Propagating Neuronal Discharges in Neocortical Slices: Computational and Experimental Study

David Golomb and Yael Amitai

Zlotowski Center for Neuroscience and Department of Physiology, Faculty of Health Sciences, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel

    ABSTRACT
Abstract
Introduction
Methods
Results
Discussion
References

Golomb, David and Yael Amitai. Propagating neuronal discharges in neocortical slices: computational and experimental study. J. Neurophysiol. 78: 1199-1211, 1997. We studied the propagation of paroxysmal discharges in disinhibited neocortical slices by developing and analyzing a model of excitatory regular-spiking neocortical cells with spatially decaying synaptic efficacies and by field potential recording in rat slices. Evoked discharges may propagate both in the model and in the experiment. The model discharge propagates as a traveling pulse with constant velocity and shape. The discharge shape is determined by an interplay between the synaptic driving force and the neuron's intrinsic currents, in particular the slow potassium current. In the model, N-methyl-D-aspartate (NMDA) conductance contributes much less to the discharge velocity than amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) conductance. Blocking NMDA receptors experimentally with 2-amino-5-phosphonovaleric acid (APV) has no significant effect on the discharge velocity. In both model and experiments, propagation occurs for AMPA synaptic coupling gAMPA above a certain threshold, at which the velocity is finite (non-zero). The discharge velocity grows linearly with the gAMPA for gAMPA much above the threshold. In the experiments, blocking AMPA receptors gradually by increasing concentrations of 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX) in the perfusing solution results in a gradual reduction of the discharge velocity until propagation stops altogether, thus confirming the model prediction. When discharges are terminated in the model by the slow potassium current, a network with the same parameter set may display discharges with several forms, which have different velocities and numbers of spikes; initial conditions select the exhibited pattern. When the discharge is also terminated by strong synaptic depression, there is only one discharge form for a particular parameter set; the velocity grows continuously with increased synaptic conductances. No indication for more than one discharge velocity was observed experimentally. If the AMPA decay rate increases while the maximal excitatory postsynaptic conductance (EPSC) a cell receives is kept fixed, the velocity increases by ~20% until it reaches a saturated value. Therefore the discharge velocity is determined mainly by the cells' integration time of input EPSCs. We conclude, on the basis of both the experiments and the model, that the total amount of excitatory conductance a typical cell receives in a control slice exhibiting paroxysmal discharges is only ~5 times larger than the excitatory conductance needed for raising the potential of a resting cell above its action potential threshold.

    INTRODUCTION
Abstract
Introduction
Methods
Results
Discussion
References

When synaptic gamma -aminobutyric acid-A (GABAA) inhibition is pharmacologically suppressed in neocortical slices, neuronal population discharges1 appear as responses to electrical stimulation above a certain strength. In extracellular recording, they are found to be abrupt, all-or-none field potentials (FPs). In intracellular recordings, they correspond to depolarizing shifts (DSs) in membrane potentials, above which rides a high-frequency train of action potentials. These discharges are referred to as "synchronous" or "epileptiform" (Gutnick et al. 1982), implying that adjacent cells tend to fire together.

Understanding the dynamics of discharge propagation is interesting for three reasons. First, by relating the discharge characteristics, such as the shape of the intracellularly measured potentials and the velocity of discharge propagation, to the intrinsic and synaptic neuronal properties, information about the underlying cortical structure can be revealed. For example, various cortical areas differ in the discharge velocities and shapes along the slice (Chervin et al. 1988; Neumann-Haefelin et al. 1996; Wadman and Gutnick 1993), most likely reflecting physiological and anatomic differences. Second, these paroxysmal discharges constitute in vitro models for interictal discharges in epilepsy: both events are relatively large, essentially all-or-none, reproducible, and relatively brief; have variable latencies; and seem to encompass almost all the neurons in a local cortical region in a synchronous manner (Chervin et al. 1988). Thus studying the mechanisms for the initiation and propagation of neuronal discharges could have therapeutic implications. Third, spatiotemporal propagating activities in the cortex have been observed recently in vivo (Nicolelis et al. 1996; Prechtl et al. 1996). Discharge propagation in slices is one of the simplest examples of such activity, and understanding its dynamics is a first step in the investigation of propagation effects in neuronal systems.

The paroxysmal discharges spread from the stimulus site and propagate along the slice. The experimental characteristics of the discharge propagation have been described in numerous experimental studies (e.g., Baldino et al. 1986; Chagnac-Amitai and Connors 1989a,b; Chervin et al. 1988; Connors 1984; Flint and Connors 1996; Gutnick et al. 1982; Wadman and Gutnick 1993). However, several questions regarding the dynamic mechanisms that lead to the creation, propagation, and cessation of discharges remain to be answered.

1) What are the spatiotemporal properties of the discharge propagation, and how are they related to the network architecture?

2) How do the velocity and shape of the discharge depend on the synaptic parameters, such as the strength of amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and N-methyl-D-aspartate (NMDA) synaptic efficacies, and the level of synaptic depression?

3) What is the relationship between the propagation velocity and the synaptic kinetics?

To answer these questions, we used both modeling and experimental techniques. We constructed a computational biophysical model that was consistent with the known cellular and synaptic properties of neocortical cells. The network model has a one-dimensional architecture with spatially decaying synaptic efficacies, and therefore spatiotemporal phenomena such as propagation can be investigated. The discharge shape and velocity were studied as functions of the synaptic coupling strengths, both for AMPA and NMDA, and the effect of synaptic depression on the network dynamics was investigated. Predictions raised by the model were tested experimentally in disinhibited rat cortical slices with the use of FP recordings and pharmacological manipulations. Combinations of modeling and experimental approaches have been used before in studying discharge propagation in hippocampal slices (Miles et al. 1988; Traub et al. 1993). Although in these hippocampal models the single cell is an endogenous burster (Golomb and Rinzel 1996; Pinsky and Rinzel 1994; Traub et al. 1991), the single cells in our model are regular spikers (Gutnick and Crill 1995), and the bursting activity during the discharge is a network phenomenon. A preliminary report of the present work has been published in abstract form (Golomb and Amitai 1996).

    METHODS
Abstract
Introduction
Methods
Results
Discussion
References

Single-cell models

When inhibition is completely blocked, both regular-spiking (RS) and intrinsic bursting excitatory cells participate in the collective activity and exhibit a DS during the discharge propagation (Connors and Amitai 1993; Gutnick et al. 1982). Because there is currently no reason to assume that the intrinsic bursting property in neocortex plays a major role in the activity, we treat here, for simplicity, only one type of excitatory pyramidal cortical cell, corresponding to the RS type. The single-cell dynamics is described by a single-compartment Hodgkin-Huxley type model with the use of a set of coupled differential equations (Hansel and Sompolinsky 1996)
<IT>C</IT><FR><NU>d<IT>V</IT>(<IT>x</IT>, <IT>t</IT>)</NU><DE>d<IT>t</IT></DE></FR> = −<IT>I</IT><SUB>Na</SUB><IT> − I</IT><SUB>NaP</SUB><IT> − I</IT><SUB>Kdr</SUB><IT> − I</IT><SUB>KA</SUB><IT> − I</IT><SUB>K-slow</SUB><IT></IT>
− <IT>I</IT><SUB>L</SUB><IT>− I</IT><SUB>AMPA</SUB><IT>− I</IT><SUB>NMDA</SUB><IT>+ I</IT><SUB>app</SUB> (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 side incorporates an applied current Iapp, and the following intrinsic and synaptic currents (Gutnick and Crill 1995): the fast sodium current INa, the persistent sodium current INaP, the delayed-rectifier potassium current IKdr, the A-type potassium current IKA, the slow potassium current IK-slow, the leak current IL, the AMPA current IAMPA, and the NMDA current INMDA.

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 on the order of several tens to hundreds of milliseconds (Storm 1990). These currents are either calcium dependent, such as IC and IAHP (Pinsky and Rinzel 1994; Sah 1996; Traub et al. 1991), or voltage dependent, such as IM (McCormick et al. 1993; Yamada et al. 1989). They are responsible for the adaptation in RS cortical cells. The equations and parameters of the model are given in the APPENDIX.

Synaptic models

Neocortical cells receive fast AMPA-mediated excitatory postsynaptic potentials (EPSPs) and slow NMDA-mediated EPSPs from neighboring excitatory cells (Douglas and Martin 1990; Gil and Amitai 1996a). Synaptic transmission occurs when the presynaptic cell emits a spike, i.e., when its potential rises above a certain voltage level. A gating variable s for an AMPA or an NMDA receptor, representing the fraction of open channels, is modeled according to
<FR><NU>d<IT>s</IT></NU><DE>d<IT>t</IT></DE></FR>= <IT>k<SUB>f</SUB>T</IT><SUB>Glu</SUB><IT>s</IT><SUB>∞</SUB>(<IT>V</IT><SUB>pre</SUB>)(1 − <IT>s</IT>) − <IT>k</IT><SUB><IT>r</IT></SUB><IT>s</IT> (2)
where Vpre is the presynaptic potential, TGlu is the normalized number of presynaptic glutamate vesicles, sinfinity (V) = {1 + exp[-(V - theta s)/sigma s]}-1, theta s = -20 mV is the presumed voltage threshold for synaptic release, and sigma s = 2 mV (Destexhe et al. 1994; Wang and Rinzel 1993). The presynaptic voltage is more depolarized than theta s only when the presynaptic cell emits a spike. The decay time of the fast AMPA synapses is kr = 0.2 ms-1, and that of the slow NMDA synapse is krN = 0.0067 ms-1 (Jonas and Spruston 1994; Stern et al. 1992); kf = 1 ms-1.

Synaptic depression is modeled phenomenologically as a depletion process of glutamate presynaptic vesicles (Stevens and Wang 1995)
<FR><NU>d<IT>T</IT><SUB>Glu</SUB></NU><DE>d<IT>t</IT></DE></FR>= −<IT>k<SUB>t</SUB>s</IT><SUB>∞</SUB>(<IT>V</IT><SUB>pre</SUB>)<IT>T</IT><SUB>Glu</SUB><IT>+ k</IT><SUB><IT>&cjs1726;</IT></SUB>(1 − <IT>T</IT><SUB>Glu</SUB>) (3)
When an action potential arrives, sinfinity (Vpre) is close to 1; the normalized number of presynaptic vesicles TGlu, 0 <=  TGlu <=  1, decreases with a rate kt. When there is no presynaptic activity, sinfinity (Vpre) is close to 0, and TGlu recovers with a rate k&cjs1726;, much slower than kt. If kt = 0 there is no synaptic depression.

The total synaptic conductance a neuron receives, gSin (for AMPA or NMDA) is calculated by summing the synaptic variable of each of its presynaptic neurons as described below. The AMPA current is
<IT>I</IT><SUB>AMPA</SUB><IT>= g</IT><SUB>AMPA</SUB><IT>S</IT><SUP>in</SUP><SUB>AMPA</SUB>(<IT>V − V</IT><SUB>Glu</SUB>) (4)
The NMDA current depends also on the postsynaptic voltage, and is calculated approximately by multiplying the presynaptic term by a sigmoid function fNMDA(V) of the postsynaptic voltage (Destexhe et al. 1994; Traub et al. 1991). The NMDA current is therefore
<IT>I</IT><SUB>NMDA</SUB><IT>= g</IT><SUB>NMDA</SUB><IT>S</IT><SUP>in</SUP><SUB>NMDA</SUB><IT>f</IT><SUB>NMDA</SUB>(<IT>V</IT>)(<IT>V − V</IT><SUB>Glu</SUB>) (5)
Axonal propagation delays are neglected because the propagation velocity of action potentials in axons (on the order of 1 m/s; see DISCUSSION) (Gil and Amitai 1996b; Haberly 1990) is considerably faster than the velocity of discharge propagation (~15 cm/s).

Network architecture

Rat coronal cortical slices used in experiments (e.g., Connors 1984; Gutnick et al. 1982; present experiments) are narrow in one dimension (thickness 400 µm) and extend ~1.5 × 0.2 cm in the two other dimensions. Our goal in this modeling work is to study the discharge propagation along the slice as neurons in consecutive columns are recruited, and not the recruitment of cells across the layers in one column. Therefore the layer structure is not modeled here. Instead, following previous modeling work on hippocampal slices (Miles et al. 1988; Traub et al. 1993) and simple neural networks (Ermentrout and McLeod 1993; Idiart and Abbott 1993; Wilson and Cowan 1973), we present a model of a one-dimensional system. Neurons are equally distributed along the interval 0 <= x <=  L, where L is the slice length and x is the neuron position. The number of neurons is N and the position of the ith neuron, <= i <=  N, is xi = iL/N. The interaction between neurons is assumed to decay with the distance between them (Fig. 1A). The term "synaptic footprint shape," w(x), denotes the functional dependence of the synaptic connectivity on the distance between the pre- and the postsynaptic cells. It is assumed here to be symmetrical. The coupling footprint length lambda  is the typical decay length of the synaptic footprint shape w(x). An exponential footprint shape w(x) = Ae-|x|/lambda is chosen (Fig. 1B) (Golomb et al. 1996; Traub et al. 1993), where the normalization constant A is determined below. The model is studied in the regime
λ ≪ <IT>L</IT> (6)
We assume that a neuron receives many synaptic inputs from neurons within a footprint length
<IT>N</IT>λ/<IT>L</IT>≫ 1 (7)
and therefore we use the continuum approximation. The total synaptic conductance gSin affecting a neuron is
<IT>gS</IT><SUP>in</SUP>(<IT>x</IT>, <IT>t</IT>) = <IT>g</IT><LIM><OP>∫</OP></LIM><SUP><IT>L</IT></SUP><SUB>0</SUB>d<IT>x</IT>′ρ(<IT>x</IT>′)<IT>w</IT>(<IT>x − x</IT>′)<IT>s</IT>(<IT>x</IT>′, <IT>t</IT>) (8)
where rho (x) = N/L is the cell density. The quantities g, s, and Sin have a subscript AMPA or NMDA for the two types of receptors. Neurons near the edges receive synaptic inputs from their existing neighbors according to the same rule, but do not get input from nonexisting "neurons" outside of the slice model (open boundary conditions). To make the model as independent as possible of the number of cells N and the footprint length lambda  (Golomb et al. 1996; Hansel and Sompolinsky 1996), we normalize the function w(x) such that the total synaptic conductance received by a cell distant from the edges when all the synaptic channels are open is g. Equations 8 and 6 imply that the normalization condition is
<FR><NU><IT>N</IT></NU><DE>L</DE></FR><LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>−∞</SUB>d<IT>x w</IT>(<IT>x</IT>) = 1 (9)
The integro-differential Equations 1, 4, 5, and 8, together with the differential equations for the auxiliary variables, are discretized on a grid, as described in the APPENDIX, and solved numerically.


View larger version (18K):
[in this window]
[in a new window]
 
FIG. 1. Schematic diagrams of model architecture and cortical slice. A: model has 1-dimensional architecture. Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and N-methyl-D-aspartate (NMDA) coupling strengths (maximal synaptic conductances) decay with distance between pre- and postsynaptic cells. Typical decay length lambda  is called synaptic footprint length. B: for exponential footprint shape, lambda  is spatial distance for which coupling strength reaches 1/e of its maximal value. C: schematic diagram of cortical slice with 2 recording electrodes and 1 stimulating electrode.

Initial conditions and parameter survey

Our goal is to understand the dynamics of activity discharges evoked experimentally by a brief local stimulation when the slice is at rest. Therefore we initiate our model from a state at which all the neurons, except a group at the left edge, are at their resting state. A wave is initiated by depolarizing a group of neurons at the left edge (0 <=  x <=  0.06) to 0 mV. The auxiliary variables of the neurons are chosen to be at their steady-state values with the appropriate V. The depolarized neurons may recruit resting neurons to the activity and initiate a propagating discharge.

Not all the parameters of our model can be determined experimentally. As was done in previous works (Golomb et al. 1994, 1996), we choose a reference parameter set that is consistent with many experimental results. Then we systematically change some parameters to investigate their effects on the network dynamics. There are parameter regimes in which several dynamic states occur for the same parameter set; the discharge form is selected by initial conditions-for example, the number of cells initially excited. Often, our goal is to modify one parameter value and to find the parameter interval in which one dynamic state, such as a propagating discharge with a specific number of spikes (see RESULTS), occurs. In this case, we simulate the system with a first parameter value and find the state. The values of the variables at the right quarter of the network, determined at the time the discharge has just arrived at the point L - lambda , are stored. Then, we simulate the network with the new parameter value. The initial conditions of the variables at the left quarter of the network are chosen to be the values of the variables at the right quarter that were stored previously. All the other neurons in the model start from rest. The borders of a parameter range for which a dynamic state is observed are determined with the use of the bisection algorithm (Press et al. 1992).

Experimental methods

The methods for preparing and maintaining slices of neocortex have been described (Amitai 1994). Briefly, Wistar rats (4-6 wk old, 100-150 g) were anesthetized with pentobarbital sodium and decapitated, and the brains quickly removed into a chilled (~4°C), oxygenated artificial cerebrospinal fluid (ACSF). The ACSF contained (in mM) 124.0 NaCl, 5.0 KCl, 2.0 MgSO4, 1.25 NaHPO4, 2.0 CaCl2, 26.0 NaHCO3, and 10.0 dextrose, saturated with 95% O2-5% CO2, pH 7.4. Coronal slices 400 µm thick were cut on a vibratome (Campden Instruments) from the primary somatosensory cortex and were kept in holding bottles that contained ACSF at room temperature, continuously bubbled with 95% O2-5% CO2. Recordings were performed in a fluid-gas interface chamber thermostatically controlled to 35-36°C after >= 1 h of incubation. All drugs were bath applied.

Two FP electrodes (filled with 1 M NaCl, 6-8 MOmega ) were placed ~7-8 mm apart in layers II/III. The cortex was stimulated by 0.1-ms, 0.01- to 0.05-mA pulses at 0.1 Hz delivered by a bipolar microelectrode made from sharpened tungsten wires and placed on the border of the white matter >= 0.5 mm away from the closer recording electrode (Fig. 1C). Recordings were made with an Axoprobe amplifier (Axon Instruments). Data were recorded on-line at 10 kHz and analyzed with a software program written under visual C++ (Labview, National Instruments).

    RESULTS
Abstract
Introduction
Methods
Results
Discussion
References

Single-cell dynamics

The model of the single excitatory cell includes two time scales: one (1-10 ms) of the action-potential-generating currents (INa, IKdr, and even IKA) and a second, on the order of 100 ms, of the activation of IK-slow. Adaptation is produced as an interplay between these two time scales. The neuronal behavior with IK-slow blocked is shown in Fig. 2, A and B. When a depolarized current above a critical value is injected into the cell, the cell fires tonically (Fig. 2AI). If the applied current is even higher, the neuron membrane potential may approach a depolarized plateau following a period of damped oscillatory firing (Fig. 2AII). Figure 2B shows that as the applied current Iapp increases, the neuron dynamics shifts from rest to tonic firing and then to a depolarized plateau (with a regime where both of the later 2 states are observed).


View larger version (27K):
[in this window]
[in a new window]
 
FIG. 2. Dynamics of single pyramidal cell model. A and B: slow potassium current IK-slow is blocked (no adaptation). A: when neuron is at rest, injecting applied current of 2.5 µA/cm2 shifts it into spiking mode (AI). Injecting Iapp at 7 µA/cm2 shifts neuron into depolarized plateau after period of fast damped firing (AII). B: properties of rest and oscillatory (tonic firing) states of single-neuron model with applied current Iapp. Narrow solid line: stable rest state. Narrow dashed line: unstable rest state (which cannot be reached experimentally). For oscillatory state, maximum and minimum voltages over cycle are denoted by wide solid lines (stable oscillations) and wide dashed lines (unstable oscillations). There is a bistable regime (1.58 µA/cm2 < Iapp < 6.56 µA/cm2) for which both tonic firing and depolarized plateau can be sustained. Arrows: the 2 Iapp values shown in A. C and D: IK-slow is intact. C: when neuron is at rest, injecting applied current of 2.5 µA/cm2 causes firing with interspike interval that increases with time until it reaches constant value. Injecting stronger current, such as 7 µA/cm2 (CII), decreases interspike interval both at onset of spikes and at steady-state firing. Minimum voltage during each firing cycle decreases with adaptation. D: properties of rest and oscillatory (tonic firing) states of single-neuron model at long times, after complete adaptation. Symbols as in B. Except for very narrow regime near onset of oscillations (Iapp gsim  0.33 µA/cm2), minimum voltage during firing is higher than at rest potential. Arrows: the 2 Iapp values shown in C.

The neuron dynamics is different when IK-slow is intact. When the neuron is initially at rest, a strong enough current injection causes the model neuron to fire, with the interspike interval increasing with time until it reaches a constant value (Fig. 2C). This behavior occurs because of the time of a few cycles that is needed for IK-slow to build up, hyperpolarize the neuron, and reduce its firing rate; it mimics the adaptation found experimentally in RS neurons (Gutnick and Crill 1995). Both the first interspike interval and the long-time interspike interval decrease with the injected current amplitude. The long-time dynamic behavior of the neuron is presented in Fig. 2D as a function of Iapp. A comparison of Fig. 2D with Fig. 2B reveals that including IK-slow increases the Iapp range for which the neuron fires tonically instead of going to a high plateau. Except for a tiny Iapp regime near the onset of oscillations (Iapp = 0.34 µA/cm2), the minimum voltage during the oscillation cycle is higher than the rest state with no applied current. This phenomenon, called here overshooting, is consistent with the observed behavior of RS pyramidal neurons (reviewed by Amitai and Connors 1994). This phenomenon, together with the existence of the high plateau, is relevant for the behavior of a cell within a network, as is shown below.

Discharge propagation without synaptic depression

As a first step, we investigate a case in which neuronal activity is terminated by IK-slow only, without synaptic depression. Stimulating the left edge of the slice creates activity that propagates by continuously recruiting new neighboring cells. The firing rastergram of the neurons in the network is shown in Fig. 3A, and membrane potential traces of five neurons along the network are plotted in Fig. 3B. Except for a region near the edges, the voltage trace (and the whole dynamics) of a neuron in a position x1 is identical to the voltage trace of a neuron in a position x2, but is shifted in time with a delay (x1 - x2)/&cjs1726;D. The discharge velocity &cjs1726;D is constant along the slice. Therefore the discharge is dynamically a traveling pulse2 (Cross and Hohenberg 1993)
<IT>V</IT>(<IT>x</IT>, <IT>t</IT>) = <IT><A><AC>V</AC><AC>˜</AC></A></IT>(<IT>x − </IT>&cjs1726;<SUB>D</SUB><IT>t</IT>) (10)
During the activity, each neuron fires several action potentials (7 in the present case). Time evolutions of synaptic and intrinsic auxiliary variables of one neuron within a network are plotted in Fig. 3C. The synaptic inputs SinAMPA and SinNMDA begin to rise from zero under the effect of neuronal activity to the left of the neuron. As a result, the neuron starts firing, sends synaptic output sAMPA and sNMDA to its neighbors, and helps to recruit more neurons on its right. During that time, the slow potassium current IK-slow builds up. When it is strong enough, it terminates the neuron's activity. As long as the neuron fires, its AMPA synaptic output variable sAMPA rises during each spike to about the same value and drops down between spikes. After the activity terminates, sAMPA returns rapidly to zero. The synaptic variable of the slow NMDA synapse, sNMDA, stays at a high value during the activity and decays slowly on its termination.


View larger version (12K):
[in this window]
[in a new window]
 
FIG. 3. Propagation of discharges in cortical slice model. There is no synaptic depression (kt = 0), gAMPA = 0.31 mS/cm2, gNMDA = 0.25 mS/cm2. A: rastergram of excitatory cells. Only firing times from every 4th cell are shown. Arrows at right: position of 5 cells along slice whose voltage time courses are plotted in B. Initially cells at left (x <=  0.06) are depolarized (V = 10 mV) and all others are at rest. Discharge propagates to right and far from edges maintains its shape; every neuron fires 7 spikes during discharge. C: time courses of internal variables of neuron denoted by a. From top: total AMPA and NMDA synaptic conductance, SinAMPA and SinNMDA, that this neuron receives (the neuron's "input"); AMPA and NMDA auxiliary variables, sAMPA and sNMDA, of postsynaptic synapses connecting this neuron to others (the neuron's "output"); slow potassium activation variable z. Auxiliary variables are normalized such that 1 means that all corresponding channels are open.

The input synaptic field varies more slowly than the interspike interval. To analyze the discharge shape, we use a crude approximation and separate the slow time scales of the synaptic input and of IK-slow from the faster time scale of the spike-generating currents. Except during action potentials, the synaptic driving force V - VGlu (Eq. 4) is approximately constant. Therefore the effect of the excitatory synaptic input is similar to the effect of an applied current. As shown above, the cortical cell model exhibits overshoot in response to current injection. Similarly, the minimum voltage during the active period is above the rest potential. This behavior is consistent with the experimental observation that the spikes ride on a depolarizing envelope (Gutnick et al. 1982).

The shape of a voltage trace of a neuron during the discharge is determined by an interplay between the input synaptic field and the slow potassium current IK-slow. Voltage traces for two values of AMPA conductance gAMPA and two values of gK-slow are shown in Fig. 4. For moderate synaptic input (e.g., 0.45 mS/cm2), the neuron fires a train of spikes with an amplitude that decays only slowly with time. At higher gAMPA (0.8 mS/cm2 in Fig. 4), the neuron's potential tends to converge to the high plateau. The amplitude of consecutive action potentials decreases rapidly; the peak voltage decreases and the minimum voltage (between the action potentials) increases. The slow potassium current terminates the activity in both cases. The stronger it is, the faster is the discharge termination.


View larger version (9K):
[in this window]
[in a new window]
 
FIG. 4. Voltage time traces of neuron during discharge, gNMDA = 0. AMPA conductance gAMPA: 0.45 mS/cm2 (left traces), 0.8 mS/cm2 (right traces). Slow potassium conductance gK-slow: 1.0 mS/cm2 (top traces), 2.0 mS/cm2 (bottom traces).

Dependence on AMPA synaptic strength

MODEL. We assess quantitatively how the number of spikes fired by each neuron, and the discharge velocity &cjs1726;D, vary with synaptic strength. Figure 5A presents the number of spikes and &cjs1726;D of a discharge as a function of gAMPA when NMDA receptors are blocked. Discharges can propagate only if gAMPA is above a threshold level gAMPA,c. At that threshold conductance, &cjs1726;D is finite (non-zero). For a wide gAMPA range, there are multiple discharge forms. The initial conditions, i.e., the number of initially depolarized cells, select the discharge pattern. For the same gAMPA, discharges with a larger number of spikes propagate faster. The reason for this is the fact that a neuron just in front of a discharge involving a larger number of spikes is affected by a stronger synaptic field. The increase in the velocity is more prominent at smaller spike numbers. For example, at gAMPA = 0.31 mS/cm2 (*), the velocity increase from the three-spike discharge to the four-spike discharge is 9.4%, whereas the velocity increase from the four-spike discharge to the five-spike discharge is only 3.6%. This effect is due to the exponential decay of the synaptic footprint shape, because each additional spike is spatially more distant, and contributes less to the synaptic field. The effects of adding more than seven spikes cannot be detected, and the velocity curve looks almost continuous.


View larger version (15K):
[in this window]
[in a new window]
 
FIG. 5. Number of spikes in propagating discharge (top) and discharge velocity &cjs1726;D (bottom) vs. AMPA synaptic conductance strength gAMPA (A) and NMDA synaptic conductance strength gNMDA (B). There is no synaptic depression. In A, NMDA conductance is blocked (gNMDA = 0). Asterisk: value gAMPA = 0.31 mS/cm2 (see text). In B, gAMPA = 0.31 mS/cm2. Arrow: gNMDA value of Fig. 3. Multiple discharge shapes with different numbers of spikes and velocities can occur for same parameter value.

For gAMPA much larger than gAMPA,c, the relationship between the velocity &cjs1726;D and gAMPA is close to linear. At gAMPA values larger than those shown in Fig. 5A, the membrane voltage during the discharge tends to go to a high plateau before returning to rest. In such cases, it is difficult to determine what a "spike" is, and therefore to count spikes, but &cjs1726;D continues to grow approximately linearly with gAMPA.

Simple scaling arguments based on the normalization condition (Eq. 9) show that the discharge velocity depends linearly on the synaptic footprint length lambda  (Golomb et al. 1996).

EXPERIMENT. The dependence of &cjs1726;D on the synaptic strength has never been tested experimentally. We tested the relationship between &cjs1726;D and gAMPA experimentally. We used cortical slices in which GABAA inhibition was blocked by bicuculline methiodide (10 µM), and excitatory transmission was manipulated pharmacologically. Initially, NMDA receptors were completely blocked by application of 2-amino-5phosphonovaleric acid (APV, 30-50 µM). To examine the effect of changes in gAMPA, AMPA receptors were blocked gradually by adding 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), at slowly increasing concentrations, to the perfusing solution, starting from 0.1-0.3 µM. The control (0 CNQX) &cjs1726;D varied between 13 and 19 cm/s among slices (n = 7). The first decrease in &cjs1726;D was measured at concentrations of 0.2-0.4 µM CNQX, varying between slices. Propagation continued to slow gradually with increasing CNQX concentrations until it was abolished altogether at concentrations of 0.6-0.8 µM. (Fig. 6, n = 4). The lowest velocity measured was between 0.25 and 0.5 of the control velocity. In higher CNQX concentrations, there was a change in the FP shape at both proximal and distal recording micropipettes.


View larger version (24K):
[in this window]
[in a new window]
 


View larger version (14K):
[in this window]
[in a new window]
 
FIG. 6. A: local field potential (FP) recording from 2 electrodes for several 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX) concentrations. Narrow line: FP measured by electrode closer to stimulus. Wide line: FP measured by more distant electrode. As CNQX concentration increases, latency between discharge arrival to the 2 electrodes increases. B: discharge velocity at several CNQX concentrations measured in 4 experiments. Results from each experiment are denoted by different line types (------, · · ·, - - -, and - - -) and different symbols (bullet , black-square, open circle  and black-diamond , respectively). Control values (at zero CNQX level) are not connected by lines. Error bars: SD, calculated by averaging over 10-20 discharges at each concentration. FP traces in A correspond to - - -.

The dependency of &cjs1726;D on the initial conditions was examined by changing the stimulus strength Istim at several CNQX concentrations. Propagation was observed only if Istim was above a threshold Istim,c, which increased with increasing CNQX concentrations. Close to Istim,c, small changes in &cjs1726;D resulted in pronounced changes in the latency to the proximal FP (Fig. 7) (see also Gutnick et al. 1982). Nevertheless, there was no significant change in &cjs1726;D as a function of Istim at all CNQX concentrations. Thus our experimental results confirmed two main predications of the model: 1) &cjs1726;D was found to depend on the strength of AMPA transmission in the slice and 2) a threshold level of AMPA strength was established (in particular, the propagation stops at a finite velocity). The experimental results could not confirm the multiple velocities for a given condition set predicted by the model at gAMPA near the threshold.


View larger version (17K):
[in this window]
[in a new window]
 
FIG. 7. FP traces recorded by 2 electrodes at 3 stimulus strengths. A: Istim = 0.03 mA. B: Istim = 0.04 mA. C: Istim = 0.1 mA. As stimulus strength increases, discharge arrives faster at electrodes, but latency between the 2 arrival times (triangle t) does not change significantly.

Dependence on NMDA synaptic strength

MODEL. Increasing gNMDA results in a larger number of spikes and a higher velocity &cjs1726;D. However, the effect on &cjs1726;D is weaker in comparison with that of varying gAMPA. For example, between the minimal gAMPA for which a six-spike discharge is obtained (0.33 mS/cm2) and the minimal gAMPA for which an eight-spike discharge is obtained (0.42 mS/cm2), &cjs1726;D increases by 38% (Fig. 5A). In comparison, between the corresponding minimal gNMDA values (0.06 mS/cm2 and 0.27 mS/cm2, Fig. 5B), &cjs1726;D increases by only 14%. This difference stems from the different time scales of the two receptors. NMDA excitatory postsynaptic conductances (EPSCs) decay more slowly than AMPA EPSCs, and therefore their main effect is on the number of spikes.

If gAMPA < gAMPA,c, strong NMDA conductance enables propagation of the discharges in the model. The minimal gNMDA needed for propagation is reduced as gAMPA approaches gAMPA,c from below.

EXPERIMENT. We measured the velocity &cjs1726;D in disinhibited slices before and after applying APV (30-50 µM) to the bath. The FPs became briefer after the NMDA blockade. However, no significant change was found in &cjs1726;D (n = 4), supporting the model result.

Discharge propagation with synaptic depression

Because we found no experimental evidence for the existence of multiple discharge velocities, we searched for ways to modify the model. Synaptic depression reduces the magnitude of consecutive EPSCs in response to a train of presynaptic action potentials (Markram and Tsodyks 1996; Thomson et al. 1993). The magnitude of synaptic depression is frequency dependent and varies considerably between synapses. Introducing strong synaptic depression turns out to prevent the occurrence of multiple discharges for a specific parameter set. To demonstrate the effect of this process, we set the parameter kt = 1 ms-1 (Eq. 3), which corresponds to strong depression. With this value, the amount of free vesicles TGlu goes down by >60% immediately after the first spike. Five voltage time traces of neurons along the slice are shown in Fig. 8A. The basic characteristic of the discharge propagation as a traveling pulse is maintained: V(x, t) = &Vtilde;(x - &cjs1726;Dt). Time evolutions of synaptic and intrinsic auxiliary variables of one neuron within a network are plotted in Fig. 8B. Synaptic depression mainly affects the trajectory of the fast AMPA variable sAMPA, because consecutive action potentials generate weaker EPSCs, similar to the recordings of Markram and Tsodyks (1996). The AMPA peak response to the fourth spike is only 30% of the response to the first spike. Therefore only the first three spikes are strongly observed by the postsynaptic cell. The input synaptic field a cell receives, SinAMPA, is therefore briefer than in the case without synaptic depression; this helps to terminate the discharge. The slow NMDA variable sNMDA has a trajectory similar to that in the case without synaptic depression, but of a smaller amplitude. As before, the IK-slow activation variable z builds up during the activity and terminates it when it is strong enough.


View larger version (11K):
[in this window]
[in a new window]
 


View larger version (13K):
[in this window]
[in a new window]
 
FIG. 8. Propagation of discharges in cortical slice model. Activity is terminated by IK-slow and by synaptic depression (kt = 1.0 ms-1), gAMPA = 0.9 mS/cm2, gNMDA = 0.9 mS/cm2. A: voltage time traces of 5 cells along slice at same position as in Fig. 3. Every neurons fires 6 spikes during discharge. B: time courses of internal variables of neuron denoted by a. From top: total AMPA and NMDA synaptic conductances, SinAMPA and SinNMDA, that this neuron receives; presynaptic variable TGlu; AMPA and NMDA auxiliary variables, sAMPA and sNMDA, of postsynaptic synapses connecting this neuron to others; slow potassium activation variable z.

The dependencies of the number of spikes and the discharge velocity &cjs1726;D on the synaptic conductance strengths gAMPA and gNMDA are shown in Fig. 9, A and B. Only one discharge form is obtained here for a specific parameter set. The number of spikes increases with gAMPA and gNMDA. The velocity &cjs1726;D increases continuously with the synaptic conductance, even at points where the number of spikes jumps. The only exception is a slight discontinuity in &cjs1726;D, hardly seen in Fig. 9A, at the gAMPA value for which the number of spikes increases from two to three. The reason for that continuous dependence of &cjs1726;D on the synaptic conductances is the fact that the postsynaptic neuron is affected primarily by the first three presynaptic spikes. Note that the gAMPA,c is higher here than gAMPA,c without synaptic depression (Fig. 9), to compensate for the reduction in SinAMPA due to the depression.


View larger version (15K):
[in this window]
[in a new window]
 
FIG. 9. Number of spikes in propagating discharge (top) and discharge velocity &cjs1726;D (bottom) vs. AMPA synaptic conductance strength gAMPA (A) and NMDA synaptic conductance strength gNMDA (B). Activity is terminated by IK-slow and by synaptic depression (kt = 1.0 ms-1). In A, NMDA conductance is blocked (gNMDA = 0); In B, gAMPA = 0.9 mS/cm2. Arrow: gNMDA value of Fig. 8. For each parameter value there is a unique discharge shape.

Even more than in the case without synaptic depression, fast AMPA receptors affect mainly the discharge velocity, whereas slow NMDA receptors contribute mainly to increasing the number of spikes. For example, between the minimal gAMPA for which a three-spike discharge is obtained (0.57 mS/cm2) and the minimal gAMPA for which a five-spike discharge is obtained (1.19 mS/cm2), &cjs1726;D increases by 205% (Fig. 9A). In contrast, &cjs1726;D increases by only 11% between the minimal gNMDA needed to obtain a five-spike discharge (0.41 mS/cm2) and the minimal gNMDA needed to obtain a seven-spike discharge (1.11 mS/cm2).

AMPA decay rate

Two main factors control the discharge velocity: the kinetics of excitatory synapses and the time needed by resting neurons to integrate input EPSCs and reach the firing threshold. To discriminate between the two factors, we investigate the dependence of &cjs1726;D on the AMPA decay rate kr when NMDA receptors are blocked. First, we vary kr while keeping all the other parameters, including the AMPA conductance gAMPA, fixed. The velocity &cjs1726;D decreases continuously with kr (Fig. 10A) because higher kr means that the synaptic conductances are active for a briefer period and therefore the total input synaptic field is weaker.


View larger version (11K):
[in this window]
[in a new window]
 
FIG. 10. Dependence of discharge shape and velocity on AMPA decay rate kr. NMDA receptors are blocked, kt = 1 ms-1. A: gAMPA = 1.0 mS/cm2 is fixed. Velocity &cjs1726;D decreases with kr. B: we start from reference parameter values gAMPA = 1 mS/cm2 and kr = 0.2 ms-1 and vary gAMPA with kr to keep effective synaptic input cell receives (gAMPA,eff, Eq. 12) fixed. Velocity &cjs1726;D grows with kr and reaches saturated value at high kr.

To study the "true" effect of kr, we vary it while keeping fixed the effective synaptic conductance gAMPA,eff, defined as the following. The maximum synaptic input SinAMPA,max a cell receives as the discharge passes it is
<IT>S</IT><SUP>in</SUP><SUB>AMPA,max</SUB>= <LIM><OP>max</OP><LL><SUB><IT>t</IT></SUB></LL></LIM><IT>S</IT><SUP>in</SUP><SUB>AMPA</SUB>(<IT>t</IT>) (11)
Note that this magnitude is the same for all the neurons. The effective synaptic conductance is:
<IT>g</IT><SUB>AMPA,eff</SUB><IT>= g</IT><SUB>AMPA</SUB><IT>S</IT><SUP>in</SUP><SUB>AMPA,max</SUB> (12)
When kr is varied, we calculate iteratively the conductance gAMPA (for each kr value) for which the gAMPA,eff remains fixed. The dependence of &cjs1726;D on kr under this condition is presented in Fig. 10B. &cjs1726;D grows with kr and reaches saturation at high kr, where the synapses can be regarded as instantaneous. At the reference value (kr = 0.2 ms-1), &cjs1726;D is ~82% of the saturated value. Therefore the main factor that limits the discharge velocity is not the synaptic decay rate but the time needed for neuronal integration of EPSPs.

    DISCUSSION
Abstract
Introduction
Methods
Results
Discussion
References

Several conclusions emerge from our study. 1) The discharge propagates in our model as a traveling pulse. Cells with a distance x between them exhibit the same voltage time course with a time shift of x/&cjs1726;D, and this is a consequence of the geometry and the recruitment process. 2) Strong synaptic input leads to plateaulike voltage trajectories. 3) Discharges can propagate if gAMPA is larger than a threshold gAMPA,c. At the threshold, &cjs1726;D is finite and increases linearly with gAMPA at high gAMPA. 4) NMDA receptors do not contribute much to the propagation. 5) Synaptic depression prevents the appearance of multiple discharge forms. 6) The discharge velocity is mainly limited by the neuronal integration time, whereas synaptic kinetics plays a smaller role.

Architecture and dynamics

Neurons in our one-dimensional model of cortical slices are at rest without an external stimulus, whereas an initial shock that is strong enough produces a propagating discharge. A quiescent neuron in front of the discharge starts firing as a result of a recruitment process by active cells. A neuron ceases to fire because the slow potassium current builds up and, possibly, because of synaptic depression that helps in reducing the strength of input EPSCs. As a result, the network exhibits a discharge that moves with a constant velocity. Localized coupling (lambda  <<  L) is an essential condition for the existence of propagating discharges. In a small domain, the activity is almost instantly felt along the whole tissue.

Synchrony in our model means that cells with a distance x between them exhibit the same voltage time course with a time shift of x/&cjs1726;D. Dynamically, the discharge is a transient response to the initial shock. This synchronized activity is a result of the interplay between the excitatory input that induces firing and the slow potassium current IK-slow that prevents the firing from becoming a sustained state. Such interplay between network excitation and intrinsic adaptation has been shown to result in synchronous bursting in other models of hippocampal (Traub and Wong 1982) and neocortical (Hansel 1997) networks as well.

There are several factors that control and limit the velocity &cjs1726;D. The decay rate of AMPA, being finite, reduces the velocity by ~20% in comparison with the theoretical case for which the EPSCs are instantaneous but the effective synaptic coupling is maintained. The velocity is mainly limited by the time a neuron needs for integrating input EPSCs as the discharge approaches it and reaching the spiking threshold.

Single-cell and synaptic properties affect network behavior

We have related the various discharge shapes to the single-cell dynamics with IK-slow blocked and intact. In particular, the model with IK-slow blocked exhibits a high-voltage plateau under the application of a strong enough external current. Similarly, under the effect of strong synaptic input the neuron potential approaches such a plateau before being brought back to rest by IK-slow. Moderate synaptic input causes the neuron to fire at an almost constant amplitude.

IK-slow is responsible for adaptation in single RS cells (Gutnick and Crill 1995; Hansel and Sompolinsky 1996), because its activation time scale is longer than that of the other intrinsic ionic current but is not too slow. If the activation rate were, say, 10 times slower, a fast-slow separation of time scales at the single-cell level (e.g., Rinzel 1987) would be possible. In that case, the single-cell potential would fire a few spikes and then converge to the high plateau as a response to the injection of a step current pulse. Later, the slow potassium current would build up until the neuron potential goes back to rest. With the reference activation time we use (tau z = 75 ms, see APPENDIX), the high plateau is not achieved with moderate levels of applied current. Instead, the slow potassium current initially decreased the spike amplitude and later increases the interspike interval. At higher levels of applied current, IK-slow causes tonic firing, whereas blocking it switches the neuron to the high plateau.

At the network level, the strength of the slow potassium conductance gK-slow determines the discharge shape, in particular the number of spikes. However, the discharge velocity is determined mainly by the first few spikes, especially in the case of strong synaptic depression. Therefore gK-slow has only a minor effect on &cjs1726;D. Strong synaptic depression prevents the occurrence of more than one discharge form for a specific parameter set. The synaptic depression in our model is assumed to be presynaptic. Another model of synaptic depression suggests a postsynaptic mechanism (Tsodyks and Markram 1997). The network dynamics is not sensitive to the biophysical details of the mechanism.

To construct the model, we used approximations with respect to the single-cell and synaptic dynamics. The single cells were taken to be isopotential, and dendritic delays (Zador et al. 1995) are neglected. Anatomic data about the positions of intracortical boutons on the dendritic tree and multicompartmental models are needed to estimate the effects of those delays. All the cells are assumed to be regular spikers for two reasons. First, most of the cells in the slice are regular spikers (Gutnick and Crill 1995). Second, neocortical excitatory cells in a disinhibited slice exhibit a depolarized shift regardless of their intrinsic properties (Gutnick et al. 1982), and there is currently no reason to assume that the intrinsic bursting property in neocortex plays a major role in the activity.

Neglecting axonal delays is justified as long as the axonal velocity is much faster than &cjs1726;D. In our experiments, &cjs1726;D is ~16 cm/s. There is a large variation in the axonal conduction velocities measured in experiments. Estimations range from 0.15 to 32 m/s (Chervin et al. 1988; Mitani and Shimokouchi 1985; Murakoshi et al. 1993; Waxman and Swadlow 1977). In our experiments in the somatosensory cortex of adult mice (Gil and Amital 1996b), electrical stimulation of a slice 1 mm away from the recording electrode resulted in latencies to the monosynaptic EPSPs of <2 ms. Assuming a synaptic delay of 1 ms, the axonal velocity is estimated to be ~1 m/s (see also Haberly 1990) in olfactory cortex. Such conduction velocity is an order of magnitude larger that the discharge velocity, and its effects on the discharge propagation are therefore negligible.

Relations to wave propagation in other systems

Models of propagating activity in other brain areas have been built, some of them related to experiments performed on the same area. When inhibition is blocked, discharges propagate in longitudinal CA3 hippocampal slices with a velocity of ~15 cm/s, similar to our experimental results. A one-dimensional model of that network was analyzed by Miles et al. (1988). Single cells were taken to be endogenous bursters. As in the neocortical case, hippocampal discharges propagate by recruitment of new cells into the activity. The discharge velocity was calculated as a function of the footprint length lambda . Although in our case &cjs1726;D is linear with lambda , &cjs1726;D in the model of Miles et al. increases less than linearly with lambda  because of the effect of finite axonal conductance (0.5m/s in that case). Traub et al. (1993) studied the propagation of discharges along a transverse hippocampal slice. Their claim that the critical factor for propagation of the initial burst is the integration time over excitatory synaptic input is consistent with our result for neocortical models.

Activity can also propagate in thalamic slices, composed of excitatory thalamocortical neurons and inhibitory reticular thalamic neurons (Destexhe et al. 1996; Golomb et al. 1996; Kim et al. 1995). However, the mechanism of propagation is different there. Although in models of neocortical slices the discharge propagates continuously such that V(x, t) = &Vtilde;(x - &cjs1726;Dt), thalamic waves propagate in a lurching manner as a group of neurons from one type recruits a new group of neurons from the second type to the wave and vice versa. As a result of the lurching propagation, Golomb et al. (1996) showed that the wavefront velocity in thalamic slices depends approximately logarithmically on the synaptic strengths for an exponential footprint length. In contrast, we show here that in the neocortical (and hippocampal) networks &cjs1726;D grows linearly with the synaptic strength if gAMPA >> gAMPA,c. In models of thalamic slices, there is no minimal velocity, and &cjs1726;D can be close to zero. In neocortical slices, the discharge disappears at a finite velocity if gAMPA is reduced below gAMPA,c.

Propagation of activity was studied in simple models of neural networks (Ermentrout and McLoed 1993; Idiart and Abbott 1993). These works were the first to analyze rigorously propagation in integro-differential equations. In these models, the neuronal activity is represented by one variable and a moving wavefront separates a regime of high activity from a regime of low activity. Spikes and synaptic kinetics were not considered there.

Computational versus experimental results

Despite the strong heterogeneity in the intrinsic and synaptic properties of cortical slices (Chervin et al. 1988; Wadman and Gutnick 1993), several predictions of the model were confirmed experimentally. Propagation occurs if gAMPA is higher than a threshold value. At the threshold, &cjs1726;D is finite (not zero), and very slow propagation is not possible. Above it, &cjs1726;D increases with gAMPA. More quantitatively, prediction of a linear dependence of &cjs1726;D on gAMPA at strong gAMPA is hard to test experimentally, because the exact relation between the CNQX concentration and the gAMPA value in slices is not known.

Another prediction of the model is that NMDA conductance has only a small effect on &cjs1726;D. In our model, the NMDA rise time is considered to be fast. If a more elaborated NMDA model with a rise time of 10 ms (Jonas and Spruston 1994) is used, the NMDA effect on &cjs1726;D is even smaller. In the experiments, we have not found any significant change in &cjs1726;D when the NMDA receptors were completely blocked. The model predicts that NMDA conductance, being slow, will increase the DS duration and the number of spikes during the discharge. This result is also consistent with the experimental results of Baldino et al. (1986) in rat organotypic explant cultures. In contrast to the negligible effect of blocking NMDA on &cjs1726;D in cortex, Traub et al. (1993) showed experimentally and computationally that in longitudinal CA3 slices, blocking NMDA receptors slows &cjs1726;D by ~50%. The difference in the NMDA effect between cortex and hippocampus may be a result of a larger gNMDA and smaller gAMPA in the hippocampus. In addition, the resting potential of the hippocampal pyramidal cells reported by Traub et al. (1993) is -60 mV, ~14 mV more depolarized than in our model. As a result, the magnesium blockade of NMDA receptors in the hippocampal model is weaker, and the more effective NMDA synapses have a stronger effect on &cjs1726;D.

The model exhibits two functional regimes. Without synaptic depression, multiple discharge forms are observed at the same parameter values. The pattern is selected by the initial conditions, i.e., by the strength of the initial stimulus. Discharges with a larger number of spikes propagate faster; this effect is prominent when the number of spikes is small, but even then the velocity differences are on the order of a few percent. With strong synaptic depression, only one discharge form is observed for a specific parameter value. The dependence of &cjs1726;D on parameters such as the synaptic strengths is continuous. In the experiments, no significant dependence of &cjs1726;D on the stimulus strength has been observed.

Neurons in the model of a homogeneous network display the same firing pattern. In a real system, sparseness and heterogeneity lead to the appearance of various patterns displayed by neurons in the same slice (Connors and Gutnick 1984; Gutnick et al. 1982). We expect that a neuron that receives stronger excitation will fire in a "plateaulike" fashion, whereas a neuron that receives weaker excitation will fire action potentials with approximately fixed amplitude. High enough levels of sparseness and heterogeneity are expected to destroy the synchrony between neurons at the time scale of spikes but to preserve the synchrony at the time scale of the discharge. They may eliminate the existence of multiple discharge forms even in the case of no synaptic depression.

In the experiment, &cjs1726;D at the threshold for propagation is ~5 cm/s. In the model, &cjs1726;D at the threshold is found to be ~100 lambda /s for the parameter sets that we used. Combining these results together, lambda  is estimated to be on the order of 0.5 mm. This value depends mainly on the time a needed for the cell to integrate input EPSPs until its potential rises above the threshold of emitting a spike, and therefore the value is sensitive to the model's parameters, in particular the leak conductance. The fact that this estimated lambda  is consistent with the known length of local regions in cortex (Chervin et al. 1988; Wadman and Gutnick 1993) confirms that the parameter set we chose is reasonable.

Estimating the synaptic conductance needed for evoking a spike

Our experiments show that &cjs1726;D for AMPA conductance intact is 2-3 times faster than &cjs1726;D at threshold (Fig. 6). As shown in Fig. 9, such &cjs1726;D corresponds to gAMPA that is less than twice the threshold level gAMPA,c. Therefore we conclude that the amount of excitation gAMPA a typical cortical cell receives is ~2 times larger than the amount needed for propagation.

To understand what determines the threshold gAMPA,c, we consider a system of two cells, in which one cell excites the other with an AMPA synapse with a strength gAMPA. A single action potential is evoked in the presynaptic cell by injection of a step current pulse. If gAMPA is below a threshold level, the postsynaptic cell will exhibit a subthreshold depolarization; if gAMPA is above that threshold, the cell will fire a spike. The threshold for the parameter set of Fig. 9 is 0.17 mS/cm2.

Going back to Fig. 9, we find that the threshold for discharge propagation is gAMPA,c = 0.54 mS/cm2. The total synaptic efficacy from neurons left of a postsynaptic cell (subtracting the synaptic input from the right and the self-coupling term) is 0.24 mS/cm2. This means that to make discharge propagation possible with strong synaptic depression, the total synaptic conductance from the left should be only somewhat bigger than the synaptic input that evokes a spike in full synchrony. In the case of no synaptic depression the threshold can be even smaller because the postsynaptic neuron integrates the input conductance of several spikes.

Combining all these results, we estimate that the total amount of conductance gAMPA a typical cell receives in an intact slice during a paroxysmal discharge is only ~5 times larger than the AMPA conductance needed for raising the potential of a resting cell above its threshold. This ratio is expected to be higher, but of the same order of magnitude, where a two-dimensional model is concerned.

Although it is difficult to calculate the implication of this result for systems in vivo, because the ratio of intact synapses that remain after cutting the slice is hard to estimate, we believe that similar experimental and modeling techniques can be used to find the ratio of input conductance needed for spike generation in vivo.

    ACKNOWLEDGEMENTS

  We are grateful to E. Barkai, B. W. Connors, G. B. Ermentrout, I. Fleidervish, Z. Gil, M. J. Gutnick, D. Hansel, and J. Rinzel for helpful discussions and to D. Hansel and B. W. Connors for carefully reading this manuscript.

  This research was supported by a grant from the Israel Science Foundation.

    APPENDIX

Model equations and parameters

We use a Hodgkin-Huxley-like formulation for the excitatory cortical cells. The intrinsic ionic currents were modeled on the basis of existing voltage-clamp and current-clamp data. There are N cells, representing a spatial discretization of the continuous integro-differential Eq. 1. The position of the ith neuron is xi = i/N. We specify here our reference parameter set. These parameters are used throughout the paper unless stated otherwise. The maximal conductances (g's) and current kinetics are adjusted to fit action potential parameters and firing patterns of RS neocortical neurons (Connors and Gutnick 1990; Fleidervish et al. 1996; Huguenard et al. 1988; Yang et al. 1996).

Current balance equation
− <IT>I</IT><SUB>K-slow</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>, <IT>z</IT><SUB><IT>i</IT></SUB>) − <IT>I</IT><SUB>L</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>) − <IT>I</IT><SUB>AMPA</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>, {<IT>s</IT><SUB>AMPA,<IT>j</IT></SUB>})
− <IT>I</IT><SUB>NMDA</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>, {<IT>s</IT><SUB>NMDA,<IT>j</IT></SUB>}) + <IT>I</IT><SUB>app,<IT>i</IT></SUB> (A1)
C = 1µF/cm2, Iapp,j = 0.
<IT>C</IT><FR><NU>d<IT>V<SUB>i</SUB></IT></NU><DE>d<IT>t</IT></DE></FR> = −<IT>I</IT><SUB>Na</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>, <IT>h</IT><SUB><IT>i</IT></SUB>) − <IT>I</IT><SUB>NaP</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>) − <IT>I</IT><SUB>Kdr</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>, <IT>n</IT><SUB><IT>i</IT></SUB>) − <IT>I</IT><SUB>KA</SUB>(<IT>V</IT><SUB><IT>i</IT></SUB>, <IT>b</IT><SUB><IT>i</IT></SUB>)

Intrinsic currents SODIUM CURRENT INa.
<FR><NU>d<IT>h</IT></NU><DE>d<IT>t</IT></DE></FR>= [<IT>h</IT><SUB>∞</SUB>(<IT>V</IT>) − <IT>h</IT>]/τ<SUB><IT>h</IT></SUB>(<IT>V</IT>) (A3)
<IT>m</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>m</IT></SUB>)/σ<SUB><IT>m</IT></SUB>]}<SUP>−1</SUP> (A4)
<IT>h</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>h</IT></SUB>)/σ<SUB><IT>h</IT></SUB>]}<SUP>−1</SUP> (A5)
τ<SUB><IT>h</IT></SUB>(<IT>V</IT>) = 0.37 + 2.78 × {1 + exp[−(<IT>V</IT>− θ<SUB><IT>ht</IT></SUB>)/σ<SUB><IT>ht</IT></SUB>]}<SUP>−1</SUP> (A6)
gNa = 24 mS/cm2, VNa = 55 mV, theta m = -30 mV, sigma m = 9.5 mV, theta h = -53 mV, sigma h = -7 mV, theta ht -40.5 mV, sigma ht = -6 mV. Our INa model is based on the experimental results of Fleidervish et al. (1996). As in other models of cortical cells (e.g., Hansel and Sompolinsky 1996; Traub et al. 1991; Wang 1993), the function hinfinity (V) is shifted to the right on the voltage axis to enable the model to replicate the current-clamp current injection experiments. Gating kinetics has been adjusted to 36°C.
<IT>I</IT><SUB>Na</SUB>(<IT>V</IT>, <IT>h</IT>) = <IT>g</IT><SUB>Na</SUB><IT>m</IT><SUP>3</SUP><SUB>∞</SUB>(<IT>V</IT>)<IT>h</IT>(<IT>V − V</IT><SUB>Na</SUB>) (A2)

PERSISTENT SODIUM CURRENT INaP.
<IT>p</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>p</IT></SUB>)/σ<SUB><IT>p</IT></SUB>]}<SUP>−1</SUP> (A8)
gNaP = 0.07 mS/cm2, theta p = -40 mV, sigma p = 5.0 mV (Alzheimer et al. 1993; Fleidervish and Gutnick 1996).
<IT>I</IT><SUB>NaP</SUB>(<IT>V</IT>, <IT>h</IT>) = <IT>g</IT><SUB>NaP</SUB><IT>P</IT><SUB>∞</SUB>(<IT>V</IT>)(<IT>V − V</IT><SUB>Na</SUB>) (A7)

DELAYED RECTIFIER POTASSIUM CURRENT IKdr.
<FR><NU>d<IT>n</IT></NU><DE>d<IT>t</IT></DE></FR>= [<IT>n</IT><SUB>∞</SUB>(<IT>V</IT>) − <IT>h</IT>]/τ<SUB><IT>n</IT></SUB>(<IT>V</IT>) (A10)
τ<SUB><IT>n</IT></SUB>(<IT>V</IT>) = 0.37 + 1.85 × {1 + exp[−(<IT>V</IT>− θ<SUB><IT>nt</IT></SUB>)/σ<SUB><IT>nt</IT></SUB>]}<SUP>−1</SUP> (A11)
<IT>n</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>n</IT></SUB>)/σ<SUB><IT>n</IT></SUB>]}<SUP>−1</SUP> (A12)
gKdr = 3 mS/cm2, VK = -90 mV, theta n = -30 mV, sigma n = 10 mV, theta nt = -27 mV, sigma nt = -15 mV. (Fleidervish et al. 1996; Wang 1993).
<IT>I</IT><SUB>Kdr</SUB>(<IT>V</IT>, <IT>n</IT>) = <IT>g</IT><SUB>Kdr</SUB><IT>n</IT><SUP>4</SUP>(<IT>V − V</IT><SUB>K</SUB>) (A9)

A-TYPE POTASSIUM CURRENT IKA.
<FR><NU>d<IT>b</IT></NU><DE>d<IT>t</IT></DE></FR>= [<IT>b</IT><SUB>∞</SUB>(<IT>V</IT>) − <IT>b</IT>]/τ<SUB><IT>b</IT></SUB> (A14)
<IT>a</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>a</IT></SUB>)/σ<SUB><IT>a</IT></SUB>]}<SUP>−1</SUP> (A15)
<IT>b</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>b</IT></SUB>)/σ<SUB><IT>b</IT></SUB>]}<SUP>−1</SUP> (A16)
gKA = 1.4 mS/cm2, theta a = -50 mV, sigma a = 20 mV, theta b = -80 mV, sigma b = -6 mV, tau b = 15 ms (Hansel and Sompolinsky 1996; Rush and Rinzel 1995).
<IT>I</IT><SUB>KA</SUB>(<IT>V</IT>, <IT>b</IT>) = <IT>g</IT><SUB>KA</SUB><IT>a</IT><SUP>3</SUP><SUB>∞</SUB><IT>b</IT>(<IT>V − V</IT><SUB>K</SUB>) (A13)

SLOW POTASSIUM CURRENT IK-slow.
<FR><NU>d<IT>z</IT></NU><DE>d<IT>t</IT></DE></FR>= [<IT>z</IT><SUB>∞</SUB>(<IT>V</IT>) − <IT>z</IT>]/τ<SUB><IT>z</IT></SUB> (A18)
<IT>z</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>z</IT></SUB>)/σ<SUB><IT>z</IT></SUB>]}<SUP>−1</SUP> (A19)
gK-slow = 1 mS/cm2, theta z = -39 mV, sigma z = 5 mV, tau z = 75 ms (Hansel and Sompolinsky 1996).
<IT>I</IT><SUB>K-slow</SUB>(<IT>V</IT>, <IT>z</IT>) = <IT>g</IT><SUB>K-slow</SUB><IT>z</IT>(<IT>V − V</IT><SUB>K</SUB>) (A17)

LEAK CURRENTS IL. gL = 0.02 mS/cm2, VL = -70 mV.
<IT>I</IT><SUB>L</SUB><IT>= g</IT><SUB>L</SUB>(<IT>V − V</IT><SUB>L</SUB>) (A20)

Synaptic currents AMPA CURRENT IAMPA.
<FR><NU>d<IT>s</IT><SUB>AMPA,<IT>j</IT></SUB></NU><DE>d<IT>t</IT></DE></FR>= <IT>k<SUB>f</SUB>T</IT><SUB>Glu</SUB><IT>s</IT><SUB>∞</SUB>(<IT>V</IT><SUB><IT>j</IT></SUB>)(1 − <IT>s</IT><SUB>AMPA,<IT>j</IT></SUB>) − <IT>k</IT><SUB><IT>r</IT></SUB><IT>s</IT><SUB>AMPA,<IT>j</IT></SUB> (A22)
<FR><NU>d<IT>T</IT><SUB>Glu,<IT>j</IT></SUB></NU><DE>d<IT>t</IT></DE></FR>= −<IT>k<SUB>t</SUB>s</IT><SUB>∞</SUB>(<IT>V</IT><SUB><IT>j</IT></SUB>)<IT>T</IT><SUB>Glu,<IT>j</IT></SUB><IT>+ k</IT><SUB><IT>&cjs1726;</IT></SUB>(1 − <IT>T</IT><SUB>Glu,<IT>j</IT></SUB>) (A23)
<IT>s</IT><SUB>∞</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>s</IT></SUB>)/σ<SUB><IT>s</IT></SUB>]}<SUP>−1</SUP> (A24)
theta s = -20 mV, sigma s = 2 mV, kf = 1 ms-1, kr = 0.2 ms-1, kt = 0 or 1 ms-1, k&cjs1726; = 0.001 ms-1, gAMPA = 0.9 mS/cm2, VGlu = 0 mV (Golomb et al. 1994, 1996; Jonas and Spruston 1994; Stern et al. 1992).
<IT>I</IT><SUB>AMPA</SUB>(<IT>V</IT>, {<IT>s</IT><SUB>AMPA,<IT>j</IT></SUB>}) = <IT>g</IT><SUB>AMPA</SUB>(<IT>V − V</IT><SUB>Glu</SUB>) <LIM><OP>∑</OP><LL><SUB><IT>j</IT>=1</SUB></LL><UL><SUP><IT>N</IT></SUP></UL></LIM><IT><A><AC>w</AC><AC>˜</AC></A></IT>(<IT>i − j</IT>)<IT>s</IT><SUB>AMPA,<IT>j</IT></SUB> (A21)

NMDA CURRENT INMDA.
<FR><NU>d<IT>s</IT><SUB>NMDA,<IT>j</IT></SUB></NU><DE>d<IT>t</IT></DE></FR>= <IT>k<SUB>f<UP>N</UP></SUB>T</IT><SUB>Glu</SUB><IT>s</IT><SUB>∞</SUB>(<IT>V</IT><SUB><IT>j</IT></SUB>)(1 − <IT>s</IT><SUB>NMDA,<IT>j</IT></SUB>) − <IT>k</IT><SUB><IT>r</IT>N</SUB><IT>s</IT><SUB>NMDA,<IT>j</IT></SUB> (A26)
<IT>f</IT><SUB>NMDA</SUB>(<IT>V</IT>) = {1 + exp[−(<IT>V</IT>− θ<SUB><IT>fN</IT></SUB>)/σ<SUB><IT>f</IT>N</SUB>]}<SUP>−1</SUP> (A27)
kf = 1 ms-1, ktau N = 0.0067 ms-1, theta fN = -25 mV, sigma fN = 12.5 mV, gNMDA = 0.9 mS/cm2, VGlu = 0 mV (Destexhe et al. 1994; Jonas and Spruston 1994; Stern et al. 1992).
<IT>I</IT><SUB>NMDA</SUB>(<IT>V</IT>, {<IT>s</IT><SUB>NMDA,<IT>j</IT></SUB>}) = <IT>g</IT><SUB>NMDA</SUB><IT>f</IT><SUB>NMDA</SUB>(<IT>V</IT>)(<IT>V − V</IT><SUB>Glu</SUB>) <LIM><OP>∑</OP><LL><SUB><IT>j</IT>=1</SUB></LL><UL><SUP><IT>N</IT></SUP></UL></LIM><IT><A><AC>w</AC><AC>˜</AC></A></IT>(<IT>i − j</IT>)<IT>s</IT><SUB>NMDA,<IT>j</IT></SUB> (A25)

Architecture The synaptic coupling term rho w(x) (Eq. 8) is discretized and yields the discrete function &wtilde;(j), where
<IT><A><AC>w</AC><AC>˜</AC></A></IT>(<IT>j</IT>) = tanh [<IT>L</IT>/(2λ<IT>N</IT>)] exp[−‖<IT>j</IT>‖<IT>L</IT>/(λ<IT>N</IT>)] (A28)
The normalization constant tanh [L/(2lambda N)] is determined by the constraint <LIM><OP>Σ</OP><LL><SUB><IT>j</IT>=−∞</SUB></LL><UL>∞</UL></LIM> w(j) = 1. lambda  = 0.03125, N = 256, L = 1.

The simulations reported in this paper were performed with the use of the fourth-order Runge-Kutta method with time step Delta t = 0.03 ms. We simulated networks with N between 128 and 512. Increasing N has little effect on the results as long as lambda N is large enough (say >= 4). Synaptic fields (Eq. A21 and A25) were computed with the use of fast Fourier transform (Golomb et al. 1996; Press et al. 1992). Bifurcation diagrams (Fig. 2, B and D) were calculated with the use of the program AUTO (Doedel 1981) incorporated into software package XPP, written by G. B. Ermentrout (Rinzel and Ermentrout 1997).

    FOOTNOTES

1   The propagating activity is called "discharge" here. In the dynamics literature, propagating activity that is localized in both space and time is often referred to as a "pulse" (Cross and Hohenberg 1993). 2   In the mathematical sense, traveling pulses only exist in an infinite domain; the discharges that are seen here are like traveling pulses in large domains.

  Address for reprint requests: D. Golomb, Dept. of Physiology, Faculty of Health Sciences, PO Box 653, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel.

  Received 11 January 1997; accepted in final form 7 May 1997.

    REFERENCES
Abstract
Introduction
Methods
Results
Discussion
References

0022-3077/97 $5.00 Copyright ©1997 The American Physiological Society