QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

HOME  |   SEARCH  |   ARCHIVE  |   SUBSCRIBE  |   CONTACT  |   HELP

The Journal of Neuroscience, July 16, 2003, 23(15):6280-6294

This Article Abstract Full Text (PDF) A correction has been published Submit an eLetter Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Web of Science (40) Citing Articles via Google Scholar Google Scholar Articles by Pfeuty, B. Articles by Hansel, D. Search for Related Content PubMed PubMed Citation Articles by Pfeuty, B. Articles by Hansel, D. Pubmed/NCBI databases Compound via MeSH Substance via MeSH Hazardous Substances DB POTASSIUM SODIUM

 Previous Article  |  Next Article 

Electrical Synapses and Synchrony: The Role of Intrinsic Currents

Benjamin Pfeuty,¹ Germán Mato,² David Golomb,^3,4 and David Hansel^1,5

¹Laboratoire de Neurophysique et Physiologie du Système Moteur, Centre National de la Recherche Scientifique–Unité Mixte de Recherche 8119, Université René Descartes, 75270 Paris Cedex 06, France, ²Comisión Nacional de Energia Atómica and Consejo Nacional de Investigaciones Cientificas y Técnicas, Centro Atómico Bariloche and Instituto Balsiero, Universidad Nacional de Cordoba, 8400 San Carlos de Bariloche, Argentina, ³Department of Physiology and Zlotowski Center for Neuroscience, Faculty of Health Sciences, Ben Gurion University of the Negev, Be'er-Sheva 84105, Israel, ⁴Mathematical Biosciences Institute, Ohio State University, Columbus, Ohio 43210, and ⁵Interdisciplinary Center for Neural Computation, The Hebrew University, Jerusalem 91904, Israel

	Abstract

Top Abstract Introduction Materials and Methods Results Discussion Appendix References

 
Electrical synapses are ubiquitous in the mammalian CNS. Particularly in the neocortex, electrical synapses have been shown to connect low-threshold spiking (LTS) as well as fast spiking (FS) interneurons. Experiments have highlighted the roles of electrical synapses in the dynamics of neuronal networks. Here we investigate theoretically how intrinsic cell properties affect the synchronization of neurons interacting by electrical synapses. Numerical simulations of a network of conductance-based neurons randomly connected with electrical synapses show that potassium currents promote synchrony, whereas the persistent sodium current impedes it. Furthermore, synchrony varies with the firing rate in qualitatively different ways depending on the intrinsic currents. We also study analytically a network of quadratic integrate-and-fire neurons. We relate the stability of the asynchronous state of this network to the phase-response function (PRF), which characterizes the effect of small perturbations on the firing timing of the neurons. In particular, we show that the greater the skew of the PRF toward the first half of the period, the more stable the asynchronous state. Combining our simulations with our analytical results, we establish general rules to predict the dynamic state of large networks of neurons coupled with electrical synapses. Our work provides a natural explanation for surprising experimental observations that blocking electrical synapses may increase the synchrony of neuronal activity. It also suggests different synchronization properties for LTS and FS cells. Finally, we propose to further test our predictions in experiments using dynamic clamp techniques.

Key words: electrical synapses; conductance-based model; synchrony; neuronal network model; intrinsic currents; neocortex

	Introduction

Top Abstract Introduction Materials and Methods Results Discussion Appendix References

 
Electrical synapses are sites at which gap-junctions bridge the membranes of two neurons. They have long been known to exist in invertebrates (Watanabe, 1958; Furshpan and Potter, 1959), but only recently has evidence of their ubiquity been unequivocally found in the mammalian brain. Electrical synapses are present in the inferior olive (Llinas and Yarom, 1986), the hippocampus (Draghun et al., 1998; Venance et al., 2000), the cerebellum (Mann-Metzer and Yarom, 1999), the locus coereleus (Christie et al., 1989; Alvarez et al., 2002), the striatum (Kita et al., 1990), the neocortex (Galarreta and Hestrin, 1999; Gibson et al., 1999), the reticular thalamic nucleus (Landisman et al., 2002), and between motoneurons (Kiehn and Tresch, 2002).

Experiments have revealed that electrical synapses are involved in synchronizing neural activity (Draghun et al., 1998; Mann-Metzer and Yarom, 1999; Beierlein et al., 2000; Perez-Velazquez and Carlen, 2000; Tamas et al., 2000; Deans et al., 2001; Hormuzdi et al., 2001). In contrast to the results of these studies, it has been reported recently that inspiratory motoneurons may display more strongly synchronized activity in presence of carbenoxolone (CBX), a blocker of electrical synapses, than in the control situation (Bou-Flores and Berger, 2001). Therefore, in this case, electrical synapses desynchronize neural activity.

The dynamics of networks of neurons interacting via chemical synapses have been studied extensively (Golomb et al., 2001). However, only a few theoretical studies have addressed the dynamics of networks in which neurons are coupled by electrical synapses. Models for pattern generation in the lobster pyloric system have been investigated (Kepler et al., 1990; Abbott et al., 1991; Meunier, 1992). Stable antiphase locking and transitions between in-phase and antiphase locking were demonstrated in pairs of model neurons coupled by electrical synapses, and its impact on rhythmogenesis was investigated (Sherman and Rinzel, 1992; Cymbalyuk et al., 1994; Han et al., 1995). More recently, the discovery of electrical synapses in the CNS of mammals has motivated simulations of conductance-based models (Traub et al., 2001) and analytical studies in the framework of leaky integrate-and-fire (LIF) models (Chow and Kopell, 2000; Lewis and Rinzel, 2003).

Although it is clear that synchronization properties of neurons depend on their intrinsic currents, the general principles that determine this dependence are still unknown in the case of neurons interacting via electrical synapses. In this paper, we combine numerical simulations of conductance-based network models and analytical investigations of quadratic integrate-and-fire (QIF) neurons (Hansel and Mato, 2003) to discover such principles. In particular, we show that potassium currents promote synchrony, whereas sodium currents impede it. We compare our results with existing experimental data. We suggest that our findings may account for those reported by Bou-Flores and Berger (2001). We also propose experiments to further test the predictions of our work. Part of this work has been published previously in abstract form (Pfeuty et al., 2002).

	Materials and Methods

Top Abstract Introduction Materials and Methods Results Discussion Appendix References

 
The conductance-based model. The membrane potential, V, of the neuron follows the equation:

(1)

where I_L = -g_L(V - V_L) is a leak current and _ionI_ion is the sum of all voltage-dependent ionic currents. The currents I_ext and I_noise are a constant external current and a Gaussian white noise with a zero mean and SD, , respectively.

Our neuron model incorporates an inactivating sodium current, I_Na = g_Nam ³ h(V - V_Na); a delayed rectifier potassium current, I_K = g_Kn ⁴(V - V_K); a slow potassium current, I_Ks = g_Kss ⁴(V - V_K) (Erisir et al., 1999); and a noninactivating (persistent) sodium current (French et al., 1990), I_NaP = g_NaPp(V - V_Na). The kinetics of the gating variable h, n, s are given by:

(2)

with x = h, n, s and _h(V) = 0.21e ^-(V+58)/20, _h(V) = 3/(1 + e ^-(V+28)/10), _n(V) = 0.03(V + 34)/(1 - e ^-(V+34)/10), _n(V) = 0.375e ^-(V+44)/80, _s(V) = 0.07(44 + V)/(1 - e ^-(V+44)/4.6), _s(V) = 0.008e ^-(V+44)/68. The activation functions m and p are given by: m(V) = _m(V)/(_m(V) + _m(V)), where _m(V) = 0.1(V + 35)/(1 - e ^-(V+35)/10), M(V) = 4e ^-(V+60)/18, and p(V) = 1/(1 + e ^-(V+50)/6).

Throughout this work, the parameters g_Na = 35 mS/cm ², V_Na = 55 mV, V_K = -90 mV, g_L = 0.1 mS/cm ², V_L = -65 mV, and C = 1 µF/cm ² are kept constant, and we study the ways in which the network dynamics depend on the conductances g_K, g_Ks, and g_NaP.

The network architecture. The network consists of N neurons randomly connected by bidirectional electrical couplings. The average number of connections per neuron is denoted by M. We assume that all existing synapses have the same conductance, g, and we define a connectivity matrix by J_ij = g if neuron i and j, (i, j = 1..., N), are connected and J_ij = 0 otherwise. The connection between neuron i and j adds a contribution, -g(V_i - V_j), to the total synaptic current received by neuron i, where V_i and V_j are the membrane potential of neurons i and j, respectively. Therefore, to include the effect of the electrical synapses in the dynamics of the network, an additional current I_i^ES = -_jJ_ij(V_i - V_j) is added to the right side of the differential equation satisfied by the membrane potential of neuron i (see Eqs. 1, 9). In our simulations, the network size is N = 1600, the average connectivity is M = 10, and the conductance is the same for all of the synapses: g = 0.005 mS/cm ².

Definition of synchrony and measure of the synchrony level in numerical simulations. By synchrony of the activity of a pair of neurons, we mean the tendency of these neurons to fire spikes at the same time. In this sense, a pair of neurons will be said to fire in synchrony if the cross-correlation of their spike trains displays a central peak that is above chance. The degree of synchrony can be characterized by the amplitude of this peak. According to this definition, if two neurons fire action potentials periodically and in antiphase, they can be considered to fire asynchronously.

In a large network of neurons, the activity is asynchronous if at any time the number of action potentials fired in the network is the same up to some random fluctuations. This implies that in a large network, the asynchronous state is unique. When this state is unstable, the network activity is necessarily synchronous. In this case, neurons tend to fire preferentially in some windows of time.

Note that a stable asynchronous state can coexist with a stable synchronous state if the network dynamics display multistability. In this work we focus on the conditions under which asynchronous states in large networks and antiphase locking states in pairs of identical neurons become unstable because of electrical synaptic interactions.

One can characterize the degree of synchrony in a population of N neurons by measuring the temporal fluctuations of macroscopic observables, as for example, the membrane potential averaged over the population (Hansel and Sompolinsky, 1992; Golomb and Rinzel, 1994). The quantity

(3)

is evaluated over time, and the variance of its temporal fluctuations is computed, where denotes time-averaging. After normalization of _V to the average over the population of the single-cell membrane potentials, , one defines (N):

(4)

which varies between 0 and 1. The central limit theorem implies that in the limit n , (N) behaves as:

(5)

where a > 0 is a constant and O(1/N) means a term of order 1/N. In particular, (N) = 1 if the activity of the network is fully synchronized (i.e., V_i(t) = V(t) for all i), and (N) = O(1/N) if the state of the network activity is asynchronous. In the asynchronous state, () = 0. More generally, the larger (), the more synchronized the population. Note that this measure of synchrony is sensitive not only to the correlations in the spike timing but also to the correlations in the time course of the membrane potentials in the subthreshold range.

Measure of the irregularity of spike trains. The irregularity of the spike trains is characterized by the coefficient of variability (CV) (i.e., the ratio between the SD of the interspike intervals and their mean values averaged over the entire network). For periodic spike trains, CV = 0, whereas CV = 1 for spike trains randomly distributed with Poisson statistics.

Quadratic integrate-and-fire model. It can be shown that near the onset of periodic firing, the detailed subthreshold dynamics of a large class of neurons can be reduced to a simple one-dimensional model in which the dynamical variable, v_red, evolves in time (Ermentrout, 1996; Hansel and Mato, 2003):

(6)

The constants, C_red, A, V^*, I_c^red, the reduced external current, I_ext ^red, and the reduced noise, I_{noise red}, can be computed as functions of the parameters of the full model as shown by Ermentrout (1996) and Hansel and Mato (2003). For I_red > I_c^red, the solution of Equation 6 diverges in finite time. This corresponds to the firing of an action potential. If one supplements Equation 6 with the condition that the variable v_red is reset to V_r < V^* immediately after v_red reaches some threshold, V_T, from below, this yields a reduced model that accurately describes the dynamics of the neuron in the limit I_ext ^red I_c^red. In particular, the current–frequency relationship (I–F curve) of this model and the neuron behave similarly in this limit, for any value of the parameters V_r and V_T. When I_ext ^red-I_c^red is not small, the reduced model is no longer an exact description of the neuronal dynamics. However, one can fit the parameters so that the reduced model provides a good approximation of the I–F curve of the neuron. The difference, V_T - V_r, will be called the reset depth.

It is convenient to rewrite the subthreshold dynamics of the reduced model in terms of the dimensionless variables _red and ^red defined by:

(7)

(8)

where ₀ has the dimension of a time. This yields:

(9)

The variable _red is reset to _r = A(V_r - V^*)₀ /C_red whenever it reaches the threshold value: _T = A(V_T - V^*)₀ /C_red. Note that _r < 0. A similar model has been used to study networks of excitatory and inhibitory neurons (Latham et al., 2000; Hansel and Mato, 2003).

The subthreshold dynamics (Eq. 9) need to be supplemented with a model for the suprathreshold part of the membrane potential time course. Assuming that the width of the action potentials is much smaller than the interspikes, we represent their time course by a -function at each time a spike is fired. Therefore, the (dimensionless) reduced membrane potential of the neuron can be written:

(10)

where measures the integral over time of the suprathreshold part of an action potential.

The model defined by Equations 9 and 10 will be called the QIF model. To simplify notations of the reduced model, we will drop the index ("red") and the tildes.

Phase reduction in the weak coupling limit. Let us consider a neuron firing action potentials periodically with an interspike T. A small and instantaneous perturbation applied at time t after a spike induces a small change in the timing of the subsequent spikes. This change, which depends on t, or equivalently on = 2t/T, can be characterized by a function, Z(), which measures the delay or the advance induced in the firing times after the perturbation. A positive value of Z() indicates that the perturbation advances the subsequent spikes. A negative value of Z() indicates a delay. The effect of noninstantaneous weak perturbations (such as spikelets) can be estimated by convolving the phase-response function (PRF), Z, with the perturbation. It can be shown that the dynamic behavior of a network of weakly interacting neurons can be completely described in terms of the response function in the framework of the phase reduction approach (Kuramoto, 1984; Ermentrout and Kopell, 1986; Hansel et al., 1993, 1995; Golomb and Hansel, 2000; Neltner et al., 2000; Lewis and Rinzel, 2003). In this approach, a phase variable is associated with each neuron in the network. This variable, _i(i = 1,..., N), measures the time elapsed since the last action potential fired by neuron i. For a network of identical neurons coupled with electrical synapses, one can show (see Appendix) that the phase variables follow a set of n first-order coupled differential equations:

(11)

(_i - _j) the phase coupling between neurons i and j given by:

(12)

J_ij is the connectivity matrix, and _i(t) is a white noise with zero mean and variance

with , the SD of the noise of the nonreduced dynamics.

Numerical integration. In the simulations of the conductance-based model, the differential equations were integrated using the second-order Runge-Kutta scheme with fixed-time step: _t = 0.01 msec. Averaged quantities (firing rate, CV, ) were computed over a time period of 1 sec after discarding a transient of 500 msec.

	Results

Top Abstract Introduction Materials and Methods Results Discussion Appendix References

 
Single neuron and coupling properties of the conductance-based model
Firing properties of single neurons
In the absence of persistent sodium and slow potassium currents (the "control model"), our conductance-based model neuron fires tonically for large enough external currents, I_ext > I_c = 0.16 µA/cm². The frequency of the discharge in response to a step of current is plotted as a function of the step amplitude (I–F curve) in Figure 1A. For I_ext larger than but close to I_c, the minimum value of the external current required for the neuron to fire, the firing frequency can be arbitrarily small, because action potentials appear at I_c through a saddle-node bifurcation [for a definition of a saddle-node bifurcation, see Strogatz (1994) and Rinzel and Ermentrout (1998)]. Decreasing g_K does not significantly modify the resting potential and the rheobase of the neuron, but it increases the gain of the I–F curve (Fig. 1B). The slow potassium (I_Ks) and persistent sodium (I_NaP) currents modify the onset of firing, because these currents are activated at rest. As should be expected, the persistent sodium current increases the excitability of the neuron and reduces its rheobase, I_c (Fig. 1C). For g_NaP > 0.08 mS/cm², the neuron is spontaneously active. As shown in Figure 1D, the main effects of adding I_Ks on the I–F curve are a translation to the right, a linearization, and, if g_Ks is large enough, a suppression of the ability of the neuron to fire at low rates (because for large enough g_Ks, a subcritical Hopf bifurcation occurs at the onset of firing [for a definition of a subcritical Hopf bifurcation, see Strogatz (1994) and Rinzel and Ermentrout (1998)].

View larger version (26K): [in this window] [in a new window]   Figure 1. I–F curves of the conductance-based model neuron. Each panel corresponds to a different set of intrinsic conductance values. The firing rate was computed from the steady-state response of the neuron to steps of constant external current of different amplitudes. A, Control case: gK = 9 mS/cm 2, gKs = gNaP = 0. B, gK = 2.5 mS/cm 2, gKs = gNaP = 0. Note that the scale of the x-axis is one-half the scale used in A. The main effect on the I–F curve of the reduction of gK is multiplicative (change in the gain). C, gK = 9 mS/cm 2, gNaP = 0.2 mS/cm 2, gKs = 0. The main effect on the I–F curve of the persistent sodium is subtractive (change in the rheobase). D, gK = 2.5 mS/cm 2, gKs = 0.2 mS/cm 2, gNaP = 0. The slow potassium current reduces the gain of the neuron and prevents low-frequency firing. It also linearizes the I–Fcurve (compare with B).

 

The response of the neuron to a step of external current, whose value is adjusted to give a firing frequency of 50 Hz, is plotted in Figure 2 for different combinations of the intrinsic currents. The action potentials are shown with a higher temporal resolution in Figure 3C. In the control case (Fig. 2A; and Fig. 3C, solid line), action potentials are narrow and are followed by a strong afterhyperpolarization (AHP) that brings the membrane potential 15 mV below its resting value. When g_K is decreased (Fig. 2B), the resting potential of the neuron does not change, but the AHP is reduced and the membrane potential always remains above rest. This conductance also controls the width of the spike and the depth of the AHP [Fig. 3C, compare the solid line (g_K = 9 mS/cm²) with the dotted line (g_K = 2.5 mS/cm²)].

View larger version (27K): [in this window] [in a new window]   Figure 2. The response of the neuron to a step of current. A step of current of 200 msec was injected into the neuron. No noise was included. The membrane potential before and after the current injection is indicated to the left of each panel. In all cases, the external current during the step was chosen to obtain a discharge rate of 50 Hz. A, Control case. External current, Iext, varies from 0 to 1.10 µA/cm 2. Resting potential is -63 mV. B, gK = 2.5 mS/cm 2, gKs = gNaP = 0. Iext varies from 0 to 0.48 µA/cm 2. Note that the concavity of the membrane potential time course is always upward in contrast to A. C, gK = 9 mS/cm 2, gNaP = 0.2 mS/cm 2, gKs = 0. The neuron is hyperpolarized by the injection of a negative current, Iext = -1.55 µA/cm 2, before the current step, to prevent spontaneous firing. The step amplitude is 1.0 µA/cm 2. D, gK = 2.5 mS/cm 2, gKs = 0.2 mS/cm 2, gNaP = 0. The external current varies from 0 to 4.88 µA/cm 2. The resting membrane potential is more hyperpolarized than in B, because of the slow potassium current that is activated at rest. The (shifted) time course of the membrane potential during the discharge is very similar to the trace in B. Note the transient hyperpolarization right after the switching off of the current, which is caused by the slow relaxation of IKs back to its resting value.

 

View larger version (24K): [in this window] [in a new window]   Figure 3. The spikelets induced for different combinations of ionic currents. Solid line, gK = 9 mS/cm 2, gKs = gNaP = 0 (control case). Dotted line, gK = 2.5 mS/cm 2, gKs = gNaP = 0. Dashed line, gK = 9 mS/cm 2, gKs = 0, gNaP = 0.2 mS/cm 2. Dashed-dotted line, gK = 2.5 mS/cm 2, gKs = 0.2 mS/cm 2, gNaP = 0. A, An action potential generated by a short and strong current pulse. The duration and the amplitude A of the pulse are the same for all four traces: = 1 msec and A = 50 µA/cm 2. The neuron was hyperpolarized to prevent firing when persistent sodium was present and to obtain the same initial membrane potential for all four parameter sets. The solid line and the dashed line overlap. B, Spikelet induced by a presynaptic neuron firing an action potential as in A. The synaptic conductance is g = 0.005 mS/cm 2, which corresponds to a coupling coefficient (Amitai et al., 2002) CC5%. The postsynaptic neuron has the same intrinsic properties as the presynaptic neuron. C, A constant external current is injected to make the neuron fire at 50 Hz. The external current is Iext = 1.10 µA/cm 2 for the solid line, Iext = 0.48 µA/cm 2 for the dotted line, Iext = 4.88 µA/cm 2 for the dashed-dotted line, and Iext = -0.55 µA/cm 2 for the dashed line. The subthreshold voltage of the neuron when the persistent sodium is added (dashed line) is below the control (solid line), although the persistent sodium current accelerates the depolarization of the neuron, because here we compare the firing patterns for the same firing rate. The external current is therefore smaller when the persistent sodium is present than in the control case. D, Spikelets induced by a presynaptic neuron firing tonically, as in C. The postsynaptic and the presynaptic neurons have the same intrinsic properties.

 

The persistent sodium and the slow potassium currents significantly modify the resting potential and the external current required to adjust the firing rate. In contrast, they only slightly affect the depth of the repolarization, measured from threshold, after the action potential (Fig. 2C,D). The size and width of the action potentials remain almost unchanged when the conductances of these currents are varied [Fig. 3C, compare the dotted line with the dashed-dotted line (effect of g_Ks) and the solid line with the dashed line (effect of g_NaP)].

Properties of the spikelets
When a neuron fires an action potential, it induces a depolarization, called a spikelet, in the neurons that are connected to it. The spikelets generated by a neuron firing one isolated spike in response to a very brief but strong pulse of current (Fig. 3A) or firing tonically at 50 Hz (Fig. 3C), are shown in Figure 3, B and D, respectively, for different parameters of the ionic currents.

In all four cases displayed, the width of the spikelets is much larger than the width of the presynaptic spikes that generate them. This is because the spikelets are a filtered version of the presynaptic membrane potential profile.

The size of the spikelet is practically unaffected by the presence of the persistent sodium (Fig. 3B,D, compare the solid line with the dashed line), because this current does not greatly affect the shape of the action potential. The only noticeable effect of the persistent sodium on the spikelet is to slow down its decay time course. This is because it reduces the input conductance of the postsynaptic neuron (by 50% at -76 mV).

The modulation of the spikelets by the delayed rectifier potassium current is primarily caused by changes in the presynaptic neuron voltage time course. When g_K is reduced, the action potential is broader, and this induces an increase in the amplitude and the width of the spikelets. Although the slow potassium current does not greatly affect the width of the action potential, it contributes to bringing the membrane of the presynaptic neuron transiently below its holding potential. This explains why in Figure 3B the spikelet is narrower and displays a faster and deeper repolarization in the presence of I_Ks.

The spikelets generated by tonically firing neurons at low firing rates are similar to those in Figure 3B (data not shown). For high-enough firing rates, the slow potassium saturates, affecting the dynamics primarily as would a constant hyperpolarizing current. Therefore, its effect on the voltage traces and on the spikelets becomes less pronounced (Fig. 3, compare B and D).

Synchrony in the conductance-based network model
The goal of this section is to show how synchrony properties of our conductance-based network model depend on the intrinsic currents of the neurons. Further understanding of the effects described in this section will be provided in the next section by investigating analytically a simplified and more abstract model.

Synchrony is modified when intrinsic conductances are changed A pair of conductance-based neurons
The effect of intrinsic currents on neuronal network dynamics can first be demonstrated by studying the dynamics of a pair of neurons. Here, as an example, we show that a persistent sodium current tends to promote antiphase locking, whereas a potassium current promotes in-phase locking. In Figure 4A, the traces of two identical neurons (g_K = 9 mS/cm²; g_Ks = g_NaP = 0) are plotted over a time interval of 100 sec. Both neurons receive a noisy input that makes them fire at 50 Hz. These traces suggest that the neurons fire action potentials in synchrony. This is confirmed by the cross-correlation of the traces computed over a longer time interval (100 sec), which displays a strong peak centered around t = 0 msec, as shown in Figure 4B. In contrast, when the persistent sodium conductance is large enough (g_NaP = 0.4 mS/cm²), neurons tend to fire in antiphase (Fig. 4C), and the peak of the cross-correlogram is shifted to 10 msec (Fig. 4D). If a slow potassium current is added (g_Ks = 0.15 mS/cm²) while keeping the persistent sodium conductance the same, antiphase locking becomes unstable and neurons tend to fire in-phase (Fig. 4E, F). Therefore, a slow potassium current is a promoter of synchronous activity.

View larger version (28K): [in this window] [in a new window]   Figure 4. Effect of persistent sodium and potassium currents on the synchronization of a pair of neurons: phase versus antiphase stability. Results from numerical simulations. The synaptic conductance is g = 0.005 mS/cm 2. The SD of the noise is = 0.3 mV/msec1/2. In all of the simulations, neurons were firing in antiphase at the beginning of the simulation (data not shown). Spike trains were computed with a time bin of 1 msec. The spike cross-correlograms (B, D, F) were calculated from spike trains recorded over 100 sec (after discarding a transient of 5 sec) and were normalized to the expected number of spikes during this time interval. Neurons fire at an average firing rate of 50 Hz. A, B, Traces of the voltage of the two neurons and the cross-correlograms of the spike trains for the control parameters and Iext = 1.08 µA/cm 2. Neurons tend to fire in-phase. C, D, Persistent sodium current stabilizes antiphase locking. Parameters are gK = 9 mS/cm 2, gNaP = 0.4 mS/cm 2, and gKs = 0. Iext = -1.38 µA/cm 2. E, F, Slow potassium conductance destabilizes antiphase locking. Parameters are gK = 9 mS/cm 2, gNaP = 0.4 mS/cm 2, and gKs = 0.15 mS/cm 2. Iext = 0.80 µA/cm 2. Neurons tend to fire in-phase.

 

Large networks of conductance-based neurons
In Figure 5 the synchrony measure, , is plotted against the conductances g_K, g_Ks, and g_NaP. Examples of raster plots of network activity are also displayed in this figure, for particular values of the intrinsic conductances (Fig. 5, insets). In each of the panels, the external input (average deviation and SD) is kept constant. Therefore, when the intrinsic conductances vary, the average firing rate of the neurons and the CV of their interspike intervals change. In particular, the firing rate decreases (respectively increases) when the potassium (respectively the persistent sodium) conductances increase (top right figures in each of the panels in Figure 5).

View larger version (29K): [in this window] [in a new window]   Figure 5. Effect of intrinsic currents on the synchronization of a large network. Results from numerical simulations. Details on the simulation and the network parameters are given in Materials and Methods. The synaptic conductance is g = 0.005 mS/cm 2. In all of the simulations, the SD of the noise was set at = 0.6 mV/msec 1/2. The synchrony level, (left), the average frequency, and the average CV of the neuronal discharges are plotted as a function of an intrinsic conductance that is varied. A, The conductance, gK, varies; gKs = gNaP = 0. The external current is Iext = 0.8 µA/cm 2. For gK > g*K, increases and saturates. Inset, Raster plot representing the spiking times of 100 neurons for gK = 9 mS/cm 2 (neurons are arrayed vertically, and the x-axis corresponds to time where the scale denotes 25 msec). Panels on the right indicate that the average frequency of the neurons decreases when gK increases, because the gain of the I–F curve decreases. The variability of the neuronal discharge varies slightly and nonmonotonically in the investigated range of gK. B, The conductance, gKs, varies; gK = 2.5 mS/cm 2, gNaP = 0. The external current is Iext = 2 µA/cm 2. Inset, Raster plots are for gKs = 0 and for gKs = 0.1 mS/cm 2. At large values of gKs, a slight decrease occurs, but the activity of the network remains strongly synchronized. Panels on the right indicate that the average firing rate of the neurons is strongly reduced when gKs varies, primarily because IKs has a strong subtractive effect on the I–F curve of the neurons (the rheobase increases with gKs). C, The conductance, gNaP, varies; gK = 9 mS/cm 2, gKs = 0. The external current is Iext = 0.8 µA/cm 2. The synchrony is reduced when gNaP increases. Inset, Raster plot for gNaP = 0.2 µA/cm 2. Panels on the right indicate that the average firing rate of the neurons increases with gNaP, primarily because of the subtractive effect of INaP on the I–F curve of the neurons (the rheobase is reduced when gNaP increases).

 

In Figure 5A, the effect of g_K on the synchrony level is shown (g_Ks = g_NaP = 0). When g_K is small enough (g_K < 4.5 mS/cm²), is very close to zero. A detailed study shows that in this region, is on the order of 1/N, where N is the number of neurons in the network, and vanishes in the limit of a very large network (data not shown): non-zero values of are attributable to finite size effects, and the network is asynchronous. For g_K of >4.5 mS/cm², increases rapidly up to 0.35. For example, in the raster plot shown in the right inset of Figure 5A (control situation, g_K = 9 mS/cm²), = 0.34. This corresponds to a state in which the neurons fire together within time windows of approximately one-third of the period.

Figure 5B shows the effect of the slow potassium current for g_NaP = 0 and g_K = 2.5 mS/cm². For g_Ks = 0, the neurons fire asynchronously. At g_Ks 0.06 mS/cm², starts to increase. For large g_Ks, saturates at a value, 0.55, that is substantially larger than in the control situation. This value of corresponds to a tight synchrony of the action potentials fired by the neurons (Fig. 5B inset on the right).

In Figure 5C, is plotted against g_NaP for g_K and g_Ks at their control values (g_K = 9 mS/cm², g_Ks = 0). Clearly, is a decreasing function of g_NaP. When g_NaP is large enough (g_NaP > 0.1 mS/cm²), the neurons fire asynchronously.

In the results presented in Figure 5, the firing rate of the neurons is not controlled. Therefore, the origin of the variation of with the intrinsic current conductances is not clear. It may be attributable to the fact that, in general, synchrony depends on the firing rate for given intrinsic currents. It might also be a result of the fact that dynamic properties of neurons with different intrinsic currents may be different even if they have the same firing rate.

Dependence of synchrony on intrinsic currents at fixed firing rates
To clarify further how potassium and sodium currents affect the synchrony of the network, we performed another set of numerical simulations in which the external input and the noise level were tuned to keep constant both the average firing rate of the neurons and the CV of their spike trains when the conductances of the intrinsic currents were changed. The results are shown in Figure 6 (three panels on the left). In all of the simulations, the average firing rate of the network was 50 Hz and the CV 0.12 (Fig. 6, panels on the right).

View larger version (24K): [in this window] [in a new window]   Figure 6. The synchrony level depends on the intrinsic currents when the frequency is kept constant. Results from numerical simulations. Details on the parameters are given in Materials and Methods. The synaptic conductance is g = 0.005 mS/cm 2. The external current, Iext, and the SD of the noise, , were changed to control the average firing rate and the average CV of the neurons. In all of the simulations, the average firing rate is 50 Hz ± 10%, and the CV 0.12. In each of the panels, the corresponding values of Iext and are plotted versus the conductance that is varied (panels to the right). A, The conductance, gK, varies; gKs = gNaP = 0. B, The conductance, gKs, varies; gK = 2.5 mS/cm 2, gNaP = 0. C, The conductance, gNaP, varies; gK = 9 mS/cm 2, gKs = 0.

 

The level of synchrony increases with the potassium conductances. For small-enough g_K or g_Ks, the activity of the network is weakly correlated. A detailed study of the dependence of with the network size reveals that for g_K (respectively g_Ks) smaller than g^*_K 4 mS/cm² (respectively g_Ks smaller than g^*_Ks 0.075 mS/cm²), is on the order of 1/N (data not shown). Therefore in this range of conductances, the network is in the asynchronous state. At g^*_K and g^*_Ks, a transition to a synchronous state occurs, and beyond these values the network settles into a synchronous state. This occurs although the level of noise has been increased.

In contrast to the potassium conductances, which promote synchrony, the persistent sodium current impedes it, as shown in Figure 6C. For large-enough g_NaP, synchrony is even destroyed. A sharp transition to the asynchronous state occurs for g_NaP 0.15 mS/cm².

Relying on the way intrinsic currents affect spikelet modulation (Fig. 3), one might naively expect that potassium currents that reduce the size of the spikelets should also reduce synchrony, and that persistent sodium current, which increases this size, should promote it. As a matter of fact, we found a trend that was exactly the opposite.

Firing rate affects synchrony differently when different intrinsic currents are involved
We also investigated in what ways the network dynamic state depends on the firing rate of the neurons. The firing rate was changed by varying the average external input, and was plotted as a function of the frequency of the neurons for several combinations of potassium and sodium conductances.

Figure 7 plots the synchrony measure, , versus the average firing rate of the neurons for four sets of conductance parameters and two levels of noise. The qualitative behavior of as a function of the firing rate depends on the intrinsic currents involved.

View larger version (25K): [in this window] [in a new window]   Figure 7. Synchrony varies with the firing rate in a way that depends on the intrinsic currents. Results from numerical simulations. Details on the simulation parameters are given in Materials and Methods. The synaptic conductance is g = 0.005 mS/cm 2. The average firing rate of the neurons is varied by changing the external current. In each of the panels, the results are displayed for two noise levels. Circles, = 0.6 mV/(msec)1/2. Squares, = 0.3 mV/(msec)1/2. A, gK = 9 mS/cm 2, gKs = gNaP = 0. The external current was varied between 0.15 and 2.7 µA/cm 2 (100 Hz). B, gK = 3.5 mS/cm 2, gKs = gNaP = 0. The external current, Iext, was varied between 0.11 and 1.4 µA/cm 2. C, gK = 9 mS/cm 2, gKs = 0, gNaP = 0.15 mS/cm 2. Iext was varied between -1.17 and 0.8 µA/cm 2. D, gK = 3.5 mS/cm 2, gKs = 0.15 mS/cm 2, gNaP = 0. Iext was varied between 2.3 and 7.2 µA/cm 2. The firing rate cannot be reduced to <20 Hz because of the presence of the slow potassium current (see also Fig. 1 D).

 

As expected, is always larger for the smaller level of noise ( = 0.3 mV/msec^1/2) (Fig. 7, circles) than for the larger one ( = 0.6 mV/msec^1/2) (Fig. 7, squares). In particular, for = 0.6 mV/msec^1/2, the noise is strong enough to prevent synchrony in the entire range of firing frequencies for g_K = 3.5 mS/cm² (Fig. 7B) and for g_NaP = 0.15 mS/cm² (Fig. 7C). This level of noise is sufficient to destroy synchrony as well in the control situation, but only at <20 Hz (Fig. 7A).

In the control case, increases monotonously with the firing rate (except for a small decrease followed by a slight increase between 50 and 100 Hz for = 0.3 mV/msec^1/2). In contrast, for smaller g_K (e.g., g_K = 3.5 mS/cm²), we find a nonmonotonous variation of synchrony: increases at low firing rates but decreases again for firing rates of >20 Hz. The network state becomes asynchronous for firing rates of >60 Hz (Fig. 7B, squares). If the slow potassium current is added (g_Ks = 0.15 mS/cm²), again varies monotonously with the firing rate (Fig. 7D). Finally, with enough persistent sodium current and not overly strong noise (g_K = 9 mS/cm², with g_NaP = 0.15 mS/cm², g_Ks = 0) (Fig. 7C, squares), varies nonmonotonously with the firing rate, as in the reduced g_K case (Fig. 7B).

Synchrony in networks of quadratic integrate-and-fire neurons depends on the phase response function
The PRF, which characterizes how neurons respond to small perturbations, is a key concept to understand the relationship between intrinsic properties of neurons and their collective dynamics (Kuramoto, 1984; Hansel et al., 1993; van Vreeswijk et al., 1994; Kopell and Ermentrout, 2002). This function depends on the excitability properties of the neurons and therefore is determined by the intrinsic currents involved in their dynamics. One expects that the differences in the dynamic behavior of our conductance-based network model for different sets of parameters reflect the changes in the excitability properties of the neurons. However, relying on numerical simulations alone, it would very difficult to establish general principles to relate the excitability properties and the neuronal PRF to the synchronization properties. Below, we consider a network of QIF neurons (Eqs. 9, 10) fully connected by electrical synapses. As we show below, the qualitative properties of the excitability of the QIF neurons crucially depend on the parameters (the threshold and the reset voltages and the external current). Thanks to its relative simplicity, this model can be investigated using analytical techniques. It reveals how, in the framework of this model, one can relate the stability of antiphase locking of a pair of neurons and the stability of the asynchronous state of a large network to the shape of the PRF of the neurons. Subsequently, we show that similar rules can be applied to conductance-based models, and that they provide a unified framework explaining the diversity of behavior found in our numerical simulations of the conductance-based model.

The membrane potential and the PRF of QIF neurons
In the Appendix, it is shown that between two successive spikes the subthreshold membrane potential of a QIF neuron reads:

(13)

where t measures the time elapsed since the first spike. The firing period, T, is determined by the condition: v(T) = V_T. Therefore:

(14)

In the weak coupling limit, the dynamics of the network are completely determined by the phase-response function, Z() (see Materials and Methods), which depends on three parameters, namely, the firing frequency of the neurons (or equivalently the external current I_ext), the threshold V_T, and the reset potential V_r. The response function of the QIF model can be derived analytically as shown in the Appendix. One finds:

(15)

Changing the ratio -V_r/V_T for a fixed reset depth, V_T - V_r has a strong influence on the shape of the trajectory of the phase-response function of the neuron, as shown in Figure 8. In Figure 8A, the membrane potential of a QIF neuron firing at 50 Hz is plotted for three values of the ratio V_r/V_T and a fixed reset depth: V_T - V_r = 3. For -V_r/V_T = 0.1, the concavity of the subthreshold trajectory is always directed upward. In contrast, for -V_r/V_T = 10, it is always downward. For -V_r/V_T = 1, the concavity changes from upward to downward.

View larger version (21K): [in this window] [in a new window]   Figure 8. Voltage traces and phase-response functions of the QIF neuron. The voltage traces and the phase-response functions were computed using Equations 13 and 15, respectively. 0 = 10 msec. In A and B, the external current is such that the firing rate is 50 Hz. A, The voltage trace of the neuron in response to a constant injected current for three values of the ratio -Vr/VT. The reset depth is the same in all three cases: VT - Vr = 3. Note the changes in the concavity of the voltage trace between the spikes when -Vr/VT increases. The external current is Iext = 0.97 for -Vr/VT = 1 and for -Vr/VT = 10, and Iext = 0.66 for -Vr/VT = 0.1. B, The phase-response function, Z, for the same cases as in A. C, Same as in B, but for a firing rate of 10 Hz. The phase-response changes less than in B when -Vr/VT varies.

 

It is easy to see that Z() > 0 for all , and that it is with one maximum. The location of the maximum of Z depends on V_r, V_T, and I_ext. This is depicted in Figure 8, where the voltage trace and the response function are plotted for different values of these parameters. For -V_r/V_T = 1, Z() is symmetric and its maximum is always at = . For -V_r/V_T < 1, the maximum of Z is located in the first half of the firing period, whereas for -V_r/V_T > 1, it is located in the second half. The value of Z at the maximum also depends on V_r, V_T. It has smaller values when -V_r/V_T = 1. The larger the firing rate, the stronger the dependence of Z on -V_r/V_T (Fig. 8, compare B, for 50 Hz, and C, for 10 Hz).

A pair of identical QIF neurons without noise
When, in the absence of noise, two identical neurons interact in a symmetric way, they reach, at large time, a phase-locked state in which the two neurons fire action potentials with a fixed phase shift, . In the Appendix, it is shown that is a solution to the equation:

(16)

where the function _-() 1/2(() - (-)), and () is the effective phase coupling (see Materials and Methods). Because of the symmetry of the system and the 2 periodicity of the function _-, there are always at least two solutions to this equation: = 0 and = . The first solution corresponds to the state in which the two neurons fire in-phase (in-phase locking), the second one to the state in which they fire in antiphase (antiphase locking). Other solutions with 0 < < may also exist. However, only solutions that are stable can be reached at large time. Therefore, out of all of the solutions to Equation 16, only those that satisfy the condition (see Appendix):

(17)

correspond to phase-locked states that the two neurons can eventually reach starting from appropriate initial conditions. In the following, we focus on the stability of the antiphase state of a pair of QIF neurons.

Stability conditions of antiphase locking
The voltage trajectory consists of three parts: (1) a subthreshold part, during which the membrane potential increases with time from a value V_r to a threshold value V_T; (2) a spike, modeled by a -function of amplitude , and (3) a resetting, where potential is instantaneously brought back to V_r. Then, as explained in the Appendix, antiphase locking of a pair of QIF neurons is stable if:

(18)

with:

(19)

(20)

(21)

The first term, S_sp, is the contribution of the presynaptic spikes to the destabilization of antiphase locking. The second term, S_r, is the contribution of the instantaneous reset of the membrane potential. The last term, S_sub, corresponds to the effect of the coupling between the two neurons when the presynaptic neuron is subthreshold. These three terms are plotted in Figure 9A as a function of the ratio -V_r/V_T, for = 1, V_T - V_r = 3, and a frequency v = 50 Hz.

View larger version (29K): [in this window] [in a new window]   Figure 9. Stability of the asynchronous state of QIF networks at weak coupling in the absence of noise. Results from the analytical study of the QIF model. Time constant0 = 10 msec. The size of the spikes is = 1. A, A pair of identical neurons. The quantity '() is plotted (solid line). The three contributions to this quantity (Eq. 21) are also plotted: Ssub (dotted line), Sr (dashed-dotted line), and Ssp (dashed line). The double-dotted-dashed line corresponds to Sr + Ssub. All of these quantities were computed from Equations 21, 13, and 15. Note the logarithmic scale on the x-axis. The neurons fire at 50 Hz. B, A large network of all-to-all weakly coupled identical neurons: the quantity µ1. The asynchronous state is stable with respect to the mode n = 1 if µ1 < 0. The quantities µ' 1 (dashed line) and µ ' 1 (double-dotted-dashed line) are also plotted. The firing rate of the neurons is 50 Hz. C, The phase diagram of antiphase locking for a pair of identical neurons in the weak coupling limit. Note the logarithmic scale on the x-axis. The solid lines separate the region in which antiphase locking is unstable (in the middle) from the regions in which it is stable (on the left and in the small region on the right). The parameters Vr and VT are varied, whereas the reset depth is kept constant: VT - Vr = 3. The boundary lines are defined by the condition: '() = 0. D, The phase diagram for the stability of the asynchronous state in a large network of all-to-all coupled neurons. The solid line is the boundary between the domains in which the asynchronous state is stable (to the left) and unstable (to the right).

 

The terms S_r and S_sub (Fig. 9, dotted line and dashed-dotted line, respectively) are invariant under the transformation -V_r/V_T V_T/-V_r. That is why the corresponding curves in Figure 8B are symmetric around -V_r/V_T = 1 (note the logarithmic scale of the x-axis). S_r is always positive and therefore destabilizes antiphase locking. In contrast, S_sub is always negative and leads to stabilized antiphase locking. The sum of the two terms is also plotted in Figure 9A (double-dotted-dashed line) to show that their overall contribution is destabilizing except for very small or very large values of -V_r/V_T.

Because the maximum of Z moves from the first part of the period to the second part of the period when -V_r/V_T crosses 1 (Fig. 8), the sign of the derivative of Z at = also changes from negative to positive, and the term S_sp (dashed line) tends to destabilize antiphase locking for -V_r/V_T > 1. The sum of all three terms of the right side of Equation 18 is plotted in Figure 9A (solid line). Antiphase locking is stable only if the phase at which Z reaches its maximum is small enough.

The phase diagram for the stability of antiphase locking as a function of -V_r/V_T and the firing frequencies is displayed in Figure 9C. The stability of antiphase locking requires that -V_r/V_T be small enough. Moreover, the range of values of -V_r/V_T where it is stable increases with the frequency of the neurons. This is because the contribution of the term S_sp in the right side of Equation 18 grows with the firing rate and tends to dominate the contribution of the two other terms. In contrast, at low frequencies, the response function weakly depends on -V_r/V_T (Fig. 8), and S_sp, which is proportional to the frequency, is negligible. Moreover, S_r + S_sub is positive except for very small or very large values of -V_r/V_T. Therefore, for small frequencies, the domain of stability of antiphase locking is very narrow.

This analysis shows that for a pair of neurons, the stability of antiphase locking depends critically on the shape of the phase-response function and in particular on the sign of the derivative of the response function at the mid-period. The more skewed toward the right (the second part of the firing period) the response function of the neurons, the more unstable antiphase locking becomes.

Large networks of weakly coupled QIF neurons
The stability analysis of the asynchronous state of a large network of neurons coupled all-to-all receiving a noisy input is also simplified if one assumes weak coupling. Using the phase-reduction approach, it can be shown (see Appendix) that the asynchronous state is stable if:

(22)

For all integers, n > 0, with Z_n and v_n the nth-Fourier components of the function Z() and v(), respectively (Eqs. 13, 15) defined by and a similar equation for v_n. If for some n, µ_n is negative (respectively positive), it is said that the mode of order n is stable (respectively unstable). For the asynchronous state to be stable, the modes at any order must be stable. It is unstable when at least one mode is unstable.

The first two terms in Equation 22 correspond to the effect of the interaction. The first term represents the effect of the spikes and can be written:

(23)

The second term, µ ' _n = ng Im(nZ_nv_-n), corresponds to the combined effect of the reset of the membrane potential and its subsequent subthreshold evolution.

The sign of their sum depends on the parameters, V_r, V_T and on the frequency of the neurons, v. The last term in Equation 22 corresponds to the effect of the noise. It is always positive. Therefore, as should be expected, noise increases the stability of the asynchronous state (because of the negative sign in front of this term in Eq. 22). The stability of the asynchronous state depends on the competition between the first two terms and the last term. In particular, the sign of µ_n depends on the ratio ²/g. Note that because of the factor n² in the last term of Equation 22, the stability of the modes increases rapidly with their order.

We first consider the case of a noiseless network (i.e., with = 0). In Figure 9B, we plotted µ' ₁ (dashed line) and µ ' ₁ (dashed-dotted line) against -V_r/V_T. The qualitative behavior of these quantities is similar to those of S_sp and S_sub + S_r, respectively, for a pair of neurons. In particular µ ' ₁ is symmetric (in a semilogarithmic scale) around -V_r/V_T = 1. Moreover, µ' ₁ increases monotonically from negative values to positive values when -V_r/V_T increases, and it changes sign at -V_r/V_T = 1. This can be easily explained. Indeed, for -V_r/V_T > 1 (respectively, -V_r/V_T < 1), the function Z() is skewed toward < (respectively > ) where sin() > 0 [respectively sin() < 0]; therefore, for n = 1, the integral in Equation 23 is negative (respectively positive) and µ' ₁ < 0 (respectively µ' ₁ > 0). In particular, because for -V_r/V_T = 1, Z()is symmetric around = , the integral in Equation 23 vanishes for n = 1 [because sin( + ) = -sin()]. The sign of µ₁ (solid line) is primarily determined by the sign of µ' ₁. It is negative for small -V_r/V_T and becomes positive for -V_r/V_T 1. Therefore the mode n = 1 is stable only if -V_r/V_T is small enough. Similar analysis can be performed for the modes n > 1. It reveals that if -V_r/V_T is not too large, the mode n = 1 determines the stability of the asynchronous state.

The phase diagram for the stability of the asynchronous state as a function of -V_r/V_T and the firing frequency is plotted in Figure 9D. It is very similar to the phase diagram for the stability of the antiphase state for a pair of neurons (Fig. 9C). The only qualitative difference between these two phase diagrams is that for a pair of neurons, antiphase locking is stable at low firing rates and very large -V_r/V_T, whereas for a large network, the asynchronous state is unstable in that limit. Therefore, we can conclude that neurons coupled with electrical synapses will be more easily synchronized if their phase-response functions are skewed to the right than to the left.

The phase diagram for the stability of the asynchronous state is plotted in Figure 10 for different levels of noise. For large-enough firing rates, the instability lines for ²/g = 0 and ²/g = 0.01 are very close to each other. The distance between the two lines increases when the frequency decreases. At some value, v = v^*, the two lines separate completely. At v^*, the line for ²/g = 0.01 changes direction and continues toward the right of the phase diagram (large -V_r/V_T). In contrast, the line for = 0 continues toward the left of the diagram. In particular, for ²/g = 0.01, the asynchronous state is stable for any value of -V_r/V_T for small-enough firing rates. As a matter of fact, one can show analytically that the asynchronous state is stable for any value of -V_r/V_T if the firing rate is smaller than a critical value, v^*( ²/g), which vanishes with ²/g. The size of the region in which the asynchronous state is stable increases with the noise level, as expected. In particular, v^*, and the frequency range in which the asynchronous state is stable for all -V_r/V_T, increases with ²/g.

View larger version (19K): [in this window] [in a new window]   Figure 10. Stability of the asynchronous (AS) state in a large network of all-to-all weakly coupled QIF neurons in the presence of noise. Time constant 0 = 10 msec. The size of the spikes is = 1. The solid line separates the region in which the asynchronous state is stable from the region in which it is unstable for 2/g = 0.05. The dashed line is for 2/g = 0.01. For comparison, the instability boundary is also plotted for the noiseless case (dotted line). The parameters Vr and VT are varied, but the reset depth is kept constant: VT - Vr = 3. Note the logarithmic scale on the x-axis.

 

How synchrony depends on the ratio -V_r/V_T at fixed frequencies can be deduced from the phase diagram in Figure 10. For low and moderate noise levels (e.g., ²/g = 0.01 and ²/g = 0.05), we find three generic behaviors as a function of -V_r/V_T: (1) At low firing rates, the asynchronous state is stable for all -V_r/V_T. (2) At high firing rates, the asynchronous state is stable for small -V_r/V_T and becomes unstable when -V_r/V_T increases. (3) In some intermediate range of firing rates, the stability of the asynchronous state varies nonmonotonously with -V_r/V_T: when -V_r/V_T increases from 0, the asynchronous state is first stable, then unstable, and stable again for large -V_r/V_T. The size of the intermediate region in which the asynchronous state is unstable decreases when the noise level increases.

Figure 10 also allows us to predict how the stability of the asynchronous state depends on the frequency for fixed -V_r/V_T. For small-enough -V_r/V_T, the asynchronous state is stable at all firing rates. The range of -V_r/V_T where this happens is larger when the noise is stronger. For large-enough -V_r/V_T, the asynchronous state is stable at low firing rates and becomes unstable when the firing rate is large enough. If the noise is not too strong, the stability of the asynchronous state varies nonmonotonously in some intermediate range of -V_r/V_T < 1. In this domain, the asynchronous state is stable at low rates, loses stability when the rate increases, but becomes stable again at high firing rates.

Interpretation of the simulation results of the conductance-based network
Our analysis of the synchrony properties of QIF networks shows that the stability of the asynchronous states depends crucially on the shape of the phase-response function, Z(), and in particular on the location of its maximum. This suggests that the desynchronization effect of the persistent sodium and the synchronization effect of the potassium currents revealed by our numerical simulations may be related to the way these currents shape the phase-response functions of the neurons.

The function, Z(), for our conductance-based neuron model is plotted in Figure 11 for different values of g_K, g_Ks, and g_NaP. In Figure 11A–C, Z() > 0 for all , except in a very tiny region after the spike peak. In Figure 11D, the region of negative Z() is larger. However, in all four cases, negative values of Z() are much smaller than positive values. This means that most of the time, a small and brief depolarizing perturbation advances subsequent firing of action potentials.

View larger version (20K): [in this window] [in a new window]   Figure 11. The phase-response functions depend on the intrinsic currents. The conductance-based model is defined in Materials and Methods. The phase-response function was computed for the same four sets of intrinsic conductances as in Figure 1 using XPPAUT software (Ermentrout, 2002). For each set of parameters, the external current is such that the firing frequency of the neuron is 50 Hz. A, Control cases: gK = 9 mS/cm 2, gNaP = gKs = 0. External current is Iext = 1.10 µA/cm 2. B, The conductance of the delayed rectifier is reduced to gK = 2.5 mS/cm 2 (gNaP = 0, gKs = 0 as in the control case). External current is Iext = 0.48 µA/cm 2. C, gNaP = 0.2 mS/cm 2 [other currents as in the control (gK = 9 mS/cm 2, gKs = 0)]. Externalcurrent isIext = -0.55 µA/cm 2. D, The parameters are as inB, butgKs = 0.2 mS/cm 2. External current is Iext = 4.88 µA/cm 2.

 

In all of the cases presented in Figure 11, the function Z() has only one maximum, whose location depends on the intrinsic currents. Comparison of Figure 11A with Figure 11B shows that the maximum of Z is shifted to the right when g_K increases, and that the global shape of the phase-response function is more skewed to the right (respectively to the left) for large (respectively small) values of g_K. The slow potassium current has an effect similar to that seen for the delayed rectifier current. It distorts the phase-response function toward the second half of the firing period. Therefore, increasing the potassium conductances has a similar effect on the phase-response function of the conductance-based neuron as increasing -V_r/V_T does for the QIF neuron. In contrast, the persistent sodium current shifts the peak of the response to the left and skews the shape of the phase-response function in that direction, as seen by comparing Figure 11A with Figure 11C. Hence, increasing g_NaP in the conductance-based model and decreasing -V_r/V_T in the QIF model have the same qualitative effect.

We are now able to understand the results of our numerical simulations. In Figure 6, A and B, the network state is asynchronous for small potassium conductances (both delay-rectifier and slow-potassium) and becomes synchronous when these conductances are increased. A similar behavior is found in the QIF network for v > v_*(): the asynchronous state is stable when -V^r/V^T is small, and the asynchronous state is destabilized when -V_r/V_T increases enough. This stems from the fact that the qualitative changes in the phase-response function of the conductance-based neuron when g_K increases and the QIF neuron when -V_r/V_T increases are qualitatively similar. Conversely, increasing the conductance of the persistent sodium current has the same effect on the phase-response function of the neuron as a decrease of -V_r/V_T. This explains the results of Figure 6C.

The phase diagram of Figure 10 also explains the different behaviors we found in our simulations for fixed-current conductances when the frequency varies (Fig. 7) in presence of noise. For instance, in the control case, the response function is skewed to the right. This is equivalent to -V_r/V_T > 1 for the QIF model. This explains the monotonous transition from the asynchronous state at low firing rates to synchrony when the firing rate increases, as shown in Figure 7, A and D. The nonmonotonous behavior in Figure 7, B and C, can also be explained. In Figure 7C, g_K is smaller than in the control model. The corresponding effective value of -V_r/V_T is in the intermediate range, where the asynchronous state changes stability nonmonotonously with the frequency. A similar argument can account for the results in Figure 7B, but here the change in the phase-response function is attributable to the persistent sodium current.

	Discussion

Top Abstract Introduction Materials and Methods Results Discussion Appendix References

 
How intrinsic currents affect synchrony in networks of neurons connected by electrical synapses
In this work, we have derived explicit rules for the stability of the asynchronous state in networks of neurons interacting via electrical synapses. We have shown analytically how the shape of the PRF, which depends on the model parameters, determines the stability of the asynchronous state in QIF networks. We have extended these results to provide a unified framework explaining the diversity of behavior found in our numerical simulations when potassium and sodium conductances are changed.

One can understand intuitively how potassium and sodium currents shape the PRF. The current, I_K, increases the refractoriness of the neuron; therefore, it reduces its responsiveness to small depolarization after a spike. A similar effect, but stronger and more lasting, occurs with I_Ks. This intuitively explains why we found that larger conductances of these currents skew the PRF toward the second half of the firing period. Similar effects have been found by Ermentrout et al. (2001) for M current and AHP potassium current. The current I_NaP is an inward current. It is already activated near rest, and depolarizing perturbations amplify this activation. Hence I_NaP increases the responsiveness of the neuron after a spike and shifts the maximum of the response function toward the first half of the period, as we have found here.

We predict that potassium currents, like I_K and I_Ks, tend to promote synchrony of neurons coupled via electrical synapses and that, in contrast, I_NaP tends to oppose it. Our other predictions concern the dependency of the synchrony level with the firing rate. We expect that these conclusions do not depend on the particular models we have chosen for the intrinsic currents.

Assumption of weak coupling
The stability of the asynchronous state of a large network of all-to-all connected QIF neurons can be calculated analytically in the absence of noise for any coupling strength using the population density method (Abbott and van Vreeswijk, 1993; Hansel and Mato, 2003). An alternative approach is the phase-reduction method. Although it is only mathematically justified for weak interactions, it leads to nontrivial results that remain valid in a reasonable range of interaction strengths. Moreover, it is straightforward to study the effect of noise in this framework. These were the motivations to use the phase-reduction method in the framework of which all of the analytical results of this paper were derived.

At finite coupling strength, an additional instability of the asynchronous state appears at low firing rates and small -V_r/V_T. It corresponds to the impossibility of controlling low firing rates when the coupling is too strong and the AHP is too small. This is because for a small AHP, the recurrent electrical synapse interactions tend to depolarize the neurons on the average, and this increases their firing rates. This instability is similar to the rate instability in networks of excitatory neurons (Hansel and Mato, 2003).

A detailed analysis of the QIF network at finite coupling strength shows that in a broad range of coupling strengths, the phase diagrams for two coupled neurons and for large networks are similar to those derived here, except for the additional rate instability line. A full report of this analysis will be published elsewhere.

Simulation results were performed here for a coupling conductance, g = 0.005 mS/cm². This corresponds to coupling coefficients (Amitai et al., 2002) CC 5% and to a spikelet size of 0.5 mV (Fig. 3). We have verified that similar results have been found for conductances four times larger (data not shown). Therefore the conclusions of the present work are relevant for electrical synapse conductances in the physiological range (Galarreta and Hestrin, 2001, 2002; Amitai et al., 2002).

Related works
Stable antiphase locking has been demonstrated in previous studies for a pair of neurons coupled via electrical synapses. Sherman and Rinzel (1992) used numerical simulations to show that this nonintuitive phenomenon occurs in a simple conductance-based model. However, they did not relate it to the intrinsic properties of their model. Antiphase locking was found by Han et al. (1995) in simulations of coupled Morris–Lecar oscillators. For this model, they computed the effective phase interaction, , and showed that their simulation results could be predicted from it. Our work further clarifies the conditions of stable antiphase locking in that it gives general qualitative criteria that the PRF must satisfy, which allows us to predict how it depends on intrinsic currents.

Phase locking of neurons connected by electrical synapses has been also investigated analytically (Chow and Kopell, 2000; Lewis and Rinzel, 2003). In both studies, neurons were modeled in the subthreshold range as leaky integrators, but the spikes were described in different ways. Chow and Kopell (2000) used a phenomenological model to describe the time course of spikes, whereas Lewis and Rinzel (2003), like us, modeled them with a -function. Both studies found stable antiphase locking at low firing rates and destabilization for large-enough firing rates. This can be understood using Equations 18–21 as follows. For a passive integrator, Z(), Z() = 1/[I_extT] exp(T/2_m), where T is the firing rate, I_ext is the external current, and _m is the neuronal membrane time constant (Kuramoto, 1991; Hansel et al., 1995). The value of Z diverges exponentially with the period at all , as do S_sub and S_r. Detailed analysis reveals that the most divergent of these terms is S_sub < 0. Hence, for low-enough firing rates, '() < 0. The term S_sp is positive (Z'() > 0) and decreases with T, whereas S_sub and S_r increase. Hence, for a large-enough frequency, '() changes sign. Therefore for passive integrators, antiphase locking is stable at low rates and loses stability as the firing rate increases. In contrast, we have found that the QIF model behaves differently than the LIF except in the limit of very large -V_r/V_T. Therefore our results show that the predictions of Chow and Kopell (2000) and Lewis and Rinzel (2003) may be relevant only for neurons with strong potassium currents.

Relationship to experiments
Electrical synapses connect fast spiking (FS) as well as low-threshold spiking (LTS) interneurons (Kawaguchi and Kubota, 1997; Gibson et al., 1999). Recently, the synchrony level of pairs of LTS and FS cells (characterized by the amplitude of the peak of their spikes cross-correlation) has been studied as a function of the firing rate v (Mancilla et al., 2002). It has been found that for FS cells, synchrony increases with v but decreases for LTS cells. Our results may explain this difference if one assumes that I_NaP is stronger in LTS than in FS cells and/or that potassium conductances are stronger in FS than in LTS cells. This assumption is compatible with the fact that FS and LTS cells have different firing patterns: in response to current injection, the action potentials of FS but not LTS cells are followed by a pronounced AHP, which may reflect stronger potassium conductances in FS than in LTS cells. A test of this explanation would be to show that the PRFs of these LTS and FS cells have different shapes.

The dynamic clamp technique (Sharp et al., 1993) makes it possible to artificially manipulate the intrinsic currents of the neurons and to couple two neurons via artificial electrical synapses. It provides a nice way of further testing our predictions. Therefore, it enables systematic verification of the effects of sodium and potassium currents on the synchronization properties of a pair of neurons (Fig. 4) and how they correlate with their PRF.

Several studies in the mouse and the rat have reported the presence of electrical synapses in hypoglossal and other inspiratory brainstem and phrenic motoneurons (Mazza et al., 1992; Rekling and Feldman, 1998). Recently, it has been found that inspiratory motoneuron synchrony is modulated by electrical synapses: blocking these synapses increases their level of synchrony (Bou-Flores and Berger, 2001). An explanation of the surprising findings of Bou-Flores and Berger (2001) is suggested by our work. Indeed, since strong I_NaP (Powers and Binder, 2003) is present in these neurons, electrical synapses can reduce the level of synchrony of their activity compared with the case in which electrical synapses are blocked. Synchrony in the latter situation would be attributable, for example, to GABAergic synaptic interactions between the neurons or to a spatially correlated time-dependent external input. As a matter of fact, we have verified that this effect occurs in numerical simulations of our network model with GABAergic synapses added. A detailed study of this phenomenon will be reported elsewhere.

	Appendix

Top Abstract Introduction Materials and Methods Results Discussion Appendix References

 
The I–F curve, the membrane potential, and the phase-response function of the QIF neuron
The subthreshold membrane potential, v(t), of a QIF neuron satisfies the dynamic equation:

(24)

Integrating this differential equation, one finds that its general solution is for I_ext > I_c:

(25)

where is a constant of integration. The condition that at t = 0 the membrane potential of the neuron is at its reset value, V_r, determines . One finds:

(26)

The function v(t) increases monotonously to infinity. Therefore, after some time T, v(t) reaches the threshold, V_T. At that time, the neuron fires an action potential, and v(t) is instantaneously reset to V_r. Subsequently, v(t) starts again to increase until it again reaches the threshold after another time T. Therefore, the neuron is firing periodically. The condition: v(T) = V_T determines the firing period. One finds:

(27)

The I–F curve of the QIF neuron is given by v = 1/T(I_ext). When the subthreshold membrane potential, v(t) reaches V_T, the neuron fires an action potential represented by a -function:

(28)

where measures the integral over time of the suprathreshold part of an action potential.

The response function Z is the phase resetting curve of the neuron (Kuramoto, 1984; Hansel et al., 1993; Rinzel and Ermentrout, 1998) in the limit of vanishing small perturbations of the membrane potential. For the case of one-dimensional models, it can be written as (Kuramoto, 1991; Hansel et al., 1995). Using Equation 24, one finds for the QIF model:

(29)

Phase interaction for electrical synapses
The phase coupling function, (), between two identical neurons, i and j, interacting synaptically depends on the phase difference, _i - _j, of the two neurons and is given by (Kuramoto, 1984; Hansel et al., 1993, 1995; Kopell and Ermentrout, 2002):

(30)

where Z is the phase-response function of the neurons and I_syn(_i, _j) is the synaptic current coming from the presynaptic neuron j to the postsynaptic neuron i. For electrical synapses, this current is:

(31)

Phase locking of a pair of weakly interacting neurons
When, in the absence of noise, two identical neurons interact weakly in a symmetric way, they reach a phase-locked state at large time. This is easily shown in the phase-reduction framework (Hansel et al., 1993, 1995; Van Vreeswick et al., 1994; Kopell and Ermentrout, 2002). The phases, ₁ and ₂, of the two neurons satisfy the two coupled equations:

(32)

(33)

Subtracting these equations, one finds:

(34)

where is the phase shift between the two neurons, = ₁ - ₂, and _-() = ¹₂(() = - (-)). At large time, the phase shift, , reaches a fixed point (i.e., it is such that d/dt = 0) and the two neurons become locked with a phase shift that satisfies the condition:

(35)

Because of the symmetry of the system and the 2 periodicity of the function _-, there are always at least two solutions to this equation. One is the in-phase solution, = 0, and the other is the antiphase solution: = . Other solutions with 0 < < may also exist. However, only those that are stable can be reached at large time. Linearizing Equation 34 around these solutions, one finds that the stability condition is:

(36)

where combining Equations 30 and 31:

(37)

Differentiating this function with respect to , one finds:

(38)

The membrane potential V can be expanded in its components (spike, reset, and subthreshold component) as:

(39)

which allows us to rewrite '():

(40)

Integrating by parts the first term and using the properties of the -function, we finally obtain:

(41)

n neurons: stability of the asynchronous state for the QIF network with noise
The stability analysis of the asynchronous state can be investigated for a network of phase oscillators in the presence of white noise (zero mean and variance D). For a general phase-coupling interaction, (), one can show (Kuramoto, 1984) that the asynchronous state is stable if for all n > 0:

(42)

where we have defined:

(43)

with _n, the nth-Fourier component of the 2 periodic function ():

(44)

For a network of QIF neurons interacting with electrical synapses (Eqs. 24, 28) the nth-Fourier component of phase interaction has the form:

(45)

where Z_n and v_n are the nth-Fourier components of the functions Z() and v() defined in Equations 29 and 25. The variance of the noise, D, depends on the average of the phase-response function, and is given by D = ²Z₀² (Kuramoto, 1984), where ² is the variance of the white noise in the full QIF model. This yields the stability condition for the asynchronous state:

(46)

Note that Z₀ depends on the parameters V_r, V_T and on the firing frequency of the neurons. Therefore, the desynchronizing effect of the noise depends on these parameters. In particular, because Z₀ decreases when the firing rate increases, the noise is more efficient at stabilizing the asynchronous state at low firing rates than at high firing rates.

	Footnotes

 
Received Mar. 13, 2003; revised Apr. 21, 2003; accepted Apr. 22, 2003.

This work was supported by North Atlantic Treaty Organization Physical and Engineering Science and Technology Collaborative Linkage Grant 977683, Les Programmes Internationaux de Coopération Scientifique–Centre National de la Recherche Scientifique (number 837), L'Action Concertée Incitative "Neurosciences integratives et computationnelles" (Ministère de la Recherche, France), and project A99E01 from the Comisión Asesora Científica de la Cooperación Argentino-Francesa. Research by D.G. was supported by Israel Science Foundation Grant 657/01 and by the National Science Foundation, under Agreement 0112050. We thank Y. Loewenstein, C. van Vreeswijk, B. Connors, and Y. Yarom for stimulating discussions and B. Connors, Y. Loewenstein, and C. Meunier for careful and critical reading of this manuscript.

Correspondence should be addressed to Dr. David Hansel, Laboratoire de Neurophysique et Physiologie du Système Moteur, 45 Rue des Saints-Pères, 75270 Paris Cedex 06, France. E-mail: david.hansel{at}biomedicale.univparis5.fr.

Copyright © 2003 Society for Neuroscience 0270-6474/03/236280-15$15.00/0

	References

Top Abstract Introduction Materials and Methods Results Discussion Appendix References

 

Abbott LF, van Vreeswijk C (1993) Asynchronous state in networks of pulse-coupled oscillators. Phys Rev E 48: 1483–1490.

Abbott LF, Hooper SL, Marder E (1991) Oscillating networks: control of burst duration by electrically coupled neurons. Neural Comput 3: 487–497.

Alvarez LF, Chow C, van Bockstaele EJ, Williams JT (2002) Frequency-dependent synchrony in locus coeruleus: role of electrotonic coupling. Proc Natl Acad Sci USA 99: 4032–4036.[Abstract/Free Full Text]

Amitai Y, Gibson JR, Patrick A, Ho B, Connors BW, Golomb D (2002) Spatial organization of electrically coupled network of interneurons in neocortex. J Neurosci 22: 4142–4152.[Abstract/Free Full Text]

Beierlein M, Gibson JR, Connors BW (2000) A network of electrically coupled interneurons drives synchronized inhibition in neocortex. Nat Neurosci 3: 904–910.[Web of Science][Medline]

Bou-Flores C, Berger AJ (2001) Gap junctions and inhibitory synapses modulate inspiratory motoneuron synchronization. J Neurophysiol 85: 1543–1551.[Abstract/Free Full Text]

Chow CC, Kopell N (2000) Dynamics of spiking neurons with electrical coupling. Neural Comput 12: 1643–1678.[Web of Science][Medline]

Christie MJ, Williams JT, North RA (1989) Electrical coupling synchronizes subthreshold activity in locus coeruleus neurons in vitro from neonatal rats. J Neurosci 9: 3584–3589.[Abstract]

Cymbalyuk GS, Nikolaev EV, Borisyuk RM (1994) In-phase and antiphase self-oscillations in a model of two electrically coupled pacemakers. Biol Cybern 71: 153–160.[Web of Science][Medline]

Deans MR, Gibson JR, Sellitto C, Connors BW, Paul DL (2001) Synchronous activity of inhibitory networks in neocortex requires electrical synapses containing connexin36. Neuron 31: 477–485.[Web of Science][Medline]

Draghun A, Traub RD, Schmitz D, Jefferys JGR (1998) Electrical coupling underlies high-frequency oscillations in the hippocampus in vitro. Nature 394: 189–192.[Medline]

Erisir A, Lau D, Rudy B, Leonard CS (1999) Function of specific K⁺ channels in sustained high-frequency firing of fast-spiking neocortical interneurons. J Neurophysiol 82: 2476–2489.[Abstract/Free Full Text]

Ermentrout B (1996) Type I membranes, phase resetting curves and synchrony. Neural Comput 8: 979–1001.[Web of Science][Medline]

Ermentrout B (2002) Simulating, analyzing, and animating dynamical systems: a guide to Xppaut for researchers and students (software, environments, tools). Philadelphia: Society for Industrial and Applied Mathematics.

Ermentrout B, Kopell N (1986) Parabolic bursting in an excitable system coupled with a slow oscillation. SIAM J Appl Math 46: 233–253.

Ermentrout B, Pascal M, Gutkin B (2001) The effect of spike frequency adaptation and negative feedback on the synchronization of neural oscillators. Neural Comput 13: 1285–1310.[Web of Science][Medline]

French CR, Sah P, Buckett KJ, Gage PW (1990) A voltage-dependent persistent sodium current in mammalian hippocampal neurons. J Gen Physiol 95: 1139–1157.[Abstract/Free Full Text]

Furshpan EJ, Potter DD (1959) Transmission at the giant motor synapses of the crayfish. J Physiol (Lond) 145: 289–325.

Galarreta M, Hestrin S (1999) A network of fast spiking cells in the neocortex connected by electrical synapses. Nature 402: 72–75.[Medline]

Galarreta M, Hestrin S (2001) Electrical synapses between GABA-releasing neurons. Nat Neurosci 2: 425–433.[Web of Science][Medline]

Galarreta M, Hestrin S (2002) Electrical and chemical synapses among parvalbumin fast-spiking GABAergic interneurons in adult mouse neocortex. Proc Natl Acad Sci USA 19: 12438–12443.

Gibson JR, Beierlein M, Connors B (1999) Two networks of inhibitory neurons electrically coupled. Nature 402: 75–79.[Medline]

Golomb D, Hansel D (2000) The number of synaptic inputs and the synchrony of large sparse neuronal networks. Neural Comput 12: 1095–1139.[Web of Science][Medline]

Golomb D, Rinzel J (1994) Clustering in globally coupled inhibitory neurons. Physica D 72: 259–282.

Golomb D, Hansel D, Mato G (2001) Theory of synchrony of neuronal activity. In: Handbook of biological physics, Vol 4: Neuro-informatics and neural modeling (Gielen S, Moss M, eds), pp 887–968. Amsterdam: Elsevier.

Han SK, Kurrer C, Kuramoto Y (1995) Dephasing and bursting in coupled neural oscillators. Phys Rev Lett 75: 3190–3193.[Web of Science][Medline]

Hansel D, Mato G (2003) Asynchronous states and the emergence of synchrony in large networks of interacting excitatory and inhibitory neurons. Neural Comput 15: 1–56.[Web of Science][Medline]

Hansel D, Sompolinsky H (1992) Synchrony and computation in a chaotic neural network. Phys Rev Lett 68: 718–721.[Web of Science][Medline]

Hansel D, Mato G, Meunier C (1993) Phase dynamics for weakly coupled Hodgkin-Huxley neurons. Europhys Lett 23: 367–372.

Hansel D, Mato G, Meunier C (1995) Synchrony in excitatory neural networks. Neural Comput 7: 307–337.[Web of Science][Medline]

Hormuzdi SG, Pais I, Lebeau FEN, Towers SK, Rozov A, Buhl E, Whittington MA, Monyer H (2001) Impaired electrical signaling disrupts gamma frequency oscillations in connexin 36-deficient mice. Neuron 31: 487–495.[Web of Science][Medline]

Kawaguchi Y, Kubota Y (1997) GABAergic cell subtypes and their synaptic connections in rat frontal cortex. Cereb Cortex 7: 476–486.[Abstract/Free Full Text]

Kepler TB, Marder E, Abbott LF (1990) The effect of electrical coupling on the frequency of model neuronal oscillators. Science 248: 83–85.[Abstract/Free Full Text]

Kiehn O, Tresch MC (2002) Gap junctions and motor behavior. Trends Neurosci 25: 108–115.[Web of Science][Medline]

Kita H, Kosaka T, Heizmann CW (1990) Parvalbumin-immunoreactive neurons in the rat neostriatum: a light and electron microscopic study. Brain Res 536: 1–15.[Web of Science][Medline]

Kopell N, Ermentrout GB (2002) Mechanisms of phase-locking and frequency control in pairs of coupled neural oscillators. In: Handbook of dynamical systems, Vol 2 (Fiedler B, ed), pp 5–55. Amsterdam: Elsevier.

Kuramoto Y (1984) Chemical oscillations, waves and turbulence. New York: Springer.

Landisman CE, Long MA, Beierlein M, Deans MR, Paul DL, Connors BW (2002) Electrical synapses in the thalamic reticular nucleus. J Neurosci 22: 1002–1009.[Abstract/Free Full Text]

Latham PE, Richmond BJ, Nelson PG, Nirenberg S (2000) Intrinsic dynamics in neuronal networks. I. Theory. J Neurophysiol 83: 808–827.[Abstract/Free Full Text]

Lewis T, Rinzel J (2003) Dynamics of spiking neurons connected by both inhibitory and electrical coupling. J Comput Neurosci 14: 283–309.[Web of Science][Medline]

Llinas R, Yarom Y (1986) Oscillatory properties of guinea-pig inferior olivary neurones and their pharmacological modulation: an in vitro study. J Physiol (Lond) 376: 163–182.[Abstract/Free Full Text]

Mancilla JG, Lewis T, Pinto DJ, Rinzel J, Connors BW (2002) Firing dynamics of single and coupled pairs of inhibitory interneurons in neocortex. 28: 840.13.

Mann-Metzer P, Yarom Y (1999) Electrical coupling interacts with intrinsic properties to generate synchronized activity in cerebellar networks of inhibitory interneurons. J Neurosci 19: 3298–3306.[Abstract/Free Full Text]

Mazza A, Nunez-Abades PA, Spielmann JM, Cameron WE (1992) Anatomical and electronical coupling in developing genioglossal motoneurons of the rat. Brain Res 598: 127–137.[Web of Science][Medline]

Meunier C (1992) The electrical coupling of two simple oscillators: load and acceleration effects. Biol Cybern 67: 155–164.[Web of Science][Medline]

Neltner L, Hansel D, Mato G, Meunier C (2000) Synchrony in heterogeneous networks of spiking neurons. Neural Comput 12: 1607–1642.[Web of Science][Medline]

Perez-Velazquez JL, Carlen PL (2000) Gap junctions, synchrony and seizures. Trends Neurosci 23: 68–74.[Web of Science][Medline]

Pfeuty B, Golomb D, Mato G, Hansel D (2002) Potassium currents modulate synchronization in networks of neurons connected with electrical synapses. 28: 840.15.

Powers RK, Binder (2003) Persistent sodium and calcium currents in rat hypoglossal motoneurons. J Neurophysiol 89: 615–624.[Abstract/Free Full Text]

Rekling JC, Feldman JL (1998) Bidirectional electrical coupling between inspiratory motoneurons in the newborn mouse nucleus ambigus. J Neurophysiol 78: 385–405.

Rinzel J, Ermentrout B (1998) Analysis of neural excitability and oscillations. In: Methods in neuronal modeling: from ions to networks, Ed 2 (Koch C, Segev I, eds), pp 251–292. Cambridge, MA: MIT.

Sharp AA, O'Neil MB, Abbott LF, Marder E (1993) The dynamic clamp: artificial conductances in biological neurons. Trends Neurosci 16: 389–394.[Web of Science][Medline]

Sherman A, Rinzel J (1992) Rhythmogenic effects of weak electrotonic coupling in neuronal models. Proc Natl Acad Sci USA 89: 2471–2474.[Abstract/Free Full Text]

Strogatz S (1994) Nonlinear dynamics and chaos. Reading, MA: Addison-Wesley.

Tamas G, Buhl EH, Lorincz A, Somogyi P (2000) Proximally targeted GABAergic synapses and gap-junctions precisely synchronize cortical interneurons. Nat Neurosci 3: 366–371.[Web of Science][Medline]

Traub RD, Kopell N, Bibbig A, Buhl E, Lebeau FEN, Whittington M (2001) Gap junction between interneuron dendrites can enhance synchrony of oscillations in distributed networks. J Neurosci 21: 9478–9486.[Abstract/Free Full Text]

Van Vreeswijk C, Abbott LF, Ermentrout GB (1994) When inhibition not excitation synchronizes neural firing. J Comput Neurosci 1: 313–321.[Medline]

Venance L, Rozov A, Blatow M, Burnashev N, Feldmeyer D, Monyer H (2000) Connexin expression in electrically coupled postnatal rat brain neurons. Proc Natl Acad Sci USA 97: 10260–10265.[Abstract/Free Full Text]

Watanabe A (1958) The interaction of electrical activity among neurons of lobster cardiac ganglion. Jpn J Physiol 8: 305–318.[Medline]

This article has been cited by other articles:

				J. Hjorth, K. T. Blackwell, and J. Hellgren Kotaleski Gap Junctions between Striatal Fast-Spiking Interneurons Regulate Spiking Activity and Synchronization as a Function of Cortical Activity J. Neurosci., April 22, 2009; 29(16): 5276 - 5286. [Abstract] [Full Text] [PDF]

				S. A. Prescott, S. Ratte, Y. De Koninck, and T. J. Sejnowski Pyramidal Neurons Switch From Integrators In Vitro to Resonators Under In Vivo-Like Conditions J Neurophysiol, December 1, 2008; 100(6): 3030 - 3042. [Abstract] [Full Text] [PDF]

				M. L. Veruki, L. Oltedal, and E. Hartveit Electrical Synapses Between AII Amacrine Cells: Dynamic Range and Functional Consequences of Variation in Junctional Conductance J Neurophysiol, December 1, 2008; 100(6): 3305 - 3322. [Abstract] [Full Text] [PDF]

				M. Vandecasteele, J. Glowinski, J.-M. Deniau, and L. Venance Chemical transmission between dopaminergic neuron pairs PNAS, March 25, 2008; 105(12): 4904 - 4909. [Abstract] [Full Text] [PDF]

				J. Gansert, J. Golowasch, and F. Nadim Sustained Rhythmic Activity in Gap-Junctionally Coupled Networks of Model Neurons Depends on the Diameter of Coupled Dendrites J Neurophysiol, December 1, 2007; 98(6): 3450 - 3460. [Abstract] [Full Text] [PDF]

				J. G. Mancilla, T. J. Lewis, D. J. Pinto, J. Rinzel, and B. W. Connors Synchronization of Electrically Coupled Pairs of Inhibitory Interneurons in Neocortex J. Neurosci., February 21, 2007; 27(8): 2058 - 2073. [Abstract] [Full Text] [PDF]

				J. A. Goldberg, C. A. Deister, and C. J. Wilson Response Properties and Synchronization of Rhythmically Firing Dendritic Neurons J Neurophysiol, January 1, 2007; 97(1): 208 - 219. [Abstract] [Full Text] [PDF]

				F. Saraga, L. Ng, and F. K. Skinner Distal Gap Junctions and Active Dendrites Can Tune Network Dynamics J Neurophysiol, March 1, 2006; 95(3): 1669 - 1682. [Abstract] [Full Text] [PDF]

				V. Zsiros and G. Maccaferri Electrical Coupling between Interneurons with Different Excitable Properties in the Stratum Lacunosum-Moleculare of the Juvenile CA1 Rat Hippocampus J. Neurosci., September 21, 2005; 25(38): 8686 - 8695. [Abstract] [Full Text] [PDF]

				B. S. Gutkin, G. B. Ermentrout, and A. D. Reyes Phase-Response Curves Give the Responses of Neurons to Transient Inputs J Neurophysiol, August 1, 2005; 94(2): 1623 - 1635. [Abstract] [Full Text] [PDF]

				P. Rabbah, J. Golowasch, and F. Nadim Effect of Electrical Coupling on Ionic Current and Synaptic Potential Measurements J Neurophysiol, July 1, 2005; 94(1): 519 - 530. [Abstract] [Full Text] [PDF]

				E. B. Merriam, T. I. Netoff, and M. I. Banks Bistable Network Behavior of Layer I Interneurons in Auditory Cortex J. Neurosci., June 29, 2005; 25(26): 6175 - 6186. [Abstract] [Full Text] [PDF]

				J. R. Gibson, M. Beierlein, and B. W. Connors Functional Properties of Electrical Synapses Between Inhibitory Interneurons of Neocortical Layer 4 J Neurophysiol, January 1, 2005; 93(1): 467 - 480. [Abstract] [Full Text] [PDF]

				S. Hidaka, Y. Akahori, and Y. Kurosawa Dendrodendritic Electrical Synapses between Mammalian Retinal Ganglion Cells J. Neurosci., November 17, 2004; 24(46): 10553 - 10567. [Abstract] [Full Text] [PDF]

				N. Kopell and B. Ermentrout Chemical and electrical synapses perform complementary roles in the synchronization of interneuronal networks PNAS, October 26, 2004; 101(43): 15482 - 15487. [Abstract] [Full Text] [PDF]

				N. Schweighofer, K. Doya, H. Fukai, J. V. Chiron, T. Furukawa, and M. Kawato Chaos may enhance information transmission in the inferior olive PNAS, March 30, 2004; 101(13): 4655 - 4660. [Abstract] [Full Text] [PDF]

				S. Schreiber, J.-M. Fellous, P. Tiesinga, and T. J. Sejnowski Influence of Ionic Conductances on Spike Timing Reliability of Cortical Neurons for Suprathreshold Rhythmic Inputs J Neurophysiol, January 1, 2004; 91(1): 194 - 205. [Abstract] [Full Text]

This Article Abstract Full Text (PDF) A correction has been published Submit an eLetter Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Web of Science (40) Citing Articles via Google Scholar Google Scholar Articles by Pfeuty, B. Articles by Hansel, D. Search for Related Content PubMed PubMed Citation Articles by Pfeuty, B. Articles by Hansel, D. Pubmed/NCBI databases Compound via MeSH Substance via MeSH Hazardous Substances DB POTASSIUM SODIUM

Home  |   Search  |   Archive  |   Subscribe  |   Contact  |   Help