Mini‐ review: Recent insights into the relative timing of myosin’s powerstroke and release of phosphate

A novel kinetic model to demonstrate the independent effects of ATP and ADP/Pi concentrations on sarcomere function

Understanding muscle contraction mechanisms is a standing challenge, and one of the approaches has been to create models of the sarcomere-the basic contractile unit of striated muscle. While these models have been successful in elucidating many aspects of muscle contraction, they fall short in explaining the energetics of functional phenomena, such as rigor, and in particular, their dependence on the concentrations of the biomolecules involved in the cross-bridge cycle. Our hypothesis posits that the stochastic time delay between ATP adsorption and ADP/Pi release in the cross-bridge cycle necessitates a modeling approach where the rates of these two reaction steps are controlled by two independent parts of the total free energy change of the hydrolysis reaction. To test this hypothesis, we built a two-filament, stochastic-mechanical half-sarcomere model that separates the energetic roles of ATP and ADP/Pi in the cross-bridge cycle's free energy landscape. Our results clearly demonstrate that there is a nontrivial dependence of the cross-bridge cycle's kinetics on the independent concentrations of ATP, ADP, and Pi. The simplicity of the proposed model allows for analytical solutions of the more basic systems, which provide novel insight into the dominant mechanisms driving some of the experimentally observed contractile phenomena.

Human cardiac β-myosin powerstroke energetics: Thin filament, Pi displacement, and mutation effects

The powerstroke of human cardiac β-myosin is an important stage of the cross-bridge cycle that generates force for muscle contraction. However, the starting structure of this process has never been resolved, and the relative timing of the powerstroke and inorganic phosphate (Pi) release is still controversial. In this study, we generated an atomistic model of myosin on the thin filament and utilized metadynamics simulations to predict the absent starting structure of the powerstroke. We demonstrated that the displacement of Pi from the active site during the powerstroke is likely necessary, reducing the energy barrier of the conformation change. The effects of the presence of the thin filament, the hypertrophic cardiomyopathy mutation R712L, and the binding of mavacamten on the powerstroke process were also investigated.

Study of the Myosin Relay Helix Peptide by Molecular Dynamics Simulations, Pump-Probe and 2D Infrared Spectroscopy

The Davydov model was conjectured to describe how an amide I excitation created during ATP hydrolysis in myosin might be significant in providing energy to drive myosin's chemomechanical cycle. The free energy surfaces of the myosin relay helix peptide dissolved in 2,2,2-trifluoroethanol (TFE), determined by metadynamics simulations, demonstrate local minima differing in free energy by only ~2 kT, corresponding to broken and stabilized hydrogen bonds, respectively. Experimental pump-probe and 2D infrared spectroscopy were performed on the peptide dissolved in TFE. The relative heights of two peaks seen in the pump-probe data and the corresponding relative volumes of diagonal peaks seen in the 2D-IR spectra at time delays between 0.5 ps and 1 ps differ noticeably from what is seen at earlier or later time delays or in the linear spectrum, indicating that a vibrational excitation may influence the conformational state of this helix. Thus, it is possible that the presence of an amide I excitation may be a direct factor in the conformational state taken on by the myosin relay helix following ATP hydrolysis in myosin.

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.

Engines of change: Nonmuscle myosin II in mechanobiology

The emergence of mechanobiology has unveiled complex mechanisms by which cells adjust intracellular force production to their needs. Most communicable intracellular forces are generated by myosin II, an actin-associated molecular motor that transforms adenosine triphosphate (ATP) hydrolysis into contraction in nonmuscle and muscle cells. Myosin II-dependent force generation is tightly regulated, and deregulation is associated with specific pathologies. Here, we focus on the role of myosin II (nonmuscle myosin II, NMII) in force generation and mechanobiology. We outline the regulation and molecular mechanism of force generation by NMII, focusing on the actual outcome of contraction, that is, force application to trigger mechanosensitive events or the building of dissipative structures. We describe how myosin II-generated forces drive two major types of events: modification of the cellular morphology and/or triggering of genetic programs, which enhance the ability of cells to adapt to, or modify, their microenvironment. Finally, we address whether targeting myosin II to impair or potentiate its activity at the motor level is a viable therapeutic strategy, as illustrated by recent examples aimed at modulating cardiac myosin II function in heart disease.

Modeling thick filament activation suggests a molecular basis for force depression

Multiscale models aiming to connect muscle's molecular and cellular function have been difficult to develop, in part due to a lack of self-consistent multiscale data. To address this gap, we measured the force response from single, skinned rabbit psoas muscle fibers to ramp shortenings and step stretches performed on the plateau region of the force-length relationship. We isolated myosin from the same muscles and, under similar conditions, performed single-molecule and ensemble measurements of myosin's ATP-dependent interaction with actin using laser trapping and in vitro motility assays. We fit the fiber data by developing a partial differential equation model that includes thick filament activation, whereby an increase in force on the thick filament pulls myosin out of an inhibited state. The model also includes a series elastic element and a parallel elastic element. This parallel elastic element models a titin-actin interaction proposed to account for the increase in isometric force after stretch (residual force enhancement). By optimizing the model fit to a subset of our fiber measurements, we specified seven unknown parameters. The model then successfully predicted the remainder of our fiber measurements and also our molecular measurements from the laser trap and in vitro motility. The success of the model suggests that our multiscale data are self-consistent and can serve as a testbed for other multiscale models. Moreover, the model captures the decrease in isometric force observed in our muscle fibers after active shortening (force depression), suggesting a molecular mechanism for force depression, whereby a parallel elastic element combines with thick filament activation to decrease the number of cycling cross-bridges.

Mechanosensitive accumulation of non-muscle myosin IIB during mitosis requires its translocation activity

Non-muscle myosin II (NMII) is a force-generating mechanosensitive enzyme that responds to mechanical forces. NMIIs mechanoaccumulate at the cell cortex in response to mechanical forces. It is essential for cells to mechanically adapt to the physical environment, failure of which results in mitotic defects when dividing in confined environment. Much less is known about how NMII mechanoaccumulation is regulated during mitosis. We show that mitotic cells respond to compressive stress by promoting accumulation of active RhoA at the cell cortex as in interphase cells. RhoA mechanoresponse during mitosis activates and stabilizes NMIIB via ROCK signaling, leading to NMIIB mechanoaccumulation at the cell cortex. Using disease-related myosin II mutations, we found that NMIIB mechanoaccumulation requires its motor activity that translocates actin filaments, but not just its actin-binding function. Thus, the motor activity coordinates structural movement and nucleotide state changes to fine-tune actin-binding affinity optimal for NMIIs to generate and respond to forces.

Modeling Thick Filament Activation Suggests a Molecular Basis for Force Depression

Multiscale models aiming to connect muscle's molecular and cellular function have been difficult to develop, in part, due to a lack of self-consistent multiscale data. To address this gap, we measured the force response from single skinned rabbit psoas muscle fibers to ramp shortenings and step stretches performed on the plateau region of the force-length relationship. We isolated myosin from the same muscles and, under similar conditions, performed single molecule and ensemble measurements of myosin's ATP-dependent interaction with actin using laser trapping and in vitro motility assays. We fit the fiber data by developing a partial differential equation model that includes thick filament activation, whereby an increase in force on the thick filament pulls myosin out of an inhibited state. The model also includes a series elastic element and a parallel elastic element. This parallel elastic element models a titin-actin interaction proposed to account for the increase in isometric force following stretch (residual force enhancement). By optimizing the model fit to a subset of our fiber measurements, we specified seven unknown parameters. The model then successfully predicted the remainder of our fiber measurements and also our molecular measurements from the laser trap and in vitro motility. The success of the model suggests that our multiscale data are self-consistent and can serve as a testbed for other multiscale models. Moreover, the model captures the decrease in isometric force observed in our muscle fibers after active shortening (force depression), suggesting a molecular mechanism for force depression, whereby a parallel elastic element combines with thick filament activation to decrease the number of cycling cross-bridges.

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.

Revolving hexameric ATPases as asymmetric motors to translocate double-stranded DNA genome along one strand

DsDNA translocation through nanoscale pores is generally accomplished by ATPase biomotors. The discovery of the revolving dsDNA translocation mechanism, as opposed to rotation, in bacteriophage phi29 elucidated how ATPase motors move dsDNA. Revolution-driven, hexameric dsDNA motors have been reported in herpesvirus, bacterial FtsK, Streptomyces TraB, and T7 phage. This review explores the common relationship between their structure and mechanisms. Commonalities include moving along the 5'→3' strand, inchworm sequential action leading to an asymmetrical structure, channel chirality, channel size, and 3-step channel gating for controlling motion direction. The revolving mechanism and contact with one of the dsDNA strands addresses the historic controversy of dsDNA packaging using nicked, gapped, hybrid, or chemically modified DNA. These controversies surrounding dsDNA packaging activity using modified materials can be answered by whether the modification was introduced into the 3'→5' or 5'→3' strand. Perspectives concerning solutions to the controversy of motor structure and stoichiometry are also discussed.

A mutation in switch I alters the load-dependent kinetics of myosin Va

Myosin Va is the molecular motor that drives intracellular vesicular transport, powered by the transduction of chemical energy from ATP into mechanical work. The coupling of the powerstroke and phosphate (Pi) release is key to understanding the transduction process, and crucial details of this process remain unclear. Therefore, we determined the effect of elevated Pi on the force-generating capacity of a mini-ensemble of myosin Va S1 (WT) in a laser trap assay. By increasing the stiffness of the laser trap we determined the effect of increasing resistive loads on the rate of Pi-induced detachment from actin, and quantified this effect using the Bell approximation. We observed that WT myosin generated higher forces and larger displacements at the higher laser trap stiffnesses in the presence of 30 mM Pi, but binding event lifetimes decreased dramatically, which is most consistent with the powerstroke preceding the release of Pi from the active site. Repeating these experiments using a construct with a mutation in switch I of the active site (S217A) caused a seven-fold increase in the load-dependence of the Pi-induced detachment rate, suggesting that the S217A region of switch I may help mediate the load-dependence of Pi-rebinding.

Opportunities and Challenges for the Development of MRCK Kinases Inhibitors as Potential Cancer Chemotherapeutics

Cytoskeleton organization and dynamics are rapidly regulated by post-translational modifications of key target proteins. Acting downstream of the Cdc42 GTPase, the myotonic dystrophy-related Cdc42-binding kinases MRCKα, MRCKβ, and MRCKγ have recently emerged as important players in cytoskeleton regulation through the phosphorylation of proteins such as the regulatory myosin light chain proteins. Compared with the closely related Rho-associated coiled-coil kinases 1 and 2 (ROCK1 and ROCK2), the contributions of the MRCK kinases are less well characterized, one reason for this being that the discovery of potent and selective MRCK pharmacological inhibitors occurred many years after the discovery of ROCK inhibitors. The disclosure of inhibitors, such as BDP5290 and BDP9066, that have marked selectivity for MRCK over ROCK, as well as the dual ROCK + MRCK inhibitor DJ4, has expanded the repertoire of chemical biology tools to study MRCK function in normal and pathological conditions. Recent research has used these novel inhibitors to establish the role of MRCK signalling in epithelial polarization, phagocytosis, cytoskeleton organization, cell motility, and cancer cell invasiveness. Furthermore, pharmacological MRCK inhibition has been shown to elicit therapeutically beneficial effects in cell-based and in vivo studies of glioma, skin, and ovarian cancers.

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.