Geometric analysis of transient bursts

Bursting in cerebellar stellate cells induced by pharmacological agents: Non-sequential spike adding

Cerebellar stellate cells (CSCs) are spontaneously active, tonically firing (5-30 Hz), inhibitory interneurons that synapse onto Purkinje cells. We previously analyzed the excitability properties of CSCs, focusing on four key features: type I excitability, non-monotonic first-spike latency, switching in responsiveness and runup (i.e., temporal increase in excitability during whole-cell configuration). In this study, we extend this analysis by using whole-cell configuration to show that these neurons can also burst when treated with certain pharmacological agents separately or jointly. Indeed, treatment with 4-Aminopyridine (4-AP), a partial blocker of delayed rectifier and A-type K⁺ channels, at low doses induces a bursting profile in CSCs significantly different than that produced at high doses or when it is applied at low doses but with cadmium (Cd²⁺), a blocker of high voltage-activated (HVA) Ca²⁺ channels. By expanding a previously revised Hodgkin–Huxley type model, through the inclusion of Ca²⁺-activated K⁺ (K(Ca)) and HVA currents, we explain how these bursts are generated and what their underlying dynamics are. Specifically, we demonstrate that the expanded model preserves the four excitability features of CSCs, as well as captures their bursting patterns induced by 4-AP and Cd²⁺. Model investigation reveals that 4-AP is potentiating HVA, inducing square-wave bursting at low doses and pseudo-plateau bursting at high doses, whereas Cd²⁺ is potentiating K(Ca), inducing pseudo-plateau bursting when applied in combination with low doses of 4-AP. Using bifurcation analysis, we show that spike adding in square-wave bursts is non-sequential when gradually changing HVA and K(Ca) maximum conductances, delayed Hopf is responsible for generating the plateau segment within the active phase of pseudo-plateau bursts, and bursting can become “chaotic” when HVA and K(Ca) maximum conductances are made low and high, respectively. These results highlight the secondary effects of the drugs applied and suggest that CSCs have all the ingredients needed for bursting.

Excitable cells, including neurons, fire action potentials (APs) in their membrane voltage that allow them to communicate with each other and to serve certain physiological purposes. They do so either tonically by firing APs periodically, or episodically by repeatedly firing clusters of APs (called bursts) separated by quiescent periods. Each one of those firing patterns can be neuron-specific and dependent on synaptic inputs and/or their physiological environment. Cerebellar stellate cells (CSCs) that synapse onto Purkinje cells, the sole output of the cerebellum responsible for motor control, are spontaneously active inhibitory interneurons that fire APs tonically. We previously studied the excitability properties of these neurons and showed that they possess several important key features, including type I excitability, runup, non-monotonic first spike latency and switching in responsiveness. In this study, we show that CSCs can also exhibit two modes of burst firing, called square-wave and pseudo-plateau, when treated with certain pharmacological agents. Using bifurcation theory, we demonstrate that spike adding in the square-wave burst is non-sequential, changing by several spikes when certain conductances are altered gradually. This study thus sheds lights onto the overall effects of the pharmacological agents and highlights the ability of CSCs to burst in certain biological conditions.

Homoclinic organization in the Hindmarsh-Rose model: A three parameter study

Bursting phenomena are found in a wide variety of fast-slow systems. In this article, we consider the Hindmarsh-Rose neuron model, where, as it is known in the literature, there are homoclinic bifurcations involved in the bursting dynamics. However, the global homoclinic structure is far from being fully understood. Working in a three-parameter space, the results of our numerical analysis show a complex atlas of bifurcations, which extends from the singular limit to regions where a fast-slow perspective no longer applies. Based on this information, we propose a global theoretical description. Surfaces of codimension-one homoclinic bifurcations are exponentially close to each other in the fast-slow regime. Remarkably, explained by the specific properties of these surfaces, we show how the Hindmarsh-Rose model exhibits isolas of homoclinic bifurcations when appropriate two-dimensional slices are considered in the three-parameter space. On the other hand, these homoclinic bifurcation surfaces contain curves corresponding to parameter values where additional degeneracies are exhibited. These codimension-two bifurcation curves organize the bifurcations associated with the spike-adding process and they behave like the "spines-of-a-book," gathering "pages" of bifurcations of periodic orbits. Depending on how the parameter space is explored, homoclinic phenomena may be absent or far away, but their organizing role in the bursting dynamics is beyond doubt, since the involved bifurcations are generated in them. This is shown in the global analysis and in the proposed theoretical scheme.

Understanding Epileptiform After-Discharges as Rhythmic Oscillatory Transients

Electro-cortical activity in patients with epilepsy may show abnormal rhythmic transients in response to stimulation. Even when using the same stimulation parameters in the same patient, wide variability in the duration of transient response has been reported. These transients have long been considered important for the mapping of the excitability levels in the epileptic brain but their dynamic mechanism is still not well understood. To investigate the occurrence of abnormal transients dynamically, we use a thalamo-cortical neural population model of epileptic spike-wave activity and study the interaction between slow and fast subsystems. In a reduced version of the thalamo-cortical model, slow wave oscillations arise from a fold of cycles (FoC) bifurcation. This marks the onset of a region of bistability between a high amplitude oscillatory rhythm and the background state. In vicinity of the bistability in parameter space, the model has excitable dynamics, showing prolonged rhythmic transients in response to suprathreshold pulse stimulation. We analyse the state space geometry of the bistable and excitable states, and find that the rhythmic transient arises when the impending FoC bifurcation deforms the state space and creates an area of locally reduced attraction to the fixed point. This area essentially allows trajectories to dwell there before escaping to the stable steady state, thus creating rhythmic transients. In the full thalamo-cortical model, we find a similar FoC bifurcation structure. Based on the analysis, we propose an explanation of why stimulation induced epileptiform activity may vary between trials, and predict how the variability could be related to ongoing oscillatory background activity. We compare our dynamic mechanism with other mechanisms (such as a slow parameter change) to generate excitable transients, and we discuss the proposed excitability mechanism in the context of stimulation responses in the epileptic cortex.

Introduction to focus issue: rhythms and dynamic transitions in neurological disease: modeling, computation, and experiment

Rhythmic neuronal oscillations across a broad range of frequencies, as well as spatiotemporal phenomena, such as waves and bumps, have been observed in various areas of the brain and proposed as critical to brain function. While there is a long and distinguished history of studying rhythms in nerve cells and neuronal networks in healthy organisms, the association and analysis of rhythms to diseases are more recent developments. Indeed, it is now thought that certain aspects of diseases of the nervous system, such as epilepsy, schizophrenia, Parkinson's, and sleep disorders, are associated with transitions or disruptions of neurological rhythms. This focus issue brings together articles presenting modeling, computational, analytical, and experimental perspectives about rhythms and dynamic transitions between them that are associated to various diseases.