Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience: Dynamical Systems in Neuroscience:

ib.cnea.gov.ar
from ib.cnea.gov.ar More from this publisher
12.07.2015 Views

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

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!