Multiple high-level QM/MM reaction paths demonstrate transition-state stabilization in chorismate mutase: correlation of barrier height with transition-state stabilization

Multiple profiles for the reaction from chorismate to prephenate in the enzyme chorismate mutase calculated with hybrid density functional combined quantum mechanics/molecular mechanics methods (B3LYP/6-31G(d)-CHARMM27) agree well with experiment, and provide direct evidence of transition-state stabilization by this important enzyme, which is at the centre of current debates about the nature of enzyme catalysis.     

All electron quantum chemical calculation of the entire enzyme system confirms a collective catalytic device in the chorismate mutase reaction

To elucidate the catalytic power of enzymes, we analyzed the reaction profile of Claisen rearrangement of Bacillus subtilis chorismate mutase (BsCM) by all electron quantum chemical calculations using the fragment molecular orbital (FMO) method. To the best of our knowledge, this is the first report of ab initio-based quantum chemical calculations of the entire enzyme system, where we provide a detailed analysis of the catalytic factors that accomplish transition-state stabilization (TSS). FMO calculations deliver an ab initio-level estimate of the intermolecular interaction between the substrate and the amino acid residues of the enzyme. To clarify the catalytic role of Arg90, we calculated the reaction profile of the wild-type BsCM as well as Lys90 and Cit90 mutant BsCMs. Structural refinement and the reaction path determination were performed at the ab initio QM/MM level, and FMO calculations were applied to the QM/MM refined structures. Comparison between three types of reactions established two collective catalytic factors in the BsCM reaction: (1) the hydrogen bonds connecting the Glu78-Arg90-substrate cooperatively control the stability of TS relative to the ES complex and (2) the positive charge on Arg90 polarizes the substrate in the TS region to gain more electrostatic stabilization.     

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.     

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.     

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     

A transition path sampling study of the reaction catalyzed by the enzyme chorismate mutase

The study of the chemical steps in enzyme-catalyzed reactions represents a challenge for molecular simulation techniques. One concern is how to calculate paths for the reaction. Common techniques include the definition of a reaction coordinate in terms of a small set of (normally) geometrical variables or the determination of minimum energy paths on the potential energy surface of the reacting system. Both have disadvantages, the former because it presupposes knowledge of which variables are likely to be important for reaction and the latter because it provides a static picture and dynamical effects are ignored. In this paper, we employ the transition path sampling method developed by Chandler and co-workers, which overcomes some of these limitations. The reaction that we have chosen is the chorismate-mutase-catalyzed conversion of chorismate into prephenate, which has become something of a test case for simulation studies of enzyme mechanisms. We generated an ensemble of approximately 1000 independent transition paths for the reaction in the enzyme and another approximately 500 for the corresponding reaction in solution. A large variety of analyses of these paths was performed, but we have concentrated on characterizing the transition state ensemble, particularly the flexibility of its structures with respect to other ligands of the enzyme and the time evolution of various geometrical and energetic properties as the reaction proceeds. We have also devised an approximate technique for locating transition state structures along the paths.     

Conformational effects in enzyme catalysis: QM/MM free energy calculation of the 'NAC' contribution in chorismate mutase

The controversial 'near attack conformation'(NAC) effect in the important model enzyme chorismate mutase is calculated to be 3.8-4.6 kcal mol(-1) by QM/MM free energy perturbation molecular dynamics methods, showing that the NAC effect by itself does not account for catalysis in this enzyme.     

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.     

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.     

Exploring the quantum mechanical/molecular mechanical replica path method: a pathway optimization of the chorismate to prephenate Claisen rearrangement catalyzed by chorismate mutase

 A replica path method has been developed and extended for use in complex systems involving hybrid quantum/classical (quantum mechanical/molecular mechanical) coupled potentials. This method involves the definition of a reaction path via replication of a set of macromolecular atoms. An “important” subset of these replicated atoms is restrained with a penalty function based on weighted root-mean-square rotation/translation best-fit distances between adjacent (i±1) and next adjacent (i±2) pathway steps. An independent subset of the replicated atoms may be treated quantum mechanically using the computational engine Gamess-UK. This treatment can be performed in a highly parallel manner in which many dozens of processors can be efficiently employed. Computed forces may be projected onto a reference pathway and integrated to yield a potential of mean force (PMF). This PMF, which does not suffer from large errors associated with calculated potential-energy differences, is extremely advantageous. As an example, the QM/MM replica path method is applied to the study of the Claisen rearrangement of chorismate to prephenate which is catalyzed by the Bacillus subtilis isolated, chorismate mutase. Results of the QM/MM pathway minimizations yielded an activation enthalpy ΔH ^†† of 14.9 kcal/mol and a reaction enthalpy of −19.5 kcal/mol at the B3LYP/6-31G(d) level of theory. The resultant pathway was compared and contrasted with one obtained using a forced transition approach based on a reaction coordinate constrained repeated walk procedure (ΔH ^†† =20.1 kcal/mol, ΔH _rxn = −20.1 kcal/mol, RHF/4-31G). The optimized replica path results compare favorably to the experimental activation enthalpy of 12.7±0.4 kcal/mol. Received: 16 December 2001 / Accepted: 6 September 2002 / Published online: 8 April 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 22nd National Meeting of the American Chemical Society, 2001. Correspondence to: H.L. Woodcock e-mail: hlwood@ccqc.uga.edu Acknowledgements. The authors thank Eric Billings, Xiongwu Wu, and Stephen Bogusz for helpful discussions and related work. The authors also show grateful appreciation to The National Institutes of Health and The National Science Foundation for support of the current research.     

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.     

Electrostatic effects in enzyme catalysis: a quantum mechanics/molecular mechanics study of the nucleophilic substitution reaction in haloalkane dehalogenase

Combined quantum mechanics/molecular mechanics molecular dynamics simulations have been carried out to study the cleavage of the carbon–chlorine bond in 1,2-dichloroethane catalysed by haloalkane dehalogenase from Xanthobacter Autotrophicus GJ10. The process has been compared with an adequate counterpart in aqueous solution, the nucleophilic attack of acetate anion on 1,2-dichloroethane. Within the limitations of the model, mainly due to the use of a semiempirical Hamiltonian, our results reproduce the magnitude and characteristics of the catalytic effect. Comparisons of the enzymatic and in solution potentials of mean force reveal that, irrespective of the reference state, the enzyme shows a larger affinity for the transition state. The origin of this increased affinity is found in the differences in the electrostatic pattern created by the environment in aqueous solution and in the enzyme.Proceedings of the 11th International Congress of Quantum Chemistry satellite meeting in honor of Jean-Louis Rivail     

Decomposition and isomerization of 1,2-benzisoxazole: single-pulse shock-tube experiments, quantum chemical and transition-state theory calculations

Isomerization and decomposition of 1,2-benzisoxazole were studied behind reflected shock waves in a pressurized driver, single-pulse shock tube. It isomerizes to o-hydroxybenzonitrile, and no fragmentation is observed up to a temperature where the isomerization is almost complete (approximately 1040 K at 2 ms reaction time). The isomerization experiments in this investigation covered the temperature range 900-1040 K. The lack of fragmentation is in complete contrast to the thermal behavior of isoxazole, where no isomerization was observed and the main decomposition products over the same temperature range were carbon monoxide and acetonitrile. In a series of experiments covering the temperature range 1190-1350 K, a plethora of fragmentation products appear in the post shock samples of 1,2-benzisoxazole. The product distribution is exactly the same regardless of whether the starting material is 1,2-benzisoxazole or o-hydroxybenzonitrile, indicating that over this temperature range the 1,2-benzisoxazole has completely isomerized to o-hydroxybenzonitrile prior to fragmentation. Two potential energy surfaces that lead to the isomerization were evaluated by quantum chemical calculations. One surface with one intermediate and two transition states has a high barrier and does not contribute to the process. The second surface is more complex. It has three intermediates and four transition states, but it has a lower overall barrier and yields the isomerization product o-hydroxybenzonitrile at a much higher rate. The unimolecular isomerization rate constants kinfinity at a number of temperatures in the range of 900-1040 K were calculated from the potential energy surface using transition-state theory and then expressed in an Arrhenius form. The value obtained is kfirst=4.15x10(14) exp(-51.7x10(3)/RT) s-1, where R is expressed in units of cal/(K mol). The calculated value is somewhat higher than the one obtained from the experimental results. When it is expressed in terms of energy difference it corresponds of ca. 2 kcal/mol.     

Water clusters on graphite: methodology for quantum chemical a priori prediction of reaction rate constants

The properties, interactions, and reactions of cyclic water clusters (H(2)O)(n=1-5) on model systems for a graphite surface have been studied using pure B3LYP, dispersion-augmented density functional tight binding (DFTB-D), and integrated ONIOM(B3LYP:DFTB-D) methods. Coronene C(24)H(12) as well as polycircumcoronenes C(96)H(24) and C(216)H(36) in monolayer, bilayer, and trilayer arrangements were used as model systems to simulate ABA bulk graphite. Structures, binding energies, and vibrational frequencies of water clusters on mono- and bilayer graphite models have been calculated, and structural changes and frequency shifts due to the water cluster-graphite interactions are discussed. ONIOM(B3LYP:DFTB-D) with coronene and water in the high level and C(96)H(24) in the low level mimics the effect of extended graphite pi-conjugation on the water-graphite interaction very reasonably and suggests that water clusters only weakly interact with graphite surfaces, as suggested by the fact that water is an excellent graphite lubricant. We use the ONIOM(B3LYP:DFTB-D) method to predict rate constants for model pathways of water dissociative adsorption on graphite. Quantum chemical molecular dynamics (QM/MD) simulations of water clusters and water addition products on the C(96)H(24) graphite model are presented using the DFTB-D method. A three-stage strategy is devised for a priori investigations of high temperature corrosion processes of graphite surfaces due to interaction with water molecules and fragments.     

Facile S(N)2 reaction in protic solvent: quantum chemical analysis

We study the effects of protic solvent (water, methanol, ethanol, and tert-butyl alcohol) and cation (Na+, K+, Cs+) on the unsymmetrical SN2 reaction X- + RY --> RX + Y- (X = F, Br; R = CH3,C3H7;Y = Cl, OMs). We describe a series of calculations for the S(N)2 reaction mechanism under the influence of cation and protic solvent, presenting the structures of pre- and postreaction complexes and transition states and the magnitude of the activation barrier. An interesting mechanism is proposed, in which the protic solvent molecules that are shielded from the nucleophile by the intervening cation act as a Lewis base to reduce the unfavorable Coulombic influence of the cation on the nucleophile. We predict that the reaction barrier for the S(N)2 reaction is significantly lowered by the cooperative effects of cation and protic solvent. We show that the cation and protic solvent, each of which has been considered to retard the SN2 reactivity of the nucleophile, can accelerate the reaction tremendously when they interact with the fluoride ion in an intricate, combined fashion. This alternative S(N)2 mechanism is discussed in relation to the recently observed phenomenal efficiency of fluorination in tert-alcohol media [Kim, D. W.; et al. J. Am. Chem. Soc. 2006, 128, 16394].