Modulation of post-powerstroke dynamics in myosin II by 2'-deoxy-ADP

Interacting myosin head 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 inhibited from performing 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 can pathologically disrupt IHM structure and the mechanical performance of muscle tissue. Thus, the IHM state has become a target for therapeutic intervention. An ATP analog called 2'-deoxy-ATP (dATP) is a potent myosin activator that 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 of the IHM state containing ADP.Pi in both nucleotide binding pockets revealed dynamic motions of the blocked head-free head interface, light chain binding domain, and S2 in this "inactive" state of myosin. Replacement of ADP.Pi by dADP.Pi triggered a series of structural changes that increased heterogeneity among residue contact pairs at the blocked head-free head interface and a 14% decrease in the interaction energy at the interface. Dynamic changes to this interface were accompanied by dynamics in the light chain binding region. A comparative analysis of these dynamics predicted new structural sites that may affect IHM stability.

Multiscale modeling shows how 2'-deoxy-ATP rescues ventricular function in heart failure

2'-deoxy-ATP (dATP) improves cardiac function by increasing the rate of crossbridge cycling and Ca[Formula: see text] transient decay. However, the mechanisms of these effects and how therapeutic responses to dATP are achieved when dATP is only a small fraction of the total ATP pool remain poorly understood. Here, we used a multiscale computational modeling approach to analyze the mechanisms by which dATP improves ventricular function. We integrated atomistic simulations of prepowerstroke myosin and actomyosin association, filament-scale Markov state modeling of sarcomere mechanics, cell-scale analysis of myocyte Ca[Formula: see text] dynamics and contraction, organ-scale modeling of biventricular mechanoenergetics, and systems level modeling of circulatory dynamics. Molecular and Brownian dynamics simulations showed that dATP increases the actomyosin association rate by 1.9 fold. Markov state models predicted that dATP increases the pool of myosin heads available for crossbridge cycling, increasing steady-state force development at low dATP fractions by 1.3 fold due to mechanosensing and nearest-neighbor cooperativity. This was found to be the dominant mechanism by which small amounts of dATP can improve contractile function at myofilament to organ scales. Together with faster myocyte Ca[Formula: see text] handling, this led to improved ventricular contractility, especially in a failing heart model in which dATP increased ejection fraction by 16% and the energy efficiency of cardiac contraction by 1%. This work represents a complete multiscale model analysis of a small molecule myosin modulator from single molecule to organ system biophysics and elucidates how the molecular mechanisms of dATP may improve cardiovascular function in heart failure with reduced ejection fraction.

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.

Incomplete-penetrant hypertrophic cardiomyopathy MYH7 G256E mutation causes hypercontractility and elevated mitochondrial respiration

Determining the pathogenicity of hypertrophic cardiomyopathy-associated mutations in the β-myosin heavy chain (MYH7) can be challenging due to its variable penetrance and clinical severity. This study investigates the early pathogenic effects of the incomplete-penetrant MYH7 G256E mutation on myosin function that may trigger pathogenic adaptations and hypertrophy. We hypothesized that the G256E mutation would alter myosin biomechanical function, leading to changes in cellular functions. We developed a collaborative pipeline to characterize myosin function across protein, myofibril, cell, and tissue levels to determine the multiscale effects on structure-function of the contractile apparatus and its implications for gene regulation and metabolic state. The G256E mutation disrupts the transducer region of the S1 head and reduces the fraction of myosin in the folded-back state by 33%, resulting in more myosin heads available for contraction. Myofibrils from gene-edited MYH7WT/G256E human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) exhibited greater and faster tension development. This hypercontractile phenotype persisted in single-cell hiPSC-CMs and engineered heart tissues. We demonstrated consistent hypercontractile myosin function as a primary consequence of the MYH7 G256E mutation across scales, highlighting the pathogenicity of this gene variant. Single-cell transcriptomic and metabolic profiling demonstrated upregulated mitochondrial genes and increased mitochondrial respiration, indicating early bioenergetic alterations. This work highlights the benefit of our multiscale platform to systematically evaluate the pathogenicity of gene variants at the protein and contractile organelle level and their early consequences on cellular and tissue function. We believe this platform can help elucidate the genotype-phenotype relationships underlying other genetic cardiovascular diseases.

Atomistic Simulations of Sarcomere Proteins

Concerted atomic motions are requisite for sarcomere protein function and may become disrupted in HCM pathologies. Computational approaches such as molecular dynamics simulation can resolve such dynamics with unrivalled spatial and temporal resolution. This chapter describes methods to model structural and dynamical changes in biomolecules with HCM-associated perturbations.

Multi-scale models reveal hypertrophic cardiomyopathy MYH7 G256E mutation drives hypercontractility and elevated mitochondrial respiration

Rationale

Over 200 mutations in the sarcomeric protein β-myosin heavy chain (MYH7) have been linked to hypertrophic cardiomyopathy (HCM). However, different mutations in MYH7 lead to variable penetrance and clinical severity, and alter myosin function to varying degrees, making it difficult to determine genotype-phenotype relationships, especially when caused by rare gene variants such as the G256E mutation.

Objective

This study aims to determine the effects of low penetrant MYH7 G256E mutation on myosin function. We hypothesize that the G256E mutation would alter myosin function, precipitating compensatory responses in cellular functions.

Methods

We developed a collaborative pipeline to characterize myosin function at multiple scales (protein to myofibril to cell to tissue). We also used our previously published data on other mutations to compare the degree to which myosin function was altered.

Results

At the protein level, the G256E mutation disrupts the transducer region of the S1 head and reduces the fraction of myosin in the folded-back state by 50.9%, suggesting more myosins available for contraction. Myofibrils isolated from hiPSC-CMs CRISPR-edited with G256E (MYH7 WT/G256E ) generated greater tension, had faster tension development and slower early phase relaxation, suggesting altered myosin-actin crossbridge cycling kinetics. This hypercontractile phenotype persisted in single-cell hiPSC-CMs and engineered heart tissues. Single-cell transcriptomic and metabolic profiling demonstrated upregulation of mitochondrial genes and increased mitochondrial respiration, suggesting altered bioenergetics as an early feature of HCM.

Conclusions

MYH7 G256E mutation causes structural instability in the transducer region, leading to hypercontractility across scales, perhaps from increased myosin recruitment and altered crossbridge cycling. Hypercontractile function of the mutant myosin was accompanied by increased mitochondrial respiration, while cellular hypertrophy was modest in the physiological stiffness environment. We believe that this multi-scale platform will be useful to elucidate genotype-phenotype relationships underlying other genetic cardiovascular diseases.

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.

Structural OFF/ON transitions of myosin in relaxed porcine myocardium predict calcium-activated force

Contraction in striated muscle is initiated by calcium binding to troponin complexes, but it is now understood that dynamic transition of myosin between resting, ordered OFF states on thick filaments and active, disordered ON states that can bind to thin filaments is critical in regulating muscle contractility. These structural OFF to ON transitions of myosin are widely assumed to correspond to transitions from the biochemically defined, energy-sparing, super-relaxed (SRX) state to the higher ATPase disordered-relaxed (DRX) state. Here we examined the effect of 2'-deoxy-ATP (dATP), a naturally occurring energy substrate for myosin, on the structural OFF to ON transitions of myosin motors in porcine cardiac muscle thick filaments. Small-angle X-ray diffraction revealed that titrating dATP in relaxation solutions progressively moves the myosin heads from ordered OFF states on the thick filament backbone to disordered ON states closer to thin filaments. Importantly, we found that the structural OFF to ON transitions are not equivalent to the biochemically defined SRX to DRX transitions and that the dATP-induced structural OFF to ON transitions of myosin motors in relaxed muscle are strongly correlated with submaximal force augmentation by dATP. These results indicate that structural OFF to ON transitions of myosin in relaxed muscle can predict the level of force attained in calcium-activated cardiac muscle. Computational modeling and stiffness measurements suggest a final step in the OFF to ON transition may involve a subset of DRX myosins that form weakly bound cross-bridges prior to becoming active force-producing cross-bridges.