Molecular modeling of disease causing mutations in domain C1 of cMyBP-C

Functional Characterization of the MYO6 Variant p.E60Q in Non-Syndromic Hearing Loss Patients

Hereditary hearing loss (HHL) is a common genetic disorder accounting for at least 60% of pre-lingual deafness in children, of which 70% is inherited in an autosomal recessive pattern. The long tradition of consanguinity among the Qatari population has increased the prevalence of HHL, which negatively impacts the quality of life. Here, we functionally validated the pathogenicity of the c.178G>C, p.E60Q mutation in the MYO6 gene, which was detected previously in a Qatari HHL family, using cellular and animal models. In vitro analysis was conducted in HeLa cells transiently transfected with plasmids carrying MYO6WT or MYO6p.E60Q, and a zebrafish model was generated to characterize the in vivo phenotype. Cells transfected with MYO6WT showed higher expression of MYO6 in the plasma membrane and increased ATPase activity. Modeling the human MYO6 variants in zebrafish resulted in severe otic defects. At 72 h post-injection, MYO6p.E60Q embryos demonstrated alterations in the sizes of the saccule and utricle. Additionally, zebrafish with MYO6p.E60Q displayed super-coiled and bent hair bundles in otic hair cells when compared to control and MYO6WT embryos. In conclusion, our cellular and animal models add support to the in silico prediction that the p.E60Q missense variant is pathogenic and damaging to the protein. Since the c.178G>C MYO6 variant has a 0.5% allele frequency in the Qatari population, about 400 times higher than in other populations, it could contribute to explaining the high prevalence of hearing impairment in Qatar.

The mechanics of the heart: zooming in on hypertrophic cardiomyopathy and cMyBP-C

Hypertrophic cardiomyopathy (HCM), a disease characterized by cardiac muscle hypertrophy and hypercontractility, is the most frequently inherited disorder of the heart. HCM is mainly caused by variants in genes encoding proteins of the sarcomere, the basic contractile unit of cardiomyocytes. The most frequently mutated among them is MYBPC3, which encodes cardiac myosin-binding protein C (cMyBP-C), a key regulator of sarcomere contraction. In this review, we summarize clinical and genetic aspects of HCM and provide updated information on the function of the healthy and HCM sarcomere, as well as on emerging therapeutic options targeting sarcomere mechanical activity. Building on what is known about cMyBP-C activity, we examine different pathogenicity drivers by which MYBPC3 variants can cause disease, focussing on protein haploinsufficiency as a common pathomechanism also in nontruncating variants. Finally, we discuss recent evidence correlating altered cMyBP-C mechanical properties with HCM development.

Clinical, Genetic and Functional Characterization of a Novel AVPR2 Missense Mutation in a Woman with X-Linked Recessive Nephrogenic Diabetes Insipidus

Nephrogenic diabetes insipidus (NDI) is a rare disorder characterized by renal unresponsiveness to the hormone vasopressin, leading to excretion of large volumes of diluted urine. Mutations in the arginine vasopressin receptor-2 (AVPR2) gene cause congenital NDI and have an X-linked recessive inheritance. The disorder affects almost exclusively male family members, but female carriers occasionally present partial phenotypes due to skewed inactivation of the X-chromosome. Here, we report a rare case of a woman affected with X-linked recessive NDI, presenting an average urinary output of 12 L/day. Clinical and biochemical studies showed incomplete responses to water deprivation and vasopressin stimulation tests. Genetic analyses revealed a novel heterozygous missense mutation (c.493G > C, p.Ala165Pro) in the AVPR2 gene. Using a combination of in-silico protein modeling with human cellular models and molecular phenotyping, we provide functional evidence for phenotypic effects. The mutation destabilizes the helical structure of the AVPR2 transmembrane domains and disrupts its plasma membrane localization and downstream intracellular signaling pathways upon activation with its agonist vasopressin. These defects lead to deficient aquaporin 2 (AQP2) membrane translocation, explaining the inability to concentrate urine in this patient.

Functional characterization of human myosin-binding protein C3 variants associated with hypertrophic cardiomyopathy reveals exon-specific cardiac phenotypes in zebrafish model

Myosin-binding protein C 3 (MYBPC3) variants are the most common cause of hypertrophic cardiomyopathy (HCM). HCM is a complex cardiac disorder due to its significant genetic and clinical heterogeneity. MYBPC3 variants genotype-phenotype associations remain poorly understood. We investigated the impact of two novel human MYBPC3 splice-site variants: V1: c.654+2_654+4dupTGG targeting exon 5 using morpholino MOe5i5; and V2: c.772+1G>A targeting exon 6 using MOe6i6; located within C1 domain of cMyBP-C protein, known to be critical in regulating sarcomere structure and contractility. Zebrafish MOe5i5 and MOe6i6 morphants recapitulated typical characteristics of human HCM with cardiac phenotypes of varying severity, including reduced cardiomyocyte count, thickened ventricular myocardial wall, a drastic reduction in heart rate, stroke volume, and cardiac output. Analysis of all cardiac morphological and functional parameters demonstrated that V2 cardiac phenotype was more severe than V1. Coinjection with synthetic human MYBPC3 messenger RNA (mRNA) partially rescued disparate cardiac phenotypes in each zebrafish morphant. While human MYBPC3 mRNA partially restored the decreased heart rate in V1 morphants and displayed increased percentages of ejection fraction, fractional shortening, and area change, it failed to revert the V1 ventricular myocardial thickness. These results suggest a possible V1 impact on cardiac contractility. In contrast, attempts to rescue V2 morphants only restored the ventricular myocardial wall hypertrophy phenotype but had no significant effect on impaired heart rate, suggesting a potential V2 impact on the cardiac structure. Our study provides evidence of an association between MYBPC3 exon-specific cardiac phenotypes in the zebrafish model providing important insights into how these genetic variants contribute to HCM disease.

Multiscale Models of Cardiac Muscle Biophysics and Tissue Remodeling in Hypertrophic Cardiomyopathies

Myocardial hypertrophy is the result of sustained perturbations to the mechanical and/or neurohormonal homeostasis of cardiac cells and is driven by integrated, multiscale biophysical and biochemical processes that are currently not well defined. In this brief review, we highlight recent computational and experimental models of cardiac hypertrophy that span mechanisms from the molecular level to the tissue level. Specifically, we focus on: (i) molecular-level models of the structural dynamics of sarcomere proteins in hypertrophic hearts, (ii) cellular-level models of excitation-contraction coupling and mechanosensitive signaling in disease-state myocytes, and (iii) organ-level models of myocardial growth kinematics and predictors thereof. Finally, we discuss how spanning these scales and combining multiple experimental/computational models will provide new information about the processes governing hypertrophy and potential methods to prevent or reverse them.

In silico and in vivo models for Qatari-specific classical homocystinuria as basis for development of novel therapies

Homocystinuria is a rare inborn error of methionine metabolism caused by cystathionine β-synthase (CBS) deficiency. The prevalence of homocystinuria in Qatar is 1:1,800 births, mainly due to a founder Qatari missense mutation, c.1006C>T; p.R336C (p.Arg336Cys). We characterized the structure-function relationship of the p.R336C-mutant protein and investigated the effect of different chemical chaperones to restore p.R336C-CBS activity using three models: in silico, ΔCBS yeast, and CRISPR/Cas9 p.R336C knock-in HEK293T and HepG2 cell lines. Protein modeling suggested that the p.R336C induces severe conformational and structural changes, perhaps influencing CBS activity. Wild-type CBS, but not the p.R336C mutant, was able to restore the yeast growth in ΔCBS-deficient yeast in a complementation assay. The p.R336C knock-in HEK293T and HepG2 cells decreased the level of CBS expression and reduced its structural stability; however, treatment of the p.R336C knock-in HEK293T cells with betaine, a chemical chaperone, restored the stability and tetrameric conformation of CBS, but not its activity. Collectively, these results indicate that the p.R336C mutation has a deleterious effect on CBS structure, stability, and activity, and using the chemical chaperones approach for treatment could be ineffective in restoring p.R336C CBS activity.

Hypertrophic cardiomyopathy-linked variants of cardiac myosin-binding protein C3 display altered molecular properties and actin interaction

The most common inherited cardiac disorder, hypertrophic cardiomyopathy (HCM), is characterized by thickening of heart muscle, for which genetic mutations in cardiac myosin-binding protein C3 (c-MYBPC3) gene, is the leading cause. Notably, patients with HCM display a heterogeneous clinical presentation, onset and prognosis. Thus, delineating the molecular mechanisms that explain how disparate c-MYBPC3 variants lead to HCM is essential for correlating the impact of specific genotypes on clinical severity. Herein, five c-MYBPC3 missense variants clinically associated with HCM were investigated; namely V1 (R177H), V2 (A216T), V3 (E258K), V4 (E441K) and double mutation V5 (V3 + V4), all located within the C1 and C2 domains of MyBP-C, a region known to interact with sarcomeric protein, actin. Injection of the variant complementary RNAs in zebrafish embryos was observed to recapitulate phenotypic aspects of HCM in patients. Interestingly, V3- and V5-cRNA injection produced the most severe zebrafish cardiac phenotype, exhibiting increased diastolic/systolic myocardial thickness and significantly reduced heart rate compared with control zebrafish. Molecular analysis of recombinant C0-C2 protein fragments revealed that c-MYBPC3 variants alter the C0-C2 domain secondary structure, thermodynamic stability and importantly, result in a reduced binding affinity to cardiac actin. V5 (double mutant), displayed the greatest protein instability with concomitant loss of actin-binding function. Our study provides specific mechanistic insight into how c-MYBPC3 pathogenic variants alter both functional and structural characteristics of C0-C2 domains leading to impaired actin interaction and reduced contractility, which may provide a basis for elucidating the disease mechanism in HCM patients with c- MYBPC3 mutations.

N-Terminal Domains of Cardiac Myosin Binding Protein C Cooperatively Activate the Thin Filament

Muscle contraction relies on interaction between myosin-based thick filaments and actin-based thin filaments. Myosin binding protein C (MyBP-C) is a key regulator of actomyosin interactions. Recent studies established that the N'-terminal domains (NTDs) of MyBP-C can either activate or inhibit thin filaments, but the mechanism of their collective action is poorly understood. Cardiac MyBP-C (cMyBP-C) harbors an extra NTD, which is absent in skeletal isoforms of MyBP-C, and its role in regulation of cardiac contraction is unknown. Here we show that the first two domains of human cMyPB-C (i.e., C0 and C1) cooperate to activate the thin filament. We demonstrate that C1 interacts with tropomyosin via a positively charged loop and that this interaction, stabilized by the C0 domain, is required for thin filament activation by cMyBP-C. Our data reveal a mechanism by which cMyBP-C can modulate cardiac contraction and demonstrate a function of the C0 domain.

Next-generation sequencing identified SPATC1L as a possible candidate gene for both early-onset and age-related hearing loss

Hereditary hearing loss (HHL) and age-related hearing loss (ARHL) are two major sensory diseases affecting millions of people worldwide. Despite many efforts, additional HHL-genes and ARHL genetic risk factors still need to be identified. To fill this gap a large genomic screening based on next-generation sequencing technologies was performed. Whole exome sequencing in a 3-generation Italian HHL family and targeted re-sequencing in 464 ARHL patients were performed. We detected three variants in SPATC1L: a nonsense allele in an HHL family and a frameshift insertion and a missense variation in two unrelated ARHL patients. In silico molecular modelling of all variants suggested a significant impact on the structural stability of the protein itself, likely leading to deleterious effects and resulting in truncated isoforms. After demonstrating Spatc1l expression in mice inner ear, in vitro functional experiments were performed confirming the results of the molecular modelling studies. Finally, a candidate-gene population-based statistical study in cohorts from Caucasus and Central Asia revealed a statistically significant association of SPATC1L with normal hearing function at low and medium hearing frequencies. Overall, the amount of different genetic data presented here (variants with early-onset and late-onset hearing loss in addition to genetic association with normal hearing function), together with relevant functional evidence, likely suggest a role of SPATC1L in hearing function and loss.

Targeted sequencing identifies novel variants involved in autosomal recessive hereditary hearing loss in Qatari families

Hereditary hearing loss is characterized by a very high genetic heterogeneity. In the Qatari population the role of GJB2, the worldwide HHL major player, seems to be quite limited compared to Caucasian populations. In this study we analysed 18 Qatari families affected by non-syndromic hearing loss using a targeted sequencing approach that allowed us to analyse 81 genes simultaneously. Thanks to this approach, 50% of these families (9 out of 18) resulted positive for the presence of likely causative alleles in 6 different genes: CDH23, MYO6, GJB6, OTOF, TMC1 and OTOA. In particular, 4 novel alleles were detected while the remaining ones were already described to be associated to HHL in other ethnic groups. Molecular modelling has been used to further investigate the role of novel alleles identified in CDH23 and TMC1 genes demonstrating their crucial role in Ca2+ binding and therefore possible functional role in proteins. Present study showed that an accurate molecular diagnosis based on next generation sequencing technologies might largely improve molecular diagnostics outcome leading to benefits for both genetic counseling and definition of recurrence risk.

The Role of Cardiac Myosin Binding Protein C3 in Hypertrophic Cardiomyopathy-Progress and Novel Therapeutic Opportunities

Hypertrophic cardiomyopathy (HCM) is a common autosomal dominant genetic cardiovascular disorder marked by genetic and phenotypic heterogeneity. Mutations in the gene encodes the cardiac myosin-binding protein C, cMYBPC3 is amongst the various sarcomeric genes that are associated with HCM. These mutations produce mutated mRNAs and truncated cMyBP-C proteins. In this review, we will discuss the implications and molecular mechanisms involved in MYBPC3 different mutations. Further, we will highlight the novel targets that can be developed into potential therapeutics for the treatment of HMC. J. Cell. Physiol. 232: 1650-1659, 2017. © 2016 Wiley Periodicals, Inc.

Two patients with Canavan disease and structural modeling of a novel mutation

Canavan disease (CD) is a rare fatal childhood neurological autosomal recessive genetic disease caused by mutations in the ASPA gene, which lead to catalytic deficiency of the ASPA enzyme, which catalyzes the hydrolysis of N-acetyl-L-aspartate (NAA) into aspartate and acetate. CD occurs frequently among Ashkenazi Jewish population, however it has been reported in many other ethnic groups with significantly lower frequency. Here, we report on two Egyptian patients diagnosed with CD, the first patient harbors five missense mutations (c.427 A > G; p. I143V, c.502C > T; p. R168C, c.530 T > C; p. I177T, c.557 T > C; p. V186D c.548C > T; p. P183L) and a silent mutation (c.693 C > T; p. Y231Y). The second patient was found to be homozygous for two missense mutations (c.427 A > G; p. I143V and c.557 T > A; p. V186D). Furthermore, molecular modeling of the novel mutation p. P183L provides an instructive explanation of the mutational impact on the protein structure that can affect the function of the ASPA. Here, the clinical, radiological, and biochemical profile of the two patients are reviewed in details.

Structural modeling of p.V31F variant in the aspartoacylase gene

Aspartoacylase (ASPA) is an abundant enzyme in the brain, which catalyzes the conversion of N-acetylaspartate into acetate and aspartate, deficiency in its activity leads to degeneration of the white matter of the brain and is a recognized cause of Canavan disease (CD), which affect children. Although genotype-phenotype correlation have been reported for Canavan disease patients, this relationships is still not straightforward. In this communication, we use molecular modeling to address the structural consequences resulting from the missense variant p.V31F in the ASPA enzyme, which we previously reported in a homozygous form in an Egyptian patient with infantile CD. This modeling suggests that this variant brings significant changes to the catalytic core by introducing structural flexibility through neighbouring key residues. In particular, it provides a molecular explanation for the pathogenic effect of this variant and provides a meaningful genotype-phonotype relationships. The mutational impact appears to have an influence on the function of the protein and initiates molecular event for the mechanism of the disease.

Genetics of hypertrophic cardiomyopathy: advances and pitfalls in molecular diagnosis and therapy

Hypertrophic cardiomyopathy (HCM) is a primary disease of the cardiac muscle that occurs mainly due to mutations (>1,400 variants) in genes encoding for the cardiac sarcomere. HCM, the most common familial form of cardiomyopathy, affecting one in every 500 people in the general population, is typically inherited in an autosomal dominant pattern, and presents variable expressivity and age-related penetrance. Due to the morphological and pathological heterogeneity of the disease, the appearance and progression of symptoms is not straightforward. Most HCM patients are asymptomatic, but up to 25% develop significant symptoms, including chest pain and sudden cardiac death. Sudden cardiac death is a dramatic event, since it occurs without warning and mainly in younger people, including trained athletes. Molecular diagnosis of HCM is of the outmost importance, since it may allow detection of subjects carrying mutations on HCM-associated genes before development of clinical symptoms of HCM. However, due to the genetic heterogeneity of HCM, molecular diagnosis is difficult. Currently, there are mainly four techniques used for molecular diagnosis of HCM, including Sanger sequencing, high resolution melting, mutation detection using DNA arrays, and next-generation sequencing techniques. Application of these methods has proven successful for identification of mutations on HCM-related genes. This review summarizes the features of these technologies, highlighting their strengths and weaknesses. Furthermore, current therapeutics for HCM patients are correlated with clinically observed phenotypes and are based on the alleviation of symptoms. This is mainly due to insufficient knowledge on the mechanisms involved in the onset of HCM. Tissue engineering alongside regenerative medicine coupled with nanotherapeutics may allow fulfillment of those gaps, together with screening of novel therapeutic drugs and target delivery systems.

Decade in review--cardiomyopathies: Cardiomyopathy on the move

Since Wallace Brigden first used the term 'cardiomyopathy' in 1952, this group of diseases has continued to attract the interest of clinicians, researchers, and importantly, patients. The past decade has seen a substantial accumulation of knowledge relating to various cardiomyopathies, which has partially lifted the mystery surrounding this topic.

Insights from molecular dynamics simulations: structural basis for the V567D mutation-induced instability of zebrafish alpha-dystroglycan and comparison with the murine model

A missense amino acid mutation of valine to aspartic acid in 567 position of alpha-dystroglycan (DG), identified in dag1-mutated zebrafish, results in a reduced transcription and a complete absence of the protein. Lacking experimental structural data for zebrafish DG domains, the detailed mechanism for the observed mutation-induced destabilization of the DG complex and membrane damage, remained unclear. With the aim to contribute to a better clarification of the structure-function relationships featuring the DG complex, three-dimensional structural models of wild-type and mutant (V567D) C-terminal domain of alpha-DG from zebrafish were constructed by a template-based modelling approach. We then ran extensive molecular dynamics (MD) simulations to reveal the structural and dynamic properties of the C-terminal domain and to evaluate the effect of the single mutation on alpha-DG stability. A comparative study has been also carried out on our previously generated model of murine alpha-DG C-terminal domain including the I591D mutation, which is topologically equivalent to the V567D mutation found in zebrafish. Trajectories from MD simulations were analyzed in detail, revealing extensive structural disorder involving multiple beta-strands in the mutated variant of the zebrafish protein whereas local effects have been detected in the murine protein. A biochemical analysis of the murine alpha-DG mutant I591D confirmed a pronounced instability of the protein. Taken together, the computational and biochemical analysis suggest that the V567D/I591D mutation, belonging to the G beta-strand, plays a key role in inducing a destabilization of the alpha-DG C-terminal Ig-like domain that could possibly affect and propagate to the entire DG complex. The structural features herein identified may be of crucial help to understand the molecular basis of primary dystroglycanopathies.

Cardiac myosin binding protein-C: a structurally dynamic regulator of myocardial contractility

Cardiac myosin binding protein-C (cMyBP-C) is a modular protein anchored to the thick filament through interactions mediated by its C-terminal region. The N-terminal region of cMyBPC-C regulates myocardial contractility by modifying actin-myosin association. Phosphorylation of the N-terminal region diminishes cMyBP-C's capacity to regulate actin-myosin function. Despite a substantial body of literature, many issues remain unclear regarding the structural and functional roles of cMyBP-C. While no high-resolution structures of the intact protein exist, crystallographic and nuclear magnetic resonance (NMR) structures of isolated N-terminal domains provide important molecular details regarding cMyBP-C's role in controlling contractility. In this review, we summarize the emerging structural understanding of cMyBP-C with a particular emphasis placed on describing how its dynamic molecular interactions with both thin and thick filament proteins likely contribute to contractile regulation. Furthermore, we discuss the future directions and strategies by which we may improve the mechanistic understanding of its role in modulating cardiac muscle contraction.