12.07.2015 Views

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

150 Conductance-Based Modelshave comparable time constants or <strong>in</strong>activation is much faster than activation.Even though none of the experimentally measured A-currents show fast <strong>in</strong>activationand a relatively slower activation, this case is still <strong>in</strong>terest<strong>in</strong>g from the pure theoreticalpo<strong>in</strong>t of view, s<strong>in</strong>ce it shows how a s<strong>in</strong>gle K + current can give rise to oscillations.Assum<strong>in</strong>g <strong>in</strong>stantaneous <strong>in</strong>activation and us<strong>in</strong>g h = h ∞ (V ) <strong>in</strong> the voltage equation, weobta<strong>in</strong> a two-dimensional systemC ˙V = I −leak I L{ }} {g L (V − E L ) −ṁ = (m ∞ (V ) − m)/τ m (V ) ,I A , <strong>in</strong>st. <strong>in</strong>activation{ }} {g A m h ∞ (V )(V − E K ) ,whose nullcl<strong>in</strong>es can easily be found:m = I − g L(V −E L )g A h ∞ (V )(V − E K )(V -nullcl<strong>in</strong>e)andm = m ∞ (V ) (m-nullcl<strong>in</strong>e) .Two typical cases are depicted <strong>in</strong> Fig. 5.14a and b. We start with the simpler case <strong>in</strong>Fig. 5.14b.Figure 5.14b depicts nullcl<strong>in</strong>es when the A-current has low activation threshold.There is only one <strong>in</strong>tersection of the nullcl<strong>in</strong>es, hence there is only one equilibrium,which is a stable focus when the <strong>in</strong>jected dc-current I is not strong enough (Fig. 5.14b,top). Increas<strong>in</strong>g I makes the equilibrium lose stability via a supercritical Andronov-Hopf bifurcation that gives birth to a small amplitude limit cycle attractor (not shown<strong>in</strong> the figure). A further <strong>in</strong>crease of I <strong>in</strong>creases the amplitude of oscillations, and e.g.,when I = 10 (middle of Fig. 5.14b), the attractor corresponds to periodic fir<strong>in</strong>g of actionpotentials. When I = 10.5, the attractor disappears and the equilibrium becomesstable (via Andronov-Hopf bifurcation) aga<strong>in</strong>. The model, however, becomes excitable.A small hyperpolarization does not change significantly the A-current, and the voltagereturns back to rest result<strong>in</strong>g <strong>in</strong> a “subthreshold response”. A sufficiently largehyperpolarization de<strong>in</strong>activates enough I A to open the K + current and hyperpolarizethe membrane even further. This regenerative process produces the downstroke andbr<strong>in</strong>gs V close to E K . Dur<strong>in</strong>g the state of hyperpolarization, the A-current deactivates(m → 0), and the dc-current I br<strong>in</strong>gs V back to rest. Notice that the action potentialis directed downward.In Fig. 5.14a we consider the I A -model with exactly the same parameters exceptthat we shift the half-voltage activation V 1/2 of I A by 10 mV so that the A-currenthas higher activation threshold. This does not affect much the behavior of the systemwhen I is small. However, when I ≈ 10.7 the spik<strong>in</strong>g limit cycle attractor undergoesa new k<strong>in</strong>d of bifurcation — saddle-node bifurcation — result<strong>in</strong>g <strong>in</strong> the appearance oftwo new equilibria: a stable node and a saddle. If the reader looks at Fig. 5.14a upsidedown,he will notice that this figure resembles figures 5.4a, 5.6a, or 5.12a, with all the

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

Saved successfully!

Ooh no, something went wrong!