1 / 34

Neuronal excitability from a dynamical systems perspective

Mike Famulare CSE/NEUBEH 528 Lecture April 14, 2009. Neuronal excitability from a dynamical systems perspective. Outline. Analysis of biophysical models Hodgkin's classification of neurons by response to steady input currents Introduction to dynamical systems

aglaia
Download Presentation

Neuronal excitability from a dynamical systems perspective

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Mike Famulare CSE/NEUBEH 528 Lecture April 14, 2009 Neuronal excitability from a dynamical systems perspective

  2. Mike Famulare, CSE/NEUBEH 528 Outline Analysis of biophysical models Hodgkin's classification of neurons by response to steady input currents Introduction to dynamical systems Phase portraits and some bifurcation theory quadratic-integrate-and-fire model Fitzhugh-Nagmuo model General simple neuron models

  3. Mike Famulare, CSE/NEUBEH 528 Biophysical Modeling Neurons can be modelled with a set of nonlinear differential equations (Hodgkin-Huxley)‏

  4. Mike Famulare, CSE/NEUBEH 528 Hodgkin-Huxley Model

  5. What causes the spike in the HH model? • In response to a step current at t=5 ms: • fast inward current followed by slower outward current • sodium channel (m) activates • more slowly, potassium channel (n) activates and sodium (h) deinactivates. Mike Famulare, CSE/NEUBEH 528

  6. Mike Famulare, CSE/NEUBEH 528 Hodgkin-Huxley f-I curve Rate coding: firing rate response (f) to input current (I), steady state There is a minimum firing rate (58 Hz)‏ Can you infer the f-I curve from the Hodgkin-Huxley equations?

  7. Mike Famulare, CSE/NEUBEH 528 Connor-Stevens Model Model of a neuron in anisodoris (AKA the sea lemon nudibranch)‏

  8. Mike Famulare, CSE/NEUBEH 528 What causes a Connor-Stevens spike?

  9. Connor-Stevens f-I curve • Does not have a minimum firing frequency. Mike Famulare, CSE/NEUBEH 528

  10. Classifying Neurons by f-I type cortical pyramidal brainstem mesV Hodgkin's Classification of Neuronal Excitability Class 1: shows a continuous f-I curve (like Connor-Stevens)‏ Class 2: shows a discontinuous f-I curve (like Hodgkin-Huxley)‏ Class 3: shows no persistent firing (as can be found in auditory brainstem)‏ Mike Famulare, CSE/NEUBEH 528

  11. Mike Famulare, CSE/NEUBEH 528 Model-by-model is not the way to go! We want to understand why neurons are excitable. We want to understand what makes different neurons behave in different ways. Going model-by-model is difficult and not at all general: there are hundreds of channel types in nature any cell expresses a few or ten or so of them What do we do with cells whose response is measurable but for which we don't have a model?

  12. Mike Famulare, CSE/NEUBEH 528 Dynamical systems and simple models The models we've talked about are dynamical systems. What's a dynamical system? We can analyze dynamical systems to understand: equilibria (resting potentials)‏ “unstable manifolds” (spike thresholds)‏ bifurcations (f-I class, root of subthreshold behavior)‏

  13. Mike Famulare, CSE/NEUBEH 528 System: dynamical variables , control variable I. Fixed points : Linear response near a fixed point: Stability analysis: what are the eigenvalues of Bifurcation: change in the qualitative behavior of the system as the “control variable”, I, is changed. Overview of stability and bifurcation analysis

  14. Mike Famulare, CSE/NEUBEH 528 A note about the leaky-integrate-and-fire Leaky-integrate-and-fire (LIF) model:vo = resting potential, vr = reset voltage, vth = threshold voltage This model is not a spiking model in the sense that it doesn't have any dynamics for the spike itself. Piecewise “spike” is not dynamically similar to any real neuron Useful, but not for what we want to do today.

  15. Mike Famulare, CSE/NEUBEH 528 Quadratic-Integrate-and-Fire Model (QIF)‏ • Simplest model with dynamical spikes:

  16. Mike Famulare, CSE/NEUBEH 528 QIF fixed points Fixed points: or Do the fixed point exist for all I? for , the fixed points no longer exist(we'll come back to what this means soon!)

  17. Mike Famulare, CSE/NEUBEH 528 Phase plane analysis of the QIF: fixed points phase portrait for various external currents

  18. Mike Famulare, CSE/NEUBEH 528 Phase plane analysis of the QIF: stability Phase portrait of QIF: τm=1, α=1, and Iext=0

  19. Mike Famulare, CSE/NEUBEH 528 Stability Analysis of QIF linear response and stability at each fixed point for v- stable! v- is the resting potential for v+ unstable! v+ is the threshold voltage

  20. Mike Famulare, CSE/NEUBEH 528 Saddle-node bifurcation in the QIF Loss of stability via a saddle-node bifurcation: two fixed points “annihilate” each other also, we see that the QIF is an integrator!

  21. Mike Famulare, CSE/NEUBEH 528 Two types of saddle-node bifurcation “Saddle-node on invariant circle” (SNIC)‏ reset below “ghost” of fixed point arbitrarily low firing rate—Type I Saddle-node (SN)‏ reset above “ghost” fp slow first spike finite minimum firing rate—Type II

  22. Mike Famulare, CSE/NEUBEH 528 Summary of saddle-node bifurcations Saddle-node bifurcations occur when two fixed points disappear in response to a changing input Systems showing an SN bifurcation will be act as integrators For neurons, depending on details of the nonlinear spike return mechanism, SN bifurcators can be Type I (continuous f-I curve) or Type II (discontinuous f-I curve)‏ The Connor-Stevens model shows a SNIC bifurcation

  23. Mike Famulare, CSE/NEUBEH 528 What about resonating neurons? The saddle-node bifurcation type is only one of the very simple (“codimension 1”) bifurcations Hodgkin-Huxley does not show a saddle node bifurcation one of the many ways to see this is that the HH model cannot show an arbitrarily-long delayed first spike to step current With one dynamical variable, the saddle-node is the only possible continuous bifurcation, so we need two variables now!

  24. Mike Famulare, CSE/NEUBEH 528 Fitzhugh-Nagumo Model The Fitzhugh-Nagumo (FN) model is Hodgkin-Huxley like. Equations: Has two dynamical variables (is a two-dimensional dynamical system)‏ a voltage variable, v a recovery variable, w

  25. Mike Famulare, CSE/NEUBEH 528 Fitzhugh-Nagumo phase portrait For standard parameters, the FN has one fixed point that exists for all I. dx/dt = 0 lines are known as “nullclines”

  26. Mike Famulare, CSE/NEUBEH 528 Stability of the fixed point in the FN model finding the critical current bifurcation happens when intersection of nullclines is at the local minimum: fixed point nearby critical point

  27. Mike Famulare, CSE/NEUBEH 528 Linear Response of FN near critical point Linear response of 2D model: eigenvalues are a complex-conjugate pair stable when

  28. Mike Famulare, CSE/NEUBEH 528 FN model shows a Hopf bifurcation Hopf bifurcation: stable, oscillatory fixed point becomes an unstable, oscillatory fixed point there is a non-zero minimum firing rate controlled by the linear response frequency at the critical input current Type II excitability

  29. Mike Famulare, CSE/NEUBEH 528 Visualizing Hopf Dynamics phase portrait of a Hopf-bifurcating model (not the FN) for currents below the critical current looping around = subthreshold oscillation

  30. Mike Famulare, CSE/NEUBEH 528 Two types of Hopf bifurcation There are two types of Hopf bifurcation: supercritical (like Fitzhugh-Nagumo and HH)‏ subcritical see “Dynamical Systems in Neuroscience” by Izhikevich or Scholarpedia for details For real single neurons, the difference has never been found to be experimentally important (Izhikevich 2007) Both are Type II

  31. Mike Famulare, CSE/NEUBEH 528 Bifurcation Theory Review Bifurcation: a change in the qualitative behavior in response to a changing control parameter With only one control parameter (e.g. current), there are only two types of equilibrium bifurcations (“codimension one”): There is a lot more to this bifurcation business! what if you've got more inputs (drugs, hormones)? how do you fit a simple model (“normal form”, ”canonical model”) to a more complex model or real data?

  32. Mike Famulare, CSE/NEUBEH 528 Simple Models can cover a lot of ground Saddle-node and Hopf bifurcations are very common and can describe the single-spike properties of the spike-generating mechanisms of most neurons One model to do a lot (Izhikevich 2003)‏

  33. Mike Famulare, CSE/NEUBEH 528 Izhikevich's simple model (adaptive QIF)‏

  34. Mike Famulare, CSE/NEUBEH 528 References Dayan and Abbott “Dynamical Systems in Neuroscience” by E.M. Izhikevich “Spiking Neuron Models” by Gerstner and Kistler “Nonlinear Dynamics and Chaos” by Strogatz Scholarpedia

More Related