420 likes | 436 Views
Basic Models in Theoretical Neuroscience. Oren Shriki 2010. Integrate and Fire and Conductance Based Neurons. References about neurons as electrical circuits:. Koch, C. Biophysics of Computation, Oxford Univ. Press, 1998.
E N D
Basic Models in Theoretical Neuroscience Oren Shriki 2010 Integrate and Fire and Conductance Based Neurons
References about neurons as electrical circuits: • Koch, C. Biophysics of Computation, Oxford Univ. Press, 1998. • Tuckwell, HC. Introduction to Theoretical Neurobiology, I&II, Cambridge UP, 1988.
Generation of Electric Potential on Nerve Cell Membranes • Chief factors that determine the resting membrane potential: • The relative permeability of the membrane to different ions • Differences in ionic concentrations • Ion pumps – Maintain the concentration gradient by actively moving ions against the gradient using metabolic resources. • Ion channels – “Holes” that allow the passage of ions in the direction of the concentration gradient. Some channels are selective for specific ions and some are not selective.
The Neuron as an Electric Circuit • Differences in ionic concentrations Battery • Cell membrane Capacitor • Ionic channels Resistors
The Neuron as an Electric Circuit Extracellular Intracellular
RC circuits • R – Resistance (in Ohms) • C – Capacitance (in Farads) Current source C I R
RC circuits • The dynamical equation is: C I R
RC circuits • Defining: • We obtain: • The general solution is:
RC circuit • Response to a step current:
RC circuit • Response to a step current:
The Integrate-and-Fire Neuron inside • R – Membrane Resistance (1/conductance) • C – Membrane Capacitance (in Farads) EL C I Threshold mechanism R outside
Integrate-and-Fire Neuron • If we define: • The dynamical equation will be: • To simplify, we define: • Thus:
Integrate-and-Fire Neuron • The threshold mechanism: • For V<θ the cell obeys its passive dynamics • For V=θ the cell fires a spike and the voltage resets to 0. • After voltage reset there is a refractory period, τR.
Integrate-and-Fire Neuron • Response to a step current: IR<θ: V t
Integrate-and-Fire Neuron • Response to a step current: IR>θ: V t τR τR τR T
Integrate-and-Fire Neuron • Finding the firing rate as a function of the applied current:
Each gating variable obeys the following dynamics: - Time constant - Represents the effect of temperature The Hodgkin & Huxley Framework
The Hodgkin & Huxley Framework The current through each channel has the form: - Maximal conductance (when all channels are open) - Fraction of open channels (can depend on several W variables).
The Temperature Parameter Φ • Allows for taking into account different temperatures. • Increasing the temperature accelerates the kinetics of the underlying processes. • However, increasing the temperature does not necessarily increase the excitability. Both increasing and decreasing the temperature can cause the neuron to stop firing. • A phenomenological model for Φ is:
Repetitive Firing in Hodgkin–Huxley Model A: Voltage time courses in response to a step of constant depolarizing current. from bottom to top: Iapp= 5, 15, 50, 100, 200 in μamp/cm2). Scale bar is 10 msec. B: f-I curves for temperatures of 6.3,18.5, 26◦C, as marked. Dotted curves show frequency of the unstable periodic orbits.
Fast-Slow Dissection of the Action Potential • n and h are slow compared to m and V. • Based on this observation, the system can be dissected into two time-scales. • This simplifies the analysis. • For details see:Borisyuk A & Rinzel J. Understanding neuronal dynamics by geometrical dissection of minimal models. In, Chow et al, eds: Models and Methods in Neurophysics (Les Houches Summer School 2003), Elsevier, 2005: 19-72.
Correlation between n and h • During the action potential the variables n and h vary together. • Using this correlation one can construct a reduced model. • The first to observe this was Fitzhugh.
Simplified Versions of the HH Model • Models that generate action potentials can be constructed with fewer dynamic variables. • These models are more amenable for analysis and are useful for learning the basic principles of neuronal excitability. • We will focus on the model developed by Morris and Lecar.
The Morris-Lecar Model (1981) • Developed for studying the barnacle muscle. • Model equations:
Morris-Lecar Model • The model contains K and Ca currents. • The variable w represents the fraction of open K channels. • The Ca conductance is assumed to behave in an instantaneous manner.
Morris-Lecar Model • A set of parameters for example:
Morris-Lecar Model • Voltage dependence of the various parameters (at long times):
Conductance-Based Models of Cortical Neurons • Cortical neurons behave differently than the squid axon that Hodgkin and Huxley investigated. • Over the years, people developed several variations of the HH model that are more appropriate for describing cortical neurons. • We will now see an example of a simple model which will later be useful in network simulations. • The model was developed by playing with the parameters such that its f-I curve is similar to that of cortical neurons.
Frequency-Current Responses of Cortical Neurons Excitatory Neuron: Ahmed et. al., Cerebral Cortex 8, 462-476, 1998 Inhibitory Neurons: Azouz et. al., Cerebral Cortex 7, 534-545, 1997
Frequency-Current Responses of Cortical Neurons The experimental findings show what f-I curves of cortical neurons are: • Continuous – starting from zero frequency. • Semi-Linear – above the threshold current the curve is linear on a wide range. How can we reconstruct this behavior in a model?
An HH Neuron with a Linear f-I Curve Shriki et al., Neural Computation 15, 1809–1841 (2003)
Linearization of the f-I Curve • We start with an HH neuron that has a continuous f-I curve (type I, saddle-node bifurcation). • The linearization is made possible by the addition of a certain K-current called A-current. • The curve becomes linear only when the time constant of the A-current is slow enough (~20 msec). • There are other mechanisms for linearizing f-I curves.
Model Equations: Shriki et al., Neural Computation 15, 1809–1841 (2003)