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

UNC-45 assisted myosin folding depends on a conserved FX3HY motif implicated in Freeman Sheldon Syndrome

Myosin motors are critical for diverse motility functions, ranging from cytokinesis and endocytosis to muscle contraction. The UNC-45 chaperone controls myosin function mediating the folding, assembly, and degradation of the muscle protein. Here, we analyze the molecular mechanism of UNC-45 as a hub in myosin quality control. We show that UNC-45 forms discrete complexes with folded and unfolded myosin, forwarding them to downstream chaperones and E3 ligases. Structural analysis of a minimal chaperone:substrate complex reveals that UNC-45 binds to a conserved FX3HY motif in the myosin motor domain. Disrupting the observed interface by mutagenesis prevents myosin maturation leading to protein aggregation in vivo. We also show that a mutation in the FX3HY motif linked to the Freeman Sheldon Syndrome impairs UNC-45 assisted folding, reducing the level of functional myosin. These findings demonstrate that a faulty myosin quality control is a critical yet unexplored cause of human myopathies.

A weak coupling mechanism for the early steps of the recovery stroke of myosin VI: A free energy simulation and string method analysis

Myosin motors use the energy of ATP to produce force and directed movement on actin by a swing of the lever-arm. ATP is hydrolysed during the off-actin re-priming transition termed recovery stroke. To provide an understanding of chemo-mechanical transduction by myosin, it is critical to determine how the reverse swing of the lever-arm and ATP hydrolysis are coupled. Previous studies concluded that the recovery stroke of myosin II is initiated by closure of the Switch II loop in the nucleotide-binding site. Recently, we proposed that the recovery stroke of myosin VI starts with the spontaneous re-priming of the converter domain to a putative pre-transition state (PTS) intermediate that precedes Switch II closing and ATPase activation. Here, we investigate the transition from the pre-recovery, post-rigor (PR) state to PTS in myosin VI using geometric free energy simulations and the string method. First, our calculations rediscover the PTS state agnostically and show that it is accessible from PR via a low free energy transition path. Second, separate path calculations using the string method illuminate the mechanism of the PR to PTS transition with atomic resolution. In this mechanism, the initiating event is a large movement of the converter/lever-arm region that triggers rearrangements in the Relay-SH1 region and the formation of the kink in the Relay helix with no coupling to the active site. Analysis of the free-energy barriers along the path suggests that the converter-initiated mechanism is much faster than the one initiated by Switch II closure, which supports the biological relevance of PTS as a major on-pathway intermediate of the recovery stroke in myosin VI. Our analysis suggests that lever-arm re-priming and ATP hydrolysis are only weakly coupled, so that the myosin recovery stroke is initiated by thermal fluctuations and stabilised by nucleotide consumption via a ratchet-like mechanism.

Insights into the Mechanism of the Cardiac Drug Omecamtiv Mecarbil─A Computational Study

Omecamtiv mecarbil (OM) is a positive inotrope that is thought to bind directly to an allosteric site of the β-cardiac myosin. The drug is under investigation for the treatment of systolic heart failure. The drug is classified as a cardiac myosin modulator and has been observed to affect multiple vital steps of the cross-bridge cycle including the recovery stroke and the chemical step. We explored the free-energy surface of the recovery stroke of the human cardiac β-myosin in the presence of OM to determine its influence on this process. We also investigated the effects of OM on the recovery stroke in the presence of genetic cardiomyopathic mutations R712L, F764L, and P710R using metadynamics. We also utilized the method of transition path sampling to generate an unbiased ensemble of reactive trajectories for the ATP hydrolysis step in the presence of OM that were able to provide insight into the differences observed due to OM in the dynamics and mechanism of the decomposition of ATP to ADP and HPO₄^2-, a central part of the power generation in cardiac muscle. We studied chemistry in the presence of the same three mutations to further elucidate the effect of OM, and its use in the treatment of cardiac disease.

Investigation of the Recovery Stroke and ATP Hydrolysis and Changes Caused Due to the Cardiomyopathic Point Mutations in Human Cardiac β Myosin

Human cardiac β myosin undergoes the cross-bridge cycle as part of the force-generating mechanism of cardiac muscle. The recovery stroke is considered one of the key steps of the kinetic cycle as it is the conformational rearrangement required to position the active site residues for hydrolysis of ATP and interaction with actin. We explored the free-energy surface of the transition and investigated the effect of the genetic cardiomyopathy causing mutations R453C, I457T, and I467T on this step using metadynamics. This work extends previous studies on Dictyostelium myosin II with engineered mutations. Here, like previously, we generated an unbiased thermodynamic ensemble of reactive trajectories for the chemical step using transition path sampling. Our methodologies were able to predict the changes to the dynamics of the recovery stroke as well as predict the pathway of breakdown of ATP to ADP and HPO₄^2- with the stabilization of the metaphosphate intermediate. We also observed clear differences between the Dictyostelium myosin II and human cardiac β myosin for ATP hydrolysis as well as predict the effect of the mutation I467T on the chemical step.

Mechanochemical Function of Myosin II: Investigation into the Recovery Stroke and ATP Hydrolysis

Myosin regulates muscle function through a complex cycle of conformational rearrangements coupled with the hydrolysis of adenosine triphosphate (ATP). The recovery stroke reorganizes the myosin active site to hydrolyze ATP and cross bridge with the thin filament to produce muscle contraction. Engineered mutations K84M and R704E in Dictyostelium myosin have been designed to specifically inhibit the recovery stroke and have been shown to indirectly affect the ATPase activity of myosin. We investigated these mutagenic perturbations to the recovery stroke and generated thermodynamically correct and unbiased trajectories for native ATP hydrolysis with computationally enhanced sampling methods. Our methodology was able to resolve experimentally observed changes to kinetic and equilibrium dynamics for the recovery stroke with the correct prediction in the severity of these changes. For ATP hydrolysis, the sequential nature along with the stabilization of a metaphosphate intermediate was observed in agreement with previous studies. However, we observed glutamate 459 being utilized as a proton abstractor to prime the attacking water instead of a lytic water, a phenomenon not well categorized in myosin but has in other ATPases. Both rare event methodologies can be extended to human myosin to investigate isoformic differences from Dictyostelium and scan cardiomyopathic mutations to see differential perturbations to kinetics of other conformational changes in myosin such as the power stroke.

Regulation and Plasticity of Catalysis in Enzymes: Insights from Analysis of Mechanochemical Coupling in Myosin

The mechanism of ATP hydrolysis in the myosin motor domain is analyzed using a combination of DFTB3/CHARMM simulations and enhanced sampling techniques. The motor domain is modeled in the pre-powerstroke state, in the post-rigor state, and as a hybrid based on the post-rigor state with a closed nucleotide-binding pocket. The ATP hydrolysis activity is found to depend on the positioning of nearby water molecules, and a network of polar residues facilitates proton transfer and charge redistribution during hydrolysis. Comparison of the observed hydrolysis pathways and the corresponding free energy profiles leads to detailed models for the mechanism of ATP hydrolysis in the pre-powerstroke state and proposes factors that regulate the hydrolysis activity in different conformational states. In the pre-powerstroke state, the scissile P_γ-O_3β bond breaks early in the reaction. Proton transfer from the lytic water to the γ-phosphate through active site residues is an important part of the kinetic bottleneck; several hydrolysis pathways that feature distinct proton transfer routes are found to have similar free energy barriers, suggesting a significant degree of plasticity in the hydrolysis mechanism. Comparison of hydrolysis in the pre-powerstroke state and the closed post-rigor model suggests that optimization of residues beyond the active site for electrostatic stabilization and preorganization is likely important to enzyme design.

Poorly understood aspects of striated muscle contraction

Muscle contraction results from cyclic interactions between the contractile proteins myosin and actin, driven by the turnover of adenosine triphosphate (ATP). Despite intense studies, several molecular events in the contraction process are poorly understood, including the relationship between force-generation and phosphate-release in the ATP-turnover. Different aspects of the force-generating transition are reflected in the changes in tension development by muscle cells, myofibrils and single molecules upon changes in temperature, altered phosphate concentration, or length perturbations. It has been notoriously difficult to explain all these events within a given theoretical framework and to unequivocally correlate observed events with the atomic structures of the myosin motor. Other incompletely understood issues include the role of the two heads of myosin II and structural changes in the actin filaments as well as the importance of the three-dimensional order. We here review these issues in relation to controversies regarding basic physiological properties of striated muscle. We also briefly consider actomyosin mutation effects in cardiac and skeletal muscle function and the possibility to treat these defects by drugs.

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

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

Alternative relay and converter domains tune native muscle myosin isoform function in Drosophila

Myosin isoforms help define muscle-specific contractile and structural properties. Alternative splicing of myosin heavy chain gene transcripts in Drosophila melanogaster yields muscle-specific isoforms and highlights alternative domains that fine-tune myosin function. To gain insight into how native myosin is tuned, we expressed three embryonic myosin isoforms in indirect flight muscles lacking endogenous myosin. These isoforms differ in their relay and/or converter domains. We analyzed isoform-specific ATPase activities, in vitro actin motility and myofibril structure/stability. We find that dorsal acute body wall muscle myosin (EMB-9c11d) shows a significant increase in MgATPase V(max) and actin sliding velocity, as well as abnormal myofibril assembly compared to cardioblast myosin (EMB-11d). These properties differ as a result of alternative exon-9-encoded relay domains that are hypothesized to communicate signals among the ATP-binding pocket, actin-binding site and the converter domain. Further, EMB-11d shows significantly reduced levels of basal Ca- and MgATPase as well as MgATPase V(max) compared to embryonic body wall muscle isoform (EMB) (expressed in a multitude of body wall muscles). EMB-11d also induces increased actin sliding velocity and stabilizes myofibril structure compared to EMB. These differences arise from exon-11-encoded alternative converter domains that are proposed to reposition the lever arm during the power and recovery strokes. We conclude that relay and converter domains of native myosin isoforms fine-tune ATPase activity, actin motility and muscle ultrastructure. This verifies and extends previous studies with chimeric molecules and indicates that interactions of the relay and converter during the contractile cycle are key to myosin-isoform-specific kinetic and mechanical functions.

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.

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.

Structural mechanism of the ATP-induced dissociation of rigor myosin from actin

Myosin is a true nanomachine, which produces mechanical force from ATP hydrolysis by cyclically interacting with actin filaments in a four-step cycle. The principle underlying each step is that structural changes in separate regions of the protein must be mechanically coupled. The step in which myosin dissociates from tightly bound actin (the rigor state) is triggered by the 30 Å distant binding of ATP. Large conformational differences between the crystal structures make it difficult to perceive the coupling mechanism. Energetically accessible transition pathways computed at atomic detail reveal a simple coupling mechanism for the reciprocal binding of ATP and actin.

Predicting order of conformational changes during protein conformational transitions using an interpolated elastic network model

The decryption of sequence of structural events during protein conformational transitions is essential to a detailed understanding of molecular functions of various biological nanomachines. Coarse-grained models have proven useful by allowing highly efficient simulations of protein conformational dynamics. By combining two coarse-grained elastic network models constructed based on the beginning and end conformations of a transition, we have developed an interpolated elastic network model to generate a transition pathway between the two protein conformations. For validation, we have predicted the order of local and global conformational changes during key ATP-driven transitions in three important biological nanomachines (myosin, F(1) ATPase and chaperonin GroEL). We have found that the local conformational change associated with the closing of active site precedes the global conformational change leading to mechanical motions. Our finding is in good agreement with the distribution of intermediate experimental structures, and it supports the importance of local motions at active site to drive or gate various conformational transitions underlying the workings of a diverse range of biological nanomachines.

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.

Atomically detailed simulation of the recovery stroke in myosin by Milestoning

Myosin II is a molecular motor that converts chemical to mechanical energy and enables muscle operations. After a power stroke, a recovery transition completes the cycle and returns the molecular motor to its prestroke state. Atomically detailed simulations in the framework of the Milestoning theory are used to calculate kinetics and mechanisms of the recovery stroke. Milestoning divides the process into transitions between hyper-surfaces (Milestones) along a reaction coordinate. Decorrelation of dynamics between sequential Milestones is assumed, which speeds up the atomically detailed simulations by a factor of millions. Two hundred trajectories of myosin with explicit water solvation are used to sample transitions between sequential pairs of Milestones. Collective motions of hundreds of atoms are described at atomic resolution and at the millisecond time scale. The experimentally measured transition time of about a millisecond is in good agreement with the computed time. The simulations support a sequential mechanism. In the first step the P-loop and switch 2 close on the ATP and in the second step the mechanical relaxation is induced via the relay and the SH1 helices. We propose that the entropy of switch 2 helps to drive the power stroke. Secondary structure elements are progressing through a small number of discrete states in a network of activated transitions and are assisted by side chain flips between rotameric states. The few-state sequential mechanism is likely to enhance the efficiency of the relaxation reducing the probability of off-pathway intermediates.

Analysis of the interaction of the nucleotide base with myosin and the effect on substrate efficacy

A wide variety of purine- and pyrimidine-based nucleotides can serve as a substrate for actomyosin mechanics, but with varying effectiveness. To understand the myosin-ATP interaction and in particular, the interactions with the base, we have used molecular dynamics simulations to model the interactions of myosin with ATP, CTP, UTP, aza-ATP, ITP, and GTP (in decreasing order of effectiveness as a substrate for the generation of motility) docked at the active site. The simulations with ATP, and x-ray structures, show a triad of conserved amino acids lining the nucleotide site that form a cyclical chain of nucleotide-protein hydrogen bonding interactions: ATP --> Y135 --> Y116 --> N188 --> ATP. Mechanical efficacy of a substrate correlates with its ability to maintain this coordination. Simulations modeling the active site of other myosin isoforms with different amino acids in the triad likewise imply that the amino acid composition at the nucleotide site could modulate function. The modeling has predictive power. In silico mutation experiments suggest mutations that would enhance GTP as a substrate for myosin while simultaneously making ATP a less effective substrate.

Genetic and clinical profile of Indian patients of idiopathic restrictive cardiomyopathy with and without hypertrophy

Both idiopathic restrictive cardiomyopathy (IRCM) and hypertrophic cardiomyopathy (HCM) are part of the same disease spectrum and are due to sarcomeric gene mutations. A patient with restrictive physiology without left ventricular hypertrophy (LVH) would be diagnosed as IRCM, while one with LVH would be diagnosed as HCM with restrictive physiology. We studied a group of patients with restrictive physiology for mutations in beta-myosin heavy chain (MYH7) and troponin I (TNNI3) gene. Consecutive probands in the HCM and IRCM cohort over a 4-year period were considered for this study. These included 10 IRCM and 102 HCM patients. All were Asian Indians. Among the 17 patients who had restrictive physiology 10 were IRCM patients and seven were HCM patients. Of the HCM patients, seven (6.9%) had restrictive physiology. Mean age of these 17 patients was 40.1 +/- 19.2 years (range: 15-67 ), six (35.3%) were males. Maximal left ventricular wall thickness of the seven HCM probands was 20.7 +/- 5.2 mm (range: 16-31), while it was normal in the IRCM probands. Ten probands (58.8%) were in NYHA class III or IV. Seven patients (41.2%) had atrial fibrillation. All the probands were screened for mutations in selected exons of MYH7 and TNNI3 genes. One IRCM patient was found to have p.Arg721Lys mutation in the MYH7 gene. She died due to progressive congestive cardiac failure at the age of 47 years. One HCM proband with a maximal left ventricular wall thickness of 17 mm had p.Arg192His mutation in the TNNI3 gene. She had features consistent with restrictive physiology. Her father and sister had died of restrictive cardiomyopathy. IRCM and HCM with restrictive physiology, both are part of the clinical expression of MYH7 and TNNI3 mutations and lead to worse clinical onset and progression of the disease.

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.

Genotype phenotype correlations of cardiac beta-myosin heavy chain mutations in Indian patients with hypertrophic and dilated cardiomyopathy

The aim of the current study was to determine the frequency of mutations in the beta-myosin heavy chain gene (MYH7) in a cohort of hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM) and their families, and to investigate correlations between genotype and phenotype. About 130 consecutive patients diagnosed with HCM or DCM (69 with HCM and 61 with DCM) attending the cardiology clinic of Post Graduate Institute of Medical Education and Research were screened for mutations in the MYH7 gene. The control group for genetic studies consisted of 100 healthy subjects. We report 14 mutations in 6 probands (5 probands in HCM and 1 proband in DCM) and their family members. Out of these 6 mutations, 3 are new and are being reported for the first time. One known mutation (p.Gly716Arg) was found to be "de novo" which resulted in severe asymmetric septal hypertrophy (31 mm) and resulted in the sudden cardiac death (SCD) of the proband at the age of 21 years. Further, a DCM causing novel mutation p.Gly377Ser was identified which resulted in the milder phenotype. The present study shows that there is genetic and phenotypic heterogeneity of cardiomyopathies in Indian population. Further, the location and type of mutation in a given sarcomeric gene determines the severity and phenotypic plasticity in cardiomyopathies.

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 A away.

Can a synthetic thread act as an electrochemically switchable molecular device?

Here we report that at room temperature in acetonitrile after the reduction of the naphthalimide-site, a synthetic molecular thread undergoes a complete conformational change which makes possible an efficient conversion of chemical energy into mechanical work; such results point out the ability of the thread to act as a molecular device under electrochemical control.

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.

The mechanism of the reverse recovery step, phosphate release, and actin activation of Dictyostelium myosin II

The rate-limiting step of the myosin basal ATPase (i.e. in absence of actin) is assumed to be a post-hydrolysis swinging of the lever arm (reverse recovery step), that limits the subsequent rapid product release steps. However, direct experimental evidence for this assignment is lacking. To investigate the binding and the release of ADP and phosphate independently from the lever arm motion, two single tryptophan-containing motor domains of Dictyostelium myosin II were used. The single tryptophans of the W129+ and W501+ constructs are located at the entrance of the nucleotide binding pocket and near the lever arm, respectively. Kinetic experiments show that the rate-limiting step in the basal ATPase cycle is indeed the reverse recovery step, which is a slow equilibrium step (k(forward) = 0.05 s(-1), k(reverse) = 0.15 s(-1)) that precedes the phosphate release step. Actin directly activates the reverse recovery step, which becomes practically irreversible in the actin-bound form, triggering the power stroke. Even at low actin concentrations the power stroke occurs in the actin-attached states despite the low actin affinity of myosin in the pre-power stroke conformation.

The post-rigor structure of myosin VI and implications for the recovery stroke

Myosin VI has an unexpectedly large swing of its lever arm (powerstroke) that optimizes its unique reverse direction movement. The basis for this is an unprecedented rearrangement of the subdomain to which the lever arm is attached, referred to as the converter. It is unclear at what point(s) in the myosin VI ATPase cycle rearrangements in the converter occur, and how this would effect lever arm position. We solved the structure of myosin VI with an ATP analogue (ADP.BeF3) bound in its nucleotide-binding pocket. The structure reveals that no rearrangement in the converter occur upon ATP binding. Based on previously solved myosin structures, our structure suggests that no reversal of the powerstroke occurs during detachment of myosin VI from actin. The structure also reveals novel features of the myosin VI motor that may be important in maintaining the converter conformation during detachment from actin, and other features that may promote rapid rearrangements in the structure following actin detachment that enable hydrolysis of ATP.

Myosin V movement: lessons from molecular dynamics studies of IQ peptides in the lever arm

Myosin V moves along actin filaments by an arm-over-arm motion, known as the lever mechanism. Each of its arms is composed of six consecutive IQ peptides that bind light chain proteins, such as calmodulin or calmodulin-like proteins. We have employed a multistage approach in order to investigate the mechanochemical structural basis of the movement of myosin V from the budding yeast Saccharomyces cerevisiae. For that purpose, we previously carried out molecular dynamics simulations of the Mlc1p-IQ2 and the Mlc1p-IQ4 protein-peptide complexes, and the present study deals with the structures of the IQ peptides when stripped from the Mlc1p protein. We have found that the crystalline structure of the IQ2 peptide retains a stable rodlike configuration in solution, whereas that of the IQ4 peptide grossly deviates from its X-ray conformation exhibiting an intrinsic tendency to curve and bend. The refolding process of the IQ4 peptide is initially driven by electrostatic interactions followed by nonpolar stabilization. Its bending appears to be affected by the ionic strength, when ionic strength higher than approximately 300 mM suppresses it from flexing. Considering that a poly-IQ sequence is the lever arm of myosin V, we suggest that the arm may harbor a joint, localized within the IQ4 sequence, enabling the elasticity of the neck of myosin V. Given that a poly-IQ sequence is present at the entire class of myosin V and the possibility that the yeast's myosin V molecule can exist either as a nonprocessive monomer or as a processive dimer depending on conditions (Krementsova, E. B., Hodges, A. R., Lu, H., and Trybus, K. M. (2006) J. Biol. Chem. 281, 6079-6086), our observations may account for a general structural feature for the myosins' arm embedded flexibility.

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

The Structural Coupling between ATPase Activation and Recovery Stroke in the Myosin II Motor

Before the myosin motor head can perform the next power stroke, it undergoes a large conformational transition in which the converter domain, bearing the lever arm, rotates approximately 65 degrees . Simultaneous with this "recovery stroke," myosin activates its ATPase function by closing the Switch-2 loop over the bound ATP. This coupling between the motions of the converter domain and of the 40 A-distant Switch-2 loop is essential to avoid unproductive ATP hydrolysis. The coupling mechanism is determined here by finding a series of optimized intermediates between crystallographic end structures of the recovery stroke (Dictyostelium discoideum), yielding movies of the transition at atomic detail. The successive formation of two hydrogen bonds by the Switch-2 loop is correlated with the successive see-saw motions of the relay and SH1 helices that hold the converter domain. SH1 helix and Switch-2 loop communicate via a highly conserved loop that wedges against the SH1-helix upon Switch-2 closing.

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

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.