A Model Metabolic Cycle Simulated with the Mathematica Program


A system of differential equations has been used to simulate a model metabolic cycle, containing only one initial substrate and 8 irreversible enzymes. The metabolic course of the intermediates (products/substrates) was also explored in different situations: changes in the kinetic constants of one of the enzymes of the cycle; introduction of one reversible step; consecutive increments in the Km values of each enzyme of the cycle and the presence of two initial substrates. A decrease or an increase in the Km value of one of the enzymes of the cycle promotes a decrease or an increase in the steady state level of its own substrate, respectively; by the contrary, a decrease or an increase in the Vmax value promotes, respectively, an increase or a decrease in the stationary level of the corresponding substrate. A comparison between a linear and a cyclic metabolic pathway is also presented.

Share and Cite:

García-Herrero, V. , López-Cánovas, F. and Sillero, A. (2014) A Model Metabolic Cycle Simulated with the Mathematica Program. Journal of Biomedical Science and Engineering, 7, 286-295. doi: 10.4236/jbise.2014.75031.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] López-Cánovas, F.J., Gomes, P.J. and Sillero, A. (2013) Mathematica Program: Its Use to Simulate Metabolic Irreversible Pathways and Inhibition of the First Enzyme of a Metabolic Pathway as Visuaized with the Reservoir Model. Computers in Biology and Medicine, 42, 853-864.
[2] Voet, D. and Voet, J.G. (2004) Biochemistry, XV +1591, 482-493.
[3] Cleland, W.W. (1970) Steady State Kinetics in the Enzymes. 3rd Edition, Academic Press, New York, 1-65.
[4] Krauss, M., Schaller, S., Borchers, S., Findelsen, R., Lippert, J. and Kuepfer, L. (2012) Integrating Cellular Metabolism into a Multiscale Whole-Body Model. PLOS Computational Biology, 8, e1002750.
[5] Fell, D. (1997) Understanding the Control of Metabolism. Portland Press; Distributed by Ashgate Pub. Co. in North America, London; Miami Brookfild, VT.
[6] Segel, I.H. (1975) Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady State Enzyme Systems. Wiley, New York.
[7] Klipp, R., Herwig, R., Kowald, A., Wierling, C. and Lehrach, H. (2005) System Biology in Practice. Wiley-VCH Verlag GmbH &Co, Winheim (FRG).
[8] Alberty, R.A. (2011) Enzyme Kinetics. Rapid-Equilibrium Applications of Mathematica, Hoboken.
[9] Alon, U. (2006) An Introduction to System Biology. Design Principles of Biological Circuits. CRC Pres, Taylor & Francis Group, London.
[10] Wolfram, S. (2013) Wolfram Mathematica Tutorial Collection.
[11] Bisswanger, H. (2004) Enzyme Kinetics. Principle and Methods, Wiley-VCH, Weinheim (FRG).
[12] Albe, K.R. and Wright, B.E. (1992) Systems Analysis of the Tricarboxylic Acid Cycle in Dictyostelium Discoideum. II. Control Analysis. Journal of Biological Chemistry, 267, 3106-3114.
[13] Oliveira, J.S., Bailey, C.G., Jones-Oliveira, J.B., Dixon, D.A., Gull, D.W. and Chandler, M.L. (2003) A Computational Model for the Identification of Biochemical Pathways in the Krebs Cycle. Journal of Computational Biology, 10, 57-82.
[14] Wu, F., Yang, F., Vinnakota, K.C. and Beard, D.A. (2007) Computer Modeling of Mitochondrial Tricarboxylic Acid Cycle, Oxidative Phosphorylation, Metabolite Transport, and Electrophysiology. Journal of Biological Chemistry, 282, 24525-24537.
[15] Nazaret, C., Heiske, M., Thurley, K. and Mazat, J.P. (2009) Mitochondrial Energetic Metabolism: A Simplified Model of TCA Cycle with ATP Production. Journal of Theoretical Biology, 258, 455-464.
[16] Amador-Noguez, D., Feng, X.J., Fan, J., Roquet, N., Rabitz, H. and Rabinowitz, J.D. (2010) Systems-Level Metabolic Flux Profiling Elucidates a Complete, Bifurcated Tricarboxylic Acid Cycle in Clostridium Acetobutylicum. Journal of Bacteriology, 192, 4452-4461.
[17] Smith, A.C. and Robinson, A.J. (2011) A Metabolic Model of the Mitochondrion and Its Use in Modelling Diseases of the Tricarboxylic Acid Cycle. BMC Systems Biology, 5, 102.
[18] Bachmann, C. and Colombo, J.P. (1981) Computer Simulation of the Urea Cycle: Trials for an Appropriate Model. Enzyme, 26, 259-264.
[19] Morris Jr., S.M. (2002) Regulation of Enzymes of the Urea Cycle and Arginine Metabolism. Annual Review of Nutrition, 22, 87-105.
[20] Maher, A.D., Kuchel, P.W., Ortega, F., de Atauri, P., Centelles, J. and Cascante, M. (2003) Mathematical Modelling of the Urea Cycle. A Numerical Investigation into Substrate Channelling. European Journal of Biochemistry, 270, 3953-3961.
[21] Sillero, A., Selivanov, V.A. and Cascante, M. (2006) Pentose Phosphate and Calvin Cycles: Similarities and Three-Dimensional Views. Biochemistry and Molecular Biology Education, 34, 275-277.
[22] Bogorad, I.W., Lin, T.-S. and Liao, J.C. (2013) Synthetic Non-Oxidative Glycolysis Enables Complete Carbon Conservation. Journal of Biological Chemistry, 202, 693-697.

Copyright © 2022 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.