# Quantitative models of the action potential

In neurophysiology, several mathematical models of the action potential have been developed, which fall into two basic types. The first type seeks to model the experimental data quantitatively, i.e., to reproduce the measurements of current and voltage exactly. The renowned Hodgkin–Huxley model of the axon from the *Loligo* squid exemplifies such models.^{[1]} Although qualitatively correct, the H-H model does not describe every type of excitable membrane accurately, since it considers only two ions (sodium and potassium), each with only one type of voltage-sensitive channel. However, other ions such as calcium may be important and there is a great diversity of channels for all ions.^{[2]} As an example, the cardiac action potential illustrates how differently shaped action potentials can be generated on membranes with voltage-sensitive calcium channels and different types of sodium/potassium channels. The second type of mathematical model is a simplification of the first type; the goal is not to reproduce the experimental data, but to understand qualitatively the role of action potentials in neural circuits. For such a purpose, detailed physiological models may be unnecessarily complicated and may obscure the "forest for the trees". The Fitzhugh-Nagumo model is typical of this class, which is often studied for its entrainment behavior.^{[3]} Entrainment is commonly observed in nature, for example in the synchronized lighting of fireflies, which is coordinated by a burst of action potentials;^{[4]} entrainment can also be observed in individual neurons.^{[5]} Both types of models may be used to understand the behavior of small biological neural networks, such as the central pattern generators responsible for some automatic reflex actions.^{[6]} Such networks can generate a complex temporal pattern of action potentials that is used to coordinate muscular contractions, such as those involved in breathing or fast swimming to escape a predator.^{[7]}

## Contents

## Hodgkin–Huxley model

In 1952 Alan Lloyd Hodgkin and Andrew Huxley developed a set of equations to fit their experimental voltage-clamp data on the axonal membrane.^{[1]}^{[8]} The model assumes that the membrane capacitance *C* is constant; thus, the transmembrane voltage *V* changes with the total transmembrane current *I*_{tot} according to the equation

where *I*_{Na}, *I*_{K}, and *I*_{L} are currents conveyed through the local sodium channels, potassium channels, and "leakage" channels (a catch-all), respectively. The initial term *I*_{ext} 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 conductance is zero, whereas if open, its conductance is some constant value *g*. Hence, the net current through an ion channel depends on two variables: the probability *p*_{open} of the channel being open, and the difference in voltage from that ion's equilibrium voltage, *V* − *V*_{eq}. For example, the current through the potassium channel may be written as

which is equivalent to Ohm's law. By definition, no net current flows (*I*_{K} = 0) when the transmembrane voltage equals the equilibrium voltage of that ion (when *V* = *E*_{K}).

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 *p*_{n} 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., *p*_{open, K} = *n*^{4}. 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, *p*_{open, Na} = *m*^{3}*h*. The probabilities for each gate are assumed to obey first-order kinetics

where both the equilibrium value *m*_{eq} 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 *m*_{eq}; however, if *V* changes more quickly, then *m* will lag behind *m*_{eq}. 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 ordinary differential equations in four independent variables: the transmembrane voltage *V*, and the probabilities *m*, *h* and *n*.^{[10]} No general solution of these equations has been discovered. A less ambitious but generally applicable method for studying such non-linear dynamical systems is to consider their behavior in the vicinity of a fixed point.^{[11]} This analysis shows that the Hodgkin–Huxley system undergoes a transition from stable quiescence to bursting oscillations as the stimulating current *I*_{ext} is gradually increased; remarkably, the axon becomes stably quiescent again as the stimulating current is increased further still.^{[12]} A more general study of the types of qualitative behavior of axons predicted by the Hodgkin–Huxley equations has also been carried out.^{[10]}

## 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).^{[3]}^{[21]} All of these models are capable of exhibiting entrainment, which is commonly observed in nervous systems.^{[3]}

## 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.^{[22]} Such models are helpful in interpreting data from extracellular electrodes, which were common prior to the invention of the glass pipette electrode that allowed intracellular recording.^{[23]} The extracellular medium may be modeled as a normal isotropic ionic solution; in such solutions, the current follows the electric field lines, according to the continuum form of Ohm's Law

where **j** and **E** are vectors representing the current density and electric field, respectively, and where σ is the conductivity. Thus, **j** can be found from **E**, which in turn may be found using Maxwell's equations. Maxwell's equations can be reduced to a relatively simple problem of electrostatics, since the ionic concentrations change too slowly (compared to the speed of light) for magnetic effects to be important. The electric potential φ(**x**) at any extracellular point **x** can be solved using Green's identities^{[22]}

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.^{[24]}

## See also

## References

- ↑
^{1.0}^{1.1}^{1.2}Hodgkin AL, Huxley AF, Katz B (1952). "Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo".*Journal of Physiology*.**116**(4): 424–448. doi:10.1113/jphysiol.1952.sp004717. PMC 1392213. PMID 14946713.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Hodgkin AL, Huxley AF (1952). "Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo".*Journal of Physiology*.**116**(4): 449–472. doi:10.1113/jphysiol.1952.sp004717. PMC 1392213. PMID 14946713.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Hodgkin AL, Huxley AF (1952). "The components of membrane conductance in the giant axon of Loligo".*J Physiol*.**116**(4): 473–496. doi:10.1113/jphysiol.1952.sp004718. PMC 1392209. PMID 14946714.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Hodgkin AL, Huxley AF (1952). "The dual effect of membrane potential on sodium conductance in the giant axon of Loligo".*J Physiol*.**116**(4): 497–506. doi:10.1113/jphysiol.1952.sp004719. PMC 1392212. PMID 14946715.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Hodgkin AL, Huxley AF (1952). "A quantitative description of membrane current and its application to conduction and excitation in nerve".*J Physiol*.**117**(4): 500–544. doi:10.1113/jphysiol.1952.sp004764. PMC 1392413. PMID 12991237.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Yamada WM, Koch C, Adams PR (1989). "Multiple Channels and Calcium Dynamics". In C. Koch, I Segev (ed.).
*Methods in Neuronal Modeling: From Synapses to Networks*. Cambridge, Massachusetts: Bradford Book, The MIT Press. pp. 97–133. ISBN 0-262-11133-0.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑
^{3.0}^{3.1}^{3.2}^{3.3}Hoppensteadt FC (1986).*An introduction to the mathematics of neurons*. Cambridge: Cambridge University Press. ISBN 0-521-31574-3.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Hanson, F.E.; Case, J.F.; Buck, E.; Buck, J. (1971). "Synchrony and Flash Entrainment in a New Guinea Firefly".
*Science*.**174**(4005): 161–164. Bibcode:1971Sci...174..161H. doi:10.1126/science.174.4005.161. PMID 17742039. Retrieved 2008-05-05.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Guttman R, Feldman L, Jacobsson E (1980). "Frequency entrainment of squid axon".
*J. Membr. Biol*.**56**(1): 9–18. doi:10.1007/BF01869347. PMID 7441721.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Getting PA (1989). "Reconstruction of Small Neural Networks". In C Koch and I Segev (ed.).
*Methods in Neuronal Modeling: From Synapses to Networks*. Cambridge, Massachusetts: Bradford Book, The MIT Press. pp. 171–194. ISBN 0-262-11133-0.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Hooper, Scott L. "Central Pattern Generators."
*Embryonic ELS*(1999) http://www.els.net/elsonline/figpage/I0000206.html (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 (ed.).
*The Book of GENESIS: Exploring Realistic Neural Models with the GEneral NEural SImulation System*. New York: Springer Verlag. pp. 29–49.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Armstrong CM, Bezanilla F, Rojas E (1973). "Destruction of Sodium Conductance Inactivation in Squid Axons Perfused with Pronase".
*J. Gen. Physiol*.**48**(4): 375–391. doi:10.1085/jgp.62.4.375. PMC 2226121. PMID 4755846.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Rojas E, Rudy B (1976). "Destruction of the sodium conductance inactivation by a specific protease in perfused nerve fibres from*Loligo*".*J. Physiol*.**262**: 501–531.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑
^{10.0}^{10.1}Sato S, Fukai H, Nomura T, Doi S (2005). "Bifurcation Analysis of the Hodgkin–Huxley Equations". In Reeke GN, Poznanski RR, Lindsay KA, Rosenberg JR, Sporns O (ed.).*Modeling in the Neurosciences: From Biological Systems to Neuromimetic Robotics*(2nd ed.). Boca Raton: CRC Press. pp. 459–478. ISBN 978-0-415-32868-5.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Guckenheimer J, Holmes P (1986).
*Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields*(2nd printing, revised and corrected ed.). New York: Springer Verlag. pp. 12–16. ISBN 0-387-90819-6.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Sabah NH, Spangler RA (1970). "Repetitive response of the Hodgkin-Huxley model for the squid giant axon".
*Journal of Theoretical Biology*.**29**(2): 155–171. doi:10.1016/0022-5193(70)90017-2. PMID 5500466.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Evans JW (1972). "Nerve axon equations. I. Linear approximations".*Indiana U. Math. Journal*.**21**(9): 877–885. doi:10.1512/iumj.1972.21.21071.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Evans JW, Feroe J (1977). "Local stability theory of the nerve impulse".*Math. Biosci*.**37**: 23–50. doi:10.1016/0025-5564(77)90076-1.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ FitzHugh R (1960). "Thresholds and Plateaus in the Hodgkin-Huxley Nerve Equations".
*J. Gen. Physiol*.**43**(5): 867–896. doi:10.1085/jgp.43.5.867. PMC 2195039. PMID 13823315.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Kepler TB, Abbott LF, Marder E (1992). "Reduction of conductance-based neuron models".*Biological Cybernetics*.**66**(5): 381–387. doi:10.1007/BF00197717. PMID 1562643.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ FitzHugh R (1961). "Impulses and Physiological States in Theoretical Models of Nerve Membrane".
*Biophysical Journal*.**1**(6): 445–466. Bibcode:1961BpJ.....1..445F. doi:10.1016/S0006-3495(61)86902-6. PMC 1366333. PMID 19431309.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Nagumo J, Arimoto S, Yoshizawa S (1962). "An active pulse transmission line simulating nerve axon".
*Proceedings of the IRE*.**50**(10): 2061–2070. doi:10.1109/JRPROC.1962.288235.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ FitzHugh R (1969). "Mathematical models of axcitation and propagation in nerve". In HP Schwann (ed.).
*Biological Engineering*. New York: McGraw-Hill. pp. 1–85.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Bonhoeffer KF (1948). "ACTIVATION OF PASSIVE IRON AS A MODEL FOR THE EXCITATION OF NERVE".
*J. Gen. Physiol*.**32**(1): 69–91. doi:10.1085/jgp.32.1.69. PMC 2213747. PMID 18885679.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Bonhoeffer KF (1953). "Modelle der Nervenerregung".*Naturwissenschaften*.**40**(11): 301–311. Bibcode:1953NW.....40..301B. doi:10.1007/BF00632438.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

van der Pol B (1926). "On relaxation-oscillations".*Philosophical Magazine*.**2**: 978–992.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

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.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

van der Pol B, 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.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Keener JP (1983). "Analogue circuitry for the van der Pol and FitzHugh-Nagumo equations".
*IEEE Trans. on Systems, Man and Cybernetics*.**13**: 1010–1014. doi:10.1109/TSMC.1983.6313098.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Morris C, Lecar H (1981). "Voltage oscillations in the barnacle giant muscle fiber".
*Biophysical Journal*.**35**(1): 193–213. Bibcode:1981BpJ....35..193M. doi:10.1016/S0006-3495(81)84782-0. PMC 1327511. PMID 7260316.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Rinzel J, Ermentrout GB (1989). "Analysis of Neural Excitability and Oscillations". In C. Koch, I Segev (ed.).
*Methods in Neuronal Modeling: From Synapses to Networks*. Cambridge, Massachusetts: Bradford Book, The MIT Press. pp. 135–169. ISBN 0-262-11133-0.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Keener JP, Hoppensteadt FC, Rinzel J (1981). "Integrate-and-fire models of nerve membrane response to oscillatory input".
*SIAM Journal on Applied Mathematics*.**41**(3): 503–517. doi:10.1137/0141042.CS1 maint: multiple names: authors list (link)<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑
^{22.0}^{22.1}Stevens, Charles F. (1966).*Neurophysiology: A Primer*. New York: John Wiley and Sons. pp. 161–173. LCCN 66015872.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Ling G, Gerard RW (1949). "The normal membrane potential of frog sartorius fibers".
*J. Cell. Comp. Physiol*.**34**(3): 383–396. doi:10.1002/jcp.1030340304. PMID 15410483.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles> - ↑ Lorente de No R (1947). "A Study of Nerve Physiology".
*Stud. Rockefeller Inst. Med. Research*.**132**: Chap. 16.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Mauro A (1960). "Properties of thin generators pertaining to electrophysiological potentials in volume conductors".*J. Neurophysiol*.**23**: 132–143.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

Woodbury JW (1965). "Chapter 3: Potentials in a volume conductor". In TC Ruch, HD Patton (ed.).*Physiology and Biophysics*. Philadelphia: W. B. Saunders Co.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>

## Further reading

- Glass L, Mackey MC (1988).
*From Clocks to Chaos: The Rhythms of Life*. Princeton, New Jersey: Princeton University. ISBN 0-691-08496-3.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>