OpenFLUX: efficient modelling software for 13C-based metabolic flux analysis

BackgroundThe quantitative analysis of metabolic fluxes, i.e., in vivo activities of intracellular enzymes and pathways, provides key information on biological systems in systems biology and metabolic engineering. It is based on a comprehensive approach combining (i) tracer cultivation on 13C substrates, (ii) 13C labelling analysis by mass spectrometry and (iii) mathematical modelling for experimental design, data processing, flux calculation and statistics. Whereas the cultivation and the analytical part is fairly advanced, a lack of appropriate modelling software solutions for all modelling aspects in flux studies is limiting the application of metabolic flux analysis.ResultsWe have developed OpenFLUX as a user friendly, yet flexible software application for small and large scale 13C metabolic flux analysis. The application is based on the new Elementary Metabolite Unit (EMU) framework, significantly enhancing computation speed for flux calculation. From simple notation of metabolic reaction networks defined in a spreadsheet, the OpenFLUX parser automatically generates MATLAB-readable metabolite and isotopomer balances, thus strongly facilitating model creation. The model can be used to perform experimental design, parameter estimation and sensitivity analysis either using the built-in gradient-based search or Monte Carlo algorithms or in user-defined algorithms. Exemplified for a microbial flux study with 71 reactions, 8 free flux parameters and mass isotopomer distribution of 10 metabolites, OpenFLUX allowed to automatically compile the EMU-based model from an Excel file containing metabolic reactions and carbon transfer mechanisms, showing it's user-friendliness. It reliably reproduced the published data and optimum flux distributions for the network under study were found quickly (<20 sec).ConclusionWe have developed a fast, accurate application to perform steady-state 13C metabolic flux analysis. OpenFLUX will strongly facilitate and enhance the design, calculation and interpretation of metabolic flux studies. By providing the software open source, we hope it will evolve with the rapidly growing field of fluxomics.

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

[2]  T. Szyperski.  Biosynthetically Directed Fractional 13C‐labeling of Proteinogenic Amino Acids , 1995 .

[3]  T. Szyperski.  Biosynthetically directed fractional 13C-labeling of proteinogenic amino acids. An efficient analytical tool to investigate intermediary metabolism. , 1995, European journal of biochemistry.

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

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

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

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

[8]  U. Sauer,et al.  Bioreaction network topology and metabolic flux ratio analysis by biosynthetic fractional 13C labeling and two-dimensional NMR spectroscopy. , 1999, Metabolic engineering.

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

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

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

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

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

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

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

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

[17]  C. Wittmann,et al.  Metabolic flux analysis using mass spectrometry. , 2002, Advances in biochemical engineering/biotechnology.

[18]  Preben Krabben,et al.  Metabolic flux and metabolic network analysis of Penicillium chrysogenum using 2D [13C, 1H] COSY NMR measurements and cumulative bondomer simulation. , 2003, Biotechnology and bioengineering.

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

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

[21]  C. Wittmann,et al.  In-Depth Profiling of Lysine-Producing Corynebacterium glutamicum by Combined Analysis of the Transcriptome, Metabolome, and Fluxome , 2004, Journal of bacteriology.

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

[23]  C. Des Rosiers,et al.  Profiling substrate fluxes in the isolated working mouse heart using 13C-labeled substrates: focusing on the origin and fate of pyruvate and citrate carbons. , 2004, American journal of physiology. Heart and circulatory physiology.

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

[25]  J. Schwender,et al.  Understanding flux in plant metabolic networks. , 2004, Current opinion in plant biology.

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

[27]  Sang Yup Lee,et al.  Systems biotechnology for strain improvement. , 2005, Trends in biotechnology.

[28]  Nicola Zamboni,et al.  FiatFlux – a software for metabolic flux analysis from 13C-glucose experiments , 2005, BMC Bioinformatics.

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

[30]  B. Palsson,et al.  Isotopomer analysis of cellular metabolism in tissue culture: A comparative study between the pathway and network-based methods , 2007, Metabolomics.

[31]  Visakan Kadirkamanathan,et al.  Markov Chain Monte Carlo Algorithm based metabolic flux distribution analysis on Corynebacterium glutamicum , 2006, Bioinform..

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

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

[34]  G. Stephanopoulos,et al.  Metabolic flux analysis in a nonstationary system: fed-batch fermentation of a high yielding strain of E. coli producing 1,3-propanediol. , 2007, Metabolic engineering.

[35]  Ralf Takors,et al.  Metabolic flux analysis at ultra short time scale: isotopically non-stationary 13C labeling experiments. , 2007, Journal of biotechnology.

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

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

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

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

[40]  Jamey D. Young,et al.  An elementary metabolite unit (EMU) based method of isotopically nonstationary flux analysis , 2008, Biotechnology and bioengineering.

[41]  Juho Rousu,et al.  An analytic and systematic framework for estimating metabolic flux ratios from 13C tracer experiments , 2008, BMC Bioinformatics.

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