Catalytic strategy used by the myosin motor to hydrolyze ATP

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.

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.

Polyphosphate Nanoparticles: Balancing Energy Requirements in Tissue Regeneration Processes

Nanoparticles of a particular, evolutionarily old inorganic polymer found across the biological kingdoms have attracted increasing interest in recent years not only because of their crucial role in metabolism but also their potential medical applicability: it is inorganic polyphosphate (polyP). This ubiquitous linear polymer is composed of 10-1000 phosphate residues linked by high-energy anhydride bonds. PolyP causes induction of gene activity, provides phosphate for bone mineralization, and serves as an energy supplier through enzymatic cleavage of its acid anhydride bonds and subsequent ATP formation. The biomedical breakthrough of polyP came with the development of a successful fabrication process, in depot form, as Ca- or Mg-polyP nanoparticles, or as the directly effective polymer, as soluble Na-polyP, for regenerative repair and healing processes, especially in tissue areas with insufficient blood supply. Physiologically, the platelets are the main vehicles for polyP nanoparticles in the circulating blood. To be biomedically active, these particles undergo coacervation. This review provides an overview of the properties of polyP and polyP nanoparticles for applications in the regeneration and repair of bone, cartilage, and skin. In addition to studies on animal models, the first successful proof-of-concept studies on humans for the healing of chronic wounds are outlined.

Does muscle-type myosin have ADPase activity?

Adenosine diphosphate (ADP) is a nucleotide that is structurally very similar to ATP but lacks one of the two high-energy bonds due to hydrolysis. In muscle studies, ADP is usually considered exclusively as a product formed during myosin cross-bridge cycling and is not otherwise involved in this process. In our study, we question the widely held view of ADP as a final product formed during muscle contraction. Using biophysical and biochemical methods, we managed to show that ADP can act as a substrate for myosins in at least three types of muscles: smooth and striated adductor muscles of bivalves (Mytilidae and Pectinidae), and also vertebrate skeletal muscles. According to our data, the differences in the effect of ATP and ADP on the optical, biochemical, and structural properties of actomyosins are exclusively quantitative. We explain the previous ideas about ADP as a compound capable of inhibiting the ATPase activity of actomyosin by the ability of ATP and ADP to depolymerize the polymeric myosin when the concentration in the medium reaches more than 0.3 mM.

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.

Structures and mechanisms of actin ATP hydrolysis

A variety of cellular functions are driven by actin, which undergoes cyclic transitions between the monomeric G-form and the filamentous F-form. To gain insights into actin dynamics, the mechanism by which the energy is supplied by the ATP hydrolysis reaction in the F-form actin must be elucidated. This has been hampered by the lack of actin filament structures at atomic resolutions. Here, we have crystallized actin molecules trapped in the F-form without forming filaments, and based upon these structures we determined the reaction path by quantum mechanics calculations. The results are consistent with previous biochemical data. Remarkably, the hydrolysis reaction mechanism is essentially identical to those of motor proteins, while the process of Pi release is distinct.

The major cytoskeleton protein actin undergoes cyclic transitions between the monomeric G-form and the filamentous F-form, which drive organelle transport and cell motility. This mechanical work is driven by the ATPase activity at the catalytic site in the F-form. For deeper understanding of the actin cellular functions, the reaction mechanism must be elucidated. Here, we show that a single actin molecule is trapped in the F-form by fragmin domain-1 binding and present their crystal structures in the ATP analog-, ADP-Pi-, and ADP-bound forms, at 1.15-Å resolutions. The G-to-F conformational transition shifts the side chains of Gln137 and His161, which relocate four water molecules including W1 (attacking water) and W2 (helping water) to facilitate the hydrolysis. By applying quantum mechanics/molecular mechanics calculations to the structures, we have revealed a consistent and comprehensive reaction path of ATP hydrolysis by the F-form actin. The reaction path consists of four steps: 1) W1 and W2 rotations; 2) P^G–O^3B bond cleavage; 3) four concomitant events: W1–PO₃⁻ formation, OH⁻ and proton cleavage, nucleophilic attack by the OH⁻ against P^G, and the abstracted proton transfer; and 4) proton relocation that stabilizes the ADP-Pi–bound F-form actin. The mechanism explains the slow rate of ATP hydrolysis by actin and the irreversibility of the hydrolysis reaction. While the catalytic strategy of actin ATP hydrolysis is essentially the same as those of motor proteins like myosin, the process after the hydrolysis is distinct and discussed in terms of Pi release, F-form destabilization, and global conformational changes.

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.

Phosphorus Compounds of Natural Origin: Prebiotic, Stereochemistry, Application

Organophosphorus compounds play a vital role as nucleic acids, nucleotide coenzymes, metabolic intermediates and are involved in many biochemical processes. They are part of DNA, RNA, ATP and a number of important biological elements of living organisms. Synthetic compounds of this class have found practical application as agrochemicals, pharmaceuticals, bioregulators, and othrs. In recent years, a large number of phosphorus compounds containing P-O, P-N, P-C bonds have been isolated from natural sources. Many of them have shown interesting biological properties and have become the objects of intensive scientific research. Most of these compounds contain asymmetric centers, the absolute configurations of which have a significant effect on the biological properties of the products of their transformations. This area of research on natural phosphorus compounds is still little-studied, that prompted us to analyze and discuss it in our review. Moreover natural organophosphorus compounds represent interesting models for the development of new biologically active compounds, and a number of promising drugs and agrochemicals have already been obtained on their basis. The review also discusses the history of the development of ideas about the role of organophosphorus compounds and stereochemistry in the origin of life on Earth, starting from the prebiotic period, that allows us in a new way to consider this most important problem of fundamental science.

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.

miR-29c&b2 encourage extramedullary infiltration resulting in the poor prognosis of acute myeloid leukemia

Extramedullary infiltration (EMI), as a concomitant symptom of acute myeloid leukemia (AML), is associated with low complete remission and poor prognosis in AML. However, the mechanism of EMI remains indistinct. Clinical trials showed that increased miR-29s were associated with a poor overall survival in AML [14]. Nevertheless, they were proved to work as tumor suppressor genes by encouraging apoptosis and inhibiting proliferation in vitro. These contradictory results led us to the hypothesis that miR-29s may play a notable role in the prognosis of AML rather than leukemogenesis. Thus, we explored the specimens of AML patients and addressed this issue into miR-29c&b2 knockout mice. As a result, a poor overall survival and invasive blast cells were observed in high miR-29c&b2-expression patients, and the wildtype mice presented a shorter survival with heavier leukemia infiltration in extramedullary organs. Subsequently, we found that the miR-29c&b2 inside leukemia cells promoted EMI, but not the one in the microenvironment. The analysis of signal pathway revealed that miR-29c&b2 could target HMG-box transcription factor 1 (Hbp1) directly, then reduced Hbp1 bound to the promoter of non-muscle myosin IIB (Myh10) as a transcript inhibitor. Thus, increased Myh10 encouraged the migration of leukemia cells. Accordingly, AML patients with EMI were confirmed to have high miR-29c&b2 and MYH10 with low HBP1. Therefore, we identify that miR-29c&b2 contribute to the poor prognosis of AML patients by promoting EMI, and related genes analyses are prospectively feasible in assessment of AML outcome.

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.

COVID-19: the CaMKII-like system of S protein drives membrane fusion and induces syncytial multinucleated giant cells

The SARS-CoV-2 S protein on the membrane of infected cells can promote receptor-dependent syncytia formation, relating to extensive tissue damage and lymphocyte elimination. In this case, it is challenging to obtain neutralizing antibodies and prevent them through antibodies effectively. Considering that, in the current study, structural domain search methods are adopted to analyze the SARS-CoV-2 S protein to find the fusion mechanism. The results show that after the EF-hand domain of S protein bound to calcium ions, S2 protein had CaMKII protein activities. Besides, the CaMKII_AD domain of S2 changed S2 conformation, facilitating the formation of HR1-HR2 six-helix bundles. Apart from that, the Ca²⁺-ATPase of S2 pumped calcium ions from the virus cytoplasm to help membrane fusion, while motor structures of S drove the CaATP_NAI and CaMKII_AD domains to extend to the outside and combined the viral membrane and the cell membrane, thus forming a calcium bridge. Furthermore, the phospholipid-flipping-ATPase released water, triggering lipid mixing and fusion and generating fusion pores. Then, motor structures promoted fusion pore extension, followed by the cytoplasmic contents of the virus being discharged into the cell cytoplasm. After that, the membrane of the virus slid onto the cell membrane along the flowing membrane on the gap of the three CaATP_NAI. At last, the HR1-HR2 hexamer would fall into the cytoplasm or stay on the cell membrane. Therefore, the CaMKII_like system of S protein facilitated membrane fusion for further inducing syncytial multinucleated giant cells.

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.

Inorganic Polyphosphates As Storage for and Generator of Metabolic Energy in the Extracellular Matrix

Inorganic polyphosphates (polyP) consist of linear chains of orthophosphate residues, linked by high-energy phosphoanhydride bonds. They are evolutionarily old biopolymers that are present from bacteria to man. No other molecule concentrates as much (bio)chemically usable energy as polyP. However, the function and metabolism of this long-neglected polymer are scarcely known, especially in higher eukaryotes. In recent years, interest in polyP experienced a renaissance, beginning with the discovery of polyP as phosphate source in bone mineralization. Later, two discoveries placed polyP into the focus of regenerative medicine applications. First, polyP shows morphogenetic activity, i.e., induces cell differentiation via gene induction, and, second, acts as an energy storage and donor in the extracellular space. Studies on acidocalcisomes and mitochondria provided first insights into the enzymatic basis of eukaryotic polyP formation. In addition, a concerted action of alkaline phosphatase and adenylate kinase proved crucial for ADP/ATP generation from polyP. PolyP added extracellularly to mammalian cells resulted in a 3-fold increase of ATP. The importance and mechanism of this phosphotransfer reaction for energy-consuming processes in the extracellular matrix are discussed. This review aims to give a critical overview about the formation and function of this unique polymer that is capable of storing (bio)chemically useful energy.

Mechanistic study of the ATP hydrolysis reaction in dynein motor protein

Dynein, a large and complex motor protein, harnesses energy from adenosine triphosphate (ATP) hydrolysis to regulate essential cellular activities. The ATP hydrolysis mechanism for the dynein motor is still shrouded in mystery. Herein, molecular dynamics simulations of a dynein motor disclosed that two water molecules are present close to the γ-phosphate of ATP and Glu1742 at the AAA1 site of dynein. We have proposed three possible mechanisms for the ATP hydrolysis. We divulge by using a quantum mechanics/molecular mechanics (QM/MM) study that two water molecules and Glu1742 are crucial for facilitating the ATP hydrolysis reaction in dynein. Moreover, the ATP hydrolysis step is initiated by the activation of lytic water (W1) by Glu1742 through relay proton transfers with the help of auxiliary water (W2) yielding HPO₄^2- and ADP, as a product. In the next step, a proton is shifted back from Glu1742 to generate inorganic phosphate (H₂PO₄^-) via another relay proton transfer event. The overall activation barrier for the Glu1742 assisted ATP hydrolysis is found to be the most favourable pathway compared to other plausible pathways. We also unearthed that ATP hydrolysis in dynein follows a so-called associative-like pathway in its rate-limiting step. Our study ascertained the important indirect roles of the two amino acids (such as Arg2109, Asn1792) and Mg²⁺ ion in the ATP hydrolysis of dynein. Additionally, multiple sequence alignment of the different organisms of dynein motors has conveyed the evolutionary importance of the Glu1742, Asn1742, and Arg2109 residues, respectively. As similar mechanisms are also prevalent in other motors, and GTPase and ATPase enzymes, the present finding spells out the definitive requirement of a proton relay process through an extended water-chain as one of the key components in an enzymatic ATP hydrolysis reaction.

Using Graphs of Dynamic Hydrogen-Bond Networks To Dissect Conformational Coupling in a Protein Motor

DExD/H-box proteins are soluble enzymes that couple binding and hydrolysis of adenosine triphosphate (ATP) with reactions involving RNA metabolism or bind and push newly synthesized proteins across bacterial cell membranes. Knowledge of the reaction mechanism of these enzymes could help the development of new therapeutics. In order to explore the mechanism of long-distance conformational coupling in SecA, the DEAD-box motor of the Sec protein secretion in bacteria, we implemented algorithms that provide simplified graph representations of the protein's dynamic hydrogen-bond networks. We find that mutations near the nucleotide-binding site or changes of the nucleotide-binding state of SecA associate with altered dynamics at the preprotein binding domain and identify extended networks of hydrogen bonds that connect the active site of SecA to the region where SecA binds newly synthesized secretory proteins. Water molecules participate in hydrogen-bonded water chains that bridge functional domains of SecA and could contribute to long-distance conformational coupling.

Dynamic water bridging and proton transfer at a surface carboxylate cluster of photosystem II

A hydrogen-bond cluster at a negatively-charged protein interface with a bound protein and long-lived waters might be a proton storage site.

Mapping Free Energy Pathways for ATP Hydrolysis in the E. coli ABC Transporter HlyB by the String Method

HlyB functions as an adenosine triphosphate (ATP)-binding cassette (ABC) transporter that enables bacteria to secrete toxins at the expense of ATP hydrolysis. Our previous work, based on potential energy profiles from combined quantum mechanical and molecular mechanical (QM/MM) calculations, has suggested that the highly conserved H-loop His residue H662 in the nucleotide binding domain (NBD) of E. coli HlyB may catalyze the hydrolysis of ATP through proton relay. To further test this hypothesis when entropic contributions are taken into account, we obtained QM/MM minimum free energy paths (MFEPs) for the HlyB reaction, making use of the string method in collective variables. The free energy profiles along the MFEPs confirm the direct participation of H662 in catalysis. The MFEP simulations of HlyB also reveal an intimate coupling between the chemical steps and a local protein conformational change involving the signature-loop residue S607, which may serve a catalytic role similar to an Arg-finger motif in many ATPases and GTPases in stabilizing the phosphoryl-transfer transition state.

Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations

We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.

Mechanism by which water and protein electrostatic interactions control proton transfer at the active site of channelrhodopsin

Channelrhodopsins are light-sensitive ion channels whose reaction cycles involve conformation-coupled transfer of protons. Understanding how channelrhodopsins work is important for applications in optogenetics, where light activation of these proteins triggers changes in the transmembrane potential across excitable membranes. A fundamental open question is how the protein environment ensures that unproductive proton transfer from the retinal Schiff base to the nearby carboxylate counterion is avoided in the resting state of the channel. To address this question, we performed combined quantum mechanical/molecular mechanical proton transfer calculations with explicit treatment of the surrounding lipid membrane. The free energy profiles computed for proton transfer to the counterion, either via a direct jump or mediated by a water molecule, demonstrate that, when retinal is all-trans, water and protein electrostatic interactions largely favour the protonated retinal Schiff base state. We identified a conserved lysine group as an essential structural element for the proton transfer energetics in channelrhodopsins.

Thiacalix[4]arenes Remove the Inhibitory Effects of Zn Cations on the Myosin ATPase Activity

Numerous female reproductive abnormalities are caused by uterine smooth muscle (myometrium) disorders. Heavy metals have an adverse effect on the contractility of the uterine smooth muscle. Although zinc is an essential biogenic element for most of the organisms, high doses of this element are toxic. The study of 0.5-5 mM Zn2+ effect on myosin S1 ATPase activity from the uterus found that 5 mM Zn2+ cations have the most pronounced inhibitory effect. The calculation of the kinetic parameters (Km and Vmax, ATP) revealed that the apparent maximum velocity of the hydrolysis ATP catalyzed by myosin in the presence of 5 mM Zn2+ decreased by 1.6 times. The value of Кm for ATP hydrolysis by myosin S1 in the presence of Zn2+ does not change statistically, although it tends to decrease. It was determined that uterine myosin S1 ATPase activity does not depend on the concentration of Mg2+ in the presence of 5 mM Zn2+. Also, it was demonstrated that tetrahydroxythiacalix[4]arene-tetrasulfosphonate (C-798) and tetrahydroxythiacalix[4]arene-tetraphosphonate (C-800) restored myosin S1 ATPase activity to the control level in the presence of 5 mM Zn2+. One of the most probable mechanisms of restoring the action of these thiacalix[4]arenes protective effect is based on its ability to chelate heavy metal cations from the incubation medium. The molecular docking of C-798 and C-800 into the myosin S1 region showed that these thiacalix[4]arenes could interact with Zn cation bond by myosin amino acid residues near the ATPase active site. Therefore, thiacalix[4]arenes may weaken the interaction between this cation and myosin S1. It was speculated that the obtained results could be used for further research with the aim of using this thiacalix[4]arenes as pharmacological compounds in the case of poisoning with high concentrations of zinc.

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.

A QM/MM approach on the structural and stereoelectronic factors governing glycosylation by GTF-SI fromStreptococcus mutans

This manuscript contains novel insights into the reaction mechanism catalyzed by GTF-SI. Structural and electronic features of the system are revealed, such as the strong hydrogen bond depicted above.

Theoretical Study of the Phosphoryl Transfer Reaction from ATP to Dha Catalyzed by DhaK from Escherichia coli

Protein kinases, representing one of the largest protein families involved in almost all aspects of cell life, have become one of the most important targets for the development of new drugs to be used in, for instance, cancer treatments. In this article an exhaustive theoretical study of the phosphoryl transfer reaction from adenosine triphosphate (ATP) to dihydroxyacetone (Dha) catalyzed by DhaK from Escherichia coli (E. coli) is reported. Two different mechanisms, previously proposed for the phosphoryl transfer from ATP to the hydroxyl side chain of specific serine, threonine, or tyrosine residues, have been explored based on the generation of free energy surfaces (FES) computed with hybrid QM/MM potentials. The results suggest that the substrate-assisted phosphoryl and proton-transfer mechanism is kinetically more favorable than the mechanism where an aspartate would be activating the Dha. Although the details of the mechanisms appear to be dramatically dependent on the level of theory employed in the calculations (PM3/MM, B3LYP:PM3/MM, or B3LYP/MM), the transition states (TSs) for the phosphoryl transfer step appear to be described as a concerted step with different degrees of synchronicity in the breaking and forming bonds process in both explored mechanisms. Residues of the active site belonging to different subunits of the protein, such as Gly78B, Thr79A, Ser80A, Arg178B, and one Mg²⁺ cation, would be stabilizing the transferred phosphate in the TS. Asp109A would have a structural role by posing the Dha and other residues of the active site in the proper orientation. The information derived from our calculations not only reveals the role of the enzyme and the particular residues of its active site, but it can assist in the rational design of new more specific inhibitors.

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.

Metal Fluorides: Tools for Structural and Computational Analysis of Phosphoryl Transfer Enzymes

The phosphoryl group, PO₃^-, is the dynamic structural unit in the biological chemistry of phosphorus. Its transfer from a donor to an acceptor atom, with oxygen much more prevalent than nitrogen, carbon, or sulfur, is at the core of a great majority of enzyme-catalyzed reactions involving phosphate esters, anhydrides, amidates, and phosphorothioates. The serendipitous discovery that the phosphoryl group could be labeled by "nuclear mutation," by substitution of PO₃^- by MgF₃^- or AlF₄^-, has underpinned the application of metal fluoride (MF _x ) complexes to mimic transition states for enzymatic phosphoryl transfer reactions, with sufficient stability for experimental analysis. Protein crystallography in the solid state and ¹⁹F NMR in solution have enabled direct observation of ternary and quaternary protein complexes embracing MF _x transition state models with precision. These studies have underpinned a radically new mechanistic approach to enzyme catalysis for a huge range of phosphoryl transfer processes, as varied as kinases, phosphatases, phosphomutases, and phosphohydrolases. The results, without exception, have endorsed trigonal bipyramidal geometry (tbp) for concerted, "in-line" stereochemistry of phosphoryl transfer. QM computations have established the validity of tbp MF _x complexes as reliable models for true transition states, delivering similar bond lengths, coordination to essential metal ions, and virtually identical hydrogen bond networks. The emergence of protein control of reactant orbital overlap between bond-forming species within enzyme transition states is a new challenging theme for wider exploration.

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.

The GTPase hGBP1 converts GTP to GMP in two steps via proton shuttle mechanisms

GTPases play a crucial role in the regulation of many biological processes by catalyzing the hydrolysis of GTP into GDP. The focus of this work is on the dynamin-related large GTPase human guanine nucleotide binding protein-1 (hGBP1) which is able to hydrolyze GTP even to GMP. Here, we studied the largely unknown mechanisms of both GTP and GDP hydrolysis steps utilizing accelerated ab initio QM/MM metadynamics simulations to compute multi-dimensional free energy landscapes. We find an indirect substrate-assisted catalysis (SAC) mechanism for GTP hydrolysis involving transfer of a proton from the water nucleophile to a nonbridging phosphoryl oxygen via a proton relay pathway where the rate-determining first step is concerted-dissociative nature. A "composite base" consisting of Ser73, Glu99, a bridging water molecule, and GTP was found to activate the nucleophilic water, thus disclosing the complex nature of the general base in hGBP1. A nearly two-fold reduction in the free energy barrier was obtained for GTP hydrolysis in the enzyme in comparison to bulk solvent. The subsequent GDP hydrolysis in hGBP1 was also found to follow a water-mediated proton shuttle mechanism. It is expected that the proton shuttle mechanisms unravelled for hGBP1 apply to many classes of GTPases/ATPases that possess an optimally-arranged hydrogen bonding network, which connects the catalytic water to a proton acceptor.

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.

A Networks Approach to Modeling Enzymatic Reactions

Modeling enzymatic reactions is a demanding task due to the complexity of the system, the many degrees of freedom involved and the complex, chemical, and conformational transitions associated with the reaction. Consequently, enzymatic reactions are not determined by precisely one reaction pathway. Hence, it is beneficial to obtain a comprehensive picture of possible reaction paths and competing mechanisms. By combining individually generated intermediate states and chemical transition steps a network of such pathways can be constructed. Transition networks are a discretized representation of a potential energy landscape consisting of a multitude of reaction pathways connecting the end states of the reaction. The graph structure of the network allows an easy identification of the energetically most favorable pathways as well as a number of alternative routes.

Theoretical studies of the function switch and mechanism of AceK as a highly active ATPase

As a multi-function enzyme, AceK integrates kinase, phosphatase and ATPase activities in a single active site and these functions are delicately regulated..

Mechanism of conformational coupling in SecA: Key role of hydrogen-bonding networks and water interactions

SecA uses the energy yielded by the binding and hydrolysis of adenosine triphosphate (ATP) to push secretory pre-proteins across the plasma membrane in bacteria. Hydrolysis of ATP occurs at the nucleotide-binding site, which contains the conserved carboxylate groups of the DEAD-box helicases. Although crystal structures provide valuable snapshots of SecA along its reaction cycle, the mechanism that ensures conformational coupling between the nucleotide-binding site and the other domains of SecA remains unclear. The observation that SecA contains numerous hydrogen-bonding groups raises important questions about the role of hydrogen-bonding networks and hydrogen-bond dynamics in long-distance conformational couplings. To address these questions, we explored the molecular dynamics of SecA from three different organisms, with and without bound nucleotide, in water. By computing two-dimensional hydrogen-bonding maps we identify networks of hydrogen bonds that connect the nucleotide-binding site to remote regions of the protein, and sites in the protein that respond to specific perturbations. We find that the nucleotide-binding site of ADP-bound SecA has a preferred geometry whereby the first two carboxylates of the DEAD motif bridge via hydrogen-bonding water. Simulations of a mutant with perturbed ATP hydrolysis highlight the water-bridged geometry as a key structural element of the reaction path.

Advances in quantum simulations of ATPase catalysis in the myosin motor

During its contraction cycle, the myosin motor catalyzes the hydrolysis of ATP. Several combined quantum/classical mechanics (QM/MM) studies of this step have been published, which substantially contributed to our thinking about the catalytic mechanism. The methodological difficulties encountered over the years in the simulation of this complex reaction are now understood: (a) Polarization of the protein peptide groups surrounding the highly charged ATP(4-) cannot be neglected. (b) Some unsuspected protein groups need to be treated QM. (c) Interactions with the γ-phosphate versus the β-phosphate favor a concurrent versus a sequential mechanism, respectively. Thus, these practical aspects strongly influence the computed mechanism, and should be considered when studying other catalyzed phosphor-ester hydrolysis reactions, such as in ATPases or GTPases.