## Abstract

A method is presented that allows one to estimate transmitter release rates from fluctuations of postsynaptic current records under conditions of stationary or slowly varying release. For experimental applications, we used the calyx of Held, a glutamatergic synapse, in which “residual current,” i.e., current attributable to residual glutamate in the synaptic cleft, is present. For a characterization of synaptic transmission, several postsynaptic parameters, such as the mean amplitude of the miniature postsynaptic current and an apparent single channel conductance, have to be known. These were obtained by evaluating variance and two more higher moments of the current fluctuations. In agreement with Fesce et al. (1986), we found both by simulations and by analyzing experimental records that high-pass filtering of postsynaptic currents renders the estimates remarkably tolerant against nonstationarities. We also found that release rates and postsynaptic parameters can be reliably obtained when release rates are low (∼10 events/msec). Furthermore, during a long-lasting stimulus, the transmitter release at the calyx of Held was found to decay to a low, stationary rate of 10 events/msec after depletion of the “releasable pool” of synaptic vesicles. This stationary release rate is compatible with the expected rate of recruitment of new vesicles to the release-ready pool of vesicles. MiniatureEPSC (mEPSC) size is estimated to be similar to the value of spontaneously occurring mEPSC under this condition.

Quantitative analysis of transmitter release is crucial for the understanding of mechanisms underlying various forms of plasticity, such as paired-pulse facilitation and short-term depression. However, estimation of release rates requires that the size and time course of the synaptic quantum [the miniature EPSC (mEPSC) for the case of excitatory transmission] are known. Furthermore, analysis of postsynaptic currents at glutamatergic synapses is complicated by the fact that the transmitter is cleared only slowly from the synaptic cleft, such that “residual glutamate” accumulates and activates “residual current” (Trussell et al., 1993; Barbour et al., 1994;Mennerick and Zorumski, 1995; Otis et al., 1996a; Kinney et al., 1997; Carter and Regehr, 2000). Neher and Sakaba (2001) adapted the method of deconvolution to allow for such residual current. By combining this method with the analysis of the variance of postsynaptic current records, they showed that the mEPSC amplitude decreases during prolonged stimulation attributable to desensitization of postsynaptic AMPA-type receptors. This confirmed the results of previous studies at various glutamatergic synapses, demonstrating both desensitization (Trussell et al., 1993; Otis et al., 1996b) and saturation (Clements et al., 1992; Jonas et al., 1993; Tang et al., 1994; Liu et al., 1999) of postsynaptic receptors during intense stimulation. Deconvolution was found to be very useful for studying large and rapid changes in transmitter release rate (Sakaba and Neher, 2001). However, the method is not very accurate during low, continuous release, in the presence of residual current. The analysis of EPSCs in such a situation is necessary not only at the calyx of Held but also at other synapses at which asynchronous release plays a major role in transmission (Parsons et al., 1994; Rieke and Schwartz, 1994; von Gersdorff and Matthews, 1994;Lu and Trussell, 2000). Therefore, it is desirable to explore an alternative method that is reliable under such conditions.

Fluctuation analysis has been used extensively in electrophysiology to study the properties of elementary events that underlie cellular signals, either channel currents or synaptic quanta (Katz and Miledi, 1972; Anderson and Stevens, 1973; Neher and Stevens, 1977; Colquhoun and Hawkes, 1981; Heinemann and Sigworth, 1993). However, the application of fluctuation analysis to the study of synaptic processes has been primarily restricted to either the statistical analysis of fluctuations in the amplitudes of evoked responses (for review, see Zucker et al., 1999;McLachlan, 1978) or the study of channel properties (Traynelis and Jaramillo, 1998). A few studies, however, have estimated release rates and mEPSC properties under stationary or slowly varying conditions (Segal et al., 1985;Fesce et al., 1986; Martin and Finger, 1988; Rossi et al., 1994). They demonstrated that such an analysis, based on an extension of Campbell's theorem, can be very powerful in elucidating mEPSC properties, particularly if high-pass-filtered records are used.

Here, we adapt this method to the specific case of a large glutamatergic synapse, by calculating the variance and cumulants of higher order (skew and fourth cumulant) (Courtney, 1978; Segal et al., 1985; Fesce, 1990; Heinemann and Sigworth, 1991) of the postsynaptic current. By combining information from these quantities, we allow for the contribution of the residual current to noise and obtain reliable estimates for the amplitudes and rates of mEPSCs under moderately nonstationary conditions. In particular, we estimate the release rate and mEPSC amplitude under the condition of deep depression (long-lasting strong stimulation) in which releasable vesicle pools should be empty and remaining release should be representing the delivery of new vesicles to the release-ready pool. We find a release rate that is compatible with the time course of recovery of the release-ready pool after stimulation, and we find normal mEPSC size. This latter finding indicates that, even after prolonged stimulation, partially filled vesicles do not contribute appreciably to the EPSC, as might be expected if vesicles recycle rapidly (Pyle et al., 2000). It also indicates that, even under these strong stimulation conditions, no vesicles contribute that are located more remote from postsynaptic densities.

## THEORY

Following Segal et al. (1985) and Fesce (1990), we assume that the postsynaptic current is a sum of randomly occurring mEPSCs, each of which has a mean time course, h·F(t). Here, h is an amplitude factor, and F(t) is the time course of the mEPSC scaled to a peak amplitude of 1. Additionally, current fluctuations around the mean, attributable to the random opening and closing of channels, have to be considered. For simplicity, we assume that these fluctuations are independent of the mean mEPSC time course, have a Gaussian distribution, and are characterized by a variance, V_{c}, that is proportional to the total current. This assumption is problematic for the rising phase of the mEPSC, when channel openings are highly synchronized; however, for any given mEPSC, the variance attributable to the mean time course is large compared with the channel variance around its peak. The assumption is well warranted in the late decay phase of an mEPSC, when channels flicker is uncorrelated, and also for the residual current generated by slowly varying accumulations of neurotransmitter in the synaptic cleft. Below, we will demonstrate by simulation that the assumption is valid for the purpose of this study.

Campbell's theorem (Campbell, 1909; Rice, 1944), which is the basis of this technique, states that mean, variance, and the higher semi-invariants (or cumulants) of the noise signal generated by a stream of randomly occurring elementary events are proportional to the frequency of occurrence of elementary events, ξ, and to the integral over the n th power of the elementary signal waveform:
Equation 1Here λ_{n} is the n th cumulant of the fluctuating signal, ξ is the frequency of occurrence of elementary events, h is the peak amplitude of the elementary signal (mEPSC, see above), and I_{n} is a short-hand notation for the integral over the n th power of the normalized mEPSC time course [F(t)]. Three requirements have to be satisfied for Equation 1 to hold: (1) ξ, the rate of independently and randomly occurring mEPSCs, has to be reasonably constant within the observation interval, (2) F(t) can be different from zero only within an interval short with respect to the observation interval, and (3) the total current is a linear superposition of randomly occurring elementary events (Rice, 1944; Segal et al., 1985).

When these conditions are satisfied, Equation 1, applied to the cases n = 2 (variance) and n = 1 (mean), readily leads to the equation used for standard noise analysis (h = λ_{2}/λ_{1}·(I_{1}/I_{2})). However, in the presence of a residual current, standard noise analysis is not applicable, because the elementary event is very long lasting (including a long “tail” of its contribution to the residual current) and the ratio I_{2}/I_{1}cannot be evaluated because of the ill-defined and variable nature of the residual current. Also, residual current is most likely a nonlinear function of residual glutamate, which violates the assumption of linear superposition. However, a high-pass-filtered version of the original current more readily satisfies the requirement of Campbell's theorem because high-pass filtering turns the elementary events into short, spike-like signals and eliminates the slowly-varying nonlinear components. Additionally, high-pass filtering allows more accurate measurement of the skew and higher moments (Segal et al., 1985). After high-pass filtering, however, the mean of the signal is zero and can no longer be used for calculation of ξ. For the case that the variance originating from channel flickering, V_{c}, is small after high-pass filtering, an estimate for ξ can be obtained from the ratio of skew and variance instead. Unfortunately, we found that, at the calyx of Held, channel variance, V_{c}, very often makes a substantial contribution to the total variance, especially when cyclothiazide (CTZ) is used to block glutamate receptor desensitization (Trussell et al., 1993; Yamada and Tang, 1993). However, V_{c} is expected to be proportional to the mean current (as measured before high-pass filtering). If, as argued above, V_{c} is statistically independent of mEPSC variance λ′_{2}, the total variance after filtering can be written as a sum of the two contributions:
Equation 2where i′ is a proportionality constant representing the filtered single channel current, I′_{2}is the integral over the square of the filtered EPSC time course, and I_{p} is the mean postsynaptic current before high-pass filtering. The primes (′) in this and the following equations denote the fact that they refer to the filtered record.

The skew (λ′_{3}) of the filtered postsynaptic current does not contain contributions from the residual current (the skew of a Gaussian distribution is zero) and is, therefore, given by the following:
Equation 3where I′_{3} is the integral of the filtered mEPSC raised to the third power.

From Equations 2 and 3, we obtain the following:
Equation 4and
Equation 5Because λ′_{2}, λ′_{3}, and I_{p} can readily be measured and both I′_{2} and I′_{3} can be calculated from the mEPSC time course (see legend of Fig. 1), these equations allow the amplitude, h, and the rate, ξ, of mEPSCs to be calculated from postsynaptic current traces. i′, the proportionality constant between channel variance and mean current, has to be determined independently as described in Results. Alternatively, it can be calculated from the cumulants (including information from the fourth cumulant, λ_{4}; see below). The fourth cumulant λ_{4} is given by the following:
Equation 6where μ_{4} is the fourth central moment. Like the other cumulants, λ_{4} has the property that, when derived from a sum of two statistically independent signals, it is equal to the sum of the two cumulants. Furthermore, both the skew and the fourth cumulant of a random signal with Gaussian distribution are zero, such that the channel variance does not contribute to either of them. We therefore obtain, in analogy to Equations 4 and 5, the following:
Equation 7and
Equation 8In practice, however, λ′_{4} can be determined with reasonable accuracy only for long stretches of data; therefore, as a rule, we determine h and ξ from Equations 4 and 5.

In some of our records, however, we encounter suitable stretches of data with relatively large residual current and little mEPSC activity. From these, i′ can be calculated, because the variance is dominated by residual current. mEPSC-induced fluctuations, in this case, constitute a minor contribution, which can be taken into account by combining Equations 2, 3, and 6 for the calculation of i′:
Equation 9It can be seen that i′ is given by the ratio λ′_{2}/ I_{p} (which is the expression for pure channel-derived fluctuations) and a correction term involving the ratio of third and fourth cumulants. We thus can calculate i′ on a suitable segment and determine h and ξ according to Equations 4 and 5 for the rest of the trace unless there is reason to believe that i′ varies along that trace. It will be shown below that h can be calculated at much better resolution than ξ. Thus, it is advantageous to calculate h from Equation 4 and insert this (averaged over longer stretches of data) together with i′ into Equation 2, for calculation of ξ.

An additional complication arises from the fact that, for most synapses, the amplitude, h, is not a constant but rather a stochastic variable with a certain distribution. Fortunately, this amplitude distribution of mEPSCs can be measured at the calyx of Held (Borst and Sakmann, 1996; Chuhma and Ohmori, 1998; Schneggenburger et al., 1999). Assuming that the different size classes of mEPSCs occur statistically independently and that no correlation exists between size and release probability, the quantities h, h^{2}, h^{3}, and h^{4} in Equations 2, 3, and 6 can be replaced by their expectation values (Fesce, 1990), and we arrive at the following:
Equation 10and
Equation 11Here, we have introduced the quantities h′_{s} and ξ′_{s}, which, respectively, represent estimates proportional to h and ξ from the skew of the filtered records. Both h′_{s} and ξ′_{s} have to be multiplied by calibration factors H′_{s} and Z′_{s}, which, according to Equations10 and 11, are given by the following:
Equation 12and
Equation 13Equivalent equations and definitions (H′_{4}, Z′_{4}; see legend of Table 1) can be introduced to replace Equations 7-9 for obtaining estimates from the fourth cumulant and the skew:
Equation 14
Equation 15
Equation 16Calibration factors H′_{s}, Z′_{s}, H′_{4}, and Z′_{4} were calculated from the average time course of the mEPSCs (Fig. 1*B*, Table 1) and from the mEPSC amplitude distribution. The moments of the distribution were as follows:〈*h*〉 = 31.1 pA (mean); 〈*h*^{2}〉 = 1182 pA^{2}; 〈*h*^{3}〉 = 5.4 · 10^{−32}A^{3}; and 〈*h*^{4}〉 = 2.91 · 10^{−42}A^{4}. The fact that the calibration factors contain the quantity of interest, 〈*h*〉, is of little concern, because they do so in combination with the higher moments in a way that they are nondimensional and invariant, when an amplitude distribution uniformly shrinks or expands. Because we are mainly interested in using Equation9 to study possible postsynaptic desensitization, and, assuming that desensitization scales all of the mEPSC amplitudes by a constant factor, we expect H′_{s} to be constant during our recordings.

In practice, the calibration factors H′_{s}and Z′_{s} are readily calculated from the equations above, provided that the time course of an individual mEPSC and its amplitude distribution is known (see legends of Fig. 1 and Table 1). Alternatively, they can be determined from simulated data, using the measured amplitude distribution for the simulation (see below) and applying exactly the same routines, as used for the analysis of experimental postsynaptic currents. In two previous papers (Neher and Sakaba, 2001; Sakaba and Neher, 2001), we used this latter method to obtain relative values of mEPSCs amplitudes. In the present work, we provide calculated values throughout, using mainly Equations 10-16.

## MATERIALS AND METHODS

*Analysis.* Postsynaptic, as well as presynaptic, currents were sampled at rates of 20 kHz after low-pass filtering at either 2.9 or 6 kHz and using a partial hardware series resistance compensation (see below for recording conditions). The remaining series resistance (typically 2–3 MΩ) was compensated by software using a routine similar to that of Traynelis (1998). In this manner, 100% compensation was readily obtained. Analysis was performed with custom-written macros in the environment of IGOR-Pro (WaveMetrics Inc., Lake Oswego, OR). The deconvolution method was performed as described by Neher and Sakaba (2001). The Igor macros used are available on our departmental homepage (www.mpibpc.gwdg.de/abteilungen/140/software).

A first type of fluctuation analysis was performed on single traces of postsynaptic currents by calculating sliding averages of the mean, variance, skew, and fourth cumulant (averaged over a specified analysis window). For the latter three quantities, the data were band-pass filtered before calculating the moments. Variance and skew were set equal to the second and third moment. The fourth cumulant was calculated according to Equation 6 from the fourth moment and the variance after smoothing these quantities by the sliding window. This order of processing is important, because Equation 6 involves a nonlinear operation. Likewise, the quantities h′_{s} and ξ′_{s} (from Eqs. 10 and 11) or the corresponding quantities h′_{4} and ξ′_{4} (estimates from the fourth cumulant and skew) were calculated after smoothing. For this type of analysis, the sliding windows had to be relatively large and the traces had to be reasonably stationary within the length of a window (see Results for an analysis of the effects of nonstationarity).

For most types of experiments, however, we wanted to analyze records with quite pronounced nonstationarities. Typically, we had 5–10 such traces with very similar mean time courses available. It was then possible to improve the analysis by a method similar to “ensemble noise analysis,” introduced by Sigworth (1980) for nonstationary records of voltage-dependent currents. This technique analyzes fluctuations around the mean of an ensemble. The most effective way of removing long-term trends is to subtract subsequent records from one another. The variance of such difference records is twice the variance around the ensemble mean (Heinemann and Sigworth, 1993). For the skew, such simple subtraction is not possible, because the expected skew of a difference between two similar records is zero. However, the difference between an individual record and the sample mean over N such records is quite suitable for analysis. Such a difference record is given by the following:
Equation 17and the expectation value for its skew, λ_{3D,i} (given statistical independence among traces), is as follows:
Equation 18
where λ_{3} is the skew around the true ensemble mean. When N is greater than 5, the multiplier of λ_{3} is larger than 0.5, and, thus, a substantial fraction of the skew is recovered. The corresponding equations for variance (λ_{2D}) and fourth cumulant (λ_{4D}) are as follows:
Equation 19and
Equation 20The advantage of this type of ensemble analysis is that it tolerates substantial nonstationarities of the original records. To efficiently reduce slow trends between records, we subtracted scaled and shifted versions of the sample mean record. Scalers and shift parameters for individual records were determined by minimizing the mean square deviation between a given record and the ensemble mean over a suitable time window. Thereby, trends and transient artifacts are eliminated but results are not modified in any other way, provided that scalers stay within the range of 0.8–1.2 (see below). The time window for determining the scalers should be selected from a region of the mean trace in which the largest nonstationarities occur (or else the nonstationary region of interest), and it should be long enough that it contains a large number of random events in any one trace. The latter requirement derives from the fact that the least-square optimization procedure reduces the number of degrees of freedom of the fluctuations by two.

The combined effect of high-pass filtering of the current records and of subtracting the weighted mean made the analysis quite independent of nonstationarities, as will be documented by simulations below.

*Filtering.* Digital filtering was performed using combinations of smoothing routines supplied by IGOR, as follows. The rationale for the particular combinations will be given in Results. For low-pass filtering, a sliding average using the “Box-Smooth” function with a specified averaging window T_{1}(usually 0.3 msec) was used. To remove the secondary maxima of the resulting sinx/x-type filtering characteristics, the first smoothing was followed by a second smoothing with a window 0.8 × T_{1}, resulting in a filter function with a −3dD point at f_{o} = 0.5/T_{1} (= 1670 Hz for T_{1} = 0.3 msec), which dropped to values lower than −60 db at ∼1.6·f_{o}. Also, secondary maxima were lower than −50 db. For high-pass filtering, a low-pass-filtered record (see above) was subtracted from unfiltered data. Two passes of such high-pass filtering were performed on all records, one using a window T_{h} and a second one using 8·T_{h}. In the first pass, the low-pass-filtered record was shifted by T_{h}/2 before subtraction (to preserve the asymmetry of the waveform), and, in the second pass, it was shifted by −4T_{h}. For band-pass filtering, low-pass filtering and high-pass filtering were applied sequentially. The choice of the smoothing windows T_{1} and T_{h} for low- and high-pass filtering, respectively, will be explained in the context of the simulations described below. The filter characteristics used for most of the variance and skew analyses is shown in Fig. 1. It was calculated with T_{1} = T_{h} = 0.3 msec. Note that a different filter procedure was used byNeher and Sakaba (2001) and Sakaba and Neher (2001).

*Simulation.* The purpose of the simulation was to provide traces for a given set of model assumptions and a specified time course of the release rate r(t). The corresponding IGOR macros were written such that either experimental release time courses [as derived from deconvolution (Neher and Sakaba, 2001)] or appropriate test functions (this paper) could be specified. The macros build the output trace as a sum of miniature currents progressing from beginning to end. For each sample interval, Δt, the number of release events, n_{R}, was determined as the output of a random number generator with a Poisson distribution and the mean number of events equal to r(t)·Δt. If n_{R} was greater than 0, an mEPSC-type waveform (as specified by model parameters) was added to the output trace with its origin at the given sample point. The amplitude of that addition was determined as a sum of n_{R} random numbers drawn from a random number generator with the distribution of the mEPSC amplitudes. In effect, this procedure merged the mEPSCs originating from one sample interval into a single one. This simplification is considered to be valid, if the sample interval (typically 50 μsec) is much shorter than the mEPSC rise time (200–300 μsec). The simulation also tracked the mean residual glutamate concentration, C_{D}, in the synaptic cleft. This was assumed to be a convolution of the released amount equal to the following:
Equation 21where r(t′) is the release rate and c_{r}(t) is a diffusion-type kernel, given by (Crank, 1975):
Equation 22The residual current, I_{D}, which was added to the sum of mEPSC currents, was assumed to be as follows:
Equation 23In a previous publication (Neher and Sakaba, 2001), it has been shown that this formalism allows one to fit glutamatergic postsynaptic currents over a wide range of experimental situations. The choice was explained in that paper; here we just list the parameters and give typical values. The peak amplitude of mEPSC was calculated according to the measured distribution with a mean of 32.1 pA. Its time course was assumed to be a double or triple exponential sum with the following: τ_{r} of 0.2 msec, rise time constant; τ_{1} of 2 msec, decay time constant; α of 0 or 0.2, fraction of a second slow component of decay; and τ_{2}of 10 msec time constant of slow decay phase. These values are typical for recordings in the presence of 100 μm cyclothiazide (which was used in most recordings to prevent glutamate receptor desensitization).

The residual current was described by an amplitude factor, selected such that residual current correctly reproduced experimental currents at long times after release. The parameters were as follows: n = 1.2 exponent of the power law of glutamate channel activation; r_{D} = 0.76 μm diffusional distance; n_{D} = 0.9 exponent of the diffusion law; and D = 30 μm^{2}/sec, diffusion coefficient of the transmitter.

For the analysis of simulated records, the resulting sum of mEPSC-derived and residual currents was filtered and processed in the same way as experimental traces. It should be noted that most of the simulation did not include noise from channel open–close fluctuations and background recording. If required, such “channel noise” was introduced by adding to the simulated trace a real current trace obtained from an experiment in which glutamate-activated current was elicited by bath application of 100 μmS-AMPA.

*Simulation of mEPSCs as a superposition of single channel currents.* For some of the simulations, we constructed simple mEPSCs with mean time courses rising with τ_{0} and decaying exponentially with τ_{1} by activating channels at a rate that has an initial value of A_{0} and decays exponentially with time constant τ_{0}. After initial activation, channels fluctuate between an open state and a short-lived closed state, C_{sl}, according to the following scheme:
The rates and k_{−1} were selected such that bursts of openings were obtained with a mean burst duration of 2 msec, consisting of flickers with a mean duration of 0.5 msec and mean gaps of 0.2 msec. The corresponding parameters are as follows: k_{−1} = 1590 sec^{−1}, β = 3409 sec^{−1}, and α = 2000 sec^{−1} (Colquhoun and Hawkes, 1981). We chose the amplitude of the single channel current such that the mean current during the burst was 1.5 pA, and we chose A_{o} such that the mean mEPSC amplitude was 32 pA. This simple model is not intended to reproduce the complicated and controversial open–close kinetics of AMPA-activated currents (Partin et al., 1996; Rosenmund et al., 1998) but only to estimate the errors made by assuming statistical independence of channel gating fluctuations and release statistics (see Results).

*Physiological recordings.* Presynaptic and postsynaptic recordings at the calyx of Held were performed in the slice preparation of the rat brainstem as has been described previously (Neher and Sakaba, 2001). Briefly, 8- to 10-d-old Wistar rats were decapitated without anaesthesia according to local guidelines. The brainstem was immersed in ice-cold, low-calcium saline, which contained (in mm): 125 NaCl, 2.5 KCl, 0.1 CaCl_{2}, 3 MgCl_{2}, 25 glucose, 1.25 NaHCO_{3}, 0.4 ascorbic acid, 3 myo-inositol, and 2 Na-pyruvate, pH 7.3–7.4 (320 mOsm; bubbled with 95% O_{2} and 5% CO_{2}). Transverse slices of the brainstem (150–200 μm thick) were cut using a vibratome. Slices were incubated in the chamber for at least 30 min at 36°C in normal extracellular solution while being continuously bubbled with 95% O_{2} and 5% CO_{2}. Normal extracellular solution was the same as the low-Ca^{2+}saline, except that 2.0 mm CaCl_{2} and 1.0 mm MgCl_{2} were used and 100 μm CTZ and/or 1 mm kynurenic acid (Kyn) was added, as indicated. Experiments were done within 4 hr after preparation of the slices. All recordings were done at room temperature (∼21–24°C). A presynaptic terminal and a postsynaptic target were simultaneously clamped at −80 mV with patch pipettes. The presynaptic pipette (4–6 MΩ) was filled with a solution containing (in mm): 125–130 Cs-gluconate, 20 tetraethylammonium (TEA)-Cl, 10 HEPES, 5 Na_{2}-phosphocreatine, 4 MgATP, 0.3 GTP, 0.5 EGTA, pH 7.2 with CsOH (310 mOsm). The postsynaptic pipette (2–3.5 MΩ) was filled with the same solution as the presynaptic pipette, except that the concentration of EGTA was increased to 5 mm.

During recordings, 0.5–1 μm TTX, 10 mmTEA-Cl, 0.1 mm 3,4-diaminopyridine, and 50 μmd-AP-5 were added to the normal extracellular solution to isolate the presynaptic calcium current and to block NMDA receptors. This way, AMPA receptor-mediated EPSCs were isolated to monitor quantal release. TTX was purchased from Alomone Labs (Jerusalem, Israel). *S*-AMPA, d-AP-5, and CTZ were from Tocris Cookson (Köln, Germany). Other drugs were from Sigma (Deisenhofen, Germany).

Both presynaptic and postsynaptic cells were whole-cell clamped to a holding potential of −80 mV, using EPC9/2 amplifiers controlled by the Pulse program (Heka Elektronik, Lambrecht, Germany). No liquid junction potential correction was applied. Thirty to 70% of the presynaptic series resistance (R_{s} of 8–30 MΩ, typically 15 MΩ) was compensated. The postsynaptic series resistance (3–8 MΩ, typically 5 MΩ) was compensated so that the uncompensated series resistance was ∼2–3 MΩ. This remaining series resistance was software-compensated off-line, as described above.

## RESULTS

### Simulations

#### Optimizing the filtering procedure for high signal-to-noise ratio and tolerance against nonstationaries

Optimal filtering of postsynaptic currents depends on a number of criteria, such as suppression of background noise (including channel noise), optimization of the ratio between the mean of the estimated cumulants and their intrinsic statistical variation, suppression of trends and nonstationarities in the data, and preservation of asymmetry of the original signal (for skew). The dependence of noise estimates on bandwidth can be derived analytically, and, for many reasons, high-pass filtering with a high cutoff frequency is beneficial (Fesce, 1990). However, at very high frequencies, the noise power originating from the mEPSC time course decreases rapidly, whereas background noise, such as amplifier noise or ion channel noise, increases or stays constant. Therefore, an optimal frequency band for filtering has to be found. In the following, we will first discuss some theoretical expectations and then search for optimal parameters by simulation.

We are interested in determining the relative accuracy of our estimates of the cumulants and want to optimize the filtering procedure so that the coefficient of variation (CV), ς_{λ′v}/λ′_{v}, is as small as possible for a given observation interval. Neher and Sakaba (2001) pointed out that the coefficient of variation of the variance is independent of the frequency of elementary events and is proportional to
, where n is the number of independent samples that can be obtained during the observation interval T. This, in turn, depends on the filtering, given the requirement that the observation for obtaining such a sample has to be longer than the duration of the filtered elementary event. For this reason, the filtering should be done such that the filtered elementary event has the shortest possible half-width. This requirement for short half-width is even more serious for the skew, as pointed out by Segal et al. (1985). The argument calls for high-pass filtering (which eliminates the residual current and reduces the slow decay phase of mEPSCs) with a corner frequency as high as possible. Additionally, high-pass filtering eliminates trends and other nonstationarities in the record. However, when examining noise power spectra of stationary records derived from either a stream of mEPSCs or currents evoked by superfusion of AMPA, it is noticed that channel and instrumentation noise may very well dominate noise spectral density at frequencies above 1 kHz. Also, the spectral density of the mEPSC-derived signal at such high frequencies critically depends on the mEPSC rise time, which, in turn, may be dependent on clamp speed and series resistance of the postsynaptic recording. Therefore, we restricted the band pass on the high end to ∼1.6 kHz, which corresponds to a smoothing window T_{1} of 0.3 msec (see Materials and Methods), or 1.5–2 times the estimated mEPSC rise time (see paragraph on filtering, above). With this restriction, a relatively narrow bandpass (implemented by software, as shown in Fig.1A), appeared to be optimal. We also restricted the filter kernel to as few points as possible to avoid contamination of those parts of the record, which are relatively stationary, by large transients of current in neighboring segments. This, too, is a property of the software filter described in Materials and Methods, as shown in Fig. 1B. The filter converts an mEPSC into a pulse with ∼0.5 msec half-width, preceded and followed by shallow depressions that extend over ∼5 msec. Fig. 1B also shows that the square and the cube of the filtered signal are narrow pulses with hardly any components preceding and following. In the following simulations, we explore the properties of fluctuation analysis of records filtered in such a way.

#### Averaging over 500 msec of a rather wavy record gives useful information on the variance, skew, and fourth cumulant

Postsynaptic current traces were simulated as described in Materials and Methods section. For the first series of simulations, we used the mEPSC parameter set (described in Materials and Methods) without a second decay component (α of 0) and without a residual current (β of 0). These parameters are appropriate for recordings in the presence of 100 μm CTZ.

First, we simulated stretches (500 msec each) of superimposed mEPSCs at constant release rates of 0.5, 1, 2, 8, 12, and 24 events/msec. Example excerpts from such simulated traces before and after filtering are shown in Figure 2, A and B. We calculated the variance, skew, and fourth moment of the filtered records (mean of 0) and formed averages over individual records (500 msec each). We repeated this procedure 50 times and obtained mean values and SDs (given in Table1, top seven rows). Also given in Table 1are the estimates for mEPSC amplitude and release rates as derived from the skew and variance (according to Eqs. 10 and 11) using quantities H′_{s} and Z′_{s}(given by Eqs. 12 and 13). The same estimates based on the fourth cumulant and skew are also shown (Eqs. 14 and 15). Estimates of the mean mEPSC amplitudes (32.1 pA), as well as the release rates, are reproduced quite accurately by the analysis based on skew and variance. In contrast, estimates deviate for rates above 2 events/msec when the analysis is instead based on the fourth cumulant (Table 1, Fig.3A). This is because at higher frequency, according to the central limit theorem, the distribution of current values approaches a Gaussian, which has a relatively small fourth cumulant.

Coefficients of variation, calculated from the means and SDs of Table1, are plotted in Figure 3, B and C. As expected from the consideration above, the estimates for the coefficient of variation of the variance are almost independent of the mEPSC frequency for rates greater than 2 events/msec. The same is also observed for estimates of the fourth moment. However, the accuracy of the fourth cumulant and of estimates derived from it deteriorate at rates greater than 2 events/msec. The skew is significantly less well resolved than the variance, and its resolution decreases for rates greater than 5 events/msec (Segal et al., 1985). Both the variance and (even more so) the fourth moment are less well resolved at very small rates, probably because their distributions are far from Gaussian, even if averaged over long time intervals. The coefficients of variation of rates and amplitudes reflect the properties of the moments from which they are derived (Fig. 3C) (see also the theoretical predictions; Fesce, 1990). The rate estimate derived from the fourth moment is without any significance for rates greater than 2 events/msec. Amplitudes and rates estimated from the skew and variance are optimally resolved at rates between 1 and 2 events/msec and are quite accurate over the entire frequency range studied.

Performing the same analysis on segments of different lengths or averaging values from several similar traces readily demonstrates that the coefficient of variation changes with 1/ T(Neher and Sakaba, 2001), where T is the total recording time over which the cumulants are averaged. Thus, for each estimate, the total recording time necessary for a desired accuracy can be calculated. This recording time is 70 msec for variance when 20% accuracy is desired (with release rates of 2 events/msec) and much longer for the other estimates. For example, an estimation of the mEPSC amplitude with 20% accuracy requires only 110 msec, significantly less than a similar estimate for the release rate (280 msec). This difference is expected, realizing that higher powers of cumulants enter the equation for release rate (Eq. 11) than that for amplitude (Eq. 10).

Second, we wanted to know how sensitive such estimates would be to slow variations in the release rate. To this end, we simulated traces with a mean rate of 2 events/msec that increased and decreased sinusoidally by ±50% with a period of 50 msec. The result for the mean cumulants are shown (Table 1) for comparison with the values for a constant release process. Within the statistical error it is seen that the values agree with those corresponding to a constant rate of 2 events/msec. Thus, we confirm that bandpass filtering effectively eliminates the effects of nonstationarities (Fesce, 1990), if the nonstationarities are sufficiently slower than the low frequency cutoff of the filter. Below, we will present in more detail the effects of more rapid transitions in release rates.

Third, we simulated traces with a constant rate of 2 events/msec and a fourfold-reduced cutoff frequency of the high-pass branch of the filter. Variance and higher cumulants were increased sevenfold to 30-fold (Table 1, compare row two, row nine), whereas the estimates of rate and amplitude were still correct, although significantly more noisy. This analysis shows that high-pass filtering should be done with the highest corner frequency that is compatible with other constraints.

### Contributions by slow decay components, residual current, and channel noise

So far, the simulations were performed with a single component mEPSC, in which fluctuations attributable to opening and closing of channels were not considered. We now proceed to simulations in which we include a slow component of mEPSC decay, as well as channel fluctuations. We demonstrate that these additions lead only to small changes in the mean values, provided that their contributions are slow relative to the bandpass of the filter.

When a slow decay component of 20% amplitude and 10 msec time constant was added in simulations at a release rate of two per millisecond, the mean current at steady state was about twice as large as that for a single component mEPSC. If we included a residual current component, using parameters typical for the calyx of Held (see Materials and Methods) (Neher and Sakaba, 2001), the steady-state currents were increased approximately fivefold. Nevertheless, the variance was unchanged and the skew increased only insignificantly (see Table 1 for simulations at 0.5, 2, 8, and 24 events/msec). Correspondingly, amplitudes and rates estimated from the skew and variance were correct. A similar result was obtained when we added a slow negative component, producing an undershoot.

We did two kinds of simulations to examine the influence of channel noise. In the first simulation, we simply added current traces recorded during S-AMPA bath application (at 100 μm) from a postsynaptic neuron to the simulated traces at 0.5, 2, 8, and 24 events/msec. To do so, we selected experimental traces that had the same mean current as the simulated ones. As expected, the variance of the summed records was increased (Table 1), and the skew and fourth cumulant were changed only insignificantly. More importantly, the SD of all estimates changed very little. The contribution of channel noise to most of the traces was ∼15% of the total variance. When the variance of the AMPA-induced current was measured directly (after bandpass filtering) and plotted against mean current, a straight line was usually obtained. However, the slopes of such lines (corresponding to the parameter i′ in Eqs. 2-10) varied between 12 and 24 fA for different cells (Neher and Sakaba, 2001). When AMPA-induced currents were added to simulated ones, the variance of the superposition was very close to the sum of the individual variances of the AMPA currents and the EPSC simulations. Thus, accurate estimates for amplitudes and rates were obtained when channel variance was accounted for by Equations 10 and 11 (Table 1, last four rows).

In the second simulation, we examined the error that might arise by the assumption of statistical independence between channel noise and the noise originating from the random superposition of mEPSCs. Specifically, we simulated simple mEPSCs explicitly as superpositions of channels (see Materials and Methods), superimposed such mEPSCs randomly at a mean rate of 2 events/msec, and compared the resulting total noise with that of a random superposition of mEPSCs with a fixed time course. Time courses of the latter mEPSCs were similar to that of the mean of the simulated ones. For this simulation, we used a constant mEPSC amplitude, unlike all other simulations. If our assumption about the independence of channel variance, V_{c}, and mEPSCs variance (Eq. 2) were correct, we would expect the difference in variance between the two cases to be given by V_{c}, according to Equation 2. This is exactly what we found. Variance values in the two simulations differed only by ∼20%. This difference (1.9·10^{−23}A^{2}) was very close to the variance of another simulation (1.7·10^{−23} A^{2}), in which the same type of channels was superimposed randomly at a rate that resulted in the same mean current (176 pA). In the control simulation, the total variance (9·10^{−23}A^{2}) was, however, smaller than that of the standard simulation at the same mEPSC rate (Table 1, row three) because, in the latter case, variance is increased by the dispersion of mEPSC amplitudes.

#### Influence of the dispersion of mEPSC amplitudes

The heterogeneity of mEPSC amplitudes has a pronounced influence on the estimates for size and rate. The effect on mean values are taken care of by the correction factors described in Equations 10-13. To investigate whether this dispersion influences the coefficients of variation of the estimates and to test the overall correctness of the analysis (with regard to the expectation values), we compared simulations based on the experimental amplitude distribution (described above) with two additional simulations in which the mean amplitudes were close to the mean of the experimental distribution. In one case, the amplitude was fixed at this value (31.1 pA), and, in the other, it had a very large dispersion, with half of the mEPSCs having an amplitude of 12.5 pA and the other half 52.5 pA. In both cases, amplitudes and rates were estimated accurately and with the same coefficients of variation as in control cases (analyses performed at 2 events/msec).

#### Time resolution of noise estimates

So far, we have shown that high-pass filtering of the records renders noise estimates relatively insensitive to slow variations in release rate. When we analyze experiments in which the release process changes slowly (such as during slowly decaying, asynchronous release after a strong stimulus), the question arises what the “time resolution” of the estimates will be. We already showed above that a total recording time of 70 msec is required to estimate the variance with an accuracy of 20%. Likewise, we can ask what recording time is required to achieve an SD of the variance estimate that is equal to its mean. It is reasonable to define this time as the time resolution of the variance estimate, because this duration is the shortest time interval, T, for which the estimate bears any meaning. Because the SD of the noise estimates at high enough mEPSC rates varies with 1/
T (when the estimates have a Gaussian distribution), we can calculate the time resolution, T_{r}, from values found in Table 1 from the following equation:
Equation 24where CV_{0} is the coefficient of variation (SD/mean), and T_{0} is the time window over which the values in Table 1 were obtained (500 msec). Table2 gives some of the calculated values for the following quantities: variance, AmpS (amplitude, as estimated from skew and variance), RateS (rate from skew, and variance), and the corresponding values for the estimates based on the fourth cumulant and skew (AmpF and RateF). We see that the resolution for variance is as short as 1–3 msec, depending on the mEPSC rate. This demonstrates that an order of magnitude estimate for variance can be obtained within surprising short times. The time resolutions for amplitude and rate estimates are within 10 and 100 msec, with better values at lower mEPSC rates. It should be noted that the values given in Table 2 are the total recording times necessary. Thus, if there is an ensemble of n similar traces available, the time resolution for the average estimate will be approximately T_{r}/n or else the coefficient of variation for the time interval T_{r} will be 1/
n. We performed simulations at 0.5 events/msec over time windows equal to the T values predicted from Equation 24. In approximately half of the cases shown in Table 2, these calculations confirmed Equation 24. However, for estimates involving the fourth cumulant, particularly at large release rates, T_{r} values turned out to be up to two times larger than those extrapolated from the numbers obtained with a 500 msec analysis window (Eq. 24 and Table 1). The reason for the discrepancy is most likely a violation of the assumption (underlying Eq. 24) that the estimates have a Gaussian distribution.

#### Tolerance towards transient changes in release rate

The above consideration addressed only the question of how long a recording interval must be to yield meaningful estimates. In addition, we wanted to find out how the estimates are influenced by trends and transient changes in release rates. To this end, we simulated postsynaptic current traces with abruptly changing release rates, subjected these to the same analysis as used above, and tested how well the resulting rate estimates reproduce such transient changes. We averaged estimates over time windows of 5 msec, which is expected to result in a coefficient of variation of ∼0.6 for the variance trace. The parameters for the mEPSCs were selected to resemble those found in recordings from the calyx of Held, and a residual current was included. In agreement with a previous study, in which variance was analysed only (Neher and Sakaba, 2001), we found that variance and skew are capable of tracking changes in release rate, provided that those are not excessive. In fact, the resolution (both in time and amplitude) is quite remarkable, especially if a few similar traces of a given protocol are available for averaging. Single traces provide order of magnitude estimates for variance and release rates with a time resolution of 5–10 msec. Correspondingly, for averages over 10 similar traces time resolution approaches 1.5–3 msec. As expected from Figure3, the coefficient of variation of the skew was larger than that of variance, and that of the amplitude estimate (as derived from skew and variance) was in between. The estimate of the release rate, on the other hand, was not as well resolved (as expected from Fig. 3 and on the basis of theory; Fesce, 1990). Thus, the best strategy might be to analyze variance and skew together only for mEPSC amplitude. In case that mEPSC amplitude is found to be constant or only slowly varying with time, the rate can be inferred from variance alone by inverting Equation 2 and inserting the known h: Equation 25Averaged simulation results were quite insensitive to variations between individual traces. Thus, results did not change when release rates used for simulations fluctuated by up to 20% between traces. Transients in the variance or skew during abrupt changes in release rate were greatly reduced, when mean subtraction was applied in an ensemble of records (described in Materials and Methods). Figure4 compares simulation results with and without mean subtraction. A series of release episodes (each with a square pulse-like release rate at increasing amplitudes between 5 and 20 events/msec) was simulated, and averages over 10 records are displayed: first without (Fig. 4A) and then after mean subtraction (Fig. 4C). The mean trace (average of 10 records), as well as an example of a difference trace, is shown in Figure 4B. The skew is particularly sensitive to transients. In Figure 4A, it is seen that its estimate is out of range for two of the larger transitions, which is not the case after mean subtraction (Fig. 4C). Variance estimates, on the other hand, were more tolerant to transient changes. An example is given in Figure 5, in which the variance estimate for an abrupt transition from 2 to 100 events/msec and back again is shown. It is free of artifacts, although the rates of individual traces were allowed to fluctuate by up to 20%.

### Experiments on the calyx of Held

#### Estimation of mEPSCs amplitudes using higher moments

We performed experiments on the calyx of Held synapse to determine the range of stimulation conditions over which the amplitude estimate for mEPSCs, calculated as described above, remains constant and agrees with the directly measured mEPSC. From the results of Neher and Sakaba (2001), it is expected that this is the case under standard conditions only for quite mild, short-lasting stimulation. Any stimulus releasing a substantial fraction of the releasable vesicle pool, however, is expected to cause desensitization and a decrease in the mEPSC size. Although Neher and Sakaba (2001) showed that desensitization is largely blocked by CTZ, their analysis of mEPSC amplitudes was not very reliable when release rates were small (<10 msec^{−1}) and stationary in the presence of a large residual current. Therefore, the possibility remains that mEPSC amplitudes were reduced under CTZ during extended stimulation, when release rates dropped toward very low levels and depletion of the release-ready pool of vesicles occurred (Meyer et al., 2001). Then, the size of the release-ready pool of vesicles might in reality be unlimited, which, of course, would imply the breakdown of the vesicle pool concept at the calyx of Held (Schneggenburger et al., 1999; Wu and Borst, 1999; Sakaba and Neher, 2001).

To address these issues, we designed two stimulation protocols. In a first set of experiments the mEPSC amplitude was estimated under mild stimulation conditions, to learn about the applicability of the present method to the calyx synapse (Fig.6A). A second stimulation protocol was designed to estimate the mEPSC amplitude during an episode of steady-state release after a strong bout of exocytosis (Fig.6B). Similarly to Neher and Sakaba (2001), we tested these protocols under three conditions: in the absence of CTZ (control condition), in the presence of CTZ (100 μm), and in the presence of CTZ plus Kyn (1 mm). Unlike the previous study, the present method does not depend on deconvolution but is based entirely on the analysis of fluctuations.

First, we did a set of experiments under control conditions. We used a voltage-clamp protocol for the presynaptic terminal as shown in Figure6A, which is quite similar to the “early noise” protocol of Neher and Sakaba (2001), except that we aimed at lower release rates. The presynaptic terminal was depolarized from −80 to +80 mV. This activated Ca^{2+} channels, but the high positive potential primarily prevented Ca^{2+}influx. Then, the terminal was repolarized to +55 mV for 100 msec. The clamp potential was adjusted so that release rates were low (<5 msec^{−1}), which is suitable for the analysis described here (Fig. 3, Table 2). During this period, a slowly rising EPSC, recorded at the postsynaptic holding potential of −80 mV, was evoked, which was paralleled by an increase in postsynaptic current noise. Finally, the terminal was held at 0 mV for 20 msec to deplete the releasable pool of synaptic vesicles completely (the depleting pulse), and then the terminal was repolarized to −80 mV. Using 5–10 similar traces (obtained at intervals of 10–15 sec), variance (λ′_{2}), skewness (λ′_{3}), and fourth cumulant (λ′_{4}) were calculated (Fig. 6A) after mean subtraction and bandpass filtering, as described in Materials and Methods (Eqs. 18-20). To estimate mEPSC amplitudes (and release rates) from skewness and variance, the AMPA channel variance (i′) must be subtracted from the total variance. Measuring λ′_{2} together with λ′_{3}, λ′_{4}, and mean current I_{p} allows one to solve Equations 2-4 for i′, h, and ξ according to Equation 16. We estimated i′ using the EPSC episode, which was recorded between 100 and 200 msec, after the depleting pulse. At this time, release rates should be small, but the residual current component (attributable to delayed clearance of glutamate) can still be observed. The lack of asynchronous release during this period is confirmed by the small amount of skew and fourth cumulant (Fig. 6A). In our experiments, the apparent channel size i′ was between 20 and 30 fA, which is similar to the estimate from AMPA-induced current noise (Neher and Sakaba, 2001). It did not differ systematically between control and the various conditions of drugs applied. Measurement of i′ is not necessary if the analysis is based on the fourth cumulant and skewness.

In the second protocol, a 6 msec depolarizing pulse to 0 mV was applied presynaptically, leading to a short EPSC-like current of 15 nA amplitude (Fig. 6B). This stimulus is strong enough to deplete 50% of the release-ready pool (Wu and Borst, 1999; Sakaba and Neher, 2001). Subsequently, the terminal was held at 20 mV for 100 msec to slowly deplete the remaining pool. This prolonged depolarization led to a steady postsynaptic inward current of only 300 pA (Fig. 6B), which indicates that postsynaptic receptors were strongly desensitized during this whole period of release (see below). Finally, the terminal was held at 0 mV for 20 msec at which maximal I_{Ca} was elicited, to confirm vesicle pool depletion.

Figure 7A shows EPSCs and mEPSC amplitudes estimated from skewness and variance during episodes of moderate release, calculated according to Equation 10 after averaging skewness and variance over 20 msec segments. The mEPSC amplitude was initially approximately −30 pA and then decreased gradually to −10 pA. This trend was not always seen; in some cell pairs, mEPSC amplitudes were constant throughout episodes of moderate release. The CV of the amplitude estimate of this cell, which is based on 200 msec (10 × 20 msec), should be ∼0.2, according to the simulation data of Figure 3; the experimental results confirm this expectation. Figure 7B shows EPSC and mEPSC amplitudes after strong release in the second type of protocol (Fig. 6B), and the mEPSC amplitude was only 10 pA in this case, suggesting desensitization of postsynaptic AMPA receptors.

Second, we performed similar experiments in the presence of 100 μmCTZ (Figs. 8, 9), an agent known to suppress or retard desensitization in several glutamatergic synapses (Trussell et al., 1993; Yamada and Tang, 1993). In Figure 8, it is seen that substantial residual currents develop when a large fraction of the releasable pool of vesicles is released. Figure 9 demonstrates that mEPSC estimates are constant under CTZ both during the moderate release episode (A) and after strong release (B).

Third, we estimated mEPSC amplitudes in the presence of CTZ and Kyn (1 mm). Because evoked EPSCs are reduced to 50% by 1 mm Kyn, we expected that the underlying mEPSCs should decrease similarly. Because mEPSCs are too small to be resolved directly under this condition (Neher and Sakaba, 2001), this experiment is important for estimating the mEPSC size accurately without any special assumptions. Figure10 shows an example of recordings in such conditions with a protocol similar to the “early noise” protocol of Figure 6A. Variance and skewness (Fig. 10) under Kyn are smaller than those in its absence, consistent with reductions in mEPSC amplitudes. Nevertheless, it was still possible to estimate mEPSC amplitudes from skewness and variance. In the case of Figure 10, the mEPSC amplitude as calculated from Equation 10 was ∼13 pA. Unfortunately, we were not able to estimate mEPSC amplitudes from the fourth moment and skewness, possibly because the signal of the fourth cumulant is too noisy. It is also difficult to estimate mEPSC amplitudes under CTZ and Kyn, when the residual current is large (data not shown). This is possibly the result of an increased noise level caused by the residual current. Thus, we did no additional measurement under this condition (CTZ plus Kyn).

Figure 11 shows a summary of mEPSC amplitudes estimated from higher moments under three conditions (control, CTZ, and CTZ plus Kyn). Figure 11A compares mEPSC estimates from skewness and variance during moderate release. These were 25.9 ± 1.7 pA (control; n = 4 cell pairs), 28.2 ± 3.6 pA (CTZ; n = 6 cell pairs), and 11.8 ± 2.6 pA (CTZ plus Kyn; n = 5 cell pairs). Amplitudes under control and CTZ conditions compared well with those of directly sampled mEPSCs. The similarity of mEPSC amplitudes in the presence and absence of CTZ suggests that postsynaptic receptors did not desensitize significantly during moderate release. Under CTZ plus Kyn, mEPSCs were 42% of the value obtained in the absence of Kyn. This reduction factor is similar to that of the evoked EPSC (∼50%) and to the value (47%) obtained from an analysis using variance and deconvolution-based release rates (Neher and Sakaba, 2001).

In Figure 11B, mEPSC amplitudes estimated during moderate release are compared with those derived from episodes after strong release in the absence of CTZ. In each case, both the estimates from skewness and variance (*3rd/2nd*) and fourth cumulant and skewness (*4th/3rd*) are given. We obtained mEPSC estimates of 25.9 ± 1.7 and 35.2 ± 5.3 pA (n = 4 cell pairs) when using skewness and variance or fourth cumulant and skewness, respectively. However, the CV of the fourth cumulant is expected to be quite large (Fig. 3). Therefore, it was necessary to average fourth cumulants over the whole episode of moderate release (0.07 to 0.15 sec in Figs. 6A and 7A) before calculating an mEPSC amplitude estimate. The amplitude estimates based on skew and variance were smaller, whereas those based on the fourth cumulant were larger than the mean of the mEPSC amplitude distribution. On the basis of the simulations, however, we consider the former as more reliable. During steady-state release after an initial bout of strong release, a mean mEPSC amplitude of 9.1 ± 1.3 pA was estimated from skewness and variance. From the fourth cumulant and skewness, a value of 16.9 ± 1.3 pA was obtained (Fig.11B). Compared with the amplitude estimate during moderate release, the mEPSC amplitude was depressed to values between 35 and 50%.

Figure 11C shows summary results in the presence of CTZ. Six experiments, in which i′ and cumulants were calculated, resulted in mean mEPSC size of 28.2 ± 3.6 pA during moderate release periods. Mean values of 21.6 ± 4.1 pA were estimated after strong release using skewness and variance for estimation. When the same experiments were analyzed for fourth cumulant and skewness, we estimated an mEPSC amplitude of 29.5 ± 4.1 pA for moderate release periods and 36.4 ± 12.2 pA after strong release. The comparison shows that mEPSC amplitudes are similar or reduced by not more than 20% after the strong bout of exocytosis in the presence of CTZ. This implies that desensitization of postsynaptic receptors contributes to a decrease in mEPSC amplitudes after strong transmitter release if CTZ is not present (Meyer et al., 2001;Neher and Sakaba, 2001).

### Estimation of release rates using higher moments

Neher and Sakaba (2001) developed a deconvolution method for calculating transmitter release rates that takes into account the residual current component attributable to delayed clearance of glutamate from the synaptic cleft. By using fluctuation analysis, they showed that release rates can be estimated when desensitization was blocked. However, this method has relatively large errors when release rates are small (<10 events/msec^{−1}), and the residual current component contributes most of the EPSC. In that case, a small error in the estimation of the residual current causes serious overall errors. Furthermore, the method is accurate only when mEPSC amplitude is constant. In contrast, rate estimates using higher moments are reliable when release rates are low (Fig. 3) (Fesce, 1990). In principle, the method is able to simultaneously estimate rate and mEPSC amplitude, although with only limited time resolution. In this section, we therefore estimate release rates using higher moments and compare these with the deconvolution-based rates. For the latter, we assumed an mEPSC amplitude, as given by the mean of its distribution. Experiments were performed in the presence of CTZ, because the constancy of mEPSC amplitude, which is the requirement for deconvolution, can be guaranteed by this drug with an error of ∼20% (Fig. 11) (Neher and Sakaba, 2001). No kynurenic acid was added, however, to preserve large fluctuations.

A first type of experiment was performed to test whether the analysis is reliable under conditions of low release rate (Fig.12A). Release rates (ξ) were estimated during moderate release periods by holding the presynaptic terminal close to +55 mV. The protocol was the same as has been described in Figure 6A. As before, 5–10 similar traces were available form each cell pair. Rate estimates, averaged over 10 msec windows, were well defined, and the time course was similar to that of release rates derived from deconvolution (Fig. 12A). The release rates estimated from skewness and variance usually were slightly smaller than those from the deconvolution (Fig.12B). To estimate release rates using fourth cumulants and skewness, much more time averaging was required to reduce CVs of these estimates (Fig. 3). Therefore, rates were estimated after averaging fourth cumulants and skewness over the whole period of moderate release (70–130 msec) (Fig. 12B). Averaged over six cell pairs, a mean release rate of 1.8 ± 0.3 events/msec was obtained from skewness and variance. This compares with 2.7 ± 0.9 events/msec as obtained from fourth cumulant and skewness and with 2.8 ± 0.2 events/msec using the deconvolution method (Fig. 12B). Estimates by using higher moments were slightly smaller than those using deconvolution. However, considering that deconvolution is unreliable when release rates are small (as mentioned above), these differences may not be significant.

Next, we performed experiments to see whether consistent estimates for the release rate are obtained under conditions of steady-state release (Fig. 13). We did this under CTZ with strong continuous stimulation, including a short episode of maximum I_{Ca}, which completely depleted the release-ready pool of vesicles, using the protocol as has been described in Figures6B and 8B. Release rates under these conditions were low (which was appropriate for estimation using higher moments) despite the strong continuous stimulus, because a large fraction of the release-ready pool had already been released earlier in the record. Low rates of release continued after this depolarization, probably because [Ca^{2+}]_{i} was elevated globally. Release should represent mainly the rapid exocytosis of newly recruited vesicles. Therefore, this measurement also constitutes a test of the concept of the releasable pool and its refilling, as described in the recent work of Wu and Borst (1999),Schneggenburger et al. (1999), and Sakaba and Neher (2001). The expectation from this concept is that release rates toward the end of a strong stimulus should be approximately similar to that predicted from the time course of recovery from synaptic depression.

We analyzed release rates during continuous depolarization (100 msec) to +20 mV and the episodes after (0.1–0.28 sec in Fig.13A,B). As seen from the deconvolution results, release rates were quite low during this period (Fig. 13A,B) and hardly changed when I_{Ca} was switched to maximum, which is consistent with the notion that the releasable pool of vesicles is depleted. Release rates during this period were estimated using skewness and variance, which showed that rates were on the order of 5–10 msec^{−1} (Fig. 13B). The rate estimates from skewness and variance were somewhat lower than those from deconvolution in the case of the experiment shown in Figure13B. However, average values from six cell pairs were similar when both methods were compared (average from 0.1 to 0.25 sec in Fig. 13A,B) (12.9 ± 2.6 events/msec from skewness and variance and 15.2 ± 2.9 events/msec from the deconvolution). Estimates from fourth cumulant and skewness showed large fluctuations and were not meaningful, because release rates during this period were higher than those required for a reliable measurement (Fig. 3). From deconvolution, we estimated a mean size of the releasable pool of 3640 ± 692 vesicles (Sakaba and Neher, 2001). Therefore, a release rate of ∼10–15 events/msec, which reflects the recruitment of synaptic vesicles to the releasable pool, would result in a recovery time constant of 300–400 msec. This is very similar to the recovery time constant of the pool observed by Wu and Borst (1999) after strong synaptic depression.

## DISCUSSION

Postsynaptic currents fluctuate strongly whenever presynaptic [Ca^{2+}]_{i} rises for periods longer than the duration of a single mEPSC. At glutamatergic synapses, such fluctuations depend on four parameters: the presynaptic transmitter release rate, the size of the postsynaptic quantum (i.e., the amplitude of the mEPSC), the amplitude of the residual current (i.e., the current attributable to glutamate persisting in the synaptic cleft), and the amplitude of the single channel response. Single channel fluctuations contribute to the amplitude dispersion of mEPSCs, and they render the residual current noisy. We provide a theory that allows to calculate the four parameters from the postsynaptic current, its variance, skewness, and the fourth cumulant. For this analysis, postsynaptic records are required, which are either stationary or vary slowly. By analyzing simulated currents, we confirm the expectation (Fesce, 1990) that high-pass filtering of the currents makes the analysis quite immune against nonstationarities, such that valid results could be obtained under a surprisingly wide range of stimulation conditions (Figs. 1-5). Using simulated records, we determine the accuracy of such noise estimates, and we show, for instance, that a 20% accurate estimate of variance can be obtained from a record as short as 50 msec (Fig. 3, Table 2). Finally, we apply the method to records from the calyx of Held synapse and show that accurate estimates for mEPSC amplitudes are obtained (Figs. 6-11). For very strong, long-lasting stimulation of the presynaptic terminal, we found steady-state release rates on the order of magnitude of 10 events/msec, well compatible with the rate of recruitment of new vesicles to a completely depleted releasable pool (Fig. 13). Thus, our present studies have two important implications: (1) they provide convenient tools to study low rates of quantal release without complicated assumptions, and (2) they constitute a critical test of the concept of the releasable pool of vesicles at the calyx of Held (Schneggenburger et al., 1999; Wu and Borst, 1999; Sakaba and Neher, 2001).

The method presented here is an extension of previous work on end plate potentials at the neuromuscular junction (Segal et al., 1985; Fesce et al., 1986; Fesce, 1990), allowing one to handle the additional complications attributable to residual current. It has improved time resolution compared with previous studies by taking advantage of optimal filtering to extract relevant information. An essential element in our analysis is high-pass filtering of the postsynaptic currents (Fig. 1) before subjecting them to fluctuation analysis. This has been proposed previously (Segal et al., 1985; Fesce et al., 1986) to be very beneficial for eliminating trends and for improving resolution. We show by simulation that, indeed, the resolution of all of the estimates is far less than optimal when a relatively low high-pass-filtering frequency is used. With the simple filtering tools of our analysis platform (IGOR; WaveMetrics Inc.), a suitable software bandpass filter can be constructed by performing smoothing operations with a smoothing window of 0.3 msec, which is approximately 1.5–2 times longer than the rise time of the mEPSC (see Materials and Methods). The resolution might be improved by choosing a bandpass centered at even higher frequencies. However, then the result is expected to depend critically on clamp speed and on changes in mEPSC rise time, which may result from insufficient clamp.

The method is particularly well suited for low release rates, complementing the recently introduced method of Neher and Sakaba (2001), which combines measurement of variance with deconvolution. The latter method is inaccurate at low release rates (see below); however, it is well suited for measuring rapidly changing release rates under conditions in which the mEPSC amplitude can be expected to be constant. Both of these methods complement the traditional (for review, see McLachlan, 1978;Fesce, 1990; Walmsley, 1993) and more recent forms of variance mean analysis of synaptic transmission (Silver et al., 1998; Clements and Silver, 2000; Oleskevich et al., 2000; Meyer et al., 2001; Scheuss and Neher, 2001). Contrary to these methods, which analyse trial-to-trial fluctuations in highly synchronized (action potential induced) responses, our technique is applicable to slowly changing streams of mEPSCs. It has been difficult to estimate such low release rates at CNS synapses because of the overlapping current components attributable to spillover of released transmitter (Barbour and Häusser, 1997). Whenever such spillover currents can be expected to change slowly with respect to the mEPSC rise time, they should not influence the results of our analysis. Furthermore, the method should not be sensitive to small changes in clamp speed when the low-pass section of the bandpass filter (set by a smoothing time window of 0.3 msec in our study) is selected to be somewhat slower than the slowest mEPSC rise time. Therefore, we think that in cases in which the underlying quantal release is relatively homogeneous (and characterized by a relatively well defined time course and amplitude distribution), it should be possible to apply the method to study EPSCs in other synapses. These might be evoked, for example, by focal stimulation (Stevens and Wang, 1995) or by pharmacological manipulations, such as using sucrose and latrotoxin (Rosenmund and Stevens, 1996;Auger and Marty, 1997). Our method should also be useful for analysing asynchronous release (del Castillo and Katz, 1954; Rahamimoff and Yaari, 1973;Goda and Stevens, 1994; Atluri and Regehr, 1998). Especially, the method will be helpful for nonspiking synapses (photoreceptors, Ashmore and Copenhagen, 1980;Schnapf and Copenhagen, 1982; Rieke and Schwartz, 1994, 1996;DeVries and Schwartz, 1999) (retinal bipolar cells,Mittman et al., 1990; von Gersdorff and Matthews, 1994; Lagnado et al., 1996; Matsui et al., 1998) (retinal amacrine cells, Borges et al., 1995) (hair cells, Furukawa et al., 1978;Yamashita and Ohmori, 1990; Parsons et al., 1994; Moser and Beutner, 2000) and GABAergic terminals (Lu and Trussell, 2000) in which asynchronous release plays an important role in transmission. Furthermore, the method might give some hints on dissecting membrane potential fluctuations observed at the neuronal soma in which various inputs converge, and such fluctuations are important determinants for action potential timing (Mainen and Sejnowski, 1995;Stevens and Zador, 1998).

When applied to the calyx of Held, the method allows one to obtain an answer to important questions of synaptic physiology, such as on the contribution of postsynaptic factors that determine mEPSC amplitudes and on the validity of the vesicle pool concept.

Regarding the first issue, we have shown previously that strong excitation at the calyx of Held synapse leads to desensitization of postsynaptic receptors (Meyer et al., 2001;Neher and Sakaba, 2001). We confirm this with the new approach based simply on the fluctuations of the signal (Figs. 6, 7). The finding is consistent with results from other calyx-type synapses (Trussell et al., 1993; Otis et al., 1996b; Oleskevich et al., 2000). CTZ partially abolished the reduction in the mEPSC size, but possibly mEPSCs are still reduced by ∼20% (Figs. 8, 9), which may be attributable to saturation of postsynaptic receptors or slowly developing desensitization (Partin et al., 1994; Geiger et al., 1995). We conclude that the simple use of the amplitude of EPSCs as an index of transmitter release is definitely not a reliable way of estimating quantal release in the absence of CTZ when strong stimuli are applied. There might be less of a problem, however, if one restricts the analysis to low rates of release (<20 msec^{−1}) because mEPSCs are not reduced under conditions of mild stimulation (Figs. 6, 8). On the other hand, such a restriction prevents a better understanding of synaptic function, given the fact that changes in EPSC may be attributable to change in release probability or a consequence of changes in vesicle pool size, the latter of which can be only estimated by applying strong stimuli that can deplete the whole vesicle pool (Wu and Borst, 1999).

In many types of synapses, one can define a so-called releasable pool of vesicles consisting of those vesicles that can be released within a very short time (10–100 msec) by a strong stimulus, such as a depolarization to maximum I_{Ca} (Wu and Borst, 1999; Sakaba and Neher, 2001), a volley of action potentials (Dobrunz and Stevens, 1997;Schneggenburger et al., 1999), or a hyperosmotic stimulus (Rosenmund and Stevens, 1996). The size of this pool is usually calculated from some measure of the cumulative response divided by the unit response. This approach assumes the unit response to be constant. It has been shown, however, that mEPSC amplitudes decrease during strong stimulation because of receptor desensitization (Trussell et al., 1993; Otis et al., 1996b). This can lead to erroneous estimates of pool size (Sakaba and Neher, 2001). Indeed, no proof has been available so far that the final decay of the postsynaptic current in response to a pool-depleting stimulus is not primarily attributable to desensitization. In that case, the pool concept would be invalidated.Neher and Sakaba (2001) showed, by combining variance analysis with deconvolution, that, under CTZ, the mEPSC size stays relatively constant up to ∼80% depletion of the releasable pool. However, the method used in this previous investigation is not very accurate at low release rates and when a large residual current is present, two conditions that become more and more serious as the pool depletes.

The method that we present here is optimally suited for low release rates, as are measured after complete pool depletion. We applied it to this condition under CTZ and, indeed, found a release rate of ∼10 events/msec (Fig. 13). This rate most likely reflects the refilling rate of synaptic vesicles to the releasable pool, the value being consistent with the pool recovery time course estimated from dual-pulse experiment by Wu and Borst (1999). This study relied on the assumption of constant mEPSCs during the stimulus, which could not be proven by the method used. Furthermore, the time course of recovery from depression may be complicated by contributions of mechanisms other than a depletion–refilling process (Hsu et al., 1996;Dobrunz et al., 1997; Bellingham and Walmsley, 1999). Nevertheless, the agreement between the recovery time constant by Wu and Borst (1999) and the refilling rates demonstrated here favors a simple mechanism of pool depletion and refilling.

We also showed that mEPSC size is nearly constant after full vesicle depletion (Fig. 11). This indicates that vesicles only partially filled with transmitter cannot contribute substantially to the EPSC, which might be expected if vesicles recycle rapidly during strong stimulation (Pyle et al., 2000). In addition, a relatively normal mEPSC size even after strong stimulation suggests that not many synaptic vesicles fuse at locations remote from postsynaptic densities; one might expect this to happen if globally elevated Ca^{2+} starts to trigger release of vesicles in parts of the terminal, which are remote from postsynaptic densities.

In summary, we provided the theoretical basis and demonstrated experimental usefulness of higher moments of current fluctuations for the analysis of transmitter release. We hope this method to be helpful for a detailed study of synaptic function in various synapses.

## Footnotes

This work was supported in part by a grant from the Deutsche Forschungsgemeinschaft (Sonderforschungsbereich 406) to E.N. and by an Alexander von Humbolt fellowship to T.S. We thank C. Pouzat for providing a macro for the detection of mEPSCs and L. de Felice, S. Pyott, R. Schnegenburger, and F. Sigworth for providing valuable advice regarding this manuscript.

Correspondence should be addressed to Erwin Neher, Max Planck Institute for Biophysical Chemistry, Department of Membrane Biophysics, D-37077 Göttingen, Germany. E-mail:eneher{at}gwdg.de.