Structure of the native myosin filament in the relaxed cardiac sarcomere.

The thick filament is a key component of sarcomeres, the basic units of striated muscle1. Alterations in thick filament proteins are associated with familial hypertrophic cardiomyopathy and other heart and muscle diseases2. Despite the central importance of the thick filament, its molecular organization remains unclear. Here we present the molecular architecture of native cardiac sarcomeres in the relaxed state, determined by cryo-electron tomography. Our reconstruction of the thick filament reveals the three-dimensional organization of myosin, titin and myosin-binding protein C (MyBP-C). The arrangement of myosin molecules is dependent on their position along the filament, suggesting specialized capacities in terms of strain susceptibility and force generation. Three pairs of titin-α and titin-ß chains run axially along the filament, intertwining with myosin tails and probably orchestrating the length-dependent activation of the sarcomere. Notably, whereas the three titin-α chains run along the entire length of the thick filament, titin-ß chains do not. The structure also demonstrates that MyBP-C bridges thin and thick filaments, with its carboxy-terminal region binding to the myosin tails and directly stabilizing the OFF state of the myosin heads in an unforeseen manner. These results provide a foundation for future research investigating muscle disorders involving sarcomeric components.

Cryo-EM structure of the human cardiac myosin filament.

Pumping of the heart is powered by filaments of the motor protein myosin that pull on actin filaments to generate cardiac contraction. In addition to myosin, the filaments contain cardiac myosin-binding protein C (cMyBP-C), which modulates contractility in response to physiological stimuli, and titin, which functions as a scaffold for filament assembly1. Myosin, cMyBP-C and titin are all subject to mutation, which can lead to heart failure. Despite the central importance of cardiac myosin filaments to life, their molecular structure has remained a mystery for 60 years2. Here we solve the structure of the main (cMyBP-C-containing) region of the human cardiac filament using cryo-electron microscopy. The reconstruction reveals the architecture of titin and cMyBP-C and shows how myosin's motor domains (heads) form three different types of motif (providing functional flexibility), which interact with each other and with titin and cMyBP-C to dictate filament architecture and function. The packing of myosin tails in the filament backbone is also resolved. The structure suggests how cMyBP-C helps to generate the cardiac super-relaxed state3; how titin and cMyBP-C may contribute to length-dependent activation4; and how mutations in myosin and cMyBP-C might disturb interactions, causing disease5,6. The reconstruction resolves past uncertainties and integrates previous data on cardiac muscle structure and function. It provides a new paradigm for interpreting structural, physiological and clinical observations, and for the design of potential therapeutic drugs.

Binding pocket dynamics along the recovery stroke of human ß-cardiac myosin.

The druggability of small-molecule binding sites can be significantly affected by protein motions and conformational changes. Ligand binding, protein dynamics and protein function have been shown to be closely interconnected in myosins. The breakthrough discovery of omecamtiv mecarbil (OM) has led to an increased interest in small molecules that can target myosin and modulate its function for therapeutic purposes (myosin modulators). In this work, we use a combination of computational methods, including steered molecular dynamics, umbrella sampling and binding pocket tracking tools, to follow the evolution of the OM binding site during the recovery stroke transition of human ß-cardiac myosin. We found that steering two internal coordinates of the motor domain can recapture the main features of the transition and in particular the rearrangements of the binding site, which shows significant changes in size, shape and composition. Possible intermediate conformations were also identified, in remarkable agreement with experimental findings. The differences in the binding site properties observed along the transition can be exploited for the future development of conformation-selective myosin modulators.

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.

Myosin 7b is a regulatory long noncoding RNA (lncMYH7b) in the human heart.

Myosin heavy chain 7b (MYH7b) is an ancient member of the myosin heavy chain motor protein family that is expressed in striated muscles. In mammalian cardiac muscle, MYH7b RNA is expressed along with two other myosin heavy chains, ß-myosin heavy chain (ß-MyHC) and α-myosin heavy chain (α-MyHC). However, unlike ß-MyHC and α-MyHC, which are maintained in a careful balance at the protein level, the MYH7b locus does not produce a full-length protein in the heart due to a posttranscriptional exon-skipping mechanism that occurs in a tissue-specific manner. Whether this locus has a role in the heart beyond producing its intronic microRNA, miR-499, was unclear. Using cardiomyocytes derived from human induced pluripotent stem cells as a model system, we found that the noncoding exon-skipped RNA (lncMYH7b) affects the transcriptional landscape of human cardiomyocytes, independent of miR-499. Specifically, lncMYH7b regulates the ratio of ß-MyHC to α-MyHC, which is crucial for cardiac contractility. We also found that lncMYH7b regulates beat rate and sarcomere formation in cardiomyocytes. This regulation is likely achieved through control of a member of the TEA domain transcription factor family (TEAD3, which is known to regulate ß-MyHC). Therefore, we conclude that this ancient gene has been repurposed by alternative splicing to produce a regulatory long-noncoding RNA in the human heart that affects cardiac myosin composition.

The Local Environment of Loop Switch 1 Modulates the Rate of ATP-Induced Dissociation of Human Cardiac Actomyosin.

Two isoforms of human cardiac myosin, alpha and beta, share significant sequence similarities but show different kinetics. The alpha isoform is a faster motor; it spends less time being strongly bound to actin during the actomyosin cycle. With alpha isoform, actomyosin dissociates faster upon ATP binding, and the affinity of ADP to actomyosin is weaker. One can suggest that the isoform-specific actomyosin kinetics is regulated at the nucleotide binding site of human cardiac myosin. Myosin is a P-loop ATPase; the nucleotide-binding site consists of P-loop and loops switch 1 and 2. All three loops position MgATP for successful hydrolysis. Loops sequence is conserved in both myosin isoforms, and we hypothesize that the isoform-specific structural element near the active site regulates the rate of nucleotide binding and release. Previously we ran molecular dynamics simulations and found that loop S291-E317 near loop switch 1 is more compact and exhibits larger fluctuations of the position of amino acid residues in beta isoform than in alpha. In alpha isoform, the loop forms a salt bridge with loop switch 1, the bridge is not present in beta isoform. Two isoleucines I303 and I313 of loop S291-E317 are replaced with valines in alpha isoform. We introduced a double mutation I303V:I313V in beta isoform background and studied how the mutation affects the rate of ATP binding and ADP dissociation from actomyosin. We found that ATP-induced actomyosin dissociation occurs faster in the mutant, but the rate of ADP release remains the same as in the wild-type beta isoform. Due to the proximity of loop S291-E317 and loop switch 1, a faster rate of ATP-induced actomyosin dissociation indicates that loop S291-E317 affects structural dynamics of loop switch 1, and that loop switch 1 controls ATP binding to the active site. A similar rate of ADP dissociation from actomyosin in the mutant and wild-type myosin constructs indicates that loop switch 1 does not control ADP release from actomyosin.

Mavacamten stabilizes an autoinhibited state of two-headed cardiac myosin.

We used transient biochemical and structural kinetics to elucidate the molecular mechanism of mavacamten, an allosteric cardiac myosin inhibitor and a prospective treatment for hypertrophic cardiomyopathy. We find that mavacamten stabilizes an autoinhibited state of two-headed cardiac myosin not found in the single-headed S1 myosin motor fragment. We determined this by measuring cardiac myosin actin-activated and actin-independent ATPase and single-ATP turnover kinetics. A two-headed myosin fragment exhibits distinct autoinhibited ATP turnover kinetics compared with a single-headed fragment. Mavacamten enhanced this autoinhibition. It also enhanced autoinhibition of ADP release. Furthermore, actin changes the structure of the autoinhibited state by forcing myosin lever-arm rotation. Mavacamten slows this rotation in two-headed myosin but does not prevent it. We conclude that cardiac myosin is regulated in solution by an interaction between its two heads and propose that mavacamten stabilizes this state.

Heart failure drug changes the mechanoenzymology of the cardiac myosin powerstroke.

Omecamtiv mecarbil (OM), a putative heart failure therapeutic, increases cardiac contractility. We hypothesize that it does this by changing the structural kinetics of the myosin powerstroke. We tested this directly by performing transient time-resolved FRET on a ventricular cardiac myosin biosensor. Our results demonstrate that OM stabilizes myosin's prepowerstroke structural state, supporting previous measurements showing that the drug shifts the equilibrium constant for myosin-catalyzed ATP hydrolysis toward the posthydrolysis biochemical state. OM slowed the actin-induced powerstroke, despite a twofold increase in the rate constant for actin-activated phosphate release, the biochemical step in myosin's ATPase cycle associated with force generation and the conversion of chemical energy into mechanical work. We conclude that OM alters the energetics of cardiac myosin's mechanical cycle, causing the powerstroke to occur after myosin weakly binds to actin and releases phosphate. We discuss the physiological implications for these changes.

Dilated cardiomyopathy myosin mutants have reduced force-generating capacity.

Dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM) can cause arrhythmias, heart failure, and cardiac death. Here, we functionally characterized the motor domains of five DCM-causing mutations in human ß-cardiac myosin. Kinetic analyses of the individual events in the ATPase cycle revealed that each mutation alters different steps in this cycle. For example, different mutations gave enhanced or reduced rate constants of ATP binding, ATP hydrolysis, or ADP release or exhibited altered ATP, ADP, or actin affinity. Local effects dominated, no common pattern accounted for the similar mutant phenotype, and there was no distinct set of changes that distinguished DCM mutations from previously analyzed HCM myosin mutations. That said, using our data to model the complete ATPase contraction cycle revealed additional critical insights. Four of the DCM mutations lowered the duty ratio (the ATPase cycle portion when myosin strongly binds actin) because of reduced occupancy of the force-holding A·M·D complex in the steady state. Under load, the A·M·D state is predicted to increase owing to a reduced rate constant for ADP release, and this effect was blunted for all five DCM mutations. We observed the opposite effects for two HCM mutations, namely R403Q and R453C. Moreover, the analysis predicted more economical use of ATP by the DCM mutants than by WT and the HCM mutants. Our findings indicate that DCM mutants have a deficit in force generation and force-holding capacity due to the reduced occupancy of the force-holding state.

Multidimensional structure-function relationships in human ß-cardiac myosin from population-scale genetic variation.

Myosin motors are the fundamental force-generating elements of muscle contraction. Variation in the human ß-cardiac myosin heavy chain gene (MYH7) can lead to hypertrophic cardiomyopathy (HCM), a heritable disease characterized by cardiac hypertrophy, heart failure, and sudden cardiac death. How specific myosin variants alter motor function or clinical expression of disease remains incompletely understood. Here, we combine structural models of myosin from multiple stages of its chemomechanical cycle, exome sequencing data from two population cohorts of 60,706 and 42,930 individuals, and genetic and phenotypic data from 2,913 patients with HCM to identify regions of disease enrichment within ß-cardiac myosin. We first developed computational models of the human ß-cardiac myosin protein before and after the myosin power stroke. Then, using a spatial scan statistic modified to analyze genetic variation in protein 3D space, we found significant enrichment of disease-associated variants in the converter, a kinetic domain that transduces force from the catalytic domain to the lever arm to accomplish the power stroke. Focusing our analysis on surface-exposed residues, we identified a larger region significantly enriched for disease-associated variants that contains both the converter domain and residues on a single flat surface on the myosin head described as the myosin mesa. Notably, patients with HCM with variants in the enriched regions have earlier disease onset than patients who have HCM with variants elsewhere. Our study provides a model for integrating protein structure, large-scale genetic sequencing, and detailed phenotypic data to reveal insight into time-shifted protein structures and genetic disease.