Previous studies have shown that oscillations of the metabolism can occur in cardiomyocytes under conditions simulating ischemia/reperfusion. It is not known whether they can also occur during real ischemia with near-anoxic oxygen tension. Here, using oxygen clamp in on-chip picochambers, we exposed single resting cardiomyocytes to near anoxia (pO2 < 0.1 mm Hg). We show that at near anoxia, the mitochondrial membrane potential (ΔΨ) was kept by the F1F0-ATPase reversal, using glycolytic adenosine triphosphate (ATP). In many cells, activation of current through sarcolemmal KATP channels (IKATP) started after a delay with one or several oscillations (frequency of 0.044 ± 0.002 Hz). These oscillations were time correlated with oscillations of ΔΨ. Metabolic oscillations at near anoxia are driven by glycolysis because (a) they were inhibited when glycolysis was blocked, (b) they persisted in cells treated with cytoplasmic reactive oxygen species scavengers, and (c) the highest rate of ATP synthesis during an oscillation cycle was associated with the generation of reducing equivalents. Glycolytic oscillations could be initiated upon rapid, but not slow, transition to near anoxia, indicating that the speed of ATP/ADP ratio drop is a determinant of their occurrence. At enhanced oxidative stress, the rate of ATP consumption was increased as indicated by rapid IKATP activation with large-scale oscillations. These results show that metabolic oscillations occur in cardiomyocytes at near anoxia and are driven by glycolysis and modulated by mitochondria through the rate of ATP hydrolysis, which, in turn, can be accelerated by oxidative stress.
INTRODUCTION
When cardiomyocytes are exposed to substrate removal (O’Rourke et al., 1994; Romashko et al., 1998), chemical ischemia (Ryu et al., 2005; Yang et al., 2008), oxidative stress induced by glutathione depletion (Aon et al., 2007; Slodzinski et al., 2008), or photooxidation (Aon et al., 2003), metabolic oscillations are frequently observed. These oscillations are manifested either by cyclic depolarization of the mitochondrial membrane potential (ΔΨ) or cyclic changes of the cytoplasmic ATP/ADP ratio, which, in turn, is reflected by changes in IKATP and cell excitability (O’Rourke et al., 1994; Aon et al., 2003; Ryu et al., 2005; Yang et al., 2008). Metabolic oscillations have been attributed to synchronization of the mitochondrial network of coupled oscillators by cytoplasmic reactive oxygen species (ROS) (Aon et al., 2006). The trigger for oscillations is thought to be oxidative stress due to radicals produced by laser flash–induced photooxidation (Zorov et al., 2000; Aon et al., 2003; Brady et al., 2004) that is supposed to simulate an elevated ROS production by the respiratory chain. During prolonged ischemia, various components of the respiratory chain involved in ROS production become damaged (Chen et al., 2008). ROS accumulation above some “critical” threshold is followed by a subsequent ROS release into the cytosol through a postulated ROS-induced ROS release mechanism, involving the mitochondrial inner membrane anion channel (Aon et al., 2003) or the cyclosporine A–sensitive permeability transition pore (PTP) (Zorov et al., 2000; Brady et al., 2004).
Metabolic oscillations were also observed during chemical ischemia in cardiomyocytes with largely eliminated mitochondrial function, for example, when the mitochondrial “clamping” of ATP/ADP levels was prevented by blocking oxidative phosphorylation with cyanide and rotenone in the presence of a creatine kinase inhibitor (Ryu et al., 2005; Yang et al., 2008). Oscillations were not affected by cytoplasmic ROS scavengers. They were attributed to oscillatory glycolysis under conditions of ADP depletion in the phosphofructokinase reaction (Yang et al., 2008). Finally, the possibility of a cyclic activation of mitochondrial KATP channels has also been suggested (Ryu et al., 2005).
What all studies of metabolic oscillations in cardiomyocytes have in common is that they were performed at unrestricted access of atmospheric oxygen, e.g., with oxygen concentrations >10 times the physiological oxygen tension (Wilson et al., 1988; Rumsey et al., 1990; Gnaiger, 2001). Whether metabolic oscillations also occur at near anoxia is not known. However, the oxygen tension is very important for evaluating the putative role of ROS in the genesis of metabolic oscillations. During ischemia, the oxygen in the affected tissue is consumed rapidly and the oxygen tension drops to near anoxia. Because the mitochondrial ROS production depends on both the degree of reduction of relevant electron carriers of the respiratory chain and the local oxygen availability (Turrens, 2003; Murphy, 2009), it can substantially vary in cardiomyocytes exposed to ischemia with near-anoxic oxygen tensions. ROS production is expected to be increased at more reduced states of the electron carriers in the respiratory chain (Murphy, 2009; Stowe and Camara, 2009), whereas a decreased oxygen availability should result in a decreased ROS production due to less substrate accepting electrons. Indeed, total ROS production of isolated mitochondria was reported to decrease under near anoxia (Turrens, 2003; Hoffman et al., 2007; Korge et al., 2008; Starkov, 2008), despite the fact that the proportion of electrons leaking to oxygen increases with slower respiration at lower oxygen tension (Hoffman and Brookes, 2009). Diminished ROS production at near anoxia most likely contributes to the result that the main ischemic injury happens upon reperfusion, when the oxygen availability is increased again and a burst of the ROS production occurs (Yellon and Hausenloy, 2007).
In this study, we used oxygen clamp in picochambers to address the question of whether or not metabolic oscillations appear in isolated cardiomyocytes exposed to near anoxia, and to specify the anticipated role of ROS in the genesis of these oscillations. We show that metabolic oscillations occur at near anoxia in ischemic cardiomyocytes at times when mitochondrial function is still preserved, i.e., during early ischemia, before ΔΨ collapse and rigor develop. They are driven by glycolysis but could be affected by mitochondria through the changes in the rate of ATP hydrolysis leading to changes of ATP/ADP ratio. During early ischemia, at steady-state near anoxia ROS does not importantly contribute. However, under conditions favoring oxidative stress, cellular radical formation can contribute to this control by affecting the mitochondrial ATP hydrolysis rate and, hence, the ATP/ADP ratio. Glycolytic oscillations could be triggered by rapid, but not slow, transition to near anoxia, indicating that the rate of oxygen consumption in the ischemic heart can be a determinant of their occurrence.
MATERIALS AND METHODS
Preparation of mouse ventricular cardiomyocytes
Animals were housed and handled according to the guidelines of the German Animal Protection law, and the experiments were conducted in accordance with the Guide for the Care and Use of Laboratory Animals (1996. National Academy of Sciences, Washington, D.C.). Mouse ventricular cardiomyocytes were prepared as described previously (Ganitkevich et al., 2006).
Picochambers and experimental setup
The experimental setup was described in detail previously (Ganitkevich et al., 2006). In brief, picochambers were constructed on glass chips by coating one side with a 32-µm thick layer of a negative photoresist in which cuboid-like pits of 150-µm length, 40-µm width, and 32-µm depth were created, resulting in a geometric volume of 192 pL. The glass chips formed the bottom of the experimental chamber. In the headspace of the experimental chamber, an atmosphere of humidified argon was established to insulate the picochambers from the atmospheric air while having free access for glass pipettes. pO2 in the solution of the picochamber could be controlled by setting the argon flow. When changing the argon flow rate, pO2 at the cell could be changed rapidly in a step-like fashion (<5 s) in the range between <0.1 and 10 mm Hg, which is the critical range for control of oxidative metabolism (Wilson et al., 1988; Rumsey et al., 1990; Gnaiger, 2001). Hence, the argon flow rate defined the oxygen tension at the bottom of the experimental chamber. If not otherwise mentioned, the measurements were performed at the two argon flow rates of 30 and 100 ml/min. With an argon flow rate of 30 ml/min, the pO2 at the bottom of the chamber was 10 mm Hg, as measured with either a Clark electrode or Pt-porphyrin quenching (MFPF-100-1 with FOSPOR-R sensor; Ocean Optics). With an argon flow rate progressively higher, the pO2 levels became, respectively, lower (Fig. S1 A). Because below pO2 = 2 mm Hg the Clark electrode became progressively erroneous, the pO2 at the argon flow rate of 100 ml/min was preferably measured by Pt-porphyrin quenching, resulting in a value between 0.04 and 0.10 mm Hg. The value of 0.04 mm Hg was obtained as follows. A further increase of the argon flow rate to 200 ml/min led to a further slowing of the phosphorescence lifetime, indicating a further decrease of pO2 (Fig. S1 B), whereas no further increase of the lifetime was observed when the sensor was exposed to 100% nitrogen outside the setup. This indicates that any further reduction of pO2 was below the resolution of the sensor of 0.04 mm Hg. Hence, with an argon flow rate of 100 ml/min, the cardiomyocyte in a picochamber was exposed to pO2 between 0.04 and 0.1 mm Hg. Routinely, pO2 was monitored with the Clark electrode because the measurement with Pt-porphyrin could interfere with the fluorescence measurement of ΔΨ. We conclude that changing the argon flow rate between 30 and 100 ml/min varied the oxygen tension at the cell between 10 and <0.1 mm Hg, respectively. Rapid (5 s) and slow (minutes) changes of pO2 were achieved by applying step and ramp voltages to the flow controller (Smart- Track-100; Sierra Instruments).
Patch clamp technique
Transmembrane currents of cardiomyocytes were recorded with the patch clamp technique in the whole cell configuration as described previously (Ganitkevich et al., 2006). If not otherwise noted, the holding potential was −80 mV and the cells were depolarized to +40 mV for 100 ms at 2 Hz. All experiments were performed at room temperature (22 ± 1°C).
Fluorescence optical measurements
ΔΨ was measured optically (Ganitkevich et al., 2006) by using the fluorescence dye tetramethylrhodamine methyl (or ethyl) ester (TMRM or TMRE). Most experiments were performed in the de-quench mode (Nicholls and Ward, 2000; Duchen et al., 2003). The cells were incubated with a KB solution containing 2 µM of either TMRM or TMRE at room temperature for 15 min. Under the experimental conditions used, the two dyes provided indistinguishable results. For most experiments, TMRM was used in the de-quench mode, which provided two advantages. First, the sensitivity is ∼100 times higher compared with the nonquenching mode (Voronina et al., 2004). Second, because the photomultiplier (PMT) signal was collected from the whole picochamber, where all dye was trapped, the diffusional loss via the pipette was only minimal, avoiding redistribution problems with the dye, as present in experiments with large solution volumes (Nicholls and Ward, 2000). However, because of concerns about a high matrix concentration of the dye, some control experiments were done with low nanomolar concentrations. For these experiments in the nonquench mode, the cells were incubated in Tyrode’s solution containing 20 nM TMRM for 15 min and imaged every 5 s with an Apotome microscope (Carl Zeiss, Inc.). “Pseudoconfocal” images were processed with the ImageJ software (National Institutes of Health). To avoid bleaching of the dye in long-lasting experiments, the development of rigor was accelerated by the inclusion of low concentrations (50 or 100 nM) of the uncoupler FCCP in the patch pipette solution. To estimate the cytosolic [ATP], free [Mg2+] was measured with MgGreen as described previously (Leyssens et al., 1996). The cells were loaded with 5 µM MgGreen-AM for 30 min. 2′,7′-dichlorodihydrofluorescein (H2DCF) oxidation was measured as described previously (Ganitkevich et al., 2006). Cells were loaded with 5 µM CM-H2DCFDA for 20 min. NAD(P)H fluorescence could not be measured in these experiments because UV illumination destroyed the photoresist used for chip production. Flavoprotein fluorescence signals were too small to be usable as metabolic parameter during metabolic oscillations.
Data acquisition and analysis
Measurements were controlled and data were collected and analyzed with the ISO2/ISO3 software and hardware (12-bit resolution; MFK) running on a PC. The sampling rate was 1 kHz. The results are presented as mean ± SE.
Chemicals and solutions
All chemicals were of analytical grade. Tyrode’s solution contained (in mM): 140 NaCl, 5.4 KCl, 1.2 CaCl2, 0.5 MgCl2, 0 glucose, and 10 HEPES, pH 7.4 (NaOH). Intracellular solution in the patch pipette contained (in mM): 140 KCl, 10 HEPES, and 5 EGTA, pH 7.3 (KOH). The cells were stored for experimental use in KB medium containing (in mM): 30 KCl, 30 KH2PO4, 50 glutamine acid, 20 taurine, 20 HEPES, 10 glucose, 3 MgSO4, and 0.5 EGTA, pH 7.3 (KOH). Drugs were applied by flooding the picochamber with the appropriate solution. Normally, the PMT was switched off during solution exchanges to avoid the recording of artifacts. Recording was continued after the solution surrounding the cardiomyocyte was reduced to the picochamber volume again.
Online supplemental material
Fig. S1 shows the dependence of pO2 at the bottom of the experimental chamber on the rate of Ar flow. Oxygen quenching of Pt-porphyrin luminescence (FOSPOR-R sensor with MFPF-100-1; Ocean Optics) was measured at different Ar flow rates. Fig. S2 demonstrates that glycolytic ATP is required to sustain ΔΨ and to keep KATP channels closed at near anoxia. Glycolysis was impaired with iodoacetic acid (IAA) or 2-deoxyglucose (2-DG), whereas oligomycin or bongkrekic acid (BKA) was used to inhibit the F1F0-ATPase.
RESULTS
ΔΨ is maintained during IKATP activation in resting cardiomyocytes exposed to near anoxia
Resting cardiomyocytes (5 mM EGTA in the patch pipette) responded to a change from atmospheric pO2 (160 mm Hg) to near anoxia (pO2 < 0.1 mm Hg) with only small depolarization of ΔΨ, which was regularly <10% of the signal obtained with the uncoupler FCCP (Fig. 1 A). After a variable delay (usually 5–10 min), IKATP became activated, suggesting that the cytoplasmic ATP/ADP ratio could not be maintained. This was confirmed by measurements of the cytosolic [Mg2+] with MgGreen. At near anoxia, the free [Mg2+] became increased when IKATP was activated. Time courses of IKATP and MgGreen signal dissociated, which is likely due to a progressive decrease of pH at near anoxia (Ganitkevich et al., 2006), affecting fluorescence properties of the indicator. Reoxygenation, which restored ATP synthesis by oxidative phosphorylation, rapidly removed IKATP and decreased the [Mg2+] (n = 6; Fig. 1 B). Most notably, ΔΨ was kept largely polarized under conditions when the cytosolic [ATP] significantly dropped, as indicated by the activation of IKATP (Fig. 1 A). However, in the de-quench mode, the relation of the TMRM fluorescence signal to ΔΨ is presumably not simple. We therefore performed experiments in the nonquench mode, in which the fluorescence signals are likely to be more directly related to ΔΨ (Duchen et al., 2003). Activation of IKATP at near anoxia was not associated with a measurable signal (Fig. 1 C) until rigor contracture occurred and ΔΨ collapsed (n = 9). Notably, the mitochondria located under the patch pipette kept some ΔΨ, most likely due to an oxygen leak from the patch pipette and accumulated dye that was released from more remote mitochondria where ΔΨ collapsed. This shows that any residual ATP production by these mitochondria did not prevent development of the rigor contracture (Fig. 1 C).
To prove that glycolytic ATP is needed for keeping ΔΨ polarized at near anoxia, glycolysis was manipulated. The glycolytic ATP production was either inhibited by IAA or the cytoplasmic ATP consumption was accelerated by 2-DG. In the presence of IAA, ΔΨ was sustained at pO2 = 10 mm Hg, but rapidly and almost completely collapsed at near anoxia (Fig. S2 A). The results with 2-DG were similar to those obtained with IAA (Fig. S2 B).
To prove that the reversal of the F1F0-ATPase keeps ΔΨ at near anoxia at the expense of glycolytic ATP, cells were treated with oligomycin. In the presence of oligomycin, IKATP appeared not only at near anoxia, but after a delay also at pO2 = 10 mm Hg, suggesting that ATP synthesis was inhibited. At the same time, ΔΨ was largely sustained because of an intact respiratory chain. After activation of IKATP, stepping to near anoxia led to the collapse of ΔΨ (Fig. S2 C). Finally, the adenine nucleotide translocator (ANT) inhibitor BKA was applied to diminish ATP hydrolysis by the F1F0-ATPase (Nicholls and Ferguson, 2002). In the presence of BKA, the ability of the cell to sustainΔΨ at near anoxia became progressively reduced (Fig. S2 D).
Collectively, these results show that the activation of IKATP indicates a progressive imbalance between glycolytic ATP synthesis and hydrolysis of ATP by the reverse action of the F1F0-ATPase.
At near anoxia, oscillations of IKATP and ΔΨ occur conjointly
Activation of IKATP started in many cells (34% of 198) with one or more oscillations (Fig. 2 A). They occurred always at the beginning of IKATP activation, ceased when the current activated further, were paralleled by small oscillations of ΔΨ, and had a frequency of 0.044 ± 0.002 Hz (n = 37). IKATP oscillations also occurred in cells not loaded with TMRM (n = 8; Fig. 2 B). Comparison of the time courses of IKATP and ΔΨ during oscillations revealed an intimate temporal correlation (Fig. 2, C and D). Because at near anoxia both IKATP and ΔΨ are linked to the cytoplasmic ATP/ADP ratio, these results suggest that the cytoplasmic ATP/ADP ratio oscillates, resulting in oscillations of both IKATP and the ability of the F1F0-ATPase to keep ΔΨ. These results also suggest that in resting cardiomyocytes with an intact creatine kinase system, no significant spatial gradients in the ATP/ADP ratio were present at the time scale of metabolic oscillations.
ROS scavengers and modulators of the inner mitochondrial membrane conductance do not prevent metabolic oscillations
Only low ROS levels are expected at near anoxia. To test an involvement of ROS in the control of metabolic oscillations, cells were pretreated with redox-active compounds. Oscillatory behavior of IKATP and ΔΨ was still observed in cells preloaded with various cytoplasmic ROS scavengers or superoxide dismutase mimetics like TEMPOL, N-(2-mercaptopropionyl)-glycine (MPG), 3-methyl-1-phenyl-2-pyrazolin-5-one (MCI-186), and Mn(III)-tetrakis-(4-benzoic-acid)-porphyrin (MnTBAP; all 5 mM; ≥60-min pretreatment; mean frequency 0.039 ± 0.003 Hz; n = 22). Hence, at near anoxia, ROS are not likely to contribute significantly to metabolic oscillations.
To test for an involvement of mitochondrial KATP channel and PTP, appropriate ligands were used. Metabolic oscillations persisted in cells pretreated with 100 µM diazoxide (mitoKATP channel opener), 100 µM 5-hydroxydecanoate (5-HD; mitoKATP channel blocker), or 1 µM cyclosporine A (PTP blocker; mean frequency 0.035 ± 0.002 Hz; n = 12). These findings make the involvement of mitoKATP channels and the PTP in the control of metabolic oscillations at this early stage of ischemia unlikely.
Anion channel blockers (4,4’-diisothiocyanostilbene-2,2’-disulfonic acid [DIDS] and 4′-chlorodiazepam), which were previously shown to reduce reperfusion injury (Akar et al., 2005), induced a rapid and complete collapse of ΔΨ at near anoxia (unpublished data), consistent with reports that inhibition of voltage-dependent anion-selective channel and/or ANT inhibits the ATP transport into the mitochondrial matrix, preventing ATP hydrolysis by the F1F0-ATPase required to keep ΔΨ (Báthori et al., 2006).
The results so far show that metabolic oscillations do occur in cardiomyocytes at near anoxia, and they suggest that glycolysis is the key player. They do not reveal, however, whether and how mitochondria are involved.
Mitochondrial state affects metabolic oscillations at near anoxia
To gain further insight into the mechanisms underlying the metabolic oscillations at near anoxia, the development of IKATP was interrupted by short (30-s) intervals of reoxygenation (pO2 = 10 mm Hg). Reoxygenation allows the cells to respire, hyperpolarize ΔΨ, and produce ATP by oxidative phosphorylation, thus “clamping” the ATP/ADP ratio to a high level. Return to near anoxia again inhibits respiration and depolarizes ΔΨ, but at a higher ATP/ADP level, because ATP is produced during the preceding reoxygenation interval.
Most notably, shortly after a rapid (5-s) return to near anoxia, IKATP rapidly activated, and to a higher level than before reoxygenation, despite some ATP that must have been synthesized during the preceding reoxygenation interval (Fig. 3 A). Moreover, when the protocol was applied shortly before measurable IKATP activation, a burst of IKATP was recorded (Fig. 3 B). Repetitive application of this protocol during activation of IKATP induced progressively larger bursts of IKATP, which switched over to oscillations (Fig. 3 B). When ΔΨ was recorded together with IKATP, the IKATP burst was paralleled by a transient depolarization of ΔΨ (Fig. 3 C). The conjoint oscillations of both ΔΨ and IKATP could be evoked (Fig. 3 D, after first reoxygenation interval) or, if already present, potentiated (Fig. 3 D, after second reoxygenation interval) by a rapid return to near anoxia. ROS scavengers did not prevent potentiation effects of reoxygenation intervals (unpublished data).
These results suggest that shortly after a rapid return to near anoxia after a reoxygenation interval, the hydrolysis of ATP transiently increases. Because the effect was not sensitive to ROS scavengers, we hypothesize that it could be due to a switch of the F1F0-ATPase from ATP synthesis to ATP hydrolysis as the result of rapid mitochondrial depolarization.
ΔΨ is a determinant of the ATP hydrolysis rate at near anoxia
To prove that ΔΨ determines the ATP hydrolysis rate by the F1F0-ATPase, glycolysis was blocked by 1 mM IAA and the F1F0-ATPase was reversed by near anoxia, leading to cytoplasmic ATP consumption until the development of rigor and collapse of ΔΨ. As in cells with intact glycolysis (Fig. 2 A), ΔΨ was depolarized only slightly upon transition to near anoxia (Fig. 4 A). Within less than a minute, however, IKATP was completely activated and ΔΨ collapsed (n = 5), suggesting that with glycolytic ATP supply blocked, ATP was rapidly consumed by the F1F0-ATPase (compare Fig. 2 A). Reoxygenation restored ΔΨ and removed IKATP after a delay, when oxidative phosphorylation had produced enough ATP (Fig. 4 A). This shows that the respiratory chain was still functional. Notably, at near anoxia, ΔΨ became depolarized at an increasingly higher rate with a maximum close to time of complete activation of IKATP (Fig. 4 A). Hence, depolarized mitochondria reduce the ATP/ADP ratio at a higher rate. The rate of ΔΨ depolarization decreased again when the ATP/ADP ratio fell further, as indicated by the collapse of ΔΨ and the development of rigor. Reoxygenation restored ΔΨ, indicating that at this time, the respiratory chain is not significantly damaged (Fig. 4 A). At blocked glycolysis, metabolic oscillations were never observed (n = 5).
To directly demonstrate that the rate of ATP hydrolysis by mitochondria affects metabolic oscillations, we applied a low concentration of an uncoupler FCCP to cells with intact glycolysis to reduce ΔΨ kept at near anoxia and thus to accelerate ATP hydrolysis by the F1F0-ATPase.
In the presence of 20 nM FCCP, transition to near anoxia induced large depolarization of ΔΨ and an immediate maximal activation of IKATP with large-scale oscillations (n = 12; Fig. 4 B). Similar to cells not treated with an uncoupler, these oscillations of IKATP were potentiated after intervals of reoxygenation (Fig. 4 B). However, the amplitude of IKATP oscillations was much larger compared with that usually recorded in untreated cells (Figs. 2 and 3). Oscillations of IKATP were also paralleled by oscillations of ΔΨ, until its collapse due to the fall of the ATP/ADP ratio. Nevertheless, the respiratory chain was still able to polarize ΔΨ upon reoxygenation (Fig. 4 B).
These results further confirm the finding that glycolysis is responsible for metabolic oscillations at near anoxia, that mitochondria are suppliers of ADP for the phosphofructokinase, and that this supply depends on ΔΨ kept at near-anoxic conditions. Moreover, because ΔΨ is a determinant of the ATP hydrolysis rate, the speed of the fall of the ATP/ADP ratio should also depend on the rate of ΔΨ depolarization. This could be verified experimentally.
The rate of oxygen drop determines mitochondrial influence on metabolic oscillations
We induced both a rapid (5-s) and a slow drop of oxygen to near anoxia in the same experiment. With the rapid drop, oscillations of both IKATP and ΔΨ were observed (Fig. 5 A; n = 5). Reoxygenation restored ΔΨ and removed IKATP. When the pO2 was reduced slowly to near anoxia within 2 min, IKATP was activated again and ΔΨ depolarized at a slower rate, but oscillations of ΔΨ and IKATP were largely absent (Fig. 5 A). Similar results were obtained when protocols were applied in the opposite order (n = 5). Similarly, a stepwise reduction of pO2 to near anoxia within 3 min depolarized ΔΨ and activated IKATP without inducing oscillations (n = 5; Fig. 5 B). Reoxygenation restored ΔΨ and abolished IKATP; subsequent rapid transition to near anoxia again evoked oscillations of both IKATP and ΔΨ.
These results suggest that the speed of transition to near anoxia, i.e., the rate of ΔΨ depolarization and fall of the ATP/ADP ratio is an important factor for the genesis of glycolytic oscillations.
Oxidative stress modulates metabolic oscillations
We found no role of ROS at near anoxia, as expected with little ROS production at a very low pO2. This conclusion is based on the failure of scavengers to prevent metabolic oscillations. Instead, all results are in support of the idea that glycolysis is a critical determinant of the metabolic oscillations. Nevertheless, ROS might be important under conditions of myocardial ischemia because they are expected to be produced during reperfusion, causing damage of the respiratory chain. We therefore studied whether and how cellular ROS production affects metabolic oscillations. Doing this, we were aware that measurement of ROS in cells is a challenging task because no really suitable tools are available (Dikalov et al., 2007; Wardman, 2007). Commonly used H2DCF is oxidized in the cell to the fluorescent compound DCF, and the increase of green fluorescence is used as an indicator of ROS production. However, it was shown in several insightful studies that upon illumination, DCF could act as photosensitizer for oxidation of cellular reductants (reduced glutathione [GSH] and NADH), thus producing radicals itself (Marchesi et al., 1999; Rota et al., 1999a,b). Subsequent electron transfer from these radicals to molecular oxygen produces superoxide (Bonini et al., 2006; O’Neill and Wardman, 2009). Moreover, as DCF is a two-electron oxidation product of H2DCF, an intermediate radical is formed during oxidation, which reacts with oxygen to produce superoxide directly (Bonini et al., 2006; Wardman, 2007; Wrona et al., 2008). Therefore, DCF itself induces oxidative stress in a self-amplifying oxygen-consuming chain reaction, which depends on light, oxygen concentration, and the redox status of the cell (Wardman, 2007). We used these properties of DCF, in combination with substances having well-established effects on mitochondrial ROS production, to assess how different degrees of oxidative stress affect metabolic oscillations, reflected by IKATP.
Because in an experiment the DCF fluorescence continuously increased upon illumination, as expected for a self-amplifying chain reaction, we used the rate of fluorescence increase (dFDCF/dt) as the measure of oxidative stress (Fig. 6).
When the cells were treated with the complex III Qo site inhibitor stigmatellin (1 µM; n = 7), which blocks ROS production (Nicholls and Ferguson, 2002; Turrens, 2003), the rate of H2DCF oxidation increased only very slowly and IKATP became activated slowly and showed oscillatory behavior (Fig. 6 A). This result suggests that blocking mitochondrial ROS production at complex III did not prevent glycolytic oscillations.
When cells were treated with the complex III Qi site inhibitor antimycin (1 µM; n = 5), which is known to greatly enhance ROS production by the respiratory chain (Nicholls and Ferguson, 2002; Turrens, 2003; Starkov, 2008), the rate of H2DCF oxidation increased rapidly and continuously (Fig. 6 B). IKATP activated much faster compared with the activation in the presence of stigmatellin and showed oscillatory behavior. Rapid transition from pO2 = 10 mm Hg to near anoxia induced a transient increase of the H2DCF oxidation rate lasting ∼20 s (Fig. 6 B). It was shown before that this acceleration of H2DCF oxidation is due to an electron transfer from thiyl radicals formed during preceding oxidative stress, which was suppressed in the presence of oxygen (Korge et al., 2008; Wrona et al., 2008). This is supported by the result that transient acceleration of H2DCF oxidation was nearly absent with stigmatellin with mitochondrial ROS production blocked, indicating a much lower degree of oxidative stress in the presence of stigmatellin.
When the cells were treated with the complex I inhibitor rotenone (1 µM; n = 5), which facilitates the ROS production at the FMN site of this complex (Nicholls and Ferguson, 2002; Turrens, 2003; Starkov, 2008), the rate of H2DCF oxidation increased continuously, and an H2DCF oxidation burst similar to that observed with antimycin was seen upon transition to near anoxia (Fig. 6 C). IKATP activated quickly and showed oscillatory behavior. The degree of oxidative stress, estimated from the rate of the H2DCF oxidation at pO2 = 10 mm Hg after the addition of the blocker, was highest with antimycin, followed by rotenone. The results with stigmatellin were not significantly different from control (Fig. 6 D).
Collectively, our experiments demonstrate that elevated ROS levels are not necessary for the oscillatory behavior of IKATP, through enhanced radical formation results in a faster fall of the ATP/ADP ratio, i.e., in an acceleration of ATP hydrolysis. Our results also demonstrate the problems of using DCF within cells. Thiyl (as well as other) radicals use oxygen as an electron sink to form superoxide, which is then dismutated in its usual way (Winterbourn, 1993, 2008; Winterbourn and Hampton, 2008), so that the lowering of pO2 results in a redistribution of unpaired electrons between other acceptors. Therefore, interpretation of DCF fluorescence signals from cells is difficult, especially when pO2 is varied. We therefore restricted our measurements with DCF to a constant low pO2 of <0.1 mm Hg. We aimed to relate changes of the ATP/ADP ratio, reflected by IKATP to the cellular redox state, during metabolic oscillations at near anoxia.
Oscillations of the H2DCF oxidation rate and IKATP at near anoxia are closely time related
At near anoxia, little superoxide is produced (Turrens, 2003; Korge et al., 2008; Muller, 2009). In the absence of oxygen, H2DCF oxidation arises from thiyl radicals (Korge et al., 2008; Wrona et al., 2008) formed by photoreduction of DCF (Marchesi et al., 1999). Like other fluorescent dyes, DCF is activated by visible light to a triplet state that is reactive toward reducing agents like GSH, generating thiyl radicals (Marchesi et al., 1999). Notably, NAD(P)H can substitute for GSH and photoreduce the dye even much faster (Marchesi et al., 1999). Because we could not measure NAD(P)H fluorescence directly (see Materials and methods), we used this peculiarity of DCF to compare IKATP and the H2DCF oxidation rate, which should be related to [NAD(P)H] and hence to the cellular redox state. Because the glycolytic ATP flux is accompanied with the production of NAD(P)H, we could assume that that oscillations of the ATP/ADP ratio are also reflected by DCF fluorescence.
Upon transition to near anoxia, a small burst of H2DCF oxidation was observed, induced by thiyl radicals formed during preceding hyperoxia (160 mm Hg). IKATP started to oscillate after a delay with a frequency of 0.043 ± 0.002 Hz (n = 6; Fig. 7 A). These oscillations were accompanied by small changes of the H2DCF oxidation rate, indicating that the cellular redox state also oscillated. When IKATP became further activated, oscillations of both IKATP and the H2DCF oxidation rate ceased (Fig. 7 A). Comparison of time courses of IKATP and H2DCF oxidation during metabolic oscillations at higher time resolutions revealed that the fastest H2DCF oxidation rate was paralleled by the highest rate of IKATP activation and, vice versa, the lowest H2DCF oxidation rate was paralleled by the highest rate of KATP channel closure (Fig. 7 B). These results suggest that at near anoxia, the ATP/ADP ratio is inversely related to the cellular redox state, i.e., ATP synthesis is accompanied with the generation of reducing equivalents. This further suggests that ROS are not casually involved in the generation of metabolic oscillations at near anoxia.
DISCUSSION
Metabolic oscillations in isolated cardiomyocytes were observed under conditions simulating ischemia (Ryu et al., 2005; Yang et al., 2008) or reperfusion (Aon et al., 2003), but always at atmospheric oxygen. They were attributed either to oscillatory behavior of glycolysis during chemical ischemia, without obvious involvement of ROS (Yang et al., 2008), or to oxidative stress (ROS) -induced synchronization of metabolism during reperfusion (Akar et al., 2005; Aon et al., 2006). The results of the present study demonstrate that metabolic oscillations also occur in ischemic cardiomyocytes exposed to near anoxia, that they are driven by glycolysis, and that ROS are not essentially involved at this stage of ischemia. Because our experiments were performed under near-anoxic conditions, the metabolic state of a cardiomyocyte at near anoxia, in particular that of the mitochondria, is considered first.
Both glycolysis and mitochondria control the metabolic state of a cardiomyocyte at near anoxia
At normoxia, when ΔΨ is high, ATP synthesis by the F1F0-ATPase decreases ΔΨ by using the proton electrochemical potential. At near anoxia, the respiratory chain is not able to keep ΔΨ, resulting in a depolarization (Fig. 1 A). When the proton electrochemical potential, which is mainly ΔΨ, becomes lower than the free energy of the phosphorylation potential, the F1F0-ATPase switches from ATP synthesis to ATP hydrolysis, resulting in proton pumping. Together with the ANT, which translocates ATP4− into and ADP3− out of the matrix (net negative charge in), proton pumping by the F1F0-ATPase results in maintenance of ΔΨ at a polarized level (Fig. 1 A). This level of ΔΨ is maintained as long as the glycolytic ATP flux, the only source of ATP at near anoxia, supplies enough ATP to the mitochondria to keep the free energy of ATP hydrolysis below the electrochemical proton potential (Fig. 2 A). One determinant of the latter is the leak in the inner mitochondrial membrane. Because ATP4− must enter the matrix against the ΔΨ to become hydrolyzed, mitochondria that are more depolarized consume ATP at a higher rate (Fig. 4 A). Thus, ATP is hydrolyzed at the highest rate by uncoupled mitochondria (Nicholls and Ferguson, 2002).
Mitochondria are involved in metabolic oscillations at near anoxia
The results of this study demonstrate that at the beginning of IKATP activation, oscillatory phenomena occur (Fig. 2 A). Because activation of IKATP represents an imbalance between the glycolytic ATP flux and ATP consumption by the F1F0-ATPase, oscillations of the ATP/ADP ratio can be controlled in principle at the level of both glycolysis and mitochondria. An oscillatory behavior requires a feedback. One such feedback is the link between ΔΨ and the ATP hydrolysis rate. When the ATP concentration in the cytoplasm falls, the proton-pumping ability of the F1F0-ATPase decreases and ΔΨ gets more depolarized. If the glycolytic ATP flux cannot balance this, the result would be a further ATP fall and ΔΨ depolarization; i.e., there is a positive feedback between depolarization of ΔΨ and the ATP consumption rate (Fig. 4 A). This positive feedback is likely to be responsible for the rapid initial ATP/ADP ratio drop with the resulting activation of IKATP (Fig. 1 A). Accelerated ATP hydrolysis stimulates the glycolytic ATP flux, resulting in a slowing down of the further ATP fall (Fig. 1 A) or even transient recovery of the ATP/ADP ratio and ΔΨ (Fig. 2 C). Hence, as long as glycolysis supplies enough ATP, the positive feedback of ΔΨ on the ATP hydrolysis rate can be interrupted by a repetitive increase of the glycolytic ATP flux. Hence, as long as ΔΨ oscillates, the ATP consumption rate should also oscillate, resulting in oscillations of the ATP/ADP ratio. When oscillations of ΔΨ are small, it is likely that this feedback is of lesser significance for the control of glycolytic oscillations. Hence, it is apparently important only for the initiation of the oscillations, as it is responsible for the acceleration of the ATP/ADP ratio fall.
The other, more important positive feedback leading to oscillations of the cytoplasmic ATP/ADP ratio is driven by the effect of free ADP on the phosphofructokinase activity (Yang et al., 2008). In this study, metabolic oscillations were recorded at atmospheric oxygen in the presence of cyanide or rotenone, i.e., under conditions simulating ischemia at high oxygen. It was shown that in cardiomyocytes, when the capacity of oxidative phosphorylation and the creatine kinase system to buffer the cellular ATP/ADP ratio is suppressed, glycolysis could cause periodic oscillations in cellular ATP levels. As in the present study, metabolic oscillations were suppressed by IAA, persisted in cardiomyocytes pretreated with ROS scavengers, and showed an anti-phase relationship to redox oscillations (Yang et al., 2008). Notably, the frequency of metabolic oscillations in cardiomyocytes at 36°C was with 0.02–0.067 Hz, comparable to the respective frequency in our experiments at 22°C (0.044 Hz), suggesting a low Q10 value for the frequency-controlling process. Alternatively, the frequency-controlling process at high and low oxygen might be different.
The rate of oxygen drop affects metabolic oscillations through the rate of ATP hydrolysis
A key finding of our report is that the mitochondrial effect on metabolic oscillations depends on the speed of the oxygen drop (Fig. 5). As oxygen tension decreases, the respiratory chain slows. Hence, the speed of oxygen drop is the determinant of the rate of mitochondrial depolarization. In turn, the ΔΨ is a determinant of the ATP hydrolysis rate (Fig. 4 A). When the mitochondria are depolarized at a fast rate (seconds), metabolic oscillations were either elicited or potentiated (Figs. 3 and 4 B). When mitochondria are depolarized at a slow rate (minutes), the oscillatory phenomena were much less prominent because glycolysis has time to compensate for the slower fall of the ATP/ADP ratio (Fig. 5). When the blood flow is interrupted during acute ischemia in the heart, the local oxygen tension drops rapidly because oxygen is consumed by cells. Hence, the rapid (seconds) drop of oxygen in our experiments is likely to mirror acute ischemia in vivo. These results suggest that the rate of oxygen consumption in the ischemic heart is a determinant for the occurrence of glycolytic oscillations.
Lack of significant involvement of ROS and mitoKATP in metabolic oscillations during early ischemia
In theory, near-anoxic conditions, i.e., highly reduced conditions, favor ROS generation by complex I of the respiratory chain (Murphy, 2009). Superoxide generated at the matrix side (Turrens, 2003; Starkov, 2008) was shown to induce permeabilization of the inner mitochondrial membrane, leading to depolarization of ΔΨ (Echtay et al., 2002a), possibly by activating uncoupling proteins (Echtay et al., 2002b). Induction of oxidative stress did accelerate an ATP/ADP ratio drop (Fig. 6). However, at near anoxia, we did not find evidence for ROS contributing importantly at this early stage of ischemia, i.e., when ΔΨ is still maintained by glycolysis and the respiratory chain is not noticeably damaged. Depolarization of ΔΨ due to the opening of mitochondrial KATP channels was shown to be cardioprotective (Korge et al., 2002). However, we did not observe any effects of pretreatment of cardiomyocytes with both opener and blocker of mitochondrial KATP channels in our experiments. However, mitoKATP channels might contribute importantly during ischemic events of longer duration and ROS during hypoxic ischemia and reperfusion (Robin et al., 2007; Korge et al., 2008).
In conclusion, we show the occurrence of glycolytic oscillations during early ischemia under near-anoxic conditions. Mitochondria affect these oscillations through the speed of ATP utilization by the backwardly running F1F0-ATPase, and the rate of ATP hydrolysis can be modulated by uncoupling or oxidative stress. In the ischemic heart, the speed of oxygen consumption by the tissue determines, through the speed of mitochondrial depolarization, the rate of ATP hydrolysis, which in turn determines the occurrence of oscillations in metabolism and excitability.
Acknowledgments
We thank G. Ditze, G. Sammler, F. Horn, and S. Bernhardt for excellent technical assistance.
This work was supported by the Deutsche Forschungsgemeinschaft (BE 1250/15-3) to K. Benndorf.
Angus C. Nairn served as editor.