Protein-protein interactions in actin-myosin binding and structural effects of R405Q mutation: a molecular dynamics study

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.

Loop 2 of myosin is a force-dependent inhibitor of the rigor bond

Myosin's actin-binding loop (loop 2) carries a charge opposite to that of its binding site on actin and is thought to play an important role in ionic interactions between the two molecules during the initial binding step. However, no subsequent role has been identified for loop 2 in actin-myosin binding. We used an optical trap to measure bond formation and bond rupture between actin and rigor heavy meromyosin when loaded perpendicular to the filament axis. We studied HMM with intact or proteolytically cleaved loop 2 at low and physiologic ionic strength. Here we show that the presence of intact loop 2 allows actomyosin bonds to form quickly and that they do so in a short-lived bound state. Increasing tensile load causes the transition to a long-lived state-the distinguishing behavior of a catch bond. When loop 2 was cleaved catch bond behavior was abrogated leaving only a long-lived state. These data suggest that in addition to its role in locating binding sites on actin, loop 2 is also a force-dependent inhibitor of the long-lived actomyosin complex. This may be important for reducing the duty ratio and increasing the shortening velocity of actomyosin at low forces.

All-atom molecular dynamics simulations of actin-myosin interactions: a comparative study of cardiac α myosin, β myosin, and fast skeletal muscle myosin

Myosins are a superfamily of actin-binding motor proteins with significant variations in kinetic properties (such as actin binding affinity) between different isoforms. It remains unknown how such kinetic variations arise from the structural and dynamic tuning of the actin-myosin interface at the amino acid residue level. To address this key issue, we have employed molecular modeling and simulations to investigate, with atomistic details, the isoform dependence of actin-myosin interactions in the rigor state. By combining electron microscopy-based docking with homology modeling, we have constructed three all-atom models for human cardiac α and β and rabbit fast skeletal muscle myosin in complex with three actin subunits in the rigor state. Starting from these models, we have performed extensive all-atom molecular dynamics (MD) simulations (total of 100 ns per system) and then used the MD trajectories to calculate actin-myosin binding free energies with contributions from both electrostatic and nonpolar forces. Our binding calculations are in good agreement with the experimental finding of isoform-dependent differences in actin binding affinity between these myosin isoforms. Such differences are traced to changes in actin-myosin electrostatic interactions (i.e., hydrogen bonds and salt bridges) that are highly dynamic and involve several flexible actin-binding loops. By partitioning the actin-myosin binding free energy to individual myosin residues, we have also identified key myosin residues involved in the actin-myosin interactions, some of which were previously validated experimentally or implicated in cardiomyopathy mutations, and the rest make promising targets for future mutational experiments.

Automating the search of molecular motor templates by evolutionary methods

Biological molecular motors are nanoscale devices capable of transforming chemical energy into mechanical work, which are being researched in many scientific disciplines. From a computational point of view, the characteristics and dynamics of these motors are studied at multiple time scales, ranging from very detailed and complex molecular dynamics simulations spanning a few microseconds, to extremely simple and coarse-grained theoretical models of their working cycles. However, this research is performed only in the (relatively few) instances known from molecular biology. In this work, results from elastic network analysis and behaviour-finding methods are applied to explore a subset of the configuration space of template molecular structures that are able to transform chemical energy into directed movement, for a fixed instance of working cycle. While using methods based on elastic networks limits the scope of our results, it enables the implementation of computationally lightweight methods, in a way that evolutionary search techniques can be applied to discover novel molecular motor templates. The results show that molecular motion can be attained from a variety of structural configurations, when a functional working cycle is provided. Additionally, these methods enable a new computational way to test hypotheses about molecular motors.

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.

Entropic potential field formed for a linear-motor protein near a filament: Statistical-mechanical analyses using simple models

We report a new progress in elucidating the mechanism of the unidirectional movement of a linear-motor protein (e.g., myosin) along a filament (e.g., F-actin). The basic concept emphasized here is that a potential field is entropically formed for the protein on the filament immersed in solvent due to the effect of the translational displacement of solvent molecules. The entropic potential field is strongly dependent on geometric features of the protein and the filament, their overall shapes as well as details of the polyatomic structures. The features and the corresponding field are judiciously adjusted by the binding of adenosine triphosphate (ATP) to the protein, hydrolysis of ATP into adenosine diphosphate (ADP)+Pi, and release of Pi and ADP. As the first step, we propose the following physical picture: The potential field formed along the filament for the protein without the binding of ATP or ADP+Pi to it is largely different from that for the protein with the binding, and the directed movement is realized by repeated switches from one of the fields to the other. To illustrate the picture, we analyze the spatial distribution of the entropic potential between a large solute and a large body using the three-dimensional integral equation theory. The solute is modeled as a large hard sphere. Two model filaments are considered as the body: model 1 is a set of one-dimensionally connected large hard spheres and model 2 is a double helical structure formed by two sets of connected large hard spheres. The solute and the filament are immersed in small hard spheres forming the solvent. The major findings are as follows. The solute is strongly confined within a narrow space in contact with the filament. Within the space there are locations with sharply deep local potential minima along the filament, and the distance between two adjacent locations is equal to the diameter of the large spheres constituting the filament. The potential minima form a ringlike domain in model 1 while they form a pointlike one in model 2. We then examine the effects of geometric features of the solute on the amplitudes and asymmetry of the entropic potential field acting on the solute along the filament. A large aspherical solute with a cleft near the solute-filament interface, which mimics the myosin motor domain, is considered in the examination. Thus, the two fields in our physical picture described above are qualitatively reproduced. The factors to be taken into account in further studies are also discussed.

The actin-myosin interface

In order to understand the mechanism of muscle contraction at the atomic level, it is necessary to understand how myosin binds to actin in a reversible way. We have used a novel molecular dynamics technique constrained by an EM map of the actin-myosin complex at 13-A resolution to obtain an atomic model of the strong-binding (rigor) actin-myosin interface. The constraining force resulting from the EM map during the molecular dynamics simulation was sufficient to convert the myosin head from the initial weak-binding state to the strong-binding (rigor) state. Our actin-myosin model suggests extensive contacts between actin and the myosin head (S1). S1 binds to two actin monomers. The contact surface between actin and S1 has increased dramatically compared with previous models. A number of loops in S1 and actin are involved in establishing the interface. Our model also suggests how the loop carrying the critical Arg 405 Glu mutation in S1 found in a familial cardiomyopathy might be functionally involved.

Multiscale modeling of structural dynamics underlying force generation and product release in actomyosin complex

To decrypt the mechanistic basis of myosin motor function, it is essential to probe the conformational changes in actomyosin with high spatial and temporal resolutions. In a computational effort to meet this challenge, we have performed a multiscale modeling of the allosteric couplings and transition pathway of actomyosin complex by combining coarse-grained modeling of the entire complex with all-atom molecular dynamics simulations of the active site. Our modeling of allosteric couplings at the pre-powerstroke state has pinpointed key actin-activated couplings to distant myosin parts which are critical to force generation and the sequential release of phosphate and ADP. At the post-powerstroke state, we have identified isoform-dependent couplings which underlie the reciprocal coupling between actin binding and nucleotide binding in fast Myosin II, and load-dependent ADP release in Myosin V. Our modeling of transition pathway during powerstroke has outlined a clear sequence of structural events triggered by actin binding, which lead to subsequent force generation, twisting of central beta-sheet, and the sequential release of phosphate and ADP. Finally we have performed atomistic simulations of active-site dynamics based on an on-path "transition-state" myosin conformation, which has revealed significantly weakened coordination of phosphate by Switch II, and a disrupted key salt bridge between Switch I and II. Meanwhile, the coordination of MgADP by Switch I and P loop is less perturbed. As a result, the phosphate can be released prior to MgADP. This study has shed new lights on the controversy over the structural mechanism of actin-activated phosphate release and force generation in myosin motor.

Myosin individualized: single nucleotide polymorphisms in energy transduction

BACKGROUND

Myosin performs ATP free energy transduction into mechanical work in the motor domain of the myosin heavy chain (MHC). Energy transduction is the definitive systemic feature of the myosin motor performed by coordinating in a time ordered sequence: ATP hydrolysis at the active site, actin affinity modulation at the actin binding site, and the lever-arm rotation of the power stroke. These functions are carried out by several conserved sub-domains within the motor domain. Single nucleotide polymorphisms (SNPs) affect the MHC sequence of many isoforms expressed in striated muscle, smooth muscle, and non-muscle tissue. The purpose of this work is to provide a rationale for using SNPs as a functional genomics tool to investigate structurefunction relationships in myosin. In particular, to discover SNP distribution over the conserved sub-domains and surmise what it implies about sub-domain stability and criticality in the energy transduction mechanism.

RESULTS

An automated routine identifying human nonsynonymous SNP amino acid missense substitutions for any MHC gene mined the NCBI SNP data base. The routine tested 22 MHC genes coding muscle and non-muscle isoforms and identified 89 missense mutation positions in the motor domain with 10 already implicated in heart disease and another 8 lacking sequence homology with a skeletal MHC isoform for which a crystallographic model is available. The remaining 71 SNP substitutions were found to be distributed over MHC with 22 falling outside identified functional sub-domains and 49 in or very near to myosin sub-domains assigned specific crucial functions in energy transduction. The latter includes the active site, the actin binding site, the rigid lever-arm, and regions facilitating their communication. Most MHC isoforms contained SNPs somewhere in the motor domain.

CONCLUSIONS

Several functional-crucial sub-domains are infiltrated by a large number of SNP substitution sites suggesting these domains are engineered by evolution to be too-robust to be disturbed by otherwise intrusive sequence changes. Two functional sub-domains are SNP-free or relatively SNP-deficient but contain many disease implicated mutants. These sub-domains are apparently highly sensitive to any missense substitution suggesting they have failed to evolve a robust sequence paradigm for performing their function.

The myosin C-loop is an allosteric actin contact sensor in actomyosin

Actin and myosin form the molecular motor in muscle. Myosin is the enzyme performing ATP hydrolysis under the allosteric control of actin such that actin binding initiates product release and force generation in the myosin power stroke. Non-equilibrium Monte Carlo molecular dynamics simulation of the power stroke suggested that a structured surface loop on myosin, the C-loop, is the actin contact sensor initiating actin activation of the myosin ATPase. Previous experimental work demonstrated C-loop binds actin and established the forward and reverse allosteric link between the C-loop and the myosin active site. Here, smooth muscle heavy meromyosin C-loop chimeras were constructed with skeletal (sCl) and cardiac (cCl) myosin C-loops substituted for the native sequence. In both cases, actin-activated ATPase inhibition is indicated mainly by the lower V(max). In vitro motility was also inhibited in the chimeras. Motility data were collected as a function of myosin surface density, with unregulated actin, and with skeletal and cardiac isoforms of tropomyosin-bound actin for the wild type, cCl, and sCl. Slow and fast subpopulations of myosin velocities in the wild-type species were discovered and represent geometrically unfavorable and favorable actomyosin interactions, respectively. Unfavorable interactions are detected at all surface densities tested. Favorable interactions are more probable at higher myosin surface densities. Cardiac tropomyosin-bound actin promotes the favorable actomyosin interactions by lowering the inhibiting geometrical constraint barriers with a structural effect on actin. Neither higher surface density nor cardiac tropomyosin-bound actin can accelerate motility velocity in cCl or sCl, suggesting the element initiating maximal myosin activation by actin resides in the C-loop.

Dynamic conformational changes due to the ATP hydrolysis in the motor domain of myosin: 10-ns molecular dynamics simulations

Muscle contraction is caused by directed movement of myosin heads along actin filaments. This movement is triggered by ATP hydrolysis, which occurs within the motor domain of myosin. The mechanism for this intramolecular process remains unknown owing to a lack of ways to observe the detailed motions of each atom in the myosin molecule. We carried out 10-ns all-atom molecular dynamics simulations to investigate the types of dynamic conformational changes produced in the motor domain by the energy released from ATP hydrolysis. The results revealed that the thermal fluctuations modulated by perturbation of ATP hydrolysis are biased in one direction that is relevant to directed movement of the myosin head along the actin filament.

Actomyosin interaction: mechanical and energetic properties in different nucleotide binding states

The mechanics of the actomyosin interaction is central in muscle contraction and intracellular trafficking. A better understanding of the events occurring in the actomyosin complex requires the examination of all nucleotide-dependent states and of the energetic features associated with the dynamics of the cross-bridge cycle. The aim of the present study is to estimate the interaction strength between myosin in nucleotide-free, ATP, ADP·Pi and ADP states and actin monomer. The molecular models of the complexes were constructed based on cryo-electron microscopy maps and the interaction properties were estimated by means of a molecular dynamics approach, which simulate the unbinding of the complex applying a virtual spring to the core of myosin protein. Our results suggest that during an ATP hydrolysis cycle the affinity of myosin for actin is modulated by the presence and nature of the nucleotide in the active site of the myosin motor domain. When performing unbinding simulations with a pulling rate of 0.001 nm/ps, the maximum pulling force applied to the myosin during the experiment is about 1nN. Under these conditions the interaction force between myosin and actin monomer decreases from 0.83 nN in the nucleotide-free state to 0.27 nN in the ATP state, and increases to 0.60 nN after ATP hydrolysis and Pi release from the complex (ADP state).

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.

Molecular modeling of the cytoskeleton

Molecular modeling techniques have truly come of age in recent decades, and here we cover several of the most commonly used techniques, namely molecular dynamics, Brownian dynamics, and molecular docking. In each case, we explain the physical basis and limitations of the various techniques and then illustrate their application to various problems related to the cytoskeleton. This set of studies covers a relatively wide range of examples and is comprehensive enough to clearly see how these techniques could be applied to other systems. Finally, we cover several related methodologies that expand on these basic techniques to allow for more detailed and specific simulation and analysis.

Conformational changes in protein function

Conformational changes are the hallmarks of protein dynamics and are often intimately related to protein functions. Molecular dynamics (MD) simulation is a powerful tool to study the time-resolved properties of protein structure in atomic details. In this chapter, we discuss the various applications of MD simulation to the study of protein conformational changes, and introduce several selected advanced techniques that may significantly increase the sampling efficiencies, including locally enhanced sampling (LES), and grow-to-fit molecular dynamics (G2FMD).

Myosin dynamics on the millisecond time scale

Myosin is a motor protein associating with actin and ATP. It translates along actin filaments against a force by transduction of free energy liberated with ATP hydrolysis. Various myosin crystal structures define time points during ATPase showing the protein undergoes large conformation change during transduction over a cycle with approximately 10 ms periodicity. The protein conformation trajectory between two intermediates in the cycle is surmised by non-equilibrium Monte Carlo simulation utilizing free-energy minimization. The trajectory shows myosin transduction of free energy to mechanical work giving evidence for: (i) a causal relationship between product release and work production in the native isoform that is correctly disrupted in a chemically modified protein, (ii) the molecular basis of ATP-sensitive tryptophan fluorescence enhancement and acrylamide quenching, (iii) an actin-binding site peptide containing the free-energy barrier to ATPase product release defining the rate limiting step and, (iv) a scenario for actin-activation of myosin ATPase.

Understanding antibody–antigen associations by molecular dynamics simulations: Detection of important intra- and inter-molecular salt bridges

1 NSec molecular dynamics (MD) simulation of anti-hen egg white antibody, HyHEL63 (HH63), complexed with HEL reveals important molecular interactions, not revealed in its X-ray crystal structure. These molecular interactions were predicted to be critical for the complex formation, based on structure-function studies of this complex and 3-other anti-HEL antibodies, HH8, HH10 and HH26, HEL complexes. All four antibodies belong to the same structural family, referred to here as HH10 family. Ala scanning results show that they recognize 'coincident epitopes'. 1 NSec explicit, with periodic boundary condition, MD simulation of HH63- HEL reveals the presence of functionally important saltbridges. Around 200 ps in vacuo and an additional 20 ps explicit simulation agree with the observations from 1 Nsec simulation. Intra-molecular salt-bridges predicted to play significant roles in the complex formation, were revealed during MD simulation. A very stabilizing saltbridge network, and another intra-molecular salt-bridge, at the binding site of HEL, revealed during the MD simulation, is proposed to predipose binding site geometry for specific binding. All the revealed saltbridges are present in one or more of the other three complexes and/or involve \"hot-spot\" epitope and paratope residues. Most of these charged epitope residues make large contribution to the binding free energy. The "hot spot" epitope residue Lys97Y, which significantly contributes to the free energy of binding in all the complexes, forms an intermolecular salt-bridge in several MD conformers. Our earlier computations have shown that this inter-molecular salt-bridge plays a significant role in determining specificity and flexibility of binding in the HH8-HEL and HH26-HEL complexes. Using a robust criterion of salt-bridge detection, this intermolecular salt-bridge was detected in the native structures of the HH8-HEL and HH26-HEL complexes, but was not revealed in the crystal structure of HH63-HEL complex. The electrostatic strength of this revealed saltbridge was very strong. During 1 Nsec MD simulation this salt-bridge networks with another inter-molecular salt-bridge to form an inter-molecular salt-bridge triad. Participation of Lys97Y in the formation of inter-molecular triad further validates the functional importance of Lys97Y in HH63-HEL associations. These results demonstrate that many important structural details of biomolecular interactions can be better understood when studied in a dynamic environment, and that MD simulations can complement and expand information obtained from static X-ray structure. This study also highlights "hot-spot" molecular interactions in HyHEL63-HEL complex.

The structural determinants of macrolide-actin binding: in silico insights

By the use of x-ray structures and flexible docking, we have developed the first in silico ligand-based view of the structural determinants of the binding of small molecule mimics of gelsolin, natural products bound to actin. Our technique highlights those residues on the actin binding site forming important hydrophobic and hydrogen-bonding interactions with the ligands. Significantly, through the flexible docking of toxin fragments, we have also identified potential residues on the actin binding site that have yet to be exploited. Guided by these observations, we have demonstrated that kabiramide C can be modified to produce a structure with a predicted binding energy increased by 20% while the molecular mass is reduced by 20%, clearly indicating the potential for future elaboration of structures targeting this important component of the cytoskeleton.

Myosin Surface Loop 4 Modulates Inhibition of Actomyosin 1b ATPase Activity by Tropomyosin

Structural studies of the class I myosin, MyoE, led to the predictions that loop 4, a surface loop near the actin-binding region that is longer in class I myosins than in other myosin subclasses, might limit binding of myosins I to actin when actin-binding proteins, like tropomyosin, are present, and might account for the exclusion of myosin I from stress fibers. To test these hypotheses, mutant molecules of the related mammalian class I myosin, Myo1b, in which loop 4 was truncated (from an amino acid sequence of RMNGLDES to NGLD) or replaced with the shorter and distinct loop 4 found in Dictyostelium myosin II (GAGEGA), were expressed in vitro and their interaction with actin and with actin-tropomyosin was tested. Saturating amounts of expressed fibroblast tropomyosin-2 resulted in a decrease in the maximum actin-activated Mg2+-ATPase activity of wild-type Myo1b but had little or no effect on the actin-activated Mg2+-ATPase activity of the two mutants. In motility assays, few actin filaments bound tightly to Myo1b-WT-coated cover slips when tropomyosin-2 was present, whereas actin filaments both bound and were translocated by Myo1b-NGLD or Myo1b-GAGEGA in both the presence and absence of tropomyosin-2. When expressed in mammalian cells, like the wild type, the mutant myosins were largely excluded from tropomyosin-containing actin filaments, indicating that in the cell additional factors besides loop 4 determine targeting of myosins I to specific subpopulations of actin filaments.

Conformational sampling with implicit solvent models: application to the PHF6 peptide in tau protein

Implicit solvent models approximate the effects of solvent through a potential of mean force and therefore make solvated simulations computationally efficient. Yet despite their computational efficiency, the inherent approximations made by implicit solvent models can sometimes lead to inaccurate results. To test the accuracy of a number of popular implicit solvent models, we determined whether implicit solvent simulations can reproduce the set of potential energy minima obtained from explicit solvent simulations. For these studies, we focus on a six-residue amino-acid sequence, referred to as the paired helical filament 6 (PHF6), which may play an important role in the formation of intracellular aggregates in patients with Alzheimer's disease. Several implicit solvent models form the basis of this work--two based on the generalized Born formalism, and one based on a Gaussian solvent-exclusion model. All three implicit solvent models generate minima that are in good agreement with minima obtained from simulations with explicit solvent. Moreover, free-energy profiles generated with each implicit solvent model agree with free-energy profiles obtained with explicit solvent. For the Gaussian solvent-exclusion model, we demonstrate that a straightforward ranking of the relative stability of each minimum suggests that the most stable structure is extended, a result in excellent agreement with the free-energy profiles. Overall, our data demonstrate that for some peptides like PHF6, implicit solvent can accurately reproduce the set of local energy minimum arising from quenched dynamics simulations with explicit solvent. More importantly, all solvent models predict that PHF6 forms extended beta-structures in solution, a finding consistent with the notion that PHF6 initiates neurofibrillary tangle formation in patients with Alzheimer's disease.

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.