Emerging ensembles of kinetic parameters to characterize observed metabolic phenotypes

BackgroundDetermining the value of kinetic constants for a metabolic system in the exact physiological conditions is an extremely hard task. However, this kind of information is of pivotal relevance to effectively simulate a biological phenomenon as complex as metabolism.ResultsTo overcome this issue, we propose to investigate emerging properties of ensembles of sets of kinetic constants leading to the biological readout observed in different experimental conditions. To this aim, we exploit information retrievable from constraint-based analyses (i.e. metabolic flux distributions at steady state) with the goal to generate feasible values for kinetic constants exploiting the mass action law. The sets retrieved from the previous step will be used to parametrize a mechanistic model whose simulation will be performed to reconstruct the dynamics of the system (until reaching the metabolic steady state) for each experimental condition. Every parametrization that is in accordance with the expected metabolic phenotype is collected in an ensemble whose features are analyzed to determine the emergence of properties of a phenotype. In this work we apply the proposed approach to identify ensembles of kinetic parameters for five metabolic phenotypes of E. Coli, by analyzing five different experimental conditions associated with the ECC2comp model recently published by Hädicke and collaborators.ConclusionsOur results suggest that the parameter values of just few reactions are responsible for the emergence of a metabolic phenotype. Notably, in contrast with constraint-based approaches such as Flux Balance Analysis, the methodology used in this paper does not require to assume that metabolism is optimizing towards a specific goal.

[1]  L. Alberghina,et al.  Systems Biology: Definitions and Perspectives , 2005 .

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

[3]  Jing Zhao,et al.  Complex networks theory for analyzing metabolic networks , 2006, q-bio/0603015.

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

[5]  Giancarlo Mauri,et al.  Constraining Mechanism Based Simulations to Identify Ensembles of Parametrizations to Characterize Metabolic Features , 2016, CIBB.

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

[7]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[8]  E. Klipp,et al.  Modelling signalling pathways-A yeast approach , 2005 .

[9]  Jong Myoung Park,et al.  Constraints-based genome-scale metabolic simulation for systems metabolic engineering. , 2009, Biotechnology advances.

[10]  S. Henry,et al.  Revising the Representation of Fatty Acid, Glycerolipid, and Glycerophospholipid Metabolism in the Consensus Model of Yeast Metabolism. , 2013, Industrial biotechnology.

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

[12]  Ram Krishnamurthy,et al.  ECMDB: The E. coli Metabolome Database , 2012, Nucleic Acids Res..

[13]  Marco S. Nobile,et al.  Computational Strategies for a System-Level Understanding of Metabolism , 2014, Metabolites.

[14]  Giancarlo Mauri,et al.  An ensemble evolutionary constraint-based approach to understand the emergence of metabolic phenotypes , 2014, Natural Computing.

[15]  Erwin P. Gianchandani,et al.  The application of flux balance analysis in systems biology , 2010, Wiley interdisciplinary reviews. Systems biology and medicine.

[16]  Giancarlo Mauri,et al.  Zooming-in on cancer metabolic rewiring with tissue specific constraint-based models , 2016, Comput. Biol. Chem..

[17]  Jan-Hendrik S. Hofmeyr,et al.  Modelling cellular systems with PySCeS , 2005, Bioinform..

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

[19]  Giancarlo Mauri,et al.  cuTauLeaping: A GPU-Powered Tau-Leaping Stochastic Simulator for Massive Parallel Analyses of Biological Systems , 2014, PloS one.

[20]  Giancarlo Mauri,et al.  GPU-accelerated simulations of mass-action kinetics models with cupSODA , 2014, The Journal of Supercomputing.

[21]  Christopher R. Myers,et al.  Universally Sloppy Parameter Sensitivities in Systems Biology Models , 2007, PLoS Comput. Biol..

[22]  L. Petzold Automatic Selection of Methods for Solving Stiff and Nonstiff Systems of Ordinary Differential Equations , 1983 .

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

[24]  M. Reuss,et al.  In vivo analysis of metabolic dynamics in Saccharomyces cerevisiae : I. Experimental observations. , 1997, Biotechnology and bioengineering.

[25]  S. Klamt,et al.  EColiCore2: a reference network model of the central metabolism of Escherichia coli and relationships to its genome-scale parent model , 2017, Scientific Reports.