17O Electron Nuclear Double Resonance Analysis of Compound I: Inverse Correlation between Oxygen Spin Population and Electron Donation

Regulation of ferryl reactivity by the cytochrome P450 decarboxylase OleT

The cytochrome P450 OleT catalyzes the decarboxylation of long-chain fatty acid substrates to produce terminal alkenes using hydrogen peroxide as a co-substrate. The facile activation of peroxide to form Compound I in the first step of the reaction, and subsequent CC bond cleavage mediated by Compound II, provides a unique opportunity to visualize both ferryl intermediates using transient kinetic approaches. Analysis of the Arrhenius behavior yields activation barriers of ∼6 kcal/mol and ∼ 18 kcal/mol for the decay of Compound I and Compound II respectively. The influence of the secondary coordination sphere, probed through site-directed mutagenesis approaches, suggests that restriction of the donor-acceptor distance contributes to the reactivity of Compound I. The reactivity of Compound II was further probed using kinetic solvent isotope effect approaches, confirming that the large barrier owes to a proton-gated mechanism in the decarboxylation reaction coordinate. Hydrogen-bonding to an active-site histidine (H85) in the distal pocket plays a key role in this process.

Importance of the Ferryl Quintet State in Determining the Electronic Properties of P450 Compound I

We previously reported a selenolate-ligated P450 compound I intermediate (SeP450-I) to be more reactive toward C-H bonds than its thiolate-ligated counterpart. To gain insight into how the selenolate axial ligand influences the reactivity of compound I, we have investigated the electronic structure of the SeP450-I intermediate using variable temperature Mössbauer (VTM) spectroscopy. The VTM data indicate that electronic spin relaxation rates are significantly slower in SeP450-I than in P450-I. Analyses of these data provide Δ, the energy spacing between the two lowest electronic energy levels in compound I. This spacing is typically determined by the zero-field splitting of the ferryl moiety, D, and the exchange coupling, J, between the iron(IV)oxo unit and the ligand-based radical. However, the systems examined are antiferromagnetically coupled with |J/D| > 1. As a result, Δ ∼ (3/2) J, and measurements of Δ provide J (to within ∼5%). These measurements reveal that the sign and magnitude of J track with the reactivity of compound I toward C-H bonds. Efforts to analyze these and other data highlight the inadequacy of the standard ligand field model that is often used to explain the electronic properties of compound I. Additional analyses combining our data with state energies from a previous theoretical investigation support predictions of a low-lying quintet state within the iron(IV)oxo unit. We discuss these findings in light of computational studies that suggest that access to excited states, particularly those of a high-spin nature, can promote metal-oxo mediated C-H bond cleavage.

Mechanistic Insights into CYP199A4-Catalyzed α-Hydroxyketone Formation and Hydrogen Bond-Assisted C-C Bond Cleavage Catalyzed by the CYP199A4 F182L Mutant

CYP199A4 is a cytochrome P450 and can catalyze the hydroxylation of 4-propionylbenzoic acid (4-pIBA) to generate α-hydroxyketone with high stereoselectivity. The F182L mutant of CYP199A4 (F182L-CYP199A4) has been shown to support the cleavage of the C-C bond between the carbonyl and hydroxyl groups of α-hydroxyketone, whereas wild-type CYP199A4 cannot. To uncover how the Phe182 regulates substrate reactivity, we conducted classical molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) MD simulations on these systems. The results predicted that the formation of α-hydroxyketone preferentially led to the (S)-enantiomer. Moreover, the findings revealed that the F182L-CYP199A4 facilitated the formation of a hydrogen bond between the α-hydroxyketone and the reactive peroxoanion (POA) species. This interaction stabilized the α-hydroxyketone near POA and promoted the subsequent C-C bond cleavage. The mechanism of α-hydroxyketone formation and the subsequent C-C bond cleavage were elucidated by employing the hybrid density functional theory (DFT). The α-hydroxyketone formation mechanism involved C-H hydroxylation of 4-pIBA with a rate-limiting energy barrier of 17.1 kcal/mol. The C-C bond cleavage of α-hydroxyketone catalyzed by F182L-CYP199A4 occurred via a radical attack mechanism.

Mechanism of the Oxidative Ring-Closure Reaction during Gliotoxin Biosynthesis by Cytochrome P450 GliF

During gliotoxin biosynthesis in fungi, the cytochrome P450 GliF enzyme catalyzes an unusual C-N ring-closure step while also an aromatic ring is hydroxylated in the same reaction cycle, which may have relevance to drug synthesis reactions in biotechnology. However, as the details of the reaction mechanism are still controversial, no applications have been developed yet. To resolve the mechanism of gliotoxin biosynthesis and gain insight into the steps leading to ring-closure, we ran a combination of molecular dynamics and density functional theory calculations on the structure and reactivity of P450 GliF and tested a range of possible reaction mechanisms, pathways and models. The calculations show that, rather than hydrogen atom transfer from the substrate to Compound I, an initial proton transfer transition state is followed by a fast electron transfer en route to the radical intermediate, and hence a non-synchronous hydrogen atom abstraction takes place. The radical intermediate then reacts by OH rebound to the aromatic ring to form a biradical in the substrate that, through ring-closure between the radical centers, gives gliotoxin products. Interestingly, the structure and energetics of the reaction mechanisms appear little affected by the addition of polar groups to the model and hence we predict that the reaction can be catalyzed by other P450 isozymes that also bind the same substrate. Alternative pathways, such as a pathway starting with an electrophilic attack on the arene to form an epoxide, are high in energy and are ruled out.

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

Mechanistic insights on the epoxidation of alkenes by high-valent non-heme Fe(iv) and Fe(v) oxidants: a comparative theoretical study

This work is based on the screening of better high-valent oxidants, and also includes a mechanistic study during oxygen atom transfer reactions.

The role of basicity in selective C-H bond activation by transition metal-oxidos

The development of (bio)catalysts capable of selectively activating strong C-H bonds is a continuing challenge in modern chemistry. In both metalloenzymes and synthetic systems capable of activating C-H bonds, transition metal-oxido intermediates serve as the active species for reactivity whose thermodynamic properties influence the bond strengths they are capable of activating. In this Frontier article, we present current ideas of how the basicity of transition metal-oxidos impacts their reactivity with C-H bonds and present new opportunities within this field. We highlight recent insights into the role basicity plays in the activation process and its influence on mechanism, as well as the important role that secondary coordination sphere effects, such as hydrogen bonds, in tuning the basicity of the metal-oxido species is discussed.

Quantifiable polarity match effect on C-H bond cleavage reactivity and its limits in reaction design

When oxidants favour cleaving a strong C-H bond at the expense of weaker ones, which are otherwise inherently preferred due to their favourable reaction energy, reactivity factors such as the polarity match effect are often invoked. Polarity match follows the intuition of electrophilic (nucleophilic) oxidants reacting faster with nucleophilic (electrophilic) C-H bonds. Nevertheless, this concept is purely qualitative and is best suited for a posteriori rationalization of experimental observations. Here, we propose and inspect two methods to quantify polar effects in C-H cleavage reactions, one by computation via the difference of atomic charges (Δq) of reacting atoms, and one amenable to experimental measurement through asynchronicity factors, η. By their application to three case studies, we observe that both Δq and η faithfully capture the notion of polarity match. The polarity match model, however, proves insufficient as a predictor of H-atom abstraction reactivity and we discourage its use as a standalone variable in reaction design. Besides this caveat, η and Δq (through its mapping on η) allow the implementation of polarity match into a Marcus-type model of reactivity, alleviating its shortcomings and making reaction planning feasible.

Comparative oxidative ability of mononuclear and dinuclear high-valent iron-oxo species towards the activation of methane: does the axial/bridge atom modulate the reactivity?

Over the years, mononuclear Fe^IVO species have been extensively studied, but the presence of dinuclear Fe^IVO species in soluble methane monooxygenase (sMMO) has inspired the development of biomimic models that could activate inert substrates such as methane. There are some successful attempts; particularly the [(Por)(m-CBA) Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- species has been reported to activate methane and yield decent catalytic turnover numbers and therefore regarded as the closest to the sMMO enzyme functional model, as no mononuclear Fe^IVO analogues could achieve this feat. In this work, we have studied a series of mono and dinuclear models using DFT and ab initio DLPNO-CCSD(T) calculations to probe the importance of nuclearity in enhancing the reactivity. We have probed the catalytic activities of four complexes: [(HO)Fe^IV(O)(Por)]^- (1), [(HO)Fe^IV(O)(Por˙⁺)] (2), μ-oxo dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-O)Fe^IV(O) (Por˙⁺)]^- (3) and N-bridged dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- (4) towards the activation of methane. Additionally, calculations were performed on the mononuclear models [(X)Fe^IV(O)(Por˙⁺)]ⁿ {X = N 4a (n = -2), NH 4b (n = -1) and NH₂4c (n = 0)} to understand the role of nuclearity in the reactivity. DFT calculations performed on species 1-4 suggest an interesting variation among them, with species 1-3 possessing an intermediate spin (S = 1) as a ground state and species 4 possessing a high-spin (S = 2) as a ground state. Furthermore, the two Fe^IV centres in species 3 and 4 are antiferromagnetically coupled, yielding a singlet state with a distinct difference in their electronic structure. On the other hand, species 2 exhibits a ferromagnetic coupling between the Fe^IV and the Por˙⁺ moiety. Our calculations suggest that the higher barriers for the C-H bond activation of methane and the rebound step for species 1 and 3 are very high in energy, rendering them unreactive towards methane, while species 2 and 4 have lower barriers, suggesting their reactivity towards methane. Studies on the system reveal that model 4a has multiple FeN bonds facilitating greater reactivity, whereas the other two models have longer Fe-N bonds and less radical character with steeper barriers. Strong electronic cooperativity is found to be facilitated by the bridging nitride atom, and this cooperativity is suppressed by substituents such as oxygen, rendering them inactive. Thus, our study unravels that apart from enhancing the nuclearity, bridging atoms that facilitate strong cooperation between the metals are required to activate very inert substrates such as methane, and our results are broadly in agreement with earlier experimental findings.