Abstract
Short-term facilitation is a widely observed form of synaptic enhancement that is not well understood. Although presynaptic calcium has long been implicated in this process, its role is unclear, particularly at synapses in the mammalian brain. We tested the role of presynaptic residual free calcium ([Ca]res) in facilitation of synapses between granule cells and Purkinje cells in rat cerebellar slices. Paired-pulse facilitation of synaptic currents resulted in an approximately 2.5-fold enhancement that decayed with a time constant of ∼200 msec, as assessed by voltage-clamp recordings. Measurements of [Ca]res using fluorescent calcium-sensitive indicators revealed that [Ca]resdecayed more rapidly than did facilitation. Manipulation of [Ca]res dynamics by introducing EGTA into presynaptic terminals sped the decays of [Ca]res and facilitation in a dose-dependent manner. When [Ca]res was reduced to a brief impulse lasting several milliseconds, facilitation was still present, although reduced in amplitude and duration. Facilitation decayed with an intrinsic time constant of ∼40 msec. These results suggest that facilitation at this synapse is produced by a calcium-driven process with a high affinity and a slow effective off-rate. A combination of [Ca]res dynamics and the properties of a calcium-driven reaction determine the time course and amplitude of facilitation.
- synaptic transmission
- parallel fiber
- presynaptic terminal
- short-term enhancement
- calcium
- magnesium green
- cerebellum
At many synapses, when two action potentials depolarize a presynaptic bouton in rapid succession, the second action potential releases more neurotransmitter than the first (Magleby, 1987;Zucker, 1989). This short-term enhancement of release that persists for tens to hundreds of milliseconds after a conditioning pulse is known as facilitation. Facilitation may play an important role in the synaptic coding of information contained in the temporal pattern of spike trains (Buonomano and Merzenich, 1995). In addition, characterizing facilitation has become an important way of studying synaptic transmission in the mammalian brain, and changes in the properties of paired-pulse facilitation are often taken to indicate a presynaptic change in neurotransmitter release (McNaughton, 1982; Schulz et al., 1994).
We are particularly interested in elucidating the mechanism of facilitation of synapses in the mammalian CNS. Residual free calcium ([Ca]res) that persists in the presynaptic terminal after termination of release has long been implicated in facilitation at synapses with large presynaptic terminals, such as the neuromuscular junction (NMJ) and the squid giant synapse (Katz and Miledi, 1968;Magleby, 1987; Bittner, 1989; Zucker, 1989; Liu and Stanley, 1995). Synapses in the mammalian CNS, at the NMJ, and in squid display facilitation with many common properties, which suggests that similar mechanisms are used by these synapses, but diversity in [Ca]res dynamics suggests that there may be important differences (Delaney et al., 1989; Swandulla et al., 1991; Delaney and Tank, 1994; Regehr et al., 1994; Regehr and Atluri, 1995; Feller et al., 1996). In large presynaptic terminals, buffered calcium diffusion plays an important role in determining [Ca]res at the release site on the time scale of facilitation (Yamada and Zucker, 1992; Winslow et al., 1994). In contrast, in small presynaptic boutons (∼1 μm in diameter) typical of the mammalian CNS, [Ca]res dynamics is controlled primarily by calcium extrusion (Regehr and Atluri, 1995). This suggested to us the possibility that the action of calcium in facilitation and the properties of calcium-binding molecules involved in facilitation might well be tailored to match the fast, extrusion-limited [Ca]res dynamics characteristic of small presynaptic terminals associated with synapses in the mammalian brain.
The rapid spatial equilibration of calcium in these small boutons confers an important technical advantage for studying the role of calcium in facilitation: the ability to measure [Ca]resdynamics on a tens of milliseconds time scale using calcium-sensitive fluorescent indicators. This has not been possible at the NMJ and squid synapses, because direct measurement of the time course of [Ca]res cannot be made until calcium has spatially equilibrated by buffered diffusion. A conservative estimate of the equilibration time for calcium in a presynaptic bouton of radiusr can be obtained using the equation for the characteristic diffusion time of τch = r2/6D, and using a rather small effective diffusion coefficient for calcium, D, of 1 × 10−7 cm2/sec. We estimate the equilibration time of calcium to be ∼3 msec for a 0.8-μm-diameter bouton, such as at the granule cell→Purkinje cell synapse, and 200 msec for a 7-μm-diameter terminal, such as at the crayfish NMJ.
Using modeling and indirect measures of [Ca]res dynamics, several hypotheses concerning the role of calcium in facilitation have emerged from the study of the NMJ and squid synapses (Charlton et al., 1982; Stanley, 1986; Bain and Quastel, 1992; Tanabe and Kijima, 1992;Yamada and Zucker, 1992; Blundon et al., 1993; Delaney and Tank, 1994;Kamiya and Zucker, 1994; Winslow et al., 1994; Liu and Stanley, 1995). Although there is agreement that facilitation is produced by calcium, there is uncertainty concerning the relative importance of calcium dynamics and the effective kinetics of calcium-activated processes in determining the time course of facilitation.
According to the simplest form of the residual calcium hypothesis, in facilitation calcium acts at the same low-affinity calcium binding sites involved in triggering release. In this model, the very high (>100 μm) submembrane calcium levels that trigger release dissipate by diffusion into the terminal, and calcium unbinds rapidly from the low-affinity trigger sites to terminate release of neurotransmitter (Zucker, 1989). Calcium entering during a subsequent action potential summates with [Ca]res, enhancing release. If release were just proportional to the calcium concentration detected at the trigger, then the time course of facilitation should track the time course of decay of [Ca]res; however, because of the power relation between calcium concentration and release (Dodge and Rahamimoff, 1967; Heidelberger et al., 1994; Heinemann et al., 1994), the decay of facilitation predicted by such a mechanism should be faster than the decay of [Ca]res. Calcium may also remain bound to these sites, as in facilitation of calcium-dependent bioluminescence in Obelia (Naranjo et al., 1994).
Another possibility is that [Ca]res acts at a site with a higher calcium affinity than those involved in triggering release (Stanley, 1986; Delaney et al., 1989; Swandulla et al., 1991; Yamada and Zucker, 1992; Liu and Stanley, 1995). Some studies suggest that this site has rapid kinetics, with the time course of facilitation determined primarily by the time course of [Ca]res(Kamiya and Zucker, 1994). Others favor a model with slow effective kinetics, so that the main determinant of the time course of facilitation is the kinetics of the calcium unbinding reaction or some consequence thereof (Blundon et al., 1993; Winslow et al., 1994).
Less is known about the role of calcium in setting the time course of facilitation at synapses with small presynaptic boutons. In frog tectal synapses, the magnitude but not the time course has been shown to be influenced by changes in [Ca]res (Feller et al., 1996). Studies of synapses in the mammalian CNS have raised many questions regarding the role of [Ca]res in facilitation (Wu and Saggau, 1994) (see Discussion).
Here we examine the issue of what sets the time course of paired-pulse facilitation at the synapse between cerebellar granule cells and Purkinje cells. By using fast low-affinity calcium indicators, we are able to track changes in [Ca]res on a time scale appropriate for the study of facilitation (Regehr and Atluri, 1995). In addition, we manipulated [Ca]restransients to test further the role of calcium in facilitation. These studies establish that at this synapse the time course of facilitation is dictated in part by [Ca]res dynamics and in part by the properties of a calcium-driven process with slow effective kinetics that is activated by modest calcium levels.
MATERIALS AND METHODS
Transverse slices (300 μm thick) were cut from the cerebellar vermis of 11- to 14-d-old rats, as described previously. Experiments were conducted at 24 ± 0.5°C, except for the experiments associated with Figure 3B, which were performed at 34 ± 0.5°C. The external solution (2 ml/min flow rate) consisted of (in mm): 125 NaCl, 2.5 KCl, 2 CaCl2, 1 MgCl2, 26 NaHCO3, 1.25 NaH2PO4, 25 glucose, and 0.02 bicuculline, bubbled with 95% O2 and 5% CO2.
For EGTA wash-in experiments, a 100 mm stock solution of EGTA-AM in dimethyl sulfoxide (DMSO) was aliquoted and frozen. Immediately before use, aliquots were diluted in external saline to final EGTA concentrations (in μm) of 100, 20, or 1, and DMSO concentrations (%) of 0.1, 0.02, and 0.001. Wash-ins were 15 min long, with flow rates of 2 ml/min. To control for the effects of DMSO on facilitation, control experiments were performed with wash-in of DMSO alone. DMSO had no significant effect on the time course of facilitation in these experiments. Some experiments were also performed in which DMSO was present in the external solution, bathing the slice both before and after application of EGTA.
Detecting presynaptic calcium transients. Parallel fibers, made up of granule cell axons and presynaptic terminals, were labeled by local application of a solution containing the membrane-permeant forms of calcium indicators, as described previously (Regehr and Tank, 1991; Regehr and Atluri, 1995). The loading time was 3–5 min for mag-fura-5 AM, 3 min for fura-2 AM, and 8–10 min for magnesium green-AM (all from Molecular Probes, Eugene, OR). Experiments commenced 2 hr after dye-loading. Parallel fiber tracts were stimulated extracellularly with an electrode placed in the molecular layer near the fill site. Fluorescence changes were measured in a 150-μm-diameter spot 400–700 μm away from the stimulus site.
The filter set that was used for fura-2 and mag-fura-5 was a 380 HT 15 excitation, a 430DCLP02 dichroic, and a 510WB40 emission filter (Omega Optical, Brattleboro, VT), and for magnesium green, a 450–490 excitation, an FT510 dichroic, and an LP520 emission filter (Zeiss, Thornwood, NY). Fluorescence was detected with a photomultiplier tube (Hamamatsu HC124-06MOD; Hamamatsu, Bridgewater, NJ). Because calcium increases produce a decrease in fluorescence of fura-2 and mag-fura-5 with 380 nm excitation, ΔF/F traces in Figure 2have been inverted for these dyes.
Correction of fura-2 traces. As a consequence of its high affinity, fura-2 is known to distort the large calcium transients that occur in skeletal muscle (Konishi et al., 1988) and in small structures, such as parallel fiber presynaptic terminals (Regehr and Atluri, 1995). The degree of fura-2 saturation has been used to estimate the change in [Ca]res produced by a single action potential (Regehr and Atluri, 1995; Sabatini and Regehr, 1995). In a straightforward extension of the reasoning of these previous studies, it is also possible to correct for the distortion of the time course of ΔF/F decay. The curve relating the number of spikes in a train, Nspikes, to peak ΔF/F changes produced by 100 Hz trains, (ΔF/F)peak, was fit to a function of the form: Equation 1where K* and (ΔF/F)max are constants. The measured ΔF/F signal was then corrected for saturation of the indicator using the constants determined above and the following equation: Equation 2where c is a constant.
The similarity of the time course of calcium decay for the corrected fura-2 transients and the low-affinity calcium indicators (see Fig. 2) supports the view that the time course of [Ca]res decay is reported faithfully by mag-fura-5 and magnesium green. Although the procedure described above corrects for the distortion produced by steady-state saturation of the indicator, it does not correct for dye kinetics (Regehr and Atluri, 1995).
Measuring synaptic currents. Whole-cell recordings of Purkinje neurons were obtained as described previously (Regehr and Mintz, 1994), using 0.9–1.5 MΩ glass pipettes containing an internal solution of (in mm): 35 CsF, 100 CsCl, 10 EGTA, 10 HEPES, and 0.1 D600, pH 7.3 with CsOH. The access resistance (<5 MΩ after series resistance compensation) and leak current (0 pA to −200 pA) were monitored continuously.
Parallel fibers were stimulated extracellularly with a bipolar electrode placed in the molecular layer several hundred micrometers from the recording electrode. The resulting excitatory postsynaptic current (EPSC) decayed with a time constant of 5–7 msec at 24°C and 3–5 msec at 34°C. Low stimulus intensities were used to keep synaptic currents small and the resulting voltage errors arising from uncompensated series resistance below 4 mV.
Facilitation was computed for each trial according to the equation 100 × (A2 − A1)/A1, whereA1 and A2 are the peak amplitudes of the unfacilitated and facilitated synaptic currents, respectively. Similar results were obtained when averages of trials were used. At short interpulse intervals (≤30 msec), the test EPSC was superimposed on the conditioning EPSC. For these trials, the waveform of the conditioning EPSC was subtracted from the superimposed test EPSC before determination of the maximum amplitude of the test EPSC.
Data acquisition and analysis. Outputs of the photomultiplier tube and Axopatch 200A were filtered at 1 kHz with a model 900C9L8L 8-pole Bessel filter (Frequency Devices, Haverhill, MA) and digitized with a 16-bit D/A converter (Instrutech, Great Neck, NY), Pulse Control software (Herrington and Bookman, 1995), and an Apple Macintosh Centris 650 computer. Analysis was done on- and off-line with Igor Pro software (Wavemetrics, Lake Oswego, OR). Exponential fits to facilitation were performed between 10 and 1000 msec. The amplitude of facilitation was calculated from this fit for Δt = 10 msec.
RESULTS
Facilitation and presynaptic calcium for control conditions
The synapse between cerebellar granule cells and Purkinje cells exhibits robust paired-pulse facilitation, as shown in Figure1. Synaptic currents evoked by extracellular stimulation of parallel fibers were measured with whole-cell recordings from the soma of cerebellar Purkinje cells (Perkel et al., 1990; Llano et al., 1991; Mintz et al., 1995). As shown in the inset of Figure 1, stimulation with two pulses separated by 50 msec resulted in a facilitation of ∼180%, as calculated by the formula 100 × (A2–A1)/A1, where A1 is the amplitude of the unfacilitated synaptic current and A2 is the amplitude of the facilitated current. In Figure 1, the percentage facilitation as a function of interpulse interval (Δt) can be approximated by a single exponential decay with an amplitude of 194% and τ = 190 msec (solid curve). For a series of such experiments (n = 15), the amplitude of facilitation was 153 ± 11% (mean ± SEM), and the time constant of decay was 203 ± 18 msec.
To investigate the role of [Ca]res in facilitation, we sought to obtain a measure of [Ca]res in parallel fiber presynaptic terminals on the time scale of facilitation, namely in the 1 sec interval after stimulation. Parallel fiber tracts were labeled with calcium indicators as described in Materials and Methods. We have shown previously that indicators with low calcium affinities, such as mag-fura-5 (Claflin et al., 1994) and magnesium green (Zhao et al., 1996), faithfully follow the large and rapid [Ca]restransients produced by parallel fiber activation (Regehr and Atluri, 1995; Sabatini and Regehr, 1995). Like many of the other indicators that were developed originally to detect magnesium, these fluorophores have proven to be very useful in detecting large calcium transients without significant interference from magnesium transients (Konishi et al., 1991; Regehr and Atluri, 1995; Zhao et al., 1996). As shown in Figure 2, the time courses of ΔF/F signals for mag-fura-5 and magnesium green are similar, but fura-2 ΔF/F transients decay more slowly. The slower ΔF/F decay with the high-affinity dye fura-2 is expected for large [Ca]restransients, which begin to saturate the dye (Regehr and Atluri, 1995;Feller et al., 1996). It is possible, however, to obtain a better estimate of the time course of [Ca]res transients by correcting for most of the distortion of fura-2 (see Materials and Methods), as shown in Figure 2C. After applying this procedure, the decay of the corrected fura-2 signal resembles more closely those of the low-affinity dyes. For the rest of the studies described here we used magnesium green to detect [Ca]res. It is less sensitive to calcium (KD ∼ 5 μm) than is fura-2 (KD ∼ 200 nm) (Grynkiewicz et al., 1985) but more sensitive than is mag-fura-5 (KD ∼ 23 μm) (Delbono and Stefani, 1993); this allows it to report calcium changes accurately, but with a better signal-to-noise ratio than mag-fura-5.
On the basis of anatomical considerations, such ΔF/F signals obtained using low-affinity indicators provide a good measure of calcium transients in parallel fiber boutons that synapse onto Purkinje cells, with minimal contamination from other structures. Parallel fibers make a series ofen passant synapses as they course through the molecular layer. It is estimated that ∼94% of these synapses are onto Purkinje cell dendrites (Palkovits et al., 1971). Presynaptic contacts are made via closely spaced varicosities, which average 0.8 μm in length and 0.18 μm3 in volume (Palay and Chan-Palay, 1974). They are separated by thin segments of axon 0.15 μm in diameter, with an average length of 1.7 μm and an average volume of ∼0.03 μm3. We estimate that these presynaptic varicosities comprise ∼86% of the volume of a parallel fiber, with the remainder made up of thin axon.
We also considered the possibility of differential distribution of calcium channels and extrusion mechanisms on the axon relative to the presynaptic terminals. If calcium channels and extrusion mechanisms were both located preferentially on axonal segments, then axonal signals could contribute a rapid component to the fluorescence transient, which would make the fluorescence changes we detect difficult to interpret. This is unlikely, because calcium channels and extrusion mechanisms are generally distributed preferentially in presynaptic terminals relative to axons (Cohen et al., 1991; Delaney et al., 1991; Kortje et al., 1991; Luther et al., 1992; Westenbroek et al., 1992, 1995; Smith et al., 1993). This is also true in parallel fibers for the α1A calcium channel subunit (Westenbroek et al., 1995). Thus axonal signals are unlikely to contribute significantly to the fluorescence changes we detect.
To compare the time course of calcium decay to the time course of facilitation, it seemed most natural to compare fits to exponential decays. Fits of the fluorescence transients to single exponential decays provide a crude measure of the decay time of calcium, with a control time constant of ∼90 msec over the range 10 msec to 1 sec; however, the decay of calcium was not very well approximated by a single exponential, and the decay could not be fit reliably by the sum of two distinct exponentials. Half-decay times (t1/2), the time for [Ca]res to decay to 50% of peak elevation, are useful measures of the decay time of calcium; however, they do not fully describe calcium dynamics and cannot be directly related to time constants of exponential decay. Therefore, to avoid assumptions about the form of the decays, we directly compared the time course of facilitation and residual calcium simply by plotting them on the same graph, as in Figure3.
After single stimuli, [Ca]res decayed more rapidly than did facilitation. Figure 3A shows that at 24°C the decay of facilitation, fit by a single exponential with τfac = 184 msec (solid curve), lags behind the decay of [Ca]res. As shown in Figure 3B, raising the temperature to 34°C speeds the decay of both facilitation and [Ca]res, but here too facilitation outlasts the elevation of [Ca]res. All subsequent experiments were conducted at 24°C.
These results have important implications for the role of calcium in facilitation. They are inconsistent with the simplest form of the residual calcium hypothesis for paired-pulse facilitation, in which the calcium sensor that triggers release also mediates facilitation. That model predicts that facilitation decays at the same rate as or faster than [Ca]res. Instead, we see that calcium decays more rapidly than facilitation.
Altering presynaptic calcium transients with EGTA
To pursue further the role of [Ca]res in paired-pulse facilitation, we decided to manipulate the time course of [Ca]res and test its effect on facilitation. As described in Materials and Methods and the (see Fig. 10), by introducing EGTA into presynaptic boutons we were able to speed the decay of [Ca]res (Adler et al., 1991; Winslow et al., 1994; Feller et al., 1996). Because EGTA is a slow calcium chelator, however, it had relatively small effects on the large, brief calcium increases that trigger release.
Figure 4 shows a series of experiments in which successively larger concentrations of EGTA-AM, the membrane-permeant form of EGTA, were bath-applied for 15 min to parallel fibers. The time courses and amplitudes of [Ca]res elevations were monitored with magnesium green (Fig. 4, left panels), with averages of fluorescence transients before and after application of EGTA-AM shown in the right panels; 1 μm EGTA-AM accelerated the t1/2 of [Ca]restransients from 36 msec to 16 msec with no significant effect on peak [Ca]res levels (Fig. 4A). Increasing the concentration of presynaptic EGTA with 20 μm EGTA-AM led to a more pronounced acceleration of the t1/2 of [Ca]res (32 msec to 6 msec), with little if any effect on peak [Ca]res levels (Fig. 4B). Higher concentrations of EGTA-AM (100 μm) sped thet1/2 from 33 msec to 2 msec (Fig. 4C) and reduced peak ΔF/F values by 33%.
A summary of the effects of 15 min applications of EGTA-AM on parallel fiber [Ca]res transients is shown in Figure5. The t1/2 for control, 1 μm, 20 μm, and 100 μm were 39 ± 1 msec (n = 20), 23 ± 2 msec (n = 5), 10 ± 1 msec (n = 5), and 2 ± 0.2 msec (n = 5), respectively; 20 μm EGTA-AM caused a small 8 ± 4% decrease in peak amplitude, whereas 100 μm EGTA-AM decreased peak amplitude by 45 ± 4%. These experiments demonstrate that EGTA accelerates the decay of [Ca]res transients in a dose-dependent manner and at higher concentrations begins to attenuate the peak of the [Ca]res transient.
Effect of rapid presynaptic calcium dynamics on facilitation
We then assessed the effect of altered calcium dynamics on facilitation. Parallel fibers were stimulated every 15 sec while the interpulse interval Δt was varied and conditioning and test EPSCs were recorded. Figure 6 shows facilitation before and after 15 min applications of successively higher concentrations of EGTA-AM. A solution of 0 EGTA-AM in 0.1% DMSO did not significantly change the time course of facilitation (Fig.6A). Application of 1 μm EGTA-AM accelerated the decay of facilitation without greatly affecting the amplitude of facilitation (Fig. 6B). Higher concentrations of EGTA-AM had more pronounced effects on the rate of decay of facilitation and also decreased the amplitude of facilitation, as shown in Figure6C for 20 μm EGTA-AM and in Figure6D for 100 μm EGTA-AM.
Although EGTA-AM will load all structures within the slice, including presynaptic terminals and Purkinje cells, it is highly likely that the effect of EGTA-AM on facilitation is a consequence of altered granule cell presynaptic calcium dynamics. Measurements of synaptic currents were made with whole-cell patch pipettes, which dialyzed the contents of the Purkinje cells with 10 mm EGTA; thus, EGTA-AM loading is unlikely to change significantly calcium buffering within recorded Purkinje cells.
Experiments such as those presented in Figure 6 demonstrate the effects of EGTA-AM on recordings from single cells. Because of the length of these experiments (≥1 hr), we were unable to avoid small increases in access resistance. To minimize the effects of changes in access resistance, we used very low resistance electrodes (0.9–1.5 MΩ) and adjusted the stimulus intensity to keep the amplitudes of unfacilitated EPSCs small. Control experiments on untreated cells recorded for >1 hr (not shown) showed that the time course of facilitation decay was insensitive to access-resistance changes in this range.
To quantitate better the effect of EGTA-AM on facilitation, we also recorded facilitation curves in cells that had been exposed to EGTA-AM, without measuring the control facilitation curves in the same cells. Such experiments had the advantage that more trials could be performed at a constant series resistance, allowing a better estimate of the time course of facilitation. To overcome cell-to-cell variability, we recorded facilitation curves for many cells, with averages shown in Figure 7. The amplitude of facilitation went from 153 ± 11% (n = 15) in control conditions to 159 ± 25% (n = 14), 74 ± 8% (n = 9), and 70 ± 10% (n = 10), respectively, for applications of 1 μm, 20 μm, or 100 μm EGTA-AM. As is seen most clearly in Figure 7B, the time constant of decay of facilitation was also affected, declining from a control value of 203 ± 18 msec to 138 ± 19 msec, 85 ± 11 msec, and 52 ± 4 msec for 1 μm, 20 μm, or 100 μm EGTA-AM, respectively.
Figure 8 compares the time course of [Ca]res transients with the time course of facilitation for control conditions and after altering [Ca]resdynamics. For all of our experimental conditions, [Ca]resdecayed more rapidly than did facilitation (Fig. 8A). As shown in Figure 8B, increasing the decay rate of [Ca]res in turn speeds the decay of facilitation.
Simulations of calcium-activated facilitation
In considering the mechanisms underlying facilitation, we explored the following general reaction scheme, in which calcium binds to a receptor, X, to produce CaX, whose concentration is directly proportional to facilitation, where k+ andk− are the forward and reverse rate constants: Equation 3The relationship between [Ca]res and facilitation could be accounted for by a simple model (Fig. 9,solid traces) that assumes facilitation is the result of a second-order calcium-dependent reaction with one-to-one binding, ak+ of 1.5 × 108m−1 sec−1, and an effectivek− of 25 sec−1, corresponding to a dissociation constant of 167 nm.
The purpose of these simulations was to examine the feasibility of the involvement in facilitation of a high-affinity calcium binding site, and to test the range of parameters that adequately relate the observed calcium transients and facilitation. The strategy was to impose a [Ca]res transient with an amplitude of 200 nm for control conditions (as estimated previously) (Regehr and Atluri, 1995), and a time course determined experimentally, as shown in Figure 5. We then tried to approximate the experimental facilitation curves by varying k+ andk−. k− was estimated from the facilitation observed after exposure to 100 μmEGTA-AM as k− ≈ (τfac)−1, which is ∼25 sec−1. We then tried various values of k+ and found that none of them were consistent with the data. To see why no value ofk+ proved to be satisfactory, consider the two constraints that had to be satisfied simultaneously. First, facilitation is observed to peak <20 msec after the conditioning pulse. For a model of this sort, the rise time of facilitation can be approximated by τrise = (k+[Ca] + k−)−1. Fork− = 25 sec−1 and [Ca] = 200 nm, k+ must be >4 × 108m−1 sec−1. On the other hand, KD must be >150 nm to be consistent with the observed amplitudes of facilitation. BecauseKD = k−/k+,k+ must be <2 × 108m−1 sec−1, a conclusion inconsistent with the constraint imposed by the rapid rise of facilitation. Another complication is that the extremely brief and rather small calcium signals observed for 100 μm EGTA-AM experiments are predicted to produce much less facilitation than the larger, longer-lasting calcium signals of the 20 μmEGTA-AM experiments, yet this is not observed.
Considerations of models of calcium transients within cells suggested a likely explanation for the inadequacies of the simulations described above. In models that take into account the spatial gradients of calcium on rapid time scales, it is apparent that calcium levels can be extremely high near open calcium channels (Simon and Llinas, 1985;Augustine and Neher, 1992; Yamada and Zucker, 1992). Fluorometric measurements such as those used here, however, are not very sensitive to the high local calcium transients near open calcium channels. It seemed to us that the most straightforward explanation of our results is that in addition to the calcium transients we can measure, there is a brief pulse of calcium that produces some facilitation by binding to its receptor X, but which we are unable to detect. The amplitude of the calcium pulse will depend on the location of X relative to the open calcium channels. Very near calcium channels this value could approach 100 μm (Roberts et al., 1990), but at more distant sites this value will be lower. Because we do not know the location of X, we fixed the duration of this pulse at 1 msec and varied its amplitude. This approach provided a good approximation to the observed facilitation for various different presynaptic calcium transients, as shown in Figure 9. These simulations used a 1.5 μmcalcium pulse of 1 msec duration combined with calcium transients with time courses and relative amplitudes determined on the basis of Figure5. This second component increased [Ca]res by 200 nm for control traces.
DISCUSSION
Facilitation at the granule cell→Purkinje cell synapse involves a high-affinity calcium binding site with slow effective kinetics
Our experiments indicate that [Ca]res plays an important role in facilitation. Speeding the decay of [Ca]res levels with EGTA-AM (Figs. 4, 5) accelerated the decay of facilitation (Figs. 6, 7) in a dose-dependent manner. The application of 1 μm EGTA-AM increased the decay rate of both [Ca]res and facilitation without affecting their peak levels. Higher concentrations of EGTA-AM (20 μm) also left peak [Ca]res largely unaffected, whereas they further accelerated the decay of [Ca]res. For these conditions the peak amplitude and duration of facilitation both decreased. We have shown previously that a single stimulus results in a peak increase in [Ca]res of 200–300 nm(Regehr and Atluri, 1995). These results indicate that the processes responsible for facilitation are sensitive to the small elevations in [Ca]res (tens to hundreds of nanomolar) remaining in these terminals in the first second after action-potential invasion, which supports the existence of a high-affinity calcium receptor whose occupancy is related to maintenance of facilitation.
There is also a component of facilitation that persists after [Ca]res has returned to resting levels. This is most apparent after treating presynaptic terminals with high (100 μm) concentrations of EGTA-AM, which speeds the decay of [Ca]res to such an extent that it becomes a brief impulse lasting only a few milliseconds. Facilitation is still present for such conditions, and persists for tens of milliseconds after a conditioning pulse.
Taken together, our results establish that the time course of facilitation is set in part by the decay of [Ca]res and in part by the slow kinetics of either calcium binding or of the calcium-activated process that produces facilitation.
Simulations of facilitation
The simulations of facilitation (Fig. 9) demonstrate that a model for facilitation based on a calcium-driven reaction with second-order kinetics provides an adequate description of facilitation at this synapse. One particularly interesting aspect of these simulations is that it is possible to think in terms of two components of facilitation, one driven by residual calcium and another driven by locally high levels of calcium that our calcium measurement techniques cannot detect (see Materials and Methods). Facilitation curves in conditions in which residual calcium decays extremely rapidly (Fig.8A) (100 μm EGTA-AM experiments) allow us to estimate the intrinsic time constant of facilitation for our experimental conditions to be ∼40 msec. According to reaction scheme (3) above, this intrinsic time constant reflects the reverse rate constant of calcium binding, which is given by (τfac)−1 = 25 sec−1. According to this model, facilitation at this synapse involves a site with rapid forward kinetics (1.5 × 108m−1 sec−1) but with a high affinity for calcium (KD = 167 nm).
In this model, facilitation decay lags behind [Ca]resdecay for two reasons. The first reason is the slow kinetics of calcium dissociation, which determines the time course of facilitation when residual calcium decays very rapidly. The second reason to expect facilitation to lag behind Ca is saturation of its receptor X, which significantly slows the decay of facilitation in control conditions. When X is nearly saturated, an increase in [Ca]resproduces a very small increase in [CaX]. Conversely, for a given decrease in [Ca]res, a much smaller decrease in [CaX] occurs if X is nearly saturated than would occur if X were far from saturation. If the KD of X is comparable in magnitude to the 200–300 nm Ca elevations produced in these terminals by single action potentials, then the saturation of X may contribute to the observed differences between the time course of facilitation and that of [Ca]res transients. In this model, the slow decay of facilitation owing to saturation of X is analogous to the slow decay of the fluorescence transients reported by high-affinity dyes such as fura-2 (Fig. 2C and Materials and Methods) (Regehr and Atluri, 1995).
According to the model presented at the end of the Results section, there is significant calcium binding to the receptor X even at resting calcium levels, because the KD for the calcium binding site is just 167 nm. In the simulations of Figure 9, which assumed a resting calcium of 40 nm, [CaX]/[Xtotal] went from 0.19 to 0.48, an increase of 150%. The amplitude of this increase is about the same as the magnitude of facilitation for these conditions, suggesting that unfacilitated release may also depend on calcium binding to X.
It must be stressed that the working model based on a calcium-driven reaction with simple kinetics (Fig. 9) is still preliminary and that additional experiments will be required to unravel some of the details of the role of calcium in facilitation. One deficiency of the model is that for control conditions elevations of [Ca]res seem to outlive facilitation (Fig. 3A for interpulse intervals >500 msec). This may reflect small contributions of additional processes, such as depression, to use-dependent changes in synaptic strength. In addition, other models, which are slightly more complex, can also describe our data. For example, consider the following general reaction scheme, in which calcium first binds to its receptor X to produce CaX, which then undergoes an activation reaction to produce CaX*, a species directly proportional to facilitation: Equation 4where k+1,k−1, k+2, andk−2 are forward and backward rate constants for steps 1 and 2, respectively. According to this model, when [Ca]res decay has been greatly accelerated by 100 μm EGTA-AM treatment, the duration of facilitation could be dominated by either the kinetics of calcium dissociation,k−1, or the kinetics of CaX* returning to CaX, k−2.
Comparison to other synapses
It is instructive to compare the role of [Ca]res in facilitation at the granule cell to Purkinje cell synapse with that described previously at other synapses. It seems that a high-affinity calcium-binding site distinct from that involved in triggering release is involved in facilitation for the crayfish NMJ and frog tectal synapses (Yamada and Zucker, 1992; Blundon et al., 1993; Delaney and Tank, 1994; Winslow et al., 1994). As we have shown here, this is true also for a synapse in mammalian brain, suggesting that this may be a universal feature of facilitation that extends from invertebrates to mammals. Similarly, delayed release of neurotransmitter, a process long thought to be closely related to facilitation, may involve a high-affinity calcium binding site (Goda and Stevens, 1994).
Important features of facilitation, however, are not universal. The relative importance of calcium dynamics and the kinetics of the calcium-activated process in determining the time course of facilitation depend on the unique properties of a synapse. In one extreme, calcium-activated reaction kinetics dictate the time course of facilitation when they are much slower than calcium dynamics. This is consistent with the conclusion that [Ca]res dynamics does not play a prominent role in determining the time course of facilitation at the crayfish NMJ (Blundon et al., 1993; Delaney and Tank, 1994; Winslow et al., 1994) and at frog tectal synapses (Feller et al., 1996). This also seems to hold in our experiments when 100 μm EGTA-AM has been used to speed the decay of [Ca]res. Conversely, calcium dynamics much slower than the effective reaction kinetics would dictate the time course of facilitation, as suggested by experiments with caged calcium chelators (Kamiya and Zucker, 1994). More generally, both calcium dynamics and effective reaction kinetics combine to determine the time course of facilitation, as is apparent here for the granule cell to Purkinje cell synapse.
Our results differ greatly from the only other direct examination of calcium dynamics and facilitation in mammalian brain of which we are aware (Wu and Saggau, 1994). They reported that at the hippocampal CA3–CA1 synapse, [Ca]res and facilitation have the same decay rate, suggesting that the time course of facilitation closely follows the [Ca]res time course. They also reported that the calcium influx evoked by the second of two closely spaced stimuli is reduced in amplitude, suggesting the paradoxical result that facilitation occurs when there is less calcium entering to trigger release. Interpretation of their experiments may have been confounded, however, by their use of a high-affinity dye to measure [Ca]res. In preliminary studies of [Ca]resat the CA3–CA1 synapse using low- and high-affinity fluorescent calcium indicators, we found that fura-2 begins to saturate and overestimates the decay time of [Ca]res (W. Regehr, unpublished observations) (Fig. 2). We conclude that [Ca]res decays more rapidly than facilitation at the CA3–CA1 synapse, much as shown here for the granule cell to Purkinje cell synapse. It seems likely, therefore, that [Ca]resalso plays a similar role in facilitation at both synapses. Furthermore, we have found that the calcium influx evoked by the second of two closely spaced stimuli remains constant, both in granule cell presynaptic terminals (Regehr and Atluri, 1995) and in CA3 pyramidal cell presynaptic terminals.
Comparison to the role of calcium in augmentation and PTP
The role of [Ca]res in facilitation as described here is reminiscent of that proposed for [Ca]res in augmentation and PTP. These forms of use-dependent synaptic enhancement, which last for tens of seconds, are observed during and after periods of elevated firing that increase [Ca]res. In augmentation and PTP, modest elevations of [Ca]res(tens to hundreds of nanomolar) enhance release by binding to high-affinity sites (Delaney et al., 1989; Swandulla et al., 1991;Delaney and Tank, 1994; Kamiya and Zucker, 1994; Regehr and Mintz, 1994). It seems that facilitation, augmentation, and PTP all involve calcium-activated processes with effective affinities and kinetics suited to their task.
Here we have used a stimulus protocol that emphasizes a rapidly decaying form of enhancement. With longer stimulus trains, additional forms of use-dependent plasticity such as PTP are also present at this synapse, but as yet the details of the action of calcium in such processes have not been described.
Implications for the properties of molecules involved in facilitation
Our findings constrain the properties of the molecule or molecules that govern the time course of facilitation at the granule cell to Purkinje cell synapse. The calcium sensor for facilitation must have a high calcium affinity and must be activated rapidly (i.e., have a fast on-rate). Furthermore, the molecular mechanism responsible for facilitation introduces a delay between calcium and facilitation, either through slow kinetics of calcium binding or of a long-lived product of a calcium-driven reaction. Calmodulin and synaptotagmin III, both of which are present in parallel fibers and sensitive to moderate levels of calcium (Ullrich et al., 1994; Li et al., 1995), are candidate calcium sensors involved in facilitation (Llinas et al., 1991; Rosahl et al., 1995; Schweizer et al., 1995; Sudhof, 1995). A comparison of the kinetics of calcium-triggered processes involving such proteins with the kinetics observed here at the granule cell to Purkinje cell synapse would help to identify the molecules underlying paired-pulse facilitation.
APPENDIX
Simulations of the effect of EGTA on calcium transients
We sped the decay of [Ca]res by introducing EGTA to presynaptic terminals (Figs. 4, 5, 6, 7, 8). Bath-applied EGTA-AM, a membrane-permeant form of the chelator, freely entered cells, where it was deesterified and rendered cell-impermeant, and accumulated within parallel fibers (Tsien, 1981). In this way, achieved levels of intracellular EGTA can greatly exceed the concentration of the bath-applied AM form.
To understand how EGTA speeds the decay of calcium, we simulated calcium dynamics (Fig. 10A) with a single-compartment model using methods similar to those described previously (Tank et al., 1995; Feller et al., 1996). It is thought that in control conditions, when calcium enters a presynaptic terminal, most of the calcium binds rapidly to fast endogenous calcium buffers, and then [Ca]res decays as calcium is extruded from the cytoplasm (Neher and Augustine, 1992; Tank et al., 1995). For our simulations we assumed that the terminal contained a fast endogenous buffer (total [B] = 1 mm, KB = 40 μm, k+ = 5 × 108m−1 sec−1), and calcium was detected with an indicator with properties similar to those of magnesium green (total [Indicator] = 1 μm,KD = 5 μm,k+ = 5 × 108m−1 sec−1). Levels of indicator in our imaging experiments are likely higher than 1 μm. This value was chosen to allow examination of the effects of EGTA without the complication of the additional buffer. The effect of EGTA was also observed when higher concentrations of indicator are present in the terminal. A Michaelis-Menton calcium extrusion pump (KP = 10 μm,Vmax = 4.8 × 10−3m sec−1) was added so that the control ΔF/F transient had at1/2 of 40 msec, as shown in Figure10A. Appropriate total concentrations of EGTA (KEGTA = 200 nm,K+EGTA = 1.5 × 106m−1 sec−1) were used to give calcium decay times similar to those observed in our experiments (compare Figs. 10A and 5B). As we observed experimentally, the addition of concentrations of EGTA sufficient to speed the t1/2 to 2 msec attenuated peak calcium levels.
The key to understanding the actions of EGTA on this time scale lies in realizing that because it is a slow buffer, EGTA takes some time to equilibrate with free calcium. After calcium entry and rapid equilibration with the fast endogenous buffer, calcium begins to bind to EGTA, resulting in a reduction of free calcium levels. The time course of [Ca]res decay, until free calcium equilibrates with EGTA, is determined mostly by the time course of calcium binding to EGTA (compare Fig. 10, A and B). Extending reasoning similar to that of Tank et al. (1995) to the nonequilibrium case, the rapid phase of decay of [Ca]res transients is approximated by an exponential decay with a time constant given by: Equation 5Once calcium has equilibrated with both the endogenous buffer and EGTA, and calcium has been reduced to low levels, the small remaining [Ca]res decays slowly with a time constant approximated by the following equation (Neher and Augustine, 1992; Tank et al., 1995): Equation 6It should be noted that although the control [Ca]restransients we detect experimentally are not perfect monoexponentials, the model described above seems to capture the major features of the speeding of [Ca]res transients by EGTA.
Note added in proof:Bertram et al. (1996) have recently proposed a model of facilitation in which the time course of facilitation is set both by free calcium dynamics and by the kinetics of one or more high-affinity calcium binding sites.
Footnotes
This work was supported by National Institutes of Health Grant R01-NS32405-01, a McKnight Scholars Award, and a Klingenstein Fellowship Award in the Neurosciences to W.R. We thank Bernardo Sabatini, Jeremy Dittman, A. Vyshedskiy, and J.-W. Lin for comments on this manuscript.
Correspondence should be addressed to Dr. Wade G. Regehr, Department of Neurobiology, Harvard Medical School, 220 Longwood Avenue, Boston, MA 02115.