Phase resetting and bifurcation in the ventricular myocardium.

With the dynamic differential equations of Beeler, G. W., and H. Reuter (1977, J. Physiol. [Lond.]. 268:177-210), we have studied the oscillatory behavior of the ventricular muscle fiber stimulated by a depolarizing applied current I app. The dynamic solutions of BR equations revealed that as I app increases, a periodic repetitive spiking mode appears above the subthreshold I app, which transforms to a periodic spiking-bursting mode of oscillations, and finally to chaos near the suprathreshold I app (i.e., near the termination of the periodic state). Phase resetting and annihilation of repetitive firing in the ventricular myocardium were demonstrated by a brief current pulse of the proper magnitude applied at the proper phase. These phenomena were further examined by a bifurcation analysis. A bifurcation diagram constructed as a function of I app revealed the existence of a stable periodic solution for a certain range of current values. Two Hopf bifurcation points exist in the solution, one just above the lower periodic limit point and the other substantially below the upper periodic limit point. Between each periodic limit point and the Hopf bifurcation, the cell exhibited the coexistence of two different stable modes of operation; the oscillatory repetitive firing state and the time-independent steady state. As in the Hodgkin-Huxley case, there was a low amplitude unstable periodic state, which separates the domain of the stable periodic state from the stable steady state. Thus, in support of the dynamic perturbation methods, the bifurcation diagram of the BR equation predicts the region where instantaneous perturbations, such as brief current pulses, can send the stable repetitive rhythmic state into the stable steady state.     

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     

Ionic channel density of excitable membranes can act as a bifurcation parameter

As the maximal K⁺-conductance (or K⁺-channel density) of the Hodgkin-Huxley equations is reduced, the stable resting membrane potential bifurcates at a subcritical Hopf bifurcation into small amplitude unstable oscillations. These small amplitude solutions jump to large amplitude periodic solutions that correspond to a repetitive discharge of action potentials. Thus the specific channel density can act as a bifurcation parameter, and can control the excitability and autorhythmicity of excitable membranes.     

Understanding bursting oscillations as periodic slow passages through bifurcation and limit points

We consider a biochemical system consisting of two allosteric enzyme reactions coupled in series. The system has been modeled by Decroly and Goldbeter (J. Theor. Biol. 124, 219 (1987)) and is described by three coupled, first-order, nonlinear, differential equations. Bursting oscillations correspond to a succession of alternating active and silent phases. The active phase is characterized by rapid oscillations while the silent phase is a period of quiescence. We propose an asymptotic analysis of the differential equations which is based on the limit of large allosteric constants. This analysis allows us to construct a time-periodic bursting solution. This solution is jumping periodically between a slowly varying steady state and a slowly varying oscillatory state. Each jump follows a slow passage through a bifurcation or limit point which we analyze in detail. Of particular interest is the slow passage through a supercritical Hopf bifurcation. The transition is from an oscillatory solution to a steady state solution. We show that the transition is delayed considerably and characterize this delay by estimating the amplitude of the oscillations at the Hopf bifurcation point.     

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.     

Hopf bifurcations in multiple-parameter space of the Hodgkin-Huxley equations II. Singularity theoretic approach and highly degenerate bifurcations

In the Hodgkin-Huxley equations (HH), we have identified the parameter regions in which either two stable periodic solutions with different amplitudes and periods and an equilibrium point or two stable periodic solutions coexist. The global structure of bifurcations in the multiple-parameter space in the HH suggested that the bistabilities of the periodic solutions are associated with the degenerate Hopf bifurcation points by which several qualitatively different behaviors are organized. In this paper, we clarify this by analyzing the details of the degenerate Hopf bifurcations using the singularity theory approach which deals with local bifurcations near a highly degenerate fixed point. Received: 23 April 1999 / Accepted in revised form: 24 September 1999     

Fluctuation-driven rhythmogenesis in an excitatory neuronal network with slow adaptation

We study an excitatory all-to-all coupled network of N spiking neurons with synaptically filtered background noise and slow activity-dependent hyperpolarization currents. Such a system exhibits noise-induced burst oscillations over a range of values of the noise strength (variance) and level of cell excitability. Since both of these quantities depend on the rate of background synaptic inputs, we show how noise can provide a mechanism for increasing the robustness of rhythmic bursting and the range of burst frequencies. By exploiting a separation of time scales we also show how the system dynamics can be reduced to low-dimensional mean field equations in the limit N → ∞. Analysis of the bifurcation structure of the mean field equations provides insights into the dynamical mechanisms for initiating and terminating the bursts.     

Analysis of bursting in a thalamic neuron model

We extend a quantitative model for low-voltage, slow-wave excitability based on the T-type calcium current (Wang et al. 1991) by juxtaposing it with a Hodgkin-Huxley-like model for fast sodium spiking in the high voltage regime to account for the distinct firing modes of thalamic neurons. We employ bifurcation analysis to illustrate the stimulus-response behavior of the full model under both voltage regimes. The model neuron shows continuous sodium spiking when depolarized sufficiently from rest. Depending on the parameters of calcium current inactivation, there are two types of low-voltage responses to a hyperpolarizing current step: a single rebound low threshold spike (LTS) upon release of the step and periodic LTSs. Bursting is seen as sodium spikes ride the LTS crest. In both cases, we analyze the LTS burst response by projecting its trajectory into a fast/slow phase plane. We also use phase plane methods to show that a potassium A-current shifts the threshold for sodium spikes, reducing the number of fast sodium spikes in an LTS burst. It can also annihilate periodic bursting. We extend the previous work of Rose and Hindmarsh (1989a–c) for a thalamic neuron and propose a simpler model for thalamic activity. We consider burst modulation by using a neuromodulator-dependent potassium leakage conductance as a control parameter. These results correspond with experiments showing that the application of certain neurotransmitters can switch firing modes. Received: 18 July 1993/Accepted in revised form: 22 January 1994     

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.     

Exploring how extracellular electric field modulates neuron activity through dynamical analysis of a two-compartment neuron model

To investigate how extracellular electric field modulates neuron activity, a reduced two-compartment neuron model in the presence of electric field is introduced in this study. Depending on neuronal geometric and internal coupling parameters, the behaviors of the model have been studied extensively. The neuron model can exist in quiescent state or repetitive spiking state in response to electric field stimulus. Negative electric field mainly acts as inhibitory stimulus to the neuron, positive weak electric field could modulate spiking frequency and spike timing when the neuron is already active, and positive electric fields with sufficient intensity could directly trigger neuronal spiking in the absence of other stimulations. By bifurcation analysis, it is observed that there is saddle-node on invariant circle bifurcation, supercritical Hopf bifurcation and subcritical Hopf bifurcation appearing in the obtained two parameter bifurcation diagrams. The bifurcation structures and electric field thresholds for triggering neuron firing are determined by neuronal geometric and coupling parameters. The model predicts that the neurons with a nonsymmetric morphology between soma and dendrite, are more sensitive to electric field stimulus than those with the spherical structure. These findings suggest that neuronal geometric features play a crucial role in electric field effects on the polarization of neuronal compartments. Moreover, by determining the electric field threshold of our biophysical model, we could accurately distinguish between suprathreshold and subthreshold electric fields. Our study highlights the effects of extracellular electric field on neuronal activity from the biophysical modeling point of view. These insights into the dynamical mechanism of electric field may contribute to the investigation and development of electromagnetic therapies, and the model in our study could be further extended to a neuronal network in which the effects of electric fields on network activity may be investigated.     

Functional significance of the A-current

This work considers the response to simulated synaptic inputs of an excitable membrane model. The model is essentially of the Hodgkin-Huxley type, but contains an A-current in addition to sodium and delayed-rectifier potassium channels. The results were compared with previous simulations in which the stimulus was an injected current. These two types of stimuli give somewhat different results because synaptic stimuli directly change the membrane resistance, whereas injected current does not. The results of synaptic stimulation were similar to injected current in that very low frequencies of action potentials were elicited only where the stimulus was slightly above threshold. For most of the range of synaptic inputs that produced oscillatory behavior, the A-current had little effect on oscillation frequency. With synaptic stimuli as with injected current, the model membrane's spiking behavior does not begin immediately when an excitatory stimulus is imposed on a quiescent state. The delay before spiking is closely related to the inactivation time of the A-current. The synaptic results were different from the injected current results in that when substantial inhibition was present, the ability to produce very-low-frequency spiking was absent, even just above the excitatory threshold. The higher the degree of inhibition, the narrower the range of spike frequencies that could be elicited by excitation. At very high inhibition, no degree of excitation could elicit spiking.     

The potassium A-current,low firing rates and rebound excitation in Hodgkin-Huxley models

It is widely believed, following the work of Connor and Stevens (1971,J. Physiol. Lond. 214, 31–53) that the ability to fire action potentials over a wide frequency range, especially down to very low rates, is due to the transient, potassium A-current (I _A). Using a reduction of the classical Hodgkin-Huxley model, we study the effects ofI _A on steady firing rate, especially in the near-threshold regime for the onset of firing. A minimum firing rate of zero corresponds to a homoclinic bifurcation of periodic solutions at a critical level of stimulating current. It requires that the membrane's steady-state current-voltage relation be N-shaped rather than monotonic. For experimentally based genericI _A parameters, the model does not fire at arbitrarily low rates, although it can for the more atypicalI _A parameters given by Connor and Stevens for the crab axon. When theI _A inactivation rate is slow, we find that the transient potassium current can mediate more complex firing patterns, such as periodic bursting in some parameter regimes. The number of spikes per burst increases asg _A decreases and as inactivation rate decreases. We also study howI _A affects properties of transient voltage responses, such as threshold and firing latency for anodal break excitation. We provide mathematical explanations for several of these dynamic behaviors using bifurcation theory and averaging methods.     

Reduction of conductance-based neuron models

We present a scheme for systematically reducing the number of differential equations required for biophysically realistic neuron models. The techniques are general, are designed to be applicable to a large set of such models and retain in the reduced system as high a degree of fidelity to the original system as possible. As examples, we provide reductions of the Hodgkin-Huxley system and the A-current model of Connor et al. (1977).     

Slow Passage Through a Hopf Bifurcation in Excitable Nerve Cables: Spatial Delays and Spatial Memory Effects

It is well established that in problems featuring slow passage through a Hopf bifurcation (dynamic Hopf bifurcation) the transition to large-amplitude oscillations may not occur until the slowly changing parameter considerably exceeds the value predicted from the static Hopf bifurcation analysis (temporal delay effect), with the length of the delay depending upon the initial value of the slowly changing parameter (temporal memory effect). In this paper we introduce new delay and memory effect phenomena using both analytic (WKB method) and numerical methods. We present a reaction–diffusion system for which slowly ramping a stimulus parameter (injected current) through a Hopf bifurcation elicits large-amplitude oscillations confined to a location a significant distance from the injection site (spatial delay effect). Furthermore, if the initial current value changes, this location may change (spatial memory effect). Our reaction–diffusion system is Baer and Rinzel’s continuum model of a spiny dendritic cable; this system consists of a passive dendritic cable weakly coupled to excitable dendritic spines. We compare results for this system with those for nerve cable models in which there is stronger coupling between the reactive and diffusive portions of the system. Finally, we show mathematically that Hodgkin and Huxley were correct in their assertion that for a sufficiently slow current ramp and a sufficiently large cable length, no value of injected current would cause their model of an excitable cable to fire; we call this phenomenon “complete accommodation.”     

Stability of periodic solutions for a model of genetic repression with delays

A technique is discussed for locating the Hopf bifurcation of an n-dimensional system of delay differential equations which arises from a model for control of protein biosynthesis. Certain parameter values are shown to allow a Hopf bifurcation to periodic orbits. At the Hopf bifurcation the periodic orbits are shown to be stable either analytically or numerically depending on the parameter values.On leave from North Carolina State University.Supported in part by N.S.F. Grant # MCS 81-02828     

Dynamics from Seconds to Hours in Hodgkin-Huxley Model with Time-Dependent Ion Concentrations and Buffer Reservoirs

The classical Hodgkin-Huxley (HH) model neglects the time-dependence of ion concentrations in spiking dynamics. The dynamics is therefore limited to a time scale of milliseconds, which is determined by the membrane capacitance multiplied by the resistance of the ion channels, and by the gating time constants. We study slow dynamics in an extended HH framework that includes time-dependent ion concentrations, pumps, and buffers. Fluxes across the neuronal membrane change intra- and extracellular ion concentrations, whereby the latter can also change through contact to reservoirs in the surroundings. Ion gain and loss of the system is identified as a bifurcation parameter whose essential importance was not realized in earlier studies. Our systematic study of the bifurcation structure and thus the phase space structure helps to understand activation and inhibition of a new excitability in ion homeostasis which emerges in such extended models. Also modulatory mechanisms that regulate the spiking rate can be explained by bifurcations. The dynamics on three distinct slow times scales is determined by the cell volume-to-surface-area ratio and the membrane permeability (seconds), the buffer time constants (tens of seconds), and the slower backward buffering (minutes to hours). The modulatory dynamics and the newly emerging excitable dynamics corresponds to pathological conditions observed in epileptiform burst activity, and spreading depression in migraine aura and stroke, respectively.     

Qualitative study of a dynamical system for metrazol-induced paroxysmal depolarization shifts

The convulsant pentylenetetrazol (PTZ) was used to trigger spike bursts and paroxysmal discharges inAplysia neurons. Voltage clamp experiments showed that PTZ induced a slow voltage-dependent potassium current and a persistent inward current. These currents are incorporated into a membrane model together with modified spike-generating Hodgkin-Huxley equations. From these data a metaphoric model is constructed and represented by a slow-fast dynamical system defined inR ⁴. With some values of the main physiological parameters, the system might have limit cycles for the fast dynamic. A qualitative study of the system shows that it satisfactorily reproduces the various observed patterns produced by PTZ.     

Bifurcations in a white-blood-cell production model

We study the dynamics of a model of white-blood-cell (WBC) production. The model consists of two compartmental differential equations with two discrete delays. We show that from normal to pathological parameter values, the system undergoes supercritical Hopf bifurcations and saddle-node bifurcations of limit cycles. We characterize the steady states of the system and perform a bifurcation analysis. Our results indicate that an increase in apoptosis rate of either hematopoietic stem cells or WBC precursors induces a Hopf bifurcation and an oscillatory regime takes place. These oscillations are seen in some hematological diseases.     

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.     

Critical Slowing Down Governs the Transition to Neuron Spiking

Many complex systems have been found to exhibit critical transitions, or so-called tipping points, which are sudden changes to a qualitatively different system state. These changes can profoundly impact the functioning of a system ranging from controlled state switching to a catastrophic break-down; signals that predict critical transitions are therefore highly desirable. To this end, research efforts have focused on utilizing qualitative changes in markers related to a system’s tendency to recover more slowly from a perturbation the closer it gets to the transition—a phenomenon called critical slowing down. The recently studied scaling of critical slowing down offers a refined path to understand critical transitions: to identify the transition mechanism and improve transition prediction using scaling laws.Here, we outline and apply this strategy for the first time in a real-world system by studying the transition to spiking in neurons of the mammalian cortex. The dynamical system approach has identified two robust mechanisms for the transition from subthreshold activity to spiking, saddle-node and Hopf bifurcation. Although theory provides precise predictions on signatures of critical slowing down near the bifurcation to spiking, quantitative experimental evidence has been lacking. Using whole-cell patch-clamp recordings from pyramidal neurons and fast-spiking interneurons, we show that 1) the transition to spiking dynamically corresponds to a critical transition exhibiting slowing down, 2) the scaling laws suggest a saddle-node bifurcation governing slowing down, and 3) these precise scaling laws can be used to predict the bifurcation point from a limited window of observation. To our knowledge this is the first report of scaling laws of critical slowing down in an experiment. They present a missing link for a broad class of neuroscience modeling and suggest improved estimation of tipping points by incorporating scaling laws of critical slowing down as a strategy applicable to other complex systems.