Biophysical attributes that affect CaMKII activation deduced with a novel spatial stochastic simulation approach
Authors: Ximing Li, William R Holmes
Journal: PLOS Computational Biology
Accepted: 2017-12-19T00:00:00Z


Calcium/calmodulin-dependent protein kinase II (CaMKII) holoenzymes play a critical role in decoding Ca2+ signals in neurons. Understanding how this occurs has been the focus of numerous studies including many that use models. However, CaMKII is notoriously difficult to simulate in detail because of its multi-subunit nature, which causes a combinatorial explosion in the number of species that must be modeled. To study the Ca2+-calmodulin-CaMKII reaction network with detailed kinetics while including the effect of diffusion, we have customized an existing stochastic particle-based simulator, Smoldyn, to manage the problem of combinatorial explosion. With this new method, spatial and temporal aspects of the signaling network can be studied without compromising biochemical details. We used this new method to examine how calmodulin molecules, both partially loaded and fully loaded with Ca2+, choose pathways to interact with and activate CaMKII under various Ca2+ input conditions. We found that the dependence of CaMKII phosphorylation on Ca2+ signal frequency is intrinsic to the network kinetics and the activation pattern can be modulated by the relative amount of Ca2+ to calmodulin and by the rate of Ca2+ diffusion. Depending on whether Ca2+ influx is saturating or not, calmodulin molecules could choose different routes within the network to activate CaMKII subunits, resulting in different frequency dependence patterns. In addition, the size of the holoenzyme produces a subtle effect on CaMKII activation. The more extended the subunits are organized, the easier for calmodulin molecules to access and activate the subunits. The findings suggest that particular intracellular environmental factors such as crowding and calmodulin availability can play an important role in decoding Ca2+ signals and can give rise to distinct CaMKII activation patterns in dendritic spines, Ca2+ channel nanodomains and cytoplasm.
Author summary: Ca2+ signals are commonly used by cells for various types of activities. In neurons, Ca2+ can regulate gene expression and dendritic spine enlargement, strengthen synaptic connectivity, and promote neural growth or even death. One important Ca2+ binding protein is calmodulin, which has a wide range of downstream effectors that relay and interpret Ca2+ signals. These downstream proteins have different calmodulin binding kinetics and can interact with calmodulin differently depending on Ca2+ signaling patterns and in particular, on signaling frequencies. In this study, we focus on one particular pathway that is thought to be important for learning and memory and involves calmodulin and the effector protein, Ca2+/calmodulin-dependent protein kinase II (CaMKII). We have developed a novel computational method that avoids the problem of combinatorial explosion inherent with modeling this pathway, and the method allows us to study the pathway in more detail than conventional methods. Experiments have shown that CaMKII activity is sensitive to Ca2+ signal frequency and our models demonstrate how this frequency dependence relies on the amount of Ca2+ input, calmodulin availability and the Ca2+ diffusion rate.