Reactivity of Phosphate Monoester Monoanions in Aqueous Solution. 2. A Theoretical Study of the Elusive Zwitterion Intermediates RO+(H)PO32-

Mechanistic Insights into Promoted Hydrolysis of Phosphoester Bonds by MoO2Cl2(DMF)2

A phosphoester bond is a crucial structural block in biological systems, whose occurrence is regulated by phosphatases. Molybdenum compounds have been reported to be active in phosphate ester hydrolysis of model phosphates. Specifically, MoO₂Cl₂(DMF)₂ is active in the hydrolysis of para-nitrophenyl phosphate (pNPP), leading to heteropolyoxometalate structures. We use density functional theory (DFT) to clarify the mechanism by which these species promote the hydrolysis of the phosphoester bond. The present calculations give insight into several key aspects of this reaction: (i) the speciation of this complex prior to interaction with the phosphate (DMF release, Mo-Cl hydrolysis, and pH influence on the speciation), (ii) the competition between phosphate addition and the molybdate nucleation process, (iii) and the mechanisms by which some plausible active species promote this hydrolysis in different conditions. We described thoroughly two different pathways depending on the nucleation possibilities of the molybdenum complex: one mononuclear mechanism, which is preferred in conditions in which very low complex concentrations are used, and another dinuclear mechanism, which is preferred at higher concentrations.

Inductive and External Electric Field Effects in Pentacoordinated Phosphorus Compounds

Pentacoordination at phosphorus is associated with a nucleophilic displacement reaction at tetracoordinated phosphorus compounds and shows a great variability in what respects their geometrical and energetic features. By means of a systematic theoretical study on a series of elementary model compounds, we have analyzed the bonding features. The pentacoordinated phosphorus compounds are held together by dative bonds, and the geometry and stability depends on the inductive effects originated by different substitutes at phosphorus. We show also that an external electric field can modify the geometrical features and the reactivity of the nucleophilic substitution reactions. This issue may have great interest in biological reactions involving pentacoordinated phosphorus where the electric field originated by the folded protein could influence the catalytic process. We report also additional calculations on the geometry and NMR spectra on three triphenyl phosphonium ylide derivatives, and our results compare well with the experimental data.

Benchmarking of DFT Functionals for the Hydrolysis of Phosphodiester Bonds

Phosphodiester bonds are an important chemical component of biological systems, and their hydrolysis and formation reactions are involved in major steps throughout metabolic pathways of all organisms. In this work, we applied dimethylphosphate as a model for this kind of bonds and calculated the potential energy surface for its hydrolysis at the approximated CCSD(T)/CBS//B3LYP/6-311++G(2d,2p) level. By varying the nucleophile (water or hydroxide) and the medium (vacuum or aqueous implicit solvent) we obtained and described four reaction paths. These structures were then used in a DFT functional benchmarking in which we tested a total of 52 functionals. Furthermore, the performances of HF, MP2, MP3, MP4, and CCSD were also evaluated. This benchmarking showed that MPWB1K, MPW1B95, and PBE1PBE are the more accurate functionals to calculate the energies of dimethylphosphate hydrolysis as far as activation and reaction energies are concerned. If considering only the activation energies, MPWB1K, MPW1B95, and B1B95 give the lowest errors when comparing to CCSD(T). A basis set benchmarking on the same system shows that 6-311+G(2d,2p) is the best basis set concerning the relationship between computational time and accuracy. We believe that our results will be of great help to further studies on related phosphodiester systems. This includes not only pure chemical problems but also biochemical studies in which DNA, RNA, and phospholipids are required to be depicted at a quantum level.

An implicit solvent model for SCC-DFTB with Charge-Dependent Radii

Motivated by the need of rapidly exploring the potential energy surface of chemical reactions that involve highly charged species, we have developed an implicit solvent model for the approximate density functional theory, SCC-DFTB. The solvation free energy is calculated using the popular model that employs Poisson-Boltzmann for electrostatics and a surface-area term for non-polar contributions. To balance the treatment of species with different charge distributions, we make the atomic radii that define the dielectric boundary and solute cavity depend on the solute charge distribution. Specifically, the atomic radii are assumed to be linearly dependent on the Mulliken charges and solved self-consistently together with the solute electronic structure. Benchmark calculations indicate that the model leads to solvation free energies of comparable accuracy to the SM6 model (especially for ions), which requires much more expensive DFT calculations. With analytical first derivatives and favorable computational speed, the SCC-DFTB based solvation model can be effectively used, in conjunction with high-level QM calculations, to explore the mechanism of solution reactions. This is illustrated with a brief analysis of the hydrolysis of mono-methyl mono-phosphate ester (MMP) and tri-methyl mono-phosphate ester (TMP). Possible future improvements are also briefly discussed.

Hydrolysis of aspartic acid phosphoramidate nucleotides: a comparative quantum chemical study

L-Aspartic acid has recently been found to be a good leaving group during HIV reverse transcriptase catalyzed incorporation of deoxyadenosine monophosphate (dAMP) in DNA. This showed that L-Asp is a good mimic for the pyrophosphate moiety of deoxyadenosine triphosphate. The present work explores the thermochemistry and mechanism for hydrolysis of several models for L-aspartic-dAMP using B3LYP/DGDZVP, MP2/6-311++G** and G3MP2 level of theory. The effect of the new compound is gradually investigated: starting from a simple methyl amine leaving group up to the aspartic acid leaving group. The enzymatic environment was mimicked by involving two Mg(2+) ions and some important active site residues in the reaction. All reactions are compared to the corresponding O-coupled leaving group, which is methanol for methyl amine and malic acid for aspartic acid. With methyl amine as a leaving group a tautomeric associative or tautomeric dissociative mechanism is preferred and the barrier is lower than the comparable mechanism with methanol as a leaving group. The calculations on the aspartic acid in the enzymatic environment show that qualitatively the mechanism is the same as for triphosphate but the barrier for hydrolysis by the associative mechanism is higher for L-aspartic-dAMP than for L-malic-dAMP and pyrophosphate.

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.

Description of phosphate hydrolysis reactions with the Self-Consistent-Charge Density-Functional-Tight-Binding (SCC-DFTB) theory. 1. Parameterization

Phosphate chemistry is involved in many key biological processes yet the underlying mechanism often remains unclear. For theoretical analysis to effectively complement experimental mechanistic analysis, it is essential to develop computational methods that can capture the complexity of the underlying potential energy surface and allow for sufficient sampling of the configurational space. To this end, we report the parameterization of an approximate density functional theory, Self-Consistent-Charge Density-Functional Tight-Binding (SCC-DFTB) method for systems containing phosphorus. Compared to high-level density functional theory and ab initio (MP2 and G3B3) results, the standard second-order parameterization is shown to give reliable structures for a diverse set of phosphate compounds but inaccurate energetics. With the on-site third-order terms included, referred to as SCC-DFTBPA, calculated proton affinities of phosphate compounds are substantially improved, although it remains difficult to obtain reliable proton affinity for both phosphates and compounds that do not contain phosphorus, indicating that further improvement in the formulation of SCC-DFTB is still a challenge to meet. To make SCC-DFTB applicable to phosphate reactions in the current (on-site-third-order-only) formulation, a "reaction-specific" parameterization, referred to as SCC-DFTBPR, is developed based on hydrolysis reactions of model phosphate species. Benchmark calculations in both the gas-phase and solution-phase indicate that SCC-DFTBPR gives reliable structural properties and semi-quantitative energetics for phosphate hydrolysis reactions. Since the number of reaction-specific parameters is small, it is likely that SCC-DFTBPR is applicable to a broad set of phosphate species. Indeed, for 56 reaction exothermicities and 47 energy barriers related to RNA catalysis model reactions collected from the QCRNA database, which involve molecules rather different from those used to parameterize SCC-DFTBPR, the corresponding root-mean-square difference between SCC-DFTBPR and high-level DFT results is only 5.3 kcal/mol. We hope that the parameterized SCC-DFTB models will complement NDDO based reaction-specific models (e.g., AM1-d/PhoT) and high-level ab initio QM/MM methods in better understanding the mechanism of phosphate chemistry in condensed phase, particularly biological systems.

Implications of protonation and substituent effects for C-O and O-P bond cleavage in phosphate monoesters

A recent study of phosphate monoesters that broke down exclusively through C-O bond cleavage and whose reactivity was unaffected by protonation of the nonbridging oxygens (Byczynski et al. J. Am. Chem. Soc. 2003, 125, 12541) raised several questions about the reactivity of phosphate monoesters, R-O-P(i). Potential catalytic strategies, particularly with regard to selectively promoting C-O or O-P bond cleavage, were investigated computationally through simple alkyl and aryl phosphate monoesters. Both C-O and O-P bonds lengthened upon protonating the bridging oxygen, R-O(H(+))-P(i), and heterolytic bond dissociation energies, DeltaH(C)(-)(O) and DeltaH(O)(-)(P), decreased. Which bond will break depends on the protonation state of the phosphoryl moiety, P(i), and the identity of the organosubstituent, R. Protonating the bridging oxygen when the nonbridging oxygens were already protonated favored C-O cleavage, while protonating the bridging oxygen of the dianion form, R-O-PO(3)(2)(-), favored O-P cleavage. Alkyl R groups capable of forming stable cations were more prone to C-O bond cleavage, with tBu > iPr > F(2)iPr > Me. The lack of effect on the C-O cleavage rate from protonating nonbridging oxygens could arise from two precisely offsetting effects: Protonating nonbridging oxygens lengthens the C-O bond, making it more reactive, but also decreases the bridging oxygen proton affinity, making it less likely to be protonated and, therefore, less reactive. The lack of effect could also arise without bridging oxygen protonation if the ratio of rate constants with different protonation states precisely matched the ratio of acidity constants, K(a). Calculations used hybrid density functional theory (B3PW91/6-31++G) methods with a conductor-like polarizable continuum model (CPCM) of solvation. Calculations on Me-phosphate using MP2/aug-cc-pVDZ and PBE0/aug-cc-pVDZ levels of theory, and variations on the solvation model, confirmed the reproducibility with different computational models.

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.

Supermolecule density functional calculations suggest a key role for solvent in alkaline hydrolysis of p-nitrophenyl phosphate

Supermolecule density functional theory calculations show that solvent is responsible for the concerted transition state in alkaline hydrolysis of p-nitrophenyl phosphate suggested by heavy atom kinetic isotope effects.

Mechanism of Triphosphate Hydrolysis in Aqueous Solution: QM/MM Simulations in Water Clusters

The mechanism of the hydrolysis reaction of the unprotonated methyl triphosphate (MTP) ester in water clusters has been modeled. The effective fragment potential based quantum mechanical-molecular mechanical (QM/MM) approach has been applied in the simulations. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step dissociative-type process similar to the case of the guanosine triphosphate (GTP) hydrolysis in the Ras-GAP protein complex (Grigorenko, B. L.; Nemukhin, A. V.; Topol, I. A.; Cachau, R. E.; Burt, S. K. Proteins: Struct., Funct., Bioinf. 2005, 60, 495). At the first stage, a unified action of environmental molecular groups and the catalytic water molecule leads to a substantial spatial separation of the gamma-phosphate group from the rest of the molecule. At the second stage, inorganic phosphate H2PO4- is formed from water and the metaphosphate anion PO3- through the chain of proton transfers along hydrogen bonds. The estimated activation barriers for MTP in aqueous solution at both stages (20 and 14 kcal/mol) are substantially higher than the corresponding barriers for the GTP hydrolysis in the protein.

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.

Theoretical Studies of Dissociative Phosphoryl Transfer in Interconversion of Phosphoenolpyruvate to Phosphonopyruvate: Solvent Effects, Thio Effects, and Implications for Enzymatic Reactions

The conversion of phosphoenolpyruvate (PEP) to phosphonopyruvate (P-pyr) is catalyzed by PEP mutase via a dissociative mechanism. In this work, we investigate the uncatalyzed reaction using ab initio methods, density functional theory, and the semiempirical MNDO/d method. Comparisons of geometries and relative energies of stationary points (minima and transition states) with density functional results indicate that the semiempirical method is reasonably accurate. Solvent effects are examined using implicit solvent models, including the recently extended smooth conductor-like screening model. Due to the large negative charge carried by the system, solvation is found to drastically alter the location and energy of stationary points along the dissociative reaction pathways. The influence of substituting a nonbridging phosphoryl oxygen by sulfur (thio effects) was also investigated. Implications of these results for the enzymatic reaction are discussed.

Quantum chemical modeling of the GTP hydrolysis by the RAS-GAP protein complex

We present results of the modeling for the hydrolysis reaction of guanosine triphosphate (GTP) in the RAS-GAP protein complex using essentially ab initio quantum chemistry methods. One of the approaches considers a supermolecular cluster composed of 150 atoms at a consistent quantum level. Another is a hybrid QM/MM method based on the effective fragment potential technique, which describes interactions between quantum and molecular mechanical subsystems at the ab initio level of the theory. Our results show that the GTP hydrolysis in the RAS-GAP protein complex can be modeled by a substrate-assisted catalytic mechanism. We can locate a configuration on the top of the barrier corresponding to the transition state of the hydrolysis reaction such that the straightforward descents from this point lead either to reactants GTP+H(2)O or to products guanosine diphosphate (GDP)+H(2)PO(4)(-). However, in all calculations such a single-step process is characterized by an activation barrier that is too high. Another possibility is a two-step reaction consistent with formation of an intermediate. Here the Pgamma-O(Pbeta) bond is already broken, but the lytic water molecule is still in the pre-reactive state. We present arguments favoring the assumption that the first step of the GTP hydrolysis reaction in the RAS-GAP protein complex may be assigned to the breaking of the Pgamma-O(Pbeta) bond prior to the creation of the inorganic phosphate.

Probing the Transition-State Structure of Dual-Specificity Protein Phosphatases Using a Physiological Substrate Mimic

Dual-specificity phosphatases (DSPs) belong to the large family of protein tyrosine phosphatases that contain the active-site motif (H/V)CxxGxxR(S/T), but unlike the tyrosine-specific enzymes, DSPs are able to catalyze the efficient hydrolysis of both phosphotyrosine and phosphoserine/threonine found on signaling proteins, as well as a variety of small-molecule aryl and alkyl phosphates. It is unclear how DSPs accomplish similar reaction rates for phosphoesters, whose reactivity (i.e., pK(a) of the leaving group) can vary by more than 10(8). Here, we utilize the alkyl phosphate m-nitrobenzyl phosphate (mNBP), leaving-group pK(a) = 14.9, as a physiological substrate mimic to probe the mechanism and transition state of the DSP, Vaccinia H1-related (VHR). Detailed pH and kinetic isotope effects of the V/K value for mNBP indicates that VHR reacts with the phosphate dianion of mNBP and that the nonbridge phosphate oxygen atoms are unprotonated in the transition state. (18)O and solvent isotope effects indicate differences in the respective timing of the proton transfer to the leaving group and P-O fission; with the alkyl ester substrate, protonation is ahead of P-O fission, while with the aryl substrate, the two processes are more synchronous. Kinetic analysis of the general-acid mutant D92N with mNBP was consistent with the requirement of Asp-92 in protonating the ester oxygen, either in a step prior to significant P-O bond cleavage or in a concerted but asynchronous mechanism in which protonation is ahead of P-O bond fission. Collectively, the data indicate that VHR and likely all DSPs can match leaving-group potential with the timing of the proton transfer to the ester oxygen, such that diverse aryl and alkyl phosphoesters are turned over with similar catalytic efficiency.

Transition State Differences in Hydrolysis Reactions of Alkyl versus Aryl Phosphate Monoester Monoanions

Although aryl phosphates have been the subject of numerous experimental studies, far less data bearing on the mechanism and transition states for alkyl phosphate reactions have been presented. Except for esters with very good leaving groups such as 2,4-dinitrophenol, the monoanion of phosphate esters is more reactive than the dianion. Several mechanisms have been proposed for the hydrolysis of the monoanion species. (18)O kinetic isotope effects in the nonbridging oxygen atoms and in the P-O(R) ester bond, and solvent deuterium isotope effects, have been measured for the hydrolysis of m-nitrobenzyl phosphate. The results rule out a proposed mechanism in which the phosphoryl group deprotonates water and then undergoes attack by hydroxide. The results are most consistent with a preequilibrium proton transfer from the phosphoryl group to the ester oxygen atom, followed by rate-limiting P-O bond fission, as originally proposed by Kirby and co-workers in 1967. The transition state for m-nitrobenzyl phosphate (leaving group pK(a) 14.9) exhibits much less P-O bond fission than the reaction of the more labile p-nitrophenyl phosphate (leaving group pK(a) = 7.14). This seemingly anti-Hammond behavior results from weakening of the P-O(R) ester bond resulting from protonation, an effect which calculations have shown is much more pronounced for aryl phosphates than for alkyl ones.