Theoretical studies on the hydrolysis of mono-phosphate and tri-phosphate in gas phase and aqueous solution

Prebiotic chemical reactivity in solution with quantum accuracy and microsecond sampling using neural network potentials

While RNA appears as a good candidate for the first autocatalytic systems preceding the emergence of modern life, the synthesis of RNA oligonucleotides without enzymes remains challenging. Because the uncatalyzed reaction is extremely slow, experimental studies bring limited and indirect information on the reaction mechanism, the nature of which remains debated. Here, we develop neural network potentials (NNPs) to study the phosphoester bond formation in water. While NNPs are becoming routinely applied to nonreactive systems or simple reactions, we demonstrate how they can systematically be trained to explore the reaction phase space for complex reactions involving several proton transfers and exchanges of heavy atoms. We then propagate at moderate computational cost hundreds of nanoseconds of a variety of enhanced sampling simulations with quantum accuracy in explicit solvent conditions. The thermodynamically preferred reaction pathway is a concerted, dissociative mechanism, with the transient formation of a metaphosphate transition state and direct participation of water solvent molecules that facilitate the exchange of protons through the nonbridging phosphate oxygens. Associative-dissociative pathways, characterized by a much tighter pentacoordinated phosphate, are higher in free energy. Our simulations also suggest that diprotonated phosphate, whose reactivity is never directly assessed in the experiments, is significantly less reactive than the monoprotonated species, suggesting that it is probably never the reactive species in normal pH conditions. These observations rationalize unexplained experimental results and the temperature dependence of the reaction rate, and they pave the way for the design of more efficient abiotic catalysts and activating groups.

Molecular Mechanisms of Phosphoester Bond Formation in Water Using Tight-Binding Ab Initio Molecular Dynamics

Phosphate groups are ubiquitous in biomolecules and are usually incorporated through phosphoester bonds between alcohol groups and orthophosphate. The formation of this bond is exceptionally difficult, with associated barriers of 30-45 kcal/mol in the absence of catalysts. In abiotic conditions, polymerizing nucleic acids without enzymes remains very challenging and is still a partly unsolved problem that severely questions the RNA World hypothesis for the origins of life. Offering a solution to this problem would involve a detailed knowledge of the reaction energetics and mechanisms, yet these remain not fully understood at a molecular level, especially because of the very slow reaction rates that represent a significant challenge for the experiments. The number of involved reaction coordinates and the possible role of the solvent in assisting the reaction are challenging for computational studies. Here, we use extensive ab initio molecular dynamics simulations using semiempirical tight-binding methods and enhanced sampling to address these issues. We first show that the choice of the tight-binding method is greatly limited by the instability of the water liquid phase for most DFTB generations and parameter sets that are widely available. We then focus on a model reaction involving methanol and orthophosphate, for which the two protonation states (mono- and dianionic) that are dominant around neutral pH are considered. We compare different proton coordinates that enable (or not) the participation of solvent water molecules. Our simulations suggest that in all cases, a dissociative associative mechanism, with an intermediate metaphosphate, is favored. The main difference between the two phosphate species is that reaction with the monoanion is assisted by the substrate, while that with the dianion involves solvent water molecules. Our results are in agreement with early experimental measurements, but the reaction barriers are underestimated in our framework. We believe that our approach provides an interesting perspective on how to sample the reaction phase space efficiently, but it calls for future studies using more accurate descriptions of chemical reactivity.

Phosphoryl-Transfer Reaction in RNA under Alkaline Conditions

The phosphoryl-transfer reaction in RNA under alkaline conditions by exploring the influence of several solvents theoretically was studied. The calculations were carried out by using the M06-2X functional and the solvents were taken as a continuum by using the solvent model density (SMD) method. The main findings show that the O2'-P-O5' angle in the reactants, the free activation energies, and the reaction mechanism are clearly dependent on the dielectric constant of the environment, thus showing that the electrostatic term is the determining factor for this chemical system with two negative charges. Our study seems to indicate that water, the solvent with the greatest dielectric constant, would be the solvent that increases the reaction rate the most. As this outcome was not the case in enzymatic catalysis, one has to conclude that, in the case of proteins as well as for ribozymes, the enzymatic catalysis is not mainly due to the solvent reaction field, but to local electrical fields as a result of enzyme preorganization.

How Mg2+ ions lower the SN2@P barrier in enzymatic triphosphate hydrolysis

Our quantum chemical activation strain analyses demonstrate how Mg2+ lowers the barrier of the enzymatic triphosphate hydrolysis through two distinct mechanisms: (a) weakening of the leaving-group bond, thereby decreasing activation strain; and (b) transition state (TS) stabilization through enhanced electrophilicity of the triphosphate PPP substrate, thereby strengthening the interaction with the nucleophile.

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 Insights on the Mechanism of the GTP Hydrolysis Catalyzed by the Elongation Factor Tu (EF-Tu)

The purpose of this work is to have a better understanding of the mechanism of GTP hydrolysis catalyzed by the elongation factor Tu. Two main aspects are being discussed in the literature: the associative or dissociative character of the process and the nature of nucleophile activation. The calculations of the QM subsystem have been done by means of the M06-2X density functional and the split valence triple-ζ 6-311+G(d,p) basis set. The environmental effect has been introduced through the continuum SMD method. We have studied three models of increasing complexity in order to analyze the different factors that intervene in the catalytic action. The results obtained in this paper confirm that the protonated His84 plays a fundamental role in the catalytic mechanism, but we have also found that the crystallographic sodium ion has a notable effect in the catalysis. So, our work has permitted a new insight, complementary to those obtained with QM/MM calculations, into this very complex process.

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.

What is “hypermobile” water?: detected in alkali halide, adenosine phosphate, and F-actin solutions by high-resolution microwave dielectric spectroscopy

Experimental observation by high-resolution microwave dielectric spectroscopy of hydration properties of alkali halide ions, adenosine phosphate ions, and F-actin revealed the existence of hypermobile water (HMW) molecules around those solutes. To understand the molecular process of HMW, two theoretical approaches are reviewed here. One is based on a statistical mechanical approach to analyze the rotational freedom of water molecules around a charged particle. Another approach reports direct calculation of dielectric relaxation process of water molecules around an ion. Experimentally observed HMW molecules are theoretically explained with the significance of multi-correlations among an ion and water molecules.

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.

Addressing open questions about phosphate hydrolysis pathways by careful free energy mapping

The nature and mechanism of phosphate hydrolysis reactions are of great interest in view of the crucial role of these reactions in key biological processes. Although it is becoming clearer that the ultimate way of resolving mechanistic controversies must involve reliable theoretical studies, it is not widely realized that such studies cannot be performed at present by using most existing automated ways and that only careful systematic studies can lead to meaningful conclusions. The present work clarifies the above point by considering the hydrolysis of phosphate monoesters. The clarification starts by defining the actual issues that should be addressed in careful studies and by highlighting the problems with studies that ignore the need for unique mechanistic definitions (e.g., works that confuse associative and dissociative pathways). We then focus on the analysis of the proton transfer (PT) pathways in phosphate hydrolysis and on recent suggestions that PT involves more than one water molecule. Here we point out that most of the studies that found a proton transfer through several water molecules have not involved a sufficient systematic search of the relevant reaction coordinates. This includes both energy minimization approaches as well as a recent metadynamics (MTD) simulation study. To illustrate the crucial need of exploring the potential surfaces reliably, rather than relying on automated approaches, we present here a very careful study of the free energy landscape along a 3D reaction coordinate (RC) exploring both the standard 2D RC, comprised of the attacking and leaving group reaction coordinates, as well as of the proton transfer (PT) coordinate. Our study points out that QM/MM minimization or MTD studies that concluded that the hydrolysis of phosphate monoesters involves a PT through several water molecules, have not explored carefully the single water (1W) path (that involves a direct PT form the attacking water molecule to the phosphate oxygen). Furthermore, we identified the most likely reason for the difficulty in finding the 1W path by QM/MM minimization methods, as well as by the current MTD simulations. We also discuss the problems with current studies that challenge the phosphate as a base mechanism and emphasize that all recent studies found associative/concerted paths (although many have not realized the meaning of their results). Finally, although we clearly do not have the last word about the 1W versus 2W paths we believe that we illustrated that the crucial mechanistic problems with alternative pathways should not be resolved by just running black box search approaches.

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.

Hydration properties of adenosine phosphate series as studied by microwave dielectric spectroscopy

Hydration properties of adenine nucleotides and orthophosphate (Pi) in aqueous solutions adjusted to pH=8 with NaOH were studied by high-resolution microwave dielectric relaxation (DR) spectroscopy at 20 °C. The dielectric spectra were analyzed using a mixture theory combined with a least-squares Debye decomposition method. Solutions of Pi and adenine nucleotides showed qualitatively similar dielectric properties described by two Debye components. One component was characterized by a relaxation frequency (f(c)=18.8-19.7 GHz) significantly higher than that of bulk water (17 GHz) and the other by a much lower f(c) (6.4-7.6 GHz), which are referred to here as hyper-mobile water and constrained water, respectively. By contrast, a hydration shell of only the latter type was found for adenosine (f(c)~6.7 GHz). The present results indicate that phosphoryl groups are mostly responsible for affecting the structure of the water surrounding the adenine nucleotides by forming one constrained water layer and an additional three or four layers of hyper-mobile water.

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.

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.

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.

Associative versus dissociative mechanisms of phosphate monoester hydrolysis: on the interpretation of activation entropies

Phosphate monoester and anhydride hydrolysis is ubiquitous in biology, being involved in, amongst other things, signal transduction, energy production, and the regulation of protein function. Therefore, this reaction has understandably been the focus of intensive research. Nevertheless, the precise mechanism by which phosphate monoester hydrolysis proceeds remains controversial. Traditionally, it has been assumed and frequently implied that a near-zero activation entropy is indicative of a dissociative pathway. Herein, we examine free-energy surfaces for the hydrolysis of the methyl phosphate dianion and the methyl pyrophosphate trianion in aqueous solution. In both cases, the reaction can proceed through either compact or expansive concerted (A(N)D(N)) transition states, with fairly similar barriers. We have evaluated the activation entropies for each transition state and demonstrate that both associative and dissociative transition states have near-zero entropies of activation that are in good agreement with experimental values. Therefore, we believe that the activation entropy alone is not a useful diagnostic tool, as it depends not only on bond orders at the transition state, but also on other issues that include (but are not limited to) steric factors determining the configurational volumes available to reactants during the reaction, solvation and desolvation effects that may be associated with charge redistribution upon approaching the transition state and entropy changes associated with intramolecular degrees of freedom as the transition state is approached.

Mechanism of the chemical step for the guanosine triphosphate (GTP) hydrolysis catalyzed by elongation factor Tu

Elongation factor Tu (EF-Tu), the protein responsible for delivering aminoacyl-tRNAs (aa-tRNAs) to ribosomal A site during translation, belongs to the group of guanosine-nucleotide (GTP/GDP) binding proteins. Its active 'on'-state corresponds to the GTP-bound form, while the inactive 'off'-state corresponds to the GDP-bound form. In this work we focus on the chemical step, GTP+H(2)O-->GDP+Pi, of the hydrolysis mechanism. We apply molecular modeling tools including molecular dynamics simulations and the combined quantum mechanical-molecular mechanical calculations for estimates of reaction energy profiles for two possible arrangements of switch II regions of EF-Tu. In the first case we presumably mimic binding of the ternary complex EF-Tu.GTP.aa-tRNA to the ribosome and allow the histidine (His85) side chain of the protein to approach the reaction active site. In the second case, corresponding to the GTP hydrolysis by EF-Tu alone, the side chain of His85 stays away from the active site, and the chemical reaction GTP+H(2)O-->GDP+Pi proceeds without participation of the histidine but through water molecules. In agreement with the experimental observations which distinguish rate constants for the fast chemical reaction in EF-Tu.GTP.aa-tRNA.ribosome and the slow spontaneous GTP hydrolysis in EF-Tu, we show that the activation energy barrier for the first scenario is considerably lower compared to that of the second case.

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.

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.

Gas-Phase Chemistry of Diphosphate Anions as a Tool To Investigate the Intrinsic Requirements of Phosphate Ester Enzymatic Reactions: The [M1M2HP2O7]− Ions

Experimental studies on gaseous inorganic phosphate ions are practically nonexistent, yet they can prove helpful for a better understanding of the mechanisms of phosphate ester enzymatic processes. The present contribution extends our previous investigations on the gas-phase ion chemistry of diphosphate species to the [M(1)M(2)HP(2)O(7)](-) ions where M(1) and M(2) are the same or different and correspond to the Li, Na, K, Cs, and Rb cations. The diphosphate ions are formed by electrospray ionization of 10(-4) M solutions of Na(5)P(3)O(10) in CH(3)CN/H(2)O (1/1) and MOH bases or M salts as a source of M(+) cations. The joint application of mass spectrometric techniques and quantum-mechanical calculations makes it possible to characterize the gaseous [M(1)M(2)HP(2)O(7)](-) ions as a mixed ionic population formed by two isomeric species: linear diphosphate anion coordinated to two M(+) cations (group I) and [PO(3)M(1)M(2)HPO(4)](-) clusters (group II). The relative gas-phase stabilities and activation barriers for the isomerization I-->II, which depend on the nature of the M(+) cations, highlight the electronic susceptibility of P-O-P bond breaking in the active site of enzymes. The previously unexplored gas-phase reactivity of [M(1)M(2)HP(2)O(7)](-) ions towards alcohols of different acidity was investigated by Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR/MS). The reaction proceeds by addition of the alcohol molecule followed by elimination of a water molecule.

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.

Reactivity of Metaphosphate and Thiometaphosphate in Water: A DFT Study

Metaphosphate is known to be highly reactive to water, whereas thiometaphosphate is relatively stable in aqueous solution. The difference in their reactivity has important mechanistic implications in interpreting the "thio effect" in phosphoryl transfer reactions. In this work, density functional theory is used to investigate the reactivity of both metaphosphate and its thio-substitute in their complexes with one, two, and three waters, and in aqueous solution. Barrier heights for converting metaphosphate to orthophosphate have been determined by geometry optimization. The results confirm that metaphosphate is consistently more reactive than thiometaphosphate and the activation free energy for both species decreases with the number of water molecules. The relative stability of thiometaphosphate is attributed to its less positively charged phosphorus atom.

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.

Polarization spectroscopy of gaseous tropolone in a strong electric field

We report studies of polarization spectroscopy of gaseous tropolone in a strong electric field using resonantly enhanced multiphoton ionization. The electric field induces localization of the tunneling proton between the two equivalent oxygen atoms. As a result, the C2v symmetry of the molecular frame is broken, and the parity selection rule is violated. The field induced transitions are type A with transition dipoles perpendicular to those under field free conditions. The polarization ratios, i.e., the ratios of the overall excitation yield under different polarizations of the resonant laser, thus deviate from those of a pure type B transition. In a field of 60 kV/cm, the experimental polarization ratio implies an essentially equal mixture of type B and type A transitions. Moreover, the induced transitions overlap with the two field-free subbands, and the resulting intensity ratios between the two subbands demonstrate dependence on the applied electric field. These observations can be qualitatively modeled using a quantum mechanical approach by assuming a two level system. A puzzling result is the magnitude of the transition dipole of the induced transition, which is proven to be essentially linearly dependent on the applied electric field.

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.