Dendritic Voltage Recordings Explain Paradoxical Synaptic Plasticity: A Modeling Study

Experiments have shown that the same stimulation pattern that causes Long-Term Potentiation in proximal synapses, will induce Long-Term Depression in distal ones. In order to understand these, and other, surprising observations we use a phenomenological model of Hebbian plasticity at the location of the synapse. Our model describes the Hebbian condition of joint activity of pre- and postsynaptic neurons in a compact form as the interaction of the glutamate trace left by a presynaptic spike with the time course of the postsynaptic voltage. Instead of simulating the voltage, we test the model using experimentally recorded dendritic voltage traces in hippocampus and neocortex. We find that the time course of the voltage in the neighborhood of a stimulated synapse is a reliable predictor of whether a stimulated synapse undergoes potentiation, depression, or no change. Our computational model can explain the existence of different -at first glance seemingly paradoxical- outcomes of synaptic potentiation and depression experiments depending on the dendritic location of the synapse and the frequency or timing of the stimulation.

Introduction

How are memories encoded in the brain? In 1949, Donald Hebb postulated that a synapse connecting two neurons strengthens if both neurons are active together (Hebb, 1949). Numerous experiments have confirmed the interaction of pre- and postsynaptic neuronal activity during the induction of synaptic plasticity (see Levy and Stewart, 1983; Bliss and Collingridge, 1993; Sjöström et al., 2001; Wang et al., 2005; Markram et al., 2011). The critical postsynaptic signal for plasticity induction might be related to voltage (Artola et al., 1990), calcium (Cormier et al., 2001), or backpropagating action potentials (Markram et al., 1997b). If changes in subthreshold voltage or calcium concentration are the critical signals on the postsynaptic side, then plasticity does not require the postsynaptic neuron to fire a somatic spike. On the other hand, if backpropagating action potentials are critical, then synaptic plasticity outcomes can be completely described by the relative timing of pre- and postsynaptic spikes in the form of a Spike-Timing Dependent Plasticity rule (STDP rule, Abbott and Nelson, 2000).

A first, and fundamental, challenge for all STDP models is the existence of subthreshold plasticity in the absence of somatic spikes (Ngezahayo et al., 2000; Golding et al., 2002; Lisman and Spruston, 2005; Brandalise and Gerber, 2014). Another challenge for some (Gerstner et al., 1996; Kistler and van Hemmen, 2000; Song et al., 2000), but not all (Senn et al., 2001; Pfister and Gerstner, 2006; Clopath and Gerstner, 2010; Graupner and Brunel, 2012) STDP models is the interaction of frequency and spike-timing dependence so that long-term potentiation (LTP) for pre-before-post timing disappears at low frequencies (Markram et al., 1997b; Sjöström et al., 2001) whereas long-term-depression (LTD) for post-before-pre timing does not (Sjöström et al., 2001). Finally, an important finding that challenges classical models of STDP (Gerstner et al., 1996; Song et al., 2000; Van Rossum et al., 2000) is the observation that plasticity rules depend on synapse location: whereas normally a protocol of presynaptic stimulation followed by postsynaptic activity induces potentiation, it was found to induce depression in distal synapses (Froemke et al., 2005; Letzkus et al., 2006). We refer to the above challenges as paradoxical effects of STDP and ask whether a single phenomenological model can account for all of these.

The observations that plasticity depends on dendritic synapse location and does not require somatic spikes hint at dendritic effects that are not accounted for by standard STDP models. Indeed, dendritic spikes have been shown to play a key role for the induction of plasticity in various brain regions (Holthoff et al., 2004; Kampa et al., 2006, 2007; Remy and Spruston, 2007; Gambino et al., 2014). Dendritic events are linked to active channel properties which can vary along the dendritic tree (Spruston, 2008). Local dendritic non-linearities could thus explain why different learning rules can be obtained with similar protocols in different brain regions or even within the same cell as a function of synapse location.

However, it is difficult to translate such an insight into a concrete biophysical model because it would require as a starting point a valid, and broadly accepted, model of local dendritic non-linearities as well as a biophysically plausible synaptic plasticity model – but neither of these are readily available. While a first step in this direction has been taken recently (Ebner et al., 2019), most of the biophysical and phenomenological plasticity rules proposed over the years have in practice been tested using simplified point neuron models. And even if a biophysically detailed non-linear dendrite model with active zones (Hay et al., 2011) were to be used, the location, and composition of ion channels in such active zones, might not be exactly the one encountered in the specific neuron recorded in an experiment.

In this paper we propose a voltage-based plasticity model that can be fitted to experiments without the need of fine-tuning any biophysical neuron model. Our model can be seen as a variation of earlier phenomenological voltage-based (Brader et al., 2007; Clopath et al., 2010) and calcium-based plasticity models (Shouval et al., 2002; Rubin et al., 2005; Graupner and Brunel, 2012). Our model has a set of plasticity parameters that need to be tuned. However, tuning of additional neuronal parameters is not necessary simply because we do not use any biophysical neuron model but work directly with the experimentally measured time course of the voltage in the neighborhood of the synapse.

In this paper we investigate whether such a phenomenological model of synaptic plasticity in which the arrival of neurotransmitter is paired with the postsynaptic voltage at the location of the synapse can explain the aforementioned paradoxical experimental results on excitatory synapses: (i) LTP in the absence of somatic spikes; (ii) the interaction of spike-timing and spike frequency; and (iii) inversion of a plasticity rule as a function of dendritic location. We focus on three experiments where the time course of dendritic voltage was measured during the application of a plasticity-inducing protocol (in neocortex, Letzkus et al., 2006 and in hippocampus, Brandalise and Gerber, 2014; Brandalise et al., 2016). We show that the time-course of the postsynaptic voltage in the neighborhood of the synapse, in combination with presynaptic signaling, is a reliable predictor of synaptic plasticity and is sufficient to explain the outcome of the experiments. In brief, our voltage-based model replicates plasticity behaviors of synapses across various dendritic locations in neocortex and hippocampus.

Results

Voltage Dependence of Plasticity

Our model combines ideas from phenomenological models of voltage-based plasticity (Brader et al., 2007; Clopath et al., 2010) with the ‘veto’ concept of Rubin et al. (2005). As described in the Methods section, each presynaptic spike leaves, in our model, a trace $\overline{x}$ at the synapse; analogously, the activity of the postsynaptic neuron also leaves two traces at the synapse, described as two low-pass filtered version $\overline{\text{u}}$₊ and $\overline{\text{u}}$_– of the dendritic voltage u. Potentiation can occur if the variable $\overline{\text{u}}$₊ (i.e., the voltage filtered with time constant τ₊) is above some threshold θ₊. Similarly, depression can occur if $\overline{\text{u}}$_– (i.e., the voltage filtered with time constant τ_–) is above some threshold θ_–. In both cases, the amount of change depends on the momentary value of the trace $\overline{x}$ left by a presynaptic spike (Figures 1A–C). Importantly, to translate the competition between the molecular actors involved in LTP and LTD (phosphatase vs. kinase, see Bhalla and Iyengar, 1999; Xia and Storm, 2005; Herring and Nicoll, 2016) into mathematical equations, we introduce into our model a ‘veto’ concept: a potentiation signal overwrites LTD that would occur otherwise (Cho et al., 2001; O’Connor et al., 2005; Rubin et al., 2005). In our model, the veto mechanism is implemented by a dynamic LTP-dependent increase of the LTD-threshold θ_– that is characterized by parameters b_θ and τ_θ.

FIGURE 1

Figure 1. Voltage-dependent plasticity model. (A) The activity X(t) of the presynaptic neuron induces local dendritic voltage changes u(t) in the postsynaptic neuron. The change in synaptic strength (LTP or LTD induction) depends on the timing of the presynaptic spike, and the voltage close to the synapse. (B) The presynaptic spike X(t) leaves a trace $\overline{x}$(t) at the synapse. The voltage u is low-pass filtered with a time constant τ₊ (for the variable $\overline{\text{u}}$₊) or τ_– (for $\overline{\text{u}}$_–). The amount of LTP is proportional to $\overline{x}$ multiplied by $\overline{\text{u}}$₊, while $\overline{\text{u}}$₊ is above a threshold θ₊. (C) Similarly, the amount of LTD is proportional to $\overline{x}$ multiplied by $\overline{\text{u}}$_–, while $\overline{\text{u}}$_– is above a threshold θ_– which is lower than θ₊ and increases when LTP occurs. (D–G) Plasticity in hippocampal model cells: extracellular afferent stimulation is paired with voltage-clamp of the postsynaptic neuron at different potentials (see Ngezahayo et al., 2000). 100 brief extracellular afferent stimulations are done at 2 different frequencies: 2 (full line) and 40 Hz (dotted line). (E,F) Synaptic strength w in percentage of its initial value as a function of voltage with respect to resting potential, in mV, for two sets of parameters. (G) The presynaptic trace $\overline{x}$ (blue) and the voltage u superimposed with its filtered versions $\overline{\text{u}}$₊ and $\overline{\text{u}}$_– (black) for 3 different values of clamped voltage (8, 20 and 30 mV) and different stimulation frequencies. The thresholds are indicated by the dashed green (θ₊) and dashed purple (θ_–) lines. The parameters are: (E,G) τ_x = 5ms, τ₊ = 6ms, τ_– = 15ms, θ₊ = 10mV, θ₀ = 5mV, A_LTP = 0.0001mV^{– 1}.ms^{– 1}, A_LTD = 0.0001mV^{– 1}.ms^{– 1}, b_θ = 31000mV.ms, τ_θ = 14ms. (F) τ_x = 5ms, τ₊ = 7ms, τ_– = 15ms, θ₊ = 13mV, θ₀ = 7mV, A_LTP = 0.0001mV^{– 1}.ms^{– 1}, A_LTD = 0.0001mV^{– 1}.ms^{– 1}, b_θ = 45000mV.ms, τ_θ = 5ms.

If we pair presynaptic stimulations with a constant voltage at the location of the synapse, our model shows three regimes (Figures 1D–G): (i) for hyperpolarization or voltage close to rest, synapses do not show any plasticity; (ii) for voltages above a first threshold θ₀, presynaptic stimulation leads to a depression of the synapses; (iii) for voltages above a second threshold θ₁ the synapses exhibit potentiation. Depending on the parameters, the voltage-plasticity relationship can be linear (Figure 1E) or non-linear (Figure 1F). Our model is consistent with experimental results of Ngezahayo et al. (2000) who paired 2 Hz presynaptic stimulations with constant postsynaptic depolarizations (voltage clamp) and determined the stationary voltage dependence of LTP and LTD induction. For a wide range of parameter choices, our model qualitatively reproduces this voltage dependence, if we assume that during clamping the voltage u at the dendrite is equal to the somatic voltage.

In more realistic experiments, the voltage at the location of the synapse is not constant but changes as a function of time. In the following, we investigate if our model can reproduce the experimentally measured plasticity observed with various LTP- or LTD-inducing protocols. We focus on experimental paradigms where the dendritic voltage was recorded close to the stimulated synapse during plasticity induction (Letzkus et al., 2006; Brandalise and Gerber, 2014; Brandalise et al., 2016). In brief, the methodological approach is as follows: we feed dendritic voltage time courses and corresponding presynaptic spike trains into our plasticity model. The variable u(t) is thus the experimentally recorded voltage, and not a simulated one. Parameters of our model are then fitted so that the plasticity predicted by our model matches as closely as possible the experimental plasticity values. To this end, we use an optimization algorithm which minimizes the mean-squared error (squared difference between theoretical and experimental plasticity outcomes, see Methods).

Subthreshold Plasticity in the Hippocampus

We first used experimental data from rat hippocampus where Brandalise and Gerber, 2014 investigated plasticity at the CA3 recurrent synapses using a subthreshold protocol. In this protocol, an excitatory postsynaptic potential (EPSP) of a few millivolts induced by stimulation of the CA3 recurrent pathway was paired with a subthreshold mossy fiber (MF) stimulation that, if stimulated separately, also led to an EPSP of a few millivolts (Figure 2A). Dendritic voltage recordings (data kindly shared by F. Brandalise) were used as the variable u(t) in our plasticity model (see Materials and Methods).

FIGURE 2

Figure 2. Subthreshold and spike-timing dependent plasticity at CA3 synapses in the hippocampus. (A) Experimental setup. Stimulation of CA3 recurrent inputs (blue electrode) was paired with a subthreshold stimulation of mossy fiber inputs (MF, brown electrode). The pairing is repeated 60 times at 0.1 Hz (Brandalise and Gerber, 2014). (B–D) Experimental voltage traces (black, middle panel) caused by subthreshold stimulations with a 10ms (B), 0 ms (C), or –40 ms (D) interval. Blue traces (top panels, $\overline{x}$(t)), green and purple (middle panels: $\overline{\text{u}}$₊ green full line; $\overline{\text{u}}$_– purple full line; θ₊ green dashed line; and θ_– purple dashed line), and black (bottom panels, w) show the time course of selected model variables during the simulated experiments. Two different cases are illustrated in (B), because for the same stimulation two types of voltage responses were recorded in the dendrite (black electrode in A): linear (left) and supralinear (right) ones. The supralinear responses correspond to the occurrence of dendritic spikes. (E) STDP protocol (Brandalise et al., 2016): stimulations of CA3 recurrent inputs were paired 50 times with brief somatic current injections (2 ms; 4 nA) which evoked action potentials (APs). (F) When 3 APs were evoked at a frequency of 200 Hz, dendritic spikes (black, right panels) occurred in 60% of the trials. In the remaining 40%, a linear response was generated (left); color of other traces as in B and C. (G) Responses were always linear, if only one AP was paired with stimulation of CA3 recurrent inputs. (B–D,F,G) Time of presynaptic stimulation is set to 0. (H) Plasticity outcome using different plasticity protocols. STDP protocol (see E): presynaptic stimulation paired either with 3 APs (see F), with 3 APs and concomitant application of a hyperpolarizing pulse, with 3 APs generated at a lower frequency (50 Hz) or with 1 AP (see G). Subthreshold protocol: CA3 stimulation only, no pairing; pairings with different time intervals (0 ms, 10 ms, –40 ms, see B–D) and, for the 10 ms time interval, simultaneously blocking the occurrence of supralinear events (10ms, block). Filled circles represent data from individual cells. Black error bars represent experimental mean ± SD. Red crosses represent simulations using the parameters obtained with the best fit. All data points in (H) are fitted with a single set of parameters (see text and Table 1). Differences in plasticity for the same protocol (e.g. 10 ms) arise due to differences in experimental voltage traces. Data shared by F. Brandalise. (I) Squared Error (SE) of the best fit subtracted from the SE obtained after increasing (upper panel) or decreasing (lower panel) each parameter by 5%, one parameter at a time (ΔSE).

In agreement with the experimental results of Brandalise and Gerber (2014), no plasticity was induced in our model synapses when there was no MF stimulation (CA3 alone) or when the MF and CA3 stimulations occurred at the same time (0 ms). LTD was observed when the MF stimulation preceded the CA3 stimulation with a 40 ms interval (noted as −40 ms in Figures 2D,H). LTP was observed when the MF stimulation followed the CA3 stimulation with a 10 ms time interval (+10 ms).

FIGURE 3

Figure 3. Distance-dependent STDP at synapses between layer 2/3 and layer 5 pyramidal neurons in somatosensory cortex. (A–C) Voltage traces at distal (middle) and proximal (bottom) synapses. Postsynaptic bursts (3 action potentials, APs at 200 Hz) are paired with presynaptic action potentials (± 10 ms time interval, pairing frequency of 1 Hz, 150 repetitions). Experimental voltage traces u (black line) are redrawn from Letzkus et al. (2006). At distal synapses, dendritic spikes are generated. Glutamate trace $\overline{x}$ (blue line, top), filtered versions of the voltage $\overline{\text{u}}$₊ and $\overline{\text{u}}$_– (solid green and purple lines), and synaptic weight w (number: value of w after one pairing; dashed line: b_θ is set to 0, no veto) as a function of time. (A) During + 10 ms pairings, the value of the presynaptic trace $\overline{x}$ has already decreased significantly before $\overline{\text{u}}$₊ reaches θ₊. The amount of LTP is not high enough for the veto to have a significant impact on LTD induction. (B) In contrast, for –10 ms pairings, $\overline{x}$ switches from 0 to its maximal value 1 at around the time when $\overline{\text{u}}$₊ reaches its maximal value far above θ₊ and $\overline{\text{u}}$_– is slightly above θ_–. Therefore, the amount of LTP induced is high and significantly reduces LTD via an increase of the LTD threshold θ_–. (C) However, in the presence of NiCl₂ (blocks a subtype of voltage-gated calcium channels), the difference between $\overline{\text{u}}$₊ and θ₊ is significantly reduced compared to the control case in B, leading to no change in the synaptic strength. (D,E) Experimental and predicted change in synaptic weight. Crosses and dots or crosses with error bars represent plasticity from Letzkus et al. (2006) and red lines simulations. (D) Plasticity along the dendrite for the protocols described in A-C: 1pre-3 APs (left) or 3 APs-1 pre (right). EPSP rise time at the soma is a proxy of the distance between the plastic synapse and the soma (a distance of 110 μm, 330 μm and 660 μm correspond to a rise time of 1.8 ms, 2.5 ms and 3.5 ms, respectively, see Letzkus et al., 2006). See text for more details. (E) Plasticity at distal synapses in the presence of NiCl₂.

FIGURE 4

Figure 4. Pairing and timing-dependence of plasticity at neocortical synapses. (A) Two synaptically connected L5 neurons were stimulated with different time intervals (–10, 0, 10 and 25 ms) at different pairing repetition frequencies: 0.1 Hz, 10 Hz, 20 Hz and 40 Hz. (B) Simulated dendritic (black), somatic (orange) and experimentally recorded somatic (blue) voltage time course for + 10 ms time interval (number: peak value). The experimental voltage trace is redrawn from Sjöström et al. (2001); inset: EPSP time course. (C) Plasticity as a function of spike timing. Each panel represents one pairing repetition frequency. LTP is induced at high frequencies. Black errorbars represent data from Sjöström et al. (2001) and red squares represent our plasticity model. (D) Presynaptic trace $\overline{x}$ (blue), voltage u (black) and its filtered versions $\overline{\text{u}}$₊ (green) and $\overline{\text{u}}$_– (purple) for + 10 ms time interval: 10 Hz (left) or 40 Hz (right) repetition frequency.

Importantly, in the experiments of Brandalise and Gerber (2014), repetitive pairings at +10 ms were found to induce two types of dendritic voltage trajectories (Figure 2B): either a nearly linear addition of the EPSPs caused by MF and CA3 stimulation or (in 31 ± 4% of cases) a strongly supralinear voltage response (Brandalise and Gerber, 2014). For example, a given cell could exhibit during one paired stimulation a linear response, and the same cell could then show in the next stimulation a supralinear response. Since our model uses the experimental voltage time course, we could predict the plasticity outcome for a given cell (Figure 2H) based on the dendritic voltage recorded for that specific cell across different conditions (CA3 alone, 0 ms and +10 ms protocols, see Materials and Methods). Hence variations in the amount of plasticity are explained in the model by differences in the voltage recordings – without any tuning of parameters between different cells. In particular, when the occurrence of supralinear events was completely blocked, no plasticity was observed in our model synapse (10 ms, block, Figure 2H), in agreement with experiments (Brandalise and Gerber, 2014).

To understand how the model works, let us focus on a few examples (Figures 2B–D). During the −40 ms protocol, the low-pass filtered voltage trace $\overline{\text{u}}$₊ did not reach the threshold θ₊ for LTP induction, whereas the voltage $\overline{\text{u}}$_– filtered with a larger time constant reached θ_–, inducing LTD (Figure 2D). With the +10 ms protocol, $\overline{\text{u}}$₊ reached θ₊ only during trials in which a supralinear event occurred. During linear events, $\overline{\text{u}}$₊ and $\overline{\text{u}}$_– did not reach their respective thresholds θ₊ and θ_– (Figure 2B). This was also the case during the 0 ms protocol (Figure 2C).

STDP Protocol in the Hippocampus and Cross-Validation Procedure

In a further set of experiments using a burst STDP protocol, Brandalise et al. (2016) paired a recurrent CA3 EPSP with 3 action potentials (APs) at 200 Hz (10 ms time interval, see Figure 2E). This stimulation led in more than half of the trials to the generation of a dendritic spike (Figure 2F), unless a hyperpolarizing step current was applied in the dendrite during the brief somatic injections triggering the APs. Similarly, pairing the CA3 EPSP with 3 APs at 50 Hz or with a single AP did not generate a dendritic spike (Figure 2G). Representative dendritic voltage traces u(t) measured by Brandalise et al. (2016) were used in our plasticity model which was able to reproduce the plasticity outcomes of the experiments (Figure 2H).

We emphasize that our plasticity model with a fixed set of parameters (Table 1) could reproduce the outcome of all the STDP experiments as well as that of all the earlier subthreshold protocols (Figure 2H). The set of parameters in Table 1 was obtained with all available voltage traces corresponding to 15 plasticity outcomes.

TABLE 1

Table 1. Parameters minimizing the error (see section “Materials and Methods”).

Since our model has 9 free parameters, the question arises whether the model is overfitting the available data points or whether it would correctly generalize to novel data. In order to check the model’s predictive power, we used an additional, independent, optimization procedure (leave-one-out cross-validation): we fitted the model parameters on plasticity outcomes for 14 voltage traces by minimizing the mean-squared error and predicted the plasticity outcome on the remaining trace (see Table 2 for the statistics over all 15 leave-one-out experiments). Even though the median error after testing the plasticity outcome on the excluded voltage traces was (as expected) larger than the median training error (Table 2), its actual value of 1.5 ^∗ 10^–3 was comparable to the normalized error of 6.2 ^∗ 10^–4 observed in the direct fitting approach of Table 1. Furthermore, we found that most parameter values are consistent across the 15 leave-one-out experiments as indicated by a small standard deviation of the parameter value compared to its mean value (Table 2); exceptions were the veto parameters b_θ and τ_θ which showed rather large standard deviations. A sensitivity analysis further confirmed that the exact values of these two parameters was not critical (Figure 2I and Supplementary Figure 1). Thus, cross-validation and sensitivity analysis confirm that the model has predictive power.

TABLE 2

Table 2. Cross-validation results.

Location-Dependent Plasticity in Neocortical Apical Dendrites

We next tested our model on data recorded at synapses between layer 2/3 and layer 5 pyramidal neurons in slices from rat somatosensory cortex (Letzkus et al., 2006). We fed our plasticity model with representative dendritic voltage traces recorded close to synapses located 100 μm (proximal), 330 μm or 660 μm (distal) away from the soma (Figures 3A–C; Letzkus et al., 2006). Using the fitting procedure described above, we found that our voltage-based plasticity model with a single set of parameters could account for the plasticity results obtained with a burst STDP protocol (see Figure 3D and Table 1). Consistent with experimental results, the observed plasticity results varied in the model depending on synapse location on the dendritic tree. Proximal EPSPs were potentiated during pairings with somatic bursts of APs occurring 10 ms after the onset of the EPSP (+10 ms) and depressed when the postsynaptic bursts occured 10 ms before the EPSP (−10 ms); at distal synapses, however, the pattern was reversed and EPSPs depressed during +10 ms pairings and were potentiated for −10 ms pairings (Letzkus et al., 2006).

Moreover, the model with the same set of parameters could also explain why distal EPSPs no longer potentiated after pairings at −10ms but still depressed during pairings at +10 ms, if the amplitude of the dendritic spikes evoked by the AP bursts decreased due to the presence of NiCl₂ (a blocker of a subtype of voltage-gated calcium channels, Figures 3C,E, and Letzkus et al., 2006).

To understand the workings of our model, we observed different model variables as a function of time. At distal synapses, during +10 ms pairings, the value of the presynaptic trace $\overline{x}$ had already decreased significantly when $\overline{\text{u}}$₊ reached the threshold θ₊ (Figure 3A). The amount of LTP was therefore not high enough for the veto to have a significant impact on LTD induction. As a result, weak LTD occurs (Figure 3D). In contrast, for −10 ms pairings, $\overline{x}$ switched from 0 to its maximal value 1 at a moment when $\overline{\text{u}}$₊ was close to its maximal value well above θ₊ (Figure 3B). Therefore, the amount of LTP induced was high. The large LTP signal vetoed the induction of LTD as manifested by an increase in the LTD threshold θ_–. As a result, LTP dominates, in agreement with experiments (Figure 3D). However, for −10 ms pairings, in the presence of NiCl₂ or at proximal synapses, the difference between $\overline{\text{u}}$₊ and θ₊ was significantly reduced compared to what was observed at distal synapses, leading to an absence of synaptic potentiation (Figures 3B,C). Moreover, blocking of the veto-mechanism reduces the quality of the fit (Table 3).

TABLE 3

Table 3. Parameters minimizing the error with and without the veto term.

Thus, the results of our voltage-based plasticity model support the idea that differences in the voltage traces can explain the spatial differences in the learning rule, as suggested by Letzkus et al. (2006). Importantly, all the above experiments (Letzkus et al., 2006) are explained by the same model with the same set of parameters (Table 1).

High-Frequency Pairings in Neocortical Basal Dendrites

We have until now focused on plasticity results obtained after repeated pairings of pre and postsynaptic activities at a low frequency (0.1 and 1 Hz). Yet, an important feature of synaptic plasticity is its frequency-dependence. Different amounts of plasticity are obtained by repeating the same pairings at different frequencies. Unfortunately, experimental dendritic recordings do not exist for these types of experiments. Results have been obtained among others at L5-L5 synapses of rat neocortical neurons (Sjöström et al., 2001), which are well-characterized. More than half of the synaptic contacts between L5 neurons are made on basal dendrites (80 ± 35 μm from the soma in young rats, Markram et al., 1997a), for which a detailed model exists (Nevian et al., 2007). We simulated dendritic voltage around 80 μm from the soma using the aforementioned neuron model at four different frequencies (0.1, 10, 20 Hz and 40 Hz) and with different time intervals between the presynaptic and the postsynaptic stimulation (−10, 0, +10 and +25 ms). As shown in Figure 4, our model with fixed set of parameters can reproduce both the frequency-dependence and spike-timing dependence of plasticity.

At low frequencies, the time between two pairings is long enough so that the membrane potential u repolarizes back to its resting value. As a consequence, $\overline{\text{u}}$₊ is close to zero when the next pairing occurs. This is not the case at high frequencies: the residual depolarization between postsynaptic spikes allows $\overline{\text{u}}$₊ to reach the threshold θ_+, leading to LTP induction (see Figure 4D). Similarly, a correlation between the amount of residual depolarization and the amount of LTP has been found in Sjöström et al. (2001). Interestingly, the spike after-depolarization of the neuron model by Nevian et al. (2007) seems to have a shorter time constant than the one recorded in Sjöström et al. (2001), see Figure 4B. This explains why in our model no LTP is induced for pairing frequencies below 20Hz. At high frequencies, the veto mechanism tunes down LTD (see Figure 4D). Blocking of the veto-mechanism reduces the quality of the fit (Table 3).

Since we have a validated model of the L5 basal dendrites, we can predict the plasticity outcome for plasticity protocols with triplets of spikes at L5-L5 synapses. During triplet experiments, a presynaptic [postsynaptic] spike is triggered between the occurrence of two postsynaptic [presynaptic] spikes. As demonstrated experimentally by Wang et al. (2005), the outcome of triplet experiments is not a linear sum of the outcome of each pair of pre- and post-synaptic spikes taken separately. Among others, while a post-pre pair with a 5 ms interval can trigger LTD, addition of a postsynaptic spike 5 ms after the presynaptic spike leads to LTP (see Table 4). Note that the data of Wang et al. (2005) is from hippocampal neurons and that our model predicts that plasticity behaves somewhat differently in basal dendrites of L5 neurons compared to hippocampus.

TABLE 4

Table 4. Predicted plasticity for pairs and triplets of spikes.

To summarize, the same voltage-based plasticity model can account for three different series of experiments corresponding to four publications (Sjöström et al., 2001; Letzkus et al., 2006; Brandalise and Gerber, 2014; Brandalise et al., 2016). Importantly, the model parameters are slightly different for different synapse types, but each series of experiments from one synapse type is explained by a single set of model parameters (Table 1). In other words, model parameters are kept fixed across all experimental results in a given experimental preparation.

Discussion

Long-term potentiation or long-term depression are induced through the combined action of the presynaptic and postsynaptic activities. We showed that a single phenomenological voltage-based model could explain results using various synaptic plasticity protocols: experiments (i) with voltage clamp (Figure 1); (ii) with variable time interval between presynaptic and postsynaptic spikes (Figures 2–4); (iii) with variable pairing frequency (Figure 4); (iv) with multiple postsynaptic spikes (Figures 2–4); (v) with subthreshold plasticity (Figure 2) and (vi) with location-dependence (Figure 3).

Comparison With Other Plasticity Models

The model proposed here, as well as other voltage-based and calcium based models (Abarbanel et al., 2002; Shouval et al., 2002; Rubin et al., 2005; Brader et al., 2007; Clopath and Gerstner, 2010; Graupner and Brunel, 2012), is a phenomenological one since it does not aim to describe the full mechanistic signaling chain from presynaptic spike arrival to a change in the number of AMPA receptors or presynaptic release probabilities. Rather it should be considered as a ‘black-box’ model that summarizes a large range of experimental results in the form of a compressed ‘learning rule’ with only a few variables.

In calcium-based models (Shouval et al., 2002; Rubin et al., 2005; Graupner and Brunel, 2012), calcium concentration acts as a summary variable that includes effects of both pre- and postsynaptic activity (influx of calcium through NMDA channels and other voltage-gated ion channels e.g.). The level (Shouval et al., 2002) or time course (Rubin et al., 2005; Graupner and Brunel, 2012) of the calcium concentration in the simulated model is then compared with threshold variables in order to predict occurrence of LTP or LTD. Rubin et al. (2005) proposed to add a “veto” of LTP on LTD when a relatively high calcium threshold was reached which inspired the veto mechanism in the present model.

In voltage-based models (Brader et al., 2007; Clopath and Gerstner, 2010), and similarly in our model, presynaptic activity leaves a filtered trace at the synapse. This trace can be interpreted for excitatory synapses as the amount of glutamate bound to postsynaptic receptors which induces calcium influx into the neuron (and/or, in the case of LTD, conformational changes of the receptors initiating intracellular signaling). It is this glutamate trace that interacts either with the postsynaptic voltage directly or with a low-pass filtered version thereof. The comparison of the voltage variables with several thresholds allows to predict the induction of LTP or LTD of those synapses that have been presynaptically stimulated (Brader et al., 2007; Clopath and Gerstner, 2010). Thus voltage-based models jump over the biophysics of calcium dynamics and connect the presynaptic stimulation in combination with the time course of the postsynaptic voltage directly with the outcome of plasticity experiments. We note that LTD could have a different mechanism (e.g., non-ionotropic, Nabavi et al., 2013) compared to LTP (ionotropic), while still consistent with our phenomenological model.

Our model and the model of Clopath et al. (2010) are different in some respects.

First, in Clopath et al. (2010), LTD is triggered by the joint action of an instantaneous spike event, rather than the glutamate trace and postsynaptic voltage. Using a glutamate trace (as opposed to a presynaptic spike event that covers a much shorter moment in time) is in our hands the only way to make LTD possible for pre-before-post pairings at low pairing frequencies (see Letzkus et al., 2006 as an example). Also, an extended glutamate trace looks biologically more plausible than a “point-like event” assumed in some classic STDP models (Kistler and van Hemmen, 2000; Song et al., 2000).

Furthermore, in Clopath et al. (2010), LTP can occur only if two conditions are met: the momentary voltage u(t) and the low-pass filtered voltage $\overline{\text{u}}$₊ need to be above a threshold θ₊ and θ_–, respectively. In other words, the membrane must already be depolarized before a spike occurs (see Sjöström et al., 2001). The two conditions together imply a quadratic dependence on voltage in the LTP inducing term (Clopath and Gerstner, 2010). Instead of a quadratic voltage term for LTP induction, our model works with a linear dependence on the thresholded, low-pass filtered voltage in combination with a veto-mechanism similar to the one suggested by Rubin et al. (2005).

Previous models were able to quantitatively fit the frequency dependence of STDP (Senn et al., 2001) as well as triplet and quadruplet effects of STDP protocols (Pfister and Gerstner, 2006; Cai et al., 2007; Clopath and Gerstner, 2010; Graupner and Brunel, 2012). The model of Graupner and Brunel (2012) also indicated how changes of STDP rules as a function of synaptic location on the dendrite could be qualitatively accounted for by changes of model parameters; in the absence of dendritic recordings and an appropriate dendrite model, a quantitative fit was not to be expected. A recent study by Ebner et al. (2019) used a detailed multicompartmental model of a neocortical neuron to simulate postsynaptic voltage at different dendritic locations and combined this voltage with a phenomenological model for four pathways of LTP and LTD induction. However, our voltage-based model is probably the first one to directly link dendritic voltage recordings with plasticity outcome, bypassing the need for a biophysically correct dendrite model.

In order to stabilize plasticity, the papers of Pfister and Gerstner (2006) and Clopath et al. (2010) suggest a form of metaplasticity implemented by a slow adjustment of the coefficients of LTD driven by the average mean firing rate. For network simulations we suggest to either implement an analogous slow adjustment of our LTD term as in Clopath et al. (2010) or keep parameters fixed and replace metaplasticity by heterosynaptic plasiticy (Zenke et al., 2015; Zenke and Gerstner, 2017). The Clopath model also uses hard bounds for the weights which is recommended for network simulations.

Role of Dendritic Spikes

Letzkus et al. (2006) showed that at distal locations, the peak amplitude of isolated backpropagating action potentials was half the size than that at proximal locations. Furthermore, postsynaptic bursts at the soma generated dendritic calcium spikes at distal locations. The two observations suggest that the somatic spike is less important for plasticity in distal dendrites than localized depolarizations at the location of the synapse. Similarly, in the hippocampal experiments of Brandalise and Gerber (2014) and Brandalise et al. (2016), dendritic NMDA spikes were generated: they resulted from high frequency bursting during the STDP protocol, and from broad and long mossy-fiber evoked EPSPs during the subthreshold protocol.

Both Letzkus et al. (2006) and Brandalise and Gerber (2014); Brandalise et al. (2016) showed that LTP was abolished when dendritic spikes were blocked (pharmacologically or by hyperpolarizing the cell). In our model, the voltage time course at the location of the synapse determines whether or not LTP (or LTD) is induced at stimulated synapses. If the low-pass filtered voltage $\overline{\text{u}}$₊ does not reach a threshold θ₊, then potentiation is impossible. Intrinsic dendritic non-linearities of the postsynaptic neuron can boost voltage and explain the existence of different – at a first glance seemingly paradoxical - outcomes of plasticity experiments (Mishra et al., 2016). Since we paste the experimentally measured voltage traces directly into our plasticity model, the biophysical source of the depolarization does not matter.

Predictions

Since our model is a phenomenological one (as opposed to a biophysical model that attempts to describe the full signal induction chain, e.g., Bhalla and Iyengar, 1999; Lisman and Zhabotinsky, 2001; Badoual et al., 2006; Saudargiene and Graham, 2015), it cannot be used as a predictive tool in cases where specific biochemical molecules are manipulated without affecting the voltage time course. However, one interesting qualitative prediction follows from the interaction of the veto-concept in our voltage based model. We predict a voltage-dependence of LTP induction (Figures 1E,F) that depends on the stimulation frequency of glutamate pulses. Since presynaptic vesicles are likely to deplete rapidly, we propose an experiment where presynaptic spike arrivals are replaced by glutamate puffs of standardized size while the postsynaptic voltage is clamped at a constant voltage. The prediction from our simple voltage-based model is that the voltage dependence of LTP induction becomes steeper at higher stimulation frequencies – even if the number of pulses is kept constant.

A second prediction concerns the shape of the dendritic voltage time course. Suppose that via dendritic voltage clamp, we artificially impose the postsynaptic voltage to follow a square-wave of amplitude Δu and duration T. The plasticity behavior will depend on both Δu and T, as shown on Figure 5. When the presynaptic neuron spikes in the middle of a square pulse of duration T = 5 ms, the amount of LTP induced increases with Δu. For T = 15 ms, plasticity will follow an ‘inverted u-shape’ as a function of voltage amplitude. If now the presynaptic spike is delivered 10 ms after the end of the square pulse, the synapse undergoes LTD.

FIGURE 5

Figure 5. Non-linear voltage-dependence. (A) The dendritic voltage is clamped for a fixed duration T and varying amplitudes Δu. The resulting squared voltage pulse is paired with a presynaptic spike X arriving 10 ms before the start of the pulse (dashed blue), 10 ms after the end of the pulse (dash-dotted blue), or in the center of the pulse (full blue). (B) Plasticity as a function of voltage amplitude Δu for T = 5 ms (left panels), T = 15 ms (middle panels) or T = 25 ms (right panels), using two sets of parameters (Letzkus: 100 or 10 pairings, top and middle panels respectively, or Brandalise: 10 pairings, bottom panels, see Table 1).

Conclusion

We do not claim that elevated voltage in combination with neurotransmitter release is the direct cause of induction of LTP or LTD. Rather our philosophy is that the voltage time course, if experimentally available, is a very good indicator of whether or not synaptic changes are induced in those synapses that have been presynaptically stimulated. In other words, our model describes the Hebbian condition of joint activity of pre- and postsynaptic neuron in a compact form as the interaction of the glutamate trace left by a presynaptic spike with the time course of the postsynaptic voltage. This philosophy does not exclude that a pharmacological block of later steps in the signaling chain could interrupt the LTP/LTD induction or that a direct experimental manipulation of postsynaptic calcium could induce synaptic plasticity in the absence of presynaptic spike arrival or postsynaptic depolarization. Rather our intuition is that, under physiological conditions, the time course of the voltage in the neighborhood of a stimulated synapse is a reliable indicator of the likelihood of that synapse to undergo plasticity. Our leave-one-out cross-validation results (Table 2) show that this intuition can be transformed into a working model to predict the outcome of future plasticity induction experiments given the voltage trace.

Materials and Methods

Voltage-Based Model of Synaptic Plasticity

The plasticity model (Figure 1) is a combination of earlier voltage-based models (Brader et al., 2007; Clopath et al., 2010) and the veto concept of Rubin et al. (2005).

Plastic changes of a synapse are caused by potentiation (LTP) or depression (LTD) of the synaptic weight w and add up to a total weight change

$$\frac{d}{dt}w\left(t\right)=\frac{d}{dt}{w}_{\text{LTP}}-\frac{d}{dt}{w}_{\text{LTD}}.$$

Potentiation or depression of the weight is induced by a Hebbian combination of presynaptic and postsynaptic activity. Postsynaptic activity is represented by the (low-pass filtered) voltage at the location of the synapse. Presynaptic activity is represented by the spike train X(t) (a sequence of Dirac delta-pulses) arriving at the synapse. The spike train is low-pass filtered and gives rise to a ‘trace’

$${\mathrm{\tau }}_{\text{x}}\frac{d}{dt}\overline{x}\left(t\right)=-\overline{x}\left(t\right)+X(t)$$

where $\overline{x}$ can be thought of as the amount of neurotransmitter bound to the postsynaptic receptors. The value of $\overline{x}$ increases at the arrival of a spike and decays exponentially with a time constant τ_x during the interval between spike arrivals (see Figure 1B).

Depression (LTD) is induced if a low-pass filtered version $\overline{u}$_– of the postsynaptic voltage is above a threshold θ_–and the “trace of presynaptic activity” $\overline{x}$ does not tend to zero,

$$\frac{d}{dt}{w}_{\text{LTD}}\left(t\right)=A_{\text{LTD}}\overline{x}\left(t\right){\left[{\overline{u}}_{-}\left(t\right)-{\mathrm{\theta }}_{-}(t)\right]}_{+}$$

where $\overline{u}$_– is defined as

$${\mathrm{\tau }}_{-}\frac{d}{dt}{\overline{u}}_{-}\left(t\right)=-{\overline{u}}_{-}\left(t\right)+u(t)$$

with time constant τ_–. The amplitude parameter A_LTD characterizes the magnitude of LTD. [y₊ equals y if y > 0, 0 otherwise.

Potentiation (LTP) is induced if another low-pass filtered version $\overline{u}$₊ of the voltage is above a threshold θ₊ and the “trace of presynaptic activity” $\overline{x}$ does not tend to zero,

$$\frac{d}{dt}{w}_{\text{LTP}}\left(t\right)=A_{\text{LTP}}\overline{x}\left(t\right){\left[{\overline{u}}_{+}\left(t\right)-{\mathrm{\theta }}_{+}\right]}_{+}$$

where $\overline{u}$₊ is defined as

$${\mathrm{\tau }}_{+}\frac{d}{dt}{\overline{u}}_{+}\left(t\right)=-{\overline{u}}_{+}\left(t\right)+u(t)$$

with time constant τ₊. The amplitude parameter A_LTP characterizes the magnitude of LTP.

Finally, depression and potentiation compete. If potentiation occurs, the threshold θ_– increases. The value of θ_– is determined by the following equation:

$${\mathrm{\theta }}_{-}(t)={\mathrm{\theta }}_0+\mathrm{\theta }(t)$$

with a fixed part θ₀ and a variable part θ(t) that follows the equation

$${\mathrm{\tau }}_{\mathrm{\theta }}\frac{d\mathrm{\theta }}{dt}=-\mathrm{\theta }+b_{\mathrm{\theta }}\frac{d}{dt}{w}_{\text{LTP}}$$

with time constant τ_θ and interaction parameter b_θ. This interaction of LTD and LTP parallels the ‘veto’ concept of Rubin et al. (2005).

We assume that the plasticity framework defined by the above set of equations is generic for glutamatergic NMDA synapses whereas the specific choice of parameters for amplitudes, thresholds and time constants depends on the specific neuron and synapse type as well as on temperature and ion concentrations in the bath of the experimental slice preparation.

Postsynaptic Voltage Trace

In the above plasticity model, the value of the postsynaptic voltage at the location of the synapse plays a crucial role. We have access to three experimental datasets where voltage has been measured at a dendritic location close to the synapse (Letzkus et al., 2006; Brandalise and Gerber, 2014; Brandalise et al., 2016 and data kindly shared by the authors). Thus, for these plasticity experiments, we do not need to use a neuron model to generate voltage traces; rather, we directly insert a representative experimental voltage trace u(t) into the equations of our plasticity model.

From the Brandalise and Gerber (2014) dataset, we had access to three cells which had the mean amplitude of their supralinear events, calculated in relation to linear event amplitude increased by at least a factor of two (i.e., difference 100% or more, Table 5, column 3). This indicates that the dendritic recording electrode was close enough to the stimulated CA3 recurrent synapses to pick up such a non-linear effect (see Brandalise and Gerber, 2014). For these cells, the voltage time course combined with plasticity outcome was available for 3 different protocols, which were performed in the following order:

TABLE 5

Table 5. Characteristics of the +10ms pairing protocol.

1. no MF stimulation (CA3 alone).

2. the MF and CA3 stimulations occurred at the same time (0 ms).

3. the MF stimulation followed the CA3 stimulation with a 10 ms time interval (+10 ms). The percentage of supralinear events that occurred during this protocol for the 3 individual cells is given in Table 5.

For the remaining protocols (10 ms block, −40 ms and STDP), we used representative voltage time courses and averaged plasticity values, as dendritic recordings and plasticity measurements were done in two different set of cells.

We also model results from Sjöström et al. (2001). In this case, we only had access to representative voltage traces measured at the soma. Since we need for our plasticity model voltage traces in the neighborhood of a synapse, we used the model of L5 basal dendrites from Nevian et al. (2007), available on ModelDB (#124394) to mimic dendritic voltage traces. The multicompartmental model was simulated in NEURON. Action potentials were generated by a 5 ms step current of 3 nA in the somatic compartment and backpropagated through Hodgkin-Huxley-like sodium and potassium channels located on the soma and dendrite. EPSPs were generated by an EPSP-like current injection (double exponential current: 0.5 ms rise time constant, 10 ms decay time constant and peak conductance of 0.1 nS).

The resting potential of all voltage traces (experimental ones and simulation-based ones) has been shifted to 0. This shift allows us to counteract any discrepancies in absolute voltage arising from the electrophysiological recording system or from differences in resting membrane potential across different brain regions and neuron types.

Parameter Optimization

Our model only defines a mathematical framework whereas specific parameter values may depend on neuron type, synapse type, brain region, as well as details of slice preparations. Therefore, we use different sets of parameters, depending on the experiments we want to model. We take (experimental or simulated) voltage traces as input to our model. Differential equations were solved using forward Euler and with an integration time step of 0.1 ms. Synaptic weights w were initialized at w_i = 0.5 and at the end of the simulation we read out the final value w_f.

The 9 parameters of our model were fitted to the outcome of different experiments using the Matlab function fmincon (interior-point algorithm). We fixed θ₊ > θ₀ and defined some upper and lower bounds for the parameters (see Table 6). Time constants are in milliseconds with lower bounds always at 2 ms and upper bounds below 100 ms. In order to mitigate the problem of local minima, we used 25 predefined combinations of parameters as initial points for the optimization algorithm (all inside the bounds). We calculated the least squared error (LSE), which minimizes the quantity SE

TABLE 6

Table 6. Lower and upper bound used during the fmincon search (same units as in Table 1).

$$SE=\sum _{pp}{({\left[\frac{w_f-w_i}{w_i}\right]}_{pp}-\mathrm{\Delta }{\text{w}}_{pp}^{exp})}^2$$

where $\mathrm{\Delta }{\text{w}}_{pp}^{exp}$ is the experimental plasticity value measured during the protocol pp. Since we are interested in the optimal set of parameters, we report in the paper always the parameters from the optimization run which yielded the smallest LSE. We checked that an automatic generation of initial points did not alter the results (Matlab function GlobalSearch).

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://github.com/clairemb90/Voltage-based-model.

Author Contributions

CM-B, WG, and LL designed research. CM-B performed research and analyzed data. MCT contributed unpublished analytic tools. CM-B and WG wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the Swiss National Science Foundation (no. 200020_184615) and by the European Union Horizon 2020 Framework Program under grant agreement no. 785907 (Human Brain Project, SGA2).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Federico Brandalise and Friedemann Zenke for careful reading and critical comments on the manuscript. We also thank Federico Brandalise and Johannes Letzkus for providing additional voltage traces. This manuscript has been released as a pre-print at https://arxiv.org/abs/2001.03614, Meissner-Bernard et al. (2020).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnsyn.2020.585539/full#supplementary-material

Supplementary Figure 1 | Variation of the squared error (SE) as a function of parameter change. (A) SE (vertical axis) when all the parameters (see Table 1) were increased or decreased (horizontal axis) by a fixed percentage. (B) Only the parameter τ_x is changed (horizontal axis) by a fixed percentage. (C) Filled contour plot of the SE while 2 parameters are increased or decreased by a given percentage: τ₊ & θ₊ (C1), τ_– & θ₀ (C2), A_LTP & A_LTD (C3) b_θ & τ_θ (C4). (D) Plasticity value in 9 different conditions (see Figure 2). Black circles and error bars represent experimental data. Red crosses represent simulations using the parameters obtained with the best fit. Gray symbols in D represent simulations using the parameters obtained with the best fit except a few which were changed by a certain percentage or when all parameters were changed by a fixed percentage (compare symbols in A,C1–C3): hexagon in (A) (−2%), upwards-pointing triangle in (C1) (−5 and +4%), rectangle in (C2) (+3 and −4%), cross in (C3) (+7 and +7%).

References

Abarbanel, H., Huerta, R., and Rabinovich, M. (2002). Dynamical model of long-term synaptic plasticity. Proc. Natl. Acad. Sci. U S A. 99, 10132–10137. doi: 10.1073/pnas.132651299

PubMed Abstract | CrossRef Full Text | Google Scholar

Abbott, L. F., and Nelson, S. B. (2000). Synaptic plasticity: taming the beast. Nature 3, 1178–1183. doi: 10.1038/81453

PubMed Abstract | CrossRef Full Text | Google Scholar

Artola, A., Bröcher, S., and Singer, W. (1990). Different voltage dependent thresholds for inducing long-term depression and long-term potentiation in slices of rat visual cortex. Nature 347, 69–72. doi: 10.1038/347069a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Badoual, M., Zou, Q., Davison, A. P., Rudolph, M., Bal, T., Fregnac, Y., et al. (2006). Biophysical and phenomenological models of multiple spike interactions in spike-timing dependent plasticity. Int. J. Neural. Systems 16, 79–97. doi: 10.1142/s0129065706000524

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhalla, U. S., and Iyengar, R. (1999). Emergent properties of networks of biological signaling pathways. Science 283, 381–387. doi: 10.1126/science.283.5400.381

PubMed Abstract | CrossRef Full Text | Google Scholar

Bliss, T. V. P., and Collingridge, G. L. (1993). A synaptic model of memory: long-term potentiation in the hippocampus. Nature 361, 31–39. doi: 10.1038/361031a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Brader, J., Senn, W., and Fusi, S. (2007). Learning real-world stimuli in a neural network with spike-driven synaptic dynamics. Neural. Comput. 19, 2881–2912. doi: 10.1162/neco.2007.19.11.2881

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandalise, F., Carta, S., Helmchen, F., Lisman, J., and Gerber, U. (2016). Dendritic NMDA spikes are necessary for timing-dependent associative LTP in CA3 pyramidal cells. Nat. Comm. 7:13480.

Google Scholar

Brandalise, F., and Gerber, U. (2014). Mossy fiber-evoked subthreshold responses induce timing-dependent plasticity at hippocampal CA3 recurrent synapses. PNAS 111, 4303–4308. doi: 10.1073/pnas.1317667111

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, Y., Gavornik, J. P., Cooper, L. N., Yeung, L. C., and Shouval, H. Z. (2007). Effect of stochastic synaptic, and dendritic dynamics on synaptic plasticity in visual cortex, and hippocampus. J. Neurophysiol. 97, 375–386. doi: 10.1152/jn.00895.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, K., Aggleton, J. P., Brown, M. W., and Bashir, Z. I. (2001). An experimental test of the role of postsynaptic calcium levels in determining synaptic strength using perirhinal cortex of rat. J. Physiol. 532(pt 2), 459–466. doi: 10.1111/j.1469-7793.2001.0459f.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Clopath, C., Busing, L., Vasilaki, E., and Gerstner, W. (2010). Connectivity reflects coding: a model of voltage based spike-timing-dependent plasticity with homeostasis. Nat. Neurosci. 13, 344–352. doi: 10.1038/nn.2479

PubMed Abstract | CrossRef Full Text | Google Scholar

Clopath, C., and Gerstner, W. (2010). Voltage and spike timing interact in STDP a unified model. Front. Synaptic Neurosci. 2:25. doi: 10.3389/fnsyn.2010.00025

PubMed Abstract | CrossRef Full Text | Google Scholar

Cormier, R. J., Greenwood, A. C., and Connor, J. A. (2001). Bidirectional synaptic plasticity correlated with the magnitude of dendritic calcium transients above a threshold. J. Neurophysiol. 85, 399–406. doi: 10.1152/jn.2001.85.1.399

PubMed Abstract | CrossRef Full Text | Google Scholar

Ebner, C., Jedlicka, P., and Cuntz, H. (2019). Unifying long-term plasticity rules for excitatory synapses by modeling dendrites of cortical pyramidal neurons. Cell Rep. 29, 4295–4307. doi: 10.1016/j.celrep.2019.11.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Froemke, R. C., Poo, M. M., and Dan, Y. (2005). Spike-timing-dependent synaptic plasticity depends on dendritic location. Nature 434, 221–225. doi: 10.1038/nature03366

PubMed Abstract | CrossRef Full Text | Google Scholar

Gambino, F., Pages, S., Kehayas, V., Baptista, D., Tatti, R., Carleton, A., et al. (2014). Sensory-evoked LTP driven by dendritic plateau potentials in vivo. Nature 515, 116–119. doi: 10.1038/nature13664

CrossRef Full Text | Google Scholar

Gerstner, W., Kempter, R., van Hemmen, J. L., and Wagner, H. (1996). A neuronal learning rule for sub-millisecond temporal coding. Nature 386, 76–78. doi: 10.1038/383076a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Golding, N. L., Staff, N. P., and Spruston, N. (2002). Dendritic spikes as a mechanism for cooperative long-term potentiation. Nature 418, 326–331. doi: 10.1038/nature00854

PubMed Abstract | CrossRef Full Text | Google Scholar

Graupner, M., and Brunel, N. (2012). Calcium-based plasticity model explains sensitivity of synaptic changes to spike pattern, rate, and dendritic location. Proc. Natl. Acad. Sci. 109, 3991–3996. doi: 10.1073/pnas.1109359109

PubMed Abstract | CrossRef Full Text | Google Scholar

Hay, E., Hill, S., Schürmann, F., Markram, H., and Segev, I. (2011). Models of neocortical layer 5b pyramidal cells capturing a wide range of dendritic and perisomatic active properties. PLoS Comput. Biol. 7:e1002107. doi: 10.1371/journal.pcbi.1002107

PubMed Abstract | CrossRef Full Text | Google Scholar

Hebb, D. O. (1949). The Organisation of Behavior. New York, NY: Wiley.

Google Scholar

Herring, B. E., and Nicoll, R. A. (2016). Long-term potentiation: from CaMKII to AMPA receptor trafficking. Ann. Rev. Physiol. 78, 351–365. doi: 10.1146/annurev-physiol-021014-071753

PubMed Abstract | CrossRef Full Text | Google Scholar

Holthoff, K., Kovalchuk, Y., Yuste, R., and Konnerth, A. (2004). Single-shock LTD by local dendritic spikes in pyramidal neurons of mouse visual cortex: dendritic spikes and LTD. J. Physiol. 560, 27–36. doi: 10.1113/jphysiol.2004.072678

PubMed Abstract | CrossRef Full Text | Google Scholar

Kampa, B., Letzkus, J., and Stuart, G. (2006). Requirement of dendritic calcium spikes for induction of spike-timing-dependent synaptic plasticity. J. Physiol. 574. 1, 283–290. doi: 10.1113/jphysiol.2006.111062

PubMed Abstract | CrossRef Full Text | Google Scholar

Kampa, B. M., Letzkus, J. J., and Stuart, G. J. (2007). Dendritic mechanisms controlling spike-timing dependent synaptic plasticity. Trends Neurosci. 30:456. doi: 10.1016/j.tins.2007.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kistler, W. M., and van Hemmen, J. L. (2000). Modeling synaptic plasticity in conjunction with the timing of pre- and postsynaptic action potentials. Neural. Comput. 12:385. doi: 10.1162/089976600300015844

PubMed Abstract | CrossRef Full Text | Google Scholar

Letzkus, J., Kampa, B., and Stuart, G. (2006). Learning rules for spike timing-dependent plasticity depend on dendritic synapse location. J. Neurosci. 26, 10420–10429. doi: 10.1523/jneurosci.2650-06.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, W. B., and Stewart, D. (1983). Temporal contiguity requirements for long-term associative potentiation/depression in hippocampus. Neuroscience 8, 791–797. doi: 10.1016/0306-4522(83)90010-6

CrossRef Full Text | Google Scholar

Lisman, J., and Spruston, N. (2005). Postsynaptic depolarization requirements for LTP and LTD: a critique of spike timing dependent plasticity. Nat. Neurosci. 8:839. doi: 10.1038/nn0705-839

PubMed Abstract | CrossRef Full Text | Google Scholar

Lisman, J., and Zhabotinsky, A. M. (2001). A model of synaptic memory: a CaMKII/PP1 switch that potentiates transmission by organizing an AMPA receptor anchoring assembly. Neuron 31, 191–2001.

Google Scholar

Markram, H., Lübke, J., Frotscher, M., Roth, A., and Sakmann, B. (1997a). Physiology and anatomy of synaptic connections between thick tufted pyramidal neurons in the developing rat neocortex. J. Physiol. 500, 409–440. doi: 10.1113/jphysiol.1997.sp022031

PubMed Abstract | CrossRef Full Text | Google Scholar

Markram, H., Lübke, J., Frotscher, M., and Sakmann, B. (1997b). Regulation of synaptic efficacy by coincidence of postsynaptic AP and EPSP. Science 275, 213–215. doi: 10.1126/science.275.5297.213

PubMed Abstract | CrossRef Full Text | Google Scholar

Markram, H., Sjostrom, J., and Gerstner, W. (2011). A history of spike-timing-dependent plasticity. Front. Syn. Neurosci. 3:4. doi: 10.3389/fnsyn.2011.00004

PubMed Abstract | CrossRef Full Text | Google Scholar

Meissner-Bernard, C., Tsai, M., Logiaco, L., and Gerstner, W. (2020). Paradoxical results of long-term potentiation explained by voltage-based plasticity rule. arXiv. doi: arXiv 2001.03614.

CrossRef Full Text | Google Scholar

Mishra, R. K., Kim, S., Guzman, S. J., and Jonas, P. (2016). Symmetric spike timing-dependent plasticity at CA3–CA3 synapses optimizes storage and recall in autoassociative networks. Nat. Commun. 7:11552.

Google Scholar

Nabavi, S., Kessels, H. W., Alfonso, S., Aow, J., Fox, R., and Malinow, R. (2013). Metabotropic NMDA receptor function is required for NMDA receptor-dependent long-term depression. Proc. Natl. Acad. Sci. U S A. 110, 4027–4032. doi: 10.1073/pnas.1219454110

PubMed Abstract | CrossRef Full Text | Google Scholar

Nevian, T., Larkum, M. E., Polsky, A., and Schiller, J. (2007). Properties of basal dendrites of layer 5 pyramidal neurons: a direct patch-clamp recording study. Nat. Neurosci. 10, 206–214. doi: 10.1038/nn1826

PubMed Abstract | CrossRef Full Text | Google Scholar

Ngezahayo, A., Schachner, M., and Artola, A. (2000). Synaptic activity modulates the induction of bidirectional synaptic changes in adult mouse hippocampus. J. Neurosci. 20, 2451–2458. doi: 10.1523/jneurosci.20-07-02451.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Connor, D. H., Wittenberg, G. M., and Wang, S. S. H. (2005). Dissection of bidirectional synaptic plasticity into saturable unidirectional processes. J. Neurophysiol. 94, 1565–1573. doi: 10.1152/jn.00047.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfister, J.-P., and Gerstner, W. (2006). Triplets of spikes in a model of spike timing-dependent plasticity. J. Neurosci. 26. 38, 9673–9682. doi: 10.1523/jneurosci.1425-06.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Remy, S., and Spruston, N. (2007). Dendritic spikes induce single-burst long-term potentiation. Proc. Natl. Acad. Sci. 104, 17192–17197. doi: 10.1073/pnas.0707919104

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubin, J. E., Gerkin, R. C., Bi, G., and Chow, C. C. (2005). Calcium time course as a signal for spike-timing dependent plasticity. J. Neurophysiol. 93, 2600–2613. doi: 10.1152/jn.00803.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Saudargiene, A., and Graham, B. P. (2015). Inhibitory control of site-specific synaptic plasticity in a model CA1 pyramidal neuron. BioSystems 130, 37–50. doi: 10.1016/j.biosystems.2015.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Senn, W., Tsodyks, M., and Markram, H. (2001). An algorithm for modifying neurotransmitter release probability based on pre- and postsynaptic spike timing. Neural. Comput. 13, 35–67. doi: 10.1162/089976601300014628

PubMed Abstract | CrossRef Full Text | Google Scholar

Shouval, H. Z., Bear, M. F., and Cooper, L. N. (2002). A unified model of NMDA receptor-dependent bidirectional synaptic plasticity. Proc. Natl. Acad. Sci. 99, 10831–10836. doi: 10.1073/pnas.152343099

PubMed Abstract | CrossRef Full Text | Google Scholar

Sjöström, P., Turrigiano, G., and Nelson, S. (2001). Rate, timing, and cooperativity jointly determine cortical synaptic plasticity. Neuron 32, 1149–1164. doi: 10.1016/s0896-6273(01)00542-6

CrossRef Full Text | Google Scholar

Song, S., Miller, K. D., and Abbott, L. F. (2000). Competitive Hebbian learning through spike-time-dependent synaptic plasticity. Nat. Neurosci. 3, 919–926. doi: 10.1038/78829

PubMed Abstract | CrossRef Full Text | Google Scholar

Spruston, N. (2008). Pyramidal neurons: dendritic structure and synaptic integration. Nat. Rev. Neurosci. 9, 206–221. doi: 10.1038/nrn2286

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Rossum, M. C. W. V., Bi, G.-Q., and Turrigiano, G. G. (2000). Stable Hebbian learning from spike timing-dependent plasticity. J. Neurosci. 20, 8812–8821. doi: 10.1523/jneurosci.20-23-08812.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H.-X., Gerkin, R. C., Nauen, D. W., and Wang, G.-Q. (2005). Coactivation and timing-dependent integration of synaptic potentiation and depression. Nat. Neurosci. 8, 187–193. doi: 10.1038/nn1387

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, Z., and Storm, D. R. (2005). The role of calmodulin as a signal integrator for synaptic plasticity. Nat. Rev. 6:267276.

Google Scholar

Zenke, F., Agnes, E., and Gerstner, W. (2015). Diverse synaptic plasticity mechanisms orchestrated to form and retrieve memories in spiking neural networks. Nat. Comm. 6:6922.

Google Scholar

Zenke, F., and Gerstner, W. (2017). Hebbian plasticity requires compensatory processes on multiple timescales. Phil. Trans. R. Soc. B 372:20160259. doi: 10.1098/rstb.2016.0259

PubMed Abstract | CrossRef Full Text | Google Scholar