Dynamical Systems in Neuroscience:
Dynamical Systems in Neuroscience: Dynamical Systems in Neuroscience:
148 Conductance-Based Models1I=651I=700.80.8V-nullclineV-nullclineK + activation, n0.60.40.2n-nullcline0 50K + activation, n0.60.40.2n-nullcline000 50-80 -60 -40 -20 0 20membrane voltage, V-80 -60 -40 -20 0 20membrane voltage, V1I=681I=73V-nullcline0.80.8V-nullclineK + activation, n0.60.40.2n-nullcline0 100K + activation, n0.60.40.2n-nullcline000 100-80 -60 -40 -20 0 20membrane voltage, Va-80 -60 -40 -20 0 20membrane voltage, VbFigure 5.12: Possible intersections of nullclines in the I K +I Kir -model. Parameters:E K = −80 mV, g Kir = 20, g K = 2. Instantaneous I Kir with V 1/2 = −80 mV andk = −12. Slower I K with k = 5, τ(V ) = 5 ms, and V 1/2 = −40 mV (in a) orV 1/2 = −55 mV (in b).potential. If we increase I, the node and the saddle approach, coalesce, and annihilateeach other via a saddle-node bifurcation, and the model start to fire action potentialsperiodically.We see that I K +I Kir -model has essentially the same dynamic repertoire as the moreconventional I Na,p +I K -model or I Na,t -model, despite the fact that it is based on a ratherbizarre ionic mechanism for excitability and spiking.5.1.7 I A -modelThe last minimal voltage-gated model has only one transient K + current, often referredto as being A-current I A , yet it can also generate sustained oscillations. In some sense,the model is similar to the I Na,t -model. Indeed, each consists of only one transient
Conductance-Based Models 149de-inactivationof I Ahyperpolarizationdeactivationof I Ainjecteddc-currentactivation +inactivationof I Aresting potentialFigure 5.13: Mechanism of generation of sustainedvoltage oscillations in the I A -model.current and an Ohmic leak current. The only difference is that A-current is outward,and as a result, the action potentials are fired downward; see Fig. 5.14 and Fig. 5.15below.The A-current has activation and inactivation variables m and h, respectively, andthe model has the formC ˙V = I −leak I L{ }} {g L (V − E L ) −ṁ = (m ∞ (V ) − m)/τ m (V )ḣ = (h ∞ (V ) − h)/τ h (V ) .I A{ }} {g A mh(V − E K )The mechanism of generation of downward action potentials is summarized in Fig. 5.13.Due to a strong injected dc-current, the rest state is in the depolarizing voltage range,and it corresponds to the balance of the partially activated, partially inactivated A-current, leak outward current, and the injected dc-current. A small hyperpolarizationcan simultaneously deactivate and deinactivate the A-current, i.e., decrease variablem and increase variable h. Depending on the relationship between the activation andinactivation time constants, this may result in an increase of the A-current conductance,which is proportional to the product mh. More outward current produces morehyperpolarization and even more outward current. As a result of this regenerativeprocess, the membrane voltage produces a sudden downstroke. While hyperpolarized,the A-current deactivates (variable m → 0), and the injected dc-current slowly bringsthe membrane potential toward the resting state, resulting in a slow upstroke. Fastdownstroke and a slower upstroke from a depolarized resting state look like an actionpotential pointing downwards.If activation kinetics is much faster than the inactivation kinetics, we can substitutem = m ∞ (V ) into the voltage equation above and reduce the I A -model to atwo-dimensional system, which hopefully would have the right kind of nullclines anda limit cycle attractor. After all, this is what we have done with previous minimalmodels and it always worked. As the reader is asked to prove in Ex. 1, the I A -modelcannot have a limit cycle attractor when the A-current activation kinetics is instantaneous.Oscillations are possible only when the activation and inactivation kinetics
- Page 108 and 109: 98 Two-Dimensional Systems(f(x(t),y
- Page 110 and 111: 100 Two-Dimensional Systems0.70acti
- Page 112 and 113: 102 Two-Dimensional Systems0.70.6K
- Page 114 and 115: 104 Two-Dimensional Systemsy1.510.5
- Page 116 and 117: 106 Two-Dimensional SystemsFor exam
- Page 118 and 119: 108 Two-Dimensional SystemsIn gener
- Page 120 and 121: 110 Two-Dimensional Systemsv 2 v 2v
- Page 122 and 123: 112 Two-Dimensional SystemsV-nullcl
- Page 124 and 125: 114 Two-Dimensional SystemsV-nullcl
- Page 126 and 127: 116 Two-Dimensional Systemsheterocl
- Page 128 and 129: 118 Two-Dimensional Systems0.60.5n-
- Page 130 and 131: 120 Two-Dimensional Systemseigenval
- Page 132 and 133: 122 Two-Dimensional SystemsK + acti
- Page 134 and 135: 124 Two-Dimensional Systemsexcitati
- Page 136 and 137: 126 Two-Dimensional SystemsReview o
- Page 138 and 139: 128 Two-Dimensional SystemsabcdFigu
- Page 140 and 141: 130 Two-Dimensional SystemsFigure 4
- Page 142 and 143: 132 Two-Dimensional Systems
- Page 144 and 145: 134 Conductance-Based Models• If
- Page 146 and 147: 136 Conductance-Based Modelsinward(
- Page 148 and 149: 138 Conductance-Based Models1I=01I=
- Page 150 and 151: 140 Conductance-Based Modelsleak cu
- Page 152 and 153: 142 Conductance-Based Modelshyperpo
- Page 154 and 155: 144 Conductance-Based Models0I=0ina
- Page 156 and 157: 146 Conductance-Based Models0I=9.75
- Page 160 and 161: 150 Conductance-Based Modelshave co
- Page 162 and 163: 152 Conductance-Based Models1 sec10
- Page 164 and 165: 154 Conductance-Based Modelsvoltage
- Page 166 and 167: 156 Conductance-Based Models1h=0.89
- Page 168 and 169: 158 Conductance-Based Models5.2.2 E
- Page 170 and 171: 160 Conductance-Based Models1recove
- Page 172 and 173: 162 Conductance-Based Modelscurrent
- Page 174 and 175: 164 Conductance-Based Modelsmodels
- Page 176 and 177: 166 Conductance-Based Models
- Page 178 and 179: 168 Bifurcationssaddle-node bifurca
- Page 180 and 181: 170 Bifurcationssaddle-node bifurca
- Page 182 and 183: 172 Bifurcations0.5v 2V-nullclineK
- Page 184 and 185: 174 BifurcationsBoth types of the b
- Page 186 and 187: 176 BifurcationsV-nullcline0.5K + g
- Page 188 and 189: 178 Bifurcationsrstable limit cycle
- Page 190 and 191: 180 BifurcationsnVstable limit cycl
- Page 192 and 193: ¡ ¡¡ ¡¡ ¡ ¡¡ ¡ ¡182 Bifur
- Page 194 and 195: ¡ ¡ ¡¡ ¡ ¡¡ ¡ ¡¡ ¡ ¡184
- Page 196 and 197: 186 Bifurcationsmembrane potential,
- Page 198 and 199: 188 BifurcationsBifurcation of a li
- Page 200 and 201: 190 Bifurcationsstableunstablelimit
- Page 202 and 203: 192 Bifurcationsamplitude (max-min)
- Page 204 and 205: 194 Bifurcationshomoclinic orbithom
- Page 206 and 207: 196 Bifurcations0.80.60.4stable lim
Conductance-Based Models 149de-<strong>in</strong>activationof I Ahyperpolarizationdeactivationof I A<strong>in</strong>jecteddc-currentactivation +<strong>in</strong>activationof I Arest<strong>in</strong>g potentialFigure 5.13: Mechanism of generation of susta<strong>in</strong>edvoltage oscillations <strong>in</strong> the I A -model.current and an Ohmic leak current. The only difference is that A-current is outward,and as a result, the action potentials are fired downward; see Fig. 5.14 and Fig. 5.15below.The A-current has activation and <strong>in</strong>activation variables m and h, respectively, andthe model has the formC ˙V = I −leak I L{ }} {g L (V − E L ) −ṁ = (m ∞ (V ) − m)/τ m (V )ḣ = (h ∞ (V ) − h)/τ h (V ) .I A{ }} {g A mh(V − E K )The mechanism of generation of downward action potentials is summarized <strong>in</strong> Fig. 5.13.Due to a strong <strong>in</strong>jected dc-current, the rest state is <strong>in</strong> the depolariz<strong>in</strong>g voltage range,and it corresponds to the balance of the partially activated, partially <strong>in</strong>activated A-current, leak outward current, and the <strong>in</strong>jected dc-current. A small hyperpolarizationcan simultaneously deactivate and de<strong>in</strong>activate the A-current, i.e., decrease variablem and <strong>in</strong>crease variable h. Depend<strong>in</strong>g on the relationship between the activation and<strong>in</strong>activation time constants, this may result <strong>in</strong> an <strong>in</strong>crease of the A-current conductance,which is proportional to the product mh. More outward current produces morehyperpolarization and even more outward current. As a result of this regenerativeprocess, the membrane voltage produces a sudden downstroke. While hyperpolarized,the A-current deactivates (variable m → 0), and the <strong>in</strong>jected dc-current slowly br<strong>in</strong>gsthe membrane potential toward the rest<strong>in</strong>g state, result<strong>in</strong>g <strong>in</strong> a slow upstroke. Fastdownstroke and a slower upstroke from a depolarized rest<strong>in</strong>g state look like an actionpotential po<strong>in</strong>t<strong>in</strong>g downwards.If activation k<strong>in</strong>etics is much faster than the <strong>in</strong>activation k<strong>in</strong>etics, we can substitutem = m ∞ (V ) <strong>in</strong>to the voltage equation above and reduce the I A -model to atwo-dimensional system, which hopefully would have the right k<strong>in</strong>d of nullcl<strong>in</strong>es anda limit cycle attractor. After all, this is what we have done with previous m<strong>in</strong>imalmodels and it always worked. As the reader is asked to prove <strong>in</strong> Ex. 1, the I A -modelcannot have a limit cycle attractor when the A-current activation k<strong>in</strong>etics is <strong>in</strong>stantaneous.Oscillations are possible only when the activation and <strong>in</strong>activation k<strong>in</strong>etics