Insights into enzyme catalysis from QM/MM modelling: transition state stabilization in chorismate mutase

Chorismate mutase provides an important test of theories of enzyme catalysis, and of modelling methods. The Claisen rearrangement of chorismate to prephenate in the enzyme has been modelled here by a combined quantum mechanics/molecular mechanics (QM/MM) method. Several pathways have been calculated. The sensitivity of the results to details of model preparation and pathway calculation is tested, and the results are compared in detail to previous similar studies and experiments. The potential energy barrier for the enzyme reaction is estimated at 24.5-31.6 kcal mol(-1) (AM1/CHARMM), and 2.7-11.9 kcal mol(-1) with corrections (e.g. B3LYP/6-31+G(d)). In agreement with previous studies, the present analysis of the calculated paths provides unequivocal evidence of significant transition state stabilization by the enzyme, indicating that this is central to catalysis by the enzyme. The active site is exquisitely complementary to the transition state, stabilizing it more than the substrate, so reducing the barrier to reaction. A number of similar pathways for reaction exist in the protein, as expected. Small structural differences give rise to differences in energetic contributions. Major electrostatic contributions to transition state stabilization come in all cases from Arg90, Arg7, one or two water molecules, and Glu78 (Glu78 destabilizes the transition state less than the substrate), while Arg63 contributes significantly in one model.