Biophysical attributes that affect CaMKII activation deduced with a novel spatial stochastic simulation approach

Ca2+ dysregulation in cardiac stromal cells sustains fibro-adipose remodeling in Arrhythmogenic Cardiomyopathy and can be modulated by flecainide

Background

Cardiac mesenchymal stromal cells (C-MSC) were recently shown to differentiate into adipocytes and myofibroblasts to promote the aberrant remodeling of cardiac tissue that characterizes arrhythmogenic cardiomyopathy (ACM). A calcium (Ca2+) signaling dysfunction, mainly demonstrated in mouse models, is recognized as a mechanism impacting arrhythmic risk in ACM cardiomyocytes. Whether similar mechanisms influence ACM C-MSC fate is still unknown.

Thus, we aim to ascertain whether intracellular Ca2+ oscillations and the Ca2+ toolkit are altered in human C-MSC obtained from ACM patients, and to assess their link with C-MSC-specific ACM phenotypes.

Methods and results

ACM C-MSC show enhanced spontaneous Ca2+ oscillations and concomitant increased Ca2+/Calmodulin dependent kinase II (CaMKII) activation compared to control cells. This is manly linked to a constitutive activation of Store-Operated Ca2+ Entry (SOCE), which leads to enhanced Ca2+ release from the endoplasmic reticulum through inositol-1,4,5-trisphosphate receptors. By targeting the Ca2+ handling machinery or CaMKII activity, we demonstrated a causative link between Ca2+ oscillations and fibro-adipogenic differentiation of ACM C-MSC. Genetic silencing of the desmosomal gene PKP2 mimics the remodelling of the Ca2+ signalling machinery occurring in ACM C-MSC. The anti-arrhythmic drug flecainide inhibits intracellular Ca2+ oscillations and fibro-adipogenic differentiation by selectively targeting SOCE.

Conclusions

Altogether, our results extend the knowledge of Ca2+ dysregulation in ACM to the stromal compartment, as an etiologic mechanism of C-MSC-related ACM phenotypes. A new mode of action of flecainide on a novel mechanistic target is unveiled against the fibro-adipose accumulation in ACM.

pSpatiocyte: a high-performance simulator for intracellular reaction-diffusion systems

Background

Studies using quantitative experimental methods have shown that intracellular spatial distribution of molecules plays a central role in many cellular systems. Spatially resolved computer simulations can integrate quantitative data from these experiments to construct physically accurate models of the systems. Although computationally expensive, microscopic resolution reaction-diffusion simulators, such as Spatiocyte can directly capture intracellular effects comprising diffusion-limited reactions and volume exclusion from crowded molecules by explicitly representing individual diffusing molecules in space. To alleviate the steep computational cost typically associated with the simulation of large or crowded intracellular compartments, we present a parallelized Spatiocyte method called pSpatiocyte.

Results

The new high-performance method employs unique parallelization schemes on hexagonal close-packed (HCP) lattice to efficiently exploit the resources of common workstations and large distributed memory parallel computers. We introduce a coordinate system for fast accesses to HCP lattice voxels, a parallelized event scheduler, a parallelized Gillespie’s direct-method for unimolecular reactions, and a parallelized event for diffusion and bimolecular reaction processes. We verified the correctness of pSpatiocyte reaction and diffusion processes by comparison to theory. To evaluate the performance of pSpatiocyte, we performed a series of parallelized diffusion runs on the RIKEN K computer. In the case of fine lattice discretization with low voxel occupancy, pSpatiocyte exhibited 74% parallel efficiency and achieved a speedup of 7686 times with 663552 cores compared to the runtime with 64 cores. In the weak scaling performance, pSpatiocyte obtained efficiencies of at least 60% with up to 663552 cores. When executing the Michaelis-Menten benchmark model on an eight-core workstation, pSpatiocyte required 45- and 55-fold shorter runtimes than Smoldyn and the parallel version of ReaDDy, respectively. As a high-performance application example, we study the dual phosphorylation-dephosphorylation cycle of the MAPK system, a typical reaction network motif in cell signaling pathways.

Conclusions

pSpatiocyte demonstrates good accuracies, fast runtimes and a significant performance advantage over well-known microscopic particle methods in large-scale simulations of intracellular reaction-diffusion systems. The source code of pSpatiocyte is available at https://spatiocyte.org.

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

Ca2+/calmodulin-dependent protein kinase II (CaMKII) accounts for up to 2 percent of all brain protein and is essential to memory function. CaMKII activity is known to regulate dynamic shifts in the size and signaling strength of neuronal connections, a process known as synaptic plasticity. Increasingly, computational models are used to explore synaptic plasticity and the mechanisms regulating CaMKII activity. Conventional modeling approaches may exclude biophysical detail due to the impractical number of state combinations that arise when explicitly monitoring the conformational changes, ligand binding, and phosphorylation events that occur on each of the CaMKII holoenzyme's subunits. To manage the combinatorial explosion without necessitating bias or loss in biological accuracy, we use a specialized syntax in the software MCell to create a rule-based model of a twelve-subunit CaMKII holoenzyme. Here we validate the rule-based model against previous experimental measures of CaMKII activity and investigate molecular mechanisms of CaMKII regulation. Specifically, we explore how Ca2+/CaM-binding may both stabilize CaMKII subunit activation and regulate maintenance of CaMKII autophosphorylation. Noting that Ca2+/CaM and protein phosphatases bind CaMKII at nearby or overlapping sites, we compare model scenarios in which Ca2+/CaM and protein phosphatase do or do not structurally exclude each other's binding to CaMKII. Our results suggest a functional mechanism for the so-called "CaM trapping" phenomenon, wherein Ca2+/CaM may structurally exclude phosphatase binding and thereby prolong CaMKII autophosphorylation. We conclude that structural protection of autophosphorylated CaMKII by Ca2+/CaM may be an important mechanism for regulation of synaptic plasticity.

From membrane receptors to protein synthesis and actin cytoskeleton: Mechanisms underlying long lasting forms of synaptic plasticity

Synaptic plasticity, the activity dependent change in synaptic strength, forms the molecular foundation of learning and memory. Synaptic plasticity includes structural changes, with spines changing their size to accomodate insertion and removal of postynaptic receptors, which are correlated with functional changes. Of particular relevance for memory storage are the long lasting forms of synaptic plasticity which are protein synthesis dependent. Due to the importance of spine structural plasticity and protein synthesis, this review focuses on the signaling pathways that connect synaptic stimulation with regulation of protein synthesis and remodeling of the actin cytoskeleton. We also review computational models that implement novel aspects of molecular signaling in synaptic plasticity, such as the role of neuromodulators and spatial microdomains, as well as highlight the need for computational models that connect activation of memory kinases with spine actin dynamics.

Ca2+ Requirements for Long-Term Depression Are Frequency Sensitive in Purkinje Cells

Cerebellar long-term depression (LTD) is a form of synaptic plasticity dependent on postsynaptic Ca²⁺ changes. One fundamental question is how LTD is selectively induced by specific numbers of Ca²⁺ pulses and which are the frequency and duration of this train of pulses required for LTD induction. The molecular mechanism which leads the integration of postsynaptic Ca²⁺ pulses in the LTD signaling network has not been elucidated either. Recent publications have shown that Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) is required for LTD induction. Additionally, protein kinase C (PKC), CaMKII, and MAPK play an important role to transduce the frequency of Ca²⁺ pulses into their enzymatic activity levels; however, it is still unknown which enzymes are involved in decoding Ca²⁺ pulses in LTD. We have extended a stochastic model of LTD by adding the molecular network regulating CaMKII activity and its activation. We solved this model with stochastic engine for pathway simulation to include the effect of biochemical noise in LTD. We systematically investigated the dependence of LTD induction on stimulus frequencies, and we found that LTD is selectively induced by a specific number of Ca²⁺ spikes at different frequencies. We observed that CaMKII is essential to induce LTD, and LTD is only weakly induced when its Thr286 phosphorylation site has been deleted. We found that CaMKII decodes the frequency of Ca²⁺ spikes into different amounts of kinase activity during LTD induction. In addition, PKC and ERK enzyme activity is highly sensitive to the frequency and the number of Ca²⁺ pulses and this sensitivity has an important effect on LTD activation. This research predicts the postsynaptic Ca²⁺ requirements to induce LTD using a typical synaptic activation sequence and explains how LTD is selectively induced by specific number of Ca²⁺ pulses at different frequencies.