# Isoprenaline: A Potential Contributor in Sick Sinus Syndrome—Insights from a Mathematical Model of the Rabbit Sinoatrial Node

Authors: Xiang Li; Ji-qian Zhang; Jian-wei Shuai
Journal: The Scientific World Journal
 (2014)

---

## Abstract

The mechanism of isoprenaline exerting its effects on cardiac pacemaking and driving in sick sinus syndrome is controversial and unresolved. In this paper, mathematical models for rabbit sinoatrial node cells were modified by incorporating equations for the known dose-dependent actions of isoprenaline on various ionic channel currents, the intracellular Ca2+ transient, and i
Na changes induced by SCN5A gene mutations; the cell models were also incorporated into an intact SAN-atrium model of the rabbit heart that is based on both heterogeneities of the SAN electrophysiology and histological structure. Our results show that, in both central and peripheral cell models, isoprenaline could not only shorten the action potential duration, but also increase the amplitude of action potential. The mutation impaired the SAN pacemaking. Simulated vagal nerve activity amplified the bradycardic effects of the mutation. However, in tissue case, the pacemaker activity may show temporal, spatial, or even spatiotemporal cessation caused by the mutation. Addition of isoprenaline could significantly diminish the bradycardic effect of the mutation and the SAN could restart pacing and driving the surrounding tissue. Positive effects of isoprenaline may primarily be attributable to an increase in i
Na and i
Ca,
T which were reduced by the mutation.


---

## Body

## 1. Introduction

Sick sinus syndrome (SSS) comprises a variety of conditions involving sinus node dysfunction (SND) which occurs as a result of anatomical damage to the sinoatrial node (SAN) of the heart. Abnormalities encompassed in this syndrome include sinus bradycardia, sinus arrest or exit block, combinations of sinoatrial and atrioventricular nodal conduction disturbances, and atrial tachyarrhythmias [1]. Patients may experience syncope, presyncope, palpitations, or dizziness who have to be fitted with an electronic pacemaker [1, 2]. While the syndrome can occur in elderly and pediatric patients, it can also occur in healthy people without any evident structural heart disease but with genetic defects. Recent studies have shown that familial SSS is linked to loss-of-function mutations of the SCN5A gene [3], which result in decreased inward Na+ current [4, 5].Currently, no medications are routinely used to treat symptomatic SND. However, acute treatment with the anticholinergic agent atropine and the adrenergic agonist isoprenaline (Iso) may be warranted. Studies have shown thatβ-adrenergic stimulation could increase the heart rate through accelerating the spontaneous activity of the pacemaker of the heart [6]. It is also believed that this occurs through β-adrenoceptor mediated modulation of ionic currents that contribute to pacemaker activity [7]. Iso is a sympathomimetic β-adrenergic agonist medication and known to accelerate the sinus node and enhance AV nodal conduction [8]. Experiments have shown that Iso increases the L-type calcium current [9–11] and the delayed rectifier potassium current [10, 11] and shifts voltage-dependent i
f activation curve to more positive potentials [12–15]. Furthermore, a number of clinical studies have shown that Iso exerts positive inotropic and chronotropic effects on the heart.Whereas, to our knowledge, many researchers have investigated the effect of Iso mainly at the cellular level, its effect on the abnormal activity at the tissue level has been rarely studied before, particularly for the abnormal activity associated with sick sinus syndrome. Thus, an intriguing problem is whether patients with SSS could recover with the help of Iso and do not need cardiac pacing. In this study, we modified models of SAN to investigate theoretically the effect of Iso on the SSS and discussed the inherent mechanism between them.


## 2. Methods

### 2.1. Model Description

In order to address the above issues, mathematical model of the rabbit central and peripheral SAN cells proposed by Kurata et al. [16] was used in this paper. The action potentials generated by using the models of a central cell and a peripheral cell were shown in Figure 1(c). Besides, a 2D anatomic model of the intact SAN-atrium tissue developed by Butters et al. [17] was also used to study the effects of Iso on SAN pacemaking and driving. This 2D intact SAN-atrium tissue model was based on histologically reconstructed geometry of a single slice of the rabbit right atrium, which was cut through the atrial muscle of the crista terminalis (CT) and the intercaval region with central and peripheral SAN areas. And an experimentally observed [18, 19] nonconductive region (block zone) next to the SAN toward the atrial septum was also incorporated in this model (see Figure 1(a)). Figure 1(b) shows the initiation and conduction of pacemaker activity computed from the 2D model. Snapshots of membrane potentials across the 2D lattice were shown at 50 ms (Figure 1(b)(i)), 70 ms (Figure 1(b)(ii)), and 90 ms (Figure 1(b)(iii)) after the initial conduction. The pacemaking action was first initiated in the center (Figure 1(b)(i)) and then propagated toward the surrounding tissues (Figure 1(b)(iii)).Figure 1
Model of the SAN and surrounding atrial tissue. (a) Color-coded distribution of cell types throughout the 2D tissue slice. (b) Snapshots of initiation and conduction in 2D anatomical model at various times after initial configuration; (i) 50 ms; (ii) 70 ms; (iii) 90 ms. (c) Action potentials in one cell of the central and peripheral SAN. (d) Action potential profiles during conduction through the slice. As, atrial septum; SA-C, SAN center; AM, atrial muscle (CT).The mechanism underlying the genesis of SND (clinically known as SSS) is believed to be associated with changes in the intrinsic properties of the SAN with age [20]. Intercellular electrical coupling strength (D) is considered to reflect the changes in the intrinsic properties of the SAN. To incorporate the effects of SCN5A mutation and acetylcholine (Ach) into the models, (i) parameters of the Na+ current were changed according to the study of Butters et al. [17] and (ii) an Ach-activated K+ current (i
K
,
Ach), the inhibited L-type Ca2+ current (i
Ca
,
L), and the shift in hyperpolarization-activated current (i
f) activation curve were introduced.The model equations used for the cells, 2D tissue, modified Na+ channel current, and the channel currents relative to Ach are presented in Table 1. A detailed description of parameters of these models are listed in Table 3 and can be found in [16, 19–21].Table 1
Model equations of the SAN single cell, 2D tissue, and  some channel currents.
Single cell model of SAN
d
V
/
d
t
=
-
(
i
Na
+
i
Ca,
L
+
i
Ca,
T
+
i
K,
r
+
i
K,
s
+
i
st
+
i
to
+
i
sus
+
i
b
,Na
+
i
b
,K
+
i
NaK
+
i
NaCa
)
/
C
m
2D intact SAN-atrium tissue model
∂
V
∂
t
=
∇
·
(
D
∇
V
)
-
i
tot
C
m
Changes of sodium current induced by mutation
i
Na
=
(
S
Na
×
g
Na
)
m
3
h
[
Na
+
]
o
F
2
R
T
e
(
V
-
E
Na
)
F
/
R
T
-
1
e
V
F
/
R
T
-
1
V
τ
h
1
=
S
τ
h
1
(
3.717
×
1
0
-
6
e
-
0.2815
(
V
+
17.11
)
1
+
3.732
×
1
0
-
3
e
-
0.3426
(
V
+
37.76
)
+
5.977
×
1
0
-
4
)
τ
h
2
=
S
τ
h
2
(
3.186
×
1
0
-
8
e
-
0.6219
(
V
+
18.8
)
1
+
7.189
×
1
0
-
5
e
-
0.6683
(
V
+
34.07
)
+
3.556
×
1
0
-
3
)
h
1
∞
=
1
1
+
e
(
V
-
V
*
S
h
+
66.1
)
/
6.4
Effects of Ach on currents
i
K,Ach
=
g
K,Ach,max
j
k
[
Ach
]
n
K,Ach
K
0.5
,
K,Ach
n
K,Ach
+
[
Ach
]
n
K,Ach
(
[
K
+
]
o
10
+
[
K
+
]
o
)
·
(
V
-
E
K
1
+
e
(
V
-
E
K
-
140
)
F
/
2.5
R
T
)
i
Ca,
L
=
(
1
-
b
max
⁡
[
Ach
]
K
0.5
,
Ca
+
[
Ach
]
)
×
g
Ca
,
L
(
d
L
f
L
+
0.006
1
+
e
-
(
V
+
14.1
)
/
6.0
)
(
V
-
E
Ca,
L
)
i
f
=
y
[
g
f
,
Na
(
V
-
E
Na
)
+
g
f
,
K
(
V
-
E
K
)
]
d
y
d
t
=
y
∞
-
y
τ
y
y
∞
=
e
-
(
V
+
78.91
-
s
)
/
26.2
e
-
(
V
+
78.91
-
s
)
/
26.2
+
e
(
V
+
75.131
-
s
)
/
21.25
τ
y
=
1000.0
e
-
(
V
+
78.91
-
s
)
/
26.2
+
e
(
V
+
75.131
-
s
)
/
21.25
s
=
s
max
⁡
[
Ach
]
n
f
K
0.5
,
f
n
f
+
[
Ach
]
n
f


### 2.2. Effects of Iso on Ionic Currents

In order to simulate the chronotropic effect of Iso on the pacemaker activity, Zhang et al. have incorporated a set of equations of the known actions of Iso on related ion channels into two established action potential models of rabbit SAN recently [22]. According to the study of Zhang et al., we modified the models mentioned above as follows.(1)
Based on the voltage clamp experimental data on the kinetics of ionic channels and regional differences in the ionic current densities, we incorporated equations for concentration-dependent actions of Iso on the macroscopic conductance ofi
Ca
,
L, i
K
,
r, i
st, and i
K
,
s and to describe changes in the kinetics of i
f and i
K
,
r by the same amount as observed experimentally [23–25].
(2)
Because Iso was also found to increase the amplitude and minimal diastolic level of[
Ca
2
+
]
i in mammalian [26, 27] pacemaker cells and to alter SR Ca2+ uptake and release by stimulation of calmodulin kinase II [26, 28], in our simulations, we adopted the approach of Kharche et al. [29] to modify the Ca2+ handling equations to increase the amplitude and minimal diastolic level of [
Ca
2
+
]
i, as observed in experimental studies [30, 31] by increasing the maximal SR Ca2+ release (by 20%) and reducing the SR Ca2+ release (by 20%) fluxes.A general remodeling of these ionic channels is given in Table2. Using the models, the simulated effects of Iso are shown in Figure 2, in which the percentage of decrease of pacemaking cycle length (CL) is plotted against Iso concentration. The simulated data (solid line) from the central (Figure 2(a)) and the peripheral (Figure 2(b)) cell models were compared with the experimental data obtained from rabbit isolated SAN cells by Zaza et al. [23] (open circle) and Lei et al. [24, 25] (open diamond).Table 2
Actions of Iso on currents.
Increase ofL-type Ca2+ current
i
Ca,
L
=
[
(
1
+
f
Ca,max
[
Iso
]
K
0.5
,
Ca
′
+
[
Iso
]
)
×
g
Ca,
L
]
[
f
L
d
L
+
0.006
1
+
e
-
(
V
+
14.1
)
/
6
]
(
V
-
E
Ca,
L
)6
Shift of the activation curve of hyperpolarization-activated current
y
∞
=
e
-
(
V
+
78.91
-
S
f
)
/
26.2
e
-
(
V
+
78.91
-
S
f
)
/
26.2
+
e
(
V
+
75.131
-
S
f
)
/
21.25,  τ
y
=
1000.0
e
-
(
V
+
78.91
-
S
f
)
/
26.2
+
e
(
V
+
75.131
-
S
f
)
/
21.25
d
y
d
t
=
y
∞
-
y
τ
y,  S
f
=
S
f
,
max
⁡
[
Iso
]
n
f
′
K
0.5
,
f
′
n
f
′
+
[
Iso
]
n
f
′
Actions on rapid delayed rectifying K+ current
i
K
,
r
=
[
(
1
+
f
K
,
max
⁡
[
Iso
]
K
0.5
,
g
K
+
[
Iso
]
)
×
g
K
,
r
]
[
(
1
-
F
K
,
r
)
p
a
,
f
+
F
K
,
r
p
a
,
s
]
p
i
(
V
-
E
K
)
d
p
a
,
f
d
t
=
1
/
(
1
+
e
-
(
V
+
14.2
-
S
K
)
/
10.6
)
-
p
a
,
f
τ
p
a
,
f,  d
p
a
,
s
d
t
=
1
/
(
1
+
e
-
(
V
+
14.2
-
S
K
)
/
10.6
)
-
p
a
,
s
τ
p
a
,
s
S
K
=
S
K
,
max
⁡
[
Iso
]
K
0.5
,
Kacti
+
[
Iso
],  d
p
i
d
t
=
1
/
(
1
+
e
(
V
+
18.6
)
/
10.1
)
-
p
i
(
1
-
d
τ
K,
r
,
max
⁡
(
[
Iso
]
/
K
0.5
,
τ
K,
r
+
[
Iso
]
)
)
×
τ
p
i
Increase of slow delayed rectifying K+ current
i
K
,
s
=
[
(
1
+
f
K
,
max
⁡
[
Iso
]
K
0.5
,
g
K
+
[
Iso
]
)
×
g
K
,
s
]
x
s
2
(
V
-
E
K
,
s
)
Increase of sustained inward current
i
st
=
[
(
1
+
f
st,max
[
Iso
]
K
0.5
,
st
+
[
Iso
]
)
×
g
st
]
q
n
q
i
(
V
-
E
st
)
Actions on Ca2+ handling
j
up
=
(
1
-
β
Iso
)
×
P
up
(
[
C
a
2+
]
i
/
K
mf
)
n
up
-
(
[
C
a
2+
]
up
/
K
mr
)
n
up
1
+
(
[
C
a
2+
]
i
/
K
mf
)
n
up
-
(
[
C
a
2+
]
up
/
K
mr
)
n
up
j
rel
⁡
=
(
1
+
α
Iso
)
×
k
s
O
(
[
C
a
2+
]
rel
⁡
-
[
C
a
2+
]
sub
)Table 3
Model parameter values.
Glossary
Central SAN model
Peripheral SAN model
C
m
Cell membrane capacitance
20 pF
65 pF
g
Na
Maximumi
Na conductance
0μS/pF
1.85
×
1
0
-
8
μ
S/pF
g
K
,
Ach,max
Maximumi
K
,
Ach conductance
3.53
×
1
0
-
10 μS/pF
1.218
×
1
0
-
9
μ
S/pF
g
Ca,
L
Maximumi
Ca,
L conductance
2.9
×
1
0
-
4 μS/pF
1.0
×
1
0
-
3
μ
S/pF
g
f
,Na
Maximumi
f
,
Na conductance
0.27
×
1
0
-
4 μS/pF
1.05
×
1
0
-
4
μ
S/pF
g
f
,
K
Maximumi
f
,
K conductance
0.27
×
1
0
-
4 μS/pF
1.05
×
1
0
-
4
μ
S/pF
g
K
,
r
Maximumi
K
,
r conductance
3.99
×
1
0
-
5 μS/pF
2.46
×
1
0
-
4
μ
S/pF
g
K
,
s
Maximumi
K
,
s conductance
2.59
×
1
0
-
5 μS/pF
1.6
×
1
0
-
4
μ
S/pF
g
st
Maximumi
st conductance
0.75
×
10
-
5 μS/pF
0μS/pF
E
Ca,
L
Apparent reversal potential fori
Ca,
L
46.4 mV
46.4 mV
E
st
Apparent reversal potential fori
st
37.4 mV
37.4 mV
[
Na
+
]
o
Extracellular Na+ concentration
140 mM
140 mM
[
K
+
]
o
Extracellular K+ concentration
5.4 mM
5.4 mM
K
0.5
,
K
,
Ach
[Ach] that produces a half-maximal activation ofg
K
,
Ach,max
0.28μM
0.28μM
K
0.5
,
Ca
[Ach] that produces a half-maximal block ofi
Ca,
L
0.09μM
0.09μM
K
0.5
,
Ca
′
[Iso] that produces a half-maximal increase ofi
Ca,
L
7 nM
7 nM
K
0.5
,
f
[Ach] that produces a half-maximal shift ofi
f activation curve
1.26
×
10
-
2 μM
1.26
×
10
-
2 μM
K
0.5
,
f
′
[Iso] that produces a half-maximal shift ofi
f activation curve
13.5 nM
13.5 nM
K
0.5
,
g
K
[Iso] that produces a half-maximal increase ofg
K
,
r
19 nM
19 nM
K
0.5
,
Kacti
[Iso] that produces a half-maximal shift ofi
K
,
r activation curve
7.5 nM
7.5 nM
K
0.5
,
τ
K,
r
[Iso] that produces a half-maximal decrease ofτ
K,
r
24 nM
24 nM
K
0.5
,
st
[Iso] that produces a half-maximal Increase ofi
st
33 nM
33 nM
P
up
Rate constant of Ca2+ uptake by j
up of the network SR
0.01 mM/ms
0.02 mM/ms
n
K
,
Ach
i
K,Ach affected by [Ach] and Hill coefficient
1.5
1.5
n
f
i
f affected by [Ach] and Hill coefficient
0.69
0.69
n
up
SR Ca2+ uptake and Hill coefficient
2
2
n
f
′
i
f affected by [Iso] and Hill coefficient
0.392
0.392
S
Na
Percentage change ofi
Na conductance induced by DelF1617 mutation
0.38
0.38
S
τ
h
1
Percentage change of the fast-inactivation time constants ofi
Na induced by DelF1617 mutation
2.18
2.18
S
τ
h
2
Percentage change of the slow-inactivation time constants ofi
Na induced by DelF1617 mutation
2.75
2.75
S
h
Parameter of the shift in the inactivation curve ofi
Na induced by DelF1617 mutation
0.128
0.128
b
max
⁡
Maximum fraction ofi
Ca,
L block caused by Ach
0.56
0.56
s
max
⁡
Maximum shift ofi
f activation curve caused by Ach
−7.2 mV
−7.2 mV
f
Ca
,
max
⁡
Maximum percentage increase ofi
Ca,
L caused by Iso
0.54
0.54
f
K
,
max
⁡
Maximum percentage increase ofi
K
,
r caused by Iso
1.87
1.87
f
st.max
Maximum percentage increase ofi
st caused by Iso
1.0
1.0
S
f
,max
Maximum activation curve shift ofi
f caused by Iso
9.62 mV
9.62 mV
S
K
,
max
⁡
Maximum activation curve shift ofi
K
,
r caused by Iso
−15 mV
−15 mV
d
τ
K,
r,
max
⁡
Maximum decrease ofτ
K,
r, caused by Iso
3.0
3.0
α
Iso
Percentage increase of the SR Ca2+ release fluxes caused by Iso
0.2
0.2
β
Iso
Percentage decrease of the SR Ca2+ uptake caused by Iso
0.2
0.2Concentration-dependent percentage decrease of pacemaking CL by Iso. (a) Central cell model and relevant experimental data. (b) Peripheral cell model and relevant experimental data. Open circles: experimental data of Zaza et al. from rabbit SAN cells. Open diamonds: data of Lei et al. from rabbit SAN cells.
(a)
(b)


## 2.1. Model Description

In order to address the above issues, mathematical model of the rabbit central and peripheral SAN cells proposed by Kurata et al. [16] was used in this paper. The action potentials generated by using the models of a central cell and a peripheral cell were shown in Figure 1(c). Besides, a 2D anatomic model of the intact SAN-atrium tissue developed by Butters et al. [17] was also used to study the effects of Iso on SAN pacemaking and driving. This 2D intact SAN-atrium tissue model was based on histologically reconstructed geometry of a single slice of the rabbit right atrium, which was cut through the atrial muscle of the crista terminalis (CT) and the intercaval region with central and peripheral SAN areas. And an experimentally observed [18, 19] nonconductive region (block zone) next to the SAN toward the atrial septum was also incorporated in this model (see Figure 1(a)). Figure 1(b) shows the initiation and conduction of pacemaker activity computed from the 2D model. Snapshots of membrane potentials across the 2D lattice were shown at 50 ms (Figure 1(b)(i)), 70 ms (Figure 1(b)(ii)), and 90 ms (Figure 1(b)(iii)) after the initial conduction. The pacemaking action was first initiated in the center (Figure 1(b)(i)) and then propagated toward the surrounding tissues (Figure 1(b)(iii)).Figure 1
Model of the SAN and surrounding atrial tissue. (a) Color-coded distribution of cell types throughout the 2D tissue slice. (b) Snapshots of initiation and conduction in 2D anatomical model at various times after initial configuration; (i) 50 ms; (ii) 70 ms; (iii) 90 ms. (c) Action potentials in one cell of the central and peripheral SAN. (d) Action potential profiles during conduction through the slice. As, atrial septum; SA-C, SAN center; AM, atrial muscle (CT).The mechanism underlying the genesis of SND (clinically known as SSS) is believed to be associated with changes in the intrinsic properties of the SAN with age [20]. Intercellular electrical coupling strength (D) is considered to reflect the changes in the intrinsic properties of the SAN. To incorporate the effects of SCN5A mutation and acetylcholine (Ach) into the models, (i) parameters of the Na+ current were changed according to the study of Butters et al. [17] and (ii) an Ach-activated K+ current (i
K
,
Ach), the inhibited L-type Ca2+ current (i
Ca
,
L), and the shift in hyperpolarization-activated current (i
f) activation curve were introduced.The model equations used for the cells, 2D tissue, modified Na+ channel current, and the channel currents relative to Ach are presented in Table 1. A detailed description of parameters of these models are listed in Table 3 and can be found in [16, 19–21].Table 1
Model equations of the SAN single cell, 2D tissue, and  some channel currents.
Single cell model of SAN
d
V
/
d
t
=
-
(
i
Na
+
i
Ca,
L
+
i
Ca,
T
+
i
K,
r
+
i
K,
s
+
i
st
+
i
to
+
i
sus
+
i
b
,Na
+
i
b
,K
+
i
NaK
+
i
NaCa
)
/
C
m
2D intact SAN-atrium tissue model
∂
V
∂
t
=
∇
·
(
D
∇
V
)
-
i
tot
C
m
Changes of sodium current induced by mutation
i
Na
=
(
S
Na
×
g
Na
)
m
3
h
[
Na
+
]
o
F
2
R
T
e
(
V
-
E
Na
)
F
/
R
T
-
1
e
V
F
/
R
T
-
1
V
τ
h
1
=
S
τ
h
1
(
3.717
×
1
0
-
6
e
-
0.2815
(
V
+
17.11
)
1
+
3.732
×
1
0
-
3
e
-
0.3426
(
V
+
37.76
)
+
5.977
×
1
0
-
4
)
τ
h
2
=
S
τ
h
2
(
3.186
×
1
0
-
8
e
-
0.6219
(
V
+
18.8
)
1
+
7.189
×
1
0
-
5
e
-
0.6683
(
V
+
34.07
)
+
3.556
×
1
0
-
3
)
h
1
∞
=
1
1
+
e
(
V
-
V
*
S
h
+
66.1
)
/
6.4
Effects of Ach on currents
i
K,Ach
=
g
K,Ach,max
j
k
[
Ach
]
n
K,Ach
K
0.5
,
K,Ach
n
K,Ach
+
[
Ach
]
n
K,Ach
(
[
K
+
]
o
10
+
[
K
+
]
o
)
·
(
V
-
E
K
1
+
e
(
V
-
E
K
-
140
)
F
/
2.5
R
T
)
i
Ca,
L
=
(
1
-
b
max
⁡
[
Ach
]
K
0.5
,
Ca
+
[
Ach
]
)
×
g
Ca
,
L
(
d
L
f
L
+
0.006
1
+
e
-
(
V
+
14.1
)
/
6.0
)
(
V
-
E
Ca,
L
)
i
f
=
y
[
g
f
,
Na
(
V
-
E
Na
)
+
g
f
,
K
(
V
-
E
K
)
]
d
y
d
t
=
y
∞
-
y
τ
y
y
∞
=
e
-
(
V
+
78.91
-
s
)
/
26.2
e
-
(
V
+
78.91
-
s
)
/
26.2
+
e
(
V
+
75.131
-
s
)
/
21.25
τ
y
=
1000.0
e
-
(
V
+
78.91
-
s
)
/
26.2
+
e
(
V
+
75.131
-
s
)
/
21.25
s
=
s
max
⁡
[
Ach
]
n
f
K
0.5
,
f
n
f
+
[
Ach
]
n
f


## 2.2. Effects of Iso on Ionic Currents

In order to simulate the chronotropic effect of Iso on the pacemaker activity, Zhang et al. have incorporated a set of equations of the known actions of Iso on related ion channels into two established action potential models of rabbit SAN recently [22]. According to the study of Zhang et al., we modified the models mentioned above as follows.(1)
Based on the voltage clamp experimental data on the kinetics of ionic channels and regional differences in the ionic current densities, we incorporated equations for concentration-dependent actions of Iso on the macroscopic conductance ofi
Ca
,
L, i
K
,
r, i
st, and i
K
,
s and to describe changes in the kinetics of i
f and i
K
,
r by the same amount as observed experimentally [23–25].
(2)
Because Iso was also found to increase the amplitude and minimal diastolic level of[
Ca
2
+
]
i in mammalian [26, 27] pacemaker cells and to alter SR Ca2+ uptake and release by stimulation of calmodulin kinase II [26, 28], in our simulations, we adopted the approach of Kharche et al. [29] to modify the Ca2+ handling equations to increase the amplitude and minimal diastolic level of [
Ca
2
+
]
i, as observed in experimental studies [30, 31] by increasing the maximal SR Ca2+ release (by 20%) and reducing the SR Ca2+ release (by 20%) fluxes.A general remodeling of these ionic channels is given in Table2. Using the models, the simulated effects of Iso are shown in Figure 2, in which the percentage of decrease of pacemaking cycle length (CL) is plotted against Iso concentration. The simulated data (solid line) from the central (Figure 2(a)) and the peripheral (Figure 2(b)) cell models were compared with the experimental data obtained from rabbit isolated SAN cells by Zaza et al. [23] (open circle) and Lei et al. [24, 25] (open diamond).Table 2
Actions of Iso on currents.
Increase ofL-type Ca2+ current
i
Ca,
L
=
[
(
1
+
f
Ca,max
[
Iso
]
K
0.5
,
Ca
′
+
[
Iso
]
)
×
g
Ca,
L
]
[
f
L
d
L
+
0.006
1
+
e
-
(
V
+
14.1
)
/
6
]
(
V
-
E
Ca,
L
)6
Shift of the activation curve of hyperpolarization-activated current
y
∞
=
e
-
(
V
+
78.91
-
S
f
)
/
26.2
e
-
(
V
+
78.91
-
S
f
)
/
26.2
+
e
(
V
+
75.131
-
S
f
)
/
21.25,  τ
y
=
1000.0
e
-
(
V
+
78.91
-
S
f
)
/
26.2
+
e
(
V
+
75.131
-
S
f
)
/
21.25
d
y
d
t
=
y
∞
-
y
τ
y,  S
f
=
S
f
,
max
⁡
[
Iso
]
n
f
′
K
0.5
,
f
′
n
f
′
+
[
Iso
]
n
f
′
Actions on rapid delayed rectifying K+ current
i
K
,
r
=
[
(
1
+
f
K
,
max
⁡
[
Iso
]
K
0.5
,
g
K
+
[
Iso
]
)
×
g
K
,
r
]
[
(
1
-
F
K
,
r
)
p
a
,
f
+
F
K
,
r
p
a
,
s
]
p
i
(
V
-
E
K
)
d
p
a
,
f
d
t
=
1
/
(
1
+
e
-
(
V
+
14.2
-
S
K
)
/
10.6
)
-
p
a
,
f
τ
p
a
,
f,  d
p
a
,
s
d
t
=
1
/
(
1
+
e
-
(
V
+
14.2
-
S
K
)
/
10.6
)
-
p
a
,
s
τ
p
a
,
s
S
K
=
S
K
,
max
⁡
[
Iso
]
K
0.5
,
Kacti
+
[
Iso
],  d
p
i
d
t
=
1
/
(
1
+
e
(
V
+
18.6
)
/
10.1
)
-
p
i
(
1
-
d
τ
K,
r
,
max
⁡
(
[
Iso
]
/
K
0.5
,
τ
K,
r
+
[
Iso
]
)
)
×
τ
p
i
Increase of slow delayed rectifying K+ current
i
K
,
s
=
[
(
1
+
f
K
,
max
⁡
[
Iso
]
K
0.5
,
g
K
+
[
Iso
]
)
×
g
K
,
s
]
x
s
2
(
V
-
E
K
,
s
)
Increase of sustained inward current
i
st
=
[
(
1
+
f
st,max
[
Iso
]
K
0.5
,
st
+
[
Iso
]
)
×
g
st
]
q
n
q
i
(
V
-
E
st
)
Actions on Ca2+ handling
j
up
=
(
1
-
β
Iso
)
×
P
up
(
[
C
a
2+
]
i
/
K
mf
)
n
up
-
(
[
C
a
2+
]
up
/
K
mr
)
n
up
1
+
(
[
C
a
2+
]
i
/
K
mf
)
n
up
-
(
[
C
a
2+
]
up
/
K
mr
)
n
up
j
rel
⁡
=
(
1
+
α
Iso
)
×
k
s
O
(
[
C
a
2+
]
rel
⁡
-
[
C
a
2+
]
sub
)Table 3
Model parameter values.
Glossary
Central SAN model
Peripheral SAN model
C
m
Cell membrane capacitance
20 pF
65 pF
g
Na
Maximumi
Na conductance
0μS/pF
1.85
×
1
0
-
8
μ
S/pF
g
K
,
Ach,max
Maximumi
K
,
Ach conductance
3.53
×
1
0
-
10 μS/pF
1.218
×
1
0
-
9
μ
S/pF
g
Ca,
L
Maximumi
Ca,
L conductance
2.9
×
1
0
-
4 μS/pF
1.0
×
1
0
-
3
μ
S/pF
g
f
,Na
Maximumi
f
,
Na conductance
0.27
×
1
0
-
4 μS/pF
1.05
×
1
0
-
4
μ
S/pF
g
f
,
K
Maximumi
f
,
K conductance
0.27
×
1
0
-
4 μS/pF
1.05
×
1
0
-
4
μ
S/pF
g
K
,
r
Maximumi
K
,
r conductance
3.99
×
1
0
-
5 μS/pF
2.46
×
1
0
-
4
μ
S/pF
g
K
,
s
Maximumi
K
,
s conductance
2.59
×
1
0
-
5 μS/pF
1.6
×
1
0
-
4
μ
S/pF
g
st
Maximumi
st conductance
0.75
×
10
-
5 μS/pF
0μS/pF
E
Ca,
L
Apparent reversal potential fori
Ca,
L
46.4 mV
46.4 mV
E
st
Apparent reversal potential fori
st
37.4 mV
37.4 mV
[
Na
+
]
o
Extracellular Na+ concentration
140 mM
140 mM
[
K
+
]
o
Extracellular K+ concentration
5.4 mM
5.4 mM
K
0.5
,
K
,
Ach
[Ach] that produces a half-maximal activation ofg
K
,
Ach,max
0.28μM
0.28μM
K
0.5
,
Ca
[Ach] that produces a half-maximal block ofi
Ca,
L
0.09μM
0.09μM
K
0.5
,
Ca
′
[Iso] that produces a half-maximal increase ofi
Ca,
L
7 nM
7 nM
K
0.5
,
f
[Ach] that produces a half-maximal shift ofi
f activation curve
1.26
×
10
-
2 μM
1.26
×
10
-
2 μM
K
0.5
,
f
′
[Iso] that produces a half-maximal shift ofi
f activation curve
13.5 nM
13.5 nM
K
0.5
,
g
K
[Iso] that produces a half-maximal increase ofg
K
,
r
19 nM
19 nM
K
0.5
,
Kacti
[Iso] that produces a half-maximal shift ofi
K
,
r activation curve
7.5 nM
7.5 nM
K
0.5
,
τ
K,
r
[Iso] that produces a half-maximal decrease ofτ
K,
r
24 nM
24 nM
K
0.5
,
st
[Iso] that produces a half-maximal Increase ofi
st
33 nM
33 nM
P
up
Rate constant of Ca2+ uptake by j
up of the network SR
0.01 mM/ms
0.02 mM/ms
n
K
,
Ach
i
K,Ach affected by [Ach] and Hill coefficient
1.5
1.5
n
f
i
f affected by [Ach] and Hill coefficient
0.69
0.69
n
up
SR Ca2+ uptake and Hill coefficient
2
2
n
f
′
i
f affected by [Iso] and Hill coefficient
0.392
0.392
S
Na
Percentage change ofi
Na conductance induced by DelF1617 mutation
0.38
0.38
S
τ
h
1
Percentage change of the fast-inactivation time constants ofi
Na induced by DelF1617 mutation
2.18
2.18
S
τ
h
2
Percentage change of the slow-inactivation time constants ofi
Na induced by DelF1617 mutation
2.75
2.75
S
h
Parameter of the shift in the inactivation curve ofi
Na induced by DelF1617 mutation
0.128
0.128
b
max
⁡
Maximum fraction ofi
Ca,
L block caused by Ach
0.56
0.56
s
max
⁡
Maximum shift ofi
f activation curve caused by Ach
−7.2 mV
−7.2 mV
f
Ca
,
max
⁡
Maximum percentage increase ofi
Ca,
L caused by Iso
0.54
0.54
f
K
,
max
⁡
Maximum percentage increase ofi
K
,
r caused by Iso
1.87
1.87
f
st.max
Maximum percentage increase ofi
st caused by Iso
1.0
1.0
S
f
,max
Maximum activation curve shift ofi
f caused by Iso
9.62 mV
9.62 mV
S
K
,
max
⁡
Maximum activation curve shift ofi
K
,
r caused by Iso
−15 mV
−15 mV
d
τ
K,
r,
max
⁡
Maximum decrease ofτ
K,
r, caused by Iso
3.0
3.0
α
Iso
Percentage increase of the SR Ca2+ release fluxes caused by Iso
0.2
0.2
β
Iso
Percentage decrease of the SR Ca2+ uptake caused by Iso
0.2
0.2Concentration-dependent percentage decrease of pacemaking CL by Iso. (a) Central cell model and relevant experimental data. (b) Peripheral cell model and relevant experimental data. Open circles: experimental data of Zaza et al. from rabbit SAN cells. Open diamonds: data of Lei et al. from rabbit SAN cells.
(a)
(b)


## 3. Results and Discussion

### 3.1. Iso Effects on SAN Cells

We discussed the effects of Iso on the dynamics of action potential in cells from the modified model. The simulated results of the action potentials under control conditions (i.e., in the absence of Iso) and in the presence of Iso (i.e., [Iso] = 0.1μM/L) are depicted in Figure 3. One can find that, with both the central and peripheral cell models, Iso shortened the action potential duration, increased the amplitude of action potentials, and caused an increase in the rate of spontaneous action potentials (Figures 3(a)(i) and 3(b)(i), dashed line). These results are consistent with experimental observations on rabbit SAN cells [23, 24].Effect of Iso on the SAN pacemaking rate. Simulated voltage from (a) central cell models and ((b) and (c)) peripheral cell models. Action potentials are shown under control conditions and in the presence 0.1μM/L Iso, 0.01 μM/L Ach, and the DelF1617 mutant channel.
(a)
(b)
(c)And then, we introduced the effect of Ach into the modified models. Application of Ach slowed down pacemaking rates in both central and peripheral SAN cells, with larger effect on the central cells [32]. Simulating effects of 0.01 μM/L Ach resulted in the CL increasing from 406 ms to 660 ms for the central SAN cell (Figure 3(a)(ii), solid line) and from 280 ms to 460 ms for the peripheral ones (Figure 3(b)(ii), solid line). When Iso was added (i.e, [Iso] = 0.1 μM/L), the negative chronotropic effect of Ach was diminished for both types of cell models (Figures 3(a)(ii) and 3(b)(ii), dashed line).It has been shown that the mutations slowed down pacemaking rate in peripheral, but not in central, SAN cells that control the heart rhythm [17]. Thus, we simulated the effect of Iso on pacemaking action potentials with the DelF1617 mutant channel for the peripheral SAN cells, and the results are shown in Figures 3(c)(i) and 3(c)(ii). One notices that, when only the effects of the DelF1617 mutation were introduced into the peripheral cells, the CL increased from 280 ms to 344 ms (Figure 3(c)(i), solid line). When Iso was added and its concentration was adjusted to be 0.1 μM/L, the CL decreased by 14.2% (Figure 3(c)(i), dashed line). Furthermore, the negative chronotropic effect of Ach (concentration is set to be 0.01 μM/L) was amplified by the mutation; that is, the cell with the wild type channel exhibited pacemaking activity (Figure 3(b)(ii), solid line), but the cell with the mutant channel became quiescent (Figure 3(c)(ii), solid line), and in this case, addition of Iso (i.e., [Iso] = 0.1 μM/L) could recover the pacemaking activity of the cell with the mutant channel (Figure 3(c)(ii), dashed line). The above adjustive behavior of Iso is similar to the clinical behavior observed in patients with a longer CL [33].


### 3.2. Iso Effects on Two-Dimensional Tissues

We investigated further the functional influences of Iso on action potential conduction across the intact SAN-atrium under control and other conditions associated with the symptoms of SSS (i.e., DelF1617 mutation; [Ach] = 0.03μM/L; D
=
0.1 mm2/ms). Results are shown in Figure 4, which presents the spatial and temporal profiles of action potentials recorded from representative cells across the SAN-atrium in the model (Figure 1(a)). And the conduction in the normal tissue (with the wild type channel, D
=
1.0 mm2/ms, no Iso and no Ach) is presented in Figure 1(d).Effects of Iso on action potential conduction. Action potential profiles in the 2D tissue with the DelF1617 mutant channel under different conditions are shown. (a) Mutant channel alone; (b) [Iso] = 0.5μM/L; (c) [Ach] = 0.03 μM/L; (d) [Ach] = 0.03 μM/L and D
=
0.1; (e) [Ach] = 0.03 μM/L, D
=
0.1, and [Iso] = 0.5 μM/L; (f) [Ach] = 0.03 μM/L, D
=
0.1, and [Iso] = 2.0 μM/L.
(a)
(b)
(c)
(d)
(e)
(f)Figure3(a) shows the action potential profile with the DelF1617 mutation. The CL was increased by 25.8% compared with the normal tissue, and the ability of the SAN to drive atrium was not impaired. However, with Ach (i.e., [Ach] = 0.03 μM/L), the mutation significantly increased the CL (Figure 4(c)) and produced a conduction block in the direction toward the atrial septum; the action potential conduction in the direction toward the CT was sustained. As we know, electrotonic interaction between the SAN and the surrounding tissues played an important role in determining the SAN pacing and driving. When the coupling strength D was fixed to be small (i.e., D
=
0.1 mm2/ms), the SAN was able to generate spontaneous activity but failed to drive the surrounding atrium; conduction exit block occurred in both directions (Figure 4(d)).When Iso addition was simulated at tissue level, three noteworthy results can be addressed. First, pacemaking rate with the DelF1617 mutation was significantly increased. The measured CL was decreased from 472.5 ms (Figure4(a)) to 414.2 ms (Figure 4(b)) in the presence of 0.5 μM/L Iso. Second, the spatial cessation of pacemaker activity in the direction toward the atrial septum cannot be recovered by addition of Iso (Figure 4(c)). Third, a conduction block in the direction toward the CT can be avoided. Through the application of 0.5 μM/L Iso, the action potential was initiated in the center and propagated toward the CT, but not continually (Figure 4(e)); however, with the Iso concentration increasing to 2.0 μM/L, the spatial cessation of pacemaker activity can be recovered completely (Figure 4(f)).To characterize the effects of Iso, a series of simulations were performed with [Iso] changing systematically from 0.0001μM/L to 10 μM/L. The results are represented in Figure 5(a), in which the computed CL was plotted against [Iso] under control conditions (solid circles), DelF1617 mutation (solid triangles), DelF1617 mutation together with [Ach] = 0.03 μM/L (solid inverted triangles), and D
=
0.1 mm2/ms (open circles). Under both control and DelF1617 mutation conditions, the SAN could pace and drive the atrium successfully. However, in the presence of 0.03 μM/L Ach, the measured CL decreased with [Iso] in a similar way as under both control and DelF1617 mutation conditions; though the SAN paces, it does not drive the atrial septum (Figure 5(a), solid inverted triangles, dashed line).(a) Measured pacemaking CL against [Iso] computed from tissue model. Solid circles: controls; solid triangle: DelF1617 mutation; solid inverted triangles: DelF1617 and [Ach] = 0.03μM/L; open circles: DelF1617, [Ach] = 0.03 μM/L, and D
=
0.1; solid line: without conduction block; dashed line: conduction block occurs in AS. (b) Distribution of CL in controls and patients (i) before and (ii) after Iso injection. Open circles: clinical data of Vallin and Edhag from controls; solid circles: data of Vallin and Edhag from patients.
(a)
(b)When coupling strength was considered, another important phenomenon was that, in the range of 10−4 μM/L ≤ [Iso] ≤ 10−2 μM/L, the SAN failed to drive the surrounding tissues in both sides (Figure 5(a), open circles without line). When [Iso] is over 10−1 μM/L, the SAN could drive the atrium continually (Figure 5(a), open circles, dashed line); however, conduction block in the atrial septum could not be avoided. Such transition happened because Iso has shown its positive effects on the pacemaker activity in the SAN. With [Iso] ≥ 10 μM/L, the pacemaker activity of the SAN was abolished completely, which suggested that high concentrations of Iso are considered toxic. These results are consistent with clinical observations by Vallin and Edhag [33] who injected Iso in 30 patients with symptoms of sinus node disease and 18 control subjects (Figure 5(b)).The above simulations suggested that Iso plays an important role in determining the initiation and conduction of the pacemaker activity in the SAN, especially under the conditions of mutation, high [Ach], and weak coupling strength.


### 3.3. Role ofi
Na and i
Ca
,
T

Due to the functional role of several major ionic currents, it is necessary to study the corresponding deterministic kinetics of these channel currents. Figure6 presents the time series of the ionic currents, in which some useful clues can be found to infer the mechanism responsible for the positive effect of Iso on the pacemaker activity.Time series ofi
Na and i
Ca
,
T recording from a SAN peripheral cell under different conditions. (a) i
Na; (b) i
Ca
,
T. Tissues with normal conditions (i), DelF1617 mutant channel (ii), and Iso with mutant channel (iii) are shown.
(a)
(b)Figures6(a)(i), (6(b)(i)), 6(a)(ii), and (6(b)(ii)) show the time series of i
Na (i
Ca
,
T) under normal condition and with the DelF1617 mutation. Compared with Figures 6(a)(i) and (6(b)(i)), i
Na (i
Ca
,
T) was reduced or even zero in Figures 6(a)(ii) and (6(b)(ii)). Besides, other major underlying currents (such as i
Ca
,
L and i
K
,
r) were not greatly affected (data are not shown). Thus, we can conclude that the slowing of the pacemaking was primarily attributable to the mutation-induced decrease in i
Na and i
Ca
,
T (Figures 6(a)(ii) and 6(b)(ii)).However, in this case, through the application of Iso (i.e., [Iso] = 0.1μM/L), i
Na was significantly increased as well as i
Ca
,
T (Figures 6(a)(iii) and 6(b)(iii)), and both of the currents can be even increased to a normal level. This means that proper concentration of Iso can increase i
Na and i
Ca
,
T which were reduced by the mutation. As a result of changes in the two currents, such abnormal pacemaking behavior could be eliminated.


## 3.1. Iso Effects on SAN Cells

We discussed the effects of Iso on the dynamics of action potential in cells from the modified model. The simulated results of the action potentials under control conditions (i.e., in the absence of Iso) and in the presence of Iso (i.e., [Iso] = 0.1μM/L) are depicted in Figure 3. One can find that, with both the central and peripheral cell models, Iso shortened the action potential duration, increased the amplitude of action potentials, and caused an increase in the rate of spontaneous action potentials (Figures 3(a)(i) and 3(b)(i), dashed line). These results are consistent with experimental observations on rabbit SAN cells [23, 24].Effect of Iso on the SAN pacemaking rate. Simulated voltage from (a) central cell models and ((b) and (c)) peripheral cell models. Action potentials are shown under control conditions and in the presence 0.1μM/L Iso, 0.01 μM/L Ach, and the DelF1617 mutant channel.
(a)
(b)
(c)And then, we introduced the effect of Ach into the modified models. Application of Ach slowed down pacemaking rates in both central and peripheral SAN cells, with larger effect on the central cells [32]. Simulating effects of 0.01 μM/L Ach resulted in the CL increasing from 406 ms to 660 ms for the central SAN cell (Figure 3(a)(ii), solid line) and from 280 ms to 460 ms for the peripheral ones (Figure 3(b)(ii), solid line). When Iso was added (i.e, [Iso] = 0.1 μM/L), the negative chronotropic effect of Ach was diminished for both types of cell models (Figures 3(a)(ii) and 3(b)(ii), dashed line).It has been shown that the mutations slowed down pacemaking rate in peripheral, but not in central, SAN cells that control the heart rhythm [17]. Thus, we simulated the effect of Iso on pacemaking action potentials with the DelF1617 mutant channel for the peripheral SAN cells, and the results are shown in Figures 3(c)(i) and 3(c)(ii). One notices that, when only the effects of the DelF1617 mutation were introduced into the peripheral cells, the CL increased from 280 ms to 344 ms (Figure 3(c)(i), solid line). When Iso was added and its concentration was adjusted to be 0.1 μM/L, the CL decreased by 14.2% (Figure 3(c)(i), dashed line). Furthermore, the negative chronotropic effect of Ach (concentration is set to be 0.01 μM/L) was amplified by the mutation; that is, the cell with the wild type channel exhibited pacemaking activity (Figure 3(b)(ii), solid line), but the cell with the mutant channel became quiescent (Figure 3(c)(ii), solid line), and in this case, addition of Iso (i.e., [Iso] = 0.1 μM/L) could recover the pacemaking activity of the cell with the mutant channel (Figure 3(c)(ii), dashed line). The above adjustive behavior of Iso is similar to the clinical behavior observed in patients with a longer CL [33].


## 3.2. Iso Effects on Two-Dimensional Tissues

We investigated further the functional influences of Iso on action potential conduction across the intact SAN-atrium under control and other conditions associated with the symptoms of SSS (i.e., DelF1617 mutation; [Ach] = 0.03μM/L; D
=
0.1 mm2/ms). Results are shown in Figure 4, which presents the spatial and temporal profiles of action potentials recorded from representative cells across the SAN-atrium in the model (Figure 1(a)). And the conduction in the normal tissue (with the wild type channel, D
=
1.0 mm2/ms, no Iso and no Ach) is presented in Figure 1(d).Effects of Iso on action potential conduction. Action potential profiles in the 2D tissue with the DelF1617 mutant channel under different conditions are shown. (a) Mutant channel alone; (b) [Iso] = 0.5μM/L; (c) [Ach] = 0.03 μM/L; (d) [Ach] = 0.03 μM/L and D
=
0.1; (e) [Ach] = 0.03 μM/L, D
=
0.1, and [Iso] = 0.5 μM/L; (f) [Ach] = 0.03 μM/L, D
=
0.1, and [Iso] = 2.0 μM/L.
(a)
(b)
(c)
(d)
(e)
(f)Figure3(a) shows the action potential profile with the DelF1617 mutation. The CL was increased by 25.8% compared with the normal tissue, and the ability of the SAN to drive atrium was not impaired. However, with Ach (i.e., [Ach] = 0.03 μM/L), the mutation significantly increased the CL (Figure 4(c)) and produced a conduction block in the direction toward the atrial septum; the action potential conduction in the direction toward the CT was sustained. As we know, electrotonic interaction between the SAN and the surrounding tissues played an important role in determining the SAN pacing and driving. When the coupling strength D was fixed to be small (i.e., D
=
0.1 mm2/ms), the SAN was able to generate spontaneous activity but failed to drive the surrounding atrium; conduction exit block occurred in both directions (Figure 4(d)).When Iso addition was simulated at tissue level, three noteworthy results can be addressed. First, pacemaking rate with the DelF1617 mutation was significantly increased. The measured CL was decreased from 472.5 ms (Figure4(a)) to 414.2 ms (Figure 4(b)) in the presence of 0.5 μM/L Iso. Second, the spatial cessation of pacemaker activity in the direction toward the atrial septum cannot be recovered by addition of Iso (Figure 4(c)). Third, a conduction block in the direction toward the CT can be avoided. Through the application of 0.5 μM/L Iso, the action potential was initiated in the center and propagated toward the CT, but not continually (Figure 4(e)); however, with the Iso concentration increasing to 2.0 μM/L, the spatial cessation of pacemaker activity can be recovered completely (Figure 4(f)).To characterize the effects of Iso, a series of simulations were performed with [Iso] changing systematically from 0.0001μM/L to 10 μM/L. The results are represented in Figure 5(a), in which the computed CL was plotted against [Iso] under control conditions (solid circles), DelF1617 mutation (solid triangles), DelF1617 mutation together with [Ach] = 0.03 μM/L (solid inverted triangles), and D
=
0.1 mm2/ms (open circles). Under both control and DelF1617 mutation conditions, the SAN could pace and drive the atrium successfully. However, in the presence of 0.03 μM/L Ach, the measured CL decreased with [Iso] in a similar way as under both control and DelF1617 mutation conditions; though the SAN paces, it does not drive the atrial septum (Figure 5(a), solid inverted triangles, dashed line).(a) Measured pacemaking CL against [Iso] computed from tissue model. Solid circles: controls; solid triangle: DelF1617 mutation; solid inverted triangles: DelF1617 and [Ach] = 0.03μM/L; open circles: DelF1617, [Ach] = 0.03 μM/L, and D
=
0.1; solid line: without conduction block; dashed line: conduction block occurs in AS. (b) Distribution of CL in controls and patients (i) before and (ii) after Iso injection. Open circles: clinical data of Vallin and Edhag from controls; solid circles: data of Vallin and Edhag from patients.
(a)
(b)When coupling strength was considered, another important phenomenon was that, in the range of 10−4 μM/L ≤ [Iso] ≤ 10−2 μM/L, the SAN failed to drive the surrounding tissues in both sides (Figure 5(a), open circles without line). When [Iso] is over 10−1 μM/L, the SAN could drive the atrium continually (Figure 5(a), open circles, dashed line); however, conduction block in the atrial septum could not be avoided. Such transition happened because Iso has shown its positive effects on the pacemaker activity in the SAN. With [Iso] ≥ 10 μM/L, the pacemaker activity of the SAN was abolished completely, which suggested that high concentrations of Iso are considered toxic. These results are consistent with clinical observations by Vallin and Edhag [33] who injected Iso in 30 patients with symptoms of sinus node disease and 18 control subjects (Figure 5(b)).The above simulations suggested that Iso plays an important role in determining the initiation and conduction of the pacemaker activity in the SAN, especially under the conditions of mutation, high [Ach], and weak coupling strength.


## 3.3. Role ofi
Na and i
Ca
,
T

Due to the functional role of several major ionic currents, it is necessary to study the corresponding deterministic kinetics of these channel currents. Figure6 presents the time series of the ionic currents, in which some useful clues can be found to infer the mechanism responsible for the positive effect of Iso on the pacemaker activity.Time series ofi
Na and i
Ca
,
T recording from a SAN peripheral cell under different conditions. (a) i
Na; (b) i
Ca
,
T. Tissues with normal conditions (i), DelF1617 mutant channel (ii), and Iso with mutant channel (iii) are shown.
(a)
(b)Figures6(a)(i), (6(b)(i)), 6(a)(ii), and (6(b)(ii)) show the time series of i
Na (i
Ca
,
T) under normal condition and with the DelF1617 mutation. Compared with Figures 6(a)(i) and (6(b)(i)), i
Na (i
Ca
,
T) was reduced or even zero in Figures 6(a)(ii) and (6(b)(ii)). Besides, other major underlying currents (such as i
Ca
,
L and i
K
,
r) were not greatly affected (data are not shown). Thus, we can conclude that the slowing of the pacemaking was primarily attributable to the mutation-induced decrease in i
Na and i
Ca
,
T (Figures 6(a)(ii) and 6(b)(ii)).However, in this case, through the application of Iso (i.e., [Iso] = 0.1μM/L), i
Na was significantly increased as well as i
Ca
,
T (Figures 6(a)(iii) and 6(b)(iii)), and both of the currents can be even increased to a normal level. This means that proper concentration of Iso can increase i
Na and i
Ca
,
T which were reduced by the mutation. As a result of changes in the two currents, such abnormal pacemaking behavior could be eliminated.


## 4. Conclusions

In this study, by using a 2D anatomical model of the intact SAN-atrium, we have investigated the effects of Iso and the DelF1617 mutation on initiation and conduction of the pacemaker activity. The main findings of this study are as follows. (1) Iso could accelerate the pacemaking rate and affect the shape of action potentials in both the central and peripheral SAN cells, which are consistent with experimental observations [23, 24]. (2) The DelF1617 mutation slowed down pacemaking rate. Simulated addition of Ach amplified the bradycardic effect of the mutation and compromised the ability of the SAN to pace and drive the atrial septum. (3) When combined effects of Ach and the changes in the intrinsic properties of the SAN with the mutation were considered, conduction block can occur in both directions toward the atrium and CT, leading to a higher probability of SAN exit block and sinus arrest than with the mutation or Ach alone. This finding agrees with that of Smits et al. [34] and may provide further insight into the mechanism underlying high risk of cardiac arrest in SSS patients at night, when Ach concentration is high [35, 36]. (4) Abnormal pacemaking or even cessation state can be eliminated by proper concentration of Iso. Functional effects of Iso on the other mutations were similar (i.e., T220I, P1298L, and E161 K). What is more, Iso would be also useful for the treatment of bradycardia and improve atrial mechanical function [37]. It can be seen that these results are in good accordance with the clinical findings reported by Rokseth and Hatle [27], who provided the results of a 4-year prospective study regarding the occurrence, clinical pattern, and treatment of sinus arrest in patients with acute myocardial infarction. Taken together, we conclude that Iso could improve the pacemaking ability of the SAN cells as one of the fundamental medicines for SSS.



---

*Source: 540496-2014-01-21.xml*