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

Sparse Firing in a Hybrid Central Pattern Generator for Spinal Motor Circuits

Central pattern generators are circuits generating rhythmic movements, such as walking. The majority of existing computational models of these circuits produce antagonistic output where all neurons within a population spike with a broad burst at about the same neuronal phase with respect to network output. However, experimental recordings reveal that many neurons within these circuits fire sparsely, sometimes as rarely as once within a cycle. Here we address the sparse neuronal firing and develop a model to replicate the behavior of individual neurons within rhythm-generating populations to increase biological plausibility and facilitate new insights into the underlying mechanisms of rhythm generation. The developed network architecture is able to produce sparse firing of individual neurons, creating a novel implementation for exploring the contribution of network architecture on rhythmic output. Furthermore, the introduction of sparse firing of individual neurons within the rhythm-generating circuits is one of the factors that allows for a broad neuronal phase representation of firing at the population level. This moves the model toward recent experimental findings of evenly distributed neuronal firing across phases among individual spinal neurons. The network is tested by methodically iterating select parameters to gain an understanding of how connectivity and the interplay of excitation and inhibition influence the output. This knowledge can be applied in future studies to implement a biologically plausible rhythm-generating circuit for testing biological hypotheses.

Bistability and Chaos Emergence in Spontaneous Dynamics of Astrocytic Calcium Concentration

In this work, we consider a mathematical model describing spontaneous calcium signaling in astrocytes. Based on biologically relevant principles, this model simulates experimentally observed calcium oscillations and can predict the emergence of complicated dynamics. Using analytical and numerical analysis, various attracting sets were found and investigated. Employing bifurcation theory analysis, we examined steady state solutions, bistability, simple and complicated periodic limit cycles and also chaotic attractors. We found that astrocytes possess a variety of complex dynamical modes, including chaos and multistability, that can further provide different modulations of neuronal circuits, enhancing their plasticity and flexibility.

Temperature compensation in a small rhythmic circuit

Temperature affects the conductances and kinetics of the ionic channels that underlie neuronal activity. Each membrane conductance has a different characteristic temperature sensitivity, which raises the question of how neurons and neuronal circuits can operate robustly over wide temperature ranges. To address this, we employed computational models of the pyloric network of crabs and lobsters. We produced multiple different models that exhibit a triphasic pyloric rhythm over a range of temperatures and explored the dynamics of their currents and how they change with temperature. Temperature can produce smooth changes in the relative contributions of the currents to neural activity so that neurons and networks undergo graceful transitions in the mechanisms that give rise to their activity patterns. Moreover, responses of the models to deletions of a current can be different at high and low temperatures, indicating that even a well-defined genetic or pharmacological manipulation may produce qualitatively distinct effects depending on the temperature.

Propensity for Bistability of Bursting and Silence in the Leech Heart Interneuron

The coexistence of neuronal activity regimes has been reported under normal and pathological conditions. Such multistability could enhance the flexibility of the nervous system and has many implications for motor control, memory, and decision making. Multistability is commonly promoted by neuromodulation targeting specific membrane ionic currents. Here, we investigated how modulation of different ionic currents could affect the neuronal propensity for bistability. We considered a leech heart interneuron model. It exhibits bistability of bursting and silence in a narrow range of the leak current parameters, conductance (g_leak) and reversal potential (E_leak). We assessed the propensity for bistability of the model by using bifurcation diagrams. On the diagram (g_leak, E_leak), we mapped bursting and silent regimes. For the canonical value of E_leak we determined the range of g_leak which supported the bistability. We use this range as an index of propensity for bistability. We investigated how this index was affected by alterations of ionic currents. We systematically changed their conductances, one at a time, and built corresponding bifurcation diagrams in parameter planes of the maximal conductance of a given current and the leak conductance. We found that conductance of only one current substantially affected the index of propensity; the increase of the maximal conductance of the hyperpolarization-activated cationic current increased the propensity index. The second conductance with the strongest effect was the conductance of the low-threshold fast Ca²⁺ current; its reduction increased the propensity index although the effect was about two times smaller in magnitude. Analyzing the model with both changes applied simultaneously, we found that the diagram (g_leak, E_leak) showed a progressively expanded area of bistability of bursting and silence.

Coexistence between silent and bursting states in a biophysical Hodgkin-Huxley-type of model

Classification of the dynamical mechanisms that support bistability between bursting oscillations and silence has not yet been clarified in detail. The purpose of this paper is to demonstrate that the coexistence of a stable equilibrium point with a state of continuous bursting can occur in a slightly modified, biophysical model that describe the dynamics of pancreatic beta-cells. To realize this form of coexistence, we have introduced an additional voltage-dependent potassium current that is activated in the region around the original, unstable equilibrium point. It is interesting to note that this modification also leads the model to display a blue-sky catastrophe in the transition region between chaotic and bursting states.

Selective monostability in multi-stable systems

We propose a robust method that allows a periodic or a chaotic multi-stable system to be transformed to a monostable system at an orbit with dominant frequency of any of the coexisting attractors. Our approach implies the selection of a particular attractor by periodic external modulation with frequency close to the dominant frequency in the power spectrum of a desired orbit and simultaneous annihilation of all other coexisting states by positive feedback, both applied to one of the system parameters. The method does not require any preliminary knowledge of the system dynamics and the phase space structure. The efficiency of the method is demonstrated in both a non-autonomous multi-stable laser with coexisting periodic orbits and an autonomous Rössler-like oscillator with coexisting chaotic attractors. The experiments with an erbium-doped fibre laser provide evidence for the robustness of the proposed method in making the system monostable at an orbit with dominant frequency of any preselected attractor.

A codimension-2 bifurcation controlling endogenous bursting activity and pulse-triggered responses of a neuron model

The dynamics of individual neurons are crucial for producing functional activity in neuronal networks. An open question is how temporal characteristics can be controlled in bursting activity and in transient neuronal responses to synaptic input. Bifurcation theory provides a framework to discover generic mechanisms addressing this question. We present a family of mechanisms organized around a global codimension-2 bifurcation. The cornerstone bifurcation is located at the intersection of the border between bursting and spiking and the border between bursting and silence. These borders correspond to the blue sky catastrophe bifurcation and the saddle-node bifurcation on an invariant circle (SNIC) curves, respectively. The cornerstone bifurcation satisfies the conditions for both the blue sky catastrophe and SNIC. The burst duration and interburst interval increase as the inverse of the square root of the difference between the corresponding bifurcation parameter and its bifurcation value. For a given set of burst duration and interburst interval, one can find the parameter values supporting these temporal characteristics. The cornerstone bifurcation also determines the responses of silent and spiking neurons. In a silent neuron with parameters close to the SNIC, a pulse of current triggers a single burst. In a spiking neuron with parameters close to the blue sky catastrophe, a pulse of current temporarily silences the neuron. These responses are stereotypical: the durations of the transient intervals-the duration of the burst and the duration of latency to spiking-are governed by the inverse-square-root laws. The mechanisms described here could be used to coordinate neuromuscular control in central pattern generators. As proof of principle, we construct small networks that control metachronal-wave motor pattern exhibited in locomotion. This pattern is determined by the phase relations of bursting neurons in a simple central pattern generator modeled by a chain of oscillators.

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.

It is often not appreciated that different activity regimes can coexist with each other in a given neuron so that a transient stimulus can cause a persistent change of activity. Such multistability of the neuronal dynamics has in fact been shown in a variety of neurons and can play either a functional role or present a substrate for neurological diseases. We explored the propensity for multistability in a database of a leech heart interneuron model, testing each case (parameter set) in a database for multistability. We found a large proportion of multistable cases, especially the coexistence of silent and bursting regimes. This was a surprising result, since these cells pace the heartbeat of the leech, and the coexistence of silence and bursting could disrupt the functional pattern, threatening the viability of the leech. Analysis of networks of mutually inhibitory multistable neurons, however, showed robustness in maintaining functional activity, suggesting that the mutually inhibitory coupling can act as a protective mechanism against failures induced by multistability.