Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9

To elucidate enzyme catalysis through computer simulation, a prerequisite is to reliably compute free energy barriers for both enzyme and solution reactions. By employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the ab initio quantum mechanical/molecular mechanical approach and the umbrella sampling method, we have determined free energy profiles for the methyl-transfer reaction catalyzed by the histone lysine methyltransferase SET7/9 and its corresponding uncatalyzed reaction in aqueous solution, respectively. Our calculated activation free energy barrier for the enzyme catalyzed reaction is 22.5 kcal/mol, which agrees very well with the experimental value of 20.9 kcal/mol. The difference in potential of mean force between a corresponding prereaction state and the transition state for the solution reaction is computed to be 30.9 kcal/mol. Thus, our simulations indicate that the enzyme SET7/9 plays an essential catalytic role in significantly lowering the barrier for the methyl-transfer reaction step. For the reaction in solution, it is found that the hydrogen bond network near the reaction center undergoes a significant change, and there is a strong shift in electrostatic field from the prereaction state to the transition state, whereas for the enzyme reaction, such an effect is much smaller and the enzyme SET7/9 is found to provide a preorganized electrostatic environment to facilitate the methyl-transfer reaction. Meanwhile, we find that the transition state in the enzyme reaction is a little more dissociative than that in solution.     

Kinetic isotope effects of L-Dopa decarboxylase

A mixed centroid path integral and free energy perturbation method (PI-FEP/UM) has been used to investigate the primary carbon and secondary hydrogen kinetic isotope effects (KIEs) in the amino acid decarboxylation of L-Dopa catalyzed by the enzyme L-Dopa decarboxylase (DDC) along with the corresponding uncatalyzed reaction in water. DDC is a pyridoxal 5'-phosphate (PLP) dependent enzyme. The cofactor undergoes an internal proton transfer between the zwitterionic protonated Schiff base configuration and the neutral hydroxyimine tautomer. It was found that the cofactor PLP makes significant contributions to lowering the decarboxylation barrier, while the enzyme active site provides further stabilization of the transition state. Interestingly, the O-protonated configuration is preferred both in the Michaelis complex and at the decarboxylation transition state. The computed kinetic isotope effects (KIE) on the carboxylate C-13 are consistent with that observed on decarboxylation reactions of other PLP-dependent enzymes, whereas the KIEs on the α carbon and secondary proton, which can easily be validated experimentally, may be used as a possible identification for the active form of the PLP tautomer in the active site of DDC.     

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.     

吡咯喹啉醌模型化合物与氨亲核加成的理论探讨

为了揭示辅酶PQQ结构与反应性的关系,在B3LYP/D95(d, p)水平上对一系列PQQ模型化合物及其类似物与氨的亲核加成进行了理论计算.结果表明:对单羰基体系,羰基碳的亲电性对反应能垒有重要的影响;对双羰基体系,过渡态中邻位羰基氧与亲核试剂氨上的H形成的氢键对反应的活化能起着关键的作用;稠合芳香环本身对反应的能垒影响不大,但当稠合杂环的1-位为可提供氢键受体的N原子时,由于N1与氨上H原子间可形成氢键而进一步降低反应的活化能.发现过渡态中被进攻羰基与氨上N原子之间形成的夹角(OCN)与活化能有良好的线性关系.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Transannular 1,5-hydride shift in 5-hydroxycyclooctanone: an experimental and theoretical investigation

1H NMR measurements have shown that eight out of twelve methylene hydrogen atoms of 5-hydroxycyclooctanone (1HK) are exchanged for deuterium atoms under acidic and basic conditions. For the reaction in 7.9 M DCl/D2O, an activation energy Ea = 19.2 +/- 0.4 kcal mol(-1) is found. In order to explain these findings, a degenerate transannular 1,5-hydride shift is essential, and this reaction has been analyzed by quantum chemical calculations. Hydride transfer takes place via a tight transition state with a six-membered ring. The activation barrier is lowest in the presence of base. The molecular conformation of the eight-membered ring in the transition state resembles that of the starting structure. It is unlikely that solvent molecules such as water participate in the formation of the transition state.     

Hydride transfer in liver alcohol dehydrogenase: quantum dynamics, kinetic isotope effects, and role of enzyme motion.

The quantum dynamics of the hydride transfer reaction catalyzed by liver alcohol dehydrogenase (LADH) are studied with real-time dynamical simulations including the motion of the entire solvated enzyme. The electronic quantum effects are incorporated with an empirical valence bond potential, and the nuclear quantum effects of the transferring hydrogen are incorporated with a mixed quantum/classical molecular dynamics method in which the transferring hydrogen nucleus is represented by a three-dimensional vibrational wave function. The equilibrium transition state theory rate constants are determined from the adiabatic quantum free energy profiles, which include the free energy of the zero point motion for the transferring nucleus. The nonequilibrium dynamical effects are determined by calculating the transmission coefficients with a reactive flux scheme based on real-time molecular dynamics with quantum transitions (MDQT) surface hopping trajectories. The values of nearly unity for these transmission coefficients imply that nonequilibrium dynamical effects such as barrier recrossings are not dominant for this reaction. The calculated deuterium and tritium kinetic isotope effects for the overall rate agree with experimental results. These simulations elucidate the fundamental nature of the nuclear quantum effects and provide evidence of hydrogen tunneling in the direction along the donor-acceptor axis. An analysis of the geometrical parameters during the equilibrium and nonequilibrium simulations provides insight into the relation between specific enzyme motions and enzyme activity. The donor-acceptor distance, the catalytic zinc-substrate oxygen distance, and the coenzyme (NAD(+)/NADH) ring angles are found to strongly impact the activation free energy barrier, while the donor-acceptor distance and one of the coenzyme ring angles are found to be correlated to the degree of barrier recrossing. The distance between VAL-203 and the reactive center is found to significantly impact the activation free energy but not the degree of barrier recrossing. This result indicates that the experimentally observed effect of mutating VAL-203 on the enzyme activity is due to the alteration of the equilibrium free energy difference between the transition state and the reactant rather than nonequilibrium dynamical factors. The promoting motion of VAL-203 is characterized in terms of steric interactions involving THR-178 and the coenzyme.     

On the mechanisms of oxidation of organic sulfides by H2O2 in aqueous solutions

The mechanism of oxidation of organic sulfides in aqueous solutions by hydrogen peroxide was investigated via ab initio calculations. Specifically, two reactions, hydrogen transfer of hydrogen peroxide to form water oxide and the oxidation of dimethyl sulfide (DMS) by hydrogen peroxide to form dimethyl sulfoxide, were studied as models of these processes in general. Solvent effects are included both via including explicitly water molecules and via the polarizable continuum model. The former was found to have a much more significant effect than the latter. When explicit water molecules are included, a mechanism different from those proposed in the literature was found. Specific interactions including hydrogen bonding with 2-3 water molecules can provide enough stabilization for the charge separation of the activation complex. The energy barrier of the oxidation of DMS by hydrogen peroxide was estimated to be 12.7 kcal/mol, within the experimental range of the oxidation of analogous compounds (10-20 kcal/mol). The major reaction coordinates of the reaction are the breaking of the O-O bond of H2O2 and the formation of the S-O bond, the transfer of hydrogen to the distal oxygen of hydrogen peroxide occurring after the system has passed the transition state. Reaction barriers of the hydrogen transfer of H2O2 are an average of 10 kcal/mol or higher than the reaction barriers of the oxidation of DMS. Therefore, a two-step oxidation mechanism in which, first, the transfer of a hydrogen atom occurs to form water oxide and, second, the transfer of oxygen to the substrate occurs is unlikely to be correct. Our proposed oxidation mechanism does not suggest a pH dependence of oxidation rate within a moderate range around neutral pH (i.e., under conditions in which hydronium and hydroxide ions do not participate directly in the reaction), and it agrees with experimental observations over moderate pH values. Also, without including a protonated solvent molecule, it has activation energies that correspond to measured activation energies.     

Water catalysis and anticatalysis in photochemical reactions: observation of a delayed threshold effect in the reaction quantum yield

The possible catalysis of photochemical reactions by water molecules is considered. Using theoretical simulations, we investigate the HF-elimination reaction of fluoromethanol in small water clusters initiated by the overtone excitation of the hydroxyl group. The reaction occurs in competition with the process of water evaporation that dissipates the excitation and quenches the reaction. Although the transition state barrier is stabilized by over 20 kcal/mol through hydrogen bonding with water, the quantum yield versus energy shows a pronounced delayed threshold that effectively eliminates the catalytic effect. It is concluded that the quantum chemistry calculations of barrier lowering are not sufficient to infer water catalysis in some photochemical reactions, which instead require dynamical modeling.     

(H_2O)_n(n=0～2)催化HS和HOCl反应机理的理论研究

在CCSD(T)/6-311+G(3df,2p)//M06-2X/6-311+G(3df,2p)水平上研究了(H_2O)n(n=0~2)催化HS和HOCl的反应机理.结果表明,HS与HOCl反应中HS夺取HOCl上的H原子形成产物H_2S和ClO.在无水催化时,该反应存在2种不同的路径(分别经过过渡态TS1和TS2,二者互为顺反结构),对应的能垒分别为100.28和100.91kJ/mol,到达产物(H_2S+ClO)需吸收18.99kJ/mol能量,反应不易发生;在单个水分子参与时,水分子可通过形成弱相互作用或者作为H原子转移桥梁影响反应机理,获得了4种水催化路径,能垒(间于53.97~92.39kJ/mol之间)均低于无水催化过程.同时发现,在反应到达产物前,水分子可以与产物形成中间体IM,IM相对能仅为0.46kJ/mol,有利于产物形成;有2个水分子参与反应时,找到了3条催化路径,最优反应路径过渡态TS7的能垒为45.05kJ/mol,低于无水催化过程,相比单个水分子最优路径能垒(53.97kJ/mol)并无显著降低.     

Reaction force analyses of nitro-aci tautomerizations of trinitromethane, the elusive trinitromethanol, picric acid and 2,4-dinitro-1H-imidazole

We have analyzed computationally, in terms of the reaction force, the nitro → aci tautomerizations of trinitromethane, trinitromethanol, picric acid and 2,4-dinitro-1H-imidazole. These processes involve intramolecular transfer of a hydrogen to an NO₂ oxygen, forming the aci tautomer (a nitronic acid). The reaction force naturally and unambiguously divides an activation barrier into two components: (1) the energy required for initial structural changes in the reactant(s), and (2) the energy associated with the first portion of the transition to product(s). In each of these tautomerizations, the first component is dominant. For trinitromethane, it is so large that the resulting total activation barrier makes C–NO₂ bond scission energetically preferable. On the other hand, trinitromethanol—which appears to be unknown—readily undergoes fragmentation in conjunction with hydrogen transfer. Picric acid has the interesting feature that the reaction is almost complete after the first portion of the activation process, marked by the minimum of the reaction force. In all four reactions, the properties of the systems at the force minimum, transition state and force maximum are consistent with the concept of a “transition” region in a chemical reaction versus simply a transition state.     

Characterizing the mechanism of the double proton transfer in the formamide dimer

The double proton transfer in the formamide dimer is characterized computationally by combining density functional theory and ab initio methods. The intrinsic reaction coordinate (IRC) is obtained at the B3LYP level of theory. Energies of several points along the IRC are treated by the more rigorous focal point method to test the validity of the B3LYP functional. The reaction mechanism is examined in terms of the energy profile, the reaction force, the chemical potential, and the reaction electronic flux. The energy profile for the activation process of the formamide dimer to the imino ether product obtained with the B3LYP functional is in agreement with the results of the focal point method. Together with the reaction force analysis and the reaction electronic flux a precise assignment of the structural and electronic contributions to the activation barrier becomes possible. The results show that the reaction starts with a structural rearrangement, where the two dimers approach each other, and is followed by electronic changes before the system reaches the transition state. This electronic contribution to the activation barrier steers the activation process. After the transition state is reached, deviations of the B3LYP functional from the more accurate focal point energies become apparent, where the errors may be rationalized in terms of the treatment of exchange. The inconsistency could be assigned to the incapacity of the functional to describe delocalization effects over the whole system.     

Theoretical study of the gas-phase reaction of hydrogen isocyanate with water

Transition states for the two probable pathways for the gas-phase hydrolysis of hydrogen isocyanate have been determined using the MINDO /3 method. Activation barriers obtained showed that the one-step mechanism is preferred to the two-step mechanism involving an intermediate. It was shown that the reaction of polymeric water has lower activation barrier than the reaction of monomeric water. Energy and charge decomposition analysis showed that in the former less energy is required in the deformation of molecules for the transition state (TS ) formation due to the cyclic flow of electronic charges in the TS .     

甲醇和氨对腺嘌呤水解脱氨机理的影响(英文)

用密度泛函理论在B3LYP/6-311G(d,p)水平上研究了甲醇和氨辅助的腺嘌呤水解脱氨机理,结果表明该反应首先是在水的亲核进攻下形成一个四面体结构中间体,然后该中间体通过构象变化得到两个不同构象,从而找到两条不同反应路径.在路径a中,辅助分子参与了过渡态的形成,起到转移氢原子的作用;而在路径b中,辅助分子仅作为介质,没有参与过渡态的形成.在氨辅助的情况下,腺嘌呤在亲核反应前发生了胺-亚胺异构化,而在甲醇辅助机理中则未发生该异构化.能量结果表明甲醇辅助腺嘌呤脱氨反应具有与水辅助类似的势垒,而氨辅助反应的势垒则比水辅助的高.     

Effect of electrostatic interaction on the mechanism of dehalogenation catalyzed by haloalkane dehalogenase

Theoretical calculation has been carried out for the nucleophilic displacement reaction of 1,2‐dichloroethane catalyzed by haloalkane dehalogenase. The results indicate that different hydrogen bond patterns of the oxyanion hole and the halide‐stabilizing residues play an important role in the dehalogenation reaction. They cause concertedly an earlier transition state (TS) with the activation barrier of 16.60 kcal/mol. The stabilization effect of Trp125 and Trp175 on chlorine atom in the TS is larger than that of the reactant complex by 15.67 kcal/mol so that, they make contribution to the stabilization of the TS. Moreover, the reaction shows the enzymatic action can be attributed to a combination of reactant‐state destabilization and transition‐state electrostatic stabilization. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.     

A first-principles study of molecular oxygen dissociation at an electrode surface: a comparison of potential variation and coadsorption effects

Influences of coadsorbed sodium and water, aqueous solvent, and electrode potential on the kinetics of O(2) dissociation over Pt(111) are systematically investigated using density functional theory models of vacuum and electrochemical interfaces. Na coadsorption alters the electronic states of Pt to stabilize the reactant (O(2)*), transition, and product (2O*) states by facilitating electron donation to oxygen, causing a more exothermic reaction energy (-0.84 eV for Na and O(2), -0.81 eV for isolated O(2)) and a decrease in dissociation barrier (0.39 eV for Na and O(2), 0.57 eV for isolated O(2)). Solvation decreases the reaction energy (-0.67 eV) due to enhanced hydrogen bond stabilization of O(2)* compared to 2O*. The influence of Na is less pronounced at the solvated interface (barrier decreases by only 0.11 eV) because H(2)O screens Na charge-donation. In the electrochemical model system, the dissociation energy becomes more exothermic and the barrier decreases toward more positive potentials. Potential-dependent behavior results from changes in interfacial dipole moment and polarizability between O(2)*, the dissociation transition state, and 2O*; each are influenced by changes in adsorption and hydrogen bonding. Coadsorption of Na in the solvated system dampens the dipole moment change between O(2)* and 2O* and significantly increases the polarizability at the dissociation transition state and for 2O*; the combination causes little change in the reaction energy but reduces the activation barrier by 0.08 eV at 0 V versus NHE. The potential-dependent behavior contrasts that determined at a constant surface charge or from an applied electric field, illustrating the importance of considering the electrochemical potential at the fully-solvated interface in determining reaction energetics, even for non-redox reactions.     

Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions

Chorismate mutase is at the centre of current controversy about fundamental features of biological catalysts. Some recent studies have proposed that catalysis in this enzyme does not involve transition state (TS) stabilization but instead is due largely to the formation of a reactive conformation of the substrate. To understand the origins of catalysis, it is necessary to compare equivalent reactions in different environments. The pericyclic conversion of chorismate to prephenate catalysed by chorismate mutase also occurs (much more slowly) in aqueous solution. In this study we analyse the origins of catalysis by comparison of multiple quantum mechanics/molecular mechanics (QM/MM) reaction pathways at a reliable, well tested level of theory (B3LYP/6-31G(d)/CHARMM27) for the reaction (i) in Bacillus subtilis chorismate mutase (BsCM) and (ii) in aqueous solvent. The average calculated reaction (potential energy) barriers are 11.3 kcal mol(-1) in the enzyme and 17.4 kcal mol(-1) in water, both of which are in good agreement with experiment. Comparison of the two sets of reaction pathways shows that the reaction follows a slightly different reaction pathway in the enzyme than in it does in solution, because of a destabilization, or strain, of the substrate in the enzyme. The substrate strain energy within the enzyme remains constant throughout the reaction. There is no unique reactive conformation of the substrate common to both environments, and the transition state structures are also different in the enzyme and in water. Analysis of the barrier heights in each environment shows a clear correlation between TS stabilization and the barrier height. The average differential TS stabilization is 7.3 kcal mol(-1) in the enzyme. This is significantly higher than the small amount of TS stabilization in water (on average only 1.0 kcal mol(-1) relative to the substrate). The TS is stabilized mainly by electrostatic interactions with active site residues in the enzyme, with Arg90, Arg7 and Glu78 generally the most important. Conformational effects (e.g. strain of the substrate in the enzyme) do not contribute significantly to the lower barrier observed in the enzyme. The results show that catalysis is mainly due to better TS stabilization by the enzyme.     

Theoretical studies on the hydrolysis of phosphate diesters in the gas phase,solution, and RNase A

Density functional theory, polarizable continuum models and semiempirical hybrid quantum mechanical/molecular mechanical (QM/MM) calculations were applied to the hydrolysis of phosphate diesters in the gas phase, in solution, and in the enzyme RNase A. Neutralization of the negative charge of the pentacovalent phosphorane intermediates provides a substantial stabilization of the transition‐state structures in the gas phase. Inclusion of solvent effects on the phosphate/phosphorane species was critical to reproducing the trends in reactivity observed experimentally. Finally, the catalytic mechanism for the hydrolysis of uridine 2′,3′‐cyclic phosphate by RNase A was studied by QM/MM calculations. Our results suggest that the rate‐limiting transition state of the reaction corresponds to the approach of a water molecule to the phosphate and its activation by His119. Thus, His119 acts as a generalized base for the reaction. The water attack leads to a pentacovalent phosphorane transition state of formal charge ?2; this excess of negative charge in the transition state is stabilized by a number of positively charged residues including His12 and Lys41. In the second stage of the reaction, the phosphorane is converted into products. This part of the reaction proceeds without a detectable barrier, and it is facilitated by a proton transfer from Lys41 to the departing O_2′. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001