An analytic approximation of the feasible space of metabolic networks

Assuming a steady-state condition within a cell, metabolic fluxes satisfy an underdetermined linear system of stoichiometric equations. Characterizing the space of fluxes that satisfy such equations along with given bounds (and possibly additional relevant constraints) is considered of utmost importance for the understanding of cellular metabolism. Extreme values for each individual flux can be computed with linear programming (as flux balance analysis), and their marginal distributions can be approximately computed with Monte Carlo sampling. Here we present an approximate analytic method for the latter task based on expectation propagation equations that does not involve sampling and can achieve much better predictions than other existing analytic methods. The method is iterative, and its computation time is dominated by one matrix inversion per iteration. With respect to sampling, we show through extensive simulation that it has some advantages including computation time, and the ability to efficiently fix empirically estimated distributions of fluxes.

[1]  Enzo Marinari,et al.  A Scalable Algorithm to Explore the Gibbs Energy Landscape of Genome-Scale Metabolic Networks , 2012, PLoS Comput. Biol..

[2]  Alfredo Braunstein,et al.  The space of feasible solutions in metabolic networks , 2008 .

[3]  Rick L. Stevens,et al.  High-throughput generation, optimization and analysis of genome-scale metabolic models , 2010, Nature Biotechnology.

[4]  Ronan M. T. Fleming,et al.  Reconstruction and Use of Microbial Metabolic Networks: the Core Escherichia coli Metabolic Model as an Educational Guide. , 2010, EcoSal Plus.

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

[6]  E. Jaynes.  Information Theory and Statistical Mechanics , 1957 .

[7]  Ole Winther,et al.  Gaussian Processes for Classification: Mean-Field Algorithms , 2000, Neural Computation.

[8]  Tom Minka,et al.  Expectation Propagation for approximate Bayesian inference , 2001, UAI.

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

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

[11]  Jan Schellenberger,et al.  Use of Randomized Sampling for Analysis of Metabolic Networks* , 2009, Journal of Biological Chemistry.

[12]  Fabrizio Capuani,et al.  Growth against entropy in bacterial metabolism: the phenotypic trade-off behind empirical growth rate distributions in E. coli , 2016, Physical biology.

[13]  Francesco Alessandro Massucci,et al.  A Novel Methodology to Estimate Metabolic Flux Distributions in Constraint-Based Models , 2013, Metabolites.

[14]  Philip Miller,et al.  BiGG Models: A platform for integrating, standardizing and sharing genome-scale models , 2015, Nucleic Acids Res..

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

[16]  Monica L. Mo,et al.  Global reconstruction of the human metabolic network based on genomic and bibliomic data , 2007, Proceedings of the National Academy of Sciences.

[17]  B. Palsson,et al.  Uniform sampling of steady-state flux spaces: means to design experiments and to interpret enzymopathies. , 2004, Biophysical journal.

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

[19]  Roberto Mulet,et al.  Fast inference of ill-posed problems within a convex space , 2016, ArXiv.

[20]  V. Turchin.  On the Computation of Multidimensional Integrals by the Monte-Carlo Method , 1971 .

[21]  Francesco Alessandro Massucci,et al.  A weighted belief-propagation algorithm for estimating volume-related properties of random polytopes , 2012, 1208.1295.

[22]  Martin E. Dyer,et al.  On the Complexity of Computing the Volume of a Polyhedron , 1988, SIAM J. Comput..

[23]  Alfredo Braunstein,et al.  Estimating the size of the solution space of metabolic networks , 2007, BMC Bioinformatics.

[24]  M. Opper,et al.  Adaptive and self-averaging Thouless-Anderson-Palmer mean-field theory for probabilistic modeling. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[26]  B. Palsson,et al.  Large-scale in silico modeling of metabolic interactions between cell types in the human brain , 2010, Nature Biotechnology.

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

[28]  S. Stenholm.  Information, Physics and Computation, by Marc Mézard and Andrea Montanari , 2010 .

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

[30]  John T. Sauls,et al.  Cell-Size Control and Homeostasis in Bacteria , 2015, Current Biology.

[31]  Robert L. Smith,et al.  Efficient Monte Carlo Procedures for Generating Points Uniformly Distributed over Bounded Regions , 1984, Oper. Res..

[32]  Matteo Mori,et al.  Uniform Sampling of Steady States in Metabolic Networks: Heterogeneous Scales and Rounding , 2013, PloS one.

[33]  M. Huynen,et al.  optGpSampler: An Improved Tool for Uniformly Sampling the Solution-Space of Genome-Scale Metabolic Networks , 2014, PloS one.

[34]  Philippe Nghe,et al.  Individuality and universality in the growth-division laws of single E. coli cells. , 2014, Physical review. E.