The Relationship between Changes in MYBPC3 Single-Nucleotide Polymorphism-Associated Metabolites and Elite Athletes' Adaptive Cardiac Function

Athletic performance is a multifactorial trait influenced by a complex interaction of environmental and genetic factors. Over the last decades, understanding and improving elite athletes' endurance and performance has become a real challenge for scientists. Significant tools include but are not limited to the development of molecular methods for talent identification, personalized exercise training, dietary requirements, prevention of exercise-related diseases, as well as the recognition of the structure and function of the genome in elite athletes. Investigating the genetic markers and phenotypes has become critical for elite endurance surveillance. The identification of genetic variants contributing to a predisposition for excellence in certain types of athletic activities has been difficult despite the relatively high genetic inheritance of athlete status. Metabolomics can potentially represent a useful approach for gaining a thorough understanding of various physiological states and for clarifying disorders caused by strength-endurance physical exercise. Based on a previous GWAS study, this manuscript aims to discuss the association of specific single-nucleotide polymorphisms (SNPs) located in the MYBPC3 gene encoding for cardiac MyBP-C protein with endurance athlete status. MYBPC3 is linked to elite athlete heart remodeling during or after exercise, but it could also be linked to the phenotype of cardiac hypertrophy (HCM). To make the distinction between both phenotypes, specific metabolites that are influenced by variants in the MYBPC3 gene are analyzed in relation to elite athletic performance and HCM. These include theophylline, ursodeoxycholate, quinate, and decanoyl-carnitine. According to the analysis of effect size, theophylline, quinate, and decanoyl carnitine increase with endurance while decreasing with cardiovascular disease, whereas ursodeoxycholate increases with cardiovascular disease. In conclusion, and based on our metabolomics data, the specific effects on athletic performance for each MYBPC3 SNP-associated metabolite are discussed.

Phosphorylation-dependent interactions of myosin-binding protein C and troponin coordinate the myofilament response to protein kinase A

PKA-mediated phosphorylation of sarcomeric proteins enhances heart muscle performance in response to β-adrenergic stimulation and is associated with accelerated relaxation and increased cardiac output for a given preload. At the cellular level, the latter translates to a greater dependence of Ca2+ sensitivity and maximum force on sarcomere length (SL), that is, enhanced length-dependent activation. However, the mechanisms by which PKA phosphorylation of the most notable sarcomeric PKA targets, troponin I (cTnI) and myosin-binding protein C (cMyBP-C), lead to these effects remain elusive. Here, we specifically altered the phosphorylation level of cTnI in heart muscle cells and characterized the structural and functional effects at different levels of background phosphorylation of cMyBP-C and with two different SLs. We found Ser22/23 bisphosphorylation of cTnI was indispensable for the enhancement of length-dependent activation by PKA, as was cMyBP-C phosphorylation. This high level of coordination between cTnI and cMyBP-C may suggest coupling between their regulatory mechanisms. Further evidence for this was provided by our finding that cardiac troponin (cTn) can directly interact with cMyBP-C in vitro, in a phosphorylation- and Ca2+-dependent manner. In addition, bisphosphorylation at Ser22/Ser23 increased Ca2+ sensitivity at long SL in the presence of endogenously phosphorylated cMyBP-C. When cMyBP-C was dephosphorylated, bisphosphorylation of cTnI increased Ca2+ sensitivity and decreased cooperativity at both SLs, which may translate to deleterious effects in physiological settings. Our results could have clinical relevance for disease pathways, where PKA phosphorylation of cTnI may be functionally uncoupled from cMyBP-C phosphorylation due to mutations or haploinsufficiency.

Cardiac Myosin Filaments are Maintained by Stochastic Protein Replacement

Myosin and myosin-binding protein C are exquisitely organized into giant filamentous macromolecular complexes within cardiac muscle sarcomeres, yet these proteins must be continually replaced to maintain contractile fidelity. The overall hypothesis that myosin filament structure is dynamic and allows for the stochastic replacement of individual components was tested in vivo, using a combination of mass spectrometry- and fluorescence-based proteomic techniques. Adult mice were fed a diet that marked all newly synthesized proteins with a stable isotope-labeled amino acid. The abundance of unlabeled and labeled proteins was quantified by high-resolution mass spectrometry over an 8-week period. The rates of change in the abundance of these proteins were well described by analytical models in which protein synthesis defined stoichiometry and protein degradation was governed by the stochastic selection of individual molecules. To test whether the whole myosin filaments or the individual components were selected for replacement, cardiac muscle was chemically skinned to remove the cellular membrane and myosin filaments were solubilized with ionic solutions. The composition of the filamentous and soluble fractions was quantified by mass spectrometry, and filament depolymerization was visualized by real-time fluorescence microscopy. Myosin molecules were preferentially extracted from ends of the filaments in the presence of the ionic solutions, and there was only a slight bias in the abundance of unlabeled molecules toward the innermost region on the myosin filaments. These data demonstrate for the first time that the newly synthesized myosin and myosin-binding protein C molecules are randomly mixed into preexisting thick filaments in vivo and the rate of mixing may not be equivalent along the length of the thick filament. These data collectively support a new model of cardiac myosin filament structure, with the filaments being dynamic macromolecular assemblies that allow for replacement of their components, rather than rigid bodies.

Microscale thermophoresis suggests a new model of regulation of cardiac myosin function via interaction with cardiac myosin-binding protein C

The cardiac isoform of myosin-binding protein C (cMyBP-C) is a key regulatory protein found in cardiac myofilaments that can control the activation state of both the actin-containing thin and myosin-containing thick filaments. However, in contrast to thin filament–based mechanisms of regulation, the mechanism of myosin-based regulation by cMyBP-C has yet to be defined in detail. To clarify its function in this process, we used microscale thermophoresis to build an extensive interaction map between cMyBP-C and isolated fragments of β-cardiac myosin. We show here that the regulatory N-terminal domains (C0C2) of cMyBP-C interact with both the myosin head (myosin S1) and tail domains (myosin S2) with micromolar affinity via phosphorylation-independent and phosphorylation-dependent interactions of domain C1 and the cardiac-specific m-motif, respectively. Moreover, we show that the interaction sites with the highest affinity between cMyBP-C and myosin S1 are localized to its central domains, which bind myosin with submicromolar affinity. We identified two separate interaction regions in the central C2C4 and C5C7 segments that compete for the same binding site on myosin S1, suggesting that cMyBP-C can crosslink the two myosin heads of a single myosin molecule and thereby stabilize it in the folded OFF state. Phosphorylation of the cardiac-specific m-motif by protein kinase A had no effect on the binding of either the N-terminal or the central segments to the myosin head domain, suggesting this might therefore represent a constitutively bound state of myosin associated with cMyBP-C. Based on our results, we propose a new model of regulation of cardiac myosin function by cMyBP-C.

Stress-dependent activation of myosin in the heart requires thin filament activation and thick filament mechanosensing

The efficiency of the heart as a pump depends on an autoregulatory mechanism, the Frank–Starling law of the heart, that potentiates the strength of contraction in response to an increase in ventricular filling. Disruption of this mechanism compromises the ability of the heart to pump blood, potentially leading to heart failure. We used fluorescent probes on myosin in heart muscle cells to investigate the molecular basis of the Frank–Starling mechanism. Our results show that the stronger contraction of heart muscle at longer lengths is due to a calcium-dependent interfilament signaling pathway that links stress sensing in the myosin-containing filaments with calcium activation of the actin-containing filaments. This pathway can potentially be targeted for treating heart failure.

Myosin-based regulation in the heart muscle modulates the number of myosin motors available for interaction with calcium-regulated thin filaments, but the signaling pathways mediating the stronger contraction triggered by stretch between heartbeats or by phosphorylation of the myosin regulatory light chain (RLC) remain unclear. Here, we used RLC probes in demembranated cardiac trabeculae to investigate the molecular structural basis of these regulatory pathways. We show that in relaxed trabeculae at near-physiological temperature and filament lattice spacing, the RLC-lobe orientations are consistent with a subset of myosin motors being folded onto the filament surface in the interacting-heads motif seen in isolated filaments. The folded conformation of myosin is disrupted by cooling relaxed trabeculae, similar to the effect induced by maximal calcium activation. Stretch or increased RLC phosphorylation in the physiological range have almost no effect on RLC conformation at a calcium concentration corresponding to that between beats. These results indicate that in near-physiological conditions, the folded myosin motors are not directly switched on by RLC phosphorylation or by the titin-based passive tension at longer sarcomere lengths in the absence of thin filament activation. However, at the higher calcium concentrations that activate the thin filaments, stretch produces a delayed activation of folded myosin motors and force increase that is potentiated by RLC phosphorylation. We conclude that the increased contractility of the heart induced by RLC phosphorylation and stretch can be explained by a calcium-dependent interfilament signaling pathway involving both thin filament sensitization and thick filament mechanosensing.

Myosin filament-based regulation of the dynamics of contraction in heart muscle

Cardiovascular disease continues to be the leading cause of death worldwide, and is frequently associated with heart failure. Efforts to develop better therapeutics for heart failure have been held back by limited understanding of the normal control of contraction on the timescale of the heartbeat. We used synchrotron X-ray diffraction to determine the dynamic structural changes in the myosin motors that drive contraction in the heart muscle, and show that myosin filament-based control mechanisms determine the time course and strength of contraction, allowing those mechanisms to be targeted for developing new therapies for heart disease.

Myosin-based mechanisms are increasingly recognized as supplementing their better-known actin-based counterparts to control the strength and time course of contraction in both skeletal and heart muscle. Here we use synchrotron small-angle X-ray diffraction to determine the structural dynamics of local domains of the myosin filament during contraction of heart muscle. We show that, although myosin motors throughout the filament contribute to force development, only about 10% of the motors in each filament bear the peak force, and these are confined to the filament domain containing myosin binding protein-C, the “C-zone.” Myosin motors in domains further from the filament midpoint are likely to be activated and inactivated first in each contraction. Inactivated myosin motors are folded against the filament core, and a subset of folded motors lie on the helical tracks described previously. These helically ordered motors are also likely to be confined to the C-zone, and the associated motor conformation reforms only slowly during relaxation. Myosin filament stress-sensing determines the strength and time course of contraction in conjunction with actin-based regulation. These results establish the fundamental roles of myosin filament domains and the associated motor conformations in controlling the strength and dynamics of contraction in heart muscle, enabling those structures to be targeted to develop new therapies for heart disease.

Noncoding RNAs versus Protein Biomarkers in Cardiovascular Disease

The development of more sensitive protein biomarker assays results in continuous improvements in detectability, extending the range of clinical applications to the detection of subclinical cardiovascular disease (CVD). However, these efforts have not yet led to improvements in risk assessment compared with existing risk scores. Noncoding RNAs (ncRNAs) have been assessed as biomarkers, and miRNAs have attracted most attention. More recently, other ncRNA classes have been identified, including long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs). Here, we compare emerging ncRNA biomarkers in the cardiovascular field with protein biomarkers for their potential in clinical application, focusing on myocardial injury.

Remodeling of repolarization and arrhythmia susceptibility in a myosin-binding protein C knockout mouse model

Hypertrophic cardiomyopathy (HCM) is one of the most common genetic cardiac diseases and among the leading causes of sudden cardiac death (SCD) in the young. The cellular mechanisms leading to SCD in HCM are not well known. Prolongation of the action potential (AP) duration (APD) is a common feature predisposing hypertrophied hearts to SCD. Previous studies have explored the roles of inward Na⁺ and Ca²⁺ in the development of HCM, but the role of repolarizing K⁺ currents has not been defined. The objective of this study was to characterize the arrhythmogenic phenotype and cellular electrophysiological properties of mice with HCM, induced by myosin-binding protein C (MyBPC) knockout (KO), and to test the hypothesis that remodeling of repolarizing K⁺ currents causes APD prolongation in MyBPC KO myocytes. We demonstrated that MyBPC KO mice developed severe hypertrophy and cardiac dysfunction compared with wild-type (WT) control mice. Telemetric electrocardiographic recordings of awake mice revealed prolongation of the corrected QT interval in the KO compared with WT control mice, with overt ventricular arrhythmias. Whole cell current- and voltage-clamp experiments comparing KO with WT mice demonstrated ventricular myocyte hypertrophy, AP prolongation, and decreased repolarizing K⁺ currents. Quantitative RT-PCR analysis revealed decreased mRNA levels of several key K⁺ channel subunits. In conclusion, decrease in repolarizing K⁺ currents in MyBPC KO ventricular myocytes contributes to AP and corrected QT interval prolongation and could account for the arrhythmia susceptibility.NEW & NOTEWORTHY Ventricular myocytes isolated from the myosin-binding protein C knockout hypertrophic cardiomyopathy mouse model demonstrate decreased repolarizing K⁺ currents and action potential and QT interval prolongation, linking cellular repolarization abnormalities with arrhythmia susceptibility and the risk for sudden cardiac death in hypertrophic cardiomyopathy.

Lack of Phenotypic Differences by Cardiovascular Magnetic Resonance Imaging in MYH7 (β-Myosin Heavy Chain)- Versus MYBPC3 (Myosin-Binding Protein C)-Related Hypertrophic Cardiomyopathy

BACKGROUND

The 2 most commonly affected genes in hypertrophic cardiomyopathy (HCM) are MYH7 (β-myosin heavy chain) and MYBPC3 (β-myosin-binding protein C). Phenotypic differences between patients with mutations in these 2 genes have been inconsistent. Scarce data exist on the genotype-phenotype association as assessed by tomographic imaging using cardiac magnetic resonance imaging.

METHODS AND RESULTS

Cardiac magnetic resonance imaging was performed on 358 consecutive genotyped hypertrophic cardiomyopathy probands at 5 tertiary hypertrophic cardiomyopathy centers. Genetic testing revealed a pathogenic mutation in 159 patients (44.4%). The most common genes identified were MYH7 (n=53) and MYBPC3 (n=75); 33.1% and 47% of genopositive patients, respectively. Phenotypic characteristics by cardiac magnetic resonance imaging of these 2 groups were similar, including left ventricular volumes, mass, maximal wall thickness, morphology, left atrial volume, and mitral valve leaflet lengths (all P=non-significant). The presence of late gadolinium enhancement (65% versus 64%; P=0.99) and the proportion of total left ventricular mass (%late gadolinium enhancement; 10.4±13.2% versus 8.5±8.5%; P=0.44) were also similar.

CONCLUSIONS

This multicenter multinational study shows lack of phenotypic differences between MYH7- and MYBPC3-associated hypertrophic cardiomyopathy when assessed by cardiac magnetic resonance imaging. Postmutational mechanisms appear more relevant to thick-filament disease expression and outcome than the disease-causing variant per se.

Sarcomeric protein modification during adrenergic stress enhances cross-bridge kinetics and cardiac output

Molecular adaptations to chronic neurohormonal stress, including sarcomeric protein cleavage and phosphorylation, provide a mechanism to increase ventricular contractility and enhance cardiac output, yet the link between sarcomeric protein modifications and changes in myocardial function remains unclear. To examine the effects of neurohormonal stress on posttranslational modifications of sarcomeric proteins, mice were administered combined α- and β-adrenergic receptor agonists (isoproterenol and phenylephrine, IPE) for 14 days using implantable osmotic pumps. In addition to significant cardiac hypertrophy and increased maximal ventricular pressure, IPE treatment accelerated pressure development and relaxation (74% increase in dP/dt_max and 14% decrease in τ), resulting in a 52% increase in cardiac output compared with saline (SAL)-treated mice. Accelerated pressure development was maintained when accounting for changes in heart rate and preload, suggesting that myocardial adaptations contribute to enhanced ventricular contractility. Ventricular myocardium isolated from IPE-treated mice displayed a significant reduction in troponin I (TnI) and myosin-binding protein C (MyBP-C) expression and a concomitant increase in the phosphorylation levels of the remaining TnI and MyBP-C protein compared with myocardium isolated from saline-treated control mice. Skinned myocardium isolated from IPE-treated mice displayed a significant acceleration in the rate of cross-bridge (XB) detachment (46% increase) and an enhanced magnitude of XB recruitment (43% increase) at submaximal Ca²⁺ activation compared with SAL-treated mice but unaltered myofilament Ca²⁺ sensitivity of force generation. These findings demonstrate that sarcomeric protein modifications during neurohormonal stress are molecular adaptations that enhance in vivo ventricular contractility through accelerated XB kinetics to increase cardiac output.NEW & NOTEWORTHY Posttranslational modifications to sarcomeric regulatory proteins provide a mechanism to modulate cardiac function in response to stress. In this study, we demonstrate that neurohormonal stress produces modifications to myosin-binding protein C and troponin I, including a reduction in protein expression within the sarcomere and increased phosphorylation of the remaining protein, which serve to enhance cross-bridge kinetics and increase cardiac output. These findings highlight the importance of sarcomeric regulatory protein modifications in modulating ventricular function during cardiac stress.

Molecular Screen Identifies Cardiac Myosin-Binding Protein-C as a Protein Kinase G-Iα Substrate

BACKGROUND

Pharmacological activation of cGMP-dependent protein kinase G I (PKGI) has emerged as a therapeutic strategy for humans with heart failure. However, PKG-activating drugs have been limited by hypotension arising from PKG-induced vasodilation. PKGIα antiremodeling substrates specific to the myocardium might provide targets to circumvent this limitation, but currently remain poorly understood.

METHODS AND RESULTS

We performed a screen for myocardial proteins interacting with the PKGIα leucine zipper (LZ)-binding domain to identify myocardial-specific PKGI antiremodeling substrates. Our screen identified cardiac myosin-binding protein-C (cMyBP-C), a cardiac myocyte-specific protein, which has been demonstrated to inhibit cardiac remodeling in the phosphorylated state, and when mutated leads to hypertrophic cardiomyopathy in humans. GST pulldowns and precipitations with cGMP-conjugated beads confirmed the PKGIα-cMyBP-C interaction in myocardial lysates. In vitro studies demonstrated that purified PKGIα phosphorylates the cMyBP-C M-domain at Ser-273, Ser-282, and Ser-302. cGMP induced cMyBP-C phosphorylation at these residues in COS cells transfected with PKGIα, but not in cells transfected with LZ mutant PKGIα, containing mutations to disrupt LZ substrate binding. In mice subjected to left ventricular pressure overload, PKGI activation with sildenafil increased cMyBP-C phosphorylation at Ser-273 compared with untreated mice. cGMP also induced cMyBP-C phosphorylation in isolated cardiac myocytes.

CONCLUSIONS

Taken together, these data support that PKGIα and cMyBP-C interact in the heart and that cMyBP-C is an anti remodeling PKGIα kinase substrate. This study provides the first identification of a myocardial-specific PKGIα LZ-dependent antiremodeling substrate and supports further exploration of PKGIα myocardial LZ substrates as potential therapeutic targets for heart failure.

Phosphorylation of cardiac Myosin-binding protein-C is a critical mediator of diastolic function

BACKGROUND

Heart failure (HF) with preserved ejection fraction (HFpEF) accounts for ≈50% of all cases of HF and currently has no effective treatment. Diastolic dysfunction underlies HFpEF; therefore, elucidation of the mechanisms that mediate relaxation can provide new potential targets for treatment. Cardiac myosin-binding protein-C (cMyBP-C) is a thick filament protein that modulates cross-bridge cycling rates via alterations in its phosphorylation status. Thus, we hypothesize that phosphorylated cMyBP-C accelerates the rate of cross-bridge detachment, thereby enhancing relaxation to mediate diastolic function.

METHODS AND RESULTS

We compared mouse models expressing phosphorylation-deficient cMyBP-C(S273A/S282A/S302A)-cMyBP-C(t3SA), phosphomimetic cMyBP-C(S273D/S282D/S302D)-cMyBP-C(t3SD), and wild-type-control cMyBP-C(tWT) to elucidate the functional effects of cMyBP-C phosphorylation. Decreased voluntary running distances, increased lung/body weight ratios, and increased brain natriuretic peptide levels in cMyBP-C(t3SA) mice demonstrate that phosphorylation deficiency is associated with signs of HF. Echocardiography (ejection fraction and myocardial relaxation velocity) and pressure/volume measurements (-dP/dtmin, pressure decay time constant τ-Glantz, and passive filling stiffness) show that cMyBP-C phosphorylation enhances myocardial relaxation in cMyBP-C(t3SD) mice, whereas deficient cMyBP-C phosphorylation causes diastolic dysfunction with HFpEF in cMyBP-C(t3SA) mice. Simultaneous force and [Ca(2+)]i measurements on intact papillary muscles show that enhancement of relaxation in cMyBP-C(t3SD) mice and impairment of relaxation in cMyBP-C(t3SA) mice are not because of altered [Ca(2+)]i handling, implicating that altered cross-bridge detachment rates mediate these changes in relaxation rates.

CONCLUSIONS

cMyBP-C phosphorylation enhances relaxation, whereas deficient phosphorylation causes diastolic dysfunction and phenotypes resembling HFpEF. Thus, cMyBP-C is a potential target for treatment of HFpEF.

Compound heterozygosity deteriorates phenotypes of hypertrophic cardiomyopathy with founder MYBPC3 mutation: evidence from patients and zebrafish models

Although most founder mutation carriers of hypertrophic cardiomyopathy (HCM), such as the cardiac myosin-binding protein C gene (MYBPC3), arose from a common ancestor exhibit favorable clinical phenotypes, there still remain small fractions of these carriers associated with increased cardiovascular events. However, few data exist regarding the defining factors that modify phenotypes of these patients, particularly in terms of multiple gene mutations. Therefore, we assessed genotype-phenotype correlations and investigated factors that contribute to phenotypic diversities of mutation carriers from 488 unrelated HCM probands. A prevalent founder mutation (Val762Asp) in MYBPC3 was identified in 33 subjects from 19 families. Among them, 28 carriers harbored an isolated Val762Asp mutation and exhibited a late onset of overt HCM compared with other MYBPC3 mutation carriers (62.8 ± 3.0 vs 50.1 ± 2.6 yr, P < 0.05). In contrast, the remaining five carriers had additional sarcomere gene mutations (3 carriers in MYBPC3 and 2 carriers in the cardiac troponin T gene). Of these five carriers, two carriers showed early disease onset and one carrier exhibited end-stage HCM. These phenotypes were recapitulated in zebrafish models; injection of MYBPC3 Val762Asp alone did not alter ventricular size or function, but ventricular dimension was significantly increased when MYBPC3 Val762Asp mRNA was coinjected with MYBPC3 Arg820Gln mRNA. These results demonstrate that MYBPC3 Val762Asp may be associated with unfavorable HCM phenotypes in some cases when combined with another MYBPC3 mutation.

AT1 blockade abolishes left ventricular hypertrophy in heterozygous cMyBP-C null mice: role of FHL1

This research investigated the impact of angiotensin AT1 receptor (Agtr1) blockade on left ventricular (LV) hypertrophy in a mouse model of human hypertrophic cardiomyopathy (HCM), which carries one functional allele of Mybpc3 gene coding cardiac myosin-binding protein C (cMyBP-C). Five-month-old heterozygous cMyBP-C knockout (Het-KO) and wild-type mice were treated with irbesartan (50 mg/kg/day) or vehicle for 8 weeks. Arterial blood pressure was measured by tail cuff plethysmography. LV dimension and function were accessed by echocardiography. Myocardial gene expression was evaluated using RT-qPCR. Compared with wild-type littermates, Het-KO mice had greater LV/body weight ratio (4.0 ± 0.1 vs. 3.3 ± 0.1 mg/g, P < 0.001), thicker interventricular septal wall (0.70 ± 0.02 vs. 0.65 ± 0.01 mm, P < 0.02), lower Mybpc3 mRNA level (-43%, P < 0.02), higher four-and-a-half LIM domains 1 (Fhl1, +110%, P < 0.01), and angiotensin-converting enzyme 1 (Ace1, +67%, P < 0.05), but unchanged Agtr1 mRNA levels in the septum. Treatment with irbesartan had no effect in wild-type mice but abolished septum-predominant LV hypertrophy and Fhl1 upregulation without changes in Ace1 but with an increased Agtr1 (+42%) in Het-KO mice. Thus, septum-predominant LV hypertrophy in Het-KO mice is combined with higher Fhl1 expression, which can be abolished by AT1 receptor blockade, indicating a role of the renin-angiotensin system and Fhl1 in cMyBP-C-related HCM.