A Ca2+-Based Computational Model for NMDA Receptor-Dependent Synaptic Plasticity at Individual Post-Synaptic Spines in the Hippocampus

Associative synaptic plasticity is synapse specific and requires coincident activity in pre-synaptic and post-synaptic neurons to activate NMDA receptors (NMDARs). The resultant Ca2+ influx is the critical trigger for the induction of synaptic plasticity. Given its centrality for the induction of synaptic plasticity, a model for NMDAR activation incorporating the timing of pre-synaptic glutamate release and post-synaptic depolarization by back-propagating action potentials could potentially predict the pre- and post-synaptic spike patterns required to induce synaptic plasticity. We have developed such a model by incorporating currently available data on the timecourse and amplitude of the post-synaptic membrane potential within individual spines. We couple this with data on the kinetics of synaptic NMDARs and then use the model to predict the continuous spine [Ca2+] in response to regular or irregular pre- and post-synaptic spike patterns. We then incorporate experimental data from synaptic plasticity induction protocols by regular activity patterns to couple the predicted local peak [Ca2+] to changes in synaptic strength. We find that our model accurately describes [Ca2+] in dendritic spines resulting from NMDAR activation during pre-synaptic and post-synaptic activity when compared to previous experimental observations. The model also replicates the experimentally determined plasticity outcome of regular and irregular spike patterns when applied to a single synapse. This model could therefore be used to predict the induction of synaptic plasticity under a variety of experimental conditions and spike patterns.

to vary Ca 2+ influx through NMDARs by fixing the post-synaptic membrane potential during low frequency synaptic stimulation (Isaac et al., 1995;Daw et al., 2000).
Spike timing-dependent plasticity (STDP) is a form of Hebbian synaptic plasticity that incorporates a temporal specificity to coincident pre-and post-synaptic activity. In the hippocampus, STDP was originally thought to be induced by single pairs of pre-and post-synaptic action potentials such that if the pre-synaptic action potential occurs before the post-synaptic action potential LTP is induced whereas if the order of action potentials is reversed then LTD is induced (Bi and Poo, 1998;Debanne et al., 1998;Nishiyama et al., 2000;Campanac and Debanne, 2008;Kwag and Paulsen, 2009). Other data have proposed this model should include a requirement for bursts of post-synaptic action potentials for the induction of LTP although LTD may be induced by single pairs (Pike et al., 1999;Wittenberg and Wang, 2006;Buchanan and Mellor, 2007) reviewed in . This is a divergence from the situation at cortical synapses where single pairs of action potentials can induce both LTP and LTD (Sjostrom et al., 2001;Sjostrom and Nelson, 2002;Froemke et al., 2006, but see Nevian and Sakmann, 2006).

IntroductIon
Hebbian synaptic plasticity is the cellular and molecular correlate of associative learning in the brain. During presentation of information that needs to be retained for future use, specific synapses are subjected to activity patterns that induce a long-term change in synaptic strength. For Hebbian synaptic plasticity at Schaffer collateral synapses in the hippocampus, these patterns require coincident activity in pre-and post-synaptic neurons to activate NMDA receptors (NMDARs) present on the membrane of the post-synaptic dendritic spine. The resulting Ca 2+ influx through NMDARs is the critical trigger for induction of synapse specific plasticity (Lisman, 1989).
Classically, high frequency synaptic stimulation induces longterm potentiation (LTP) whereas low frequency stimulation induces long-term depression (LTD) suggesting that brief high concentrations of Ca 2+ in the post-synaptic spine induce LTP whereas prolonged lower concentrations of Ca 2+ induce LTD (Bear et al., 1987;Hansel et al., 1996). This hypothesis is supported by measurements of Ca 2+ concentration during plasticity induction (Hansel et al., 1997;Cho et al., 2001;Cormier et al., 2001;Ismailov et al., 2004;Gall et al., 2005) and by plasticity induction protocols designed Since Ca 2+ influx through NMDARs is pivotal for LTP and LTD, this suggests the induction of synaptic plasticity can be predicted by NMDAR opening kinetics in response to pre-synaptic glutamate release and post-synaptic depolarization. This approach has been adopted for the modeling of post-synaptic calcium dynamics in response to synaptic stimulation or back-propagating action potentials (Franks et al., 2002;Grunditz et al., 2008;Keller et al., 2008) and to STDP induction protocols (Shouval et al., 2002;Rubin et al., 2005;Graupner and Brunel, 2007;Helias et al., 2008;Urakubo et al., 2008;Castellani et al., 2009). However, these STDP models are limited by the experimental data used to determine their parameters and, in addition, ought to accurately predict the plasticity outcomes of a variety of induction protocols. Recent advances in dendritic spine imaging provide data on spine depolarization and Ca 2+ concentrations in response to pre-and post-synaptic action potentials (Sabatini et al., 2002;Nevian and Sakmann, 2006;Bloodgood and Sabatini, 2007;Canepari et al., 2007;Palmer and Stuart, 2009) that potentially greatly increase the accuracy of such models of plasticity induction.
We have developed a computational model of synaptic plasticity induction based on one originally described by Shouval et al. Our model incorporates the latest experimental data on dendritic spine depolarization and Ca 2+ dynamics. We also test the predictive power of our model on many plasticity induction protocols by calculating continuous Ca 2+ concentrations during long induction periods. We find that our model accurately predicts the experimental data tested and we hypothesize that it can thus be used to search for instances of synaptic plasticity induction during continuous activity at Hebbian synapses in the hippocampus.

MaterIals and Methods
We use a physiologically plausible model based on intracellular Ca 2+ dynamics caused by NMDAR activation during the induction of synaptic plasticity to predict the plasticity outcome of any set of pre-and post-synaptic activity patterns that occur at the Schaffer collateral synapse in the hippocampus. Since we are interested in studying experimental spike trains we modify a model originally proposed by Shouval et al. (2002) to allow us to carry out such analysis. We make a number of critical modifications to analyze the Ca 2+ dynamics in individual dendritic spines during long periods of irregular spiking activity. This is illustrated using short epochs of overlapping hippocampal place cell activity (Figure 1) (Isaac et al., 2009).
Essential components of the experimental spike trains for the activation of NMDARs and therefore the induction of synaptic plasticity are (i) the pre-synaptic release of glutamate that dictates the binding of glutamate to the NMDARs and (ii) the postsynaptic membrane potential that determines the relative blockade of NMDARs by Mg 2+ . For the purposes of this model, we have assumed the two events that determine the post-synaptic membrane potential within an individual dendritic spine are excitatory post-synaptic potentials (EPSPs) and back-propagating action potentials (BPAPs). We start by modeling the BPAPs as follows: BPAP ,  where V max bs is the maximum depolarization due to the BPAP, I f bs , and I s bs are the relative magnitudes of the fast and slow components of the BPAP, respectively, that sum to one, and the integration time step δ is 0.1 ms. Due to the slower (and much smaller) after-depolarizing  I f and I s are the relative magnitudes of the fast and slow component of the NMDAR current as a result of glutamate binding, respectively, that sum to one, and Θ is the Heaviside (unit) step function. The voltage dependence of the current that takes into account Mg 2+ block of the receptor (Jahr and Stevens, 1990) is represented by is the driving force determined by the reversal potential, V r1 (0 mV), and We then calculate the spine membrane potential as the summation of BPAP and EPSP AMPA and EPSP NMDA : where V rest is set at −65 mV unless otherwise stated. An example of the predicted spine voltage can be seen for a sample epoch of overlapping place cell activity in Figure 1. Since NMDARs provide the major source of Ca 2+ influx into post-synaptic dendritic spines (Bloodgood and Sabatini, 2007), we incorporate in our model the Ca 2+ current through NMDAR that takes the following form (Shouval et al., 2002): This is similar to Eq. 3 except for the terms P 0 and G NMDA that represent the open channel probability and NMDAR Ca 2+ conductance respectively and V r2 is the reversal potential for calcium (130 mV).
Next, the rate of change of the [Ca 2+ ] inside the post-synaptic spine is governed by: where α is a factor that converts current to flux and τ Ca is the calcium passive decay time constant. An example of the Ca 2+ current flow through NMDARs and the resulting predicted [Ca 2+ ] in the spine can be seen in Figure 1.
Finally, we assume that spine [Ca 2+ ] is the trigger for synaptic strength change. For the purposes of our study the continuous model for synaptic strength used in Shouval et al. is modified to act as a Ca 2+ -gated function based on local peaks in [Ca 2+ ] as follows: potential, if two spikes happen near enough to each other that the first spike is still decaying, the effect of the BPAPs is additive. Since we are modeling the BPAP at the spine V max bs is set at 67 mV in line with experimental data measuring membrane potential in spines with voltage-sensitive dyes (Canepari et al., 2007;Palmer and Stuart, 2009). This is smaller than the maximum BPAP amplitude found at the soma used by Shouval et al. An example of the modeled BPAP during place cell activity can be seen in Figure 1.
The equation that governs the behavior of AMPAR-mediated EPSPs in the model is similar to (1) having a slow and a fast exponential component: where the parameter N a reflects the maximum effect that a single AMPAR-mediated EPSP can have. The value of N a can vary depending on the number of synapses activated. Activation of a single synapse results in a membrane depolarization in the spine of approximately 10 mV (Palmer and Stuart, 2009). Again this deviates from the value of 1 mV recorded at the soma and used by Shouval et al. Assuming that the maximum depolarization that a single EPSP can generate is 10 mV we define N a in the following way: EPSP AMPA also depends on the membrane potential, V m . This dependence is represented by the term (V m −V r1 )/V rest where V r1 is the reversal potential for AMPARs (0 mV) and V rest is the resting membrane potential (−65 mV).
The equation that governs the behavior of NMDA-mediated EPSPs in the model has the following form: where the parameter N n reflects the maximum effect of the NMDAR-mediated component of the EPSP. This is calculated in a similar fashion to N a for EPSP AMPA using a value of 5 mV for the NMDAR-mediated EPSP at −65 mV in the absence of Mg 2+ measured by dendritic recordings (Fernandez de Sevilla et al., 2007).
Having validated our model for the observed Ca 2+ influx at dendritic spines we next asked the question if the model could replicate experimental data for the induction of synaptic plasticity using a variety of protocols. We have restricted our model to comparisons with experimental data from the Schaffer collateral synapse of the hippocampus and not considered other synapses in other brain regions.

spIke tIMIng-dependent plastIcIty wIth paIrs of pre-and post-synaptIc spIkes
To model STDP with pairs of pre-and post-synaptic spikes we initially assumed single synaptic activation and varied ∆t between −20 and + 100 ms at intervals of 0.1 ms measuring the peak [Ca 2+ ] at each value of ∆t ( Figure 3A). [Ca 2+ ] rose from its baseline of 72 nM (the peak [Ca 2+ ] attained for a single EPSP in isolation) to a peak of 230 nM at ∆t ≈10 ms ( Figure 3B). Experiments such as these have been shown to generate no significant synaptic plasticity (Buchanan and Mellor, 2007) whereas those using larger amplitude EPSPs have been shown to generate LTD (Wittenberg and Wang, 2006). We estimated the activation of multiple synapses at the same time would increase the depolarization within a single dendritic spine during an EPSP from 10 to 20 mV based on experimental predictions (Palmer and Stuart, 2009). Thus we have estimated that the activation of other spines will contribute an additional 10 mV of depolarization within an activated spine above and beyond the experimentally determined 10 mV for activation of a single synapse. This doubling of the EPSP amplitude resulted in an increase in peak [Ca 2+ ] at all values of ∆t with a peak of 279 nM occurring at ∆t ≈10 ms ( Figure 3C).
It has also been shown that the frequency of spike pairing is important for the induction of plasticity such that at higher frequencies (>5-10 Hz) LTP can be induced (Wittenberg and Wang, 2006;Buchanan and Mellor, 2007). We varied the frequency of spike pairings in our model for 10 mV EPSPs over a range of frequencies from 1 to 100 Hz ( Figure 3D). Summation of Ca 2+ transients was found to occur at frequencies greater than ∼5 Hz indicating that increasing the frequency will shift the STDP protocol towards larger [Ca 2+ ] and therefore LTP in line with the experimental data. where The critical target for Ca 2+ influx through NMDARs is the enzyme CAMKII. Due to its ability to autophosphorylate, the activation of this molecule can be long lived and the level is determined by local peak [Ca 2+ ]. Thus, synaptic weight change is determined at local peak [Ca 2+ ] (Miller et al., 2005;Graupner and Brunel, 2007;Helias et al., 2008;Urakubo et al., 2008;Castellani et al., 2009). Since there is no noise associated with our model, these peaks are measured instantaneously without smoothing. Experimentally, increases in synaptic weight tend towards saturation as synaptic weight increases. In addition, decreases require synaptic weight to always be >0. These constraints explain the form of Eq. 7.
Numerical integration was performed using forward Euler method implemented in MATLAB.
The parameter values used in the simulations are given for completeness in Table 1.

results
Our starting point for developing a model for the induction of synaptic plasticity was to incorporate the most recent and accurate measurements of voltage changes within dendritic spines using data from measurements of voltage-dependent dyes (Canepari et al., 2007;Palmer and Stuart, 2009). We model the membrane potential at the spine rather than the soma because this is the site of the NMDARs critical for the induction of synaptic plasticity. This shifts the determination of membrane depolarization away from BPAPs and towards EPSPs since the former attenuate as they pass along the dendrite and the latter are now measured at their site of origin. This is a departure from previous models that used values for BPAPs and EPSPs recorded at the soma (Shouval et al., 2002). With this change, our model predicts that an EPSP resulting from the activation of a single synapse is sufficient to cause a significant Ca 2+ influx through NMDARs (Figure 2A) in line with experimentally observed data (Bloodgood and Sabatini, 2007;Canepari et al., 2007;Sobczyk and Svoboda, 2007). The pairing of a BPAP with a single EPSP with a time delay of 10 ms produces 3-4 times the Ca 2+ influx (Figure 2A) that again agrees qualitatively with experimentally observed data (Bloodgood and Sabatini, 2007). For comparison we changed the maximal EPSP and BPAP amplitudes to those known to occur at the soma (∼1 and ∼100 mV respectively). With these parameters, a single EPSP produces limited Ca 2+ influx whereas pairing an EPSP with a BPAP produces a large Ca 2+ influx ( Figure 2B).

Rackham et al.
Model for hippocampal synaptic plasticity plasticity is induced (Buchanan and Mellor, 2007). When we used theta burst stimulation with five stimuli to only the pre-synaptic input, the model predicted peak [Ca 2+ ] within the spine to be 325 nM ( Figure 5) and with four stimuli 250 nM. The value for four stimuli is more physiologically relevant since the probability of neurotransmitter release at any one Schaffer collateral synapse is considerably less than 1. Therefore it is highly unlikely that an experimental theta burst will ever generate five EPSPs at an individual synapse.
paIrIng post-synaptIc depolarIzatIon wIth pre-synaptIc stIMulatIon Other common synaptic plasticity induction protocols have dispensed with the need for post-synaptic spikes altogether and use voltage clamp to depolarize the post-synaptic membrane and allow NMDAR activation. This technique neatly demonstrates the bidirectional nature of NMDAR-dependent plasticity since depolarization to moderate levels (−40 mV) produces LTD whereas higher depolarization (0 mV) produces LTP (Isaac et al., 1995;Daw et al., 2000). We tested this with our model by clamping the membrane potential (V m ) at either −40 or 0 mV (Figure 6). Peak [Ca 2+ ] in response to EPSPs were 336 nM and 2.43 μM respectively, which when compared to peak [Ca 2+ ] produced by other protocols would be expected to induce LTD and LTP respectively in agreement with experimental data.

the ca 2+ hypothesIs can explaIn prevIous experIMental data
The Ca 2+ hypothesis states that brief high concentrations of Ca 2+ in the post-synaptic spine induce LTP whereas prolonged lower concentrations of Ca 2+ induce LTD (Bear et al., 1987;Hansel et al., 1996). This is expressed graphically in Figure 7. Points are indicated representing the predicted [Ca 2+ ] from our model for specific plasticity inducing protocols. STDP with single pairs of BPAPs and small EPSPs do not induce plasticity (Buchanan and Mellor, 2007) but when large EPSPs are used LTD is induced (Wittenberg and Wang, 2006) and STDP with triplets of single EPSPs and bursts of BPAPs produces LTP (Pike et al., 1999;Wittenberg and Wang, 2006;Buchanan and Mellor, 2007). When the post-synaptic membrane potential is set at −40 mV during pre-synaptic stimulation LTD is induced (Daw et al., 2000) whereas at 0 mV LTP is induced (Isaac et al., 1996). Theta burst pairing also induces LTP (Frick et al.,

spIke tIMIng-dependent plastIcIty wIth trIplets of spIkes
Post-synaptic burst firing has been shown to be important for the induction of LTP at Schaffer collateral synapses in the hippocampus where burst firing in this instance refers to any number of spikes greater than one (Pike et al., 1999;Wittenberg and Wang, 2006;Buchanan and Mellor, 2007). We tested this on our model using triplets of spikes composed of one pre-synaptic spike and two postsynaptic spikes where ∆t is the time between the pre-synaptic spike and the first post-synaptic spike and ∆s is the delay between the two post-synaptic spikes ( Figure 4A). We first used 10 mV EPSPs with a constant ∆s of 10 ms and varied ∆t between −20 and +100 ms. This produced a peak [Ca 2+ ] of 420 nM at ∆t = 4 ms which increased to a peak [Ca 2+ ] of 475 nM when 20 mV EPSPs were used ( Figure 4B) confirming that spike triplets produce higher peak [Ca 2+ ] than spike pairs and therefore are more likely to induce LTP. We next varied ∆s whilst maintaining ∆t constant at 10 ms for both 10 and 20 mV EPSPs revealing a decrease in peak [Ca 2+ ] as ∆s increases ( Figure 4C). Finally, we varied the frequency of triplets for 10 mV EPSPs over a range of frequencies from 1 to 100 Hz whilst keeping ∆t and ∆s constant at 10 ms each ( Figure 4D). Summation of Ca 2+ transients was found to occur at frequencies greater than ∼4 Hz.

theta burst plastIcIty
We now moved away from STDP to look at other common synaptic plasticity induction protocols. The theta burst protocol was developed to mimic the activity patterns believed to occur at hippocampal synapses during learning and consists of bursts of four or five spikes at 100 Hz with an interburst interval of 200 ms. These can either be applied to the pre-and post-synaptic neuron coincidentally (Frick et al., 2004) or to just the pre-synaptic neuron (Larson et al., 1986). The latter then leads to post-synaptic spikes through EPSP summation if the initial EPSP amplitude is sufficiently large (Buchanan and Mellor, 2007). We used our model to mimic coincident theta burst activity in pre-and post-synaptic neurons using 10 mV EPSPs and found that this type of synaptic stimulation produces very large peak [Ca 2+ ] within dendritic spines (Figure 5) indicating that this protocol is very efficient at producing LTP in agreement with experimental data. Experimental data also shows when theta burst stimulation is given to only the pre-synaptic neuron without initiating action potentials then no

InductIon of synaptIc plastIcIty by hIppocaMpal place cell fIrIng patterns
One of the main purposes for developing a model that is capable of continuously monitoring [Ca 2+ ] in spines and therefore predicts changes in synaptic strength is to scan long periods of neuronal activity for epochs that would be expected to induce plasticity without having to directly measure synaptic strength. To test if the model could perform this task we used data from experiments where long sections of hippocampal place cell activity were replayed into single hippocampal synapses to test the plasticity outcome (Isaac et al., 2009(Isaac et al., ). 2004Buchanan and Mellor, 2007) whereas theta burst to only pre-synaptic inputs does not (Buchanan and Mellor, 2007). In this instance the absolute [Ca 2+ ] values predicted by the model are not as important as the relative magnitudes between plasticity induction protocols. However, it is interesting to note that the absolute [Ca 2+ ] values predicted by the model broadly agree with those measured experimentally for the induction of synaptic plasticity at Schaffer collateral synapses on CA1 pyramidal neurons (Cormier et al., 2001). Thus the predictions from our model support the Ca 2+ hypothesis for synaptic plasticity induction.  (Figures 9A-D). We also tested two further pairs of place cells with non-overlapping or adjacent place fields (1 A and 1 C , 2 E and 2 D ) and found the model predicted only a small LTD for the non-overlapping pair and a small LTP for the adjacent pair (Figures 9E-F).
Our original experimental data also tested a pair of place cells that had an asymmetric cross-correlation such that cell 1 A preferentially fired just before cell 1 B . Because classical STDP rules state that the temporal order of pre-and post-synaptic spikes controls the direction of synaptic plasticity (Bi and Poo, 1998;Song et al., 2000), the existence of this asymmetry suggested that when cell 1 A was pre-synaptic and cell 1 B post-synaptic then LTP would be induced but if the cells were reversed then LTD would be induced.
We first took an ∼16-min period of activity from a pair of place cells (1 A and 1 B ) that had overlapping place fields and therefore would be expected to fire at approximately the same time (Isaac et al., 2009). During the ∼16-min period short coincident bursts of activity could be seen in the two place cells that the model predicted would produce large [Ca 2+ ] sufficient to induce LTP (Figures 1  and 8). This LTP was initiated in the first few minutes of activity and eventually reached a plateau.
We tested a set of four further pairs of place cells (2 A and 2 B , 2 C and 2 D , 3 A and 3 B , 4 A and 4 B ) with overlapping place fields but with strikingly different spiking characteristics [for a full description of the place cell spike pattern characteristics and plasticity outcomes see Isaac et al. (2009)] and found in each case the model predicted  (Isaac et al., 2009). Finally we have compared the experimentally determined plasticity outcome from nine pairs of place cells with the outcome predicted by our model. We find that the correlation between the predicted and observed values is significant (Figure 10C, r 2 = 0.58, P < 0.05 by linear regression) and therefore conclude that the model successfully predicts the induction of synaptic plasticity by irregular activity patterns.

dIscussIon
The model described in this study incorporates two important components of Ca 2+ dynamics in dendritic spines that are necessary for the induction of synaptic plasticity. Firstly, our model is capable of analyzing Ca 2+ influx and concentration continuously and therefore it can determine the plasticity outcome of multiple synaptic events that occur in vivo in an irregular pattern. Secondly, [Ca 2+ ] is modeled at the synapse in dendritic spines rather than at the soma. This is important since the critical Ca 2+ signal for the induction of synaptic plasticity occurs at the spine. It also changes the relative importance of EPSP vs BPAP depolarization which has major implications for the predicted induction of STDP. This approach is validated by comparison of the predicted vs observed Ca 2+ transients in response to either a single EPSP or coupled with a BPAP (Figure 2) (Bloodgood and Sabatini, 2007).
The absolute values for [Ca 2+ ] within the dendrite required for the induction of synaptic plasticity have been estimated as 150-500 nM for LTD and >500 nM for LTP (Cormier et al., 2001). However, other researchers have estimated [Ca 2+ ] within a spine  However, when we reversed the place cell firing patterns such that cell 1 B was pre-synaptic and cell 1 A post-synaptic the model predicted LTP (Figure 10A) that corroborates the experimental results and closely reproduces the experimentally determined timecourse of LTP development (Isaac et al., 2009). We also manipulated the spike patterns in cell 1 B to remove all spikes that occurred less than 100 ms after a spike in cell 1 A leaving only spikes that occurred before any spike in cell 1 A . Classical STDP rules would again predict that if cell 1 A was pre-synaptic and the modified cell 1 B was  in response to a single EPSP at 700 nM and a much higher 12 μM during pairing of post-synaptic depolarization with synaptic stimulation (Sabatini et al., 2002). This discrepancy could be explained in a number of ways.  Figure 7 that predicts the plasticity outcome. Inhibitory synaptic transmission has a major role regulating the induction of synaptic plasticity in the hippocampus. The transient depression of inhibition induced by activation of pre-synaptic cannabinoid or GABA B receptors facilitates the induction of LTP (Davies et al., 1991;Chevaleyre and Castillo, 2004). This modulation of synaptic plasticity is not included in our current model but incorporation of the hyperpolarizing effects of GABAergic transmission would be an important future improvement and might, for example, contribute to the frequency dependence of STDP induction.
NMDARs are not the only sources of Ca 2+ within dendritic spines but are certainly the most important for the induction of synaptic plasticity. A role has also been demonstrated for Ca 2+ stores present in the endoplasmic reticulum in dendrites and spines (linked to Ca 2+ influx through NMDARs or mGluRs) and also voltage-dependent   , 1 B (B). (C) A comparison between the induced plasticity predicted by the model and the observed plasticity from experimental data (Isaac et al., 2009). receptor activation no plasticity is induced but in the presence of a muscarinic agonist we obtain the results described above. This suggests that muscarinic receptor activation is permissive for the induction of synaptic plasticity under conditions where single synapses are activated and EPSP amplitude is small. A similar role has also been suggested for other G-protein coupled receptors such as mGluRs (Bashir et al., 1993) and β-adrenoreceptors (Seol et al., 2007). The observation that muscarinic receptor activation is not required when multiple synapses are activated and EPSP amplitude is large also demonstrates the critical role that multiple synaptic activation plays in the induction of synaptic plasticity. The model described in this study does not currently include a role for G-protein-mediated modulation of synaptic plasticity and this is an important consideration for future modifications and improvements. For example, it would be potentially possible to incorporate a negative feedback on Ca 2+ influx under normal conditions that is inhibited by activation of G-protein coupled receptors or overridden by larger EPSPs.
Technical considerations make the visualization or recording of single synapses during behavior extremely difficult. This makes determination of the events that induce synaptic plasticity in the hippocampus almost impossible. These limitations could be circumvented by using the model described in this study to predict the plasticity outcome at a single synapse for any pattern of pre-and post-synaptic activity. Potentially, this will enable the determination of when synaptic plasticity is induced during behavior by examining the spike patterns of pairs of cells recorded during a memory task. This represents an exciting future application for this model. acknowledgMents Jack R. Mellor is supported by the Wellcome Trust and by the European Union ENI-NET. The Bristol Centre for Complexity Sciences is funded by the EPSRC, UK (EP/E501214/1) which supports OJLR and AG. calcium channels present on the membrane of dendritic spines that generate local dendritic Ca 2+ spikes (Golding et al., 2002;Remy and Spruston, 2007) and can contribute to some forms of LTD (Nevian and Sakmann, 2006). The contribution of non-NMDAR-mediated Ca 2+ increases to synaptic plasticity induction is specifically omitted from our model but would be a useful future addition. There are also potential alternative mechanisms for the induction of synaptic plasticity that do not involve Ca 2+ such as mGluRs (Moult et al., 2006) or muscarinic acetylcholine receptors (Dickinson et al., 2009) that are also specifically omitted from our model. Future experimental data may, however, require their addition.
There is a wealth of data on the spiking requirements for plasticity at various synapses in the cortex (Dan and Poo, 2006). However, these requirements have been shown to be somewhat different to those found in the hippocampus perhaps due to variations in dendritic architecture and NMDAR subunit expression and therefore we have restricted our model specifically to the Schaffer collateral synapse in the hippocampus. Even within the hippocampus there is some controversy surrounding the precise spiking requirements for STDP  for example whether single pairs of spikes at low frequency can induce LTP (Bi and Poo, 1998;Debanne et al., 1998;Nishiyama et al., 2000;Wittenberg and Wang, 2006;Buchanan and Mellor, 2007;Campanac and Debanne, 2008;Kwag and Paulsen, 2009). These discrepancies can potentially be explained by variations in experimental conditions [For example, size of EPSPs, age of animals, frequency of stimulus presentation, presence of inhibitory inputs, etc. discussed in ]. To reduce these discrepancies our model is based on data from acute hippocampal slices but future modifications could enable it to accurately replicate data from all systems including cortical synapses if one assumes that the general mechanism for the induction of synaptic plasticity is similar for all synapses that have NMDAR-dependent synaptic plasticity.
Our experimental data has revealed a critical role for muscarinic acetylcholine receptors in facilitating LTP induced by place cell firing patterns (Isaac et al., 2009). In the absence of muscarinic