Oxygen Atom Transfer Using an Iron(IV)-Oxo Embedded in a Tetracyclic N-Heterocyclic Carbene System: How Does the Reactivity Compare to Cytochrome P450 Compound I?

Cyclic iron tetra N-heterocyclic carbenes: synthesis, properties, reactivity, and catalysis

Cyclic iron tetracarbenes are an emerging class of macrocyclic iron N-heterocyclic carbene (NHC) complexes. They can be considered as an organometallic compound class inspired by their heme analogs, however, their electronic properties differ, e.g. due to the very strong σ-donation of the four combined NHCs in equatorial coordination. The ligand framework of iron tetracarbenes can be readily modified, allowing fine-tuning of the structural and electronic properties of the complexes. The properties of iron tetracarbene complexes are discussed quantitatively and correlations are established. The electronic nature of the tetracarbene ligand allows the isolation of uncommon iron(III) and iron(IV) species and reveals a unique reactivity. Iron tetracarbenes are successfully applied in C-H activation, CO₂ reduction, aziridination and epoxidation catalysis and mechanisms as well as decomposition pathways are described. This review will help researchers evaluate the structural and electronic properties of their complexes and target their catalyst properties through ligand design.

Product Distributions of Cytochrome P450 OleTJE with Phenyl-Substituted Fatty Acids: A Computational Study

There are two types of cytochrome P450 enzymes in nature, namely, the monooxygenases and the peroxygenases. Both enzyme classes participate in substrate biodegradation or biosynthesis reactions in nature, but the P450 monooxygenases use dioxygen, while the peroxygenases take H₂O₂ in their catalytic cycle instead. By contrast to the P450 monooxygenases, the P450 peroxygenases do not require an external redox partner to deliver electrons during the catalytic cycle, and also no external proton source is needed. Therefore, they are fully self-sufficient, which affords them opportunities in biotechnological applications. One specific P450 peroxygenase, namely, P450 OleT_JE, reacts with long-chain linear fatty acids through oxidative decarboxylation to form hydrocarbons and, as such, has been implicated as a suitable source for the biosynthesis of biofuels. Unfortunately, the reactions were shown to produce a considerable amount of side products originating from C^α and C^β hydroxylation and desaturation. These product distributions were found to be strongly dependent on whether the substrate had substituents on the C^α and/or C^β atoms. To understand the bifurcation pathways of substrate activation by P450 OleT_JE leading to decarboxylation, C^α hydroxylation, C^β hydroxylation and C^α–C^β desaturation, we performed a computational study using 3-phenylpropionate and 2-phenylbutyrate as substrates. We set up large cluster models containing the heme, the substrate and the key features of the substrate binding pocket and calculated (using density functional theory) the pathways leading to the four possible products. This work predicts that the two substrates will react with different reaction rates due to accessibility differences of the substrates to the active oxidant, and, as a consequence, these two substrates will also generate different products. This work explains how the substrate binding pocket of P450 OleT_JE guides a reaction to a chemoselectivity.

Modeling the hydroxylation of estragole via human liver cytochrome P450

Natural compounds derived from plants are generally regarded safe and devoid of adverse effects. However, there are individual ingredients that possess toxic, genotoxic, and carcinogenic activities. These compounds when exposed at specific level become hazardous to health. Estragole (1-allyl-4-methoxybenzene) is a common component of spice plants. Its toxicity gets activated with the hydroxylation at benzylic carbon (C1') position by P450 enzymes present in the human liver. The present study grounds to explore the reaction mechanism of conversion of estragole to hydroxylated metabolite using computational methodology. Density functional theory (DFT)-based calculations were employed to explore the cytochrome P450-catalyzed mechanism at C1 position aliphatic hydroxylation of estragole. Overall reaction energy profile, electronic configuration, and 3D structure of all intermediates, transition states, and product complexes formed during the reaction along with their free energies were tried to be investigated.

Catalytic Mechanism of Aromatic Nitration by Cytochrome P450 TxtE: Involvement of a Ferric-Peroxynitrite Intermediate

The cytochromes P450 are heme-dependent enzymes that catalyze many vital reaction processes in the human body related to biodegradation and biosynthesis. They typically act as mono-oxygenases; however, the recently discovered P450 subfamily TxtE utilizes O₂ and NO to nitrate aromatic substrates such as L-tryptophan. A direct and selective aromatic nitration reaction may be useful in biotechnology for the synthesis of drugs or small molecules. Details of the catalytic mechanism are unknown, and it has been suggested that the reaction should proceed through either an iron(III)-superoxo or an iron(II)-nitrosyl intermediate. To resolve this controversy, we used stopped-flow kinetics to provide evidence for a catalytic cycle where dioxygen binds prior to NO to generate an active iron(III)-peroxynitrite species that is able to nitrate l-Trp efficiently. We show that the rate of binding of O₂ is faster than that of NO and also leads to l-Trp nitration, while little evidence of product formation is observed from the iron(II)-nitrosyl complex. To support the experimental studies, we performed density functional theory studies on large active site cluster models. The studies suggest a mechanism involving an iron(III)-peroxynitrite that splits homolytically to form an iron(IV)-oxo heme (Compound II) and a free NO₂ radical via a small free energy of activation. The latter activates the substrate on the aromatic ring, while compound II picks up the ipso-hydrogen to form the product. The calculations give small reaction barriers for most steps in the catalytic cycle and, therefore, predict fast product formation from the iron(III)-peroxynitrite complex. These findings provide the first detailed insight into the mechanism of nitration by a member of the TxtE subfamily and highlight how the enzyme facilitates this novel reaction chemistry.

Theoretical Identification of the Factors Governing the Reactivity of C-H Bond Activation by Non-Heme Iron(IV)-Oxo Complexes

Selective functionalization of C-H bonds provides a straightforward approach to a large variety of well-defined derivatives. High-valent mononuclear iron(IV)-oxo complexes are proposed to carry out these C-H activation reactions in enzymes or in biomimetic syntheses. In this Minireview, we aim to highlight the features that delineate the distinct reactivity of non-heme oxo-iron(IV) motifs to cleave strong C-H bonds in hydrocarbons, primarily focusing on the hydrogen atom transfer (HAT) process. We describe how the structural and electronic properties of supporting ligands modulate the oxidative property of the iron(IV)-oxo complexes. Furthermore, we highlight the decisive role played by spin-state in these biomimetic reactions. We also discuss how tunneling and external perturbations like electric field influence the transfer of hydrogen atoms. Lastly, we emphasize how computations could work as a practical guide to sketch and develop synthetic models with greater efficacy.

Epoxidation Catalyzed by the Nonheme Iron(II)- and 2-Oxoglutarate-Dependent Oxygenase, AsqJ: Mechanistic Elucidation of Oxygen Atom Transfer by a Ferryl Intermediate

Mechanisms of enzymatic epoxidation via oxygen atom transfer (OAT) to an olefin moiety is mainly derived from the studies on thiolate-heme containing epoxidases, such as cytochrome P450 epoxidases. The molecular basis of epoxidation catalyzed by nonheme-iron enzymes is much less explored. Herein, we present a detailed study on epoxidation catalyzed by the nonheme iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase, AsqJ. The native substrate and analogues with different para substituents ranging from electron-donating groups (e.g., methoxy) to electron-withdrawing groups (e.g., trifluoromethyl) were used to probe the mechanism. The results derived from transient-state enzyme kinetics, Mössbauer spectroscopy, reaction product analysis, X-ray crystallography, density functional theory calculations, and molecular dynamic simulations collectively revealed the following mechanistic insights: (1) The rapid O₂ addition to the AsqJ Fe(II) center occurs with the iron-bound 2OG adopting an online-binding mode in which the C1 carboxylate group of 2OG is trans to the proximal histidine (His134) of the 2-His-1-carboxylate facial triad, instead of assuming the offline-binding mode with the C1 carboxylate group trans to the distal histidine (His211); (2) The decay rate constant of the ferryl intermediate is not strongly affected by the nature of the para substituents of the substrate during the OAT step, a reactivity behavior that is drastically different from nonheme Fe(IV)-oxo synthetic model complexes; (3) The OAT step most likely proceeds through a stepwise process with the initial formation of a C(benzylic)-O bond to generate an Fe-alkoxide species, which is observed in the AsqJ crystal structure. The subsequent C3-O bond formation completes the epoxide installation.

Spin state and reactivity of iron(iv)oxido complexes with tetradentate bispidine ligands

The iron(iv)oxido complex [(bispidine)Fe^IVO(Cl)]⁺is shown by experiment and high-level DLPNO-CCSD(T) quantum-chemical calculations to be an extremely short-lived and very reactive intermediate-spin (S= 1) species.

Sluggish reactivity by a nonheme iron(iv)-tosylimido complex as compared to its oxo analogue

A comparative spectroscopic and computational study of reactivity between ferryl-tosylimido and ferryl-oxo complexes of two biomimetic model systems. The Fe(iv)-tosylimido complex was found to be sluggish in comparison to its fellow oxo counterpart.

Flavonol biosynthesis by nonheme iron dioxygenases: A computational study into the structure and mechanism

Plants produce flavonol compounds for vital functions regarding plant growth, fruit and flower colouring as well as fruit ripening processes. Several of these biosynthesis steps are stereo- and regioselective and are being carried out by nonheme iron enzymes. Using density functional theory calculations on a large active site model complex of flavanone-3β-hydroxylase (FHT), we established the mechanism for conversion of naringenin to its dihydroflavonol, which is a key step in the mechanism of flavonol biosynthesis. The reaction starts with dioxygen binding to the iron(II) centre and a reaction with α-ketoglutarate co-substrate gives succinate, an iron(IV)-oxo species and CO₂ with large exothermicity and small reaction barriers. The rate-determining reaction step in the mechanism; however, is hydrogen atom abstraction of an aliphatic CH bond by the iron(IV)-oxo species. We identify a large kinetic isotope effect for the replacement of the transferring hydrogen atom by deuterium. In a final step the OH and substrate radicals combine to form the alcohol product with a barrier of several kcal mol^-1. We show that the latter is the result of geometric constraints in the active site pocket. Furthermore, the calculations show that a weak tertiary CH bond is shielded from the iron(IV)-oxo species in the substrate binding position and therefore the enzyme is able to activate a stronger CH bond. As such, the flavanone-3β-hydroxylase enzyme reacts regioselectively with one specific CH bond of naringenin by avoiding activation of weaker bonds through tight substrate and oxidant positioning.

Mechanistic Studies of Fatty Acid Activation by CYP152 Peroxygenases Reveal Unexpected Desaturase Activity

The majority of cytochrome P450 enzymes (CYPs) predominantly operate as monooxygenases, but recently a class of P450 enzymes was discovered, that can act as peroxygenases (CYP152). These enzymes convert fatty acids through oxidative decarboxylation, yielding terminal alkenes, and through α- and β-hydroxylation to yield hydroxy-fatty acids. Bioderived olefins may serve as biofuels, and hence understanding the mechanism and substrate scope of this class of enzymes is important. In this work, we report on the substrate scope and catalytic promiscuity of CYP OleT_JE and two of its orthologues from the CYP152 family, utilizing α-monosubstituted branched carboxylic acids. We identify α,β-desaturation as an unexpected dominant pathway for CYP OleT_JE with 2-methylbutyric acid. To rationalize product distributions arising from α/β-hydroxylation, oxidative decarboxylation, and desaturation depending on the substrate’s structure and binding pattern, a computational study was performed based on an active site complex of CYP OleT_JE containing the heme cofactor in the substrate binding pocket and 2-methylbutyric acid as substrate. It is shown that substrate positioning determines the accessibility of the oxidizing species (Compound I) to the substrate and hence the regio- and chemoselectivity of the reaction. Furthermore, the results show that, for 2-methylbutyric acid, α,β-desaturation is favorable because of a rate-determining α-hydrogen atom abstraction, which cannot proceed to decarboxylation. Moreover, substrate hydroxylation is energetically impeded due to the tight shape and size of the substrate binding pocket.

Does Substrate Positioning Affect the Selectivity and Reactivity in the Hectochlorin Biosynthesis Halogenase?

In this work we present the first computational study on the hectochlorin biosynthesis enzyme HctB, which is a unique three-domain halogenase that activates non-amino acid moieties tethered to an acyl-carrier, and as such may have biotechnological relevance beyond other halogenases. We use a combination of small cluster models and full enzyme structures calculated with quantum mechanics/molecular mechanics methods. Our work reveals that the reaction is initiated with a rate-determining hydrogen atom abstraction from substrate by an iron (IV)-oxo species, which creates an iron (III)-hydroxo intermediate. In a subsequent step the reaction can bifurcate to either halogenation or hydroxylation of substrate, but substrate binding and positioning drives the reaction to optimal substrate halogenation. Furthermore, several key residues in the protein have been identified for their involvement in charge-dipole interactions and induced electric field effects. In particular, two charged second coordination sphere amino acid residues (Glu223 and Arg245) appear to influence the charge density on the Cl ligand and push the mechanism toward halogenation. Our studies, therefore, conclude that nonheme iron halogenases have a chemical structure that induces an electric field on the active site that affects the halide and iron charge distributions and enable efficient halogenation. As such, HctB is intricately designed for a substrate halogenation and operates distinctly different from other nonheme iron halogenases.

A Comparative Review on the Catalytic Mechanism of Nonheme Iron Hydroxylases and Halogenases

Enzymatic halogenation and haloperoxidation are unusual processes in biology; however, a range of halogenases and haloperoxidases exist that are able to transfer an aliphatic or aromatic C–H bond into C–Cl/C–Br. Haloperoxidases utilize hydrogen peroxide, and in a reaction with halides (Cl−/Br−), they react to form hypohalides (OCl−/OBr−) that subsequently react with substrate by halide transfer. There are three types of haloperoxidases, namely the iron-heme, nonheme vanadium, and flavin-dependent haloperoxidases that are reviewed here. In addition, there are the nonheme iron halogenases that show structural and functional similarity to the nonheme iron hydroxylases and form an iron(IV)-oxo active species from a reaction of molecular oxygen with α-ketoglutarate on an iron(II) center. They subsequently transfer a halide (Cl−/Br−) to an aliphatic C–H bond. We review the mechanism and function of nonheme iron halogenases and hydroxylases and show recent computational modelling studies of our group on the hectochlorin biosynthesis enzyme and prolyl-4-hydroxylase as examples of nonheme iron halogenases and hydroxylases. These studies have established the catalytic mechanism of these enzymes and show the importance of substrate and oxidant positioning on the stereo-, chemo- and regioselectivity of the reaction that takes place.

High-Oxidation-State 3d Metal (Ti-Cu) Complexes with N-Heterocyclic Carbene Ligation

High-oxidation-state 3d metal species have found a wide range of applications in modern synthetic chemistry and materials science. They are also implicated as key reactive species in biological reactions. These applications have thus prompted explorations of their formation, structure, and properties. While the traditional wisdom regarding these species was gained mainly from complexes supported by nitrogen- and oxygen-donor ligands, recent studies with N-heterocyclic carbenes (NHCs), which are widely used for the preparation of low-oxidation-state transition metal complexes in organometallic chemistry, have led to the preparation of a large variety of isolable high-oxidation-state 3d metal complexes with NHC ligation. Since the first report in this area in the 1990s, isolable complexes of this type have been reported for titanium(IV), vanadium(IV,V), chromium(IV,V), manganese(IV,V), iron(III,IV,V), cobalt(III,IV,V), nickel(IV), and copper(II). With the aim of providing an overview of this intriguing field, this Review summarizes our current understanding of the synthetic methods, structure and spectroscopic features, reactivity, and catalytic applications of high-oxidation-state 3d metal NHC complexes of titanium to copper. In addition to this progress, factors affecting the stability and reactivity of high-oxidation-state 3d metal NHC species are also presented, as well as perspectives on future efforts.

Quantum Mechanics/Molecular Mechanics Studies on the Relative Reactivities of Compound I and II in Cytochrome P450 Enzymes

The cytochromes P450 are drug metabolizing enzymes in the body that typically react with substrates through a monoxygenation reaction. During the catalytic cycle two reduction and protonation steps generate a high-valent iron (IV)-oxo heme cation radical species called Compound I. However, with sufficient reduction equivalents present, the catalytic cycle should be able to continue to the reduced species of Compound I, called Compound II, rather than a reaction of Compound I with substrate. In particular, since electron transfer is usually on faster timescales than atom transfer, we considered this process feasible and decided to investigate the reaction computationally. In this work we present a computational study using density functional theory methods on active site model complexes alongside quantum mechanics/molecular mechanics calculations on full enzyme structures of cytochrome P450 enzymes. Specifically, we focus on the relative reactivity of Compound I and II with a model substrate for O⁻H bond activation. We show that generally the barrier heights for hydrogen atom abstraction are higher in energy for Compound II than Compound I for O⁻H bond activation. Nevertheless, for the activation of such bonds, Compound II should still be an active oxidant under enzymatic conditions. As such, our computational modelling predicts that under high-reduction environments the cytochromes P450 can react with substrates via Compound II but the rates will be much slower.

Mechanistic Insight on the Activity and Substrate Selectivity of Nonheme Iron Dioxygenases

Nonheme iron dioxygenases catalyze vital reactions for human health particularly related to aging processes. They are involved in the biosynthesis of amino acids, but also the biodegradation of toxic compounds. Typically they react with their substrate(s) through oxygen atom transfer, although often with the assistance of a co-substrate like α-ketoglutarate that is converted to succinate and CO₂ . Many reaction processes catalyzed by the nonheme iron dioxygenases are stereoselective or regiospecific and hence understanding the mechanism and protein involvement in the selectivity is important for the design of biotechnological applications of these enzymes. To this end, I will review recent work of our group on nonheme iron dioxygenases and include background information on their general structure and catalytic cycle. Examples of stereoselective and regiospecific reaction mechanisms we elucidated are for the AlkB repair enzyme, prolyl-4-hydroxylase and the ergothioneine biosynthesis enzyme. Finally, I cover an example where we bioengineered S-p-hydroxymandelate synthase into the R-p-hydroxymandelate synthase.

Dramatic rate-enhancement of oxygen atom transfer by an iron(iv)-oxo species by equatorial ligand field perturbations

The reactivity and characterization of a novel iron(iv)-oxo species is reported that gives enhanced reactivity as a result of second-coordination sphere perturbations of the ligand system.

Nitrogen Reduction to Ammonia on a Biomimetic Mononuclear Iron Centre: Insights into the Nitrogenase Enzyme

Nitrogenases catalyse nitrogen fixation to ammonia on a multinuclear Fe‐Mo centre, but their mechanism and particularly the order of proton and electron transfer processes that happen during the catalytic cycle is still unresolved. Recently, a unique biomimetic mononuclear iron model was developed using tris(phosphine)borate (TPB) ligands that was shown to convert N₂ into NH₃. Herein, we present a computational study on the [(TPB)FeN₂]⁻ complex and describe its conversion into ammonia through the addition of electrons and protons. In particular, we tested the consecutive proton transfer on only the distal nitrogen atom or alternated protonation of the distal/proximal nitrogen. It is found that the lowest energy pathway is consecutive addition of three protons to the same site, which forms ammonia and an iron‐nitrido complex. In addition, the proton transfer step of complexes with the metal in various oxidation and spin states were tested and show that the pK_a values of biomimetic mononuclear nitrogenase intermediates vary little with iron oxidation states. As such, the model gives several possible NH₃ formation pathways depending on the order of electron/proton transfer, and all should be physically accessible in the natural system. These results may have implications for enzymatic nitrogenases and give insight into the catalytic properties of mononuclear iron centres.

Keto-Enol Tautomerization Triggers an Electrophilic Aldehyde Deformylation Reaction by a Nonheme Manganese(III)-Peroxo Complex

Oxygen atom transfer by high-valent enzymatic intermediates remains an enigma in chemical catalysis. In particular, manganese is an important first-row metal involved in key biochemical processes, including the biosynthesis of molecular oxygen (through the photosystem II complex) and biodegradation of toxic superoxide to hydrogen peroxide by superoxide dismutase. Biomimetic models of these biological systems have been developed to gain understanding on the structure and properties of short-lived intermediates but also with the aim to create environmentally benign oxidants. In this work, we report a combined spectroscopy, kinetics and computational study on aldehyde deformylation by two side-on manganese(III)-peroxo complexes with bispidine ligands. Both manganese(III)-peroxo complexes are characterized by UV-vis and mass spectrometry techniques, and their reactivity patterns with aldehydes was investigated. We find a novel mechanism for the reaction that is initiated by a hydrogen atom abstraction reaction, which enables a keto-enol tautomerization in the substrate. This is an essential step in the mechanism that makes an electrophilic attack on the olefin bond possible as the attack on the aldehyde carbonyl is too high in energy. Kinetics studies determine a large kinetic isotope effect for the replacement of the transferring hydrogen atom by deuterium, while replacing the transferring hydrogen atom by a methyl group makes the substrate inactive and hence confirm the hypothesized mechanism. Our new mechanism is confirmed with density functional theory modeling on the full mechanism and rationalized through valence bond and thermochemical cycles. Our unprecedented new mechanism may have relevance to biological and biomimetic chemistry processes in general and gives insight into the reactivity patterns of metal-peroxo and metal-hydroperoxo intermediates in general.

How Are Substrate Binding and Catalysis Affected by Mutating Glu127 and Arg161 in Prolyl-4-hydroxylase? A QM/MM and MD Study

Prolyl-4-hydroxylase is a vital enzyme for human physiology involved in the biosynthesis of 4-hydroxyproline, an essential component for collagen formation. The enzyme performs a unique stereo- and regioselective hydroxylation at the C4 position of proline despite the fact that the C5 hydrogen atoms should be thermodynamically easier to abstract. To gain insight into the mechanism and find the origin of this regioselectivity, we have done a quantum mechanics/molecular mechanics (QM/MM) study on wildtype and mutant structures. In a previous study (Timmins et al., 2017) we identified several active site residues critical for substrate binding and positioning. In particular, the Glu127 and Arg161 were shown to form multiple hydrogen bonding and ion-dipole interactions with substrate and could thereby affect the regio- and stereoselectivity of the reaction. In this work, we decided to test that hypothesis and report a QM/MM and molecular dynamics (MD) study on prolyl-4-hydroxylase and several active site mutants where Glu127 or Arg161 are mutated for Asp, Gln, or Lys. Thus, the R161D and R161Q mutants give very high barriers for hydrogen atom abstraction from any proline C-H bond and therefore will be inactive. The R161K mutant, by contrast, sees the regio- and stereoselectivity of the reaction change but still is expected to hydroxylate proline at room temperature. By contrast, the Glu127 mutants E127D and E127Q show possible changes in regioselectivity with the former being more probable to react compared to the latter.

Understanding How Prolyl-4-hydroxylase Structure Steers a Ferryl Oxidant toward Scission of a Strong C-H Bond

Prolyl-4-hydroxylase (P4H) is a non-heme iron hydroxylase that regio- and stereospecifically hydroxylates proline residues in a peptide chain into R-4-hydroxyproline, which is essential for collagen cross-linking purposes in the human body. Surprisingly, in P4H, a strong aliphatic C-H bond is activated, while thermodynamically much weaker aliphatic C-H groups, that is, at the C³ and C⁵ positions, are untouched. Little is known on the origins of the high regio- and stereoselectivity of P4H and many non-heme and heme enzymes in general, and insight into this matter may be relevant to Biotechnology as well as Drug Development. The active site of the protein contains two aromatic residues (Tyr₁₄₀ and Trp₂₄₃) that we expected to be crucial for guiding the regioselectivity of the reaction. We performed a detailed quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) study on wild-type and mutant structures. The work shows that Trp₂₄₃ is involved in key protein loop-loop interactions that affect the shape and size of the substrate binding pocket and its mutation has major long-range effects. By contrast, the Tyr₁₄₀ residue is shown to guide the regio- and stereoselectivity by holding the substrate and ferryl oxidant in a specific orientation through hydrogen bonding and π-stacking interactions. Compelling evidence is found that the Tyr₁₄₀ residue is involved in expelling the product from the binding pocket after the reaction is complete. It is shown that mutations where the hydrogen bonding network that involves the Tyr₁₄₀ and Trp₂₄₃ residues is disrupted lead to major changes in folding of the protein and the size and shape of the substrate binding pocket. Specifically, the Trp₂₄₃ residue positions the amino acid side chains of Arg₁₆₁ and Glu₁₂₇ in specific orientations with substrate. As such, the P4H enzyme is a carefully designed protein with a subtle and rigid secondary structure that enables the binding of substrate, guides the regioselectivity, and expels product efficiently.

Nonclassical Single-State Reactivity of an Oxo-Iron(IV) Complex Confined to Triplet Pathways

C-H bond activation mediated by oxo-iron (IV) species represents the key step of many heme and nonheme O₂-activating enzymes. Of crucial interest is the effect of spin state of the Fe^IV(O) unit. Here we report the C-H activation kinetics and corresponding theoretical investigations of an exclusive tetracarbene ligated oxo-iron(IV) complex, [L^NHCFe^IV(O)(MeCN)]²⁺ (1). Kinetic traces using substrates with bond dissociation energies (BDEs) up to 80 kcal mol^-1 show pseudo-first-order behavior and large but temperature-dependent kinetic isotope effects (KIE 32 at -40 °C). When compared with a topologically related oxo-iron(IV) complex bearing an equatorial N-donor ligand, [L^TMCFe^IV(O) (MeCN)]²⁺ (A), the tetracarbene complex 1 is significantly more reactive with second order rate constants k'₂ that are 2-3 orders of magnitude higher. UV-vis experiments in tandem with cryospray mass spectrometry evidence that the reaction occurs via formation of a hydroxo-iron(III) complex (4) after the initial H atom transfer (HAT). An extensive computational study using a wave function based multireference approach, viz. complete active space self-consistent field (CASSCF) followed by N-electron valence perturbation theory up to second order (NEVPT2), provided insight into the HAT trajectories of 1 and A. Calculated free energy barriers for 1 reasonably agree with experimental values. Because the strongly donating equatorial tetracarbene pushes the Fe-d_x²-y² orbital above d_z², 1 features a dramatically large quintet-triplet gap of ∼18 kcal/mol compared to ∼2-3 kcal/mol computed for A. Consequently, the HAT process performed by 1 occurs on the triplet surface only, in contrast to complex A reported to feature two-state-reactivity with contributions from both triplet and quintet states. Despite this, the reactive Fe^IV(O) units in 1 and A undergo the same electronic-structure changes during HAT. Thus, the unique complex 1 represents a pure "triplet-only" ferryl model.

Reactivity Patterns of (Protonated) Compound II and Compound I of Cytochrome P450: Which is the Better Oxidant?

The cytochromes P450 are versatile enzymes in human physiology that perform substrate hydroxylation reactions extremely efficiently. In this work, we present results of a computational study on the reactivity patterns of Compound I, Compound II, and protonated Compound II with model substrates, and we address the question of which of these compounds is the most effective oxidant? All calculations, regardless of the substrate, implicated that Compound I is the superior oxidant of the three. However, Compound II and protonated Compound II were found to react with free energies of activation that are only a few kcal mol^-1 higher in energy than those obtained with Compound I. Therefore, Compound II and protonated Compound II should be able to react with aliphatic groups with moderate C-H bond strengths. We have analysed all results in detail and have given electronic, thermochemical, valence bond, and molecular orbital rationalizations on the reactivity differences and explained experimental product distributions. Overall, the findings implied that alternative oxidants could operate alongside Compound I in complex reaction mechanisms of enzymatic and synthetic iron porphyrinoid complexes.