Am J Physiol Heart Circ Physiol 291: H1036-H1046, 2006.
First published April 7, 2006; doi:10.1152/ajpheart.01382.2005
0363-6135/06 $8.00
TRANSLATIONAL PHYSIOLOGY
Regulation of myocardial substrate metabolism during increased energy expenditure: insights from computational studies
Lufang Zhou,1,5
Marco E. Cabrera,1,2,3,5
Isidore C. Okere,2,5
Naveen Sharma,4,5 and
William C. Stanley2,4,5
Departments of Biomedical Engineering,1 Physiology and Biophysics,2 Pediatrics,3 Nutrition,4 and the Center for Modeling Integrated Metabolic Systems,5 Case Western Reserve University, Cleveland, Ohio
Submitted 29 December 2005
; accepted in final form 27 March 2006
 |
ABSTRACT
|
|---|
In response to exercise, the heart increases its metabolic rate severalfold while maintaining energy species (e.g., ATP, ADP, and Pi) concentrations constant; however, the mechanisms that regulate this response are unclear. Limited experimental studies show that the classic regulatory species NADH and NAD+ are also maintained nearly constant with increased cardiac power generation, but current measurements lump the cytosol and mitochondria and do not provide dynamic information during the early phase of the transition from low to high work states. In the present study, we modified our previously published computational model of cardiac metabolism by incorporating parallel activation of ATP hydrolysis, glycolysis, mitochondrial dehydrogenases, the electron transport chain, and oxidative phosphorylation, and simulated the metabolic responses of the heart to an abrupt increase in energy expenditure. Model simulations showed that myocardial oxygen consumption, pyruvate oxidation, fatty acids oxidation, and ATP generation were all increased with increased energy expenditure, whereas ATP and ADP remained constant. Both cytosolic and mitochondrial NADH/NAD+ increased during the first minutes (by 40% and 20%, respectively) and returned to the resting values by 1015 min. Furthermore, model simulations showed that an altered substrate selection, induced by either elevated arterial lactate or diabetic conditions, affected cytosolic NADH/NAD+ but had minimal effects on the mitochondrial NADH/NAD+, myocardial oxygen consumption, or ATP production. In conclusion, these results support the concept of parallel activation of metabolic processes generating reducing equivalents during an abrupt increase in cardiac energy expenditure and suggest there is a transient increase in the mitochondrial NADH/NAD+ ratio that is independent of substrate supply.
diabetes; exercise; heart; lactate; mitochondria; modeling
CARDIAC PUMP FUNCTION is fueled by ATP hydrolysis, which is precisely matched by ATP formation, primarily in the mitochondria (42). In the transition from rest to intense exercise, there is a three- to sixfold increase in the rate of cardiac power generation, myocardial oxygen consumption (M
O2), and ATP turnover (25). Nevertheless, at high work states the myocardial ATP and ADP concentrations are maintained at a relatively constant level (2, 4, 36). There is rapid activation of NADH generation in the mitochondria, flux through the electron transport chain (ETC), and ATP production by oxidative phosphorylation to match exactly ATP breakdown in the cytosol. Studies in isolated mitochondria show that the generation of NADH from carbon substrates, oxygen consumption, and oxidative phosphorylation are turned on by feedback from an increase in ADP concentration (9); however, the regulatory mechanisms in vivo are still unclear. Experimental studies suggest that the mitochondrial NADH and NAD+ are also maintained near steady-state values at high rates of energy expenditure due to the matching of NADH oxidation by complex I of the ETC with the rapid activation of substrate metabolism and NADH formation by dehydrogenases [e.g., pyruvate dehydrogenase (PDH), fatty acids oxidation, and the Krebs tricarboxylic acid cycle] (22, 36). The stability of these regulatory species in the heart demonstrate that metabolic control models based solely on changes in the levels of these intermediates to control metabolic flux rates cannot fully explain the observed responses (5, 2729).
An alternative regulatory hypothesis proposes that there is simultaneous parallel activation of ATP utilization and production, glycolysis, and mitochondrial dehydrogenases by signaling molecules, e.g., Ca2+ (6, 7, 31, 43). Simulations from an integrated computational model of cardiac mitochondrial energy metabolism demonstrate that an increase in Ca2+ concentration in the mitochondrial matrix activates NADH generation and oxidative phosphorylation (13). Indeed, Korzeniewski et al. (2729) found additional support for the parallel-activation mechanism in the skeletal muscle and the heart by using a dynamic computer model of oxidative phosphorylation; however, this model did not incorporate metabolic pathways for substrate metabolism and NADH generation. Because current experimental measurements of myocardial metabolism lump the cytosol and mitochondria, it is important to distinguish these subcellular compartments in computational models to provide better insight into the metabolic responses during the transition to increase cardiac power.
The plasma concentration of substrates effect the metabolic response to increased cardiac work such as with the rise in arterial lactate concentration during physical exercise. Myocardial lactate uptake during exercise is increased as a function of the arterial lactate concentration (16, 23, 38), which is dependent on the rate of systemic lactate production and clearance (39). In addition, studies in resting dogs (37) and perfused rat hearts (15) demonstrate that an elevation in arterial lactate concentration reduces the relative contribution of fatty acids to M
O2 (37). The mechanism for this effect is not clear but could be due to a lactate-induced increase in cytosolic and/or mitochondrial NADH/NAD+ and inhibition of mitochondrial fatty acid oxidation through NADH inhibition of
-oxidation. Plasma substrate levels are also altered with diabetes, as are the expression and activity of key enzymes and transporters. The metabolic adaptations to diabetes are particularly complex because there is both altered substrate availability (e.g., increased plasma ketone bodies and glucose) and altered expression of key proteins (e.g., phosphofructokinase, pyruvate dehydrogenase kinase, and glucose transporters) that effect myocardial substrate selection. During increased cardiac power we observed that diabetic pigs have decreased PDH activity and lactate uptake and enhanced
-hydroxybutyrate uptake compared with healthy animals (19). It is generally postulated that impaired pyruvate oxidation in the diabetic heart is partially due to an increase in mitochondrial NADH/NAD+ (19, 32, 33, 40); however, evidence for this is lacking due to the inability to measure mitochondrial NADH and NAD+ in the intact heart.
The overall goal of the present study was to 1) assess the role of parallel activation of multiple metabolic steps in the regulation of cytosolic and mitochondrial NADH/NAD+ during the dynamic transition from resting conditions to high rates of cardiac energy expenditure, as observed with exercise; 2) investigate the effects of elevated arterial lactate concentration on mitochondrial pyruvate and fatty acid oxidation and the cytosolic and mitochondrial NADH/NAD+ ratios; and 3) determine whether impaired pyruvate oxidation under diabetic conditions requires an increase in mitochondrial NADH/NAD+. We modified our recently published computational model of cardiac metabolism (8, 46, 47) by incorporating multistep activation mechanisms (29) and simulated the metabolic responses of the heart to an abrupt increase in energy expenditure. We hypothesized that in response to an abrupt increased cardiac energy expenditure, the simultaneous activation of ATP hydrolysis, ATP synthesis, glycogen breakdown, NADH generation, and electron transport chain activity would maintain constant ATP and ADP levels but cause a sharp increase in the cytosolic NADH/NAD+ ratio, with less effect on mitochondrial NADH/NAD+. On the other hand, elevation in arterial lactate concentration simultaneously with the increase in energy expenditure would increase the cytosolic NADH/NAD+ ratio and mitochondrial pyruvate oxidation and reduce fatty acid oxidation, with little effect on mitochondrial NADH/NAD+ ratio. Finally, we hypothesized that simulation of diabetic conditions would reduce pyruvate oxidation but would not require an increase in mitochondrial NADH/NAD+.
 |
METHODS
|
|---|
Model development.
To quantitatively investigate the effects of increased energy expenditure on cardiac energy metabolism and potential regulatory mechanisms, we modified our previously published cardiac model of ischemia (48) by including multistep activation of specific metabolic processes. As in our ischemia model (47, 48), the dynamic mass balance equation for species j in a spatially lumped capillary blood domain (accounting for fractional mixing between the arterial and venous blood) was described as:
 | (1) |
where Caj is the arterial concentration of species j, Cbj is the blood concentration, Vb is the domain volume, Q is blood flow, 0 < F
1 is the mixing factor, and Jb-c,j is the mass transport rate between blood and cytosol.
The dynamic mass balance of species j in the spatially lumped cytosolic domain was described by the involved metabolic reactions in the domain and mass transfer across the cellular membrane and mitochondrial membrane:
 | (2) |
where Ccj is the concentration, Vcj is the effective volume of species j in the cytosol, Rcj is the metabolic reaction rate of species j, and Jcm,j is the mass transport rate between cytosol and mitochondria.
In the mitochondrial domain, the dynamic mass balance equation for species j is
 | (3) |
where Cmj is the concentration, Vmj is the effective volume of species j in the mitochondrial domain, and Rmj is the net reaction rate.
The dynamic mass balance for protons in the mitochondrial intermembrane space is
 | (4) |
where Vi is the effective volume of intermembrane space, Jfm
i,H is the proton efflux from mitochondrial matrix to intermembrane space, Jfi
m,H is the proton influx, and Jleak is the proton leak flux through the mitochondrial inner membrane. The proton fluxes values under steady state and the corresponding parameters can be found in our recently published paper (47).
The reaction and transport kinetics terms in the mass balance equations (e.g., glycolysis, oxygen consumption, ATP synthesis, glucose uptake, etc.) were described by Michaelis-Menten equations as in our ischemia model. To simulate increased energy demand in the model of cardiac cellular metabolism, we assumed simultaneous activations of multistep metabolic fluxes in both cytosol and mitochondria (Fig. 1), as well as increase in myocardial blood flow (MBF). Consequently, simultaneous activation occurred with glucose uptake, glycogen breakdown, glycolysis, ATP hydrolysis, dehydrogenase, Krebs tricarboxylic acid cycle, oxygen consumption (complex I), and ATP synthesis (complex V) (35). This activation was implemented in the model by increasing the maximal velocity (Vmax) of metabolic fluxes n times through the Michaelis-Menten coefficients for corresponding processes that are proportional to the increase in energy expenditure:
 | (5) |
The results of this multistep activation mechanism will be evident in model simulations, which correspond closely with published experimental data in pigs (36) involving oxygen consumption, glucose oxidation, and fatty acid oxidation with increased cardiac work induced by a simultaneous injection of atropine, infusion of dobutamine, and constriction of the proximal thoracic aorta.

View larger version (30K):
[in this window]
[in a new window]
|
Fig. 1. Myocardium metabolic pathway. This complex network incorporates most key biochemical reactions and pathways involved in cardiac metabolism, including glycolysis, pyruvate oxidation, fatty oxidation, the tricarboxylic acid cycle, and oxidative phosphorylation. It also includes species transport across cellular membrane and mitochondrial membrane, as well as proton influx and efflux between mitochondrial matrix and inter-membrane space. GLU, glucose; G6P, glucose-6-phosphate; GLY, glycogen; GAP, glyceraldyhe-3-phosphate; BPG, 1,3-bisphosphate-glycerate; PYR, pyruvate; LAC, lactate; TG, triglyceride; GLR, glycerol; FFA, free FA; FAC, fatty acyl-CoA; PCr, phosphocreatine; Cr, creatine; CIT, citrate; -KG, -ketoglutarate; SCA, succinyl-CoA; SUC, succinate; MAL, malate; OAA, oxaloacetate; ACoA, acetyl-CoA; bADP, bound ADP; ADP, free ADP; BHB, -hydroxybutyrate. *Activation of metabolic process with onset of increase energy expenditure.
|
|
To simulate the dynamic responses of the diabetic heart to increased energy expenditure, we incorporated
-hydroxybutyrate in the present model (Fig. 1) because in the diabetic heart, there is an elevated uptake and oxidation of ketone bodies (primarily
-hydroxybutyrate), which is concentration dependent (19, 41). We assumed that the concentration of
-hydroxybutyrate under normal conditions is very low, and it increases with diabetes, therefore, the rate of
-hydroxybutyrate oxidation increases.
Parameter estimation and simulation strategy.
The initial values of the reaction flux rates, transport flux rates, and species concentrations of this model were modified so that they matched the published experimental values collected in pigs during the baseline period (26, 36). The species concentrations (in the blood, cytosol, and mitochondria) were obtained by taking into account the relative distribution of these species in the specific subcellular domains as in our previous models and are listed in Table 1. Values for the reaction and transport rates and associated parameter values that were updated in this model are listed in Tables 2 and 3; the remainder of the parameter values were the same as in our previously published model (8, 46, 47).
Simulations were oriented to investigate the metabolic responses and the regulation of metabolic processes in the cytosol and mitochondria during the transition from resting to conditions of a fourfold exponential increase (
= 0.3 min) in energy expenditure and MBF. These simulations emphasized the time profiles of key metabolite concentrations (e.g., ATP, ADP, NADH, NAD+, glycogen, lactate, etc.) and metabolic fluxes (e.g., oxygen consumption, glycolysis, ATP hydrolysis, FAO, etc.) and were compared with previously published experimental data in pigs. The model was then used to simulate the metabolic responses of the heart to the increased energy expenditure under different metabolic milieus: 1) constant arterial levels of substrates [reflecting the baseline conditions in our recent studies in anesthetized pigs (26, 36)]; 2) elevated arterial lactate concentrations [0.7 mM (normal), 2.1 mM, and 4.2 mM] initiated with the onset of increased myocardial energy expenditure; and 3) diabetic conditions (high
-hydroxybutyrate + high arterial glucose + inhibited glucose uptake). Computer simulations were obtained by solving the model equations by using DLSODE, a robust implicit integrator for stiff and sparse systems (CASC, 2003; http://www.llnl.gov/CASC/odepack/software/dlsode.f).
 |
SIMULATION RESULTS
|
|---|
Elevated metabolic fluxes during increased energy expenditure.
With the increase of cardiac ATP hydrolysis, there were significant and immediate increases in M
O2 and ATP generation through oxidative phosphorylation (Fig. 2). M
O2 increased 3.8-fold (from 2.7 to 10.1 µmol·g1·min1) and then decreased slightly to 9.9 µmol·g1·min1 at 15 min. Consistent with increased M
O2, ATP generation rate increased approximately fourfold from 14.8 to 58.8 µmol·g1·min1 and then decreased to 57.5 µmol·g1·min1. ATP hydrolysis rate dynamics followed similar pattern (increased from 15.6 to 63.7 µmol·g1·min1 then decreased to 61.4 µmol·g1·min1) in response to increased cardiac energy expenditure.
The increase in energy expenditure also stimulated glycolysis, pyruvate oxidation, and
-oxidation (Fig. 2). The glycolytic rate increased from its resting state value of 0.12 µmol·g1·min1 to a peak value of 0.91 µmol·g1·min1 (7.6-fold) then decreased to 0.63 µmol·g1·min1 (5.3-fold) as glycogen was depleted (Fig. 2C). The pyruvate oxidation rate increased to a peak value of 1.88 µmol·g1·min1 and then decreased to 1.38 µmol µmol·g1·min1 (Fig. 2C). Fatty acid
-oxidation increased from 0.07 to 0.24 µmol·g1·min1at 15 min of increased energy expenditure (Fig. 2C).
Species concentration dynamics in response to increased cardiac energy expenditure.
The myocardial glycogen concentration decreased exponentially from its initial value of 36 to 15.9 µmol/g at 15 min of increased energy expenditure, whereas the lactate concentration increased from 2.6 to 5.9 µmol/g (Fig. 3A). Although the ATP concentration in the mitochondria slightly decreased, the cytosolic ATP concentration remained constant (Fig. 3).

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 3. Model simulated dynamics of glycogen and lactate (A), cytosolic ATP and ADP (B), and mitochondrial ATP and ADP (C) concentrations in response to fourfold increased energy expenditure. Solid circles and squares represent experimental data on lactate and glycogen concentrations, respectively.
|
|
There was a small transient increase in NADH and decrease in NAD+ in both the cytosol and mitochondria. The changes were greater in the cytosolic than in the mitochondrial compartment, with the peak increase in the NADH/NAD+ ratio being a 40% increase in the cytosol and 20% in the mitochondria, with a return to baseline values by 15 min of increased cardiac energy expenditure (Fig. 4). To compare simulation results with experimental data, we calculated the tissue NADH and NAD+ concentrations from their cytosolic and mitochondrial values assuming the volume ratio of cytosol to mitochondria to be 70:17 (46). The result shows that the simulated tissue NADH and NAD+ concentrations and NADH/NAD+ ratio matched the tissue measurements made at 15 min of increased cardiac energy expenditure; however, sampling late in the stimulation period missed the predicted transient increase in NADH/NAD+ ratio at the onset of increased cardiac energy expenditure (Fig. 4C).

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 4. Model simulated dynamics of NADH and NAD+ concentrations in the cytosol (A), mitochondria (B), and whole tissue (C) in response to increased energy expenditure. Solid triangle and circles represent experimental data on tissue NAD+ and NADH concentrations, respectively. Insets: NADH/NAD+ ratios plotted on the same time scale.
|
|
Effects of pathway activation on metabolic homeostasis.
The mechanisms responsible for the maintenance of regulatory species (high-energy phosphates and nicotinamide adenine dinucleotides) were further investigated by systematically eliminating the activation of specific metabolic processes with the increase in cardiac energy expenditure (Fig. 5). Specifically, we simulated the cardiac responses to increased energy expenditure with 1) activation of all pathways as in Fig. 1; 2) activation of all pathways except glycolysis and glycogen breakdown; 3) activation of all pathways except mitochondrial dehydrogenases; 4) activation of all pathways except complex I; 5) activation of all pathways except complex V; and 6) activation of all pathways except transferring of NADH-NAD+ and ATP-ADP between the cytosol and mitochondria (denoted as "communication" in Fig. 5).

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 5. Effects of pathway activation on metabolic homeostasis in response to increased energy expenditure. A: mitochondrial NADH/NAD+; B: cytosolic NADH/NAD+; C: mitochondrial ATP/ADP; D: cytosolic ATP/ADP.
|
|
Prevention of the activation of glycolysis resulted in a major reduction in the cytosolic redox and the phosphorylation potential but had more modest effects in the mitochondria (Fig. 5). Preventing the activation of glycolysis resulted in a decrease in the cytosolic and mitochondrial NADH/NAD+ ratio when energy expenditure was increased instead of the increase that was observed with parallel activation (Fig. 5). There was a modest decrease in the ATP/ADP ratio in both compartments when the activation of glycolysis was prevented and a greater effect in the cytosol than in the mitochondria (Fig. 5). In addition, the steady-state M
O2 was decreased from 10.1 to 9.1 µmol·g1·min1 when activation of glycolysis was blocked (data not shown).
Blockade of the activation of mitochondrial dehydrogenases had profound effects on both mitochondrial and cytosolic redox, the phosphorylation potential, and the increase in M
O2, suggesting that this process is critical for maintaining metabolic homeostasis in the transition from normal to high rates of energy expenditure. Simulation results show that preventing the activation of mitochondrial dehydrogenases resulted in a
50% decrease in both cytosolic and mitochondrial NADH/NAD+ and ATP/ADP ratios (Fig. 5). The increase from baseline in steady-state M
O2 was reduced by 74% (4.6 µmol·g1·min1) when the activation of mitochondrial dehydrogenases was prevented.
Preventing the activation of complex I caused a large increase in the cytosolic and mitochondrial NADH/NAD+ ratio, a 40% decrease in the mitochondrial ATP/ADP ratio, and a 57% reduction the increase from baseline in steady-state M
O2 (5.8 µmol·g1·min1). In contrast to blockade of the activation of dehydrogenases, preventing the activation of complex I caused modest fall in the cytosolic ATP/ADP ratio. Blocking the activation of complex V had minor effects on metabolic homeostasis (Fig. 5) and had no effect on the increase in M
O2 (data not shown). Blocking activation of complex V had little effect on the mitochondrial and cytosolic NADH/NAD+ or ATP/ADP ratios compared with parallel activation, suggesting that this step is less important than activation of mitochondrial dehydrogenases or complex I for maintaining metabolic homeostasis during the transition for normal to high rates of energy expenditure.
The role of communication between the cytosol and mitochondrial was investigated by preventing the activation of the exchange of ATP ADP, and NADH NAD+ between the cytosol and mitochondria. This resulted in a slight decrease in mitochondrial NADH/NAD+ and a moderate decrease in cytosolic ATP/ADP but a huge increase in mitochondrial ATP/ADP ratio and cytosolic NADH/NAD+. The steady-state increase in M
O2 was lowered by only 6% when the activation of the exchange of ATP ADP and NADH NAD+ was prevented.
Effects of elevated arterial lactate.
In silico experiments were performed to investigate the role of elevated arterial lactate concentration on cardiac metabolism during the transition from resting to increased energy expenditure conditions (Fig. 6). Simulations were performed with a constant arterial lactate concentration of 0.7 mM (as in Figs. 25) and with an increase in arterial lactate to either 2.1 or 4.2 mM (Fig. 6) to mimic the increases in arterial lactate observed in human exercise. Arterial lactate concentration had no effect on M
O2, the rate of oxidative phosphorylation, and the concentrations of ATP, ADP, or phosphocreatine (data not shown). Increased arterial lactate concentration increased pyruvate oxidation and reduced myocardial fatty acid oxidation (Fig. 6C). Cytosolic NADH/NAD+ ratios (
4090% at 2.1 mM and
90%-140% at 4.2 mM) were increased as a function of elevations in arterial lactate concentration, reflecting a major effect on lactate uptake and the lactate dehydrogenase reaction (Fig. 6). On the other hand, the mitochondrial NADH/NAD+ ratio was not affected significantly by elevated arterial lactate concentration during the first 5 min but was 1015% higher with elevated arterial lactate at 16 min of increased work, reflecting a higher cytosolic NADH/NAD+ ratio and elevated pyruvate oxidation compared with low arterial lactate concentrations.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 6. Model simulated dynamics of arterial lactate (LAC) concentration (A), myocardial lactate uptake (B), fatty acid oxidation (C), pyruvate oxidation (D), cytosolic NADH/NAD+ ratio (E), and mitochondrial NADH/NAD+ (F) in response to increased energy expenditure.
|
|
Dynamic responses of diabetic heart to increased energy expenditure.
Conditions in the diabetic myocardium were simulated by increasing the arterial concentrations of
-hydroxybutyrate from 0.03 to 0.8 mM and glucose from 4.75 to 9.5 mM and reducing the maximal velocities for glucose transport by 50% and PDH activity by 25% with normal energy expenditure (Fig. 7, phase I). There was a decrease in pyruvate oxidation (from 0.37 to 0.31 µmol·g1·min1) and a slight increase in fatty acids oxidation (from 0.07 to 0.078 µmol·g1·min1), and the heart began to utilize
-hydroxybutyrate as a fuel (Fig. 7). In response to a fourfold increase in energy expenditure, the peak value of pyruvate oxidation under diabetic conditions was 1.2 µmol·g1·min1, which was lower than that in the normal heart (1.9 µmol·g1·min1) (Fig. 7B, phase II). On the other hand, the rate of fatty acids oxidation under diabetic conditions was increased compared with the nondiabetic conditions (Fig. 7C), and
-hydroxybutyrate oxidation rate increased dramatically to 0.18 µmol·g1·min1 (Fig. 7A). Diabetes decreased the cytosolic NADH/NAD+ ratio and had no effect on the mitochondrial NADH/NAD+ ratio (Fig. 7D) nor on the rates of ATP generation or M
O2 (data not shown).
 |
DISCUSSION
|
|---|
The results of the present investigation add further support to the concept that an abrupt increase in cardiac energy expenditure results in a simultaneous activation of metabolic processes in both cytosol (ATP hydrolysis, glycogen breakdown, and glycolysis) and mitochondria (dehydrogenase activity, NADH oxidation, oxidative phosphorylation of ADP) (4, 21, 28, 29, 43, 45). This parallel activation has been demonstrated to be essential for maintaining metabolic homeostasis (e.g., ATP, ADP, and Pi) in response to increased energy demand, as shown in Fig. 5, supporting previous findings by other investigators (3, 27). The results also show the novel finding of a transient increase in the NADH/NAD+ ratio in both the cytosolic and mitochondrial compartments, with a peak increase of a 40% increase in the cytosol over baseline values but only a 20% increase in the mitochondria. On the other hand, we found that an elevation in arterial lactate concentration with the onset of the increase in energy expenditure increased the cytosolic NADH/NAD+ ratio and mitochondrial pyruvate oxidation and reduced fatty acid oxidation but had little effect on the mitochondrial NADH/NAD+ ratio. Finally, our simulation results under diabetic conditions showed that the dynamics of the mitochondrial NADH/NAD+ ratio was not altered during increase energy expenditure despite enhanced oxidation of
-hydroxybutyrate and fatty acids and reduced pyruvate oxidation. It is important to point out that in the present model, the total ATP+ADP+AMP and NADH+NAD+ pools were conserved; however, there were no constraints to maintain constant ATP/ADP and NADH/NAD+ ratios.
Model validation.
The findings of the present study cannot be obtained using current experimental methods, which clearly demonstrates the utility of computational models in the investigation of the regulatory mechanism of cellular energy metabolism. We can, however, validate our computational model by comparing some of the simulation results with data we recently obtained in anesthetized pigs subjected to a fourfold increase in myocardial workload for 17 min (26, 36). Our model simulations corresponded closely with available experimental data on M
O2 (36) (Fig. 2A) and the concentrations of glycogen, lactate (Fig. 3A), tissue NADH, and NAD+ (Fig. 4C). Accompanying the increased energy expenditure was an increase in pyruvate oxidation (4.7-fold) (Fig. 5A) and fatty acid oxidation (3.2-fold) (Fig. 2C), even in the absence of an increase in arterial lactate concentration, which is similar to the relative increase we previously measure for oxidation of exogenous glucose and fatty acids in pigs subjected to a three- to fourfold increase in M
O2 (20, 26, 36). Thus the heart increases its reliance on carbohydrate oxidation with increased cardiac work, as demonstrated in previous experimental studies in humans and pigs (16, 30, 36) and in the isolated perfused rat heart (12, 17). Moreover, model simulations predicted a greater rate of lactate uptake and pyruvate oxidation as a function of increased arterial lactate concentration, which resulted in inhibition of fatty acid oxidation. The greater rate of pyruvate oxidation with elevated lactate corresponded with an increase in the cytosolic NADH/NAD+ ratio, without significant effect on the mitochondrial NADH/NAD+ ratio (Fig. 6), suggesting that alterations in mitochondrial NADH/NAD+ are not essential for lactate inhibition of fatty acid oxidation at high rates of energy expenditure.
Parallel activation of metabolic pathways in response to increased energy expenditure.
There is growing evidence that the fine regulation of myocardial ATP concentration during the transition from resting conditions to increased cardiac workload is due to the parallel activation of multiple metabolic steps concomitant with the increase in contractile power and ATP hydrolysis (28, 29). ATP formation is driven by the activity of complex V, availability of ADP and Pi, the proton motive force across the mitochondrial inner membrane, the activity of the complexes of the ETC, and the delivery of reducing equivalents (NADH and FADH2) to complexes I and II. The results of the present investigation are consistent with evidence from previous experimental studies (3, 4, 6, 34, 44) and computational models (28, 29). The mechanisms for the parallel activation remains somewhat unclear but can be at least partially attributed to Ca2+ activation of glycolysis, mitochondrial dehydrogenases, and complex V (31, 34, 42, 44). One possible signaling source between aerobic metabolism and myocardial work is Ca2+ (3) because it not only regulates myosin ATPase and sarcoplasmic reticulum Ca2+-ATPase but also directly regulates cytosolic ATP hydrolysis, mitochondrial dehydrogenases (3, 6, 31, 43, 44), and complex V (34, 43). Furthermore, in vitro studies show that Ca2+ can also activate PDH phosphatase (therefore activates PDH) and the tricarboxylic acid cycle enzymes isocitrate dehydrogenase and
-ketoglutarate dehydrogenase, stimulating NADH generation (31). Moreover, it has been demonstrated that there is an increase in mitochondrial Ca2+ when there is an increase in extramitochondrial Ca2+, such as that occurs with adrenergic stimulation (31). An additional activator is an increase in free Pi, which can activate complex I and oxidative phosphorylation during the transition to high rates of energy expenditure (6). Another possibility is that there is direct allosteric activation of complex I, as suggested by the observation that complex I in heart mitochondria is phosphorylated in a cAMP-dependent manner (11).
To investigate quantitatively the extent of control exerted by various components of "parallel activation" in maintaining the metabolic homeostasis observed in the present study and other experimental studies, we simulated the cardiac responses to increased energy expenditure without activation of specific pathways (Fig. 5). We found that only when all the components of "parallel activation" investigated were activated simultaneously and to the same extent, the ADP and ATP concentration in both cytosol and mitochondria remained nearly constant. Without activating glycolysis and mitochondrial dehydrogenases, the NADH/NAD+ ratios in both cytosol and mitochondria decreased, which makes sense because these pathways are the major producers of reducing equivalents. This also decreased the ATP/ADP ratio in the cytosol and mitochondria because, as we discussed above, mitochondrial ATP resynthesis relies on a steady-state supply of reducing equivalents. Hindrance of complex I blocked the consumption of NADH causing an increase in NADH and a decrease in ATP due to lowering the H+ gradient across the inner mitochondrial membrane and inhibiting ATP synthesis (Fig. 5, C and D). Activation of complex V did not affect the NADH/NAD+ and ATP/ADP ratios significantly, indicating that this component is not the key site for the regulation of metabolic homeostasis under these conditions. Whereas we have found that with ischemia the metabolic communication between cytosol and mitochondria is impaired (47), in the present investigation we show that active metabolic communication is essential for maintaining constant concentrations of energy transfer species during the transition from normal condition to increased energy expenditure. Without the transferring systems activated, mitochondrial ATP/ADP increased significantly, whereas cytosolic ATP/ADP decreased. This is not surprising, considering that more than 96% of ATP is produced inside the mitochondria, whereas almost all ATP is hydrolyzed inside the cytosol to support cardiomyocytes contraction. From these results, we concluded that 1) mitochondrial ATP/ADP and cytosolic NADH/NAD+ are mostly affected by the transferring of ATP-ADP and NADH-NAD+, 2) mitochondrial NADH/NAD+ is mainly regulated by complex I and mitochondrial dehydrogenases, and 3) cytosolic ATP/ADP is regulated by both mitochondrial dehydrogenases and complex I.
Substrate selection in diabetic heart in response to increased energy expenditure.
One established feature of the diabetic heart is altered myocardial substrate utilization [e.g., increased fatty acids and/or ketone body utilization, and decreased glucose uptake and pyruvate oxidation, and utilization (41)]; however, the effect of increased workload on substrate metabolism and the corresponding regulatory mechanism in diabetic hearts are not well understood. We observed enhanced
-hydroxybutyrate and fatty acids oxidation and depressed pyruvate oxidation under conditions of increased energy expenditure (Fig. 7). Previous studies in diabetic hearts show reduced PDH activity and pyruvate oxidation, greater ketone body oxidation, and normal or enhanced fatty acid oxidation (1, 10, 14, 18). The results of the present investigation provide evidence that impaired pyruvate oxidation in the diabetic heart is not due to an increase in mitochondrial NADH/NAD+ (Fig. 7), as previously suggested by measurements in isolated mitochondria from diabetics rats (24).
Model limitations and future directions.
Whereas the current model is able to predict the cardiac responses to increased energy expenditure and provide useful insight into metabolic regulatory mechanism, there are several limitations that need to be addressed. The present model does not link substrate metabolism and ATP hydrolysis with mechanical power generation, cardiac myocytes extraction-contraction coupling, and Ca2+ transients. Incorporation of mechanical power generation would allow for the evaluation of mechanical efficiency of cardiac muscle during increased energy expenditure. Although the model has included a multistep activation mechanism, the activation was implemented in the model by increasing the rate coefficients and not through a direct link to the dynamics of Ca2+ in the cytosol and mitochondria in response to increased workload. Finally, the experimental validation data does not provide the dynamic information during the transition to increased energy expenditure. Model simulation shows that at the very beginning of increased workload, the cytosolic and mitochondrial NADH and NAD+ are slightly changed, which may be caused by the activation of glycogen breakdown and the increased NADH generation from pyruvate and fatty acid oxidation. Additional experimental studies are needed to validate these transient dynamics. Finally, a systematic sensitivity analysis is needed to quantitatively analyze the effect of each reaction or transport process on the dynamics of metabolites in response to metabolic stress.
In summary, our in silico studies suggest that with the onset of an abrupt increase in cardiac energy expenditure there is a parallel activation of glycolysis, pyruvate oxidation, fatty acid oxidation, and ATP generation, resulting in constant concentrations of ATP and ADP in the cytosol and mitochondria. On the other hand, we observed that both cytosolic and mitochondrial NADH/NAD+ increased during the first minutes (by 40% and 20%, respectively) and returned to the resting values by 15 min. Furthermore, model simulations showed that an elevation in arterial lactate concentration greatly increased the cytosolic NADH/NAD+ ratio but had a lesser effect on the mitochondrial NADH/NAD+ ratio. On the other hand, simulation of diabetic conditions reduced pyruvate oxidation and the cytosolic NADH/NAD+ ratio but did not effect mitochondrial NADH/NAD+ ratio, illustrating that alterations in mitochondrial substrate selection can occur independent of large changes in the mitochondrial NADH/NAD+ ratio.
 |
ACKNOWLEDGMENTS
|
|---|
This work was support by National Institutes of Health (NIH) Grant GM-66309 that established the Center for Modeling Integrated Metabolic Systems (MIMS) and by NIH Grant HL-074237. N. Sharma was supported by a predoctoral fellowship from the American Heart Association, Ohio Valley Affiliate.
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: W. C. Stanley, Dept. of Physiology and Biophysics, School of Medicine, Case Western Reserve Univ., 10900 Euclid Ave., Cleveland, OH 44106-4970 (e-mail: wcs4{at}case.edu)
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
 |
REFERENCES
|
|---|
- Avogaro A, Nosadini R, Doria A, Fioretto P, Velussi M, Vigorito C, Sacca L, Toffolo G, Cobelli C, and Trevisan R. Myocardial metabolism in insulin-deficient diabetic humans without coronary artery disease. Am J Physiol Endocrinol Metab 258: E606E618, 1990.[Abstract/Free Full Text]
- Balaban RS. Regulation of oxidative phosphorylation in the mammalian cell. Am J Physiol Cell Physiol 258: C377C389, 1990.[Abstract/Free Full Text]
- Balaban RS, Bose S, French SA, and Territo PR. Role of calcium in metabolic signaling between cardiac sarcoplasmic reticulum and mitochondria in vitro. Am J Physiol Cell Physiol 284: C285C293, 2003.[Abstract/Free Full Text]
- Balaban RS, Kantor HL, Katz LA, and Briggs RW. Relation between work and phosphate metabolite in the in vivo paced mammalian heart. Science 232: 11211123, 1986.[Abstract/Free Full Text]
- Beard DA. A biophysical model of the mitochondrial respiratory system and oxidative phosphorylation. PLoS Comput Biol 1: e36, 2005.
- Bose S, French S, Evans FJ, Joubert F, and Balaban RS. Metabolic network control of oxidative phosphorylation: multiple roles of inorganic phosphate. J Biol Chem 278: 3915539165, 2003.[Abstract/Free Full Text]
- Bunger R, Permanetter B, and Yaffe S. Energy utilization and pyruvate as determinants of pyruvate dehydrogenase in norepinephrine-stimulated heart. Pflügers Arch 397: 214219, 1983.[CrossRef][ISI][Medline]
- Cabrera ME, Zhou L, Stanley WC, and Saidel GM. Regulation of cardiac energetics: role of redox state and cellular compartmentation during ischemia. Ann NY Acad Sci 1047: 259270, 2005.[CrossRef][ISI][Medline]
- Chance B and Williams GR. The respiratory chain and oxidative phosphorylation. Adv Enzymol Relat Subj Biochem 17: 65134, 1956.[ISI][Medline]
- Chatham JC, Gao ZP, Bonen A, and Forder JR. Preferential inhibition of lactate oxidation relative to glucose oxidation in the rat heart following diabetes. Cardiovasc Res 43: 96106, 1999.[Abstract/Free Full Text]
- Chen R, Fearnley IM, Peak-Chew SY, and Walker JE. The phosphorylation of subunits of complex I from bovine heart mitochondria. J Biol Chem 279: 2603626045, 2004.[Abstract/Free Full Text]
- Collins-Nakai RL, Noseworthy D, and Lopaschuk GD. Epinephrine increases ATP production in hearts by preferentially increasing glucose metabolism. Am J Physiol Heart Circ Physiol 267: H1862H1871, 1994.[Abstract/Free Full Text]
- Cortassa S, Aon MA, Marban E, Winslow RL, and O'Rourke B. An integrated model of cardiac mitochondrial energy metabolism and calcium dynamics. Biophys J 84: 27342755, 2003.[Abstract/Free Full Text]
- Doria A, Nosadini R, Avogaro A, Fioretto P, and Crepaldi G. Myocardial metabolism in type 1 diabetic patients without coronary artery disease. Diabet Med 8: S104S107, 1991.
- Forsey RG, Reid K, and Brosnan JT. Competition between fatty acids and carbohydrate or ketone bodies as metabolic fuels for the isolated perfused heart. Can J Physiol Pharmacol 65: 401406, 1987.[ISI][Medline]
- Gertz EW, Wisneski JA, Stanley WC, and Neese RA. Myocardial substrate utilization during exercise in humans dual carbon-labeled carbohydrate isotope experiments. J Clin Invest 82: 20172025, 1988.[ISI][Medline]
- Goodwin GW, Taylor CS, and Taegtmeyer H. Regulation of energy metabolism of the heart during acute increase in heart work. J Biol Chem 273: 2953029539, 1998.[Abstract/Free Full Text]
- Hall JL, Lopaschuk GD, Barr A, Bringas J, Pizzurro RD, and Stanley WC. Increased cardiac fatty acid uptake with dobutamine infusion in swine is accompanied by a decrease in maloanyl CoA levels. Circ Res 32: 879885, 1996.
- Hall JL, Stanley WC, Lopaschuk GD, Wisneski JA, Pizzurro RD, Hamilton CD, and McCormack JG. Impaired pyruvate osidation but normal glucose uptake in diabetic pig heart during dobutamine induced work. Am J Physiol Heart Circ Physiol 271: H2320H2329, 1996.[Abstract/Free Full Text]
- Hall JL, Van Wylen DG, Pizzurro RD, Hamilton CD, Reiling CM, and Stanley WC. Myocardial interstitial purine metabolites and lactate with increased work in swine. Cardiovasc Res 30: 351356, 1995.[CrossRef][ISI][Medline]
- Heineman FW and Balaban RS. Control of mitochondrial respiration in the heart in vivo. Annu Rev Physiol 52: 523542, 1990.[CrossRef][ISI][Medline]
- Heineman FW and Balaban RS. Effects of afterload and heart rate on NAD(P)H redox state in the isolated rabbit heart. Am J Physiol Heart Circ Physiol 264: H433H440, 1993.[Abstract/Free Full Text]
- Heiss HW, Barmeyer J, Wink K, Hell G, Cerny FJ, Keul J, and Reindell H. Studies on the regulation of myocardial blood flow in man. I: Training effects on blood flow and metabolism of the healthy heart at rest and during standardized heavy exercise. Basic Res Cardiol 71: 658675, 1976.[CrossRef][ISI][Medline]
- Kerbey AL, Radcliffe PM, and Randle PJ. Diabetes and the control of pyruvate dehydrogenase in rat heart mitochondria by concentration ratios of adenosine triphosphate/adenosine diphosphate, of reduced/oxidized nicotinamide-adenine dinucleotide and of acetyl-coenzyme A/coenzyme A. Biochem J 164: 509519, 1977.[ISI][Medline]
- Khouri EM, Gregg DE, and Rayford CR. Effect of exercise on cardiac output, left coronary flow and myocardial metabolism in the unanesthetized dog. Circ Res 17: 427437, 1965.[Abstract/Free Full Text]
- King KL, Okere IC, Sharma N, Dyck JR, Reszko AE, McElfresh TA, Kerner J, Chandler MP, Lopaschuk GD, and Stanley WC. Regulation of cardiac malonyl-CoA content and fatty acid oxidation during increased cardiac power. Am J Physiol Heart Circ Physiol 289: H1033H1037, 2005.[Abstract/Free Full Text]
- Korzeniewski B. Regulation of ATP supply in mammalian skeletal muscle during resting stateintensive work transition. Biophys Chem 83: 1934, 2000.[CrossRef][ISI][Medline]
- Korzeniewski B. Regulation of oxidative phosphorylation in different muscles and various experimental conditions. Biochem J 375: 799804, 2003.[CrossRef][ISI][Medline]
- Korzeniewski B, Noma A, and Matsuoka S. Regulation of oxidative phosphorylation in intact mammalian heart in vivo. Biophys Chem 116: 145157, 2005.[CrossRef][ISI][Medline]
- Massie BM, Schwartz GG, Garcia J, Wisneski JA, Weiner MW, and Owens T. Myocardial metabolism during increased work states in the porcine left ventricle in vivo. Circ Res 74: 6473, 1994.[Abstract/Free Full Text]
- McCormack JG, Halestrap AP, and Denton RM. Role of calcium ions in regulation of mammalian intramitochondrial metabolism. Physiol Rev 70: 391425, 1990.[Free Full Text]
- Newsholme EA and Randle PJ. Regulation of glucose uptake by muscle. Effects of fatty acids, ketone bodies and pyruvate, and of alloxan-diabetes, starvation, hypophysectomy and adrenalectomy, on the concentrations of hexose phosphates, nucleotides and inorganic phosphate in perfused rat heart. Biochem J 93: 641651, 1964.[ISI][Medline]
- Randle PJ. Fuel selection in animals. Biochim Soc Trans 14: 799806, 1986.
- Scholz TD and Balaban RS. Mitochondrial F1-ATPase activity of canine myocardium: effects of hypoxia and stimulation. Am J Physiol Heart Circ Physiol 266: H2396H2403, 1994.[Abstract/Free Full Text]
- Sharma N and Stanley WC. Power and control: transferring energy when there is work to be done. Physiology 59: 2930, 2005.[Medline]
- Sharma N, Okere IC, Brunengraber DZ, McElfresh TA, King KL, Sterk JP, Huang H, Chandler MP, and Stanley WC. Regulation of pyruvate dehydrogenase activity and citric acid cycle intermediates during high cardiac power generation. J Physiol 562: 593603, 2005.[Abstract/Free Full Text]
- Spitzer JJ. Effect of lactate infusion on canine myocardial free fatty acid metabolism in vivo. Am J Physiol 226: 213217, 1974.[Free Full Text]
- Stanley WC. Myocardial lactate metabolism during exercise. Med Sci Sports Exerc 23: 920924, 1991.[ISI][Medline]
- Stanley WC, Gertz EW, Wisneski JA, Morris DL, Neese RA, and Brooks GA. Systemic lactate kinetics during graded exercise in man. Am J Physiol Endocrinol Metab 249: E595E602, 1985.[Abstract/Free Full Text]
- Stanley WC, Lopaschuk GD, and McCormack JG. Regulation of energy substrate metabolism in the diabetic heart. Cardiovasc Res 34: 2533, 1997.[Free Full Text]
- Stanley WC, Meadows SR, Kivilo KM, Roth BA, and Lopaschuk GD.
-Hydroxybutyrate inhibits myocardial fatty acid oxidation in vivo independent of changes in malonyl-CoA content. Am J Physiol Heart Circ Physiol 285: H1626H1631, 2003.[Abstract/Free Full Text] - Stanley WC, Recchia FA, and Lopaschuk GD. Myocardial substrate metabolism in the normal and failing heart. Physiol Rev 85: 10931129, 2005.[Abstract/Free Full Text]
- Territo PR, French SA, and Balaban RS. Simulation of cardiac work transitions, in vitro: effects of simultaneous Ca2+ and ATPase additions on isolated porcine heart mitochondria. Cell Calcium 30: 1927, 2001.[CrossRef][ISI][Medline]
- Territo PR, Mootha VK, French SA, and Balaban RS. Ca2+ activation of heart mitochondrial oxidative phosphorylation: role of the F(0)/F(1)-ATPase. Am J Physiol Cell Physiol 278: C423C435, 2000.[Abstract/Free Full Text]
- Wan B, Doumen C, Duszynski J, Salama G, Vary TC, and LaNoue KF. Effects of cardiac work on electrical potential gradient across mitochondrial membrane in perfused rat hearts. Am J Physiol Heart Circ Physiol 265: H453H460, 1993.[Abstract/Free Full Text]
- Zhou L, Salem JE, Saidel GM, Stanley WC, and Cabrera ME. Mechanistic model of cardiac energy metabolism predicts localization of glycolysis to cytosolic subdomain during ischemia. Am J Physiol Heart Circ Physiol 288: H2400H2411, 2005.[Abstract/Free Full Text]
- Zhou L, Stanley WC, Saidel GM, Yu X, and Cabrera ME. Regulation of lactate production at the onset of ischaemia is independent of mitochondrial NADH/NAD+: insights from in silico studies. J Physiol 569: 925937, 2005.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
Y. Li, R. K. Dash, J. Kim, G. M. Saidel, and M. E. Cabrera
Role of NADH/NAD+ transport activity and glycogen store on skeletal muscle energy metabolism during exercise: in silico studies
Am J Physiol Cell Physiol,
January 1, 2009;
296(1):
C25 - C46.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
L. Zhou, M. E. Cabrera, H. Huang, C. L. Yuan, D. K. Monika, N. Sharma, F. Bian, and W. C. Stanley
Parallel activation of mitochondrial oxidative metabolism with increased cardiac energy expenditure is not dependent on fatty acid oxidation in pigs
J. Physiol.,
March 15, 2007;
579(3):
811 - 821.
[Abstract]
[Full Text]
[PDF]
|
 |
|