Inhibiting Na+/K+ ATPase Can Impair Mitochondrial Energetics and Induce Abnormal Ca2+ Cycling and Automaticity in Guinea Pig Cardiomyocytes

Qince Li; Steven M. Pogwizd; Sumanth D. Prabhu; Lufang Zhou

doi:10.1371/journal.pone.0093928

Abstract

Cardiac glycosides have been used for the treatment of heart failure because of their capabilities of inhibiting Na⁺/K⁺ ATPase (NKA), which raises [Na⁺]_i and attenuates Ca²⁺ extrusion via the Na⁺/Ca²⁺ exchanger (NCX), causing [Ca²⁺]_i elevation. The resulting [Ca²⁺]_i accumulation further enhances Ca²⁺-induced Ca²⁺ release, generating the positive inotropic effect. However, cardiac glycosides have some toxic and side effects such as arrhythmogenesis, confining their extensive clinical applications. The mechanisms underlying the proarrhythmic effect of glycosides are not fully understood. Here we investigated the mechanisms by which glycosides could cause cardiac arrhythmias via impairing mitochondrial energetics using an integrative computational cardiomyocyte model. In the simulations, the effect of glycosides was mimicked by blocking NKA activity. Results showed that inhibiting NKA not only impaired mitochondrial Ca²⁺ retention (thus suppressed reactive oxygen species (ROS) scavenging) but also enhanced oxidative phosphorylation (thus increased ROS production) during the transition of increasing workload, causing oxidative stress. Moreover, concurrent blocking of mitochondrial Na⁺/Ca²⁺ exchanger, but not enhancing of Ca²⁺ uniporter, alleviated the adverse effects of NKA inhibition. Intriguingly, NKA inhibition elicited Ca²⁺ transient and action potential alternans under more stressed conditions such as severe ATP depletion, augmenting its proarrhythmic effect. This computational study provides new insights into the mechanisms underlying cardiac glycoside-induced arrhythmogenesis. The findings suggest that targeting both ion handling and mitochondria could be a very promising strategy to develop new glycoside-based therapies in the treatment of heart failure.

Citation: Li Q, Pogwizd SM, Prabhu SD, Zhou L (2014) Inhibiting Na⁺/K⁺ ATPase Can Impair Mitochondrial Energetics and Induce Abnormal Ca²⁺ Cycling and Automaticity in Guinea Pig Cardiomyocytes. PLoS ONE 9(4): e93928. https://doi.org/10.1371/journal.pone.0093928

Editor: Lai-Hua Xie, Rutgers-New Jersey Medical School, United States of America

Received: December 23, 2013; Accepted: March 11, 2014; Published: April 10, 2014

Copyright: © 2014 Li et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work is supported by NIH/NHLBI R00HL095648 and University of Alabama at Birmingham Comprehensive Cardiovascular Center. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Glycosides are capable of inhibiting sarcolemmal Na⁺/K⁺-ATPase (NKA), which blocks the extrusion of Na⁺ and results in cytosolic Na⁺ accumulation [1], [2]. Elevation of Na⁺ consequently suppresses Na⁺/Ca²⁺ exchanger (NCX), the primary Ca²⁺ efflux pathway in cardiac myocytes, leading to Ca²⁺ overload and increased sarcoplasmic reticulum (SR) Ca²⁺ uptake. The resultant greater Ca²⁺-induced Ca²⁺ release (CICR) allows for more powerful contractions by cross-bridge cycling in response to stimulations [3]. Because of their positive inotropic effects, cardiac glycosides, such as digoxin, have been widely used in the treatment of congestive heart failure [2], [4]. However, recently the use of glycoside treatment in HF patients has been largely supplanted by other drugs (e.g., angiotensin-converting enzyme (ACE) inhibitors, β-blockers and aldosterone antagonists) and cardiac resynchronization therapies [5], [6]. The diminished use of glycosides in the clinic was partially due to their well-known side effects such as cardiac arrhythmias, gastrointestinal symptoms, and central nervous system abnormalities [2], [5], [7], [8]. Whereas the proarrhythmic effect of glycosides largely confines their clinic applications, the detailed underlying molecular mechanisms are not completely understood.

A classical hypothesis on the proarrhythmic effect of glycosides is that when SR Ca²⁺ stores become too high, some Ca²⁺ might be released spontaneously through ryanodine receptors (RyRs), causing early or delayed afterdepolarizations or triggered activity [1], [8], [9]. However, emerging evidence suggests that SR Ca²⁺ overload is not the only root cause of glycoside-induced cardiac arrhythmias [5]. A series of work by Dr. O’Rourke’s group has suggested that glycosides (e.g. ouabain) may impair mitochondrial energy metabolism and raise oxidative stress in guinea pig cardiomyocytes, especially those at increased workload [1], [10], [11]. Particularly, their studies showed that ouabain-induced cytosolic Na⁺ accumulation caused mitochondrial Ca²⁺ deficiency, NADH imbalance, and increased reactive oxygen species (ROS) accumulation. Since many ion channels/exchangers underlying the action potential (AP) or involved in Ca²⁺ handling (e.g., the fast Na⁺ channels, RyRs, and SR Ca²⁺ ATPase) are redox and/or ATP sensitive, ouabain-induced mitochondrial dysfunction can disturb Ca²⁺ cycling and elicit erratic action potentials. Indeed, Liu et al. have demonstrated that ouabain caused mitochondrial oxidative stress and DADs in guinea pig ventricular myocytes [1]. They also showed that concurrent application of CGP-37157 (a mitochondrial Na⁺/Ca²⁺ exchanger, a.k.a. mNCE, inhibitor) with ouabain retained mitochondrial Ca²⁺ and NADH levels, suppressed ROS production and prevented ouabain-induced DADs [1], [10]. These findings highlight the important roles of mitochondrial Ca²⁺ and NADH homeostasis in glycoside-induced oxidative stress and cardiac arrhythmogenesis.

In cardiac cells, mitochondrial Ca²⁺ is regulated not only by mNCE but also by mitochondrial Ca²⁺ uniporters (MCU), the primary pathway of mitochondrial Ca²⁺ uptake. While the molecular identity of MCU has been revealed recently [12], [13], the mechanisms governing MCU Ca²⁺ uptake are still incompletely understood [14]–[16]. Consequently, the role of MCU in regulating Ca²⁺ cycling and excitation-contraction (E-C) coupling remains controversial. The major restriction of mitochondrial Ca²⁺ uptake is attributed to the low affinity of MCU to Ca²⁺. Specifically, the Ca²⁺ concentration for half-maximal mitochondrial Ca²⁺ uptake via MCU was reported as ∼10–20 μM in isolated mitochondria [17], [18], which substantially exceeds the cytosolic Ca²⁺ concentration during E-C coupling (0.5–1.5 μM). Therefore, it has been argued that mitochondria may have a minor role in regulating cellular Ca²⁺ cycling and E-C coupling [14], [19], [20]. On the contrary, other evidence suggests that mitochondria are functionally and physically tethered to SR by a mitochondrial fusion protein (namely mitofusin 2, a.k.a mfn2) [21]–[23]. The proximity between mitochondria and SR [24]–[28] can form a high Ca²⁺ microdomain, facilitating dynamic interorganellar coupling (e.g. rapid, beat-to-beat MCU Ca²⁺ uptake [10], [15], [29]) and the modulation of SR Ca²⁺ release by mitochondria [30], [31]. For instance, several studies have shown that the impairment of mitochondrial function, such as dissipation of mitochondrial membrane potential or inhibition of ATP synthesis, promotes proarrhythmic Ca²⁺ alternans [32]–[34]. How glycoside-induced mitochondrial dysfunction affects Ca²⁺ handling and how the altered Ca²⁺ cycling influences energy metabolism through feedback mechanisms are not well understood and difficult to examine experimentally due to the complex and nonlinear interactions formed by the participating subcellular components.

In this study, we incorporated a detailed antioxidant system into a recently published computational model of the cardiomyocyte to quantitatively investigate the mechanisms underlying cardiac glycoside-induced mitochondrial energetic dysfunction and glycoside- mediated abnormal Ca²⁺ cycling and automaticity. Our simulation work showed that NKA inhibition caused cytosolic Na⁺ and Ca²⁺ accumulations that adversely affected mitochondrial Ca²⁺ retention and energetics (such as reduced NADH and ATP production) similar to experimental observations. In addition, the impaired mitochondrial function could be significantly alleviated by inhibiting mNCE but barely ameliorated by enhancing MCU. Finally both Ca²⁺ and AP alternans could be elicited when severe ATP depletion occurred as the result of high frequency pacing and concurrent NKA inhibition, implicating a novel mechanism underlying the arrhythmogenic effect of cardiac glycosides.

Methods

Model Development

A multiscale guinea pig cardiomyocyte model was developed in order to systematically examine the effects of glycoside-induced high cytosolic Na⁺ ([Na⁺]_i) on the interaction between mitochondrial energetics and cellular electrophysiology. The model was based on our recently published ECME-RIRR model [35] and consisted of major ion channels/transporters underlying the action potential, E-C coupling and Ca²⁺ handling, as well as mitochondrial energetics, ROS-induced ROS release (RIRR) and intramitochondrial Ca²⁺ regulation. To replicate mitochondrial Ca²⁺ and Na⁺ handling, a refined mitochondrial model [36] consisting of Na⁺, H⁺ and P_i⁺ channels was incorporated. Moreover, a detailed antioxidant subsystem developed by Kembro et al. [37] was modified and added to account for the regulation of mitochondrial energetics (e.g. NADH) on ROS production and scavenging. Specifically, the expanded ROS metabolism pathway is comprised of: (i) glutathione peroxidase (GPX), Mn superoxide dismutase (MnSOD), glutaredoxin (GRX), and thioredoxin (Trx) subsystems in mitochondrial matrix; (ii) GPX, catalase (CAT), and Cu,Zn superoxide dismutase (Cu,ZnSOD) subsystems in the extra-mitochondrial space; and (iii) transhydrogenase (THD) on mitochondrial inner membrane. For simplicity, the pathway of isocitrate dehydrogenase 2 (IDH2) and diffusion of reduced glutathione (GSH) in the Kembro model [37] were not considered in the current work. The general model scheme is shown in Figure 1 and the complete model equations and parameters are listed in the Supplemental Information (File S1 and S2, respectively).

Download:

Figure 1. General scheme of guinea pig cardiomyocyte model.

The electrophysiological module describes major ion channels underlying the action potential, E-C coupling and Ca²⁺ handling. The mitochondrial module accounts for major components of mitochondrial energetics such as TCA cycle. The ROS-induced ROS release (RIRR) module describes ROS production, transport and scavenging. The mitochondrial energetics and RIRR are linked to Ca²⁺ handling and cellular electrophysiology through signaling ions (e.g., Ca²⁺ and Na⁺) and metabolic intermediates (e.g. ATP). In the figure, shunt represents the fraction of ROS production from the electron transport chain, and p1 and p2 represent the fraction of Ca²⁺ that MCU take up from the cytosol and the mitochondria-SR microdomain, respectively.

https://doi.org/10.1371/journal.pone.0093928.g001

Modeling Mitochondrial Ca²⁺ Uptake through MCU

It is still debated how mitochondria take up Ca²⁺ in the cardiomyocyte. Some studies suggest that mitochondria take up Ca²⁺ solely from the cytoplasm [20], while others propose that MCU channels are tightly coupled to the SR and thus the main uptake is from the SR-mitochondria Ca²⁺ microdomain [10], [16]. We hypothesized that MCU distribute non-uniformly along the mitochondrial membrane, with expression abundant in areas facing the microdomain and scarce on membrane away from the SR. Consequently, MCU were modeled to draw Ca²⁺ from both the cytosol and the microdomain. The kinetics formulas and parameters of Ca²⁺ uptake remained the same as those in the ECME-RIRR model [35], except (the maximal rate of MCU) was split as (representing uptake from the cytosol) and (representing uptake from the microdomain). For details see equations E168.1 and E168.2 in the Supplemental Information. The values of and were adjusted accordingly to ensure that the total steady-state mitochondrial Ca²⁺ uptake flux equals the flux in the ECME-RIRR model: , where V_uni, V_uni,1 and V_uni,2 are the flux of the total Ca²⁺ uptake, the Ca²⁺ uptake from the cytosol, and the Ca²⁺ uptake from the microdomain, respectively. p1 and p2 represent the percentages of Ca²⁺ uptake flux from each pool, and their values are determined by the distribution of MCU proteins on the mitochondrial membrane.

Model Simulation Protocol

The formulas of other processes, such as membrane potential, ion channels and metabolic reactions, and model parameters were the same as those in the ECME-RIRR model [35] (see the Supplemental Information). The code of the whole cell model was written in C⁺⁺ (Visual Studio; Microsoft, Redmond, WA). The nonlinear ordinary differential equations were integrated numerically with CVODE as described previously [35], [38].

The cardiomyocyte was stimulated at 0.25 Hz until the steady state was reached. The steady state values were then used as initial conditions for all simulations. In the model, the effect of cardiac glycosides was mimicked by inhibiting NKA activity by 50%, which caused a 2.4-fold increase in cytosolic Na⁺ concentration in 3 mins. This amount of Na⁺ increase is comparable to that caused by 1 μM ouabain in the experimental study [1]. The cell was paced at 0.25 Hz for 1 min and then at an increasing rate (e.g., 2 Hz or 4 Hz) for 3 mins with or without NKA inhibition. We first simulated the dynamics of mitochondrial Ca²⁺ uptake, energy production and ROS metabolism during the transition of increasing workload. Various p1:p2 ratios (ranging from 1∶0 to 0∶1) were evaluated to match the experimental data [1]. The defined p1:p2 ratio was then used in the subsequent simulations. We then examined the effect of NKA inhibition-induced mitochondrial dysfunction on Ca²⁺ transients and AP under low (i.e. the fraction of electron transport chain O₂⁻ production, shunt, = 1%) and relatively high (i.e. shunt = 2.5%) O₂⁻ production conditions. We also explored whether blocking mNCE or enhancing MCU alleviates NKA inhibition-induced mitochondrial energetics impairment. Finally we examined whether NKA inhibition and associated mitochondrial dysfunction are capable of triggering Ca²⁺ alternans in a paced cardiomyocyte. The simulation results were post-processed and plotted using Origin 8.6 (OriginLab, Northampton, MA).

Results

Effects of MCU Ca²⁺ Pools and NKA Inhibition on Intracellular Ca²⁺ Dynamics during the Transition of Increasing Workload

We first examined the impacts of different subcellular mitochondrial Ca²⁺ uptake pools (i.e. cytosol vs. mitochondria-SR microdomain or p1:p2 ratio) on Ca²⁺ dynamics during the transition of increasing workload with or without NKA inhibition. As expected, increasing pacing frequency (from 0.25 Hz to 2 Hz) caused significant increase in both cytosolic and mitochondrial [Ca²⁺] (Fig. 2A and 2C, respectively). Interestingly, these elevations were differentially influenced by the p1:p2 ratio. Particularly, when MCU were modeled to take up Ca²⁺ completely from the cytoplasmic pool (i.e. p1 = 1 & p2 = 0), [Ca²⁺]_i had the smallest elevation (e.g. systolic [Ca²⁺]_i increased by 62.5%) in response to the increasing workload (Fig. 2A). It is clear that shifting the portion of Ca²⁺ uptake from the cytosol to the microdomain (or reducing p1:p2 ratio) enhanced [Ca²⁺]_i accumulation. When p1 = 0 & p2 = 1 (i.e. completely from the microdomain), the increase of [Ca²⁺]_i was the greatest (Fig. 2A). The effect of p1:p2 ratio on mitochondrial Ca²⁺ ([Ca²⁺]_m) was similar to that of [Ca²⁺]_i, except when [Ca²⁺]_m was completely taken from the microdomain (i.e. p1 = 0 & p2 = 1), in which case [Ca²⁺]_m increased first and then decreased slightly (Fig. 2C).

Download:

Figure 2. Model simulated effects of MCU Ca²⁺ uptake pool on [Ca²⁺]_i and [Ca²⁺]_m under control (A and C) and NKA inhibition (B and D) conditions during the transition of increasing pacing (from 0.25 Hz to 2 Hz).

Different colors represent different p1:p2 ratios: red (1∶0); green (3∶1); black (1∶1); blue (1∶3); and pink (0∶1). (E). Summarized effect of p1:p2 ratio on [Ca²⁺]_m accumulation in response to increased energy demand. (F) Experimental data of the effect of NKA inhibition on [Ca²⁺]_m accumulation (Reproduced from [1] with permission). The pacing frequency was increased to 1 Hz at 60 s.

https://doi.org/10.1371/journal.pone.0093928.g002

NKA inhibition (by 50%) monotonically enhanced the accumulation of [Ca²⁺]_i and amplified the influence of MCU Ca²⁺ uptake pool (Fig. 2B). The effect of blocking NKA on pacing-induced [Ca²⁺]_m elevation was complex and dependent upon the p1:p2 ratio (Fig. 2D). As summarized in Figure 1E, NKA inhibition potentiated [Ca²⁺]_m accumulation when mitochondria were modeled to sequestrate all Ca²⁺ from the cytosol (i.e. p1 = 1). NKA inhibition started to suppress the pacing-induced [Ca²⁺]_m elevation with the pool shifting to the microdomain (i.e. reducing p1:p2 ratio). Specifically, when p1:p2 = 1∶3, NKA inhibition caused about 15% reduction in [Ca²⁺]_m accumulation, which was comparable to the experimental data (18% reduction in systolic [Ca²⁺]_m and 16% in diastolic [Ca²⁺]_m) (Fig. 2F, reproduced from [1]). Therefore, we proposed that mitochondria take up Ca²⁺ ions mainly, if not exclusively, from mitochondrial-SR microdomain. The p1:p2 ratio of 1∶3 was used in all subsequent simulations. It is worth pointing out that in the experiments [1] the cardiomyocytes were treated with isoproterenol in addition to 1 Hz pacing. Thus, in our simulations, a higher (2 Hz) pacing frequency was used to produce comparable energy demand. In addition, the percentage of NKA activity inhibition caused by ouabain was not measured in the experiment, so the comparison was more qualitative than quantitative to some extent.

Effect of NKA Inhibition on Intracellular Ion Homeostasis and Mitochondrial Energetics

Next, the effects of NKA inhibition on ionic dynamics and mitochondrial energetics during the transition of increasing workload were simulated under the physiological ROS production conditions (e.g. ROS production fraction, shunt = 1%). Increasing pacing frequency from 0.25 Hz to 2 Hz caused a 1.6-fold elevation in [Na⁺]_i and a 1.5-fold increase in mitochondrial Na⁺ concentration ([Na⁺]_m), respectively (Fig. 3A and 3B). The NADH concentration decreased transiently and then returned to the basal level (Fig. 3E), which is consistent with the experimental data when workload is acutely increased [10], [39]. Increasing pacing rate alone caused a slight increase in ROS (Fig. 3F) and had minor effects on mitochondrial membrane potential (ΔΨ_m) (Fig. 3G), [ATP]_i and [ATP]_m (Fig. S1). Inhibiting NKA led to further increases of [Na⁺]_i and [Na⁺]_m (by 2.4-fold and 2.1-fold, respectively, Fig. 3A and 3B) and enhanced [Ca²⁺]_i accumulation and [Ca²⁺]_m attenuation (Fig. 3C and 3D). The elevation of [Na⁺]_i was comparable to that caused by 1 μM ouabain in the experimental study [1]. The NADH level decreased by 9% and did not return to the basal level (Fig. 3E), along with a large (∼30- fold) increase of ROS (Fig. 3F) and small loss of ΔΨ_m (Fig. 3G) and ATP (Fig. S1). These simulations were comparable to experimental results (Inserts in Fig. 3; reproduced from [1] with permission).

Download:

Figure 3. Effects of blocking NKA on ion homeostasis and mitochondrial energetics under low O₂⁻ production conditions with or without concurrent mNCE inhibition.

The cell was paced at 0.25: control; Red: 50% NKA inhibition; Dark green: 50% NKA inhibition+60% mNCE inhibition. (A): [Na⁺]_i; (B): [Na⁺]_m; (C): [Ca²⁺]_i; (D): [Ca²⁺]_m; (E): NADH; (F): ROS; and (G): mitochondrial membrane potential (ΔΨ_m). In this simulation, shunt = 1% and p1:p2 = 1∶3. Inserts: Experimental data of effects of ouabain (red) and CGP-37157 (dark green) on diastolic [Ca²⁺]_m (Insert D), NADH (Insert E), and oxidative stress (Insert F) (Reproduced from [1]). In insert F, ouabain and CGP-37157 were applied at 1 min and pacing frequency was increased to 1 Hz at 5 min.

https://doi.org/10.1371/journal.pone.0093928.g003

As heart failure often associates with increased ROS production, we simulated the effect of NKA inhibition on mitochondrial energetics under a relatively higher ROS production condition (e.g. shunt = 2.5%). As shown in Figure 4, the small increase of shunt per se did not significantly alter ion concentrations and mitochondrial energetics under control condition. However, simultaneously blocking NKA (by 50%) and increasing energy demand (from 0.25 Hz to 2 Hz) caused detriment on mitochondrial energetics under the higher ROS production condition. Specifically, NKA inhibition triggered sustained mitochondrial oscillations, including NADH (Fig. 4E), ROS production (Fig. 4F) and ΔΨ_m (Fig. 4G). The oscillations also included Ca²⁺ and Na⁺ ions in both cytosol and mitochondria (Fig. 4A–4D).

Download:

Figure 4. Effects of blocking NKA on ion homeostasis and mitochondrial energetics under higher O₂⁻ production conditions with or without concurrent mNCE inhibition.

Black: control; Red: 50% NKA inhibition; Blue grey: 50% NKA inhibition+60% mNCE inhibition; Dark green: 50% NKA inhibition+90% mNCE inhibition. (A): [Na⁺]_i; (B): [Na⁺]_m; (C): [Ca²⁺]_i; (D): [Ca²⁺]_m; (E): NADH; (F): ROS; and (G): mitochondrial membrane potential (ΔΨ_m). In this simulation, shunt = 2.5% and p1:p2 = 1∶3.

https://doi.org/10.1371/journal.pone.0093928.g004

The Effect of Modulating mNCE or MCU on NKA Inhibition-induced Mitochondrial Dysfunction

In cardiomyocyte, mitochondrial energetics largely relies on mitochondrial Ca²⁺ retention and NADH homeostasis mediated by that Ca²⁺ retention, especially during the transition of increasing energy demand [38], [40], [41]. Therefore, we examined whether modulating mitochondrial Ca²⁺ handling channels, namely mNCE and MCU, could influence the effects of NKA inhibition on mitochondrial energetics. Our results showed that the concurrent inhibition of mNCE (e.g. by 60%) significantly ameliorated NKA inhibition-induced mitochondrial dysfunction, without abolishing the inotropic effect of NKA inhibition (Fig. 3). Specifically, NKA inhibition-induced [Na⁺]_i and [Na⁺]_m elevations were slightly attenuated and [Ca²⁺]_i accumulation was reduced by 13%. mNCE inhibition also overturned the blunted [Ca²⁺]_m retention and restored NADH concentration, as wells as suppressed ROS production and preserved ΔΨ_m (Fig. 3, dark green curves) and ATP (Fig. S1). The beneficial effect of mNCE inhibition was consistent with previous experimental findings that application of CGP, an mNCE inhibitor, prevented ouabain-induced mitochondrial dysfunction and resultant oxidative stress (Fig. 3 Inserts) [1], [10]. While moderate inhibition of mNCE was beneficial, too much mNCE blocking (e.g. by 90%) completely suppressed NKA inhibition-induced [Ca²⁺]_i elevation (and thus the inotropic effect) (data not shown). At the higher ROS production rate condition (e.g. shunt = 2.5%), blocking mNCE by 60% significantly delayed NKA inhibition-induced mitochondrial oscillations (Fig. 4, sky grey curves). Again, 90% mNCE inhibition counteracted the inotropic effect of glycoside, although it successfully prevented NKA inhibition-induced oxidative stress and mitochondrial oscillations (Fig. 4, dark green curves).

Another way to compensate NKA inhibition-induced mitochondrial Ca²⁺ decline is to increase mitochondrial Ca²⁺ uptake via MCU. Therefore, we examined whether enhancing MCU can alleviate NKA inhibition-induced mitochondrial energetics dysfunction similarly to mNCE inhibition. As expected, increasing MCU rate (e.g. by 100%) significantly increased [Ca²⁺]_m (Fig. 5D). However, its effects on NKA inhibition-induced [Na⁺]_i, [Na⁺]_m, and [Ca²⁺]_i alterations were not significant (Fig. 5A–5C). Enhancing MCU barely alleviated NKA inhibition-induced mitochondrial energetic impairments. For instance, NADH level slightly increased but did not completely recover (Fig. 5E), ROS still accumulated (Fig. 5F), and ΔΨ_m still dropped (Fig. 5G). Further enhancing MCU (e.g. by 4-fold) could increase mitochondrial Ca²⁺ accumulation to a level comparable to that caused by mNCE inhibition, but did not completely abolish NKA inhibition-induced mitochondrial energetics dysfunction (data not shown).

Download:

Figure 5. Effect of enhancing MCU on NKA inhibition-induced mitochondrial energetic dysfunction.

Black: control; Red: 50% NKA inhibition; Dark green: 50% NKA inhibition+100% MCU enhancement. (A): [Na⁺]_i; (B): [Na⁺]_m; (C): [Ca²⁺]_i; (D): [Ca²⁺]_m; (E): NADH; (F): ROS; and (G): mitochondrial membrane potential (ΔΨ_m). In this simulation, shunt = 1% and p1:p2 = 1∶3.

https://doi.org/10.1371/journal.pone.0093928.g005

In addition to blocking mNCE and enhancing MCU, we also examined the effects of concurrently enhancing mNCE or inhibiting MCU with NKA inhibition. Increasing mNCE (100%) slightly enhanced NKA inhibition-induced [Ca²⁺]_i elevation (Fig. S2A) and significantly exacerbated [Ca²⁺]_m decline (Fig. S2B), NADH reduction (Fig. S2C) and ROS elevation (Fig. S2D). Concurrently blocking MCU (60%) with NKA inhibition had similar effects except it didn’t substantially promote ROS elevation (Fig. S2D). Further blocking MCU (100%) induced a more profound ROS elevation (data not shown).

Effect of NKA Inhibition on Ca²⁺ Cycling and AP under more Stressed Conditions

Finally, we investigated the effect of NKA inhibition on mitochondrial energetics and cellular electrophysiology under more stressed conditions, such as more severe NKA inhibition and/or higher energy demand. Increasing pacing rate from 0.25 Hz to 4 Hz resulted in mitochondrial depolarization and oscillations when the shunt was set to 2.5% (Fig. 6A), as well as cyclic ATP depletions (Fig. 6B). Similar results have been reported in our previous study [35]. Concurrent induction of NKA inhibition by 50% prompted mitochondrial depolarization (Fig. 6A red curve) and accelerated ATP consumption, resulting in further decline of ATP (Fig. 6B red curve). Figures 6C–6F present the dynamics of SERCA Ca²⁺ uptake (J_up), SR Ca²⁺ content, Ca²⁺ transients and action potential (AP) at the last second (i.e. 249 s) of (the 4 Hz) simulations. It is clear that the depletion of ATP significantly impaired SERCA Ca²⁺ uptake (Fig. 6C) and caused imbalance of SR Ca²⁺ cycling, leading to reducing and alternating SR Ca²⁺ content (Fig. 6D), as well as Ca²⁺ transient alternans with significantly increased amplitude of the larger Ca²⁺ transients within the large-small Ca²⁺ alternations (Fig. 6E). AP alternans with elevated resting potential and suppressed amplitude was also observed (Fig. 6F). More severe NKA inhibition (e.g. by 90%) accelerated ATP depletion, resulting in earlier occurrence of alternans (data not shown). Concurrent inhibition of mNCE (by 90%) retarded mitochondrial oscillations and mitigated cytosolic ATP depletion (Fig. 6A and 6B, green curve), which consequently augmented SR Ca²⁺ uptake (Fig. 6C), leading to a higher SR Ca²⁺ concentration (Fig. 6D). The stabilization and restoration of SR Ca²⁺ content abolished Ca²⁺ and AP alternans, producing a smaller Ca²⁺ transient and a larger AP amplitude (Fig. 6E and 6F).

Download:

Figure 6. Action potential (AP) and Ca²⁺ alternans induced by NKA inhibition.

The cell was paced at 0.25(by 50%). AP and Ca²⁺ alternans were observed at the end of pacing. Black: control; Red: 50% NKA inhibition; Dark green: 50% NKA inhibition+90% mNCE inhibition. (A): mitochondrial membrane potential (ΔΨ_m); (B): ATP_i; (C): SERCA Ca²⁺ uptake; (D): [Ca²⁺]_SR; (E): [Ca²⁺]_i and (F): AP. In this simulation, shunt = 2.5% and p1:p2 = 1∶3.

https://doi.org/10.1371/journal.pone.0093928.g006

To further explore the development of NKA inhibition-induced Ca²⁺ and AP alternans, the dynamics of ATP, [Ca²⁺]_i, [Ca²⁺]_SR and AP before, during and after the 5^th and 9^th mitochondrial depolarization (indicated by arrows in Fig. 6A) were analyzed. Within the first 4 oscillations, the morphology of APs did not change significantly, although their amplitudes dropped gradually (data not shown). At the 5^th mitochondrial depolarization, [ATP]_i decreased dramatically from 4.8 mM to 2.1 mM (Fig. 7A). Consequently, SERCA Ca²⁺ uptake was significantly impaired and the SR Ca²⁺ release-uptake cycle was disturbed, resulting in a large [Ca²⁺]_SR depletion and [Ca²⁺]_i accumulation (Fig. 7C and 7E). The elevated intracellular Ca²⁺ stimulated Na⁺/Ca²⁺ exchanger and caused APD prolongation (Fig. 7G). As a result, the AP could not fully repolarize and was still in the refractory period when the next electrical stimulation fired, resulting in EAD-like abnormal APs (Fig. 7G). Nonetheless, when ΔΨ_m recovered and [ATP]_i rebounded to a relative higher level (Fig. 7A, dark green curve), the disturbed AP and Ca²⁺ transients could return to normal rhythms (Fig. 7E and 7G, dark green curves).

Download:

Figure 7. Comparisons of ATP_i (A and B), [Ca²⁺]_SR (C and D), [Ca²⁺]_i (E and F), and action potential (AP) (G and H) before (Black), during (Red) and after (Dark green) the 5^th (A,C,E,G) and 9^th (B,D,F,H) mitochondrial depolarization (marked with arrows in Fig. 6A).

In this simulation, shunt = 2.5%, p1:p2 = 1∶3, and NKA inhibition was 50%.

https://doi.org/10.1371/journal.pone.0093928.g007

As mitochondrial oscillations continued, ATP (Fig. 7B) and SR Ca²⁺ content further declined (Fig. 7D), the morphologies of Ca²⁺ transients and AP became more distorted. For instance, right before the 9^th depolarization, [ATP]_i and basal [Ca²⁺]_SR were already dramatically low (1.6 mM and 0.36 mM, respectively) (Fig. 7B and 7D, black curves). The subsequent ΔΨ_m depolarization caused extra decreases of [ATP]_i and SR Ca²⁺ (Fig. 7B and 7D, red curves) and a much larger Ca²⁺ transient that spanned several cycle lengths (Fig. 7F, red curve). The APs also displayed erratic morphologies, with a significantly prolonged AP followed by several ones with much smaller AP amplitudes (Fig. 7H, red curve). Interestingly, the Ca²⁺ transients and AP developed alternating behaviors (i.e., alternans) rather than recovering to normal rhythm (Fig. 7F and 7H, dark green curves) when ΔΨ_m repolarized and ATP restored (Fig. 7B).

Discussion

The present study provides, for the first time, a computational framework to quantitatively examine the mechanisms underlying mitochondrial dysfunction and cardiac arrhythmias induced by cardiac glycosides that are used in the treatment of heart failure. The main findings are: (i) NKA inhibition simultaneously raises [Ca²⁺]_i and blunts [Ca²⁺]_m retention, causing ATP depletion and ROS accumulation that disturb SR Ca²⁺ handling and elicit abnormal APs; (ii) NKA inhibition-induce adverse effects on mitochondrial energetics and Ca²⁺ cycling can be ameliorated by blocking mNCE, but not by enhancing MCU; (iii) under more severe conditions, NKA inhibition can cause proarrhythmic Ca²⁺ and AP alternans; and (iv) mitochondrial uniporter takes up Ca²⁺ mainly from the mitochondria-SR microdomain instead of the cytoplasm.

Glycosides and Mitochondrial Energetics Dysfunction: New Links Revealed by Model Simulations

In cardiomyocytes, mitochondria are the major sites of ROS production, and NADH plays a vital role in regulating both the production and scavenging of ROS. Previous studies showed that inhibiting mitochondrial Ca²⁺ uptake impaired NADH production and promoted ROS accumulation [1], [10]. Our simulations were able to reproduce these experimental observations. In addition, our data demonstrated for the first time that the blunted [Ca²⁺]_m uptake (and resultant reduced NADH production) was not the sole reason causing oxidative stress during NKA inhibition, because ROS accumulation could occur even when NADH did not decline significantly (Fig. 8). It is known that increasing pacing frequency not only increases energy consumption but also causes ADP accumulation, which directly stimulates oxidative phosphorylation and the associated ROS production. Since NKA inhibition exacerbates [Ca²⁺]_i overload and enhances ATP hydrolysis, the resulting further ADP_m elevation can activate ATP synthesis and respiration, resulting in additional ROS production. Therefore, we propose that NKA inhibition such as by cardiac glycosides causes mitochondrial oxidative stress through two separate but concomitant pathways: (i) NKA inhibition blunting [Ca²⁺]_m accumulation, which reduces NADH production and therefore ROS removal, and (ii) NKA inhibition increasing [Ca²⁺]_i and ATP hydrolysis, which produces a large amount of ADP that stimulates mitochondrial respiration and therefore ROS production.

Download:

Figure 8. Summary of [Ca²⁺]_i and mitochondrial energetics under 0.25 Hz, 2 Hz, 2 Hz + NKA inhibition, and 2 Hz + NKA inhibition + mNCE inhibition conditions.

(A): [Ca²⁺]_i, NADH and ROS_m; and (B): ADP_m, complex V activity (V_ATPase) and respiration rate (V_O2). The data was recorded at the end point of simulations (i.e. 3 mins after increasing pacing). In this simulation shunt = 1%, p1:p2 = 1∶3, and NKA inhibition was 50%.

https://doi.org/10.1371/journal.pone.0093928.g008

While the role of NKA inhibition in regulating mitochondrial ROS is well appreciated, the effect of ouabain on mitochondrial ATP synthesis remains elusive. A recent study showed that ouabain suppressed glucose-induced mitochondrial ATP production in rat pancreatic islets [42]. This inhibitory effect, caused by enhanced ROS production that impaired electron transport chain, could be observed only at high glucose concentration but not at the basal level of glucose. Similarly, another study reported that ouabain caused ATP reduction, a dramatic increase of intracellular Ca²⁺ (up to 40-fold) and enhanced ATP consumption in cultured chick ventricular cells [43]. The ATP reduction induced by ouabain was much less than that induced by cyanide and 2-deoxyglucosein, so the ATP loss observed in this experiment was likely caused by the increase of ATP consumption rather than the impairment of mitochondrial ATP synthesis. On the contrary, it has been shown that ouabain directly enhanced the rate of electron transport chain and the synthesis of ATP in cardiomyocytes via increasing the synthesis of ubiquinone [7]. These contradictory results might be due to differences in experimental conditions, animal species or tissue variations. Our results suggest that NKA inhibition differentially influence ATP production, depending on the severity of stress: (i) concurrent NKA inhibition and moderate increase of workload stimulates oxidative phosphorylation and ATP production, and (ii) under more severe conditions such as higher frequency pacing, NKA inhibition caused dramatic ROS accumulation, mitochondrial depolarization, and tremendous ATP depletion.

Glycosides-induced Ca²⁺ Alternans: A New Mitochondria-originated Arrhythmic Substrate

In previous animal and clinical studies, the arrhythmogenic effects of glycosides have been mainly attributed to NKA inhibition-induced Ca²⁺ overload and associated abnormal RyR Ca²⁺ release [2], [5], [7], [8], [44]. More recently, studies by Liu et al. suggested that mitochondrial dysfunction, particularly ouabain-induced ROS production, might also contribute to the arrhythmogenic effect of glycosides [1]. To better understand the mechanisms underlying glycoside-induced cardiac arrhythmias, we analyzed the dynamics of Ca²⁺ in the cytosol, mitochondria and SR, as well as the profiles of metabolic intermediates such as ATP and NADH. Our results showed that NKA inhibition caused a significant ROS accumulation during the transition of moderate workload increment. The mitochondria-derived ROS may lead to abnormal Ca²⁺ handling and erratic APs, as demonstrated in [1]. More importantly, the model revealed a novel phenomenon that has not reported before: under more stressed conditions such as higher frequency pacing, NKA inhibition causes Ca²⁺ and AP alternans by reducing ATP. It is worthy to point out that although in our simulations Ca²⁺ alternans could be induced only at higher pacing rate (i.e. 4 Hz), the mechanism underlying the genesis of this type of alternans is different from that of classical alternans induced by increasing pacing rate. Actually, in our hands pacing at 4 Hz per se did not induce Ca²⁺ alternans. Further analysis showed that the genesis of Ca²⁺ alternans was attributed to NKA inhibition-induced mitochondrial dysfunction. Particularly, severe ATP depletion impaired SERCA activity, reduced the SR Ca²⁺ content, and thereby generated SR release and uptake mismatch, causing intracellular Ca²⁺ alternans. A similar mitochondrial dysfunction-induced Ca²⁺ alternans has been reported by several experimental studies [45], [46]. The role of mitochondrial dysfunction in the genesis of Ca²⁺ alternans was further supported by the mNCE inhibition simulation. As shown in Figure 5, blocking mNCE ameliorated NKA inhibition-induced mitochondrial dysfunction, mitigating the decline of SR Ca²⁺ and suppressing the genesis of alternans. Thus, our in silico study revealed a new plausible mechanism responsible for glycoside-induced cardiac arrhythmias, which is different from the traditional paradigm of SR Ca²⁺ overload and associated abnormal RyR Ca²⁺ release. Our model also suggested that the increase of workload could potentiate the mismatch between energy supply and demand and increase the risk of NKA inhibition-mediated arrhythmias. Thus, digitalis toxic effects could be more problematic at increased heart rate.

Can Modulating Mitochondrial Ca²⁺ Handling Alleviate Glycoside-induced Arrhythmogenesis?

Studies have shown that CGP-17157, a mNCE inhibitor, mitigated ouabain-induced adverse effects on mitochondrial function and lowered the occurrence of cardiac arrhythmias in guinea pigs [1], [10]. The beneficial effect of CGP was thought due to its capability of retaining [Ca²⁺]_m accumulation that preserves mitochondrial NADH homeostasis and reduces ROS production. Our simulations are quantitatively comparable to these experimental results. Moreover, we showed that the mNCE inhibition attenuated cytosolic Ca²⁺ accumulation, which contributed to the anti-arrhythmic effect of CGP-17157. In another word, the therapeutic effect of mNCE inhibitors was attributed to not only preserving [Ca²⁺]_m retention and mitochondrial function, but also attenuating glycoside-induced intracellular Ca²⁺ overload. One should be aware though that it is the [Ca²⁺]_i elevation that is primarily responsible for the (beneficial) inotropic effect of glycosides treatment. Therefore, the dosage of the mNCE inhibitor is likely to be critical to counteracting the adverse effects of glycosides. Actually, we showed that too much mNCE inhibition, such as by 90%, completely abolished NKA inhibition-induced Ca²⁺ accumulation, although it significantly improved mitochondrial function and suppressed oxidative stress.

While inhibiting mNCE efficiently relieved NKA inhibition-induced adverse influence on mitochondrial energetics, the effect of increasing mitochondrial Ca²⁺ uniporters (MCU) activity was not evident, which might be caused by the substrate compartmentation of MCU. Our model analysis showed that MCU need to take up most of Ca²⁺ (>70%) from the microdomain in order to reproduce NKA inhibition-induced [Ca²⁺]_m attenuation observed in the experimental studies [1], [10]. These results indicate that the majority of Ca²⁺ taken up by MCU is from the microdomain instead of cytoplasm. Thus, enhancing MCU activity increased [Ca²⁺]_m uptake and retention, but it did not substantially reduce NKA inhibition-induced [Ca²⁺]_i overload and ROS production. These in silico results also imply that MCU are not uniformly distributed across the mitochondrial membrane, but rather are concentrated at mitochondria-SR contacting sites. The non-homogenous localized expression of ion channels on the sarcoplasmic or intracellular membranes has been reported previously. For instance, studies showed that L-type Ca²⁺ channels are concentrated in the neighborhood of the junctional SR at T-tubules [47], and this unique microstructure allows a correlation between ionic channel function and distribution. Similarly, abundant expression of MCU at the mitochondria-SR contacting sites can facilitate fast mitochondrial Ca²⁺ uptake in response to increased energy demand. With the molecular identity of MCU emerging [13], experimentally verifying the non-homogenous distribution of MCU in cardiac mitochondria is no longer far-fetched.

Study Limitations

The model developed in this study incorporates subcellular Ca²⁺ compartmentation for MCU, therefore providing advantages over existing cardiomyocyte models [48]–[51]. However, it also has some limitations that need be addressed. (i) The current model lacks an isoproterenol signaling pathway, which impedes the direct comparison between model simulations and experimental data of Liu et al. [1]. However, this limitation would not alter our conclusion on the NKA inhibition-induced mitochondrial dysfunction. (ii) Previous studies have suggested that Ca²⁺ wave propagation plays a critical role in intracellular Ca²⁺ regulation and genesis of Ca²⁺ alternans [45], [46]. The “mitochondria-SR Ca²⁺ microdomain” notion proposed in this study (and by others elsewhere) may be involved in regulating Ca²⁺ wave propagation. Incorporating spatial resolution of local Ca²⁺ distribution in the microdomain would allow examination of its role in regulating cellular electrophysiology and mitochondrial energetics. (iii) The current model did not incorporate the direct link between mitochondrial ROS and redox sensitive SR Ca²⁺ proteins. It has been shown that oxidative stress can modulate SR Ca²⁺ handling proteins and cause abnormal electrical activities in cardiomyocytes [52]–[54]. The NKA inhibition-induced ROS production would exacerbate glycoside-induced abnormal cellular electrophysiology [1]. Adding mdROS modulation of redox sensitive ion channels/transporters is ongoing in our laboratory. (iv) In the current model, the equations of complexes of the electron transport chain are lumped, which makes it hard to investigate the specific role of the complexes of electron transport chain (ETC) in mitochondrial dysfunction. A more detailed ETC model [55] can be incorporated into the ECME-RIRR model to examine the effect of ETC complexes on mitochondrial function. (v) The present model did not consider the effect of energy depletion on ATP-sensitive potassium channels (K_ATP), as our main focus is on mitochondrial dysfunction and Ca²⁺ handling. It is expected that incorporating K_ATP would exacerbate the deleterious effect of NKA inhibition-induced mitochondrial dysfunction on cardiac action potentials. As shown before, depletion of ATP activates K_ATP currents, resulting in AP shortening or even cellular inexcitability [35], [56].

In summary, the present study provides a novel computational tool to quantitatively investigate the proarrhythmic effects of glycosides (and resultant NKA inhibition) in the cardiomyocyte. The simulations reveal that NKA inhibition can cause mitochondrial oxidative stress through two separate and concomitant pathways, which involve an intricate interplay between cardiac electrophysiology and mitochondrial energetics. Moreover, the model suggests that under more severe conditions, NKA inhibition can induce Ca²⁺ alternans, contributing to its proarrhythmic effect. Finally, the study verifies the role of mitochondria-SR Ca²⁺ microdomain in regulating mitochondrial Ca²⁺ uptake and energetics. Taken together, this work underscores the importance of targeting both intracellular ion handling and mitochondrial energetics in developing new glycoside-based heart failure therapies.

Supporting Information

Figure S1.

Effects of blocking NKA on cytosolic and mitochondrial ATP under low O₂⁻ production conditions with or without concurrent mNCE inhibition. The cell was paced at 0.25 Hz for 1 min then at 2 Hz for 3 mins. Black: control; Red: 50% NKA inhibition; Dark green: 50% NKA inhibition+60% mNCE inhibition. (A): [ATP]_i; (B): [ATP]_m. shunt = 1% and p1:p2 = 1∶3.

https://doi.org/10.1371/journal.pone.0093928.s001

(TIF)

Figure S2.

Effect of blocking MCU or enhancing mNCE on NKA inhibition-induced Ca²⁺ and mitochondrial energetics changes. Black: control; Red: 50% NKA inhibition; Dark green: 50% NKA inhibition+60% MCU inhibition; Blue: 50% NKA inhibition+100% mNCE enhancement. (A): [Ca²⁺]_i; (B): [Ca²⁺]_m; (C): NADH; (D): ROS. shunt = 1% and p1:p2 = 1∶3.

https://doi.org/10.1371/journal.pone.0093928.s002

(TIF)

File S1.

Supplemental Information model equations.

https://doi.org/10.1371/journal.pone.0093928.s003

(DOCX)

File S2.

Supplemental Information model parameters.

https://doi.org/10.1371/journal.pone.0093928.s004

(DOCX)

Acknowledgments

The authors thank Dr. X. Margaret Liu for reading and insightful criticisms.

Author Contributions

Conceived and designed the experiments: QL SMP SDP LZ. Performed the experiments: QL. Analyzed the data: QL SMP LZ. Contributed reagents/materials/analysis tools: QL LZ. Wrote the paper: QL SMP SDP LZ.

References

1. Liu T, Brown DA, O’Rourke B (2010) Role of mitochondrial dysfunction in cardiac glycoside toxicity. J Mol Cell Cardiol 49: 728–736.
- View Article
- Google Scholar
2. Rahimtoola SH (2004) Digitalis therapy for patients in clinical heart failure. Circulation 109: 2942–2946.
- View Article
- Google Scholar
3. Bers DM, Despa S (2009) Na+ transport in cardiac myocytes; Implications for excitation-contraction coupling. IUBMB Life 61: 215–221.
- View Article
- Google Scholar
4. Hashim T, Elbaz S, Patel K, Morgan CJ, Fonarow GC, et al.. (2013) Digoxin and 30-day All-cause Hospital Admission in Older Patients with Chronic Diastolic Heart Failure. The American journal of medicine.
5. Gonano LA, Sepulveda M, Rico Y, Kaetzel M, Valverde CA, et al. (2011) Calcium-calmodulin kinase II mediates digitalis-induced arrhythmias. Circ Arrhythm Electrophysiol 4: 947–957.
- View Article
- Google Scholar
6. Rathore SS, Curtis JP, Wang Y, Bristow MR, Krumholz HM (2003) Association of serum digoxin concentration and outcomes in patients with heart failure. JAMA 289: 871–878.
- View Article
- Google Scholar
7. Campia I, Sala V, Kopecka J, Leo C, Mitro N, et al. (2012) Digoxin and ouabain induce the efflux of cholesterol via liver X receptor signalling and the synthesis of ATP in cardiomyocytes. Biochem J 447: 301–311.
- View Article
- Google Scholar
8. Ferrier GR (1977) Digitalis arrhythmias: role of oscillatory afterpotentials. Prog Cardiovasc Dis 19: 459–474.
- View Article
- Google Scholar
9. Roberts BN, Christini DJ (2011) NHE Inhibition Does Not Improve Na+ or Ca2+ Overload During Reperfusion: Using Modeling to Illuminate the Mechanisms Underlying a Therapeutic Failure. Plos Computational Biology 7.
10. Maack C, Cortassa S, Aon MA, Ganesan AN, Liu T, et al. (2006) Elevated cytosolic Na+ decreases mitochondrial Ca2+ uptake during excitation-contraction coupling and impairs energetic adaptation in cardiac myocytes. Circ Res 99: 172–182.
- View Article
- Google Scholar
11. Liu T, O’Rourke B (2008) Enhancing mitochondrial Ca2+ uptake in myocytes from failing hearts restores energy supply and demand matching. Circ Res 103: 279–288.
- View Article
- Google Scholar
12. De Stefani D, Raffaello A, Teardo E, Szabo I, Rizzuto R (2011) A forty-kilodalton protein of the inner membrane is the mitochondrial calcium uniporter. Nature 476: 336–340.
- View Article
- Google Scholar
13. Baughman JM, Perocchi F, Girgis HS, Plovanich M, Belcher-Timme CA, et al. (2011) Integrative genomics identifies MCU as an essential component of the mitochondrial calcium uniporter. Nature 476: 341–345.
- View Article
- Google Scholar
14. Dedkova EN, Blatter LA (2008) Mitochondrial Ca2+ and the heart. Cell Calcium 44: 77–91.
- View Article
- Google Scholar
15. Liu T, O’Rourke B (2009) Regulation of mitochondrial Ca2+ and its effects on energetics and redox balance in normal and failing heart. J Bioenerg Biomembr 41: 127–132.
- View Article
- Google Scholar
16. O’Rourke B, Blatter LA (2009) Mitochondrial Ca2+ uptake: tortoise or hare? J Mol Cell Cardiol 46: 767–774.
- View Article
- Google Scholar
17. Denton RM, McCormack JG (1985) Ca2+ transport by mammalian mitochondria and its role in hormone action. Am J Physiol 249: E543–554.
- View Article
- Google Scholar
18. Huser J, Blatter LA, Sheu SS (2000) Mitochondrial calcium in heart cells: beat-to-beat oscillations or slow integration of cytosolic transients? J Bioenerg Biomembr 32: 27–33.
- View Article
- Google Scholar
19. Bassani JW, Bassani RA, Bers DM (1994) Relaxation in rabbit and rat cardiac cells: species-dependent differences in cellular mechanisms. J Physiol 476: 279–293.
- View Article
- Google Scholar
20. Bassani RA, Bassani JW, Bers DM (1992) Mitochondrial and sarcolemmal Ca2+ transport reduce [Ca2+]i during caffeine contractures in rabbit cardiac myocytes. J Physiol 453: 591–608.
- View Article
- Google Scholar
21. Chen Y, Csordas G, Jowdy C, Schneider TG, Csordas N, et al. (2012) Mitofusin 2-containing mitochondrial-reticular microdomains direct rapid cardiomyocyte bioenergetic responses via interorganelle Ca2+ crosstalk. Circ Res 111: 863–875.
- View Article
- Google Scholar
22. Csordas G, Hajnoczky G (2009) SR/ER-mitochondrial local communication: calcium and ROS. Biochim Biophys Acta 1787: 1352–1362.
- View Article
- Google Scholar
23. Csordás G, Hajnóczky G (2009) SR/ER–mitochondrial local communication: Calcium and ROS. Biochimica et Biophysica Acta (BBA) - Bioenergetics 1787: 1352–1362.
- View Article
- Google Scholar
24. Sharma VK, Ramesh V, Franzini-Armstrong C, Sheu SS (2000) Transport of Ca2+ from sarcoplasmic reticulum to mitochondria in rat ventricular myocytes. J Bioenerg Biomembr 32: 97–104.
- View Article
- Google Scholar
25. Berridge MJ (2006) Calcium microdomains: organization and function. Cell Calcium 40: 405–412.
- View Article
- Google Scholar
26. Franzini-Armstrong C (2007) ER-mitochondria communication. How privileged? Physiology (Bethesda) 22: 261–268.
- View Article
- Google Scholar
27. Rizzuto R, Pozzan T (2006) Microdomains of intracellular Ca2+: molecular determinants and functional consequences. Physiol Rev 86: 369–408.
- View Article
- Google Scholar
28. Ramesh V, Sharma VK, Sheu SS, Franzini-Armstrong C (1998) Structural proximity of mitochondria to calcium release units in rat ventricular myocardium may suggest a role in Ca2+ sequestration. Ann NY Acad Sci 853: 341–344.
- View Article
- Google Scholar
29. Rizzuto R, Duchen MR, Pozzan T (2004) Flirting in little space: the ER/mitochondria Ca2+ liaison. Sci STKE 2004: re1.
- View Article
- Google Scholar
30. Yan Y, Liu J, Wei C, Li K, Xie W, et al. (2008) Bidirectional regulation of Ca2+ sparks by mitochondria-derived reactive oxygen species in cardiac myocytes. Cardiovasc Res 77: 432–441.
- View Article
- Google Scholar
31. Zhou L, Aon MA, Liu T, O’Rourke B (2011) Dynamic modulation of Ca(2+) sparks by mitochondrial oscillations in isolated guinea pig cardiomyocytes under oxidative stress. Journal of molecular and cellular cardiology 51: 632–639.
- View Article
- Google Scholar
32. Florea SM, Blatter LA (2010) The role of mitochondria for the regulation of cardiac alternans. Front Physiol 1: 141.
- View Article
- Google Scholar
33. Florea SM, Blatter LA (2012) Regulation of cardiac alternans by beta-adrenergic signaling pathways. Am J Physiol Heart Circ Physiol 303: H1047–1056.
- View Article
- Google Scholar
34. Smith RM, Visweswaran R, Talkachova I, Wothe JK, Tolkacheva EG (2013) Uncoupling the mitochondria facilitates alternans formation in the isolated rabbit heart. Am J Physiol Heart Circ Physiol 305: H9–H18.
- View Article
- Google Scholar
35. Zhou L, Cortassa S, Wei AC, Aon MA, Winslow RL, et al. (2009) Modeling cardiac action potential shortening driven by oxidative stress-induced mitochondrial oscillations in guinea pig cardiomyocytes. Biophys J 97: 1843–1852.
- View Article
- Google Scholar
36. Wei AC, Aon MA, O’Rourke B, Winslow RL, Cortassa S (2011) Mitochondrial energetics, pH regulation, and ion dynamics: a computational-experimental approach. Biophys J 100: 2894–2903.
- View Article
- Google Scholar
37. Kembro JM, Aon MA, Winslow RL, O’Rourke B, Cortassa S (2013) Integrating mitochondrial energetics, redox and ROS metabolic networks: a two-compartment model. Biophys J 104: 332–343.
- View Article
- Google Scholar
38. Cortassa S, Aon MA, O’Rourke B, Jacques R, Tseng HJ, et al. (2006) A computational model integrating electrophysiology, contraction, and mitochondrial bioenergetics in the ventricular myocyte. Biophys J 91: 1564–1589.
- View Article
- Google Scholar
39. Brandes R, Bers DM (1999) Analysis of the mechanisms of mitochondrial NADH regulation in cardiac trabeculae. Biophys J 77: 1666–1682.
- View Article
- Google Scholar
40. Brandes R, Bers DM (1997) Intracellular Ca2+ increases the mitochondrial NADH concentration during elevated work in intact cardiac muscle. Circ Res 80: 82–87.
- View Article
- Google Scholar
41. Brandes R, Bers DM (2002) Simultaneous measurements of mitochondrial NADH and Ca(2+) during increased work in intact rat heart trabeculae. Biophys J 83: 587–604.
- View Article
- Google Scholar
42. Kajikawa M, Tsuura Y, Mukai E, Takeda T, Hamamoto Y, et al. (2002) Ouabain Suppresses Glucose-Induced Mitochondrial ATP Production and Insulin Release by Generating Reactive Oxygen Species in Pancreatic Islets. Diabetes 51: 2522–2529.
- View Article
- Google Scholar
43. Ishida H, Kohmoto O, Bridge JH, Barry WH (1988) Alterations in cation homeostasis in cultured chick ventricular cells during and after recovery from adenosine triphosphate depletion. J Clin Invest 81: 1173–1181.
- View Article
- Google Scholar
44. Fedida D, Noble D, Rankin AC, Spindler AJ (1987) The arrhythmogenic transient inward current iTI and related contraction in isolated guinea-pig ventricular myocytes. J Physiol 392: 523–542.
- View Article
- Google Scholar
45. Li Q, O’Neill SC, Tao T, Li Y, Eisner D, et al. (2012) Mechanisms by which cytoplasmic calcium wave propagation and alternans are generated in cardiac atrial myocytes lacking T-tubules-insights from a simulation study. Biophys J 102: 1471–1482.
- View Article
- Google Scholar
46. Mackenzie L, Roderick HL, Berridge MJ, Conway SJ, Bootman MD (2004) The spatial pattern of atrial cardiomyocyte calcium signalling modulates contraction. J Cell Sci 117: 6327–6337.
- View Article
- Google Scholar
47. Bers DM (2001) Excitation–Contraction Coupling and Cardiac Contractile Force edn 2 Kluwer Academic, Dordrecht, Netherlands.
48. Magnus G, Keizer J (1998) Model of β-cell mitochondrial calcium handling and electrical activity. I.Cytoplasmic variables. Am J Physiol Cell Physiol 274: C1158–1173.
- View Article
- Google Scholar
49. Matsuoka S, Sarai N, Jo H, Noma A (2004) Simulation of ATP metabolism in cardiac excitation-contraction coupling. Prog Biophys Mol Biol 85: 279–299.
- View Article
- Google Scholar
50. Greenstein JL, Hinch R, Winslow RL (2006) Mechanisms of excitation-contraction coupling in an integrative model of the cardiac ventricular myocyte. Biophys J 90: 77–91.
- View Article
- Google Scholar
51. Gauthier LD, Greenstein JL, Winslow RL (2012) Toward an integrative computational model of the Guinea pig cardiac myocyte. Front Physiol 3: 244.
- View Article
- Google Scholar
52. Xie LH, Chen F, Karagueuzian HS, Weiss JN (2009) Oxidative-stress-induced afterdepolarizations and calmodulin kinase II signaling. Circ Res 104: 79–86.
- View Article
- Google Scholar
53. Barrington PL, Meier CF Jr, Weglicki WB (1988) Abnormal electrical activity induced by free radical generating systems in isolated cardiocytes. J Mol Cell Cardiol 20: 1163–1178.
- View Article
- Google Scholar
54. Horackova M, Ponka P, Byczko Z (2000) The antioxidant effects of a novel iron chelator salicylaldehyde isonicotinoyl hydrazone in the prevention of H(2)O(2) injury in adult cardiomyocytes. Cardiovasc Res 47: 529–536.
- View Article
- Google Scholar
55. Gauthier LD, Greenstein JL, Cortassa S, O’Rourke B, Winslow RL (2013) A computational model of reactive oxygen species and redox balance in cardiac mitochondria. Biophys J 105: 1045–1056.
- View Article
- Google Scholar
56. Akar FG, Aon MA, Tomaselli GF, O’Rourke B (2005) The mitochondrial origin of postischemic arrhythmias. J Clin Invest 115: 3527–3535.
- View Article
- Google Scholar

[ref1] 1. Liu T, Brown DA, O’Rourke B (2010) Role of mitochondrial dysfunction in cardiac glycoside toxicity. J Mol Cell Cardiol 49: 728–736.
View Article
Google Scholar

[2] View Article

[3] Google Scholar

[ref2] 2. Rahimtoola SH (2004) Digitalis therapy for patients in clinical heart failure. Circulation 109: 2942–2946.
View Article
Google Scholar

[5] View Article

[6] Google Scholar

[ref3] 3. Bers DM, Despa S (2009) Na+ transport in cardiac myocytes; Implications for excitation-contraction coupling. IUBMB Life 61: 215–221.
View Article
Google Scholar

[8] View Article

[9] Google Scholar

[ref4] 4. Hashim T, Elbaz S, Patel K, Morgan CJ, Fonarow GC, et al.. (2013) Digoxin and 30-day All-cause Hospital Admission in Older Patients with Chronic Diastolic Heart Failure. The American journal of medicine.

[ref5] 5. Gonano LA, Sepulveda M, Rico Y, Kaetzel M, Valverde CA, et al. (2011) Calcium-calmodulin kinase II mediates digitalis-induced arrhythmias. Circ Arrhythm Electrophysiol 4: 947–957.
View Article
Google Scholar

[12] View Article

[13] Google Scholar

[ref6] 6. Rathore SS, Curtis JP, Wang Y, Bristow MR, Krumholz HM (2003) Association of serum digoxin concentration and outcomes in patients with heart failure. JAMA 289: 871–878.
View Article
Google Scholar

[15] View Article

[16] Google Scholar

[ref7] 7. Campia I, Sala V, Kopecka J, Leo C, Mitro N, et al. (2012) Digoxin and ouabain induce the efflux of cholesterol via liver X receptor signalling and the synthesis of ATP in cardiomyocytes. Biochem J 447: 301–311.
View Article
Google Scholar

[18] View Article

[19] Google Scholar

[ref8] 8. Ferrier GR (1977) Digitalis arrhythmias: role of oscillatory afterpotentials. Prog Cardiovasc Dis 19: 459–474.
View Article
Google Scholar

[21] View Article

[22] Google Scholar

[ref9] 9. Roberts BN, Christini DJ (2011) NHE Inhibition Does Not Improve Na+ or Ca2+ Overload During Reperfusion: Using Modeling to Illuminate the Mechanisms Underlying a Therapeutic Failure. Plos Computational Biology 7.

[ref10] 10. Maack C, Cortassa S, Aon MA, Ganesan AN, Liu T, et al. (2006) Elevated cytosolic Na+ decreases mitochondrial Ca2+ uptake during excitation-contraction coupling and impairs energetic adaptation in cardiac myocytes. Circ Res 99: 172–182.
View Article
Google Scholar

[25] View Article

[26] Google Scholar

[ref11] 11. Liu T, O’Rourke B (2008) Enhancing mitochondrial Ca2+ uptake in myocytes from failing hearts restores energy supply and demand matching. Circ Res 103: 279–288.
View Article
Google Scholar

[28] View Article

[29] Google Scholar

[ref12] 12. De Stefani D, Raffaello A, Teardo E, Szabo I, Rizzuto R (2011) A forty-kilodalton protein of the inner membrane is the mitochondrial calcium uniporter. Nature 476: 336–340.
View Article
Google Scholar

[31] View Article

[32] Google Scholar

[ref13] 13. Baughman JM, Perocchi F, Girgis HS, Plovanich M, Belcher-Timme CA, et al. (2011) Integrative genomics identifies MCU as an essential component of the mitochondrial calcium uniporter. Nature 476: 341–345.
View Article
Google Scholar

[34] View Article

[35] Google Scholar

[ref14] 14. Dedkova EN, Blatter LA (2008) Mitochondrial Ca2+ and the heart. Cell Calcium 44: 77–91.
View Article
Google Scholar

[37] View Article

[38] Google Scholar

[ref15] 15. Liu T, O’Rourke B (2009) Regulation of mitochondrial Ca2+ and its effects on energetics and redox balance in normal and failing heart. J Bioenerg Biomembr 41: 127–132.
View Article
Google Scholar

[40] View Article

[41] Google Scholar

[ref16] 16. O’Rourke B, Blatter LA (2009) Mitochondrial Ca2+ uptake: tortoise or hare? J Mol Cell Cardiol 46: 767–774.
View Article
Google Scholar

[43] View Article

[44] Google Scholar

[ref17] 17. Denton RM, McCormack JG (1985) Ca2+ transport by mammalian mitochondria and its role in hormone action. Am J Physiol 249: E543–554.
View Article
Google Scholar

[46] View Article

[47] Google Scholar

[ref18] 18. Huser J, Blatter LA, Sheu SS (2000) Mitochondrial calcium in heart cells: beat-to-beat oscillations or slow integration of cytosolic transients? J Bioenerg Biomembr 32: 27–33.
View Article
Google Scholar

[49] View Article

[50] Google Scholar

[ref19] 19. Bassani JW, Bassani RA, Bers DM (1994) Relaxation in rabbit and rat cardiac cells: species-dependent differences in cellular mechanisms. J Physiol 476: 279–293.
View Article
Google Scholar

[52] View Article

[53] Google Scholar

[ref20] 20. Bassani RA, Bassani JW, Bers DM (1992) Mitochondrial and sarcolemmal Ca2+ transport reduce [Ca2+]i during caffeine contractures in rabbit cardiac myocytes. J Physiol 453: 591–608.
View Article
Google Scholar

[55] View Article

[56] Google Scholar

[ref21] 21. Chen Y, Csordas G, Jowdy C, Schneider TG, Csordas N, et al. (2012) Mitofusin 2-containing mitochondrial-reticular microdomains direct rapid cardiomyocyte bioenergetic responses via interorganelle Ca2+ crosstalk. Circ Res 111: 863–875.
View Article
Google Scholar

[58] View Article

[59] Google Scholar

[ref22] 22. Csordas G, Hajnoczky G (2009) SR/ER-mitochondrial local communication: calcium and ROS. Biochim Biophys Acta 1787: 1352–1362.
View Article
Google Scholar

[61] View Article

[62] Google Scholar

[ref23] 23. Csordás G, Hajnóczky G (2009) SR/ER–mitochondrial local communication: Calcium and ROS. Biochimica et Biophysica Acta (BBA) - Bioenergetics 1787: 1352–1362.
View Article
Google Scholar

[64] View Article

[65] Google Scholar

[ref24] 24. Sharma VK, Ramesh V, Franzini-Armstrong C, Sheu SS (2000) Transport of Ca2+ from sarcoplasmic reticulum to mitochondria in rat ventricular myocytes. J Bioenerg Biomembr 32: 97–104.
View Article
Google Scholar

[67] View Article

[68] Google Scholar

[ref25] 25. Berridge MJ (2006) Calcium microdomains: organization and function. Cell Calcium 40: 405–412.
View Article
Google Scholar

[70] View Article

[71] Google Scholar

[ref26] 26. Franzini-Armstrong C (2007) ER-mitochondria communication. How privileged? Physiology (Bethesda) 22: 261–268.
View Article
Google Scholar

[73] View Article

[74] Google Scholar

[ref27] 27. Rizzuto R, Pozzan T (2006) Microdomains of intracellular Ca2+: molecular determinants and functional consequences. Physiol Rev 86: 369–408.
View Article
Google Scholar

[76] View Article

[77] Google Scholar

[ref28] 28. Ramesh V, Sharma VK, Sheu SS, Franzini-Armstrong C (1998) Structural proximity of mitochondria to calcium release units in rat ventricular myocardium may suggest a role in Ca2+ sequestration. Ann NY Acad Sci 853: 341–344.
View Article
Google Scholar

[79] View Article

[80] Google Scholar

[ref29] 29. Rizzuto R, Duchen MR, Pozzan T (2004) Flirting in little space: the ER/mitochondria Ca2+ liaison. Sci STKE 2004: re1.
View Article
Google Scholar

[82] View Article

[83] Google Scholar

[ref30] 30. Yan Y, Liu J, Wei C, Li K, Xie W, et al. (2008) Bidirectional regulation of Ca2+ sparks by mitochondria-derived reactive oxygen species in cardiac myocytes. Cardiovasc Res 77: 432–441.
View Article
Google Scholar

[85] View Article

[86] Google Scholar

[ref31] 31. Zhou L, Aon MA, Liu T, O’Rourke B (2011) Dynamic modulation of Ca(2+) sparks by mitochondrial oscillations in isolated guinea pig cardiomyocytes under oxidative stress. Journal of molecular and cellular cardiology 51: 632–639.
View Article
Google Scholar

[88] View Article

[89] Google Scholar

[ref32] 32. Florea SM, Blatter LA (2010) The role of mitochondria for the regulation of cardiac alternans. Front Physiol 1: 141.
View Article
Google Scholar

[91] View Article

[92] Google Scholar

[ref33] 33. Florea SM, Blatter LA (2012) Regulation of cardiac alternans by beta-adrenergic signaling pathways. Am J Physiol Heart Circ Physiol 303: H1047–1056.
View Article
Google Scholar

[94] View Article

[95] Google Scholar

[ref34] 34. Smith RM, Visweswaran R, Talkachova I, Wothe JK, Tolkacheva EG (2013) Uncoupling the mitochondria facilitates alternans formation in the isolated rabbit heart. Am J Physiol Heart Circ Physiol 305: H9–H18.
View Article
Google Scholar

[97] View Article

[98] Google Scholar

[ref35] 35. Zhou L, Cortassa S, Wei AC, Aon MA, Winslow RL, et al. (2009) Modeling cardiac action potential shortening driven by oxidative stress-induced mitochondrial oscillations in guinea pig cardiomyocytes. Biophys J 97: 1843–1852.
View Article
Google Scholar

[100] View Article

[101] Google Scholar

[ref36] 36. Wei AC, Aon MA, O’Rourke B, Winslow RL, Cortassa S (2011) Mitochondrial energetics, pH regulation, and ion dynamics: a computational-experimental approach. Biophys J 100: 2894–2903.
View Article
Google Scholar

[103] View Article

[104] Google Scholar

[ref37] 37. Kembro JM, Aon MA, Winslow RL, O’Rourke B, Cortassa S (2013) Integrating mitochondrial energetics, redox and ROS metabolic networks: a two-compartment model. Biophys J 104: 332–343.
View Article
Google Scholar

[106] View Article

[107] Google Scholar

[ref38] 38. Cortassa S, Aon MA, O’Rourke B, Jacques R, Tseng HJ, et al. (2006) A computational model integrating electrophysiology, contraction, and mitochondrial bioenergetics in the ventricular myocyte. Biophys J 91: 1564–1589.
View Article
Google Scholar

[109] View Article

[110] Google Scholar

[ref39] 39. Brandes R, Bers DM (1999) Analysis of the mechanisms of mitochondrial NADH regulation in cardiac trabeculae. Biophys J 77: 1666–1682.
View Article
Google Scholar

[112] View Article

[113] Google Scholar

[ref40] 40. Brandes R, Bers DM (1997) Intracellular Ca2+ increases the mitochondrial NADH concentration during elevated work in intact cardiac muscle. Circ Res 80: 82–87.
View Article
Google Scholar

[115] View Article

[116] Google Scholar

[ref41] 41. Brandes R, Bers DM (2002) Simultaneous measurements of mitochondrial NADH and Ca(2+) during increased work in intact rat heart trabeculae. Biophys J 83: 587–604.
View Article
Google Scholar

[118] View Article

[119] Google Scholar

[ref42] 42. Kajikawa M, Tsuura Y, Mukai E, Takeda T, Hamamoto Y, et al. (2002) Ouabain Suppresses Glucose-Induced Mitochondrial ATP Production and Insulin Release by Generating Reactive Oxygen Species in Pancreatic Islets. Diabetes 51: 2522–2529.
View Article
Google Scholar

[121] View Article

[122] Google Scholar

[ref43] 43. Ishida H, Kohmoto O, Bridge JH, Barry WH (1988) Alterations in cation homeostasis in cultured chick ventricular cells during and after recovery from adenosine triphosphate depletion. J Clin Invest 81: 1173–1181.
View Article
Google Scholar

[124] View Article

[125] Google Scholar

[ref44] 44. Fedida D, Noble D, Rankin AC, Spindler AJ (1987) The arrhythmogenic transient inward current iTI and related contraction in isolated guinea-pig ventricular myocytes. J Physiol 392: 523–542.
View Article
Google Scholar

[127] View Article

[128] Google Scholar

[ref45] 45. Li Q, O’Neill SC, Tao T, Li Y, Eisner D, et al. (2012) Mechanisms by which cytoplasmic calcium wave propagation and alternans are generated in cardiac atrial myocytes lacking T-tubules-insights from a simulation study. Biophys J 102: 1471–1482.
View Article
Google Scholar

[130] View Article

[131] Google Scholar

[ref46] 46. Mackenzie L, Roderick HL, Berridge MJ, Conway SJ, Bootman MD (2004) The spatial pattern of atrial cardiomyocyte calcium signalling modulates contraction. J Cell Sci 117: 6327–6337.
View Article
Google Scholar

[133] View Article

[134] Google Scholar

[ref47] 47. Bers DM (2001) Excitation–Contraction Coupling and Cardiac Contractile Force edn 2 Kluwer Academic, Dordrecht, Netherlands.

[ref48] 48. Magnus G, Keizer J (1998) Model of β-cell mitochondrial calcium handling and electrical activity. I.Cytoplasmic variables. Am J Physiol Cell Physiol 274: C1158–1173.
View Article
Google Scholar

[137] View Article

[138] Google Scholar

[ref49] 49. Matsuoka S, Sarai N, Jo H, Noma A (2004) Simulation of ATP metabolism in cardiac excitation-contraction coupling. Prog Biophys Mol Biol 85: 279–299.
View Article
Google Scholar

[140] View Article

[141] Google Scholar

[ref50] 50. Greenstein JL, Hinch R, Winslow RL (2006) Mechanisms of excitation-contraction coupling in an integrative model of the cardiac ventricular myocyte. Biophys J 90: 77–91.
View Article
Google Scholar

[143] View Article

[144] Google Scholar

[ref51] 51. Gauthier LD, Greenstein JL, Winslow RL (2012) Toward an integrative computational model of the Guinea pig cardiac myocyte. Front Physiol 3: 244.
View Article
Google Scholar

[146] View Article

[147] Google Scholar

[ref52] 52. Xie LH, Chen F, Karagueuzian HS, Weiss JN (2009) Oxidative-stress-induced afterdepolarizations and calmodulin kinase II signaling. Circ Res 104: 79–86.
View Article
Google Scholar

[149] View Article

[150] Google Scholar

[ref53] 53. Barrington PL, Meier CF Jr, Weglicki WB (1988) Abnormal electrical activity induced by free radical generating systems in isolated cardiocytes. J Mol Cell Cardiol 20: 1163–1178.
View Article
Google Scholar

[152] View Article

[153] Google Scholar

[ref54] 54. Horackova M, Ponka P, Byczko Z (2000) The antioxidant effects of a novel iron chelator salicylaldehyde isonicotinoyl hydrazone in the prevention of H(2)O(2) injury in adult cardiomyocytes. Cardiovasc Res 47: 529–536.
View Article
Google Scholar

[155] View Article

[156] Google Scholar

[ref55] 55. Gauthier LD, Greenstein JL, Cortassa S, O’Rourke B, Winslow RL (2013) A computational model of reactive oxygen species and redox balance in cardiac mitochondria. Biophys J 105: 1045–1056.
View Article
Google Scholar

[158] View Article

[159] Google Scholar

[ref56] 56. Akar FG, Aon MA, Tomaselli GF, O’Rourke B (2005) The mitochondrial origin of postischemic arrhythmias. J Clin Invest 115: 3527–3535.
View Article
Google Scholar

[161] View Article

[162] Google Scholar

Figures

Abstract

Introduction

Methods

Model Development

Modeling Mitochondrial Ca2+ Uptake through MCU

Model Simulation Protocol

Results

Effects of MCU Ca2+ Pools and NKA Inhibition on Intracellular Ca2+ Dynamics during the Transition of Increasing Workload

Effect of NKA Inhibition on Intracellular Ion Homeostasis and Mitochondrial Energetics

The Effect of Modulating mNCE or MCU on NKA Inhibition-induced Mitochondrial Dysfunction

Effect of NKA Inhibition on Ca2+ Cycling and AP under more Stressed Conditions

Discussion

Glycosides and Mitochondrial Energetics Dysfunction: New Links Revealed by Model Simulations

Glycosides-induced Ca2+ Alternans: A New Mitochondria-originated Arrhythmic Substrate

Can Modulating Mitochondrial Ca2+ Handling Alleviate Glycoside-induced Arrhythmogenesis?

Study Limitations

Supporting Information

Figure S1.

Figure S2.

File S1.

File S2.

Acknowledgments

Author Contributions

References

Modeling Mitochondrial Ca²⁺ Uptake through MCU

Effects of MCU Ca²⁺ Pools and NKA Inhibition on Intracellular Ca²⁺ Dynamics during the Transition of Increasing Workload

Effect of NKA Inhibition on Ca²⁺ Cycling and AP under more Stressed Conditions

Glycosides-induced Ca²⁺ Alternans: A New Mitochondria-originated Arrhythmic Substrate

Can Modulating Mitochondrial Ca²⁺ Handling Alleviate Glycoside-induced Arrhythmogenesis?