Quantitative models of the action potential
In
Hodgkin–Huxley model
In 1952
where INa, IK, and IL are currents conveyed through the local sodium channels, potassium channels, and "leakage" channels (a catch-all), respectively. The initial term Iext represents the current arriving from external sources, such as excitatory postsynaptic potentials from the dendrites or a scientist's electrode.
The model further assumes that a given ion channel is either fully open or closed; if closed, its
which is equivalent to Ohm's law. By definition, no net current flows (IK = 0) when the transmembrane voltage equals the equilibrium voltage of that ion (when V = EK).
To fit their data accurately, Hodgkin and Huxley assumed that each type of ion channel had multiple "gates", so that the channel was open only if all the gates were open and closed otherwise. They also assumed that the probability of a gate being open was independent of the other gates being open; this assumption was later validated for the inactivation gate.[9] Hodgkin and Huxley modeled the voltage-sensitive potassium channel as having four gates; letting pn denote the probability of a single such gate being open, the probability of the whole channel being open is the product of four such probabilities, i.e., popen, K = n4. Similarly, the probability of the voltage-sensitive sodium channel was modeled to have three similar gates of probability m and a fourth gate, associated with inactivation, of probability h; thus, popen, Na = m3h. The probabilities for each gate are assumed to obey first-order kinetics
where both the equilibrium value meq and the relaxation time constant τm depend on the instantaneous voltage V across the membrane. If V changes on a time-scale more slowly than τm, the m probability will always roughly equal its equilibrium value meq; however, if V changes more quickly, then m will lag behind meq. By fitting their voltage-clamp data, Hodgkin and Huxley were able to model how these equilibrium values and time constants varied with temperature and transmembrane voltage.[1] The formulae are complex and depend exponentially on the voltage and temperature. For example, the time constant for sodium-channel activation probability h varies as 3(θ−6.3)/10 with the Celsius temperature θ, and with voltage V as
In summary, the Hodgkin–Huxley equations are complex, non-linear
FitzHugh–Nagumo model
Because of the complexity of the Hodgkin–Huxley equations, various simplifications have been developed that exhibit qualitatively similar behavior.[3][13] The FitzHugh–Nagumo model is a typical example of such a simplified system.[14][15] Based on the tunnel diode, the FHN model has only two independent variables, but exhibits a similar stability behavior to the full Hodgkin–Huxley equations.[16] The equations are
where g(V) is a function of the voltage V that has a region of negative slope in the middle, flanked by one maximum and one minimum (Figure FHN). A much-studied simple case of the FitzHugh–Nagumo model is the Bonhoeffer-van der Pol nerve model, which is described by the equations[17]
where the coefficient ε is assumed to be small. These equations can be combined into a second-order differential equation
This van der Pol equation has stimulated much research in the mathematics of nonlinear dynamical systems. Op-amp circuits that realize the FHN and van der Pol models of the action potential have been developed by Keener.[18]
A hybrid of the Hodgkin–Huxley and FitzHugh–Nagumo models was developed by Morris and Lecar in 1981, and applied to the muscle fiber of barnacles.[19] True to the barnacle's physiology, the Morris–Lecar model replaces the voltage-gated sodium current of the Hodgkin–Huxley model with a voltage-dependent calcium current. There is no inactivation (no h variable) and the calcium current equilibrates instantaneously, so that again, there are only two time-dependent variables: the transmembrane voltage V and the potassium gate probability n. The bursting, entrainment and other mathematical properties of this model have been studied in detail.[20]
The simplest models of the action potential are the "flush and fill" models (also called "integrate-and-fire" models), in which the input signal is summed (the "fill" phase) until it reaches a threshold, firing a pulse and resetting the summation to zero (the "flush" phase).
Extracellular potentials and currents
Whereas the above models simulate the transmembrane voltage and current at a single patch of membrane, other mathematical models pertain to the voltages and currents in the ionic solution surrounding the neuron.
where j and E are vectors representing the
where the integration is over the complete surface of the membrane; is a position on the membrane, σinside and φinside are the conductivity and potential just within the membrane, and σoutside and φoutside the corresponding values just outside the membrane. Thus, given these σ and φ values on the membrane, the extracellular potential φ(x) can be calculated for any position x; in turn, the electric field E and current density j can be calculated from this potential field.[25]
See also
- Biological neuron models
- GHK current equation
- Models of neural computation
- Saltatory conduction
- Bioelectronics
- Cable theory
References
- ^ PMID 12991237.
- ISBN 978-0-262-11133-1.)
{{cite book}}
: CS1 maint: multiple names: authors list (link - ^ ISBN 978-0-521-31574-6.
- S2CID 41760422.
- S2CID 10775478.
- ISBN 978-0-262-11133-1.
- ^ Hooper, Scott L. "Central Pattern Generators." Embryonic ELS (1999) http://www.els.net/elsonline/figpage/I0000206.html[permanent dead link] (2 of 2) [2/6/2001 11:42:28 AM] Online: Accessed 27 November 2007 [1].
- ^ Nelson ME, Rinzel J (1994). "The Hodgkin–Huxley Model" (PDF). In Bower J, Beeman D (eds.). The Book of GENESIS: Exploring Realistic Neural Models with the GEneral NEural SImulation System. New York: Springer Verlag. pp. 29–49.
- PMID 994046.
- ^ ISBN 978-0-415-32868-5.
- ISBN 978-0-387-90819-9.
- .
- S2CID 6789007.
- PMID 19431309.
- S2CID 51648050.
- ^ FitzHugh R (1969). "Mathematical models of axcitation and propagation in nerve". In HP Schwann (ed.). Biological Engineering. New York: McGraw-Hill. pp. 1–85.
- S2CID 19149460., van der Mark J (1929). "The heartbeat considered as a relaxation oscillation, and an electrical model of the heart". Arch. Neerl. Physiol. 14: 418–443.
van der Pol B (1926). "On relaxation-oscillations". Philosophical Magazine. 2: 978–992.
van der Pol B, van der Mark J (1928). "The heartbeat considered as a relaxation oscillation, and an electrical model of the heart". Philosophical Magazine. 6: 763–775.
van der Pol B - S2CID 20077648.
- PMID 7260316.
- ISBN 978-0-262-11133-1.
- doi:10.1137/0141042.
- PMID 32355185.
- ^ LCCN 66015872.
- PMID 15410483.
- .
Woodbury JW (1965). "Chapter 3: Potentials in a volume conductor". In TC Ruch; HD Patton (eds.). Physiology and Biophysics. Philadelphia: W. B. Saunders Co.
Further reading
- Glass L, Mackey MC (1988). From Clocks to Chaos: The Rhythms of Life. Princeton, New Jersey: Princeton University. ISBN 978-0-691-08496-1.