The Role of Equatorial and Axial Ligands in Promoting the Activity of Non-Heme Oxidoiron(IV) Catalysts in Alkane Hydroxylation

New Frontiers in Nonheme Enzymatic Oxyferryl Species

Non-heme mononuclear iron dependent (NHM-Fe) enzymes exhibit exceedingly diverse catalytic reactivities. Despite their catalytic versatilities, the mononuclear iron centers in these enzymes show a relatively simple architecture, in which an iron atom is ligated with 2-4 amino acid residues, including histidine, aspartic or glutamic acid. In the past two decades, a common high-valent reactive iron intermediate, the S=2 oxyferryl (Fe(IV)-oxo or Fe(IV)=O) species, has been repeatedly discovered in NHM-Fe enzymes containing a 2-His-Fe or 2-His-1-carboxylate-Fe center. However, for 3-His/4-His-Fe enzymes, no common reactive intermediate has been identified. Recently, we have spectroscopically characterized the first S=1 Fe(IV) intermediate in a 3-His-Fe containing enzyme, OvoA, which catalyzes a novel oxidative carbon-sulfur bond formation. In this review, we summarize the broad reactivities demonstrated by S=2 Fe(IV)-oxo intermediates, the discovery of the first S=1 Fe(IV) intermediate in OvoA and the mechanistic implication of such a discovery, and the intrinsic reactivity differences of the S=2 and the S=1 Fe(IV)-oxo species. Finally, we postulate the possible reasons to utilize an S=1 Fe(IV) species in OvoA and their implications to other 3-His/4-His-Fe enzymes.

A tale of two topological isomers: Uptuning [FeIV(O)(Me4cyclam)]2+ for olefin epoxidation

TMC-anti and TMC-syn, the two topological isomers of [FeIV(O)(TMC)(CH3CN)]2+ (TMC = 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane, or Me4cyclam), differ in the orientations of their FeIV=O units relative to the four methyl groups of the TMC ligand framework. The FeIV=O unit of TMC-anti points away from the four methyl groups, while that of TMC-syn is surrounded by the methyl groups, resulting in differences in their oxidative reactivities. TMC-syn reacts with HAT (hydrogen atom transfer) substrates at 1.3- to 3-fold faster rates than TMC-anti, but the reactivity difference increases dramatically in oxygen-atom transfer reactions. R2S substrates are oxidized into R2S=O products at rates 2-to-3 orders of magnitude faster by TMC-syn than TMC-anti. Even more remarkably, TMC-syn epoxidizes all the olefin substrates in this study, while TMC-anti reacts only with cis-cyclooctene but at a 100-fold slower rate. Comprehensive quantum chemical calculations have uncovered the key factors governing such reactivity differences found between these two topological isomers.

Aqueous Iron(IV)-Oxo Complex: An Emerging Powerful Reactive Oxidant Formed by Iron(II)-Based Advanced Oxidation Processes for Oxidative Water Treatment

High-valent iron(IV)-oxo complexes are of great significance as reactive intermediates implicated in diverse chemical and biological systems. The aqueous iron(IV)-oxo complex (Fe_aq^IVO²⁺) is the simplest but one of the most powerful ferryl ion species, which possesses a high-spin state, high reduction potential, and long lifetime. It has been well documented that Fe_aq^IVO²⁺ reacts with organic compounds through various pathways (hydrogen-atom, hydride, oxygen-atom, and electron transfer as well as electrophilic addition) at moderate reaction rates and show selective reactivity toward inorganic ions prevailing in natural water, which single out Fe_aq^IVO²⁺ as a superior candidate for oxidative water treatment. This review provides state-of-the-art knowledge on the chemical properties and oxidation mechanism and kinetics of Fe_aq^IVO²⁺, with special attention to the similarities and differences to two representative free radicals (hydroxyl radical and sulfate radical). Moreover, the prospective role of Fe_aq^IVO²⁺ in Fe_aq²⁺ activation-initiated advanced oxidation processes (AOPs) has been intensively investigated over the past 20 years, which has significantly challenged the conventional recognition that free radicals dominated in these AOPs. The latest progress in identifying the contribution of Fe_aq^IVO²⁺ in Fe_aq²⁺-based AOPs is thereby reviewed, highlighting controversies on the nature of the reactive oxidants formed in several Fe_aq²⁺ activated peroxide and oxyacid processes. Finally, future perspectives for advancing the evaluation of Fe_aq^IVO²⁺ reactivity from an engineering viewpoint are proposed.

Catalytic properties of the ferryl ion in the solid state: a computational review

This review summarises the last findings in the emerging field of heterogeneous catalytic oxidation of light alkanes by ferryl species supported on solid-state systems such as the conversion of methane into methanol by FeO-MOF74.

A Pseudotetrahedral Terminal Oxoiron(IV) Complex: Mechanistic Promiscuity in C-H Bond Oxidation Reactions

S=2 oxoiron(IV) species act as reactive intermediates in the catalytic cycle of nonheme iron oxygenases. The few available synthetic S=2 FeIV =O complexes known to date are often limited to trigonal bipyramidal and very rarely to octahedral geometries. Herein we describe the generation and characterization of an S=2 pseudotetrahedral FeIV =O complex 2 supported by the sterically demanding 1,4,7-tri-tert-butyl-1,4,7-triazacyclononane ligand. Complex 2 is a very potent oxidant in hydrogen atom abstraction (HAA) reactions with large non-classical deuterium kinetic isotope effects, suggesting hydrogen tunneling contributions. For sterically encumbered substrates, direct HAA is impeded and an alternative oxidative asynchronous proton-coupled electron transfer mechanism prevails, which is unique within the nonheme oxoiron community. The high reactivity and the similar spectroscopic parameters make 2 one of the best electronic and functional models for a biological oxoiron(IV) intermediate of taurine dioxygenase (TauD-J).

Unveiling the catalytic potential of the Fe(iv)oxo species for the oxidation of hydrocarbons in the solid state

We have investigated the three steps in the conversion of methane into methanol by Fe(iv)Ooxo species supported in MOF-74. We use ab initio MD and static approximations to predict the reaction barriers using enthalpy ΔH and free energy ΔG.

Deciphering the origin of million-fold reactivity observed for the open core diiron [HO-FeIII-O-FeIV[double bond, length as m-dash]O]2+ species towards C-H bond activation: role of spin-states, spin-coupling, and spin-cooperation

High-valent metal-oxo species have been characterised as key intermediates in both heme and non-heme enzymes that are found to perform efficient aliphatic hydroxylation, epoxidation, halogenation, and dehydrogenation reactions. Several biomimetic model complexes have been synthesised over the years to mimic both the structure and function of metalloenzymes. The diamond-core [Fe2(μ-O)2] is one of the celebrated models in this context as this has been proposed as the catalytically active species in soluble methane monooxygenase enzymes (sMMO), which perform the challenging chemical conversion of methane to methanol at ease. In this context, a report of open core [HO(L)FeIII-O-FeIV(O)(L)]2+ (1) gains attention as this activates C-H bonds a million-fold faster compared to the diamond-core structure and has the dual catalytic ability to perform hydroxylation as well as desaturation with organic substrates. In this study, we have employed density functional methods to probe the origin of the very high reactivity observed for this complex and also to shed light on how this complex performs efficient hydroxylation and desaturation of alkanes. By modelling fifteen possible spin-states for 1 that could potentially participate in the reaction mechanism, our calculations reveal a doublet ground state for 1 arising from antiferromagnetic coupling between the quartet FeIV centre and the sextet FeIII centre, which regulates the reactivity of this species. The unusual stabilisation of the high-spin ground state for FeIV[double bond, length as m-dash]O is due to the strong overlap of with the orbital, reducing the antibonding interactions via spin-cooperation. The electronic structure features computed for 1 are consistent with experiments offering confidence in the methodology chosen. Further, we have probed various mechanistic pathways for the C-H bond activation as well as -OH rebound/desaturation of alkanes. An extremely small barrier height computed for the first hydrogen atom abstraction by the terminal FeIV[double bond, length as m-dash]O unit was found to be responsible for the million-fold activation observed in the experiments. The barrier height computed for -OH rebound by the FeIII-OH unit is also smaller suggesting a facile hydroxylation of organic substrates by 1. A strong spin-cooperation between the two iron centres also reduces the barrier for second hydrogen atom abstraction, thus making the desaturation pathway competitive. Both the spin-state as well as spin-coupling between the two metal centres play a crucial role in dictating the reactivity for species 1. By exploring various mechanistic pathways, our study unveils the fact that the bridged μ-oxo group is a poor electrophile for both C-H activation as well for -OH rebound. As more and more evidence is gathered in recent years for the open core geometry of sMMO enzymes, the idea of enhancing the reactivity via an open-core motif has far-reaching consequences.

Surface-orientation- and ligand-dependent quenching of the spin magnetic moment of Co porphyrins adsorbed on Cu substrates

Porphyrin molecules are particularly interesting candidates for spintronic applications due to their bonding flexibility, which allows to modify their properties substantially by the addition or transformation of ligands. Here, we investigate the electronic and magnetic properties of cobalt octaethylporphyrin (CoOEP), deposited on copper substrates with two distinct crystallographic surface orientations, Cu(100) and Cu(111), with X-ray absorption spectroscopy (XAS) and X-ray magnetic circular dichroism (XMCD). A significant magnetic moment is present in the Co ions of the molecules deposited on Cu(100), but it is completely quenched on Cu(111). Heating the molecules on both substrates to 500 K induces a ring-closure reaction with cobalt tetrabenzoporphyrin (CoTBP) as reaction product. In these molecules, the magnetic moment is quenched on both surfaces. Our XMCD and XAS measurements suggest that the filling of the d_z² orbital leads to a non-integer valence state and causes the quench of the spin moments on all samples except CoOEP/Cu(100), where the molecular conformation induces variations to the ligand field that lift the quench. We further employ density functional theory calculations, supplemented with on-site Coulomb correlations (DFT+U), to study the adsorption of these spin-bearing molecules on the Cu substrates. Our calculations show that charge transfer from the Cu substrates as well as charge redistribution within the Co 3d orbitals lead to the filling of the Co minority spin d_z² orbital, causing a 'turning off' of the exchange splitting and quenching of the spin moment at the Co magnetic centers. Our investigations suggest that, by this mechanism, molecule-substrate interactions can be used to control the quenching of the magnetic moments of the adsorbed molecules.

Density-functional theory models of Fe(iv)O reactivity in metal-organic frameworks: self-interaction error, spin delocalisation and the role of hybrid exchange

We study the reactivity of Fe(iv)O moieties supported by a metal-organic framework (MOF-74) in the oxidation reaction of methane to methanol using all-electron, periodic density-functional theory calculations. We compare results concerning the electronic properties and reactivity obtained using two hybrid (B3LYP and sc-BLYP) and two standard generalised gradient corrected (PBE and BLYP) semi-local density functional approximations. The semi-local functionals are unable to reproduce the expected reaction profiles and yield a qualitatively incorrect representation of the reactivity. Non-local hybrid functionals provide a substantially more reliable description and predict relatively modest (ca. 60 kJ mol^-1) reaction energy barriers for the H-atom abstraction reaction from CH₄ molecules. We examine the origin of these differences and we highlight potential means to overcome the limitations of standard semi-local functionals in reactivity calculations in solid-state systems.

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.

Mechanism and stereoselectivity of benzylic C-H hydroxylation by Ru-porphyrin: a computational study

The mechanism and origin of the stereoselectivity of asymmetric benzylic C-H hydroxylation by Ru-porphyrin were elucidated with density functional theory calculations. The reaction proceeds via a hydrogen-atom abstraction/oxygen-rebound pathway, wherein a high-valent ruthenium-oxo species abstracts a hydrogen atom from ethylbenzene to generate a radical pair intermediate, followed by the oxygen-rebound process to form 1-phenylethanol. The hydrogen-atom abstraction step is the rate- and stereoselectivity-determining step. Based on the mechanistic model, the computed stereoselectivity is in agreement with the experimental observations. Analysis of the distortion/interaction model suggests that stereoselectivity is determined by both the distortion energy of the ethylbenzene and the interaction energy between the ethylbenzene and the chiral Ru-porphyrin. The steric repulsion between the phenyl group of ethylbenzene and the bulky substituent of Ru-porphyrin is the leading cause of chiral induction.

Advancing Fenton and photo-Fenton water treatment through the catalyst design

The review is devoted to modern Fenton, photo-Fenton, as well as Fenton-like and photo-Fenton-like reactions with participation of iron species in liquid phase and as heterogeneous catalysts. Mechanisms of these reactions were considered that include hydroxyl radical and oxoferryl species as the reactive intermediates. The barriers in the way of application of these reactions to wastewater treatment were discussed. The following fundamental problems need further research efforts: inclusion of more mechanism steps and quantum calculations of all rate constants lacking in the literature, checking the outer sphere electron transfer contribution, determination of the causes for the key changes in the homogeneous Fenton reaction mechanism with a change in the reagents concentration. The key advances for Fenton reactions implementation for the water treatment are related to tremendous hydrodynamical effects on the catalytic activity, design of ligands for high rate and completeness of mineralization in short time, and design of highly active heterogeneous catalysts. While both homogeneous and heterogeneous Fenton and photo-Fenton systems are open for further improvements, heterogeneous photo-Fenton systems are most promising for practical applications because of the inherent higher catalyst stability. Modern methods of quantum chemistry are expected to play a continuously increasing role in development of such catalysts.

Electronic structure and reactivity of Fe(iv)oxo species in metal-organic frameworks

We investigate the potential use of Fe(iv)oxo species supported on a metal-organic framework in the catalytic hydroxylation of methane to produce methanol. We use periodic density-functional theory calculations at the 6-31G**/B3LYP level of theory to study the electronic structure and chemical reactivity in the hydrogen abstraction reaction from methane in the presence of Fe(iv)O(oxo) supported on MOF-74. Our results indicate that the Fe(iv)O moiety in MOF-74 is characterised by a highly reactive (quintet) ground-state, with a distance between Fe(iv) and O(oxo) of 1.601 Å, consistent with other high-spin Fe(iv)O inorganic complexes in the gas phase and in aqueous solution. Similar to the latter systems, the highly electrophilic character (and thus the reactivity) of Fe(iv)O in MOF-74 is determined by the presence of a low-lying anti-bonding virtual orbital (3σ*), which acts as an electron acceptor in the early stages of the hydrogen atom abstraction from methane. We estimate an energy barrier for hydrogen abstraction of 50.77 kJ mol-1, which is comparable to the values estimated in other gas-phase and hydrated Fe(iv)O-based complexes with the ability to oxidise methane. Our findings therefore suggest that metal-organic frameworks can provide suitable supports to develop new solid-state catalysts for organic oxidation reactions.

The enzyme-like catalytic hydrogen abstraction reaction mechanisms of cyclic hydrocarbons with magnesium-diluted Fe-MOF-74

Enzymatic heme and non-heme Fe(iv)–O species usually play an important role in hydrogen abstraction of biocatalytic reactions, yet duplicating the reactivity in biomimicry remains a great challenge. Based on Xiao et al.'s experimental work [Nat. Chem., 2014, 6(7), 590], we theoretically found that in the presence of the oxidant N₂O, the enzyme-like metal organic framework, i.e., magnesium-diluted Fe-MOF-74 [Fe/(Mg)-MOF-74] can activate the C–H bonds of 1,4-cyclohexadiene (CHD) into benzene with a two-step hydrogen abstraction mechanism based on the density functional theory (DFT) level. It is shown that the first transition state about the cleavage of the N–O bond of N₂O to form the Fe(iv)–O species is the rate-determining step with activation enthalpy of 19.4 kcal mol⁻¹ and the complete reaction is exothermic by 62.8 kcal mol⁻¹ on quintet rather than on triplet PES. In addition, we proposed a rebound mechanism of cyclic cyclohexane (CHA) hydroxylation to cyclohexanol which has not been studied experimentally. Note that the activation enthalpies on the first hydrogen abstraction for both cyclic CHD and cyclohexane are just 8.1 and 3.5 kcal mol⁻¹, respectively, which are less than that of 13.9 kcal mol⁻¹ for chained ethane. Most importantly, for the hydrogen abstraction of methane catalyzed by M/(Mg)-MOF-74 (M = Cu, Ni, Fe, and Co), we found that the activation enthalpies versus the C–H bond length of methane of TSs, NPA charge of the reacting oxyl atom have linear relationships with different slopes, i.e., shorter C–H bond and less absolute value of NPA charge of oxyl atom are associated with lower activation enthalpy; while for the activation of methane, ethane, propane and CHD catalyzed by Fe/(Mg)-MOF-74, there also exists positive correlations between activation enthalpies, bond dissociation energies (BDEs) and C–H bond lengths in TSs, respectively. We hope the present theoretical study may provide the guideline to predict the performance of MOFs in C–H bond activation reactions.

The enzyme-like catalytic hydrogen abstraction reaction of cyclic hydrocarbons.

Electronic Structure Contributions of Non-Heme Oxo-Iron(V) Complexes to the Reactivity

Oxo-iron(V) species have been implicated in the catalytic cycle of the Rieske dioxygenase. Their synthetic analog, [Fe^V(O)(OC(O)CH₃)(PyNMe₃)]²⁺ (1, PyNMe₃ = 3,6,9,15-tetraazabicyclo[9.3.1]pentadeca-1(15),11,13-triene-3,6,9-trimethyl), derived from the O-O bond cleavage of its acetylperoxo iron(III) precursor, has been shown experimentally to perform regio- and stereoselective C-H and C═C bond functionalization. However, its structure-activity relation is poorly understood. Herein we present a detailed electronic-structure and spectroscopic analysis of complex 1 along with well-characterized oxo-iron(V) complexes, [Fe^V(O)(TAML)]^- (2, TAML = tetraamido macrocyclic ligand), [Fe^V(O)(TMC)(NC(O)CH₃)]⁺ (4, TMC = 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane), and [Fe^V(O)(TMC)(NC(OH)CH₃)]²⁺ (4-H⁺), using wave function-based multireference complete active-space self-consistent field calculations. Our results reveal that the x/ y anisotropy of the ⁵⁷Fe A-matrix is not a reliable spectroscopic marker to identify oxo-iron(V) species and that the drastically different A _x and A _y values determined for complexes 1, 4, and 4-H⁺ have distinctive origins compared to complex 2, a genuine oxo-iron(V) species. Complex 1, in fact, has a dominant character of [Fe^IV(O···OC(O)CH₃)^2-•]²⁺, i.e., an S_Fe = 1 iron(IV) center antiferromagnetically coupled to an O-O σ* radical, where the O-O bond has not been completely broken. Complex 4 is best described as a triplet ferryl unit that strongly interacts with the trans acetylimidyl radical in an antiferromagnetic fashion, [Fe^IV(O)(^•N═C(O^-)CH₃)]⁺. Complex 4-H⁺ features a similar electronic structure, [Fe^IV(O)(^•N═C(OH)CH₃)]²⁺. Owing to the remaining approximate half σ-bond in the O-O moiety, complex 1 can arrange two electron-accepting orbitals (α σ*_O-O and β Fe-d _xz) in such a way that both orbitals can simultaneously interact with the doubly occupied electron-donating orbitals (σ_C-H or π_C-C). Hence, complex 1 can promote a concerted yet asynchronous two-electron oxidation of the C-H and C═C bonds, which nicely explains the stereospecificity observed for complex 1 and the related species.

What drives the H-abstraction reaction in bio-mimetic oxoiron-bTAML complexes? A computational investigation

Monomeric iron-oxo units have been confirmed as intermediates involved in the C-H bond activation in various metallo-enzymes. Biomimetic oxoiron complexes of the biuret modified tetra-amido macrocyclic ligand (bTAML) have been demonstrated to oxidize a wide variety of unactivated C-H bonds. In the current work, density functional theory (DFT) has been employed to investigate the hydrogen abstraction (HAT) reactivity differences across a series of bTAML complexes. The cause for the differences in the HAT energy barriers has been found to be the relative changes in the energy of the frontier molecular orbitals (FMOs) induced by electronic perturbation.

A Highly Reactive Oxoiron(IV) Complex Supported by a Bioinspired N3 O Macrocyclic Ligand

The sluggish oxidants [FeIV (O)(TMC)(CH3 CN)]2+ (TMC=1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane) and [FeIV (O)(TMCN-d12 )(OTf)]+ (TMCN-d12 =1,4,7,11-tetra(methyl-d3 )-1,4,7,11-tetraazacyclotetradecane) are transformed into the highly reactive oxidant [FeIV (O)(TMCO)(OTf)]+ (1; TMCO=4,8,12-trimethyl-1-oxa-4,8,12-triazacyclotetradecane) upon replacement of an NMe donor in the TMC and TMCN ligands by an O atom. A rate enhancement of five to six orders of magnitude in both H atom and O atom transfer reactions was observed upon oxygen incorporation into the macrocyclic ligand. This finding was explained in terms of the higher electrophilicity of the iron center and the higher availability of the more reactive S=2 state in 1. This rationalizes nature's preference for using O-rich ligand environments for the hydroxylation of strong C-H bonds in enzymatic reactions.

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.

Spinning around in Transition-Metal Chemistry

The great diversity and richness of transition metal chemistry, such as the features of an open d-shell, opened a way to numerous areas of scientific research and technological applications. Depending on the nature of the metal and its environment, there are often several energetically accessible spin states, and the progress in accurate theoretical treatment of this complicated phenomenon is presented in this Account. The spin state energetics of a transition metal complex can be predicted theoretically on the basis of density functional theory (DFT) or wave function based methodology, where DFT has advantages since it can be applied routinely to medium-to-large-sized molecules and spin-state consistent density functionals are now available. Additional factors such as the effect of the basis set, thermochemical contributions, solvation, relativity, and dispersion, have been investigated by many researchers, but challenges in unambiguous assignment of spin states still remain. The first DFT studies showed intrinsic spin-state preferences of hybrid functionals for high spin and early generalized gradient approximation functionals for low spin. Progress in the development of density functional approximations (DFAs) then led to a class of specially designed DFAs, such as OPBE, SSB-D, and S12g, and brought a very intriguing and fascinating observation that the spin states of transition metals and the S_N2 barriers of organic molecules are somehow intimately linked. Among the many noteworthy results that emerged from the search for the appropriate description of the complicated spin state preferences in transition metals, we mainly focused on the examination of the connection between the spin state and the structures or coordination modes of the transition metal complexes. Changes in spin states normally lead only to changes in the metal-ligand bond lengths, but to the best of our knowledge, the dapsox ligand showed the first example of a transition-metal complex where a change in spin state leads also to changes in the coordination, switching between pentagonal-bipyramidal and capped-octahedron. Moreover, we have summarized the results of the thorough study that corrected the experimental assignment of the nature of the recently synthesized Sc³⁺ adduct of [Fe^IV(O)(TMC)]²⁺ (TMC = 1,4,8,11-tetramethylcyclam) and firmly established that the Sc³⁺-capped iron-oxygen complex corresponds to high-spin Fe^III. Last, but not least, we have provided deeper insight and rationalization of the observation that unlike in metalloenzymes, where the Fe^IV-oxo is usually observed with high spin, biomimetic Fe^IV-oxo complexes typically have a intermediate spin state. Energy decomposition analyses on the trigonal-bypiramidal (TBP) and octahedral model systems with ammonia ligands have revealed that the interaction energy of the prepared metal ion in the intermediate spin state is much smaller for the TBP structure. This sheds light on the origin of the intermediate spin state of the biomimetic TBP Fe^IV-oxo complexes.

Toward the synthesis of more reactive S = 2 non-heme oxoiron(IV) complexes

2003 marked a banner year in the bioinorganic chemistry of mononuclear non-heme iron enzymes. The first non-heme oxoiron(IV) intermediate (called J) was trapped and characterized by Bollinger and Krebs in the catalytic cycle of taurine dioxygenase (TauD), and the first crystal structure of a synthetic non-heme oxoiron(IV) complex was reported by Münck, Nam, and Que. These results stimulated inorganic chemists to synthesize related oxoiron(IV) complexes to shed light on the electronic structures and spectroscopic properties of these novel intermediates and gain mechanistic insights into their function in biology. All of the biological oxoiron(IV) intermediates discovered since 2003 have an S = 2 ground spin state, while over 90% of the 60 or so synthetic oxoiron(IV) complexes reported to date have an S = 1 ground spin state. This difference in electronic structure has fueled an interest to more accurately model these enzymatic intermediates and synthesize S = 2 oxoiron(IV) complexes.

This Account follows up on a previous Account (Acc. Chem. Res. 2007, 40, 493) that provided a perspective on the early developments in this field up to 2007 and details our group’s efforts in the development of synthetic strategies to obtain oxoiron(IV) complexes with an S = 2 ground state. Upon inspection of a qualitative d-orbital splitting diagram for a d⁴ metal–oxo center, it becomes evident that the key to achieving an S = 2 ground state is to decrease the energy gap between the d_x²–y² and d_xy orbitals. Described below are two different synthetic strategies we used to accomplish this goal.

The first strategy took advantage of the realization that the d_x²–y² and d_xy orbitals become degenerate in a C₃-symmetric ligand environment. Thus, by employing bulky tripodal ligands, trigonal-bipyramidal S = 2 oxoiron(IV) complexes were obtained. However, substrate access to the oxoiron(IV) center was hindered by the bulky ligands, and the complexes showed limited ability to cleave substrate C–H bonds. The second strategy entailed introducing weaker-field equatorial ligands in six-coordinate oxoiron(IV) complexes to decrease the d_x²–y²/d_xy energy gap to the point where the S = 2 ground state is favored. These pseudo-octahedral S = 2 oxoiron(IV) complexes exhibit high H-atom transfer reactivity relative to their S = 1 counterparts and shed light on the role that the spin state may play in these reactions. Among these complexes is a highly reactive species that to date represents the closest electronic and functional model of the enzymatic intermediate, TauD-J.

Nonheme Fe(IV) Oxo Complexes of Two New Pentadentate Ligands and Their Hydrogen-Atom and Oxygen-Atom Transfer Reactions

Two new pentadentate {N5} donor ligands based on the N4Py (N4Py = N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine) framework have been synthesized, viz. [N-(1-methyl-2-benzimidazolyl)methyl-N-(2-pyridyl)methyl-N-(bis-2-pyridyl methyl)amine] (L(1)) and [N-bis(1-methyl-2-benzimidazolyl)methyl-N-(bis-2-pyridylmethyl)amine] (L(2)), where one or two pyridyl arms of N4Py have been replaced by corresponding (N-methyl)benzimidazolyl-containing arms. The complexes [Fe(II)(CH3CN)(L)](2+) (L = L(1) (1); L(2) (2)) were synthesized, and reaction of these ferrous complexes with iodosylbenzene led to the formation of the ferryl complexes [Fe(IV)(O)(L)](2+) (L = L(1) (3); L(2) (4)), which were characterized by UV-vis spectroscopy, high resolution mass spectrometry, and Mössbauer spectroscopy. Complexes 3 and 4 are relatively stable with half-lives at room temperature of 40 h (L = L(1)) and 2.5 h (L = L(2)). The redox potentials of 1 and 2, as well as the visible spectra of 3 and 4, indicate that the ligand field weakens as ligand pyridyl substituents are progressively substituted by (N-methyl)benzimidazolyl moieties. The reactivities of 3 and 4 in hydrogen-atom transfer (HAT) and oxygen-atom transfer (OAT) reactions show that both complexes exhibit enhanced reactivities when compared to the analogous N4Py complex ([Fe(IV)(O)(N4Py)](2+)), and that the normalized HAT rates increase by approximately 1 order of magnitude for each replacement of a pyridyl moiety; i.e., [Fe(IV)(O)(L(2))](2+) exhibits the highest rates. The second-order HAT rate constants can be directly related to the substrate C-H bond dissociation energies. Computational modeling of the HAT reactions indicates that the reaction proceeds via a high spin transition state.

Mechanism of Oxidation of Ethane to Ethanol at Iron(IV)-Oxo Sites in Magnesium-Diluted Fe2(dobdc)

The catalytic properties of the metal-organic framework Fe2(dobdc), containing open Fe(II) sites, include hydroxylation of phenol by pure Fe2(dobdc) and hydroxylation of ethane by its magnesium-diluted analogue, Fe0.1Mg1.9(dobdc). In earlier work, the latter reaction was proposed to occur through a redox mechanism involving the generation of an iron(IV)-oxo species, which is an intermediate that is also observed or postulated (depending on the case) in some heme and nonheme enzymes and their model complexes. In the present work, we present a detailed mechanism by which the catalytic material, Fe0.1Mg1.9(dobdc), activates the strong C-H bonds of ethane. Kohn-Sham density functional and multireference wave function calculations have been performed to characterize the electronic structure of key species. We show that the catalytic nonheme-Fe hydroxylation of the strong C-H bond of ethane proceeds by a quintet single-state σ-attack pathway after the formation of highly reactive iron-oxo intermediate. The mechanistic pathway involves three key transition states, with the highest activation barrier for the transfer of oxygen from N2O to the Fe(II) center. The uncatalyzed reaction, where nitrous oxide directly oxidizes ethane to ethanol is found to have an activation barrier of 280 kJ/mol, in contrast to 82 kJ/mol for the slowest step in the iron(IV)-oxo catalytic mechanism. The energetics of the C-H bond activation steps of ethane and methane are also compared. Dehydrogenation and dissociation pathways that can compete with the formation of ethanol were shown to involve higher barriers than the hydroxylation pathway.

Modeling TauD-J: a high-spin nonheme oxoiron(IV) complex with high reactivity toward C-H bonds

High-spin oxoiron(IV) species are often implicated in the mechanisms of nonheme iron oxygenases, their C-H bond cleaving properties being attributed to the quintet spin state. However, the few available synthetic S = 2 Fe(IV)═O complexes supported by polydentate ligands do not cleave strong C-H bonds. Herein we report the characterization of a highly reactive S = 2 complex, [Fe(IV)(O)(TQA)(NCMe)](2+) (2) (TQA = tris(2-quinolylmethyl)amine), which oxidizes both C-H and C═C bonds at -40 °C. The oxidation of cyclohexane by 2 occurs at a rate comparable to that of the oxidation of taurine by the TauD-J enzyme intermediate after adjustment for the different temperatures of measurement. Moreover, compared with other S = 2 complexes characterized to date, the spectroscopic properties of 2 most closely resemble those of TauD-J. Together these features make 2 the best electronic and functional model for TauD-J to date.

How does tunneling contribute to counterintuitive H-abstraction reactivity of nonheme Fe(IV)O oxidants with alkanes?

This article addresses the intriguing hydrogen-abstraction (H-abstraction) and oxygen-transfer (O-transfer) reactivity of a series of nonheme [Fe(IV)(O)(TMC)(Lax)](z+) complexes, with a tetramethyl cyclam ligand and a variable axial ligand (Lax), toward three substrates: 1,4-cyclohexadiene, 9,10-dihydroanthracene, and triphenyl phosphine. Experimentally, O-transfer-reactivity follows the relative electrophilicity of the complexes, whereas the corresponding H-abstraction-reactivity generally increases as the axial ligand becomes a better electron donor, hence exhibiting an antielectrophilic trend. Our theoretical results show that the antielectrophilic trend in H-abstraction is affected by tunneling contributions. Room-temperature tunneling increases with increase of the electron donation power of the axial-ligand, and this reverses the natural electrophilic trend, as revealed through calculations without tunneling, and leads to the observed antielectrophilic trend. By contrast, O-transfer-reactivity, not being subject to tunneling, retains an electrophilic-dependent reactivity trend, as revealed experimentally and computationally. Tunneling-corrected kinetic-isotope effect (KIE) calculations matched the experimental KIE values only if all of the H-abstraction reactions proceeded on the quintet state (S = 2) surface. As such, the present results corroborate the initially predicted two-state reactivity (TSR) scenario for these reactions. The increase of tunneling with the electron-releasing power of the axial ligand, and the reversal of the "natural" reactivity pattern, support the "tunneling control" hypothesis (Schreiner et al., ref 19). Should these predictions be corroborated, the entire field of C-H bond activation in bioinorganic chemistry would lay open to reinvestigation.

Does a higher metal oxidation state necessarily imply higher reactivity toward H-atom transfer? A computational study of C-H bond oxidation by high-valent iron-oxo and -nitrido complexes

In this work, the reactions of C-H bond activation by two series of iron-oxo ( (Fe(IV)), (Fe(V)), (Fe(VI))) and -nitrido model complexes ( (Fe(IV)), (Fe(V)), (Fe(VI))) with a nearly identical coordination geometry but varying iron oxidation states ranging from iv to vi were comprehensively investigated using density functional theory. We found that in a distorted octahedral coordination environment, the iron-oxo species and their isoelectronic nitrido analogues feature totally different intrinsic reactivities toward C-H bond cleavage. In the case of the iron-oxo complexes, the reaction barrier monotonically decreases as the iron oxidation state increases, consistent with the gradually enhanced electrophilicity across the series. The iron-nitrido complex is less reactive than its isoelectronic iron-oxo species, and more interestingly, a counterintuitive reactivity pattern was observed, i.e. the activation barriers essentially remain constant independent of the iron oxidation states. The detailed analysis using the Polanyi principle demonstrates that the different reactivities between these two series originate from the distinct thermodynamic driving forces, more specifically, the bond dissociation energies (BDEE-Hs, E = O, N) of the nascent E-H bonds in the FeE-H products. Further decomposition of the BDEE-Hs into the electron and proton affinity components shed light on how the oxidation states modulate the BDEE-Hs of the two series.

Alternative Pathway for the Reaction Catalyzed by DNA Dealkylase AlkB from Ab Initio QM/MM Calculations

AlkB is the title enzyme of a family of DNA dealkylases that catalyze the direct oxidative dealkylation of nucleobases. The conventional mechanism for the dealkylation of N1-methyl adenine (1-meA) catalyzed by AlkB after the formation of FeIV-oxo is comprised by a reorientation of the oxo moiety, hydrogen abstraction, OH rebound from the Fe atom to the methyl adduct, and the dissociation of the resulting methoxide to obtain the repaired adenine base and formaldehyde. An alternative pathway with hydroxide as a ligand bound to the iron atom is proposed and investigated by QM/MM simulations. The results show OH- has a small impact on the barriers for the hydrogen abstraction and OH rebound steps. The effects of the enzyme and the OH- ligand on the hydrogen abstraction by the FeIV-oxo moiety are discussed in detail. The new OH rebound step is coupled with a proton transfer to the OH- ligand and results in a novel zwitterion intermediate. This zwitterion structure can also be characterized as Fe-O-C complex and facilitates the formation of formaldehyde. In contrast, for the pathway with H2O bound to iron, the hydroxyl product of the OH rebound step first needs to unbind from the metal center before transferring a proton to Glu136 or other residue/substrate. The consistency between our theoretical results and experimental findings is discussed. This study provides new insights into the oxidative repair mechanism of DNA repair by nonheme FeII and α-ketoglutarate (α-KG) dependent dioxygenases and a possible explanation for the substrate preference of AlkB.

Large equatorial ligand effects on C-H bond activation by nonheme iron(IV)-oxo complexes

In this article, we present density functional theory (DFT) calculations on the iron(IV)-oxo catalyzed methane C-H activation reactions for complexes in which the Fe(IV)═O core is surrounded by five negatively charged ligands. We found that it follows a hybrid pathway that mixes features of the classical σ- and π-pathways in quintet surfaces. These calculations show that the Fe-O-H arrangement in this hybrid pathway is bent in sharp contrast to the collinear character as observed for the classical quintet σ-pathways before. The calculations have also shown that it is the equatorial ligands that play key roles in tuning the reactivity of Fe(IV)═O complexes. The strong π-donating equatorial ligands employed in the current study cause a weak π(FeO) bond and thereby shift the electronic accepting orbitals (EAO) from the vertically orientated O pz orbital to the horizontally orientated O px. In addition, all the equatorial ligands are small in size and would therefore be expected have small steric effects upon substrate horizontal approaching. Therefore, for the small and strong π-donating equatorial ligands, the collinear Fe-O-H arrangement is not the best choice for the quintet reactivity. This study adds new element to iron(IV)-oxo catalyzed C-H bond activation reactions.

Dichotomous hydrogen atom transfer vs proton-coupled electron transfer during activation of X-H bonds (X = C, N, O) by nonheme iron-oxo complexes of variable basicity

We describe herein the hydrogen-atom transfer (HAT)/proton-coupled electron-transfer (PCET) reactivity for Fe(IV)-oxo and Fe(III)-oxo complexes (1-4) that activate C-H, N-H, and O-H bonds in 9,10-dihydroanthracene (S1), dimethylformamide (S2), 1,2-diphenylhydrazine (S3), p-methoxyphenol (S4), and 1,4-cyclohexadiene (S5). In 1-3, the iron is pentacoordinated by tris[N'-tert-butylureaylato)-N-ethylene]aminato ([H3buea](3-)) or its derivatives. These complexes are basic, in the order 3 ≫ 1 > 2. Oxidant 4, [Fe(IV)N4Py(O)](2+) (N4Py: N,N-bis(2-pyridylmethyl)bis(2-pyridyl)methylamine), is the least basic oxidant. The DFT results match experimental trends and exhibit a mechanistic spectrum ranging from concerted HAT and PCET reactions to concerted-asynchronous proton transfer (PT)/electron transfer (ET) mechanisms, all the way to PT. The singly occupied orbital along the O···H···X (X = C, N, O) moiety in the TS shows clearly that in the PCET cases, the electron is transferred separately from the proton. The Bell-Evans-Polanyi principle does not account for the observed reactivity pattern, as evidenced by the scatter in the plot of calculated barrier vs reactions driving forces. However, a plot of the deformation energy in the TS vs the respective barrier provides a clear signature of the HAT/PCET dichotomy. Thus, in all C-H bond activations, the barrier derives from the deformation energy required to create the TS, whereas in N-H/O-H bond activations, the deformation energy is much larger than the corresponding barrier, indicating the presence of a stabilizing interaction between the TS fragments. A valence bond model is used to link the observed results with the basicity/acidity of the reactants.

Hydrogen-bonding effects on the reactivity of [X-Fe(III)-O-Fe(IV)═O] (X = OH, F) complexes toward C-H bond cleavage

Complexes 1-OH and 1-F are related complexes that share similar [X-Fe(III)-O-Fe(IV)═O](3+) core structures with a total spin S of ½, which arises from antiferromagnetic coupling of an S = 5/2 Fe(III)-X site and an S = 2 Fe(IV)═O site. EXAFS analysis shows that 1-F has a nearly linear Fe(III)-O-Fe(IV) core compared to that of 1-OH, which has an Fe-O-Fe angle of ~130° due to the presence of a hydrogen bond between the hydroxo and oxo groups. Both complexes are at least 1000-fold more reactive at C-H bond cleavage than 2, a related complex with a [OH-Fe(IV)-O-Fe(IV)═O](4+) core having individual S = 1 Fe(IV) units. Interestingly, 1-F is 10-fold more reactive than 1-OH. This raises an interesting question about what gives rise to the reactivity difference. DFT calculations comparing 1-OH and 1-F strongly suggest that the H-bond in 1-OH does not significantly change the electrophilicity of the reactive Fe(IV)═O unit and that the lower reactivity of 1-OH arises from the additional activation barrier required to break its H-bond in the course of H-atom transfer by the oxoiron(IV) moiety.

A theory for bioinorganic chemical reactivity of oxometal complexes and analogous oxidants: the exchange and orbital-selection rules

Over the past decades metalloenzymes and their synthetic models have emerged as an area of increasing research interest. The metalloenzymes and their synthetic models oxidize organic molecules using oxometal complexes (OMCs), especially oxoiron(IV)-based ones. Theoretical studies have helped researchers to characterize the active species and to resolve mechanistic issues. This activity has generated massive amounts of data on the relationship between the reactivity of OMCs and the transition metal's identity, oxidation state, ligand sphere, and spin state. Theoretical studies have also produced information on transition state (TS) structures, reaction intermediates, barriers, and rate-equilibrium relationships. For example, the experimental-theoretical interplay has revealed that nonheme enzymes carry out H-abstraction from strong C-H bonds using high-spin (S = 2) oxoiron(IV) species with four unpaired electrons on the iron center. However, other reagents with higher spin states and more unpaired electrons on the metal are not as reactive. Still other reagents carry out these transformations using lower spin states with fewer unpaired electrons on the metal. The TS structures for these reactions exhibit structural selectivity depending on the reactive spin states. The barriers and thermodynamic driving forces of the reactions also depend on the spin state. H-Abstraction is preferred over the thermodynamically more favorable concerted insertion into C-H bonds. Currently, there is no unified theoretical framework that explains the totality of these fascinating trends. This Account aims to unify this rich chemistry and understand the role of unpaired electrons on chemical reactivity. We show that during an oxidative step the d-orbital block of the transition metal is enriched by one electron through proton-coupled electron transfer (PCET). That single electron elicits variable exchange interactions on the metal, which in turn depend critically on the number of unpaired electrons on the metal center. Thus, we introduce the exchange-enhanced reactivity (EER) principle, which predicts the preferred spin state during oxidation reactions, the dependence of the barrier on the number of unpaired electrons in the TS, and the dependence of the deformation energy of the reactants on the spin state. We complement EER with orbital-selection rules, which predict the structure of the preferred TS and provide a handy theory of bioinorganic oxidative reactions. These rules show how EER provides a Hund's Rule for chemical reactivity: EER controls the reactivity landscape for a great variety of transition-metal complexes and substrates. Among many reactivity patterns explained, EER rationalizes the abundance of high-spin oxoiron(IV) complexes in enzymes that carry out bond activation of the strongest bonds. The concepts used in this Account might also be applicable in other areas such as in f-block chemistry and excited-state reactivity of 4d and 5d OMCs.