AJP - Heart BIOPAC complete lab solutions
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 291: H1466-H1474, 2006. First published May 5, 2006; doi:10.1152/ajpheart.00004.2006
0363-6135/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
291/3/H1466    most recent
00004.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Korzeniewski, B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Korzeniewski, B.

Oxygen consumption and metabolite concentrations during transitions between different work intensities in heart

Bernard Korzeniewski

Faculty of Biotechnology, Jagiellonian University, Kraków, Poland

Submitted 3 January 2006 ; accepted in final form 20 April 2006


    ABSTRACT
 TOP
 ABSTRACT
 THEORETICAL PROCEDURES
 THEORETICAL RESULTS
 DISCUSSION
 REFERENCES
 
Steady-state metabolite (ADP, ATP, Pi, PCr, and NADH) concentrations usually differ little between different workloads with significantly different oxygen consumption rates in the heart. However, during transitions between steady states, metabolite concentrations may in some cases change transiently, exhibiting a significant overshoot or undershoot, whereas in other cases they approach near-exponentially new steady-state values. Oxygen consumption rate usually reaches the new steady-state value very quickly (within a few seconds). The present in silico studies, performed using a previously developed computer model of oxidative phosphorylation in the heart, demonstrate that such a behavior of the oxidative phosphorylation system can be reproduced only under the assumption that ATP usage, substrate dehydrogenation, and (particular steps of) oxidative phosphorylation are directly activated to a similar extend by some cytosolic factor/mechanism during transition from low work to high work (the so-called parallel-activation mechanism). Computer simulations show that some differences observed between different experimental systems can be explained by a slightly different balance of the activation of particular components of the system and/or by a delay in time of the activation/inactivation of substrate dehydrogenation and oxidative phosphorylation during low-to-high and high-to-low work transitions. Thus the presented theoretical approach offers a general idea that is able to unify, at least semiquantitatively, different experimental data available in the literature.

parallel activation; oxidative phosphorylation; transition time; computer model


OXIDATIVE PHOSPHORYLATION in the mitochondria constitutes the main source of energy in the form of ATP in the heart. In cardiac myocytes, ATP is utilized by actomyosin-ATPase, sarcoplasmic reticulum Ca2+-ATPase, and sarcolemmal Na+-K+-ATPase and by basic processes keeping the cell alive. The rate of respiration and ATP turnover in the heart can increase several times during transition from a low beating frequency (workload) to a high beating frequency (workload) and decrease to its initial value during the opposite transition. In the intact heart in vivo, essentially no differences in oxidative phosphorylation metabolite [ADP, ATP, Pi, phosphocreatine (PCr), and NADH] concentrations between different steady states with different work intensities are observed (1, 10, 12, 26). In the perfused heart, [ADP] and [Pi] may be higher, whereas [PCr] may be lower at high work than at low work (8, 29, 30). In isolated trabeculae, [NADH] may either drop or increase significantly when electrical stimulation frequency is increased (2, 3). Additionally, in this system there is an initial undershoot/overshoot of [NADH] during transitions between steady states with different work intensities set by varying electrical stimulation frequencies: a significant transient decrease/increase in [NADH] takes place immediately after an increase/decrease in stimulation frequency, respectively (2, 3). A rather similar pattern of behavior is observed if the work intensity in isolated trabeculae is varied by external calcium concentration or sarcomere length (2). A recent study (11) showed a similar time course of [NADH] in isolated cardiac myocytes. On the other hand, in paced isolated, perfused hearts, no overshoots or undershoots were observed in [PCr] and [Pi] when [Pi] increased/decreased and [PCr] decreased/increased after an increase/decrease in pacing frequency, respectively (7, 29). Finally, under physiological conditions (intact heart, 37°C, normoxia, and lack of inhibitors), the oxygen consumption (VO2) in the heart rises very quickly during transition from low work to high work: the half-transition time is usually shorter than 10 s and can be as low as 4 s (7, 29). The half-transition time for [PCr] can be even shorter –2.5 s (7, 29). In contrast, the half-transition time for VO2 and [PCr] in skeletal muscle at the onset of exercise is much greater: usually between 20 and 65 s (31).

The rate of ATP production by oxidative phosphorylation must be adjusted to the current energy demand to match exactly the rate of ATP consumption. Three main mechanisms of adjusting the rate of ATP supply by oxidative phosphorylation have been proposed in the literature. Because of the first mechanism, which can be called the output-activation mechanism, only the ATP usage block (output of the system) is directly activated by calcium ions (Ca2+) during elevated muscle work, whereas oxidative phosphorylation is activated only indirectly through the negative feedback involving an increase in [ADP] (and [Pi]) (4). The input/output activation mechanism involves the direct activation of substrate dehydrogenation (input of the system), especially the "key" tricarboxylic acid (TCA) cycle dehydrogenases: pyruvate dehydrogenase, isocitrate dehydrogenase, and 2-oxoglutarate dehydrogenase by Ca2+ in parallel with the activation of ATP usage (9, 22). Finally, because of the each-step-activation mechanism or parallel-activation mechanisms, some cytosolic factor (e.g., Ca2+) directly activates all oxidative phosphorylation enzymes, substrate dehydrogenation block, and ATP usage block (1317). In particular, it has been postulated recently (17) that only the each-step-activation mechanism is able to account for the almost perfect stability of several oxidative phosphorylation metabolites (ADP, PCr, Pi, and NADH) between steady states with very different work intensities and VO2 rates in the heart. It was also shown earlier that this mechanism is able to explain many, apparently independent, kinetic properties of oxidative phosphorylation in the skeletal muscle (1316, 18).

However, in the previous paper (17) concerning the regulation of oxidative phosphorylation in intact mammalian hearts in vivo, only the steady-state properties of the system were analyzed. Therefore, the purpose of the present article is to study in silico the kinetic behavior of oxidative phosphorylation during low-to-high and high-to-low work transitions, assuming different mechanisms of the regulation of oxidative phosphorylation and to confront the theoretical predictions with experimental data obtained in the intact heart and isolated trabeculae under various conditions.


    THEORETICAL PROCEDURES
 TOP
 ABSTRACT
 THEORETICAL PROCEDURES
 THEORETICAL RESULTS
 DISCUSSION
 REFERENCES
 
Computer model. The model of oxidative phosphorylation in the intact mammalian heart published previously (17) was used in the present in silico study. Because the VO2 and ATP turnover per unit of heart mass depend on the animal size (both being greater in smaller animals) (see Ref. 17), it is important to stress that this model was developed for middle-size mammals (dogs, pigs, and humans). The discussed model takes into account explicitly the following enzymes/processes/metabolic blocks: substrate dehydrogenation (glycolysis, TCA cycle, fatty acid beta-oxidation, etc.), complex I, complex III, complex IV (cytochrome-c oxidase), proton leak, ATP synthase, ATP/ADP carrier, phosphate carrier, adenylate kinase, creatine kinase, and ATP usage system. The time variations of the metabolite concentrations that constitute independent variables (NADH, ubiquinol, reduced form of cytochrome c, O2, internal protons, internal ATP, internal Pi, external ATP, external ADP, external Pi, external protons, and PCr) are expressed in the form of a set of ordinary differential equations. The set of differential equations is integrated numerically using the Gear procedure. The simulation program is written in FORTRAN. The complete description of the model of oxidative phosphorylation in the intact heart is given in Ref. 17 and is located on the website http://awe.mol.uj.edu.pl/~benio.

Computer simulations. In the present paper it was assumed that ATP usage at high work is five times greater than at low work (12, 17). The VO2 rate equals 2.5 mM/min at low work and 11.2 mM at high work (17). The ratio 11.2/2.5 is slightly lower than 5 because some VO2 is due to the proton leak across the inner mitochondrial membrane and not the ATP turnover (within the model the proton leak accounts for 17% and 4% of VO2 at low and high work, respectively). The above values of VO2 correspond to the minimal and maximal work intensity in the paced intact dog heart in vivo measured in (12), although even higher values of the maximal VO2 equal to 12–15 mM/min were found in epinephrine- and phenylephrine-stimulated dog hearts (12). Therefore, it can be estimated that at least a three- to sixfold maximal increase in VO2 during low-to-high work transition takes place. However, the present article refers to different experiments performed in different systems (heart in vivo, perfused heart, isolated trabeculae) on hearts of different mammals stimulated in different ways (pacing, isoprenaline administration, adrenaline administration, Ca2+ administration, and increase in sarcomere length). Therefore, the theoretical predictions generated by computer simulations may be in many cases regarded as only semiquantitative ones. Nevertheless, they are useful even when they are able to explain only a given type of behavior of the system.

The fivefold (or, in some cases, twofold) activation/inactivation of ATP usage in computer simulations during low-to-high/high-to-low work transition was equivalent to a fivefold (twofold) increase/decrease in the rate constant of ATP usage. Also in some cases, other components of the system (substrate dehydrogenation and/or particular steps of oxidative phosphorylation) were directly activated/inactivated n times, which corresponded to an n-fold increase/decrease in the appropriate rate constants (or maximal velocities). It was proposed previously that in skeletal muscle the activation/inactivation of the ATP-producing system during rest-to-work/work-to-rest transition is not instant but occurs with some delay (16, 20). An analogous delay, caused by calcium uptake by mitochondria, was proposed for the activation of TCA cycle dehydrogenases in isolated trabeculae (2, 5). Therefore, in some simulations in the present study, the activation/inactivation of substrate dehydrogenation and oxidative phosphorylation during low-to-high/high-to-low work transition was not instant but delayed in time. This delay was expressed by the following equation during low-to-high work transient (16, 20):

Formula 1(1)
and high-to-low work transient (16, 20):

Formula 2(2)
where m is the current (at time t) relative activity of a given step expressed as a multiplicity of the low work activity, t is the time after the onset of on- or off-transition, n is the relative activity of this step at the high-work steady state expressed as a multiplicity of the low work activity, whereas {tau}(ON) and {tau}(OFF) are characteristic delay times of activation/inactivation of this step during low-to-high/high-to-low work transition, respectively. Therefore, according to Eqs. 1 and 2, the degree of activation m increases exponentially from 1 to n or decreases exponentially from n to 1 during low-to-high and high-to-low work transition, respectively.

For the purposes of the present in silico studies, the oxidative phosphorylation system was divided into five blocks (subsystems) that are presented in Fig. 1. All components of a given block were always (in)activated to the same extent in computer simulations. These blocks comprise the following: substrate dehydrogenation (DH), oxidation subsystem (OX) (complexes I, III and IV of the respiratory chain), proton leak (LK), phosphorylation subsystem (PH) (ATP synthase, ATP/ADP carrier, Pi carrier), and ATP utilization (UT). DH produces NADH that is oxidized by OX, which is coupled with building up the proton motive force {Delta}p (related to the proton gradient across the inner mitochondrial membrane). Some {Delta}p is dissipated by LK, but most of it is used for the synthesis (and transport) of ATP, which is in turn used by UT. In different simulations, different blocks were directly activated/inactivated in parallel with the activation/inactivation of ATP usage. The only exception was a proton leak that was never activated during transition to high work, because such an activation would be equivalent to a net energy waste.


Figure 1
View larger version (4K):
[in this window]
[in a new window]
 
Fig. 1. General scheme of oxidative phosphorylation in heart. Blocks distinguished in the present study are shown. DH, substrate dehydrogenation (tricarboxylic acid cycle, glycolysis, fatty acid beta-oxidation, and so on); OX, oxidative block (respiratory chain: complexes I, III, and IV); LK, proton leak; PH, phosphorylation block (ATP synthase, ATP/ADP carrier, Pi carrier); UT, ATP utilization block (actomyosin-ATPase, Ca2+-ATPase, Na+-K+-ATPase, etc.); {Delta}p, protonmotive force; X, factor(s)/mechanism(s) activating directly particular components of the system.

 
Generally, five modes of computer simulations were used in the present study. They are presented in Table 1.


View this table:
[in this window]
[in a new window]
 
Table 1. Characterization of the simulation modes used in the present study

 
Mode A (output activation) corresponds to the output-activation mechanism. Only ATP usage (UT) is (instantly) activated/inactivated here (5 times) during low-to-high/high-to-low work transition, whereas other blocks are not (in)activated [1-fold (in)activation].

Mode B (input/output activation) represents the input/output activation where ATP usage and substrate dehydrogenation are (instantly) (in)activated to the same extent (5 times), whereas the activity of the remaining blocks remains constant.

Mode C (perfectly balanced each-step activation) represents a perfectly balanced parallel activation of all blocks (each-step activation). This does not mean that all blocks are activated/inactivated to exactly the same extent during low-to-high/high-to-low work transition, but that the concentrations of all metabolites remain constant. Because the pathway is branched at the proton motive force {Delta}p (some external protons are not used by ATP synthase but by proton leak LK) (compare Fig. 1), PH and UT must be (in)activated slightly more than DH and OX to keep intermediate metabolite concentrations constant. Therefore, a compensation factor (comp) of 0.89 was introduced, and DH and OX were (in)activated in mode C 5·comp (= 4.45) times, whereas PH and UT were (in)activated five times. It must be emphasized that mode C represents an idealized case in which no oxygen diffusion limitations and an instant activation/inactivation of oxidative phosphorylation during on- and off-transient are assumed. This results in an instant increase/decrease in Formula 2O2 during low-to-high and high-to-low work transition (see below, Fig. 4). This assumption is rather unrealistic, and therefore simulations performed in mode C demonstrate only that the transitions of Formula 2O2 may be potentially very quick, but not instant, in the real system.


Figure 4
View larger version (12K):
[in this window]
[in a new window]
 
Fig. 4. Simulated low-to-high and high-to-low work transition in mode C (perfectly balanced parallel-activation mechanism). A: time courses of Formula 4O2, [ADP], and [NADH]. B: time courses of [ATP], [PCr], and [Pi]. C: time courses of relative activities of UT, PH, OX, and DH.

 
Mode D (disturbed each-step activation) is intended to represent some imbalance in perfect parallel activation because ATP production is activated slightly less than ATP consumption. Therefore, all blocks in the ATP production system (apart from proton leak) were less (in)activated by an imbalance factor imb = 0.8 than in mode C. Thus PH is (in)activated 4 times (5 x imb), whereas DH and OX are (in)activated 3.56 times (5 x comp x imb).

In mode E (disturbed and delayed each-step activation), the extent of (in)activation of particular blocks is exactly the same as in mode D, but the (in)activation of DH, OX, and PH is delayed in time. For each of these blocks, the characteristic delay times for low-to-high and high-to-low work transition {tau}(ON) and {tau}(OFF) were equal to 20 and 6 s, respectively. These values were selected to match approximately the experimental data.

In all simulations, the low-work steady-state was maintained for 1 min, energy demand (rate constant of ATP usage) was increased for 2 min to the high-work value, and then energy demand was decreased back to its low-work value.


    THEORETICAL RESULTS
 TOP
 ABSTRACT
 THEORETICAL PROCEDURES
 THEORETICAL RESULTS
 DISCUSSION
 REFERENCES
 
Figure 2 presents the computer simulation performed in mode A (output-activation mechanism). It can be seen that an activation of only ATP usage during low-to-high work transition results in two general properties of the kinetic behavior of the system. First, changes in intermediate metabolite (ADP, Pi, PCr, NADH) concentrations are huge. During transition from low work to high work, [ADP] increases many times from 32 µM to over 1,500 µM. A significant decrease in [ATP] related to a great increase in [ADP] and [AMP] takes place. [PCr] drops almost to zero, whereas [Pi] increases several times. [NADH] decreases significantly. Of course, during high-to-low work transition, opposite changes take place. Second, the changes in Formula 2O2 and metabolite concentrations are not very quick but occur with some half-transition time t1/2. For instance, this time for both Formula 2O2 and [PCr] is equal to 14 s during low-to-high work transition. As discussed above, in the simulation presented in Fig. 2, the ATP usage (energy demand) was five times higher at high work than at low work, which is quite a large difference. When a twofold difference in ATP demand was assumed in computer simulations (which resulted in an 1.8-fold difference in Formula 2O2; compare the above discussion concerning proton leak), the changes in metabolite concentrations were significantly smaller, but still large (for instance [ADP] increased twice), whereas the half-transition times for Formula 2O2 and [PCr] remained essentially the same (14 s) (computer simulations not shown). Therefore, the transition time for Formula 2O2 seems to be essentially independent on the work intensity in the stimulated heart.


Figure 2
View larger version (14K):
[in this window]
[in a new window]
 
Fig. 2. Simulated low-to-high and high-to-low work transition in mode A (output-activation mechanism). A: time courses of oxygen consumption (Formula 4O2), [ADP], and [NADH]. B: time courses of [ATP], phosphocreatine ([PCr]), and [Pi]. C: time courses of relative activities of UT, PH, OX, and DH.

 
The simulation performed in mode B (input/output activation), presented in Fig. 3, is generally similar to that performed in mode A, although changes in metabolite concentrations are slightly smaller. The main difference is that [NADH] slightly increases during low-to-high work transition, exhibiting however a large overshoot, whereas during the opposite transition, it slightly decreases with a large undershoot. Depending on the proportional extend of the direct activation/inactivation of ATP production and ATP consumption, the steady-state [NADH] may decrease/increase more than it is shown in Fig. 3 when ATP production is activated/inactivated less than ATP consumption, or it can even increase/decrease when ATP production is activated/inactivated more than ATP consumption during low-to-high/high-to-low work transition, respectively. The half-transition times for Formula 2O2 and [PCr] are slightly shorter than in mode A and equal 11 s.


Figure 3
View larger version (13K):
[in this window]
[in a new window]
 
Fig. 3. Simulated low-to-high and high-to-low work transition in mode B (input/output-activation mechanism). A: time courses of Formula 4O2, [ADP], and [NADH]. B: time courses of [ATP], [PCr], and [Pi]. C: time courses of relative activities of UT, PH, OX, and DH.

 
A completely different behavior of the system can be observed in mode C (perfectly balanced parallel activation). In this case, as it can be seen in Fig. 4, no changes in metabolite concentrations can be observed during low-to-high and high-to-low work transitions, whereas Formula 2O2 increases or decreases instantly after transition. However, as it is discussed throughout this study, the last property is most probably some idealization of the real system where, e.g., some oxygen diffusion limitations are likely to take place.

When the direct activation of the ATP-producing system was slightly decreased in relation to mode C in the simulation performed in mode D (disturbed parallel activation), some very moderate changes in metabolite concentrations took place, whereas changes in Formula 2O2 were very quick, although not instant, during low-to-high and high-to-low work transition (Fig. 5). Of course, when the imbalance of parallel activation was increased in mode D (imbalance factor imb was <0.8), the differences between the behavior of the system in mode C and mode D were greater (theoretical results not shown).


Figure 5
View larger version (12K):
[in this window]
[in a new window]
 
Fig. 5. Simulated low-to-high and high-to-low work transition in mode D (disturbed parallel-activation mechanism). A: time courses of Formula 4O2, [ADP], and [NADH]. B: time courses of [ATP], [PCr], and [Pi]. C: time courses of relative activities of UT, PH, OX, and DH.

 
Figure 6 presents the computer simulation performed in mode E (disturbed parallel activation with time delay). Here the delay in the activation/inactivation of substrate dehydrogenation and oxidative phosphorylation (OX + PH) resulted in significant overshoots/undershoots in metabolite concentrations. [ADP] and [Pi] experienced overshoot and [NADH] and [PCr] experienced undershoot during low-to-high work transition, whereas the opposite behavior of the system took place during high-to-low work transition.


Figure 6
View larger version (13K):
[in this window]
[in a new window]
 
Fig. 6. Simulated low-to-high and high-to-low work transition in mode E (disturbed and delayed parallel-activation mechanism). A: time courses of Formula 4O2, [ADP], and [NADH]. B: time courses of [ATP], [PCr], and [Pi]. C: time courses of relative activities of UT, PH, OX, and DH.

 

    DISCUSSION
 TOP
 ABSTRACT
 THEORETICAL PROCEDURES
 THEORETICAL RESULTS
 DISCUSSION
 REFERENCES
 
The present in silico study shows that the each-step-activation mechanism, postulated previously (13, 16, 17) according to which all oxidative phosphorylation complexes are directly activated in parallel with ATP usage during elevated work intensity, is able to explain the regulation of oxidative phosphorylation during low-to-high and high-to-low work transition in the heart. In particular, it offers a general theory that can unify, at least semiquantitatively, the kinetic behavior of the system during work transitions in the intact heart in vivo, in the perfused heart, and in isolated trabeculae.

Computer simulations show that both the output-activation mechanism (mode A, Fig. 2) and input/output activation mechanism (mode B, Fig. 3) generate large changes in metabolite (ADP, PCr, Pi, and NADH) concentrations and half-transition times t1/2 for Formula 2O2 and [PCr] longer than 10 s. On the other hand, intermediate metabolite concentrations are essentially constant in intact mammalian heart in vivo (1, 10, 12, 26), whereas in perfused hearts and isolated trabeculae relatively moderate changes in metabolite concentrations take place (2, 8, 29, 30). Additionally, the half-transition time for Formula 2O2 is usually very short in the intact heart. In most cases it equals 4–8 s [and sometimes 11 s (7)] under "normal" physiological conditions (37°C, normoxia, lack of inhibitors) (see Table 2 in Ref. 29), although the transition time may be longer at 20°C or 28°C, in hypoxia or in the presence of inhibitors (7, 29). The last finding is logical because low temperature slows down all biological processes, whereas hypoxia and inhibitors lower the activity of oxidative phosphorylation. Therefore, the computer simulations performed in the present study are relevant to the heart working under "physiological" conditions. In addition, the t1/2 for [PCr] may be even much shorter (about 2.5 s) (29). Therefore, these experimental results cannot be matched by computer simulations performed in modes A and B.

The perfectly balanced parallel-activation (each-step-activation) mechanism (mode C, Fig. 4) predicts no changes in [ADP], [PCr], [Pi], and [NADH]. Therefore, as it was shown previously (17), it is able to explain the perfect stability of these metabolite concentrations at different work intensities encountered in the intact heart in vivo (1, 10, 12, 26). This mechanism also implies that, in the absence of any delay in the activation/inactivation of ATP production, Formula 2O2 may increase/decrease very quickly (in an idealized case instantly) to its new steady-state value during low-to-high and high-to-low work transition (Fig. 4A), and thus t1/2 for Formula 2O2 may be very close to zero. However, most probably Formula 2O2 does not increase/decrease instantly in the real system because of oxygen diffusion limitations. In fact, it was postulated (29) that diffusion limitations are responsible for the significantly shorter t1/2 for [PCr] than for Formula 2O2. Unfortunately, the time course of Formula 2O2 has not been measured in the intact heart in vivo. All existing experimental data, especially those discussed in Refs. 7 and 29, concern perfused heart where some changes in metabolite concentrations take place and therefore correspond to mode D. However, the parallel-activation mechanism offers, at least potentially, a possibility of a very quick Formula 2O2 adjustment during transitions between different steady states with different workloads.

The disturbed parallel-activation mechanism (mode D, Fig. 5) predicts moderate changes in intermediate metabolite concentrations and a short (a few seconds) half-transition time for Formula 2O2 and [PCr]. Such a kinetic behavior of the system reflects well, at least semiquantitatively, the situation prevailing in many cases in the perfused heart (7, 8, 29, 30). The extent of the changes of metabolite concentrations between low work and high work may depend significantly on the composition of the perfusion medium, for instance, on the respiratory substrates and hormones present in it (8). Therefore, these factors probably determine, at least partly, the size of the imbalance of the parallel activation in the perfused heart.

Finally, computer simulations involving the disturbed parallel-activation mechanism with delayed activation/inactivation of ATP production (mode E, Fig. 6) predict significant overshoots/undershoots in [NADH], [ADP], [PCr], and [Pi]. A very similar [NADH] undershoot during low-to-high work transition and [NADH] overshoot during high-to-low work transition are observed in isolated trabeculae (2) and isolated heart cells (11). Therefore, in this case it seems likely that the activation/inactivation of ATP production is significantly delayed in time. This delay may be due to the transport of calcium through the inner mitochondrial membrane (2, 5) and/or to the kinetics of the activation/inactivation of TCA cycle and oxidative phosphorylation enzymes (16, 20). It is important to emphasize that in isolated trabeculae the steady-state [NADH] differs much more between a steady state with very low, unphysiological stimulation frequency (e.g., 0.25 Hz) and a steady state with a higher, more physiological stimulation frequency (e.g., 1 Hz) than between two steady states with higher stimulation frequencies (e.g., 1 and 3 Hz) (2), as in the simulation presented in Fig. 6. In fact, when the highest stimulation frequencies (4–6 Hz) are applied, NADH level even increases (2). Such a behavior of the system can be reproduced by the model if different degrees of activation of the NADH-producing and NADH-consuming subsystem are assumed (not shown). Unfortunately, changes over time in [ADP], [PCr], and [Pi] during transitions between different stimulation frequencies have not been measured in the discussed experimental system (isolated trabeculae). Significant overshoots/undershoots in these metabolite concentrations, similar to those presented in Fig. 6, were observed in a isoprenaline-stimulated heart (28). However, in this case, they were probably due to a transiently increased work output. The mode E of computer simulations predicts a small overshoot in Formula 2O2 during low-to-high work transition (see Fig. 6A). However, again, the time course of mitochondrial Formula 2O2 was not measured in isolated trabeculae, and therefore it is not possible to test this theoretical prediction [it also depends on the characteristic delay time {tau}(ON), not shown].

Some experimental results seem to represent intermediate cases between some of the modes defined in the present article. For instance, in electrically paced rat heart cells at 25°C, there is a small [NADH] undershoot during low-to-high work transition and a small overshoot during the opposite transition (Fig. 2 in Ref. 32). Such a behavior is intermediate between modes D and E. In the presence of deoxyglucose there is even some overshoot of [NADH] during low-to-high work transition (Fig. 2 in Ref. 32), whereas steady-state [NADH] remains essentially constant, which suggests a case intermediate between modes B and C (but closer to mode C).

Of course, the five modes of simulations defined in the present study do not exhaust all potential possibilities and combinations of the kinetic properties of the system. For instance, one can imagine mode C (perfectly balanced parallel activation) with some time delay of the activation/inactivation of different steps of ATP production. Additionally, as discussed above, mode C constitutes a certain approximation of the reality and most probably should be supplemented with oxygen diffusion limitations. Maybe the future experimental studies will cause the need of introducing additional modes. However, essentially all experimental data known presently can be reproduced, at least semiquantitatively, within the frame of modes C, D, and E.

Generally, different versions of the parallel-activation mechanism (perfectly balanced, disturbed, disturbed, and delayed) seem to be able to explain the changes in Formula 2O2 and intermediate metabolite concentrations observed during transitions among different work intensities in the intact heart, perfused hearts, and isolated trabeculae studied in different conditions. Therefore, the idea of parallel activation supplemented with some modifications (some imbalance and/or delay in the activation/inactivation of ATP production) offers a unified theory of the regulation of oxidative phosphorylation during transitions between different work intensities in the heart. In particular, this idea is necessary to explain the short transition times for Formula 2O2 and [PCr] and very moderate (if any) changes in metabolite concentrations accompanying substantial changes in Formula 2O2 and work intensity. It seems that, roughly, mode C (supplemented with oxygen diffusion limitations) represents the intact heart in vivo, mode D represents the perfused heart, and mode E represents isolated trabeculae. Therefore, it can be said that the further from the physiological conditions prevailing in vivo, the more disturbed/delayed the parallel activation is.

The half-transition time t1/2 is proportional to the difference in metabolite concentrations between different steady states and inversely proportional to the metabolite flux through the system. Easterby (6) defined the transition time {tau}A->B between the initial steady-state A with flux JA and final steady-state B with flux JB in the following way:

Formula 3(3)
where [Ij]ss is the steady-state concentration of the intermediate metabolite j. However, as it was discussed previously (20), the involvement of only the final flux JB in the above expression causes a great asymmetry between the transition times during on-transient (rest-to-work transition) and off-transient (work-to-rest transition) in skeletal muscle. Therefore, the following expression seems to be more appropriate:

Formula 4(4)
where (JB + JA)/2 is the average transition flux. This equation helps to understand why the transition times for Formula 4O2 and [PCr] are generally greater in skeletal muscle [about 20–65 s (31)] than in the heart working under physiological conditions (37°C, normoxia, lack of inhibitors) (2.5–11 s, see Refs. 7 and 29). Namely, changes in intermediate metabolite concentrations during transitions between different work intensities in the heart are smaller than in skeletal muscle, whereas fluxes are higher. Additionally, the computer simulation presented in Fig. 6 (mode E) demonstrates that the transition time depends not only on the changes in metabolite concentrations, but also on the time course of the activation/inactivation of particular components of a considered system. As it was discussed before, the half-transition time t1/2 is not identical with {tau}A->B. However, {tau}A->B defined in Eq. 4 is proportional to t1/2.

Cortassa and coworkers (5) developed a computer model of oxidative phosphorylation and TCA cycle involving the calcium-activation of three "key" TCA cycle dehydrogenases (pyruvate dehydrogenase, isocitrate dehydrogenase, and 2-oxoglutarate dehydrogenase). They used their model to study in silico the time course of [NADH] in isolated trabeculae during low-to-high and high-to-low work transition (caused by varying electrical stimulation frequency) (2). They were able to reproduce, at least semiquantitatively, the undershoot in [NADH] during the low-to-high work transition and the overshoot in [NADH] during the opposite transition (5). However, an eightfold increase in stimulation frequency, and thus in ATP usage for contraction and calcium transport, was accompanied in this model by only a very small relative increase in the respiration rate. Formula 4O2 increased by 20% if it was assumed that the TCA cycle keeps most of the control over the Formula 4O2 flux. However, it has been demonstrated experimentally that the respiration rate is mostly controlled (under conditions where ATP usage has no control) by oxidative phosphorylation (21), whereas the flux control coefficient for TCA cycle equals about 0.1 (10% of metabolic control). Under this more realistic assumption that the control over Formula 4O2 is exerted predominantly by oxidative phosphorylation, the respiration rate in the model developed by Cortassa and coworkers increases by only 10% (5). This example shows very clearly that, first, many different models can be adjusted by manipulating with assumptions and free parameter values to some simple set of experimental data, and therefore, second, any model should be thoroughly validated by comparison of its predictions with possibly a great number of different parameter values and system properties. The computer simulation shown in Fig. 6 demonstrates that both the undershoot/overshoot in [NADH] and a substantial increase/decrease in Formula 4O2 during low-to-high/high-to-low work transition can be predicted by a computer model when the delayed-and-disturbed parallel-activation mechanism (mode E) and not the delayed input/output-activation mechanism used by Cortassa and coworkers is applied (in fact, in the model by Cortassa and coworkers there is no ATP usage but only a pulse of ADP introduced ad hoc to force the expected behavior of the system).

The parallel activation of different steps of the ATP-producing system seems to develop gradually as an animal matures, because in young sheep, changes in metabolite concentrations between different work intensities are significantly greater than in adult sheep (24). It was shown earlier that the each-step-activation mechanism is able to explain many, apparently unrelated to each other, kinetic properties of oxidative phosphorylation in the skeletal muscle (13–16, 18). These properties comprise 1) the relative stability of [ADP], [Pi], PCr/Cr, NADH/NAD+, and {Delta}p during rest-to-work transition (13, 16); 2) the much steeper phenomenological Formula 4O2/[ADP] relationship in the intact skeletal muscle than in isolated muscle mitochondria (13, 16); 3) the two to four (to six) times greater maximum oxygen consumption in intact skeletal muscle than in isolated mitochondria, skinned fibers, or muscle homogenate when scaled for the amount of mitochondrial proteins (18, 25, 27); 4) the increase in the relative slope of the phenomenological Formula 4O2/[ADP] relationship as a result of muscle training (19); 5) the greater Formula 4O2 at a given [ADP] and amount of mitochondrial protein in trained muscle than in untrained muscle and in untrained muscle than in hypothyroid muscle (16); 6) the asymmetry of the half-transition time t1/2 for [PCr] between the on-transient (rest-to-work transition) and off-transient (work-to-rest transition) (16); and 7) the PCr recovery overshoot (16) and the variability of the kinetic properties of oxidative phosphorylation in different muscles and various experimental conditions (16). Therefore, the idea of the each-step activation in both heart and skeletal muscle seems to be well founded. The differences in the regulation of oxidative phosphorylation between these tissues are discussed in Ref. 17.

The parallel-activation mechanism implies that the potential capacity of oxidative phosphorylation for ATP synthesis in the intact heart is significantly higher than the maximal ATP synthesis rate by isolated mitochondria (when scaled for the amount of mitochondrial proteins). In intact skeletal muscle the maximal Formula 4O2 is two to six times greater than in isolated mitochondria (18, 25, 27). On the other hand, Mootha and coworkers (23) reported that the hearts of the exercising dog and pig use only up to 80–90% of the maximal oxidative phosphorylation capacity in isolated heart mitochondria. However, the maximal Formula 4O2 in mitochondria was measured at saturating concentrations of ADP, Pi, oxygen, and respiratory substrates. On the other hand, [ADP] in intact heart in vivo is around the Km value of oxidative phosphorylation for ADP and is essentially constant even during transition to very high work intensities (12). Also [Pi] is constant and not saturating (12). Furthermore, it is very likely that at high work intensities in the intact heart, oxygen and respiratory substrate supply by blood becomes limiting for Formula 4O2 and ATP synthesis by oxidative phosphorylation. Taken together, the above properties of the system in vivo suggest that the "true" (potential) capacity of oxidative phosphorylation in the intact heart is at least twice (and perhaps much more) greater than in isolated mitochondria, in accordance with the idea of parallel activation.

The physical nature of the factor(s)/mechanism(s) X responsible for the direct parallel activation of (particular components of) the oxidative phosphorylation system during low-to-high work transition remains unclear. Several evidences suggest some involvement of Ca2+. However, it is not certain whether a high calcium concentration is the (only) factor looked for. It has been proposed (14) that the relevant factor may be the frequency of calcium oscillations. This frequency could be integrated over time by some protein, analogous to calmodulin, that is lost or inactivated during isolation of mitochondria. On the other hand, the fact that ruthenium red, an inhibitor of Ca2+ uptake by mitochondria, only slightly lengthens the low-to-high work transition but does not affect much the steady-state metabolite concentrations (28) suggests that either this inhibitor only delays but does not block calcium entry into mitochondria, or there exists another factor/mechanism that acts in parallel with Ca2+. There is also a possibility that, whereas TCA cycle dehydrogenases are activated by matrix calcium, the oxidative phosphorylation complexes, being transmembrane proteins located in the inner mitochondrial membrane, are activated by cytosolic calcium and therefore do not need the Ca2+ transport to the matrix for their activation.

Of course, the computer model used in the present study contains several simplifications. For instance, as discussed above, it does not take into account the possible oxygen diffusion limitations that can delay Formula 4O2 transitions in the real system. Additionally, the NADH-producing system is described by a simple phenomenological equation. Theoretical studies using a model containing a more sophisticated description of this system demonstrated that rapid activation of glycolysis and cytosolic NADH formation may be important in the transition to high work states in the heart (33, 34).

In summary, the present theoretical study demonstrates that the parallel-activation mechanism of the adjustment of ATP production to the present energy demand is able to explain, at least in principle, different kinds of the dynamic behavior of the oxidative phosphorylation system (changes over time of Formula 4O2 and intermediate metabolite concentrations) encountered in different experimental conditions during low-to-high and high-to-low work transitions. Generally, under physiological conditions (intact heart in vivo) the direct activation of particular elements of the ATP-producing system (substrate dehydrogenation and oxidative phosphorylation) in parallel with the activation of ATP usage during transition from low work to high work seems to be perfectly balanced and not significantly delayed in time. In less physiological conditions (perfused heart, isolated trabeculae), this parallel activation is likely to be disturbed (ATP production is activated slightly less than ATP consumption) and/or delayed (it takes some time to activate/inactivate TCA cycle and oxidative phosphorylation). However, in all the analyzed cases a direct activation/inactivation of particular oxidative phosphorylation complexes during low-to-high/high-to-low work transition is needed to reproduce the experimental data.


    FOOTNOTES
 

Address for reprint requests and other correspondence: B. Korzeniewski, Faculty of Biotechnology, Jagiellonian Univ., ul. Gronostajowa 7, 30-387 Kraków, Poland (e-mail: benio{at}mol.uj.edu.pl)

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
 TOP
 ABSTRACT
 THEORETICAL PROCEDURES
 THEORETICAL RESULTS
 DISCUSSION
 REFERENCES
 

  1. Balaban RS, Kantor HL, Katz LA, and Briggs RW. Relation between work and phosphate metabolite in the in vivo paced mammalian heart. Science 232: 1121–1123, 1986.[Abstract/Free Full Text]
  2. Brandes R and Bers DM. Intracellular Ca2+ increases the mitochondrial NADH concentration during elevated work in intact cardiac muscle. Circ Res 80: 82–87, 1997.[Abstract/Free Full Text]
  3. Brandes R and Bers DM. Increased work in cardiac trabeculae causes decreased mitochondrial NADH fluorescence followed by slow recovery. Biophys J 71: 1024–1035, 1996.[Web of Science][Medline]
  4. Chance B and Williams GR. The respiratory chain and oxidative phosphorylation. Adv Enzymol Relat Subj Biochem 17: 65–134, 1956.[Web of Science][Medline]
  5. Cortassa S, Aon MA, Marban E, Winslow RL, and O’Rourke B. An integrated model of cardiac energy metabolism and calcium dynamics. Biophys J 84: 2734–2755, 2003.[Web of Science][Medline]
  6. Easterby JS. A general theory of the transition time for sequential enzyme reactions. Biochem J 199: 155–161, 1981.[Web of Science][Medline]
  7. Eijgelshoven MHJ, van Beek JHGM, Mottet I, Nederhoff MGJ, Vanechteld CJA, and Westerhof N. Cardiac high-energy phosphates adapt faster than oxygen-consumption to changes in heart-rate. Circ Res 75: 751–759, 1994.[Abstract/Free Full Text]
  8. From AHL, Zimmer SD, Michurski SP, Mohanakrishnan P, Ulstad VK, Thoma WJ, and Ugurbil K. Regulation of the oxidative phosphorylation rate in the intact cell. Biochemistry 29: 3731–3743, 1990.[CrossRef][Medline]
  9. Hansford RG. Control of mitochondrial substrate oxidation. Curr Top Bioenerg 10: 217–277, 1980.
  10. 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: H433–H440, 1993.[Abstract/Free Full Text]
  11. Jo H, Noma A, and Matsuoka S. Calcium-mediated coupling between mitochondrial substrate dehydrogenation and cardiac workload in single guinea-pig venticular myocytes. J Mol Cell Cardiol 40:394–404, 2006.
  12. Katz LA, Swain JA, Portman MA, and Balaban RS. Relation between phosphate metabolites and oxygen consumption of heart in vivo. Am J Physiol Heart Circ Physiol 256: H265–H274, 1989.[Abstract/Free Full Text]
  13. Korzeniewski B. Regulation of ATP supply during muscle contraction: theoretical studies. Biochem J 330: 1189–1195, 1998.
  14. Korzeniewski B. Regulation of ATP supply in mammalian skeletal muscle during resting state -> intensive work transition. Biophys Chem 83: 19–34, 2000.[CrossRef][Web of Science][Medline]
  15. Korzeniewski B. Theoretical studies on the regulation of oxidative phosphorylation in intact tissues. Biochim Biophys Acta 1504: 31–45, 2001.[Medline]
  16. Korzeniewski B. Regulation of oxidative phosphorylation in different muscles and various experimental conditions. Biochem J 375: 799–804, 2003.[CrossRef][Web of Science][Medline]
  17. Korzeniewski B, Noma A, and Matsuoka S. Regulation of oxidative phosphorylation in intact mammalian heart in vivo. Biophys Chem 116: 145–157, 2005.[CrossRef][Web of Science][Medline]
  18. Korzeniewski B and Zoladz JA. A model of oxidative phosphorylation in mammalian skeletal muscle. Biophys Chem 92: 17–34, 2001.[CrossRef][Web of Science][Medline]
  19. Korzeniewski B and Zoladz JA. Training-induced adaptation of oxidative phosphorylation in skeletal muscle. Biochem J 374: 37–40, 2003.[CrossRef][Web of Science][Medline]
  20. Korzeniewski B and Zoladz JA. Factors determining the oxygen consumption rate (Formula 4O2) on-kinetics in skeletal muscle. Biochem J 379: 703–710, 2004.[CrossRef][Web of Science][Medline]
  21. Letellier T, Malgat M, and Mazat JP. Control of oxidative-phosphorylation in rat muscle mitochondria: implications for mitochondrial myopathies. Biochim Biophys Acta 1141: 58–63, 1993.[Medline]
  22. McCormack JG, Halestrap AP, and Denton RM. Role of calcium ions in regulation of mammalian intramitochondrial metabolism. Physiol Rev 70: 391–425, 1990.[Free Full Text]
  23. Mootha VK, Arai AE, and Balaban RS. Maximum oxidative phosphorylation capacity of the mammalian heart. Am J Physiol Heart Circ Physiol 272: H769–H775, 1997.[Abstract/Free Full Text]
  24. Portman NA, Heineman FW, and Balaban RS. Developmental changes in the relation between phosphate metabolites and oxygen consumption in the sheep heart in vivo. J Clin Invest 83: 456–464, 1989.[Web of Science][Medline]
  25. Richardson RS, Knight DR, Poole DC, Kurack SS, Hogan MC, Grassi B, and Wagner PD. Determinations of maximal exercise Formula 4O2 during single leg knee-extensor exercise in humans. Am J Physiol Heart Circ Physiol 268: H1453–H1461, 1995.[Abstract/Free Full Text]
  26. 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: 593–603, 2005.[Abstract/Free Full Text]
  27. Tonkonogi M and Sahlin K. Rate of oxidative phosphorylation in isolated mitochondria from human skeletal muscle: effect of training status. Acta Physiol Scand 161: 345–353, 1997.[CrossRef][Web of Science][Medline]
  28. Unitt JF, McCormack JG, Reid D, MacLahlan LK, and England PJ. Direct evidence for a role of intramitochondrial Ca2+ in the regulation of oxidative phosphorylation in the stimulated rat heart. Biochem J 262: 293–301, 1989.[Web of Science][Medline]
  29. Van Beek JGHM, Tian X, Zuurbier CJ, de Groot B, van Echteld CJA, Eijgelshoven MHJ, and Hak JB. The dynamic regulation of myocardial oxidative phosphorylation: analysis of the response time of oxygen consumption. Mol Cell Biochem 184: 321–344, 1998.[CrossRef][Web of Science][Medline]
  30. 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 hearts. Am J Physiol Heart Circ Physiol 265: H453–H460, 1993.[Abstract/Free Full Text]
  31. Whipp BJ, Rossiter HB, and Ward S. Exertional oxygen uptake kinetics: a stamen of stamina? Biochem Soc Trans 30: 237–247, 2002.[CrossRef][Web of Science][Medline]
  32. White RL and Wittenberg BA. Mitochondrial NAD(P)H, ADP, oxidative phosphorylation and concentration in isolated heart cells. Am J Physiol Heart Circ Physiol 279: H1849–H1857, 2000.[Abstract/Free Full Text]
  33. Zhou LF, 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: H2400–H2411, 2005.[Abstract/Free Full Text]
  34. Zhou LF, 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: 925–937, 2005.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Am. J. Physiol. Cell Physiol.Home page
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]


Home page
J. Physiol.Home page
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]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
291/3/H1466    most recent
00004.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Korzeniewski, B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Korzeniewski, B.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2006 by the American Physiological Society.