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.