Quantum Mechanical/Molecular Mechanical Investigation of the Mechanism of C−H Hydroxylation of Camphor by Cytochrome P450cam: Theory Supports a Two-State Rebound Mechanism

Investigation on spin-flip reaction of Re + CH3CN by relativistic density functional theory

To explore the integrated reaction mechanisms for Re atom with acetonitrile theoretically, density functional theory with zero-order regular approximation (ZORA) relativistic corrections has been employed at the BP86/TZ2P level. There have been three adiabatic potential energy surfaces in the study along sextet, quartet and doublet spin states. However, the detailed minimum energy reaction pathway altogether contains six stationary states () to (), five transition states (), and two intersystem crossings with spin inversion (marked by ⇒): (6)Re + CH3CN → η(1)-ReNCCH3 () → ⇒ η(2)-Re(NC)CH3 () → → η(3)-HRe(NCCH2) () → → CH3-ReNC () → → CH2[double bond, length as m-dash]Re(H)NC () ⇒ → CH[triple bond, length as m-dash]Re(H)2NC (). Thereinto, the lowest energy crossing points (LECP) have been determined by the DFT fractional-occupation-number (FON) approach. The first spin inversion has transferred the potential energy surfaces from high-spin sextet to the quartet intermediate () with the subsequent C-C bond breakage. The second one from the quartet to the low-spin doublet state accompanies the C-H activation, decreasing the transition barrier by 157 kJ mol(-1). The overall reaction could be exothermic by about 210 kJ mol(-1). Harmonic vibration frequencies and NBO, WBO analysis are also applied to verified the experimental observed information.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

Unified view of oxidative C-H bond cleavage and sulfoxidation by a nonheme iron(IV)-oxo complex via Lewis acid-promoted electron transfer

Oxidative C-H bond cleavage of toluene derivatives and sulfoxidation of thioanisole derivatives by a nonheme iron(IV)-oxo complex, [(N4Py)Fe(IV)(O)](2+) (N4Py = N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine), were remarkably enhanced by the presence of triflic acid (HOTf) and Sc(OTf)3 in acetonitrile at 298 K. All the logarithms of the observed second-order rate constants of both the oxidative C-H bond cleavage and sulfoxidation reactions exhibit remarkably unified correlations with the driving forces of proton-coupled electron transfer (PCET) and metal ion-coupled electron transfer (MCET) in light of the Marcus theory of electron transfer when the differences in the formation constants of precursor complexes between PCET and MCET were taken into account, respectively. Thus, the mechanisms of both the oxidative C-H bond cleavage of toluene derivatives and sulfoxidation of thioanisole derivatives by [(N4Py)Fe(IV)(O)](2+) in the presence of HOTf and Sc(OTf)3 have been unified as the rate-determining electron transfer, which is coupled with binding of [(N4Py)Fe(IV)(O)](2+) by proton (PCET) and Sc(OTf)3 (MCET). There was no deuterium kinetic isotope effect (KIE) on the oxidative C-H bond cleavage of toluene via the PCET pathway, whereas a large KIE value was observed with Sc(OTf)3, which exhibited no acceleration of the oxidative C-H bond cleavage of toluene. When HOTf was replaced by DOTf, an inverse KIE (0.4) was observed for PCET from both toluene and [Ru(II)(bpy)3](2+) (bpy =2,2'-bipyridine) to [(N4Py)Fe(IV)(O)](2+). The PCET and MCET reactivities of [(N4Py)Fe(IV)(O)](2+) with Brønsted acids and various metal triflates have also been unified as a single correlation with a quantitative measure of the Lewis acidity.

Effect of protein environment within cytochrome P450cam evaluated using a polarizable-embedding QM/MM method

Metalloenzymes accommodate cofactors and substrates in their active sites, thereby exerting powerful catalytic effects. Understanding the key elements of the mechanism via which such binding is accomplished using a number of atoms in a protein is a fundamental challenge. To address this issue computationally, here we used mechanical-embedding (ME), electronic-embedding (EE), and polarizable-embedding (PE) hybrid quantum mechanics and molecular mechanics (QM/MM) methods and performed an energy decomposition analysis (EDA) of the nonbonding protein environmental effect in the "compound I" intermediate state of cytochrome P450cam. The B3LYP and AMBER99/QP302 methods were used to deal with the QM and MM subsystems, respectively, and the nonbonding interaction energy between these subsystems was decomposed into electrostatic, van der Waals, and polarization contributions. The PE-QM/MM calculation was performed using polarizable force fields that were capable of describing induced dipoles within the MM subsystem, which arose in response to the electric field generated by QM electron density, QM nuclei, and MM point charges. The present QM/MM EDA revealed that the electrostatic term constituted the largest stabilizing interaction between the QM and MM subsystems. When proper adjustment was made for the point charges of the MM atoms located at the QM-MM boundary, EE-QM/MM and PE-QM/MM calculations yielded similar QM electron density distributions, indicating that the MM polarization effect does not have a significant influence on the extent of QM polarization in this particular enzyme system.

QM/MM study of the C-C coupling reaction mechanism of CYP121, an essential cytochrome p450 of Mycobacterium tuberculosis

Among 20 p450s of Mycobacterium tuberculosis (Mt), CYP121 has received an outstanding interest, not only due to its essentiality for bacterial viability but also because it catalyzes an unusual carbon-carbon coupling reaction. Based on the structure of the substrate bound enzyme, several reaction mechanisms were proposed involving first Tyr radical formation, second Tyr radical formation, and C-C coupling. Key and unknown features, being the nature of the species that generate the first and second radicals, and the role played by the protein scaffold each step. In the present work we have used classical and quantum based computer simulation methods to study in detail its reaction mechanism. Our results show that substrate binding promotes formation of the initial oxy complex, Compound I is the responsible for first Tyr radical formation, and that the second Tyr radical is formed subsequently, through a PCET reaction, promoted by the presence of key residue Arg386. The final C-C coupling reaction possibly occurs in bulk solution, thus yielding the product in one oxygen reduction cycle. Our results thus contribute to a better comprehension of MtCYP121 reaction mechanism, with direct implications for inhibitor design, and also contribute to our general understanding of these type of enzymes.

Use of density functional theory in drug metabolism studies

INTRODUCTION

The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes.

AREAS COVERED

We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom.

EXPERT OPINION

DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.

Spin-flip reactions of Zr + C2H6 researched by relativistic density functional theory

Density functional theory (DFT) with relativistic corrections of zero-order regular approximation (ZORA) has been applied to explore the reaction mechanisms of ethane dehydrogenation by Zr atom with triplet and singlet spin-states. Among the complicated minimum energy reaction path, the available states involves three transition states (TS), and four stationary states (1) to (4) and one intersystem crossing with spin-flip (marked by -->): (3) Zr + C 2 H 6 → (3) Zr-CH 3 -CH 3 ((3)1) → (3)TS 1/2 → (3) ZrH-CH 2 -CH 3 ((3)2) → (3) TS 2/3 --> (1) ZrH2-CH2 = CH2 ((1) 3) → (1) TS 3/4 → (1) ZrH 3 -CH = CH 2 ((1)4). The minimum energy crossing point is determined with the help of the DFT fractional-occupation-number (FON) approach. The spin inversion leads the reaction pathway transferring from the triplet potential energy surface (PES) to the singlet's accompanying with the activation of the second C-H bond. The overall reaction is calculated to be exothermic by about 231 kJ mol(-1). Frequency and NBO analysis are also applied to confirm with the experimental observed data.

Water complexes of cytochrome P450: insights from energy decomposition analysis

Water is a small molecule that nevertheless perturbs, sometimes significantly, the electronic properties of an enzyme’s active site. In this study, interactions of a water molecule with the ferric heme and the compound I (Cpd I) intermediate of cytochrome P450 are studied. Energy decomposition analysis (EDA) schemes are used to investigate the physical origins of these interactions. Localized molecular orbital EDA (LMOEDA) implemented in the quantum chemistry software GAMESS and the EDA method implemented in the ADF quantum chemistry program are used. EDA reveals that the electrostatic and polarization effects act as the major driving force in both of these interactions. The hydrogen bonding in the Cpd I•••H₂O complex is similar to that in the water dimer; however, the relative importance of the electrostatic effect is somewhat larger in the water dimer.

Quantum mechanics/molecular mechanics modeling of regioselectivity of drug metabolism in cytochrome P450 2C9

Cytochrome P450 enzymes (P450s) are important in drug metabolism and have been linked to adverse drug reactions. P450s display broad substrate reactivity, and prediction of metabolites is complex. QM/MM studies of P450 reactivity have provided insight into important details of the reaction mechanisms and have the potential to make predictions of metabolite formation. Here we present a comprehensive study of the oxidation of three widely used pharmaceutical compounds (S-ibuprofen, diclofenac, and S-warfarin) by one of the major drug-metabolizing P450 isoforms, CYP2C9. The reaction barriers to substrate oxidation by the iron-oxo species (Compound I) have been calculated at the B3LYP-D/CHARMM27 level for different possible metabolism sites for each drug, on multiple pathways. In the cases of ibuprofen and warfarin, the process with the lowest activation energy is consistent with the experimentally preferred metabolite. For diclofenac, the pathway leading to the experimentally observed metabolite is not the one with the lowest activation energy. This apparent inconsistency with experiment might be explained by the two very different binding modes involved in oxidation at the two competing positions. The carboxylate of diclofenac interacts strongly with the CYP2C9 Arg108 side chain in the transition state for formation of the observed metabolite—but not in that for the competing pathway. We compare reaction barriers calculated both in the presence and in the absence of the protein and observe a marked improvement in selectivity prediction ability upon inclusion of the protein for all of the substrates studied. The barriers calculated with the protein are generally higher than those calculated in the gas phase. This suggests that active-site residues surrounding the substrate play an important role in controlling selectivity in CYP2C9. The results show that inclusion of sampling (particularly) and dispersion effects is important in making accurate predictions of drug metabolism selectivity of P450s using QM/MM methods.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

How Do Perfluorinated Alkanoic Acids Elicit Cytochrome P450 to Catalyze Methane Hydroxylation? An MD and QM/MM Study

Recent experimental studies show that usage of perfluoro decanoic acid (PFDA), as a dummy substrate, can elicit P450_BM3 to perform hydroxylation of small alkanes, such as methane (ref. 17) and propane (ref. 17 and ref. 18). To comprehend the mechanism whereby PFDA operates to potentiate P450_BM3 to catalyze the hydroxylation of small alkanes, we used molecular dynamics (MD) and hybrid quantum mechanical / molecular mechanical (QM/MM) calculations. The MD results show that without the PFDA, methane escapes the active site, while the presence of PFDA can potentially induce a productive Cpd I-Methane juxtaposition for rapid oxidation. Nevertheless, when only a single methane molecule is present near the PFDA, it still escapes the pocket within less than a nanosecond. However, when three methane molecules are present in the pocket, they alternate quasi-periodically such that at all times (within 10 ns), a molecule of methane is always present in the proximity of Cpd I in a reactive conformation. Our results further demonstrate that the PFDA does not exert any electrostatic catalysis, whether the PFDA is in the protonated or deprotonated forms. Taken together, we conclude that methane hydroxylation requires, in addition to PFDA, a high partial pressure of methane that will cause a high methane concentration in the active site. Further study of ethane and propane hydroxylations demonstrates that higher alkane concentration is helpful for all the three small alkanes. Thus for the smallest alkane, methane, at least three molecules are necessary whereas for the larger ethane, two molecules are needed to force one ethane to be closer to Cpd I. Finally, for propane a second molecule is helpful but not absolutely necessary; for this molecule the PFDA may well be sufficient to keep propane close to Cpd I for efficient oxidation. We therefore propose that high alkane pressure should assist small alkane hydroxylation by P450 in a manner inversely proportional to the size of the alkanes.

Effects of Dispersion in Density Functional Based Quantum Mechanical/Molecular Mechanical Calculations on Cytochrome P450 Catalyzed Reactions

Density functional theory (DFT) based quantum mechanical/molecular mechanical (QM/MM) calculations have provided valuable insight into the reactivity of the cytochrome P450 family of enzymes (P450s). A failure of commonly used DFT methods, such as B3LYP, is the neglect of dispersion interactions. An empirical dispersion correction has been shown to improve the accuracy of gas phase DFT calculations of P450s. The current work examines the effect of the dispersion correction in QM/MM calculations on P450s. The hydrogen abstraction from camphor, and hydrogen abstraction and C-O addition of cyclohexene and propene by P450cam have been modeled, along with the addition of benzene to Compound I in CYP2C9, at the B3LYP-D2/CHARMM27 level of theory. Single point energy calculations were also performed at the B3LYP-D3//B3LYP-D2/CHARMM27 level. The dispersion corrections lower activation energy barriers significantly (by ∼5 kcal/mol), as seen for gas phase calculations, but has a small effect on optimized geometries.These effects are likely to be important in modeling reactions catalyzed by other enzymes also. Given the low computational cost of including such dispersion corrections, we recommend doing so in all B3LYP based QM/MM calculations.

Synthesis and catalytic properties of a series of cobalt porphyrins as cytochrome P450 model: the effect of substituents on the catalytic activity

A series of cobalt porphyrins derived from hemin was prepared as cytochrome P450 models. Effects of substituents at the cobalt deuteroporphyrin-propionate side chains are investigated in oxidation of toluene with air to benzaldehyde and benzyl alcohol without the use of solvent and sacrificial co-reductant. The catalytic activity of cobalt porphyrins depends on the type of substituents. When the electron-withdrawing groups like -Cl, -Br, were introduced into the double propionate side chains, they can increase the catalyst stability and selectivity to benzaldehyde. In comparison with these electron-withdrawing groups, the electron-donor groups, such as -CH₃, -S-S- and -NH₂ groups, can improve their catalytic activities. Moreover, the electron-donor group containing an unpaired electron (such as -S-S-, -NH₂) is benefit for improving its catalytic efficiency and promoting the electron delivery. It can be concluded that the double propionate side chains in the deuteroporphyrin complex may participate in oxidation process and effect electron transfer from the high-valent metalloporphyrin species to the substrate.

Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms

Metabolism of xenobiotics remains a central challenge for the discovery and development of drugs, cosmetics, nutritional supplements, and agrochemicals. Metabolic transformations are frequently related to the incidence of toxic effects that may result from the emergence of reactive species, the systemic accumulation of metabolites, or by induction of metabolic pathways. Experimental investigation of the metabolism of small organic molecules is particularly resource demanding; hence, computational methods are of considerable interest to complement experimental approaches. This review provides a broad overview of structure- and ligand-based computational methods for the prediction of xenobiotic metabolism. Current computational approaches to address xenobiotic metabolism are discussed from three major perspectives: (i) prediction of sites of metabolism (SOMs), (ii) elucidation of potential metabolites and their chemical structures, and (iii) prediction of direct and indirect effects of xenobiotics on metabolizing enzymes, where the focus is on the cytochrome P450 (CYP) superfamily of enzymes, the cardinal xenobiotics metabolizing enzymes. For each of these domains, a variety of approaches and their applications are systematically reviewed, including expert systems, data mining approaches, quantitative structure–activity relationships (QSARs), and machine learning-based methods, pharmacophore-based algorithms, shape-focused techniques, molecular interaction fields (MIFs), reactivity-focused techniques, protein–ligand docking, molecular dynamics (MD) simulations, and combinations of methods. Predictive metabolism is a developing area, and there is still enormous potential for improvement. However, it is clear that the combination of rapidly increasing amounts of available ligand- and structure-related experimental data (in particular, quantitative data) with novel and diverse simulation and modeling approaches is accelerating the development of effective tools for prediction of in vivo metabolism, which is reflected by the diverse and comprehensive data sources and methods for metabolism prediction reviewed here. This review attempts to survey the range and scope of computational methods applied to metabolism prediction and also to compare and contrast their applicability and performance.

Hydroxylation catalysis by mononuclear and dinuclear iron oxo catalysts: a methane monooxygenase model system versus the Fenton reagent Fe(IV)O(H2O)5(2+)

Hydroxylation of aliphatic C-H bonds is a chemically and biologically important reaction, which is catalyzed by the oxidoiron group FeO(2+) in both mononuclear (heme and nonheme) and dinuclear complexes. We investigate the similarities and dissimilarities of the action of the FeO(2+) group in these two configurations, using the Fenton-type reagent [FeO(2+) in a water solution, FeO(H(2)O)(5)(2+)] and a model system for the methane monooxygenase (MMO) enzyme as representatives. The high-valent iron oxo intermediate MMOH(Q) (compound Q) is regarded as the active species in methane oxidation. We show that the electronic structure of compound Q can be understood as a dimer of two Fe(IV)O(2+) units. This implies that the insights from the past years in the oxidative action of this ubiquitous moiety in oxidation catalysis can be applied immediately to MMOH(Q). Electronically the dinuclear system is not fundamentally different from the mononuclear system. However, there is an important difference of MMOH(Q) from FeO(H(2)O)(5)(2+): the largest contribution to the transition state (TS) barrier in the case of MMOH(Q) is not the activation strain (which is in this case the energy for the C-H bond lengthening to the TS value), but it is the steric hindrance of the incoming CH(4) with the ligands representing glutamate residues. The importance of the steric factor in the dinuclear system suggests that it may be exploited, through variation in the ligand framework, to build a synthetic oxidation catalyst with the desired selectivity for the methane substrate.

A highly reactive mononuclear non-heme manganese(IV)-oxo complex that can activate the strong C-H bonds of alkanes

A mononuclear non-heme manganese(IV)-oxo complex has been synthesized and characterized using various spectroscopic methods. The Mn(IV)-oxo complex shows high reactivity in oxidation reactions, such as C-H bond activation, oxidations of olefins, alcohols, sulfides, and aromatic compounds, and N-dealkylation. In C-H bond activation, the Mn(IV)-oxo complex can activate C-H bonds as strong as those in cyclohexane. It is proposed that C-H bond activation by the non-heme Mn(IV)-oxo complex does not occur via an oxygen-rebound mechanism. The electrophilic character of the non-heme Mn(IV)-oxo complex is demonstrated by a large negative ρ value of -4.4 in the oxidation of para-substituted thioanisoles.

Quantum mechanically derived AMBER-compatible heme parameters for various states of the cytochrome P450 catalytic cycle

Molecular mechanics (MM) methods are computationally affordable tools for screening chemical libraries of novel compounds for sites of P450 metabolism. One challenge for MM methods has been the absence of a consistent and transferable set of parameters for the heme within the P450 active site. Experimental data indicate that mammalian P450 enzymes vary greatly in the size, architecture, and plasticity of their active sites. Thus, obtaining X-ray-based geometries for the development of accurate MM parameters for the major classes of hepatic P450 remains a daunting task. Our previous work with preliminary gas-phase quantum mechanics (QM)-derived atomic partial charges greatly improved the accuracy of docking studies of raloxifene to CYP3A4. We have therefore developed and tested a consistent set of transferable MM parameters based on gas-phase QM calculations of two model systems of the heme-a truncated (T-HM) and a full (F-HM) for four states of the P450 catalytic cycle. Our results indicate that the use of the atomic partial charges from the F-HM further improves the accuracy of docked predictions for raloxifene to CYP3A4. Different patterns for substrate docking are also observed depending on the choice of heme model and state. Newly parameterized heme models are tested in implicit and explicitly solvated MD simulations in the absence and presence of enzyme structures, for CYP3A4, and appear to be stable on the nanosecond simulation timescale. The new force field for the various heme states may aid the community for simulations of P450 enzymes and other heme-containing enzymes.

The rate-limiting step in P450 hydroxylation of hydrocarbons a direct comparison of the "somersault" versus the "consensus" mechanism involving compound I

Model theoretical quantum mechanical (QM) calculations are described for the P-450 hydroxylation of methane, isobutane, and camphor that compare the concerted somersault H-abstraction mechanism with the oxidation step involving Cpd I. Special emphasis has been placed on maintaining a balanced basis set in the oxidation step. QM calculations, employing the 6-311+G(d,p) basis set on the Fe atom and all of the key surrounding atoms involved in the C-H abstraction step, reaffirm a mechanism involving rearrangement of the iron hydroperoxide group (FeO-OH --> FeO...HO(*)) in concert with hydrogen abstraction from the C-H bond of the substrate by the incipient bound hydroxyl radical HO(*). The barrier for the somersault rearrangement of model Cpd 0 (FeO-OH) is calculated to be 21.4 kcal/mol in the absence of substrate. The overall activation energy for the oxidation of camphor involving the somersault motion of the FeO-OH group of P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] --> [Por(SH)Fe(III)-O....HO(-)] in concert with hydrogen abstraction is DeltaE(++) = 12.4 kcal/mol. The corresponding abstraction of the hydrogen atom from the C-H bond of camphor by Cpd I has an activation barrier of 17.6 kcal/mol. Arguments are presented that the somersault rearrangement is induced by steric compression at the active site. Kinetic isotope effect data are discussed that provides compelling evidence for a rate-limiting step involving C-H bond cleavage.

What factors influence the rate constant of substrate epoxidation by compound I of cytochrome P450 and analogous iron(IV)-oxo oxidants?

The cytochromes P450 are a versatile range of mono-oxygenase enzymes that catalyze a variety of different chemical reactions, of which the key reactions include aliphatic hydroxylation and C=C double bond epoxidation. To establish the fundamental factors that govern substrate epoxidation by these enzymes we have done a systematic density functional theory study on substrate epoxidation by the active species of P450 enzymes, namely the iron(IV)-oxo porphyrin cation radical oxidant or Compound I. We show here, for the first time, that the rate constant of substrate epoxidation, and hence the activation energy, correlates with the ionization potential of the substrate as well as with intrinsic electronic properties of the active oxidant such as the polarizability volume. To explain these findings we present an electron-transfer model for the reaction mechanism that explains the factors that determine the barrier heights and developed a valence bond (VB) curve crossing mechanism to rationalize the observed trends. In addition, we have found a correlation for substrate epoxidation reactions catalyzed by a range of heme and nonheme iron(IV)-oxo oxidants with the strength of the O-H bond in the iron-hydroxo complex, i.e. BDE(OH), which is supported by the VB model. Finally, the fundamental factors that determine the regioselectivity change between substrate hydroxylation and epoxidation are discussed. It is shown that the regioselectivity of aliphatic hydroxylation versus double bond epoxidation is not influenced by the choice of the oxidant but is purely substrate dependent.