Trends in predicted chemoselectivity of cytochrome P450 oxidation: B3LYP barrier heights for epoxidation and hydroxylation reactions.

Prediction of epoxide formation in drug metabolism is a difficult but important task, as epoxide formation is linked to drug toxicity. A comparison of the energy barriers for cytochrome P450 mediated epoxidation of alkenes to the barriers for the hydroxylation of an aliphatic carbon atom next to a double bond has been performed using B3LYP and B3LYP-D3. Relevant experimental data on oxidation selectivity has also been assessed. The results show that density functional theory, when using B3LYP-D3, does well in reproducing the experimental trends. Considering that the comparison involves chemical steps with quite different features this is remarkable. We also find that B3LYP consistently underestimates the hydrogen abstraction barriers relative to the epoxidation barriers, and that including a dispersion correction reduces this problem.