Can a Mononuclear Iron(III)‐Superoxo Active Site Catalyze the Decarboxylation of Dodecanoic Acid in UndA to Produce Biofuels?

Effect of monodentate heterocycle co-ligands on the μ-1,2-peroxo-diiron(III) mediated aldehyde deformylation reactions

Peroxo-diiron(III) species are present in the active sites of many metalloenzymes that carry out challenging organic transformations. The reactivity of these species is influenced by various factors, such as the structure and topology of the supporting ligands, the identity of the axial and equatorial co-ligands, and the oxidation states of the metal ion(s). In this study, we aim to diversify the importance of equatorial ligands in controlling the reactivity of peroxo-diiron(III) species. As a model compound, we chose the previously published and fully characterized [(PBI)2(CH3CN)FeIII(μ-O2)FeIII(CH3CN)(PBI)2]4+ complex, where the steric effect of the four PBI ligands is minimal, so the labile CH3CN molecules easily can be replaced by different monodentate co-ligands (substituted pyridines and N-donor heterocyclic compounds). Thus, their effect on the electronic and spectral properties of peroxo-divas(III) intermediates could be easily investigated. The relationship between structure and reactivity was also investigated in the stoichiometric deformylation of PPA mediated by peroxo-diiron(III) complexes. It was found that the deformylation rates are influenced by the Lewis acidity and redox properties of the metal centers, and showed a linear correlation with the FeIII/FeII redox potentials (in the range of 197 to 415 mV).

DFT mechanistic study of biomimetic diiron complex catalyzed dehydrogenation: Unexpected Fe(III)Fe(III)-1,1-μ-hydroperoxy active species for hydride abstraction

The diiron active site is pivotal in catalyzing transformations in both biological and chemical systems. Recently, a range of biomimetic diiron catalysts have been synthesized, drawing inspiration from the active architecture of soluble methane monooxygenase (sMMO). These catalysts have been successfully deployed for the dehydrogenation of indolines, marking a significant advancement in the field. Using density functional theory (DFT) calculations, we have identified a novel mechanistic pathway that governs the dehydrogenation of indolines catalyzed by a biomimetic diiron complex. Specifically, this reaction is facilitated by the transfer of a hybrid atom from the C1 position of the substrate to the distal oxygen atom of the Fe(III)Fe(III)-1,1-μ-hydroperoxy active species. This transfer serves as the rate-limiting step for the heterolytic cleavage of the OO bond, ultimately generating the substrate cation. The mechanism we propose aligns well with mechanistic investigations incorporating both kinetic isotope effect (KIE) measurements and evaluations of stereochemical selectivity. This research contributes to the broader scientific understanding of catalysis involving biomimetic diiron complexes and offers valuable insights into the catalytic behaviors of non-heme diiron metalloenzymes.

Unmasking the Iron-Oxo Bond of the [(Ligand)Fe-OIAr]2+/+ Complexes

ArIO (ArI = 2-(^tBuSO₂)C₆H₄I) is an oxidant used to oxidize Fe^II species to their Fe^IV-oxo state, enabling hydrogen-atom transfer (HAT) and oxygen-atom transfer (OAT) reactions at low energy barriers. ArIO, as a ligand, generates masked Feⁿ═O species of the type Fe^(n-2)-OIAr. Herein, we used gas-phase ion-molecule reactions and DFT calculations to explore the properties of masked iron-oxo species and to understand their unmasking mechanisms. The theory shows that the I-O bond cleavage in [(TPA)Fe^IVO(ArIO)]²⁺ (1²⁺, TPA = tris(2-pyridylmethyl)amine)) is highly endothermic; therefore, it can be achieved only in collision-induced dissociation of 1²⁺ leading to the unmasked iron(VI) dioxo complex. The reduction of 1²⁺ by HAT leads to [(TPA)Fe^IIIOH(ArIO)]²⁺ with a reduced energy demand for the I-O bond cleavage but is, however, still endothermic. The exothermic unmasking of the Fe═O bond is predicted after one-electron reduction of 1²⁺ or after OAT reactivity. The latter leads to the 4e^- oxidation of unsaturated hydrocarbons: The initial OAT from [(TPA)Fe^IVO(ArIO)]²⁺ leads to the epoxidation of an alkene and triggers the unmasking of the second Fe═O bond still within one collisional complex. The second oxidation step starts with HAT from a C-H bond and follows with the rebound of the C-radical and the OH group. The process starting with the one-electron reduction could be studied with [(TQA)Fe^IVO(ArIO)]²⁺ (2²⁺, TQA = tris(2-quinolylmethyl)amine)) because it has a sufficient electron affinity for electron transfer with alkenes. Accordingly, the reaction of 2²⁺ with 2-carene leads to [(TQA)Fe^IIIO(ArIO)]²⁺ that exothermically eliminates ArI and unmasks the reactive Fe^V-dioxo species.

What Drives Radical Halogenation versus Hydroxylation in Mononuclear Nonheme Iron Complexes? A Combined Experimental and Computational Study

Nonheme iron halogenases are unique enzymes in nature that selectively activate an aliphatic C–H bond of a substrate to convert it into C–X (X = Cl/Br, but not F/I). It is proposed that they generate an Fe^III(OH)(X) intermediate in their catalytic cycle. The analogous Fe^III(OH) intermediate in nonheme iron hydroxylases transfers OH^• to give alcohol product, whereas the halogenases transfer X^• to the carbon radical substrate. There remains significant debate regarding what factors control their remarkable selectivity of the halogenases. The reactivity of the complexes Fe^III(BNPA^Ph2O)(OH)(X) (X = Cl, Br) with a secondary carbon radical (R^•) is described. It is found that X^• transfer occurs with a secondary carbon radical, as opposed to OH^• transfer with tertiary radicals. Comprehensive computational studies involving density functional theory were carried out to examine the possible origins of this selectivity. The calculations reproduce the experimental findings, which indicate that halogen transfer is not observed for the tertiary radicals because of a nonproductive equilibrium that results from the endergonic nature of these reactions, despite a potentially lower reaction barrier for the halogenation pathway. In contrast, halogen transfer is favored for secondary carbon radicals, for which the halogenated product complex is thermodynamically more stable than the reactant complex. These results are rationalized by considering the relative strengths of the C–X bonds that are formed for tertiary versus secondary carbon centers. The computational analysis also shows that the reaction barrier for halogen transfer is significantly dependent on secondary coordination sphere effects, including steric and H-bonding interactions.

Which is the real oxidant in the competitive ligand self-hydroxylation and substrate oxidation, a biomimetic iron(II)-hydroperoxo species or an oxo-iron(IV)-hydroxy one?

Nonheme iron(II)-hydroperoxo species (FeII-(η2-OOH)) 1 and the concomitant oxo-iron(IV)-hydroxyl one 2 are proposed as the key intermediates of a large class of 2-oxoglutarate dependent dioxygenases (e.g., isopenicillin N synthase). Extensive...

Decoding Regioselective Reaction Mechanism of the Gentisic Acid Catalyzed by Gentisate 1,2-Dioxygenase Enzyme

Gentisate 1,2-dioxygenase (GDO), a ring-fission non-heme dioxygenase enzyme, displays a unique regioselective reaction of gentisic acid (GTQ) in the presence of molecular oxygen. GTQ is an important intermediate in the...

Electrostatic perturbations in the substrate-binding pocket of taurine/α-ketoglutarate dioxygenase determine its selectivity

Taurine/α‐ketoglutarate dioxygenase is an important enzyme that takes part in the cysteine catabolism process in the human body and selectively hydroxylates taurine at the C¹‐position. Recent computational studies showed that in the gas‐phase the C²−H bond of taurine is substantially weaker than the C¹−H bond, yet no evidence exists of 2‐hydroxytaurine products. To this end, a detailed computational study on the selectivity patterns in TauD was performed. The calculations show that the second‐coordination sphere and the protonation states of residues play a major role in guiding the enzyme to the right selectivity. Specifically, a single proton on an active site histidine residue can change the regioselectivity of the reaction through its electrostatic perturbations in the active site and effectively changes the C¹−H and C²−H bond strengths of taurine. This is further emphasized by many polar and hydrogen bonding interactions of the protein cage in TauD with the substrate and the oxidant that weaken the pro‐R C¹−H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently.

BesC Initiates C-C Cleavage through a Substrate-Triggered and Reactive Diferric-Peroxo Intermediate

BesC catalyzes the iron- and O₂-dependent cleavage of 4-chloro-l-lysine to form 4-chloro-l-allylglycine, formaldehyde, and ammonia. This process is a critical step for a biosynthetic pathway that generates a terminal alkyne amino acid which can be leveraged as a useful bio-orthogonal handle for protein labeling. As a member of an emerging family of diiron enzymes that are typified by their heme oxygenase-like fold and a very similar set of coordinating ligands, recently termed HDOs, BesC performs an unusual type of carbon-carbon cleavage reaction that is a significant departure from reactions catalyzed by canonical dinuclear-iron enzymes. Here, we show that BesC activates O₂ in a substrate-gated manner to generate a diferric-peroxo intermediate. Examination of the reactivity of the peroxo intermediate with a series of lysine derivatives demonstrates that BesC initiates this unique reaction trajectory via cleavage of the C4-H bond; this process represents the rate-limiting step in a single turnover reaction. The observed reactivity of BesC represents the first example of a dinuclear-iron enzyme that utilizes a diferric-peroxo intermediate to capably cleave a C-H bond as part of its native function, thus circumventing the formation of a high-valent intermediate more commonly associated with substrate monooxygenations.

Biodegradation of herbicides by a plant nonheme iron dioxygenase: mechanism and selectivity of substrate analogues

Aryloxyalkanoate dioxygenases are unique herbicide biodegrading nonheme iron enzymes found in plants and hence, from environmental and agricultural point of view they are important and valuable. However, they often are substrate specific and little is known on the details of the mechanism and the substrate scope. To this end, we created enzyme models and calculate the mechanism for 2,4-dichlorophenoxyacetic acid biodegradation and 2-methyl substituted analogs by density functional theory. The work shows that the substrate binding is tight and positions the aliphatic group close to the metal center to enable a chemoselective reaction mechanism to form the C 2 -hydroxy products, whereas the aromatic hydroxylation barriers are well higher in energy. Subsequently, we investigated the metabolism of R - and S -methyl substituted inhibitors and show that these do not react as efficiently as 2,4-dichlorophenoxyacetic acid substrate due to stereochemical clashes in the active site and particularly for the R -isomer give high rebound barriers.

Structure and Functional Differences of Cysteine and 3-Mercaptopropionate Dioxygenases: A Computational Study

Thiol dioxygenases are important enzymes for human health; they are involved in the detoxification and catabolism of toxic thiol-containing natural products such as cysteine. As such, these enzymes have relevance to the development of Alzheimer's and Parkinson's diseases in the brain. Recent crystal structure coordinates of cysteine and 3-mercaptopropionate dioxygenase (CDO and MDO) showed major differences in the second-coordination spheres of the two enzymes. To understand the difference in activity between these two analogous enzymes, we created large, active-site cluster models. We show that CDO and MDO have different iron(III)-superoxo-bound structures due to differences in ligand coordination. Furthermore, our studies show that the differences in the second-coordination sphere and particularly the position of a positively charged Arg residue results in changes in substrate positioning, mobility and enzymatic turnover. Furthermore, the substrate scope of MDO is explored with cysteinate and 2-mercaptosuccinic acid and their reactivity is predicted.

Density Functional Theory Study into the Reaction Mechanism of Isonitrile Biosynthesis by the Nonheme Iron Enzyme ScoE

The nonheme iron enzyme ScoE catalyzes the biosynthesis of an isonitrile substituent in a peptide chain. To understand details of the reaction mechanism we created a large active site cluster model of 212 atoms that contains substrate, the active oxidant and the first- and second-coordination sphere of the protein and solvent. Several possible reaction mechanisms were tested and it is shown that isonitrile can only be formed through two consecutive catalytic cycles that both use one molecule of dioxygen and α-ketoglutarate. In both cycles the active species is an iron(IV)-oxo species that in the first reaction cycle reacts through two consecutive hydrogen atom abstraction steps: first from the N–H group and thereafter from the C–H group to desaturate the NH-CH2 bond. The alternative ordering of hydrogen atom abstraction steps was also tested but found to be higher in energy. Moreover, the electronic configurations along that pathway implicate an initial hydride transfer followed by proton transfer. We highlight an active site Lys residue that is shown to donate charge in the transition states and influences the relative barrier heights and bifurcation pathways. A second catalytic cycle of the reaction of iron(IV)-oxo with desaturated substrate starts with hydrogen atom abstraction followed by decarboxylation to give isonitrile directly. The catalytic cycle is completed with a proton transfer to iron(II)-hydroxo to generate the iron(II)-water resting state. The work is compared with experimental observation and previous computational studies on this system and put in a larger perspective of nonheme iron chemistry.

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.

What Determines the Selectivity of Arginine Dihydroxylation by the Nonheme Iron Enzyme OrfP?

The nonheme iron enzyme OrfP reacts with l-Arg selectively to form the 3R,4R-dihydroxyarginine product, which in mammals can inhibit the nitric oxide synthase enzymes involved in blood pressure control. To understand the mechanisms of dioxygen activation of l-Arg by OrfP and how it enables two sequential oxidation cycles on the same substrate, we performed a density functional theory study on a large active site cluster model. We show that substrate binding and positioning in the active site guides a highly selective reaction through C³ -H hydrogen atom abstraction. This happens despite the fact that the C³ -H and C⁴ -H bond strengths of l-Arg are very similar. Electronic differences in the two hydrogen atom abstraction pathways drive the reaction with an initial C³ -H activation to a low-energy ⁵ σ-pathway, while substrate positioning destabilizes the C⁴ -H abstraction and sends it over the higher-lying ⁵ π-pathway. We show that substrate and monohydroxylated products are strongly bound in the substrate binding pocket and hence product release is difficult and consequently its lifetime will be long enough to trigger a second oxygenation cycle.

Computational studies of DNA base repair mechanisms by nonheme iron dioxygenases: selective epoxidation and hydroxylation pathways

DNA base repair mechanisms of alkylated DNA bases is an important reaction in chemical biology and particularly in the human body. It is typically catalyzed by an α-ketoglutarate-dependent nonheme iron dioxygenase named the AlkB repair enzyme. In this work we report a detailed computational study into the structure and reactivity of AlkB repair enzymes with alkylated DNA bases. In particular, we investigate the aliphatic hydroxylation and C[double bond, length as m-dash]C epoxidation mechanisms of alkylated DNA bases by a high-valent iron(iv)-oxo intermediate. Our computational studies use quantum mechanics/molecular mechanics methods on full enzymatic structures as well as cluster models on active site systems. The work shows that the iron(iv)-oxo species is rapidly formed after dioxygen binding to an iron(ii) center and passes a bicyclic ring structure as intermediate. Subsequent cluster models explore the mechanism of substrate hydroxylation and epoxidation of alkylated DNA bases. The work shows low energy barriers for substrate activation and consequently energetically feasible pathways are predicted. Overall, the work shows that a high-valent iron(iv)-oxo species can efficiently dealkylate alkylated DNA bases and return them into their original form.

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.

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.

Mechanism of fatty acid decarboxylation catalyzed by a non-heme iron oxidase (UndA): a QM/MM study

UndA is a non-heme iron enzyme that was recognized to catalyze the decarboxylation of medium chain (C10-C14) fatty acids to produce trace amounts of 1-alkenes. Owing to the electron imbalance during the oxidative decarboxylation of the substrate and the reduction of O₂, only single turnover reactions were obtained in UndA in vitro assays. Unlike the general non-heme iron enzymes, the catalytic efficiency of UndA is quite low. According to the previous proposal, both Fe^III-OO˙^- and Fe^IV[double bond, length as m-dash]O complexes may abstract the β-H of fatty acids to trigger the oxidative decarboxylation reaction. Herein, on the basis of the crystal structures of UndA in complex with the substrate analogues, we constructed a series of computational models and performed quantum mechanics/molecular mechanics (QM/MM) calculations to explore the UndA-catalyzed decarboxylation using lauric acid as the substrate. Our calculation results reveal that only the Fe^III-OO˙^- complex can initiate the decarboxylation, and the substrate (lauric acid) should monodentately coordinate to the Fe center to facilitate the β-H abstraction. In addition, the monodentate coordination corresponds to higher relative energy than the bidentate mode, which may explain the low efficiency of UndA. It is also revealed that as long as the β-H is extracted by the Fe^III-OO˙^-, the decarboxylation of the substrate radical is quite easy, and an electron transfer from the substrate to the iron center is the prerequisite. For the Fe^IV[double bond, length as m-dash]O complex, since the β-H is far from the O_Fe atom and the angle of ∠Fe-O-H is 53.1°, the H-abstraction is calculated to be difficult.