Effect of Porphyrin Ligands on the Regioselective Dehydrogenation versus Epoxidation of Olefins by Oxoiron(IV) Mimics of Cytochrome P450

What is the Origin of the Regioselective C3-Hydroxylation of L-Arg by the Nonheme Iron Enzyme Capreomycin C?

The nonheme iron dioxygenase capreomycin C (CmnC) hydroxylates a free L-arginine amino acid regio- and stereospecifically at the C3-position as part of the capreomycin antibiotics biosynthesis. Little is known on its structure, catalytic cycle and substrate specificity and, therefore, a comprehensive computational study was performed. A large QM cluster model of CmnC was created of 297 atoms and the mechanisms for C3-H, C4-H and C5-H hydroxylation and C3-C4 desaturation were investigated. All low-energy pathways correspond to radical reaction mechanisms with an initial hydrogen atom abstraction followed by OH rebound to form alcohol product complexes. The work is compared to alternative L-Arg hydroxylating nonheme iron dioxygenases and the differences in active site polarity are compared. We show that a tight hydrogen bonding network in the substrate binding pocket positions the substrate in an ideal orientation for C3-H activation, whereby the polar groups in the substrate binding pocket induce an electric field effect that guides the selectivity.

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.

Computational study on the endocrine-disrupting metabolic activation of Benzophenone-3 catalyzed by cytochrome P450 1A1: A QM/MM approach

Predicting the metabolic activation mechanism and potential hazardous metabolites of environmental endocrine-disruptors is a challenging and significant task in risk assessment. Here the metabolic activation mechanism of benzophenone-3 catalyzed by P450 1A1 was investigated by using Molecular Dynamics, Quantum Mechanics/Molecular Mechanics and Density Functional Theory approaches. Two elementary reactions involved in the metabolic activation of BP-3 with P450 1A1: electrophilic addition and hydrogen abstraction reactions were both discussed. Further conversion reactions of epoxidation products, ketone products and the formaldehyde formation reaction were investigated in the non-enzymatic environment based on previous experimental reports. Binding affinities analysis of benzophenone-3 and its metabolites to sex hormone binding globulin indirectly demonstrates that they all exhibit endocrine-disrupting property. Toxic analysis shows that the eco-toxicity and bioaccumulation values of the benzophenone-3 metabolites are much lower than those of benzophenone-3. However, the metabolites are found to have skin-sensitization effects. The present study provides a deep insight into the biotransformation process of benzophenone-3 catalyzed by P450 1A1 and alerts us to pay attention to the adverse effects of benzophenone-3 and its metabolites in human livers.

Mechanistic and predictive studies on the oxidation of furans by cytochrome P450: A DFT study

Furan-containing compounds distribute widely in food, herbal medicines, industrial synthetic products, and environmental media. These compounds can undergo oxidative metabolism catalyzed by cytochrome P450 enzymes (CYP450) within organisms, which may produce reactive products, possibly reacting with biomolecules to induce toxic effects. In this work, we performed DFT calculations to investigate the CYP450-mediated metabolic mechanism of furan-ring oxidation using 2-methylfuran as a model substrate, meanwhile, we studied the regioselective competition of another hydroxylation reaction involving methyl group of 2-methylfuran. As a result, we found the toxicological-relevant cis-enedione product can be produced from O-addition directly via a concerted manner without formation of an epoxide intermediate as traditionally believed. Moreover, our calculations demonstrate the kinetic and thermodynamic feasibility of both furan-ring oxidation and methyl hydroxylation pathways, although the former pathway is a bit more favorable. We then constructed a linear model to predict the rate-limiting activation energies (ΔE*) of O-addition with 11 diverse furan substates based on their adiabatic ionization potentials (AIPs) and condensation Fukui functions (CFFs). The results show a good predictive ability (R2=0.94, Q2CV=0.87). Therefore, AIP and CFF with clear physichem meanings relevant to the mechanism, emerge as pivotal molecular descriptors to enable the fast prediction of furan-ring oxidation reactivities for quick insight into the toxicological risk of furans, using just ground-state calculations.

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,β-Dioxygenase

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways.

Enhanced Reactivity through Equatorial Sulfur Coordination in Nonheme Iron(IV)-Oxo Complexes: Insights from Experiment and Theory

Sulfur ligation in metalloenzymes often gives the active site unique properties, whether it is the axial cysteinate ligand in the cytochrome P450s or the equatorial sulfur/thiol ligation in nonheme iron enzymes. To understand sulfur ligation to iron complexes and how it affects the structural, spectroscopic, and intrinsic properties of the active species and the catalysis of substrates, we pursued a systematic study and compared sulfur with amine-ligated iron(IV)-oxo complexes. We synthesized and characterized a biomimetic N4S-ligated iron(IV)-oxo complex and compared the obtained results with an analogous N5-ligated iron(IV)-oxo complex. Our work shows that the amine for sulfur replacement in the equatorial ligand framework leads to a rate enhancement for oxygen atom and hydrogen atom transfer reactions. Moreover, the sulfur-ligated iron(IV)-oxo complex reacts through a different reaction mechanism as compared to the N5-ligated iron(IV)-oxo complex, where the former reacts through hydride transfer with the latter reacting via radical pathways. We show that the reactivity differences are caused by a dramatic change in redox potential between the two complexes. Our studies highlight the importance of implementing a sulfur ligand into the equatorial ligand framework of nonheme iron(IV)-oxo complexes and how it affects the physicochemical properties of the oxidant and its reactivity.

Insight into chemically reactive metabolites of aliphatic amine pollutants: A de novo prediction strategy and case study of sertraline

The uncommon metabolic pathways of organic pollutants are easily overlooked, potentially leading to idiosyncratic toxicity. Prediction of their biotransformation associated with the toxic effects is the very purpose that this work focuses, to develop a de novo method to mechanistically predict the reactive toxicity pathways of uncommon metabolites from start aliphatic amine molecules, which employed sertraline triggered by CYP450 enzymes as a model system, as there are growing concerns about the effects on human health posed by antidepressants in the aquatic environment. This de novo prediction strategy combines computational and experimental methods, involving DFT calculations upon sequential growth, in vitro and in vivo assays, dissecting chemically reactive mechanism relevant to toxicity, and rationalizing the fundamental factors. Significantly, desaturation and debenzylation-aromatization as the emerging metabolic pathways of sertraline have been elucidated, with the detection of DNA adducts of oxaziridine metabolite in mice, highlighting the potential reactive toxicity. Molecular orbital analysis supports the reactivity preference for toxicological-relevant C-N desaturation over N-hydroxylation of sertraline, possibly extended to several other aliphatic amines based on the Bell-Evans-Polanyi principle. It was further validated toward some other wide-concerned aliphatic amine pollutants involving atrazine, ε-caprolactam, 6PPD via in silico and in vitro assays, thereby constituting a complete path for de novo prediction from case study to general applications.

Catalytic divergencies in the mechanism of L-arginine hydroxylating nonheme iron enzymes

Many enzymes in nature utilize a free arginine (L-Arg) amino acid to initiate the biosynthesis of natural products. Examples include nitric oxide synthases, which generate NO from L-Arg for blood pressure control, and various arginine hydroxylases involved in antibiotic biosynthesis. Among the groups of arginine hydroxylases, several enzymes utilize a nonheme iron(II) active site and let L-Arg react with dioxygen and α-ketoglutarate to perform either C3-hydroxylation, C4-hydroxylation, C5-hydroxylation, or C4-C5-desaturation. How these seemingly similar enzymes can react with high specificity and selectivity to form different products remains unknown. Over the past few years, our groups have investigated the mechanisms of L-Arg-activating nonheme iron dioxygenases, including the viomycin biosynthesis enzyme VioC, the naphthyridinomycin biosynthesis enzyme NapI, and the streptothricin biosynthesis enzyme OrfP, using computational approaches and applied molecular dynamics, quantum mechanics on cluster models, and quantum mechanics/molecular mechanics (QM/MM) approaches. These studies not only highlight the differences in substrate and oxidant binding and positioning but also emphasize on electronic and electrostatic differences in the substrate-binding pockets of the enzymes. In particular, due to charge differences in the active site structures, there are changes in the local electric field and electric dipole moment orientations that either strengthen or weaken specific substrate C-H bonds. The local field effects, therefore, influence and guide reaction selectivity and specificity and give the enzymes their unique reactivity patterns. Computational work using either QM/MM or density functional theory (DFT) on cluster models can provide valuable insights into catalytic reaction mechanisms and produce accurate and reliable data that can be used to engineer proteins and synthetic catalysts to perform novel reaction pathways.

Melatonin Activation by Human Cytochrome P450 Enzymes: A Comparison between Different Isozymes

Cytochrome P450 enzymes in the human body play a pivotal role in both the biosynthesis and the degradation of the hormone melatonin. Melatonin plays a key role in circadian rhythms in the body, but its concentration is also linked to mood fluctuations as well as emotional well-being. In the present study, we present a computational analysis of the binding and activation of melatonin by various P450 isozymes that are known to yield different products and product distributions. In particular, the P450 isozymes 1A1, 1A2, and 1B1 generally react with melatonin to provide dominant aromatic hydroxylation at the C6-position, whereas the P450 2C19 isozyme mostly provides O-demethylation products. To gain insight into the origin of these product distributions of the P450 isozymes, we performed a comprehensive computational study of P450 2C19 isozymes and compared our work with previous studies on alternative isozymes. The work covers molecular mechanics, molecular dynamics and quantum mechanics approaches. Our work highlights major differences in the size and shape of the substrate binding pocket amongst the different P450 isozymes. Consequently, substrate binding and positioning in the active site varies substantially within the P450 isozymes. Thus, in P450 2C19, the substrate is oriented with its methoxy group pointing towards the heme, and therefore reacts favorably through hydrogen atom abstraction, leading to the production of O-demethylation products. On the other hand, the substrate-binding pockets in P450 1A1, 1A2, and 1B1 are tighter, direct the methoxy group away from the heme, and consequently activate an alternative site and lead to aromatic hydroxylation instead.

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.

Oriented External Electric Fields Regurating the Reaction Mechanism of CH4 Oxidation Catalyzed by Fe(IV)-Oxo-Corrolazine: Insight from Density Functional Calculations

Methane is the simplest alkane and can be used as an alternative energy source for oil and coal, but the greenhouse effect caused by its leakage into the air is not negligible, and its conversion into liquid methanol not only facilitates transportation, but also contributes to carbon neutrality. In order to find an efficient method for converting methane to methanol, CH4 oxidation catalyzed by Fe(IV)-Oxo-corrolazine (Fe(IV)-Oxo-Cz) and its reaction mechanism regulation by oriented external electric fields (OEEFs) are systematically studied by density functional calculations. The calculations show that Fe(IV)-Oxo-Cz can abstract one H atom from CH4 to form the intermediate with OH group connecting on the corrolazine ring, with the energy barrier of 25.44 kcal mol-1. And then the product methanol is formed through the following rebound reaction. Moreover, the energy barrier can be reduced to 20.72 kcal mol-1 through a two-state reaction pathway. Furthermore, the effect of OEEFs on the reaction is investigated. We found that OEEFs can effectively regulate the reaction by adjusting the stability of the reactant and the transition state through the interaction of electric field-molecular dipole moment. When the electric field is negative, the energy barrier of the reaction decreases with the increase of electric intensity. Moreover, the OEEF aligned along the intrinsic Fe‒O reaction axis can effectively regulate the ability of forming the OH on the corrolazine ring by adjusting the charges of O and H atoms. When the electric field intensity is -0.010 a.u., the OH can be directly rebounded to the CH3· before it is connecting on the corrolazine ring, thus forming the product directly from the transition state without passing through the intermediate with only an energy barrier of 17.34 kcal mol-1, which greatly improves the selectivity of the reaction.

Heme compound II models in chemoselectivity and disproportionation reactions

Heme compound II models bearing electron-deficient and -rich porphyrins, [FeIV(O)(TPFPP)(Cl)]- (1a) and [FeIV(O)(TMP)(Cl)]- (2a), respectively, are synthesized, spectroscopically characterized, and investigated in chemoselectivity and disproportionation reactions using cyclohexene as a mechanistic probe. Interestingly, cyclohexene oxidation by 1a occurs at the allylic C-H bonds with a high kinetic isotope effect (KIE) of 41, yielding 2-cyclohexen-1-ol product; this chemoselectivity is the same as that of nonheme iron(iv)-oxo intermediates. In contrast, as observed in heme compound I models, 2a yields cyclohexene oxide product with a KIE of 1, demonstrating a preference for C[double bond, length as m-dash]C epoxidation. The latter result is interpreted as 2a disproportionating to form [FeIV(O)(TMP+˙)]+ (2b) and FeIII(OH)(TMP), and 2b becoming the active oxidant to conduct the cyclohexene epoxidation. In contrast to 2a, 1a does not disproportionate under the present reaction conditions. DFT calculations confirm that compound II models prefer C-H bond hydroxylation and that disproportionation of compound II models is controlled thermodynamically by the porphyrin ligands. Other aspects, such as acid and base effects on the disproportionation of compound II models, have been discussed as well.

Enthalpy-Entropy Compensation Effect in Oxidation Reactions by Manganese(IV)-Oxo Porphyrins and Nonheme Iron(IV)-Oxo Models

"Enthalpy-Entropy Compensation Effect" (EECE) is ubiquitous in chemical reactions; however, such an EECE has been rarely explored in biomimetic oxidation reactions. In this study, six manganese(IV)-oxo complexes bearing electron-rich and -deficient porphyrins are synthesized and investigated in various oxidation reactions, such as hydrogen atom transfer (HAT), oxygen atom transfer (OAT), and electron-transfer (ET) reactions. First, all of the six Mn(IV)-oxo porphyrins are highly reactive in the HAT, OAT, and ET reactions. Interestingly, we have observed a reversed reactivity in the HAT and OAT reactions by the electron-rich and -deficient Mn(IV)-oxo porphyrins, depending on reaction temperatures, but not in the ET reactions; the electron-rich Mn(IV)-oxo porphyrins are more reactive than the electron-deficient Mn(IV)-oxo porphyrins at high temperature (e.g., 0 °C), whereas at low temperature (e.g., -60 °C), the electron-deficient Mn(IV)-oxo porphyrins are more reactive than the electron-rich Mn(IV)-oxo porphyrins. Such a reversed reactivity between the electron-rich and -deficient Mn(IV)-oxo porphyrins depending on reaction temperatures is rationalized with EECE; that is, the lower is the activation enthalpy, the more negative is the activation entropy, and vice versa. Interestingly, a unified linear correlation between the activation enthalpies and the activation entropies is observed in the HAT and OAT reactions of the Mn(IV)-oxo porphyrins. Moreover, from the previously reported HAT reactions of nonheme Fe(IV)-oxo complexes, a linear correlation between the activation enthalpies and the activation entropies is also observed. To the best of our knowledge, we report the first detailed mechanistic study of EECE in the oxidation reactions by synthetic high-valent metal-oxo complexes.

Nonheme Iron Imido Complexes Bearing a Non-Innocent Ligand: A Synthetic Chameleon Species in Oxidation Reactions

High-valent iron-imido complexes can perform C-H activation and sulfimidation reactions, but are far less studied than the more ubiquitous iron-oxo species. As case studies, we have looked at a recently published iron(V)-imido ligand π-cation radical complex, which is formally an iron(VI)-imido complex [Fe^V (NTs)(TAML^+. )] (1; NTs=tosylimido), and an iron(V)-imido complex [Fe^V (NTs)(TAML)]^- (2). Using a theoretical approach, we found that they have multiple energetically close-lying electromers, sometimes even without changing spin states, reminiscent of the so-called Compound I in Cytochrome P450. When studying their reactivity theoretically, it is indeed found that their electronic structures may change to perform efficient oxidations, emulating the multi-spin state reactivity in Fe^IV O systems. This is actually in contrast to the known [Fe^V (O)(TAML)]^- species (3), where the reactions occur only on the ground spin state. We also looked into the whole reaction pathway for the C-H bond activation of 1,4-cyclohexadiene by these intermediates to reproduce the experimentally observed products, including steps that usually attract no interest (neither theoretically nor experimentally) due to their non-rate-limiting status and fast reactivity. A new "clustering non-rebound mechanism" is presented for this C-H activation reaction.

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.

Reactive Metabolites from Thiazole-Containing Drugs: Quantum Chemical Insights into Biotransformation and Toxicity

Drugs containing thiazole and aminothiazole groups are known to generate reactive metabolites (RMs) catalyzed by cytochrome P450s (CYPs). These RMs can covalently modify essential cellular macromolecules and lead to toxicity and induce idiosyncratic adverse drug reactions. Molecular docking and quantum chemical hybrid DFT study were carried out to explore the molecular mechanisms involved in the biotransformation of thiazole (TZ) and aminothiazole (ATZ) groups leading to RM epoxide, S-oxide, N-oxide, and oxaziridine. The energy barrier required for the epoxidation is 13.63 kcal/mol, that is lower than that of S-oxidation, N-oxidation, and oxaziridine formation (14.56, 17.90, and 20.20, kcal/mol respectively). The presence of the amino group in ATZ further facilitates all the metabolic pathways, for example, the barrier for the epoxidation reaction is reduced by ∼2.5 kcal/mol. Some of the RMs/their isomers are highly electrophilic and tend to form covalent bonds with nucleophilic amino acids, finally leading to the formation of metabolic intermediate complexes (MICs). The energy profiles of these competitive pathways have also been explored.

Mechanism-Based Inactivation of Cytochrome P450 Enzymes: Computational Insights

Mechanism-based inactivation (MBI) refers to the metabolic bioactivation of a xenobiotic by cytochrome P450s to a highly reactive intermediate which subsequently binds to the enzyme and leads to the quasi-irreversible or irreversible inhibition. Xenobiotics, mainly drugs with specific functional units, are the major sources of MBI. Two possible consequences of MBI by medicinal compounds are drug-drug interaction and severe toxicity that are observed and highlighted by clinical experiments. Today almost all of these latent functional groups (e.g., thiophene, furan, alkylamines, etc.) are known, and their features and mechanisms of action, owing to the vast experimental and theoretical studies, are determined. In the past decade, molecular modeling techniques, mostly density functional theory, have revealed the most feasible mechanism that a drug undergoes by P450 enzymes to generate a highly reactive intermediate. In this review, we provide a comprehensive and detailed picture of computational advances toward the elucidation of the activation mechanisms of various known groups with MBI activity. To this aim, we briefly describe the computational concepts to carry out and analyze the mechanistic investigations, and then, we summarize the studies on compounds with known inhibition activity including thiophene, furan, alkylamines, terminal acetylene, etc. This study can be reference literature for both theoretical and experimental (bio)chemists in several different fields including rational drug design, the process of toxicity prevention, and the discovery of novel inhibitors and catalysts.

Metalloporphyrins in Medicine: From History to Recent Trends

The history of metalloporphyrins dates back more than 200 years ago. Metalloporphyrins are excellent catalysts, capable of forming supramolecular systems, participate in oxygen photosynthesis, transport, and used as contrast agents or superoxide dismutase mimetics. Today, metalloporphyrins represent complexes of conjugated π-electron system and metals from the entire periodic system. However, the effect of these compounds on living systems has not been fully understood, and researchers are exploring the properties of metalloporphyrins thereby extending their further application. This review provides an overview of the variety of metalloporphyrins that are currently used in different medicine fields and how metalloporphyrins became the subject of scientists' interest. Currently, metalloporphyrins utilization has expanded significantly, which gave us an opprotunuty to summarize recent progress in metalloporphyrins derivatives and prospects of their application in the treatment and diagnosis of different diseases.

Oxo-Free Hydrocarbon Oxidation by an Iron(III)-Isoporphyrin Complex

Metal-halides that perform proton coupled electron-transfer (PCET) oxidation are an important new class of high-valent oxidant. In investigating metal-dihalides, we reacted [Fe^III(Cl)(T(OMe)PP)] (1, T(OMe)PP = meso-tetra(4-methoxyphenyl)porphyrinyl) with (dichloroiodo)benzene. An Fe^III-meso-chloro-isoporphyrin complex [Fe^III(Cl)₂(T(OMe)PP-Cl)] (2) was obtained. 2 was characterized by electronic absorption, ¹H NMR, EPR, and X-ray absorption spectroscopies and mass spectrometry with support from computational analyses. 2 was reacted with a series of hydrocarbon substrates. The measured kinetic data exhibited a nonlinear behavior, whereby the oxidation followed a hydrogen-atom-transfer (HAT) PCET mechanism. The meso-chlorine atom was identified as the HAT agent. In one case, a halogenated product was identified by mass spectrometry. Our findings demonstrate that oxo-free hydrocarbon oxidation with heme systems is possible and show the potential for iron-dihalides in oxidative hydrocarbon halogenation.

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.

Theoretical Study on the Structural-Function Relationship of Manganese(III)-Iodosylarene Adducts

Metal-iodosylarene complexes have been recently viewed as a second oxidant alongside of the well-known high-valent metal-oxo species. Extensive efforts have been exerted to unveil the structure-function relationship of various metal-iodosylarene complexes. In the present manuscript, density functional theoretical calculations were employed to investigate such relationship of a specific manganese-iodosylbenzene complex [Mn^III(TBDAP)(PhIO)(OH)]²⁺ (1). Our results fit the experimental observations and revealed new mechanistic findings. 1 acts as a stepwise 1e+1e oxidant in sulfoxidation reactions. Surprisingly, C-H bond activation of 9,10-dihydroanthracene (DHA) by 1 proceeds via a novel ionic hydride transfer/proton transfer (HT/PT) mechanism. As a comparison to 1, the electrophilicity of an iodosylbenzene monomer PhIO was investigated. PhIO performs concerted 2e-oxidations both in sulfoxidation and C-H activation. Hydroxylation of DHA by PhIO was found to proceed via a novel ionic and concerted proton-transfer/hydroxyl-rebound mechanism involving 2e-oxidation to form a transient carbonium species.

Reactivity patterns of vanadium(iv/v)-oxo complexes with olefins in the presence of peroxides: a computational study

Vanadium porphyrin complexes are naturally occurring substances found in crude oil and have been shown to have medicinal properties as well. Little is known on their activities with substrates; therefore, we decided to perform a detailed density functional theory study on the properties and reactivities of vanadium(iv)- and vanadium(v)-oxo complexes with a TPPCl8 or 2,3,7,8,12,13,17,18-octachloro-meso-tetraphenylporphyrinato ligand system. In particular, we investigated the reactivity of [VV(O)(TPPCl8)]+ and [VIV(O)(TPPCl8)] with cyclohexene in the presence of H2O2 or HCO4-. The work shows that vanadium(iv)-oxo and vanadium(v)-oxo are sluggish oxidants by themselves and react with olefins slowly. However, in the presence of hydrogen peroxide, these metal-oxo species can be transformed into a side-on vanadium-peroxo complex, which reacts with substrates more efficiently. Particularly with anionic axial ligands, the side-on vanadium-peroxo and vanadium-oxo complexes produced epoxides from cyclohexene via small barrier heights. In addition to olefin epoxidation, we investigated aliphatic hydroxylation mechanisms by the same oxidants and some oxidants show efficient and viable cyclohexene hydroxylation mechanisms. The work implies that vanadium-oxo and vanadium-peroxo complexes can react with double bonds through epoxidation, and under certain conditions also undergo hydroxylation, but the overall reactivity is highly dependent on the equatorial ligand, the local environment and the presence or absence of anionic axial ligands.

Heterogenization of manganese porphyrin via hydrogen bond in zeolite imidazolate framework-8 matrix, a host–guest interaction, as catalytic system for olefin epoxidation

A heterogenized meso-tetrakis(2,3-dihydroxyphenyl)porphyrinatomanganese(III) acetate at zeolite imidazolate framework-8 (T(2,3-OHP)PorMn@ZIF-8) is investigated for the catalytic olefin epoxidation reactions at room temperature. Heterogenization is accomplished through a non-classical hydrogen bond proposed between T(2,3-OHP)PorMn bearing O–H groups and C–H of the 2-methylimidazolate linkers in the ZIF-8 structure. The aforementioned compound is characterized by X-ray powder diffraction (XRD), atomic absorption spectroscopy (AAS), nitrogen adsorption−desorption, FT-IR spectroscopy, field emission scanning electron microscopy (FE-SEM), transmission electron microscopy (TEM) and Raman spectroscopy. The catalytic system with rather high potential of reusability is proposed as a fairly efficient epoxidation catalyst compared to reports in homogeneous media.

Quantum Mechanics/Molecular Mechanics Studies on the Relative Reactivities of Compound I and II in Cytochrome P450 Enzymes

The cytochromes P450 are drug metabolizing enzymes in the body that typically react with substrates through a monoxygenation reaction. During the catalytic cycle two reduction and protonation steps generate a high-valent iron (IV)-oxo heme cation radical species called Compound I. However, with sufficient reduction equivalents present, the catalytic cycle should be able to continue to the reduced species of Compound I, called Compound II, rather than a reaction of Compound I with substrate. In particular, since electron transfer is usually on faster timescales than atom transfer, we considered this process feasible and decided to investigate the reaction computationally. In this work we present a computational study using density functional theory methods on active site model complexes alongside quantum mechanics/molecular mechanics calculations on full enzyme structures of cytochrome P450 enzymes. Specifically, we focus on the relative reactivity of Compound I and II with a model substrate for O⁻H bond activation. We show that generally the barrier heights for hydrogen atom abstraction are higher in energy for Compound II than Compound I for O⁻H bond activation. Nevertheless, for the activation of such bonds, Compound II should still be an active oxidant under enzymatic conditions. As such, our computational modelling predicts that under high-reduction environments the cytochromes P450 can react with substrates via Compound II but the rates will be much slower.

Mechanistic Insight on the Activity and Substrate Selectivity of Nonheme Iron Dioxygenases

Nonheme iron dioxygenases catalyze vital reactions for human health particularly related to aging processes. They are involved in the biosynthesis of amino acids, but also the biodegradation of toxic compounds. Typically they react with their substrate(s) through oxygen atom transfer, although often with the assistance of a co-substrate like α-ketoglutarate that is converted to succinate and CO₂ . Many reaction processes catalyzed by the nonheme iron dioxygenases are stereoselective or regiospecific and hence understanding the mechanism and protein involvement in the selectivity is important for the design of biotechnological applications of these enzymes. To this end, I will review recent work of our group on nonheme iron dioxygenases and include background information on their general structure and catalytic cycle. Examples of stereoselective and regiospecific reaction mechanisms we elucidated are for the AlkB repair enzyme, prolyl-4-hydroxylase and the ergothioneine biosynthesis enzyme. Finally, I cover an example where we bioengineered S-p-hydroxymandelate synthase into the R-p-hydroxymandelate synthase.

Reactivity Patterns of (Protonated) Compound II and Compound I of Cytochrome P450: Which is the Better Oxidant?

The cytochromes P450 are versatile enzymes in human physiology that perform substrate hydroxylation reactions extremely efficiently. In this work, we present results of a computational study on the reactivity patterns of Compound I, Compound II, and protonated Compound II with model substrates, and we address the question of which of these compounds is the most effective oxidant? All calculations, regardless of the substrate, implicated that Compound I is the superior oxidant of the three. However, Compound II and protonated Compound II were found to react with free energies of activation that are only a few kcal mol^-1 higher in energy than those obtained with Compound I. Therefore, Compound II and protonated Compound II should be able to react with aliphatic groups with moderate C-H bond strengths. We have analysed all results in detail and have given electronic, thermochemical, valence bond, and molecular orbital rationalizations on the reactivity differences and explained experimental product distributions. Overall, the findings implied that alternative oxidants could operate alongside Compound I in complex reaction mechanisms of enzymatic and synthetic iron porphyrinoid complexes.

Tunneling Effect That Changes the Reaction Pathway from Epoxidation to Hydroxylation in the Oxidation of Cyclohexene by a Compound I Model of Cytochrome P450

The rate constants of the C═C epoxidation and the C-H hydroxylation (i.e., allylic C-H bond activation) in the oxidation of cyclohexene by a high-valent iron(IV)-oxo porphyrin π-cation radical complex, [(TMP^•+)Fe^IV(O)(Cl)] (1, TMP = meso-tetramesitylporphyrin dianion), were determined at various temperatures by analyzing the overall rate constants and the products obtained in the cyclohexene oxidation by 1, leading us to conclude that reaction pathway changes from the C═C epoxidation to C-H hydroxylation by decreasing reaction temperature. When cyclohexene was replaced by deuterated cyclohexene (cyclohexene-d₁₀), the epoxidation pathway dominated irrespective of the reaction temperature. The temperature dependence of the rate constant of the C-H hydroxylation pathway in the reactions of cyclohexene and cyclohexene-d₁₀ by 1 suggests that there is a significant tunneling effect on the hydrogen atom abstraction of allylic C-H bonds of cyclohexene by 1, leading us to propose that the tunneling effect is a determining factor for the switchover of the reaction pathway from the C═C epoxidation pathway to the C-H hydroxylation pathway by decreasing reaction temperature. By performing density functional theory (DFT) calculations, the reaction energy barriers of the C═C epoxidation and C-H bond activation reactions by 1 were found to be similar, supporting the notion that small environmental changes, such as the reaction temperature, can flip the preference for one reaction to another.

Oxygen Atom Transfer Using an Iron(IV)-Oxo Embedded in a Tetracyclic N-Heterocyclic Carbene System: How Does the Reactivity Compare to Cytochrome P450 Compound I?

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

Substrate Sulfoxidation by an Iron(IV)-Oxo Complex: Benchmarking Computationally Calculated Barrier Heights to Experiment

High-valent metal-oxo oxidants are common reactive species in synthetic catalysts as well as heme and nonheme iron enzymes. In general, they efficiently react with substrates through oxygen atom transfer, and for a number of cases, experimental rate constants have been determined. However, because these rate constants are generally measured in a polar solution, it has been found difficult to find computational methodologies to reproduce experimental trends and reactivities. In this work, we present a detailed computational study into para-substituted thioanisole sulfoxidation by a nonheme iron(IV)-oxo complex. A range of density functional theory methods and basis sets has been tested for their suitability to describe the reaction mechanism and compared with experimentally obtained free energies of activation. It is found that the enthalpy of activation is reproduced well, but all methods overestimate the entropy of activation by about 50%, for which we recommend a correction factor. The effect of solvent and dispersion on the barrier heights is explored both at the single-point level and also through inclusion in geometry optimizations, and particularly, solvent is seen as highly beneficial to reproduce experimental free energies of activation. Interestingly, in general, experimental trends and Hammett plots are reproduced well with almost all methods and procedures, and only a systematic error seems to apply for these chemical systems. Very good agreement between experiment and theory is found for a number of different methods, including B3LYP and PBE0, and procedures that are highlighted in the paper.

Toxicity Originating from Thiophene Containing Drugs: Exploring the Mechanism using Quantum Chemical Methods

Drug metabolism of thiophene containing substrates by cytochrome P450s (CYP450) leads to toxic side effects, for example, nephrotoxicity (suprofen, ticlopidine), hepatotoxicity (tienilic acid), thrombotic thrombocytopenic purpura (clopidogrel), and aplastic anemia (ticlopidine). The origin of toxicity in these cases has been attributed to two different CYP450 mediated metabolic reactions: S-oxidation and epoxidation. In this work, the molecular level details of the bioinorganic chemistry associated with the generation of these competitive reactions are reported. Density functional theory was utilized (i) to explore the molecular mechanism for S-oxidation and epoxidation using the radical cationic center Cpd I [(iron(IV)-oxo-heme porphine system with SH(-) as the axial ligand, to mimic CYP450s] as the model oxidant, (ii) to establish the 3D structures of the reactants, transition states, and products on both the metabolic pathways, and (iii) to examine the potential energy (PE) profile for both the pathways to determine the energetically preferred toxic metabolite formation. The energy barrier required for S-oxidation was observed to be 14.75 kcal/mol as compared to that of the epoxidation reaction (13.23 kcal/mol) on the doublet PE surface of Cpd I. The formation of the epoxide metabolite was found to be highly exothermic (-23.24 kcal/mol), as compared to S-oxidation (-8.08 kcal/mol). Hence, on a relative scale the epoxidation process was observed to be thermodynamically and kinetically more favorable. The energy profiles associated with the reactions of the S-oxide and epoxide toxic metabolites were also explored. This study helps in understanding the CYP450-catalyzed toxic reactions of drugs containing the thiophene ring at the atomic level.

Site-selective formation of an iron(iv)-oxo species at the more electron-rich iron atom of heteroleptic μ-nitrido diiron phthalocyanines

A combination of MS and computation on μ-nitrido bridged diiron complexes reveals H₂O₂ binding to the complex and generates an oxidant capable of oxidizing methane.

Iron(iv)–oxo species have been identified as the active intermediates in key enzymatic processes, and their catalytic properties are strongly affected by the equatorial and axial ligands bound to the metal, but details of these effects are still unresolved. In our aim to create better and more efficient oxidants of H-atom abstraction reactions, we have investigated a unique heteroleptic diiron phthalocyanine complex. We propose a novel intramolecular approach to determine the structural features that govern the catalytic activity of iron(iv)–oxo sites. Heteroleptic μ-nitrido diiron phthalocyanine complexes having an unsubstituted phthalocyanine (Pc1) and a phthalocyanine ligand substituted with electron-withdrawing alkylsulfonyl groups (PcSO₂R) were prepared and characterized. A reaction with terminal oxidants gives two isomeric iron(iv)–oxo and iron(iii)–hydroperoxo species with abundances dependent on the equatorial ligand. Cryospray ionization mass spectrometry (CSI-MS) characterized both hydroperoxo and diiron oxo species in the presence of H₂O₂. When m-CPBA was used as the oxidant, the formation of diiron oxo species (PcSO₂R)FeNFe(Pc1) Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O was also evidenced. Sufficient amounts of these transient species were trapped in the quadrupole region of the mass-spectrometer and underwent a CID-MS/MS fragmentation. Analyses of fragmentation patterns indicated a preferential formation of hydroperoxo and oxo moieties at more electron-rich iron sites of both heteroleptic μ-nitrido complexes. DFT calculations show that both isomers are close in energy. However, the analysis of the iron(iii)–hydroperoxo bond strength reveals major differences for the (Pc1)FeN(PcSO₂R)Fe^IIIOOH system as compared to (PcSO₂R)FeN(Pc1)Fe^IIIOOH system, and, hence binding of a terminal oxidant will be preferentially on more electron-rich sides. Subsequent kinetics studies showed that these oxidants are able to even oxidize methane to formic acid efficiently.

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

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

Drug metabolism by cytochrome p450 enzymes: what distinguishes the pathways leading to substrate hydroxylation over desaturation?

Cytochrome P450 enzymes are highly versatile biological catalysts in our body that react with a broad range of substrates. Key functions in the liver include the metabolism of drugs and xenobiotics. One particular metabolic pathway that is poorly understood relates to the P450 activation of aliphatic groups leading to either hydroxylation or desaturation pathways. A DFT and QM/MM study has been carried out on the factors that determine the regioselectivity of aliphatic hydroxylation over desaturation of compounds by P450 isozymes. The calculations establish multistate reactivity patterns, whereby the product distributions differ on each of the spin-state surfaces; hence spin-selective product formation was found. The electronic and thermochemical factors that determine the bifurcation pathways were analysed and a model that predicts the regioselectivity of aliphatic hydroxylation over desaturation pathways was established from valence bond and molecular orbital theories. Thus, the difference in energy of the OH versus the OC bond formed and the π-conjugation energy determines the degree of desaturation products. In addition, environmental effects of the substrate binding pocket that affect the regioselectivities were identified. These studies imply that bioengineering P450 isozymes for desaturation reactions will have to include modifications in the substrate binding pocket to restrict the hydroxylation rebound reaction.

A comprehensive test set of epoxidation rate constants for iron(iv)-oxo porphyrin cation radical complexes

Cytochrome P450 enzymes are heme based monoxygenases that catalyse a range of oxygen atom transfer reactions with various substrates, including aliphatic and aromatic hydroxylation as well as epoxidation reactions. The active species is short-lived and difficult to trap and characterize experimentally, moreover, it reacts in a regioselective manner with substrates leading to aliphatic hydroxylation and epoxidation products, but the origin of this regioselectivity is poorly understood. We have synthesized a model complex and studied it with low-pressure Fourier transform-ion cyclotron resonance (FT-ICR) mass spectrometry (MS). A novel approach was devised using the reaction of [FeIII(TPFPP)]+ (TPFPP = meso-tetrakis(pentafluorophenyl)porphinato dianion) with iodosylbenzene as a terminal oxidant which leads to the production of ions corresponding to [FeIV(O)(TPFPP+˙)]+. This species was isolated in the gas-phase and studied in its reactivity with a variety of olefins. Product patterns and rate constants under Ideal Gas conditions were determined by FT-ICR MS. All substrates react with [FeIV(O)(TPFPP+˙)]+ by a more or less efficient oxygen atom transfer process. In addition, substrates with low ionization energies react by a charge-transfer channel, which enabled us to determine the electron affinity of [FeIV(O)(TPFPP+˙)]+ for the first time. Interestingly, no hydrogen atom abstraction pathways are observed for the reaction of [FeIV(O)(TPFPP+˙)]+ with prototypical olefins such as propene, cyclohexene and cyclohexadiene and also no kinetic isotope effect in the reaction rate is found, which suggests that the competition between epoxidation and hydroxylation - in the gas-phase - is in favour of substrate epoxidation. This notion further implies that P450 enzymes will need to adapt their substrate binding pocket, in order to enable favourable aliphatic hydroxylation over double bond epoxidation pathways. The MS studies yield a large test-set of experimental reaction rates of iron(iv)-oxo porphyrin cation radical complexes, so far unprecedented in the gas-phase, providing a benchmark for calibration studies using computational techniques. Preliminary computational results presented here confirm the observed trends excellently and rationalize the reactivities within the framework of thermochemical considerations and valence bond schemes.