OpenFLUX2: 13C-MFA modeling software package adjusted for the comprehensive analysis of single and parallel labeling experiments

BackgroundSteady-state 13C-based metabolic flux analysis (13C-MFA) is the most powerful method available for the quantification of intracellular fluxes. These analyses include concertedly linked experimental and computational stages: (i) assuming the metabolic model and optimizing the experimental design; (ii) feeding the investigated organism using a chosen 13C-labeled substrate (tracer); (iii) measuring the extracellular effluxes and detecting the 13C-patterns of intracellular metabolites; and (iv) computing flux parameters that minimize the differences between observed and simulated measurements, followed by evaluating flux statistics. In its early stages, 13C-MFA was performed on the basis of data obtained in a single labeling experiment (SLE) followed by exploiting the developed high-performance computational software. Recently, the advantages of parallel labeling experiments (PLEs), where several LEs are conducted under the conditions differing only by the tracer(s) choice, were demonstrated, particularly with regard to improving flux precision due to the synergy of complementary information. The availability of an open-source software adjusted for PLE-based 13C-MFA is an important factor for PLE implementation.ResultsThe open-source software OpenFLUX, initially developed for the analysis of SLEs, was extended for the computation of PLE data. Using the OpenFLUX2, in silico simulation confirmed that flux precision is improved when 13C-MFA is implemented by fitting PLE data to the common model compared with SLE-based analysis. Efficient flux resolution could be achieved in the PLE-mediated analysis when the choice of tracer was based on an experimental design computed to minimize the flux variances from different parts of the metabolic network. The analysis provided by OpenFLUX2 mainly includes (i) the optimization of the experimental design, (ii) the computation of the flux parameters from LEs data, (iii) goodness-of-fit testing of the model’s adequacy, (iv) drawing conclusions concerning the identifiability of fluxes and construction of a contribution matrix reflecting the relative contribution of the measurement variances to the flux variances, and (v) precise determination of flux confidence intervals using a fine-tunable and convergence-controlled Monte Carlo-based method.ConclusionsThe developed open-source OpenFLUX2 provides a friendly software environment that facilitates beginners and existing OpenFLUX users to implement LEs for steady-state 13C-MFA including experimental design, quantitative evaluation of flux parameters and statistics.

[1]  U. Sauer,et al.  GC‐MS Analysis of Amino Acids Rapidly Provides Rich Information for Isotopomer Balancing , 2000, Biotechnology progress.

[2]  Tae Hoon Yang,et al.  13C-based metabolic flux analysis: fundamentals and practice. , 2013, Methods in molecular biology.

[3]  M. Dauner.  From fluxes and isotope labeling patterns towards in silico cells. , 2010, Current opinion in biotechnology.

[4]  B. Christensen,et al.  Isotopomer analysis using GC-MS. , 1999, Metabolic engineering.

[5]  G. Stephanopoulos,et al.  Systematic quantification of complex metabolic flux networks using stable isotopes and mass spectrometry. , 2003, European journal of biochemistry.

[6]  Nicola Zamboni,et al.  13C metabolic flux analysis in complex systems. , 2011, Current opinion in biotechnology.

[7]  Taghi M. Khoshgoftaar,et al.  Controlling Overfitting in Classification-Tree Models of Software Quality , 2001, Empirical Software Engineering.

[8]  J J Heijnen,et al.  A priori analysis of metabolic flux identifiability from (13)C-labeling data. , 2001, Biotechnology and bioengineering.

[9]  J E Bailey,et al.  Metabolic flux analysis with a comprehensive isotopomer model in Bacillus subtilis. , 2001, Biotechnology and bioengineering.

[10]  Christopher P. Long,et al.  Quantifying biomass composition by gas chromatography/mass spectrometry. , 2014, Analytical chemistry.

[11]  Andrej Pázman,et al.  Foundations of Optimum Experimental Design , 1986 .

[12]  Jean-Charles Portais,et al.  IsoDesign: a software for optimizing the design of 13C-metabolic flux analysis experiments. , 2014, Biotechnology and bioengineering.

[13]  Anthony S. Bryk,et al.  Hierarchical Linear Models: Applications and Data Analysis Methods , 1992 .

[14]  Wolfgang Wiechert,et al.  Metabolic isotopomer labeling systems. Part II: structural flux identifiability analysis. , 2003, Mathematical biosciences.

[15]  Christoph Wittmann,et al.  Amplified Expression of Fructose 1,6-Bisphosphatase in Corynebacterium glutamicum Increases In Vivo Flux through the Pentose Phosphate Pathway and Lysine Production on Different Carbon Sources , 2005, Applied and Environmental Microbiology.

[16]  L. Eggeling,et al.  Handbook of Corynebacterium glutamicum , 2005 .

[17]  Gregory Stephanopoulos,et al.  Accurate assessment of amino acid mass isotopomer distributions for metabolic flux analysis. , 2007, Analytical chemistry.

[18]  W. Wiechert,et al.  In Vivo Quantification of Parallel and Bidirectional Fluxes in the Anaplerosis of Corynebacterium glutamicum * , 2000, The Journal of Biological Chemistry.

[19]  Stephan Noack,et al.  Improved L‐lysine production with Corynebacterium glutamicum and systemic insight into citrate synthase flux and activity , 2012, Biotechnology and bioengineering.

[20]  A. D. de Graaf,et al.  Flux partitioning in the split pathway of lysine synthesis in Corynebacterium glutamicum. Quantification by 13C- and 1H-NMR spectroscopy. , 1993, European journal of biochemistry.

[21]  P. Verheijen,et al.  Cumulative bondomers: a new concept in flux analysis from 2D [13C,1H] COSY NMR data. , 2002, Biotechnology and bioengineering.

[22]  Albert A. de Graaf,et al.  Metabolic Flux Analysis of Corynebacterium glutamicum , 2000 .

[23]  H. Sahm,et al.  Identification of glyA (Encoding Serine Hydroxymethyltransferase) and Its Use Together with the Exporter ThrE To Increase l-Threonine Accumulation by Corynebacterium glutamicum , 2002, Applied and Environmental Microbiology.

[24]  U. Sauer,et al.  Large-scale 13C-flux analysis reveals distinct transcriptional control of respiratory and fermentative metabolism in Escherichia coli , 2011, Molecular systems biology.

[25]  Wolfgang M. Hartmann,et al.  Computing the Moore-Penrose Inverse for the Covariance Matrix in Constrained Nonlinear Estimation , 1996, SIAM J. Optim..

[26]  Bernhard O. Palsson,et al.  Predicting outcomes of steady-state 13C isotope tracing experiments using Monte Carlo sampling , 2012, BMC Systems Biology.

[27]  Gregory Stephanopoulos,et al.  Quantifying Reductive Carboxylation Flux of Glutamine to Lipid in a Brown Adipocyte Cell Line* , 2008, Journal of Biological Chemistry.

[28]  Scott B. Crown,et al.  Parallel labeling experiments and metabolic flux analysis: Past, present and future methodologies. , 2013, Metabolic engineering.

[29]  Christoph Wittmann,et al.  Comparative Metabolic Flux Analysis of Lysine-Producing Corynebacterium glutamicum Cultured on Glucose or Fructose , 2004, Applied and Environmental Microbiology.

[30]  S. Billings,et al.  Metabolic flux distribution analysis by 13C-tracer experiments using the Markov chain-Monte Carlo method. , 2005, Biochemical Society transactions.

[31]  Maciek R Antoniewicz,et al.  Resolving the TCA cycle and pentose-phosphate pathway of Clostridium acetobutylicum ATCC 824: Isotopomer analysis, in vitro activities and expression analysis. , 2011, Biotechnology journal.

[32]  Maciek R Antoniewicz,et al.  Selection of tracers for 13C-metabolic flux analysis using elementary metabolite units (EMU) basis vector methodology. , 2012, Metabolic engineering.

[33]  L. Foster.  Rank and null space calculations using matrix decomposition without column interchanges , 1986 .

[34]  M. Antoniewicz.  Using multiple tracers for 13C metabolic flux analysis. , 2013, Methods in molecular biology.

[35]  Christoph Wittmann,et al.  Metabolic responses to pyruvate kinase deletion in lysine producing Corynebacterium glutamicum , 2008, Microbial cell factories.

[36]  Yonina C. Eldar,et al.  Fluxomers: a new approach for 13C metabolic flux analysis , 2011, BMC Systems Biology.

[37]  W. Wiechert,et al.  Bidirectional reaction steps in metabolic networks: I. Modeling and simulation of carbon isotope labeling experiments. , 1997, Biotechnology and bioengineering.

[38]  B. Pfeifer,et al.  Metabolic flux analysis and pharmaceutical production. , 2010, Metabolic engineering.

[39]  G. Stephanopoulos,et al.  Application of macroscopic balances to the identification of gross measurement errors , 1983, Biotechnology and bioengineering.

[40]  Gregory Stephanopoulos,et al.  Determination of confidence intervals of metabolic fluxes estimated from stable isotope measurements. , 2006, Metabolic engineering.

[41]  Oliver Frick,et al.  Hybrid optimization for 13C metabolic flux analysis using systems parametrized by compactification , 2008, BMC Systems Biology.

[42]  A. D. de Graaf,et al.  Analysis of carbon metabolism in Escherichia coli strains with an inactive phosphotransferase system by (13)C labeling and NMR spectroscopy. , 2002, Metabolic Engineering.

[43]  T Szyperski,et al.  13C-NMR, MS and metabolic flux balancing in biotechnology research , 1998, Quarterly Reviews of Biophysics.

[44]  M. Antoniewicz,et al.  COMPLETE-MFA: complementary parallel labeling experiments technique for metabolic flux analysis. , 2013, Metabolic engineering.

[45]  Madhukar S. Dasika,et al.  Metabolic flux elucidation for large-scale models using 13C labeled isotopes. , 2007, Metabolic engineering.

[46]  B. Palsson,et al.  Metabolic Flux Balancing: Basic Concepts, Scientific and Practical Use , 1994, Bio/Technology.

[47]  Christoph Wittmann,et al.  Correcting mass isotopomer distributions for naturally occurring isotopes. , 2002, Biotechnology and bioengineering.

[48]  L. Quek,et al.  OpenFLUX: efficient modelling software for 13C-based metabolic flux analysis , 2009, Microbial cell factories.

[49]  G. Sriram,et al.  Designer labels for plant metabolism: statistical design of isotope labeling experiments for improved quantification of flux in complex plant metabolic networks. , 2013, Molecular bioSystems.

[50]  H. Shimizu,et al.  OpenMebius: An Open Source Software for Isotopically Nonstationary 13C-Based Metabolic Flux Analysis , 2014, BioMed research international.

[51]  W Wiechert,et al.  Metabolic isotopomer labeling systems. Part I: global dynamic behavior. , 2001, Mathematical biosciences.

[52]  Scott B. Crown,et al.  Rational design of 13C-labeling experiments for metabolic flux analysis in mammalian cells , 2012, BMC Systems Biology.

[53]  Elmar Heinzle,et al.  Numerical bias estimation for mass spectrometric mass isotopomer analysis. , 2009, Analytical biochemistry.

[54]  U. Sauer,et al.  Article number: 62 REVIEW Metabolic networks in motion: 13 C-based flux analysis , 2022 .

[55]  Elmar Heinzle,et al.  Metabolic flux analysis in eukaryotes. , 2010, Current opinion in biotechnology.

[56]  Christoph Wittmann,et al.  Genealogy Profiling through Strain Improvement by Using Metabolic Network Analysis: Metabolic Flux Genealogy of Several Generations of Lysine-Producing Corynebacteria , 2002, Applied and Environmental Microbiology.

[57]  G. Stephanopoulos,et al.  Elementary metabolite units (EMU): a novel framework for modeling isotopic distributions. , 2007, Metabolic engineering.

[58]  George Stephanopoulos,et al.  Modeling of Isotope Distributions and Intracellular Fluxes in Metabolic Networks Using Atom Mapping Matrices , 1994 .

[59]  Kazuyuki Shimizu,et al.  An improved method for statistical analysis of metabolic flux analysis using isotopomer mapping matrices with analytical expressions. , 2003, Journal of biotechnology.

[60]  W Wiechert,et al.  A universal framework for 13C metabolic flux analysis. , 2001, Metabolic engineering.

[61]  Jorge Nocedal,et al.  An Interior Point Algorithm for Large-Scale Nonlinear Programming , 1999, SIAM J. Optim..

[62]  Maciek R Antoniewicz,et al.  Tandem mass spectrometry for measuring stable-isotope labeling. , 2013, Current opinion in biotechnology.

[63]  C. Wittmann,et al.  Modeling and experimental design for metabolic flux analysis of lysine-producing Corynebacteria by mass spectrometry. , 2001, Metabolic engineering.

[64]  M. Antoniewicz,et al.  Parallel labeling experiments with [U-13C]glucose validate E. coli metabolic network model for 13C metabolic flux analysis. , 2012, Metabolic engineering.

[65]  J. Heijnen,et al.  Cumulative Bondomers: A New Concept in Flux Analysis from 2D ( 13 C, 1 H) COSY , 2002 .

[66]  C. Maranas,et al.  Identification of optimal measurement sets for complete flux elucidation in metabolic flux analysis experiments. , 2008, Biotechnology and bioengineering.

[67]  Patrick F Suthers,et al.  Construction of an E. Coli genome‐scale atom mapping model for MFA calculations , 2011, Biotechnology and bioengineering.

[68]  G. Stephanopoulos.  CHAPTER 5 – Regulation of Metabolic Pathways , 1998 .

[69]  J. Villadsen,et al.  Modeling isotopomer distributions in biochemical networks using isotopomer mapping matrices. , 1997, Biotechnology and bioengineering.

[70]  Gregory Stephanopoulos,et al.  Kinetic isotope effects significantly influence intracellular metabolite [superscript 13]C labeling patterns and flux determination , 2013 .

[71]  Yinjie J. Tang,et al.  Central metabolic responses to the overproduction of fatty acids in Escherichia coli based on 13C‐metabolic flux analysis , 2014, Biotechnology and bioengineering.

[72]  Jean-Charles Portais,et al.  Influx_s: Increasing Numerical Stability and Precision for Metabolic Flux Analysis in Isotope Labelling Experiments , 2012, Bioinform..

[73]  U. Sauer,et al.  Metabolic flux profiling of Escherichia coli mutants in central carbon metabolism using GC-MS. , 2003, European journal of biochemistry.

[74]  Yohei Yamada,et al.  Metabolic flux analysis in biotechnology processes , 2008, Biotechnology Letters.

[75]  Christoph Wittmann,et al.  Metabolic network analysis of lysine producing Corynebacterium glutamicum at a miniaturized scale , 2004, Biotechnology and bioengineering.

[76]  G. Stephanopoulos,et al.  Metabolic Engineering: Principles And Methodologies , 1998 .

[77]  Christoph Wittmann,et al.  Metabolic fluxes and beyond—systems biology understanding and engineering of microbial metabolism , 2010, Applied Microbiology and Biotechnology.

[78]  W. Wiechert,et al.  Bidirectional reaction steps in metabolic networks: II. Flux estimation and statistical analysis. , 1997, Biotechnology and bioengineering.

[79]  Maciek R Antoniewicz,et al.  Publishing 13C metabolic flux analysis studies: a review and future perspectives. , 2013, Metabolic engineering.

[80]  Wolfgang Wiechert,et al.  Collisional fragmentation of central carbon metabolites in LC‐MS/MS increases precision of 13C metabolic flux analysis , 2012, Biotechnology and bioengineering.

[81]  Wolfgang Wiechert,et al.  13CFLUX2—high-performance software suite for 13C-metabolic flux analysis , 2012, Bioinform..

[82]  S. Noack,et al.  Comparative 13C Metabolic Flux Analysis of Pyruvate Dehydrogenase Complex-Deficient, l-Valine-Producing Corynebacterium glutamicum , 2011, Applied and Environmental Microbiology.

[83]  W. Wiechert.  13C metabolic flux analysis. , 2001, Metabolic engineering.

[84]  W. Wiechert,et al.  Bidirectional reaction steps in metabolic networks: III. Explicit solution and analysis of isotopomer labeling systems. , 1999, Biotechnology and bioengineering.

[85]  Juho Rousu,et al.  Planning optimal measurements of isotopomer distributions for estimation of metabolic fluxes , 2006, Bioinform..

[86]  U. Sauer,et al.  High-throughput metabolic flux analysis based on gas chromatography-mass spectrometry derived 13C constraints. , 2004, Analytical biochemistry.

[87]  C. Wittmann,et al.  From zero to hero--design-based systems metabolic engineering of Corynebacterium glutamicum for L-lysine production. , 2011, Metabolic engineering.

[88]  W Wiechert,et al.  Bidirectional reaction steps in metabolic networks: IV. Optimal design of isotopomer labeling experiments. , 1999, Biotechnology and bioengineering.

[89]  Christoph Wittmann,et al.  Fluxome analysis using GC-MS , 2007, Microbial cell factories.

[90]  Christoph Wittmann,et al.  Increased lysine production by flux coupling of the tricarboxylic acid cycle and the lysine biosynthetic pathway--metabolic engineering of the availability of succinyl-CoA in Corynebacterium glutamicum. , 2013, Metabolic engineering.

[91]  U. Sauer,et al.  13C-based metabolic flux analysis , 2009, Nature Protocols.