Natural variant frequencies across domains from different sarcomere proteins cross-correlate to identify inter-protein contacts associated with cardiac muscle function and disease

Coordinated sarcomere proteins produce contraction force for muscle shortening. In human ventriculum they include the cardiac myosin motor (βmys), repetitively converting ATP free energy into work, and myosin binding protein C (MYBPC3) that in complex with βmys is regulatory. Single nucleotide variants (SNVs) causing hereditary heart diseases frequently target this protein pair. The βmys/MYBPC3 complex models a regulated motor and is used here to study how the proteins couple. SNVs in βmys or MYBPC3 survey human populations worldwide. Their protein expression modifies domain structure affecting phenotype and pathogenicity outcomes. When the SNV modified domain locates to inter-protein contacts it could affect complex coordination. Domains involved, one in βmys the other in MYBPC3, form coordinated domains (co-domains). Co-domain bilateral structure implies the possibility for a shared impact from SNV modification in either domain suggesting a correlated response to a common perturbation could identify their location. Genetic divergence over human populations is proposed to perturb SNV probability coupling that is detected by cross-correlation in 2D correlation genetics (2D-CG). SNV probability data and 2D-CG identify three critical sites, two in MYBPC3 with links to several domains across the βmys motor, and, one in βmys with links to the MYBPC3 regulatory domain. MYBPC3 sites are hinges sterically enabling regulatory interactions with βmys. The βmys site is the actin binding C-loop (residues 359-377). The C-loop is a trigger for actin-activated myosin ATPase and a contraction velocity modulator. Co-domain identification implies their spatial proximity suggesting a novel approach for in vivo protein complex structure determination.

Reconnoitering the Role of Long-Noncoding RNAs in Hypertrophic Cardiomyopathy: A Descriptive Review

Hypertrophic cardiomyopathy (HCM) is the most common form of hereditary cardiomyopathy. It is characterized by an unexplained non-dilated hypertrophy of the left ventricle with a conserved or elevated ejection fraction. It is a genetically heterogeneous disease largely caused by variants of genes encoding for cardiac sarcomere proteins, including MYH7, MYBPC3, ACTC1, TPM1, MYL2, MYL3, TNNI3, and TNNT23. Preclinical evidence indicates that the enhanced calcium sensitivity of the myofilaments plays a key role in the pathophysiology of HCM. Notably, this is not always a direct consequence of sarcomeric variations but may also result from secondary mutation-driven alterations. Long non-coding RNAs (lncRNAs) are a large class of transcripts ≥200 nucleotides in length that do not encode proteins. Compared to coding mRNAs, most lncRNAs are not as well-annotated and their functions are greatly unexplored. Nevertheless, increasing evidence shows that lncRNAs are involved in a variety of biological processes and diseases including HCM. Accumulating evidence has indicated that lncRNAs are dysregulated in HCM, and closely related to sarcomere construction, calcium channeling and homeostasis of mitochondria. In this review, we have summarized the known regulatory and functional roles of lncRNAs in HCM.

Investigating the correlation of muscle function tests and sarcomere organization in C. elegans

Background Caenorhabditis elegans has been widely used as a model to study muscle structure and function. Its body wall muscle is functionally and structurally similar to vertebrate skeletal muscle with conserved molecular pathways contributing to sarcomere structure, and muscle function. However, a systematic investigation of the relationship between muscle force and sarcomere organization is lacking. Here, we investigate the contribution of various sarcomere proteins and membrane attachment components to muscle structure and function to introduce C. elegans as a model organism to study the genetic basis of muscle strength. Methods We employ two recently developed assays that involve exertion of muscle forces to investigate the correlation of muscle function to sarcomere organization. We utilized a microfluidic pillar-based platform called NemaFlex that quantifies the maximum exertable force and a burrowing assay that challenges the animals to move in three dimensions under a chemical stimulus. We selected 20 mutants with known defects in various substructures of sarcomeres and compared the physiological function of muscle proteins required for force generation and transmission. We also characterized the degree of sarcomere disorganization using immunostaining approaches. Results We find that mutants with genetic defects in thin filaments, thick filaments, and M-lines are generally weaker, and our assays are successful in detecting the functional changes in response to each sarcomere location tested. We find that the NemaFlex and burrowing assays are functionally distinct informing on different aspects of muscle physiology. Specifically, the burrowing assay has a larger bandwidth in phenotyping muscle mutants, because it could pick ten additional mutants impaired while exerting normal muscle force in NemaFlex. This enabled us to combine their readouts to develop an integrated muscle function score that was found to correlate with the score for muscle structure disorganization. Conclusions Our results highlight the suitability of NemaFlex and burrowing assays for evaluating muscle physiology of C. elegans. Using these approaches, we discuss the importance of the studied sarcomere proteins for muscle function and structure. The scoring methodology we have developed enhances the utility of C. elegans as a genetic model to study muscle function.

Investigating the correlation of muscle function tests and sarcomere organization in C. elegans

BackgroundCaenorhabditis elegans has been widely used as a model to study muscle structure and function due to many genes having human homologs. Its body wall muscle is functionally and structurally similar to vertebrate skeletal muscle with conserved molecular pathways contributing to sarcomere structure, and muscle function. However, a systematic investigation of the relationship between muscle force and sarcomere organization is lacking. Here, we investigate the contribution of various sarcomere proteins and membrane attachment components to muscle structure and function to introduce C. elegans as a model organism to study the genetic basis of muscle strength.MethodsWe employ two recently developed assays that involve exertion of muscle forces to investigate the correlation of muscle function to sarcomere organization. We utilized a microfluidic pillar-based platform called NemaFlex that quantifies the maximum exertable force and a burrowing assay that challenges the animals to move in three dimensions under a chemical stimulus. We selected 20 mutants with known defects in various substructures of sarcomeres and compared the physiological function of muscle proteins required for force generation and transmission. We also characterized the degree of sarcomere disorganization using immunostaining approaches.ResultsWe find that mutants with genetic defects in thin filaments, thick filaments and M-lines are generally weaker, and our assays are successful in detecting the functional changes in response to each sarcomere location tested. We find that the NemaFlex and burrowing assays are functionally distinct informing on different aspects of muscle physiology. Specifically, the burrowing assay has a larger bandwidth in phenotyping muscle mutants, because it could pick ten additional mutants impaired while exerting normal muscle force in NemaFlex. This enabled us to combine their readouts to develop an integrated muscle function score that was found to correlate with the score for muscle structure disorganization.ConclusionsOur results highlight the suitability of NemaFlex and burrowing assays for evaluating muscle physiology of C. elegans. Using these approaches, we discuss the importance of the studied sarcomere proteins for muscle function and structure. The scoring methodology we have developed lays the foundation for investigating the contribution of conserved sarcomere proteins and membrane attachment components to human muscle function and strength.

Multiscale modeling of twitch contractions in cardiac trabeculae

Understanding the dynamics of a cardiac muscle twitch contraction is complex because it requires a detailed understanding of the kinetic processes of the Ca2+ transient, thin-filament activation, and the myosin–actin cross-bridge chemomechanical cycle. Each of these steps has been well defined individually, but understanding how all three of the processes operate in combination is a far more complex problem. Computational modeling has the potential to provide detailed insight into each of these processes, how the dynamics of each process affect the complexity of contractile behavior, and how perturbations such as mutations in sarcomere proteins affect the complex interactions of all of these processes. The mechanisms involved in relaxation of tension during a cardiac twitch have been particularly difficult to discern due to nonhomogeneous sarcomere lengthening during relaxation. Here we use the multiscale MUSICO platform to model rat trabecular twitches. Validation of computational models is dependent on being able to simulate different experimental datasets, but there has been a paucity of data that can provide all of the required parameters in a single experiment, such as simultaneous measurements of force, intracellular Ca2+ transients, and sarcomere length dynamics. In this study, we used data from different studies collected under similar experimental conditions to provide information for all the required parameters. Our simulations established that twitches either in an isometric sarcomere or in fixed-length, multiple-sarcomere trabeculae replicate the experimental observations if models incorporate a length–tension relationship for the nonlinear series elasticity of muscle preparations and a scheme for thick-filament regulation. The thick-filament regulation assumes an off state in which myosin heads are parked onto the thick-filament backbone and are unable to interact with actin, a state analogous to the super-relaxed state. Including these two mechanisms provided simulations that accurately predict twitch contractions over a range of different conditions.

Unraveling the Genotype‐Phenotype Relationship in Hypertrophic Cardiomyopathy: Obesity‐Related Cardiac Defects as a Major Disease Modifier

Hypertrophic cardiomyopathy (HCM) is the most common inherited cardiomyopathy and is characterized by asymmetric septal thickening and diastolic dysfunction. More than 1500 mutations in genes encoding sarcomere proteins are associated with HCM. However, the genotype‐phenotype relationship in HCM is incompletely understood and involves modification by additional disease hits. Recent cohort studies identify obesity as a major adverse modifier of disease penetrance, severity, and clinical course. In this review, we provide an overview of these clinical findings. Moreover, we explore putative mechanisms underlying obesity‐induced sensitization and aggravation of the HCM phenotype. We hypothesize obesity‐related stressors to impact on cardiomyocyte structure, metabolism, and homeostasis. These may impair cardiac function by directly acting on the primary mutation‐induced myofilament defects and by independently adding to the total cardiac disease burden. Last, we address important clinical and pharmacological implications of the involvement of obesity in HCM disease modification.

Cardiomyocyte Contractile Impairment in Heart Failure Results from Reduced BAG3-mediated Sarcomeric Protein Turnover

ABSTRACTThe association between reduced myofilament force-generating capacity (Fmax) and heart failure (HF) is clear, however the underlying molecular mechanisms are poorly understood. Here, we show the Fmax decrease arises from impaired BAG3-mediated sarcomere turnover. Myofilament BAG3 decreased in human HF and predicted the extent of Fmax decrease. This relationship was confirmed using BAG3+/- mice, which had reduced Fmax and increased myofilament ubiquitination, suggesting impaired protein turnover. We show cardiac BAG3 operates via the chaperone-assisted selective autophagy complex (CASA), conserved from skeletal muscle, and confirm sarcomeric CASA localization is BAG3/proteotoxic stress-dependent. To determine if increasing BAG3 expression in HF would restore sarcomere proteostasis/Fmax, HF mice were treated with AAV9/BAG3. Gene therapy fully restored Fmax after four weeks and decreased ubiquitination. Using mass spectrometry, we identified several sarcomere proteins with increased ubiquitination in HF and four that decreased with AAV9/BAG3. Our findings indicate BAG3-mediated sarcomere turnover is required for myofilament functional maintenance.

Protein Quality Control Activation and Microtubule Remodeling in Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is the most common inherited cardiac disorder. It is mainly caused by mutations in genes encoding sarcomere proteins. Mutant forms of these highly abundant proteins likely stress the protein quality control (PQC) system of cardiomyocytes. The PQC system, together with a functional microtubule network, maintains proteostasis. We compared left ventricular (LV) tissue of nine donors (controls) with 38 sarcomere mutation-positive (HCMSMP) and 14 sarcomere mutation-negative (HCMSMN) patients to define HCM and mutation-specific changes in PQC. Mutations in HCMSMP result in poison polypeptides or reduced protein levels (haploinsufficiency, HI). The main findings were 1) several key PQC players were more abundant in HCM compared to controls, 2) after correction for sex and age, stabilizing heat shock protein (HSP)B1, and refolding, HSPD1 and HSPA2 were increased in HCMSMP compared to controls, 3) α-tubulin and acetylated α-tubulin levels were higher in HCM compared to controls, especially in HCMHI, 4) myosin-binding protein-C (cMyBP-C) levels were inversely correlated with α-tubulin, and 5) α-tubulin levels correlated with acetylated α-tubulin and HSPs. Overall, carrying a mutation affects PQC and α-tubulin acetylation. The haploinsufficiency of cMyBP-C may trigger HSPs and α-tubulin acetylation. Our study indicates that proliferation of the microtubular network may represent a novel pathomechanism in cMyBP-C haploinsufficiency-mediated HCM.

Fibre-type-specific and Mitochondrial Biomarkers of Muscle Damage after Mountain Races

Consequences of running mountain races on muscle damage were investigated by analysing serum muscle enzymes and fibre-type-specific sarcomere proteins. We studied 10 trained amateur and 6 highly trained runners who ran a 35 km and 55 km mountain trail race (MTR), respectively. Levels of creatine kinase (CK), CK-MB isoform (CK-MB), sarcomeric mitochondrial CK (sMtCK), transaminases (AST and ALT), cardiac troponin I (cTnI) and fast (FM) and slow myosin (SM) isoforms, were assessed before, 1 h, 24 h and 48 h after the beginning of MTR. Significant SM increases were found at 24 h in the 55 km group. Levels of CK, CK-MB, AST and cTnI were significantly elevated in both groups following MTR, but in the 55 km group they tended to stabilize in at 48 h. Using pooled data, time-independent serum peaks of SM and CK-MB were significantly correlated. Moreover, concentration of sMtCK was significantly elevated at 1 and 24 h after the race in the 35 km group. Although training volume could confer protection on the mitochondria, the increase in serum CK-MB and SM in the 55 km group might be related to damage to the contractile apparatus type I fibres. Competing in long-distance MTRs might be related to deeper type I muscle fibre damage, even in highly trained individuals