Quantum Mechanical/Molecular Mechanical Study on the Mechanism of the Enzymatic Baeyer–Villiger Reaction
- 25 January 2012
- journal article
- research article
- Published by American Chemical Society (ACS) in Journal of the American Chemical Society
- Vol. 134 (5) , 2732-2741
Abstract
We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer–Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP+, while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme–reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP+ cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP+ cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.Keywords
This publication has 92 references indexed in Scilit:
- Snapshots of Enzymatic Baeyer-Villiger CatalysisJournal of Biological Chemistry, 2011
- Baeyer−Villiger Monooxygenases: More Than Just Green ChemistryChemical Reviews, 2011
- Induced allostery in the directed evolution of an enantioselective Baeyer–Villiger monooxygenaseProceedings of the National Academy of Sciences, 2010
- CHARMM: The biomolecular simulation programJournal of Computational Chemistry, 2009
- H++: a server for estimating pKas and adding missing hydrogens to macromoleculesNucleic Acids Research, 2005
- All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of ProteinsThe Journal of Physical Chemistry B, 1998
- Recombinant Baker's Yeast as a Whole-Cell Catalyst for Asymmetric Baeyer−Villiger OxidationsJournal of the American Chemical Society, 1998
- DL_POLY_2.0: A general-purpose parallel molecular dynamics simulation packageJournal of Molecular Graphics, 1996
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983
- Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozymeJournal of Molecular Biology, 1976