Theoretical examination of Mg(2+)-mediated hydrolysis of a phosphodiester linkage as proposed for the hammerhead ribozyme

The hammerhead ribozyme is an RNA molecule capable of self-cleavage at a unique site within its sequence. Hydrolysis of this phosphodiester linkage has been proposed to occur via an in-line attack geometry for nucleophilic displacement by the 2'-hydroxyl on the adjoining phosphorus to generate a 2',3'-cyclic phosphate ester with elimination of the 5'-hydroxyl group, requiring a divalent metal ion under physiological conditions. The proposed S(N)2(P) reaction mechanism was investigated using density functional theory calculations incorporating the hybrid functional B3LYP to study this metal ion-dependent reaction with a tetraaquo magnesium (II)-bound hydroxide ion. For the Mg(2+)-catalyzed reaction, the gas-phase geometry optimized calculations predict two transition states with a kinetically insignificant, yet clearly defined, pentacoordinate intermediate. The first transition state located for the reaction is characterized by internal nucleophilic attack coupled to proton transfer. The second transition state, the rate-determining step, involves breaking of the exocyclic P-O bond where a metal-ligated water molecule assists in the departure of the leaving group. These calculations demonstrate that the reaction mechanism incorporating a single metal ion, serving as a Lewis acid, functions as a general base and can afford the necessary stabilization to the leaving group by orienting a water molecule for catalysis.     

A DFT study on the formation of a phosphohistidine intermediate in prostatic acid phosphatase

Histidine phosphatases are a class of enzymes that are characterized by the presence of a conserved RHGXRXP motif. This motif contains a catalytic histidine that is being phosphorylated in the course of a dephosphorylation reaction catalyzed by these enzymes. Prostatic acid phosphatase (PAP) is one such enzyme. The dephosphorylation of phosphotyrosine by PAP is a two-step process. The first step involves the transfer of a phosphate group from the substrate to the histidine (His12). The present study reports on the details of the first step of this reaction, which was investigated using a series of quantum chemistry calculations. A number of quantum models were constructed containing various residues that were thought to play a role in the mechanism. In all these models, the transition state displayed an associative character. The transition state is stabilized by three active site arginines (Arg11, Arg15, and Arg79), two of which belong to the aforementioned conserved motif. The work also demonstrated that His12 could act as a nucleophile. The enzyme is further characterized by a His257-Asp258 motif. The role of Asp258 has been elusive. In this work, we propose that Asp258 acts as a proton donor which becomes protonated when the substrate enters the binding pocket. Evidence is also obtained that the transfer of a proton from Asp258 to the leaving group is possibly mediated by a water molecule in the active site. The work also underlines the importance of His257 in lowering the energy barrier for the nucleophilic attack.     

DNA polymerase beta catalysis: are different mechanisms possible?

DNA polymerases are crucial constituents of the complex cellular machinery for replicating and repairing DNA. Discerning mechanistic pathways of DNA polymerase on the atomic level is important for revealing the origin of fidelity discrimination. Mammalian DNA polymerase beta (pol beta), a small (39 kDa) member of the X-family, represents an excellent model system to investigate polymerase mechanisms. Here, we explore several feasible low-energy pathways of the nucleotide transfer reaction of pol beta for correct (according to Watson-Crick hydrogen bonding) G:C basepairing versus the incorrect G:G case within a consistent theoretical framework. We use mixed quantum mechanics/molecular mechanics (QM/MM) techniques in a constrained energy minimization protocol to effectively model not only the reactive core but also the influence of the rest of the enzymatic environment and explicit solvent on the reaction. The postulated pathways involve initial proton abstraction from the terminal DNA primer O3'H group, nucleophilic attack that extends the DNA primer chain, and elimination of pyrophosphate. In particular, we analyze several possible routes for the initial deprotonation step: (i) direct transfer to a phosphate oxygen O(Palpha) of the incoming nucleotide, (ii) direct transfer to an active site Asp group, and (iii) transfer to explicit water molecules. We find that the most probable initial step corresponds to step (iii), involving initial deprotonation to water, which is followed by proton migration to active site Asp residues, and finally to the leaving pyrophosphate group, with an activation energy of about 15 kcal/mol. We argue that initial deprotonation steps (i) and (ii) are less likely as they are at least 7 and 11 kcal/mol, respectively, higher in energy. Overall, the rate-determining step for both the correct and the incorrect nucleotide cases is the initial deprotonation in concert with nucleophilic attack at the phosphate center; however, the activation energy we obtain for the mismatched G:G case is 5 kcal/mol higher than that of the matched G:C complex, due to active site structural distortions. Taken together, our results support other reported mechanisms and help define a framework for interpreting nucleotide specificity differences across polymerase families, in terms of the concept of active site preorganization or the so-called "pre-chemistry avenue".     

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.     

Mechanism of intramolecular catalysis in the hydrolysis of alkyl monoesters of 1,8-naphthalic acid

Hydrolysis of alkyl 1,8-naphthalic acid monoesters 1a-d is subject to highly efficient intramolecular nucleophilic catalysis by the neighboring COOH group. The reactivity for the COOH reaction depends on the leaving group pK(a), with values of β(LG) of -0.50, consistent with a mechanism involving rate determining breakdown of tetrahedral addition intermediates. The release of the steric strain of the peri-substitiuents in the highly reactive alkyl 1,8-naphthalic acid monoesters is fundamental to understand the observed special reactivity in this intramolecular reaction. DFT calculations show how the proton transfers involved in the cleavage of the neutral ester can be catalyzed by solvent water, thus facilitating the departure of poor alkoxide leaving groups.     

Water-assisted reaction mechanism of monozinc beta-lactamases

Hybrid Car-Parrinello QM/MM calculations are used to investigate the reaction mechanism of hydrolysis of a common beta-lactam substrate (cefotaxime) by the monozinc beta-lactamase from Bacillus cereus (BcII). The calculations suggest a fundamental role for an active site water in the catalytic mechanism. This water molecule binds the zinc ion in the first step of the reaction, expanding the zinc coordination number and providing a proton donor adequately oriented for the second step. The free energy barriers of the two reaction steps are similar and consistent with the available experimental data. The conserved hydrogen bond network in the active site, defined by Asp120, Cys221, and His263, not only contributes to orient the nucleophile (as already proposed), but it also guides the second catalytic water molecule to the zinc ion after the substrate is bound. The hydrolysis reaction in water has a relatively high free energy barrier, which is consistent with the stability of cefotaxime in water solution. The modeled Michaelis complexes for other substrates are also characterized by the presence of an ordered water molecule in the same position, suggesting that this mechanism might be general for the hydrolysis of different beta-lactam substrates.     

New insights into the mechanism of the Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase from S. enterica: a theoretical study

The reaction pathway of Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase (DHQD) from S. enterica has been studied by performing molecular dynamics (MD) simulations and density functional theory (DFT) calculations and the corresponding potential energy profile has also been identified. On the basis of the results, the catalytic hydrolysis process for the wild-type enzyme consists of three major reaction steps, including nucleophilic attack on the carbon atom involved in the carbon-nitrogen double bond of the Schiff base intermediate by a water molecule, deprotonation of the His143 residue, and dissociation between the product and the Lys170 residue of the enzyme. The remarkable difference between this and the previously proposed reaction mechanism is that the second step here, absent in the previously proposed reaction mechanism, plays an important role in facilitating the reaction through a key proton transfer by the His143 residue, resulting in a lower energy barrier. Comparison with our recently reported results on the Schiff base formation and dehydration processes clearly shows that the Schiff base hydrolysis is rate-determining in the overall reaction catalyzed by type I DHQD, consistent with the experimental prediction, and the calculated energy barrier of ～16.0 kcal mol(-1) is in good agreement with the experimentally derived activation free energy of ～14.3 kcal mol(-1). When the imidazole group of His143 residue is missing, the Schiff base hydrolysis is initiated by a hydroxide ion in the solution, rather than a water molecule, and both the reaction mechanism and the kinetics of Schiff base hydrolysis have been remarkably changed, clearly elucidating the catalytic role of the His143 residue in the reaction. The new mechanistic insights obtained here will be valuable for the rational design of high-activity inhibitors of type I DHQD as non-toxic antimicrobials, anti-fungals, and herbicides.     

Dinuclear Zn(II) complex catalyzed phosphodiester cleavage proceeds via a concerted mechanism: a density functional theory study

Density functional theory (DFT) calculations were used to study the mechanism for the cleavage reaction of the RNA analogue HpPNP (HpPNP = 2-hydroxypropyl-4-nitrophenyl phosphate) catalyzed by the dinuclear Zn(II) complex of 1,3-bis(1,4,7-triazacyclonon-1-yl)-2-hydroxypropane (Zn(2)(L(2)O)). We present a binding mode in which each terminal phosphoryl oxygen atom binds to one zinc center, respectively, and the nucleophilic 2-hydroxypropyl group coordinates to one of the zinc ions, while the hydroxide from deprotonation of a water molecule coordinates to the other zinc ion. Our calculations found a concerted mechanism for the HpPNP cleavage with a 16.5 kcal/mol reaction barrier. An alternative proposed stepwise mechanism through a pentavalent oxyphosphorane dianion reaction intermediate for the HpPNP cleavage was found to be less feasible with a significantly higher energy barrier. In this stepwise mechanism, the deprotonation of the nucleophilic 2-hydroxypropyl group is accompanied with nucleophilic attack in the rate-determining step. Calculations of the nucleophile (18)O kinetic isotope effect (KIE) and leaving (18)O KIE for the concerted mechanism are in reasonably good agreement with the experimental values. Our results indicate a specific-base catalysis mechanism takes place in which the deprotonation of the nucleophilic 2-hydroxypropyl group occurs in a pre-equilibrium step followed by a nucleophilic attack on the phosphorus center. Detailed comparison of the geometric and electronic structure for the HpPNP cleavage reaction mechanisms in the presence/absence of catalyst revealed that the catalyst significantly altered the determining-step transition state to become far more associative or tight, that is, bond formation to the nucleophile was remarkably more advanced than leaving group bond fission in the catalyzed mechanism. Our results are consistent with and provide a reliable interpretation for the experimental observations that suggest the reaction occurs by a concerted mechanism (see Humphry, T.; Iyer, S.; Iranzo, O.; Morrow, J. R.; Richard, J. P.; Paneth, P.; Hengge, A. C. J. Am. Chem. Soc. 2008, 130, 17858-17866) and has a specific-base catalysis character (see Yang, M.-Y.; Iranzo, O.; Richard, J. P.; Morrow, J. R. J. Am. Chem. Soc. 2005, 127, 1064-1065).     

Reaction mechanism of deoxyribonucleotidase: a theoretical study

The reaction mechanism of human deoxyribonucleotidase (dN) is studied using high-level quantum-chemical methods. dN catalyzes the dephosphorylation of deoxyribonucleoside monophosphates (dNMPs) to their nucleoside form in human cells. Large quantum models are employed (99 atoms) based on a recent X-ray crystal structure. The calculations support the proposed mechanism in which Asp41 performs a nucleophilic attack on the phosphate to form a phospho-enzyme intermediate. Asp43 acts in the first step as an acid, protonating the leaving nucleoside, and in the second step as a base, deprotonating the lytic water. No pentacoordinated intermediates could be located.     

Reaction of imidazole with toluene-4-sulfonate salts of substituted phenyl N-methylpyridinium-4-carboxylate esters: special base catalysis by imidazole

The reaction of imidazole in aqueous solution with toluene-4-sulfonate salts of substituted phenyl N-methylpyridinium-4-carboxylate esters obeys the rate law: k(obs) - k(background) = k2[Im] + k3[Im]2 where [Im] is the imidazole concentration present as free base. The parameters k2 and k3 fit Br?nsted type free energy correlations against the pKa of the leaving phenol with betaLg values of -0.65 and -0.42 respectively. The imidazolysis is insensitive to catalysis by general bases and yet k3 for the 3-cyanophenyl ester possesses a deuterium oxide solvent isotope effect of 4.43 consistent with rate limiting proton transfer. A special catalytic function is proposed for decomposition of the tetrahedral addition intermediate (T+/-) via k3 whereby the catalytic imidazole interacts electrophilically with the leaving phenolate ion and removes a proton from the nitrogen in the rate limiting step with subsequent non-rate limiting ArO-C bond fission. This is consistent with the change in effective charge on the leaving oxygen in the transition structure of k3 which is more positive (-0.42) than that expected (-0.60) for the equilibrium formation of the zwitterion intermediate. The catalytic function at the leaving oxygen is likely to be an electrophilic role of the NH as a hydrogen bond donor. In the k2 step the deuterium oxide solvent isotope effect of 1.51 for the 3-cyanophenyl ester and the betaLg of -0.65 are consistent with rate limiting expulsion of the phenolate ion from the T+/- intermediate. The absence of general base catalysis of imidazolysis rules out the established mechanism for aminolysis of esters where T+/- is stabilised by a standard rate limiting proton transfer. The kinetically equivalent term for k3 where T- reacts with the imidazolium ion as an acid catalyst would require this step to be rate limiting and involve proton transfer not consistent with departure of the good aryl oxide leaving group.     

Antibiotic deactivation by a dizinc beta-lactamase: mechanistic insights from QM/MM and DFT studies

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods and density functional theory (DFT) were used to investigate the initial ring-opening step in the hydrolysis of moxalactam catalyzed by the dizinc L1 beta-lactamase from Stenotrophomonas maltophilia. Anchored at the enzyme active site via direct metal binding as suggested by a recent X-ray structure of an enzyme-product complex (Spencer, J.; et al. J. Am. Chem. Soc. 2005, 127, 14439), the substrate is well aligned with the nucleophilic hydroxide that bridges the two zinc ions. Both QM/MM and DFT results indicate that the addition of the hydroxide nucleophile to the carbonyl carbon in the substrate lactam ring leads to a metastable intermediate via a dominant nucleophilic addition barrier. The potential of mean force obtained by SCC-DFTB/MM simulations and corrected by DFT/MM calculations yields a reaction free energy barrier of 23.5 kcal/mol, in reasonable agreement with the experimental value of 18.5 kcal/mol derived from kcat of 0.15 s(-1). It is further shown that zinc-bound Asp120 plays an important role in aligning the nucleophile, but accepts the hydroxide proton only after the nucleophilic addition. The two zinc ions are found to participate intimately in the catalysis, consistent with the proposed mechanism. In particular, the Zn(1) ion is likely to serve as an "oxyanion hole" in stabilizing the carbonyl oxygen, while the Zn(2) ion acts as an electrophilic catalyst to stabilize the anionic nitrogen leaving group.     

Spontaneous dehydration mechanism of aromatic aldehyde reactions with hydroxyl and non‐hydroxyl amines

The plot of rate constants vs. pH for the dehydration step of the reaction between furfural and 5‐nitrofurfural with hydroxylamine, N‐methylhydroxylamine, and O‐methylhydroxylamine, shows two regions corresponding to the oxonium ion‐catalyzed and spontaneous dehydration. The oxonium ion‐catalyzed dehydration region of the reaction of furfural with the above mentioned hydroxylamines exhibits general acid catalysis with excellent Brønsted correlation (Brønsted coefficients: 0.76 (r = 0.986), 0.68 (r = 0.987), and 0.67 (r = 0.993) respectively). However, the rate constants of the spontaneous dehydration of these hydroxylamines, where water is considered the general acid catalyst, exhibit a large positive deviation from the Brønsted line. This fact was not observed in the reaction of non‐hydroxyl amines with different aromatic aldehydes by other authors, thus supporting that the spontaneous dehydration steps for these reactions proceed by intramolecular catalysis. The mechanism of intramolecular catalysis might be stepwise. First, a zwitterionic intermediate is formed. It can then evolve in the second step by loss of water, or follow a concerted pathway, with the transference of a proton through a five‐membered ring (general intramolecular acid catalysis). In the case of non‐hydroxyl amines, data suggested the possibility of a mechanism of intramolecular proton transfer through one or two water molecules, from the nitrogen of the amine to the leaving hydroxide ion. © 2002 Wiley Periodicals, Inc. Int J Chem Kinet 34: 685–692, 2002     

The Catalytic Mechanism of HIV-1 Integrase for DNA 3'-End Processing Established by QM/MM Calculations

The development of HIV-1 integrase (INT) inhibitors has been hampered by incomplete structural and mechanistic information. Despite the efforts made to overcome these limitations, only one compound has been approved for clinical use so far. In this work, we have used all experimental information available for INT and similar enzymes, to build a model of the holo-integrase:DNA complex that includes an entire central core domain, a ssDNA GCAGT substrate, and two magnesium ions. Subsequently, we used a large array of computational techniques, which included molecular dynamics, thermodynamic integration, and high-level quantum mechanics/molecular mechanics (QM/MM) calculations to study the possible pathways for the mechanism of 3' end processing catalyzed by INT. We found that the only viable mechanism to hydrolyze the DNA substrate is a nucleophilic attack of an active site water molecule to the phosphorus atom of the scissile phosphoester bond, with the attacking water being simultaneously deprotonated by an Mg(2+)-bound hydroxide ion. The unstable leaving oxoanion is protonated by an Mg(2+)-bound water molecule within the same elementary reaction step. This reaction has an activation free energy of 15.4 kcal/mol, well within the limits imposed by the experimental turnover. This work significantly improves the fundamental knowledge on the integrase chemistry. It can also contribute to the discovery of leads against HIV-1 infection as it provides, for the first time, accurate transition states structures that can be successfully used as templates for high-throughput screening of new INT inhibitors.     

Proton storage site in bacteriorhodopsin: new insights from quantum mechanics/molecular mechanics simulations of microscopic pK(a) and infrared spectra

Identifying the group that acts as the proton storage/loading site is a challenging but important problem for understanding the mechanism of proton pumping in biomolecular proton pumps, such as bacteriorhodopsin (bR) and cytochrome c oxidase. Recent experimental studies of bR propelled the idea that the proton storage/release group (PRG) in bR is not an amino acid but a water cluster embedded in the protein. We argue that this idea is at odds with our knowledge of protein electrostatics, since invoking the water cluster as the PRG would require the protein to raise the pK(a) of a hydronium by almost 11 pK(a) units, which is difficult considering known cases of pK(a) shifts in proteins. Our recent quantum mechanics/molecular mechanics (QM/MM) simulations suggested an alternative "intermolecular proton bond" model in which the stored proton is shared between two conserved Glu residues (194 and 204). Here we show that this model leads to microscopic pK(a) values consistent with available experimental data and the functional requirement of a PRG. Extensive QM/MM simulations also show that, independent of a number of technical issues, such as the influence of QM region size, starting X-ray structure, and nuclear quantum effects, the "intermolecular proton bond" model is qualitatively consistent with available spectroscopic data. Potential of mean force calculations show explicitly that the stored proton strongly prefers the pair of Glu residues over the water cluster. The results and analyses help highlight the importance of considering protein electrostatics and provide arguments for why the "intermolecular proton bond" model is likely applicable to the PRG in biomolecular proton pumps in general.     

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.     

Intracomplex general acid/base catalyzed cleavage of RNA phosphodiester bonds: the leaving group effect

The general acid/base catalyzed cleavage of a number of alkyl esters of uridine-3'- (and -5'-)phosphate has been studied by utilizing a cleaving agent, in which the catalytic moiety (a substituted 1,3,5-triazine) is tethered to an anchoring Zn(II):cyclen moiety. Around pH 7, formation of a strong ternary complex between uracil, Zn(II) and cyclen brings the general acid/base catalyst close to the scissile phosphodiester linkage, resulting in rate acceleration of 1-2 orders of magnitude with the uridine-3'-phosphodiesters. Curiously, no acceleration was observed with their 5'-counterparts. A β(lg) value of -0.7 has been determined for the general acid/base catalyzed cleavage, consistent with a proton transfer to the leaving group in the rate-limiting step.     

The role of the putative catalytic base in the phosphoryl transfer reaction in a protein kinase: first-principles calculations

Protein kinases are important enzymes controlling the majority of cellular signaling events via a transfer of the gamma-phosphate of ATP to a target protein. Even after many years of study, the mechanism of this reaction is still poorly understood. Among many factors that may be responsible for the 1011-fold rate enhancement due to this enzyme, the role of the conserved aspartate (Asp166) has been given special consideration. While the essential presence of Asp166 has been established by mutational studies, its function is still debated. The general base catalyst role assigned to Asp166 on the basis of its position in the active site has been brought into question by the pH dependence of the reaction rate, isotope measurements, and pre-steady-state kinetics. Recent semiempirical calculations have added to the controversy surrounding the role of Asp166 in the catalytic mechanism. No major role for Asp166 has been found in these calculations, which have predicted the reaction process consisting of an early transfer of a substrate proton onto the phosphate group. These conclusions were inconsistent with experimental observations. To address these differences between experimental results and theory with a more reliable computational approach and to provide a theoretical platform for understanding catalysis in this important enzyme family, we have carried out first-principles structural and dynamical calculations of the reaction process in cAPK kinase. To preserve the essential features of the reaction, representations of all of the key conserved residues (82 atoms) were included in the calculation. The structural calculations were performed using the local basis density functional (DFT) approach with both hybrid B3LYP and PBE96 generalized gradient approximations. This kind of calculation has been shown to yield highly accurate structural information for a large number of systems. The optimized reactant state structure is in good agreement with X-ray data. In contrast to semiempirical methods, the lowest energy product state places the substrate proton on Asp166. First-principles molecular dynamics simulations provide additional support for the stability of this product state. The latter also demonstrate that the proton transfer to Asp166 occurs at a point in the reaction where bond cleavage at the PO bridging position is already advanced. This mechanism is further supported by the calculated structure of the transition state in which the substrate hydroxyl group is largely intact. A metaphoshate-like structure is present in the transition state, which is consistent with the X-ray structures of transition state mimics. On the basis of the calculated structure of the transition state, it is estimated to be 85% dissociative. Our analysis also indicates an increase in the hydrogen bond strength between Asp166 and substrate hydroxyl and a small decrease in the bond strength of the latter in the transition state. In summary, our calculations demonstrate the importance of Asp166 in the enzymatic mechanism as a proton acceptor. However, the proton abstraction from the substrate occurs late in the reaction process. Thus, in the catalytic mechanism of cAPK protein kinase, Asp166 plays a role of a "proton trap" that locks the transferred phosphoryl group to the substrate. These results resolve prior inconsistencies between theory and experiment and bring new understanding of the role of Asp166 in the protein kinase catalytic mechanism.     

Hydrolysis of phosphotriesters: a theoretical analysis of the enzymatic and solution mechanisms

A theoretical study on the alkaline hydrolysis of paraoxon, one of the most popular organophosphorus pesticides, in aqueous solution and in the active site of Pseudomonas diminuta phosphotriesterase (PTE) is presented. Simulations by means of hybrid quantum mechanics/molecular mechanics (QM/MM) potentials show that the hydrolysis of paraoxon takes place through an A(N)D(N) or associative mechanism both in solution and in the active site of PTE. The results correctly reproduce the magnitude of the activation free energies and can be used to rationalize the observed kinetic isotope effects (KIEs) for the hydrolysis of paraoxon in both media. Enzymatic hydrolysis of O,O-diethyl p-chlorophenyl phosphate, a phosphotriester having a leaving group with higher pK(a) than paraoxon, was also simulated. Hydrolysis of this phosphotriester by PTE follows a A(N)+D(N) mechanism with a pentacoordinate intermediate. Moreover, the leaving group of this new substrate coordinates to one of the zinc ions of the bimetallic active site in order to stabilize the large negative charge developed on the oxygen atom of the leaving group when the P-O bond is broken in the products state. To accommodate this new ligand in the coordination shell, carbamylated Lys169 must be displaced from one zinc ion to the other, which in turn affects the acidity of Asp301, a residue originally bound to the second zinc ion. This ability to displace some of the ligands of the coordination shell of the zinc centers would explain the promiscuity of this enzyme, which is capable of catalyzing hydrolysis of different substrate by means of different mechanisms.     

A theoretical study of the mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II

The mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II has been investigated using molecular mechanics and quantum mechanics methods.Molecular dynamics(MD) simulations were carried out using the TIP3 water model and generalized solvent boundary potential(GSBP) by CHARMM based on the X-ray crystal structure.Two models of the ternary elongation complex were constructed based on CHARMM MD calculations.All the species including reactants,transition states,intermediates,and products were optimized using the DFT-PBE method coupled with the basis set DZVP and the auxiliary basis set GEN-A2.Three pathways were explored using the DFT method.The most favorable reaction pathway involves indirect proton migration from the RNA primer 3’-OH to the oxygen atom of-phosphate via a solvent water molecule,proton rotation from the oxygen atom of-phosphate to the-phosphate side,the RNA primer 3’-O nucleophilic attack on the-phosphorus atom,and P-O bond breakage.The corresponding reaction potential profile was obtained.The rate limiting step,with a barrier height of 21.5 kcal/mol,is the RNA primer 3’-O nucleophilic attack,rather than the commonly considered proton transfer process.A high-resolution crystal structure including crystallographic water molecules is required for further studies.