Molecular dynamics analysis of structural factors influencing back door pi release in myosin

Mechanisms of myosin II force generation: insights from novel experimental techniques and approaches

Myosin II is a molecular motor that converts chemical energy derived from ATP hydrolysis into mechanical work. Myosin II isoforms are responsible for muscle contraction and a range of cell functions relying on the development of force and motion. When the motor attaches to actin, ATP is hydrolyzed and inorganic phosphate (Pi) and ADP are released from its active site. These reactions are coordinated with changes in the structure of myosin, promoting the so-called "power stroke" that causes the sliding of actin filaments. The general features of the myosin-actin interactions are well accepted, but there are critical issues that remain poorly understood, mostly due to technological limitations. In recent years, there has been a significant advance in structural, biochemical, and mechanical methods that have advanced the field considerably. New modeling approaches have also allowed researchers to understand actomyosin interactions at different levels of analysis. This paper reviews recent studies looking into the interaction between myosin II and actin filaments, which leads to power stroke and force generation. It reviews studies conducted with single myosin molecules, myosins working in filaments, muscle sarcomeres, myofibrils, and fibers. It also reviews the mathematical models that have been used to understand the mechanics of myosin II in approaches focusing on single molecules to ensembles. Finally, it includes brief sections on translational aspects, how changes in the myosin motor by mutations and/or posttranslational modifications may cause detrimental effects in diseases and aging, among other conditions, and how myosin II has become an emerging drug target.

Single-molecule tracking reveals dual front door/back door inhibition of Cel7A cellulase by its product cellobiose

Degrading cellulose is a key step in the processing of lignocellulosic biomass into bioethanol. Cellobiose, the disaccharide product of cellulose degradation, has been shown to inhibit cellulase activity, but the mechanisms underlying product inhibition are not clear. We combined single-molecule imaging and biochemical investigations with the goal of revealing the mechanism by which cellobiose inhibits the activity of Trichoderma reesei Cel7A, a well-characterized exo-cellulase. We find that cellobiose slows the processive velocity of Cel7A and shortens the distance moved per encounter; effects that can be explained by cellobiose binding to the product release site of the enzyme. Cellobiose also strongly inhibits the binding of Cel7A to immobilized cellulose, with a Ki of 2.1 mM. The isolated catalytic domain (CD) of Cel7A was also inhibited to a similar degree by cellobiose, and binding of an isolated carbohydrate-binding module to cellulose was not inhibited by cellobiose, suggesting that cellobiose acts on the CD alone. Finally, cellopentaose inhibited Cel7A binding at micromolar concentrations without affecting the enzyme's velocity of movement along cellulose. Together, these results suggest that cellobiose inhibits Cel7A activity both by binding to the "back door" product release site to slow activity and to the "front door" substrate-binding tunnel to inhibit interaction with cellulose. These findings point to strategies for engineering cellulases to reduce product inhibition and enhance cellulose degradation, supporting the growth of a sustainable bioeconomy.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

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

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

Embedding Quantum Statistical Excitations in a Classical Force Field

Quantum-mechanically driven charge polarization and charge transfer are ubiquitous in biomolecular systems, controlling reaction rates, allosteric interactions, ligand-protein binding, membrane transport, and dynamically driven structural transformations. Molecular dynamics (MD) simulations of these processes require quantum mechanical (QM) information in order to accurately describe their reactive dynamics. However, current techniques-empirical force fields, subsystem approaches, ab initio MD, and machine learning-vary in their ability to achieve a consistent chemical description across multiple atom types, and at scale. Here we present a physics-based, atomistic force field, the ensemble DFT charge-transfer embedded-atom method, in which QM forces are described at a uniform level of theory across all atoms, avoiding the need for explicit solution of the Schrödinger equation or large, precomputed training data sets. Coupling between the electronic and atomistic length scales is effected through an ensemble density functional theory formulation of the embedded-atom method originally developed for elemental materials. Charge transfer is expressed in terms of ensembles of ionic state basis densities of individual atoms, and charge polarization, in terms of atomic excited-state basis densities. This provides a highly compact yet general representation of the force field, encompassing both local and system-wide effects. Charge rearrangement is realized through the evolution of ensemble weights, adjusted at each dynamical time step via chemical potential equalization.

Long-range coupling between ATP-binding and lever-arm regions in myosin via dielectric allostery

A protein molecule is a dielectric substance, so the binding of a ligand is expected to induce dielectric response in the protein molecule, considering that ligands are charged or polar in general. We previously reported that binding of adenosine triphosphate (ATP) to molecular motor myosin actually induces such a dielectric response in myosin due to the net negative charge of ATP. By this dielectric response, referred to as "dielectric allostery," spatially separated two regions in myosin, the ATP-binding region and the actin-binding region, are allosterically coupled. In this study, from the statistically stringent analyses of the extensive molecular dynamics simulation data obtained in the ATP-free and the ATP-bound states, we show that there exists the dielectric allostery that transmits the signal of ATP binding toward the distant lever-arm region. The ATP-binding-induced electrostatic potential change observed on the surface of the main domain induced a movement of the converter subdomain from which the lever arm extends. The dielectric response was found to be caused by an underlying large-scale concerted rearrangement of the electrostatic bond network, in which highly conserved charged/polar residues are involved. Our study suggests the importance of the dielectric property for molecular machines in exerting their function.

Physical driving force of actomyosin motility based on the hydration effect

We propose a driving force hypothesis based on previous thermodynamics, kinetics and structural data as well as additional experiments and calculations presented here on water-related phenomena in the actomyosin systems. Although Szent-Györgyi pointed out the importance of water in muscle contraction in 1951, few studies have focused on the water science of muscle because of the difficulty of analyzing hydration properties of the muscle proteins, actin, and myosin. The thermodynamics and energetics of muscle contraction are linked to the water-mediated regulation of protein-ligand and protein-protein interactions along with structural changes in protein molecules. In this study, we assume the following two points: (1) the periodic electric field distribution along an actin filament (F-actin) is unidirectionally modified upon binding of myosin subfragment 1 (M or myosin S1) with ADP and inorganic phosphate Pi (M.ADP.Pi complex) and (2) the solvation free energy of myosin S1 depends on the external electric field strength and the solvation free energy of myosin S1 in close proximity to F-actin can become the potential force to drive myosin S1 along F-actin. The first assumption is supported by integration of experimental reports. The second assumption is supported by model calculations utilizing molecular dynamics (MD) simulation to determine solvation free energies of a small organic molecule and two small proteins. MD simulations utilize the energy representation method (ER) and the roughly proportional relationship between the solvation free energy and the solvent-accessible surface area (SASA) of the protein. The estimated driving force acting on myosin S1 is as high as several piconewtons (pN), which is consistent with the experimentally observed force.

Modulation of Skeletal Muscle Contraction by Myosin Phosphorylation

The striated muscle sarcomere is a highly organized and complex enzymatic and structural organelle. Evolutionary pressures have played a vital role in determining the structure-function relationship of each protein within the sarcomere. A key part of this multimeric assembly is the light chain-binding domain (LCBD) of the myosin II motor molecule. This elongated "beam" functions as a biological lever, amplifying small interdomain movements within the myosin head into piconewton forces and nanometer displacements against the thin filament during the cross-bridge cycle. The LCBD contains two subunits known as the essential and regulatory myosin light chains (ELC and RLC, respectively). Isoformic differences in these respective species provide molecular diversity and, in addition, sites for phosphorylation of serine residues, a highly conserved feature of striated muscle systems. Work on permeabilized skeletal fibers and thick filament systems shows that the skeletal myosin light chain kinase catalyzed phosphorylation of the RLC alters the "interacting head motif" of myosin motor heads on the thick filament surface, with myriad consequences for muscle biology. At rest, structure-function changes may upregulate actomyosin ATPase activity of phosphorylated cross-bridges. During activation, these same changes may increase the Ca2+ sensitivity of force development to enhance force, work, and power output, outcomes known as "potentiation." Thus, although other mechanisms may contribute, RLC phosphorylation may represent a form of thick filament activation that provides a "molecular memory" of contraction. The clinical significance of these RLC phosphorylation mediated alterations to contractile performance of various striated muscle systems are just beginning to be understood. © 2017 American Physiological Society. Compr Physiol 7:171-212, 2017.

Dielectric Allostery of Protein: Response of Myosin to ATP Binding

Protein uses allostery to execute biological function. The physical mechanism underlying the allostery has long been studied, with the focus on the mechanical response by ligand binding. Here, we highlight the electrostatic response, presenting an idea of "dielectric allostery". We conducted molecular dynamics simulations of myosin, a motor protein with allostery, and analyzed the response to ATP binding which is a crucial step in force-generating function, forcing myosin to unbind from the actin filament. We found that the net negative charge of ATP causes a large-scale, anisotropic dielectric response in myosin, altering the electrostatic potential in the distant actin-binding region and accordingly retracting a positively charged actin-binding loop. A large-scale rearrangement of electrostatic bond network was found to occur upon ATP binding. Since proteins are dielectric and ligands are charged/polar in general, the dielectric allostery might underlie a wide spectrum of functions by proteins.

How Myosin Generates Force on Actin Filaments

How myosin interacts with actin to generate force is a subject of considerable controversy. The major debate centers on understanding at what point in force generation the inorganic phosphate is released with respect to the lever arm swing, or powerstroke. Resolving the controversy is essential for understanding how force is produced as well as the mechanisms underlying disease-causing mutations in myosin. Recent structural insights into the powerstroke have come from a high-resolution structure of myosin in a previously unseen state and from an electron cryomicroscopy (cryo-EM) 3D reconstruction of the actin-myosin-MgADP complex. Here, we argue that seemingly contradictory data from time-resolved fluorescence resonance energy transfer (FRET) studies can be reconciled, and we put forward a model for myosin force generation on actin.

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

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

Orientation of the N- and C-terminal lobes of the myosin regulatory light chain in cardiac muscle

The orientations of the N- and C-terminal lobes of the cardiac isoform of the myosin regulatory light chain (cRLC) in the fully dephosphorylated state in ventricular trabeculae from rat heart were determined using polarized fluorescence from bifunctional sulforhodamine probes. cRLC mutants with one of eight pairs of surface-accessible cysteines were expressed, labeled with bifunctional sulforhodamine, and exchanged into demembranated trabeculae to replace some of the native cRLC. Polarized fluorescence data from the probes in each lobe were combined with RLC crystal structures to calculate the lobe orientation distribution with respect to the filament axis. The orientation distribution of the N-lobe had three distinct peaks (N1–N3) at similar angles in relaxation, isometric contraction, and rigor. The orientation distribution of the C-lobe had four peaks (C1–C4) in relaxation and isometric contraction, but only two of these (C2 and C4) remained in rigor. The N3 and C4 orientations are close to those of the corresponding RLC lobes in myosin head fragments bound to isolated actin filaments in the absence of ATP (in rigor), but also close to those of the pair of heads folded back against the filament surface in isolated thick filaments in the so-called J-motif conformation. The N1 and C1 orientations are close to those expected for actin-bound myosin heads with their light chain domains in a pre-powerstroke conformation. The N2 and C3 orientations have not been observed previously. The results show that the average change in orientation of the RLC region of the myosin heads on activation of cardiac muscle is small; the RLC regions of most heads remain in the same conformation as in relaxation. This suggests that the orientation of the dephosphorylated RLC region of myosin heads in cardiac muscle is primarily determined by an interaction with the thick filament surface.

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

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

Phosphate release coupled to rotary motion of F1-ATPase

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

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

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

Muscle contraction mechanism based on single molecule measurements

Single molecule measurements have shown that a muscle myosin step is driven by biased Brownian movement. Furthermore, they have also demonstrated that in response to strain in the backward direction a detached myosin head preferentially attaches to the forward direction due to an accelerated transition from a weak binding to strong binding state. Because they are consistent with the original Huxley model for muscle contraction, we have built a model that describes macroscopic muscle characteristics based on these single molecule results.

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

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

Minimum energy reaction profiles for ATP hydrolysis in myosin

The minimum energy reaction profiles corresponding to two possible reaction mechanisms of adenosine triphosphate (ATP) hydrolysis in myosin are computed in this work within the framework of the quantum mechanics-molecular mechanics (QM/MM) method by using the same partitioning of the model system to the QM and MM parts and the same computational protocol. On the first reaction route, one water molecule performs nucleophilic attack at the phosphorus center P(γ) from ATP while the second water molecule in the closed protein cleft serves as a catalytic base assisted by the Glu residue from the myosin salt bridge. According to the present QM/MM calculations consistent with the results of kinetic studies this reaction pathway is characterized by a low activation energy barrier about 10 kcal/mol. The computed activation energy barrier for the second mechanism, which assumes the penta-coordinated oxyphosphorane transition state upon involvement of single water molecule in the reaction, is considerably higher than that for the two-water mechanism.

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.

Proteomic analysis of the enterocyte brush border

The brush border domain at the apex of intestinal epithelial cells is the primary site of nutrient absorption in the intestinal tract and the primary surface of interaction with microbes that reside in the lumen. Because the brush border is positioned at such a critical physiological interface, we set out to create a comprehensive list of the proteins that reside in this domain using shotgun mass spectrometry. The resulting proteome contains 646 proteins with diverse functions. In addition to the expected collection of nutrient processing and transport components, we also identified molecules expected to function in the regulation of actin dynamics, membrane bending, and extracellular adhesion. These results provide a foundation for future studies aimed at defining the molecular mechanisms underpinning brush border assembly and function.

Pi release from myosin: a simulation analysis of possible pathways

The release of phosphate (Pi) is an important element in actomyosin function and has been shown to be accelerated by the binding of myosin to actin. To provide information about the structural elements important for Pi release, possible escape pathways from various isolated myosin II structures have been determined by molecular dynamics simulations designed for studying such slow processes. The residues forming the pathways were identified and their role was evaluated by mutant simulations. Pi release is slow in the pre-powerstroke structure, an important element in preventing the powerstroke prior to actin binding, and is much more rapid for Pi modeled into the post-rigor and rigor-like structures. The previously proposed backdoor route is dominant in the pre-powerstroke and post-rigor states, whereas a different path is most important in the rigor-like state. This finding suggests a mechanism for the actin-activated acceleration of Pi release.

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.

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.

Models of protein linear molecular motors for dynamic nanodevices

Protein molecular motors are natural nano-machines that convert the chemical energy from the hydrolysis of adenosine triphosphate into mechanical work. These efficient machines are central to many biological processes, including cellular motion, muscle contraction and cell division. The remarkable energetic efficiency of the protein molecular motors coupled with their nano-scale has prompted an increasing number of studies focusing on their integration in hybrid micro- and nanodevices, in particular using linear molecular motors. The translation of these tentative devices into technologically and economically feasible ones requires an engineering, design-orientated approach based on a structured formalism, preferably mathematical. This contribution reviews the present state of the art in the modelling of protein linear molecular motors, as relevant to the future design-orientated development of hybrid dynamic nanodevices.

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.

Characterization of the pre-force-generation state in the actomyosin cross-bridge cycle

Myosin is an actin-based motor protein that generates force by cycling between actin-attached (strong binding: ADP or rigor) and actin-detached (weak binding: ATP or ADP.P(i)) states during its ATPase cycle. However, it remains unclear what specific conformational changes in the actin binding site take place on binding to actin, and how these structural changes lead to product release and the production of force and motion. We studied the dynamics of the actin binding region of myosin V by using fluorescence resonance energy transfer (FRET) to monitor conformational changes in the upper-50-kDa domain of the actin binding cleft in the weak and strong actin binding states. Steady-state and lifetime data monitoring the FRET signal suggest that the cleft is in a more open conformation in the weak actin binding states. Transient kinetic experiments suggest that a rapid conformational change occurs, which is consistent with cleft closure before actin-activated phosphate release. Our results have identified a pre-force-generation actomyosin ADP.P(i) state, and suggest force generation may occur from a state not yet seen by crystallography in which the actin binding cleft and the nucleotide binding pocket are closed. Computational modeling uncovers dramatic changes in the rigidity of the upper-50-kDa domain in different nucleotide states, which suggests that the intrinsic flexibility of this domain allows myosin motors to accomplish simultaneous tight nucleotide binding (closed nucleotide binding pocket) and high-affinity actin binding (closed actin binding cleft).

Acetylcholinesterase: how is structure related to function?

In accordance with its biological role, termination of neurotransmission at cholinergic synapses by rapid hydrolysis of the neurotransmitter, acetylcholine, acetylcholinesterase is one of nature's most efficient enzymes. Solution of its three-dimensional structure revealed that its active site is located at the bottom of a deep and narrow gorge. Such an architecture was unanticipated in view of its high turnover number. The present review examines how the highly specialized structure of acetylcholinesterase, with its sequestered active site, contributes to its catalytic efficacy, and discusses how the traffic of substrate and products to and from the active site is controlled.

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.

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.

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.

Human deafness mutation E385D disrupts the mechanochemical coupling and subcellular targeting of myosin-1a

Missense mutations in the membrane-binding actin-based motor protein, myosin-1a (Myo1a), have recently been linked to sensorineural deafness in humans. One of these mutations, E385D, impacts a residue in the switch II region of the motor domain that is present in virtually all members of the myosin superfamily. We sought to examine the impact of E385D on the function of Myo1a, both in terms of mechanochemical activity and ability to target to actin-rich microvilli in polarized epithelial cells. While E385D-Myo1a demonstrated actin-activated ATPase activity, the V(MAX) was reduced threefold relative to wild-type. Despite maintaining an active mechanochemical cycle, E385D-Myo1a was unable to move actin in the sliding filament assay. Intriguingly, when an enhanced-green-fluorescent-protein-tagged form of E385D-Myo1a was stably expressed in polarized epithelial cells, this mutation abolished the microvillar targeting normally demonstrated by wild-type Myo1a. Notably, these data are the first to suggest that mechanical activity is essential for proper localization of Myo1a in microvilli. These studies also provide a unique example of how even the most mild substitution of invariant switch II residues can effectively uncouple enzymatic and mechanical activity of the myosin motor domain.

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

Binding modes of decavanadate to myosin and inhibition of the actomyosin ATPase activity

Decavanadate, a vanadate oligomer, is known to interact with myosin and to inhibit the ATPase activity, but the putative binding sites and the mechanism of inhibition are still to be clarified. We have previously proposed that the decavanadate (V(10)O(28)(6-)) inhibition of the actin-stimulated myosin ATPase activity is non-competitive towards both actin and ATP. A likely explanation for these results is that V(10) binds to the so-called back-door at the end of the Pi-tube opposite to the nucleotide-binding site. In order to further investigate this possibility, we have carried out molecular docking simulations of the V(10) oligomer on three different structures of the myosin motor domain of Dictyostelium discoideum, representing distinct states of the ATPase cycle. The results indicate a clear preference of V(10) to bind at the back-door, but only on the "open" structures where there is access to the phosphate binding-loop. It is suggested that V(10) acts as a "back-door stop" blocking the closure of the 50-kDa cleft necessary to carry out ATP-gamma-phosphate hydrolysis. This provides a simple explanation to the non-competitive behavior of V(10) and spurs the use of the oligomer as a tool to elucidate myosin back-door conformational changes in the process of muscle contraction.

High-resolution cryo-EM maps show the nucleotide binding pocket of KIF1A in open and closed conformations

Kinesin is an ATP-driven microtubule (MT)-based motor fundamental to organelle transport. Although a number of kinesin crystal structures have been solved, the structural evidence for coupling between the bound nucleotide and the conformation of kinesin is elusive. In addition, the structural basis of the MT-induced ATPase activity of kinesin is not clear because of the absence of the MT in the structure. Here, we report cryo-electron microscopy structures of the monomeric kinesin KIF1A-MT complex in two nucleotide states at about 10 A resolution, sufficient to reveal the secondary structure. These high-resolution maps visualized clear structural changes that suggest a mechanical pathway from the nucleotide to the neck linker via the motor core rotation. In addition, new nucleotide binding pocket conformations are observed that are different from X-ray crystallographic structures; it is closed in the 5'-adenylyl-imidodiphosphate state, but open in the ADP state. These results suggest a structural model of biased diffusion movement of monomeric kinesin motor.

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

Detailed residue-wise interactions involved in the binding of myosin to actin in the rigor conformation without nucleotides have been examined using molecular dynamics simulations of the chicken skeletal myosin head complexed with two actin monomers, based on the cryo-microscopic model of Holmes et al. (Nature 2003;425:423-427). The overall interaction is largely electrostatic in nature, because of the charged residues in the four loops surrounding the central primary binding site. The 50k/20k loop, disordered in crystal structures and in simulations of free myosin in solution, was found to be in a conformation stabilized with 1 - 2 internal salt bridges. The cardiomyopathy loop forms 2 - 3 interprotein salt bridges with actin monomers upon binding, whereas its Arg405 residue, the mutation site associated with the hypertrophic cardiomyopathy, forms a strong salt bridge with Glu605 in the neighboring helix away from actin in the actin-bound myosin. The myopathy loop of the R405Q mutant maintains a high degree of two-strand beta-sheet character when bound to actin with the corresponding salt bridges broken.

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.

Dynamics of the upper 50-kDa domain of myosin V examined with fluorescence resonance energy transfer

The upper 50-kDa region of myosin may be critical for coupling between the nucleotide- and actin-binding regions. We introduced a tetracysteine motif in the upper 50-kDa domain (residues 292-297) of myosin V containing a single IQ domain (MV 1IQ), allowing us to label this site with the fluorescein biarscenical hairpin-binding dye (FlAsH) (MV 1IQ FlAsH). The enzymatic properties of MV 1IQ FlAsH were similar to those of unlabeled MV 1IQ except for a 3-fold reduced ADP-release rate. MV 1IQ FlAsH was also capable of moving actin filaments in the in vitro motility assay. To examine rotation of the upper 50-kDa region, we determined the difference in the degree of energy transfer from N-methylanthraniloyl (mant)-labeled nucleotides to FlAsH in both steady-state and transient kinetic experiments. The energy transfer efficiency was higher with mant-ATP (0.65 +/- 0.02) compared with mant-ADP (0.55 +/- 0.02) in the absence of actin. Stopped-flow measurements suggested that the energy transfer efficiency decreased with phosphate release (0.04 s(-1)) in the absence of actin. In contrast, upon mixing MV 1IQ FlAsH in the ADP.P(i) state with actin, a decrease in the energy transfer signal was observed at a rate of 13 s(-1), similar to the ADP release rate. Our results demonstrate there was no change in the energy transfer signal upon actin-activated phosphate release and suggest that actin binding alters the dynamics of the upper 50-kDa region, which may be critical for the ability of myosin to bind tightly to both ADP and actin.

The motor mechanism of myosin V: insights for muscle contraction

It is 50 years since the sliding of actin and myosin filaments was proposed as the basis of force generation and shortening in striated muscle. Although this is now generally accepted, the detailed molecular mechanism of how myosin uses adenosine triphosphate to generate force during its cyclic interaction with actin is only now being unravelled. New insights have come from the unconventional myosins, especially myosin V. Myosin V is kinetically tuned to allow movement on actin filaments as a single molecule, which has led to new kinetic, mechanical and structural data that have filled in missing pieces of the actomyosin-chemo-mechanical transduction puzzle.

Molecular dynamics simulations of evolved collective motions of atoms in the myosin motor domain upon perturbation of the ATPase pocket

A crucial point for mechanical force generation in actomyosin systems is how the energy released by ATP hydrolysis in the myosin motor domain gives rise to the movement of the myosin head along the actin filament. We assumed the signal of the ATP hydrolysis to be transmitted as modulated atomic vibrations from the nucleotide-binding site throughout the myosin head, and carried out 1-ns all-atom molecular dynamics simulations for that signal transmission. We distributed the released energy to atoms located around the ATPase pocket as kinetic energies and examined how the effect of disturbance extended throughout the motor domain. The result showed that the disturbance signal extended over the motor domain in 150 ps and induced slowly varying collective motions of atoms at the actin-binding site and the junction with the neck, both of which are relevant to the movement of the myosin head along the actin filament. We also performed a principal component analysis of thermal atomic motions for the motor domain, and the first principal component was consistent with the response to the disturbance given to the ATPase pocket.