Modelling the dynamics of CaMKII–NMDAR complex related to memory formation in synapses: The possible roles of threonine 286 autophosphorylation of CaMKII in long term potentiation
Introduction
Recent research shows that memory is stored in a set of synapses as synaptic strength, which strongly correlates to the magnitude of the postsynaptic response to an environmental stimulus (see supplementary material (SI) for the details of synaptic structure and synaptic transmission). The memory is retrieved by recognising the stored synaptic response patterns that match or are most sensitive to the patterns of the environmental stimulus (Bear et al., 2007, Martin et al., 2000). The synaptic adaptation to the external stimulus is experimentally linked to the long term potentiation (LTP): LTP is a synaptic process that long-lastingly increases synaptic strength, through increasing the excitatory postsynaptic potential (EPSP) or receptor mediated current (EPSC) in vivo, in response to high frequency stimulations (Bliss and Collingridge, 1993, Bliss and Lømo, 1973; and see SI). LTP and memory formation are strongly correlated: (1) both the memory formation and LTP trigger formation of new synapses (Engert and Bonhoeffer, 1999, Morgado-Bernal, 2011); (2) blocking LTP diminishes spatial memory formation (Davis et al., 1992, Morris, 1989, Tsien et al., 1996); and (3) the disruption of LTP by mutations impairs memory formation (Abel et al., 1997, Bourtchuladze et al., 1994, Chang et al., 1999, Grant and Silva, 1994, Huang and Lee, 1995, Silva et al., 1992a, Silva et al., 1992b). Therefore, the elucidation of the mechanisms and the pathways underlying LTP is imperative to understand synaptic adaptation and memory formation. In this paper, we investigate a potential pathway and the associated mechanisms to adapt synaptic strength using a novel mathematical model. The model is centred around the conformational state transitions (STs) of the “memory molecule”, Ca2+/Calmodulin (CaM) dependent protein kinase II, CaMKII (Lisman et al., 2012).
LTP contains two distinct phases: (1) the early phase (E-LTP) occurs immediately after the transient (a few seconds) induction signal and rapidly potentiates synaptic strength; and (2) the late phase (L-LTP) follows E-LTP to maintain the potentiation by triggering gene expression pathways and may change the synaptic structure (Abel et al., 1997, Bliss and Collingridge, 1993, Frey et al., 1993). Our focus here is modelling mechanisms related to E-LTP as the mechanisms related to L-LTP are not well understood.
The consensus of the experimental studies of LTP in CA1 region of hippocampus (Bear et al., 2007, Bliss and Collingridge, 1993, Kandel, 2009) is that LTP is induced by a postsynaptic protein cascade (Isaac et al., 1995, Malenka and Bear, 2004, Sheng and Kim, 2002, Xia and Storm, 2005; see SI for the description of the cascade) initiated by N-methyl-d-aspartate receptor (NMDAR). Once the high frequency stimulation is applied, NMDAR triggers a Ca2+ influx to initiate the protein cascade, which leads to the changes in A-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor (AMPAR)-mediated EPSC (Citri and Malenka, 2007, Kandel, 2009, Mayford et al., 2012).
CaMKII, a holoenzyme with 12 subunits (Lisman et al., 2012), phosphorylates AMPAR (see SI) to increase: (1) the conductance of the embedded ion channel of AMPAR (Derkach et al., 1999), and (2) the anchoring of AMPAR into the postsynaptic density (PSD) (Opazo et al., 2010, Tomita et al., 2005), where AMPAR becomes active; both of these increases will raise the AMPAR-mediated EPSC. Moreover, CaMKII can stay in PSD for more than 30 min (Fig. 2B, Bayer et al., 2006) that may be critical to cause the long-lasting increase in the EPSC and to link the two phases of LTP.
The binding between CaMKII and NMDAR serves a number of functions: (1) anchoring CaMKII in PSD (Gardoni et al., 1998, Leonard et al., 1999, Strack and Colbran, 1998); (2) targeting CaMKII to AMPAR since there is a linkage between NMDAR and AMPAR via A-kinase anchoring protein (Colledge et al., 2000, Robertson et al., 2009); and (3) longer sustaining periods of CaMKII in PSD (Bayer et al., 2006). Further, the formation of CaMKII–NMDAR complex is an essential component of LTP and memory formation (Barria and Malinow, 2005, Halt et al., 2012, Sanhueza et al., 2011, Otmakhov et al., 2004, Vest et al., 2007, Zhang et al., 2008, Zhou et al., 2007). Disparate experimental evidence exists on the significance of the autophosphorylation of CaMKII for the formation of CaMKII–NMDAR complex and for the induction of E-LTP. Some studies state that the formation of CaMKII–NMDAR complex does not require the autophosphorylation (Barcomb et al., 2013, Bayer et al., 2001, Bayer et al., 2006): T286A mutant (Bayer et al., 2006), which lacks the autophosphorylation, and the autophosphorylation repressors (Barcomb et al., 2013) do not block the formation of CaMKII–NMDAR complex, although the strength of CaMKII and NMDAR binding decreases. On the other hand, another set of studies show that the autophosphorylation is required for LTP and spatial memory formation (Giese et al., 1998, Irvine et al., 2006), especially during E-LTP (Buard et al., 2010). Furthermore, a study reports that the autophosphorylation is irreversible in PSD that may have a critical role in LTP (Mullasseril et al., 2007).
Therefore, we hope to understand the role of the autophosphorylation and the relationships among the autophosphorylation, CaMKII–NMDAR complex and the LTP by developing a dynamic model of the processes involved in the STs of CaMKII and its binding to NMDAR in response to Ca2+ signals. The main contributions of this model development are the inclusion of the STs based on the conformation of CaMKII, the inclusion of translocation of the holoenzyme between dendritic spine and PSD, and the development of a probabilistic framework for the binding between NMDAR and CaMKII holoenzyme. The previous models either considered only the STs of single subunits of CaMKII or did not include the binding to NMDAR (Chao et al., 2011, Chiba et al., 2008, Dosemeci and Albers, 1996, Dupont et al., 2003, Hayer and Bhalla, 2005, Kubota and Bower, 2001, Lisman, 1989, Lisman and Zhabotinsky, 2001, Miller et al., 2005, Pepke et al., 2010, Zhabotinsky, 2000). To the best of our knowledge, this is the first attempt to model the holoenzyme state transitions (HSTs) of CaMKII, including the binding to NMDAR. The modelling of HSTs is important to understand the dynamic nature of the formation of CaMKII–NMDAR complex as the experiments show (Strack et al., 2000) that the holoenzyme structure of CaMKII is crucial for its binding to NMDAR.
Our model provides insights into the formation of CaMKII–NMDAR complex which is regulated not only by the binding affinity between CaMKII and NMDAR as expected, but also by the rate of the CaMKII translocation into PSD. The model predicts well the alteration on the CaMKII translocation by T286A mutant in comparison to the experimental observation. In addition, the model predicts the formation of CaMKII–NMDAR complex under T286A mutant well: in particular, the model predicts that T286A mutant is not able to distinguish the frequency of tetanus; this prediction is consistent with the experimental observation. With regard to the autophosphorylation, the model shows that (1) the autophosphorylation amplifies the postsynaptic responses in reaction to the stimulation signals; (2) the autophosphorylation decodes the frequency of the stimulation to trigger appropriate postsynaptic response as reported frequently by previous modelling studies of CaMKII (Dosemeci and Albers, 1996, Kubota and Bower, 2001); and (3) the autophosphorylation can couple multiple trains of tetanus that are separated by relatively long inter-train intervals.
The paper is organised as follows: we briefly discuss the biochemical background needed for the model development in Section 2; Section 3 is about the previous models of ISTs; Section 4 is on the model development; and we discuss our results in Section 5, and we make concluding remarks in Section 6.
Section snippets
Biochemical background for modelling
A CaMKII subunit consists of a kinase domain, a regulatory domain, and a hub domain (see Lisman et al., 2002, Lisman et al., 2012 for an in-depth review and Chao et al., 2011 for the structural images of CaMKII). Twelve subunits are folded into two hexametric rings and each ring contains six subunits with their hub domains attached one by one to form a central hub (Fig. 1A) (Chao et al., 2011). The functionality of CaMKII relies on the following critical sites (Fig. 1B): (1) the substrate
Previous models
A number of earlier mathematical models describe the CaM dynamics (Holmes, 2000, Pepke et al., 2010) and the successive ISTs (Chao et al., 2011, Chiba et al., 2008). We use Holmes׳s model of the CaM dynamics (Holmes, 2000) and its successor (Chiba et al., 2008) to model ISTs of CaMKII in response to the Ca2+ influx. The deciding factor for this choice is the relative simplicity of Holmes׳s approach without compromising the agreement to the experimental data (Chiba et al., 2008). Holmes׳s model
A new model of HSTs of CaMKII
This section gives the details of a new mathematical model of HSTs of CaMKII, including the assumptions on which the model is based, the conceptual models, the equations and the estimated parameter values.
The autophosphorylation and CaM number in CaMKII translocation
As shown in Fig. 7A, the time course of the CaMKII translocation into PSD agrees well with the wild type data (no T286A mutant data found); the simulated time courses of the CaMKII dissociation from PSD have reasonable agreements to both the wild type and T286A mutant data (Shen and Meyer, 1999, Shen et al., 2000, Lee et al., 2009). The wild type CaMKII takes about 1 min to have complete CaM dissociation from CaMKII subunit (Fig. 7B) which is consistent with the experimentally reported time
Summary and conclusions
We develop a theoretical model of the formation of CaMKII–NMDAR complex which has the following novel contributions in comparison to the previously published studies: (1) the conformations of CaMKII subunits are expressed as ratios of the total number of CaMKII subunits and the compositions of CaMKII are formulated by binomial distribution based on these ratios; (2) the probability of the binding between CaMKII and NMDAR is dynamically changing with respect to the STs of CaMKII; (3) the time
References (97)
- et al.
Genetic demonstration of a role for PKA in the late phase of LTP and in hippocampus-based long-term memory
Cell
(1997) - et al.
Local calcium signaling in neurons
Neuron
(2003) - et al.
NMDA receptor subunit composition controls synaptic plasticity by regulating binding to CaMKII
Neuron
(2005) - et al.
Deficient long-term memory in mice with a targeted mutation of the cAMP-responsive element-binding protein
Cell
(1994) - et al.
A mechanism for tunable autoinhibition in the structure of a human Ca2+/calmodulin-dependent kinase II holoenzyme
Cell
(2011) - et al.
A simulation study on the activation of cardiac CaMKII< i> δ</i>-isoform and its regulation by phosphatases
Biophys. J.
(2008) - et al.
Targeting of PKA to glutamate receptors through a MAGUK–AKAP complex
Neuron
(2000) - et al.
A mechanism for synaptic frequency detection through autophosphorylation of CaM kinase II
Biophys. J.
(1996) - et al.
Sensitivity of CaM kinase II to the frequency of Ca2+ oscillations: a simple model
Cell Calcium
(2003) - et al.
Inhibitory autophosphorylation of CaMKII controls PSD association, plasticity, and learning
Neuron
(2002)
Studies of the regulatory mechanism of Ca2+/calmodulin-dependent protein kinase II. Mutation of threonine 286 to alanine and aspartate
J. Biol. Chem.
Targeting learning
Trends Neurosci.
Role of hippocampal nitric oxide in memory retention in rats
Pharmacol. Biochem. Behavior
αCaMKII autophosphorylation: a fast track to memory
Trends Neurosci.
Evidence for silent synapses: implications for the expression of LTP
Neuron
Robustness of G1/S checkpoint pathways in cell cycle regulation based on probability of DNA-damaged cells passing as healthy cells
Biosystems
A model of synaptic memory: a CaMKII/PP1 switch that potentiates transmission by organizing an AMPA receptor anchoring assembly
Neuron
LTP and LTD: an embarrassment of riches
Neuron
A methodology for performing global uncertainty and sensitivity analysis in systems biology
J. Theor. Biol.
Regulation of brain Type II Ca2+ calmodulin-dependent protein kinase by autophosphorylation: a Ca2+-triggered molecular switch
Cell
Sequences of autophosphorylation sites in neuronal type II CaM kinase that control Ca2+-independent activity
Neuron
Learning and memory consolidation: linking molecular and behavioral data
Neuroscience
CaMKII triggers the diffusional trapping of surface AMPARs through phosphorylation of stargazin
Neuron
CaMKIIβ functions as an F-actin targeting module that localizes CaMKIIα/β heterooligomers to dendritic spines
Neuron
Serine/threonine phosphatases in the nervous system
Curr. Opin. Neurobiol.
Autophosphorylation-dependent targeting of calcium/calmodulin-dependent protein kinase II by the NR2B subunit of theN-methyl-d-aspartate receptor
J. Biol. Chem.
Mechanism and regulation of calcium/calmodulin-dependent protein kinase II targeting to the NR2B subunit of the N-methyl-d-aspartate receptor
J. Biol. Chem.
Structural studies on the regulation of Ca2+/calmodulin dependent protein kinase II
Curr. Opin. Struct. Biol.
α-and βCaMKII: inverse regulation by neuronal activity and opposing effects on synaptic strength
Neuron
Bidirectional synaptic plasticity regulated by phosphorylation of stargazin-like TARPs
Neuron
The essential role of hippocampal CA1 NMDA receptor-dependent synaptic plasticity in spatial memory
Cell
Structural examination of autoregulation of multifunctional calcium/calmodulin-dependent protein kinase II
J. Biol. Chem.
Regulation of NMDA receptor subunit expression and its implications for LTD, LTP, and metaplasticity
Neuropharmacology
Bistability in the Ca2+/calmodulin-dependent protein kinase-phosphatase system
Biophys. J.
Enzymatic activity of CaMKII is not required for its interaction with the glutamate receptor subunit GluN2B
Mol. Pharmacol.
Interaction with the NMDA receptor locks CaMKII in an active conformation
Nature
Transition from reversible to persistent binding of CaMKII to postsynaptic sites and NR2B
J. Neurosci.
Neuroscience
Still Looking for the Memories: Molecules and Synaptic Plasticity. 20 Years of Computational Neuroscience
Long-lasting potentiation of synaptic transmission in the dentate area of the anaesthetized rabbit following stimulation of the perforant path
J. Physiol.
A synaptic model of memory: long-term potentiation in the hippocampus
Nature
CaMKII “autonomy” is required for initiating but not for maintaining neuronal long-term information storage
J. Neurosci.
Impaired memory retention and decreased long-term potentiation in integrin-associated protein-deficient mice
Learn. Memory
Synaptic plasticity: multiple forms, functions, and mechanisms
Neuropsychopharmacology
A significant but rather mild contribution of T286 autophosphorylation to Ca2+/CaM-stimulated CaMKII activity
PloS ONE
The NMDA receptor antagonist D-2-amino-5-phosphonopentanoate (D-AP5) impairs spatial learning and LTP in vivo at intracerebral concentrations comparable to those that block LTP in vitro
J. Neurosci.
Ca2+/calmodulin-kinase II enhances channel conductance of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionate type glutamate receptors
Proc. Natl. Acad. Sci.
Cited by (12)
Changed expressions of N-methyl-D-aspartate receptors in the brains of rats and primary neurons exposed to high level of fluoride
2018, Journal of Trace Elements in Medicine and BiologyNeural changes in Alzheimer's disease from circuit to molecule: Perspective of optogenetics
2017, Neuroscience and Biobehavioral ReviewsGlobal sensitivity analysis of a model related to memory formation in synapses: Model reduction based on epistemic parameter uncertainties and related issues
2017, Journal of Theoretical BiologyCitation Excerpt :If we are to reduce STM to a simpler one, we could use the information provided by the GSA and PRCC index in a biologically meaningful manner. Many important processes are involved in the mechanisms that connect the Ca2+ influx to the formation of CaMKII-NMDAR complex: Ca2+/CaM complex formation; activation of CaMKII subunit; autophosphorylation of CaMKII subunit; translocation of CaMKII holoenzyme into PSD; CaMKII holoenzyme binding to NMDAR at (S site of) a CaM bound subunit; CaMKII holoenzyme binding to NMDAR at (S site of) an autophosphorylated subunit; and S site bound CaMKII transferring to T site bound (see Figs. 2 and 4 in He et al. (2015)). These processes form building blocks of STM model.
Alzheimer's disease: Biology, biophysics and computational models
2022, Alzheimer's Disease: Biology, Biophysics And Computational Models