Mechanochemical coupling in the myosin motor domain. I. Insights from equilibrium active-site simulations

Although the major structural transitions in molecular motors are often argued to couple to the binding of Adenosine triphosphate (ATP), the recovery stroke in the conventional myosin has been shown to be dependent on the hydrolysis of ATP. To obtain a clearer mechanistic picture for such "mechanochemical coupling" in myosin, equilibrium active-site simulations with explicit solvent have been carried out to probe the behavior of the motor domain as functions of the nucleotide chemical state and conformation of the converter/relay helix. In conjunction with previous studies of ATP hydrolysis with different active-site conformations and normal mode analysis of structural flexibility, the results help establish an energetics-based framework for understanding the mechanochemical coupling. It is proposed that the activation of hydrolysis does not require the rotation of the lever arm per se, but the two processes are tightly coordinated because both strongly couple to the open/close transition of the active site. The underlying picture involves shifts in the dominant population of different structural motifs as a consequence of changes elsewhere in the motor domain. The contribution of this work and the accompanying paper [] is to propose the actual mechanism behind these "population shifts" and residues that play important roles in the process. It is suggested that structural flexibilities at both the small and large scales inherent to the motor domain make it possible to implement tight couplings between different structural motifs while maintaining small free-energy drops for processes that occur in the detached states, which is likely a feature shared among many molecular motors. The significantly different flexibility of the active site in different X-ray structures with variable level arm orientations supports the notation that external force sensed by the lever arm may transmit into the active site and influence the chemical steps (nucleotide hydrolysis and/or binding).     

Extensive conformational transitions are required to turn on ATP hydrolysis in myosin

Conventional myosin is representative of biomolecular motors in which the hydrolysis of adenosine triphosphate (ATP) is coupled to large-scale structural transitions both in and remote from the active site. The mechanism that underlies such “mechanochemical coupling,” especially the causal relationship between hydrolysis and allosteric structural changes, has remained elusive despite extensive experimental and computational analyses. In this study, using combined quantum mechanical and molecular mechanical simulations and different conformations of the myosin motor domain, we provide evidence to support that regulation of ATP hydrolysis activity is not limited to residues in the immediate environment of the phosphate. Specifically, we illustrate that efficient hydrolysis of ATP depends not only on the proper orientation of the lytic water but also on the structural stability of several nearby residues, especially the Arg238-Glu459 salt bridge (the numbering of residues follows myosin II in Dictyostelium discoideum) and the water molecule that spans this salt bridge and the lytic water. More importantly, by comparing the hydrolysis activities in two motor conformations with very similar active-site (i.e., Switches I and II) configurations, which distinguished this work from our previous study, the results clearly indicate that the ability of these residues to perform crucial electrostatic stabilization relies on the configuration of residues in the nearby N-terminus of the relay helix and the “wedge loop.” Without the structural support from those motifs, residues in a closed active site in the post-rigor motor domain undergo subtle structural variations that lead to consistently higher calculated ATP hydrolysis barriers than in the pre-powerstroke state. In other words, starting from the post-rigor state, turning on the ATPase activity requires not only displacement of Switch II to close the active site but also structural transitions in the N-terminus of the relay helix and the “wedge loop,” which have been proposed previously to be ultimately coupled to the rotation of the converter subdomain 40 Å away.     

Computational Tool for Ensemble Averaging of Single-Molecule Data

Molecular motors couple chemical transitions to conformational changes that perform mechanical work in a wide variety of biological processes. Disruption of this coupling can lead to diseases, and therefore there is a need to accurately measure mechanochemical coupling in motors in both health and disease. Optical tweezers with nanometer spatial and millisecond temporal resolution have provided valuable insights into these processes. However, fluctuations due to Brownian motion can make it difficult to precisely resolve these conformational changes. One powerful analysis technique that has improved our ability to accurately measure mechanochemical coupling in motor proteins is ensemble averaging of individual trajectories. Here, we present a user-friendly computational tool, Software for Precise Analysis of Single Molecules (SPASM), for generating ensemble averages of single-molecule data. This tool utilizes several conceptual advances, including optimized procedures for identifying single-molecule interactions and the implementation of a change-point algorithm, to more precisely resolve molecular transitions. Using both simulated and experimental data, we demonstrate that these advances allow for accurate determination of the mechanics and kinetics of the myosin working stroke with a smaller set of data. Importantly, we provide our open-source MATLAB-based program with a graphical user interface that enables others to readily apply these advances to the analysis of their own data.     

The principal motions involved in the coupling mechanism of the recovery stroke of the myosin motor

Muscle contraction is driven by a cycle of conformational changes in the myosin II head. After myosin binds ATP and releases from the actin fibril, myosin prepares for the next power stroke by rotating back the converter domain that carries the lever arm by 60 degrees . This recovery stroke is coupled to the activation of myosin ATPase by a mechanism that is essential for an efficient motor cycle. The mechanics of this coupling have been proposed to occur via two distinct and successive motions of the two helices that hold the converter domain: in a first phase a seesaw motion of the relay helix, followed by a piston-like motion of the SH1 helix in a second phase. To test this model, we have determined the principal motions of these structural elements during equilibrium molecular dynamics simulations of the crystallographic end states of the recovery-stroke by using principal component analysis. This reveals that the only principal motions of these two helices that make a large-amplitude contribution towards the conformational change of the recovery stroke are indeed the predicted seesaw and piston motions. Moreover, the results demonstrate that the seesaw motion of the relay helix dominates in the dynamics of the pre-recovery stroke structure, but not in the dynamics of the post-recovery stroke structure, and vice versa for the piston motion of the SH1 helix. This is consistent with the order of the proposed two-phase model for the coupling mechanism of the recovery stroke. Molecular movies of these principal motions are available at http://www.iwr.uni-heidelberg.de/groups/biocomp/fischer.     

Mobility of Molecular Motors Regulates Contractile Behaviors of Actin Networks

Cells generate mechanical forces primarily from interactions between F-actin, cross-linking proteins, myosin motors, and other actin-binding proteins in the cytoskeleton. To understand how molecular interactions between the cytoskeletal elements generate forces, a number of in vitro experiments have been performed but are limited in their ability to accurately reproduce the diversity of motor mobility. In myosin motility assays, myosin heads are fixed on a surface and glide F-actin. By contrast, in reconstituted gels, the motion of both myosin and F-actin is unrestricted. Because only these two extreme conditions have been used, the importance of mobility of motors for network behaviors has remained unclear. In this study, to illuminate the impacts of motor mobility on the contractile behaviors of the actin cytoskeleton, we employed an agent-based computational model based on Brownian dynamics. We find that if motors can bind to only one F-actin like myosin I, networks are most contractile at intermediate mobility. In this case, less motor mobility helps motors stably pull F-actins to generate tensile forces, whereas higher motor mobility allows F-actins to aggregate into larger clustering structures. The optimal intermediate motor mobility depends on the stall force and affinity of motors that are regulated by mechanochemical rates. In addition, we find that the role of motor mobility can vary drastically if motors can bind to a pair of F-actins. A network can exhibit large contraction with high motor mobility because motors bound to antiparallel pairs of F-actins can exert similar forces regardless of their mobility. Results from this study imply that the mobility of molecular motors may critically regulate contractile behaviors of actin networks in cells.     

Making sense of melanosome dynamics in mouse melanocytes

Molecular motors drive most if not all organelle movements in Eukaryotic cells. These proteins are thought to bind to the organelle surface and, through the action of their mechanochemical domains, to translocate the organelle along a cytoskeletal track. In the case of the myosin family of molecular motors, the cytoskeletal track is filamentous actin. Microtubules serve as the cytoskeletal track for the kinesins and dyneins. While a considerable amount is known about the motors and tracks responsible for the bi-directional movement of pigment granules in fish and frog melanophores, relatively little is known about how melanosomes in mammalian melanocytes are transported out the cells dendritic arbor, accumulated at the ends of these dendrites, and transferred to keratinocytes. In this short review, we focus on the use of video microscopy to address these questions in mouse melanocytes, and we describe how an analysis of melanosome dynamics within wild type and dilute melanocytes shaped our thinking regarding the role of an unconventional myosin in melanosome transport and distribution.     

Exploration of the conformational space of myosin recovery stroke via molecular dynamics

Muscle contractions are driven by cyclic conformational changes of myosin, whose molecular mechanisms of operation are being elucidated by recent advances in crystallographic studies and single molecule experiments. To complement such structural studies and consider the energetics of the conformational changes of myosin head, umbrella sampling molecular dynamics (MD) simulations were performed with the all-atom model of the scallop myosin sub-fragment 1 (S1) with a bound ATP in solution in explicit water using the crystallographic near-rigor and transition state conformations as two references. The constraints on RMSD reaction coordinates used for the umbrella sampling were found to steer the conformational changes efficiently, and relatively close correlations have been observed between the set of characteristic structural changes including the lever arm rotation and the closing of the nucleotide binding pocket. The lever arm angle and key residue interaction distances in the nucleotide binding pocket and the relay helix show gradual changes along the recovery stroke reaction coordinate, consistent with previous crystallographic and computational minimum energy studies. Thermal fluctuations, however, appear to make the switch-2 coordination of ATP more flexible than suggested by crystal structures. The local solvation environment of the fluorescence probe, Trp 507 (scallop numbering), also appears highly mobile in the presence of thermal fluctuations.     

Early stages of the recovery stroke in myosin II studied by molecular dynamics simulations

The recovery stroke is a key step in the functional cycle of muscle motor protein myosin, during which pre-recovery conformation of the protein is changed into the active post-recovery conformation, ready to exersice force. We study the microscopic details of this transition using molecular dynamics simulations of atomistic models in implicit and explicit solvent. In more than 2 μs of aggregate simulation time, we uncover evidence that the recovery stroke is a two-step process consisting of two stages separated by a time delay. In our simulations, we directly observe the first stage at which switch II loop closes in the presence of adenosine triphosphate at the nucleotide binding site. The resulting configuration of the nucleotide binding site is identical to that detected experimentally. Distribution of inter-residue distances measured in the force generating region of myosin is in good agreement with the experimental data. The second stage of the recovery stroke structural transition, rotation of the converter domain, was not observed in our simulations. Apparently it occurs on a longer time scale. We suggest that the two parts of the recovery stroke need to be studied using separate computational models.     

A deterministic mechanism producing the loose coupling phenomenon observed in an actomyosin system

Myosins are molecular motors that convert the chemical energy of ATP into mechanical work called a power stroke. Class II myosin engaged in muscle contraction is reported to show a "loose coupling phenomenon", in which the number of power strokes is greater than the number of ATP hydrolyses. This phenomenon cannot be explained by the lever-arm hypothesis, which is currently accepted as a standard theory for myosin motility. In this paper, a model is proposed to reproduce the loose coupling phenomenon. The model is based on a mechanochemical process called "Driven by Detachment (DbD)" mechanism, which assumes that the energy of the power strokes originates from the potential energy generated by the attractive force between myosin and actin. During the docking process, the potential energy is converted into an intramolecular strain in a myosin molecule, which drives the power stroke after the myosin is firmly attached to an actin filament. The energy of ATP is used to temporarily reduce the attractive force and to increase the potential energy. Therefore, it is not directly linked to the power strokes. When myosin molecules work as an aggregate, the sliding movement of a myosin filament driven by the power strokes of some myosin heads makes other myosin heads that have completed their power strokes detach from the actin without consuming ATP. Under the DbD mechanism, these passively detached myosins can be again engaged in power strokes after the next attachment to actin. As a result, the number of power strokes becomes greater than the number of ATP hydrolyses, and the loose coupling phenomenon will be observed. A theoretical analysis indicates that the efficiency of converting the potential energy into intramolecular elastic energy determines the number of power strokes per each ATP hydrolysis. Computer simulations showed that the DbD mechanism actually produced the loose coupling phenomenon. A critical requirement for this mechanism is that ATP must preferentially facilitate the detachment of myosins that have completed their power strokes, but are still strongly attached to the actin. This requirement may be fulfilled by ATP hydrolysis tightly depending on the conformation of a myosin molecule.     

Cardiac ventricular myosin and slow skeletal myosin exhibit dissimilar chemomechanical properties despite bearing the same myosin heavy chain isoform

The myosin II motors are ATP-powered force-generating machines driving cardiac and muscle contraction. Myosin II heavy chain isoform-beta (β-MyHC) is primarily expressed in the ventricular myocardium and in slow-twitch muscle fibers, such as M. soleus. M. soleus–derived myosin II (SolM-II) is often used as an alternative to the ventricular β-cardiac myosin (βM-II); however, the direct assessment of biochemical and mechanical features of the native myosins is limited. By employing optical trapping, we examined the mechanochemical properties of native myosins isolated from the rabbit heart ventricle and soleus muscles at the single-molecule level. We found purified motors from the two tissue sources, despite expressing the same MyHC isoform, displayed distinct motile and ATPase kinetic properties. We demonstrate βM-II was approximately threefold faster in the actin filament–gliding assay than SolM-II. The maximum actomyosin (AM) detachment rate derived in single-molecule assays was also approximately threefold higher in βM-II, while the power stroke size and stiffness of the “AM rigor” crossbridge for both myosins were comparable. Our analysis revealed a higher AM detachment rate for βM-II, corresponding to the enhanced ADP release rates from the crossbridge, likely responsible for the observed differences in the motility driven by these myosins. Finally, we observed a distinct myosin light chain 1 isoform (MLC1sa) that associates with SolM-II, which might contribute to the observed kinetics differences between βM-II and SolM-II. These results have important implications for the choice of tissue sources and justify prerequisites for the correct myosin heavy and light chains to study cardiomyopathies.     

Direct observation of the myosin Va recovery stroke that contributes to unidirectional stepping along actin

Myosins are ATP-driven linear molecular motors that work as cellular force generators, transporters, and force sensors. These functions are driven by large-scale nucleotide-dependent conformational changes, termed "strokes"; the "power stroke" is the force-generating swinging of the myosin light chain-binding "neck" domain relative to the motor domain "head" while bound to actin; the "recovery stroke" is the necessary initial motion that primes, or "cocks," myosin while detached from actin. Myosin Va is a processive dimer that steps unidirectionally along actin following a "hand over hand" mechanism in which the trailing head detaches and steps forward ～72 nm. Despite large rotational Brownian motion of the detached head about a free joint adjoining the two necks, unidirectional stepping is achieved, in part by the power stroke of the attached head that moves the joint forward. However, the power stroke alone cannot fully account for preferential forward site binding since the orientation and angle stability of the detached head, which is determined by the properties of the recovery stroke, dictate actin binding site accessibility. Here, we directly observe the recovery stroke dynamics and fluctuations of myosin Va using a novel, transient caged ATP-controlling system that maintains constant ATP levels through stepwise UV-pulse sequences of varying intensity. We immobilized the neck of monomeric myosin Va on a surface and observed real time motions of bead(s) attached site-specifically to the head. ATP induces a transient swing of the neck to the post-recovery stroke conformation, where it remains for ～40 s, until ATP hydrolysis products are released. Angle distributions indicate that the post-recovery stroke conformation is stabilized by ≥ 5 k(B)T of energy. The high kinetic and energetic stability of the post-recovery stroke conformation favors preferential binding of the detached head to a forward site 72 nm away. Thus, the recovery stroke contributes to unidirectional stepping of myosin Va.     

Coupling of Lever Arm Swing and Biased Brownian Motion in Actomyosin

An important unresolved problem associated with actomyosin motors is the role of Brownian motion in the process of force generation. On the basis of structural observations of myosins and actins, the widely held lever-arm hypothesis has been proposed, in which proteins are assumed to show sequential structural changes among observed and hypothesized structures to exert mechanical force. An alternative hypothesis, the Brownian motion hypothesis, has been supported by single-molecule experiments and emphasizes more on the roles of fluctuating protein movement. In this study, we address the long-standing controversy between the lever-arm hypothesis and the Brownian motion hypothesis through in silico observations of an actomyosin system. We study a system composed of myosin II and actin filament by calculating free-energy landscapes of actin-myosin interactions using the molecular dynamics method and by simulating transitions among dynamically changing free-energy landscapes using the Monte Carlo method. The results obtained by this combined multi-scale calculation show that myosin with inorganic phosphate (P_i) and ADP weakly binds to actin and that after releasing P_i and ADP, myosin moves along the actin filament toward the strong-binding site by exhibiting the biased Brownian motion, a behavior consistent with the observed single-molecular behavior of myosin. Conformational flexibility of loops at the actin-interface of myosin and the N-terminus of actin subunit is necessary for the distinct bias in the Brownian motion. Both the 5.5–11 nm displacement due to the biased Brownian motion and the 3–5 nm displacement due to lever-arm swing contribute to the net displacement of myosin. The calculated results further suggest that the recovery stroke of the lever arm plays an important role in enhancing the displacement of myosin through multiple cycles of ATP hydrolysis, suggesting a unified movement mechanism for various members of the myosin family.     

Hydration studies on the archaeal protein Sso7d using NMR measurements and MD simulations

Background

Cytoplasmic class XI myosins are the fastest processive motors known. This class functions in high-velocity cytoplasmic streaming in various plant cells from algae to angiosperms. The velocities at which they process are ten times faster than its closest class V homologues.

Results

To provide sequence determinants and structural rationale for the molecular mechanism of this fast pace myosin, we have compared the sequences from myosin class V and XI through Evolutionary Trace (ET) analysis. The current study identifies class-specific residues of myosin XI spread over the actin binding site, ATP binding site and light chain binding neck region. Sequences for ET analysis were accumulated from six plant genomes, using literature based text search and sequence searches, followed by triple validation viz. CDD search, string-based searches and phylogenetic clustering. We have identified nine myosin XI genes in sorghum and seven in grape by sequence searches. Both the plants possess one gene product each belonging to myosin type VIII as well. During this process, we have re-defined the gene boundaries for three sorghum myosin XI genes using fgenesh program.

Conclusion

Molecular modelling and subsequent analysis of putative interactions involving these class-specific residues suggest a structural basis for the molecular mechanism behind high velocity of plant myosin XI. We propose a model of a more flexible switch I region that contributes to faster ADP release leading to high velocity movement of the algal myosin XI.     

Predicting allosteric communication in myosin via a pathway of conserved residues

We present a computational method that predicts a pathway of residues that mediate protein allosteric communication. The pathway is predicted using only a combination of distance constraints between contiguous residues and evolutionary data. We applied this analysis to find pathways of conserved residues connecting the myosin ATP binding site to the lever arm. These pathway residues may mediate the allosteric communication that couples ATP hydrolysis to the lever arm recovery stroke. Having examined pre-stroke conformations of Dictyostelium, scallop, and chicken myosin II as well as Dictyostelium myosin I, we observed a conserved pathway traversing switch II and the relay helix, which is consistent with the understood need for allosteric communication in this conformation. We also examined post-rigor and rigor conformations across several myosin species. Although initial residues of these paths are more heterogeneous, all but one of these paths traverse a consistent set of relay helix residues to reach the beginning of the lever arm. We discuss our results in the context of structural elements and reported mutational experiments, which substantiate the significance of the pre-stroke pathways. Our method provides a simple, computationally efficient means of predicting a set of residues that mediate allosteric communication. We provide a refined, downloadable application and source code (on https://simtk.org) to share this tool with the wider community (https://simtk.org/home/allopathfinder).     

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.     

Coupling between catalytic site and collective dynamics: a requirement for mechanochemical activity of enzymes

Growing evidence supports the view that enzymatic activity results from a subtle interplay between chemical kinetics and molecular motions. A systematic analysis is performed here to delineate the type and level of coupling between catalysis and conformational mechanics. The dynamics of a set of 98 enzymes representative of different EC classes are analyzed with the Gaussian network model (GNM) and compared with experimental data. In more than 70% of the examined enzymes, the global hinge centers predicted by the GNM are found to be colocalized with the catalytic sites experimentally identified. Low translational mobility (< 7%) is observed for the catalytic residues, consistent with the fine-tuned design of enzymes to achieve precise mechanochemical activities. Ligand binding sites, while closely neighboring catalytic sites, enjoy a moderate flexibility to accommodate the ligand binding. These findings could serve as additional criteria for assessing drug binding residues and could lessen the computational burden of substrate docking searches.     

Influence of modulated structural dynamics on the kinetics of alpha-chymotrypsin catalysis. Insights through chemical glycosylation, molecular dynamics and domain motion analysis

Although the chemical nature of the catalytic mechanism of the serine protease alpha-chymotrypsin (alpha-CT) is largely understood, the influence of the enzyme's structural dynamics on its catalysis remains uncertain. Here we investigate whether alpha-CT's structural dynamics directly influence the kinetics of enzyme catalysis. Chemical glycosylation [Solá RJ & Griebenow K (2006) FEBS Lett 580, 1685-1690] was used to generate a series of glycosylated alpha-CT conjugates with reduced structural dynamics, as determined from amide hydrogen/deuterium exchange kinetics (k(HX)). Determination of their catalytic behavior (K(S), k(2), and k(3)) for the hydrolysis of N-succinyl-Ala-Ala-Pro-Phe p-nitroanilide (Suc-Ala-Ala-Pro-Phe-pNA) revealed decreased kinetics for the catalytic steps (k(2) and k(3)) without affecting substrate binding (K(S)) at increasing glycosylation levels. Statistical correlation analysis between the catalytic (DeltaG( not equal)k(i)) and structurally dynamic (DeltaG(HX)) parameters determined revealed that the enzyme acylation and deacylation steps are directly influenced by the changes in protein structural dynamics. Molecular modelling of the alpha-CT glycoconjugates coupled with molecular dynamics simulations and domain motion analysis employing the Gaussian network model revealed structural insights into the relation between the protein's surface glycosylation, the resulting structural dynamic changes, and the influence of these on the enzyme's collective dynamics and catalytic residues. The experimental and theoretical results presented here not only provide fundamental insights concerning the influence of glycosylation on the protein biophysical properties but also support the hypothesis that for alpha-CT the global structural dynamics directly influence the kinetics of enzyme catalysis via mechanochemical coupling between domain motions and active site chemical groups.     

Structural basis of mechanochemical coupling in a hexameric molecular motor

The P4 protein of bacteriophage phi12 is a hexameric molecular motor closely related to superfamily 4 helicases. P4 converts chemical energy from ATP hydrolysis into mechanical work, to translocate single-stranded RNA into a viral capsid. The molecular basis of mechanochemical coupling, i.e. how small approximately 1 A changes in the ATP-binding site are amplified into nanometer scale motion along the nucleic acid, is not understood at the atomic level. Here we study in atomic detail the mechanochemical coupling using structural and biochemical analyses of P4 mutants. We show that a conserved region, consisting of superfamily 4 helicase motifs H3 and H4 and loop L2, constitutes the moving lever of the motor. The lever tip encompasses an RNA-binding site that moves along the mechanical reaction coordinate. The lever is flanked by gamma-phosphate sensors (Asn-234 and Ser-252) that report the nucleotide state of neighboring subunits and control the lever position. Insertion of an arginine finger (Arg-279) into the neighboring catalytic site is concomitant with lever movement and commences ATP hydrolysis. This ensures cooperative sequential hydrolysis that is tightly coupled to mechanical motion. Given the structural conservation, the mutated residues may play similar roles in other hexameric helicases and related molecular motors.     

Selective perturbation of the myosin recovery stroke by point mutations at the base of the lever arm affects ATP hydrolysis and phosphate release

After ATP binding the myosin head undergoes a large structural rearrangement called the recovery stroke. This transition brings catalytic residues into place to enable ATP hydrolysis, and at the same time it causes a swing of the myosin lever arm into a primed state, which is a prerequisite for the power stroke. By introducing point mutations into a subdomain interface at the base of the myosin lever arm at positions Lys(84) and Arg(704), we caused modulatory changes in the equilibrium constant of the recovery stroke, which we could accurately resolve using the fluorescence signal of single tryptophan Dictyostelium myosin II constructs. Our results shed light on a novel role of the recovery stroke: fine-tuning of this reversible equilibrium influences the functional properties of myosin through controlling the effective rates of ATP hydrolysis and phosphate release.     

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.