Impact of adaptation currents on synchronization of coupled exponential integrate-and-fire neurons

The Roles of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slow-wave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking (RS) cells with frequency adaptation and do not exhibit bursts in current-clamp experiments (in vitro). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings (gsyn) and external current inputs (I). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to the bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold (IL) and low-threshold (IT) ion channels play a role in increasing and decreasing the parameter conditions (gsyn and I) in which bistable dynamics occur, respectively. For high values of IL conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance IT increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

The Role of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slowwave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking cells (RS) with frequency adaptation and do not exhibit bursts in current-clamp experiments ( in vitro ). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings ( g syn ) and external current inputs ( I ). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold ( I L ) and low-threshold ( I T ) ion channels play a role in increasing and decreasing the parameter conditions ( g syn and I ) in which bistable dynamics occur, respectively. For high values of I L conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance I T increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

Synchronization of relaxation oscillators with adaptive thresholds and application to automated guided vehicles

The present paper proposes an adaptive control law for inducing in-phase and antiphase synchronization in a pair of relaxation oscillators. We analytically show that the phase dynamics of the oscillators coupled by the control law is equivalent to that of Kuramoto phase oscillators and then extend the results for a pair of oscillators to three or more oscillators. We also provide a systematic procedure for designing the controller parameters for oscillator networks with all-to-all and ring topologies. Our numerical simulations demonstrate that these analytical results can be used to solve a dispatching problem encountered by automated guided vehicles (AGVs) in factories. AGV congestion can be avoided and the peak value of the amount of materials or parts in buffers can be suppressed.

Interspike interval correlations in neuron models with adaptation and correlated noise

The generation of neural action potentials (spikes) is random but nevertheless may result in a rich statistical structure of the spike sequence. In particular, contrary to the popular renewal assumption of theoreticians, the intervals between adjacent spikes are often correlated. Experimentally, different patterns of interspike-interval correlations have been observed and computational studies have identified spike-frequency adaptation and correlated noise as the two main mechanisms that can lead to such correlations. Analytical studies have focused on the single cases of either correlated (colored) noise or adaptation currents in combination with uncorrelated (white) noise. For low-pass filtered noise or adaptation, the serial correlation coefficient can be approximated as a single geometric sequence of the lag between the intervals, providing an explanation for some of the experimentally observed patterns. Here we address the problem of interval correlations for a widely used class of models, multidimensional integrate-and-fire neurons subject to a combination of colored and white noise sources and a spike-triggered adaptation current. Assuming weak noise, we derive a simple formula for the serial correlation coefficient, a sum of two geometric sequences, which accounts for a large class of correlation patterns. The theory is confirmed by means of numerical simulations in a number of special cases including the leaky, quadratic, and generalized integrate-and-fire models with colored noise and spike-frequency adaptation. Furthermore we study the case in which the adaptation current and the colored noise share the same time scale, corresponding to a slow stochastic population of adaptation channels; we demonstrate that our theory can account for a nonmonotonic dependence of the correlation coefficient on the channel's time scale. Another application of the theory is a neuron driven by network-noise-like fluctuations (green noise). We also discuss the range of validity of our weak-noise theory and show that by changing the relative strength of white and colored noise sources, we can change the sign of the correlation coefficient. Finally, we apply our theory to a conductance-based model which demonstrates its broad applicability.

Transition to synchronization in heterogeneous inhibitory neural networks with structured synapses

Inhibitory neurons form an extensive network involved in the development of different rhythms in the cerebral cortex. A transition from an incoherent state, where all inhibitory neurons fire unrelated to each other, to a synchronized or locked state, where all or most neurons define a tight firing pattern, is maybe the most salient process to analyze when considering neuronal rhythms. In this work, we analyzed whether different patterns of effective synaptic connectivity may support a first-order-like transition in this path to synchronization. Such an "explosive" phenomenon may be relevant in neural processes, as normal cognitive processing in different tasks and some neurological disorders manifest an increased power in many neuronal rhythms, supported by an extended concerted spiking activity and an abrupt change to this state. Furthermore, we built an adaptive mechanism that supports the generation of this kind of network, which rapidly creates the underlying structure based on the ongoing firing statistics.

Synchronization of Electrically Coupled Resonate-and-Fire Neurons

Electrical coupling between neurons is broadly present across brain areas and is typically assumed to synchronize network activity. However, intrinsic properties of the coupled cells can complicate this simple picture. Many cell types with electrical coupling show a diversity of post-spike subthreshold fluctuations, often linked to subthreshold resonance, which are transmitted through electrical synapses in addition to action potentials. Using the theory of weakly coupled oscillators, we explore the effect of both subthreshold and spike-mediated coupling on synchrony in small networks of electrically coupled resonate-and-fire neurons, a hybrid neuron model with damped subthreshold oscillations and a range of post-spike voltage dynamics. We calculate the phase response curve using an extension of the adjoint method that accounts for the discontinuous post-spike reset rule. We find that both spikes and subthreshold fluctuations can jointly promote synchronization. The subthreshold contribution is strongest when the voltage exhibits a significant post-spike elevation in voltage, or plateau potential. Additionally, we show that the geometry of trajectories approaching the spiking threshold causes a "reset-induced shear" effect that can oppose synchrony in the presence of network asymmetry, despite having no effect on the phase-locking of symmetrically coupled pairs.

Integrate-and-fire network model of activity propagation from thalamus to cortex

The thalamus plays a crucial role in modulating the cortical activity underlying sensory and cognitive processes. In particular, recent experimental findings highlighted that the thalamus does not merely act as a binary gate for sensory stimuli, but rather participates to the processing of sensory information. Clarifying such thalamic influence on cortical dynamics is also important as the thalamus is the target of therapies such as DBS for Tourette patients. In this perspective, various computational models have been proposed in the last decades. However, a detailed description of the propagation of thalamic activity to the cortex is missing. Here we present a simple computational model of thalamocortical connectivity accounting for the propagation of activity from the thalamus to the cortex. The model includes both the single-neuron scale and the mesoscopic level of Local Field Potential (LFP) signals. Numerical simulations at both levels reproduce typical thalamocortical dynamics which are consistent with experimental measurements and robust to parameters changes. In particular, our model correctly reproduces locally generated rhythms as spindle oscillations in the thalamus and gamma oscillations in the cortex. Our model paves the way to deeper investigations of the thalamic influence on cortical dynamics, with and without sensory inputs or therapeutic electrical stimulation.

Contrasting the effects of adaptation and synaptic filtering on the timescales of dynamics in recurrent networks

Neural activity in awake behaving animals exhibits a vast range of timescales that can be several fold larger than the membrane time constant of individual neurons. Two types of mechanisms have been proposed to explain this conundrum. One possibility is that large timescales are generated by a network mechanism based on positive feedback, but this hypothesis requires fine-tuning of the strength or structure of the synaptic connections. A second possibility is that large timescales in the neural dynamics are inherited from large timescales of underlying biophysical processes, two prominent candidates being intrinsic adaptive ionic currents and synaptic transmission. How the timescales of adaptation or synaptic transmission influence the timescale of the network dynamics has however not been fully explored. To address this question, here we analyze large networks of randomly connected excitatory and inhibitory units with additional degrees of freedom that correspond to adaptation or synaptic filtering. We determine the fixed points of the systems, their stability to perturbations and the corresponding dynamical timescales. Furthermore, we apply dynamical mean field theory to study the temporal statistics of the activity in the fluctuating regime, and examine how the adaptation and synaptic timescales transfer from individual units to the whole population. Our overarching finding is that synaptic filtering and adaptation in single neurons have very different effects at the network level. Unexpectedly, the macroscopic network dynamics do not inherit the large timescale present in adaptive currents. In contrast, the timescales of network activity increase proportionally to the time constant of the synaptic filter. Altogether, our study demonstrates that the timescales of different biophysical processes have different effects on the network level, so that the slow processes within individual neurons do not necessarily induce slow activity in large recurrent neural networks.

Pulse-coupled mixed-mode oscillators: Cluster states and extreme noise sensitivity

Motivated by rhythms in the olfactory system of the brain, we investigate the synchronization of all-to-all pulse-coupled neuronal oscillators exhibiting various types of mixed-mode oscillations (MMOs) composed of sub-threshold oscillations (STOs) and action potentials ("spikes"). We focus particularly on the impact of the delay in the interaction. In the weak-coupling regime, we reduce the system to a Kuramoto-type equation with non-sinusoidal phase coupling and the associated Fokker-Planck equation. Its linear stability analysis identifies the appearance of various cluster states. Their type depends sensitively on the delay and the width of the pulses. Interestingly, long delays do not imply slow population rhythms, and the number of emerging clusters only loosely depends on the number of STOs. Direct simulations of the oscillator equations reveal that for quantitative agreement of the weak-coupling theory the coupling strength and the noise have to be extremely small. Even moderate noise leads to significant skipping of STO cycles, which can enhance the diffusion coefficient in the Fokker-Planck equation by two orders of magnitude. Introducing an effective diffusion coefficient extends the range of agreement significantly. Numerical simulations of the Fokker-Planck equation reveal bistability and solutions with oscillatory order parameters that result from nonlinear mode interactions. These are confirmed in simulations of the full spiking model.

Self-Consistent Scheme for Spike-Train Power Spectra in Heterogeneous Sparse Networks

Recurrent networks of spiking neurons can be in an asynchronous state characterized by low or absent cross-correlations and spike statistics which resemble those of cortical neurons. Although spatial correlations are negligible in this state, neurons can show pronounced temporal correlations in their spike trains that can be quantified by the autocorrelation function or the spike-train power spectrum. Depending on cellular and network parameters, correlations display diverse patterns (ranging from simple refractory-period effects and stochastic oscillations to slow fluctuations) and it is generally not well-understood how these dependencies come about. Previous work has explored how the single-cell correlations in a homogeneous network (excitatory and inhibitory integrate-and-fire neurons with nearly balanced mean recurrent input) can be determined numerically from an iterative single-neuron simulation. Such a scheme is based on the fact that every neuron is driven by the network noise (i.e., the input currents from all its presynaptic partners) but also contributes to the network noise, leading to a self-consistency condition for the input and output spectra. Here we first extend this scheme to homogeneous networks with strong recurrent inhibition and a synaptic filter, in which instabilities of the previous scheme are avoided by an averaging procedure. We then extend the scheme to heterogeneous networks in which (i) different neural subpopulations (e.g., excitatory and inhibitory neurons) have different cellular or connectivity parameters; (ii) the number and strength of the input connections are random (Erdős-Rényi topology) and thus different among neurons. In all heterogeneous cases, neurons are lumped in different classes each of which is represented by a single neuron in the iterative scheme; in addition, we make a Gaussian approximation of the input current to the neuron. These approximations seem to be justified over a broad range of parameters as indicated by comparison with simulation results of large recurrent networks. Our method can help to elucidate how network heterogeneity shapes the asynchronous state in recurrent neural networks.

Nanodiamonds-induced effects on neuronal firing of mouse hippocampal microcircuits

Fluorescent nanodiamonds (FND) are carbon-based nanomaterials that can efficiently incorporate optically active photoluminescent centers such as the nitrogen-vacancy complex, thus making them promising candidates as optical biolabels and drug-delivery agents. FNDs exhibit bright fluorescence without photobleaching combined with high uptake rate and low cytotoxicity. Focusing on FNDs interference with neuronal function, here we examined their effect on cultured hippocampal neurons, monitoring the whole network development as well as the electrophysiological properties of single neurons. We observed that FNDs drastically decreased the frequency of inhibitory (from 1.81 Hz to 0.86 Hz) and excitatory (from 1.61 to 0.68 Hz) miniature postsynaptic currents, and consistently reduced action potential (AP) firing frequency (by 36%), as measured by microelectrode arrays. On the contrary, bursts synchronization was preserved, as well as the amplitude of spontaneous inhibitory and excitatory events. Current-clamp recordings revealed that the ratio of neurons responding with AP trains of high-frequency (fast-spiking) versus neurons responding with trains of low-frequency (slow-spiking) was unaltered, suggesting that FNDs exerted a comparable action on neuronal subpopulations. At the single cell level, rapid onset of the somatic AP (“kink”) was drastically reduced in FND-treated neurons, suggesting a reduced contribution of axonal and dendritic components while preserving neuronal excitability.

Dichotomous Dynamics in E-I Networks with Strongly and Weakly Intra-connected Inhibitory Neurons

The interconnectivity between excitatory and inhibitory neural networks informs mechanisms by which rhythmic bursts of excitatory activity can be produced in the brain. One such mechanism, Pyramidal Interneuron Network Gamma (PING), relies primarily upon reciprocal connectivity between the excitatory and inhibitory networks, while also including intra-connectivity of inhibitory cells. The causal relationship between excitatory activity and the subsequent burst of inhibitory activity is of paramount importance to the mechanism and has been well studied. However, the role of the intra-connectivity of the inhibitory network, while important for PING, has not been studied in detail, as most analyses of PING simply assume that inhibitory intra-connectivity is strong enough to suppress subsequent firing following the initial inhibitory burst. In this paper we investigate the role that the strength of inhibitory intra-connectivity plays in determining the dynamics of PING-style networks. We show that networks with weak inhibitory intra-connectivity exhibit variations in burst dynamics of both the excitatory and inhibitory cells that are not obtained with strong inhibitory intra-connectivity. Networks with weak inhibitory intra-connectivity exhibit excitatory rhythmic bursts with weak excitatory-to-inhibitory synapses for which classical PING networks would show no rhythmic activity. Additionally, variations in dynamics of these networks as the excitatory-to-inhibitory synaptic weight increases illustrates the important role that consistent pattern formation in the inhibitory cells serves in maintaining organized and periodic excitatory bursts. Finally, motivated by these results and the known diversity of interneurons, we show that a PING-style network with two inhibitory subnetworks, one strongly intra-connected and one weakly intra-connected, exhibits organized and periodic excitatory activity over a larger parameter regime than networks with a homogeneous inhibitory population. Taken together, these results serve to better articulate the role of inhibitory intra-connectivity in generating PING-like rhythms, while also revealing how heterogeneity amongst inhibitory synapses might make such rhythms more robust to a variety of network parameters.

Inhibition of Metabotropic Glutamate Receptor Subtype 1 Alters the Excitability of the Commissural Pyramidal Neuron in the Rat Anterior Cingulate Cortex after Chronic Constriction Injury to the Sciatic Nerve

Background

Inhibition of the metabotropic glutamate receptor subtype 1 in the anterior cingulate cortex has an analgesic effect during sustained nociceptive hypersensitivity. However, the specific changes in different subtypes of anterior cingulate cortex layer 5 pyramidal neurons, as well as the distinct effect of metabotropic glutamate receptor subtype 1 inhibition on different neuronal subtypes, have not been well studied.

Methods

Retrograde labeling combined with immunofluorescence, whole cell clamp recording, and behavioral tests combined with RNA interference were performed in a rat model of chronic constriction injury to the sciatic nerve.

Results

Commissural layer 5 pyramidal neurons (projecting to the contralateral cortex) existed in the anterior cingulate cortex. The voltage-gated potassium channel subunit 2–mediated current in these neurons were substantially reduced after chronic constriction injury (current densities at +30 mV for the sham, and chronic constriction injury neurons were [mean ± SD] 10.22 ± 3.42 pA/pF vs. 5.58 ± 2.71 pA/pF, respectively; n = 11; P < 0.01), which increased the spike width and fast afterhyperpolarization potential, resulting in hyperexcitability. Inhibition of metabotropic glutamate receptor subtype 1 alleviated the down-regulation of voltage-gated potassium channel subunit 2 currents (current density increased by 8.11 ± 3.22 pA/pF; n = 7; P < 0.01). Furthermore, knockdown of voltage-gated potassium channel subunit 2 current in the commissural neurons attenuated the analgesic effect of metabotropic glutamate receptor subtype 1 inhibition (n = 6 rats; P < 0.05).

Conclusions

The effect of metabotropic glutamate receptor subtype 1 inhibition on commissural anterior cingulate cortex layer 5 pyramidal neurons is likely different with the modification of previously studied hyperpolarization-activated/cyclic nucleotide-gated channel-dependent neurons but relies on the alteration of voltage-gated potassium channel subunit 2 currents. These results will contribute to a better understanding of the therapeutic role of metabotropic glutamate receptor subtype 1 in chronic pain.

Effect of Phase Response Curve Shape and Synaptic Driving Force on Synchronization of Coupled Neuronal Oscillators

The role of the phase response curve (PRC) shape on the synchrony of synaptically coupled oscillating neurons is examined. If the PRC is independent of the phase, because of the synaptic form of the coupling, synchrony is found to be stable for both excitatory and inhibitory coupling at all rates, whereas the antisynchrony becomes stable at low rates. A faster synaptic rise helps extend the stability of antisynchrony to higher rates. If the PRC is not constant but has a profile like that of a leaky integrate-and-fire model, then, in contrast to the earlier reports that did not include the voltage effects, mutual excitation could lead to stable synchrony provided the synaptic reversal potential is below the voltage level the neuron would have reached in the absence of the interaction and threshold reset. This level is controlled by the applied current and the leakage parameters. Such synchrony is contingent on significant phase response (that would result, for example, by a sharp PRC jump) occurring during the synaptic rising phase. The rising phase, however, does not contribute significantly if it occurs before the voltage spike reaches its peak. Then a stable near-synchronous state can still exist between type 1 PRC neurons if the PRC shows a left skewness in its shape. These results are examined comprehensively using perfect integrate-and-fire, leaky integrate-and-fire, and skewed PRC shapes under the assumption of the weakly coupled oscillator theory applied to synaptically coupled neuron models.

Modulation of neuronal dynamic range using two different adaptation mechanisms

The capability of neurons to discriminate between intensity of external stimulus is measured by its dynamic range. A larger dynamic range indicates a greater probability of neuronal survival. In this study, the potential roles of adaptation mechanisms (ion currents) in modulating neuronal dynamic range were numerically investigated. Based on the adaptive exponential integrate-and-fire model, which includes two different adaptation mechanisms, i.e. subthreshold and suprathreshold (spike-triggered) adaptation, our results reveal that the two adaptation mechanisms exhibit rather different roles in regulating neuronal dynamic range. Specifically, subthreshold adaptation acts as a negative factor that observably decreases the neuronal dynamic range, while suprathreshold adaptation has little influence on the neuronal dynamic range. Moreover, when stochastic noise was introduced into the adaptation mechanisms, the dynamic range was apparently enhanced, regardless of what state the neuron was in, e.g. adaptive or non-adaptive. Our model results suggested that the neuronal dynamic range can be differentially modulated by different adaptation mechanisms. Additionally, noise was a non-ignorable factor, which could effectively modulate the neuronal dynamic range.

Extending Integrate-and-Fire Model Neurons to Account for the Effects of Weak Electric Fields and Input Filtering Mediated by the Dendrite

Transcranial brain stimulation and evidence of ephaptic coupling have recently sparked strong interests in understanding the effects of weak electric fields on the dynamics of brain networks and of coupled populations of neurons. The collective dynamics of large neuronal populations can be efficiently studied using single-compartment (point) model neurons of the integrate-and-fire (IF) type as their elements. These models, however, lack the dendritic morphology required to biophysically describe the effect of an extracellular electric field on the neuronal membrane voltage. Here, we extend the IF point neuron models to accurately reflect morphology dependent electric field effects extracted from a canonical spatial “ball-and-stick” (BS) neuron model. Even in the absence of an extracellular field, neuronal morphology by itself strongly affects the cellular response properties. We, therefore, derive additional components for leaky and nonlinear IF neuron models to reproduce the subthreshold voltage and spiking dynamics of the BS model exposed to both fluctuating somatic and dendritic inputs and an extracellular electric field. We show that an oscillatory electric field causes spike rate resonance, or equivalently, pronounced spike to field coherence. Its resonance frequency depends on the location of the synaptic background inputs. For somatic inputs the resonance appears in the beta and gamma frequency range, whereas for distal dendritic inputs it is shifted to even higher frequencies. Irrespective of an external electric field, the presence of a dendritic cable attenuates the subthreshold response at the soma to slowly-varying somatic inputs while implementing a low-pass filter for distal dendritic inputs. Our point neuron model extension is straightforward to implement and is computationally much more efficient compared to the original BS model. It is well suited for studying the dynamics of large populations of neurons with heterogeneous dendritic morphology with (and without) the influence of weak external electric fields.

How extracellular electric fields—as generated endogenously or through transcranial brain stimulation—affect the dynamics of neuronal populations is of great interest but not well understood. To study neuronal activity at the network level single-compartment neuron models have been proven very successful, because of their computational efficiency and analytical tractability. Unfortunately, these models lack the dendritic morphology to biophysically account for the effects of electric fields, and for changes in synaptic integration due to morphology alone. Here, we consider a canonical, spatially extended model neuron and characterize its responses to fluctuating synaptic input as well as an oscillatory, weak electric field. In order to accurately reproduce these responses we analytically derive an extension for the popular integrate-and-fire point neuron models. We show that the dendritic cable acts as a filter for the synaptic input current, which depends on the input location, and that an electric field modulates the neuronal spike rate strongest at a certain (preferred) field frequency. These phenomena can be successfully reproduced using integrate-and-fire models, extended by a small number of components that are straightforward to implement. The extended point models are thus well suited for studying populations of coupled neurons with different morphology, exposed to extracellular electric fields.

The Impact of Structural Heterogeneity on Excitation-Inhibition Balance in Cortical Networks

Models of cortical dynamics often assume a homogeneous connectivity structure. However, we show that heterogeneous input connectivity can prevent the dynamic balance between excitation and inhibition, a hallmark of cortical dynamics, and yield unrealistically sparse and temporally regular firing. Anatomically based estimates of the connectivity of layer 4 (L4) rat barrel cortex and numerical simulations of this circuit indicate that the local network possesses substantial heterogeneity in input connectivity, sufficient to disrupt excitation-inhibition balance. We show that homeostatic plasticity in inhibitory synapses can align the functional connectivity to compensate for structural heterogeneity. Alternatively, spike-frequency adaptation can give rise to a novel state in which local firing rates adjust dynamically so that adaptation currents and synaptic inputs are balanced. This theory is supported by simulations of L4 barrel cortex during spontaneous and stimulus-evoked conditions. Our study shows how synaptic and cellular mechanisms yield fluctuation-driven dynamics despite structural heterogeneity in cortical circuits.

•

Structural heterogeneity threatens the dynamic balance of excitation and inhibition

•

Reconstruction of cortical networks reveals significant structural heterogeneity

•

Spike-frequency adaptation can act locally to facilitate global balance

•

Inhibitory homeostatic plasticity can compensate for structural imbalance

Intrinsic Cellular Properties and Connectivity Density Determine Variable Clustering Patterns in Randomly Connected Inhibitory Neural Networks

The plethora of inhibitory interneurons in the hippocampus and cortex play a pivotal role in generating rhythmic activity by clustering and synchronizing cell firing. Results of our simulations demonstrate that both the intrinsic cellular properties of neurons and the degree of network connectivity affect the characteristics of clustered dynamics exhibited in randomly connected, heterogeneous inhibitory networks. We quantify intrinsic cellular properties by the neuron's current-frequency relation (IF curve) and Phase Response Curve (PRC), a measure of how perturbations given at various phases of a neurons firing cycle affect subsequent spike timing. We analyze network bursting properties of networks of neurons with Type I or Type II properties in both excitability and PRC profile; Type I PRCs strictly show phase advances and IF curves that exhibit frequencies arbitrarily close to zero at firing threshold while Type II PRCs display both phase advances and delays and IF curves that have a non-zero frequency at threshold. Type II neurons whose properties arise with or without an M-type adaptation current are considered. We analyze network dynamics under different levels of cellular heterogeneity and as intrinsic cellular firing frequency and the time scale of decay of synaptic inhibition are varied. Many of the dynamics exhibited by these networks diverge from the predictions of the interneuron network gamma (ING) mechanism, as well as from results in all-to-all connected networks. Our results show that randomly connected networks of Type I neurons synchronize into a single cluster of active neurons while networks of Type II neurons organize into two mutually exclusive clusters segregated by the cells' intrinsic firing frequencies. Networks of Type II neurons containing the adaptation current behave similarly to networks of either Type I or Type II neurons depending on network parameters; however, the adaptation current creates differences in the cluster dynamics compared to those in networks of Type I or Type II neurons. To understand these results, we compute neuronal PRCs calculated with a perturbation matching the profile of the synaptic current in our networks. Differences in profiles of these PRCs across the different neuron types reveal mechanisms underlying the divergent network dynamics.

Transition between Functional Regimes in an Integrate-And-Fire Network Model of the Thalamus

The thalamus is a key brain element in the processing of sensory information. During the sleep and awake states, this brain area is characterized by the presence of two distinct dynamical regimes: in the sleep state activity is dominated by spindle oscillations (7 − 15 Hz) weakly affected by external stimuli, while in the awake state the activity is primarily driven by external stimuli. Here we develop a simple and computationally efficient model of the thalamus that exhibits two dynamical regimes with different information-processing capabilities, and study the transition between them. The network model includes glutamatergic thalamocortical (TC) relay neurons and GABAergic reticular (RE) neurons described by adaptive integrate-and-fire models in which spikes are induced by either depolarization or hyperpolarization rebound. We found a range of connectivity conditions under which the thalamic network composed by these neurons displays the two aforementioned dynamical regimes. Our results show that TC-RE loops generate spindle-like oscillations and that a minimum level of clustering (i.e. local connectivity density) in the RE-RE connections is necessary for the coexistence of the two regimes. We also observe that the transition between the two regimes occurs when the external excitatory input on TC neurons (mimicking sensory stimulation) is large enough to cause a significant fraction of them to switch from hyperpolarization-rebound-driven firing to depolarization-driven firing. Overall, our model gives a novel and clear description of the role that the two types of neurons and their connectivity play in the dynamical regimes observed in the thalamus, and in the transition between them. These results pave the way for the development of efficient models of the transmission of sensory information from periphery to cortex.

Information Coding through Adaptive Gating of Synchronized Thalamic Bursting

It has been posited that the regulation of burst/tonic firing in the thalamus could function as a mechanism for controlling not only how much but what kind of information is conveyed to downstream cortical targets. Yet how this gating mechanism is adaptively modulated on fast timescales by ongoing sensory inputs in rich sensory environments remains unknown. Using single-unit recordings in the rat vibrissa thalamus (VPm), we found that the degree of bottom-up adaptation modulated thalamic burst/tonic firing as well as the synchronization of bursting across the thalamic population along a continuum for which the extremes facilitate detection or discrimination of sensory inputs. Optogenetic control of baseline membrane potential in thalamus further suggests that this regulation may result from an interplay between adaptive changes in thalamic membrane potential and synaptic drive from inputs to thalamus, setting the stage for an intricate control strategy upon which cortical computation is built.

Bud extracts from Tilia tomentosa Moench inhibit hippocampal neuronal firing through GABAA and benzodiazepine receptors activation

ETHNOPHARMACOLOGICAL RELEVANCE

Tilia tomentosa Moench bud extracts (TTBEs) is used in traditional medicine for centuries as sedative compound. Different plants belonging to the Tilia genus have shown their efficacy in the treatment of anxiety but still little is known about the mechanism of action of their bud extracts.

AIM OF THE STUDY

To evaluate the action of TTBEs as anxiolytic and sedative compound on in vitro hippocampal neurons.

MATERIAL AND METHODS

The anxiolytic effect of TTBEs was assayed by testing the effects of these compounds on GABAA receptor-activated chloride current of hippocampal neurons by means of the patch-clamp technique and microelectrode-arrays (MEAs).

RESULTS

TTBEs acutely administered on mouse hippocampal neurons, activated a chloride current comparable to that measured in the presence of GABA (100 µM). Bicuculline (100 µM) and picrotoxin (100 µM) blocked about 90% of this current, while the remaining 10% was blocked by adding the benzodiazepine (BDZ) antagonist flumazenil (30 µM). Flumazenil alone blocked nearly 60% of the TTBEs activated current, suggesting that TTBEs binds to both GABAA and BDZ receptor sites. Application of high-doses of TTBEs on spontaneous active hippocampal neurons grown for 3 weeks on MEAs blocked the synchronous activity of these neurons. The effects were mimicked by GABA and prevented by picrotoxin (100µM) and flumazenil (30 µM). At minimal doses, TTBEs reduced the frequency of synchronized bursts and increased the cross-correlation index of synchronized neuronal firing.

CONCLUSIONS

Our data suggest that TTBEs mimics GABA and BDZ agonists by targeting hippocampal GABAergic synapses and inhibiting network excitability by increasing the strength of inhibitory synaptic outputs. Our results contribute toward the validation of TTBEs as effective sedative and anxiolytic compound.

Cav1.3 Channels as Key Regulators of Neuron-Like Firings and Catecholamine Release in Chromaffin Cells

Neuronal and neuroendocrine L-type calcium channels (Cav1.2, Cav1.3) open readily at relatively low membrane potentials and allow Ca(2+) to enter the cells near resting potentials. In this way, Cav1.2 and Cav1.3 shape the action potential waveform, contribute to gene expression, synaptic plasticity, neuronal differentiation, hormone secretion and pacemaker activity. In the chromaffin cells (CCs) of the adrenal medulla, Cav1.3 is highly expressed and is shown to support most of the pacemaking current that sustains action potential (AP) firings and part of the catecholamine secretion. Cav1.3 forms Ca(2+)-nanodomains with the fast inactivating BK channels and drives the resting SK currents. These latter set the inter-spike interval duration between consecutive spikes during spontaneous firing and the rate of spike adaptation during sustained depolarizations. Cav1.3 plays also a primary role in the switch from "tonic" to "burst" firing that occurs in mouse CCs when either the availability of voltage-gated Na channels (Nav) is reduced or the β2 subunit featuring the fast inactivating BK channels is deleted. Here, we discuss the functional role of these "neuron-like" firing modes in CCs and how Cav1.3 contributes to them. The open issue is to understand how these novel firing patterns are adapted to regulate the quantity of circulating catecholamines during resting condition or in response to acute and chronic stress.

Effect of phase response curve skew on synchronization with and without conduction delays

A central problem in cortical processing including sensory binding and attentional gating is how neurons can synchronize their responses with zero or near-zero time lag. For a spontaneously firing neuron, an input from another neuron can delay or advance the next spike by different amounts depending upon the timing of the input relative to the previous spike. This information constitutes the phase response curve (PRC). We present a simple graphical method for determining the effect of PRC shape on synchronization tendencies and illustrate it using type 1 PRCs, which consist entirely of advances (delays) in response to excitation (inhibition). We obtained the following generic solutions for type 1 PRCs, which include the pulse-coupled leaky integrate and fire model. For pairs with mutual excitation, exact synchrony can be stable for strong coupling because of the stabilizing effect of the causal limit region of the PRC in which an input triggers a spike immediately upon arrival. However, synchrony is unstable for short delays, because delayed inputs arrive during a refractory period and cannot trigger an immediate spike. Right skew destabilizes antiphase and enables modes with time lags that grow as the conduction delay is increased. Therefore, right skew favors near synchrony at short conduction delays and a gradual transition between synchrony and antiphase for pairs coupled by mutual excitation. For pairs with mutual inhibition, zero time lag synchrony is stable for conduction delays ranging from zero to a substantial fraction of the period for pairs. However, for right skew there is a preferred antiphase mode at short delays. In contrast to mutual excitation, left skew destabilizes antiphase for mutual inhibition so that synchrony dominates at short delays as well. These pairwise synchronization tendencies constrain the synchronization properties of neurons embedded in larger networks.

Patterns of interval correlations in neural oscillators with adaptation

Neural firing is often subject to negative feedback by adaptation currents. These currents can induce strong correlations among the time intervals between spikes. Here we study analytically the interval correlations of a broad class of noisy neural oscillators with spike-triggered adaptation of arbitrary strength and time scale. Our weak-noise theory provides a general relation between the correlations and the phase-response curve (PRC) of the oscillator, proves anti-correlations between neighboring intervals for adapting neurons with type I PRC and identifies a single order parameter that determines the qualitative pattern of correlations. Monotonically decaying or oscillating correlation structures can be related to qualitatively different voltage traces after spiking, which can be explained by the phase plane geometry. At high firing rates, the long-term variability of the spike train associated with the cumulative interval correlations becomes small, independent of model details. Our results are verified by comparison with stochastic simulations of the exponential, leaky, and generalized integrate-and-fire models with adaptation.

How adaptation currents change threshold, gain, and variability of neuronal spiking

Many types of neurons exhibit spike rate adaptation, mediated by intrinsic slow K(+) currents, which effectively inhibit neuronal responses. How these adaptation currents change the relationship between in vivo like fluctuating synaptic input, spike rate output, and the spike train statistics, however, is not well understood. In this computational study we show that an adaptation current that primarily depends on the subthreshold membrane voltage changes the neuronal input-output relationship (I-O curve) subtractively, thereby increasing the response threshold, and decreases its slope (response gain) for low spike rates. A spike-dependent adaptation current alters the I-O curve divisively, thus reducing the response gain. Both types of an adaptation current naturally increase the mean interspike interval (ISI), but they can affect ISI variability in opposite ways. A subthreshold current always causes an increase of variability while a spike-triggered current decreases high variability caused by fluctuation-dominated inputs and increases low variability when the average input is large. The effects on I-O curves match those caused by synaptic inhibition in networks with asynchronous irregular activity, for which we find subtractive and divisive changes caused by external and recurrent inhibition, respectively. Synaptic inhibition, however, always increases the ISI variability. We analytically derive expressions for the I-O curve and ISI variability, which demonstrate the robustness of our results. Furthermore, we show how the biophysical parameters of slow K(+) conductances contribute to the two different types of an adaptation current and find that Ca(2+)-activated K(+) currents are effectively captured by a simple spike-dependent description, while muscarine-sensitive or Na(+)-activated K(+) currents show a dominant subthreshold component.

Firing properties of entorhinal cortex neurons and early alterations in an Alzheimer's disease transgenic model

The entorhinal cortex (EC) is divided into medial (MEC) and lateral (LEC) anatomical areas, and layer II neurons of these two regions project to granule cells of the dentate gyrus through the medial and lateral perforant pathways (MPP and LPP), respectively. Stellate cells (SCs) represent the main neurons constituting the MPP inputs, while fan cells (FCs) represent the main LPP inputs. Here, we first characterized the excitability properties of SCs and FCs in adult wild-type (WT) mouse brain. Our data indicate that, during sustained depolarization, action potentials (APs) generated by SCs exhibit increased fast afterhyperpolarization and overshoot, making them able to fire at higher frequencies and to exhibit higher spike frequency adaptation (SFA) than FCs. Since the EC is one of the earliest brain regions affected during Alzheimer's disease (AD) progression, we compared SCs and FCs firing in 4-month-old WT and transgenic Tg2576 mice, a well-established AD mouse model. Tg2576-SCs displayed a slight increase in firing frequency during mild depolarization but otherwise normal excitability properties during higher stimulations. On the contrary, Tg2576-FCs exhibited a decreased firing frequency during mild and higher depolarizations, as well as an increased SFA. Our data identify the FCs as a neuronal population particularly sensitive to early pathological effects of chronic accumulation of APP-derived peptides, as it occurs in Tg2576 mice. As FCs represent the major input of sensory information to the hippocampus during memory acquisition, early alterations in their excitability profile could significantly contribute to the onset of cognitive decline in AD.

Interaction function of oscillating coupled neurons

Large scale simulations of electrically coupled neuronal oscillators often employ the phase coupled oscillator paradigm to understand and predict network behavior. We study the nature of the interaction between such coupled oscillators using weakly coupled oscillator theory. By employing piecewise linear approximations for phase response curves and voltage time courses and parametrizing their shapes, we compute the interaction function for all such possible shapes and express it in terms of discrete Fourier modes. We find that reasonably good approximation is achieved with four Fourier modes that comprise of both sine and cosine terms.

Adaptation controls synchrony and cluster states of coupled threshold-model neurons

We analyze zero-lag and cluster synchrony of delay-coupled nonsmooth dynamical systems by extending the master stability approach, and apply this to networks of adaptive threshold-model neurons. For a homogeneous population of excitatory and inhibitory neurons we find (i) that subthreshold adaptation stabilizes or destabilizes synchrony depending on whether the recurrent synaptic excitatory or inhibitory couplings dominate, and (ii) that synchrony is always unstable for networks with balanced recurrent synaptic inputs. If couplings are not too strong, synchronization properties are similar for very different coupling topologies, i.e., random connections or spatial networks with localized connectivity. We generalize our approach for two subpopulations of neurons with nonidentical local dynamics, including bursting, for which activity-based adaptation controls the stability of cluster states, independent of a specific coupling topology.

Spike width and frequency alter stability of phase-locking in electrically coupled neurons

The stability of phase-locked states of electrically coupled type-1 phase response curve neurons is studied using piecewise linear formulations for their voltage profile and phase response curves. We find that at low frequency and/or small spike width, synchrony is stable, and antisynchrony unstable. At high frequency and/or large spike width, these phase-locked states switch their stability. Increasing the ratio of spike width to spike height causes the antisynchronous state to transition into a stable synchronous state. We compute the interaction function and the boundaries of stability of both these phase-locked states, and present analytical expressions for them. We also study the effect of phase response curve skewness on the boundaries of synchrony and antisynchrony.

Effect of sharp jumps at the edges of phase response curves on synchronization of electrically coupled neuronal oscillators

We study synchronization phenomenon of coupled neuronal oscillators using the theory of weakly coupled oscillators. The role of sudden jumps in the phase response curve profiles found in some experimental recordings and models on the ability of coupled neurons to exhibit synchronous and antisynchronous behavior is investigated, when the coupling between the neurons is electrical. The level of jumps in the phase response curve at either end, spike width and frequency of voltage time course of the coupled neurons are parameterized using piecewise linear functional forms, and the conditions for stable synchrony and stable antisynchrony in terms of those parameters are computed analytically. The role of the peak position of the phase response curve on phase-locking is also investigated.

How adaptation shapes spike rate oscillations in recurrent neuronal networks

Neural mass signals from in-vivo recordings often show oscillations with frequencies ranging from <1 to 100 Hz. Fast rhythmic activity in the beta and gamma range can be generated by network-based mechanisms such as recurrent synaptic excitation-inhibition loops. Slower oscillations might instead depend on neuronal adaptation currents whose timescales range from tens of milliseconds to seconds. Here we investigate how the dynamics of such adaptation currents contribute to spike rate oscillations and resonance properties in recurrent networks of excitatory and inhibitory neurons. Based on a network of sparsely coupled spiking model neurons with two types of adaptation current and conductance-based synapses with heterogeneous strengths and delays we use a mean-field approach to analyze oscillatory network activity. For constant external input, we find that spike-triggered adaptation currents provide a mechanism to generate slow oscillations over a wide range of adaptation timescales as long as recurrent synaptic excitation is sufficiently strong. Faster rhythms occur when recurrent inhibition is slower than excitation and oscillation frequency increases with the strength of inhibition. Adaptation facilitates such network-based oscillations for fast synaptic inhibition and leads to decreased frequencies. For oscillatory external input, adaptation currents amplify a narrow band of frequencies and cause phase advances for low frequencies in addition to phase delays at higher frequencies. Our results therefore identify the different key roles of neuronal adaptation dynamics for rhythmogenesis and selective signal propagation in recurrent networks.

Ca(V)1.3-driven SK channel activation regulates pacemaking and spike frequency adaptation in mouse chromaffin cells

Mouse chromaffin cells (MCCs) fire spontaneous action potentials (APs) at rest. Ca(v)1.3 L-type calcium channels sustain the pacemaker current, and their loss results in depolarized resting potentials (V(rest)), spike broadening, and remarkable switches into depolarization block after BayK 8644 application. A functional coupling between Ca(v)1.3 and BK channels has been reported but cannot fully account for the aforementioned observations. Here, using Ca(v)1.3(-/-) mice, we investigated the role of Ca(v)1.3 on SK channel activation and how this functional coupling affects the firing patterns induced by sustained current injections. MCCs express SK1-3 channels whose tonic currents are responsible for the slow irregular firing observed at rest. Percentage of frequency increase induced by apamin was found inversely correlated to basal firing frequency. Upon stimulation, MCCs build-up Ca(v)1.3-dependent SK currents during the interspike intervals that lead to a notable degree of spike frequency adaptation (SFA). The major contribution of Ca(v)1.3 to the subthreshold Ca(2+) charge during an AP-train rather than a specific molecular coupling to SK channels accounts for the reduced SFA of Ca(v)1.3(-/-) MCCs. Low adaptation ratios due to reduced SK activation associated with Ca(v)1.3 deficiency prevent the efficient recovery of Na(V) channels from inactivation. This promotes a rapid decline of AP amplitudes and facilitates early onset of depolarization block following prolonged stimulation. Thus, besides serving as pacemaker, Ca(v)1.3 slows down MCC firing by activating SK channels that maintain Na(V) channel availability high enough to preserve stable AP waveforms, even upon high-frequency stimulation of chromaffin cells during stress responses.