Full system bifurcation analysis of endocrine bursting models

Plateau bursting is typical of many electrically excitable cells, such as endocrine cells that secrete hormones and some types of neurons that secrete neurotransmitters. Although in many of these cell types the bursting patterns are regulated by the interplay between voltage-gated calcium channels and calcium-sensitive potassium channels, they can be very different. We investigate so-called square-wave and pseudo-plateau bursting patterns found in endocrine cell models that are characterized by a super- or subcritical Hopf bifurcation in the fast subsystem, respectively. By using the polynomial model of Hindmarsh and Rose (Proceedings of the Royal Society of London B 221 (1222) 87-102), which preserves the main properties of the biophysical class of models that we consider, we perform a detailed bifurcation analysis of the full fast-slow system for both bursting patterns. We find that both cases lead to the same possibility of two routes to bursting, that is, the criticality of the Hopf bifurcation is not relevant for characterizing the route to bursting. The actual route depends on the relative location of the full-system's fixed point with respect to a homoclinic bifurcation of the fast subsystem. Our full-system bifurcation analysis reveals properties of endocrine bursting that are not captured by the standard fast-slow analysis.     

A geometric understanding of how fast activating potassium channels promote bursting in pituitary cells

The electrical activity of endocrine pituitary cells is mediated by a plethora of ionic currents and establishing the role of a single channel type is difficult. Experimental observations have shown however that fast-activating voltage- and calcium-dependent potassium (BK) current tends to promote bursting in pituitary cells. This burst promoting effect requires fast activation of the BK current, otherwise it is inhibitory to bursting. In this work, we analyze a pituitary cell model in order to answer the question of why the BK activation must be fast to promote bursting. We also examine how the interplay between the activation rate and conductance of the BK current shapes the bursting activity. We use the multiple timescale structure of the model to our advantage and employ geometric singular perturbation theory to demonstrate the origin of the bursting behaviour. In particular, we show that the bursting can arise from either canard dynamics or slow passage through a dynamic Hopf bifurcation. We then compare our theoretical predictions with experimental data using the dynamic clamp technique and find that the data is consistent with a burst mechanism due to a slow passage through a Hopf.     

State-Dependent Effects of Na Channel Noise on Neuronal Burst Generation

We explore the effects of stochastic sodium (Na) channel activation on the variability and dynamics of spiking and bursting in a model neuron. The complete model segregates Hodgin-Huxley-type currents into two compartments, and undergoes applied current-dependent bifurcations between regimes of periodic bursting, chaotic bursting, and tonic spiking. Noise is added to simulate variable, finite sizes of the population of Na channels in the fast spiking compartment.During tonic firing, Na channel noise causes variability in interspike intervals (ISIs). The variance, as well as the sensitivity to noise, depend on the model's biophysical complexity. They are smallest in an isolated spiking compartment; increase significantly upon coupling to a passive compartment; and increase again when the second compartment also includes slow-acting currents. In this full model, sufficient noise can convert tonic firing into bursting.During bursting, the actions of Na channel noise are state-dependent. The higher the noise level, the greater the jitter in spike timing within bursts. The noise makes the burst durations of periodic regimes variable, while decreasing burst length duration and variance in a chaotic regime. Na channel noise blurs the sharp transitions of spike time and burst length seen at the bifurcations of the noise-free model. Close to such a bifurcation, the burst behaviors of previously periodic and chaotic regimes become essentially indistinguishable.We discuss biophysical mechanisms, dynamical interpretations and physiological implications. We suggest that noise associated with finite populations of Na channels could evoke very different effects on the intrinsic variability of spiking and bursting discharges, depending on a biological neuron's complexity and applied current-dependent state. We find that simulated channel noise in the model neuron qualitatively replicates the observed variability in burst length and interburst interval in an isolated biological bursting neuron.     

Dynamics of in-phase and anti-phase bursting in the coupled pre-Bötzinger complex cells

Activity of neurons in the pre-Bötzinger complex within the mammalian brain stem has an important role in the generation of respiratory rhythms. Previous experimental results have shown that the dynamics of sodium and calcium within each cell may be responsible for various bursting mechanisms. In this paper, we study the bursting dynamics of the two-coupled pre-Bötzinger complex neurons. Using a combination of fast-slow decomposition and two-parameter bifurcation analysis, we explore the possible forms of dynamics that the model network can produce as well the transitions of in-phase and anti-phase bursting respectively.     

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     

神经起步点自发放电节律及节律转化的分岔规律

在神经起步点的实验中观察到了复杂多样的神经放电([Ca^2 ]o)节律模式,如周期簇放电、周期峰放电、混沌簇放电、混沌峰放电以及随机放电节律等。随着细胞外钙离子浓度的降低,神经放电节律从周期l簇放电,经过复杂的分岔过程(包括经倍周期分岔到混沌簇放电、混沌簇放电经激变到混沌峰放电、以及混沌峰放电经逆倍周期分岔到周期峰放电)转化为周期l峰放电。在神经放电理论模型——Chay模型中,调节与实验相关的参数(Ca^2 平衡电位),可以获得与实验相似的神经放电节律和节律转换规律。这表明复杂的神经放电节律之间存在着一定的分岔规律,它们是理解神经元信息编码的基础。     

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.     

Period doubling of calcium spike firing in a model of a Purkinje cell dendrite

Recordings from cerebellar Purkinje cell dendrites have revealed that in response to sustained current injection, the cell firing pattern can move from tonic firing of Ca²⁺ spikes to doublet firing and even to quadruplet firing or more complex firing. These firing patterns are not modified substantially if Na⁺ currents are blocked. We show that the experimental results can be viewed as a slow transition of the neuronal dynamics through a period-doubling bifurcation. To further support this conclusion and to understand the underlying mechanism that leads to doublet firing, we develop and study a simple, one-compartment model of Purkinje cell dendrite. The neuron can also exhibit quadruplet and chaotic firing patterns that are similar to the firing patterns that some of the Purkinje cells exhibit experimentally. The effects of parameters such as temperature, applied current, and potassium reversal potential in the model resemble their effects in experiments. The model dynamics involve three time scales. Ca²⁺- dependent K⁺ currents, with intermediate time scales, are responsible for the appearance of doublet firing, whereas a very slow hyperpolarizing current transfers the neuron from tonic to doublet firing. We use the fast-slow analysis to separate the effects of the three time scales. Fast-slow analysis of the neuronal dynamics, with the activation variable of the very slow, hyperpolarizing current considered as a parameter, reveals that the transitions occurs via a cascade of period-doubling bifurcations of the fast and intermediate subsystem as this slow variable increases. We carry out another analysis, with the Ca²⁺ concentration considered as a parameter, to investigate the conditions for the generation of doublet firing in systems with one effective variable with intermediate time scale, in which the rest state of the fast subsystem is terminated by a saddle-node bifurcation. We find that the scenario of period doubling in these systems can occur only if (1) the time scale of the intermediate variable (here, the decay rate of the calcium concentration) is slow enough in comparison with the interspike interval of the tonic firing at the transition but is not too slow and (2) there is a bistability of the fast subsystem of the spike-generating variables.     

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.     

Coexistence of Tonic Spiking Oscillations in a Leech Neuron Model

The leech neuron model studied here has a remarkable dynamical plasticity. It exhibits a wide range of activities including various types of tonic spiking and bursting. In this study we apply methods of the qualitative theory of dynamical systems and the bifurcation theory to analyze the dynamics of the leech neuron model with emphasis on tonic spiking regimes. We show that the model can demonstrate bi-stability, such that two modes of tonic spiking coexist. Under a certain parameter regime, both tonic spiking modes are represented by the periodic attractors. As a bifurcation parameter is varied, one of the attractors becomes chaotic through a cascade of period-doubling bifurcations, while the other remains periodic. Thus, the system can demonstrate co-existence of a periodic tonic spiking with either periodic or chaotic tonic spiking. Pontryagins averaging technique is used to locate the periodic orbits in the phase space.     

Nonlinear electronic circuit with neuron like bursting and spiking dynamics

It is difficult to design electronic nonlinear devices capable of reproducing complex oscillations because of the lack of general constructive rules, and because of stability problems related to the dynamical robustness of the circuits. This is particularly true for current analog electronic circuits that implement mathematical models of bursting and spiking neurons. Here we describe a novel, four-dimensional and dynamically robust nonlinear analog electronic circuit that is intrinsic excitable, and that displays frequency adaptation bursting and spiking oscillations. Despite differences from the classical Hodgkin–Huxley (HH) neuron model, its bifurcation sequences and dynamical properties are preserved, validating the circuit as a neuron model. The circuit's performance is based on a nonlinear interaction of fast–slow circuit blocks that can be clearly dissected, elucidating burst's starting, sustaining and stopping mechanisms, which may also operate in real neurons. Our analog circuit unit is easily linked and may be useful in building networks that perform in real-time.     

Symmetric bursting behaviors in the generalized FitzHugh–Nagumo model

In the current paper, we have investigated the generalized FitzHugh–Nagumo model. We have shown that symmetric bursting behaviors of different types could be observed in this model with an appropriate recovery term. A modified version of this system is used to construct bursting activities. Furthermore, we have shown some numerical examples of delayed Hopf bifurcation and canard phenomenon in the symmetric bursting of super-Hopf/homoclinic type near its super-Hopf and homoclinic bifurcations, respectively.     

Intrinsic and network rhythmogenesis in a reduced traub model for CA3 neurons

We have developed a two-compartment, eight-variable model of a CA3 pyramidal cell as a reduction of a complex 19-compartment cable model [Traub et al, 1991]. Our reduced model segregates the fast currents for sodium spiking into a proximal, soma-like, compartment and the slower calcium and calcium-mediated currents into a dendrite-like compartment. In each model periodic bursting gives way to repetitive soma spiking as somatic injected current increases. Steady dendritic stimulation can produce periodic bursting of significantly higher frequency (8–20 Hz) than can steady somatic input (<8 Hz). Bursting in our model occurs only for an intermediate range of electronic coupling conductance. It depends on the segregation of channel types and on the coupling current that flows back-and-forth between compartments. When the soma and dendrite are tightly coupled electrically, our model reduces to a single compartment and does not burst. Network simulations with our model using excitatory AMPA and NMDA synapses (without inhibition) give results similar to those obtained with the complex cable model [Traub et al, 1991; Traub et al, 1992]. Brief stimulation of a single cell in a resting network produces multiple synchronized population bursts, with fast AMPA synapses providing the dominant synchronizing mechanism. The number of bursts increases with the level of maximal NMDA conductance. For high enough maximal NMDA conductance synchronized bursting repeats indefinitely. We find that two factors can cause the cells to desynchronize when AMPA synapses are blocked: heterogeneity of properties amongst cells and intrinsically chaotic burst dynamics. But even when cells are identical, they may synchronize only approximately rather than exactly. Since our model has a limited number of parameters and variables, we have studied its cellular and network dynamics computationally with relative ease and over wide parameter ranges. Thereby, we identify some qualitative features that parallel or are distinguished from those of other neuronal systems; e.g., we discuss how bursting here differs from that in some classical models.     

Modulation of neuronal excitability by intracellular calcium buffering: from spiking to bursting

We have investigated the detailed regulation of neuronal firing pattern by the cytosolic calcium buffering capacity using a combination of mathematical modeling and patch-clamp recording in acute slice. Theoretical results show that a high calcium buffer concentration alters the characteristic regular firing of cerebellar granule cells and that a transition to various modes of oscillations occurs, including bursting. Using bifurcation analysis, we show that this transition from spiking to bursting is a consequence of the major slowdown of calcium dynamics. Patch-clamp recordings on cerebellar granule cells loaded with a high concentration of the fast calcium buffer BAPTA (15 mM) reveal dramatic alterations in their excitability as compared to cells loaded with 0.15 mM BAPTA. In high calcium buffering conditions, granule cells exhibit all bursting behaviors predicted by the model whereas bursting is never observed in low buffering. These results suggest that cytosolic calcium buffering capacity can tightly modulate neuronal firing patterns leading to generation of complex patterns and therefore that calcium-binding proteins may play a critical role in the non-synaptic plasticity and information processing in the central nervous system.     

From spikers to bursters via coupling: Help from heterogeneity

It has been shown previously that identical spiking cells, incapable of bursting by themselves, may burst under weak diffusive coupling conditions. With stronger coupling, the coupled system reverts to bursting can be recovered by introducing heterogeneity in the model parameters. For a two-cell system, we explain the phenomenon with bifurcation analysis. For larger clusters of coupled cells, we demonstrate by numerical simulation that the phenomenon carries over. In addition, we use a perturbation analysis to deduce the dependence of the heterogeneity parameter used in the bifurcation analysis on the original heterogeneity parameters and the coupling strength. Implications of the phenomenon of emergent bursting are discussed in the context of electrical activity in pancreatic beta cells.     

Spontaneous secondary spiking in excitable cells

Kepler & Marder (1993, Biol. Cybern.68, 209-214) proposed a model describing the electrical activity of a crab neuron in which a train of directly induced action potentials is sometimes followed by one or more spontaneous action potentials, referred to as spontaneous secondary spikes. We reduce their five-dimensional model to three dimensions in two different ways in order to gain insight into the mechanism underlying the spontaneous spikes. We then treat a slowly varying current as a parameter in order to give a qualitative explanation of the phenomenon using phase-plane and bifurcation analysis. We demonstrate that a three-dimensional model, consisting of a two-dimensional excitable system plus a slow inward current, is sufficient to produce the behaviour observed in the original model. The exact dynamics of the excitable system are not important, but the relative time constant and amplitude of the slow inward current are crucial. Using the numerical bifurcation analysis package AUTO (Doedel & Kernevez, 1986, AUTO: Software for Continuation and Bifurcation Problems in Ordinary Differential Equations. California Institute of Technology), we compute bifurcation diagrams using the maximum amplitude of the slow inward current as the bifurcation parameter. The full and reduced models have a stable resting potential for all values of the bifurcation parameter. At a critical value of the bifurcation parameter, a stable tonic firing mode arises via a saddle-node of periodics bifurcation. Whether or not the models can exhibit transient or continuous spontaneous spiking depends on their position in parameter space relative to this saddle-node of periodics.     

The Onset and Extinction of Neural Spiking: A Numerical Bifurcation Approach

We study the onset of neural spiking when the equilibrium rest state loses stability by the change of a critical parameter, the applied current. In the case of the well-known Morris-Lecar model, we start from a complete numerical study of the bifurcation diagram in the most relevant two-parameter range. This diagram includes all equilibrium and limit cycle bifurcations, thus correcting and completing earlier studies.We discuss and classify the behavior of the spiking orbits, when increasing or decreasing the applied current. A complete classification can be extracted from the complete bifurcation diagram. It is based on three components: bifurcation type of the equilibrium at the loss of stability, subcritical behavior in the limit of decreasing the applied current and supercritical behavior in the limit of increasing the applied current.     

神经放电节律加周期分岔序列中的簇放电到峰放电的转迁

含快慢子系统的神经元数学模型仿真预期,神经放电节律经历加周期分岔序列,可以进一步表现激变,并通过逆倍周期分岔级联进入周期1峰放电。实验调节胞外钙离子浓度,观察到从周期1簇放电开始的带有随机节律的加周期分岔到簇内有多个峰的簇放电,再经激变转迁到峰放电节律的分岔序列,提供了这种分岔序列模式实验证据。实验所见之激变表现为簇放电节律的休止期消失,放电节律变为混沌峰放电和周期峰放电。作者利用随机Chay模型更加逼真地仿真再现了实验所见的分岔序列。该实验结果验证了以前的确定性数学模型的理论预期,并利用随机理论模型仿真了其在现实神经系统的表现;揭示了一类完整的神经放电节律的转换规律。     

High Prevalence of Multistability of Rest States and Bursting in a Database of a Model Neuron

Flexibility in neuronal circuits has its roots in the dynamical richness of their neurons. Depending on their membrane properties single neurons can produce a plethora of activity regimes including silence, spiking and bursting. What is less appreciated is that these regimes can coexist with each other so that a transient stimulus can cause persistent change in the activity of a given neuron. Such multistability of the neuronal dynamics has been shown in a variety of neurons under different modulatory conditions. It can play either a functional role or present a substrate for dynamical diseases. We considered a database of an isolated leech heart interneuron model that can display silent, tonic spiking and bursting regimes. We analyzed only the cases of endogenous bursters producing functional half-center oscillators (HCOs). Using a one parameter (the leak conductance ()) bifurcation analysis, we extended the database to include silent regimes (stationary states) and systematically classified cases for the coexistence of silent and bursting regimes. We showed that different cases could exhibit two stable depolarized stationary states and two hyperpolarized stationary states in addition to various spiking and bursting regimes. We analyzed all cases of endogenous bursters and found that 18% of the cases were multistable, exhibiting coexistences of stationary states and bursting. Moreover, 91% of the cases exhibited multistability in some range of . We also explored HCOs built of multistable neuron cases with coexisting stationary states and a bursting regime. In 96% of cases analyzed, the HCOs resumed normal alternating bursting after one of the neurons was reset to a stationary state, proving themselves robust against this perturbation.