Abstract

Large-scale kinetic models are used for designing, predicting, and understanding the metabolic responses of living cells. Kinetic models are particularly attractive for the biosynthesis of target molecules in cells as they are typically better than other types of models at capturing the complex cellular biochemistry. Using simpler stoichiometric models as scaffolds, kinetic models are built around a steady-state flux profile and a metabolite concentration vector that are typically determined via optimization. However, as the underlying optimization problem is underdetermined, even after incorporating available experimental omics data, one cannot uniquely determine the operational configuration in terms of metabolic fluxes and metabolite concentrations. As a result, some reactions can operate in either the forward or reverse direction while still agreeing with the observed physiology. Here, we analyze how the underlying uncertainty in intracellular fluxes and concentrations affects predictions of constructed kinetic models and their design in metabolic engineering and systems biology studies. To this end, we integrated the omics data of optimally grown Escherichia coli into a stoichiometric model and constructed populations of non-linear large-scale kinetic models of alternative steady-state solutions consistent with the physiology of the E. coli aerobic metabolism. We performed metabolic control analysis (MCA) on these models, highlighting that MCA-based metabolic engineering decisions are strongly affected by the selected steady state and appear to be more sensitive to concentration values rather than flux values. To incorporate this into future studies, we propose a workflow for moving towards more reliable and robust predictions that are consistent with all alternative steady-state solutions. This workflow can be applied to all kinetic models to improve the consistency and accuracy of their predictions. Additionally, we show that, irrespective of the alternative steady-state solution, increased activity of phosphofructokinase and decreased ATP maintenance requirements would improve cellular growth of optimally grown E. coli.

Highlights

  • Over the last decades, advances in genome editing technologies have allowed the redirection of carbon flow within the organism towards specialty products of interest and desired physiologies (Nielsen, 2017)

  • This work studied the impact of alternative concentration and flux steady states on the conclusions derived from the metabolic control analysis (MCA) outputs of the non-linear kinetic models built around them using the physiology of optimally grown E. coli

  • We show that different flux directionality profiles (FDPs) can lead to distinct metabolic engineering conclusions when analyzing output flux control coefficients (FCCs) of the non-linear models

Read more

Summary

Introduction

Advances in genome editing technologies have allowed the redirection of carbon flow within the organism towards specialty products of interest and desired physiologies (Nielsen, 2017). Available high-throughput sequencing data has enabled the construction of stoichiometric genome-scale metabolic models (GEMs) that describe mathematically the balanced metabolic fluxes within an organism (Thiele and Palsson, 2010) Metabolic models such as these GEMs have been extensively used to characterize overall network behavior of organisms, which can provide guidance about the genes that can be modified to improve a desired product biosynthesis. Improved guidance for metabolic engineering and basic biology will be achieved with kinetic models of the reactions/networks in GEMs. The construction of a kinetic model of metabolism requires knowledge of steady states and/or dynamics of metabolic fluxes and metabolite concentrations that can be used to estimate the unknown kinetic parameters that describe these data.

Methods
Results
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call