Skip to main content

A whole-body mathematical model of cholesterol metabolism and its age-associated dysregulation



Global demographic changes have stimulated marked interest in the process of aging. There has been, and will continue to be, an unrelenting rise in the number of the oldest old ( >85 years of age). Together with an ageing population there comes an increase in the prevalence of age related disease. Of the diseases of ageing, cardiovascular disease (CVD) has by far the highest prevalence. It is regarded that a finely tuned lipid profile may help to prevent CVD as there is a long established relationship between alterations to lipid metabolism and CVD risk. In fact elevated plasma cholesterol, particularly Low Density Lipoprotein Cholesterol (LDL-C) has consistently stood out as a risk factor for having a cardiovascular event. Moreover it is widely acknowledged that LDL-C may rise with age in both sexes in a wide variety of groups. The aim of this work was to use a whole-body mathematical model to investigate why LDL-C rises with age, and to test the hypothesis that mechanistic changes to cholesterol absorption and LDL-C removal from the plasma are responsible for the rise. The whole-body mechanistic nature of the model differs from previous models of cholesterol metabolism which have either focused on intracellular cholesterol homeostasis or have concentrated on an isolated area of lipoprotein dynamics. The model integrates both current and previously published data relating to molecular biology, physiology, ageing and nutrition in an integrated fashion.


The model was used to test the hypothesis that alterations to the rate of cholesterol absorption and changes to the rate of removal of LDL-C from the plasma are integral to understanding why LDL-C rises with age. The model demonstrates that increasing the rate of intestinal cholesterol absorption from 50% to 80% by age 65 years can result in an increase of LDL-C by as much as 34 mg/dL in a hypothetical male subject. The model also shows that decreasing the rate of hepatic clearance of LDL-C gradually to 50% by age 65 years can result in an increase of LDL-C by as much as 116 mg/dL.


Our model clearly demonstrates that of the two putative mechanisms that have been implicated in the dysregulation of cholesterol metabolism with age, alterations to the removal rate of plasma LDL-C has the most significant impact on cholesterol metabolism and small changes to the number of hepatic LDL receptors can result in a significant rise in LDL-C. This first whole-body systems based model of cholesterol balance could potentially be used as a tool to further improve our understanding of whole-body cholesterol metabolism and its dysregulation with age. Furthermore, given further fine tuning the model may help to investigate potential dietary and lifestyle regimes that have the potential to mitigate the effects aging has on cholesterol metabolism.


Lipid metabolism has a key role to play in human longevity and healthy ageing. This has been emphasized by recent genetic studies examining the lipoprotein phenotype in individuals with exceptional longevity[1, 2]. It is highly probable such subjects with exceptional longevity and favourable lipoprotein profiles have avoided CVD. CVD is the primary cause of mortality in developed countries, with almost 40% of males and 30% of females in the United Kingdom over the age of 85 years living with the condition[3]. Of the components of lipid metabolism, elevated LDL-C has consistently stood out as a risk factor for CVD[4]. This cholesterol sub-fraction has been connected to atherosclerosis, a process regarded as the underlying pathogenesis for coronary heart disease (CHD) and stroke, the leading forms of CVD[5]. Studies have repeatedly demonstrated that regardless of physical activity levels and nutritional status, LDL-C has been shown to rise with age in both males and females in a diverse range of groups[6, 7]. Understanding why LDL-C rises with age is complex; nevertheless rodent studies have indicated that increases in intestinal cholesterol absorption or possibly a decrease in the plasma clearance rate of LDL-C may have a mechanistic role to play[810]. In this paper we examine the hypothesis that these two mechanisms are central to understanding why LDL-C increases with age. To test this hypothesis a number of steps were required. Firstly, it was necessary to study both intestinal cholesterol absorption and LDL-C plasma clearance within an integrative framework which incorporated the other fundamental biological components of this complex system[11, 12]. This required incorporating the interactions of the individual elements of this system. Secondly, it was important to investigate changes both to intestinal cholesterol absorption and LDL-C clearance over an extended period of time. Clearly it would have been difficult to test our hypothesis within this framework using conventional in vivo or in vitro techniques as such approaches can be resource intensive, expensive, time consuming, unpractical and potentially unethical[13]. Furthermore, whole-body cholesterol metabolism is a complex system with a variety of non-linear interactions among its various components, including both positive and negative feedback and complex crosstalk between elements such as cholesterol synthesis and absorption[14].

Utilizing a mechanistic mathematical model offered an alternative means of overcoming these difficulties and presented a cheap, ethical and practical way of investigating our hypothesis[13]. However, as existing computational models of cholesterol metabolism were either of an intracellular nature, lipoprotein kinetic focused, or compartmental in nature, they failed to incorporate the ageing process or address the holistic nature of our question and were thus determined to be unsuitable[15, 16]. Consequently we constructed a whole-body mathematical model of cholesterol metabolism which was used to explore changes to both the rate of intestinal cholesterol absorption and the hepatic rate of clearance of LDL-C from the plasma. The model showed that of these two mechanisms, changes to the rate of LDL-C removal from the plasma with age had the most significant effect on cholesterol metabolism. The model was constructed using a series of coupled ordinary differential equations (ODEs). Additionally, the model was coded in the systems biology markup language (SBML) format and submitted to the Biomodels database to facilitate its updating and future exchange[17, 18]. To build the model a wide variety of data was used including previously published data from a range of fields such as molecular biology, nutrition, physiology and biochemistry.

Whole-body cholesterol metabolism

Cholesterol has a vital role to play in the human body. It is a key constituent of all cell membranes being involved in membrane fluidity; it is the precursor of steroid hormones which control a range of physiological functions, and bile salts, which are necessary for the intestinal absorption of cholesterol, fats and lipid soluble vitamins[19, 20]. The mechanisms underpinning cholesterol metabolism interact to preserve the balance of cholesterol in the body. This balance is maintained by coordinated interactions between cholesterol absorption, excretion and synthesis.

Intake, absorption and excretion

Cholesterol from both the diet and bile is taken up by the small intestine daily. Cholesterol absorption is an inefficient process and can vary significantly from person to person. The percentage of cholesterol absorbed in healthy subjects (both male and female) is in the range of 29-80%[21]. The heterogeneity in cholesterol absorption has been attributed to its complexity, which involves a large number of enzymes and transport proteins in a multi-step process, however the way these mechanisms interact to regulate absorption efficiency remains unknown[22]. It is clear that during digestion, bile acids are released from the liver and gall bladder into the intestine and as cholesterol is practically insoluble in aqueous environments, the bile acids serve to create lipid micelles[23]. The micelles are then transported to the brush border of jejunal enterocytes, where the cholesterol is transferred into the enterocytes. At this point the story becomes nebulous, as the mechanism(s) by which micellar cholesterol is absorbed through the brush border membranes independent of bile salt uptake remains a mystery. A long standing hypothesis suggests that cholesterol absorption occurs by passive diffusion down a concentration gradient[24, 25]. However, recent evidence supports the hypothesis that protein facilitated mechanisms are involved in cholesterol uptake by the enterocyte[26, 27]. When cholesterol is inside enterocytes, it has been suggested that it can be transported back to the intestinal lumen by a class of membrane proteins known as adenosine triphosphate binding cassette transporters. It has been proposed that these transporters may serve to efflux cholesterol from the enterocyte back into the intestinal lumen for excretion[28]. Although the precise mechanisms of cholesterol absorption are not yet fully understood, it is known that cholesterol not transported back to the intestinal lumen is esterified in the enterocyte by acyl-CoA-cholesterol acyltransferase (ACAT), assembled and packaged, together with triglycerides, into chylomicrons. Chylomicrons are just one of many lipid carriers, known as lipoproteins, found in the circulation of mammals. According to their density, lipoproteins are categorised as chylomicrons, very low density lipoproteins (VLDL), intermediate density lipoproteins (IDL), low density lipoproteins (LDL) and high density lipoproteins (HDL)[29]. Chylomicrons transport absorbed cholesterol via the lymph to the liver[30].

The liver and whole-body cholesterol metabolism

The liver is the central organ involved in cholesterol metabolism. It is actively involved in the uptake of cholesterol from lipoproteins, is an important site for de novo synthesis and is capable of storing cholesterol as esters after esterification by ACAT. Furthermore, it is involved in the secretion of cholesterol-containing lipoproteins and is the sole organ capable of removing excess cholesterol from the body, either by secretion into bile or by conversion into bile acids[31]. Bile acids are synthesised from cholesterol in the liver and facilitate the solubilisation of cholesterol, a prerequisite for cholesterol absorption. Bile acids are released postprandially and will not solubilise dietary lipids unless above a critical concentration. The concentration gradient is generated by two factors. Firstly, bile acids are powerful acids that are impermeable to cell membranes. Secondly, the majority of bile acids are reabsorbed from the small intestine and return to the liver via the portal vein, where they are taken up by hepatocytes and re-secreted into bile[32]. The conservation of bile acids by active absorption from the small intestine results in a hepatic pool of bile acids that cycles several times with each meal. In humans, the bile acids circulate between six and 10 times per day, while around 400 mg of bile acids is lost daily through faecal excretion. A decreased return of bile acids to the liver is compensated for by increased de novo synthesis from cholesterol in order to maintain the bile salt pool[33]. Cholesterol not converted to bile acids or secreted into bile can be released from the liver into the circulation and transported to peripheral tissue. The liver secretes VLDLs and on entering the circulation VLDLs are hydrolysed by lipoprotein lipase (LPL)[34], leading to the formation of VLDL remnants and IDLs. IDLs are either taken up by the liver or further hydrolysed to LDLs, which are the main cholesterol carrier in the blood[35]. LDL is taken up by the liver or by peripheral cells, either independently or via the LDL-receptor (LDLR). Expression is high in the liver but LDLRs are also expressed in peripheral tissue[36, 37]. Hepatic LDLR (HLDLR) expression is transcriptionally regulated in response to intracellular cholesterol levels[38].

Reverse cholesterol transport

The flux of cholesterol from peripheral tissue to the liver is known as reverse cholesterol transport (RCT). This pathway plays a vital role in maintaining cholesterol balance, with its action resulting in the movement of cholesterol from peripheral tissue to the liver[39]. Consequently RCT presents the only route for excess cholesterol generated in peripheral tissue to be eliminated from the body. HDL plays a crucial part in RCT[40, 41]. HDL synthesis takes place in the liver and intestine, commencing with the generation of nascent particles, small discoidal structures lacking cholesterol. HDL acquires free cholesterol from the cell membranes of several tissues, which is esterified by the action of lecithin:cholesterol acyltransferase (LCAT), resulting in cholesteryl ester-rich, mature HDL particles. The transport of high density lipoprotein cholesterol (HDL-C) to the liver may follow several routes. In the presence of the enzyme cholesteryl ester transport protein (CETP), a portion of cholesterol esters from the HDL particle can be redistributed to other lipoproteins (e.g. VLDL and LDL) in exchange for triglyceride[42]. This cholesterol is delivered back to the liver via the LDLR. An alternative route involves the mature HDL particle binding directly to scavenger receptors BI (SR-BI) on the liver[43].

Cholesterol synthesis

Almost all the tissues in the body are capable of synthesizing cholesterol from acetyl-Co enzyme A (CoA) in a series of enzyme-mediated steps that are mainly restricted to the endoplasmic reticulum (ER). The major rate-controlling enzyme in this pathway is 3-hydroxy-3-methylglutaryl CoA (HMG-CoA) reductase. Cholesterol homeostasis is regulated in two ways, firstly, when the content of unesterified cholesterol in cells increases, the expression of the LDLR protein decreases. Furthermore, two key enzymes involved in cholesterol synthesis are repressed, thus, any subsequent increase in cellular cholesterol diminishes. When cholesterol levels fall, these pathways are reactivated[38, 44, 45].

Dysregulation of cholesterol metabolism with age

The interaction between cholesterol absorption, synthesis, and excretion maintains whole- body cholesterol metabolism. A change to any of them can significantly influence the others. For example, numerous studies have shown that inhibition of cholesterol synthesis results in increased intestinal cholesterol absorption[46, 47]. The close interaction between absorption and synthesis in maintaining cholesterol balance is consolidated by studies, which have shown that interference with cholesterol absorption is associated with increased cholesterol synthesis[47]. Dysregulation of cholesterol metabolism can lead to high levels of LDL-C, while population studies have consistently demonstrated that LDL-C rises with age in both males and females[7]. The reason for an increase in LDL-C in so many individuals across both genders remains unknown. Furthermore, the intrinsic mechanisms that contribute to the rise in LDL-C with age remain incomplete; however a picture is beginning to emerge of potential candidates. A significant finding was the discovery that there is a gradual reduction in the rate of clearance of LDL-C from the circulation with age[4850]. Adding further to this is the evidence that the number of hepatic LDLRs (HLDLRs) diminish with age in certain species[10, 50]. Additionally, it has been reported that there is an increase in intestinal cholesterol absorption efficiency with age in rats[8, 51]. However, to date no single intrinsic mechanism or combination of intrinsic mechanisms has conclusively been attributed to the rise in LDL-C with age[9]. Consequently our whole-body mathematical model investigated the hypothesis that age-related changes to the rate of intestinal cholesterol absorption and the rate of removal of LDL-C from the plasma could be responsible for the rise in LDL-C with age in humans.


Results overview

Using the values listed in tables one and two, a sensitivity analysis was conducted, followed by a number of simulations that examined the key questions associated with this research. The model was able to replicate many of the integral features of whole-body cholesterol balance and displays qualitative behaviour similar to that observed in experimental studies[52]. For example, plasma LDL-C levels rose in response to increases in the intake of dietary cholesterol, however this was somewhat more sensitive than that predicted by a meta-analysis of the experimental literature, as for every 100mg/day increase in dietary cholesterol there was a ≈10mg/dL increase in LDL-C (Figure1). Changes to cholesterol absorption were also investigated and the model showed that for every 10% increment in the rate of cholesterol absorption, this resulted in a 12.5 mg/dL increase in LDL-C. Finally, the model was used to investigate the hypothesis that age related changes to cholesterol absorption/a decrease in the clearance of LDL-C from the circulation result in the rise of LDL-C with age. The model showed that of these two mechanisms, changes to the rate of removal of LDL-C from the circulation had the most significant impact on LDL-C levels with age.

Figure 1

(A-F) The response of the model to changes in the intake of dietary cholesterol compared to a published meta-analysis. A) The relationship between mean change in dietary cholesterol (mg/day) and mean change in plasma cholesterol levels (mg/dL) in 167 cholesterol feeding studies published between 1960 and 1999 (McNamara 2000). B-D) The response of the model to changes in dietary cholesterol over a range of time periods. E) Steady-state levels of LDL-C for various different intakes of dietary cholesterol. F) The change in LDL-C levels in response to various different values of the parameter k 1 (cholesterol ingestion). Results of both E&F were generated from a sensitivity analysis of the model using the software tool Copasi.

Predictive capability of the model

The model demonstrated that increasing cholesterol absorption from 50% to 80% and running a simulation from age 20 years to age 65 years resulted in an increase of plasma LDL-C by as much as 34 mg/dL (Figure2, graph A). However, reducing the number of hepatic LDL receptors had a profound impact on the system, as an increase of 116 mg/dL in LDL-C was witnessed by age 65 years in response to a reduction in the number of these receptors by 50% (Figure2 graph B). These rises were compared to that found in population studies (Figure2 graphs C and D) with changes to the rate of clearance of LDL-C clearly having the more significant impact on LDL-C levels.

Figure 2

(A-E) The relationship between LDL-C and age. The relationship between LDL-C and age A) for males and B) for females. Graphs A and b were both generated using data from Abbott, Garrison et al. (1983). C) the response of LDL-C in the model to increased cholesterol absorption. D) The response of LDL-C in the model to decreases in the rate of clearance of LDL-C by hepatic LDLRs. E) Parameter scan for k hrs (daily rate of synthesis of hepatic cholesterol receptors) and its impact on LDL-C levels.

Sensitivity analysis-cholesterol absorption

The initial concentrations of the various species are summarised in Table1. Where quantitative data was available, these values reflected the numbers found in the literature. A sensitivity analysis of the model was also conducted as the parameter values differ considerably in terms of their magnitude (Table2) and we wanted to establish how changes to these parameters would affect the concentrations of key species. Cholesterol absorption (k 6 ) is an example of one of the key parameters that was altered. The initial value of this parameter was set so that 50% of cholesterol entering the small intestine would be absorbed. However, although cholesterol absorption is tightly regulated, it is a process that can vary significantly from one individual to the next[21]. This range was used to test the sensitivity of the model to changes in k 6 . The impact variations in cholesterol absorption has on LDL-C levels in the model is shown in (Figure3).

Table 1 List of species, their abbreviations and their initial values
Table 2 Summary of the reactions and parameter values used in the model
Figure 3

(A-E) Changes to cholesterol absorption in the range 30-85%. A-D) simulations of the model using MathSBML to show the response of LDL-C to changes in the rate of efficiency of cholesterol absorption in the range 30-80% over a number of different time periods. E) Steady-state levels of LDL-C for various different percentages of cholesterol absorption. This was generated with the software tool Copasi.

Sensitivity analysis-dietary cholesterol

A literature review using PubMed revealed the existence of almost 200 publications containing information relating to the response of plasma cholesterol to dietary cholesterol feeding. These investigations were performed in both young and old individuals. Importantly among these publications, a meta-analysis by McNamara (2000) revealed that in response to increased levels of dietary cholesterol, plasma cholesterol plateaus at high levels of cholesterol intake[52]. The meta-analysis by McNamara also reported that the total plasma cholesterol response to dietary cholesterol is 0.023 mg/dL per mg/day cholesterol increase (Figure1, graph A). Of this, 0.019 mg/dL is in the LDL-C fraction. Therefore, for each 100 mg/day increase in dietary cholesterol intake in the model, a 1.9 mg/dL increase in LDL-C was expected. Importantly, McNamara reported that this increase is independent of other types of lipid in the diet and the baseline plasma cholesterol level. Therefore, the model was exposed to a decrease of 200 mg/day and increases of up to 700 mg/day of dietary cholesterol (Figure1, graphs A-C). A ≈10 mg/dL increase in LDL-C was observed for every 100mg increase in dietary cholesterol which was 8mg greater than the value predicted in the meta-analysis by McNamara. We also conducted a sensitivity analysis of the parameter k 1 to highlight the different steady-state levels of LDL-C for various intakes of dietary cholesterol (Figure1, graphs E&F and Additional file1). No increase in HDL-C was observed in response to changes in dietary cholesterol.

Hypothesis testing- ageing and cholesterol absorption/LDL-C plasma removal

Aging results in a rise in LDL-C in both sexes (Figure2, graphs C and D)[7]. The reason(s) for an increase in LDL-C in so many individuals across both genders remains unknown and the issue is the subject of debate as nutritional status and physical activity levels may differ within the population sample. The intrinsic mechanism(s) that contribute to the rise in LDL-C with age remain incomplete; however a picture is beginning to emerge of potential candidates. As mentioned previously, rodent studies have shown that cholesterol absorption efficiency increases markedly with aging[8, 21]. It was also mentioned that there is a gradual reduction in the rate of clearance of LDL-C from the circulation with age in rodents. Adding further to this is the evidence that the number of hepatic LDLRs diminish with age in certain species[29]. When cholesterol absorption efficiency was increased in 10% increments between 50% and 80% by age 65 years, 12, 24 and 34mg/dL increases respectively, in LDL-C were observed (Figure2, graph A). Next, reducing the number of hepatic LDL receptors (HLDLRs) was investigated. The number of HLDLRs was reduced in the range 90–50 by age 65 years. This was done by gradually decreasing the rate of synthesis k hrs . This had a significant effect on the model as LDL-C raised by 116,76,49,26 and 11mg/dL respectively (Figure2, graph B). A parameter scan of k hrs was also conducted (Figure2, graph E).


We have constructed a whole-body mathematical model of cholesterol metabolism using data from a wide variety of sources and integrated this data within a series of coupled ODEs. Previous computational models have focused on the intracellular regulation of cholesterol metabolism[16, 53] or have focused on the metabolic fate of a particular lipoprotein or their receptor mediated endocytosis[54, 55], while other mathematical models of lipid metabolism have centred on compartmental approaches, where by lipoproteins are represented by compartments[56]. Our model attempts to include all the major components of whole-body cholesterol balance and to our knowledge is the first model of its type built using such a framework. Our model provides an insight into the complex interplay of cholesterol metabolism with the aging process. The model has many areas that require further development, however, when the model progresses in the future it would be worthwhile to investigate how changes to different combinations of parameters affect the overall behaviour of the system, as the model indicates that the dysregulation of cholesterol metabolism with age may involve perturbations to several components of this system. In tandem with investigating changes to combinations of parameters, the model may also be used to investigate combination therapy that may help to mitigate the effects aging has on cholesterol metabolism. For example, it would be straightforward to include hypothetical interventions that investigate different dietary regimes and also include the effects of variations in physical activity. For example, based on available literature the simple assumption could be made that consuming 3 g/day of fibre reduces cholesterol absorption by 15%[57, 58]. Therefore, the parameters could be adjusted accordingly to reflect this. Similarly the effect of consuming plant sterols, which are known to reduce cholesterol absorption[59, 60] could easily be investigated, as consumption of 1.8 to 2.0 g/day of plant sterols has been shown to lower both total and LDL-C concentrations by 10% to 15%, respectively in a variety of different population groups[6163]. Thus, this could be investigated by making alterations to the rate of cholesterol absorption in the model.

As mentioned previously the model is by no means the finished product and several assumptions were made during the building of the model, which were necessary as considerable uncertainty still surrounds cholesterol metabolism. For example, in the majority of cases the functional form that best describes a reaction remains unknown. We did however include a number of feed-back and feed-forward equations to describe the behaviour of cholesterol and bile salt synthesis as these are known to be subject to this type of action. At steady state, parameter values attempted to reflect a generic normolipididemic 20 year old male. Although parameter values were chosen for the model from available published literature, the ranges of parameter values are quite broad; however, given the nature and diversity of the literature used this was unavoidable. On certain occasions, parameter values simply did not exist, therefore reasonable values were chosen based on the published literature. An area of the model that could be developed further in the future is the kinetic parameters of the various enzymes involved in a number of the reactions. Presently, most of these enzymes have arbitrary values of 100 assigned to them to indicate that they have ‘normal’ enzymatic activity. In the future it would be worthwhile to use a resource such as BRENDA to amend these enzymes to include km and kcat values[64]. Despite these limitations it is hoped that in the future the model will prove to be a useful tool for testing nutritional and lifestyle interventions that would be difficult, challenging or perhaps unethical to conduct using convention means.

Building the model has highlighted priorities for future experimental work, for example the need for experimentalists to work closely with computational modellers in order to generate appropriate quantitative data that is of benefit to systems biology models such as this one. Such data could be used to expand and enhance our model leading to a better understanding of cholesterol metabolism, particularly the interrelationship between cholesterol synthesis and absorption. This improved model could be used to gain additional insights into the factors associated with cholesterol metabolism that help to confer increased longevity and healthy aging. Appropriate dietary intervention strategies based on such models could be used to decrease the risk of CVD and prolong healthy aging. Models could be tailored for an individual’s dietary, biochemical, genetic and social circumstances and as such models can be coded in SBML they would be relatively easy to extend, enhance and improve as more quantitative data becomes available. Furthermore, it would be worthwhile combining our model with existing models, such as that of atherosclerosis progression[65]. The portable nature of SBML also makes this a very realistic possibility, especially if models that are being merged have been designed using SBML[66, 67]. Moreover, the model could be used to investigate crosstalk between cholesterol metabolism and fatty acid metabolism.


From the evidence presented in this paper it is clear that cholesterol metabolism is a complex multi-component system. Changes to this multi-component system as a result of the ageing process can lead to high levels of plasma cholesterol particularly LDL-C. There is an indisputable link between elevated levels of LDL-C and the risk of developing atherosclerosis, with population studies indicating that LDL-C rises with age in both sexes in a large number of individuals[5, 68, 69]. It is therefore unsurprising that recent genetic studies have associated lipid metabolism with longevity as it is apparent that maintaining cholesterol balance is vital to cardiovascular health[1, 2, 70, 71].

We developed a whole-body computational model to investigate the hypothesis that an increase in the rate of cholesterol absorption and a reduction in the rate of removal of LDL-C from the plasma were integral to understanding the dysregulation of cholesterol metabolism with age. We were of the opinion that existing models of cholesterol metabolism were limited and could not offer a sufficiently deep understanding of how an entire biological system changes with age, thus in this work a conscious effort was made to represent the interactions throughout the body. Results from the sensitivity analysis of the model indicated that the model requires further fine-tuning in the future before it could be used as a comprehensive tool for predicting changes to LDL-C in response to dietary cholesterol/fats and intrinsic aging. None the less we used the model to investigate aging and it demonstrated that changes to the rate of removal of LDL-C to half its original value at age 65 years resulted in a 116 mg/dL rise in this cholesterol sub fraction, which is intimately connected with CVD. The model contrasts with other models of cholesterol metabolism in that it is implemented in a whole-body fashion and provides a template for building a quantitative systems level understanding of cholesterol metabolism and its interaction with aging. In the future the model could contribute to a better understanding of cholesterol metabolism so that the wider population could benefit in the same way as those individuals with exceptional longevity have from an altered lipid profile.


Deciding on model structure and the rationale for a whole-body approach

Historically, mathematical modelling has been used to investigate various aspects of both cholesterol and more generally lipid metabolism; however to our knowledge no model to date has been implemented to investigate cholesterol metabolism mechanistically within a whole- body framework, with the goal of understanding how intrinsic age-related biological changes affect this crucial biological system. In order to appreciate the rationale underpinning the whole-body approach that we adopted, it is worthwhile examining some of the models of cholesterol metabolism that have been developed previously. Additionally, it is worth exploring why these were unsuitable for this work. In general these models can be divided into three types; intracellular models of the cholesterol homeostasis genetic regulatory pathway; models of lipoprotein dynamics, and compartmental models. If we take the gene network approach, a recent version involved a boolean network model that assigned 33 different components of intracellular cholesterol homeostasis a boolean value. Simulations then led to the representation of intracellular cholesterol homeostasis as a boolean vector with each coordinate of the system denoting a biological species of the pathway[16]. It is easy to appreciate the utility of this approach as a means of better understanding the functioning of a gene regulatory network; however this strategy is not a mechanistic approach. It is primarily a state based representation of intracellular gene activity; therefore it is limited as it does not include the different biological compartments involved in cholesterol metabolism, such as absorption, and hepatic/peripheral cholesterol synthesis. Nor does it account for the synthesis of the various lipoproteins such as VLDL, IDL, LDL and their subsequent interactions with one another and with receptors both hepatically and peripherally. Finally, and significantly for this work it is difficult to imagine such models being able to test a biological hypothesis that centres on three factors, the aging process, cholesterol absorption and the removal of LDL-C from the plasma.

If mathematical models focused on lipoprotein kinetics are examined; one discovers a wide variety of models that have represented varying degrees of cholesterol metabolism and have centred to a large extent on the dynamics of the LDL receptor. Examples of such models go back as far as the early 1980s when Goldstein et al. presented a theoretical study of the interaction of LDL receptors with coated pits[72, 73]. These were followed by other models of the same process[73] or more recent slight adaptations of this process[74]. These models have been beneficial from the point of view of elucidating the underlying dynamics which describe the interaction of LDL with its receptor and the associated underlying kinetics of this process. However, as with the intracellular models, these models focus on one particular, isolated aspect of cholesterol metabolism and do not portray the entire picture which is that cholesterol metabolism involves the coordinated action of several different biological systems operating in unison throughout the entire body.

The final approach that will be discussed is the compartmental methodology. Compartmental models, as the name suggests, contain a number of compartments, each containing a well mixed material. Each compartment has a number of connections leading into it and out of it. Biological material can flow from one compartment to another, and it can be added from the outside through a source, or it can be removed through a sink. Such systems exchange material with each other following certain rules[75]. These models have advantages such as being able to examine in some detail the interaction of lipoproteins and their derivatives. The major limitation of these models is that they are not based on underlying biological mechanisms; rather exchange rates are based on experimental inferences. Therefore, as with the previous approaches this method is also unsuitable for representing a mechanistic whole-body of cholesterol metabolism.

In order to overcome the limitations of the approaches that have been discussed, it was decided to represent cholesterol metabolism with the aid of biological compartments that represented their real-world biological counterparts. The compartments contained quantitative details of the different forms of cholesterol. These forms are present in very high numbers (e.g. intestinal, hepatic and peripheral cholesterol is present in grams). This means that a stochastic approach would not offer significant benefits, as the intrinsic stochasticity in the system is negligible, and so a deterministic model was considered more appropriate to represent these species and their interactions. In order to implement this deterministic solution, we created a series of ordinary differential equations that were coupled and numerically integrated. This type of approach has previously been shown to be eminently suitable both for representing mechanistic biological reactions and for dealing with the time dimension associated with modelling the aging process[13]. In order to build the ODEs, we first needed to concisely represent the compartments and the fundamental components of whole-body cholesterol metabolism diagrammatically.

Construction of a network diagram

Whole-body cholesterol metabolism is a difficult process to reason about directly; therefore the first step in constructing a model of this nature was to assemble a network diagram (Figure4) using Systems Biology Graphical Notation (SBGN). SBGN is a standardized means of graphical representation of biological systems that facilitated the drafting of a network diagram of cholesterol metabolism in a clear and unambiguous fashion[76]. This network diagram explicitly outlines each of the compartments and the biological species involved in whole-body cholesterol metabolism which we have discussed. Each of the components of Figure4 was then abbreviated and used to construct a series of ODEs (Table1 contains a summary of these abbreviations).

Figure 4

Network diagram of the cholesterol metabolism model. The model is laid out in 6 compartments, 1) intake, 2) intestinal tissue, 3) excretion, 4) plasma, 5) hepatic and 6) peripheral tissue. The arrows represent the flow of cholesterol around the system into its different forms. Enzymes are represented by blue spheres and their catalytic influence on the reaction is indicated by a round arrow head coming into contact with a reaction arrow. Synthesis is represented by the Greek symbol theta, while inhibition is represented by T -shaped arrows. In summary cholesterol from the diet and bile is formed into micelles in the small intestine. Absorbed cholesterol is then transported to the liver where it is exported into the plasma via VLDL. VLDL is in turn catabolised to LDL. Excess cholesterol from peripheral tissue is transferred to the liver via HDL.

Model assembly

Rationale for reaction forms

The assumption was made that a large number of the reactions are of a first-order nature, where the rate of the reaction is directly proportional to the concentration of one of the reactants. This assumption was made to facilitate the evaluation of the model based on steady states and to allow for the determination of parameters when the model was compared to experimental data, moreover the reaction form that best describes a large number of these remains unknown. ODEs were constructed in an incremental fashion to meet certain steady states, eventually giving rise to a unified series of coupled ODEs and ultimately a whole-body mathematical model representing a generic twenty year old male. The model was built using MathSBML[77]. MathSBML is a software package that has been designed to work with Mathematica (version 5.2)[78], a pre-existing commercial software package developed to perform numeric, algebraic, graphical and many other tasks. In addition the software tool Copasi was used to perform some of the parameter scans and steady state examinations of the model[79]. The next section outlines the assembly of the ODEs and the biological functions that lie at their core. Additionally, the reasoning for the inclusion of each of the species in Figure4 will be detailed along with the quantitative data that was used to make inferences about each of the steady states that these species gave rise to.

Cholesterol intake

As whole-body cholesterol metabolism was being dealt with, it was important first to include the dietary intake of cholesterol. According to Henderson et al. the average amount of cholesterol ingested daily in the UK by a male is 304 mg/day[80]. Cholesterol intake was assigned the rate k1 and set at an initial value of 1. While dietary cholesterol [DC] had its initial condition set to a fixed value of 304.

Cholesterol synthesis

Dietary cholesterol joins synthesised intestinal cholesterol. Additionally, cholesterol is also synthesised throughout the body[81]. To determine the daily rates of synthesis in each compartment, whole-body cholesterol synthesis was first calculated. According to Dietschy et al. humans synthesise ≈10 mg/kg of cholesterol daily, therefore a 70 kg man synthesises ≈700 mg/day[81]. The amount of cholesterol synthesis that can be attributed to the liver, intestine and remaining peripheral tissue is uncertain[81]. Therefore, whole-body synthesis of cholesterol for a 70 kg man was based on the assumption that ≈70% of cholesterol synthesis takes place in the peripheral tissue and that 10% of this occurs in the intestine[82] (Table3). This data was then used to make inferences about inter-compartmental cholesterol synthesis and functional relationships were derived to describe cholesterol synthesis in each compartment. Firstly, intestinal cholesterol synthesis was represented with a negative feedback function (equation 1) were ICSmax represents the maximum amount of cholesterol that can be synthesised by the intestine in a day, IC represents intestinal cholesterol. ICt represents the intestinal cholesterol synthesis threshold and IS is a fitting parameter that represents the sensitivity of the negative feedback. It was assumed that the maximum rate of intestinal cholesterol synthesis would be 100 mg/day (Figure5, graph A).

Intestinal cholesterol Synthesis = I C S max 1 + I C / I C t IS
Table 3 Breakdown of whole-body cholesterol synthesis
Figure 5

(A-D) Functions and simulations associated with model building. A) Graph of concentration of intestinal cholesterol versus cholesterol synthesis. C and D) Graphs of the functions for hepatic and peripheral cholesterol synthesis respectively. D) Graph of concentration of hepatic cholesterol versus biliary cholesterol release state values being reached for IDL-C and VLDL-C, respectively.

We also derived similar functions for both hepatic and peripheral synthesis (equations 2 &3, respectively). HCSmax and PCSmax represent the maximum amount of cholesterol that can be synthesised daily in the hepatic and peripheral compartments respectively. HCSt and PPCt represent hepatic and peripheral daily rates of synthesis thresholds respectively and HS and PCSS are fitting parameters that represent the sensitivity of the feedback in the two respective compartments. Based on experimental data it was assumed that the maximum daily rate of cholesterol synthesis in the liver would be 500 mg/day, while the maximum rate of peripheral cholesterol synthesis was also estimated to be 500 mg/day (Figure5, graphs B and C).

Hepatic cholesterol synthesis= H C S max 1 + H F C / H C S t HS
Peripheral cholesterol synthesis= P C S max 1 + P F C / P P C t PCSS

Hepatic biliary cholesterol release

According to Grundy and Metzger dietary cholesterol also mixes intestinally with ≈ 1200 mg/day of biliary cholesterol[83]. Therefore, a functional relationship was derived to describe the release of biliary cholesterol. This was done using a feed-forward equation (equation 2) as by nature this biological interaction is feed-forward. BCRmax represents the maximum daily release of biliary cholesterol. BCRt represents the threshold of release of biliary cholesterol and HFC represents the hepatic pool of free cholesterol. Based on experimental data it was assumed that the maximum rate of release of biliary cholesterol would be 2000 mg/day (Figure5, graph D)[83].

Release of biliary cholesterol= B C R max 1 + B C R t / H F C BS

Estimating the steady-state levels of cholesterol in each compartment

It was necessary to estimate the steady-state level of cholesterol in each compartment. According to Soars et al. in mammals this usually falls within the range of 100-200 mg/100 g of intestine. Also according to Soras et al. humans contain 30 g of intestine per kg of body weight[84]. The concentration of cholesterol in the intestine was determined based on several assumptions; firstly it was assumed that there is 150 mg of cholesterol/100 g of intestine and secondly the assumption was made that the hypothetical individual weighed ≈70 kg. Details of this calculation are provided in Table4.

Table 4 Calculating the steady-state level of intestinal cholesterol

Bile salts- hepatic synthesis and enterohepatic circulation

Bile salts are synthesised from the hepatic pool of free cholesterol [HFC] at a rate of ≈400 mg/day, while the hepatic bile salt pool [HBS] contains bile salt in the region of ≈4000 mg. This pool is released on average four to six times per day into the small intestine. Thus, a total of ≈24000 mg of bile salts enter the small intestine daily. The synthesis of bile salts is subject to negative feedback as the bile salt pool increases. Based on this information a reciprocal type function was derived to represent the production of bile salts (equation 5). Where k 5 is the rate constant for the hepatic production of bile salts.

Hepatic Production of Bile Salts= k 5 HFC HBS

Almost all bile salts are reabsorbed (≈23600 mg/day) and return to the liver, while the remainder (≈400 mg/day) are excreted[85]. Bile salt release was assigned the rate constant k 2 , while the hepatic return of bile salts was assigned the rate constant k 3 . This rate was dependent on the concentration of intestinal bile salts [IBS] and the concentration of [IC]. The excretion of bile salts [EBS] was assigned the rate constant k 4 . Putting this information together, ODEs 1 and 2 were assembled to describe how hepatic and intestinal bile salts change with time.

d H B S d t = k 5 H F C H B S + k 3 I B S I C - k 2 H B S ODE 1
d I B S d t = k 2 B S P - k 3 I B S I C - k 4 I B S ODE 2

Cholesterol absorption and excretion

The rate of cholesterol absorption and excretion depend on the concentration of intestinal bile salts. Over the years there has been controversy as to whether there is a difference between the absorption of dietary cholesterol and endogenously derived cholesterol. It is now accepted that they form an indistinguishable intestinal pool[54]. An assumption was made based on experimental data that 50% of cholesterol in the small intestine whether it originates from de novo synthesis, diet or bile will be absorbed in a normal individual each day while 50% would be excreted[21, 86]. Table5 presents a summary of the estimated flux of cholesterol into and out of the small intestine on a day to day basis. The daily rates of cholesterol absorption and excretion were represented by the rate constants k6 and k7 respectively, and thus we were able to derive ODEs 3 and 4 that represent the change with time of IC and excreted cholesterol EC, respectively.

d I C d t = k 1 D C + B C R max 1 + B C R t / H F C BS - k 6 I B S I C - k 7 I B S I C + I C S max 1 + I C / I C t IS ODE 3
d E C d t = + k 7 I B S I C ODE 4
Table 5 Flux of cholesterol into and out of the small intestine on a daily basis

Simulations of the mini model

At this stage it was critical to create a mini-model, as determination of parameter values becomes increasingly complicated as models gain in size and complexity. The parameters for the mini model were determined using the baseline data outlined and the ODEs were numerically solved using MathSBML. Figure6, graph A, shows the amount of cholesterol and bile salts excreted over a 100 day period. This output from the mini-system confirmed that the intestinal and intake compartments, and part of the hepatic compartment were behaving in a biologically realistic manner, thus further compartments and species were added to the system. It is important to note that the hepatic pool of cholesterol was fixed at a value of 100 during these simulations to facilitate parameterisation.

Figure 6

(A-E) Outputs and steady states associated with model building. A) Graph of the output from a simulation of the mini model that was constructed initially. It shows the faecal excretion of both cholesterol and bile salts over a 100 day period. B and C) Graphs showing the steady state values being reached for both the hepatic cholesterol pool and hepatic cholesterol esters respectively. D and E) Graphs showing steady state levels of LDL-C, IDL-C and VLDL-C.

Expansion of the mini model- further additions to the hepatic compartment

It was now necessary to estimate the steady state level of hepatic cholesterol. The amount of cholesterol in the liver was determined using data from Sahlin et al., the molecular weight of cholesterol and the weight of an average human liver[87]. Calculations are outlined in Table6. According to Cook humans weighing ≈70 kg contain ≈ 140 g of cholesterol[88]. Based on the calculations in Table6 the assumption was made that half the cholesterol in man resides in the liver, while the remainder resides in peripheral tissue. The assumption was also made that 10 g of cholesterol is stored as esters in each compartment. The hepatic synthesis of cholesterol was previously estimated to be 210 mg/day. When the hepatic pool of cholesterol increases, there is a rise in the conversion of cholesterol to esterified cholesterol. This reaction is catalysed by ACAT and the forward reaction was assigned the rate constant k 9 . The reverse reaction is catalysed by cholesterol ester hydrolases (CEH) and was assigned the rate constant k 10 . Both ACAT and CEH were given fixed vales of 100 to represent ‘normal’ enzymatic activity. Thus, based on this information we were able to derive ODEs 5 and 6 to describe the change in both the [HFC] and hepatic cholesterol esters [HCE] with time.

d H F C d t = k 19 + k 10 H C E C E H - B C R max 1 + B C R t / H F C BS
- k 12 H F C - k 5 H F C B S P - k 9 H A C A T H F C
H C S max 1 + H F C / H C S t HS ODE 5
d H C E d t = k 9 A C A T H F C - k 10 H C E C E H ODE 6
Table 6 Outlining how the concentration of cholesterol in the liver was determined (Non obese subjects (n=19))

These ODEs were added to the system and the model was re-parameterised by sequential adjustment of the parameters in order to bring both hepatic cholesterol esters and the hepatic pool of free cholesterol into a steady-state (Figure6, graphs B and C).

Hepatic release of very low density lipoprotein cholesterol

According to Glomset et al. the amount of cholesterol secreted hepatically in VLDL is considered to be in the region of 800 mg/day and could be as high as 3000 mg/day[90]. Therefore, it was assumed that in the steady state ≈1000 mg/day of cholesterol was incorporated into VLDL. The rate of incorporation of cholesterol into VLDL was assumed to be proportional to the concentration of hepatic cholesterol and this reaction was assigned the rate constant k12.

The turnover of hepatic LDL receptors

According to Dietschy, the rate of lipoprotein removal from the circulation by the liver is dependent on the concentration of the lipoprotein and hepatic LDL receptors (HLDLRs)[82]. Brown and Goldstein established that the synthesis of HLDLRs is subject to inhibition by the intracellular concentration of free cholesterol[44, 91]. Thus represented the synthesis of HLDLRs with equation 6, where khrs is the rate constant for the synthesis of HLDLRs.

Synthesis of hepatic LDL receptors = k hrs H L D L R s H F C

It is challenging to determine the actual number of HLDLRs present on a human organ such as the liver. Although, the number has been estimated at a cellular level, it was not feasible to try and determine it for a whole-body model of cholesterol metabolism. Thus, an arbitrary value of 100 was assigned to represent a fully active complement of HLDLRs and the rate constant k 13 was assigned to represent the rate of HLDLRs degradation. Thus, ODE 7 was derived to describe the turnover of HLDLRs. It is important to note that this type of approach has been used previously to model this process, whereby the regulation of receptor synthesis is handled via reciprocal feedback due to cholesterol concentration[92].

d HLDLRs d t = k hrs HLDLRs H F C - k 13 HLDLRs ODE 7

It is known that a certain amount of VLDL-C is removed via HLDLRs[82]. It was difficult to obtain quantitative data about this process; however the literature indicated that this is not the major route by which cholesterol re-enters the liver. It was therefore assumed that only a negligible amount of VLDL-C is taken up via this avenue (≈100 mg/day) and was assigned the rate constant k14.

Hepatic LDL receptor mediated removal of VLDL-C from the plasma

The enzyme LPL acts on VLDLs to reduce them to IDLs[93]. It was assumed that this process is dependent on the concentration of LPL and of VLDL-C and was assigned the rate constant k15. Thus, ODE 8 was derived to represent the change with time of VLDL-C.

d V L D L C d t = k 12 H F C k 14 HLDLRs VLDLC k 15 LPL][VLDLC] (ODE 8

Similar to VLDL-C, a certain amount of IDL-C re-enters the liver via HLDLRs and the assumption was made that this is ≈ 10% of the overall particle composition. IDLs are further catabolised to LDLs by the enzyme hormone sensitive lipase (HSL). The majority of LDL-C (≈ 75%) is taken up by the liver, while the remainder is removed by peripheral tissue. Based on the literature it was assumed that receptor independent uptake of LDL-C is of minor importance and accounts for only 50 mg/day of LDL-C removed from the tissue[9496]. This value was split equally between the liver and peripheral tissue. Based on this information ODEs 9 and 10 were derived to describe the change in IDL-C and LDL-C with time. k 16 is the rate constant for removal of IDLC via the hepatic receptor and k 17 is the rate constant for the formation of LDL-C, while k 18 is the rate constant for HLDLR uptake of LDL-C and k 19 is the rate constant for the hepatic receptor independent uptake of LDL-C. Peripheral LDL receptor (PLDLRs) and non-peripheral receptor uptake are represented by the rate constants k 20 and k21, respectively.

d I D L C d t = k 15 L P L V L D L C - k 16 H L D L R I D L C - k 17 H S L I D L C (ODE 9 )
d L D L C d t = k 17 H S L I D L C - k 18 L D L C H L D L R s - k 19 L D L C - k 20 L D L C P L D L R s - k 21 ODE 10

Parameterisation to establish the steady states of LDL-C, IDL-C VLDL-C

As the system had now increased in size significantly, further simulations were completed to establish the steady-state levels of the various species. Again the model was re-parameterized using the data outlined and the system numerically solved with MathSBML. After adjustment of the parameters LDL-C, IDL-C and VLDL-C entered steady states based on known literature values (Figures6, graphs D and E). It is important to note that the values of both the peripheral pool of free cholesterol (PFC) and PLDLRs were held at fixed arbitrary values of 100 to facilitate this process.

Further additions to the peripheral compartment

Once the first three compartments demonstrated reliable biological outputs, the final phase of model building involved the addition of the peripheral compartment. Table7 outlines how the amount of cholesterol in the peripheral tissue was calculated based on previous calculations for the amount of cholesterol in the liver and intestine. Previously, the rate of peripheral cholesterol synthesis was estimated to be 441 mg/day. This was completed in exactly the same manner as for the hepatic pool, where the rate constants k prs represents the synthesis of PLDLRs, k 22 represents the degradation of PLDLRs, k 23 represents the formation of peripheral cholesterol esters (PCE) and k 24 is the release of stored peripheral cholesterol. Additionally, according to Myant, a certain amount of cholesterol can be converted to steroid hormones. For example, a normal man may excrete 25-50 mg/day of total adrenocortical and gonadal hormones[62]. As this loss has to be replaced, this was included in the model by the addition of the rate constant k 25 which represents steroid hormone synthesis (PSS), which was in essence a sink species. ODEs 11 and 12 were then derived. Again at this point the model was re-parameterised and steady-states produced.

d P C E d t = k 23 A C A T P F C - k 24 PCE CEH ODE 11
d P F C d t = k 21 - k 25 + k 24 P C E C E H - k 23 A C A T P F C P C S max 1 + P F C / P F C t PCSS + k 20 L D L C [ P L D L R s (ODE 12
Table 7 Calculation of cholesterol in the peripheral tissue

Reverse cholesterol transport and the final steps to a whole-body mathematical model

Based on current literature it was assumed that the production of HDL-C is dependent on three factors; the concentration of peripheral free cholesterol, the population of nascent HDL particles and the activity of the enzyme LCAT which has a role to play in the conversion of free peripheral cholesterol into cholesteryl esters which are then sequestered into the core of a lipoprotein particle, eventually making the newly synthesized HDL particle. The enzyme is bound to HDLs and LDLs in the blood plasma.

Firstly, the population of nascent HDLs had to be represented. Nascent HDLs are synthesised in the intestine. As extensive literature searches revealed little quantitative information related to this process, an arbitrary value of 100 mg/day was assigned for the synthesis of nascent HDL. This value was then split evenly between its two points of origin; the intestine and liver. ODE 13 was then derived where k8 represents the intestinal rate of nascent HDL (NHDL) synthesis and k11 represents the hepatic synthesis of NHDL.

d N H D L d t = k 8 + k 11 ODE 13

Reverse cholesterol transport is important as it represents the only route for excess cholesterol generated in the peripheral tissue to be removed from the body, either by secretion into bile or by conversion into bile acids[55]. In the steady state, RCT should equal the rate of synthesis of cholesterol in the peripheral tissue[41]. Previously the synthesis of cholesterol in the peripheral tissue was calculated to equal 441 mg/day. Therefore, 441 mg of cholesterol that originated in peripheral tissue should enter the liver via HDL and this return of cholesterol via HDL-C to the liver may follow one of several routes. In the presence of CETP, a portion of the HDL cholesterol is transferred to either VLDL or LDL and ultimately returned to the liver via HLDLRs; therefore it was assumed that 10% of the cholesterol that is scavenged by HDL-C makes its way to the liver via this path. Alternatively, HDL-C is delivered to the liver via scavenger receptor class B1 [SRB1]. Thus based on this information the change in HDL-C with time was represented by ODE 14, where k 26 is the rate constant for the enzymatic dependent scavenging of cholesterol from the peripheral compartment, k 27 is the rate constant for the CETP dependent transfer of cholesterol to VLDL and k 28 is the rate constant for the CETP mediated transfer of cholesterol to LDL. Finally, k 29 is the rate constant for reverse cholesterol transport.

d HDLC dt = k 26 LCAT NHDL P F C - k 27 C E T P H D L C - k 28 C E T P H D L C - k 29 H D L C S R B 1 ODE 14

Updating of ODEs 5 and ODE6

ODE 5 and ODE 6 were updated to include the changes described above. Again the system was re-parameterized and brought into a steady state to reflect these new changes. HDL-C reached a steady state which was based on the literature for normolipidemic male. All reactions and parameters are summarised in Table2.

d H F C d t = k 19 + k 10 H C E C E H - B C R max 1 + B C R t / H F C BS H C S max 1 + H F C / H C S t HS + k 6 I B S I C + k 16 H L D L R s I D L C + k 18 H L D L R s L D L C + k 29 H D L C S R B 1 k 14 H L D L R s V L D L C - k 12 H F C - k 5 H F C B S P - k 9 A C A T H F C Updated ODE 5 d P F C d t = k 21 - k 25 + k 24 P C E C E H - k 26 L C A T N H D L P F C - k 23 A C A T P F C + P C S max 1 + P F C / P F C t PCSS + k 20 L D L C [ P L D L R s (Updated ODE 6

Model exchangeability, SBML and submission to the Biomodels database

The model was coded into SBML using an export function in MathSBML (Additional file2). Unlike computer languages such as C, C++ and Java, SBML was not designed to be coded manually; instead a number of tools are available that automatically generate SBML. The code can then be exchanged between tools supporting SBML[17]. In this case the model code was transferred into the software tool Copasi to facilitate parameter scanning[79]. Please note that when the ODEs are integrated in Copasi using the initial conditions, LDL-C enters a slightly higher steady-state than it does with MathSBML. In addition the model was submitted to the Biomodels database to facilitate its updating and future exchange (MODEL1206010000).


  1. 1.

    Barzilai N, Atzmon G, Derby CA, Bauman JM, Lipton RB: A genotype of exceptional longevity is associated with preservation of cognitive function. Neurology. 2006, 67 (12): 2170-2175. 10.1212/01.wnl.0000249116.50854.65.

    CAS  Article  Google Scholar 

  2. 2.

    Barzilai N, Atzmon G, Schechter C, Schaefer EJ, Cupples AL, Lipton R, Cheng S, Shuldiner AR: Unique lipoprotein phenotype and genotype associated with exceptional longevity. JAMA. 2003, 290 (15): 2030-2040. 10.1001/jama.290.15.2030.

    CAS  Article  Google Scholar 

  3. 3.

    Capewell S, Allender S, Critchley J, Lloyd-Williams F, O’Flaherty M, Rayner M, Scarborough P: Cardio & Vascular Coalition Modelling the UK burden of Cardiovascular Disease to 2020:A Research Report for the Cardio & Vascular Coalition and the British Heart Foundation. British Heart Foundation. 2008

    Google Scholar 

  4. 4.

    Castelli WP, Garrison RJ, Wilson PW, Abbott RD, Kalousdian S, Kannel WB: Incidence of coronary heart disease and lipoprotein cholesterol levels. The Framingham Study. JAMA. 1986, 256 (20): 2835-2838. 10.1001/jama.1986.03380200073024.

    CAS  Article  Google Scholar 

  5. 5.

    Berliner J, Leitinger N, Watson A, Huber J, Fogelman A, Navab M: Oxidized lipids in atherogenesis: formation, destruction and action. Thromb Haemost. 1997, 78 (1): 195-199.

    CAS  Google Scholar 

  6. 6.

    Heiss G, Tamir I, Davis CE, Tyroler HA, Rifkand BM, Schonfeld G, Jacobs D, Frantz ID: Lipoprotein-cholesterol distributions in selected North American populations: the lipid research clinics program prevalence study. Circulation. 1980, 61 (2): 302-315. 10.1161/01.CIR.61.2.302.

    CAS  Article  Google Scholar 

  7. 7.

    Abbott RD, Garrison RJ, Wilson PW, Epstein FH, Castelli WP, Feinleib M, LaRue C: Joint distribution of lipoprotein cholesterol classes. The Framingham study. Arteriosclerosis. 1983, 3 (3): 260-272. 10.1161/01.ATV.3.3.260.

    CAS  Article  Google Scholar 

  8. 8.

    Hollander D, Morgan D: Increase in cholesterol intestinal absorption with aging in the rat. Exp Gerontol. 1979, 14 (4): 201-204. 10.1016/0531-5565(79)90020-2.

    CAS  Article  Google Scholar 

  9. 9.

    Wang DQ: Aging per se is an independent risk factor for cholesterol gallstone formation in gallstone susceptible mice. J Lipid Res. 2002, 43 (11): 1950-1959. 10.1194/jlr.M200078-JLR200.

    CAS  Article  Google Scholar 

  10. 10.

    Grundy SM, Vega GL, Bilheimer DW: Kinetic mechanisms determining variability in low density lipoprotein levels and rise with age. Arteriosclerosis. 1985, 5 (6): 623-630. 10.1161/01.ATV.5.6.623.

    CAS  Article  Google Scholar 

  11. 11.

    Hood L: Systems biology: integrating technology, biology, and computation. Mech Ageing Dev. 2003, 124 (1): 9-16. 10.1016/S0047-6374(02)00164-1.

    Article  Google Scholar 

  12. 12.

    Wang K, Lee I, Carlson G, Hood L, Galas D: Systems biology and the discovery of diagnostic biomarkers. Dis Markers. 2010, 28 (4): 199-207.

    CAS  Article  Google Scholar 

  13. 13.

    Kirkwood TBL, Boys RJ, Gillespie CS, Proctor CJ, Shanley DP, Wilkinson DJ: (2005) Computer Modeling in the Study of Aging, in S. N. Austad, E. J. Masoro (eds.): Computer Modeling in the Study of Aging. Handbook of the Biology of Aging. Academic. 2005, 334-357.

    Chapter  Google Scholar 

  14. 14.

    Homan R, Krause BR: Established and emerging strategies for inhibition of cholesterol absorption. Curr PharmDes. 1997, 3: 29-44.

    CAS  Google Scholar 

  15. 15.

    Goodman DS, Smith FR, Seplowitz AH, Ramakrishnan R, Dell RB: Prediction of the parameters of whole body cholesterol metabolism in humans. J Lipid Res. 1980, 21 (6): 699-713.

    CAS  Google Scholar 

  16. 16.

    Kervizic G, Corcos L: Dynamical modeling of the cholesterol regulatory pathway with Boolean networks. BMC Syst Biol. 2008, 2: 99-10.1186/1752-0509-2-99.

    Article  Google Scholar 

  17. 17.

    Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, et al.: The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics. 2003, 19 (4): 524-531. 10.1093/bioinformatics/btg015.

    CAS  Article  Google Scholar 

  18. 18.

    Li C, Donizelli M, Rodriguez N, Dharuri H, Endler L, Chelliah V, Li L, He E, Henry A, Stefan MI, et al.: BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models. BMC Syst Biol. 2010, 4: 92-10.1186/1752-0509-4-92.

    Article  Google Scholar 

  19. 19.

    Havel RJ: Biology of cholesterol, lipoproteins and atherosclerosis. Clin Exp Hypertens A. 1989, 11 (5–6): 887-900.

    CAS  Google Scholar 

  20. 20.

    Brown AJ: Atherosclerosis: cell biology and lipoproteins: cholesterol absorption inhibitors: gateway therapy for hypercholesterolaemia. Curr Opin Lipidol. 2002, 13 (6): 701-703. 10.1097/00041433-200212000-00016.

    CAS  Article  Google Scholar 

  21. 21.

    Bosner MS, Lange LG, Stenson WF, Ostlund RE: Percent cholesterol absorption in normal women and men quantified with dual stable isotopic tracers and negative ion mass spectrometry. J Lipid Res. 1999, 40 (2): 302-308.

    CAS  Google Scholar 

  22. 22.

    Hui DY, Howles PN: Molecular mechanisms of cholesterol absorption and transport in the intestine. Semin Cell Dev Biol. 2005, 16 (2): 183-192. 10.1016/j.semcdb.2005.01.003.

    CAS  Article  Google Scholar 

  23. 23.

    Hernell O, Staggers JE, Carey MC: Physical-chemical behavior of dietary and biliary lipids during intestinal digestion and absorption. 2. Phase analysis and aggregation states of luminal lipids during duodenal fat digestion in healthy adult human beings. Biochemistry. 1990, 29 (8): 2041-2056. 10.1021/bi00460a012.

    CAS  Article  Google Scholar 

  24. 24.

    Westergaard H, Dietschy JM: The mechanism whereby bile acid micelles increase the rate of fatty acid and cholesterol uptake into the intestinal mucosal cell. J Clin Invest. 1976, 58 (1): 97-108. 10.1172/JCI108465.

    CAS  Article  Google Scholar 

  25. 25.

    Chijiiwa K, Linscheer WG: Distribution and monomer activity of cholesterol in micellar bile salt: effect of cholesterol level. Am J Physiol. 1987, 252 (3 Pt 1): G309-G314.

    CAS  Google Scholar 

  26. 26.

    Hauser H, Dyer JH, Nandy A, Vega MA, Werder M, Bieliauskaite E, Weber FE, Compassi S, Gemperli A, Boffelli D, et al.: Identification of a receptor mediating absorption of dietary cholesterol in the intestine. Biochemistry. 1998, 37 (51): 17843-17850. 10.1021/bi982404y.

    CAS  Article  Google Scholar 

  27. 27.

    Detmers PA, Patel S, Hernandez M, Montenegro J, Lisnock JM, Pikounis B, Steiner M, Kim D, Sparrow C, Chao YS, et al.: A target for cholesterol absorption inhibitors in the enterocyte brush border membrane. Biochim Biophys Acta. 2000, 1486 (2–3): 243-252.

    CAS  Article  Google Scholar 

  28. 28.

    Voshol PJ, Schwarz M, Rigotti A, Krieger M, Groen AK, Kuipers F: Down-regulation of intestinal scavenger receptor class B, type I (SR-BI) expression in rodents under conditions of deficient bile delivery to the intestine. Biochem J. 2001, 356 (Pt 2): 317-325.

    CAS  Article  Google Scholar 

  29. 29.

    Mahley RW, Innerarity TL: Lipoprotein receptors and cholesterol homeostasis. Biochim Biophys Acta. 1983, 737 (2): 197-222. 10.1016/0304-4157(83)90001-1.

    CAS  Article  Google Scholar 

  30. 30.

    Havel RJ: Postprandial lipid metabolism: an overview. Proc Nutr Soc. 1997, 56 (2): 659-666. 10.1079/PNS19970065.

    CAS  Article  Google Scholar 

  31. 31.

    Spady DK, Turley SD, Dietschy JM: Rates of low density lipoprotein uptake and cholesterol synthesis are regulated independently in the liver. J Lipid Res. 1985, 26 (4): 465-472.

    CAS  Google Scholar 

  32. 32.

    St-Pierre MV, Kullak-Ublick GA, Hagenbuch B, Meier PJ: Transport of bile acids in hepatic and non-hepatic tissues. J Exp Biol. 2001, 204 (Pt 10): 1673-1686.

    CAS  Google Scholar 

  33. 33.

    Hofmann AF: Bile Acids: The Good, the Bad, and the Ugly. News Physiol Sci. 1999, 14: 24-29.

    CAS  Google Scholar 

  34. 34.

    Olivecrona G, Olivecrona T: Triglyceride lipases and atherosclerosis. Curr Opin Lipidol. 1995, 6 (5): 291-305. 10.1097/00041433-199510000-00009.

    CAS  Article  Google Scholar 

  35. 35.

    Packard CJ, Shepherd J: Lipoprotein heterogeneity and apolipoprotein B metabolism. Arterioscler Thromb Vasc Biol. 1997, 17 (12): 3542-3556. 10.1161/01.ATV.17.12.3542.

    CAS  Article  Google Scholar 

  36. 36.

    Hussain MM, Strickland DK, Bakillah A: The mammalian low-density lipoprotein receptor family. Annu Rev Nutr. 1999, 19: 141-172. 10.1146/annurev.nutr.19.1.141.

    CAS  Article  Google Scholar 

  37. 37.

    Willnow TE, Nykjaer A, Herz J: Lipoprotein receptors: new roles for ancient proteins. Nat Cell Biol. 1999, 1 (6): E157-E162. 10.1038/14109.

    CAS  Article  Google Scholar 

  38. 38.

    Brown MS, Goldstein JL: Sterol regulatory element binding proteins (SREBPs): controllers of lipid synthesis and cellular uptake. Nutr Rev. 1998, 56 (2 Pt 2): S1-S3. discussion S54-75

    CAS  Google Scholar 

  39. 39.

    Small DM: Mechanisms of reversed cholesterol transport. Agents Actions Suppl. 1988, 26: 135-146.

    CAS  Google Scholar 

  40. 40.

    Fielding CJ, Fielding PE: Molecular physiology of reverse cholesterol transport. J Lipid Res. 1995, 36 (2): 211-228.

    CAS  Google Scholar 

  41. 41.

    Groen AK, Oude Elferink RP, Verkade HJ, Kuipers F: The ins and outs of reverse cholesterol transport. Ann Med. 2004, 36 (2): 135-145. 10.1080/07853890310020635.

    CAS  Article  Google Scholar 

  42. 42.

    Johnson WJ, Mahlberg FH, Rothblat GH, Phillips MC: Cholesterol transport between cells and high-density lipoproteins. Biochim Biophys Acta. 1991, 1085 (3): 273-298. 10.1016/0005-2760(91)90132-2.

    CAS  Article  Google Scholar 

  43. 43.

    Lewis GF, Rader DJ: New insights into the regulation of HDL metabolism and reverse cholesterol transport. Circ Res. 2005, 96 (12): 1221-1232. 10.1161/01.RES.0000170946.56981.5c.

    CAS  Article  Google Scholar 

  44. 44.

    Brown MS, Goldstein JL: A receptor-mediated pathway for cholesterol homeostasis. Science. 1986, 232 (4746): 34-47. 10.1126/science.3513311.

    CAS  Article  Google Scholar 

  45. 45.

    Brown MS, Goldstein JL: A proteolytic pathway that controls the cholesterol content of membranes, cells, and blood. Proc Natl Acad Sci U S A. 1999, 96 (20): 11041-11048. 10.1073/pnas.96.20.11041.

    CAS  Article  Google Scholar 

  46. 46.

    Miettinen TA, Gylling H: Synthesis and absorption markers of cholesterol in serum and lipoproteins during a large dose of statin treatment. Eur J Clin Invest. 2003, 33 (11): 976-982. 10.1046/j.1365-2362.2003.01229.x.

    CAS  Article  Google Scholar 

  47. 47.

    Watts GF, Chan DC, Barrett PH, O'Neill FH, Thompson GR: Effect of a statin on hepatic apolipoprotein B-100 secretion and plasma campesterol levels in the metabolic syndrome. Int J Obes Relat Metab Disord. 2003, 27 (7): 862-865. 10.1038/sj.ijo.0802287.

    CAS  Article  Google Scholar 

  48. 48.

    Grundy SM, Vega GL, Bilheimer DW: Influence of combined therapy with mevinolin and interruption of bile-acid reabsorption on low density lipoproteins in heterozygous familial hypercholesterolemia. Ann Intern Med. 1985, 103 (3): 339-343.

    CAS  Article  Google Scholar 

  49. 49.

    Ericsson S, Eriksson M, Vitols S, Einarsson K, Berglund L, Angelin B: Influence of age on the metabolism of plasma low density lipoproteins in healthy males. J Clin Invest. 1991, 87 (2): 591-596. 10.1172/JCI115034.

    CAS  Article  Google Scholar 

  50. 50.

    Millar JS, Lichtenstein AH, Cuchel M, Dolnikowski GG, Hachey DL, Cohn JS, Schaefer EJ: Impact of age on the metabolism of VLDL, IDL, and LDL apolipoprotein B-100 in men. J Lipid Res. 1995, 36 (6): 1155-1167.

    CAS  Google Scholar 

  51. 51.

    Uchida K, Chikai T, Takase H, Nomura Y, Seo S, Nakao H, Takeuchi N: Age-related changes of bile acid metabolism in rats. Arch Gerontol Geriatr. 1990, 10 (1): 37-48. 10.1016/0167-4943(90)90042-5.

    CAS  Article  Google Scholar 

  52. 52.

    McNamara DJ: The impact of egg limitations on coronary heart disease risk: do the numbers add up?. J Am Coll Nutr. 2000, 19 (5 Suppl): 540S-548S.

    CAS  Article  Google Scholar 

  53. 53.

    Knoblauch H, Schuster H, Luft FC, Reich J: A pathway model of lipid metabolism to predict the effect of genetic variability on lipid levels. J Mol Med (Berl). 2000, 78 (9): 507-515. 10.1007/s001090000156.

    CAS  Article  Google Scholar 

  54. 54.

    Chun PW, Espinosa AJ, Lee CW, Shireman RB, Brumbaugh EE: Low density lipoprotein receptor regulation. Kinetic models. Biophys Chem. 1985, 21 (3–4): 185-209.

    CAS  Article  Google Scholar 

  55. 55.

    Wattis JA, O'Malley B, Blackburn H, Pickersgill L, Panovska J, Byrne HM, Jackson KG: Mathematical model for low density lipoprotein (LDL) endocytosis by hepatocytes. Bull Math Biol. 2008, 70 (8): 2303-2333. 10.1007/s11538-008-9347-9.

    CAS  Article  Google Scholar 

  56. 56.

    Adiels M, Packard C, Caslake MJ, Stewart P, Soro A, Westerbacka J, Wennberg B, Olofsson SO, Taskinen MR, Boren J: A new combined multicompartmental model for apolipoprotein B-100 and triglyceride metabolism in VLDL subfractions. J Lipid Res. 2005, 46 (1): 58-67.

    CAS  Article  Google Scholar 

  57. 57.

    Erkkila AT, Lichtenstein AH: Fiber and cardiovascular disease risk: how strong is the evidence?. J Cardiovasc Nurs. 2006, 21 (1): 3-8.

    Article  Google Scholar 

  58. 58.

    Brown L, Rosner B, Willett WW, Sacks FM: Cholesterol-lowering effects of dietary fiber: a meta-analysis. Am J Clin Nutr. 1999, 69 (1): 30-42.

    CAS  Google Scholar 

  59. 59.

    Lees AM, Mok HY, Lees RS, McCluskey MA, Grundy SM: Plant sterols as cholesterol-lowering agents: clinical trials in patients with hypercholesterolemia and studies of sterol balance. Atherosclerosis. 1977, 28 (3): 325-338. 10.1016/0021-9150(77)90180-0.

    CAS  Article  Google Scholar 

  60. 60.

    Mattson FH, Volpenhein RA, Erickson BA: Effect of plant sterol esters on the absorption of dietary cholesterol. J Nutr. 1977, 107 (7): 1139-1146.

    CAS  Google Scholar 

  61. 61.

    Law M: Plant sterol and stanol margarines and health. BMJ. 2000, 320 (7238): 861-864. 10.1136/bmj.320.7238.861.

    CAS  Article  Google Scholar 

  62. 62.

    Katan MB, Grundy SM, Jones P, Law M, Miettinen T, Paoletti R: Efficacy and safety of plant stanols and sterols in the management of blood cholesterol levels. Mayo Clin Proc. 2003, 78 (8): 965-978.

    CAS  Article  Google Scholar 

  63. 63.

    St-Onge MP, Jones PJ: Phytosterols and human lipid metabolism: efficacy, safety, and novel foods. Lipids. 2003, 38 (4): 367-375. 10.1007/s11745-003-1071-3.

    CAS  Article  Google Scholar 

  64. 64.

    Scheer M, Grote A, Chang A, Schomburg I, Munaretto C, Rother M, Sohngen C, Stelzer M, Thiele J, Schomburg D: BRENDA, the enzyme information system in 2011. Nucleic Acids Res. 2011, 39 (Database issue): D670-D676.

    CAS  Article  Google Scholar 

  65. 65.

    Cobbold CA, Sherratt JA, Maxwell SR: Lipoprotein oxidation and its significance for atherosclerosis: a mathematical approach. Bull Math Biol. 2002, 64 (1): 65-95. 10.1006/bulm.2001.0267.

    CAS  Article  Google Scholar 

  66. 66.

    Krause F, Uhlendorf J, Lubitz T, Schulz M, Klipp E, Liebermeister W: Annotation and merging of SBML models with semanticSBML. Bioinformatics. 2010, 26 (3): 421-422. 10.1093/bioinformatics/btp642.

    CAS  Article  Google Scholar 

  67. 67.

    Schulz M, Uhlendorf J, Klipp E, Liebermeister W: SBMLmerge, a system for combining biochemical network models. Genome Inform. 2006, 17 (1): 62-71.

    CAS  Google Scholar 

  68. 68.

    Witztum JL, Steinberg D: Role of oxidized low density lipoprotein in atherogenesis. J Clin Invest. 1991, 88 (6): 1785-1792. 10.1172/JCI115499.

    CAS  Article  Google Scholar 

  69. 69.

    Steinberg D: Low density lipoprotein oxidation and its pathobiological significance. J Biol Chem. 1997, 272 (34): 20963-20966. 10.1074/jbc.272.34.20963.

    CAS  Article  Google Scholar 

  70. 70.

    Heijmans BT, Beekman M, Houwing-Duistermaat JJ, Cobain MR, Powell J, Blauw GJ, van der Ouderaa F, Westendorp RG, Slagboom PE: Lipoprotein particle profiles mark familial and sporadic human longevity. PLoS Med. 2006, 3 (12): e495-10.1371/journal.pmed.0030495.

    Article  Google Scholar 

  71. 71.

    Rahilly-Tierney CR, Spiro A, Vokonas P, Gaziano JM: Relation between high-density lipoprotein cholesterol and survival to age 85 years in men (from the VA normative aging study). Am J Cardiol. 2011, 107 (8): 1173-1177. 10.1016/j.amjcard.2010.12.015.

    CAS  Article  Google Scholar 

  72. 72.

    Goldstein B, Wofsy C, Bell G: Interactions of low density lipoprotein receptors with coated pits on human fibroblasts: estimate of the forward rate constant and comparison with the diffusion limit. Proc Natl Acad Sci U S A. 1981, 78 (9): 5695-5698. 10.1073/pnas.78.9.5695.

    CAS  Article  Google Scholar 

  73. 73.

    Keizer J, Ramirez J, Peacock-Lopez E: The effect of diffusion on the binding of membrane-bound receptors to coated pits. Biophys J. 1985, 47 (1): 79-87. 10.1016/S0006-3495(85)83879-0.

    CAS  Article  Google Scholar 

  74. 74.

    Gao H, Shi W, Freund LB: Mechanics of receptor-mediated endocytosis. Proc Natl Acad Sci U S A. 2005, 102 (27): 9469-9474. 10.1073/pnas.0503879102.

    CAS  Article  Google Scholar 

  75. 75.

    Pont F, Duvillard L, Verges B, Gambert P: Development of compartmental models in stable-isotope experiments: application to lipid metabolism. Arterioscler Thromb Vasc Biol. 1998, 18 (6): 853-860. 10.1161/01.ATV.18.6.853.

    CAS  Article  Google Scholar 

  76. 76.

    Le Novere N, Hucka M, Mi H, Moodie S, Schreiber F, Sorokin A, Demir E, Wegner K, Aladjem MI, Wimalaratne SM, et al.: The systems biology graphical notation. Nat Biotechnol. 2009, 27 (8): 735-741. 10.1038/nbt.1558.

    CAS  Article  Google Scholar 

  77. 77.

    Shapiro BE, Hucka M, Finney A, Doyle J: MathSBML: a package for manipulating SBML-based biological models. Bioinformatics. 2004, 20 (16): 2829-2831. 10.1093/bioinformatics/bth271.

    CAS  Article  Google Scholar 

  78. 78.

    Wolfram Research I: Mathematica Edition: Version 5.2. 2005, Champaign, Illinois: Wolfram Research Inc

    Google Scholar 

  79. 79.

    Mendes P, Hoops S, Sahle S, Gauges R, Dada J, Kummer U: Computational modeling of biochemical networks using COPASI. Methods Mol Biol. 2009, 500: 17-59. 10.1007/978-1-59745-525-1_2.

    CAS  Article  Google Scholar 

  80. 80.

    Henderson LG, Gregory J, Irving K, Swan G: National diet and Nutrtion Survey (NDNS) of adults aged 19 to 64. 2003, 70-2

    Google Scholar 

  81. 81.

    Dietschy JM, Turley SD, Spady DK: Role of liver in the maintenance of cholesterol and low density lipoprotein homeostasis in different animal species, including humans. J Lipid Res. 1993, 34 (10): 1637-1659.

    CAS  Google Scholar 

  82. 82.

    Dietschy JM: Theoretical considerations of what regulates low-density-lipoprotein and high-density-lipoprotein cholesterol. Am J Clin Nutr. 1997, 65 (5 Suppl): 1581S-1589S.

    CAS  Google Scholar 

  83. 83.

    Grundy SM, Metzger AL: A physiological method for estimation of hepatic secretion of biliary lipids in man. Gastroenterology. 1972, 62 (6): 1200-1217.

    CAS  Google Scholar 

  84. 84.

    Soars MG, Burchell B, Riley RJ: In vitro analysis of human drug glucuronidation and prediction of in vivo metabolic clearance. J Pharmacol Exp Ther. 2002, 301 (1): 382-390. 10.1124/jpet.301.1.382.

    CAS  Article  Google Scholar 

  85. 85.

    Redinger RN: The coming of age of our understanding of the enterohepatic circulation of bile salts. Am J Surg. 2003, 185 (2): 168-172. 10.1016/S0002-9610(02)01212-6.

    CAS  Article  Google Scholar 

  86. 86.

    Sudhop T, Lutjohann D, Kodal A, Igel M, Tribble DL, Shah S, Perevozskaya I, von Bergmann K: Inhibition of intestinal cholesterol absorption by ezetimibe in humans. Circulation. 2002, 106 (15): 1943-1948. 10.1161/01.CIR.0000034044.95911.DC.

    CAS  Article  Google Scholar 

  87. 87.

    Sahlin S, Granstrom L, Gustafsson U, Stahlberg D, Backman L, Einarsson K: Hepatic esterification rate of cholesterol and biliary lipids in human obesity. J Lipid Res. 1994, 35 (3): 484-490.

    CAS  Google Scholar 

  88. 88.

    Cook R: Distribution of sterols in organisms and in tissues. Chemistry, biochemistry, and pathology. 1958, New York: Academic Press Inc

    Google Scholar 

  89. 89.

    Erickson SK, Bosterling B: Cholesterol 7 alpha-hydroxylase from human liver: partial purification and reconstruction into defined phospholipid-cholesterol vesicles. J Lipid Res. 1981, 22 (5): 872-876.

    CAS  Google Scholar 

  90. 90.

    Glomset JA, Norum KR: The metabolic role of lecithin: cholesterol acyltransferase: perspectives form pathology. Adv Lipid Res. 1973, 11: 1-65.

    CAS  Google Scholar 

  91. 91.

    Brown MS, Goldstein JL: Receptor-mediated control of cholesterol metabolism. Science. 1976, 191 (4223): 150-154. 10.1126/science.174194.

    CAS  Article  Google Scholar 

  92. 92.

    August E, Parker KH, Barahona M: A dynamical model of lipoprotein metabolism. Bull Math Biol. 2007, 69 (4): 1233-1254. 10.1007/s11538-006-9132-6.

    CAS  Article  Google Scholar 

  93. 93.

    Yamada N, Shames DM, Stoudemire JB, Havel RJ: Metabolism of lipoproteins containing apolipoprotein B-100 in blood plasma of rabbits: heterogeneity related to the presence of apolipoprotein E. Proc Natl Acad Sci U S A. 1986, 83 (10): 3479-3483. 10.1073/pnas.83.10.3479.

    CAS  Article  Google Scholar 

  94. 94.

    Spady DK, Bilheimer DW, Dietschy JM: Rates of receptor-dependent and -independent low density lipoprotein uptake in the hamster. Proc Natl Acad Sci U S A. 1983, 80 (11): 3499-3503. 10.1073/pnas.80.11.3499.

    CAS  Article  Google Scholar 

  95. 95.

    Spady DK, Meddings JB, Dietschy JM: Kinetic constants for receptor-dependent and receptor-independent low density lipoprotein transport in the tissues of the rat and hamster. J Clin Invest. 1986, 77 (5): 1474-1481. 10.1172/JCI112460.

    CAS  Article  Google Scholar 

  96. 96.

    Spady DK, Huettinger M, Bilheimer DW, Dietschy JM: Role of receptor-independent low density lipoprotein transport in the maintenance of tissue cholesterol balance in the normal and WHHL rabbit. J Lipid Res. 1987, 28 (1): 32-41.

    CAS  Google Scholar 

Download references


We acknowledge Fotios Drenos, Daryl Shanley and Carole Proctor for their advice throughout the duration of this work. We also acknowledge Chris Packard for his advice on cholesterol metabolism. MT Mc Auley was funded by an EPSRC CASE studentship with Unilever PLC. The research was also supported by the National Institute for Health Research (NIHR) Newcastle Biomedical Research Centre based at Newcastle upon Tyne Hospitals NHS Foundation Trust and Newcastle University.

Author information



Corresponding author

Correspondence to Mark T Mc Auley.

Additional information

Competing interests

The authors declared that they have no competing interests.

Authors' contributions

MTMcA was responsible for the computational mathematical modeling and sensitivity analysis of the system. MTMcA was also responsible for data mining and acquisition of variable and parameter data. TBLK provided insight into the effects of ageing on the system, while DJW and JLJ advised on the modeling of the system. All authors read and approved the final manuscript.

Mark T Mc Auley, Darren J Wilkinson, Janette JL Jones and Thomas BL Kirkwood contributed equally to this work.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Mc Auley, M.T., Wilkinson, D.J., Jones, J.J. et al. A whole-body mathematical model of cholesterol metabolism and its age-associated dysregulation. BMC Syst Biol 6, 130 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Cholesterol
  • Bile Salt
  • High Density Lipoprotein
  • Cholesterol Metabolism
  • System Biology Markup Language