JN Fuel your research with LabChart
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 ISI 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 ISI Web of Science (6)
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 2003