Original Research ARTICLE
Depotentiation from Potentiated Synaptic Strength in a Tristable System of Coupled Phosphatase and Kinase
- 1School of Physics and Information Technology, Shaanxi Normal University, Xi'an, China
- 2Key Laboratory of Modern Teaching Technology, Ministry of Education, Shaanxi Normal University, Xi'an, China
Long-term potentiation (LTP) of synaptic strength is strongly implicated in learning and memory. On the other hand, depotentiation, the reversal of synaptic strength from potentiated LTP state to the pre-LTP level, is required in extinction of the obsolete memory. A generic tristable system, which couples the phosphatase and kinase switches, exclusively explains how moderate and high elevation of intracellular calcium concentration triggers long-term depression (LTD) and LTP, respectively. The present study, introducing calcium influx and calcium release from internal store into the tristable system, further show that significant elevation of cytoplasmic calcium concentration switches activation of both kinase and phosphatase to their basal states, thereby depotentiate the synaptic strength. A phase-plane analysis of the combined model was employed to explain the previously reported depotentiation in experiments and predict a threshold-like effect with calcium concentration. The results not only reveal a mechanism of NMDAR- and mGluR-dependent depotentiation, but also predict further experiments about the role of internal calcium store in induction of depotentiation and extinction of established memories.
Over the past decades, ample evidences support that long-lasting modifications of synaptic strength, including long-term potentiation (LTP) and long-term depression (LTD), are strongly implicated in learning and memory. Meanwhile, depotentiation, which refers to reversal of synaptic strength from potentiated LTP state to its pre-LTP level, is required by extinction of the obsolete memory. In the amygdala, cued fear-conditioning induces widespread synaptic strengthening (McKernan and Shinnick-Gallagher, 1997; Rogan et al., 1997; Rumpel et al., 2005). There is growing evidence to support that fear memories are transiently susceptible to erasure due to the ability of the same experience to reverse this form of synaptic strengthening (Kim et al., 2007; Clem and Huganir, 2010; Díaz-Mataix et al., 2011). A similar role of depotentiation is also found in hippocampal-dependent memory extinction (Zhang et al., 2011; Wang and Zhang, 2012). Novelty acquisition is not only able to induce LTP (Li et al., 2003; Gruart et al., 2006; Whitlock et al., 2006; Madroñal et al., 2007), but also enables induction of depotentiation from previously established LTP (Xu et al., 1998; Manahan-Vaughan and Braunewell, 1999; Abraham et al., 2002; Straube et al., 2003; Collingridge et al., 2010). These findings indicate that experience or stimulation protocol resulting in LTP can further induce depotentiation. However, mechanisms for depotentiation induced by further application of stimulations for LTP induction still remain unclear.
There have been a number of substantial theoretical investigations of the multi-stable mechanism of phosphatase and kinase switches which could account for LTP and LTD. Early simulation works proposed that the bistable switch formed by kinase, Ca2+/CaM protein kinase II (CaMKII), autophosphorylation underlies the stable transition between the basal and LTP states (Lisman, 1989). Recently, the authors further show that autodephosphorylation of protein phosphatase 2A (PP2A), one of the few serine/threonine-specific phosphatases, also forms a bistable switch, which could account for the stable transition between the basal and LTD states. Even more importantly, coupling of the phosphatase and kinase switches produces a tristable system, which exclusively explains how moderate and high elevation of intracellular calcium concentration triggers LTD and LTP, respectively. A phase-plane analysis was employed to understand the model graphically, in which nullclines at different calcium concentrations were plotted to show how the system moves during LTD and LTP induction and how it is finally stabilized after the induction, respectively (Pi and Lisman, 2008). This tristable system provides a generic theoretical framework for understanding dynamic mechanisms for LTP and LTD. Depotentiation induced by LTP induction protocol from the established LTP state to its basal state, which may be triggered by even higher elevation of intracellular concentration, also needs to be studied theoretically.
In a number of experiments, it has been shown that biochemical machinery for depotentiation is also formed by the molecular processes responsible for LTP, including calcium influxes, activation of specific enzymes, and glutamate receptor trafficking. Blockage of depotentiation when the induction stimulation was delivered in the presence of an antagonist of N-methyl-D-aspartic acid (NMDA) receptors (NMDARs) suggests that activation of NMDA receptors is required (Christie et al., 1995; Abraham and Huggett, 1997). Activation of metabotropic glutamate receptors (mGluRs) has also been shown required for depotentiation (Bikbaev et al., 2008; Qi et al., 2013). Besides membrane permeability (activation of Ca2+ channels), the amount of Ca2+ influx depends on the calcium ion concentration difference in the intracellular vs. extracellular spaces, which causes diffusion of calcium ion from high to low concentrations. So, in experiment, lowering extracellular calcium ion is applied as one way to reduce influx and intracellular calcium ion elevation. Stimulation is no longer able to induce depotentiation when it is delivered in a condition in which extracellular Ca2+ concentration is lower than normal (Abraham and Huggett, 1997), suggesting elevation of intracellular Ca2+ concentration is also required. In addition, phosphatase inhibitor can block depotentiation, showing involvement of protein phosphatase pathways in depotentiation induction (Kang-Park et al., 2003). The above experimentations strongly suggest that, after the establishment of LTP, further increase of intracellular calcium concentration, which could be initiated by calcium influxes mediated by NMDARs and mGluRs, may regulate phosphatase and kinase pathways and then switch synaptic strength back to its pre-LTP level.
In the present study, calcium influxes through NMDARs and calcium release from internal calcium stores by activation of mGluRs are introduced to the tristable system to simulate the experimentally observed depotentiation by using LTP-induction protocol. An phase-plane analysis of the dynamics for transition from LTP to its basal level is also provided.
Model and Methods
Our model contains two primary ingredients as schematized in Figure 1. The first is a mechanism through which the plasticity of synaptic strength can be depicted. Here, we adopt the tristable biochemical network described by Pi and Lisman (2008), which describes coupling of the CaMKII switch and PP2A switch in the postsynaptic site of synapses. The system yields three different stable states, corresponding to the basal, the LTP, and the LTD levels of synaptic strength, respectively. Transient moderate elevation of calcium concentration can induce LTD, while high elevation induces LTP. In this model, the values of phosphorylated CaMKII and dephosphorylated PP2A govern α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) receptors (AMPARs) trafficking which determines the number of AMPAR inserted in the postsynaptic membrane and hence the synaptic efficacy. A basal state (with basal AMPARs level) is determined by low concentrations of both phosphorylated CaMKII and dephosphorylated PP2A, a LTP state (with high AMPARs level) by high phosphorylated CaMKII and low dephosphorylated PP2A, and a LTD state (with low AMPARs level) by low phosphorylated CaMKII and high dephosphorylated PP2A, respectively. Switching between these states is controlled by cytoplasmic calcium concentration in the post synaptic site. These states are self-sustaining; once the system moves into one of these states, it remains there, although calcium concentration returns to the basal level.
Figure 1. Schematic diagram of the full model. Gray shading represents spine structure and light yellow shading represents internal stores [endoplasmic reticulum (ER)]. Pointed and circled arrows indicate the activation and the inhibition pathways, respectively. An agonist such as neurotransmitter binds to NMDARs and mGluRs. The latter one initiates a series of reactions, linked through G-protein, that ends in the production of the second messenger inositol trisphosphate (IP3), which diffuses through the cytoplasm and binds to IP3 receptors (IP3Rs) locating at ER membrane. IP3Rs are calcium channels binding with IP3 and mediating calcium release from the ER into cytoplasm. NMDARs are membrane calcium channels binding with glutamate mediating calcium influx. Cytoplasmic Ca2+ activates CaMKII (K) and protein phosphatase (P) in a manner that depends on the concentration of cytoplasmic Ca2+. CaMKII and phosphatase can self-activate themselves and inhibit each other. Active kinase and phosphatase control AMPARs insertion and removal, respectively.
The second key part of our model is a description of the cytoplasmic calcium dynamics in the post synaptic site, including both calcium accumulation and calcium extrusion. Calcium accumulation in cytoplasm is mainly dependent on calcium influxes through NMDARs and calcium release from internal calcium store [for example, in endoplasmic reticulum (ER)].
Postsynaptic calcium dynamics in cytoplasm and in stores are expressed by Kusters et al. (2005):
Here, INMDAR is the calcium current through NMDARs, JIP3R denotes calcium flux from the ER to the cytoplasmic space through IP3Rs, JSERCA is the calcium flux pumped from the cytoplasmic space into the ER, JIKER is the leak of calcium ions from the ER to cytoplasmic space, ZCa is calcium ion valence, F is Faraday's constant, VolCyt is the volume of synapses, VolER is the volume of stores, AER is the surface area of the store exposed in synapses, and are calcium concentrations of ER and cytoplasmic region, respectively, is the basal concentration of calcium in cytoplasmic region and the decay constant is τCa. The values of the different parameters in equations (1, 2) using in this paper are furnished in the Table 1.
INMDAR is given by Shouval et al. (2002)
where V is the postsynaptic potential and glutamate is released at a time t = tglut. Here P0 is the fraction of NMDARs that change to the open state after glutamate binding and GNMDAR is the conductance of the open NMDARs from calcium ions. θ1 is a step function: θ1(t) = 1 when t ≥ tglut, and θ1(t) = 0 when t < tglut. Glutamate unbinding is characterized by two decay constants and , and the fast and slow unbinding components and . The term H(V) describes Mg2+ unblock due to changes in the postsynaptic potential and is modeled by
where ECa is the reversal potential of calcium, V is the postsynaptic membrane potential, and [Mg2+] is the external magnesium ion concentration. All parameter values used for simulation of INMDAR are listed in Table 2.
The expressions for each calcium fluxes from the store are
And w is the fraction of activated IP3Rs; its dynamics is modeled by
IP3 is used to mimic activation of mGluRs, we describe the change in IP3 by production and degradation as follows (De Young and Keizer, 1992):
Where [IP3] is concentration of IP3 in cytoplasm, Ir is the rate constant for degradation of IP3, Ip is the maximal production generated by each activation of mGluRs, and f(t) controls the timing of activation of mGluRs and is described by a step function: If tglut + 2 ms > t > tglut, f(t) = 1; otherwise f(t) = 0. Details of the parameters are listed in Table 3.
Postsynaptic potential change has been modeled by (Shouval et al., 2002)
where Vrest is the resting membrane potential, Vrest = −65 mV. EPSP is the net depolarized voltage generated by binding glutamate to receptors in post synapse. The EPSP is:
where s is the maximal potential of depolarization, s = 50 mV, τ1 and τ2 are membrane potential decay constants, τ1 = 50 ms and τ2 = 5 ms, and tglut is the time that glutamate released.
Pi and Lisman Model of the CaMKII/PP2A Reaction Network
In the present modeling, the CaMKII/PP2A reaction network exactly follows the work of Pi and Lisman (2008). The cytoplasmic calcium dynamics in the post synaptic site affects the rates of autocatalytic reactions that convert CaMKII and PP2A from their non-phosphorylated to phosphorylated forms. For CaMKII the concentration of the phosphorylated form, which is the active form, is denoted by K*, while for PP2A the non-phosphorylated form is active, and its concentration is denoted by P*. All equations describing the evolution of K* and P* were taken directly from Pi and Lisman (2008):
where k1, k2, k3, k4, k11, k12, k13, and k14 are rate constants for each reaction. Km1, Km2, Km, Km11, and Km12 are equilibrium constants for kinase and phosphatase. K0 and P0 indicate the basal concentration of active kinase and phosphatase. Ktot and Ptot indicate the total amount of kinase and phosphatase.
Dependence of AMPARs on K* and P* is described as:
where A and Atot indicate the AMPARs on the synaptic membrane and the total amount. c1 and c2 are scaling factors, c3 and c4 are rate constants for processes independent of kinase and phosphatase activities. All parameters of CaMKII/PP2A reaction network are listed in Table 4.
All the calculations were implemented in the C programming language, with the differential equations expressed as different equations, and solved by using the fourth order Runge–Kutta method in steps of 0.01 ms. According to the related experimental procedures, theta-burst stimulation (TBS) was used as triggers for elevation of calcium concentration in cytoplasm. One unit of TBS is consisted of 3 pulses (100 Hz) delivered at 200 ms intervals. For application of each specific units of TBS, the variation of cytoplasmic calcium concentration in the post-synaptic site, the values of concentration of active kinase and phosphatase, and synaptic efficacy were calculated as functions of time, respectively.
Transitions of Synaptic Strength Induced by CaMKII/PP2A Network Switching
To better understand the transition of synaptic strength, it is useful to consider how the original Pi and Lisman model switches in simpler situations. We begin with a study on the responses of the original CaMKII/PP2A network to a simple calcium pulse (PCa) with different amplitude and the same duration (100 ms). Simulation results are shown in Figure 2. The results for PCa = 4.5 μM are consistent with those reported by Pi and Lisman (2008), changes of phosphatase and kinase activities are shown in Figure 2A. Both enzymes were nearly inactive in the basal state. During elevation of PCa, the kinase and phosphatase become strongly activated. At the end of the induction period, PCa returns to the basal level, phosphatase activation returns back to its basal level subsequently but the kinase activation remains at a high level. In this way a LTP is induced.
Figure 2. Simulation results for original pi and lisman model of coupled kinase and phosphatase swithes in response to different simple calcium pulse. K* and P* denote concentration of phosphorylated kinase and dephosphorylated phosphatase, respectively. (A) Changes of K* and P* in response to application of 4.5 μM calcium pulse (PCa) for 100 ms. Before PCa stimulation, K* and P* both stay at basal level. During PCa application, K* and P* rise, and K* is dominate. After PCa is over, P* returns back to basal level while K* stays at high level, which finally leads to LTP. (B) Changes of K* and P* in response to application of 6.5 μM calcium pulse for 100 ms. Both K* and P* are transiently increased by PCa and return back to basal level after Ca2+ is removed. Therefore, synaptic efficacy remains unchanged. Red line represents K* and blue line for P*, respectively. (C) Dependence of synaptic efficacy on amplitudes of PCa. As the level of calcium increase from moderate to high, synaptic efficacy changes from LTD to LTP, then return to basal state when amplitude of PCa increases further. (D) Reversals of LTP. LTP is induced by the first stimulation (4.5 μM, 100 ms) and reversed to basal level (depotentiation) by the second stimulation (6.5 μM, 100 ms). Changes in PCa is shown with red line and synaptic efficacy with black line, respectively.
However, further transient elevations of calcium concentration failed to induce a LTP state. The simulation result of changes in phosphatase and kinase activities to a transient increase of PCa to 6.5 μM is shown in Figure 2B. Similarly, two enzymes are inactive in the basal state. During PCa elevation, the kinase and phosphatase also become strongly activated. At the end of the induction period, PCa returns to its basal level and both phosphatase and kinase return to their basal states subsequently. So there is no sustained change in synaptic strength.
To further characterize the transitions between different states, we systematically varied amplitude of PCa with the duration fixed at 100 ms. For PCa between 0.1 and 1.3 μM, the system is stabled at basal state. For PCa in a range between 1.4 and 2.4 μM, the system is switched from the basal state to LTD states. For PCa between 3.0 μM and 6.2 μM, it is switched to LTP states. However, in conditions with PCa > 6.2 μM, neither LTP nor LTD is inducible (Figure 2C), the system returns to its basal state, suggesting that significant elevation of calcium concentration may depotentiate synaptic strength to basal levels.
Depotentiation from a LTP state was also tested by inputting a large PCa at a LTP level. As shown in Figure 2D, after a LTP is induced from the basal state by a proper PCa (4.5 μM for 100 ms), a large PCa (6.5 μM for 100 ms) reverses the potentiated state to its basal level, this mechanism of depotentiation is different from the ones induced by small elevations of PCa (Pi and Lisman, 2008), which induce LTD from a basal synaptic strength.
The Amount of TBS Trains Determines Changes in Cytoplasmic Calcium Concentration, K*, and P*
In order to directly simulate the experimental induction of depotentiation from a LTP state, we applied the induction procedure used in the experiments to the modified tristable system (including calcium dynamics) in simulation. A postsynaptic cytoplasmic calcium concentration dynamics produced by NMDARs activation and calcium release from internal calcium stores was considered and coupled in the original system. We then observed very different responses induced by application of different amount of TBS.
We first consider the responses to a brief TBS. A train of corresponding elevations of cytoplasmic calcium concentration is induced, and the amplitudes of these elevations are moderate (Figure 3A1). Details of the result are presented in Figure 3A2. Before application of the stimulation, no transmitters are released, thus NMDARs are closed, while the basal concentration of IP3 is very low. The IP3Rs are also closed, therefore, there are neither Ca2+ influx through NMDARs nor Ca2+ release through IP3Rs. During TBS delivery, NMDARs are activated to mediate a series of Ca2+ influx through NMDARs (Figure 3A2 top), while IP3Rs is still closed (Figure 3A2 bottom) as the concentration of IP3 is not high enough, although it shows growing pulsatile elevations (Figure 3A2 middle). The results show that, a comparatively high elevation in cytoplasmic calcium concentration is induced by brief TBS, which is mainly mediated by activation of NMDARs.
Figure 3. The amount of TBS trains determines the changes in cytoplasmic calcium concentration, K*, and P*. Cytoplasmic calcium concentration in post-synapse induced by a brief TBS (4 trains) (A1) and a prolonged TBS (24 trains) (B1). Details of each calcium fluxes induced by the brief (A2) and the prolonged TBS (B2). Enzyme activation induced by the brief TBS (4 trains) (A3) and the prolonged TBS (B3). Before the brief TBS, both K* and P* are nearly inactive. During stimulation, both K* and P* are transiently increased. After calcium removal, K* stays at high and P* returns back to basal level. Therefore, LTP is induced. On the other hand, the prolonged TBS induces a depotentiation.
As shown in Figure 3A3, before elevation of cytoplasmic calcium concentration, both enzymes are nearly inactive in the basal state. During the elevation of cytoplasmic calcium concentration, kinase and phosphates become strongly activated. Then cytoplasmic calcium concentration decline to the basal level, phosphates returns back to the basal level subsequently, while the kinase still stays at the activated state. Therefore, synaptic efficacy is switched to a LTP state.
Then, we turn to consider the responses to a prolonged TBS. Correspondingly, a series of elevations of cytoplasmic calcium concentration is induced by this prolonged TBS, and the amplitudes of the elevations gradually increase to significantly high levels (Figure 3B1). We then analyzed in details the sources of Ca2+. Before the stimulation, both NMDARs and IP3Rs are closed and there is no Ca2+ influx. As this prolonged TBS is delivered, NMDARs are activated to mediate stronger Ca2+ influx (Figure 3B2 top), meanwhile, IP3Rs are gradually activated (Figure 3B2 bottom) as the concentration of IP3 is accumulated to a higher lever (Figure 3B2 middle). Therefore, a more significant elevation in cytoplasmic calcium concentration is eventually formed by activation of both NMDARs and IP3Rs.
Changes in responses of phosphates and kinase caused by this significant elevation in cytoplasmic calcium concentration are shown in Figure 3B3. Both enzymes are nearly inactive in the basal state. During this significant elevation in cytoplasmic calcium concentration, the kinase and phosphates become strongly activated. However, both of them return back to their basal states subsequently as cytoplasmic calcium concentration declines to basal levels. Therefore, synaptic efficacy remains at the basal state.
Reversal of an Established LTP
Previous behavioral experiments have also shown that brief novelty acquisition was insufficient to reverse LTP, but a prolonged novelty exposure did reverse an established LTP (Qi et al., 2013). Similar phenomena were also observed in hippocampus slices by using electrical stimulations. In the CA1 region of hippocampus slices from young adult rats, TBS, which mimics a physiologically relevant frequency of neuronal activity exhibited in the hippocampus of behaving animals (Kandel and Spencer, 1961; Ranck, 1973; Bikbaev et al., 2008), induces long-lasting reversal of previously induced LTP (Barr et al., 1995). To test for whether a prolonged TBS can reverse LTP while a brief TBS cannot, we first simulated the effect of a brief TBS on established LTP. As shown in Figure 4A, activation of kinase was initially set at an activated state and phosphates at a basal state, corresponding to a LTP state of synaptic strength. When an elevation of cytoplasmic calcium concentration is induced by a brief TBS, kinase becomes slightly inactivated and phosphatase strongly activated. After cytoplasmic calcium concentration declines to the basal state, kinase returns to its initial activated state and phosphatase returns to its basal state. In other words, before and after brief TBS, activation levels of both kinase and phosphatase are unchanged. Hence, a brief TBS is not able to induce a transition from the established LTP state to the basal state.
Figure 4. Effect of different TBS on an established LTP. Changes of kinase (K*) and phosphatase (P*) (upper) and calcium dynamics (bottom) induced by the brief TBS (A) and the prolonged TBS (B).
Then, we simulated the effect of a prolonged TBS on an established LTP. The results are presented in Figure 4B. Before application of TBS, activation of kinase was still set at an activated state and phosphatase at a basal state. During a significant elevation by a prolonged TBS, kinase becomes strongly inactivated but phosphatase becomes strongly activated. When cytoplasmic calcium concentration declines to its basal level, activated forms of both kinase and phosphatase are switched to their basal states. The results clearly show that a prolonged TBS can reverse an established LTP by switching kinase from its activated state to its basal state, which is induced by a significant elevation of calcium concentration at the post synaptic site.
Phase Space Analysis
In order to show the dynamics of depotentiation graphically, the nullclines of CaMKII/PP2A network at different Ca2+ concentrations is given in Figure 5. Firstly, we studied the case with unchanged synaptic efficacy after action of Ca2+ pulses. In this case, the system stays at the basal state before application of Ca2+ pulses. At the basal Ca2+ level, nullclines (Figure 5A, left panel) intersects at five points; three are stable points (filled circles at top left, bottom left, and bottom right) and two are unstable (open circles). In the basal state, the system is stable at the bottom left point. During a high Ca2+ elevation, the nullclines are deformed to create only one stable state (Figure 5A, middle panel, light blue and light red lines, and gray dot), and the system finally stays at the new unique stable state (Figure 5A, middle panel, gray trace). After action of Ca2+ pulses, the three stable states restore, but the system stays at the top left stable point (LTP) (Figure 5A, top right panel). However, during a significant elevation of Ca2+ concentration, nullclines are deformed to create another unique stable state (Figure 5A, middle panel, blue and red lines, and black dot), and the system moves from this stable state (Figure 5A, middle panel, black trace) to and stays at the bottom left stable point (basal state), when the original three stable points restores after action of Ca2+ pulses (Figure 5A, bottom right panel).
Figure 5. Phase space analysis of the system's dynamics when initial state is at basal state (A) and at LTP (B), respectively. Left panel, red and blue curves indicate nullclines for kinase and phosphatase, respectively. Nullclines at basal Ca2+ concentration (0.1 μM) create five steady-states. Open and filled circles denote unstable and stable steady-states, the filled circle with black open circle indicates the basal state (A) and LTP (B), K* and P* denote concentration of active kinase and phosphatase, respectively. Middle panels, nullclines during Ca2+ influx. Ca2+ elevation deforms the nullclines to create one stable state. Red and blue lines indicate nullclines for K* and P* and black dot denotes the new stable state under 6.5 μM Ca2+ elevation. Light red and light blue lines indicate nullclines for K* and P* and gray dot denotes the new stable state under 4.5 μM Ca2+ elevation. The system becomes unstable and moves to the new stable states (gray trace for 4.5 μM and black trace for 6.5 μM) from basal state (A) and LTP (B), respectively. Right panels, nullclines after Ca2+ removal. Nullclines form three stable states again in the basal Ca2+ level. The state moves to the closest stable state, LTP (gray trace, top), or basal state (black trace, bottom), respectively.
Then, we studied reversal of the established LTP by application of Ca2+ pulses. At the LTP level, the system was stable at the top left stable point (Figure 5B, left panel). During action of high Ca2+ pulses, high Ca2+ elevations deform the nullclines to create one stable state (Figure 5B, middle panel, light blue and light red lines, and gray dot), and the system moves to this new stable state (Figure 5B, middle panel, gray trace). After action of Ca2+ pulses, nullclines restore the three stable states and the system returns to the top left stable point (LTP) (Figure 5B, top right panel). During action of significant Ca2+ pulses, large Ca2+ elevations deform nullclines to create the unique stable state (Figure 5B, middle panel, blue and red lines, and black dot), and the system moves to this new stable state (Figure 5B, middle panel, black trace). After action of Ca2+ pulses, nullclines restore the three stable states, but the system moves to and stays at the bottom left stable point (basal state) (Figure 5B, bottom right panel).
Attractive basins of each stable state were depicted in the K* − P* plane when calcium concentration was set at the basal level (Figure 6). We have known that nullclines form three stable states (filled circles at top left, bottom left, and bottom right of left panels in Figures 5A,B) at basal Ca2+ level. The K* − P* plane is then divided into three regions; points located within the red region will evolve to the stable point at top left (named LTP basin), those in the green region to the stable point at bottom left (named basal state basin), and those in the blue region to the one point at bottom right (named LTD basin), respectively.
Figure 6. Attractive basins with basal level Ca2+ and evolvement of the stable state induced by adjusting Ca2+ elevation. Red region indicates the attractive basin of LTP, Green for basal state, and blue for LTD. Black dots indicate evolvement of the stable state. Two white dots represent the stable states induced by 4.5 μM (bottom left) and 6.5 μM (top right), respectively (also shown in the insertion).
Then we turned to observe the distribution of the stable points created by Ca2+ elevations in the K* − P* plane. As the amplitude of Ca2+ elevation increases from 3 to 10 μM gradually, the stable points induced by Ca2+ elevation moves from bottom left to top right, from the LTP basin to the basal state basin. The stable point induced by 4.5 μM Ca2+ elevation is in the LTP basin (white dot in left), while the stable point induced by 6.5 μM Ca2+ elevation is in the basal state basin (white dot in right, graph inserted in Figure 6). After the action of Ca2+ elevation, nullclines restore the three stable states, and the system will evolve from its new initial position (the position determined by Ca2+ elevation) to one of the three stable states, according to the location of this initial position in the K* − P* plane. The new initial position determined by high Ca2+ elevation is in the LTP basin and the system evolves to LTP. The new initial position determined by significant Ca2+ elevation is in the basal state basin, therefore, the system evolves to the basal state and can manifest in this way the depotentiation from a LTP state. Just as in the original CaMKII/PP2A network model, Ca2+ elevation is still the primary controller, which controls the synaptic strength through activation of kinase and phosphatase. In experimentations synaptic strength is not directly manipulated by Ca2+ elevation but by presynaptic simulation, including TBS. The above analysis provides a link between experiments and the generic dynamics revealed by the tristable system (Pi and Lisman, 2008).
Although there has been substantial investigation on how a tristable system composed of coupled kinase and phosphatase switches could account for stable transition of synaptic strength from basal state to both LTP and LTD, depotentiation from LTP to basal state has not been theoretically studied. Based on a series of experimental observation, the present work introduced two Ca2+ resources with different temporal dynamics, mediated by influx through NMDARs and release from internal calcium store, respectively, in the generic tristable system and further show that significant elevation of cytoplasmic calcium concentration may dynamically switch activation of both kinase and phosphatase to their basal states, thereby depotentiate the synaptic strength. The generation of this significant elevation of cytoplasmic calcium concentration depends on the length of stimulation, a prolonged stimulation activates both NMDARs and internal Ca2+ stores and then produces a significant Ca2+ elevation, while a brief stimulation only activates NMDARs and induces a high Ca2+ elevation. When initial condition of synaptic strength is at basal state, a prolonged TBS dynamically firstly induces a high Ca2+ elevation which might result in LTP, and then induces a large Ca2+ elevation and depotentiates the possible LTP.
Involvement of internal calcium store is required for induction of depotentiation. Both experiments and simulations have shown that experience exposure or stimulation protocol elevating cytoplasmic calcium concentration by activating cellular membrane calcium channels can switch the synaptic strength from basal level to LTP. We have shown in this paper that repeated experience exposure or prolonged stimulation could further activate postsynaptic group I mGluRs, produce more IP3 to activate IP3Rs on endoplasmic reticulum, finally lead to calcium release from internal calcium store. Relatively slow dynamics of this process requires repeated experience or prolonged stimulation. Recent experiment results show that in synapses containing endoplasmic reticulum, presynaptic stimulation induced two successive elevation peaks of calcium concentration: The first one is mediated by NMDARs and the second by IP3Rs activation, respectively (Holbro et al., 2009; Sheridan et al., 2014). These experiments support the simulation results of the present paper, not only by showing the necessity of a significant calcium elevation, but also by an equivalent time course for induction of depotentiation. The results in our modeling also suggest that depotentiation should have a threshold-like effect with calcium concentration. The prediction of this threshold and its changes with changes of internal calcium store needs to be demonstrated in further experiments.
In an interesting experiment, mutant mice lacking IP3Rs display a significantly greater magnitude of LTP induced by tetanus stimulation (Fujii et al., 2000). A similar mutant strain lacking ryanodine receptor type 3 (RYR3), another type of calcium channel on internal calcium store, exhibits in hippocampus CA1 pyramidal neurons facilitated LTP as well (Futatsugi et al., 1999). These results are controversial to some expectations for release of calcium from internal store facilitates LTP. Our simulation work could provide an explanation for these experiments, by indicating that lack of IP3Rs or RYR3 may impair calcium release, and thus prevent depotentioation, finally facilitate LTP.
Compared with the calcium elevation inducing LTP, more significant and prolonged elevation of cytoplasmic calcium concentration is required for depotentiation induction. It has been proposed that moderate elevation in [Ca2+]cyt that is produced during LTD induction may preferentially activate calcinerin, the Ca2+/calmodulin (CaM)-dependent protein phosphatase. Calcineurin then can induce inhibition of protein 1 (I-1), the endogenous inhibitor of serine/threonine protein phosphatase-1 (PP1). This allows the activation of PP1 (Oliver and Shenolikar, 1998). PP1 dephosphorylates CaMKII and other proteins (Shields et al., 1985), including glutamate receptors, to promote LTD. However, LTP-inducing stimuli are associated with higher increase in cytoplasmic calcium concentration. This can activate CaMKII, which plays a key role in LTP induction through phosphorylation of the GluR1 subunit of AMPARs (Cormier et al., 2001). In addition, high elevation of cytoplasmic calcium concentration can activate PKA, PKA can phosphorylate and activate I-1, which then can inhibit PP1 (Barria et al., 1997). Thus, LTP is induced and LTD is prevented.
PKA, which negatively regulates PP1, is modulated by the activation of Ca2+/CaM-dependent adenylyl cyclases (AC) (Piascik et al., 1980; Potter et al., 1980; Ahlijanian and Cooper, 1988). And the activity of AC, which shows a bell-shaped activity curve relative to increasing calcium concentrations (Piascik et al., 1980; Potter et al., 1980; Ahlijanian and Cooper, 1988), may be decreased by large rise in cytoplasmic calcium concentration induced by prolonged TBS. The resulting low PKA activity is accompanied by a disinhibition of PP1, which inactivates CaMKII, leading occurrence of depotentiation (Lisman, 1989).
When certain memory becomes obsolete, effective extinction of the previously established memory is essential for animals to adapt to the changing environment. An exposure-based therapy for human traumatic fear memories has long been established (McNally, 2007). Fear memories present as defensive reactions to the neutral cue for a period of up to many months after learning. However, when conditioned subjects repeatedly encounter a neutral cue without a reinforcing unpleasant event, a large rise in cytoplasmic calcium concentration decreases kinase activity and subsequently induces depotentiation by AMPA receptors endocytosis (Clem and Huganir, 2010). The present simulation results may account for the mechanism of these behavioral investigations.
It has been studied and clinically applied that agonist of NMDA receptors, D-cycloserine, could facilitate fear extinction when given systematically or locally into the amygdale (Davis, 2010). Targeting mGluRs and calcium pump may provide other mechanistic insights into potential treatments for fear extinction. Moreover, significant elevation of cytoplasmic calcium concentration resulting in depotentiation is achieved by enzyme reactions, hence modulations on the activities of enzymes become to the second kind of approaches. It is observed that ACs is required for spatial memory extinction (Zhang et al., 2011) and mice over-expressing type 1 ACs show enhancing spatial memory extinction (Zhang and Wang, 2013). Overall, the present findings uncover a potentially important time- and state-dependent mechanism of NMDAR- and mGluR-dependent depotentiation. Furthermore, the engagement of such depotentiation may be involved in the activity-dependent erasure of recently stored information in the hippocampus (Zhang et al., 2011; Wang and Zhang, 2012; Zhang and Wang, 2013), showing strong parallels with similar mechanisms engaged in the amygdala (Kim et al., 2007).
MC prepared the methods of feature construction and conducted the experiments. WR and MC prepared the manuscript. WR and XW supervised all aspects of the work. All authors discussed the results and commented on the manuscript.
This work was supported by the Fundamental Research Funds for the Central Universities of China (No. GK201503027 and GK201601001) and the National Natural Science Foundation of China (No. 11375109).
Conflict of Interest Statement
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.
Abraham, W. C., and Huggett, A. (1997). Induction and reversal of long-term potentiation by repeated high-frequency stimulation in rat hippocampal slices. Hippocampus 7, 137–145. doi: 10.1002/(SICI)1098-1063(1997)7:2<137::AID-HIPO3>3.0.CO;2-K
Abraham, W. C., Logan, B., Greenwood, J. M., and Dragunow, M. (2002). Induction and experience-dependent consolidation of stable long-term potentiation lasting months in the hippocampus. J. Neurosci. 22, 9626–9634.
Ahlijanian, M. K., and Cooper, D. M. (1988). Distinct interactions between Ca2+/calmodulin and neurotransmitter stimulation of adenylate cyclase in striatum and hippocampus. Cell. Mol. Neurobiol. 8, 459–469. doi: 10.1007/BF00711229
Barr, D. S., Lambert, N. A., Hoyt, K. L., Moore, S. D., and Wilson, W. A. (1995). Induction and reversal of long-term potentiation by low- and high-intensity theta pattern stimulation. J. Neurosci. 15, 5402–5410.
Barria, A., Muller, D., Derkach, V., Griffith, L. C., and Soderling, T. R. (1997). Regulatory phosphorylation of AMPA-type glutamate receptors by CaM-KII during long-term potentiation. Science 276, 2042–2045. doi: 10.1126/science.276.5321.2042
Bikbaev, A., Neyman, S., Ngomba, R. T., Conn, P. J., Nicoletti, F., and Manahan-Vaughan, D. (2008). MGluR5 mediates the interaction between late-LTP, network activity, and learning. PLoS ONE 3:e2155. doi: 10.1371/journal.pone.0002155
Christie, B. R., Stellwagen, D., and Abraham, W. C. (1995). Reduction of the threshold for long-term potentiation by prior theta-frequency synaptic activity. Hippocampus 5, 52–59. doi: 10.1002/hipo.450050107
Cormier, R. J., Greenwood, A. C., and Connor, J. A. (2001). Bidirectional synaptic plasticity correlated with the magnitude of dendritic calcium transients above a threshold. J. Neurophysiol. 85, 399–406.
De Young, G. W., and Keizer, J. (1992). A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca2+ concentration. Proc. Natl. Acad. Sci. U.S.A. 89, 9895–9899. doi: 10.1073/pnas.89.20.9895
Díaz-Mataix, L., Debiec, J., LeDoux, J. E., and Doyère, V. (2011). sensory-specific associations stored in the lateral amygdala allow for selective alteration of fear memories. J. Neurosci. 31, 9538–9543. doi: 10.1523/JNEUROSCI.5808-10.2011
Fujii, S., Matsumoto, M., Igarashi, K., Kato, H., and Mikoshiba, K. (2000). Synaptic plasticity in hippocampal CA1 neurons of mice lacking type 1 inositol-1,4,5-trisphosphate receptors. Learn. Mem. 7, 312–320. doi: 10.1101/lm.34100
Futatsugi, A., Kato, K., Ogura, H., Li, S. T., Nagata, E., Kuwajima, G., et al. (1999). Facilitation of NMDAR-independent LTP and spatial learning in mutant mice lacking ryanodine receptor type 3. Neuron 24, 701–713. doi: 10.1016/S0896-6273(00)81123-X
Gruart, A., Muñoz, M. D., and Delgado-García, J. M. (2006). Involvement of the CA3-CA1 synapse in the acquisition of associative learning in behaving mice. J. Neurosci. 26, 1077–1087. doi: 10.1523/JNEUROSCI.2834-05.2006
Holbro, N., Grunditz, A., and Oertner, T. G. (2009). Differential distribution of endoplasmic reticulum controls metabotropic signaling and plasticity at hippocampal synapses. Proc. Natl. Acad. Sci. U.S.A. 106, 15055–15060. doi: 10.1073/pnas.0905110106
Kang-Park, M. H., Sarda, M. A., Jones, K. H., Moore, S. D., Shenolikar, S., Clark, S., et al. (2003). Protein phosphatases mediate depotentiation induced by high-intensity theta-burst stimulation. J. Neurophysiol. 89, 684–690. doi: 10.1152/jn.01041.2001
Kusters, J. M., Dernison, M. M., van Meerwijk, W. P., Ypey, D. L., Theuvenet, A. P., and Gielen, C. C. (2005). Stabilizing role of calcium store-dependent plasma membrane calcium channels in action-potential firing and intracellular calcium oscillations. Biophys. J. 89, 3741–3756. doi: 10.1529/biophysj.105.062984
Li, S., Cullen, W. K., Anwyl, R., and Rowan, M. J. (2003). Dopamine-dependent facilitation of LTP induction in hippocampal CA1 by exposure to spatial novelty. Nat. Neurosci. 6, 526–531. doi: 10.1038/nn1049
Madroñal, N., Delgado-García, J. M., and Gruart, A. (2007). Differential effects of long-term Potentiation evoked at the CA3-CA1 synapse before, during, and after the acquisition of classical eyeblink conditioning in behaving mice. J. Neurosci. 27, 12139–12146. doi: 10.1523/JNEUROSCI.3397-07.2007
Manahan-Vaughan, D., and Braunewell, K. H. (1999). Novelty acquisition is associated with induction of hippocampal long-term depression. Proc. Natl. Acad. Sci. U.S.A. 96, 8739–8744. doi: 10.1073/pnas.96.15.8739
Pi, H. J., and Lisman, J. E. (2008). Coupled phosphatase and kinase switches produce the tristability required for long-term potentiation and long-term depression. J. Neurosci. 28, 13132–13138. doi: 10.1523/JNEUROSCI.2348-08.2008
Potter, J. D., Piascik, M. T., Wisler, P. L., Robertson, S. P., and Johnson, C. L. (1980). Calcium dependent regulation of brain and cardiac muscle adenylate cyclase. Ann. N. Y. Acad. Sci. 356, 220–231. doi: 10.1111/j.1749-6632.1980.tb29613.x
Qi, Y. J., Hu, N. W., and Rowan, M. J. (2013). Switching off LTP: mGlu and NMDA receptor-dependent novelty exploration-induced depotentiation in the rat hippocampus. Cereb. Cortex. 23, 932–939. doi: 10.1093/cercor/bhs086
Sheridan, G. K., Moeendarbary, E., Pickering, M., O'Connor, J. J., and Murphy, K. J. (2014). Theta-burst stimulation of hippocampal slices induces network-level calcium oscillations and activates analogous gene transcription to spatial learning. PLoS ONE 9:e100546. doi: 10.1371/journal.pone.0100546
Shields, S. M., Ingebritsen, T. S., and Kelly, P. T. (1985). Identification of protein phosphatase 1 in synaptic junctions: dephosphorylation of endogenous calmodulin-dependent kinase II and synapse-enriched phosphoproteins. J. Neurosci. 5, 3414–3422.
Shouval, H. Z., Bear, M. F., and Cooper, L. N. (2002). A unified model of NMDA receptor-dependent bidirectional synaptic plasticity. Proc. Natl. Acad. Sci. U.S.A. 99, 10831–10836. doi: 10.1073/pnas.152343099
Straube, T., Korz, V., and Frey, J. U. (2003). Bidirectional modulation of long-term potentiation by novelty-exploration in rat dentate gyrus. Neurosci. Lett. 344, 5–8. doi: 10.1016/S0304-3940(03)00349-5
Zhang, M., Storm, D. R., and Wang, H. B. (2011). Bidirectional synaptic plasticity and spatial memory flexibility require Ca2+-stimulated adenylyl cyclases. J. Neurosci. 31, 10174–10183. doi: 10.1523/JNEUROSCI.0009-11.2011
Keywords: long-term potentiation, depotentiation, tristability, kinase, phosphatase
Citation: Chen M, Ren W and Wang X (2016) Depotentiation from Potentiated Synaptic Strength in a Tristable System of Coupled Phosphatase and Kinase. Front. Comput. Neurosci. 10:104. doi: 10.3389/fncom.2016.00104
Received: 28 June 2016; Accepted: 26 September 2016;
Published: 19 October 2016.
Edited by:David Holcman, École Normale Supérieure, France
Copyright © 2016 Chen, Ren and Wang. 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) or licensor 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.