The Substrate-Assisted General Base Catalysis Model for Phosphate Monoester Hydrolysis: Evaluation Using Reactivity Comparisons

Solvent and solvation effects on reactivities and mechanisms of phospho group transfers from phosphate and phosphinate esters to nucleophiles

Organophosphorus esters fulfil many industrial, agricultural, and household roles. Nature has deployed phosphates and their related anhydrides as energy carriers and reservoirs, as constituents of genetic materials in the form of DNA and RNA, and as intermediates in key biochemical conversions. The transfer of the phosphoryl (PO₃) group is thus a ubiquitous biological process that is involved in a variety of transformations at the cellular level such as bioenergy and signals transductions. Significant attention has been paid in the last seven decades to understanding the mechanisms of uncatalyzed (solution) chemistry of the phospho group transfer because of the notion that enzymes convert the dissociative transition state structures in the uncatalyzed reactions into associative ones in the biological processes. In this regard, it has also been proposed that the rate enhancements enacted by enzymes result from the desolvation of the ground state in the hydrophobic active site environments, although theoretical calculations seem to disagree with this position. As a result, some attention has been paid to the study of the effects of solvent change, from water to less polar solvents, in uncatalyzed phospho transfer reactions. Such changes have consequences on the stabilities of the ground and the transition states of reactions which affect reactivities and, sometimes, the mechanisms of reactions. This review seeks to collate and evaluate what is known about solvent effects in this domain, especially their effects on rates of reactions of different classes of organophosphorus esters. The outcome of this exercise shows that a systematized study of solvent effects needs to be undertaken to fully understand the physical organic chemistry of the transfer of phosphates and related molecules from aqueous to substantially hydrophobic environments, since significant knowledge gaps exist.

Synthesis, Stability, and Kinetics of Hydrogen Sulfide Release of Dithiophosphates

The development of chemicals to slowly release hydrogen sulfide would aid the survival of plants under environmental stressors as well as increase harvest yields. We report a series of dialkyldithiophosphates and disulfidedithiophosphates that slowly degrade to release hydrogen sulfide in the presence of water. Kinetics of the degradation of these chemicals were obtained at 85 °C and room temperature, and it was shown that the identity of the alkyl or sulfide group had a large impact on the rate of hydrolysis, and the rate constant varied by more than 10⁴×. For example, using tert-butanol as the nucleophile yielded a dithiophosphate (8) that hydrolyzed 13,750× faster than the dithiophosphate synthesized from n-butanol (1), indicating that the rate of hydrolysis is structure-dependent. The rates of hydrolysis at 85 °C varied from a low value of 6.9 × 10^-4 h^-1 to a high value of 14.1 h^-1. Hydrogen sulfide release in water was also quantified using a hydrogen sulfide-sensitive electrode. Corn was grown on an industrial scale and dosed with dibutyldithiophosphate to show that these dithiophosphates have potential applications in agriculture. At a loading of 2 kg per acre, a 6.4% increase in the harvest yield of corn was observed.

Challenges and advances in the computational modeling of biological phosphate hydrolysis

Phosphate ester hydrolysis is fundamental to many life processes, and has been the topic of substantial experimental and computational research effort. However, even the simplest of phosphate esters can be hydrolyzed through multiple possible pathways that can be difficult to distinguish between, either experimentally, or computationally. Therefore, the mechanisms of both the enzymatic and non-enzymatic reactions have been historically controversial. In the present contribution, we highlight a number of technical issues involved in reliably modeling these computationally challenging reactions, as well as proposing potential solutions. We also showcase examples of our own work in this area, discussing both the non-enzymatic reaction in aqueous solution, as well as insights obtained from the computational modeling of organophosphate hydrolysis and catalytic promiscuity amongst enzymes that catalyze phosphoryl transfer.

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.

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.

Catalytic Mechanism of the Maltose Transporter Hydrolyzing ATP

We use quantum mechanical and molecular mechanical (QM/MM) simulations to study ATP hydrolysis catalyzed by the maltose transporter. This protein is a prototypical member of a large family that consists of ATP-binding cassette (ABC) transporters. The ABC proteins catalyze ATP hydrolysis to perform a variety of biological functions. Despite extensive research efforts, the precise molecular mechanism of ATP hydrolysis catalyzed by the ABC enzymes remains elusive. In this work, the reaction pathway for ATP hydrolysis in the maltose transporter is evaluated using a QM/MM implementation of the nudged elastic band method without presuming reaction coordinates. The potential of mean force along the reaction pathway is obtained with an activation free energy of 19.2 kcal/mol in agreement with experiments. The results demonstrate that the reaction proceeds via a dissociative-like pathway with a trigonal bipyramidal transition state in which the cleavage of the γ-phosphate P-O bond occurs and the O-H bond of the lytic water molecule is not yet broken. Our calculations clearly show that the Walker B glutamate as well as the switch histidine stabilizes the transition state via electrostatic interactions rather than serving as a catalytic base. The results are consistent with biochemical and structural experiments, providing novel insight into the molecular mechanism of ATP hydrolysis in the ABC proteins.

Catalytic Mechanism of Mammalian Adenylyl Cyclase: A Computational Investigation

Adenylyl cyclase (AC) catalyzes the synthesis of cyclic AMP, an important intracellular regulatory molecule, from ATP. We propose a catalytic mechanism for class III mammalian AC based on density functional theory calculations. We employ a model of the AC active site derived from a crystal structure of mammalian AC activated by Gα·GTP and forskolin at separate allosteric sites. We compared the calculated activation free energies for 13 possible reaction sequences involving proton transfer, nucleophilic attack, and elimination of pyrophosphate. The proposed most probable mechanism is initiated by deprotonation of 3'OH and water-mediated transfer of the 3'H to the γ-phosphate. Proton transfer is followed by changes in coordination of the two magnesium ion cofactors and changes in the conformation of ATP to enhance the role of 3'O as a nucleophile and to bring 3'O close to Pα. The subsequent phosphoryl transfer step is concerted and rate-limiting. Comparison of the enzyme-catalyzed and nonenzymatic reactions reveals that the active site residues lower the free energy barrier for both phosphoryl transfer and proton transfer and significantly shift the proton transfer equilibrium. Calculations for mutants K1065A and R1029A demonstrate that K1065 plays a significant role in shifting the proton transfer equilibrium, whereas R1029 is important for making the transition state of concerted phosphoryl transfer tight rather than loose.

Structure and reactivity of phosphate diesters. Dependence on the nonleaving group

The hydrolytic reactivity of simple phosphate diesters with very good leaving groups is known to be practically independent of the nonleaving group at 100 °C. Calculations on the (too slow to measure) hydrolysis at 25 °C of a series of p-nitrophenyl diesters ROPO2––OpNP with a wide range of nonleaving group OR indicate a small but significant effect at the lower temperature, making the R = methyl ester the most reactive. This is in the opposite sense to the much larger effect observed for the reactions of triesters and consistent with a reaction driven primarily by leaving group departure. The calculations use a continuum model with up to four discrete water molecules: two or three waters give the best results.

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.

A DFT study on proton transfers in hydrolysis reactions of phosphate dianion and sulfate monoanion

B3LYP calculations were carried out on hydrolysis reactions of monosubstituted(R) phosphate dianion and sulfate monoanion. In the reacting system, water clusters (H2O)22 and (H2O)35 are included to trace reaction paths. For both P and S substrates with R = methyl group, elementary processes were calculated. While the phosphate undergoes the substitution at the phosphorus, the sulfate does at the methyl carbon. For the S substrate with R = neopentyl group, the product tert-amyl alcohol was found to be formed via a dyotropic rearrangement from the neopentyl alcohol intermediate. For R = aryl groups, transition-state geometries were calculated to be similar between P and S substrates. Calculated activation energies are in good agreement with experimental values. After the rate-determining transition state of the substitution, the hydronium ion H3O(+) is formed at the third water molecule. It was suggested that alkyl and aryl substrates are of the different reactivity of the hydrolysis.

Catalytic mechanism of bacteriophage T4 Rad50 ATP hydrolysis

Spontaneous double-strand breaks (DSBs) are one of the most deleterious forms of DNA damage, and their improper repair can lead to cellular dysfunction. The Mre11 and Rad50 proteins, a nuclease and an ATPase, respectively, form a well-conserved complex that is involved in the initial processing of DSBs. Here we examine the kinetic and catalytic mechanism of ATP hydrolysis by T4 Rad50 (gp46) in the presence and absence of Mre11 (gp47) and DNA. Single-turnover and pre-steady state kinetics on the wild-type protein indicate that the rate-limiting step for Rad50, the MR complex, and the MR-DNA complex is either chemistry or a conformational change prior to catalysis. Pre-steady state product release kinetics, coupled with viscosity steady state kinetics, also supports that the binding of DNA to the MR complex does not alter the rate-limiting step. The lack of a positive deuterium solvent isotope effect for the wild type and several active site mutants, combined with pH-rate profiles, implies that chemistry is rate-limiting and the ATPase mechanism proceeds via an asymmetric, dissociative-like transition state. Mutation of the Walker A/B and H-loop residues also affects the allosteric communication between Rad50 active sites, suggesting possible routes for cooperativity between the ATP active sites.

Phosphoryl transfers of the phospholipase D superfamily: a quantum mechanical theoretical study

The HKD-containing Phospholipase D superfamily catalyzes the cleavage of the headgroup of phosphatidylcholine to produce phosphatidic acid and choline. The mechanism of this cleavage process is studied theoretically. The geometric basis of our models is the X-ray crystal structure of the five-coordinate phosphohistidine intermediate from Streptomyces sp . Strain PMF (PDB Code = 1V0Y ). Hybrid ONIOM QM:QM methodology with Density Functional Theory (DFT) and semiempirical PM6 (DFT:PM6) is used to acquire thermodynamic and kinetic data for the initial phosphoryl transfer, subsequent hydrolysis, and finally, the formation of the experimentally observed ″dead-end″ phosphohistidine product (PDB Code = 1V0W ). The model contains nineteen amino acid residues (including the two highly conserved HKD-motifs), four explicit water molecules, and the substrate. Via computations, the persistence of the short-lived five-coordinate phosphorane intermediate on the minutes times scale is rationalized. This five-coordinate phosphohistidine intermediate energetically exists between the hydrolysis event and ″substrate reorganization″ (the reorganization of the in vitro model substrate within the active site). Computations directly support the thermodynamic favorability of the in vitro four-coordinate phosphohistidine product. In vivo, the activation energy of substrate reorganization is too high, perhaps due to a combination of substrate immobility when embedded in the lipid bilayer, as well as its larger steric bulk compared to the compound used in the in vitro substrate soaks. On this longer time scale, the enzyme will migrate along the lipid membrane toward its next substrate target, rather than promote the formation of the dead-end product.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Efficient, crosswise catalytic promiscuity among enzymes that catalyze phosphoryl transfer

The observation that one enzyme can accelerate several chemically distinct reactions was at one time surprising because the enormous efficiency of catalysis was often seen as inextricably linked to specialization for one reaction. Originally underreported, and considered a quirk rather than a fundamental property, enzyme promiscuity is now understood to be important as a springboard for adaptive evolution. Owing to the large number of promiscuous enzymes that have been identified over the last decade, and the increased appreciation for promiscuity's evolutionary importance, the focus of research has shifted to developing a better understanding of the mechanistic basis for promiscuity and the origins of tolerant or restrictive specificity. We review the evidence for widespread crosswise promiscuity amongst enzymes that catalyze phosphoryl transfer, including several members of the alkaline phosphatase superfamily, where large rate accelerations between 10(6) and 10(17) are observed for both native and multiple promiscuous reactions. This article is part of a Special Issue entitled: Chemistry and mechanism of phosphatases, diesterases and triesterases.

Mechanistic insights into the hydrolysis of a nucleoside triphosphate model in neutral and acidic solution

Nucleoside triphosphate hydrolysis is an essential component of all living systems. Despite extensive research, the exact modus and mechanism of this ubiquitous reaction still remain elusive. In this work, we examined the detailed hydrolysis mechanisms of a model nucleoside triphosphate in acidic and neutral solution by means of ab initio simulations. The timescale of the reaction was accessed through use of an accelerated sampling method, metadynamics. Both hydrolyses were found to proceed via different mechanisms; the acidic system reacted by means of concerted general acid catalysis (found to be a so-called D(N)A(N)A(H)D(xh) mechanism), whereas the neutral system reacted by way of a different mechanism (namely, D(N)*A(N)D(xh)A(H)). A neighboring water molecule took on the role of a general base in both systems, which has not been seen before but is a highly plausible reaction path, meaning that substrate-assisted catalysis was not observed in the bulk water environment.

Examining the promiscuous phosphatase activity of Pseudomonas aeruginosa arylsulfatase: a comparison to analogous phosphatases

Pseudomonas aeruginosa arylsulfatase (PAS) is a bacterial sulfatase capable of hydrolyzing a range of sulfate esters. Recently, it has been demonstrated to also show very high proficiency for phosphate ester hydrolysis. Such proficient catalytic promiscuity is significant, as promiscuity has been suggested to play an important role in enzyme evolution. Additionally, a comparative study of the hydrolyses of the p-nitrophenyl phosphate and sulfate monoesters in aqueous solution has demonstrated that despite superficial similarities, the two reactions proceed through markedly different transition states with very different solvation effects, indicating that the requirements for the efficient catalysis of the two reactions by an enzyme will also be very different (and yet they are both catalyzed by the same active site). This work explores the promiscuous phosphomonoesterase activity of PAS. Specifically, we have investigated the identity of the most likely base for the initial activation of the unusual formylglycine hydrate nucleophile (which is common to many sulfatases), and demonstrate that a concerted substrate-as-base mechanism is fully consistent with the experimentally observed data. This is very similar to other related systems, and suggests that, as far as the phosphomonoesterase activity of PAS is concerned, the sulfatase behaves like a "classical" phosphatase, despite the fact that such a mechanism is unlikely to be available to the native substrate (based on pK(a) considerations and studies of model systems). Understanding such catalytic versatility can be used to design novel artificial enzymes that are far more proficient than the current generation of designer enzymes.

Insights into the phosphoryl transfer mechanism of cyclin-dependent protein kinases from ab initio QM/MM free-energy studies

Phosphorylation reactions catalyzed by kinases and phosphatases play an indispensible role in cellular signaling, and their malfunctioning is implicated in many diseases. A better understanding of the catalytic mechanism will help design novel and effective mechanism-based inhibitors of these enzymes. In this work, ab initio quantum mechanical/molecular mechanical studies are reported for the phosphoryl transfer reaction catalyzed by a cyclin-dependent kinase, CDK2. Our results suggest that an active-site Asp residue, rather than ATP as previously proposed, serves as the general base to activate the Ser nucleophile. The corresponding transition state features a dissociative, metaphosphate-like structure, stabilized by the Mg(2+) ion and several hydrogen bonds. The calculated free-energy barrier is consistent with experimental values. Implications of our results in this and other protein kinases are discussed.

Density functional calculations on the effect of sulfur substitution for 2'-hydroxypropyl-p-nitrophenyl phosphate: C-O vs. P-O bond cleavage

Density functional theory calculations have been used to investigate the intra-molecular attack of 2'-hydroxypropyl-p-nitrophenyl phosphate (HPpNP) and its analogous compound 2-thiouridyl-p-nitrophenyl phosphate (s-2'pNP). Bulk solvent effect has been tested at the geometry optimization level with the polarized continuum model. It is found that the P-path involving the intra-molecular attack at the phosphorus atom and C-path involving the attack at the beta carbon atom proceed through the S(N)2-type mechanism for HPpNP and s-2'pNP. The calculated results indicate that the P-path with the free energy barrier of about 11 kcal/mol is more accessible than the C-path for the intra-molecular attack of HPpNP, which favors the formation of the five-membered phosphate diester. While for s-2'pNP, the C-path with the free energy barrier of about 21 kcal/mol proceeds more favorably than the P-path. The calculated energy barriers of the favorable pathways for HPpNP and s-2'pNP are both in agreement with the experimental results.

Protein kinase biochemistry and drug discovery

Protein kinases are fascinating biological catalysts with a rapidly expanding knowledge base, a growing appreciation in cell regulatory control, and an ascendant role in successful therapeutic intervention. To better understand protein kinases, the molecular underpinnings of phosphoryl group transfer, protein phosphorylation, and inhibitor interactions are examined. This analysis begins with a survey of phosphate group and phosphoprotein properties which provide context to the evolutionary selection of phosphorylation as a central mechanism for biological regulation of most cellular processes. Next, the kinetic and catalytic mechanisms of protein kinases are examined with respect to model aqueous systems to define the elements of catalysis. A brief structural biology overview further delves into the molecular basis of catalysis and regulation of catalytic activity. Concomitant with a prominent role in normal physiology, protein kinases have important roles in the disease state. To facilitate effective kinase drug discovery, classic and emerging approaches for characterizing kinase inhibitors are evaluated including biochemical assay design, inhibitor mechanism of action analysis, and proper kinetic treatment of irreversible inhibitors. As the resulting protein kinase inhibitors can modulate intended and unintended targets, profiling methods are discussed which can illuminate a more complete range of an inhibitor's biological activities to enable more meaningful cellular studies and more effective clinical studies. Taken as a whole, a wealth of protein kinase biochemistry knowledge is available, yet it is clear that a substantial extent of our understanding in this field remains to be discovered which should yield many new opportunities for therapeutic intervention.

The effect of leaving group on mechanistic preference in phosphate monoester hydrolysis

We present 2-dimensional potential energy surfaces and optimised transition states (TS) for water attack on a series of substituted phosphate monoester monoanions at the DFT level of theory, comparing a standard 6-31++g(d,p) basis set with a larger triple-zeta (augmented cc-pVTZ) basis set. Small fluorinated model compounds are used to simulate increasing leaving group stability without adding further geometrical complexity to the system. We demonstrate that whilst changing the leaving group causes little qualitative change in the potential energy surfaces (with the exception of the system with the most electron withdrawing leaving group, CF(3)O(-), in which the associative pathway changes from a stepwise A(N) + D(N) pathway to a concerted A(N)D(N) pathway), there is a quantitative change in relative gas-phase and solution barriers for the two competing pathways. In line with previous studies, in the case of OCH(3), the barriers for the associative and dissociative pathways are similar in solution, and the two pathways are equally viable and indistinguishable in solution. However, significantly increasing the stability of the leaving group (decreasing proton affinity, PA) results in the progressive favouring of a stepwise dissociative, D(N) + A(N), mechanism over associative mechanisms.

Alkaline hydrolysis of ethylene phosphate: an ab initio study by supermolecule model and polarizable continuum approach

The alkaline hydrolysis reaction of ethylene phosphate (EP) has been investigated using a supermolecule model, in which several explicit water molecules are included. The structures and single-point energies for all of the stationary points are calculated in the gas phase and in solution at the B3LYP/6-31++G(df,p) and MP2/6-311++G(df,2p) levels. The effect of water bulk solvent is introduced by the polarizable continuum model (PCM). Water attack and hydroxide attack pathways are taken into account for the alkaline hydrolysis of EP. An associative mechanism is observed for both of the two pathways with a kinetically insignificant intermediate. The water attack pathway involves a water molecule attacking and a proton transfer from the attacking water to the hydroxide in the first step, followed by an endocyclic bond cleavage to the leaving group. While in the first step of the hydroxide attack pathway the nucleophile is the hydroxide anion. The calculated barriers in aqueous solution for the water attack and hydroxide attack pathways are all about 22 kcal/mol. The excellent agreement between the calculated and observed values demonstrates that both of the two pathways are possible for the alkaline hydrolysis of EP.

Nucleotide recognition and phosphate linkage hydrolysis at a lipid cubic interface

Mononucleotides, when entrapped within a mono-olein-based cubic Ia3d liquid crystalline phase, have been found to undergo hydrolysis at the sugar-phosphate ester bond in spite of their natural inertness toward hydrolysis. Here, kinetics of the hydrolysis reaction and interactions between the lipid matrix and the mononucleotide adenosine 5'-monophosphate disodium salt (AMP) and its 2'-deoxy derivative (dAMP) are thoroughly investigated in order to shed some light on the mechanism of the nucleotide recognition and phosphate ester hydrolysis. Experiments evidenced that molecular recognition occurs essentially through the sn-2 and the sn-3 alcoholic OH groups of mono-olein. As deduced from the apparent activation energies, the mechanism underlying the hydrolysis reaction is the same for AMP and dAMP. Nevertheless, the reaction proceeds slower for the latter, highlighting a substantial difference in the chemical behavior of the two nucleotides. A model that explains the hydrolysis reaction is presented. Remarkably, the hydrolysis mechanism appears to be highly specific for the Ia3d phase.

Are mixed explicit/implicit solvation models reliable for studying phosphate hydrolysis? A comparative study of continuum, explicit and mixed solvation models

Phosphate hydrolysis is ubiquitous in biology. However, despite intensive research on this class of reactions, the precise nature of the reaction mechanism remains controversial. Herein, we have examined the hydrolysis of three homologous phosphate diesters. The solvation free energy was simulated by means of either an implicit solvation model (COSMO), hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) or a mixed solvation model in which N water molecules were explicitly included in the ab initio description of the reacting system (where N=1-3), with the remainder of the solvent being implicitly modelled as a continuum. Here, both COSMO and QM/MM-FEP reproduce DeltaG(obs) within an error of about 1 kcal mol(-1). However, we demonstrate that in order to obtain any kind of reliable results from a mixed model, it is essential to carefully select the explicit water molecules from short QM/MM runs that act as a model for the true infinite system. Additionally, the mixed models tend to be increasingly unstable and miss larger entropic contributions as more explicit water molecules are placed into the system. Thus, our analysis indicates that this approach provides an unreliable way for modelling phosphate hydrolysis in solution.

Does water relay play an important role in phosphoryl transfer reactions? Insights from theoretical study of a model reaction in water and tert-butanol

To investigate whether water relay plays an important role in phosphoryl transfer reactions, we have used several theoretical approaches to compare key properties of uridine 3'-m-nitrobenzyl phosphate (UNP) in aqueous and tert-butanol solutions. Previous kinetic experiments found that the isomerization reaction of UNP is abolished in tert-butanol, which was interpreted as the direct evidence that supports the role of water relay in phosphoryl transfer. We have analyzed solute flexibility and solvent structure near the solute using equilibrium molecular dynamics simulations and a combined quantum mechanical/molecular mechanism (QM/MM) potential function for the solute. Snapshots from the simulations are then used in minimum energy path calculations to compare the energetics of direct nucleophilic attack and water-mediated nucleophilic attack pathways. QM/MM simulations are also used to compare the pseudorotation barriers for the pentavalent intermediate formed following the nucleophilic attack, another key step for the isomerization reaction. Combined results from these calculations suggest that water relay does not offer any significant energetic advantage over the direct nucleophilic attack. Unfortunately, the lack of isomerization in tert-butanol solution cannot be straightforwardly explained based on the results we have obtained here and therefore requires additional analysis. This study, nevertheless, has provided new insights into several most commonly discussed possibilities.

Dineopentyl phosphate hydrolysis: evidence for stepwise water attack

Phosphate ester hydrolysis is ubiquitous in biology, playing a central role in energy production, signaling, biosynthesis, and the regulation of protein function among other things. Although the mechanism of action of the enzymes regulating this reaction has been the focus of intensive research in the past few decades, the correct description of this apparently simple reaction remains controversial. A clear understanding of the mechanism that takes place in solution is crucial to be able to evaluate whether proposals for the enzyme-catalyzed mechanisms are reasonable. For the pH-independent hydrolysis of phosphate diesters, several kinetically equivalent mechanisms are plausible, including hydroxide attack on the neutral phosphate. However, it is very difficult to measure the rate of this reaction directly by experimental methods, so it has been evaluated by examining the rate of hydrolysis of neutral phosphate triesters, where a methyl group has replaced a proton. This may not be an accurate model of the neutral phosphate diester and does not provide information about a reaction pathway that is concerted with nucleophilic attack to generate a similar phosphorane. We have carefully mapped out free energy surfaces for both hydroxide and water attack on the dineopentyl phosphate anion and for water attack on the neutral diester. In doing so, we have accurately reproduced existing experimental data and demonstrate that water attack proceeds through an associative mechanism with proton transfer to the phosphate to generate a phosphorane intermediate. Our data show that the substrate-as-base mechanism is viable for phosphate ester hydrolysis, which may have important implications for the studies of phosphate ester hydrolysis by enzymes.

The case for the intermediacy of monomeric metaphosphate analogues during oxidation of H-phosphonothioate, H-phosphonodithioate, and H-phosphonoselenoate monoesters: mechanistic and synthetic studies

Studies on the reaction of H-phosphonothioate, H-phosphonodithioate, and H-phosphonoselenoate monoesters with iodine in the presence of a base led to identification of a unique oxidation pathway, which consists of the initial oxidation of the sulfur or selenium atom in these compounds, followed by oxidative elimination of hydrogen iodide to generate the corresponding metaphosphate analogues. The intermediacy of the latter species during oxidation of the investigated H-phosphonate monoester derivatives with iodine was supported by various diagnostic experiments. The scope and limitation of these oxidative transformations for the purpose of the synthesis of nucleoside phosphorothioate, nucleoside phosphorodithioate, and nucleoside phosphoroselenoate diesters was also investigated.

On the interpretation of the observed linear free energy relationship in phosphate hydrolysis: a thorough computational study of phosphate diester hydrolysis in solution

The hydrolysis of phosphate esters is crucially important to biological systems, being involved in, among other things, signaling, energy transduction, biosynthesis, and the regulation of protein function. Despite this, there are many questions that remain unanswered in this important field, particularly with regard to the preferred mechanism of hydrolysis of phosphate esters, which can proceed through any of multiple pathways that are either associative or dissociative in nature. Previous comparisons of calculated and observed linear free energy relationships (LFERs) for phosphate monoester dianions with different leaving groups showed that the TS character gradually changes from associative to dissociative with the increasing acidity of the leaving group, while reproducing the experimental LFER. Here, we have generated ab initio potential energy surfaces for the hydrolysis of phosphate diesters in solution, with a variety of leaving groups. Once again, the reaction changes from a compact concerted pathway to one that is more expansive in character when the acidity of the leaving group increases. When such systems are examined in solution, it is essential to take into consideration the contribution of solute to the overall activation entropy, which remains a major computational challenge. The popular method of calculating the entropy using a quasi-harmonic approximation appears to markedly overestimate the configurational entropy for systems with multiple occupied energy wells. We introduce an improved restraint release approach for evaluating configurational entropies and apply this approach to our systems. We demonstrate that when this factor is taken into account, it is possible to reproduce the experimental LFER for this system with reasonable accuracy.

Theoretical study of specific hydrogen-bonding effects on the bridging P-OR bond strength of phosphate monoester dianions

It has been proposed that the driving force for the initial phosphoryl transfer step of protein tyrosine phosphatases (PTPases) could be activation of the substrate ROPO32- by means of an enforced hydrogen-bonding interaction between an aspartic general acid and the bridging oxygen atom O (Zhang et al. Biochemistry 1995, 34, 16088-16096). The potential catalytic effect of this type of interaction, with regard to P-OR bond cleavage, was investigated computationally through simple model systems in which an efficient intramolecular hydrogen bond can take place between a H-bond donor group and the bridging oxygen atom of the dianionic phosphate. The dielectric effect of the environment (epsilon = 1, 4, and 78) was also explored. The results indicate that this interaction causes significant lengthenings of the scissile P-OR bond in all media but with more extreme effects observed in the low dielectric fields epsilon = 1 and epsilon = 4. It is interesting that, in all cases examined, this interaction actually contributes to stabilize the reactant state while causing its P-OR bond to lengthen. Overall, our results support the idea that this specific hydrogen-bonding situation might well be used by PTPases as an important driving force for promoting phosphoryl transfer reactions through highly dissociative transition states.

Kinetics and Mechanism of the Aminolysis of Aryl Phenyl Chlorothiophosphates with Anilines

Kinetic studies of the reactions of aryl phenyl chlorothiophosphates (1) and aryl 4-chlorophenyl chlorothiophosphates (2) with substituted anilines in acetonitrile at 55.0 degrees C are reported. The negative values of the cross-interaction constant rhoXY (rhoXY = -0.22 and -0.50 for 1 and 2, respectively) between substituents in the nucleophile (X) and substrate (Y) indicate that the reactions proceed by concerted SN2 mechanism. The primary kinetic isotope effects (kH/kD = 1.11-1.13 and 1.10-1.46 for 1 and 2, respectively) involving deuterated aniline nucleophiles are obtained. Front- and back-side nucleophilic attack on the substrates is proposed mainly on the basis of the primary kinetic isotope effects. A hydrogen-bonded, four-center-type transition state is suggested for a front-side attack, while the trigonal bipyramidal pentacoordinate transition state is suggested for a back-side attack. The MO theoretical calculations of the model reactions of dimethyl chlorothiophosphate (1') and dimethyl chlorophosphate (3') with ammonia are carried out. Considering the specific solvation effect, the front-side nucleophilic attack can occur competitively with the back-side attack in the reaction of 1'.

Theoretical Evaluation of the Substrate-Assisted Catalysis Mechanism for the Hydrolysis of Phosphate Monoester Dianions

Quantum chemistry methods coupled with a continuum solvation model have been applied to evaluate the substrate-assisted catalysis (SAC) mechanism recently proposed for the hydrolysis of phosphate monoester dianions. The SAC mechanism, in which a proton from the nucleophile is transferred to a nonbridging phosphoryl oxygen atom of the substrate prior to attack, has been proposed in opposition to the widely accepted mechanism of direct nucleophilic reaction. We have assessed the SAC proposal for the hydrolysis of three representative phosphate monoester dianions (2,4-dinitrophenyl phosphate, phenyl phosphate, and methyl phosphate) by considering the reactivity of the hydroxide ion toward the phosphorus center of the corresponding singly protonated monoesters. The reliability of the calculations was verified by comparing the calculated and the observed values of the activation free energies for the analogous S(N)2(P) reactions of F- with the monoanion of the monoester 2,4-dinitrophenyl phosphate and its diester analogue, methyl 2,4-dinitrophenyl phosphate. It was found that the orientation of the phosphate hydrogen atom has important implications with regard to the nature of the transition state. Hard nucleophiles such as OH- and F- can attack the phosphorus atom of a singly protonated phosphate monoester only if the phosphate hydrogen atom is oriented toward the leaving-group oxygen atom. As a result of this proton orientation, the SAC mechanism in solution is characterized by a small Brønsted coefficient value (beta(lg)=-0.25). This mechanism is unlikely to apply to aryl phosphates, but becomes a likely possibility for alkyl phosphate esters. If oxyanionic nucleophiles of pK(a)<11 are involved, as in alkaline phosphatase, then the S(N)2(P) reaction may proceed with the phosphate hydrogen atom oriented toward the nucleophile. In this situation, a large negative value of beta(lg) (-0.95) is predicted for the substrate-assisted catalysis mechanism.

On the mechanism of hydrolysis of phosphate monoesters dianions in solutions and proteins

The nature of the hydrolysis of phosphate monoester dianions in solutions and in proteins is a problem of significant current interest. The present work explores this problem by systematic calculations of the potential surfaces of the reactions of a series of phosphate monoesters with different leaving groups. These calculations involve computational studies ranging from ab initio calculations with implicit solvent models to ab initio QM/MM free energy calculations. The calculations reproduce the observed linear free energy relationship (LFER) for the solution reaction and thus are consistent with the overall experimental trend and can be used to explore the nature of the transition state (TS) region, which is not accessible to direct experimental studies. It is found that the potential surface for the associative and dissociative paths is very flat and that the relative height of the associative and dissociative TS is different in different systems. In general, the character of the TS changes from associative to dissociative upon decrease in the pKa of the leaving group. It is also demonstrated that traditional experimental markers such as isotope effects and the LFER slope cannot be used in a conclusive way to distinguish between the two classes of transition states. In addition it is found that the effective charges of the TS do not follow the previously assumed simple rule. Armed with that experience we explore the free energy surface for the GTPase reaction of the RasGap system. In this case it is found that the surface is flat but that the lowest TS is associative. The present study indicates that the nature of the potential surfaces for the phosphoryl transfer reactions in solution and proteins is quite complicated and cannot be determined in a conclusive way without the use of careful theoretical studies that should, of course, reproduce the available experimental information.

Theoretical Studies of the Hydroxide-Catalyzed P−O Cleavage Reactions of Neutral Phosphate Triesters and Diesters in Aqueous Solution: Examination of the Changes Induced by H/Me Substitution

DFT calculations and dielectric continuum methods have been employed to map out the lowest activation free-energy profiles for the alkaline hydrolysis of representative phosphate triesters and diesters, including trimethyl phosphate (TMP), dimethyl 4-nitrophenyl phosphate (DMNPP), dimethyl hydrogen phosphate (DMHP), and the dimethyl phosphate anion (DMP-). The reliability of the calculations is supported by the excellent agreement observed between the calculated and the experimentally determined activation enthalpies for phosphate triesters with poor (TMP) and good (DMNPP) leaving groups. The results obtained for the OH- + DMHP and OH- + DMP- reactions are also consistent with all the available experimental information concerning the hydrolysis reaction of dimethyl phosphate anion at pH > 5. By performing geometry optimizations in the dielectric field (epsilon = 78.39), we found that OH- can attack the phosphorus atom of DMHP without capturing its proton only if the O-H bond of DMHP is oriented opposite the attacking OH- group. In these conditions, the rate for OH- attack on DMHP was found to be approximately 10(3)-fold faster than that for OH- attack on TMP. The calculated rate acceleration induced by the phosphoryl proton corresponds to the maximum rate effect expected from kinetic studies. Overall, our calculations performed on the dimethyl phosphate ester predict that, contrary to what is generally observed for RNA and aryl phosphodiesters, the water-promoted P-O cleavage reaction of DNA should dominate the base-catalyzed reaction at pH 7. These results are suggestive that nucleases may be less proficient as catalysts than has recently been suspected.

Protein kinase structure and function analysis with chemical tools

Protein kinases are the largest enzyme superfamily involved in cell signal transduction and represent therapeutic targets for a range of diseases. There have been intensive efforts from many labs to understand their catalytic mechanisms, discover inhibitors and discern their cellular functions. In this review, we will describe two approaches developed to analyze protein kinases: bisubstrate analog inhibition and phosphonate analog utilization. Both of these methods have been used in combination with the protein semisynthesis method expressed protein ligation to advance our understanding of kinase-substrate interactions and functional elucidation of phosphorylation. Previous work on the nature of the protein kinase mechanism suggests it follows a dissociative transition state. A bisubstrate analog was designed against the insulin receptor kinase to mimic the geometry of a dissociative transition state reaction coordinate distance. This bisubstrate compound proved to be a potent inhibitor against the insulin receptor kinase and occupied both peptide and nucleotide binding sites. Bisubstrate compounds with altered hydrogen bonding potential as well as varying spacers between the adenine and the peptide demonstrate the importance of the original design features. We have also shown that related bisubstrate analogs can be used to potently block serine/threonine kinases including protein kinase A. Since many protein kinases recognize folded protein substrates for efficient phosphorylation, it was advantageous to incorporate the peptide-ATP conjugates into protein structures. Using expressed protein ligation, a Src-ATP conjugate was produced and shown to be a high affinity ligand for the Csk tyrosine kinase. Nonhydrolyzable mimics of phosphoSer/phosphoTyr can be useful in examining the functionality of phosphorylation events. Using expressed protein ligation, we have employed phosphonomethylene phenylalanine and phosphonomethylene alanine to probe the phosphorylation of Tyr and Ser, respectively. These tools have permitted an analysis of the SH2-phosphatases (SHP1 and SHP2), revealing a novel intramolecular stimulation of catalytic activity mediated by the corresponding phosphorylation events. They have also been used to characterize the cellular regulation of the melatonin rhythm enzyme by phosphorylation.

Gaseous H5P2O8? Ions: A Theoretical and Experimental Study on the Hydrolysis and Synthesis of Diphosphate Ion

The structure and reactivity of gaseous H5P2O8- ions obtained from the chemical ionization (CI) of an H4P2O7/H2O mixture and from electrospray ionization (ESI) of CH3CN/H2O/H4P2O7 solutions were investigated by Fourier transform ion cyclotron (FTICR) and triple quadrupole mass spectrometry. Theoretical calculations performed at the B3LYP/6-31+G* level of theory and collisionally activated dissociation (CAD) mass spectrometric results allowed the ionic population obtained in the CI conditions to be structurally characterized as a mixture of gaseous [H3P2O7...H2O]-, [H3PO4...H2PO4]-, and [PO3...H3PO4...H2O]- clusters. The energy profile emerging from theoretical calculations affords insight into the mechanism of diphosphate ion hydrolysis and synthesis.

On the role of the conserved aspartate in the hydrolysis of the phosphocysteine intermediate of the low molecular weight tyrosine phosphatase

The usual rate-determining step in the catalytic mechanism of the low molecular weight tyrosine phosphatases involves the hydrolysis of a phosphocysteine intermediate. To explain this hydrolysis, general base-catalyzed attack of water by the anion of a conserved aspartic acid has sometimes been invoked. However, experimental measurements of solvent deuterium kinetic isotope effects for this enzyme do not reveal a rate-limiting proton transfer accompanying dephosphorylation. Moreover, base activation of water is difficult to reconcile with the known gas-phase proton affinities and solution phase pK(a)'s of aspartic acid and water. Alternatively, hydrolysis could proceed by a direct nucleophilic attack by a water molecule. To understand the hydrolysis mechanism, we have used high-level density functional methods of quantum chemistry combined with continuum electrostatics models of the protein and the solvent. Our calculations do not support a catalytic activation of water by the aspartate. Instead, they indicate that the water oxygen directly attacks the phosphorus, with the aspartate residue acting as a H-bond acceptor. In the transition state, the water protons are still bound to the oxygen. Beyond the transition state, the barrier to proton transfer to the base is greatly diminished; the aspartate can abstract a proton only after the transition state, a result consistent with experimental solvent isotope effects for this enzyme and with established precedents for phosphomonoester hydrolysis.

Kinetic isotope effects in Ras-catalyzed GTP hydrolysis: evidence for a loose transition state

A remote labeling method has been developed to determine (18)O kinetic isotope effects (KIEs) in Ras-catalyzed GTP hydrolysis. Substrate mixtures consist of (13)C-depleted GTP and [(18)O,(13)C]GTP that contains (18)O at phosphoryl positions of mechanistic interest and (13)C at all carbon positions of the guanosine moiety. Isotope ratios of the nonvolatile substrates and products are measured by using a chemical reaction interface/isotope ratio mass spectrometer. The isotope effects are 1.0012 (0.0026) in the gamma nonbridge oxygens, 1.0194 (0.0025) in the leaving group oxygens (the beta-gamma oxygen and the two beta nonbridge oxygens), and 1.0105 (0.0016) in the two beta nonbridge oxygens. The KIE in the beta-gamma bridge oxygen was computed to be 1.0116 or 1.0088 by two different methods. The significant KIE in the leaving group reveals that chemistry is largely rate-limiting whereas the KIEs in the gamma nonbridge oxygens and the leaving group indicate a loose transition state that approaches a metaphosphate. The KIE in the two beta nonbridge oxygens is roughly equal to that in the beta-gamma bridge oxygen. This indicates that, in the transition state, Ras shifts one-half of the negative charge that arises from P(gamma)-O(beta-gamma) fission from the beta-gamma bridge oxygen to the two beta nonbridge oxygens. The KIE effects, interpreted in light of structural and spectroscopic data, suggest that Ras promotes a loose transition state by stabilizing negative charge in the beta-gamma bridge and beta nonbridge oxygens of GTP.

Why does the Ras switch "break" by oncogenic mutations?

The elucidation of the structure of the RasGAP complex provides what is perhaps the most detailed link between protein structure and cancer causing mutations. In particular, it is known that mutations of Gln 61 destroy the GTPase activity of the complex, locks the cell in its ON state and thus, can cause cancer. It is entirely unclear however, why this specific mutation is so important. The present work uncovers the elusive role of Gln 61 by computer simulation of the GTPase reaction in Ras, RasGAP and of their mutants. Simulations of the effects of mutations of Gln 61 reproduce the corresponding observed changes in activation energies and allow us to analyze the energy contributions to these effects. It is found that Gln 61 does not operate in a direct chemical way nor by a direct electrostatic or steric interaction with the transition state (TS). Instead, oncogenic mutations of Gln 61 lead to the destruction of the exquisitely preorganized catalytic configuration of the active site of the RasGAP complex. This "allosteric" effect causes a major reduction in the electrostatic stabilization of the TS. Our findings have general relevance to other proteins that control signal transduction processes.