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     

Routes to chaos in a model of a bursting neuron.

Chaotic regimens have been observed experimentally in neurons as well as in deterministic neuronal models. The R15 bursting cell in the abdominal ganglion of Aplysia has been the subject of extensive mathematical modeling. Previously, the model of Plant and Kim has been shown to exhibit both bursting and beating modes of electrical activity. In this report, we demonstrate (a) that a chaotic regime exists between the bursting and beating modes of the model, and (b) that the model approaches chaos from both modes by a period doubling cascade. The bifurcation parameter employed is the external stimulus current. In addition to the period doubling observed in the model-generated trajectories, a period three "window" was observed, power spectra that demonstrate the approaches to chaos were generated, and the Lyaponov exponents and the fractal dimension of the chaotic attractors were calculated. Chaotic regimes have been observed in several similar models, which suggests that they are a general characteristic of cells that exhibit both bursting and beating modes.     

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.     

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.     

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     

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.     

Channel density regulation of firing patterns in a cortical neuron model

Modifying the density and distribution of ion channels in a neuron (by natural up- and downregulation or by pharmacological intervention or by spontaneous mutations) changes its activity pattern. In this investigation we analyzed how the impulse patterns are regulated by the density of voltage-gated channels in a neuron model based on voltage-clamp measurements of hippocampal interneurons. At least three distinct oscillatory patterns, associated with three distinct regions in the Na-K channel density plane, were found. A stability analysis showed that the different regions are characterized by saddle-node, double-orbit, and Hopf-bifurcation threshold dynamics, respectively. Single, strongly graded action potentials occur in an area outside the oscillatory regions, but less graded action potentials occur together with repetitive firing over a considerable range of channel densities. The relationship found here between channel densities and oscillatory behavior may partly explain the difference between the principal spiking patterns previously described for crab axons (class 1 and 2) and cortical neurons (regular firing and fast spiking).     

Modulation of endogenous bursting pacemaker activity in a Helix pomatia neuron

Considerable membrane depolarization was shown to arise periodically, at intervals of up to a few minutes, in the PPa1 bursting neuron ofHelix pomatia. Pulses of slow depolarizing current were found by the voltage clamping method. The frequency of the pulses was independent of the holding potential. The equilibrium potential for the slow depolarizing current was about 45 mV. During development of the depolarizing current a region of negative conductivity was observed on the steady-state voltage-current characteristic curve of the membrane. It is suggested that the pulses of slow depolarizing current are associated with the presence of secretory connections between the molluscan neurons.A. A. Bogomolets Institute of Physiology, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Neirofiziologiya, Vol. 9, No. 6, pp. 606–612, November–December, 1977.     

The interspike interval of a cable model neuron with white noise input

The firing time of a cable model neuron in response to white noise current injection is investigated with various methods. The Fourier decomposition of the depolarization leads to partial differential equations for the moments of the firing time. These are solved by perturbation and numerical methods, and the results obtained are in excellent agreement with those obtained by Monte Carlo simulation. The convergence of the random Fourier series is found to be very slow for small times so that when the firing time is small it is more efficient to simulate the solution of the stochastic cable equation directly using the two different representations of the Green's function, one which converges rapidly for small times and the other which converges rapidly for large times. The shape of the interspike interval density is found to depend strongly on input position. The various shapes obtained for different input positions resemble those for real neurons. The coefficient of variation of the interspike interval decreases monotonically as the distance between the input and trigger zone increases. A diffusion approximation for a nerve cell receiving Poisson input is considered and input/output frequency relations obtained for different input sites. The cases of multiple trigger zones and multiple input sites are briefly discussed.     

The influences of somatic and dendritic inhibition on bursting patterns in a neuronal circuit model

The balance between inhibition and excitation plays a crucial role in the generation of synchronous bursting activity in neuronal circuits. In human and animal models of epilepsy, changes in both excitatory and inhibitory synaptic inputs are known to occur. Locations and distribution of these excitatory and inhibitory synaptic inputs on pyramidal cells play a role in the integrative properties of neuronal activity, e.g., epileptiform activity. Thus the location and distribution of the inputs onto pyramidal cells are important parameters that influence neuronal activity in epilepsy. However, the location and distribution of inhibitory synapses converging onto pyramidal cells have not been fully studied. The objectives of this study are to investigate the roles of the relative location of inhibitory synapses on the dendritic tree and soma in the generation of bursting activity. We investigate influences of somatic and dendritic inhibition on bursting activity patterns in several paradigms of potential connections using a simplified multicompartmental model. We also investigate the effects of distribution of fast and slow components of GABAergic inhibition in pyramidal cells. Interspike interval (ISI) analysis is used for examination of bursting patterns. Simulations show that the inhibitory interneuron regulates neuronal bursting activity. Bursting behavior patterns depend on the synaptic weight and delay of the inhibitory connection as well as the location of the synapse. When the inhibitory interneuron synapses on the pyramidal neuron, inhibitory action is stronger if the inhibitory synapse is close to the soma. Alterations of synaptic weight of the interneuron can be compensatory for changes in the location of synaptic input. The relative changes in these parameters exert a considerable influence on whether synchronous bursting activity is facilitated or reduced. Additional simulations show that the slow GABAergic inhibitory component is more effective than the fast component in distal dendrites. Taken together, these findings illustrate the potential for GABAergic inhibition in the soma and dendritic tree to play an important modulatory role in bursting activity patterns.     

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.     

Isolation of a peptide initiating bursting activity in an identifiedHelix pomatia neuron

A peptide initiating bursting activity when applied to the soma of identified neuron RPal ofHelix pomatia (if such activity was absent) or increasing the amplitude of waves of membrane potential (if it was low), was isolated from the water-soluble fraction of brain homogenate. Application of the peptide and of the original material for its isolation (the water-soluble fraction of snail brain homogenate) evoked identical changes in the character of electrical activity of neuron RPal. It is concluded from the experimental results that the isolated protein possesses a specific action, qualitatively different from that of other known peptides, on this neuron and is the active principle of the modulating factor described previously. It is postulated that under natural conditions the modulating factor is secreted by an unidentified peptidergic interneuron in the region of axo-somatic synapses, evoking bursting activity in neuron RPal.A. A. Bogomolets Institute of Physiology, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Neirofiziologiya, Vol. 16, No. 4, pp. 488–492, July–August, 1984.     

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.     

Propagation effects of current and conductance noise in a model neuron with subthreshold oscillations

We have examined the effects of current and conductance noise in a single-neuron model which can generate a variety of physiologically important impulse patterns. Current noise enters the membrane equation directly while conductance noise is propagated through the activation variables. Additive Gaussian white noise which is implemented as conductance noise appears in the voltage equations as an additive and a multiplicative term. Moreover, the originally white noise is turned into colored noise. The noise correlation time is a function of the system's control parameters which may explain the different effects of current and conductance noise in different dynamic states. We have found the most significant, qualitative differences between different noise implementations in a pacemaker-like, tonic firing regime at the transition to chaotic burst discharges. This reflects a dynamic state of high physiological relevance.     

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.