JN AJP: Gastrointestinal and Liver Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 95: 1049-1067, 2006. First published October 19, 2005; doi:10.1152/jn.00932.2005
0022-3077/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
95/2/1049    most recent
00932.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (10)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Golomb, D.
Right arrow Articles by Ermentrout, G. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Golomb, D.
Right arrow Articles by Ermentrout, G. B.

Persistent Synchronized Bursting Activity in Cortical Tissues With Low Magnesium Concentration: A Modeling Study

David Golomb1, Anat Shedmi1, Rodica Curtu2,3 and G. Bard Ermentrout2

1Department of Physiology and Zlotowski Center for Neuroscience, Faculty of Health Sciences, Ben-Gurion University, Be'er-Sheva, Israel; 2Department of Mathematics, University of Pittsburgh, Pennsylvania; and 3Department of Mathematics, University of Transilvania Brasov, Brasov, Romania

Submitted 7 September 2005; accepted in final form 10 October 2005


 ABSTRACT
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We explore the mechanism of synchronized bursting activity with frequency of ~10 Hz that appears in cortical tissues at low extracellular magnesium concentration [Mg2+]o. We hypothesize that this activity is persistent, namely coexists with the quiescent state and depends on slow N-methyl-D-aspartate (NMDA) conductances. To explore this hypothesis, we construct and investigate a conductance-based model of excitatory cortical networks. Population bursting activity can persist for physiological values of the NMDA decay time constant (~100 ms). Neurons are synchronized at the time scale of bursts but not of single spikes. A reduced model of a cell coupled to itself can encompass most of this highly synchronized network behavior and is analyzed using the fast-slow method. Synchronized bursts appear for intermediate values of the NMDA conductance gNMDA if NMDA conductances are not too fast. Regular spiking activity appears for larger gNMDA. If the single cell is a conditional burster, persistent synchronized bursts become more robust. Weakly synchronized states appear for zero AMPA conductance gAMPA. Enhancing gAMPA increases both synchrony and the number of spikes within bursts and decreases the bursting frequency. Too strong gAMPA, however, prevents the activity because it enhances neuronal intrinsic adaptation. When [Mg2+]o is increased, higher gNMDA values are needed to maintain bursting activity. Bursting frequency decreases with [Mg2+]o, and the network is silent with physiological [Mg2+]o. Inhibition weakly decreases the bursting frequency if inhibitory cells receive enough NMDA-mediated excitation. This study explains the importance of conditional bursters in layer V in supporting epileptiform activity at low [Mg2+]o.


 INTRODUCTION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
A persistent state of a neuronal network is a state in which periods of increased neuronal discharge dynamically coexist with silent, or an almost silent, epochs. Such neural activity has been demonstrated in the cortex during "working memory" tasks (Funahashi et al. 1989Go); reviewed in (McCormick et al. 2003Go). A special type of persistent state, asynchronous firing, was studied using theoretical and computational methods (Amit and Brunel 1997Go; Hansel and Mato 2001Go; Wang 1999Go, 2001Go). Experiments in monkeys, however, have revealed partially synchronized activity during working memory tasks that have oscillatory (Cardoso de Oliveira et al. 2001Go) or nonoscillatory (Vaadia et al. 1995Go) characteristics. The mechanisms that generate persistent synchronized states in cortical networks are still unknown.

Here we study the mechanism for a persistent synchronized state in a specific neocortical system that was studied in vitro: discharges of rhythmic bursting activity at low extracellular Mg2+ concentration ([Mg2+]o) solution (Kawaguchi 2001Go; Sutor and Hablitz 1989Go). These synchronized oscillations depend on NMDA receptors and appear in tissue segments containing layer V only (Silva et al. 1991Go). In rats, this layer contains, among other cell types, pyramidal neurons that are quiescent at rest and burst rhythmically when constant, depolarizing current is applied, suggesting that this intrinsic bursting property could facilitate the generation of these waves (Steriade 2004Go). The oscillation frequency in the somatosensory cortex is ~8–12 Hz. Blocking GABAA-mediated inhibition considerably increases the amplitude of the local field potentials and the duration of the oscillatory episode (up to ~3 s) and only slightly increases the oscillation frequency, namely the peak of the power spectrum increases by ~10–20% (Flint and Connors 1996Go). With inhibition intact, some inhibitory neurons burst in synchrony with the excitatory cells, whereas other fire almost continuously (Kawaguchi 2001Go). Recently, properties of these waves were studied using a combination of optical, voltage-sensitive dye imaging and local field potential techniques. Episodes of bursting oscillations can be evoked or can start spontaneously at a few confined initiation foci (Tsau et al. 1998Go). Correlation between local field potentials measured at two points did not decrease considerably with distance, but there were large phase shifts at a large distance (Fig. 5 in Wu et al. 1999Go). The phase relationship between the starting times of bursting period was not constant in time (Figs. 4 and 7 in Wu et al. 1999Go), hinting that the activity after the steady-state activity is a result of collective dynamics and not of periodic stimulation of the slice by a group of pacemakers located at a certain location.


Figure 5
View larger version (28K):
[in this window]
[in a new window]
 
FIG. 5. Weakly synchronized persistent waves occur in the cortical slice model in response to a brief, local stimulation. Parameters: gAMPA = 0, gNMDA = 0.045 mS/cm2. Formula 5. A: rastergrams; Formula 5. Only firing times from every 32th cell are shown. Left: response to a local stimulation. Initially cells at left (x ≤ {sigma}) are depolarized and all others are at rest. Discharge propagates to the right. Right: system activity at steady state. B: voltage time courses of 2 adjacent neighboring cells at Formula 5. These, as well as the rastergrams, show that the bursting activity of one neuron is shifted in time in comparison with the bursting activity of the 2nd one. C: synchrony measure {chi} as a function of 1/{rho}, where {rho} varies from 64 to 1,024; {chi} is computed for a local population of 2{sigma}{rho} + 1 neurons in the middle of the slice model.

 

Figure 4
View larger version (22K):
[in this window]
[in a new window]
 
FIG. 4. Dependence of the number of spikes within a burst Ns, the bursting frequency f, the burst duration Td, and the synchrony measure {chi} on gNMDA for Formula 5(A and B) and on {tau}NMDA for Formula 5(C and D). Calculations are carried out for Formula 5(A and C) and for gAMPA = 0.08 mS/cm2 (B and D). Red lines denote results of simulations of the full network model with Formula 5. Black lines denote results of simulating the reduced model of 1 cell coupled to itself.

 

Figure 7
View larger version (23K):
[in this window]
[in a new window]
 
FIG. 7. Dependence of the number of spikes within a burst Ns, the bursting frequency f, the burst duration Td, and the synchrony measure {chi} on gAMPA (A) and gNaP (B and C). A, Formula 5; Formula 5; C, gAMPA = 0.08 mS/cm2. Red lines denote results of simulations of the full network model with Formula 5. Black lines denote results of simulating the reduced model of 1 cell coupled to itself.

 
In this paper, we investigate the simplest scenario of evoked propagating waves in a slice model without noise or disorder. We focus on dynamics of order of a few hundred milliseconds or less and do not consider the slow processes of the order of seconds or more that eventually terminate the episode of activity (therefore in reality, the activity is not strictly persistent at long times). First, we study disinhibited networks because inhibitory neurons do not have a major effect on the basic discharge properties (Flint and Connors 1996Go) of the population spiking activity. Then we explore the effect of inhibition. We hypothesize that the synchronized waves at low [Mg2+]o appear as synchronized persistent bursting activity. This hypothesis raises several questions.

First, the bursting time period of the population (field) is of the order of 100 ms. Because in each burst the active (spiking) phase of the burst eventually terminates, there is a process that terminates it that is slow in comparison with the inter-spike interval within the burst. To generate a subsequent burst, the slow N-methyl-D-aspartate (NMDA) conductance should overcome this slow process and initiate the spiking activity. For that, the NMDA-mediated conductance should decay slowly enough. How small can this decay time constant be to still support the next active phase? Is this value within the physiological range?

Second, should the single cells have a propensity to burst to obtain network bursting or can bursting be obtained as a network effect with the single cells only being regular spiking cells?

Third, does the fast, AMPA-mediated excitation support or disrupt persistent bursting? How does it affect the properties of this activity, such as bursting duration, frequency and level of synchronization?

Fourth, how large should [Mg2+]o be to prevent the persistent activity? How do the burst properties depend on [Mg2+]o in the parameter regime where they exist?

And fifth, how does intracortical inhibition affect the discharge frequency? Does the inhibitory effect depend on the type and conductance strength of excitatory receptors (AMPA and NMDA) on inhibitory cells?

To answer these questions, we construct and analyze a conductance-based model of one-dimensional networks with spatially decaying connectivity.


 METHODS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Single-cell model

EXCITATORY CELLS.  We model excitatory cortical cells that can fire either in trains of single action potentials (regular spiking, RS) or in fast bursts of action potentials (intrinsically bursting, IB) (Connors et al. 1982Go; McCormick et al. 1985Go). To evaluate the importance of the intrinsic bursting property, we use a model that can be transformed, by modifying a single parameter, from RS to IB with an increasing number of spikes within a burst (Golomb and Amitai 1997Go; D. Golomb and Y. Yaari, unpublished data). The single-cell dynamics is described by a single-compartment Hodgkin-Huxley type model with the use of a set of coupled differential equations

Formula 1(1)
Where C is the membrane capacitance and V(x,t) is the membrane potential of a neuron at a position x and time t. The right-hand side incorporates an applied current Iapp and the following intrinsic and synaptic currents: the transient sodium current INa, the persistent sodium current INaP, the delayed-rectifier potassium current IKdr, the slow potassium current IK-slow, the leak current IL, the AMPA current IAMPA, the NMDA current INMDA, and the GABAA current IGABAAIE. The currents INa and IKdr are needed for spike generation. The slow potassium current IK-slow in our model represents potassium currents with kinetics slower than the action potential time scale, i.e., with activation time scale of several tens to hundreds of milliseconds (Storm 1990Go). These currents are either calcium dependent, such as IAHP (Pinsky and Rinzel 1994Go; Sah 1996Go), or voltage dependent, such as IM (Bibbig et al. 2001Go; Mainen and Sejnowski 1996Go; Yue and Yaari 2004Go). The current IK-slow contributes the slow variable that is needed for adaptation or bursting in RS cortical cells (Rinzel and Ermentrout 1998Go). If the kinetics of INa and IKdr are fast enough, and with strong enough persistent sodium conductance (gNaP), the model may produce IB behavior in response to prolonged current injection. For example, for Formula 1, the neuron bursts if gNaP > 0.091 mS/cm2. Neurons from this type are called "conditional intrinsically bursting" (cIB). Hence, the conductance gNaP can be viewed as an "intrinsic bursting parameter." The equations and parameters of the model are given in APPENDIX A.

This single-cell model for cortical pyramidal neurons was chosen for two reasons. First, it is a minimal model that can show both regular spiking and bursting behavior. Second, the role of gNaP in the intrinsic bursting property in the model is consistent with experimental observation showing the importance of this conductance for both single-cell bursting in cortical structures (e.g., Su et al. 2001Go) and for bursting in neocortical networks bathed in low [Mg2+]o (Castro-Alamancos and Rigas 2004Go).

INHIBITORY CELLS.  We model inhibitory cortical cells using the Wang-Buzsáki model of fast spiking (FS) neurons (Wang and Buzsáki 1996Go). The single-cell dynamics is described by a single-compartment model

Formula 2(2)
The equations and parameters of the model are given in APPENDIX A.

Synaptic models

Neocortical cells receive fast, AMPA-mediated, and slow, NMDA-mediated excitatory postsynaptic potentials (EPSPs) from neighboring excitatory cells (Gil and Amitai 1996Go), and GABAA-mediated inhibitory postsynaptic potentials (IPSPs) from neighboring inhibitory cells (Hansel and Mato 2003Go; Wang and Buzsáki 1996Go). A gating variable sAMPA for an AMPA receptor, representing the fraction of open channels, is modeled according to

Formula 3(3)
where Vpre is the presynaptic potential, Formula 3Formula 3 and Formula 3(Golomb and Amitai 1997Go; Wang and Rinzel 1993Go), Formula 3. The rise time of the AMPA synapses is assumed to be instantaneous, and decay time of those synapse is Formula 3(Stern et al. 1992Go).

The rise time of NMDA receptors cannot be neglected, and therefore two differential equations are needed to model those synapses. We use the phenomenological equations of Golomb et al. (1996)Go

Formula 4(4)

Formula 5(5)
where Formula 5, Formula 5, Formula 5. The decay time of NMDA-mediated excitatory postsynaptic conductances (EPSCs) in cortex varies from 65 to 108 ms (Kumar and Huguenard 2003Go). For our reference parameter set, we use the value {tau}NMDA =100 ms, corresponding to the value of intracortically evoked NMDA-mediated EPSCs (Kumar and Huguenard 2003Go). The effect of this parameter is assessed in the following text. For simplicity, effects of synaptic depression and facilitation (Golomb and Amitai 1997Go; Tsodyks and Markram 1997Go) are not included in the model.

A gating variable sGABAA for a GABAA receptor, is modeled according to

Formula 6(6)
where Formula 6, Formula 6(Hansel and Mato 2003Go; Wang and Buzsaki 1996Go).

The total synaptic conductance a neuron receives, gSin (for AMPA, NMDA, or GABAA) is calculated by summing the synaptic variable of each of its presynaptic neurons as described in the following text. The AMPA current is

Formula 7(7)
where VGlu is the reversal potential of glutamatergic synapses. If [Mg2+]o is larger than 0, the NMDA current also depends on the postsynaptic voltage and is calculated by multiplying the presynaptic term by a sigmoid function fNMDA(V) of the postsynaptic voltage (Destexhe et al. 1994Go; Traub et al. 1991Go). The NMDA current is therefore

Formula 8(8)
The half-maximum voltage of fNMDA(V) is {theta}NMDA. The value of {theta}NMDA is –{infty} for Formula 8 and is –31 mV for NMDA-mediated intracortical synapses at physiological levels of [Mg2+]o (Kumar and Huguenard 2003Go); it increases logarithmically with [Mg2+]o (Jahr and Stevens 1990Go). Using the data for intracortical connections, we find that {theta}NMDA = 10.5 mV x ln([Mg2+]o/38.3 mM).

GABAA current is

Formula 9(9)
where VGABAA is the reversal potential of GABAA synapses.

Axonal propagation delay and synaptic delays are neglected because the time scale of bursting oscillations is in the order of 100 ms, much larger than the time scale of those delays, which are in the order of a few milliseconds.

Network architecture

Propagation of activity along cortical slices has been modeled before in one-dimensional networks (Ermentrout 1998Go; Golomb 1998Go; Golomb and Amitai 1997Go; Golomb and Ermentrout 1999Go, 2002Go; Golomb et al. 2001Go), where the extensive intralaminar connectivity is neglected, assuming cortical columns are rapidly and fully recruited. In the same manner, our model consists of a one-dimensional system. Excitatory and inhibitory neurons are equally distributed along the interval 0 ≤ x ≤ L, where L is the slice length and x is the neuron position. The cell density is {rho}, and the number of neurons from each type is Formula 9. The position of the ith neuron, either excitatory or inhibitory, 1 ≤ i ≤ N, is Formula 9. The interaction between neurons is assumed to decay with the distance between them. The "synaptic footprint shape," w(x), denotes the functional dependence of the synaptic connectivity on the distance between the preand postsynaptic cells. It is assumed here to be exponential with a characteristic delay length ("footprint length") {sigma}, namely Formula 9. The model is studied in the parameter regime 1/{rho} << {sigma} << L. The total synaptic conductance gSin affecting a neuron is

Formula 10(10)
The quantities g, s, and Sin have a subscript AMPA, NMDA, or GABAA for the three types of receptors. They also have superscripts {alpha}beta, where {alpha} and beta and denote the pre- and postsynaptic populations, respectively, that can be E (excitatory) or I (inhibitory). The coupling among neurons in the inhibitory population, either excitatory or inhibitory (Beierlein et al. 2003Go), is neglected here. To shorten the notation, we omit the superscripts if they are EE. We use open boundary conditions. The integro-differential equations (110), together with the differential equations for the auxiliary variables, are discretized on a grid, as described in APPENDIX A.

Initial conditions

As in our previous work (Golomb 1998Go; Golomb and Amitai 1997Go), we initiate our model from a state at which all the neurons, except a group at the left edge, are in their resting state. A wave is initiated by depolarizing a group of neuron, both excitatory and inhibitory, within a length {sigma} or {sigma}/2 at the left edge, such that they generate action potentials. The firing neurons may recruit resting neurons through their synaptic connections to initiate a propagating discharge.

Definitions of states

We define the state of the system based on its long-time behavior. Therefore a state of transient activity, in which the network is silent after brief transient activity, is considered to be quiescent. A transient bursting state followed by a tonic state is considered to be tonic.

For a single-cell model, or a state of a disinhibited network that can be represented by a single cell coupled to itself, the notion of "bursting cell" means here "a spiking cell that is not firing periodically (tonically)."1 To determine whether the voltage time course of an isolated neuron is tonic or bursting, we run a simulation for a long time Tm after a transient time Ttransient. We find the minimal inter-spike interval Tisi,min and the maximal inter-spike interval Tisi,max. The cell is considered to be in a tonic mode if Tisi,min/Tisi,max > 0.9. Otherwise, the cell is in a bursting mode.

The situation is more complicated in a network, where neurons can fire in an irregular and aperiodic manner. Therefore a neuron in a network is defined to be "bursting" if there is a clear distinction between brief and prolonged inter-spike interval, namely Tisi,min/Tisi,max < 0.33. The properties of the network are determined by considering the middle group of excitatory neurons with length L/2 and ignoring the two edges of length L/4. A network is considered to be in a tonic mode if all those middle excitatory neurons fire tonically. Borders between regimes of activity are determined using the bisection method (Golomb and Ermentrout 1999Go).

For both single neurons and networks, we compute the average number over the excitatory population of spikes within a burst, Ns, and the bursting frequency, f. The burst duration Td is defined as the average time between the first spike and the last spike in a burst.

Correlation calculation.

The average membrane potential of a neuron is

Formula 11(11)
The cross-correlation between the membrane potentials of excitatory neurons at positions x and x' is

Formula 12(12)
Note that we subtract the product of the average voltages.

Synchrony measure

For networks with all-to-all coupling, the synchrony measure {chi} of the excitatory population is defined as the fluctuation of the population-average membrane potential, normalized by the average fluctuations of the membrane potentials of single neurons (Golomb and Rinzel 1993Go, 1994Go; Golomb et al. 2001Go). We extend this definition for networks with spatially decaying synaptic coupling by calculating the population average over excitatory neurons in the neighborhood of a specific neuron in the middle of the chain, as explained in APPENDIX B. The measure {chi} varies between 0 (asynchronized state) and 1 (fully synchronized state).

Fast-slow analysis

We used the fast-slow method (Bertram et al. 1995Go; Hoppensteadt and Izhikevich 1997Go; Izhikevich 2000Go; Rinzel and Ermentrout 1998Go) to study the bursting mechanism of a single excitatory cell coupled to itself. This method has been applied successfully to analyze periodic bursting of various biophysical models of neurons and networks (e.g., Jung et al. 1996Go; Mandelblat et al. 2001Go; Tabak et al. 2000Go). The method makes use of basic dynamical systems theory (Strogatz 1994Go), which separates the variables of the system into two subsystems, "fast" and "slow," according to their overall timescales. In our system (the reduced model of an excitatory cell coupled to itself), the variables V (the membrane potential), h (inactivation of INa), n (activation of IKdr), and sAMPA (the fraction of open AMPA channels), with time scales on the order of 1 ms, are considered to be fast. The variables z (activation of IK-slow) and sNMDA (the fraction of open NMDA channels), with time scales on the order of 100 ms, are considered to be slow.2 In the first stage of the analysis, as will be described later, the variable sNMDA is considered to be constant, and the fast subsystem includes only the variable z. The bifurcation diagram of the fast subsystem is computed with the slow variable z considered as a parameter. Then the dynamics of the slow subsystem are computed using the time-averaged values of the fast subsystem. This is done by plotting two more curves on the bifurcation diagram. One is the activation curve of IM, Formula 12(Eq. A15). The second curve is the equivalent voltage during the oscillatory state, Vequiv (Bertram et al. 1995Go; Mandelblat et al. 2001Go), which is computed to determine the slow evolution of z in the full model. The value of Vequiv for a periodic cycle with a specific z is the voltage that, for a rest state (fixed point), yielded a value of dz/dt that is equal to the value of dz/dt averaged over the cycle. It was defined implicitly by the equation

Formula 13(13)
where V(t) is the voltage as a function of time along the cycle calculated for a specific value of z and Tz is the time period of the cycle that depends z. If, for a certain value of z, Vequiv(z) > z{infty}–1(z) (where z{infty}–1 (z) means the inverse function of z{infty}), than z in the full dynamical system (fast + slow together) increases slowly with a rate of order {tau}z. If Vequiv(z) < z{infty}–1(z), z decreases slowly.

The type of bursting we study in this paper belongs to the class of "square-wave bursters" (Rinzel and Ermentrout 1998Go), where the fast subsystem exhibits bistability of a rest state and a periodic firing state. During the rest state, the slow variable z decreases slowly, until the rest state does not exist anymore, and the fast subsystem jumps to the periodic firing state. During this state, the slow variable z increases slowly, until the firing state does not exist anymore and the fast subsystem jumps back to its rest state. This mechanism generates periodic bursting.

The bursting dynamics in our system is different from the classical square wave bursting mechanism because the slow variable sNMDA is not really a constant. We examine the effects of this second slow variable below.


 RESULTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We start our investigation of the persistent synchronized bursting states by considering disinhibited networks at Formula 13. In a second step, we study how the properties of the states vary when [Mg2+]o is elevated. Finally, we explore the role of inhibitory interneurons.

Synchronized persistent bursting for [Mg2+]o = 0

The disinhibited slice model with the reference parameter set Formula 13 exhibits persistent synchronized bursting activity in response to stimulation at the "left." The initial burst propagates as a traveling pulse at a constant velocity, similar to propagating single bursts in models of slices in which inhibition was reduced (Golomb and Amitai 1997Go) (Fig. 1A, left). Subsequent bursts propagate in a less orderly manner, especially near the edges of the slice model, and do not maintain a constant velocity of the left-to-right propagation. At long times after the initiation of the activity, cells burst with a frequency of Formula 13, and the bursting activity is synchronized, although not fully synchronized (Fig. 1A, right). Neighboring cells tend to burst at similar times (Fig. 1B, left), but spikes within bursts are only weakly synchronized, if at all (Fig. 1B, right). Summarizing, the model with reference parameter set generates persistent synchronized behavior with properties similar to those seen experimentally (Flint and Connors 1996Go; Silva et al. 1991Go; Sutor and Hablitz 1989Go; Tsau et al. 1998Go; Wu et al. 1999Go).


Figure 1
View larger version (32K):
[in this window]
[in a new window]
 
FIG. 1. Synchronized persistent waves occur in the disinhibited cortical slice model in response to a brief, local stimulation. The reference parameter set is used with Formula 5. A: rastergrams. Only firing times from every 32th cell are shown. Each spike is represented by a solid circle. Because the circles representing spikes within a burst are adjacent, each burst of spikes looks like a continuous line in the rastergram. Left: response to a local stimulation. Initially cells at left (x ≤ {sigma}/2) are depolarized and all others are at rest. Discharge propagates to the right. At each point, the periodic bursts follow the initial 1. Right: system activity at steady state. B: voltage time courses of 2 neighboring cells at Formula 5 and Formula 5, located in the middle of the slice model. The traces on the left, as well as the rastergrams, show that bursts of neighboring neurons are synchronized, although not fully. The same traces plotted at a magnified time scale (right) demonstrates that spikes within the bursts are synchronized very weakly, if at all. The vertical dotted lines are plotted at the times when the upper neuron is firing.

 
To further quantify the synchronization properties of the bursting activity, we compute the auto-correlation Cx,x({tau}) of the membrane potential V(t) of one neuron (Eq. 12) in the middle of the slice (Formula 13, Fig. 2A) and the cross-correlation Cx,x'({tau}) of this neuron with its neighbor at Formula 13(Fig. 2B). Both the auto- and cross-correlations oscillate on the bursting time scale, reflecting a periodic firing pattern and synchronization at the bursting level. The auto-correlation shows peaks at the spiking time scale for {tau} near zero (Fig. 2A, top right) but almost no peaks at the spiking time scale (and only one wide peak at the bursting time scale) for {tau} value around other integer multiplications of the time period Formula 13(Fig. 2A, bottom right). The cross-correlation shows almost no peaks at the spiking time scale around Formula 13(Fig. 2B, top right) and no peaks at the spiking time scale near other integer multiplications of Tper (Fig. 2B, bottom right). For all neurons within a footprint length {sigma} around the neuron at Formula 13, the amplitude of the cross-correlation is almost constant with the distance between the neurons, and the phase of the cross-correlation is zero. We conclude that the neurons are synchronized on the bursting time scale, at least within distances of order {sigma}, but not on the spiking time scale. Note that noise, heterogeneity, and sparseness are expected to cause the cross-correlations to decay with space and time (Golomb 1998Go).


Figure 2
View larger version (30K):
[in this window]
[in a new window]
 
FIG. 2. A: auto-correlation (Eq. 12) of the voltage time course of 1 neuron located at Formula 5. B: cross-correlation (Eq. 12) between the voltage time courses of the 2 neurons shown in Fig. 1B, located at Formula 5 and Formula 5. In both A and B, sections of the graphs within the dotted boxes around the middle and 1st right peaks are magnified on the right. Parameters: Formula 5.

 
The persistent activity in the model does not stop. In the experiments, however, bursting activity in disinhibited slices is terminated after ~3 s (Flint and Connors 1996Go). This termination can be explained by the effects of slow processes, such as slow inactivation of Na+ channels (Fleidervish and Gutnick 1996Go; Fleidervish et al. 1996Go), that are not included in the model.

Dependence of activity regimes on {tau}NMDA and gNMDA

The bursting activity of each neuron is terminated by the slow potassium current IK-slow, with activation kinetics constant Formula 13. Once this current has decayed, a strong enough slow NMDA-mediated excitation is responsible for the generation of a new burst. Therefore the decay time constant of the NMDA-mediated synapses, {tau}NMDA, is an important parameter for generating a persistent synchronized bursting state. We expect that this state can appear only if {tau}NMDA is large enough in comparison with {tau}z. Furthermore, one may expect that the {tau}NMDA value that allows persistent bursting should be large enough in comparison with Tper and that increasing gNMDA may allow persistent bursting even for smaller {tau}NMDA. To examine the roles of {tau}NMDA and gNMDA, we computed the regimes in the {tau}NMDAgNMDA plane where stimulation of the slice led to one of three states: persistent bursting activity, persistent tonic activity, or nothing but the quiescent state. The calculation was carried out without (Fig. 3A) and with (B) the inclusion of AMPA-mediated synaptic conductances (gAMPA = 0 and 0.08 mS/cm2 respectively). For both cases, bursting activity appear for large enough {tau}NMDA and moderate values of gNMDA. Interestingly, bursting activity is found even for {tau}NMDA values that are smaller than {tau}z, especially for Formula 13. For larger values of gNMDA and for {tau}NMDA that is not too small, the activity is tonic, whereas for small gNMDA or small {tau}NMDA, the persistent activity does not prevail. The gNMDA interval where bursting activity exists increases with {tau}NMDA and extends toward low-gNMDA values. Increasing gAMPA shifts the boundary of the bursting regime toward larger {tau}NMDA values. This means that for moderate values of {tau}NMDA, increasing the fast, AMPA-mediated excitation may abolish the persistent bursting behavior. This network behavior stems from the fact that increasing gAMPA increases the number of spikes within a burst (see following text), and therefore increases the activation of IK-slow that terminates the burst and prevents the next one. Similarly, the transition from the tonic regime to the bursting regime occurs for larger gNMDA values as gAMPA increases, because the stronger IK-slow due to the increased number of spikes terminate the spiking activity and generates a quiescent period before the next firing period begins.


Figure 3
View larger version (40K):
[in this window]
[in a new window]
 
FIG. 3. Regimes of activity as functions of {tau}NMDA and gNMDA. A and B: regimes obtained by simulating the 1-dimensional network with {rho} = 8,L = 16{sigma}. C and D: regimes obtained by simulating the reduced model of a cell coupled to itself. The value of gAMPA is 0 for (A and C) and 0.08 mS/cm2 for (B and D). In the regime denoted "quiescent" (white), the quiescent state is the only state at long times; it exists and is stable for all parameter values. A "bursting" state (dark gray) or a "tonic" state (light gray) coexists with a quiescent state. In the "tristable" regime (diagonal gray bands), the tonic, bursting and quiescent states are stable states. The tristable regime has a tongue-like shape. The numbers written in each tongue denote the number of spikes within a burst.

 
Reduced model: cell coupled to itself

The large synchronization on the bursting time scale among neighboring neurons in cases such as those shown in Figs. 1 and 2 allows the use of a simple approximation for the large dynamical system: a neuron coupled to itself with the same gNMDA and gAMPA values as in the network model (van Vreeswijk and Hansel 2001Go). This reduced model is amenable to mathematical analysis and is much faster to simulate to produce phase diagrams describing the bursting and tonic regimes as functions of various parameters. This simplified model is only an approximation, however, because the spikes within the bursts are not (or only weakly) synchronized. If the bursts themselves are less synchronized in the full network model, this approximation becomes even worse. Therefore we compare between the behavioral regimes of the reduced model (Fig. 3, C and D) and the full network model (Fig. 3, A and B). Qualitatively, the behavioral regimes of the two models are similar. Two main differences are found, however. First, for Formula 13, the bursting regime extends toward low {tau}NMDA in the full model in comparison with the reduced one. This effect is caused by the appearance of states with low synchronization in the full model as will be shown in the following text. Second, for Formula 13, there is a "tristable" regime where the bursting, tonic and quiescent states are stable states. Interestingly, the tristable regime has a tongue-like shape. Each tongue corresponds to a certain number of spikes in the burst, as written in Fig. 3D. Because the transition between a certain number of spikes to the subsequent one occurs through chaotic behavior (Terman 1992Go), it is not surprising that the tongue structure is complex. In our network simulations, the state is selected by the initial conditions, and therefore tristability is not seen. Overall, the reduced model can be used to describe the dynamical regimes of the whole network unless the network exhibits bursts that are not highly synchronized.

Dependence of bursting patterns on {tau}NMDA and gNMDA

We next examine how the main properties of the wave: the number of spikes within a burst Ns, the bursting frequency f, the burst duration Td, and the synchrony measure {chi} depend on the properties of NMDA synapses: gNMDA (Fig. 4, A and B) and {tau}NMDA (Fig. 4, C and D). The examination is carried out for parameters that lead to bursting, without (Formula 13) and with (Formula 13) AMPA excitation. Computation was performed for both the full network model (red) and the reduced model of one cell coupled to itself (black), to further compare their behavior. In the regimes where the two models yield persistent synchronized bursting, the bursting properties Ns, f, and Td, are, in general, very similar, showing that the reduced model mimics the full network model well. The number of spikes Ns is almost independent of gNMDA and {tau}NMDA. The bursting frequency f increases with gNMDA but only weakly with {tau}NMDA. To obtain physiological values of f, ~12 Hz (Flint and Connors 1996Go), gNMDA should not be too large. The burst duration Td increases weakly with gNMDA and is almost independent of {tau}NMDA. A main determinant of the burst duration is the time constant of the slow potassium current, {tau}z (see the section on the fast-slow analysis).

States with low synchrony

For Formula 13, the full model, but not the reduced model, shows bursting activity for low gNMDA values (Fig. 4A). These bursting states have low synchrony: {chi} is <0.2 in comparison to values >0.5 that are obtained for larger gNMDA values for which the bursts are highly synchronized. We find that for the value Formula 13, {chi} increases sharply. For this same gNMDA, f of the full model shows a small deflection downward, and the reduced model starts to fire bursts instead of being quiescent. Similarly, states with low bursting synchrony are found for Formula 13, the reference value of gNMDA and low {tau}NMDA, <60 ms (Fig. 4C). For Formula 13, such low-synchrony states do not appear, and all the bursting states are highly synchronized.

To explore the nature of this state further, we present a rastergram obtained in response to stimulation of the slice model with Formula 13 and Formula 13 in Fig. 5A. Initially, the pattern of activity (Fig. 5A, left) resembles that of the reference parameter set (Fig. 1A, left). The firing pattern in the steady state, however, is more complicated (Fig. 5A, right). The bursting time of neurons is widely distributed as exhibited both in the rastergram and in the voltage time courses of two bursting neurons (Fig. 5B). To determine whether the network is asynchronous, we computed the synchrony measure {chi} for a local population of 2{sigma}{rho} + 1 neurons, as a function of 1/{rho} where {rho} varies from 64 to 1,024 (Fig. 5C). The value of {chi} remains roughly constant near the value of 0.17, indicating that the state shows a low level of synchronization and is not asynchronized.

Dependence of activity regimes and bursting patterns on gAMPA and gNaP

The parameters gNMDA and {tau}NMDA control the slow dynamics of the system, on the order of 100 ms. Bursting behavior, however, depends also on the fast dynamics of the system (Bertram et al. 1995Go; Izhikevich 2000Go; Rinzel and Ermentrout 1998Go). In particular, fast, AMPA-mediated excitation was shown to induce bursting activity in networks of excitatory neurons with spike adaptation (van Vreeswijk and Hansel 2001Go). Even without AMPA excitation, our model neuron can generate endogenous bursting behavior in response to constant applied current, and the parameter gNaP represents the intrinsic conditional bursting property of the single neuron. Without gNaP, the neuron fires only tonically in response to a step current pulse. It responds by bursting to this current pulse for large enough gNaP, and the number of spikes within each burst increases with increasing gNaP (Golomb and Yaari, unpublished).

We study how persistent bursting in the network depends on gNaP and gAMPA by determining the various behavioral regimes as functions of these parameters using the reduced model. The phase diagram is plotted in Fig. 6 for two values of gNMDA: the reference value Formula 13(A) and a value that is 50% larger, Formula 13(B). At the lower gNMDA value, tonic firing is obtained for small values of gNaP and gAMPA (but not for Formula 13). Bursts are obtained for moderate gNaP values and gAMPA values that are not too large. The behavior for the higher gNMDA value, 0.105 mS/cm2 shares some similarity with the behavior for Formula 13 but differs in the following aspects. First, the tonic regime extends toward larger gNaP values. Second, the bursting regime extends toward larger gAMPA values. Third, bursts are generated even for Formula 13, though in a restricted regime. Fourth, a tristable regime emerges for moderate values of gAMPA and gNaP. For both values of gNMDA, increasing gAMPA or gNaP transfers the activity from the bursting to quiescent states. Reducing gAMPA may transfer the activity from bursting to tonic state but not to quiescence (except in the small areas near the tongues). Reducing gNaP, however, mostly transfers the activity from bursting to quiescent state. The intrinsic conductance gNaP and the synaptic conductance gAMPA display therefore some similarities, but also some differences, in their effects on the bursting activity.


Figure 6
View larger version (23K):
[in this window]
[in a new window]
 
FIG. 6. Regimes of activity of the reduced model of a cell coupled to itself as functions of gNaP and gAMPA. The parameter gNMDA is 0.07 mS/cm2 in A and 0.105 mS/cm2 in B. The color code is: bursting, dark gray; tonic, light gray; quiescent, white; tristable, diagonal gray bands.

 
The dependence of Ns, f, Td, and {chi} on gAMPA is shown in Fig. 7A for Formula 13. Increasing gAMPA causes Ns to increase and f to decrease. These two phenomena are related because an increase in Ns leads to an increase in the activation variables of IK-slow, resulting in an elongation of the inter-burst interval. The burst duration Td increases only weakly with gAMPA. Larger Ns values that result from increasing gAMPA indeed prolong the burst duration, but when Ns is kept constant, the burst duration even decreases with gAMPA. The bursts are synchronized ({chi} > 0.5), and {chi} gradually increases with gAMPA with many local peaks corresponding to minima in Td.

The effect of Ns and f on gNaP is similar to the effect of gAMPA, whether gAMPA is zero (Fig. 7B) or gAMPA = 0.08 mS/cm2 (Fig. 7C). The frequency f, however, depends only weakly on gNaP for Formula 13, and Td is almost constant with some fluctuations. For Formula 13, {chi} fluctuates and may reach low values representing states of low synchrony (Fig. 7B, bottom). Such low synchrony states do not exist for larger gAMPA (Fig. 7C, bottom). We conclude that the network exhibits persistent burst synchronization with substantial AMPA excitation. Without it, bursting states of low synchrony appear in many parameter regimes.

Fast-slow analysis of bursting activity

DESCRIPTION OF THE ANALYSIS.  The reduced model of one cell coupled to itself describes the properties of the network discharges very well (Fig. 3, 4, and 7). Therefore we use this simple model to analyze the properties of the persistent synchronized activity, using the fast-slow method (Bertram et al. 1995Go; Izhikevich 2000Go; Rinzel and Ermentrout 1998Go). We separate the variables of the system (the reduced model) into two subsystems, fast and slow, according to their overall time scales. The fast subsystem includes the variables V, h, n, and sAMPA, and the slow subsystem includes the variables z and sNMDA. We analyze the behavior of the single neuron coupled to itself for the reference parameter set, gNMDA = 0.07 mS/cm2, gAMPA = 0.08 mS/cm2 (Fig. 8, AC). To describe the effect of the two types of synapses on the burst properties, we carry out the analysis also for two more cases: AMPA excitation is blocked (gNMDA = 0.07 mS/cm2, gAMPA = 0; Fig. 8, DF) and NMDA excitation is enhanced (gNMDA = 0.105 mS/cm2, gAMPA = 0.08 mS/cm2; Fig. 8, GI). As a first step, we plot the voltage time course, V(t), and the time courses of the slow variables z(t) and sNMDA(t), computed by simulating the reduced model during one time period Tper, in Fig. 8, A, D, and G, for those three parameter sets.


Figure 8
View larger version (31K):
[in this window]
[in a new window]
 
FIG. 8. Fast-slow analysis of the model of 1 neuron coupled to itself. Synaptic conductances (in mS/cm2): A–C: gNMDA = 0.07, gAMPA = 0.08; DF: gNMDA = 0.07, gAMPA = 0; GI: gNMDA = 0.105, gAMPA = 0.08. A, D, and G, the time courses of V (top), z (middle), and sNMDA (bottom) in the full system (fast and slow) during one firing time period Tper are shown. B, E, and H, bifurcation diagram of the fast subsystem as a function of z. The value sNMDA is fixed at an average value it obtains during the cycle in the full system: 0.85, 0.92 and 0.93 for B, E, and H, respectively. The various line types are: thin solid line, stable rest; thin dotted line, unstable rest; thick solid line, stable firing (oscillatory) state; thick dotted line, unstable firing (oscillatory) state; green line, z{infty}(V) (Eq. A15); blue line, Vequiv (Eq. 13). HB, Hopf bifurcation (violet solid circle). C, F, and I: dynamical regime of the fast subsystem in the sNMDA-z plane. The black line denotes a saddle-node bifurcation of the fast subsystem, and the red line denotes the bifurcation in which the stable firing state ceases to exist (via a saddle-node of periodic bifurcation). Below the black line, the fast subsystem is in a firing mode. Above the red line, it is at rest. Between the 2 lines, there is bistability and both firing and rest are stable state. The green line denotes the projection of the trajectory of the full dynamical system on the sNMDA-z plane. The blue arrows point to the direction of this trajectory.

 
When a cell bursts, its state alternates between an active firing state and a silent rest state. The slow hyperpolarizing variable z and the slow depolarizing variable sNMDA increase during the firing phase of the burst and decrease during the silent phase. Their role in the burst generation, however, is opposite. As z increases, it terminates the firing phase of burst, and as it decreases, it may allow firing to resume again. In contrast, sNMDA delays the termination of the firing phase. If sNMDA decreases too much during the silent phase, the neuron does not resume firing. Furthermore, increasing {tau}NMDA >100 ms affects the bursting properties only weakly, namely the behavior for physiological value of {tau}NMDA is quite similar to the behavior for very large {tau}NMDA. Therefore as a first step, we study the system in the limit {tau}NMDA>> {tau}z and fix the slow variable sNMDA at a constant, average value it obtains during the cycle: 0.85, 0.92, and 0.93 for the three parameter sets (Fig. 8, A, D, and G, respectively, bottom) and consider the slow variable z to be a parameter. We compute the bifurcation diagrams of the fast subsystem as a function of z (Fig. 8, B, E, and H). The steady state (fixed point; thin black line) is stable for large z. This stable rest state coalesces with an unstable state and ceases to exist in a saddle-node bifurcation. The rest state, this time a high plateau, is stable again for small values of z. At the z value where the high rest state gains its stability (a Hopf bifurcation), an oscillatory state (limit cycle) emerges, corresponding to periodic firing. This oscillatory firing state extends toward the right. We plot the minimal and maximal values of the membrane potential during the firing state (thick black lines). There is a narrow regime in which the continuity of the stable oscillatory state is disturbed by an unstable oscillatory state, and for large values of z the amplitude of the oscillatory state that gains its stability is larger than the oscillation amplitude for smaller z. The oscillatory state coalesces with another unstable oscillator state (saddle-node of periodics bifurcation) that immediately terminates by intersecting the unstable fixed point (homoclinic, or saddle-loop, bifurcation).

Bistability of firing (oscillatory) and quiescent states in the dynamics of the fast subsystem is necessary, but not sufficient, for allowing bursting activity with one slow variable. To further determine the conditions for bursting in the full (fast and slow) model of a cell coupled to itself, we plot two more curves: the activation curve of the current IK-slow, z{infty}(V) (Eq. A15), plotted in green, and the equivalent voltage during the oscillatory state of the fast subsystem, Vequiv (Eq. 13), plotted in blue. For values of Vequiv above this green line, z slowly increases, and for values of Vequiv below this green line, z slowly decreases. For the three parameter sets we investigate (Fig. 8, B, E, and H), the solid (stable) blue line is above the green line, and therefore z increases during the firing state until it terminates the oscillations. The thin black line is below the green line, and therefore z decreases during the silent state until the rest state disappears (at the saddle-node bifurcation). The full dynamical system exhibits, therefore periodic bursting ("square-wave bursting") (Rinzel and Ermentrout 1998Go) for the three parameter sets shown in Fig. 8.

To complete the analysis, the two variables should be considered as slow variables with similar time scale. The activity regimes: rest, bistability, and firing are shown in the sNMDA-z plane (Fig. 8, C, F, and I). The boundaries of those regimes are diagonal lines, and their z coordinate increases with their sNMDA coordinate. The projection of the trajectory of the bursting state is shown in green. It oscillates through the bistable regime between the silent, rest state and the firing state. Both intrinsic neuronal properties (Fig. 8E) and AMPA-mediated excitation (Fig. 8B) contribute to the generation of this bistable regime. Using this analysis, we now explain the dynamical behavior of the model.

EXPLAINING THE EFFECTS OF VARYING gAMPA.  We first explain why Ns increases with gAMPA and f decreases with it, as shown in Fig. 7A. Comparison between the two values of gAMPA (Fig. 8, B and E) shows that the bistable regime is wider for larger gAMPA. The low rest states (both stable and unstable) are not affected by gAMPA because the AMPA activation curve is substantially larger than 0 only for very depolarized potentials (the dependencies of the rest voltages on z are slightly different for the two values of gAMPA because of the different values of the constant sNMDA). The conductance gAMPA, however, affects the periodic action potentials by increasing their frequency and reducing their amplitude. This happens because after the self-coupled neuron fires an action potential, the depolarizing conductance advances the firing of the subsequent one, the sodium current INa is still partially inactivated, and therefore the action potential overshoot is lower, the activation of IKdr is smaller, and the afterhyperpolarization that follows the action potential is less pronounced. Hence, for the same z, the minimal voltage during the action potential is more depolarized for larger gAMPA. The value of z, for which the minimal value of the action potential nearly coincides with the voltage of the unstable rest state, is larger and the bistable regime is wider when gAMPA increases. As a result, Ns is larger and f is smaller for larger gAMPA.

The intrinsic conductance gNaP also increases the bistable regime (not shown), and therefore Ns increases with gNaP (Fig. 7, B and C).

EXPLAINING THE EFFECTS OF VARYING {tau}NMDA.  We now turn to explain the dependence of Ns on {tau}NMDA and why activity stops at small {tau}NMDA values. The analysis becomes somewhat more complicated by the fact that the value of {tau}NMDA is not much larger than {tau}z. If the kinetics of sNMDA were very slow, the direction of the projection of the neuronal dynamics onto this plane would be parallel to the ordinate (z-axis; Fig. 8F). Because these kinetics are not much slower than that of z, the value of both z and sNMDA increase during the firing phase and decrease during the silent phase. Because the boundaries of the bistable regime are diagonal, the z interval that the neuronal trajectory covers during the bursting cycle is larger for faster NMDA kinetics. As a result, the number of spikes decreases somewhat with {tau}NMDA, for small {tau}NMDA values, as shown in Fig. 4, C and D.

If {tau}NMDA becomes smaller in comparison to {tau}z, the trajectory of the dynamical system in the sNMDA-z plane becomes more and more parallel to the abscissa. After firing several spikes, the system is stuck at the rest state (on the left), the bursting process is terminated, sNMDA decreases to 0 and z decreases to a very small value. This explains why persistent bursts cannot occur for too small {tau}NMDA (Fig. 3).

EXPLAINING THE EFFECTS OF VARYING gNMDA.  The numerical simulations (Figs. 4, A and B, and 3) show that, first, in the persistent bursts regime, f increases with gNMDA, whereas Ns and Td are almost independent of gNMDA; second, for large enough gNMDA, the activity becomes tonic instead of bursting. To explain these two numerical observations, we note that the slow currents IK-slow and INMDA for the self-coupled neuron at Formula 13 are (Eqs. A13 and A20)

Formula 14(14)

Formula 15(15)
If z and sNMDA are considered constants, varying gK-slow x z or gNMDA x sNMDA is roughly like varying the external current Iapp. There is, however, a difference that stems from the driving force term VVrev, where Vrev, the reversal potential of K+ or glutamate. Because of this driving force, increasing z is more powerful for depolarizing potentials, whereas increasing sNMDA is more powerful for hyperpolarizing potentials. This means that when gNMDA x sNMDA increases, the bifurcation diagram should move to the right (toward larger z values), and it does it more for hyperpolarizing currents to compensate for the smaller driving force. This is why both the steady-state curve and the firing curves are shifted to the right in Fig. 8H in comparison with Fig. 8B, but the steady-state curve is shifted more to the right and the bistable regime is smaller.

If gNMDA is large enough, the solid blue curve, representing the effective potential Vequiv for the firing state, is shifted enough to the right such that it intersects with the green curve. In this case, the firing state is stable (in the full system) and the system does not burst.

If gNMDA is large, but still below the value for the tonic behavior of the full system, the neuron still bursts but the bistable regime is narrower. The value of dz/dt in the full system depends on [z{infty}(V) – z)/{tau}z (Eq. A14), namely on the distance between the thin solid black line and the green line for rest states, or the solid blue line and the green line for the tonic states. As shown in Fig. 8H, the impact of the shift to the right of the steady state curve is weak, because the value z{infty}(V) – z remains about 0.1. The bistable regime is narrower, and therefore the silent phase of the burst decreases as gNMDA increases. The shift of the solid blue curve to the right is much more important because the solid blue curve and the green curve almost intersect. Therefore the value of dz/dt during the firing phase of the burst is smaller for larger gNMDA. This compensates for the decrease in the z regime where the system is bistable, and therefore the burst duration Td, as well as Ns, are almost independent of gNMDA.

Effects of nonzero [Mg2+]o

Raising [Mg2+]o leads to several physiological effect. The main effect is on the response of NMDA receptors. Under physiological conditions, when the membrane potential is close to its resting values, NMDA receptors are blocked by Mg2+ ions. Increasing [Mg2+]o corresponds to shifting the sigmoid function fNMDA(V), the voltage dependency of the NMDA synaptic conductance on the postsynaptic voltage, toward the right. This means that {theta}NMDA, the half-maximum value of fNMDA(V), is depolarized as [Mg2+]o increases from –{infty} for Formula 15 to around –30 mV for physiological value (Jahr and Stevens 1990Go; Kumar and Huguenard 2003Go). In addition, elevating [Mg2+]o, like elevating the concentration of other divalent cations, decreases synaptic transmission (Johnston and Wu 1995Go) and increase the action potential threshold (Frankenhaeuser and Hodgkin 1957Go). We analyze the consequences of varying {theta}NMDA and then address the two other effects.

EFFECTS OF VARYING {theta}NMDA.  We computed the dynamical regimes in the {theta}NMDA-gNMDA plane for Formula 15(Fig. 9A) and for Formula 15(Fig. 9B). The bursting regime has a "backward L" shape. Whereas at hyperpolarized {theta}NMDA values, it gradually shifts to larger gNMDA values, as {theta}NMDA increases, it turns toward very large gNMDA values for {theta}NMDA above –60 mV. The firing patterns at more depolarized {theta}NMDA and large gNMDA are characterized by bursting with a large number of fast spikes in a burst (see traces on the right of Fig. 9, A and B). No bursting is obtained for gNMDA values that are not extremely large and {theta}NMDA above –40 mV for Formula 15 or –55 mV for Formula 15. Thus the model predicts that there is no persistent synchronized activity under physiological conditions for Formula 15(Kumar and Huguenard 2003Go). Moreover, raising gAMPA is expected to shift the bursting regime toward more negative {theta}NMDA values.


Figure 9
View larger version (30K):
[in this window]
[in a new window]
 
FIG. 9. A and B: regimes of activity as functions of {theta}NMDA and gNMDA. Regimes obtained by simulating the 1-dimensional network with Formula 5. The color code is: bursting, dark gray; tonic, light gray; quiescent, white. The value of gAMPA is 0 for A and 0.08 mS/cm2 for B. The 2 voltage time courses on the right of each panel were simulated for the parameters denoted by I and II in the diagram on the left. The dotted line denotes –90 mV. C: dependence of the number of spikes within a burst Ns, the bursting frequency f, the burst duration Td, and the synchrony measure {chi} on {theta}NMDA. Calculations are carried out for the network model with gAMPA = 0.08 mS/cm2, gNMDA = 0.105 mS/cm2, {rho} = 64, L = 16{sigma}.

 
The dependence of Ns, f, Td, and {chi} on {theta}NMDA is shown in Fig. 9C for Formula 15. The number of spikes Ns increases and the frequency f decreases with {theta}NMDA. This means that, for nonzero [Mg2+]o values, neurons fire low-frequency bursts with large number of spikes. The burst duration Td depends only weakly on {theta}NMDA, and the system is in a highly synchronized state ({chi} above 0.5).

EXPLAINING THE EFFECTS OF VARYING {theta}NMDA.  We now explain, using the fast-slow analysis, the dependence of Ns, f, Td, and {chi}. We perform this analysis of the reduced model of one cell coupled to itself for Formula 15, and Formula 15(Fig. 10) and compare this parameter set with a similar one with Formula 15(Fig. 7, G–I). Because the frequency f is lower for more depolarized {theta}NMDA, sNMDA decreases to lower values (Fig. 10A, bottom), and hence the bifurcation diagram (Fig. 10B) is computed for an average value of Formula 15. The structure of the bifurcationdiagram for Formula 15 resembles that for {theta}NMDA = –{infty}, but there are important quantitative differences. Because raising {theta}NMDA has a stronger impact on the NMDA current at hyperpolarized voltages (the effect of Mg2+ blockade on NMDA-mediated excitation is more effective at hyperpolarized potentials), the bifurcation diagram is shifted more to the left at more negative values of V. This means that both the rest state and the firing state extend to lower values of z, but the rest state extends more, and therefore the bistable regime is wider for more depolarized {theta}NMDA. In addition, the rest state potential is more hyperpolarized because the depolarizing effect of the NMDA excitation is weaker because of the Mg2+ blockade. The widening of the bistable state and its shift toward lower values of z are not just a consequence of the smaller value of sNMDA for more depolarized {theta}NMDA, as shown in Fig. 10C (compare with Fig. 7H). Because the bistable regime is wider and the rest state is more hyperpolarized, the silent phase of the burst is more prolonged and reaches lower values of V. Hence, f decreases with {theta}NMDA. The neuron fires more spikes with {theta}NMDA because the bistable regime extends toward the z values for which the firing state has lower amplitude and higher firing frequencies. In contrast to the duration of the silent phase of the burst, the duration of the firing state Td remains more or less the same despite the widening of the bistable regime. As in the case of increasing gNMDA, this is a result of the fact that shifting the firing state to the right increases the distance between the blue and green lines, and therefore the value of dz/dt during the firing phase (Eq. 13) is larger for {theta}NMDA = –70 (Fig. 10B) than for Formula 15(Fig. 7H). The dynamics of z during the firing phase of the burst is therefore faster and compensates for the widening of the bistable regime, such that Td remains almost constant.


Figure 10
View larger version (15K):
[in this window]
[in a new window]
 
FIG. 10. Fast-slow analysis of the model of 1 neuron coupled to itself for non-0 [Mg2+]o. Parameters: gNMDA = 0.105 mS/cm2, gAMPA = 0.08 mS/cm2, {theta}NMDA = –70 mV. A: the time courses of V (top), z (middle), and sNMDA (bottom) in the full system (fast and slow) during one time period Tper. B: bifurcation diagram of the fast subsystem as a function of z. The value sNMDA is fixed at 0.76, the average value it obtains during the cycle in the full system. The various line types are the same as in Fig. 8B. HB, Hopf bifurcation. C: dynamical regime of the fast subsystem in the sNMDA-z plane. Symbols are as in Fig. 8C.

 
DIVALENT CATION EFFECTS.  Elevating [Mg2+]o increases firing threshold. We mimic this effect phenomenologically in the model by decreasing the rest potential of the neuron, namely by hyperpolarizing VL. This causes a small increase in Ns and a decrease in f; Td, and {chi} are almost unchanged (Fig. 11). If VL is too depolarized, the persistent activity ceases to exist. Hence, hyperpolarizing VL affects the activity in a similar way to depolarizing {theta}NMDA.


Figure 11
View larger version (16K):
[in this window]
[in a new window]
 
FIG. 11. Dependence of the number of spikes within a burst Ns, the bursting frequency f, the burst duration Td, and the synchrony measure {chi} on VL. Calculations are carried out for the network model with gAMPA = 0.08 mS/cm2, gNMDA = 0.07 mS/cm2, {rho} = 64, L = 16{sigma}.

 
Decreasing synaptic transmission as a result of elevating [Mg2+]o leads to a decrease in gAMPA and gNMDA. Decreasing gNMDA causes f to decrease, almost does not affect Ns, and eventually terminates the activity (Fig. 4, A and B). Its effect is therefore similar to the effect of depolarizing {theta}NMDA. Decreasing gAMPA has an opposite effect: it causes f to increase and Ns to decrease (Fig. 7, A and B) and extends the regime in which persistent bursting exist toward larger {theta}NMDA values (Fig. 9A). Note that persist bursting activity cannot occur for {theta}NMDA = –31 mV even for gAMPA = 0. Near the parameter regime in which bursting disappear, the bursting frequency is large even for Formula 15.

Networks with inhibition

Most of the experimental recordings in cortical slices with low [Mg2+]o were carried out with inhibition intact. We therefore examine the effect of including an inhibitory population on the dynamics of the persistent synchronized bursting. An example is shown in Fig. 12, A and B, where the excitatory cells with the reference parameter set (Fig. 1) are inhibited by GABAA-mediated receptors with conductance strength Formula 15. The inhibitory cells are excited by AMPA-mediated receptors with Formula 15. Rastergrams of the excitatory and inhibitory neurons are shown in Fig. 12A, and voltage time traces of one excitatory neuron and one inhibitory neuron in the middle of the slice model are shown in Fig. 12B. The firing pattern of the excitatory cells is similar to the firing pattern without inhibition (compare Fig. 12, A and B, with Fig. 1). The inhibitory cells burst as well, and their bursting times follow those of the adjacent excitatory cells. The firing patterns of the inhibitory cells, but not those of excitatory cells, are different when the inhibitory cells are also excited by NMDA-mediated receptors with Formula 15= 0.05 mS/cm2 (Fig. 12, C and D). The sustained NMDA-mediated excitation causes the inhibitory cell to fire even during the silent periods of the excitatory cells, although the firing rate of the inhibitory cells is enhanced during the firing period of the excitatory cells in response to the fast AMPA-mediated excitation. The number of spikes Ns fired by excitatory cells during a burst decreases with increasing the excitatory-to-inhibitory synaptic conductances (Fig. 13, A and B) or the inhibitory-to-excitatory synaptic conductances (C and D) because the recurrent inhibition contributes to the burst termination. This decrease in Ns causes an increase in f when Formula 15 because the slow K+ current is less activated (Fig. 13, A and C). When gNMDAEI is large enough and the inhibitory neurons fire during the whole cycle, they hyperpolarize the excitatory cells and reduce f (Fig. 13, B and D). To conclude, our modeling results show that GABAA-mediated inhibition causes the f to increases gradually if the inhibitory neurons are excited by AMPA-mediated synaptic conductance only and to decrease gradually if the inhibitory neurons are also excited by strong enough NMDA-mediated synaptic conductance.


Figure 12
View larger version (52K):
[in this window]
[in a new window]
 
FIG. 12. Synchronized persistent waves occur in a cortical slice model, composed of excitatory and inhibitory populations, in response to a brief, local stimulation. The reference parameter set is used with gAMPA = 0.08 mS/cm2, gNMDA = 0.07 mS/cm2, {rho} = 32, L = 32{sigma}, n = 1,024, gGABAAIE = 0.05 mS/cm2, gAMPAEI = 0.2 mS/cm2. The value gNMDAEI is 0 in A and B and 0.05 mS/cm2 in C and D. A and C: rastergrams of excitatory cells (top) and inhibitory cells (bottom). Only firing times from every 32th cell are shown. Each spike is represented by a solid circle. Because the circles representing spikes within a burst are adjacent, bursts of spikes of excitatory cells often look like continuous lines in the rastergram. C and D: voltage time courses of an excitatory cell (top) and an inhibitory cell (bottom) located in the middle of the slice model (Formula 5).

 

Figure 13
View larger version (26K):
[in this window]
[in a new window]
 
FIG. 13. Dependence of the number of spikes within a burst Ns and the bursting frequency (f) on gAMPAEI for Formula 5(A and B) and on gGABAAIE for Formula 5(C and D). Calculations are carried out for Formula 5(A and C) and for Formula 5gAMPAEI (B and D) and are based on simulations of the network model with Formula 5Formula 5.

 

 DISCUSSION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We have shown that a network of excitatory neurons coupled by NMDA-mediated synapses can show persistent synchronized bursting activity. Our model proposes two main physiologically relevant conditions for the emergence of such activity: 1) the presence of enough neurons with intrinsic bursting properties in the population, or enough AMPA excitation; and 2) relatively open NMDA receptors, such that may exist if {theta}NMDA is hyperpolarized or the neuronal population is depolarized enough. This persistent activity is robust for physiological values of the {tau}NMDA (decay time constant of NMDA-mediated synapses) of ~100 ms. In addition, the conductance gNMDA should be moderate; large gNMDA supports tonic firing. AMPA conductance often shifts tonic to bursting activity. Surprisingly, gAMPA that is too large abolishes all persistent activity in the model. Bursts are often highly synchronized, whereas the spikes within the burst are not (Pinsky and Rinzel 1994Go; van Vreeswijk and Hansel 2001Go). Without AMPA conductance, the system may settle into a weakly synchronized state. Elevating [Mg2+]o decreases the bursting frequency and eventually terminates the activity. Even strong increases in gNMDA cannot generate synchronized persistent activity in the cortical slice for physiological [Mg2+]o. If inhibitory cells do not receive NMDA-mediated excitation, inhibition increases f; if they receive strong NMDA-mediated excitation, inhibition decreases f.

Consequences of the specific model

We have analyzed a bursting neuron model from the "square-wave" type (Bertram et al. 1995Go; Izhikevich 2000Go; Rinzel and Ermentrout 1998Go), which is based on bistability of the fast subsystem of variables. By applying the fast-slow analysis, we have shown that the results are a consequence of the bursting mechanism and not just of the specific choice of intrinsic and synaptic parameters. We expect most of the main qualitative results of this article pertaining to the persistent synchronized bursts to be valid also for other bursting mechanisms. For every burster, {tau}NMDA should be large enough to overcome the hyperpolarizing current(s) that terminate a burst and start a new one. The exact minimal value of {tau}NMDA is parameter-dependent, but roughly it should be at least around the time period, 100 ms. This value is similar to the decay rate of NMDA-mediated EPSCs in intracortical synapses (Kumar and Huguenard 2003Go). If gNMDA is too large, the activity is not terminated and neurons fire tonically. Increasing AMPA is expected to increase the number of spikes within a burst and therefore terminate the activity completely in models where a burst is terminated by activation of hyperpolarizing currents.

One necessary condition for this model to generate bursts, however, may not be needed in other bursting mechanisms. In our model, the kinetics of the spike-generating currents INa and IKdr are fast enough to allow bistability in the model of a cell coupled to itself (Bertram et al. 1995Go; Rinzel and Ermentrout 1998Go). In other bursting models, for example, those currents based on ping-ponging between soma and dendrite (Mainen and Sejnowski 1996Go; Pinsky and Rinzel 1994Go) may have slower kinetics.

Roles of gAMPA and gNaP

The parameter gNaP controls the tendency of the isolated neuron to burst in response to external depolarizing current. According to our model, for small or moderate values of gAMPA and moderate values of gNMDA, moderate values of gNaP are needed to obtain network bursting. This result is consistent with the experimental observation that drugs that suppress INaP abolish the 10-Hz oscillations in the agranular neocortex (Castro-Alamancos and Rigas 2004Go).

The AMPA conductance gAMPA may support persistent bursting behavior even when Formula 15, but very large gAMPA values terminate the activity (Fig. 6). This conductance shifts the bursting regime toward larger gNMDA values (Fig. 3). In addition, enhancing gAMPA increases the number of spikes within a burst, decreases the bursting frequency, enhances bursting synchronization, and prevents the appearance of weakly synchronized bursting states (Fig. 7).

In the reduced model of a cell coupled to itself, the synaptic AMPA conductance and the persistent Na+ (NaP) conductance have a similar, although somewhat different, properties. Their kinetics are fast, but that of the NaP conductance is much faster and is considered here to be instantaneous. They are activated at depolarized membrane potential, but the NaP conductance is activated above Formula 15, whereas the AMPA conductance is activated above Formula 15(namely, when the presynaptic neuron fires an action potential). As a result, the dynamical roles of the two conductances, while similar, have some differences. Increasing gNaP, like increasing gAMPA, extends the bistable regime in the bifurcation diagram of the fast subsystem and therefore causes an increases in the number of spikes within a burst and, in general, a reduction of bursting frequency. Because the NaP conductances are activated above –47 mV, increasing gNaP, but not gAMPA, can support a high plateau or fast spiking riding on a high plateau. Decreasing gNaP, but not gAMPA, can transfer a bursting state into a quiescent state (Fig. 6). This is a result of the fact that the NaP conductance is partially activated in subthreshold voltages and can support the firing of a new spike, whereas the AMPA conductance is activated much above firing threshold.

The reduced model cannot, of course, address the issue of synchronization in large networks. Indeed, the AMPA conductance supports synchronization and prevents weakly synchronized states, whereas such states may exist with substantial gNaP values.

Contribution to the theory of synchronization and comparison with other models

Most theoretical works on synchronization deal with synchrony of spikes (reviewed in Golomb et al. 2001Go). Only a few models deal with synchronization of bursts of spikes, and in most of them, the activity is not persistent, namely there is no stable rest state that coexists with the firing state. Our model differs from those models by creating an active state that coexists with a quiescent state. Still, it is possible to compare the bursting mechanism of our model with the bursting mechanism in those models by noticing that the applied current there plays a similar role to the role of NMDA synaptic current here, which is roughly constant over the duration of a burst cycle. Previous work (van Vreeswijk and Hansel 2001Go) describes bursting activity that emerges from an excitatory network of neurons as a result of adaptation in neurons that are firing tonically when isolated. Similarly to our model, bistability of the fast subsystem (without the activation variable of the adaptation current, which is similar to z) is needed for generating bursting activity. Such bistability is generated in the model of van Vreeswijk and Hansel (2001)Go by the effect gAMPA only, whereas in our model, it can be generated by gAMPA or gNaP. The active phase of the burst in the conductance-based model of van Vreeswijk and Hansel (2001)Go is characterized by the strong undershoot of the membrane potential, which is more hyperpolarized than the membrane potential during the silent phase of the burst (Fig. 10B there). In our model, the minimal membrane potential during the active phase of the burst is more depolarized than the membrane potential during the silent phase of the burst (Figs. 1, 5, 8, and 9) as found in intracellular recording from cortical neurons during the NMDA-induced activity (Kawaguchi 2001Go; Sutor and Hablitz 1989Go).

Like our model and the model of van Vreeswijk and Hansel (2001)Go, highly, but not fully, synchronized bursts can be obtained in networks of spontaneous bursters (Pinsky and Rinzel 1994Go). In addition, we find states with low, but not zero, synchrony that occur with slow excitation only.

A model of epileptiform activity induced by low Mg2+ in the rat hippocampal slices (Traub et al. 1994Go) showed bursting activity for hundreds of milliseconds. NMDA receptors in that model decayed very slowly, with the time scale of desensitization (Fig. 5 in Traub et al. 1994Go). In our work, we study the roles of NMDA kinetics, synaptic and intrinsic strengths in determining bursting mechanism, and level of synchrony, issues that were not addressed there. Persistent tonic activity was found in a model of NMDA-coupled neurons without Mg2+ blockade (Ermentrout 2003Go). That model does not show bursting because the neurons lacked AMPA excitation and the necessary intrinsic properties, for example, slow adaptation.

Bursting activity can be generated in networks of excitatory and inhibitory neurons without adaptation (Hansel and Mato 2003Go) as a result of destabilization of an asynchronous firing state. Those synchronized bursts are not persistent. It can be shown, however, that the persistent synchronized bursting states can occur in the model of Hansel and Mato (2003)Go if excitatory synapse are slow enough (D. Hansel, private communication). The inhibitory population in that network model takes the role of adaptation in our model in terminating the firing phase of the burst. Unlike our model, such a model will not be consistent with the experimental observation that persistent synchronized states occur in disinhibited networks.

Slow NMDA-mediated synapses without Mg2+ blockade were shown to be important for generating persistent asynchronized activity in excitatory networks (Compte et al. 2000Go; Wang 1999Go) but not necessarily when an inhibitory population is added to the model (Hansel and Mato 2001Go). Here we have shown that when neurons possess an adaptation current and the have the proper intrinsic and/or synaptic properties (namely, bistability of the fast subsystem), NMDA-mediated synapses can lead to synchronized bursting activity. Mg2+ blockade of NMDA synapses prevents the appearance of persistent bursting for physiological [Mg2+]o. It will be interesting to systematically examine the effects of [Mg2+]o on persistent asynchronized states.

Endogenously active cells may spontaneously initiate a bursting episode (Latham et al. 2000Go). When a neuronal system shows bistability of a bursting state and a quiescent state, as described here, an episode of bursting activity can also be evoked by an external stimulus. Bursting activity, without bistability with a quiescent state, is observed in cultures of cortical neurons receiving stochastic input (Giugliano et al. 2004Go).

Comparison with experimental results.

The bursting frequency f in disinhibited slices is ~10–14 Hz (Flint and Connors 1996Go). Our model predicts that f cannot be much less than that because in that case the NMDA conductance would decay and nothing would start a new burst after the previous one has terminated. The results of our model (Figs. 3 and 4) show that decreasing f (or, equivalently, decreasing {tau}NMDA) cannot be compensated by increasing gNMDA because tonic, and not bursting, activity is obtained for large gNMDA. With AMPA-mediated synapses blocked and GABAA inhibition intact, Flint and Connors (1996)Go recorded only two transient population bursts within a time difference of ~300 ms. This activity is classified as transient and not as persistent according to the classification of this paper because it does not last for many cycles.

Local field potential recording and optical imaging reveal that the bursting activity is synchronized locally (Silva et al. 1991Go) and that phase shifts can develop as the distance between neurons grows (Wu et al. 1999Go). Similar behavior is obtained in our model (Figs. 1 and 2). The phase shifts are expected to grow if the network is heterogeneous. The fact that the activity in our model is persistent and does not result from entrainment by a group of pacemakers is consistent with the optical imaging results that the phase relationship between starting times of bursts vary with time (Wu et al. 1999Go).

Our modeling results indicate that inhibitory neurons fire bursts in synchrony with their neighboring excitatory cells if they receive mainly AMPA-mediated excitation; they fire almost continuously if they receive strong enough NMDA-mediated excitation. Intracellular recordings from fast-spiking inhibitory interneurons show examples for the two firing patterns (Figs. 4D and 8A in Kawaguchi 2001Go). Even inhibitory cells that fire continuously tend to fire more spikes during the bursting period of the excitatory cells (Fig. 8 in Kawaguchi 2001Go), in agreement with our simulations (Fig. 12). Our theory suggests that blocking GABAA-mediated inhibition decreases the frequency f of the network bursting if inhibitory neurons are excited mainly by AMPA-mediated receptors, but this blockade increases f if inhibitory neurons are substantially excited by NMDA-mediated receptors (Fig. 13). Because experiments carried out in the somatosensory cortex revealed that the bursting frequency increases weakly with blocking inhibition (Flint and Connors 1996Go), we predict that inhibitory cells in this cortical area should receive NMDA-mediated excitation from neighboring excitatory cells.

Our model takes into account kinetic processes that are faster than ~100 ms. Slow processes, such as slow inactivation of Na+ currents (Fleidervish and Gutnick 1996Go; Fleidervish et al. 1996Go) or slow synaptic depression, are neglected. As a result, the synchronized bursting activity is persistent and does not stop. In the experimental systems, those slow processes and others finally terminate the activity after <1 s if inhibition is intact and after ~3 s if inhibition is blocked (Flint and Connors 1996Go). The epileptic-like activity generated in response to lowering [Mg2+]o leads also to accumulation of extracellular Ca2+ and K+ ion and the alternation of the K+-Cl electrogenic pump. Such slow processes are not considered here.

Short-term synaptic processes, depression and facilitation, are neglected here for simplicity (Golomb and Amitai 1997Go; Markram and Tsodyks 1996Go). Because NMDA-mediated receptors in this model are close to saturation after a brief burst of a few spikes, the effects of depression and facilitation on the NMDA-mediated excitation are expected to be weak. Depression effects on AMPA-mediated excitation can be mostly compensated for by increasing gAMPA. The residual effect may lead to a moderate decrease in Ns.

The present paper describes a deterministic and ordered network. Noise, heterogeneity, and sparseness are expected to cause the cross-correlations to decay with space and time and to eliminate any synchrony at the spiking time scale (Fig. 2). Only large levels of sparseness are expected to desynchronize bursts of neighboring neurons (Golomb 1998Go).

The most pronounced effect of elevating [Mg2+]o is a shift of the activation curve fNMDA(V) (Eq. 8) toward more depolarized values, namely making {theta}NMDA more depolarized. As a result, f decrease and Ns increases, and for large enough {theta}NMDA the persistent activity cannot exist (Fig. 9). Elevating [Mg2+]o has several other effects, although the strength of these effects in cortical neurons is not known quantitatively. It causes an increase in firing threshold and decrease of synaptic transmission. The effects of raising firing threshold and reducing NMDA-mediated synaptic conductance enhance the effects of depolarizing {theta}NMDA. Reducing AMPA-mediated synaptic conductance reduces the effects of depolarizing {theta}NMDA on f but cannot prevent the termination of activity for depolarizing enough {theta}NMDA. Near the {theta}NMDA value where the persistent activity ceases to exist, f is low even when gAMPA is reduced (Fig. 9).

Prediction from the model

Our analysis yields several predictions that can be tested in experimental systems such as cortical slices. First, gradual blockade of AMPA-mediated excitation is expected to decrease the number of spikes in a burst Ns and increases the bursting frequency f, or even transfer a synchronized bursting state into a tonic, asynchronous state. It will be interesting to look for weakly synchronized states when gAMPA is blocked. Second, Ns is decreased when the persistent Na+ current is blocked. Third, gradual blockade of NMDA-mediated excitation is expected to decrease f and but not to change Ns considerably. Fourth, raising [Mg2+]o gradually will decrease f and increase Ns, until eventually the activity will stop. Fifth, the excitatory-to-inhibitory conductance in the somatosensory cortex includes an NMDA-mediated component.

Functional implications

This work shows that slow NMDA synapses can support persistent synchronized bursts as long as they are not highly blocked by Mg2+ ions. This mechanism is relevant for epileptiform activity (Telfeian and Connors 1999Go). It may also be relevant for cortical tissues in vivo provided that the Mg2+ blockade is weak, namely {theta}NMDA is hyperpolarized (Fleidervish et al. 1998Go) or the neuron is depolarized, such as during the "up state" (McCormick 2005Go; Steriade et al. 1993Go). Models of such networks with weak Mg2+ blockade of NMDA-mediated excitation were shown to exhibit persistent asynchronized states (Compte et al. 2000Go). We show that under the certain intrinsic and synaptic conditions we define, persistent synchronized bursting state states (Cardoso de Oliveira et al. 2001Go) can also occur in such networks.

Intrinsic bursters in cortex appear in layer V, and the low [Mg2+]o waves propagate in layer V alone but not in cortical slices without layer V (Flint and Connors 1996Go; Silva et al. 1991Go). The model shows that the intrinsic bursting property is important (although not necessary) for obtaining robust persistent network bursting activity at low [Mg2+]o. Therefore the cellular properties of layer V neurons together with the results of this model, taken together, can explain why this type of epileptiform propagates only in layer V. Other factors, such as the existence of long-range cortico-cortical connections in layer V (Connors and Amitai 1995Go), can also support those waves.

NMDA receptors underlying intracortical connections in young rats (P13–P21) are composed of the NR2B subunit, and the decay rate of their EPSCs is ~108 ms (Kumar and Huguenard 2003Go). The decay rate of EPSCs of the NR2A subunits is faster, ~65 ms. The NR2A subunit becomes abundant in cortex with respect to the NR2B subunit as a result of development (Quinlan et al. 1999Go) or rule learning (Quinlan et al. 2004Go). Our model predicts that NMDA-mediated persistent synchronized activity at ~10 Hz can barely exist if most of the NMDA receptors are composed of the fast NR2A subunits. Faster synchronized activity, however, such that is found in vivo (Cardoso de Oliveira et al. 2001Go), can be supported even by NMDA receptors with NR2A subunits.


 APPENDIX A
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Model equations and parameters

There are N excitatory cells and N inhibitory cells, representing a spatial discretization of the continuous integro-differential equations (Eqs. 1 and 2). The position of the ith neuron is xi = i/{rho}. We specify here our reference parameter set, that is used throughout the paper unless stated otherwise.

Excitatory cells

We use the Hodgkin-Huxley-like formulation for the excitatory cortical cells based on a previous model (Golomb and Amitai 1997Go).

Current balance equation


Formula A1(A1)
Formula A1

Intrinsic Currents.

Sodium current INa 

Formula A2(A2)

Formula A3(A3)

Formula A4(A4)

Formula A5(A5)

Formula A6(A6)
gNa = 35 mS/cm2, VNa = 55 mV, {theta}m = – 30 mV, {sigma}m = 9.5 mV, {theta}h =–45 mV, {sigma}h = –7 mV, {theta}th = –40.5 mV, {sigma}th = –6 mV.

Persistent sodium current INaP 

Formula A7(A7)

Formula A8(A8)
Formula A8.

Delayed rectifier potassium current IKdr 

Formula A9(A9)

Formula A10(A10)

Formula A11(A11)

Formula A12(A12)
gKdr = 3 mS/cm2, VK = –90 mV, {theta}n = –33 mV, {sigma}n = 10 mV, {theta}tn = –27 mV,{sigma}tn = –15 mV. The kinetics INa and IKdr here are faster than those of Golomb and Amitai (1997)Go.

Slow potassium current IK-slow 

Formula A13(A13)

Formula A14(A14)

Formula A15(A15)
Formula A15.

Leak current IL 

Formula A16(A16)
Formula A16.

Synaptic Currents.

AMPA current IAMPA 

Formula A17(A17)

Formula A18(A18)

Formula A19(A19)
Formula A19 gAMPA=Formula A19.

NMDA current INMDA 

Formula A20(A20)

Formula A21(A21)

Formula A22(A22)

Formula A23(A23)
Formula A23Formula A23(Kumar and Huguenard 2003Go). The value of {theta}NMDA is –{infty} for Formula A23, and increases logarithmically with [Mg2+]o (Jahr and Stevens 1990Go).

GABAA current IGABAAIE


Formula A24(A24)
Formula A24(for disinhibited networks) or 0.05 mS/cm2 (for networks with inhibition), Formula A24.

Inhibitory cells

We use the Wang-Buzsáki model (Wang and Buzsáki 1996Go).

Current Balance Equation


Formula A25(A25)
Formula A25.

Intrinsic Currents.

Sodium current INaI 

Formula A26(A26)

Formula A27(A27)

Formula A28(A28)

Formula A29(A29)

Formula A30(A30)

Formula A31(A31)

Formula A32(A32)
Formula A32.

Delayed Rectifier Potassium Current IKdrI


Formula A33(A33)

Formula A34(A34)

Formula A35(A35)

Formula A36(A36)
Formula A36.

Leak Current ILI


Formula A37(A37)
Formula A37.

Synaptic Currents.

AMPA current IAMPAEI 

Formula A38(A38)
Formula A38.

NMDA current INMDAEI 

Formula A39(A39)
Formula A39.

GABAA variables 

Formula A40(A40)
Formula A40.

Architecture

The excitatory-to-excitatory synaptic footprint shape (Eqs. 10, A17, and A20) is used with the discrete function w(j), where

Formula A41(A41)
The basic unit length is the footprint length Formula A41, corresponding to a typical width of a cortical column, ~0.5 mm. Values of L/{sigma} vary from 16 to 32, and {rho} varies from 8 to 64. The excitatory-to-inhibitory synaptic footprint shape (Eqs. A38 and A39) is Formula A41. The inhibitory-to-excitatory synaptic footprint shape (Eq. A24) is wIE (j) =Formula A41] with Formula A41(Golomb and Ermentrout 2002Go).

Numerical methods

Simulations were performed using the fourth-order Runge-Kutta method with time step Formula A41. Varying {rho} above 8 has little effect on the results. Synaptic fields (Eqs. A17, A20, A24, A38, and A39) were computed using fast Fourier transform. Bifurcation diagrams (Figs. 8 and 10) were calculated with the software package XPPAUT (Ermentrout 2002Go).


 APPENDIX B
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Definition of local synchrony measure

The synchrony measure {chi} quantifies the normalized average voltage fluctuations, where the average is taken over a local population of excitatory neurons within a certain distance of order {sigma} from a specific neuron. Here, we choose this neuron to be in the middle of the chain and compute averages over all the neurons within a distance {sigma} from it. The population average voltage is

Formula 1(B1)
The variance of the time fluctuation of V(t) is

Formula 2(B2)
where Formula 2dt . . . denotes time average over a large time, Tm. After normalization of {sigma}V to the average over the population of the single-cell membrane potentials

Formula 3(B3)
we define a synchrony measure {chi}({rho}) for the activity of a system with cell density {rho} by

Formula 4(B4)
This synchrony measure, {chi}({rho}), is between 0 and 1. In the limit {sigma}{rho}->{infty} it behaves as

Formula 5(B5)
where a is a constant. In particular, Formula 5 if the system is fully synchronized [i.e., Formula 5 for all i], and Formula 5 if the state of the system is asynchronous. Asynchronous states are unambiguously characterized only when the number {sigma}{rho} of neurons is infinite (Golomb and Hansel 2000Go). To test the level of synchrony of a neuronal system in simulations, one needs to calculate {chi}({rho}) for various values of {rho} and determine the asymptotic value for very large {sigma}{rho} (Hansel and Sompolinsky 1996Go).


 GRANTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The research was supported by United States-Israel Binational Science Foundation Grant 2003019 to D. Golomb. G. B. Ermentrout is funded by National Science Foundation and National Institute on Mental Health.


 ACKNOWLEDGMENTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank Y. Grossman for helpful discussions and Y. Amitai and D. Hansel for careful reading of the manuscript.


 FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1 As a control parameter varies, such bursting states are obtained from tonic, periodic state via a period doubling or a torus (Hopf of periodics) bifurcation (Mandelblat et al. 2001Go). Back

2 The auxiliary variables xNMDA is used just for determining sAMPA, and is not included in the analysis. Back

Address for reprint requests and other correspondence: D. Golomb, Dept. of Physiology, Faculty of Health Sciences, P.O. Box 653, Ben-Gurion University, Be'er-Sheva, Israel 84105 (E-mail: golomb{at}bgu.ac.il)


 REFERENCES
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX A
 APPENDIX B
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Amit DJ and Brunel N. Model of global spontaneous activity and local structured activity during delay periods in the cerebral cortex. Cereb Cortex 7: 237–252, 1997.[Abstract/Free Full Text]

Beierlein M, Gibson JR, and Connors BW. Two dynamically distinct inhibitory networks in layer 4 of the neocortex. J Neurophysiol 90: 2987–3000, 2003.[Abstract/Free Full Text]

Bertram R, Butte MJ, Kiemel T, and Sherman A. Topological and phenomenological classification of bursting oscillations. Bull Math Biol 57: 413–439, 1995.[Web of Science][Medline]

Bibbig A, Faulkner HJ, Whittington MA, and Traub RD. Self-organized synaptic plasticity contributes to the shaping of gamma and beta oscillations in vitro. J Neurosci 21: 9053–9067, 2001.[Abstract/Free Full Text]

Cardoso de Oliveira S, Gribova A, Donchin O, Bergman H, and Vaadia E. Neural interactions between motor cortical hemispheres during bimanual and unimanual arm movements. Eur J Neurosci 14: 1881–1896, 2001.[CrossRef][Web of Science][Medline]

Castro-Alamancos MA and Rigas P. Cellular mechanisms of 7–14 Hz oscillations in the motor cortex. Soc Neurosci Abstr 30: 641.3, 2004.

Compte A, Brunel N, Goldman-Rakic PS, and Wang XJ. Synaptic mechanisms and network dynamics underlying spatial working memory in a cortical network model. Cereb Cortex 10: 910–923, 2000.[Abstract/Free Full Text]

Connors BW and Amitai Y. Synchronized excitation and inhibition driven by intrinsically bursting neurons in neocortex. In: The Cortical Neuron, edited by Gutnick MJ and Mody I. New York: Cambridge Univ. Press, 1995.

Connors BW, Gutnick MJ, and Prince DA. Electrophysiological properties of neocortical neurons in vitro. J Neurophysiol 48: 1302–1320, 1982.[Abstract/Free Full Text]

Destexhe A, Mainen ZF, and Sejnowski TJ. Synthesis of models for excitable membranes, synaptic transmission and neuromodulation using a common kinetic formalism. J Comput Neurosci 1: 195–230, 1994.[CrossRef][Medline]

Ermentrout B. The analysis of synaptically generated traveling waves. J Comput Neurosci 5: 191–208, 1998.[CrossRef][Web of Science][Medline]

Ermentrout B. Simulating, Analazing, and Animating Dynamical Systems: A Guide to XPPAUT for Researchers and Students (software, environment, tools). Philadelphia, PA: Society for Industrial and Applied Mathematics, 2002.

Ermentrout B. Dynamical consequences of fast-rising, slow-decaying synapses in neuronal networks. Neural Comput 15: 2483–2522, 2003.[CrossRef][Web of Science][Medline]

Fleidervish IA, Binshtok AM, and Gutnick MJ. Functionally distinct NMDA receptors mediate horizontal connectivity within layer 4 of mouse barrel cortex. Neuron 21: 1055–1065, 1998.[CrossRef][Web of Science][Medline]

Fleidervish IA, Friedman A, and Gutnick MJ. Slow inactivation of Na+ current and slow cumulative spike adaptation in mouse and guinea pig neocortical neurones in slices. J Physiol 493: 83–97, 1996.[Abstract/Free Full Text]

Fleidervish IA and Gutnick MJ. Kinetics of slow inactivation of persistent sodium current in layer V neurons of mouse neocortical slices. J Neurophysiol 76: 2125–2130, 1996.[Abstract/Free Full Text]

Flint AC and Connors BW. Two types of network oscillations in neocortex mediated by distinct glutamate receptor subtypes and neuronal populations. J Neurophysiol 75: 951–957, 1996.[Abstract/Free Full Text]

Frankenhaeuser B and Hodgkin AL. The action of calcium on the electrical properties of squid axons. J Physiol 137: 218–244, 1957.[Free Full Text]

Funahashi S, Bruce CJ, and Goldman-Rakic PS. Mnemonic coding of visual space in the monkey's dorsolateral prefrontal cortex. J Neurophysiol 61: 331–349, 1989.[Abstract/Free Full Text]

Gil Z and Amitai Y. Properties of convergent thalamocortical and intracortical synaptic potentials in single neurons of neocortex. J Neurosci 16: 6567–6578, 1996.[Abstract/Free Full Text]

Giugliano M, Darbon P, Arsiero M, Luscher HR, and Streit J. Single-neuron discharge properties and network activity in dissociated cultures of neocortex. J Neurophysiol 92: 977–996, 2004.[Abstract/Free Full Text]

Golomb D. Models of neuronal transient synchrony during propagation of activity through neocortical circuitry. J Neurophysiol 79: 1–12, 1998.[Abstract/Free Full Text]

Golomb D and Amitai Y. Propagating neuronal discharges in neocortical slices: computational and experimental study. J Neurophysiol 78: 1199–1211, 1997.[Abstract/Free Full Text]

Golomb D and Ermentrout GB. Continuous and lurching traveling pulses in neuronal networks with delay and spatially decaying connectivity. Proc Natl Acad Sci USA 96: 13480–13485, 1999.[Abstract/Free Full Text]

Golomb D and Ermentrout GB. Slow excitation supports propagation of slow pulses in networks of excitatory and inhibitory populations. Phys Rev E 65: 061911, 2002.[CrossRef]

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

Golomb D, Hansel D, and Mato G. Mechanisms of synchrony of neural activity in large networks. In: Handbook of Biological Physics, edited by Moss F and Gielen S: Elsevier Science, 2001, p. 887–968.

Golomb D and Rinzel J. Dynamics of globally coupled inhibitory neurons with heterogeneity. Phys Rev E 48: 4810–4814, 1993.[CrossRef]

Golomb D and Rinzel J. Clustering in globally coupled inhibitory neurons. Physica D 72: 259–282, 1994.[CrossRef]

Golomb D, Wang XJ, and Rinzel J. Propagation of spindle waves in a thalamic slice model. J Neurophysiol 75: 750–769, 1996.[Abstract/Free Full Text]

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

Hansel D and Mato G. Existence and stability of persistent states in large neuronal networks. Phys Rev Lett 86: 4175–4178, 2001.[CrossRef][Web of Science][Medline]

Hansel D and Sompolinsky H. Chaos and synchrony in a model of a hypercolumn in visual cortex. J Comput Neurosci 3: 7–34, 1996.[CrossRef][Web of Science][Medline]

Hoppensteadt FC and Izhikevich EM. Weakly Connected Neural Networks. New York: Springer-Verlag, 1997.

Izhikevich EM. Neural excitability, spiking and bursting. Int J Bifurcation Chaos 10: 1171–1266, 2000.[CrossRef][Web of Science]

Jahr CE and Stevens CF. Voltage dependence of NMDA-activated macroscopic conductances predicted by single-channel kinetics. J Neurosci 10: 3178–3182, 1990.[Abstract]

Johnston D and Wu SM-S. Foundations of Cellular Neurophysiology. Cambridge, MA: MIT Press, 1995.

Jung R, Kiemel T, and Cohen AH. Dynamic behavior of a neural network model of locomotor control in the lamprey. J Neurophysiol 75: 1074–1086, 1996.[Abstract/Free Full Text]

Kawaguchi Y. Distinct firing patterns of neuronal subtypes in cortical synchronized activities. J Neurosci 21: 7261–7272, 2001.[Abstract/Free Full Text]

Kumar SS and Huguenard JR. Pathway-specific differences in subunit composition of synaptic NMDA receptors on pyramidal neurons in neocortex. J Neurosci 23: 10074–10083, 2003.[Abstract/Free Full Text]

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

Mainen ZF and Sejnowski TJ. Influence of dendritic structure on firing pattern in model neocortical neurons. Nature 382: 363–366, 1996.[CrossRef][Medline]

Mandelblat Y, Etzion Y, Grossman Y, and Golomb D. Period doubling of calcium spike firing in a model of a Purkinje cell dendrite. J Comput Neurosci 11: 43–62, 2001.[CrossRef][Web of Science][Medline]

Markram H and Tsodyks M. Redistribution of synaptic efficacy between neocortical pyramidal neurons. Nature 382: 807–810, 1996.[CrossRef][Medline]

McCormick DA. Neuronal networks: flip-flops in the brain. Curr Biol 15: R294–296, 2005.[CrossRef][Web of Science][Medline]

McCormick DA, Connors BW, Lighthall JW, and Prince DA. Comparative electrophysiology of pyramidal and sparsely spiny stellate neurons of the neocortex. J Neurophysiol 54: 782–806, 1985.[Abstract/Free Full Text]

McCormick DA, Shu Y, Hasenstaub A, Sanchez-Vives M, Badoual M, and Bal T. Persistent cortical activity: mechanisms of generation and effects on neuronal excitability. Cereb Cortex 13: 1219–1231, 2003.[Abstract/Free Full Text]

Pinsky PF and Rinzel J. Intrinsic and network rhythmogenesis in a reduced Traub model for CA3 neurons. J Comput Neurosci 1: 39–60, 1994.[CrossRef][Medline]

Quinlan EM, Lebel D, Brosh I, and Barkai E. A molecular mechanism for stabilization of learning-induced synaptic modifications. Neuron 41: 185–192, 2004.[CrossRef][Web of Science][Medline]

Quinlan EM, Olstein DH, and Bear MF. Bidirectional, experience-dependent regulation of N-methyl-D-aspartate receptor subunit composition in the rat visual cortex during postnatal development. Proc Natl Acad Sci USA 96: 12876–12880, 1999.[Abstract/Free Full Text]

Rinzel J and Ermentrout GB.Analysis of neural excitability and oscillations. In: Methods in Neuronal Modeling: From Ions to Networks (2nd ed.), edited by Koch C and Segev I. Cambridge. MA: MIT Press, 1998, p. 251–291.

Sah P. Ca2+-activated K+ currents in neurones: types, physiological roles and modulation. Trends Neurosci 19: 150–154, 1996.[CrossRef][Web of Science][Medline]

Silva LR, Amitai Y, and Connors BW. Intrinsic oscillations of neocortex generated by layer 5 pyramidal neurons. Science 251: 432–435, 1991.[Abstract/Free Full Text]

Steriade M. Neocortical cell classes are flexible entities. Nat Rev Neurosci 5: 121–134, 2004.[CrossRef][Web of Science][Medline]

Steriade M, Nunez A, and Amzica F. A novel slow (< 1 Hz) oscillation of neocortical neurons in vivo: depolarizing and hyperpolarizing components. J Neurosci 13: 3252–3265, 1993.[Abstract]

Stern P, Edwards FA, and Sakmann B. Fast and slow components of unitary EPSCs on stellate cells elicited by focal stimulation in slices of rat visual cortex. J Physiol 449: 247–278, 1992.[Abstract/Free Full Text]

Storm JF. Potassium currents in hippocampal pyramidal cells. Prog Brain Res 83: 161–187, 1990.[Web of Science][Medline]

Strogatz SH. Nonlinear Dynamics and Chaos, with Applications to Physics, Biology, Chemistry and Engineering. Reading: Addison-Wesley, 1994.

Su H, Alroy G, Kirson ED, and Yaari Y. Extracellular calcium modulates persistent sodium current-dependent burst-firing in hippocampal pyramidal neurons. J Neurosci 21: 4173–4182, 2001.[Abstract/Free Full Text]

Sutor B and Hablitz JJ. EPSPs in rat neocortical neurons in vitro. II. Involvement of N-methyl-D-aspartate receptors in the generation of EPSPs. J Neurophysiol 61: 621–634, 1989.[Abstract/Free Full Text]

Tabak J, Senn W, O'Donovan MJ, and Rinzel J. Modeling of spontaneous activity in developing spinal cord using activity-dependent depression in an excitatory network. J Neurosci 20: 3041–3056, 2000.[Abstract/Free Full Text]

Telfeian AE and Connors BW. Epileptiform propagation patterns mediated by NMDA and non-NMDA receptors in rat neocortex. Epilepsia 40: 1499–1506, 1999.[CrossRef][Web of Science][Medline]

Terman D. The transition from bursting to continuous spiking in excitable membrane model. J Nonlinear Sci 2: 135–182, 1992.

Traub RD, Jefferys JG, and Whittington MA. Enhanced NMDA conductance can account for epileptiform activity induced by low Mg2+ in the rat hippocampal slice. J Physiol 478: 379–393, 1994.[Abstract/Free Full Text]

Traub RD, Wong RK, Miles R, and Michelson H. A model of a CA3 hippocampal pyramidal neuron incorporating voltage-clamp data on intrinsic conductances. J Neurophysiol 66: 635–650, 1991.[Abstract/Free Full Text]

Tsau Y, Guan L, and Wu JY. Initiation of spontaneous epileptiform activity in the neocortical slice. J Neurophysiol 80: 978–982, 1998.[Abstract/Free Full Text]

Tsodyks MV and Markram H. The neural code between neocortical pyramidal neurons depends on neurotransmitter release probability. Proc Natl Acad Sci USA 94: 719–723, 1997.[Abstract/Free Full Text]

Vaadia E, Haalman I, Abeles M, Bergman H, Prut Y, Slovin H, and Aertsen A. Dynamics of neuronal interactions in monkey cortex in relation to behavioural events. Nature 373: 515–518, 1995.[CrossRef][Medline]

van Vreeswijk C, and Hansel D. Patterns of synchrony in neural networks with spike adaptation. Neural Comput 13: 959–992, 2001.[CrossRef][Web of Science][Medline]

Wang XJ. Synaptic basis of cortical persistent activity: the importance of NMDA receptors to working memory. J Neurosci 19: 9587–9603, 1999.[Abstract/Free Full Text]

Wang XJ. Synaptic reverberation underlying mnemonic persistent activity. Trends Neurosci 24: 455–463, 2001.[CrossRef][Web of Science][Medline]

Wang XJ and Buzsáki G. Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network model. J Neurosci 16: 6402–6413, 1996.[Abstract/Free Full Text]

Wang XJ and Rinzel J. Spindle rhythmicity in the reticularis thalami nucleus: synchronization among mutually inhibitory neurons. Neuroscience 53: 899–904, 1993.[CrossRef][Web of Science][Medline]

Wu JY, Guan L, and Tsau Y. Propagating activation during oscillations and evoked responses in neocortical slices. J Neurosci 19: 5005–5015, 1999.[Abstract/Free Full Text]

Yue C and Yaari Y. KCNQ/M channels control spike afterdepolarization and burst generation in hippocampal neurons. J Neurosci 24: 4614–4624, 2004.[Abstract/Free Full Text]




This article has been cited by other articles:


Home page
Cereb CortexHome page
D. Durstewitz and T. Gabriel
Dynamical Basis of Irregular Spiking in NMDA-Driven Prefrontal Cortex Neurons
Cereb Cortex, April 1, 2007; 17(4): 894 - 908.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
M. A. Castro-Alamancos, P. Rigas, and Y. Tawara-Hirata
Resonance (~10 Hz) of excitatory networks in motor cortex: effects of voltage-dependent ion channel blockers
J. Physiol., January 1, 2007; 578(1): 173 - 191.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
95/2/1049    most recent
00932.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (10)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Golomb, D.
Right arrow Articles by Ermentrout, G. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Golomb, D.
Right arrow Articles by Ermentrout, G. B.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2006 by the The American Physiological Society.