Global structure, robustness, and modulation of neuronal models

Degenerate neuronal and circuit mechanisms important for generating rhythmic motor patterns

In 1996, we published a review article (Marder E, Calabrese RL. Physiol Rev 76: 687-717, 1996) describing the state of knowledge about the structure and function of the central pattern-generating circuits important for producing rhythmic behaviors. Although many of the core questions persist, much has changed since 1996. Here, we focus on newer studies that reveal ambiguities that complicate understanding circuit dynamics, despite the enormous technical advances of the recent past. In particular, we highlight recent studies of animal-to-animal variability and our understanding that circuit rhythmicity may be supported by multiple state-dependent mechanisms within the same animal and that robustness and resilience in the face of perturbation may depend critically on the presence of modulators and degenerate circuit mechanisms. Additionally, we highlight the use of computational models to ask whether there are generalizable principles about circuit motifs that can be found across rhythmic motor systems in different animal species.

The Brain's Best Kept Secret Is Its Degenerate Structure

Degeneracy is defined as multiple sets of solutions that can produce very similar system performance. Degeneracy is seen across phylogenetic scales, in all kinds of organisms. In neuroscience, degeneracy can be seen in the constellation of biophysical properties that produce a neuron's characteristic intrinsic properties and/or the constellation of mechanisms that determine circuit outputs or behavior. Here, we present examples of degeneracy at multiple levels of organization, from single-cell behavior, small circuits, large circuits, and, in cognition, drawing conclusions from work ranging from bacteria to human cognition. Degeneracy allows the individual-to-individual variability within a population that creates potential for evolution.

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.

Milking a spherical cow: Toy models in neuroscience

There are many different kinds of models, and they play many different roles in the scientific endeavour. Neuroscience, and biology more generally, has understandably tended to emphasise empirical models that are grounded in data and make specific, experimentally testable predictions. Meanwhile, strongly idealised or 'toy' models have played a central role in the theoretical development of other sciences such as physics. In this paper, we examine the nature of toy models and their prospects in neuroscience.

Circuit function is more robust to changes in synaptic than intrinsic conductances

Circuit function results from both intrinsic conductances of network neurons and the synaptic conductances that connect them. In models of neural circuits, different combinations of maximal conductances can give rise to similar activity. We compared the robustness of a neural circuit to changes in their intrinsic versus synaptic conductances. To address this, we performed a sensitivity analysis on a population of conductance-based models of the pyloric network from the crustacean stomatogastric ganglion (STG). The model network consists of three neurons with nine currents: a sodium current (Na), three potassium currents (Kd, KCa, A-type), two calcium currents (CaS and CaT), a hyperpolarization-activated current (H), a non-voltage-gated leak current (leak), and a neuromodulatory current (MI). The model cells are connected by seven synapses of two types, glutamatergic and cholinergic. We produced one hundred models of the pyloric network that displayed similar activities with values of maximal conductances distributed over wide ranges. We evaluated the robustness of each model to changes in their maximal conductances. We found that individual models have different sensitivities to changes in their maximal conductances, both in their intrinsic and synaptic conductances. As expected the models become less robust as the extent of the changes increase. Despite quantitative differences in their robustness, we found that in all cases, the model networks are more sensitive to the perturbation of their intrinsic conductances than their synaptic conductances.

Ion-channel degeneracy and heterogeneities in the emergence of signature physiological characteristics of dentate gyrus granule cells

Complex systems are neither fully determined nor completely random. Biological complex systems, including single neurons, manifest intermediate regimes of randomness that recruit integration of specific combinations of functionally specialized subsystems. Such emergence of biological function provides the substrate for the expression of degeneracy, the ability of disparate combinations of subsystems to yield similar function. Here, we present evidence for the expression of degeneracy in morphologically realistic models of dentate gyrus granule cells (GCs) through functional integration of disparate ion-channel combinations. We performed a 45-parameter randomized search spanning 16 active and passive ion channels, each biophysically constrained by their gating kinetics and localization profiles, to search for valid GC models. Valid models were those that satisfied 17 sub- and suprathreshold cellular-scale electrophysiological measurements from rat GCs. A vast majority (>99%) of the 15,000 random models were not electrophysiologically valid, demonstrating that arbitrarily random ion-channel combinations would not yield GC functions. The 141 valid models (0.94% of 15,000) manifested heterogeneities in and cross-dependencies across local and propagating electrophysiological measurements, which matched with their respective biological counterparts. Importantly, these valid models were widespread throughout the parametric space and manifested weak cross-dependencies across different parameters. These observations together showed that GC physiology could neither be obtained by entirely random ion-channel combinations nor is there an entirely determined single parametric combination that satisfied all constraints. The complexity, the heterogeneities in measurement and parametric spaces, and degeneracy associated with GC physiology should be rigorously accounted for while assessing GCs and their robustness under physiological and pathological conditions.NEW & NOTEWORTHY A recent study from our laboratory had demonstrated pronounced heterogeneities in a set of 17 electrophysiological measurements obtained from a large population of rat hippocampal granule cells. Here, we demonstrate the manifestation of ion-channel degeneracy in a heterogeneous population of morphologically realistic conductance-based granule cell models that were validated against these measurements and their cross-dependencies. Our analyses show that single neurons are complex entities whose functions emerge through intricate interactions among several functionally specialized subsystems.

Similarity in evoked responses does not imply similarity in macroscopic network states

It is commonplace in neuroscience to assume that if two tasks activate the same brain areas in the same way, then they are recruiting the same underlying networks. Yet computational theory has shown that the same pattern of activity can emerge from many different underlying network representations. Here we evaluated whether similarity in activation necessarily implies similarity in network architecture by comparing region-wise activation patterns and functional correlation profiles from a large sample of healthy subjects (N = 242). Participants performed two executive control tasks known to recruit nearly identical brain areas, the color-word Stroop task and the Multi-Source Interference Task (MSIT). Using a measure of instantaneous functional correlations, based on edge time series, we estimated the task-related networks that differed between incongruent and congruent conditions. We found that the two tasks were much more different in their network profiles than in their evoked activity patterns at different analytical levels, as well as for a wide range of methodological pipelines. Our results reject the notion that having the same activation patterns means two tasks engage the same underlying representations, suggesting that task representations should be independently evaluated at both node and edge (connectivity) levels.

Interactions between circuit architecture and plasticity in a closed-loop cerebellar system

Determining the sites and directions of plasticity underlying changes in neural activity and behavior is critical for understanding mechanisms of learning. Identifying such plasticity from neural recording data can be challenging due to feedback pathways that impede reasoning about cause and effect. We studied interactions between feedback, neural activity, and plasticity in the context of a closed-loop motor learning task for which there is disagreement about the loci and directions of plasticity: vestibulo-ocular reflex learning. We constructed a set of circuit models that differed in the strength of their recurrent feedback, from no feedback to very strong feedback. Despite these differences, each model successfully fit a large set of neural and behavioral data. However, the patterns of plasticity predicted by the models fundamentally differed, with the direction of plasticity at a key site changing from depression to potentiation as feedback strength increased. Guided by our analysis, we suggest how such models can be experimentally disambiguated. Our results address a long-standing debate regarding cerebellum-dependent motor learning, suggesting a reconciliation in which learning-related changes in the strength of synaptic inputs to Purkinje cells are compatible with seemingly oppositely directed changes in Purkinje cell spiking activity. More broadly, these results demonstrate how changes in neural activity over learning can appear to contradict the sign of the underlying plasticity when either internal feedback or feedback through the environment is present.

Oscillatory network spontaneously recovers both activity and robustness after prolonged removal of neuromodulators

Robustness of neuronal activity is a property necessary for a neuronal network to withstand perturbations, which may otherwise disrupt or destroy the system. The robustness of complex systems has been shown to depend on a number of features of the system, including morphology and heterogeneity of the activity of the component neurons, size of the networks, synaptic connectivity, and neuromodulation. The activity of small networks, such as the pyloric network of the crustacean stomatogastric nervous system, appears to be robust despite some of the factors not being consistent with the expected properties of complex systems, e.g., small size and homogeneity of the synaptic connections. The activity of the pyloric network has been shown to be stable and robust in a neuromodulatory state-dependent manner. When neuromodulatory inputs are severed, activity is initially disrupted, losing both stability and robustness. Over the long term, however, stable activity homeostatically recovers without the restoration of neuromodulatory input. The question we address in this study is whether robustness can also be restored as the network reorganizes itself to compensate for the loss of neuromodulatory input and recovers the lost activity. Here, we use temperature changes as a perturbation to probe the robustness of the network's activity. We develop a simple metric of robustness, i.e., the variances of the network phase relationships, and show that robustness is indeed restored simultaneously along with its stable network activity, indicating that, whatever the reorganization of the network entails, it is deep enough also to restore this important property.

The enigmatic HCN channels: A cellular neurophysiology perspective

What physiological role does a slow hyperpolarization-activated ion channel with mixed cation selectivity play in the fast world of neuronal action potentials that are driven by depolarization? That puzzling question has piqued the curiosity of physiology enthusiasts about the hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, which are widely expressed across the body and especially in neurons. In this review, we emphasize the need to assess HCN channels from the perspective of how they respond to time-varying signals, while also accounting for their interactions with other co-expressing channels and receptors. First, we illustrate how the unique structural and functional characteristics of HCN channels allow them to mediate a slow negative feedback loop in the neurons that they express in. We present the several physiological implications of this negative feedback loop to neuronal response characteristics including neuronal gain, voltage sag and rebound, temporal summation, membrane potential resonance, inductive phase lead, spike triggered average, and coincidence detection. Next, we argue that the overall impact of HCN channels on neuronal physiology critically relies on their interactions with other co-expressing channels and receptors. Interactions with other channels allow HCN channels to mediate intrinsic oscillations, earning them the "pacemaker channel" moniker, and to regulate spike frequency adaptation, plateau potentials, neurotransmitter release from presynaptic terminals, and spike initiation at the axonal initial segment. We also explore the impact of spatially non-homogeneous subcellular distributions of HCN channels in different neuronal subtypes and their interactions with other channels and receptors. Finally, we discuss how plasticity in HCN channels is widely prevalent and can mediate different encoding, homeostatic, and neuroprotective functions in a neuron. In summary, we argue that HCN channels form an important class of channels that mediate a diversity of neuronal functions owing to their unique gating kinetics that made them a puzzle in the first place.

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.

A biophysical perspective on the resilience of neuronal excitability across timescales

Neuronal membrane excitability must be resilient to perturbations that can take place over timescales from milliseconds to months (or even years in long-lived animals). A great deal of attention has been paid to classes of homeostatic mechanisms that contribute to long-term maintenance of neuronal excitability through processes that alter a key structural parameter: the number of ion channel proteins present at the neuronal membrane. However, less attention has been paid to the self-regulating 'automatic' mechanisms that contribute to neuronal resilience by virtue of the kinetic properties of ion channels themselves. Here, we propose that these two sets of mechanisms are complementary instantiations of feedback control, together enabling resilience on a wide range of temporal scales. We further point to several methodological and conceptual challenges entailed in studying these processes - both of which involve enmeshed feedback control loops - and consider the consequences of these mechanisms of resilience.

Gating of homeostatic regulation of intrinsic excitability produces cryptic long-term storage of prior perturbations

Neurons and neuronal circuits must maintain their function throughout the life of the organism despite changing environments. Previous theoretical and experimental work suggests that neurons monitor their activity using intracellular calcium concentrations to regulate their intrinsic excitability. Models with multiple sensors can distinguish among different patterns of activity, but previous work using models with multiple sensors produced instabilities that lead the models' conductances to oscillate and then to grow without bound and diverge. We now introduce a nonlinear degradation term that explicitly prevents the maximal conductances to grow beyond a bound. We combine the sensors' signals into a master feedback signal that can be used to modulate the timescale of conductance evolution. Effectively, this means that the negative feedback can be gated on and off according to how far the neuron is from its target. The modified model recovers from multiple perturbations. Interestingly, depolarizing the models to the same membrane potential with current injection or with simulated high extracellular K+ produces different changes in conductances, arguing that caution must be used in interpreting manipulations that serve as a proxy for increased neuronal activity. Finally, these models accrue traces of prior perturbations that are not visible in their control activity after perturbation but that shape their responses to subsequent perturbations. These cryptic or hidden changes may provide insight into disorders such as posttraumatic stress disorder that only become visible in response to specific perturbations.

Larval zebrafish as a model for studying individual variability in translational neuroscience research

The larval zebrafish is a popular model for translational research into neurological and psychiatric disorders due to its conserved vertebrate brain structures, ease of genetic and experimental manipulation and small size and scalability to large numbers. The possibility of obtaining in vivo whole-brain cellular resolution neural data is contributing important advances into our understanding of neural circuit function and their relation to behavior. Here we argue that the larval zebrafish is ideally poised to push our understanding of how neural circuit function relates to behavior to the next level by including considerations of individual differences. Understanding variability across individuals is particularly relevant for tackling the variable presentations that neuropsychiatric conditions frequently show, and it is equally elemental if we are to achieve personalized medicine in the future. We provide a blueprint for investigating variability by covering examples from humans and other model organisms as well as existing examples from larval zebrafish. We highlight recent studies where variability may be hiding in plain sight and suggest how future studies can take advantage of existing paradigms for further exploring individual variability. We conclude with an outlook on how the field can harness the unique strengths of the zebrafish model to advance this important impending translational question.

Slow negative feedback enhances robustness of square-wave bursting

Square-wave bursting is an activity pattern common to a variety of neuronal and endocrine cell models that has been linked to central pattern generation for respiration and other physiological functions. Many of the reduced mathematical models that exhibit square-wave bursting yield transitions to an alternative pseudo-plateau bursting pattern with small parameter changes. This susceptibility to activity change could represent a problematic feature in settings where the release events triggered by spike production are necessary for function. In this work, we analyze how model bursting and other activity patterns vary with changes in a timescale associated with the conductance of a fast inward current. Specifically, using numerical simulations and dynamical systems methods, such as fast-slow decomposition and bifurcation and phase-plane analysis, we demonstrate and explain how the presence of a slow negative feedback associated with a gradual reduction of a fast inward current in these models helps to maintain the presence of spikes within the active phases of bursts. Therefore, although such a negative feedback is not necessary for burst production, we find that its presence generates a robustness that may be important for function.

Neural learning rules for generating flexible predictions and computing the successor representation

The predictive nature of the hippocampus is thought to be useful for memory-guided cognitive behaviors. Inspired by the reinforcement learning literature, this notion has been formalized as a predictive map called the successor representation (SR). The SR captures a number of observations about hippocampal activity. However, the algorithm does not provide a neural mechanism for how such representations arise. Here, we show the dynamics of a recurrent neural network naturally calculate the SR when the synaptic weights match the transition probability matrix. Interestingly, the predictive horizon can be flexibly modulated simply by changing the network gain. We derive simple, biologically plausible learning rules to learn the SR in a recurrent network. We test our model with realistic inputs and match hippocampal data recorded during random foraging. Taken together, our results suggest that the SR is more accessible in neural circuits than previously thought and can support a broad range of cognitive functions.

Neuronal morphology enhances robustness to perturbations of channel densities

Biological neurons show significant cell-to-cell variability but have the striking ability to maintain their key firing properties in the face of unpredictable perturbations and stochastic noise. Using a population of multi-compartment models consisting of soma, neurites, and axon for the lateral pyloric neuron in the crab stomatogastric ganglion, we explore how rebound bursting is preserved when the 14 channel conductances in each model are all randomly varied. The coupling between the axon and other compartments is critical for the ability of the axon to spike during bursts and consequently determines the set of successful solutions. When the coupling deviates from a biologically realistic range, the neuronal tolerance of conductance variations is lessened. Thus, the gross morphological features of these neurons enhance their robustness to perturbations of channel densities and expand the space of individual variability that can maintain a desired output pattern.

Machine learning provides insight into models of heterogeneous electrical activity in human beta-cells

Understanding how heterogeneous cellular responses emerge from cell-to-cell variations in expression and function of subcellular components is of general interest. Here, we focus on human insulin-secreting beta-cells, which are believed to constitute a population in which heterogeneity is of physiological importance. We exploit recent single-cell electrophysiological data that allow biologically realistic population modeling of human beta-cells that accounts for cellular heterogeneity and correlation between ion channel parameters. To investigate how ion channels influence the dynamics of our updated mathematical model of human pancreatic beta-cells, we explore several machine learning techniques to determine which model parameters are important for determining the qualitative patterns of electrical activity of the model cells. As expected, K^＋ channels promote absence of activity, but once a cell is active, they increase the likelihood of having action potential firing. HERG channels were of great importance for determining cell behavior in most of the investigated scenarios. Fast bursting is influenced by the time scales of ion channel activation and, interestingly, by the type of Ca^2＋ channels coupled to BK channels in BK-CaV complexes. Slow, metabolically driven oscillations are promoted mostly by K(ATP) channels. In summary, combining population modeling with machine learning analysis provides insight into the model and generates new hypotheses to be investigated both experimentally, via simulations and through mathematical analysis.

Energy-efficient network activity from disparate circuit parameters

Neural circuits can produce similar activity patterns from vastly different combinations of channel and synaptic conductances. These conductances are tuned for specific activity patterns but might also reflect additional constraints, such as metabolic cost or robustness to perturbations. How do such constraints influence the range of permissible conductances? Here we investigate how metabolic cost affects the parameters of neural circuits with similar activity in a model of the pyloric network of the crab Cancer borealis . We present a machine learning method that can identify a range of network models that generate activity patterns matching experimental data and find that neural circuits can consume largely different amounts of energy despite similar circuit activity. Furthermore, a reduced but still significant range of circuit parameters gives rise to energy-efficient circuits. We then examine the space of parameters of energy-efficient circuits and identify potential tuning strategies for low metabolic cost. Finally, we investigate the interaction between metabolic cost and temperature robustness. We show that metabolic cost can vary across temperatures but that robustness to temperature changes does not necessarily incur an increased metabolic cost. Our analyses show that despite metabolic efficiency and temperature robustness constraining circuit parameters, neural systems can generate functional, efficient, and robust network activity with widely disparate sets of conductances.

Voltage-Gated Ion Channels and the Variability in Information Transfer

The prerequisites for neurons to function within a circuit and be able to contain and transfer information efficiently and reliably are that they need to be homeostatically stable and fire within a reasonable range, characteristics that are governed, among others, by voltage-gated ion channels (VGICs). Nonetheless, neurons entail large variability in the expression levels of VGICs and their corresponding intrinsic properties, but the role of this variability in information transfer is not fully known. In this study, we aimed to investigate how this variability of VGICs affects information transfer. For this, we used a previously derived population of neuronal model neurons, each with the variable expression of five types of VGICs, fast Na+, delayed rectifier K+, A-type K+, T-type Ca++, and HCN channels. These analyses showed that the model neurons displayed variability in mutual information transfer, measured as the capability of neurons to successfully encode incoming synaptic information in output firing frequencies. Likewise, variability in the expression of VGICs caused variability in EPSPs and IPSPs amplitudes, reflected in the variability of output firing frequencies. Finally, using the virtual knockout methodology, we show that among the ion channels tested, the A-type K+ channel is the major regulator of information processing and transfer.

Scaling and Benchmarking an Evolutionary Algorithm for Constructing Biophysical Neuronal Models

Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.

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.

Predicting Spike Features of Hodgkin-Huxley-Type Neurons With Simple Artificial Neural Network

Hodgkin-Huxley (HH)-type model is the most famous computational model for simulating neural activity. It shows the highest accuracy in capturing neuronal spikes, and its model parameters have definite physiological meanings. However, HH-type models are computationally expensive. To address this problem, a previous study proposed a spike prediction module (SPM) to predict whether a spike will take place 1 ms later based on three voltage values with intervals of 1 ms. Although SPM does well, it fails to evaluate the informative features of the spike. In this study, the feature prediction module (FPM) based on simple artificial neural network (ANN) was proposed to predict spike features including maximum voltage, minimum voltage, and dropping interval. Nine different HH-type models were adopted whose firing patterns cover most of the firing behaviors observed in the brain. Voltage and spike feature samples under constant external input current were collected for training and testing. Experiment results illustrated that the combination of SPM and FPM can accurately predict the spiking part of different HH-type models and can generalize to unseen types of input current. The combination of SPM and FPM may offer a possible way to simulate the action potentials of biological neurons with high accuracy and efficiency.

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.

A transcriptional rheostat couples past activity to future sensory responses

Animals traversing different environments encounter both stable background stimuli and novel cues, which are thought to be detected by primary sensory neurons and then distinguished by downstream brain circuits. Here, we show that each of the ∼1,000 olfactory sensory neuron (OSN) subtypes in the mouse harbors a distinct transcriptome whose content is precisely determined by interactions between its odorant receptor and the environment. This transcriptional variation is systematically organized to support sensory adaptation: expression levels of more than 70 genes relevant to transforming odors into spikes continuously vary across OSN subtypes, dynamically adjust to new environments over hours, and accurately predict acute OSN-specific odor responses. The sensory periphery therefore separates salient signals from predictable background via a transcriptional rheostat whose moment-to-moment state reflects the past and constrains the future; these findings suggest a general model in which structured transcriptional variation within a cell type reflects individual experience.

Neuromodulation and Individuality

Within populations, individuals show a variety of behavioral preferences, even in the absence of genetic or environmental variability. Neuromodulators affect these idiosyncratic preferences in a wide range of systems, however, the mechanism(s) by which they do so is unclear. I review the evidence supporting three broad mechanisms by which neuromodulators might affect variability in idiosyncratic behavioral preference: by being a source of variability directly upstream of behavior, by affecting the behavioral output of a circuit in a way that masks or accentuates underlying variability in that circuit, and by driving plasticity in circuits leading to either homeostatic convergence toward a given behavior or divergence from a developmental setpoint. I find evidence for each of these mechanisms and propose future directions to further understand the complex interplay between individual variability and neuromodulators.

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.

Modeling spinal locomotor circuits for movements in developing zebrafish

Many spinal circuits dedicated to locomotor control have been identified in the developing zebrafish. How these circuits operate together to generate the various swimming movements during development remains to be clarified. In this study, we iteratively built models of developing zebrafish spinal circuits coupled to simplified musculoskeletal models that reproduce coiling and swimming movements. The neurons of the models were based upon morphologically or genetically identified populations in the developing zebrafish spinal cord. We simulated intact spinal circuits as well as circuits with silenced neurons or altered synaptic transmission to better understand the role of specific spinal neurons. Analysis of firing patterns and phase relationships helped to identify possible mechanisms underlying the locomotor movements of developing zebrafish. Notably, our simulations demonstrated how the site and the operation of rhythm generation could transition between coiling and swimming. The simulations also underlined the importance of contralateral excitation to multiple tail beats. They allowed us to estimate the sensitivity of spinal locomotor networks to motor command amplitude, synaptic weights, length of ascending and descending axons, and firing behavior. These models will serve as valuable tools to test and further understand the operation of spinal circuits for locomotion.

Does Differential Receptor Distribution Underlie Variable Responses to a Neuropeptide in the Lobster Cardiac System?

Central pattern generators produce rhythmic behaviors independently of sensory input; however, their outputs can be modulated by neuropeptides, thereby allowing for functional flexibility. We investigated the effects of C-type allatostatins (AST-C) on the cardiac ganglion (CG), which is the central pattern generator that controls the heart of the American lobster, Homarus americanus, to identify the biological mechanism underlying the significant variability in individual responses to AST-C. We proposed that the presence of multiple receptors, and thus differential receptor distribution, was at least partly responsible for this observed variability. Using transcriptome mining and PCR-based cloning, we identified four AST-C receptors (ASTCRs) in the CG; we then characterized their cellular localization, binding potential, and functional activation. Only two of the four receptors, ASTCR1 and ASTCR2, were fully functional GPCRs that targeted to the cell surface and were activated by AST-C peptides in our insect cell expression system. All four, however, were amplified from CG cDNAs. Following the confirmation of ASTCR expression, we used physiological and bioinformatic techniques to correlate receptor expression with cardiac responses to AST-C across individuals. Expression of ASTCR1 in the CG showed a negative correlation with increasing contraction amplitude in response to AST-C perfusion through the lobster heart, suggesting that the differential expression of ASTCRs within the CG is partly responsible for the specific physiological response to AST-C exhibited by a given individual lobster.

Interrogating theoretical models of neural computation with emergent property inference

A cornerstone of theoretical neuroscience is the circuit model: a system of equations that captures a hypothesized neural mechanism. Such models are valuable when they give rise to an experimentally observed phenomenon -- whether behavioral or a pattern of neural activity -- and thus can offer insights into neural computation. The operation of these circuits, like all models, critically depends on the choice of model parameters. A key step is then to identify the model parameters consistent with observed phenomena: to solve the inverse problem. In this work, we present a novel technique, emergent property inference (EPI), that brings the modern probabilistic modeling toolkit to theoretical neuroscience. When theorizing circuit models, theoreticians predominantly focus on reproducing computational properties rather than a particular dataset. Our method uses deep neural networks to learn parameter distributions with these computational properties. This methodology is introduced through a motivational example of parameter inference in the stomatogastric ganglion. EPI is then shown to allow precise control over the behavior of inferred parameters and to scale in parameter dimension better than alternative techniques. In the remainder of this work, we present novel theoretical findings in models of primary visual cortex and superior colliculus, which were gained through the examination of complex parametric structure captured by EPI. Beyond its scientific contribution, this work illustrates the variety of analyses possible once deep learning is harnessed towards solving theoretical inverse problems.

Modeling of sustained spontaneous network oscillations of a sexually dimorphic brainstem nucleus: the role of potassium equilibrium potential

Intrinsic oscillators in the central nervous system play a preeminent role in the neural control of rhythmic behaviors, yet little is known about how the ionic milieu regulates their output patterns. A powerful system to address this question is the pacemaker nucleus of the weakly electric fish Apteronotus leptorhynchus. A neural network comprised of an average of 87 pacemaker cells and 20 relay cells produces tonic oscillations, with higher frequencies in males compared to females. Previous empirical studies have suggested that this sexual dimorphism develops and is maintained through modulation of buffering of extracellular K⁺ by a massive meshwork of astrocytes enveloping the pacemaker and relay cells. Here, we constructed a model of this neural network that can generate sustained spontaneous oscillations. Sensitivity analysis revealed the potassium equilibrium potential, E_K (as a proxy of extracellular K⁺ concentration), and corresponding somatic channel conductances as critical determinants of oscillation frequency and amplitude. In models of both the pacemaker nucleus network and isolated pacemaker and relay cells, the frequency increased almost linearly with E_K, whereas the amplitude decreased nonlinearly with increasing E_K. Our simulations predict that this frequency increase is largely caused by a shift in the minimum K⁺ conductance over one oscillation period. This minimum is close to zero at more negative E_K, converging to the corresponding maximum at less negative E_K. This brings the resting membrane potential closer to the threshold potential at which voltage-gated Na⁺ channels become active, increasing the excitability, and thus the frequency, of pacemaker and relay cells.

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.

Ion-channel regulation of response decorrelation in a heterogeneous multi-scale model of the dentate gyrus

Heterogeneities in biological neural circuits manifest in afferent connectivity as well as in local-circuit components such as neuronal excitability, neural structure and local synaptic strengths. The expression of adult neurogenesis in the dentate gyrus (DG) amplifies local-circuit heterogeneities and guides heterogeneities in afferent connectivity. How do neurons and their networks endowed with these distinct forms of heterogeneities respond to perturbations to individual ion channels, which are known to change under several physiological and pathophysiological conditions? We sequentially traversed the ion channels-neurons-network scales and assessed the impact of eliminating individual ion channels on conductance-based neuronal and network models endowed with disparate local-circuit and afferent heterogeneities. We found that many ion channels differentially contributed to specific neuronal or network measurements, and the elimination of any given ion channel altered several functional measurements. We then quantified the impact of ion-channel elimination on response decorrelation, a well-established metric to assess the ability of neurons in a network to convey complementary information, in DG networks endowed with different forms of heterogeneities. Notably, we found that networks constructed with structurally immature neurons exhibited functional robustness, manifesting as minimal changes in response decorrelation in the face of ion-channel elimination. Importantly, the average change in output correlation was dependent on the eliminated ion channel but invariant to input correlation. Our analyses suggest that neurogenesis-driven structural heterogeneities could assist the DG network in providing functional resilience to molecular perturbations.

•

Perturbations at one scale result in a cascading impact on physiology across scales.

•

Heterogeneous multi-scale models used to assess the impact of ion-channel deletion.

•

Mapping of structural components to functional outcomes is many-to-many.

•

Differential & variable impact of ion channel deletion on response decorrelation.

•

Neurogenesis-induced structural heterogeneity confers resilience to perturbations.

Generalizable Machine Learning in Neuroscience Using Graph Neural Networks

Although a number of studies have explored deep learning in neuroscience, the application of these algorithms to neural systems on a microscopic scale, i.e. parameters relevant to lower scales of organization, remains relatively novel. Motivated by advances in whole-brain imaging, we examined the performance of deep learning models on microscopic neural dynamics and resulting emergent behaviors using calcium imaging data from the nematode C. elegans. As one of the only species for which neuron-level dynamics can be recorded, C. elegans serves as the ideal organism for designing and testing models bridging recent advances in deep learning and established concepts in neuroscience. We show that neural networks perform remarkably well on both neuron-level dynamics prediction and behavioral state classification. In addition, we compared the performance of structure agnostic neural networks and graph neural networks to investigate if graph structure can be exploited as a favourable inductive bias. To perform this experiment, we designed a graph neural network which explicitly infers relations between neurons from neural activity and leverages the inferred graph structure during computations. In our experiments, we found that graph neural networks generally outperformed structure agnostic models and excel in generalization on unseen organisms, implying a potential path to generalizable machine learning in neuroscience.

Intrinsic plasticity and birdsong learning

Although information processing and storage in the brain is thought to be primarily orchestrated by synaptic plasticity, other neural mechanisms such as intrinsic plasticity are available. While a number of recent studies have described the plasticity of intrinsic excitability in several types of neurons, the significance of non-synaptic mechanisms in memory and learning remains elusive. After reviewing plasticity of intrinsic excitation in relation to learning and homeostatic mechanisms, we focus on the intrinsic properties of a class of basal-ganglia projecting song system neurons in zebra finch, how these related to each bird's unique learned song, how these properties change over development, and how they are maintained dynamically to rapidly change in response to auditory feedback perturbations. We place these results in the broader theme of learning and changes in intrinsic properties, emphasizing the computational implications of this form of plasticity, which are distinct from synaptic plasticity. The results suggest that exploring reciprocal interactions between intrinsic and network properties will be a fruitful avenue for understanding mechanisms of birdsong learning.

The many bits of positional information

Half a century after Lewis Wolpert's seminal conceptual advance on how cellular fates distribute in space, we provide a brief historical perspective on how the concept of positional information emerged and influenced the field of developmental biology and beyond. We focus on a modern interpretation of this concept in terms of information theory, largely centered on its application to cell specification in the early Drosophila embryo. We argue that a true physical variable (position) is encoded in local concentrations of patterning molecules, that this mapping is stochastic, and that the processes by which positions and corresponding cell fates are determined based on these concentrations need to take such stochasticity into account. With this approach, we shift the focus from biological mechanisms, molecules, genes and pathways to quantitative systems-level questions: where does positional information reside, how it is transformed and accessed during development, and what fundamental limits it is subject to?

Lateral habenula dysfunctions in Tm4sf2-/y mice model for neurodevelopmental disorder

Mutations in the TM4SF2 gene, which encodes TSPAN7, cause a severe form of intellectual disability (ID) often comorbid with autism spectrum disorder (ASD). Recently, we found that TM4SF2 loss in mice affects cognition. Here, we report that Tm4sf2-/y mice, beyond an ID-like phenotype, display altered sociability, increased repetitive behaviors, anhedonic- and depressive-like states. Cognition relies on the integration of information from several brain areas. In this context, the lateral habenula (LHb) is strategically positioned to coordinate the brain regions involved in higher cognitive functions. Furthermore, in Tm4sf2-/y mice we found that LHb neurons present hypoexcitability, aberrant neuronal firing pattern and altered sodium and potassium voltage-gated ion channels function. Interestingly, we also found a reduced expression of voltage-gated sodium channel and a hyperactivity of the PKC-ERK pathway, a well-known modulator of ion channels activity, which might explain the functional phenotype showed by Tm4sf2-/y mice LHb neurons. These findings support Tm4sf2-/y mice as useful in modeling some ASD-like symptoms. Additionally, we can speculate that LHb functional alteration in Tm4sf2-/y mice might play a role in the disease pathophysiology.

Bayesian inference for biophysical neuron models enables stimulus optimization for retinal neuroprosthetics

While multicompartment models have long been used to study the biophysics of neurons, it is still challenging to infer the parameters of such models from data including uncertainty estimates. Here, we performed Bayesian inference for the parameters of detailed neuron models of a photoreceptor and an OFF- and an ON-cone bipolar cell from the mouse retina based on two-photon imaging data. We obtained multivariate posterior distributions specifying plausible parameter ranges consistent with the data and allowing to identify parameters poorly constrained by the data. To demonstrate the potential of such mechanistic data-driven neuron models, we created a simulation environment for external electrical stimulation of the retina and optimized stimulus waveforms to target OFF- and ON-cone bipolar cells, a current major problem of retinal neuroprosthetics.

An Understandable, Extensible, and Reusable Implementation of the Hodgkin-Huxley Equations Using Modelica

The Hodgkin-Huxley model of the squid giant axon has been used for decades as the basis of many action potential models. These models are usually communicated using just a list of equations or a circuit diagram, which makes them unnecessarily complicated both for novices and for experts. We present a modular version of the Hodgkin-Huxley model that is more understandable than the usual monolithic implementations and that can be easily reused and extended. Our model is written in Modelica using software engineering concepts, such as object orientation and inheritance. It retains the electrical analogy, but names and explains individual components in biological terms. We use cognitive load theory to measure understandability as the amount of items that have to be kept in working memory simultaneously. The model is broken down into small self-contained components in human-readable code with extensive documentation. Additionally, it features a hybrid diagram that uses biological symbols in an electrical circuit and that is directly tied to the model code. The new model design avoids many redundancies and reduces the cognitive load associated with understanding the model by a factor of 6. Extensions can be easily applied due to an unifying interface and inheritance from shared base classes. The model can be used in an educational context as a more approachable introduction to mathematical modeling in electrophysiology. Additionally the modeling approach and the base components can be used to make complex Hodgkin-Huxley-type models more understandable and reusable.

Activity-dependent compensation of cell size is vulnerable to targeted deletion of ion channels

In many species, excitable cells preserve their physiological properties despite significant variation in physical size across time and in a population. For example, neurons in crustacean central pattern generators generate similar firing patterns despite several-fold increases in size between juveniles and adults. This presents a biophysical problem because the electrical properties of cells are highly sensitive to membrane area and channel density. It is not known whether specific mechanisms exist to sense membrane area and adjust channel expression to keep a consistent channel density, or whether regulation mechanisms that sense activity alone are capable of compensating cell size. We show that destabilising effects of growth can be specifically compensated by feedback mechanism that senses average calcium influx and jointly regulate multiple conductances. However, we further show that this class of growth-compensating regulation schemes is necessarily sensitive to perturbations that alter the expression of subsets of ion channel types. Targeted perturbations of specific ion channels can trigger a pathological response of the regulation mechanism and a failure of homeostasis. Our findings suggest that physiological regulation mechanisms that confer robustness to growth may be specifically vulnerable to deletions or mutations that affect subsets of ion channels.

Training deep neural density estimators to identify mechanistic models of neural dynamics

Mechanistic modeling in neuroscience aims to explain observed phenomena in terms of underlying causes. However, determining which model parameters agree with complex and stochastic neural data presents a significant challenge. We address this challenge with a machine learning tool which uses deep neural density estimators-trained using model simulations-to carry out Bayesian inference and retrieve the full space of parameters compatible with raw data or selected data features. Our method is scalable in parameters and data features and can rapidly analyze new data after initial training. We demonstrate the power and flexibility of our approach on receptive fields, ion channels, and Hodgkin-Huxley models. We also characterize the space of circuit configurations giving rise to rhythmic activity in the crustacean stomatogastric ganglion, and use these results to derive hypotheses for underlying compensation mechanisms. Our approach will help close the gap between data-driven and theory-driven models of neural dynamics.

A Physiological Instability Displayed in Hippocampal Neurons Derived From Lithium-Nonresponsive Bipolar Disorder Patients

BACKGROUND

We recently reported a hyperexcitability phenotype displayed in dentate gyrus granule neurons derived from patients with bipolar disorder (BD) as well as a hyperexcitability that appeared only in CA3 pyramidal hippocampal neurons that were derived from patients with BD who responded to lithium treatment (lithium responders) and not in CA3 pyramidal hippocampal neurons that were derived from patients with BD who did not respond to lithium (nonresponders).

METHODS

Here we used our measurements of currents in neurons derived from 4 control subjects, 3 patients with BD who were lithium responders, and 3 patients with BD who were nonresponders. We changed the conductances of simulated dentate gyrus and CA3 hippocampal neurons according to our measurements to derive a numerical simulation for BD neurons.

RESULTS

The computationally simulated BD dentate gyrus neurons had a hyperexcitability phenotype similar to the experimental results. Only the simulated BD CA3 neurons derived from lithium responder patients were hyperexcitable. Interestingly, our computational model captured a physiological instability intrinsic to hippocampal neurons that were derived from nonresponder patients that we also observed when re-examining our experimental results. This instability was caused by a drastic reduction in the sodium current, accompanied by an increase in the amplitude of several potassium currents. These baseline alterations caused nonresponder BD hippocampal neurons to drastically shift their excitability with small changes to their sodium currents, alternating between hyperexcitable and hypoexcitable states.

CONCLUSIONS

Our computational model of BD hippocampal neurons that was based on our measurements reproduced the experimental phenotypes of hyperexcitability and physiological instability. We hypothesize that the physiological instability phenotype strongly contributes to affective lability in patients with BD.

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.

Ionic channel blockage in stochastic Hodgkin-Huxley neuronal model driven by multiple oscillatory signals

Ionic channel blockage and multiple oscillatory signals play an important role in the dynamical response of pulse sequences. The effects of ionic channel blockage and ionic channel noise on the discharge behaviors are studied in Hodgkin-Huxley neuronal model with multiple oscillatory signals. It is found that bifurcation points of spontaneous discharge are altered through tuning the amplitude of multiple oscillatory signals, and the discharge cycle is changed by increasing the frequency of multiple oscillatory signals. The effects of ionic channel blockage on neural discharge behaviors indicate that the neural excitability can be suppressed by the sodium channel blockage, however, the neural excitability can be reversed by the potassium channel blockage. There is an optimal blockage ratio of potassium channel at which the electrical activity is the most regular, while the order of neural spike is disrupted by the sodium channel blockage. In addition, the frequency of spike discharge is accelerated by increasing the ionic channel noise, the firing of neuron becomes more stable if the ionic channel noise is appropriately reduced. Our results might provide new insights into the effects of ionic channel blockages, multiple oscillatory signals, and ionic channel noises on neural discharge behaviors.

Rapid adaptation to elevated extracellular potassium in the pyloric circuit of the crab, Cancer borealis

Elevated potassium concentration ([K⁺]) is often used to alter excitability in neurons and networks by shifting the potassium equilibrium potential (E_K) and, consequently, the resting membrane potential. We studied the effects of increased extracellular [K⁺] on the well-described pyloric circuit of the crab Cancer borealis. A 2.5-fold increase in extracellular [K⁺] (2.5×[K⁺]) depolarized pyloric dilator (PD) neurons and resulted in short-term loss of their normal bursting activity. This period of silence was followed within 5-10 min by the recovery of spiking and/or bursting activity during continued superfusion of 2.5×[K⁺] saline. In contrast, when PD neurons were pharmacologically isolated from pyloric presynaptic inputs, they exhibited no transient loss of spiking activity in 2.5×[K⁺], suggesting the presence of an acute inhibitory effect mediated by circuit interactions. Action potential threshold in PD neurons hyperpolarized during an hour-long exposure to 2.5×[K⁺] concurrent with the recovery of spiking and/or bursting activity. Thus the initial loss of activity appears to be mediated by synaptic interactions within the network, but the secondary adaptation depends on changes in the intrinsic excitability of the pacemaker neurons. The complex sequence of events in the responses of pyloric neurons to elevated [K⁺] demonstrates that electrophysiological recordings are necessary to determine both the transient and longer term effects of even modest alterations of K⁺ concentrations on neuronal activity.NEW & NOTEWORTHY Solutions with elevated extracellular potassium are commonly used as a depolarizing stimulus. We studied the effects of high potassium concentration ([K⁺]) on the pyloric circuit of the crab stomatogastric ganglion. A 2.5-fold increase in extracellular [K⁺] caused a transient loss of activity that was not due to depolarization block, followed by a rapid increase in excitability and recovery of spiking within minutes. This suggests that changing extracellular potassium can have complex and nonstationary effects on neuronal circuits.

Degeneracy in the emergence of spike-triggered average of hippocampal pyramidal neurons

Hippocampal pyramidal neurons are endowed with signature excitability characteristics, exhibit theta-frequency selectivity - manifesting as impedance resonance and as a band-pass structure in the spike-triggered average (STA) - and coincidence detection tuned for gamma-frequency inputs. Are there specific constraints on molecular-scale (ion channel) properties in the concomitant emergence of cellular-scale encoding (feature detection and selectivity) and excitability characteristics? Here, we employed a biophysically-constrained unbiased stochastic search strategy involving thousands of conductance-based models, spanning 11 active ion channels, to assess the concomitant emergence of 14 different electrophysiological measurements. Despite the strong biophysical and physiological constraints, we found models that were similar in terms of their spectral selectivity, operating mode along the integrator-coincidence detection continuum and intrinsic excitability characteristics. The parametric combinations that resulted in these functionally similar models were non-unique with weak pair-wise correlations. Employing virtual knockout of individual ion channels in these functionally similar models, we found a many-to-many relationship between channels and physiological characteristics to mediate this degeneracy, and predicted a dominant role for HCN and transient potassium channels in regulating hippocampal neuronal STA. Our analyses reveals the expression of degeneracy, that results from synergistic interactions among disparate channel components, in the concomitant emergence of neuronal excitability and encoding characteristics.

Robust emergence of sharply tuned place-cell responses in hippocampal neurons with structural and biophysical heterogeneities

Hippocampal pyramidal neurons sustain propagation of fast electrical signals and are electrotonically non-compact structures exhibiting cell-to-cell variability in their complex dendritic arborization. In this study, we demonstrate that sharp place-field tuning and several somatodendritic functional maps concomitantly emerge despite the presence of geometrical heterogeneities in these neurons. We establish this employing an unbiased stochastic search strategy involving thousands of models that spanned several morphologies and distinct profiles of dispersed synaptic localization and channel expression. Mechanistically, employing virtual knockout models (VKMs), we explored the impact of bidirectional modulation in dendritic spike prevalence on place-field tuning sharpness. Consistent with the prior literature, we found that across all morphologies, virtual knockout of either dendritic fast sodium channels or N-methyl-D-aspartate receptors led to a reduction in dendritic spike prevalence, whereas A-type potassium channel knockouts resulted in a non-specific increase in dendritic spike prevalence. However, place-field tuning sharpness was critically impaired in all three sets of VKMs, demonstrating that sharpness in feature tuning is maintained by an intricate balance between mechanisms that promote and those that prevent dendritic spike initiation. From the functional standpoint of the emergence of sharp feature tuning and intrinsic functional maps, within this framework, geometric variability was compensated by a combination of synaptic democracy, the ability of randomly dispersed synapses to yield sharp tuning through dendritic spike initiation, and ion-channel degeneracy. Our results suggest electrotonically non-compact neurons to be endowed with several degrees of freedom, encompassing channel expression, synaptic localization and morphological microstructure, in achieving sharp feature encoding and excitability homeostasis.

Non-synaptic Plasticity in Leech Touch Cells

The role of Na⁺/K⁺-pumps in activity-dependent synaptic plasticity has been described in both vertebrates and invertebrates. Here, we provide evidence that the Na⁺/K⁺-pump is also involved in activity-dependent non-synaptic cellular plasticity in leech sensory neurons. We show that the resting membrane potential (RMP) of T cells hyperpolarizes in response to repeated somatic current injection, while at the same time their spike count (SC) and the input resistance (IR) increase. Our Hodgkin–Huxley-type neuron model, adjusted to physiological T cell properties, suggests that repetitive action potential discharges lead to increased Na⁺/K⁺-pump activity, which then hyperpolarizes the RMP. In consequence, a slow, non-inactivating current decreases, which is presumably mediated by voltage-dependent, low-threshold potassium channels. Closing of these putative M-type channels due to hyperpolarization of the resting potential increases the IR of the cell, leading to a larger number of spikes. By this mechanism, the response behavior switches from rapidly to slowly adapting spiking. These changes in spiking behavior also effect other T cells on the same side of the ganglion, which are connected via a combination of electrical and chemical synapses. An increased SC in the presynaptic T cell results in larger postsynaptic responses (PRs) in the other T cells. However, when the number of elicited presynaptic spikes is kept constant, the PR does not change. These results suggest that T cells change their responses in an activity-dependent manner through non-synaptic rather than synaptic plasticity. These changes might act as a gain-control mechanism. Depending on the previous activity, this gain could scale the relative impacts of synaptic inputs from other mechanoreceptors, versus the spike responses to tactile skin stimulation. This multi-tasking ability, and its flexible adaptation to previous activity, might make the T cell a key player in a preparatory network, enabling the leech to perform fast behavioral reactions to skin stimulation.

Degeneracy in hippocampal physiology and plasticity

Degeneracy, defined as the ability of structurally disparate elements to perform analogous function, has largely been assessed from the perspective of maintaining robustness of physiology or plasticity. How does the framework of degeneracy assimilate into an encoding system where the ability to change is an essential ingredient for storing new incoming information? Could degeneracy maintain the balance between the apparently contradictory goals of the need to change for encoding and the need to resist change towards maintaining homeostasis? In this review, we explore these fundamental questions with the mammalian hippocampus as an example encoding system. We systematically catalog lines of evidence, spanning multiple scales of analysis that point to the expression of degeneracy in hippocampal physiology and plasticity. We assess the potential of degeneracy as a framework to achieve the conjoint goals of encoding and homeostasis without cross-interferences. We postulate that biological complexity, involving interactions among the numerous parameters spanning different scales of analysis, could establish disparate routes towards accomplishing these conjoint goals. These disparate routes then provide several degrees of freedom to the encoding-homeostasis system in accomplishing its tasks in an input- and state-dependent manner. Finally, the expression of degeneracy spanning multiple scales offers an ideal reconciliation to several outstanding controversies, through the recognition that the seemingly contradictory disparate observations are merely alternate routes that the system might recruit towards accomplishment of its goals.