Periodic and non-periodic responses of a periodically forced Hodgkin-Huxley oscillator

Membrane potential responses of a Hodgkin-Huxley oscillator to an externally-applied sinusoidal current were numerically calculated with relation to bifurcation parameters of the amplitude and the frequency of the stimulating current. The Hodgkin-Huxley oscillator, or the Hodgkin-Huxley axon in the state of self-sustained oscillation of action potentials, was realized by immersing the axon in calcium-deficient sea water. The forced oscillations were analysed by the stroboscopic plots and/or the Lorenz plots. The results show that the periodically forced Hodgkin-Huxley oscillator exhibits not only periodic motions (harmonic or sub-harmonic synchronization) but also non-periodic motions (quasi-periodic or chaotic oscillation), that the motions were determined by the amplitude and the frequency of the stimulating current, and that the characteristic motions obtained in the present study were in reasonable agreement with those of our previous results, found experimentally in squid giant axons. Also, two kinds of routes to the chaotic oscillations were found; successive period-doubling bifurcations and formation of the intermittently chaotic oscillation from sub-harmonic synchronization.     

Dynamical mechanism for coexistence of dispersing species.

Dispersal of organisms may play an essential role in the coexistence of species. Recent studies of the evolution of dispersal in temporally varying environments suggest that clones differing in dispersal rates can coexist indefinitely. In this work, we explore the mechanism permitting such coexistence for a model of dispersal in a patchy environment, where temporal heterogeneity arises from endogenous chaotic dynamics. We show that coexistence arises from an extreme type of intermittent behavior, namely the phenomenon known as on-off intermittency. In effect, coexistence arises because of an alternation between synchronized and de-synchronized dynamical behaviors. Our analysis of the dynamical mechanism for on-off intermittency lends strong credence to the proposition that chaotic synchronism may be a general feature of species coexistence, where competing species differ only in dispersal rate.     

A fully coupled transient excited state model for the sodium channel

Summary The axon membrane is simulated by standard Hodgkin-Huxley leakage and potassium channels plus a coupled transient excited state kinetic scheme for the sodium channel. This scheme for the sodium channel is as proposed previously by the author. Simultations are presented showing the form of the action potential, threshold behavior, accommodation, and repetitive firing. It is seen that the form of the individual action potential, its all-or-none nature, and its refractory period are well simulated by this model, as they are by the standard Hodgkin-Huxley model. However, the model differs markedly from the Hodgkin-Huxley model with respect to repetitive firing and accommodation to stimulating currents of slowly rising intensity, in ways that are anomn to be related to those features of the sodium inactivation which are anomalous to the H-H model. The tendency for repetitive firing is highly dependent on that parameter which primarily determintes the existence of the inactivation shift in voltage clamp experiments, in such a way that the more pronounced the inactivation shift, the less the tendency for repetitive firing,. The tendency for accommodation is highly dependent on that parameter which primarily determines the “τ_c − τ_h” separation, in such a way that the greater the separation the greater the tendency for the membrane to accommodate without firing action potentials to a slowly rising current.     

Action Potential Initiation in the Hodgkin-Huxley Model

A recent paper of B. Naundorf et al. described an intriguing negative correlation between variability of the onset potential at which an action potential occurs (the onset span) and the rapidity of action potential initiation (the onset rapidity). This correlation was demonstrated in numerical simulations of the Hodgkin-Huxley model. Due to this antagonism, it is argued that Hodgkin-Huxley-type models are unable to explain action potential initiation observed in cortical neurons in vivo or in vitro. Here we apply a method from theoretical physics to derive an analytical characterization of this problem. We analytically compute the probability distribution of onset potentials and analytically derive the inverse relationship between onset span and onset rapidity. We find that the relationship between onset span and onset rapidity depends on the level of synaptic background activity. Hence we are able to elucidate the regions of parameter space for which the Hodgkin-Huxley model is able to accurately describe the behavior of this system.     

Effects of transient receptor potential-like current on the firing pattern of action potentials in the Hodgkin-Huxley neuron during exposure to sinusoidal external voltage

Transient receptor potential vanilloid-1 (TRPV1) channels play a role in several inflammatory and nociceptive processes. Previous work showed that magnetic electrical field-induced antinociceptive [corrected] action is mediated by activation of capsaicin-sensitive sensory afferents. In this study, a modified Hodgkin-Huxley model, in which TRP-like current (ITRP) was incorporated, was implemented to predict the firing behavior of action potentials (APs), as the model neuron was exposed to sinusoidal changes in externally-applied voltage. When model neuron is exposed to low-frequency sinusoidal voltage, increased maximal conductance of ITRP can enhance repetitive bursts of APs accompanied by a shortening of inter-spike interval (ISI) in AP firing. The change in ISIs with number of interval is periodic with the phase-locking. In addition, increased maximal conductance of ITRP can abolish chaotic pattern of AP firing in model neuron during exposure to high-frequency voltage. The ISI pattern is converted from irregular to constant, as maximal conductance of ITRP is increased under such high-frequency voltage. Our simulation results suggest that modulation of TRP-like channels functionally expressed in small-diameter peripheral sensory neurons should be an important mechanism through which it can contribute to the firing pattern of APs.     

Ghostbursting: A Novel Neuronal Burst Mechanism

Pyramidal cells in the electrosensory lateral line lobe (ELL) of weakly electric fish have been observed to produce high-frequency burst discharge with constant depolarizing current (Turner et al., 1994). We present a two-compartment model of an ELL pyramidal cell that produces burst discharges similar to those seen in experiments. The burst mechanism involves a slowly changing interaction between the somatic and dendritic action potentials. Burst termination occurs when the trajectory of the system is reinjected in phase space near the ghost of a saddle-node bifurcation of fixed points. The burst trajectory reinjection is studied using quasi-static bifurcation theory, that shows a period doubling transition in the fast subsystem as the cause of burst termination. As the applied depolarization is increased, the model exhibits first resting, then tonic firing, and finally chaotic bursting behavior, in contrast with many other burst models. The transition between tonic firing and burst firing is due to a saddle-node bifurcation of limit cycles. Analysis of this bifurcation shows that the route to chaos in these neurons is type I intermittency, and we present experimental analysis of ELL pyramidal cell burst trains that support this model prediction. By varying parameters in a way that changes the positions of both saddle-node bifurcations in parameter space, we produce a wide gallery of burst patterns, which span a significant range of burst time scales.     

Degenerate Hopf bifurcation and isolated periodic solutions of the Hodgkin-Huxley model with varying sodium ion concentration

Points of degenerate Hopf bifurcation in the Hodgkin-Huxley model are found as parameters temperature T and voltage level of sodium VNa are varied. Local techniques of degenerate Hopf bifurcation analysis are used to show the existence of families of periodic solutions of the model: isolated branches of periodic solutions (i.e. branches not connected to the stationary branch) are found in addition to Hopf branches. Purely numerical techniques are used to show that the isolas persist for VNa up to a value slightly greater than 114 mV. Under some conditions there are multiple stable periodic solutions, so "jumping" between action potentials of different amplitudes might be observed.     

Bifurcations in the decremental propagation of a spike train in the Hodgkin-Huxley model of low excitability

Response of a nerve fiber of low excitability to periodic stimulus pulses is studied with computer simulation of the Hodgkin-Huxley model. The excitability of the Hodgkin-Huxley model is reduced by decreasing the equilibrium potential for the sodium ion and by increasing the temperature, so that the decremental propagation of spikes occurs in the refractory period. It is shown that, as the period of stimulus pulses is decreased, the propagation length of the spikes is continuously changed, and period-doubling bifurcations occur. The response of a nerve fiber of low excitability is then qualitatively different from that of a normal fiber. Received: 6 December 1996 / Accepted in revised form: 12 June 1998     

Dissection and reduction of a modeled bursting neuron

An 11-variable Hodgkin-Huxley type model of a bursting neuron was investigated using numerical bifurcation analysis and computer simulations. The results were applied to develop a reduced model of the underlying subthreshold oscillations (slow-wave) in membrane potential. Two different low-order models were developed: one 3-variable model, which mimicked the slow-wave of the full model in the absence of action potentials and a second 4-variable model, which included expressions accounting for the perturbational effects of action potentials on the slow-wave. The 4-variable model predicted more accurately the activity mode (bursting, beating, or silence) in response to application of extrinsic stimulus current or modulatory agents. The 4-variable model also possessed a phase-response curve that was very similar to that of the original 11-variable model. The results suggest that low-order models of bursting cells that do not consider the effects of action potentials may erroneously predict modes of activity and transient responses of the full model on which the reductions are based. These results also show that it is possible to develop low-order models that retain many of the characteristics of the activity of the higher-order system.     

Computation of Impulse Initiation and Saltatory Conduction in a Myelinated Nerve Fiber

A mathematical model of the electrical properties of a myelinated nerve fiber is given, consisting of the Hodgkin-Huxley ordinary differential equations to represent the membrane at the nodes of Ranvier, and a partial differential cable equation to represent the internodes. Digital computer solutions of these equations show an impulse arising at a stimulating electrode and being propagated away, approaching a constant velocity. Action potential curves plotted against distance show discontinuities in slope, proportional to the nodal action currents, at the nodes. Action potential curves plotted against time, at the nodes and in the internodes, show a marked difference in steepness of the rising phase, but little difference in peak height. These results and computed action current curves agree fairly accurately with published experimental data from frog and toad fibers.     

Slow inward Ca current in frog heart: theoretical evidence against a voltage-dependent inactivation

The validity of a Hodgkin-Huxley type voltage-dependent inactivation of slow inward Ca current (Isi) was tested in frog heart using a computer simulation. The time course of Isi was calculated during the development of a frog atrial action potential (AP). With a time constant of inactivation (tauf) of 55 ms at a membrane potential (Em) of -15 mV, the variation of Isi was biphasic: after a transient increase followed by a decrease to zero, Isi partially "reactivated" (at the beginning of the AP repolarization phase) and then fully deactivated. The "reactivation" phase of Isi developed whether tauf was an increasing, decreasing, U-shaped, or bell-shaped function of Em. The addition of an independent and slower process responsible for the recovery from inactivation only partly suppressed the "reactivation" phase. However, until now there was no experimental evidence supporting such a biphasic variation of Isi during AP repolarization. Thus our results indicate that the Hodgkin-Huxley type model of the voltage-dependence of Isi-inactivation process may not correctly represent the actual behavior of frog cardiac muscle.     

A Statistical Model for In Vivo Neuronal Dynamics

Single neuron models have a long tradition in computational neuroscience. Detailed biophysical models such as the Hodgkin-Huxley model as well as simplified neuron models such as the class of integrate-and-fire models relate the input current to the membrane potential of the neuron. Those types of models have been extensively fitted to in vitro data where the input current is controlled. Those models are however of little use when it comes to characterize intracellular in vivo recordings since the input to the neuron is not known. Here we propose a novel single neuron model that characterizes the statistical properties of in vivo recordings. More specifically, we propose a stochastic process where the subthreshold membrane potential follows a Gaussian process and the spike emission intensity depends nonlinearly on the membrane potential as well as the spiking history. We first show that the model has a rich dynamical repertoire since it can capture arbitrary subthreshold autocovariance functions, firing-rate adaptations as well as arbitrary shapes of the action potential. We then show that this model can be efficiently fitted to data without overfitting. We finally show that this model can be used to characterize and therefore precisely compare various intracellular in vivo recordings from different animals and experimental conditions.     

Giant squid-hidden canard: the 3D geometry of the Hodgkin–Huxley model

This work is motivated by the observation of remarkably slow firing in the uncoupled Hodgkin-Huxley model, depending on parameters tau( h ), tau( n ) that scale the rates of change of the gating variables. After reducing the model to an appropriate nondimensionalized form featuring one fast and two slow variables, we use geometric singular perturbation theory to analyze the model's dynamics under systematic variation of the parameters tau( h ), tau( n ), and applied current I. As expected, we find that for fixed (tau( h ), tau( n )), the model undergoes a transition from excitable, with a stable resting equilibrium state, to oscillatory, featuring classical relaxation oscillations, as I increases. Interestingly, mixed-mode oscillations (MMO's), featuring slow action potential generation, arise for an intermediate range of I values, if tau( h ) or tau( n ) is sufficiently large. Our analysis explains in detail the geometric mechanisms underlying these results, which depend crucially on the presence of two slow variables, and allows for the quantitative estimation of transitional parameter values, in the singular limit. In particular, we show that the subthreshold oscillations in the observed MMO patterns arise through a generalized canard phenomenon. Finally, we discuss the relation of results obtained in the singular limit to the behavior observed away from, but near, this limit.     

Two and three dimensional reductions of the Hodgkin-Huxley system: separation of time scales and bifurcation schemes

We study two different two-dimensional reductions of the Hodgkin-Huxley equations. We show that they display the same qualitative bifurcation scheme as the original equations but overestimate the current range where periodic emission occurs. This is essentially due to the assumption that the evolution of the sodium activation variable m is instantaneous with respect to the dynamics of the variables h and n, an hypothesis that breaks down at high values of the injected current. To prove this point we compare the current-amplitude relation, the current-frequency relation, and the shapes of individual spikes for the two reduced models to the results obtained for the original Hodgkin-Huxley model and for a three dimensional model with instantaneous sodium activation. We show that a more satisfying agreement with the original Hodgkin-Huxley equations is obtained if we modify the evolution equation for the potential by incorporating the prominent features of the dynamics of m.     

Birhythmicity, trirhythmicity and chaos in bursting calcium oscillations

We have analyzed various types of complex calcium oscillations. The oscillations are explained with a model based on calcium-induced calcium release (CICR). In addition to the endoplasmic reticulum as the main intracellular Ca2+ store, mitochondrial and cytosolic Ca2+ binding proteins are also taken into account. This model was previously proposed for the study of the physiological role of mitochondria and the cytosolic proteins in gene rating complex Ca2+ oscillations [1]. Here, we investigated the occurrence of different types of Ca2+ oscillations obtained by the model, i.e. simple oscillations, bursting, and chaos. In a bifurcation diagram, we have shown that all these various modes of oscillatory behavior are obtained by a change of only one model parameter, which corresponds to the physiological variability of an agonist. Bursting oscillations were studied in more detail because they express birhythmicity, trirhythmicity and chaotic behavior. Two different routes to chaos are observed in the model: in addition to the usual period doubling cascade, we also show intermittency. For the characterization of the chaotic behavior, we made use of return maps and Lyapunov exponents. The potential biological role of chaos in intracellular signaling is discussed.     

The relation of Vmax to INa, GNa, and h infinity in a model of the cardiac Purkinje fiber.

The inward sodium current in cardiac muscle is difficult to study by voltage clamp methods, so various indirect experimental measures have been used to obtain insight into its characteristics. These methods depend on the relationship between maximal upstroke velocity of the action potential (Vmax) and the sodium current (INa), usually defined in terms of the Hodgkin-Huxley model. These relationships were explored using an adaptation of this model to cardiac Purkinje fibers. In general Vmax corresponded to INa, and it could be used to determine the relationship of membrane potential to GNa, and h infinity. The results, however, depended on the method of stimulation of the action potential, and an optimal stimulation method was determined. A commonly used experimental technique called "membrane responsiveness" was shown to distort seriously the properties of steady-state gating inactivation that is supposed to measure. Estimation of the changes in maximal sodium conductance, such as those produced by tetrodotoxin (TTX), would be accurately measured. Some experimental results have indicated a voltage-dependent effect of TTX. Characteristics of the measures of TTX effect under those conditions were illustrated. In summary, calculations with a model of the cardiac Purkinje fiber action potential provide insight into the accuracy of certain experimental methods using maximal upstroke velocity as a measure of INa, and cast doubt on other experimental methods, such as membrane responsiveness.     

An analysis of the mammalian ventricular action potential

Summary A Hodgkin-Huxley model for ventricular excitation is abstracted from electrophysiological data. A singular perturbation analysis of the 8-dimensional phase portrait of the model characterizes the role of calcium during the plateau phase of the ventricular action potential and demonstrates how the calcium refractory period prevents tetanization. Supported in part by the Undergraduate Research Opportunities Program, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA     

Stochastic differential equation model for cerebellar granule cell excitability

Neurons in the brain express intrinsic dynamic behavior which is known to be stochastic in nature. A crucial question in building models of neuronal excitability is how to be able to mimic the dynamic behavior of the biological counterpart accurately and how to perform simulations in the fastest possible way. The well-established Hodgkin-Huxley formalism has formed to a large extent the basis for building biophysically and anatomically detailed models of neurons. However, the deterministic Hodgkin-Huxley formalism does not take into account the stochastic behavior of voltage-dependent ion channels. Ion channel stochasticity is shown to be important in adjusting the transmembrane voltage dynamics at or close to the threshold of action potential firing, at the very least in small neurons. In order to achieve a better understanding of the dynamic behavior of a neuron, a new modeling and simulation approach based on stochastic differential equations and Brownian motion is developed. The basis of the work is a deterministic one-compartmental multi-conductance model of the cerebellar granule cell. This model includes six different types of voltage-dependent conductances described by Hodgkin-Huxley formalism and simple calcium dynamics. A new model for the granule cell is developed by incorporating stochasticity inherently present in the ion channel function into the gating variables of conductances. With the new stochastic model, the irregular electrophysiological activity of an in vitro granule cell is reproduced accurately, with the same parameter values for which the membrane potential of the original deterministic model exhibits regular behavior. The irregular electrophysiological activity includes experimentally observed random subthreshold oscillations, occasional spontaneous spikes, and clusters of action potentials. As a conclusion, the new stochastic differential equation model of the cerebellar granule cell excitability is found to expand the range of dynamics in comparison to the original deterministic model. Inclusion of stochastic elements in the operation of voltage-dependent conductances should thus be emphasized more in modeling the dynamic behavior of small neurons. Furthermore, the presented approach is valuable in providing faster computation times compared to the Markov chain type of modeling approaches and more sophisticated theoretical analysis tools compared to previously presented stochastic modeling approaches.