Abstract:Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13C-Metabolic Flux Analysis (13C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated.Microbial production hosts are used for production of a wide range of commodity chemicals. Improving the conversion efficiency of microbial strains is critical to the economic viability and the continued push towards the use of environmentally neutral bioprocesses as a means for producing the chemicals society depends on. Metabolic models have played a key role in helping us to predict metabolic behavior in response to environmental and genetic perturbation that can maximize efficiency. Recently, kinetic models of metabolism have re-emerged as a means for characterizing metabolism, offering improvements over their stoichiometric counterparts in both the type of information that can be gleaned from them, and in prediction accuracy. Despite recent developments, a lack of raw experimental data needed for flux elucidation and, subsequently, kinetic parameterization, and high computation cost have prevented the development of a uniform workflow for construction of the most informative kinetic models. Here, we have incorporated raw 13C-isotopic labeling data and a computationally inexpensive parameterization algorithm into a kinetic parameterization pipeline to ensure that the resulting kinetic model (k-ecoli74) conforms to experimental data. We show how the use of an identical metabolic network for flux elucidation and kinetic parameterization influences predictive capabilities.

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

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

Optimal re-design of primary metabolism in Escherichia coli using linlog kinetics

From Escherichia coli mutant 13C labeling data to a core kinetic model: A kinetic model parameterization pipeline

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Complexity analysis and parameter estimation of dynamic metabolic systems

Elucidating yeast glycolytic dynamics at steady state growth and glucose pulses through kinetic metabolic modeling

Towards Stable Kinetics of Large Metabolic Networks: Nonequilibrium Potential Function Approach

Application of dynamic metabolic flux analysis for process modeling: Robust flux estimation with regularization, confidence bounds, and selection of elementary modes

Building kinetic models for metabolic engineering

Dynamic optimization of metabolic networks coupled with gene expression

Towards kinetic modeling of metabolic networks with incomplete parameters

Structural Thermokinetic Modelling

Omics data for sampling thermodynamically feasible kinetic models

Kriging-Based Parameter Estimation Algorithm for Metabolic Networks Combined with Single-Dimensional Optimization and Dynamic Coordinate Perturbation

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

Parameter Estimation and Identifiability in Kinetic Flux Profiling Models of Metabolism

Global Optimization Approach for Parameter Estimation in Stochastic Dynamic Models of Biosystems

Enhancing the accuracy of genome-scale metabolic models with kinetic information

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

Dynamic modeling of the central carbon metabolism of Escherichia coli