A hybrid model of anaerobic E. coli GJT001: combination of elementary flux modes and cybernetic variables

Abstract

Flux balance analysis (FBA) in combination with the decomposition of metabolic networks into elementary modes has provided a route to modeling cellular metabolism. It is dependent, however, on the availability of external fluxes such as substrate uptake or growth rate before estimates can become available of intracellular fluxes. The framework classically does not allow modeling of metabolic regulation or the formulation of dynamic models except through dynamic measurement of external fluxes. The cybernetic modeling approach of Ramkrishna and coworkers provides a dynamic framework for modeling metabolic systems because of its focus on describing regulatory processes based on cybernetic arguments and hence has the capacity to describe both external and internal fluxes. In this article, we explore the alternative of developing hybrid models combining cybernetic models for the external fluxes with the flux balance approach for estimation of the internal fluxes. The approach has the merit of the simplicity of the early cybernetic models and hence computationally facile while also providing detailed information on intracellular fluxes. The hybrid model of this article is based on elementary mode decomposition of the metabolic network. The uptake rates for the various elementary modes are combined using global cybernetic variables based on maximizing substrate uptake rates. Estimation of intracellular metabolism is based on its stoichiometric coupling with the external fluxes under the assumption of (pseudo-) steady state conditions. The set of parameters of the hybrid model was estimated with the aid of nonlinear optimization routine, by fitting simulations with dynamic experimental data on concentrations of biomass, substrate, and fermentation products. The hybrid model estimations were tested with FBA (based on measured substrate uptake rate) for two different metabolic networks (one is a reduced network which fixes ATP contribution to the biomass and maintenance requirement of ATP, and the other network is a more complex network which has a separate reaction for maintenance.) for the same experiment involving anaerobic growth of E. coli GJT001. The hybrid model estimated glucose consumption and all fermentation byproducts to better than 10%. The FBA makes similar estimations of fermentation products, however, with the exception of succinate. The simulation results show that the global cybernetic variables alone can regulate the metabolic reactions obtaining a very satisfactory fit to the measured fermentation byproducts. In view of the hybrid model's ability to predict biomass growth and fermentation byproducts of anaerobic E. coli GJT001, this reduced order model offers a computationally efficient alternative to more detailed models of metabolism and hence useful for the simulation of bioreactors.