Valence bond modeling of trends in hydrogen abstraction barriers and transition states of hydroxylation reactions catalyzed by cytochrome P450 enzymes

The paper outlines the fundamental factors that govern the mechanisms of alkane hydroxylation by cytochrome P450 and the corresponding barrier heights during the hydrogen abstraction and radical rebound steps of the process. This is done by a combination of density functional theory calculations for 11 alkanes and valence bond (VB) modeling of the results. The energy profiles and transition states for the various steps are reconstructed using VB diagrams (Shaik, S. S. J. Am. Chem. Soc. 1981, 103, 3692-3701. Shaik, S.; Shurki, A. Angew. Chem. Int. Ed. 1999, 38, 586-625.) and the DFT barriers are reproduced by the VB model from raw data based on C-H bond energies. The model explains a variety of other features of P450 hydroxylations: (a) the nature of the polar effect during hydrogen abstraction, (b) the difference between the activation mechanisms leading to the Fe-IV vs the Fe-III electromers, (c) the difference between the gas phase and the enzymatic reaction, and (d) the dependence of the rebound barrier on the spin state. the VB mechanism shows that the active species of the enzyme performs a complex reaction that involves multiple bond making and breakage mechanisms by utilizing an intermediate VB structure that cuts through the high barrier of the principal transformation between reactants and products, thereby mediating the process at a low energy cost The correlations derived in this paper create order and organize the data for a process of a complex and important enzyme. This treatment can be generalized to the reactivity patterns of nonheme systems and synthetic iron-oxo porphyrin reagents.