Qualitative changes in phase-response curve and synchronization at the saddle-node-loop bifurcation

The Na+/K+-ATPase generically enables deterministic bursting in class I neurons by shearing the spike-onset bifurcation structure

Slow brain rhythms, for example during slow-wave sleep or pathological conditions like seizures and spreading depolarization, can be accompanied by oscillations in extracellular potassium concentration. Such slow brain rhythms typically have a lower frequency than tonic action-potential firing. They are assumed to arise from network-level mechanisms, involving synaptic interactions and delays, or from intrinsically bursting neurons. Neuronal burst generation is commonly attributed to ion channels with slow kinetics. Here, we explore an alternative mechanism generically available to all neurons with class I excitability. It is based on the interplay of fast-spiking voltage dynamics with a one-dimensional slow dynamics of the extracellular potassium concentration, mediated by the activity of the Na+/K+-ATPase. We use bifurcation analysis of the complete system as well as the slow-fast method to reveal that this coupling suffices to generate a hysteresis loop organized around a bistable region that emerges from a saddle-node loop bifurcation-a common feature of class I excitable neurons. Depending on the strength of the Na+/K+-ATPase, bursts are generated from pump-induced shearing the bifurcation structure, spiking is tonic, or cells are silenced via depolarization block. We suggest that transitions between these dynamics can result from disturbances in extracellular potassium regulation, such as glial malfunction or hypoxia affecting the Na+/K+-ATPase activity. The identified minimal mechanistic model outlining the sodium-potassium pump's generic contribution to burst dynamics can, therefore, contribute to a better mechanistic understanding of pathologies such as epilepsy syndromes and, potentially, inform therapeutic strategies.

How neuronal morphology impacts the synchronisation state of neuronal networks

The biophysical properties of neurons not only affect how information is processed within cells, they can also impact the dynamical states of the network. Specifically, the cellular dynamics of action-potential generation have shown relevance for setting the (de)synchronisation state of the network. The dynamics of tonically spiking neurons typically fall into one of three qualitatively distinct types that arise from distinct mathematical bifurcations of voltage dynamics at the onset of spiking. Accordingly, changes in ion channel composition or even external factors, like temperature, have been demonstrated to switch network behaviour via changes in the spike onset bifurcation and hence its associated dynamical type. A thus far less addressed modulator of neuronal dynamics is cellular morphology. Based on simplified and anatomically realistic mathematical neuron models, we show here that the extent of dendritic arborisation has an influence on the neuronal dynamical spiking type and therefore on the (de)synchronisation state of the network. Specifically, larger dendritic trees prime neuronal dynamics for in-phase-synchronised or splayed-out activity in weakly coupled networks, in contrast to cells with otherwise identical properties yet smaller dendrites. Our biophysical insights hold for generic multicompartmental classes of spiking neuron models (from ball-and-stick-type to anatomically reconstructed models) and establish a connection between neuronal morphology and the susceptibility of neural tissue to synchronisation in health and disease.

How to incorporate biological insights into network models and why it matters

Due to the staggering complexity of the brain and its neural circuitry, neuroscientists rely on the analysis of mathematical models to elucidate its function. From Hodgkin and Huxley's detailed description of the action potential in 1952 to today, new theories and increasing computational power have opened up novel avenues to study how neural circuits implement the computations that underlie behaviour. Computational neuroscientists have developed many models of neural circuits that differ in complexity, biological realism or emergent network properties. With recent advances in experimental techniques for detailed anatomical reconstructions or large-scale activity recordings, rich biological data have become more available. The challenge when building network models is to reflect experimental results, either through a high level of detail or by finding an appropriate level of abstraction. Meanwhile, machine learning has facilitated the development of artificial neural networks, which are trained to perform specific tasks. While they have proven successful at achieving task-oriented behaviour, they are often abstract constructs that differ in many features from the physiology of brain circuits. Thus, it is unclear whether the mechanisms underlying computation in biological circuits can be investigated by analysing artificial networks that accomplish the same function but differ in their mechanisms. Here, we argue that building biologically realistic network models is crucial to establishing causal relationships between neurons, synapses, circuits and behaviour. More specifically, we advocate for network models that consider the connectivity structure and the recorded activity dynamics while evaluating task performance.

Gap junctions desynchronize a neural circuit to stabilize insect flight

Insect asynchronous flight is one of the most prevalent forms of animal locomotion used by more than 600,000 species. Despite profound insights into the motor patterns¹, biomechanics^2,3 and aerodynamics underlying asynchronous flight^4,5, the architecture and function of the central-pattern-generating (CPG) neural network remain unclear. Here, on the basis of an experiment-theory approach including electrophysiology, optophysiology, Drosophila genetics and mathematical modelling, we identify a miniaturized circuit solution with unexpected properties. The CPG network consists of motoneurons interconnected by electrical synapses that, in contrast to doctrine, produce network activity splayed out in time instead of synchronized across neurons. Experimental and mathematical evidence support a generic mechanism for network desynchronization that relies on weak electrical synapses and specific excitability dynamics of the coupled neurons. In small networks, electrical synapses can synchronize or desynchronize network activity, depending on the neuron-intrinsic dynamics and ion channel composition. In the asynchronous flight CPG, this mechanism translates unpatterned premotor input into stereotyped neuronal firing with fixed sequences of cell activation that ensure stable wingbeat power and, as we show, is conserved across multiple species. Our findings prove a wider functional versatility of electrical synapses in the dynamic control of neural circuits and highlight the relevance of detecting electrical synapses in connectomics.

Temperature elevations can induce switches to homoclinic action potentials that alter neural encoding and synchronization

Almost seventy years after the discovery of the mechanisms of action potential generation, some aspects of their computational consequences are still not fully understood. Based on mathematical modeling, we here explore a type of action potential dynamics - arising from a saddle-node homoclinic orbit bifurcation - that so far has received little attention. We show that this type of dynamics is to be expected by specific changes in common physiological parameters, like an elevation of temperature. Moreover, we demonstrate that it favours synchronization patterns in networks - a feature that becomes particularly prominent when system parameters change such that homoclinic spiking is induced. Supported by in-vitro hallmarks for homoclinic spikes in the rodent brain, we hypothesize that the prevalence of homoclinic spikes in the brain may be underestimated and provide a missing link between the impact of biophysical parameters on abrupt transitions between asynchronous and synchronous states of electrical activity in the brain.

A dynamic clamp protocol to artificially modify cell capacitance

Dynamics of excitable cells and networks depend on the membrane time constant, set by membrane resistance and capacitance. Whereas pharmacological and genetic manipulations of ionic conductances of excitable membranes are routine in electrophysiology, experimental control over capacitance remains a challenge. Here, we present capacitance clamp, an approach that allows electrophysiologists to mimic a modified capacitance in biological neurons via an unconventional application of the dynamic clamp technique. We first demonstrate the feasibility to quantitatively modulate capacitance in a mathematical neuron model and then confirm the functionality of capacitance clamp in in vitro experiments in granule cells of rodent dentate gyrus with up to threefold virtual capacitance changes. Clamping of capacitance thus constitutes a novel technique to probe and decipher mechanisms of neuronal signaling in ways that were so far inaccessible to experimental electrophysiology.

Noise-induced quasiperiod and period switching

We employ a typical genetic circuit model to explore how noise can influence dynamic structure. With the increase of a key interactive parameter, the model will deterministically go through two bifurcations and three dynamic structure regions. We find that a quasiperiodic component, which is not allowed by deterministic dynamics, will be generated by noise inducing in the first two regions, and this quasiperiod will be more and more stable along with the increase in noise. In particular, in the second region the quasiperiod will compete with a stable limit cycle and perform a new transient rhythm. Furthermore, we ascertain the entropy production rate and the heat dissipation rate, and discover a minimal value with theoretical elucidation. In the end, we unveil the mechanism of the formation of quasiperiods, and show a practical biological example. We expect this work to be helpful in solving some biological or ecological problems, such as the genetic origin of periodical cicadas and population dynamics with fluctuation.

Biophysical models of intrinsic homeostasis: Firing rates and beyond

In view of ever-changing conditions both in the external world and in intrinsic brain states, maintaining the robustness of computations poses a challenge, adequate solutions to which we are only beginning to understand. At the level of cell-intrinsic properties, biophysical models of neurons permit one to identify relevant physiological substrates that can serve as regulators of neuronal excitability and to test how feedback loops can stabilize crucial variables such as long-term calcium levels and firing rates. Mathematical theory has also revealed a rich set of complementary computational properties arising from distinct cellular dynamics and even shaping processing at the network level. Here, we provide an overview over recently explored homeostatic mechanisms derived from biophysical models and hypothesize how multiple dynamical characteristics of cells, including their intrinsic neuronal excitability classes, can be stably controlled.

Activity-mediated accumulation of potassium induces a switch in firing pattern and neuronal excitability type

During normal neuronal activity, ionic concentration gradients across a neuron’s membrane are often assumed to be stable. Prolonged spiking activity, however, can reduce transmembrane gradients and affect voltage dynamics. Based on mathematical modeling, we investigated the impact of neuronal activity on ionic concentrations and, consequently, the dynamics of action potential generation. We find that intense spiking activity on the order of a second suffices to induce changes in ionic reversal potentials and to consistently induce a switch from a regular to an intermittent firing mode. This transition is caused by a qualitative alteration in the system’s voltage dynamics, mathematically corresponding to a co-dimension-two bifurcation from a saddle-node on invariant cycle (SNIC) to a homoclinic orbit bifurcation (HOM). Our electrophysiological recordings in mouse cortical pyramidal neurons confirm the changes in action potential dynamics predicted by the models: (i) activity-dependent increases in intracellular sodium concentration directly reduce action potential amplitudes, an effect typically attributed solely to sodium channel inactivation; (ii) extracellular potassium accumulation switches action potential generation from tonic firing to intermittently interrupted output. Thus, individual neurons may respond very differently to the same input stimuli, depending on their recent patterns of activity and/or the current brain-state.

Ionic concentrations in the brain are not constant. We show that during intense neuronal activity, they can change on the order of seconds and even switch neuronal spiking patterns under identical stimulation from a regular firing mode to an intermittently interrupted one. Triggered by an accumulation of extracellular potassium, such a transition is caused by a specific, qualitative change in of the neuronal voltage dynamics—a so-called bifurcation—which affects crucial features of action-potential generation and bears consequences for how information is encoded and how neurons behave together in the network. Also, changes in intracellular sodium can induce measurable effects, like a reduction of spike amplitude that occurs independently of the fast amplitude effects attributed to sodium channel inactivation. Taken together, our results demonstrate that a neuron can respond very differently to the same stimulus, depending on its previous activity or the current brain state. This finding may be particularly relevant when other regulatory mechanisms of ionic homeostasis are challenged, for example, during pathological states of glial impairment or oxygen deprivation. Finally, categorization of cortical neurons as intrinsically bursting or regular spiking may be biased by the ionic concentrations at the time of the observation, highlighting the non-static nature of neuronal dynamics.

Firing statistics in the bistable regime of neurons with homoclinic spike generation

Neuronal voltage dynamics of regularly firing neurons typically has one stable attractor: either a fixed point (like in the subthreshold regime) or a limit cycle that defines the tonic firing of action potentials (in the suprathreshold regime). In two of the three spike onset bifurcation sequences that are known to give rise to all-or-none type action potentials, however, the resting-state fixed point and limit cycle spiking can coexist in an intermediate regime, resulting in bistable dynamics. Here, noise can induce switches between the attractors, i.e., between rest and spiking, and thus increase the variability of the spike train compared to neurons with only one stable attractor. Qualitative features of the resulting spike statistics depend on the spike onset bifurcations. This paper focuses on the creation of the spiking limit cycle via the saddle-homoclinic orbit (HOM) bifurcation and derives interspike interval (ISI) densities for a conductance-based neuron model in the bistable regime. The ISI densities of bistable homoclinic neurons are found to be unimodal yet distinct from the inverse Gaussian distribution associated with the saddle-node-on-invariant-cycle bifurcation. It is demonstrated that for the HOM bifurcation the transition between rest and spiking is mainly determined along the downstroke of the action potential-a dynamical feature that is not captured by the commonly used reset neuron models. The deduced spike statistics can help to identify HOM dynamics in experimental data.

Predicting critical ignition in slow-fast excitable models

Linearization around unstable traveling waves in excitable systems can be used to approximate strength-extent curves in the problem of initiation of excitation waves for a family of spatially confined perturbations to the rest state. This theory relies on the knowledge of the unstable traveling wave solution as well as the leading left and right eigenfunctions of its linearization. We investigate the asymptotics of these ingredients, and utility of the resulting approximations of the strength-extent curves, in the slow-fast limit in two-component excitable systems of FitzHugh-Nagumo type and test those on four illustrative models. Of these, two are with degenerate dependence of the fast kinetic on the slow variable, a feature which is motivated by a particular model found in the literature. In both cases, the unstable traveling wave solution converges to a stationary "critical nucleus" of the corresponding one-component fast subsystem. We observe that in the full system, the asymptotics of the left and right eigenspaces are distinct. In particular, the slow component of the left eigenfunction corresponding to the translational symmetry does not become negligible in the asymptotic limit. This has a significant detrimental effect on the critical curve predictions. The theory as formulated previously uses an heuristic to address a difficulty related to the translational invariance. We describe two alternatives to that heuristic, which do not use the misbehaving eigenfunction component. These new heuristics show much better predictive properties, including in the asymptotic limit, in all four examples.

How to correctly quantify neuronal phase-response curves from noisy recordings

At the level of individual neurons, various coding properties can be inferred from the input-output relationship of a cell. For small inputs, this relation is captured by the phase-response curve (PRC), which measures the effect of a small perturbation on the timing of the subsequent spike. Experimentally, however, an accurate experimental estimation of PRCs is challenging. Despite elaborate measurement efforts, experimental PRC estimates often cannot be related to those from modeling studies. In particular, experimental PRCs rarely resemble the characteristic theoretical PRC expected close to spike initiation, which is indicative of the underlying spike-onset bifurcation. Here, we show for conductance-based model neurons that the correspondence between theoretical and measured phase-response curve is lost when the stimuli used for the estimation are too large. In this case, the derived phase-response curve is distorted beyond recognition and takes on a generic shape that reflects the measurement protocol and masks the spike-onset bifurcation. We discuss how to identify appropriate stimulus strengths for perturbation and noise-stimulation methods, which permit to estimate PRCs that reliably reflect the spike-onset bifurcation - a task that is particularly difficult if a lower bound for the stimulus amplitude is dictated by prominent intrinsic neuronal noise.