Theoretical investigation of C--H hydroxylation by (N4Py)Fe(IV)=O(2+): an oxidant more powerful than P450?

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.

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?

N-Heterocyclic carbenes (NHC) are commonly featured as ligands in transition metal catalysis. Recently, a cyclic system containing four NHC groups with a central iron atom was synthesized and its iron(IV)-oxo species, [Fe^IV (O)(cNHC₄ )]²⁺ , was characterized. This tetracyclic NHC ligand system may give the iron(IV)-oxo species unique catalytic properties as compared to traditional non-heme and heme iron ligand systems. Therefore, we performed a computational study on the structure and reactivity of the [Fe^IV (O)(cNHC₄ )]²⁺ complex in substrate hydroxylation and epoxidation reactions. The reactivity patterns are compared with cytochrome P450 Compound I and non-heme iron(IV)-oxo models and it is shown that the [Fe^IV (O)(cNHC₄ )]²⁺ system is an effective oxidant with oxidative power analogous to P450 Compound I. Unfortunately, in polar solvents, a solvent molecule will bind to the sixth ligand position and decrease the catalytic activity of the oxidant. A molecular orbital and valence bond analysis provides insight into the origin of the reactivity differences and makes predictions of how to further exploit these systems in chemical catalysis.

Computational analysis of site differences in selective aliphatic C–H hydroxylation by nonheme iron–oxo complexes

The CH₃CN-solvent influences selective C–H hydroxylation by nonheme iron complexes due to interactions in the H-bonding formation in H-abstraction.

Proton-Coupled Electron Transfer in Organic Synthesis: Fundamentals, Applications, and Opportunities

Proton-coupled electron transfers (PCETs) are unconventional redox processes in which both protons and electrons are exchanged, often in a concerted elementary step. While PCET is now recognized to play a central a role in biological redox catalysis and inorganic energy conversion technologies, its applications in organic synthesis are only beginning to be explored. In this chapter, we aim to highlight the origins, development, and evolution of the PCET processes most relevant to applications in organic synthesis. Particular emphasis is given to the ability of PCET to serve as a non-classical mechanism for homolytic bond activation that is complimentary to more traditional hydrogen atom transfer processes, enabling the direct generation of valuable organic radical intermediates directly from their native functional group precursors under comparatively mild catalytic conditions. The synthetically advantageous features of PCET reactivity are described in detail, along with examples from the literature describing the PCET activation of common organic functional groups.

Vanadium bromoperoxidase (VBrPO) mimics: synthesis, structure and a comparative account of the catalytic activity of newly synthesized oxidovanadium and oxido-peroxidovanadium complexes

The catalytic peroxidative bromination (VBrPO mimic) of organic substances by oxido and oxido-peroxidovanadium complexes suggests the superiority of the oxido-peroxido complex over the oxido one.

What factors influence the reactivity of C-H hydroxylation and C=C epoxidation by [Fe(IV)(L(ax))(1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane)(O)](n+)

Density functional theory is used to investigate geometric structures and mechanisms for hydroxylation and epoxidation from propene for four non-heme iron complexes, [Fe(IV)(L(ax))(TMC)(O)](n+), which are the inverted isomers of [Fe(IV)(O)(TMC)(Lax)](n+) (Lax = acetonitrile (AN), monoanionic trifluoroacetate (TF), azide (N3), thiolate (SR)). The Fe(IV)O unit is found to be sterically less hindered in [Fe(IV)(L(ax))(TMC)(O)](n+) than that in [Fe(IV)(O)(TMC)(L(ax))](n+). Becke, three-parameter, Lee-Yang-Parr (B3LYP) calculations show that hydroxylation and epoxidation proceed via a two-state-reactivity on competing triplet and quintet spin surfaces; and the reactions have been invariably mediated by the S = 2 state. The reaction pathways computed reveal that 2-AN is the most reactive in the four [Fe(IV)(L(ax))(TMC)(O)](n+) complexes; along the reaction pathway, the axial ligand moves away from the iron center, and thus, the energy of the LUMO decreases. The anionic axial ligand, which is more electron releasing than neutral AN, shows a strong overlap of iron orbitals. Thus, the anionic ligand does not move away from the iron center. The H-abstraction is affected by the tunneling contribution, the more electron donation power of the axial ligand, the more effect of the tunneling contribution. Adding the tunneling correction, the relative reactivity of the hydroxylation follows the trend: 2-AN > 2-SR ≈ 2-N3 > 2-TF. However, for the epoxidation, the reactivity is in the following order of 2-AN > 2-TF > 2-N3 > 2-SR. Except for 2-AN, 2-X (L(ax) = TF, N3, SR) complexes chemoselectively hydroxylate even in the presence of a C=C double bond.

Mechanistic elucidation of C-H oxidation by electron rich non-heme iron(IV)-oxo at room temperature

Non-heme iron(IV)-oxo species form iron(III) intermediates during hydrogen atom abstraction (HAA) from the C-H bond. While synthesizing a room temperature stable, electron rich, non-heme iron(IV)-oxo compound, we obtained iron(III)-hydroxide, iron(III)-alkoxide and hydroxylated-substrate-bound iron(II) as the detectable intermediates. The present study revealed that a radical rebound pathway was operative for benzylic C-H oxidation of ethylbenzene and cumene. A dissociative pathway for cyclohexane oxidation was established based on UV-vis and radical trap experiments. Interestingly, experimental evidence including O-18 labeling and mechanistic study suggested an electron transfer mechanism to be operative during C-H oxidation of alcohols (e.g. benzyl alcohol and cyclobutanol). The present report, therefore, unveils non-heme iron(IV)-oxo promoted substrate-dependent C-H oxidation pathways which are of synthetic as well as biological significance.

Oxoiron(IV) Complex of the Ethylene-Bridged Dialkylcyclam Ligand Me2EBC

We report herein the first example of an oxoiron(IV) complex of an ethylene-bridged dialkylcyclam ligand, [Fe(IV)(O)(Me2EBC)(NCMe)](2+) (2; Me2EBC = 4,11-dimethyl-1,4,8,11-tetraazabicyclo[6.6.2]hexadecane). Complex 2 has been characterized by UV-vis, (1)H NMR, resonance Raman, Mössbauer, and X-ray absorption spectroscopy as well as electrospray ionization mass spectrometry, and its properties have been compared with those of the closely related [Fe(IV)(O)(TMC)(NCMe)](2+) (3; TMC = 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane), the intensively studied prototypical oxoiron(IV) complex of the macrocyclic tetramethylcyclam ligand. Me2EBC has an N4 donor set nearly identical with that of TMC but possesses an ethylene bridge in place of the 1- and 8-methyl groups of TMC. As a consequence, Me2EBC is forced to deviate from the trans-I configuration typically found for Fe(IV)(O)(TMC) complexes and instead adopts a folded cis-V stereochemistry that requires the MeCN ligand to coordinate cis to the Fe(IV)═O unit in 2 rather than in the trans arrangement found in 3. However, switching from the trans geometry of 3 to the cis geometry of 2 did not significantly affect their ground-state electronic structures, although a decrease in ν(Fe═O) was observed for 2. Remarkably, despite having comparable Fe(IV/III) reduction potentials, 2 was found to be significantly more reactive than 3 in both oxygen-atom-transfer (OAT) and hydrogen-atom-transfer (HAT) reactions. A careful analysis of density functional theory calculations on the HAT reactivity of 2 and 3 revealed the root cause to be the higher oxyl character of 2, leading to a stronger O---H bond specifically in the quintet transition state.

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.

Properties and reactivities of nonheme iron(IV)-oxo versus iron(V)-oxo: long-range electron transfer versus hydrogen atom abstraction

Recent work of Nam and co-workers [J. Yoon, S. A. Wilson, Y. K. Jang, M. S. Seo, K. Nehru, B. Hedman, K. O. Hodgson, E. Bill, E. I. Solomon and W. Nam, Angew. Chem., Int. Ed., 2009, 48, 1257] on a biomimetic iron complex implicated a mixture of iron(IV)-oxo and iron(V)-oxo intermediates but the latter could not be spectroscopically characterized, hence its involvement was postulated. To gain insight into the relative activity of these iron(IV)-oxo versus iron(V)-oxo intermediates, we have performed an extensive density functional theory (DFT) study on the chemical properties of the chemical system of Nam et al., namely [Fe(O)(BQEN)(NCCH3)](2+/3+) with BQEN = N,N'-dimethyl-N,N'-bis(8-quinolyl)ethane-1,2-diamine and their reactivity in hydrogen atom abstraction from ethylbenzene. We show that the perceived iron(V)-oxo species actually is an iron(IV)-oxo ligand cation radical, similar to cytochrome P450 compound I. Moreover, this intermediate has an extremely large electron affinity and therefore can abstract electrons from substrates readily. In our particular system, this means that prior to the hydrogen atom abstraction, an electron is abstracted to form an iron(IV)-oxo species, which subsequently abstracts a hydrogen atom from the substrate. Thus, our calculations show for the first time how some nonheme iron complexes react by long-range electron transfer and others directly via hydrogen atom abstraction. We have rationalized our results with detailed thermochemical cycles that explain the observed reactivity patterns.

ortho-Hydroxylation of aromatic acids by a non-heme Fe(V)=O species: how important is the ligand design?

There is a growing interest in probing the mechanism of catalytic transformations effected by non-heme iron-oxo complexes as these reactions set a platform for understanding the relevant enzymatic reactions. The ortho-hydroxylation of aromatic compounds is one such reaction catalysed by iron-oxo complexes. Experimentally [Fe(II)(BPMEN)(CH3CN)2](2+) (1) and [Fe(II)(TPA)(CH3CN)2](2+) (2) (where TPA = tris(2-pyridylmethyl)amine and BPMEN = N,N′-dimethyl-N,N′-bis(2-pyridylmethyl)ethane-1,2-diamine) complexes containing amino pyridine ligands along with H2O2 are employed to carry out these transformations where complex 1 is found to be more reactive than complex 2. Herein, using density functional methods employing B3LYP and dispersion corrected B3LYP (B3LYP-D) functionals, we have explored the mechanism of this reaction to reason out the importance of ligand design in fine-tuning the reactivity of such catalytic transformations. Dispersion corrected B3LYP is found to be superior to B3LYP in predicting the correct ground state of these species and also yields lower barrier heights than the B3LYP functional. Starting the reaction from the Fe(III)–OOH species, both homolytic and heterolytic cleavage of the O···O bond is explored leading to the formation of the transient Fe(IV)=O and Fe(V)=O species. For both the ligand systems, heterolytic cleavage was energetically preferable and our calculations suggest that both the reactions are catalyzed by an elusive high-valent Fe(V)=O species. The Fe(V)=O species undergoes the reaction via an electrophilic attack of the benzene ring to effect the ortho-hydroxylation reaction. The reactivity pattern observed for 1 and 2 are reflected in the computed barrier heights for the ortho-hydroxylation reaction. Electronic structure analysis reveals that the difference in reactivity between the ligand architectures described in complex 1 and 2 arise due to orientation of the pyridine ring(s) parallel or perpendicular to the Fe(V)=O bond. The parallel orientation of the pyridine ring is found to mix with the (πFe(dyz)–O(py))* orbital of the Fe-oxo bond leading to a reduction in the electrophilicity of the ferryl oxygen atom. Our calculations highlight the importance of ligand design in this chemistry and suggest that this concept can be used to (i) stabilize high-valent intermediates which can be trapped and thoroughly characterized (ii) enhance the reactivity and efficiency of the oxidants by increasing the electrophilicity of the ferryl oxygen containing FeVO species. Our computed results are in general agreement with the experimental results.

Determination of Spin Inversion Probability, H-Tunneling Correction, and Regioselectivity in the Two-State Reactivity of Nonheme Iron(IV)-Oxo Complexes

We show by experiments that nonheme Fe(IV)O species react with cyclohexene to yield selective hydrogen atom transfer (HAT) reactions with virtually no C═C epoxidation. Straightforward DFT calculations reveal, however, that C═C epoxidation on the S = 2 state possesses a low-energy barrier and should contribute substantially to the oxidation of cyclohexene by the nonheme Fe(IV)O species. By modeling the selectivity of this two-site reactivity, we show that an interplay of tunneling and spin inversion probability (SIP) reverses the apparent barriers and prefers exclusive S = 1 HAT over mixed HAT and C═C epoxidation on S = 2. The model enables us to derive a SIP value by combining experimental and theoretical results.

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.

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

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

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.

Quantum mechanics/molecular mechanics study on the oxygen binding and substrate hydroxylation step in AlkB repair enzymes

AlkB repair enzymes are important nonheme iron enzymes that catalyse the demethylation of alkylated DNA bases in humans, which is a vital reaction in the body that heals externally damaged DNA bases. Its mechanism is currently controversial and in order to resolve the catalytic mechanism of these enzymes, a quantum mechanics/molecular mechanics (QM/MM) study was performed on the demethylation of the N(1) -methyladenine fragment by AlkB repair enzymes. Firstly, the initial modelling identified the oxygen binding site of the enzyme. Secondly, the oxygen activation mechanism was investigated and a novel pathway was found, whereby the catalytically active iron(IV)-oxo intermediate in the catalytic cycle undergoes an initial isomerisation assisted by an Arg residue in the substrate binding pocket, which then brings the oxo group in close contact with the methyl group of the alkylated DNA base. This enables a subsequent rate-determining hydrogen-atom abstraction on competitive σ- and π-pathways on a quintet spin-state surface. These findings give evidence of different locations of the oxygen and substrate binding channels in the enzyme and the origin of the separation of the oxygen-bound intermediates in the catalytic cycle from substrate. Our studies are compared with small model complexes and the effect of protein and environment on the kinetics and mechanism is explained.

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.

Reactivity of compound II: electronic structure analysis of methane hydroxylation by oxoiron(IV) porphyrin complexes

The methane hydroxylation reaction by a Compound II (Cpd II) mimic PorFe(IV)=O and its hydrosulfide-ligated derivative [Por(SH)Fe(IV)=O](-) is investigated by density functional theory (DFT) calculations on the ground triplet and excited quintet spin-state surfaces. On each spin surface both the σ- and π-channels are explored. H-abstraction is invariably the rate-determining step. In the case of PorFe(IV)=O the H-abstraction reaction can proceed either through the classic π-channel or through the nonclassical σ-channel on the triplet surface, but only through the classic σ-mechanism on the quintet surface. The barrier on the quintet σ-pathway is much lower than on the triplet channels so the quintet surface cuts through the triplet surfaces and a two state reactivity (TSR) mechanism with crossover from the triplet to the quintet surface becomes a plausible scenario for C-H bond activation by PorFe(IV)=O. In the case of the hydrosulfide-ligated complex the H-abstraction follows a π-mechanism on the triplet surface: the σ* is too high in energy to make a σ-attack of the substrate favorable. The σ- and π-channels are both feasible on the quintet surface. As the quintet surface lies above the triplet surface in the entrance channel of the oxidative process and is highly destabilized on both the σ- and π-pathways, the reaction can only proceed on the triplet surface. Insights into the electron transfer process accompanying the H-abstraction reaction are achieved through a detailed electronic structure analysis of the transition state species and the reactant complexes en route to the transition state. It is found that the electron transfer from the substrate σ(CH) into the acceptor orbital of the catalyst, the Fe-O σ* or π*, occurs through a rather complex mechanism that is initiated by a two-orbital four-electron interaction between the σ(CH) and the low-lying, oxygen-rich Fe-O σ-bonding and/or Fe-O π-bonding orbitals of the catalyst.

Intramolecular gas-phase reactions of synthetic nonheme oxoiron(IV) ions: proximity and spin-state reactivity rules

The intramolecular gas-phase reactivity of four oxoiron(IV) complexes supported by tetradentate N(4) ligands (L) has been studied by means of tandem mass spectrometry measurements in which the gas-phase ions [Fe(IV)(O)(L)(OTf)](+) (OTf = trifluoromethanesulfonate) and [Fe(IV) (O)(L)](2+) were isolated and then allowed to fragment by collision-induced decay (CID). CID fragmentation of cations derived from oxoiron(IV) complexes of 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane (tmc) and N,N'-bis(2-pyridylmethyl)-1,5-diazacyclooctane (L(8)Py(2)) afforded the same predominant products irrespective of whether they were hexacoordinate or pentacoordinate. These products resulted from the loss of water by dehydrogenation of ethylene or propylene linkers on the tetradentate ligand. In contrast, CID fragmentation of ions derived from oxoiron(IV) complexes of linear tetradentate ligands N,N'-bis(2-pyridylmethyl)-1,2-diaminoethane (bpmen) and N,N'-bis(2-pyridylmethyl)-1,3-diaminopropane (bpmpn) showed predominant oxidative N-dealkylation for the hexacoordinate [Fe(IV)(O)(L)(OTf)](+) cations and predominant dehydrogenation of the diaminoethane/propane backbone for the pentacoordinate [Fe(IV)(O)(L)](2+) cations. DFT calculations on [Fe(IV)(O)(bpmen)] ions showed that the experimentally observed preference for oxidative N-dealkylation versus dehydrogenation of the diaminoethane linker for the hexa- and pentacoordinate ions, respectively, is dictated by the proximity of the target C-H bond to the oxoiron(IV) moiety and the reactive spin state. Therefore, there must be a difference in ligand topology between the two ions. More importantly, despite the constraints on the geometries of the TS that prohibit the usual upright σ trajectory and prevent optimal σ(CH)-σ*(z2) overlap, all the reactions still proceed preferentially on the quintet (S = 2) state surface, which increases the number of exchange interactions in the d block of iron and leads thereby to exchange enhanced reactivity (EER). As such, EER is responsible for the dominance of the S = 2 reactions for both hexa- and pentacoordinate complexes.