The Significance of Reduced Conductance in the Sino-Atrial Node for Cardiac Rythmicity Will Nesse May 14, 2005
Sino-Atrial Node (SA node): The Pacemaker of the Heart 1
The SA node and the Atria: How the pacemaker works Both the SA node and the Atria are excitable, spiking cells The SA node is more depolarized and has autonomous oscillations The Atria is more hyperpolarized and needs a stimulus to drive an action potential An action potential wave emanates from the SA node out through the Atria The Vagus nerve and the sympathetic cardiac nerve modulates the frequency of oscillation 2
The SA node Figure from Dobrzynski et al. 2005 3
The SA node Non contractile cells Small compared to the Atria Figure from Dobrzynski et al. 2005 4
The Atria Contractile cells Large compared to the SA node Figure from Dobrzynski et al. 2005t 5
The SA node is distinct from the Atria Figure from Dobrzynski et al. 2005 6
Electrical Conductance in the SA node- Atria The SA node expresses only low conductance gap junctions: 30-pS in SA node versus 160-pS in Atria Guess: The SA node has a lower electrical conductance 7
Question What is the Functional Significance of reduced conductance in the SA node? Dr Keener s result: There is a critical SA node size to achieve an oscillatory solution with homogeneous conductance Numerical studies (Joyner 1986, Winslow et al. 1995) suggest reduced conductance could help insulate the SA node from the Atria, reducing electrotonic loading from the Atria, preserving SA node oscillations My work: Reduced SA node conductance lowers the critical SA node size for an oscillatory solution 8
Why is a smaller SA node an advantage? Leaves more room for Atrial muscle Smaller SA node will have a more coherent oscillation Fewer cells Cells are more connected More robust to noise and intrinsic variation Easier to modulate 9
The Plan Use an analytically tractable 1D continuum model of excitable cardiac media: Piecewise linear FitzHugh-Nagumo model Model SA node as a small, excited part of the Atrial media Important parameters Excitability of SA node, Atria Conductance of SA node, Atria Size of SA node Specific aim: Find Hopf bifurcation condition for SA node- Atrial media relating SA node size and conductance 10
The Model The Continuum FitzHugh-Nagumo Model x [0, ) v = t x (k(x) v x ) + f(v) w w = ɛ(v γw α(x)) t k(x) = varying conductance α(x) = varying excitability level 11
Piecewise linear nonlinearity f(v) = v, v 1 4 v 1 2, 1 4 < v < 3 4 1 v, v 3 4 0.4 The FitzHugh Nagumo Piecewise Nonlinearity 0.3 0.2 f(v) 0.1 0 0.1 0.2 0.3 0.4 0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 v 12
Oscillatory solution to the uncoupled FitzHugh- Nagumo model 0.4 The FitzHugh Nagumo Equation 0.3 0.2 Recovery (w) 0.1 0 0.1 0.2 0.3 0.4 0.2 0 0.2 0.4 0.6 0.8 1 1.2 Voltage (V) 13
Voltage trace of uncoupled FitzHugh-Nagumo Model 1.2 The FitzHugh Nagumo Equation 1 0.8 Voltage (V) 0.6 0.4 0.2 0 0.2 0.4 0 5 10 15 20 25 30 35 40 45 50 Time 14
Piecewise constant excitability and conductance α(x) = { α s α a, x σ, x > σ k(x) = { k s k a, x σ, x > σ s = SA node a = Atria Where α s > α a and k s k a = 1 σ is the SA node boundary 15
Steady state: v t = w t = 0 0.4 SA Node Steady State Solution 0.35 0.3 Steady State Voltage 0.25 0.2 0.15 Region I II Region III 0.1 ξ σ 0.05 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x ξ is the boundary between excited (v > 1 4 ) and unexcited (v 1 4 ) states 16
Steady state solution V(x), W(x) The boundary and continuity conditions v I (0) = 0 v I (ξ) = v II (ξ) v I (ξ) = v II (ξ) v II (σ) = v III (σ) k s v II (σ) = k av III (σ) reduce to the consistency relation F (ξ, k s, σ) = 0 ξ(σ, k s ) 17
Steady state solution as a function of k s Steady State Solution as a Function of k s 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 0.2 0.4 0.6 k s 0.8 1 1 0.9 0.8 0.7 0.6 x 0.5 0.4 0.3 0.2 0.1 0 18
There is a larger excited region ξ as k s is reduced The Boundary of the Excited state For Different Values of k s 0.46 0.44 0.42 0.4 0.38 ξ 0.36 0.34 0.32 0.3 0.28 0.26 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 k s 19
Stability of steady state V (x), W (x) Linearize: v = V (x) + φ(x)e λt, w = W (x) + ψ(x)e λt Obtain Schrödinger equation µφ = x (k φ x ) + f (V (x))φ µ = λ + ɛ λ + ɛγ µ R and λ is pure imaginary (Hopf bifurcation) when µ = ɛγ 20
Schrödinger potential function f (V (x)) = Parameterized by ξ 1 x ξ 1 x > ξ 0.4 The FitzHugh Nagumo Piecewise Nonlinearity 0.3 0.2 f(v) 0.1 0 0.1 0.2 0.3 0.4 0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 v 21
Schrödinger equation solution Break solution up into the linear pieces, solve them and using the boundary data and continuity constraints φ I (0) = 0 φ I (ξ) = φ II (ξ) φ I (ξ) = φ II (ξ) φ II (σ) = φ III (σ) k s φ II (σ) = k aφ III (σ) lim III(x) x bounded generate transcendental consistency equation Φ(µ, ξ(k s, σ), k s ) = 0 22
Hopf bifurcation problem Set µ = ɛγ to find the bifurcation manifold for σ and k s : Φ(ɛγ, ξ(k s, σ), k s ) = 0 1.6 σ k s Bifurcation Manifold 1.4 1.2 1 Stable Steady State k s 0.8 0.6 0.4 0.2 Oscillatory Solution 0 0 0.5 1 1.5 2 2.5 3 σ 23
Future directions Investigate parameter effects on SA node oscillation frequency Investigate traveling wave solutions Will SA node emit a traveling pulse? How does it depend on parameters? Parameter-dependent dispersion relation of the Atrial media. Will it support high frequencies? Higher dimensions Explore numerically the existence of bifurcation manifold in more biologically realistic models 24