Hypertrophic cardiomyopathy disease results from disparate impairments of cardiac myosin function and auto-inhibition

Select autosomal dominant DFNA11 deafness mutations activate Myo7A in epithelial cells

Myosin-7A (Myo7A) is a motor protein crucial for the organization and function of stereocilia, specialized actin-rich protrusions on the surface of inner ear hair cells that mediate hearing. Mutations in Myo7A cause several forms of genetic hearing loss, including autosomal dominant DFNA11 deafness. Despite its importance, the structural elements of Myo7A that control its motor activity within cells are not well understood. In this study, we used cultured kidney epithelial cells to screen for mutations that activate the motor-dependent targeting of Myo7A to the tips of apical microvilli on these cells. Our findings reveal that Myo7A is regulated by specific IQ motifs within its lever arm, and that this regulation can function at least partially independent of its tail sequence. Importantly, we demonstrate that many of the DFNA11 deafness mutations reported in patients activate Myo7A targeting, providing a potential explanation for the autosomal dominant genetics of this form of deafness.

Dilated cardiomyopathy mutation in beta-cardiac myosin enhances actin activation of the power stroke and phosphate release

Inherited mutations in human beta-cardiac myosin (M2β) can lead to severe forms of heart failure. The E525K mutation in M2β is associated with dilated cardiomyopathy (DCM) and was found to stabilize the interacting heads motif (IHM) and autoinhibited super-relaxed (SRX) state in dimeric heavy meromyosin. However, in monomeric M2β subfragment 1 (S1) we found that E525K enhances (threefold) the maximum steady-state actin-activated ATPase activity (k cat) and decreases (eightfold) the actin concentration at which ATPase is one-half maximal (K ATPase). We also found a twofold to fourfold increase in the actin-activated power stroke and phosphate release rate constants at 30 μM actin, which overall enhanced the duty ratio threefold. Loaded motility assays revealed that the enhanced intrinsic motor activity translates to increased ensemble force in M2β S1. Glutamate 525, located near the actin binding region in the so-called activation loop, is highly conserved and predicted to form a salt bridge with another conserved residue (lysine 484) in the relay helix. Enhanced sampling molecular dynamics simulations predict that the charge reversal mutation disrupts the E525-K484 salt bridge, inducing conformations with a more flexible relay helix and a wide phosphate release tunnel. Our results highlight a highly conserved allosteric pathway associated with actin activation of the power stroke and phosphate release and suggest an important feature of the autoinhibited IHM is to prevent this region of myosin from interacting with actin. The ability of the E525K mutation to stabilize the IHM likely overrides the enhanced intrinsic motor properties, which may be key to triggering DCM pathogenesis.

An atomistic model of myosin interacting heads motif dynamics and their modification by 2'-deoxy-ADP

The contraction of striated muscle is driven by cycling myosin motor proteins embedded within the thick filaments of sarcomeres. In addition to cross-bridge cycling with actin, these myosin proteins can enter an inactive, sequestered state in which the globular S1 heads rest along the thick filament surface and are unable to perform motor activities. Structurally, this state is called the interacting heads motif (IHM) and is a critical conformational state of myosin that regulates muscle contractility and energy expenditure. Structural perturbation of the sequestered state via missense mutations can pathologically disrupt the mechanical performance of muscle tissue. Thus, the IHM state has become a target for therapeutic intervention. An ATP analogue called 2'-deoxy-ATP (dATP) is a potent myosin activator which destabilizes the IHM. Here we use molecular dynamics simulations to study the molecular mechanisms by which dATP modifies the structure and dynamics of myosin in a sequestered state. Simulations with IHM containing ADP.Pi in both nucleotide binding pockets revealed residual dynamics in an otherwise 'inactive' and 'sequestered' state of a motor protein. Replacement of ADP.Pi by dADP.Pi triggered a series of structural changes that modify the protein-protein interface that stabilizes the sequestered state, and changes to this interface were accompanied by allosteric changes in remote regions of the protein complex. A comparative analysis of these dynamics predicted new structural sites that may affect IHM stability.

Omecamtiv mecarbil and Mavacamten target the same myosin pocket despite opposite effects in heart contraction

Inherited cardiomyopathies are common cardiac diseases worldwide, leading in the late stage to heart failure and death. The most promising treatments against these diseases are small molecules directly modulating the force produced by β-cardiac myosin, the molecular motor driving heart contraction. Omecamtiv mecarbil and Mavacamten are two such molecules that completed phase 3 clinical trials, and the inhibitor Mavacamten is now approved by the FDA. In contrast to Mavacamten, Omecamtiv mecarbil acts as an activator of cardiac contractility. Here, we reveal by X-ray crystallography that both drugs target the same pocket and stabilize a pre-stroke structural state, with only few local differences. All-atom molecular dynamics simulations reveal how these molecules produce distinct effects in motor allostery thus impacting force production in opposite way. Altogether, our results provide the framework for rational drug development for the purpose of personalized medicine.

Exploring the Super-Relaxed State of Human Cardiac β-Myosin by Molecular Dynamics Simulations

Human β-cardiac myosin plays a critical role in generating the mechanical forces necessary for cardiac muscle contraction. This process relies on a delicate dynamic equilibrium between the disordered relaxed state (DRX) and the super-relaxed state (SRX) of myosin. Disruptions in this equilibrium due to mutations can lead to heart diseases. However, the structural characteristics of SRX and the molecular mechanisms underlying pathogenic mutations have remained elusive. To bridge this gap, we conducted molecular dynamics simulations and free energy calculations to explore the conformational changes in myosin. Our findings indicate that the size of the phosphate-binding pocket can serve as a valuable metric for characterizing the transition from the DRX to SRX state. Importantly, we established a global dynamic coupling network within the myosin motor head at the residue level, elucidating how the pathogenic mutation E483K impacts the equilibrium between SRX and DRX through allosteric effects. Our work illuminates molecular details of SRX and offers valuable insights into disease treatment through the regulation of SRX.

Enhanced myofilament calcium sensitivity aggravates abnormal calcium handling and diastolic dysfunction in patient-specific induced pluripotent stem cell-derived cardiomyocytes with MYH7 mutation

Hypertrophic cardiomyopathy (HCM), the most common inherited heart disease, is frequently caused by mutations in the β-cardiac myosin heavy chain gene (MYH7). Abnormal calcium handling and diastolic dysfunction are archetypical features of HCM caused by MYH7 gene mutations. However, the mechanism of how MYH7 mutations leads to these features remains unclear, which inhibits the development of effective therapies. Initially, cardiomyocytes were generated from induced pluripotent stem cells from an eight-year-old girl diagnosed with HCM carrying a MYH7(C.1063 G>A) heterozygous mutation(mutant-iPSC-CMs) and mutation-corrected isogenic iPSCs(control-iPSC-CMs) in the present study. Next, we compared phenotype of mutant-iPSC-CMs to that of control-iPSC-CMs, by assessing their morphology, hypertrophy-related genes expression, calcium handling, diastolic function and myofilament calcium sensitivity at days 15 and 40 respectively. Finally, to better understand increased myofilament Ca2+ sensitivity as a central mechanism of central pathogenicity in HCM, inhibition of calcium sensitivity with mavacamten can improveed cardiomyocyte hypertrophy. Mutant-iPSC-CMs exhibited enlarged areas, increased sarcomere disarray, enhanced expression of hypertrophy-related genes proteins, abnormal calcium handling, diastolic dysfunction and increased myofilament calcium sensitivity at day 40, but only significant increase in calcium sensitivity and mild diastolic dysfunction at day 15. Increased calcium sensitivity by levosimendan aggravates cardiomyocyte hypertrophy phenotypes such as expression of hypertrophy-related genes, abnormal calcium handling and diastolic dysfunction, while inhibition of calcium sensitivity significantly improves cardiomyocyte hypertrophy phenotypes in mutant-iPSC-CMs, suggesting increased myofilament calcium sensitivity is the primary mechanisms for MYH7 mutations pathogenesis. Our studies have uncovered a pathogenic mechanism of HCM caused by MYH7 gene mutations through which enhanced myofilament calcium sensitivity aggravates abnormal calcium handling and diastolic dysfunction. Correction of the myofilament calcium sensitivity was found to be an effective method for treating the development of HCM phenotype in vitro.

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.

Omecamtiv mecarbil and Mavacamten target the same myosin pocket despite antagonistic effects in heart contraction

Inherited cardiomyopathies are amongst the most common cardiac diseases worldwide, leading in the late-stage to heart failure and death. The most promising treatments against these diseases are small-molecules directly modulating the force produced by β-cardiac myosin, the molecular motor driving heart contraction. Two of these molecules that produce antagonistic effects on cardiac contractility have completed clinical phase 3 trials: the activator Omecamtiv mecarbil and the inhibitor Mavacamten. In this work, we reveal by X-ray crystallography that both drugs target the same pocket and stabilize a pre-stroke structural state, with only few local differences. All atoms molecular dynamics simulations reveal how these molecules can have antagonistic impact on the allostery of the motor by comparing β-cardiac myosin in the apo form or bound to Omecamtiv mecarbil or Mavacamten. Altogether, our results provide the framework for rational drug development for the purpose of personalized medicine.

Dilated cardiomyopathy mutation in beta-cardiac myosin enhances actin activation of the power stroke and phosphate release

Inherited mutations in human beta-cardiac myosin (M2β) can lead to severe forms of heart failure. The E525K mutation in M2β is associated with dilated cardiomyopathy (DCM) and was found to stabilize the interacting heads motif (IHM) and autoinhibited super-relaxed (SRX) state in dimeric heavy meromyosin. However, in monomeric M2β subfragment 1 (S1) we found that E525K enhances (3-fold) the maximum steady-state actin-activated ATPase activity (kcat) and decreases (6-fold) the actin concentration at which ATPase is one-half maximal (KATPase). We also found a 3 to 4-fold increase in the actin-activated power stroke and phosphate release rate constants at 30 μM actin, which overall enhanced the duty ratio 3-fold. Loaded motility assays revealed that the enhanced intrinsic motor activity translates to increased ensemble force in M2β S1. Glutamate 525, located near the actin binding region in the so-called activation loop, is highly conserved and predicted to form a salt-bridge with another conserved residue (lysine 484) in the relay helix. Enhanced sampling molecular dynamics simulations predict that the charge reversal mutation disrupts the E525-K484 salt-bridge, inducing conformations with a more flexible relay helix and a wide phosphate release tunnel. Our results highlight a highly conserved allosteric pathway associated with actin activation of the power stroke and phosphate release and suggest an important feature of the autoinhibited IHM is to prevent this region of myosin from interacting with actin. The ability of the E525K mutation to stabilize the IHM likely overrides the enhanced intrinsic motor properties, which may be key to triggering DCM pathogenesis.

How myosin VI traps its off-state, is activated and dimerizes

Myosin VI (Myo6) is the only minus-end directed nanomotor on actin, allowing it to uniquely contribute to numerous cellular functions. As for other nanomotors, the proper functioning of Myo6 relies on precise spatiotemporal control of motor activity via a poorly defined off-state and interactions with partners. Our structural, functional, and cellular studies reveal key features of myosin regulation and indicate that not all partners can activate Myo6. TOM1 and Dab2 cannot bind the off-state, while GIPC1 binds Myo6, releases its auto-inhibition and triggers proximal dimerization. Myo6 partners thus differentially recruit Myo6. We solved a crystal structure of the proximal dimerization domain, and show that its disruption compromises endocytosis in HeLa cells, emphasizing the importance of Myo6 dimerization. Finally, we show that the L926Q deafness mutation disrupts Myo6 auto-inhibition and indirectly impairs proximal dimerization. Our study thus demonstrates the importance of partners in the control of Myo6 auto-inhibition, localization, and activation.

Whole-Exome Sequencing Reveals Mutational Signature of Hypertrophic Cardiomyopathy

Background

Hypertrophic cardiomyopathy (HCM) is an extremely insidious and lethal disease caused by genetic variation. It has been studied for nearly 70 years since its discovery, but its cause of the disease remains a mystery. This study is aimed to explore the genetic pathogenesis of HCM in order to provide new insight for the diagnosis and treatment of HCM.

Methods

Patients with HCM at 4 hospitals from January 1, 2020, to December 31, 2021, were collected. Peripheral blood of these patients was collected for whole exome sequencing. Moreover, data on the HCM transcriptome were analyzed in the GEO database.

Results

Totally, 14 patients were enrolled, and 6 single-nucleotide variation (SNV) mutant genes represented by MUC12 were observed. Most of the gene mutations in HCM patients were synonymous and non-synonymous, and the types of base mutations were mainly C > T and G > A. Copy number variants (CNVs) predominantly occurred on chromosome 1 in HCM patients. Furthermore, we found that the only ATP2A2 gene was differentially expressed in 3 groups of transcriptome data in GEO database, and the presence of ATP2A2 mutation in 10 samples was observed in this study.

Conclusion

In summary, 7 mutated genes represented by MUC12 and ATP2A2 were found in this study, which may provide novel insights into the pathogenic mechanism of HCM.

New paradigms in actomyosin energy transduction: Critical evaluation of non-traditional models for orthophosphate release

Release of the ATP hydrolysis product ortophosphate (Pi) from the active site of myosin is central in chemo-mechanical energy transduction and closely associated with the main force-generating structural change, the power-stroke. Despite intense investigations, the relative timing between Pi-release and the power-stroke remains poorly understood. This hampers in depth understanding of force production by myosin in health and disease and our understanding of myosin-active drugs. Since the 1990s and up to today, models that incorporate the Pi-release either distinctly before or after the power-stroke, in unbranched kinetic schemes, have dominated the literature. However, in recent years, alternative models have emerged to explain apparently contradictory findings. Here, we first compare and critically analyze three influential alternative models proposed previously. These are either characterized by a branched kinetic scheme or by partial uncoupling of Pi-release and the power-stroke. Finally, we suggest critical tests of the models aiming for a unified picture.

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.

Cryo-EM structure of the folded-back state of human β-cardiac myosin

To save energy and precisely regulate cardiac contractility, cardiac muscle myosin heads are sequestered in an 'off' state that can be converted to an 'on' state when exertion is increased. The 'off' state is equated with a folded-back structure known as the interacting-heads motif (IHM), which is a regulatory feature of all class-2 muscle and non-muscle myosins. We report here the human β-cardiac myosin IHM structure determined by cryo-electron microscopy to 3.6 Å resolution, providing details of all the interfaces stabilizing the 'off' state. The structure shows that these interfaces are hot spots of hypertrophic cardiomyopathy mutations that are thought to cause hypercontractility by destabilizing the 'off' state. Importantly, the cardiac and smooth muscle myosin IHM structures dramatically differ, providing structural evidence for the divergent physiological regulation of these muscle types. The cardiac IHM structure will facilitate development of clinically useful new molecules that modulate IHM stability.

Cryo-EM structure of the folded-back state of human β-cardiac myosin

During normal levels of exertion, many cardiac muscle myosin heads are sequestered in an off-state even during systolic contraction to save energy and for precise regulation. They can be converted to an on-state when exertion is increased. Hypercontractility caused by hypertrophic cardiomyopathy (HCM) myosin mutations is often the result of shifting the equilibrium toward more heads in the on-state. The off-state is equated with a folded-back structure known as the interacting head motif (IHM), which is a regulatory feature of all muscle myosins and class-2 non-muscle myosins. We report here the human β-cardiac myosin IHM structure to 3.6 Å resolution. The structure shows that the interfaces are hot spots of HCM mutations and reveals details of the significant interactions. Importantly, the structures of cardiac and smooth muscle myosin IHMs are dramatically different. This challenges the concept that the IHM structure is conserved in all muscle types and opens new perspectives in the understanding of muscle physiology. The cardiac IHM structure has been the missing puzzle piece to fully understand the development of inherited cardiomyopathies. This work will pave the way for the development of new molecules able to stabilize or destabilize the IHM in a personalized medicine approach. *This manuscript was submitted to Nature Communications in August 2022 and dealt efficiently by the editors. All reviewers received this version of the manuscript before 9 208 August 2022. They also received coordinates and maps of our high resolution structure on the 18 208 August 2022. Due to slowness of at least one reviewer, this contribution was delayed for acceptance by Nature Communications and we are now depositing in bioRxiv the originally submitted version written in July 2022 for everyone to see. Indeed, two bioRxiv contributions at lower resolution but adding similar concepts on thick filament regulation were deposited this week in bioRxiv, one of the contributions having had access to our coordinates. We hope that our data at high resolution will be helpful for all readers that appreciate that high resolution information is required to build accurate atomic models and discuss implications for sarcomere regulation and the effects of cardiomyopathy mutations on heart muscle function.

Muscle Mechanics and Thick Filament Activation: An Emerging Two-Way Interaction for the Vertebrate Striated Muscle Fine Regulation

Contraction in striated muscle is classically described as regulated by calcium-mediated structural changes in the actin-containing thin filaments, which release the binding sites for the interaction with myosin motors to produce force. In this view, myosin motors, arranged in the thick filaments, are basically always ready to interact with the thin filaments, which ultimately regulate the contraction. However, a new “dual-filament” activation paradigm is emerging, where both filaments must be activated to generate force. Growing evidence from the literature shows that the thick filament activation has a role on the striated muscle fine regulation, and its impairment is associated with severe pathologies. This review is focused on the proposed mechanical feedback that activates the inactive motors depending on the level of tension generated by the active ones, the so-called mechanosensing mechanism. Since the main muscle function is to generate mechanical work, the implications on muscle mechanics will be highlighted, showing: (i) how non-mechanical modulation of the thick filament activation influences the contraction, (ii) how the contraction influences the activation of the thick filament and (iii) how muscle, through the mechanical modulation of the thick filament activation, can regulate its own mechanics. This description highlights the crucial role of the emerging bi-directional feedback on muscle mechanical performance.

Conformational changes linked to ADP release from human cardiac myosin bound to actin-tropomyosin

Following binding to the thin filament, β-cardiac myosin couples ATP-hydrolysis to conformational rearrangements in the myosin motor that drive myofilament sliding and cardiac ventricular contraction. However, key features of the cardiac-specific actin-myosin interaction remain uncertain, including the structural effect of ADP release from myosin, which is rate-limiting during force generation. In fact, ADP release slows under experimental load or in the intact heart due to the afterload, thereby adjusting cardiac muscle power output to meet physiological demands. To further elucidate the structural basis of this fundamental process, we used a combination of cryo-EM reconstruction methodologies to determine structures of the human cardiac actin-myosin-tropomyosin filament complex at better than 3.4 Å-resolution in the presence and in the absence of Mg2+·ADP. Focused refinements of the myosin motor head and its essential light chains in these reconstructions reveal that small changes in the nucleotide-binding site are coupled to significant rigid body movements of the myosin converter domain and a 16-degree lever arm swing. Our structures provide a mechanistic framework to understand the effect of ADP binding and release on human cardiac β-myosin, and offer insights into the force-sensing mechanism displayed by the cardiac myosin motor.

Effects of Mavacamten on Measures of Cardiopulmonary Exercise Testing Beyond Peak Oxygen Consumption: A Secondary Analysis of the EXPLORER-HCM Randomized Trial

Importance

Mavacamten, a cardiac myosin inhibitor, improved peak oxygen uptake (pVO2) in patients with symptomatic obstructive hypertrophic cardiomyopathy (HCM) in the EXPLORER-HCM study. However, the full extent of mavacamten's effects on exercise performance remains unclear.

Objective

To investigate the effect of mavacamten on exercise physiology using cardiopulmonary exercise testing (CPET).

Design, Setting, and Participants

Exploratory analyses of the data from the EXPLORER-HCM study, a randomized, double-blind, placebo-controlled, phase 3 trial that was conducted in 68 cardiovascular centers in 13 countries. In total, 251 patients with symptomatic obstructive HCM were enrolled.

Interventions

Patients were randomly assigned in a 1:1 ratio to mavacamten or placebo.

Main Outcomes and Measures

The following prespecified exploratory cardiovascular and performance parameters were assessed with a standardized treadmill or bicycle ergometer test protocol at baseline and week 30: carbon dioxide output (VCO2), minute ventilation (VE), peak VE/VCO2 ratio, ventilatory efficiency (VE/VCO2 slope), peak respiratory exchange ratio (RER), peak circulatory power, ventilatory power, ventilatory threshold, peak metabolic equivalents (METs), peak exercise time, partial pressure of end-tidal carbon dioxide (PETCO2), and VO2/workload slope.

Results

Two hundred fifty-one patients were enrolled. The mean (SD) age was 58.5 (11.9) years and 59% of patients were male. There were significant improvements with mavacamten vs placebo in the following peak-exercise CPET parameters: peak VE/VCO2 ratio (least squares [LS] mean difference, -2.2; 95% CI, -3.05 to -1.26; P < .001), peak METs (LS mean difference, 0.4; 95% CI, 0.17-0.60; P < .001), peak circulatory power (LS mean difference, 372.9 mL/kg/min × mm Hg; 95% CI, 153.12-592.61; P = .001), and peak PETCO2 (LS mean difference, 2.0 mm Hg; 95% CI, 1.12-2.79; P < .001). Mavacamten also improved peak exercise time compared with placebo (LS mean difference, 0.7 minutes; 95% CI, 0.13-1.24; P = .02). There was a significant improvement in nonpeak-exercise CPET parameters, such as VE/VCO2 slope (LS mean difference, -2.6; 95% CI, -3.58 to -1.52; P < .001) and ventilatory power (LS mean difference, 0.6 mm Hg; 95% CI, 0.29-0.90; P < .001) favoring mavacamten vs placebo.

Conclusions and Relevance

Mavacamten improved a range of CPET parameters beyond pVO2, indicating consistent and broad benefits on maximal exercise capacity. Although improvements in peak-exercise CPET parameters are clinically meaningful, the favorable effects of mavacamten on submaximal exertional tolerance provide further insights into the beneficial impact of mavacamten in patients with obstructive HCM.

Trial Registration

ClinicalTrials.gov Identifier: NCT03470545.

Integrating comparative modeling and accelerated simulations reveals conformational and energetic basis of actomyosin force generation

Muscle contraction is performed by arrays of contractile proteins in the sarcomere. Serious heart diseases, such as cardiomyopathy, can often be results of mutations in myosin and actin. Direct characterization of how small changes in the myosin-actin complex impact its force production remains challenging. Molecular dynamics (MD) simulations, although capable of studying protein structure-function relationships, are limited owing to the slow timescale of the myosin cycle as well as a lack of various intermediate structures for the actomyosin complex. Here, employing comparative modeling and enhanced sampling MD simulations, we show how the human cardiac myosin generates force during the mechanochemical cycle. Initial conformational ensembles for different myosin-actin states are learned from multiple structural templates with Rosetta. This enables us to efficiently sample the energy landscape of the system using Gaussian accelerated MD. Key myosin loop residues, whose substitutions are related to cardiomyopathy, are identified to form stable or metastable interactions with the actin surface. We find that the actin-binding cleft closure is allosterically coupled to the myosin motor core transitions and ATP-hydrolysis product release from the active site. Furthermore, a gate between switch I and switch II is suggested to control phosphate release at the prepowerstroke state. Our approach demonstrates the ability to link sequence and structural information to motor functions.

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.

Comparative study of binding pocket structure and dynamics in cardiac and skeletal myosin

The development of small molecule myosin modulators has seen an increased effort in recent years due to their possible use in the treatment of cardiac and skeletal myopathies. Omecamtiv mecarbil (OM) is the first-in-class cardiac myotrope and the first to enter clinical trials. Its selectivity toward slow/beta-cardiac myosin lies at the heart of its function; however, little is known about the underlying reasons for selectivity to this isoform as opposed to other closely related ones such as fast-type skeletal myosins. In this work, we compared the structure and dynamics of the OM binding site in cardiac and in fasttype IIa skeletal myosin to identify possible reasons for OM selectivity. We found that the different shape, size, and composition of the binding pocket in skeletal myosin directly affects the binding mode and related affinity of OM, which is potentially a result of weaker interactions and less optimal molecular recognition. Moreover, we identified a side pocket adjacent to the OM binding site that shows increased accessibility in skeletal myosin compared with the cardiac isoform. These findings could pave the way to the development of skeletal-selective compounds that can target this region of the protein and potentially be used to treat congenital myopathies where muscle weakness is related to myosin loss of function.

Insights into the Mechanism of the Cardiac Drug Omecamtiv Mecarbil─A Computational Study

Omecamtiv mecarbil (OM) is a positive inotrope that is thought to bind directly to an allosteric site of the β-cardiac myosin. The drug is under investigation for the treatment of systolic heart failure. The drug is classified as a cardiac myosin modulator and has been observed to affect multiple vital steps of the cross-bridge cycle including the recovery stroke and the chemical step. We explored the free-energy surface of the recovery stroke of the human cardiac β-myosin in the presence of OM to determine its influence on this process. We also investigated the effects of OM on the recovery stroke in the presence of genetic cardiomyopathic mutations R712L, F764L, and P710R using metadynamics. We also utilized the method of transition path sampling to generate an unbiased ensemble of reactive trajectories for the ATP hydrolysis step in the presence of OM that were able to provide insight into the differences observed due to OM in the dynamics and mechanism of the decomposition of ATP to ADP and HPO₄^2-, a central part of the power generation in cardiac muscle. We studied chemistry in the presence of the same three mutations to further elucidate the effect of OM, and its use in the treatment of cardiac disease.

Dilated cardiomyopathy mutation E525K in human beta-cardiac myosin stabilizes the interacting-heads motif and super-relaxed state of myosin

The auto-inhibited, super-relaxed (SRX) state of cardiac myosin is thought to be crucial for regulating contraction, relaxation, and energy conservation in the heart. We used single ATP turnover experiments to demonstrate that a dilated cardiomyopathy (DCM) mutation (E525K) in human beta-cardiac myosin increases the fraction of myosin heads in the SRX state (with slow ATP turnover), especially in physiological ionic strength conditions. We also utilized FRET between a C-terminal GFP tag on the myosin tail and Cy3ATP bound to the active site of the motor domain to estimate the fraction of heads in the closed, interacting-heads motif (IHM); we found a strong correlation between the IHM and SRX state. Negative stain electron microscopy and 2D class averaging of the construct demonstrated that the E525K mutation increased the fraction of molecules adopting the IHM. Overall, our results demonstrate that the E525K DCM mutation may reduce muscle force and power by stabilizing the auto-inhibited SRX state. Our studies also provide direct evidence for a correlation between the SRX biochemical state and the IHM structural state in cardiac muscle myosin. Furthermore, the E525 residue may be implicated in crucial electrostatic interactions that modulate this conserved, auto-inhibited conformation of myosin.

Hypertrophic cardiomyopathy: Mutations to mechanisms to therapies

Hypertrophic cardiomyopathy (HCM) affects more than 1 in 500 people in the general population with an extensive burden of morbidity in the form of arrhythmia, heart failure, and sudden death. More than 25 years since the discovery of the genetic underpinnings of HCM, the field has unveiled significant insights into the primary effects of these genetic mutations, especially for the myosin heavy chain gene, which is one of the most commonly mutated genes. Our group has studied the molecular effects of HCM mutations on human β-cardiac myosin heavy chain using state-of-the-art biochemical and biophysical tools for the past 10 years, combining insights from clinical genetics and structural analyses of cardiac myosin. The overarching hypothesis is that HCM-causing mutations in sarcomere proteins cause hypercontractility at the sarcomere level, and we have shown that an increase in the number of myosin molecules available for interaction with actin is a primary driver. Recently, two pharmaceutical companies have developed small molecule inhibitors of human cardiac myosin to counteract the molecular consequences of HCM pathogenesis. One of these inhibitors (mavacamten) has recently been approved by the FDA after completing a successful phase III trial in HCM patients, and the other (aficamten) is currently being evaluated in a phase III trial. Myosin inhibitors will be the first class of medication used to treat HCM that has both robust clinical trial evidence of efficacy and that targets the fundamental mechanism of HCM pathogenesis. The success of myosin inhibitors in HCM opens the door to finding other new drugs that target the sarcomere directly, as we learn more about the genetics and fundamental mechanisms of this disease.

Myosin Transducer Inter-Strand Communication Is Critical for Normal ATPase Activity and Myofibril Structure

The R249Q mutation in human β-cardiac myosin results in hypertrophic cardiomyopathy. We previously showed that inserting this mutation into Drosophila melanogaster indirect flight muscle myosin yields mechanical and locomotory defects. Here, we use transgenic Drosophila mutants to demonstrate that residue R249 serves as a critical communication link within myosin that controls both ATPase activity and myofibril integrity. R249 is located on a β-strand of the central transducer of myosin, and our molecular modeling shows that it interacts via a salt bridge with D262 on the adjacent β-strand. We find that disrupting this interaction via R249Q, R249D or D262R mutations reduces basal and actin-activated ATPase activity, actin in vitro motility and flight muscle function. Further, the R249D mutation dramatically affects myofibril assembly, yielding abnormalities in sarcomere lengths, increased Z-line thickness and split myofibrils. These defects are exacerbated during aging. Re-establishing the β-strand interaction via a R249D/D262R double mutation restores both basal ATPase activity and myofibril assembly, indicating that these properties are dependent upon transducer inter-strand communication. Thus, the transducer plays an important role in myosin function and myofibril architecture.

Hypertrophic cardiomyopathy mutations in the pliant and light chain-binding regions of the lever arm of human β-cardiac myosin have divergent effects on myosin function

Mutations in the lever arm of β-cardiac myosin are a frequent cause of hypertrophic cardiomyopathy, a disease characterized by hypercontractility and eventual hypertrophy of the left ventricle. Here, we studied five such mutations: three in the pliant region of the lever arm (D778V, L781P, and S782N) and two in the light chain-binding region (A797T and F834L). We investigated their effects on both motor function and myosin subfragment 2 (S2) tail-based autoinhibition. The pliant region mutations had varying effects on the motor function of a myosin construct lacking the S2 tail: overall, D778V increased power output, L781P reduced power output, and S782N had little effect on power output, while all three reduced the external force sensitivity of the actin detachment rate. With a myosin containing the motor domain and the proximal S2 tail, the pliant region mutations also attenuated autoinhibition in the presence of filamentous actin but had no impact in the absence of actin. By contrast, the light chain-binding region mutations had little effect on motor activity but produced marked reductions in autoinhibition in both the presence and absence of actin. Thus, mutations in the lever arm of β-cardiac myosin have divergent allosteric effects on myosin function, depending on whether they are in the pliant or light chain-binding regions.

Targeting the sarcomere in inherited cardiomyopathies

Variants in >12 genes encoding sarcomeric proteins can cause various cardiomyopathies. The two most common are hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM). Current therapeutics do not target the root causes of these diseases, but attempt to prevent disease progression and/or to manage symptoms. Accordingly, novel approaches are being developed to treat the cardiac muscle dysfunction directly. Challenges to developing therapeutics for these diseases include the diverse mechanisms of pathogenesis, some of which are still being debated and defined. Four small molecules that modulate the myosin motor protein in the cardiac sarcomere have shown great promise in the settings of HCM and DCM, regardless of the underlying genetic pathogenesis, and similar approaches are being developed to target other components of the sarcomere. In the setting of HCM, mavacamten and aficamten bind to the myosin motor and decrease the ATPase activity of myosin. In the setting of DCM, omecamtiv mecarbil and danicamtiv increase myosin activity in cardiac muscle (but omecamtiv mecarbil decreases myosin activity in vitro). In this Review, we discuss the therapeutic strategies to alter sarcomere contractile activity and summarize the data indicating that targeting one protein in the sarcomere can be effective in treating patients with genetic variants in other sarcomeric proteins, as well as in patients with non-sarcomere-based disease.

Critical Evaluation of Current Hypotheses for the Pathogenesis of Hypertrophic Cardiomyopathy

Hereditary hypertrophic cardiomyopathy (HCM), due to mutations in sarcomere proteins, occurs in more than 1/500 individuals and is the leading cause of sudden cardiac death in young people. The clinical course exhibits appreciable variability. However, typically, heart morphology and function are normal at birth, with pathological remodeling developing over years to decades, leading to a phenotype characterized by asymmetric ventricular hypertrophy, scattered fibrosis and myofibrillar/cellular disarray with ultimate mechanical heart failure and/or severe arrhythmias. The identity of the primary mutation-induced changes in sarcomere function and how they trigger debilitating remodeling are poorly understood. Support for the importance of mutation-induced hypercontractility, e.g., increased calcium sensitivity and/or increased power output, has been strengthened in recent years. However, other ideas that mutation-induced hypocontractility or non-uniformities with contractile instabilities, instead, constitute primary triggers cannot yet be discarded. Here, we review evidence for and criticism against the mentioned hypotheses. In this process, we find support for previous ideas that inefficient energy usage and a blunted Frank-Starling mechanism have central roles in pathogenesis, although presumably representing effects secondary to the primary mutation-induced changes. While first trying to reconcile apparently diverging evidence for the different hypotheses in one unified model, we also identify key remaining questions and suggest how experimental systems that are built around isolated primarily expressed proteins could be useful.

Molecular mechanisms of cardiac actomyosin transforming from rigor state to post-rigor state

Sudden cardiac death contributed to half of all deaths from cardiovascular diseases. The mechanism of the kinetic cycle of cardiac myosin is crucial for heart protection and drug development. The state change in the myosin kinetic cycle from the rigor state to the post-rigor state is fundamental to explain binding and dissociation. Here, we used β-cardiac myosin in the rigor and post-rigor states to model the actomyosin complexes. Molecular dynamics simulations, electrostatic analysis, and energetic analysis of actomyosin complexes were performed in this work. The results showed that there are fewer interactions and lower electrostatic binding strength in the post-rigor state than in the rigor state. In the post-rigor state, there were higher free binding energy, fewer salt bridges, and fewer hydrogen bonds. The results showed a lower binding affinity in the post-rigor state than in the rigor state. The decrease in the binding affinity provided important conditions for dissociation of the myosin from the actin filament. Although previous studies focused mostly on the binding process, this study provides evidence of dissociation, which is even more important in the myosin kinetic cycle. This research on the mechanism of myosin kinetic cycles provides a novel direction for future genetic disease studies.

Cryo-EM structure of the autoinhibited state of myosin-2

We solved the near-atomic resolution structure of smooth muscle myosin-2 in the autoinhibited state (10S) using single-particle cryo–electron microscopy. The 3.4-Å structure reveals the precise molecular architecture of 10S and the structural basis for myosin-2 regulation. We reveal the position of the phosphorylation sites that control myosin autoinhibition and activation by phosphorylation of the regulatory light chain. Further, we present a previously unidentified conformational state in myosin-2 that traps ADP and Pi produced by the hydrolysis of ATP in the active site. This noncanonical state represents a branch of the myosin enzyme cycle and explains the autoinhibition of the enzyme function of 10S along with its reduced affinity for actin. Together, our structure defines the molecular mechanisms that drive 10S formation, stabilization, and relief by phosphorylation of the regulatory light chain.

Two Classes of Myosin Inhibitors, Para-nitroblebbistatin and Mavacamten, Stabilize β-Cardiac Myosin in Different Structural and Functional States

In addition to a conventional relaxed state, a fraction of myosins in the cardiac muscle exists in a low-energy consuming super-relaxed (SRX) state, which is kept as a reserve pool that may be engaged under sustained increased cardiac demand. The conventional relaxed and the super-relaxed states are widely assumed to correspond to a structure where myosin heads are in an open configuration, free to interact with actin, and a closed configuration, inhibiting binding to actin, respectively. Disruption of the myosin SRX population is an emerging model in different heart diseases, such as hypertrophic cardiomyopathy, which results in excessive muscle contraction, and stabilizing them using myosin inhibitors is budding as an attractive therapeutic strategy. Here we examined the structure-function relationships of two myosin ATPase inhibitors, mavacamten and para-nitroblebbistatin, and found that binding of mavacamten at a site different than para-nitroblebbistatin populates myosin into the SRX state. Para-nitroblebbistatin, binding to a distal pocket to the myosin lever arm near the nucleotide-binding site, does not affect the usual myosin SRX state but instead appears to render myosin into a new, perhaps much more inhibited, 'ultra-relaxed' state. X-ray scattering-based rigid body modeling shows that both mavacamten and para-nitroblebbistatin induce novel conformations in human β-cardiac heavy meromyosin that diverge significantly from the hypothetical open and closed states, and furthermore, mavacamten treatment causes greater compaction than para-nitroblebbistatin. Taken together, we conclude that mavacamten and para-nitroblebbistatin stabilize myosin in different structural states, and such states may give rise to different functional energy-sparing states.

Cardiac myosin contraction and mechanotransduction in health and disease

Cardiac myosin is the molecular motor that powers heart contraction by converting chemical energy from ATP hydrolysis into mechanical force. The power output of the heart is tightly regulated to meet the physiological needs of the body. Recent multiscale studies spanning from molecules to tissues have revealed complex regulatory mechanisms that fine-tune cardiac contraction, in which myosin not only generates power output but also plays an active role in its regulation. Thus, myosin is both shaped by and actively involved in shaping its mechanical environment. Moreover, these studies have shown that cardiac myosin-generated tension affects physiological processes beyond muscle contraction. Here, we review these novel regulatory mechanisms, as well as the roles that myosin-based force generation and mechanotransduction play in development and disease. We describe how key intra- and intermolecular interactions contribute to the regulation of myosin-based contractility and the role of mechanical forces in tuning myosin function. We also discuss the emergence of cardiac myosin as a drug target for diseases including heart failure, leading to the discovery of therapeutics that directly tune myosin contractility. Finally, we highlight some of the outstanding questions that must be addressed to better understand myosin's functions and regulation, and we discuss prospects for translating these discoveries into precision medicine therapeutics targeting contractility and mechanotransduction.

Direct detection of the myosin super-relaxed state and interacting-heads motif in solution

The interacting-heads motif (IHM) is a structure of myosin that has been proposed to modulate cardiac output by occluding myosin molecules from undergoing the force-generating cycle. It is hypothesized to be the structural basis for the super-relaxed state (SRX), a low-ATPase kinetic state thought to be cardioprotective. The goal of the present study was to test this hypothesis by determining directly and quantitatively the fractions of myosin in the IHM and SRX under the same conditions in solution. To detect the structural IHM, we used time-resolved fluorescence resonance energy transfer to quantitate two distinct populations. One population was observed at a center distance of 2.0 nm, whereas the other was not detectable by fluorescence resonance energy transfer, implying a distance greater than 4 nm. We confirmed the IHM assignment to the 2.0-nm population by applying the same cross-linking protocol used previously to image the IHM by electron microscopy. Under the same conditions, we also measured the fraction of myosin in the SRX using stopped-flow kinetics. Our results show that the populations of SRX and IHM myosin were similar, unless treated with mavacamten, a drug that recently completed phase III clinical trials to treat hypertrophic cardiomyopathy and is proposed to act by stabilizing both the SRX and IHM. However, we found that mavacamten had a much greater effect on the SRX (55% increase) than on the IHM (4% increase). We conclude that the IHM structure is sufficient but not necessary to produce the SRX kinetic state.

Hypertrophic cardiomyopathy β-cardiac myosin mutation (P710R) leads to hypercontractility by disrupting super relaxed state

Hypertrophic cardiomyopathy (HCM) is the most common inherited form of heart disease, associated with over 1,000 mutations, many in β-cardiac myosin (MYH7). Molecular studies of myosin with different HCM mutations have revealed a diversity of effects on ATPase and load-sensitive rate of detachment from actin. It has been difficult to predict how such diverse molecular effects combine to influence forces at the cellular level and further influence cellular phenotypes. This study focused on the P710R mutation that dramatically decreased in vitro motility velocity and actin-activated ATPase, in contrast to other MYH7 mutations. Optical trap measurements of single myosin molecules revealed that this mutation reduced the step size of the myosin motor and the load sensitivity of the actin detachment rate. Conversely, this mutation destabilized the super relaxed state in longer, two-headed myosin constructs, freeing more heads to generate force. Micropatterned human induced pluripotent derived stem cell (hiPSC)-cardiomyocytes CRISPR-edited with the P710R mutation produced significantly increased force (measured by traction force microscopy) compared with isogenic control cells. The P710R mutation also caused cardiomyocyte hypertrophy and cytoskeletal remodeling as measured by immunostaining and electron microscopy. Cellular hypertrophy was prevented in the P710R cells by inhibition of ERK or Akt. Finally, we used a computational model that integrated the measured molecular changes to predict the measured traction forces. These results confirm a key role for regulation of the super relaxed state in driving hypercontractility in HCM with the P710R mutation and demonstrate the value of a multiscale approach in revealing key mechanisms of disease.

Hypertrophic cardiomyopathy β-cardiac myosin mutation (P710R) leads to hypercontractility by disrupting super relaxed state

Hypertrophic cardiomyopathy (HCM) is the most common inherited form of heart disease, associated with over 1,000 mutations, many in β-cardiac myosin (MYH7). Molecular studies of myosin with different HCM mutations have revealed a diversity of effects on ATPase and load-sensitive rate of detachment from actin. It has been difficult to predict how such diverse molecular effects combine to influence forces at the cellular level and further influence cellular phenotypes. This study focused on the P710R mutation that dramatically decreased in vitro motility velocity and actin-activated ATPase, in contrast to other MYH7 mutations. Optical trap measurements of single myosin molecules revealed that this mutation reduced the step size of the myosin motor and the load sensitivity of the actin detachment rate. Conversely, this mutation destabilized the super relaxed state in longer, two-headed myosin constructs, freeing more heads to generate force. Micropatterned human induced pluripotent derived stem cell (hiPSC)-cardiomyocytes CRISPR-edited with the P710R mutation produced significantly increased force (measured by traction force microscopy) compared with isogenic control cells. The P710R mutation also caused cardiomyocyte hypertrophy and cytoskeletal remodeling as measured by immunostaining and electron microscopy. Cellular hypertrophy was prevented in the P710R cells by inhibition of ERK or Akt. Finally, we used a computational model that integrated the measured molecular changes to predict the measured traction forces. These results confirm a key role for regulation of the super relaxed state in driving hypercontractility in HCM with the P710R mutation and demonstrate the value of a multiscale approach in revealing key mechanisms of disease.

Investigation of the Recovery Stroke and ATP Hydrolysis and Changes Caused Due to the Cardiomyopathic Point Mutations in Human Cardiac β Myosin

Human cardiac β myosin undergoes the cross-bridge cycle as part of the force-generating mechanism of cardiac muscle. The recovery stroke is considered one of the key steps of the kinetic cycle as it is the conformational rearrangement required to position the active site residues for hydrolysis of ATP and interaction with actin. We explored the free-energy surface of the transition and investigated the effect of the genetic cardiomyopathy causing mutations R453C, I457T, and I467T on this step using metadynamics. This work extends previous studies on Dictyostelium myosin II with engineered mutations. Here, like previously, we generated an unbiased thermodynamic ensemble of reactive trajectories for the chemical step using transition path sampling. Our methodologies were able to predict the changes to the dynamics of the recovery stroke as well as predict the pathway of breakdown of ATP to ADP and HPO₄^2- with the stabilization of the metaphosphate intermediate. We also observed clear differences between the Dictyostelium myosin II and human cardiac β myosin for ATP hydrolysis as well as predict the effect of the mutation I467T on the chemical step.

Identification of sequence changes in myosin II that adjust muscle contraction velocity

The speed of muscle contraction is related to body size; muscles in larger species contract at slower rates. Since contraction speed is a property of the myosin isoform expressed in a muscle, we investigated how sequence changes in a range of muscle myosin II isoforms enable this slower rate of muscle contraction. We considered 798 sequences from 13 mammalian myosin II isoforms to identify any adaptation to increasing body mass. We identified a correlation between body mass and sequence divergence for the motor domain of the 4 major adult myosin II isoforms (β/Type I, IIa, IIb, and IIx), suggesting that these isoforms have adapted to increasing body mass. In contrast, the non-muscle and developmental isoforms show no correlation of sequence divergence with body mass. Analysis of the motor domain sequence of β-myosin (predominant myosin in Type I/slow and cardiac muscle) from 67 mammals from 2 distinct clades identifies 16 sites, out of 800, associated with body mass (padj < 0.05) but not with the clade (padj > 0.05). Both clades change the same small set of amino acids, in the same order from small to large mammals, suggesting a limited number of ways in which contraction velocity can be successfully manipulated. To test this relationship, the 9 sites that differ between human and rat were mutated in the human β-myosin to match the rat sequence. Biochemical analysis revealed that the rat-human β-myosin chimera functioned like the native rat myosin with a 2-fold increase in both motility and in the rate of ADP release from the actin-myosin crossbridge (the step that limits contraction velocity). Thus, these sequence changes indicate adaptation of β-myosin as species mass increased to enable a reduced contraction velocity and heart rate.

Cardiomyopathic mutations in essential light chain reveal mechanisms regulating the super relaxed state of myosin

In this study, we assessed the super relaxed (SRX) state of myosin and sarcomeric protein phosphorylation in two pathological models of cardiomyopathy and in a near-physiological model of cardiac hypertrophy. The cardiomyopathy models differ in disease progression and severity and express the hypertrophic (HCM-A57G) or restrictive (RCM-E143K) mutations in the human ventricular myosin essential light chain (ELC), which is encoded by the MYL3 gene. Their effects were compared with near-physiological heart remodeling, represented by the N-terminally truncated ELC (Δ43 ELC mice), and with nonmutated human ventricular WT-ELC mice. The HCM-A57G and RCM-E143K mutations had antagonistic effects on the ATP-dependent myosin energetic states, with HCM-A57G cross-bridges fostering the disordered relaxed (DRX) state and the RCM-E143K model favoring the energy-conserving SRX state. The HCM-A57G model promoted the switch from the SRX to DRX state and showed an ∼40% increase in myosin regulatory light chain (RLC) phosphorylation compared with the RLC of normal WT-ELC myocardium. On the contrary, the RCM-E143K–associated stabilization of the SRX state was accompanied by an approximately twofold lower level of myosin RLC phosphorylation compared with the RLC of WT-ELC. Upregulation of RLC phosphorylation was also observed in Δ43 versus WT-ELC hearts, and the Δ43 myosin favored the energy-saving SRX conformation. The two disease variants also differently affected the duration of force transients, with shorter (HCM-A57G) or longer (RCM-E143K) transients measured in electrically stimulated papillary muscles from these pathological models, while no changes were displayed by Δ43 fibers. We propose that the N terminus of ELC (N-ELC), which is missing in the hearts of Δ43 mice, works as an energetic switch promoting the SRX-to-DRX transition and contributing to the regulation of myosin RLC phosphorylation in full-length ELC mice by facilitating or sterically blocking RLC phosphorylation in HCM-A57G and RCM-E143K hearts, respectively.

Filamentous tangles with nemaline rods in MYH2 myopathy: a novel phenotype

The MYH2 gene encodes the skeletal muscle myosin heavy chain IIA (MyHC-IIA) isoform, which is expressed in the fast twitch type 2A fibers. Autosomal dominant or recessive pathogenic variants in MYH2 lead to congenital myopathy clinically featured by ophthalmoparesis and predominantly proximal weakness. MYH2-myopathy is pathologically characterized by loss and atrophy of type 2A fibers. Additional myopathological abnormalities have included rimmed vacuoles containing small p62 positive inclusions, 15–20 nm tubulofilaments, minicores and dystrophic changes. We report an adult patient with late-pediatric onset MYH2-myopathy caused by two heterozygous pathogenic variants: c.3331C>T, p.Gln1111* predicted to result in truncation of the proximal tail region of MyHC-IIA, and c.1546T>G, p.Phe516Val, affecting a highly conserved amino acid within the highly conserved catalytic motor head relay loop. This missense variant is predicted to result in a less compact loop domain and in turn could affect the protein affinity state. The patient’s genotype is accompanied by a novel myopathological phenotype characterized by centralized large myofilamentous tangles associated with clusters of nemaline rods, and ring fibers, in addition to the previously reported rimmed vacuoles, paucity and atrophy of type 2A fibers. Electron microscopy demonstrated wide areas of disorganized myofibrils which were oriented in various planes of direction and entrapped multiple nemaline rods, as corresponding to the large tangles with rods seen on light microscopy. Nemaline rods were rarely observed also in nuclei. We speculate that the mutated MyHC-IIA may influence myofibril disorganization. While nemaline rods have been described in myopathies caused by pathogenic variants in genes encoding several sarcomeric proteins, to our knowledge, nemaline rods have not been previously described in MYH2-myopathy.

MYH1 is a candidate gene for recurrent rhabdomyolysis in humans

Rhabdomyolysis is a serious medical condition characterized by muscle injury, and there are recognized genetic causes especially in recurrent forms. The majority of these cases, however, remain unexplained. Here, we describe a patient with recurrent rhabdomyolysis in whom extensive clinical testing failed to identify a likely etiology. Whole-exome sequencing revealed a novel missense variant in MYH1, which encodes a major adult muscle fiber protein. Structural biology analysis revealed that the mutated residue is extremely well conserved and is located in the actin binding cleft. Furthermore, immediately adjacent mutations in that cleft in other myosins are pathogenic in humans. Our results are consistent with the finding that MYH1 is mutated in rhabdomyolysis in horses and suggest that this gene should be investigated in cases with recurrent rhabdomyolysis.

Myosin with hypertrophic cardiac mutation R712L has a decreased working stroke which is rescued by omecamtiv mecarbil

Hypertrophic cardiomyopathies (HCMs) are the leading cause of acute cardiac failure in young individuals. Over 300 mutations throughout β-cardiac myosin, including in the motor domain, are associated with HCM. A β-cardiac myosin motor mutation (R712L) leads to a severe form of HCM. Actin-gliding motility of R712L-myosin is inhibited, despite near-normal ATPase kinetics. By optical trapping, the working stroke of R712L-myosin was decreased 4-fold, but actin-attachment durations were normal. A prevalent hypothesis that HCM mutants are hypercontractile is thus not universal. R712 is adjacent to the binding site of the heart failure drug omecamtiv mecarbil (OM). OM suppresses the working stroke of normal β-cardiac myosin, but remarkably, OM rescues the R712L-myosin working stroke. Using a flow chamber to interrogate a single molecule during buffer exchange, we found OM rescue to be reversible. Thus, the R712L mutation uncouples lever arm rotation from ATPase activity and this inhibition is rescued by OM.

Cardiac myosin super relaxation (SRX): a perspective on fundamental biology, human disease and therapeutics

The fundamental basis of muscle contraction 'the sliding filament model' (Huxley and Niedergerke, 1954; Huxley and Hanson, 1954) and the 'swinging, tilting crossbridge-sliding filament mechanism' (Huxley, 1969; Huxley and Brown, 1967) nucleated a field of research that has unearthed the complex and fascinating role of myosin structure in the regulation of contraction. A recently discovered energy conserving state of myosin termed the super relaxed state (SRX) has been observed in filamentous myosins and is central to modulating force production and energy use within the sarcomere. Modulation of myosin function through SRX is a rapidly developing theme in therapeutic development for both cardiovascular disease and infectious disease. Some 70 years after the first discoveries concerning muscular function, modulation of myosin SRX may bring the first myosin targeted small molecule to the clinic, for treating hypertrophic cardiomyopathy (Olivotto et al., 2020). An often monogenic disease HCM afflicts 1 in 500 individuals, and can cause heart failure and sudden cardiac death. Even as we near therapeutic translation, there remain many questions about the governance of muscle function in human health and disease. With this review, we provide a broad overview of contemporary understanding of myosin SRX, and explore the complexities of targeting this myosin state in human disease.This article has an associated Future Leaders to Watch interview with the authors of the paper.

Hypertrophic cardiomyopathy associated E22K mutation in myosin regulatory light chain decreases calcium-activated tension and stiffness and reduces myofilament Ca2+ sensitivity

We investigated the mechanisms associated with E22K mutation in myosin regulatory light chain (RLC), found to cause hypertrophic cardiomyopathy (HCM) in humans and mice. Specifically, we characterized the mechanical profiles of papillary muscle fibers from transgenic mice expressing human ventricular RLC wild-type (Tg-WT) or E22K mutation (Tg-E22K). Because the two mouse models expressed different amounts of transgene, the B6SJL mouse line (NTg) was used as an additional control. Mechanical experiments were carried out on Ca²⁺ - and ATP-activated fibers and in rigor. Sinusoidal analysis was performed to elucidate the effect of E22K on tension and stiffness during activation/rigor, tension-pCa, and myosin cross-bridge (CB) kinetics. We found significant reductions in active tension (by 54%) and stiffness (active by 40% and rigor by 54%). A decrease in the Ca²⁺ sensitivity of tension (by ∆pCa ~ 0.1) was observed in Tg-E22K compared with Tg-WT fibers. The apparent (=measured) rate constant of exponential process B (2πb: force generation step) was not affected by E22K, but the apparent rate constant of exponential process C (2πc: CB detachment step) was faster in Tg-E22K compared with Tg-WT fibers. Both 2πb and 2πc were smaller in NTg than in Tg-WT fibers, suggesting a kinetic difference between the human and mouse RLC. Our results of E22K-induced reduction in myofilament stiffness and tension suggest that the main effect of this mutation was to disturb the interaction of RLC with the myosin heavy chain and impose structural abnormalities in the lever arm of myosin CB. When placed in vivo, the E22K mutation is expected to result in reduced contractility and decreased cardiac output whereby leading to HCM.

SUB-DISCIPLINE

Bioenergetics.

DATABASE

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

ANIMAL PROTOCOL

BK20150353 (Soochow University).

RESEARCH GOVERNANCE

School of Nursing: Hua-Gang Hu: seuboyh@163.com; Soochow University: Chen Ge chge@suda.edu.cn.

Hybrid method for representing ions in implicit solvation calculations

Fast and accurate calculations of the electrostatic features of highly charged biomolecules such as DNA, RNA, and highly charged proteins are crucial and challenging tasks. Traditional implicit solvent methods calculate the electrostatic features quickly, but these methods are not able to balance the high net biomolecular charges effectively. Explicit solvent methods add unbalanced ions to neutralize the highly charged biomolecules in molecular dynamic simulations, which require more expensive computing resources. Here we report developing a novel method, Hybridizing Ions Treatment (HIT), which hybridizes the implicit solvent method with an explicit method to realistically calculate the electrostatic potential for highly charged biomolecules. HIT utilizes the ionic distribution from an explicit method to predict the bound ions. The bound ions are then added in the implicit solvent method to perform the electrostatic potential calculations. In this study, two training sets were developed to optimize parameters for HIT. The performance on the testing set demonstrates that HIT significantly improves the electrostatic calculations. Results on molecular motors myosin and kinesin reveal some mechanisms and explain some previous experimental findings. HIT can be widely used to study highly charged biomolecules, including DNA, RNA, molecular motors, and other highly charged biomolecules. The HIT package is available at http://compbio.utep.edu/static/downloads/download_hit.zip.

High-Resolution Cryo-EM Structure of the Cardiac Actomyosin Complex

Heart contraction depends on a complicated array of interactions between sarcomeric proteins required to convert chemical energy into mechanical force. Cyclic interactions between actin and myosin molecules, controlled by troponin and tropomyosin, generate the sliding force between the actin-based thin and myosin-based thick filaments. Alterations in this sophisticated system due to missense mutations can lead to cardiovascular diseases. Numerous structural studies proposed pathological mechanisms of missense mutations at the myosin-myosin, actin-tropomyosin, and tropomyosin-troponin interfaces. However, despite the central role of actomyosin interactions a detailed structural description of the cardiac actomyosin interface remained unknown. Here, we report a cryo-EM structure of a cardiac actomyosin complex at 3.8 Å resolution. The structure reveals the molecular basis of cardiac diseases caused by missense mutations in myosin and actin proteins.

Synthetic thick filaments: A new avenue for better understanding the myosin super-relaxed state in healthy, diseased, and mavacamten-treated cardiac systems

A hallmark feature of myosin-II is that it can spontaneously self-assemble into bipolar synthetic thick filaments (STFs) in low-ionic-strength buffers, thereby serving as a reconstituted in vitro model for muscle thick filaments. Although these STFs have been extensively used for structural characterization, their functional evaluation has been limited. In this report, we show that myosins in STFs mirror the more electrostatic and cooperative interactions that underlie the energy-sparing super-relaxed (SRX) state, which are not seen using shorter myosin subfragments, heavy meromyosin (HMM) and myosin subfragment 1 (S1). Using these STFs, we show several pathophysiological insults in hypertrophic cardiomyopathy, including the R403Q myosin mutation, phosphorylation of myosin light chains, and an increased ADP:ATP ratio, destabilize the SRX population. Furthermore, WT myosin containing STFs, but not S1, HMM, or STFs-containing R403Q myosin, recapitulated the ADP-induced destabilization of the SRX state. Studies involving a clinical-stage small-molecule inhibitor, mavacamten, showed that it is more effective in not only increasing myosin SRX population in STFs than in S1 or HMM but also in increasing myosin SRX population equally well in STFs made of healthy and disease-causing R403Q myosin. Importantly, we also found that pathophysiological perturbations such as elevated ADP concentration weakens mavacamten's ability to increase the myosin SRX population, suggesting that mavacamten-bound myosin heads are not permanently protected in the SRX state but can be recruited into action. These findings collectively emphasize that STFs serve as a valuable tool to provide novel insights into the myosin SRX state in healthy, diseased, and therapeutic conditions.

Alpha and beta myosin isoforms and human atrial and ventricular contraction

Human atrial and ventricular contractions have distinct mechanical characteristics including speed of contraction, volume of blood delivered and the range of pressure generated. Notably, the ventricle expresses predominantly β-cardiac myosin while the atrium expresses mostly the α-isoform. In recent years exploration of the properties of pure α- & β-myosin isoforms have been possible in solution, in isolated myocytes and myofibrils. This allows us to consider the extent to which the atrial vs ventricular mechanical characteristics are defined by the myosin isoform expressed, and how the isoform properties are matched to their physiological roles. To do this we Outline the essential feature of atrial and ventricular contraction; Explore the molecular structural and functional characteristics of the two myosin isoforms; Describe the contractile behaviour of myocytes and myofibrils expressing a single myosin isoform; Finally we outline the outstanding problems in defining the differences between the atria and ventricles. This allowed us consider what features of contraction can and cannot be ascribed to the myosin isoforms present in the atria and ventricles.

Heart Failure Drug Modifies the Intrinsic Dynamics of the Pre-Power Stroke State of Cardiac Myosin

Omecamtiv mecarbil (OM), currently investigated for the treatment of heart failure, is the first example of a new class of drugs (cardiac myotropes) that can modify muscle contractility by directly targeting sarcomeric proteins. Using atomistic molecular dynamics simulations, we show that the binding of OM to the pre-power stroke state of cardiac myosin inhibits the functional motions of the protein and potentially affects Pi release from the nucleotide binding site. We also show that the changes in myosin ATPase activity induced by a set of OM analogues can be predicted from their relative affinity to the pre-power stroke state compared to the near rigor one, indicating that conformational selectivity plays an important role in determining the activity of these compounds.