Abstract:Abstract Background It is quite important to simulate the metabolic changes of a cell in response to the change in culture environment and/or specific gene knockouts particularly for the purpose of application in industry. If this could be done, the cell design can be made without conducting exhaustive experiments, and one can screen out the promising candidates, proceeded by experimental verification of a select few of particular interest. Although several models have so far been proposed, most of them focus on the specific metabolic pathways. It is preferred to model the whole of the main metabolic pathways in Escherichia coli , allowing for the estimation of energy generation and cell synthesis, based on intracellular fluxes and that may be used to characterize phenotypic growth. Results In the present study, we considered the simulation of the main metabolic pathways such as glycolysis, TCA cycle, pentose phosphate (PP) pathway, and the anapleorotic pathways using enzymatic reaction models of E. coli . Once intracellular fluxes were computed by this model, the specific ATP production rate, the specific CO 2 production rate, and the specific NADPH production rate could be estimated. The specific ATP production rate thus computed was used for the estimation of the specific growth rate. The CO 2 production rate could be used to estimate cell yield, and the specific NADPH production rate could be used to determine the flux of the oxidative PP pathway. The batch and continuous cultivations were simulated where the changing patterns of extracellular and intra-cellular metabolite concentrations were compared with experimental data. Moreover, the effects of the knockout of such pathways as Ppc, Pck and Pyk on the metabolism were simulated. It was shown to be difficult for the cell to grow in Ppc mutant due to low concentration of OAA, while Pck mutant does not necessarily show this phenomenon. The slower growth rate of the Ppc mutant was properly estimated by taking into account the lower specific ATP production rate. In the case of Pyk mutant, the enzyme level regulation was made clear such that Pyk knockout caused PEP concentration to be up-regulated and activated Ppc, which caused the increase in MAL concentration and backed up reduced PYR through Mez, resulting in the phenotypic growth characteristics similar to the wild type. Conclusions It was shown to be useful to simulate the main metabolism of E. coli for understanding metabolic changes inside the cell in response to specific pathway gene knockouts, considering the whole main metabolic pathways. The comparison of the simulation result with the experimental data indicates that the present model could simulate the effect of the specific gene knockouts to the changes in the metabolisms to some extent.

Dynamic simulation of an in vitro multi-enzyme system

Kinetic profiling of metabolic specialists demonstrates stability and consistency of in vivo enzyme turnover numbers

Analysis of l-glutamic acid fermentation by using a dynamic metabolic simulation model of Escherichia coli

Dynamic modeling of the central carbon metabolism of Escherichia coli

Modeling and simulation of the main metabolism in Escherichia coli and its several single-gene knockout mutants with experimental verification

Interrogating the effect of enzyme kinetics on metabolism using differentiable constraint-based models

Population-based Optimization for Kinetic Parameter Identification in Glycolytic Pathway in Saccharomyces cerevisiae

Combining systems and synthetic biology for in vivo enzymology

Towards Kinetic Modeling of Global Metabolic Networks： Methylobacterium Extorquens AM1 Growth As Validation

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

GENERIC ENZYMATIC RATE EQUATION UNDER LIVING CONDITIONS

A method for estimation of elasticities in metabolic networks using steady state and dynamic metabolomics data and linlog kinetics

Modeling, Simulating, and Parameter Fitting of Biochemical Kinetic Experiments

Kinetic modelling of the cellular metabolic responses underpinning in vitro glycolysis assays

Beyond the Michaelis-Menten equation: Accurate and efficient estimation of enzyme kinetic parameters

A Generic Rate Equation for modeling Enzymatic Reactions under Living Conditions

Omics data for sampling thermodynamically feasible kinetic models

Multi-Scale Kinetics Modeling for Cell Culture Process with Metabolic State Transition

Quantitative Online Monitoring of an Immobilized Enzymatic Network by Ion Mobility-Mass Spectrometry

Towards Kinetic Modeling of Global Metabolic Networks with Incomplete Experimental Input on Kinetic Parameters

Development of a Digital Twin for Enzymatic Hydrolysis Processes