The Ability of Flux Balance Analysis to Predict Evolution of Central Metabolism Scales with the Initial Distance to the Optimum

The most powerful genome-scale framework to model metabolism, flux balance analysis (FBA), is an evolutionary optimality model. It hypothesizes selection upon a proposed optimality criterion in order to predict the set of internal fluxes that would maximize fitness. Here we present a direct test of the optimality assumption underlying FBA by comparing the central metabolic fluxes predicted by multiple criteria to changes measurable by a 13C-labeling method for experimentally-evolved strains. We considered datasets for three Escherichia coli evolution experiments that varied in their length, consistency of environment, and initial optimality. For ten populations that were evolved for 50,000 generations in glucose minimal medium, we observed modest changes in relative fluxes that led to small, but significant decreases in optimality and increased the distance to the predicted optimal flux distribution. In contrast, seven populations evolved on the poor substrate lactate for 900 generations collectively became more optimal and had flux distributions that moved toward predictions. For three pairs of central metabolic knockouts evolved on glucose for 600–800 generations, there was a balance between cases where optimality and flux patterns moved toward or away from FBA predictions. Despite this variation in predictability of changes in central metabolism, two generalities emerged. First, improved growth largely derived from evolved increases in the rate of substrate use. Second, FBA predictions bore out well for the two experiments initiated with ancestors with relatively sub-optimal yield, whereas those begun already quite optimal tended to move somewhat away from predictions. These findings suggest that the tradeoff between rate and yield is surprisingly modest. The observed positive correlation between rate and yield when adaptation initiated further from the optimum resulted in the ability of FBA to use stoichiometric constraints to predict the evolution of metabolism despite selection for rate.

[1]  J. Collins,et al.  DIVERSITY-BASED, MODEL-GUIDED CONSTRUCTION OF SYNTHETIC GENE NETWORKS WITH PREDICTED FUNCTIONS , 2009, Nature Biotechnology.

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

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

[4]  Richard E. Lenski,et al.  Experimental Tests for an Evolutionary Trade‐Off between Growth Rate and Yield in E. coli , 2006, The American Naturalist.

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

[6]  U. Sauer,et al.  Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli , 2007, Molecular systems biology.

[7]  Gabriela Kalna,et al.  Haem oxygenase is synthetically lethal with the tumour suppressor fumarate hydratase , 2011, Nature.

[8]  Jeffrey E. Barrick,et al.  Genome evolution and adaptation in a long-term experiment with Escherichia coli , 2009, Nature.

[9]  S. Gould,et al.  The spandrels of San Marco and the Panglossian paradigm: a critique of the adaptationist programme , 1979, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[10]  U. Sauer,et al.  Multidimensional Optimality of Microbial Metabolism , 2012, Science.

[11]  Tsvi Tlusty,et al.  Cross-species analysis traces adaptation of Rubisco toward optimality in a low-dimensional landscape , 2010, Proceedings of the National Academy of Sciences.

[12]  J. Reed,et al.  Synergy between (13)C-metabolic flux analysis and flux balance analysis for understanding metabolic adaptation to anaerobiosis in E. coli. , 2011, Metabolic engineering.

[13]  Richard E. Lenski,et al.  Evolution of competitive fitness in experimental populations of E. coli: What makes one genotype a better competitor than another? , 1998, Antonie van Leeuwenhoek.

[14]  Wim Soetaert,et al.  Effect of iclR and arcA knockouts on biomass formation and metabolic fluxes in Escherichia coli K12 and its implications on understanding the metabolism of Escherichia coli BL21 (DE3) , 2011, BMC Microbiology.

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

[16]  R. Lenski,et al.  Long-Term Experimental Evolution in Escherichia coli. I. Adaptation and Divergence During 2,000 Generations , 1991, The American Naturalist.

[17]  William R. Harcombe,et al.  Multiple long-term, experimentally-evolved populations of Escherichia coli acquire dependence upon citrate as an iron chelator for optimal growth on glucose , 2012, BMC Evolutionary Biology.

[18]  Balázs Papp,et al.  Systems-biology approaches for predicting genomic evolution , 2011, Nature Reviews Genetics.

[19]  Achille Messac,et al.  Integrated Energy and Flux Balance Based Multiobjective Framework for Large-Scale Metabolic Networks , 2007, Annals of Biomedical Engineering.

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

[21]  Richard E. Lenski,et al.  Long-Term Experimental Evolution in Escherichia coli. II. Changes in Life-History Traits During Adaptation to a Seasonal Environment , 1994, The American Naturalist.

[22]  Andrew R. Joyce,et al.  Metabolic Characterization of Escherichia coli Strains Adapted to Growth on Lactate , 2007, Applied and Environmental Microbiology.

[23]  Bernhard O. Palsson,et al.  Constraint-based analysis of metabolic capacity of Salmonella typhimurium during host-pathogen interaction , 2009, BMC Systems Biology.

[24]  J. Liao,et al.  Metabolic ensemble modeling for strain engineers , 2012, Biotechnology journal.

[25]  Matthew D. Rolfe,et al.  Network analysis of the transcriptional pattern of young and old cells of Escherichia coli during lag phase , 2009, BMC Systems Biology.

[26]  A. Kempf,et al.  Optimization Models and the Structure of the Genetic Code , 2009, Journal of Molecular Evolution.

[27]  Jay D Keasling,et al.  Model-Driven Engineering of RNA Devices to Quantitatively Program Gene Expression , 2011, Science.

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

[29]  Adam M. Feist,et al.  The biomass objective function. , 2010, Current opinion in microbiology.

[30]  U. Sauer,et al.  Determination of metabolic flux ratios from 13C-experiments and gas chromatography-mass spectrometry data: protocol and principles. , 2007, Methods in molecular biology.

[31]  Tom M. Conrad,et al.  Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models , 2010, Molecular systems biology.

[32]  U. Alon,et al.  Optimality and evolutionary tuning of the expression level of a protein , 2005, Nature.

[33]  Bernhard O Palsson,et al.  Latent Pathway Activation and Increased Pathway Capacity Enable Escherichia coli Adaptation to Loss of Key Metabolic Enzymes* , 2006, Journal of Biological Chemistry.

[34]  Adam M. Feist,et al.  A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information , 2007, Molecular systems biology.

[35]  B. Palsson,et al.  Parallel adaptive evolution cultures of Escherichia coli lead to convergent growth phenotypes with different gene expression states. , 2005, Genome research.

[36]  Bas Teusink,et al.  Understanding the Adaptive Growth Strategy of Lactobacillus plantarum by In Silico Optimisation , 2009, PLoS Comput. Biol..

[37]  James J Bull,et al.  TESTING OPTIMALITY WITH EXPERIMENTAL EVOLUTION: LYSIS TIME IN A BACTERIOPHAGE , 2007, Evolution; international journal of organic evolution.

[38]  U. Sauer,et al.  Metabolic Flux Responses to Pyruvate Kinase Knockout in Escherichia coli , 2002, Journal of bacteriology.

[39]  H. Kacser,et al.  The control of flux. , 1995, Biochemical Society transactions.

[40]  Richard E. Lenski,et al.  Long-Term Experimental Evolution in Escherichia coli. XIII. Phylogenetic History of a Balanced Polymorphism , 2005, Journal of Molecular Evolution.

[41]  Barbara M. Bakker,et al.  Unraveling the complexity of flux regulation: A new method demonstrated for nutrient starvation in Saccharomyces cerevisiae , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[42]  R. Lenski,et al.  Historical contingency and the evolution of a key innovation in an experimental population of Escherichia coli , 2008 .