# Role of CaMKII and PKA in Early Afterdepolarization of Human Ventricular Myocardium Cell: A Computational Model Study

Authors: Ling Dai; Yunliang Zang; Dingchang Zheng; Ling Xia; Yinglan Gong
Journal: Computational and Mathematical Methods in Medicine
 (2016)

---

## Abstract

Early afterdepolarization (EAD) plays an important role in arrhythmogenesis. Many experimental studies have reported that Ca2+/calmodulin-dependent protein kinase II (CaMKII) and β-adrenergic signaling pathway are two important regulators. In this study, we developed a modified computational model of human ventricular myocyte to investigate the combined role of CaMKII and β-adrenergic signaling pathway on the occurrence of EADs. Our simulation results showed that (1) CaMKII overexpression facilitates EADs through the prolongation of late sodium current’s (INaL) deactivation progress; (2) the combined effect of CaMKII overexpression and activation of β-adrenergic signaling pathway further increases the risk of EADs, where EADs could occur at shorter cycle length (2000 ms versus 4000 ms) and lower rapid delayed rectifier K+ current (IKr) blockage (77% versus 85%). In summary, this study computationally demonstrated the combined role of CaMKII and β-adrenergic signaling pathway on the occurrence of EADs, which could be useful for searching for therapy strategies to treat EADs related arrhythmogenesis.


---

## Body

## 1. Introduction

Early afterdepolarizations (EADs) are triggered before the completion of repolarization [1] and associated with polymorphic ventricular tachyarrhythmia for long QT syndrome patients [2]. Prolongation of action potential duration (APD) and recovery of L-type Ca2+ current have been reported as two important factors for the occurrence of EADs [3]. It is also known that the increase of inward currents (e.g., ICaL and late sodium current, INaL) or the decrease of outward currents (e.g., rapid delayed rectifier K+ current, IKr and slow delayed rectifier K+ current, IKs) at plateau membrane voltage could increase the probability of EADs events. Therefore, any factors that could change the intensity or time sequence of these currents may lead to the occurrence of EADs [4–8].Ca2+/calmodulin-dependent protein kinase II (CaMKII) is a key kinase in tuning cardiac excitation-contraction coupling. Its substrates include ion channels, transporters, and accessory proteins [9]. It has been reported that CaMKII phosphorylates ICaL, leading to increased amplitude and APD prolongation and facilitating the occurrence of EADs [10, 11]. CaMKII can also alter INaL, transient outward K current (Ito), SR Ca2+-ATPase (SERCA) [12], and ryanodine receptor (RyR) channels [10]. It would therefore be useful to understand and quantify these regulatory roles. However, it is very difficult for the laboratory experiments to achieve this. Computer modelling approaches provide alternative ways, allowing us to distinguish the most effective phosphorylation target of arrhythmogenesis, which would ultimately provide useful tool in searching for antiarrhythmia therapy.It has been known thatβ-adrenergic signaling pathway regulates Ca2+ cycling partly via phosphorylation of ICaL and phospholamban (PLB) [13]. ICaL elevates intracellular Ca2+ ([Ca]i) and increases spontaneous Ca2+ release via the SERCA inhibition by PLB. The broken balance of Ca2+ cycling may contribute to the occurrence of EADs [14–16]. Volders et al. investigated the ionic mechanisms of β-adrenergic on the occurrence of EADs in canine ventricular myocyte and concluded that cellular Ca2+ overload and spontaneous SR Ca2+ release played important roles in EADs [17]. On the contrary, other studies reported that β-adrenergic agonists activate protein kinase A (PKA), which phosphorylates ICaL, RyR, PLB, SERCA, and IKs, resulting in delayed afterdepolarization (DADs) [18]. These numerous targets and different temporal characteristics of phosphorylation effects complicate the mechanism analysis of β-adrenergic agonists in relation to EADs. Recently, different computational models have been used to investigate these complex interactions. Xie et al. developed a biophysically detailed rabbit model and found that the faster time course of ICaL versus IKs increased ISO-induced transient EADs [19] and emphasized the importance of understanding the nonsteady state of kinetics in meditating β-adrenergic-induced EADs and arrhythmia. However, other targets including INaL have not been investigated thoroughly in their model. It is noted that, although the computational studies of EAD mechanisms have been widely taken, the majority of these published modelling studies have been developed based on nonhuman myocyte. Furthermore, to the best of our knowledge, there is no modelling study considering the combined effect of CaMKII overexpression and β-adrenergic agonists on EADs.This study aimed to develop a modified computational model of human ventricular myocyte that integrates CaMKII andβ-adrenergic signaling networks into a modified ORd’s dynamic model [20], with which the combined role of CaMKII and β-adrenergic signaling pathway on the occurrence of EADs would be investigated.


## 2. Methods

### 2.1. Integration of CaMKII

Since there is a lack of experimental measurements of human ventricle CaMKII pathway, O’Hara et al. used the Hund-Decker-Rudy’s dog model to describe CaMKII kinetics [21, 22]. Our model was developed from the O’Hara-Rudy dynamic model (ORd model) to integrate CaMKII pathway [20]. The following equations describe the CaMKII kinetic of human ventricle:(1)CaMKbound=CaMK0·1-CaMKtrap1+KmCaM/Ca2+ss,CaMKactive=CaMKbound+CaMKtrap,dCaMKtrapdt=αCaMK·CaMKbound·CaMKbound+CaMKtrap-βCaMK·CaMKtrap,αCaMK=0.05ms-1,βCaMK=0.00068ms-1,CaMK0=0.05,KmCaM=0.0015mM.The fraction of active CaMKII binding sites at equilibrium state (CaMK0) was set to 0.05 at the control state. CaMK0 of 0.12 was used to simulate CaMKII overexpression according to the study from Kohlhaas et al. [23].


### 2.2. Integration ofβ-Adrenergic Signaling Networks

The detailed description ofβ-adrenergic signaling networks could be found from the published study by Soltis and Saucerman [24]. PKA has been reported to phosphorylate ICaL, PLB, troponin I, RyR, myosin binding protein-C, protein phosphates Inhibitor-I [25], and IKs [24]. ICaL, PLB, and IKs were the three key factors in this study to model the inotropic effect of β-adrenergic related to potential EADs occurrence. Although the PKA phosphorylation of Na+/K+ ATPase current (INaK) has been previously described in the computational models [19, 26], its effect was not included in this study partially because there is a lack of direct measurements of INaK for the normal human ventricle [20]. The phosphorylation by PKA to the three targets (ICaL, PLB, and IKs) is described as follows:(2)favail=0.017·LCCbPKAPfracLCCbp0+0.983,(3)dss=11.0+exp⁡-V+vshift/4.230,(4)ICaL=ICaL¯·d·1-ϕICaL,CaMK·f·1-n+fCa·n·jCa+ICaL,CaMK¯·d·ϕICaL,CaMK·fCaMK·1-n+fCa,CaMK·n·jCa·favail,(5)ΦICaL,CaMK=11+Km,CaMK/CaMKactive,(6)td=1.2·0.6+1exp⁡-0.05·V+6.0+exp⁡0.09·V+14.0.Equation (2) shows the coefficient that represents the effect of PKA on ICaL phosphorylation. In (3), the value of “vshift” was increased from 3.94 (no ISO application) to 10.0 (for saturated ISO application), which means that the steady state activation curve of ICaL was moved left by 6.06 mV. The permeability of ion Ca2+ was increased by 10% with saturated ISO application. Equation (4) describes the augmentation of ICaL amplitude by multiplying “favail” with the value without PKA phosphorylation, and this represents how PKA regulates ICaL. In (6), the time constant of activation gate (td) was extended by 20% with saturated ISO application (7)fracIKsavail=0.49·IKsPKAPfracIKsp0+0.51,(8)GKs=fracIKsavail·GKs,(9)Xs05=11.60·fracIKsavail,(10)XS1,∞=11+exp⁡-V+Xs05/8.932.Equation (7) describes the factor of phosphorylation to IKs by PKA, which was used to alter maximum conductance of IKs (GKs) in (8). Meanwhile, IKs state steady activation curve was adjusted by time dependent gate value through the factor “fracIKsavail” in (9) and (10).


### 2.3. Combination of CaMKII andβ-Adrenergic Signaling Networks

(11) f P K A P L B = P L B _ P K A n f r a c P K A P L B 0 · 1 4 + 3 4 , (12) f J u p p = 1.0 1.0 + K m , C a M K · f P K A P L B / C a M K a c t i v e , (13) J u p = 1.0 - f J u p p · J u p n p + f J u p p · J u p p - J l e a k .CaM K active that is affected by CaMK0 as shown in (1) would influence ICaL in (4) via the fraction of ICaL channels phosphorylated by CaMKII (ΦICaL,CaMK). The fraction of SERCAs phosphorylated by CaMKII in (12) was affected by “fPKAPLB” and “CaMKactive”, representing the effects of PKA and CaMKII on SERCAs, respectively. As shown in (13), SERCAs were separated into nonphosphorylated populations and CaMKII phosphorylated populations. Therefore, the total Ca2+ uptake via SERCAs was adjusted by these two networks simultaneously.


### 2.4. Simulation Strategy

In order to determine the potential targets in CaMKII-induced EADs, CaMKII was solely overexpressed by assigning the CaMK0of 0.12 to specific targets (including INaL, ICaL, ICaK, INaCa, and Ito), respectively, while β-adrenergic signaling was maintained inactive. Next, CaMKII overexpression was applied with ISO administration. 1 μM ISO was applied in this study to simulate its effect on myocyte action potential and ion currents.The cycle length of 2000 ms was used in this study since it was closer to normal human beat rhythm than the length of 4000 ms used in the experiments from Guo et al. [27]. In Guo et al.’s work, EADs appeared with IKr blockage of about 85% [27], which was used for comparison in this study. 500 cycles were performed when the simulation reached steady state. In each individual cycle length, the upper bound of solver step size was set 2 ms.


## 2.1. Integration of CaMKII

Since there is a lack of experimental measurements of human ventricle CaMKII pathway, O’Hara et al. used the Hund-Decker-Rudy’s dog model to describe CaMKII kinetics [21, 22]. Our model was developed from the O’Hara-Rudy dynamic model (ORd model) to integrate CaMKII pathway [20]. The following equations describe the CaMKII kinetic of human ventricle:(1)CaMKbound=CaMK0·1-CaMKtrap1+KmCaM/Ca2+ss,CaMKactive=CaMKbound+CaMKtrap,dCaMKtrapdt=αCaMK·CaMKbound·CaMKbound+CaMKtrap-βCaMK·CaMKtrap,αCaMK=0.05ms-1,βCaMK=0.00068ms-1,CaMK0=0.05,KmCaM=0.0015mM.The fraction of active CaMKII binding sites at equilibrium state (CaMK0) was set to 0.05 at the control state. CaMK0 of 0.12 was used to simulate CaMKII overexpression according to the study from Kohlhaas et al. [23].


## 2.2. Integration ofβ-Adrenergic Signaling Networks

The detailed description ofβ-adrenergic signaling networks could be found from the published study by Soltis and Saucerman [24]. PKA has been reported to phosphorylate ICaL, PLB, troponin I, RyR, myosin binding protein-C, protein phosphates Inhibitor-I [25], and IKs [24]. ICaL, PLB, and IKs were the three key factors in this study to model the inotropic effect of β-adrenergic related to potential EADs occurrence. Although the PKA phosphorylation of Na+/K+ ATPase current (INaK) has been previously described in the computational models [19, 26], its effect was not included in this study partially because there is a lack of direct measurements of INaK for the normal human ventricle [20]. The phosphorylation by PKA to the three targets (ICaL, PLB, and IKs) is described as follows:(2)favail=0.017·LCCbPKAPfracLCCbp0+0.983,(3)dss=11.0+exp⁡-V+vshift/4.230,(4)ICaL=ICaL¯·d·1-ϕICaL,CaMK·f·1-n+fCa·n·jCa+ICaL,CaMK¯·d·ϕICaL,CaMK·fCaMK·1-n+fCa,CaMK·n·jCa·favail,(5)ΦICaL,CaMK=11+Km,CaMK/CaMKactive,(6)td=1.2·0.6+1exp⁡-0.05·V+6.0+exp⁡0.09·V+14.0.Equation (2) shows the coefficient that represents the effect of PKA on ICaL phosphorylation. In (3), the value of “vshift” was increased from 3.94 (no ISO application) to 10.0 (for saturated ISO application), which means that the steady state activation curve of ICaL was moved left by 6.06 mV. The permeability of ion Ca2+ was increased by 10% with saturated ISO application. Equation (4) describes the augmentation of ICaL amplitude by multiplying “favail” with the value without PKA phosphorylation, and this represents how PKA regulates ICaL. In (6), the time constant of activation gate (td) was extended by 20% with saturated ISO application (7)fracIKsavail=0.49·IKsPKAPfracIKsp0+0.51,(8)GKs=fracIKsavail·GKs,(9)Xs05=11.60·fracIKsavail,(10)XS1,∞=11+exp⁡-V+Xs05/8.932.Equation (7) describes the factor of phosphorylation to IKs by PKA, which was used to alter maximum conductance of IKs (GKs) in (8). Meanwhile, IKs state steady activation curve was adjusted by time dependent gate value through the factor “fracIKsavail” in (9) and (10).


## 2.3. Combination of CaMKII andβ-Adrenergic Signaling Networks

(11) f P K A P L B = P L B _ P K A n f r a c P K A P L B 0 · 1 4 + 3 4 , (12) f J u p p = 1.0 1.0 + K m , C a M K · f P K A P L B / C a M K a c t i v e , (13) J u p = 1.0 - f J u p p · J u p n p + f J u p p · J u p p - J l e a k .CaM K active that is affected by CaMK0 as shown in (1) would influence ICaL in (4) via the fraction of ICaL channels phosphorylated by CaMKII (ΦICaL,CaMK). The fraction of SERCAs phosphorylated by CaMKII in (12) was affected by “fPKAPLB” and “CaMKactive”, representing the effects of PKA and CaMKII on SERCAs, respectively. As shown in (13), SERCAs were separated into nonphosphorylated populations and CaMKII phosphorylated populations. Therefore, the total Ca2+ uptake via SERCAs was adjusted by these two networks simultaneously.


## 2.4. Simulation Strategy

In order to determine the potential targets in CaMKII-induced EADs, CaMKII was solely overexpressed by assigning the CaMK0of 0.12 to specific targets (including INaL, ICaL, ICaK, INaCa, and Ito), respectively, while β-adrenergic signaling was maintained inactive. Next, CaMKII overexpression was applied with ISO administration. 1 μM ISO was applied in this study to simulate its effect on myocyte action potential and ion currents.The cycle length of 2000 ms was used in this study since it was closer to normal human beat rhythm than the length of 4000 ms used in the experiments from Guo et al. [27]. In Guo et al.’s work, EADs appeared with IKr blockage of about 85% [27], which was used for comparison in this study. 500 cycles were performed when the simulation reached steady state. In each individual cycle length, the upper bound of solver step size was set 2 ms.


## 3. Results

### 3.1. Effect of CaMKII Overexpression on Ion Currents

#### 3.1.1. Late Sodium Current (INaL)

As shown in Table1 and Figure 1, with the overexpressed CaMK0 value to INaL of 0.12 and the normal CaMK0 value of 0.05 to other targets, the alternated EADs occur from our simulation with the cycle length (CL) of 2000 ms and IKr blockage of 85%.Table 1
Current increment and EADs occurrence when different targets were phosphorylated by CaMKII independently.
Cycle length (ms) I K r blockage level (%) CaMKII target (CaMK0 = 0.12) Target current variation Deactivation time (ms) EADs 2000 85 I N a L Decreased by 24%(from 0.25 to 0.19 µA/µF) 1297 versus 761 ms Alternated I C a L Increased by 4.2%(from 1.67 to 1.74 µA/µF)  No I C a K Increased by 6.5%(from 0.62 to 0.66 µA/µF)  No I C a N a Increased by 8.1%(from 0.37 to 0.40 µA/µF)  No I t o Increased by 3.2%(from 0.95 to 0.98 µA/µF)  NoFigure 1
(a) No EAD was produced when CaMKII phosphorylation level was in control (CaMK0 = 0.05). (b) Corresponding INaL when no EADs occurred in (a). (c) Alternated EADs were produced when INaL phosphorylation by CaMKII was enhanced with CaMK0 of 0.12 and other targets’ phosphorylation levels were in control (CaMK0 = 0.05). (d) Corresponding INaL when alternated EADs occurred in (c). Under these conditions, CL was 2000 ms and IKr was blocked by 85%.As shown in Figure1, INaL amplitudes alternated with overexpressed CaMK0. In beats with EADs, INaL amplitude was about 24% smaller (0.19 μA/μF versus 0.25 μA/μF) than these with normal CaMK0 value to INaL and in beats without EADs, INaL amplitude was also reduced by 32% (0.17 μA/μF versus 0.25 μA/μF). In Figure 1(d), when EADs occurred, INaL deactivated, which was about 168% of normal beats (1297 ms versus 773 ms) and about 170% of the control situation in Figure 1(b) (1297 ms versus 761 ms), indicating that overexpressed CaMKII phosphorylation level of INaL reduces its amplitude and prolongs INaL deactivation process. Furthermore, the results also indicate that the delayed deactivation of INaL, rather than its amplitude variation, contributes to the formation of EADs.


#### 3.1.2. L-Type Calcium Current (ICaL)

Similarly, as shown in Table1, with the overexpressed CaMK0 to ICaL of 0.12 and the normal CaMK0 value of 0.05 to other targets, no EADs occur at CL of 2000 ms and IKr blockage of 85%, suggesting that the probability of EADs might have little relation with amplitude of ICaL. With the enhanced CaMK0 value to ICaL, the amplitudes of ICaL only increased by 4.2% (1.74 μA/μF versus 1.67 μA/μF).


#### 3.1.3. K+ Current through the L-Type Ca2+ Channel (ICaK), Na+ Current through the L-Type Ca2+ Channel (ICaNa), and Transient Outward K+ Current (Ito)

As above, CL was set to 2000 ms andIKr was blocked by 85%. As shown in Table 1, with enhanced CaMKII phosphorylation level (CaMK0 = 0.12) to different targets, ICaK increased only by 6.5% (0.66 μA/μF versus 0.62 μA/μF), ICaNa by 8.1% (0.40 μA/μF versus 0.37 μA/μF), and Ito by 3.2% (0.98 μA/μF versus 0.95 μA/μF). No EADs occurred under all these conditions.


### 3.2. Combined Effect of CaMKII Overexpression andβ-Adrenergic Agonist

#### 3.2.1. Normal CaMKII and 1μM ISO

The action potentials with differentIKr blockage level are shown in Figure 2. A fixed CL =2000 ms was used and CaMKII phosphorylation level to all targets was kept control (CaMK0 = 0.05). In Figure 2(a), with IKr blocked by 85%, stable EADs occurred with the application of ISO, indicating that β-adrenergic agonist facilitates EADs. The IKr blockage level decreased gradually from 85%, and, when it was decreased to 77%, as shown in Figure 2(b), these EADs disappeared, indicating that 77% was the threshold value for EADs disappearance in this setting. Therefore, 77% was used in following simulation to compare with previously published level of 85% [20]. These results are listed from row 2 to row 3 in Table 2.Table 2
Combined effect of CaMKII overexpression andβ-adrenergic agonist on EADs.
Cycle length (ms) ISO application CaMKII target (CaMK0 = 0.12) I K r blockage level (%) EADs 2000 None None 85 None 1µM None 85 Yes None 77 No I N a L 85 Yes I N a L 77 Yes I C a L 85 Yes I C a L 77 NoFigure 2
(a) EADs were induced when 1μM ISO was applied and IKr was blocked by 85%. In (b), EADs disappeared when IKr blockage was reduced to 77%. Cycle length was set 2000 ms and CaMKII phosphorylation level to all targets was kept control (CaMK0 = 0.05).
(a) (b)


#### 3.2.2. Enhanced CaMKII toINaL and 1 μM ISO

As shown in Figure3, with enhanced CaMKII to INaL and 1 μM ISO, when IKr was blocked by 85%, EADs were induced. The EADs were still observed when IKr was blocked by 77%, suggesting that INaL phosphorylation by CaMKII and ISO application together increased the probability of EADs. These results are listed from row 4 to row 5 in Table 2.Figure 3
EADs occurred whenIKr was blocked by 85% (a) and 77% (c); corresponding INaL when IKr was blocked by 85% (b) and 77% (d). CL was set 2000 ms, 1 μM ISO was applied and CaMK0 for INaL was set 0.12 independently.


#### 3.2.3. Enhanced CaMKII toICaL and 1 μM ISO

As shown in Figure4, with enhanced CaMKII to ICaL and 1 μM ISO, when IKr was blocked by 85%, EADs were induced, but when IKr blockage was reduced to 77%, EADs disappeared. These results are listed from row 6 to row 7 in Table 2.Figure 4
1μM ISO was applied and cycle length was 2000 ms. CaMK0 for ICaL was 0.12 but CaMK0 for other targets was 0.05. EADs occurred when IKr was blocked by 85% in (a), but when IKr blockage was reduced to 77%, EADs vanished in (c). (b) ICaL figures when EADs existed. (d) ICaL figures when EADs vanished.


## 3.1. Effect of CaMKII Overexpression on Ion Currents

### 3.1.1. Late Sodium Current (INaL)

As shown in Table1 and Figure 1, with the overexpressed CaMK0 value to INaL of 0.12 and the normal CaMK0 value of 0.05 to other targets, the alternated EADs occur from our simulation with the cycle length (CL) of 2000 ms and IKr blockage of 85%.Table 1
Current increment and EADs occurrence when different targets were phosphorylated by CaMKII independently.
Cycle length (ms) I K r blockage level (%) CaMKII target (CaMK0 = 0.12) Target current variation Deactivation time (ms) EADs 2000 85 I N a L Decreased by 24%(from 0.25 to 0.19 µA/µF) 1297 versus 761 ms Alternated I C a L Increased by 4.2%(from 1.67 to 1.74 µA/µF)  No I C a K Increased by 6.5%(from 0.62 to 0.66 µA/µF)  No I C a N a Increased by 8.1%(from 0.37 to 0.40 µA/µF)  No I t o Increased by 3.2%(from 0.95 to 0.98 µA/µF)  NoFigure 1
(a) No EAD was produced when CaMKII phosphorylation level was in control (CaMK0 = 0.05). (b) Corresponding INaL when no EADs occurred in (a). (c) Alternated EADs were produced when INaL phosphorylation by CaMKII was enhanced with CaMK0 of 0.12 and other targets’ phosphorylation levels were in control (CaMK0 = 0.05). (d) Corresponding INaL when alternated EADs occurred in (c). Under these conditions, CL was 2000 ms and IKr was blocked by 85%.As shown in Figure1, INaL amplitudes alternated with overexpressed CaMK0. In beats with EADs, INaL amplitude was about 24% smaller (0.19 μA/μF versus 0.25 μA/μF) than these with normal CaMK0 value to INaL and in beats without EADs, INaL amplitude was also reduced by 32% (0.17 μA/μF versus 0.25 μA/μF). In Figure 1(d), when EADs occurred, INaL deactivated, which was about 168% of normal beats (1297 ms versus 773 ms) and about 170% of the control situation in Figure 1(b) (1297 ms versus 761 ms), indicating that overexpressed CaMKII phosphorylation level of INaL reduces its amplitude and prolongs INaL deactivation process. Furthermore, the results also indicate that the delayed deactivation of INaL, rather than its amplitude variation, contributes to the formation of EADs.


### 3.1.2. L-Type Calcium Current (ICaL)

Similarly, as shown in Table1, with the overexpressed CaMK0 to ICaL of 0.12 and the normal CaMK0 value of 0.05 to other targets, no EADs occur at CL of 2000 ms and IKr blockage of 85%, suggesting that the probability of EADs might have little relation with amplitude of ICaL. With the enhanced CaMK0 value to ICaL, the amplitudes of ICaL only increased by 4.2% (1.74 μA/μF versus 1.67 μA/μF).


### 3.1.3. K+ Current through the L-Type Ca2+ Channel (ICaK), Na+ Current through the L-Type Ca2+ Channel (ICaNa), and Transient Outward K+ Current (Ito)

As above, CL was set to 2000 ms andIKr was blocked by 85%. As shown in Table 1, with enhanced CaMKII phosphorylation level (CaMK0 = 0.12) to different targets, ICaK increased only by 6.5% (0.66 μA/μF versus 0.62 μA/μF), ICaNa by 8.1% (0.40 μA/μF versus 0.37 μA/μF), and Ito by 3.2% (0.98 μA/μF versus 0.95 μA/μF). No EADs occurred under all these conditions.


## 3.1.1. Late Sodium Current (INaL)

As shown in Table1 and Figure 1, with the overexpressed CaMK0 value to INaL of 0.12 and the normal CaMK0 value of 0.05 to other targets, the alternated EADs occur from our simulation with the cycle length (CL) of 2000 ms and IKr blockage of 85%.Table 1
Current increment and EADs occurrence when different targets were phosphorylated by CaMKII independently.
Cycle length (ms) I K r blockage level (%) CaMKII target (CaMK0 = 0.12) Target current variation Deactivation time (ms) EADs 2000 85 I N a L Decreased by 24%(from 0.25 to 0.19 µA/µF) 1297 versus 761 ms Alternated I C a L Increased by 4.2%(from 1.67 to 1.74 µA/µF)  No I C a K Increased by 6.5%(from 0.62 to 0.66 µA/µF)  No I C a N a Increased by 8.1%(from 0.37 to 0.40 µA/µF)  No I t o Increased by 3.2%(from 0.95 to 0.98 µA/µF)  NoFigure 1
(a) No EAD was produced when CaMKII phosphorylation level was in control (CaMK0 = 0.05). (b) Corresponding INaL when no EADs occurred in (a). (c) Alternated EADs were produced when INaL phosphorylation by CaMKII was enhanced with CaMK0 of 0.12 and other targets’ phosphorylation levels were in control (CaMK0 = 0.05). (d) Corresponding INaL when alternated EADs occurred in (c). Under these conditions, CL was 2000 ms and IKr was blocked by 85%.As shown in Figure1, INaL amplitudes alternated with overexpressed CaMK0. In beats with EADs, INaL amplitude was about 24% smaller (0.19 μA/μF versus 0.25 μA/μF) than these with normal CaMK0 value to INaL and in beats without EADs, INaL amplitude was also reduced by 32% (0.17 μA/μF versus 0.25 μA/μF). In Figure 1(d), when EADs occurred, INaL deactivated, which was about 168% of normal beats (1297 ms versus 773 ms) and about 170% of the control situation in Figure 1(b) (1297 ms versus 761 ms), indicating that overexpressed CaMKII phosphorylation level of INaL reduces its amplitude and prolongs INaL deactivation process. Furthermore, the results also indicate that the delayed deactivation of INaL, rather than its amplitude variation, contributes to the formation of EADs.


## 3.1.2. L-Type Calcium Current (ICaL)

Similarly, as shown in Table1, with the overexpressed CaMK0 to ICaL of 0.12 and the normal CaMK0 value of 0.05 to other targets, no EADs occur at CL of 2000 ms and IKr blockage of 85%, suggesting that the probability of EADs might have little relation with amplitude of ICaL. With the enhanced CaMK0 value to ICaL, the amplitudes of ICaL only increased by 4.2% (1.74 μA/μF versus 1.67 μA/μF).


## 3.1.3. K+ Current through the L-Type Ca2+ Channel (ICaK), Na+ Current through the L-Type Ca2+ Channel (ICaNa), and Transient Outward K+ Current (Ito)

As above, CL was set to 2000 ms andIKr was blocked by 85%. As shown in Table 1, with enhanced CaMKII phosphorylation level (CaMK0 = 0.12) to different targets, ICaK increased only by 6.5% (0.66 μA/μF versus 0.62 μA/μF), ICaNa by 8.1% (0.40 μA/μF versus 0.37 μA/μF), and Ito by 3.2% (0.98 μA/μF versus 0.95 μA/μF). No EADs occurred under all these conditions.


## 3.2. Combined Effect of CaMKII Overexpression andβ-Adrenergic Agonist

### 3.2.1. Normal CaMKII and 1μM ISO

The action potentials with differentIKr blockage level are shown in Figure 2. A fixed CL =2000 ms was used and CaMKII phosphorylation level to all targets was kept control (CaMK0 = 0.05). In Figure 2(a), with IKr blocked by 85%, stable EADs occurred with the application of ISO, indicating that β-adrenergic agonist facilitates EADs. The IKr blockage level decreased gradually from 85%, and, when it was decreased to 77%, as shown in Figure 2(b), these EADs disappeared, indicating that 77% was the threshold value for EADs disappearance in this setting. Therefore, 77% was used in following simulation to compare with previously published level of 85% [20]. These results are listed from row 2 to row 3 in Table 2.Table 2
Combined effect of CaMKII overexpression andβ-adrenergic agonist on EADs.
Cycle length (ms) ISO application CaMKII target (CaMK0 = 0.12) I K r blockage level (%) EADs 2000 None None 85 None 1µM None 85 Yes None 77 No I N a L 85 Yes I N a L 77 Yes I C a L 85 Yes I C a L 77 NoFigure 2
(a) EADs were induced when 1μM ISO was applied and IKr was blocked by 85%. In (b), EADs disappeared when IKr blockage was reduced to 77%. Cycle length was set 2000 ms and CaMKII phosphorylation level to all targets was kept control (CaMK0 = 0.05).
(a) (b)


### 3.2.2. Enhanced CaMKII toINaL and 1 μM ISO

As shown in Figure3, with enhanced CaMKII to INaL and 1 μM ISO, when IKr was blocked by 85%, EADs were induced. The EADs were still observed when IKr was blocked by 77%, suggesting that INaL phosphorylation by CaMKII and ISO application together increased the probability of EADs. These results are listed from row 4 to row 5 in Table 2.Figure 3
EADs occurred whenIKr was blocked by 85% (a) and 77% (c); corresponding INaL when IKr was blocked by 85% (b) and 77% (d). CL was set 2000 ms, 1 μM ISO was applied and CaMK0 for INaL was set 0.12 independently.


### 3.2.3. Enhanced CaMKII toICaL and 1 μM ISO

As shown in Figure4, with enhanced CaMKII to ICaL and 1 μM ISO, when IKr was blocked by 85%, EADs were induced, but when IKr blockage was reduced to 77%, EADs disappeared. These results are listed from row 6 to row 7 in Table 2.Figure 4
1μM ISO was applied and cycle length was 2000 ms. CaMK0 for ICaL was 0.12 but CaMK0 for other targets was 0.05. EADs occurred when IKr was blocked by 85% in (a), but when IKr blockage was reduced to 77%, EADs vanished in (c). (b) ICaL figures when EADs existed. (d) ICaL figures when EADs vanished.


## 3.2.1. Normal CaMKII and 1μM ISO

The action potentials with differentIKr blockage level are shown in Figure 2. A fixed CL =2000 ms was used and CaMKII phosphorylation level to all targets was kept control (CaMK0 = 0.05). In Figure 2(a), with IKr blocked by 85%, stable EADs occurred with the application of ISO, indicating that β-adrenergic agonist facilitates EADs. The IKr blockage level decreased gradually from 85%, and, when it was decreased to 77%, as shown in Figure 2(b), these EADs disappeared, indicating that 77% was the threshold value for EADs disappearance in this setting. Therefore, 77% was used in following simulation to compare with previously published level of 85% [20]. These results are listed from row 2 to row 3 in Table 2.Table 2
Combined effect of CaMKII overexpression andβ-adrenergic agonist on EADs.
Cycle length (ms) ISO application CaMKII target (CaMK0 = 0.12) I K r blockage level (%) EADs 2000 None None 85 None 1µM None 85 Yes None 77 No I N a L 85 Yes I N a L 77 Yes I C a L 85 Yes I C a L 77 NoFigure 2
(a) EADs were induced when 1μM ISO was applied and IKr was blocked by 85%. In (b), EADs disappeared when IKr blockage was reduced to 77%. Cycle length was set 2000 ms and CaMKII phosphorylation level to all targets was kept control (CaMK0 = 0.05).
(a) (b)


## 3.2.2. Enhanced CaMKII toINaL and 1 μM ISO

As shown in Figure3, with enhanced CaMKII to INaL and 1 μM ISO, when IKr was blocked by 85%, EADs were induced. The EADs were still observed when IKr was blocked by 77%, suggesting that INaL phosphorylation by CaMKII and ISO application together increased the probability of EADs. These results are listed from row 4 to row 5 in Table 2.Figure 3
EADs occurred whenIKr was blocked by 85% (a) and 77% (c); corresponding INaL when IKr was blocked by 85% (b) and 77% (d). CL was set 2000 ms, 1 μM ISO was applied and CaMK0 for INaL was set 0.12 independently.


## 3.2.3. Enhanced CaMKII toICaL and 1 μM ISO

As shown in Figure4, with enhanced CaMKII to ICaL and 1 μM ISO, when IKr was blocked by 85%, EADs were induced, but when IKr blockage was reduced to 77%, EADs disappeared. These results are listed from row 6 to row 7 in Table 2.Figure 4
1μM ISO was applied and cycle length was 2000 ms. CaMK0 for ICaL was 0.12 but CaMK0 for other targets was 0.05. EADs occurred when IKr was blocked by 85% in (a), but when IKr blockage was reduced to 77%, EADs vanished in (c). (b) ICaL figures when EADs existed. (d) ICaL figures when EADs vanished.


## 4. Discussion

This study developed a modified computational model of human ventricular myocardium cell based on the ORd human model with the integration of regulation mechanism by CaMKII and PKA [20], with which their effects on EADs have been investigated.EADs often occur during bradycardia under the condition of reduced repolarization reserve. O’Hara’s group successfully elicited EADs with cycle length reduced to 4000 ms andIKr blocked by 85% [20]. In our simulation, with the cycle length halved (2000 ms) and IKr blocked by 85%, alternated EADs occurred with the sole effect of CaMKII overexpression to INaL, suggesting that ventricular myocardial cell with CaMKII overexpression is more susceptible to EADs in normal HR range (CL = 2000 ms). Additionally, previous work has shown that ICaL plays an important role in the occurrence of EADs [6]. Our simulation showed that CaMKII overexpression slightly increased ICaL amplitude, but this effect alone did not induce EADs, and CaMKII did not alter the overlap region of ICaL steady state activation and reactivation curves. Zaza et al. reported that [28], when repolarization is suitably slow, channel reactivation within the overlap region may break the current balance and support the possibility of autoregenerative depolarization. Other inward currents such as Na+-Ca2+ exchange current (INaCa) may play certain roles in triggering EADs if amplitude is augmented and falls into the overlap region. Our simulation results suggest that CaMKII enhances these currents but does not induce EADs via this effect individually. Therefore, when CaMKII is overexpressed, the susceptibility to EADs is mainly originated from INaL variation.Our study also demonstrated the behavior of ventricular myocardial cell with the integration of CaMKII overexpression and ISO application. With the ISO application, cyclic AMP (cAMP) is formed throughβ-adrenergic mediated activation of adenylyl cyclase, which activates PKA, a well-described mediator with targets that promote myocardial performance. In the case where PKA took effect independently, EADs occurred with shorter cycle length (2000 ms versus 4000 ms), indicating that the precondition of EADs is relaxed. Our results were different from that from Xie et al.’s study, where the APD shortening after ISO application was observed without EADs in steady state [19]. Xie et al.’s results could be caused by the simulation of transient ICaL recovery and the prevented spontaneous SR Ca2+ release, limiting the INaCa in forward mode. Therefore, when the cells step into steady state with small inward INaCa, the shortening of APD could be reasonable. In our work, SR Ca2+ release was enhanced by ICaL amplification. With the integration with IKr blockage, INaCa was more likely to work as inward current, contributing to the prolongation of APD and occurrence of EADs. Additionally, the shortening of APD was also obtained without IKr blockage when β-adrenergic pathway was activated. Simulation results showed that when cycle length was chosen at 500 ms, 1000 ms, 1500 ms, 2000 ms, and 4000 ms, the corresponding APD90 shortening was 31.6 ms (233.0 versus 201.4 ms), 30.5 ms (268.6 versus 238.1 ms), 22.3 ms (282.6 versus 260.3 ms), 17.3 ms (289.1 versus 271.8 ms), and 18.6 ms (305.6 versus 287.0 ms). These results were consistent with the experiment measurements from Volders et al. [29]. In the case that CaMKII overexpression acted on INaL alone with ISO application, EADs occur at CL = 2000 ms and IKr blockage by 77%, suggesting that the combination of PKA and CaMKII overexpression on INaL relaxes the precondition further (IKr blockage 77% versus 85%) and increases the probability of EADs. With the CaMKII overexpressed on ICaL alone, PKA induced EADs occurrence with 85%  IKr blockage, not 77%  IKr blockage, suggesting that, with the ISO application, INaL phosphorylation by CaMKII has more effect on EADs than ICaL. The steady state of CaMKactive was simulated with and without ISO application when CL = 2000 ms, and their corresponding maximum values were 0.0469 and 0.0421, respectively. This suggests that CaMKII activation increased with the application of ISO.Our simulation results have shown that our proposed model is useful for exploring the interaction ofβ-adrenergic receptor signaling and CaMKII in formation of EADs; some potential limitations need to be addressed. Firstly, the role of CaMKII on RyR function has not been considered. There is growing evidence about its role in modulating RyR function. Increasing Ca2+ leak via RyR would limit SR Ca2+ load and disturb Ca2+ cycling. This should be incorporated into an improved model in a future study. Secondly, CaMKII overexpression may downregulate inward rectifier potassium current (IK1) and increase baseline IK1 amplitude [30]. We did not include any CaMKII effect on IK1 in our model. Thirdly, our model only describes the short-term effect of CaMKII on different current targets. It has been suggested that short-term and chronic effects of CaMKII are different [30–33]. Short-term (milliseconds to hours) CaMKII overexpression may slow Ito inactivation and accelerate recovery, but chronic overexpression may downregulate Ito,fast and upregulate Ito,slow [34]. Therefore, there is a scope to improve our model by taking chronic effects of CaMKII into consideration in future. Fourthly, there is a lack of an accurate measurement about how ISO regulates INaK with dynamic calcium change during the beat cycle, although it has been published by Gao et al. [35] that ISO regulated INaK increased (when calcium is 1.4 μM) or decreased in pump current (when calcium is 0.15 μM) depending on the intracellular Ca2+ concentrations. In our model the resting calcium concentration was 0.12 μM and the peak calcium concentration was 1.1 μM. Similar to the simulation work of Heijman et al. [26], by simply increasing the pump current from 17% to 33%, no significant change of INaK in triggering EADs (results not shown here) has been observed. However, when the continuous experimental data from human ventricular cells is available, the regulation of ISO on INaK could be incorporated into the model to reconfirm this nonsignificant effect. Fifthly, it has been reported that local regulation of cAMP and substrate phosphorylation play important roles in β-adrenergic receptor signaling [26, 36–38], so it could be useful to incorporate local control mechanism in a future study. Lastly, there are four CaMKII isoforms (α,β,γ,δ) with different distributions, kinetics, and roles in physiological and pathological adjustments. Until now, these differences have not been fully understood. Developing a model with detailed CaMKII isoforms information could be useful when experiment data are available.In conclusion, our simulation results computationally demonstrated a better understanding of the combinational effect of CaMKII and ISO stimulus on the occurrence of EADs in human ventricular myocyte, which may provide useful tool to research therapeutic methods for the treatment of arrhythmia.



---

*Source: 4576313-2016-12-08.xml*