A method for estimation of elasticities in metabolic networks using steady state and dynamic metabolomics data and linlog kinetics

BackgroundDynamic modeling of metabolic reaction networks under in vivo conditions is a crucial step in order to obtain a better understanding of the (dis)functioning of living cells. So far dynamic metabolic models generally have been based on mechanistic rate equations which often contain so many parameters that their identifiability from experimental data forms a serious problem. Recently, approximative rate equations, based on the linear logarithmic (linlog) format have been proposed as a suitable alternative with fewer parameters.ResultsIn this paper we present a method for estimation of the kinetic model parameters, which are equal to the elasticities defined in Metabolic Control Analysis, from metabolite data obtained from dynamic as well as steady state perturbations, using the linlog kinetic format. Additionally, we address the question of parameter identifiability from dynamic perturbation data in the presence of noise. The method is illustrated using metabolite data generated with a dynamic model of the glycolytic pathway of Saccharomyces cerevisiae based on mechanistic rate equations. Elasticities are estimated from the generated data, which define the complete linlog kinetic model of the glycolysis. The effect of data noise on the accuracy of the estimated elasticities is presented. Finally, identifiable subset of parameters is determined using information on the standard deviations of the estimated elasticities through Monte Carlo (MC) simulations.ConclusionThe parameter estimation within the linlog kinetic framework as presented here allows the determination of the elasticities directly from experimental data from typical dynamic and/or steady state experiments. These elasticities allow the reconstruction of the full kinetic model of Saccharomyces cerevisiae, and the determination of the control coefficients. MC simulations revealed that certain elasticities are potentially unidentifiable from dynamic data only. Addition of steady state perturbation of enzyme activities solved this problem.

[1]  B. Palsson,et al.  Metabolic capabilities of Escherichia coli: I. synthesis of biosynthetic precursors and cofactors. , 1993, Journal of theoretical biology.

[2]  J. Heijnen,et al.  A new framework for the estimation of control parameters in metabolic pathways using lin-log kinetics. , 2004, European journal of biochemistry.

[3]  J. Heijnen,et al.  Determination of elasticities, concentration and flux control coefficients from transient metabolite data using linlog kinetics. , 2005, Metabolic engineering.

[4]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[5]  M A Aon,et al.  Fluxes of carbon, phosphorylation, and redox intermediates during growth of saccharomyces cerevisiae on different carbon sources , 1995, Biotechnology and bioengineering.

[6]  M. Reuss,et al.  In vivo analysis of metabolic dynamics in Saccharomyces cerevisiae : I. Experimental observations. , 1997, Biotechnology and bioengineering.

[7]  J. Nielsen,et al.  In vivo dynamics of galactose metabolism in Saccharomyces cerevisiae: metabolic fluxes and metabolite levels. , 2001, Biotechnology and bioengineering.

[8]  M C Mossing,et al.  Variability of the intracellular ionic environment of Escherichia coli. Differences between in vitro and in vivo effects of ion concentrations on protein-DNA interactions and gene expression. , 1987, The Journal of biological chemistry.

[9]  M. Reuss,et al.  In VivoDynamics of the Pentose Phosphate Pathway inSaccharomyces cerevisiae , 1999 .

[10]  L. Rédei.  CHAPTER I – SET-THEORETICAL PRELIMINARIES , 1967 .

[11]  D Weuster-Botz,et al.  Automated sampling device for monitoring intracellular metabolite dynamics. , 1999, Analytical biochemistry.

[12]  J. Heijnen,et al.  Analysis of in vivo kinetics of glycolysis in aerobic Saccharomyces cerevisiae by application of glucose and ethanol pulses , 2004, Biotechnology and bioengineering.

[13]  I. Birol,et al.  Metabolic control analysis under uncertainty: framework development and case studies. , 2004, Biophysical journal.

[14]  M. Savageau.  Biochemical Systems Analysis: A Study of Function and Design in Molecular Biology , 1976 .

[15]  J. Bailey,et al.  Fermentation pathway kinetics and metabolic flux control in suspended and immobilized Saccharomyces cerevisiae , 1990 .

[16]  James E. Bailey,et al.  Ethanol Production in Baker's Yeast: Application of Experimental Perturbation Techniques for Model Development and Resultant Changes in Flux Control Analysis , 1994 .

[17]  Antje Chang,et al.  BRENDA, enzyme data and metabolic information , 2002, Nucleic Acids Res..

[18]  J. Heijnen,et al.  A metabolic network stoichiometry analysis of microbial growth and product formation , 1995, Biotechnology and bioengineering.

[19]  V. Hatzimanikatis,et al.  Metabolic engineering under uncertainty--II: analysis of yeast metabolism. , 2006, Metabolic engineering.

[20]  G. Church,et al.  Analysis of optimality in natural and perturbed metabolic networks , 2002 .

[21]  J. Heijnen,et al.  Rapid sampling for analysis of in vivo kinetics using the BioScope: a system for continuous-pulse experiments. , 2002, Biotechnology and bioengineering.

[22]  J. Heijnen,et al.  The mathematics of metabolic control analysis revisited. , 2002, Metabolic engineering.

[23]  G Zacchi,et al.  Influence of experimental errors on the determination of flux control coefficients from transient metabolite concentrations. , 1996, The Biochemical journal.

[24]  G. T. Tsao,et al.  Are Microbes Optimal Strategists? , 1987 .

[25]  Sonia Cortassa,et al.  Metabolic control analysis of glycolysis and branching to ethanol production in chemostat cultures of Saccharomyces cerevisiae under carbon, nitrogen, or phosphate limitations , 1994 .

[26]  Mark M. Meerschaert,et al.  Mathematical Modeling , 2014, Encyclopedia of Social Network Analysis and Mining.

[27]  M. Reuss,et al.  In vivo dynamics of the pentose phosphate pathway in Saccharomyces cerevisiae. , 1999, Metabolic engineering.

[28]  W. Knorre,et al.  M. A. Savageau, Biochemical Systems Analysis. A Study of Function and Design in Molecular Biology. 396 S., 115 Abb., 14 Tab. Reading, Mass. 1976. Addison‐Wesley Pbl. Co./Advanced Book Program. £ 26,50 , 1979 .

[29]  Matthias Reuss,et al.  Optimal re-design of primary metabolism in Escherichia coli using linlog kinetics. , 2004, Metabolic engineering.

[30]  A. Cornish-Bowden,et al.  Extending double modulation: combinatorial rules for identifying the modulations necessary for determining elasticities in metabolic pathways. , 1996, Journal of theoretical biology.

[31]  M. Reuss,et al.  In vivo analysis of metabolic dynamics in Saccharomyces cerevisiae: II. Mathematical model. , 1997, Biotechnology and bioengineering.

[32]  Barbara M. Bakker,et al.  Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry. , 2000, European journal of biochemistry.

[33]  L. Wackett.  An annotated selection of World Wide Web sites relevant to the topics in Microbial Biotechnology , 2013, Microbial biotechnology.

[34]  Vassily Hatzimanikatis,et al.  Metabolic engineering under uncertainty. I: framework development. , 2006, Metabolic engineering.

[35]  J. Bailey,et al.  Toward a science of metabolic engineering , 1991, Science.

[36]  A. Cornish-Bowden,et al.  Co-response analysis: a new experimental strategy for metabolic control analysis. , 1996, Journal of theoretical biology.

[37]  M. Reuss,et al.  In vivo analysis of glucose-induced fast changes in yeast adenine nucleotide pool applying a rapid sampling technique. , 1993, Analytical biochemistry.

[38]  J. Varner,et al.  Metabolic Engineering from a Cybernetic Perspective. 2. Qualitative Investigation of Nodal Architechtures and Their Response to Genetic Perturbation , 1999, Biotechnology progress.

[39]  B. Palsson,et al.  Metabolic Capabilities of Escherichia coli II. Optimal Growth Patterns , 1993 .

[40]  J. Liao,et al.  Experimental determination of flux control distribution in biochemical systems: In vitro model to analyze transient metabolite concentrations , 1993, Biotechnology and bioengineering.

[41]  J. Heijnen.  Approximative kinetic formats used in metabolic network modeling , 2005, Biotechnology and bioengineering.

[42]  H. Westerhoff,et al.  Transcriptome meets metabolome: hierarchical and metabolic regulation of the glycolytic pathway , 2001, FEBS letters.

[43]  J J Heijnen,et al.  Metabolic flux control analysis of branch points: an improved approach to obtain flux control coefficients from large perturbation data. , 2004, Metabolic engineering.

[44]  Christoph Giersch.  Determining Elasticities from Multiple Measurements of Steady-state Flux Rates and Metabolite Concentrations: Theory , 1994 .

[45]  Eberhard O. Voit,et al.  Computational Analysis of Biochemical Systems: A Practical Guide for Biochemists and Molecular Biologists , 2000 .

[46]  J. Heijnen,et al.  Metabolome dynamic responses of Saccharomyces cerevisiae to simultaneous rapid perturbations in external electron acceptor and electron donor. , 2007, FEMS yeast research.

[47]  D. Ramkrishna,et al.  Metabolic Engineering from a Cybernetic Perspective. 1. Theoretical Preliminaries , 1999, Biotechnology progress.

[48]  J. Heijnen,et al.  Quantitative analysis of the microbial metabolome by isotope dilution mass spectrometry using uniformly 13C-labeled cell extracts as internal standards. , 2005, Analytical biochemistry.

[49]  Jonas S. Almeida,et al.  Decoupling dynamical systems for pathway identification from metabolic profiles , 2004, Bioinform..

[50]  H. Kacser,et al.  MOlecular democracy: who shares the controls? , 1979, Biochemical Society transactions.

[51]  J. Heijnen,et al.  Dynamic simulation and metabolic re-design of a branched pathway using linlog kinetics. , 2003, Metabolic engineering.

[52]  Dynamic Metabolic Control Theory , 1992 .