Simultaneous parameters identifiability and estimation of an E. coli metabolic network model

In silico Determination of Some Conditions Leading to Glycolytic Oscillations and Their Interference With Some Other Processes in E. coli Cells

Autonomous oscillations of species levels in the glycolysis express the self-control of this essential cellular pathway belonging to the central carbon metabolism (CCM), and this phenomenon takes place in a large number of bacteria. Oscillations of glycolytic intermediates in living cells occur according to the environmental conditions and to the cell characteristics, especially the adenosine triphosphate (ATP) recovery system. Determining the conditions that lead to the occurrence and maintenance of the glycolytic oscillations can present immediate practical applications. Such a model-based analysis allows in silico (model-based) design of genetically modified microorganisms (GMO) with certain characteristics of interest for the biosynthesis industry, medicine, etc. Based on our kinetic model validated in previous works, this paper aims to in silico identify operating parameters and cell factors leading to the occurrence of stable glycolytic oscillations in the Escherichia coli cells. As long as most of the glycolytic intermediates are involved in various cellular metabolic pathways belonging to the CCM, evaluation of the dynamics and average level of its intermediates is of high importance for further applicative analyses. As an example, by using a lumped kinetic model for tryptophan (TRP) synthesis from literature, and its own kinetic model for the oscillatory glycolysis, this paper highlights the influence of glycolytic oscillations on the oscillatory TRP synthesis through the PEP (phosphoenolpyruvate) glycolytic node shared by the two oscillatory processes. The numerical analysis allows further TRP production maximization in a fed-batch bioreactor (FBR).

Structured model and parameter estimation in plant cell cultures of Thevetia peruviana

In this work, a mechanistic model for predicting the dynamic behavior of extracellular and intracellular nutrients, biomass production, and the main metabolites involved in the central carbon metabolism in plant cell cultures of Thevetia peruviana is presented. The proposed model is the first mechanistic model implemented for plant cell cultures of this species, and includes 28 metabolites, 33 metabolic reactions, and 61 parameters. Given the over-parametrization of the model, its nonlinear nature and the strong correlation among the effects of the parameters, a parameter estimation routine based on identifiability analysis was implemented. This routine reduces the parameter's search space by selecting the most sensitive and linearly independent parameters. Results have shown that only 19 parameters are identifiable. Finally, the model was used for analyzing the fluxes distribution in plant cell cultures of T. peruviana. This analysis shows high uptake of phosphates and parallel uptake of glucose and fructose. Furthermore, it has pointed out the main central carbon metabolism routes for promoting biomass production in this cell culture.

Development of an accurate kinetic model for the central carbon metabolism of Escherichia coli

Background

A kinetic model provides insights into the dynamic response of biological systems and predicts how their complex metabolic and gene regulatory networks generate particular functions. Of many biological systems, Escherichia coli metabolic pathways have been modeled extensively at the enzymatic and genetic levels, but existing models cannot accurately reproduce experimental behaviors in a batch culture, due to the inadequate estimation of a specific cell growth rate and a large number of unmeasured parameters.

Results

In this study, we developed a detailed kinetic model for the central carbon metabolism of E. coli in a batch culture, which includes the glycolytic pathway, tricarboxylic acid cycle, pentose phosphate pathway, Entner-Doudoroff pathway, anaplerotic pathway, glyoxylate shunt, oxidative phosphorylation, phosphotransferase system (Pts), non-Pts and metabolic gene regulations by four protein transcription factors: cAMP receptor, catabolite repressor/activator, pyruvate dehydrogenase complex repressor and isocitrate lyase regulator. The kinetic parameters were estimated by a constrained optimization method on a supercomputer. The model estimated a specific growth rate based on reaction kinetics and accurately reproduced the dynamics of wild-type E. coli and multiple genetic mutants in a batch culture.

Conclusions

This model overcame the intrinsic limitations of existing kinetic models in a batch culture, predicted the effects of multilayer regulations (allosteric effectors and gene expression) on central carbon metabolism and proposed rationally designed fast-growing cells based on understandings of molecular processes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0511-x) contains supplementary material, which is available to authorized users.