|
|
||||||||
1Department of Neurobiology and Anatomy, University of Texas Medical School at Houston, Houston; and 2Department of Biomedical Engineering, University of Texas, Austin, Texas
Submitted 12 May 2004; accepted in final form 21 September 2004
|
|
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
One form of NMDA receptordependent synaptic plasticity, spike timedependent plasticity (STDP), is induced by repetitively firing pre- and postsynaptic neurons with a fixed time lag of
t between the pre- and postsynaptic spikes. Positive values of
t (pre-post) result in long-term potentiation (LTP), whereas negative values of
t (post-pre) result in long-term depression (LTD) (Bi and Poo 1998
; Feldman 2000
; Froemke and Dan 2002
; Markram et al. 1997
; Sjostrom et al. 2001
; Zhang et al. 1998
).
We previously proposed a model of calcium-dependent synaptic plasticity that can account for various induction protocols of synaptic plasticity (Shouval et al. 2002b
). This theory is based on three fundamental assumptions. 1) Calcium influx into spines determines the sign, magnitude, and rate of synaptic plasticity; a moderate elevation in calcium results in LTD, whereas a large elevation in calcium levels results in LTP (Fig. 1A). We also assume that the rate of plasticity is a monotonically increasing function of calcium (Eqs. 1 and 3). 2) Calcium influx through NMDA receptors is sufficient to account for the various forms of synaptic plasticity. 3) The back-propagating action potentials in the postsynaptic neuron have a slow afterdepolarization. Related models were proposed by other groups (Karmarkar and Buonomano 2002
; Kitajima and Hara 2000
). These calcium-based theories make the testable prediction that, at positive values of
t, larger than those that produce LTP, LTD will be induced (Fig. 1B). This prediction is consistent with results by Nishiyama et al. (2000)
in hippocampal slices. Other experiments have not explored this region of
t extensively; however, from the available data, there is no strong support for a universal form of pre-post LTD (Bi and Poo 1998
; Feldman 2000
; Froemke and Dan 2002
; Markram et al. 1997
; Sjostrom et al. 2001
).
|
t, assuming a simplified stochastic model of synaptic transmission (Yeung et al. 2004
t for
t > 0 (Fig. 1D).
Previous theoretical predictions about the existence of a pre-post form of LTD were obtained using deterministic models of synaptic transmission (Karmarkar and Buonomano 2002
; Kitajima and Hara 2000
; Shouval et al. 2002b
). In this paper, we take into account the stochastic nature of synaptic transmission and show that the shape of the STDP curve can be altered by properties of stochastic synaptic transmission. The finding that the shape of plasticity curves can critically depend on the stochastic nature of synaptic transmission is relevant to other nonlinear biophysical theories of synaptic plasticity as well.
|
|
METHODS |
|---|
|
The calcium-dependent synaptic plasticity equation used here is from Shouval et al. 2002b
![]() | (1) |
is the decay constant, set here to 1, and the
function, as depicted in Fig. 1A, has the form
![]() | (2) |
) = exp(
x)/[1 + exp(
x)], and
1 = 0.4,
2 = 0.65, and
1 =
2 = 30. The calcium-dependent learning rate function,
, has the form
![]() | (3) |
function is based qualitatively on the notion that a moderate rise in calcium produces LTD, whereas a large rise in calcium produces LTP. The actual parameters were set arbitrarily to obtain reasonable plasticity curves, given the calcium transients assumed here. The sigmoidal calcium-dependent learning rate is essential for potentiation during STDP because and equal rate of depression and potentiation would result in oscillations and not robust potentiation. In addition, this monotonically increasing leaning rate is consistent with observations that induction of LTP is faster that induction of LTD (Dudek and Bear 1992
The NMDA receptor kinetics, which we have assumed, are characterized by a double exponential curve of the form
![]() |
f = 50 ms, and
s = 150 ms such that the mean time constant is
av = If x
f + (1 If) x
s = 75 ms, the same time constant as assumed below in the calculations of the CV of calcium transient in the stochastic case. The parameter GNMDA is not fixed as in the stochastic case, but is chosen, as explained below, from a distribution to mimic the stochastic calcium dynamics. When we simulate the deterministic case (Fig. 1B), we used a fixed value GNMDA= 1/325 µM/(ms · V). The voltage dependence of NMDA receptor calcium conductance is as in Shouval et al. 2002b. Calcium concentration is calculated by a first-order differential equation of the form
![]() | (4) |
Ca = 25 ms (Sabatini et al. 2002
The back-propagating action potential has a double exponential form (Shouval et al. 2002b
): BPAP(t) = Vfbpexp(t/
fbp) + Vsbpexp(t/
sbp), here we chose Vfbp = 60 mV, Vsbp = 25 mV,
fbp = 2 ms, and
sbp = 60 ms.
Stochastic calcium dynamics
In previous papers, calcium dynamics were simulated as described in the section above, assuming deterministic calcium dynamics. However, it is better to model synaptic transmission as a stochastic model since given a presynaptic action potential the presynaptic terminal may fail to release neurotransmitter, and if it does release neurotransmitter, both the number of receptors that open and their duration in the open state can vary. In Fig. 1C, we show schematically a histogram of peak calcium transients. The peak at zero is due to release failures from the presynaptic terminal, and the distribution at positive calcium values is due to variability of calcium transients given that there was a release. We have previously calculated analytically both the means and variability of calcium transients invoked by different values of
t, assuming a simplified model of NMDA receptors (Yeung et al. 2004
). These calculations take into account presynaptic variability due to release failures and postsynaptic variability modeled by a simple Markov model of the NMDA receptor. As the number of NMDA receptors increases, the relative postsynaptic contribution to this variability decreases. In Fig. 1D, we show the calculated CV of the peak calcium transients, given that there was a presynaptic release, for 10 postsynaptic NMDA receptors. Anatomical experimental evidence of single postsynaptic spines indicates that the number of postsynaptic NMDA receptors is of the order 1020 (Racca et al. 2000
). Here we assume that the time constant of the NMDA receptor is 75 ms. It is important to note that the Markov model we have used for this analysis is significantly simplified with respect to more realistic models that can account for most of the complexity of NMDA receptor dynamics. However, simulations of calcium dynamics using more complex models of NMDA receptors yield qualitatively similar results (see web APPENDIX A1 ). Intuitively, this phenomenon arises because the number of NMDA receptors that are still open when the back propagating action potential (BPAP) arrives at the synapse and falls as the time between the presynaptic release and the postsynaptic spike (
t) increases, resulting in an effectively smaller number of NMDA receptors as
t increases. This fall in the effective number of NMDA receptors as
t increases is therefore likely to hold for reasonable models of NMDA receptors. At values of
t larger than those depicted in Fig. 1D, the value of the CV drops abruptly to levels equal to those found at large negative values of
t. Therefore at large values of positive
t, larger than those simulated here, the linear fit (Fig. 1D, line) fails.
The shape of calcium transients in this paper is determined by the deterministic calcium dynamics outlined above. However, the magnitude of the transients is modulated by a random process by which we choose the maximal chord conductance parameter GNMDA. The variance of this random process is determined for each
t by the results of the theoretical analysis. The stochastic process for choosing GNMDA has two phases. First, we emulate the effect of failures of presynaptic release by randomly setting GNMDA to zero, with a probability Pfailure= 1 µ = 0.5, where µ is the probability of release. It must be noted that, in our analysis, we do not differentiate between a release failure or a failure to bind postsynaptic receptors when there was a release. Next, if there is a release, we chose the maximal chord conductance, GNMDA randomly from a gamma distribution with a mean of 1/325 µM/(ms · mV), and a variance chosen differently for every
t and the number of NMDA receptors, Z, according to linear fit to the theoretical derivation. The linear fits used in this paper are CV(
t > 0) = 0.095 + 0.0045
t and CV(
t
0) = 0.095 0.00067
t. We also limit the maximal value of GNMDA to be smaller than the mean NMDA receptor conductance times the number of NMDA receptors Z; this has a small effect on the results for small Z. When this procedure is carried out, the CV of the peak calcium transients matches the linear fit to the calculated CV (Fig. 1D, line). Note that we do not change the conductance as a function of Z, only the relative variance. This implies that we are scaling the single NMDA receptor conductance parameter inversely to Z to approximately preserve the form of the plasticity curves, without changing other parameters.
|
|
RESULTS |
|---|
|
t, the CV significantly and monotonically increases as
t increases. This is quite surprising since two different values of
t, one positive and one negative, that produce the same mean peak transient result in significantly different variability. An intuitive explanation for this result is that as
t increases (for
t > 0), the number of NMDA receptors that are still bound to the neurotransmitter when the BPAP arrives decreases; thus the effective number of NMDA receptors is smaller, and the relative level of fluctuations is larger. This intuitive explanation suggests that the qualitative form of the dependence of CV on
t does not depend on the exact details of the Markov model used for the NMDA receptors, and indeed, simulations using more complex Markov models of NMDA receptors produce qualitatively similar results (See web APPENDIX A1). The CV of the calcium transients is inversely proportional to the square root of the number of NMDA receptors at each spine (CV
1/
), and as the number of NMDA receptors increases, the relative postsynaptic contribution to the variance of calcium transients decreases. We have used these analytical results about the dependence of CV on
t and the number of postsynaptic NMDA receptors at each spine (Z) as a basis for a stochastic implementation of the unified plasticity model.
Simulations of STDP with the unified model (Shouval et al. 2002b
), using deterministic calcium dynamics, result in STDP curves that exhibit a pre-post form of LTD (Fig. 1B). In Fig. 2, we show the effect of stochastic synaptic transmission on the shape of the STDP curve. For a small number of NMDA receptors (Z = 10, Fig. 2A) the pre-post form of LTD is nearly eliminated, whereas the post-pre form of LTD is not. This happens because the CV at values of
t, which produce pre-post LTD in the deterministic model (dashed black line), are significantly larger than for values of
t that produce post-pre LTD. For example: CV(60)/CV(10), given that there is a release, is 3.6. Since the omega function (Fig. 1A) is nonmonotonic in the LTD range, a large variability of calcium transients implies that some of the transients can actually produce LTP. Due to the calcium dependent learning rate
, these instances of LTP have a more significant effect than the instances of LTD. Experimental evidence (Racca et al. 2000
) indicates that the number of NMDA receptors in postsynaptic spines is 1020. To compare the STDP curve for Z = 10 to experimental results, we have fit these data points to exponential curves; we obtain time constants
= 14 ms for
t > 0 and
= 57 ms for
t < 0. These results are comparable with experimental results (Feldman 2000
; Froemke and Dan 2002
), although we have not made an attempt to use biophysical parameters appropriate for each system.
|
|
|
DISCUSSION |
|---|
|
t (Karmarkar and Buonomano 2002
t increases, the magnitude of calcium influx through NMDA receptors monotonically and continuously decreases, thus it seems that it must pass through the LTD region before plasticity is completely eliminated at large values of
t. However, these calcium-based models assume deterministic calcium dynamics that represent average values of calcium transients. We found that, due to stochastic properties of NMDA receptors, the pre-post form of LTD predicted by these calcium-dependent models is not a robust prediction of such theories. Instead we propose that, in systems that have a large number of NMDA receptors, and consequently a smaller relative variability of calcium transients, the pre-post form of LTD should be observed. However, if the number of postsynaptic NMDA receptors is small and the variability large, the pre-post form of LTD should have a significantly reduced magnitude.
Our prediction that the CV of peak calcium transients increases as a function of
t (for
t > 0) could be tested in experiments in which calcium transients in spines are imaged (Koester and Sakmann 1998
; Nevian and Sakmann 2004
; Nimchinsky et al. 2004
; Sabatini et al. 2002
). Since the interesting quantity to be experimentally examined is the variability of calcium transients, results would be meaningful if the variability of the method of measurement is small, the number of trials for each synapse large, and estimated variability should be over a single synapse rather than over several synapses, because that will obscure the single spine variability by adding between spine variability.
In the extreme cases (Z < 10) where the pre-post form of LTD is nearly eliminated, most pairs of spikes still fall within the LTD region of the
function. However, the effect of those spikes that fall in the LTP region is larger because we have assumed that the learning rate increases with intracellular calcium levels. This assumption arises from lower level biophysical models (Castellani et al. 2001
; Shouval et al. 2002a
) and is necessary for avoiding oscillations during the induction of spike timedependent LTP. It is also consistent with several experimental observations, first, rate-based induction of LTD is typically induced by 900 (Dudek and Bear 1992
) spikes over 15 min, whereas LTP is much faster and requires less spikes; furthermore, when synaptic plasticity is obtained by direct control of intracellular calcium levels (Yang et al. 1999
), LTP is obtained by a brief large increase in calcium levels, whereas LTD is induced by a moderate and sustained increase in calcium levels.
In the original deterministic model, we assumed calcium influx only through NMDA receptors. This was done because influx through NMDA receptors was sufficient to qualitatively account for the different induction mechanisms (Shouval et al. 2002b
), not because this is the only calcium source. However, imaging of calcium transients when the postsynaptic cell is depolarized (Sabatini et al. 2002
) shows that NMDA receptors are the major source for calcium in spines. Furthermore, when the voltage dependence of NMDA receptors is eliminated using a magnesium-free solution, the difference in calcium transients between pre-post and post-pre conditions was eliminated (Nevian and Sakmann 2004
). Nevertheless, inclusion of other calcium sources such as voltage-gated channels or calcium-induced calcium release from internal stores will have a quantitative effect on results of the model.
The parameters assumed in this paper are generic parameters and are not based on measurements in a specific experimental system. Therefore the consequences of the model should not be taken as quantitative predictions. For example, the LTP and LTD time windows in Fig. 2A depend on specific assumptions about the NMDA receptor and the BPAP time constants, and therefore should not be taken as quantitative predictions of what these time constants should be. This theory, if shown to be qualitatively correct, could be taken to the next quantitative level by using measured parameters of a specific system and making quantitative predictions for that specific system.
We have based our estimate of the number of postsynaptic NMDA receptors of anatomical measurements (Racca et al. 2000
). Recent physiological measurements (Nimchinsky et al. 2004
) indicate that the median number of NMDA receptors open during synaptic transmission is
3, given that there is release. In our analysis, we have not distinguished between probability of release and probability of receptors binding; however, for 10 postsynaptic NMDA receptors (Z = 10) and probability of release of 0.5, our model assumes a mean number of open receptors of 5, above the physiological estimate of 3. The physiological estimate of three bound receptors would imply even larger fluctuations than we have assumed, which would result in an even larger difference between the deterministic and stochastic plasticity curves (see Web APPENDIX B1).
Although the consequences of the unified plasticity model are qualitatively consistent with experimental results, recent experiments indicate that the basic assumptions of this theory might not be. It has been shown, in different systems, that spike timedependent LTD requires additional receptors such as cannabinoid receptors in neocortex (Sjostrom et al. 2003
) and metabotropic glutamate receptors in hippocampus (Nishiyama et al. 2000
). These findings might imply that the basic assumptions of the unified plasticity model need to be revised.
The central finding of this paper is that the stochastic properties of synaptic transmission not only add variability to resulting plasticity curves, but can also alter their shape. The significance of this result is independent of whether the exact assumptions of the unified plasticity model are correct, and it potentially applies to all nonlinear biophysical models of synaptic plasticity.
|
|
GRANTS |
|---|
|
|
|
ACKNOWLEDGMENTS |
|---|
|
|
|
FOOTNOTES |
|---|
1 The Supplementary Material for this article (Appendix A and Appendix B is available online at http://jn.physiology.org/cgi/content/full/00504.2004/DC1. ![]()
Address for reprint requests and other correspondence: H. Shouval, Rm. 7.046, Dept. of Neurobiology. Univ. of Texas Medical School at Houston, 6431 Fannin St., Houston, TX 77030 (E-mail: harel.shouval{at}uth.tmc.edu)
|
|
REFERENCES |
|---|
|
Castellani GC, Quinlan EM, Cooper LN, and Shouval HZ. A biophysical model of bidirectional synaptic plasticity: dependence on AMPA and NMDA receptors. Proc Natl Acad Sci USA 98: 1277212777, 2001.
Cho K, Aggleton JP, Brown MW, and Bashir ZI. An experimental test of the role of postsynaptic calcium levels in determining synaptic strength using perirhinal cortex of rat. J Physiol 532: 459466, 2001.
Cormier RJ, Greenwood AC, and Connor JA. Bidirectional synaptic plasticity correlated with the magnitude of dendritic calcium transients above a threshold. J Neurophysiol 85: 399406, 2001.
Cummings JA, Mulkey RM, Nicoll RA, and Malenka RC. Ca2+ signaling requirements for long-term depression in the hippocampus. Neuron 16: 825833, 1996.[CrossRef][Web of Science][Medline]
Dudek SM and Bear MF. Homosynaptic long-term depression in area CA1 of hippocampus and effects of N-methyl-D-aspartate receptor blockade. Proc Natl Acad Sci USA 89: 43634367, 1992.
Feldman DE. Timing-based LTP and LTD at vertical inputs to layer II/III pyramidal cells in rat barrel cortex. Neuron 27: 4556, 2000.[CrossRef][Web of Science][Medline]
Froemke RC and Dan Y. Spike-timing-dependent synaptic modification induced by natural spike trains. Nature 416: 433438, 2002.[CrossRef][Medline]
Jahr CE and Stevens CF. Voltage dependence of NMDA-activated macroscopic conductances predicted by single-channel kinetics. J Neurosci 10: 31783182, 1990.[Abstract]
Johnston D and Wu MS. Foundation of Cellular Neurophysiology. Cambridge, MA: MIT Press, 1995.
Karmarkar UR and Buonomano DV. A model of spike-timing dependent plasticity: one or two coincidence detectors? J Neurophysiol 88: 507513, 2002.
Kitajima T and Hara K. A generalized Hebbian rule for activity-dependent synaptic modifications. Neural Netw 13: 445454, 2000.[CrossRef][Web of Science][Medline]
Koester HJ and Sakmann B. Calcium dynamics in single spines during coincident pre- and postsynaptic activity depend on relative timing of back-propagating action potentials and subthreshold excitatory postsynaptic potentials. Proc Natl Acad Sci USA 95: 95969601, 1998.
Lisman J. A mechanism for the Hebb and the anti-Hebb processes underlying learning and memory. Proc Natl Acad Sci USA 86: 95749578, 1989.
Markram H, Lubke J, Frotscher M, and Sakmann B. Regulation of synaptic efficacy by coincidence of postsynaptic APs and EPSPs. Science 275: 213215, 1997.
Nevian T and Sakmann B. Single spine Ca2+ signals evoked by coincident EPSPs and backpropagating action potentials in spiny stellate cells of layer 4 in the juvenile rat somatosensory barrel cortex. J Neurosci 24: 16891699, 2004.
Nimchinsky EA, Yasuda R, Oertner TG, and Svoboda K. The number of glutamate receptors opened by synaptic stimulation in single hippocampal spines. J Neurosci 24: 20542064, 2004.
Nishiyama M, Hong K, Mikoshiba K, Poo MM, and Kato K. Calcium stores regulate the polarity and input specificity of synaptic modification. Nature 408: 584588, 2000.[CrossRef][Medline]
Racca C, Stephenson FA, Streit P, Roberts JD, and Somogyi P. NMDA receptor content of synapses in stratum radiatum of the hippocampal CA1 area. J Neurosci 20: 25122522, 2000.
Sabatini BL, Oertner TG, and Svoboda K. The life cycle of Ca2+ ions in dendritic spines. Neuron 33: 439452, 2002.
Shouval H, Castellani GC, Yeung LC, Blais BS, and Cooper LN. Converging evidence for a simplified biophysical model of synaptic plasticity. Biol Cybern 87: 383391, 2002a.[CrossRef][Web of Science][Medline]
Shouval HZ, Bear MF, and Cooper LN. A unified model of NMDA receptor-dependent bidirectional synaptic plasticity. Proc Natl Acad Sci USA 99: 1083110836, 2002b.
Sjostrom PJ, Turrigiano GG, and Nelson SB. Rate, timing, and cooperativity jointly determine cortical synaptic plasticity. Neuron 32: 11491164, 2001.[CrossRef][Web of Science][Medline]
Sjostrom PJ, Turrigiano GG, and Nelson SB. Neocortical LTD via coincident activation of presynaptic NMDA and cannabinoid receptors. Neuron 39: 641654, 2003.[CrossRef][Web of Science][Medline]
Yang SN, Tang YG, and Zucker RS. Selective induction of LTP and LTD by postsynaptic [Ca2+]i elevation. J Neurophysiol 81: 781787, 1999.
Yeung LC, Castellani GC, and Shouval HZ. Analysis of the interspinal calcium dynamics and its implications for the plasticity of spiking neurons. Physiol Rev E 69: 011907011912, 2004.[CrossRef]
Zhang LI, Tao HW, Holt CE, Harris WA, and Poo M. A critical window for cooperation and competition among developing retinotectal synapses. Nature 395: 3744, 1998.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
R. Corlew, D. J. Brasier, D. E. Feldman, and B. D. Philpot Presynaptic NMDA Receptors: Newly Appreciated Roles in Cortical Synaptic Function and Plasticity Neuroscientist, December 1, 2008; 14(6): 609 - 625. [Abstract] [PDF] |
||||
![]() |
X. Yu, H. Z. Shouval, and J. J. Knierim A Biophysical Model of Synaptic Plasticity and Metaplasticity Can Account for the Dynamics of the Backward Shift of Hippocampal Place Fields J Neurophysiol, August 1, 2008; 100(2): 983 - 992. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Urakubo, M. Honda, R. C. Froemke, and S. Kuroda Requirement of an Allosteric Kinetics of NMDA Receptors for Spike Timing-Dependent Plasticity J. Neurosci., March 26, 2008; 28(13): 3310 - 3323. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. C. Gerkin, P.-M. Lau, D. W. Nauen, Y. T. Wang, and G.-Q. Bi Modular Competition Driven by NMDA Receptor Subtypes in Spike-Timing-Dependent Plasticity J Neurophysiol, April 1, 2007; 97(4): 2851 - 2862. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Cai, J. P. Gavornik, L. N. Cooper, L. C. Yeung, and H. Z. Shouval Effect of Stochastic Synaptic and Dendritic Dynamics on Synaptic Plasticity in Visual Cortex and Hippocampus J Neurophysiol, January 1, 2007; 97(1): 375 - 386. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Nevian and B. Sakmann Spine Ca2+ Signaling in Spike-Timing-Dependent Plasticity J. Neurosci., October 25, 2006; 26(43): 11001 - 11013. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. A. Bender, K. J. Bender, D. J. Brasier, and D. E. Feldman Two coincidence detectors for spike timing-dependent plasticity in somatosensory cortex. J. Neurosci., April 19, 2006; 26(16): 4166 - 4177. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. C. Froemke, I. A. Tsay, M. Raad, J. D. Long, and Y. Dan Contribution of Individual Spikes in Burst-Induced Long-Term Synaptic Modification J Neurophysiol, March 1, 2006; 95(3): 1620 - 1629. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. Rubin, R. C. Gerkin, G.-Q. Bi, and C. C. Chow Calcium Time Course as a Signal for Spike-Timing-Dependent Plasticity J Neurophysiol, May 1, 2005; 93(5): 2600 - 2613. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |