Energetics of subdomain movements and fluorescence probe solvation environment change in ATP-bound myosin

Energetics of conformational changes experienced by an ATP-bound myosin head detached from actin was studied by all-atom explicit water umbrella sampling simulations. The statistics of coupling between large scale domain movements and smaller scale structural features were examined, including the closing of the ATP binding pocket, and a number of key hydrogen bond formations shown to play roles in structural and biochemical studies. The statistics for the ATP binding pocket open/close transition show an evolution of the relative stability from the open state in the early stages of the recovery stroke to the stable closed state after the stroke. The change in solvation environment of the fluorescence probe Trp507 (scallop numbering; 501 in Dictyostelium discoideum) indicates that the probe faithfully reflects the closing of the binding pocket as previously shown in experimental studies, while being directly coupled to roughly the early half of the overall large scale conformational change of the converter domain rotation. The free energy change of this solvation environment change, in particular, is −1.3 kcal/mol, in close agreement with experimental estimates. In addition, our results provide direct molecular level data allowing for interpretations of the fluorescence experiments of myosin conformational change in terms of the de-solvation of Trp side chain.     

Myosin flexibility: structural domains and collective vibrations

The movement of the myosin motor along an actin filament involves a directed conformational change within the cross-bridge formed between the protein and the filament. Despite the structural data that has been obtained on this system, little is known of the mechanics of this conformational change. We have used existing crystallographic structures of three conformations of the myosin head, containing the motor domain and the lever arm, for structural comparisons and mechanical studies with a coarse-grained elastic network model. The results enable us to define structurally conserved domains within the protein and to better understand myosin flexibility. Notably they point to the role of the light chains in rigidifying the lever arm and to changes in flexibility as a consequence of nucleotide binding.     

Conformational selection during weak binding at the actin and myosin interface

The molecular mechanism of the powerstroke in muscle is examined by resonance energy transfer techniques. Recent models suggesting a pre-cocking of the myosin head involving an enormous rotation between the lever arm and the catalytic domain were tested by measuring separation distances among myosin subfragment-2, the nucleotide site, and the regulatory light chain in the presence of nucleotide transition state analogs. Only small changes (<0.5 nm) were detected that are consistent with internal conformational changes of the myosin molecule, but not with extreme differences in the average lever arm position suggested by some atomic models. These results were confirmed by stopped-flow resonance energy transfer measurements during single ATP turnovers on myosin. To examine the participation of actin in the powerstroke process, resonance energy transfer between the regulatory light chain on myosin subfragment-1 and the C-terminus of actin was measured in the presence of nucleotide transition state analogs. The efficiency of energy transfer was much greater in the presence of ADP-AlF(4), ADP-BeF(x), and ADP-vanadate than in the presence of ADP or no nucleotide. These data detect profound differences in the conformations of the weakly and strongly attached cross-bridges that appear to result from a conformational selection that occurs during the weak binding of the myosin head to actin.     

Crystal structure of scallop Myosin s1 in the pre-power stroke state to 2.6 a resolution: flexibility and function in the head

We have extended the X-ray structure determination of the complete scallop myosin head in the pre-power stroke state to 2.6 A resolution, allowing an atomic comparison of the three major (weak actin binding) states of various myosins. We can now account for conformational differences observed in crystal structures in the so-called "pliant region" at the motor domain-lever arm junction between scallop and vertebrate smooth muscle myosins. A hinge, which may contribute to the compliance of the myosin crossbridge, has also been identified for the first time within the regulatory light-chain domain of the lever arm. Analysis of temperature factors of key joints of the motor domain, especially the SH1 helix, provides crystallographic evidence for the existence of the "internally uncoupled" state in diverse isoforms. The agreement between structural and solution studies reinforces the view that the unwinding of the SH1 helix is a part of the cross-bridge cycle in many myosins.     

How myosin VI coordinates its heads during processive movement

A processive molecular motor must coordinate the enzymatic state of its two catalytic domains in order to prevent premature detachment from its track. For myosin V, internal strain produced when both heads of are attached to an actin track prevents completion of the lever arm swing of the lead head and blocks ADP release. However, this mechanism cannot work for myosin VI, since its lever arm positions are reversed. Here, we demonstrate that myosin VI gating is achieved instead by blocking ATP binding to the lead head once it has released its ADP. The structural basis for this unique gating mechanism involves an insert near the nucleotide binding pocket that is found only in class VI myosin. Reverse strain greatly favors binding of ADP to the lead head, which makes it possible for myosin VI to function as a processive transporter as well as an actin-based anchor. While this mechanism is unlike that of any other myosin superfamily member, it bears remarkable similarities to that of another processive motor from a different superfamily--kinesin I.     

Mutation of a conserved glycine in the SH1-SH2 helix affects the load-dependent kinetics of myosin

The ATP hydrolysis rate and shortening velocity of muscle are load-dependent. At the molecular level, myosin generates force and motion by coupling ATP hydrolysis to lever arm rotation. When a laser trap was used to apply load to single heads of expressed smooth muscle myosin (S1), the ADP release kinetics accelerated with an assistive load and slowed with a resistive load; however, ATP binding was mostly unaffected. To investigate how load is communicated within the motor, a glycine located at the putative fulcrum of the lever arm was mutated to valine (G709V). In the absence of load, stopped-flow and laser trap studies showed that the mutation significantly slowed the rates of ADP release and ATP binding, accounting for the ~270-fold decrease in actin sliding velocity. The load dependence of the mutant's ADP release rate was the same as that of wild-type S1 (WT) despite the slower rate. In contrast, load accelerated ATP binding by ~20-fold, irrespective of loading direction. Imparting mechanical energy to the mutant motor partially reversed the slowed ATP binding by overcoming the elevated activation energy barrier. These results imply that conformational changes near the conserved G709 are critical for the transmission of mechanochemical information between myosin's active site and lever arm.     

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).     

Myosin motor domain lever arm rotation is coupled to ATP hydrolysis

We have investigated coupling of lever arm rotation to the ATP binding and hydrolysis steps for the myosin motor domain. In several current hypotheses of the mechanism of force production by muscle, the primary mechanical feature is the rotation of a lever arm that is a subdomain of the myosin motor domain. In these models, the lever arm rotates while the myosin motor domain is free, and then reverses the rotation to produce force while it is bound to actin. These mechanical steps are coupled to steps in the ATP hydrolysis cycle. Our hypothesis is that ATP hydrolysis induces lever arm rotation to produce a more compact motor domain that has stored mechanical energy. Our approach is to use transient electric birefringence techniques to measure changes in hydrodynamic size that result from lever arm rotation when various ligands are bound to isolated skeletal muscle myosin motor domain in solution. Results for ATP and CTP, which do support force production by muscle fibers, are compared to those of ATPgammaS and GTP, which do not. Measurements are also made of conformational changes when the motor domain is bound to NDP's and PP(i) in the absence and presence of the phosphate analogue orthovanadate, to determine the roles the nucleoside moieties of the nucleotides have on lever arm rotation. The results indicate that for the substrates investigated, rotation does not occur upon substrate binding, but is coupled to the NTP hydrolysis step. The data are consistent with a model in which only substrates that produce a motor domain-NDP-P(i) complex as the steady-state intermediate make the motor domain more compact, and only those substrates support force production.     

Simulations of the myosin II motor reveal a nucleotide-state sensing element that controls the recovery stroke

During the recovery stroke, the myosin motor is primed for the next power stroke by a 60° rotation of its lever arm. This reversible motion is coupled to the activation of the ATPase function of myosin through conformational changes along the relay helix, which runs from the Switch-2 loop near the ATP to the converter domain carrying the lever arm. Via a hydrogen bond between the side-chain of Asn475 on the relay helix and the Gly457/Ser456 peptide group on the Switch-2, the rotation of the converter domain is coupled to the formation of a hydrogen bond between Gly457 and γ-phosphate that is essential for ATP hydrolysis. Here, molecular dynamics simulations of Dictyostelium discoideum myosin II in the two end conformations of the recovery stroke with different nucleotide states (ATP, ADP·Pi, ADP) reveal that the side-chain of Asn475 breaks away from Switch-2 upon ATP hydrolysis to make a hydrogen bond with Tyr573. This sensing of the nucleotide state is achieved by a small displacement of the cleaved γ-phosphate towards Gly457 which in turn pushes Asn475 away. The sensing plays a dual role by (i) preventing the wasteful reversal of the recovery stroke while the nucleotide is in the ADP·Pi state, and (ii) decoupling the relay helix from Switch-2, thus allowing the power stroke to start upon initial binding to actin while Gly457 of Switch-2 keeps interacting with the Pi (known to be released only later after tight actin binding). A catalytically important salt bridge between Arg238 (on Switch-1) and Glu459 (on Switch-2), which covers the hydrolysis site, is seen to form rapidly when ATP is added to the pre-recovery stroke conformer and remains stable after the recovery stroke, indicating that it has a role in shaping the ATP binding site by induced fit.     

Kinetics of structural changes in the relay loop and SH3 domain of myosin

The intrinsic fluorescence of smooth muscle myosin signals conformational changes associated with different catalytic states of the ATPase cycle. To elucidate this relationship, we have examined the pre-steady-state kinetics of nucleotide binding, hydrolysis, and product release in motor domain-essential light chain mutants containing a single endogenous tryptophan, either residue 512 in the rigid relay loop or residue 29 adjacent to the SH3 domain. The intrinsic fluorescence of W512 is sensitive to both nucleotide binding and hydrolysis, and appears to report structural changes at the active site, presumably through a direct connection with switch II. The intrinsic fluorescence of W29 is sensitive to nucleotide binding but not hydrolysis, and does not appear to be tightly linked with structural changes occurring at the active site. We propose that the SH3 domain may be sensitive to conformational changes in the lever arm through contacts with the essential light chain.     

The structural basis of muscle contraction

The myosin cross-bridge exists in two conformations, which differ in the orientation of a long lever arm. Since the lever arm undergoes a 60 degree rotation between the two conformations, which would lead to a displacement of the myosin filament of about 11 nm, the transition between these two states has been associated with the elementary 'power stroke' of muscle. Moreover, this rotation is coupled with changes in the active site (CLOSED to OPEN), which probably enable phosphate release. The transition CLOSED to OPEN appears to be brought about by actin binding. However, kinetics shows that the binding of myosin to actin is a two-step process which affects both ATP and ADP affinity and vice versa. The structural basis of these effects is only partially explained by the presently known conformers of myosin. Therefore, additional states of the myosin cross-bridge should exist. Indeed, cryoelectron microscopy has revealed other angles of the lever arm induced by ADP binding to a smooth muscle actin-myosin complex.     

Detailed tuning of structure and intramolecular communication are dispensable for processive motion of myosin VI

Dimeric myosin VI moves processively hand-over-hand along actin filaments. We have characterized the mechanism of this processive motion by measuring the impact of structural and chemical perturbations on single-molecule processivity. Processivity is maintained despite major alterations in lever arm structure, including replacement of light chain binding regions and elimination of the medial tail. We present kinetic models that can explain the ATP concentration-dependent processivities of myosin VI constructs containing either native or artificial lever arms. We conclude that detailed tuning of structure and intramolecular communication are dispensable for processive motion, and further show theoretically that one proposed type of nucleotide gating can be detrimental rather than beneficial for myosin processivity.     

Switch I closure simultaneously promotes strong binding to actin and ADP in smooth muscle myosin

The motor protein myosin uses energy derived from ATP hydrolysis to produce force and motion. Important conserved components (P-loop, switch I, and switch II) help propagate small conformational changes at the active site into large scale conformational changes in distal regions of the protein. Structural and biochemical studies have indicated that switch I may be directly responsible for the reciprocal opening and closing of the actin and nucleotide-binding pockets during the ATPase cycle, thereby aiding in the coordination of these important substrate-binding sites. Smooth muscle myosin has displayed the ability to simultaneously bind tightly to both actin and ADP, although it is unclear how both substrate-binding clefts could be closed if they are rigidly coupled to switch I. Here we use single tryptophan mutants of smooth muscle myosin to determine how conformational changes in switch I are correlated with structural changes in the nucleotide and actin-binding clefts in the presence of actin and ADP. Our results suggest that a closed switch I conformation in the strongly bound actomyosin-ADP complex is responsible for maintaining tight nucleotide binding despite an open nucleotide-binding pocket. This unique state is likely to be crucial for prolonged tension maintenance in smooth muscle.     

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.     

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).     

Experimental investigation of the seesaw mechanism of the relay region that moves the myosin lever arm

A seesaw-like movement of the relay region upon the recovery step of myosin was recently simulated in silico. In this model the relay helix tilts around its pivoting point formed by a phenylalanine cluster (Phe(481), Phe(482), and Phe(652)), which moves the lever arm of myosin. To study the effect of the elimination of the proposed pivoting point, these phenylalanines were mutated to alanines in two Dictyostelium myosin II motor domain constructs (M(F481A, F482A) and M(F652A)). The relay movement was followed by the fluorescence change of Trp(501) located in the relay region. The steady-state and transient kinetic fluorescence experiments showed that the lack of the phenylalanine fulcrum perturbs the formation of the "up" lever arm state, and only moderate effects were found in the nucleotide binding, the formation of the "down" lever arm position, and the ATP hydrolysis steps. We conclude that the lack of the fulcrum decouples the distal part of the relay from the nucleotide binding site upon the recovery step. Our molecular dynamics simulations also showed that the conformation of the motor is not perturbed by the mutation in the down lever arm state, however, the lack of the pivoting point rearranges the dynamic pattern of the kink region of the relay helix.     

Conformational sampling and nucleotide-dependent transitions of the GroEL subunit probed by unbiased molecular dynamics simulations

GroEL is an ATP dependent molecular chaperone that promotes the folding of a large number of substrate proteins in E. coli. Large-scale conformational transitions occurring during the reaction cycle have been characterized from extensive crystallographic studies. However, the link between the observed conformations and the mechanisms involved in the allosteric response to ATP and the nucleotide-driven reaction cycle are not completely established. Here we describe extensive (in total long) unbiased molecular dynamics (MD) simulations that probe the response of GroEL subunits to ATP binding. We observe nucleotide dependent conformational transitions, and show with multiple 100 ns long simulations that the ligand-induced shift in the conformational populations are intrinsically coded in the structure-dynamics relationship of the protein subunit. Thus, these simulations reveal a stabilization of the equatorial domain upon nucleotide binding and a concomitant "opening" of the subunit, which reaches a conformation close to that observed in the crystal structure of the subunits within the ADP-bound oligomer. Moreover, we identify changes in a set of unique intrasubunit interactions potentially important for the conformational transition.     

The effect of F-actin on the relay helix position of myosin II, as revealed by tryptophan fluorescence, and its implications for mechanochemical coupling

The fluorescence properties of Dictyostelium discoideum (Dd) myosin II constructs containing a single tryptophan residue have revealed detailed information regarding nucleotide binding and hydrolysis steps. Here we extend these studies to investigate the influence of actin on nucleotide-induced fluorescence transients. The fluorescence from native actin tryptophan residues is not significantly perturbed on binding to myosin, although an apparent signal is detected as a consequence of a light scatter artifact. Actin has a minor effect on the response of W129, located at the entrance to the nucleotide-binding pocket, and reduces the forward rate constants for the isomerization(s) associated with binding of ATP, ATPgammaS, and ADP by 3-fold or less. The isomerization detected by W129 clearly precedes the dissociation of actin in the case of ADP and ATPgammaS binding. The fluorescence from the conserved W501 residue, located at the distal end of the relay helix, is very sensitive to the switch 2 and/or lever arm disposition. Consequently, the observed fluorescence emission intensity can be used to estimate the equilibrium constant between the pre- and post-power stroke conformations. Actin modulates this equilibrium by no more than 2-fold in the presence of nucleoside triphosphate. These data have implications for the mechanism of product release and suggest that actin activates another process in the mechanism, such as switch 1 movement and Pi release, rather than influencing the switch 2 equilibrium and lever arm position directly.     

Three myosin V structures delineate essential features of chemo-mechanical transduction

The molecular motor, myosin, undergoes conformational changes in order to convert chemical energy into force production. Based on kinetic and structural considerations, we assert that three crystal forms of the myosin V motor delineate the conformational changes that myosin motors undergo upon detachment from actin. First, a motor domain structure demonstrates that nucleotide-free myosin V adopts a specific state (rigor-like) that is not influenced by crystal packing. A second structure reveals an actomyosin state that favors rapid release of ADP, and differs from the rigor-like state by a P-loop rearrangement. Comparison of these structures with a third structure, a 2.0 angstroms resolution structure of the motor bound to an ATP analog, illuminates the structural features that provide communication between the actin interface and nucleotide-binding site. Paramount among these is a region we name the transducer, which is composed of the seven-stranded beta-sheet and associated loops and linkers. Reminiscent of the beta-sheet distortion of the F1-ATPase, sequential distortion of this transducer region likely controls sequential release of products from the nucleotide pocket during force generation.     

Glycine 699 is pivotal for the motor activity of skeletal muscle myosin

Myosin couples ATP hydrolysis to the translocation of actin filaments to power many forms of cellular motility. A striking feature of the structure of the muscle myosin head domain is a 9-nm long "lever arm" that has been postulated to produce a 5-10-nm power stroke. This motion must be coupled to conformational changes around the actin and nucleotide binding sites. The linkage of these sites to the lever arm has been analyzed by site-directed mutagenesis of a conserved glycine residue (G699) found in a bend joining two helices containing the highly reactive and mobile cysteine residues, SH1 and SH2. Alanine mutagenesis of this glycine (G699A) dramatically alters the motor activity of skeletal muscle myosin, inhibiting the velocity of actin filament movement by > 100-fold. Analysis of the defect in the G699A mutant myosin is consistent with a marked slowing of the transition within the motor domain from a strong binding to a weak binding interaction with actin. This result is interpreted in terms of the role of this residue (G699) as a pivot point for motion of the lever arm. The recombinant myosin used in these experiments has been produced in a unique expression system. A shuttle vector containing a regulated muscle-specific promoter has been developed for the stable expression of recombinant myosin in C2C12 cells. The vector uses the promoter/enhancer region, the first two and the last five exons of an embryonic rat myosin gene, to regulate the expression of an embryonic chicken muscle myosin cDNA. Stable cell lines transfected with this vector express the unique genetically engineered myosin after differentiation into myotubes. The myosin assembles into myofibrils, copurifies with the endogenous myosin, and contains a complement of muscle-specific myosin light chains. The functional activity of the recombinant myosin is readily analyzed with an in vitro motility assay using a species-specific anti-S2 mAb to selectively assay the recombinant protein. This expression system has facilitated manipulation and analysis of the skeletal muscle myosin motor domain and is also amenable to a wide range of structure-function experiments addressing questions unique to the muscle-specific cytoarchitecture and myosin isoforms.