Origin of the Enhanced Reactivity of μ-Nitrido-Bridged Diiron(IV)-Oxo Porphyrinoid Complexes over Cytochrome P450 Compound I

Biomimetic [MFe3S4]3+ Cubanes (M = V/Mo) as Catalysts for a Fischer-Tropsch-like Hydrocarbon Synthesis─A Computational Study

Nitrogenase is the enzyme primarily responsible for reducing atmospheric nitrogen to ammonia. There are three general forms of nitrogenase based on the metal ion present in the cofactor binding site, namely, molybdenum-dependent nitrogenases with the iron-molybdenum cofactor (FeMoco), the vanadium-dependent nitrogenases with FeVco, and the iron-only nitrogenases. It has been shown that the vanadium-dependent nitrogenases tend to have a lesser efficacy in reducing dinitrogen but a higher efficacy in binding and reducing carbon monoxide. In biomimetic chemistry, [MFe3S4] (M = Mo/V) cubanes have been synthesized, studied, and shown to be promising mimics of some of the geometric and electronic properties of the nitrogenase cofactors. In this work, a density functional theory (DFT) study is presented on Fischer-Tropsch catalysis by these cubane complexes by studying CO binding and reduction to hydrocarbons. Our work implies that molybdenum has stronger binding interactions with the iron-sulfur framework of the cubane, which results in easier reduction of substrates like N2H4. However, this inhibits the binding and activation of CO, and hence, the molybdenum-containing complexes are less suitable for Fischer-Tropsch catalysis than vanadium-containing complexes.

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

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

Origin of Unprecedented Formation and Reactivity of FeIV═O Species via Oxygen Activation: Role of Noncovalent Interactions and Magnetic Coupling

Emulating the capabilities of the soluble methane monooxygenase (sMMO) enzymes, which effortlessly activate oxygen at diiron(II) centers to form a reactive diiron(IV) intermediate Q, which then performs the challenging oxidation of methane to methanol, poses a significant challenge. Very recently, one of us reported the mononuclear complex [(cyclam)FeII(CH3CN)2]2+ (1), which performed a rare bimolecular activation of the molecule of O2 to generate two molecules of FeIV═O without the requirement of external proton or electron sources, similar to sMMO. In the present study, we employed the density functional theory (DFT) calculations to investigate this unique mechanism of O2 activation. We show that secondary hydrogen-bonding interactions between ligand N-H groups and O2 play a vital role in reducing the energy barrier associated with the initial O2 binding at 1 and O-O bond cleavage to form the FeIV═O complex. Further, the unique reactivity of FeIV═O species toward simultaneous C-H and O-H bond activation process has been demonstrated. Our study unveils that the nature of the magnetic coupling between the diiron centers is also crucial. Given that the influence of magnetic coupling and noncovalent interactions in catalysis remains largely unexplored, this unexplored realm presents numerous avenues for experimental chemists to develop novel structural and functional analogues of sMMO.

Effect of porphyrin ligands on the catalytic CH4 oxidation activity of monocationic μ-nitrido-bridged iron porphyrinoid dimers by using H2O2 as an oxidant

The μ-nitrido-bridged iron phthalocyanine homodimer is a potent molecule-based CH4 oxidation catalyst that can effectively oxidize chemically stable CH4 under mild reaction conditions in an acidic aqueous solution including an oxidant such as H2O2. The reactive intermediate is a high-valent iron-oxo species generated upon reaction with H2O2. However, a detailed comparison of the CH4 oxidation activity of the μ-nitrido-bridged iron phthalocyanine dimer with those of μ-nitrido-bridged iron porphyrinoid dimers containing one or two porphyrin ring(s) has not been yet reported, although porphyrins are the most important class of porphyrinoids. Herein, we compare the catalytic CH4 and CH3CH3 oxidation activities of a monocationic μ-nitrido-bridged iron porphyrin homodimer and a monocationic μ-nitrido-bridged heterodimer of an iron porphyrin and an iron phthalocyanine with those of a monocationic μ-nitrido-bridged iron phthalocyanine homodimer in an acidic aqueous solution containing H2O2 as an oxidant. It was demonstrated that the CH4 oxidation activities of monocationic μ-nitrido-bridged iron porphyrinoid dimers containing porphyrin ring(s) were much lower than that of a monocationic μ-nitrido-bridged iron phthalocyanine homodimer. These findings suggested that the difference in the electronic structure of the porphyrinoid rings of monocationic μ-nitrido-bridged iron porphyrinoid dimers strongly affected their catalytic light alkane oxidation activities.

Unravelling the Effect of Acid-Driven Electron Transfer in High-Valent FeIV =O/MnIV =O Species and Its Implications for Reactivity

The electron transfer (ET) step is one of the crucial processes in biochemical redox reactions that occur in nature and has been established as a key step in dictating the reactivity of high-valent metal-oxo species. Although metalloenzymes possessing metal-oxo units at their active site are typically associated with outer-sphere electron transfer (OSET) processes, biomimetic models, in contrast, have been found to manifest either an inner-sphere electron transfer (ISET) or OSET mechanism. This distinction is clearly illustrated through the behaviour of [(N4Py)MnIV (O)]2+ (1) and [(N4Py)FeIV (O)]2+ (2) complexes, where complex 1 showcases an OSET mechanism, while complex 2 exhibits an ISET mechanism, especially evident in their reactions involving C-H bond activation and oxygen atom transfer reactions in the presence of a Lewis/Bronsted acid. However, the precise reason for this puzzling difference remains elusive. This work unveils the origin of the perplexing inner-sphere vs outer-sphere electron transfer process (ISET vs OSET) in [(N4Py)MnIV (O)]2+ (1) and [(N4Py)FeIV (O)]2+ (2) species in the presence of Bronsted acid. The calculations indicate that when the substrate (toluene) approaches both 1 and 2 that is hydrogen bonded with two HOTf molecules (denoted as 1-HOTf and 2-HOTf, respectively), proton transfer from one of the HOTf molecules to the metal-oxo unit is triggered and a simultaneous electron transfer occurs from toluene to the metal centre. Interestingly, the preference for OSET by 1-HOTf is found to originate from the choice of MnIV =O centre to abstract spin-down (β) electron from toluene to its δ(dxy ) orbital. On the other hand, in 2-HOTf, a spin state inversion from triplet to quintet state takes place during the proton (from HOTf) coupled electron transfer (from toluene) preferring a spin-up (α) electron abstraction to its σ* (dz 2 ) orbital mediated by HOTf giving rise to ISET. In addition, 2-HOTf was calculated to possess a larger reorganisation energy, which facilitates the ISET process via the acid. The absence of spin-inversion and smaller reorganisation energy switch the mechanism to OSET for 1-HOTf. Therefore, for the first time, the significance of spin-state and spin-inversion in the electron transfer process has been identified and demonstrated within the realm of high-valent metal-oxo chemistry. This discovery holds implications for the potential involvement of high-valent Mn-oxo species in performing similar transformative processes within Photosystem II.

Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides

Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations.

Bioinspired microenvironment modulation of metal-organic framework-based catalysts for selective methane oxidation

Inspiration from natural enzymes enabling creationary catalyst design is appealing yet remains extremely challenging for selective methane (CH4) oxidation. This study presents the construction of a biomimetic catalyst platform for CH4 oxidation, which is constructed by incorporating Fe-porphyrin into a robust metal-organic framework, UiO-66, furnished with saturated monocarboxylic fatty acid bearing different long alkyl chains. The catalysts demonstrate the high efficiency in the CH4 to methanol (CH3OH) conversion at 50 °C. Moreover, the selectivity to CH3OH can be effectively regulated and promoted through a fine-tuned microenvironment by hydrophobic modification around the Fe-porphyrin. The long-chain fatty acids anchored on the Zr-oxo cluster of UiO-66 can not only tune the electronic state of the Fe sites to improve CH4 adsorption, but also restrict the amount of H2O2 around the Fe sites to reduce the overoxidation. This behavior resembles the microenvironment regulation in methane monooxygenase, resulting in high CH3OH selectivity.

Melatonin Activation by Cytochrome P450 Isozymes: How Does CYP1A2 Compare to CYP1A1?

Cytochrome P450 enzymes are versatile enzymes found in most biosystems that catalyze mono-oxygenation reactions as a means of biosynthesis and biodegradation steps. In the liver, they metabolize xenobiotics, but there are a range of isozymes with differences in three-dimensional structure and protein chain. Consequently, the various P450 isozymes react with substrates differently and give varying product distributions. To understand how melatonin is activated by the P450s in the liver, we did a thorough molecular dynamics and quantum mechanics study on cytochrome P450 1A2 activation of melatonin forming 6-hydroxymelatonin and N-acetylserotonin products through aromatic hydroxylation and O-demethylation pathways, respectively. We started from crystal structure coordinates and docked substrate into the model, and obtained ten strong binding conformations with the substrate in the active site. Subsequently, for each of the ten substrate orientations, long (up to 1 μs) molecular dynamics simulations were run. We then analyzed the orientations of the substrate with respect to the heme for all snapshots. Interestingly, the shortest distance does not correspond to the group that is expected to be activated. However, the substrate positioning gives insight into the protein residues it interacts with. Thereafter, quantum chemical cluster models were created and the substrate hydroxylation pathways calculated with density functional theory. These relative barrier heights confirm the experimental product distributions and highlight why certain products are obtained. We make a detailed comparison with previous results on CYP1A1 and identify their reactivity differences with melatonin.

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

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

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

DFT Quantum-chemical prediction of molecular structure of iron(VI) macrocyclic complex with phthalocyanine and two oxo ligands

The possibility of the existence of the unknown iron complex containing phthalocyanine and two oxygen atoms in the inner coordination sphere was predicted using DFT quantum chemical calculation (OPBE/TZVP and B3PW91/TZVP levels). The structural parameters, the ground state multiplicity, the NBO analysis results, and the standard thermodynamic parameters for complex formation (standard enthalpy [Formula: see text], entropy [Formula: see text] and Gibbs’s energy [Formula: see text] are represented.

Organometallic μ-Nitridodiiron Complexes in Oxidation States Ranging from (III/III) to (IV/IV)

Iron complexes with nitrido ligands are of interest as molecular analogues of key intermediates during N₂-to-NH₃ conversion in industrial or enzymatic processes. Dinuclear iron complexes with a bridging nitrido unit are mostly known in relatively high oxidation states (III/IV or IV/IV), originating from the decomposition of azidoiron precursors via high-valent Fe≡N intermediates. The use of a tetra-NHC macrocyclic scaffold ligand (NHC = N-heterocyclic carbene) has now allowed for the isolation of a series of organometallic μ-nitridodiiron complexes ranging from the mid-valent Fe^III-N-Fe^III (1) via mixed-valent Fe^III-N-Fe^IV (type 4) to the high-valent Fe^IV-N-Fe^IV (type 5) species that are interconverted at moderate potentials, accompanied by axial ligand binding at the Fe^IV sites. Magnetic measurements and electron paramagnetic resonance spectroscopy showed the homovalent complexes to be diamagnetic and the mixed-valent system to feature an S = 1/2 ground state due to very strong antiferromagnetic coupling. The bonding in the Fe-N-Fe moiety has been further probed by crystallographic structure determination, ⁵⁷Fe Mössbauer and UV-vis spectroscopies, as well as density functional theory computations, which revealed high covalency and nearly identical Fe-N distances across this redox series. The latter has been rationalized in terms of the nonbonding nature of the combination of Fe d_z² atomic orbitals from which electrons are successively removed upon oxidation, and these redox processes are best described as being metal-centered. The tetra-NHC-ligated μ-nitridodiiron series complements a set of related complexes with single-atom μ-oxido and μ-phosphido bridges, but the Fe-N-Fe core exhibits a comparatively high stability over several oxidation states. This promises interesting applications in view of the manifold catalytic uses of μ-nitridodiiron complexes based on macrocyclic N-donor porphinato(2-) or phthalocyaninato(2-) ligands.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65, which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

Activation and catalytic transformation of methane under mild conditions

In the last few decades, worldwide scientists have been motivated by the promising production of chemicals from the widely existing methane (CH₄) under mild conditions for both chemical synthesis with low energy consumption and climate remediation. To achieve this goal, a whole library of catalytic chemistries of transforming CH₄ to various products under mild conditions is required to be developed. Worldwide scientists have made significant efforts to reach this goal. These significant efforts have demonstrated the feasibility of oxidation of CH₄ to value-added intermediate compounds including but not limited to CH₃OH, HCHO, HCOOH, and CH₃COOH under mild conditions. The fundamental understanding of these chemical and catalytic transformations of CH₄ under mild conditions have been achieved to some extent, although currently neither a catalyst nor a catalytic process can be used for chemical production under mild conditions at a large scale. In the academic community, over ten different reactions have been developed for converting CH₄ to different types of oxygenates under mild conditions in terms of a relatively low activation or catalysis temperature. However, there is still a lack of a molecular-level understanding of the activation and catalysis processes performed in extremely complex reaction environments under mild conditions. This article reviewed the fundamental understanding of these activation and catalysis achieved so far. Different oxidative activations of CH₄ or catalytic transformations toward chemical production under mild conditions were reviewed in parallel, by which the trend of developing catalysts for a specific reaction was identified and insights into the design of these catalysts were gained. As a whole, this review focused on discussing profound insights gained through endeavors of scientists in this field. It aimed to present a relatively complete picture for the activation and catalytic transformations of CH₄ to chemicals under mild conditions. Finally, suggestions of potential explorations for the production of chemicals from CH₄ under mild conditions were made. The facing challenges to achieve high yield of ideal products were highlighted and possible solutions to tackle them were briefly proposed.

Synthesis of a monocationic μ-nitrido-bridged iron porphycene dimer and its methane oxidation activity

Herein, we report the synthesis of a monocationic μ-nitrido-bridged iron porphycene dimer, a structural analogue of a monocationic μ-nitrido-bridged iron phthalocyanine dimer, which is known to be one of the most potent molecule-based catalysts for methane oxidation. ¹H-NMR and single-crystal X-ray structural analyses showed that the porphycene complex includes two Fe(IV) ions, and the structure around the Fe-NFe core is quite similar to that of the monocationic μ-nitrido-bridged iron phthalocyanine dimer. Although methane was oxidized into MeOH, HCHO, and HCOOH in the presence of a silica-supported catalyst of this monocationic μ-nitrido-bridged iron porphycene dimer in an acidic aqueous solution containing excess H₂O₂, its reactive intermediate was not a high-valence iron-oxo species, as in the case of a monocationic μ-nitrido-bridged iron phthalocyanine dimer, but ˙OH. It is suggested that the high-valent iron-oxo species of the μ-nitrido-bridged iron porphycene dimer was gradually decomposed under these reaction conditions, and the decomposed compound catalyzed a Fenton-type reaction. This result indicates that the stability of the oxo-species is indispensable for achieving high catalytic methane oxidation activity using a μ-nitrido-bridged iron porphyrinoid dimer with an Fe-NFe core as a catalyst.

A Noncanonical Tryptophan Analogue Reveals an Active Site Hydrogen Bond Controlling Ferryl Reactivity in a Heme Peroxidase

Nature employs high-energy metal-oxo intermediates embedded within enzyme active sites to perform challenging oxidative transformations with remarkable selectivity. Understanding how different local metal-oxo coordination environments control intermediate reactivity and catalytic function is a long-standing objective. However, conducting structure-activity relationships directly in active sites has proven challenging due to the limited range of amino acid substitutions achievable within the constraints of the genetic code. Here, we use an expanded genetic code to examine the impact of hydrogen bonding interactions on ferryl heme structure and reactivity, by replacing the N-H group of the active site Trp51 of cytochrome c peroxidase by an S atom. Removal of a single hydrogen bond stabilizes the porphyrin π-cation radical state of CcP W191F compound I. In contrast, this modification leads to more basic and reactive neutral ferryl heme states, as found in CcP W191F compound II and the wild-type ferryl heme-Trp191 radical pair of compound I. This increased reactivity manifests in a >60-fold activity increase toward phenolic substrates but remarkably has negligible effects on oxidation of the biological redox partner cytc. Our data highlight how Trp51 tunes the lifetimes of key ferryl intermediates and works in synergy with the redox active Trp191 and a well-defined substrate binding site to regulate catalytic function. More broadly, this work shows how noncanonical substitutions can advance our understanding of active site features governing metal-oxo structure and reactivity.

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 H2O2 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 OleTJE, 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 OleTJE 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 OleTJE guides a reaction to a chemoselectivity.

Electrostatic Perturbations from the Protein Affect C-H Bond Strengths of the Substrate and Enable Negative Catalysis in the TmpA Biosynthesis Enzyme

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C¹ -position of the substrate as part of the biosynthesis of glycine betaine in bacteria and carnitine in humans. To understand how the enzyme avoids breaking the weak C² -H bond in favor of C¹ -hydroxylation, we set up a cluster model of 242 atoms representing the first and second coordination sphere of the metal center and substrate binding pocket, and investigated possible reaction mechanisms of substrate activation by an iron(IV)-oxo species by density functional theory methods. In agreement with experimental product distributions, the calculations predict a favorable C¹ -hydroxylation pathway. The calculations show that the selectivity is guided through electrostatic perturbations inside the protein from charged residues, external electric fields and electric dipole moments. In particular, charged residues influence and perturb the homolytic bond strength of the C¹ -H and C² -H bonds of the substrate, and strongly strengthens the C² -H bond in the substrate-bound orientation.

High catalytic methane oxidation activity of monocationic μ-nitrido-bridged iron phthalocyanine dimer with sixteen methyl groups

Herein, we report the highly potent catalytic methane oxidation activity of a monocationic μ-nitrido-bridged iron phthalocyanine dimer with 16 peripheral methyl groups. It was confirmed that this complex oxidized methane stably into MeOH, HCHO, and HCOOH in a catalytic manner in an acidic aqueous solution containing excess H2O2 at 60 °C. The total turnover number of the reaction reached 135 after 12 h, which is almost seven times higher than that of a monocatinoic μ-nitrido-bridged iron phthalocyanine dimer with no peripheral substituents. This suggests that the increased number of peripheral electron-donating substituents could have facilitated the generation of a reactive high-valent iron-oxo species as well as hydrogen abstraction from methane by the reactive iron-oxo species.

A comprehensive insight into aldehyde deformylation: mechanistic implications from biology and chemistry

Aldehyde deformylation is an important reaction in biology, organic chemistry and inorganic chemistry and the process has been widely applied and utilized. For instance, in biology, the aldehyde deformylation reaction has wide differences in biological function, whereby cyanobacteria convert aldehydes into alkanes or alkenes, which are used as natural products for, e.g., defense mechanisms. By contrast, the cytochromes P450 catalyse the biosynthesis of hormones, such as estrogen, through an aldehyde deformylation reaction step. In organic chemistry, the aldehyde deformylation reaction is a common process for replacing functional groups on a molecule, and as such, many different synthetic methods and procedures have been reported that involve an aldehyde deformylation step. In bioinorganic chemistry, a variety of metal(iii)-peroxo complexes have been synthesized as biomimetic models and shown to react efficiently with aldehydes through deformylation reactions. This review paper provides an overview of the various aldehyde deformylation reactions in organic chemistry, biology and biomimetic model systems, and shows a broad range of different chemical reaction mechanisms for this process. Although a nucleophilic attack at the carbonyl centre is the consensus reaction mechanism, several examples of an alternative electrophilic reaction mechanism starting with hydrogen atom abstraction have been reported as well. There is still much to learn and to discover on aldehyde deformylation reactions, as deciphered in this review paper.

Selective carbene transfer to amines and olefins catalyzed by ruthenium phthalocyanine complexes with donor substituents

Electron-rich ruthenium phthalocyanine complexes were evaluated in carbene transfer reactions from ethyl diazoacetate (EDA) to aromatic and aliphatic olefins as well as to a wide range of aromatic, heterocyclic and aliphatic amines for the first time. It was revealed that the ruthenium octabutoxyphthalocyanine carbonyl complex [(BuO)8Pc]Ru(CO) is the most efficient catalyst converting electron-rich and electron-poor aromatic olefins to cyclopropane derivatives with high yields (typically 80-100%) and high TON (up to 1000) under low catalyst loading and nearly equimolar substrate/EDA ratio. This catalyst shows a rare efficiency in the carbene insertion into amine N-H bonds. Using a 0.05 mol% catalyst loading, a high amine concentration (1 M) and 1.1 eq. of EDA, a number of structurally divergent amines were selectively converted to mono-substituted glycine derivatives with up to quantitative yields and turnover numbers reaching 2000. High selectivity, large substrate scope, low catalyst loading and practical reaction conditions place [(BuO)8Pc]Ru(CO) among the most efficient catalysts for the carbene insertion into amines.

Oxidation of Organic Compounds with Peroxides Catalyzed by Polynuclear Metal Compounds

The review describes articles that provide data on the synthesis and study of the properties of catalysts for the oxidation of alkanes, olefins, and alcohols. These catalysts are polynuclear complexes of iron, copper, osmium, nickel, manganese, cobalt, vanadium. Such complexes for example are: [Fe2(HPTB)(m-OH)(NO3)2](NO3)2·CH3OH·2H2O, where HPTB-¼N,N,N0,N0-tetrakis(2-benzimidazolylmethyl)-2-hydroxo-1,3-diaminopropane; complex [(PhSiO1,5)6]2[CuO]4[NaO0.5]4[dppmO2]2, where dppm-1,1-bis(diphenylphosphino)methane; (2,3-η-1,4-diphenylbut-2-en-1,4-dione)undecacarbonyl triangulotriosmium; phenylsilsesquioxane [(PhSiO1.5)10(CoO)5(NaOH)]; bi- and tri-nuclear oxidovanadium(V) complexes [{VO(OEt)(EtOH)}2(L2)] and [{VO(OMe)(H2O)}3(L3)]·2H2O (L2 = bis(2-hydroxybenzylidene)terephthalohydrazide and L3 = tris(2-hydroxybenzylidene)benzene-1,3,5-tricarbohydrazide); [Mn2L2O3][PF6]2 (L = 1,4,7-trimethyl-1,4,7-triazacyclononane). For comparison, articles are introduced describing catalysts for the oxidation of alkanes and alcohols with peroxides, which are simple metal salts or mononuclear metal complexes. In many cases, polynuclear complexes exhibit higher activity compared to mononuclear complexes and exhibit increased regioselectivity, for example, in the oxidation of linear alkanes. The review contains a description of some of the mechanisms of catalytic reactions. Additionally presented are articles comparing the rates of oxidation of solvents and substrates under oxidizing conditions for various catalyst structures, which allows researchers to conclude about the nature of the oxidizing species. This review is focused on recent works, as well as review articles and own original studies of the authors.

Neighbouring effects on catalytic epoxidation by Fe-cyclam in M2-PDIxCy complexes

The unsymmetric PDIeCy ligand, featuring pyridinediimine and cylam sites, can be selectively metalated. Complementing the diiron compound, we have synthesized two heterobimetallic isomers, [ZnPDIFeCy(PDIeCy)(OTf)4] (3) and [FePDIZnCy(PDIeCy)(OTf)4] (4), and a dizinc complex, [Zn2(PDIeCy)(OTf)4] (5). Olefin epoxidation by the series of complexes was investigated. The M-PDI site influences the reactivity of the M-cyclam, resulting in increased activity toward enones.

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.

Hydroxyl Transfer to Carbon Radicals by Mn(OH) vs Fe(OH) Corrole Complexes

The transfer of •OH from metal-hydroxo species to carbon radicals (R•) to give hydroxylated products (ROH) is a fundamental process in metal-mediated heme and nonheme C-H bond oxidations. This step, often referred to as the hydroxyl "rebound" step, is typically very fast, making direct study of this process challenging if not impossible. In this report, we describe the reactions of the synthetic models M(OH)(ttppc) (M = Fe (1), Mn (3); ttppc = 5,10,15-tris(2,4,6-triphenyl)phenyl corrolato^3-) with a series of triphenylmethyl carbon radical (R•) derivatives ((4-X-C₆H₄)₃C•; X = OMe, tBu, Ph, Cl, CN) to give the one-electron reduced M^III(ttppc) complexes and ROH products. Rate constants for 3 for the different radicals ranged from 11.4(1) to 58.4(2) M^-1 s^-1, as compared to those for 1, which fall between 0.74(2) and 357(4) M^-1 s^-1. Linear correlations for Hammett and Marcus plots for both Mn and Fe were observed, and the small magnitudes of the slopes for both correlations imply a concerted •OH transfer reaction for both metals. Eyring analyses of reactions for 1 and 3 with (4-X-C₆H₄)₃C• (X = tBu, CN) also give good linear correlations, and a comparison of the resulting activation parameters highlight the importance of entropy in these •OH transfer reactions. Density functional theory calculations of the reaction profiles show a concerted process with one transition state for all radicals investigated and help to explain the electronic features of the OH rebound process.

Lignin Biodegradation by a Cytochrome P450 Enzyme: A Computational Study into Syringol Activation by GcoA

A recently characterized cytochrome P450 isozyme GcoA activates lignin components through a selective O-demethylation or alternatively an acetal formation reaction. These are important reactions in biotechnology and, because lignin is readily available; it being the main component in plant cell walls. In this work we present a density functional theory study on a large active site model of GcoA to investigate syringol activation by an iron(IV)-oxo heme cation radical oxidant (Compound I) leading to hemiacetal and acetal products. Several substrate-binding positions were tested and full energy landscapes calculated. The study shows that substrate positioning determines the product distributions. Thus, with the phenol group pointing away from the heme, an O-demethylation is predicted, whereas an initial hydrogen-atom abstraction of the weak phenolic O-H group would trigger a pathway leading to ring-closure to form acetal products. Predictions on how to engineer P450 GcoA to get more selective product distributions are given.

How Does Replacement of the Axial Histidine Ligand in Cytochrome c Peroxidase by Nδ-Methyl Histidine Affect Its Properties and Functions? A Computational Study

Heme peroxidases have important functions in nature related to the detoxification of H₂O₂. They generally undergo a catalytic cycle where, in the first stage, the iron(III)-heme-H₂O₂ complex is converted into an iron(IV)-oxo-heme cation radical species called Compound I. Cytochrome c peroxidase Compound I has a unique electronic configuration among heme enzymes where a metal-based biradical is coupled to a protein radical on a nearby Trp residue. Recent work using the engineered N_δ-methyl histidine-ligated cytochrome c peroxidase highlighted changes in spectroscopic and catalytic properties upon axial ligand substitution. To understand the axial ligand effect on structure and reactivity of peroxidases and their axially N_δ-methyl histidine engineered forms, we did a computational study. We created active site cluster models of various sizes as mimics of horseradish peroxidase and cytochrome c peroxidase Compound I. Subsequently, we performed density functional theory studies on the structure and reactivity of these complexes with a model substrate (styrene). Thus, the work shows that the N_δ-methyl histidine group has little effect on the electronic configuration and structure of Compound I and little changes in bond lengths and the same orbital occupation is obtained. However, the N_δ-methyl histidine modification impacts electron transfer processes due to a change in the reduction potential and thereby influences reactivity patterns for oxygen atom transfer. As such, the substitution of the axial histidine by N_δ-methyl histidine in peroxidases slows down oxygen atom transfer to substrates and makes Compound I a weaker oxidant. These studies are in line with experimental work on N_δ-methyl histidine-ligated cytochrome c peroxidases and highlight how the hydrogen bonding network in the second coordination sphere has a major impact on the function and properties of the enzyme.

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.

Bioengineering of Cytochrome P450 OleTJE: How Does Substrate Positioning Affect the Product Distributions?

The cytochromes P450 are versatile enzymes found in all forms of life. Most P450s use dioxygen on a heme center to activate substrates, but one class of P450s utilizes hydrogen peroxide instead. Within the class of P450 peroxygenases, the P450 OleT_JE isozyme binds fatty acid substrates and converts them into a range of products through the α-hydroxylation, β-hydroxylation and decarboxylation of the substrate. The latter produces hydrocarbon products and hence can be used as biofuels. The origin of these product distributions is unclear, and, as such, we decided to investigate substrate positioning in the active site and find out what the effect is on the chemoselectivity of the reaction. In this work we present a detailed computational study on the wild-type and engineered structures of P450 OleT_JE using a combination of density functional theory and quantum mechanics/molecular mechanics methods. We initially explore the wild-type structure with a variety of methods and models and show that various substrate activation transition states are close in energy and hence small perturbations as through the protein may affect product distributions. We then engineered the protein by generating an in silico model of the double mutant Asn242Arg/Arg245Asn that moves the position of an active site Arg residue in the substrate-binding pocket that is known to form a salt-bridge with the substrate. The substrate activation by the iron(IV)-oxo heme cation radical species (Compound I) was again studied using quantum mechanics/molecular mechanics (QM/MM) methods. Dramatic differences in reactivity patterns, barrier heights and structure are seen, which shows the importance of correct substrate positioning in the protein and the effect of the second-coordination sphere on the selectivity and activity of enzymes.

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.

Cross-linking of aromatic phenolate groups by cytochrome P450 enzymes: a model for the biosynthesis of vancomycin by OxyB

Aromatic cross-linking by cytochrome P450 enzymes was studied computationally. P450 Compound I rapidly abstracts two weak phenolic H-atoms that link up via a rate-determining C–O bond formation.

Iron- and cobalt-catalyzed C(sp3)–H bond functionalization reactions and their application in organic synthesis

This review highlights the developments in iron and cobalt catalyzed C(sp³)–H bond functionalization reactions with emphasis on their applications in organic synthesis, i.e. natural products and pharmaceuticals synthesis and/or modification.

Properties and reactivity of μ-nitrido-bridged dimetal porphyrinoid complexes: how does ruthenium compare to iron?

Methane hydroxylation by metal-oxo oxidants is one of the Holy Grails in biomimetic and biotechnological chemistry. The only enzymes known to perform this reaction in Nature are iron-containing soluble methane monooxygenase and copper-containing particulate methane monooxygenase. Furthermore, few biomimetic iron-containing oxidants have been designed that can hydroxylate methane efficiently. Recent studies reported that μ-nitrido-bridged diiron(IV)-oxo porphyrin and phthalocyanine complexes hydroxylate methane to methanol efficiently. To find out whether the reaction rates are enhanced by replacing iron by ruthenium, we performed a detailed computational study. Our work shows that the μ-nitrido-bridged diruthenium(IV)-oxo reacts with methane via hydrogen atom abstraction barriers that are considerably lower in energy (by about 5 kcal mol^‒1) as compared to the analogous diiron(IV)-oxo complex. An analysis of the electronic structure implicates similar spin and charge distributions for the diiron(IV)-oxo and diruthenium(IV)-oxo complexes, but the strength of the O‒H bond formed during the reaction is much stronger for the latter. As such a larger hydrogen atom abstraction driving force for the Ru complex than for the Fe complex is found, which should result in higher reactivity in the oxidation of methane.

Mechanism of Oxidative Activation of Fluorinated Aromatic Compounds by N-Bridged Diiron-Phthalocyanine: What Determines the Reactivity?

The biodegradation of compounds with C-F bonds is challenging due to the fact that these bonds are stronger than the C-H bond in methane. In this work, results on the unprecedented reactivity of a biomimetic model complex that contains an N-bridged diiron-phthalocyanine are presented; this model complex is shown to react with perfluorinated arenes under addition of H₂ O₂ effectively. To get mechanistic insight into this unusual reactivity, detailed density functional theory calculations on the mechanism of C₆ F₆ activation by an iron(IV)-oxo active species of the N-bridged diiron phthalocyanine system were performed. Our studies show that the reaction proceeds through a rate-determining electrophilic C-O addition reaction followed by a 1,2-fluoride shift to give the ketone product, which can further rearrange to the phenol. A thermochemical analysis shows that the weakest C-F bond is the aliphatic C-F bond in the ketone intermediate. The oxidative defluorination of perfluoroaromatics is demonstrated to proceed through a completely different mechanism compared to that of aromatic C-H hydroxylation by iron(IV)-oxo intermediates such as cytochrome P450 Compound I.

Electronic Structure of N-Bridged High-Valent Diiron-Oxo

Density functional theory (DFT) and an advanced ab initio technique based on density matrix renormalization group (DMRG-CASPT2) were employed to investigate a reactive N-bridged high-valent diiron-oxo species involved in H-abstraction reactions. We studied in detail two important doublet states, the ground state with two iron(IV) centers and a mixed valence Fe^V -Fe^IV excited state. We found that the latter state is low-lying. Furthermore, its electronic structure and spin density imply that it has significantly higher H-abstraction reactivity than the ground state. This low-lying excited state might be the reason behind the high oxidation reactivity of this diiron-oxo species towards methane.

DFT study on the effect of proximal residues on the Mycobacterium tuberculosis catalase-peroxidase (katG) heme compound I intermediate and its bonding interaction with isoniazid

Isoniazid (INH) is converted into isonicotinyl radical through its interaction with the catalase-peroxidase (katG) enzyme present in the cells of Mycobacterium tuberculosis (M. tb.), the bacteria that causes the tuberculosis disease. This process is important because resistance of M. tb. cells to INH treatment has been associated with the failure of completion of this process. However, this process is poorly understood and there are a variety of conflicting theories about the details of the mechanism of this process. One theory suggests that INH binds to katG and transfers a single electron to the heme while the heme is in its two electron oxidized state, compound I [Fe(iv)Por˙+] (CpdI). In this study, DFT calculations at the UB3LYP/6-31g(d)-LANL2DZ level of theory are used to study the M. tb. katG CpdI molecule. Different models of the M. tb. CpdI molecule were prepared and the calculations revealed the impact of Trp321 on the electronic properties of the heme. Without Trp321 the heme assumed a triradical state with single electrons on the πxy and πyz orbitals of Fe and another on the a2u orbital of the porphyrin ring that can either be coupled with the first two, to assume a quartet state, or decoupled to form a doublet state. With Trp321, however, a transfer of an electron from the πTrp orbital to a2u porphyrin orbital leads to loss of radical character of the porphyrin and leaves the Trp321 group with a radical character. INH was observed to have strong interaction with CpdI and the absence of Trp321 significantly decreased the binding energy by 2 kcal mol-1 explaining the importance of Trp321 in the binding of INH. The results in this study show the importance of Trp321 in the binding of INH and its effect on its electronic properties, which is consistent with previous experimental findings that mutation of Trp321 results in an increase in drug resistance.