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

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.

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

CO2 Reduction by an Iron(I) Porphyrinate System: Effect of Hydrogen Bonding on the Second Coordination Sphere

Transforming CO2 into valuable materials is an important reaction in catalysis, especially because CO2 concentrations in the atmosphere have been growing steadily due to extensive fossil fuel usage. From an environmental perspective, reduction of CO2 to valuable materials should be catalyzed by an environmentally benign catalyst and avoid the use of heavy transition-metal ions. In this work, we present a computational study into a novel iron(I) porphyrin catalyst for CO2 reduction, namely, with a tetraphenylporphyrin ligand and analogues. In particular, we investigated iron(I) tetraphenylporphyrin with one of the meso-phenyl groups substituted with o-urea, p-urea, or o-2-amide groups. These substituents can provide hydrogen-bonding interactions in the second coordination sphere with bound ligands and assist with proton relay. Furthermore, our studies investigated bicarbonate and phenol as stabilizers and proton donors in the reaction mechanism. Potential energy landscapes for double protonation of iron(I) porphyrinate with bound CO2 are reported. The work shows that the bicarbonate bridges the urea/amide groups to the CO2 and iron center and provides a tight bonding pattern with strong hydrogen-bonding interactions that facilitates easy proton delivery and reduction of CO2. Specifically, bicarbonate provides a low-energy proton shuttle mechanism to form CO and water efficiently. Furthermore, the o-urea group locks bicarbonate and CO2 in a tight orientation and helps with ideal proton transfer, while there is more mobility and lesser stability with an o-amide group in that position instead. Our calculations show that the o-urea group leads to reduction in proton-transfer barriers, in line with experimental observation. We then applied electric-field-effect calculations to estimate the environmental effects on the two proton-transfer steps in the reaction. These calculations describe the perturbations that enhance the driving forces for the proton-transfer steps and have been used to make predictions about how the catalysts can be further engineered for more enhanced CO2 reduction processes.

Mechanism of CO2 Reduction to Methanol with H2 on an Iron(II)-scorpionate Catalyst

CO2 utilization is an important process in the chemical industry with great environmental power. In this work we show how CO2 and H2 can be reacted to form methanol on an iron(II) center and highlight the bottlenecks for the reaction and what structural features of the catalyst are essential for efficient turnover. The calculations predict the reactions to proceed through three successive reaction cycles that start with heterolytic cleavage of H2 followed by sequential hydride and proton transfer processes. The H2 splitting process is an endergonic process and hence high pressures will be needed to overcome this step and trigger the hydrogenation reaction. Moreover, H2 cleavage into a hydride and proton requires a metal to bind hydride and a nearby source to bind the proton, such as an amide or pyrazolyl group, which the scorpionate ligand used here facilitates. As such the computations highlight the non-innocence of the ligand scaffold through proton shuttle from H2 to substrate as an important step in the reaction mechanism.

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.

Mechanism of Nitrogen Reduction to Ammonia in a Diiron Model of Nitrogenase

Nitrogenase is a fascinating enzyme in biology that reduces dinitrogen from air to ammonia through stepwise reduction and protonation. Despite it being studied in detail by experimental and computational groups, there are still many unknown factors in the catalytic cycle of nitrogenase, especially related to the addition of protons and electrons and their order. A recent biomimetic study characterized a potential dinitrogen-bridged diiron cluster as a synthetic model of nitrogenase. Using strong acid and reductants, the dinitrogen was converted into ammonia molecules, but details of the mechanism remains unknown. In particular, it was unclear from the experimental studies whether the proton and electron transfer steps are sequential or alternating. Moreover, the work failed to establish what the function of the diiron core is and whether it split into mononuclear iron fragments during the reaction. To understand the structure and reactivity of the biomimetic dinitrogen-bridged diiron complex [(P2P'PhFeH)2(μ-N2)] with triphenylphosphine ligands, we performed a density functional theory study. Our computational methods were validated against experimental crystal structure coordinates, Mössbauer parameters, and vibrational frequencies and show excellent agreement. Subsequently, we investigated the alternating and consecutive addition of electrons and protons to the system. The calculations identify a number of possible reaction channels, namely, same-site protonation, alternating protonation, and complex dissociation into mononuclear iron centers. The calculations show that the overall mechanism is not a pure sequential set of electron and proton transfers but a mixture of alternating and consecutive steps. In particular, the first reaction steps will start with double proton transfer followed by an electron transfer, while thereafter, there is another proton transfer and a second electron transfer to give a complex whereby ammonia can split off with a low energetic barrier. The second channel starts with alternating protonation of the two nitrogen atoms, whereafter the initial double proton transfer, electrons and protons are added sequentially to form a hydrazine-bound complex. The latter split off ammonia spontaneously after further protonation. The various reaction channels are analyzed with valence bond and orbital diagrams. We anticipate the nitrogenase enzyme to operate with mixed alternating and consecutive protonation and electron transfer steps.

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.

Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study

Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho- or para-position. The barriers of para-rebound are well lower in energy than those for ortho-rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C⁴-aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.

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.

Charge-Transfer Spectroscopy of Bisaxially Coordinated Iron(II) Phthalocyanines through the Prism of the Lever's EL Parameters Scale, MCD Spectroscopy, and TDDFT Calculations

The position of the experimentally observed (in the UV-vis and magnetic circular dichroism (MCD) spectra) low-energy metal-to-ligand charge-transfer (MLCT) band in low-spin iron(II) phthalocyanine complexes of general formula PcFeL₂, PcFeL'L″, and [PcFeX₂]^2- (L, L', or L″ are neutral and X^- is an anionic axial ligand) was correlated with the Lever's electrochemical E_L scale values for the axial ligands. The time-dependent density functional theory (TDDFT)-predicted UV-vis spectra are in very good agreement with the experimental data for all complexes. In the majority of compounds, TDDFT predicts that the first degenerate MLCT band that correlates with the MCD A-term observed between 360 and 480 nm is dominated by an e_g (Fe, d_π) → b_1u (Pc, π*) single-electron excitation (in traditional D_4h point group notation) and agrees well with the previous assignment discussed by Stillman and co-workers[ Inorg. Chem. 1994, 33, 573-583]. The TDDFT calculations also suggest a small energy gap for b_1u/b_2u (Pc, π*) orbital splitting and closeness of the MLCT₁ e_g (Fe, d_π) → b_1u (Pc, π*) and MLCT₂ e_g (Fe, d_π) → b_2u (Pc, π*) transitions. In the case of the PcFeL₂ complexes with phosphines as the axial ligands, additional degenerate charge-transfer transitions were observed between 450 and 500 nm. These transitions are dominated by a_2u (Pc + L, π) → e_g (Pc, π*) single-electron excitations and are unique for the PcFe(PR₃)₂ complexes. The energy of the phthalocyanine-based a_2u orbital has large axial ligand dependency and is the reason for a large energy deviation for B1 a_2u (Pc + L, π) → e_g (Pc, π*) transition. The energies of the axial ligand-to-iron, axial ligand-to-phthalocyanine, iron-to-axial ligand, and phthalocyanine-to-axial ligand charge-transfer transitions were discussed on the basis of TDDFT calculations.

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.

Oxidative dehalogenation of halophenols by high-valent nonheme iron(IV)-oxo intermediates

Mononuclear high-valent iron(IV)-oxo intermediates are excellent oxidants towards oxygenation reactions by heme and nonheme metalloenzymes and their model systems. One of the most important functions of these intermediates in nature is to detoxify various environmental pollutants. Organic substrates, such as halogenated phenols, are known to be water pollutants which can be degraded to their less hazardous forms through an oxidation reaction by iron(IV)-oxo complexes. Metalloproteins in nature utilize various types of second-coordination sphere interactions to anchor the substrate in the vicinity of the active site. This concept of substrate-binding is well-known for natural enzymes, but is elusive for the relevant biomimetic model systems. Herein, we report the oxidative reactivity patterns of an iron(IV)-oxo intermediate, [Fe^IV(O)(2PyN2Q)]²⁺, (2PyN2Q = 1,1-di(pyridin-2yl)-N,N-bis(quinolin-2-ylmethyl)methanamine) with a series of mono-, di- and tri-halophenols. A detailed experimental study shows that the dehalogenation reactions of the halophenols by such iron(IV)-oxo intermediates proceed via an initial hydrogen atom abstraction from the phenolic O-H group. Furthermore, based on the size and nucleophilicity of the halophenol, an intermediate substrate-bound species forms that is a phenolate adduct to the ferric species, which thereafter leads to the formation of the corresponding products.

Application of Lever's EL Parameter Scale toward Fe(II)/Fe(III) versus Pc(2-)/Pc(1-) Oxidation Process Crossover Point in Axially Coordinated Iron(II) Phthalocyanine Complexes

The electronic structures and, particularly, the nature of the HOMO in a series of PcFeL₂, PcFeL'L″, and [PcFeX₂]^2- complexes (Pc = phthalocyaninato(2-) ligand; L = NH₃, n-BuNH₂, imidazole (Im), pyridine (Py), PMe₃, PBu₃, t-BuNC, P(OBu)₃, and DMSO; L' = CO; L″ = NH₃ or n-BuNH₂; X = NCO^-, NCS^-, CN^-, imidazolate (Im^-), or 1,2,4-triazolate(Tz^-)) were probed by electrochemical, spectroelectrochemical, and chemical oxidation as well as theoretical (density functional theory, DFT) studies. In general, energies of the metal-centered occupied orbitals in various six-coordinate iron phthalocyanine complexes correlate well with Lever Electrochemical Parameter E_L and intercross the phthalocyanine-centered a_1u orbital in several compounds with moderate-to-strong π-accepting axial ligands. In these cases, an oxidation of the phthalocyanine macrocycle (Pc(2-)/Pc(1-)) rather than the central metal ion (Fe(II)/Fe(III)) was theoretically predicted and experimentally confirmed.

pH Changes That Induce an Axial Ligand Effect on Nonheme Iron(IV) Oxo Complexes with an Appended Aminopropyl Functionality

Nonheme iron enzymes often utilize a high-valent iron(IV) oxo species for the biosynthesis of natural products, but their high reactivity often precludes structural and functional studies of these complexes. In this work, a combined experimental and computational study is presented on a biomimetic nonheme iron(IV) oxo complex bearing an aminopyridine macrocyclic ligand and its reactivity toward olefin epoxidation upon changes in the identity and coordination ability of the axial ligand. Herein, we show a dramatic effect of the pH on the oxygen-atom-transfer (OAT) reaction with substrates. In particular, these changes have occurred because of protonation of the axial-bound pendant amine group, where its coordination to iron is replaced by a solvent molecule or anionic ligand. This axial ligand effect influences the catalysis, and we observe enhanced cyclooctene epoxidation yields and turnover numbers in the presence of the unbound protonated pendant amine group. Density functional theory studies were performed to support the experiments and highlight that replacement of the pendant amine with a neutral or anionic ligand dramatically lowers the rate-determining barriers of cyclooctene epoxidation. The computational work further establishes that the change in OAT is due to electrostatic interactions of the pendant amine cation that favorably affect the barrier heights.

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.

Formation Mechanism of Cofactor Cys-Tyr in the Cysteine Dioxygenases (CDO and F2-CDO) and Its Influence on Catalysis: A QM/MM Study

Cysteine dioxygenase (CDO) is a nonheme mononuclear iron enzyme, which catalyzes the oxidation of cysteine to cysteine sulfinic acid. Crystal structure studies of mammalian CDO showed that there is a cross-linked cysteine-tyrosine (Cys-Tyr) cofactor in its active site. Moreover, the formation of the Cys-Tyr cofactor requires the metal cofactor (Fe²⁺) and O₂, and it was previously considered to substantially enhance the catalytic efficiency and half-life of CDO. Recently, a pure human CDO (F₂-CDO) without including the Cys-Tyr cofactor was crystalized by the site-directed mutagenesis approach in the anaerobic condition. In this work, to gain insights into the formation mechanism of the Cys-Tyr cofactor and whether it can really promote the catalytic reactivity of CDO, a series of computational models have been constructed, and quantum mechanical/molecular mechanical (QM/MM) calculations have been performed. Our calculation results reveal that WT-CDO and F₂-CDO follow different mechanisms for the formation of the Cys-Tyr cofactor. In F₂-CDO, the cofactor formation contains the H-abstraction, C-S bond formation, intramolecular F migration, and aromatization of the residue F₂Y157, in which the Fe-coordinate dioxygen can be recovered after the formation cofactor; however, in the WT-CDO, the cofactor formation shows some differences. During the reaction, hydrogen peroxide is generated, and the final aromatization requires the assistance of one water molecule. Furthermore, the overall barriers of cofactor formation are always higher than l-cysteine oxidation for both WT-CDO and F₂-CDO irrespective of the absence or presence of the cofactor. Thus, we can theoretically confirm that the Cys-Tyr cofactor is not essential for the oxidation activity of CDO, and cofactor formation is just an accompanying reaction but not a prerequisite for the oxidation reaction. These results may provide useful information for understanding the catalysis of CDO.

Oxidative degradation of toxic organic pollutants by water soluble nonheme iron(iv)-oxo complexes of polydentate nitrogen donor ligands

The ability of four mononuclear nonheme iron(iv)-oxo complexes supported by polydentate nitrogen donor ligands to degrade organic pollutants has been investigated. The water soluble iron(ii) complexes upon treatment with ceric ammonium nitrate (CAN) in aqueous solution are converted into the corresponding iron(iv)-oxo complexes. The hydrogen atom transfer (HAT) ability of iron(iv)-oxo species has been exploited for the oxidation of halogenated phenols and other toxic pollutants with weak X-H (X = C, O, S, etc.) bonds. The iron-oxo oxidants can oxidize chloro- and fluorophenols with moderate to high yields under stoichiometric as well as catalytic conditions. Furthermore, these oxidants perform selective oxidative degradation of several persistent organic pollutants (POPs) such as bisphenol A, nonylphenol, 2,4-D (2,4-dichlorophenoxyacetic acid) and gammaxene. This work demonstrates the utility of water soluble iron(iv)-oxo complexes as potential catalysts for the oxidative degradation of a wide range of toxic pollutants, and these oxidants could be considered as an alternative to conventional oxidation methods.

Accurate Prediction of Mössbauer Hyperfine Parameters in Bis-Axially Coordinated Iron(II) Phthalocyanines Using Density Functional Theory Calculations: A Story of a Single Orbital Revealed by Natural Bond Orbital Analysis

Density Functional Theory (DFT) calculations coupled with several exchange-correlation functionals were used for the prediction of Mössbauer hyperfine parameters of 36 bis-axially coordinated iron(II) phthalocyanine complexes with the general formulas PcFeL₂, PcFeL'L″, and [PcFeX₂]^2-, including four new compounds. Both gas-phase and PCM calculations using BPW91 and MN12L exchange-correlation functionals were found to accurately predict both Mössbauer quadrupole splittings and the correct trends in experimentally observed isomer shifts. In comparison, hybrid exchange-correlation functionals underestimated quadrupole splittings, while still accurately predicted isomer shifts. Out of ∼40 exchange-correlation functionals tested, only MN12L was found to correctly reproduce quadrupole splitting trends in the PcFeL₂ complexes coordinated with phosphorus-donor axial ligands (i.e., P(OnBu)₃ ≈ P(OEt)₃ < PMe₃ < P[(CH₂O)₂CH₂]-p-C₆H₄NO₂ < PEt₃ ≈ PnBu₃). Natural Bond Orbital (NBO) analysis was successfully used to explain the general trends in the observed quadrupole splitting for all compounds of interest. In particular, the general trends in the quadrupole splitting correlate well with the axial ligand dependent, NBO-predicted population of the 3d_z2 orbital of the Fe ion and are reflective of the hypothesis proposed by Ohya and co-workers ( Inorg. Chem., 1984, 23, 1303) on the adaptability of the phthalocyanine's π-system toward Fe-L_ax interactions. The first X-ray crystal structure of a PcFeL₂ complex with axial phosphine ligands is also reported.

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.

Metallophthalocyanines as Catalysts in Aerobic Oxidation

The first remarkable property associated to metallophthalocyanines (MPcs) was their chemical “inertness”, which made and make them very attractive as stable and durable industrial dyes. Nevertheless, their rich redox chemistry was also explored in the last decades, making available a solid and detailed knowledge background for further studies on the suitability of MPcs as redox catalysts. An overlook of MPcs and their catalytic activity with dioxygen as oxidants will be discussed here with a special emphasis on the last decade. The mini-review begins with a short introduction to phthalocyanines, from their structure to their main features, going then through the redox chemistry of metallophthalocyanines and their catalytic activity in aerobic oxidation reactions. The most significant systems described in the literature comprise the oxidation of organosulfur compounds such as thiols and thiophenes, the functionalization of alkyl arenes, alcohols, olefins, among other substrates.

Mechanistic Insights into the ortho-Defluorination-Hydroxylation of 2-Halophenolates Promoted by a Bis(μ-oxo)dicopper(III) Complex

C-F bonds are one of the most inert functionalities. Nevertheless, some [Cu₂O₂]²⁺ species are able to defluorinate-hydroxylate ortho-fluorophenolates in a chemoselective manner over other ortho-halophenolates. Albeit it is known that such reactivity is promoted by an electrophilic attack of a [Cu₂O₂]²⁺ core over the arene ring, the crucial details of the mechanism that explain the chemo- and regioselectivity of the reaction remain unknown, and it has not being determined either if Cu^II₂(η²:η²-O₂) or Cu^III₂(μ-O)₂ species are responsible for the initial attack on the arene. Herein, we present a combined theoretical and experimental mechanistic study to unravel the origin of the chemoselectivity of the ortho-defluorination-hydroxylation of 2-halophenolates by the [Cu₂(O)₂(DBED)₂]²⁺ complex (DBED = N,N'-di-tert-butylethylenediamine). Our results show that the equilibria between (side-on)peroxo (P) and bis(μ-oxo) (O) isomers plays a key role in the mechanism, with the latter being the reactive species. Furthermore, on the basis of quantum-mechanical calculations, we were able to rationalize the chemoselective preference of the [Cu₂(O)₂(DBED)₂]²⁺ catalyst for the C-F activation over C-Cl and C-H activations.

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.

Carbon-fluorine bond cleavage mediated by metalloenzymes

Fluorochemicals are a widely distributed class of compounds and have been utilized across a wide range of industries for decades. Given the environmental toxicity and adverse health threats of some fluorochemicals, the development of new methods for their decomposition is significant to public health. However, the carbon-fluorine (C-F) bond is among the most chemically robust bonds; consequently, the degradation of fluorinated hydrocarbons is exceptionally difficult. Here, metalloenzymes that catalyze the cleavage of this chemically challenging bond are reviewed. These enzymes include histidine-ligated heme-dependent dehaloperoxidase and tyrosine hydroxylase, thiolate-ligated heme-dependent cytochrome P450, and four nonheme oxygenases, namely, tetrahydrobiopterin-dependent aromatic amino acid hydroxylase, 2-oxoglutarate-dependent hydroxylase, Rieske dioxygenase, and thiol dioxygenase. While much of the literature regarding the aforementioned enzymes highlights their ability to catalyze C-H bond activation and functionalization, in many cases, the C-F bond cleavage has been shown to occur on fluorinated substrates. A copper-dependent laccase-mediated system representing an unnatural radical defluorination approach is also described. Detailed discussions on the structure-function relationships and catalytic mechanisms provide insights into biocatalytic defluorination, which may inspire drug design considerations and environmental remediation of halogenated contaminants.

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.

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.

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.