Robust Analysis of Fluxes in Genome-Scale Metabolic Pathways

Constraint-based optimization, such as flux balance analysis (FBA), has become a standard systems-biology computational method to study cellular metabolisms that are assumed to be in a steady state of optimal growth. The methods are based on optimization while assuming (i) equilibrium of a linear system of ordinary differential equations, and (ii) deterministic data. However, the steady-state assumption is biologically imperfect, and several key stoichiometric coefficients are experimentally inferred from situations of inherent variation. We propose an approach that explicitly acknowledges heterogeneity and conducts a robust analysis of metabolic pathways (RAMP). The basic assumption of steady state is relaxed, and we model the innate heterogeneity of cells probabilistically. Our mathematical study of the stochastic problem shows that FBA is a limiting case of our RAMP method. Moreover, RAMP has the properties that: A) metabolic states are (Lipschitz) continuous with regards to the probabilistic modeling parameters, B) convergent metabolic states are solutions to the deterministic FBA paradigm as the stochastic elements dissipate, and C) RAMP can identify biologically tolerable diversity of a metabolic network in an optimized culture. We benchmark RAMP against traditional FBA on genome-scale metabolic reconstructed models of E. coli, calculating essential genes and comparing with experimental flux data.

[1]  S. Lee,et al.  Metabolic engineering of Escherichia coli for the production of l-valine based on transcriptome analysis and in silico gene knockout simulation , 2007, Proceedings of the National Academy of Sciences.

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

[3]  Allen Holder,et al.  An Introduction to Systems Biology for Mathematical Programmers , 2008 .

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

[5]  Costas D. Maranas,et al.  Optimization Methods in Metabolic Networks: Maranas/Optimization Methods in Metabolic Networks , 2016 .

[6]  Albert-László Barabási,et al.  The Activity Reaction Core and Plasticity of Metabolic Networks , 2005, PLoS Comput. Biol..

[7]  Bernhard O. Palsson,et al.  Context-Specific Metabolic Networks Are Consistent with Experiments , 2008, PLoS Comput. Biol..

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

[9]  B. Palsson.  Systems Biology: Constraint-based Reconstruction and Analysis , 2015 .

[10]  Arkadi Nemirovski,et al.  Robust Convex Optimization , 1998, Math. Oper. Res..

[11]  Ronan M. T. Fleming,et al.  Do genome-scale models need exact solvers or clearer standards? , 2015, Molecular systems biology.

[12]  S. Henderson,et al.  Robust optimization for intensity modulated radiation therapy treatment planning under uncertainty , 2005, Physics in medicine and biology.

[13]  Ali Navid,et al.  Genome-level transcription data of Yersinia pestis analyzed with a New metabolic constraint-based approach , 2012, BMC Systems Biology.

[14]  B. Palsson,et al.  Biochemical production capabilities of escherichia coli , 1993, Biotechnology and bioengineering.

[15]  B. Palsson,et al.  Characterization of Metabolism in the Fe(III)-Reducing Organism Geobacter sulfurreducens by Constraint-Based Modeling , 2006, Applied and Environmental Microbiology.

[16]  Zachary A. King,et al.  Constraint-based models predict metabolic and associated cellular functions , 2014, Nature Reviews Genetics.

[17]  B. Palsson,et al.  A protocol for generating a high-quality genome-scale metabolic reconstruction , 2010 .

[18]  S. Lee,et al.  Systems metabolic engineering of Escherichia coli for L-threonine production , 2007, Molecular systems biology.

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

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

[21]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[22]  M. Ataman,et al.  Heading in the right direction: thermodynamics-based network analysis and pathway engineering. , 2015, Current opinion in biotechnology.

[23]  Michael M. Zavlanos,et al.  Robust flux balance analysis of metabolic networks , 2011, Proceedings of the 2011 American Control Conference.

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

[25]  C. Maranas,et al.  Succinate Overproduction: A Case Study of Computational Strain Design Using a Comprehensive Escherichia coli Kinetic Model , 2015, Front. Bioeng. Biotechnol..

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

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

[28]  U. Sauer,et al.  Impact of Global Transcriptional Regulation by ArcA, ArcB, Cra, Crp, Cya, Fnr, and Mlc on Glucose Catabolism in Escherichia coli , 2005, Journal of bacteriology.

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

[30]  Eytan Ruppin,et al.  A computational study of the Warburg effect identifies metabolic targets inhibiting cancer migration , 2014, Molecular Systems Biology.

[31]  Bonnie Berger,et al.  An exact arithmetic toolbox for a consistent and reproducible structural analysis of metabolic network models , 2014, Nature Communications.

[32]  Tomer Shlomi,et al.  Prediction of Microbial Growth Rate versus Biomass Yield by a Metabolic Network with Kinetic Parameters , 2012, PLoS Comput. Biol..

[33]  Arkadi Nemirovski,et al.  Robust solutions of uncertain linear programs , 1999, Oper. Res. Lett..

[34]  Arkadi Nemirovski,et al.  Robust Truss Topology Design via Semidefinite Programming , 1997, SIAM J. Optim..

[35]  A. Burgard,et al.  Metabolic engineering of Escherichia coli for direct production of 1,4-butanediol. , 2011, Nature chemical biology.

[36]  Oscar P. Kuipers,et al.  Phenotypic variation in bacteria: the role of feedback regulation , 2006, Nature Reviews Microbiology.

[37]  B. Palsson,et al.  Regulation of gene expression in flux balance models of metabolism. , 2001, Journal of theoretical biology.

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

[39]  Tetsuya Yomo,et al.  Dynamic clustering of bacterial population , 1994 .

[40]  Annik Nanchen,et al.  Nonlinear Dependency of Intracellular Fluxes on Growth Rate in Miniaturized Continuous Cultures of Escherichia coli , 2006, Applied and Environmental Microbiology.

[41]  F. Neidhardt,et al.  Physiology of the bacterial cell : a molecular approach , 1990 .

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

[43]  Adam M. Feist,et al.  A comprehensive genome-scale reconstruction of Escherichia coli metabolism—2011 , 2011, Molecular systems biology.

[44]  W. Martin,et al.  Cell individuality : the bistability of competence development , 2005 .

[45]  Bonnie Berger,et al.  Reply to “Do genome-scale models need exact solvers or clearer standards?” , 2015, Molecular systems biology.

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

[47]  B. Palsson,et al.  Constraining the metabolic genotype–phenotype relationship using a phylogeny of in silico methods , 2012, Nature Reviews Microbiology.

[48]  M. A. de Menezes,et al.  Intracellular crowding defines the mode and sequence of substrate uptake by Escherichia coli and constrains its metabolic activity , 2007, Proceedings of the National Academy of Sciences.