BioMed Research International
Volume 2020 (2020), Article ID 1597012, 13 pages
https://doi.org/10.1155/2020/1597012

Research Article

Role of Oxidation-Dependent CaMKII Activation in the Genesis of Abnormal Action Potentials in Atrial Cardiomyocytes: A Simulation Study

Na Zhao,¹ Qince Li,1,2 Haibo Sui,¹ and Henggui Zhang1,2,3

¹School of Computer Science and Technology, Harbin Institute of Technology, Harbin 150000, China
²Peng Cheng Laboratory, Shenzhen 518000, China
³School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK

Correspondence should be addressed to Qince Li; qinceli@hit.edu.cn and Henggui Zhang; henggui.zhang@manchester.ac.uk

Received 25 October 2019; Revised 20 May 2020; Accepted 2 June 2020; Published 22 June 2020

Guest Editor: Roland E. Akhigbe

Copyright © 2020 Na Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Atrial fibrillation is a common cardiac arrhythmia with an increasing incidence rate. Particularly for the aging population, understanding the underlying mechanisms of atrial arrhythmia is important in designing clinical treatment. Recently, experiments have shown that atrial arrhythmia is associated with oxidative stress. In this study, an atrial cell model including oxidative-dependent Ca²⁺/calmodulin- (CaM-) dependent protein kinase II (CaMKII) activation was developed to explore the intrinsic mechanisms of atrial arrhythmia induced by oxidative stress. The simulation results showed that oxidative stress caused early afterdepolarizations (EADs) of action potentials by altering the dynamics of transmembrane currents and intracellular calcium cycling. Oxidative stress gradually elevated the concentration of calcium ions in the cytoplasm by enhancing the L-type Ca²⁺ current and sarcoplasmic reticulum (SR) calcium release. Owing to increased intracellular calcium concentration, the inward Na⁺/Ca²⁺ exchange current was elevated which slowed down the repolarization of the action potential. Thus, the action potential was prolonged and the L-type Ca²⁺ current was reactivated, resulting in the genesis of EAD. Furthermore, based on the atrial single-cell model, a two-dimensional (2D) ideal tissue model was developed to explore the effect of oxidative stress on the electrical excitation wave conduction in 2D tissue. Simulation results demonstrated that, under oxidative stress conditions, EAD hindered the conduction of electrical excitation and caused an unstable spiral wave, which could disrupt normal cardiac rhythm and cause atrial arrhythmia. This study showed the effects of excess reactive oxygen species on calcium cycling and action potential in atrial myocytes and provided insights regarding atrial arrhythmia induced by oxidative stress.

1. Introduction

Atrial fibrillation (AF) is the most common cardiac arrhythmia [1–3]. To design optimal treatment of AF, the mechanisms underlying AF need to be better understood. Both reactive oxygen species (ROS) and Ca²⁺/calmodulin- (CaM-) dependent protein kinase II (CaMKII) have been shown to be associated with the development of cardiac arrhythmias [4, 5]. The kinase CaMKII is ubiquitously expressed in the cardiomyocytes [4, 6]. It is involved in numerous cellular signaling cascades, such as phosphorylation of L-type Ca²⁺ channels [7–9], Na⁺ channels [10, 11], ryanodine receptors (RyRs) [12–18], and phospholamban (PLB) [12, 17, 19, 20]. Overexpression of CaMKII increases fractional sarcoplasmic reticulum (SR) Ca²⁺ release [13] and SR Ca²⁺ leakage [16], enhancing the activation of RyRs during both systole and diastole. As mentioned, it has been reported that CaMKII can also phosphorylate PLB [20]. In its unphosphorylated state, PLB acts as an endogenous inhibitor of sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA). Therefore, phosphorylation of PLB by CaMKII will enhance SERCA activity [21]. In addition to SR Ca²⁺ dynamics, many ion channels show CaMKII-dependent phosphorylation. Phosphorylation of L-type Ca²⁺ channels induced by CaMKII is associated with Ca²⁺-dependent Ca²⁺ current facilitation [22, 23]; CaMKII-dependent phosphorylation of Na⁺ channels may cause an increase in the late sodium current, which predisposes cardiomyocytes to arrhythmias [10, 11].

Conventionally, CaMKII is activated by Ca²⁺-bound calmodulin (Ca-CaM). Recently, a novel mechanism of CaMKII activation, which is ROS dependent, has been revealed [19, 24]. Previous studies have demonstrated that oxidative stress is closely associated with cardiac arrhythmias through alteration of the electrical activity and intracellular calcium dynamics of cardiac myocytes [25–27]. Oxidative stress is the main manifestation of cell metabolism disorders with excessive ROS. This oxidation-induced CaMKII activation is in relation to apoptosis [28], sinus node dysfunction [29], heart injury [30], and arrhythmias [31]. Although both ROS and CaMKII are associated with cardiac arrhythmias, the role of oxidative-dependent CaMKII activation in the development of atrial arrhythmias is not yet well understood.

In previous studies of cardiac arrhythmias, in vivo and in vitro experiments provide insights into mechanisms underlying arrhythmogenesis [5, 31]. However, these approaches have several limitations. For example, it is difficult to record the electrophysiological properties at different physical scales (from subcellular level to organ level) at the same time in one experiment. The electrophysiological properties recorded from different experiments, such as patch clamp and optical mapping, might be affected by gradual changes in cell or tissue properties. In addition to in vivo and in vitro experiments, a common method widely used in cardiac arrhythmia studies is in silico modeling. In this study, to simulate the effects of ROS-dependent CaMKII activation on regulations of intracellular Ca²⁺ and ionic currents, we first developed a Markov chain model of CaMKII, including both the autophosphorylation and oxidation pathways. Then, the effect of CaMKII on related proteins and ion channels was integrated to establish a computational model of the atrial cell to simulate changes in the electrophysiology under oxidative stress conditions, including ion currents, ion concentrations, calcium cycling, and transmembrane potential. Finally, the two-dimensional (2D) spiral wave was induced and the electrical excitation propagation was analyzed under normal and oxidative stress conditions.

2. Methods

2.1. CaMKII Model, including Both Autophosphorylation- and Oxidation-Dependent Activation

A six-state Markov chain model of CaMKII developed by our previous study of Zhang et al. [32] was incorporated into a human atrial model developed by Grandi et al. [33] to simulate autophosphorylation- and oxidation-dependent CaMKII activation.

The CaMKII monomer consists of three domains (an association domain, a regulatory domain, and a catalytic domain). When CaMCa4 (one calmodulin bound with four Ca²⁺) binds to the regulatory domain, the catalytic domain will be exposed, leading to CaMKII activation. In this state, CaMKII can be further autophosphorylated or oxidized, producing a long-lasting activation, even if CaMCa4 dissociates from CaMKII. Finally, CaMKII is fully deactivated by dephosphorylation with protein phosphatases, or reduction with methionine sulfoxide reductases (MsrA). Based on the four-state model of CaMKII developed by Chiba et al. [34], two oxidized states, with and without CaMCa4 binding with CaMKII, were added (see Figure S1 in the Supplementary Material for illustration), corresponding to the conformational change of oxidation-dependent CaMKII activation. Since there is no evidence to show that oxidation and autophosphorylation can occur at the same time, in this study, autophosphorylation and oxidation of CaMKII were treated as two independent processes [24]. The CaMKII autophosphorylation was based on the model developed by Chiba et al. [34]. The parameters of CaMKII oxidation were fitted to the experimental data recorded by Erickson et al. [24]. The detailed parameters of the CaMKII model are listed in Table S1 in the Supplementary Material.

2.2. Effects of CaMKII on Ion Channels

Previous studies showed that the ion channels influenced by CaMKII included the fast Na⁺ current () and the L-type Ca²⁺ current () [7–11]. According to the method developed by O’Hara et al. [35], the current of each ion channel () was divided into two parts:

where is the part of ion channel current not affected by CaMKII, is the part of ion channel current affected by CaMKII, and is the proportion affected by CaMKII using the following equations:

where is the fraction of CaMKII activation and is the Michaelis constant with the same value in O’Hara et al. [35]. For , the time constant of gate was slowed down by 1.46-fold [32]. And for , experimental data showed that CaMKII activation produced an increased amplitude and a slowed inactivation of [36]. Therefore, in this study, the time constant of gate was slowed down by 1.5-fold and the part of affected by CaMKII was increased by shown in Equation (3).

2.3. Effects of CaMKII on Ca²⁺ Cycling

Both RyR and PLB are primary regulatory proteins, controlling SR Ca²⁺ release and uptake; these are crucial processes, maintaining the balance of intracellular Ca²⁺. CaMKII activation-induced RyR phosphorylation can increase RyR opening probability and SR Ca²⁺ release. CaMKII activation-induced PLB phosphorylation can reduce SERCA inhibition, increasing SR Ca²⁺ uptake during diastole. The RyR and PLB models developed by Soltis and Saucerman [37] were used to simulate the effects of CaMKII activation on intracellular Ca²⁺ cycling. The rate constants controlling RyR opening () and SR leakage () were modified (Equations (4) and (5)) under CaMKII activation to increase RyR opening in both systolic and diastolic conditions. Meanwhile, the half maximal saturation constant of SERCA () was modified (Equation (6)) to mimic the phenomenon of increased SERCA pump calcium sensitivity, which is induced by Thr17 phosphorylation of PLB under CaMKII activation.

2.4. Single Atrial Cell Model

The cell membrane of an atrial myocyte was mimicked as an electrical circuit. The AP in a human atrial cell was calculated using the following ordinary differential equation [33]:

where is the cell capacitance, is the transmembrane voltage, is time, is the total transmembrane ionic current, and is the stimulus current. was calculated as

where is the fast Na⁺ current, is the background Na⁺ current, is the transient outward K⁺ current, is the ultrarapid delayed rectifier K⁺ current, is the rapid activating K⁺ current, is the slowly activating K⁺ current, is the inward rectifier K⁺ current, is the Na⁺/K⁺ pump current, is the L-type Ca²⁺ current, is the background Ca²⁺ current, is the sarcolemmal Ca²⁺ pump current, is the Na⁺/Ca²⁺ exchange current, is the Ca²⁺-activated Cl⁻ current, and is the background Cl⁻ current.

In this study, 0.5, 1, and 2 Hz stimulation frequencies were used to investigate the frequency dependency of CaMKII activation. Action potentials of atrial cells were produced by applying a stimulus current with an amplitude of −12.5 pA/pF and a duration of 5 ms. The time step used in the simulation was 0.1 ms using the CVODE solver of the SUite of Nonlinear and DIfferential/ALgebraic equation Solvers (SUNDIALS) to solve initial value problems for ordinary differential equation systems. These simulation results were consistent with the ones using forward Euler with a time step of 0.005 ms. Under the control condition, the value of H₂O₂ in the CaMKII model was 0 μM, and the one was 200 μM under the oxidative stress condition. To ensure that the model reached a quasistable steady state, the simulations were carried out for more than 50 s under the control condition and 150 s under the oxidative stress condition.

2.5. Excitation Wave Conduction in 2D Tissue

In this study, a monodomain model was used for 2D simulations in an ideal square tissue. The excitation wave propagation was simulated using the following equation:

where is the diffusion tensor describing the conductivity of the tissue and is the spatial gradient operator.

In 2D simulations, the 2D ideal isotropic tissue was constructed as a sheet of nodes to investigate the stability of the spiral waves, as used in [38]. The time step was set at 0.005 ms, and the space step was 0.1 mm along both direction axes, as used in [39]. The S1–S2 stimulation protocol was used to produce spiral waves. An S1 stimulus was applied to the five column nodes on the left side to induce a plane wave (Figure 1(a)). After the nodes in the middle line (the dashed line in Figure 1(b)) reached the end point of the refractory period, an S2 stimulus was applied in the left lower quadrant to induce a spiral wave (Figure 1(b)).

(a)

(b)

Figure 1: The S1–S2 protocol used in 2D simulations. The red bar in (a) and the red square in (b) represent the regions where S1 and S2 stimuli were applied.

The conduction velocity (CV) was measured using the following equation:

where and are the positions and and are the times when the excitation reaches the two recording points. In 2D simulation, the tissue was isotropic and the diffusion tensor () along both direction axes was set to 0.029 mm²/ms to produce a CV of 69 cm/s in control conditions, which was consistent with the CV (70 cm/s) previously used for human atrial tissue simulation [40–42].

3. Results

3.1. Electrophysiological Properties of a Single Atrial Cell under Control Conditions

The action potential (AP) in a single-cell model was generated by applying a series of 1 Hz stimuli. The time trace of the simulated AP is presented in Figure 2, along with other AP traces from previous mathematical models [33, 43, 44]. The resting membrane potential (RMP), maximum upstroke velocity (), action potential amplitude (APA), and action potential duration at 90% repolarization (APD₉₀) were measured as −75.2 mV, 160 mV/ms, 110 mV, and 312 ms, respectively. These parameters were consistent with previous models developed by Grandi et al. [33], Courtemanche et al. [43], and Nygren et al. [44], and APD₉₀ was within the scope of previously reported experimental data [45–48], as listed in Table 1.

Figure 2: Action potential traces produced by our simulation and previous mathematical models and experiments at 1 Hz. Traces labeled “Courtemanche et al.,” “Nygren et al.,” and “Grandi et al.” were reproduced from mathematical models [33, 43, 44].

Table 1: Comparison of action potential characteristics.

Model and experimental data
RMP
(mV)

(mV/ms)
APA
(mV)
APD₉₀
(ms)

Courtemanche et al. [43]
−79.5
216
106
297

Nygren et al. [44]
−76.4
116
103.8
245

Grandi et al. [33]
−74.5
135
105
294

Our model
−75.2
160
110
312

Dawodu et al. [45]
—
—
—

Katoh et al. [46]
—
—
—

Bosch et al. [47]
—
—
—

Kim et al. [48]
—
—
—

To further validate the model, the ionic currents and intracellular calcium cycling process were investigated and compared with those from the model of Grandi et al. [33] at 1 Hz, as shown in Figure 3. These results demonstrated an elevation in the repolarization phase of the AP and, therefore, a longer action potential duration (APD) in our model (Figure 3(a)). This was attributed to the change in . In our model, showed a larger amplitude and a slowed inactivation (Figure 3(b)), which was consistent with the effect of CaMKII activation on [36]. In addition to the AP, the variation in altered intracellular calcium regulation, as was the main influx of intracellular Ca²⁺. The increase in ICaL caused a larger Ca²⁺ influx, leading to a larger SR Ca²⁺ release (Figure 3(c)) and thus a larger intracellular Ca²⁺ concentration (Figure 3(d)). The accumulation of intracellular Ca²⁺ finally caused an elevation in the concentration of SR Ca²⁺ (Figure 3(e)). Interestingly, the peak inward in our model declined even when the intracellular Ca²⁺ concentration increased (Figure 3(f)). This may be explained by the fact that the elevation of the membrane potential during the repolarization phase suppressed the inward Na⁺ flux and decreased .

Figure 3: Traces of AP, ionic currents, and calcium cycling produced by our model (red) and the model of Grandi et al. [33] (black). (a) AP. (b) . (c) SR Ca²⁺ release (). (d) intracellular Ca²⁺ concentration (). (e) SR Ca²⁺ concentration ([Ca²⁺]SR). (f) . (g) SR Ca²⁺ reuptake (). (h) SR Ca²⁺ leakage (). (i) . (j) .

To justify the increase of intracellular Ca²⁺ concentration in our model, the Ca²⁺ concentration was compared with previous studies. The calcium transient amplitude in our model was ~0.53 μM, which was consistent with the studies of Courtemanche et al. and Colman et al. (0.45–0.55 μM) [43, 49], although the Ca²⁺ transient amplitude in the study of Nygren et al. was reported to be even larger (~1.2 μM) [44].

The variations in and intracellular calcium cycling were caused by adding the effect of CaMKII activation in our model. Previous studies demonstrated that CaMKII activation was significantly influenced by pacing frequencies. In our simulation, CaMKII activation and intracellular Ca²⁺ dynamics were investigated at 0.5, 1, and 2 Hz. Figure 4(a)shows that the fraction of CaMKII activation increased with increasing pacing rate. Meanwhile, the time required for CaMKII activation was also decreased, indicating a faster activation at a higher pacing rate (Figure 4(a)). The frequency-dependent activation of CaMKII was mainly due to the intracellular Ca²⁺ concentration. With increasing pacing rate, both RyR Ca²⁺ release and intracellular Ca²⁺ concentration increased (Figure 4(b)), causing an increase in CaMCa4 and, therefore, augmenting CaMKII activation.

(a)

(b)

Figure 4: Frequency-dependent activation of CaMKII. (a) CaMKII activation curves at 0.5, 1, and 2 Hz in 80 s. (b) Normalized CaMKII activation, RyR calcium release (RyR Activity), and intracellular Ca²⁺ transient (CaT) to 2 Hz at 0.5, 1, and 2 Hz.

3.2. Electrophysiological Properties of a Single Atrial Cell under Oxidative Stress Conditions

The simulation results showed that an EAD was induced in the AP under conditions of oxidative stress at a pacing rate of 1 Hz (Figure 5(a)). In this case, the elevated ROS concentration (0.2 mM) induced a significant increase in CaMKII activation (Figure 5(b)). ROS-enhanced CaMKII activation further enlarged , , fraction of RyR phosphorylation, and fraction of PLB phosphorylation (Figures 5(c)–5(f)). Under oxidative stress conditions, the intracellular Ca²⁺ dynamics were remarkably changed. Enhanced RyR and PLB phosphorylation caused significant increases in , , and (Figures 5(g)–5(i)), which further induced a dramatic increase in intracellular Ca²⁺ concentration (Figure 5(j)) and a decrease in minimum SR Ca²⁺ (Figure 5(k)). However, there was no obvious change in the intracellular Na⁺ concentration (Figure 5(l)).

Figure 5: Traces of AP, ionic currents, and calcium cycling under normal (black) and oxidative stress (red) conditions. (a) AP. (b) CaMKII activity. (c) . (d) . (e) Fraction of RyR phosphorylation. (f) Fraction of PLB phosphorylation. (g) . (h) . (i) . (j) . (k) [Ca²⁺]SR. (l) Intracellular Na⁺ concentration ().

As CaMKII activation showed frequency-dependent behavior, ROS-induced EADs are rate dependent as well. Under oxidative stress conditions, the occurrence of EADs in 50 s became more prominent with decreasing pacing rate, as shown in Figure 6. Under control conditions, the CaMKII activation (without ROS-induced activation) gradually decreased with decreasing pacing rate (Figure 4(b)). This meant that more CaMKII could be activated by ROS under oxidative conditions at low pacing rates. Consequently, the CaMKII activation induced by ROS gradually increased with decreasing pacing rate (green bars in Figure 6). Therefore, under oxidative stress conditions, more EADs were induced at low pacing rates.

Figure 6: Number of EAD occurring in the same duration (orange) and the fraction of CaMKII activation induced by ROS (green) at different pacing rates under oxidative stress conditions.

3.3. Mechanisms Underlying the Genesis of EAD under Oxidative Stress Conditions

In the case of EAD, and were the main currents significantly influenced by oxidative stress. Blocking at the point of reactivation abolished EAD (data not shown), implying that was a main factor contributing to the genesis of EAD. However, blocking did not guarantee elimination of EAD. Figure 7 indicated the effects of blocking at different times from the point of reactivation (70 ms) to the time of the peak of EAD (350 ms). With the delay in blocking (Figure 7(a)), the reactivation of became more prominent (Figure 7(b)) and the APD gradually increased (Figure 7(c)). When blocking in the interval 70–210 ms, the reactivation of did not induce AP depolarization and therefore abolished EAD. On blocking during the interval 210–280 ms, the enlarged reactivation of started to transfer AP repolarization to depolarization. However, in this time interval, AP depolarization was not obvious and no EAD was induced. When blocking during the interval 280–350 ms, the reactivation of was able to induce AP depolarization. Consequently, in this time interval, blocking could not eliminate EAD (Figure 7(c)). Therefore, played an important role in reactivation and thus in triggering EAD.

Figure 7: Effects of blocking at different times on the genesis of EAD under oxidative stress condition: (a) . (b) . (c) AP. The gray box represents the time interval during which EAD can occur while blocking .

In our model, was not directly regulated by ROS. As the main efflux of intracellular Ca²⁺, was regulated by intracellular Ca²⁺ concentration; Ca²⁺ cycling was another important process affected by ROS-induced CaMKII activation. Three factors associated with intracellular cycling in our model were directly regulated by ROS-induced CaMKII activation, namely, , , and . Under oxidative stress conditions, intracellular Ca²⁺ concentration significantly increased (Figure 5(j)). To investigate the role of regulation in the genesis of EAD, decreased , reduced , and enhanced were independently applied to reduce the under oxidative stress conditions, as shown in Figure 8. Figure 8(a) shows that decreasing by 10% abolished EAD. This was due not only to decline but also to a decrease in reactivation. Reducing SR Ca²⁺ release postponed the occurrence of EAD but did not successfully eliminate EAD (Figure 8(b)). Although a reduction in SR Ca²⁺ release temporarily reduced , it caused Ca²⁺ accumulation in SR, which finally induced elevation (data not shown) and EAD. Increasing SR Ca²⁺ uptake by 10% had a similar effect to that of reducing SR calcium release (Figure 8(c)). Nonetheless, further increasing by 20%–30% abolished EAD. In these two cases, the SR Ca²⁺ uptake and release reached a new balance and maintained and [Ca²⁺]SR at a steady state (data not shown).

Figure 8: Effect of (a) reducing by 10%, (b) reducing by 10%, 20%, and 30%, and (c) increasing by 10%, 20%, and 30% on the genesis of EAD under oxidative stress conditions.

3.4. Effect of Oxidative Stress on Excitation Wave Propagation in 2D Tissue

The excitation wave propagation was investigated in 2D ideal tissue under control and oxidative stress conditions. Spiral waves were induced using the S1–S2 protocol. Under control conditions with the S1–S2 interval of 550 ms, a stationary spiral wave was produced, with its tip anchored in the central region of the square tissue, as shown in Figure 9. Under oxidative stress conditions, the ROS concentration was set to 0.2 mM for all nodes in the tissue. In this case with the S1–S2 interval of 730 ms, a nonstationary spiral wave was induced. As the EAD was induced under this condition, the region to which the S2 stimulus was applied showed a prolonged depolarization (Figure 10, 1125 ms) and a second wavefront was induced at 1450 ms, as shown in Figure 10. When the wavefronts reached the region to which the S2 stimulus was applied (Figure 10, 1650 ms), the nodes were not fully repolarized, owing to the prolonged APD induced by EAD. Therefore, the excitation wave propagation was suppressed (Figure 10, 1650 ms). Finally, the tip of the spiral wave wandered across the tissue and produced a nonstationary spiral wave. When oxidative stress only occurred in the right third of the tissue, the spatial heterogeneity further aggravated the instability of the nonstationary spiral wave and gave rise to a breakup of the spiral wave (data not shown).

Figure 9: Spiral wave induced under control conditions using the S1–S2 protocol. The number in each panel represents the recording time.

Figure 10: Spiral wave induced under oxidative stress conditions using the S1–S2 protocol. The number in each panel represents the recording time.

4. Discussion

In this study, the effects of ROS-dependent CaMKII activation on regulations of intracellular Ca²⁺ and ionic currents were investigated in silico using the updated human atrial cell model including the CaMKII model with autophosphorylation- and oxidation-dependent activation at the single-cell level and 2D tissue model. Our major findings follow: (i) A new Markov chain model of CaMKII, including both the autophosphorylation and oxidation pathways, was developed, and the effects of CaMKII on ion channels and Ca²⁺ cycling were incorporated into the computational model of the atrial myocyte. (ii) The mechanisms of oxidative stress-induced EADs in atrial cells were thoroughly examined which are helpful to further understand or investigate the mechanisms underlying oxidative stress-induced AF. The simulation results at the single-cell level illustrated that oxidative stress resulted in EADs of the AP at the normal pacing rate of 1 Hz. It was contributed by reactivation of and intracellular Ca²⁺ elevation induced by CaMKII activation under oxidative stress conditions. (iii) The 2D simulations provide insights into reentry in atrial tissue under the oxidative stress condition, which plays significant roles in AF mechanisms. In 2D simulation, oxidative stress aggravated the instability of the excitation wave and gave rise to a nonstationary spiral wave, owing to the EAD induced by oxidative stress hindering the electrical conduction.

4.1. Model Development of ROS-Dependent CaMKII Activation

Previous experimental studies have revealed that both ROS and CaMKII are associated with the development of atrial arrhythmia [4, 5]. In addition, it has recently been reported that CaMKII can keep persistent activity by oxidation [24]. However, the role of ROS-dependent CaMKII activation in the genesis of atrial fibrillation is not yet well understood. In this study, we developed a computational model of human atrial cell including both ROS-dependent and autophosphorylation-dependent CaMKII activation to investigate effects of CaMKII activation on atrial electrophysiology under oxidative stress condition. In this CaMKII model, a pathway of ROS-dependent CaMKII activation was introduced which was different from previous CaMKII models only including the autophosphorylation-dependent CaMKII activation [34, 35, 50]. Christensen et al. [39] also developed a model of CaMKII activity including oxidation and autophosphorylation activation pathways. However, different from their model, the oxidation and autophosphorylation activation pathways in our six-state CaMKII are mutual independence, implying that there is no such a state with both oxidation and autophosphorylation activation simultaneously. This fact is consistent with the experimental observations reported in [24].

The effect of CaMKII on , , RyR, and PLB was integrated into the human atrial myocyte model, producing action potential characteristics (such as RMP, , APA, APD₉₀, and the calcium transient amplitude) consistent with previous models and experimental data of human atrial myocytes [33, 43–49] (Figure 2 and Table 1).

4.2. Mechanisms Underlying the Genesis of EAD under Oxidative Stress Conditions

The generation of EADs associated with ectopic (triggered) activity can be contributed by CaMKII activation [51–55]. Animal experiments have shown that oxidative stress can promote EADs in the ventricle of guinea pigs and rabbits [56]. Therefore, this study explored the intrinsic mechanisms of EADs associated with atrial arrhythmia induced by oxidative stress, using the developed CaMKII model.

First, simulation results demonstrated that the characteristics of AP and calcium cycling generated by our model were consistent with previous computational and experimental studies [43–49]. After model validation, the role of oxidation-dependent CaMKII activation on the AP was investigated in a single atrial cell model. The simulation results demonstrated that, under oxidative stress conditions, increasing ROS concentrations in the cytoplasm enhanced CaMKII activation and consequently augmented , RyR phosphorylation, and PLB phosphorylation (Figures 5(c), 5(e), and (f)), inducing an enhanced Ca²⁺ influx, a larger SR Ca²⁺ release, and a promoted SR Ca²⁺ uptake (Figures 5(g) and 5(h)). These effects generated a remarkable elevation in intracellular Ca²⁺ concentration and therefore enlarged via promoted Ca²⁺ extrusion (Figure 5(d)). augmentation provided a depolarization component to counterbalance the repolarization reserve and prolong the APD. The elevated and prolonged AP in phase 3 gradually caused reactivation of (Figure 5(c)), which ultimately induced action potential depolarization and produced EAD (Figure 5(a)), consistent with previous studies [57–59].

In this process, as was the direct depolarization current inducing EAD, blocking completely abolished EAD. However, blocking at different time points had different effects on EAD. If was blocked before reactivation could induce AP depolarization, EAD was eliminated (Figure 7). This phenomenon suggested that augmentation acted as a trigger of reactivation. EAD can be induced only when the trigger is large enough. Under oxidative stress condition, ROS elevation did not affect directly and the increase of primarily accounted for augmentation. Therefore, the effects of ROS-induced CaMKII activation on intracellular calcium cycling also played a crucial role in the genesis of EAD.

To examine the role of ROS-induced elevation in the genesis of EAD, Ca²⁺ influx and SR Ca²⁺ release and uptake were modified to reduce ROS-induced elevation. First, was reduced by 10% to inhibit Ca²⁺ influx, causing a decline of . Together with the decreased reactivation of , the decline caused by reducing abolished EAD (Figure 8(a)). Second, partial inhibition of SR Ca²⁺ release (by 10%–30%) temporarily reduced . However, inhibition of SR Ca²⁺ release gradually caused SR Ca²⁺ accumulation and ultimately gave rise to further elevation, which accounted for the fact that partial inhibition of SR Ca²⁺ release could only postpone the occurrence of EAD but could not abolish EAD (Figure 8(b)). Increasing SR Ca²⁺ uptake by 10% had a similar effect to partial inhibition of SR Ca²⁺ release. Interestingly, on further increasing SR Ca²⁺ uptake by 20% and 30%, declined and EAD was abolished. This phenomenon might be attributed to a new balance of increased SR Ca²⁺ uptake and ROS-induced release (Figure 8(c)). Therefore, the dynamic balance between SR Ca²⁺ release and uptake is the key factor in maintaining at a steady state.

4.3. Effect of Rate-Dependent CaMKII Activation on EAD under Oxidative Stress Conditions

Under control conditions, CaMKII activation gradually decreased with a decrease in the pacing rate, as Ca²⁺ concentration gradually decreased at low pacing rates. However, under oxidative stress conditions, once CaMKII was activated, the high ROS concentration tended to maintain CaMKII activation. Therefore, under oxidative stress conditions, fractions of CaMKII activation at different pacing rates were similar and maintained at a high level (about 75%). This indicated that, at low pacing rates, ROS-induced CaMKII activation played a more important role. At low pacing rates, intracellular Ca²⁺ concentration could be dramatically increased by ROS-induced CaMKII activation and therefore could induce more EADs (Figure 6). Although the heart rate is not so slow such as 0.5 Hz and 0.33 Hz under normal conditions, it may be such slow under some pathological conditions, such as bradycardia [60]. The purpose of exploring up to such slow pacing frequencies is to investigate whether the occurrence of EAD is linear with the decrease of pacing frequency or not. The results showed that the occurrence of EAD was more dependent on ROS-induced CaMKII activation rather than the pacing frequencies.

4.4. Effect of Oxidative Stress on Excitation Wave Propagation in 2D Tissue

The spiral wave is a well-known approach to understanding reentry in cardiac tissue, which plays significant roles in AF mechanisms [61]. This study investigated the stability of the spiral wave in 2D ideal tissue under control and oxidative stress conditions. In tissue with a uniform distribution of oxidative stress (Figure 10), the region to which the S2 stimulus was applied had a prolonged depolarization phase, owing to the ROS-induced EAD, forming a spiral wave barrier. In addition, the ROS-induced EAD produced a second excitation wavefront during wave propagation. These effects finally produced a nonstationary spiral wave. In the case of nonuniform distribution of oxidative stress, a breakup of the spiral wave was induced, which might disrupt normal cardiac rhythm and cause atrial arrhythmia [62–64].

4.5. Model Limitations and Future Directions

The present model was based on models of Grandi et al. [33] and O’Hara et al. [35] and thus inherited the same limitations of both models. The main limitation of the current model, which will be addressed in the future, is that the ROS concentration is constant. In future versions, pathways of ROS production and scavenging should be added. Another issue for model development is to expand the model into real 2D heart tissue and 3D organ models, to explore electrical wave propagation induced by abnormal action potentials, and electrocardiography in real heart geometry.

Although the S1–S2 interval in 2D simulation to produce the spiral wave was dependent on the tissue size, it is reasonable to investigate the difference in the stability of the spiral wave using the same tissue size under control and oxidative stress conditions. Further work is required to measure the range of S1–S2 intervals for producing the spiral wave (vulnerable window) varying from control to oxidative stress conditions to further determine the proarrhythmic effect of oxidative stress.

5. Conclusions

In this study, we have investigated the role of oxidation-dependent CaMKII activation in the genesis of abnormal action potentials in atria. It was shown that, at the atrial cell level, oxidation-dependent CaMKII activation-induced reactivation and elevation contributed to EAD generation. Moreover, at the 2D tissue level, oxidative stress-induced EAD contributed to the instability of excitation waves, facilitating atrial arrhythmia. This study investigated the role of ROS-dependent CaMKII activation in the development of atrial arrhythmias, shedding light on the genesis of atrial arrhythmias under oxidative stress conditions.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Qince Li, Henggui Zhang, and Na Zhao conceived the study and designed the simulations. Na Zhao and Haibo Sui developed the model and conducted the simulations. The paper was written by Qince Li. Qince Li and Henggui Zhang analyzed the simulation results. Qince Li and Na Zhao contributed equally to this work. All the authors read and approved the final manuscript.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (grant numbers 61601143, 81770328, and 61572152) and the Heilongjiang and China Postdoctoral Science Foundation (grant number 2015M581448).

Supplementary Materials

Figure S1: The CaMKII model. B, CaMKII-CaMCa4 (bound state); BO, CaMKIIOx-CaMCa4 (oxidized and bound state); BP, CaMKIIP-CaMCa4 (phosphorylated and bound state); I, CaMKII (inactive state); O, CaMKIIOx (oxidized state); P, CaMKIIP (phosphorylated state). Table S1: Parameters of CaMKII model (Supplementary Materials)

References

C. R. Wyndham, “Atrial fibrillation: the most common arrhythmia,” Texas Heart Institute Journal, vol. 27, no. 3, pp. 257–267, 2000.

M. C. E. F. Wijffels, C. J. H. J. Kirchhof, R. Dorland, J. Power, and M. A. Allessie, “Electrical remodeling due to atrial fibrillation in chronically instrumented conscious goats,” Circulation, vol. 96, no. 10, pp. 3710–3720, 1997.

M. A. Allessie, “Atrial electrophysiologic remodeling: another vicious circle?” Journal of Cardiovascular Electrophysiology, vol. 9, no. 12, pp. 1378–1393, 1998.

S. Wagner, A. G. Rokita, M. E. Anderson, and L. S. Maier, “Redox regulation of sodium and calcium handling,” Antioxidants & Redox Signaling, vol. 18, no. 9, pp. 1063–1077, 2013.

P. D. Swaminathan, A. Purohit, T. J. Hund, and M. E. Anderson, “Calmodulin-dependent protein kinase II: linking heart failure and arrhythmias,” Circulation Research, vol. 110, no. 12, pp. 1661–1677, 2012.

T. Zhang and J. H. Brown, “Role of Ca²⁺/calmodulin-dependent protein kinase II in cardiac hypertrophy and heart failure,” Cardiovascular Research, vol. 63, no. 3, pp. 476–486, 2004.

A. Hudmon, H. Schulman, J. Kim, J. M. Maltez, R. W. Tsien, and G. S. Pitt, “CaMKII tethers to L-type Ca²⁺ channels, establishing a local and dedicated integrator of Ca²⁺ signals for facilitation,” The Journal of Cell Biology, vol. 171, no. 3, pp. 537–547, 2005.

I. Dzhura, Y. Wu, R. J. Colbran, J. R. Balser, and M. E. Anderson, “Calmodulin kinase determines calcium-dependent facilitation of L-type calcium channels,” Nature Cell Biology, vol. 2, no. 3, pp. 173–177, 2000.

A. Blaich, A. Welling, S. Fischer et al., “Facilitation of murine cardiac L-type Cav1.2 channel is modulated by calmodulin kinase II-dependent phosphorylation of S1512 and S1570,” Proceedings of the National Academy of Sciences, vol. 107, no. 22, pp. 10285–10289, 2010.

S. Wagner, N. Dybkova, E. C. L. Rasenack et al., “Ca²⁺/calmodulin-dependent protein kinase II regulates cardiac Na⁺ channels,” The Journal of Clinical Investigation, vol. 116, no. 12, pp. 3127–3138, 2006.

S. Wagner, H. M. Ruff, S. L. Weber et al., “Reactive oxygen species–activated Ca/calmodulin kinase IIδ is required for LateINaAugmentation leading to cellular Na and Ca overload,” Circulation Research, vol. 108, no. 5, pp. 555–565, 2011.

L. S. Maier and D. M. Bers, “Role of Ca²⁺/calmodulin-dependent protein kinase (CaMK) in excitation–contraction coupling in the heart,” Cardiovascular Research, vol. 73, no. 4, pp. 631–640, 2007.

L. S. Maier, T. Zhang, L. Chen, J. DeSantiago, J. H. Brown, and D. M. Bers, “Transgenic CaMKIIδCOverexpression uniquely alters cardiac myocyte Ca²⁺Handling,” Circulation Research, vol. 92, no. 8, pp. 904–911, 2003.

X. H. T. Wehrens, S. E. Lehnart, S. R. Reiken, and A. R. Marks, “Ca²⁺/calmodulin-dependent protein kinase II phosphorylation regulates the cardiac ryanodine receptor,” Circulation Research, vol. 94, no. 6, pp. e61–e70, 2004.

T. Guo, T. Zhang, R. Mestril, and D. M. Bers, “Ca²⁺/calmodulin-dependent protein kinase II phosphorylation of ryanodine receptor does affect calcium sparks in mouse ventricular myocytes,” Circulation Research, vol. 99, no. 4, pp. 398–406, 2006.

M. Kohlhaas, T. Zhang, T. Seidler et al., “Increased sarcoplasmic reticulum calcium leak but unaltered contractility by acute CaMKII overexpression in isolated rabbit cardiac myocytes,” Circulation Research, vol. 98, no. 2, pp. 235–244, 2006.

C. M. Sag, H. A. Wolff, K. Neumann et al., “Ionizing radiation regulates cardiac Ca handling via increased ROS and activated CaMKII,” Basic research in cardiology, vol. 108, no. 6, p. 385, 2013.

H.-T. Ho, B. Liu, J. S. Snyder et al., “Ryanodine receptor phosphorylation by oxidized CaMKII contributes to the cardiotoxic effects of cardiac glycosides,” Cardiovascular Research, vol. 101, no. 1, pp. 165–174, 2014.

J. R. Erickson, B. J. He, I. M. Grumbach, and M. E. Anderson, “CaMKII in the cardiovascular system: sensing redox states,” Physiological Reviews, vol. 91, no. 3, pp. 889–915, 2011.

H. K. Simmerman, J. H. Collins, J. L. Theibert, A. D. Wegener, and L. R. Jones, “Sequence analysis of phospholamban. Identification of phosphorylation sites and two major structural domains,” Journal of Biological Chemistry, vol. 261, no. 28, pp. 13333–13341, 1986.

A. G. Brittsan and E. G. Kranias, “Phospholamban and cardiac contractile function,” Journal of Molecular and Cellular Cardiology, vol. 32, no. 12, pp. 2131–2139, 2000.

M. E. Anderson, A. P. Braun, H. Schulman, and B. A. Premack, “Multifunctional Ca²⁺/calmodulin-dependent protein kinase mediates Ca(2+)-induced enhancement of the L-type Ca²⁺ current in rabbit ventricular myocytes,” Circulation Research, vol. 75, no. 5, pp. 854–861, 1994.

W. Yuan and D. M. Bers, “Ca-dependent facilitation of cardiac Ca current is due to Ca-calmodulin-dependent protein kinase,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 267, no. 3, pp. H982–H993, 1994.

J. R. Erickson, M.-l. A. Joiner, X. Guan et al., “A dynamic pathway for calcium-independent activation of CaMKII by methionine oxidation,” Cell, vol. 133, no. 3, pp. 462–474, 2008.

Q. Li, D. Su, B. O'Rourke, S. M. Pogwizd, and L. Zhou, “Mitochondria-derived ROS bursts disturb Ca²⁺ cycling and induce abnormal automaticity in guinea pig cardiomyocytes: a theoretical study,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 308, no. 6, pp. H623–H636, 2015.

Q. Li, S. M. Pogwizd, S. D. Prabhu, and L. Zhou, “Inhibiting Na⁺/K⁺ ATPase can impair mitochondrial energetics and induce abnormal Ca²⁺ cycling and automaticity in guinea pig cardiomyocytes,” PLoS One, vol. 9, no. 4, p. e93928, 2014.

D. B. Zorov, M. Juhaszova, and S. J. Sollott, “Mitochondrial ROS-induced ROS release: an update and review,” Biochimica et Biophysica Acta (BBA)-Bioenergetics, vol. 1757, no. 5–6, pp. 509–517, 2006.

J. Palomeque, O. V. Rueda, L. Sapia et al., “Angiotensin II–induced oxidative stress resets the Ca²⁺ dependence of Ca²⁺–calmodulin protein kinase II and promotes a death pathway conserved across different species,” Circulation Research, vol. 105, no. 12, pp. 1204–1212, 2009.

P. D. Swaminathan, A. Purohit, S. Soni et al., “Oxidized CaMKII causes cardiac sinus node dysfunction in mice,” The Journal of Clinical Investigation, vol. 121, no. 8, pp. 3277–3288, 2011.

B. J. He, M.-l. A. Joiner, M. V. Singh et al., “Oxidation of CaMKII determines the cardiotoxic effects of aldosterone,” Nature Medicine, vol. 17, no. 12, pp. 1610–1618, 2011.

L.-H. Xie, F. Chen, H. S. Karagueuzian, and J. N. Weiss, “Oxidative stress–induced afterdepolarizations and calmodulin kinase II signaling,” Circulation Research, vol. 104, no. 1, pp. 79–86, 2009.

S. Zhang, Q. Li, L. Zhou, K. Wang, and H. Zhang, “Development of a novel Markov chain model for oxidative-dependent CaMKIIδ activation,” in 2015 Computing in Cardiology Conference (CinC), pp. 881–884, Nice, France, September, 2015, IEEE.

E. Grandi, S. V. Pandit, N. Voigt et al., “Human atrial action potential and Ca²⁺ model,” Circulation Research, vol. 109, no. 9, pp. 1055–1066, 2011.

H. Chiba, N. S. Schneider, S. Matsuoka, and A. Noma, “A Simulation Study on the Activation of Cardiac CaMKII _δ_ -Isoform and Its Regulation by Phosphatases,” Biophysical Journal, vol. 95, no. 5, pp. 2139–2149, 2008.

T. O'Hara, L. Virág, A. Varró, and Y. Rudy, “Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation,” PLoS computational biology, vol. 7, no. 5, article e1002061, 2011.

T. Christ, P. Boknik, S. Wöhrl et al., “L-type Ca²⁺ current downregulation in chronic human atrial fibrillation is associated with increased activity of protein phosphatases,” Circulation, vol. 110, no. 17, pp. 2651–2657, 2004.

A. R. Soltis and J. J. Saucerman, “Synergy between CaMKII substrates and β-adrenergic signaling in regulation of cardiac myocyte Ca²⁺ handling,” Biophysical Journal, vol. 99, no. 7, pp. 2038–2047, 2010.

K. H. W. J. ten Tusscher and A. V. Panfilov, “Alternans and spiral breakup in a human ventricular tissue model,” American Journal of Physiology-heart and Circulatory Physiology, vol. 291, no. 3, pp. H1088–H1100, 2006.

M. D. Christensen, W. Dun, P. A. Boyden, M. E. Anderson, P. J. Mohler, and T. J. Hund, “Oxidized calmodulin kinase II regulates conduction following myocardial infarction: a computational analysis,” PLOS Computational Biology, vol. 5, no. 12, article e1000583, 2009.

H. Ni, D. G. Whittaker, W. Wang, W. R. Giles, S. M. Narayan, and H. Zhang, “Synergistic anti-arrhythmic effects in human atria with combined use of sodium blockers and acacetin,” Frontiers in Physiology, vol. 8, 2017.

H. Ni, I. Adeniran, and H. Zhang, “_In-silico_ investigations of the functional impact of _KCNA5_ mutations on atrial mechanical dynamics,” Journal of Molecular and Cellular Cardiology, vol. 111, pp. 86–95, 2017.

D. G. Whittaker, M. A. Colman, H. Ni, J. C. Hancox, and H. Zhang, “Human atrial arrhythmogenesis and sinus bradycardia in KCNQ1-linked short QT syndrome: insights from computational modelling,” Frontiers in Physiology, vol. 9, article 1402, 2018.

M. Courtemanche, R. J. Ramirez, and S. Nattel, “Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 275, no. 1, pp. H301–H321, 1998.

A. Nygren, C. Fiset, L. Firek et al., “Mathematical model of an adult human atrial cell,” Circulation research., vol. 82, no. 1, pp. 63–81, 1998.

A. A. Dawodu, F. Monti, K. Iwashiro, M. Schiariti, R. Chiavarelli, and P. E. Puddu, “The shape of human atrial action potential accounts for different frequency- related changes in vitro,” International Journal of Cardiology, vol. 54, no. 3, pp. 237–249, 1996.

H. Katoh, T. Shinozaki, S. Baba et al., “Monophasic action potential duration at the crista terminalis in patients with sinus node disease,” Circulation, vol. 69, no. 11, pp. 1361–1367, 2005.

R. F. Bosch, X. Zeng, J. B. Grammer, K. Popovic, C. Mewis, and V. Kühlkamp, “Ionic mechanisms of electrical remodeling in human atrial fibrillation,” Cardiovascular Research, vol. 44, no. 1, pp. 121–131, 1999.

B.-S. Kim, Y.-H. Kim, G.-S. Hwang et al., “Action potential duration restitution kinetics in human atrial fibrillation,” Journal of the American College of Cardiology, vol. 39, no. 8, pp. 1329–1336, 2002.

M. A. Colman, O. V. Aslanidi, S. Kharche et al., “Pro-arrhythmogenic effects of atrial fibrillation-induced electrical remodelling: insights from the three-dimensional virtual human atria,” The Journal of Physiology, vol. 591, no. 17, pp. 4249–4272, 2013.

T. J. Hund and Y. Rudy, “Rate dependence and regulation of action potential and calcium transient in a canine cardiac ventricular cell model,” Circulation, vol. 110, no. 20, pp. 3168–3174, 2004.

M. E. Anderson, A. P. Braun, Y. Wu et al., “KN-93, an inhibitor of multifunctional Ca⁺⁺/calmodulin-dependent protein kinase, decreases early afterdepolarizations in rabbit heart,” Journal of Pharmacology and Experimental Therapeutics, vol. 287, no. 3, pp. 996–1006, 1999.

Y. Wu, L. B. MacMillan, R. B. McNeill, R. J. Colbran, and M. E. Anderson, “CaM kinase augments cardiac L-type Ca²⁺ current: a cellular mechanism for long Q-T arrhythmias,” American Journal of Physiology-heart and Circulatory Physiology, vol. 276, no. 6, pp. H2168–H2178, 1999.

Y. Wu, J. Temple, R. Zhang et al., “Calmodulin kinase II and arrhythmias in a mouse model of cardiac hypertrophy,” Circulation, vol. 106, no. 10, pp. 1288–1293, 2002.

D. M. Bers and S. Morotti, “Ca²⁺ current facilitation is CaMKII-dependent and has arrhythmogenic consequences,” Frontiers in Pharmacology, vol. 5, article 144, 2014.

J. Heijman, N. Voigt, X. H. T. Wehrens, and D. Dobrev, “Calcium dysregulation in atrial fibrillation: the role of CaMKII,” Frontiers in Pharmacology, vol. 5, article 30, 2014.

Y. Song, J. C. Shryock, S. Wagner, L. S. Maier, and L. Belardinelli, “Blocking late sodium current reduces hydrogen peroxide-induced arrhythmogenic activity and contractile dysfunction,” Journal of Pharmacology and Experimental Therapeutics, vol. 318, no. 1, pp. 214–222, 2006.

X. Y. Qi, Y.-H. Yeh, D. Chartier et al., “The calcium/calmodulin/kinase system and arrhythmogenic afterdepolarizations in bradycardia-related acquired long-QT syndrome,” Circulation. Arrhythmia and Electrophysiology, vol. 2, no. 3, pp. 295–304, 2009.

A. Burashnikov and C. Antzelevitch, “Reinduction of atrial fibrillation immediately after termination of the arrhythmia is mediated by late phase 3 early afterdepolarization–induced triggered activity,” Circulation, vol. 107, no. 18, pp. 2355–2360, 2003.

E. Patterson, R. Lazzara, B. Szabo et al., “Sodium-calcium exchange initiated by the Ca2+ transient,” Journal of the American College of Cardiology, vol. 47, no. 6, pp. 1196–1206, 2006.

F. M. Kusumoto, M. H. Schoenfeld, C. Barrett et al., “2018 ACC/AHA/HRS guideline on the evaluation and management of patients with bradycardia and cardiac conduction delay: a report of the American College of Cardiology/American Heart Association task force on clinical practice guidelines and the Heart Rhythm Society,” Journal of the American College of Cardiology, vol. 74, no. 7, pp. e51–e156, 2019.

S. Nattel, F. Xiong, and M. Aguilar, “Demystifying rotors and their place in clinical translation of atrial fibrillation mechanisms,” Nature Reviews Cardiology, vol. 14, no. 9, pp. 509–520, 2017.

J. Jalife, O. Berenfeld, and M. Mansour, “Mother rotors and fibrillatory conduction: a mechanism of atrial fibrillation,” Cardiovascular Research, vol. 54, no. 2, pp. 204–216, 2002.

J. Kneller, R. Zou, E. J. Vigmond, Z. Wang, L. J. Leon, and S. Nattel, “Cholinergic atrial fibrillation in a computer model of a two-dimensional sheet of canine atrial cells with realistic ionic properties,” Circulation Research, vol. 90, no. 9, pp. E73–E87, 2002.

J. Kneller, J. Kalifa, R. Zou et al., “Mechanisms of atrial fibrillation termination by pure sodium channel blockade in an ionically-realistic mathematical model,” Circulation Research, vol. 96, no. 5, pp. e35–e47, 2005.