Abstract
The presynaptic Ca2+ signal is a key determinant of transmitter release at chemical synapses. In cortical synaptic terminals, however, little is known about the kinetic properties of the presynaptic Ca2+ channels. To investigate the timing and magnitude of the presynaptic Ca2+ inflow, we performed whole-cell patch-clamp recordings from mossy fiber boutons (MFBs) in rat hippocampus. MFBs showed large high-voltage-activated Ca2+ currents, with a maximal amplitude of ∼100 pA at a membrane potential of 0 mV. Both activation and deactivation were fast, with time constants in the submillisecond range at a temperature of ∼23°C. An MFB action potential (AP) applied as a voltage-clamp command evoked a transient Ca2+ current with an average amplitude of ∼170 pA and a half-duration of 580 μsec. A prepulse to +40 mV had only minimal effects on the AP-evoked Ca2+ current, indicating that presynaptic APs open the voltage-gated Ca2+ channels very effectively. On the basis of the experimental data, we developed a kinetic model with four closed states and one open state, linked by voltage-dependent rate constants. Simulations of the Ca2+ current could reproduce the experimental data, including the large amplitude and rapid time course of the current evoked by MFB APs. Furthermore, the simulations indicate that the shape of the presynaptic AP and the gating kinetics of the Ca2+ channels are tuned to produce a maximal Ca2+ influx during a minimal period of time. The precise timing and high efficacy of Ca2+ channel activation at this cortical glutamatergic synapse may be important for synchronous transmitter release and temporal information processing.
- mossy fiber boutons
- presynaptic Ca2+inflow
- presynaptic Ca2+ channels
- hippocampus
- kinetic model
- glutamatergic synapse
The mossy fiber synapse on CA3 pyramidal neurons is a key synapse that mediates information flow in the trisynaptic circuit of the hippocampus (Jonas et al., 1993; Jung and McNaughton, 1993; Lisman, 1999; Henze et al., 2000). Glutamate release at this synapse is initiated by a complex series of events. These include the invasion of the action potential (AP) into the presynaptic terminal (Geiger and Jonas, 2000), the activation of presynaptic P/Q-, N-, and R-type Ca2+channels (Castillo et al., 1994; Gasparini et al., 2001), and the inflow of Ca2+ ions into the presynaptic terminal. The presynaptic Ca2+ inflow finally triggers the fusion of synaptic vesicles but also contributes to the regulation of synaptic strength (Salin et al., 1996; Tong et al., 1996; Sabatini and Regehr, 1999).
The spatiotemporal profile of the presynaptic Ca2+ transient is of critical importance for the efficacy, timing, and regulation of transmitter release at the mossy fiber synapse. A brief, local transient may trigger fast and highly synchronized glutamate release (Brown and Johnston, 1983; Jonas et al., 1993). In contrast, a more long-lasting, global Ca2+ increase is necessary to evoke the release of peptides, e.g., dynorphin (Verhage et al., 1991; Derrick and Martinez, 1996; Williams and Johnston, 1996). Finally, changes in synaptic strength, such as facilitation, post-tetanic potentiation, and long-term potentiation, are thought to be primarily determined by the residual Ca2+ (Salin et al., 1996). Short-term plasticity is completely suppressed by the slow exogenous Ca2+ buffer EGTA, whereas fast transmitter release is only partly reduced (Salin et al., 1996). This suggests that Ca2+ sensors for release and plasticity sense Ca2+ signals with different spatiotemporal profiles.
The shape and duration of the presynaptic Ca2+ signal is determined by several molecular and structural properties of the presynaptic element, including the duration of the spike, the number of Ca2+ channels per release site, the gating properties of the presynaptic channels, and the concentration of endogenous Ca2+ buffers (Neher, 1998). The calyx of Held is the only mammalian synapse in which these factors have been determined systematically (Forsythe, 1994;Borst and Sakmann, 1996). In cortical glutamatergic synapses, however, these properties are unknown.
To directly examine the presynaptic Ca2+ channels at a cortical glutamatergic synapse, we made patch-clamp recordings from hippocampal mossy fiber boutons (MFBs) in acute brain slices (Geiger and Jonas, 2000; Geiger et al., 2002). We investigated the time course of presynaptic Ca2+ inflow and the gating properties of the presynaptic Ca2+ channels using whole-cell voltage-clamp recordings. Our approach revealed several fundamental functional parameters of synaptic transmission at the mossy fiber synapse, such as the proportion of Ca2+ channels that open during a presynaptic AP and the number of Ca2+ ions entering per spike.
Part of the results have been published previously in abstract form (Bischofberger et al., 2001).
MATERIALS AND METHODS
Slice preparation. Transverse 300-μm-thick slices were cut from the hippocampus of 20- to 26-d-old Wistar rats with a custom-built (Geiger et al., 2002) or a commercially availableVibratome (DTK-1000, Dosaka, Japan). The animals were killed by decapitation, in accordance with national and institutional guidelines. Slices were kept at 35°C for 30 min after slicing and then at room temperature. For the dissection and storage of slices, we used a saline containing (in mm): 64 NaCl, 25 NaHCO3, 10 glucose, 120 sucrose, 2.5 KCl, 1.25 NaH2PO4, 0.5 CaCl2, and 7 MgCl2. For the experiments, the slices were transferred into a recording chamber and superfused with a physiological extracellular solution containing (in mm): 125 NaCl, 25 NaHCO3, 25 glucose, 2.5 KCl, 1.25 NaH2PO4, 2 CaCl2, and 1 MgCl2(equilibrated with 95% O2/5% CO2 gas mixture).
Electrophysiology. Mossy fiber boutons in stratum lucidum of the hippocampal CA3 region were visually identified by their size and location (3–4 μm apparent diameter) (see Fig.1A) (Geiger and Jonas, 2000) using infrared differential interference contrast videomicroscopy (Stuart et al., 1993). Identification was confirmed by the small capacitance (∼1 pF whole-cell capacitance readout of the amplifier) and the high input resistance (>1 GΩ). A subset of boutons was filled with biocytin for 30 min in the whole-cell recording configuration. Slices were fixed in 4% paraformaldehyde, stained with Fluorescein Avidin D (Vector Laboratories, Burlingame, CA), and embedded with ProLong Antifade (Molecular Probes, Eugene, OR). Images were taken using a confocal laser scanning microscope (LSM 510, Zeiss, Göttingen, Germany) (see Fig. 1A).
Patch pipettes were pulled from borosilicate glass tubing (2.0 mm outer diameter, 0.7 mm wall thickness; Hilgenberg, Malsfeld, Germany). A modified Axopatch 200A amplifier (Axon Instruments, Foster City, CA) was used for current-clamp (I-Clamp fast) and voltage-clamp recording. The amplifier included a bridge-balance circuit for compensation of series resistance in the current-clamp mode, similar to that available in the Axopatch 200B. Current and voltage signals were filtered at 5 and 10 kHz, respectively, with a four-pole low-pass Bessel filter and stored online using a CED 1401plus interface (Cambridge Electronic Design, Cambridge, UK) connected to a personal computer. The sampling frequency was two to four times the filter frequency. Pulses were generated using commercial (CED) or self-made (P. Jonas) programs; the latter allowed us to apply previously recorded APs as voltage-clamp commands.
For whole-cell current-clamp recordings, we used 7–10 MΩ patch pipettes filled with internal solution containing (in mm): 150 KMeSO4, 0.2 EGTA, 2 MgCl2, 2 Na2ATP, 0.3 NaGTP, and 10 HEPES (the pH was adjusted to 7.3 with KOH). Bridge balance was used to compensate the series resistance. In addition, capacitance compensation was used to decrease the charging time of the pipette to <50 μsec. Resting potentials were between −65 and −80 mV. Membrane potentials were set to −80 mV by applying a small hyperpolarizing holding current if necessary. Stimulation of the mossy fiber tract was made with a patch pipette located in the stratum lucidum in subregion CA3c. A stimulus isolation unit (List, Darmstadt, Germany) was used to generate pulses with a duration of 200 μsec and amplitude of −5 to −10 V.
For whole-cell voltage-clamp recordings, we used 5–7 MΩ pipettes filled with a solution containing (in mm): 135 CsCl, 10 EGTA, 4 MgCl2, 4 Na2ATP, 0.3 NaGTP, 10 Na2-phosphocreatine, and 10 HEPES (the pH was adjusted to 7.3 with CsOH). The bath solution contained (in mm): 105 NaCl, 25 NaHCO3, 25 glucose, 2.5 KCl, 1.25 NaH2PO4, 2 CaCl2, 1 MgCl2, 1 μm tetrodotoxin, 20 tetraethylammonium chloride, and 5 4-aminopyridine to block voltage-gated Na+and K+ channels, respectively. Leak and capacitive currents were subtracted using two P/−4 sequences within each protocol. Current traces shown in the figures represent the averages of two to five sweeps. Series resistance compensation was enabled in all experiments (compensation 80–90%; lag 20–30 μsec).
The series resistance was typically between 20 and 50 MΩ. Assuming a membrane capacitance of ∼1 pF (Geiger and Jonas, 2000), this would correspond to a voltage-clamp time constant τ =RsCmof 20–50 μsec. If the series resistance exceeded a value of 50 MΩ, the recordings were aborted. In some recordings, unclamped tail currents with very slow time course were apparent (these currents may have been generated in neighboring boutons). These recordings were discarded. All chemicals were obtained from Sigma (St. Louis, MO), Merck (Darmstadt, Germany), Riedel-de Haën, or Gerbu. Recordings were made at 23 ± 1°C.
Data analysis. The activation time course of the Ca2+ current was fitted with a monoexponential function with delayed onset: Equation 1where I0 is the amplitude, τ is the activation time constant, and δ is the delay of Ca2+ channel activation. To account for the small delay introduced by the lag of series resistance compensation (20–30 μsec) and by electronic filtering (mostly at 5 kHz), the current trace was shifted in time by 100 μsec; with this shift, the end of the rectangular current pulse corresponded to the time when the tail current rose to ∼50% of its peak amplitude. The time course of deactivation was fitted with a monoexponential function starting from 50–100 μsec after the peak of the tail currents to minimize contamination by residual capacitive artifacts. The number of Ca2+ ions was calculated from the total charge q as n = qNA/z F with the Avogadro number NA = 6.02 1023/mol, the Faraday constantF = 96,485 C/mol, and z = 2. Time course fitting was performed using programs written in Pascal (P. Jonas) based on a modified Gauss–Newton algorithm.
The steady-state activation curve was fitted with a Boltzmann function of the form f(V) = 1/(1 + exp[(Vhalf −V)/Vslope]) withVhalf = −3.9 mV andVslope = 7.1 mV as best-fit parameters. The values for Vhalf andVslope were then used to fit theI–V relationship with a modified Goldman–Hodgkin–Katz equation of the form: Equation 2with the voltage V, an amplitude factor P, a steepness factor C, and a parameter Ddetermining current rectification and reversal potential (Sala, 1991). The best-fit parameters were P = −3.003 pA/mV,C = 80.36 mV, and D = 0.3933. The voltage dependence of activation time constants was fitted with the function: Equation 3where k1 arek2 are coefficients andV1 andV2 are steepness factors. The voltage dependence of current amplitudes and time constants was fitted using Prism 3.0 (GraphPad, San Diego, CA).
Amplitude and duration at half-maximal amplitude of APs were measured from baseline potential. Membrane potentials are given without correction for liquid junction potentials. The values indicate mean ± SEM. The error bars also represent SEM. The significance of differences was assessed by a two-tailed Wilcoxon–Mann–Whitney test (Prism 3.0) at the significance level (P) indicated.
Modeling. The simulations were performed with Mathematica 4.0 (Wolfram Research, Champaign, IL). A set of linear differential equations was defined for a kinetic model with four closed and one open state with voltage-dependent microscopic rate constants for the forward [αi(V)] and backward [βi(V)] transitions. The differential equations were defined using the Q-matrix approach ofColquhoun and Hawkes (1977) and numerically solved using Mathematica's NDSolve routine. The analysis of simulated data was identical to that of experimental data; the simulated activation and deactivation time course of the open state were fitted with a monoexponential function with delayed onset and an exponential decay, respectively. The values for activation τ, delay, deactivation τ, and steady-state activation were compared with experimental data by calculating the total sum of the squares of the differences. This sum was minimized by using the FindMinimum routine of Mathematica.
The microscopic rate constants were exponentially dependent on voltage as described by the equations: Equation 4with the best-fit parameters α1o = 4.04 msec−1, β1o = 2.88 msec−1,V1 = 49.14 mV, α2o = 6.70 msec−1, β2o = 6.30 msec−1,V2 = 42.08 mV, α3o = 4.39 msec−1, β3o = 8.16 msec−1,V3 = 55.31 mV, α4o = 17.33 msec−1, β4o = 1.84 msec−1, andV4 = 26.55 mV.
RESULTS
Whole-cell recordings were made from MFBs in acute rat hippocampal slices (Fig. 1). MFBs were identified by their small size, their small capacitance (∼1 pF), and their high input resistance (>1 GΩ) (Geiger and Jonas, 2000). A subset of the boutons was filled with biocytin and stained with Fluorescein Avidin D. Figure 1A shows a confocal image of a biocytin-filled MFB, which illustrates the staining of the parent axon and an adjacent MFB. The image also shows filopodial extensions originating from the boutons, which are highly characteristic structural properties of mossy fiber terminals (Acsády et al., 1998).
The presynaptic AP is the primary event that activates presynaptic Ca2+ channels and thereby triggers transmitter release. Similar to other presynaptic elements (Borst and Sakmann, 1996; Sabatini and Regehr, 1996; Bischofberger and Jonas, 1997), the MFBs were shown to have APs with very rapid time course (Geiger and Jonas, 2000). Figure 1B shows a current-clamp recording from an MFB. A small current injection of 20–40 pA evoked an AP with an amplitude of 122.9 ± 2.4 mV and a half-duration of 915 ± 44 μsec (range, 540–1140 μsec;n = 13). Similarly, APs could be evoked by extracellular stimulation of the mossy fiber tract (Fig.1C). The amplitude (122.8 ± 2.7 mV) and half-duration (852 ± 62 μsec; n = 6) of the AP evoked by remote stimulation were not significantly different from those evoked by direct current injection (p > 0.4). In contrast, the half-duration of APs recorded from granule cell somata was much longer, as shown in Figure 1D (2.06 ± 0.09 msec; n = 9; p < 0.01). Thus both locally evoked and orthodromically propagated APs in MFBs showed a fast time course that was significantly faster than APs in granule cell somata.
MFB Ca2+ channels have steep voltage dependence and rapid activation kinetics
To determine the proportion of Ca2+channels opened by the fast presynaptic AP waveform, we analyzed the gating properties of the presynaptic Ca2+channels in MFBs (Figs. 2-4). Figure2A shows whole-cell voltage-clamp recordings of pharmacologically isolated Ca2+ currents evoked by 20 msec voltage steps from −70 mV to +60 mV, and Figure2B illustrates the corresponding current–voltage (I–V) relationship. The Ca2+ currents were activated above −40 mV and reached a maximal current amplitude of 98 ± 11 pA at 0 mV (n = 16). At more positive potentials the inward currents decreased, resulting from a reduction in the driving force (Augustine et al., 1985; Borst and Sakmann, 1998).
To obtain the steady-state activation curve, we plotted the integral of the tail currents against the voltage of the preceding pulse (Fig.2C). Because the tail currents were measured at the same potential, the integral gives a precise measure of the number of channels open at the end of the preceding pulse. The steady-state activation data were adequately fitted by a Boltzmann function with a slope factor of 7.1 mV and a midpoint potential of −3.9 mV. We then fitted the bell-shaped I–V relation in Figure2B with the product of a Boltzmann function representing channel activation and a modified Goldman–Hodgkin–Katz equation describing the driving force for ion flow (Sala, 1991). The parameters of the Boltzmann function were constrained to the previously determined values. From the fitted curve, the extrapolated reversal potential was determined to be 75.0 mV (Fig. 2B). In conclusion, presynaptic Ca2+ channels in MFBs show a high activation threshold and a steep voltage dependence.
We next examined the activation and deactivation kinetics of presynaptic Ca2+ channels. The time course of activation could be well fitted by a monoexponential function with delayed onset relative to the onset of the pulse (Fig.3A,B). The activation time constant τ and the delay are plotted in Figure3C as a function of voltage. The activation τ was fast and voltage dependent, ranging from τ = 842 ± 137 μsec at −10 mV to 201 ± 35 μsec at 50 mV (n = 6). The voltage dependence could be fitted by Equation 3 with coefficientsk1 = 1.12 msec−1 andk2 = 0.14 msec−1 as well as slope factors V1 = 31.5 mV andV2 = 8.6 mV. In contrast, the delay was independent of voltage, with a mean value of ∼200 μsec. Linear regression revealed that the slope of the fitted line was not significantly different from 0 (p > 0.4).
Because the amount of Ca2+ inflow depends on how long the channels remain open during the presynaptic AP, we further analyzed the deactivation kinetics after a voltage step from 0 mV to negative potentials (−10 mV to −60 mV) (Fig.4A). As shown in Figure4B, the deactivation could be adequately fitted by a monoexponential function. The deactivation τ was fast, ranging from 813 ± 90 μsec at −10 mV to 148 ± 11 μsec at −60 mV (n = 6). The voltage dependence could be described with an exponential function with an e-fold change per 15.0 mV (Fig. 4C). Thus both activation and deactivation time constants were fast and highly voltage dependent.
Presynaptic Ca2+ channels are effectively opened by MFB APs
To study the presynaptic Ca2+ inflow during natural stimuli, we used previously recorded AP waveforms as voltage-clamp commands (Scroggs and Fox, 1992) (Fig.5). A propagated AP evoked a large and brief Ca2+ inflow starting at the peak of the AP (Fig. 5A). During the repolarization phase the current rapidly increased, reaching a peak amplitude of 187 ± 36 pA (n = 9). This current was substantially larger (177 ± 17%) than the maximum peak current amplitude during rectangular voltage pulses to 0 mV in the same MFBs (102 ± 16 pA;p < 0.01). The half-duration of the AP-evoked Ca2+ current was very short (540 ± 12 μsec; n = 9), most probably because of the fast deactivation of the channels at negative potentials. Similar results were obtained with a locally evoked AP waveform (Fig. 5C,dashed line). The relative peak amplitude and half-duration were 175 ± 10% and 590 ± 19 μsec (n = 18), very similar to the respective values for the propagated AP. If data with both waveforms were pooled, a peak amplitude of 171 ± 20 pA (corresponding to a relative amplitude of 173 ± 9%), a half-duration of 579 ± 15 μsec, and a total charge of 119 ± 13 fC were obtained (n = 24). This corresponds to an inflow of ∼370,000 Ca2+ ions per single AP.
To estimate the relative open probability (Popen) of the Ca2+ channels during the AP, we added a prepulse to +40 mV to the AP waveform (Fig.5B,D) (Borst and Sakmann, 1998). According to the steady-state activation curve (Fig. 2C), this pulse should activate the Ca2+channels with a maximal Popen. The peak amplitude of the evoked Ca2+currents, however, was only slightly larger as compared with control (124 ± 3% of control; p < 0.01;n = 9). This indicates that the presynaptic Ca2+ channels are opened very effectively by the AP waveform and that the peakPopen with a natural stimulus reaches ∼80% of the maximal value.
To further elucidate the relation between AP waveform and Ca2+ channel activation, we compared the Ca2+ current evoked by a fast MFB AP with that evoked by a slow somatic AP waveform (Fig.5C,D). The somatic AP evoked a Ca2+ current with a half-duration of 1232 ± 51 μsec, significantly larger than that evoked by the brief AP (p < 0.01). However, the peak amplitude was reduced (84.6 ± 1.1% of control; p< 0.01; n = 5). This decrease in peak current amplitude appears to be fundamentally different from the increase in current during AP broadening reported previously for both somatic (Scroggs and Fox, 1992; Wheeler et al., 1996) and presynaptic Ca2+ channels (Augustine, 1990; Sabatini and Regehr, 1997). In conclusion, the data suggest that MFB Ca2+ channels are effectively opened by the presynaptic APs. To understand the complex dynamic process of activation and deactivation during AP waveforms, we developed a computational model of Ca2+ channel gating.
Ca2+ channel gating in MFBs can be described by a five-state kinetic model
The pore-forming α1 subunit of voltage-gated Ca2+ channels consists of four domains, each having a positively charged S4 segment, which is believed to be involved in the voltage-dependent gating of the channels (Hille, 2001). Therefore, a model with four voltage-dependent transitions between four closed states and one open state was used (Fig. 6A). The microscopic rate constants αi and βi were assumed to be exponentially dependent on voltage (Eq. 4) (Fig. 6B). For a given parameter set, the occupancy of the open state during rectangular voltage pulses was simulated, and the resulting traces were fitted with exponential functions, similar to the experimental data. The parameters for the rate constants αi and βi were changed until the best approximation of the experimental values for activation τ, delay, deactivation τ, and steady-state activation were achieved. Figure 6, C andD, shows a comparison of the experimental data with the predictions of the final model (continuous lines). In the final set of rate constants, the last transition shows the steepest voltage dependence (Vi = 49, 42, 55, and 27 mV), suggesting cooperativity in channel gating (Fig.6B) (see also Materials and Methods).
We then simulated the Ca2+ inflow during APs, using the established model of Ca2+channel gating and the driving force from theI–V relation fitted with Equation 2 (Fig.7). If the brief MFB AP was applied as a voltage command, the simulated Ca2+current had a peak amplitude of 161 pA and a half-duration of 596 μsec. The peak amplitude corresponded to 176% of the current amplitude during a rectangular voltage pulse to 0 mV, very similar to the experimental data (relative current amplitude 173%; half-duration 579 μsec) (Fig. 5). As shown in Figure 7A, the channels started to open during the rising phase of the AP and reached an open probability of 51% at the peak of the AP. The maximal open probability occurred during the early repolarization phase at a potential of 15 mV (Popen = 90%). The peak Ca2+ inward current occurred slightly later, at a potential of −20 mV (Popen = 68%). Subsequently, the channels closed very rapidly and deactivated before complete repolarization to the resting potential. Consistent with the effective opening of the Ca2+ channels during the MFB AP, the peak current amplitude was only slightly increased by a 5 msec prepulse to 40 mV (Fig. 7B).
To examine the mechanisms underlying the reduction in peak amplitude of the Ca2+ current during the slower somatic AP (Fig. 5C), we used simulations to compare Ca2+ channel activation by fast and slow APs. Similar to the experimental data, the simulated Ca2+ current during a broad somatic AP had a prolonged time course and a reduced amplitude (Fig. 7C). The simulation shows that the reduction in peak amplitude is caused by a reduction in driving force at this point (−10 vs−20 mV). By contrast, the open probability of the channels at the peak of the Ca2+ current was very similar. Figure7D shows a comparison of the relative peak amplitude and half-duration of simulated Ca2+ currents evoked by different waveforms. The simulated values for the MFB AP were in close agreement with experimentally determined values (Fig. 5). The results of experiments and modeling converge toward the conclusion that a brief AP waveform leads to efficient and rapid activation of the Ca2+ channels, resulting in a high peak open probability.
Open probability of presynaptic Ca2+ channels is highly sensitive to changes in AP waveform
The shape of the presynaptic AP in MFBs is not constant but regulated in an activity-dependent manner (Geiger and Jonas, 2000). Thus we systematically examined the effects of modified AP waveforms, scaled in amplitude and shape in comparison with the original MFB AP (Fig. 8). Figure 8Ashows the effect of amplitude scaling on the evoked Ca2+ inflow, and Figure8B illustrates plots of peak amplitude and charge as a function of AP amplitude. The simulations show that the control AP waveform (vertical lines) is a very effective stimulus, producing a presynaptic Ca2+ signal with near-maximal peak amplitude and total charge (89 and 98% of the maximal value, respectively). By contrast, a reduction of the AP amplitude by Na+-channel inactivation during high-frequency activity or steady-state depolarization, for example by presynaptic kainate receptors, will substantially reduce the presynaptic Ca2+ inflow (Kamiya and Ozawa, 2000; Schmitz et al., 2000).
Similarly, we modified the half-duration of the AP waveform by slowing the time course of repolarization in comparison with the MFB AP (Fig.8C,D). The peak amplitude was exponentially dependent on the half-duration (Fig. 8D, top graph). In contrast, the total charge could be fitted by linear regression (Fig. 8D, bottom graph), yielding a slope of 73 fC/msec. Thus the model predicts a total Ca2+ inflow of ∼77–121 fC for the experimental half-durations ranging from 540 to 1140 μsec. Furthermore, it explains the previous observation that dynamic AP broadening leads to a decrease in peak amplitude but an increase in total charge of presynaptic Ca2+ inflow (Geiger and Jonas, 2000).
In conclusion, the simulations show that the natural AP is a highly effective stimulus for the activation of the presynaptic Ca2+ channels. Additionally, the Ca2+ signal can be regulated via the half-duration of the presynaptic AP waveform.
Ca2+ channel gating is optimal for large Ca2+ inflow
Similarly, we examined the effects of alterations in channel gating on the amplitude and charge of Ca2+inflow (Fig. 9). Figure 9Ashows the effects of slower gating, generated by multiplication of all rates with constant factors in the range 1–0.1, and Figure9B shows the effects of acceleration of gating, using factors in the range 1–10. In all simulations, the natural MFB AP waveform was used as a stimulus. Plots of peak amplitude and total charge against the scaling factor indicate that gating kinetics is a critical determinant of the total Ca2+inflow (Fig. 9C,D). Surprisingly, both slowing down and speeding up the rate constants reduced the Ca2+ peak current and the integral. The peak amplitude and total charge produced by the native channels (scaling factor of 1) were 99 and 97% of the maximal values, which were reached with factors slightly below 1. Furthermore, both slowing down and speeding up of the rates lead to prolonged time course of Ca2+ inflow. In particular, the acceleration of the rates by a factor of >2 leads to a biphasic Ca2+ current, increasing the total time interval during which Ca2+ inflow occurs (Fig. 9B). In conclusion, these simulations unequivocally show that the gating of presynaptic Ca2+channels in MFBs is optimal for producing a transient Ca2+ inflow with maximal amplitude and minimal duration. Whether this conclusion also holds true at physiological temperatures remains to be established. However, because the Q10 values for AP duration (∼2.2) (Geiger and Jonas, 2000; this paper) and Ca2+channel activation kinetics are similar (∼2.5) (Swandulla and Armstrong, 1988), this appears to be very likely.
DISCUSSION
The present results provide the first description of Ca2+ channel gating in a cortical presynaptic terminal. The major finding was that the fast activation kinetics of the channels leads to a high open probability and a remarkably large Ca2+ inflow of ∼370,000 Ca2+ ions during a single AP. The gating of presynaptic Ca2+ channels appears to be precisely tuned to achieve a maximal Ca2+influx during a minimal period of time.
Rapid gating of presynaptic Ca2+ channels
The gating of presynaptic Ca2+channels could be adequately described by a serial gating model with four closed states and one open state. In contrast, neither alternative models with two or three closed states nor Hodgkin–Huxley models with an integer number of independent gating particles provided a sufficient description of the data (J. Bischofberger, unpublished data). In the models with two or three closed states, the delay of activation was smaller than that observed experimentally. In Hodgkin–Huxley models, deactivation is faster than activation at the same potential, with a ratio that approaches the number of gating particles (Kay and Wong, 1987; Zidanic and Fuchs, 1995; Mennerick and Matthews, 1998). Thus the functional properties of presynaptic Ca2+channels in MFBs are inconsistent with the assumption of independently moving gating particles.
Our model suggests that the voltage dependence differs among the transition rates. This has to be interpreted as a different movement of the S4 segments, the putative voltage sensors, because the number of positive charges per S4 segment is comparable among domains (approximately five) (Soong et al., 1993). For the first three transitions, the steepness factors range from 42 to 55 mV, which would correspond to an S4 segment movement by 9–12% along the total electric field. In contrast, for the final transition leading to the open state, the steepness factor is 27 mV, which would correspond to an S4 segment movement by 19%. Interdomain differences in voltage-dependent movement of S4 segments were experimentally demonstrated for voltage-dependent Na+channels (Cha et al., 1999). An analogous biophysical approach will reveal whether similar interdomain differences exist in voltage-gated Ca2+ channels, as our results would suggest.
Efficiency and timing of Ca2+ channel opening during single APs
The activation kinetics of the presynaptic Ca2+ channels in MFBs appear to be much faster than the activation of Ca2+channels in invertebrate presynaptic terminals (Llinás et al., 1981; Wright et al., 1996). For the squid giant synapse it was suggested that the relative open probability of Ca2+ channels during a single AP was ∼10% (Pumplin et al., 1981; Augustine, 1990), markedly less than the value estimated for the MFB (∼90%). The gating is also faster than that described for somatodendritic high-voltage-activated Ca2+ channels in CA1 pyramidal neurons (Kay and Wong, 1987). This may be correlated to the slower time course of somatodendritic voltage signals, e.g., synaptic potentials or backpropagating APs (Magee and Johnston, 1995; Normann et al., 2000).
The activation kinetics of the Ca2+channels in MFBs, however, is similar to that reported for presynaptic Ca2+ channels in giant vertebrate synaptic terminals (Zidanic and Fuchs, 1995; Borst and Sakmann, 1998; Mennerick and Matthews, 1998). In the calyx of Held, single APs evoke a Ca2+ current with a relative peak open probability of ∼70% and an average half-duration of ∼400 μsec (24°C) (Borst and Sakmann, 1998), comparable with the values reported here. Thus, although the two types of synapses differ in morphological properties, subunit composition of postsynaptic receptors, and functional role in the neuronal circuit, the efficacy and timing of presynaptic Ca2+ inflow appear to be similar (Von Gersdorff and Borst, 2002).
In the MFB, the Ca2+ channels open very early during the rising phase of the AP with a relativePopen of ∼50% at the peak of the AP; however, no substantial current can be detected at this point, because of the small driving force. The Ca2+ inflow is restricted to the repolarization phase of the AP (Figs. 5A, 7A). In contrast, for parallel fiber synapses in the cerebellum, it was reported that rapid Ca2+ channel gating leads to substantial Ca2+ current during the upstroke and peak of the presynaptic AP (Sabatini and Regehr, 1996,1997). This situation appears to be very different from that in the MFB because the experimentally determined I–Vrelationship for the Ca2+ currents in the MFB is obviously inconsistent with a large driving force at the peak of the AP (at ∼45 mV) (Fig. 2B). If these discrepancies are not caused by technical difficulties in deriving the voltage and Ca2+ current time course from fluorescence signals (Sabatini and Regehr, 1996, 1997), they may suggest a synapse-specific difference in presynaptic ion flow.
How many Ca2+ ions are necessary for transmitter release?
The present results provide quantitative estimates of fundamental parameters of the early presynaptic events that are relevant for both efficacy and timing of synaptic transmission. We found that the average peak amplitude of the presynaptic Ca2+current evoked by a natural AP is 170 pA. With a single-channel current of 0.2 pA at physiological external Ca2+concentrations (Gollasch et al., 1992), this corresponds to 850 Ca2+ channels open at the peak in a presynaptic terminal. The total Ca2+inflow evoked by a single AP was 120 fC, corresponding to ∼370,000 Ca2+ ions per bouton.
MFBs have multiple release sites (up to 37 in an incompletely reconstructed large MFB) (Chicurel and Harris, 1992). Thus the measured presynaptic Ca2+ current would approximately correspond to 23 Ca2+channels open at the peak at a single release site, producing an inflow of ∼10,000 Ca2+ ions per site. Although these are only rough estimates, the results are consistent with previous suggestions at boutons of neocortical pyramidal cells (∼40 channels per bouton, ∼5500 Ca2+ ions) (Koester and Sakmann, 2000) and the calyx of Held (∼20 channels per site, ∼13000 Ca2+ ions) (Borst and Sakmann, 1996; Sätzler et al., 2002). Furthermore, they are consistent with the cooperativity between different types of Ca2+ channels in triggering transmitter release at the mossy fiber–CA3 synapse (Castillo et al., 1994) and other cortical glutamatergic synapses (Luebke et al., 1993; Wheeler et al., 1994; Dunlap et al., 1995). Although we cannot exclude a subpopulation of release sites endowed with a single Ca2+ channel (Stanley, 1997), our results suggest that most of the sites are controlled by multiple channels. This may increase the reliability of excitation–secretion coupling in central glutamatergic synapses.
Functional significance for synaptic transmission at MFB–CA3 pyramidal neuron synapse
The quantitative properties of the presynaptic Ca2+ inflow are also critical for the timing of transmitter release. A large Ca2+ signal may help to minimize synaptic delay. Transmitter release requires the binding of multiple Ca2+ ions to the Ca2+ sensor (Augustine, 2001). As shown for the calyx of Held, the sequential binding of multiple Ca2+ ions introduces a significant delay between the Ca2+ inflow and the transmitter release (Schneggenburger and Neher, 2000). This delay was strongly dependent on the Ca2+concentration. Thus, a large Ca2+ signal will provide a rapid onset of neurotransmitter release. Additionally, the fast deactivation of the presynaptic Ca2+ channels may be important for the termination of exocytosis. The rapid Ca2+signaling will be a critical determinant of the precise timing of transmitter release at the MFB to CA3–pyramidal cell synapse, where the half-width of the first latency distribution is ∼500 μsec (Jonas et al., 1993). Similar rules may apply for synapses formed by the mossy fiber system onto other types of target cells, e.g., at the granule cell to basket cell synapse, formed by mossy fiber collaterals on inhibitory interneurons (Geiger et al., 1997).
The precise timing during single APs may be relevant for the temporal encoding of information in the hippocampus (Lisman, 1999). It was shown that granule cells and pyramidal cells fire APs in a specific temporal relationship to the theta rhythm when the animal passes through the place field of the cell (Jung and McNaughton, 1993; O'Keefe and Recce, 1993). When the animal is in the center of the place field of a given neuron, this neuron fires APs earlier in the theta cycle than neurons with adjacent place fields (Skaggs and McNaughton, 1996; Skaggs et al., 1996). For temporal coding with high fidelity, neurons with the same place field should fire exactly at the same time (Lisman, 1999). Because the mossy fiber tract is the major excitatory pathway that connects the associative network in the dentate gyrus/hilus with that in the CA3 region, the precise timing of presynaptic Ca2+ inflow and transmitter release might be important for the correct temporal sequence of place cell firing and, more generally, for the temporal coding of episodic memory.
Footnotes
J.B. was supported by grants from the Deutsche Forschungsgemeinschaft (Bi 642/1-2 and SFB 505/C9). We thank Dr. U. Kraushaar, Dr. S. Hefft, and C. Schmidt-Hieber for critically reading this manuscript, F. Heyde for secretarial help, and A. Blomenkamp and K. Winterhalter for technical assistance.
Correspondence should be addressed to Dr. Josef Bischofberger, Physiologisches Institut, Universität Freiburg, Hermann-Herder-Strasse 7, D-79104 Freiburg, Germany. E-mail:Josef.Bischofberger{at}physiologie.uni-freiburg.de.