A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data

Ali Khodayari, Ali R. Zomorrodi, James C. Liao, Costas D. Maranas

Research output: Contribution to journalArticle

78 Citations (Scopus)

Abstract

In contrast to stoichiometric-based models, the development of large-scale kinetic models of metabolism has been hindered by the challenge of identifying kinetic parameter values and kinetic rate laws applicable to a wide range of environmental and/or genetic perturbations. The recently introduced ensemble modeling (EM) procedure provides a promising remedy to address these challenges by decomposing metabolic reactions into elementary reaction steps and incorporating all phenotypic observations, upon perturbation, in its model parameterization scheme. Here, we present a kinetic model of Escherichia coli core metabolism that satisfies the fluxomic data for wild-type and seven mutant strains by making use of the EM concepts. This model encompasses 138 reactions, 93 metabolites and 60 substrate-level regulatory interactions accounting for glycolysis/gluconeogenesis, pentose phosphate pathway, TCA cycle, major pyruvate metabolism, anaplerotic reactions and a number of reactions in other parts of the metabolism. Parameterization is performed using a formal optimization approach that minimizes the discrepancies between model predictions and flux measurements. The predicted fluxes by the model are within the uncertainty range of experimental flux data for 78% of the reactions (with measured fluxes) for both the wild-type and seven mutant strains. The remaining flux predictions are mostly within three standard deviations of reported ranges. Converting the EM-based parameters into a Michaelis-Menten equivalent formalism revealed that 35% of Km and 77% of kcat parameters are within uncertainty range of the literature-reported values. The predicted metabolite concentrations by the model are also within uncertainty ranges of metabolomic data for 68% of the metabolites. A leave-one-out cross-validation test to evaluate the flux prediction performance of the model showed that metabolic fluxes for the mutants located in the proximity of mutations used for training the model can be predicted more accurately. The constructed model and the parameterization procedure presented in this study pave the way for the construction of larger-scale kinetic models with more narrowly distributed parameter values as new metabolomic/fluxomic data sets are becoming available for E. coli and other organisms.

Original languageEnglish (US)
Pages (from-to)50-62
Number of pages13
JournalMetabolic engineering
Volume25
DOIs
StatePublished - Jan 1 2014

Fingerprint

Enzyme kinetics
Metabolism
Escherichia coli
Fluxes
Uncertainty
Metabolomics
Metabolites
Parameterization
Pentose Phosphate Pathway
Gluconeogenesis
Glycolysis
Pyruvic Acid
Mutation
Pentoses
Kinetics
Kinetic parameters
Phosphates

All Science Journal Classification (ASJC) codes

  • Biotechnology
  • Bioengineering
  • Applied Microbiology and Biotechnology

Cite this

@article{aa7297c888164eefb9206e4f9e4a1562,
title = "A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data",
abstract = "In contrast to stoichiometric-based models, the development of large-scale kinetic models of metabolism has been hindered by the challenge of identifying kinetic parameter values and kinetic rate laws applicable to a wide range of environmental and/or genetic perturbations. The recently introduced ensemble modeling (EM) procedure provides a promising remedy to address these challenges by decomposing metabolic reactions into elementary reaction steps and incorporating all phenotypic observations, upon perturbation, in its model parameterization scheme. Here, we present a kinetic model of Escherichia coli core metabolism that satisfies the fluxomic data for wild-type and seven mutant strains by making use of the EM concepts. This model encompasses 138 reactions, 93 metabolites and 60 substrate-level regulatory interactions accounting for glycolysis/gluconeogenesis, pentose phosphate pathway, TCA cycle, major pyruvate metabolism, anaplerotic reactions and a number of reactions in other parts of the metabolism. Parameterization is performed using a formal optimization approach that minimizes the discrepancies between model predictions and flux measurements. The predicted fluxes by the model are within the uncertainty range of experimental flux data for 78{\%} of the reactions (with measured fluxes) for both the wild-type and seven mutant strains. The remaining flux predictions are mostly within three standard deviations of reported ranges. Converting the EM-based parameters into a Michaelis-Menten equivalent formalism revealed that 35{\%} of Km and 77{\%} of kcat parameters are within uncertainty range of the literature-reported values. The predicted metabolite concentrations by the model are also within uncertainty ranges of metabolomic data for 68{\%} of the metabolites. A leave-one-out cross-validation test to evaluate the flux prediction performance of the model showed that metabolic fluxes for the mutants located in the proximity of mutations used for training the model can be predicted more accurately. The constructed model and the parameterization procedure presented in this study pave the way for the construction of larger-scale kinetic models with more narrowly distributed parameter values as new metabolomic/fluxomic data sets are becoming available for E. coli and other organisms.",
author = "Ali Khodayari and Zomorrodi, {Ali R.} and Liao, {James C.} and Maranas, {Costas D.}",
year = "2014",
month = "1",
day = "1",
doi = "10.1016/j.ymben.2014.05.014",
language = "English (US)",
volume = "25",
pages = "50--62",
journal = "Metabolic Engineering",
issn = "1096-7176",
publisher = "Academic Press Inc.",

}

A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data. / Khodayari, Ali; Zomorrodi, Ali R.; Liao, James C.; Maranas, Costas D.

In: Metabolic engineering, Vol. 25, 01.01.2014, p. 50-62.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data

AU - Khodayari, Ali

AU - Zomorrodi, Ali R.

AU - Liao, James C.

AU - Maranas, Costas D.

PY - 2014/1/1

Y1 - 2014/1/1

N2 - In contrast to stoichiometric-based models, the development of large-scale kinetic models of metabolism has been hindered by the challenge of identifying kinetic parameter values and kinetic rate laws applicable to a wide range of environmental and/or genetic perturbations. The recently introduced ensemble modeling (EM) procedure provides a promising remedy to address these challenges by decomposing metabolic reactions into elementary reaction steps and incorporating all phenotypic observations, upon perturbation, in its model parameterization scheme. Here, we present a kinetic model of Escherichia coli core metabolism that satisfies the fluxomic data for wild-type and seven mutant strains by making use of the EM concepts. This model encompasses 138 reactions, 93 metabolites and 60 substrate-level regulatory interactions accounting for glycolysis/gluconeogenesis, pentose phosphate pathway, TCA cycle, major pyruvate metabolism, anaplerotic reactions and a number of reactions in other parts of the metabolism. Parameterization is performed using a formal optimization approach that minimizes the discrepancies between model predictions and flux measurements. The predicted fluxes by the model are within the uncertainty range of experimental flux data for 78% of the reactions (with measured fluxes) for both the wild-type and seven mutant strains. The remaining flux predictions are mostly within three standard deviations of reported ranges. Converting the EM-based parameters into a Michaelis-Menten equivalent formalism revealed that 35% of Km and 77% of kcat parameters are within uncertainty range of the literature-reported values. The predicted metabolite concentrations by the model are also within uncertainty ranges of metabolomic data for 68% of the metabolites. A leave-one-out cross-validation test to evaluate the flux prediction performance of the model showed that metabolic fluxes for the mutants located in the proximity of mutations used for training the model can be predicted more accurately. The constructed model and the parameterization procedure presented in this study pave the way for the construction of larger-scale kinetic models with more narrowly distributed parameter values as new metabolomic/fluxomic data sets are becoming available for E. coli and other organisms.

AB - In contrast to stoichiometric-based models, the development of large-scale kinetic models of metabolism has been hindered by the challenge of identifying kinetic parameter values and kinetic rate laws applicable to a wide range of environmental and/or genetic perturbations. The recently introduced ensemble modeling (EM) procedure provides a promising remedy to address these challenges by decomposing metabolic reactions into elementary reaction steps and incorporating all phenotypic observations, upon perturbation, in its model parameterization scheme. Here, we present a kinetic model of Escherichia coli core metabolism that satisfies the fluxomic data for wild-type and seven mutant strains by making use of the EM concepts. This model encompasses 138 reactions, 93 metabolites and 60 substrate-level regulatory interactions accounting for glycolysis/gluconeogenesis, pentose phosphate pathway, TCA cycle, major pyruvate metabolism, anaplerotic reactions and a number of reactions in other parts of the metabolism. Parameterization is performed using a formal optimization approach that minimizes the discrepancies between model predictions and flux measurements. The predicted fluxes by the model are within the uncertainty range of experimental flux data for 78% of the reactions (with measured fluxes) for both the wild-type and seven mutant strains. The remaining flux predictions are mostly within three standard deviations of reported ranges. Converting the EM-based parameters into a Michaelis-Menten equivalent formalism revealed that 35% of Km and 77% of kcat parameters are within uncertainty range of the literature-reported values. The predicted metabolite concentrations by the model are also within uncertainty ranges of metabolomic data for 68% of the metabolites. A leave-one-out cross-validation test to evaluate the flux prediction performance of the model showed that metabolic fluxes for the mutants located in the proximity of mutations used for training the model can be predicted more accurately. The constructed model and the parameterization procedure presented in this study pave the way for the construction of larger-scale kinetic models with more narrowly distributed parameter values as new metabolomic/fluxomic data sets are becoming available for E. coli and other organisms.

UR - http://www.scopus.com/inward/record.url?scp=84904317199&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84904317199&partnerID=8YFLogxK

U2 - 10.1016/j.ymben.2014.05.014

DO - 10.1016/j.ymben.2014.05.014

M3 - Article

C2 - 24928774

AN - SCOPUS:84904317199

VL - 25

SP - 50

EP - 62

JO - Metabolic Engineering

JF - Metabolic Engineering

SN - 1096-7176

ER -