Entrainment and Chaos in the Hodgkin-Huxley Oscillator Kevin K. Lin http://www.cims.nyu.edu/ klin Courant Institute, New York University Mostly Biomath - 2005.4.5 p.1/42
Overview (1) Goal: Show that the Hodgkin-Huxley neuron model, driven by a periodic impulse train, can exhibit entrainment, transient chaos, and fully chaotic behavior. Mostly Biomath - 2005.4.5 p.2/42
Overview (2) Why? 1. Suggested by general, rigorous perturbation theory of kicked oscillators (Qiudong Wang & Lai-Sang Young). 2. Hodgkin-Huxley is a paradigm for excitable biological systems where pulse-like inputs and outputs are natural. Mostly Biomath - 2005.4.5 p.3/42
Overview (3) Results: 1. Entrainment and chaos are readily observable in the pulse-driven Hodgkin-Huxley system. 2. The pulse-driven Hodgkin-Huxley system prefers entrainment. 3. Strong expansion is caused by invariant structures. Mostly Biomath - 2005.4.5 p.4/42
Outline Classical Hodgkin-Huxley neuron model Kicked nonlinear oscillators & Wang-Young theory Pulse-driven Hodgkin-Huxley neuron model Mostly Biomath - 2005.4.5 p.5/42
Squid giant axon http://hermes.mbl.edu/publications/loligo/squid Mostly Biomath - 2005.4.5 p.6/42
Schematic (rest state) References: Abbott and Dayan, Theoretical Neuroscience, MIT Press 2001 Cronin, Mathematical Aspects of Hodgkin-Huxley Neural Theory, Cambridge University Press 1987 Mostly Biomath - 2005.4.5 p.7/42
Hodgkin-Huxley equations v = C 1 [ I ḡ K n 4 (v v K ) ḡ Na m 3 h(v v Na ) ḡ leak (v v leak ) ] ṁ = αm(v)(1 m) βm(v)m ṅ = αn(v)(1 n) βn(v)n ḣ = α h (v)(1 h) β h (v)h αm(v) = Ψ ( ) v+25 10, β m(v) = 4 exp (v/18), αn(v) = 0.1Ψ ( ) v+10 10, β n(v) = 0.125 exp (v/80), 1 α h (v) = 7 exp (v/20), β h (v) = 1+exp( v+30 10 ), v Ψ(v) = exp(v) 1. Mostly Biomath - 2005.4.5 p.8/42
Equivalent circuit [ ] v = C 1 I ḡ K n 4 (v v K ) ḡ Na m 3 h(v v Na ) ḡ leak (v v leak ) http://www.syssim.ecs.soton.ac.uk/vhdl-ams/examples/hodhuxneu/hh2.htm Mostly Biomath - 2005.4.5 p.9/42
Parameters In this study: All parameters take on original squid values except the injected current I This guarantees stable oscillations Mostly Biomath - 2005.4.5 p.10/42
Parameters (cont d) v Na = 115 mv, ḡ Na = 120 mω 1 /cm 2, v K = +12 mv, ḡ K = 36 mω 1 /cm 2, v leak = 10.613 mv, ḡ leak = 0.3 mω 1 /cm 2, C = 1 µf/cm 2, I = 14.2211827403 Mostly Biomath - 2005.4.5 p.11/42
Parameters (cont d) 20 0-20 Stable fixed point Unstable fixed point V -40-60 Unstable cycle -80 Limit cycle -100 6 7 8 9 10 11 12 13 14 -I Mostly Biomath - 2005.4.5 p.12/42
Dynamics without kicks 0.7-2 0.65 v -4 n 0.6 0.55-6 0.5-8 0.45 1 2 t 3 4-8 -6-4 v -2 Mostly Biomath - 2005.4.5 p.13/42
Outline Classical Hodgkin-Huxley neuron model Kicked nonlinear oscillators & Wang-Young theory Pulse-driven Hodgkin-Huxley neuron model Mostly Biomath - 2005.4.5 p.14/42
Kicked oscillators A stable, nonlinear oscillator is a flow with a limit cycle γ (period=t γ ) and basin of attraction U. A kick instantaneously changes the system s state: Mostly Biomath - 2005.4.5 p.15/42
Examples of kicked oscillators Circadian rhythm, phase reset experiments (Winfree). Possible approach to disrupting synchronous firing of neuron (Tass). Mostly Biomath - 2005.4.5 p.16/42
Simple Example 1.0 y -1.0-1.0 x 1.0 2.0 ṙ = (µ αr 2 )r+ 1 2 sin(3θ) δ(t nt) n Z θ = ω + βr 2 Mostly Biomath - 2005.4.5 p.17/42
Effect of Kick-and-Flow on Phase Space 1.0 0.8 1.0 0.6 0.5 0.4 0.5 0.2-0.2-0.5-0.4-0.5-0.6-1.0-0.8-1.0-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1.0-1.5-1.0-0.5 0.5 1.0-1.0-0.5 0.5 1.0 t = 0 t = 0 (after kick) t = 1 1.0 1.0 1.0 0.5 0.5 0.5-0.5-0.5-1.0-0.5-1.0-1.5-1.0-1.0-0.5 0.5 1.0-1.0 1.0-1.0-0.5 0.5 1.0 t = 2 t = 2 (after kick) t = 4 Mostly Biomath - 2005.4.5 p.18/42
Discrete time map Define F T : R 4 R 4 by F T (x) = Φ T (K(x)), where K(x) represents kicks Φ T (x) = flow map T = period of kicks. Continuous time Discrete time Entrainment F T has sinks Chaos F T chaotic Mostly Biomath - 2005.4.5 p.19/42
Reduction to 1-D Wang and Young start with F T and 1. Reduces from the map F T on R n to a circle map f T : lim F T+nT γ (x) n induces a map f T on γ S 1. We refer to f T as the singular limit or the phase resetting curve. 2. Analyze f T and infer properties of F T. Mostly Biomath - 2005.4.5 p.20/42
Wang-Young Conditions If 1. Kicks do not send limit cycle to bad places, i.e. K(γ) does not go outside the basin of γ 2. Kicks are in the right directions (e.g. not along W ss (x)) to take advantage of shear Mostly Biomath - 2005.4.5 p.21/42
Wang-Young Consequences Then for different kick amplitude A & kick period T the discrete-time system F T can have 1. Rotation-like behavior (small A) 2. Sinks and sources (for all A large enough) 3. Transient chaos / horseshoes (for large A & T) 4. Strange attractor & chaos (for large A & T 1) Mostly Biomath - 2005.4.5 p.22/42
Wang-Young Theory (cont d) Notes: 1. The conditions are satisfied for the simple example. 2. For Hodgkin-Huxley there is not too much choice if we want to stay close to physical interpretation of model. Mostly Biomath - 2005.4.5 p.23/42
Lyapunov exponents The largest Lyapunov exponent λ of F T is useful for distinguishing different scenarios numerically: Rotations λ = 0 Sinks λ < 0 Chaos λ > 0 Mostly Biomath - 2005.4.5 p.24/42
Outline Classical Hodgkin-Huxley neuron model Kicked nonlinear oscillators & Wang-Young theory Pulse-driven Hodgkin-Huxley neuron model Mostly Biomath - 2005.4.5 p.25/42
Hodgkin-Huxley equations v = C 1 [ I ḡ K n 4 (v v K ) ḡ Na m 3 h(v v Na ) ḡ leak (v v leak ) ] +A n Z δ (t nt) ṁ = αm(v)(1 m) βm(v)m ṅ = αn(v)(1 n) βn(v)n ḣ = α h (v)(1 h) β h (v)h Prior work: Winfree, Best on null space and degree-change. Mostly Biomath - 2005.4.5 p.26/42
Entrainment A = 10, T = 81.0-2 v -4-6 -8 5 t 10 15 6 4 2 v -2-4 -6 50 t 100 150 Mostly Biomath - 2005.4.5 p.27/42
Entrainment (cont d) -2 A = 10, T = 81.0 v -4-6 -8 5 t 10 15 10 5 v1(t)-v2(t) -5-10 50 t 100 150 Mostly Biomath - 2005.4.5 p.28/42
Entrainment (cont d) Time-T map: F T = Φ T K 2 v1(n)-v2(n) -2-4 1 Multiple of T (n) 2 Mostly Biomath - 2005.4.5 p.29/42
Chaos v 10 8 6 4 2-2 -4-6 -8 A = 10, T = 80.8 200 400 600 t 2.0 Log10(dist) -2.0-4.0-6.0 200 400 600 t Mostly Biomath - 2005.4.5 p.30/42
λ (F T ) versus T Largest Lyapunov exponent of F_T -2.0-4.0-6.0-8.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 T / T_gamma Mostly Biomath - 2005.4.5 p.31/42
λ (F T ) versus A 1.0 0.8 ROTATIONS 0.6 SINKS 0.4 0.2 CHAOS UNKNOWN 1 2 3 4 Drive amplitude Chaos: Prob( ˆλ > 3ǫ) Rotations: Prob( ˆλ < ǫ/3) Sinks: Prob( ˆλ < 3ǫ) Unknown: everything else Mostly Biomath - 2005.4.5 p.32/42
Phase resetting curves ( f T ) A=5, T=101.5 A=10, T=80.8 15.0 2 15.0 1 1 5.0 5.0 2.5 5.0 7.5 1 12.5 2.5 5.0 7.5 1 12.5 A=20, T=101.5 15.0 1 5.0 2.5 5.0 7.5 1 12.5 Mostly Biomath - 2005.4.5 p.33/42
Plateau and fixed points The first return map R ft to the interval [4, 10] (enclosing the plateau), for A = 10 and T = 17.6. 1 9.0 8.0 7.0 6.0 5.0 4.0 4.0 5.0 6.0 7.0 8.0 9.0 1 Mostly Biomath - 2005.4.5 p.34/42
Plateau and fixed points 1 9.0 8.0 7.0 6.0 5.0 4.0 4.0 5.0 6.0 7.0 8.0 9.0 1 Drive amplitude A Probability of sink near plateau 5 41% 10 58% 20 68% 30 76% Mostly Biomath - 2005.4.5 p.35/42
Zooming into the kink 16.0 14.0 12.0 1 8.0 6.0 4.0 2.0 0.2 0.4 0.6 0.8 1.0 Mostly Biomath - 2005.4.5 p.36/42
Why the kink? 0.38 0.36 0.45 0.44 h 0.34 n 0.43 0.32 0.42-14.0-12.0-1-8.0-6.0-4.0-2.0 0.1 0.12 0.14 0.16 0.18 v m But Hodgkin-Huxley lives in R 4... Mostly Biomath - 2005.4.5 p.37/42
Why the kink? (cont d) Approaching critical A crit 13.5895...: A=13.58 A=13.589 A=13.5895 3 4 5 4 3 2 3 2 2 1 1 1 9.0 9.5 1 10.5 9.0 9.5 1 10.5 9.0 9.5 1 10.5 A=13.5896 A=13.59 A=13.6 4 3 2 3 2 2 1 1 1 9.0 9.5 1 10.5 9.0 9.5 1 10.5 9.0 9.5 1 10.5 Mostly Biomath - 2005.4.5 p.38/42
Why the plateau? Graph of f T for A = 10, around plateau. 9.0 8.0 f_t 7.0 6.0 5.0 5.0 6.0 7.0 phase 8.0 9.0-1.0-2.0-3.0-4.0-5.0 5.0 6.0 7.0 phase 8.0 9.0 Black: log 10 f Blue: log 10 (Ess(K(γ(θ)), γ(θ))) Mostly Biomath - 2005.4.5 p.39/42
Finding horseshoes Horseshoes can produce transient chaos: A = 10, T = 81 12.0 1 8.0 6.0 4.0 2.0 0.3 0.4 0.5 0.6 0.7 Mostly Biomath - 2005.4.5 p.40/42
Summary Can observe entrainment and chaos in the pulse-driven Hodgkin-Huxley neuron model. The pulse-driven Hodgkin-Huxley model prefers entrainment. This can be explained. Complex phase response can arise from kicks going near invariant structures. Mostly Biomath - 2005.4.5 p.41/42
References 1. Eric N. Best, Null space in the Hodgkin-Huxley equations, Biophys. J. 27 (1979) 2. Kevin K. Lin, Entrainment and chaos in the Hodgkin-Huxley oscillator, in preparation 3. Qiudong Wang and Lai-Sang Young, Strange attractors in periodically-kicked limit cycles and Hopf bifurcations, Comm. Math. Phys. 240 (2003) 4. Arthur Winfree, The Geometry of Biological Time, 2nd Edition, Springer-Verlag (2000) Acknowledgements I am grateful to Lai-Sang Young for her help with this work, and to Eric Brown, Adi Rangan, Alex Barnett, and Toufic Suidan for critical comments. Many thanks to Charlie Peskin for the invitation. This work is supported by the National Science Foundation. Mostly Biomath - 2005.4.5 p.42/42