Mechanisms of firing patterns in fast-spiking cortical interneurons

Basolateral amygdala oscillations enable fear learning in a biophysical model

The basolateral amygdala (BLA) is a key site where fear learning takes place through synaptic plasticity. Rodent research shows prominent low theta (~3-6 Hz), high theta (~6-12 Hz), and gamma (>30 Hz) rhythms in the BLA local field potential recordings. However, it is not understood what role these rhythms play in supporting the plasticity. Here, we create a biophysically detailed model of the BLA circuit to show that several classes of interneurons (PV, SOM, and VIP) in the BLA can be critically involved in producing the rhythms; these rhythms promote the formation of a dedicated fear circuit shaped through spike-timing-dependent plasticity. Each class of interneurons is necessary for the plasticity. We find that the low theta rhythm is a biomarker of successful fear conditioning. The model makes use of interneurons commonly found in the cortex and, hence, may apply to a wide variety of associative learning situations.

Existence and Stability Criteria for Global Synchrony and for Synchrony in two Alternating Clusters of Pulse-Coupled Oscillators Updated to Include Conduction Delays

Phase Response Curves (PRCs) have been useful in determining and analyzing various phase-locking modes in networks of oscillators under pulse-coupling assumptions, as reviewed in Mathematical Biosciences, 226:77-96, 2010. Here, we update that review to include progress since 2010 on pulse coupled oscillators with conduction delays. We then present original results that extend the derivation of the criteria for stability of global synchrony in networks of pulse-coupled oscillators to include conduction delays. We also incorporate conduction delays to extend previous studies that showed how an alternating firing pattern between two synchronized clusters could enforce within cluster synchrony, even for clusters unable to synchronize themselves in isolation. To obtain these results, we used self-connected neurons to represent clusters. These results greatly extend the applicability of the stability analyses to networks of pulse-coupled oscillators since conduction delays are ubiquitous and strongly impact the stability of synchrony. Although these analyses only strictly apply to identical oscillators with identical connections to other oscillators, the principles are general and suggest how to promote or impede synchrony in physiological networks of neurons, for example. Heterogeneity can be interpreted as a form of frozen noise, and approximate synchrony can be sustained despite heterogeneity. The pulse-coupled oscillator model can not only be used to describe biological neuronal networks but also cardiac pacemakers, lasers, fireflies, artificial neural networks, social self-organization, and wireless sensor networks.

AMS Subject Classification

37N25, 39A06, 39A30, 92B25, 92C20.

Utilising activity patterns of a complex biophysical network model to optimise intra-striatal deep brain stimulation

A large-scale biophysical network model for the isolated striatal body is developed to optimise potential intrastriatal deep brain stimulation applied to, e.g. obsessive-compulsive disorder. The model is based on modified Hodgkin-Huxley equations with small-world connectivity, while the spatial information about the positions of the neurons is taken from a detailed human atlas. The model produces neuronal spatiotemporal activity patterns segregating healthy from pathological conditions. Three biomarkers were used for the optimisation of stimulation protocols regarding stimulation frequency, amplitude and localisation: the mean activity of the entire network, the frequency spectrum of the entire network (rhythmicity) and a combination of the above two. By minimising the deviation of the aforementioned biomarkers from the normal state, we compute the optimal deep brain stimulation parameters, regarding position, amplitude and frequency. Our results suggest that in the DBS optimisation process, there is a clear trade-off between frequency synchronisation and overall network activity, which has also been observed during in vivo studies.

Data-driven and equation-free methods for neurological disorders: analysis and control of the striatum network

The striatum as part of the basal ganglia is central to both motor, and cognitive functions. Here, we propose a large-scale biophysical network for this part of the brain, using modified Hodgkin-Huxley dynamics to model neurons, and a connectivity informed by a detailed human atlas. The model shows different spatio-temporal activity patterns corresponding to lower (presumably normal) and increased cortico-striatal activation (as found in, e.g., obsessive-compulsive disorder), depending on the intensity of the cortical inputs. By applying equation-free methods, we are able to perform a macroscopic network analysis directly from microscale simulations. We identify the mean synaptic activity as the macroscopic variable of the system, which shows similarity with local field potentials. The equation-free approach results in a numerical bifurcation and stability analysis of the macroscopic dynamics of the striatal network. The different macroscopic states can be assigned to normal/healthy and pathological conditions, as known from neurological disorders. Finally, guided by the equation-free bifurcation analysis, we propose a therapeutic close loop control scheme for the striatal network.

Persistent Interruption in Parvalbumin-Positive Inhibitory Interneurons: Biophysical and Mathematical Mechanisms

Persistent activity in excitatory pyramidal cells (PYRs) is a putative mechanism for maintaining memory traces during working memory. We have recently demonstrated persistent interruption of firing in fast-spiking parvalbumin-expressing interneurons (PV-INs), a phenomenon that could serve as a substrate for persistent activity in PYRs through disinhibition lasting hundreds of milliseconds. Here, we find that hippocampal CA1 PV-INs exhibit type 2 excitability, like striatal and neocortical PV-INs. Modeling and mathematical analysis showed that the slowly inactivating potassium current KV1 contributes to type 2 excitability, enables the multiple firing regimes observed experimentally in PV-INs, and provides a mechanism for robust persistent interruption of firing. Using a fast/slow separation of times scales approach with the KV1 inactivation variable as a bifurcation parameter shows that the initial inhibitory stimulus stops repetitive firing by moving the membrane potential trajectory onto a coexisting stable fixed point corresponding to a nonspiking quiescent state. As KV1 inactivation decays, the trajectory follows the branch of stable fixed points until it crosses a subcritical Hopf bifurcation (HB) and then spirals out into repetitive firing. In a model describing entorhinal cortical PV-INs without KV1, interruption of firing could be achieved by taking advantage of the bistability inherent in type 2 excitability based on a subcritical HB, but the interruption was not robust to noise. Persistent interruption of firing is therefore broadly applicable to PV-INs in different brain regions but is only made robust to noise in the presence of a slow variable, KV1 inactivation.

Ketamine can produce oscillatory dynamics by engaging mechanisms dependent on the kinetics of NMDA receptors

Ketamine is an N-methyl-D-aspartate (NMDA)-receptor antagonist that produces sedation, analgesia, and dissociation at low doses and profound unconsciousness with antinociception at high doses. At high and low doses, ketamine can generate gamma oscillations (>25 Hz) in the electroencephalogram (EEG). The gamma oscillations are interrupted by slow-delta oscillations (0.1 to 4 Hz) at high doses. Ketamine's primary molecular targets and its oscillatory dynamics have been characterized. However, how the actions of ketamine at the subcellular level give rise to the oscillatory dynamics observed at the network level remains unknown. By developing a biophysical model of cortical circuits, we demonstrate how NMDA-receptor antagonism by ketamine can produce the oscillatory dynamics observed in human EEG recordings and nonhuman primate local field potential recordings. We have identified how impaired NMDA-receptor kinetics can cause disinhibition in neuronal circuits and how a disinhibited interaction between NMDA-receptor-mediated excitation and GABA-receptor-mediated inhibition can produce gamma oscillations at high and low doses, and slow-delta oscillations at high doses. Our work uncovers general mechanisms for generating oscillatory brain dynamics that differs from ones previously reported and provides important insights into ketamine's mechanisms of action as an anesthetic and as a therapy for treatment-resistant depression.

Parvalbumin interneuron activity in autism underlies susceptibility to PTSD-like memory formation

A rising concern in autism spectrum disorder (ASD) is the heightened sensitivity to trauma, the potential consequences of which have been overlooked, particularly upon the severity of the ASD traits. We first demonstrate a reciprocal relationship between ASD and post-traumatic stress disorder (PTSD) and reveal that exposure to a mildly stressful event induces PTSD-like memory in four mouse models of ASD. We also establish an unanticipated consequence of stress, as the formation of PTSD-like memory leads to the aggravation of core autistic traits. Such a susceptibility to developing PTSD-like memory in ASD stems from hyperactivation of the prefrontal cortex and altered fine-tuning of parvalbumin interneuron firing. Traumatic memory can be treated by recontextualization, reducing the deleterious effects on the core symptoms of ASD in the Cntnap2 KO mouse model. This study provides a neurobiological and psychological framework for future examination of the impact of PTSD-like memory in autism.

Ketamine can produce oscillatory dynamics by engaging mechanisms dependent on the kinetics of NMDA receptors

Ketamine is an NMDA-receptor antagonist that produces sedation, analgesia and dissociation at low doses and profound unconsciousness with antinociception at high doses. At high and low doses, ketamine can generate gamma oscillations (>25 Hz) in the electroencephalogram (EEG). The gamma oscillations are interrupted by slow-delta oscillations (0.1-4 Hz) at high doses. Ketamine's primary molecular targets and its oscillatory dynamics have been characterized. However, how the actions of ketamine at the subcellular level give rise to the oscillatory dynamics observed at the network level remains unknown. By developing a biophysical model of cortical circuits, we demonstrate how NMDA-receptor antagonism by ketamine can produce the oscillatory dynamics observed in human EEG recordings and non-human primate local field potential recordings. We have discovered how impaired NMDA-receptor kinetics can cause disinhibition in neuronal circuits and how a disinhibited interaction between NMDA-receptor-mediated excitation and GABA-receptor-mediated inhibition can produce gamma oscillations at high and low doses, and slow-delta oscillations at high doses. Our work uncovers general mechanisms for generating oscillatory brain dynamics that differs from ones previously reported, and provides important insights into ketamine's mechanisms of action as an anesthetic and as a therapy for treatment-resistant depression.

Persistent Interruption in Parvalbumin Positive Inhibitory Interneurons: Biophysical and Mathematical Mechanisms

Persistent activity in principal cells is a putative mechanism for maintaining memory traces during working memory. We recently demonstrated persistent interruption of firing in fast-spiking parvalbumin-expressing interneurons (PV-INs), a phenomenon which could serve as a substrate for persistent activity in principal cells through disinhibition lasting hundreds of milliseconds. Here, we find that hippocampal CA1 PV-INs exhibit type 2 excitability, like striatal and neocortical PV-INs. Modelling and mathematical analysis showed that the slowly inactivating potassium current Kv1 contributes to type 2 excitability, enables the multiple firing regimes observed experimentally in PV-INs, and provides a mechanism for robust persistent interruption of firing. Using a fast/slow separation of times scales approach with the Kv1 inactivation variable as a bifurcation parameter shows that the initial inhibitory stimulus stops repetitive firing by moving the membrane potential trajectory onto a co-existing stable fixed point corresponding to a non-spiking quiescent state. As Kv1 inactivation decays, the trajectory follows the branch of stable fixed points until it crosses a subcritical Hopf bifurcation then spirals out into repetitive firing. In a model describing entorhinal cortical PV-INs without Kv1, interruption of firing could be achieved by taking advantage of the bistability inherent in type 2 excitability based on a subcritical Hopf bifurcation, but the interruption was not robust to noise. Persistent interruption of firing is therefore broadly applicable to PV-INs in different brain regions but is only made robust to noise in the presence of a slow variable.

Circuit-based intervention corrects excessive dentate gyrus output in the fragile X mouse model

Abnormal cellular and circuit excitability is believed to drive many core phenotypes in fragile X syndrome (FXS). The dentate gyrus is a brain area performing critical computations essential for learning and memory. However, little is known about dentate circuit defects and their mechanisms in FXS. Understanding dentate circuit dysfunction in FXS has been complicated by the presence of two types of excitatory neurons, the granule cells and mossy cells. Here we report that loss of FMRP markedly decreased excitability of dentate mossy cells, a change opposite to all other known excitability defects in excitatory neurons in FXS. This mossy cell hypo-excitability is caused by increased Kv7 function in Fmr1 knockout (KO) mice. By reducing the excitatory drive onto local hilar interneurons, hypo-excitability of mossy cells results in increased excitation/inhibition ratio in granule cells and thus paradoxically leads to excessive dentate output. Circuit-wide inhibition of Kv7 channels in Fmr1 KO mice increases inhibitory drive onto granule cells and normalizes the dentate output in response to physiologically relevant theta-gamma coupling stimulation. Our study suggests that circuit-based interventions may provide a promising strategy in this disorder to bypass irreconcilable excitability defects in different cell types and restore their pathophysiological consequences at the circuit level.

Neuronal Spike Shapes (NSS): A straightforward approach to investigate heterogeneity in neuronal excitability states

The mammalian brain exhibits a remarkable diversity of neurons, contributing to its intricate architecture and functional complexity. The analysis of multimodal single-cell datasets enables the investigation of cell types and states heterogeneity. In this study, we introduce the Neuronal Spike Shapes (NSS), a straightforward approach for the exploration of excitability states of neurons based on their Action Potential (AP) waveforms. The NSS method describes the AP waveform based on a triangular representation complemented by a set of derived electrophysiological (EP) features. To support this hypothesis, we validate the proposed approach on two datasets of murine cortical neurons, focusing it on GABAergic neurons. The validation process involves a combination of NSS-based clustering analysis, features exploration, Differential Expression (DE), and Gene Ontology (GO) enrichment analysis. Results show that the NSS-based analysis captures neuronal excitability states that possess biological relevance independently of cell subtype. In particular, Neuronal Spike Shapes (NSS) captures, among others, a well-characterized fast-spiking excitability state, supported by both electrophysiological and transcriptomic validation. Gene Ontology Enrichment Analysis reveals voltage-gated potassium (K+) channels as specific markers of the identified NSS partitions. This finding strongly corroborates the biological relevance of NSS partitions as excitability states, as the expression of voltage-gated K+ channels regulates the hyperpolarization phase of the AP, being directly implicated in the regulation of neuronal excitability.

Circuit-based intervention corrects excessive dentate gyrus output in the Fragile X mouse model

Abnormal cellular and circuit excitability is believed to drive many core phenotypes in fragile X syndrome (FXS). The dentate gyrus is a brain area performing critical computations essential for learning and memory. However, little is known about dentate circuit defects and their mechanisms in FXS. Understanding dentate circuit dysfunction in FXS has been complicated by the presence of two types of excitatory neurons, the granule cells and mossy cells. Here we report that loss of FMRP markedly decreased excitability of dentate mossy cells, a change opposite to all other known excitability defects in excitatory neurons in FXS. This mossy cell hypo-excitability is caused by increased Kv7 function in Fmr1 KO mice. By reducing the excitatory drive onto local hilar interneurons, hypo-excitability of mossy cells results in increased excitation/inhibition ratio in granule cells and thus paradoxically leads to excessive dentate output. Circuit-wide inhibition of Kv7 channels in Fmr1 KO mice increases inhibitory drive onto granule cells and normalizes the dentate output in response to physiologically relevant theta-gamma coupling stimulation. Our study suggests that circuit-based interventions may provide a promising strategy in this disorder to bypass irreconcilable excitability defects in different cell types and restore their pathophysiological consequences at the circuit level.

Brief synaptic inhibition persistently interrupts firing of fast-spiking interneurons

Neurons perform input-output operations that integrate synaptic inputs with intrinsic electrical properties; these operations are generally constrained by the brevity of synaptic events. Here, we report that sustained firing of CA1 hippocampal fast-spiking parvalbumin-expressing interneurons (PV-INs) can be persistently interrupted for several hundred milliseconds following brief GABAAR-mediated inhibition in vitro and in vivo. A single presynaptic neuron could interrupt PV-IN firing, occasionally with a single action potential (AP), and reliably with AP bursts. Experiments and computational modeling reveal that the persistent interruption of firing maintains neurons in a depolarized, quiescent state through a cell-autonomous mechanism. Interrupted PV-INs are strikingly responsive to Schaffer collateral inputs. The persistent interruption of firing provides a disinhibitory circuit mechanism favoring spike generation in CA1 pyramidal cells. Overall, our results demonstrate that neuronal silencing can far outlast brief synaptic inhibition owing to the well-tuned interplay between neurotransmitter release and postsynaptic membrane dynamics, a phenomenon impacting microcircuit function.

Interacting rhythms enhance sensitivity of target detection in a fronto-parietal computational model of visual attention

Even during sustained attention, enhanced processing of attended stimuli waxes and wanes rhythmically, with periods of enhanced and relatively diminished visual processing (and subsequent target detection) alternating at 4 or 8 Hz in a sustained visual attention task. These alternating attentional states occur alongside alternating dynamical states, in which lateral intraparietal cortex (LIP), the frontal eye field (FEF), and the mediodorsal pulvinar (mdPul) exhibit different activity and functional connectivity at α, β, and γ frequencies-rhythms associated with visual processing, working memory, and motor suppression. To assess whether and how these multiple interacting rhythms contribute to periodicity in attention, we propose a detailed computational model of FEF and LIP. When driven by θ-rhythmic inputs simulating experimentally-observed mdPul activity, this model reproduced the rhythmic dynamics and behavioral consequences of observed attentional states, revealing that the frequencies and mechanisms of the observed rhythms allow for peak sensitivity in visual target detection while maintaining functional flexibility.

Altered excitatory and decreased inhibitory transmission in the prefrontal cortex of male mice with early developmental disruption to the ventral hippocampus

Ventral hippocampal (vHPC)-prefrontal cortical (PFC) pathway dysfunction is a core neuroimaging feature of schizophrenia. However, mechanisms underlying impaired connectivity within this pathway remain poorly understood. The vHPC has direct projections to the PFC that help shape its maturation. Here, we wanted to investigate the effects of early developmental vHPC perturbations on long-term functional PFC organization. Using whole-cell recordings to assess PFC cellular activity in transgenic male mouse lines, we show early developmental disconnection of vHPC inputs, by excitotoxic lesion or cell-specific ablations, impairs pyramidal cell firing output and produces a persistent increase in excitatory and decrease in inhibitory synaptic inputs onto pyramidal cells. We show this effect is specific to excitatory vHPC projection cell ablation. We further identify PV-interneurons as a source of deficit in inhibitory transmission. We find PV-interneurons are reduced in density, show a reduced ability to sustain high-frequency firing, and show deficits in excitatory inputs that emerge over time. We additionally show differences in vulnerabilities to early developmental vHPC disconnection, wherein PFC PV-interneurons but not pyramidal cells show deficits in NMDA receptor-mediated current. Our results highlight mechanisms by which the PFC adapts to early developmental vHPC perturbations, providing insights into schizophrenia circuit pathology.

Interneuronal network model of theta-nested fast oscillations predicts differential effects of heterogeneity, gap junctions and short term depression for hyperpolarizing versus shunting inhibition

Theta and gamma oscillations in the hippocampus have been hypothesized to play a role in the encoding and retrieval of memories. Recently, it was shown that an intrinsic fast gamma mechanism in medial entorhinal cortex can be recruited by optogenetic stimulation at theta frequencies, which can persist with fast excitatory synaptic transmission blocked, suggesting a contribution of interneuronal network gamma (ING). We calibrated the passive and active properties of a 100-neuron model network to capture the range of passive properties and frequency/current relationships of experimentally recorded PV+ neurons in the medial entorhinal cortex (mEC). The strength and probabilities of chemical and electrical synapses were also calibrated using paired recordings, as were the kinetics and short-term depression (STD) of the chemical synapses. Gap junctions that contribute a noticeable fraction of the input resistance were required for synchrony with hyperpolarizing inhibition; these networks exhibited theta-nested high frequency oscillations similar to the putative ING observed experimentally in the optogenetically-driven PV-ChR2 mice. With STD included in the model, the network desynchronized at frequencies above ~200 Hz, so for sufficiently strong drive, fast oscillations were only observed before the peak of the theta. Because hyperpolarizing synapses provide a synchronizing drive that contributes to robustness in the presence of heterogeneity, synchronization decreases as the hyperpolarizing inhibition becomes weaker. In contrast, networks with shunting inhibition required non-physiological levels of gap junctions to synchronize using conduction delays within the measured range.

Cell-Type Specific Inhibition Controls the High-Frequency Oscillations in the Medial Entorhinal Cortex

The medial entorhinal cortex (mEC) plays a critical role for spatial navigation and memory. While many studies have investigated the principal neurons within the entorhinal cortex, much less is known about the inhibitory circuitries within this structure. Here, we describe for the first time in the mEC a subset of parvalbumin-positive (PV+) interneurons (INs)-stuttering cells (STUT)-with morphological, intrinsic electrophysiological, and synaptic properties distinct from fast-spiking PV+ INs. In contrast to the fast-spiking PV+ INs, the axon of the STUT INs also terminated in layer 3 and showed subthreshold membrane oscillations at gamma frequencies. Whereas the synaptic output of the STUT INs was only weakly reduced by a μ-opioid agonist, their inhibitory inputs were strongly suppressed. Given these properties, STUT are ideally suited to entrain gamma activity in the pyramidal cell population of the mEC. We propose that activation of the μ-opioid receptors decreases the GABA release from the PV+ INs onto the STUT, resulting in disinhibition of the STUT cell population and the consequent increase in network gamma power. We therefore suggest that the opioid system plays a critical role, mediated by STUT INs, in the neural signaling and oscillatory network activity within the mEC.

Juvenile social isolation immediately affects the synaptic activity and firing property of fast-spiking parvalbumin-expressing interneuron subtype in mouse medial prefrontal cortex

A lack of juvenile social experience causes various behavioral impairments and brain dysfunction, especially in the medial prefrontal cortex (mPFC). Our previous studies revealed that juvenile social isolation for 2 weeks immediately after weaning affects the synaptic inputs and intrinsic excitability of fast-spiking parvalbumin-expressing (FSPV) interneurons as well as a specific type of layer 5 (L5) pyramidal cells, which we termed prominent h-current (PH) cells, in the mPFC. However, since these changes were observed at the adult age of postnatal day 65 (P65), the primary cause of these changes to neurons immediately after juvenile social isolation (postnatal day 35) remains unknown. Here, we investigated the immediate effects of juvenile social isolation on the excitability and synaptic inputs of PH pyramidal cells and FSPV interneurons at P35 using whole-cell patch-clamp recording. We observed that excitatory inputs to FSPV interneurons increased immediately after juvenile social isolation. We also found that juvenile social isolation increases the firing reactivity of a subtype of FSPV interneurons, whereas only a fractional effect was detected in PH pyramidal cells. These findings suggest that juvenile social isolation primarily disturbs the developmental rebuilding of circuits involving FSPV interneurons and eventually affects the circuits involving PH pyramidal cells in adulthood.

Huygens synchronization of medial septal pacemaker neurons generates hippocampal theta oscillation

Episodic learning and memory retrieval are dependent on hippocampal theta oscillation, thought to rely on the GABAergic network of the medial septum (MS). To test how this network achieves theta synchrony, we recorded MS neurons and hippocampal local field potential simultaneously in anesthetized and awake mice and rats. We show that MS pacemakers synchronize their individual rhythmicity frequencies, akin to coupled pendulum clocks as observed by Huygens. We optogenetically identified them as parvalbumin-expressing GABAergic neurons, while MS glutamatergic neurons provide tonic excitation sufficient to induce theta. In accordance, waxing and waning tonic excitation is sufficient to toggle between theta and non-theta states in a network model of single-compartment inhibitory pacemaker neurons. These results provide experimental and theoretical support to a frequency-synchronization mechanism for pacing hippocampal theta, which may serve as an inspirational prototype for synchronization processes in the central nervous system from Nematoda to Arthropoda to Chordate and Vertebrate phyla.

Deep brain stimulation in the subthalamic nucleus for Parkinson's disease can restore dynamics of striatal networks

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is highly effective in alleviating movement disability in patients with Parkinson’s disease (PD). However, its therapeutic mechanism of action is unknown. The healthy striatum exhibits rich dynamics resulting from an interaction of beta, gamma, and theta oscillations. These rhythms are essential to selection and execution of motor programs, and their loss or exaggeration due to dopamine (DA) depletion in PD is a major source of behavioral deficits. Restoring the natural rhythms may then be instrumental in the therapeutic action of DBS. We develop a biophysical networked model of a BG pathway to study how abnormal beta oscillations can emerge throughout the BG in PD and how DBS can restore normal beta, gamma, and theta striatal rhythms. Our model incorporates STN projections to the striatum, long known but understudied, found to preferentially target fast-spiking interneurons (FSI). We find that DBS in STN can normalize striatal medium spiny neuron activity by recruiting FSI dynamics and restoring the inhibitory potency of FSIs observed in normal conditions. We also find that DBS allows the reexpression of gamma and theta rhythms, thought to be dependent on high DA levels and thus lost in PD, through cortical noise control. Our study highlights that DBS effects can go beyond regularizing BG output dynamics to restoring normal internal BG dynamics and the ability to regulate them. It also suggests how gamma and theta oscillations can be leveraged to supplement DBS treatment and enhance its effectiveness.

Imbalanced Subthreshold Currents Following Sepsis and Chemotherapy: A Shared Mechanism Offering a New Therapeutic Target?

Both sepsis and treatment of cancer with chemotherapy are known to cause neurologic dysfunction. The primary defects seen in both groups of patients are neuropathy and encephalopathy; the underlying mechanisms are poorly understood. Analysis of preclinical models of these disparate conditions reveal similar defects in ion channel function contributing to peripheral neuropathy. The defects in ion channel function extend to the central nervous system where lower motoneurons are affected. In motoneurons the defect involves ion channels responsible for subthreshold currents that convert steady depolarization into repetitive firing. The inability to correctly translate depolarization into steady, repetitive firing has profound effects on motor function, and could be an important contributor to weakness and fatigue experienced by both groups of patients. The possibility that disruption of function, either instead of, or in addition to neurodegeneration, may underlie weakness and fatigue leads to a novel approach to therapy. Activation of serotonin (5HT) receptors in a rat model of sepsis restores the normal balance of subthreshold currents and normal motoneuron firing. If an imbalance of subthreshold currents also occurs in other central nervous system neurons, it could contribute to encephalopathy. We hypothesize that pharmacologically restoring the proper balance of subthreshold currents might provide effective therapy for both neuropathy and encephalopathy in patients recovering from sepsis or treatment with chemotherapy.

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.

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

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

Initiation of migraine-related cortical spreading depolarization by hyperactivity of GABAergic neurons and NaV1.1 channels

Spreading depolarizations (SDs) are involved in migraine, epilepsy, stroke, traumatic brain injury, and subarachnoid hemorrhage. However, the cellular origin and specific differential mechanisms are not clear. Increased glutamatergic activity is thought to be the key factor for generating cortical spreading depression (CSD), a pathological mechanism of migraine. Here, we show that acute pharmacological activation of NaV1.1 (the main Na+ channel of interneurons) or optogenetic-induced hyperactivity of GABAergic interneurons is sufficient to ignite CSD in the neocortex by spiking-generated extracellular K+ build-up. Neither GABAergic nor glutamatergic synaptic transmission were required for CSD initiation. CSD was not generated in other brain areas, suggesting that this is a neocortex-specific mechanism of CSD initiation. Gain-of-function mutations of NaV1.1 (SCN1A) cause familial hemiplegic migraine type-3 (FHM3), a subtype of migraine with aura, of which CSD is the neurophysiological correlate. Our results provide the mechanism linking NaV1.1 gain of function to CSD generation in FHM3. Thus, we reveal the key role of hyperactivity of GABAergic interneurons in a mechanism of CSD initiation, which is relevant as a pathological mechanism of Nav1.1 FHM3 mutations, and possibly also for other types of migraine and diseases in which SDs are involved.

Conventional measures of intrinsic excitability are poor estimators of neuronal activity under realistic synaptic inputs

Activity-dependent regulation of intrinsic excitability has been shown to greatly contribute to the overall plasticity of neuronal circuits. Such neuroadaptations are commonly investigated in patch clamp experiments using current step stimulation and the resulting input-output functions are analyzed to quantify alterations in intrinsic excitability. However, it is rarely addressed, how such changes translate to the function of neurons when they operate under natural synaptic inputs. Still, it is reasonable to expect that a strong correlation and near proportional relationship exist between static firing responses and those evoked by synaptic drive. We challenge this view by performing a high-yield electrophysiological analysis of cultured mouse hippocampal neurons using both standard protocols and simulated synaptic inputs via dynamic clamp. We find that under these conditions the neurons exhibit vastly different firing responses with surprisingly weak correlation between static and dynamic firing intensities. These contrasting responses are regulated by two intrinsic K-currents mediated by Kv1 and Kir channels, respectively. Pharmacological manipulation of the K-currents produces differential regulation of the firing output of neurons. Static firing responses are greatly increased in stuttering type neurons under blocking their Kv1 channels, while the synaptic responses of the same neurons are less affected. Pharmacological blocking of Kir-channels in delayed firing type neurons, on the other hand, exhibit the opposite effects. Our subsequent computational model simulations confirm the findings in the electrophysiological experiments and also show that adaptive changes in the kinetic properties of such currents can even produce paradoxical regulation of the firing output.

Neural excitability increases with axonal resistance between soma and axon initial segment

The position of the axon initial segment (AIS) is thought to play a critical role in neuronal excitability. Previous experimental studies have found that a distal shift in AIS position correlates with a reduction in excitability. Yet theoretical work has suggested the opposite, because of increased electrical isolation. A distal shift in AIS position corresponds to an elevation of axial resistance R _a We therefore examined how changes in R _a at the axon hillock impact the voltage threshold (V_th) of the somatic action potential in L5 pyramidal neurons. Increasing R _a by mechanically pinching the axon between the soma and the AIS was found to lower V_th by ∼6 mV. Conversely, decreasing R _a by substituting internal ions with higher mobility elevated V_th All R _a -dependent changes in V_th could be reproduced in a Hodgkin-Huxley compartmental model. We conclude that in L5 pyramidal neurons, excitability increases with axial resistance and therefore with a distal shift of the AIS.

Modeling NaV1.1/SCN1A sodium channel mutations in a microcircuit with realistic ion concentration dynamics suggests differential GABAergic mechanisms leading to hyperexcitability in epilepsy and hemiplegic migraine

Loss of function mutations of SCN1A, the gene coding for the voltage-gated sodium channel Na_V1.1, cause different types of epilepsy, whereas gain of function mutations cause sporadic and familial hemiplegic migraine type 3 (FHM-3). However, it is not clear yet how these opposite effects can induce paroxysmal pathological activities involving neuronal networks’ hyperexcitability that are specific of epilepsy (seizures) or migraine (cortical spreading depolarization, CSD). To better understand differential mechanisms leading to the initiation of these pathological activities, we used a two-neuron conductance-based model of interconnected GABAergic and pyramidal glutamatergic neurons, in which we incorporated ionic concentration dynamics in both neurons. We modeled FHM-3 mutations by increasing the persistent sodium current in the interneuron and epileptogenic mutations by decreasing the sodium conductance in the interneuron. Therefore, we studied both FHM-3 and epileptogenic mutations within the same framework, modifying only two parameters. In our model, the key effect of gain of function FHM-3 mutations is ion fluxes modification at each action potential (in particular the larger activation of voltage-gated potassium channels induced by the Na_V1.1 gain of function), and the resulting CSD-triggering extracellular potassium accumulation, which is not caused only by modifications of firing frequency. Loss of function epileptogenic mutations, on the other hand, increase GABAergic neurons’ susceptibility to depolarization block, without major modifications of firing frequency before it. Our modeling results connect qualitatively to experimental data: potassium accumulation in the case of FHM-3 mutations and facilitated depolarization block of the GABAergic neuron in the case of epileptogenic mutations. Both these effects can lead to pyramidal neuron hyperexcitability, inducing in the migraine condition depolarization block of both the GABAergic and the pyramidal neuron. Overall, our findings suggest different mechanisms of network hyperexcitability for migraine and epileptogenic Na_V1.1 mutations, implying that the modifications of firing frequency may not be the only relevant pathological mechanism.

The voltage-gated sodium channel Na_V1.1 is a major target of human mutations implicated in different pathologies. In particular, mutations identified in certain types of epilepsy cause loss of function of the channel, whereas mutations identified in certain types of migraine (in which spreading depolarizations of the cortical circuits of the brain are involved) cause instead gain of function. Here, we study dysfunctions induced by these differential effects in a two-neuron (GABAergic and pyramidal) conductance-based model with dynamic ion concentrations. We obtain results that can be related to experimental findings in both situations. Namely, extracellular potassium accumulation induced by the activity of the GABAergic neuron in the case of CSD, and higher propensity of the GABAergic neuron to depolarization block in the epileptogenic scenario, without significant modifications of its firing frequency prior to it. Both scenarios can induce hyperexcitability of the pyramidal neuron, leading in the migraine condition to depolarization block of both the GABAergic and the pyramidal neuron. Our results are successfully confronted to experimental data and suggest that modification of firing frequency is not the only key mechanism in these pathologies of neuronal excitability.

Altered corticostriatal synchronization associated with compulsive-like behavior in APP/PS1 mice

Mild behavioral impairment (MBI), which can include compulsive behavior, is an early sign of Alzheimer's disease (AD), but its underlying neural mechanisms remain unclear. Here, we show that 3-5-month-old APP/PS1 mice display obsessive-compulsive disorder (OCD)-like behavior. The number of parvalbumin-positive (PV) interneurons and level of high gamma (γ_high) oscillation are significantly decreased in the striatum of AD mice. This is accompanied by enhanced β-γ_high coupling and firing rates of putative striatal projection neurons (SPNs), indicating decorrelation between PV interneurons and SPNs. Local field potentials (LFPs) simultaneously recorded in prefrontal cortex (PFC) and striatum (Str) demonstrate a decrease in γ_high-band coherent activity and spike-field coherence in corticostriatal circuits of APP/PS1 mice. Furthermore, levels of GABA_B receptor (GABA_BR), but not GABA_A receptor (GABA_AR), and glutamatergic receptors, were markedly reduced, in line with presymptomatic AD-related behavioral changes. These findings suggest that MBI occurs as early as 3-5 months in APP/PS1 mice and that altered corticostriatal synchronization may play a role in mediating the behavioral phenotypes observed.

Kcns3 deficiency disrupts Parvalbumin neuron physiology in mouse prefrontal cortex: Implications for the pathophysiology of schizophrenia

The unique fast spiking (FS) phenotype of cortical parvalbumin-positive (PV) neurons depends on the expression of multiple subtypes of voltage-gated potassium channels (Kv). PV neurons selectively express Kcns3, the gene encoding Kv9.3 subunits, suggesting that Kcns3 expression is critical for the FS phenotype. KCNS3 expression is lower in PV neurons in the neocortex of subjects with schizophrenia, but the effects of this alteration are unclear, because Kv9.3 subunit function is poorly understood. Therefore, to assess the role of Kv9.3 subunits in PV neuron function, we combined gene expression analyses, computational modeling, and electrophysiology in acute slices from the cortex of Kcns3-deficient mice. Kcns3 mRNA levels were ~ 50% lower in cortical PV neurons from Kcns3-deficient relative to wildtype mice. While silent per se, Kv9.3 subunits are believed to amplify the Kv2.1 current in Kv2.1-Kv9.3 channel complexes. Hence, to assess the consequences of reducing Kv9.3 levels, we simulated the effects of decreasing the Kv2.1-mediated current in a computational model. The FS cell model with reduced Kv2.1 produced spike trains with irregular inter-spike intervals, or stuttering, and greater Na+ channel inactivation. As in the computational model, PV basket cells (PVBCs) from Kcns3-deficient mice displayed spike trains with strong stuttering, which depressed PVBC firing. Moreover, Kcns3 deficiency impaired the recruitment of PVBC firing at gamma frequency by stimuli mimicking synaptic input observed during cortical UP states. Our data indicate that Kv9.3 subunits are critical for PVBC physiology and suggest that KCNS3 deficiency in schizophrenia could impair PV neuron firing, possibly contributing to deficits in cortical gamma oscillations in the illness.

Hyperexcitable Neurons Enable Precise and Persistent Information Encoding in the Superficial Retrosplenial Cortex

The retrosplenial cortex (RSC) is essential for memory and navigation, but the neural codes underlying these functions remain largely unknown. Here, we show that the most prominent cell type in layers 2/3 (L2/3) of the mouse granular RSC is a hyperexcitable, small pyramidal cell. These cells have a low rheobase (LR), high input resistance, lack of spike frequency adaptation, and spike widths intermediate to those of neighboring fast-spiking (FS) inhibitory neurons and regular-spiking (RS) excitatory neurons. LR cells are excitatory but rarely synapse onto neighboring neurons. Instead, L2/3 is a feedforward, not feedback, inhibition-dominated network with dense connectivity between FS cells and from FS to LR neurons. Biophysical models of LR but not RS cells precisely and continuously encode sustained input from afferent postsubicular head-direction cells. Thus, the distinct intrinsic properties of LR neurons can support both the precision and persistence necessary to encode information over multiple timescales in the RSC.

40-Hz Auditory Steady-State Response (ASSR) as a Biomarker of Genetic Defects in the SHANK3 Gene: A Case Report of 15-Year-Old Girl with a Rare Partial SHANK3 Duplication

SHANK3 encodes a scaffold protein involved in postsynaptic receptor density in glutamatergic synapses, including those in the parvalbumin (PV)+ inhibitory neurons—the key players in the generation of sensory gamma oscillations, such as 40-Hz auditory steady-state response (ASSR). However, 40-Hz ASSR was not studied in relation to SHANK3 functioning. Here, we present a 15-year-old girl (SH01) with previously unreported duplication of the first seven exons of the SHANK3 gene (22q13.33). SH01’s electroencephalogram (EEG) during 40-Hz click trains of 500 ms duration binaurally presented with inter-trial intervals of 500–800 ms were compared with those from typically developing children (n = 32). SH01 was diagnosed with mild mental retardation and learning disabilities (F70.88), dysgraphia, dyslexia, and smaller vocabulary than typically developing (TD) peers. Her clinical phenotype resembled the phenotype of previously described patients with 22q13.33 microduplications (≈30 reported so far). SH01 had mild autistic symptoms but below the threshold for ASD diagnosis and microcephaly. No seizures or MRI abnormalities were reported. While SH01 had relatively preserved auditory event-related potential (ERP) with slightly attenuated P1, her 40-Hz ASSR was totally absent significantly deviating from TD’s ASSR. The absence of 40-Hz ASSR in patients with microduplication, which affected the SHANK3 gene, indicates deficient temporal resolution of the auditory system, which might underlie language problems and represent a neurophysiological biomarker of SHANK3 abnormalities.

Nucleus accumbens fast-spiking interneurons in motivational and addictive behaviors

The development of drug addiction is associated with functional adaptations within the reward circuitry, within which the nucleus accumbens (NAc) is anatomically positioned as an interface between motivational salience and behavioral output. The functional output of NAc is profoundly altered after exposure to drugs of abuse, and some of the functional changes continue to evolve during drug abstinence, contributing to numerous emotional and motivational alterations related drug taking, seeking, and relapse. As in most brain regions, the functional output of NAc is critically dependent on the dynamic interaction between excitation and inhibition. One of the most prominent sources of inhibition within the NAc arises from fast-spiking interneurons (FSIs). Each NAc FSI innervates hundreds of principal neurons, and orchestrates population activity through its powerful and sustained feedforward inhibition. While the role of NAc FSIs in the context of drug addiction remains poorly understood, emerging evidence suggests that FSIs and FSI-mediated local circuits are key targets for drugs of abuse to tilt the functional output of NAc toward a motivational state favoring drug seeking and relapse. In this review, we discuss recent findings and our conceptualization about NAc FSI-mediated regulation of motivated and cocaine-induced behaviors. We hope that the conceptual framework proposed in this review may provide a useful guidance for ongoing and future studies to determine how FSIs influence the function of NAc and related reward circuits, ultimately leading to addictive behaviors.

Characterization of the Synergistic Inhibition of IK(erg) and IK(DR) by Ribociclib, a Cyclin-Dependent Kinase 4/6 Inhibitor

Ribociclib (RIB, LE011, Kisqali^®), an orally administered inhibitor of cyclin-dependent kinase-4/6 (CDK-4/6) complex, is clinically effective for the treatment of several malignancies, including advanced breast cancer. However, information regarding the effects of RIB on membrane ion currents is limited. In this study, the addition of RIB to pituitary tumor (GH₃) cells decreased the peak amplitude of erg-mediated K⁺ current (I_K(erg)), which was accompanied by a slowed deactivation rate of the current. The IC₅₀ value for RIB-perturbed inhibition of deactivating I_K(erg) in these cells was 2.7 μM. In continued presence of μM RIB, neither the subsequent addition of 17β-estradiol (30 μM), phorbol 12-myristate 13-acetate (10 μM), or transforming growth factor-β (1 μM) counteracted the inhibition of deactivating I_K(erg). Its presence affected the decrease in the degree of voltage-dependent hysteresis for I_K(erg) elicitation by long-duration triangular ramp voltage commands. The presence of RIB differentially inhibited the peak or sustained component of delayed rectifier K⁺ current (I_K(DR)) with an effective IC₅₀ of 28.7 or 11.4 μM, respectively, while it concentration-dependently decreased the amplitude of M-type K⁺ current with IC₅₀ of 13.3 μM. Upon 10-s long membrane depolarization, RIB elicited a decrease in the I_K(DR) amplitude, which was concomitant with an accelerated inactivation time course. However, the inability of RIB (10 μM) to modify the magnitude of the hyperpolarization-activated cation current was disclosed. The mean current–voltage relationship of I_K(erg) present in HL-1 atrial cardiomyocytes was inhibited in the presence of RIB (10 μM). Collectively, the hyperpolarization-activated cation current was observed. RIB-mediated perturbations in ionic currents presented herein are upstream of its suppressive action on cytosolic CDK-4/6 activities and partly participates in its modulatory effects on the functional activities of pituitary tumor cells (e.g., GH₃ cells) or cardiac myocytes (e.g., HL-1 cells).

Chronic defects in intraspinal mechanisms of spike encoding by spinal motoneurons following chemotherapy

Chemotherapy-induced sensorimotor disabilities, including gait and balance disorders, as well as physical fatigue often persist for months and sometimes years into disease free survival from cancer. While associated with impaired sensory function, chronic sensorimotor disorders might also depend on chemotherapy-induced defects in other neuron types. In this report, we extend consideration to motoneurons, which, if chronically impaired, would necessarily degrade movement behavior. The present study was undertaken to determine whether motoneurons qualify as candidate contributors to chronic sensorimotor disability independently from sensory impairment. We tested this possibility in vivo from rats 5 weeks following human-scaled treatment with one of the platinum-based compounds, oxaliplatin, widely used in chemotherapy for a variety of cancers. Action potential firing of spinal motoneurons responding to different fixed levels of electrode-current injection was measured in order to assess the neurons' intrinsic capacity for stimulus encoding. The encoding of stimulus duration and intensity corroborated in untreated control rats was severely degraded in oxaliplatin treated rats, in which motoneurons invariably exhibited erratic firing that was unsustained, unpredictable from one stimulus trial to the next, and unresponsive to changes in current strength. Direct measurements of interspike oscillations in membrane voltage combined with computer modeling pointed to aberrations in subthreshold conductances as a plausible contributor to impaired firing behavior. These findings authenticate impaired spike encoding as a candidate contributor to, in the case of motoneurons, deficits in mobility and fatigue. Aberrant firing also becomes a deficit worthy of testing in other CNS neurons as a potential contributor to perceptual and cognitive disorders induced by chemotherapy in patients.

Auditory Long-Range Parvalbumin Cortico-Striatal Neurons

Previous studies have shown that cortico-striatal pathways link auditory signals to action-selection and reward-learning behavior through excitatory projections. Only recently it has been demonstrated that long-range GABAergic cortico-striatal somatostatin-expressing neurons in the auditory cortex project to the dorsal striatum, and functionally inhibit the main projecting neuronal population, the spiny projecting neuron. Here we tested the hypothesis that parvalbumin-expressing neurons of the auditory cortex can also send long-range projections to the auditory striatum. To address this fundamental question, we took advantage of viral and non-viral anatomical tracing approaches to identify cortico-striatal parvalbumin neurons (CS-Parv inhibitory projections → auditory striatum). Here, we describe their anatomical distribution in the auditory cortex and determine the anatomical and electrophysiological properties of layer 5 CS-Parv neurons. We also analyzed their characteristic voltage-dependent membrane potential gamma oscillation, showing that intrinsic membrane mechanisms generate them. The inherent membrane mechanisms can also trigger intermittent and irregular bursts (stuttering) of the action potential in response to steps of depolarizing current pulses.

Shunting Inhibition Improves Synchronization in Heterogeneous Inhibitory Interneuronal Networks with Type 1 Excitability Whereas Hyperpolarizing Inhibition Is Better for Type 2 Excitability

All-to-all homogeneous networks of inhibitory neurons synchronize completely under the right conditions; however, many modeling studies have shown that biological levels of heterogeneity disrupt synchrony. Our fundamental scientific question is “how can neurons maintain partial synchrony in the presence of heterogeneity and noise?” A particular subset of strongly interconnected interneurons, the PV+ fast-spiking (FS) basket neurons, are strongly implicated in γ oscillations and in phase locking of nested γ oscillations to theta. Their excitability type apparently varies between brain regions: in CA1 and the dentate gyrus they have type 1 excitability, meaning that they can fire arbitrarily slowly, whereas in the striatum and cortex they have type 2 excitability, meaning that there is a frequency thresh old below which they cannot sustain repetitive firing. We constrained the models to study the effect of excitability type (more precisely bifurcation type) in isolation from all other factors. We use sparsely connected, heterogeneous, noisy networks with synaptic delays to show that synchronization properties, namely the resistance to suppression and the strength of theta phase to γ amplitude coupling, are strongly dependent on the pairing of excitability type with the type of inhibition. Shunting inhibition performs better for type 1 and hyperpolarizing inhibition for type 2. γ Oscillations and their nesting within theta oscillations are thought to subserve cognitive functions like memory encoding and recall; therefore, it is important to understand the contribution of intrinsic properties to these rhythms.

A biophysical model of striatal microcircuits suggests gamma and beta oscillations interleaved at delta/theta frequencies mediate periodicity in motor control

Striatal oscillatory activity is associated with movement, reward, and decision-making, and observed in several interacting frequency bands. Local field potential recordings in rodent striatum show dopamine- and reward-dependent transitions between two states: a "spontaneous" state involving β (∼15-30 Hz) and low γ (∼40-60 Hz), and a state involving θ (∼4-8 Hz) and high γ (∼60-100 Hz) in response to dopaminergic agonism and reward. The mechanisms underlying these rhythmic dynamics, their interactions, and their functional consequences are not well understood. In this paper, we propose a biophysical model of striatal microcircuits that comprehensively describes the generation and interaction of these rhythms, as well as their modulation by dopamine. Building on previous modeling and experimental work suggesting that striatal projection neurons (SPNs) are capable of generating β oscillations, we show that networks of striatal fast-spiking interneurons (FSIs) are capable of generating δ/θ (ie, 2 to 6 Hz) and γ rhythms. Under simulated low dopaminergic tone our model FSI network produces low γ band oscillations, while under high dopaminergic tone the FSI network produces high γ band activity nested within a δ/θ oscillation. SPN networks produce β rhythms in both conditions, but under high dopaminergic tone, this β oscillation is interrupted by δ/θ-periodic bursts of γ-frequency FSI inhibition. Thus, in the high dopamine state, packets of FSI γ and SPN β alternate at a δ/θ timescale. In addition to a mechanistic explanation for previously observed rhythmic interactions and transitions, our model suggests a hypothesis as to how the relationship between dopamine and rhythmicity impacts motor function. We hypothesize that high dopamine-induced periodic FSI γ-rhythmic inhibition enables switching between β-rhythmic SPN cell assemblies representing the currently active motor program, and thus that dopamine facilitates movement in part by allowing for rapid, periodic shifts in motor program execution.

Frequency-dependent entrainment of striatal fast-spiking interneurons

Striatal fast-spiking interneurons (FSIs) fire in variable-length runs of action potentials at 20-200 spikes/s separated by pauses. In vivo, or with fluctuating applied current, both runs and pauses become briefer and more variable. During runs, spikes are entrained specifically to gamma-frequency components of the input fluctuations. We stimulated parvalbumin-expressing striatal FSIs in mouse brain slices with broadband noise currents added to direct current steps and measured spike entrainment across all frequencies. As the constant current level was increased, FSIs produced longer runs and showed sharper frequency tuning, with best entrainment at the stimulus frequency matching their intrarun firing rate. We separated the contributions of previous spikes from that of the fluctuating stimulus, revealing a strong contribution of previous action potentials to gamma-frequency entrainment. In contrast, after subtraction of the effect inherited from the previous spike, the remaining stimulus contribution to spike generation was less sharply tuned, showing a larger contribution of lower frequencies. The frequency specificity of entrainment within a run was reproduced with a phase resetting model based on experimentally measured phase resetting curves of the same FSIs. In the model, broadly tuned phase entrainment for the first spike in a run evolved into sharply tuned gamma entrainment over the next few spikes. The data and modeling results indicate that for FSIs firing in brief runs and pauses firing within runs is entrained by gamma-frequency components of the input, whereas the onset timing of runs may be sensitive to a wider range of stimulus frequency components.NEW & NOTEWORTHY Specific types of neurons entrain their spikes to particular oscillation frequencies in their synaptic input. This entrainment is commonly understood in terms of the subthreshold voltage response, but how this translates to spiking is not clear. We show that in striatal fast-spiking interneurons, entrainment to gamma-frequency input depends on rhythmic spike runs and is explained by the phase resetting curve, whereas run initiation can be triggered by a broad range of input frequencies.

Biophysical Modeling Suggests Optimal Drug Combinations for Improving the Efficacy of GABA Agonists after Traumatic Brain Injuries

Traumatic brain injuries (TBI) lead to dramatic changes in the surviving brain tissue. Altered ion concentrations, coupled with changes in the expression of membrane-spanning proteins, create a post-TBI brain state that can lead to further neuronal loss caused by secondary excitotoxicity. Several GABA receptor agonists have been tested in the search for neuroprotection immediately after an injury, with paradoxical results. These drugs not only fail to offer neuroprotection, but can also slow down functional recovery after TBI. Here, using computational modeling, we provide a biophysical hypothesis to explain these observations. We show that the accumulation of intracellular chloride ions caused by a transient upregulation of Na+-K+-2Cl- (NKCC1) co-transporters as observed following TBI, causes GABA receptor agonists to lead to excitation and depolarization block, rather than the expected hyperpolarization. The likelihood of prolonged, excitotoxic depolarization block is further exacerbated by the extremely high levels of extracellular potassium seen after TBI. Our modeling results predict that the neuroprotective efficacy of GABA receptor agonists can be substantially enhanced when they are combined with NKCC1 co-transporter inhibitors. This suggests a rational, biophysically principled method for identifying drug combinations for neuroprotection after TBI.

Ionic mechanisms underlying tonic and burst firing behavior in subfornical organ neurons: a combined experimental and modeling study

Subfornical organ (SFO) neurons exhibit heterogeneity in current expression and spiking behavior, where the two major spiking phenotypes appear as tonic and burst firing. Insight into the mechanisms behind this heterogeneity is critical for understanding how the SFO, a sensory circumventricular organ, integrates and selectively influences physiological function. To integrate efficient methods for studying this heterogeneity, we built a single-compartment, Hodgkin-Huxley-type model of an SFO neuron that is parameterized by SFO-specific in vitro patch-clamp data. The model accounts for the membrane potential distribution and spike train variability of both tonic and burst firing SFO neurons. Analysis of model dynamics confirms that a persistent Na⁺ and Ca²⁺ currents are required for burst initiation and maintenance and suggests that a slow-activating K⁺ current may be responsible for burst termination in SFO neurons. Additionally, the model suggests that heterogeneity in current expression and subsequent influence on spike afterpotential underlie the behavioral differences between tonic and burst firing SFO neurons. Future use of this model in coordination with single neuron patch-clamp electrophysiology provides a platform for explaining and predicting the response of SFO neurons to various combinations of circulating signals, thus elucidating the mechanisms underlying physiological signal integration within the SFO. NEW & NOTEWORTHY Our understanding of how the subfornical organ (SFO) selectively influences autonomic nervous system function remains incomplete but theoretically results from the electrical responses of SFO neurons to physiologically important signals. We have built a computational model of SFO neurons, derived from and supported by experimental data, which explains how SFO neurons produce different electrical patterns. The model provides an efficient system to theoretically and experimentally explore how changes in the essential features of SFO neurons affect their electrical activity.

Firing rate equations require a spike synchrony mechanism to correctly describe fast oscillations in inhibitory networks

Recurrently coupled networks of inhibitory neurons robustly generate oscillations in the gamma band. Nonetheless, the corresponding Wilson-Cowan type firing rate equation for such an inhibitory population does not generate such oscillations without an explicit time delay. We show that this discrepancy is due to a voltage-dependent spike-synchronization mechanism inherent in networks of spiking neurons which is not captured by standard firing rate equations. Here we investigate an exact low-dimensional description for a network of heterogeneous canonical Class 1 inhibitory neurons which includes the sub-threshold dynamics crucial for generating synchronous states. In the limit of slow synaptic kinetics the spike-synchrony mechanism is suppressed and the standard Wilson-Cowan equations are formally recovered as long as external inputs are also slow. However, even in this limit synchronous spiking can be elicited by inputs which fluctuate on a time-scale of the membrane time-constant of the neurons. Our meanfield equations therefore represent an extension of the standard Wilson-Cowan equations in which spike synchrony is also correctly described.

Population models describing the average activity of large neuronal ensembles are a powerful mathematical tool to investigate the principles underlying cooperative function of large neuronal systems. However, these models do not properly describe the phenomenon of spike synchrony in networks of neurons. In particular, they fail to capture the onset of synchronous oscillations in networks of inhibitory neurons. We show that this limitation is due to a voltage-dependent synchronization mechanism which is naturally present in spiking neuron models but not captured by traditional firing rate equations. Here we investigate a novel set of macroscopic equations which incorporate both firing rate and membrane potential dynamics, and that correctly generate fast inhibition-based synchronous oscillations. In the limit of slow-synaptic processing oscillations are suppressed, and the model reduces to an equation formally equivalent to the Wilson-Cowan model.

Increasing motor neuron excitability to treat weakness in sepsis

OBJECTIVE

Weakness induced by critical illness (intensive care unit acquired weakness) is a major cause of disability in patients and is currently untreatable. We recently identified a defect in repetitive firing of lower motor neurons as a novel contributor to intensive care unit acquired weakness. To develop therapy for intensive care unit acquired weakness, it was necessary to determine the mechanism underlying the defect in repetitive firing.

METHODS

Both computer simulation and in vivo dynamic voltage clamp of spinal motor neurons in septic rats were employed to explore potential mechanisms underlying defective repetitive firing.

RESULTS

Our results suggest alteration in subthreshold voltage-activated currents might be the mechanism underlying defective repetitive firing. It has been shown previously that pharmacologic activation of serotonin receptors on motor neurons increases motor neuron excitability, in part by enhancing subthreshold voltage-activated inward currents. Administration of a U.S. Food and Drug Administration-approved serotonin agonist (lorcaserin) to septic rats greatly improved repetitive firing and motor unit force generation.

INTERPRETATION

Our findings suggest activation of serotonin receptors with lorcaserin may provide the first ever therapy for intensive care unit acquired weakness in patients. Ann Neurol 2017;82:961-971.

Pallidostriatal Projections Promote β Oscillations in a Dopamine-Depleted Biophysical Network Model

In the basal ganglia, focused rhythmicity is an important feature of network activity at certain stages of motor processing. In disease, however, the basal ganglia develop amplified rhythmicity. Here, we demonstrate how the cellular architecture and network dynamics of an inhibitory loop in the basal ganglia yield exaggerated synchrony and locking to β oscillations, specifically in the dopamine-depleted state. A key component of this loop is the pallidostriatal pathway, a well-characterized anatomical projection whose function has long remained obscure. We present a synaptic characterization of this pathway in mice and incorporate these data into a computational model that we use to investigate its influence over striatal activity under simulated healthy and dopamine-depleted conditions. Our model predicts that the pallidostriatal pathway influences striatal output preferentially during periods of synchronized activity within GPe. We show that, under dopamine-depleted conditions, this effect becomes a key component of a positive feedback loop between the GPe and striatum that promotes synchronization and rhythmicity. Our results generate novel predictions about the role of the pallidostriatal pathway in shaping basal ganglia activity in health and disease.

SIGNIFICANCE STATEMENT

This work demonstrates that functional connections from the globus pallidus externa (GPe) to striatum are substantially stronger onto fast-spiking interneurons (FSIs) than onto medium spiny neurons. Our circuit model suggests that when GPe spikes are synchronous, this pallidostriatal pathway causes synchronous FSI activity pauses, which allow a transient window of disinhibition for medium spiny neurons. In simulated dopamine-depletion, this GPe-FSI activity is necessary for the emergence of strong synchronization and the amplification and propagation of β oscillations, which are a hallmark of parkinsonian circuit dysfunction. These results suggest that GPe may play a central role in propagating abnormal circuit activity to striatum, which in turn projects to downstream basal ganglia structures. These findings warrant further exploration of GPe as a target for interventions for Parkinson's disease.

The shaping of intrinsic membrane potential oscillations: positive/negative feedback, ionic resonance/amplification, nonlinearities and time scales

The generation of intrinsic subthreshold (membrane potential) oscillations (STOs) in neuronal models requires the interaction between two processes: a relatively fast positive feedback that favors changes in voltage and a slower negative feedback that opposes these changes. These are provided by the so-called resonant and amplifying gating variables associated to the participating ionic currents. We investigate both the biophysical and dynamic mechanisms of generation of STOs and how their attributes (frequency and amplitude) depend on the model parameters for biophysical (conductance-based) models having qualitatively different types of resonant currents (activating and inactivating) and an amplifying current. Combinations of the same types of ionic currents (same models) in different parameter regimes give rise to different types of nonlinearities in the voltage equation: quasi-linear, parabolic-like and cubic-like. On the other hand, combinations of different types of ionic currents (different models) may give rise to the same type of nonlinearities. We examine how the attributes of the resulting STOs depend on the combined effect of these resonant and amplifying ionic processes, operating at different effective time scales, and the various types of nonlinearities. We find that, while some STO properties and attribute dependencies on the model parameters are determined by the specific combinations of ionic currents (biophysical properties), and are different for models with different such combinations, others are determined by the type of nonlinearities and are common for models with different types of ionic currents. Our results highlight the richness of STO behavior in single cells as the result of the various ways in which resonant and amplifying currents interact and affect the generation and termination of STOs as control parameters change. We make predictions that can be tested experimentally and are expected to contribute to the understanding of how rhythmic activity in neuronal networks emerge from the interplay of the intrinsic properties of the participating neurons and the network connectivity.

Stochastic and deterministic dynamics of intrinsically irregular firing in cortical inhibitory interneurons

Most cortical neurons fire regularly when excited by a constant stimulus. In contrast, irregular-spiking (IS) interneurons are remarkable for the intrinsic variability of their spike timing, which can synchronize amongst IS cells via specific gap junctions. Here, we have studied the biophysical mechanisms of this irregular spiking in mice, and how IS cells fire in the context of synchronous network oscillations. Using patch-clamp recordings, artificial dynamic conductance injection, pharmacological analysis and computational modeling, we show that spike time irregularity is generated by a nonlinear dynamical interaction of voltage-dependent sodium and fast-inactivating potassium channels just below spike threshold, amplifying channel noise. This active irregularity may help IS cells synchronize with each other at gamma range frequencies, while resisting synchronization to lower input frequencies.

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

Neurons send information to other neurons in the brain by generating fast electrical pulses called action potentials (or spikes). When stimulated by input signals of a constant size, neurons generally respond with regular patterns of spiking leading to rhythmical brain activity. However, neurons known as irregular spiking interneurons are unique: the relationship between the input they receive and whether or not they produce a spike appears to be random. The molecular mechanism behind this phenomenon is not clear.

Mendonça et al. set out to investigate whether irregular spiking is truly random, or whether there is some degree of predictability. The experiments used genetically modified mice in which irregular spiking interneurons were specifically labeled with a fluorescent protein, which made them easier to find to record their electrical activity. Sophisticated statistical analyses showed that these neurons are not firing at random. Instead, there is a pattern to the timings of the spikes they produce.

It was previously known that electrical spikes in neurons are generated by sodium ions and potassium ions moving across the membrane that surrounds each cell. Proteins called ion channels provide routes for these ions to pass through the membrane. Mendonça et al. show that compared to other neurons, irregular spiking interneurons have larger numbers of a specific type of potassium ion channel. Mimicking the effect of increasing the number of these potassium ion channels in the interneurons made the firing pattern of these neurons more irregular, while decreasing the number of these channels made the firing patterns more regular and predictable.

A computer model of an irregular spiking interneuron showed that the activity of these potassium ion channels and a type of sodium ion channel plays a key role in producing irregular electrical spiking. Further analysis showed that irregular spiking interneurons can synchronize their activity with fast, but not slow, rhythms in brain activity.

The findings of Mendonça et al. suggest that irregular spiking interneurons can disrupt slow regular electrical activity in the brain. Rhythms in brain activity vary depending on whether we are awake or asleep, and are altered in diseases such as epilepsy and schizophrenia. Now that we have a better understanding of how irregular spiking interneurons work, it should be possible to find out how they coordinate their activity with each other, and what effect they have on animal behavior.

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

Hierarchical spike clustering analysis for investigation of interneuron heterogeneity

Action potentials represent the output of a neuron. Especially interneurons display a variety of discharge patterns ranging from regular action potential firing to prominent spike clustering or stuttering. The mechanisms underlying this heterogeneity remain incompletely understood. We established hierarchical cluster analysis of spike trains as a measure of spike clustering. A clustering index was calculated from action potential trains recorded in the whole-cell patch clamp configuration from hippocampal (CA1, stratum radiatum) and entorhinal (medial entorhinal cortex, layer 2) interneurons in acute slices and simulated data. Prominent, region-dependent, but also variable spike clustering was detected using this measure. Further analysis revealed a strong positive correlation between spike clustering and membrane potentials oscillations but an inverse correlation with neuronal resonance. Furthermore, clustering was more pronounced when the balance between fast-activating K(+) currents, assessed by the spike repolarisation time, and hyperpolarization-activated currents, gauged by the size of the sag potential, was shifted in favour of fast K(+) currents. Simulations of spike clustering confirmed that variable ratios of fast K(+) and hyperpolarization-activated currents could underlie different degrees of spike clustering and could thus be crucial for temporally structuring interneuron spike output.

Optogenetic stimulation reveals distinct modulatory properties of thalamostriatal vs corticostriatal glutamatergic inputs to fast-spiking interneurons

Parvalbumin-containing fast-spiking interneurons (FSIs) exert a powerful feed-forward GABAergic inhibition on striatal medium spiny neurons (MSNs), playing a critical role in timing striatal output. However, how glutamatergic inputs modulate their firing activity is still unexplored. Here, by means of a combined optogenetic and electrophysiological approach, we provide evidence for a differential modulation of cortico- vs thalamo-striatal synaptic inputs to FSIs in transgenic mice carrying light-gated ion channels channelrhodopsin-2 (ChR2) in glutamatergic fibers. Corticostriatal synapses show a postsynaptic facilitation, whereas thalamostriatal synapses present a postsynaptic depression. Moreover, thalamostriatal synapses exhibit more prominent AMPA-mediated currents than corticostriatal synapses, and an increased release probability. Furthermore, during current-evoked firing activity, simultaneous corticostriatal stimulation increases bursting activity. Conversely, thalamostriatal fiber activation shifts the canonical burst-pause activity to a more prolonged, regular firing pattern. However, this change in firing pattern was accompanied by a significant rise in the frequency of membrane potential oscillations. Notably, the responses to thalamic stimulation were fully abolished by blocking metabotropic glutamate 1 (mGlu1) receptor subtype, whereas both acetylcholine and dopamine receptor antagonists were ineffective. Our findings demonstrate that cortical and thalamic glutamatergic input differently modulate FSIs firing activity through specific intrinsic and synaptic properties, exerting a powerful influence on striatal outputs.

Postnatal development attunes olfactory bulb mitral cells to high-frequency signaling

Mitral cells (MCs) are a major class of principal neurons in the vertebrate olfactory bulb, conveying odor-evoked activity from the peripheral sensory neurons to olfactory cortex. Previous work has described the development of MC morphology and connectivity during the first few weeks of postnatal development. However, little is known about the postnatal development of MC intrinsic biophysical properties. To understand stimulus encoding in the developing olfactory bulb, we have therefore examined the development of MC intrinsic biophysical properties in acute slices from postnatal day (P)7-P35 mice. Across development, we observed systematic changes in passive membrane properties and action potential waveforms consistent with a developmental increase in sodium and potassium conductances. We further observed developmental decreases in hyperpolarization-evoked membrane potential sag and firing regularity, extending recent links between MC sag heterogeneity and firing patterns. We then applied a novel combination of statistical analyses to examine how the evolution of these intrinsic biophysical properties specifically influenced the representation of fluctuating stimuli by MCs. We found that immature MCs responded to frozen fluctuating stimuli with lower firing rates, lower spike-time reliability, and lower between-cell spike-time correlations than more mature MCs. Analysis of spike-triggered averages revealed that these changes in spike timing were driven by a developmental shift from broad integration of inputs to more selective detection of coincident inputs. Consistent with this shift, generalized linear model fits to MC firing responses demonstrated an enhanced encoding of high-frequency stimulus features by mature MCs.

Interneurons and oligodendrocyte progenitors form a structured synaptic network in the developing neocortex

NG2 cells, oligodendrocyte progenitors, receive a major synaptic input from interneurons in the developing neocortex. It is presumed that these precursors integrate cortical networks where they act as sensors of neuronal activity. We show that NG2 cells of the developing somatosensory cortex form a transient and structured synaptic network with interneurons that follows its own rules of connectivity. Fast-spiking interneurons, highly connected to NG2 cells, target proximal subcellular domains containing GABA_A receptors with γ2 subunits. Conversely, non-fast-spiking interneurons, poorly connected with these progenitors, target distal sites lacking this subunit. In the network, interneuron-NG2 cell connectivity maps exhibit a local spatial arrangement reflecting innervation only by the nearest interneurons. This microcircuit architecture shows a connectivity peak at PN10, coinciding with a switch to massive oligodendrocyte differentiation. Hence, GABAergic innervation of NG2 cells is temporally and spatially regulated from the subcellular to the network level in coordination with the onset of oligodendrogenesis.

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

Neurons are outnumbered in the brain by cells called glial cells. The brain contains various types of glial cells that perform a range of different jobs, including the supply of nutrients and the removal of dead neurons. The role of glial cells called oligodendrocytes is to produce a material called myelin: this is an electrical insulator that, when wrapped around a neuron, increases the speed at which electrical impulses can travel through the nervous system.

Neurons communicate with one another through specialized junctions called synapses, and at one time it was thought that only neurons could form synapses in the brain. However, this view had to be revised when researchers discovered synapses between neurons and glial cells called NG2 cells, which go on to become oligodendrocytes. These neuron-NG2 cell synapses have a lot in common with neuron–neuron synapses, but much less is known about them.

Orduz, Maldonado et al. have now examined these synapses in unprecedented detail by analyzing individual synapses between a type of neuron called an interneuron and an NG2 cell in mice aged only a few weeks. Interneurons can be divided into two major classes based on how quickly they fire, and Orduz, Maldonado et al. show that both types of interneuron form synapses with NG2 cells. However, these two types of interneuron establish synapses on different parts of the NG2 cell, and these synapses involve different receptor proteins.

Together, the synapses give rise to a local interneuron-NG2 cell network that reaches a peak of activity roughly two weeks after birth, after which the network is disassembled. This period of peak activity is accompanied by a sudden increase in the maturation of NG2 cells into oligodendrocytes. Further experiments are needed to test the possibility that activity in the interneuron-NG2 cell network acts as the trigger for the NG2 cells to turn into oligodendrocytes, which then supply myelin for the developing brain.

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