Restricted access
Articles
ArticlesRobust modelling, measurement and analysis of human and animal metabolic systems
Published:28 May 2009https://doi.org/10.1098/rsta.2008.0305
Abstract
Modelling human and animal metabolism is impeded by the lack of accurate quantitative parameters and the large number of biochemical reactions. This problem may be tackled by: (i) study of modules of the network independently; (ii) ensemble simulations to explore many plausible parameter combinations; (iii) analysis of ‘sloppy’ parameter behaviour, revealing interdependent parameter combinations with little influence; (iv) multiscale analysis that combines molecular and whole network data; and (v) measuring metabolic flux (rate of flow) in vivo via stable isotope labelling. For the latter method, carbon transition networks were modelled with systems of ordinary differential equations, but we show that coloured Petri nets provide a more intuitive graphical approach. Analysis of parameter sensitivities shows that only a few parameter combinations have a large effect on predictions. Model analysis of high-energy phosphate transport indicates that membrane permeability, inaccurately known at the organellar level, can be well determined from whole-organ responses. Ensemble simulations that take into account the imprecision of measured molecular parameters contradict the popular hypothesis that high-energy phosphate transport in heart muscle is mostly by phosphocreatine. Combining modular, multiscale, ensemble and sloppy modelling approaches with in vivo flux measurements may prove indispensable for the modelling of the large human metabolic system.
References
Antoniewicz M., Kelleher J.& Stephanopoulos G. . 2007Elementary metabolite units (EMU): a novel framework for modeling isotopic distributions. Metab. Eng. 9, 68–86.doi:10.1016/j.ymben.2006.09.001. . Crossref, PubMed, ISI, Google ScholarBessman S.P.& Geiger P.J. . 1981Transport of energy in muscle: the phosphorylcreatine shuttle. Science. 211, 448–452.doi:10.1126/science.6450446. . Crossref, PubMed, ISI, Google ScholarBinsl T.W., Mullen K.M., van Stokkum I.H.M., Heringa J.& van Beek J.H.G.M. . 2007FluxSimulator: an R package to simulate isotopomer distributions in metabolic networks. J. Stat. Softw. 18, 1–17. Crossref, ISI, Google ScholarBrown K.S., Hill C.C., Calero G.A., Myers C.R., Lee K.H., Sethna J.P.& Cerione R.A. . 2004The statistical mechanics of complex signaling networks: nerve growth factor signaling. Phys. Biol. 1, 184–195.doi:10.1088/1478-3967/1/3/006. . Crossref, PubMed, ISI, Google Scholar- Chouikha, M. & Schnieder, E. 1998 Modelling continuous–discrete systems with hybrid Petri nets. In IEEE-MC Multiconf. Comput. Eng. Syst. Appl., pp. 606–612. Google Scholar
Duarte N.C., Becker S.A., Jamshidi N., Thiele I., Mo M.L., Vo T.D., Srivas R.& Palsson B.Ø. . 2007Global reconstruction of the human metabolic network based on genomic and bibliomic data. Proc. Natl Acad. Sci. USA. 104, 1777–1782.doi:10.1073/pnas.0610772104. . Crossref, PubMed, ISI, Google ScholarFunahashi A., Tanimura N., Morohashi M.& Kitano H. . 2003CellDesigner: a process diagram editor for gene-regulatory and biochemical networks. BIOSILICO. 1, 159–162.doi:10.1016/S1478-5382(03)02370-9. . Crossref, Google ScholarGokulakrishnan P., Lawrence A.D., McLellan P.J.& Grandmaison E.W. . 2006A functional-PCA approach for analyzing and reducing complex chemical mechanisms. Comput. Chem. Eng. 30, 1093–1101.doi:10.1016/j.compchemeng.2006.02.007. . Crossref, ISI, Google Scholar- Gutenkunst, R. N., Atlas, J. C., Casey, F. P., Kuczenski, R. S., Waterfall, J. J., Myers, C. R. & Sethna, J. 2007aSloppyCell documentation. See http://sloppycell.sourceforge.net/. Google Scholar
Gutenkunst R.N., Waterfall J., Casey F., Brown K., Myers C.& Sethna J. Universally sloppy parameter sensitivities in systems biology models. PLoS Comput. Biol. 3, 2007b1–10.doi:10.1371/journal.pcbi.0030189. . Crossref, ISI, Google Scholar- Jensen, K. 1998 A brief introduction to coloured Petri nets. In Proc. Workshop on the Applicability of Formal Models, pp. 55–58. Google Scholar
Kelleher J. . 2001Flux estimation using isotopic tracers: common ground for metabolic physiology and metabolic engineering. Metab. Eng. 3, 100–110.doi:10.1006/mben.2001.0185. . Crossref, PubMed, ISI, Google Scholar- Kongas, O. & van Beek, J. H. G. M. 2007 Creatine kinase in energy metabolic signaling in muscle. Nat. Precedings. (doi:10.1038/npre.2007.1317.1) (Publication accompanying computational model published in the BioModels database. ID: BIOMD0000000041-Kongas2001_creatine; see http://www.ebi.ac.uk/compneur-srv/biomodels-main/publ-model.do?mid=BIOMD0000000041. See also in the CellML model repository http://www.cellml.org/models/kongas_vanbeek_2001_version01). Google Scholar
Le Novère N., 2006BioModels database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems. Nucleic Acids Res. 34, Database issue, D689–D691.doi:10.1093/nar/gkj092. . Crossref, PubMed, ISI, Google ScholarMa H., Sorokin A., Mazein A., Selkov A., Selkov E., Demin O.& Goryanin I. . 2007The Edinburgh human metabolic network reconstruction and its functional analysis. Mol. Syst. Biol. 3, 135doi:10.1038/msb4100177. . Crossref, PubMed, ISI, Google ScholarMurata T. . 1989Petri nets: properties, analysis and applications. Proc. IEEE. 77, 541–580.doi:10.1109/5.24143. . Crossref, ISI, Google ScholarNelder J.A.& Mead R. . 1965A simplex method for function minimization. Comput. J. 17, 308–313. Crossref, ISI, Google Scholar- R Development Core Team 2008 R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. See http://www.R-project.org. Google Scholar
Sherry A., Jeffrey F.& Malloy C. . 2004Analytical solutions for 13C isotopomer analysis of complex metabolic conditions: substrate oxidation, multiple pyruvate cycles, and gluconeogenis. Metab. Eng. 6, 12–24.doi:10.1016/j.ymben.2003.10.007. . Crossref, PubMed, ISI, Google ScholarTeusink B., 2000Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry. Eur. J. Biochem. 267, 5313–5329.doi:10.1046/j.1432-1327.2000.01527.x. . Crossref, PubMed, Google Scholarvan Beek J.H.G.M. . 2007Adenine nucleotide–creatine–phosphate module in myocardial metabolic system explains fast phase of dynamic regulation of oxidative phosphorylation. Am. J. Physiol. Cell Physiol. 293, C815–C829.doi:10.1152/ajpcell.00355.2006. . Crossref, PubMed, ISI, Google Scholarvan Beek J.H.G.M. . 2008Multiscale and modular analysis of cardiac energy metabolism. Repairing the broken interfaces of isolated system components. Ann. NY Acad. Sci. 1123, 155–168.doi:10.1196/annals.1420.018. . Crossref, PubMed, ISI, Google Scholarvan Beek J.H.G.M., Csont T., de Kanter F.& Bussemaker J. . 1998Simple model analysis of 13C NMR spectra to measure oxygen consumption using frozen tissue samples. Adv. Exp. Med. Biol. 454, 475–485. Crossref, PubMed, ISI, Google Scholarvan Beek J.H.G.M., van Mil H., King R., de Kanter F., Alders D.& Bussemaker J. . 1999A 13C NMR double-labeling method to quantitate local myocardial O2 consumption using frozen tissue samples. Am. J. Physiol. Heart Circ. Physiol. 277, H1630–H1640. Crossref, ISI, Google ScholarWiechert W. . 200113C metabolic flux analysis. Metab. Eng. 3, 195–206.doi:10.1006/mben.2001.0187. . Crossref, PubMed, ISI, Google ScholarWishart D.S., 2007HMDB: the Human Metabolome Database. Nucleic Acids Res. 35, Database issue, D521–D526.doi:10.1093/nar/gkl923. . Crossref, PubMed, ISI, Google Scholar
