An adaptable in silico ensemble model of the arachidonic acid cascade.

Eicosanoids are a family of bioactive lipids, including derivatives of the ubiquitous fatty acid arachidonic acid (AA). The intimate involvement of eicosanoids in inflammation motivates the development of predictive in silico models for a systems-level exploration of disease mechanisms, drug development and replacement of animal models. Using an ensemble modelling strategy, we developed a computational model of the AA cascade. This approach allows the visualisation of plausible and thermodynamically feasible predictions, overcoming the limitations of fixed-parameter modelling. A quality scoring method was developed to quantify the accuracy of ensemble predictions relative to experimental data, measuring the overall uncertainty of the process. Monte Carlo ensemble modelling was used to quantify the prediction confidence levels. Model applicability was demonstrated using mass spectrometry mediator lipidomics to measure eicosanoids produced by HaCaT epidermal keratinocytes and 46BR.1N dermal fibroblasts, treated with stimuli (calcium ionophore A23187), (ultraviolet radiation, adenosine triphosphate) and a cyclooxygenase inhibitor (indomethacin). Experimentation and predictions were in good qualitative agreement, demonstrating the ability of the model to be adapted to cell types exhibiting differences in AA release and enzyme concentration profiles. The quantitative agreement between experimental and predicted outputs could be improved by expanding network topology to include additional reactions. Overall, our approach generated an adaptable, tuneable ensemble model of the AA cascade that can be tailored to represent different cell types and demonstrated that the integration of in silico and in vitro methods can facilitate a greater understanding of complex biological networks such as the AA cascade.

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

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Identification of McbR as transcription regulator of aecD and genes involved in methionine and cysteine biosynthesis in Corynebacterium jeikeium K411.

The C-S lyase aecD (MetC) from skin corynebacteria plays an important role in body odour formation by releasing odoriferous sulfanylalkanols from cysteine conjugates in human axilla secretions. The expression of the aecD gene from Corynebacterium jeikeium K411, a strain originally isolated from the human axilla, was down-regulated in cells grown in minimal medium supplemented with methionine. A candidate transcription regulator binding in front of the aecD coding region was detected by DNA affinity chromatography and identified as McbR by peptide mass fingerprinting. A 16-bp McbR-binding site was localized in the mapped promoter region of the aecD gene. The binding of purified McbR protein to the 16-bp sequence motif was demonstrated by DNA band shift assays. Comparative DNA microarray hybridizations and bioinformatic motif searches revealed the gene composition of the McbR regulon from C. jeikeium, including 28 genes that are organized in 16 transcription units. The McbR protein from C. jeikeium K411 directly regulates genes involved in methionine uptake and biosynthesis, in cysteine biosynthesis and sulfate reduction, and in the biosynthesis of amino acids belonging to the aspartate family.

A mathematical model of aging-related and cortisol induced hippocampal dysfunction.

BACKGROUND: The hippocampus is essential for declarative memory synthesis and is a core pathological substrate for Alzheimer's disease (AD), the most common aging-related dementing disease. Acute increases in plasma cortisol are associated with transient hippocampal inhibition and retrograde amnesia, while chronic cortisol elevation is associated with hippocampal atrophy. Thus, cortisol levels could be monitored and managed in older people, to decrease their risk of AD type hippocampal dysfunction. We generated an in silicomodel of the chronic effects of elevated plasma cortisol on hippocampal activity and atrophy, using the systems biology mark-up language (SBML). We further challenged the model with biologically based interventions to ascertain if cortisol associated hippocampal dysfunction could be abrogated. RESULTS: The in silicoSBML model reflected the in vivoaging of the hippocampus and increased plasma cortisol and negative feedback to the hypothalamic pituitary axis. Aging induced a 12% decrease in hippocampus activity (HA), increased to 30% by acute and 40% by chronic elevations in cortisol. The biological intervention attenuated the cortisol associated decrease in HA by 2% in the acute cortisol simulation and by 8% in the chronic simulation. CONCLUSION: Both acute and chronic elevations in cortisol secretion increased aging-associated hippocampal atrophy and a loss of HA in the model. We suggest that this first SMBL model, in tandem with in vitroand in vivostudies, may provide a backbone to further frame computational cortisol and brain aging models, which may help predict aging-related brain changes in vulnerable older people.

Duality, thermodynamics, and the linear programming problem in constraint-based models of metabolism.

It is shown that the dual to the linear programming problem that arises in constraint-based models of metabolism can be given a thermodynamic interpretation in which the shadow prices are chemical potential analogues, and the objective is to minimize free energy consumption given a free energy drain corresponding to growth. The interpretation is distinct from conventional nonequilibrium thermodynamics, although it does satisfy a minimum entropy production principle. It can be used to motivate extensions of constraint-based modeling, for example, to microbial ecosystems.