Metabolic Flux Estimation-A Self-Adaptive Evolutionary Algorithm with Singular Value Decomposition

Metabolic flux analysis is important for metabolic system regulation and intracellular pathway identification. A popular approach for intracellular flux estimation involves using 13C tracer experiments to label states that can be measured by nuclear magnetic resonance spectrometry or gas chromatography mass spectrometry. However, the bilinear balance equations derived from 13C tracer experiments and the noisy measurements require a nonlinear optimization approach to obtain the optimal solution. In this paper, the flux quantification problem is formulated as an error-minimization problem with equality and inequality constraints through the 13C balance and stoichiometric equations. The stoichiometric constraints are transformed to a null space by singular value decomposition. Self-adaptive evolutionary algorithms are then introduced for flux quantification. The performance of the evolutionary algorithm is compared with ordinary least squares estimation by the simulation of the central pentose phosphate pathway. The proposed algorithm is also applied to the central metabolism of Corynebacterium glutamicum under lysine-producing conditions. A comparison between the results from the proposed algorithm and data from the literature is given. The complexity of a metabolic system with bidirectional reactions is also investigated by analyzing the fluctuations in the flux estimates when available measurements are varied.

[1]  E. R. Andrew,et al.  Nuclear Magnetic Resonance , 1955 .

[2]  John H. Holland,et al.  Outline for a Logical Theory of Adaptive Systems , 1962, JACM.

[3]  David H. Anderson.  Compartmental Modeling and Tracer Kinetics , 1983 .

[4]  G. Stephanopoulos,et al.  Computer‐aided synthesis of biochemical pathways , 1990, Biotechnology and bioengineering.

[5]  David B. Fogel,et al.  Meta-evolutionary programming , 1991, [1991] Conference Record of the Twenty-Fifth Asilomar Conference on Signals, Systems & Computers.

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

[7]  Heinz Mühlenbein,et al.  Predictive Models for the Breeder Genetic Algorithm I. Continuous Parameter Optimization , 1993, Evolutionary Computation.

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

[9]  G. Stephanopoulos,et al.  Metabolic flux distributions in Corynebacterium glutamicum during growth and lysine overproduction , 2000, Biotechnology and bioengineering.

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

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

[12]  S. Isaacs,et al.  AN EVOLUTIONARY ALGORITHM FOR INITIAL STATE AND PARAMETER ESTIMATION IN COMPLEX BIOCHEMICAL MODELS , 1995 .

[13]  J Villadsen,et al.  Metabolic flux distributions in Penicillium chrysogenum during fed‐batch cultivations , 1995, Biotechnology and bioengineering.

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

[15]  D. Fogel,et al.  A comparison of methods for self-adaptation in evolutionary algorithms. , 1995, Bio Systems.

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

[17]  H Sahm,et al.  Determination of the fluxes in the central metabolism of Corynebacterium glutamicum by nuclear magnetic resonance spectroscopy combined with metabolite balancing , 1996, Biotechnology and bioengineering.

[18]  H. Bonarius,et al.  Flux analysis of underdetermined metabolic networks: the quest for the missing constraints. , 1997 .

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

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

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

[22]  B. Palsson,et al.  The underlying pathway structure of biochemical reaction networks. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[23]  David B. Fogel,et al.  Inductive reasoning and bounded rationality reconsidered , 1999, IEEE Trans. Evol. Comput..

[24]  J Villadsen,et al.  Quantification of intracellular metabolic fluxes from fractional enrichment and 13C-13C coupling constraints on the isotopomer distribution in labeled biomass components. , 1999, Metabolic engineering.

[25]  R. Heinrich,et al.  Metabolic Pathway Analysis: Basic Concepts and Scientific Applications in the Post‐genomic Era , 1999, Biotechnology progress.

[26]  D. Fell,et al.  Detection of elementary flux modes in biochemical networks: a promising tool for pathway analysis and metabolic engineering. , 1999, Trends in biotechnology.

[27]  J. Nielsen,et al.  Quantitative analysis of metabolic fluxes in Escherichia coli, using two-dimensional NMR spectroscopy and complete isotopomer models. , 1999, Journal of biotechnology.

[28]  D. Fell,et al.  A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks , 2000, Nature Biotechnology.

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

[30]  Elena Marchiori,et al.  How to handle constraints with evolutionary algorithms , 2001 .

[31]  A. Eiben.  Evolutionary algorithms and constraint satisfaction: definitions, survey, methodology, and research directions , 2001 .

[32]  Reinhart Heinrich,et al.  Evolutionary optimization of metabolic pathways. Theoretical reconstruction of the stoichiometry of ATP and NADH producing systems , 2001, Bulletin of mathematical biology.

[33]  Jason A. Papin,et al.  Determination of redundancy and systems properties of the metabolic network of Helicobacter pylori using genome-scale extreme pathway analysis. , 2002, Genome research.

[34]  Ertugrul M. Ozbudak,et al.  Regulation of noise in the expression of a single gene , 2002, Nature Genetics.

[35]  P. Swain,et al.  Stochastic Gene Expression in a Single Cell , 2002, Science.

[36]  T. Strovas,et al.  Quantification of central metabolic fluxes in the facultative methylotroph methylobacterium extorquens AM1 using 13C‐label tracing and mass spectrometry , 2003, Biotechnology and bioengineering.

[37]  Sorin Draghici,et al.  Predicting HIV drug resistance with neural networks , 2003, Bioinform..

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

[39]  A. E. Eiben,et al.  Introduction to Evolutionary Computing , 2003, Natural Computing Series.

[40]  David Corne,et al.  Evolutionary Computation In Bioinformatics , 2003 .

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

[42]  Jason A. Papin,et al.  Analysis of metabolic capabilities using singular value decomposition of extreme pathway matrices. , 2003, Biophysical journal.

[43]  Gregory Stephanopoulos,et al.  Ion-trap mass spectrometry used in combination with gas chromatography for high-resolution metabolic flux determination. , 2003, BioTechniques.

[44]  Mads Kærn,et al.  Noise in eukaryotic gene expression , 2003, Nature.

[45]  Visakan Kadirkamanathan,et al.  Self-adaptive evolutionary algorithm based methods for quantification in metabolic systems , 2004, 2004 Symposium on Computational Intelligence in Bioinformatics and Computational Biology.

[46]  R. Mahadevan,et al.  Using metabolic flux data to further constrain the metabolic solution space and predict internal flux patterns: the Escherichia coli spectrum , 2004, Biotechnology and bioengineering.

[47]  Hans-Paul Schwefel,et al.  Evolution strategies – A comprehensive introduction , 2002, Natural Computing.

[48]  Phillip C Wright,et al.  Novel approach for peptide quantitation and sequencing based on 15N and 13C metabolic labeling. , 2005, Journal of proteome research.

[49]  G. Stephanopoulos,et al.  Strain improvement by metabolic engineering: lysine production as a case study for systems biology. , 2005, Current opinion in biotechnology.

[50]  Takashi Gojobori,et al.  Comparative study of flux redistribution of metabolic pathway in glutamate production by two coryneform bacteria. , 2005, Metabolic engineering.

[51]  Stefan Schuster,et al.  YANA – a software tool for analyzing flux modes, gene-expression and enzyme activities , 2005, BMC Bioinformatics.

[52]  Robert Urbanczik,et al.  An improved algorithm for stoichiometric network analysis: theory and applications , 2005, Bioinform..

[53]  Yun S. Song,et al.  Efficient computation of close lower and upper bounds on the minimum number of recombinations in biological sequence evolution , 2005, ISMB.