Quantifying the mechanism of phosphate monoester hydrolysis in aqueous solution by evaluating the relevant ab initio QM/MM free-energy surfaces

Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1

Collagen hydrolysis, catalyzed by Zn(II)-dependent matrix metalloproteinases (MMPs), is a critical physiological process. Despite previous computational investigations into the catalytic mechanisms of MMP-mediated collagenolysis, a significant knowledge gap in understanding remains regarding the influence of conformational sampling and entropic contributions at physiological temperature on enzymatic collagenolysis. In our comprehensive multilevel computational study, employing quantum mechanics/molecular mechanics (QM/MM) metadynamics (MetD) simulations, we aimed to bridge this gap and provide valuable insights into the catalytic mechanism of MMP-1. Specifically, we compared the full enzyme-substrate complex in solution, clusters in solution, and gas-phase to elucidate insights into MMP-1-catalyzed collagenolysis. Our findings reveal significant differences in the catalytic mechanism when considering thermal effects and the dynamic evolution of the system, contrasting with conventional static potential energy surface QM/MM reaction path studies. Notably, we observed a significant stabilization of the critical tetrahedral intermediate, attributed to contributions from conformational flexibility and entropy. Moreover, we found that protonation of the scissile bond nitrogen occurs via proton transfer from a Zn(II)-coordinated hydroxide rather than from a solvent water molecule. Following C-N bond cleavage, the C-terminus remains coordinated to the catalytic Zn(II), while the N-terminus forms a hydrogen bond with a solvent water molecule. Subsequently, the release of the C-terminus is facilitated by the coordination of a water molecule. Our study underscores the pivotal role of protein conformational dynamics at physiological temperature in stabilizing the transition state of the rate-limiting step and key intermediates, compared to the corresponding reaction in solution. These fundamental insights into the mechanism of collagen degradation provide valuable guidance for the development of MMP-1-specific inhibitors.

Mechanism of Phosphate Release from Actin Filaments

After ATP-actin monomers assemble filaments, the ATP's γ-phosphate is hydrolyzed within seconds and dissociates over minutes. We used all-atom molecular dynamics simulations to sample the release of phosphate from filaments and study residues that gate release. Dissociation of phosphate from Mg2+ is rate limiting and associated with an energy barrier of 20 kcal/mol, consistent with experimental rates of phosphate release. Phosphate then diffuses in an internal cavity toward a gate formed by R177 suggested in prior computational studies and cryo-EM structures. The gate is closed when R177 hydrogen bonds with N111 and is open when R177 forms a salt bridge with D179. Most of the time interactions of R177 with other residues occludes the phosphate release pathway. Machine learning analysis reveals that the occluding interactions fluctuate rapidly, underscoring the secondary role of backdoor gate opening in Pi release, in contrast with the previous hypothesis that gate opening is the primary event.

Density functional theory explorations of parathion and paraoxon hydrolysis as a function of the underlying alkaline environment

Parathion, a once commonly used pesticide known for its potential toxicity, can follow several degradation mechanisms in the environment. Given the species stability and persistence, parathion can be washed into waterways from rain, and therefore an atomistic perspective of the hydrolysis of parathion, and its byproduct paraoxon, is required in order to understand its fate in the environment. Experimental studies have determined that pH plays an important role in the calculated hydrolysis rate constants of parathion degradation. In this work, the degradation of parathion into either paraoxon or 4-nitrophenol, and the degradation of paraoxon to 4-nitrophenol are explored through density functional theory using the M06-2X functional. How the level of basicity affects the reaction mechanism is explored through two different hydroxide/water environments. Our calculations support the anticipated mechanisms determined by previous experimental work that the formation of 4-nitrophenol is the predominant pathway in hydrolysis of parathion.

Multiscale computation delivers organophosphorus reactivity and stereoselectivity to immunoglobulin scavengers

Quantum mechanics/molecular mechanics (QM/MM) maturation of an immunoglobulin (Ig) powered by supercomputation delivers novel functionality to this catalytic template and facilitates artificial evolution of biocatalysts. We here employ density functional theory-based (DFT-b) tight binding and funnel metadynamics to advance our earlier QM/MM maturation of A17 Ig-paraoxonase (WTIgP) as a reactibody for organophosphorus toxins. It enables regulation of biocatalytic activity for tyrosine nucleophilic attack on phosphorus. The single amino acid substitution l-Leu47Lys results in 340-fold enhanced reactivity for paraoxon. The computed ground-state complex shows substrate-induced ionization of the nucleophilic l-Tyr37, now H-bonded to l-Lys47, resulting from repositioning of l-Lys47. Multiple antibody structural homologs, selected by phenylphosphonate covalent capture, show contrasting enantioselectivities for a P-chiral phenylphosphonate toxin. That is defined by crystallographic analysis of phenylphosphonylated reaction products for antibodies A5 and WTIgP. DFT-b analysis using QM regions based on these structures identifies transition states for the favored and disfavored reactions with surprising results. This stereoselection analysis is extended by funnel metadynamics to a range of WTIgP variants whose predicted stereoselectivity is endorsed by experimental analysis. The algorithms used here offer prospects for tailored design of highly evolved, genetically encoded organophosphorus scavengers and for broader functionalities of members of the Ig superfamily, including cell surface-exposed receptors.

Insights into Catalytic Gas-Phase Hydrolysis of Organophosphate Chemical Warfare Agents by MOF-Supported Bimetallic Metal-Oxo Clusters

Zirconium-based metal-organic frameworks (Zr-MOFs) have been reported to be efficient catalysts for the hydrolysis of organophosphate chemical warfare agents (CWAs) in buffered solutions. However, for the gas-phase reaction, which is more relevant to the situation in a battlefield gas mask application, the kinetics of Zr-MOF catalysts may be severely hindered by strong product inhibition. To improve the catalytic performance, we computationally screened a series of synthetically accessible Zr-MOF-supported bimetallic metal-oxo clusters in which the metal-oxygen-metal active motif is preserved, aiming to find catalysts that have lower binding affinities to the hydrolysis product. For the promising catalyst Al₂O₂(OH)₂@NU-1000 identified from the screening using density functional theory, we mapped out the full reaction pathway of gas-phase dimethyl p-nitrophenolphosphate (DMNP) hydrolysis and analyzed the free energy profile as well as the turnover frequency (TOF). We found that the catalytic mechanism on the new catalyst is slightly different from the one on NU-1000, which also led to a different TOF-limiting step. Additional factors that can affect the overall catalytic performance in practical application, such as the amount of ambient moisture and the existence of acid gases that may poison the catalyst, have also been evaluated.

Diversity of mechanisms in Ras-GAP catalysis of guanosine triphosphate hydrolysis revealed by molecular modeling

The mechanism of the deceptively simple reaction of guanosine triphosphate (GTP) hydrolysis catalyzed by the cellular protein Ras in complex with the activating protein GAP is an important issue because of the significance of this reaction in cancer research. We show that molecular modeling of GTP hydrolysis in the Ras-GAP active site reveals a diversity of mechanisms of the intrinsic chemical reaction depending on molecular groups at position 61 in Ras occupied by glutamine in the wild-type enzyme. First, a comparison of reaction energy profiles computed at the quantum mechanics/molecular mechanics (QM/MM) level shows that an assignment of the Gln61 side chain in the wild-type Ras either to QM or to MM parts leads to different scenarios corresponding to the glutamine-assisted or the substrate-assisted mechanisms. Second, replacement of Gln61 by the nitro-analog of glutamine (NGln) or by Glu, applied in experimental studies, results in two more scenarios featuring the so-called two-water and the concerted-type mechanisms. The glutamine-assisted mechanism in the wild-type Ras-GAP, in which the conserved Gln61 plays a decisive role, switching between the amide and imide tautomer forms, is consistent with the known experimental results of structural, kinetic and spectroscopy studies. The results emphasize the role of the Ras residue Gln61 in Ras-GAP catalysis and explain the retained catalytic activity of the Ras-GAP complex towards GTP hydrolysis in the Gln61NGln and Gln61Glu mutants of Ras.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Mapping the working route of phosphate monoester hydrolysis catalyzed by copper based models with special emphasis on the role of oxoanions by experimental and theoretical studies

The mechanistic pathway of phosphate-ester bond hydrolysis with special emphasis on the role of oxoanions was explored by experimental and theoretical study.

EF-Tu and EF-G are activated by allosteric effects

Many cellular processes are controlled by GTPases, and gaining quantitative understanding of the activation of such processes has been a major challenge. In particular, it is crucial to obtain reliable free-energy surfaces for the relevant reaction paths both in solution and in GTPases active sites. Here, we revisit the energetics of the activation of EF-G and EF-Tu by the ribosome and explore the nature of the catalysis of the GTPase reaction. The comparison of EF-Tu to EF-G allows us to explore the impact of possible problems with the available structure of EF-Tu. Additionally, mutational effects are used for a careful validation of the emerging conclusions. It is found that the reaction may proceed by both a two-water mechanism and a one-water (GTP as a base) mechanism. However, in both cases, the activation involves a structural allosteric effect, which is likely to be a general-activation mechanism for all GTPases.

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.

Mechanistic Insights on Human Phosphoglucomutase Revealed by Transition Path Sampling and Molecular Dynamics Calculations

Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted S_N 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg²⁺ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM.

Phosphoryl Transfer Reaction in RNA: Is the Substrate-Assisted Catalysis a Possible Mechanism in Certain Solvents?

A proton shuttle mechanism for the phosphoryl transfer reaction in RNA, in which a proton is transferred from the nucleophile to the leaving group through a nonbridged oxygen atom of the phosphate, was explored using the MO6-2X density functional method and the solvent continuum model. This reaction is the initial step of the RNA hydrolysis. We used different solvents characterized by their dielectric constant, and, for each of them, we studied the nuclear and electronic relaxation, produced by the solvent reaction field, for the stationary points. Given that RNA has a poor leaving group, the bond breaking corresponds to the rate-determining step. If the O atom is substituted by a S atom, the leaving group is now good, and the rate-determining step is now the nucleophilic attack concerted with the proton transfer. The most relevant result we found is that none of the solvents we studied has a free energy of activation that is smaller than the one in water. This suggests that the enzyme catalysis following this mechanism must be due to the permanent electric field that is created by a preorganized charge distribution but not to the solvent reaction field.

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.

The effect of magnesium ions on triphosphate hydrolysis

The role of metal ions in catalyzing phosphate ester hydrolysis has been the subject of much debate, both in terms of whether they change the transition state structure or mechanistic pathway. Understanding the impact of metal ions on these biologically critical reactions is central to improving our understanding of the role of metal ions in the numerous enzymes that facilitate them. In the present study, we have performed density functional theory studies of the mechanisms of methyl triphosphate and acetyl phosphate hydrolysis in aqueous solution to explore the competition between solvent- and substrate-assisted pathways, and examined the impact of Mg2+ on the energetics and transition state geometries. In both cases, we observe a clear preference for a more dissociative solvent-assisted transition state, which is not significantly changed by coordination of Mg2+. The effect of Mg2+ on the transition state geometries for the two pathways is minimal. While our calculations cannot rule out a substrate-assisted pathway as a possible solution for biological phosphate hydrolysis, they demonstrate that a significantly higher energy barrier needs to be overcome in the enzymatic reaction for this to be an energetically viable reaction pathway.

Exploring the Dependence of QM/MM Calculations of Enzyme Catalysis on the Size of the QM Region

Although QM/MM calculations are the primary current tool for modeling enzymatic reactions, the reliability of such calculations can be limited by the size of the QM region. Thus, we examine in this work the dependence of QM/MM calculations on the size of the QM region, using the reaction of catechol-O-methyl transferase (COMT) as a test case. Our study focuses on the effect of adding residues to the QM region on the activation free energy, obtained with extensive QM/MM sampling. It is found that the sensitivity of the activation barrier to the size of the QM is rather limited, while the dependence of the reaction free energy is somewhat larger. Of course, the results depend on the inclusion of the first solvation shell in the QM regions. For example, the inclusion of the Mg²⁺ ion can change the activation barrier due to charge transfer effects. However, such effects can easily be included in semiempirical approaches by proper parametrization. Overall, we establish that QM/MM calculations of activation barriers of enzymatic reactions are not highly sensitive to the size of the QM region, beyond the immediate region that describes the reacting atoms.

Exploring the mechanism of DNA polymerases by analyzing the effect of mutations of active site acidic groups in Polymerase β

Elucidating the catalytic mechanism of DNA polymerase is crucial for a progress in the understanding of the control of replication fidelity. This work tries to advance the mechanistic understanding by analyzing the observed effect of mutations of the acidic groups in the active site of Polymerase β as well as the pH effect on the rate constant. The analysis involves both empirical valence bond (EVB) free energy calculations and considerations of the observed pH dependence of the reaction. The combined analysis indicates that the proton transfer (PT) from the nucleophilic O3' has two possible pathways, one to D256 and the second to the bulk. We concluded based on calculations and the experimental pH profile that the most likely path for the wild-type (WT) and the D256E and D256A mutants is a PT to the bulk, although the WT may also use a PT to Asp 256. Our analysis highlights the need for very extensive sampling in the calculations of the activation barrier and also clearly shows that ab initio QM/MM calculations that do not involve extensive sampling are unlikely to give a clear quantitative picture of the reaction mechanism. Proteins 2016; 84:1644-1657. © 2016 Wiley Periodicals, Inc.

The Competing Mechanisms of Phosphate Monoester Dianion Hydrolysis

Despite the numerous experimental and theoretical studies on phosphate monoester hydrolysis, significant questions remain concerning the mechanistic details of these biologically critical reactions. In the present work we construct a linear free energy relationship for phosphate monoester hydrolysis to explore the effect of modulating leaving group pK_a on the competition between solvent- and substrate-assisted pathways for the hydrolysis of these compounds. Through detailed comparative electronic-structure studies of methyl phosphate and a series of substituted aryl phosphate monoesters, we demonstrate that the preferred mechanism is dependent on the nature of the leaving group. For good leaving groups, a strong preference is observed for a more dissociative solvent-assisted pathway. However, the energy difference between the two pathways gradually reduces as the leaving group pK_a increases and creates mechanistic ambiguity for reactions involving relatively poor alkoxy leaving groups. Our calculations show that the transition-state structures vary smoothly across the range of pK_as studied and that the pathways remain discrete mechanistic alternatives. Therefore, while not impossible, a biological catalyst would have to surmount a significantly higher activation barrier to facilitate a substrate-assisted pathway than for the solvent-assisted pathway when phosphate is bonded to good leaving groups. For poor leaving groups, this intrinsic preference disappears.

Enhancing Paradynamics for QM/MM Sampling of Enzymatic Reactions

Despite the enormous increase in computer power, it is still extremely challenging to obtain computationally converging sampling of ab initio QM/MM (QM(ai)/MM) free energy surfaces in condensed phases. The sampling problem can be significantly reduced by the use of the reference potential paradynamics (PD) approach, but even this approach still requires major computer time in studies of enzymatic reactions. To further reduce the sampling problem we developed here a new PD version where we use an empirical valence bond reference potential that has a minimum rather than a maximum at the transition state region of the target potential (this is accomplished conveniently by shifting the EVB of the product state). Hence, we can map the TS region in a more efficient way. Here, we introduce and validate the inverted EVB PD approach. The validation involves the study of the S(N)2 step of the reaction catalyzed by haloakene dehalogenase (DhlA) and the GTP hydrolysis in the RasGAP system. In addition, we have also studied the corresponding reaction in water for each of the systems described here and the reaction involving trimethylsulfonium and dimethylamine in solution. The results are encouraging and the new strategy appears to provide a powerful way of evaluating QM(ai)/MM activation free energies.

Phosphorane lifetime and stereo-electronic effects along the alkaline hydrolysis of phosphate esters

Ab initio molecular dynamics simulations revealed that phosphorane, an important intermediate in the hydrolysis of phosphate diesters, has a lifetime of ∼1 ps in aqueous solution. QTAIM and EDA analyses along the reaction coordinate show that the hydrolysis reaction of phosphate esters is driven mainly by electrostatic interactions.

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

A computational study of the phosphoryl transfer reaction between ATP and Dha in aqueous solution

Phosphoryl transfer reactions are ubiquitous in biology, being involved in processes ranging from energy and signal transduction to the replication genetic material. Dihydroxyacetone phosphate (Dha-P), an intermediate of the synthesis of pyruvate and a very important building block in nature, can be generated by converting free dihydroxyacetone (Dha) through the action of the dihydroxyacetone kinase enzyme. In this paper the reference uncatalyzed reaction in solution has been studied in order to define the foundations of the chemical reaction and to determine the most adequate computational method to describe this electronically complex reaction. In particular, the phosphorylation reaction mechanism between adenosine triphosphate (ATP) and Dha in aqueous solution has been studied by means of quantum mechanics/molecular mechanics (QM/MM) Molecular Dynamics (MD) simulations with the QM subset of atoms described with semi-empirical and DFT methods. The results appear to be strongly dependent on the level of calculation, which will have to be taken into account for future studies of the reaction catalyzed by enzymes. In particular, PM3/MM renders lower free energy barriers and a less endergonic process than AM1d/MM and PM6/MM methods. Nevertheless, the concerted pathway was not located with the former combination of potentials.

Linear free energy relationships in RNA transesterification: theoretical models to aid experimental interpretations

RNA cleavage transesterification is of fundamental reaction in biology that is catalyzed by both protein and RNA enzymes. In this work, a series of RNA transesterification model reactions with a wide range of leaving groups are investigated with density-functional calculations in an aqueous solvation environment in order to study linear free energy relationships (LFERs) and their connection to transition state structure and bonding. Overall, results obtained from the polarizable continuum solvation model with UAKS radii produce the best linear correlations and closest overall agreement with experimental results. Reactions with a poor leaving group are predicted to proceed via a stepwise mechanism with a late transition state that is rate controlling. As leaving group becomes more acidic and labile, the barriers of both early and late transition states decrease. LFERs for each transition state are computed, with the late transition state barrier showing greater sensitivity to leaving group pKa. For sufficiently enhanced leaving groups, the reaction mechanism transits to a concerted mechanism characterized by a single early transition state. Further linear relationships were derived for bond lengths and bond orders as a function of leaving group pKa and rate constant values that can be used for prediction. This work provides important benchmark linear free energy data that allows a molecular-level characterization of the structure and bonding of the transition states for this important class of phosphoryl transfer reactions. The relations reported herein can be used to aid in the interpretation of data obtained from experimental studies of non-catalytic and catalytic mechanisms.

Computational characterization of the chemical step in the GTP hydrolysis by Ras-GAP for the wild-type and G13V mutated Ras

The free energy profiles for the chemical reaction of the guanosine triphosphate hydrolysis GTP + H2O → GDP + Pi by Ras-GAP for the wild-type and G13V mutated Ras were computed by using molecular dynamics protocols with the QM(ab initio)/MM potentials. The results are consistent with the recent measurements of reaction kinetics in Ras-GAP showing about two-order reduction of the rate constant upon G13V mutation in Ras: the computed activation barrier on the free energy profile is increased by 3 kcal/mol upon the G13V replacement. The major reason for a higher energy barrier is a shift of the "arginine finger" (R789 from GAP) from the favorable position in the active site. The results of simulations provide support for the mechanism of the reference reaction according to which the Q61 side chain directly participates in chemical transformations at the proton transfer stage.

QM/MM investigation of ATP hydrolysis in aqueous solution

Adenosine-5'-triphosphate (ATP) hydrolysis represents a most important reaction in biology. Despite extensive research efforts, the mechanism for ATP hydrolysis in aqueous solution still remains under debate. Previous theoretical studies often predefined reaction coordinates to characterize the mechanism for ATP hydrolysis in water with Mg(2+) by evaluating free energy profiles through these preassumed reaction paths. In the present work, a nudged elastic band method is applied to identify the minimum energy path calculated with a hybrid quantum mechanics and molecular mechanics approach. Along the reaction path, the free energy profile was obtained to have a single transition state and the activation energy of 32.5 kcal/mol. This transition state bears a four-centered structure that describes a concerted nature of the reaction. In the More-O'Ferrall-Jencks diagram, the results show that the reaction proceeds through a concerted path before the system reaches the transition state and along an associative path after the transition state. In addition, the calculated reaction free energy is -7.0 kcal/mol, in good agreement with experiment, capturing the exothermic feature of MgATP(2-) hydrolysis in aqueous solution, whereas the reaction was often shown to be endothermic in the previous theoretical studies. As Mg(2+) is required for ATP hydrolysis in cells, its role in the reaction is also elucidated.

A QM/MM study of Kemptide phosphorylation catalyzed by protein kinase A. The role of Asp166 as a general acid/base catalyst

In this work a theoretical study of the γ-phosphoryl group transfer from ATP to Ser17 of the synthetic substrate Kemptide (LRRASLG) in protein kinase A (PKA) has been carried out with a solvated model of the PKA-Mg2ATP-Kemptide system based on the X-ray crystallographic structure. We have used high levels (B3LYP/MM and MP2/MM) of theory to determine the overall reaction paths of the so-called concerted loose mechanism trying to clarify some aspects of that mechanism still under debate. Our calculations demonstrate for the first time in a complete model of the ternary system the viability of the final step of the catalytic mechanism in which the protonation of the phosphokemptide product by Asp166 takes place. Asp166 is a base catalyst that abstracts the HγSer17 of Kemptide thus facilitating the phosphoryl transfer, but it also acts as an acid catalyst by donating the proton just accepted from Ser17 to the O2γATP atom of the phosphoryl group.

Finding Chemical Reaction Paths with a Multilevel Preconditioning Protocol

Finding transition paths for chemical reactions can be computationally costly owing to the level of quantum-chemical theory needed for accuracy. Here, we show that a multilevel preconditioning scheme that was recently introduced (Tempkin et al. J. Chem. Phys.2014, 140, 184114) can be used to accelerate quantum-chemical string calculations. We demonstrate the method by finding minimum-energy paths for two well-characterized reactions: tautomerization of malonaldehyde and Claissen rearrangement of chorismate to prephanate. For these reactions, we show that preconditioning density functional theory (DFT) with a semiempirical method reduces the computational cost for reaching a converged path that is an optimum under DFT by several fold. The approach also shows promise for free energy calculations when thermal noise can be controlled.

Resolving apparent conflicts between theoretical and experimental models of phosphate monoester hydrolysis

Understanding phosphoryl and sulfuryl transfer is central to many biochemical processes. However, despite decades of experimental and computational studies, a consensus concerning the precise mechanistic details of these reactions has yet to be reached. In this work we perform a detailed comparative theoretical study of the hydrolysis of p-nitrophenyl phosphate, methyl phosphate and p-nitrophenyl sulfate, all of which have served as key model systems for understanding phosphoryl and sulfuryl transfer reactions, respectively. We demonstrate the existence of energetically similar but mechanistically distinct possibilities for phosphate monoester hydrolysis. The calculated kinetic isotope effects for p-nitrophenyl phosphate provide a means to discriminate between substrate- and solvent-assisted pathways of phosphate monoester hydrolysis, and show that the solvent-assisted pathway dominates in solution. This preferred mechanism for p-nitrophenyl phosphate hydrolysis is difficult to find computationally due to the limitations of compressing multiple bonding changes onto a 2-dimensional energy surface. This problem is compounded by the need to include implicit solvation to at least microsolvate the system and stabilize the highly charged species. In contrast, methyl phosphate hydrolysis shows a preference for a substrate-assisted mechanism. For p-nitrophenyl sulfate hydrolysis there is only one viable reaction pathway, which is similar to the solvent-assisted pathway for phosphate hydrolysis, and the substrate-assisted pathway is not accessible. Overall, our results provide a unifying mechanistic framework that is consistent with the experimentally measured kinetic isotope effects and reconciles the discrepancies between theoretical and experimental models for these biochemically ubiquitous classes of reaction.

Modeling the role of G12V and G13V Ras mutations in the Ras-GAP-catalyzed hydrolysis reaction of guanosine triphosphate

Cancer-associated point mutations in Ras, in particular, at glycine 12 and glycine 13, affect the normal cycle between inactive GDP-bound and active GTP-bound states. In this work, the role of G12V and G13V replacements in the GAP-stimulated intrinsic GTP hydrolysis reaction in Ras is studied using molecular dynamics (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials. A model molecular system was constructed by motifs of the relevant crystal structure (Protein Data Bank entry 1WQ1 ). QM/MM optimization of geometry parameters in the Ras-GAP-GTP complex and QM/MM-MD simulations were performed with a quantum subsystem comprising a large fraction of the enzyme active site. For the system with wild-type Ras, the conformations fluctuated near the structure ready to be involved in the efficient chemical reaction leading to the cleavage of the phosphorus-oxygen bond in GTP upon approach of the properly aligned catalytic water molecule. Dynamics of the system with the G13V mutant is characterized by an enhanced flexibility in the area occupied by the γ-phosphate group of GTP, catalytic water, and the side chains of Arg789 and Gln61, which should somewhat hinder fast chemical steps. Conformational dynamics of the system with the G12V mutant shows considerable displacement of the Gln61 side chain and catalytic water from their favorable arrangement in the active site that may lead to a marked reduction in the reaction rate. The obtained computational results correlate well with the recent kinetic measurements of the Ras-GAP-catalyzed hydrolysis of GTP.

Simulating the catalytic effect of a designed mononuclear zinc metalloenzyme that catalyzes the hydrolysis of phosphate triesters

One of the greatest challenges in biotechnology and in biochemistry is the ability to design efficient enzymes. In fact, such an ability would be one of the most convincing manifestations of a full understanding of the origin of enzyme catalysis. Despite some progress on this front, most of the advances have been made by placing the reacting fragments in the proper places rather than by optimizing the preorganization of the environment, which is the key factor in enzyme catalysis. A rational improvement of the preorganization and a consistent assessment of the effectiveness of different design options require approaches capable of evaluating reliably the actual catalytic effect. In this work we examine the ability of the empirical valence bond (EVB) to reproduce the results of directed evolution improvements of the catalysis of diethyl 7-hydroxycoumarinyl by a designed mononuclear zinc metalloenzyme. Encouragingly, our study reproduced the catalytic effect obtained by directed evolution and offers a good start for further studies of this system.

Ribosome-induced tuning of GTP hydrolysis by a translational GTPase

GTP hydrolysis by elongation factor Tu (EF-Tu), a translational GTPase that delivers aminoacyl-tRNAs to the ribosome, plays a crucial role in decoding and translational fidelity. The basic reaction mechanism and the way the ribosome contributes to catalysis are a matter of debate. Here we use mutational analysis in combination with measurements of rate/pH profiles, kinetic solvent isotope effects, and ion dependence of GTP hydrolysis by EF-Tu off and on the ribosome to dissect the reaction mechanism. Our data suggest that--contrary to current models--the reaction in free EF-Tu follows a pathway that does not involve the critical residue H84 in the switch II region. Binding to the ribosome without a cognate codon in the A site has little effect on the GTPase mechanism. In contrast, upon cognate codon recognition, the ribosome induces a rearrangement of EF-Tu that renders GTP hydrolysis sensitive to mutations of Asp21 and His84 and insensitive to K(+) ions. We suggest that Asp21 and His84 provide a network of interactions that stabilize the positions of the γ-phosphate and the nucleophilic water, respectively, and thus play an indirect catalytic role in the GTPase mechanism on the ribosome.

Unraveling the mystery of ATP hydrolysis in actin filaments

Actin performs its myriad cellular functions by the growth and disassembly of its filamentous form. The hydrolysis of ATP in the actin filament has been shown to modulate properties of the filament, thus making it a pivotal regulator of the actin life cycle. Actin has evolved to selectively hydrolyze ATP in the filamentous form, F-actin, with an experimentally observed rate increase over the monomeric form, G-actin, of 4.3 × 10⁴. The cause of this dramatic increase in rate is investigated in this paper using extensive QM/MM simulations of both G- and F-actin. To compute the free energy of hydrolysis in both systems, metadynamics is employed along two collective variables chosen to describe the reaction coordinates of hydrolysis. F-actin is modeled as a monomer with restraints applied to coarse-grained variables enforced to keep it in a filament-like conformation. The simulations reveal a barrier height reduction for ATP hydrolysis in F-actin as compared to G-actin of 8 ± 1 kcal/mol, in good agreement with the experimentally measured barrier height reduction of 7 ± 1 kcal/mol. The barrier height reduction is influenced by an enhanced rotational diffusion of water in F-actin as compared to G-actin and shorter water wires between Asp154 and the nucleophilic water in F-actin, leading to more rapid proton transport.

QM/MM simulation (B3LYP) of the RNase A cleavage-transesterification reaction supports a triester A(N) + D(N) associative mechanism with an O2' H internal proton transfer

The mechanism of the backbone cleavage-transesterification step of the RNase A enzyme remains controversial even after 60 years of study. We report quantum mechanics/molecule mechanics (QM/MM) free energy calculations for two optimized reaction paths based on an analysis of all structural data and identified by a search for reaction coordinates using a reliable quantum chemistry method (B3LYP), equilibrated structural optimizations, and free energy estimations. Both paths are initiated by nucleophilic attack of the ribose O2' oxygen on the neighboring diester phosphate bond, and both reach the same product state (PS) (a O3'-O2' cyclic phosphate and a O5' hydroxyl terminated fragment). Path 1, resembles the widely accepted dianionic transition-state (TS) general acid (His119)/base (His12) classical mechanism. However, this path has a barrier (25 kcal/mol) higher than that of the rate-limiting hydrolysis step and a very loose TS. In Path 2, the proton initially coordinating the O2' migrates to the nonbridging O1P in the initial reaction path rather than directly to the general base resulting in a triester (substrate as base) AN + DN mechanism with a monoanionic weakly stable intermediate. The structures in the transition region are associative with low barriers (TS1 10, TS2 7.5 kcal/mol). The Path 2 mechanism is consistent with the many results from enzyme and buffer catalyzed and uncatalyzed analog reactions and leads to a PS consistent with the reactive state for the following hydrolysis step. The differences between the consistently estimated barriers in Path 1 and 2 lead to a 10(11) difference in rate strongly supporting the less accepted triester mechanism.

Quantitative exploration of the molecular origin of the activation of GTPase

GTPases play a major role in cellular processes, and gaining quantitative understanding of their activation demands reliable free energy surfaces of the relevant mechanistic paths in solution, as well as the interpolation of this information to GTPases. Recently, we generated ab initio quantum mechanical/molecular mechanical free energy surfaces for the hydrolysis of phosphate monoesters in solution, establishing quantitatively that the barrier for the reactions with a proton transfer (PT) step from a single attacking water (1 W) is higher than the one where the PT is assisted by a second water (2 W). The implication of this finding on the activation of GTPases is quantified here, by using the ab initio solution surfaces to calibrate empirical valence bond surfaces and then exploring the origin of the activation effect. It is found that, although the 2 W PT path is a new element, this step is not rate determining, and the catalytic effect is actually due to the electrostatic stabilization of the pre-PT transition state and the subsequent plateau. Thus, the electrostatic catalytic effect found in our previous studies of the Ras GTPase activating protein (RasGAP) and the elongation factor-Tu (EF-Tu) with a 1 W mechanism is still valid for the 2 W path. Furthermore, as found before, the corresponding activation appears to involve a major allosteric effect. Overall, we believe that our finding is general to both GTPases and ATPases. In addition to the biologically relevant finding, we also provide a critical discussion of the requirements from reliable surfaces for enzymatic reactions.