The spatial arrangement of Ca2+ channels and vesicles remains unknown for most CNS synapses, despite of the crucial importance of this geometrical parameter for the Ca2+ control of transmitter release. At a large model synapse, the calyx of Held, transmitter release is controlled by several Ca2+ channels in a" domain overlap" mode, at least in young animals. To study the geometrical constraints of Ca2+ channel placement in domain overlap control of release, we used stochastic MCell modelling, at active zones for which the position of docked vesicles was derived from electron microscopy (EM). We found that random placement of Ca2+ channels was unable to produce high slope values between release and presynaptic Ca2+ entry, a hallmark of domain overlap, and yielded excessively large release probabilities. The simple assumption that Ca2+ channels can be located anywhere at active zones, except below a critical distance of ~ 30 nm away from docked vesicles (" exclusion zone" ), rescued high slope values and low release probabilities. Alternatively, high slope values can also be obtained by placing all Ca2+ channels into a single supercluster, which however results in significantly higher heterogeneity of release probabilities. We also show experimentally that high slope values, and the sensitivity to the slow Ca2+ chelator EGTA-AM, are maintained with developmental maturation of the calyx synapse. Taken together, domain overlap control of release represents a highly organized active zone architecture in which Ca2+ channels must obey a certain distance to docked vesicles. Furthermore, domain overlap can be employed by near-mature, fast-releasing synapses.
Transmitter release at CNS synapses happens at active zones of sub-micrometer dimensions, which harbor docked vesicles and vesicle fusion proteins, as well as presynaptic scaffold proteins and voltage-gated Ca2+ channels [1]. The number, and distance of Ca2+ channels to readily-releasable vesicles are crucial for determining release probability, because the Ca2+ signal generated by a single open Ca2+ channel drops off steeply with distance [2–5]. There is, however, only sparse morphological information on the co-localization of individual vesicles and Ca2+ channels at synapses. For this reason, Ca2+ channel—vesicle coupling distances have often been inferred from functional measurements. One approach to functionally assess the number of Ca2+ channels controlling release utilizes the high intrinsic Ca2+ cooperativity of release [6–8]. In experiments in which the number of open Ca2+ channels is varied in presynaptic voltage-clamp experiments, the slope value in plots of transmitter release vs integral Ca2+ influx in double-logarithmic coordinates (which we will call the" Ca2+ current release cooperativity" or simply" slope value" ) can inform about the number of Ca2+ channels involved in release control [9–13]. A large Ca2+ current release cooperativity close to the intrinsic Ca2+ cooperativity indicates that Ca2+ signals from individual channels mix to produce a graded Ca2+ signal at each docked vesicle. Conversely, a slope value close to one indicates that a few, or a single Ca2+ channel control the release of a given docked vesicle [9–13]. The diffusional distance between Ca2+ channels and vesicles can also be probed by the slowly binding Ca2+ chelator EGTA, which affects Ca2+ signals only at some distance from the Ca2+ source [14–17]. Experiments with these functional approaches have shown different coupling regimes at different CNS synapses. At the calyx of Held synapse which is amenable to direct presynaptic patch-clamp experiments, release is sensitive to EGTA at low millimolar concentrations, and high slope values are observed, indicating domain overlap control of release [18,19]. However, Ca2+ channel—release coupling was later shown to be developmentally regulated, such that brief presynaptic Ca2+ currents become more effective in causing transmitter release [12,20–22]. Electron microscopy (EM) has shown that active zones of various types of CNS synapses are small (~ 0. 1 μm2), and contain a high density of docked vesicles [23–26]. The development of EM freeze fracture replica labeling techniques has enabled the visualization of ion channels and receptors in en-face views at nanometer resolution [27]. Using this method, multiple stripe—like clusters of P/Q—type Ca2+ channels (CaV2. 1 subunits) have been observed at hippocampal synapses [28]. However, with freeze fracture replica labeling, the position of docked vesicles cannot be visualized, and antibody-labelled particles do not necessarily overlap with the entire population of functionally active Ca2+ channels. Therefore, despite important advances in ultrastructural methods, the functionally relevant positions of Ca2+ channels and vesicles on the nanometer scale are still uncertain. Here, we use stochastic MCell modelling of Ca2+ influx through individual Ca2+ channels at EM—reconstructed active zones, to explore how Ca2+ channels might be positioned to produce domain overlap control of release. The position of docked vesicles was fixed by previous EM reconstructions of single active zones [29]. Other model parameters, like presynaptic Ca2+ channel gating and—conductance [18,19,30–32], Ca2+ buffering [33–35], and the intracellular Ca2+—sensitivity of vesicle fusion [7,8, 22,36–39] were constrained by previous measurements at the calyx of Held and at other large model synapses. We find that random placement of Ca2+ channels leads to an excessively high release probability and low slope values. In order to enable domain overlap control of release, we rather need to assume that Ca2+ channels are kept at some distance from vesicles. We also show experimentally that the Ca2+ current—release cooperativity stays high after the onset of hearing in mice, suggesting that despite a characteristic developmental tightening in the Ca2+ channel vesicle co-localization [12,20–22], fast release in more mature calyx synapses continues to be controlled by several Ca2+ channels.
We wished to explore the possible placement of Ca2+ channels compatible with domain overlap control of release. For this, it is important to know the distribution, and density of docked vesicles typically observed at active zones. To determine these parameters, we analyzed a sample of n = 15 reconstructed active zones of calyx synapses of a P11 wild-type mouse [29] (Fig 1A). In this sample, the average surface area of active zones was 0. 07 ± 0. 03 μm2, with an average density of docked vesicles of 110 ± 40 ves / μm2 (n = 15). Across individual active zones, the number of docked vesicles and active zone surface were correlated (r = 0. 74; Fig 1B), indicating that large active zones tend to harbor more docked vesicles [20,23]. Visual inspection of the active zone maps with their docked vesicles (Fig 1A) suggests that some active zones show sub-areas where docked vesicles are sparse. We will show below that in domain overlap control scenarios, such areas might be preferentially occupied by Ca2+ channels. To investigate whether such void spaces could arise randomly, we compared experimentally observed vesicle distributions with the same parameters derived from random x-y placements of vesicles at a given active zone (Fig 1C, 1D, 1E, 1F, 1G and 1H). Analyzing nearest neighbor distances, and the largest hole radii capable to fill empty spaces within active zones suggested that the observed vesicle positions were close to the random case (Fig 1C and 1D, Fig 1E and 1F, respectively). Similarly, calculating the summed area extending 30 nm away from the edge of docked vesicles resulted in similar areas for randomly distributed vesicles, and for the docked vesicles found in the data set (Fig 1G and 1H). Finally, we used Ripley’s K function test, a metric for detecting deviations from spatial homogeneity [40] (see S1 Text). It gave values within the 99% confidence interval for all distances tested (Fig 1I), indicating no significant vesicle clustering in the reconstructed active zones. Thus, there was overall no strong vesicle clustering in the reconstructed active zones. We next investigated whether random placement of Ca2+ channels in a realistic active zone with several docked vesicles can predict release control by several channels. We developed an MCell based model which incorporated individual Ca2+ channels with realistic gating and permeation properties [19,31,41] (Fig 2A). The model contained several intracellular Ca2+ buffers, and each vesicle had a Ca2+ sensor for vesicle fusion according to a highly non-linear five-site ion binding mechanism [7,8] (Fig 2A). The parameters of the Ca2+ sensor model were determined in several independent Ca2+ uncaging studies [7,8, 22,36,37], and were fixed to the values reported in one previous study [37]. Several other parameters, including the width of the presynaptic AP (Fig 2B) were also tightly constrained by previous biophysical experiments at the calyx of Held synapse (see Materials and Methods, and S2 Text). The Ca2+ channel density was set at a value of 280 /μm2 (ref. [28,41]). For most simulations, we used single representative active zones drawn from the sample of reconstructed active zones. We started with active zone #5 (see Fig 1A, star symbol) which has an area of 0. 05 μm2 (slightly smaller than the mean), and a vesicle docking density of 120 ves / μm2, representative of the sample (Fig 1B, vertical arrow). We placed n = 14 Ca2+ channels randomly in this active zone (Fig 2C), corresponding to a Ca2+ channel density of 280 /μm2. We simulated a" Ca2+ current—release cooperativity" experiment by driving release with APs of different widths (Fig 2D), and plotted the resulting average vesicle release probability, pves, as a function of Ca2+ entry in double-logarithmic coordinates (Fig 2E). This yielded a slope value of only 1. 3 ± 0. 2 (n = 5 independent random spatial samples), far smaller than a value of ~ 3–3. 5 expected for domain overlap control of release (see below). In addition, the average pves in response to the standard AP was 0. 66 ± 0. 18 (Fig 2F), much higher than experimental estimates (~ 0. 1; ref. [8]). As a control, we placed n = 16 Ca2+ channels in a tight cluster at 100 nm from a single vesicle. This resulted in a slope value of 2. 8 (Fig 2G), close to the maximal value of 3. 5 which can be achieved for kinetic models with 5 Ca2+ binding sites [11,42]. Thus, these simulations suggest that a random placement of Ca2+ channels is unable to produce high Ca2+ current—release cooperativities close to the experimentally observed values. It has been shown that a large number of Ca2+ channels (> 10–20) is needed to produce high Ca2+ current—release cooperativity [11]. If such a large number of channels controls release of a given vesicle with roughly equal strength, then it seems plausible to assume that the channels must be located at some distance to vesicles. A simple implementation of such a rule in a field of randomly placed docked vesicles at relatively high density (Fig 1) is to assume an exclusion zone around each docked vesicle; a zone into which Ca2+ channels cannot enter. We implemented this exclusion zone rule by assuming that Ca2+ channels can be located anywhere at the active zone, but not below a distance of 30 nm (Fig 3A). When the exclusion zone model was driven with a standard AP, brief local [Ca2+]i transients with amplitudes of 10–20 μM resulted at individual docking sites (Fig 3B), similar to results based on back-calculation from Ca2+ uncaging data [7,8]. The vesicular release probabilities pves were variable across individual docked vesicles, with a realistic average value of ~ 0. 1 (Fig 3C, upper panel). The distribution of the number of released vesicles showed that most trials (~ 60%) led to release failures at an active zone (Fig 3C), again in agreement with previous estimates [43]. We then simulated the Ca2+ current—release cooperativity experiment using different AP widths, and found that the exclusion zone model with a distance of 30 nm produced high values of Ca2+ current—release cooperativity (n = 2. 8 ± 0. 2; n = 4 independent spatial samples; Fig 3D). We attribute this to the fact that multiple Ca2+ channels (n = 6–14) should influence the release of a given vesicle. Thus, assuming that Ca2+ channels can be placed anywhere except at very close distances to docked vesicles (~ 30 nm) predicts a high Ca2+ current—release cooperativity, a hallmark of release control by several Ca2+ channels [9–13,22]. We also simulated the sensitivity of release to the slow Ca2+ buffer EGTA in the exclusion zone model (Fig 3E). For this, release in response to single standard APs was modelled in the presence of EGTA, added to the Ca2+ buffers present in the standard model (see Materials and Methods). These simulations showed that 0. 5 mM EGTA suppressed release by ~ half, with an overall biphasic concentration—dependence of EGTA (Fig 3E), similar to what was found at the calyx of Held synapse and at cortical synapses of young rats [15,18]. Finally, we tested the influence of the exclusion zone width on release probability pves, and found an inverse relation of pves with exclusion zone distance, with a value for pves of 0. 096 reached at 30 nm (Fig 3F, horizontal line). We conclude that the exclusion zone model can reproduce many features of domain overlap control of release reported in previous functional studies. In a previous model, release control by several Ca2+ channels was achieved by placing all channels in a single cluster [11], an arrangement which we will call" supercluster" . We next tested this arrangement, by first using the docked vesicle distribution of active zone #5. The distance of the supercluster to the nearest vesicle was 30 nm (Fig 4A), which yielded an acceptable average pves of 0. 08 (Fig 4B) and a high Ca2+ current—release cooperativity (2. 8 ± 0. 1; n = 4 independent spatial samples; Fig 4C). A characteristic property of the supercluster arrangement is, however, a very large heterogeneity between pves values for individual vesicles (range, 0. 005 to 0. 24; c. v. = 1. 1; Fig 4B), because far—away vesicles experienced a much lower Ca2+ signal. We next compared the exclusion zone—and supercluster channel arrangements at another example active zone with larger size (active zone # 11; see above Fig 1A, triangle symbol). This active zone also had a representative density of docked vesicles (see Fig 1B, horizontal arrow). In this active zone, 25 Ca2+ channels were used to again yield a Ca2+ channel density of 280 /μm2. Assuming an exclusion zone with a distance of 30 nm (Fig 5A), we found a somewhat higher average pves (0. 17; Fig 5B) as compared to the smaller active zone # 5 (see above; Fig 3; pves ~ 0. 1). This finding, which is compatible with recent experimental findings [41], might suggest that at comparable Ca2+ channel densities, large active zones use Ca2+ more efficiently because of smaller diffusional loss of Ca2+ away from the active zone. However, further modelling studies are needed to investigate this mechanism in more detail. In the simulations of Fig 5A–5C using an exclusion zone distance of 30 nm, the individual pves values ranged from 0. 02 to 0. 36 with a c. v. of 0. 7 (Fig 5B), and the Ca2+ current—release cooperativity was 2. 7. We next modelled Ca2+—release control at this large active zone assuming that all Ca2+ channels (n = 25) are located within a single supercluster. We placed the supercluster into a void space in the bottom half of the active zone, and varied its exact position until we achieved a reasonable average pves of 0. 1 (Fig 5D and 5E). The supercluster model again yielded a high slope value (3. 3; Fig 5F). However, the heterogeneity of pves was very large, with a c. v. value of 1. 4. Indeed, visual inspection of this relatively large active zone shows many vesicles that are far away from the cluster (e. g. vesicles 1,2, 7 and 8; Fig 5D). Correspondingly, these vesicles showed extremely low pves values (0,0. 0001,0. 0129, and 0. 012, respectively). As expected, there was an inverse relation between pves and distance of the docked vesicle from the Ca2+ channel supercluster (Fig 5G). We note, therefore, that the supercluster model produces a high heterogeneity of pves, especially in large active zones. Since the two arrangements of Ca2+ channels predict distinct heterogeneities of pves, experimental data on the heterogeneity of pves within individual active zones would be necessary to verify the supercluster versus the exclusion zone arrangement. Unfortunately, the distribution of pves within individual active zones following AP stimulation has, to our knowledge, not been estimated experimentally; a previous theoretical study has highlighted the role of pves heterogeneity for synaptic depression [44]. On the other hand, studies at the calyx using prolonged presynaptic depolarizations have revealed that release occurs in distinct kinetic phases, with time constants of ~ 2 ms and ~ 20 ms called fast- and slowly releasable pool (FRP and SRP, respectively) [45–47]. Some studies have interpreted the different release kinetics as a result of different distances of vesicles to Ca2+ channels (the" positional" model; [46,48,49]), but other studies concluded that SRP release results from intrinsically slower release [47,50] (see Discussion). To further validate the exclusion zone model and the supercluster model, we modelled release in response to a prolonged presynaptic depolarization to 0 mV, using the large active zone with either the supercluster—or the exclusion zone arrangement of Ca2+ channels (Fig 6). To our initial surprise, both models produced fast release, without obvious SRP component (Fig 6A, 6B, and 6C). Exponential fits of the average cumulative release rates yielded time constants of 1. 52 ms and of 0. 78 ms for the supercluster model and the exclusion zone model, respectively (Fig 6B, red fit lines). In each case, a double-exponential function did not improve the fit (not shown). When we plotted release from the individual sites separately, it became apparent that the spread of release delays and 20–80% release times was larger for the supercluster model (Fig 6C and 6D; left) than for the exclusion zone model (Fig 6C and 6D; right). This is expected, since the heterogeneity of the Ca2+ channel—vesicle distances is larger in the supercluster model than in the exclusion zone arrangement (see above; Fig 5). The faster average release time constant in the exclusion zone model as compared to the supercluster model (0. 8 ms versus 1. 5 ms) corresponds to the lower degree of pves heterogeneity in the exclusion zone arrangement. Indeed, far-away vesicles in the supercluster arrangement have release probabilities near zero during AP stimulation (e. g. vesicles #1 and # 2; see above Fig 5D and 5G), and these vesicles are released with an additional delay of ~ 2 ms upon prolonged depolarization to 0 mV (Fig 6C and 6D). Taken together, simulating the response to prolonged presynaptic depolarizations does not give further clues to distinguish between the two Ca2+ channel arrangements (exclusion zone versus supercluster), since both models were unable to predict SRP release with time constants of ~ 20 ms, as observed experimentally in immature calyx of Held synapses [45,47]. We conclude that to explain SRP release, one would either need to postulate Ca2+ channel—vesicle distances longer than 400 nm. Alternatively, SRP release could be caused by intrinsically slower release kinetics as suggested by previous Ca2+ uncaging experiments [22,47,51]. We next wanted to investigate whether transmitter release in more mature calyx of Held synapses might still be controlled in a domain overlap fashion. Previous work has shown a developmental tightening in the spatial co-localization between Ca2+ channels and docked vesicles, as demonstrated by an increased efficiency of small Ca2+ charges in inducing release [20–22], and a reduced efficiency of EGTA in suppressing release and reduced Ca2+ current—cooperativity values [12]. The question therefore arises whether release in more mature calyx synapses is converted into a single channel control mode or else, whether several Ca2+ channels still control release of a given vesicle, albeit at somewhat shorter distance. To investigate this issue, we made measurements of Ca2+ current release cooperativity at postnatal day (P) P8–P11, and in P15–P16 old mice, the oldest age group amenable to paired recordings under our conditions. We varied the number of open Ca2+ channels in Ca2+" tail" current experiments with voltage steps of different lengths, and measured the corresponding EPSC amplitudes as a proxy of release (Fig 7A). Plotting the EPSC amplitudes versus the integral Ca2+ charge in double logarithmic coordinates then allowed us to measure the Ca2+ current—release cooperativity (slope value), and the Ca2+ charge (QCa) necessary to evoke an EPSC of 2 nA. The slope values were 4. 75 ± 0. 60 (n = 6 pair recordings) for mice at P8–P11, and 3. 10 ± 0. 11 (n = 5 pair recordings) for P15–P16 mice (Fig 7C). Thus, although there was a significant decrease of the slope value with development (p < 0. 05), the slope of ~ 3 was still significantly higher than 1 (p < 0. 001; one-sample two-tailed t-test), the value expected for the extreme case of single-channel vesicle release control [9,13]. In addition, we observed that the Ca2+ charge (QCa) necessary to evoke an interpolated EPSC of 2 nA was shifted to smaller values (Fig 7B and 7C; p < 0. 05), similarly as reported before for the rat calyx of Held [22]. Both findings indicate a tighter coupling in the more mature mice [12,20,21]. Despite the tighter coupling, the Ca2+ current—release cooperativity was still high (~ 3), which indicates that several Ca2+ channels control the release of a given vesicle also in more mature synapses. Finally, we wished to investigate the EGTA sensitivity of transmitter release, and its dependence on the developmental state of the calyx of Held synapse. We used the membrane-permeable EGTA analog, EGTA-AM, which allowed us to investigate synapses from adult mice (P60–P100) not amenable to paired recordings. We first investigated immature synapses (P8–P11), and found a gradual suppression of EPSC amplitudes upon acute application of 200 μM EGTA-AM (Fig 8A1 and 8A2). After applications times of 15–20 minutes, EPSC amplitudes were suppressed to 16 ± 7% of their control value (n = 2). Using lower concentrations of EGTA-AM (100, and 50 μM), the suppression of EPSCs was smaller, which demonstrates a concentration-dependent effect of EGTA-AM with an apparent half-maximal effective concentration of 87 μM (Fig 8B and 8C). These results agree with earlier findings showing that low millimolar concentrations of EGTA, or 200 μM EGTA-AM suppress transmitter release at the immature rat [18,52] and mouse [12] calyx of Held. We then investigated adult synapses from P60–P100 mice using the same experimental approach. Acute application of 200 μM EGTA-AM again strongly suppressed EPSC amplitudes, with 20. 1 ± 5% of the control EPSC amplitude remaining (Fig 8B1 and 8B2). This value was slightly larger than the remaining fraction of the EPSC at the young calyx synapse, but the difference did not reach statistical significance (p = 0. 6). The dose-dependency of the steady-state block of EPSCs by EGTA-AM indicated a half-maximal effective concentration of 92 μM (Fig 8F), similar to the value found at young calyx synapses. These results show that overall, ultrafast transmitter release at the calyx synapses remains sensitive to EGTA-AM, although its efficiency might be slightly smaller in mature mice, in agreement with recent findings [53]. Therefore, the measurements of Ca2+ current release cooperativity (Fig 7) and of the EGTA-AM sensitivity (Fig 8) suggest that despite a developmental tightening of the Ca2+ channel—release coupling (Fig 7; [12,20–22], domain overlap remains relevant for the mature calyx of Held synapse. This also suggests that while some fast-releasing synapses employ single channel release control [54–56], other fast-releasing synapses like the calyx of Held can rely on release control by several Ca2+ channels.
The supercluster model and the exclusion zone model predicted values of Ca2+ current—release cooperativity in the range of 2. 8–3. 3 (Figs 3–5). The maximal value of this parameter predicted by Ca2+ sensor models with 5 Ca2+ ion binding sites [7,8] is ~ 3. 5, as shown by previous simulations with AP-like stimuli [11,42]. Therefore, both models preserve the high intrinsic Ca2+ cooperativity nearly completely when the number of open Ca2+ channels is varied, as expected for domain overlap control of release [9–13]. Since the Ca2+ current—release cooperativity was therefore not a distinguishing feature between the two Ca2+ channel arrangements, our further analysis concentrated on the heterogeneity of pves between vesicles. Our model allowed us to track the individual release probabilities of multiple docked vesicles, an approach which received little attention in previous modeling studies. We assumed that all docked vesicles belonged to the readily-releasable pool consistent with experimental considerations (see discussion in [29]). Furthermore, we assumed that no mechanism exists which would limit release to single fusion events at a given active zone (thus, no" lateral inhibition" mechanism as proposed in refs. [58,59]). For simplicity, we also assumed that all vesicles have the same Ca2+ sensitivity; thus, our model would be unable to describe the biphasic release response upon Ca2+ uncaging [47]. In the exclusion zone model with a distance of 30 nm which yielded a realistic average pves of ~ 0. 1, we observed mostly single or no release events at the small active zone, whereas multiple fusion events (" multivesicular release" ) were rare (Fig 3C, bottom). Therefore, we think that no special depression mechanism needs to be postulated to explain that under normal conditions, there is a low degree of multivesicular release [60]. We found that the two Ca2+ channel arrangements could mainly be distinguished by the degree of heterogeneity of release probability. The large heterogeneity of pves produced by the supercluster model is expected, because placing all Ca2+ channels in a small spot inevitably produces strong gradients of local [Ca2+]i across an active zone, since active zones can extend over a few hundred nanometers. Furthermore, docked vesicles likely create a hindrance for Ca2+ diffusion towards the more distant vesicles (Figs 4A and 5D). Our modelling approach included such a" shading" effect of Ca2+ diffusion by docked vesicles (see Materials and Methods). To further distinguish between the two Ca2+ channel arrangements, more insights into the heterogeneity of pves amongst vesicles within individual active zones must be obtained. Earlier work using optical methods has shown a substantial variability of release probabilities [61], but this heterogeneity likely represents differences across active zones. At the calyx synapse, prolonged presynaptic depolarizations have been used to estimate release probability heterogeneity [45–47]. Some studies concluded that SRP release with time constants of ~ 20 ms is caused by longer distances of SRP vesicles to Ca2+ channels (the" positional" model; ref. [45,46]). We modelled the response of the supercluster- and exclusion zone models to prolonged voltage-clamp depolarizations (Fig 6). As expected, the supercluster arrangement produced a larger degree of heterogeneity in the release times than the exclusion zone model. However, even in the supercluster arrangement at the large active zone, the spread in release delays was in the range of only 2–3 ms, too small to explain SRP release with time constants of ~ 20 ms [45–47]. We conclude that the separation of FRP and SRP release cannot be used to inform about the validity of the supercluster versus the exclusion zone arrangement. The inability of realistic active zone models to predict slow SRP release with time constants ~ 20 ms under positional assumptions is consistent with the earlier view that the slow release phase is caused by an intrinsic mechanism [47]. To distinguish between the supercluster and exclusion zone arrangements, it will also be necessary to obtain more ultrastructural data on the exact localization of Ca2+ channels at active zones of CNS synapses. At hippocampal synapses, freeze-fracture replica labeling has shown several loosely organized, stripe-like clusters of particles detected by an antibody against CaV2. 1 (P/Q-type Ca2+ channel) (see Figs 6–8 in ref. [28]). The hypothetical distributions of Ca2+ channels produced by the exclusion zone (Fig 3A and 5A) could correspond to such stripe-like clusters of CaV2. 1 particles observed in the previous study [28]. A recent study obtained evidence for the localization of P/Q-type Ca2+ channels from freeze-fracture labeling at the calyx of Held synapse. Combined with functional findings on the inhibition of release by EGTA, a release model for the calyx synapse was derived [53]. The study found clusters of Ca2+ channels; a cluster was defined if particles were located within 100 nm distance. Example images show that CaV2. 1 particles can be located over quite long distances (~ 100–300nm; see e. g. Fig 1C and 6H of ref. [53]). Modelling then suggested that docked vesicles must be localized at a distance of 20–30 nm from a somewhat arbitrary Ca2+ channel cluster of ~ 100 nm diameter, to explain the EGTA-sensitivity of release [53]. Our finding of an exclusion zone distance of 30 nm is consistent with this previous distance estimate. Nevertheless, such distance estimates from modelling are subject to parameters which are only incompletely known, including the Ca2+ channel density [41] and the kinetics, concentration and mobility of endogenous Ca2+ buffers (see S2 Text), and should therefore be taken with care. Based on our computational model, and on the placement of Ca2+ channels at active zones shown so far in EM studies [28,53], we find it unlikely that all Ca2+ channels should be localized in a single cluster (" supercluster" ) within an active zone, as predicted by a previous modelling study [11]. Rather, we regard it as more likely that Ca2+ channels are placed either randomly, or else in several sub-clusters in the void spaces in-between docked vesicles, but with a certain minimal distance to docked vesicles. Ion channels are mobile in membranes as has been found for postsynaptic glutamate receptors [62]; however, the mobile properties of presynaptic Ca2+ channels are less well known [63]. In addition, docked vesicles might also undergo lateral movement at synapses. An exclusion zone of several tens of nanometer around docked vesicles might be caused by a hindered mobility of Ca2+ channels close to the docked vesicle, maybe caused by the multitude of vesicle-near proteins which make up the docking complex and the vesicle fusion machinery [64,65]. In addition, specific proteins likely regulate the minimal distance between vesicles and Ca2+ channels. At the neuromuscular synapse, rib-like proteins have been found which likely hold Ca2+ channels at a certain distance from vesicles [66], and filamentous structures playing similar roles might also be present at CNS synapses [25]. One candidate for such a structure at CNS active zones is the filamentous protein Septin-5. When this protein was genetically inactivated, an increased release probability and a tighter spatial coupling resulted [67], which suggested that Septin-5 normally keeps Ca2+ channels at a certain distance from docked vesicles. During developmental maturation, the Ca2+ channel—vesicle coupling distance becomes tighter (Fig 7) [12,20–22], which probably also requires protein re-arrangements to specifically regulate the minimal Ca2+ channel—vesicle distance with a precision of ~ 10 nm or less. Scaffold proteins at the active zone which regulate Ca2+ channel vesicle coupling distance should be investigated in future studies.
Procedures of mouse breeding, handling, and the sacrification of mice before slice preparation were approved by the Veterinary Office of the Canton of Vaud, Switzerland (authorization no. 2063. 3). The sample of the location of docked vesicles at n = 15 active zones is the same as that previously published by [29] (n = 15 active zones from a P11 control mouse). In brief, serial transmission EM images were taken from 10–20 subsequent thin sections (50 nm), which were made following standard fixation and resin embedding procedures. The images were aligned, and pre- and postsynaptic membranes and vesicles were drawn in each section using the TrackEM2 program [68] of the FIJI software [69], and the extent of the PSD was annotated. Only active zones that were completely contained in the series were maintained in the final data set. Each active zone was 3D- reconstructed by the FIJI software, and the nearest distance between vesicle membrane and active zone membrane was measured in 3D, and plotted. Vesicles at 10 nm or less showed a distinct peak in these distance distributions [29], and were regarded as the pool of" docked" vesicles. In order to map the x-y location of docked vesicles, we unfurled the curved surface of the active zone with its attached docked vesicles, using a Matlab program. The shortest distance from the vesicle membrane to the active zone membrane was calculated in the 3D model, and the corresponding point in the flat surface was taken as the remapped vesicle position. We used the Monte Carlo simulator MCell [70] to track the fate of individual Ca2+ ions that entered through various individual Ca2+ channels. The simulator was run on a Blue Gene/Q system, and the MCell based model is available on ModelDB (http: //senselab. med. yale. edu/modeldb/). In general, 2000–4000 repetitions were run to determine an average" release response" for a given simulated stimulus, which is the release rate for all individual vesicles in response to one stimulus; see e. g. Fig 3B, bottom). This simulation was then repeated for each AP width (usually, n = 9 different AP widths), to construct a given dose-response curve of release rates versus Ca2+ influx (e. g. Fig 2E, bottom). The dose-response curve was fitted with a line in double-logarithmic coordinates to obtain the" Ca2+ current—release cooperativity" (see main text for definition). The highest datapoints were usually excluded from the fit range, since they deviated from the steepest slopes, as expected from the beginning saturation of the Ca2+ sensor for vesicle fusion > 10 μM [Ca2+]i. In addition, the lowermost data points at very low release probabilities were often inaccurate despite the number of individual repetitions (2000–4000, see above); in this case, the lowermost data points were not included in the fit range. In order to evaluate the statistical robustness of the slope estimates, the entire simulation was repeated 3–5 times, using different random seeds for the spatial arrangements of Ca2+ channels; these repetitions are referred to as" independent spatial seeds" in the main text. The simulated compartment measured 600 nm by 600 nm and had a height of 1200 nm. Ca2+ was reflected off the walls; the size of the compartment was chosen based on published nearest-neighbor distances of active zones at the immature rat calyx of Held [71]. Models contained Ca2+ channels, endogenous Ca2+ buffers and Ca2+ sensors for release. Docked vesicles were present as spheres with 45 nm diameter inaccessible to Ca2+ ions. Ca2+ channels opened and closed stochastically during the AP according to a two-state model (Fig 2A). Ca2+ channels had a realistic single channel current of 0. 12 pA at 0 mV [31,41]. The standard AP had a half-width of 0. 49 ms; the peak open probability of Ca2+ channels was ~ 0. 6 (Fig 2D). Ca2+ diffusion away from individual channels was modelled stochastically with MCell. Two mobile Ca2+ buffers, ATP and Parvalbumin, and an immobile (fixed) buffer were present in the simulations, with parameters as given in S1 Table. The Ca2+ sensor was a kinetic model of Ca2+ binding and vesicle fusion, with parameters which reflect extensive previous Ca2+ uncaging studies at the calyx of Held synapse. The specific kinetic model was a 5-site model (Fig 2A), with parameters as reported in ref. [37]. See S2 Text for further justification of the model parameters. Transverse brain slices on the level of the medial nucleus of the trapezoid body were prepared from C57Bl6 mice of three different age groups: P8–P11; P15–P16, and P60–P100 (" adult" ). Mice were killed by decapitation, sometimes following a brief isoflurane anesthesia in the oldest age group, in a protocol approved by the Veterinary office of the Canton of Vaud, Switzerland. Transverse brain slices of 200 μm thickness were made with a vibratome (VT1000 or VT1200; Leica Microsystems, Wetzlar, Germany), and stored in a submerged keeping chamber containing a standard bicarbonate buffered solution (for composition, see below), bubbled with 95% O2 and 5% CO2. For paired pre- and postsynaptic recordings (Fig 6), the extracellular solution contained 50 μM D-APV, 100 μM cyclothiazide (CTZ), 1 μM tetrodotoxin (TTX; all from BIOTREND, Wangen, Switzerland) and 10 mM TEA-Cl (Sigma Aldrich/Fluka, Buchs, Switzerland), to suppress unwanted current components and AMPA-receptor desensitization (CTZ). Series resistances (Rs) were in the range of 3–10 MΩ and 5–25 MΩ for post- and presynaptic recordings, and were compensated by up to 80 and 50% for post- and presynaptic recordings respectively, using the patch-clamp amplifier (EPC10/2; HEKA Electronik, Lambrecht/Pfalz, Germany). For postsynaptic recordings, the remaining Rs error was compensated offline. Brief steps to +60 mV of various lengths were applied from a holding potential of -70 mV, resulting in the recruitment of an increasing number of Ca2+ channels giving rise to Ca2+" tail" currents during the repolarization phase (Fig 6A). In fiber stimulation experiments (Fig 8), afferent fibers were stimulated by a custom-made bipolar stimulation electrode placed close to the midline of the slice, using 0. 2 ms long pulses of 1–10 V amplitude from an isolated stimulation unit (ISO-STIM01D-100, NPI electronic, Tamm, Germany). EGTA-AM (Life Technologies, Zug, Switzerland) was applied by bath perfusion, using stock solutions of 200 mM (in DMSO, Sigma), at the indicated final concentrations (50,100 or 200 μM). Stocks of EGTA-AM solutions in DMSO were kept desiccated at–20°C for a maximum of 3 months. The extracellular solution contained (in mM): 125 NaCl, 2. 5 KCl, 25 NaHCO3,1. 25 NaH2PO4,25 glucose, 1 MgCl2,2 CaCl2,0. 4 ascorbic acid, 3 myo-inositol and 2 Na-pyruvate, continuously bubbled with 95% O2/5% CO2 (pH 7. 4). The intracellular (pipette) solution for pre- and postsynaptic recordings contained (in mM): 130 Cs-gluconate, 20 TEA-Cl, 10 HEPES, 5 Na2-phosphocreatine, 4 MgATP, 0. 3 Na2GTP (pH 7. 2 adjusted with CsOH). For pre- or postsynaptic recordings, this solution was supplemented with 0. 1 or 5 mM EGTA, respectively. All chemicals were from Sigma Aldrich/Fluka (Buchs, Switzerland) unless indicated.
Ca2+ channels provide the rise in intracellular Ca2+ concentration necessary to initiate the membrane fusion of transmitter—filled vesicles at synapses. Because Ca2+ diffuses away from Ca2+ channels, the distance between Ca2+ channels and vesicles on the range of tens of nanometers is a crucial determinant of the vesicle fusion probability. However, there is still little experimental evidence on how Ca2+ channels and vesicles co-localize in the nanospace of a single synapse. We show by computational modelling that the channels should be located at some distance to vesicles (~ 30 nm), to allow for release control by several channels, a release mechanism found at many synapses. In realistic synapses with a high density of docked vesicles, this translates into a likely localization of Ca2+ channels at membrane sites not occupied by docked vesicles. Thus, we present a computational model of how Ca2+ channels can be localized in an active zone with several docked vesicles, to enable control of release by several Ca2+ channels.
An Exclusion Zone for Ca2+ Channels around Docked Vesicles Explains Release Control by Multiple Channels at a CNS Synapse