PLoS Computational Biology

PLoS Comput Biol. 14(6): e1006232

Transmission of temporally correlated spike trains through synapses with short-term depression

http://orcid.org/0000-0003-4895-8253Alex D. Bird, Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing¹²³*, Magnus J. E. Richardson, Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing⁴*

Farzan Nadim, Editor

1. Warwick Systems Biology Centre, University of Warwick, Coventry, United Kingdom

2. Ernst Strüngmann Institute for Neuroscience, Max Planck Society, Frankfurt, Germany

3. Frankfurt Institute for Advanced Studies, Frankfurt, Germany

4. Warwick Mathematics Institute, University of Warwick, Coventry, United Kingdom

New Jersey Institute of Technology / Rutgers Univ-Newark, UNITED STATES

The authors have declared that no competing interests exist.

* E-mail: alex.neurosci@gmail.com (ADB); magnus.richardson@warwick.ac.uk (MJER)

© 2018 Bird, Richardson

DOI: 10.1371/journal.pcbi.1006232

Published online: 22 June 2018

Abstract

Short-term synaptic depression, caused by depletion of releasable neurotransmitter, modulates the strength of neuronal connections in a history-dependent manner. Quantifying the statistics of synaptic transmission requires stochastic models that link probabilistic neurotransmitter release with presynaptic spike-train statistics. Common approaches are to model the presynaptic spike train as either regular or a memory-less Poisson process: few analytical results are available that describe depressing synapses when the afferent spike train has more complex, temporally correlated statistics such as bursts. Here we present a series of analytical results—from vesicle release-site occupancy statistics, via neurotransmitter release, to the post-synaptic voltage mean and variance—for depressing synapses driven by correlated presynaptic spike trains. The class of presynaptic drive considered is that fully characterised by the inter-spike-interval distribution and encompasses a broad range of models used for neuronal circuit and network analyses, such as integrate-and-fire models with a complete post-spike reset and receiving sufficiently short-time correlated drive. We further demonstrate that the derived post-synaptic voltage mean and variance allow for a simple and accurate approximation of the firing rate of the post-synaptic neuron, using the exponential integrate-and-fire model as an example. These results extend the level of biological detail included in models of synaptic transmission and will allow for the incorporation of more complex and physiologically relevant firing patterns into future studies of neuronal networks.

Author summary

Synapses between neurons transmit signals with strengths that vary with the history of their activity, over scales from milliseconds to decades. Short-term changes in synaptic strength modulate and sculpt ongoing neuronal activity, whereas long-term changes underpin memory formation. Here we focus on changes of strength over timescales of less than a second caused by transitory depletion of the neurotransmitters that carry signals across the synapse. Neurotransmitters are stored in small vesicles that release their contents, with a certain probability, when the presynaptic neuron is active. Once a vesicle has been used it is replenished after a variable delay. There is therefore a complex interaction between the pattern of incoming signals to the synapse and the probablistic release and restock of packaged neurotransmitter. Here we extend existing models to examine how correlated synaptic activity is transmitted through synapses and affects the voltage fluctuations and firing rate of the target neuron. Our results provide a framework that will allow for the inclusion of biophysically realistic synaptic behaviour in studies of neuronal circuits.

Introduction

Variability in synaptic function arises from stochasticity in processes ranging in scale from the transitory opening and closing of ion channels to probabilistic neurotransmitter release and vesicle restock [1–4]. The transmission of signals between neurons is therefore inherently stochastic and, moreover, will interact in a history-dependent manner with the patterns in the incoming presynaptic drive [5–9]. A common approach to treating this stochasticity analytically assumes that neuronal firing is uncorrelated, with a Poisson process typically used to model spike times [10–12]. However, non-Poissonian activity is regularly observed in vivo [13–16] and models suggest that, even in the absence of short-term synaptic dynamics, it can have a substantial effect on the propagation of activity [17–21]. It can be expected that the impact of non-Poissonian presynaptic activity will be further complicated when combined with vesicle-depletion depression, in which synaptic transmission becomes weaker and less reliable as stores of available neurotransmitter are depleted and yet to be restocked [22–24].

Transmission through plastic synapses has been shown to decorrelate input spike trains [5, 25, 26], typically increasing the computational power [27, 28] and efficiency [29] of a neuronal network. While average rate effects under the influence of depression have been extensively studied [30–32], a compact set of analytical results for correlated spike trains and stochastic quantal vesicle release from multiple sites has remained elusive (for a full discussion of existing results, see the Discussion). Recently, it has been shown [33] that correlated firing patterns more regular than Poissonian can increase the rate of vesicle release, thereby enhancing the fidelity and efficiency of signal transmission, whilst more irregular spike trains can lead to a decrease in neurotransmitter release [25, 27].

To further analyse how the interaction of correlated presynaptic drive and short-term depression affect quantal synaptic transmission, here we derive a number of analytic results for renewal processes, for which the incoming spike train is fully characterised by the inter-spike-interval (ISI) distribution. This type of presynaptic drive includes that generated by the leaky, quadratic and exponential integrate-and-fire models driven by white-noise [34–36], which are models commonly used to fit experimental data [37–39]. It can be noted that these models will also generate ISIs that are well approximated by renewal processes when the correlations of their incoming synaptic drive are much shorter in time than the typical outgoing ISIs.

We derive equations for the occupancy and the temporal structure of release events when presynaptic cells have spiking patterns fully characterised by their inter-spike-interval distribution. We then show how these can be used to calculate the post-synaptic voltage mean and variance when the presynaptic neurons make multiple independent contacts. These results are illustrated using gamma-distributed ISIs and presynaptic integrate-and-fire neurons. We show that the results allow for a good estimation of post-synaptic firing rates thereby opening the way for the analysis of the role synapses with stochastic short-term depression play in feed-forward and recurrent neuronal networks in which presynaptic spike-trains are typically non-Poissonian.

Materials and methods

Synaptic release sites

A quantal model of synaptic dynamics is used where the binary variable x represents the occupancy (x = 1) of a release site by a neurotransmitter-filled vesicle, with x = 0 otherwise. On the arrival of a presynaptic action potential, neurotransmitter is released with probability p if a vesicle is present at the site. For brevity, the probability that a present vesicle is not released, 1 − p, is written as q. Sites that are empty are then restocked at memoryless rate λ (Poisson process). Note that because x is a binary variable taking values of 0 and 1, x² ≡ x and so Var(x) = 〈x〉(1 − 〈x〉) where the notation 〈X〉 is used as the expectation of any quantity X over all stochastic processes (spike times, release and restock events). Values of parameters used to generate figures, unless otherwise stated, are given in Table 1.

See full tableTable 1. Parameters and their values used for figures.Parameters are grouped together as: those pertaining to occupancy and release statistics; those required for calculation of the post-synaptic response; and those used to generated ISI statistics. Certain neuronal parameters, such as τ and μ, are used for both pre and post-synaptic integrate-fire models.

ISI distribution

The presynaptic spike train is modelled as a renewal process characterised by the ISI distribution f(t) where the firing rate r is the reciprocal of the mean ISI. Though there are no correlations between successive ISIs, the spike train itself will in general be non-Poissonian and can range between bursting and regular extremes: bursting neurons have positively autocorrelated spike trains at short times whereas more periodically firing neurons generate trains that are negatively autocorrelated over short time frames. Many results in this paper will involve expectations over the ISI distribution of exponential functions

[Image: pcbi.1006232.e001.jpg]⟨e-zt⟩=∫0∞f(t)e-ztdt=L(z)(1)

where z can be complex (with restrictions on its real part related to the specific functional form the ISI distribution). This can also be interpreted as the Laplace transform of the ISI distribution, which we will denote as [Image: pcbi.1006232.e002.jpg]ℒ(z). For Laplace transforms of other quantities, say X(t), we will use the notation [Image: pcbi.1006232.e003.jpg]ℒX(t). Note that the Fourier transform [Image: pcbi.1006232.e004.jpg]f^(ω) of the ISI distribution can also be interpreted as an expectation

[Image: pcbi.1006232.e005.jpg]f^(ω)=∫-∞∞f(t)e-iωtdt=⟨e-iωt⟩=L(iω)(2)

and is directly related to its Laplace transform, remembering that f(t) = 0 for t < 0. This will allow many existing results from the literature on integrate-and-fire ISI distributions in the frequency domain to be used in the subsequent analyses.

Gamma-distributed ISI distributions

Gamma distributed ISIs provide a useful illustration of the results presented here as a single shape parameter α continuously varies the train between bursting for α < 1 and more regular α > 1. The ISI distribution f(t) is given by

[Image: pcbi.1006232.e006.jpg]f(t)=(αr)αΓ(α)tα-1e-αrt(3)

for positive t and zero otherwise, where [Image: pcbi.1006232.e007.jpg]Γ(α)=∫0∞tα-1e-tdt is the gamma function. When α = 1 the ISI distribution becomes exponential and the presynaptic spike train is a Poisson process. The expectation of an exponential function (Eq 1) over this class of ISI distribution is

[Image: pcbi.1006232.e008.jpg]⟨e-zt⟩=(αrz+αr)α(4)

where z can be a complex constant as long as the real part of z + αr is greater than zero.

Data availability

Together with this paper we provide JULIA code in the Jupyter Notebook environment for generating each of the figures shown in the paper. All five scripts are published under the GNU General Public Licence, Version 3 (http://www.gnu.org/copyleft/gpl.html).

Results

After writing down some general statements for arbitrary spike times, we focus on deriving exact results for the case when correlated presynaptic spiking is a renewal process and fully described by the ISI distribution. We present formulae for the steady-state vesicle-release site occupancy averaged over time as well as its mean value just before the arrival of a presynaptic action potential. We then derive an integral equation for the occupancy at some later time given a release event at an earlier time: this will allow us to calculate the autocovariance of release events which in turn leads to an analytical form for the postsynaptic voltage variance. We then generalise this result to a scenario in which each presynaptic cell makes multiple independent contacts. Because of the shared presynaptic drive across these contacts, an additional level of correlation is generated in the input to the post-synaptic cell. We characterise these correlations through the cross-covariance of release events and extend the formula for the post-synaptic voltage to multiple contacts. These results are illustrated using gamma-distributed ISIs for the presynaptic trains. It is then further shown how all results can be derived exactly for presynaptic LIF models or numerically for other classes of integrate-and-fire model, such as the EIF model that are themselves all driven by Gaussian white-noise drive. Finally, we consider a straightforward extension to consider biophysically realistic post-synaptic potentials.

Average and pre-spike mean release-site occupancy

The statistics of a binary occupancy variable x(t), where x = 1 if the site is occupied and is zero otherwise, are first considered. We first write down an obvious steady-state result which links two averages of this quantity: 〈x〉, the occupancy averaged over time, and 〈x〉∞, the mean occupancy just before the arrival of a presynaptic spike. It should be noted that these two quantities are only the same for (memory-less) Poisson processes. Given that the total restock rate must equal the total release rate in the steady state, the following balance equation holds

[Image: pcbi.1006232.e009.jpg]λ(1 - ⟨x⟩) = pr⟨x⟩∞(5)

where λ is the restock rate given no vesicle is present, p is the probability of release given a vesicle is present and r is the firing rate of the presynaptic neuron. As will be seen later, it is the quantity 〈x〉∞ that is required for analysing the effect on the post-synaptic cell.

To derive 〈x〉∞ it is convenient to first consider a less complete expectation of x, denoted by [Image: pcbi.1006232.e010.jpg]x¯, which implies the expectation of x for a fixed pattern of spike times {t₁, t₂, …tm} but averaged over all possible patterns of restock and release events. Consider the expected occupancy [Image: pcbi.1006232.e011.jpg]x¯m immediately before the mth spike at time tm as a function of the expected occupancy [Image: pcbi.1006232.e012.jpg]x¯m-1 immediately before the (m − 1)th spike: this obeys the recursion equation

[Image: pcbi.1006232.e013.jpg]x¯m=x¯m-1qe-λ(tm-tm-1)+(1-e-λ(tm-tm-1))(6)

where q = 1 − p. This can be solved, with the initial condition [Image: pcbi.1006232.e014.jpg]x¯1=1, to give

[Image: pcbi.1006232.e015.jpg]x¯m=1-pq∑k=1m-1qke-λ(tm-tm-k).(7)

Taking expectations over all realisations of presynaptic spike times, as m → ∞, gives the expected occupancy before the arrival of a presynaptic action potential in the steady-state

[Image: pcbi.1006232.e016.jpg]⟨x⟩∞=1-pq∑k=1∞qk⟨e-λTk⟩(8)

where Tk is the sum of the last k ISIs. This result is quite general and holds for arbitrarily correlated spike trains; however, we now consider the specific case of spike-times generated by a renewal process. In this case the ISIs will be independent so that the expectation over the exponential term factorises [Image: pcbi.1006232.e017.jpg]〈e-λTk〉 = 〈e-λt〉k and the sum can be evaluated for 〈x〉∞ to give

[Image: pcbi.1006232.e018.jpg]⟨x⟩∞=1-⟨e-λt⟩1-q⟨e-λt⟩and⟨x⟩=1-pr(1-⟨e-λt⟩)λ(1-q⟨e-λt⟩).(9)

The corresponding form for 〈x〉 was found from Eq (5). The expectation 〈e−λt〉 is straightforward to evaluate for many classes of renewal process and is directly related to the ISI Laplace or Fourier transform (Eqs 1 and 2) [Image: pcbi.1006232.e019.jpg]ℒ(λ) = 〈e-λt〉.

Release-site occupancy for gamma-distributed ISIs

To illustrate the theoretical results in this paper we consider ISIs that have a distribution given by a gamma function with shape factor α and mean rate r (Eq 3). This is a convenient choice because a value α < 1 generates a spike train that is bursty, α = 1 is Poissonian firing and α > 1 is regular firing, thereby allowing for a range of behaviours to be examined as a single parameter is varied (at fixed presynaptic rate). Illustrations of the behaviour for three different values of α, from bursting to regular, are provided in Fig 1A with the corresponding ISI distribution (Eq 3) and its cumulation given in Fig 1B and 1C. For this class of ISI distribution we have

[Image: pcbi.1006232.e020.jpg]⟨e-λt⟩=(αrαr+λ)α(10)

from Eq (4). This result can be used to calculate the forms of the average occupancy and mean pre-spike occupancy given in equation pair (9) and plotted as a function of α in Fig 1D. As α increases, the pre-spike mean 〈x〉∞ increases and the overall mean 〈x〉 decreases. The two means take the same value when α = 1 and the input spike train is an uncorrelated Poisson process. It is interesting to consider the two extreme limits of the quantity in Eq (10)

[Image: pcbi.1006232.e021.jpg]⟨e-λt⟩→1-αlog(λ/αr)whenα→0and⟨e-λt⟩→e-λrwhenα→∞(11)

for highly bursty and highly regular presynaptic firing, respectively. In the limit α → 0 of very bursty presynaptic firing the mean occupancies tend to

[Image: pcbi.1006232.e022.jpg]⟨x⟩∞→αplog(λ/αr)and⟨x⟩→1-rλαlog(λ/αr).(12)

Hence, in the bursty limit the dominant behaviour of 〈x〉∞ is as [Image: pcbi.1006232.e023.jpg]-αplog(α), which is dependent on p and α only, and the probability of neurotransmitter release p〈x〉∞ is a function of α only. Both these quantities and the variance of the prespike occupancy vanish in the extreme bursty limit α → 0, a trend that can be seen in Fig 1D, 1E and 1F.

Fig 1View larger versionFig 1. Vesicle occupancy and neurotransmitter release rate are a function of the ISI distribution.(A) Presynaptic spike train, vesicle occupancy and neurotransmitter release time-course for gamma-distributed ISIs with α as marked (in panel B) for the same release probabilty p = 0.6, presynaptic rate r = 5Hz and postsynaptic rate λ = 2Hz. (B) Corresponding ISI distributions for the three α values. (C) The cumulative distribution. (D) The mean 〈x〉 (grey) and the pre-spike 〈x〉∞ (black) release-site occupancies (Eqs 9 and 10). Note that 〈x〉∞ increases with increasing regularity, larger α. (E) Variance of the two occupancies exhibiting non-monotonic behaviour. (F) The mean neurotransmitter release rate 〈χ〉 is directly proportional to 〈x〉∞ and so shares its qualitative dependence on α. The code used to generate this figure is provided in the Supporting Material.

As the presynaptic firing becomes increasingly regular, for large α, the quantity 〈e−λt〉 begins to lose its α dependency. For the particular case of regular firing when λ/r ≪ 1, meaning that the restock rate is low relative to the presynaptic firing rate and the vesicles are in a depleted state the limit is 〈x〉∞ → (1/p)(λ/r) and so in this case the probability of neurotransmitter release p〈x〉∞ → λ/r is independent of both p and α and that the total mean rate of vesicle release rp〈x〉∞ is a function of λ only. This observation, that when r ≫ λ, the release rate loses its dependency on the presynaptic rate has been previously highlighted [32] for Poissonian processes (when α = 1). However, it is interesting to note that for bursty firing the release rate does not lose its dependency on the presynaptic rate so readily in this limit, as argued above.

Vesicle release average and autocovariance

We now consider the statistics of the release of neurotransmitter and, in particular, the autocovariance of the release. This quantity will be required to calculate the variance of the postsynaptic voltage. We define the release events at a single site as a series of Dirac-delta events

[Image: pcbi.1006232.e024.jpg]χ(t)=∑{tk}δ(t-tk)(13)

where here {tk} are the times of the neurotransmitter release events. These occur with probability p only when a presynaptic spike arrives and a vesicle is present, so that the steady-state mean of this quantity is 〈χ〉 = pr〈x〉∞, as already observed in the previous section.

We now consider the steady-state autocovariance that, because of its time-translation invariance, can be written 〈χ(t)χ(0)〉 − 〈χ〉². For t > 0 we note that 〈χ(t)χ(0)〉 = 〈χ(t|0)〉〈χ〉 where 〈χ(t|0)〉 is understood to be the probability density of a vesicle being released at time t given that one was released previously at time 0. The autocovariance can therefore be written

[Image: pcbi.1006232.e025.jpg]AutoCov(χ) = ⟨χ⟩(δ(t) + ⟨χ(t|0)⟩ - ⟨χ⟩)(14)

for t ≥ 0 (the function is even in time, thus specifying the t < 0 component) with the Dirac delta function coming from the zero-time contribution.

The quantity 〈χ(t|0)〉 itself is the rate of release, given a release at t = 0 and can be written as pG(t) where G(t) is the probability density of presynaptic spike arriving while a vesicle is present at time t given that initially (at t = 0) a presynaptic spike arrived and immediately after there was no vesicle present. We also introduce a related quantity H(t) which has the same conditionality but is the density of a presynaptic spike arriving while there is no vesicle present at a time t. Note that F(t) = G(t) + H(t) where F(t) is the probability density of an action potential arriving at t given there was one at t = 0 (the conditionality is the same because the arrival of a spike at t does not depend on whether a release site was stocked or not just before an earlier spike). This last quantity can be directly related to the ISI distribution f(t) via a convolution

[Image: pcbi.1006232.e026.jpg]F(t)=f(t)+∫0tdsF(s)f(t-s)(15)

that can be solved in terms of integral transforms. It is straightforward to derive similar formulae for G(t) and H(t) via the introduction of the quantities

[Image: pcbi.1006232.e027.jpg]g(t)=f(t)(1-e-λt)andh(t)=f(t)e-λt(16)

where g(t) and h(t) are the probability densities that the release site (initially empty at t = 0 following a presynaptic spike) are either restocked or not, respectively, by the time the next spike arrives; note also that f(t) = g(t) + h(t). To derive a self-consistent integral equation for G(t) we need to decompose it into the various histories that start with a vesicle absent at t = 0 following a presynaptic spike and end with a vesicle present at t when a spike arrives. There are four distinct contributions that need to be accounted for. The first is straightforward as it arises from the first spike arriving at t giving a contribution of g(t). The other contributions imply that the penultimate spike arrives at an intermediate time s that needs to be integrated over (as in Eq 15). The second contribution has a vesicle present before the intermediate time s and no release and so contributes G(s)qf(t − s). The third contribution has a vesicle present before the intermediate time s and there is a release and so contributes G(s)pg(t − s). The final contribution has no vesicle present before the intermediate time s and then there is a restock, contributing H(s)g(t − s). Combining these four contributions and simplifying, using H = F − G and h = f − g, results in the following integral equation

[Image: pcbi.1006232.e028.jpg]G(t)=g(t)+∫0tdsF(s)g(t-s)+q∫0tdsG(s)h(t-s).(17)

The equations for F(t) and G(t) can either be solved numerically using iterative procedures or, alternatively, solved using integral transforms (see the next section). Noting that in the limit t → ∞ the conditional quantity G(t) converges to r〈x〉∞ allows the autocovariance

[Image: pcbi.1006232.e029.jpg]AutoCov(χ(t)) = pr⟨x⟩∞(δ(t) + p(G(|t|) - r⟨x⟩∞))(18)

of the neurotransmitter release time series χ(t) to be written in terms of G(t) and the occupancy 〈x〉∞.

Vesicle-release train in the Laplace domain and its power spectrum

As will be seen later, the Laplace transform of the densities F(t) and G(t) are required for the calculation of the post-synaptic voltage variance. The Laplace transforms of g(t) and h(t) can be written in terms of the ISI-distribution Laplace transform

[Image: pcbi.1006232.e030.jpg]Lg(z)=L(z)-L(z+λ)andLh(z)=L(z+λ).(19)

With these results and the convolution theorem, it is straightforward to solve the integral Eqs (15) and (17) in the Laplace domain.

[Image: pcbi.1006232.e031.jpg]LF(z)=L(z)1-L(z)andLG(z)=1(1-L(z))L(z)-L(z+λ)(1-qL(z+λ)).(20)

These results will be of use for the later calculation of the voltage variance. It is also useful to state the solution in the Fourier domain, as this is required for the power spectra of the spike train and the release process. Because both F(t) and G(t) tend to a constant as t → ∞ and so the Fourier transform will diverge at zero frequency, we separate the solutions into non-zero and zero frequency components to give the Fourier transforms as [Image: pcbi.1006232.e032.jpg]F^(ω)=LF(iω)+rπδ(ω) and [Image: pcbi.1006232.e033.jpg]G^(ω)=LG(iω)+r〈x〉∞πδ(ω). These results can be used [40] to extend the autocorrelation of reference [29] to account for the biophysically important case of unreliable vesicle release. The power spectrum [40] of the presynaptic action potentials [Image: pcbi.1006232.e034.jpg]S^ϱ as well as that for the release events [Image: pcbi.1006232.e035.jpg]S^χ follows as

[Image: pcbi.1006232.e036.jpg]S^ϱ(ω)=r(1+2ℜLF(iω))andS^χ(ω)=⟨χ⟩(1+2pℜLG(iω)))(21)

where 〈χ〉 = pr〈x〉∞.

Auto-covariance and power spectrum for gamma-distributed ISIs

For the particular case of gamma-distributed ISIs, the integral equation (Eq 15) for F(t) can be found by iteratively substituting for F(s) under the integral to provide a series of multiple integrals over products of f(s). On substituting for the gamma-distribution form for f(t), and using standard results for the normalisation of Dirichlet distributions, the following series solution is found

[Image: pcbi.1006232.e037.jpg]F(t)=e-αrtt∑m=1∞(αrt)mαΓ(mα)(22)

which converges relatively rapidly. We were unable to find a similarly compact series solution for G(t), given in Eq (17); however, it is straightforward to develop a numerical scheme that solves the integral equation by iterating forward in time on a grid of time points. Fig 2A plots the presynaptic spike-triggered average rate F(t) together with the corresponding release-triggered average rate pG(t) for three choices of α ranging from bursty to regular (the same as those used in Fig 1). The late-time asymptote for F(t) is the presynaptic rate whereas for pG(t) it is the average release rate 〈χ〉. The corresponding finite component of the autocovariance for these processes is plotted in Fig 2B for positive time (they are even in time). Note that though the presynaptic rate is positively correlated for bursty spike trains, its autocovariance becomes negative at short times (like for the Poissonian and regular trains) once it has passed through synapses exhibiting short-term depression.

Fig 2View larger versionFig 2. Temporal and spectral statistics of the presynaptic spike-train and neurotransmitter release required for post-synaptic voltage mean and variance.(A) The spike-train (Ai) and neurotransmitter release (Aii) event-triggered rate for three cases ranging from bursty to regular with α as marked (see Eq 22 and surrounding text). (B) The equivalent autocovariances of the spike train (Bi) and synaptic release events (Bii). For bursty presynaptic trains (small α = 0.4) the autocovariance is positive; however all cases have negative autocovariances for the release events themselves (Eq 18) because of the filtering by the depressing synapses. (C) The corresponding power spectra (Eq 21) of the spike train (Ci) and synaptic release events (Cii). (D) The postsynaptic voltage mean, standard deviation (std) and coefficient of variation std/mean (Eqs 24 and 29). Though the mean and std both increase with increasing presynaptic regularity, mirroring the behaviour of 〈x〉∞ (Fig 1D), the voltage CV itself decreases with increasing regularity. For panels A-C parameters and color are same as Fig 1, and in panel D parameter α is varied with other parameters provided in Table 1. The code used to generate this figure is provided in the Supporting Material.

The postsynaptic response

In the previous section the pre-spike occupancy of a release site was calculated and the release-train autocovariance derived. We now use these results to derive the voltage mean and variance of the postsynaptic neuron. The subthreshold mean and variance of the postsynaptic neuron are the key quantities necessary to estimate the output firing rate of a neuron [18, 36, 41–45] and hence its resultant effect on the network.

It is assumed that the presynaptic neurons are uncorrelated and each fires at a rate r with the same ISI distribution. For simplicity we assume that each neurotransmitter-release event causes the postsynaptic membrane to increase by a fixed voltage a (this restriction is for simplicity and can be relaxed). In this section we consider presynaptic cells that make contacts with only one vesicle release site, leaving multiple release sites to a later section. The postsynaptic voltage therefore follows the equation

[Image: pcbi.1006232.e038.jpg]τdvdt=μ-v+aτ∑i=1Nχi(t)(23)

where τ is the membrane time constant, μ the resting potential in absence of synaptic drive and there are N presynaptic neurons, with the ith having a release train χi(t).

The steady-state mean post-synaptic voltage is found using the result 〈χ〉 = pr〈x〉∞ so that

[Image: pcbi.1006232.e039.jpg]⟨v⟩ = μ + aτNpr⟨x⟩∞.(24)

This is an increasing function of the prespike occupancy 〈x〉∞ and therefore also increases with the regularity of the presynaptic spike train (see Fig 1D for the dependence of 〈x〉∞ on the burstiness of the presynaptic spike train).

To find the post-synaptic voltage variance we first solve differential Eq (23) formally as an integral over the release trains

[Image: pcbi.1006232.e040.jpg]v(t)=⟨v⟩+a∫-∞tdt′e-(t-t′)/τ∑i(χi(t′)-⟨χ⟩)(25)

where we included the component of the mean stemming from the synaptic input in the summation. The voltage variance can therefore be written as a double integral over the autocovariance of χ(t) as follows

[Image: pcbi.1006232.e041.jpg]Var(v)=Na2⟨χ⟩∫0∞ds∫0∞ds′e-s/τe-s′/τ(δ(s-s′)+p(G(|s-s′|)-r⟨x⟩∞))(26)

where the fact that the χi(t) from different presynaptic neurons are uncorrelated has been used. The Dirac-delta component is straightforward to evaluate and the other components are symmetric around s − s′ so that

[Image: pcbi.1006232.e042.jpg]Var(v)=τNa22⟨χ⟩+2Na2⟨χ⟩p∫0∞ds′e-s′/τ∫s′∞dse-s/τ(G(|s-s′|)-r⟨x⟩∞)(27)

which, on a change of the inner integration variable z = s − s′, allows for the outer integral over s′ to be performed resulting in

[Image: pcbi.1006232.e043.jpg]Var(v)=τNa22⟨χ⟩(1+2p∫0∞dze-z/τ(G(z)-r⟨x⟩∞)).(28)

Part of the integral in the above equation is simply the Laplace transform of G(t), which was already provided in the second of equation pair (20). On substitution, this gives a compact form for the postsynaptic voltage variance in terms of [Image: pcbi.1006232.e044.jpg]ℒ(z) = 〈e-zt〉 the Laplace transform of the ISI distribution

[Image: pcbi.1006232.e045.jpg]Var(v)=τNa22⟨χ⟩(1+2p(L(1/τ)-L(1/τ+λ)(1-L(1/τ))(1-qL(1/τ+λ))-τr⟨x⟩∞)).(29)

This is a central result of this paper and is general for neurons that receive synapses with single release sites from presynaptic neurons that fire as a renewal process. We now go on to examine the postsynaptic voltage behaviour for gamma-distributed presynaptic ISIs, and consider the case for presynaptic integrate-and-fire models in a later subsection.

Post-synaptic voltage statistics for gamma-distributed ISIs

The expectations appearing in the voltage mean and variance equation are of the form of Eq (4) for the case of gamma-distributed ISIs. The mean post-synaptic voltage is a linear function of 〈x〉∞ and therefore also of 〈χ〉 (Fig 1D and 1F), which in turn take their α dependence from Eq (10). The voltage mean therefore increases monotonically with increasing regularity of the presynaptic train, as can be seen in Fig 2D. For the calculation of the variance (Eq 29) two additional expectations are required

[Image: pcbi.1006232.e046.jpg]L(1/τ)=⟨e-t/τ⟩=(ατr1+ατr)αandL(1/τ+λ)=⟨e-(1/τ+Rr)T⟩=(ατr1+τλ+ατr)α.(30)

Using these results the voltage standard-deviation (std) is plotted in Fig 2D along with the coefficient of variation (CV) of the post-synaptic voltage (std/mean). While the variance increase with increasing regularity of the presynaptic train passing through depressing synapses, it can be noted that the voltage CV itself decreases (last panel of Fig 2D).

Multiple release sites sharing the same presynaptic neuron

A single presynaptic neuron will typically make contacts with a postsynaptic cell that result in multiple vesicle release sites, with estimates of this parameter varying from 1 to as much as 100 [46] for neocortical layer-5 pyramidal cells. Here we consider a case where each of the N presynaptic cells makes connections with n indepedent release sites, a scenario illustrated in Fig 3A for a case N = 1 and n = 3. We assume that all processes (such as restock and release) are statistically independent at each of these sites, but those sharing the same presynaptic neuron receive the same spike train. The postsynaptic-voltage dynamics now take the form

[Image: pcbi.1006232.e047.jpg]τdvdt=μ-v+aτ∑i=1N∑j=1nχij(t)(31)

where here χij is the release time course of the jth contact on the ith neuron. Only the χij that share the same presynaptic neuron will be correlated. The steady-state mean voltage is straightforward to derive using the result 〈χ〉 = pr〈x〉∞ so that

[Image: pcbi.1006232.e048.jpg]⟨v⟩ = μ + aτNnpr⟨x⟩∞.(32)

However, to calculate the variance correlations between release sites need to be accounted for, which requires the solution of an additional integral equation.

Fig 3View larger versionFig 3. Correlated occupancy, release and post-synaptic voltage statistics when there are multiple contacts per presynaptic cell.(A) Illustration of a case with one presynaptic cell N = 1 making three contacts n = 3. Shown in descending order: presynaptic spike train (α = 0.4, r = 5Hz), occupancy of the three release sites (restock λ = 2Hz), the corresponding release-event time courses (p = 0.6), and the voltage time course (see Table 1 for parameters). (B) The occupancy correlation 〈xz〉∞ for a pair of sites receiving the same presynaptic spike train and the covariance [Image: pcbi.1006232.e049]〈xz〉∞-〈x〉∞2 (see Eqs 34 and 35). (C) The cross-covariance of release events from a pair of sites sharing the same presynaptic cell (Eq 39). (D) Post-synaptic voltage standard deviation (top: from Eq 42) and CV (bottom) as a function of α for four combinations of N, n as marked. With each of these choices, where Nn = 1000, the mean voltage (Eq 32) is the same as in Fig 2D. The code used to generate this figure is provided in the Supporting Material.

Expectation of joint vesicle occupancy for release sites

The jth contact of presynaptic cell i will have a release train χij(t) that is correlated with that of other contacts from the same presynaptic cell due to receiving the same spike train. To take this into account, the first quantity to consider is the joint prespike occupancy 〈xz〉∞ for two release sites (labelled x(t) and z(t)) sharing the same presynaptic cell. Given that restock and release events are uncorrelated, and that only the spike times are, it is possible to simply multiply the result Eq (6) by the equivalent for [Image: pcbi.1006232.e050.jpg]z¯m and take expectations

[Image: pcbi.1006232.e051.jpg]⟨xz⟩m = ⟨xz⟩m-1q²⟨e-2λt⟩ + 2⟨x⟩m-1q⟨e-λt(1 - e-λt)⟩ + ⟨(1-e-λt)²⟩(33)

where it should be noted that first-order expectations 〈x〉m = 〈z〉m are identical for the two release sites as they are statistically indistinguishable. In the steady-state limit m → ∞ the difference Eq (33) results in

[Image: pcbi.1006232.e052.jpg]⟨xz⟩∞=2q⟨x⟩∞⟨(1-e-λt)e-λt⟩+⟨(1-e-λt)2⟩1-q2⟨e-2λt⟩(34)

which gives the covariance as

[Image: pcbi.1006232.e053.jpg]⟨xz⟩∞-⟨x⟩∞2=p2⟨e-2λt⟩-⟨e-λt⟩2(1-q2⟨e-2λt⟩)(1-q⟨e-λt⟩)2(35)

where the expectations are Laplace transforms of the ISI distribution [Image: pcbi.1006232.e054.jpg]〈e-mλt〉 = ℒ(mλ). It will be necessary in the next section to find the steady-state probability that site z is still stocked given a vesicle was just released from site x and also the probability that site z is unstocked with the same conditionality. These probabilities are

[Image: pcbi.1006232.e055.jpg]q⟨xz⟩∞/⟨x⟩∞and1-q⟨xz⟩∞/⟨x⟩∞(36)

respectively.

Cross-covariance function for release sites

To compute the cross-covariance for two release sites sharing the same presynaptic cell requires the introduction of a quantity G′(t), which is the probability density that a spike arrives while a vesicle is present at time t given that the site was occupied immediately after a spike at time 0. Note that this differs from G(t) defined earlier only by the initial condition on the occupancy of a vesicle release site. In analogy, a quantity H′(t) can also be introduced which has the conditionality as G′(t) but is the probability density that a spike arrives while a vesicle is absent. Note that F(t) = G′(t) + H′(t). To derive an integral equation for G′(t), involving the previously defined quantities f(t), g(t) and h(t), we consider the four distinct histories that contribute to it. The first involves no intermediate spike between 0 and t, with the first spike arriving at t, and is therefore simply f(t). The other three contributions involve at least one intermediate spike, the penultimate one arriving at s and being integrated over. The first of these three contributions involves the penultimate spike arriving at a time s when a vesicle is present and there being no release, to give G′(s)qf(t − s). The second involves the penultimate spike arriving at a time s when a vesicle is present, there being a release and then the empty site being restocked, thereby contributing G′(s)pg(t − s). The final contribution involves there being no vesicle present at the penultimate spike time s but then being restocked, to give H′(s)g(t − s). Combining these four terms and then simplifying results in the following integral equation

[Image: pcbi.1006232.e056.jpg]G′(t)=f(t)+∫0tdsF(s)g(t-s)+q∫0tdsG′(s)h(t-s).(37)

It will be seen later that this quantity typically only appears with G(t) subtracted from it. Eq (17) includes one of the integrals seen in G′(t) and hence the difference of these quantites obeys

[Image: pcbi.1006232.e057.jpg]G′(t)-G(t)=h(t)+q∫0tds(G′(s)-G(s))h(t-s)(38)

which is an integral equation with a similar structure to Eq (15) for the spike-triggered presynaptic rate F(t) and is simpler to treat analytically and numerically. The cross-covariance in vesicle release is a sum of G′(t) and G(t) weighted by the probabilities that the second site is stocked or unstocked immediately after the first releases (these are given in Eq 36) so that the cross-covariance can ultimately be written in the form

[Image: pcbi.1006232.e058.jpg]CrossCov(χx,χz)=p2r⟨xz⟩∞(δ(t)+⟨x⟩∞⟨xz⟩∞(G(|t|)-r⟨x⟩∞)+q(G′(|t|)-G(|t|)))(39)

where release trains χx and χz are from two distinct release sites that share the same presynaptic cell. To derive the post-synaptic voltage variance the Laplace transform of the difference G′(t) − G(t) is required. From Eq (38) this is

[Image: pcbi.1006232.e059.jpg]LG′(z)-LG(z)=L(z+λ)1-qL(z+λ)soLG′(z)=L(z)(1-L(z+λ)(1-L(z))(1-qL(z+λ)).(40)

For completeness, the Fourier transform is [Image: pcbi.1006232.e060.jpg]G^′(ω)=LG′(iω)+r〈x〉∞πδ(ω).

Postsynaptic voltage variance for multiple release sites

The voltage Eq (31) has a solution analgous to Eq (25) but with a double sum over the N presynaptic neurons and their n contacts. On squaring this and taking expectations there are two groups of non-zero terms amongst the (Nn)² components of the double sum: Nn terms of the auto-covariance (Eq 18) and Nn(n − 1) terms of the cross-covariance (Eq 39). The double integral over the cross-covariance can be performed in the same way as for the auto-covariance, and when combined results in the following form for the post-synaptic voltage variance

[Image: pcbi.1006232.e061.jpg]Var(v)=Nnτa22pr⟨x⟩∞(1+2p(LG(1/τ)-τr⟨x⟩∞))+Nn(n-1)τa22p2r⟨xx′⟩∞(1+2⟨x⟩∞⟨xx′⟩∞(LG(1/τ)-τr⟨x⟩∞)+2q(LG′(1/τ)-LG(1/τ))).(41)

The first line arises from the autocorrelation of χ and the second line from the cross-correlation between χ and χ′. This form can be simplified and the equations for [Image: pcbi.1006232.e062.jpg]ℒG(z) and [Image: pcbi.1006232.e063.jpg]ℒG′(z) used to get

[Image: pcbi.1006232.e064.jpg]Var(v)=τa2Nn2⟨χ⟩(1+2pn(L(1/τ)-L(1/τ+λ)(1-L(1/τ))(1-qL(1/τ+λ))-τr⟨x⟩∞))+τa2Nn2⟨χ⟩(n-1)p⟨xx′⟩∞⟨x⟩∞(1+2qL(1/τ+λ)1-qL(1/τ+λ)).(42)

This formula represents the extension of the analytical forms for the voltage variance to include biophysical details such as: stochastic transmission, quantal effects, short-term depression, multiple contacts, all with non-Poissonian input.

Post-synaptic voltage for gamma-distributed ISIs with multiple contacts

Extending the calculation of the variance of the post-synaptic voltage for gamma-distributed ISIs to the case of multiple contacts is straightforward as it still requires only knowledge of the Laplace transform of the ISI distribution. The temporal form of the cross-covariance, however, requires G′(t) which is given by an integral equation that is awkard to solve numerically. However, the equation for the difference (Eq 38) is easier to analyse and it is the difference that is required for the cross-covariance of the release events between sites sharing a presynaptic neuron. This quantity can be solved in the form of a series in much the same was as was done for F(t) in Eq (22), giving in this case

[Image: pcbi.1006232.e065.jpg]G′(t)-G(t)=e-(αr+λ)tqt∑m=1∞qm(αrt)mαΓ(mα)(43)

which also converges rapidly. Fig 3B shows the correlation (Eq 34) and covariance (Eq 35) between two release sites x and z that share the same presynaptic neuron: note the contrasting behaviour as α scans from bursty to regular at constant presynaptic rate. Fig 3C shows the cross-covariance (Eq 39) for the same three α values used in Fig 2Bii. In Fig 3D the standard-deviation of the voltage and its CV are plotted for three pairs of N and n, such that their product is Nn = 1000. For this choice the voltage mean is the same as in Fig 2D. Note that the standard deviations have a qualitatively different dependence on the regularity of the presynaptic train for different ratios of contact number and presynaptic neuron number: for many contacts the std is highest for bursting presynaptic cells because the effects of the large voltage deviations from multiple release events exceed the effect of mean lower release rate for smaller values of α.

ISIs generated by presynaptic integrate-and-fire neurons

In previous sections analytical forms for the pre-spike occupation 〈x〉∞, autocovariance of the release-train χ and the resulting voltage moments were derived with results illustrated using gamma-distributed ISIs (Eq 3). Integrate-and-fire models, such as the Leaky, Quadratic or Exponential IF models driven by white noise that, following a spike, retain no memory of their previous state will generate spike trains with uncorrelated ISIs. In this case the ISI distribution is identical to the first-passage-time density of the steady-state dynamics. Because of this, all the general results derived thus far are applicable when the presynaptic population is comprised of integrate-and-fire neurons. For the LIF the Fourier transform of the first-passage-time density is available analytically, whereas for non-linear IF models like the QIF or EIF it can be straightforwardly obtained numerically. We now consider examples of these two cases.

Presynaptic leaky integrate-and-fire neurons

We first consider a presynaptic population comprised of LIF neurons that are each driven by an input that has a constant and (independent) fluctuating component. The voltage evolution obeys the equation

[Image: pcbi.1006232.e066.jpg]τv˙=μ-v+σ2τξ(t)(44)

where τ is the membrane time constant, μ is the voltage mean and σ its standard deviation in the absence of a threshold, and ξ(t) is zero-mean 〈ξ(t)〉 = 0, delta-correlated Gaussian white noise such that 〈ξ(t)ξ(t′)〉 = δ(t − t′). With a threshold at vth and reset at vre the neuron will fire at an average rate r, which can be written as a single integral [41]

[Image: pcbi.1006232.e067.jpg]1rτ=∫0∞dyye-y2/2(eyyth-eyyre)(45)

where yth = (vth − μ)/σ with an analogous definition for yre. The Fourier-transform [Image: pcbi.1006232.e068.jpg]f^(ω) of the first-passage-time density, or ISI distribution has been derived in terms of cylinder functions [47] and can also be written as a ratio of integrals [42] of similar form to Eq (45)

[Image: pcbi.1006232.e069.jpg]f^(ω)=∫0∞dyyiωτddy[eyyre-y2/2]∫0∞dyyiωτddy[eyyth-y2/2].(46)

Using the relation (Eq 2) between the expectation of an exponential, and the Fourier and Laplace transform of the ISI distribution, the key quantity required for all the main analytical results derived in this paper is

[Image: pcbi.1006232.e070.jpg]L(z)=⟨e-zt⟩=∫0∞dyyτz-1e-y2/2eyyre∫0∞dyyτz-1e-y2/2eyyth(47)

for ISIs generated by LIF neurons (note that a partial integration step has been performed between Eqs 46 and 47). For example, the mean release-site occupation (Eq 9) just before the arrival of a presynaptic spike takes the form

[Image: pcbi.1006232.e071.jpg]⟨x⟩∞=∫0∞dyyτλ-1e-y22(eyyth-eyyre)∫0∞dyyτλ-1e-y22(eyyth-qeyyre).(48)

for presynaptic LIF neurons. Though not explicitly addressed in this paper, the corresponding formulae for LIF neurons driven by Poissonian shot-noise [42] take a similar form, with 〈x〉∞ again expressible as a ratio of integrals.

Fig 4 examines the release-site occupancy and post-synaptic voltage statistics when the presynaptic spike trains are generated by LIF models. At a fixed rate, a reset near threshold with strong noise will lead to bursting behaviour [43] whereas a low reset with weak noise noise will lead to regular spiking (varying the steady component μ allows the rate to be kept the same). In Fig 4A the presynaptic rates versus μ for three pairs of vth − vre and σ (values given in the figure) ranging from bursting (blue), intermediate (green) and regular (red) are plotted. Fig 4B shows example presynaptic-voltage time courses for these three cases when the presynaptic rate is 10Hz (symbols on Fig 4A). We can simultaneously vary vre and σ linearly, from their values on the blue curve to those on the red curve (with μ compensating so that the rate is constant) to cover the range from bursty to regular firing. The resulting occupancy (Eq 9) and post-synaptic voltage mean, std and CV (Eqs 24 and 29) are plotted in Fig 4C and demonstrate qualitatively similar behaviour to that seen for gamma-distributed ISIs.

Fig 4View larger versionFig 4. Occupancy and post-synaptic voltage statistics from a presynaptic pool of LIF neurons.(A) Presynaptic rate (Eq 45) as a function of drive μ for three pairs of (vth − vre, σ) as marked leading to bursty (blue), intermediate (green) and regular (red) firing statistics. Lines of constant rate are at 5, 10 and 20Hz are marked (dotted lines). (B) Example presynaptic voltage time courses (Eq 44) with spikes marked (black) for the parameters marked with symbols in panel A (μ is varied so they are all at rate 5Hz). (C) The parameters vre and σ were co-varied linearly (with μ compensating so that rates were constant at 5, 10 and 20Hz: see dotted lines in panel A, with same color coding used) and the occupancy (Eq 9), post-synaptic voltage mean (Eq 24), standard deviation (from Eq 29) and CV plotted (against vth − vre). The behaviour seen is qualitatively the same as for the gamma-generated ISIs. In this figure N = 1000, n = 1 and vth = 10mV, with other parameters given in Table 1. The code used to generate this figure is provided in the Supporting Material.

ISI distributions from exponential integrate-and-fire neurons

A number of extensions of the LIF neuron model have been proposed to better capture the non-linearities at the onset of the spike. These include the Quadratic Integrate-Fire model [35] and, more recently, the Exponential Integrate-and-Fire (EIF) model [36] which has been shown to be in good agreement with experimental data [39]. The EIF model takes the form

[Image: pcbi.1006232.e072.jpg]τv˙=μ-v+δTe(v-vT)/δT+σ2τξ(t)(49)

where the parameter δT sets the voltage range over which the spike initiates, vT sets the spike-onset threshold and the other parameters have been previously defined earlier in the context of the LIF model. For the threshold the limit vth → ∞ is typically taken but, practically, as long as vth is significantly above vT the behaviour is insensitive to its exact value. The reset has the same function as for the LIF model. Because of the non-linearity of the model many of the neuronal input-output functions are unavailable in analytical form; however, it is straightforward to solve the underlying differential equations numerically in the steady state [48] and for the first-passage-time density in the Fourier domain [49] (see Section 3.2 of that paper). For the latter calculation, replacing iω with z allows for all of the exponential expectations (Eq 1) required for the analytical results of this paper to be straightforwardly derived. Using this methodology, in Fig 5A–5E the results given in Fig 4 for the LIF are generalised to the EIF model, with the additional complexity of multiple contacts per presynaptic neuron included.

Fig 5View larger versionFig 5. Approximate rate of a postsynaptic EIF neuron driven by a pool of presynaptic EIF neurons, using the matched-variance approximation.(A) Firing rate of presynaptic neurons as a function of μ for three pairs (vT − vre, σ) representing bursting (−3mV, 2.0mV; blue), intermediate (1.0mV, 1.45mV; green) and regular (10mV, 0.2mV; red) firing. (B) Example presynaptic-voltage time courses (Eq 49) corresponding to symbols in panel A. In panels C-F parameters vre and σ were simultaneously varied linearly between the values of the blue and red curves in panel A (μ adjusted so the presynaptic rate remained 10Hz). Hence, vT − vre ranges from bursting −3mV to regular 10mV. (C) Shows the mean pre-spike vesicle occupancy (Eq 9), (D) the mean post-synaptic voltage (Eq 32) and (E) its standard deviation (From Eq 42). In (F) the post-synaptic rate was approximated using the matched-variance approximation, in which the voltage mean and variance (exactly calculated for this filtered drive) are used in the white-noise EIF firing rate calculation. For comparison, the symbols are simulations of the true firing rate. In Panels D-F Nn = 1000 with (N, n) taking the values (1000, 1: grey), (100, 10: yellow), (50, 20: orange) and (25, 40: purple). The reset for the post-synaptic neuron was held constant at vre = 5mV and other parameters are given in Table 1. The code used to generate this figure is provided in the Supporting Material.

From voltage moments to firing rate of IF neurons

It has been shown [44] that the postsynaptic firing rate of EIF neurons is surprisingly insensitive to (positive) temporal correlations when driven by coloured Gaussian noise for a given voltage mean and variance. This allows for a matched-variance approximation to be used in which the firing rate of a coloured-noise driven EIF neuron is approximated by its white noise equivalent but using the same subthreshold mean and variance. Eqs (32) and (42) can be used to find the voltage mean and variance and the EIF rate for white-noise drive is given in [36]. To test whether this approximation has validity, we consider the firing rate of an EIF neuron driven by depressing, stochastic synapses from a presynaptic population of EIF neurons. Following the same approach used in Fig 4 for LIF neurons, we covaried vre and σ to get a range of spiking statistics in the presynaptic EIF population (see Eq 49 for the EIF model definition). In Fig 5A the firing rate of an EIF neuron as a function of the constant drive μ is shown for three different combinations of vT − vre and σ, with example time courses given in Fig 5B each having a rate 10Hz. In Fig 5C–5F the occupancy, voltage statistics and post-synaptic rate are plotted, as vT − vre and σ are simultaneosly linearly varied from their values on the blue and red curves in Fig 5A (with μ adapted so that the presynaptic rate is always 10Hz). A range of connectivity is considered by four combinations of N and n (as marked). As can be seen, the simulations agree well with the post-synaptic firing rate over a range of presynaptic firing patterns and forward connectivity choices. Hence, the matched variance approximation provides a fair account of the firing rate even in cases where the incoming fluctuations are negatively correlated, extending previous results [44].

Generalisation to filtered synapses

It is worth noting it is fairly straightforward to generalise the voltage-variance calculations Eqs (29 and 42) to neurons receiving more biophysically shaped excitatory post-synaptic potentials (EPSPs). For example, if an isolated release event generates an EPSP of the form [Image: pcbi.1006232.e073.jpg]ℰ(t) then, under the assumption of additivity,

[Image: pcbi.1006232.e074.jpg]v=⟨v⟩+∫-∞tdt′E(t-t′)∑i(χi(t′)-⟨χ⟩)(50)

becomes the generalisation of Eq (25) where the voltage mean is

[Image: pcbi.1006232.e075.jpg]⟨v⟩=μ+N⟨χ⟩∫0∞dtE(t).(51)

Following the same approach that led to Eq (29), the voltage variance in this case can be written

[Image: pcbi.1006232.e076.jpg]Var(v)=N⟨χ⟩(∫0∞dtE(t)2)(1+2p∫0∞dz(G(z)-r⟨x⟩∞)∫0∞dtE(t)E(t+z)∫0∞dtE(t)2).(52)

If the form of the EPSP is modelled as a sum of multiple exponentials then the variance can again be expressed by Laplace transforms of G(t). For example, a two-exponential model for an EPSP

[Image: pcbi.1006232.e077.jpg]E(t)=aτ1+τ2τ1-τ2(e-t/τ1-e-t/τ2)(53)

has the integrals

[Image: pcbi.1006232.e078.jpg]∫0∞dtE(t)2=a22(τ1+τ2)and∫0∞dtE(t)E(t+z)∫0∞dtE(t)2=τ1e-z/τ1-τ2e-z/τ2τ1-τ2.(54)

On substitution into the equation for the variance, and after identifying any Laplace transforms of G(t), this gives

[Image: pcbi.1006232.e079.jpg]Var(v)=(τ1+τ2)Na22⟨χ⟩(1+2p(τ1LG(1/τ1)τ1-τ2+τ2LG(1/τ2)τ2-τ1-r⟨x⟩∞(τ1+τ2))).(55)

The forms [Image: pcbi.1006232.e080.jpg]ℒG(z) are provided in the second equation of the pair (20) in terms of the Laplace transform of the first-passage-time density.

Discussion

We have presented a series of novel analytical results that extend the level of biophysical detail incorporated in models of synaptic transmission. Non-Poissonian activity is commonly seen in vivo and can have a substantial effect on neuronal activity; relating this to synaptic dynamics allows for a more comprehensive understanding of the typical behaviour of plastic synapses. We have derived the exact spike-triggered mean vesicle occupancy, noting that it takes a particularly compact form when the input spike train is a renewal process, and highlighted the relationship between the spike-triggered mean and the overall level of vesicle occupancy, confirming the numerical results of Matveev and Wang (2000, [25]), de la Rocha and Parga (2005, [8]) and Reich and Rosenbaum (2013, [33]). We have derived the autocorrelations in vesicle release in terms of integral transforms, confirming the numerical results of Goldman et al (2002, [26]) and extending the analytical results of Goldman (2004, [29]) to account for the biophysically important case of probabalistic vesicle release.

The exact subthreshold voltage variance calculated from the neurotransmitter release autocorrelations is a potentially useful result and incorporates many biophysical details of autocorrelated input spikes, quantal effects, stochastic and cross-correlated vesicle release, and short-term plasticity. The relative effects of these biophysically relevant phenomena can now be quantitatively analysed [27, 50], providing insights into how changes in release-site number n arising from long-term plasticity [51] or in release probability p due to developmental changes [52–54] reshape the postsynaptic response to correlated spike trains.

The results for synaptic transmission when the presynaptic integrate-and-fire neuron is driven by short-time correlated noise processes are another important application, expanding the utility of the model to allow for study of synaptic dynamics alongside other known results concerning the firing rate [41, 42] and correlation structure [18, 43] of such neurons.

Context

Many previous studies have examined how plastic, probabilistic and quantal synapses affect the statistics of patterned transmission through the synapse, and we now provide a selective overview. Vere-Jones (1966) [5] examined a model of a quantal, probabilistic synapse, finding that the release of neurotransmitter is more Poissonian than the afferent activity. Maass and Zador (1999) [27] considered the binary output of a single vesicle release site, investigating how triplets of incoming spikes corresponded to different release patterns under varying synaptic parameters. In particular they showed that dynamic synapses transmit bursts of spikes more reliably than static synapses and have enhanced computational power. Matveev and Wang (2000, [25]) numerically studied the effects of naturalistic presynaptic firing patterns on vesicle release, both with bursts generated by a two-state Markov model and long-time correlated trains. They found that spikes within a burst are suppressed by synaptic depression compared to isolated spikes, and, like Vere-Jones, that neurotransmitter release is more Poissonian than the incoming spikes. Goldman et al (2002, [26]) examined the transmission of doubly-stochastic Poissonian spike trains, constructed to reflect experimentally recorded neuronal bursting, finding again that dynamic synapses decorrelate afferent spike trains and so reduce coding redundancy across a broad range of synaptic parameters. As part of this study, the autocovariance of neurotransmitter release in response to temporally structured spike trains was calculated numerically. Goldman (2004, [29]) derived the information transmission efficiency of a depressing synapse analytically, finding the autocovariance of neurotransmitter release under the assumption of a synapse that reliably releases neurotransmitter whenever a vesicle is present. Using a model comprising a single neurotransmitter release site, de la Rocha et al (2002, [55]) showed that dynamic synapses were more effective transmitters of afferent signals only when the input is non-Poisson, analytically describing the distribution of synaptic release events when the input is a renewal process. They later [8] numerically studied the impact of temporal correlations on synapses containing multiple release sites, showing that bursty stimuli elicited fewer releases of neurotransmitter but that there could be a non-monotonic relationship between presynaptic and postsynaptic firing rates in the presence of input correlations and synaptic dynamics. Fuhrmann et al (2002, [56]) developed the stochastic quantal model used in this paper, capturing the same processes as the continuous phenomenological Tsodyks-Markram model of short-term plasticity [30], but focussed the initial analysis on Poisson spike trains. Ly and Tranchina (2009, [57]) considered numerically the transmission of temporally correlated spike trains across stochastic, but not dynamic, synapses and plotted the autocovariances in vesicle release, as well as the postsynaptic firing rate for renewal process inputs. Rosenbaum et al (2012, [12]) studied information transmission for Poissonian inputs, finding that the incorporation of stochastic quantal effects differentially affected information transmitted at different presynaptic rates. This paper approximated the auto- and cross covariances in neurotransmitter release in response to Poisson drive and these results were shown to be exact by Bird and Richardson (2014, [9]). Reich and Rosenbaum (2013, [33]) studied models of presynaptic spiking both more and less regular than a Poisson process, showing numerically that more regular firing patterns can increase the rate of vesicle release, thereby enhancing the fidelity and efficiency of signal transmission, whilst more irregular spike trains can lead to a decrease in neurotransmitter release. Zhang and Peskin (2015, [50]) developed the results of [12] on information transfer with unreliable dynamic synapses using a slightly simpler model of vesicle recovery, analytically studying the effects of a more general model of presynaptic spiking on neurotransmitter release rates and numerically simulating the effects on the postsynaptic membrane.

Extensions

The matched-variance firing rate approximation in Fig 5 does not constitute a complete framework for treating recurrent networks of neurons with stochastic depressing synapses, because only the rate and not the full ISI statistics of the post-synaptic neuron were derived. However, it does suggest that some approximation scheme that goes beyond the first-order statistics of the ISI distribution might be used to analyse recurrent networks. This is currently a problem of great interest [21, 43, 58] given the strong effects correlated spike trains are acknowledged to have even across static synapses [18]. To include the components of vesicle-release autocorrelation arising from short-term depression, as modelled here, would increase physiological relevance and bring studies of output firing patterns into line with much of the literature on neuronal networks [31, 32, 59].

Another interesting extension is to account for the effect of spike-frequency adaptation currents on presynaptic firing. Adaptation is present across the nervous system [60] and can modulate responses to persistent activity by high-pass filtering and response selectivity [61–63]. These functional roles overlap with those attributed to synaptic depression [64], and there have been a number of recent studies on the interactions of short-term synaptic plasticity with slow adaptation mechanisms [65, 66]. A key feature of adaptation currents is the creation of correlations between interspike intervals [67], generating non-renewal spike trains. These correlations have recently been shown to take the form of a geometric series [68]. Eq (8) presents a way of deriving approximate results for the synaptic transmission for weakly correlated ISIs and suggests a promising avenue of research to go beyond renewal processes and study the effect of these two key short-term adaptive processes in neural circuits.

Supplementary Material

Open In Web Browser

Acknowledgments

We would like to thank Robert Rosenbaum for suggesting the extension to filtered synaptic drive.

Articles from PLoS Computational Biology are provided here courtesy of Public Library of Science

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Verveen A, Derksen H, and Schick K. Voltage fluctuations of neural membrane. Nature. 1967 11; 216: 588–589. doi: 10.1038/216588a0 5583506

White J, Rubinstein J, and Kay A. Channel noise in neurons. Trends Neurosci. 2000 3; 23(3): 131–137. doi: 10.1016/S0166-2236(99)01521-0 10675918

Faisal A, Selen L, and Wolpert D. Noise in the nervous system. Nat Rev Neuro. 2008 4; 9(4): 292–303. doi: 10.1038/nrn2258

Ribrault C, Sekimoto K and Triller A. From the stochasticity of molecular processes to the variability of synaptic transmission. Nat Rev Neuro. 2011 6; 12(7): 375–387. doi: 10.1038/nrn3025

Vere-Jones D. Simple stochastic models for the release of quanta of transmitter from a nerve terminal. Aust J Stat 8:53–63 (1996) doi: 10.1111/j.1467-842X.1966.tb00164.x

Quastel DMJ. The binomial model in fluctuation analysis of quantal neurotransmitter release. Biophys J 72:728–753 (1997) doi: 10.1016/S0006-3495(97)78709-5 9017200

Matveev V and Wang X-J. Implications of all-or-none synaptic transmission and short-term depression beyond vesicle depletion: A computational study. J Neurosci 20:1575–1588 (2000) doi: 10.1523/JNEUROSCI.20-04-01575.2000 10662847

de la Rocha J and Parga N. Short-term synaptic depression causes a non-monotonic response to correlated stimuli. J Neurosci. 2005 9; 25(37): 8416–8431. (2005) doi: 10.1523/JNEUROSCI.0631-05.2005 16162924

Bird A and Richardson M. Long-term plasticity determines the postsynaptic response to correlated afferents with multivesicular short-term synaptic depression. Front Comp Neuro. 2014 1; 8(1)

Stein R. A theoretical analysis of neuronal variability. Biophys J. 1965 3; 5: 173–194 doi: 10.1016/S0006-3495(65)86709-1 14268952

Kuhn A, Aertsen S, and Rotter S. Higher-order statistics of input ensembles and the response of simple model neurons. Neural Comp. 2003 1; 15(1): 67–101 doi: 10.1162/089976603321043702

Rosenbaum R, Rubin J, and Doiron B. Short term synaptic depression imposes a frequency dependent filter on synaptic information transfer. PLoS Comp Bio. 2012 6; 8(6): e1002557 doi: 10.1371/journal.pcbi.1002557

Fellous J, Rudolph M, Destexhe A and Sejnowski T. Synaptic background noise controls the input/output charactersitics of single cells in an in vitro model of in vivo activity. Neuroscience. 2003 8; 122(2): 811–829 doi: 10.1016/j.neuroscience.2003.08.027 14622924

Reyes A. Synchrony-dependent propagation of firing rate in iteratively constructed networks in vitro. Nat Neurosci. 2003 6; 6(6): 593–599 doi: 10.1038/nn1056 12730700

Buzsaki G and Draguhn A. Neuronal oscillations in cortical networks. Science. 2004 6; 304(5679): 1926–1929 doi: 10.1126/science.1099745 15218136

Shinomoto S, Kim H, Shimokawa T, Matsuno N, Funahashi S, Shima K et al Relating neuronal firing patterns to functional differentiation of cerebral cortex. PLoS Comp Bio. 2009 7; 5(7): e1000433 doi: 10.1371/journal.pcbi.1000433

Wilbur W and Rinzel J. A theoretical basis for large coefficient of variation and bimodality in neuronal interspike interval distributions. J Theor Bio. 1983 11; 105(2): 345–368 doi: 10.1016/S0022-5193(83)80013-76656286

Câteau H and Reyes A. Relation between single neuron and population spiking statistics and effects on network activity. Phys Rev Letts. 2006 2; 96(5): e058101

Lindner B. Superposition of many independent spike trains is generally not a Poisson process. Phys Rev E. 2006 2; 73(2): e022901 doi: 10.1103/PhysRevE.73.022901

Pipa G, Grun S, and van Vreeswijk C. Impact of spike-train auto-structure on probability distribution of joint-spike events. Neural Comput. 2013 5; 25(5): 1123–1163 doi: 10.1162/NECO_a_00432 23470124

Dummer B, Wieland S, and Lindner B. Self-consistent determination of the spike-train power spectrum in a neural network with sparse connectivity. Front Comp Neuro. 2014 9; 8: e00104

Eccles J, Katz B and Kuffler S. Nature of the endplate potential in curarized muscle. J Neurophysiol. 1941 1; 4: 362–387 doi: 10.1152/jn.1941.4.5.362

Zucker R and Regehr W. Short-term synaptic plasticity. Ann Rev Physiol. 2002 3; 64: 355–405 doi: 10.1146/annurev.physiol.64.092501.11454711826273

Sudhof T. The synaptic vesicle cycle. Ann Rev Neurosc. 2004 7; 27: 509–547 doi: 10.1146/annurev.neuro.26.041002.131412

Matveev V and Wang X-J. Differential short-term plasticity and transmission of complex spike trains: to depress or to facilitate? Cerebral Cortex 10:1143–1153 (2000) doi: 10.1093/cercor/10.11.1143 11053234

Goldman MS, Maldonado P and Abbott LF. Redundancy reduction and sustained firing with stochastic depressing synapses. J Neurosci 22:584–591 (2002) doi: 10.1523/JNEUROSCI.22-02-00584.2002 11784806

Maass W and Zador A. Dynamic stochastic synapses as computational units. Neural Comp 11:903–917 (1999) doi: 10.1162/089976699300016494

Natschläger T, Maass W and Zador A. Efficient temporal processing with biologically realistic dynamic synapses. Network 12:75–87 (2001)11254083

Goldman MS. Enhancement of information transmission efficiency with unreliable synapses. Neural Computation 16:1137–1162 (2004)15130245

Tsodyks M, and Markram H. The neural code between neocortical pyramidal neurons depends on neurotransmitter release probability. PNAS. 1997 1; 94(2): 719–723 doi: 10.1073/pnas.94.2.719 9012851

Abbott L, Varela J, Sen K, and Nelson S. Synaptic depression and cortical gain control. Science. 275: 220–224 (1997) doi: 10.1126/science.275.5297.221 8985017

Tsodyks M, Pawelzik K and Markram H. Neural networks with dynamic synapses. Neural Comp. 1998 5; 10(4): 821–835 doi: 10.1162/089976698300017502

Reich S and Rosenbaum R. The impact of short term synaptic depression and stochastic vesicle dynamics on neuronal variability. J Comp Neuro. 2013 8; 35(1): 39–51 doi: 10.1007/s10827-012-0438-0

Lapicque L. Recherches quantitatives sur l’excitation électrique des nerfs traiteé comme une polarization. J Physiol Pathol Gen. 1907; 9:620–635

Ermentrout GB and Kopell N. Parabolic bursting in an excitable system coupled with a slow oscillation. SIAM J. Applied Mathematics 1986; 46: 233–253 doi: 10.1137/0146017

Fourcaud-Trocmé N, Hansel D, van Vreeswijk C, and Brunel N. How spike generation mechanisms determine the neuronal response to fluctuating inputs. J Neurosci. 2003 12; 23(37): 11628–11640 doi: 10.1523/JNEUROSCI.23-37-11628.2003 14684865

Gerstein G and Mandelbrot B. Random walk models for the spike activity of a single neuron. Biophys J. 1964 1; 4: 41–68 doi: 10.1016/S0006-3495(64)86768-0 14104072

Jolivet R, Lewis T, and Gerstner W. Generalized integrate-and-fire models of neuronal activity approximate spike trains of a detailed model to a high degree of accuracy. J Neurophysiol. 2004 8; 92(2): 959–976 doi: 10.1152/jn.00190.2004 15277599

Badel L, Lefort S, Brette R, Petersen CCH, Gerstner W, and Richardson MJE. Dynamic I-V curves are reliable predictors of naturalistic pyramidal-neuron voltage traces. J. Neurophysiol. 99: 656–666 (2008) doi: 10.1152/jn.01107.2007 18057107

Franklin J and Bair W. The effect of a refractory period on the power spectrum of neuronal discharge. SIAM Journal of Applied Mathematics 55: 1074–1093 (1995) doi: 10.1137/S0036139993258691

Brunel N and Hakim V. Fast global oscillations in networks of integrate-and-fire neurons with low firing rates. Neural Comp. 1999 10; 11(7): 1621–1671 doi: 10.1162/089976699300016179

Richardson M and Swarbrick R. Firing-rate response of a neuron receiving excitatory and inhibitory synaptic shot noise. Phys Rev Letts 105: e178102 (2010) doi: 10.1103/PhysRevLett.105.178102

Ostojic S. Interspike interval distributions of spiking neurons driven by fluctuating inputs. J Neurophysiol. 2011 7; 106(1): 361–373. doi: 10.1152/jn.00830.2010 21525364

Alijani A and Richardson M. Rate-response of neurons subject to fast or frozen noise: from stochastic and homogeneous to deterministic and heterogeneous populations. Phys Rev E. 2011 7; 84(1): e011919 doi: 10.1103/PhysRevE.84.011919

Renart A, Moreno-Bote R, Wang XJ, and Parga N Mean-driven and fluctuation driven persistent activity in recurrent networks. Neural Comp. 2007 1; 19(1):1–46. doi: 10.1162/neco.2007.19.1.1

Loebel A, Silberberg G, Helbig D, Markram H, Tsodyks M, and Richardson M. Multiquantal release underlies the distribution of synaptic efficacies in the neocortex. Front Comp Neuro. 2009 11; 3(27).

Lindner B, Schimansky-Geier L and Longtin A. Maximizing spike train coherence or incoherence in the leaky integrate-and-fire model. Phys. Rev. E 66 031916 (2002) doi: 10.1103/PhysRevE.66.031916

Richardson MJE. Firing-rate response of linear and nonlinear integrate-and-fire neurons to modulated current-based and conductance-based synaptic drive. Phys Rev E 76: article-no 021919 (2007) doi: 10.1103/PhysRevE.76.021919

Richardson MJE. Spike-train spectra and network response functions for non-linear integrate-and-fire neurons. Biological Cybernetics 99: 381–392 (2008) doi: 10.1007/s00422-008-0244-y 19011926

Zhang C and Peskin CS. Improved signaling as a result of randomness in synaptic vesicle release. PNAS 112: 14954–14959 (2015) doi: 10.1073/pnas.1513160112 26627245

Loebel A, Le Be J, Richardson M, Markram H, and Herz A. Matched pre- and post-synaptic changes underlie synaptic plasticity over long time scales. J Neurosci. 2013 4; 33(15): 6257–6266. doi: 10.1523/JNEUROSCI.3740-12.2013 23575825

Reyes A and Sakmann B. Developmental switch in the short-term modification of unitary EPSPs evoked in Layer 2/3 and Layer 5 pyramidal neurons of rat neocortex. J Neurosci. 1999 5; 19(10): 3827–3835. doi: 10.1523/JNEUROSCI.19-10-03827.1999 10234015

Frick A, Feldmeyer D, and Sakmann B. Postnatal development of synaptic transmission in local networks of L5A pyramidal neurons in rat somatosensory cortex. J Physiol. 2007 11; 585(1): 103–116. doi: 10.1113/jphysiol.2007.141788 17916610

Kerr M, Wall M, and Richardson M. Adenosine a1-receptor activation mediates the developmental shift at Layer-5 pyramidal-cell synapses and is a determinant of mature synaptic strength. J Physiol. 2013 7; 591(13): 3371–3380. doi: 10.1113/jphysiol.2012.244392 23613526

de la Rocha J, Nevado A, and Parga N. Information transmission by stochastic synapses with short-term depression: Neural coding and optimization. Neurocomputing. 44-46:85–90. (2002) doi: 10.1016/S0925-2312(02)00362-4

Fuhrmann G, Segev I, Markram H, and Tsodyks M. Coding of temporal information by activity-dependent synapses. J Neurophysiol. 87:140–148. (2002) doi: 10.1152/jn.00258.2001 11784736

Ly C and Tranchina D. Spike train statistics and dynamics with synaptic input from any renewal process: A population density approach. Neural Comp 21:360–396 (2009) doi: 10.1162/neco.2008.03-08-743

Schwalger T, Droste F, and Lindner B. Statistical structure of neural spiking under non-Poissonian or other non-white stimulation. J Comp Neurosci. 2015 5; 39: 29–51. doi: 10.1007/s10827-015-0560-x

Lindner B, Gangloff D, Longtin A and Lewis J. Broadband coding with dynamic synapses. J Neurosci. 2009 2; 29(7): 2076–2088. doi: 10.1523/JNEUROSCI.3702-08.2009 19228961

Benda J and Herz A. A universal model of spike frequency adaptation. Neural Comp. 2003 11; 15(11): 2523–2564. doi: 10.1162/089976603322385063

Benda J, Longtin A, and Maler L. Spike-frequency adaptation separates transient communication signals from background oscillations. J Neurosci. 2005 3; 25(9): 2312–2321. doi: 10.1523/JNEUROSCI.4795-04.2005 15745957

Peron S and Gabbiani F. Role of spike-frequency adaptation in shaping neuronal response to dynamic stimuli. Biol Cybern. 2009 6; 100(6): 505–520. doi: 10.1007/s00422-009-0304-y 19381681

Wang X. Calcium coding and adaptive temporal computation in cortical pyramidal neurons. J Neurophysiol. 1998 3; 79(3): 1549–1566. doi: 10.1152/jn.1998.79.3.1549 9497431

Abbott L and Regehr W. Synaptic computation. Nature. 2004 10; 431: 796–803. doi: 10.1038/nature03010 15483601

Moezzi B, Iannella N, and McDonnell M. Modelling the influence of short term depression in vesicle release and stochastic calcium channel gating on auditory nerve spontaneous firing statistics. Front Comp Neuro. 2014 12; 8(163).

Nagel K, Hong E, and Wilson R. Synaptic and circuit mechanisms promoting broadband transmission of olfactory stimulus dynamics. Nat Neurosci. 2014 12; 18: 56–65. doi: 10.1038/nn.3895 25485755

Schwalger T, Fisch K, Benda J, and Lindner B. How noisy adaptation of neurons shapes interspike interval histograms and correlations. PLoS Comp Bio. 2010 12; 6(12): e1001026 doi: 10.1371/journal.pcbi.1001026

Schwalger T and Lindner B. Patterns of interval correlations in neural oscillators with adaptation. Front Comp Neuro. 2013 11; 7(164).

[Back]

Table 1. Parameters and their values used for figures.

Parameters are grouped together as: those pertaining to occupancy and release statistics; those required for calculation of the post-synaptic response; and those used to generated ISI statistics. Certain neuronal parameters, such as τ and μ, are used for both pre and post-synaptic integrate-fire models.

Parameter | Interpretation | Value
p | probability of release on arrival of a spike, given vesicle is present | 0.6
λ | rate empty release sites are restocked | 2Hz
r | rate of presynaptic spiking | varies
τ | membrane time constant | 20ms
μ | voltage resting potential in absence of input | varies
a | amplitude of EPSP induced by fusion and release of single vesicle | 0.3mV
N | number of presynaptic neurons | varies
n | vesicle release sites per presynaptic neuron | varies
α | shape parameter for gamma-generated ISIs | varies
σ | voltage standard deviation for presynaptic integrate-and-fire neurons | varies
v re | voltage reset after an action potential | varies
v th | voltage threshold at which a spike is registered | 10mV (LIF) 15mV (EIF)
δ T | spike onset range for EIF model | 1.5mV
v T | spike onset threshold for EIF model | 10mV

DOI: 10.1371/journal.pcbi.1006232.t001

[Back]

[Back]

Fig 1. Vesicle occupancy and neurotransmitter release rate are a function of the ISI distribution.

(A) Presynaptic spike train, vesicle occupancy and neurotransmitter release time-course for gamma-distributed ISIs with α as marked (in panel B) for the same release probabilty p = 0.6, presynaptic rate r = 5Hz and postsynaptic rate λ = 2Hz. (B) Corresponding ISI distributions for the three α values. (C) The cumulative distribution. (D) The mean 〈x〉 (grey) and the pre-spike 〈x〉∞ (black) release-site occupancies (Eqs 9 and 10). Note that 〈x〉∞ increases with increasing regularity, larger α. (E) Variance of the two occupancies exhibiting non-monotonic behaviour. (F) The mean neurotransmitter release rate 〈χ〉 is directly proportional to 〈x〉∞ and so shares its qualitative dependence on α. The code used to generate this figure is provided in the Supporting Material.

[Back]

[Back]

Fig 2. Temporal and spectral statistics of the presynaptic spike-train and neurotransmitter release required for post-synaptic voltage mean and variance.

(A) The spike-train (Ai) and neurotransmitter release (Aii) event-triggered rate for three cases ranging from bursty to regular with α as marked (see Eq 22 and surrounding text). (B) The equivalent autocovariances of the spike train (Bi) and synaptic release events (Bii). For bursty presynaptic trains (small α = 0.4) the autocovariance is positive; however all cases have negative autocovariances for the release events themselves (Eq 18) because of the filtering by the depressing synapses. (C) The corresponding power spectra (Eq 21) of the spike train (Ci) and synaptic release events (Cii). (D) The postsynaptic voltage mean, standard deviation (std) and coefficient of variation std/mean (Eqs 24 and 29). Though the mean and std both increase with increasing presynaptic regularity, mirroring the behaviour of 〈x〉∞ (Fig 1D), the voltage CV itself decreases with increasing regularity. For panels A-C parameters and color are same as Fig 1, and in panel D parameter α is varied with other parameters provided in Table 1. The code used to generate this figure is provided in the Supporting Material.

[Back]

[Back]

Fig 3. Correlated occupancy, release and post-synaptic voltage statistics when there are multiple contacts per presynaptic cell.

(A) Illustration of a case with one presynaptic cell N = 1 making three contacts n = 3. Shown in descending order: presynaptic spike train (α = 0.4, r = 5Hz), occupancy of the three release sites (restock λ = 2Hz), the corresponding release-event time courses (p = 0.6), and the voltage time course (see Table 1 for parameters). (B) The occupancy correlation 〈xz〉∞ for a pair of sites receiving the same presynaptic spike train and the covariance [Image: pcbi.1006232.e049]〈xz〉∞-〈x〉∞2 (see Eqs 34 and 35). (C) The cross-covariance of release events from a pair of sites sharing the same presynaptic cell (Eq 39). (D) Post-synaptic voltage standard deviation (top: from Eq 42) and CV (bottom) as a function of α for four combinations of N, n as marked. With each of these choices, where Nn = 1000, the mean voltage (Eq 32) is the same as in Fig 2D. The code used to generate this figure is provided in the Supporting Material.

[Back]

[Back]

Fig 4. Occupancy and post-synaptic voltage statistics from a presynaptic pool of LIF neurons.

(A) Presynaptic rate (Eq 45) as a function of drive μ for three pairs of (vth − vre, σ) as marked leading to bursty (blue), intermediate (green) and regular (red) firing statistics. Lines of constant rate are at 5, 10 and 20Hz are marked (dotted lines). (B) Example presynaptic voltage time courses (Eq 44) with spikes marked (black) for the parameters marked with symbols in panel A (μ is varied so they are all at rate 5Hz). (C) The parameters vre and σ were co-varied linearly (with μ compensating so that rates were constant at 5, 10 and 20Hz: see dotted lines in panel A, with same color coding used) and the occupancy (Eq 9), post-synaptic voltage mean (Eq 24), standard deviation (from Eq 29) and CV plotted (against vth − vre). The behaviour seen is qualitatively the same as for the gamma-generated ISIs. In this figure N = 1000, n = 1 and vth = 10mV, with other parameters given in Table 1. The code used to generate this figure is provided in the Supporting Material.

[Back]

[Back]

Fig 5. Approximate rate of a postsynaptic EIF neuron driven by a pool of presynaptic EIF neurons, using the matched-variance approximation.

(A) Firing rate of presynaptic neurons as a function of μ for three pairs (vT − vre, σ) representing bursting (−3mV, 2.0mV; blue), intermediate (1.0mV, 1.45mV; green) and regular (10mV, 0.2mV; red) firing. (B) Example presynaptic-voltage time courses (Eq 49) corresponding to symbols in panel A. In panels C-F parameters vre and σ were simultaneously varied linearly between the values of the blue and red curves in panel A (μ adjusted so the presynaptic rate remained 10Hz). Hence, vT − vre ranges from bursting −3mV to regular 10mV. (C) Shows the mean pre-spike vesicle occupancy (Eq 9), (D) the mean post-synaptic voltage (Eq 32) and (E) its standard deviation (From Eq 42). In (F) the post-synaptic rate was approximated using the matched-variance approximation, in which the voltage mean and variance (exactly calculated for this filtered drive) are used in the white-noise EIF firing rate calculation. For comparison, the symbols are simulations of the true firing rate. In Panels D-F Nn = 1000 with (N, n) taking the values (1000, 1: grey), (100, 10: yellow), (50, 20: orange) and (25, 40: purple). The reset for the post-synaptic neuron was held constant at vre = 5mV and other parameters are given in Table 1. The code used to generate this figure is provided in the Supporting Material.

[Back]