Probabilistic thermodynamic analysis of metabolic networks

Random sampling of metabolic fluxes can provide an unbiased 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. 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 E. 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. C++ source code with MATLAB interface available at https://gitlab.com/csb.ethz/pta.

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

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

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

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

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

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

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

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

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

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

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

[12]  H. P. Williams,et al.  Model Building in Mathematical Programming , 1979 .

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

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

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

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

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

[18]  Cedric Chauve,et al.  Evolution of genes neighborhood within reconciled phylogenies: an ensemble approach , 2015, bioRxiv.

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

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

[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]  Santosh S. Vempala,et al.  CHRR: coordinate hit-and-run with rounding for uniform sampling of constraint-based models , 2017, Bioinform..

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

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

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

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

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

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

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

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

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

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

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