Ionic current correlations underlie the global tuning of large numbers of neuronal activity attributes

Alterations in network robustness upon simultaneous temperature and pH perturbations

Nervous systems have evolved to function consistently in the face of the normal environmental fluctuations experienced by animals. The stomatogastric nervous system (STNS) of the crab, Cancer borealis, produces a motor output that has been studied for its remarkable robustness in response to single global perturbations. Changes in environments, however, are often complex and multifactorial. Therefore, we studied the robustness of the pyloric network in the stomatogastric ganglion (STG) in response to simultaneous perturbations of temperature and pH. We compared the effects of elevated temperatures on the pyloric rhythm at control, acid, or base pHs. In each pH recordings were made at 11°C, and then the temperature was increased until the rhythms became disorganized ("crashed"). Pyloric burst frequencies and phase relationships showed minor differences between pH groups until reaching close to the crash temperatures. However, the temperatures at which the rhythms were disrupted were lower in the two extreme pH conditions. This indicates that one environmental stress can make an animal less resilient to a second stressor.NEW & NOTEWORTHY Resilience to environmental fluctuations is important for all animals. It is common that animals encounter multiple stressful events at the same time, the cumulative impacts of which are largely unknown. This study examines the effects of temperature and pH on the nervous system of crabs that live in the fluctuating environments of the Northern Atlantic Ocean. The ranges of tolerance to one perturbation, temperature, are reduced under the influence of a second, pH.

Heterogeneous off-target impact of ion-channel deletion on intrinsic properties of hippocampal model neurons that self-regulate calcium

How do neurons that implement cell-autonomous self-regulation of calcium react to knockout of individual ion-channel conductances? To address this question, we used a heterogeneous population of 78 conductance-based models of hippocampal pyramidal neurons that maintained cell-autonomous calcium homeostasis while receiving theta-frequency inputs. At calcium steady-state, we individually deleted each of the 11 active ion-channel conductances from each model. We measured the acute impact of deleting each conductance (one at a time) by comparing intrinsic electrophysiological properties before and immediately after channel deletion. The acute impact of deleting individual conductances on physiological properties (including calcium homeostasis) was heterogeneous, depending on the property, the specific model, and the deleted channel. The underlying many-to-many mapping between ion channels and properties pointed to ion-channel degeneracy. Next, we allowed the other conductances (barring the deleted conductance) to evolve towards achieving calcium homeostasis during theta-frequency activity. When calcium homeostasis was perturbed by ion-channel deletion, post-knockout plasticity in other conductances ensured resilience of calcium homeostasis to ion-channel deletion. These results demonstrate degeneracy in calcium homeostasis, as calcium homeostasis in knockout models was implemented in the absence of a channel that was earlier involved in the homeostatic process. Importantly, in reacquiring homeostasis, ion-channel conductances and physiological properties underwent heterogenous plasticity (dependent on the model, the property, and the deleted channel), even introducing changes in properties that were not directly connected to the deleted channel. Together, post-knockout plasticity geared towards maintaining homeostasis introduced heterogenous off-target effects on several channels and properties, suggesting that extreme caution be exercised in interpreting experimental outcomes involving channel knockouts.

Motor neurons within a network use cell-type specific feedback mechanisms to constrain relationships among ion channel mRNAs

Recently, activity has been proposed as a primary feedback mechanism used by continuously bursting neurons to coordinate ion channel mRNA relationships that underlie stable output. However, some neuron types only have intermittent periods of activity and so may require alternative mechanisms that induce and constrain the appropriate ion channel profile in different states of activity. To address this, we used the pyloric dilator (PD; constitutively active) and the lateral gastric (LG; periodically active) neurons of the stomatogastric ganglion (STG) of the crustacean Cancer borealis. We experimentally stimulated descending inputs to the STG to cause release of neuromodulators known to elicit the active state of LG neurons and quantified the mRNA abundances and pairwise relationships of 11 voltage-gated ion channels in active and silent LG neurons. The same stimulus does not significantly alter PD activity. Activation of LG upregulated ion channel mRNAs and lead to a greater number of positively correlated pairwise channel mRNA relationships. Conversely, this stimulus did not induce major changes in ion channel mRNA abundances and relationships of PD cells, suggesting their ongoing activity is sufficient to maintain channel mRNA relationships even under changing modulatory conditions. In addition, we found that ion channel mRNA correlations induced by the active state of LG are influenced by a combination of activity- and neuromodulator-dependent feedback mechanisms. Interestingly, some of these same correlations are maintained by distinct mechanisms in PD, suggesting that these motor networks use distinct feedback mechanisms to coordinate the same mRNA relationships across neuron types.NEW & NOTEWORTHY Neurons use various feedback mechanisms to adjust and maintain their output. Here, we demonstrate that different neurons within the same network can use distinct signaling mechanisms to regulate the same ion channel mRNA relationships.

Homeostatic regulation of neuronal function: importance of degeneracy and pleiotropy

Neurons maintain their average firing rate and other properties within narrow bounds despite changing conditions. This homeostatic regulation is achieved using negative feedback to adjust ion channel expression levels. To understand how homeostatic regulation of excitability normally works and how it goes awry, one must consider the various ion channels involved as well as the other regulated properties impacted by adjusting those channels when regulating excitability. This raises issues of degeneracy and pleiotropy. Degeneracy refers to disparate solutions conveying equivalent function (e.g., different channel combinations yielding equivalent excitability). This many-to-one mapping contrasts the one-to-many mapping described by pleiotropy (e.g., one channel affecting multiple properties). Degeneracy facilitates homeostatic regulation by enabling a disturbance to be offset by compensatory changes in any one of several different channels or combinations thereof. Pleiotropy complicates homeostatic regulation because compensatory changes intended to regulate one property may inadvertently disrupt other properties. Co-regulating multiple properties by adjusting pleiotropic channels requires greater degeneracy than regulating one property in isolation and, by extension, can fail for additional reasons such as solutions for each property being incompatible with one another. Problems also arise if a perturbation is too strong and/or negative feedback is too weak, or because the set point is disturbed. Delineating feedback loops and their interactions provides valuable insight into how homeostatic regulation might fail. Insofar as different failure modes require distinct interventions to restore homeostasis, deeper understanding of homeostatic regulation and its pathological disruption may reveal more effective treatments for chronic neurological disorders like neuropathic pain and epilepsy.

Pareto optimality, economy-effectiveness trade-offs and ion channel degeneracy: improving population modelling for single neurons

Neurons encounter unavoidable evolutionary trade-offs between multiple tasks. They must consume as little energy as possible while effectively fulfilling their functions. Cells displaying the best performance for such multi-task trade-offs are said to be Pareto optimal, with their ion channel configurations underpinning their functionality. Ion channel degeneracy, however, implies that multiple ion channel configurations can lead to functionally similar behaviour. Therefore, instead of a single model, neuroscientists often use populations of models with distinct combinations of ionic conductances. This approach is called population (database or ensemble) modelling. It remains unclear, which ion channel parameters in the vast population of functional models are more likely to be found in the brain. Here we argue that Pareto optimality can serve as a guiding principle for addressing this issue by helping to identify the subpopulations of conductance-based models that perform best for the trade-off between economy and functionality. In this way, the high-dimensional parameter space of neuronal models might be reduced to geometrically simple low-dimensional manifolds, potentially explaining experimentally observed ion channel correlations. Conversely, Pareto inference might also help deduce neuronal functions from high-dimensional Patch-seq data. In summary, Pareto optimality is a promising framework for improving population modelling of neurons and their circuits.

Minimal requirements for a neuron to co-regulate many properties and the implications for ion channel correlations and robustness

Neurons regulate their excitability by adjusting their ion channel levels. Degeneracy – achieving equivalent outcomes (excitability) using different solutions (channel combinations) – facilitates this regulation by enabling a disruptive change in one channel to be offset by compensatory changes in other channels. But neurons must coregulate many properties. Pleiotropy – the impact of one channel on more than one property – complicates regulation because a compensatory ion channel change that restores one property to its target value often disrupts other properties. How then does a neuron simultaneously regulate multiple properties? Here, we demonstrate that of the many channel combinations producing the target value for one property (the single-output solution set), few combinations produce the target value for other properties. Combinations producing the target value for two or more properties (the multioutput solution set) correspond to the intersection between single-output solution sets. Properties can be effectively coregulated only if the number of adjustable channels (n_in) exceeds the number of regulated properties (n_out). Ion channel correlations emerge during homeostatic regulation when the dimensionality of solution space (n_in − n_out) is low. Even if each property can be regulated to its target value when considered in isolation, regulation as a whole fails if single-output solution sets do not intersect. Our results also highlight that ion channels must be coadjusted with different ratios to regulate different properties, which suggests that each error signal drives modulatory changes independently, despite those changes ultimately affecting the same ion channels.

Reciprocally inhibitory circuits operating with distinct mechanisms are differently robust to perturbation and modulation

Reciprocal inhibition is a building block in many sensory and motor circuits. We studied the features that underly robustness in reciprocally inhibitory two neuron circuits. We used the dynamic clamp to create reciprocally inhibitory circuits from pharmacologically isolated neurons of the crab stomatogastric ganglion by injecting artificial graded synaptic (I_Syn) and hyperpolarization-activated inward (I_H) currents. There is a continuum of mechanisms in circuits that generate antiphase oscillations, with ‘release’ and ‘escape’ mechanisms at the extremes, and mixed mode oscillations between these extremes. In release, the active neuron primarily controls the off/on transitions. In escape, the inhibited neuron controls the transitions. We characterized the robustness of escape and release circuits to alterations in circuit parameters, temperature, and neuromodulation. We found that escape circuits rely on tight correlations between synaptic and H conductances to generate bursting but are resilient to temperature increase. Release circuits are robust to variations in synaptic and H conductances but fragile to temperature increase. The modulatory current (I_MI) restores oscillations in release circuits but has little effect in escape circuits. Perturbations can alter the balance of escape and release mechanisms and can create mixed mode oscillations. We conclude that the same perturbation can have dramatically different effects depending on the circuits’ mechanism of operation that may not be observable from basal circuit activity.

Degeneracy in the nervous system: from neuronal excitability to neural coding

Degeneracy is ubiquitous across biological systems where structurally different elements can yield a similar outcome. Degeneracy is of particular interest in neuroscience too. On the one hand, degeneracy confers robustness to the nervous system and facilitates evolvability: Different elements provide a backup plan for the system in response to any perturbation or disturbance. On the other, a difficulty in the treatment of some neurological disorders such as chronic pain is explained in light of different elements all of which contribute to the pathological behavior of the system. Under these circumstances, targeting a specific element is ineffective because other elements can compensate for this modulation. Understanding degeneracy in the physiological context explains its beneficial role in the robustness of neural circuits. Likewise, understanding degeneracy in the pathological context opens new avenues of discovery to find more effective therapies.

Ion Channel Degeneracy, Variability, and Covariation in Neuron and Circuit Resilience

The large number of ion channels found in all nervous systems poses fundamental questions concerning how the characteristic intrinsic properties of single neurons are determined by the specific subsets of channels they express. All neurons display many different ion channels with overlapping voltage- and time-dependent properties. We speculate that these overlapping properties promote resilience in neuronal function. Individual neurons of the same cell type show variability in ion channel conductance densities even though they can generate reliable and similar behavior. This complicates a simple assignment of function to any conductance and is associated with variable responses of neurons of the same cell type to perturbations, deletions, and pharmacological manipulation. Ion channel genes often show strong positively correlated expression, which may result from the molecular and developmental rules that determine which ion channels are expressed in a given cell type.

Circuit Stability to Perturbations Reveals Hidden Variability in the Balance of Intrinsic and Synaptic Conductances

Neurons and circuits each with a distinct balance of intrinsic and synaptic conductances can generate similar behavior but sometimes respond very differently to perturbation. Examining a large family of circuit models with non-identical neurons and synapses underlying rhythmic behavior, we analyzed the circuits' response to modifications in single and multiple intrinsic conductances in the individual neurons. To summarize these changes over the entire range of perturbed parameters, we quantified circuit output by defining a global stability measure. Using this measure, we identified specific subsets of conductances that when perturbed generate similar behavior in diverse individuals of the population. Our unbiased clustering analysis enabled us to quantify circuit stability when simultaneously perturbing multiple conductances as a nonlinear combination of single conductance perturbations. This revealed surprising conductance combinations that can predict the response to specific perturbations, even when the remaining intrinsic and synaptic conductances are unknown. Therefore, our approach can expose hidden variability in the balance of intrinsic and synaptic conductances of the same neurons across different versions of the same circuit solely from the circuit response to perturbations. Developed for a specific family of model circuits, our quantitative approach to characterizing high-dimensional degenerate systems provides a conceptual and analytic framework to guide future theoretical and experimental studies on degeneracy and robustness.

SIGNIFICANCE STATEMENT Neural circuits can generate nearly identical behavior despite neuronal and synaptic parameters varying several-fold between individual instantiations. Yet, when these parameters are perturbed through channel deletions and mutations or environmental disturbances, seemingly identical circuits can respond very differently. What distinguishes inconsequential perturbations that barely alter circuit behavior from disruptive perturbations that drastically disturb circuit output remains unclear. Focusing on a family of rhythmic circuits, we propose a computational approach to reveal hidden variability in the intrinsic and synaptic conductances in seemingly identical circuits based solely on circuit output to different perturbations. We uncover specific conductance combinations that work similarly to maintain stability and predict the effect of changing multiple conductances simultaneously, which often results from neuromodulation or injury.

Neuronal Homeostasis: Voltage Brings It All Together

Neurons generate consistent patterns of activity combining a large set of ionic currents and variable conductance levels. How they regulate this variability so that it does not run out of control seems to depend most prominently on neuronal activity itself.

Membrane Voltage Is a Direct Feedback Signal That Influences Correlated Ion Channel Expression in Neurons

The number and type of ion channels present in the membrane determines the electrophysiological function of every neuron. In many species, stereotyped output of neurons often persists for years [1], and ion channel dysregulation can change these properties to cause severe neurological disorders [2-4]. Thus, a fundamental question is how do neurons coordinate channel expression to uphold their firing patterns over long timescales [1, 5]? One major hypothesis purports that neurons homeostatically regulate their ongoing activity through mechanisms that link membrane voltage to expression relationships among ion channels [6-10]. However, experimentally establishing this feedback loop for the control of expression relationships has been a challenge: manipulations that aim to test for voltage feedback invariably disrupt trophic signaling from synaptic transmission and neuromodulation in addition to activity [9, 11, 12]. Further, neuronal activity often relies critically on these chemical mediators, obscuring the contribution of voltage activity of the membrane per se in forming the channel relationships that determine neuronal output [6, 13]. To resolve this, we isolated an identifiable neuron in crustaceans and then either kept this neuron silent or used a long-term voltage clamp protocol to artificially maintain activity. We found that physiological voltage activity-independent of all known forms of synaptic and neuromodulatory feedback-maintains most channel mRNA relationships, while metabotropic influences may play a relatively smaller role. Thus, ion channel relationships likely needed to maintain neuronal identity are actively and continually regulated at least in part at the level of channel mRNAs through feedback by membrane voltage.

Modulator-Gated, SUMOylation-Mediated, Activity-Dependent Regulation of Ionic Current Densities Contributes to Short-Term Activity Homeostasis

Neurons operate within defined activity limits, and feedback control mechanisms dynamically tune ionic currents to maintain this optimal range. This study describes a novel, rapid feedback mechanism that uses SUMOylation to continuously adjust ionic current densities according to changes in activity. Small ubiquitin-like modifier (SUMO) is a peptide that can be post-translationally conjugated to ion channels to influence their surface expression and biophysical properties. Neuronal activity can regulate the extent of protein SUMOylation. This study on the single, unambiguously identifiable lateral pyloric neuron (LP), a component of the pyloric network in the stomatogastric nervous system of male and female spiny lobsters (Panulirus interruptus), focused on dynamic SUMOylation in the context of activity homeostasis. There were four major findings: First, neuronal activity adjusted the balance between SUMO conjugation and deconjugation to continuously and bidirectionally fine-tune the densities of two opposing conductances: the hyperpolarization activated current (I_h) and the transient potassium current (I_A). Second, tonic 5 nm dopamine (DA) gated activity-dependent SUMOylation to permit and prevent activity-dependent regulation of I_h and I_A, respectively. Third, DA-gated, activity-dependent SUMOylation contributed to a feedback mechanism that restored the timing and duration of LP activity during prolonged modulation by 5 μm DA, which initially altered these and other activity features. Fourth, DA modulatory and metamoduatory (gating) effects were tailored to simultaneously alter and stabilize neuronal output. Our findings suggest that modulatory tone may select a subset of rapid activity-dependent mechanisms from a larger menu to achieve homeostasis under varying conditions.SIGNIFICANCE STATEMENT Post-translational SUMOylation of ion channel subunits controls their interactions. When subunit SUMOylation is dysregulated, conductance densities mediated by the channels are distorted, leading to nervous system disorders, such as seizures and chronic pain. Regulation of ion channel SUMOylation is poorly understood. This study demonstrated that neuronal activity can regulate SUMOylation to reconfigure ionic current densities over minutes, and this regulation was gated by tonic nanomolar dopamine. Dynamic SUMOylation was necessary to maintain specific aspects of neuronal output while the neuron was being modulated by high (5 μm) concentrations of dopamine, suggesting that the gating function may ensure neuronal homeostasis during extrinsic modulation of a circuit.

The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow

Every neuron is part of a network, exerting its function by transforming multiple spatiotemporal synaptic input patterns into a single spiking output. This function is specified by the particular shape and passive electrical properties of the neuronal membrane, and the composition and spatial distribution of ion channels across its processes. For a variety of physiological or pathological reasons, the intrinsic input/output function may change during a neuron's lifetime. This process results in high variability in the peak specific conductance of ion channels in individual neurons. The mechanisms responsible for this variability are not well understood, although there are clear indications from experiments and modeling that degeneracy and correlation among multiple channels may be involved. Here, we studied this issue in biophysical models of hippocampal CA1 pyramidal neurons and interneurons. Using a unified data-driven simulation workflow and starting from a set of experimental recordings and morphological reconstructions obtained from rats, we built and analyzed several ensembles of morphologically and biophysically accurate single cell models with intrinsic electrophysiological properties consistent with experimental findings. The results suggest that the set of conductances expressed in any given hippocampal neuron may be considered as belonging to two groups: one subset is responsible for the major characteristics of the firing behavior in each population and the other is responsible for a robust degeneracy. Analysis of the model neurons suggests several experimentally testable predictions related to the combination and relative proportion of the different conductances that should be expressed on the membrane of different types of neurons for them to fulfill their role in the hippocampus circuitry.

Origin of heterogeneous spiking patterns from continuously distributed ion channel densities: a computational study in spinal dorsal horn neurons

KEY POINTS

Distinct spiking patterns may arise from qualitative differences in ion channel expression (i.e. when different neurons express distinct ion channels) and/or when quantitative differences in expression levels qualitatively alter the spike generation process. We hypothesized that spiking patterns in neurons of the superficial dorsal horn (SDH) of spinal cord reflect both mechanisms. We reproduced SDH neuron spiking patterns by varying densities of KV 1- and A-type potassium conductances. Plotting the spiking patterns that emerge from different density combinations revealed spiking-pattern regions separated by boundaries (bifurcations). This map suggests that certain spiking pattern combinations occur when the distribution of potassium channel densities straddle boundaries, whereas other spiking patterns reflect distinct patterns of ion channel expression. The former mechanism may explain why certain spiking patterns co-occur in genetically identified neuron types. We also present algorithms to predict spiking pattern proportions from ion channel density distributions, and vice versa.

ABSTRACT

Neurons are often classified by spiking pattern. Yet, some neurons exhibit distinct patterns under subtly different test conditions, which suggests that they operate near an abrupt transition, or bifurcation. A set of such neurons may exhibit heterogeneous spiking patterns not because of qualitative differences in which ion channels they express, but rather because quantitative differences in expression levels cause neurons to operate on opposite sides of a bifurcation. Neurons in the spinal dorsal horn, for example, respond to somatic current injection with patterns that include tonic, single, gap, delayed and reluctant spiking. It is unclear whether these patterns reflect five cell populations (defined by distinct ion channel expression patterns), heterogeneity within a single population, or some combination thereof. We reproduced all five spiking patterns in a computational model by varying the densities of a low-threshold (KV 1-type) potassium conductance and an inactivating (A-type) potassium conductance and found that single, gap, delayed and reluctant spiking arise when the joint probability distribution of those channel densities spans two intersecting bifurcations that divide the parameter space into quadrants, each associated with a different spiking pattern. Tonic spiking likely arises from a separate distribution of potassium channel densities. These results argue in favour of two cell populations, one characterized by tonic spiking and the other by heterogeneous spiking patterns. We present algorithms to predict spiking pattern proportions based on ion channel density distributions and, conversely, to estimate ion channel density distributions based on spiking pattern proportions. The implications for classifying cells based on spiking pattern are discussed.

Mechanisms of generation of membrane potential resonance in a neuron with multiple resonant ionic currents

Neuronal membrane potential resonance (MPR) is associated with subthreshold and network oscillations. A number of voltage-gated ionic currents can contribute to the generation or amplification of MPR, but how the interaction of these currents with linear currents contributes to MPR is not well understood. We explored this in the pacemaker PD neurons of the crab pyloric network. The PD neuron MPR is sensitive to blockers of H- (I_H) and calcium-currents (I_Ca). We used the impedance profile of the biological PD neuron, measured in voltage clamp, to constrain parameter values of a conductance-based model using a genetic algorithm and obtained many optimal parameter combinations. Unlike most cases of MPR, in these optimal models, the values of resonant- (f_res) and phasonant- (f_{ϕ = 0}) frequencies were almost identical. Taking advantage of this fact, we linked the peak phase of ionic currents to their amplitude, in order to provide a mechanistic explanation the dependence of MPR on the I_Ca gating variable time constants. Additionally, we found that distinct pairwise correlations between I_Ca parameters contributed to the maintenance of f_res and resonance power (Q_Z). Measurements of the PD neuron MPR at more hyperpolarized voltages resulted in a reduction of f_res but no change in Q_Z. Constraining the optimal models using these data unmasked a positive correlation between the maximal conductances of I_H and I_Ca. Thus, although I_H is not necessary for MPR in this neuron type, it contributes indirectly by constraining the parameters of I_Ca.

Many neuron types exhibit membrane potential resonance (MPR) in which the neuron produces the largest response to oscillatory input at some preferred (resonant) frequency and, in many systems, the network frequency is correlated with neuronal MPR. MPR is captured by a peak in the impedance vs. frequency curve (Z-profile), which is shaped by the dynamics of voltage-gated ionic currents. Although neuron types can express variable levels of ionic currents, they may have a stable resonant frequency. We used the PD neuron of the crab pyloric network to understand how MPR emerges from the interplay of the biophysical properties of multiple ionic currents, each capable of generating resonance. We show the contribution of an inactivating current at the resonant frequency in terms of interacting time constants. We measured the Z-profile of the PD neuron and explored possible combinations of model parameters that fit this experimentally measured profile. We found that the Z-profile constrains and defines correlations among parameters associated with ionic currents. Furthermore, the resonant frequency and amplitude are sensitive to different parameter sets and can be preserved by co-varying pairs of parameters along their correlation lines. Furthermore, although a resonant current may be present in a neuron, it may not directly contribute to MPR, but constrain the properties of other currents that generate MPR. Finally, constraining model parameters further to those that modify their MPR properties to changes in voltage range produces maximal conductance correlations.

A balance of outward and linear inward ionic currents is required for generation of slow-wave oscillations

Regenerative inward currents help produce slow oscillations through a negative-slope conductance region of their current-voltage relationship that is well approximated by a linear negative conductance. We used dynamic-clamp injections of a linear current with such conductance, I_NL, to explore why some neurons can generate intrinsic slow oscillations whereas others cannot. We addressed this question in synaptically isolated neurons of the crab Cancer borealis after blocking action potentials. The pyloric network consists of a distinct pacemaker and follower neurons, all of which express the same complement of ionic currents. When the pyloric dilator (PD) neuron, a member of the pacemaker group, was injected with I_NL with dynamic clamp, it consistently produced slow oscillations. In contrast, all follower neurons failed to oscillate with I_NL To understand these distinct behaviors, we compared outward current levels of PD with those of follower lateral pyloric (LP) and ventral pyloric (VD) neurons. We found that LP and VD neurons had significantly larger high-threshold potassium currents (I_HTK) than PD and LP had lower-transient potassium current (I_A). Reducing I_HTK pharmacologically enabled both LP and VD neurons to produce I_NL-induced oscillations, whereas modifying I_A levels did not affect I_NL-induced oscillations. Using phase-plane and bifurcation analysis of a simplified model cell, we demonstrate that large levels of I_HTK can block I_NL-induced oscillatory activity whereas generation of oscillations is almost independent of I_A levels. These results demonstrate the general importance of a balance between inward pacemaking currents and high-threshold K⁺ current levels in determining slow oscillatory activity.NEW & NOTEWORTHY Pacemaker neuron-generated rhythmic activity requires the activation of at least one inward and one outward current. We have previously shown that the inward current can be a linear current (with negative conductance). Using this simple mechanism, here we demonstrate that the inward current conductance must be in relative balance with the outward current conductances to generate oscillatory activity. Surprisingly, an excess of outward conductances completely precludes the possibility of achieving such a balance.

Activation mechanism of a neuromodulator-gated pacemaker ionic current

The neuromodulator-gated current (I_MI) found in the crab stomatogastric ganglion is activated by neuromodulators that are essential to induce the rhythmic activity of the pyloric network in this system. One of these neuromodulators is also known to control the correlated expression of voltage-gated ionic currents in pyloric neurons, as well as synaptic plasticity and strength. Thus understanding the mechanism by which neuromodulator receptors activate I_MI should provide insights not only into how oscillations are initiated but also into how other processes, and currents not directly activated by them, are regulated. To determine what specific signaling molecules are implicated in this process, we used a battery of agonists and antagonists of common signal transduction pathways. We found that the G protein inhibitor GDPβS and the G protein activator GTPγS significantly affect I_MI amplitude, suggesting that its activation is mediated by G proteins. Interestingly, when using the more specific G protein blocker pertussis toxin, we observed the expected inhibition of I_MI amplitude but, unexpectedly, in a calcium-dependent fashion. We also found that antagonists of calcium- and calmodulin-associated signaling significantly reduce I_MI amplitude. In contrast, we found little evidence for the role of cyclic nucleotide signaling, phospholipase C (PLC), or kinases and phosphatases, except two calmodulin-dependent kinases. In sum, these results suggest that proctolin-induced I_MI is mediated by a G protein whose pertussis toxin sensitivity is altered by external calcium concentration and appears to depend on intracellular calcium, calmodulin, and calmodulin-activated kinases. In contrast, we found no support for I_MI being mediated by PLC signaling or cyclic nucleotides.NEW & NOTEWORTHY Neuronal rhythmic activity is generated by either network-based or cell-autonomous mechanisms. In the pyloric network of decapod crustaceans, the activation of a neuromodulator-gated pacemaker current is crucial for the generation of rhythmic activity. This current is activated by several neuromodulators, including peptides and acetylcholine, presumably via metabotropic receptors. We have previously demonstrated a novel extracellular calcium-sensitive voltage-dependence mechanism of this current. We presently report that the activation mechanism depends on intracellular and extracellular calcium-sensitive components.

Genetic perturbations suggest a role of the resting potential in regulating the expression of the ion channels of the KCNA and HCN families in octopus cells of the ventral cochlear nucleus

Low-voltage-activated K⁺ (g_KL) and hyperpolarization-activated mixed cation conductances (g_h) mediate currents, I_KL and I_h, through channels of the Kv1 (KCNA) and HCN families respectively and give auditory neurons the temporal precision required for signaling information about the onset, fine structure, and time of arrival of sounds. Being partially activated at rest, g_KL and g_h contribute to the resting potential and shape responses to even small subthreshold synaptic currents. Resting g_KL and g_h also affect the coupling of somatic depolarization with the generation of action potentials. To learn how these important conductances are regulated we have investigated how genetic perturbations affect their expression in octopus cells of the ventral cochlear nucleus (VCN). We report five new findings: First, the magnitude of g_h and g_KL varied over more than two-fold between wild type strains of mice. Second, average resting potentials are not different in different strains of mice even in the face of large differences in average g_KL and g_h. Third, I_KL has two components, one being α-dendrotoxin (α-DTX)-sensitive and partially inactivating and the other being α-DTX-insensitive, tetraethylammonium (TEA)-sensitive, and non-inactivating. Fourth, the loss of Kv1.1 results in diminution of the α-DTX-sensitive I_KL, and compensatory increased expression of an α-DTX-insensitive, tetraethylammonium (TEA)-sensitive I_KL. Fifth, I_h and I_KL are balanced at the resting potential in all wild type and mutant octopus cells even when resting potentials vary in individual cells over nearly 10 mV, indicating that the resting potential influences the expression of g_h and g_KL. The independence of resting potentials on g_KL and g_h shows that g_KL and g_h do not, over days or weeks, determine the resting potential but rather that the resting potential plays a role in regulating the magnitude of either or both g_KL and g_h.

Consequences of degeneracy in network function

Often distinct elements serve similar functions within a network. However, it is unclear whether this network degeneracy is beneficial, or merely a reflection of tighter regulation of overall network performance relative to individual neuronal properties. We review circumstances where data strongly suggest that degeneracy is beneficial in that it makes network function more robust. Importantly, network degeneracy is likely to have functional consequences that are not widely appreciated. This is likely to be true when network activity is configured by modulators with persistent actions, and the history of network activity potentially impacts subsequent functioning. Data suggest that degeneracy in this context may be important for the creation of latent memories, and for state-dependent task switching.

Dynamic compensation mechanism gives rise to period and duty-cycle level sets in oscillatory neuronal models

Rhythmic oscillation in neurons can be characterized by various attributes, such as the oscillation period and duty cycle. The values of these features depend on the amplitudes of the participating ionic currents, which can be characterized by their maximum conductance values. Recent experimental and theoretical work has shown that the values of these attributes can be maintained constant for different combinations of two or more ionic currents of varying conductances, defining what is known as level sets in conductance space. In two-dimensional conductance spaces, a level set is a curve, often a line, along which a particular oscillation attribute value is conserved. In this work, we use modeling, dynamical systems tools (phase-space analysis), and numerical simulations to investigate the possible dynamic mechanisms responsible for the generation of period and duty-cycle levels sets in simplified (linearized and FitzHugh-Nagumo) and conductance-based (Morris-Lecar) models of neuronal oscillations. A simplistic hypothesis would be that the tonic balance between ionic currents with the same or opposite effective signs is sufficient to create level sets. According to this hypothesis, the dynamics of each ionic current during a given cycle are well captured by some constant quantity (e.g., maximal conductances), and the phase-plane diagrams are identical or are almost identical (e.g., cubic-like nullclines with the same maxima and minima) for different combinations of these maximal conductances. In contrast, we show that these mechanisms are dynamic and involve the complex interaction between the nonlinear voltage dependencies and the effective time scales at which the ionic current's dynamical variables operate.

The complexity of small circuits: the stomatogastric nervous system

The crustacean stomatogastric nervous system is a long-standing test bed for studies of circuit dynamics and neuromodulation. We give a brief update on the most recent work on this system, with an emphasis on the broader implications for understanding neural circuits. In particular, we focus on new findings underlining that different levels of dynamics taking place at different time scales all interact in multiple ways. Dynamics due to synaptic and intrinsic neuronal properties, neuromodulation, and long-term gene expression-dependent regulation are not independent, but influence each other. Extensive research on the stomatogastric system shows that these dynamic interactions convey robustness to circuit operation, while facilitating the flexibility of producing multiple circuit outputs.

Computational implications of biophysical diversity and multiple timescales in neurons and synapses for circuit performance

Despite advances in experimental and theoretical neuroscience, we are still trying to identify key biophysical details that are important for characterizing the operation of brain circuits. Biological mechanisms at the level of single neurons and synapses can be combined as 'building blocks' to generate circuit function. We focus on the importance of capturing multiple timescales when describing these intrinsic and synaptic components. Whether inherent in the ionic currents, the neuron's complex morphology, or the neurotransmitter composition of synapses, these multiple timescales prove crucial for capturing the variability and richness of circuit output and enhancing the information-carrying capacity observed across nervous systems.

Monoaminergic tone supports conductance correlations and stabilizes activity features in pattern generating neurons of the lobster, Panulirus interruptus

Experimental and computational studies demonstrate that different sets of intrinsic and synaptic conductances can give rise to equivalent activity patterns. This is because the balance of conductances, not their absolute values, defines a given activity feature. Activity-dependent feedback mechanisms maintain neuronal conductance correlations and their corresponding activity features. This study demonstrates that tonic nM concentrations of monoamines enable slow, activity-dependent processes that can maintain a correlation between the transient potassium current (I_A) and the hyperpolarization activated current (I_h) over the long-term (i.e., regulatory change persists for hours after removal of modulator). Tonic 5 nM DA acted through an RNA interference silencing complex (RISC)- and RNA polymerase II-dependent mechanism to maintain a long-term positive correlation between I_A and I_h in the lateral pyloric neuron (LP) but not in the pyloric dilator neuron (PD). In contrast, tonic 5 nM 5HT maintained a RISC-dependent positive correlation between I_A and I_h in PD but not LP over the long-term. Tonic 5 nM OCT maintained a long-term negative correlation between I_A and I_h in PD but not LP; however, it was only revealed when RISC was inhibited. This study also demonstrated that monoaminergic tone can also preserve activity features over the long-term: the timing of LP activity, LP duty cycle and LP spike number per burst were maintained by tonic 5 nM DA. The data suggest that low-level monoaminergic tone acts through multiple slow processes to permit cell-specific, activity-dependent regulation of ionic conductances to maintain conductance correlations and their corresponding activity features over the long-term.

Feedback control of variability in the cycle period of a central pattern generator

We address how feedback to a bursting biological pacemaker with intrinsic variability in cycle length can affect that variability. Specifically, we examine a hybrid circuit constructed of an isolated crab anterior burster (AB)/pyloric dilator (PD) pyloric pacemaker receiving virtual feedback via dynamic clamp. This virtual feedback generates artificial synaptic input to PD with timing determined by adjustable phase response dynamics that mimic average burst intervals generated by the lateral pyloric neuron (LP) in the intact pyloric network. Using this system, we measure network period variability dependence on the feedback element's phase response dynamics and find that a constant response interval confers minimum variability. We further find that these optimal dynamics are characteristic of the biological pyloric network. Building upon our previous theoretical work mapping the firing intervals in one cycle onto the firing intervals in the next cycle, we create a theoretical map of the distribution of all firing intervals in one cycle to the distribution of firing intervals in the next cycle. We then obtain an integral equation for a stationary self-consistent distribution of the network periods of the hybrid circuit, which can be solved numerically given the uncoupled pacemaker's distribution of intrinsic periods, the nature of the network's feedback, and the phase resetting characteristics of the pacemaker. The stationary distributions obtained in this manner are strongly predictive of the experimentally observed distributions of hybrid network period. This theoretical framework can provide insight into optimal feedback schemes for minimizing variability to increase reliability or maximizing variability to increase flexibility in central pattern generators driven by pacemakers with feedback.

Dynamic Input Conductances Shape Neuronal Spiking

Reliable neuron activity is ensured by a tight regulation of the ion channels that resides in the neuron’s membrane. Understanding the causal mechanisms that relate this regulation to physiological and pathological neuronal activity is a necessary step for developing efficient therapies for neurological diseases associated with abnormal nervous system activity.

Assessing the role of biophysical parameter variations in neuronal activity is critical to the understanding of modulation, robustness, and homeostasis of neuronal signalling. The paper proposes that this question can be addressed through the analysis of dynamic input conductances. Those voltage-dependent curves aggregate the concomitant activity of all ion channels in distinct timescales. They are shown to shape the current−voltage dynamical relationships that determine neuronal spiking. We propose an experimental protocol to measure dynamic input conductances in neurons. In addition, we provide a computational method to extract dynamic input conductances from arbitrary conductance-based models and to analyze their sensitivity to arbitrary parameters. We illustrate the relevance of the proposed approach for modulation, compensation, and robustness studies in a published neuron model based on data of the stomatogastric ganglion of the crab Cancer borealis.

Tonic 5nM DA stabilizes neuronal output by enabling bidirectional activity-dependent regulation of the hyperpolarization activated current via PKA and calcineurin

Volume transmission results in phasic and tonic modulatory signals. The actions of tonic dopamine (DA) at type 1 DA receptors (D1Rs) are largely undefined. Here we show that tonic 5nM DA acts at D1Rs to stabilize neuronal output over minutes by enabling activity-dependent regulation of the hyperpolarization activated current (I _h). In the presence but not absence of 5nM DA, I _h maximal conductance (G _max) was adjusted according to changes in slow wave activity in order to maintain spike timing. Our study on the lateral pyloric neuron (LP), which undergoes rhythmic oscillations in membrane potential with depolarized plateaus, demonstrated that incremental, bi-directional changes in plateau duration produced corresponding alterations in LP I _hG _max when preparations were superfused with saline containing 5nM DA. However, when preparations were superfused with saline alone there was no linear correlation between LP I _hG_max and duty cycle. Thus, tonic nM DA modulated the capacity for activity to modulate LP I _h G _max; this exemplifies metamodulation (modulation of modulation). Pretreatment with the Ca2+-chelator, BAPTA, or the specific PKA inhibitor, PKI, prevented all changes in LP I _h in 5nM DA. Calcineurin inhibitors blocked activity-dependent changes enabled by DA and revealed a PKA-mediated, activity-independent enhancement of LP I _hG _max. These data suggested that tonic 5nM DA produced two simultaneous, PKA-dependent effects: a direct increase in LP I _h G _max and a priming event that permitted calcineurin regulation of LP I _h. The latter produced graded reductions in LP I _hG _max with increasing duty cycles. We also demonstrated that this metamodulation preserved the timing of LP’s first spike when network output was perturbed with bath-applied 4AP. In sum, 5nM DA permits slow wave activity to provide feedback that maintains spike timing, suggesting that one function of low-level, tonic modulation is to stabilize specific features of a dynamic output.

Robust circuit rhythms in small circuits arise from variable circuit components and mechanisms

Small central pattern generating circuits found in invertebrates have significant advantages for the study of the circuit mechanisms that generate brain rhythms. Experimental and computational studies of small oscillatory circuits reveal that similar rhythms can arise from disparate mechanisms. Animal-to-animal variation in the properties of single neurons and synapses may underly robust circuit performance, and can be revealed by perturbations. Neuromodulation can produce altered circuit performance but also ensure reliable circuit function.

Cell types, network homeostasis, and pathological compensation from a biologically plausible ion channel expression model

How do neurons develop, control, and maintain their electrical signaling properties in spite of ongoing protein turnover and perturbations to activity? From generic assumptions about the molecular biology underlying channel expression, we derive a simple model and show how it encodes an "activity set point" in single neurons. The model generates diverse self-regulating cell types and relates correlations in conductance expression observed in vivo to underlying channel expression rates. Synaptic as well as intrinsic conductances can be regulated to make a self-assembling central pattern generator network; thus, network-level homeostasis can emerge from cell-autonomous regulation rules. Finally, we demonstrate that the outcome of homeostatic regulation depends on the complement of ion channels expressed in cells: in some cases, loss of specific ion channels can be compensated; in others, the homeostatic mechanism itself causes pathological loss of function.

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.

Ionic Current Variability and Functional Stability in the Nervous System

Identified neurons in different animals express ionic currents at highly variable levels (population variability). If neuronal identity is associated with stereotypical function, as is the case in genetically identical neurons or in unambiguously identified individual neurons, this variability poses a conundrum: How is activity the same if the components that generate it-ionic current levels-are different? In some cases, ionic current variability across similar neurons generates an output gradient. However, many neurons produce very similar output activity, despite substantial variability in ionic conductances. It appears that, in many such cells, conductance levels of one ionic current vary in proportion to the conductance levels of another current. As a result, in a population of neurons, these conductances appear to be correlated. Here, I review theoretical and experimental work that suggests that neuronal ionic current correlation can reduce the global ionic current variability and can contribute to functional stability.

Criticality and degeneracy in injury-induced changes in primary afferent excitability and the implications for neuropathic pain

Neuropathic pain remains notoriously difficult to treat despite numerous drug targets. Here, we offer a novel explanation for this intractability. Computer simulations predicted that qualitative changes in primary afferent excitability linked to neuropathic pain arise through a switch in spike initiation dynamics when molecular pathologies reach a tipping point (criticality), and that this tipping point can be reached via several different molecular pathologies (degeneracy). We experimentally tested these predictions by pharmacologically blocking native conductances and/or electrophysiologically inserting virtual conductances. Multiple different manipulations successfully reproduced or reversed neuropathic changes in primary afferents from naïve or nerve-injured rats, respectively, thus confirming the predicted criticality and its degenerate basis. Degeneracy means that several different molecular pathologies are individually sufficient to cause hyperexcitability, and because several such pathologies co-occur after nerve injury, that no single pathology is uniquely necessary. Consequently, single-target-drugs can be circumvented by maladaptive plasticity in any one of several ion channels.

DOI:http://dx.doi.org/10.7554/eLife.02370.001

Although the pain associated with an injury is unpleasant, it normally serves an important purpose: to make you avoid its source. However, some pain appears to arise from nowhere. Frustratingly, this type of pain, known as neuropathic pain, does not respond to common painkillers and is thus very difficult to treat.

The neurons that transmit pain and other sensory information do so using electrical signals. In response to a stimulus, ions travel through channels in the membrane of a neuron, which leads to a change in the electrical potential of the membrane. When this change is large enough, a voltage spike is produced: this signal is ultimately transmitted to the brain.

When certain neurons fire too easily or too often, neuropathic pain can arise. This hyperexcitability can make something painful feel even worse, or it can make things hurt that shouldn’t. To prevent this, extensive research has been devoted to identify drugs that target particular types of ion channels and block them. However, despite the discovery of many promising drugs, those drugs have been frustratingly ineffective in clinical trials.

Using simulations and experiments, Ratté et al. have examined the behavior of a type of neuron that normally conducts information about touch, but the brain sometimes misinterprets this information as pain. Increasing the flow of ions through the cell membrane in these simulations eventually causes a ‘tipping point’ to be crossed, which triggers a dramatic, discontinuous change in spiking pattern. However, as several different types of ion channels contribute to the current, there are several different ways in which the tipping point can be crossed.

This ability to produce the same result by multiple means is a common feature of complex systems. Known as degeneracy, it makes systems more robust, as a given result can still be achieved if one particular attempt to achieve this result fails. The work of Ratté et al. helps to explain why drugs that target just one type of ion channel may fail to relieve neuropathic pain: maladaptive changes in any one of several other ion channels may circumvent the therapeutic effect.

DOI:http://dx.doi.org/10.7554/eLife.02370.002

The integrative role of the sigh in psychology, physiology, pathology, and neurobiology

"Sighs, tears, grief, distress" expresses Johann Sebastian Bach in a musical example for the relationship between sighs and deep emotions. This review explores the neurobiological basis of the sigh and its relationship with psychology, physiology, and pathology. Sighs monitor changes in brain states, induce arousal, and reset breathing variability. These behavioral roles homeostatically regulate breathing stability under physiological and pathological conditions. Sighs evoked in hypoxia evoke arousal and thereby become critical for survival. Hypoarousal and failure to sigh have been associated with sudden infant death syndrome. Increased breathing irregularity may provoke excessive sighing and hyperarousal, a behavioral sequence that may play a role in panic disorders. Essential for generating sighs and breathing is the pre-Bötzinger complex. Modulatory and synaptic interactions within this local network and between networks located in the brainstem, cerebellum, cortex, hypothalamus, amygdala, and the periaqueductal gray may govern the relationships between physiology, psychology, and pathology. Unraveling these circuits will lead to a better understanding of how we balance emotions and how emotions become pathological.

Correlated conductance parameters in leech heart motor neurons contribute to motor pattern formation

Neurons can have widely differing intrinsic membrane properties, in particular the density of specific conductances, but how these contribute to characteristic neuronal activity or pattern formation is not well understood. To explore the relationship between conductances, and in particular how they influence the activity of motor neurons in the well characterized leech heartbeat system, we developed a new multi-compartmental Hodgkin-Huxley style leech heart motor neuron model. To do so, we evolved a population of model instances, which differed in the density of specific conductances, capable of achieving specific output activity targets given an associated input pattern. We then examined the sensitivity of measures of output activity to conductances and how the model instances responded to hyperpolarizing current injections. We found that the strengths of many conductances, including those with differing dynamics, had strong partial correlations and that these relationships appeared to be linked by their influence on heart motor neuron activity. Conductances that had positive correlations opposed one another and had the opposite effects on activity metrics when perturbed whereas conductances that had negative correlations could compensate for one another and had similar effects on activity metrics.