Homeostasis of functional maps in active dendrites emerges in the absence of individual channelostasis

Network motifs in cellular neurophysiology

Concepts from network science and graph theory, including the framework of network motifs, have been frequently applied in studying neuronal networks and other biological complex systems. Network-based approaches can also be used to study the functions of individual neurons, where cellular elements such as ion channels and membrane voltage are conceptualized as nodes within a network, and their interactions are denoted by edges. Network motifs in this context provide functional building blocks that help to illuminate the principles of cellular neurophysiology. In this review we build a case that network motifs operating within neurons provide tools for defining the functional architecture of single-neuron physiology and neuronal adaptations. We highlight the presence of such computational motifs in the cellular mechanisms underlying action potential generation, neuronal oscillations, dendritic integration, and neuronal plasticity. Future work applying the network motifs perspective may help to decipher the functional complexities of neurons and their adaptation during health and disease.

Dendritic effects of tDCS insights from a morphologically realistic model neuron

Transcranial direct current stimulation (tDCS) induces subcellular compartmental-dependent polarization, maximal in the distal portions of axons and dendrites. Using a morphologically realistic neuron model, we simulated tDCS-induced membrane polarization of the apical dendrite. Thus, we investigated the differential dendritic effects of anodal and cathodal tDCS on membrane potential polarization along the dendritic structure and its subsequent effects on dendritic membrane resistance, excitatory postsynaptic potential amplitude, backpropagating action potential amplitude, input/output relations, and long-term synaptic plasticity. We further showed that the effects of anodal and cathodal tDCS on the backpropagating action potential were asymmetric, and explained this asymmetry. Additionally, we showed that the effects on input/output relations were rather weak and limited to the low-mid range of stimulation frequencies, and that synaptic plasticity effects were mostly limited to the distal portion of the dendrite. Thus, we demonstrated how tDCS modifies dendritic physiology due to the dendrite's unique morphology and composition of voltage-gated ion channels.

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.

Simulations predict differing phase responses to excitation vs. inhibition in theta-resonant pyramidal neurons

Rhythmic activity is ubiquitous in neural systems, with theta-resonant pyramidal neurons integrating rhythmic inputs in many cortical structures. Impedance analysis has been widely used to examine frequency-dependent responses of neuronal membranes to rhythmic inputs, but it assumes that the neuronal membrane is a linear system, requiring the use of small signals to stay in a near-linear regime. However, postsynaptic potentials are often large and trigger nonlinear mechanisms (voltage-gated ion channels). The goals of this work were to 1) develop an analysis method to evaluate membrane responses in this nonlinear domain and 2) explore phase relationships between rhythmic stimuli and subthreshold and spiking membrane potential (Vmemb) responses in models of theta-resonant pyramidal neurons. Responses in these output regimes were asymmetrical, with different phase shifts during hyperpolarizing and depolarizing half-cycles. Suprathreshold theta-rhythmic stimuli produced nonstationary Vmemb responses. Sinusoidal inputs produced "phase retreat": action potentials occurred progressively later in cycles of the input stimulus, resulting from adaptation. Sinusoidal current with increasing amplitude over cycles produced "phase advance": action potentials occurred progressively earlier. Phase retreat, phase advance, and subthreshold phase shifts were modulated by multiple ion channel conductances. Our results suggest differential responses of cortical neurons depending on the frequency of oscillatory input, which will play a role in neuronal responses to shifts in network state. We hypothesize that intrinsic cellular properties complement network properties and contribute to in vivo phase-shift phenomena such as phase precession, seen in place and grid cells, and phase roll, also observed in hippocampal CA1 neurons.NEW & NOTEWORTHY We augmented electrical impedance analysis to characterize phase shifts between large-amplitude current stimuli and nonlinear, asymmetric membrane potential responses. We predict different frequency-dependent phase shifts in response excitation vs. inhibition, as well as shifts in spike timing over multiple input cycles, in theta-resonant pyramidal neurons. We hypothesize that these effects contribute to navigation-related phenomena such as phase precession and phase roll. Our neuron-level hypothesis complements, rather than falsifies, prior network-level explanations of these phenomena.

Degeneracy in epilepsy: multiple routes to hyperexcitable brain circuits and their repair

Due to its complex and multifaceted nature, developing effective treatments for epilepsy is still a major challenge. To deal with this complexity we introduce the concept of degeneracy to the field of epilepsy research: the ability of disparate elements to cause an analogous function or malfunction. Here, we review examples of epilepsy-related degeneracy at multiple levels of brain organisation, ranging from the cellular to the network and systems level. Based on these insights, we outline new multiscale and population modelling approaches to disentangle the complex web of interactions underlying epilepsy and to design personalised multitarget therapies.

Synergies between synaptic and HCN channel plasticity dictates firing rate homeostasis and mutual information transfer in hippocampal model neuron

Homeostasis is a precondition for any physiological system of any living organism. Nonetheless, models of learning and memory that are based on processes of synaptic plasticity are unstable by nature according to Hebbian rules, and it is not fully clear how homeostasis is maintained during these processes. This is where theoretical and computational frameworks can help in gaining a deeper understanding of the various cellular processes that enable homeostasis in the face of plasticity. A previous simplistic single compartmental model with a single synapse showed that maintaining input/output response homeostasis and stable synaptic learning could be enabled by introducing a linear relationship between synaptic plasticity and HCN conductance plasticity. In this study, we aimed to examine whether this approach could be extended to a more morphologically realistic model that entails multiple synapses and gradients of various VGICs. In doing so, we found that a linear relationship between synaptic plasticity and HCN conductance plasticity was able to maintain input/output response homeostasis in our morphologically realistic model, where the slope of the linear relationship was dependent on baseline HCN conductance and synaptic permeability values. An increase in either baseline HCN conductance or synaptic permeability value led to a decrease in the slope of the linear relationship. We further show that in striking contrast to the single compartment model, here linear relationship was insufficient in maintaining stable synaptic learning despite maintaining input/output response homeostasis. Additionally, we showed that homeostasis of input/output response profiles was at the expense of decreasing the mutual information transfer due to the increase in noise entropy, which could not be fully rescued by optimizing the linear relationship between synaptic and HCN conductance plasticity. Finally, we generated a place cell model based on theta oscillations and show that synaptic plasticity disrupts place cell activity. Whereas synaptic plasticity accompanied by HCN conductance plasticity through linear relationship maintains the stability of place cell activity. Our study establishes potential differences between a single compartmental model and a morphologically realistic model.

Heterogeneous stochastic bifurcations explain intrinsic oscillatory patterns in entorhinal cortical stellate cells

Stellate cells (SC) in the medial entorhinal cortex manifest intrinsic membrane potential oscillatory patterns. Although different theoretical frameworks have been proposed to explain these patterns, a robust unifying framework that jointly accounts for intrinsic heterogeneities and stochasticity is missing. Here, we first performed in vitro patch-clamp electrophysiological recordings from rat SCs and found pronounced cell-to-cell variability in their characteristic physiological properties, including peri-threshold oscillatory patterns. We demonstrate that noise introduced into two independent populations (endowed with deterministic or stochastic ion-channel gating kinetics) of heterogeneous biophysical models yielded activity patterns that were qualitatively similar to electrophysiological peri-threshold oscillatory activity in SCs. We developed spectrogram-based quantitative metrics for the identification of valid oscillations and confirmed that these metrics reliably captured the variable-amplitude and arhythmic oscillatory patterns observed in electrophysiological recordings. Using these quantitative metrics, we validated activity patterns from both heterogeneous populations of SC models, with each model assessed with multiple trials of different levels of noise at distinct membrane depolarizations. Our analyses unveiled the manifestation of stochastic resonance (detection of the highest number of valid oscillatory traces at an optimal level of noise) in both heterogeneous populations of SC models. Finally, we show that a generalized network motif comprised of a slow negative feedback loop amplified by a fast positive feedback loop manifested stochastic bifurcations and stochastic resonance in the emergence of oscillations. Together, through a unique convergence of the degeneracy and stochastic resonance frameworks, our unifying framework centered on heterogeneous stochastic bifurcations argues for state-dependent emergence of SC oscillations.

Efficient information coding and degeneracy in the nervous system

Efficient information coding (EIC) is a universal biological framework rooted in the fundamental principle that system responses should match their natural stimulus statistics for maximizing environmental information. Quantitatively assessed through information theory, such adaptation to the environment occurs at all biological levels and timescales. The context dependence of environmental stimuli and the need for stable adaptations make EIC a daunting task. We argue that biological complexity is the principal architect that subserves deft execution of stable EIC. Complexity in a system is characterized by several functionally segregated subsystems that show a high degree of functional integration when they interact with each other. Complex biological systems manifest heterogeneities and degeneracy, wherein structurally different subsystems could interact to yield the same functional outcome. We argue that complex systems offer several choices that effectively implement EIC and homeostasis for each of the different contexts encountered by the system.

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.

The voltage and spiking responses of subthreshold resonant neurons to structured and fluctuating inputs: persistence and loss of resonance and variability

We systematically investigate the response of neurons to oscillatory currents and synaptic-like inputs and we extend our investigation to non-structured synaptic-like spiking inputs with more realistic distributions of presynaptic spike times. We use two types of chirp-like inputs consisting of (i) a sequence of cycles with discretely increasing frequencies over time, and (ii) a sequence having the same cycles arranged in an arbitrary order. We develop and use a number of frequency-dependent voltage response metrics to capture the different aspects of the voltage response, including the standard impedance (Z) and the peak-to-trough amplitude envelope ([Formula: see text]) profiles. We show that Z-resonant cells (cells that exhibit subthreshold resonance in response to sinusoidal inputs) also show [Formula: see text]-resonance in response to sinusoidal inputs, but generally do not (or do it very mildly) in response to square-wave and synaptic-like inputs. In the latter cases the resonant response using Z is not predictive of the preferred frequencies at which the neurons spike when the input amplitude is increased above subthreshold levels. We also show that responses to conductance-based synaptic-like inputs are attenuated as compared to the response to current-based synaptic-like inputs, thus providing an explanation to previous experimental results. These response patterns were strongly dependent on the intrinsic properties of the participating neurons, in particular whether the unperturbed Z-resonant cells had a stable node or a focus. In addition, we show that variability emerges in response to chirp-like inputs with arbitrarily ordered patterns where all signals (trials) in a given protocol have the same frequency content and the only source of uncertainty is the subset of all possible permutations of cycles chosen for a given protocol. This variability is the result of the multiple different ways in which the autonomous transient dynamics is activated across cycles in each signal (different cycle orderings) and across trials. We extend our results to include high-rate Poisson distributed current- and conductance-based synaptic inputs and compare them with similar results using additive Gaussian white noise. We show that the responses to both Poisson-distributed synaptic inputs are attenuated with respect to the responses to Gaussian white noise. For cells that exhibit oscillatory responses to Gaussian white noise (band-pass filters), the response to conductance-based synaptic inputs are low-pass filters, while the response to current-based synaptic inputs may remain band-pass filters, consistent with experimental findings. Our results shed light on the mechanisms of communication of oscillatory activity among neurons in a network via subthreshold oscillations and resonance and the generation of network resonance.

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.

Integrative Models of Brain Structure and Dynamics: Concepts, Challenges, and Methods

The anatomical architecture of the brain constrains the dynamics of interactions between various regions. On a microscopic scale, neural plasticity regulates the connections between individual neurons. This microstructural adaptation facilitates coordinated dynamics of populations of neurons (mesoscopic scale) and brain regions (macroscopic scale). However, the mechanisms acting on multiple timescales that govern the reciprocal relationship between neural network structure and its intrinsic dynamics are not well understood. Studies empirically investigating such relationships on the whole-brain level rely on macroscopic measurements of structural and functional connectivity estimated from various neuroimaging modalities such as Diffusion-weighted Magnetic Resonance Imaging (dMRI), Electroencephalography (EEG), Magnetoencephalography (MEG), and functional Magnetic Resonance Imaging (fMRI). dMRI measures the anisotropy of water diffusion along axonal fibers, from which structural connections are estimated. EEG and MEG signals measure electrical activity and magnetic fields induced by the electrical activity, respectively, from various brain regions with a high temporal resolution (but limited spatial coverage), whereas fMRI measures regional activations indirectly via blood oxygen level-dependent (BOLD) signals with a high spatial resolution (but limited temporal resolution). There are several studies in the neuroimaging literature reporting statistical associations between macroscopic structural and functional connectivity. On the other hand, models of large-scale oscillatory dynamics conditioned on network structure (such as the one estimated from dMRI connectivity) provide a platform to probe into the structure-dynamics relationship at the mesoscopic level. Such investigations promise to uncover the theoretical underpinnings of the interplay between network structure and dynamics and could be complementary to the macroscopic level inquiries. In this article, we review theoretical and empirical studies that attempt to elucidate the coupling between brain structure and dynamics. Special attention is given to various clinically relevant dimensions of brain connectivity such as the topological features and neural synchronization, and their applicability for a given modality, spatial or temporal scale of analysis is discussed. Our review provides a summary of the progress made along this line of research and identifies challenges and promising future directions for multi-modal neuroimaging analyses.

Active Dendrites and Local Field Potentials: Biophysical Mechanisms and Computational Explorations

Neurons and glial cells are endowed with membranes that express a rich repertoire of ion channels, transporters, and receptors. The constant flux of ions across the neuronal and glial membranes results in voltage fluctuations that can be recorded from the extracellular matrix. The high frequency components of this voltage signal contain information about the spiking activity, reflecting the output from the neurons surrounding the recording location. The low frequency components of the signal, referred to as the local field potential (LFP), have been traditionally thought to provide information about the synaptic inputs that impinge on the large dendritic trees of various neurons. In this review, we discuss recent computational and experimental studies pointing to a critical role of several active dendritic mechanisms that can influence the genesis and the location-dependent spectro-temporal dynamics of LFPs, spanning different brain regions. We strongly emphasize the need to account for the several fast and slow dendritic events and associated active mechanisms - including gradients in their expression profiles, inter- and intra-cellular spatio-temporal interactions spanning neurons and glia, heterogeneities and degeneracy across scales, neuromodulatory influences, and activitydependent plasticity - towards gaining important insights about the origins of LFP under different behavioral states in health and disease. We provide simple but essential guidelines on how to model LFPs taking into account these dendritic mechanisms, with detailed methodology on how to account for various heterogeneities and electrophysiological properties of neurons and synapses while studying LFPs.

Neural circuit function redundancy in brain disorders

Redundancy is a ubiquitous property of the nervous system. This means that vastly different configurations of cellular and synaptic components can enable the same neural circuit functions. However, until recently, very little brain disorder research has considered the implications of this characteristic when designing experiments or interpreting data. Here, we first summarise the evidence for redundancy in healthy brains, explaining redundancy and three related sub-concepts: sloppiness, dependencies and multiple solutions. We then lay out key implications for brain disorder research, covering recent examples of redundancy effects in experimental studies on psychiatric disorders. Finally, we give predictions for future experiments based on these concepts.

Ion-channel degeneracy: Multiple ion channels heterogeneously regulate intrinsic physiology of rat hippocampal granule cells

Degeneracy, the ability of multiple structural components to elicit the same characteristic functional properties, constitutes an elegant mechanism for achieving biological robustness. In this study, we sought electrophysiological signatures for the expression of ion-channel degeneracy in the emergence of intrinsic properties of rat hippocampal granule cells. We measured the impact of four different ion-channel subtypes-hyperpolarization-activated cyclic-nucleotide-gated (HCN), barium-sensitive inward rectifier potassium (Kir ), tertiapin-Q-sensitive inward rectifier potassium, and persistent sodium (NaP) channels-on 21 functional measurements employing pharmacological agents, and report electrophysiological data on two characteristic signatures for the expression of ion-channel degeneracy in granule cells. First, the blockade of a specific ion-channel subtype altered several, but not all, functional measurements. Furthermore, any given functional measurement was altered by the blockade of many, but not all, ion-channel subtypes. Second, the impact of blocking each ion-channel subtype manifested neuron-to-neuron variability in the quantum of changes in the electrophysiological measurements. Specifically, we found that blocking HCN or Ba-sensitive Kir channels enhanced action potential firing rate, but blockade of NaP channels reduced firing rate of granule cells. Subthreshold measures of granule cell intrinsic excitability (input resistance, temporal summation, and impedance amplitude) were enhanced by blockade of HCN or Ba-sensitive Kir channels, but were not significantly altered by NaP channel blockade. We confirmed that the HCN and Ba-sensitive Kir channels independently altered sub- and suprathreshold properties of granule cells through sequential application of pharmacological agents that blocked these channels. Finally, we found that none of the sub- or suprathreshold measurements of granule cells were significantly altered upon treatment with tertiapin-Q. Together, the heterogeneous many-to-many mapping between ion channels and single-neuron intrinsic properties emphasizes the need to account for ion-channel degeneracy in cellular- and network-scale physiology.

Spatial information transfer in hippocampal place cells depends on trial-to-trial variability, symmetry of place-field firing, and biophysical heterogeneities

The relationship between the feature-tuning curve and information transfer profile of individual neurons provides vital insights about neural encoding. However, the relationship between the spatial tuning curve and spatial information transfer of hippocampal place cells remains unexplored. Here, employing a stochastic search procedure spanning thousands of models, we arrived at 127 conductance-based place-cell models that exhibited signature electrophysiological characteristics and sharp spatial tuning, with parametric values that exhibited neither clustering nor strong pairwise correlations. We introduced trial-to-trial variability in responses and computed model tuning curves and information transfer profiles, using stimulus-specific (SSI) and mutual (MI) information metrics, across locations within the place field. We found spatial information transfer to be heterogeneous across models, but to reduce consistently with increasing levels of variability. Importantly, whereas reliable low-variability responses implied that maximal information transfer occurred at high-slope regions of the tuning curve, increase in variability resulted in maximal transfer occurring at the peak-firing location in a subset of models. Moreover, experience-dependent asymmetry in place-field firing introduced asymmetries in the information transfer computed through MI, but not SSI, and the impact of activity-dependent variability on information transfer was minimal compared to activity-independent variability. We unveiled ion-channel degeneracy in the regulation of spatial information transfer, and demonstrated critical roles for N-methyl-d-aspartate receptors, transient potassium and dendritic sodium channels in regulating information transfer. Our results demonstrate that trial-to-trial variability, tuning-curve shape and biological heterogeneities critically regulate the relationship between the spatial tuning curve and spatial information transfer in hippocampal place 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.

Efficient phase coding in hippocampal place cells

Neural codes have been postulated to build efficient representations of the external world. The hippocampus, an encoding system, employs neuronal firing rates and spike phases to encode external space. Although the biophysical origin of such codes is at a single neuronal level, the role of neural components in efficient coding is not understood. The complexity of this problem lies in the dimensionality of the parametric space encompassing neural components, and is amplified by the enormous biological heterogeneity observed in each parameter. A central question that spans encoding systems therefore is how neurons arrive at efficient codes in the face of widespread biological heterogeneities. To answer this, we developed a conductance-based spiking model for phase precession, a phase code of external space exhibited by hippocampal place cells. Our model accounted for several experimental observations on place cell firing and electrophysiology: the emergence of phase precession from exact spike timings of conductance-based models with neuron-specific ion channels and receptors; biological heterogeneities in neural components and excitability; the emergence of subthreshold voltage ramp, increased firing rate, enhanced theta power within the place field; a signature reduction in extracellular theta frequency compared to its intracellular counterpart; and experience-dependent asymmetry in firing-rate profile. We formulated phase-coding efficiency, using Shannon's information theory, as an information maximization problem with spike phase as the response and external space within a single place field as the stimulus. We employed an unbiased stochastic search spanning an 11-dimensional neural space, involving thousands of iterations that accounted for the biophysical richness and neuron-to-neuron heterogeneities. We found a small subset of models that exhibited efficient spatial information transfer through the phase code, and investigated the distinguishing features of this subpopulation at the parametric and functional scales. At the parametric scale, which spans the molecular components that defined the neuron, several nonunique parametric combinations with weak pairwise correlations yielded models with similar high phase-coding efficiency. Importantly, placing additional constraints on these models in terms of matching other aspects of hippocampal neural responses did not hamper parametric degeneracy. We provide quantitative evidence demonstrating this parametric degeneracy to be a consequence of a many-to-one relationship between the different parameters and phase-coding efficiency. At the functional scale, involving the cellular-scale neural properties, our analyses revealed an important higher-order constraint that was exclusive to models exhibiting efficient phase coding. Specifically, we found a counterbalancing negative correlation between neuronal gain and the strength of external synaptic inputs as a critical functional constraint for the emergence of efficient phase coding. These observations implicate intrinsic neural properties as important contributors in effectuating such counterbalance, which can be achieved by recruiting nonunique parametric combinations. Finally, we show that a change in afferent statistics, manifesting as input asymmetry onto these neuronal models, induced an adaptive shift in the phase code that preserved its efficiency. Together, our analyses unveil parametric degeneracy as a mechanism to harness widespread neuron-to-neuron heterogeneity towards accomplishing stable and efficient encoding, provided specific higher-order functional constraints on the relationship of neural gain to external inputs are satisfied.

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.

Heterogeneities in intrinsic excitability and frequency-dependent response properties of granule cells across the blades of the rat dentate gyrus

The dentate gyrus (DG), the input gate to the hippocampus proper, is anatomically segregated into three different sectors, namely, the suprapyramidal blade, the crest region, and the infrapyramidal blade. Although there are well-established differences between these sectors in terms of neuronal morphology, connectivity patterns, and activity levels, differences in electrophysiological properties of granule cells within these sectors have remained unexplored. Here, employing somatic whole cell patch-clamp recordings from the rat DG, we demonstrate that granule cells in these sectors manifest considerable heterogeneities in their intrinsic excitability, temporal summation, action potential characteristics, and frequency-dependent response properties. Across sectors, these neurons showed positive temporal summation of their responses to inputs mimicking excitatory postsynaptic currents and showed little to no sag in their voltage responses to pulse currents. Consistently, the impedance amplitude profile manifested low-pass characteristics and the impedance phase profile lacked positive phase values at all measured frequencies and voltages and for all sectors. Granule cells in all sectors exhibited class I excitability, with broadly linear firing rate profiles, and granule cells in the crest region fired significantly fewer action potentials compared with those in the infrapyramidal blade. Finally, we found weak pairwise correlations across the 18 different measurements obtained individually from each of the three sectors, providing evidence that these measurements are indeed reporting distinct aspects of neuronal physiology. Together, our analyses show that granule cells act as integrators of afferent information and emphasize the need to account for the considerable physiological heterogeneities in assessing their roles in information encoding and processing.NEW & NOTEWORTHY We employed whole cell patch-clamp recordings from granule cells in the three subregions of the rat dentate gyrus to demonstrate considerable heterogeneities in their intrinsic excitability, temporal summation, action potential characteristics, and frequency-dependent response properties. Across sectors, granule cells did not express membrane potential resonance, and their impedance profiles lacked inductive phase leads at all measured frequencies. Our analyses also show that granule cells manifest class I excitability characteristics, categorizing them as integrators of afferent information.

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.

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.

Simple models of quantitative firing phenotypes in hippocampal neurons: Comprehensive coverage of intrinsic diversity

Patterns of periodic voltage spikes elicited by a neuron help define its dynamical identity. Experimentally recorded spike trains from various neurons show qualitatively distinguishable features such as delayed spiking, spiking with or without frequency adaptation, and intrinsic bursting. Moreover, the input-dependent responses of a neuron not only show different quantitative features, such as higher spike frequency for a stronger input current injection, but can also exhibit qualitatively different responses, such as spiking and bursting under different input conditions, thus forming a complex phenotype of responses. In previous work, the comprehensive knowledge base of hippocampal neuron types Hippocampome.org systematically characterized various spike pattern phenotypes experimentally identified from 120 neuron types/subtypes. In this paper, we present a complete set of simple phenomenological models that quantitatively reproduce the diverse and complex phenotypes of hippocampal neurons. In addition to point-neuron models, we created compact multi-compartment models with up to four compartments, which will allow spatial segregation of synaptic integration in network simulations. Electrotonic compartmentalization observed in our compact multi-compartment models is qualitatively consistent with experimental observations. The models were created using an automated pipeline based on evolutionary algorithms. This work maps 120 neuron types/subtypes in the rodent hippocampus to a low-dimensional model space and adds another dimension to the knowledge accumulated in Hippocampome.org. Computationally efficient representations of intrinsic dynamics, along with other pieces of knowledge available in Hippocampome.org, provide a biologically realistic platform to explore the large-scale interactions of various neuron types at the mesoscopic level.

Frequency-dependent responses of neuronal models to oscillatory inputs in current versus voltage clamp

Action potential generation in neurons depends on a membrane potential threshold and therefore on how subthreshold inputs influence this voltage. In oscillatory networks, for example, many neuron types have been shown to produce membrane potential ([Formula: see text]) resonance: a maximum subthreshold response to oscillatory inputs at a nonzero frequency. Resonance is usually measured by recording [Formula: see text] in response to a sinusoidal current ([Formula: see text]), applied at different frequencies (f), an experimental setting known as current clamp (I-clamp). Several recent studies, however, use the voltage clamp (V-clamp) method to control [Formula: see text] with a sinusoidal input at different frequencies [[Formula: see text]] and measure the total membrane current ([Formula: see text]). The two methods obey systems of differential equations of different dimensionality, and while I-clamp provides a measure of electrical impedance [[Formula: see text]], V-clamp measures admittance [[Formula: see text]]. We analyze the relationship between these two measurement techniques. We show that, despite different dimensionality, in linear systems the two measures are equivalent: [Formula: see text]. However, nonlinear model neurons produce different values for Z and [Formula: see text]. In particular, nonlinearities in the voltage equation produce a much larger difference between these two quantities than those in equations of recovery variables that describe activation and inactivation kinetics. Neurons are inherently nonlinear, and notably, with ionic currents that amplify resonance, the voltage clamp technique severely underestimates the current clamp response. We demonstrate this difference experimentally using the PD neurons in the crab stomatogastric ganglion. These findings are instructive for researchers who explore cellular mechanisms of neuronal oscillations.

Disparate forms of heterogeneities and interactions among them drive channel decorrelation in the dentate gyrus: Degeneracy and dominance

The ability of a neuronal population to effectuate channel decorrelation, which is one form of response decorrelation, has been identified as an essential prelude to efficient neural encoding. To what extent are diverse forms of local and afferent heterogeneities essential in accomplishing channel decorrelation in the dentate gyrus (DG)? Here, we incrementally incorporated four distinct forms of biological heterogeneities into conductance-based network models of the DG and systematically delineate their relative contributions to channel decorrelation. First, to effectively incorporate intrinsic heterogeneities, we built physiologically validated heterogeneous populations of granule (GC) and basket cells (BC) through independent stochastic search algorithms spanning exhaustive parametric spaces. These stochastic search algorithms, which were independently constrained by experimentally determined ion channels and by neurophysiological signatures, revealed cellular-scale degeneracy in the DG. Specifically, in GC and BC populations, disparate parametric combinations yielded similar physiological signatures, with underlying parameters exhibiting significant variability and weak pair-wise correlations. Second, we introduced synaptic heterogeneities through randomization of local synaptic strengths. Third, in including adult neurogenesis, we subjected the valid model populations to randomized structural plasticity and matched neuronal excitability to electrophysiological data. We assessed networks comprising different combinations of these three local heterogeneities with identical or heterogeneous afferent inputs from the entorhinal cortex. We found that the three forms of local heterogeneities were independently and synergistically capable of mediating significant channel decorrelation when the network was driven by identical afferent inputs. However, when we incorporated afferent heterogeneities into the network to account for the divergence in DG afferent connectivity, the impact of all three forms of local heterogeneities was significantly suppressed by the dominant role of afferent heterogeneities in mediating channel decorrelation. Our results unveil a unique convergence of cellular- and network-scale degeneracy in the emergence of channel decorrelation in the DG, whereby disparate forms of local and afferent heterogeneities could synergistically drive input discriminability.

Membrane potential resonance in non-oscillatory neurons interacts with synaptic connectivity to produce network oscillations

Several neuron types have been shown to exhibit (subthreshold) membrane potential resonance (MPR), defined as the occurrence of a peak in their voltage amplitude response to oscillatory input currents at a preferred (resonant) frequency. MPR has been investigated both experimentally and theoretically. However, whether MPR is simply an epiphenomenon or it plays a functional role for the generation of neuronal network oscillations and how the latent time scales present in individual, non-oscillatory cells affect the properties of the oscillatory networks in which they are embedded are open questions. We address these issues by investigating a minimal network model consisting of (i) a non-oscillatory linear resonator (band-pass filter) with 2D dynamics, (ii) a passive cell (low-pass filter) with 1D linear dynamics, and (iii) nonlinear graded synaptic connections (excitatory or inhibitory) with instantaneous dynamics. We demonstrate that (i) the network oscillations crucially depend on the presence of MPR in the resonator, (ii) they are amplified by the network connectivity, (iii) they develop relaxation oscillations for high enough levels of mutual inhibition/excitation, and (iv) the network frequency monotonically depends on the resonators resonant frequency. We explain these phenomena using a reduced adapted version of the classical phase-plane analysis that helps uncovering the type of effective network nonlinearities that contribute to the generation of network oscillations. We extend our results to networks having cells with 2D dynamics. Our results have direct implications for network models of firing rate type and other biological oscillatory networks (e.g, biochemical, genetic).

Active dendrites regulate the spatiotemporal spread of signaling microdomains

Microdomains that emerge from spatially constricted spread of biochemical signaling components play a central role in several neuronal computations. Although dendrites, endowed with several voltage-gated ion channels, form a prominent structural substrate for microdomain physiology, it is not known if these channels regulate the spatiotemporal spread of signaling microdomains. Here, we employed a multiscale, morphologically realistic, conductance-based model of the hippocampal pyramidal neuron that accounted for experimental details of electrical and calcium-dependent biochemical signaling. We activated synaptic N-Methyl-d-Aspartate receptors through theta-burst stimulation (TBS) or pairing (TBP) and assessed microdomain propagation along a signaling pathway that included calmodulin, calcium/calmodulin-dependent protein kinase II (CaMKII) and protein phosphatase 1. We found that the spatiotemporal spread of the TBS-evoked microdomain in phosphorylated CaMKII (pCaMKII) was amplified in comparison to that of the corresponding calcium microdomain. Next, we assessed the role of two dendritically expressed inactivating channels, one restorative (A-type potassium) and another regenerative (T-type calcium), by systematically varying their conductances. Whereas A-type potassium channels suppressed the spread of pCaMKII microdomains by altering the voltage response to TBS, T-type calcium channels enhanced this spread by modulating TBS-induced calcium influx without changing the voltage. Finally, we explored cross-dependencies of these channels with other model components, and demonstrated the heavy mutual interdependence of several biophysical and biochemical properties in regulating microdomains and their spread. Our conclusions unveil a pivotal role for dendritic voltage-gated ion channels in actively amplifying or suppressing biochemical signals and their spatiotemporal spread, with critical implications for clustered synaptic plasticity, robust information transfer and efficient neural coding.

The spatiotemporal spread of biochemical signals in neurons and other cells regulate signaling specificity, tuning of signal propagation, along with specificity and clustering of adaptive plasticity. Theoretical and experimental studies have demonstrated a critical role for cellular morphology and the topology of signaling networks in regulating this spread. In this study, we add a significantly complex dimension to this narrative by demonstrating that voltage-gated ion channels on the plasma membrane could actively amplify or suppress the strength and spread of downstream signaling components. Given the expression of different ion channels with wide-ranging heterogeneity in gating kinetics, localization and density, our results point to an increase in complexity of and degeneracy in signaling spread, and unveil a powerful mechanism for regulating biochemical-signaling pathways across different cell types.

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.

Spatially dispersed synapses yield sharply-tuned place cell responses through dendritic spike initiation

KEY POINTS

The generation of dendritic spikes and the consequent sharp tuning of neuronal responses are together attainable even when iso-feature synapses are randomly dispersed across the dendritic arbor. Disparate combinations of channel conductances with distinct configurations of randomly dispersed place field synapses concomitantly yield similar sharp tuning profiles and similar functional maps of several intrinsic properties. Targeted synaptic plasticity converts silent cells to place cells for specific place fields in models with disparate channel combinations that receive dispersed synaptic inputs from multiple place field locations. Dispersed localization of iso-feature synapses is a strong candidate for achieving sharp feature selectivity in neurons across sensory-perceptual systems, with several degrees of freedom in relation to synaptic locations. Quantitative evidence for the possibility that degeneracy (i.e. the ability of disparate structural components to yield similar functional outcomes) could act as a broad framework that effectively accomplishes the twin goals of input-feature encoding and homeostasis of intrinsic properties without cross interferences.

ABSTRACT

A prominent hypothesis spanning several sensory-perceptual systems implicates spatially clustered synapses in the generation of dendritic spikes that mediate sharply-tuned neuronal responses to input features. In this conductance-based morphologically-precise computational study, we tested this hypothesis by systematically analysing the impact of distinct synaptic and channel localization profiles on sharpness of spatial tuning in hippocampal pyramidal neurons. We found that the generation of dendritic spikes, the emergence of an excitatory ramp in somatic voltage responses, the expression of several intrinsic somatodendritic functional maps and sharp tuning of place-cell responses were all attainable even when iso-feature synapses are randomly dispersed across the dendritic arbor of models with disparate channel combinations. Strikingly, the generation and propagation of dendritic spikes, reliant on dendritic sodium channels and N-methyl-d-asparate receptors, mediated the sharpness of spatial tuning achieved with dispersed synaptic localization. To ensure that our results were not artefacts of narrow parametric choices, we confirmed these conclusions with independent multiparametric stochastic search algorithms spanning thousands of unique models for each synaptic localization scenario. Next, employing virtual knockout models, we demonstrated a vital role for dendritically expressed voltage-gated ion channels, especially the transient potassium channels, in maintaining sharpness of place-cell tuning. Importantly, we established that synaptic potentiation targeted to afferents from one specific place field was sufficient to impart place field selectivity even when intrinsically disparate neurons received randomly dispersed afferents from multiple place field locations. Our results provide quantitative evidence for disparate combinations of channel and synaptic localization profiles to concomitantly yield similar tuning and similar intrinsic properties.

Degeneracy in the robust expression of spectral selectivity, subthreshold oscillations, and intrinsic excitability of entorhinal stellate cells

Biological heterogeneities are ubiquitous and play critical roles in the emergence of physiology at multiple scales. Although neurons in layer II (LII) of the medial entorhinal cortex (MEC) express heterogeneities in channel properties, the impact of such heterogeneities on the robustness of their cellular-scale physiology has not been assessed. Here, we performed a 55-parameter stochastic search spanning nine voltage- or calcium-activated channels to assess the impact of channel heterogeneities on the concomitant emergence of 10 in vitro electrophysiological characteristics of LII stellate cells (SCs). We generated 150,000 models and found a heterogeneous subpopulation of 449 valid models to robustly match all electrophysiological signatures. We employed this heterogeneous population to demonstrate the emergence of cellular-scale degeneracy in SCs, whereby disparate parametric combinations expressing weak pairwise correlations resulted in similar models. We then assessed the impact of virtually knocking out each channel from all valid models and demonstrate that the mapping between channels and measurements was many-to-many, a critical requirement for the expression of degeneracy. Finally, we quantitatively predict that the spike-triggered average of SCs should be endowed with theta-frequency spectral selectivity and coincidence detection capabilities in the fast gamma-band. We postulate this fast gamma-band coincidence detection as an instance of cellular-scale-efficient coding, whereby SC response characteristics match the dominant oscillatory signals in LII MEC. The heterogeneous population of valid SC models built here unveils the robust emergence of cellular-scale physiology despite significant channel heterogeneities, and forms an efficacious substrate for evaluating the impact of biological heterogeneities on entorhinal network function.

NEW & NOTEWORTHY We assessed the impact of heterogeneities in channel properties on the robustness of cellular-scale physiology of medial entorhinal cortical stellate neurons. We demonstrate that neuronal models with disparate channel combinations were endowed with similar physiological characteristics, as a consequence of the many-to-many mapping between channel properties and the physiological characteristics that they modulate. We predict that the spike-triggered average of stellate cells should be endowed with theta-frequency spectral selectivity and fast gamma-band coincidence detection capabilities.

T2N as a new tool for robust electrophysiological modeling demonstrated for mature and adult-born dentate granule cells

Compartmental models are the theoretical tool of choice for understanding single neuron computations. However, many models are incomplete, built ad hoc and require tuning for each novel condition rendering them of limited usability. Here, we present T2N, a powerful interface to control NEURON with Matlab and TREES toolbox, which supports generating models stable over a broad range of reconstructed and synthetic morphologies. We illustrate this for a novel, highly detailed active model of dentate granule cells (GCs) replicating a wide palette of experiments from various labs. By implementing known differences in ion channel composition and morphology, our model reproduces data from mouse or rat, mature or adult-born GCs as well as pharmacological interventions and epileptic conditions. This work sets a new benchmark for detailed compartmental modeling. T2N is suitable for creating robust models useful for large-scale networks that could lead to novel predictions. We discuss possible T2N application in degeneracy studies.

Spiking resonances in models with the same slow resonant and fast amplifying currents but different subthreshold dynamic properties

The generation of spiking resonances in neurons (preferred spiking responses to oscillatory inputs) requires the interplay of the intrinsic ionic currents that operate at the subthreshold voltage level and the spiking mechanisms. Combinations of the same types of ionic currents in different parameter regimes may give rise to different types of nonlinearities in the voltage equation (e.g., parabolic- and cubic-like), generating subthreshold (membrane potential) oscillations patterns with different properties. These nonlinearities are not apparent in the model equations, but can be uncovered by plotting the voltage nullclines in the phase-plane diagram. We investigate the spiking resonant properties of conductance-based models that are biophysically equivalent at the subthreshold level (same ionic currents), but dynamically different (parabolic- and cubic-like voltage nullclines). As a case study we consider a model having a persistent sodium and a hyperpolarization-activated (h-) currents, which exhibits subthreshold resonance in the theta frequency band. We unfold the concept of spiking resonance into evoked and output spiking resonance. The former focuses on the input frequencies that are able to generate spikes, while the latter focuses on the output spiking frequencies regardless of the input frequency that generated these spikes. A cell can exhibit one or both types of resonances. We also measure spiking phasonance, which is an extension of subthreshold phasonance (zero-phase-shift response to oscillatory inputs) to the spiking regime. The subthreshold resonant properties of both types of models are communicated to the spiking regime for low enough input amplitudes as the voltage response for the subthreshold resonant frequency band raises above threshold. For higher input amplitudes evoked spiking resonance is no longer present in these models, but output spiking resonance is present primarily in the parabolic-like model due to a cycle skipping mechanism (involving mixed-mode oscillations), while the cubic-like model shows a better 1:1 entrainment. We use dynamical systems tools to explain the underlying mechanisms and the mechanistic differences between the resonance types. Our results demonstrate that the effective time scales that operate at the subthreshold regime to generate intrinsic subthreshold oscillations, mixed-mode oscillations and subthreshold resonance do not necessarily determine the existence of a preferred spiking response to oscillatory inputs in the same frequency band. The results discussed in this paper highlight both the complexity of the suprathreshold responses to oscillatory inputs in neurons having resonant and amplifying currents with different time scales and the fact that the identity of the participating ionic currents is not enough to predict the resulting patterns, but additional dynamic information, captured by the geometric properties of the phase-space diagram, is needed.

Theta-frequency selectivity in the somatic spike-triggered average of rat hippocampal pyramidal neurons is dependent on HCN channels

The ability to distill specific frequencies from complex spatiotemporal patterns of afferent inputs is a pivotal functional requirement for neurons residing in networks receiving frequency-multiplexed inputs. Although the expression of theta-frequency subthreshold resonance is established in hippocampal pyramidal neurons, it is not known if their spike initiation dynamics manifest spectral selectivity, or if their intrinsic properties are tuned to process gamma-frequency inputs. Here, we measured the spike-triggered average (STA) of rat hippocampal pyramidal neurons through electrophysiological recordings and quantified spectral selectivity in their spike initiation dynamics and their coincidence detection window (CDW). Our results revealed strong theta-frequency selectivity in the STA, which was also endowed with gamma-range CDW, with prominent neuron-to-neuron variability that manifested distinct pairwise dissociations and correlations with different intrinsic measurements. Furthermore, we demonstrate that the STA and its measurements substantially adapted to the state of the neuron defined by its membrane potential and to the statistics of its afferent inputs. Finally, we tested the effect of pharmacologically blocking the hyperpolarization-activated cyclic-nucleotide-gated (HCN) channels on the STA and found that the STA characteristic frequency reduced significantly to the delta-frequency band after HCN channel blockade. This delta-frequency selectivity in the STA emerged in the absence of subthreshold resonance, which was abolished by HCN channel blockade, thereby confirming computational predictions on the dissociation between these two forms of spectral selectivity. Our results expand the roles of HCN channels to theta-frequency selectivity in the spike initiation dynamics, apart from underscoring the critical role of interactions among different ion channels in regulating neuronal physiology.NEW & NOTEWORTHY We had previously predicted, using computational analyses, that the spike-triggered average (STA) of hippocampal neurons would exhibit theta-frequency (4-10 Hz) spectral selectivity and would manifest coincidence detection capabilities for inputs in the gamma-frequency band (25-150 Hz). Here, we confirmed these predictions through direct electrophysiological recordings of STA from rat CA1 pyramidal neurons and demonstrate that blocking HCN channels reduces the frequency of STA spectral selectivity to the delta-frequency range (0.5-4 Hz).

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.

Resonance modulation, annihilation and generation of anti-resonance and anti-phasonance in 3D neuronal systems: interplay of resonant and amplifying currents with slow dynamics

Subthreshold (membrane potential) resonance and phasonance (preferred amplitude and zero-phase responses to oscillatory inputs) in single neurons arise from the interaction between positive and negative feedback effects provided by relatively fast amplifying currents and slower resonant currents. In 2D neuronal systems, amplifying currents are required to be slave to voltage (instantaneously fast) for these phenomena to occur. In higher dimensional systems, additional currents operating at various effective time scales may modulate and annihilate existing resonances and generate antiresonance (minimum amplitude response) and antiphasonance (zero-phase response with phase monotonic properties opposite to phasonance). We use mathematical modeling, numerical simulations and dynamical systems tools to investigate the mechanisms underlying these phenomena in 3D linear models, which are obtained as the linearization of biophysical (conductance-based) models. We characterize the parameter regimes for which the system exhibits the various types of behavior mentioned above in the rather general case in which the underlying 2D system exhibits resonance. We consider two cases: (i) the interplay of two resonant gating variables, and (ii) the interplay of one resonant and one amplifying gating variables. Increasing levels of an amplifying current cause (i) a response amplification if the amplifying current is faster than the resonant current, (ii) resonance and phasonance attenuation and annihilation if the amplifying and resonant currents have identical dynamics, and (iii) antiresonance and antiphasonance if the amplifying current is slower than the resonant current. We investigate the underlying mechanisms by extending the envelope-plane diagram approach developed in previous work (for 2D systems) to three dimensions to include the additional gating variable, and constructing the corresponding envelope curves in these envelope-space diagrams. We find that antiresonance and antiphasonance emerge as the result of an asymptotic boundary layer problem in the frequency domain created by the different balances between the intrinsic time constants of the cell and the input frequency f as it changes. For large enough values of f the envelope curves are quasi-2D and the impedance profile decreases with the input frequency. In contrast, for f ≪ 1 the dynamics are quasi-1D and the impedance profile increases above the limiting value in the other regime. Antiresonance is created because the continuity of the solution requires the impedance profile to connect the portions belonging to the two regimes. If in doing so the phase profile crosses the zero value, then antiphasonance is also generated.

Increased transient Na+ conductance and action potential output in layer 2/3 prefrontal cortex neurons of the fmr1-/y mouse

KEY POINTS

Layer 2/3 neurons of the prefrontal cortex display higher gain of somatic excitability, responding with a higher number of action potentials for a given stimulus, in fmr1^-/y mice. In fmr1^-/y L2/3 neurons, action potentials are taller, faster and narrower. Outside-out patch clamp recordings revealed that the maximum Na⁺ conductance density is higher in fmr1^-/y L2/3 neurons. Measurements of three biophysically distinct K⁺ currents revealed a depolarizing shift in the activation of a rapidly inactivating (A-type) K⁺ conductance. Realistic neuronal simulations of the biophysical observations recapitulated the elevated action potential and repetitive firing phenotype.

ABSTRACT

Fragile X syndrome is the most common form of inherited mental impairment and autism. The prefrontal cortex is responsible for higher order cognitive processing, and prefrontal dysfunction is believed to underlie many of the cognitive and behavioural phenotypes associated with fragile X syndrome. We recently demonstrated that somatic and dendritic excitability of layer (L) 5 pyramidal neurons in the prefrontal cortex of the fmr1^-/y mouse is significantly altered due to changes in several voltage-gated ion channels. In addition to L5 pyramidal neurons, L2/3 pyramidal neurons play an important role in prefrontal circuitry, integrating inputs from both lower brain regions and the contralateral cortex. Using whole-cell current clamp recording, we found that L2/3 pyramidal neurons in prefrontal cortex of fmr1^-/y mouse fired more action potentials for a given stimulus compared with wild-type neurons. In addition, action potentials in fmr1^-/y neurons were significantly larger, faster and narrower. Voltage clamp of outside-out patches from L2/3 neurons revealed that the transient Na⁺ current was significantly larger in fmr1^-/y neurons. Furthermore, the activation curve of somatic A-type K⁺ current was depolarized. Realistic conductance-based simulations revealed that these biophysical changes in Na⁺ and K⁺ channel function could reliably reproduce the observed increase in action potential firing and altered action potential waveform. These results, in conjunction with our prior findings on L5 neurons, suggest that principal neurons in the circuitry of the medial prefrontal cortex are altered in distinct ways in the fmr1^-/y mouse and may contribute to dysfunctional prefrontal cortex processing in fragile X syndrome.

Strings on a Violin: Location Dependence of Frequency Tuning in Active Dendrites

Strings on a violin are tuned to generate distinct sound frequencies in a manner that is firmly dependent on finger location along the fingerboard. Sound frequencies emerging from different violins could be very different based on their architecture, the nature of strings and their tuning. Analogously, active neuronal dendrites, dendrites endowed with active channel conductances, are tuned to distinct input frequencies in a manner that is dependent on the dendritic location of the synaptic inputs. Further, disparate channel expression profiles and differences in morphological characteristics could result in dendrites on different neurons of the same subtype tuned to distinct frequency ranges. Alternately, similar location-dependence along dendritic structures could be achieved through disparate combinations of channel profiles and morphological characteristics, leading to degeneracy in active dendritic spectral tuning. Akin to strings on a violin being tuned to different frequencies than those on a viola or a cello, different neuronal subtypes exhibit distinct channel profiles and disparate morphological characteristics endowing each neuronal subtype with unique location-dependent frequency selectivity. Finally, similar to the tunability of musical instruments to elicit distinct location-dependent sounds, neuronal frequency selectivity and its location-dependence are tunable through activity-dependent plasticity of ion channels and morphology. In this morceau, we explore the origins of neuronal frequency selectivity, and survey the literature on the mechanisms behind the emergence of location-dependence in distinct forms of frequency tuning. As a coda to this composition, we present some future directions for this exciting convergence of biophysical mechanisms that endow a neuron with frequency multiplexing capabilities.

Degeneracy in the regulation of short-term plasticity and synaptic filtering by presynaptic mechanisms

KEY POINTS

We develop a new biophysically rooted, physiologically constrained conductance-based synaptic model to mechanistically account for short-term facilitation and depression, respectively through residual calcium and transmitter depletion kinetics. We address the specific question of how presynaptic components (including voltage-gated ion channels, pumps, buffers and release-handling mechanisms) and interactions among them define synaptic filtering and short-term plasticity profiles. Employing global sensitivity analyses (GSAs), we show that near-identical synaptic filters and short-term plasticity profiles could emerge from disparate presynaptic parametric combinations with weak pairwise correlations. Using virtual knockout models, a technique to address the question of channel-specific contributions within the GSA framework, we unveil the differential and variable impact of each ion channel on synaptic physiology. Our conclusions strengthen the argument that parametric and interactional complexity in biological systems should not be viewed from the limited curse-of-dimensionality standpoint, but from the evolutionarily advantageous perspective of providing functional robustness through degeneracy.

ABSTRACT

Information processing in neurons is known to emerge as a gestalt of pre- and post-synaptic filtering. However, the impact of presynaptic mechanisms on synaptic filters has not been quantitatively assessed. Here, we developed a biophysically rooted, conductance-based model synapse that was endowed with six different voltage-gated ion channels, calcium pumps, calcium buffer and neurotransmitter-replenishment mechanisms in the presynaptic terminal. We tuned our model to match the short-term plasticity profile and band-pass structure of Schaffer collateral synapses, and performed sensitivity analyses to demonstrate that presynaptic voltage-gated ion channels regulated synaptic filters through changes in excitability and associated calcium influx. These sensitivity analyses also revealed that calcium- and release-control mechanisms were effective regulators of synaptic filters, but accomplished this without changes in terminal excitability or calcium influx. Next, to perform global sensitivity analysis, we generated 7000 randomized models spanning 15 presynaptic parameters, and computed eight different physiological measurements in each of these models. We validated these models by applying experimentally obtained bounds on their measurements, and found 104 (∼1.5%) models to match the validation criteria for all eight measurements. Analysing these valid models, we demonstrate that analogous synaptic filters emerge from disparate combinations of presynaptic parameters exhibiting weak pairwise correlations. Finally, using virtual knockout models, we establish the variable and differential impact of different presynaptic channels on synaptic filters, underlining the critical importance of interactions among different presynaptic components in defining synaptic physiology. Our results have significant implications for protein-localization strategies required for physiological robustness and for degeneracy in long-term synaptic plasticity profiles.

Active dendrites regulate the impact of gliotransmission on rat hippocampal pyramidal neurons

An important consequence of gliotransmission, a signaling mechanism that involves glial release of active transmitter molecules, is its manifestation as N-methyl-d-aspartate receptor (NMDAR)-dependent slow inward currents in neurons. However, the intraneuronal spatial dynamics of these events or the role of active dendrites in regulating their amplitude and spatial spread have remained unexplored. Here, we used somatic and/or dendritic recordings from rat hippocampal pyramidal neurons and demonstrate that a majority of NMDAR-dependent spontaneous slow excitatory potentials (SEP) originate at dendritic locations and are significantly attenuated through their propagation across the neuronal arbor. We substantiated the astrocytic origin of SEPs through paired neuron-astrocyte recordings, where we found that specific infusion of inositol trisphosphate (InsP3) into either distal or proximal astrocytes enhanced the amplitude and frequency of neuronal SEPs. Importantly, SEPs recorded after InsP3 infusion into distal astrocytes exhibited significantly slower kinetics compared with those recorded after proximal infusion. Furthermore, using neuron-specific infusion of pharmacological agents and morphologically realistic conductance-based computational models, we demonstrate that dendritically expressed hyperpolarization-activated cyclic-nucleotide-gated (HCN) and transient potassium channels play critical roles in regulating the strength, kinetics, and compartmentalization of neuronal SEPs. Finally, through the application of subtype-specific receptor blockers during paired neuron-astrocyte recordings, we provide evidence that GluN2B- and GluN2D-containing NMDARs predominantly mediate perisomatic and dendritic SEPs, respectively. Our results unveil an important role for active dendrites in regulating the impact of gliotransmission on neurons and suggest astrocytes as a source of dendritic plateau potentials that have been implicated in localized plasticity and place cell formation.

Transient potassium channels augment degeneracy in hippocampal active dendritic spectral tuning

Hippocampal pyramidal neurons express an intraneuronal map of spectral tuning mediated by hyperpolarization-activated cyclic-nucleotide-gated nonspecific-cation channels. Modeling studies have predicted a critical regulatory role for A-type potassium (KA) channels towards augmenting functional robustness of this map. To test this, we performed patch-clamp recordings from soma and dendrites of rat hippocampal pyramidal neurons, and measured spectral tuning before and after blocking KA channels using two structurally distinct pharmacological agents. Consistent with computational predictions, we found that blocking KA channels resulted in a significant reduction in resonance frequency and significant increases in input resistance, impedance amplitude and action-potential firing frequency across the somato-apical trunk. Furthermore, across all measured locations, blocking KA channels enhanced temporal summation of postsynaptic potentials and critically altered the impedance phase profile, resulting in a significant reduction in total inductive phase. Finally, pair-wise correlations between intraneuronal percentage changes (after blocking KA channels) in different measurements were mostly weak, suggesting differential regulation of different physiological properties by KA channels. Our results unveil a pivotal role for fast transient channels in regulating theta-frequency spectral tuning and intrinsic phase response, and suggest that degeneracy with reference to several coexisting functional maps is mediated by cross-channel interactions across the active dendritic arbor.

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.

Molecular variability elicits a tunable switch with discrete neuromodulatory response phenotypes

Recent single cell studies show extensive molecular variability underlying cellular responses. We evaluated the impact of molecular variability in the expression of cell signaling components and ion channels on electrophysiological excitability and neuromodulation. We employed a computational approach that integrated neuropeptide receptor-mediated signaling with electrophysiology. We simulated a population of neurons in which expression levels of a neuropeptide receptor and multiple ion channels were simultaneously varied within a physiological range. We analyzed the effects of variation on the electrophysiological response to a neuropeptide stimulus. Our results revealed distinct response patterns associated with low versus high receptor levels. Neurons with low receptor levels showed increased excitability and neurons with high receptor levels showed reduced excitability. These response patterns were separated by a narrow receptor level range forming a separatrix. The position of this separatrix was dependent on the expression levels of multiple ion channels. To assess the relative contributions of receptor and ion channel levels to the response profiles, we categorized the responses into six phenotypes based on response kinetics and magnitude. We applied several multivariate statistical approaches and found that receptor and channel expression levels influence the neuromodulation response phenotype through a complex though systematic mapping. Our analyses extended our understanding of how cellular responses to neuromodulation vary as a function of molecular expression. Our study showed that receptor expression and biophysical state interact with distinct relative contributions to neuronal excitability.

Variability in State-Dependent Plasticity of Intrinsic Properties during Cell-Autonomous Self-Regulation of Calcium Homeostasis in Hippocampal Model Neurons

How do neurons reconcile the maintenance of calcium homeostasis with perpetual switches in patterns of afferent activity? Here, we assessed state-dependent evolution of calcium homeostasis in a population of hippocampal pyramidal neuron models, through an adaptation of a recent study on stomatogastric ganglion neurons. Calcium homeostasis was set to emerge through cell-autonomous updates to 12 ionic conductances, responding to different types of synaptically driven afferent activity. We first assessed the impact of theta-frequency inputs on the evolution of ionic conductances toward maintenance of calcium homeostasis. Although calcium homeostasis emerged efficaciously across all models in the population, disparate changes in ionic conductances that mediated this emergence resulted in variable plasticity to several intrinsic properties, also manifesting as significant differences in firing responses across models. Assessing the sensitivity of this form of plasticity, we noted that intrinsic neuronal properties and the firing response were sensitive to the target calcium concentration and to the strength and frequency of afferent activity. Next, we studied the evolution of calcium homeostasis when afferent activity was switched, in different temporal sequences, between two behaviorally distinct types of activity: theta-frequency inputs and sharp-wave ripples riding on largely silent periods. We found that the conductance values, intrinsic properties, and firing response of neurons exhibited differential robustness to an intervening switch in the type of afferent activity. These results unveil critical dissociations between different forms of homeostasis, and call for a systematic evaluation of the impact of state-dependent switches in afferent activity on neuronal intrinsic properties during neural coding and homeostasis.

Active dendrites mediate stratified gamma-range coincidence detection in hippocampal model neurons

KEY POINTS

Quantitative metrics for the temporal window of integration/coincidence detection, based on the spike-triggered average, were employed to assess the emergence and dependence of gamma-range coincidence detection in hippocampal pyramidal neurons on various ion channel combinations. The presence of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels decreased the coincidence detection window (CDW) of the neuronal compartment to the gamma frequency range. Interaction of HCN channels with T-type calcium channels and persistent sodium channels further reduced the CDW, whereas interaction with A-type potassium channels broadened the CDW. When multiple channel gradients were co-expressed, the high density of resonating conductances in the distal dendrites led to a slow gamma CDW in the proximal dendrites and a fast-gamma CDW in the distal dendrites. The presence of resonating and spike-generating conductances serve as a mechanism underlying the emergence of stratified gamma-range coincidence detection in the dendrites of CA1 pyramidal neurons, enabling them to perform behaviour- and state-dependent gamma frequency multiplexing.

ABSTRACT

Hippocampal pyramidal neurons exhibit gamma-phase preference in their spikes, selectively route inputs through gamma frequency multiplexing and are considered part of gamma-bound cell assemblies. How do these neurons exhibit gamma-frequency coincidence detection capabilities, a feature that is essential for the expression of these physiological observations, despite their slow membrane time constant? In this conductance-based modelling study, we developed quantitative metrics for the temporal window of integration/coincidence detection based on the spike-triggered average (STA) of the neuronal compartment. We employed these metrics in conjunction with quantitative measures for spike initiation dynamics to assess the emergence and dependence of coincidence detection and STA spectral selectivity on various ion channel combinations. We found that the presence of resonating conductances (hyperpolarization-activated cyclic nucleotide-gated or T-type calcium), either independently or synergistically when expressed together, led to the emergence of spectral selectivity in the spike initiation dynamics and a significant reduction in the coincidence detection window (CDW). The presence of A-type potassium channels, along with resonating conductances, reduced the STA characteristic frequency and broadened the CDW, but persistent sodium channels sharpened the CDW by strengthening the spectral selectivity in the STA. Finally, in a morphologically precise model endowed with experimentally constrained channel gradients, we found that somatodendritic compartments expressed functional maps of strong theta-frequency selectivity in spike initiation dynamics and gamma-range CDW. Our results reveal the heavy expression of resonating and spike-generating conductances as the mechanism underlying the robust emergence of stratified gamma-range coincidence detection in the dendrites of hippocampal and cortical pyramidal neurons.

Analogous synaptic plasticity profiles emerge from disparate channel combinations

An open question within the Bienenstock-Cooper-Munro theory for synaptic modification concerns the specific mechanism that is responsible for regulating the sliding modification threshold (SMT). In this conductance-based modeling study on hippocampal pyramidal neurons, we quantitatively assessed the impact of seven ion channels (R- and T-type calcium, fast sodium, delayed rectifier, A-type, and small-conductance calcium-activated (SK) potassium and HCN) and two receptors (AMPAR and NMDAR) on a calcium-dependent Bienenstock-Cooper-Munro-like plasticity rule. Our analysis with R- and T-type calcium channels revealed that differences in their activation-inactivation profiles resulted in differential impacts on how they altered the SMT. Further, we found that the impact of SK channels on the SMT critically depended on the voltage dependence and kinetics of the calcium sources with which they interacted. Next, we considered interactions among all the seven channels and the two receptors through global sensitivity analysis on 11 model parameters. We constructed 20,000 models through uniform randomization of these parameters and found 360 valid models based on experimental constraints on their plasticity profiles. Analyzing these 360 models, we found that similar plasticity profiles could emerge with several nonunique parametric combinations and that parameters exhibited weak pairwise correlations. Finally, we used seven sets of virtual knock-outs on these 360 models and found that the impact of different channels on the SMT was variable and differential. These results suggest that there are several nonunique routes to regulate the SMT, and call for a systematic analysis of the variability and state dependence of the mechanisms underlying metaplasticity during behavior and pathology.

HCN channels enhance spike phase coherence and regulate the phase of spikes and LFPs in the theta-frequency range

What are the implications for the existence of subthreshold ion channels, their localization profiles, and plasticity on local field potentials (LFPs)? Here, we assessed the role of hyperpolarization-activated cyclic-nucleotide-gated (HCN) channels in altering hippocampal theta-frequency LFPs and the associated spike phase. We presented spatiotemporally randomized, balanced theta-modulated excitatory and inhibitory inputs to somatically aligned, morphologically realistic pyramidal neuron models spread across a cylindrical neuropil. We computed LFPs from seven electrode sites and found that the insertion of an experimentally constrained HCN-conductance gradient into these neurons introduced a location-dependent lead in the LFP phase without significantly altering its amplitude. Further, neurons fired action potentials at a specific theta phase of the LFP, and the insertion of HCN channels introduced large lags in this spike phase and a striking enhancement in neuronal spike-phase coherence. Importantly, graded changes in either HCN conductance or its half-maximal activation voltage resulted in graded changes in LFP and spike phases. Our conclusions on the impact of HCN channels on LFPs and spike phase were invariant to changes in neuropil size, to morphological heterogeneity, to excitatory or inhibitory synaptic scaling, and to shifts in the onset phase of inhibitory inputs. Finally, we selectively abolished the inductive lead in the impedance phase introduced by HCN channels without altering neuronal excitability and found that this inductive phase lead contributed significantly to changes in LFP and spike phase. Our results uncover specific roles for HCN channels and their plasticity in phase-coding schemas and in the formation and dynamic reconfiguration of neuronal cell assemblies.