## Abstract

Although the relationship between exocytosis and calcium is fundamental both to synaptic and nonneuronal secretory function, analysis is problematic because of the temporal and spatial properties of calcium, and the fact that vesicle transport, priming, retrieval, and recycling are coupled. By analyzing the kinetics of sea urchin egg secretory vesicle exocytosis in vitro*,* the final steps of exocytosis are resolved. These steps are modeled as a three-state system: activated, committed, and fused, where interstate transitions are given by the probabilities that an active fusion complex commits (α) and that a committed fusion complex results in fusion, *p*. The number of committed complexes per vesicle docking site is Poisson distributed with mean *n*. Experimentally, *p* and *n* increase with increasing calcium, whereas α and the *p**n* ratio remain constant, reducing the kinetic description to only one calcium-dependent, controlling variable, *n*. On average, the calcium dependence of the maximum rate (R_{max}) and the time to reach R_{max} (T_{peak}) are described by the calcium dependence of *n*. Thus, the nonlinear relationship between the free calcium concentration and the rate of exocytosis can be explained solely by the calcium dependence of the distribution of fusion complexes at vesicle docking sites.

## Introduction

The exquisite physiological regulation of local calcium concentration at the active zone of a presynaptic terminal supports the long-debated hypothesis that synaptic plasticity depends heavily on the modulation of the relationship between calcium concentration and transmitter release. This fundamental relationship is observed in temporally diverse secretory systems including neurons, other secretory cells, and eggs. However, a single analytical description for the action of calcium at the synapse that includes the many sites of modulation and control has yet to be developed. Such a description would be equally valuable in synthesizing and testing the molecular basis of synaptic function, since calcium must interact with ensembles of proteins to effect release, and kinetic parameters need to be assigned to the underlying molecular events. Here, we present an analytical description for the final steps of calcium-triggered exocytosis that relies, on average, only upon the calcium dependence of fusion complexes at vesicle docking sites and may be applicable in describing calcium-regulated exocytosis over the entire temporal range observed from eggs to neurons.

To understand the modulation of calcium-regulated release, it is essential to fully understand the origin of the nonlinear relationship between calcium concentration (external, internal, and local) and calcium-triggered exocytosis (amplitude and rate of response). Dodge and Rahamimoff 1967 postulated that the cooperative action of four calcium ions is involved in transmitter release at the frog neuromuscular junction. This observation was based on the calcium concentration in the external bathing solution. Efforts to relate the external, internal, and local calcium concentrations have focused attention on multiple steps in the exocytotic process that involve calcium. Calcium entry through channels complicates the relationship between internal and external calcium concentrations and, ultimately, the interpretation of power-law dependencies (Llinas et al. 1976, Llinas et al. 1981; Charlton et al. 1982; Augustine et al. 1985; Augustine and Charlton 1986; Schneggenburger et al. 1999). Calcium diffusion and the dynamics and density of channels, and their location relative to the fusion site are predicted to influence the exocytotic response (Chad and Eckert 1984; Simon and Llinas 1985; Zucker and Fogelson 1986; Yamada and Zucker 1992; Bertram et al. 1996; Klingauf and Neher 1997). In fact, direct measurements of free calcium concentrations reveal nonlinear relationships in a variety of calcium-triggered systems including neurosecretory cells and neurons support the existence of power-law dependencies (Zucker et al. 1991; Augustine and Neher 1992; Thomas et al. 1993; Heidelberger et al. 1994; Heinemann et al. 1994; Huang and Neher 1996; Tse et al. 1997; Schneggenburger et al. 1999; Bollman et al. 2000; Schneggenburger and Neher 2000). Models of the exocytotic pathway that incorporate multiple binding reactions at one or more key steps describe many of the kinetic features observed in these preparations (Bittner and Holz 1992; Heinemann et al. 1993, Heinemann et al. 1994; Neher and Zucker 1993; Thomas et al. 1993; von Ruden and Neher 1993). The fraction of vesicles available for release and the calcium dependence and dynamics of this pool are important aspects of these models (von Ruden and Neher 1993; Heinemann et al. 1993, Heinemann et al. 1994; Gillis et al. 1996; Rosenmund and Stevens 1996; Murthy et al. 1997; Stevens and Wesseling 1998; Schneggenburger et al. 1999; Wu and Borst 1999). Overall, the nonlinear relationship between calcium concentration and exocytosis has proven to be a consistent finding whether the calcium is measured as an external calcium concentration, a channel-dependent calcium influx, or an intracellular free calcium concentration.

Preparations that isolate selected steps in the exocytotic pathway reduce the complexity of the overall process, and are useful in developing an analytical description of the secretory process. Cortical vesicle exocytosis in the isolated planar cortex of the sea urchin egg, an example of calcium-triggered exocytosis, is an attractive system for directly studying the relationship between calcium and exocytosis. The final fusion steps of exocytosis are preserved in isolation from kinetically confounding processes, such as the temporal and spatial properties of the calcium signal, vesicle transport and priming, and membrane retrieval and recycling. The magnitude of the free calcium concentration is the single controlling variable determining which fraction of vesicles fuse (Blank et al. 1998a,Blank et al. 1998b). Consequently, the kinetics of exocytosis in the sea urchin cortex reflects the relationship between calcium triggering and the final fusion steps of the exocytotic process.

Previously, we investigated the consequences of having multiple fusion complexes at fusion sites and proposed that the calcium-triggered exocytotic behavior observed in this preparation can be explained by the hypothesis that an increase in the free calcium concentration increases the average number of participating fusion complexes (Vogel et al. 1996). For simplicity, fusion complex activation was assumed to be instantaneous, but this assumption led to limitations regarding the calcium dependence of the fusion probability, and a major difference between the observed and predicted initial time courses. In this paper, the analysis assumes that, upon calcium-dependent activation, the following exists: (1) a probability in time that an active fusion complex enters a committed state; and (2) a probability in time that a committed fusion complex enters a fusion-competent state. Either of these probabilities may be a function of the calcium concentration, but neither is assumed to be a function of either time or vesicle location. The probability of mediating a calcium-triggered fusion event (the fusion efficacy) and the relationship between this probability and the number of fusion complexes has been determined here. The fusion efficacy is linearly correlated with the number of fusion complexes. This correlation supports the hypothesis that, in the absence of other kinetic processes, the calcium dependence of multiple fusion complexes at vesicle docking sites is sufficient to explain the nonlinear relationship between the calcium concentration and the rate of exocytosis. Our analysis will show that Poisson-distributed fusion complexes may be operant in neuronal systems under conditions in which large numbers of vesicles in an ensemble of boutons are released from the most readily releasable pool.

## Materials And Methods

All materials and methods related to the study of exocytosis in the sea urchin, *Strongylocentrotus purpuratus* are described in Blank et al. 1998a. Curve fitting was performed using the Levenburg-Marquardt algorithm with Pearson minimization as implemented in Table Curve2D software (Windows v2.0; Jandel Scientific). Fusion curves were fit using the three-parameter kinetic model. Several criteria were established for evaluating goodness of fit upon minimization. These criteria are parameter standard errors, magnitude and distribution of the residuals, and convergence from different initial values in parameter space. The residuals were randomly distributed about zero (zero mean Gaussian) and, typically, showed <1% deviation (Fitted − Observed)/Observed. The mean coefficients of variation for fitting the three parameters *n*, α, and *p* were 0.4 ± 0.1%, 1.2 ± 0.2%, and 3.1 ± 0.5% (*n* = 73, mean ± SEM), respectively. To further test the goodness of fit upon minimization a small sample, resampling analysis was performed. In this procedure, 90% of the data was randomly selected without replacement and the process was repeated 10 times. These 10 datasets were then individually fit for the three parameters, and parameter values and error estimates were calculated from the mean and standard deviation of the 10 values obtained for each parameter. The estimated errors for all parameters were, both from each separate trial and calculated from the 10 trials, smaller than that obtained in the original dataset and strongly argues for a global minimum in the fitting. The results described were obtained from fitting 73 different fusion curves (*n* = 3, 19, 25, 20, and 6 for pCa 4.12, 4.46, 4.62, 4.84, and 4.98; and *n* = 10, 11, and 3 for pCa 4.46, 4.62, and 4.84, for single and double challenge experiments, respectively), a subset of the 111 egg preparations studied previously (Blank et al. 1998a). Experiments involving calcium concentrations below threshold, in which fusion did not occur, were not fit (*n* = 12). Of the 99 fusion curves available for fitting, 26 are not included in this analysis. These fusion curves were not amenable to fitting because the data were too noisy, or corrupted by the transient lifting of the plasma membrane or passage of air bubbles during solution exchange. Derivative curves were calculated from the percent fusion versus time using third order forward differences of data smoothed using a low-pass filter as implemented in SigmaPlot (Windows v2.0; Jandel Scientific) or Savitsky-Golay smoothing (1–5%) as implemented in TableCurve2D software. Absence of a calcium-dependent trend in α and the *p**n* ratio was evaluated using a one-factor analysis of variance (ANOVA) on both the data and a log transformation of the data. Differences between the parameters describing single and double challenge kinetics were evaluated using a two-sample *t* test done at each calcium level; significance was evaluated at the *P* = 0.05 level. Population estimates are presented as the mean ± SEM for *n* experiments. Note, the error (SEM) of the population estimates for the parameters *n*, α, and *p*, represents the observed biological variability and is larger than the fitting errors associated with the determination of parameter values from a single experiment.

## Results

### A Three-parameter Kinetic Analysis of Exocytosis

We have developed a three-parameter kinetic description for calcium-triggered exocytosis by including an explicit transition between a calcium-triggered precursor and a committed state, to address the major incongruity between the observed initial time course of fusion and the initial time course predicted by our previous analysis (Vogel et al. 1996). The active fusion complex was redefined as the calcium-triggered precursor of the committed fusion complex, where the committed fusion complex is defined as the minimal set of components that can cause the fusion of a vesicle with the plasma membrane. Triggered exocytosis of a population of vesicles was treated as a survival process in which the population of vesicles decreases by fusion. Let *t* equal the time in an arbitrary dimensionless time interval that, for example, could be scaled to 1 s, *n* is an integer representing the number of committed fusion complexes per docking site, α is the probability per time interval that an active fusion complex becomes a committed fusion complex, and *p* is the probability per time interval that a committed fusion complex mediates a fusion event between a vesicle and the plasma membrane. The survival of a single vesicle having a time dependent appearance of *n* committed fusion complexes isS*t,n*,α,*t*,*p*=1−*p*^{t}^{nα,t}.

In general, *n*(α, *t*) = *n* · L(α, *t*), where L(α, *t*) allows for the introduction of a lag phase signifying that the appearance of committed fusion complexes does not occur instantaneously. The lag phase is a characteristic feature of exocytotic kinetics. The deterministic function L(α, *t*) = [1 − (1 − α)* ^{t}*] was used because it is the simplest function that described the observed transition. If higher temporal resolution data provides evidence for multiple transitions in the lag phase, then a more detailed function describing these transitions would be necessary. Both the appearance of

*n*committed fusion complexes and the disappearance (fusion) of a vesicle have been formulated as Bernoulli trials (Papoulis 1984). The Poisson-weighted sum of individual vesicle survival curves is:

*n* represents the average number of committed complexes per docking site. Fusion, F*t*,*n*,α,*p*=1−S*t*,*n*,α,*p*. The rate of fusion, R*t,n*,α,*p* is:

As described previously, the extent of fusion is:

In the limit of small α and *p* (≤ 0.1), (1 − α)^{t} and (1 − *p*)^{t} approach exp(−α*t*) and exp(−*pt*), respectively. In this limit, an alternative expression for the Poisson-weighted sum of individual vesicle survival curves that is continuous in time is:

Now, α and *p* represent rate constants for the specified transition. For the same transition probability (for example, 0.1) but different time intervals (millisecond and second) the characteristic time constants for the transitions τ_{α} = 1/α and τ_{P} = 1/*p*, are 10 ms and 10 s, respectively. Using a discrete time parameter, represented by a sequence of independent trials, is analogous to approximating a one-dimensional random walk using interval transition probabilities to describe movement to the left or right. The discrete and continuous expressions for the Poisson-weighted sum of individual vesicle survival curves describe the kinetics with statistically indistinguishable parameters. Both analytical expressions eliminate the need to fit data to exponential sums, a notoriously difficult analysis requiring data of high precision over a wide range (Istratov and Vyvenko 1999).

### Kinetic Analysis of the Calcium Dependence of Exocytosis

We have used multiple calcium challenges to investigate submaximal exocytotic responses of isolated planar cortices of the sea urchin *S. purpuratus* at suboptimal calcium concentrations (Blank et al. 1998a,Blank et al. 1998b; Coorssen et al. 1998). In a sequence of two challenges, one to an intermediate level of calcium followed by a second higher concentration (double challenge), the same extent of fusion as a single challenge at the higher concentration was observed (Blank et al. 1998a). Our analysis describes the time course of cortical vesicle exocytosis triggered by calcium using multiple calcium challenges (Fig. 1). The calcium dependence of the model parameters *p*, *n*, and α is shown in Fig. 2. There is no statistical difference between the parameters describing either single or double challenge kinetics: this analysis describes both responses equally well (Fig. 1). Note, in the present study, *n* is estimated from a fitting parameter and not calculated from the extent of fusion. Over the ranges of calcium concentrations studied, the *p**n* ratio and the parameter α were found to be independent of pCa (−log [Ca^{2}+]_{free}) (Fig. 2). Both the *p**n* ratio and the log transformation of *p**n*, log (*p**n*), are independent of pCa with *p**n* = 0.029 ± 0.002 (mean ± SEM, *n* = 73 experiments; F = 0.514, df = 5, 67, *P* = 0.765, and F = 0.671, df = 5, 67, *P* = 0.647, ANOVA). The invariance of α with [Ca^{2}+]_{free} (α = 0.113 ± 0.013 [mean ± SEM], *n* = 73 experiments; F = 1.722, df = 5, 67, *P* = 0.141, ANOVA; Fig. 2 C) is consistent with either calcium-independent conversion of activated to committed complexes or calcium-dependent conversion saturating at ∼5 μM calcium. In particular, as neither the ANOVA nor the within subset, two-sample *t* tests were significant, the hypothesis that no significant trend was detected in either α or the *p**n* ratio for both single and double challenges is supported. The correlation between *p* and *n* is consistent with a rate constant proportional to *n*. Since *p* represents a bounded quantity in the discrete expression, the invariance of *p**n* with [Ca^{2}+]_{free} (Fig. 2 D) suggests that this is a limiting relationship of the form *p**n* = κ, a characteristic constant. This was tested by fitting all pairs of *p* and *n* with *p* = 1 − exp(−κ*n*), which approaches κ*n* and 1 for small and large values of the exponent, respectively. The value of κ was identical to that obtained by ANOVA of *p**n* with [Ca^{2}+]_{free}. For 0 < *n* < 0.2, >98% of vesicles have zero or one activated complex, suggesting that at low [Ca^{2}+]_{free}, *p* has a limiting value whose distribution is presently unknown. Substituting α = 0.113 and *p* = 0.029 *n* ignores the distribution of α and κ, but reduces the analysis to one parameter that should, on average, describe the kinetics of fusion. This was tested using a reduced kinetic analysis with only one controlling variable (pCa) derived using the relationship between the extent of fusion and *n* and the empirical relationship (log-normal cumulative distribution with midpoint M = 18.2 μM and width W = 0.23), which describes the extent of fusion as a function of [Ca^{2}+]_{free} (Blank et al. 1998a):

The maximum rate and the time to reach the maximum rate were calculated as a function of pCa and compared with measured values (Fig. 3). Note this is not a fit to the average data, but rather it is the predicted behavior following reduction to one controlling variable, the calcium concentration. The maximum rate is a quadratic function of *n*, (R_{max} ∼ *n*^{a}; a = 1.93 ± 0.25); this is exactly the dependence predicted when *p* is proportional to *n*. The relationship between the maximum rate and *n* suggests that *n* represents the local concentration of fusion complexes. Furthermore, these predicted characteristic relationships strongly support the hypothesis that understanding the properties of the local and stochastic concentration of fusion complexes is important for a full kinetic description of calcium-triggered exocytosis. The nonlinear relationship between the maximum rate of exocytosis and the free calcium concentration can be described by the calcium dependence of *n* and two characteristic constants of the system, α and κ.

### Extension of the Analysis to the Millisecond Time Scale

Variations in the model parameters can capture kinetic features observed in both “egglike” and “neuronal-like” calcium-triggered systems (Fig. 4). Calcium concentrations eliciting *n* ∼ 5 will trigger fusion of >99% of all vesicles present. Changing the time constants from the second to millisecond time scale and decreasing *n* results in limited depletion of the available vesicle pool with fusion occurring on a millisecond time scale. Both the extent of vesicle release and the maximum rate of release can be approximated by a [Ca^{2}+]^{4} power-law dependence (Fig. 5A and Fig. B). Calcium-dependent variations in the lag (time to first fusion event) and time to peak rate (Fig. 5 C) are similar to those observed in neurons (Bollman et al. 2000; Schneggenburger and Neher 2000). Poisson-distributed fusion complexes, or reserve capacity (*n*), coupled with variations in the kinetic time constants and distribution of calcium thresholds is therefore sufficient to explain the wide range in rates and extents of fusion observed in other calcium-triggered exocytotic systems.

## Discussion

What factors contribute to the observed nonlinear relationship between exocytosis and calcium concentration? We present an analytical description for the final steps of calcium-triggered exocytosis that relies, on average, only upon the calcium dependence of Poisson-distributed fusion complexes at vesicle docking sites. The calcium dependence of distributed fusion complexes is sufficient to explain the nonlinear relationship observed in sea urchin cortical vesicle exocytosis. This analysis may be applicable in describing experimental features of calcium-triggered exocytosis that are observed in diverse systems covering a wide temporal range.

### Kinetic Properties of Exocytosis in Eggs Are Similar to Those of Other Calcium-triggered Systems

How do the rates of exocytosis in neurons and other secretory cells compare with those observed in the egg? An 80-μm-diam *S*. *purpuratus* egg has a density of ∼0.9 vesicles/μm^{2} based on ∼1.8 × 10^{4} vesicles/egg (Schroeder 1979). An 80-μm-diam planar cortex is expected to have ∼4.5 × 10^{3} vesicles, which is in agreement with measured values (Zimmerberg et al. 1985). Although kinetically unresolved in this study, fusion rates, calculated from the initial slope, plateau at ∼200% fusion/s (∼10^{4} vesicles/s) when challenged with saturating calcium concentrations (>1 mM) that activate all fusion complexes. As a function of [Ca^{2}+]_{free}, the maximum rates calculated from the analysis of the resolved fusion kinetics correspond to ∼20–2,700 vesicles/s. In the urchin cortex, such rates are primarily a consequence of the large number of releasable vesicles rather than the underlying kinetic rate constant (less than ∼10 s^{−1} maximum based on the time to peak rate occurring in 0.1–0.2 s at saturating [Ca^{2}+]_{free}). The releasable pool in other preparations is estimated to be ∼10–1,000 vesicles which, at the highest rates, are released on a ∼1–10-ms time scale (for review see Zucker 1996). Upon fertilization, the intact egg releases ∼10^{4} vesicles at a maximum rate of ∼2,000 vesicles/s (Vogel et al. 1996). Cortical vesicle release rates (events per second) are therefore comparable to rates reported in other calcium-triggered systems such as melanotrophs (Thomas et al. 1993), chromaffin cells (Augustine and Neher 1992; Heinemann et al. 1994; Parsons et al. 1995; Voets et al. 1999), neutrophils (Nusse et al. 1998), and bipolar neurons (von Gersdorff and Matthews 1994; although Heidelberger et al. 1994 report higher rates). However, estimates of the maximal time constant describing the final fusion steps vary over three orders of magnitude depending upon the preparation. For example, in orders of magnitude, the rate limiting time constant in egg is ∼100 ms, in melanotrophs and neutrophils ∼10 ms (Thomas et al. 1993; Nusse et al. 1998), in chromaffin cells ∼1 ms (Heinemann et al. 1994), and in neurons ∼0.1 ms (Heidelberger et al. 1994). Both the transition to the committed state (α) and the coupling between the fusion probability and the number of committed complexes (fusion efficacy [κ]) must be sites of modulation if an underlying conserved mechanism, captured in this analysis, is responsible for the final fusion pathway in such temporally diverse systems.

### Are Poisson-distributed Fusion Complexes with Distributed Calcium Thresholds a Conserved Feature of Calcium-triggered Exocytosis?

One major feature of sea urchin cortical vesicle exocytosis is submaximal secretion at suboptimal stimulation, also observed in a wide variety of permeabilized preparations (for reviews see Knight and Scrutton 1986; Knight and Baker 1987). This has been reported more recently in the high Ca^{2}+ affinity granule population of human neutrophils (Nusse et al. 1998), PC12 cells (Kasai et al. 1996), and mast cells (stimulated with compound 48/80; Hide et al. 1993). Neuronal preparations also have been shown to produce submaximal responses to suboptimal stimulation. However, these results are consistent with an underlying kinetic process that scales with the number of vesicles undergoing release, and this number is a function of the effective calcium concentration. In cultured CA1 hippocampal neurons, reducing the release probability 20-fold did not affect the time course of the averaged AMPA receptor EPSC (Diamond and Jahr 1995). In calyces of Held, the time course of the averaged EPSC at low and physiological calcium concentrations is similar, whereas the amplitude depended strongly upon [Ca^{2}+]_{free} (Borst and Sakman 1996). If synaptic vesicles have reserve capacity then the prediction of our analysis is that the calcium concentration under these conditions did not reach levels sufficient to activate more than one complex per synaptic vesicle docking site. This does not imply that multiple complexes do not exist at synaptic vesicle docking sites, but rather that the reserve capacity was not used under these stimulus conditions. For example, in the urchin, the maximum number of activated complexes (*n*_{Max}) is approximately seven to nine (Vogel et al. 1996; Coorssen et al. 1998), whereas the [Ca^{2}+]_{free} present during fertilization activates only approximately four to five complexes (Vogel et al. 1996). Calcium-triggered systems requiring fast cyclic patterns of release may be designed to work in a regime where only one complex is activated so as to prevent massive depletion of releasable vesicles, and thereby preventing vesicle transport, maturation, and pool refilling from becoming rate limiting. Since the release of synaptic vesicles is part of a cyclic process, comparison with the urchin system requires that the release probability of the most readily releasable pool be considered.

The calcium dependence of this readily releasable pool was determined in calyces of Held where 20% of 700 quanta (vesicles) are released, on average, by a single presynaptic action potential under physiological conditions (Schneggenburger et al. 1999). A releasable pool of ∼20% suggests that when stimulation is low, synaptic transmission operates with *n* < 0.2; vesicles have only one activated complex. We predict that manipulations decreasing the available calcium decrease *n* further and results in a kinetic response that normalizes to a common time course in which the number of participating quanta varies with [Ca^{2}+]_{free}. Do synaptic vesicles utilize multiple fusion complexes with distributed calcium thresholds under conditions of higher calcium concentrations then those observed during a single action potential?

A single presynaptic action potential is capable of depleting the readily releasable pool in calyces of Held, provided the external calcium concentration is high enough (Schneggenburger et al. 1999). In the same preparation, rapid photolysis of caged calcium producing a [Ca^{2}+]_{free} of ∼10 μM released ∼100% of the releasable pool while the release rate continued to increase with free calcium concentrations greater than the minimum concentrations required to release ∼100% of the releasable pool (Bollman et al. 2000; Schneggenburger and Neher 2000). The variability in the sensitivity to calcium and the increased rate with manipulations that increase [Ca^{2}+]_{free} are observations consistent with synaptic vesicles using multiple fusion complexes with distributed calcium thresholds. In addition, a decrease in the release probability of recently recruited releasable vesicles was observed in the same preparation (Wu and Borst 1999). Although this decrease in release probability of recently recruited vesicles may be due to vesicle location relative to the calcium channel, these results are consistent also with a population of synaptic vesicles having a decreased sensitivity to calcium. Calcium-dependent modulation of the release probability may occur through both the distribution of sensitivities and the number of complexes activated. However, the neuronal data are also consistent with other kinetic schemes; multiple calcium binding steps predict similar behavior (Heinemann et al. 1993, Heinemann et al. 1994). Although synaptic transmission may operate in the regime where less than one complex, on average, is activated, the system may have the capacity to operate in the regime where multiple complexes are activated. This may represent another site for synaptic modulation. As *n* increases, corresponding to the appearance of vesicles with two or more activated complexes, additional terms in the kinetics begin to appear and may manifest as poorer fits with a single exponential. The failure to rigorously describe synaptic release kinetics with one scalable process is predicted to increase with increasing concentrations of calcium and may provide evidence for the existence of reserve capacity on synaptic vesicles.

### Does *n*_{Max} Vary in Different Calcium-triggered Systems?

The size of an exocytotic vesicle may limit the maximum average number of complexes *n*_{Max} present at vesicle docking sites. For example, *n*_{Max} on synaptic vesicles may be fractional such that the proportion of vesicles with two or more complexes is small. One consequence of *n*_{Max} = 0.5–0.1 is that only 8% to <1% of the vesicles present will have two or more complexes, and that 61–90% of the vesicles present will be incapable of fusion at any calcium concentration. The pool of inactive vesicles may represent inefficiency in the process of creating the readily releasable/recycling pool and may not be detected using techniques that rely solely on vesicle release. Such a large pool of inactive vesicles is consistent with the observations that, even under strong stimulus conditions, only ∼10% of the total synaptic vesicle population in a terminal participates in neurotransmitter release and recycling (Murthy et al. 1997; Murthy and Stevens 1999). Our analysis remains valid at fractional values of *n*_{Max} because of the hypothesis that the overall population of complexes is described by an underlying distribution of thresholds defined by the derivative of the calcium activity curve (Blank et al. 1998a,Blank et al. 1998b). In this regime, vesicles either fuse under the action of one committed complex, or they are inactive. An alternative analysis, that fixes the number of complexes on every vesicle at one, would require calcium-dependent rate constants and a threshold function defined by the calcium activity curve to be consistent with the existing data. With α and *p* representing the described previously rate constants, and Thresh(pCa) describing the normalized calcium activity curve, an alternative formulation for the fusion process is:

Although this expression is not appropriate for the urchin, where reserve capacity is well established (Vogel et al. 1996; Blank et al. 1998a,Blank et al. 1998b; Coorssen et al. 1998), it may be descriptive of systems where the number of complexes is fixed at one and the population of complexes has an underlying distribution of calcium thresholds. In a sequential, homogeneous model based on a single, linear reaction, calcium is postulated to influence the rate of exocytosis via calcium-dependent rate constants; both the number and binding characteristics of the calcium-dependent transitions are explicitly specified. Although these models can be highly parameterized and may not statistically differentiate between the numbers and characteristics of the calcium-binding site, the data appear to be best described with three or four calcium-binding steps (Heidelberger et al. 1994; Heinemann et al. 1994). In the present analysis, the detailed relationship between calcium binding and the conversion of inactive to active fusion complexes is unspecified. However, even in the absence of a specified calcium activation process, the reduced kinetic analysis describes the same dependence on calcium as that approximated by a [Ca^{2}+]^{4} power function for R_{max} and extent of fusion, respectively.

### Issues in Identifying the Molecular Basis of Calcium-triggered Exocytosis

Understanding the interactions of calcium with ensembles of proteins is a key step in identifying the molecular basis of calcium-triggered exocytosis. The conserved multigene families of SNARE proteins, VAMP, SNAP-25, and syntaxin, have homologues present in the cortical membranes of sea urchin eggs (Avery et al. 1997; Conner et al. 1997; Tahara et al. 1998). These homologues form the same complexes as those observed in other systems (Coorssen et al. 1998; Tahara et al. 1998). Although the exact roles for these and other proteins have yet to be elucidated, the correspondence and conservation among proteins present from yeast, to urchin, to man suggests a conserved pathway for exocytosis; specific steps in this pathway can be directly probed in sea urchin egg preparations (Vogel et al. 1996; Blank et al. 1998a,Blank et al. 1998b; Coorssen et al. 1998; Zimmerberg et al. 1999, Zimmerberg et al. 2000). It may be premature to couple kinetic properties to identified proteins without knowing the relationship between calcium, triggering, and fusion. The kinetic coupling between multiple identified and unknown steps in cyclic processes can result in the incorrect association of a protein or ensemble of proteins with specific kinetic transitions. The stage specific preparations of the sea urchin simplify this difficult problem (Zimmerberg et al. 1999, Zimmerberg et al. 2000).

### Conclusion

The nonlinear relationship between [Ca^{2}+]_{free} and the rate of exocytosis is a characteristic feature of calcium-triggered secretory systems. The final steps of calcium-triggered vesicle exocytosis are described by an analysis using three parameters: (1) the interval probability that an active fusion complex becomes committed (α); (2) the average number of committed complexes (*n*); and (3) the interval probability that a committed fusion complex results in fusion (*p*). This analysis is based on the hypothesis that exocytotic vesicle docking sites have multiple fusion complexes (reserve capacity), and that increasing the free calcium concentration increases the average number of participating fusion complexes. In the sea urchin planar cortex, the average kinetic behavior of submaximal exocytosis at suboptimal free calcium concentrations, including maximum (peak) rate and time to reach maximum rate, is described by the calcium dependence of the participating fusion complexes. The reserve capacity may be limited in other systems when the rapid release of only a small number of vesicles is required. In contrast, the egg requires that essentially every cortical vesicle fuse after fertilization and utilizes reserve capacity to ensure that this level of calcium-triggered exocytosis occurs.

## Acknowledgments

We thank Drs. Jens R. Coorssen and Nevin A. Lambert for stimulating discussions and critically reading our manuscript.

S.S. Vogel has received support from the National Institutes of Health grant No. NS41055.