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     

A two-variable model of somatic-dendritic interactions in a bursting neuron

We present a two-variable delay-differential-equation model of a pyramidal cell from the electrosensory lateral line lobe of a weakly electric fish that is capable of burst discharge. It is a simplification of a six-dimensional ordinary differential equation model for such a cell whose bifurcation structure has been analyzed (Doiron et al., J. Comput. Neurosci., 12, 2002). We have modeled the effects of back-propagating action potentials by a delay, and use an integrate-and-fire mechanism for action potential generation. The simplicity of the model presented here allows one to explicitly derive a two-dimensional map for successive interspike intervals, and to analytically investigate the effects of time-dependent forcing on such a model neuron. Some of the effects discussed include ‘burst excitability’, the creation of resonance tongues under periodic forcing, and stochastic resonance. We also investigate the effects of changing the parameters of the model.     

Multiple timescale mixed bursting dynamics in a respiratory neuron model

Experimental results in rodent medullary slices containing the pre-Bötzinger complex (pre-BötC) have identified multiple bursting mechanisms based on persistent sodium current (I _NaP) and intracellular Ca²⁺. The classic two-timescale approach to the analysis of pre-BötC bursting treats the inactivation of I _NaP, the calcium concentration, as well as the Ca²⁺-dependent inactivation of IP ₃ as slow variables and considers other evolving quantities as fast variables. Based on its time course, however, it appears that a novel mixed bursting (MB) solution, observed both in recordings and in model pre-BötC neurons, involves at least three timescales. In this work, we consider a single-compartment model of a pre-BötC inspiratory neuron that can exhibit both I _NaP and Ca²⁺ oscillations and has the ability to produce MB solutions. We use methods of dynamical systems theory, such as phase plane analysis, fast-slow decomposition, and bifurcation analysis, to better understand the mechanisms underlying the MB solution pattern. Rather surprisingly, we discover that a third timescale is not actually required to generate mixed bursting solutions. Through our analysis of timescales, we also elucidate how the pre-BötC neuron model can be tuned to improve the robustness of the MB solution.     

The role of a transient potassium current in a bursting neuron model

Presented here is a biophysical cell model which can exhibit low-frequency repetitive activity and bursting behavior. The model is developed from previous models (Av-Ron et al. 1991, 1993) for excitability, oscillations and bursting. A stepwise development of the present model shows the contribution of a transient potassium current (I _A ) to the overall dynamics. By changing a limited set of model parameters one can describe different firing patterns; oscillations with frequencies ranging from 2–200 Hz and a wide range of bursting behaviors in terms of the durations of bursting and quiescence, peak firing frequency and rate of change of the firing frequency.     

Variability of bursting patterns in a neuron model in the presence of noise

Spiking and bursting patterns of neurons are characterized by a high degree of variability. A single neuron can demonstrate endogenously various bursting patterns, changing in response to external disturbances due to synapses, or to intrinsic factors such as channel noise. We argue that in a model of the leech heart interneuron existing variations of bursting patterns are significantly enhanced by a small noise. In the absence of noise this model shows periodic bursting with fixed numbers of interspikes for most parameter values. As the parameter of activation kinetics of a slow potassium current is shifted to more hyperpolarized values of the membrane potential, the model undergoes a sequence of incremental spike adding transitions accumulating towards a periodic tonic spiking activity. Within a narrow parameter window around every spike adding transition, spike alteration of bursting is deterministically chaotic due to homoclinic bifurcations of a saddle periodic orbit. We have found that near these transitions the interneuron model becomes extremely sensitive to small random perturbations that cause a wide expansion and overlapping of the chaotic windows. The chaotic behavior is characterized by positive values of the largest Lyapunov exponent, and of the Shannon entropy of probability distribution of spike numbers per burst. The windows of chaotic dynamics resemble the Arnold tongues being plotted in the parameter plane, where the noise intensity serves as a second control parameter. We determine the critical noise intensities above which the interneuron model generates only irregular bursting within the overlapped windows.     

Dynamical analysis of periodic bursting in piece-wise linear planar neuron model

A piece-wise linear planar neuron model, namely, two-dimensional McKean model with periodic drive is investigated in this paper. Periodical bursting phenomenon can be observed in the numerical simulations. By assuming the formal solutions associated with different intervals of this non-autonomous system and introducing the generalized Jacobian matrix at the non-smooth boundaries, the bifurcation mechanism for the bursting solution induced by the slowly varying periodic drive is presented. It is shown that, the discontinuous Hopf bifurcation occurring at the non-smooth boundaries, i.e., the bifurcation taking place at the thresholds of the stimulation, leads the alternation between the rest state and spiking state. That is, different oscillation modes of this non-autonomous system convert periodically due to the non-smoothness of the vector field and the slow variation of the periodic drive as well.     

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.     

A study of the connection between the interneuron initiating pacemaker activity in a bursting neuron and the bursting neuron of the snailHelix pomatia

The connection between an interneuron initiating pacemaker activity in the bursting RPa1 neuron and the bursting neuron itself (Pin and Gola, 1983) has been analyzed in the snail Helix pomatia. Prolonged depolarization of the interneuronal membrane produced in it a series of action potentials as well as a parallel initiation or enhancement of bursting activity in the RPa1 neuron. If the discharge in the interneuron was evoked by short current pulses of threshold amplitude, no bursting activity was seen in the RPa1 neuron. However, short stimuli delivered on the background of subthreshold depolarization of the interneuronal membrane produced bursting activity in the RPa1 neuron. Under voltage-clamp conditions a slow inward current could be recorded in the RPa1 neuronal membrane after stimulation of the interneuron with a latency of about 2 sec. Short shifts of the holding potential in the hyperpolarizing direction at the maximum of this current produced a transient outward current. Replacement of extracellular Ca2+ by Mg2+ ions, as well as addition of 1 mM CdCl2 to the external solution, prevented the response to the interneuronal stimulation in the RPa1 neuron. Electron microscopic investigation of the interneuron has shown the abundance of Golgi complexes in its cytoplasm with electron-dense granules in their vicinity. It is concluded that the connection between the interneuron and the bursting neuron is of chemical origin, based on secretion by the former of some substances which activate at least two types of ionic channels in the membrane of the RPa1 neuron.     

Method for constructing the boundary of the bursting oscillations region in the neuron model

We examine the problem of constructing the boundary of bursting oscillations on a parameter plane for the system of equations describing the electrical behaviour of the membrane neuron arising from the interaction of fast oscillations of the cytoplasma membrane potential and slow oscillations of the intracellular calcium concentration. As the boundary point on the parameter plane we consider the values at which the limit cycle of the slow subsystem is tangent to the Hopf bifurcation curve of the fast subsystem. The method suggested for determining the boundary is based on the dissection of the system variables into slow and fast. The strong point of the method is that it requires the integration of the slow subsystem only. An example of the application of the method for the stomatogastric neuron model [Guckenheimer J, Gueron S, Harris-Warrick RM (1993) Philos Trans R Soc Lond B 341: 345–359] is given. Received: 31 May 1999 / Accepted in revised form: 19 November 1999     

Firing patterns in a conductance-based neuron model: bifurcation, phase diagram, and chaos

Responding to various stimuli, some neurons either remain resting or can fire several distinct patterns of action potentials, such as spiking, bursting, subthreshold oscillations, and chaotic firing. In particular, Wilson’s conductance-based neocortical neuron model, derived from the Hodgkin–Huxley model, is explored to understand underlying mechanisms of the firing patterns. Phase diagrams describing boundaries between the domains of different firing patterns are obtained via extensive numerical computations. The boundaries are further studied by standard instability analyses, which demonstrates that the chaotic neural firing could develop via period-doubling and/or period- adding cascades. Sequences of the firing patterns often observed in many neural experiments are also discussed in the phase diagram framework developed. Our results lay the groundwork for wider use of the model, especially for incorporating it into neural field modeling of the brain.     

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.     

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.     

A pulse-type hardware neuron model with beating, bursting excitation and plateau potential

We proposed a pulse-type hardware neuron model. It could reproduce simple excitations, beating and bursting discharges as well as an action potential with a plateau potential observed in living membranes. The model exhibited one of these dynamics depending on parameter values of the model's circuit. They include resistance, capacitance and externally injected DC current intensity. We studied the model's dynamics based on hardware experiments and mathematical analyses. Our results showed that two inward currents introduced into the model and differences in their operating time scales determined dynamics of the model. In particular, we illustrated a mechanism of the bursting discharges generation in terms of bifurcation theory and time-dependent changes in the form of instantaneous current voltage characteristics of the model.     

Dominant ionic mechanisms explored in spiking and bursting using local low-dimensional reductions of a biophysically realistic model neuron

The large number of variables involved in many biophysical models can conceal potentially simple dynamical mechanisms governing the properties of its solutions and the transitions between them as parameters are varied. To address this issue, we extend a novel model reduction method, based on “scales of dominance,” to multi-compartment models. We use this method to systematically reduce the dimension of a two-compartment conductance-based model of a crustacean pyloric dilator (PD) neuron that exhibits distinct modes of oscillation—tonic spiking, intermediate bursting and strong bursting. We divide trajectories into intervals dominated by a smaller number of variables, resulting in a locally reduced hybrid model whose dimension varies between two and six in different temporal regimes. The reduced model exhibits the same modes of oscillation as the 16 dimensional model over a comparable parameter range, and requires fewer ad hoc simplifications than a more traditional reduction to a single, globally valid model. The hybrid model highlights low-dimensional organizing structure in the dynamics of the PD neuron, and the dependence of its oscillations on parameters such as the maximal conductances of calcium currents. Our technique could be used to build hybrid low-dimensional models from any large multi-compartment conductance-based model in order to analyze the interactions between different modes of activity.     

Modification of bursting in a Helix neuron by drugs influencing intracellular regulation of calcium level

The effect of ruthenium red, caffein and EGTA (ethyleneglycol tetraacetic acid) influencing intracellular Ca2+ level as well as that of pH-lowering was investigated on identified RPal neuron of Helix pomatia characterized by bimodal pacemaker (bursting) activity. Drugs were applied both extracellularly and intracellularly. Intracellular injection was performed from micropipettes by pressure. It was found that intracellular injection of ruthenium red, caffein, EGTA and pH-lowering caused immediate short hyperpolarization and suspension of bursting. The effect of caffein and lowering of pH was biphasic, hyperpolarization was followed by an increase of spiking. Following EGTA injection the amplitudes of interburst hyperpolarizing waves decreased, and prolongation of spikes occurred. Extracellular application of ruthenium red caused slight depolarization, while caffein produced mainly effects that were similar to those of the intracellular injection. Adding EGTA into the bath resulted in cessation of bursting, and later on also spike generation was blocked. All these effects could be eliminated by washing. It is concluded that Ca-influx during spiking cannot be considered as a single factor in maintaining bursting activity, nevertheless, intracellular binding and liberation of Ca depending on the cell metabolism should also be taken into consideration as a possible mechanism of burst regulation.     

Dissection of a model for neuronal parabolic bursting

We have obtained new insight into the mechanisms for bursting in a class of theoretical models. We study Plant's model [24] for Aplysia R-15 to illustrate our view of these so-called parabolic bursters, which are characterized by low spike frequency at the beginning and end of a burst. By identifying and analyzing the fast and slow processes we show how they interact mutually to generate spike activity and the slow wave which underlies the burst pattern. Our treatment is essentially the first step of a singular perturbation approach presented from a geometrical viewpoint and carried out numerically with AUTO [12]. We determine the solution sets (steady state and oscillatory) of the fast subsystem with the slow variables treated as parameters. These solutions form the slow manifold over which the slow dynamics then define a burst trajectory. During the silent phase of a burst, the solution trajectory lies approximately on the steady state branch of the slow manifold and during the active phase of spiking, the trajectory sweeps through the oscillation branch. The parabolic nature of bursting arises from the (degenerate) homoclinic transition between the oscillatory branch and the steady state branch. We show that, for some parameter values, the trajectory remains strictly on the steady state branch (to produce a resting steady state or a pure slow wave without spike activity) or strictly in the oscillatory branch (continuous spike activity without silent phases). Plant's model has two slow variables: a calcium conductance and the intracellular free calcium concentration, which activates a potassium conductance. We also show how bursting arises from an alternative mechanism in which calcium inactivates the calcium current and the potassium conductance is insensitive to calcium. These and other biophysical interpretations are discussed.     

Complete chaos in a simple epidemiological model

A first-order difference equation which constitutes a simple model for a lethal parasite-host interaction is studied. Completing a study initiated by May and Anderson, the dynamics are shown to be completely chaotic.     

Transient responses of a modeled bursting neuron: analysis with equilibrium and averaged nullclines

We utilized a state-space approach to study the dynamics of a modeled bursting neuron consisting of 11 state variables. Such an approach may be used on a high-order system when a small number of variables are rate-limiting and dominate the dynamics of the model. Calculation of equilibrium and averaged nullclines and saddle-node bifurcations of the full and reduced models provided measures that indicated the transition between silence and spiking and the dynamics of the system during both the silent and spiking phases of the burst cycle. The relative stability of tonic beating solutions in the presence and absence of 5-HT was calculated in the state-space of the slow variables and related to specific biophysical mechanisms. The results were compared with similar simulations performed in Butera et al. (1995) which utilized a current-voltage (I-V)-based method for analysis. While the state-space method is sometimes more difficult to link to specific biophysical mechanisms, it offers a wider portrait of the dynamics of the system. In contrast, the use of I-V plots offers a direct relationship to biophysical processes, but provides no information on the dynamics of non-voltage-dependent processes such as Ca. Received: 6 December 1996 / Accepted in revised form: 1 July 1997     

Homoclinic chaos in a model of natural selection

 We modify a simple mathematical model for natural selection originally formulated by Robert M. May in 1983 by permitting one homozygote to have a larger selective advantage when rare than the other, and show that the new model exhibits dynamical chaos. We determine an open region of parameter space associated with homoclinic points, and prove that there are infinite sequences of period-doubling bifurcations along selected paths through parameter space. We also discuss the possibility of chaos arising from imbalance in the homozygote fitnesses in more realistic biological situations, beyond the constraints of the model. Received 3 February 1995; received in revised form 1 November 1995