Coincidence detection in the Hodgkin-Huxley equations

Some of the cochlear nuclei in the auditory pathway are specialized for the sound localization. They compute the interaural time difference. The difference in sound timing is transduced by the dedicated neuronal circuit into a labeled line difference. The detector neurons along the delay line fire only when synaptic inputs reflecting signals from both cars arrive within a short time window. It was therefore called coincidence detection. We show, (1) what are the limits of coincidence detection in the leaky integrator model, which is a linear system, (2) how should the ideal coincidence detector based on the Hodkin-Huxley equations from real neurons look like, (3) what are the properties and physical limits in the real coincidence detection system. The conclusion is that the neuron with the Hodgkin Huxley dynamics has a fixed precision for the coincidence detection. The limits of the sound localization precision are set by the frequency of the sound and, therefore, by the vector strength of spike trains generated in the neuronal circuit in response to the sound.     

Isolated periodic solutions of the Hodgkin-Huxley equations

Periodic solutions of the current clamped Hodgkin-Huxley equations (Hodgkin & Huxley, 1952 J. Physiol. 117, 500) that arise by degenerate Hopf bifurcation were studied recently by Labouriau (1985 SIAM J. Math. Anal. 16, 1121, 1987 Degenerate Hopf Bifurcation and Nerve Impulse (Part II), in press). Two parameters, temperature T and sodium conductance gNa were varied from the original values obtained by Hodgkin & Huxley. Labouriau's work proved the existence of small amplitude periodic solution branches that do not connect locally to the stationary solution branch, and had not been previously computed. In this paper we compute these solution branches globally. We find families of isolas of periodic solutions (i.e. branches not connected to the stationary branch). For values of gNa in the range measured by Hodgkin & Huxley, and for physically reasonable temperatures, there are isolas containing orbitally asymptotically stable solutions. The presence of isolas of periodic solutions suggests that in certain current space clamped membrane experiments, action potentials could be observed even though the stationary state is stable for all current stimuli. Once produced, such action potentials will disappear suddenly if the current stimulus is either increased or decreased past certain values. Under some conditions, "jumping" between action potentials of different amplitudes might be observed.     

Canards for a reduction of the Hodgkin-Huxley equations

This paper shows that canards, which are periodic orbits for which the trajectory follows both the attracting and repelling part of a slow manifold, can exist for a two-dimensional reduction of the Hodgkin-Huxley equations. Such canards are associated with a dramatic change in the properties of the periodic orbit within a very narrow interval of a control parameter. By smoothly connecting stable and unstable manifolds in an asymptotic limit, we predict with great accuracy the parameter value at which the canards exist for this system. This illustrates the power of using singular perturbation theory to understand the dynamical properties of realistic biological systems.     

Modeling of time disparity detection by the Hodgkin-Huxley equations

Phase-sensitive neurons in the electrosensory lateral line lobe in the electrosensory pathway of the wave-type electric fish, Gymnarchus niloticus, are specialized for sensing the time disparity between sensory inputs at different parts of the body surface that is necessary for an electrical behavior, jamming avoidance response. These neurons are sensitive to time disparity in the microsecond range between synaptic inputs that represent occurrence times of electrosensory signals at different areas on the body surface. We showed that an ideal Hodgkin-Huxley equation may serve as a time disparity detector that fits physiological precision, and the precision for the time disparity detection is largely regulated by the maximal g(K) conductance in the Hodgkin-Huxley equations.     

Spontaneous synchronization of coupled circadian oscillators

In mammals, the circadian pacemaker, which controls daily rhythms, is located in the suprachiasmatic nucleus (SCN). Circadian oscillations are generated in individual SCN neurons by a molecular regulatory network. Cells oscillate with periods ranging from 20 to 28 h, but at the tissue level, SCN neurons display significant synchrony, suggesting a robust intercellular coupling in which neurotransmitters are assumed to play a crucial role. We present a dynamical model for the coupling of a population of circadian oscillators in the SCN. The cellular oscillator, a three-variable model, describes the core negative feedback loop of the circadian clock. The coupling mechanism is incorporated through the global level of neurotransmitter concentration. Global coupling is efficient to synchronize a population of 10,000 cells. Synchronized cells can be entrained by a 24-h light-dark cycle. Simulations of the interaction between two populations representing two regions of the SCN show that the driven population can be phase-leading. Experimentally testable predictions are: 1), phases of individual cells are governed by their intrinsic periods; and 2), efficient synchronization is achieved when the average neurotransmitter concentration would dampen individual oscillators. However, due to the global neurotransmitter oscillation, cells are effectively synchronized.     

Complex nonlinear dynamics of the Hodgkin-Huxley equations induced by time scale changes

 The Hodgkin–Huxley equations with a slight modification are investigated, in which the inactivation process (h) of sodium channels or the activation process of potassium channels (n) is slowed down. We show that the equations produce a variety of action potential waveforms ranging from a plateau potential, such as in heart muscle cells, to chaotic bursting firings. When h is slowed down – differently from the case of n variable being slow – chaotic bursting oscillations are observed for a wide range of parameter values although both variables cause a decrease in the membrane potential. The underlying nonlinear dynamics of various action potentials are analyzed using bifurcation theory and a so-called slow–fast decomposition analysis. It is shown that a simple topological property of the equilibrium curves of slow and fast subsystems is essential to the production of chaotic oscillations, and this is the cause of the large difference in global firing characteristics between the h-slow and n-slow cases. Received: 9 August 2000 / Accepted in revised form: 10 January 2001     

Noise controlled synchronization in potassium coupled neural models

The paper applies biologically plausible models to investigate how noise input to small ensembles of neurons, coupled via the extracellular potassium concentration, can influence their firing patterns. Using the noise intensity and the volume of the extracellular space as control parameters, we show that potassium induced depolarization underlies the formation of noise-induced patterns such as delayed firing and synchronization. These phenomena are associated with the appearance of new time scales in the distribution of interspike intervals that may be significant for the spatio-temporal oscillations in neuronal ensembles.     

Homomorphism on a physical system of the Hodgkin-Huxley equations for ion conductance in nerve

A homomorphism on a physical system of the Hodgkin-Huxley equations for ion conductance in nerve is derived. It is pointed out that a homomorphism can correct the Cole-Moore discrepancy in delay of conductance for voltage clamp data with initial hyperpolarization. The voltage dependence of the rate constants can also be removed. Curves are presented to compare the representation of the nerve conductances by the Hodgkin-Huxley equations and the new homomorphism.     

Pinning synchronization of coupled inertial delayed neural networks

The paper is devoted to the investigation of synchronization for an array of linearly and diffusively coupled inertial delayed neural networks (DNNs). By placing feedback control on a small fraction of network nodes, the entire coupled DNNs can be synchronized to a common objective trajectory asymptotically. Two different analysis methods, including matrix measure strategy and Lyapunov–Krasovskii function approach, are employed to provide sufficient criteria for the synchronization control problem. Comparisons of these two techniques are given at the end of the paper. Finally, an illustrative example is provided to show the effectiveness of the obtained theoretical results.     

A note on the asymptotic reduction of the Hodgkin-Huxley equations for nerve impulses

The (standard) FitzHugh reduction of the Hodgkin-Huxley equations for the propagation of nerve impulses ignores the dynamics of the activation gates. This assumption is invalid and leads to an over-estimation of the wave speed by a factor of 5 and the wrong dependence of wave speed on sodium channel conductance. The error occurs because a non-dimensional parameter, which is assumed to be small in the FitzHugh reduction, is in fact large (≈18). We analyse the Hodgkin-Huxley equations for propagating nerve impulses in the limit that this non-dimensional parameter is large, and show that the analytical results are consistent with numerical simulations of the Hodgkin-Huxley equations.     

Finite-time synchronization of coupled neural networks via discontinuous controllers

This paper investigates finite-time synchronization of an array of coupled neural networks via discontinuous controllers. Based on Lyapunov function method and the discontinuous version of finite-time stability theory, some sufficient criteria for finite-time synchronization are obtained. Furthermore, we propose switched control and adaptive tuning parameter strategies in order to reduce the settling time. In addition, pinning control scheme via a single controller is also studied in this paper. With the hypothesis that the coupling network topology contains a directed spanning tree and each of the strongly connected components is detail-balanced, we prove that finite-time synchronization can be achieved via pinning control. Finally, some illustrative examples are given to show the validity of the theoretical results.     

Intercellular synchronization of diffusively coupled Ca2+ oscillators

We examine the synchrony in the dynamics of localized [Ca^2 +]_i oscillations among a group of cells exhibiting such complex Ca^2 + oscillations, connected in the form of long chain, via diffusing coupling where cytosolic Ca^2 + and inositol 1,4,5-triphosphate are coupling molecules. Based on our numerical results, we could able to identify three regimes, namely desynchronized, transition and synchronized regimes in the (T − k_e) (time period-coupling constant) and (A − k_e) (amplitude-coupling constant) spaces which are supported by phase plots (Δϕ verses time) and recurrence plots, respectively. We further show the increase of synchronization among the cells as the number of coupling molecules increases in the (T − k_e) and (A − k_e) spaces.     

Hopf bifurcations in multiple-parameter space of the Hodgkin-Huxley equations I. Global organization of bistable periodic solutions

 The Hodgkin-Huxley equations (HH) are parameterized by a number of parameters and shows a variety of qualitatively different behaviors depending on the parameter values. We explored the dynamics of the HH for a wide range of parameter values in the multiple-parameter space, that is, we examined the global structure of bifurcations of the HH. Results are summarized in various two-parameter bifurcation diagrams with I _ext (externally applied DC current) as the abscissa and one of the other parameters as the ordinate. In each diagram, the parameter plane was divided into several regions according to the qualitative behavior of the equations. In particular, we focused on periodic solutions emerging via Hopf bifurcations and identified parameter regions in which either two stable periodic solutions with different amplitudes and periods and a stable equilibrium point or two stable periodic solutions coexist. Global analysis of the bifurcation structure suggested that generation of these regions is associated with degenerate Hopf bifurcations. Received: 23 April 1999 / Accepted in revised form: 24 September 1999     

Application of the Hodgkin-Huxley equations to an electric field model for interaction between excitable cells

We previously described a model for the electrical transfer of excitation from one cell to the next which utilized the electric potential generated in the junctional cleft between the cells. Low-resistance connections between the cells were not used in the model, and it was assumed that the junctional membranes were excitable. This model was analyzed for the static case without capacitances and for the dynamic case in which capacitances were part of the circuit elements. For simplicity, the Na⁺ resistance (R_Na), after a threshold potential was exceeded, was allowed to decrease exponentially (to 1% of its initial value) within 0·25–1·0 ms, and possible changes in the K⁺ resistance were ignored. In this paper, we have incorporated the Hodgkin-Huxley equations into the operation of the lumped membrane units for the electrical equivalent circuit of the cell membrane. The parameters varied are the membrane capacitances, resistances, maximum Na⁺ conductance ($\text{g}{}_{\mathrm{<mtext>Na</mtext>}}$), and the radial cleft resistance (R_jc). We demonstrated that our model worked very well, i.e. the successful transfer of action potentials was achieved, with the membrane units following Hodgkin-Huxley dynamics for changes in g_Na and g_K. The calculations indicate that transmission is facilitated when the junctional units have a higher $\text{g}{}_{\mathrm{<mtext>Na</mtext>}}$ and a lower capacitance and when R_jc is elevated. Lowering the resistance of the junctional membrane units several fold, relative to the surface membrane units, also facilitated transmission; however, the absolute resistance of the junctional membrane was still well above the maximum value that would allow sufficient local-circuit current to flow to effect transmission. Thus, the electric field model provides an alternative means of cell-to-cell propagation between myocardial cells which is electrical in nature but does not require the presence of low-resistance connections between cells.     

Mathematical description of a bursting pacemaker neuron by a modification of the Hodgkin-Huxley equations.

Modifications based on experimental results reported in the literature are made to the Hodgkin-Huxley equations to describe the electrophysiological behavior of the Aplysia abdominal ganglion R15 cell. The system is then further modified to describe the effects with the application of the drug tetrodotoxin (TTX) to the cells' bathing medium. Methods of the qualitative theory of differential equations are used to determine the conditions necessary for such a system of equations to have an oscillatory solution. A model satisfying these conditions is shown to preduct many experimental observations of R15 cell behavior. Numerical solutions are obtained for differential equations satisfying the conditions of the model. These solutions are shown to have a form similar to that of the bursting which is characteristic of this cell, and to preduct many results of experiments conducted on this cell. The physiological implications of the model are discussed.     

Power-rate synchronization of coupled genetic oscillators with unbounded time-varying delay

In this paper, a new synchronization problem for the collective dynamics among genetic oscillators with unbounded time-varying delay is investigated. The dynamical system under consideration consists of an array of linearly coupled identical genetic oscillators with each oscillators having unbounded time-delays. A new concept called power-rate synchronization, which is different from both the asymptotical synchronization and the exponential synchronization, is put forward to facilitate handling the unbounded time-varying delays. By using a combination of the Lyapunov functional method, matrix inequality techniques and properties of Kronecker product, we derive several sufficient conditions that ensure the coupled genetic oscillators to be power-rate synchronized. The criteria obtained in this paper are in the form of matrix inequalities. Illustrative example is presented to show the effectiveness of the obtained results.