Myosin-Catalyzed ATP Hydrolysis in the Presence of Disease-Causing Mutations: Mavacamten as a Way to Repair Mechanism.

Hypertrophic cardiomyopathy is one of the most common forms of genetic cardiomyopathy. Mavacamten is a first-in-class myosin modulator that was identified via activity screening on the wild type, and it is FDA-approved for the treatment of obstructive hypertrophic cardiomyopathy (HCM). The drug selectively binds to the cardiac ß-myosin, inhibiting myosin function to decrease cardiac contractility. Though the drug is thought to affect multiple steps of the myosin cross-bridge cycle, its detailed mechanism of action is still under investigation. Individual steps in the overall cross-bridge cycle must be queried to elucidate the full mechanism of action. In this study, we utilize the rare-event method of transition path sampling to generate reactive trajectories to gain insights into the action of the drug on the dynamics and rate of the ATP hydrolysis step for human cardiac ß-myosin. We study three known HCM causative myosin mutations: R453C, P710R, and R712L to observe the effect of the drug on the alterations caused by these mutations in the chemical step. Since the crystal structure of the drug-bound myosin was not available at the time of this work, we created a model of the drug-bound system utilizing a molecular docking approach. We find a significant effect of the drug in one case, where the actual mechanism of the reaction is altered from the wild type by mutation. The drug restores both the rate of hydrolysis to the wildtype level and the mechanism of the reaction. This is a way to check the effect of the drug on untested mutations.

The HCM - Linked Mutation Arg92Leu in TNNT2 Allosterically Alters the cTnC - cTnI Interface and Disrupts the PKA-mediated Regulation of Myofilament Relaxation.

Background: Impaired left ventricular relaxation, high filling pressures, and dysregulation of Ca 2+ homeostasis are common findings contributing to diastolic dysfunction in hypertrophic cardiomyopathy (HCM). Studies have shown that impaired relaxation is an early observation in the sarcomere-gene-positive preclinical HCM cohort which suggests potential involvement of myofilament regulators of relaxation. Yet, a molecular level understanding of mechanism(s) at the level of the myofilament is lacking. We hypothesized that mutation-specific, allosterically mediated, changes to the cardiac troponin C-cardiac troponin I (cTnC-cTnI) interface can account for the development of early-onset diastolic dysfunction via decreased PKA accessibility to cTnI. Methods: HCM mutations R92L-cTnT (Arg92Leu) and Δ160E-cTnT (Glu160 deletion) were studied in vivo , in vitro, and in silico via 2D echocardiography, western blotting, ex vivo hemodynamics, stopped-flow kinetics, time resolved fluorescence resonance energy transfer (TR-FRET), and molecular dynamics simulations. Results: The HCM-causative mutations R92L-cTnT and Δ160E-cTnT result in different time-of-onset of diastolic dysfunction. R92L-cTnT demonstrated early-onset diastolic dysfunction accompanied by a localized decrease in phosphorylation of cTnI. Constitutive phosphorylation of cTnI (cTnI-D 23 D 24 ) was sufficient to recover diastolic function to Non-Tg levels only for R92L-cTnT. Mutation-specific changes in Ca 2+ dissociation rates associated with R92L-cTnT reconstituted with cTnI-D 23 D 24 led us to investigate potential involvement of structural changes in the cTnC-cTnI interface as an explanation for these observations. We probed the interface via TR-FRET revealing a repositioning of the N-terminus of cTnI, closer to cTnC, and concomitant decreases in distance distributions at sites flanking the PKA consensus sequence. Implementing TR-FRET distances as constraints into our atomistic model identified additional electrostatic interactions at the consensus sequence. Conclusion: These data indicate that the early diastolic dysfunction observed in a subset of HCM is likely attributable to structural changes at the cTnC-cTnI interface that impair accessibility of PKA thereby blunting ß-adrenergic responsiveness and identifying a potential molecular target for therapeutic intervention.

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.

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 HPO42-, 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.

Thin filament cardiomyopathies: A review of genetics, disease mechanisms, and emerging therapeutics.

All muscle contraction occurs due to the cyclical interaction between sarcomeric thin and thick filament proteins within the myocyte. The thin filament consists of the proteins actin, tropomyosin, Troponin C, Troponin I, and Troponin T. Mutations in these proteins can result in various forms of cardiomyopathy, including hypertrophic, restrictive, and dilated phenotypes and account for as many as 30% of all cases of inherited cardiomyopathy. There is significant evidence that thin filament mutations contribute to dysregulation of Ca2+ within the sarcomere and may have a distinct pathomechanism of disease from cardiomyopathy associated with thick filament mutations. A number of distinct clinical findings appear to be correlated with thin-filament mutations: greater degrees of restrictive cardiomyopathy and relatively less left ventricular (LV) hypertrophy and LV outflow tract obstruction than that seen with thick filament mutations, increased morbidity associated with heart failure, increased arrhythmia burden and potentially higher mortality. Most therapies that improve outcomes in heart failure blunt the neurohormonal pathways involved in cardiac remodeling, while most therapies for hypertrophic cardiomyopathy involve use of negative inotropes to reduce LV hypertrophy or septal reduction therapies to reduce LV outflow tract obstruction. None of these therapies directly address the underlying sarcomeric dysfunction associated with thin-filament mutations. With mounting evidence that thin filament cardiomyopathies occur through a distinct mechanism, there is need for therapies targeting the unique, underlying mechanisms tailored for each patient depending on a given mutation.

Structure and Dynamics of the Flexible Cardiac Troponin T Linker Domain in a Fully Reconstituted Thin Filament.

The structural analysis of large protein complexes has been greatly enhanced through the application of electron microscopy techniques. One such multiprotein complex, the cardiac thin filament (cTF), has cyclic interactions with thick filament proteins to drive contraction of the heart that has recently been the subject of such studies. As important as these studies are, they provide limited or no information on highly flexible regions that in isolation would be characterized as inherently disordered. One such region is the extended cardiac troponin T (cTnT) linker between the regions of cTnT which have been labeled TNT1 and TNT2. It comprises a hinge region (residues 158-166) and a highly flexible region (residues 167-203). Critically, this region modulates the troponin/tropomyosin complex's position across the actin filament. Thus, the cTnT linker structure and dynamics are central to the regulation of the function of cardiac muscles, but up to now, it was ill-understood. To establish the cTnT linker structure, we coupled an atomistic computational cTF model with time-resolved fluorescence resonance energy transfer measurements in both ±Ca2+ conditions utilizing fully reconstituted cTFs. We mapped the cTnT linker's positioning across the actin filament, and by coupling the experimental results to computation, we found mean structures and ranges of motion of this part of the complex. With this new insight, we can now address cTnT linker structural dynamics in both myofilament activation and disease.

Free-Energy Surfaces of Two Cardiac Thin Filament Conformational Changes during Muscle Contraction.

The troponin core is an important regulatory complex in cardiac sarcomeres. Contraction is initiated by a calcium ion binding to cardiac troponin C (cTnC), initiating a conformational shift within the protein, altering its interactions with cardiac troponin I (cTnI). The change in cTnC-cTnI interactions prompts the C-terminal domain of cTnI to dissociate from actin, allowing tropomyosin to reveal myosin-binding sites on actin. Each of the concerted movements in the cardiac thin filament (CTF) is crucial for allowing the contraction of cardiomyocytes, yet little is known about the free energy associated with each transition, which is vital for understanding contraction on a molecular level. Using metadynamics, we calculated the free-energy surface of two transitions in the CTF: cTnC opening in the presence and absence of Ca2+ and cTnI dissociating from actin with both open and closed cTnC. These results not only provide the free-energy surface of the transitions but will also be shown to determine if the order of transitions in the contraction cycle is important. From our calculations, we found that the calcium ion helps stabilize the open conformation of cTnC and that the C-terminus of cTnI is stabilized by cTnC in the open conformation when dissociating from the actin surface.

Sex-Related Differences in Genetic Cardiomyopathies.

Cardiomyopathies are a heterogeneous collection of diseases that have in common primary functional and structural abnormalities of the heart muscle, often genetically determined. The most effective categorization of cardiomyopathies is based on the presenting phenotype, with hypertrophic, dilated, arrhythmogenic, and restrictive cardiomyopathy as the prototypes. Sex modulates the prevalence, morpho-functional manifestations and clinical course of cardiomyopathies. Aspects as diverse as ion channel expression and left ventricular remodeling differ in male and female patients with myocardial disease, although the reasons for this are poorly understood. Moreover, clinical differences may also result from complex societal/environmental discrepancies between sexes that may disadvantage women. This review provides a state-of-the-art appraisal of the influence of sex on cardiomyopathies, highlighting the many gaps in knowledge and open research questions.

Myosin modulators: emerging approaches for the treatment of cardiomyopathies and heart failure.

Myosin modulators are a novel class of pharmaceutical agents that are being developed to treat patients with a range of cardiomyopathies. The therapeutic goal of these drugs is to target cardiac myosins directly to modulate contractility and cardiac power output to alleviate symptoms that lead to heart failure and arrhythmias, without altering calcium signaling. In this Review, we discuss two classes of drugs that have been developed to either activate (omecamtiv mecarbil) or inhibit (mavacamten) cardiac contractility by binding to ß-cardiac myosin (MYH7). We discuss progress in understanding the mechanisms by which the drugs alter myosin mechanochemistry, and we provide an appraisal of the results from clinical trials of these drugs, with consideration for the importance of disease heterogeneity and genetic etiology for predicting treatment benefit.

Computational and biophysical determination of pathogenicity of variants of unknown significance in cardiac thin filament.

Point mutations within sarcomeric proteins have been associated with altered function and cardiomyopathy development. Difficulties remain, however, in establishing the pathogenic potential of individual mutations, often limiting the use of genotype in management of affected families. To directly address this challenge, we utilized our all-atom computational model of the human full cardiac thin filament (CTF) to predict how sequence substitutions in CTF proteins might affect structure and dynamics on an atomistic level. Utilizing molecular dynamics calculations, we simulated 21 well-defined genetic pathogenic cardiac troponin T and tropomyosin variants to establish a baseline of pathogenic changes induced in computational observables. Computational results were verified via differential scanning calorimetry on a subset of variants to develop an experimental correlation. Calculations were performed on 9 independent variants of unknown significance (VUS), and results were compared with pathogenic variants to identify high-resolution pathogenic signatures. Results for VUS were compared with the baseline set to determine induced structural and dynamic changes, and potential variant reclassifications were proposed. This unbiased, high-resolution computational methodology can provide unique structural and dynamic information that can be incorporated into existing analyses to facilitate classification both for de novo variants and those where established approaches have provided conflicting information.

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 HPO42- 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.

A Proposed Mechanism for the Initial Myosin Binding Event on the Cardiac Thin Filament: A Metadynamics Study.

The movement of tropomyosin over filamentous actin regulates the cross-bridge cycle of the thick with thin filament of cardiac muscle by blocking and revealing myosin binding sites. Tropomyosin exists in three, distinct equilibrium states with one state blocking myosin-actin interactions (blocked position) and the remaining two allowing for weak (closed position) and strong myosin binding (open position). However, experimental information illuminating how myosin binds to the thin filament and influences tropomyosin's transition across the actin surface is lacking. Using metadynamics, we mimic the effect of a single myosin head binding by determining the work required to pull small segments of tropomyosin toward the open position in several distinct regions of the thin filament. We find differences in required work due to the influence of cardiac troponin T lead to preferential binding sites and determine the mechanism of further myosin head recruitment.

Complexity in genetic cardiomyopathies and new approaches for mechanism-based precision medicine.

Genetic cardiomyopathies have been studied for decades, and it has become increasingly clear that these progressive diseases are more complex than originally thought. These complexities can be seen both in the molecular etiologies of these disorders and in the clinical phenotypes observed in patients. While these disorders can be caused by mutations in cardiac genes, including ones encoding sarcomeric proteins, the disease presentation varies depending on the patient mutation, where mutations even within the same gene can cause divergent phenotypes. Moreover, it is challenging to connect the mutation-induced molecular insult that drives the disease pathogenesis with the various compensatory and maladaptive pathways that are activated during the course of the subsequent progressive, pathogenic cardiac remodeling. These inherent complexities have frustrated our ability to understand and develop broadly effective treatments for these disorders. It has been proposed that it might be possible to improve patient outcomes by adopting a precision medicine approach. Here, we lay out a practical framework for such an approach, where patient subpopulations are binned based on common underlying biophysical mechanisms that drive the molecular disease pathogenesis, and we propose that this function-based approach will enable the development of targeted therapeutics that ameliorate these effects. We highlight several mutations to illustrate the need for mechanistic molecular experiments that span organizational and temporal scales, and we describe recent advances in the development of novel therapeutics based on functional targets. Finally, we describe many of the outstanding questions for the field and how fundamental mechanistic studies, informed by our more nuanced understanding of the clinical disorders, will play a central role in realizing the potential of precision medicine for genetic cardiomyopathies.

Mechanochemical Function of Myosin II: Investigation into the Recovery Stroke and ATP Hydrolysis.

Myosin regulates muscle function through a complex cycle of conformational rearrangements coupled with the hydrolysis of adenosine triphosphate (ATP). The recovery stroke reorganizes the myosin active site to hydrolyze ATP and cross bridge with the thin filament to produce muscle contraction. Engineered mutations K84M and R704E in Dictyostelium myosin have been designed to specifically inhibit the recovery stroke and have been shown to indirectly affect the ATPase activity of myosin. We investigated these mutagenic perturbations to the recovery stroke and generated thermodynamically correct and unbiased trajectories for native ATP hydrolysis with computationally enhanced sampling methods. Our methodology was able to resolve experimentally observed changes to kinetic and equilibrium dynamics for the recovery stroke with the correct prediction in the severity of these changes. For ATP hydrolysis, the sequential nature along with the stabilization of a metaphosphate intermediate was observed in agreement with previous studies. However, we observed glutamate 459 being utilized as a proton abstractor to prime the attacking water instead of a lytic water, a phenomenon not well categorized in myosin but has in other ATPases. Both rare event methodologies can be extended to human myosin to investigate isoformic differences from Dictyostelium and scan cardiomyopathic mutations to see differential perturbations to kinetics of other conformational changes in myosin such as the power stroke.

Modulating the tension-time integral of the cardiac twitch prevents dilated cardiomyopathy in murine hearts.

Dilated cardiomyopathy (DCM) is often associated with sarcomere protein mutations that confer reduced myofilament tension-generating capacity. We demonstrated that cardiac twitch tension-time integrals can be targeted and tuned to prevent DCM remodeling in hearts with contractile dysfunction. We employed a transgenic murine model of DCM caused by the D230N-tropomyosin (Tm) mutation and designed a sarcomere-based intervention specifically targeting the twitch tension-time integral of D230N-Tm hearts using multiscale computational models of intramolecular and intermolecular interactions in the thin filament and cell-level contractile simulations. Our models predicted that increasing the calcium sensitivity of thin filament activation using the cardiac troponin C (cTnC) variant L48Q can sufficiently augment twitch tension-time integrals of D230N-Tm hearts. Indeed, cardiac muscle isolated from double-transgenic hearts expressing D230N-Tm and L48Q cTnC had increased calcium sensitivity of tension development and increased twitch tension-time integrals compared with preparations from hearts with D230N-Tm alone. Longitudinal echocardiographic measurements revealed that DTG hearts retained normal cardiac morphology and function, whereas D230N-Tm hearts developed progressive DCM. We present a computational and experimental framework for targeting molecular mechanisms governing the twitch tension of cardiomyopathic hearts to counteract putative mechanical drivers of adverse remodeling and open possibilities for tension-based treatments of genetic cardiomyopathies.

Modifications of Sarcoplasmic Reticulum Function Prevent Progression of Sarcomere-Linked Hypertrophic Cardiomyopathy Despite a Persistent Increase in Myofilament Calcium Response.

Hypertrophic cardiomyopathy (HCM) is a genetic disorder caused by mutations in different genes mainly encoding myofilament proteins and therefore called a "disease of the sarcomere." Despite the discovery of sarcomere protein mutations linked to HCM almost 30 years ago, the cellular mechanisms responsible for the development of this disease are not completely understood and likely vary among different mutations. Moreover, despite many efforts to develop effective treatments for HCM, these have largely been unsuccessful, and more studies are needed to better understand the cellular mechanisms of the disease. In experiments reported here, we investigated a mouse model expressing the mutant cTnT-R92Q, which is linked to HCM and induces an increase in myofilament Ca2+ sensitivity and diastolic dysfunction. We found that early correction of the diastolic dysfunction by phospholamban knockout (PLNKO) was able to prevent the development of the HCM phenotype in troponin T (TnT)-R92Q transgenic (TG) mice. Four groups of mice in FVB/N background were generated and used for the experiments: (1) non-transgenic (NTG)/PLN mice, which express wild-type TnT and normal level of PLN; (2) NTG/PLNKO mice, which express wild-type TnT and no PLN; (3) TG/PLN mice, which express TnT-R92Q and normal level of PLN; (4) TG/PLNKO mice, which express TnT-R92Q and no PLN. Cardiac function was determined using both standard echocardiographic parameters and speckle tracking strain measurements. We found that both atrial morphology and diastolic function were altered in TG/PLN mice but normal in TG/PLNKO mice. Histological analysis showed a disarray of myocytes and increased collagen deposition only in TG/PLN hearts. We also observed increased Ca2+/calmodulin-dependent protein kinase II (CaMKII) phosphorylation only in TG/PLN hearts but not in TG/PLNKO hearts. The rescue of the HCM phenotype was not associated with differences in myofilament Ca2+ sensitivity between TG/PLN and TG/PLNKO mice. Moreover, compared to standard systolic echo parameters, such as ejection fraction (EF), speckle strain measurements provided a more sensitive approach to detect early systolic dysfunction in TG/PLN mice. In summary, our results indicate that targeting diastolic dysfunction through altering Ca2+ fluxes with no change in myofilament response to Ca2+ was able to prevent the development of the HCM phenotype and should be considered as a potential additional treatment for HCM patients.

FRET-based analysis of the cardiac troponin T linker region reveals the structural basis of the hypertrophic cardiomyopathy-causing Δ160E mutation.

Mutations in the cardiac thin filament (TF) have highly variable effects on the regulatory function of the cardiac sarcomere. Understanding the molecular-level dysfunction elicited by TF mutations is crucial to elucidate cardiac disease mechanisms. The hypertrophic cardiomyopathy-causing cardiac troponin T (cTnT) mutation Δ160Glu (Δ160E) is located in a putative "hinge" adjacent to an unstructured linker connecting domains TNT1 and TNT2. Currently, no high-resolution structure exists for this region, limiting significantly our ability to understand its role in myofilament activation and the molecular mechanism of mutation-induced dysfunction. Previous regulated in vitro motility data have indicated mutation-induced impairment of weak actomyosin interactions. We hypothesized that cTnT-Δ160E repositions the flexible linker, altering weak actomyosin electrostatic binding and acting as a biophysical trigger for impaired contractility and the observed remodeling. Using time-resolved FRET and an all-atom TF model, here we first defined the WT structure of the cTnT-linker region and then identified Δ160E mutation-induced positional changes. Our results suggest that the WT linker runs alongside the C terminus of tropomyosin. The Δ160E-induced structural changes moved the linker closer to the tropomyosin C terminus, an effect that was more pronounced in the presence of myosin subfragment (S1) heads, supporting previous findings. Our in silico model fully supported this result, indicating a mutation-induced decrease in linker flexibility. Our findings provide a framework for understanding basic pathogenic mechanisms that drive severe clinical hypertrophic cardiomyopathy phenotypes and for identifying structural targets for intervention that can be tested in silico and in vitro.