|Home | About | Journals | Submit | Contact Us | Français|
During short bursts of neuronal activity, changes in the efficacy of neurotransmitter release are governed primarily by two counteracting processes: (1) Ca2+-dependent elevations of vesicle release probability and (2) depletion of synaptic vesicles. The dynamic interplay of both processes contributes to the expression of activity-dependent synaptic plasticity. Here, we exploited various facets of short-term plasticity at the Drosophila neuromuscular junction to dissect these two processes. This enabled us to rigorously analyze different models of synaptic vesicle pools in terms of their size and mobilization properties. Independent of the specific model, we estimate ~300 readily releasable vesicles with an average release probability of ~50% in 1 mM extracellular calcium (~5% in 0.4 mM extracellular calcium) under resting conditions. The models also helped interpreting the altered short-term plasticity of the previously reported mutant of the active zone component Bruchpilot (BRP). Finally, our results were independently confirmed through fluctuation analysis. Our data reveal that the altered short-term plasticity observed in BRP mutants cannot be accounted for by delocalized Ca2+ channels alone and thus suggest an additional role of BRP in short-term plasticity.
Synapses are specialized intercellular contact sites, which serve as the communication link between neurons and their partner cells. At chemical synapses, calcium-ion (Ca2+) influx triggers the fusion of transmitter-laden vesicles with the presynaptic membrane at a specific subcellular region termed the active zone. Transmitter substances released by this process then diffuse across the synaptic cleft and are sensed by the postsynaptic reception apparatus to convey signal transduction. The discovery that synaptic transmission is not only highly dynamic but that its strength also crucially depends on the recent history of presynaptic activity, has attracted decades of considerable scientific interest, and has identified synapses as key regulators of complex brain processes. Activity-dependent changes in synaptic function can last from milliseconds to hours or days, and this plasticity enables synapses to filter, modify, or integrate information. In appropriate neural networks, synaptic plasticity gives rise to processes such as sensory adaptation, rhythm generation, or learning and memory.
During short-term plasticity, the postsynaptic response to brief bursts of presynaptic activity may increase (facilitation) or decrease (depression) over time. Besides postsynaptic changes that contribute to these phenomena, the main presynaptic influences appear to be exerted by two counteracting processes: hand-in-hand with the proposal of the quantal nature of transmitter release, went the suggestion that depression was brought about by a depletion of transmitter “parcels,” whereas facilitation was caused by an elevated probability of their release (del Castillo and Katz, 1954a, 1954b). Derived from this, the postsynaptic response can be expressed as the product of quantal size (q, usually taken as the postsynaptic sensitivity), the number of synaptic vesicles available for release (N, or alternatively: the number of available release sites multiplied by their probability of occupancy), and the probability of transmitter release from a vesicle (p or pr).
Regarding p, the calcium hypothesis proposed that the release of neurotransmitter is triggered by elevations of the Ca2+ concentration in the presynaptic terminal (Katz and Miledi, 1965). This in turn, led to the suggestion that facilitation of release during closely spaced stimuli was caused by the intracellular build-up of residual Ca2+ (Katz and Miledi, 1968).
To address N and to describe the depletion of the store of transmitter during synaptic depression, the concept of functionally distinct vesicle pools, with respect to their release and mobilization properties, was adopted (Liley and North, 1953; Elmqvist and Quastel, 1965). In particular, the idea that a “readily releasable pool” (RRP, Birks and Macintosh, 1961) of synaptic vesicles is supplied by a larger storage pool has been used to explain short-term plasticity phenomena, such as paired-pulse facilitation or depression, in terms of either a change in the size of the RRP or a change in the release probability of the vesicles comprising this pool.
The various expression patterns of short-term plasticity observed at different synapses, in various Ca2+ concentrations, and during different frequencies of activity, appear to mainly stem from distinct changes in the influence of p and N on release. And in this respect, we are essentially still faced with the difficulty first formulated by del Castillo and Katz (1954a). To accurately dissect changes in p and N, we need to know “[…] more about the molecular nature of the reaction whose probability we are considering.”
In the present study, we addressed this challenge, by making use of recent advances in our understanding of molecular aspects of presynaptic function. The work was carried out at the larval neuromuscular junction (NMJ) of Drosophila. This synaptic system is highly accessible to genetic manipulations, which can be employed to give insight into molecular mechanisms of synaptic communication. Our approach was based on applying different stimulation protocols in several extracellular Ca2+ concentrations ([Ca2+]e) to sample a broad range of synaptic plasticity. We then utilized several previously reported, independent models of activity-dependent facilitation (Tsodyks and Markram, 1997; Trommershäuser et al., 2003). By focusing mainly on a mechanistic description of the relationship between vesicle release probability and Ca2+ dynamics (Trommershäuser et al., 2003), information was gained regarding p. This was then used to examine the behavior of various vesicle pool concepts during plasticity and to resolve N. Based on their ability to faithfully reproduce the experimental data, the pool concepts were evaluated and their implications for active zone function at the Drosophila NMJ were interpreted. Additionally, the resting parameters determined in this manner were independently supported by nonstationary fluctuation analysis.
Finally, the models were applied to analyze the physiology of synapses lacking the active zone component Bruchpilot (BRP, Kittel et al., 2006b; Wagh et al., 2006). Our results indicate that BRP mutant active zones are not only affected by a reduced release probability but also by impaired trafficking of synaptic vesicles.
A broad range of short-term plasticity (STP) was evoked through prototypical stimulation protocols and sampled via the peak amplitude of excitatory postsynaptic currents (EPSCs) in voltage clamped muscles. Under conditions of low release probability (0.4 mM [Ca2+]e) paired-pulse stimulation [Fig. [Fig.1A]1A] led to prominent facilitation, which subsided at higher pr (1.0 mM [Ca2+]e) due to simultaneously occurring depression. Assuming that facilitation is based on a build-up of residual Ca2+ in the nerve terminal, the time course of paired-pulse facilitation (PPF) in low [Ca2+]e (i.e., with little contribution by depletion) should correspond to the decay of the residual Ca2+ concentration. A mono-exponential fit to the PPF revealed a time constant of 57 ms, which is in good agreement with the previously recorded decay of residual Ca2+ (τ=60 ms, Macleod et al., 2002). The models of facilitation used in this study therefore assume that residual Ca2+ decays with τ=60 ms.
Next, we measured the amplitude maintained under steady-state conditions at different stimulation frequencies [Fig. [Fig.1B].1B]. In high [Ca2+]e depression became dominant at increased frequencies while under conditions of low pr high-frequency stimulation was capable of increasing the steady-state level (SSL) of EPSC amplitudes via elevated facilitation.
Finally, a complex stimulation protocol was investigated, consisting of a train of 100 pulses at 60 Hz followed by test pulses at increasing intervals [Wu et al., 2005; Fig. Fig.1C].1C]. During the 60 Hz, the tetanus residual glutamate likely accumulates in the synaptic cleft, thereby mediating a “tonic” current [see, e.g., Sakaba (2006)]. In this study we concentrated on the “phasic” component and used the peak EPSC amplitude as a measure for synchronous release. We determined time courses for the facilitation and depression during the train and for the recovery thereof. EPSC amplitudes in high and low [Ca2+]e were similar after the first ~10 pulses [Fig. [Fig.1C].1C]. Recovery of EPSCs after the train critically depends on [Ca2+]e. In 1.0 mM [Ca2+]e, a biphasic recovery process was discernible (Wu et al., 2005), characterized by a fast initial component during which EPSC amplitudes increased almost twofold within the first ~100 ms and a slow subsequent phase lasting ~50 s. In contrast, low [Ca2+]e led EPSC amplitudes to decrease during the first 100 ms after the train, which is most likely attributable to the decay of residual Ca2+ in the terminal and with that, cessation of facilitation.
The observed STP results from the complex interplay of facilitation and depression, which in turn are both subject to specific constraints. Thus for example, a model of facilitation must be capable of describing the twofold PPF at short intervals in low [Ca2+]e [Fig. [Fig.1A].1A]. We focused on implementing a mechanistic model for facilitation proposed by Neher and co-workers (Neher, 1998; Trommershäuser et al., 2003), which is based on detailed investigations of Ca2+ fluxes and Ca2+ buffer saturation in the vicinity of Ca2+ channels. This work provides estimations of Ca2+ concentrations by assuming that microdomains arising from individual Ca2+ channels add up linearly (Neher, 1998). While several parameters of this model were experimentally determined at the calyx of Held synapse, here, two (or three) free parameters were implemented to allow faithful descriptions of the experimental data: (1) the parameter α is related to the distance between Ca2+ channels and the Ca2+ sensor for vesicle fusion and thereby critically determines pr (a large value for α corresponds to a high pr). (2) The parameter γ characterizes the saturation of Ca2+ buffers and thus sets the strength of facilitation. (3) If heterogeneity in pr is assumed, the parameter η defines the difference in pr between vesicle populations (Trommershäuser et al., 2003).
The model of facilitation was then combined with several frameworks describing the kinetics of vesicle trafficking to active zone release sites. Six of such “vesicle pool models” are illustrated in Fig. Fig.2A.2A. The free parameters of these concepts (including the free parameters of the facilitation model) were now systematically varied to simultaneously describe PPF, SSL, and STP during and after the train in both [Ca2+]e. The predictions of the models are superimposed on the experimentally recorded data in Figs. Figs.2B,2B, ,2C,2C, ,2D2D and the best-fit parameters are given in Supplementary Table S1 (see Supplementary Figs. S1 and S3 for more details and further models). The deviation of the model predictions is expressed as the difference between the data and the predictions (χ2) and is plotted in Fig. Fig.2E2E for all sets of stimulation protocols. As illustrated in Fig. Fig.2F,2F, vesicle pool model 6 (black) delivered the best description of the experimental data using only six free parameters.
One RR-pool: We started with a model based on the idea that a number of vesicles are readily releasable (RR) and that after release, new RR vesicles are supplied from a reserve pool of vesicles (Fig. (Fig.2,2, model 1, blue; Liley and North, 1953; Elmqvist and Quastel, 1965; Betz, 1970). The release probability of a vesicle in the RR-pool changes during the train according to the implemented model of facilitation (Trommershäuser et al., 2003). The best description of the data was obtained with 262 RR vesicles (NRR, see Supplementary Table S1) with a pr of 0.41 in 1.0 mM [Ca2+]e, a recovery time constant of 83 ms (τRR), and a facilitation parameter, γ, of to 2.74 μM−1. This resulted in an increase in pr from 0.41 to 0.97 based on the rise of residual Ca2+ from 0.1 μM to ~1 μM during the train (see Supplementary Fig. S1 and Table S1 for model rates). After the initial pulses of the train, the number of vesicles released per action potential equals the number of vesicles newly supplied in between stimuli. On average 30 RR vesicles exist in this steady-state condition (see Supplementary Fig. S1A).
One RR-pool+CDR: It has been suggested that the supply of vesicles is accelerated after stimulation (Worden et al., 1997; Wang and Kaczmarek, 1998; Sakaba and Neher, 2001a; Hosoi et al., 2007). Especially the rate of vesicle mobilization from the reserve pool to the RR-pool was assumed to be accelerated when residual Ca2+ is increased (Weis et al., 1999; Dittman et al., 2000; Zucker and Regehr, 2002; Hosoi et al., 2007). Consistently, introducing such Ca2+ dependent recovery (CDR) led to an improved description of the STP found at the Drosophila NMJ (Fig. (Fig.2,2, model 2, pink). The best description was obtained with 386 RR vesicles that recover with a time constant of 5.2 s under resting conditions (see Supplementary Table S1). During the train, residual Ca2+ accumulates and the mobilization rate increases, elevating the steady-state level. After the train, the mobilization rate drops back to its resting value, resulting in slow refilling of the RR-pool.
One RR-pool+a finite supply pool: The implemented infinite reserve pool is used as an approximation of a large pool of vesicles whose depletion does not limit the supply of vesicles to the release sites. Abandoning this simplification leads directly to model 3, which incorporates an intermediate supply pool (Fig. (Fig.2,2, green). This model was recently applied to cerebellar mossy fiber bottoms (Saviane and Silver, 2006). At the Drosophila NMJ the best description of the data was obtained with 224 RR vesicles, which are supplied with a time constant τRR of 20 ms by 6941 supply vesicles (Nsuppl), which in turn, are refilled with a time constant τsuppl of 9 s from the reserve pool. The delayed depression of EPSC amplitudes during the train, which is especially apparent in logarithmic plots, can be interpreted with this model. Since each vesicle in the supply pool is able to refill a depleted RR vesicle, a decreased number of vesicles in the supply pool leads to a proportionally slowed supply of vesicles to the RR-pool. In between stimuli this results in fewer RR vesicles available per action potential. As the intervals between stimuli increase after the end of the train, the RR-pool is refilled quite rapidly from the only partially depleted supply pool. In summary while the depression during the train is captured accurately, the fast monophasic recovery predicted by this model is obviously incorrect.
Two RR-pools: A possible interpretation of the observed biphasic recovery could be that two populations of RR vesicles exist, which recover with different kinetics (Fig. (Fig.2,2, model 4, red). In fact, the calyx of Held synapse appears to harbor one population of vesicles that posses a high pr but recover slowly and another population that has a low pr and recovers rapidly (Wu and Borst, 1999; Neher and Sakaba, 2001; Sakaba and Neher, 2001b; Sakaba and Neher, 2001a). At the Drosophila NMJ the data are described best by a population of 181 vesicles that recover with a time constant of 44 ms and a population of 163 vesicles that recover with a time constant of 3.7 s. During the train, the slowly recovering vesicles are completely depleted and transmission is sustained by the rapidly recovering pool. However, the delayed depression during the train is not captured by this model and, in contrast to the expectation derived from the findings at the calyx, the data demand similar release probabilities (pr=0.44 and pr=0.36) for both populations of vesicles.
Two RR-pools+a finite supply pool: Our aforementioned results indicated that heterogeneous vesicle populations could explain the biphasic recovery after the train and that a finite supply pool of vesicles captured the build-up of depression during the train. We therefore investigated a combination of both mechanisms (Fig. (Fig.2,2, model 5, brown). It turned out that while this model was capable of describing all of the experimental findings (see also Supplementary Fig. S1 and Table S1), some of its properties seemed suspicious. (1) The model predicts a very small population of 82 RR vesicles with a high pr. Estimates from anatomical investigations suggest that about 500 synapses mediate the EPSCs studied here (Atwood et al., 1993), which implies that only 16% of the synapses would posses a vesicle from this population. (2) The predicted pr (~0.80) of these vesicles is surprisingly high and the remaining vesicles have a pr of ~0.40, which is similar to the estimated pr in the other models (see Supplementary Table S1). According to the results obtained at the calyx of Held, we would expect the high pr pool to posses a release probability of ~0.5 and that these vesicles should dominate the results acquired through fluctuation analysis under resting conditions (Fig. (Fig.5).5). The low pr pool was expected to have a much lower release probability (pr<0.1) and to operate preferentially during the train. We therefore addressed whether simpler models might also be able to explain our data appropriately.
One RR-pool in fast equilibrium with a finite supply pool: It emerged that our synaptic recordings could be described by a model consisting of an RR-pool, which is in fast equilibrium with a finite supply pool (Fig. (Fig.2,2, model 6, black, see also Zucker and Regehr, 2002). The most satisfactory fit to the data was obtained with 299 RR vesicles, 11 075 supply vesicles, and an equilibrium time constant of 47 ms. Both the time course of depression during the train, attributable to depletion of the supply pool, and the biphasic recovery thereafter are adequately captured. The initial fast recovery stems from the fast equilibrium between RR vesicles and the supply pool. Since vesicles are able to rapidly leave the RR-pool, its replenishment, in turn, depends on the refilling of the supply pool. Thus, the later phase of recovery follows the slow refilling kinetics of the supply pool.
In summary, the parameters of diverse models were quantified and several of the investigated models were able to reproduce the facets of STP recorded in high and low [Ca2+]e. Notably, several parameters appeared comparably robust despite the substantially different mechanisms assumed by the models (compare, e.g., NRR, Nsuppl, and pr of different models in Supplementary Table S1).
Next, we investigated several assumptions underlying the data presented up to this point. First, the implemented model of facilitation, which is based on detailed mechanistic interpretations of calcium dynamics in the terminal (Neher, 1998; Trommershäuser et al., 2003), uses parameters determined at the calyx of Held to describe facilitation at the Drosophila NMJ. In addition, we therefore applied an alternative, phenomenological model of facilitation developed at neocortical pyramidal neurons (Tsodyks and Markram, 1997; Markram et al., 1998). Using such an entirely different description of facilitation (see Materials and methods section) to drive vesicle pool models delivered essentially identical results (see Supplementary Table S1B).
Second, we tested how postsynaptic mechanisms of STP would influence our conclusions. To this end, a strong but realistic degree of postsynaptic depression was assumed with kinetic properties based on previous investigations. Again, similar results were obtained (Supplementary Table S1).
The analysis of STP critically depended on the information gained from the time course of recovery after depression [e.g., Fig. Fig.2D].2D]. The different pool models predict very different recovery behavior and only those models with a supply pool (models 3, 5, and 6) would predict an altered recovery after trains with more stimuli. Fig. Fig.3A3A shows the average EPSC amplitudes during and following trains with 1000 stimuli at 60 Hz (dark gray) and, for comparison, the 100 pulse train (light gray). Fitting model 6 to the train of 1000 stimuli including the recovery (black line) gives very similar parameters as estimated using the aforementioned protocols (see legend of Fig. Fig.33 for values). The more rapid recovery of the supply pool may be due to CDR of the supply pool from the reserve pool (Zucker and Regehr, 2002). Notably, model 6 accurately describes the mono-exponential nature of the recovery from 1000 stimuli, whereas model 4, for instance, would predict an unaltered biphasic recovery independent of the train duration.
As an additional approach to gain information about vesicle pools and their mobilization rates, a linear extrapolation of the cumulative EPSC amplitude was performed [Schneggenburger et al., 1999; Fig. Fig.3B].3B]. Focus on the first 30 EPSCs of a 60 Hz train suggested ~300 RR vesicles. However, when all 1000 EPSCs of a long 60 Hz train were taken into account ~5000 vesicles were estimated. Besides indicating that back-extrapolation of cumulative EPSC amplitudes should be interpreted cautiously (as the results depend heavily on the identified SSL), these data also clearly highlight that complex mechanisms (as, e.g., several pools of synaptic vesicles) are involved in the depression during trains of stimuli.
Finally, we addressed how well the parameters of a model were defined depending on the number of stimuli within the train. Focusing on model 6, EPSC amplitudes during 60 Hz trains of 10–1000 stimuli were simulated and subsequently used as data to extract the parameters of the model. The calculated errors of the parameters (square roots of the diagonal elements of the covariance matrix; Colquhoun et al., 2003) were expressed as coefficients of variation (CV) and indicate how well each parameter could be estimated. The values for the supply pool (especially τsupply) could be estimated only poorly from the depression during the train [Fig. [Fig.3C].3C]. However, when the same calculations were repeated with trains of stimuli including 14 pulses during recovery [as in the experiments lasting ~100 s; compare with Fig. Fig.1C]1C] all parameters could be determined very reliably and almost independently of the number of stimuli within the train (CV<100%; data not shown). These results indicate that the parameters of the model can be estimated reliably as long as the information contained in the time course of recovery is included.
Bruchpilot (BRP) is an integral component of the Drosophila active zone cytomatrix where it helps clustering presynaptic Ca2+ channels at vesicle fusion sites to ensure adequate release probability (Kittel et al., 2006a, 2006b; Wagh et al., 2006; Fouquet et al., 2009). The modeling study presented here was partially motivated by findings made at BRP null mutant synapses. In particular, the strongly altered recovery of BRP mutants after a train [Fig. [Fig.4A]4A] prompted us to develop quantitative models of active zone function that would allow us to interpret the complex observations.
Figures Figures4A,4A, ,4B,4B, ,4C4C (right column) shows the PPF, the SSL, and the train plus recovery recorded at BRP mutant neuromuscular synapses in 1.0 mM [Ca2+]e (blue). The train experiments were also performed in 2 mM [Ca2+]e [light gray, Fig. Fig.4A]4A] in order to give roughly the same pr as experienced by control animals in 1 mM [Ca2+]e [compare initial values of trains in Fig. Fig.4A].4A]. By applying model 6 to the STP data of BRP mutants and by optimizing all six free parameters, we were able to obtain good fits to the experimental observations [dark blue and dark gray lines in Figs. Figs.4A,4A, ,4B,4B, ,4C].4C]. Figure Figure4E4E summarizes the parameters used by the model to fit both controls (black dotted lines) and mutants (blue dotted lines).
Importantly, the modeling approach captured the previously reported low pr of BRP mutants (Kittel et al., 2006a). However, while the number of RR vesicles appeared unaltered, BRP mutant active zones were described as possessing a significantly smaller supply pool with slower refilling kinetics [compare with Fig. Fig.4E4E].
To confirm the reliability of the extracted parameters, we fitted single experiments of the train with subsequent recovery and estimated the variability of the parameters form one experiment to the next (Figs. (Figs.4D,4D, ,4E4E and Supplementary Fig. S2). The data were acquired separately under conditions of low and high release probability (0.4 mM and 1 mM [Ca2+]e for controls; 1 mM and 2 mM [Ca2+]e for BRP mutants) and model 6 was utilized to fit the results. The available information was now strongly reduced compared to the previous approach of simultaneously fitting the averages of PPF, SSL, and train plus recovery in both Ca2+ concentrations. Nevertheless, in 1 mM [Ca2+]e the single trains sufficed to extract most of the model parameters. The parameter γ critically depends on PPF and was therefore fixed to the value estimated previously (2.74 μM−1). In addition, vesicle release in 0.4 mM [Ca2+]e is so moderate that it was not significantly limited by the size and the kinetics of the supply pool. For single experiments in low [Ca2+]e the number of supply vesicles and their refilling kinetics were therefore fixed to the previously obtained values (Fig. (Fig.4E,4E, see Supplementary Table S1). The reproducibility of the individual parameters [open and filled symbols in Fig. Fig.4E]4E] supports the validity of the investigated model (Colquhoun and Sigworth, 1995).
In addition, we tested model 2 (Wang and Kaczmarek, 1998; Hosoi et al., 2007) and assumed that the miniature amplitude was decreased by 1% per action potential with a recovery time constant of τ=1.5 s, resembling, e.g., postsynaptic receptor saturation, desensitization, or Ca2+-dependent Ca2+ channel inactivation (Forsythe et al., 1998; Hennig et al., 2008; see Supplementary Fig. S3). Although this model assumes different mechanisms it also highlighted that the BRP mutant phenotype cannot be accounted for by lower release probability alone (compare the changes in k and k1,b in Supplementary Fig. S3)
Fluctuation analysis can be used to determine the three quantal parameters N, p, and q under resting conditions (Silver et al., 1996; Quastel, 1997; Meyer et al., 2001; Scheuss and Neher, 2001) and independently of the modeling approach. To increase the signal to noise ratio, all EPSCs used forfluctuation analysis were recorded with a focal electrode [Fig. [Fig.5B],5B], whereas the two-electrode voltage clamp method had been applied for STP protocols to prevent muscle contractions [Fig. [Fig.5A].5A]. Elevating [Ca2+]e increased the peak amplitude of EPSCs such that similar values were reached by BRP mutants and controls in high [Ca2+]e [Figs. [Figs.5C,5C, ,5D].5D]. The variance of EPSC amplitudes rose with increase in [Ca2+]e, reached a maximum at intermediate values, and then decreased again. To obtain the number of RR vesicles (N), the quantal size (q), and the release probability (pr), mean-variance plots were generated and stationary fluctuation analysis was used [Fig. [Fig.5E].5E]. Consistent with earlier work (Kittel et al., 2006b) and confirming the modeling results, pr was reduced to 0.08±0.02 in BRP mutants compared to 0.29±0.06 in controls (P<0.01) while N was unchanged [Fig. [Fig.5F].5F]. The absolute value obtained for N is smaller in focal recordings because the two-electrode voltage clamp technique (TEVC), applied for STP protocols samples a larger number of synapses (Fig. (Fig.5;5; see also Materials and methods section). Furthermore, the absolute values obtained for pr with fluctuation analysis are smaller compared to the values determined through STP experiments because the fluctuation analysis was performed at 18 °C (to prevent muscle contraction in the focal recording configuration) and not at 22 °C as for the TEVC recordings. Assuming a temperature coefficient Q10 of 4 (Hille, 2001), a pr of 0.29 at 18 °C would correspond to a pr of 0.50 at 22 °C. This is in good agreement with the model estimates based on TEVC recordings (see Supplementary Table S1).
To gain insight into the mechanisms underlying short-term plasticity at neuromuscular active zones of Drosophila, we compared several mechanistic models and quantified their parameters. Through utilization of these models we were then able to further analyze the function of the active zone protein BRP. Finally, the basic parameters under resting conditions were confirmed independently for both control animals and BRP mutants via fluctuation analysis.
In this study, we focused on the peak EPSC amplitude and related it linearly to the number of vesicles released per action potential. While linear summation of quanta has been demonstrated over a wide range of release probabilities at, e.g., cerebellar mossy fiber synapses (Sargent et al., 2005), a number of mechanisms could corrupt this assumption. For instance, increased asynchronous release may change the time course but not necessarily the peak amplitude of EPSCs (see, e.g., Sakaba 2006). However, deconvolving EPSCs during 60 Hz trains with templates of recorded miniature EPSCs (single quanta) resulted in very similar time courses of vesicle release per stimulus (data not shown). Also, the similar release probabilities estimated for both RR-pools of models 4 and 5 is consistent with recent data, which indicates that asynchronous release contributes little to the peak EPSC amplitude at the calyx of Held (Sakaba, 2006).
Furthermore, Ca2+-dependent Ca2+ channel modulation (Forsythe et al., 1998; Hennig et al., 2008; Mochida et al., 2008) and action potential broadening have been reported to influence STP (Geiger and Jonas, 2000). It is thus also conceivable that the increased pr during the train may be influenced by such mechanisms in addition to buffer saturation (Trommershäuser et al., 2003). Postsynaptic depression is an obvious mechanism that contributes to STP (Scheuss et al., 2002; Zucker and Regehr, 2002; Taschenberger et al., 2005; Scheuss et al., 2007). Therefore, both fluctuation analysis in high [Ca2+]e (Fig. (Fig.5)5) and investigations of SSLs and trains (Figs. (Figs.1112223334)4) could be distorted by, e.g., postsynaptic receptor desensitization or saturation. This would underestimate the number of RR vesicles and thereby overestimate pr. However, implementing strong postsynaptic depression in our simulations did not affect the presented conclusions (see Supplementary Table S1C, Materials and Methods Section). Since the time course of recovery differs after 100 versus 1000 pulses [Fig. [Fig.3A]3A] it appears unlikely that postsynaptic mechanisms elicit the biphasic recovery after 100 pulses (as by this time a steady-state situation concerning the glutamate concentration in the cleft and receptor saturation and desensitization should already have been attained (Adelsberger et al., 1997; Heckmann and Dudel, 1997; Pawlu et al., 2004). Instead, it appears more likely that a steady-state is not reached after 100 stimuli because the depletion of vesicles has not yet been completed.
Previous work has described variations in pr across active zones formed on the same axon or innervating the same target (Wojtowicz et al., 1994; Murthy et al., 1997; Marrus and DiAntonio, 2004; Wairkar et al., 2009). Similarly, our estimations of pr most likely reflect a mean value given by a broad distribution. As the neuromuscular preparation is innervated by two motoneurons differing in their physiological properties (Kurdyak et al., 1994), it will be interesting to examine whether these contribute heterogeneously to vesicle pool models.
Taking into account how many active zones (AZs) are present in this preparation (~500, as estimated from serial electron-microscopic images of larval longitudinal muscle 6 (Atwood et al., 1993), our estimated rate of vesicle recruitment during sustained release is ~15/s per AZ (Nsupply×kRR=10,000×0.7/s per 500 AZs (model 6) and k1,b+ΔCa2+×k=0.2/s+1 μM×15/(μM s) per AZ [modified model 2 in Supplementary Fig. S3)]. Hence, an AZ requires about 70 ms to recruit a vesicle during sustained release. This is faster than previous calculations at the Drosophila NMJ (500 ms, Delgado et al., 2000; Kidokoro et al., 2004) but within the range of recent estimates at other synapses (12 ms, Saviane and Silver, 2006), 250 ms (Hosoi et al., 2007). Our model 6 assumes a fast equilibrium between RR and supply vesicles. Such a mass-action model was previously suggested for the Drosophila NMJ in order to account for the delayed depression during the train (Li and Schwarz, 1999). However, the estimated undocking rate of RR vesicles at central synapses (Murthy and Stevens, 1999) is lower than the high transition rate of vesicles from the RR to the supply pool described by model 6 (21.1 s−1).
While these differences may be reconciled by distinct functional specializations present at different synapses, the aim of our modeling approach was to generate experimentally testable hypotheses for further analysis. For each investigated model a set of parameters was determined and several of these estimates appeared relatively independent of the underlying vesicle pool model (compare, e.g., NRR, Nsuppl, and pr of different models in Supplementary Table S1). Together with additional experimental confirmations, such as through fluctuation analysis, this increases the confidence in the parameter estimates. Finally, both models 6 and 2 identified defective trafficking of vesicles to BRP mutant active zones (Fig. (Fig.44 and Supplementary Fig. S3). Examining the precise molecular mechanisms underlying this supply defect will be of great interest.
As control animals we used the genotype w1 and the genotype of Bruchpilot null mutants was brp69/df(2R)BSC29 (Kittel et al., 2006b).
TEVC recordings of EPSCs were obtained at room temperature from late third instar male Drosophila larvae (ventral longitudinal muscle 6, segments A2 and A3, and Fig. Fig.5),5), essentially as previously described (Kittel et al., 2006b). For the TEVC experiments, both mutants and control animals carried a copy of either elav-GAL4 or ok6-GAL4 to ensure comparability with a previous study (Kittel et al., 2006b). The composition of the extracellular hemolymph-like saline (HL-3, Stewart et al., 1994) was (in mM): NaCl 70, KCl 5, MgCl2 20, NaHCO3 10, trehalose 5, sucrose 115, HEPES 5, CaCl2 as indicated, and pH adjusted to 7.4. Recordings were made from cells with an initial Vm between −50 mV and −70 mV (holding potential at −60 mV) using intracellular electrodes with resistances of 10–32 MΩ, filled with 3 M KCl. Train stimulation protocols consisted of either 100 or 1,000 pulses applied at 60 Hz. The recovery was assessed by evoking APs at (in ms following the last pulse in the train): 25, 50, 100, 200, 500, 1,000, 2,000, 5,000, 10,000, 20,000, 50,000, and 100,000 (Wu et al., 2005). EPSCs reflect the compound response to stimulation of both motoneurons innervating muscle 6, and care was therefore taken to ensure their stable recruitment. Infrequently observed recruitment failures were linearly interpolated. Only cells that recovered at least 70% of their initial EPSC amplitude following tetanic stimulation were included in the analysis. The recordings were analyzed with pClamp 9 (now Molecular Devices, Sunnyvale, CA, USA) and the peak amplitude was determined as the difference between the peak value of the EPSC and the baseline value before onset of that EPSC. The amplitude of the second response in 10 ms interpulse recordings was measured from the peak to the point of interception with the extrapolated first response.
EPSCs for fluctuation analysis were recorded with a focal electrode from third instar larvae (Fig. (Fig.5),5), in HL-3 essentially as reported (Pawlu et al., 2004). For each [Ca2+]e, more than ~20 EPSC were recorded at 0.2 Hz, after a steady-state amplitude was reached. The bath temperature was kept within 18 °C±0.5 °C with a Peltier element (27 W, Conrad Electronic, Hirschau, Germany) glued to the bath inflow with heat-conductive paste (Fischer Elektronik, Lüdenscheid, Germany). Focally recorded EPSCs were analyzed with Igor Pro 5.04 (WaveMetrics, Lake Oswego, OR). EPSCs were digitally filtered at 3 kHz (Gaussian filter), baseline subtracted, and the average of all failures was subtracted from the currents. Whether action potentials had been initiated was assessed by the presence of the current resulting from action potential invasion of the nerve terminal (Dudel, 1981). The peak amplitude of the EPCSs was determined within a 5-ms-window after the presynaptic action potential. The average peak amplitude of the failures was subtracted from the peak amplitude of the EPSCs.
To plot the relationship between EPSC amplitudes and [Ca2+]e, [Fig. [Fig.5D],5D], the EPSC amplitudes of each experiment were first normalized to the amplitude in 1 mM [Ca2+]e. The normalized values were then averaged across experiments and subsequently the average was rescaled with the average amplitude of all experiments in 1 mM [Ca2+]e. The data were fitted with a Hill equation according to I([Ca2+]e)=Imax[1+(EC50/[Ca2+]e)slope]−1, where Imax=10.6 and 9.5 nA, EC50=1.5 and 2.2 mM, and slope=2.1 and 2.8 for controls and BRP mutants, respectively.
To calculate the variance (and the variance of the variance) we followed the procedures described by Neher and co-workers (Meyer et al., 2001; Scheuss and Neher, 2001; Scheuss et al., 2002). For a set of N peak amplitudes xi (i=1,…,N) recorded at a certain [Ca2+]e, the sample variance was calculated as the average of the sample variance of overlapping segments of size two. The variance of the variance was calculated as the standard error of the mean (SEM) of all sample variances (Meyer et al., 2001).
Stationary EPSC fluctuation analysis was performed as formerly reported (Meyer et al., 2001; Scheuss and Neher, 2001; Silver, 2003). Based on a binomial model for synaptic transmission, the variance of the EPSC amplitude depends on the quantal size (q), the number of release sites (N), and the probability of release at each site (pr), resulting in a parabolic dependence of the EPSC variance on the average EPSC amplitude. The variance-mean plots were fitted with and weighted with the reciprocals of Var(Var) and constrained to pass the origin. To correct for the variability in q, we corrected for intrasite variance (CVQI) and intersite variance (CVQII). The apparent and estimates are related to the true parameters by (Frerking and Wilson, 1996; Meyer et al., 2001; Scheuss and Neher, 2001; Silver, 2003)
The experimentally determined variance of q at theDrosophila NMJ [Pawlu et al., 2004; CVQ=0.30=(CVQI2+CVQII2)0.5] in combination with the assumption of CVQI=CVQII (Frerking and Wilson, 1996; Meyer et al., 2001; Scheuss and Neher, 2001; Silver, 2003) results in and .
The simulations were based on a biophysical interpretation of facilitation described in Trommershäuser et al. (2003). In brief, a linear build up of residual calcium (ΔCa2+) was assumed, where each AP leads to a constant increase in ΔCa2+ by x0=0.4 μM and decays mono-exponentially with τx=60 ms. Instead of 100 ms (Helmchen et al., 1997; Trommershäuser et al., 2003) we used 60 ms since (1) recent data indicate τx≈30 ms for undialysed Calyces (Müller et al., 2007); (2) measurements at the Drosophila NMJ indicated τx=60 ms (Macleod et al., 2002); and (3) a mono-exponential fit of the paired pulse facilitation versus interval duration recorded in 0.4 mM extracellular Ca2+ concentration ([Ca2+]e) gave τ=57 ms.
Two types of readily releasable vesicles can be considered (Trommershäuser et al., 2003): vesicles close to Ca2+ channels (pool2) and vesicles further away from Ca2+ channels (pool1). The corresponding peak local Ca2+ concentration at a release site of pool j, [Ca2+]j, can be expressed as (Trommershäuser et al., 2003):
where δj,2 is the Kronecker symbol, which has the value 1 for j=2 and 0 for j=1. The local [Ca2+] at the release site is related to the probability of release according to a fourth power relationship with a half-maximal release at K1/2=42.5 μM (Trommershäuser et al., 2003). The global Ca2+ concentration ([Ca2+]gl) is elevated above the resting Ca2+ concentration ([Ca2+]r) by the residual Ca2+ concentration. The Ca2+ influx (J) depends on [Ca2+]e as described by a Michaelis–Menten saturation equation
with half-maximal Ca2+-influx at an EC50 of 2.61 mM [Ca2+]e and a dimensionless quantity Jmax=2.31 (Schneggenburger et al., 1999; Trommershäuser et al., 2003). In addition, the increase per AP (x0) was scaled according to
The parameter α is proportional to the inverse of the distance between Ca2+ channels and the vesicles (Trommershäuser et al., 2003) and defines the resting release probability. The parameter η determines the difference in release probability between pool1 and pool2. The parameter γ quantifies facilitation.
The free parameters (α, η, and γ) were used to fit the experimental data of control animals including paired-pulse plasticity, frequency-dependent steady-state levels, and the response to and recovery from a train of 100 pulses at 60 Hz, in 0.4 mM and 1.0 mM [Ca2+]e (Fig. (Fig.2).2). In models with only one pool of readily releasable vesicles (e.g., model 6), only the vesicles close to Ca2+ channels (pool2) were simulated (compare Sakaba, 2006; Hosoi et al., 2007) and η was fixed to the value in Trommershäuser et al. (2003). Due to the degree of facilitation required to reproduce the twofold paired pulse facilitation at short intervals in low [Ca2+]e, γ was (2.74 μM−1). When individual experiments (train plus recovery) were fitted [Fig. [Fig.4D4D and Supplementary Fig. S2], γ was fixed to 2.85 μM−1 because the paired-pulse ratio in low Ca2+ critically constrained the degree of facilitation. For fitting individual train-experiments in 0.4 mM [Ca2+]e, in addition, Nsupply, τsupply, and τRR (model 6, Fig. Fig.4D,4D, and Supplementary Fig. S2) and α, k1.b, and k (model 2, see Supplementary Fig. S3D) were fixed to the values obtained from the average data [Fig. [Fig.4E4E and Supplementary Fig. S2].
Starting with analytically calculated steady-state values, the number of vesicles in each pool was numerically calculated with a fifth-order Cash–Karp Runge–Kutta method with adaptive stepsize control (Cash and Krarp, 1990). For computational speed, global calcium [Ca2+]gl(t) was calculated as an array with a temporal resolution of 100 μs. For steady-state EPSC amplitude, 10, 30, 50, 70, 100 and 150 stimuli were simulated at a frequency of 0.1 Hz, 0.3 Hz, 1 Hz, 3 Hz, 10 Hz, and 60 Hz, comparable to the experiments. The last EPSC of the simulated sweep was taken as the steady-state level at the corresponding frequency. The predicted number of released vesicles was converted to a current assuming a miniature EPSC amplitude of 0.63 nA for controls and 0.87 nA for BRP mutants (based on 0.84 nA for controls and 1.17 nA for BRP mutants recorded at −80 mV holding potential (Kittel et al., 2006b) and scaled to −60 mV, assuming 0 mV reversal potential and a linear current voltage relation).
For parameter optimization, all free parameters were varied with a simplex algorithm as described (Press et al., 2002). The deviations of the model predictions from the experimental data were measured as the sum of squared differences, χ2, as follows: first, the number of data points obtained from the train and recovery experiments (100 during the train+13 from the recovery) was reduced by exponential sampling to 28 values (15 during the train+13 from the recovery). Second, the six values from the steady-state and the five values from paired-pulse plasticity were weighted fivefold. Finally, the χ2 values from the steady-state and train plus recovery were scaled by the inverse of the squared average current amplitude under resting conditions. This resulted in six χ2 values of similar size for the paired-pulse plasticity, the steady-state level, and the train plus recovery at both 0.4 mM and 1.0 mM [Ca2+]e [compare with Fig. Fig.2E2E].
Subsequently, the model predictions and the experimental data were plotted superimposed using MATHEMATICA 5.0 (Wolfram Research, Champaign, IL, USA) for visual inspection. Simulation and plotting lasted only several seconds on a conventional 2.7 GHz dual core computer, which allowed us to gain a better understanding of the reliability and the constraints of the parameters. After the fit was completed, internal estimates of the errors of the best-fit parameters were calculated as the square roots of the diagonals of the covariance matrix [Colquhoun et al., 2003; compare with Fig. Fig.3C].3C]. To confirm that the global minimum was reached, the best-fit parameters were shown to be independent of the starting values within a plausible range. In addition, the reproducibility of the best-fit parameters to individual experiments, as performed in the analysis of single channel recordings (Colquhoun et al., 2003; Hallermann et al., 2005), indicated that the global minimum was reached (Colquhoun and Sigworth, 1995).
The nonparametric Mann–Whitney rank sum test was used for statistical analysis. The data are reported as mean±SEM, n indicates the sample number and P denotes the significance (*P<0.05, **P<0.01, and ***P<0.001).
We thank Jens Eilers and Stephan J. Sigrist for scientific exchange. This work was supported by grants from the DFG to M.H. (Contract Nos. HE 2621/4-2 and TP B27/SFB 581) and by formel.1 grants to S.H. and R.J.K. from the Medical Faculty of the University of Leipzig. R.J.K. is supported by the DFG Emmy–Noether Program.