Estimating the Distribution of Synaptic Reliabilities

Emily P. Huang1 and Charles F. Stevens2

1 Neuroscience Department, University of California San Diego, La Jolla 92093; and 2 Howard Hughes Medical Institute, The Salk Institute, La Jolla, California 92037

    ABSTRACT
Abstract
Introduction
Methods
Results
Discussion
References

Huang, Emily P. and Charles F. Stevens. Estimating the distribution of synaptic reliabilities. J. Neurophysiol. 78: 2870-2880, 1997. Using whole cell recording from CA1 hippocampal pyramidal neurons in slices, we examined the progressive decrease of N-methyl-D-aspartate receptor-mediated synaptic responses in the presence of the open-channel blocker MK-801. Previous studies analyzing this decrease have proposed that hippocampal synapses fall into two distinct classes of release probabilities, whereas studies based on other methods indicate a broad distribution of synaptic reliabilities exists. Here we derive the theoretical relationship between the MK-801-mediated decrease in excitatory postsynaptic current amplitudes and the underlying distribution of synaptic reliabilities. We find that the MK-801 data are consistent with a continuous distribution of synaptic reliabilities, in agreement with studies examining individual synapses. In addition, changes in the MK-801-mediated decrease in response size as a consequence of altering release probability are consistent with this continuous distribution of synaptic reliabilities.

    INTRODUCTION
Abstract
Introduction
Methods
Results
Discussion
References

Synapses release neurotransmitter probabilistically (Katz 1969). Release probability---the probability that at least one exocytotic event will occur when a nerve impulse invades a presynaptic terminal---is regulated by the Ca2+-dependent processes governing synaptic vesicle exocytosis and is generally much less than one in central synapses (Allen and Stevens 1994; Hessler et al. 1993; Jack et al. 1981; Korn et al. 1986; Raastad 1995; Raastad et al. 1992; Rosenmund et al. 1993; Walmsley et al. 1988).

Synaptic release probability (Pr) plays an important role in determining the efficacy, or strength, of a given synapse: the average size of the postsynaptic response is proportional to Pr. Thus the strength of a synapse can be altered by modification of its release probability (Bekkers and Stevens 1990; Bolshakov and Siegelbaum 1995; Magleby 1987; Malinow and Tsien 1990; Manabe et al. 1993; Stevens and Wang 1994; Zucker 1989). Because synaptic strength is so important for the function of neuronal circuits, a significant question is how release probabilities are distributed in central synapses. In other words, we want to know the density function that describes the occurrence frequency of synapses with a given Pr.

Two studies (Hessler et al. 1993; Rosenmund at al. 1993) addressed this question in large populations of hippocampal synapses in culture and in brain slices. These studies analyzed the progressive decrease of N-methyl-D-aspartate (NMDA) receptor-mediated synaptic current amplitudes in the presence of the open-channel NMDA receptor blocker MK-801. MK-801 blocks receptors only when they are opened, and because receptors open only when transmitter is released, the progressive block of receptors in MK-801 provides an indirect measure of the transmitter release probability.

The conclusion from this analysis was that populations of hippocampal synapses are nonuniform in release probability. More specifically, it was proposed that hippocampal synapses fall into two classes with about sixfold difference in release probability; we call these classes unreliable and very unreliable. In this scenario, more than half of the synapses belong to the very unreliable class. This conclusion, however, is in conflict with the results of studies using other, more direct methods to examine Pr of synapses (Allen and Stevens 1994; Murthy et al. 1997); these suggest a wide range of release probabilities exists. Indeed, Rosenmund et al. (1993) note the possibility that more than two classes of synaptic Pr are present, although they do not pursue this analytically.

Murthy et al. (1997) measured the endocytotic uptake of the styryl dye FM-143 at stimulated synaptic terminals in culture. The amount of dye taken up at each synapse is directly proportional to its release probability (a synapse with high Pr releases synaptic vesicles more often in response to stimulation and therefore undergoes more vesicle endocytosis). Murthy et al. thus could calculate the Pr for each of a large number of synapses. When a frequency histogram of synaptic Pr was plotted, the distribution appeared to be a continuous function skewed heavily toward low Pr values; in contrast to the conclusions of the MK-801 studies, this directly observed Pr distribution cannot be fitted as a bimodal function. These results are corroborated by data from minimal stimulation experiments (Allen and Stevens 1994), which indicate synapses have Prs ranging continuously from <0.1 to 1.0.

In the following analysis, we demonstrate that the progressive block of NMDA-mediated synaptic responses in the presence of MK-801 is in fact consistent with a continuous release probability distribution. First, we analyze the MK-801 method and derive the mathematical relationship between the progressive block of synaptic responses and the underlying release probability distribution. Then we show that a continuous distribution fits the data as well as the bimodal (two-class) distribution postulated in previous studies. Finally, we demonstrate that the same continuous distribution fits robustly the progressive block of excitatory postsynaptic current (EPSC) amplitudes in MK-801 under two conditions that alter the synaptic release probability: paired pulse facilitation and changing the external Ca2+ concentration.

    METHODS
Abstract
Introduction
Methods
Results
Discussion
References

Slice preparation and physiology

Transverse hippocampal slices (350-µm thick) were prepared from 2- to 3-wk-old Long-Evans rats using standard procedures (Stevens and Wang 1993). Slices were stored at room temperature (24°C) in saline bubbled with 95% O2-5% CO2 for >= 1.5 h before recording; the saline contained (in mM) 120 NaCl, 3.5 KCl, 1.25 NaH2PO4, 26 NaHCO3, 1.3 MgCl2, and 2.5 CaCl2. After incubation, slices were mounted in a recording chamber at room temperature and superfused with external solution identical to the storage solution described above with the addition of 50 µM picrotoxin (a cut was made between area CA3 and CA1 in the slices to prevent epileptiform activity). Flow rate was 2 ml/min. Whole cell recordings were obtained in CA1 pyramidal cells using standard procedures as described by Stevens and Wang 1993. Recording pipettes were 3-4 MOmega and were filled with internal solution containing (in mM) 130 Cs gluconate, 5 CsCl, 5 NaCl, 10 N-2-hydroxyethylpiperazine-N'-2-ethanesulfonic acid, 0.5 ethylene glycol-bis(beta -aminoethyl ether)-N,N,N',N'-tetraacetic acid, 1 MgCl2, 2 Mg-ATP, and 0.2 Li-guanosine 5'-triphosphate.

MK-801 experiments

Whole cell recordings were obtained as above in standard recording solution containing 5 µM 6,7-dinitroquinoxaline-2,3-dione (DNQX). NMDA-mediated EPSCs were evoked by stimulating Schaffer collateral fibers with a bipolar tungsten electrode while clamping the recorded cell at -40 mV; the stimulus repetition interval was 8 s. Recording solution containing MK-801 was added as specified in the text. DNQX and MK-801 were obtained from Research Biochemicals.

Two-pathway experiments were performed as follows: two bipolar tungsten electrodes were placed in the Schaffer collaterals, each ~200 µm to either side of the recorded cell. The electrodes were stimulated alternately, evoking two independent populations of synapses on the cell. Pathway independence was ensured by confirming there was no paired pulse facilitation when one pathway was stimulated shortly after the other (at a 50-ms interval). In two-pathway paired pulse facilitation (PPF) experiments, the PPF interval was 50 ms; the stimulation patterns (repetition = 8 s) were arranged so that the same number of pulses were delivered to PPF and control pathways in the presence of MK-801, as follows
—‖——‖——‖——‖——‖— Control
—∥—————∥————∥— PPF
In Ca2+ experiments, the external Ca2+ concentration was established by perfusing slices for >= 15 min at the specified [Ca2+] before beginning the experiment.

Minimal stimulation

Minimal stimulation experiments were performed using previously described methods (Dobrunz and Stevens 1997; Raastad et al. 1992; Stevens and Wang 1995). Whole cell recordings were obtained in standard recording solution, and EPSCs were evoked by stimulating the Schaffer collaterals with a tungsten bipolar electrode while clamping the cell at -70 mV. Stimulation rate was 0.1 Hz. Stimulation intensity was lowered until a significant number of failures were seen among the responses and the presence of a putative single synapse was confirmed as described in Dobrunz and Stevens 1997. The criteria for single synapse selection can be summarized as follows: the stimulus intensity is just high enough to elicit responses; small changes in stimulus intensity (± 5%) do not alter failure probability or average response amplitude; there is no variation in epsc shape or latency; and the response amplitude histogram does not change when the synapse is paired-pulse facilitated. Failure rate for the synapse was calculated as the fraction of responses that are failures; success rate (or Pr) was one minus the failure rate.

PPF at minimal stimulation was evoked by stimulating in pairs with interstimulus intervals of 19-25 ms (pairs were presented at a rate of 0.1 Hz). The PPF ratio was calculated as the success rate on the second pulse divided by success rate on the first pulse, whereas initial Pr was just the success rate on the first pulse.

Analysis

Signals were filtered at 2 kHz, digitized at 5 kHz, and stored for analysis. Analysis of epscs was performed using programs written in AXOBASIC; relative EPSC size was measured by integrating the current in a 40-ms window around the peak. For minimal stimulation experiments, responses were inspected visually to determine failures.

All fits of data with mathematical models were carried out in Mathcad by least-square fitting procedures. Differential equations describing the NMDA-receptor kinetic scheme also were solved in Mathcad.

    RESULTS
Abstract
Introduction
Methods
Results
Discussion
References

Evaluation of the MK-801 method

Whole cell recordings were obtained from CA1 pyramidal cells in hippocampal slices, and NMDA receptor-mediated EPSCs were evoked by stimulating Schaffer collateral fibers in the presence of 5 mM DNQX while clamping the cells at -40 mV. After a stable baseline was obtained, stimulation was stopped and the NMDA-receptor open-channel blocker MK-801 was added to the bathing medium for 10 min. When stimulation was resumed in the presence of MK-801, the NMDA receptor-mediated EPSCs decreased progressively in amplitude with each stimulus event (Fig. 1A). The shape and time course of the amplitude decrease in 40 µm MK-801 plotted against stimulus number agrees generally with previous studies (Hessler et al. 1993; Kullmann et al. 1996; Manabe and Nicoll 1994; Rosenmund et al. 1993) and is fairly consistent from experiment to experiment (Fig. 1B). In subsequent discussion, we shall call curves such as the one illustrated in Fig. 1B, blocking functions.


View larger version (19K):
[in this window]
[in a new window]
 
FIG. 1. Decrease of N-methyl-D-aspartate (NMDA) receptor-mediated excitatory postsynaptic current (EPSC) amplitudes in MK-801. A: representative experiment showing EPSC amplitude plotted against stimulation number. After a stable baseline is obtained, stimulation is stopped and 40 µM MK-801 is perfused as indicated; stimulation is resumed after 10 min. EPSC amplitudes are normalized to the first response in the presence of MK-801. Insets: average EPSCs (n = 8 and 5, respectively) at indicated points in the experiment. B: summary graph of the decrease of EPSC amplitudes in the presence of 40 µM MK-801, 2.5 mM [Ca2+], and 1.3 mM [Mg2+]. EPSC amplitudes are plotted against stimulus number and are normalized to the first response. Error bars represent SE (n = 6).

Our goal is to determine the relationship between this EPSC amplitude decrease in MK-801, the blocking function, and the underlying release probability distribution for a population of synapses. MK-801 blocks NMDA receptors only when they are open; the fraction of total NMDA receptors that gets blocked at each stimulus, then, depends on synaptic Pr (the higher the Pr, the more likely a synapse will release and have open receptors to block) as well as the effectiveness of receptor block by MK-801 when a release event does occur. The MK-801 blocking function therefore provides us with a transform of the synaptic Pr distribution, but interpreting the underlying Pr distribution from these experiments presupposes an understanding of the nature of receptor block by MK-801. Thus a proper evaluation of the MK-801 blocking function must involve an analysis of channel block by MK-801 under our particular experimental conditions. Toward this end, we obtained blocking functions for different concentrations of MK-801 (10, 20, and 40 µM) and also examined the shape changes of the individual EPSCs as a function of concentration and/or time.

Scaling of NMDA receptor-mediated EPSCs so that their initial peaks match in the presence and absence of 40 µM MK-801 (Fig. 2A) confirmed that the decay of the EPSC is faster when MK-801 is present; this faster decay represents the block of open channels by MK-801 (Hestrin et al. 1990; Huettner and Bean 1988). EPSCs were fitted with a simple four-state kinetic model (Fig. 2B), in which channels make the indicated transitions between the unbound, closed, open, or blocked states. An impulse of transmitter is assumed. For any given experiment, control EPSCs were first fitted with the blocking rate set to zero; EPSCs in the presence of MK-801 then were fitted using the same rate constants as control and allowing the blocking rate to vary. In general, an "unblocking" rate was unnecessary to fit the EPSCs except at low concentrations of MK-801 (Jahr 1992). When unblocking was used, the unblocking rate was at least an order of magnitude smaller than the blocking rate.


View larger version (15K):
[in this window]
[in a new window]
 
FIG. 2. Kinetic model of EPSCs in the presence and absence of MK-801. A: EPSCs in the presence (ii) and absence (i) of 40 µM MK-801 from a representative experiment. i: average of 8 EPSCs from the baseline before the perfusion of MK-801. ii: average of the 1st 5 EPSCs evoked after the perfusion of MK-801. EPSCs have been scaled so that their initial peaks match to show the faster decay of ii. B: EPSCs in the absence (open circle ) and presence (down-triangle) of MK-801 fitted by a 4-state kinetic model (inset). On the left are EPSCs taken from an experiment using 10 µM MK-801, and on the right are EPSCs taken from an experiment using 40 µM MK-801. open circle , average of 8 EPSCs from the baseline; down-triangle, average of the 1st 5 EPSCs in the presence of MK-801. Inset: kinetic scheme; U, unbound state; C, bound and closed; O, bound and open; B, blocked; koff is unbinding rate; ko, opening rate; kc, closing rate; and kb, blocking rate. Average koff, ko, kc, and kb for 10 µM MK-801 experiments are 3.6, 14, 43.3, and 5.2 s-1, respectively (n = 3); for 40 µM MK-801, they are 3.1, 8.8, 43.4, and 24.5 s-1, respectively (n = 5). C: open-channel block fraction as a function of MK-801 concentration. Error bars represent SE (n = 3 for 10 µM; n = 4 for 20 µM; and n = 5 for 40 µM MK-801). Data are fitted by the blocking curve BF = ([MK]/k1/2)/(1 + ([MK]/k1/2), where BF is block fraction and [MK] is MK-801 concentration. Best-fit k1/2 by a least-squares fitting procedure is 28.6 µM MK-801.

The estimated blocking rate increased as the concentration of MK-801 was varied from 10 to 40 µm, indicating a dose dependence. To determine channel block fraction (the fraction of open channels that get blocked by MK-801 during a single synaptic response), the differential equations for the kinetic scheme (Fig. 2B) were solved with the fitted rate constants to calculate the time dependence of the blocked state; block fraction was calculated as the fraction of opened channels that end up in the blocked state after each synaptic current. This value is plotted against MK-801 concentration in Fig. 2C, and the data are fitted by an effective dose-response curve with k1/2 = 28.63 µM MK-801. We emphasize that this constant k1/2 is not the equilibrium dissociation constant of MK-801 binding to the NMDA receptor. Rather, it is the concentration of MK-801 at which half of the receptors that open during a stimulus get blocked. This value depends on the opening and closing kinetics of the receptor as well as the specific experimental conditions. In all subsequent experiments, a concentration of 40 µM MK-801 is used, with an estimated average channel block fraction of 0.611 ± 0.068 (mean ± SE; n = 5).

One might expect that block fraction, as determined above, would influence the rate at which NMDA receptor channels are blocked with successive stimuli. That is, an increase in the efficacy of MK-801 block would be expected to yield experimental blocking functions that decrease faster as a function of stimulus number. To compare the rate of decrease of two different blocking functions, we plotted the amplitude points of one as a function of the amplitude points of the other for each stimulus number. For instance, when we plotted the amplitudes of one 40 µM MK-801 blocking function against the amplitudes of another for each stimulus number (Fig. 3A), the points fell along a straight line of slope = 1, indicating the two functions have an identical rate of decrease. When we compared a 10 µM and a 40 µM MK-801 experiment in the same way (Fig. 3B), however, the points deviated greatly from the slope = 1 line, highlighting the slower decrease of EPSC amplitudes in 10 µM MK-801. Plots for decreases in 20 µM against those in 40 µM MK-801 fell in between (plots not shown).


View larger version (14K):
[in this window]
[in a new window]
 
FIG. 3. Block fraction affects rate of decrease of EPSC amplitudes in MK-801. A: normalized amplitudes of one 40 µM MK-801 blocking function are plotted against the normalized amplitudes of another for each stimulus number. Data fall on a straight line with slope = 1.0, indicating an identical rate of decrease. B: amplitudes of a 10 µM MK-801 blocking function plotted against the amplitudes of a 40 µM MK-801 blocking function for each stimulus number. Points deviate considerably from a line of slope = 1, indicating a far slower rate of amplitude decrease in 10 µM MK-801. C: representative EPSCs at early and late stimulus numbers in the presence of 40 µM MK-801. i: average of the 1st 5 responses and average of 6 responses beginning at stimulus number 50 in the presence of MK-801. In ii, the EPSCs are scaled so that their initial peaks match. There is no significant change in the decay course of the EPSCs, indicating that block fraction remains unchanged.

As the value of block fraction does clearly affect the blocking function, we must be sure the block fraction does not change during the course of the experiment. This point can be settled by scaling and comparing EPSCs after different number of stimuli in the presence of MK-801. Figure 3C shows typical results from an experiment in 40 µm MK-801; EPSCs at later stimulation numbers scale perfectly to EPSCs at early numbers. Thus block fraction does not change with either time or stimulus number in the course of our experiments.

Decrease of EPSC amplitudes in MK-801 depends on block efficacy and synaptic release probability

Now we can derive the quantitative relationship between synaptic release probability and the decrease of NMDA-mediated EPSC amplitudes in the presence of MK-801, that is, between the synaptic reliability distribution and the blocking function. First, let us outline the assumptions used in our derivation. We assume we repetitively stimulated a fixed number of synapses (with a fixed number of postsynaptic NMDA receptors) the release probabilities of which remained constant throughout the experiment [Allen and Stevens (1994) found this to be the case for the synapses in their sample]. With each stimulus, a fraction of the available postsynaptic NMDA receptors are blocked by MK-801. This fraction (a) depends on three factors: Pr, the release probability of the stimulated synapses; theta , the fraction of receptors participating in the current given a release event; and m, the fraction of these participating receptors that get blocked by MK-801 (block fraction). If Pr varies across synapses (Allen and Stevens 1994; Hessler et al. 1993; Murthy et al. 1997; Rosenmund et al. 1993), a also will depend on rho (Pr), the function that describes the relative occurrence frequency of synapses with a given Pr.

With these initial assumptions, we write and solve a difference equation describing the fraction of synaptic receptors remaining unblocked at nth stimulus, Fn (see APPENDIX for full derivation). Our experimental blocking functions measure synaptic current, however, not receptor fraction, so we then calculate from Fn the predicted synaptic current on the nth stimulus, Sn. When this is done (see APPENDIX), we get
<IT>S<SUB>n</SUB></IT>= <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB>ρ(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP> (1)
Sn is normalized so that S0 = 1

From Eq. 1, we see that the predicted decrease of synaptic amplitudes in the presence of MK-801 depends on mtheta , the fraction of NMDA receptors blocked by MK-801 when a synapse releases transmitter, as well as rho (Pr), the Pr weighting distribution. We now are prepared to determine what release probability distributions fit our MK-801 data.

Decrease of EPSC amplitudes in MK-801 is consistent with a continuous release probability distribution

Previous studies analyzing the decrease of EPSC amplitudes in MK-801 (Hessler et al. 1993; Rosenmund et al. 1993) have suggested that synapses are distributed between two classes of release probabilities (termed high and low). As noted earlier, the conclusion that there is a bimodal distribution of release probabilities, however, contradicts more direct estimates of the Pr weighting function from studies using other methods (Allen and Stevens 1994; Murthy et al. 1997), in which release probabilities seem to follow a continuous distribution.

To resolve this conflict, we next examined whether the MK-801 blocking function could be explained as well by a continuous distribution of release probabilities as by the two-class (bimodal) Pr distribution. An ideal bimodal weighting function, as used by Hessler et al. and Rosenmund et al., is represented mathematically as
ρ<SUB><IT>b</IT></SUB>(<IT>P</IT><SUB><IT>r</IT></SUB>) = <IT>A</IT>δ(<IT>P</IT><SUB><IT>r</IT></SUB><IT> − P</IT><SUB>1</SUB>) + (1 − <IT>A</IT>)δ(<IT>P</IT><SUB><IT>r</IT></SUB><IT> − P</IT><SUB>2</SUB>)
where P1 and P2 are the Pr values of the two synaptic classes, respectively, and A is the proportion of synapses with Pr = P1. On the other hand, a continuous distribution can be characterized as a function with a negative exponential dependence on Pr (see APPENDIX), such as
ρ<SUB>c</SUB>(<IT>P</IT><SUB><IT>r</IT></SUB>) ∼ <IT>e</IT><SUP>λ<IT>P<SUB>r</SUB></IT></SUP>
where lambda  is a constant; 1/lambda is the "characteristic Pr" of the distribution.

When Eq. 1 is evaluated with the above bimodal and continuous distributions, we get two predictions for Sn, the decrease in EPSC amplitudes in the presence of MK-801
<IT>S</IT><SUP><IT>b</IT></SUP><SUB><IT>n</IT></SUB> = <IT>Ae</IT><SUP>−<IT>nb<SUB>1</SUB></IT></SUP> + (1 − <IT>A</IT>)<IT>e</IT><SUP>−<IT>nb<SUB>2</SUB></IT></SUP>
where b1 = P1mtheta , b2 = P2mtheta , 0 <=  A <=  1, and
<IT>S</IT><SUP>c</SUP><SUB><IT>n</IT></SUB> = 1/(1 + <IT>n</IT>/<IT>r</IT>)
where r = lambda /mtheta . We see that Sn for the bimodal distribution is the sum of two exponentials (the form used by Hessler et al. 1993; Rosenmund et al. 1993), whereas Sn for the continuous distribution is the sum of an infinite number of exponentials. Note that the parameter r is inversely proportional to the characteristic release probability 1/lambda .

Experimentally obtained blocking functions in the presence of 10, 20, and 40 µM MK-801 were fitted with the predicted Sn functions incorporating the bimodal and continuous distributions by a least-squares-fitting procedure (Fig. 4). For all cells (n = 15), the amplitude decreases were fitted by both models nearly indistinguishably; for 40 µm MK-801 (n = 6), mean values of A, b1, b2, and r were 0.677 ± 0.012, 0.116 ± 0.007, 0.014 ± 0.0006, and 9.78 ± 0.502, respectively. We conclude that the data from MK-801 experiments are perfectly consistent with at least one type of continuous distribution of synaptic release probabilities as well as with a bimodal Pr distribution.


View larger version (18K):
[in this window]
[in a new window]
 
FIG. 4. MK-801 data are fitted by a continuous release probability distribution. Representative blocking function in 40 µM MK-801 fitted by Sn incorporating a bimodal distribution (- - -) and a continuous distribution (------). To reduce variability and for the purposes of quantitative fitting, the EPSC amplitudes are normalized to the y intercept of a line fitted to the 1st 10 points in the decrease. ···, fit of the data by a Sn incorporating a distribution derived from minimal stimulation experiments, as described in the DISCUSSION.

To obtain a quantitative description of the above Pr distributions, we must calculate P1 (=b1/mtheta ) and P2 (=b2/mtheta ), the Pr values of the two synaptic classes in the bimodal distribution, and 1/lambda (lambda  = rmtheta ), the characteristic Pr for the continuous distribution. Calculation of these values requires an estimate of mtheta , the fraction of receptors blocked per stimulus. m is the fraction of participating receptors on a given trial that get blocked and was estimated earlier as 0.611 for 40 µM MK-801. theta is the fraction of receptors that participate in the current when transmitter is released; we take theta  at the previously estimated value 0.5 (Rosenmund et al. 1995) and at 1.0. When theta  is 0.5, P1, P2, and lambda  are 0.38, 0.05, and 2.99, respectively; when theta  is 1, P1, P2, and lambda  are 0.19, 0.025, and 5.98, respectively. Thus for the bimodal distribution, the "high" synaptic Pr may range from 0.19 to 0.38 and the low Pr from 0.025 to 0.05. For the continuous distribution, the characteristic release probability 1/lambda ranges between 0.17 and 0.33.

PPF of a continuous Pr distribution

PPF is an increase in presynaptic release probability that occurs when a synapse is stimulated to release within a few hundred milliseconds of a previous stimulation event; the increased Pr is thought to reflect residual Ca2+ in the axon terminal (see Zucker 1989). To examine the effect of PPF on the blocking function, we performed experiments in which two separate and independent pathways onto the same recorded CA1 cell were stimulated alternately. This was accomplished by stimulating the Schaffer collaterals on either side of the cell. When such pathways were alternately given single stimuli in the presence of MK-801, the EPSC amplitudes of one pathway could be plotted against the amplitudes of the other (for each stimulus number) as a straight line (Fig. 5A), with a slope of 1. 


View larger version (14K):
[in this window]
[in a new window]
 
FIG. 5. Paired pulse facilitation (PPF) changes rate of EPSC amplitude decrease. A: normalized EPSC amplitudes of 1 pathway is plotted against the normalized amplitudes of the second pathway for each stimulus number in a 2-pathway experiment. Points fall on a straight axis with slope = 1.0. B: normalized amplitudes of the first pulse of every pair in the paired pulse facilitated pathway plotted against the normalized amplitudes of every other pulse in the control pathway, for a 2-pathway PPF experiment. The points fall under the straight axis, indicating that the PPF pathway decreases more quickly than control in MK-801.

When one of the pathways was paired pulse facilitated, however, a plot of the first pulse of each pair in the PPF pathway against every other pulse in the control pathway (with stimulus number as the plot parameter) deviated considerably from a line of slope equal to one (Fig. 5B). This deviation reflects a more rapid decrease of the blocking function for the paired pulse pathway, presumably due to an increase in release probability. To see whether we could relate directly the amount of PPF enhancement to the change in amplitude decrease, we modified our continuous distribution model (see APPENDIX) to accommodate an increase in Pr for every second pulse. This increase in Pr on the second pulse is given by the function phi (Pr), defined as the amount of PPF (or PPF ratio) as a function of initial release probability. The result of this simple modification is
<IT>S<SUB>n</SUB></IT>= <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB><IT>e</IT><SUP>−λ<IT>P<SUB>r</SUB></IT></SUP><IT>e</IT><SUP>−<IT>n</IT>(1+φ(<IT>P<SUB>r</SUB></IT>))<IT>P<SUB>r</SUB>m</IT>θ</SUP> (2)
All the parameters of this equation, except for phi (Pr), are fixed by values derived from fits of data from the control pathway (lambda , m, theta ). The problem, of course, is to determine phi (Pr) because the MK-801 experiments do not provide a direct measure of the changes in release probabilities due to PPF. To estimate phi (Pr), we turned to the technique of minimal stimulation, which provides a means of activating putative single synapses on a neuron recorded in slices (Dobrunz and Stevens 1997; Raastad 1995; Stevens and Wang 1995). Synapses stimulated in this manner exhibit frequent response failures that are due to the probabilistic nature of transmitter release (Allen and Stevens 1994).

In minimal stimulation, the Pr of a given synapse is calculated as one minus the failure rate of responses to the minimal stimulus (Allen and Stevens 1994; Raastad et al. 1992; Stevens and Wang 1994). When closely spaced pairs of the minimal stimulus are given, the Pr of the second pulse is increased relative to the Pr of the first pulse, as expected with PPF. Thus we were able to measure PPF ratio and initial Pr for a number of synapses (n = 25, see METHODS). In these experiments, we found synapses with Prs ranging continuously from 0.1 to 1.0; there was no evidence of a bimodal distribution of synaptic Pr (see Fig. 6).


View larger version (9K):
[in this window]
[in a new window]
 
FIG. 6. PPF as a function of initial release probability. PPF ratio, calculated by dividing release probability of the 2nd pulse by the release probability of the 1st pulse, plotted as a function of the release probability of the 1st pulse. Data are from 25 minimal stimulation experiments with paired pulse intervals from 19 to 25 ms.

When the PPF ratio is plotted as a function of initial Pr (Fig. 6), we see a decrease in PPF as Pr goes from zero to one. We thus can approximate phi (Pr) as a straight line
φ(<IT>P<SUB>r</SUB></IT>) = PPF<SUB>max</SUB>− (PPF<SUB>max</SUB>− 1)<IT>P</IT><SUB><IT>r</IT></SUB> (3)
PPFmax is the PPF ratio at very low Pr values; phi (Pr) goes to one at Pr = 1.0. We assessed PPFmax in our MK-801 experiments by calculating average PPF ratio at late stimulus numbers; then, for ease of calculation, we further approximated phi (Pr) as a staircase function decreasing regularly from PPFmax to one in 10 steps. When this phi (Pr) was incorporated into Eq. 2, we obtained an expression describing Sn for PPF synapses that had no free parameters. For four experiments, this expression fitted the MK-801 mediated decrease of EPSC amplitudes in the PPF pathway very well (Fig. 7, A and B). The model incorporating the continuous distribution rho (Pr) e-lambda Pr is therefore robust when the release probabilities are modified, such as with PPF.


View larger version (13K):
[in this window]
[in a new window]
 
FIG. 7. Decrease of EPSC amplitudes in MK-801 for PPF synapses is fitted by a continuous Pr distribution. A: blocking function of a representative control pathway fitted by a continuous distribution model (Scn). Amplitudes are plotted against every other stimulus number; fitted r = 9.92. B: blocking function of the PPF pathway from the same experiment as A is fitted by a Sn incorporating the same Pr distribution as in A. Amplitudes of the 1st pulse of each stimulus pair are plotted against the stimulus pair number. Estimated phi (Pr) = 1.86-0.86 Pr was used to calculate the fitting function Sn.

We also calculated the Sn for PPF synapses assuming a bimodal Pr distribution
<IT>S<SUB>n</SUB></IT>= <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB>[<IT>A</IT>δ(<IT>P</IT><SUB><IT>r</IT></SUB><IT>− P</IT><SUB>1</SUB>) + (1 − <IT>A</IT>)δ(<IT>P</IT><SUB><IT>r</IT></SUB><IT>− P</IT><SUB>2</SUB>)]<IT>e</IT><SUP>−<IT>n</IT>(1+φ(<IT>P<SUB>r</SUB></IT>))<IT>P<SUB>r</SUB>m</IT>θ</SUP> (4)
Again, the parameters A, P1, P2, m, and theta  are taken from fits of data from the control pathway, and phi (Pr) is taken from the minimal stimulation data by selecting the values of Eq. 3 at P1 and P2. The PPF blocking functions are fitted acceptably with these calculated Sn, so we cannot exclude a bimodal Pr distribution from consideration on the basis of fitting the PPF data. As we noted before, however, that the phi (Pr) obtained from minimal stimulation data are itself a continuous function; in Eq. 4, it is treated artificially as a bimodal function because only two of its points are selected in the calculation of Sn. Ignoring this inherent contradiction yields reasonable fits but indicates a bimodal Pr distribution is a less appropriate choice than a continuous one.

Ca2+ alters the Pr distribution according to the Dodge-Rahamimoff equation

Transmitter release is a Ca2+-dependent process, and changes in external Ca2+ profoundly affect synaptic release probability (Dodge and Rahamimoff 1967). We examined the MK-801 blocking function at four different external Ca2+ concentrations: 1.0, 1.75, 2.5 (standard), and 6.0 mM. Figure 8A shows blocking functions under these four conditions.


View larger version (17K):
[in this window]
[in a new window]
 
FIG. 8. Blocking functions in varying external [Ca2+] follow the Dodge-Rahamimoff equation. A: representative blocking functions for 1.0 mM (square ), 1.75 mM (bullet ), 2.5 mM (triangle , standard), and 6.0 mM (black-square) external [Ca2+]. B: average r-1 plotted against [Ca2+]. Error bars represent SE; data are normalized to r-1 for [Ca2+] = 1.75 mM. Data are fitted by the Dodge-Rahamimoff equation: normalized r-1 = {[1 + (1.3/kMg) + (2.5/kCa)][Ca]/[2.5][1 + ([Mg]/kMg) + ([Ca]/kCa)]}4, where [Mg] and [Ca] are Mg2+ and Ca2+ concentrations, respectively. Best-fit kMg and kCa are 3.9 and 0.57 mM, respectively.

Blocking functions for all Ca2+ concentrations were well fitted by our model incorporating a continuous Pr distribution. The average best-fit rs for external [Ca2+] of 1.0 mM (n = 2), 1.75 mM (n = 5), 2.5 mM (n = 6), and 6.0 mM (n = 4) were 22.51 ± 4.5, 12.67 ± 0.868, 9.78 ± 0.502, and 8.81 ± 0.607, respectively.

The change in r reflects changes in the Pr distribution at different [Ca2+] concentrations. r, however, also is affected by MK-801 block efficacy (recall that r = lambda /mtheta ); because Ca2+ ions may interact with the MK-801 block site (Reynolds and Miller 1988), we examined the EPSCs from each condition to determine whether block fraction had been significantly changed from standard conditions (2.5 mM [Ca2+]). Analysis of the EPSCs revealed that MK-801 block fraction was not significantly affected by changing the external Ca2+ concentration except at 6.0 mM [Ca2+], where the average block fraction was reduced from 0.611 ± 0.068 to 0.330 ± 0.029. Because r is inversely proportional to block fraction, the fitted r for the condition of 6.0 mM [Ca2+] was scaled by a factor of 0.330/0.611 = 0.54 for the purposes of quantitative comparison.

Figure 8B shows normalized r-1 plotted against external [Ca2+]; the data are fitted by the Dodge-Rahamimoff equation (Dodge and Rahamimoff 1967). Best fit by a least-squares fitting procedure yielded kMg and kCa of 3.9 and 0.57 mM, respectively.

Blocking functions at different [Ca2+] also were fitted reasonably with a Sn incorporating a bimodal Pr distribution by varying the three required parameters for a best fit. However, the three parameters of the fitting function, A, b1, and b2, changed in a manner that we were unable to relate to the Dodge-Rahamimoff equation: for external [Ca2+] of 1.0 mM, average A, b1, and b2 were 0.649 ± 0.008, 0.039 ± 0.007, and 0.007 ± 0.002 (n = 2); for 1.75 mM, they were 0.655 ± 0.019, 0.089 ± 0.008, and 0.012 ± 0.001 (n = 5); for 2.5 mM, they were 0.677 ± 0.012, 0.116 ± 0.007, and 0.014 ± 0.0006 (n = 6); and for 6.0 mM, they were 0.749 ± 0.012, 0.224 ± 0.03, and 0.019 ± 0.004 (n = 4). The weighting factor A changed as [Ca2+] changed, implying some synapses move from low to high Pr in this model. Thus we know of no straightforward way to relate the changes of the three parameters of the bimodal distribution to the changes in external calcium concentration.

    DISCUSSION
Abstract
Introduction
Methods
Results
Discussion
References

In the work presented here, we analyzed the relationship between synaptic release probability and the decrease of NMDA-mediated synaptic current amplitudes in the presence of MK-801. Previous studies exploiting this method (Hessler et al. 1993; Rosenmund et al. 1993) have argued that the simplest fit of the MK-801 blocking function points to two classes of synapses with a sixfold difference in release probability (although Rosenmund et al. did not exclude the possibility of multiple Pr classes). The proposal of two synaptic Pr classes, however, is inconsistent with the results of studies in which synaptic Pr is directly measured (Allen and Stevens 1994; Murthy et al. 1997) (see also Fig. 6); such studies indicate synapses range over all possible release probabilities. Furthermore, studies have shown that individual synapses are capable of both increasing and decreasing their Pr by varying amounts under conditions such as long-term potentiation and long-term depression (Bolshakov and Siegelbaum 1995; Stevens and Wang 1994), implying that two Pr classes may be inadequate to describe the available Pr landscape.

We have demonstrated that the data from MK-801 can be described with a continuous Pr distribution. One attractive feature of this distribution is that it describes the MK-801 blocking function with a single parameter (r). In contrast, three parameters are needed with the bimodal distribution (A, b1, b2). Our distribution thus provides the simplest fit of the MK-801 data. Furthermore, this distribution fits the MK-801 blocking function robustly under conditions of altered release probability, such as PPF and changes in external [Ca2+]. In our PPF experiments, we can account for the change in MK-801 blocking function in terms of our distribution and the amount of observed PPF alone. In our Ca2+ experiments, the predicted change in our distribution at different Ca2+ concentrations follows the Dodge-Rahamimoff equation. In contrast, we cannot identify a consistent relationship between the parameters of the bimodal distribution and the [Ca2+].

The equation we use to describe the continuous distribution, rho (Pr) = e-lambda Pr, is, of course, a mathematical idealization. For instance, we actually expect the distribution to go to 0 at Pr = 0, as described in Murthy et al. (1997). Our idealization is acceptable because very low values of Pr are not detectable by the MK-801 method, although this method is more accurate at low Pr than minimal stimulation. We estimate the MK-801 method cannot detect synaptic Pr < 0.05; with minimal stimulation, synapses begin to be heavily underrepresented at Pr of 0.1-0.2. The FM1-43 method used by Murthy et al. (1997) is also inaccurate for low probability synapses and underrepresents synapses with release probabilities less than ~0.05-0.1.

Instead of idealizing the minimal stimulation-derived distribution as a negative exponential, we can express it exactly as a histogram (Allen and Stevens 1994). This "real" distribution is problematic in that it assuredly underestimates the relative number of synapses at the lower end of the Pr range (<0.2, as discussed earlier), so it is of limited use for analysis. Nonetheless, to see whether this distribution is consistent with the MK-801 data, we calculated its Sn and fitted it to the experimental blocking function. The calculated Sn had only two free parameters, which were varied for best fit: theta , the number of receptors that participate in the current, and y, a small constant added to the fitting equation to compensate for the underrepresentation of low Pr synapses. When these parameters were adjusted, the data were fitted quite well (Fig. 4, ···); the best fit yielded theta  = 0.54. Thus a real, albeit incomplete, distribution derived from minimal stimulation data are consistent with the MK-801 blocking function.

Although we have confirmed here that the MK-801 blocking functions reflect the distribution of synaptic reliabilities, the limitations of this method for estimating synaptic reliabilities must be stressed. We consider these limitations in the following text.

We have used a particular continuous weighting function of synaptic reliabilities [rho (Pr) = e-lambda Pr] to describe our data---this distribution is the lowest order approximation to a more general function---but we must emphasize that the MK-801 method alone cannot be used to determine the exact shape of the function that describes the synaptic reliability distribution. The problem is that the MK-801 blocking function can be fitted with sums of anywhere between two (bimodal) and an infinite (continuous) number of exponential functions. Distinguishing whether a function with even a small amount of noise is the sum of 2, 10, or more exponentials is difficult because exponential functions with similar rate constants are very similar in shape. On the basis of our analysis, the reliability distribution rho (Pr) could be a continuous exponential function, a sum of two delta functions, or any function in between (such as sums of multiple delta or Gaussian functions); with correctly chosen parameters, the predicted Sn are nearly indistinguishable. We can determine the range of rho (Pr) that is consistent with the observed blocking function, but within that range, cannot decide whether one rho (Pr) is better based on fit. To make that decision, one must use other criteria, as we have done in the preceding text.

Whatever the precise problems with the MK-801 method for determining the distribution of synaptic reliabilities, we have resolved the conflict between the observed continuous distribution of release probabilities and the original interpretation of the MK-801 blocking data as indicating the presence of two or a few classes of synapses. A firm conclusion of the present work is that synapses are generally quite unreliable and that a broad distribution of reliabilities that is weighted to the less reliable synapses is present in slices. A measure of how much the distribution is weighted to low Prs is provided by the characteristic release probability 1/lambda of our fitted continuous distribution: two-thirds of the synapses have a Pr value less than or equal to the characteristic Pr. Our estimation of 1/lambda is 0.17-0.33. We believe that any theory of brain function will have to take account of this broad distribution of synaptic reliabilities.

    ACKNOWLEDGEMENTS

  This work was supported by the Howard Hughes Medical Institute and a National Institute of Neurological Disorders and Stroke Grant NS-12961 to C. F. Stevens and by a Department of Defense Graduate Fellowship from the Office of Naval Research to E. P. Huang.

    APPENDIX

MK-801 blocking equation

Here we derive the relationship between synaptic Pr and the decrease of EPSC amplitudes in the presence of MK-801. Let us start with the simplest case where a hypothetical population of synapses all have the same release probability. In that case, the fraction Fn of receptors remaining unblocked after n stimulation trials in the presence of MK-801 can be calculated by solving the difference equation
<IT>F<SUB>n<UP>+1</UP></SUB></IT> = (1 − <IT>a</IT>)<IT>F<SUB>n</SUB></IT>
a, the fraction of receptors blocked for any given stimulus, is the product of three factors: Pr, synaptic release probability; theta , the fraction of receptors participating in the current; and m, the fraction of these participating receptors that get blocked by MK-801 (block fraction). In the text, we found m = 0.611 for an MK-801 concentration of 40 µM.

Solving the above difference equation yields
<IT>F<SUB>n</SUB></IT> = (1 − <IT>a</IT>)<SUP><IT>n</IT></SUP>
For small a, such as when release probability Pr is small, this equation becomes
<IT>F<SUB>n</SUB> ≈ e</IT><SUP>−<IT>na</IT></SUP>
Thus if release probability is uniform across a population of synapses and m and theta  are constants, the number of unblocked receptors at trial n will decline as a single exponential. Actual EPSC amplitude at trial n, Sn, will be proportional to the product of the fraction of unblocked receptors (Fn), the probability of release (Pr), and the average current per receptor (z)
<IT>S<SUB>n</SUB> ∼ zP<SUB>r</SUB>F<SUB>n</SUB></IT>
The actual value of Sn can be obtained by scaling this quantity by the total number of receptors. From the preceding text, we see that the assumption of uniform Pr yields a measured quantity Sn that declines exponentially; so for this hypothetical case, the blocking function is an exponential. However, as has been noted previously (Hessler et al. 1993; Rosenmund et al. 1993), the decline of synaptic current amplitudes in the presence of MK-801 cannot be fitted with a single exponential, indicating that release probability is not constant across the population of synapses responsible for the measured synaptic currents.

Because release probability is not constant, there must be a probability density function omega (Pr) that describes the relative distribution of synaptic release probabilities in the population; thus Fn can be weighted accordingly
<IT>F<SUB>n</SUB></IT> = <LIM><OP>∫</OP></LIM><SUP>1</SUP><SUB>0</SUB>ω(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP>d<IT>P</IT><SUB><IT>r</IT></SUB> ≈ <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>ω(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP>d<IT>P</IT><SUB><IT>r</IT></SUB>
The upper limit of the integral can be taken as infinity if omega (Pr) approaches zero rapidly as Pr approaches one, so that higher values of Pr do not contribute; in practical terms, this means there are many more low Pr synapses than high Pr synapses. Several studies (Allen and Stevens 1994; Hessler et al. 1993; Rosenmund et al. 1993) have shown this to be the case.

Again, synaptic current amplitude for trial n will be proportional to the product of Fn, Pr, and the size of the average current as a function of Pr, z(Pr)
<IT>S<SUB>n</SUB></IT> ∼ <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB><IT>z</IT>(<IT>P</IT><SUB><IT>r</IT></SUB>)ω(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>P</IT><SUB><IT>r</IT></SUB><IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP> = <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB>ρ(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP>
where rho (Pr) triple-bond  omega (Pr)z(Pr)Pr.

The actual magnitude of Sn depends on the total number of receptors, but this factor is unimportant in our analysis because we will normalize Sn to its initial value.

Motivation for rho (Pr) as a negative exponential function The choice of an exponential weighting function can be motivated as follows. First, recall
<IT>S<SUB>n</SUB></IT> = <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB>ρ(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP>
which also can be written as
<IT>S<SUB>n</SUB></IT> = <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB><IT>e</IT><SUP>lnρ(<IT>P<SUB>r</SUB></IT>)</SUP><IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP>
Expanding the function lnrho (Pr) to first order about &Pmacr;r, we find the equation
<IT>S<SUB>n</SUB></IT> = <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB><IT>e</IT><SUP>lnρ(<IT><A><AC>P</AC><AC>¯</AC></A><SUB>r</SUB></IT>)+(<IT>P<SUB>r</SUB>−<A><AC>P</AC><AC>¯</AC></A><SUB>r</SUB></IT>)<IT>k</IT>+. . .</SUP><IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP>where <IT>k</IT>is a constant &cjs0358;= <FR><NU>ρ′(<IT><A><AC>P</AC><AC>¯</AC></A><SUB>r</SUB></IT>)</NU><DE>ρ(<IT><A><AC>P</AC><AC>¯</AC></A><SUB>r</SUB></IT>)</DE></FR>&cjs0359;
Generally, rho (Pr) decreases with Pr, so k is negative. Taking the first order approximation, we get
<IT>S<SUB>n</SUB></IT> ∼ <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB><IT>e</IT><SUP>−‖<IT>k</IT>‖<IT>P<SUB>r</SUB></IT></SUP><IT>e</IT><SUP>−<IT>nP<SUB>r</SUB>m</IT>θ</SUP>
Thus the continuous distribution e-lambda Pr can be thought of as the lowest order approximation to the general function rho (Pr). This distribution reasonably describes the dominant tail portion of Pr histograms derived from the study of individual synapses (Allen and Stevens 1994; Murthy et al. 1997).

MK-801 blocking equation for PPF synapses A modified Sn for paired pulse facilitated synapses is derived next. The assumptions and analysis are similar to those used above to derive Sn for normal synapses (see Eq. 1). We write the difference equation that describes the fraction of unblocked receptors at stimulus trial k for the case in which every second pulse is facilitated
<IT>F<SUB>k<UP>+1</UP></SUB> = F<SUB>k</SUB> − r<SUB>k</SUB>F<SUB>k</SUB></IT>,
where rk = a if k is even (first pulse) and rk = alpha  if k is odd (second pulse).

Take n = k/2, so that n counts the first stimulus of each pair. The fraction of unblocked receptors on the first pulse of the nth stimulus pair is
<IT>F<SUB>n</SUB></IT> = (1 − <IT>a</IT>)<SUP><IT>n</IT></SUP>(1 − α)<SUP><IT>n</IT></SUP> ≈ <IT>e</IT><SUP>−(<IT>a</IT>+α)<IT>n</IT></SUP>
Using the same reasoning as in the text, the actual synaptic response size for the first stimulus of the nth pair is given by
<IT>S<SUB>n</SUB></IT> ∼ <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB>ρ(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>e</IT><SUP>−(<IT>a</IT>+α)<IT>n</IT></SUP>
Remember that the fraction blocked a on the first pulse equals Prmtheta . For the second (facilitated) pulse, the fraction blocked alpha  is related to a by
α = φ(<IT>P<SUB>r</SUB></IT>)<IT>a</IT>
where phi (Pr) is the PPF ratio at a given Pr. Thus
<IT>S<SUB>n</SUB></IT> = <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB>ρ(<IT>P</IT><SUB><IT>r</IT></SUB>)<IT>e</IT><SUP>−<IT>n</IT>(1+φ(<IT>P<SUB>r</SUB></IT>))<IT>P<SUB>r</SUB>m</IT>θ</SUP>
With a continuous distribution rho (Pr) = e-lambda Pr
<IT>S<SUB>n</SUB></IT> = <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB><IT>e</IT><SUP>−λ<IT>P<SUB>r</SUB></IT></SUP><IT>e</IT><SUP>−<IT>n</IT>(1+φ(<IT>P<SUB>r</SUB></IT>))<IT>P<SUB>r</SUB>m</IT>θ</SUP>
With a bimodal distribution rho (Pr) = Adelta (Pr - P1) + (1 - A)delta (Pr - P2)
<IT>S<SUB>n</SUB></IT> = <LIM><OP>∫</OP></LIM><SUP>∞</SUP><SUB>0</SUB>d<IT>P</IT><SUB><IT>r</IT></SUB>[<IT>A</IT>δ(<IT>P</IT><SUB><IT>r</IT></SUB><IT> − P</IT><SUB>1</SUB>) + (1 − <IT>A</IT>)δ(<IT>P</IT><SUB><IT>r</IT></SUB><IT> − P</IT><SUB>2</SUB>)]<IT>e</IT><SUP>−<IT>n</IT>(1+φ(<IT>P<SUB>r</SUB></IT>))<IT>P<SUB>r</SUB>m</IT>θ</SUP>

    FOOTNOTES

  Address for reprint requests: C. F. Stevens, The Salk Institute, 10010 N. Torrey Pines Rd., La Jolla, CA 92037.

  Received 24 March 1997; accepted in final form 13 August 1997.

    REFERENCES
Abstract
Introduction
Methods
Results
Discussion
References

0022-3077/97 $5.00 Copyright ©1997 The American Physiological Society