A Novel Method for Flux Distribution Computation in Metabolic Networks

In recent years, the study on metabolic networks has attracted considerable attention from the research community. Though the topological structures of genome-scale metabolic networks of some organisms have been investigated, their metabolic flux distributions still remain unclear. The understanding of flux distributions in metabolic networks, especially when it comes to the gene-knockout mutants, is helpful for suggesting potential ways to improve strain design. The traditional method of flux distribution computation, i.e., flux balance analysis (FBA) method, is based on the idea of maximizing biomass yield. However, this method overestimates the production of biomass. In this paper, we develop a novel approach to overcome the drawback of the FBA method. First, we adopt a series of extended equations to model reaction flux; Second, we build the stoichiometric matrix of a metabolic network by using a more complex but accurate model - carbon mole balance - rather than mass balance used in FBA. Computation results with real-world data of Escherichia coli show that our approach outperforms FBA in the accuracy of flux distribution computation.

[1]  B. Palsson,et al.  Saccharomyces cerevisiae phenotypes can be predicted by using constraint-based analysis of a genome-scale reconstructed metabolic network , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[2]  A. Burgard,et al.  Optimization-based framework for inferring and testing hypothesized metabolic objective functions. , 2003, Biotechnology and bioengineering.

[3]  Kenneth J. Kauffman,et al.  Advances in flux balance analysis. , 2003, Current opinion in biotechnology.

[4]  R. Mahadevan,et al.  The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. , 2003, Metabolic engineering.

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

[6]  E. Ruppin,et al.  Regulatory on/off minimization of metabolic flux changes after genetic perturbations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[7]  U. Sauer,et al.  Large-scale in vivo flux analysis shows rigidity and suboptimal performance of Bacillus subtilis metabolism , 2005, Nature Genetics.

[8]  Erwin P. Gianchandani,et al.  Flux balance analysis in the era of metabolomics , 2006, Briefings Bioinform..

[9]  J. Pronk,et al.  Role of Transcriptional Regulation in Controlling Fluxes in Central Carbon Metabolism of Saccharomyces cerevisiae , 2004, Journal of Biological Chemistry.

[10]  B. Palsson,et al.  An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR) , 2003, Genome Biology.

[11]  Costas D Maranas,et al.  OptStrain: a computational framework for redesign of microbial production systems. , 2004, Genome research.

[12]  Kazuyuki Shimizu,et al.  Metabolic flux analysis based on 13C-labeling experiments and integration of the information with gene and protein expression patterns. , 2004, Advances in biochemical engineering/biotechnology.

[13]  C. Schilling,et al.  Flux coupling analysis of genome-scale metabolic network reconstructions. , 2004, Genome research.

[14]  B. Palsson.  The challenges of in silico biology , 2000, Nature Biotechnology.

[15]  C. Maranas,et al.  An optimization framework for identifying reaction activation/inhibition or elimination candidates for overproduction in microbial systems. , 2006, Metabolic engineering.

[16]  Stephen S Fong,et al.  Metabolic gene–deletion strains of Escherichia coli evolve to computationally predicted growth phenotypes , 2004, Nature Genetics.

[17]  S. Oliver,et al.  Chance and necessity in the evolution of minimal metabolic networks , 2006, Nature.

[18]  S. Schuster,et al.  Metabolic network structure determines key aspects of functionality and regulation , 2002, Nature.

[19]  A. Burgard,et al.  Optknock: A bilevel programming framework for identifying gene knockout strategies for microbial strain optimization , 2003, Biotechnology and bioengineering.

[20]  B. Palsson,et al.  Genome-scale models of microbial cells: evaluating the consequences of constraints , 2004, Nature Reviews Microbiology.

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

[22]  B. Palsson,et al.  Reconstructing metabolic flux vectors from extreme pathways: defining the α-spectrum , 2003 .

[23]  B. Palsson,et al.  Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110 , 1994, Applied and environmental microbiology.

[24]  H. Mori,et al.  Global metabolic response of Escherichia coli to gnd or zwf gene-knockout, based on 13C-labeling experiments and the measurement of enzyme activities , 2004, Applied Microbiology and Biotechnology.

[25]  I. Grossmann,et al.  Recursive MILP model for finding all the alternate optima in LP models for metabolic networks , 2000 .

[26]  Bernhard O. Palsson,et al.  Identification of Genome-Scale Metabolic Network Models Using Experimentally Measured Flux Profiles , 2006, PLoS Comput. Biol..

[27]  M. Araúzo-Bravo,et al.  Effect of a pyruvate kinase (pykF-gene) knockout mutation on the control of gene expression and metabolic fluxes in Escherichia coli. , 2004, FEMS microbiology letters.