Synchrony of neuronal oscillations controlled by GABAergic reversal potentials

Non-monotonic effects of GABAergic synaptic inputs on neuronal firing

GABA is generally known as the principal inhibitory neurotransmitter in the nervous system, usually acting by hyperpolarizing membrane potential. However, GABAergic currents sometimes exhibit non-inhibitory effects, depending on the brain region, developmental stage or pathological condition. Here, we investigate the diverse effects of GABA on the firing rate of several single neuron models, using both analytical calculations and numerical simulations. We find that GABAergic synaptic conductance and output firing rate exhibit three qualitatively different regimes as a function of GABA reversal potential, EGABA: monotonically decreasing for sufficiently low EGABA (inhibitory), monotonically increasing for EGABA above firing threshold (excitatory); and a non-monotonic region for intermediate values of EGABA. In the non-monotonic regime, small GABA conductances have an excitatory effect while large GABA conductances show an inhibitory effect. We provide a phase diagram of different GABAergic effects as a function of GABA reversal potential and glutamate conductance. We find that noisy inputs increase the range of EGABA for which the non-monotonic effect can be observed. We also construct a micro-circuit model of striatum to explain observed effects of GABAergic fast spiking interneurons on spiny projection neurons, including non-monotonicity, as well as the heterogeneity of the effects. Our work provides a mechanistic explanation of paradoxical effects of GABAergic synaptic inputs, with implications for understanding the effects of GABA in neural computation and development.

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.

Distinct effects of heterogeneity and noise on gamma oscillation in a model of neuronal network with different reversal potential

Gamma oscillation is crucial in brain functions such as attentional selection, and is inextricably linked to both heterogeneity and noise (or so-called stochastic fluctuation) in neuronal networks. However, under coexistence of these factors, it has not been clarified how the synaptic reversal potential modulates the entraining of gamma oscillation. Here we show distinct effects of heterogeneity and noise in a population of modified theta neurons randomly coupled via GABAergic synapses. By introducing the Fokker-Planck equation and circular cumulants, we derive a set of two-cumulant macroscopic equations. In bifurcation analyses, we find a stabilizing effect of heterogeneity and a nontrivial effect of noise that results in promoting, diminishing, and shifting the oscillatory region, and is largely dependent on the reversal potential of GABAergic synapses. These findings are verified by numerical simulations of a finite-size neuronal network. Our results reveal that slight changes in reversal potential and magnitude of stochastic fluctuations can lead to immediate control of gamma oscillation, which would results in complex spatio-temporal dynamics for attentional selection and recognition.

The effects of chloride dynamics on substantia nigra pars reticulata responses to pallidal and striatal inputs

As a rodent basal ganglia (BG) output nucleus, the substantia nigra pars reticulata (SNr) is well positioned to impact behavior. SNr neurons receive GABAergic inputs from the striatum (direct pathway) and globus pallidus (GPe, indirect pathway). Dominant theories of action selection rely on these pathways' inhibitory actions. Yet, experimental results on SNr responses to these inputs are limited and include excitatory effects. Our study combines experimental and computational work to characterize, explain, and make predictions about these pathways. We observe diverse SNr responses to stimulation of SNr-projecting striatal and GPe neurons, including biphasic and excitatory effects, which our modeling shows can be explained by intracellular chloride processing. Our work predicts that ongoing GPe activity could tune the SNr operating mode, including its responses in decision-making scenarios, and GPe output may modulate synchrony and low-frequency oscillations of SNr neurons, which we confirm using optogenetic stimulation of GPe terminals within the SNr.

Functional characterization of novel bumetanide derivatives for epilepsy treatment

Temporal lobe epilepsy (TLE) is the most common type of focal epilepsies, affecting approximately 35 million people worldwide. Despite the introduction of numerous novel antiepileptic drugs during the last decades, the proportion of patients with therapy-resistant TLE is still high. As an impaired cellular chloride homeostasis appears involved in disease pathophysiology, bumetanide, an antagonist to Na-K-Cl cotransporters, gained interest as potential therapeutic option. However, bumetanide induces a strong diuretic effect and displays poor penetration across the blood-brain barrier (BBB). To reduce these unwanted effects, we modified the already described BUM690 by exchanging the allyl-into a trifluoro-ethyl group to yield BUM532. Furthermore, we exchanged the nitrogen for oxygen in the trifluoro-ethyl group to yield BUM97. In the intrahippocampal kainic acid mouse model of TLE BUM532 ± phenobarbital (PB), bumetanide ± PB and PB alone significantly reduced hippocampal paroxysmal discharges (HPDs) but not spike trains. By contrast, treatment with BUM97 suppressed HPDs as well as spike trains dose-dependently, more pronounced compared to the other tested compounds and exerted a synergistic anticonvulsant effect with PB. Moreover, at higher doses BUM97 achieved long-lasting reduction of spike trains. In pentylenetetrazole-induced acute seizures only BUM532 combined with a sub-effective dose of PB increased the seizure threshold. No diuretic effects were observed at any dose of the three derivatives. Our data demonstrate the successful optimization of the pharmacological profile of bumetanide and the potential of the improved derivative BUM97 for the treatment of therapy-resistant TLE, in particular in combinatorial drug regimens with a GABA mimetic.

Phase response theory explains cluster formation in sparsely but strongly connected inhibitory neural networks and effects of jitter due to sparse connectivity

We show how to predict whether a neural network will exhibit global synchrony (a one-cluster state) or a two-cluster state based on the assumption of pulsatile coupling and critically dependent upon the phase response curve (PRC) generated by the appropriate perturbation from a partner cluster. Our results hold for a monotonically increasing (meaning longer delays as the phase increases) PRC, which likely characterizes inhibitory fast-spiking basket and cortical low-threshold-spiking interneurons in response to strong inhibition. Conduction delays stabilize synchrony for this PRC shape, whereas they destroy two-cluster states, the former by avoiding a destabilizing discontinuity and the latter by approaching it. With conduction delays, stronger coupling strength can promote a one-cluster state, so the weak coupling limit is not applicable here. We show how jitter can destabilize global synchrony but not a two-cluster state. Local stability of global synchrony in an all-to-all network does not guarantee that global synchrony can be observed in an appropriately scaled sparsely connected network; the basin of attraction can be inferred from the PRC and must be sufficiently large. Two-cluster synchrony is not obviously different from one-cluster synchrony in the presence of noise and may be the actual substrate for oscillations observed in the local field potential (LFP) and the electroencephalogram (EEG) in situations where global synchrony is not possible. Transitions between cluster states may change the frequency of the rhythms observed in the LFP or EEG. Transitions between cluster states within an inhibitory subnetwork may allow more effective recruitment of pyramidal neurons into the network rhythm. NEW & NOTEWORTHY We show that jitter induced by sparse connectivity can destabilize global synchrony but not a two-cluster state with two smaller clusters firing alternately. On the other hand, conduction delays stabilize synchrony and destroy two-cluster states. These results hold if each cluster exhibits a phase response curve similar to one that characterizes fast-spiking basket and cortical low-threshold-spiking cells for strong inhibition. Either a two-cluster or a one-cluster state might provide the oscillatory substrate for neural computations.

Developmental Emergence of Sparse Coding: A Dynamic Systems Approach

During neocortical development, network activity undergoes a dramatic transition from largely synchronized, so-called cluster activity, to a relatively sparse pattern around the time of eye-opening in rodents. Biophysical mechanisms underlying this sparsification phenomenon remain poorly understood. Here, we present a dynamic systems modeling study of a developing neural network that provides the first mechanistic insights into sparsification. We find that the rest state of immature networks is strongly affected by the dynamics of a transient, unstable state hidden in their firing activities, allowing these networks to either be silent or generate large cluster activity. We address how, and which, specific developmental changes in neuronal and synaptic parameters drive sparsification. We also reveal how these changes refine the information processing capabilities of an in vivo developing network, mainly by showing a developmental reduction in the instability of network’s firing activity, an effective availability of inhibition-stabilized states, and an emergence of spontaneous attractors and state transition mechanisms. Furthermore, we demonstrate the key role of GABAergic transmission and depressing glutamatergic synapses in governing the spatiotemporal evolution of cluster activity. These results, by providing a strong link between experimental observations and model behavior, suggest how adult sparse coding networks may emerge developmentally.

Dopamine Neurons Change the Type of Excitability in Response to Stimuli

The dynamics of neuronal excitability determine the neuron's response to stimuli, its synchronization and resonance properties and, ultimately, the computations it performs in the brain. We investigated the dynamical mechanisms underlying the excitability type of dopamine (DA) neurons, using a conductance-based biophysical model, and its regulation by intrinsic and synaptic currents. Calibrating the model to reproduce low frequency tonic firing results in N-methyl-D-aspartate (NMDA) excitation balanced by γ-Aminobutyric acid (GABA)-mediated inhibition and leads to type I excitable behavior characterized by a continuous decrease in firing frequency in response to hyperpolarizing currents. Furthermore, we analyzed how excitability type of the DA neuron model is influenced by changes in the intrinsic current composition. A subthreshold sodium current is necessary for a continuous frequency decrease during application of a negative current, and the low-frequency "balanced" state during simultaneous activation of NMDA and GABA receptors. Blocking this current switches the neuron to type II characterized by the abrupt onset of repetitive firing. Enhancing the anomalous rectifier Ih current also switches the excitability to type II. Key characteristics of synaptic conductances that may be observed in vivo also change the type of excitability: a depolarized γ-Aminobutyric acid receptor (GABAR) reversal potential or co-activation of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) leads to an abrupt frequency drop to zero, which is typical for type II excitability. Coactivation of N-methyl-D-aspartate receptors (NMDARs) together with AMPARs and GABARs shifts the type I/II boundary toward more hyperpolarized GABAR reversal potentials. To better understand how altering each of the aforementioned currents leads to changes in excitability profile of DA neuron, we provide a thorough dynamical analysis. Collectively, these results imply that type I excitability in dopamine neurons might be important for low firing rates and fine-tuning basal dopamine levels, while switching excitability to type II during NMDAR and AMPAR activation may facilitate a transient increase in dopamine concentration, as type II neurons are more amenable to synchronization by mutual excitation.

Adaptation and shunting inhibition leads to pyramidal/interneuron gamma with sparse firing of pyramidal cells

Gamma oscillations are a prominent phenomenon related to a number of brain functions. Data show that individual pyramidal neurons can fire at rate below gamma with the population showing clear gamma oscillations and synchrony. In one kind of idealized model of such weak gamma, pyramidal neurons fire in clusters. Here we provide a theory for clustered gamma PING rhythms with strong inhibition and weaker excitation. Our simulations of biophysical models show that the adaptation of pyramidal neurons coupled with their low firing rate leads to cluster formation. A partially analytic study of a canonical model shows that the phase response curves with a near zero flat region, caused by the presence of the slow adaptive current, are the key to the formation of clusters. Furthermore we examine shunting inhibition and show that clusters become robust and generic.

A new regime for highly robust gamma oscillation with co-exist of accurate and weak synchronization in excitatory-inhibitory networks

A great number of biological experiments show that gamma oscillation occurs in many brain areas after the presentation of stimulus. The neural systems in these brain areas are highly heterogeneous. Specifically, the neurons and synapses in these neural systems are diversified; the external inputs and parameters of these neurons and synapses are heterogeneous. How the gamma oscillation generated in such highly heterogeneous networks remains a challenging problem. Aiming at this problem, a highly heterogeneous complex network model that takes account of many aspects of real neural circuits was constructed. The network model consists of excitatory neurons and fast spiking interneurons, has three types of synapses (GABAA, AMPA, and NMDA), and has highly heterogeneous external drive currents. We found a new regime for robust gamma oscillation, i.e. the oscillation in inhibitory neurons is rather accurate but the oscillation in excitatory neurons is weak, in such highly heterogeneous neural networks. We also found that the mechanism of the oscillation is a mixture of interneuron gamma (ING) and pyramidal-interneuron gamma (PING). We explained the mixture ING and PING mechanism in a consistent-way by a compound post-synaptic current, which has a slowly rising-excitatory stage and a sharp decreasing-inhibitory stage.

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.

Short-term episodes of hypoxia induce posthypoxic hyperexcitability and selective death of GABAergic hippocampal neurons

We have previously developed a rat hippocampal neuronal cell model for the registration of the preconditioning effect and posthypoxic hyperexcitability (Turovskaya et al., 2011). Repeated episodes of short-term hypoxia are reported to suppress the amplitude of Ca(2+) response to NMDA in majority of neurons, reflecting the effect of preconditioning in the culture. In addition, exposure to hypoxia causes posthypoxic hyperexcitability: this is characterized by the onset of spontaneous synchronous Ca(2+) transients in a population of neurons in a neural network during the period of reoxygenation after each hypoxic episode. The nature of this phenomenon is unknown, although it has been observed that there always exists a minority of neurons in which there is no effect of hypoxic preconditioning. In this small population of neurons, the amplitude of Ca(2+) response to NMDA is not suppressed, but rather increases after each episode of hypoxia. Here we report the type of these neurons and their role in the generation of posthypoxic hyperexcitability. We compared the effect of short-term hypoxia on the amplitude of the Ca(2+) response to NMDA and the Ca(2+) transient generation in two populations of neurons - inhibitory GABAergic and excitatory glutamatergic. We have demonstrated that the neurons in which the preconditioning effect was not observed are GABAergic. Moreover at the instant moment of the posthypoxic synchronous Ca(2+)-transient generation (during reoxygenation) there is a global increase of [Ca(2+)]i and subsequent apoptosis in some GABAergic neurons. Anti-inflammatory cytokine interleukin-10 prevents the development of posthypoxic hyperexcitability, inhibiting the spontaneous synchronous Ca(2+) transients. At the same time, interleukin-10 protects GABAergic neurons from death, by restoring the effect of hypoxic preconditioning in them. Activation of one of the signaling pathways initiated by interleukin-10 appears to be necessary for the development of hypoxic preconditioning in GABAergic neurons. Overall our results indicate that short-term episodes of hypoxia can damage GABAergic neurons and weaken the inhibitory action of GABAergic neurons in a neural network. Activation of PI3K-dependent survival signaling pathways in neurons of this type is a possible strategy to protect these cells against hypoxia.

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.

When Long-Range Zero-Lag Synchronization is Feasible in Cortical Networks

Many studies have reported long-range synchronization of neuronal activity between brain areas, in particular in the beta and gamma bands with frequencies in the range of 14–30 and 40–80 Hz, respectively. Several studies have reported synchrony with zero phase lag, which is remarkable considering the synaptic and conduction delays inherent in the connections between distant brain areas. This result has led to many speculations about the possible functional role of zero-lag synchrony, such as for neuronal communication, attention, memory, and feature binding. However, recent studies using recordings of single-unit activity and local field potentials report that neuronal synchronization may occur with non-zero phase lags. This raises the questions whether zero-lag synchrony can occur in the brain and, if so, under which conditions. We used analytical methods and computer simulations to investigate which connectivity between neuronal populations allows or prohibits zero-lag synchrony. We did so for a model where two oscillators interact via a relay oscillator. Analytical results and computer simulations were obtained for both type I Mirollo–Strogatz neurons and type II Hodgkin–Huxley neurons. We have investigated the dynamics of the model for various types of synaptic coupling and importantly considered the potential impact of Spike-Timing Dependent Plasticity (STDP) and its learning window. We confirm previous results that zero-lag synchrony can be achieved in this configuration. This is much easier to achieve with Hodgkin–Huxley neurons, which have a biphasic phase response curve, than for type I neurons. STDP facilitates zero-lag synchrony as it adjusts the synaptic strengths such that zero-lag synchrony is feasible for a much larger range of parameters than without STDP.

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

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies.

Synchronization of neuronal spiking in the brain is related to cognitive functions, such as perception, attention, and memory. It is therefore important to determine which properties of neurons influence their collective behavior in a network and to understand how. A prominent feature of many cortical neurons is spike frequency adaptation, which is caused by slow transmembrane currents. We investigated how these adaptation currents affect the synchronization tendency of coupled model neurons. Using the efficient adaptive exponential integrate-and-fire (aEIF) model and a biophysically detailed neuron model for validation, we found that increased adaptation currents promote synchronization of coupled excitatory neurons at lower spike frequencies, as long as the conduction delays between the neurons are negligible. Inhibitory neurons on the other hand synchronize in presence of conduction delays, with or without adaptation currents. Our results emphasize the utility of the aEIF model for computational studies of neuronal network dynamics. We conclude that adaptation currents provide a mechanism to generate low frequency oscillations in local populations of excitatory neurons, while faster rhythms seem to be caused by inhibition rather than excitation.

Sparse gamma rhythms arising through clustering in adapting neuronal networks

Gamma rhythms (30-100 Hz) are an extensively studied synchronous brain state responsible for a number of sensory, memory, and motor processes. Experimental evidence suggests that fast-spiking interneurons are responsible for carrying the high frequency components of the rhythm, while regular-spiking pyramidal neurons fire sparsely. We propose that a combination of spike frequency adaptation and global inhibition may be responsible for this behavior. Excitatory neurons form several clusters that fire every few cycles of the fast oscillation. This is first shown in a detailed biophysical network model and then analyzed thoroughly in an idealized model. We exploit the fact that the timescale of adaptation is much slower than that of the other variables. Singular perturbation theory is used to derive an approximate periodic solution for a single spiking unit. This is then used to predict the relationship between the number of clusters arising spontaneously in the network as it relates to the adaptation time constant. We compare this to a complementary analysis that employs a weak coupling assumption to predict the first Fourier mode to destabilize from the incoherent state of an associated phase model as the external noise is reduced. Both approaches predict the same scaling of cluster number with respect to the adaptation time constant, which is corroborated in numerical simulations of the full system. Thus, we develop several testable predictions regarding the formation and characteristics of gamma rhythms with sparsely firing excitatory neurons.

Network with shunting synapses as a non-linear frequency modulator

The role of 'noisy' excitation in synchronizing interneuron networks with shunting synapses was studied. The excitatory input was simulated as a Poisson pattern of presynaptic conductance with varying frequencies and amplitudes. We find that higher excitation frequencies induce stronger synchronisation of the network. Within the range of 1-10000 Hz, only frequencies between 20 Hz and 200 Hz affected network synchronisation. No detectable network synchronisation was found at excitation frequencies below 20 Hz, and the network's synchronisation was either almost independent of the external input or falling down to zero when the input frequency was greater than 200 Hz. Thus the network transformed the input signals with frequencies above 20 Hz into output signals with the network's synchronisation frequency. The network's synchronisation frequency in our model ranged from 20 to 68 Hz depending on the frequency of the excitatory input. We conclude that a network of interconnected interneurons is capable of converting an asynchronous excitatory input into a synchronous inhibitory output as a frequency amplifier with the amplification coefficient dependent on the number of converging excitatory inputs. Another important result of our work revealed that the external frequency may affect, in opposite ways, the frequency of the network with shunting synapses depending on the excitatory synaptic conductance and the magnitude of leak conductance.

Control of neural synchrony using channelrhodopsin-2: a computational study

In this paper, we present an optical stimulation based approach to induce 1:1 in-phase synchrony in a network of coupled interneurons wherein each interneuron expresses the light sensitive protein channelrhodopsin-2 (ChR2). We begin with a transition rate model for the channel kinetics of ChR2 in response to light stimulation. We then define "functional optical time response curve (fOTRC)" as a measure of the response of a periodically firing interneuron (transfected with ChR2 ion channel) to a periodic light pulse stimulation. We specifically consider the case of unidirectionally coupled (UCI) network and propose an open loop control architecture that uses light as an actuation signal to induce 1:1 in-phase synchrony in the UCI network. Using general properties of the spike time response curves (STRCs) for Type-1 neuron model (Ermentrout, Neural Comput 8:979-1001, 1996) and fOTRC, we estimate the (open loop) optimal actuation signal parameters required to induce 1:1 in-phase synchrony. We then propose a closed loop controller architecture and a controller algorithm to robustly sustain stable 1:1 in-phase synchrony in the presence of unknown deviations in the network parameters. Finally, we test the performance of this closed-loop controller in a network of mutually coupled (MCI) interneurons.

Phase-response curves and synchronized neural networks

We review the principal assumptions underlying the application of phase-response curves (PRCs) to synchronization in neuronal networks. The PRC measures how much a given synaptic input perturbs spike timing in a neural oscillator. Among other applications, PRCs make explicit predictions about whether a given network of interconnected neurons will synchronize, as is often observed in cortical structures. Regarding the assumptions of the PRC theory, we conclude: (i) The assumption of noise-tolerant cellular oscillations at or near the network frequency holds in some but not all cases. (ii) Reduced models for PRC-based analysis can be formally related to more realistic models. (iii) Spike-rate adaptation limits PRC-based analysis but does not invalidate it. (iv) The dependence of PRCs on synaptic location emphasizes the importance of improving methods of synaptic stimulation. (v) New methods can distinguish between oscillations that derive from mutual connections and those arising from common drive. (vi) It is helpful to assume linear summation of effects of synaptic inputs; experiments with trains of inputs call this assumption into question. (vii) Relatively subtle changes in network structure can invalidate PRC-based predictions. (viii) Heterogeneity in the preferred frequencies of component neurons does not invalidate PRC analysis, but can annihilate synchronous activity.

Synchronization of the small-world neuronal network with unreliable synapses

As is well known, synchronization phenomena are ubiquitous in neuronal systems. Recently a lot of work concerning the synchronization of the neuronal network has been accomplished. In these works, the synapses are usually considered reliable, but experimental results show that, in biological neuronal networks, synapses are usually unreliable. In our previous work, we have studied the synchronization of the neuronal network with unreliable synapses; however, we have not paid attention to the effect of topology on the synchronization of the neuronal network. Several recent studies have found that biological neuronal networks have typical properties of small-world networks, characterized by a short path length and high clustering coefficient. In this work, mainly based on the small-world neuronal network (SWNN) with inhibitory neurons, we study the effect of network topology on the synchronization of the neuronal network with unreliable synapses. Together with the network topology, the effects of the GABAergic reversal potential, time delay and noise are also considered. Interestingly, we found a counter-intuitive phenomenon for the SWNN with specific shortcut adding probability, that is, the less reliable the synapses, the better the synchronization performance of the SWNN. We also consider the effects of both local noise and global noise in this work. It is shown that these two different types of noise have distinct effects on the synchronization: one is negative and the other is positive.

The response of a classical Hodgkin-Huxley neuron to an inhibitory input pulse

A population of uncoupled neurons can often be brought close to synchrony by a single strong inhibitory input pulse affecting all neurons equally. This mechanism is thought to underlie some brain rhythms, in particular gamma frequency (30–80 Hz) oscillations in the hippocampus and neocortex. Here we show that synchronization by an inhibitory input pulse often fails for populations of classical Hodgkin–Huxley neurons. Our reasoning suggests that in general, synchronization by inhibitory input pulses can fail when the transition of the target neurons from rest to spiking involves a Hopf bifurcation, especially when inhibition is shunting, not hyperpolarizing. Surprisingly, synchronization is more likely to fail when the inhibitory pulse is stronger or longer-lasting. These findings have potential implications for the question which neurons participate in brain rhythms, in particular in gamma oscillations.

Synchrony with shunting inhibition in a feedforward inhibitory network

Recent experiments have shown that GABA(A) receptor mediated inhibition in adult hippocampus is shunting rather than hyperpolarizing. Simulation studies of realistic interneuron networks with strong shunting inhibition have been demonstrated to exhibit robust gamma band (20-80 Hz) synchrony in the presence of heterogeneity in the intrinsic firing rates of individual neurons in the network. In order to begin to understand how shunting can contribute to network synchrony in the presence of heterogeneity, we develop a general theoretical framework using spike time response curves (STRC's) to study patterns of synchrony in a simple network of two unidirectionally coupled interneurons (UCI network) interacting through a shunting synapse in the presence of heterogeneity. We derive an approximate discrete map to analyze the dynamics of synchronous states in the UCI network by taking into account the nonlinear contributions of the higher order STRC terms. We show how the approximate discrete map can be used to successfully predict the domain of synchronous 1:1 phase locked state in the UCI network. The discrete map also allows us to determine the conditions under which the two interneurons can exhibit in-phase synchrony. We conclude by demonstrating how the information from the study of the discrete map for the dynamics of the UCI network can give us valuable insight into the degree of synchrony in a larger feed-forward network of heterogeneous interneurons.

Functional phase response curves: a method for understanding synchronization of adapting neurons

Phase response curves (PRCs) for a single neuron are often used to predict the synchrony of mutually coupled neurons. Previous theoretical work on pulse-coupled oscillators used single-pulse perturbations. We propose an alternate method in which functional PRCs (fPRCs) are generated using a train of pulses applied at a fixed delay after each spike, with the PRC measured when the phasic relationship between the stimulus and the subsequent spike in the neuron has converged. The essential information is the dependence of the recovery time from pulse onset until the next spike as a function of the delay between the previous spike and the onset of the applied pulse. Experimental fPRCs in Aplysia pacemaker neurons were different from single-pulse PRCs, principally due to adaptation. In the biological neuron, convergence to the fully adapted recovery interval was slower at some phases than that at others because the change in the effective intrinsic period due to adaptation changes the effective phase resetting in a way that opposes and slows the effects of adaptation. The fPRCs for two isolated adapting model neurons were used to predict the existence and stability of 1:1 phase-locked network activity when the two neurons were coupled. A stability criterion was derived by linearizing a coupled map based on the fPRC and the existence and stability criteria were successfully tested in two-simulated-neuron networks with reciprocal inhibition or excitation. The fPRC is the first PRC-based tool that can account for adaptation in analyzing networks of neural oscillators.

Phase-resetting curves determine synchronization, phase locking, and clustering in networks of neural oscillators

Networks of model neurons were constructed and their activity was predicted using an iterated map based solely on the phase-resetting curves (PRCs). The predictions were quite accurate provided that the resetting to simultaneous inputs was calculated using the sum of the simultaneously active conductances, obviating the need for weak coupling assumptions. Fully synchronous activity was observed only when the slope of the PRC at a phase of zero, corresponding to spike initiation, was positive. A novel stability criterion was developed and tested for all-to-all networks of identical, identically connected neurons. When the PRC generated using N-1 simultaneously active inputs becomes too steep, the fully synchronous mode loses stability in a network of N model neurons. Therefore, the stability of synchrony can be lost by increasing the slope of this PRC either by increasing the network size or the strength of the individual synapses. Existence and stability criteria were also developed and tested for the splay mode in which neurons fire sequentially. Finally, N/M synchronous subclusters of M neurons were predicted using the intersection of parameters that supported both between-cluster splay and within-cluster synchrony. Surprisingly, the splay mode between clusters could enforce synchrony on subclusters that were incapable of synchronizing themselves. These results can be used to gain insights into the activity of networks of biological neurons whose PRCs can be measured.

Self-organization of feed-forward structure and entrainment in excitatory neural networks with spike-timing-dependent plasticity

Spike-timing dependent plasticity (STDP) is an organizing principle of biological neural networks. While synchronous firing of neurons is considered to be an important functional block in the brain, how STDP shapes neural networks possibly toward synchrony is not entirely clear. We examine relations between STDP and synchronous firing in spontaneously firing neural populations. Using coupled heterogeneous phase oscillators placed on initial networks, we show numerically that STDP prunes some synapses and promotes formation of a feedforward network. Eventually a pacemaker, which is the neuron with the fastest inherent frequency in our numerical simulations, emerges at the root of the feedforward network. In each oscillatory cycle, a packet of neural activity is propagated from the pacemaker to downstream neurons along layers of the feedforward network. This event occurs above a clear-cut threshold value of the initial synaptic weight. Below the threshold, neurons are self-organized into separate clusters each of which is a feedforward network.

Asynchronous response of coupled pacemaker neurons

We study a network model of two conductance-based pacemaker neurons of differing natural frequency, coupled with either mutual excitation or inhibition, and receiving shared random inhibitory synaptic input. The networks may phase lock spike to spike for strong mutual coupling. But the shared input can desynchronize the locked spike pairs by selectively eliminating the lagging spike or modulating its timing with respect to the leading spike depending on their separation time window. Such loss of synchrony is also found in a large network of sparsely coupled heterogeneous spiking neurons receiving shared input.