Roles of the creatine kinase system and myoglobin in maintaining energetic state in the working heart
© Wu and Beard. 2009
Received: 23 July 2008
Accepted: 19 February 2009
Published: 19 February 2009
Skip to main content
© Wu and Beard. 2009
Received: 23 July 2008
Accepted: 19 February 2009
Published: 19 February 2009
The heart is capable of maintaining contractile function despite a transient decrease in blood flow and increase in cardiac ATP demand during systole. This study analyzes a previously developed model of cardiac energetics and oxygen transport to understand the roles of the creatine kinase system and myoglobin in maintaining the ATP hydrolysis potential during beat-to-beat transient changes in blood flow and ATP hydrolysis rate.
The theoretical investigation demonstrates that elimination of myoglobin only slightly increases the predicted range of oscillation of cardiac oxygenation level during beat-to-beat transients in blood flow and ATP utilization. In silico elimination of myoglobin has almost no impact on the cytoplasmic ATP hydrolysis potential (ΔG ATPase). In contrast, disabling the creatine kinase system results in considerable oscillations of cytoplasmic ADP and ATP levels and seriously deteriorates the stability of ΔG ATPase in the beating heart.
The CK system stabilizes ΔG ATPase by both buffering ATP and ADP concentrations and enhancing the feedback signal of inorganic phosphate in regulating mitochondrial oxidative phosphorylation.
The working heart relies on uninterrupted supplies of oxygen and substrates to maintain its normal function under different workloads [1, 2]. However, during the systole the heart muscle contracts and the coronary blood flow is greatly reduced; while during the diastole the heart muscle relaxes, and the coronary blood flow approaches a maximum [1, 2].
Myoglobin may work as an oxygen buffer [6, 9], facilitate oxygen diffusion at low cellular oxygen tension [10–12], and/or catalyze chemical reactions (such as NO scavenging) [13–15]. Similarly, the creatine kinase system may either buffer cellular ATP levels or facilitate ATP diffusion inside myocytes [7, 8, 16].
Under normoxic conditions facilitated diffusion of oxygen by oxy-myoglobin is not expected to play a significant role in oxygen transport in the myocardium [6, 17, 18]. Therefore we consider only the oxygen storage function of myoglobin in this computational study. The importance of CK-facilitated high-energy phosphate transport depends on diffusion path length and diffusivity of ATP and ADP in cardiomyocytes . Since the myofibrils generally have small diameters (1–2 μm) and are surrounded by dense mitochondria, it is possible that CK-facilitated transport does not play a significant role in vivo either . The CK-shuttle theory of Saks et al. [8, 19, 20] requires restricted intracellular diffusivity of ADP, for which there is no unambiguous experimental evidence. In this study the buffering role of the CK system is investigated, and CK-facilitated diffusion of high-energy phosphate is not considered.
Here, the roles of myoglobin and creatine phosphate in buffering the energy state (i.e., ATP hydrolysis potential) in the working heart are investigated using a multi-scale computer model integrating cellular metabolism with oxygen transport in the cardiac tissue [21, 22]. The cellular metabolism model was adopted from a recently-published computer model of energetic metabolism of cardiac mitochondria . The metabolic model is integrated with a model of oxygen transport accounting for heterogeneous oxygenation in cardiac tissue .
To examine roles of myoglobin and the CK system, simulations were performed for the normal model (control) and models with myoglobin and CK not active. The simulations demonstrate that the effects of myoglobin on metabolite levels and cytoplasmic ATP hydrolysis potential are barely discernible in the working heart. In contrast, the CK system plays an essential role in maintaining the energetic state in the heart.
Buffer capacity of the CK system
Total creatine pool in myocardium
mmol (l cytoplasm water)-1
Cytoplasmic ATP hydrolysis potential
Average magnitude of cytoplasmic ATP hydrolysis potential
Transformed Gibbs free energy of ATP hydrolysis
Mitochondrial inner membrane potential
ml min-1 (g tissue)-1
Mitochondrial ANT transport flux
mmol sec-1 (liter mito)-1
Cytoplasmic ATP consumption rate
mmol sec-1 (liter cell)-1
Oxygen consumption rate
μmol min-1 (g tissue)-1
Partial oxygen pressure in myocardium
P 50, Mb
Half-saturation partial pressure for the oxygen-myoglobin binding
kJ mol-1 K-1
where P 50, Mb is the half-saturation partial pressure for the oxygen-myoglobin binding, 2.39 mmHg .
where = -34.89 kJ mol-1  at ionic strength I = 0.17 M, and temperature T = 310.15 K, and R is the gas constant equal to 8.314 kJ mol-1 K-1. The subscript "c" denotes cytoplasm.
Oxygen-myoglobin binding is assumed to be maintained in equilibrium in the tissue model. The kinetics of the CK reaction is governed according to the model described in Reference . Molecular diffusion of creatine, phosphocreatine, and myoglobin are not included in the model (see Reference ). Therefore putative roles of oxy-myoglobin and phosphocreatine are not explored in this work.
To understand the roles of myoglobin and the creatine kinase system as buffers in maintaining the energetic state in the heart, we simulate transient responses and steady states of cardiac energetics under different cardiac workloads in three systems: (1) a normal system with both the normal myoglobin level and creatine kinase activity; (2) a system without myoglobin (no-Mb system); (3) a system with zero creatine kinase activity (no-CK system). The following analyses demonstrate the importance of the CK system in stabilizing energetic states in the beating heart. In the normal system, the myoglobin concentration (C Mb) is set to be a physiologically reasonable value 200 μM , and the creatine kinase activity is set to be an arbitrary large value to ensure that the CK reaction is rapid enough to remain near equilibrium . C Mb is reduced to zero in the no-Mb system, and the creatine kinase activity is set to be zero in the no-CK system.
Figure 7A illustrates that the oxygenation levels in the no-CK system vary slightly more than in the control system, while average oxygenation levels in the both system are nearly identical. Figure 7B shows that [ATP]c is maintained at an almost constant level (~9.66 mM) without noticeable oscillations over the whole range of cardiac workload in the control system; in the no-CK system [ATP]c oscillates over a wider range (~± 0.2 mM at the maximum workload), and the average [ATP]c decreases from ~9.66 mM at the baseline workload to ~9.52 mM at the maximum workload in the no-CK system. Figure 7C shows that the average [CrP]c decreases from ~23.4 mM at the baseline workload to ~19.5 mM at the maximum workload in the control system. ([CrP]c is constant in the no-CK system because of zero CK activity.) Figure 7D shows that the average [Pi]c at the maximal work rate is three-folder lower in the no-CK system than in the control system. Average [Pi]c is ~0.7 mM at the maximum workload in the no-CK system compared to ~2.1 mM in the control simulation. Figure 7E shows that [ADP]c oscillates drastically in the no-CK system, where [ADP]c varies between ~85 μm and ~270 μm at the maximum workload, compared with [ADP]c varying between ~60 and ~63 μm at the maximum workload in the control system. As shown in Figure 7F, the average values of ΔG ATPase are approximately equal between the control and no-CK. However, it is apparent from the predicted ranges of oscillation that the CK system acts in stabilizing ΔG ATPase in the myocardium. Over the whole range of cardiac workload, the range of oscillations of ΔG ATPase is approximately four times higher in the no-CK system than in the control. Mitochondrial inner membrane potential (ΔΨ) also oscillates over in a wider range (~± 5 mV) in the no-CK system than in the control system (< ± 1.2 mV) as shown in Figure 7G. Figure 7H shows that in both the control and the no-CK systems the average J ANT at the steady state is equal to the average cytoplasmic ATP hydrolysis rates at different cardiac workloads, indicating the oxidative phosphorylation flux is able to match the cellular ATP demand in both cases. However, the oxidative ATP synthesis rate in the no-CK case oscillates over a much larger range in the no-CK case.
These results demonstrate two roles of the CK system in maintaining cardiac energetic state. The first role as a buffer has been widely appreciated [7, 8, 16, 34]. The second role is to provide a source of increasing inorganic phosphate with increases in work rate. This allows the cardiomyocytes to maintain relatively stable ATP and ADP concentration, as illustrated in Figure 7B and 7E. Briefly, Figure 7 clearly demonstrates significant roles of the CK system in maintaining energetic stability during the heart cycles at varying cardiac workloads.
Comparisons between Figure 4 and 5 demonstrate the oxygen-storage function of myoglobin plays minor roles in maintaining energy state of the heart in normoxic conditions. However, possible physiological roles of myoglobin may be revealed under extreme conditions. The concentration of myoglobin in the heart of most terrain mammals (e.g., rat, dog, human) is typically between 130 and 320 μM [35, 36]. In the event of severe ischemia, these concentrations can only satisfy the oxygen demand of the heart for several seconds even at relatively low work rate. Comparatively, C Mb level is high in diving birds and mammals (e.g., 4~5 mM in seals) , and may contribute to significant oxygen buffering .
First hypothesized over 30 years ago, the relative importance of Mb-facilitated oxygen transport has remained unclear, in part, due to conflicting experimental measurements on myoglobin diffusivity [10–12]. Early conclusions on Mb-facilitated oxygen diffusion were mainly drawn based on high estimates (7~23 × 10-7 cm2 s-1) of myoglobin diffusivity from in vitro measurements in dilute solution of myoglobin [4, 10–12, 39, 40]. Using two different techniques, microinjection and fluorescence recovery after photobleaching (FRAP), Papadopoulos et al. [9, 41] report that that myoglobin diffusivity is approximately 2.0 × 10-7 cm2 s-1 at 37°C in myocardiocytes. Based on this low value of myoglobin diffusivity, Jurgen et al.  and Beard and Bassingthwaighte  determine that the intracellular myoglobin diffusivity may be too low to provide significant facilitation of oxygen transport in well-oxygenated cardiac tissue. Jue and colleagues [43, 44] report higher values of myoglobin diffusion in cardiomyocytes, 7.85 × 10-7 cm2 s-1 at 35°C and 4.24 × 10-7 cm2 s-1 at 22°C, obtained from a 1H-NMR technique. However, even at these diffusivity values, myoglobin does not impact oxygen transport into the cells significantly when cellular oxygen tension is above P 50 of myoglobin [43, 44]. Similarly, Timmons et al. [45, 46] report that oxygen supply does not limit oxidative ATP synthesis in rest-work transition in skeletal muscle, also implying the minor role of myoglobin in oxygen transport in oxidative striated muscle.
Myoglobin can scavenge and catalyze degradation of NO in vivo [13, 47]. In addition, Rassaf et al.  propose that myoglobin in the heart might work as an oxygen sensor and metabolism regulator, since myoglobin can act as either an NO scavenger at high cellular P O2 and an NO producer at low P O2. Since oxygenation in cardiac tissue is maintained above the P 50 for myoglobin in normoxia, myoglobin may act as an NO scavenger under normal conditions.
The absence of myoglobin may elevate cytoplasmic NO level and consequently impacts metabolism regulation, structure, and functions in the heart. Experimental observations show that physiological performance of Mb-knockout mice is comparable to those of normal mice, but the Mb-knockout mice have significantly elevated coronary flow, coronary reserve, and capillary density compared to the control [48, 49]. These differences may be associated with elevated NO levels, resulting from reduced NO scavenging in myocardiocytes of the Mb-knockout mice. Elevated NO may also promote angiogenesis by regulating various growth factors and increase capillary density [50–52]. Flogel et al.  propose that decreased myoglobin may directly affect substrate selection. For example, elevated NO level impacts substrate utilization by effecting gene expressions of the glucose transporter (GLUT4) and the peroxisome proliferator-activated receptor (PPARα) . Hence, abnormally elevated NO levels may contribute to the physiological adaptations observed in Mb-knockout mice.
In summary, our simulations predict that oxygen buffering by myoglobin in terrestrial mammals plays an insignificant role in maintaining the energy state of the beating heart. Significant changes of structure, hemodynamics, and metabolism observed in Mb-knockout mice [49, 53] do not necessarily reflect compensation mechanisms for oxygen deficiency in cardiac tissue. These differences may be caused by elevated cellular NO levels in the knockout mice via various pathways [2, 50–56]. In aerobic muscles, such as cardiac muscle, a large amount of mitochondria are required to provide adequate and sustained energy supply, and this high oxidative capacity is usually accompanied by abundant myoglobin, perhaps functioning primarily to scavenge NO generated by mitochondria.
Simulations presented in Figure 6, 7 show considerably diminished energetic stability in the no-CK system, and thus demonstrate significant roles of the creatine kinase system in maintaining normal function in the working heart. These analyses show that the CK system not only buffers cytoplasmic ATP and ADP concentrations but also enhances the Pi feedback signal in regulating mitochondrial oxidative ATP synthesis.
One way to understand the relative importance of Pi and ADP in regulating oxidative phosphorylation is to clamp their concentrations and simulate the expected impact on the phosphate metabolites. The concentrations of Pi and ADP are clamped by setting their time derivatives to be zero in the model, respectively. Results from these simulation experiments are illustrated in Figure 8. The model predictions of the control system are plotted as solid lines, and those of the no-CK system are illustrated as dashed lines. The previously obtained baseline and maximum work rates (see Methods section), 0.36 and 1.2 mmol s-1 (liter cell)-1, respectively, are indicated by two vertical dotted lines in the plots in Figure 8. The simulations without clamped metabolite concentrations are demonstrated in Figure 8A–C. Figure 8A shows that the magnitudes of average ΔG ATPase ( ) are approximately equal between the control and no-CK over the range of workload. Figure 8B illustrates that the average [Pi]c only increases approximately two fold over the range of work rate in the no-CK system, compared to the approximately eight-fold increase in the control. Figure 8C shows that the average [ADP]c increases approximately four fold over the range of work rate, compared to the only slight increase of the average [ADP]c in the control. With [Pi]c clamped at the baseline level, the model-predictions are the same in both the control and no-CK system, as shown Figure 8D–F. Figure 8D shows that reaches the critical value , 63.5 kJ mol-1, at J ATPase = ~0.90 mmol s-1 (liter cell)-1. Figure 8F shows a rapid increase of the average [ADP]c in responses of elevated work rate. Clamping [ADP]c impacts the simulations less significantly, as shown in Figure 8G–I. Figure 8G shows that with [ADP]c clampled, reaches the critical value at J ATPase = ~1.04 mmol s-1 (liter cell)-1 in the control, and at J ATPase = ~1.08 mmol s-1 (liter cell)-1 in the no-CK system. With [ADP]c clamped (Figure 8H) and with [ADP]c free to vary (Figure 8B), Pi concentrations are higher in the control system than in the no-CK system because of a larger effective total pool of exchangeable phosphates in the control.
To summarize these findings, clamping [Pi]c has a more significant impact on energetic stability than clamping [ADP]c in both the control and no-CK systems. For both systems, when Pi feedback is removed from the model (by clamping [Pi]c) the free energy of ATP hydrolysis drops more quickly with increasing work rate than in the unclamped simulations. When [ADP]c is clamped, the impact on the predicted energetic state is less significant than when [Pi]c is clamped. Therefore, inorganic phosphate is a more important controller of oxidative phosphorylation than ADP in this model. Since this model is well validated against steady-state and transient in vivo data, it is likely that inorganic phosphate is a primary controller of oxidative phosphorylation in vivo.
Previous analysis of data from purified mitochondria suggests that Pi activates complex III of the respiratory chain . Indeed this phenomenon is incorporated into our computational model [21, 23]. When Pi level decreases in the no-CK system, the mitochondrial membrane potential (Figure 6G) is diminished, leading to increased [ADP]c and decreased |ΔG ATPase| compared to control. However, even though the Pi signal is diminished in the no-CK system, Pi remains the most important feedback signal for oxidative phosphorylation.
Saupe et al.  observed that as work rate increases in hearts isolated from mice with both mitochondrial and M-form CK knocked out, [Pi]c changes little, [ADP]c increases about three fold, and decreases ~3 kJ mol-1 in these hearts compared to the control. Thus, their concentrations on ADP are qualitatively matched by our simulations, while our model predicts that [Pi]c is lower in the no-CK system compared to the control. Since the knockout mouse cardiomyocytes show residual CK activity (40% of wild type), this animal model is not equivalent to our no-CK model. The data of Saupe et al.  show a decrease in CrP with increasing work rate, indicating that the CK system is potentially active in the knockout animals.
Feedback of ADP and Pi is also essential for matching oxidative ATP synthesis to cellular energy demand in skeletal muscle [23, 59]. However, cytoplasmic Pi can increase to 20 mM and higher at high work rates in skeletal muscle  while cytoplasmic Pi is predicted to stay below 3 mM at maximal work rate in the heart . Furthermore, the cytoplasmic Pi concentration in resting slow oxidative soleus muscle is in the range of 5 mM . Therefore even at rest, the Pi concentration is well above the predicted regulatory feedback range for cardiac muscle. We would expect that ADP acts as an important physiological feedback signal in those muscles, as has been established [63, 64].
In additional to a temporal buffering function, the CK system has been proposed to facilitate transport of energetic phosphates in myocardium [7, 19]. This hypothesis remains controversial. Meyer et al.  describe both temporal buffering and "spatial buffering" roles of creatine phosphate associated with the near-equilibrium creatine kinase (CK) reaction and a high cytoplasmic ATP-to-ADP ratio. The spatial buffering role may be negligible in the heart because of small diameters of cardiomyofibrils and abundant surrounding mitochondria . The phosphocreatine shuttle hypothesis – that the free energy of ATP hydrolysis is transported primarily by spatial gradients of CrP and Cr between mitochondria and sites of ATP hydrolysis – hinges on the existence of three critical phenomena: (1) restricted diffusion of adenine nucleotides in cardiomyocytes; (2) functional coupling (direct product-substrate channeling) between mitochondrial adenine nucleotide translocase and creatine kinase; and (3) disequilibrium of the CK reaction in cardiomyocytes [67–70]. While computational models that invoke these three phenomena are able to match data on the kinetics of oxidative phosphorylation in isolated skinned cardiomyocytes  and purified mitochondria , it remains to be demonstrated that the available data, particularly in vivo data, cannot be explained without invoking these phenomena. Indeed, recent measurements of diffusivities of labeled adenine nucleotides in cardiomyocytes demonstrate that bulk diffusion is not restricted to the degree necessary for the phosphocreatine shuttle to operate . The existence and significance of "microcompartments"  around mitochondria with restricted diffusion in cardiomyocytes remains an active subject of investigation and debate.
To determine the roles of myoglobin and the CK system in stabilizing cardiac energy state, a computational model of oxygen transport and cardiac metabolism is applied to simulate transient changes and steady states in the beating heart. The analysis suggests that myoglobin has little impact on cardiac energetics, while the CK system is important for the beating heart to maintain stable energetic state over a range of cardiac work rate. Two distinct functions of the CK system are apparent from this analysis: first, the CK system buffers changes of [ADP]c and [ATP]c and ΔG ATPase; second, the CK system enhances the feedback of Pi to match the rate of mitochondrial ATP synthesis to cellular ATP demand while maintaining relatively stable [ADP]c.
The authors thank Drs. Kalyan C. Vinnakota, Feng Qi, and Martin J. Kushmerick for helpful discussions. This work was supported by NIH grant HL072011 (DAB).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.