Email updates

Keep up to date with the latest news and content from BMC Neuroscience and BioMed Central.

This article is part of the supplement: Seventeenth Annual Computational Neuroscience Meeting: CNS*2008

Open Access Poster presentation

A novel method for modelling nonlinear dynamical systems applied to the Hodgkin-Huxley neuron

Caitriona Boushel* and Paul Curran

Author Affiliations

School of Electronic, Electrical and Mechanical Engineering, University College Dublin, Ireland

For all author emails, please log on.

BMC Neuroscience 2008, 9(Suppl 1):P39  doi:10.1186/1471-2202-9-S1-P39


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2202/9/S1/P39


Published:11 July 2008

© 2008 Boushel and Curran; licensee BioMed Central Ltd.

Poster presentation

Two dynamical systems are considered equivalent if their bifurcation diagrams (or part of their bifurcation diagrams) are topologically equivalent. Hence, the two systems display the same bifurcation behaviour. In addition to this, we propose that the input-output behaviour of the systems and the behaviour of small networks of the systems will be qualitatively similar. To support this claim, a second order system, which is equivalent to the Hodgkin-Huxley(HH) neuron, is presented. The HH neuron [1] is the standard biologically plausible model of the electrical activity in a neuron. It consists of four coupled nonlinear ordinary differential equations that relate the potential across the neuron's membrane to three ionic currents flowing through the membrane and an externally applied current. The second order system, with u and μ as the bifurcation parameters,

<a onClick="popup('http://www.biomedcentral.com/1471-2202/9/S1/P39/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2202/9/S1/P39/mathml/M1">View MathML</a>

(1)

with r ≥ 0, φ ∈ [0, 2π] and σ0(μ) = -μ2, τ0(μ) = 1-μ, σ1(μ) = -2μ and <a onClick="popup('http://www.biomedcentral.com/1471-2202/9/S1/P39/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2202/9/S1/P39/mathml/M2">View MathML</a>

can approximate the two-parameter bifurcation diagram of the HH neuron, with the input current (I) and the potassium current equilibrium potentials (E>K) as the bifurcation parameters (see Fig. 1). It also reproduces much of the input-output behaviour of the HH neuron. When single input pulses are applied to the second order neuron and the HH neuron, both neurons switch from their stable equilibria to stable limit cycles and back again at the same time. When a continuous train of input pulses or a sinusoidal input are applied, both neurons oscillate at the frequency of the input signal once its amplitude is large enough (see Fig. 2). In addition, the oscillations of two second order neurons synchronise when the neurons are coupled, replicating the behaviour of two coupled HH neurons. These results indicate that nonlinear dynamical systems, which have topologically equivalent bifurcation diagrams, will have qualitatively similar input-output behaviour and qualitatively similar network behaviour as well as the same bifurcation behaviour.

thumbnailFigure 1. Two-parameter bifurcation diagrams of the HH neuron (solid lines) and the second order neuron (dotted lines) in the u-μ parameter. The supercritical Hopf, subcritical Hopf and double limit cycle bifurcation curves are labelled sH, uH and dc.

thumbnailFigure 2. Frequency spectra of the HH neuron (upper plot) and the second order neuron (lower plot). The dotted spectra are with no input to the neurons and the solid spectra are with a continuous train of input pulses.

References

  1. Hodgkin AL, Huxley AF: A quantitative description of membrane current and its application to conduction and excitation in nerve.

    J Physiol 1952, 117:500-544. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL