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.     

Bifurcation analysis of a two-compartment hippocampal pyramidal cell model

The Pinsky-Rinzel model is a non-smooth 2-compartmental CA3 pyramidal cell model that has been used widely within the field of neuroscience. Here we propose a modified (smooth) system that captures the qualitative behaviour of the original model, while allowing the use of available, numerical continuation methods to perform full-system bifurcation and fast-slow analysis. We study the bifurcation structure of the full system as a function of the applied current and the maximal calcium conductance. We identify the bifurcations that shape the transitions between resting, bursting and spiking behaviours, and which lead to the disappearance of bursting when the calcium conductance is reduced. Insights gained from this analysis, are then used to firstly illustrate how the irregular spiking activity found between bursting and stable spiking states, can be influenced by phase differences in the calcium and dendritic voltage, which lead to corresponding changes in the calcium-sensitive potassium current. Furthermore, we use fast-slow analysis to investigate the mechanisms of bursting and show that bursting in the model is dependent on the intermediately slow variable, calcium, while the other slow variable, the activation gate of the afterhyperpolarisation current, does not contribute to setting the intraburst dynamics but participates in setting the interburst interval. Finally, we discuss how some of the described bifurcations affect spiking behaviour, during sharp-wave ripples, in a larger network of Pinsky-Rinzel cells.     

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.     

Mixed mode oscillations as a mechanism for pseudo-plateau bursting

We combine bifurcation analysis with the theory of canard-induced mixed mode oscillations to investigate the dynamics of a novel form of bursting. This bursting oscillation, which arises from a model of the electrical activity of a pituitary cell, is characterized by small impulses or spikes riding on top of an elevated voltage plateau. Oscillations with these characteristics have been called “pseudo-plateau bursting”. Unlike standard bursting, the subsystem of fast variables does not possess a stable branch of periodic spiking solutions, and in the case studied here the standard fast/slow analysis provides little information about the underlying dynamics. We demonstrate that the bursting is actually a canard-induced mixed mode oscillation, and use canard theory to characterize the dynamics of the oscillation. We also use bifurcation analysis of the full system of equations to extend the results of the singular analysis to the physiological regime. This demonstrates that the combination of these two analysis techniques can be a powerful tool for understanding the pseudo-plateau bursting oscillations that arise in electrically excitable pituitary cells and isolated pancreatic β-cells.     

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.     

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.     

Complex dynamics in a three-level trophic system with intraspecies interaction

In this paper, we present a three-level trophic food chain, including intraspecies interaction. In contrast with other analyses, we consider the effect on the third trophic level by the first-level parameters. The model shows complex, as well as, chaotic oscillations. Bifurcation diagrams show period doubling route to chaos and crises. Also from the forward and backwards sections of the bifurcation diagrams, we find hysteresis. This result implies the coexistence of attractors for the same parameter values. In particular, we consider the coexistence of a chaotic and a P1 attractors. Our results show that the regulation in the food chain is not exclusive to either a food-prey or prey-predator interaction, but to a more subtle food-prey-predator interaction, where, for some parameter values, a food-prey or a prey-predator regulation may dominate the system's dynamics. Finally, we consider the impact of the intraspecies interaction in the overall dynamics of the food chain.     

Ion concentration dynamics as a mechanism for neuronal bursting

We describe a simple conductance-based model neuron that includes intra- and extracellular ion concentration dynamics and show that this model exhibits periodic bursting. The bursting arises as the fast-spiking behavior of the neuron is modulated by the slow oscillatory behavior in the ion concentration variables and vice versa. By separating these time scales and studying the bifurcation structure of the neuron, we catalog several qualitatively different bursting profiles that are strikingly similar to those seen in experimental preparations. Our work suggests that ion concentration dynamics may play an important role in modulating neuronal excitability in real biological systems.     

Bursting, beating, and chaos in an excitable membrane model.

We have studied periodic as well as aperiodic behavior in the self-sustained oscillations exhibited by the Hodgkin-Huxley type model of Chay, T. R., and J. Keizer (Biophys. J., 1983, 42:181-190) for the pancreatic beta-cell. Numerical solutions reveal a variety of patterns as the glucose-dependent parameter kCa is varied. These include regimes of periodic beating (continuous spiking) and bursting modes and, in the transition between these modes, aperiodic responses. Such aperiodic behavior for a nonrandom system has been called deterministic chaos and is characterized by distinguishing features found in previous studies of chaos in nonbiophysical systems and here identified for an (endogenously active) excitable membrane model. To parallel the successful analysis of chaos in other physical/chemical contexts we introduce a simplified, but quantitative, one-variable, discrete-time representation of the dynamics. It describes the evolution of intracellular calcium (which activates a potassium conductance) from one spike upstroke to the next and exhibits the various modes of behavior.     

Finding complex oscillatory phenomena in biochemical systems. An empirical approach

Starting with a model for a product-activated enzymatic reaction proposed for glycolytic oscillations, we show how more complex oscillatory phenomena may develop when the basic model is modified by addition of product recycling into substrate or by coupling in parallel or in series two autocatalytic enzyme reactions. Among the new modes of behavior are the coexistence between two stable types of oscillations (birhythmicity), bursting, and aperiodic oscillations (chaos). On the basis of these results, we outline an empirical method for finding complex oscillatory phenomena in autonomous biochemical systems, not subjected to forcing by a periodic input. This procedure relies on finding in parameter space two domains of instability of the steady state and bringing them close to each other until they merge. Complex phenomena occur in or near the region where the two domains overlap. The method applies to the search for birhythmicity, bursting and chaos in a model for the cAMP signalling system of Dictyostelium discoideum amoebae.     

Six types of multistability in a neuronal model based on slow calcium current

Background

Multistability of oscillatory and silent regimes is a ubiquitous phenomenon exhibited by excitable systems such as neurons and cardiac cells. Multistability can play functional roles in short-term memory and maintaining posture. It seems to pose an evolutionary advantage for neurons which are part of multifunctional Central Pattern Generators to possess multistability. The mechanisms supporting multistability of bursting regimes are not well understood or classified.

Methodology/Principal Findings

Our study is focused on determining the bio-physical mechanisms underlying different types of co-existence of the oscillatory and silent regimes observed in a neuronal model. We develop a low-dimensional model typifying the dynamics of a single leech heart interneuron. We carry out a bifurcation analysis of the model and show that it possesses six different types of multistability of dynamical regimes. These types are the co-existence of 1) bursting and silence, 2) tonic spiking and silence, 3) tonic spiking and subthreshold oscillations, 4) bursting and subthreshold oscillations, 5) bursting, subthreshold oscillations and silence, and 6) bursting and tonic spiking. These first five types of multistability occur due to the presence of a separating regime that is either a saddle periodic orbit or a saddle equilibrium. We found that the parameter range wherein multistability is observed is limited by the parameter values at which the separating regimes emerge and terminate.

Conclusions

We developed a neuronal model which exhibits a rich variety of different types of multistability. We described a novel mechanism supporting the bistability of bursting and silence. This neuronal model provides a unique opportunity to study the dynamics of networks with neurons possessing different types of multistability.     

Bifurcation structure of a model of bursting pancreatic cells

One- and two-dimensional bifurcation studies of a prototypic model of bursting oscillations in pancreatic beta-cells reveal a squid-formed area of chaotic dynamics in the parameter plane, with period-doubling bifurcations on one side of the arms and saddle-node bifurcations on the other. The transition from this structure to the so-called period-adding structure is found to involve a subcritical period-doubling bifurcation and the emergence of type-III intermittency. The period-adding transition itself is not smooth but consists of a saddle-node bifurcation in which (n+1)-spike bursting behavior is born, slightly overlapping with a subcritical period-doubling bifurcation in which n-spike bursting behavior loses its stability.     

Scale invariance in biology: coincidence or footprint of a universal mechanism?

In this article, we present a self-contained review of recent work on complex biological systems which exhibit no characteristic scale. This property can manifest itself with fractals (spatial scale invariance), flicker noise or 1/f-noise where f denotes the frequency of a signal (temporal scale invariance) and power laws (scale invariance in the size and duration of events in the dynamics of the system). A hypothesis recently put forward to explain these scale-free phenomomena is criticality, a notion introduced by physicists while studying phase transitions in materials, where systems spontaneously arrange themselves in an unstable manner similar, for instance, to a row of dominoes. Here, we review in a critical manner work which investigates to what extent this idea can be generalized to biology. More precisely, we start with a brief introduction to the concepts of absence of characteristic scale (power-law distributions, fractals and 1/f-noise) and of critical phenomena. We then review typical mathematical models exhibiting such properties: edge of chaos, cellular automata and self-organized critical models. These notions are then brought together to see to what extent they can account for the scale invariance observed in ecology, evolution of species, type III epidemics and some aspects of the central nervous system. This article also discusses how the notion of scale invariance can give important insights into the workings of biological systems.     

Chaos and the (un)predictability of evolution in a changing environment

Among the factors that may reduce the predictability of evolution, chaos, characterized by a strong dependence on initial conditions, has received much less attention than randomness due to genetic drift or environmental stochasticity. It was recently shown that chaos in phenotypic evolution arises commonly under frequency‐dependent selection caused by competitive interactions mediated by many traits. This result has been used to argue that chaos should often make evolutionary dynamics unpredictable. However, populations also evolve largely in response to external changing environments, and such environmental forcing is likely to influence the outcome of evolution in systems prone to chaos. We investigate how a changing environment causing oscillations of an optimal phenotype interacts with the internal dynamics of an eco‐evolutionary system that would be chaotic in a constant environment. We show that strong environmental forcing can improve the predictability of evolution by reducing the probability of chaos arising, and by dampening the magnitude of chaotic oscillations. In contrast, weak forcing can increase the probability of chaos, but it also causes evolutionary trajectories to track the environment more closely. Overall, our results indicate that, although chaos may occur in evolution, it does not necessarily undermine its predictability.     

From simple to complex oscillatory behaviour: analysis of bursting in a multiply regulated biochemical system

We analyze the transition from simple to complex oscillatory behaviour in a three-variable biochemical system that consists of the coupling in series of two autocatalytic enzyme reactions. Complex periodic behaviour occurs in the form of bursting in which clusters of spikes are separated by phases of relative quiescence. The generation of such temporal patterns is investigated by a series of complementary approaches. The dynamics of the system is first cast into two different time-scales, and one of the variables is taken as a slowly-varying parameter influencing the behaviour of the two remaining variables. This analysis shows how complex oscillations develop from simple periodic behaviour and accounts for the existence of various modes of bursting as well as for the dependence of the number of spikes per period on key parameters of the model. We further reduce the number of variables by analyzing bursting by means of one-dimensional return maps obtained from the time evolution of the three-dimensional system. The analysis of a related piecewise linear map allows for a detailed understanding of the complex sequence leading from a bursting pattern with p spikes to a pattern with p + 1 spikes per period. We show that this transition possesses properties of self-similarity associated with the occurrence of more and more complex patterns of bursting. In addition to bursting, period-doubling bifurcations leading to chaos are observed, as in the differential system, when the piecewise-linear map becomes nonlinear.     

Coherence resonance for neuronal bursting with spike undershoot

Although the bursting patterns with spike undershoot are involved with the achievement of physiological or cognitive functions of brain with synaptic noise, noise induced-coherence resonance (CR) from resting state or subthreshold oscillations instead of bursting has been widely identified to play positive roles in information process. Instead, in the present paper, CR characterized by the increase firstly and then decease of peak value of power spectrum of spike trains is evoked from a bursting pattern with spike undershoot, which means that the minimal membrane potential within burst is lower than that of the subthreshold oscillations between bursts, while CR cannot be evoked from the bursting pattern without spike undershoot. With bifurcations and fast-slow variable dissection method, the bursting patterns with and without spike undershoot are classified into “Sub-Hopf/Fold” bursting and “Fold/Homoclinic” bursting, respectively. For the bursting with spike undershoot, the trajectory of the subthreshold oscillations is very close to that of the spikes within burst. Therefore, noise can induce more spikes from the subthreshold oscillations and modulate the bursting regularity, which leads to the appearance of CR. For the bursting pattern without spike undershoot, the trajectory of the quiescent state is not close to that of the spikes within burst, and noise cannot induce spikes from the quiescent state between bursts, which is cause for non-CR. The result provides a novel case of CR phenomenon and extends the scopes of CR concept, presents that noise can enhance rather than suppress information of the bursting patterns with spike undershoot, which are helpful for understanding the dynamics and the potential physiological or cognitive functions of the nerve fiber or brain neurons with such bursting patterns.     

Classification of bursting patterns: A tale of two ducks

Bursting is one of the fundamental rhythms that excitable cells can generate either in response to incoming stimuli or intrinsically. It has been a topic of intense research in computational biology for several decades. The classification of bursting oscillations in excitable systems has been the subject of active research since the early 1980s and is still ongoing. As a by-product, it establishes analytical and numerical foundations for studying complex temporal behaviors in multiple timescale models of cellular activity. In this review, we first present the seminal works of Rinzel and Izhikevich in classifying bursting patterns of excitable systems. We recall a complementary mathematical classification approach by Bertram and colleagues, and then by Golubitsky and colleagues, which, together with the Rinzel-Izhikevich proposals, provide the state-of-the-art foundations to these classifications. Beyond classical approaches, we review a recent bursting example that falls outside the previous classification systems. Generalizing this example leads us to propose an extended classification, which requires the analysis of both fast and slow subsystems of an underlying slow-fast model and allows the dissection of a larger class of bursters. Namely, we provide a general framework for bursting systems with both subthreshold and superthreshold oscillations. A new class of bursters with at least 2 slow variables is then added, which we denote folded-node bursters, to convey the idea that the bursts are initiated or annihilated via a folded-node singularity. Key to this mechanism are so-called canard or duck orbits, organizing the underpinning excitability structure. We describe the 2 main families of folded-node bursters, depending upon the phase (active/spiking or silent/nonspiking) of the bursting cycle during which folded-node dynamics occurs. We classify both families and give examples of minimal systems displaying these novel bursting patterns. Finally, we provide a biophysical example by reinterpreting a generic conductance-based episodic burster as a folded-node burster, showing that the associated framework can explain its subthreshold oscillations over a larger parameter region than the fast subsystem approach.     

Bursting synchronization dynamics of pancreatic β-cells with electrical and chemical coupling

Based on bifurcation analysis, the synchronization behaviors of two identical pancreatic β-cells connected by electrical and chemical coupling are investigated, respectively. Various firing patterns are produced in coupled cells when a single cell exhibits tonic spiking or square-wave bursting individually, irrespectively of what the cells are connected by electrical or chemical coupling. On the one hand, cells can burst synchronously for both weak electrical and chemical coupling when an isolated cell exhibits tonic spiking itself. In particular, for electrically coupled cells, under the variation of the coupling strength there exist complex transition processes of synchronous firing patterns such as “fold/limit cycle” type of bursting, then anti-phase continuous spiking, followed by the “fold/torus” type of bursting, and finally in-phase tonic spiking. On the other hand, it is shown that when the individual cell exhibits square-wave bursting, suitable coupling strength can make the electrically coupled system generate “fold/Hopf” bursting via “fold/fold” hysteresis loop; whereas, the chemically coupled cells generate “fold/subHopf” bursting. Especially, chemically coupled bursters can exhibit inverse period-adding bursting sequence. Fast–slow dynamics analysis is applied to explore the generation mechanism of these bursting oscillations. The above analysis of bursting types and the transition may provide us with better insight into understanding the role of coupling in the dynamic behaviors of pancreatic β-cells.     

EEG动力学模型中混沌现象的研究

通过对脑电图（electroencephalogram,EEG）动力学模型模拟出的EEG信号的相图、分岔图、功率谱、关联维数和Lyapunov指数的对比研究,得出如下结论：1）该模型是按周期行为与混沌现象交替出现的间歇突发通向混沌的,且该间歇性与Hopf分岔、倍周期分岔和逆分岔有关;2）支持了EEG中存在混沌运动的观点。