Myosins adjust their power outputs in response to mechanical loads in an isoform-dependent manner, resulting in their ability to dynamically adapt to a range of motile challenges. Here, we reveal the structural basis for force-sensing based on near-atomic resolution structures of one rigor and two ADP-bound states of myosin-IB (myo1b) bound to actin, determined by cryo-electron microscopy. The two ADP-bound states are separated by a 25° rotation of the lever. The lever of the first ADP state is rotated toward the pointed end of the actin filament and forms a previously unidentified interface with the N-terminal subdomain, which constitutes the upper half of the nucleotide-binding cleft. This pointed-end orientation of the lever blocks ADP release by preventing the N-terminal subdomain from the pivoting required to open the nucleotide binding site, thus revealing how myo1b is inhibited by mechanical loads that restrain lever rotation. The lever of the second ADP state adopts a rigor-like orientation, stabilized by class-specific elements of myo1b. We identify a role for this conformation as an intermediate in the ADP release pathway. Moreover, comparison of our structures with other myosins reveals structural diversity in the actomyosin binding site, and we reveal the high-resolution structure of actin-bound phalloidin, a potent stabilizer of filamentous actin. These results provide a framework to understand the spectrum of force-sensing capacities among the myosin superfamily.

containing (in mM): 30 KCl, 0.1 CaCl2, 0.5 MgCl2, 2 DTT, 20 Imidazole-HCl (pH 7.2) with a slightly modified Fiske and Subbarow method. The ATPase reaction was started after 10 min of incubation at 25 C by adding 1 mM of ATP and was terminated after 10 min by adding 2.5% trichloroacetic acid. The amount of Pi liberated was evaluated colorimetrically at 595nm. Results & discussion: The Mg2þ-ATPase activity of R403Q myosin (0.186 mmol/mg/min) was 2.8 times lower than the activity of WT myosin (0.53 mmol/mg/min) (p<0.05). R403Q HMM fragments also reported a drop in Mg2þ-ATPase (0.086 mmol/mg/min) compared to the WT HMM (0.186 mmol/mg/min) activity (p<0.05). However, there were no significant differences between the Mg2þ-ATPase activity of R403Q and WT myofibrils. The result of this study suggests that the R403Q mutation in the heart leads to a reduced contraction of the diaphragm.

Parallel replica exchange sampling is an extended ensemble technique often used to accelerate the exploration of the conformational ensemble of atomistic molecular simulations of chemical systems. Inter-process communication and coordination requirements have historically discouraged the deployment of replica exchange on distributed and heterogeneous resources. Here we describe the architecture of a software (named ASyncRE) for performing asynchronous replica exchange molecular simulations on volunteered computing grids and heterogeneous high performance clusters. The asynchronous replica exchange algorithm on which the software is based avoids centralized synchronization steps and the need for direct communication between remote processes. It allows molecular dynamics threads to progress at different rates and enables parameter exchanges among arbitrary sets of replicas independently from other replicas. ASyncRE is written in Python following a modular design conducive to extensions to various replica exchange schemes and molecular dynamics engines. Applications of the software for the modeling of association equilibria of supramolecular and macromolecular complexes on BOINC campus computational grids and on the CPU/MIC heterogeneous hardware of the XSEDE Stampede supercomputer are illustrated. They show the ability of ASyncRE to utilize large grids of desktop computers running the Windows, MacOS, and/or Linux operating systems as well as collections of high performance heterogeneous hardware devices.

Force field accuracy is still one of the “stalemates” in biomolecular modeling. Model systems with high quality experimental data are valuable instruments for the validation and improvement of effective potentials. With respect to protein–ligand binding, organic host–guest complexes have long served as models for both experimental and computational studies because of the abundance of binding affinity data available for such systems. Binding affinity data collected for cyclodextrin (CD) inclusion complexes, a popular model for molecular recognition, is potentially a more reliable resource for tuning energy parameters than hydration free energy measurements. Convergence of binding free energy calculations on CD host–guest systems can also be obtained rapidly, thus offering the opportunity to assess the robustness of these parameters. In this work, we demonstrate how implicit solvent parameters can be developed using binding affinity experimental data and the binding energy distribution analysis method (BEDAM) and validated using the Grid Inhomogeneous Solvation Theory analysis. These new solvation parameters were used to study protein–ligand binding in two drug targets against the HIV-1 virus and improved the agreement between the calculated and the experimental binding affinities. This work illustrates how benchmark sets of high quality experimental binding affinity data and physics-based binding free energy models can be used to evaluate and optimize force fields for protein–ligand systems.

Cancer cells bearing distinct KRAS mutations exhibit variable sensitivity to SHP2 inhibitors (SHP2i). Here we show that cells harboring KRAS Q61H are uniquely resistant to SHP2i, and investigate the underlying mechanisms using biophysics, molecular dynamics, and cell-based approaches. Q61H mutation impairs intrinsic and GAP-mediated GTP hydrolysis, and impedes activation by SOS1, but does not alter tyrosyl phosphorylation. Wild-type and Q61H-mutant KRAS are both phosphorylated by Src on Tyr32 and Tyr64 and dephosphorylated by SHP2, however, SHP2i does not reduce ERK phosphorylation in KRAS Q61H cells. Phosphorylation of wild-type and Gly12-mutant KRAS, which are associated with sensitivity to SHP2i, confers resistance to regulation by GAP and GEF activities and impairs binding to RAF, whereas the near-complete GAP/GEF-resistance of KRAS Q61H remains unaltered, and high-affinity RAF interaction is retained. SHP2 can stimulate KRAS signaling by modulating GEF/GAP activities and dephosphorylating KRAS, processes that fail to regulate signaling of the Q61H mutant.