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

Myosin's powerstroke transitions define atomic scale movement of cardiac thin filament tropomyosin

Dynamic interactions between the myosin motor head on thick filaments and the actin molecular track on thin filaments drive the myosin-crossbridge cycle that powers muscle contraction. The process is initiated by Ca2+ and the opening of troponin-tropomyosin-blocked myosin-binding sites on actin. The ensuing recruitment of myosin heads and their transformation from pre-powerstroke to post-powerstroke conformation on actin produce the force required for contraction. Cryo-EM-based atomic models confirm that during this process, tropomyosin occupies three different average positions on actin. Tropomyosin pivoting on actin away from a TnI-imposed myosin-blocking position accounts for part of the Ca2+ activation observed. However, the structure of tropomyosin on thin filaments that follows pre-powerstroke myosin binding and its translocation during myosin's pre-powerstroke to post-powerstroke transition remains unresolved. Here, we approach this transition computationally in silico. We used the myosin helix-loop-helix motif as an anchor to dock models of pre-powerstroke cardiac myosin to the cleft between neighboring actin subunits along cardiac thin filaments. We then performed targeted molecular dynamics simulations of the transition between pre- and post-powerstroke conformations on actin in the presence of cardiac troponin-tropomyosin. These simulations show Arg 369 and Glu 370 on the tip of myosin Loop-4 encountering identically charged residues on tropomyosin. The charge repulsion between residues causes tropomyosin translocation across actin, thus accounting for the final regulatory step in the activation of the thin filament, and, in turn, facilitating myosin movement along the filament. We suggest that during muscle activity, myosin-induced tropomyosin movement is likely to result in unencumbered myosin head interactions on actin at low-energy cost.

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.

Pharmacology of LRRK2 with type I and II kinase inhibitors revealed by cryo-EM

LRRK2 is one of the most promising drug targets for Parkinson's disease. Though type I kinase inhibitors of LRRK2 are under clinical trials, alternative strategies like type II inhibitors are being actively pursued due to the potential undesired effects of type I inhibitors. Currently, a robust method for LRRK2-inhibitor structure determination to guide structure-based drug discovery is lacking, and inhibition mechanisms of available compounds are also unclear. Here we present near-atomic-resolution structures of LRRK2 with type I (LRRK2-IN-1 and GNE-7915) and type II (rebastinib, ponatinib, and GZD-824) inhibitors, uncovering the structural basis of LRRK2 inhibition and conformational plasticity of the kinase domain with molecular dynamics (MD) simulations. Type I and II inhibitors bind to LRRK2 in active-like and inactive conformations, so LRRK2-inhibitor complexes further reveal general structural features associated with LRRK2 activation. Our study provides atomic details of LRRK2-inhibitor interactions and a framework for understanding LRRK2 activation and for rational drug design.

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 ( k cat ) and decreases (6-fold) the actin concentration at which ATPase is one-half maximal ( K ATPase ). 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.

Significance Statement

Heart disease can be caused by inherited mutations in beta-cardiac myosin, the molecular motor that powers systolic contraction in the ventricles of the heart. However, it remains unclear how these mutations lead to contractile dysfunction and pathogenic remodeling of the heart. We investigated a unique dilated cardiomyopathy mutation (E525K) that dramatically stabilizes the autoinhibited state while enhancing intrinsic motor function. Thus, we examined how this mutation impacts transient kinetic steps of the ATPase cycle, motile properties, and structural changes associated with the power stroke and phosphate release. Our results provide a kinetic and structural basis for how beta-cardiac myosin mutations may disrupt molecular-level contractile function in complex ways, which may inform the development of targeted therapeutics.