Advances in quantum simulations of ATPase catalysis in the myosin motor

Myosin-Catalyzed ATP Hydrolysis in the Presence of Disease-Causing Mutations: Mavacamten as a Way to Repair Mechanism

Hypertrophic cardiomyopathy is one of the most common forms of genetic cardiomyopathy. Mavacamten is a first-in-class myosin modulator that was identified via activity screening on the wild type, and it is FDA-approved for the treatment of obstructive hypertrophic cardiomyopathy (HCM). The drug selectively binds to the cardiac β-myosin, inhibiting myosin function to decrease cardiac contractility. Though the drug is thought to affect multiple steps of the myosin cross-bridge cycle, its detailed mechanism of action is still under investigation. Individual steps in the overall cross-bridge cycle must be queried to elucidate the full mechanism of action. In this study, we utilize the rare-event method of transition path sampling to generate reactive trajectories to gain insights into the action of the drug on the dynamics and rate of the ATP hydrolysis step for human cardiac β-myosin. We study three known HCM causative myosin mutations: R453C, P710R, and R712L to observe the effect of the drug on the alterations caused by these mutations in the chemical step. Since the crystal structure of the drug-bound myosin was not available at the time of this work, we created a model of the drug-bound system utilizing a molecular docking approach. We find a significant effect of the drug in one case, where the actual mechanism of the reaction is altered from the wild type by mutation. The drug restores both the rate of hydrolysis to the wildtype level and the mechanism of the reaction. This is a way to check the effect of the drug on untested mutations.

Mechanism of small molecule inhibition of Plasmodium falciparum myosin A informs antimalarial drug design

Malaria results in more than 500,000 deaths per year and the causative Plasmodium parasites continue to develop resistance to all known agents, including different antimalarial combinations. The class XIV myosin motor PfMyoA is part of a core macromolecular complex called the glideosome, essential for Plasmodium parasite mobility and therefore an attractive drug target. Here, we characterize the interaction of a small molecule (KNX-002) with PfMyoA. KNX-002 inhibits PfMyoA ATPase activity in vitro and blocks asexual blood stage growth of merozoites, one of three motile Plasmodium life-cycle stages. Combining biochemical assays and X-ray crystallography, we demonstrate that KNX-002 inhibits PfMyoA using a previously undescribed binding mode, sequestering it in a post-rigor state detached from actin. KNX-002 binding prevents efficient ATP hydrolysis and priming of the lever arm, thus inhibiting motor activity. This small-molecule inhibitor of PfMyoA paves the way for the development of alternative antimalarial treatments.

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.

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.

Alpha and beta myosin isoforms and human atrial and ventricular contraction

Human atrial and ventricular contractions have distinct mechanical characteristics including speed of contraction, volume of blood delivered and the range of pressure generated. Notably, the ventricle expresses predominantly β-cardiac myosin while the atrium expresses mostly the α-isoform. In recent years exploration of the properties of pure α- & β-myosin isoforms have been possible in solution, in isolated myocytes and myofibrils. This allows us to consider the extent to which the atrial vs ventricular mechanical characteristics are defined by the myosin isoform expressed, and how the isoform properties are matched to their physiological roles. To do this we Outline the essential feature of atrial and ventricular contraction; Explore the molecular structural and functional characteristics of the two myosin isoforms; Describe the contractile behaviour of myocytes and myofibrils expressing a single myosin isoform; Finally we outline the outstanding problems in defining the differences between the atria and ventricles. This allowed us consider what features of contraction can and cannot be ascribed to the myosin isoforms present in the atria and ventricles.

Effects of protein posttranslational modifications on meat quality: A review

Meat quality plays an important role in the purchase decision of consumers, affecting producers and retailers. The formation mechanisms determining meat quality are intricate, as several endogenous and exogenous factors contribute during antemortem and postmortem periods. Abundant research has been performed on meat quality; however, unexpected variation in meat quality remains an issue in the meat industry. Protein posttranslational modifications (PTMs) regulate structures and functions of proteins in living tissues, and recent reports confirmed their importance in meat quality. The objective of this review was to provide a summary of the research on the effects of PTMs on meat quality. The effects of four common PTMs, namely, protein phosphorylation, acetylation, S-nitrosylation, and ubiquitination, on meat quality were discussed, with emphasis on the effects of protein phosphorylation on meat tenderness, color, and water holding capacity. The mechanisms and factors that may affect the function of protein phosphorylation are also discussed. The current research confirms that meat quality traits are regulated by multiple PTMs. Cross talk between different PTMs and interactions of PTMs with postmortem biochemical processes need to be explored to improve our understanding on factors affecting meat quality.

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.

Molecular Mechanism of ATP Hydrolysis in an ABC Transporter

Hydrolysis of nucleoside triphosphate (NTP) plays a key role for the function of many biomolecular systems. However, the chemistry of the catalytic reaction in terms of an atomic-level understanding of the structural, dynamic, and free energy changes associated with it often remains unknown. Here, we report the molecular mechanism of adenosine triphosphate (ATP) hydrolysis in the ATP-binding cassette (ABC) transporter BtuCD-F. Free energy profiles obtained from hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations show that the hydrolysis reaction proceeds in a stepwise manner. First, nucleophilic attack of an activated lytic water molecule at the ATP γ-phosphate yields ADP + HPO₄ ^2- as intermediate product. A conserved glutamate that is located very close to the γ-phosphate transiently accepts a proton and thus acts as catalytic base. In the second step, the proton is transferred back from the catalytic base to the γ-phosphate, yielding ADP + H₂PO₄ ^-. These two chemical reaction steps are followed by rearrangements of the hydrogen bond network and the coordination of the Mg²⁺ ion. The rate constant estimated from the computed free energy barriers is in very good agreement with experiments. The overall free energy change of the reaction is close to zero, suggesting that phosphate bond cleavage itself does not provide a power stroke for conformational changes. Instead, ATP binding is essential for tight dimerization of the nucleotide-binding domains and the transition of the transmembrane domains from inward- to outward-facing, whereas ATP hydrolysis resets the conformational cycle. The mechanism is likely relevant for all ABC transporters and might have implications also for other NTPases, as many residues involved in nucleotide binding and hydrolysis are strictly conserved.

The hydrolysis of 6-phosphogluconolactone in the second step of pentose phosphate pathway occurs via a two-water mechanism

Hydrolysis reaction marks the basis of life yet the mechanism of this crucial biochemical reaction is not completely understood. We recently reported the mechanisms of hydrolysis of nucleoside triphosphate and phosphate monoester. These two reactions hydrolyze P-O-P and P-O-C linkages, respectively. Here, we present the mechanism of hydrolysis of δ-6-phosphogluconolactone, which is an important precursor in the second step of the pentose phosphate pathway. Its hydrolysis requires the cleavage of C-O-C linkage and its mechanism is hitherto unknown. We report three mechanisms of hydrolysis of δ-6-phosphogluconolactone based on density functional computations. In the energetically most favorable mechanism, two water molecules participate in the hydrolysis reaction and the mechanism is sequential, i.e., activation of the attacking water molecule (OH bond breaking) precedes that of the cleavage of the CO bond of the C-O-C linkage. The rate-limiting energy barrier of this mechanism is comparable to the reported experimental free energy barrier. This mechanism has similarities with the mechanism of triphosphate hydrolysis and that of hydrolytic cleavage of DNA in EcoRV enzyme. This two-water sequential hydrolysis mechanism could be the unified mechanism required for the hydrolysis of other hydrolysable species in living cells.

Medicinal Chemistry and Use of Myosin II Inhibitor ( S)-Blebbistatin and Its Derivatives

( S)-Blebbistatin, a chiral tetrahydropyrroloquinolinone, is a widely used and well-characterized ATPase inhibitor selective for myosin II. The central role of myosin II in many normal and pathological biological processes has been revealed with the aid of this small molecule. The first part of this manuscript provides a summary of myosin II and ( S)-blebbistatin literature from a medicinal chemist's perspective. The second part of this perspective deals with the physicochemical deficiencies that trouble the use of ( S)-blebbistatin in advanced biological settings: low potency and solubility, fluorescence interference, (photo)toxicity, and stability issues. A large toolbox of analogues has been developed in which particular shortcomings have been addressed. This perspective provides a necessary overview of these developments and presents guidelines for selecting the best available analogue for a given application. As the unmet need for high-potency analogues remains, we also propose starting points for medicinal chemists in search of nanomolar myosin II inhibitors.

Comparing the catalytic strategy of ATP hydrolysis in biomolecular motors

ATP-driven biomolecular motors utilize the chemical energy obtained from the ATP hydrolysis to perform vital tasks in living cells. Understanding the mechanism of enzyme-catalyzed ATP hydrolysis reaction has substantially progressed lately thanks to combined quantum/classical molecular mechanics (QM/MM) simulations. Here, we present a comparative summary of the most recent QM/MM results for myosin, kinesin and F1-ATPase motors. These completely different motors achieve the acceleration of ATP hydrolysis through a very similar catalytic mechanism. ATP hydrolysis has high activation energy because it involves the breaking of two strong bonds, namely the Pγ-Oβγ bond of ATP and the H-O bond of lytic water. The key to the four-fold decrease in the activation barrier by the three enzymes is that the breaking of the Pγ-Oβγ bond precedes the deprotonation of the lytic water molecule, generating a metaphosphate hydrate complex. The resulting singly charged trigonal planar PγO3(-) metaphosphate is a better electrophilic target for attack by an OaH(-) hydroxyl group. The formation of this OaH(-) is promoted by a strong polarization of the lytic water: in all three proteins, this water is forming a hydrogen-bond with a backbone carbonyl group and interacts with the carboxylate group of glutamate (either directly or via an intercalated water molecule). This favors the shedding of one proton by the attacking water. The abstracted proton is transferred to the γ-phosphate via various proton wires, resulting in a H2PγO4(-)/ADP(3-) product state. This catalytic strategy is so effective that most other nucleotide hydrolyzing enzymes adopt a similar approach, as suggested by their very similar triphosphate binding sites.

Effect of protonation on the mechanism of phosphate monoester hydrolysis and comparison with the hydrolysis of nucleoside triphosphate in biomolecular motors

Hydrolysis of phosphate groups is a crucial reaction in living cells. It involves the breaking of two strong bonds, i.e. the O_aH bond of the attacking water molecule, and the PO_l bond of the substrate (O_a and O_l stand for attacking and leaving oxygen atoms). Mechanism of the hydrolysis reaction can proceed either by a concurrent or a sequential mechanism. In the concurrent mechanism, the breaking of O_aH and PO_l bonds occurs simultaneously, whereas in the sequential mechanism, the O_aH and PO_l bonds break at different stages of the reaction. To understand how protonation affects the mechanism of hydrolysis of phosphate monoester, we have studied the mechanism of hydrolysis of protonated and deprotonated phosphate monoester at M06-2X/6-311+G**//M06-2X/6-31+G*+ZPE level of theory (where ZPE stands for zero point energy). Our calculations show that in both protonated and deprotonated cases, the breaking of the water O_aH bond occurs before the breaking of the PO_l bond. Because the two events are not separated by a stable intermediate, the mechanism can be categorized as semi-concurrent. The overall energy barrier is 41kcalmol^-1 in the unprotonated case. Most (5/6th) of this is due to the initial breaking of the water O_aH bond. This component is lowered from 34 to 25kcalmol^-1 by adding one proton to the phosphate. The rest of the overall energy barrier comes from the subsequent breaking of the PO_l bond and is not sensitive to protonation. This is consistent with previous findings about the effect of triphosphate protonation on the hydrolysis, where the equivalent protonation (on the γ-phosphate) was seen to lower the barrier of breaking the water O_aH bond and to have little effect on the PO_l bond breaking. Hydrolysis pathways of phosphate monoester with initial breaking of the PO_l bond could not be found here. This is because the leaving group in phosphate monoester cannot be protonated, unlike in triphosphate hydrolysis, where protonation of the β- and γ-phosphates had been shown to promote a mechanism where the PO_l bond breaks before the O_aH bond does. We also point out that the charge shift due to PO_l bond breaking during sequential ATP hydrolysis in bio-molecular motors onsets the week unbinding of hydrolysis product that finally leads to the product release during 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.

Review: The ATPase mechanism of myosin and actomyosin

Myosins are a large family of molecular motors that use the common P-loop, Switch 1 and Switch 2 nucleotide binding motifs to recognize ATP, to create a catalytic site than can efficiently hydrolyze ATP and to communicate the state of the nucleotide pocket to other allosteric binding sites on myosin. The energy of ATP hydrolysis is used to do work against an external load. In this short review I will outline current thinking on the mechanism of ATP hydrolysis and how the energy of ATP hydrolysis is coupled to a series of protein conformational changes that allow a myosin, with the cytoskeleton track actin, to operate as a molecular motor of distinct types; fast movers, processive motors or strain sensors. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 483-491, 2016.

Simulating the reactions of substituted pyridinio-N-phosphonates with pyridine as a model for biological phosphoryl transfer

This work provides a comprehensive model for non-enzymatic phosphoryl transfer, as a baseline for understanding biological phosphoryl transfer reactions.

ATP-dependent interplay between local and global conformational changes in the myosin motor

The ATPase active site of myosin is located at the core of the motor head. During the Lymn-Taylor actomyosin contractile cycle, small conformational changes in the active site upon ATP binding, ATP hydrolysis and ADP/P_i release are accompanied by large conformational transitions of the motor domains, such as opening and closing of the actin binding cleft and the movement of lever arm. Here, our previous computational studies of myosin are summarized in a comprehensive model at the level of atomic detail. Molecular movies show how the successive domain motions during the ATP induced actin dissociation and the recovery stroke are coupled with the precise positioning of the key catalytic groups in the active site. This leads to a precise timing of the activation of the ATPase function: it allows ATP hydrolysis only after unbinding from actin and the priming of the lever arm, both pre-requisites for an efficient functioning of the motor during the subsequent power stroke. These coupling mechanisms constitute essential principles of every myosin motor, of which the ATP-site can be seen as the central allosteric control unit. © 2016 Wiley Periodicals, Inc.

Effects of protonation on the hydrolysis of triphosphate in vacuum and the implications for catalysis by nucleotide hydrolyzing enzymes

BACKGROUND

Nucleoside triphosphate (NTP) hydrolysis is a key reaction in biology. It involves breaking two very stable bonds (one P-O bond and one O-H bond of water), in either a concurrent or a sequential way. Here, we systematically examine how protonation of the triphosphate affects the mechanism of hydrolysis.

RESULTS

The hydrolysis reaction of methyl triphosphate in vacuum is computed with protons in various numbers and position on the three phosphate groups. Protonation is seen to have a strong catalytic effect, with the reaction mechanism depending highly on the protonation pattern.

CONCLUSION

This dependence is apparently complicated, but is shown to obey a well-defined set of rules: Protonation of the α- and β-phosphate groups favors a sequential hydrolysis mechanism, whereas γ-protonation favors a concurrent mechanism, the two effects competing with each other in cases of simultaneous protonation. The rate-limiting step is always the breakup of the water molecule while it attacks the γ-phosphorus, and its barrier is lowered by γ-protonation. This step has significantly lower barriers in the sequential reactions, because the dissociated γ-metaphosphate intermediate (PγO3-) is a much better target for water attack than the un-dissociated γ-phosphate (-PγO42-). The simple chemical logic behind these rules helps to better understand the catalytic strategy used by NTPase enzymes, as illustrated here for the catalytic pocket of myosin. A set of rules was determined that describes how protonating the phosphate groups affects the hydrolysis mechanism of methyl triphosphate: Protonation of the α- and/or β- phosphate groups promotes a sequential mechanism in which P-O bond breaking precedes the breakup of the attacking water, whereas protonation of the γ-phosphate promotes a concurrent mechanism and lowers the rate-limiting barrier of water breakup. The role played by individual protein residues in the catalytic pocket of triphosphate hydrolysing enzymes can be assigned accordingly.