Probabilistic thermodynamic analysis of metabolic networks

Abstract Motivation Random sampling of metabolic fluxes can provide a comprehensive description of the capabilities of a metabolic network. However, current sampling approaches do not model thermodynamics explicitly, leading to inaccurate predictions of an organism’s potential or actual metabolic operations. Results We present a probabilistic framework combining thermodynamic quantities with steady-state flux constraints to analyze the properties of a metabolic network. It includes methods for probabilistic metabolic optimization and for joint sampling of thermodynamic and flux spaces. Applied to a model of Escherichia coli, we use the methods to reveal known and novel mechanisms of substrate channeling, and to accurately predict reaction directions and metabolite concentrations. Interestingly, predicted flux distributions are multimodal, leading to discrete hypotheses on E.coli’s metabolic capabilities. Availability and implementation Python and MATLAB packages available at https://gitlab.com/csb.ethz/pta. Supplementary information Supplementary data are available at Bioinformatics online.

[1]  Adam M. Feist,et al.  iML1515, a knowledgebase that computes Escherichia coli traits , 2017, Nature Biotechnology.

[2]  Robert N. Goldberg,et al.  Thermodynamics of enzyme-catalyzed reactions - a database for quantitative biochemistry , 2004, Bioinform..

[3]  V. Hatzimanikatis,et al.  Thermodynamics-based metabolic flux analysis. , 2007, Biophysical journal.

[4]  Urs von Stockar,et al.  Influence of uncertainties in pH, pMg, activity coefficients, metabolite concentrations, and other factors on the analysis of the thermodynamic feasibility of metabolic pathways , 2009, Biotechnology and bioengineering.

[5]  Richard C. Larson,et al.  Model Building in Mathematical Programming , 1979 .

[6]  D. Rakus,et al.  Interaction between muscle aldolase and muscle fructose 1,6-bisphosphatase results in the substrate channeling. , 2004, Biochemistry.

[7]  Daniel C. Zielinski,et al.  Temperature-dependent estimation of Gibbs energies using an updated group contribution method , 2018, bioRxiv.

[8]  Hongwu Ma,et al.  Metabolic engineering of Escherichia coli for poly(3-hydroxybutyrate) production via threonine bypass , 2015, Microbial Cell Factories.

[9]  Michael L. Mavrovouniotis,et al.  Identification of Localized and Distributed Bottlenecks in Metabolic Pathways , 1993, ISMB.

[10]  Santosh S. Vempala,et al.  CHRR: coordinate hit-and-run with rounding for uniform sampling of constraint-based models , 2017, Bioinform..

[11]  Anne Richelle,et al.  Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0 , 2019, Nature Protocols.

[12]  T. Shlomi,et al.  Metabolite concentrations, fluxes, and free energies imply efficient enzyme usage , 2016, Nature chemical biology.

[13]  Ronan M. T. Fleming,et al.  Consistent Estimation of Gibbs Energy Using Component Contributions , 2013, PLoS Comput. Biol..

[14]  Yinjie J. Tang,et al.  Comparative studies of glycolytic pathways and channeling under in vitro and in vivo modes , 2018, AIChE Journal.

[15]  Lars K. Nielsen,et al.  ll-ACHRB: a scalable algorithm for sampling the feasible solution space of metabolic networks , 2016, Bioinform..

[16]  Juho Rousu,et al.  Bayesian metabolic flux analysis reveals intracellular flux couplings , 2018, Bioinform..

[17]  Robert L. Smith,et al.  Hit-and-Run Algorithms for Generating Multivariate Distributions , 1993, Math. Oper. Res..

[18]  U. Sauer,et al.  Pseudo-transition Analysis Identifies the Key Regulators of Dynamic Metabolic Adaptations from Steady-State Data. , 2015, Cell systems.

[19]  Florian Jarre,et al.  CycleFreeFlux: efficient removal of thermodynamically infeasible loops from flux distributions , 2015, Bioinform..

[20]  Bernhard O. Palsson,et al.  Thermodynamic favorability and pathway yield as evolutionary tradeoffs in biosynthetic pathway choice , 2018, Proceedings of the National Academy of Sciences.

[21]  Ignacio Fita,et al.  A novel two-domain architecture within the amino acid kinase enzyme family revealed by the crystal structure of Escherichia coli glutamate 5-kinase. , 2007, Journal of molecular biology.

[22]  Matthew D. Jankowski,et al.  Group contribution method for thermodynamic analysis of complex metabolic networks. , 2008, Biophysical journal.

[23]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[24]  Ronan M. T. Fleming,et al.  Quantitative assignment of reaction directionality in a multicompartmental human metabolic reconstruction. , 2012, Biophysical journal.

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

[26]  Jörg Stelling,et al.  Large-scale computation of elementary flux modes with bit pattern trees , 2008, Bioinform..

[27]  B. Palsson,et al.  Theory for the systemic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective. , 2000, Journal of theoretical biology.

[28]  John K Kruschke,et al.  Bayesian data analysis. , 2010, Wiley interdisciplinary reviews. Cognitive science.

[29]  Steffen Klamt,et al.  An algorithm for the reduction of genome-scale metabolic network models to meaningful core models , 2015, BMC Systems Biology.

[30]  Cuiqing Ma,et al.  Coupling between d-3-phosphoglycerate dehydrogenase and d-2-hydroxyglutarate dehydrogenase drives bacterial l-serine synthesis , 2017, Proceedings of the National Academy of Sciences.

[31]  Joshua A. Lerman,et al.  COBRApy: COnstraints-Based Reconstruction and Analysis for Python , 2013, BMC Systems Biology.

[32]  Vassily Hatzimanikatis,et al.  Bioenergetics-based modeling of Plasmodium falciparum metabolism reveals its essential genes, nutritional requirements, and thermodynamic bottlenecks , 2017, PLoS Comput. Biol..

[33]  Paul A. Jensen,et al.  gapsplit: Efficient random sampling for non-convex constraint-based models , 2019, bioRxiv.

[34]  Jörg Stelling,et al.  Integrating -omics data into genome-scale metabolic network models: principles and challenges. , 2018, Essays in biochemistry.

[35]  Michael Binns,et al.  Sampling with poling-based flux balance analysis: optimal versus sub-optimal flux space analysis of Actinobacillus succinogenes , 2015, BMC Bioinformatics.