Computational prediction of protein subdomain stability in MYBPC3 enables clinical risk stratification in hypertrophic cardiomyopathy and enhances variant interpretation

PURPOSE: Variants in MYBPC3 causing loss of function are the most common cause of hypertrophic cardiomyopathy (HCM). However, a substantial number of patients carry missense variants of uncertain significance (VUS) in MYBPC3. We hypothesize that a structural-based algorithm, STRUM, which estimates the effect of missense variants on protein folding, will identify a subgroup of HCM patients with a MYBPC3 VUS associated with increased clinical risk. METHODS: Among 7,963 patients in the multicenter Sarcomeric Human Cardiomyopathy Registry (SHaRe), 120 unique missense VUS in MYBPC3 were identified. Variants were evaluated for their effect on subdomain folding and a stratified time-to-event analysis for an overall composite endpoint (first occurrence of ventricular arrhythmia, heart failure, all-cause mortality, atrial fibrillation, and stroke) was performed for patients with HCM and a MYBPC3 missense VUS. RESULTS: We demonstrated that patients carrying a MYBPC3 VUS predicted to cause subdomain misfolding (STRUM+, ΔΔG ≤ −1.2 kcal/mol) exhibited a higher rate of adverse events compared with those with a STRUM- VUS (hazard ratio = 2.29, P = 0.0282). In silico saturation mutagenesis of MYBPC3 identified 4,943/23,427 (21%) missense variants that were predicted to cause subdomain misfolding. CONCLUSION: STRUM identifies patients with HCM and a MYBPC3 VUS who may be at higher clinical risk and provides supportive evidence for pathogenicity

Spatial and Functional Distribution of MYBPC3 Pathogenic Variants and Clinical Outcomes in Patients with Hypertrophic Cardiomyopathy

Background - Pathogenic variants in MYBPC3, encoding cardiac MyBP-C, are the most common cause of familial hypertrophic cardiomyopathy. A large number of unique MYBPC3 variants and relatively small genotyped HCM cohorts have precluded detailed genotype-phenotype correlations. Methods - Patients with HCM and MYBPC3 variants were identified from the Sarcomeric Human Cardiomyopathy Registry (SHaRe). Variant types and locations were analyzed, morphologic severity was assessed, and time-event analysis was performed (composite clinical outcome of sudden death, class III/IV heart failure, LVAD/transplant, atrial fibrillation). For selected missense variants falling in enriched domains, myofilament localization and degradation rates were measured in vitro. Results - Among 4,756 genotyped HCM patients in SHaRe, 1,316 patients were identified with adjudicated pathogenic truncating (N=234 unique variants, 1047 patients) or non-truncating (N=22 unique variants, 191 patients) variants in MYBPC3. Truncating variants were evenly dispersed throughout the gene, and hypertrophy severity and outcomes were not associated with variant location (grouped by 5' - 3' quartiles or by founder variant subgroup). Non-truncating pathogenic variants clustered in the C3, C6, and C10 domains (18 of 22, 82%, p<0.001 vs. gnomAD common variants) and were associated with similar hypertrophy severity and adverse event rates as observed with truncating variants. MyBP-C with variants in the C3, C6, and C10 domains was expressed in rat ventricular myocytes. C10 mutant MyBP-C failed to incorporate into myofilaments and degradation rates were accelerated by ~90%, while C3 and C6 mutant MyBP-C incorporated normally with degradation rate similar to wild-type. Conclusions - Truncating variants account for 91% of MYBPC3 pathogenic variants and cause similar clinical severity and outcomes regardless of location, consistent with locus-independent loss-of-function. Non-truncating MYBPC3 pathogenic variants are regionally clustered, and a subset also cause loss-of-function through failure of myofilament incorporation and rapid degradation. Cardiac morphology and clinical outcomes are similar in patients with truncating vs. non-truncating variants

An Unbiased Screen Identified the Hsp70-BAG3 Complex as a Regulator of Myosin-Binding Protein C3.

Variants in the gene myosin-binding protein C3 (MYBPC3) account for approximately 50% of familial hypertrophic cardiomyopathy (HCM), leading to reduced levels of myosin-binding protein C3 (MyBP-C), the protein product made by gene MYBPC3. Elucidation of the pathways that regulate MyBP-C protein homeostasis could uncover new therapeutic strategies. Toward this goal, we screened a library of 2,426 bioactive compounds and identified JG98, an allosteric modulator of heat shock protein 70 that inhibits interaction with Bcl-2-associated athanogene (BAG) domain co-chaperones. JG98 reduces MyBP-C protein levels. Furthermore, genetic reduction of BAG3 phenocopies treatment with JG-98 by reducing MYBP-C protein levels.. Thus, an unbiased compound screen identified the heat shock protein 70-BAG3 complex as a regulator of MyBP-C stability

HSC70 is a chaperone for wild-type and mutant cardiac myosin binding protein C

Cardiac myosin binding protein C (MYBPC3) is the most commonly mutated gene associated with hypertrophic cardiomyopathy (HCM). Haploinsufficiency of full-length MYBPC3 and disruption of proteostasis have both been proposed as central to HCM disease pathogenesis. Discriminating the relative contributions of these 2 mechanisms requires fundamental knowledge of how turnover of WT and mutant MYBPC3 proteins is regulated. We expressed several disease-causing mutations in MYBPC3 in primary neonatal rat ventricular cardiomyocytes. In contrast to WT MYBPC3, mutant proteins showed reduced expression and failed to localize to the sarcomere. In an unbiased coimmunoprecipitation/mass spectrometry screen, we identified HSP70-family chaperones as interactors of both WT and mutant MYBPC3. Heat shock cognate 70 kDa (HSC70) was the most abundant chaperone interactor. Knockdown of HSC70 significantly slowed degradation of both WT and mutant MYBPC3, while pharmacologic activation of HSC70 and HSP70 accelerated degradation. HSC70 was expressed in discrete striations in the sarcomere. Expression of mutant MYBPC3 did not affect HSC70 localization, nor did it induce a protein folding stress response or ubiquitin proteasome dysfunction. Together these data suggest that WT and mutant MYBPC3 proteins are clients for HSC70, and that the HSC70 chaperone system plays a major role in regulating MYBPC3 protein turnover