Calmodulin mutations associated with long QT syndrome prevent inactivation of cardiac L-type Ca(2+) currents and promote proarrhythmic behavior in ventricular myocytes

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.

Calcium- and calmodulin-dependent inhibition of NMDA receptor currents

Calcium ions (Ca2+) reduce NMDA receptor currents through several distinct mechanisms. Among these, calmodulin (CaM)-dependent inhibition (CDI) accomplishes rapid, reversible, and incomplete reduction of the NMDA receptor currents in response to elevations in intracellular Ca2+. Quantitative and mechanistic descriptions of CDI of NMDA receptor-mediated signals have been marred by variability originating, in part, from differences in the conditions and metrics used to evaluate this process across laboratories. Recent ratiometric approaches to measure the magnitude and kinetics of NMDA receptor CDI have facilitated rapid insights into this phenomenon. Notably, the kinetics and magnitude of NMDA receptor CDI depend on the degree of saturation of its CaM binding sites, which represent the bona fide calcium sensor for this type of inhibition, the kinetics and magnitude of the Ca2+ signal, which depends on the biophysical properties of the NMDA receptor or of adjacent Ca2+ sources, and on the relative distribution of Ca2+ sources and CaM molecules. Given that all these factors vary widely during development, across cell types, and with physiological and pathological states, it is important to understand how NMDA receptor CDI develops and how it contributes to signaling in the central nervous system. Here, we review briefly these recent advances and highlight remaining questions about the structural and kinetic mechanisms of NMDA receptor CDI. Given that pathologies can arise from several sources, including mutations in the NMDA receptor and in CaM, understanding how CaM responds to intracellular Ca2+ signals to initiate conformational changes in NMDA receptors, and mapping the structural domains responsible will help to envision novel therapeutic strategies to neuropsychiatric diseases, which presently have limited available treatments.

Asymmetric contribution of a selectivity filter gate in triggering inactivation of CaV1.3 channels

Voltage-dependent and Ca2+-dependent inactivation (VDI and CDI, respectively) of CaV channels are two biologically consequential feedback mechanisms that fine-tune Ca2+ entry into neurons and cardiomyocytes. Although known to be initiated by distinct molecular events, how these processes obstruct conduction through the channel pore remains poorly defined. Here, focusing on ultrahighly conserved tryptophan residues in the interdomain interfaces near the selectivity filter of CaV1.3, we demonstrate a critical role for asymmetric conformational changes in mediating VDI and CDI. Specifically, mutagenesis of the domain III-IV interface, but not others, enhanced VDI. Molecular dynamics simulations demonstrate that mutations in distinct selectivity filter interfaces differentially impact conformational flexibility. Furthermore, mutations in distinct domains preferentially disrupt CDI mediated by the N- versus C-lobes of CaM, thus uncovering a scheme of structural bifurcation of CaM signaling. These findings highlight the fundamental importance of the asymmetric arrangement of the pseudotetrameric CaV pore domain for feedback inhibition.

Calmodulin Mutations in Human Disease

Calcium ions (Ca²⁺) are the basis of a unique and potent array of cellular responses. Calmodulin (CaM) is a small but vital protein that is able to rapidly transmit information about changes in Ca²⁺ concentrations to its regulatory targets. CaM plays a critical role in cellular Ca²⁺ signaling, and interacts with a myriad of target proteins. Ca²⁺-dependent modulation by CaM is a major component of a diverse array of processes, ranging from gene expression in neurons to the shaping of the cardiac action potential in heart cells. Furthermore, the protein sequence of CaM is highly evolutionarily conserved, and identical CaM proteins are encoded by three independent genes (CALM1-3) in humans. Mutations within any of these three genes may lead to severe cardiac deficits including severe long QT syndrome (LQTS) and/or catecholaminergic polymorphic ventricular tachycardia (CPVT). Research into disease-associated CaM variants has identified several proteins modulated by CaM that are likely to underlie the pathogenesis of these calmodulinopathies, including the cardiac L-type Ca²⁺ channel (LTCC) Ca_V1.2, and the sarcoplasmic reticulum Ca²⁺ release channel, ryanodine receptor 2 (RyR2). Here, we review the research that has been done to identify calmodulinopathic CaM mutations and evaluate the mechanisms underlying their role in disease.

Arrhythmia-associated calmodulin variants interact with KCNQ1 to confer aberrant membrane trafficking and function

Missense variants in calmodulin (CaM) predispose patients to arrhythmias associated with high mortality rates ("calmodulinopathy"). As CaM regulates many key cardiac ion channels, an understanding of disease mechanism associated with CaM variant arrhythmias requires elucidating individual CaM variant effects on distinct channels. One key CaM regulatory target is the KCNQ1 (KV7.1) voltage-gated potassium channel that carries the IKs current. Yet, relatively little is known as to how CaM variants interact with KCNQ1 or affect its function. Here, we take a multipronged approach employing a live-cell fluorescence resonance energy transfer binding assay, fluorescence trafficking assay, and functional electrophysiology to characterize >10 arrhythmia-associated CaM variants for effect on KCNQ1 CaM binding, membrane trafficking, and channel function. We identify one variant (G114W) that exhibits severely weakened binding to KCNQ1 but find that most other CaM variants interact with similar binding affinity to KCNQ1 when compared with CaM wild-type over physiological Ca2+ ranges. We further identify several CaM variants that affect KCNQ1 and IKs membrane trafficking and/or baseline current activation kinetics, thereby delineating KCNQ1 dysfunction in calmodulinopathy. Lastly, we identify CaM variants with no effect on KCNQ1 function. This study provides extensive functional data that reveal how CaM variants contribute to creating a proarrhythmic substrate by causing abnormal KCNQ1 membrane trafficking and current conduction. We find that CaM variant regulation of KCNQ1 is not uniform with effects varying from benign to significant loss of function, suggesting how CaM variants predispose patients to arrhythmia via the dysregulation of multiple cardiac ion channels. Classification: Biological, Health, and Medical Sciences, Physiology.

Arrhythmia-Associated Calmodulin E105A Mutation Alters the Binding Affinity of CaM to a Ryanodine Receptor 2 CaM-Binding Pocket

Calmodulin (CaM) is a small, multifunctional calcium (Ca2+)-binding sensor that binds and regulates the open probability of cardiac ryanodine receptor 2 (RyR2) at both low and high cytosolic Ca2+ concentrations. Recent isothermal titration calorimetry (ITC) studies of a number of peptides that correspond to different regions of human RyR2 showed that two regions of human RyR2 (3584-3602aa and 4255-4271aa) bind with high affinity to CaM, suggesting that these two regions might contribute to a putative RyR2 intra-subunit CaM-binding pocket. Moreover, a previously characterized de novo long QT syndrome (LQTS)-associated missense CaM mutation (E105A) which was identified in a 6-year-old boy, who experienced an aborted first episode of cardiac arrest revealed that this mutation dysregulates normal cardiac function in zebrafish by a complex mechanism that involves alterations in both CaM-Ca2+ and CaM-RyR2 interactions. Herein, to gain further insight into how the CaM E105A mutation leads to severe cardiac arrhythmia, we generated large quantities of recombinant CaMWT and CaME105A proteins. We then performed ITC experiments to investigate and compare the interactions of CaMWT and CaME105A mutant protein with two synthetic peptides that correspond to the two aforementioned human RyR2 regions, which we have proposed to contribute to the RyR2 CaM-binding pocket. Our data reveal that the E105A mutation has a significant negative effect on the interaction of CaM with both RyR2 regions in the presence and absence of Ca2+, highlighting the potential contribution of these two human RyR2 regions to an RyR2 CaM-binding pocket, which may be essential for physiological CaM/RyR2 association and thus channel regulation.

Asymmetric Contribution of a Selectivity Filter Gate in Triggering Inactivation of CaV1.3 Channels

Voltage-dependent and Ca2+-dependent inactivation (VDI and CDI, respectively) of CaV channels are two biologically consequential feedback mechanisms that fine-tune Ca2+ entry into neurons and cardiomyocytes. Although known to be initiated by distinct molecular events, how these processes obstruct conduction through the channel pore remains poorly defined. Here, focusing on ultra-highly conserved tryptophan residues in the inter-domain interfaces near the selectivity filter of CaV1.3, we demonstrate a critical role for asymmetric conformational changes in mediating VDI and CDI. Specifically, mutagenesis of the domain III-IV interface, but not others, enhanced VDI. Molecular dynamics simulations demonstrate that mutations in distinct selectivity filter interfaces differentially impact conformational flexibility. Furthermore, mutations in distinct domains preferentially disrupt CDI mediated by the N- versus C-lobes of CaM, thus uncovering a scheme of structural bifurcation of CaM signaling. These findings highlight the fundamental importance of the asymmetric arrangement of the pseudo-tetrameric CaV pore domain for feedback inhibition.

Clinical presentation of calmodulin mutations: the International Calmodulinopathy Registry

AIMS

Calmodulinopathy due to mutations in any of the three CALM genes (CALM1-3) causes life-threatening arrhythmia syndromes, especially in young individuals. The International Calmodulinopathy Registry (ICalmR) aims to define and link the increasing complexity of the clinical presentation to the underlying molecular mechanisms.

METHODS AND RESULTS

The ICalmR is an international, collaborative, observational study, assembling and analysing clinical and genetic data on CALM-positive patients. The ICalmR has enrolled 140 subjects (median age 10.8 years [interquartile range 5-19]), 97 index cases and 43 family members. CALM-LQTS and CALM-CPVT are the prevalent phenotypes. Primary neurological manifestations, unrelated to post-anoxic sequelae, manifested in 20 patients. Calmodulinopathy remains associated with a high arrhythmic event rate (symptomatic patients, n = 103, 74%). However, compared with the original 2019 cohort, there was a reduced frequency and severity of all cardiac events (61% vs. 85%; P = .001) and sudden death (9% vs. 27%; P = .008). Data on therapy do not allow definitive recommendations. Cardiac structural abnormalities, either cardiomyopathy or congenital heart defects, are present in 30% of patients, mainly CALM-LQTS, and lethal cases of heart failure have occurred. The number of familial cases and of families with strikingly different phenotypes is increasing.

CONCLUSION

Calmodulinopathy has pleiotropic presentations, from channelopathy to syndromic forms. Clinical severity ranges from the early onset of life-threatening arrhythmias to the absence of symptoms, and the percentage of milder and familial forms is increasing. There are no hard data to guide therapy, and current management includes pharmacological and surgical antiadrenergic interventions with sodium channel blockers often accompanied by an implantable cardioverter-defibrillator.

Sex and Gene Influence Arrhythmia Susceptibility in Murine Models of Calmodulinopathy

BACKGROUND

Pathogenic variants in genes encoding CaM (calmodulin) are associated with a life-threatening ventricular arrhythmia syndrome (calmodulinopathy). The in vivo consequences of CaM variants have not been studied extensively and there is incomplete understanding of the genotype-phenotype relationship for recurrent variants. We investigated effects of different factors on calmodulinopathy phenotypes using 2 mouse models with a recurrent pathogenic variant (N98S) in Calm1 or Calm2.

METHODS

Genetically engineered mice with heterozygous N98S pathogenic variants in Calm1 or Calm2 were generated. Differences between the sexes and affected genes were assessed using multiple physiological assays at the cellular and whole animal levels. Statistical significance among groups was evaluated using 1-way ANOVA or the Kruskal-Wallis test when data were not normally distributed.

RESULTS

Calm1N98S/+ (Calm1S/+) or Calm2N98S/+ (Calm2S/+) mice exhibited sinus bradycardia and were more susceptible to arrhythmias after exposure to epinephrine and caffeine. Male Calm1S/+ mice had the most severe arrhythmia phenotype with evidence of early embryonic lethality, greater susceptibility for arrhythmic events, frequent premature beats, corrected QT prolongation, and more heart rate variability after epinephrine and caffeine than females with the same genotype. Calm2 S/+ mice exhibited a less severe phenotype, with female Calm2 S/+ mice having the least severe arrhythmia susceptibility. Flecainide was not effective in preventing arrhythmias in heterozygous CaM-N98S mice. Intracellular Ca2+ transients observed in isolated ventricular cardiomyocytes from male heterozygous CaM-N98S mice had lower peak amplitudes and slower sarcoplasmic reticulum Ca2+ release following in vitro exposure to epinephrine and caffeine, which were not observed in cardiomyocytes from heterozygous female CaM-N98S mice.

CONCLUSIONS

We report heterogeneity in arrhythmia susceptibility and cardiomyocyte Ca2+ dynamics among male and female mice heterozygous for a recurrent pathogenic variant in Calm1 or Calm2, illustrating a complex calmodulinopathy phenotype in vivo. Further investigation of sex and genetic differences may help identify the molecular basis for this heterogeneity.

Long QT syndrome-associated calmodulin variants disrupt the activity of the slowly activating delayed rectifier potassium channel

Calmodulin (CaM) is a highly conserved mediator of calcium (Ca2+ )-dependent signalling and modulates various cardiac ion channels. Genotyping has revealed several CaM mutations associated with long QT syndrome (LQTS). LQTS patients display prolonged ventricular recovery times (QT interval), increasing their risk of incurring life-threatening arrhythmic events. Loss-of-function mutations to Kv7.1 (which drives the slow delayed rectifier potassium current, IKs, a key ventricular repolarising current) are the largest contributor to congenital LQTS (>50% of cases). CaM modulates Kv7.1 to produce a Ca2+ -sensitive IKs, but little is known about the consequences of LQTS-associated CaM mutations on Kv7.1 function. Here, we present novel data characterising the biophysical and modulatory properties of three LQTS-associated CaM variants (D95V, N97I and D131H). We showed that mutations induced structural alterations in CaM and reduced affinity for Kv7.1, when compared with wild-type (WT). Using HEK293T cells expressing Kv7.1 channel subunits (KCNQ1/KCNE1) and patch-clamp electrophysiology, we demonstrated that LQTS-associated CaM variants reduced current density at systolic Ca2+ concentrations (1 μm), revealing a direct QT-prolonging modulatory effect. Our data highlight for the first time that LQTS-associated perturbations to CaM's structure impede complex formation with Kv7.1 and subsequently result in reduced IKs. This provides a novel mechanistic insight into how the perturbed structure-function relationship of CaM variants contributes to the LQTS phenotype. KEY POINTS: Calmodulin (CaM) is a ubiquitous, highly conserved calcium (Ca2+ ) sensor playing a key role in cardiac muscle contraction. Genotyping has revealed several CaM mutations associated with long QT syndrome (LQTS), a life-threatening cardiac arrhythmia syndrome. LQTS-associated CaM variants (D95V, N97I and D131H) induced structural alterations, altered binding to Kv7.1 and reduced IKs. Our data provide a novel mechanistic insight into how the perturbed structure-function relationship of CaM variants contributes to the LQTS phenotype.

Calmodulin mutations affecting Gly114 impair binding to the NaV1.5 IQ-domain

Missense variants in CALM genes encoding the Ca2+-binding protein calmodulin (CaM) cause severe cardiac arrhythmias. The disease mechanisms have been attributed to dysregulation of RyR2, for Catecholaminergic Polymorphic Ventricular Tachycardia (CPVT) and/or CaV1.2, for Long-QT Syndrome (LQTS). Recently, a novel CALM2 variant, G114R, was identified in a mother and two of her four children, all of whom died suddenly while asleep at a young age. The G114R variant impairs closure of CaV1.2 and RyR2, consistent with a CPVT and/or mild LQTS phenotype. However, the children carrying the CALM2 G114R variant displayed a phenotype commonly observed with variants in NaV1.5, i.e., Brugada Syndrome (BrS) or LQT3, where death while asleep is a common feature. We therefore hypothesized that the G114R variant specifically would interfere with NaV1.5 binding. Here, we demonstrate that CaM binding to the NaV1.5 IQ-domain is severely impaired for two CaM variants G114R and G114W. The impact was most severe at low and intermediate Ca2+ concentrations (up to 4 µM) resulting in more than a 50-fold reduction in NaV1.5 binding affinity, and a smaller 1.5 to 11-fold reduction at high Ca2+ concentrations (25-400 µM). In contrast, the arrhythmogenic CaM-N98S variant only induced a 1.5-fold reduction in NaV1.5 binding and only at 4 µM Ca2+. A non-arrhythmogenic I10T variant in CaM did not impair NaV1.5 IQ binding. These data suggest that the interaction between NaV1.5 and CaM is decreased with certain CaM variants, which may alter the cardiac sodium current, INa. Overall, these results suggest that the phenotypic spectrum of calmodulinopathies may likely expand to include BrS- and/or LQT3-like traits.

Human disease-associated calmodulin mutations alter calcineurin function through multiple mechanisms

Calmodulin (CaM) is a ubiquitous, calcium-sensing protein that regulates a multitude of processes throughout the body. In response to changes in [Ca2+], CaM modifies, activates, and deactivates enzymes and ion channels, as well as many other cellular processes. The importance of CaM is highlighted by the conservation of an identical amino acid sequence in all mammals. Alterations to CaM amino acid sequence were once thought to be incompatible with life. During the last decade modifications to the CaM protein sequence have been observed in patients suffering from life-threatening heart disease (calmodulinopathy). Thus far, inadequate or untimely interaction between mutant CaM and several proteins (LTCC, RyR2, and CaMKII) have been identified as mechanisms underlying calmodulinopathy. Given the extensive number of CaM interactions in the body, there are likely many consequences for altering CaM protein sequence. Here, we demonstrate that disease-associated CaM mutations alter the sensitivity and activity of the Ca2+-CaM-enhanced serine/threonine phosphatase calcineurin (CaN). Biophysical characterization by circular dichroism, solution NMR spectroscopy, stopped-flow kinetic measurements, and MD simulations provide mechanistic insight into mutation dysfunction as well as highlight important aspects of CaM Ca2+ signal transduction. We find that individual CaM point mutations (N53I, F89L, D129G, and F141L) impair CaN function, however, the mechanisms are not the same. Specifically, individual point mutations can influence or modify the following properties: CaM binding, Ca2+ binding, and/or Ca2+kinetics. Moreover, structural aspects of the CaNCaM complex can be altered in manners that indicate changes to allosteric transmission of CaM binding to the enzyme active site. Given that loss of CaN function can be fatal, as well as evidence that CaN modifies ion channels already associated with calmodulinopathy, our results raise the possibility that altered CaN function contributes to calmodulinopathy.

Novel Calmodulin Variant p.E46K Associated With Severe Catecholaminergic Polymorphic Ventricular Tachycardia Produces Robust Arrhythmogenicity in Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes

BACKGROUND

CaM (calmodulin) is a ubiquitously expressed, multifunctional Ca²⁺ sensor protein that regulates numerous proteins. Recently, CaM missense variants have been identified in patients with malignant inherited arrhythmias, such as long QT syndrome and catecholaminergic polymorphic ventricular tachycardia (CPVT). However, the exact mechanism of CaM-related CPVT in human cardiomyocytes remains unclear. In this study, we sought to investigate the arrhythmogenic mechanism of CPVT caused by a novel variant using human induced pluripotent stem cell (iPSC) models and biochemical assays.

METHODS

We generated iPSCs from a patient with CPVT bearing CALM2 p.E46K. As comparisons, we used 2 control lines including an isogenic line, and another iPSC line from a patient with long QT syndrome bearing CALM2 p.N98S (also reported in CPVT). Electrophysiological properties were investigated using iPSC-cardiomyocytes. We further examined the RyR2 (ryanodine receptor 2) and Ca²⁺ affinities of CaM using recombinant proteins.

RESULTS

We identified a novel de novo heterozygous variant, CALM2 p.E46K, in 2 unrelated patients with CPVT accompanied by neurodevelopmental disorders. The E46K-cardiomyocytes exhibited more frequent abnormal electrical excitations and Ca²⁺ waves than the other lines in association with increased Ca²⁺ leakage from the sarcoplasmic reticulum via RyR2. Furthermore, the [³H]ryanodine binding assay revealed that E46K-CaM facilitated RyR2 function especially by activating at low [Ca²⁺] levels. The real-time CaM-RyR2 binding analysis demonstrated that E46K-CaM had a 10-fold increased RyR2 binding affinity compared with wild-type CaM which may account for the dominant effect of the mutant CaM. Additionally, the E46K-CaM did not affect CaM-Ca²⁺ binding or L-type calcium channel function. Finally, antiarrhythmic agents, nadolol and flecainide, suppressed abnormal Ca²⁺ waves in E46K-cardiomyocytes.

CONCLUSIONS

We, for the first time, established a CaM-related CPVT iPSC-CM model which recapitulated severe arrhythmogenic features resulting from E46K-CaM dominantly binding and facilitating RyR2. In addition, the findings in iPSC-based drug testing will contribute to precision medicine.

Palmitoylation of the pore-forming subunit of Ca(v)1.2 controls channel voltage sensitivity and calcium transients in cardiac myocytes

Mammalian voltage-activated L-type Ca2+ channels, such as Ca(v)1.2, control transmembrane Ca2+ fluxes in numerous excitable tissues. Here, we report that the pore-forming α1C subunit of Ca(v)1.2 is reversibly palmitoylated in rat, rabbit, and human ventricular myocytes. We map the palmitoylation sites to two regions of the channel: The N terminus and the linker between domains I and II. Whole-cell voltage clamping revealed a rightward shift of the Ca(v)1.2 current-voltage relationship when α1C was not palmitoylated. To examine function, we expressed dihydropyridine-resistant α1C in human induced pluripotent stem cell-derived cardiomyocytes and measured Ca2+ transients in the presence of nifedipine to block the endogenous channels. The transients generated by unpalmitoylatable channels displayed a similar activation time course but significantly reduced amplitude compared to those generated by wild-type channels. We thus conclude that palmitoylation controls the voltage sensitivity of Ca(v)1.2. Given that the identified Ca(v)1.2 palmitoylation sites are also conserved in most Ca(v)1 isoforms, we propose that palmitoylation of the pore-forming α1C subunit provides a means to regulate the voltage sensitivity of voltage-activated Ca2+ channels in excitable cells.

Arrhythmia-associated Calmodulin Variants Interact with KCNQ1 to Confer Aberrant Membrane Trafficking and Function

Rationale

Missense variants in calmodulin (CaM) predispose patients to arrhythmias associated with high mortality rates. As CaM regulates several key cardiac ion channels, a mechanistic understanding of CaM variant-associated arrhythmias requires elucidating individual CaM variant effect on distinct channels. One key CaM regulatory target is the KCNQ1 (K V 7.1) voltage-gated potassium channel that underlie the I Ks current. Yet, relatively little is known as to how CaM variants interact with KCNQ1 or affect its function.

Objective

To observe how arrhythmia-associated CaM variants affect binding to KCNQ1, channel membrane trafficking, and KCNQ1 function.

Methods and Results

We combine a live-cell FRET binding assay, fluorescence trafficking assay, and functional electrophysiology to characterize >10 arrhythmia-associated CaM variants effect on KCNQ1. We identify one variant (G114W) that exhibits severely weakened binding to KCNQ1 but find that most other CaM variants interact with similar binding affinity to KCNQ1 when compared to CaM wild-type over physiological Ca 2+ ranges. We further identify several CaM variants that affect KCNQ1 and I Ks membrane trafficking and/or baseline current activation kinetics, thereby contextualizing KCNQ1 dysfunction in calmodulinopathy. Lastly, we delineate CaM variants with no effect on KCNQ1 function.

Conclusions

This study provides comprehensive functional data that reveal how CaM variants contribute to creating a pro-arrhythmic substrate by causing abnormal KCNQ1 membrane trafficking and current conduction. We find that CaM variant regulation of KCNQ1 is not uniform with effects varying from benign to significant loss of function. This study provides a new approach to collecting details of CaM binding that are key for understanding how CaM variants predispose patients to arrhythmia via the dysregulation of multiple cardiac ion channels.

Role of Ca2+ in healthy and pathologic cardiac function: from normal excitation-contraction coupling to mutations that cause inherited arrhythmia

Calcium (Ca²⁺) ions are a key second messenger involved in the rhythmic excitation and contraction of cardiomyocytes throughout the heart. Proper function of Ca²⁺-handling proteins is required for healthy cardiac function, whereas disruption in any of these can cause cardiac arrhythmias. This comprehensive review provides a broad overview of the roles of Ca²⁺-handling proteins and their regulators in healthy cardiac function and the mechanisms by which mutations in these proteins contribute to inherited arrhythmias. Major Ca²⁺ channels and Ca²⁺-sensitive regulatory proteins involved in cardiac excitation-contraction coupling are discussed, with special emphasis on the function of the RyR2 macromolecular complex. Inherited arrhythmia disorders including catecholaminergic polymorphic ventricular tachycardia, long QT syndrome, Brugada syndrome, short QT syndrome, and arrhythmogenic right-ventricular cardiomyopathy are discussed with particular emphasis on subtypes caused by mutations in Ca²⁺-handling proteins.

Calmodulin variant E140G associated with long QT syndrome impairs CaMKIIδ autophosphorylation and L-type calcium channel inactivation

Long QT syndrome (LQTS) is a human inherited heart condition that can cause life-threatening arrhythmia including sudden cardiac death. Mutations in the ubiquitous Ca²⁺-sensing protein calmodulin (CaM) are associated with LQTS, but the molecular mechanism by which these mutations lead to irregular heartbeats is not fully understood. Here, we use a multidisciplinary approach including protein biophysics, structural biology, confocal imaging, and patch-clamp electrophysiology to determine the effect of the disease-associated CaM mutation E140G on CaM structure and function. We present novel data showing that mutant-regulated CaMKIIδ kinase activity is impaired with a significant reduction in enzyme autophosphorylation rate. We report the first high-resolution crystal structure of a LQTS-associated CaM variant in complex with the CaMKIIδ peptide, which shows significant structural differences, compared to the WT complex. Furthermore, we demonstrate that the E140G mutation significantly disrupted Ca_v1.2 Ca²⁺/CaM-dependent inactivation, while cardiac ryanodine receptor (RyR2) activity remained unaffected. In addition, we show that the LQTS-associated mutation alters CaM's Ca²⁺-binding characteristics, secondary structure content, and interaction with key partners involved in excitation-contraction coupling (CaMKIIδ, Ca_v1.2, RyR2). In conclusion, LQTS-associated CaM mutation E140G severely impacts the structure-function relationship of CaM and its regulation of CaMKIIδ and Ca_v1.2. This provides a crucial insight into the molecular factors contributing to CaM-mediated arrhythmias with a central role for CaMKIIδ.

Understanding Calmodulin Variants Affecting Calcium-Dependent Inactivation of L-Type Calcium Channels through Whole-Cell Simulation of the Cardiac Ventricular Myocyte

Mutations in the calcium-sensing protein calmodulin (CaM) have been linked to two cardiac arrhythmia diseases, Long QT Syndrome 14 (LQT14) and Catecholaminergic Polymorphic Ventricular Tachycardia Type 4 (CPVT4), with varying degrees of severity. Functional characterization of the CaM mutants most strongly associated with LQT14 show a clear disruption of the calcium-dependent inactivation (CDI) of the L-Type calcium channel (LCC). CPVT4 mutants on the other hand are associated with changes in their affinity to the ryanodine receptor. In clinical studies, some variants have been associated with both CPVT4 and LQT15. This study uses simulations in a model for excitation-contraction coupling in the rat ventricular myocytes to understand how LQT14 variant might give the functional phenotype similar to CPVT4. Changing the CaM-dependent transition rate by a factor of 0.75 corresponding to the D96V variant and by a factor of 0.90 corresponding to the F142L or N98S variants, in a physiologically based stochastic model of the LCC prolonger, the action potential duration changed by a small amount in a cardiac myocyte but did not disrupt CICR at 1, 2, and 4 Hz. Under beta-adrenergic simulation abnormal excitation-contraction coupling was observed above 2 Hz pacing for the mutant CaM. The same conditions applied under beta-adrenergic stimulation led to the rapid onset of arrhythmia in the mutant CaM simulations. Simulations with the LQT14 mutations under the conditions of rapid pacing with beta-adrenergic stimulation drives the cardiac myocyte toward an arrhythmic state known as Ca2+ overload. These simulations provide a mechanistic link to a disease state for LQT14-associated mutations in CaM to yield a CPVT4 phenotype. The results show that small changes to the CaM-regulated inactivation of LCC promote arrhythmia and underscore the significance of CDI in proper heart function.

Impaired CaV1.2 inactivation reduces the efficacy of calcium channel blockers in the treatment of LQT8

Mutations in the CaV1.2 L-type calcium channel can cause a profound form of long-QT syndrome known as long-QT type 8 (LQT8), which results in cardiac arrhythmias that are often fatal in early childhood. A growing number of such pathogenic mutations in CaV1.2 have been identified, increasing the need for targeted therapies. As many of these mutations reduce channel inactivation; resulting in excess Ca2+ entry during the action potential, calcium channel blockers (CCBs) would seem to represent a promising treatment option. Yet CCBs have been unsuccessful in the treatment of LQT8. Here, we demonstrate that this lack of efficacy likely stems from the impact of the mutations on CaV1.2 channel inactivation. As CCBs are known to preferentially bind to the inactivated state of the channel, mutation-dependent deficits in inactivation result in a decrease in use-dependent block of the mutant channel. Further, application of the CCB verapamil to induced pluripotent stem cell (iPSC) derived cardiomyocytes from an LQT8 patient demonstrates that this loss of use-dependent block translates to a lack of efficacy in correcting the LQT phenotype. As a growing number of channelopathic mutations demonstrate effects on channel inactivation, reliance on state-dependent blockers may leave a growing population of patients without a viable treatment option. This biophysical understanding of the interplay between inactivation deficits and state-dependent block may provide a new avenue to guide the development of improved therapies.

A CACNA1C variant associated with cardiac arrhythmias provides mechanistic insights in the calmodulation of L-type Ca2+ channels

We recently reported the identification of a de novo single nucleotide variant in exon 9 of CACNA1C associated with prolonged repolarization interval. Recombinant expression of the glycine to arginine variant at position 419 produced a gain in the function of the L-type Ca_V1.2 channel with increased peak current density and activation gating but without significant decrease in the inactivation kinetics. We herein reveal that these properties are replicated by overexpressing calmodulin (CaM) with Ca_V1.2 WT and are reversed by exposure to the CaM antagonist W-13. Phosphomimetic (T79D or S81D), but not phosphoresistant (T79A or S81A), CaM surrogates reproduced the impact of CaM WT on the function of Ca_V1.2 WT. The increased channel activity of Ca_V1.2 WT following overexpression of CaM was found to arise in part from enhanced cell surface expression. In contrast, the properties of the variant remained unaffected by any of these treatments. Ca_V1.2 substituted with the α-helix breaking proline residue were more reluctant to open than Ca_V1.2 WT but were upregulated by phosphomimetic CaM surrogates. Our results indicate that (1) CaM and its phosphomimetic analogs promote a gain in the function of Ca_V1.2 and (2) the structural properties of the first intracellular linker of Ca_V1.2 contribute to its CaM-induced modulation. We conclude that the CACNA1C clinical variant mimics the increased activity associated with the upregulation of Ca_V1.2 by Ca²⁺-CaM, thus maintaining a majority of channels in a constitutively active mode that could ultimately promote ventricular arrhythmias.

Markov state modelling reveals heterogeneous drug-inhibition mechanism of Calmodulin

Calmodulin (CaM) is a calcium sensor which binds and regulates a wide range of target-proteins. This implicitly enables the concentration of calcium to influence many downstream physiological responses, including muscle contraction, learning and depression. The antipsychotic drug trifluoperazine (TFP) is a known CaM inhibitor. By binding to various sites, TFP prevents CaM from associating to target-proteins. However, the molecular and state-dependent mechanisms behind CaM inhibition by drugs such as TFP are largely unknown. Here, we build a Markov state model (MSM) from adaptively sampled molecular dynamics simulations and reveal the structural and dynamical features behind the inhibitory mechanism of TFP-binding to the C-terminal domain of CaM. We specifically identify three major TFP binding-modes from the MSM macrostates, and distinguish their effect on CaM conformation by using a systematic analysis protocol based on biophysical descriptors and tools from machine learning. The results show that depending on the binding orientation, TFP effectively stabilizes features of the calcium-unbound CaM, either affecting the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the bound domain. The conclusions drawn from this work may in the future serve to formulate a complete model of pharmacological modulation of CaM, which furthers our understanding of how these drugs affect signaling pathways as well as associated diseases.

Calmodulin (CaM) is a calcium-sensing protein which makes other proteins dependent on the surrounding calcium concentration by binding to these proteins. Such protein-protein interactions with CaM are vital for calcium to control many physiological pathways within the cell. The antipsychotic drug trifluoperazine (TFP) inhibits CaM’s ability to bind and regulate other proteins. Here, we use molecular dynamics simulations together with Markov state modeling and machine learning to understand the structural and dynamical features by which TFP bound to the one domain of CaM prevents association to other proteins. We find that TFP encourages CaM to adopt a conformation that is like the one stabilized in absence of calcium: depending on the binding orientation of TFP, the drug indeed either affects the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the domain. Understanding TFP binding is a first step towards designing better drugs targeting CaM.

hPSC gene editing for cardiac disease therapy

Cardiovascular diseases (CVDs) are the leading cause of mortality worldwide. However, the lack of human cardiomyocytes with proper genetic backgrounds limits the study of disease mechanisms. Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) have significantly advanced the study of these conditions. Moreover, hPSC-CMs made it easy to study CVDs using genome-editing techniques. This article discusses the applications of these techniques in hPSC for studying CVDs. Recently, several genome-editing systems have been used to modify hPSCs, including zinc finger nucleases, transcription activator-like effector nucleases, and clustered regularly interspaced short palindromic repeat-associated protein 9 (CRISPR/Cas9). We focused on the recent advancement of genome editing in hPSCs, which dramatically improved the efficiency of the cell-based mechanism study and therapy for cardiac diseases.

An ensemble of parameters from a robust Markov-based model reproduces L-type calcium currents from different human cardiac myocytes

The development of modeling structures at the channel level that can integrate subcellular and cell models and properly reproduce different experimental data is of utmost importance in cardiac electrophysiology. In contrast to gate-based models, Markov Chain models are well suited to promote the integration of the subcellular level of the cardiomyocyte to the whole cell. In this paper, we develop Markov Chain models for the L-type Calcium current that can reproduce the electrophysiology of two established human models for the ventricular and Purkinje cells. In addition, instead of presenting a single set of parameters, we present a collection of set of parameters employing Differential Evolution algorithms that can properly reproduce very different protocol data. We show the importance of using an ensemble of a set of parameter values to obtain proper results when considering a second protocol that suppresses calcium inactivation and mimics a pathological condition. We discuss how model discrepancy, data availability, and parameter identifiability can influence the choice of the size of the collection. In summary, we have modified two cardiac models by proposing new Markov Chain models for the L-type Calcium. We keep the original whole-cell dynamics by reproducing the same characteristic action potential and calcium dynamics, whereas the Markov chain-based description of the L-type Calcium channels allows novel small spatial scale simulations of subcellular processes. Finally, the use of collections of parameters was crucial for addressing model discrepancy, identifiability issues, and avoiding fitting parameters overly precisely, i.e., overfitting.

Controlling ion channel function with renewable recombinant antibodies

Selective ion channel modulators play a critical role in physiology in defining the contribution of specific ion channels to physiological function and as proof of concept for novel therapeutic strategies. Antibodies are valuable research tools that have broad uses including defining the expression and localization of ion channels in native tissue, and capturing ion channel proteins for subsequent analyses. In this review, we detail how renewable and recombinant antibodies can be used to control ion channel function. We describe the different forms of renewable and recombinant antibodies that have been used and the mechanisms by which they modulate ion channel function. We highlight the use of recombinant antibodies that are expressed intracellularly (intrabodies) as genetically-encoded tools to control ion channel function. We also offer perspectives of avenues of future research that may be opened by the application of emerging technologies for engineering recombinant antibodies for enhanced utility in ion channel research. Overall, this review provides insights that may help stimulate and guide interested researchers to develop and incorporate renewable and recombinant antibodies as valuable tools to control ion channel function. Abstract figure legend Two different approaches for controlling ion channel function using renewable recombinant antibodies. On the left, an externally applied intact IgG antibody (purple) binds to an extracellular domain of an ion channel (light blue) to control ion channel function. On the right, a genetically-encoded intrabody, in this example a camelid nanobody (green) fused to an effector molecule (red) binds to an intracellular auxiliary subunit of an ion channel (dark blue) to control ion channel function. This article is protected by copyright. All rights reserved.

Novel CALM3 Variant Causing Calmodulinopathy With Variable Expressivity in a 4-Generation Family

BACKGROUND

CaM (calmodulin), encoded by 3 separate genes (CALM1, CALM2, and CALM3), is a multifunctional Ca²⁺-binding protein involved in many signal transduction events including ion channel regulation. CaM variants may present with early-onset long QT syndrome (LQTS), catecholaminergic polymorphic ventricular tachycardia, or sudden cardiac death. Most reported variants occurred de novo. We identified a novel CALM3 variant, p.Asn138Lys (N138K), in a 4-generation family segregating with LQTS. The aim of this study was to elucidate its pathogenicity and to compare it with that of p.D130G-CaM-a variant associated with a severe LQTS phenotype.

METHODS

We performed whole exome sequencing for a large, 4-generation family affected by LQTS. To assess the effect of the detected CALM3 variant, the intrinsic Ca²⁺-binding affinity was measured by stoichiometric Ca²⁺ titrations and equilibrium titrations. L-type Ca²⁺ and slow delayed rectifier potassium currents (I_CaL and I_Ks) were recorded by whole-cell patch-clamp. Cav1.2 and Kv7.1 membrane expression were determined by optical fluorescence assays.

RESULTS

We identified 14 p.N138K-CaM carriers in a family where 2 sudden deaths occurred in children. Several members were only mildly affected compared with CaM-LQTS patients to date described in literature. The intrinsic Ca²⁺-binding affinity of the CaM C-terminal domain was 10-fold lower for p.N138K-CaM compared with wild-type-CaM. I_CaL inactivation was slowed in cells expressing p.N138K-CaM but less than in p.D130G-CaM cells. Unexpectedly, a larger I_Ks current density was observed in cells expressing p.N138K-CaM, but not for p.D130G-CaM, compared with wild-type-CaM.

CONCLUSIONS

The p.N138K CALM3 variant impairs Ca²⁺-binding affinity of CaM and I_CaL inactivation but potentiates I_Ks. The variably expressed phenotype of this variant compared with previously published de novo LQTS-CaM variants is likely explained by a milder impairment of I_CaL inactivation combined with I_Ks augmentation.

Hypocalcemia-Induced QT Interval Prolongation

An 87-year-old man with a history of transcatheter aortic valve replacement, pulmonary hypertension, diastolic dysfunction with preserved systolic function, and myelofibrosis had a 12-lead ECG showed a prolonged QT interval of 508 ms with heart-rate correction placing it in the 99th percentile of the population. Reduction in the dose of furosemide and calcium supplementation increased serum calcium and shortened the QT interval. This case provides an opportunity to examine newer concepts for the understanding of the mechanisms by which hypocalcemia might induce QT prolongation. Hypocalcemia likely produces corrected QT interval prolongation primarily through a calcium-dependent inactivation (CDI) mechanism on the L-type calcium channel (LTCC). Lower extracellular calcium leads to a decreased ICaL, subsequently causing intracellular calcium to take longer to reach the critical threshold to induce CDI of the LTCC. The resulting prolonged repolarization of the ventricular myocyte can lead to early after-depolarizations and ensuing life-threatening ventricular arrhythmias. Genetic polymorphisms in Ca2+-binding protein calmodulin which can prolong QT, underscore the role for disturbances of intracellular myocardial calcium handling in arrhythmogenesis. Hypocalcemia is an under-recognized cause of QT prolongation and should be taken into careful consideration in patients presenting with incidental findings of a prolonged QT interval.

The link between abnormalities of calcium handling proteins and catecholaminergic polymorphic ventricular tachycardia

Catecholaminergic polymorphic ventricular tachycardia (CPVT), a rare autosomal dominant or recessive disease, usually results in syncope or sudden cardiac death. Most CPVT patients do not show abnormal cardiac structure and electrocardiogram features and symptoms, usually onset during adrenergically mediated physiological conditions. CPVT tends to occur at a younger age and is not easy to be diagnosed and managed. The main cause of CPVT is associated with mishandling Ca²⁺ in cardiomyocytes. Intracellular Ca²⁺ is strictly controlled by a protein located in the sarcoplasm reticulum (SR), such as ryanodine receptor, histidine-rich Ca²⁺-binding protein, triadin, and junctin. Mutation in these proteins results in misfolding or malfunction of these proteins, thereby affecting their Ca²⁺-binding affinity, and subsequently disturbs Ca²⁺ homeostasis during excitation–contraction coupling (E-C coupling). Furthermore, transient disturbance of Ca²⁺ homeostasis increases membrane potential and causes Ca²⁺ store overload-induced Ca²⁺ release, which in turn leads to delayed after depolarization and arrhythmia. Previous studies have focused on the interaction between ryanodine receptors and protein kinase or phosphatase in the cytosol. However, recent studies showed the regulation signaling for ryanodine receptor not only from the cytosol but also within the SR. The changing of Ca²⁺ concentration is critical for protein interaction inside the SR which changes protein conformation to regulate the open probability of ryanodine receptors. Thus, it influences the threshold of Ca²⁺ released from the SR, making it easier to release Ca²⁺ during E-C coupling. In this review, we briefly discuss how Ca²⁺ handling protein variations affect the Ca²⁺ handling in CPVT.

Calmodulinopathy in inherited arrhythmia syndromes

Calmodulin (CaM) is a ubiquitous intracellular calcium sensor that controls and regulates key cellular functions. In all vertebrates, three CaM genes located on separate chromosomes encode an identical 149 amino acid protein, implying an extraordinarily high level of evolutionary importance and suggesting that CaM mutations would be possibly fatal. Inherited arrhythmia syndromes comprise a spectrum of primary electrical disorders caused by mutations in genes encoding ion channels or associated proteins leading to various cardiac arrhythmias, unexplained syncope, and sudden cardiac death. CaM mutations have emerged as an independent entity among inherited arrhythmia syndromes, referred to as calmodulinopathies. The most common clinical presentation associated with calmodulinopathy is congenital long QT syndrome, followed by catecholaminergic polymorphic ventricular tachycardia, both of which significantly increase the possibility of repeated syncope, lethal arrhythmic events, and sudden cardiac death, especially in young individuals. Here, we aim to give an overview of biochemical and structural characteristics of CaM and progress toward updating current known CaM mutations and associated clinical phenotypes. We also review the possible mechanisms underlying calmodulinopathy, based on several key in vitro studies. We expect that further experimental studies are needed to explore the complexity of calmodulinopathy.

Different arrhythmia-associated calmodulin mutations have distinct effects on cardiac SK channel regulation

Calmodulin (CaM) plays a critical role in intracellular signaling and regulation of Ca²⁺-dependent proteins and ion channels. Mutations in CaM cause life-threatening cardiac arrhythmias. Among the known CaM targets, small-conductance Ca²⁺-activated K⁺ (SK) channels are unique, since they are gated solely by beat-to-beat changes in intracellular Ca²⁺. However, the molecular mechanisms of how CaM mutations may affect the function of SK channels remain incompletely understood. To address the structural and functional effects of these mutations, we introduced prototypical human CaM mutations in human induced pluripotent stem cell–derived cardiomyocyte-like cells (hiPSC-CMs). Using structural modeling and molecular dynamics simulation, we demonstrate that human calmodulinopathy-associated CaM mutations disrupt cardiac SK channel function via distinct mechanisms. CaM_D96V and CaM_D130G mutants reduce SK currents through a dominant-negative fashion. By contrast, specific mutations replacing phenylalanine with leucine result in conformational changes that affect helix packing in the C-lobe, which disengage the interactions between apo-CaM and the CaM-binding domain of SK channels. Distinct mutant CaMs may result in a significant reduction in the activation of the SK channels, leading to a decrease in the key Ca²⁺-dependent repolarization currents these channels mediate. The findings in this study may be generalizable to other interactions of mutant CaMs with Ca²⁺-dependent proteins within cardiac myocytes.

The heart arrhythmia-linked D130G calmodulin mutation causes premature inhibitory autophosphorylation of CaMKII

The Ca²⁺/calmodulin (CaM)-dependent kinase II (CaMKII) is well known for transmitting Ca²⁺-signals, which leads to a multitude of physiological responses. Its functionality is believed to involve CaMKII holoenzyme dynamics where trans-autophosphorylation of the crucial phosphorylation site, T286 occurs. Phosphorylation of this site does not occur when stimulated exclusively with the arrhythmia associated D130G mutant form of CaM in vitro. Here, we present evidence that the loss-of-CaMKII function correlates with premature phosphorylation of its inhibitory phosphosite T306 in CaMKIIα and T307 in CaMKIIδ as this site was up to 20-fold more phosphorylated in the presence of D130G CaM compared to wildtype CaM. Indeed, changing this phosphosite to a non-phosphorylatable alanine reversed the inhibitory effect of D130G both in vitro and in live cell experiments. In addition, several phosphosites with so far undescribed functions directing the Ca²⁺-sensitivity of the CaMKII sensor were also affected by the presence of the D130G mutation implicating a role of several additional autophosphosites (besides T286 and T306/T307) so far not known to regulate CaMKII Ca²⁺ sensitivity. Furthermore, we show that introducing a D130G mutation in the CALM2 gene of the P19CL6 pluripotent mouse embryonic carcinoma cell line using CRISPR/Cas9 decreased the spontaneous beat frequency compared to wildtype cells when differentiated into cardiomyocytes supporting an alteration of cardiomyocyte physiology caused by this point mutation. In conclusion, our observations shed for the first time light on how the D130G CaM mutation interferes with the function of CaMKII and how it affects the beating frequency of cardiomyocyte-like cells.

Calcium Release Channels: Structure and Function of IP3 Receptors and Ryanodine Receptors

Ca²⁺-release channels are giant membrane proteins that control the release of Ca²⁺ from the endoplasmic and sarcoplasmic reticulum. The two members, ryanodine receptors (RyRs) and inositol-1,4,5-trisphosphate Receptors (IP₃Rs), are evolutionarily related and are both activated by cytosolic Ca²⁺. They share a common architecture, but RyRs have evolved additional modules in the cytosolic region. Their massive size allows for the regulation by tens of proteins and small molecules, which can affect the opening and closing of the channels. In addition to Ca²⁺, other major triggers include IP₃ for the IP₃Rs, and depolarization of the plasma membrane for a particular RyR subtype. Their size has made them popular targets for study via electron microscopic methods, with current structures culminating near 3Å. The available structures have provided many new mechanistic insights int the binding of auxiliary proteins and small molecules, how these can regulate channel opening, and the mechanisms of disease-associated mutations. They also help scrutinize previously proposed binding sites, as some of these are now incompatible with the structures. Many questions remain around the structural effects of post-translational modifications, additional binding partners, and the higher-order complexes these channels can make in situ. This review summarizes our current knowledge about the structures of Ca²⁺-release channels and how this informs on their function.

Mechanisms underlying pathological Ca2+ handling in diseases of the heart

Cardiomyocyte contraction relies on precisely regulated intracellular Ca²⁺ signaling through various Ca²⁺ channels and transporters. In this article, we will review the physiological regulation of Ca²⁺ handling and its role in maintaining normal cardiac rhythm and contractility. We discuss how inherited variants or acquired defects in Ca²⁺ channel subunits contribute to the development or progression of diseases of the heart. Moreover, we highlight recent insights into the role of protein phosphatase subunits and striated muscle preferentially expressed protein kinase (SPEG) in atrial fibrillation, heart failure, and cardiomyopathies. Finally, this review summarizes current drug therapies and new advances in genome editing as therapeutic strategies for the cardiac diseases caused by aberrant intracellular Ca²⁺ signaling.

Structural Basis for the Modulation of Ryanodine Receptors

Historically, ryanodine receptors (RyRs) have presented unique challenges for high-resolution structural determination despite long-standing interest in their role in excitation-contraction coupling. Owing to their large size (nearly 2.2 MDa), high-resolution structures remained elusive until the advent of cryogenic electron microscopy (cryo-EM) techniques. In recent years, structures for both RyR1 and RyR2 have been solved at near-atomic resolution. Furthermore, recent reports have delved into their more complex structural associations with key modulators - proteins such as the dihydropyridine receptor (DHPR), FKBP12/12.6, and calmodulin (CaM), as well as ions and small molecules including Ca²⁺, ATP, caffeine, and PCB95. This review addresses the modulation of RyR1 and RyR2, in addition to the impact of such discoveries on intracellular Ca²⁺ dynamics and biophysical properties.

Infanticide vs. inherited cardiac arrhythmias

AIMS

In 2003, an Australian woman was convicted by a jury of smothering and killing her four children over a 10-year period. Each child died suddenly and unexpectedly during a sleep period, at ages ranging from 19 days to 18 months. In 2019 we were asked to investigate if a genetic cause could explain the children's deaths as part of an inquiry into the mother's convictions.

METHODS AND RESULTS

Whole genomes or exomes of the mother and her four children were sequenced. Functional analysis of a novel CALM2 variant was performed by measuring Ca2+-binding affinity, interaction with calcium channels and channel function. We found two children had a novel calmodulin variant (CALM2 G114R) that was inherited maternally. Three genes (CALM1-3) encode identical calmodulin proteins. A variant in the corresponding residue of CALM3 (G114W) was recently reported in a child who died suddenly at age 4 and a sibling who suffered a cardiac arrest at age 5. We show that CALM2 G114R impairs calmodulin's ability to bind calcium and regulate two pivotal calcium channels (CaV1.2 and RyR2) involved in cardiac excitation contraction coupling. The deleterious effects of G114R are similar to those produced by G114W and N98S, which are considered arrhythmogenic and cause sudden cardiac death in children.

CONCLUSION

A novel functional calmodulin variant (G114R) predicted to cause idiopathic ventricular fibrillation, catecholaminergic polymorphic ventricular tachycardia, or mild long QT syndrome was present in two children. A fatal arrhythmic event may have been triggered by their intercurrent infections. Thus, calmodulinopathy emerges as a reasonable explanation for a natural cause of their deaths.

Calmodulin mutations and life-threatening cardiac arrhythmias: insights from the International Calmodulinopathy Registry

AIMS

Calmodulinopathies are rare life-threatening arrhythmia syndromes which affect mostly young individuals and are, caused by mutations in any of the three genes (CALM 1-3) that encode identical calmodulin proteins. We established the International Calmodulinopathy Registry (ICalmR) to understand the natural history, clinical features, and response to therapy of patients with a CALM-mediated arrhythmia syndrome.

METHODS AND RESULTS

A dedicated Case Report File was created to collect demographic, clinical, and genetic information. ICalmR has enrolled 74 subjects, with a variant in the CALM1 (n = 36), CALM2 (n = 23), or CALM3 (n = 15) genes. Sixty-four (86.5%) were symptomatic and the 10-year cumulative mortality was 27%. The two prevalent phenotypes are long QT syndrome (LQTS; CALM-LQTS, n = 36, 49%) and catecholaminergic polymorphic ventricular tachycardia (CPVT; CALM-CPVT, n = 21, 28%). CALM-LQTS patients have extremely prolonged QTc intervals (594 ± 73 ms), high prevalence (78%) of life-threatening arrhythmias with median age at onset of 1.5 years [interquartile range (IQR) 0.1-5.5 years] and poor response to therapies. Most electrocardiograms (ECGs) show late onset peaked T waves. All CALM-CPVT patients were symptomatic with median age of onset of 6.0 years (IQR 3.0-8.5 years). Basal ECG frequently shows prominent U waves. Other CALM-related phenotypes are idiopathic ventricular fibrillation (IVF, n = 7), sudden unexplained death (SUD, n = 4), overlapping features of CPVT/LQTS (n = 3), and predominant neurological phenotype (n = 1). Cardiac structural abnormalities and neurological features were present in 18 and 13 patients, respectively.

CONCLUSION

Calmodulinopathies are largely characterized by adrenergically-induced life-threatening arrhythmias. Available therapies are disquietingly insufficient, especially in CALM-LQTS. Combination therapy with drugs, sympathectomy, and devices should be considered.

Complex Arrhythmia Syndrome in a Knock-In Mouse Model Carrier of the N98S Calm1 Mutation

BACKGROUND

Calmodulin mutations are associated with arrhythmia syndromes in humans. Exome sequencing previously identified a de novo mutation in CALM1 resulting in a p.N98S substitution in a patient with sinus bradycardia and stress-induced bidirectional ventricular ectopy. The objectives of the present study were to determine if mice carrying the N98S mutation knocked into Calm1 replicate the human arrhythmia phenotype and to examine arrhythmia mechanisms.

METHODS

Mouse lines heterozygous for the Calm1^N98S allele (Calm1^N98S/+) were generated using CRISPR/Cas9 technology. Adult mutant mice and their wildtype littermates (Calm1^+/+) underwent electrocardiographic monitoring. Ventricular de- and repolarization was assessed in isolated hearts using optical voltage mapping. Action potentials and whole-cell currents and [Ca²⁺]_i, as well, were measured in single ventricular myocytes using the patch-clamp technique and fluorescence microscopy, respectively. The microelectrode technique was used for in situ membrane voltage monitoring of ventricular conduction fibers.

RESULTS

Two biologically independent knock-in mouse lines heterozygous for the Calm1^N98S allele were generated. Calm1^N98S/+ mice of either sex and line exhibited sinus bradycardia, QT_c interval prolongation, and catecholaminergic bidirectional ventricular tachycardia. Male mutant mice also showed QRS widening. Pharmacological blockade and activation of β-adrenergic receptors rescued and exacerbated, respectively, the long-QT phenotype of Calm1^N98S/+ mice. Optical and electric assessment of membrane potential in isolated hearts and single left ventricular myocytes, respectively, revealed β-adrenergically induced delay of repolarization. β-Adrenergic stimulation increased peak density, slowed inactivation, and left-shifted the activation curve of I_Ca.L significantly more in Calm1^N98S/+ versus Calm1^+/+ ventricular myocytes, increasing late I_Ca.L in the former. Rapidly paced Calm1^N98S/+ ventricular myocytes showed increased propensity to delayed afterdepolarization-induced triggered activity, whereas in situ His-Purkinje fibers exhibited increased susceptibility for pause-dependent early afterdepolarizations. Epicardial mapping of Calm1^N98S/+ hearts showed that both reentry and focal mechanisms contribute to arrhythmogenesis.

CONCLUSIONS

Heterozygosity for the Calm1^N98S mutation is causative of an arrhythmia syndrome characterized by sinus bradycardia, QRS widening, adrenergically mediated QTc interval prolongation, and bidirectional ventricular tachycardia. β-Adrenergically induced I_Ca.L dysregulation contributes to the long-QT phenotype. Pause-dependent early afterdepolarizations and tachycardia-induced delayed afterdepolarizations originating in the His-Purkinje network and ventricular myocytes, respectively, constitute potential sources of arrhythmia in Calm1^N98S/+ hearts.

Impact of arrhythmogenic calmodulin variants on small conductance Ca2+ -activated K+ (SK3) channels

Calmodulin (CaM) is a ubiquitous Ca²⁺‐sensing protein regulating many important cellular processes. Several CaM‐associated variants have been identified in a small group of patients with cardiac arrhythmias. The mechanism remains largely unknown, even though a number of ion channels, including the ryanodine receptors and the L‐type calcium channels have been shown to be functionally affected by the presence of mutant CaM. CaM is constitutively bound to the SK channel, which underlies the calcium‐gated I_SK contributing to cardiac repolarization. The CaM binding to SK channels is essential for gating, correct assembly, and membrane expression. To elucidate the effect of nine different arrhythmogenic CaM variants on SK3 channel function, HEK293 cells stably expressing SK3 were transiently co‐transfected with CaM^{WT or variant} and whole‐cell patch‐clamp recordings were performed with a calculated free Ca²⁺ concentration of 400 nmol/L. MDCK cells were transiently transfected with SK3 and/or CaM^{WT or variant} to address SK3 and CaM localization by immunocytochemistry. The LQTS‐associated variants CaM^D96V, CaM^D130G, and CaM^F142L reduced I_SK,Ca compared with CaM^WT (P < 0.01, P < 0.001, and P < 0.05, respectively). The CPVT associated variant CaM^N54I also reduced the I_SK,Ca (P < 0.05), which was linked to an accumulation of SK3/CaM^N54I channel complexes in intracellular compartments (P < 0.05). The CPVT associated variants, CaM^A103V and CaM^D132E only revealed a tendency toward reduced current, while the variants CaM^F90L and CaM^N98S, causing LQTS syndrome, did not have any impact on I_SK,Ca. In conclusion, we found that the arrhythmogenic CaM variants CaM^N54I, CaM^D96V, CaM^D130G, and CaM^F142L significantly down‐regulate the SK3 channel current, but with distinct mechanism.

CaV channels reject signaling from a second CaM in eliciting Ca2+-dependent feedback regulation

Calmodulin (CaM) regulation of voltage-gated calcium (Ca_V1-2) channels is a powerful Ca²⁺-feedback mechanism to adjust channel activity in response to Ca²⁺ influx. Despite progress in resolving mechanisms of CaM-Ca_V feedback, the stoichiometry of CaM interaction with Ca_V channels remains ambiguous. Functional studies that tethered CaM to Ca_V1.2 suggested that a single CaM sufficed for Ca²⁺ feedback, yet biochemical, FRET, and structural studies showed that multiple CaM molecules interact with distinct interfaces within channel cytosolic segments, suggesting that functional Ca²⁺ regulation may be more nuanced. Resolving this ambiguity is critical as CaM is enriched in subcellular domains where Ca_V channels reside, such as the cardiac dyad. We here localized multiple CaMs to the Ca_V nanodomain by tethering either WT or mutant CaM that lack Ca²⁺-binding capacity to the pore-forming α-subunit of Ca_V1.2, Ca_V1.3, and Ca_V2.1 and/or the auxiliary β_2A subunit. We observed that a single CaM tethered to either the α or β_2A subunit tunes Ca²⁺ regulation of Ca_V channels. However, when multiple CaMs are localized concurrently, Ca_V channels preferentially respond to signaling from the α-subunit-tethered CaM. Mechanistically, the introduction of a second IQ domain to the Ca_V1.3 carboxyl tail switched the apparent functional stoichiometry, permitting two CaMs to mediate functional regulation. In all, Ca²⁺ feedback of Ca_V channels depends exquisitely on a single CaM preassociated with the α-subunit carboxyl tail. Additional CaMs that colocalize with the channel complex are unable to trigger Ca²⁺-dependent feedback of channel gating but may support alternate regulatory functions.

Regulation of cardiovascular calcium channel activity by post-translational modifications or interacting proteins

Voltage-gated calcium channels are the major pathway for Ca²⁺ influx to initiate the contraction of smooth and cardiac muscles. Alterations of calcium channel function have been implicated in multiple cardiovascular diseases, such as hypertension, atrial fibrillation, and long QT syndrome. Post-translational modifications do expand cardiovascular calcium channel structure and function to affect processes such as channel trafficking or polyubiquitination by two E3 ubiquitin ligases, Ret finger protein 2 (Rfp2) or murine double minute 2 protein (Mdm2). Additionally, biophysical property such as Ca²⁺-dependent inactivation (CDI) could be altered through binding of calmodulin, or channel activity could be modulated via S-nitrosylation by nitric oxide and phosphorylation by protein kinases or by interacting protein partners, such as galectin-1 and Rem. Understanding how cardiovascular calcium channel function is post-translationally remodeled under distinctive disease conditions will provide better information about calcium channel-related disease mechanisms and improve the development of more selective therapeutic agents for cardiovascular diseases.

The LQT-associated calmodulin mutant E141G induces disturbed Ca2+-dependent binding and a flickering gating mode of the CaV1.2 channel

Calmodulin (CaM) mutations are associated with congenital long QT (LQT) syndrome (LQTS), which may be related to the dysregulation of the cardiac-predominant Ca²⁺ channel isoform Ca_V1.2. Among various mutants, CaM-E141G was identified as a critical missense variant. However, the interaction of this CaM mutant with the Ca_V1.2 channel has not been determined. In this study, by utilizing a semiquantitative pull-down assay, we explored the interaction of CaM-E141G with CaM-binding peptide fragments of the Ca_V1.2 channel. Using the patch-clamp technique, we also investigated the electrophysiological effects of the mutant on Ca_V1.2 channel activity. We found that the maximum binding (B_max) of CaM-E141G to the proximal COOH-terminal region, PreIQ-IQ, PreIQ, IQ, and NT (an NH₂-terminal peptide) was decreased (by 17.71-59.26%) compared with that of wild-type CaM (CaM-WT). In particular, the Ca²⁺-dependent increase in B_max became slower with the combination of CaM-E141G + PreIQ and IQ but faster in the case of NT. Functionally, CaM-WT and CaM-E141G at 500 nM Ca²⁺ decreased Ca_V1.2 channel activity to 24.88% and 55.99%, respectively, compared with 100 nM Ca²⁺, showing that the inhibitory effect was attenuated in CaM-E141G. The mean open time of the Ca_V1.2 channel was increased, and the number of blank traces with no channel opening was significantly decreased. Overall, CaM-E141G exhibits disrupted binding with the Ca_V1.2 channel and induces a flickering gating mode, which may result in the dysfunction of the Ca_V1.2 channel and, thus, the development of LQTS. The present study is the first to investigate the detailed binding properties and single-channel gating mode induced by the interaction of CaM-E141G with the Ca_V1.2 channel.

Molecular and tissue mechanisms of catecholaminergic polymorphic ventricular tachycardia

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is a stress-induced cardiac channelopathy that has a high mortality in untreated patients. Our understanding has grown tremendously since CPVT was first described as a clinical syndrome in 1995. It is now established that the deadly arrhythmias are caused by unregulated 'pathological' calcium release from the sarcoplasmic reticulum (SR), the major calcium storage organelle in striated muscle. Important questions remain regarding the molecular mechanisms that are responsible for the pathological calcium release, regarding the tissue origin of the arrhythmic beats that initiate ventricular tachycardia, and regarding optimal therapeutic approaches. At present, mutations in six genes involved in SR calcium release have been identified as the genetic cause of CPVT: RYR2 (encoding ryanodine receptor calcium release channel), CASQ2 (encoding cardiac calsequestrin), TRDN (encoding triadin), CALM1, CALM2 and CALM3 (encoding identical calmodulin protein). Here, we review each CPVT subtype and how CPVT mutations alter protein function, RyR2 calcium release channel regulation, and cellular calcium handling. We then discuss research and hypotheses surrounding the tissue mechanisms underlying CPVT, such as the pathophysiological role of sinus node dysfunction in CPVT, and whether the arrhythmogenic beats originate from the conduction system or the ventricular working myocardium. Finally, we review the treatments that are available for patients with CPVT, their efficacy, and how therapy could be improved in the future.

The arrhythmogenic N53I variant subtly changes the structure and dynamics in the calmodulin N-terminal domain, altering its interaction with the cardiac ryanodine receptor

Mutations in the genes encoding the highly conserved Ca²⁺-sensing protein calmodulin (CaM) cause severe cardiac arrhythmias, including catecholaminergic polymorphic ventricular tachycardia or long QT syndrome and sudden cardiac death. Most of the identified arrhythmogenic mutations reside in the C-terminal domain of CaM and mostly affect Ca²⁺-coordinating residues. One exception is the catecholaminergic polymorphic ventricular tachycardia-causing N53I substitution, which resides in the N-terminal domain (N-domain). It does not affect Ca²⁺ coordination and has only a minor impact on binding affinity toward Ca²⁺ and on other biophysical properties. Nevertheless, the N53I substitution dramatically affects CaM's ability to reduce the open probability of the cardiac ryanodine receptor (RyR2) while having no effect on the regulation of the plasmalemmal voltage-gated Ca²⁺ channel, Ca_v1.2. To gain more insight into the molecular disease mechanism of this mutant, we used NMR to investigate the structures and dynamics of both apo- and Ca²⁺-bound CaM-N53I in solution. We also solved the crystal structures of WT and N53I CaM in complex with the primary calmodulin-binding domain (CaMBD2) from RyR2 at 1.84-2.13 Å resolutions. We found that all structures of the arrhythmogenic CaM-N53I variant are highly similar to those of WT CaM. However, we noted that the N53I substitution exposes an additional hydrophobic surface and that the intramolecular dynamics of the protein are significantly altered such that they destabilize the CaM N-domain. We conclude that the N53I-induced changes alter the interaction of the CaM N-domain with RyR2 and thereby likely cause the arrhythmogenic phenotype of this mutation.

Cardiomyocyte calcium handling in health and disease: Insights from in vitro and in silico studies

Calcium (Ca²⁺) plays a central role in cardiomyocyte excitation-contraction coupling. To ensure an optimal electrical impulse propagation and cardiac contraction, Ca²⁺ levels are regulated by a variety of Ca²⁺-handling proteins. In turn, Ca²⁺ modulates numerous electrophysiological processes. Accordingly, Ca²⁺-handling abnormalities can promote cardiac arrhythmias via various mechanisms, including the promotion of afterdepolarizations, ion-channel modulation and structural remodeling. In the last 30 years, significant improvements have been made in the computational modeling of cardiomyocyte Ca²⁺ handling under physiological and pathological conditions. However, numerous questions involving the Ca²⁺-dependent regulation of different macromolecular complexes, cross-talk between Ca²⁺-dependent regulatory pathways operating over a wide range of time scales, and bidirectional interactions between electrophysiology and mechanics remain to be addressed by in vitro and in silico studies. A better understanding of disease-specific Ca²⁺-dependent proarrhythmic mechanisms may facilitate the development of improved therapeutic strategies. In this review, we describe the fundamental mechanisms of cardiomyocyte Ca²⁺ handling in health and disease, and provide an overview of currently available computational models for cardiomyocyte Ca²⁺ handling. Finally, we discuss important uncertainties and open questions about cardiomyocyte Ca²⁺ handling and highlight how synergy between in vitro and in silico studies may help to answer several of these issues.

Calcium Handling Defects and Cardiac Arrhythmia Syndromes

Calcium ions (Ca²⁺) play a major role in the cardiac excitation-contraction coupling. Intracellular Ca²⁺ concentration increases during systole and falls in diastole thereby determining cardiac contraction and relaxation. Normal cardiac function also requires perfect organization of the ion currents at the cellular level to drive action potentials and to maintain action potential propagation and electrical homogeneity at the tissue level. Any imbalance in Ca²⁺ homeostasis of a cardiac myocyte can lead to electrical disturbances. This review aims to discuss cardiac physiology and pathophysiology from the elementary membrane processes that can cause the electrical instability of the ventricular myocytes through intracellular Ca²⁺ handling maladies to inherited and acquired arrhythmias. Finally, the paper will discuss the current therapeutic approaches targeting cardiac arrhythmias.

Calmodulin Mutations Associated with Heart Arrhythmia: A Status Report

Calmodulin (CaM) is a ubiquitous intracellular Ca²⁺ sensing protein that modifies gating of numerous ion channels. CaM has an extraordinarily high level of evolutionary conservation, which led to the fundamental assumption that mutation would be lethal. However, in 2012, complete exome sequencing of infants suffering from recurrent cardiac arrest revealed de novo mutations in the three human CALM genes. The correlation between mutations and pathophysiology suggests defects in CaM-dependent ion channel functions. Here, we review the current state of the field for all reported CaM mutations associated with cardiac arrhythmias, including knowledge of their biochemical and structural characteristics, and progress towards understanding how these mutations affect cardiac ion channel function.

Arrhythmia mutations in calmodulin can disrupt cooperativity of Ca2+ binding and cause misfolding

KEY POINTS

Mutations in the calmodulin protein (CaM) are associated with arrhythmia syndromes. This study focuses on understanding the structural characteristics of CaM disease mutants and their interactions with the voltage-gated calcium channel Ca_V 1.2. Arrhythmia mutations in CaM can lead to loss of Ca²⁺ binding, uncoupling of Ca²⁺ binding cooperativity, misfolding of the EF-hands and altered affinity for the calcium channel. These results help us to understand how different CaM mutants have distinct effects on structure and interactions with protein targets to cause disease.

ABSTRACT

Calmodulinopathies are life-threatening arrhythmia syndromes that arise from mutations in calmodulin (CaM), a calcium sensing protein whose sequence is completely conserved across all vertebrates. These mutations have been shown to interfere with the function of cardiac ion channels, including the voltage-gated Ca²⁺ channel Ca_V 1.2 and the ryanodine receptor (RyR2), in a mutation-specific manner. The ability of different CaM disease mutations to discriminate between these channels has been enigmatic. We present crystal structures of several C-terminal lobe mutants and an N-terminal lobe mutant in complex with the Ca_V 1.2 IQ domain, in conjunction with binding assays and complementary structural biology techniques. One mutation (D130G) causes a pathological conformation, with complete separation of EF-hands within the C-lobe and loss of Ca²⁺ binding in EF-hand 4. Another variant (Q136P) has severely reduced affinity for the IQ domain, and shows changes in the CD spectra under Ca²⁺ -saturating conditions when unbound to the IQ domain. Ca²⁺ binding to a pair of EF-hands normally proceeds with very high cooperativity, but we found that N98S CaM can adopt different conformations with either one or two Ca²⁺ ions bound to the C-lobe, possibly disrupting the cooperativity. An N-lobe variant (N54I), which causes severe stress-induced arrhythmia, does not show any major changes in complex with the IQ domain, providing a structural basis for why this mutant does not affect function of Ca_V 1.2. These findings show that different CaM mutants have distinct effects on both the CaM structure and interactions with protein targets, and act via distinct pathological mechanisms to cause disease.

Molecular Insights from Conformational Ensembles via Machine Learning

Biomolecular simulations are intrinsically high dimensional and generate noisy data sets of ever-increasing size. Extracting important features from the data is crucial for understanding the biophysical properties of molecular processes, but remains a big challenge. Machine learning (ML) provides powerful dimensionality reduction tools. However, such methods are often criticized as resembling black boxes with limited human-interpretable insight. We use methods from supervised and unsupervised ML to efficiently create interpretable maps of important features from molecular simulations. We benchmark the performance of several methods, including neural networks, random forests, and principal component analysis, using a toy model with properties reminiscent of macromolecular behavior. We then analyze three diverse biological processes: conformational changes within the soluble protein calmodulin, ligand binding to a G protein-coupled receptor, and activation of an ion channel voltage-sensor domain, unraveling features critical for signal transduction, ligand binding, and voltage sensing. This work demonstrates the usefulness of ML in understanding biomolecular states and demystifying complex simulations.

Role of cardiac ryanodine receptor calmodulin-binding domains in mediating the action of arrhythmogenic calmodulin N-domain mutation N54I

The Ca²⁺ -sensing protein calmodulin (CaM) inhibits cardiac ryanodine receptor (RyR2)-mediated Ca²⁺ release. CaM mutations associated with arrhythmias and sudden cardiac death have been shown to diminish CaM-dependent inhibition of RyR2, but the underlying mechanisms are not well understood. Nearly all arrhythmogenic CaM mutations identified are located in the C-domain of CaM and exert marked effects on Ca²⁺ binding to CaM and on the CaM C-domain interaction with the CaM-binding domain 2 (CaMBD2) in RyR2. Interestingly, the arrhythmogenic N-domain mutation CaM-N54I has little or no effect on Ca²⁺ binding to CaM or the CaM C-domain-RyR2 CaMBD2 interaction, unlike all CaM C-domain mutations. This suggests that CaM-N54I may diminish CaM-dependent RyR2 inhibition by affecting CaM N-domain interactions with RyR2 CaMBDs other than CaMBD2. To explore this possibility, we assessed the effects of deleting each of the four known CaMBDs in RyR2 (CaMBD1a, -1b, -2, or -3) on the CaM-dependent inhibition of RyR2-mediated Ca²⁺ release in HEK293 cells. We found that removing CaMBD1a, CaMBD1b, or CaMBD3 did not alter the effects of CaM-N54I or CaM-WT on RyR2 inhibition. On the other hand, deleting RyR2-CaMBD2 abolished the effects of both CaM-N54I and CaM-WT. Our results support that CaM-N54I causes aberrant RyR2 regulation via an uncharacterized CaMBD or less likely CaMBD2, and that RyR2 CaMBD2 is required for the actions of both N- and C-domain CaM mutations. Moreover, our results show that CaMBD1a is central to RyR2 regulation, but CaMBD1a, CaMBD1b, and CaMBD3 are not required for CaM-dependent inhibition of RyR2 in HEK293 cells.