Computational Analysis of Binding Interactions between the Ryanodine Receptor Type 2 and Calmodulin

Structural instability of ryanodine receptor 2 causes endoplasmic reticulum (ER) dysfunction as well as sarcoplasmic reticulum (SR) dysfunction

The type 2 ryanodine receptor (RyR2) is a giant Ca2+ (Ca)-releasing channel on the sarcoplasmic reticulum (SR) membrane, with subunits composed of 5000 amino acids constituting a homotetrameric channel. The N-terminal (1-220) and central (2300-2500) domain interactions (inter-subunit zipping interfaces) within RyR2 are located in close proximity to each other between different neighboring subunits and play an important "cornerstone" role in maintaining the tetrameric structure of RyR2. External stress such as oxidative stress causes Ca leak by destabilizing RyR2 (instability of the tetrameric structure) due to domain unzipping between N-terminal (1-220) and central (2300-2500) domains, followed by dissociation of calmodulin (CaM: binds to the RyR2 and stabilize RyR2) from RyR2. Ca leak from SR causes arrhythmias and myocardial dysfunction. RyR2 is also present in the endoplasmic reticulum (ER), thus it is not surprising that undesired Ca release from RyR2 on the ER is closely associated with various diseases involving ER dysfunction such as neurodegenerative diseases, diabetes, metabolic dysfunction-associated steatotic liver disease, chronic kidney disease, and autoimmune diseases. Pharmacological or genetic (point mutations within RyR2 that increase CaM-RyR2 affinity: knock-in RyR2-V3599K) RyR2 structural stabilization has shown potential therapeutic effects not only for SR failure-related diseases (malignant hyperthermia, arrhythmia, and heart failure) but also for ER failure-related disease. RyR2-stabilizers may function as a panacea for aging-related diseases.

Ca2+/calmodulin signaling in organismal aging and cellular senescence: Impact on human diseases

Molecular mechanisms of aging processes at the level of organisms and cells are in the focus of a large number of research laboratories. This research culminated in recent breakthroughs, which contributed to the better understanding of the natural aging process and aging associated malfunctions leading to age-related diseases. Ca2+ in connection with its master intracellular sensor protein calmodulin (CaM) regulates a plethora of crucial cellular processes orchestrating a wide range of signaling processes. This review focuses on the involvement of Ca2+/CaM in cellular mechanisms, which are associated with normal aging, as well as playing a role in the development of diseases connected with signaling processes during aging. We specifically highlight processes that involve inactivation of proteins, which take part in Ca2+/CaM regulatory systems by oxygen or nitrogen free radical species, during organismal aging and cellular senescence. As examples of organs where aging processes have recently been investigated, we chose to review the literature on molecular aging processes with involvement of Ca2+/CaM in heart and neuronal diseases, as well as in cancer and metabolic diseases, all deeply affected by aging. In addition, this article focuses on cellular senescence, a mechanism that may contribute to aging processes and therefore has been proposed as a target to interfere with the progression of age-associated diseases.

Identifying Key Binding Interactions Between the Cardiac L-Type Calcium Channel and Calmodulin Using Molecular Dynamics Simulations

Defects in the binding of the calcium sensing protein calmodulin (CaM) to the L-type calcium channel (CaV1.2) or to the ryanodine receptor type 2 (RyR2) can lead to dangerous cardiac arrhythmias with distinct phenotypes, such as long-QT syndrome (LQTS) and catecholaminergic ventricular tachycardia (CPVT). Certain CaM mutations lead to LQTS while other mutations lead to CPVT, but the mechanisms by which a specific mutation can lead to each disease phenotype are not well-understood. In this study, we use long, 2 μs molecular dynamics simulations and a multitrajectory approach to identify the key binding interactions between the IQ domain of CaV1.2 and CaM. Five key interactions are found between CaV1.2 and CaM in the C-lobe, 1 in the central linker, and 2 in the N-lobe. In addition, while 5 key interactions appear between residues 120-149 in the C-lobe of CaM when it interacts with CaV1.2, only 1 key interaction is found within this region of CaM when it interacts with the RyR2. We show that this difference in the distribution of key interactions correlates with the known distribution of CaM mutations that lead to LQTS or CPVT. This correlation suggests that a disruption of key binding interactions is a plausible mechanism that can lead to these two different disease phenotypes.

Life-threatening arrhythmogenic CaM mutations disrupt CaM binding to a distinct RyR2 CaM-binding pocket

Calmodulin (CaM) modulates the activity of several proteins that play a key role in excitation-contraction coupling (ECC). In cardiac muscle, the major binding partner of CaM is the type-2 ryanodine receptor (RyR2) and altered CaM binding contributes to defects in sarcoplasmic reticulum (SR) calcium (Ca²⁺) release. Many genetic studies have reported a series of CaM missense mutations in patients with a history of severe arrhythmogenic cardiac disorders. In the present study, we generated four missense CaM mutants (CaM^N98I, CaM^D132E, CaM^D134H and CaM^Q136P) and we used a CaM-RyR2 co-immunoprecipitation and a [³H]ryanodine binding assay to directly compare the relative RyR2-binding of wild type and mutant CaM proteins and to investigate the functional effects of these CaM mutations on RyR2 activity. Furthermore, isothermal titration calorimetry (ITC) experiments were performed to investigate and compare the interactions of the wild-type and mutant CaM proteins with various synthetic peptides located in the well-established RyR2 CaM-binding region (3584-3602aa), as well as another CaM-binding region (4255-4271aa) of human RyR2. Our data revealed that all four CaM mutants displayed dramatically reduced RyR2 interaction and defective modulation of [³H]ryanodine binding to RyR2, regardless of LQTS or CPVT association. Moreover, our isothermal titration calorimetry ITC data suggest that RyR2 3584-3602aa and 4255-4271aa regions interact with significant affinity with wild-type CaM, in the presence and absence of Ca²⁺, two regions that might contribute to a putative intra-subunit CaM-binding pocket. In contrast, screening the interaction of the four arrhythmogenic CaM mutants with two synthetic peptides that correspond to these RyR2 regions, revealed disparate binding properties and signifying differential mechanisms that contribute to reduced RyR2 association.

Experiment and Simulation Reveal Residue Details for How Target Binding Tunes Calmodulin's Calcium-Binding Properties

We aim to elucidate the molecular mechanism of the reciprocal relation of calmodulin's (CaM) target binding and its affinity for calcium ions (Ca²⁺), which is central to decoding CaM-dependent Ca²⁺ signaling in a cell. We employed stopped-flow experiments and coarse-grained molecular simulations that learn the coordination chemistry of Ca²⁺ in CaM from first-principle calculations. The associative memories as part of the coarse-grained force fields built on known protein structures further influence CaM's selection of its polymorphic target peptides in the simulations. We modeled the peptides from the Ca²⁺/CaM-binding domain of Ca²⁺/CaM-dependent kinase II (CaMKII), CaMKIIp (293-310) and selected distinctive mutations at the N-terminus. Our stopped-flow experiments have shown that the CaM's affinity for Ca²⁺ in the bound complex of Ca²⁺/CaM/CaMKIIp decreased significantly when Ca²⁺/CaM bound to the mutant peptide (296-AAA-298) compared to that bound to the wild-type peptide (296-RRK-298). The coarse-grained molecular simulations revealed that the 296-AAA-298 mutant peptide destabilized the structures of Ca²⁺-binding loops at the C-domain of CaM (c-CaM) due to both loss of electrostatic interactions and differences in polymorphic structures. We have leveraged a powerful coarse-grained approach to advance a residue-level understanding of the reciprocal relation in CaM, that could not be possibly achieved by other computational approaches.

Molecular Dynamics Simulations of the Cardiac Ryanodine Receptor Type 2 (RyR2) Gating Mechanism

Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to fatal cardiac arrhythmias such as catecholaminergic polymorphic ventricular tachycardia (CPVT). While many CPVT mutations are associated with an increase in Ca²⁺ leak from the sarcoplasmic reticulum, the mechanistic details of RyR2 channel gating are not well understood, and this poses a barrier in the development of new pharmacological treatments. To address this, we explore the gating mechanism of the RyR2 using molecular dynamics (MD) simulations. We test the effect of changing the conformation of certain structural elements by constructing chimera RyR2 structures that are derived from the currently available closed and open cryo-electron microscopy (cryo-EM) structures, and we then use MD simulations to relax the system. Our key finding is that the position of the S4-S5 linker (S4S5L) on a single subunit can determine whether the channel as a whole is open or closed. Our analysis reveals that the position of the S4S5L is regulated by interactions with the U-motif on the same subunit and with the S6 helix on an adjacent subunit. We find that, in general, channel gating is crucially dependent on high percent occupancy interactions between adjacent subunits. We compare our interaction analysis to 49 CPVT1 mutations in the literature and find that 73% appear near a high percent occupancy interaction between adjacent subunits. This suggests that disruption of cooperative, high percent occupancy interactions between adjacent subunits is a primary cause of channel leak and CPVT in mutant RyR2 channels.