From Escherichia coli mutant 13C labeling data to a core kinetic model: A kinetic model parameterization pipeline
Charles J. Foster,Saratram Gopalakrishnan,Maciek R. Antoniewicz,Costas D. Maranas
DOI: https://doi.org/10.1371/journal.pcbi.1007319
2019-09-10
PLoS Computational Biology
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 <sup>13</sup>C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of <em>E</em>. <em>coli</em> containing 74 reactions and 61 metabolites using <sup>13</sup>C-Metabolic Flux Analysis (<sup>13</sup>C-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 <em>E</em>. <em>coli</em> 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 <sup>13</sup>C-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 <sup>13</sup>C-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.
biochemical research methods,mathematical & computational biology