High prevalence of multistability of rest states and bursting in a database of a model neuron

Multistability of bursting rhythms in a half-center oscillator and the protective effects of synaptic inhibition

For animals to meet environmental challenges, the activity patterns of specialized oscillatory neural circuits, central pattern generators (CPGs), controlling rhythmic movements like breathing and locomotion, are adjusted by neuromodulation. As a representative example, the leech heartbeat is controlled by a CPG driven by two pairs of mutually inhibitory interneurons, heart interneuron (HN) half-center oscillators (HCO). Experiments and modeling indicate that neuromodulation of HCO navigates this CPG between dysfunctional regimes by employing a co-regulating inverted relation; reducing Na+/K+ pump current and increasing hyperpolarization-activated (h-) current. Simply reducing pump activity or increasing h-current leads to either seizure-like bursting or an asymmetric bursting dysfunctional regime, respectively. Here, we demonstrate through modeling that, alongside this coregulation path, a new bursting regime emerges. Both regimes fulfill the criteria for functional bursting activity. Although the cycle periods and burst durations of these patterns are roughly the same, the new one exhibits an intra-burst spike frequency that is twice as high as the other. This finding suggests that neuromodulation could introduce additional functional regimes with higher spike frequency, and thus more effective synaptic transmission to motor neurons. We found that this new regime co-exists with the original bursting. The HCO can be switched between them by a short pulse of excitatory or inhibitory conductance. In this domain of coexisting functional patterns, an isolated cell model exhibits only one regime, a severely dysfunctional plateau-containing, seizure-like activity. This aligns with widely reported notion that deficiency of inhibition can cause seizures and other dysfunctional neural activities. We show that along the coregulation path of neuromodulation, the high excitability of the single HNs induced by myomodulin is harnessed by mutually inhibitory synaptic interactions of the HCO into the functional bursting pattern.

Drive-specific selection in multistable mechanical networks

Systems with many stable configurations abound in nature, both in living and inanimate matter, encoding a rich variety of behaviors. In equilibrium, a multistable system is more likely to be found in configurations with lower energy, but the presence of an external drive can alter the relative stability of different configurations in unexpected ways. Living systems are examples par excellence of metastable nonequilibrium attractors whose structure and stability are highly dependent on the specific form and pattern of the energy flow sustaining them. Taking this distinctively lifelike behavior as inspiration, we sought to investigate the more general physical phenomenon of drive-specific selection in nonequilibrium dynamics. To do so, we numerically studied driven disordered mechanical networks of bistable springs possessing a vast number of stable configurations arising from the two stable rest lengths of each spring, thereby capturing the essential physical properties of a broad class of multistable systems. We found that there exists a range of forcing amplitudes for which the attractor states of driven disordered multistable mechanical networks are fine-tuned with respect to the pattern of external forcing to have low energy absorption from it. Additionally, we found that these drive-specific attractor states are further stabilized by precise matching between the multidimensional shape of their orbit and that of the potential energy well they inhabit. Lastly, we showed evidence of drive-specific selection in an experimental system and proposed a general method to estimate the range of drive amplitudes for drive-specific selection.

Initial Variability and Time-Dependent Changes of Neuronal Response Features Are Cell-Type-Specific

Different cell types are commonly defined by their distinct response features. But several studies proved substantial variability between cells of the same type, suggesting rather the appraisal of response feature distributions than a limitation to "typical" responses. Moreover, there is growing evidence that time-dependent changes of response features contribute to robust and functional network output in many neuronal systems. The individually characterized Touch (T), Pressure (P), and Retzius (Rz) cells in the medicinal leech allow for a rigid analysis of response features, elucidating differences between and variability within cell types, as well as their changes over time. The initial responses of T and P cells to somatic current injection cover a wide range of spike counts, and their first spike is generated with a high temporal precision after a short latency. In contrast, all Rz cells elicit very similar low spike counts with variable, long latencies. During prolonged electrical stimulation the resting membrane potential of all three cell types hyperpolarizes. At the same time, Rz cells reduce their spiking activity as expected for a departure from the spike threshold. In contrast, both mechanoreceptor types increase their spike counts during repeated stimulation, consistent with previous findings in T cells. A control experiment reveals that neither a massive current stimulation nor the hyperpolarization of the membrane potential is necessary for the mechanoreceptors' increase in excitability over time. These findings challenge the previously proposed involvement of slow K+-channels in the time-dependent activity changes. We also find no indication for a run-down of HCN channels over time, and a rigid statistical analysis contradicts several potential experimental confounders as the basis of the observed variability. We conclude that the time-dependent change in excitability of T and P cells could indicate a cell-type-specific shift between different spiking regimes, which also could explain the high variability in the initial responses. The underlying mechanism needs to be further investigated in more naturalistic experimental situations to disentangle the effects of varying membrane properties versus network interactions. They will show if variability in individual response features serves as flexible adaptation to behavioral contexts rather than just "randomness".

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.

Discrimination of bursts and tonic activity in multifunctional sensorimotor neural network using the extended hill-valley method

Individual neurons can exhibit a wide range of activity, including spontaneous spiking, tonic spiking, bursting, or spike-frequency adaptation, and can also transition between these activity types. Manual identification of these activity patterns can be subjective and inconsistent. The extended hill-valley (EHV) analysis discriminates tonic spiking and bursts in a spike train by detecting fluctuations in a local, history-dependent analysis signal derived from the spike train. Consequently, the EHV method is not susceptible to changes in baseline firing rate and can identify different types of activity patterns. In addition, output from the EHV method can be used to identify more complex activity patterns such as phasotonic bursting, in which a burst is immediately followed by a period of tonic spiking.NEW & NOTEWORTHY Neurons exhibit diverse spiking patterns, but automated activity classification has focused mainly on detecting bursts. The novel extended hill-valley algorithm uses a smoothed, history-dependent signal to discriminate different types of activity, such as bursts and tonic spiking.

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.

Optimizing computer models of corticospinal neurons to replicate in vitro dynamics

Corticospinal neurons (SPI), thick-tufted pyramidal neurons in motor cortex layer 5B that project caudally via the medullary pyramids, display distinct class-specific electrophysiological properties in vitro: strong sag with hyperpolarization, lack of adaptation, and a nearly linear frequency-current (F-I) relationship. We used our electrophysiological data to produce a pair of large archives of SPI neuron computer models in two model classes: 1) detailed models with full reconstruction; and 2) simplified models with six compartments. We used a PRAXIS and an evolutionary multiobjective optimization (EMO) in sequence to determine ion channel conductances. EMO selected good models from each of the two model classes to form the two model archives. Archived models showed tradeoffs across fitness functions. For example, parameters that produced excellent F-I fit produced a less-optimal fit for interspike voltage trajectory. Because of these tradeoffs, there was no single best model but rather models that would be best for particular usages for either single neuron or network explorations. Further exploration of exemplar models with strong F-I fit demonstrated that both the detailed and simple models produced excellent matches to the experimental data. Although dendritic ion identities and densities cannot yet be fully determined experimentally, we explored the consequences of a demonstrated proximal to distal density gradient of I_h, demonstrating that this would lead to a gradient of resonance properties with increased resonant frequencies more distally. We suggest that this dynamical feature could serve to make the cell particularly responsive to major frequency bands that differ by cortical layer.

NEW & NOTEWORTHY

We developed models of motor cortex corticospinal neurons that replicate in vitro dynamics, including hyperpolarization-induced sag and realistic firing patterns. Models demonstrated resonance in response to synaptic stimulation, with resonance frequency increasing in apical dendrites with increasing distance from soma, matching the increasing oscillation frequencies spanning deep to superficial cortical layers. This gradient may enable specific corticospinal neuron dendrites to entrain to relevant oscillations in different cortical layers, contributing to appropriate motor output commands.

Na(+)/K(+) pump interacts with the h-current to control bursting activity in central pattern generator neurons of leeches

The dynamics of different ionic currents shape the bursting activity of neurons and networks that control motor output. Despite being ubiquitous in all animal cells, the contribution of the Na(+)/K(+) pump current to such bursting activity has not been well studied. We used monensin, a Na(+)/H(+) antiporter, to examine the role of the pump on the bursting activity of oscillator heart interneurons in leeches. When we stimulated the pump with monensin, the period of these neurons decreased significantly, an effect that was prevented or reversed when the h-current was blocked by Cs(+). The decreased period could also occur if the pump was inhibited with strophanthidin or K(+)-free saline. Our monensin results were reproduced in model, which explains the pump's contributions to bursting activity based on Na(+) dynamics. Our results indicate that a dynamically oscillating pump current that interacts with the h-current can regulate the bursting activity of neurons and networks.

From global to local: exploring the relationship between parameters and behaviors in models of electrical excitability

Models of electrical activity in excitable cells involve nonlinear interactions between many ionic currents. Changing parameters in these models can produce a variety of activity patterns with sometimes unexpected effects. Further more, introducing new currents will have different effects depending on the initial parameter set. In this study we combined global sampling of parameter space and local analysis of representative parameter sets in a pituitary cell model to understand the effects of adding K (+) conductances, which mediate some effects of hormone action on these cells. Global sampling ensured that the effects of introducing K (+) conductances were captured across a wide variety of contexts of model parameters. For each type of K (+) conductance we determined the types of behavioral transition that it evoked. Some transitions were counterintuitive, and may have been missed without the use of global sampling. In general, the wide range of transitions that occurred when the same current was applied to the model cell at different locations in parameter space highlight the challenge of making accurate model predictions in light of cell-to-cell heterogeneity. Finally, we used bifurcation analysis and fast/slow analysis to investigate why specific transitions occur in representative individual models. This approach relies on the use of a graphics processing unit (GPU) to quickly map parameter space to model behavior and identify parameter sets for further analysis. Acceleration with modern low-cost GPUs is particularly well suited to exploring the moderate-sized (5-20) parameter spaces of excitable cell and signaling models.

Computational models in the age of large datasets

Technological advances in experimental neuroscience are generating vast quantities of data, from the dynamics of single molecules to the structure and activity patterns of large networks of neurons. How do we make sense of these voluminous, complex, disparate and often incomplete data? How do we find general principles in the morass of detail? Computational models are invaluable and necessary in this task and yield insights that cannot otherwise be obtained. However, building and interpreting good computational models is a substantial challenge, especially so in the era of large datasets. Fitting detailed models to experimental data is difficult and often requires onerous assumptions, while more loosely constrained conceptual models that explore broad hypotheses and principles can yield more useful insights.

Noise, transient dynamics, and the generation of realistic interspike interval variation in square-wave burster neurons

First return maps of interspike intervals for biological neurons that generate repetitive bursts of impulses can display stereotyped structures (neuronal signatures). Such structures have been linked to the possibility of multicoding and multifunctionality in neural networks that produce and control rhythmical motor patterns. In some cases, isolating the neurons from their synaptic network reveals irregular, complex signatures that have been regarded as evidence of intrinsic, chaotic behavior. We show that incorporation of dynamical noise into minimal neuron models of square-wave bursting (either conductance-based or abstract) produces signatures akin to those observed in biological examples, without the need for fine tuning of parameters or ad hoc constructions for inducing chaotic activity. The form of the stochastic term is not strongly constrained and can approximate several possible sources of noise, e.g., random channel gating or synaptic bombardment. The cornerstone of this signature generation mechanism is the rich, transient, but deterministic dynamics inherent in the square-wave (saddle-node and homoclinic) mode of neuronal bursting. We show that noise causes the dynamics to populate a complex transient scaffolding or skeleton in state space, even for models that (without added noise) generate only periodic activity (whether in bursting or tonic spiking mode).

Identifying crucial parameter correlations maintaining bursting activity

Recent experimental and computational studies suggest that linearly correlated sets of parameters (intrinsic and synaptic properties of neurons) allow central pattern-generating networks to produce and maintain their rhythmic activity regardless of changing internal and external conditions. To determine the role of correlated conductances in the robust maintenance of functional bursting activity, we used our existing database of half-center oscillator (HCO) model instances of the leech heartbeat CPG. From the database, we identified functional activity groups of burster (isolated neuron) and half-center oscillator model instances and realistic subgroups of each that showed burst characteristics (principally period and spike frequency) similar to the animal. To find linear correlations among the conductance parameters maintaining functional leech bursting activity, we applied Principal Component Analysis (PCA) to each of these four groups. PCA identified a set of three maximal conductances (leak current, _Leak; a persistent K current, _K2; and of a persistent Na+ current, _P) that correlate linearly for the two groups of burster instances but not for the HCO groups. Visualizations of HCO instances in a reduced space suggested that there might be non-linear relationships between these parameters for these instances. Experimental studies have shown that period is a key attribute influenced by modulatory inputs and temperature variations in heart interneurons. Thus, we explored the sensitivity of period to changes in maximal conductances of _Leak, _K2, and _P, and we found that for our realistic bursters the effect of these parameters on period could not be assessed because when varied individually bursting activity was not maintained.

Central pattern-generating networks (CPGs) must be remarkably robust, maintaining functional rhythmic activity despite fluctuations in internal and external conditions. Recent experimental evidence suggests that this robustness is achieved by the coordinated regulation of many membrane and synaptic current parameters. Experimental and computational studies showed that linearly correlated sets of such parameters allow CPG neurons to produce and maintain their rhythmic activity. However, the mechanisms that allow multiple parameters to interact, thereby producing and maintaining rhythmic single cell and network activity, are not clear. Here, we use a half-center oscillator (HCO) model that replicates the electrical activity (rhythmic alternating bursting of mutually inhibitory interneurons) of the leech heartbeat CPG to investigate potential relationships between parameters that maintain functional bursting activity in the HCOs and the isolated component neurons (bursters). We found a linearly correlated set of three maximal conductances that maintains functional bursting activity similar to the animal in burster model instances, therefore increasing robustness of bursting activity. In addition, we found that bursting activity was very sensitive to individual variation of these parameters; only correlated changes could maintain the activity type.