Dynamic elementary mode modelling of non-steady state flux data

BackgroundA novel framework is proposed to analyse metabolic fluxes in non-steady state conditions, based on the new concept of dynamic elementary mode (dynEM): an elementary mode activated partially depending on the time point of the experiment.ResultsTwo methods are introduced here: dynamic elementary mode analysis (dynEMA) and dynamic elementary mode regression discriminant analysis (dynEMR-DA). The former is an extension of the recently proposed principal elementary mode analysis (PEMA) method from steady state to non-steady state scenarios. The latter is a discriminant model that permits to identify which dynEMs behave strongly different depending on the experimental conditions. Two case studies of Saccharomyces cerevisiae, with fluxes derived from simulated and real concentration data sets, are presented to highlight the benefits of this dynamic modelling.ConclusionsThis methodology permits to analyse metabolic fluxes at early stages with the aim of i) creating reduced dynamic models of flux data, ii) combining many experiments in a single biologically meaningful model, and iii) identifying the metabolic pathways that drive the organism from one state to another when changing the environmental conditions.

[1]  C. Leão,et al.  Sugar Metabolism in Yeasts: an Overview of Aerobic and Anaerobic Glucose Catabolism , 2006 .

[2]  K. Larsson,et al.  A gene encoding sn‐glycerol 3‐phosphate dehydrogenase (NAD+) complements an osmosensitive mutant of Saccharomyces cerevisiae , 1993, Molecular microbiology.

[3]  Meike T. Wortel,et al.  Lost in Transition: Start-Up of Glycolysis Yields Subpopulations of Nongrowing Cells , 2014, Science.

[4]  Bernhard O. Palsson,et al.  Decomposing complex reaction networks using random sampling, principal component analysis and basis rotation , 2009, BMC Systems Biology.

[5]  Age K Smilde,et al.  Global test for metabolic pathway differences between conditions. , 2012, Analytica chimica acta.

[6]  Age K. Smilde,et al.  Principal Component Analysis , 2003, Encyclopedia of Machine Learning.

[7]  中尾 光輝,et al.  KEGG(Kyoto Encyclopedia of Genes and Genomes)〔和文〕 (特集 ゲノム医学の現在と未来--基礎と臨床) -- (データベース) , 2000 .

[8]  Kiyoko F. Aoki-Kinoshita,et al.  From genomics to chemical genomics: new developments in KEGG , 2005, Nucleic Acids Res..

[9]  Holger Fröhlich,et al.  Modeling ERBB receptor-regulated G1/S transition to find novel targets for de novo trastuzumab resistance , 2009, BMC Systems Biology.

[10]  A. Blomberg,et al.  Purification and characterization of two isoenzymes of DL-glycerol-3-phosphatase from Saccharomyces cerevisiae. Identification of the corresponding GPP1 and GPP2 genes and evidence for osmotic regulation of Gpp2p expression by the osmosensing mitogen-activated protein kinase signal transduction path , 1996, The Journal of biological chemistry.

[11]  Wei-Hao Wang,et al.  Studies , 1926 .

[12]  Jesús Picó,et al.  MCR-ALS on metabolic networks: Obtaining more meaningful pathways , 2015 .

[13]  F. Llaneras,et al.  Stoichiometric modelling of cell metabolism. , 2008, Journal of bioscience and bioengineering.

[14]  R. Bro.  Multiway calibration. Multilinear PLS , 1996 .

[15]  J. Heijnen,et al.  Determination of the cytosolic free NAD/NADH ratio in Saccharomyces cerevisiae under steady‐state and highly dynamic conditions , 2008, Biotechnology and bioengineering.

[16]  M. Barker,et al.  Partial least squares for discrimination , 2003 .

[17]  Romà Tauler,et al.  A graphical user-friendly interface for MCR-ALS: a new tool for multivariate curve resolution in MATLAB , 2005 .

[18]  L. Petzold.  Automatic Selection of Methods for Solving Stiff and Nonstiff Systems of Ordinary Differential Equations , 1983 .

[19]  Age K. Smilde,et al.  Double-check: validation of diagnostic statistics for PLS-DA models in metabolomics studies , 2011, Metabolomics.

[20]  A. Blomberg,et al.  Cloning and characterization of GPD2, a second gene encoding sn‐glycerol 3‐phosphate dehydrogenase (NAD+) in Saccharomyces cerevisiae, and its comparison with GPD1 , 1995, Molecular microbiology.

[21]  Fabian J Theis,et al.  Statistical methods for the analysis of high-throughput metabolomics data , 2013, Computational and structural biotechnology journal.

[22]  Diana M. Hendrickx,et al.  MetDFBA: incorporating time-resolved metabolomics measurements into dynamic flux balance analysis. , 2015, Molecular bioSystems.

[23]  Rui Oliveira,et al.  Principal elementary mode analysis (PEMA). , 2016, Molecular bioSystems.

[24]  Peter J.T. Verheijen,et al.  Construction of kinetic models for metabolic reaction networks: Lessons learned in analysing short-term stimulus response data , 2011 .

[25]  Age K. Smilde,et al.  UvA-DARE ( Digital Academic Repository ) Assessment of PLSDA cross validation , 2008 .

[26]  Mudita Singhal,et al.  COPASI - a COmplex PAthway SImulator , 2006, Bioinform..

[27]  L. Hood.  Systems biology: integrating technology, biology, and computation , 2003, Mechanisms of Ageing and Development.

[28]  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.

[29]  Susumu Goto,et al.  KEGG for representation and analysis of molecular networks involving diseases and drugs , 2009, Nucleic Acids Res..

[30]  Abel Folch-Fortuny,et al.  Metabolic Flux Understanding of Pichia pastoris Grown on Heterogenous , 2014 .

[31]  Rasmus Bro,et al.  The N-way Toolbox for MATLAB , 2000 .

[32]  Jörg Stelling,et al.  Large-scale computation of elementary flux modes with bit pattern trees , 2008, Bioinform..

[33]  F. Doyle,et al.  Dynamic flux balance analysis of diauxic growth in Escherichia coli. , 2002, Biophysical journal.

[34]  Anders Blomberg,et al.  Purification and Characterization of Two Isoenzymes of DL-Glycerol-3-phosphatase from Saccharomyces cerevisiae , 1996, The Journal of Biological Chemistry.