Functional divergence of the sarcomeric myosin, MYH7b, supports species-specific biological roles

Postprandial cardiac hypertrophy is sustained by mechanics, epigenetic, and metabolic reprogramming in pythons

Constricting pythons, known for their ability to consume infrequent, massive meals, exhibit rapid and reversible cardiac hypertrophy following feeding. Our primary goal was to investigate how python hearts achieve this adaptive response after feeding. Isolated myofibrils increased force after feeding without changes in sarcomere ultrastructure and without increasing energy cost. Ca2+ transients were prolonged after feeding with no changes in myofibril Ca2+ sensitivity. Feeding reduced titin-based tension, resulting in decreased cardiac tissue stiffness. Feeding also reduced the activity of sirtuins, a metabolically linked class of histone deacetylases, and increased chromatin accessibility. Transcription factor enrichment analysis on transposase-accessible chromatin with sequencing revealed the prominent role of transcription factors Yin Yang1 and NRF1 in postfeeding cardiac adaptation. Gene expression also changed with the enrichment of translation and metabolism. Finally, metabolomics analysis and adenosine triphosphate production demonstrated that cardiac adaptation after feeding not only increased energy demand but also energy production. These findings have broad implications for our understanding of cardiac adaptation across species and hold promise for the development of innovative approaches to address cardiovascular diseases.

The influence of biological sex in human skeletal muscle transcriptome during ageing

Sex is a crucial biological variable, and influence of biological sex on the change of gene expression in ageing skeletal muscle has not yet been fully revealed. In this study, the mRNA expression profiles were obtained from the Gene Expression Omnibus database. Key genes were identified by differential expression analysis and weighted gene co-expression network analysis. The gene set enrichment analysis software and Molecular Signatures Database were used for functional and enrichment analysis. A protein-protein interaction network was constructed using STRING and visualized in Cytoscape. The results were compared between female and male subgroups. Differentially expressed genes and enriched pathways in different sex subgroups shared only limited similarities. The pathways enriched in the female subgroup were more similar to the pathways enriched in the older groups without taking sex difference into consideration. The pathways enriched in the female subgroup were more similar to the pathways enriched in the older groups without taking sex difference into consideration. The muscle myosin filament pathways were downregulated in the both aged female and male samples whereas transforming growth factor beta pathway and extracellular matrix-related pathways were upregulated. With muscle ageing, the metabolism-related pathways, protein synthesis and degradation pathways, results of predicted immune cell infiltration, and gene cluster associated with slow-type myofibers drastically different between the female and male subgroups. This finding may indicate that changes in muscle type with ageing may differ between the sexes in vastus lateralis muscle.

Developmental, Physiological and Phylogenetic Perspectives on the Expression and Regulation of Myosin Heavy Chains in Craniofacial Muscles

This review deals with the developmental origins of extraocular, jaw and laryngeal muscles, the expression, regulation and functional significance of sarcomeric myosin heavy chains (MyHCs) that they express and changes in MyHC expression during phylogeny. Myogenic progenitors from the mesoderm in the prechordal plate and branchial arches specify craniofacial muscle allotypes with different repertoires for MyHC expression. To cope with very complex eye movements, extraocular muscles (EOMs) express 11 MyHCs, ranging from the superfast extraocular MyHC to the slowest, non-muscle MyHC IIB (nmMyH IIB). They have distinct global and orbital layers, singly- and multiply-innervated fibres, longitudinal MyHC variations, and palisade endings that mediate axon reflexes. Jaw-closing muscles express the high-force masticatory MyHC and cardiac or limb MyHCs depending on the appropriateness for the acquisition and mastication of food. Laryngeal muscles express extraocular and limb muscle MyHCs but shift toward expressing slower MyHCs in large animals. During postnatal development, MyHC expression of craniofacial muscles is subject to neural and hormonal modulation. The primary and secondary myotubes of developing EOMs are postulated to induce, via different retrogradely transported neurotrophins, the rich diversity of neural impulse patterns that regulate the specific MyHCs that they express. Thyroid hormone shifts MyHC 2A toward 2B in jaw muscles, laryngeal muscles and possibly extraocular muscles. This review highlights the fact that the pattern of myosin expression in mammalian craniofacial muscles is principally influenced by the complex interplay of cell lineages, neural impulse patterns, thyroid and other hormones, functional demands and body mass. In these respects, craniofacial muscles are similar to limb muscles, but they differ radically in the types of cell lineage and the nature of their functional demands.

Toward physics-based precision medicine: Exploiting protein dynamics to design new therapeutics and interpret variants

The goal of precision medicine is to utilize our knowledge of the molecular causes of disease to better diagnose and treat patients. However, there is a substantial mismatch between the small number of food and drug administration (FDA)-approved drugs and annotated coding variants compared to the needs of precision medicine. This review introduces the concept of physics-based precision medicine, a scalable framework that promises to improve our understanding of sequence-function relationships and accelerate drug discovery. We show that accounting for the ensemble of structures a protein adopts in solution with computer simulations overcomes many of the limitations imposed by assuming a single protein structure. We highlight studies of protein dynamics and recent methods for the analysis of structural ensembles. These studies demonstrate that differences in conformational distributions predict functional differences within protein families and between variants. Thanks to new computational tools that are providing unprecedented access to protein structural ensembles, this insight may enable accurate predictions of variant pathogenicity for entire libraries of variants. We further show that explicitly accounting for protein ensembles, with methods like alchemical free energy calculations or docking to Markov state models, can uncover novel lead compounds. To conclude, we demonstrate that cryptic pockets, or cavities absent in experimental structures, provide an avenue to target proteins that are currently considered undruggable. Taken together, our review provides a roadmap for the field of protein science to accelerate precision medicine.

Homologous mutations in human β, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in β-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman-Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known whether their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human β, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins but minimal effects in β myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing overall enzymatic (ATPase) cycle rate. In contrast, the only measured effect of R671C in β myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not β, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are a testament to myosin's highly allosteric nature.

MYH13, a superfast myosin expressed in extraocular, laryngeal and syringeal muscles

MYH13 is a unique type of sarcomeric myosin heavy chain (MYH) first detected in mammalian extraocular (EO) muscles and later also in vocal muscles, including laryngeal muscles of some mammals and syringeal muscles of songbirds. All these muscles are specialized in generating very fast contractions while producing relatively low force, a design appropriate for muscles acting against a much lower load than most skeletal muscles inserting into the skeleton. The definition of the physiological properties of muscle fibres containing MYH13 has been complicated by the mixed fibre type composition of EO muscles and the coexistence of different MYH types within the same fibre. A major advance in this area came from studies on isolated recombinant myosin motors and the demonstration that the affinity of actin-bound human MYH13 for ADP is much weaker than those of fast-type MYH1 (type 2X) and MYH2 (type 2A). This property is consistent with a very fast detachment of myosin from actin, a major determinant of shortening velocity. The MYH13 gene arose early during vertebrate evolution but was characterized only in mammals and birds and appears to have been lost in some teleost fish. The MYH13 gene is located at the 3' end of the mammalian fast/developmental gene cluster and in a similar position to the orthologous cluster in syntenic regions of the songbird genome. MYH13 gene regulation is controlled by a super-enhancer in the mammalian locus and deletion of the neighbouring fast MYH1 and MYH4 genes leads to abnormal MYH13 expression in mouse leg muscles.

Discovery of a novel cardiac-specific myosin modulator using artificial intelligence-based virtual screening

Direct modulation of cardiac myosin function has emerged as a therapeutic target for both heart disease and heart failure. However, the development of myosin-based therapeutics has been hampered by the lack of targeted in vitro screening assays. In this study we use Artificial Intelligence-based virtual high throughput screening (vHTS) to identify novel small molecule effectors of human β-cardiac myosin. We test the top scoring compounds from vHTS in biochemical counter-screens and identify a novel chemical scaffold called 'F10' as a cardiac-specific low-micromolar myosin inhibitor. Biochemical and biophysical characterization in both isolated proteins and muscle fibers show that F10 stabilizes both the biochemical (i.e. super-relaxed state) and structural (i.e. interacting heads motif) OFF state of cardiac myosin, and reduces force and left ventricular pressure development in isolated myofilaments and Langendorff-perfused hearts, respectively. F10 is a tunable scaffold for the further development of a novel class of myosin modulators.

Expression of MyHC-15 and -2x in human muscle spindles: An immunohistochemical study

To build on the existing data on the pattern of myosin heavy chain (MyHC) isoforms expression in the human muscle spindles, we aimed to verify whether the 'novel' MyHC-15, -2x and -2b isoforms are co-expressed with the other known isoforms in the human intrafusal fibres. Using a set of antibodies, we attempted to demonstrate nine isoforms (15, slow-tonic, 1, α, 2a, 2x, 2b, embryonic, neonatal) in different regions of intrafusal fibres in the biceps brachii and flexor digitorum profundus muscles. The reactivity of some antibodies with the extrafusal fibres was also tested in the masseter and laryngeal cricothyreoid muscles. In both upper limb muscles, the expression of slow-tonic isoform was a reliable marker for differentiating positive bag fibres from negative chain fibres. Generally, bag1 and bag2 fibres were distinguished in isoform 1 expression; the latter consistently expressed this isoform over their entire length. Although isoform 15 was not abundantly expressed in intrafusal fibres, its expression was pronounced in the extracapsular region of bag fibres. Using a 2x isoform-specific antibody, this isoform was demonstrated in the intracapsular regions of some intrafusal fibres, particularly chain fibres. To the best of our knowledge, this study is the first to demonstrate 15 and 2x isoforms in human intrafusal fibres. However, whether the labelling with an antibody specific for rat 2b isoform reflects the expression of this isoform in bag fibres and some extrafusal ones in the specialised cranial muscles requires further evaluation. The revealed pattern of isoform co-expression only partially agrees with the results of previous, more extensive studies. Nevertheless, it may be inferred that MyHC isoform expression in intrafusal fibres varies along their length, across different muscle spindles and muscles. Furthermore, the estimation of expression may also depend on the antibodies utilised, which may also react differently with intrafusal and extrafusal fibres.

Divergent Molecular Phenotypes in Point Mutations at the Same Residue in Beta-Myosin Heavy Chain Lead to Distinct Cardiomyopathies

In genetic cardiomyopathies, a frequently described phenomenon is how similar mutations in one protein can lead to discrete clinical phenotypes. One example is illustrated by two mutations in beta myosin heavy chain (β-MHC) that are linked to hypertrophic cardiomyopathy (HCM) (Ile467Val, I467V) and left ventricular non-compaction (LVNC) (Ile467Thr, I467T). To investigate how these missense mutations lead to independent diseases, we studied the molecular effects of each mutation using recombinant human β-MHC Subfragment 1 (S1) in in vitro assays. Both HCM-I467V and LVNC-I467T S1 mutations exhibited similar mechanochemical function, including unchanged ATPase and enhanced actin velocity but had opposing effects on the super-relaxed (SRX) state of myosin. HCM-I467V S1 showed a small reduction in the SRX state, shifting myosin to a more actin-available state that may lead to the "gain-of-function" phenotype commonly described in HCM. In contrast, LVNC-I467T significantly increased the population of myosin in the ultra-slow SRX state. Interestingly, molecular dynamics simulations reveal that I467T allosterically disrupts interactions between ADP and the nucleotide-binding pocket, which may result in an increased ADP release rate. This predicted change in ADP release rate may define the enhanced actin velocity measured in LVNC-I467T, but also describe the uncoupled mechanochemical function for this mutation where the enhanced ADP release rate may be sufficient to offset the increased SRX population of myosin. These contrasting molecular effects may lead to contractile dysregulation that initiates LVNC-associated signaling pathways that progress the phenotype. Together, analysis of these mutations provides evidence that phenotypic complexity originates at the molecular level and is critical to understanding disease progression and developing therapies.

Homologous mutations in β, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in β -cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known if their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human β , embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins, with the most dramatic in perinatal, but minimal effects in β myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing ATPase cycle rate. In contrast, the only measured effect of R671C in β myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not β , myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents the first direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are yet another testament to myosin's highly allosteric nature.

Right Ventricular Sarcomere Contractile Depression and the Role of Thick Filament Activation in Human Heart Failure With Pulmonary Hypertension

BACKGROUND

Right ventricular (RV) contractile dysfunction commonly occurs and worsens outcomes in patients with heart failure with reduced ejection fraction and pulmonary hypertension (HFrEF-PH). However, such dysfunction often goes undetected by standard clinical RV indices, raising concerns that they may not reflect aspects of underlying myocyte dysfunction. We thus sought to characterize RV myocyte contractile depression in HFrEF-PH, identify those components reflected by clinical RV indices, and uncover underlying biophysical mechanisms.

METHODS

Resting, calcium-, and load-dependent mechanics were prospectively studied in permeabilized RV cardiomyocytes isolated from explanted hearts from 23 patients with HFrEF-PH undergoing cardiac transplantation and 9 organ donor controls.

RESULTS

Unsupervised machine learning using myocyte mechanical data with the highest variance yielded 2 HFrEF-PH subgroups that in turn mapped to patients with decompensated or compensated clinical RV function. This correspondence was driven by reduced calcium-activated isometric tension in decompensated clinical RV function, whereas surprisingly, many other major myocyte contractile measures including peak power and myocyte active stiffness were similarly depressed in both groups. Similar results were obtained when subgroups were first defined by clinical indices, and then myocyte mechanical properties in each group compared. To test the role of thick filament defects, myofibrillar structure was assessed by x-ray diffraction of muscle fibers. This revealed more myosin heads associated with the thick filament backbone in decompensated clinical RV function, but not compensated clinical RV function, as compared with controls. This corresponded to reduced myosin ATP turnover in decompensated clinical RV function myocytes, indicating less myosin in a crossbridge-ready disordered-relaxed (DRX) state. Altering DRX proportion (%DRX) affected peak calcium-activated tension in the patient groups differently, depending on their basal %DRX, highlighting potential roles for precision-guided therapeutics. Last, increasing myocyte preload (sarcomere length) increased %DRX 1.5-fold in controls but only 1.2-fold in both HFrEF-PH groups, revealing a novel mechanism for reduced myocyte active stiffness and by extension Frank-Starling reserve in human heart failure.

CONCLUSIONS

Although there are many RV myocyte contractile deficits in HFrEF-PH, commonly used clinical indices only detect reduced isometric calcium-stimulated force, which is related to deficits in basal and recruitable %DRX myosin. Our results support use of therapies to increase %DRX and enhance length-dependent recruitment of DRX myosin heads in such patients.

Distinct effects of two hearing loss-associated mutations in the sarcomeric myosin MYH7b

For decades, sarcomeric myosin heavy chain proteins were assumed to be restricted to striated muscle where they function as molecular motors that contract muscle. However, MYH7b, an evolutionarily ancient member of this myosin family, has been detected in mammalian nonmuscle tissues, and mutations in MYH7b are linked to hereditary hearing loss in compound heterozygous patients. These mutations are the first associated with hearing loss rather than a muscle pathology, and because there are no homologous mutations in other myosin isoforms, their functional effects were unknown. We generated recombinant human MYH7b harboring the D515N or R1651Q hearing loss-associated mutation and studied their effects on motor activity and structural and assembly properties, respectively. The D515N mutation had no effect on steady-state actin-activated ATPase rate or load-dependent detachment kinetics but increased actin sliding velocity because of an increased displacement during the myosin working stroke. Furthermore, we found that the D515N mutation caused an increase in the proportion of myosin heads that occupy the disordered-relaxed state, meaning more myosin heads are available to interact with actin. Although we found no impact of the R1651Q mutation on myosin rod secondary structure or solubility, we observed a striking aggregation phenotype when this mutation was introduced into nonmuscle cells. Our results suggest that each mutation independently affects MYH7b function and structure. Together, these results provide the foundation for further study of a role for MYH7b outside the sarcomere.

Right Ventricular Sarcomere Contractile Depression and the Role of Thick Filament Activation in Human Heart Failure with Pulmonary Hypertension

Rationale

Right ventricular (RV) contractile dysfunction commonly occurs and worsens outcomes in heart failure patients with reduced ejection fraction and pulmonary hypertension (HFrEF-PH). However, such dysfunction often goes undetected by standard clinical RV indices, raising concerns that they may not reflect aspects of underlying myocyte dysfunction.

Objective

To determine components of myocyte contractile depression in HFrEF-PH, identify those reflected by clinical RV indices, and elucidate their underlying biophysical mechanisms.

Methods and Results

Resting, calcium- and load-dependent mechanics were measured in permeabilized RV cardiomyocytes isolated from explanted hearts from 23 HFrEF-PH patients undergoing cardiac transplantation and 9 organ-donor controls. Unsupervised machine learning using myocyte mechanical data with the highest variance yielded two HFrEF-PH subgroups that in turn mapped to patients with depressed (RVd) or compensated (RVc) clinical RV function. This correspondence was driven by reduced calcium-activated isometric tension in RVd, while surprisingly, many other major myocyte contractile measures including peak power, maximum unloaded shortening velocity, and myocyte active stiffness were similarly depressed in both groups. Similar results were obtained when subgroups were first defined by clinical indices, and then myocyte mechanical properties in each group compared. To test the role of thick-filament defects, myofibrillar structure was assessed by X-ray diffraction of muscle fibers. This revealed more myosin heads associated with the thick filament backbone in RVd but not RVc, as compared to controls. This corresponded to reduced myosin ATP turnover in RVd myocytes, indicating less myosin in a cross-bridge ready disordered-relaxed (DRX) state. Altering DRX proportion (%DRX) affected peak calcium-activated tension in the patient groups differently, depending on their basal %DRX, highlighting potential roles for precision-guided therapeutics. Lastly, increasing myocyte preload (sarcomere length) increased %DRX 1.5-fold in controls but only 1.2-fold in both HFrEF-PH groups, revealing a novel mechanism for reduced myocyte active stiffness and by extension Frank-Starling reserve in human HF.

Conclusions

While there are multiple RV myocyte contractile deficits In HFrEF-PH, clinical indices primarily detect reduced isometric calcium-stimulated force related to deficits in basal and recruitable %DRX myosin. Our results support use of therapies to increase %DRX and enhance length-dependent recruitment of DRX myosin heads in such patients.

Drug specificity and affinity are encoded in the probability of cryptic pocket opening in myosin motor domains

The design of compounds that can discriminate between closely related target proteins remains a central challenge in drug discovery. Specific therapeutics targeting the highly conserved myosin motor family are urgently needed as mutations in at least six of its members cause numerous diseases. Allosteric modulators, like the myosin-II inhibitor blebbistatin, are a promising means to achieve specificity. However, it remains unclear why blebbistatin inhibits myosin-II motors with different potencies given that it binds at a highly conserved pocket that is always closed in blebbistatin-free experimental structures. We hypothesized that the probability of pocket opening is an important determinant of the potency of compounds like blebbistatin. To test this hypothesis, we used Markov state models (MSMs) built from over 2 ms of aggregate molecular dynamics simulations with explicit solvent. We find that blebbistatin's binding pocket readily opens in simulations of blebbistatin-sensitive myosin isoforms. Comparing these conformational ensembles reveals that the probability of pocket opening correctly identifies which isoforms are most sensitive to blebbistatin inhibition and that docking against MSMs quantitatively predicts blebbistatin binding affinities (R2=0.82). In a blind prediction for an isoform (Myh7b) whose blebbistatin sensitivity was unknown, we find good agreement between predicted and measured IC50s (0.67 μM vs. 0.36 μM). Therefore, we expect this framework to be useful for the development of novel specific drugs across numerous protein targets.