Dynamics of Muscle Glycogenolysis Modeled with pH Time-Course Computation and pH Dependent Reaction Equilibria and Enzyme Kinetics
Cellular metabolites are moieties defined by their specific binding constants to H1, Mg21, and K1 or anions without ligands. As a consequence, every biochemical reaction in the cytoplasm has an associated proton stoichiometry that is generally noninteger- and pH-dependent. Therefore, with metabolic flux, pH is altered in a medium with finite buffer capacity. Apparent equilibrium constants and maximum enzyme velocities, which are functions of pH, are also altered. We augmented an earlier mathematical model of skeletal muscle glycogenolysis with pH-dependent enzyme kinetics and reaction equilibria to compute the time course of pH changes. Analysis shows that kinetics and final equilibrium states of the closed system are highly constrained by the pH-dependent parameters. This kinetic model of glycogenolysis, coupled to creatine kinase and adenylate kinase, simulated published experiments made with a cell-free enzyme mixture to reconstitute the network and to synthesize PCr and lactate in vitro. Using the enzyme kinetic and thermodynamic data in the literature, the simulations required minimal adjustments of parameters to describe the data. These results show that incorporation of appropriate physical chemistry of the reactions with accurate kinetic modeling gives a reasonable simulation of experimental data and is necessary for a physically correct representation of the metabolic network. The approach is general for modeling metabolic networks beyond the specific pathway and conditions presented here. Model directly based on: Vinnakota K, Kemp ML, Kushmerick MJ, Dynamics of Muscle Glycogenolysis Modeled with pH Time-Course Computation and pH Dependent Reaction Equilibria and Enzyme Kinetics, Biophy J (91), Aug 2006, p1264-1287
The equations for this model may be viewed by running the JSim model applet and clicking on the Source tab at the bottom left of JSim's Run Time graphical user interface. The equations are written in JSim's Mathematical Modeling Language (MML). See the Introduction to MML and the MML Reference Manual. Additional documentation for MML can be found by using the search option at the Physiome home page.
- Download JSim model MML code (text):
- Download translated SBML version of model (if available):
- Reproduce Figures 1 and 2 of 2006 paper:
- Reproduce Figures 12 and 13 of 2006 paper:
- Reproduce Figure 15 of 2006 paper:
Download Matlab files of original mode (ZIP file)
We welcome comments and feedback for this model. Please use the button below to send comments:
Rutman, R. J., and P. George. 1961. Hydrogen ion effects in highenergy phosphate reactions. Proc. Natl. Acad. Sci. USA. 47:1094-1109. George, P., and R. J. Rutan. 1960. The high energy phosphate bond. Concept. Prog. Biophys. Biophys. Chem. 10:2-53. George, P., R. C. Phillips, and R. J. Rutman. 1963. Estimates of thermodynamic data for the formation of Mg21 complexes of ATP and ADP at zero ionic strength. Biochemistry. 2:508-512. Smith, R. M., and R. A. Alberty. 1956. The apparent stability constants of ionic complexes of various adenosine phosphates with divalent cations. J. Am. Chem. Soc. 78:2376-2384. Alberty, R. A. 2004. A short history of the thermodynamics of enzymecatalyzed reactions. J. Biol. Chem. 279:27831-27836. Alberty, R. A. 2003. Thermodynamics of Biochemical Reactions. John Wiley and Sons, Hoboken, NJ. Harkema, S. J., and R. A. Meyer. 1997. Effect of acidosis on control of respiration in skeletal muscle. Am. J. Physiol. Cell Physiol. 41: C491-C500. Kushmerick, M. J. 1997. Multiple equilibria of cations with metabolites in muscle bioenergetics. Am. J. Physiol. 272:C1739-C1747. Sato, K., Y. Kashiwaya, C. Keon, N. Tsuchiya, M. T. King, G. K. Radda, B. Chance, K. Clarke, and R. L. Veech. 1995. Insulin, ketone bodies, and mitochondrial energy transduction. FASEB J. 9:651-658. Beard, D. A., E. Babson, E. Curtis, and H. Qian. 2004. Thermodynamic constraints for biochemical networks. J. Theor. Biol. 228:327-333. Mulquiney, P. J., W. A. Bubb, and P. W. Kuchel. 1999. Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: in vivo kinetic characterization of 2,3-bisphosphoglycerate synthase/phosphatase using 13C and 31P NMR. Biochem. J. 342:567-580. Mulquiney, P. J., and P. W. Kuchel. 1999. Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: computer simulation and metabolic control analysis. Biochem. J. 342:597-604. Mulquiney, P. J., and P. W. Kuchel. 1999. Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: equations and parameter refinement. Biochem. J. 342:581-596. Scopes, R. K. 1973. Studies with a reconstituted muscle glycolytic system. The rate and extent of creatine phosphorylation by anaerobic glycolysis. Biochem. J. 134:197-208. Scopes, R. K. 1974. Studies with a reconstituted muscle glycolytic system. The rate and extent of glycolysis in simulated post-mortem conditions. Biochem. J. 142:79-86. Lambeth, M. J., and M. J. Kushmerick. 2002. A computational model for glycogenolysis in skeletal muscle. Ann. Biomed. Eng. 30:808-827. Robergs, R. A., F. Ghiasvand, and D. Parker. 2004. Biochemistry of exercise-induced metabolic acidosis. Am. J. Physiol. Regul. Integr. Comp. Physiol. 287:R502R516. Needham, D. M. 1971. Machina Carnis: the Biochemistry of Muscular Contraction in its Historical Development. Cambridge University Press, London, UK. Gevers, W. 1977. Generation of protons by metabolic processes in heart cells. J. Mol. Cell. Cardiol. 9:867-874. Portner, H. O., N. Heisler, and M. K. Grieshaber. 1984. Anaerobiosis and acid-base status in marine invertebrates: a theoretical analysis of proton generation by anaerobic metabolism. J. Comp. Physiol. [B]. 155:1-12. Hochachka, P. W., and T. P. Mommsen. 1983. Protons and anaerobiosis. Science. 219:1391-1397. Achs, M. J., and D. Garfinkel. 1977. Computer simulation of rat heart metabolism after adding glucose to the perfusate. Am. J. Physiol. 232: R175-R184. Achs, M. J., and D. Garfinkel. 1977. Computer simulation of energy metabolism in anoxic perfused rat heart. Am. J. Physiol. 232:R164-R174. Achs, M. J., D. Garfinkel, and L. H. Opie. 1982. Computer simulation of metabolism of glucose-perfused rat heart in a work-jump. Am. J. Physiol. 243:R389-R399. Robinson, R. A., and R. H. Stokes. 1968. Electrolyte Solutions. Butterworths, London, UK. Dawson, R. M. C., D. C. Elliott, W. H. Elliott, and K. M. Jones. 1969. Data for Biochemical Research. Oxford University Press, London, UK. Clarke, E. C. W., and D. H. Glew. 1980. Evaluation of Debye-Huckel limiting slopes for water between 0 and 150C. J. Chem. Soc., Faraday Trans. 1. 76:1911-1916. Adams, G. R., J. M. Foley, and R. A. Meyer. 1990. Muscle buffer capacity estimated from pH changes during rest-to-work transitions. J. Appl. Physiol. 69:968-972. Martell, A. E., and R. M. Smith. 2003. NIST Standard Reference Database 46 Version 7.0: NIST Critically Selected Stability Constants of Metal Complexes. NIST Standard Reference Data, Gaithersburg, MD. Teague, W. E., Jr., E. M. Golding, and G. P. Dobson. 1996. Adjustment of K9 for the creatine kinase, adenylate kinase and ATP hydrolysis equilibria to varying temperature and ionic strength. J. Exp. Biol. 199:509-512. Cornish-Bowden, A. 1995. Fundamentals of Enzyme Kinetics. Portland Press, London, UK. Segel, I. H. 1975. Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady State Enzyme Systems. Wiley, New York. Dyson, J. E., and E. A. Noltmann. 1968. The effect of pH and temperature on the kinetic parameters of phosphoglucose isomerase. Participation of histidine and lysine in a proposed dual function mechanism. J. Biol. Chem. 243:1401-1414. Keleti, T., and J. Batke. 1967. The kinetics of the reactions catalyzed by d-glyceraldehyde-3-phosphate dehydrogenase. 3. The pH dependence of apparent Michaelis constants and maximum velocity. Enzymologia. 33:65-79. Najjar, V. A. 1948. The isolation and properties of phosphoglucomutase. J. Biol. Chem. 175:281-290. Paetkau, V., and H. A. Lardy. 1967. Phosphofructokinase. Correlation of physical and enzymatic properties. J. Biol. Chem. 242:2035-2042. Simon, L. M., M. Kotorman, B. Szajani, and L. Boross. 1985. Comparative studies on soluble and immobilized rabbit muscle pyruvate kinase. Appl. Biochem. Biotechnol. 11:195-205. Stankiewicz, P. J., and L. F. Hass. 1986. The catalytic bimodality of mammalian phosphoglycerate mutase. J. Biol. Chem. 261:12715-12721. Winer, A. D., and G. W. Schwert. 1958. Lactic dehydrogenase. IV. The influence of pH on the kinetics of the reaction. J. Biol. Chem. 231: 1065-1083. Yunis, A. A., E. H. Fischer, and E. G. Krebs. 1960. Crystallization and properties of human muscle phosphorylases a and b. J. Biol. Chem. 235:3163-3168. Abraham, M., L. Horvath, and B. Szajani. 1985. Isolation and characterization of pig muscle aldolase. A comparative study. Comp. Biochem. Physiol. B. 80:847-852. Bakker, B. M., P. A. Michels, F. R. Opperdoes, and H. V. Westerhoff. 1997. Glycolysis in bloodstream from Trypanosoma brucei can be understood in terms of the kinetics of the glycolytic enzymes. J. Biol. Chem. 272:3207-3215. Teusink, B., J. Passarge, C. A. Reijenga, E. Esgalhado, C. C. van der Weijden, M. Schepper, M. C. Walsh, B. M. Bakker, K. van Dam, H. V. Westerhoff, and J. L. Snoep. 2000. Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry. Eur. J. Biochem. 267:5313-5329. Ostro, M. J., and T. P. Fondy. 1977. Isolation and characterization of multiple molecular forms of cytosolic NAD-linked glycerol-3-phosphate dehydrogenase from normal and neoplastic rabbit tissues. J. Biol. Chem. 252:5575-5583. Hairer, E., and G. Wanner. 2004. Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. R. Bank, R. L. Graham, J. Stoer, R. Varga, and H. Yserentant, editors. Springer, New York. Meyer, R. A., and J. M. Foley. 1996. Cellular processes integrating the metabolic response to exercise. In Handbook of Physiology. L. B. Rowell and J. T. Shepherd, editors. Oxford University Press, New York. 841-869. MacDonald, J. A., and K. B. Storey. 2002. Purification and characterization of fructose bisphosphate aldolase from the ground squirrel, Spermophilus lateralis: enzyme role in mammalian hibernation. Arch. Biochem. Biophys. 408:279-285. Vinnakota, K. C., M. L. Kemp, and M. J. Kushmerick. 2005. Dynamics of pH changes and phoshoenergetics in mouse skeletal muscle during transient anoxia are modeled by physical chemical principles. Late Breaking Abstracts, XXXV International Congress of Physiological Sciences/Experimental Biology, San Diego. LB158. Cornell, N. W., M. Leadbetter, and R. L. Veech. 1979. Effects of free magnesium concentration and ionic strength on equilibrium constants for the glyceraldehyde phosphate dehydrogenase and phosphoglycerate kinase reactions. J. Biol. Chem. 254:6522-6527. Cornish-Bowden, A., and M. L. Cardenas. 2000. Irreversible reactions in metabolic simulations: how reversible is irreversible? In BTK2000: Animating the Cellular Map. J. H. S. Hofmeyr, J. H. Rohwer, and J. L. Snoep, editors. Stellenbosch University Press. Stellenbosch, South Africa. 65-71. Goldberg, R. N., and Y. B. Tewari. 1994. Thermodynamics of enzymecatalyzed reactions: part 2. Transferases. J. Phys. Chem. Ref. Data. 23: 547-617. Goldberg, R. N., Y. B. Tewari, D. Bell, K. Fazio, and E. Anderson. 1993. Thermodynamics of enzyme-catalyzed reactions: part 1. Oxidoreductases. J. Phys. Chem. Ref. Data. 22:515-582. Robitaille, P.-M. L., and P. A. Robitaille. 1969. J. Gordon Brown, G., and G. G. Brown. 1991. An analysis of the pH-dependent chemicalshift behavior of phosphorus-containing metabolites. J. Magn. Reson. 92:73-84. Beard, D. A., and H. Qian. 2005. Thermodynamic-based computational profiling of cellular regulatory control in hepatocyte metabolism. Am. J. Physiol. Endocrinol. Metab. 288:E633-E644.
Please cite https://www.imagwiki.nibib.nih.gov/physiome in any publication for which this software is used and send one reprint to the address given below:
The National Simulation Resource, Director J. B. Bassingthwaighte, Department of Bioengineering, University of Washington, Seattle WA 98195-5061.
Model development and archiving support at https://www.imagwiki.nibib.nih.gov/physiome provided by the following grants: NIH U01HL122199 Analyzing the Cardiac Power Grid, 09/15/2015 - 05/31/2020, NIH/NIBIB BE08407 Software Integration, JSim and SBW 6/1/09-5/31/13; NIH/NHLBI T15 HL88516-01 Modeling for Heart, Lung and Blood: From Cell to Organ, 4/1/07-3/31/11; NSF BES-0506477 Adaptive Multi-Scale Model Simulation, 8/15/05-7/31/08; NIH/NHLBI R01 HL073598 Core 3: 3D Imaging and Computer Modeling of the Respiratory Tract, 9/1/04-8/31/09; as well as prior support from NIH/NCRR P41 RR01243 Simulation Resource in Circulatory Mass Transport and Exchange, 12/1/1980-11/30/01 and NIH/NIBIB R01 EB001973 JSim: A Simulation Analysis Platform, 3/1/02-2/28/07.