An intermediate along the recovery stroke of myosin VI revealed by X-ray crystallography and molecular dynamics

Hidden complexity of α7 nicotinic acetylcholine receptor desensitization revealed by MD simulations and Markov state modeling

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that plays an important role in neuronal signaling throughout the nervous system. Its implication in neurological disorders and inflammation has spurred the development of numerous compounds that enhance channel activation. However, the therapeutic potential of these compounds has been limited by the characteristically fast desensitization of the α7 receptor. Using recent high-resolution structures from cryo-EM, and all-atom molecular dynamic simulations augmented by Markov state modeling, here we explore the mechanism of α7 receptor desensitization and its implication on allosteric modulation. The results provide a precise characterization of the desensitization gate and illuminate the mechanism of ion-pore opening/closing with an agonist bound. In addition, the simulations reveal the existence of a short-lived, open-channel intermediate between the activated and desensitized states that rationalizes the paradoxical pharmacology of the L247T mutant and may be relevant to type-II allosteric modulation. This analysis provides an interpretation of the signal transduction mechanism and its regulation in α7 receptors.

Uncovering the Differed Susceptibility of Fusarium oxysporum (Fo32931 and FocII5) to Fungicide Phenamacril: From Computational and Experimental Perspectives

Fo32931 and FoCII5 are two subtypes of Fusarium oxysporum (Fo), a pathogenic filamentous fungus. Phenamacril (PHA), a Fusarium-specific fungicide that targets myosin I, exhibits significant hyphal growth inhibition in Fo32931 but shows weak resistance in FocII5, despite only two amino acid differences in the PHA-binding pocket of myosin I. In this study, we aim to elucidate the molecular basis for the differential sensitivity ofF. oxysporum myosin I variants (FoMyoI32931 and FoMyoIcII5) to phenamacril through computational methods and biochemical validation. The results suggest that phenamacril functions as an allosteric inhibitor for FoMyoI32931, inhibiting the large oscillation of the converter lever domain (CLD) upon ATP binding and promoting the opening of the outer cleft, further impairing protein function. PHA significantly reduced the coupling between the CLD, especially the converter, and the catalytic center, diminishing the response of the CLD to the motor domain in FoMyoI32931. From the residue mutation experiment, we found that the S418T substitution in FoMyoIcII5 is the key to the reduced phenamacril sensitivity of FocII5. According to the microscale thermophoresis (MST) assay and pocket conformation analysis, the S418T mutation disturbs the orientation of pocket residues, especially Lys537, leading to a looser pocket and reduced interaction between Lys537 and phenamacril, which lowers the binding affinity of FoMyoIcII5 for phenamacril. These findings provide deeper insights into the reasons for the lower sensitivity of FoCII5 to phenamacril from both molecular and structural perspectives and will also guide the design of novel inhibitors against resistant Fusarium spp., like FoCII5.

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.

The Nicotinic Acetylcholine Receptor and Its Pentameric Homologs: Toward an Allosteric Mechanism of Signal Transduction at the Atomic Level

The nicotinic acetylcholine receptor has served, since its biochemical identification in the 1970s, as a model of an allosteric ligand-gated ion channel mediating signal transition at the synapse. In recent years, the application of X-ray crystallography and high-resolution cryo-electron microscopy, together with molecular dynamic simulations of nicotinic receptors and homologs, have opened a new era in the understanding of channel gating by the neurotransmitter. They reveal, at atomic resolution, the diversity and flexibility of the multiple ligand-binding sites, including recently discovered allosteric modulatory sites distinct from the neurotransmitter orthosteric site, and the conformational dynamics of the activation process as a molecular switch linking these multiple sites. The model emerging from these studies paves the way for a new pharmacology based, first, upon the occurrence of an original mode of indirect allosteric modulation, distinct from a steric competition for a single and rigid binding site, and second, the design of drugs that specifically interact with privileged conformations of the receptor such as agonists, antagonists, and desensitizers. Research on nicotinic receptors is still at the forefront of understanding the mode of action of drugs on the nervous system.

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.

Mechanism of proton-powered c-ring rotation in a mitochondrial ATP synthase

Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of adenosine triphosphate (ATP). Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.

Characterizing the Molecular Mechanism of the Lethal C423D Mutation in FgMyoI: A Molecular Perspective

The lethal mutation C423D in Fusarium graminearum myosin I (FgMyoI) occurs close to the binding pocket of the allosteric inhibitor phenamacril and causes severe inhibition on mycelial growth of F. graminearum strain PH-1. Here, based on extensive Gaussian accelerated molecular dynamics simulations and wet experiments, we elucidate the underlying molecular mechanism of the abnormal functioning of the FgMyoIC423D mutant at the atomistic level. Our results suggest that the damaging mutation C423D exhibits a synergistic allosteric inhibition mechanism similar to but more robust than that of phenamacril, including effects on the active site and actin binding. Unlike phenamacril-induced closure of Switch2, the mutation results in unfolding of the N-terminal relay helix with a partially opened Switch2 and blocks the structural rearrangement of the relay/SH1 helices, impairing the proper initiation of the recovery stroke. Due to the significant influence of C423D mutation on the function of FgMyoI, designing covalent inhibitors targeting this site holds tremendous potential.

Allosteric inhibition of myosin by phenamacril: a synergistic mechanism revealed by computational and experimental approaches

BACKGROUND

Myosin plays a crucial role in cellular processes, while its dysfunction can lead to organismal malfunction. Phenamacril (PHA), a highly species-specific and non-competitive inhibitor of myosin I (FgMyoI) from Fusarium graminearum, has been identified as an effective fungicide for controlling plant diseases caused by partial Fusarium pathogens, such as wheat scab and rice bakanae. However, the molecular basis of its action is still unclear.

RESULTS

This study used multiple computational approaches first to elucidate the allosteric inhibition mechanism of FgMyoI by PHA at the atomistic level. The results indicated the increase of adenosine triphosphate (ATP) binding affinity upon PHA binding, which might impede the release of hydrolysis products. Furthermore, simulations revealed a broadened outer cleft and a significantly more flexible interface for actin binding, accompanied by a decrease in signaling transduction from the catalytic center to the actin-binding interface. These various effects might work together to disrupt the actomyosin cycle and hinder the ability of motor to generate force. Our experimental results further confirmed that PHA reduces the enzymatic activity of myosin and its binding with actin.

CONCLUSION

Therefore, our findings demonstrated that PHA might suppress the function of myosin through a synergistic mechanism, providing new insights into myosin allostery and offering new avenues for drug/fungicide discovery targeting myosin. © 2023 Society of Chemical Industry.

Binding pocket dynamics along the recovery stroke of human β-cardiac myosin

The druggability of small-molecule binding sites can be significantly affected by protein motions and conformational changes. Ligand binding, protein dynamics and protein function have been shown to be closely interconnected in myosins. The breakthrough discovery of omecamtiv mecarbil (OM) has led to an increased interest in small molecules that can target myosin and modulate its function for therapeutic purposes (myosin modulators). In this work, we use a combination of computational methods, including steered molecular dynamics, umbrella sampling and binding pocket tracking tools, to follow the evolution of the OM binding site during the recovery stroke transition of human β-cardiac myosin. We found that steering two internal coordinates of the motor domain can recapture the main features of the transition and in particular the rearrangements of the binding site, which shows significant changes in size, shape and composition. Possible intermediate conformations were also identified, in remarkable agreement with experimental findings. The differences in the binding site properties observed along the transition can be exploited for the future development of conformation-selective myosin modulators.

Actomyosin Complex

Formation of cross-bridges between actin and myosin occurs ubiquitously in eukaryotic cells and mediates muscle contraction, intracellular cargo transport, and cytoskeletal remodeling. Myosin motors repeatedly bind to and dissociate from actin filaments in a cycle that transduces the chemical energy from ATP hydrolysis into mechanical force generation. While the general layout of surface elements within the actin-binding interface is conserved among myosin classes, sequence divergence within these motifs alters the specific contacts involved in the actomyosin interaction as well as the kinetics of mechanochemical cycle phases. Additionally, diverse lever arm structures influence the motility and force production of myosin molecules during their actin interactions. The structural differences generated by myosin's molecular evolution have fine-tuned the kinetics of its isoforms and adapted them for their individual cellular roles. In this chapter, we will characterize the structural and biochemical basis of the actin-myosin interaction and explain its relationship with myosin's cellular roles, with emphasis on the structural variation among myosin isoforms that enables their functional specialization. We will also discuss the impact of accessory proteins, such as the troponin-tropomyosin complex and myosin-binding protein C, on the formation and regulation of actomyosin cross-bridges.

High-resolution structures of the actomyosin-V complex in three nucleotide states provide insights into the force generation mechanism

The molecular motor myosin undergoes a series of major structural transitions during its force-producing motor cycle. The underlying mechanism and its coupling to ATP hydrolysis and actin binding are only partially understood, mostly due to sparse structural data on actin-bound states of myosin. Here, we report 26 high-resolution cryo-EM structures of the actomyosin-V complex in the strong-ADP, rigor, and a previously unseen post-rigor transition state that binds the ATP analog AppNHp. The structures reveal a high flexibility of myosin in each state and provide valuable insights into the structural transitions of myosin-V upon ADP release and binding of AppNHp, as well as the actomyosin interface. In addition, they show how myosin is able to specifically alter the structure of F-actin.

An Asymmetric Mechanism in a Symmetric Molecular Machine

The design of molecular architectures exhibiting functional motions is a promising area for disruptive technological development. Toward this goal, rotaxanes and catenanes, which undergo relative motions of their subunits in response to external stimuli, are prime candidates. Here, we report on the computational analysis of the contraction/extension of a bistable [c2]daisy chain rotaxane. Using free-energy calculations and transition path optimizations, we explore the free-energy landscape governing the functional motions of a prototypical molecular machine with atomic resolution. The calculations reveal a sequential mechanism in which the asynchronous gliding of each ring is preferred over the concerted movement. Analysis of the underlying free-energy surface indicates that the formation of partially rearranged intermediates entails crossing of much smaller barriers. Our findings illustrate an important design principle for molecular machines, namely that efficient exploitation of thermal fluctuations may be realized by breaking down the large-scale functional motions into smaller steps.

Step-Wise Hydration of Magnesium by Four Water Molecules Precedes Phosphate Release in a Myosin Motor

Molecular motors, such as myosin, kinesin, and dynein, convert the energy released by the hydrolysis of ATP into mechanical work, thus allowing them to undergo directional motion on cytoskeletal tracks. A pivotal step in the chemomechanical transduction in myosin motors occurs after they bind to the actin filament, which triggers the release of phosphate (P_i, product of ATP hydrolysis) and the rotation of the lever arm. Here, we investigate the mechanism of phosphate release in myosin VI using extensive molecular dynamics simulations involving multiple trajectories of several μs. Because the escape of phosphate is expected to occur on time-scales on the order of milliseconds or more in myosin VI, we observed P_i release only if the trajectories were initiated with a rotated phosphate inside the nucleotide binding pocket. We discovered that although P_i populates the traditional "back door" route, phosphate exits through various other gateways, thus establishing the heterogeneity in the escape routes. Remarkably, we observed that the release of phosphate is preceded by a stepwise hydration of the ADP-bound magnesium ion. The release of the anion occurred only after four water molecules hydrated the cation (Mg²⁺). By performing comparative structural analyses, we show that hydration of magnesium is the key step in the phosphate release in a number of ATPases and GTPases. Nature may have evolved hydration of Mg²⁺ as a general molecular switch for P_i release, which is a universal step in the catalytic cycle of many machines that share little sequence or structural similarity.

Unraveling a Force-Generating Allosteric Pathway of Actomyosin Communication Associated with ADP and Pi Release

The actomyosin system generates mechanical work with the execution of the power stroke, an ATP-driven, two-step rotational swing of the myosin-neck that occurs post ATP hydrolysis during the transition from weakly to strongly actin-bound myosin states concomitant with P_i release and prior to ADP dissociation. The activating role of actin on product release and force generation is well documented; however, the communication paths associated with weak-to-strong transitions are poorly characterized. With the aid of mutant analyses based on kinetic investigations and simulations, we identified the W-helix as an important hub coupling the structural changes of switch elements during ATP hydrolysis to temporally controlled interactions with actin that are passed to the central transducer and converter. Disturbing the W-helix/transducer pathway increased actin-activated ATP turnover and reduced motor performance as a consequence of prolonged duration of the strongly actin-attached states. Actin-triggered P_i release was accelerated, while ADP release considerably decelerated, both limiting maximum ATPase, thus transforming myosin-2 into a high-duty-ratio motor. This kinetic signature of the mutant allowed us to define the fractional occupancies of intermediate states during the ATPase cycle providing evidence that myosin populates a cleft-closure state of strong actin interaction during the weak-to-strong transition with bound hydrolysis products before accomplishing the power stroke.

Structural and Computational Insights into a Blebbistatin-Bound Myosin•ADP Complex with Characteristics of an ADP-Release Conformation along the Two-Step Myosin Power Stoke

The motor protein myosin drives a wide range of cellular and muscular functions by generating directed movement and force, fueled through adenosine triphosphate (ATP) hydrolysis. Release of the hydrolysis product adenosine diphosphate (ADP) is a fundamental and regulatory process during force production. However, details about the molecular mechanism accompanying ADP release are scarce due to the lack of representative structures. Here we solved a novel blebbistatin-bound myosin conformation with critical structural elements in positions between the myosin pre-power stroke and rigor states. ADP in this structure is repositioned towards the surface by the phosphate-sensing P-loop, and stabilized in a partially unbound conformation via a salt-bridge between Arg131 and Glu187. A 5 Å rotation separates the mechanical converter in this conformation from the rigor position. The crystallized myosin structure thus resembles a conformation towards the end of the two-step power stroke, associated with ADP release. Computationally reconstructing ADP release from myosin by means of molecular dynamics simulations further supported the existence of an equivalent conformation along the power stroke that shows the same major characteristics in the myosin motor domain as the resolved blebbistatin-bound myosin-II·ADP crystal structure, and identified a communication hub centered on Arg232 that mediates chemomechanical energy transduction.

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Molecular dynamics (MD) simulation is a rigorous theoretical tool that when used efficiently could provide reliable answers to questions pertaining to the structure-function relationship of proteins. Data collated from protein dynamics can be translated into useful statistics that can be exploited to sieve thermodynamics and kinetics crucial for the elucidation of mechanisms responsible for the modulation of biological processes such as protein-ligand binding and protein-protein association. Continuous modernization of simulation tools enables accurate prediction and characterization of the aforementioned mechanisms and these qualities are highly beneficial for the expedition of drug development when effectively applied to structure-based drug design (SBDD). In this review, current all-atom MD simulation methods, with focus on enhanced sampling techniques, utilized to examine protein structure, dynamics, and functions are discussed. This review will pivot around computer calculations of protein-ligand and protein-protein systems with applications to SBDD. In addition, we will also be highlighting limitations faced by current simulation tools as well as the improvements that have been made to ameliorate their efficiency.

Conformational distributions of isolated myosin motor domains encode their mechanochemical properties

Myosin motor domains perform an extraordinary diversity of biological functions despite sharing a common mechanochemical cycle. Motors are adapted to their function, in part, by tuning the thermodynamics and kinetics of steps in this cycle. However, it remains unclear how sequence encodes these differences, since biochemically distinct motors often have nearly indistinguishable crystal structures. We hypothesized that sequences produce distinct biochemical phenotypes by modulating the relative probabilities of an ensemble of conformations primed for different functional roles. To test this hypothesis, we modeled the distribution of conformations for 12 myosin motor domains by building Markov state models (MSMs) from an unprecedented two milliseconds of all-atom, explicit-solvent molecular dynamics simulations. Comparing motors reveals shifts in the balance between nucleotide-favorable and nucleotide-unfavorable P-loop conformations that predict experimentally measured duty ratios and ADP release rates better than sequence or individual structures. This result demonstrates the power of an ensemble perspective for interrogating sequence-function relationships.

Force Generation by Myosin Motors: A Structural Perspective

Generating force and movement is essential for the functions of cells and organisms. A variety of molecular motors that can move on tracks within cells have evolved to serve this role. How these motors interact with their tracks and how that, in turn, leads to the generation of force and movement is key to understanding the cellular roles that these motor-track systems serve. This review is focused on the best understood of these systems, which is the molecular motor myosin that moves on tracks of filamentous (F-) actin. The review highlights both the progress and the limits of our current understanding of how force generation can be controlled by F-actin-myosin interactions. What has emerged are insights they may serve as a framework for understanding the design principles of a number of types of molecular motors and their interactions with their tracks.

Structural basis for power stroke vs. Brownian ratchet mechanisms of motor proteins

Two mechanisms have been proposed for the function of motor proteins: The power stroke and the Brownian ratchet. The former refers to generation of a large downhill free energy gradient over which the motor protein moves nearly irreversibly in making a step, whereas the latter refers to biasing or rectifying the diffusive motion of the motor. Both mechanisms require input of free energy, which generally involves the processing of an ATP (adenosine 5'-triphosphate) molecule. Recent advances in experiments that reveal the details of the stepping motion of motor proteins, together with computer simulations of atomistic structures, have provided greater insights into the mechanisms. Here, we compare the various models of the power stroke and the Brownian ratchet that have been proposed. The 2 mechanisms are not mutually exclusive, and various motor proteins employ them to different extents to perform their biological function. As examples, we discuss linear motor proteins Kinesin-1 and myosin-V, and the rotary motor F₁-ATPase, all of which involve a power stroke as the essential element of their stepping mechanism.