A multi-state model of the CaMKII dodecamer suggests a role for calmodulin in maintenance of autophosphorylation

A spatial model of autophosphorylation of Ca 2+ /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 in excitatory glutamatergic synapses in the Central Nervous System (CNS). It is initiated by a bout of synaptic activation that is strong enough to contribute to production of an action potential in the axon of the postsynaptic neuron, and it results in an increase in the size of postsynaptic depolarization during subsequent activity. The first step leading to LTP is activation and autophosphorylation of an abundant postsynaptic enzyme, Ca 2+ /calmodulin-dependent protein kinase II (CaMKII). We use simulation of activation of CaMKII holoenzymes in a realistic spatial model of a spine synapse, created in MCell4, to test three hypotheses about how the autophosphorylation response of CaMKII is shaped during a repeated high-frequency stimulus. First, the simulation results indicate that autophosphorylation of CaMKII does not constitute a bistable switch under biologically realistic conditions. Instead, prolonged autophosphorylation of CaMKII may contribute to a biochemical "kinetic proof-reading" mechanism that controls induction of synaptic plasticity. Second, concentration of CaMKII near the postsynaptic membrane increases the local concentration of kinase activity. However, neither localization nor "Ca 2+ -calmodulin-trapping (CaM-trapping)" increase the proportion of autophosphorylated subunits in holoenzymes after a complex stimulus, as previously hypothesized. Finally, we show that, as hypothesized, the amplitude of autophosphorylation in the first 30 seconds after a stimulus is extremely sensitive to the level and location of PP1 activity when PP1 is present in biologically accurate amounts. We further show that prolonged steric hindrance of dephosphorylation of CaMKII, caused by CaM-trapping, can increase the amplitude of autophosphorylation after a complex stimulus. These simulation results sharpen our quantitative understanding of the early events leading to LTP at excitatory synapses.

Author Summary

Neurons in the brain are interconnected in an organized fashion by synapses that transmit neuronal activity from one neuron to another. Most of the billions of neurons in the brain have about 10,000 synapses spread over the neuronal membrane. Information is stored in the brain when the ability of specific synapses to pass along neuronal activity is strengthened resulting in formation of new networks. The increase in strength of a synapse is tightly controlled by the frequency and amplitude of its activity, and by neurohormonal signals, which, in combination, can cause long-lasting biochemical changes at the synapse that underlie learning and memory. Defects in these biochemical pathways cause mental and neurological diseases. To develop treatments, we need to understand the precise choreography of these critical biochemical changes. However, the tiny size of the synaptic compartment makes precise measurements of the biochemical reactions impossible. We have used computer simulation techniques and information gathered from experiments on purified synaptic proteins to simulate, within a single synapse, the choreography of the first biochemical step in synaptic strengthening: activation of the enzyme Ca 2+ / calmodulin-dependent protein kinase II. Our results provide insights that can be used in future studies to develop treatments for neuronal diseases.

Biophysical Modeling of Synaptic Plasticity

Dendritic spines are small, bulbous compartments that function as postsynaptic sites and undergo intense biochemical and biophysical activity. The role of the myriad signaling pathways that are implicated in synaptic plasticity is well studied. A recent abundance of quantitative experimental data has made the events associated with synaptic plasticity amenable to quantitative biophysical modeling. Spines are also fascinating biophysical computational units because spine geometry, signal transduction, and mechanics work in a complex feedback loop to tune synaptic plasticity. In this sense, ideas from modeling cell motility can inspire us to develop multiscale approaches for predictive modeling of synaptic plasticity. In this article, we review the key steps in postsynaptic plasticity with a specific focus on the impact of spine geometry on signaling, cytoskeleton rearrangement, and membrane mechanics. We summarize the main experimental observations and highlight how theory and computation can aid our understanding of these complex processes.

Real-time single-molecule imaging of CaMKII-calmodulin interactions

The binding of calcium/calmodulin (CAM) to calcium/calmodulin-dependent protein kinase II (CaMKII) initiates an ATP-driven cascade that triggers CaMKII autophosphorylation. The autophosphorylation in turn increases the CaMKII affinity for CAM. Here, we studied the ATP dependence of CAM association with the actin-binding CaMKIIβ isoform using single-molecule total internal reflection fluorescence microscopy. Rhodamine-CAM associations/dissociations to surface-immobilized Venus-CaMKIIβ were resolved with 0.5 s resolution from video records, batch-processed with a custom algorithm. CAM occupancy was determined simultaneously with spot-photobleaching measurement of CaMKII holoenzyme stoichiometry. We show the ATP-dependent increase of the CAM association requires dimer formation for both the α and β isoforms. The study of mutant β holoenzymes revealed that the ATP-dependent increase in CAM affinity results in two distinct states. The phosphorylation-defective (T287.306-307A) holoenzyme resides only in the low-affinity state. CAM association is further reduced in the T287A holoenzyme relative to T287.306-307A. In the absence of ATP, the affinity of CAM for the T287.306-307A mutant and the wild-type monomer are comparable. The affinity of the ATP-binding impaired (K43R) mutant is even weaker. In ATP, the K43R holoenzyme resides in the low-affinity state. The phosphomimetic mutant (T287D) resides only in a 1000-fold higher-affinity state, with mean CAM occupancy of more than half of the 14-mer holoenzyme stoichiometry in picomolar CAM. ATP promotes T287D holoenzyme disassembly but does not elevate CAM occupancy. Single Poisson distributions characterized the ATP-dependent CAM occupancy of mutant holoenzymes. In contrast, the CAM occupancy of the wild-type population had a two-state distribution with both low- and high-affinity states represented. The low-affinity state was the dominant state, a result different from published in vitro assays. Differences in assay conditions can alter the balance between activating and inhibitory autophosphorylation. Bound ATP could be sufficient for CaMKII structural function, while antagonistic autophosphorylations may tune CaMKII kinase-regulated action-potential frequency decoding in vivo.

Vitamin D, Calbindin, and calcium signaling: Unraveling the Alzheimer's connection

Calcium is a ubiquitous second messenger that is indispensable in regulating neurotransmission and memory formation. A precise intracellular calcium level is achieved through the concerted action of calcium channels, and calcium exerts its effect by binding to an array of calcium-binding proteins, including calmodulin (CAM), calcium-calmodulin complex-dependent protein kinase-II (CAMK-II), calbindin (CAL), and calcineurin (CAN). Calbindin orchestrates a plethora of signaling events that regulate synaptic transmission and depolarizing signals. Vitamin D, an endogenous fat-soluble metabolite, is synthesized in the skin upon exposure to ultraviolet B radiation. It modulates calcium signaling by increasing the expression of the calcium-sensing receptor (CaSR), stimulating phospholipase C activity, and regulating the expression of calcium channels such as TRPV6. Vitamin D also modulates the activity of calcium-binding proteins, including CAM and calbindin, and increases their expression. Calbindin, a high-affinity calcium-binding protein, is involved in calcium buffering and transport in neurons. It has been shown to inhibit apoptosis and caspase-3 activity stimulated by presenilin 1 and 2 in AD. Whereas CAM, another calcium-binding protein, is implicated in regulating neurotransmitter release and memory formation by phosphorylating CAN, CAMK-II, and other calcium-regulated proteins. CAMK-II and CAN regulate actin-induced spine shape changes, which are further modulated by CAM. Low levels of both calbindin and vitamin D are attributed to the pathology of Alzheimer's disease. Further research on vitamin D via calbindin-CAMK-II signaling may provide newer insights, revealing novel therapeutic targets and strategies for treatment.

A biochemical description of postsynaptic plasticity-with timescales ranging from milliseconds to seconds

Synaptic plasticity [long-term potentiation/depression (LTP/D)], is a cellular mechanism underlying learning. Two distinct types of early LTP/D (E-LTP/D), acting on very different time scales, have been observed experimentally-spike timing dependent plasticity (STDP), on time scales of tens of ms; and behavioral time scale synaptic plasticity (BTSP), on time scales of seconds. BTSP is a candidate for a mechanism underlying rapid learning of spatial location by place cells. Here, a computational model of the induction of E-LTP/D at a spine head of a synapse of a hippocampal pyramidal neuron is developed. The single-compartment model represents two interacting biochemical pathways for the activation (phosphorylation) of the kinase (CaMKII) with a phosphatase, with ion inflow through channels (NMDAR, CaV1,Na). The biochemical reactions are represented by a deterministic system of differential equations, with a detailed description of the activation of CaMKII that includes the opening of the compact state of CaMKII. This single model captures realistic responses (temporal profiles with the differing timescales) of STDP and BTSP and their asymmetries. The simulations distinguish several mechanisms underlying STDP vs. BTSP, including i) the flow of [Formula: see text] through NMDAR vs. CaV1 channels, and ii) the origin of several time scales in the activation of CaMKII. The model also realizes a priming mechanism for E-LTP that is induced by [Formula: see text] flow through CaV1.3 channels. Once in the spine head, this small additional [Formula: see text] opens the compact state of CaMKII, placing CaMKII ready for subsequent induction of LTP.

Synaptic memory and CaMKII

Ca2+/calmodulin-dependent protein kinase II (CaMKII) and long-term potentiation (LTP) were discovered within a decade of each other and have been inextricably intertwined ever since. However, like many marriages, it has had its up and downs. Based on the unique biochemical properties of CaMKII, it was proposed as a memory molecule before any physiological linkage was made to LTP. However, as reviewed here, the convincing linkage of CaMKII to synaptic physiology and behavior took many decades. New technologies were critical in this journey, including in vitro brain slices, mouse genetics, single-cell molecular genetics, pharmacological reagents, protein structure, and two-photon microscopy, as were new investigators attracted by the exciting challenge. This review tracks this journey and assesses the state of this marriage 40 years on. The collective literature impels us to propose a relatively simple model for synaptic memory involving the following steps that drive the process: 1) Ca2+ entry through N-methyl-d-aspartate (NMDA) receptors activates CaMKII. 2) CaMKII undergoes autophosphorylation resulting in constitutive, Ca2+-independent activity and exposure of a binding site for the NMDA receptor subunit GluN2B. 3) Active CaMKII translocates to the postsynaptic density (PSD) and binds to the cytoplasmic C-tail of GluN2B. 4) The CaMKII-GluN2B complex initiates a structural rearrangement of the PSD that may involve liquid-liquid phase separation. 5) This rearrangement involves the PSD-95 scaffolding protein, α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs), and their transmembrane AMPAR-regulatory protein (TARP) auxiliary subunits, resulting in an accumulation of AMPARs in the PSD that underlies synaptic potentiation. 6) The stability of the modified PSD is maintained by the stability of the CaMKII-GluN2B complex. 7) By a process of subunit exchange or interholoenzyme phosphorylation CaMKII maintains synaptic potentiation in the face of CaMKII protein turnover. There are many other important proteins that participate in enlargement of the synaptic spine or modulation of the steps that drive and maintain the potentiation. In this review we critically discuss the data underlying each of the steps. As will become clear, some of these steps are more firmly grounded than others, and we provide suggestions as to how the evidence supporting these steps can be strengthened or, based on the new data, be replaced. Although the journey has been a long one, the prospect of having a detailed cellular and molecular understanding of learning and memory is at hand.

CaMKII: a central molecular organizer of synaptic plasticity, learning and memory

Calcium-calmodulin (CaM)-dependent protein kinase II (CaMKII) is the most abundant protein in excitatory synapses and is central to synaptic plasticity, learning and memory. It is activated by intracellular increases in calcium ion levels and triggers molecular processes necessary for synaptic plasticity. CaMKII phosphorylates numerous synaptic proteins, thereby regulating their structure and functions. This leads to molecular events crucial for synaptic plasticity, such as receptor trafficking, localization and activity; actin cytoskeletal dynamics; translation; and even transcription through synapse-nucleus shuttling. Several new tools affording increasingly greater spatiotemporal resolution have revealed the link between CaMKII activity and downstream signalling processes in dendritic spines during synaptic and behavioural plasticity. These technologies have provided insights into the function of CaMKII in learning and memory.

Conformational spread drives the evolution of the calcium-calmodulin protein kinase II

The calcium calmodulin (Ca²⁺/CaM) dependent protein kinase II (CaMKII) decodes Ca²⁺ frequency oscillations. The CaMKIIα isoform is predominantly expressed in the brain and has a central role in learning. I matched residue and organismal evolution with collective motions deduced from the atomic structure of the human CaMKIIα holoenzyme to learn how its ring architecture abets function. Protein dynamic simulations showed its peripheral kinase domains (KDs) are conformationally coupled via lateral spread along the central hub. The underlying β-sheet motions in the hub or association domain (AD) were deconvolved into dynamic couplings based on mutual information. They mapped onto a coevolved residue network to partition the AD into two distinct sectors. A second, energetically stressed sector was added to ancient bacterial enzyme dimers for assembly of the ringed hub. The continued evolution of the holoenzyme after AD–KD fusion targeted the sector’s ring contacts coupled to the KD. Among isoforms, the α isoform emerged last and, it alone, mutated rapidly after the poikilotherm–homeotherm jump to match the evolution of memory. The correlation between dynamics and evolution of the CaMKII AD argues single residue substitutions fine-tune hub conformational spread. The fine-tuning could increase CaMKIIα Ca²⁺ frequency response range for complex learning functions.

Data-driven learning how oncogenic gene expression locally alters heterocellular networks

Developing drugs increasingly relies on mechanistic modeling and simulation. Models that capture causal relations among genetic drivers of oncogenesis, functional plasticity, and host immunity complement wet experiments. Unfortunately, formulating such mechanistic cell-level models currently relies on hand curation, which can bias how data is interpreted or the priority of drug targets. In modeling molecular-level networks, rules and algorithms are employed to limit a priori biases in formulating mechanistic models. Here we combine digital cytometry with Bayesian network inference to generate causal models of cell-level networks linking an increase in gene expression associated with oncogenesis with alterations in stromal and immune cell subsets from bulk transcriptomic datasets. We predict how increased Cell Communication Network factor 4, a secreted matricellular protein, alters the tumor microenvironment using data from patients diagnosed with breast cancer and melanoma. Predictions are then tested using two immunocompetent mouse models for melanoma, which provide consistent experimental results.

While mechanistic models play increasing roles in immuno-oncology, hand network curation is current practice. Here the authors use a Bayesian data-driven approach to infer how expression of a secreted oncogene alters the cellular landscape within the tumor.

Validation Through Collaboration: Encouraging Team Efforts to Ensure Internal and External Validity of Computational Models of Biochemical Pathways

Computational modelling of biochemical reaction pathways is an increasingly important part of neuroscience research. In order to be useful, computational models need to be valid in two senses: First, they need to be consistent with experimental data and able to make testable predictions (external validity). Second, they need to be internally consistent and independently reproducible (internal validity). Here, we discuss both types of validity and provide a brief overview of tools and technologies used to ensure they are met. We also suggest the introduction of new collaborative technologies to ensure model validity: an incentivised experimental database for external validity and reproducibility audits for internal validity. Both rely on FAIR principles and on collaborative science practices.

Physical interactions driving the activation/inhibition of calcium/calmodulin dependent protein kinase II

CaMKII is a protein kinase whose function is regulated by the binding of the Calcium/Calmodulin complex (Ca²⁺/CaM). It is a major player in the Long Term Potentiation process where it acts as a molecular switch, oscillating between inhibited and active conformations. The mechanism for the switching is thought to be initiated by Ca²⁺/CaM binding, which allows the trans-phosphorylation of a subunit of CaMKII by a neighboring kinase, leading to the active state of the system. A combination of all-atom and coarse-grained MD simulations with free energy calculations, led us to reveal an interplay of electrostatic forces exerted by Ca²⁺/CaM on CaMKII, which initiate the activation process. The highly electrically charged Ca²⁺/CaM neutralizes basic regions in the linker domain of CaMKII, facilitating its opening and consequent activation. The emerging picture of CaMKII's behavior highlights the preponderance of electrostatic interactions, which are modulated by the presence of Ca²⁺/CaM and the phosphorylation of key sites.

Interactions between calmodulin and neurogranin govern the dynamics of CaMKII as a leaky integrator

Calmodulin-dependent kinase II (CaMKII) has long been known to play an important role in learning and memory as well as long term potentiation (LTP). More recently it has been suggested that it might be involved in the time averaging of synaptic signals, which can then lead to the high precision of information stored at a single synapse. However, the role of the scaffolding molecule, neurogranin (Ng), in governing the dynamics of CaMKII is not yet fully understood. In this work, we adopt a rule-based modeling approach through the Monte Carlo method to study the effect of Ca2+ signals on the dynamics of CaMKII phosphorylation in the postsynaptic density (PSD). Calcium surges are observed in synaptic spines during an EPSP and back-propagating action potential due to the opening of NMDA receptors and voltage dependent calcium channels. Using agent-based models, we computationally investigate the dynamics of phosphorylation of CaMKII monomers and dodecameric holoenzymes. The scaffolding molecule, Ng, when present in significant concentration, limits the availability of free calmodulin (CaM), the protein which activates CaMKII in the presence of calcium. We show that Ng plays an important modulatory role in CaMKII phosphorylation following a surge of high calcium concentration. We find a non-intuitive dependence of this effect on CaM concentration that results from the different affinities of CaM for CaMKII depending on the number of calcium ions bound to the former. It has been shown previously that in the absence of phosphatase, CaMKII monomers integrate over Ca2+ signals of certain frequencies through autophosphorylation (Pepke et al, Plos Comp. Bio., 2010). We also study the effect of multiple calcium spikes on CaMKII holoenzyme autophosphorylation, and show that in the presence of phosphatase, CaMKII behaves as a leaky integrator of calcium signals, a result that has been recently observed in vivo. Our models predict that the parameters of this leaky integrator are finely tuned through the interactions of Ng, CaM, CaMKII, and PP1, providing a mechanism to precisely control the sensitivity of synapses to calcium signals. Author Summary not valid for PLOS ONE submissions.