Uniform sampling of steady-state flux spaces: means to design experiments and to interpret enzymopathies.

Reconstruction of genome-scale metabolic networks is now possible using multiple different data types. Constraint-based modeling is an approach to interrogate capabilities of reconstructed networks by constraining possible cellular behavior through the imposition of physicochemical laws. As a result, a steady-state flux space is defined that contains all possible functional states of the network. Uniform random sampling of the steady-state flux space allows for the unbiased appraisal of its contents. Monte Carlo sampling of the steady-state flux space of the reconstructed human red blood cell metabolic network under simulated physiologic conditions yielded the following key results: 1), probability distributions for the values of individual metabolic fluxes showed a wide variety of shapes that could not have been inferred without computation; 2), pairwise correlation coefficients were calculated between all fluxes, determining the level of independence between the measurement of any two fluxes, and identifying highly correlated reaction sets; and 3), the network-wide effects of the change in one (or a few) variables (i.e., a simulated enzymopathy or fixing a flux range based on measurements) were computed. Mathematical models provide the most compact and informative representation of a hypothesis of how a cell works. Thus, understanding model predictions clearly is vital to driving forward the iterative model-building procedure that is at the heart of systems biology. Taken together, the Monte Carlo sampling procedure provides a broadening of the constraint-based approach by allowing for the unbiased and detailed assessment of the impact of the applied physicochemical constraints on a reconstructed network.

[1]  S. Lowen.  The Biophysical Journal , 1960, Nature.

[2]  Abbreviations , 1971 .

[3]  A. Grimes.  Human Red Cell Metabolism , 1980 .

[4]  H. Holzhütter,et al.  Mathematical modelling of metabolic pathways affected by an enzyme deficiency. A mathematical model of glycolysis in normal and pyruvate-kinase-deficient red blood cells. , 1985, European journal of biochemistry.

[5]  P. Kuchel,et al.  Regulation of the human-erythrocyte hexose-monophosphate shunt under conditions of oxidative stress: a study using NMR spectroscopy, a kinetic isotope effect, a reconstituted system and computer simulation , 1985 .

[6]  H. Holzhütter,et al.  Interrelations between glycolysis and the hexose monophosphate shunt in erythrocytes as studied on the basis of a mathematical model. , 1988, Bio Systems.

[7]  Tanaka Kr,et al.  Red cell enzymopathies of the glycolytic pathway. , 1990 .

[8]  K. Tanaka,et al.  Red cell enzymopathies of the glycolytic pathway. , 1990, Seminars in hematology.

[9]  H. Holzhütter,et al.  Mathematical modelling of metabolic pathways affected by an enzyme deficiency. , 1990, Biomedica biochimica acta.

[10]  B O Palsson,et al.  Optimal selection of metabolic fluxes for in vivo measurement. I. Development of mathematical methods. , 1992, Journal of theoretical biology.

[11]  B O Palsson,et al.  Optimal selection of metabolic fluxes for in vivo measurement. II. Application to Escherichia coli and hybridoma cell metabolism. , 1992, Journal of theoretical biology.

[12]  Robert L. Smith,et al.  Improving Hit-and-Run for global optimization , 1993, J. Glob. Optim..

[13]  B. Schmeiser,et al.  Performance of the Gibbs, Hit-and-Run, and Metropolis Samplers , 1993 .

[14]  S. Schuster,et al.  ON ELEMENTARY FLUX MODES IN BIOCHEMICAL REACTION SYSTEMS AT STEADY STATE , 1994 .

[15]  H. Holzhütter,et al.  Use of mathematical models for predicting the metabolic effect of large-scale enzyme activity alterations. Application to enzyme deficiencies of red blood cells. , 1995, European journal of biochemistry.

[16]  G. Jacobasch,et al.  Chapter 3 Hemolytic anemias due to erythrocyte enzyme deficiencies , 1996 .

[17]  G. Jacobasch,et al.  Hemolytic anemias due to erythrocyte enzyme deficiencies. , 1996, Molecular aspects of medicine.

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

[19]  Robert L. Smith,et al.  Direction Choice for Accelerated Convergence in Hit-and-Run Sampling , 1998, Oper. Res..

[20]  László Lovász,et al.  Hit-and-run mixes fast , 1999, Math. Program..

[21]  P W Kuchel,et al.  Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: in vivo kinetic characterization of 2,3-bisphosphoglycerate synthase/phosphatase using 13C and 31P NMR. , 1999, The Biochemical journal.

[22]  Juan Carlos Nuño,et al.  METATOOL: for studying metabolic networks , 1999, Bioinform..

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

[24]  J. Edwards,et al.  Systems Properties of the Haemophilus influenzaeRd Metabolic Genotype* , 1999, The Journal of Biological Chemistry.

[25]  P W Kuchel,et al.  Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: computer simulation and metabolic control analysis. , 1999, The Biochemical journal.

[26]  Rui Alves,et al.  Systemic properties of ensembles of metabolic networks: application of graphical and statistical methods to simple unbranched pathways , 2000, Bioinform..

[27]  Rui Alves,et al.  Extending the method of mathematically controlled comparison to include numerical comparisons , 2000, Bioinform..

[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]  B. Palsson,et al.  Assessment of the metabolic capabilities of Haemophilus influenzae Rd through a genome-scale pathway analysis. , 2000, Journal of theoretical biology.

[30]  Rui Alves,et al.  Comparing systemic properties of ensembles of biological networks by graphical and statistical methods , 2000, Bioinform..

[31]  Carl D. Meyer,et al.  Matrix Analysis and Applied Linear Algebra , 2000 .

[32]  Bernhard O. Palsson,et al.  Metabolic flux balance analysis and the in silico analysis of Escherichia coli K-12 gene deletions , 2000, BMC Bioinformatics.

[33]  A. Burgard,et al.  Probing the performance limits of the Escherichia coli metabolic network subject to gene additions or deletions. , 2001, Biotechnology and bioengineering.

[34]  B. Palsson,et al.  In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data , 2001, Nature Biotechnology.

[35]  Bernhard O. Palsson,et al.  Dynamic simulation of the human red blood cell metabolic network , 2001, Bioinform..

[36]  Farren J. Isaacs,et al.  Computational studies of gene regulatory networks: in numero molecular biology , 2001, Nature Reviews Genetics.

[37]  A. Burgard,et al.  Minimal Reaction Sets for Escherichia coli Metabolism under Different Growth Requirements and Uptake Environments , 2001, Biotechnology progress.

[38]  Jeff Hasty,et al.  Engineered gene circuits , 2002, Nature.

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

[40]  Bernhard O Palsson,et al.  Extreme pathway analysis of human red blood cell metabolism. , 2002, Biophysical journal.

[41]  Neema Jamshidi,et al.  In silico model-driven assessment of the effects of single nucleotide polymorphisms (SNPs) on human red blood cell metabolism. , 2002, Genome research.

[42]  B. Palsson,et al.  Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth , 2002, Nature.

[43]  B. Palsson,et al.  Metabolic modelling of microbes: the flux-balance approach. , 2002, Environmental microbiology.

[44]  G. Church,et al.  Genome-Scale Metabolic Model of Helicobacter pylori 26695 , 2002, Journal of bacteriology.

[45]  Jason A. Papin,et al.  Extreme pathway lengths and reaction participation in genome-scale metabolic networks. , 2002, Genome research.

[46]  Neema Jamshidi,et al.  Description and analysis of metabolic connectivity and dynamics in the human red blood cell. , 2002, Biophysical journal.

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

[48]  Bernhard O Palsson,et al.  Network-based analysis of metabolic regulation in the human red blood cell. , 2003, Journal of theoretical biology.

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

[50]  B. Palsson,et al.  Large-scale evaluation of in silico gene deletions in Saccharomyces cerevisiae. , 2003, Omics : a journal of integrative biology.

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

[52]  P. Kuchel,et al.  Modelling metabolism with Mathematica : detailed examples including erythrocyte metabolism , 2003 .

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

[54]  B. Palsson,et al.  Thirteen Years of Building Constraint-Based In Silico Models of Escherichia coli , 2003, Journal of bacteriology.

[55]  Michael A Savageau,et al.  Quantitative evolutionary design of glucose 6-phosphate dehydrogenase expression in human erythrocytes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[56]  Jason A. Papin,et al.  Genome-scale microbial in silico models: the constraints-based approach. , 2003, Trends in biotechnology.

[57]  B. Palsson,et al.  Identifying constraints that govern cell behavior: a key to converting conceptual to computational models in biology? , 2003, Biotechnology and bioengineering.

[58]  Jason A. Papin,et al.  Metabolic pathways in the post-genome era. , 2003, Trends in biochemical sciences.

[59]  A. Barabasi,et al.  Global organization of metabolic fluxes in the bacterium Escherichia coli , 2004, Nature.

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

[61]  H. J. Greenberg,et al.  Monte Carlo sampling can be used to determine the size and shape of the steady-state flux space. , 2004, Journal of theoretical biology.

[62]  Virgilio L. Lew,et al.  Volume, pH, and ion-content regulation in human red cells: Analysis of transient behavior with an integrated model , 2005, The Journal of Membrane Biology.