Mitochondrial Contributions in the Genesis of Delayed Afterdepolarizations in Ventricular Myocytes

Mitochondria fulfill the cell’s energy demand and affect the intracellular calcium (Ca²⁺) dynamics via direct Ca²⁺ exchange, the redox effect of reactive oxygen species (ROS) on Ca²⁺ handling proteins, and other signaling pathways. Recent experimental evidence indicates that mitochondrial depolarization promotes arrhythmogenic delayed afterdepolarizations (DADs) in cardiac myocytes. However, the nonlinear interactions among the Ca²⁺ signaling pathways, ROS, and oxidized Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) pathways make it difficult to reveal the mechanisms. Here, we use a recently developed spatiotemporal ventricular myocyte computer model, which consists of a 3-dimensional network of Ca²⁺ release units (CRUs) intertwined with mitochondria and integrates mitochondrial Ca²⁺ signaling and other complex signaling pathways, to study the mitochondrial regulation of DADs. With a systematic investigation of the synergistic or competing factors that affect the occurrence of Ca²⁺ waves and DADs during mitochondrial depolarization, we find that the direct redox effect of ROS on ryanodine receptors (RyRs) plays a critical role in promoting Ca²⁺ waves and DADs under the acute effect of mitochondrial depolarization. Furthermore, the upregulation of mitochondrial Ca²⁺ uniporter can promote DADs through Ca²⁺-dependent opening of mitochondrial permeability transition pores (mPTPs). Also, due to much slower dynamics than Ca²⁺ cycling and ROS, oxidized CaMKII activation and the cytosolic ATP do not appear to significantly impact the genesis of DADs during the acute phase of mitochondrial depolarization. However, under chronic conditions, ATP depletion suppresses and enhanced CaMKII activation promotes Ca²⁺ waves and DADs.

Introduction

Delayed afterdepolarizations (DADs) are abnormal depolarizations during the diastolic phase following an action potential (AP) and could trigger cardiac arrhythmias (Rosen et al., 1984; January and Fozzard, 1988; Katra and Laurita, 2005; Qu et al., 2014; Song et al., 2017). DADs are known to be caused by spontaneous calcium (Ca²⁺) waves (Rosen et al., 1984; Marban et al., 1986; January and Fozzard, 1988), occurring due to spontaneous Ca²⁺ release from the intracellular Ca²⁺ store, sarcoplasmic reticulum (SR), via the ryanodine receptors (RyRs). Ca²⁺ waves are known to be promoted by Ca²⁺ overload under normal (Cheng et al., 1996) and pathological conditions, such as heart failure (Pogwizd and Bers, 2003; Hoeker et al., 2009), long QT syndrome (Mohler et al., 2003), ischemia (Ross and Howlett, 2009), and catecholaminergic polymorphic ventricular tachycardia (CPVT) (Watanabe et al., 2009). During a cardiac cycle, Ca²⁺ enters into the cytosol from the extracellular space mainly via L-type Ca²⁺ channels (LCCs) during membrane depolarization, which causes Ca²⁺ release from the SR, a process called Ca²⁺-induced Ca²⁺ release (CICR; Bers, 2002). Ca²⁺ is extruded from the cell mainly through the Na⁺-Ca²⁺ exchanger (NCX) and taken back to the SR through sarcoplasmic reticulum Ca²⁺-ATPase (SERCA). Meanwhile, mitochondria, as another Ca²⁺ store, are involved in intracellular Ca²⁺ cycling via mitochondrial Ca²⁺ uniporter (MCU; Baughman et al., 2011; De Stefani et al., 2011), mitochondrial Na⁺-Ca²⁺ exchangers (mNCX; Palty et al., 2010), and the mitochondrial permeability transition pore (mPTP; Hunter et al., 1976). Besides the direct Ca²⁺ exchange, mitochondria may indirectly alter the cytosolic Ca²⁺ dynamics through many different ways under abnormal conditions (Yan et al., 2008; Florea and Blatter, 2010; Zhao et al., 2013; Xie et al., 2018), thus impacting Ca²⁺ waves and DADs. Under normal conditions, the occurrence of mitochondrial depolarizations through the mPTP opening is rare (Lu et al., 2016). However, the mPTP open probability increases in abnormal conditions, resulting in a higher degree of mitochondrial depolarization in the cell. The critical consequences that affect intracellular Ca²⁺ dynamics include an increased cytosolic reactive oxygen species (ROS) level, enhanced Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) activation via oxidative stress, Ca²⁺ influx into the cytosol from the mitochondria, and a decrease in the cytosolic ATP level, etc.

Reactive oxygen species can directly affect the RyRs hyperactivity and SERCA pump strength via its redox effect (Zima and Blatter, 2006; Wagner et al., 2013) or indirectly via oxidized CaMKII signaling (Xie et al., 2009; Foteinou et al., 2015). CaMKII activation is known to increase SERCA pump through phosphorylation of phospholamban (Hund and Rudy, 2004; Mattiazzi and Kranias, 2014), make RyRs leakier (Ai et al., 2005), and modulate LCCs and other membrane ionic currents (Anderson et al., 1994; Xiao et al., 1994; Yuan and Bers, 1994; Hund and Rudy, 2004; Hund et al., 2008). Furthermore, ATP depletion impairs the SERCA pump function (Sakamoto and Tonomura, 1980). Due to their highly complex interactions, it is difficult to dissect out the individual roles of mitochondrial Ca²⁺ exchange, ROS, ATP, and CaMKII activation in the genesis of DADs during mitochondrial depolarization by experiments. We have recently developed a whole-cell ventricular myocyte model consisting of a network of intermingled Ca²⁺ release units (CRUs) and mitochondria, which contains physiological details of mitochondrial membrane potential, mitochondrial Ca²⁺ cycling, mPTP stochastic opening and closing, intracellular ROS, and oxidized CaMKII signaling. Using this model, we have investigated the underlying mechanisms of Ca²⁺ alternans and early afterdepolarizations caused by mitochondrial depolarization and dissected each of the components (Xie et al., 2018; Song et al., 2019; Pandey et al., 2021).

We used this model to investigate the underlying mechanisms of spontaneous Ca²⁺ release mediated DADs under mitochondrial depolarization due to mPTP openings in the present work. Specifically, we performed computer simulations to reveal individual contributions of the components mentioned earlier to the genesis of Ca²⁺ waves and DADs. Our previous work provided mechanistic insights of generation of Ca²⁺ alternans under mitochondrial depolarization, and we have reported that the redox effect of ROS on RyRs and SERCA pump synergistically promote alternans (Pandey et al., 2021). Here, we show that the ROS redox regulation of RyRs plays a significant role in the genesis of Ca²⁺ waves and DADs during the acute phase of mitochondrial depolarization. Also, upregulation of MCU can promote DADs through Ca²⁺-dependent openings of mPTPs. However, the changes of oxidized CaMKIIs activation and the cytosolic ATP level are at much slower time scales than the redox effects of ROS, and thereby, they do not significantly impact the genesis of DADs in a relatively short duration after mitochondrial depolarization. Whereas, under chronic conditions, ATP depletion suppresses and enhanced CaMKII activation promotes the Ca²⁺ waves and DADs.

Materials and Methods

The details of the model, including the mathematical formulations and control values of the parameters, can be found in Song et al. (2019) and Pandey et al. (2021). Here, we describe some of the essential aspects of the model for the sake of this study.

The Overall Ventricular Myocyte Model Structure

Our rabbit ventricular myocyte model consists of a 3-dimensional coupled network of CRUs and mitochondria. This network contains 21504 (64 2812) CRUs and 5376 (64 14 6) mitochondria. The membrane potential (V) of the cell is described by

$$C_m\frac{dV}{dt}=I_{Na}+I_{Na,L}+I_{Ca,L}+I_{NCX}+I_{K1}+I_{Kr}+I_{Ks}+I_{to,f}+I_{to,s}+I_{NaK}+I_{K,ATP}+I_{Ca,b}-I_{sti}(1)$$

where C_m=1 mF/cm² is the cell membrane capacitance, and I_sti is the stimulus pulse with the current density being −80 mA/cm² and the duration being 0.5 ms.

Regulations of Reactive Oxygen Species and CaMKII on Ryanodine Receptors

The oxidized CaMKII activation and the redox effect of ROS both increase the open probability of RyRs (Wehrens et al., 2004; Ai et al., 2005; Guo et al., 2006; Zima and Blatter, 2006; Wagner et al., 2013). To incorporate these effects, the close-to-open rate (k₁₂) of RyRs was modeled as follows:

$$k_{12}=k_{base}{k}_u\left(1+\mathrm{△}{k}_{CaMKII}+\mathrm{△}{k}_{ROS}\right){\left({\left[C{a}^{2+}\right]}_p\right)}^2(2)$$

$$\mathrm{△}{k}_{CaMKII}=\frac{\mathrm{\Delta }{k}_{CaMk,max}}{1+{\left(\frac{\mathrm{\Delta }{k}_{CaMRyR}}{{[CamKII]}_{act}}\right)}^{h_{CaMKIIyR}}}(3)$$

$$\mathrm{△}{k}_{ROS}=\frac{\mathrm{\Delta }{K}_{ROS,max}}{1+{\left(\frac{\mathrm{\Delta }{k}_{mROSRyR}}{{[ROS]}_{cyt}}\right)}^{h_{ROSRyR}}}(4)$$

where △k_CaMKII and △k_ROS are the CaMKII-dependent (Eq. 3) and ROS-dependent components (Eq. 4), respectively. k_base and k_u are the rate constants. [Ca^{2 +}]_p is the Ca²⁺ concentration in the dyadic space of a CRU. [CaMKII]_act and [ROS]_cyt are the CaMKII activation level and the cytosolic ROS concentration in each CRU, respectively. The increase in CaMKII activation and ROS level increase k₁₂, which in turn increases the open probability of RyRs.

Regulations of Reactive Oxygen Species and CaMKII on Sarcoplasmic Reticulum Ca²⁺-ATPase Pump

The formulation of SERCA is

$${\text{J}}_{\text{up}}=v_{up}{f}_{up,ATP}{f}_{up,ROS}\frac{{\left[C{a}^{2+}\right]}_i^2}{{\left[C{a}^{2+}\right]}_i^2+{\left(K_i-PLB\left({\left[CaMKII\right]}_{act}\right)\right)}^2}(5)$$

where f_up,ATP, and f_up,ROS are ATP and ROS-dependent functions (Song et al., 2019):

$$f_{up,ATP}=\frac{1}{1+\frac{{[ADP]}_f}{k_{i,up}^{\prime }}+\left(1+\frac{{[ADP]}_f}{k_{i,up}}\right)\frac{k_{mup}ATP}{[ATP]}}(6)$$

$$f_{up,ROS}=\frac{1}{1+{\left(\frac{{\left[ROS\right]}_{cyt}}{k_{d,ros}}\right)}^{h_{ros,SERCA}}}+\frac{0.75}{1+{\left(\frac{k_{d,ros}}{{\left[ROS\right]}_{cyt}}\right)}^{h_{ros,SERCA}}}(7)$$

v_up is the maximum SERCA strength and K_ithe half-maximum value. PLB([CaMKII]_act) is a CaMKII dependent function. [CaMKII]_act is CaMKII activation level in the cytosolic space of a CRU.

The Mitochondrial Permeability Transition Pore Model

We used a 3-state (two close states C₀ and C₁, and an open state O) Markov model to simulate the stochastic opening and closing of the mPTP. The transition rate from the C₀ state to the C₁ state, k_c0c1, is set as:

$$k_{c0c1}={\mathrm{\alpha }}_0\left(1+199*\frac{{\left[C{a}^{2+}\right]}_m^{h_{mPTP}}}{{\left[C{a}^{2+}\right]}_m^{h_{mPTP}}+{\left[C{a}^{2+}\right]}_0^{h_{mPTP}}}\right)(8)$$

where h_mPTP is the Hill coefficient, [Ca^{2 +}]_m is the mitochondrial free Ca²⁺ in the corresponding mitochondrion, and [Ca^{2 +} ]₀ is the half-maximum value. We assume that other transition rates are constants. To simulate different levels of mPTP open probability, we multiplied a factor, α_mPTP, to the transition rate from C₁ to O, $k_{c10}^0$, i.e.,

$$k_{c1o}={\mathrm{\alpha }}_{mPTP}{k}_{c1o}^0(9)$$

Results

Mitochondrial Depolarization Due to Openings of Mitochondrial Permeability Transition Pore Promotes Spontaneous Ca²⁺ Release and Delayed Afterdepolarizations

We investigated the impact of mitochondrial depolarizations on the occurrence of Ca²⁺ waves and DADs via mPTP opening. We performed simulations over a wide range of α_mPTP values at the PCL of 300 ms (Figure 1A). α_mPTP is a factor multiplied to transition rate of mPTP opening, and increasing its value results in higher mPTP opening. For each simulation, the cell was paced for 140 beats (42 s), following 3 s without pacing in order to observe Ca²⁺ waves and DADs. As shown in Figure 1A, the amplitude of DAD increases with α_mPTP, suggesting that mitochondrial depolarization due to openings of mPTP promotes spontaneous Ca²⁺ release and DADs. Also, when α_mPTP is greater than ∼50, the proarrhythmic effect appears to saturate. Under the control condition (α_mPTP=1), there is no occurrence of DADs (Figures 1A,B). The cytosolic ROS is ∼2 μM, CaMKII activation is ∼0.2%, the cytosolic ATP is ∼5 mM, and most of the mitochondria remain repolarized ($-\overline{\mathrm{\Delta }\mathrm{\psi }}=$∼180 mV, and the mPTP open probability ∼0.8%, mitochondrial Ca²⁺ amplitude is ∼0.8 μM) (Figure 1B). However, with the higher open probability of mPTP (∼42%, for α_mPTP=60), we observed DADs. The corresponding line-scan image clearly shows enhanced spontaneous Ca²⁺ release as compared to a few scattered Ca²⁺ sparks under the control condition (Figure 1C). In this case, we should note that the mitochondrial Ca²⁺ amplitude increased to ∼1.2 M, and the cytosolic ROS drastically increased to ∼86 μM. Still, the CaMKII activation and the cytosolic ATP level insignificantly changed (∼0.6% and ∼4.8 mM, respectively).

FIGURE 1

Figure 1. Mitochondrial depolarization via mPTP opening promotes DADs. (A) The amplitude of DADs vs. α_mPTP. The resting potential is –86.0 mV. Note that for the range of α_mPTP in these simulations, we observed only one DAD after stopping pacing. (B) Time traces of V, ${\overline{\left[C{a}^{2+}\right]}}_i$, ${\overline{\left[C{a}^{2+}\right]}}_{JSR}$, ${\overline{\left[C{a}^{2+}\right]}}_m$, ${\overline{P}}_{mPTP}$, $-\overline{\mathrm{△}\mathrm{\psi }}$, ${\overline{\left[ROS\right]}}_{cyt}$,${\overline{\left[CaMKII\right]}}_{act}$, ${\overline{\left[ATP\right]}}_{cyt}$ for normal control (α_mPTP=1) in black and mitochondrial depolarization conditions (α_mPTP=60) in red, respectively. The pacing cycle length is 300 ms, and we stopped pacing after 140 beats (i.e., 42 s). This pacing protocol was used throughout the whole study. (C) Linescan images of the cytosolic Ca²⁺ for normal (top) and mitochondrial depolarization (bottom) conditions as in panel (B).

Role of Reactive Oxygen Species in the Genesis of Delayed Afterdepolarizations

The concentration of the cytosolic ROS depends on the level of mitochondrial depolarization, and thus, increasing α_mPTP increases the open probability of mPTP, which, in turn, elevates the level of cytosolic ROS. To further identify the role of cytosolic ROS in the genesis of DADs, we performed simulations for a clamped ROS level at the PCL 300 ms. In the free-running ROS case as shown in Figure 1, the ROS dynamics in the model remains intact, while in the clamped ROS condition, the cytosolic ROS is clamped to a constant regardless the level of mitochondrial depolarization is. Here we clamped ROS to be 1.0 μM, which is close to the level under the normal control condition. We then measured the amplitude of DAD with different α_mPTP values for the clamped ROS (Figure 2, red) condition. We observed that when the ROS was clamped at the control level (1.0 μM), no DADs occurred, suggesting that the cytosolic ROS plays a critical role in inducing DADs during mitochondrial depolarization.

FIGURE 2

Figure 2. The Redox effect of cytosolic ROS on RyR facilitates the genesis of DADs. The mean and standard deviation of the DAD amplitude vs. α_mPTP for three different conditions: clamped cytosolic ROS level at 1.0 μM (red circles), the redox effect of ROS only exerted on RyRs (blue squares), and the redox effect of ROS only exerted on SERCA (black triangles). The removal of the redox effect of ROS on RyRs and SERCA was executed by setting △k_ROS=0 in Eq. 2, and f_up,ROS=1 in Eq. 5, respectively. Ten random trials were performed for each given α_mPTP value.

Furthermore, the cytosolic ROS is known to alter the characteristics of both SERCA and RyRs (Zima and Blatter, 2006; Wagner et al., 2013). Therefore, we investigated the redox effect of ROS on the RyRs and SERCA pump separately to dissect out its individual role. We observed that DADs disappeared when we removed the redox effect of ROS on the RyRs (Figure 2, green). However, removing the redox effect on SERCA did not significantly impact the amplitude of DADs (Figure 2, blue), suggesting that the redox effect on SERCA may not play a critical role in the genesis of DADs. In fact, the direct redox effect of ROS inhibits the SERCA pump activity. Thus, removing the redox effect on SERCA increased the SR Ca²⁺ load instead, causing higher amplitudes of DADs. For instance, at α_mPTP=20, the amplitude of DAD is 7.17 mV when the redox effect of ROS is only exerted on RyRs (Figure 2, blue), but 6.84 mV under control as shown in Figure 1A.

Role of the Mitochondrial Ca²⁺ in the Genesis of the Delayed Afterdepolarizations

Several studies have shown that under certain pathological conditions, MCU activity is enhanced (Santulli et al., 2015; Xie et al., 2018), which may elevate the mitochondrial free Ca²⁺. Our previous study (Song et al., 2019) showed that the increase of MCU up to 20-fold does not alter cytosolic Ca²⁺ markedly at the steady-state. However, we hypothesize that the higher mitochondrial Ca²⁺ due to MCU overexpression could increase the mPTP open probability and cause higher ROS production in the cytosol (Korge et al., 2011). To test this hypothesis, we performed simulations to examine the effect of MCU overexpression on the genesis of DADs.

We multiplied a factor, denoted as α_MCU, to the maximal MCU conductance. α_MCU=1 represents the control case and higher α_MCU values are used to represent the different levels of MCU activity. We plotted in Figure 3A the amplitude of DADs for different α_MCU and α_mPTP. We observed that at α_mPTP = 1, increasing the MCU activity did not result in DADs even for α_MCU=50. When α_mPTP becomes greater, the effect of MCU on promoting DADs appears to be more significant. Time traces of membrane voltage, the whole-cell averaged cytosolic Ca²⁺ and the mitochondrial free Ca²⁺ for the three marked locations in the phase map (Figure 3A) are shown in Figure 3B. These results indicate that increasing MCU activity could promote spontaneous Ca²⁺ release and DADs. The mechanism revealed in the model is that increasing MCU activity elevates the mitochondrial free Ca²⁺, which increases the open probability of mPTP, resulting in the elevation of the cytosolic ROS, which in turn promotes the spontaneous Ca²⁺ release and DADs.

FIGURE 3

Figure 3. Mitochondrial Ca²⁺ Uniporter upregulation enhances the genesis of DADs through Ca²⁺- dependent opening of mPTP. (A) Dependence of the DAD amplitude on α_MCU and α_mPTP. (B) Time traces of V, ${\overline{\left[C{a}^{2+}\right]}}_i$, and ${\overline{\left[C{a}^{2+}\right]}}_m$ for the (α_mPTP, α_MCU) coordinates as marked in panel (A), following the same color codes. The zoomed-in sections of DAD and spontaneous Ca²⁺ oscillation are shown in the corresponding insets.

Role of Oxidized CaMKII Activation and ATP in the Genesis of Delayed Afterdepolarizations

As seen in Figure 1B, the CaMKII activation and ATP appeared to change slowly during the simulations due to the slow kinetics in the governing equations of their dynamics. It is computationally cumbersome to run long simulations (up to several thousand beats) for these variables to reach their steady states. Therefore, to evaluate the individual role of CaMKII activation and ATP in the genesis of DADs, we clamped them to different constant values, respectively.

Figure 4A shows the relationship between the amplitude of DADs and the CaMKII activation level. Our results clearly show that increasing CaMKII activation promotes DADs. However, due to its slow kinetics, CaMKII activation insignificantly changes during the acute phase of mitochondrial depolarization (Figure 1B), suggesting that CaMKII activation may not play a primary role in the genesis of DADs during the short period immediately after mitochondrial depolarization. Still, it may promote DADs chronically due to its regulation on SERCA, LCC, and RyRs (Wang et al., 2020).

FIGURE 4

Figure 4. Effect of CaMKII activation and ATP depletion on the genesis of DADs. The mean and standard deviation of the DAD amplitude under three clamped CaMKII activation (A) and ATP (B) levels, respectively. α_mPTP=30. Ten random trials were performed for each clamped CaMKII activation or ATP level.

Since the SERCA pump requires ATP for its normal function, a lower level of ATP directly impairs the SERCA pump activity (Eq. 5). However, ATP depletion is a slow process, which is evident from Figure 1B. Hence, we clamped ATP at different levels for α_mPTP = 30 from the beginning of the simulations. Figure 4B plots the relationship between the amplitude of DADs and the cytosolic ATP concentration, and it shows that ATP depletion suppresses DADs during mitochondrial depolarization (α_mPTP = 30). Our results demonstrate that the cytosolic ATP level has a significant impact on the genesis of DADs. However, the depletion of the cytosolic ATP concentration during mitochondrial depolarization is a relatively slow process. Therefore, ATP depletion should not play a central role in the genesis of DADs during the acute phase of mitochondrial depolarization. Thus, similar to CaMKII activation, ATP depletion may only suppress DADs in a much longer time scale.

Discussion

We used a physiological detailed ventricular myocyte model consisting of a 3D network of coupled CRUs and mitochondria to investigate the roles of mitochondrial depolarization via mPTP opening in the genesis of DADs. The systematic investigation of individual roles, including the cytosolic ROS, mitochondrial Ca²⁺, CaMKII activation, and the cytosolic ATP, reveals that the redox effect of ROS on RyRs may play an essential role in the occurrence of DADs during the acute phase of mitochondrial depolarization. Furthermore, increasing the MCU activity could promote DADs by increasing the mPTP open probability through mitochondrial Ca²⁺ dependent kinetics of mPTP. In addition, oxidized CaMKII activation promotes, and ATP depletion suppresses DADs chronically in the condition of mitochondrial depolarization.

Role of Reactive Oxygen Species, CaMKII Activation and ATP Depletion on the Genesis of Delayed Afterdepolarizations

Experimental findings reported that the mitochondria depolarization through the application of FCCP promotes Ca²⁺ waves (Zhao et al., 2013). Furthermore, the effects of FCCP were counteracted by the application of mPTP blocker cyclosporine A (Zhao et al., 2013). Also, our previous experimental observations using cyclophilin D knockout mouse model showed attenuation of Ca²⁺ waves (Gordan et al., 2016). Elevation of cytosolic ROS during mPTP opening could be a significant factor, and experimental evidence showed that oxidative stress during mitochondrial depolarization slightly alters SR Ca²⁺ leaks (Ca²⁺ spark) amplitude but drastically increases its frequency (Yan et al., 2008; Zhou et al., 2011; Williams et al., 2013). Furthermore, ROS can oxidize CaMKII and enhance its activation. It has been shown that ROS and CaMKII activation act on the major Ca²⁺ handling proteins, such as RyRs and SERCA (Hund and Rudy, 2004; Wehrens et al., 2004; Ai et al., 2005; Guo et al., 2006; Zima and Blatter, 2006; Wagner et al., 2013). The direct redox effect of ROS increases the leakiness of RyRs and decreases the strength of SERCA (Zima and Blatter, 2006; Wagner et al., 2013), whereas CaMKII activation increases both the leakiness of RyRs and the strength of SERCA (Hund and Rudy, 2004; Ai et al., 2005; Mattiazzi and Kranias, 2014). Thus, the observed Ca²⁺ waves in experiments during mitochondrial depolarization are the consequences of the combined effects of the above factors. However, what is the primary player remains unclear. Here, by using our previously established physiologically detailed computer model, we show that our simulation results agree with the experimental observations that mitochondrial depolarization could induce spontaneous Ca²⁺ release and DADs. Furthermore, we found that it is the redox effect on RyRs that causes the DADs under the acute effect of mitochondrial depolarization, and the redox effect of ROS on reducing SERCA strength may not be sufficient to suppress DADs (Figure 2).

In addition, we indeed observed an increase in CaMKII activation due to the mPTP opening (Figure 1B), but the dynamics of CaMKII activation appeared much slower than that of the ROS. The CaMKII activation was increased from ∼0.2% under control to ∼0.6% during the mitochondrial depolarization for a 45-sec long simulation (Figure 1B). A further simulation showed that CaMKII activation could reach up to ∼40% for a much more extended duration (1200 s). These results suggest that CaMKII activation may be too slow to play an essential role in inducing DADs under the acute effect of mitochondrial depolarization. However, simulations with different clamped CaMKII activation levels reveal that CaMKII activation may play a vital role in causing spontaneous Ca²⁺ release and DADs chronically, since a higher CaMKII activation level caused a greater DAD amplitude (Figure 4A).

Similarly, we have observed a slight change in cytosolic ATP under the acute effect of mitochondrial depolarizations in our simulations (Figure 1B). Therefore, within a relatively short duration after mitochondrial depolarization, ATP may not be involved in the genesis of spontaneous Ca²⁺ release and DADs. Clamped ATP simulations showed that a lower cytosolic ATP level is linked to a smaller amplitude of DADs (Figure 4B). This is because the reduction of ATP impaired SERCA activity and suppress the DADs, which agree well with experimental evidence that ATP synthase inhibitor, oligomycin, does not promote DADs (Zhao et al., 2013). Although ATP reduction seems to suppress Ca²⁺ waves and DADs, our simulations and other’s experimental work suggest that ATP reduction could promote Ca²⁺ alternans, which is still arrhythmogenic (Hüser et al., 2000; Kockskämper et al., 2005; Zima and Blatter, 2006; Pandey et al., 2021).

Mitochondrial Ca²⁺ Uniporter Overexpression and Delayed Afterdepolarizations

Mitochondrial Ca²⁺ uptake has been reported to rise in heart failure (Santulli et al., 2015; Xie et al., 2018) and can promote EADs (Xie et al., 2018) and Ca²⁺ alternans (Pandey et al., 2021). Our previous work demonstrated MCU upregulation could promote EADs in heart failure conditions without mPTP openings (Xie et al., 2018). And MCU upregulation promotes Ca²⁺ alternans through the Ca²⁺ dependent opening of mPTPs (Pandey et al., 2021). The previous experiment by Zhao et al. (2013) reported that the mitochondrial Ca²⁺ efflux in the proximity of the junctional SR played an essential role in the regulation of Ca²⁺ waves. Furthermore, our previous study has shown that MCU overexpression can lead to Ca²⁺ overload in mitochondria (Song et al., 2019). Also, there is evidence that mitochondrial Ca²⁺ overload can cause the Ca²⁺-dependent opening of mPTP (Kwong and Molkentin, 2015), resulting in mitochondrial depolarization (Santulli et al., 2015). Here, our simulation study shows that increasing MCU activity promotes spontaneous Ca²⁺ release and DADs (Figure 3). The underlying mechanism revealed in our simulations is that increasing MCU activity enhanced mitochondrial depolarization through the Ca²⁺-dependent openings of mPTP, which resulted in spontaneous Ca²⁺ release primarily due to the direct redox effect of ROS on RyRs.

Pathophysiological and Clinical Relevance

Mitochondrial dysfunction has been associated with increased arrhythmic risk (Santulli et al., 2015; Shimizu et al., 2015; Xie et al., 2015, 2018; Gordan et al., 2016), which could account for mortality in many cardiac diseases such as cardiomyopathy, heart failure, and ischemia/reperfusion injury (IRI). We have demonstrated that the direct redox effect of ROS on RyRs plays a critical role in promoting Ca²⁺ waves and DADs under the acute effect of mitochondrial depolarization. Furthermore, the upregulation of MCU can promote DADs through Ca²⁺-dependent opening of mPTPs. These findings suggest that pharmacological interventions targeted at avoiding ROS buildup and MCU upregulation may provide novel therapeutics to prevent or treat cardiac arrhythmias.

Limitations

This detailed model coupled AP, CRUs, and mitochondria to capture excitation-contraction-metabolism coupling in ventricular myocytes. However, it has some limitations. For instance, heterogeneities in T-tubule networks and distributions of ion channels and Ca²⁺ handling proteins are a few examples (Soeller and Cannell, 1999; Baddeley et al., 2009) that should be considered in the future analysis. These heterogeneities in T-tubule networks may alter the genesis of DADs (Song et al., 2018).

We note that in Figure 3A, for α_mPTP=1, there was virtually no DAD occurring even with α_MCU=50, suggesting that the proposed mechanism of increasing MCU activity inducing DADs in this study requires a certain basal level of mPTP opening. In this model, we consider the mPTP gating kinetics only mitochondrial Ca²⁺ dependent. However, the ROS-induced ROS release mechanism also impacts the mPTP open probability (Zorov et al., 2000, 2006), which is essential for modeling mitochondrial depolarization waves (Yang et al., 2010; Zhou et al., 2010; Nivala et al., 2011). Thus, this ROS-induced ROS release mechanism may provide another critical positive feedback loop between mitochondrial and cytosolic Ca²⁺ instability. In the future, the ROS-induced ROS release will be added to our model to study the role of mitochondrial depolarization waves in the genesis of arrhythmias in diseased conditions, such as heart failure.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Author Contributions

VP performed the simulations, analyzed simulation results, and drafted the manuscript. ZS prepared the figures. All the authors conceived and designed the study, interpreted the results, and edited and revised the manuscript.

Funding

This study is supported by the National Institutes of Health grant R01 HL139829.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Ai, X., Curran, J. W., Shannon, T. R., Bers, D. M., and Pogwizd, S. M. (2005). Ca2+/calmodulin-dependent protein kinase modulates cardiac ryanodine receptor phosphorylation and sarcoplasmic reticulum Ca2+ leak in heart failure. Circ. Res. 97, 1314–1322. doi: 10.1161/01.RES.0000194329.41863.89

CrossRef Full Text | Google Scholar

Anderson, M. E., Braun, A. P., Schulman, H., and Premack, B. A. (1994). Multifunctional Ca2+ / Calmodulin-dependent protein kinase mediates Ca2+ -induced enhancement of the L-type Ca2+ current in Rabbit ventricular myocytes. Circ. Res. 75, 854–861.

Google Scholar

Baddeley, D., Jayasinghe, I. D., Lam, L., Rossberger, S., Cannell, M. B., and Soeller, C. (2009). Optical single-channel resolution imaging of the ryanodine receptor distribution in rat cardiac myocytes. Proc. Natl. Acad. Sci. U.S.A. 106, 22275–22280. doi: 10.1073/pnas.0908971106

PubMed Abstract | CrossRef Full Text | Google Scholar

Baughman, J. M., Perocchi, F., Girgis, H. S., Plovanich, M., Belcher-Timme, C. A., Sancak, Y., et al. (2011). Integrative genomics identifies MCU as an essential component of the mitochondrial calcium uniporter. Nature 476, 341–345. doi: 10.1038/nature10234

PubMed Abstract | CrossRef Full Text | Google Scholar

Bers, D. M. (2002). Cardiac excitation-contraction coupling. Nature 415, 379–383. doi: 10.1016/B978-0-12-378630-2.00221-8

CrossRef Full Text | Google Scholar

Cheng, H., Lederer, M. R., Lederer, W. J., and Cannell, M. B. (1996). Calcium sparks and [Ca2+]i waves in cardiac myocytes. Am. J. Physiol. Cell Physiol. 270(1 Pt 1), C148–C159. doi: 10.1152/ajpcell.1996.270.1.c148

PubMed Abstract | CrossRef Full Text | Google Scholar

De Stefani, D., Raffaello, A., Teardo, E., Szabó, I., and Rizzuto, R. (2011). A forty-kilodalton protein of the inner membrane is the mitochondrial calcium uniporter. Nature 476, 336–340. doi: 10.1038/nature10230

PubMed Abstract | CrossRef Full Text | Google Scholar

Florea, S. M., and Blatter, L. A. (2010). The role of mitochondria for the regulation of cardiac alternans. Front. Physiol. 1:141. doi: 10.3389/fphys.2010.00141

PubMed Abstract | CrossRef Full Text | Google Scholar

Foteinou, P. T., Greenstein, J. L., and Winslow, R. L. (2015). Mechanistic Investigation of the Arrhythmogenic Role of Oxidized CaMKII in the Heart. Biophys. J. 109, 838–849. doi: 10.1016/j.bpj.2015.06.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Gordan, R., Fefelova, N., Jk, G., and Xie, L. H. (2016). Involvement of mitochondrial permeability transition pore (mPTP) in cardiac arrhythmias: evidence from Cyclophilin D knockout mice. Cell Calcium 60, 363–372. doi: 10.1016/j.ceca.2016.09.001.Involvement

CrossRef Full Text | Google Scholar

Guo, T., Zhang, T., Mestril, R., and Bers, D. M. (2006). Ca2+/calmodulin-dependent protein kinase II phosphorylation of ryanodine receptor does affect calcium sparks in mouse ventricular myocytes. Circ. Res. 99, 398–406. doi: 10.1161/01.RES.0000236756.06252.13

CrossRef Full Text | Google Scholar

Hoeker, G. S., Katra, R. P., Wilson, L. D., Plummer, B. N., and Laurita, K. R. (2009). Spontaneous calcium release in tissue from the failing canine heart. Am. J. Physiol. Hear. Circ. Physiol. 297, 10–15. doi: 10.1152/ajpheart.01320.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Hund, T. J., Decker, K. F., Kanter, E., Mohler, P. J., Boyden, P. A., Schuessler, R. B., et al. (2008). Role of activated CaMKII in abnormal calcium homeostasis and INa remodeling after myocardial infarction: insights from mathematical modeling. J. Mol. Cell. Cardiol. 45, 420–428. doi: 10.1016/j.yjmcc.2008.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Hund, T. J., and Rudy, Y. (2004). Rate dependence and regulation of action potential and calcium transient in a canine cardiac ventricular cell model. Circulation 110, 3168–3174. doi: 10.1161/01.CIR.0000147231.69595.D3

CrossRef Full Text | Google Scholar

Hunter, D. R., Haworth, R. A., and Southard, J. H. (1976). Relationship between configuration, function, and permeability in calcium treated mitochondria. J. Biol. Chem. 251, 5069–5077. doi: 10.1016/s0021-9258(17)33220-9

CrossRef Full Text | Google Scholar

Hüser, J., Wang, Y. G., Sheehan, K. A., Cifuentes, F., Lipsius, S. L., and Blatter, L. A. (2000). Functional coupling between glycolysis and excitation-contraction coupling underlies alternans in cat heart cells. J. Physiol. 524, 795–806. doi: 10.1111/j.1469-7793.2000.00795.x

PubMed Abstract | CrossRef Full Text | Google Scholar

January, C. T., and Fozzard, H. (1988). Delayed afterdepolrizatons in Heart muscle: mechanisms and relevance. Pharmacol. Rev. 40, 219–224.

Google Scholar

Katra, R. P., and Laurita, K. R. (2005). Cellular mechanism of calcium-mediated triggered activity in the heart. Circ. Res. 96, 535–542. doi: 10.1161/01.RES.0000159387.00749.3c

CrossRef Full Text | Google Scholar

Kockskämper, J., Zima, A. V., and Blatter, L. A. (2005). Modulation of sarcoplasmic reticulum Ca2+ release by glycolysis in cat atrial myocytes. J. Physiol. 564, 697–714. doi: 10.1113/jphysiol.2004.078782

PubMed Abstract | CrossRef Full Text | Google Scholar

Korge, P., Yang, L., Yang, J. H., Wang, Y., Qu, Z., and Weiss, J. N. (2011). Protective role of transient pore openings in calcium handling by cardiac mitochondria. J. Biol. Chem. 286, 34851–34857. doi: 10.1074/jbc.M111.239921

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwong, J. Q., and Molkentin, J. D. (2015). Physiological and pathological roles of the mitochondrial permeability transition pore in the Heart. Cell Metab. 21, 206–214. doi: 10.1016/j.cmet.2014.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, X., Kwong, J. Q., Molkentin, J. D., and Bers, D. M. (2016). Individual cardiac mitochondria undergo rare transient permeability transition pore openings. Circ. Res. 118, 834–841. doi: 10.1161/CIRCRESAHA.115.308093

PubMed Abstract | CrossRef Full Text | Google Scholar

Marban, E., Robinson, S. W., and Wier, W. G. (1986). Mechanisms of arrhythmogenic delayed and early afterdepolarizations in ferret ventricular muscle. J. Clin. Invest. 78, 1185–1192. doi: 10.1172/JCI112701

PubMed Abstract | CrossRef Full Text | Google Scholar

Mattiazzi, A., and Kranias, E. G. (2014). The role of CaMKII regulation of phospholamban activity in heart disease. Front. Pharmacol. 5:5. doi: 10.3389/fphar.2014.00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mohler, P. J., Schott, J. J., Gramolini, A. O., Dilly, K. W., Guatimosim, S., DuBell, W. H., et al. (2003). Ankyrin-B mutation causes type 4 long-QT cardiac arrhythmia and sudden cardiac death. Nature 421, 634–639. doi: 10.1038/nature01335

PubMed Abstract | CrossRef Full Text | Google Scholar

Nivala, M., Korge, P., Nivala, M., Weiss, J. N., and Qu, Z. (2011). Linking flickering to waves and whole-cell oscillations in a mitochondrial network model. Biophys. J. 101, 2102–2111. doi: 10.1016/j.bpj.2011.09.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Palty, R., Silverman, W. F., Hershfinkel, M., Caporale, T., Sensi, S. L., Parnis, J., et al. (2010). NCLX is an essential component of mitochondrial Na+/Ca2+ exchange. Proc. Natl. Acad. Sci. U.S.A. 107, 436–441. doi: 10.1073/pnas.0908099107

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, V., Xie, L., Id, Z. Q., and Id, Z. S. (2021). Mitochondrial depolarization promotes calcium alternans: mechanistic insights from a ventricular myocyte model. PLoS Comput. Biol. 17:e1008624. doi: 10.1371/journal.pcbi.1008624

PubMed Abstract | CrossRef Full Text | Google Scholar

Pogwizd, S. M., and Bers, D. M. (2003). Calcium cycling in heart failure. Respir. Circ. 51, 67–76.

Google Scholar

Qu, Z., Gang, H., Garfinkel, A., and Weiss, J. N. (2014). Nonlinear and stochastic dynamics in the Heart. Phys. Rep. 543, 61–162. doi: 10.1038/jid.2014.371

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosen, M. R., Moak, J. P., and Damiano, B. (1984). The clinical relevence of afterdepolrizations. Ann. N. Y. Acad. Sci. 427, 84–93.

Google Scholar

Ross, J. L., and Howlett, S. E. (2009). ß-Adrenoceptor stimulation exacerbates detrimental effects of ischemia and reperfusion in isolated guinea pig ventricular myocytes. Eur. J. Pharmacol. 602, 364–372. doi: 10.1016/j.ejphar.2008.11.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakamoto, J., and Tonomura, Y. (1980). Order of release of ADP and Pi from phosphoenzyme with bound ADP of Ca2+-dependent ATPase from sarcoplasmic reticulum and of Na+, K+-dependent ATPase studied by ADP-inhibition patterns. J. Biochem. 87, 1721–1727. doi: 10.1093/oxfordjournals.jbchem.a132916

PubMed Abstract | CrossRef Full Text | Google Scholar

Santulli, G., Xie, W., Reiken, S. R., and Marks, A. R. (2015). Mitochondrial calcium overload is a key determinant in heart failure. Proc. Natl. Acad. Sci. U.S.A. 112, 11389–11394. doi: 10.1073/pnas.1513047112

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimizu, H., Schredelseker, J., Huang, J., Lu, K., Naghdi, S., Lu, F., et al. (2015). Mitochondrial Ca 2 + uptake by the voltage-dependent anion channel 2 regulates cardiac rhythmicity. eLife 4, 1–20. doi: 10.7554/eLife.04801

PubMed Abstract | CrossRef Full Text | Google Scholar

Soeller, C., and Cannell, M. B. (1999). Examination of the transverse tubular system in living cardiac rat myocytes by 2-photon microscopy and digital image-processing techniques. Circ. Res. 84, 266–275. doi: 10.1161/01.RES.84.3.266

CrossRef Full Text | Google Scholar

Song, Z., Liu, M. B., and Qu, Z. (2018). Transverse tubular network structures in the genesis of intracellular calcium alternans and triggered activity in cardiac cells. J. Mol. Cell. Cardiol. 114, 288–299. doi: 10.1016/j.yjmcc.2017.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Z., Qu, Z., and Karma, A. (2017). Stochastic initiation and termination of calciummediated triggered activity in cardiac myocytes. Proc. Natl. Acad. Sci. U.S.A. 114, E270–E279. doi: 10.1073/pnas.1614051114

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Z., Xie, L. H., Weiss, J. N., and Qu, Z. (2019). A spatiotemporal ventricular myocyte model incorporating mitochondrial calcium cycling. Biophys. J. 117, 2349–2360. doi: 10.1016/j.bpj.2019.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, S., Rokita, A. G., Anderson, M. E., and Maier, L. S. (2013). Redox regulation of sodium and calcium handling. Antioxidants Redox Signal. 18, 1063–1077. doi: 10.1089/ars.2012.4818

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Shen, W., Zhang, S., Luo, G., Wang, K., Xu, Y., et al. (2020). The Role of CaMKII overexpression and oxidation in atrial fibrillation—A simulation study. Front. Physiol. 11:607809. doi: 10.3389/fphys.2020.607809

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, H., Chopra, N., Laver, D., Hwang, H. S., Davies, S. S., Roach, D. E., et al. (2009). Flecainide prevents catecholaminergic polymorphic ventricular tachycardia in mice and humans. Nat. Med. 15, 380–383. doi: 10.1038/nm.1942

PubMed Abstract | CrossRef Full Text | Google Scholar

Wehrens, X. H. T., Lehnart, S. E., Reiken, S. R., and Marks, A. R. (2004). Ca2+/calmodulin-dependent protein kinase II phosphorylation regulates the cardiac ryanodine receptor. Circ. Res. 94, e61–e70. doi: 10.1161/01.res.0000125626.33738.e2

CrossRef Full Text | Google Scholar

Williams, G. S. B., Boyman, L., Chikando, A. C., Khairallah, R. J., and Lederer, W. J. (2013). Mitochondrial calcium uptake. Proc. Natl. Acad. Sci. U.S.A. 110, 10479–10486. doi: 10.1073/pnas.1300410110

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, R. P., Cheng, H., Lederer, W. J., Suzuki, T., and Lakatta, E. G. (1994). Dual regulation of Ca2+/calmodulin-dependent kinase II activity by membrane voltage and by calcium influx. Proc. Natl. Acad. Sci. U.S.A. 91, 9659–9663. doi: 10.1073/pnas.91.20.9659

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, A., Song, Z., Liu, H., Zhou, A., Shi, G., Wang, Q., et al. (2018). Mitochondrial Ca2+ influx contributes to arrhythmic risk in nonischemic cardiomyopathy. J. Am. Heart Assoc. 7, 1–20. doi: 10.1161/JAHA.117.007805

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, L. H., Fuhua, C., Karagueuzian, H. S., and Weiss, J. N. (2009). Oxidative stress-induced afterdepolarizations and calmodulin kinase II signaling. Circ Res. 104, 79–86. doi: 10.1161/CIRCRESAHA.108.183475.Oxidative

CrossRef Full Text | Google Scholar

Xie, W., Santulli, G., Reiken, S. R., Yuan, Q., Osborne, B. W., Chen, B. X., et al. (2015). Mitochondrial oxidative stress promotes atrial fibrillation. Sci. Rep. 5, 1–11. doi: 10.1038/srep11427

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Y., Liu, J., Wei, C., Li, K., Xie, W., Wang, Y., et al. (2008). Bidirectional regulation of Ca2+ sparks by mitochondria-derived reactive oxygen species in cardiac myocytes. Cardiovasc. Res. 77, 432–441. doi: 10.1093/cvr/cvm047

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Korge, P., Weiss, J. N., and Qu, Z. (2010). Mitochondrial oscillations and waves in cardiac myocytes: insights from computational models. Biophys. J. 98, 1428–1438. doi: 10.1016/j.bpj.2009.12.4300

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, W., and Bers, D. M. (1994). Ca-dependent facilitation of cardiac Ca current is due to Ca-calmodulin-dependent protein kinase. Am. J. Physiol. 267(3 Pt 2), 982–993.

Google Scholar

Zhao, Z., Gordan, R., Wen, H., Fefelova, N., Zang, W. J., and Xie, L. H. (2013). Modulation of intracellular calcium waves and triggered activities by mitochondrial Ca flux in mouse cardiomyocytes. PLoS One 8:e0080574. doi: 10.1371/journal.pone.0080574

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Aon, M. A., Almas, T., Cortassa, S., Winslow, R. L., and O’Rourke, B. (2010). A reaction-diffusion model of ROS-induced ROS release in a mitochondrial network. PLoS Comput. Biol. 6:e1000657. doi: 10.1371/journal.pcbi.1000657

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Aon, M. A., Liu, T., and O’Rourke, B. (2011). Dynamic modulation of Ca2+ sparks by mitochondrial oscillations in isolated guinea pig cardiomyocytes under oxidative stress. J. Mol. Cell. Cardiol. 51, 632–639. doi: 10.1016/j.yjmcc.2011.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Zima, A. V., and Blatter, L. A. (2006). Redox regulation of cardiac calcium channels and transporters. Cardiovasc. Res. 71, 310–321. doi: 10.1016/j.cardiores.2006.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Zorov, D. B., Filburn, C. R., Klotz, L. O., Zweier, J. L., and Sollott, S. J. (2000). Reactive oxygen species (ROS)-induced ROS release: a new phenomenon accompanying induction of the mitochondrial permeability transition in cardiac myocytes. J. Exp. Med. 192, 1001–1014. doi: 10.1084/jem.192.7.1001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zorov, D. B., Juhaszova, M., and Sollott, S. J. (2006). Mitochondrial ROS-induced ROS release: an update and review. Biochim. Biophys. Acta Bioenerg. 1757, 509–517. doi: 10.1016/j.bbabio.2006.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar