In Silico Prediction of Cytochrome P450-Mediated Biotransformations of Xenobiotics: A Case Study of Epoxidation.

Predicting the biotransformation of xenobiotics is important in toxicology; however, as more compounds are synthesized than can be investigated experimentally, powerful computational methods are urgently needed to prescreen potentially useful candidates. Cytochrome P450 enzymes (P450s) are the major enzymes involved in xenobiotic metabolism, and many substances are bioactivated by P450s to form active compounds. An example is the conversion of olefinic substrates to epoxides, which are intermediates in the metabolic activation of many known or suspected carcinogens. We have calculated the activation energies for epoxidation by the active species of P450 enzymes (an iron-oxo porphyrin cation radical oxidant, compound I) for a diverse set of 36 olefinic substrates with state-of-the-art density functional theory (DFT) methods. Activation energies can be estimated by the computationally less demanding method of calculating the ionization potentials of the substrates, which provides a useful and simple predictive model based on the reaction mechanism; however, the preclassification of these diverse substrates into weakly polar and strongly polar groups is a prerequisite for the construction of specific predictive models with good predictability for P450 epoxidation. This approach has been supported by both internal and external validations. Furthermore, the relation between the activation energies for the regioselective epoxidation and hydroxylation reactions of P450s and experimental data has been investigated. The results show that the computational method used in this work, single-point energy calculations with the B3LYP functional including zero-point energy and solvation and dispersion corrections based on B3LYP-optimized geometries, performs well in reproducing the experimental trends of the epoxidation and hydroxylation reactions.