Abstract
Plateau potentials are a critical feature of neuronal excitability, but their all-or-none behavior is not easily captured in modeling. In this study, we investigated models of plateau potentials in multi-compartment neuron models and found that including glutamate spillover provides robust all-or-none behavior. This result arises due to the prolonged duration of extrasynaptic glutamate. When glutamate spillover is not included, the all-or-none behavior is very sensitive to the steepness of the Mg2+ block. These results suggest a potentially significant role of glutamate spillover in plateau potential generation, providing a mechanism for robust all-or-none behavior across a wide range of slopes of the Mg2+ block curve. We also illustrate the importance of the all-or-none plateau potential behavior for nonlinear computation with regard to the nonlinear feature binding problem.
1. Introduction
Plateau potentials are all-or-none voltage elevations that occur in neuronal dendrites in various brain regions, including striatal projection neurons (SPNs), which last tens to hundreds of milliseconds (Plotkin et al., ; Oikonomou et al., ; Du et al., ). They are involved in processing sensory information, as well as triggering synaptic plasticity (Lavzin et al., ; Xu et al., 2012; Gambino et al., ; Kumar et al., ). Clustered synapses, which facilitate the induction of dendritic NMDA spikes and plateau potentials, have been shown to form during development and learning (see e.g., review by Kastellakis and Poirazi, ). Plateau potentials are generated on stretches of dendrites around and above 20 μm of length, where the N-methyl-D-aspartate (NMDA) receptors that mediate them must have a cumulative conductance above a threshold level (Antic et al., ). Experimentally, plateau potentials have been evoked by activating clusters of synapses by glutamate uncaging or by stimulating stretches of dendrites by glutamate iontophoresis or repeated synaptic stimulation (Milojkovic et al., ; Major et al., ; Plotkin et al., ; Oikonomou et al., ; Du et al., ; Kumar et al., ; Gao et al., ). In the cases of progressively increasing glutamate stimulus in equal increments, a supralinear, all-or-none response in the amplitude of the somatic voltage is observed (Major et al., ; Oikonomou et al., ; Gao et al., ). In this study, we deal with the question of how to model the all-or-none behavior of plateau potentials in multicompartment neuron models and discuss the implications of the model for the precise mechanism of generation of plateau potentials experimentally and in vivo.
The NMDA receptors are ionotropic receptors activated by glutamate and glycine as a co-factor that are additionally blocked by Mg2+ ions at resting membrane potentials. The amount of blockage depends on the local voltage, and it has been determined that this dependence is a sigmoidal nonlinearity (Figure 1C), where higher voltage relieves the block (Jahr and Stevens, ). The generation mechanism is a positive feedback loop that starts by partially alleviating the Mg2+ block on NMDA receptors (NMDARs) when the dendritic voltage reaches a threshold value, usually around −50 mV. This causes the flow of an inward current through the NMDA receptors, raising the voltage in a cycle that progressively alleviates more of the Mg2+ block and allows for more inward current. Thus, a sustained voltage elevation is produced, which lasts as long as there is glutamate to bind to the NMDARs. The sigmoidal shape of the voltage nonlinearity is a crucial ingredient for the all-or-none behavior of the plateau potentials since, for linear increases in voltage, it provides a “sudden” relief of Mg2+ ions from the NMDARs. The binding of neurotransmitters to the NMDAR initiates a transition from the closed receptor state to the open state. The behavior of an NMDAR in the open state is described with the diagram in Figure 1A, also called the three-state model, since it includes three receptor states: open, blocked, and closed. An NMDAR in the open state can transition to two different non-conducting states, closed and blocked. The transitions between the open and blocked states depend on Mg2+ concentration (denoted as [Mg2+]) and voltage, while the transitions from the open and blocked states to the closed state are voltage- and [Mg2+]-independent. Voltage crossing the −50 mV threshold unblocks some of the receptors, initiating the self-sustaining positive feedback loop. Neurotransmitter removal closes the receptor pore (the receptor transitions to the closed state), thus stopping the inward current and lowering the voltage. Lowering the voltage causes remaining open receptors (with bound neurotransmitter) to transition to the blocked state, which because of the sigmoidal nonlinearity rapidly lowers the voltage back to the resting state.
Figure 1
Table 1
| Transition rate | Value |
|---|---|
| a1 | exp(−0.016·V−2.91) [ms−1] |
| a2 | C · exp(−0.045·V−6.97) [μM−1 ms−1] |
| b1 | exp(0.009·V+1.22) [ms−1] |
| b2 | exp(0.017·V+0.96) [ms−1] |
| A | exp(−2.847) [ms−1] |
| B1 | exp(−0.693) [ms−1] |
| B2 | exp(−3.101) [ms−1] |
Transition rates for the four-state model reproduced from Jahr and Stevens (,).
The transition rates were derived by single-channel analysis. C represents the magnesium concentration in μM and V the voltage in mV.
1.1. NMDAR state diagrams and the gating function
The three-state model has been shown to be suitable for high Mg2+ concentrations ([Mg2+]>0.2 mM), but it can be easily expanded to take into account lower concentrations by introducing an additional, [Mg2+]-independent, blocked state (Figure 1B, the four-state model; Jahr and Stevens, ). All the transition rates of the four-state model have been derived from single-channel measurements and reported in Table 1 (reproduced from Jahr and Stevens, ). The rate of leaving the open state in the three-state model (O→B) is a, and the corresponding rate for the four-state model (B1←O→B2) is a1+a2. The transition rate from blocked to open is b in the three-state model, whereas the corresponding rate in the four-state model is b1a1/(a1+a2)+b2a2/(a1+a2) because the rate of leaving each blocked state is weighted by the probability that particular state had been entered. Similarly, the closing rate from the blocked state is B for the three-state model and B1a1/(a1+a2)+B2a2/(a1+a2). The Mg2+ block is expressed with a gating function, g(V), which gives the fraction of open NMDARs for a given voltage value. As shown in Jahr and Stevens () it is possible to express g(V) in terms of the transition rates as:
The single-channel analysis showed that b1 and b2 are much larger than B1 and B2 (also visible in Table 1); hence, the following is a good approximation:
Moreover, when [Mg2+] is greater than a few hundred micromolar, a2 becomes much larger than a1 since the former increases linearly with [Mg2+] while the latter is independent of it, and hence g(V) for physiological [Mg2+] can be further approximated as:
Finally, taking into account the transition rates in Table 1, the gating function can be written in its most commonly used form:
where α is the steepness of the curve, η is related to the voltage value V1/2 where half of the receptors are free of Mg2+ according to , [Mg2+] is the concentration of extracellular Mg2+ in mM and V is the membrane potential in mV (see Supplementary material for details and Supplementary Figure 1). As described in the Supplementary material, the values of α and η are determined by the transition rates between different states in the 4-state model, which in turn are determined by the structural properties of the NMDARs.
The gating function adequately describes the behavior not only of single channels but also of macroscopic NMDAR currents, which allows for the determination of the parameters α and η from macroscopic current measurements. Such macroscopic current measurements from NMDARs, as well as gating functions fitted to such measurements, have been reported numerous times since the work of Jahr and Stevens (), and are summarized in Table 2. For completeness, in Table 2 we have added parameter values which were not estimated in the original articles. The fitting procedure, as well as the results of the fittings, are described in the Methods and Results sections, respectively. Table 2 shows a natural variability in the parameters of the Mg2+ block gating function which depends on the NMDAR subunit composition, as well as on the concentration of other intra- and extra-cellular ions (Qian et al., ; McMenimen et al., ; Qian and Johnson, ; Retchless et al., ). Since the sigmoidal Mg2+ block is a crucial ingredient for the all-or-none behavior of plateau potentials, we investigate the effect of the parameter variability on the behavior of the plateau potentials.
Table 2
| NMDAR type | [Mg] (mM) | α (/mV) | η (/mM) | Cell type | Ionic concentration (mM) | |
|---|---|---|---|---|---|---|
| Nowak et al. () | ? | 0.5 | 0.04 | 1.33 | Mouse mesencephalic or striatal neurons | 140 [Na]o |
| Jahr and Stevens () | ? | 1 | 0.062 | 0.28 | Rat hippocampal neurons | 165 [Na]o / 150 [Cs]i |
| Chen and Huang () | ? | 0.03 | 0.05 | 0.49 | Trigeminal subnucleus caudalis neurons | 140 [Na]o / 125[Cs]i |
| Sharma and Stevens (1996) | N1/2A | 3 | 0.06 | 0.28 | HEK293 cells | 145 [Na]o / 160[Cs]i |
| McMenimen et al. () | N1/2B | 2 | 0.06 | 0.42 | Xenopus laevis oocytes | 96 [Na]o |
| N1/2B | 0.2 | 0.05 | 3.3 | 96 [Na]o | ||
| Chiu and Carter () | ? | 1 | 0.074 | 0.11 | CD1 mice | |
| Layer 2/3 pyramidal neurons | 155 [Na]o, 2 [Ca]o, 2.5 [K]o / 130 [Cs]i or 128 [K]i | |||||
| ? | 0.7 | 0.074 | 0.104 | 155 [Na]o, 1.2 [Ca]o, 4.2 [K]o / 130 [Cs]i or 128 [K]i | ||
| ? | 0.8 | 0.071 | 0.119 | 155 [Na]o, 1.3 [Ca]o, 3.8 [K]o / 130 [Cs]i or 128 [K]i |
Gating function parameter values from experimental data.
Parameter values are determined from fits to examples of experimental data. Fits for all articles except Jahr and Stevens () were done as described in Section 2.
1.2. Computational models of dendritic plateau potentials
Studies that model plateau potentials and/or NMDA spikes fall into two groups. In the first group are studies that simply include the voltage-dependent Mg2+ block of the NMDARs in an NMDA synapse model and activate a cluster of such synapses on a dendrite to generate the plateau potential. The second group comprises one recent study that includes the effect of glutamate spillover when modeling plateau potentials (Gao et al., ). Glutamate spillover is thought to occur when a cluster of synapses is stimulated more strongly or repeatedly, and the amounts of exocytosed glutamate surpass the ability of astrocytes to take it up, thus causing it to spill over from the synaptic cleft into the space surrounding the spine necks and dendritic shafts. Here it stimulates extrasynaptic NMDARs (eNMDARs), which can aid the generation of plateau potentials (Chalifoux and Carter, ). The eNMDARs have also been found to have a substantial contribution to neuronal up-states in striatal SPNs, which in turn have been hypothesized to be driven by dendritic plateau potentials (Oikonomou et al., ; Garcia-Munoz et al., ).
None of the studies have explicitly focused on modeling the all-or-none behavior of the plateau potentials. However, all studies have changed the originally reported Mg2+ block parameters by Jahr and Stevens () without explicitly stating a reason for the change. The parameter values for α, η, and [Mg2+] that have been used in the different studies are exemplified in Table 3. In addition, all of these studies except one use steeper sigmoid curves compared to the one reported by Jahr and Stevens ().
Table 3
| [Mg] (mM) | α (/mV) | η (/mM) | |
|---|---|---|---|
| Rhodes (2006) | 1 and 2 | 0.08 | 0.28 |
| Major et al. () | 1.8 | 0.08 | 0.11 |
| Farinella et al. () | 1 | 0.08 | 0.3 |
| Poleg-Polsky () | 1 | 0.08 | 0.25 |
| Doron et al. () | 1 | 0.08 | 0.28 and 1.45 |
| Du et al. () | 1 | 0.07 | 0.33 |
| Dorman et al. () | 1.4 | 0.099 | 0.055 |
| Kumar et al. () | 1 | 0.08 | 0.25 |
| Ecker et al. () | 1 | 0.062 | 0.38 |
| Gao et al. () | 1 | 0.08 | 0.25 |
Examples of studies that have modeled plateau potentials and the parameters used for the Mg2+ block.
Ecker et al. () have corrected η for the omitted junction potential in Jahr and Stevens (), and use the same value for α. Note that Major et al. () have used a different form of the equation, and we have converted their parameter values to match the form in Equation (4).
We analyze the two different types of models of generating plateau potentials and find that the all-or-none behavior is very sensitive to the steepness of the Mg2+ block when no glutamate spillover is included in the model and, conversely, very robust with respect to the steepness of the Mg2+ block when glutamate spillover is included. Without glutamate spillover, in order to obtain the all-or-none behavior of the plateau potentials, one needs to increase the steepness of the Mg2+ block since the originally reported curve is usually too shallow and provides graded rather than all-or-none voltage elevations. In light of the variability in the steepness of the gating function (Table 2), these results indicate that glutamate spillover might be the reason why plateau potentials exhibit a robust all-or-none quality in experimental recordings.
A robust all-or-none plateau potential behavior can significantly improve nonlinear computation, as we illustrate for the nonlinear feature binding problem. We conclude by discussing the implications of the glutamate spillover model in light of the many NMDAR isoforms, the experimentally observed all-or-none behavior of plateau potentials, as well as the predicted role for dendritic computations.
2. Methods
The code implementation is done in Python 3, and the simulations in NEURON (Carnevale and Hines, ) are used via the Python interface.
2.1. Neuron models
The SPN model is taken from the collection of models published in Lindroos and Hellgren Kotaleski (). Spines in the neurons are added as additional compartments consisting of a neck and head with lengths and diameters lneck = 0.5 μm, lhead = 0.5 μm, and dneck = 0.125 μm, dhead = 0.5 μm, respectively. Axial resistance in all compartments is 150 Ohm · cm, except for the spine neck, where it is 1, 130 Ohm · cm (Dorman et al., ). Spines contain the inwardly rectifying potassium channel, with a conductance equal to that of the parent dendritic shaft segment. Voltage-gated calcium channels of types R (Cav2.3), T (Cav3.2 and Cav3.3), and L (Cav1.2 and Cav1.3) are added to the spines as well, and their conductances have been manually tuned to match the relative proportions determined in Carter and Sabatini (), and Higley and Sabatini (), as well as calcium concentration amplitudes arising from backpropagating action potential (bAP) stimulation as in Figure 2 of Shindou et al. (2011).
2.2. Background synaptic noise
Background synaptic noise is implemented by adding synapses across the dendrites according to the excitatory and inhibitory synapse densities reported in Cheng et al. (). These synapses are not on spines. We only explicitly model spines where we have clustered synaptic inputs. To increase the speed of the simulations, in compartments where the synaptic density specifies adding more than one synapse, only one synapse was added with scaled-up input frequency to account for the actual number of synapses that would need to arrive in that compartment.
2.3. Synaptic inputs
Synaptic inputs consist of AMPA and NMDA synapses located on spines within clusters. Glutamate spillover is modeled by including extrasynaptic NMDA conductances on the dendritic shafts situated directly beneath the spines. For glutamate spillover, we put as many extrasynaptic NMDA synapses as the NMDA synapses on the spines. The AMPA and NMDA synapse models are taken from Gao et al. (), which are a variation of the saturating synapse models in Destexhe et al. () implemented in the NEURON simulator (Carnevale and Hines, ). These are kinetic models which operate according to two different kinetic schemes depending on the presence of neurotransmitters. When neurotransmitters arrive at the postsynaptic site, receptor dynamics are described by a kinetic scheme that models switching between closed (C) and open (O) states with rates of α and β:
where T is the transmitter concentration. After neurotransmitter has been cleared from the synaptic cleft, receptor dynamics evolves according to a kinetic scheme that describes just the closing of the receptor with a rate β:
When a presynaptic spike arrives, neurotransmitter levels are assumed to always reach a fixed saturating concentration, Tmax, in the synaptic cleft, i.e., are represented by a pulse with amplitude Tmax and duration of Tdur. A presynaptic spike that arrives while the neurotransmitter pulse is still on lengthens the pulse duration by Tdur.
Extrasynaptic receptors are modeled by the same kinetic scheme, with two differences: (i) Tdur is much longer, modeling the spillover effect of glutamate in the extrasynaptic space, and (ii) the synaptic spike there arrives with a delay of 1.5 ms after the glutamate threshold has been reached (taken from Szapiro and Barbour, 2007, as we describe in the next section). The channel properties of NMDARs and eNMDARs have been reported to be similar (Clark et al., ).
We use the model given in Gao et al. () but with modified parameter values, given in Table 4. For example, in the simulation code for the model in Gao et al. () different maximal transmitter concentration, Tmax, for the synaptic AMPA and NMDA synapses is used. Although difficult to justify physiologically, this difference would not affect the results in Gao et al. (), since the AMPA synapses are only used to activate the NMDARs and eNMDARs, and since the effect of the different concentrations can be offset with suitable values for the synaptic conductances. We have nevertheless used the same Tmax for both AMPARs and synaptic NMDARs. In addition, the model in Gao et al. () uses the same Tmax for NMDARs and eNMDARs. We have changed this, since, even though the extrasynaptic glutamate concentration is not precisely known, extrasynaptic spillover-induced currents are very sensitive to low-affinity AMPA antagonists, suggesting lower extrasynaptic glutamate concentrations (Szapiro and Barbour, 2007). Low glutamate concentrations during spillover have also been reported in Okubo et al. (). As for the maximal synaptic glutamate concentration, Tmax, it has been found to be very close to 1 mM at cultured hippocampal synapses (Clements et al., ).
Table 4
| gmax (nS) | w | Tdur (ms) | Tmax (mM) | α [/(ms mM)] | β (/ms) | Kd (μM) (nS) | |
|---|---|---|---|---|---|---|---|
| AMPA | 1.5 | 0.4 | 1 | 1 | 12.5 | 0.25 | 20 |
| NMDA | 3.5 | 0.4 | 1 | 1 | 4 | 0.01 | 2.5 |
| eNMDA | 3.5 | 0.4 | 50+200w | 0.2 | 4 | 0.01 | 2.5 |
Parameters for the synapse models with glutamate spillover.
In the cases without spillover the NMDA conductance in the spines is the same as in the cases with spillover, i.e., NMDA gmax = 3.5 nS, and eNMDA gmax = 0 nS. The synaptic weight parameter is w (described in the text).
We have used a short Tdur for NMDARs (the same as for AMPARs), and long Tdur for eNMDARs. Lastly, we have modified the α and β parameters for AMPARs to more realistic values. Also note that we have denoted transmitter duration and transmitter concentration with Tdur and Tmax, respectively, which are given with Cdur and Cmax in the simulation code in Gao et al. ().
eNMDARs have been found to have both diffuse and punctual distribution along dendritic shafts, and electrophysiological measurements have reported their conductance to be 20–60% of the total NMDAR conductance (Harris and Pettit, ; Petralia, ; Papouin and Oliet, ). Even though synapses are represented as punctual inputs in NEURON, the distribution of eNMDAR synapses along 20–30 μm of a dendrite resulted in us effectively modeling a diffuse distribution of eNMDARs. The maximal conductance of an extrasynaptic NMDA synapse in the model is the same as the maximal conductance for a synaptic NMDA synapse. The total synaptic surface area in the clustered spines is around 10% of the surface area of the dendritic shaft where a cluster is situated, making the density of eNMDARs 10 times smaller than the density of NMDARs on the spine.
2.4. Evoking plateau potentials
Plateau potentials are evoked by clustered synaptic input. Clusters ranged from 1 to 40 synapses for the model without spillover and from 1 to 20 synapses for the models with spillover; in SPNs, which are spherically symmetric, a cluster is placed on one of the dendrites, ~120–140 μm from the soma. The difference in the size of the clusters for the cases without and with glutamate spillover is in order to compare the same total NMDAR conductance in both cases. Since both cases have equal spine NMDAR conductance and the case without spillover has no eNMDARs, this requires that it possesses either a twice larger cluster size or a twice larger spine NMDAR conductance to match the total NMDAR conductance in the spillover scenarios; we opted for the former. Clustered synaptic input consists of one presynaptic spike per synapse, all arriving randomly within a 30 ms interval. When no glutamate spillover is modeled, synaptic inputs arrive only to the spines. For glutamate spillover, the same presynaptic spike is delivered to the extrasynaptic NMDA conductance placed in the dendritic shaft under that spine. The study by Gao et al. () models clusters ranging from 10 synapses upwards, which leaves the question of how to treat glutamate spillover for clusters with fewer than 10 synapses open. We have explored two scenarios, which we call thresholded spillover and accumulative spillover. In thresholded spillover we assume glutamate spills over suddenly and simultaneously to all eNMDARs once a glutamate threshold is reached, with a delay of 1.5 ms. The delay is due to experimental data indicating spillover currents arise after 1–2 ms of synaptic currents (Szapiro and Barbour, 2007). In accumulative glutamate spillover there is no threshold for glutamate, and a presynaptic spike delivered to a synaptic NMDA synapse is also delivered to the corresponding extrasynaptic NMDA conductance after a delay of 1.5 ms. This means that every synaptic activation causes spillover and activates the corresponding eNMDARs. An additional difference between the model in Gao et al. () and our study is that in Gao et al. (), larger synaptic clusters also have synapses with bigger weights (parameter w) and longer transmitter pulses (parameter Tdur), whereas in our study we have kept these parameters fixed when increasing cluster size.
The glutamate threshold in thresholded spillover is implemented in the following way. Each NMDA synapse has a weight parameter (w), which is a dimensionless number between 0 and 1 that scales the maximal conductance of the NMDAR. We have set up the model so that a cluster of at least ten synapses with a weight >0.4 needs to be activated to reach the threshold level of glutamate at which spillover occurs, and a spike is delivered to the eNMDARs. We implemented this in NEURON by adding the NMDA synaptic weight to a variable that conceptually models the concentration of glutamate each time a spike arrives at a spine. Added weights are normalized so that ten activated synapses with a weight >0.4 reach the glutamate threshold. To implement this variable, we have in fact used a NEURON integrate-and-fire cell that generates a spike to the eNMDARs in the cluster once it reaches a threshold value of 1.
The two glutamate spillover scenarios we consider are idealized scenarios, and biological glutamate spillover is most likely somewhere in between, varying among brain regions.
2.5. Calculation of phase plots for a single membrane compartment
To investigate the role of local membrane conductances in generating plateau potentials, we performed phase plot analysis for a single electrical compartment which we aimed to make similar to a dendritic region where plateau potentials are evoked. For this purpose, in the single compartment, we inserted the ionic channels Naf, Kaf, Kas, Kir, Kdr with maximal conductances taken from a dendritic segment in the model SPN at ~120 μm from the soma, where clustered inputs are located in most of the results showing voltage traces. (In NEURON nomenclature, the dendritic segment is located at position 0.6 of dendrite 3.) Calcium channels were omitted to simplify the analysis since they have been found to have a minor effect on the voltage, and SK and BK channels have also been omitted since they depend on intracellular calcium concentration, which we do not model in the single compartment. The equation describing the current in the single membrane compartment is:
or
where
In the equations above is the capacitive current, Eion, ENMDA, and Eleak are the reversal potentials and gion, , and gleak conductances of the ionic, NMDAR, and leak currents, respectively, and g(V) is the gating function describing the Mg2+ block. The construction of the phase plots consists of calculating the stationary membrane current as dependent on the voltage (as if a voltage value has been “on” for a long time), which means that the gating variables in the ion channels have their asymptotic, time-independent values. These are denoted as . For example, for the Naf channel, this variable is . Since, in addition to the NMDA current, bistable membrane behavior requires a shunting current of longer duration, such as Kir or GABAB currents (comparable to the duration of NMDA currents), and since the available K currents in this dendritic location have small conductances relative to the NMDA conductances, we have added a leak conductance in the same order of magnitude as the NMDA current to be able to obtain phase curves in the bistable regime (Lazarewicz et al., ; Shoemaker, 2011; Sanders et al., 2013). As described in the Results section, the y-axes of the phase plot show , i.e., they show the current:
2.6. Procedure for fitting the gating function parameters from experimental data
We follow the procedure presented in Jahr and Stevens () to fit the parameters α and η we list in Table 2. We first digitize the data using WebPlotDigitizer (Rohatgi, 2022) for current-voltage relationships for the available Mg2+ concentration from the figures in the articles, from which we then derive the conductance-voltage relationships. When converting the current amplitudes to conductance, the different voltage reversal potential across different experimental conditions was taken into account. Specifically, following Ohm's Law, . Conductance-voltage relationships were normalized with the maximal conductance of each curve, which is the point obtained for the highest applied voltage (usually +50 mV or above) when all of the NMDARs are free from Mg2+. Finally, the normalized conductance-voltage relationships are fitted to Equation (4) using the Python function scipy.optimize.curve_fit.
3. Results
3.1. Thresholded glutamate spillover causes robust all-or-none behavior of plateau potentials
In Figure 2, we compare the plateau potentials generated by a cluster of increasing size with the two spillover models and without glutamate spillover and with the originally reported parameters for the Mg2+ block but corrected for the value of the junction potential (i.e., the parameters found in Jahr and Stevens, but corrected as in Ecker et al., ). The somatic voltage traces in Figures 2B1, B2 show that plateau potentials generated without glutamate spillover exhibit a graded increase in amplitude, whereas those generated with thresholded glutamate spillover exhibit an all-or-none response. As described in the Methods, in order to compare the same total NMDAR conductance between the two scenarios, the model with no spillover has a twice larger cluster size than the models with glutamate spillover.
Figure 2
We also varied the steepness of the gating function by varying the parameter α in Equation (4), and the amplitude of the plateau potentials thus generated is shown as heatmaps in Figures 2C1, C2. The all-or-none behavior can be seen as a sharp jump in the colors of the heatmaps for a small increase in cluster size (i.e., stimulus strength for evoking a plateau). Without glutamate spillover the all-or-none behavior is present only for steep sigmoidal curves, whereas with thresholded glutamate spillover, it is always present. (Steeper gating functions require more excitation to overcome the Mg2+ block in both cases). As described above, we also estimated the gating function parameters for available experimental data (Table 2, sigmoid curves plotted in Supplementary Figure 2). Most of these gating functions are more shallow than the originally reported one by Jahr and Stevens (), and result in graded NMDA potentials in the model without spillover, but still provide all-or-none plateau potentials in the model with glutamate spillover. When varying the parameter η, which shifts the gating function along the x-axis, the same results regarding the all-or-none behavior are observed for the models with thresholded spillover and without glutamate spillover as when varying the parameter α (Figures 2D1, D2). Shifting the gating functions to very depolarized values disables the generation of plateau potentials, as seen from the heatmaps in Figures 2D1, D2. In Figure 2E, we compare the maximal difference between the somatic voltage amplitudes of two consecutive voltage traces, termed the “size of the jump,” δ (indicated in Figure 2B4). In our view, the size of the jump, δ, is essentially proportional to the system's nonlinearity. For example, a large sudden jump (≈10 mV) indicates all-or-none behavior, whereas small jumps (of a few mV) would indicate graded potentials. In our model, thresholded glutamate spillover produced voltage jumps of 8–12 mV, while no spillover produced jumps of ≈3 mV, typically (Figure 2B1).
In addition, the somatic amplitude of plateau potentials generated with thresholded spillover and without glutamate spillover decays with increasing the cluster distance from the soma as shown in Supplementary Figure 3 (Major et al., ).
3.2. Importance of extrasynaptic NMDAR activation for all-or-none plateau potentials: accumulative vs thresholded spillover
Physiologically relevant parallel fiber stimulation protocols cause gradual accumulation of glutamate with each stimulation pulse in cerebellar extrasynaptic space (Okubo et al., ). To explore the effects of such more gradual spillover, we performed simulations where each activation of spine NMDARs also activates the corresponding eNMDARs in the dendritic shaft under the spine after a delay, without the need to first reach a glutamate threshold. (the accumulative spillover described in Section 2).
The results of these simulations show that accumulative glutamate spillover exhibits plateau potential behavior which is between the behavior without spillover and with thresholded glutamate spillover (Figures 2B3, C3, D3). Shallow gating functions provide graded plateau potentials, but all-or-none behavior arises much sooner for increasing steepness than the in case without spillover. With thresholded spillover, clusters of up to 10 synapses only access the spine NMDAR conductance, and clusters with more than 10 synapses activate both NMDAR and eNMDAR conductances, thus producing the large voltage jump before and after spillover (Figure 2B2). With accumulative spillover, eNMDARs are activated with a short delay after each NMDAR activation, giving rise to the more graded increase in plateau potential amplitude, as well as duration (Figure 2B3). Moreover, for a larger proportion of eNMDAR conductance in the total conductance, the behavior is more all-or-none (Figures 2B4, E). This suggests that there may be particular conditions under which all-or-none plateau potentials appear and that plateau potential behavior might vary among brain regions.
3.3. Spillover induces all-or-none plateau potentials by facilitating the switch to bistable or self-triggering regime
To explain the results in the previous section, we use phase plot analysis of the membrane voltage calculated for a single membrane compartment (Figure 3). Phase plot analysis reveals the role of different membrane conductances in generating nonlinear behavior such as plateau potentials. In particular, it has been found that apart from the nonlinear, voltage-dependent excitatory NMDA current, an opposing, shunting current is also required acting at least on the longer timescale of NMDA currents. This role can be played by Kir channels, GABAB currents, a passive leak current, and/or an axial current diverting charge to the rest of the neuron (Lazarewicz et al., ; Shoemaker, 2011; Sanders et al., 2013).
Figure 3
We begin with a short description of how to read the phase plots (Figure 3A). A voltage phase plot normally contains voltage on the x-axis and the rate of voltage change (the time derivative of the voltage, ) on the y-axis. We have instead plotted the membrane current on the y-axis, which is simply a linear rescaling of with the membrane capacitance, . The intersections of the phase curve with the x-axis correspond to fixed points in the voltage dynamics (). Since our single compartment is a one-dimensional system (voltage is the only independent state variable), a point in the phase plot can move only along the phase curve. Initial voltage points that correspond to positive values of the phase curve () will get increased in time, whereas those corresponding to negative values get decreased. Voltage change stops when a point on the phase plot reaches a stable fixed point (where ). Points near unstable fixed points are repelled away from them. The trajectory from the multicompartment SPN model shown in Figure 3D1, on the other hand, can take on any values in the phase plot. This is because in the recorded dendritic compartment, the gating variables of the ion channels are also time-varying, i.e., they are additional state variables in the model. We note that the phase plot is different from the commonly used and similarly named I–V plot, which can be found, for example, in the aforementioned articles measuring macroscopic NMDA currents.
Depending on the total NMDAR conductance in a dendritic region, NMDARs can exhibit three types of nonlinear dynamics, as shown in Figure 3A: boosting, bistable, and self-triggering, of which the latter two correspond to plateau potentials (Schiller and Schiller, 2001). Boosting behavior occurs for small or moderate amounts of total NMDAR conductance which is not enough to trigger plateau potentials. In this case other excitatory inputs are boosted in amplitude by activating the longer-lasting NMDA potentials. The phase curve has one stable fixed point to the left of the diagram at a low voltage. Bistable behavior is exhibited when the phase curve has an inverted-N shape with three interceptions with the x-axis (fixed points), which happens for larger amounts of total NMDAR conductance. In this case a plateau potential is generated for an excitation that exceeds the voltage threshold determined by the middle (unstable) fixed point. Such an excitation causes a fast excursion in the phase plot to the right (higher) stable fixed point, which causes the all-or-none behavior (region 2 in the phase plot in Figure 3D1, which corresponds to the upward swing in the corresponding dendritic voltage in Figure 3D2, and animation in Supplementary data). Even larger total NMDAR conductance gives rise to the so-called self-triggering behavior. In this case the phase curve, now situated above the bistable region, is inverted-N-shaped and has one stable fixed point to the right of the diagram which corresponds to high voltage. In this case, any excitation produces a plateau potential by carrying the voltage to the single fixed point. Thus, as the total NMDAR conductance increases with successive activations of synapses in a cluster, the phase curve sweeps potentially across all three regions in the phase plot. On the other hand, closing of the NMDARs moves the phase curve back in the other direction, collapsing the inverted-N shape back toward a straight line and terminating the plateau potential with the voltage dropping to the low-voltage, single stable fixed point, which corresponds to the resting voltage.
To elicit plateau potentials, the total NMDAR conductance must be enough to situate the phase curve in the bistable or self-triggering region (the threshold level of NMDAR conductance being at the border between the boosting and the bistable region). An example of a plateau potential generated with thresholded glutamate spillover, with the corresponding phase plot trajectory, is given in Figures 3D1, D2. The gray phase curve in Figure 3D1 is the one obtained when only synaptic NMDAR conductance is activated, whereas the purple curve represents the total NMDAR conductance with spillover, when eNMDARs are activated. After the subthreshold excitation by nine synaptic inputs (region 1), the 10th input causes spillover and the dendritic voltage crosses the threshold determined by the unstable fixed point. This initiates the self-sustaining suprathreshold excitation (region 2), which is followed by the plateau region (region 3). The end of the extrasynaptic transmitter pulse starts the downward, hyper-polarizing phase (region 4).
To examine the role of the gating function parameters α and η, we constructed phase plots with varying α and η. Figures 3B1, B2 are a model of the total NMDAR conductance when both NMDARs an eNMDARs are activated during spillover, as well as the case when all of this total conductance is situated in the spines of the equivalent model without spillover (Figure 2B). Varying the steepness of the Mg2+ block shifts the phase curve vertically, and all three phase plot regimes may be obtained (Figure 3B1). Shallow sigmoidal curves have only one fixed point and are in the self-triggering region, whereas steep sigmoidal curves have three fixed points and are in the bistable region. Very steep sigmoidal curves are in the boosting region. In Figures 3C1, C2, on the other hand, the total NMDAR conductance is halved, to model only the synaptic NMDAR conductance accessible before thresholded spillover. In this case, all phase curves for the single membrane compartment are in the boosting region, regardless of the gating function steepness.
These results from the phase plots can explain the plateau potential behavior with and without spillover. Total NMDAR conductance increases with increasing cluster size—larger clusters have greater total NMDA conductance (Figure 2B, cluster size). When there is no glutamate spillover, all-or-none behavior is observed only for steep sigmoidal curves, since the corresponding phase curves are in the bistable regime, while for the shallow sigmoidal curves, all synaptic clusters of increasing size produce graded excitations as in Figure 2B1, because the phase curve is in the self-triggering region with a single, high-voltage fixed point (Figure 3B1). Graded excitations trigger NMDA spikes first and plateau potentials later, as the total NMDA conductance increases. On the other hand, in the model with glutamate spillover, activating only the synaptic NMDA conductance (extrasynaptic NMDARs are inactive) is insufficient to elicit plateau potentials (the phase curve is most likely in the boosting region as in Figure 3C1). When spillover occurs, the activation of the eNMDARs provides enough NMDA conductance to elicit a plateau potential. This shifts the phase curve from the boosting to the self-triggering region for shallow gating functions and to the bistable region for steep gating functions. The duration of the phase curve in both of these regions is significantly prolonged due to glutamate spillover. This shift (and inversion) of the phase curve occurs suddenly because, as described in the Methods section, in thresholded glutamate spillover all the eNMDAR conductance is simultaneously activated. For shallow sigmoidal curves, it is this sudden inversion that causes the all-or-none plateau potential behavior. For steep sigmoidal curves, the voltage excitation also needs to be higher than the voltage threshold for a plateau to occur. In this way, glutamate spillover causes all-or-none plateau potentials for a wide range of gating function steepness. Accumulative spillover, on the other hand, is very similar to no spillover. Each successive synaptic activation gradually increases the total NMDAR conductance, slowly moving the phase curve across the phase plot.
Varying the parameter η shifts the gating function along the x-axis. Shifting the function to higher potentials flattens (linearizes) the phase curve. While all three regimes can be obtained with various values of η, shifting the gating function to higher potentials results in only boosting behavior and disappearance of plateau potentials both with and without glutamate spillover (Figures 2D1–D3, 3C1, C2).
3.4. Long glutamate duration drives all-or-none plateau potentials
We verify that it is indeed the length of glutamate spillover that causes all-or-none plateau potential behavior by manipulating the duration of the transmitter pulse (the parameter Tdur, Figure 4). As explained in the Methods section, Tdur for the synapses onto spines is much shorter than Tdur for eNMDARs. Increasing Tdur in the case without spillover results in all-or-none plateau potential behavior. Again, increasing NMDAR conductance by increasing cluster size moves the stable fixed point further to the right in the phase plot in Figure 3A. When entering the bistable regime, the all-or-none jump appears in Figure 4A1, and lasts because the increased Tdur keeps the phase curve in the bistable regime for a longer time. On the other hand, decreasing the eNMDAR Tdur in the case with spillover produces rather graded plateau potentials with much smaller voltage jumps (Figure 4B1). The all-or-none behavior in Figure 4A1 and graded behavior in Figure 4B1 are also evident in the sharp and gradual color transitions in the respective heatmaps in Figures 4A2, B2, as well as in the comparison of the voltage jump, δ in Figure 4C.
Figure 4
3.5. Low extrasynaptic glutamate concentration can produce all-or-none plateau potentials
The extrasynaptic glutamate concentration (the parameter Tmax for eNMDARs) used in the simulations above is 200 μM, an order of magnitude higher than the measurements reported in Okubo et al. (). NMDARs have high-affinity for glutamate, in the low-micromolar range, so they should be able to respond to low glutamate concentrations. We have tested the thresholded spillover model with lower values for extracellular glutamate concentration (the parameter Tmax) and the results are shown in Figure 5. Concentrations of 10 and 20 μM, significantly higher than the NMDAR dissociation constant, Kd = 2.5 μM, produce all-or-none plateau potentials (Figures 5A, B). For concentrations of 5 μM, the size of the voltage jump is reduced, although the plateau potentials retain their all-or-none quality. Concentrations comparable to the NMDAR Kd do not activate the eNMDARs fully and produce graded plateau potential amplitude (Figures 5D, E).
Figure 5
3.6. All-or-none behavior of plateau potentials enables nonlinear computation, such as solving the nonlinear feature binding problem
Linearly non-separable tasks, an example of which is the nonlinear feature binding problem (NFBP) cannot be solved by the perceptron, a computational unit with a single nonlinearity, and require a network of such artificial neurons. As such, these tasks are traditionally used as a benchmark for the computational capabilities of a computational unit. Possessing dendritic nonlinearities in addition to the somatic nonlinearity, single biological neurons have been shown to be equivalent to at least a two-layer artificial neural network (Poirazi et al., ) and should, in principle, be able to solve such tasks (Tran-Van-Minh et al., 2015). Although it is not known whether the NFBP is solved by single neurons in vivo, and whether it is a relevant task for any brain region, human cortical layer 2/3 pyramidal neurons have been found to possess calcium-mediated dendritic spikes whose voltage dependency enables the dendrite exhibiting them to solve the exclusive OR (XOR) problem in simulations (Gidon et al., ), and the XOR problem is related to the NFBP (Cazé et al., ). These dendritic nonlinearities occur in the apical dendrites, while the SPN dendrites correspond to basal dendrites in pyramidal neurons, where NMDA nonlinearities occur. We, therefore, use the NFBP as a benchmark to demonstrate how all-or-none plateau potentials enable a single neuron to solve this task.
We illustrate the NFBP with an example from visual feature binding (Figure 6A). A stimulus has two features (shape and color), each of which can have two values (strawberry or banana for the shape, and red or yellow for the color), for a total of four possible combinations. Each stimulus excites the neuron with the same amount of excitation, on average. The task consists of responding to two of the feature combinations (the relevant stimuli) with a somatic spike and remaining silent for the other two (irrelevant stimuli). This means that the same amount of total excitation should be processed differently, and the plateau potentials offer a natural solution to the task: the synapses for the relevant stimuli should be clustered in separate dendrites (Figure 6B), and should elicit a plateau potential upon arrival of a relevant stimulus which will drive somatic spiking. Irrelevant stimuli would activate only half of the cluster and should not elicit a plateau potential (Tran-Van-Minh et al., 2015). This setup with thresholded glutamate spillover provides perfect performance on the NFBP (Figures 6C, D). Without spillover, the graded amplitudes of the plateaus cause somatic spiking for the irrelevant stimuli as well, and the neuron cannot perform this task (Figures 6C, E, F).
Figure 6
4. Discussion
In this computational study we investigated how to generate robust all-or-none dendritic plateau potentials using striatal projection neuron models as a test case. We also illustrated how such plateau potentials can enhance computation. Dendritic nonlinearities in the form of NMDA spikes and plateau potentials offer an enhanced ability for dendritic computation compared to the linear summation of excitatory inputs (Poirazi et al., ; Lavzin et al., ; Xu et al., 2012; Tran-Van-Minh et al., 2015; Kumar et al., ). By summating co-activated clustered inputs supralinearly, they can allow for the preferential detection of certain input patterns or combinations of inputs. For example, in the mouse barrel cortex, layer 4 spiny stellate neurons generate local and global multi-branch NMDA spikes that contribute substantially to the angular tuning of these neurons (Lavzin et al., ). In addition, distal dendrites of layer 5 pyramidal neurons in the barrel cortex might integrate correlated sensory and motor information via dendritic nonlinearities in order to produce a signal related to object localization during active sensing tasks (Xu et al., 2012). Dendritic plateau potentials are also essential processes for triggering synaptic plasticity both in vitro and in vivo (Gambino et al., ; Cichon and Gan, ; Brandalise et al., ). For example, in the somatosensory cortex, rhythmic sensory whisker stimulation can induce synaptic long-term potentiation (LTP) in layer 2/3 pyramidal cells, which is triggered by plateau potentials generated through the cooperative activity of the intracolumnar lemniscal and thalamocortical paralemniscal synaptic circuitry in the absence of somatic spiking (Gambino et al., ). Additionally, different motor learning tasks induce dendritic calcium spikes on different apical tuft branches of layer 5 pyramidal neurons in the mouse motor cortex (Cichon and Gan, ). These task-related, branch-specific calcium spikes cause long-lasting potentiation of postsynaptic dendritic spines active at the time of spike generation, suggesting a role for dendritic nonlinearities in storing new information without disrupting previously acquired memories (Cichon and Gan, ). Lastly, dendritic nonlinearities have been suggested to enable a neuron to solve linearly non-separable tasks such as the NFBP (Tran-Van-Minh et al., 2015).
Experimentally evoked plateau potentials in vitro display all-or-none behavior, and this feature is not easily nor robustly captured when modeling. A critical ingredient for the all-or-none property is the sigmoidal shape of the Mg2+ block voltage dependence, described by the gating function in Equation (4). The many NMDAR isoforms show different sensitivity to Mg2+, which results in a range of possible shapes of the sigmoid gating function in terms of its steepness and the voltage required to overcome the Mg2+ block. We have studied how the variability of these gating function properties affects the all-or-none plateau potential behavior. We found that the all-or-none property is very sensitive to the steepness of the gating function. Shallow gating functions are in the self-triggering region of the phase plot and result in graded plateau potential amplitudes for increasing cluster size, whereas steep gating functions are in the bistable region and result in all-or-none plateau potentials (Figures 2B1, B2, 3B1). Importantly, including glutamate spillover in the models provides a robust all-or-none quality of the plateau potentials for the whole tested range of Mg2+ block steepness.
Almost all the gating functions we fitted to available experimental data are more shallow than the originally reported sigmoidal curve by Jahr and Stevens (), which already results in graded plateau potential amplitude in the SPN model (Figure 2B; Table 2). Conversely, all of the studies that model plateau potentials, except for Ecker et al. (), that correct for the omitted junction potential in Jahr and Stevens (), have increased the steepness of the originally reported gating function without explicitly stating a reason for it. As shown in Figure 3B1, increasing the steepness of the gating function can move the phase curve in the bistable region, allowing to robustly generate all-or-none plateau potentials, thus offering an explanation for the steeper gating functions in modeling studies. Moreover, in light of the natural variability in the Mg2+ block gating function, the question remains as to how SPNs generate the all-or-none plateau behavior robustly in experiments with increasing stimulation. According to the results above, glutamate spillover is predicted to be a part of the plateau potential generation mechanism, providing robust all-or-none behavior across a wide range of slopes of the Mg2+ block gating function.
Whether, in vivo plateau potentials are graded or all-or-none is not known. However, in Gambino et al. () for example, robust, seemingly all-or-none plateau potentials are evoked with co-activation of the intracolumnar lemniscal and thalamocortical paralemniscal pathways, while lemniscal pathway activation alone significantly reduces the probability of evoking plateau potentials. In this study we have illustrated the importance of such all-or-none plateau potentials in a commonly used benchmark task, the NFBP, showing that graded plateau potentials cannot perform this task compared to the perfect performance with all-or-none plateau potentials. It is not known whether the NFBP is solved by single neurons in any brain region on a regular basis, but plateau potentials, whether used for the NFBP or in the detection/representation of single features, have been suggested to prepare the neuron to spike in a robust fashion as a response to the particular plateau-evoking input (Antic et al., ). In vitro, distally evoked plateau potentials in the basal dendrites of pyramidal neurons typically do not evoke somatic spiking by themselves (Milojkovic et al., ; Major et al., ; Gao et al., ), and neither do plateau potentials in SPNs (Plotkin et al., ; Du et al., ). This corresponds to the situation we modeled in Figures 2, 4, 5. However, in the presence of additional inputs or background noise, plateau potentials drive a robust somatic response to the particular input pattern that evoked them, as we exemplify in Figure 6. According to the predictions from this study, the prolonged activation of extrasynaptic NMDARs is critical for robust plateau potential generation.
Despite important advances such as the confirmation of glutamate spillover in vivo, the confirmation of its role in neuronal signaling in the cerebellum, and reports of extracellular glutamate concentrations in experimentally evoked spillover in brain slices in the low micromolar ranges, many important details await further studies (Szapiro and Barbour, 2007; Okubo et al., ; Papouin and Oliet, ; Rusakov and Stewart, 2021). Some prominent missing details are: (i) the lifetime of glutamate in the extrasynaptic and perisynaptic spaces; (ii) the conditions under which possible astrocytic glutamate release and/or reversal of glutamate uptake occurs; as well as (iii) control of eNMDAR activation by the co-agonist glycine or d-serine (Malarkey and Parpura, ; Papouin and Oliet, ). Assuming high NMDAR affinity for glutamate, we have shown that 10–20 μM of extracellular glutamate can produce all-or-none plateau potentials if glutamate persists in the extracellular space for long durations, indicating that physiological concentrations may play a large role in signal integration. Recent findings that LTP boosts glutamate spillover by initiating withdrawal of perisynaptic astroglial processes suggest dynamic regulation of the synaptic and perisynaptic environment, including eNMDAR activation, for the purposes of signal processing and learning (Henneberger et al., ; Rusakov and Stewart, 2021). Facilitation of the prolonged eNMDAR activation would, in turn, increase the probability that the already strengthened synapses can evoke a robust dendritic plateau potential.
One direction for future work would be to investigate how the distinct gating functions of various NMDAR isoforms affect dendritic computations in different neuron types in the brain, as well as to further quantify the contribution of extrasynaptic NMDA receptors for these computations (Zhou et al., 2015). Another direction for future work would be to develop more detailed diffusion models of glutamate concentration activating extrasynaptic NMDARs for a quantitatively more precise understanding of the phenomenon in particular synapse types. However, since NMDARs have high affinity for glutamate, even low extrasynaptic glutamate concentrations also generate robust plateau potentials for a prolonged extrasynaptic glutamate signal (as shown in Figure 5). The synaptic model that we used assumes fixed transmitter concentration for the duration of the transmitter pulse both in the synaptic cleft and in the extrasynaptic space. Glutamate in the synaptic cleft rapidly reaches its peak concentration and is also quickly cleared due to diffusion and reuptake (Clements et al., ). In contrast to that, extrasynaptic glutamate dynamics is likely more variable. Under conditions of temporary saturation of astrocytic glutamate uptake and possibly Ca2+-dependent astrocytic glutamate exocytosis and/or reverse operation of astrocytic glutamate transporters (Malarkey and Parpura, ), extrasynaptic glutamate concentration could remain elevated long enough to produce all-or-none plateau potentials. Further experimental studies of extrasynaptic glutamate are needed to clarify these issues in different brain regions.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://github.com/danieltrpevski/Spillover/tree/master and on ModelDB (McDougal et al., ).
Author contributions
DT designed the study, performed the simulations, and wrote the first draft of the manuscript. ZK contributed to the study's design and was involved in the simulations. IC performed the parameter estimation from experimental articles and wrote sections of the manuscript. JHK supervised the project. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Swedish Research Council (VR-M-2020-01652), Swedish e-Science Research Centre (SeRC), EU/Horizon 2020 no. 945539 (HBP SGA3), and KTH Digital Futures.
Acknowledgments
The simulations were enabled by resources provided by the Swedish National Infrastructure for Computing (SNIC) at PDC KTH partially funded by the Swedish Research Council through grant agreement no. 2018-05973. We acknowledge the use of Fenix Infrastructure resources, which are partially funded from the European Union's Horizon 2020 research and innovation programme through the ICEI project under the grant agreement No. 800858. We thank Dr. Robert Lindroos for helpful discussion and useful feedback.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2023.1196182/full#supplementary-material
Video S1Video animation of the plots in Figures 3D1, 3D2. The trajectory traversing the phase plot in the upper panel corresponds to the voltage timecourse in the lower panel. For details see figure legend of Figures 3D1, 3D2.
References
1
AnticS. D.HinesM.LyttonW. W. (2018). Embedded ensemble encoding hypothesis: the role of the “prepared” cell. J. Neurosci. Res. 96, 1543–1559. 10.1002/jnr.24240
2
AnticS. D.ZhouW.-L.MooreA. R.ShortS. M.IkonomuK. D. (2010). The decade of the dendritic NMDA spike. J. Neurosci. Res. 88, 2991–3001. 10.1002/jnr.22444
3
BrandaliseF.CartaS.HelmchenF.LismanJ.GerberU. (2016). Dendritic NMDA spikes are necessary for timing-dependent associative LTP in CA3 pyramidal cells. Nat. Commun. 7:13480. 10.1038/ncomms13480
4
CarnevaleN. T.HinesM. L. The NEURON Book. New York, NY: Cambridge University Press (2005).
5
CarterA. G.SabatiniB. L. (2004). State-dependent calcium signaling in dendritic spines of striatal medium spiny neurons. Neuron44, 483–493. 10.1016/j.neuron.2004.10.013
6
CazéR. D.HumphriesM.GutkinB. (2013). Passive dendrites enable single neurons to compute linearly non-separable functions. PLoS Comput. Biol. 9:e1002867. 10.1371/journal.pcbi.1002867
7
ChalifouxJ. R.CarterA. G. (2011). Glutamate spillover promotes the generation of NMDA spikes. J. Neurosci. 31, 16435–16446. 10.1523/JNEUROSCI.2777-11.2011
8
ChenL.HuangL.-Y. M. (1992). Protein kinase C reduces Mg2+ block of NMDA-receptor channels as a mechanism of modulation. Nature356, 521–523. 10.1038/356521a0
9
ChengH.-W.RafolsJ. A.GoshgarianH. G.AnaviY.TongJ.McNeillT. H. (1997). Differential spine loss and regrowth of striatal neurons following multiple forms of deafferentation: a Golgi study. Exp. Neurol. 147, 287–298. 10.1006/exnr.1997.6618
10
ChiuD. N.CarterB. C. (2022). Synaptic NMDA receptor activity at resting membrane potentials. Front. Cell. Neurosci. 16:916626. 10.3389/fncel.2022.916626
11
CichonJ.GanW.-B. (2015). Branch-specific dendritic Ca2+ spikes cause persistent synaptic plasticity. Nature520, 180–185. 10.1038/nature14251
12
ClarkB. A.FarrantM.Cull-CandyS. G. (1997). A direct comparison of the single-channel properties of synaptic and extrasynaptic NMDA receptors. J. Neurosci. 17, 107–116. 10.1523/JNEUROSCI.17-01-00107.1997
13
ClementsJ. D.LesterR. A. J.TongG.JahrC. E.WestbrookG. L. (1992). The time course of glutamate in the synaptic cleft. Science258, 1498–1501. 10.1126/science.1359647
14
DestexheA.MainenZ. F.SejnowskiT. J. (1994). An efficient method for computing synaptic conductances based on a kinetic model of receptor binding. Neural Comput. 6, 14–18. 10.1162/neco.1994.6.1.14
15
DormanD. B.Jkedrzejewska-SzmekJ.BlackwellK. T. (2018). Inhibition enhances spatially-specific calcium encoding of synaptic input patterns in a biologically constrained model. eLife7. 10.7554/eLife.38588.022
16
DoronM.ChindemiG.MullerE.MarkramH.SegevI. (2017). Timed synaptic inhibition shapes NMDA spikes influencing local dendritic processing and global I/O properties of cortical neurons. Cell Rep. 21, 1550–1561. 10.1016/j.celrep.2017.10.035
17
DuK.WuY.-W.LindroosR.LiuY.RózsaB.KatonaG.et al. (2017). Cell-type–specific inhibition of the dendritic plateau potential in striatal spiny projection neurons. Proc. Natl. Acad. Sci. U.S.A. 114, E7612—E7621. 10.1073/pnas.1704893114
18
EckerA.RomaniA.SárayS.KáliS.MiglioreM.FalckJ.et al. (2020). Data-driven integration of hippocampal ca1 synaptic physiologyin silico. Hippocampus30, 1129–1145. 10.1002/hipo.23220
19
FarinellaM.RuedtD. T.GleesonP.LanoreF.SilverR. A. (2014). Glutamate-bound NMDARs arising from in vivo-like network activity extend spatio-temporal integration in a L5 cortical pyramidal cell model. PLoS Comput. Biol. 10:e1003590. 10.1371/journal.pcbi.1003590
20
GambinoF.PagèsS.KehayasV.BaptistaD.TattiR.CarletonA.et al. (2014). Sensory-evoked LTP driven by dendritic plateau potentials in vivo. Nature515, 116–119. 10.1038/nature13664
21
GaoP. P.GrahamJ. W.ZhouW.-L.JangJ.AnguloS.Dura-BernalS.et al. (2021). Local glutamate-mediated dendritic plateau potentials change the state of the cortical pyramidal neuron. J. Neurophysiol. 125, 23–42. 10.1152/jn.00734.2019
22
Garcia-MunozM.Lopez-HuertaV. G.Carrillo-ReidL.ArbuthnottG. W. (2015). Extrasynaptic glutamate NMDA receptors: key players in striatal function. Neuropharmacology89, 54–63. 10.1016/j.neuropharm.2014.09.013
23
GidonA.ZolnikT. A.FidzinskiP.BolduanF.PapoutsiA.PoiraziP.et al. (2020). Dendritic action potentials and computation in human layer 2/3 cortical neurons. Science367, 83–87. 10.1126/science.aax6239
24
HarrisA. Z.PettitD. L. (2007). Extrasynaptic and synaptic NMDA receptors form stable and uniform pools in rat hippocampal slices. J. Physiol. 584, 509–519. 10.1113/jphysiol.2007.137679
25
HennebergerC.BardL.PanatierA.ReynoldsJ. P.KopachO.MedvedevN. I.et al. (2020). LTP induction boosts glutamate spillover by driving withdrawal of perisynaptic astroglia. Neuron108, 919–936.e11. 10.1016/j.neuron.2020.08.030
26
HigleyM. J.SabatiniB. L. (2010). Competitive regulation of synaptic Ca2+influx by D2 dopamine and A2A adenosine receptors. Nat. Neurosci. 13, 958–966. 10.1038/nn.2592
27
JahrC.StevensC. (1990a). A quantitative description of NMDA receptor-channel kinetic behavior. J. Neurosci. 10, 1830–1837. 10.1523/JNEUROSCI.10-06-01830.1990
28
JahrC.StevensC. (1990b). Voltage dependence of NMDA-activated macroscopic conductances predicted by single-channel kinetics. J. Neurosci. 10, 3178–3182. 10.1523/JNEUROSCI.10-09-03178.1990
29
KastellakisG.PoiraziP. (2019). Synaptic clustering and memory formation. Front. Mol. Neurosci. 12:300. 10.3389/fnmol.2019.00300
30
KumarA.SchiffO.BarkaiE.MelB. W.Poleg-PolskyA.SchillerJ. (2018). NMDA spikes mediate amplification of inputs in the rat piriform cortex. eLife7:e38446. 10.7554/eLife.38446.018
31
LavzinM.RapoportS.PolskyA.GarionL.SchillerJ. (2012). Nonlinear dendritic processing determines angular tuning of barrel cortex neurons in vivo. Nature490, 397–401. 10.1038/nature11451
32
LazarewiczM. T.AngC.-W.CarlsonG. C.CoulterD. A.FinkelL. H. (2006). Analysis of nmda-dependent voltage bistability in thin dendritic compartments. Neurocomputing69, 1025–1029. 10.1016/j.neucom.2005.12.038
33
LindroosR.Hellgren KotaleskiJ. (2021). Predicting complex spikes in striatal projection neurons of the direct pathway following neuromodulation by acetylcholine and dopamine. Eur. J. Neurosci. 53, 2117–2134. 10.1111/ejn.14891
34
MajorG.PolskyA.DenkW.SchillerJ.TankD. (2008). Spatiotemporally graded NMDA spike/plateau potentials in basal dendrites of neocortical pyramidal neurons. J. Neurophysiol. 99, 2584–601. 10.1152/jn.00011.2008
35
MalarkeyE. B.ParpuraV. (2008). Mechanisms of glutamate release from astrocytes. Neurochem. Int. 52, 142–154. 10.1016/j.neuint.2007.06.005
36
McDougalR. A.MorseT. M.CarnevaleT.MarencoL.WangR.MiglioreM.et al. (2017). Twenty years of modeldb and beyond: building essential modeling tools for the future of neuroscience. J. Comput. Neurosci. 42, 1–10. 10.1007/s10827-016-0623-7
37
McMenimenK. A.PeterssonE. J.LesterH. A.DoughertyD. A. (2006). Probing the mg2+ blockade site of an n-methyl-d-aspartate (NMDA) receptor with unnatural amino acid mutagenesis. ACS Chem. Biol. 1, 227–234. 10.1021/cb6000944
38
MilojkovicB. A.RadojicicM. S.Goldman-RakicP. S.AnticS. D. (2004). Burst generation in rat pyramidal neurones by regenerative potentials elicited in a restricted part of the basilar dendritic tree. J. Physiol. 558, 193–211. 10.1113/jphysiol.2004.061416
39
NowakL.BregestovskiP.AscherP.HerbetA.ProchiantzA. (1984). Magnesium gates glutamate-activated channels in mouse central neurones. Nature307, 462–465. 10.1038/307462a0
40
OikonomouK. D.ShortS. M.RichM. T.AnticS. D. (2012). Extrasynaptic glutamate receptor activation as cellular bases for dynamic range compression in pyramidal neurons. Front. Physiol. 3:334. 10.3389/fphys.2012.00334
41
OikonomouK. D.SinghM. B.SterjanajE. V.AnticS. D. (2014). Spiny neurons of amygdala striatum, and cortex use dendritic plateau potentials to detect network UP states. Front. Cell. Neurosci. 8:292. 10.3389/fncel.2014.00292
42
OkuboY.SekiyaH.NamikiS.SakamotoH.IinumaS.YamasakiM.et al. (2010). Imaging extrasynaptic glutamate dynamics in the brain. Proc. Natl. Acad. Sci. U.S.A. 107, 6526–6531. 10.1073/pnas.0913154107
43
PapouinT.OlietS. H. R. (2014). Organization, control and function of extrasynaptic nmda receptors. Philos. Trans. R. Soc. B Biol. Sci. 369:20130601. 10.1098/rstb.2013.0601
44
PetraliaR. S. (2012). Distribution of extrasynaptic NMDA receptors on neurons. Sci. World J. 2012:267120. 10.1100/2012/267120
45
PlotkinJ. L.DayM.SurmeierD. J. (2011). Synaptically driven state transitions in distal dendrites of striatal spiny neurons. Nat. Neurosci. 14, 881–888. 10.1038/nn.2848
46
PoiraziP.BrannonT.MelB. W. (2003). Pyramidal neuron as two-layer neural network. Neuron37, 989–999. 10.1016/S0896-6273(03)00149-1
47
Poleg-PolskyA. (2015). Effects of neural morphology and input distribution on synaptic processing by global and focal NMDA-spikes. PLoS ONE10:e0140254. 10.1371/journal.pone.0140254
48
QianA.BullerA. L.JohnsonJ. W. (2005). Nr2 subunit-dependence of NMDA receptor channel block by external mg2+. J. Physiol. 562, 319–331. 10.1113/jphysiol.2004.076737
49
QianA.JohnsonJ. W. (2006). Permeant ion effects on external mg2+ block of nr1/2d NMDA receptors. J. Neurosci. 26, 10899–10910. 10.1523/JNEUROSCI.3453-06.2006
50
RetchlessB. S.GaoW.JohnsonJ. W. (2012). A single GluN2 subunit residue controls NMDA receptor channel properties via intersubunit interaction. Nat. Neurosci. 15, 406–413. 10.1038/nn.3025
51
RhodesP. (2006). The properties and implications of NMDA spikes in neocortical pyramidal cells. J. Neurosci. 26, 6704–6715. 10.1523/JNEUROSCI.3791-05.2006
52
RohatgiA. (2022). Webplotdigitizer: Version 4.6.
53
RusakovD. A.StewartM. G. (2021). Synaptic environment and extrasynaptic glutamate signals: the quest continues. Neuropharmacology195:108688. 10.1016/j.neuropharm.2021.108688
54
SandersH.BerendsM.MajorG.GoldmanM. S.LismanJ. E. (2013). NMDA and gabab (kir) conductances: the “perfect couple” for bistability. J. Neurosci. 33, 424–429. 10.1523/JNEUROSCI.1854-12.2013
55
SchillerJ.SchillerY. (2001). NMDA receptor-mediated dendritic spikes and coincident signal amplification. Curr. Opin. Neurobiol. 11, 343–348. 10.1016/S0959-4388(00)00217-8
56
SharmaG.StevensC. F. (1996). A mutation that alters magnesium block of N-methyl-D-aspartate receptor channels. Proc. Natl. Acad. Sci. U.S.A. 93, 9259–9263. 10.1073/pnas.93.17.9259
57
ShindouT.Ochi-ShindouM.WickensJ. (2011). A Ca(2+) threshold for induction of spike-timing-dependent depression in the mouse striatum. J. Neurosci. 31, 13015–13022. 10.1523/JNEUROSCI.3206-11.2011
58
ShoemakerP. A. (2011). Neural bistability and amplification mediated by NMDA receptors: analysis of stationary equations. Neurocomputing74, 3058–3071. 10.1016/j.neucom.2011.04.018
59
SzapiroG.BarbourB. (2007). Multiple climbing fibers signal to molecular layer interneurons exclusively via glutamate spillover. Nat. Neurosci. 10, 735–742. 10.1038/nn1907
60
Tran-Van-MinhA.CazéR. D.AbrahamssonT.CathalaL.GutkinB. S.DiGregorioD. A. (2015). Contribution of sublinear and supralinear dendritic integration to neuronal computations. Front. Cell. Neurosci. 9:67. 10.3389/fncel.2015.00067
61
XuN.HarnettM. T.WilliamsS. R.HuberD.O'ConnorD. H.SvobodaK.et al. (2012). Nonlinear dendritic integration of sensory and motor input during an active sensing task. Nature492, 247–254. 10.1038/nature11601
62
ZhouX.ChenZ.YunW.RenJ.LiC.WangH. (2015). Revisiting the function of extrasynaptic nmda receptor in excitotoxicity running title: NMDA receptor and excitotoxicity. Neuroscientist21:337. 10.1177/1073858414548724
Summary
Keywords
glutamate spillover, plateau potentials, NMDA spikes, gating function, computational modeling, nonlinear dendritic computation, clustered synapses, magnesium block of NMDA receptors
Citation
Trpevski D, Khodadadi Z, Carannante I and Hellgren Kotaleski J (2023) Glutamate spillover drives robust all-or-none dendritic plateau potentials—an in silico investigation using models of striatal projection neurons. Front. Cell. Neurosci. 17:1196182. doi: 10.3389/fncel.2023.1196182
Received
29 March 2023
Accepted
31 May 2023
Published
29 June 2023
Volume
17 - 2023
Edited by
Srdjan D. Antic, University of Connecticut Health Center, United States
Reviewed by
William W. Lytton, Downstate Health Sciences University, United States; Pierre Apostolides, University of Michigan, United States
Updates
Copyright
© 2023 Trpevski, Khodadadi, Carannante and Hellgren Kotaleski.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniel Trpevski trpevski@kth.seJeanette Hellgren Kotaleski jeanette@kth.se
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.