Toward accurate barriers for enzymatic reactions: QM/MM case study on p-hydroxybenzoate hydroxylase

The hydroxylation reaction catalyzed by p-hydroxybenzoate hydroxylase has been investigated by quantum mechanical/molecular mechanical (QM/MM) calculations at different levels of QM theory. The solvated enzyme was modeled (approximately 23,000 atoms in total, 49 QM atoms). The geometries of reactant and transition state were optimized for ten representative pathways using semiempirical (AM1) and density functional (B3LYP) methods as QM components. Single-point calculations at B3LYP/MM optimized geometries were performed with local correlation methods [LMP2, LCCSD(T0)] and augmented triple-zeta basis sets. A careful validation of the latter approach with regard to all computational parameters indicates convergence of the QM contribution to the computed barriers to within approximately 1 kcal mol(-1). Comparison with the available experimental data supports this assessment.