Investigation of the enzymatic and nonenzymatic cope rearrangement of carbaprephenate to carbachorismate

The dimethyl esters of carbaprephenate and 4-epi-carbaprephenate were prepared by modification of published procedures. In methanol these compounds are converted quantitatively to isomeric 6-hydroxytricyclo[3.3.1.0(2,7)]non-3-en-1,3-dimethyl esters via a two-step sequence involving an initial Cope rearrangement, followed by intramolecular Diels-Alder reaction of the dimethyl carbachorismate or 4-epi-carbachorismate intermediates. Carbaprephenate and its epimer were obtained by alkaline hydrolysis of the corresponding dimethyl esters. These compounds, in contrast to their ester precursors, undergo spontaneous acid-catalyzed decarboxylation in aqueous solution. Only at high pH does the Cope rearrangement compete with decarboxylation. At pH 12 and 90 degrees C, carbaprephenate slowly rearranges to carbachorismate, which rapidly loses water to give 3-(2-carboxyallyl)benzoic acid as the major product. A small amount of the intramolecular Diels-Alder adduct derived from carbachorismate is also observed by NMR as a minor product. Carbaprephenate is not a substrate for the enzyme chorismate mutase from Bacillus subtilis (BsCM), nor does carbaprephenate inhibit the normal chorismate mutase activity of this enzyme, even when present in 200-fold excess over chorismate. Its low affinity for the enzyme-active site is presumably a consequence of placing a methylene group rather than an oxygen atom proximal to the essential cationic residue Arg90. Nevertheless, BsCM variants that lack this cation (R90G and R90A) do not accelerate the Cope rearrangement of carbaprephenate either, and a catalytic antibody 1F7, which exhibits modest chorismate mutase activity, is similarly inactive. Poor substrate binding and the relatively high barrier for the Cope compared to the Claisen rearrangement presumably account for the lack of detectable catalysis. Acceleration of this sigmatropic rearrangement apparently requires more than an active site that is complementary in shape to the reactive substrate conformer.     

A comparative study of claisen and cope rearrangements catalyzed by chorismate mutase. An insight into enzymatic efficiency: transition state stabilization or substrate preorganization?

In this work we present a detailed analysis of the activation free energies and averaged interactions for the Claisen and Cope rearrangements of chorismate and carbachorismate catalyzed by Bacillus subtilischorismate mutase (BsCM) using quantum mechanics/molecular mechanics (QM/MM) simulation methods. In gas phase, both reactions are described as concerted processes, with the activation free energy for carbachorismate being about 10-15 kcal mol(-)(1) larger than for chorismate, at the AM1 and B3LYP/6-31G levels. Aqueous solution and BsCM active site environments reduce the free energy barriers for both reactions, due to the fact that in these media the two carboxylate groups can be approached more easily than in the gas phase. The enzyme specifically reduces the activation free energy of the Claisen rearrangement about 3 kcal mol(-)(1) more than that for the Cope reaction. This result is due to a larger transition state stabilization associated to the formation of a hydrogen bond between Arg90 and the ether oxygen. When this oxygen atom is changed by a methylene group, the interaction is lost and Arg90 moves inside the active site establishing stronger interactions with one of the carboxylate groups. This fact yields a more intense rearrangement of the substrate structure. Comparing two reactions in the same enzyme, we have been able to obtain conclusions about the relative magnitude of the substrate preorganization and transition state stabilization effects. Transition state stabilization seems to be the dominant effect in this case.     

Differential transition-state stabilization in enzyme catalysis: quantum chemical analysis of interactions in the chorismate mutase reaction and prediction of the optimal catalytic field

Chorismate mutase is a key model system in the development of theories of enzyme catalysis. To analyze the physical nature of catalytic interactions within the enzyme active site and to estimate the stabilization of the transition state (TS) relative to the substrate (differential transition state stabilization, DTSS), we have carried out nonempirical variation-perturbation analysis of the electrostatic, exchange, delocalization, and correlation interactions of the enzyme-bound substrate and transition-state structures derived from ab initio QM/MM modeling of Bacillus subtilis chorismate mutase. Significant TS stabilization by approximately -23 kcal/mol [MP2/6-31G(d)] relative to the bound substrate is in agreement with that of previous QM/MM modeling and contrasts with suggestions that catalysis by this enzyme arises purely from conformational selection effects. The most important contributions to DTSS come from the residues, Arg90, Arg7, Glu78, a crystallographic water molecule, Arg116, and Arg63, and are dominated by electrostatic effects. Analysis of the differential electrostatic potential of the TS and substrate allows calculation of the catalytic field, predicting the optimal location of charged groups to achieve maximal DTSS. Comparison with the active site of the enzyme from those of several species shows that the positions of charged active site residues correspond closely to the optimal catalytic field, showing that the enzyme has evolved specifically to stabilize the TS relative to the substrate.     

Preorganization and reorganization as related factors in enzyme catalysis: the chorismate mutase case

In this paper a deeper insight into the chorismate-to prephenate-rearrangement, catalyzed by Bacillus subtilis chorismate mutase, is provided by means of a combination of statistical quantum mechanics/molecular mechanics simulation methods and hybrid potential energy surface exploration techniques. The main aim of this work is to present an estimation of the preorganization and reorganization terms of the enzyme catalytic rate enhancement. To analyze the first of these, we have studied different conformational equilibria of chorismate in aqueous solution and in the enzyme active site. Our conclusion is that chorismate mutase preferentially binds the reactive conformer of the substrate--that presenting a structure similar to the transition state of the reaction to be catalyzed--with shorter distances between the carbon atoms to be bonded and more diaxial character. With respect to the reorganization effect, an energy decomposition analysis of the potential energies of the reactive reactant and of the reaction transition state in aqueous solution and in the enzyme shows that the enzyme structure is better adapted to the transition structure. This means not only a more negative electrostatic interaction energy with the transition state but also a low enzyme deformation contribution to the energy barrier. Our calculations reveal that the structure of the enzyme is responsible for stabilizing the transition state structure of the reaction, with concomitant selection of the reactive form of the reactants. This is, the same enzymatic pattern that stabilizes the transition structure also promotes those reactant structures closer to the transition structure (i.e., the reactive reactants). In fact, both reorganization and preorganization effects have to be considered as the two faces of the same coin, having a common origin in the effect of the enzyme structure on the energy surface of the substrate.     

Isotope effects on the enzymatic and nonenzymatic reactions of chorismate

The important biosynthetic intermediate chorismate reacts thermally by two competitive pathways, one leading to 4-hydroxybenzoate via elimination of the enolpyruvyl side chain, and the other to prephenate by a facile Claisen rearrangement. Measurements with isotopically labeled chorismate derivatives indicate that both are concerted sigmatropic processes, controlled by the orientation of the enolpyruvyl group. In the elimination reaction of [4-2H]chorismate, roughly 60% of the label was found in pyruvate after 3 h at 60 degrees C. Moreover, a 1.846 +/- 0.057 2H isotope effect for the transferred hydrogen atom and a 1.0374 +/- 0.0005 18O isotope effect for the ether oxygen show that the transition state for this process is highly asymmetric, with hydrogen atom transfer from C4 to C9 significantly less advanced than C-O bond cleavage. In the competing Claisen rearrangement, a very large 18O isotope effect at the bond-breaking position (1.0482 +/- 0.0005) and a smaller 13C isotope effect at the bond-making position (1.0118 +/- 0.0004) were determined. Isotope effects of similar magnitude characterized the transformations catalyzed by evolutionarily unrelated chorismate mutases from Escherichia coli and Bacillus subtilis. The enzymatic reactions, like their solution counterpart, are thus concerted [3,3]-sigmatropic processes in which C-C bond formation lags behind C-O bond cleavage. However, as substantially larger 18O and smaller 13C isotope effects were observed for a mutant enzyme in which chemistry is fully rate determining, the ionic active site may favor a somewhat more polarized transition state than that seen in solution.     

A hybrid potential reaction path and free energy study of the chorismate mutase reaction.

We present a combination of two techniques--QM/MM statistical simulation methods and QM/MM internal energy minimizations--to get a deeper insight into the reaction catalyzed by the enzyme chorismate mutase. Structures, internal energies and free energies, taken from the paths of the reaction in solution and in the enzyme have been analyzed in order to estimate the relative importance of the reorganization and preorganization effects. The results we obtain for this reaction are in good agreement with experiment and show that chorismate mutase achieves its catalytic efficiency in two ways; first, it preferentially binds the active conformer of the substrate and, second, it reduces the free energy of activation for the reaction relative to that in solution by providing an environment which stabilizes the transition state.     

Apparent NAC effect in chorismate mutase reflects electrostatic transition state stabilization

The catalytic reaction of chorismate mutase (CM) has been the subject of major current attention. Nevertheless, the origin of the catalytic power of CM remains an open question. In particular, it has not been clear whether the enzyme works by providing electrostatic transition state stabilization (TSS), by applying steric strain, or by populating near attack conformation (NAC). The present work explores this issue by a systematic quantitative analysis. The overall catalytic effect is reproduced by the empirical valence bond (EVB) method. In addition, the binding free energy of the ground state and the transition state is evaluated, demonstrating that the enzyme works by TSS. Furthermore, the evaluation of the electrostatic contribution to the reduction of the activation energy establishes that the TSS results from electrostatic effects. It is also found that the apparent NAC effect is not the reason for the catalytic effect but the result of the TSS. It is concluded that in CM as in other enzymes the key catalytic effect is electrostatic TSS. However, since the charge distribution of the transition state and the reactant state is similar, the stabilization of the transition state leads to reduction in the distance between the reacting atoms in the reactant state.     

Comparison of formation of reactive conformers (NACs) for the Claisen rearrangement of chorismate to prephenate in water and in the E. coli mutase: the efficiency of the enzyme catalysis

The Claisen rearrangements of chorismate (CHOR) in water and at the active site of E. coli chorismate mutase (EcCM) have been compared. From a total of 33 ns molecular dynamics simulation of chorismate in water solvent, seven diaxial conformers I-VII were identified. Most of the time (approximately 99%), the side chain carboxylate of the chorismate is positioned away from the ring due to the electrostatic repulsion from the carboxylate in the ring. Proximity of the two carboxylates, as seen in conformer I, is a requirement for the formation of a near attack conformer (NAC) that can proceed to the transition state (TS). In the EcCM.CHOR complex, the two carboxylates of CHOR are tightly held by Arg28 of one subunit and Arg11* of the other subunit, resulting in the side chain C16 being positioned adjacent to C5 with their motions restricted by van der Waals contacts with methyl groups of Val35 and Ile81. With the definition of NAC as the C5...C16 distance < or =3.7 A and the attack angle < or =30 degrees, it was estimated from our MD trajectories that the free energy of NAC formation is approximately 8.4 kcal/mol above the total ground state in water, whereas in the enzyme it is only 0.6 kcal/mol above the average of the Michaelis complex EcCM.CHOR. The experimentally measured difference in the activation free energies of the water and enzymatic reactions (Delta Delta G(++)) is 9 kcal/mol. It follows that the efficiency of formation of NAC (7.8 kcal/mol) at the active site provides approximately 90% of the kinetic advantage of the enzymatic reaction as compared to the water reaction. Comparison of the EcCM.TSA (transition state analogue) and EcCM.NAC simulations suggests that the experimentally measured 100 fold tighter binding of TSA compared to CHOR does not originate from the difference between NAC and the TS binding affinities, but might be due to the free energy cost to bring the two carboxylates of CHOR together to interact with Arg28 and Arg11* at the active site. The two carboxylates of TSA are fixed by a bicyclic structure. The remaining approximately 10% of Delta Delta G(++) may be attributed to a preferential interaction of Lys39-NH(3)(+) with O13 ether oxygen in the TS.     

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.     

Just a near attack conformer for catalysis (chorismate to prephenate rearrangements in water,antibody, enzymes,and their mutants)

Standard free energies for formation of ground-state reactive conformers (DeltaGN degrees ) and transition states (DeltaG) in the conversion of chorismate to prephenate in water, B. subtilis mutase, E. coli mutase, and their mutants, as well as a catalytic antibody, are related by DeltaG = DeltaGN degrees + 16 kcal/mol. Thus, the differences in the rate constants for the water reaction and catalysts reactions reside in the mole fraction of substrate present as reactive conformers (NACs). These results, and knowledge of the importance of transition state stabilization in other cases, suggest a proposal that enzymes utilize both NAC and transition state stabilization in the mix required for the most efficient catalysis.     

Entropic and enthalpic components of catalysis in the mutase and lyase activities of Pseudomonas aeruginosa PchB

The isochorismate-pyruvate lyase from Pseudomonas aeruginosa (PchB) catalyzes two pericyclic reactions, demonstrating the eponymous activity and also chorismate mutase activity. The thermodynamic parameters for these enzyme-catalyzed activities, as well as the uncatalyzed isochorismate decomposition, are reported from temperature dependence of k(cat) and k(uncat) data. The entropic effects do not contribute to enzyme catalysis as expected from previously reported chorismate mutase data. Indeed, an entropic penalty for the enzyme-catalyzed mutase reaction (ΔS(++) = -12.1 ± 0.6 cal/(mol K)) is comparable to that of the previously reported uncatalyzed reaction, whereas that of the enzyme-catalyzed lyase reaction (ΔS(++) = -24.3 ± 0.2 cal/(mol K)) is larger than that of the uncatalyzed lyase reaction (-15.77 ± 0.02 cal/(mol K)) documented here. With the assumption that chemistry is rate-limiting, we propose that a reactive substrate conformation is formed upon loop closure of the active site and that ordering of the loop contributes to the entropic penalty for converting the enzyme substrate complex to the transition state.     

Charge optimization increases the potency and selectivity of a chorismate mutase inhibitor

The highest affinity inhibitor for chorismate mutases, a conformationally constrained oxabicyclic dicarboxylate transition state analogue, was modified as suggested by computational charge optimization methods. As predicted, replacement of the C10 carboxylate in this molecule with a nitro group yields an even more potent inhibitor of a chorismate mutase from Bacillus subtilis (BsCM), but the magnitude of the improvement (roughly 3-fold, corresponding to a DeltaDeltaG of -0.7 kcal/mol) is substantially lower than the gain of 2-3 kcal/mol binding free energy anticipated for the reduced desolvation penalty upon binding. Experiments with a truncated version of the enzyme show that the flexible C terminus, which was only partially resolved in the crystal structure and hence omitted from the calculations, provides favorable interactions with the C10 group that partially compensate for its desolvation. Although truncation diminishes the affinity of the enzyme for both inhibitors, the nitro derivative binds 1.7 kcal/mol more tightly than the dicarboxylate, in reasonable agreement with the calculations. Significantly, substitution of the C10 carboxylate with a nitro group also enhances the selectivity of inhibition of BsCM relative to a chorismate mutase from Escherichia coli (EcCM), which has a completely different fold and binding pocket, by 10-fold. These results experimentally verify the utility of charge optimization methods for improving interactions between proteins and low-molecular weight ligands.     

Contributions of conformational compression and preferential transition state stabilization to the rate enhancement by chorismate mutase

The rate enhancement provided by the chorismate mutase (CM) enzyme for the Claisen rearrangement of chorismate to prephenate has been investigated by application of the concept of near attack conformations (NACs). Using a combined QM/MM Monte Carlo/free-energy perturbation (MC/FEP) method, 82% and 100% of chorismate conformers were found to be NAC structures in water and in the CM active site, respectively. Consequently, the conversion of non-NACs to NACs does not contribute to the free energy of activation from preorganization of the substrate into NACs. The FEP calculations yielded differences in free energies of activation that well reproduce the experimental data. Additional calculations indicate that the rate enhancement by CM over the aqueous phase results primarily from conformational compression of NACs by the enzyme and that this process is enthalpically controlled. This suggests that preferential stabilization of the transition state in the enzyme environment relative to water plays a secondary role in the catalysis by CM.     

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.     

Quantitative evaluation of noncovalent chorismate mutase-inhibitor binding by ESI-MS

Electrospray time-of-flight mass spectrometry was used to quantitatively determine the dissociation constant of chorismate mutase and a transition state analogue inhibitor. This system presents a fairly complex stoichiometry because the native protein is a homotrimer with three equal and independent substrate binding sites. We can detect the chorismate mutase trimer as well as chorismate mutase-inhibitor complexes by choosing appropriate conditions in the ESI source. To verify that the protein-inhibitor complexes are specific, titration experiments with different enzyme variants and different inhibitors were performed. A plot of the number of bound inhibitors versus added inhibitor concentration revealed saturation behavior with 3:1 (inhibitor:functional trimer) stoichiometry for the TSA. The soft ESI conditions, the relatively high protein mass of 43.5 kDa, and the low charge state (high m/z) result in broad peaks, a typical problem in analyzing noncovalent protein complexes. Due to the low molecular weight of the TSA (226 Da) the peaks of the free protein and the protein with one, two or three inhibitors bound cannot be clearly resolved. For data analysis, relative peak areas of the deconvoluted spectra of chorismate mutase-inhibitor complexes were obtained by fitting appropriate peak shapes to the signals corresponding to the free enzyme and its complexes with one, two, or three inhibitor molecules. From the relative peak areas we were able to calculate a dissociation constant that agreed well with known solution-phase data. This method may be generally useful for interpreting mass spectra of noncovalent complexes that exhibit broad peaks in the high m/z range.     

The catalytic mechanism of an aspartic proteinase explored with neutron and X-ray diffraction

Hydrogen atoms play key roles in enzyme mechanism, but as this study shows, even high-quality X-ray data to a resolution of 1 A cannot directly visualize them. Neutron diffraction, however, can locate deuterium atoms even at resolutions around 2 A. Both neutron and X-ray diffraction data have been used to investigate the transition state of the aspartic proteinase endothiapepsin. The different techniques reveal a different part of the story, revealing the clearest picture yet of the catalytic mechanism by which the enzyme operates. Room temperature neutron and X-ray diffraction data were used in a newly developed joint refinement software package to visualize deuterium atoms within the active site of the enzyme when a gem-diol transition state analogue inhibitor is bound at the active site. These data were also used to estimate their individual occupancy, while analysis of the differences between the bond lengths of the catalytic aspartates was performed using atomic resolution X-ray data. The two methods are in agreement on the protonation state of the active site with a transition state analogue inhibitor bound confirming the catalytic mechanism at which the enzyme operates.     

Probing the origin of the compromised catalysis of E. coli alkaline phosphatase in its promiscuous sulfatase reaction

The catalytic promiscuity of E. coli alkaline phosphatase (AP) and many other enzymes provides a unique opportunity to dissect the origin of enzymatic rate enhancements via a comparative approach. Here, we use kinetic isotope effects (KIEs) to explore the origin of the 109-fold greater catalytic proficiency by AP for phosphate monoester hydrolysis relative to sulfate monoester hydrolysis. The primary 18O KIEs for the leaving group oxygen atoms in the AP-catalyzed hydrolysis of p-nitrophenyl phosphate (pNPP) and p-nitrophenylsulfate (pNPS) decrease relative to the values observed for nonenzymatic hydrolysis reactions. Prior linear free energy relationship results suggest that the transition states for AP-catalyzed reactions of phosphate and sulfate esters are "loose" and indistinguishable from that in solution, suggesting that the decreased primary KIEs do not reflect a change in the nature of the transition state but rather a strong interaction of the leaving group oxygen atom with an active site Zn2+ ion. Furthermore, the primary KIEs for the two reactions are identical within error, suggesting that the differential catalysis of these reactions cannot be attributed to differential stabilization of the leaving group. In contrast, AP perturbs the KIE for the nonbridging oxygen atoms in the reaction of pNPP but not pNPS, suggesting a differential interaction with the transferred group in the transition state. These and prior results are consistent with a strong electrostatic interaction between the active site bimetallo Zn2+ cluster and one of the nonbridging oxygen atoms on the transferred group. We suggest that the lower charge density of this oxygen atom on a transferred sulfuryl group accounts for a large fraction of the decreased stabilization of the transition state for its reaction relative to phosphoryl transfer.     

Pentavalent Organo-Vanadates as Transition State Analogues for Phosphoryl Transfer Reactions

Pentavalent organo-vanadates have been put forth as transition state analogues for a variety of phosphoryl transfer reactions. In particular, uridine 2',3'-cyclic vanadate (U>v) has been proposed to resemble the transition state during catalysis by ribonuclease A (RNase A). Here, this hypothesis is tested. Lys41 of RNase A is known to donate a hydrogen bond to a nonbridging phosphoryl oxygen in the transition state during catalysis. Site-directed mutagenesis and semisynthesis were used to create enzymes with natural and nonnatural amino acid residues at position 41. These variants differ by 10(5)-fold in their k(cat)/K(m) values for catalysis, but <40-fold in their K(i) values for inhibition of catalysis by U>v. Plots of logK(i) vs log(K(m)/k(cat)) for three distinct substrates [poly(cytidylic acid), uridine 3'-(p-nitrophenyl phosphate), and cytidine 2',3'-cyclic phosphate] have slopes that range from 0.25 and 0.36. These plots would have a slope of unity if U>v were a perfect transition state analogue. Values of K(i) for U>v correlate weakly with the equilibrium dissociation constant for the enzymic complexes with substrate or product, indicating that U>v bears some resemblance to the substrate and product as well as the transition state. Thus, U>v is a transition state analogue for RNase A, but only a marginal one. This finding indicates that a pentavalent organo-vanadate cannot necessarily be the basis for a rigorous analysis of the transition state for a phosphoryl transfer reaction.     

Transition structure selectivity in enzyme catalysis: a QM/MM study of chorismate mutase

Two different transition structures (TSs) have been located and characterized for the chorismate conversion to prephenate in Bacillus subtilis chorismate mutase by means of hybrid quantum-mechanical/molecular-mechanical (QM/MM) calculations. GRACE software, combined with an AM1/CHARMM24/TIP3P potential, has been used involving full gradient relaxation of the position of ca. 3300 atoms. These TSs have been connected with their respective reactants and products by the intrinsic reaction coordinate (IRC) procedure carried out in the presence of the protein environment, thus obtaining for the first time a realistic enzymatic reaction path for this reaction. Similar QM/MM computational schemes have been applied to study the chemical reaction solvated by ca. 500 water molecules. Comparison of these results together with gas phase calculations has allowed understanding of the catalytic efficiency of the protein. The enzyme stabilizes one of the TSs (TS_OHout) by means of specific hydrogen bond interactions, while the other TS (TS_OHin) is the preferred one in vacuum and in water. The enzyme TS is effectively more polarized but less dissociative than the corresponding solvent and gas phase TSs. Electrostatic stabilization and an intramolecular charge-transfer process can explain this enzymatically induced change. Our theoretical results provide new information on an important enzymatic transformation and the key factors responsible for efficient selectivity are clarified. Received: 25 March 2000 / Accepted: 7 August 2000 / Published online: 23 November 2000     

The role of the cleavage site 2'-hydroxyl in the Tetrahymena group I ribozyme reaction

BACKGROUND: The 2'-hydroxyl of U preceding the cleavage site, U(-1), in the Tetrahymena ribozyme reaction contributes 10(3)-fold to catalysis relative to a 2'-hydrogen atom. Previously proposed models for the catalytic role of this 2'-OH involve coordination of a catalytic metal ion and hydrogen-bond donation to the 3'-bridging oxygen. An additional model, hydrogen-bond donation by the 2'-OH to a nonbridging reactive phosphoryl oxygen, is also consistent with previous results. We have tested these models using atomic-level substrate modifications and kinetic and thermodynamic analyses. RESULTS: Replacing the 2'-OH with -NH(3)(+) increases the reaction rate approximately 60-fold, despite the absence of lone-pair electrons on the 2'-NH(3)(+) group to coordinate a metal ion. Binding and reaction of a modified oligonucleotide substrate with 2'-NH(2) at U(-1) are unaffected by soft-metal ions. These results suggest that the 2'-OH of U(-1) does not interact with a metal ion. The contribution of the 2'-moiety of U(-1) is unperturbed by thio substitution at either of the nonbridging oxygens of the reactive phosphoryl group, providing no indication of a hydrogen bond between the 2'-OH and the nonbridging phosphoryl oxygens. In contrast, the 10(3)-fold catalytic advantage of 2'-OH relative to 2'-H is eliminated when the 3'-bridging oxygen is replaced by sulfur. As sulfur is a weaker hydrogen-bond acceptor than oxygen, this effect suggests a hydrogen-bonding interaction between the 2'-OH and the 3'-bridging oxygen. CONCLUSIONS: These results provide the first experimental support for the model in which the 2'-OH of U(-1) donates a hydrogen bond to the neighboring 3'-bridging oxygen, thereby stabilizing the developing negative charge on the 3'-oxygen in the transition state.