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.

Mathematical Modeling of Microbial Community Dynamics: A Methodological Review

“…Modeling microbial consortia using indirect coupling can also be done by integrating EM analysis with a dynamic framework such as cybernetic models. The framework called hybrid cybernetic model (HCM) [127,128] identifies a relevant subset of EMs as metabolic options for accommodating the metabolic shift in individual species. Geng et al [129] applied HCM to model a situation where the culture medium contains four sugars (i.e., glucose, xylose, mannose, and galactose) which are competitively consumed by three yeast strains (i.e., S. cerevisiae, P. stipitis, and Kluyveromyces marxianus).…”

JuPOETs: A Constrained Multiobjective Optimization Approach to Estimate Biochemical Model Ensembles in the Julia Programming Language

“…We considered a common case in which the same extracellular measurements of A e , B e , C e and cellmass were made on four hypothetical cell types, each having the same biological connectivity but different performance. Network dynamics were modeled using the hybrid cybernetic model with elementary modes (HCM) approach of Ramkrishna and coworkers [34]. In the HCM approach, metabolic networks are first decomposed into a set of elementary modes (EMs) (chemically balanced steady-state pathways, see [35]).…”

Multiscale Modeling of Biopolymer Production in Multicellular Systems

“…In a first step single cell kinetics are described by a state of the art hybrid cybernetic model (HCM) (see Kim et al (2008); ). The HCM allows a systematic derivation of the model equations from elementary mode analysis (see Stelling et al (2002)) and has a moderate complexity, which is useful for nonlinear analysis.…”