Pi release from myosin: a simulation analysis of possible pathways

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.

Homologous mutations in human β, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in β-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman-Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known whether their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human β, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins but minimal effects in β myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing overall enzymatic (ATPase) cycle rate. In contrast, the only measured effect of R671C in β myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not β, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are a testament to myosin's highly allosteric nature.

Molecular mechanisms of inorganic-phosphate release from the core and barbed end of actin filaments

The release of inorganic phosphate (Pi) from actin filaments constitutes a key step in their regulated turnover, which is fundamental to many cellular functions. The mechanisms underlying Pi release from the core and barbed end of actin filaments remain unclear. Here, using human and bovine actin isoforms, we combine cryo-EM with molecular-dynamics simulations and in vitro reconstitution to demonstrate how actin releases Pi through a 'molecular backdoor'. While constantly open at the barbed end, the backdoor is predominantly closed in filament-core subunits and opens only transiently through concerted amino acid rearrangements. This explains why Pi escapes rapidly from the filament end but slowly from internal subunits. In a nemaline-myopathy-associated actin variant, the backdoor is predominantly open in filament-core subunits, resulting in accelerated Pi release and filaments with drastically shortened ADP-Pi caps. Our results provide the molecular basis for Pi release from actin and exemplify how a disease-linked mutation distorts the nucleotide-state distribution and atomic structure of the filament.

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.

Multistep orthophosphate release tunes actomyosin energy transduction

Muscle contraction and a range of critical cellular functions rely on force-producing interactions between myosin motors and actin filaments, powered by turnover of adenosine triphosphate (ATP). The relationship between release of the ATP hydrolysis product ortophosphate (Pi) from the myosin active site and the force-generating structural change, the power-stroke, remains enigmatic despite its central role in energy transduction. Here, we present a model with multistep Pi-release that unifies current conflicting views while also revealing additional complexities of potential functional importance. The model is based on our evidence from kinetics, molecular modelling and single molecule fluorescence studies of Pi binding outside the active site. It is also consistent with high-speed atomic force microscopy movies of single myosin II molecules without Pi at the active site, showing consecutive snapshots of pre- and post-power stroke conformations. In addition to revealing critical features of energy transduction by actomyosin, the results suggest enzymatic mechanisms of potentially general relevance.

Release of the ATP hydrolysis product orthophosphate (Pi) from the myosin active site is central in force generation but is poorly understood. Here, Moretto et al. present evidence for multistep Pi-release reconciling apparently contradictory results.

Cryo-EM structure of the autoinhibited state of myosin-2

We solved the near-atomic resolution structure of smooth muscle myosin-2 in the autoinhibited state (10S) using single-particle cryo–electron microscopy. The 3.4-Å structure reveals the precise molecular architecture of 10S and the structural basis for myosin-2 regulation. We reveal the position of the phosphorylation sites that control myosin autoinhibition and activation by phosphorylation of the regulatory light chain. Further, we present a previously unidentified conformational state in myosin-2 that traps ADP and Pi produced by the hydrolysis of ATP in the active site. This noncanonical state represents a branch of the myosin enzyme cycle and explains the autoinhibition of the enzyme function of 10S along with its reduced affinity for actin. Together, our structure defines the molecular mechanisms that drive 10S formation, stabilization, and relief by phosphorylation of the regulatory light chain.

A reverse stroke characterizes the force generation of cardiac myofilaments, leading to an understanding of heart function

Changes in the molecular properties of cardiac myosin strongly affect the interactions of myosin with actin that result in cardiac contraction and relaxation. However, it remains unclear how myosin molecules work together in cardiac myofilaments and which properties of the individual myosin molecules impact force production to drive cardiac contractility. Here, we measured the force production of cardiac myofilaments using optical tweezers. The measurements revealed that stepwise force generation was associated with a higher frequency of backward steps at lower loads and higher stall forces than those of fast skeletal myofilaments. To understand these unique collective behaviors of cardiac myosin, the dynamic responses of single cardiac and fast skeletal myosin molecules, interacting with actin filaments, were evaluated under load. The cardiac myosin molecules switched among three distinct conformational positions, ranging from pre- to post-power stroke positions, in 1 mM ADP and 0 to 10 mM phosphate solution. In contrast to cardiac myosin, fast skeletal myosin stayed primarily in the post-power stroke position, suggesting that cardiac myosin executes the reverse stroke more frequently than fast skeletal myosin. To elucidate how the reverse stroke affects the force production of myofilaments and possibly heart function, a simulation model was developed that combines the results from the single-molecule and myofilament experiments. The results of this model suggest that the reversal of the cardiac myosin power stroke may be key to characterizing the force output of cardiac myosin ensembles and possibly to facilitating heart contractions.

Step-Wise Hydration of Magnesium by Four Water Molecules Precedes Phosphate Release in a Myosin Motor

Molecular motors, such as myosin, kinesin, and dynein, convert the energy released by the hydrolysis of ATP into mechanical work, thus allowing them to undergo directional motion on cytoskeletal tracks. A pivotal step in the chemomechanical transduction in myosin motors occurs after they bind to the actin filament, which triggers the release of phosphate (P_i, product of ATP hydrolysis) and the rotation of the lever arm. Here, we investigate the mechanism of phosphate release in myosin VI using extensive molecular dynamics simulations involving multiple trajectories of several μs. Because the escape of phosphate is expected to occur on time-scales on the order of milliseconds or more in myosin VI, we observed P_i release only if the trajectories were initiated with a rotated phosphate inside the nucleotide binding pocket. We discovered that although P_i populates the traditional "back door" route, phosphate exits through various other gateways, thus establishing the heterogeneity in the escape routes. Remarkably, we observed that the release of phosphate is preceded by a stepwise hydration of the ADP-bound magnesium ion. The release of the anion occurred only after four water molecules hydrated the cation (Mg²⁺). By performing comparative structural analyses, we show that hydration of magnesium is the key step in the phosphate release in a number of ATPases and GTPases. Nature may have evolved hydration of Mg²⁺ as a general molecular switch for P_i release, which is a universal step in the catalytic cycle of many machines that share little sequence or structural similarity.

Improved Inhibitory and Absorption, Distribution, Metabolism, Excretion, and Toxicology (ADMET) Properties of Blebbistatin Derivatives Indicate That Blebbistatin Scaffold Is Ideal for drug Development Targeting Myosin-2

Blebbistatin, para-nitroblebbistatin (NBleb), and para-aminoblebbistatin (AmBleb) are highly useful tool compounds as they selectively inhibit the ATPase activity of myosin-2 family proteins. Despite the medical importance of the myosin-2 family as drug targets, chemical optimization has not yet provided a promising lead for drug development because previous structure-activity-relationship studies were limited to a single myosin-2 isoform. Here we evaluated the potential of blebbistatin scaffold for drug development and found that D-ring substitutions can fine-tune isoform specificity, absorption-distribution-metabolism-excretion, and toxicological properties. We defined the inhibitory properties of NBleb and AmBleb on seven different myosin-2 isoforms, which revealed an unexpected potential for isoform specific inhibition. We also found that NBleb metabolizes six times slower than blebbistatin and AmBleb in rats, whereas AmBleb metabolizes two times slower than blebbistatin and NBleb in human, and that AmBleb accumulates in muscle tissues. Moreover, mutagenicity was also greatly reduced in case of AmBleb. These results demonstrate that small substitutions have beneficial functional and pharmacological consequences, which highlight the potential of the blebbistatin scaffold for drug development targeting myosin-2 family proteins and delineate a route for defining the chemical properties of further derivatives to be developed. SIGNIFICANCE STATEMENT: Small substitutions on the blebbistatin scaffold have beneficial functional and pharmacological consequences, highlighting their potential in drug development targeting myosin-2 family proteins.

Heart Failure Drug Modifies the Intrinsic Dynamics of the Pre-Power Stroke State of Cardiac Myosin

Omecamtiv mecarbil (OM), currently investigated for the treatment of heart failure, is the first example of a new class of drugs (cardiac myotropes) that can modify muscle contractility by directly targeting sarcomeric proteins. Using atomistic molecular dynamics simulations, we show that the binding of OM to the pre-power stroke state of cardiac myosin inhibits the functional motions of the protein and potentially affects Pi release from the nucleotide binding site. We also show that the changes in myosin ATPase activity induced by a set of OM analogues can be predicted from their relative affinity to the pre-power stroke state compared to the near rigor one, indicating that conformational selectivity plays an important role in determining the activity of these compounds.

A restrained locally enhanced sampling method (RLES) for finding free energy minima in complex systems

We present an extension of the locally enhanced sampling method. A restraint potential is introduced to drive the many-replica system to the canonical ensemble corresponding to the physical, single-replica system. Convergence properties are demonstrated using a model rugged two-dimensional potential, for which sampling by conventional equilibrium molecular dynamics is inefficient. Restrained locally enhanced sampling (RLES) is found to explore the space of configurations with an efficiency comparable to that of temperature replica exchange. To demonstrate the potential of RLES for realistic applications, the method is used to fold the 12-residue tryptophan zipper miniprotein in explicit solvent. The RLES algorithm can be incorporated into existing LES implementations with minor code modifications.

Structural basis for power stroke vs. Brownian ratchet mechanisms of motor proteins

Two mechanisms have been proposed for the function of motor proteins: The power stroke and the Brownian ratchet. The former refers to generation of a large downhill free energy gradient over which the motor protein moves nearly irreversibly in making a step, whereas the latter refers to biasing or rectifying the diffusive motion of the motor. Both mechanisms require input of free energy, which generally involves the processing of an ATP (adenosine 5'-triphosphate) molecule. Recent advances in experiments that reveal the details of the stepping motion of motor proteins, together with computer simulations of atomistic structures, have provided greater insights into the mechanisms. Here, we compare the various models of the power stroke and the Brownian ratchet that have been proposed. The 2 mechanisms are not mutually exclusive, and various motor proteins employ them to different extents to perform their biological function. As examples, we discuss linear motor proteins Kinesin-1 and myosin-V, and the rotary motor F₁-ATPase, all of which involve a power stroke as the essential element of their stepping mechanism.

Overview of the mechanism of cytoskeletal motors based on structure

In the last two decades, a wealth of structural and functional knowledge has been obtained for the three major cytoskeletal motor proteins, myosin, kinesin and dynein, which we review here. The cytoskeletal motor proteins myosin and kinesin are structurally similar in the core architecture of their motor domains and have similar force-producing mechanisms that are coupled with the chemical cycles of ATP binding, hydrolysis, Pi release and subsequent ADP release. The force is generated through conformational changes in the motor domain during Pi release and ATP binding in myosin and kinesin, respectively, and then converted into the rotation of the lever arm or neck linker (referred to as a power stroke) through the common structural pathways. On the other hand, the dynein cytoskeletal motor is an AAA+ protein and has a different structure and power stroke mechanism from those of myosins and kinesins. The linker protruding from the AAA+ ring of dynein swings according to the ATPase states, which, presumably, generates force to carry cargos within a cell. The communication mechanism between the track-binding and ATPase domains of dynein is unique because the two helices that presumably slide with respect to each other work as coordinators for these domains. Details of the mechanism underlying the power stroke and interdomain communication were revealed through recent progress in the structural studies of myosin, kinesin and dynein.

Allosteric modulation of cardiac myosin dynamics by omecamtiv mecarbil

New promising avenues for the pharmacological treatment of skeletal and heart muscle diseases rely on direct sarcomeric modulators, which are molecules that can directly bind to sarcomeric proteins and either inhibit or enhance their activity. A recent breakthrough has been the discovery of the myosin activator omecamtiv mecarbil (OM), which has been shown to increase the power output of the cardiac muscle and is currently in clinical trials for the treatment of heart failure. While the overall effect of OM on the mechano-chemical cycle of myosin is to increase the fraction of myosin molecules in the sarcomere that are strongly bound to actin, the molecular basis of its action is still not completely clear. We present here a Molecular Dynamics study of the motor domain of human cardiac myosin bound to OM, where the effects of the drug on the dynamical properties of the protein are investigated for the first time with atomistic resolution. We found that OM has a double effect on myosin dynamics, inducing a) an increased coupling of the motions of the converter and lever arm subdomains to the rest of the protein and b) a rewiring of the network of dynamic correlations, which produces preferential communication pathways between the OM binding site and distant functional regions. The location of the residues responsible for these effects suggests possible strategies for the future development of improved drugs and the targeting of specific cardiomyopathy-related mutations.

Molecular mechanisms underlying deoxy-ADP.Pi activation of pre-powerstroke myosin

Myosin activation is a viable approach to treat systolic heart failure. We previously demonstrated that striated muscle myosin is a promiscuous ATPase that can use most nucleoside triphosphates as energy substrates for contraction. When 2-deoxy ATP (dATP) is used, it acts as a myosin activator, enhancing cross-bridge binding and cycling. In vivo, we have demonstrated that elevated dATP levels increase basal cardiac function and rescues function of infarcted rodent and pig hearts. Here we investigate the molecular mechanism underlying this physiological effect. We show with molecular dynamics simulations that the binding of dADP.Pi (dATP hydrolysis products) to myosin alters the structure and dynamics of the nucleotide binding pocket, myosin cleft conformation, and actin binding sites, which collectively yield a myosin conformation that we predict favors weak, electrostatic binding to actin. In vitro motility assays at high ionic strength were conducted to test this prediction and we found that dATP increased motility. These results highlight alterations to myosin that enhance cross-bridge formation and reveal a potential mechanism that may underlie dATP-induced improvements in cardiac function.

Review: The ATPase mechanism of myosin and actomyosin

Myosins are a large family of molecular motors that use the common P-loop, Switch 1 and Switch 2 nucleotide binding motifs to recognize ATP, to create a catalytic site than can efficiently hydrolyze ATP and to communicate the state of the nucleotide pocket to other allosteric binding sites on myosin. The energy of ATP hydrolysis is used to do work against an external load. In this short review I will outline current thinking on the mechanism of ATP hydrolysis and how the energy of ATP hydrolysis is coupled to a series of protein conformational changes that allow a myosin, with the cytoskeleton track actin, to operate as a molecular motor of distinct types; fast movers, processive motors or strain sensors. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 483-491, 2016.

Product Release Pathways in Human and Plasmodium falciparum Phosphoribosyltransferase

Atomistic molecular dynamics (MD) simulations coupled with the metadynamics technique were carried out to delineate the product (PPi.2Mg and IMP) release mechanisms from the active site of both human (Hs) and Plasmodium falciparum (Pf) hypoxanthine-guanine-(xanthine) phosphoribosyltransferase (HG(X)PRT). An early movement of PPi.2Mg from its binding site has been observed. The swinging motion of the Asp side chain (D134/D145) in the binding pocket facilitates the detachment of IMP, which triggers the opening of flexible loop II, the gateway to the bulk solvent. In PfHGXPRT, PPi.2Mg and IMP are seen to be released via the same path in all of the biased MD simulations. In HsHGPRT too, the product molecules follow similar routes from the active site; however, an alternate but minor escape route for PPi.2Mg has been observed in the human enzyme. Tyr 104 and Phe 186 in HsHGPRT and Tyr 116 and Phe 197 in PfHGXPRT are the key residues that mediate the release of IMP, whereas the motion of PPi.2Mg away from the reaction center is guided by the negatively charged Asp and Glu and a few positively charged residues (Lys and Arg) that line the product release channels. Mutations of a few key residues present in loop II of Trypanosoma cruzi (Tc) HGPRT have been shown to reduce the catalytic efficiency of the enzyme. Herein, in silico mutation of corresponding residues in loop II of HsHGPRT and PfHGXPRT resulted in partial opening of the flexible loop (loop II), thus exposing the active site to bulk water, which offers a rationale for the reduced catalytic activity of these two mutant enzymes. Investigations of the product release from these HsHGPRT and PfHGXPRT mutants delineate the role of these important residues in the enzymatic turnover.

Temperature effect on the chemomechanical regulation of substeps within the power stroke of a single Myosin II

Myosin IIs in the skeletal muscle are highly efficient nanoscale machines evolved in nature. Understanding how they function can not only bring insights into various biological processes but also provide guidelines to engineer synthetic nanoscale motors working in the vicinity of thermal noise. Though it was clearly demonstrated that the behavior of a skeletal muscle fiber, or that of a single myosin was strongly affected by the temperature, how exactly the temperature affects the kinetics of a single myosin is not fully understood. By adapting the newly developed transitional state model, which successfully explained the intriguing motor force regulation during skeletal muscle contraction, here we systematically explain how exactly the power stroke of a single myosin proceeds, with the consideration of the chemomechanical regulation of sub-steps within the stroke. The adapted theory is then utilized to investigate the temperature effect on various aspects of the power stroke. Our analysis suggests that, though swing rates, the isometric force, and the maximal stroke size all strongly vary with the temperature, the temperature can have a very small effect on the releasable elastic energy within the power stroke.

A gating mechanism for Pi release governs the mRNA unwinding by eIF4AI during translation initiation

Eukaryotic translation initiation factor eIF4AI, the founding member of DEAD-box helicases, undergoes ATP hydrolysis-coupled conformational changes to unwind mRNA secondary structures during translation initiation. However, the mechanism of its coupled enzymatic activities remains unclear. Here we report that a gating mechanism for Pi release controlled by the inter-domain linker of eIF4AI regulates the coupling between ATP hydrolysis and RNA unwinding. Molecular dynamic simulations and experimental results revealed that, through forming a hydrophobic core with the conserved SAT motif of the N-terminal domain and I357 from the C-terminal domain, the linker gated the release of Pi from the hydrolysis site, which avoided futile hydrolysis cycles of eIF4AI. Further mutagenesis studies suggested this linker also plays an auto-inhibitory role in the enzymatic activity of eIF4AI, which may be essential for its function during translation initiation. Overall, our results reveal a novel regulatory mechanism that controls eIF4AI-mediated mRNA unwinding and can guide further mechanistic studies on other DEAD-box helicases.

How actin initiates the motor activity of Myosin

Fundamental to cellular processes are directional movements driven by molecular motors. A common theme for these and other molecular machines driven by ATP is that controlled release of hydrolysis products is essential for using the chemical energy efficiently. Mechanochemical transduction by myosin motors on actin is coupled to unknown structural changes that result in the sequential release of inorganic phosphate (Pi) and MgADP. We present here a myosin structure possessing an actin-binding interface and a tunnel (back door) that creates an escape route for Pi with a minimal rotation of the myosin lever arm that drives movements. We propose that this state represents the beginning of the powerstroke on actin and that Pi translocation from the nucleotide pocket triggered by actin binding initiates myosin force generation. This elucidates how actin initiates force generation and movement and may represent a strategy common to many molecular machines.

Structural dynamics of myosin 5 during processive motion revealed by interferometric scattering microscopy

Myosin 5a is a dual-headed molecular motor that transports cargo along actin filaments. By following the motion of individual heads with interferometric scattering microscopy at nm spatial and ms temporal precision we found that the detached head occupies a loosely fixed position to one side of actin from which it rebinds in a controlled manner while executing a step. Improving the spatial precision to the sub-nm regime provided evidence for an ångstrom-level structural transition in the motor domain associated with the power stroke. Simultaneous tracking of both heads revealed that consecutive steps follow identical paths to the same side of actin in a compass-like spinning motion demonstrating a symmetrical walking pattern. These results visualize many of the critical unknown aspects of the stepping mechanism of myosin 5 including head-head coordination, the origin of lever-arm motion and the spatiotemporal dynamics of the translocating head during individual steps.

Structural dynamics of myosin 5 during processive motion revealed by interferometric scattering microscopy

Myosin 5a is a dual-headed molecular motor that transports cargo along actin filaments. By following the motion of individual heads with interferometric scattering microscopy at nm spatial and ms temporal precision we found that the detached head occupies a loosely fixed position to one side of actin from which it rebinds in a controlled manner while executing a step. Improving the spatial precision to the sub-nm regime provided evidence for an ångstrom-level structural transition in the motor domain associated with the power stroke. Simultaneous tracking of both heads revealed that consecutive steps follow identical paths to the same side of actin in a compass-like spinning motion demonstrating a symmetrical walking pattern. These results visualize many of the critical unknown aspects of the stepping mechanism of myosin 5 including head-head coordination, the origin of lever-arm motion and the spatiotemporal dynamics of the translocating head during individual steps.

Trapping the ATP binding state leads to a detailed understanding of the F1-ATPase mechanism

The rotary motor enzyme FoF1-ATP synthase uses the proton-motive force across a membrane to synthesize ATP from ADP and Pi (H2PO4(-)) under cellular conditions that favor the hydrolysis reaction by a factor of 2 × 10(5). This remarkable ability to drive a reaction away from equilibrium by harnessing an external force differentiates it from an ordinary enzyme, which increases the rate of reaction without shifting the equilibrium. Hydrolysis takes place in the neighborhood of one conformation of the catalytic moiety F1-ATPase, whose structure is known from crystallography. By use of molecular dynamics simulations we trap a second structure, which is rotated by 40° from the catalytic dwell conformation and represents the state associated with ATP binding, in accord with single-molecule experiments. Using the two structures, we show why Pi is not released immediately after ATP hydrolysis, but only after a subsequent 120° rotation, in agreement with experiment. A concerted conformational change of the α3β3 crown is shown to induce the 40° rotation of the γ-subunit only when the βE subunit is empty, whereas with Pi bound, βE serves as a latch to prevent the rotation of γ. The present results provide a rationalization of how F1-ATPase achieves the coupling between the small changes in the active site of βDP and the 40° rotation of γ.

A new mechanokinetic model for muscle contraction, where force and movement are triggered by phosphate release

The atomic structure of myosin-S1 suggests that its working stroke, which generates tension and shortening in muscle, is triggered by the release of inorganic phosphate from the active site. This mechanism is the basis of a new mechanokinetic model for contractility, using the biochemical actomyosin ATPase cycle, strain-dependent kinetics and dimeric myosins on buckling rods. In this model, phosphate-dependent aspects of contractility arise from a rapid reversible release of phosphate from the initial bound state (A.M.ADP.Pi), which triggers the stroke. Added phosphate drives bound myosin towards this initial state, and the transient tension response to a phosphate jump reflects the rate at which it detaches from actin. Predictions for the tensile and energetic properties of striated muscle as a function of phosphate level, including the tension responses to length steps and Pi-jumps, are compared with experimental data from rabbit psoas fibres at 10 °C. The phosphate sensitivity of isometric tension is maximal when the actin affinity of M.ADP.Pi is near unity. Hence variations in actin affinity modulate the phosphate dependence of isometric tension, and may explain why phosphate sensitivity is temperature-dependent or absent in different muscles.

Effects of ATP and actin-filament binding on the dynamics of the myosin II S1 domain

Actin and myosin interact with one another to perform a variety of cellular functions. Central to understanding the processive motion of myosin on actin is the characterization of the individual states along the mechanochemical cycle. We present an all-atom molecular dynamics simulation of the myosin II S1 domain in the rigor state interacting with an actin filament. We also study actin-free myosin in both rigor and post-rigor conformations. Using all-atom level and coarse-grained analysis methods, we investigate the effects of myosin binding on actin, and of actin binding on myosin. In particular, we determine the domains of actin and myosin that interact strongly with one another at the actomyosin interface using a highly coarse-grained level of resolution, and we identify a number of salt bridges and hydrogen bonds at the interface of myosin and actin. Applying coarse-grained analysis, we identify differences in myosin states dependent on actin-binding, or ATP binding. Our simulations also indicate that the actin propeller twist-angle and nucleotide cleft-angles are influenced by myosin at the actomyosin interface. The torsional rigidity of the myosin-bound filament is also calculated, and is found to be increased compared to previous simulations of the free filament.

Phosphate release coupled to rotary motion of F1-ATPase

F1-ATPase, the catalytic domain of ATP synthase, synthesizes most of the ATP in living organisms. Running in reverse powered by ATP hydrolysis, this hexameric ring-shaped molecular motor formed by three αβ-dimers creates torque on its central γ-subunit. This reverse operation enables detailed explorations of the mechanochemical coupling mechanisms in experiment and simulation. Here, we use molecular dynamics simulations to construct a first atomistic conformation of the intermediate state following the 40° substep of rotary motion, and to study the timing and molecular mechanism of inorganic phosphate (Pi) release coupled to the rotation. In response to torque-driven rotation of the γ-subunit in the hydrolysis direction, the nucleotide-free αβE interface forming the "empty" E site loosens and singly charged Pi readily escapes to the P loop. By contrast, the interface stays closed with doubly charged Pi. The γ-rotation tightens the ATP-bound αβTP interface, as required for hydrolysis. The calculated rate for the outward release of doubly charged Pi from the αβE interface 120° after ATP hydrolysis closely matches the ~1-ms functional timescale. Conversely, Pi release from the ADP-bound αβDP interface postulated in earlier models would occur through a kinetically infeasible inward-directed pathway. Our simulations help reconcile conflicting interpretations of single-molecule experiments and crystallographic studies by clarifying the timing of Pi exit, its pathway and kinetics, associated changes in Pi protonation, and changes of the F1-ATPase structure in the 40° substep. Important elements of the molecular mechanism of Pi release emerging from our simulations appear to be conserved in myosin despite the different functional motions.

1,3-Diethylurea-enhanced Mg-ATPase activity of skeletal muscle myosin with a converse effect on the sliding motility

We investigate the effects of urea and its derivatives on the ATPase activity and on the in vitro motility of chicken skeletal muscle actomyosin. Mg-ATPase rate of myosin subfragment-1 (S1) is increased by 4-fold by 0.3M 1,3-diethylurea (DEU), but it is unaffected by urea, thiourea, and 1,3-dimethylurea at ≤1M concentration. Thus, we further examine the effects of DEU in comparison to those of urea as reference. In in vitro motility assay, we find that in the presence of 0.3M DEU, the sliding speeds of actin filaments driven by myosin and heavy meromyosin (HMM) are significantly decreased to 1/16 and 1/6.6, respectively, compared with the controls. However, the measurement of the actin-activated ATPase activity of HMM shows that the maximal rate, Vmax, is almost unchanged with DEU. Thus, the myosin-driven sliding motility of actin filaments is significantly impeded in the presence of 0.3M DEU, whereas the cyclic interaction of myosin with F-actin occurs during the ATP turnover, the rate of which is close to that without DEU. In contrast to DEU, 0.3M urea exhibits only modest effects on both actin-activated ATPase and sliding motility of actomyosin. Thus, DEU has the effect of uncoupling the sliding motility of actomyosin from its ATP turnover.

The myosin start-of-power stroke state and how actin binding drives the power stroke

We propose that on binding to actin at the start of the power stroke the myosin cross-bridge takes on the rigor configuration at the actin interface. Starting from the prepower stroke state, this can be achieved by a small movement (16° rotation) of the lower 50K domain without twisting the central β-sheet or opening switch-1 or switch-2. The movement of the lower 50K domain puts a strain on the W-helix. This strain tries to twist the β-sheet, which could drive the power stroke. This would provide a coupling between actin binding and the execution of the power stroke. During the power stroke the β-sheet twists, moving the P-loop away from switch-2, which opens the nucleotide binding pocket and separates ADP from Pi . The power stroke is different from the recovery stroke because the upper and lower 50K domains are tethered in the rigor configuration.

Molecular dynamics simulations of isoleucine-release pathway in GAF domain of N-CodY from Bacillus Subtilis

The GAF domain located in the N-terminal motifs of CodY (N-CodY) is responsible for increasing the affinity of CodY to its target sites on DNA by its interaction with the branched chain amino acids (BCAAs) involving isoleucine, leucine and valine. The study of the interaction of GAF domain with isoleucine gains much attention in recent years, but the mechanism of isoleucine release still remains unclear. In this paper, a conventional molecular dynamics (MD) and force probe molecular dynamics (FPMD) simulations have been performed with the aim to understand how the isoleucine ligand escapes from the GAF domain of N-CodY from Bacillus subtilis. The MD results reveal that the ligand release is a gradual process, which is accompanied by the movement of the loop between β3 and β4. During the periods of ligand escaping from the bottom to the top of binding pocket, isoleucine forms hydrogen bonds one after another with series of residues, such as ARG61, THR96, PHE98, VAL100, GLU101 and ASN102, under the mediation of hydrophobic contacts. The FPMD results show that the easiest way to pull ligand out of the cavity is along x direction (i.e. the direction is opposite to MET62).

Human Tonic and Phasic Smooth Muscle Myosin Isoforms Are Unresponsive to the Loop 1 Insert

Smooth muscle myosin gene products include two isoforms, SMA and SMB, differing by a 7-residue peptide in loop 1 (i7) at the myosin active site where ATP is hydrolyzed. Using chicken isoforms, previous work indicated that the i7 deletion in SMA prolongs strong actin binding by inhibiting active site ingress and egress of nucleotide when compared to i7 inserted SMB. Additionally, i7 deletion inhibits Pi release associated with the switch 2 closed → open transition in actin-activated ATPase. Switch 2 is far from loop 1 indicating i7 deletion has an allosteric effect on Pi release. Chicken SMA and SMB have unknown and robust nucleotide-sensitive tryptophan (NST) fluorescence increments, respectively. Human SMA and SMB both lack NST increments while Pi release in Ca²⁺ ATPase is not impacted by i7 deletion. The NST reports relay helix movement following conformation change in switch 2 but in the open → closed transition. The NST is common to all known myosin isoforms except human smooth muscle. Other independent works on human SMA and SMB motility indicate no functional effect of i7 deletion. Smooth muscle myosin is a stunning example of species-specific myosin structure/function divergence underscoring the danger in extrapolating disease-linked mutant effects on myosin across species.

The mechanism of the converter domain rotation in the recovery stroke of myosin motor protein

Upon ATP binding, myosin motor protein is found in two alternative conformations, prerecovery state M* and postrecovery state M**. The transition from one state to the other, known as the recovery stroke, plays a key role in the myosin functional cycle. Despite much recent research, the microscopic details of this transition remain elusive. A critical step in the recovery stroke is the rotation of the converter domain from "up" position in prerecovery state to "down" position in postrecovery state that leads to the swing of the lever arm attached to it. In this work, we demonstrate that the two rotational states of the converter domain are determined by the interactions within a small structural motif in the force-generating region of the protein that can be accurately modeled on computers using atomic representation and explicit solvent. Our simulations show that the transition between the two states is controlled by a small helix (SH1) located next to the relay helix and relay loop. A small translation in the position of SH1 away from the relay helix is seen to trigger the transition from "up" state to "down" state. The transition is driven by a cluster of hydrophobic residues I687, F487, and F506 that make significant contributions to the stability of both states. The proposed mechanism agrees well with the available structural and mutational studies.

Complex intramolecular mechanics of G-actin--an elastic network study

Systematic numerical investigations of conformational motions in single actin molecules were performed by employing a simple elastic-network (EN) model of this protein. Similar to previous investigations for myosin, we found that G-actin essentially behaves as a strain sensor, responding by well-defined domain motions to mechanical perturbations. Several sensitive residues within the nucleotide-binding pocket (NBP) could be identified, such that the perturbation of any of them can induce characteristic flattening of actin molecules and closing of the cleft between their two mobile domains. Extending the EN model by introduction of a set of breakable links which become effective only when two domains approach one another, it was observed that G-actin can possess a metastable state corresponding to a closed conformation and that a transition to this state can be induced by appropriate perturbations in the NBP region. The ligands were roughly modeled as a single particle (ADP) or a dimer (ATP), which were placed inside the NBP and connected by elastic links to the neighbors. Our approximate analysis suggests that, when ATP is present, it stabilizes the closed conformation of actin. This may play an important role in the explanation why, in the presence of ATP, the polymerization process is highly accelerated.

Recent insights into muscle fatigue at the cross-bridge level

The depression in force and/or velocity associated with muscular fatigue can be the result of a failure at any level, from the initial events in the motor cortex of the brain to the formation of an actomyosin cross-bridge in the muscle cell. Since all the force and motion generated by muscle ultimately derives from the cyclical interaction of actin and myosin, researchers have focused heavily on the impact of the accumulation of intracellular metabolites [e.g., P(i), H(+) and adenosine diphoshphate (ADP)] on the function these contractile proteins. At saturating Ca(++) levels, elevated P(i) appears to be the primary cause for the loss in maximal isometric force, while increased [H(+)] and possibly ADP act to slow unloaded shortening velocity in single muscle fibers, suggesting a causative role in muscular fatigue. However the precise mechanisms through which these metabolites might affect the individual function of the contractile proteins remain unclear because intact muscle is a highly complex structure. To simplify problem isolated actin and myosin have been studied in the in vitro motility assay and more recently the single molecule laser trap assay with the findings showing that both P(i) and H(+) alter single actomyosin function in unique ways. In addition to these new insights, we are also gaining important information about the roles played by the muscle regulatory proteins troponin (Tn) and tropomyosin (Tm) in the fatigue process. In vitro studies, suggest that both the acidosis and elevated levels of P(i) can inhibit velocity and force at sub-saturating levels of Ca(++) in the presence of Tn and Tm and that this inhibition can be greater than that observed in the absence of regulation. To understand the molecular basis of the role of regulatory proteins in the fatigue process researchers are taking advantage of modern molecular biological techniques to manipulate the structure and function of Tn/Tm. These efforts are beginning to reveal the relevant structures and how their functions might be altered during fatigue. Thus, it is a very exciting time to study muscle fatigue because the technological advances occurring in the fields of biophysics and molecular biology are providing researchers with the ability to directly test long held hypotheses and consequently reshaping our understanding of this age-old question.

Novel therapies in acute and chronic heart failure

Despite past advances in the pharmacological management of heart failure, the prognosis of these patients remains poor, and for many, treatment options remain unsatisfactory. Additionally, the treatments and clinical outcomes of patients with acute decompensated heart failure have not changed substantially over the past few decades. Consequently, there is a critical need for new drugs that can improve clinical outcomes. In the setting of acute heart failure, new inotrops such as cardiac myosin activators and new vasodilators such as relaxin have been developed. For chronic heart failure with reduced ejection fraction, there are several new approaches that target multiple pathophysiological mechanism including novel blockers of the renin-angiotensin-aldosterone system (direct renin inhibitors, dual-acting inhibitors of the angiotensin II receptor and neprilysin, aldosterone synthase inhibitors), ryanodine receptor stabilizers, and SERCA activators. Heart failure with preserved ejection fraction represents a substantial therapeutic problem as no therapy has been demonstrated to improve symptoms or outcomes in this condition. Newer treatment strategies target specific structural and functional abnormalities that lead to increased myocardial stiffness. Dicarbonyl-breaking compounds reverse advanced glycation-induced cross-linking of collagen and improve the compliance of aged and/or diabetic myocardium. Modulation of titin-dependent passive tension can be achieved via phosphorylation of a unique sequence on the extensible region of the protein. This review describes the pathophysiological basis, mechanism of action, and available clinical efficacy data of drugs that are currently under development. Finally, new therapies for the treatment of heart failure complications, such as pulmonary hypertension and anemia, are discussed.

Interactions between relay helix and Src homology 1 (SH1) domain helix drive the converter domain rotation during the recovery stroke of myosin II

Myosin motor protein exists in two alternative conformations, prerecovery state M* and postrecovery state M**, on adenosine triphosphate binding. The details of the M*-to-M** transition, known as the recovery stroke to reflect its role as the functional opposite of the force-generating power stroke, remain elusive. The defining feature of the postrecovery state is a kink in the relay helix, a key part of the protein involved in force generation. In this article, we determine the interactions that are responsible for the appearance of the kink. We design a series of computational models that contain three other segments, relay loop, converter domain, and Src homology 1 (SH1) domain helix, with which relay helix interacts and determine their structure in accurate replica exchange molecular dynamics simulations in explicit solvent. By conducting an exhaustive combinatorial search among different models, we find that: (1) the converter domain must be attached to the relay helix during the transition, so it does not interfere with other parts of the protein and (2) the structure of the relay helix is controlled by SH1 helix. The kink is strongly coupled to the position of SH1 helix. It arises as a result of direct interactions between SH1 and the relay helix and leads to a rotation of the C-terminal part of the relay helix, which is subsequently transmitted to the converter domain.

Myosin-V as a mechanical sensor: an elastic network study

According to recent experiments, the molecular-motor myosin behaves like a strain sensor, exhibiting different functional responses when loads in opposite directions are applied to its tail. Within an elastic-network model, we explore the sensitivity of the protein to the forces acting on the tail and find, in agreement with experiments, that such forces invoke conformational changes that should affect filament binding and ADP release. Furthermore, conformational responses of myosin to the application of forces to individual residues in its principal functional regions are systematically investigated and a detailed sensitivity map of myosin-V is thus obtained. The results suggest that the strain-sensor behavior is involved in the intrinsic operation of this molecular motor.

Inter-domain movements in polyketide synthases: a molecular dynamics study

Insights into the structure and dynamics of modular polyketide synthases (PKS) are essential for understanding the mechanistic details of the biosynthesis of a large number of pharmaceutically important secondary metabolites. The crystal structures of the KS-AT di-domain from erythromycin synthase have revealed the relative orientation of various catalytic domains in a minimal PKS module. However, the relatively large distance between catalytic centers of KS and AT domains in the static structure has posed certain intriguing questions regarding mechanistic details of substrate transfer during polyketide biosynthesis. In order to investigate the role of inter-domain movements in substrate channeling, we have carried out a series of explicit solvent MD simulations for time periods ranging from 10 to 15 ns on the KS-AT di-domain and its sub-fragments. Analyses of these MD trajectories have revealed that both the catalytic domains and the structured inter-domain linker region remain close to their starting structures. Inter-domain movements at KS-linker and linker-AT interfaces occur around hinge regions which connect the structured linker region to the catalytic domains. The KS-linker interface was found to be more flexible compared to the linker-AT interface. However, inter-domain movements observed during the timescale of our simulations do not significantly reduce the distance between catalytic centers of KS and AT domains for facilitating substrate channeling. Based on these studies and prediction of intrinsic disorder we propose that the intrinsically unstructured linker stretch preceding the ACP domain might be facilitating movement of ACP domains to various catalytic centers.

Extension of a three-helix bundle domain of myosin VI and key role of calmodulins

The molecular motor protein myosin VI moves toward the minus-end of actin filaments with a step size of 30-36 nm. Such large step size either drastically limits the degree of complex formation between dimer subunits to leave enough length for the lever arms, or requires an extension of the lever arms' crystallographically observed structure. Recent experimental work proposed that myosin VI dimerization triggers the unfolding of the protein's proximal tail domain which could drive the needed lever-arm extension. Here, we demonstrate through steered molecular dynamics simulation the feasibility of sufficient extension arising from turning a three-helix bundle into a long α-helix. A key role is played by the known calmodulin binding that facilitates the extension by altering the strain path in myosin VI. Sequence analysis of the proximal tail domain suggests that further calmodulin binding sites open up when the domain's three-helix bundle is unfolded and that subsequent calmodulin binding stabilizes the extended lever arms.

A conformational transition in the myosin VI converter contributes to the variable step size

Myosin VI (MVI) is a dimeric molecular motor that translocates backwards on actin filaments with a surprisingly large and variable step size, given its short lever arm. A recent x-ray structure of MVI indicates that the large step size can be explained in part by a novel conformation of the converter subdomain in the prepowerstroke state, in which a 53-residue insert, unique to MVI, reorients the lever arm nearly parallel to the actin filament. To determine whether the existence of the novel converter conformation could contribute to the step-size variability, we used a path-based free-energy simulation tool, the string method, to show that there is a small free-energy difference between the novel converter conformation and the conventional conformation found in other myosins. This result suggests that MVI can bind to actin with the converter in either conformation. Models of MVI/MV chimeric dimers show that the variability in the tilting angle of the lever arm that results from the two converter conformations can lead to step-size variations of ∼12 nm. These variations, in combination with other proposed mechanisms, could explain the experimentally determined step-size variability of ∼25 nm for wild-type MVI. Mutations to test the findings by experiment are suggested.

Coarse-grained modeling of conformational transitions underlying the processive stepping of myosin V dimer along filamentous actin

To explore the structural basis of processive stepping of myosin V along filamentous actin, we have performed comprehensive modeling of its key conformational states and transitions with an unprecedented residue level of details. We have built structural models for a myosin V monomer complexed with filamentous actin at four biochemical states [adenosine diphosphate (ATP)-, adenosine diphosphate (ADP)-phosphate-, ADP-bound or nucleotide-free]. Then we have modeled a myosin V dimer (consisting of lead and rear head) at various two-head-bound states with nearly straight lever arms rotated by intramolecular strain. Next, we have performed transition pathway modeling to determine the most favorable sequence of transitions (namely, phosphate release at the lead head followed by ADP release at the rear head, while ADP release at the lead head is inhibited), which underlie the kinetic coordination between the two heads. Finally, we have used transition pathway modeling to reveal the order of structural changes during three key biochemical transitions (phosphate release at the lead head, ADP release and ATP binding at the rear head), which shed lights on the strain-dependence of the allosterically coupled motions at various stages of myosin V's work cycle. Our modeling results are in agreement with and offer structural insights to many results of kinetic, single-molecule and structural studies of myosin V.

Phosphate enhances myosin-powered actin filament velocity under acidic conditions in a motility assay

Elevated levels of inorganic phosphate (P(i)) are believed to inhibit muscular force by reversing myosin's force-generating step. These same levels of P(i) can also affect muscle velocity, but the molecular basis underlying these effects remains unclear. We directly examined the effect of P(i) (30 mM) on skeletal muscle myosin's ability to translocate actin (V(actin)) in an in vitro motility assay. Manipulation of the pH enabled us to probe rebinding of P(i) to myosin's ADP-bound state, while changing the ATP concentration probed rebinding to the rigor state. Surprisingly, the addition of P(i) significantly increased V(actin) at both pH 6.8 and 6.5, causing a doubling of V(actin) at pH 6.5. To probe the mechanisms underlying this increase in speed, we repeated these experiments while varying the ATP concentration. At pH 7.4, the effects of P(i) were highly ATP dependent, with P(i) slowing V(actin) at low ATP (<500 μM), but with a minor increase at 2 mM ATP. The P(i)-induced slowing of V(actin), evident at low ATP (pH 7.4), was minimized at pH 6.8 and completely reversed at pH 6.5. These data were accurately fit with a simple detachment-limited kinetic model of motility that incorporated a P(i)-induced prolongation of the rigor state, which accounted for the slowing of V(actin) at low ATP, and a P(i)-induced detachment from a strongly bound post-power-stroke state, which accounted for the increase in V(actin) at high ATP. These findings suggest that P(i) differentially affects myosin function: enhancing velocity, if it rebinds to the ADP-bound state, while slowing velocity, if it binds to the rigor state.

Simulations of allosteric transitions

Allosteric transitions are one of the subtler mechanisms used by nature to fine tune protein activity. Effector binding to a specific site on the protein surface induces significant activity change, and initiates a conformational transition that frequently includes domain motions and is very large. From a theoretical and biophysical perspective two problems are particularly intriguing. The first is the way in which a launching signal, which is spatially confined and includes only a few interacting atoms, is propagated to a large-scale conformational transition we frequently see in allosteric transitions. Hence, there is the question of how a small perturbation is magnified to yield motions of thousands of atoms. The second puzzle is of focus, coherence, and efficiency. The impact of the binding of the effector is sometimes extended over tens of angstroms. How the signal is transmitted and kept significant over such large distances in the 'noisy' molecular environment is another major direction of investigation. In the present review we examined different theoretical and computational attempts to solve the questions.