The nonlinear behavior of metabolic systems can arise from at least two different sources. One comes from the nonlinear kinetics of chemical reactions in metabolism and the other from nonlinearity associated with regulatory processes. Consequently, organisms at a constant growth rate (as experienced in a chemostat) could display multiple metabolic states or display complex oscillatory behavior both with potentially serious implications to process operation. This paper explores the nonlinear behavior of a metabolic model of Escherichia coli growth on mixed substrates with sufficient detail to include regulatory features through the cybernetic postulate that metabolic regulation is the consequence of a dynamic objective function ensuring the organism’s survival. The chief source of nonlinearity arises from the optimal formulation with the metabolic state determined by a convex combination of reactions contributing to the objective function. The model for anaerobic growth of E. coli was previously examined for multiple steady states in a chemostat fed by a mixture of glucose and pyruvate substrates under very specific conditions and experimentally verified. In this article, we explore the foregoing model for nonlinear behavior over the full range of parameters, γ (the fractional concentration of glucose in the feed mixture) and D (the dilution rate). The observed multiplicity is in the cybernetic variables combining elementary modes. The results show steady-state multiplicity up to seven. No Hopf bifurcation was encountered, however. Bifurcation analysis of cybernetic models is complicated by the non-differentiability of the cybernetic variables for enzyme activities. A methodology is adopted here to overcome this problem, which is applicable to more complicated metabolic networks.
References
[1]
Alzate, C.A.C.; Astudillo, I.C.P. Importance of stability study of continuous systems for ethanol production. J. Biotechnol. 2011, 151, 43–55.
[2]
Elnashaie, S.S.; Grace, J.R. Complexity, bifurcation and chaos in natural and man-made lumped and distributed systems. Chem. Eng. Sci. 2007, 62, 3295–3325, doi:10.1016/j.ces.2007.02.047.
[3]
Trinh, C.T.; Wlaschin, A.; Srienc, F. Elementary mode analysis: A useful metabolic pathway analysis tool for characterizing cellular metabolism. Appl. Microbiol. Biotechnol. 2009, 81, 813–826.
[4]
Wagner, C.; Urbanczik, R. The geometry of the flux cone of a metabolic network. Biophys. J. 2005, 89, 3837–3845, doi:10.1529/biophysj.104.055129.
[5]
Orth, J.D.; Thiele, I.; Palsson, B.O. What is flux balance analysis? Nat. Biotechnol. 2010, 28, 245–248.
[6]
Reed, J.L.; Palsson, B.O. Thirteen years of building constraint-based in silico models of Escherichia coli. J. Bacteriol. 2003, 185, 2692–2699.
[7]
Wiechert, W.; Noack, S. Mechanistic pathway modeling for industrial biotechnology: Challenging but worthwhile. Curr. Opin. Biotechnol. 2011, 22, 604–610, doi:10.1016/j.copbio.2011.01.001.
[8]
Ramkrishna, D.; Song, H.S. Dynamic models of metabolism: Review of the cybernetic approach. AIChE J. 2012, 58, 986–997, doi:10.1002/aic.13734.
[9]
Song, H.S.; DeVilbiss, F.; Ramkrishna, D. Modeling metabolic systems: The need for dynamics. Curr. Opin. Chem. Eng. 2013. in press.
[10]
Namjoshi, A.A.; Ramkrishna, D. Multiplicity and stability of steady states in continuous bioreactors: Dissection of cybernetic models. Chem. Eng. Sci. 2001, 56, 5593–5607, doi:10.1016/S0009-2509(01)00166-X.
[11]
Namjoshi, A.A.; Hu, W.S.; Ramkrishna, D. Unveiling steady-state multiplicity in hybridoma cultures: The cybernetic approach. Biotechnol. Bioeng. 2003, 81, 80–91, doi:10.1002/bit.10447.
[12]
Simpson, D.J.; Kompala, D.S.; Meiss, J.D. Discontinuity induced bifurcations in a model of Saccharomyces cerevisiae. Math. Biosci. 2009, 218, 40–49, doi:10.1016/j.mbs.2008.12.005.
[13]
Kim, J.I.; Varner, J.D.; Ramkrishna, D. A hybrid model of anaerobic E. coli GJT001: Combination of elementary flux modes and cybernetic variables. Biotechnol. Prog. 2008, 24, 993–1006, doi:10.1002/btpr.73.
[14]
Song, H.S.; Morgan, J.A.; Ramkrishna, D. Systematic development of hybrid cybernetic models: Application to recombinant yeast co-consuming glucose and xylose. Biotechnol. Bioeng. 2009, 103, 984–1002, doi:10.1002/bit.22332.
[15]
Song, H.S.; Ramkrishna, D. Reduction of a set of elementary modes using yield analysis. Biotechnol. Bioeng. 2009, 102, 554–568, doi:10.1002/bit.22062.
[16]
Kim, J.I.; Song, H.S.; Sunkara, S.R.; Lali, A.; Ramkrishna, D. Exacting predictions by cybernetic model confirmed experimentally: Steady state multiplicity in the chemostat. Biotechnol. Prog. 2012, 28, 1160–1166, doi:10.1002/btpr.1583.
[17]
Seydel, R. Practical Bifurcation and Stability Analysis; Springer Verlag: New York, NY, USA, 2009; Volume5.
[18]
Dhooge, A.; Govaerts, W.; Kuznetsov, Y.A. MATCONT: A MATLAB package for numerical bifurcation analysis of ODEs. ACM Trans. Math. Softw. (TOMS) 2003, 29, 141–164, doi:10.1145/779359.779362.
[19]
Dhooge, A.; Govaerts, W.; Kuznetsov, Y.A.; Meijer, H.; Sautois, B. New features of the software MatCont for bifurcation analysis of dynamical systems. Math. Comput. Model. Dyn. Syst. 2008, 14, 147–175, doi:10.1080/13873950701742754.
[20]
Schuster, S.; Fell, D.A.; Dandekar, T. A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks. Nat. Biotechnol. 2000, 18, 326–332, doi:10.1038/73786.
[21]
Kompala, D.S.; Ramkrishna, D.; Jansen, N.B.; Tsao, G.T. Investigation of bacterial-growth on mixed substrates—Experimental evaluation of cybernetic models. Biotechnol. Bioeng. 1986, 28, 1044–1055.
[22]
Young, J.D.; Ramkrishna, D. On the matching and proportional laws of cybernetic models. Biotechnol. Prog. 2007, 23, 83–99.
[23]
Von Kamp, A.; Schuster, S. Metatool 5.0: Fast and flexible elementary modes analysis. Bioinformatics 2006, 22, 1930–1931, doi:10.1093/bioinformatics/btl267.
[24]
Kim, J.I. A Hybrid Cybernetic Modeling for the Growth of Escherichia coli in Glucose-Pyruvate Mixtures. Ph.D. Thesis, Purdue University, West Lafayette, IN, USA, May 2008.
[25]
Simpson, D.J.W. Bifurcations in Piecewise-Smooth Continuous Systems; World Scientific Publishing Co.: Hackensack, NJ, USA, 2010.
