A multiscale mathematical model describing the growth and development of bambara groundnut.

A principal objective in agriculture is to maximise food production; this is particularly relevant with the added demands of an ever increasing population, coupled with the unpredictability that climate change brings. Further improvements in productivity can only be achieved with an increased understanding of plant and crop processes. In this respect, mathematical modelling of plants and crops plays an important role. In this paper we present a two-scale mathematical model of crop yield that accounts for plant growth and canopy interactions. A system of nonlinear ordinary differential equations (ODEs) is formulated to describe the growth of each individual plant, where equations are coupled via a term that describes plant competition via canopy-canopy interactions. A crop of greenhouse plants is then modelled via an agent based modelling approach in which the growth of each plant is described via our system of ODEs. The model is formulated for the African drought tolerant legume bambara groundnut (Vigna subterranea), which is currently being investigated as a food source in light of climate change and food insecurity challenges. Our model allows us to account for plant diversity and also investigate the effect of individual plant traits (e.g. plant canopy size and planting distance) on the yield of the overall crop. Informed with greenhouse data, model results show that plant positioning relative to other plants has a large impact on individual plant yield. Variation in physiological plant traits from genetic diversity and the environmental effects lead to experimentally observed variations in crop yield. These traits include plant height, plant carrying capacity, leaf accumulation rate and canopy spread. Of these traits plant height and ground cover growth rates are found to have the greatest impact on crop yield. We also consider a range of different planting arrangements (uniform grid, staggered grid, circular rings and random allocation) and find that the staggered grid leads to the greatest crop yield (6% more compared to uniform grid). Whilst formulated specifically for bambara groundnut, the generic formulation of our model means that with changes to certain parameter's, it may be extended to other crop species that form a canopy.

A mathematical model of the mevalonate cholesterol biosynthesis pathway.

We formulate, parameterise and analyse a mathematical model of the mevalonate pathway, a key pathway in the synthesis of cholesterol. Of high clinical importance, the pathway incorporates rate limiting enzymatic reactions with multiple negative feedbacks. In this work we investigate the pathway dynamics and demonstrate that rate limiting steps and negative feedbacks within it act in concert to tightly regulate intracellular cholesterol levels. Formulated using the theory of nonlinear ordinary differential equations and parameterised in the context of a hepatocyte, the governing equations are analysed numerically and analytically. Sensitivity and mathematical analysis demonstrate the importance of the two rate limiting enzymes 3-hydroxy-3-methylglutaryl-CoA reductase and squalene synthase in controlling the concentration of substrates within the pathway as well as that of cholesterol. The role of individual feedbacks, both global (between that of cholesterol and sterol regulatory element-binding protein 2; SREBP-2) and local internal (between substrates in the pathway) are investigated. We find that whilst the cholesterol SREBP-2 feedback regulates the overall system dynamics, local feedbacks activate within the pathway to tightly regulate the overall cellular cholesterol concentration. The network stability is analysed by constructing a reduced model of the full pathway and is shown to exhibit one real, stable steady-state. We close by addressing the biological question as to how farnesyl-PP levels are affected by CYP51 inhibition, and demonstrate that the regulatory mechanisms within the network work in unison to ensure they remain bounded.

A mathematical model of the sterol regulatory element binding protein 2 cholesterol biosynthesis pathway.

Cholesterol is one of the key constituents for maintaining the cellular membrane and thus the integrity of the cell itself. In contrast high levels of cholesterol in the blood are known to be a major risk factor in the development of cardiovascular disease. We formulate a deterministic nonlinear ordinary differential equation model of the sterol regulatory element binding protein 2 (SREBP-2) cholesterol genetic regulatory pathway in a hepatocyte. The mathematical model includes a description of genetic transcription by SREBP-2 which is subsequently translated to mRNA leading to the formation of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR), a main regulator of cholesterol synthesis. Cholesterol synthesis subsequently leads to the regulation of SREBP-2 via a negative feedback formulation. Parameterised with data from the literature, the model is used to understand how SREBP-2 transcription and regulation affects cellular cholesterol concentration. Model stability analysis shows that the only positive steady-state of the system exhibits purely oscillatory, damped oscillatory or monotic behaviour under certain parameter conditions. In light of our findings we postulate how cholesterol homeostasis is maintained within the cell and the advantages of our model formulation are discussed with respect to other models of genetic regulation within the literature.