|
|
||||||||
1Neuroscience Program, Department of Biological Sciences, Ohio University, Athens, Ohio; and 2Department of Pharmacology, Physiology and Toxicology, Marshall University School of Medicine, Huntington, West Virginia
Submitted 8 March 2006; accepted in final form 15 May 2006
|
|
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
The magnitude of potentiation with LTP is often quantified over multiple cells from field potentials measured with extracellular electrodes, but this potentiation represents an average over many synapses in many cells. LTP has also been studied in individual cells where extracellular stimulation activates many synapses (Gustafsson et al. 1987
; Sastry et al. 1986
). The average potentiation measured with these methods is typically 50100%. In other studies, minimal stimulation techniques have been used in an attempt to assess potentiation at individual synapses. Friedlander et al. (1990)
report that unitary excitatory postsynaptic potentials (EPSPs) were rarely potentiated, whereas compound EPSPs could be potentiated 5984% on average. Debanne et al. (1999)
found that, while multiunitary EPSPs were potentiated 75% on average, LTP of unitary EPSPs was quite heterogeneous; some synapses could not be potentiated (Petersen et al. 1998
), whereas others could be potentiated >10-fold (mean of 250%). Similarly, Isaac et al. (1996)
report individual synapses were potentiated by 145212%.
These average LTP percentages coupled with synapse LTP heterogeneity do not reveal what actually happens quantitatively at individual synapses. If LTP causes the addition of AMPA receptor channels, how many channels per activated synapse would have to be added to get the observed potentiation? If AMPA single channel conductance increases, how many channels per synapse would have to exhibit a change of a certain magnitude to manifest the observed potentiation? How many "silent" synapses would have to be unmasked or new connections formed? Here we combine experimental and modeling studies to address two specific issues. First, we seek to answer the above questions quantitatively to see if the required numbers make sense physiologically and to see if any of these mechanisms is sufficient by itself to explain LTP. Second, we address how closely the potentiation percentage explains the magnitude of actual changes at the synaptic level given inherent nonlinearities in synaptic and neuron level properties.
|
|
METHODS |
|---|
|
SLICE PREPARATION AND MAINTENANCE. Sprague-Dawley rats (3060 days old) were sedated by inhalation of a CO2/air mixture and decapitated. The brain was removed, placed into chilled artificial cerebrospinal fluid (ACSF; equilibrated with 95% O2-5% CO2), trimmed with a razor blade, glued to the stage of a vibratome, and sectioned into 400-µm-thick slices. After sectioning, brain slices were dissected further to free the hippocampus from surrounding structures. To prevent spontaneous bursting in the presence of the GABAA receptor antagonist bicuculline, area CA3 was removed by a small cut.
Hippocampal slices were stored in a holding chamber at the interface of ACSF and humidified 95% O2-5% CO2 at room temperature. Individual slices were transferred, as needed, to a small volume (
200 µl) interface recording chamber (35°C) where they were perfused at 11.5 ml/min with ACSF. The upper surfaces of the slices were exposed to warmed, humidified atmosphere consisting of 95% O2-5% CO2.
ELECTROPHYSIOLOGY.
Intracellular recordings were made using the "blind" whole cell recording technique originally described by Blanton et al. (1989
; also see Grover 1998
; Grover and Chen 1999
). The pipette solution was composed of 140 mM potassium gluconate, 10 mM Na-HEPES, and 3 mM MgCl2, pH 7.2, 280290 mOsm. Potentials were recorded with an Axoclamp 2B (Axon Instruments). Membrane potentials were low-pass filtered at 1 kHz and signals were fed through an Ithaco 4302 dual filter, where they were amplified (10100x), filtered (DC 3 kHz), digitized at 740 kHz, and stored on a personal computer. Access resistance (compensated using the Axoclamp bridge circuit) and membrane resistance were monitored throughout the recordings. Recordings were terminated if large or abrupt changes in access or membrane resistance occurred.
Recordings of voltage at the soma were made in response to four types of stimuli: brief (20 ms) and long (200 ms) 50 pA current injection through the recording electrode and near minimal stimulation (0.1 Hz) through two stimulating electrodes placed in mid-stratum radiatum on opposite sides of the recording site. Bipolar stimulating electrodes were constructed of teflon-insulated stainless steel, and constant voltage stimulation was used (0.1-ms duration, 0.67.0 V). One stimulating electrode was used to deliver low frequency test stimuli and high-frequency tetanic stimulation for LTP induction, and the second electrode was used to deliver test stimuli only to a control pathway. Synaptic stimulation was interleaved with intracellular current stimulation (50 pA, 20 or 200 ms). Stimulation (synaptic and intracellular) was delivered every 5 s, so that synaptic stimulation was given once every 10 s, and intracellular current stimulation was also given once every 10 s. Recordings were made in the presence of 10 µM bicuculline to block fast, GABAA receptormediated inhibition. Test stimulus intensity was initially adjusted so that failures were observed on some trials. Synaptic potentials were recorded at the initial resting membrane potential of the cell; any systematic drift in membrane potential was compensated by DC current injection. Test stimulation was delivered until a minimum 5-min period of stable recording was achieved. Responses for each of the four types of stimuli were averaged over a 5-min period. All membrane potentials were corrected for the calculated liquid junction potential (10 mV).
LTP INDUCTION AND MEASUREMENT. After achieving a stable baseline recording, LTP induction was attempted. Afferent stimulation was increased to 25 Hz for 1 s. During LTP induction, each afferent stimulus was paired with a 1.5-nA, 5-ms current injected through the whole cell pipette. This protocol resulted in LTP in >85% of the cells recorded. The same stimulation protocol was used before and after LTP induction. Responses to short and long current pulses at the soma and to EPSPs evoked in tetanized and control pathways were monitored during the posttetanus period for as long as a stable recording could be maintained. Posttetanus recording durations ranged from 10 to 30 min or more, over which time LTP magnitude was typically stable. Potentials were collected and initially analyzed using the WinWCP software (John Dempster, University of Strathclyde). Further analysis was accomplished using Origin (OriginLab) and custom software written in the Python programming language.
PHARMACOLOGICAL INHIBITION OF IH.
To eliminate possible effects of the hyperpolarization activated cation current Ih (Maccaferri et al. 1993
) on EPSP amplitude and kinetics, some cells were recorded in the presence of 2550 µM ZD-7288 (Gasparini and DiFrancesco 1997
; Harris and Constanti 1995
). Our intention in these recordings was to determine if the parameter fitting procedure described below would be easier without this current. In fact, with Ih blocked, the parameter fitting procedure found excellent fits to the data very quickly.
CELL SELECTION CRITERIA.
Cells were accepted for analysis if their resting membrane potentials were at least 60 mV (after correction for the liquid junction potential), and their input resistance estimated from the voltage change at the end of the 200-ms current pulses was
30 M
in standard ACSF or
60 M
in ZD-7288 with no large or abrupt changes during the recording periods before and after LTP induction. Some cells that satisfied these criteria were rejected later if, after LTP induction, EPSPs were increased in the control pathway as well as in the tetanized pathway. The LTP in these latter cells tended to be larger than in cells where there was no change in EPSP size in the control pathway, suggesting instability in the slice leading to nonspecific "run-up" of EPSPs. Finally, cells were excluded if population spikes were observed during EPSPs, which occurred occasionally during posttetanus recording, because this made model calibration impossible. The criteria were satisfied in eight cells in normal ACSF and six cells exposed to ZD-7288.
Modeling procedures
CELL RECONSTRUCTIONS.
Models were constructed from reconstruction data for four CA1 pyramidal cells obtained from public databases. Specifically, we used cells n123 and n400 from the Duke/Southampton Archive of Neuronal Morphology (Cannon et al. 1998
; Pyapali et al. 1998
) at http://www.neuro.soton.ac.uk/
jchad/cellArchive/cellArchive.html, cell pc2a from Attila Gulyas website (Megias et al. 2001
) at http://www.koki.hu/
gulyas/ca1cells/index.htm, and cell c20465 from the study of Ishizuka et al. (1995)
(morphology obtained directly from the author). Cell n123 was used in models by Poirazi et al. (2003)
, and cell n400 was used by Migliore et al. (1999)
. CA1 pyramidal cells have been shown to be electrotonically similar (Mainen et al. 1996
), but we used multiple cells to take into account some of the heterogeneity among cells and their reconstructions. Cells are pictured in Fig. 1.
|
30,000 (Megias et al. 2001
NEURON LEVEL MODELING AND FITTING PROCEDURES.
Simulations were done with the NEURON simulator (Hines and Carnevale 1997
). NEURON's multiple run fitter was used to fit parameter values simultaneously to three types of recorded experimental data (20-ms, 50 pA current pulse; 200-ms, 50 pA current pulse; EPSP stimulation data). The error function was the default function in the multiple run fitter, the square norm between data points and simulation results. Differences in parameter scale were handled by using the multiple run fitter option to use the logarithms of the parameter values in the fits. These details are discussed in the APPENDIX of Shen et al. (1999)
.
Parameter values were fit for each set of experimental data with models that used each of the morphological cells shown in Fig. 1. Each fit was repeated with different starting values to see if a better fit could be obtained. Fits were done separately for data obtained before and after LTP induction. Generally fits were excellent, although in a few cases, there were problems fitting the EPSP time-course. Sometimes the best fits to the EPSP data numerically were not the best fits visually; for example, the best numerical fit sometimes had a larger or smaller peak EPSP amplitude than the data to fit the decay phase more accurately. In these cases, we restricted the fitting domain to make the fitting procedure concentrate on the more relevant parts of the EPSP or else we took the more visually appealing fit even though the error was slightly larger.
For recordings made in the presence of ZD-7288 to block Ih, model parameters fitted to the data were axial resistivity, Ra, membrane resistivity, Rm, specific membrane capacitance, Cm, and number of synapses (for pre-LTP EPSP data; also for post-LTP EPSP data as appropriate). For recordings of cells in normal ACSF, two types of fits were done that differed in regard to whether the leak conductance was uniform or nonuniform. The parameters fitted were Ra, Rm, Cm, number of synapses, and a
h scaling factor (H-factor) described below.
Cells in normal ACSF were assumed to have an Ih whose maximum conductance density,
h, increased linearly with distance from the soma (Hoffman et al. 1997
; Magee 1998
) to the distal end of the s. radiatum and was constant at more distal locations. The starting value for
h before scaling by the H-factor was 0.22 mS/cm2 at the soma, rising to 1.33 mS/cm2 at distal locations. Equations for Ih kinetics were taken from Vasilyev and Barish (2002)
and are repeated here
![]() | (1) |
![]() | (2) |
hf and
hs are the maximum conductances of the fast and slow components [fast (80%) and slow (20%) of
h] and
![]() | (3) |
q given above produce values that are 34 times faster than their reported experimental values; the resulting values are similar to those reported by Magee (1998)
Passive conductances were modeled in two different ways. First, (nonuniform leak conductance case) a uniformly distributed passive leak conductance (Rm) reversing at the resting potential was modeled together with a separate passive K-leak current whose density increased with distance from the soma. The K-leak conductance density, GKleak, was chosen to make IKleak balance Ih in all of the compartments at the resting potential so that there was no net current flow at rest. The reversal potential for IKleak was 90 mV. Fits began with initial
h and GKleak values assigned and computed as just described, and these conductance densities were scaled by a parameter called H-factor in the fitting procedure. This scaling preserved the property of a stable uniform resting potential in all compartments. In a second set of simulations (uniform leak conductance case), a single passive leak conductance (Rm) was uniformly distributed, but the reversal potential was chosen to be between the resting potential and the potassium reversal potential. Instead of fitting the passive conductancereversal potential, e_pas, directly, we initially tried fitting an L-factor (values between 0 and 1) where e_pas is defined by
![]() | (4) |
SYNAPTIC INPUTS.
Synaptic input locations were chosen from among the synapses explicitly modeled at locations shown in Fig. 1. Synapses were randomly assigned to be long-thin, mushroom-shaped, or stubby according to the proportions given by Harris et al. (1992)
. The number of open AMPA and N-methyl-D-aspartate (NMDA) receptor channels at each activated synapse as a function of time was precomputed with a synapse level model and read into NEURON from a file. The voltage at individual synapses computed with the neuron level model determined the proportion of open NMDA receptor channels at each synapse that were not blocked by magnesium. Magnesium block kinetics was taken from Ascher and Nowak (1988)
.
Because this method of modeling synaptic inputs slowed the multiple run fitter considerably (the variable step method could not be used), we used the more traditional two exponential approximation to represent the synaptic conductance in the multiple run fitter. Parameter values for this approximation were obtained from fits to the time-course of the number of open AMPA and NMDA receptor channels computed with the synapse level model. In this representation, the maximum synaptic conductance, gmax, was 187, 496, and 243 pS for the AMPA conductance and 60, 99, and 76.5 pS for the NMDA conductance at long-thin, mushroom-shaped, and stubby spines, respectively. The two time constants were 0.23 and 1.75, 0.34 and 1.65, and 0.24 and 1.75 ms for the AMPA conductance at long-thin, mushroom-shaped, and stubby spines, respectively, and were 3.5 and 40 ms for the NMDA conductance at all spines.
In fits to experimental data, we guessed at the number of activated synapses, but then let the fitting procedure scale the gmax values at the activated synapses. This scaling factor was used to update our estimated number of activated synapses, and the fits were rerun. This process was repeated iteratively until the fit returned a scaling factor very close to 1.0. We fixed the scaling factor at 1.0 and ran the fit once more. Synapse locations were chosen from a list of synapses by choosing a starting location in the list and selecting every second or third synapse in the list until the desired number was selected. Finally, we went back to using the actual precomputed numbers of open AMPA and NMDA receptor channels at each synapse for final fits. Because the peak of the soma EPSPs differed slightly when precomputed numbers of receptors were modeled instead of the two exponential approximations of the conductances used in the multiple run fitter, we had to modify the number of activated synapses and/or the starting location in our synapse list to match the EPSP peaks better. These modifications were small, and because a high degree of overlap of actual activated synapses was preserved, the EPSP shape did not change noticeably.
SYNAPSE LEVEL MODELING PROCEDURE.
The synapse level model was a deterministic model previously described (Holmes 1995
), in which a vesicle of 2,000 molecules of glutamate is released into the synaptic cleft, glutamate diffuses in the cleft, binds to AMPA and NMDA receptors on the postsynaptic membrane, and is removed from the cleft by diffusion and uptake mechanisms. AMPA receptor binding kinetics came from Jonas et al. (1993)
, with rate constants increased 1.5 times to compensate for the higher temperature of our experiments, whereas the NMDA receptor kinetics was based on results from Sather et al. (1992)
. We presented these kinetic schemes in previous models (Holmes and Levy 1997
).
Within the postsynaptic density (PSD), the density of AMPA receptors was 1,000/µm2, and the density of NMDA receptors was 300/µm2. Taking average PSD areas from Harris et al. (1992)
and assuming that PSDs are circular, the PSD radius was 0.125, 0.15, and 0.3 µm for long-thin, stubby, and mushroom-shaped spines, respectively. With these densities and PSD areas, the numbers of AMPA and NMDA receptors were 49 and 15 at long-thin spines, 71 and 21 at stubby spines, and 283 and 85 at mushroom-shaped spines. These numbers are consistent with estimates from Jonas et al. (1993)
for CA3 synapses. Single channel conductance was 8 pS for AMPA receptor channels and 50 pS for NMDA receptor channels. With these parameters, maximum AMPA conductance was 187, 243, and 496 pS for long-thin, stubby, and mushroom-shaped spines, respectively, as mentioned above. Unitary EPSPs at the soma computed with our models were consistent with the range (30665 µV) and mean values (131 µV) found by Sayer et al. (1990)
and ranges (84197 and 66275 µV) reported by Larkman et al. (1991
, 1997
).
FITTING POST-LTP EXPERIMENTAL RESULTS WITH MODELS. Post-LTP experimental data were fit with models in three different ways. First, parameter values for Rm, Ra, Cm, and H-factor, along with the number of synapses activated, were simultaneously fit to short and long current pulse data and EPSP data, as described above for pre-LTP data. Every synapse included in the pre-LTP fits was included in the post-LTP fits with additional synapses added to fit the post-LTP EPSP. Above we mentioned that, for pre-LTP fits, we chose an arbitrary starting location in our list of synapses and took every second or third synapse in the list until we had the necessary number for the fit. For post-LTP fits, additional synapses were added by going up in the list from the last pre-LTP synapse selected and also down in the list from the first pre-LTP synapse selected, selecting every second or third synapse in the list until we had the necessary number for the fit. Second, given the parameter values for Rm, Ra, Cm, and H-factor just computed for the post-LTP data, the number and location of activated synapses were fixed to those in the best pre-LTP fits. Then the AMPA single channel conductance value was increased for all channels at all activated synapses to determine the value (within 0.1 pS) that provided the best fit to the peak of the post-LTP EPSP. Third, synapse level models computed new numbers of open AMPA and NMDA receptor channels for synapses at the three types of spines for 5400% increases in receptor density (in 5% increments) while keeping PSD radius constant. Given the post-LTP parameter values computed for Rm, Ra, Cm, and H-factor plus the pre-LTP synapse number and locations, the file with increased receptor density that best fit the peak of the post-LTP EPSP was found. Because density increases were computed in 5% increments, there were some cases when one density gave an EPSP that was too small and the next gave an EPSP that was too big. In these cases, the percentage increase in receptor density chosen to match the post-LTP EPSP was estimated by interpolation.
GABAB INPUT.
We showed previously that the hyperpolarization that followed minimal synaptic stimulation was caused by a GABAB current (Holmes et al. 2006
), so we added a GABAB current at the soma to our synaptic simulations. The functional form was
![]() | (5) |
|
|
RESULTS |
|---|
|
Characteristics of cells studied
Characteristics of cells studied experimentally are shown in Table 1. Resting potential was similar both before and after LTP induction averaging 70.1 mV before and 70.2 mV after LTP induction in cells in normal ACSF and 86.3 mV before and 85.6 mV after LTP in cells treated with ZD-7288.
|
(maximum displacement) or from 66.7 to 71.7 M
(steady-state displacement) in cells in normal ACSF and from 128.9 M
to 137.1 M
in cells exposed to ZD-7288. These differences between pre-LTP and post-LTP measurements were not statistically significant in either cell group. Cells exposed to ZD-7288 responded passively to the hyperpolarizing current pulse with no voltage sag, but the voltage did not always reach a steady state by the end of the 200-ms current pulse. Because input resistance was estimated from the voltage displacement at this point, our values underestimate the actual input resistance. We did not detect any change in EPSP kinetics after LTP induction in cells in normal ACSF. The time to the peak of the EPSP at the soma averaged 11.4 ms before and 11.3 ms after LTP induction and a scaled pre-LTP EPSP usually matched the post-LTP EPSP time-course well. However, if there was a change in EPSP kinetics, it would have been very difficult to detect. If more synapses contributed to the EPSP after LTP, these might increase or reduce the time to peak observed at the soma depending on whether the added synapses were more distal or more proximal to those activated before LTP induction.
Conversely, the time to peak of the soma EPSP increased after LTP for all of the cells treated with ZD-7288. The time to peak was larger in ZD-7288 treated cells than that for the cells in ACSF, as might be expected with Ih blocked, but rose from 12.3 to 13.8 ms after LTP (P = 0.016 with a sign test). One explanation for the increase is that ZD-7288 is a slow acting blocker, and it may have blocked additional H channels during the post-LTP period. This could contribute to the increase in input resistance in these cells as well.
The percentage potentiation was quite variable among these cells, perhaps because of the small number of synapses stimulated by our near minimal stimulation paradigm. Average potentiation was 60% in cells in normal ACSF and 74% in cells with ZD-7288 added. One cell not included in these calculations was potentiated over fivefold. While such cells have been observed by others, we did not include this cell in the modeling study summary analyses.
Characteristics of the cell models calibrated with experimental data
Because the experimental cell data and the detailed morphological data used in the models came from different, although presumably representative, cells, models were calibrated to make the experimental cells and the model cells similar electrotonically. Such calibration allowed us to use actual CA1 pyramidal cell structures without having to worry about reconstruction errors or the fact that the experimental cells were different. Because of inherent variability among CA1 pyramidal cells, we did not want to rely on the morphology from a single reconstruction, so we used morphologies from four different cells in our models. Parameter value fits were done with each of the four cell models to match simultaneously three types of data from each of the 14 experimental cells both before and after LTP induction.
In general, we were able to fit the experimental data quite closely, as shown in Figs. 2 and 3, where the fits with the four morphological models typically overlap each other and the experimental data. The one difficulty we had was that fits to the EPSP for cell pc2a (and sometimes for cell n400) often had a faster time-course than the experimental data. This occurred because the current pulse data and the EPSP data pulled parameter values in opposite directions, and the best overall fit usually matched the current pulse data better than the EPSP data. Fits were equally good for cells in ACSF whether we assumed that the leak conductance was uniform or nonuniform.
|
|
|
After the fits were done, we used the fitted parameter values for each set of experimental data for each cell model to calculate the size and time-course of unitary EPSPs for each synapse activated in our simulations. Our purpose was simply to verify that unitary EPSP amplitudes in the models were similar to those reported from experiments by others. To this end, we lumped results from the four modeled cells in the results summarized in Fig. 4. As might be expected, EPSP size at the soma had a tendency to be bimodal with synapses on long-thin spines and stubby spines producing small amplitude EPSPs and synapses on mushroom-shaped spines producing much larger EPSPs (Fig. 4, A and D). Amplitudes ranged from 29 to 251 µV, with a mean of 82 µV in models of cells in ACSF and 94 µV in models of cells treated with ZD-7288. EPSP amplitudes at the location of the synapse (at spine heads) ranged from 0.2 to 14 mV, with most lying between 0.2 and 4 mV (Fig. 4, B and E).
|
Quantifying the change in the EPSP after LTP
Model fits were done to quantify the change in synaptic parameters that may occur in LTP. We quantified the number of additional activated synapses, the increase in the AMPA single channel conductance (
), and the increase in the number of receptors needed to explain the increased EPSP amplitude observed with LTP. Results are given in Table 3.
|
The increase in AMPA single channel conductance necessary to match the post-LTP EPSP amplitude was 64 and 73% in models of cells in normal ACSF with uniform and nonuniform leak conductance, respectively, and 95% in models of cells exposed to ZD-7288. For these fits, the same synapses were activated in pre-LTP and post-LTP models, with the only change being in
, the AMPA receptor single channel conductance.
The increase in the number of AMPA receptors at the synapse needed to match the post-LTP EPSP amplitude was the largest for all of the proposed LTP mechanisms, being 73 and 83% in models of cells in normal ACSF with uniform and nonuniform leak conductance, respectively, and 118% in models of cells exposed to ZD-7288. Again, the same synapses were activated in the pre-LTP and post-LTP models, but here, only the numbers of receptor channels at the synapses were changed. While the numbers of AMPA and NMDA receptors in these simulations were both increased, the voltages achieved were not sufficient for NMDA receptor channels to be open and unblocked to a degree that would affect the EPSP time-course.
We note that the percentage increases in number of activated synapses, AMPA single channel conductance, or number of receptors closely tracked the percentage increase in EPSP amplitude found experimentally after LTP when the amount of LTP was near or below the average (Fig. 5). However, when LTP magnitude was larger, the number of activated synapses, the AMPA single channel conductance, or the number of receptors had to be increased by a greater percentage than the degree of LTP found experimentally because the nonlinear properties of synaptic transmission and synaptic integration began to play a more important role. This was particularly true for increases in the number of receptors.
|
|
|
DISCUSSION |
|---|
|
In our experiments, the increase in EPSP amplitude after LTP averaged 60% for cells in normal ACSF and 74% for cells in ZD-7288. The naïve prediction is that LTP could be explained by identical percentage increases in the number of activated synapses, the AMPA single channel conductance, or the number of receptors at a synapse. We tested this with detailed neuron and synapse level models. While this naïve prediction was close to what the models found when the degree of potentiation was average or less than the average, changes at the synaptic level had to be larger than the percentage increase in the EPSP after LTP when potentiation was large. This occurred because of nonlinear properties of synaptic transmission and synaptic integration.
First, we modeled the mechanism of LTP as an increase in the number of activated synapses, as could happen with an increase in the probability of vesicle release or through unsilencing postsynaptically silent synapses. If more synapses are activated, the voltage change in the dendrites is larger, and this reduces the driving force for each synapse. Spatial summation should be less than linear requiring a greater than linear increase in the number of activated synapses to achieve the observed percentage increase in EPSP size. We saw this effect in our models when potentiation was large, but not when potentiation was average or less. In our models, activated synapses were separated by
20 µm, and a larger increase in the number of activated synapses would have been required if activated synapses had been more clustered.
Second, we modeled the mechanism of LTP as an increase in the AMPA single channel conductance. Because the increased conductance occurred at the individual synapses, voltage changes at the synapses were larger causing the driving force to be reduced directly. This should make the increment in EPSP size sublinear, requiring a greater than linear increase in single channel conductance to match the observed EPSP increase with LTP. However, the voltage change at individual synapses in our models was typically small (Fig. 4), making the effect of the reduced driving force a factor only when potentiation, and hence the increase in single channel conductance, was large.
Third, we modeled the mechanism of LTP as an increase in the number of receptors at the synapse. The increase in EPSP amplitude was sublinear in this case for two reasons. An increase in receptor number had the same effect on driving force as the increase in the AMPA single channel conductance, but in addition, because the number of glutamate molecules in a vesicle was fixed in the models, a smaller percentage of the additional receptors were activated. Receptors were not saturated by the release of 2,000 glutamate molecules in our models, and with more receptors, there was more competition for the neurotransmitter. Consequently, the percentage increases in the number of receptors at each synapse had to be larger than the percentage increases in single channel conductance or number of activated synapses to match the post-LTP EPSP (73% ACSF cells; 118% ZD-7288 cells). We note that the use of 3,000 glutamate molecules per vesicle in the models would have reduced this additional sublinear effect significantly.
Implications of the quantification for mechanisms of LTP
The numbers found by the models show that the change with LTP is quite large regardless of the mechanism. What are the implications of these numbers for proposed mechanisms of LTP?
LTP AS AN INCREASE IN VESICLE RELEASE PROBABILITY.
If LTP is caused by an increase in vesicle release probability, the increase must be of such magnitude that, for our ACSF cells, the number of activated synapses would increase from 11 to 18 (assuming a maximum of one vesicle released per synapse). We saw considerable variability in the experimental EPSP size, which could be caused by a low probability of release at individual synapses. If the probability of release averaged 0.25 before LTP, it would have had to average 0.41 after LTP. If the probability of release averaged 0.10 before LTP, this would have had to increase to 0.16 after LTP. Although most investigators do not favor this hypothesis for LTP currently, there are still reasons to consider a presynaptic component for LTP (Voronin and Cherubini 2004
), and our quantitative analysis does not reject this hypothesis because these changes in the probability of release do not seem implausible.
LTP AS AN INCREASE IN AMPA SINGLE CHANNEL CONDUCTANCE.
If LTP is caused by an increase in the AMPA single channel conductance (
), our numbers suggest that an increase from 8 to 13 pS would explain our LTP data. Benke et al. (1998)
reported a similar percentage change in
, but their initial and final values were low (4.87.2 pS). Derkach et al. (1999)
expressed homomeric GluR1 AMPA channels in HEK-293 cells and found a mean increase in
from 11.5 to 19.9 pS on exposure to calcium-calmodulindependent protein kinase II (CaMKII). This increase was the result of a higher probability of channels conducting at the 21- and 28-pS states. Poncer et al. (2002)
also observed an increase in
with
-CaMKII exposure from a mean of 11 pS in control to 16.9 pS. These observed changes are entirely consistent with our model prediction of the changes needed to explain LTP. However, these studies were done with young or juvenile animals, and it is not clear that such changes in
occur in the adult. Andrasfalvy and Magee (2004)
saw increases in the number of AMPA receptors with CaMKII activity, but no change in
in the adult.
LTP AS UNSILENCING SILENT SYNAPSES.
It has been proposed that a large percentage of synapses are silent in that they have no AMPA receptors but do have NMDA receptors and that LTP induction can cause AMPA receptors to be inserted into these synapses in minutes, particularly early in development (Isaac et al. 1995
; Liao et al.1995
, 2001
; but see Friedman et al. 2000
; Groc et al. 2002
). Estimating the percentage of silent synapses is inherently problematic, but even accepting the estimates reported, there is a limit to the number of silent synapses that can be unsilenced. The percentage of silent synapses seems to decrease rapidly during development from 2560% in 2-wk-old animals (Isaac et al. 1995
; Liao et al. 1995
; Nusser et al. 1998
; Petralia et al. 1999
) to <20% in the adult. These percentages are too low, even in the young animals in most studies, to explain the LTP we observed in adult rats. Xiao et al. (2004)
suggest that typical test frequencies might silence synapses. This could inflate the number of possible synapses that could be unsilenced later with LTP. However, this study was done with neonatal synapses, and similar results in the adult have not been reported.
There are other problems with unsilencing silent synapses as the mechanism for the LTP we observed. First, the number of AMPA receptors to be inserted in a formerly silent synapse is likely to be small, even if lipid rafts of receptors (hyperslots with 1020 receptors; Lisman 2003
) are inserted at once. This means that the number of synapses that would have to be unsilenced should be much larger than predicted by the models, because these unsilenced synapses would probably have fewer receptors than our modeled long-thin or stubby spine synapses. Large numbers of these new synapses would be needed to produce a detectable change in the somatic EPSP. Second, it is not clear where the AMPA receptors come from. Labeling studies do not show intracellular pools of AMPA receptors near silent synapses waiting to be inserted, and the amount of labeling in extrasynaptic membrane is far less than the synaptic labeling (Baude et al. 1995
; Nusser et al. 1998
; Racca et al. 2000
). It does not seem that there are enough receptors available. Our quantification suggests strongly that unsilencing silent synapses is not by itself a sufficient mechanism to explain the degree of potentiation observed with LTP.
LTP AS AN INCREASE IN AMPA RECEPTOR NUMBER.
Increasing the number of AMPA receptors at a synapse would seem to be the least efficient method for expressing potentiation. First, this reduces the driving force, and second, a smaller proportion of the additional receptors will be activated (because of the lack of receptor saturation, unless there are also presynaptic changes as well, such as a concomitant increase in the number of neurotransmitter molecules in a vesicle or the number of vesicles released). Nevertheless, it has been convincingly shown that receptors are transported to synaptic membrane during LTP (Hayashi et al. 2000
; Shi et al. 1999
).
What is not clear from the experiments is how large of an increase in AMPA receptor number occurs with LTP. Our quantification shows that the percentage increase in number of receptors has to be larger than the percentage of observed LTP and considerably larger when the degree of potentiation is large. We found that, if the EPSP amplitude were increased 60% with LTP, this could be accounted with an increase in the number of receptors on synapses on long-thin, stubby, and mushroom-shaped spines from 49, 71, and 283 to 85, 123, and 490, respectively, in our ACSF cells. In other words, 36, 52, and 207 AMPA receptors would have to be added to the three types of spine synapses. Matsuzaki et al. (2004)
suggested that LTP occurs preferentially on small synapses on long-thin and stubby spines and not on the large synapses on mushroom-shaped spines. If the increase in receptor number was to occur only on long-thin or stubby spines, the number of additional receptors per synapse predicted by our models would have to be much greater, perhaps reflecting the conversion of a long-thin or stubby spine into a mushroom-shaped one. There are reports of an increase in spine size after LTP (Harris et al. 2003
; Yuste and Bonhoeffer 2001
) that may reflect this type of conversion.
The problem with increases in the number of AMPA receptors as a mechanism for LTP is similar to that for the silent synapse hypothesis. Do synapses have such large numbers of AMPA receptors available for insertion? Again the labeling studies do not show available AMPA receptors in the required numbers, either intracellularly or extrasynaptically, to explain the observed LTP (Baude et al. 1995
; Nusser et al. 1998
; Racca et al. 2000
), particularly near the smaller spines. Lu et al. (2001)
found that the cell surface percentage of AMPA receptors in cultured CA1 pyramidal neurons rose from 58 to 73% after LTP, but unless there is recruitment of receptors from other synapses, these numbers are again too small to explain the observed LTP.
Recently, Andrasfalvy and Magee (2004)
found that patches pulled from CA1 pyramidal dendrites after LTP contained 7075% more AMPA receptors than control patches. These numbers are closer to what might be required to explain the observed LTP but would still be too small if potentiation was limited to the smaller synapses. The increase in AMPA receptors found by Andrasfalvy and Magee (2004)
was in the number of extrasynaptic receptors and was restricted to the region of the cell where LTP was induced. Interestingly, given the size of the patches, the density of AMPA receptors could have been on the order of 200300/µm2. The starting density of AMPA receptors in our models was only 45 times larger than this. Given that there is a delay before the new AMPA receptors seen at the synapses after LTP induction (Hayashi et al. 2000
; Shi et al. 1999
) and that the number of available receptors before LTP induction seems insufficient, one explanation could be that receptors are synthesized locally. If this occurs, the required number of receptors could be obtained.
To summarize, our quantification suggests that unsilencing silent synapses or increasing the numbers of AMPA receptors at synapses may not be sufficient mechanisms to explain the LTP we observed. While an increase in
or an increase in probability of release is consistent with our results, experimental work has cast doubt on these mechanisms as explanations for LTP in adult animals. Our analysis does not rule out multiple mechanisms for LTP expression or mechanisms acting at different times, and in fact suggests that some combination of mechanisms is most probable. In particular, the very large values of potentiation that are observed on occasion (e.g., our cell 081704b, see Table 3; Debanne et al. 1999
) are inconsistent with a single mechanism. It has also been shown that different mechanisms are used in different proportions in different stages of development (Palmer et al. 2004
).
Notes on the models
The results of this study and its implications for the various proposed hypotheses about the mechanism of LTP depend on a number of modeling assumptions that deserve comment.
RECEPTOR NUMBERS.
The number of receptors we chose for the three types of synapses is based on reported PSD size and an assumed density of AMPA receptors. These numbers of receptors produce EPSPs at the soma in our calibrated models of actual morphology that are comparable in size with experimental estimates (Andersen 1982
, 1990
; Jonas et al. 1993
; Larkman et al. 1991
, 1997
; Sayer et al. 1989
, 1990
). The peak number of open AMPA receptor channels is 23, 30, and 62 for long-thin, stubby, and mushroom-shaped spine synapses, respectively. Others have proposed higher estimates for the synaptic conductance than we used in the models (Forti et al. 1997
; Markram et al. 1997
; Stricker et al. 1996
). In these studies, it seems likely that either the methods used are biased toward selection of particularly large synapses or there is some uncertainty about estimates of the number of synapses that may have contributed to a unitary EPSP. Conversely, Nimchinsky et al. (2004)
suggested that <10 AMPA channels are open at the peak, which is much smaller than in the models. This number would produce EPSPs in our calibrated models that would generally be 530 µV at the soma if AMPA single channel conductance is 8 pS. AMPA receptor number has been shown to be quite variable among synapses, whereas NMDA receptor number seems more stable, so it is possible that we biased our synapse types to those that would produce EPSPs detectable experimentally. However, the relatively small reported percentages of silent synapses in the adult suggest that our numbers are plausible. Our NMDA channel numbers produce one to two open NMDA receptor channels at the peak, in agreement with Nimchinsky et al. (2004)
, but these are blocked by magnesium at the potentials generated in the models.
CALIBRATED MODELS.
Typically, modeling studies use the morphology of a single cell. We used morphological data from four cells in this study. One reason is that we do not have reconstructions for the specific cells that were used in the physiological experiments, and we felt that four different cells might represent some of the heterogeneity in CA1 cells. However, second, even if we had morphological information for each of our experimental cells, measurements in morphological reconstructions are subject to a number of difficulties that can affect model results (Ambros-Ingerson and Holmes 2005
). These difficulties are manifested in the large variability in cell statistics listed in Fig. 1. For this reason, we believe that for modeling results to be robust, simulations should be done with several reconstructions of the same cell type, as we have done here. We argue elsewhere (Holmes et al. 2006
) that model calibration can overcome inherent problems with cell heterogeneity and reconstructions to make modeled cells similar to experimental cells electronically while preserving the basic structure of the cell type modeled. While such calibration can lead to model parameter values that are not traditional, as in Table 2, we believe this is the correct approach to use.
UNIQUENESS. Whenever one does fits of the type done here, it is always a concern whether or not fitted parameter values provide a unique solution. Here we presented two ways to represent the leak conductance and could not distinguish between them in the fits. However, once we settled on a means to represent the leak conductance, assigned Ra and Cm to be constant and fixed the functional form for the H conductance distribution, our final results were very consistent with a wide range of starting values. In no case did we find highly different sets of parameter values that provided equally good (best) fits to the data. Very different parameter values can give fits that may not look so bad visually, but these fits had an error several-fold larger than the best fit and were never the final fit returned by the fitting procedure. Our fits were done with averaged data and one might argue that perhaps averaged data are not truly representative or that it includes features not characteristic of the actual system. If true, the very good fits we obtained might be fitting these noncharacteristic features. We cannot guarantee that our fits are unique or that slight changes in the experimental data points would not cause very different parameter values in the fits. However, our fitted parameter values were quite consistent within a morphological model across experimental cells with the relatively few parameters we included in the fits, suggesting that a possible lack of uniqueness was not an issue for our results.
|
|
GRANTS |
|---|
|
|
|
ACKNOWLEDGMENTS |
|---|
|
|
|
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: W. R. Holmes, Neuroscience Program, Dept. of Biological Sciences, Ohio Univ., Athens, OH 45701 (E-mail: holmes{at}ohio.edu)
|
|
REFERENCES |
|---|
|
Andersen P. Synaptic control and modulation of hippocampal neurons. Adv Biol Psych 8: 4055, 1982.
Andersen P. Synaptic integration in CA1 hippocampal pyramidal cells. Prog Brain Res 83: 215222, 1990.[Web of Science][Medline]
Andrasfalvy BK and Magee JC. Changes in AMPA receptor currents following LTP induction on rat CA1 pyramidal neurons. J Physiol 559: 543554, 2004.
Ascher P and Nowak L. The role of divalent cations in N-methyl-D-aspartate responses of mouse central neurons in culture. J Physiol 399: 247266, 1988.
Baude A, Nusser Z, Molnar E, McIlhinney RAJ, and Somogyi P. High-resolution immunoigold localization of AMPA type glutamate receptor subunits at synaptic and non-synaptic sites in rat hippocampus. Neuroscience 69: 10311055, 1995.[CrossRef][Web of Science][Medline]
Benke TA, Luthi A, Isaac JT, and Collingridge GL. Modulation of AMPA receptor unitary conductance by synaptic activity. Nature 393: 793797, 1998.[CrossRef][Medline]
Blanton MG, Lo Turco J, and Kriegstein AR. Whole cell recording from neurons in slices of reptilian and mammalian cerebral cortex. J Neurosci Methods 30: 203210, 1989.[CrossRef][Web of Science][Medline]
Bliss TVP and Collingridge GL. A synaptic model of memory: long-term potentiation in the hippocampus. Nature 361: 3139, 1993.[CrossRef][Medline]
Cannon RC, Turner DA, Pyapali GK, and Wheal HV. An online archive of reconstructed hippocampal neurons. J Neurosci Methods 4: 4954, 1998.
Debanne D, Gahwiler BH, and Thompson SM. Heterogeneity of synaptic plasticity at unitary CA3-CA1 and CA3-CA3 connections in rat hippocampal slice cultures. J Neurosci 19: 1066410671, 1999.
Derkach V, Barria A, and Soderling TR. Ca2+/calmodulin-kinase II enhances channel conductance of
-amino-3-hydroxy-5-methyl-4-isoxazolepropionate type glutamate receptors. Proc Natl Acad Sci USA 96: 32693274, 1999.
Engert F and Bonhoeffer T. Dendritic spine changes associated with hippocampal long-term synaptic plasticity. Nature 399: 6670, 1999.[CrossRef][Medline]
Forti L, Bossi M, Bergamaschi A, Villa A, and Malgaroll A. Loose-patch recordings of single quanta at individual hippocampal synapses. Nature 388: 874878, 1997.[CrossRef][Medline]
Friedlander MJ, Sayer RJ, and Redman SJ. Evaluation of long-term potentiation of small compound and unitary EPSPs at the hippocampal CA3-CA1 synapse. J Neurosci 10: 814825, 1990.[Abstract]
Friedman HV, Bresler T, Garner CC, and Ziv NE. Assembly of new individual excitatory synapses: time course and temporal order of synaptic molecule recruitment. Neuron 27: 5769, 2000.[CrossRef][Web of Science][Medline]
Gasparini S and DiFrancesco D. Action of the hyperpolarization-activated current (Ih) blocker ZD 7288 in hippocampal CA1 neurons. Pflugers 435: 99106, 1997.[CrossRef]
Groc L, Gustafsson B, and Hanse E. Spontaneous unitary synaptic activity in CA1 pyramidal neurons during early postnatal development: constant contribution of AMPA and NMDA receptors. J Neurosci 22: 55525562, 2002.
Grover LM. Evidence for postsynaptic induction and expression of NMDA receptor independent LTP. J Neurophysiol 79: 11671182, 1998.
Grover LM and Chen Y. Evidence for involvement of group II/III metabotropic glutamate receptors in NMDA receptor-independent long-term potentiation. J Neurophysiol 82: 29562969, 1999.
Gustafsson B, Wigsrom H, Abraham WC, and Huang Y-Y. Long-term potentiation in the hippocampus using depolarizing current pulses as the conditioning stimulus to single volley of synaptic potentials. J Neurosci 7: 774780, 1987.[Abstract]
Harris KM, Fiala JC, and Ostroff L. Structural changes at dendritic spine synapses during long-term potentiation. Philos Trans R Soc Lond B Biol Sci 358: 745748, 2003.
Harris KM, Jensen FE, and Tsao B. Three-dimensional structure of dendritic spines and synapses in rat hippocampus (CA1) at postnatal day 15 and adult ages: implications for the maturation of synaptic physiology and long-term potentiation. J Neurosci 12: 26852705, 1992.[Abstract]
Harris NC and Constanti A. Mechanism of block by ZD 7288 of the hyperpolarization-activated inward rectifying current in guinea pig substantia nigra neurons in vitro. J Neurophysiol 74: 23662378, 1995.
Hayashi Y, Shi SH, Esteban JA, Piccini A, Poncer JC, and Malinow R. Driving AMPA receptors into synapses by LTP and CaMKII: requirement for GluR1 and PDZ domain interaction. Science 287: 22622267, 2000.
Hines M and Carnevale NT. The NEURON simulation environment. Neural Comput 9: 11791209, 1997.[CrossRef][Web of Science][Medline]
Hoffman DA, Magee JC, Colbert CM, and Johnston D. K+ channel regulation of signal propagation in dendrites of hippocampal pyramidal neurons. Nature 387: 869875, 1997.[CrossRef][Medline]
Holmes WR. Modeling the effect of diffusion and uptake on NMDA and non-NMDA receptor saturation. Biophys J 69: 17341747, 1995.[Web of Science][Medline]
Holmes WR, Ambros-Ingerson J, and Grover LM. Fitting experimental data to models that use morphological data from public databases. J Comput Neurosci 20: 349365, 2006.[CrossRef][Web of Science][Medline]
Holmes WR and Levy WB. Quantifying the role of inhibition in associative long-term potentiation in dentate granule cells with computational models. J Neurophysiol 78: 103116, 1997.
Isaac JTR, Hjelmstad GO, Nicoll RA, and Malenka RC. Long-term potentiation at single fiber inputs to hippocampal CA1 pyramidal cells. Proc Natl Acad Sci USA 93: 87108715, 1996.
Isaac JTR, Nicoll RA, and Malenka RC. Evidence for silent synapses: implications for the expression of LTP. Neuron 15: 427434, 1995.[CrossRef][Web of Science][Medline]
Ishizuka N, Cowan WM, and Amaryl DG. A quantitative analysis of the dendritic organization of pyramidal cells in the rat hippocampus. J Comp Neurol 362: 1745, 1995.[CrossRef][Web of Science][Medline]
Jonas P, Major G, and Sakmann B. Quantal components of unitary EPSCs at the mossy fibre synapse on CA3 pyramidal cells of rat hippocampus. J Physiol 472: 615663, 1993.[Web of Science][Medline]
Larkman A, Stratford K, and Jack J. Quantal analysis of excitatory synaptic action and depression in hippocampal slices. Nature 350: 344347, 1991.[CrossRef][Medline]
Larkman AU, Jack JJB, and Stratford KJ. Quantal analysis of excitatory synapses in rat hippocampal CA1 in vitro during low-frequency depression. J Physiol 505: 457471, 1997.
Liao D, Hessler NA, and Malinow R. Activation of postsynaptically silent synapses during pairing-induced LTP in CA1 region of hippocampal slice. Nature 375: 400404, 1995.[CrossRef][Medline]
Liao D, Scannevin RH, and Huganir R. Activation of silent synapses by rapid activity-dependent synaptic recruitment of AMPA receptors. J Neurosci 21: 60086017, 2001.
Lisman J. Long-term potentiation: outstanding questions and attempted synthesis. Philos Trans R Soc Lond B Biol Sci 358: 829842, 2003.
Lu W-Y, Man H-Y, Ju W, Trinble WS, MacDonald JF, and Wang YT. Activation of synaptic NMDA receptors and LTP in cultured hippocampal neurons. Neuron 29: 243254, 2001.[CrossRef][Web of Science][Medline]
Maccaferri G, Mangoni M, Lazzari A, and DiFrancesco D. Properties of the hyperpolarization-activated current in rat hippocampal CA1 pyramidal cells. J Neurophysiol 69: 21292136, 1993.
Magee JC. Dendritic hyperpolarization-activated currents modify the integrative properties of hippocampal CA1 pyramidal neurons. J Neurosci 18: 76137624, 1998.
Mainen ZF, Carnevale NT, Zador AM, Claiborne BJ, and Brown TH. Electrotonic architecture of hippocampal CA1 pyramidal neurons based on three-dimensional reconstructions. J Neurophysiol 76: 19041923, 1996.
Malenka RC and Nicoll RA. Long-term potentiationa decade of progress? Science 285: 18701874, 1999.
Markram H, Lubke J, Frotscher M, Roth A, and Sackmann B. Physiology and anatomy of synaptic connections between thick tufted pyramidal neurons in the developing rat neocortex. J Physiol 500: 409440, 1997.
Matsuzaki M, Honkura N, Ellis-Davis GCR, and Kasai H. Structural basis of long-term potentiation in single dendritic spines. Nature 429: 761766, 2004.[CrossRef][Medline]
Megias M, Emri A, Freund TF, and Gulyas AI. Total number and distribution of inhibitory and excitatory synapses on hippocampal CA1 pyramidal cells. Neuroscience 102: 527540, 2001.[CrossRef][Web of Science][Medline]
Migliore M, Hoffman DA, Magee JC, and Johnston D. Role of an A-type K+ conductance in the back-propagation of action potentials in the dendrites of hippocampal pyramidal neurons. J Comput Neurosci 7: 515, 1999.[CrossRef][Web of Science][Medline]
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.
Nusser A, Lujan R, Laube G, Roberts JDB, Molnar E, and Somogyi P. Cell type and pathway dependence of synaptic AMPA receptor number and variability in the hippocampus. Neuron 21: 545559, 1998.[CrossRef][Web of Science][Medline]
Otis TS, DeKoninck Y, and Mody I. Characterization of synaptically elicited GABAB responses using patch-clamp recordings in rat hippocampus slices. J Physiol 463: 391407, 1993.
Palmer MJ, Isaac JTR, and Collingridge GL. Multiple, developmentally regulated expression mechanisms of long-term potentiation at CA1 synapses. J Neurosci 24: 49034911, 2004.
Petersen CCH, Malenka RC, Nicoll RA, and Hopfield JJ. All-or-none potentiation at CA3-CA1 synapses. Proc Natl Acad Sci USA 95: 47324737, 1998.
Petralia RS, Esteban JA, Wang Y-X, Partridge JG, Zhao H-M, Wenthold RJ, and Malinow R. Selective acquisition of AMPA receptors over postnatal development suggests a molecular basis for silent synapses. Nat Neurosci 2: 3136, 1999.[CrossRef][Web of Science][Medline]
Poirazi P, Brannon T, and Mel BW. Arithmetic of subthreshold synaptic summation in a model CA1 pyramidal cell. Neuron 37: 977987, 2003.[CrossRef][Web of Science][Medline]
Poncer JC, Esteban JA, and Malinow R. Multiple mechanisms for the potentiation of AMPA receptor-mediated transmission by
-Ca2+/calmodulin-dependent protein kinase II. J Neurosci 22: 44064411, 2002.
Pyapali GK, Sik A, Penttonen M, Buzsaki G, and Turner DA. Dendritic properties of hippocampal CA1 pyramidal neurons in the rat: intracellular staining in vivo and in vitro. J Comp Neurol 391: 335352, 1998.[CrossRef][Web of Science][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.
Rall W, Burke RE, Holmes WR, Jack JJB, Redman SJ, and Segev I. Matching dendritic neuron models to experimental data. Physiol Rev 72: S159S186, 1992.[Web of Science][Medline]
Sastry BR, Goh JW, and Auyeung A. Associative induction of posttetanic and long-term potentiation in CA1 neurons of rat hippocampus. Science 232: 988990, 1986.
Sather W, Dieudonne S, MacDonald JF, and Ascher P. Activation and desensitization of N-methyl-D-aspartate receptors in nucleated outside-out patches from mouse neurons. J Physiol 450: 643672, 1992.
Sayer RJ, Friedlander MJ, and Redman SJ. The time course and amplitude of EPSPs evoked at synapses between pairs of CA3/CAl neurons in the hippocampal slice. J Neurosci 10: 826836, 1990.[Abstract]
Sayer RJ, Redman SJ, and Andersen P. Amplitude fluctuations in small EPSPs recorded from CA1 pyramidal cells in the guinea pig hippocampal slice. J Neurosci 9: 840850, 1989.[Abstract]
Shen GY, Chen WR, Midtgaard J, Shepherd GM, and Hines ML. Computational analysis of action potential initiation in mitral cell soma and dendrites based on dual patch recordings. J Neurophysiol 82: 30063020, 1999.
Shi SH, Hayashi Y, Petralia RS, Zaman SH, Wenthold RJ, Svoboda K, and Malinow R. Rapid spine delivery and redistribution of AMPA receptors after synaptic NMDA receptor activation. Science 284: 18111816, 1999.
Stevens CF and Wang Y. Changes in reliability of synaptic function as a mechanism for plasticity. Nature 371: 704707, 1994.[CrossRef][Medline]
Stricker C, Field AC, and Redman SJ. Statistical analysis of amplitude fluctuations in EPSCs evoked in rat CA1 pyramidal neurons. in vitro. J Physiol 490: 419441, 1996.
Vasilyev DV and Barish ME. Postnatal development of hyperpolarization-activated excitatory current Ih in mouse hippocampal pyramidal neurons. J Neurosci 22: 89929004, 2002.
Voronin LL and Cherubini E. 'Deaf, mute and whispering silent synapses: their role in synaptic plasticity. J Physiol 557: 312, 2004.
Xiao M-Y, Wasling P, Hanse E, and Gustafsson B. Creation of AMPA-silent synapses in the neonatal hippocampus. Nat Neurosci 7: 236243, 2004.[CrossRef][Web of Science][Medline]
Yuste R and Bonhoeffer T. Morphological changes in dendritic spines associated with long-term synaptic plasticity. Annu Rev Neurosci 24: 10711089, 2001.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
L. M. Grover, E. Kim, J. D. Cooke, and W. R. Holmes LTP in hippocampal area CA1 is induced by burst stimulation over a broad frequency range centered around delta Learn. Mem., January 7, 2009; 16(1): 69 - 81. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. K. Bichler, S. T. Nakanishi, Q.-B. Wang, M. J. Pinter, M. M. Rich, and T. C. Cope Enhanced Transmission at a Spinal Synapse Triggered In Vivo by an Injury Signal Independent of Altered Synaptic Activity J. Neurosci., November 21, 2007; 27(47): 12851 - 12859. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. A. Earnshaw and P. C. Bressloff Biophysical Model of AMPA Receptor Trafficking and Its Regulation during Long-Term Potentiation/Long-Term Depression. J. Neurosci., November 22, 2006; 26(47): 12362 - 12373. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |