A spatial model of autophosphorylation of Ca2+/calmodulin-dependent protein kinase II in a glutamatergic spine reveals dynamics of kinase activation in the first several seconds after a complex synaptic stimulus.

Long-term potentiation (LTP) is a biochemical process that underlies learning in excitatory glutamatergic synapses in the Central Nervous System (CNS). The critical early driver of LTP is autophosphorylation of the abundant postsynaptic enzyme, Ca2+/calmodulin-dependent protein kinase II (CaMKII). Autophosphorylation is initiated by Ca2+ flowing through NMDA receptors activated by strong synaptic activity. Its lifetime is ultimately determined by the balance of the rates of autophosphorylation and of dephosphorylation by protein phosphatase 1 (PP1). Here we have modeled the autophosphorylation and dephosphorylation of CaMKII during synaptic activity in a spine synapse using MCell4, an open source computer program for creating particle-based stochastic, and spatially realistic models of cellular microchemistry. The model integrates four earlier detailed models of separate aspects of regulation of spine Ca2+ and CaMKII activity, each of which incorporate experimentally measured biochemical parameters and have been validated against experimental data. We validate the composite model by showing that it accurately predicts previous experimental measurements of effects of NMDA receptor activation, including high sensitivity of induction of LTP to phosphatase activity in vivo, and persistence of autophosphorylation for a period of minutes after the end of synaptic stimulation. We then use the model to probe aspects of the mechanism of regulation of autophosphorylation of CaMKII that are difficult to measure in vivo. We examine the effects of "CaM-trapping," a process in which the affinity for Ca2+/CaM increases several hundred-fold after autophosphorylation. We find that CaM-trapping does not increase the proportion of autophosphorylated subunits in holoenzymes after a complex stimulus, as previously hypothesized. Instead, CaM-trapping may dramatically prolong the lifetime of autophosphorylated CaMKII through steric hindrance of dephosphorylation by protein phosphatase 1. The results provide motivation for experimental measurement of the extent of suppression of dephosphorylation of CaMKII by bound Ca2+/CaM. The composite MCell4 model of biochemical effects of complex stimuli in synaptic spines is a powerful new tool for realistic, detailed dissection of mechanisms of synaptic plasticity.