12.07.2015 Views

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

134 Conductance-Based Models• If one removes any current or gat<strong>in</strong>g variable, the model has only equilibriumattractors for any values of parameters.We refer to such a model as be<strong>in</strong>g m<strong>in</strong>imal or irreducible for spik<strong>in</strong>g. Thus, m<strong>in</strong>imalmodels can exhibit periodic activity, even if of small amplitude, but their reductionscannot. Accord<strong>in</strong>g to this def<strong>in</strong>ition, any space-clamped conductance-based model iseither a m<strong>in</strong>imal model, or could be reduced to a m<strong>in</strong>imal model or models by remov<strong>in</strong>ggat<strong>in</strong>g variables. This will be the basis for our classification of electrophysiologicalmechanisms <strong>in</strong> neurons.For example, the Hodgk<strong>in</strong>-Huxley model considered <strong>in</strong> Sect. 2.3 is not m<strong>in</strong>imal forspik<strong>in</strong>g. Recall that the model consists of three currents: leakage I L , transient sodiumI Na,t (gat<strong>in</strong>g variables m and h) and persistent potassium I K (gat<strong>in</strong>g variable n); seeFig. 5.1. Which of these currents are responsible for excitability and spik<strong>in</strong>g?We can remove the leakage current and the gat<strong>in</strong>g variable, h, for the <strong>in</strong>activationof the sodium current: The result<strong>in</strong>g I Na,p +I K -modelC ˙V = I −I K{ }} {g K n 4 (V − E K ) −ṅ = (n ∞ (V ) − n)/τ n (V ) ,ṁ = (m ∞ (V ) − m)/τ m (V ) ,I Na,p{ }} {g Na m 3 (V − E Na ) ,was considered <strong>in</strong> the previous chapter where we have shown that it could oscillate dueto the <strong>in</strong>terplay between the activation of persistent sodium and potassium currents.Alternatively, we can remove the K + current from the Hodgk<strong>in</strong>-Huxley model, yet thenew I Na,t -modelC ˙V = I −I Na,t{ }} {g Na m 3 h(V − E Na ) −ṁ = (m ∞ (V ) − m)/τ m (V ) ,ḣ = (h ∞ (V ) − h)/τ h (V ) ,I L{ }} {g L (V − E L ) ,can still oscillate via the <strong>in</strong>terplay between activation and <strong>in</strong>activation of the Na +current, as we will see later <strong>in</strong> this chapter. Both models are m<strong>in</strong>imal, because removalof any other gat<strong>in</strong>g variable results <strong>in</strong> either the I Na,p -, I K -, or I h -models, neither ofwhich can have a limit cycle attractor, as the reader is asked to prove at the end of theprevious chapter.We see that the Hodgk<strong>in</strong>-Huxley model is not m<strong>in</strong>imal, but it is a comb<strong>in</strong>ation oftwo m<strong>in</strong>imal models. M<strong>in</strong>imal models are appeal<strong>in</strong>g because they are relatively simple;each <strong>in</strong>dividual variable has an established electrophysiological mean<strong>in</strong>g, and its role<strong>in</strong> dynamics can be easily identified. As we show below, many m<strong>in</strong>imal models can

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

Saved successfully!

Ooh no, something went wrong!