Genetic Mosaicism in Calmodulinopathy

Role of Calmodulin in Cardiac Disease: Insights on Genotype and Phenotype

Calmodulin, a protein critically important for the regulation of all major cardiac ion channels, is the quintessential cellular calcium sensor and plays a key role in preserving cardiac electrical stability. Its unique importance is highlighted by the presence of 3 genes in 3 different chromosomes encoding for the same protein and by their extreme conservation. Indeed, all 3 calmodulin (CALM) genes are among the most constrained genes in the human genome, that is, the observed variants are much less than expected by chance. Not surprisingly, CALM variants are poorly tolerated and accompany significant clinical phenotypes, of which the most important are those associated with increased risk for life-threatening arrhythmias. Here, we review the current knowledge about calmodulin, its specific physiological, structural, and functional characteristics, and its importance for cardiovascular disease. Given our role in the development of this knowledge, we also share some of our views about currently unanswered questions, including the rational approaches to the clinical management of the affected patients. Specifically, we present some of the most critical information emerging from the International Calmodulinopathy Registry, which we established 10 years ago. Further progress clearly requires deep phenotypic information on as many carriers as possible through international contributions to the registry, in order to expand our knowledge about Calmodulinopathies and guide clinical management.

Molecular and functional profiling of cell diversity and identity in the lateral superior olive, an auditory brainstem center with ascending and descending projections

The lateral superior olive (LSO), a prominent integration center in the auditory brainstem, contains a remarkably heterogeneous population of neurons. Ascending neurons, predominantly principal neurons (pLSOs), process interaural level differences for sound localization. Descending neurons (lateral olivocochlear neurons, LOCs) provide feedback into the cochlea and are thought to protect against acoustic overload. The molecular determinants of the neuronal diversity in the LSO are largely unknown. Here, we used patch-seq analysis in mice at postnatal days P10-12 to classify developing LSO neurons according to their functional and molecular profiles. Across the entire sample (n = 86 neurons), genes involved in ATP synthesis were particularly highly expressed, confirming the energy expenditure of auditory neurons. Two clusters were identified, pLSOs and LOCs. They were distinguished by 353 differentially expressed genes (DEGs), most of which were novel for the LSO. Electrophysiological analysis confirmed the transcriptomic clustering. We focused on genes affecting neuronal input-output properties and validated some of them by immunohistochemistry, electrophysiology, and pharmacology. These genes encode proteins such as osteopontin, Kv11.3, and Kvβ3 (pLSO-specific), calcitonin-gene-related peptide (LOC-specific), or Kv7.2 and Kv7.3 (no DEGs). We identified 12 "Super DEGs" and 12 genes showing "Cluster similarity." Collectively, we provide fundamental and comprehensive insights into the molecular composition of individual ascending and descending neurons in the juvenile auditory brainstem and how this may relate to their specific functions, including developmental aspects.

Case report of a child with long QT syndrome type 14 caused by CALM1 gene mutation and literature review

OBJECTIVE

To analyze the clinical and genetic characteristics of a patient with long QT syndrome type 14 (long QT syndrome-14, LQT14, OMIM # 616247) caused by a de novo CALM1 mutation.

METHODS

The clinical data of the patient were collected, next-generation sequencing technology was used to determine the exome gene sequence of the patient, and the suspected pathogenic locus was verified by Sanger sequencing.

RESULTS

A 5-year and 9-month-old girl was admitted to the hospital due to a syncopal episode. During the attack, the main symptoms were loss of consciousness, cyanosis of the face and lips, and weakness of limbs. The child had multiple seizures in the past, all of which occurred after emotional excitement and activity. She was diagnosed with epilepsy for more than 3 years, but the effect of antiepileptic treatment was not satisfactory. The electrocardiogram was normal in the past. A month ago, convulsions occurred again after exercise, and the electrocardiogram showed QTc 496 ms. The treadmill test showed a significant prolongation of QTc after exercise, and the genetic results suggested a new heterozygous variant of CALM1, c.395A>G; p. (Asp132Gly). Consequently, she was diagnosed with LQT14 and treated with propranolol. During a follow-up of 15 months, there were no seizures or syncope.

CONCLUSIONS

This patient had multiple episodes of convulsions or syncope after emotional stimulation or activity, with intermittent prolongation of the QTc on routine ECG, marked prolongation of the QTc after exercise, and T-wave alternans, which differed from the LQT14 phenotype caused by the previous CALM1 mutation.

Calmodulinopathy in Japanese Children - Their Cardiac Phenotypes Are Severe and Show Early Onset in Fetal Life and Infancy

BACKGROUND

Cardiac calmodulinopathy, characterized by a life-threatening arrhythmia and sudden death in the young, is extremely rare and caused by genes encoding calmodulin, namely calmodulin 1 (CALM1), CALM2, and CALM3.

METHODS AND RESULTS

We screened 195 symptomatic children (age 0-12 years) who were suspected of inherited arrhythmias for 48 candidate genes, using a next-generation sequencer. Ten probands were identified as carrying variants in any of CALM1-3 (5%; median age 5 years), who were initially diagnosed with long QT syndrome (LQTS; n=5), catecholaminergic polymorphic ventricular tachycardia (CPVT; n=3), and overlap syndrome (n=2). Two probands harbored a CALM1 variant and 8 probands harbored 6 CALM2 variants. There were 4 clinical phenotypes: (1) documented lethal arrhythmic events (LAEs): 4 carriers of N98S in CALM1 or CALM2; (2) suspected LAEs: CALM2 p.D96G and D132G carriers experienced syncope and transient cardiopulmonary arrest under emotional stimulation; (3) critical cardiac complication: CALM2 p.D96V and p.E141K carriers showed severe cardiac dysfunction with QTc prolongation; and (4) neurological and developmental disorders: 2 carriers of CALM2 p.E46K showed cardiac phenotypes of CPVT. Beta-blocker therapy was effective in all cases except cardiac dysfunction, especially in combination with flecainide (CPVT-like phenotype) and mexiletine (LQTS-like).

CONCLUSIONS

Calmodulinopathy patients presented severe cardiac features, and their onset of LAEs was earlier in life, requiring diagnosis and treatment at the earliest age possible.

CALM3 affects the prognosis of leukemia and hemorrhoids

Leukemia is an abnormal proliferation of white blood cells in the bone marrow, resulting in a large accumulation of abnormal leukemia cells in the blood and bone marrow. Hemorrhoids are dilated and swollen veins in the rectum or anal area. However, the relationship between CALM3 and leukemia and hemorrhoids remains unclear. The hemorrhoids dataset GSE154650 and leukemia dataset GSE26294 were downloaded from GEO databases generated by GPL20301 and GPL571.The R package limma was used to screen differentially expressed genes (DEDs). Weighted gene co-expression network analysis (WGCNA) was performed. The construction and analysis of protein-protein interaction (PPI) network, functional enrichment analysis, Gene Set Enrichment Analysis (GSEA) and comparative toxicogenomics database (CTD) analysis were performed. TargetScan was used to screen miRNAs regulating central DEGs. It was verified by western blot basic cell assay. A total of 125 DEGs were co-identified. According to the GO analysis, they are mainly enriched in small molecule catabolic processes, skin development, and chemokine receptor binding. The KEGG analysis results show that the target cells are mainly enriched in the interaction of cytokines and cytokine receptors, as well as butyric acid metabolism. The GSEA analysis results indicate enrichment in small molecule catabolic processes, skin development, and chemokine receptor binding. Six core genes (CALM3, ACE2, PPARGC1A, XCR1, CFTR, PRKCA) were identified. We found that the core gene CALM3 is highly expressed in hemorrhoid samples, low in leukemia samples, and has low expression in normal samples, which may play a regulatory role in hemorrhoids and leukemia. Immunoinfiltration results showed a higher proportion of T_cells_CD4_memory_resting and a correlation with T_cells_CD8. WB experiment verified the result. CALM3 expression is low in leukemia, and the lower the expression is, the worse the prognosis is. CALM3 is highly expressed in hemorrhoids, and the higher the expression, the worse the prognosis.

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.

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.

Single-cell analysis of adult human heart across healthy and cardiovascular disease patients reveals the cellular landscape underlying SARS-CoV-2 invasion of myocardial tissue through ACE2

BACKGROUND

The distribution of ACE2 and accessory proteases (ANAD17 and CTSL) in cardiovascular tissue and the host cell receptor binding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are crucial to understanding the virus's cell invasion, which may play a significant role in determining the viral tropism and its clinical manifestations.

METHODS

We conducted a comprehensive analysis of the cell type-specific expression of ACE2, ADAM17, and CTSL in myocardial tissue from 10 patients using RNA sequencing. Our study included a meta-analysis of 2 heart single-cell RNA-sequencing studies with a total of 90,024 cells from 250 heart samples of 10 individuals. We used co-expression analysis to locate specific cell types that SARS-CoV-2 may invade.

RESULTS

Our results revealed cell-type specific associations between male gender and the expression levels of ACE2, ADAM17, and CTSL, including pericytes and fibroblasts. AGT, CALM3, PCSK5, NRP1, and LMAN were identified as potential accessory proteases that might facilitate viral invasion. Enrichment analysis highlighted the extracellular matrix interaction pathway, adherent plaque pathway, vascular smooth muscle contraction inflammatory response, and oxidative stress as potential immune pathways involved in viral infection, providing potential molecular targets for therapeutic intervention. We also found specific high expression of IFITM3 and AGT in pericytes and differences in the IFN-II signaling pathway and PAR signaling pathway in fibroblasts from different cardiovascular comorbidities.

CONCLUSIONS

Our data indicated possible high-risk groups for COVID-19 and provided emerging avenues for future investigations of its pathogenesis.

TRIAL REGISTRATION

(Not applicable).

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.

Calcium Ion Binding to the Mutants of Calmodulin: A Structure-Based Computational Predictive Model of Binding Affinity Using a Charge Scaling Approach in Molecular Dynamics Simulation

The binding of calcium ions (Ca²⁺) to the calcium-binding proteins (CBPs) controls a plethora of regulatory processes. Among the roles played by CBPs in several diseases, the onset and progress of some cardiovascular diseases are caused by mutations in calmodulin (CaM), an important member of CBPs. Rationalization and prediction of the binding affinity of Ca²⁺ ions to the CaM can play important roles in understanding the origin of cardiovascular diseases. However, there is no robust structure-based computational method for predicting the binding affinity of Ca²⁺ ions to the different forms of CBPs in general and CaM in particular. In the current work, we have devised a fast yet accurate computational technique to accurately calculate the binding affinity of Ca²⁺ to the different forms of CaM. This method combines the well-known molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method and a charge scaling approach developed by previous authors that takes care of the polarization of CaM and Ca²⁺ ions. Our detailed analysis of the different components of binding free energy shows that subtle changes in electrostatics and van der Waals contribute to the difference in the binding affinity of mutants from that of the wild type (WT), and the charge scaling approach is superior in calculating these subtle changes in electrostatics as compared to the nonpolarizable force field used in this work. A statistically significant regression model made from our binding free energy calculations gives a correlation coefficient close to 0.8 to the experimental results. This structure-based predictive model can open up a new strategy to understand and predict the binding of Ca²⁺ to the mutants of CBPs, in general.

Human Engineered Heart Tissue Models for Disease Modeling and Drug Discovery

The emergence of human induced pluripotent stem cells (hiPSCs) and efficient differentiation of hiPSC-derived cardiomyocytes (hiPSC-CMs) induced from diseased donors have the potential to recapitulate the molecular and functional features of the human heart. Although the immaturity of hiPSC-CMs, including the structure, gene expression, conduct, ion channel density, and Ca²⁺ kinetics, is a major challenge, various attempts to promote maturation have been effective. Three-dimensional cardiac models using hiPSC-CMs have achieved these functional and morphological maturations, and disease models using patient-specific hiPSC-CMs have furthered our understanding of the underlying mechanisms and effective therapies for diseases. Aside from the mechanisms of diseases and drug responses, hiPSC-CMs also have the potential to evaluate the safety and efficacy of drugs in a human context before a candidate drug enters the market and many phases of clinical trials. In fact, novel drug testing paradigms have suggested that these cells can be used to better predict the proarrhythmic risk of candidate drugs. In this review, we overview the current strategies of human engineered heart tissue models with a focus on major cardiac diseases and discuss perspectives and future directions for the real application of hiPSC-CMs and human engineered heart tissue for disease modeling, drug development, clinical trials, and cardiotoxicity tests.

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.

Potential protein markers in children with Autistic Spectrum Disorder (ASD) revealed by salivary proteomics

The lack of specific pharmacological therapy for Autistic Spectrum Disorder (ASD) and its clinical heterogeneity demand efforts directed toward the identification of biomarkers to aid in diagnosis. Proteomics offers a new perspective for studying the altered proteins associated with autism spectrum disorders (ASD) and we have saliva as an easy-to-collect biological fluid with important biomolecules for investigating biomarkers in various diseases. In this sense, saliva could be used to identify potential biomarkers of ASD. In the current work, saliva samples were collected from children with different degrees of ASD and healthy children and proteomics approaches were applied to generate data on differentially expressed proteins between groups which will serve as a basis for future validation studies as protein markers. Data are available via ProteomeXchange with identifier PXD030065. As results, 132 proteins were present in 80% of the saliva pools of all analyzed groups. Twenty-five proteins were identified as overexpressed in the group of severe and mild/moderate ASD carriers, among which, eight were identified as potential biomarkers for ASD.

Optimizing peptide inhibitors of SARS-Cov-2 nsp10/nsp16 methyltransferase predicted through molecular simulation and machine learning

Coronaviruses, including the recent pandemic strain SARS-Cov-2, use a multifunctional 2'-O-methyltransferase (2'-O-MTase) to restrict the host defense mechanism and to methylate RNA. The nonstructural protein 16 2'-O-MTase (nsp16) becomes active when nonstructural protein 10 (nsp10) and nsp16 interact. Novel peptide drugs have shown promise in the treatment of numerous diseases and new research has established that nsp10 derived peptides can disrupt viral methyltransferase activity via interaction of nsp16. This study had the goal of optimizing new analogous nsp10 peptides that have the ability to bind nsp16 with equal to or higher affinity than those naturally occurring. The following research demonstrates that in silico molecular simulations can shed light on peptide structures and predict the potential of new peptides to interrupt methyltransferase activity via the nsp10/nsp16 interface. The simulations suggest that misalignments at residues F68, H80, I81, D94, and Y96 or rotation at H80 abrogate MTase function. We develop a new set of peptides based on conserved regions of the nsp10 protein in the Coronaviridae species and test these to known MTase variant values. This results in the prediction that the H80R variant is a solid new candidate for potential new testing. We envision that this new lead is the beginning of a reputable foundation of a new computational method that combats coronaviruses and that is beneficial for new peptide drug development.

Network Pharmacology-Based Strategy and Molecular Docking to Explore the Potential Mechanism of Jintiange Capsule for Treating Osteoporosis

Background

With the advent of ageing population, osteoporosis (OP) has already become a global challenge. Jintiange capsule is extensively applied to treat OP in China. Although recent studies demonstrate that it generates significant effects on strengthening bone, the exact mechanism of the jintiange capsule for treating OP remains unknown.

Purpose

To understand the main ingredients of the jintiange capsule, predict the possible targets and the relevant signal transduction pathways, and explore the mechanism of the jintiange capsule for the treatment of OP.

Methods

Main ingredients of the jintiange capsule, drug targets, and potential disease targets for OP were obtained from public databases. Molecular biological processes and signaling pathways were determined via bioinformatic analysis, containing protein-protein interaction (PPI), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG). Subsequently, the disease-drug-ingredient-targets-pathways networks were constructed using Cytoscape. According to CytoNCA, core targets were acquired. Finally, the present study conducted molecular docking for better testing the abovementioned results.

Results

In the current work, we found that 4 main ingredients of the jintiange capsule, 33 drug targets, 4745 potential disease targets for OP, and 12 overlapping targets were identified. PPI network containing 12 nodes and 25 edges proved that there existed a complex relationship. As revealed by GO functional annotation, the intersected targets were mostly associated with BP, CC, and MF. The targets were enriched to 368 items in BP, 27 items in CC, and 42 items in MF. They mainly included calcium ion homeostasis, calcium channel complex, and calcium channel regulator activity. According to KEGG pathway analysis, the intersected targets were mostly associated with Rap 1, cGMP-PKG, Ras, cAMP, calcium pathways, and so on. Based on the analysis with CytoNCA, we acquired 4 core targets, respectively—CALR, SPARC, CALM1, and CALM2. Besides, 2 core targets, CALR and CALM1, were selected for molecular docking experiments. Molecular docking revealed that the main ingredient, calcium phosphate, had good binding with the CALR protein and CALM1 protein.

Conclusion

To conclude, the main ingredient of the jintiange capsule, particularly calcium phosphate, may interact with 2 targets, CALR and CALM1, and regulate multiple signaling pathways to treat OP. Additionally, this also benefits us in further understanding the mechanism of the jintiange capsule for treating OP.

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.

Revealing the Molecular Mechanisms of Alzheimer's Disease Based on Network Analysis

The complex pathology of Alzheimer's disease (AD) emphasises the need for comprehensive modelling of the disease, which may lead to the development of efficient treatment strategies. To address this challenge, we analysed transcriptome data of post-mortem human brain samples of healthy elders and individuals with late-onset AD from the Religious Orders Study and Rush Memory and Aging Project (ROSMAP) and Mayo Clinic (MayoRNAseq) studies in the AMP-AD consortium. In this context, we conducted several bioinformatics and systems medicine analyses including the construction of AD-specific co-expression networks and genome-scale metabolic modelling of the brain in AD patients to identify key genes, metabolites and pathways involved in the progression of AD. We identified AMIGO1 and GRPRASP2 as examples of commonly altered marker genes in AD patients. Moreover, we found alterations in energy metabolism, represented by reduced oxidative phosphorylation and ATPase activity, as well as the depletion of hexanoyl-CoA, pentanoyl-CoA, (2E)-hexenoyl-CoA and numerous other unsaturated fatty acids in the brain. We also observed that neuroprotective metabolites (e.g., vitamins, retinoids and unsaturated fatty acids) tend to be depleted in the AD brain, while neurotoxic metabolites (e.g., β-alanine, bilirubin) were more abundant. In summary, we systematically revealed the key genes and pathways related to the progression of AD, gained insight into the crucial mechanisms of AD and identified some possible targets that could be used in the treatment of AD.

Allosteric Switching of Calmodulin in a Mycobacterium smegmatis porin A (MspA) Nanopore-Trap

Recent developments concerning large protein nanopores suggest a new approach to structure profiling of native folded proteins. In this work, the large vestibule of Mycobacterium smegmatis porin A (MspA) and calmodulin (CaM), a Ca²⁺ -binding protein, were used in the direct observation of the protein structure. Three conformers, including the Ca²⁺ -free, Ca²⁺ -bound, and target peptide-bound states of CaM, were unambiguously distinguished. A disease related mutant, CaM D129G was also discriminated by MspA, revealing how a single amino acid replacement can interfere with the Ca²⁺ -binding capacity of the whole protein. The binding capacity and aggregation effect of CaM induced by different ions (Mg²⁺ /Sr²⁺ /Ba²⁺ /Ca²⁺ /Pb²⁺ /Tb³⁺ ) were also investigated and the stability of MspA in extreme conditions was evaluated. This work demonstrates the most systematic single-molecule investigation of different allosteric conformers of CaM, acknowledging the high sensing resolution offered by the MspA nanopore trap.

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.

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.

Network Pharmacology-Based Systematic Analysis of Molecular Mechanisms of Dingji Fumai Decoction for Ventricular Arrhythmia

Background

Dingji Fumai Decoction (DFD), a traditional herbal mixture, has been widely used to ventricular arrhythmia (VA) in clinical practice in China. However, research on the bioactive components and underlying mechanisms of DFD in VA is still scarce.

Methods

Components of DFD were collected from TCMSP, ETCM, and literature. The chemical structures of each component were obtained from PubChem. Next, SwissADME and SwissTargetPrediction were applied for compounds screening and targets prediction of DFD; meanwhile, targets of VA were collected from DrugBank and Online Mendelian Inheritance in Man (OMIM). Then, the H-C-T-D network and the protein-protein interaction (PPI) network were constructed based on the data obtained above. CytoNCA was utilized to filter hub genes and VarElect was used to analyze the relationship between genes and diseases. At last, Metascape was employed for systematic analysis on the potential targets of herbals against VA, and AutoDock was applied for molecular docking to verify the results.

Results

A total of 434 components were collected, 168 of which were qualified, and there were 28 shared targets between DFD and VA. Three function modules of DFD were found from the PPI network. Further systematic analysis of shared genes and function modules explained the potential mechanism of DFD in the treatment of VA; molecular docking has verified the interactions.

Conclusions

DFD could be employed for VA through mechanisms, including complex interactions between related components and targets, as predicted by network pharmacology and molecular docking. This work confirmed that DFD could apply to the treatment of VA and promoted the explanation of DFD for VA in the molecular mechanisms.

Predicting Genetic Variation Severity Using Machine Learning to Interpret Molecular Simulations

Distinct missense mutations in a specific gene have been associated with different diseases as well as differing severity of a disease. Current computational methods predict the potential pathogenicity of a missense variant but fail to differentiate between separate disease or severity phenotypes. We have developed a method to overcome this limitation by applying machine learning to features extracted from molecular dynamics simulations, creating a way to predict the effect of novel genetic variants in causing a disease, drug resistance, or another specific trait. As an example, we have applied this novel approach to variants in calmodulin associated with two distinct arrhythmias as well as two different neurodegenerative diseases caused by variants in amyloid-β peptide. The new method successfully predicts the specific disease caused by a gene variant and ranks its severity with more accuracy than existing methods. We call this method molecular dynamics phenotype prediction model.

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.

MicroRNA-129-1-3p protects cardiomyocytes from pirarubicin-induced apoptosis by down-regulating the GRIN2D-mediated Ca2+ signalling pathway

Pirarubicin (THP), an anthracycline anticancer drug, is a first‐line therapy for various solid tumours and haematologic malignancies. However, THP can cause dose‐dependent cumulative cardiac damage, which limits its therapeutic window. The mechanisms underlying THP cardiotoxicity are not fully understood. We previously showed that MiR‐129‐1‐3p, a potential biomarker of cardiovascular disease, was down‐regulated in a rat model of THP‐induced cardiac injury. In this study, we used Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway enrichment analyses to determine the pathways affected by miR‐129‐1‐3p expression. The results linked miR‐129‐1‐3p to the Ca²⁺ signalling pathway. TargetScan database screening identified a tentative miR‐129‐1‐3p‐binding site at the 3′‐UTR of GRIN2D, a subunit of the N‐methyl‐D‐aspartate receptor calcium channel. A luciferase reporter assay confirmed that miR‐129‐1‐3p directly regulates GRIN2D. In H9C2 (rat) and HL‐1 (mouse) cardiomyocytes, THP caused oxidative stress, calcium overload and apoptotic cell death. These THP‐induced changes were ameliorated by miR‐129‐1‐3p overexpression, but exacerbated by miR‐129‐1‐3p knock‐down. In addition, miR‐129‐1‐3p overexpression in cardiomyocytes prevented THP‐induced changes in the expression of proteins that are either key components of Ca²⁺ signalling or important regulators of intracellular calcium trafficking/balance in cardiomyocytes including GRIN2D, CALM1, CaMKⅡδ, RyR2‐pS2814, SERCA2a and NCX1. Together, these bioinformatics and cell‐based experiments indicate that miR‐129‐1‐3p protects against THP‐induced cardiomyocyte apoptosis by down‐regulating the GRIN2D‐mediated Ca²⁺ pathway. Our results reveal a novel mechanism underlying the pathogenesis of THP‐induced cardiotoxicity. The miR‐129‐1‐3p/Ca²⁺ signalling pathway could serve as a target for the development of new cardioprotective agents to control THP‐induced cardiotoxicity.