Functional maps within a single neuron

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.

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.

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.

Conjunctive changes in multiple ion channels mediate activity-dependent intrinsic plasticity in hippocampal granule cells

Plasticity in the brain is ubiquitous. How do neurons and networks encode new information and simultaneously maintain homeostasis in the face of such ubiquitous plasticity? Here, we unveil a form of neuronal plasticity in rat hippocampal granule cells, which is mediated by conjunctive changes in HCN, inward-rectifier potassium, and persistent sodium channels induced by theta-modulated burst firing, a behaviorally relevant activity pattern. Cooperation and competition among these simultaneous changes resulted in a unique physiological signature: sub-threshold excitability and temporal summation were reduced without significant changes in action potential firing, together indicating a concurrent enhancement of supra-threshold excitability. This form of intrinsic plasticity was dependent on calcium influx through L-type calcium channels and inositol trisphosphate receptors. These observations demonstrate that although brain plasticity is ubiquitous, strong systemic constraints govern simultaneous plasticity in multiple components-referred here as plasticity manifolds-thereby providing a cellular substrate for concomitant encoding and homeostasis in engram cells.

Network instability dynamics drive a transient bursting period in the developing hippocampus in vivo

Spontaneous correlated activity is a universal hallmark of immature neural circuits. However, the cellular dynamics and intrinsic mechanisms underlying network burstiness in the intact developing brain are largely unknown. Here, we use two-photon Ca²⁺ imaging to comprehensively map the developmental trajectories of spontaneous network activity in the hippocampal area CA1 of mice in vivo. We unexpectedly find that network burstiness peaks after the developmental emergence of effective synaptic inhibition in the second postnatal week. We demonstrate that the enhanced network burstiness reflects an increased functional coupling of individual neurons to local population activity. However, pairwise neuronal correlations are low, and network bursts (NBs) recruit CA1 pyramidal cells in a virtually random manner. Using a dynamic systems modeling approach, we reconcile these experimental findings and identify network bi-stability as a potential regime underlying network burstiness at this age. Our analyses reveal an important role of synaptic input characteristics and network instability dynamics for NB generation. Collectively, our data suggest a mechanism, whereby developing CA1 performs extensive input-discrimination learning prior to the onset of environmental exploration.

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.

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.

Local Homeostatic Regulation of the Spectral Radius of Echo-State Networks

Recurrent cortical networks provide reservoirs of states that are thought to play a crucial role for sequential information processing in the brain. However, classical reservoir computing requires manual adjustments of global network parameters, particularly of the spectral radius of the recurrent synaptic weight matrix. It is hence not clear if the spectral radius is accessible to biological neural networks. Using random matrix theory, we show that the spectral radius is related to local properties of the neuronal dynamics whenever the overall dynamical state is only weakly correlated. This result allows us to introduce two local homeostatic synaptic scaling mechanisms, termed flow control and variance control, that implicitly drive the spectral radius toward the desired value. For both mechanisms the spectral radius is autonomously adapted while the network receives and processes inputs under working conditions. We demonstrate the effectiveness of the two adaptation mechanisms under different external input protocols. Moreover, we evaluated the network performance after adaptation by training the network to perform a time-delayed XOR operation on binary sequences. As our main result, we found that flow control reliably regulates the spectral radius for different types of input statistics. Precise tuning is however negatively affected when interneural correlations are substantial. Furthermore, we found a consistent task performance over a wide range of input strengths/variances. Variance control did however not yield the desired spectral radii with the same precision, being less consistent across different input strengths. Given the effectiveness and remarkably simple mathematical form of flow control, we conclude that self-consistent local control of the spectral radius via an implicit adaptation scheme is an interesting and biological plausible alternative to conventional methods using set point homeostatic feedback controls of neural firing.

Integration of Within-Cell Experimental Data With Multi-Compartmental Modeling Predicts H-Channel Densities and Distributions in Hippocampal OLM Cells

Determining biophysical details of spatially extended neurons is a challenge that needs to be overcome if we are to understand the dynamics of brain function from cellular perspectives. Moreover, we now know that we should not average across recordings from many cells of a given cell type to obtain quantitative measures such as conductance since measures can vary multiple-fold for a given cell type. In this work we examine whether a tight combination of experimental and computational work can address this challenge. The oriens-lacunosum/moleculare (OLM) interneuron operates as a “gate” that controls incoming sensory and ongoing contextual information in the CA1 of the hippocampus, making it essential to understand how its biophysical properties contribute to memory function. OLM cells fire phase-locked to the prominent hippocampal theta rhythms, and we previously used computational models to show that OLM cells exhibit high or low theta spiking resonance frequencies that depend respectively on whether their dendrites have hyperpolarization-activated cation channels (h-channels) or not. However, whether OLM cells actually possess dendritic h-channels is unknown at present. We performed a set of whole-cell recordings of OLM cells from mouse hippocampus and constructed three multi-compartment models using morphological and electrophysiological parameters extracted from the same OLM cell, including per-cell pharmacologically isolated h-channel currents. We found that the models best matched experiments when h-channels were present in the dendrites of each of the three model cells created. This strongly suggests that h-channels must be present in OLM cell dendrites and are not localized to their somata. Importantly, this work shows that a tight integration of model and experiment can help tackle the challenge of characterizing biophysical details and distributions in spatially extended neurons. Full spiking models were built for two of the OLM cells, matching their current clamp cell-specific electrophysiological recordings. Overall, our work presents a technical advancement in modeling OLM cells. Our models are available to the community to use to gain insight into cellular dynamics underlying hippocampal function.

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.

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.

Stores, Channels, Glue, and Trees: Active Glial and Active Dendritic Physiology

Glial cells and neuronal dendrites were historically assumed to be passive structures that play only supportive physiological roles, with no active contribution to information processing in the central nervous system. Research spanning the past few decades has clearly established this assumption to be far from physiological realities. Whereas the discovery of active channel conductances and their localized plasticity was the turning point for dendritic structures, the demonstration that glial cells release transmitter molecules and communicate across the neuroglia syncytium through calcium wave propagation constituted path-breaking discoveries for glial cell physiology. An additional commonality between these two structures is the ability of calcium stores within their endoplasmic reticulum (ER) to support active propagation of calcium waves, which play crucial roles in the spatiotemporal integration of information within and across cells. Although there have been several demonstrations of regulatory roles of glial cells and dendritic structures in achieving common physiological goals such as information propagation and adaptability through plasticity, studies assessing physiological interactions between these two active structures have been few and far. This lacuna is especially striking given the strong connectivity that is known to exist between these two structures through several complex and tightly intercoupled mechanisms that also recruit their respective ER structures. In this review, we present brief overviews of the parallel literatures on active dendrites and active glial physiology and make a strong case for future studies to directly assess the strong interactions between these two structures in regulating physiology and pathophysiology of the brain.

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.

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.

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).

Modulating STDP Balance Impacts the Dendritic Mosaic

The ability for cortical neurons to adapt their input/output characteristics and information processing capabilities ultimately relies on the interplay between synaptic plasticity, synapse location, and the nonlinear properties of the dendrite. Collectively, they shape both the strengths and spatial arrangements of convergent afferent inputs to neuronal dendrites. Recent experimental and theoretical studies support a clustered plasticity model, a view that synaptic plasticity promotes the formation of clusters or hotspots of synapses sharing similar properties. We have previously shown that spike timing-dependent plasticity (STDP) can lead to synaptic efficacies being arranged into spatially segregated clusters. This effectively partitions the dendritic tree into a tessellated imprint which we have called a dendritic mosaic. Here, using a biophysically detailed neuron model of a reconstructed layer 2/3 pyramidal cell and STDP learning, we investigated the impact of altered STDP balance on forming such a spatial organization. We show that cluster formation and extend depend on several factors, including the balance between potentiation and depression, the afferents' mean firing rate and crucially on the dendritic morphology. We find that STDP balance has an important role to play for this emergent mode of spatial organization since any imbalances lead to severe degradation- and in some case even destruction- of the mosaic. Our model suggests that, over a broad range of of STDP parameters, synaptic plasticity shapes the spatial arrangement of synapses, favoring the formation of clustered efficacy engrams.

Synaptic plasticity in dendrites: complications and coping strategies

The elaborate morphology, nonlinear membrane mechanisms and spatiotemporally varying synaptic activation patterns of dendrites complicate the expression, compartmentalization and modulation of synaptic plasticity. To grapple with this complexity, we start with the observation that neurons in different brain areas face markedly different learning problems, and dendrites of different neuron types contribute to the cell's input-output function in markedly different ways. By committing to specific assumptions regarding a neuron's learning problem and its input-output function, specific inferences can be drawn regarding the synaptic plasticity mechanisms and outcomes that we 'ought' to expect for that neuron. Exploiting this assumption-driven approach can help both in interpreting existing experimental data and designing future experiments aimed at understanding the brain's myriad learning processes.

The Segregated Expression of Voltage-Gated Potassium and Sodium Channels in Neuronal Membranes: Functional Implications and Regulatory Mechanisms

Neurons are highly polarized cells with apparent functional and morphological differences between dendrites and axon. A critical determinant for the molecular and functional identity of axonal and dendritic segments is the restricted expression of voltage-gated ion channels (VGCs). Several studies show an uneven distribution of ion channels and their differential regulation within dendrites and axons, which is a prerequisite for an appropriate integration of synaptic inputs and the generation of adequate action potential (AP) firing patterns. This review article will focus on the signaling pathways leading to segmented expression of voltage-gated potassium and sodium ion channels at the neuronal plasma membrane and the regulatory mechanisms ensuring segregated functions. We will also discuss the relevance of proper ion channel targeting for neuronal physiology and how alterations in polarized distribution contribute to neuronal pathology.

Computational models of O-LM cells are recruited by low or high theta frequency inputs depending on h-channel distributions

Although biophysical details of inhibitory neurons are becoming known, it is challenging to map these details onto function. Oriens-lacunosum/moleculare (O-LM) cells are inhibitory cells in the hippocampus that gate information flow, firing while phase-locked to theta rhythms. We build on our existing computational model database of O-LM cells to link model with function. We place our models in high-conductance states and modulate inhibitory inputs at a wide range of frequencies. We find preferred spiking recruitment of models at high (4-9 Hz) or low (2-5 Hz) theta depending on, respectively, the presence or absence of h-channels on their dendrites. This also depends on slow delayed-rectifier potassium channels, and preferred theta ranges shift when h-channels are potentiated by cyclic AMP. Our results suggest that O-LM cells can be differentially recruited by frequency-modulated inputs depending on specific channel types and distributions. This work exposes a strategy for understanding how biophysical characteristics contribute to function.

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.

Reduced Glutamatergic Currents and Dendritic Branching of Layer 5 Pyramidal Cells Contribute to Medial Prefrontal Cortex Deactivation in a Rat Model of Neuropathic Pain

Multiple studies have demonstrated that neuropathic pain is associated with major reorganization in multiple brain areas. In line with the strong emotional salience of chronic pain, involvement of the limbic system appears particularly important. Within the past few years, it has become clear that the functional deactivation of the prefrontal cortex (PFC) is critical for both the cognitive/emotional and the sensory components of pain. However, at the cellular level, details of this deactivation remain in large part unclear. Here we show that 1 week after a peripheral neuropathic injury (Spared Nerve Injury model) pyramidal cells in layer 5 (L5) of the rat medial PFC show responses to excitatory glutamatergic inputs that are reduced by about 50%, as well as reduced frequency of spontaneous excitatory synaptic currents. Additionally, these cells have reduced membrane capacitance and increased input resistance. All these findings are consistent with decreased dendritic length, thus we performed a detailed morphological analysis on a subset of the recorded neurons. We found that the apical dendrites proximal to the soma (excluding the tuft) are shorter and less complex in SNI animals, in agreement with the reduced capacitance and glutamatergic input. Finally, we used in vivo microdialysis to compare the basal concentrations of glutamate and GABA in the PFC of sham and SNI rats and found that ambient glutamate is decreased in SNI rats. Taken together, these data show that impaired glutamatergic transmission contributes to the functional deactivation of the mPFC in neuropathic pain. Additionally, the reduced branching of apical dendrites of L5 pyramidal neurons may underlay the gray matter reduction in chronic pain.

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.

Visual Stimuli Evoked Action Potentials Trigger Rapidly Propagating Dendritic Calcium Transients in the Frog Optic Tectum Layer 6 Neurons

The superior colliculus in mammals or the optic tectum in amphibians is a major visual information processing center responsible for generation of orientating responses such as saccades in monkeys or prey catching avoidance behavior in frogs. The conserved structure function of the superior colliculus the optic tectum across distant species such as frogs, birds monkeys permits to draw rather general conclusions after studying a single species. We chose the frog optic tectum because we are able to perform whole-cell voltage-clamp recordings fluorescence imaging of tectal neurons while they respond to a visual stimulus. In the optic tectum of amphibians most visual information is processed by pear-shaped neurons possessing long dendritic branches, which receive the majority of synapses originating from the retinal ganglion cells. Since the first step of the retinal input integration is performed on these dendrites, it is important to know whether this integration is enhanced by active dendritic properties. We demonstrate that rapid calcium transients coinciding with the visual stimulus evoked action potentials in the somatic recordings can be readily detected up to the fine branches of these dendrites. These transients were blocked by calcium channel blockers nifedipine CdCl₂ indicating that calcium entered dendrites via voltage-activated L-type calcium channels. The high speed of calcium transient propagation, >300 μm in <10 ms, is consistent with the notion that action potentials, actively propagating along dendrites, open voltage-gated L-type calcium channels causing rapid calcium concentration transients in the dendrites. We conclude that such activation by somatic action potentials of the dendritic voltage gated calcium channels in the close vicinity to the synapses formed by axons of the retinal ganglion cells may facilitate visual information processing in the principal neurons of the frog optic tectum.

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.

Branch specific and spike-order specific action potential invasion in basal, oblique, and apical dendrites of cortical pyramidal neurons

In neocortical pyramidal neurons, action potentials (APs) propagate from the axon into the dendritic tree to influence distal synapses. Traditionally, AP backpropagation was studied in the thick apical trunk. Here, we used the principles of optical imaging developed by Cohen to investigate AP invasion into thin dendritic branches (basal, oblique, and tuft) of prefrontal cortical L5 pyramidal neurons. Multisite optical recordings from neighboring dendrites revealed a clear dichotomy between two seemingly equal dendritic branches belonging to the same cell ("sister branches"). We documented the variable efficacy of AP invasion in basal and oblique branches by revealing their AP voltage waveforms. Using fast multisite calcium imaging, we found that trains of APs are filtered differently between two apical tuft branches. Although one dendritic branch passes all spikes in an AP train, another branch belonging to the same neuron, same cortical layer, and same path distance from the cell body, experiences only one spike. Our data indicate that the vast differences in dendritic voltage and calcium transients, detected in dendrites of pyramidal neurons, arise from a nonuniform distribution of A-type [Formula: see text] conductance, an aggregate number of branch points in the path of the AP propagation and minute differences in dendritic diameter.

Dendritic atrophy constricts functional maps in resonance and impedance properties of hippocampal model neurons

A gradient in the density of hyperpolarization-activated cyclic-nucleotide gated (HCN) channels is necessary for the emergence of several functional maps within hippocampal pyramidal neurons. Here, we systematically analyzed the impact of dendritic atrophy on nine such functional maps, related to input resistance and local/transfer impedance properties, using conductance-based models of hippocampal pyramidal neurons. We introduced progressive dendritic atrophy in a CA1 pyramidal neuron reconstruction through a pruning algorithm, measured all functional maps in each pruned reconstruction, and arrived at functional forms for the dependence of underlying measurements on dendritic length. We found that, across frequencies, atrophied neurons responded with higher efficiency to incoming inputs, and the transfer of signals across the dendritic tree was more effective in an atrophied reconstruction. Importantly, despite the presence of identical HCN-channel density gradients, spatial gradients in input resistance, local/transfer resonance frequencies and impedance profiles were significantly constricted in reconstructions with dendritic atrophy, where these physiological measurements across dendritic locations converged to similar values. These results revealed that, in atrophied dendritic structures, the presence of an ion channel density gradient alone was insufficient to sustain homologous functional maps along the same neuronal topograph. We assessed the biophysical basis for these conclusions and found that this atrophy-induced constriction of functional maps was mediated by an enhanced spatial spread of the influence of an HCN-channel cluster in atrophied trees. These results demonstrated that the influence fields of ion channel conductances need to be localized for channel gradients to express themselves as homologous functional maps, suggesting that ion channel gradients are necessary but not sufficient for the emergence of functional maps within single neurons.

High-conductance states and A-type K+ channels are potential regulators of the conductance-current balance triggered by HCN channels

An increase in the hyperpolarization-activated cyclic nucleotide-gated (HCN) channel conductance reduces input resistance, whereas the consequent increase in the inward h current depolarizes the membrane. This results in a delicate and unique conductance-current balance triggered by the expression of HCN channels. In this study, we employ experimentally constrained, morphologically realistic, conductance-based models of hippocampal neurons to explore certain aspects of this conductance-current balance. First, we found that the inclusion of an experimentally determined gradient in A-type K+ conductance, but not in M-type K+ conductance, tilts the HCN conductance-current balance heavily in favor of conductance, thereby exerting an overall restorative influence on neural excitability. Next, motivated by the well-established modulation of neuronal excitability by synaptically driven high-conductance states observed under in vivo conditions, we inserted thousands of excitatory and inhibitory synapses with different somatodendritic distributions. We measured the efficacy of HCN channels, independently and in conjunction with other channels, in altering resting membrane potential (RMP) and input resistance ( Rin) when the neuron received randomized or rhythmic synaptic bombardments through variable numbers of synaptic inputs. We found that the impact of HCN channels on average RMP, Rin, firing frequency, and peak-to-peak voltage response was severely weakened under high-conductance states, with the impinging synaptic drive playing a dominant role in regulating these measurements. Our results suggest that the debate on the role of HCN channels in altering excitability should encompass physiological and pathophysiological neuronal states under in vivo conditions and the spatiotemporal interactions of HCN channels with other channels.

Activation of InsP₃ receptors is sufficient for inducing graded intrinsic plasticity in rat hippocampal pyramidal neurons

The synaptic plasticity literature has focused on establishing necessity and sufficiency as two essential and distinct features in causally relating a signaling molecule to plasticity induction, an approach that has been surprisingly lacking in the intrinsic plasticity literature. In this study, we complemented the recently established necessity of inositol trisphosphate (InsP3) receptors (InsP3R) in a form of intrinsic plasticity by asking if InsP3R activation was sufficient to induce intrinsic plasticity in hippocampal neurons. Specifically, incorporation of d-myo-InsP3 in the recording pipette reduced input resistance, maximal impedance amplitude, and temporal summation but increased resonance frequency, resonance strength, sag ratio, and impedance phase lead. Strikingly, the magnitude of plasticity in all these measurements was dependent on InsP3 concentration, emphasizing the graded dependence of such plasticity on InsP3R activation. Mechanistically, we found that this InsP3-induced plasticity depended on hyperpolarization-activated cyclic nucleotide-gated channels. Moreover, this calcium-dependent form of plasticity was critically reliant on the release of calcium through InsP3Rs, the influx of calcium through N-methyl-d-aspartate receptors and voltage-gated calcium channels, and on the protein kinase A pathway. Our results delineate a causal role for InsP3Rs in graded adaptation of neuronal response dynamics, revealing novel regulatory roles for the endoplasmic reticulum in neural coding and homeostasis.

Frequency-dependent signal processing in apical dendrites of hippocampal CA1 pyramidal cells

Depending on an animal's behavioral state, hippocampal CA1 pyramidal cells receive distinct patterns of excitatory and inhibitory synaptic inputs. The time-dependent changes in the frequencies of these inputs and the nonuniform distribution of voltage-gated channels lead to dynamic fluctuations in membrane conductance. In this study, using a whole-cell patch-clamp method, we attempted to record and analyze the frequency dependencies of membrane responsiveness in Wistar rat hippocampal CA1 pyramidal cells following noise current injection directly into dendrites and somata under pharmacological blockade of all synaptic inputs. To estimate the frequency-dependent properties of membrane potential, membrane impedance was determined from the voltage response divided by the input current in the frequency domain. The cell membrane of most neurons showed low-pass filtering properties in all regions. In particular, the properties were strongly expressed in the somata or proximal dendrites. Moreover, the data revealed nonuniform distribution of dendritic impedance, which was high in the intermediate segment of the apical dendritic shaft (∼220-260μm from the soma). The low-pass filtering properties in the apical dendrites were more enhanced by membrane depolarization than those in the somata. Coherence spectral analysis revealed high coherence between the input signal and the output voltage response in the theta-gamma frequency range, and large lags emerged in the distal dendrites in the gamma frequency range. Our results suggest that apical dendrites of hippocampal CA1 pyramidal cells integrate synaptic inputs according to the frequency components of the input signal along the dendritic segments receiving the inputs.

An Augmented Two-Layer Model Captures Nonlinear Analog Spatial Integration Effects in Pyramidal Neuron Dendrites

In pursuit of the goal to understand and eventually reproduce the diverse functions of the brain, a key challenge lies in reverse engineering the peculiar biology-based "technology" that underlies the brain's remarkable ability to process and store information. The basic building block of the nervous system is the nerve cell, or "neuron," yet after more than 100 years of neurophysiological study and 60 years of modeling, the information processing functions of individual neurons, and the parameters that allow them to engage in so many different types of computation (sensory, motor, mnemonic, executive, etc.) remain poorly understood. In this paper, we review both historical and recent findings that have led to our current understanding of the analog spatial processing capabilities of dendrites, the major input structures of neurons, with a focus on the principal cell type of the neocortex and hippocampus, the pyramidal neuron (PN). We encapsulate our current understanding of PN dendritic integration in an abstract layered model whose spatially sensitive branch-subunits compute multidimensional sigmoidal functions. Unlike the 1-D sigmoids found in conventional neural network models, multidimensional sigmoids allow the cell to implement a rich spectrum of nonlinear modulation effects directly within their dendritic trees.

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

The maintenance of ion channel homeostasis, or channelostasis, is a complex puzzle in neurons with extensive dendritic arborization, encompassing a combinatorial diversity of proteins that encode these channels and their auxiliary subunits, their localization profiles, and associated signaling machinery. Despite this, neurons exhibit amazingly stereotypic, topographically continuous maps of several functional properties along their active dendritic arbor. Here, we asked whether the membrane composition of neurons, at the level of individual ion channels, is constrained by this structural requirement of sustaining several functional maps along the same topograph. We performed global sensitivity analysis on morphologically realistic conductance-based models of hippocampal pyramidal neurons that coexpressed six well-characterized functional maps along their trunk. We generated randomized models by varying 32 underlying parameters and constrained these models with quantitative experimental measurements from the soma and dendrites of hippocampal pyramidal neurons. Analyzing valid models that satisfied experimental constraints on all six functional maps, we found topographically analogous functional maps to emerge from disparate model parameters with weak pairwise correlations between parameters. Finally, we derived a methodology to assess the contribution of individual channel conductances to the various functional measurements, using virtual knockout simulations on the valid model population. We found that the virtual knockout of individual channels resulted in variable, measurement- and location-specific impacts across the population. Our results suggest collective channelostasis as a mechanism behind the robust emergence of analogous functional maps and have significant ramifications for the localization and targeting of ion channels and enzymes that regulate neural coding and homeostasis.

Active dendrites regulate spectral selectivity in location-dependent spike initiation dynamics of hippocampal model neurons

How does the presence of plastic active dendrites in a pyramidal neuron alter its spike initiation dynamics? To answer this question, we measured the spike-triggered average (STA) from experimentally constrained, conductance-based hippocampal neuronal models of various morphological complexities. We transformed the STA computed from these models to the spectral and the spectrotemporal domains and found that the spike initiation dynamics exhibited temporally localized selectivity to a characteristic frequency. In the presence of the hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, the STA characteristic frequency strongly correlated with the subthreshold resonance frequency in the theta frequency range. Increases in HCN channel density or in input variance increased the STA characteristic frequency and its selectivity strength. In the absence of HCN channels, the STA exhibited weak delta frequency selectivity and the characteristic frequency was related to the repolarization dynamics of the action potentials and the recovery kinetics of sodium channels from inactivation. Comparison of STA obtained with inputs at various dendritic locations revealed that nonspiking and spiking dendrites increased and reduced the spectrotemporal integration window of the STA with increasing distance from the soma as direct consequences of passive filtering and dendritic spike initiation, respectively. Finally, the presence of HCN channels set the STA characteristic frequency in the theta range across the somatodendritic arbor and specific STA measurements were strongly related to equivalent transfer-impedance-related measurements. Our results identify explicit roles for plastic active dendrites in neural coding and strongly recommend a dynamically reconfigurable multi-STA model to characterize location-dependent input feature selectivity in pyramidal neurons.

The log-dynamic brain: how skewed distributions affect network operations

We often assume that the variables of functional and structural brain parameters - such as synaptic weights, the firing rates of individual neurons, the synchronous discharge of neural populations, the number of synaptic contacts between neurons and the size of dendritic boutons - have a bell-shaped distribution. However, at many physiological and anatomical levels in the brain, the distribution of numerous parameters is in fact strongly skewed with a heavy tail, suggesting that skewed (typically lognormal) distributions are fundamental to structural and functional brain organization. This insight not only has implications for how we should collect and analyse data, it may also help us to understand how the different levels of skewed distributions - from synapses to cognition - are related to each other.

Preconfigured, skewed distribution of firing rates in the hippocampus and entorhinal cortex

Despite the importance of the discharge frequency in neuronal communication, little is known about the firing-rate patterns of cortical populations. Using large-scale recordings from multiple layers of the entorhinal-hippocampal loop, we found that the firing rates of principal neurons showed a lognormal-like distribution in all brain states. Mean and peak rates within place fields of hippocampal neurons were also strongly skewed. Importantly, firing rates of the same neurons showed reliable correlations in different brain states and testing situations, as well as across familiar and novel environments. The fraction of neurons that participated in population oscillations displayed a lognormal pattern. Such skewed firing rates of individual neurons may be due to a skewed distribution of synaptic weights, which is supported by our observation of a lognormal distribution of the efficacy of spike transfer from principal neurons to interneurons. The persistent skewed distribution of firing rates implies that a preconfigured, highly active minority dominates information transmission in cortical networks.

A calcium-dependent plasticity rule for HCN channels maintains activity homeostasis and stable synaptic learning

Theoretical and computational frameworks for synaptic plasticity and learning have a long and cherished history, with few parallels within the well-established literature for plasticity of voltage-gated ion channels. In this study, we derive rules for plasticity in the hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, and assess the synergy between synaptic and HCN channel plasticity in establishing stability during synaptic learning. To do this, we employ a conductance-based model for the hippocampal pyramidal neuron, and incorporate synaptic plasticity through the well-established Bienenstock-Cooper-Munro (BCM)-like rule for synaptic plasticity, wherein the direction and strength of the plasticity is dependent on the concentration of calcium influx. Under this framework, we derive a rule for HCN channel plasticity to establish homeostasis in synaptically-driven firing rate, and incorporate such plasticity into our model. In demonstrating that this rule for HCN channel plasticity helps maintain firing rate homeostasis after bidirectional synaptic plasticity, we observe a linear relationship between synaptic plasticity and HCN channel plasticity for maintaining firing rate homeostasis. Motivated by this linear relationship, we derive a calcium-dependent rule for HCN-channel plasticity, and demonstrate that firing rate homeostasis is maintained in the face of synaptic plasticity when moderate and high levels of cytosolic calcium influx induced depression and potentiation of the HCN-channel conductance, respectively. Additionally, we show that such synergy between synaptic and HCN-channel plasticity enhances the stability of synaptic learning through metaplasticity in the BCM-like synaptic plasticity profile. Finally, we demonstrate that the synergistic interaction between synaptic and HCN-channel plasticity preserves robustness of information transfer across the neuron under a rate-coding schema. Our results establish specific physiological roles for experimentally observed plasticity in HCN channels accompanying synaptic plasticity in hippocampal neurons, and uncover potential links between HCN-channel plasticity and calcium influx, dynamic gain control and stable synaptic learning.