Can Ferric-Superoxide Act as a Potential Oxidant in P450cam? QM/MM Investigation of Hydroxylation, Epoxidation, and Sulfoxidation

Substrate Conformational Switch Enables the Stereoselective Dimerization in P450 NascB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations

P450 NascB catalyzes the coupling of cyclo-(l-tryptophan-l-proline) (1) to generate (-)-naseseazine C (2) through intramolecular C-N bond formation and intermolecular C-C coupling. A thorough understanding of its catalytic mechanism is crucial for the engineering or design of P450-catalyzed C-N dimerization reactions. By employing MD simulations, QM/MM calculations, and enhanced sampling, we assessed various mechanisms from recent works. Our study demonstrates that the most favorable pathway entails the transfer of a hydrogen atom from N7-H to Cpd I. Subsequently, there is a conformational change in the substrate radical, shifting it from the Re-face to the Si-face of N7 in Substrate 1. The Si-face conformation of Substrate 1 is stabilized by the protein environment and the π-π stacking interaction between the indole ring and heme porphyrin. The subsequent intermolecular C3-C6' bond formation between Substrate 1 radical and Substrate 2 occurs via a radical attack mechanism. The conformational switch of the Substrate 1 radical not only lowers the barrier of the intermolecular C3-C6' bond formation but also yields the correct stereoselectivity observed in experiments. In addition, we evaluated the reactivity of the ferric-superoxide species, showing it is not reactive enough to initiate the hydrogen atom abstraction from the indole NH group of the substrate. Our simulation provides a comprehensive mechanistic insight into how the P450 enzyme precisely controls both the intramolecular C-N cyclization and intermolecular C-C coupling. The current findings align with the available experimental data, emphasizing the pivotal role of substrate dynamics in governing P450 catalysis.

Deciphering the Molecular Mechanisms of Reactive Metabolite Formation in the Mechanism-Based Inactivation of Cytochrome p450 1B1 by 8-Methoxypsoralen and Assessing the Driving Effect of phe268

This study provides a comprehensive computational exploration of the inhibitory activity and metabolic pathways of 8-methoxypsoralen (8-MP), a furocoumarin derivative used for treating various skin disorders, on cytochrome P450 (P450). Employing quantum chemical DFT calculations, molecular docking, and molecular dynamics (MD) simulations analyses, the biotransformation mechanisms and the active site binding profile of 8-MP in CYP1B1 were investigated. Three plausible inactivation mechanisms were minutely scrutinized. Further analysis explored the formation of reactive metabolites in subsequent P450 metabolic processes, including covalent adduct formation through nucleophilic addition to the epoxide, 8-MP epoxide hydrolysis, and non-CYP-catalyzed epoxide ring opening. Special attention was paid to the catalytic effect of residue Phe268 on the mechanism-based inactivation (MBI) of P450 by 8-MP. Energetic profiles and facilitating conditions revealed a slight preference for the C4'=C5' epoxidation pathway, while recognizing a potential kinetic competition with the 8-OMe demethylation pathway due to comparable energy demands. The formation of covalent adducts via nucleophilic addition, particularly by phenylalanine, and the generation of potentially harmful reactive metabolites through autocatalyzed ring cleavage are likely to contribute significantly to P450 metabolism of 8-MP. Our findings highlight the key role of Phe268 in retaining 8-MP within the active site of CYP1B1, thereby facilitating initial oxygen addition transition states. This research offers crucial molecular-level insights that may guide the early stages of drug discovery and risk assessment related to the use of 8-MP.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Heme-copper and Heme O2-derived synthetic (bioinorganic) chemistry toward an understanding of cytochrome c oxidase dioxygen chemistry

Cytochrome c oxidase (CcO), also widely known as mitochondrial electron-transport-chain complex IV, is a multi-subunit transmembrane protein responsible for catalyzing the last step of the electron transport chain, dioxygen reduction to water, which is essential to the establishment and maintenance of the membrane proton gradient that drives ATP synthesis. Although many intermediates in the CcO catalytic cycle have been spectroscopically and/or computationally authenticated, the specifics regarding the IP intermediate, hypothesized to be a heme-Cu (hydro)peroxo species whose O-O bond homolysis is supported by a hydrogen-bonding network of water molecules, are largely obscured by the fast kinetics of the A (FeIII-O2•-/CuI/Tyr) → PM (FeIV=O/CuII-OH/Tyr•) step. In this review, we have focused on the recent advancements in the design, development, and characterization of synthetic heme-peroxo‑copper model complexes, which can circumvent the abovementioned limitation, for the investigation of the formation of IP and its O-O cleavage chemistry. Novel findings regarding (a) proton and electron transfer (PT/ET) processes, together with their contributions to exogenous phenol induced O-O cleavage, (b) the stereo-electronic tunability of the secondary coordination sphere (especially hydrogen-bonding) on the geometric and spin state alteration of the heme-peroxo‑copper unit, and (c) a plausible mechanism for the Tyr-His cofactor biogenesis, are discussed in great detail. Additionally, since the ferric-superoxide and the ferryl-oxo (Compound II) species are critically involved in the CcO catalytic cycle, this review also highlights a few fundamental aspects of these heme-only (i.e., without copper) species, including the structural and reactivity influences of electron-donating trans-axial ligands and Lewis acid-promoted H-bonding.

A Ferric-Superoxide Intermediate Initiates P450-Catalyzed Cyclic Dipeptide Dimerization

The cytochrome P450 (CYP) AspB is involved in the biosynthesis of the diketopiperazine (DKP) aspergilazine A. Tryptophan-linked dimeric DKP alkaloids are a large family of natural products that are found in numerous species and exhibit broad and often potent bioactivity. The proposed mechanisms for C-N bond formation by AspB, and similar C-C bond formations by related CYPs, have invoked the use of a ferryl-intermediate as an oxidant to promote substrate dimerization. Here, the parallel application of steady-state and transient kinetic approaches reveals a very different mechanism that involves a ferric-superoxide species as a primary oxidant to initiate DKP-assembly. Single turnover kinetic isotope effects and a substrate analog suggest the probable nature and site for abstraction. The direct observation of CYP-superoxide reactivity rationalizes the atypical outcome of AspB and reveals a new reaction manifold in heme enzymes.

The Role of Heme Peroxo Oxidants in the Rational Mechanistic Modeling of Nitric Oxide Synthase: Characterization of Key Intermediates and Elucidation of the Mechanism

Mammalian nitric oxide synthase (NOS) mediates the two-step O2 -dependent oxidative degradation of arginine, and has been linked to a medley of disease situations in humans. Nonetheless, its exact mechanism of action still remains unclear. This work presents the first NOS model system where biologically proposed heme superoxo and peroxo intermediates are assessed as active oxidants against oxime substrates. Markedly, heme peroxo intermediates engaged in a bioinspired oxime oxidation reaction pathway, converting oximes to ketones and nitroxyl anions (NO- ). Detailed thermodynamic, kinetic, and mechanistic interrogations all evince a rate-limiting step primarily driven by the nucleophilicity of the heme peroxo moiety. Coherent with other findings, 18 O and 15 N isotope substitution experiments herein suffice compelling evidence toward a detailed mechanism, which draw close parallels to one of the enzymatic proposals. Intriguingly, recent enzymatic studies also lend credence to these findings, and several relevant reaction intermediates have been observed during NOS turnover.

In Silico simulation of Cytochrome P450-Mediated metabolism of aromatic amines: A case study of N-Hydroxylation

Aromatic amines, the widely used raw materials in industry, cause long-term exposure to human bodies. They can be metabolized by cytochrome P450 enzymes to form active electrophilic compounds, which will potentially react with nucleophilic DNA to exert carcinogenesis. The short lifetime and versatility of the oxidant (a high-valent iron (IV)-oxo species, compound I) of P450 enzymes prompts us to use theoretical methods to investigate the metabolism of aromatic amines. In this work, the density functional theory (DFT) has been employed to simulate the hydroxylation metabolism through H-abstraction and to calculate the activation energy of this reaction for 28 aromatic amines. The results indicate that the steric effects, inductive effects and conjugative effects greatly contribute to the metabolism activity of the chemicals. The further correlation reveals that the dissociation energy of -NH₂ (BDE_N-H) can successfully predict the time-consuming calculated activation energy (R² for aromatic and heteroaromatic amines are 0.93 and 0.86, respectively), so BDE_N-H can be taken as a key parameter to characterize the relative stability of aromatic amines in P450 enzymes and further to quickly assess their potential toxicity. The validation results prove such relationship has good statistical performance (q_cv² for aromatic and heteroaromatic amines are 0.95 and 0.90, respectively) and can be used to other aromatic amines in the application domain, greatly reducing computational cost and providing useful support for experimental research.

Triphenyllead Hydroperoxide: A 1D Coordination Peroxo Polymer, Single-Crystal-to-Single-Crystal Disproportionation to a Superoxo/Hydroxo Complex, and Application in Catalysis

The synthesis, transformation, and application in catalysis of triphenyllead hydroperoxide, the first dioxygen lead complex, are described. Triphenyllead hydroperoxide is characterized by ²⁰⁷Pb nuclear magnetic resonance (NMR), Fourier transform infrared spectroscopy (FTIR), Raman spectroscopy, and single-crystal X-ray diffraction, revealing the first one-dimensional (1D) coordination peroxo polymer. Photolytic isomorphous transformation of Ph₃PbOOH yields a mixed hydroxo/superoxo crystalline structure, the first nonalkali superoxo crystalline metal salt, which is stable up to 100 °C. Upon further photolysis, another isomorphous transformation of the superoxide to hydroxide is observed. These are the first single-crystal-to-single-crystal hydroperoxide-to-superoxide and then to hydroxide transformations reported to date. Photolysis of triphenyllead hydroperoxide yields two forms of superoxide-doped crystalline structures that are distinguished by widely different characteristic relaxation times. The use of Ph₃PbOOH as an easy-to-handle solid two-electron oxidant for the highly enantioselective epoxidation of olefins is described.

Selective Photocatalyst for styrene epoxidation with atmospheric O2 using covalent organic frameworks

Fe@POG-OH was synthetized and used to photo-catalyze styrene epoxidation with high selectivity and high conversion at room temperature. O2˙− plays crucial roles in the effective and selective oxidation of styrene to styrene oxide.

Computational Study of the C5-Hydroxylation Mechanism Catalyzed by the Diiron Monooxygenase PtmU3 as Part of the Platensimycin Biosynthesis

PtmU3 is a newly identified nonheme diiron monooxygenase, which installs a C-5 β-hydroxyl group into the C-19 CoA-ester intermediate involved in the biosynthesis of unique diterpene-derived scaffolds of platensimycin and platencin. PtmU3 possesses a noncanonical diiron active site architecture of a saturated six-coordinate iron center and lacks the μ-oxo bridge. Although the hydroxylation process is a simple reaction for nonheme mononuclear iron-dependent enzymes, how PtmU3 employs the diiron center to catalyze the H-abstraction and OH-rebound is still unknown. In particular, the electronic characteristic of diiron is also unclear. To understand the catalytic mechanism of PtmU3, we constructed two reactant models in which both the Fe1^II-Fe2^III-superoxo and Fe1^II-Fe2^IV═O are considered to trigger the H-abstraction and performed a series of quantum mechanics/molecular mechanics calculations. Our calculation results reveal that PtmU3 is a special monooxygenase, that is, both atoms of the dioxygen molecule can be incorporated into two molecules of the substrate by the successive reactions. In the first-round reaction, PtmU3 uses the Fe1^II-Fe2^III-superoxo to install a hydroxyl group into the substrate, generating the high-reactive Fe1^II-Fe2^IV═O complex. In the second-round reaction, the Fe1^II-Fe2^IV═O species is responsible for the hydroxylation of another molecule of the substrate. In the diiron center, Fe2 adopts the high spin state (S = 5/2) during the catalysis, whereas for Fe1, in addition to its structural role, it may also play an assistant role for Fe1 catalysis. In the two successive OH-installing steps, the H-abstraction is always the rate-liming step. E241 and D308 not only act as bridging ligands to connect two Fe ions but also take part in the electron reorganization. Owing to the high reactivity of Fe1^II-Fe2^IV═O compared to Fe1^II-Fe2^III-superoxo, besides the C5-hydroxylation, the C3- or C18-hydroxylation was also calculated to be feasible.

Proton-coupled electron transfer reactivities of electronically divergent heme superoxide intermediates: a kinetic, thermodynamic, and theoretical study

Heme superoxides are one of the most versatile metallo-intermediates in biology, and they mediate a vast variety of oxidation and oxygenation reactions involving O2(g). Overall proton-coupled electron transfer (PCET) processes they facilitate may proceed via several different mechanistic pathways, attributes of which are not yet fully understood. Herein we present a detailed investigation into concerted PCET events of a series of geometrically similar, but electronically disparate synthetic heme superoxide mimics, where unprecedented, PCET feasibility-determining electronic effects of the heme center have been identified. These electronic factors firmly modulate both thermodynamic and kinetic parameters that are central to PCET, as supported by our experimental and theoretical observations. Consistently, the most electron-deficient superoxide adduct shows the strongest driving force for PCET, whereas the most electron-rich system remains unreactive. The pivotal role of these findings in understanding significant heme systems in biology, as well as in alternative energy applications is also discussed.

Inactivation Mechanism of Neuronal Nitric Oxide Synthase by (S)-2-Amino-5-(2-(methylthio)acetimidamido)pentanoic Acid: Chemical Conversion of the Inactivator in the Active Site

Neuronal nitric oxide synthase (nNOS) is one of the three isoforms of nitric oxide synthase (NOS). The other two isoforms include inducible NOS (iNOS) and endothelial NOS (eNOS). These three isoforms of NOS are widely present in both human and other mammals and are responsible for the biosynthesis of NO. As an essential biological molecule, NO plays an essential role in neurotransmission, immune response, and vasodilation; however, the overproduction of NO can cause a series of diseases. Thus, the selective inhibition of three isoforms of NOS has been considered to be important in treating related diseases. The active sites of the three enzymes are highly conserved, causing the selective inhibition of the three enzymes to be a great challenge. (S)-2-Amino-5-(2-(methylthio)acetimidamido)pentanoic acid (1) has been experimentally proved to be a selective and time-dependent irreversible inhibitor of nNOS, and three pathways, including sulfide oxidation, oxidative dethiolation, and oxidative demethylation, have been suggested. In this work, we performed quantum mechanics/molecular mechanics calculations to verify the chemical conversion of inactivator 1. Although we agree with the previously suggested chemical transformation process, our calculations demonstrated that there are lower energy pathways to accomplish both oxidative dethiolation and oxidative demethylation. These three branching reactions are competitive, but only dethiolation and demethylation reactions can generate inhibitory intermediates. As a powerful time-dependent irreversible inhibitor of nNOS, the key sulfur atom and middle imine are all necessary. Our calculation results not only verified the chemical reaction of inhibitor 1 occurring in the enzymatic active site but also explained the inactivation mechanism of inhibitor 1. This is also the first verified example of the heme-enzyme-catalyzed S-demethylation mechanism.

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.

Heme-FeIII Superoxide, Peroxide and Hydroperoxide Thermodynamic Relationships: FeIII-O2•- Complex H-Atom Abstraction Reactivity

Establishing redox and thermodynamic relationships between metal-ion-bound O₂ and its reduced (and protonated) derivatives is critically important for a full understanding of (bio)chemical processes involving dioxygen processing. Here, a ferric heme peroxide complex, [(F₈)Fe^III-(O₂^2-)]^- (P) (F₈ = tetrakis(2,6-difluorophenyl)porphyrinate), and a superoxide complex, [(F₈)Fe^III-(O₂^•-)] (S), are shown to be redox interconvertible. Using Cr(η-C₆H₆)₂, an equilibrium state where S and P are present is established in tetrahydrofuran (THF) at -80 °C, allowing determination of the reduction potential of S as -1.17 V vs Fc^+/0. P could be protonated with 2,6-lutidinium triflate, yielding the low-spin ferric hydroperoxide species, [(F₈)Fe^III-(OOH)] (HP). Partial conversion of HP back to P using a derivatized phosphazene base gave a P/HP equilibrium mixture, leading to the determination of pK_a = 28.8 for HP (THF, -80 °C). With the measured reduction potential and pK_a, the O-H bond dissociation free energy (BDFE) of hydroperoxide species HP was calculated to be 73.5 kcal/mol, employing the thermodynamic square scheme and Bordwell relationship. This calculated O-H BDFE of HP, in fact, lines up with an experimental demonstration of the oxidizing ability of S via hydrogen atom transfer (HAT) from TEMPO-H (2,2,6,6-tetramethylpiperdine-N-hydroxide, BDFE = 66.5 kcal/mol in THF), forming the hydroperoxide species HP and TEMPO radical. Kinetic studies carried out with TEMPO-H(D) reveal second-order behavior, k_H = 0.5, k_D = 0.08 M^-1 s^-1 (THF, -80 °C); thus, the hydrogen/deuterium kinetic isotope effect (KIE) = 6, consistent with H-atom abstraction by S being the rate-determining step. This appears to be the first case where experimentally derived thermodynamics lead to a ferric heme hydroperoxide OO-H BDFE determination, that Fe^III-OOH species being formed via HAT reactivity of the partner ferric heme superoxide complex.

Modeling Tryptophan/Indoleamine 2,3-Dioxygenase with Heme Superoxide Mimics: Is Ferryl the Key Intermediate?

Tryptophan oxidation in biology has been recently implicated in a vast array of paramount pathogenic conditions in humans, including multiple sclerosis, rheumatoid arthritis, type-I diabetes, and cancer. This 2,3-dioxygenative cleavage of the indole ring of tryptophan with dioxygen is mediated by two heme enzymes, tryptophan 2,3-dioxygenase (TDO) and indoleamine 2,3-dioxygenase (IDO), during its conversion to N-formylkynurenine in the first and rate-limiting step of kynurenine pathway. Despite the pivotal significance of this enzymatic transformation, a vivid viewpoint of the precise mechanistic events is far from complete. A heme superoxide adduct is thought to be the active oxidant in both TDO and IDO, which, following O-O bond cleavage, presumably generates a key ferryl (Fe^IV=O) reaction intermediate. This study, for the first time in model chemistry, demonstrates the potential of synthetic heme superoxide adducts to mimic the bioinorganic chemistry of indole dioxygenation by TDO and IDO, challenging the widely accepted categorization of these metal adducts as weak oxidants. Herein, an electronically divergent series of ferric heme superoxo oxidants mediates the facile conversion of an array of indole substrates into their corresponding 2,3-dioxygenated products, while shedding light on an unequivocally occurring, putative ferryl intermediate. The oxygenated indole products have been isolated in ∼31% yield, and characterized by LC-MS, ¹H and ¹³C NMR, and FT-IR methodologies, as well as by ¹⁸O_2(g) labeling experiments. Distinctly, the most electron-deficient superoxo adduct is observed to react the fastest, specifically with the most electron-rich indole substrate, underscoring the cruciality of electrophilicity of the heme superoxide moiety in facilitating the initial indole activation step. Comprehensive understanding of such mechanistic subtleties will benefit future attempts in the rational design of salient therapeutic agents, including next generation anticancer drug targets with amplified effectivity.

Crystallographic evidence for unexpected selective tyrosine hydroxylations in an aerated achiral Ru-papain conjugate

The X-ray structure of an aerated achiral Ru-papain conjugate has revealed the hydroxylation of two tyrosine residues found near the ruthenium ion. The most likely mechanism involves a ruthenium-bound superoxide as the reactive species responsible for the first hydroxylation and the resulting high valent Ru(iv)[double bond, length as m-dash]O species for the second one.

Stereoselective Metabolism of Omeprazole by Cytochrome P450 2C19 and 3A4: Mechanistic Insights from DFT Study

The efficacy of S-omeprazole as a proton pump inhibitor compared with that of its enantiomer R-omeprazole is studied using density functional theoretical calculations. The pharmacokinetic studies suggest that the efficacy of S-omeprazole presumably depends on metabolic pathway and excretion from the human body. The density functional theory calculations at SMD_water-B3LYP-D3/6-311+G(d,p)/LANL2DZ//B3LYP/6-31G(d)/LANL2DZ with triradicaloid model active species, [Por^•+Fe^IV(SH)O], of CYP2C19 enzyme with high-spin quartet and low-spin doublet states demonstrate C-H bond activation mechanism through a two-state rebound process for the hydroxylation of R-omeprazole and S-omeprazole. The calculated activation free energy barriers for the hydrogen abstraction are 15.7 and 17.5 kcal/mol for R-omeprazole and S-omeprazole, respectively. The hydroxylation of R-omeprazole and S-omeprazole is thermodynamically favored; however, the hydroxylated intermediate of S-omeprazole further disintegrates to metabolite 5- O-desmethylomeprazole with a higher kinetic barrier. We have examined the sulfoxidation of S-omeprazole to omeprazole sulfone metabolite by CYP3A4, and the observed activation free energy barrier is 9.9 kcal/mol. The computational results reveal that CYP2C19 exclusively metabolizes R-omeprazole to hydroxyomeprazole, which is hydrophilic and can easily excrete, whereas CYP3A4 metabolizes S-omeprazole to lipophilic sulfone; hence, the excretion of this metabolite would be relatively slower from the body. The spin density analysis and molecular orbital analysis performed using biorthogonalization calculations indicate that R-omeprazole favors high-spin pathway for metabolism process whereas S-omeprazole prefers the low-spin pathway.

A theoretical study on the oxidation of alkenes to aldehydes catalyzed by ruthenium porphyrins using O2 as the sole oxidant

Density functional theory (DFT) calculations were used to study the ruthenium porphyrin-catalyzed oxidation of styrene to generate an aldehyde. The results indicate that two reactive oxidants, dioxoruthenium and monooxoruthenium-superoxo porphyrins, participate in the catalytic oxidation. In the mechanism, the resultant monooxoruthenium porphyrin acts in the tandem epoxide isomerization (E-I) to selectively yield an aldehyde and generate a dioxoruthenium porphyrin, thereby triggering new oxidation reaction cycles. In this calculation, several key elements responsible for the observed oxidative ability have been established by using Frontier molecular orbital (FMO) theory, natural bond orbital (NBO) analysis, etc., which include the reaction energy, the spin exchange effect, the spin-state conversion process, and the energy level of the lowest unoccupied molecular orbitals (LUMOs) of the reactive oxidants. The comparative oxidative abilities of the ruthenium-oxo/superoxo compounds with different axial ligands are also investigated. The results suggest that the ruthenium-oxo/superoxo species featuring a chlorine axial ligand is more reactive than that substituted with oxygen. This tuneable reactivity can be understood when considering the different electronic characters of the two ligands and the effective atomic number rule (EAN).

A siloxyl bis-pocket thiolate-tailed Fe(III) porphyrin complex

We report the synthesis and characterization of a bis-pocket iron(III) porphyrin complex that has a covalently attached siloxy thiolate group as the axial ligand. As a new cytochrome P450 model compound, this siloxyl bis-pocket thiolate-tailed iron(III) porphyrin is easily synthesized and is surprisingly stable under air due to the extreme steric protection of the siloxy pockets on both faces of the porphyrin. A single-crystal XRD structure has been determined; the Fe–S bond distance is 2.237 (7) Å with Fe–N bond distances of 2.100 (8) and the Fe is 0.5 Å out of the mean porphyrin plane. The Fe–S bond distance in the siloxy thiolate-tailed iron(III) porphyrin is very similar to that in cytochrome P450 and this structure represents a very rare crystallographically-characterized five-coordinate high spin alkylthiolate-tailed ferric porphyrin. EPR spectrum of this compound showed g values and an E/D ratio very similar to those of cytochrome P450 (CYP450) enzyme, demonstrating the importance of using a very basic thiolate group as the axial ligand in P450 model studies. UV-vis studies of its reduced form with carbon monoxide shows a hyper spectrum, which is characteristic of carbonyl complexes of Fe(II)porphyrin thiolates.

Hydrogen Abstraction of Camphor Catalyzed by Cytochrome P450cam: A QM/MM Study

A combined quantum mechanics and molecular mechanics (QM/MM, QM = UB3LYP-D3, MM = AMBER) method is used to study the hydrogen abstraction reaction in P450_cam catalyzed hydroxylation of camphor in the quartet state. Compared to QM/MM calculations in the literature, this study uses larger basis sets for the most important atoms at the active site and QM/MM Hessian harmonic frequency calculations to determine the standard Gibbs free energy of activation and kinetic isotope effect. The QM/MM covalent boundary is treated with a capping hydrogen atom method, which is simple and robust. An energy barrier of 21.3 kcal/mol and a standard free energy of activation of 16.8 kcal/mol are obtained for this hydrogen abstraction reaction. These values are similar to those reported in the literature, suggesting that when a general protocol is followed, QM/MM results are reproducible. It is found that using a sufficiently large basis set is important to minimize basis set errors.

From DNA to catalysis: a thymine-acetate ligated non-heme iron(III) catalyst for oxidative activation of aliphatic C-H bonds

A non-heme, iron(III)/THA(thymine-1-acetate) catalyst together with H2O2 as an oxidant is efficient in oxidative C-H activation of alkanes. Although having a higher preference for tertiary C-H bonds, the catalyst also oxidizes aliphatic secondary C-H bonds into carbonyl compounds with good to excellent conversions. Based on the site selectivity of the catalyst and our mechanistic studies the reaction proceeds via an Fe-oxo species without long lived carbon centered radicals.

Effect of Natural Polyphenols on CYP Metabolism: Implications for Diseases

Cytochromes P450 (CYPs) are a large group of hemeproteins located on mitochondrial membranes or the endoplasmic reticulum. They play a crucial role in the metabolism of endogenous and exogenous molecules. The activity of CYP is associated with a number of factors including redox potential, protein conformation, the accessibility of the active site by substrates, and others. This activity may be potentially modulated by a variety of small molecules. Extensive experimental data collected over the past decade point at the active role of natural polyphenols in modulating the catalytic activity of CYP. Polyphenols are widespread micronutrients present in human diets of plant origin and in medicinal herbs. These compounds may alter the activity of CYP either via direct interactions with the enzymes or by affecting CYP gene expression. The polyphenol-CYP interactions may significantly alter the pharmacokinetics of drugs and thus influence the effectiveness of chemical therapies used in the treatment of different types of cancers, diabetes, obesity, and cardiovascular diseases (CVD). CYPs are involved in the oxidation and activation of external carcinogenic agents, in which case the inhibition of the CYP activity is beneficial for health. CYPs also support detoxification processes. In this case, it is the upregulation of CYP genes that would be favorable for the organism. A CYP enzyme aromatase catalyzes the formation of estrone and estradiol from their precursors. CYPs also catalyze multiple reactions leading to the oxidation of estrogen. Estrogen signaling and oxidative metabolism of estrogen are associated with the development of cancer. Thus, polyphenol-mediated modulation of the CYP's activity also plays a vital role in estrogen carcinogenesis. The aim of the present review is to summarize the data collected over the last five to six years on the following topics: (1) the mechanisms of the interactions of CYP with food constituents that occur via the direct binding of polyphenols to the enzymes and (2) the mechanisms of the regulation of CYP gene expression mediated by polyphenols. The structure-activity relationship relevant to the ability of polyphenols to affect the activity of CYP is analyzed. The application of polyphenol-CYP interactions to diseases is discussed.

Computational biotransformation profile of paracetamol catalyzed by cytochrome P450

The P450-catalyzed biotransformation of the analgesic drug paracetamol (PAR) is a long-debated topic, involving different mechanistic hypotheses as well as experimental evidence for the metabolites N-acetyl-p-benzoquinone imine (NAPQI), p-benzoquinone, acetamide, and 3-hydroxy-PAR. During the catalytic cycle of P450, a high-valent iron(IV)-oxo species known as Compound I (Cpd I) is formed as the ultimate oxidant, featuring two energetically close-lying ground states in the doublet (low-spin) and quartet (high-spin) spin states, respectively. In order to clarify the catalytic mechanism, a computational chemistry analysis has been undertaken for both the high- and low-spin routes, employing density functional theory (DFT) including PCM (polarized continuum-solvation model) that yields an approximate simulation of the bulk polarization exerted through the protein. The results demonstrate that hydrogen abstraction transfer (HAT) by the P450 oxidant Cpd I (FeO) is kinetically strongly preferred over the alternative pathways of an oxygen addition reaction (OAR) or two consecutive single-electron transfers (SET). Moreover, only the respective high-spin route yields N-acetyl-p-semiquinone imine (NAPSQI) as an intermediate that is converted to the electrophile N-acetyl-p-benzoquinone imine (NAPQI). By contrast, 3-hydroxy-PAR, acetamide, and p-benzoquinone as electrophilic and redox-active agent are formed predominantly in the low-spin state through reactions that do not involve NAPSQI. Thus, all experimentally observed PAR metabolites are in accord with an initial HAT from the phenolic oxygen, and NAPSQI should indeed be the precursor of NAPQI, both of which are generated only via the high-spin pathway.

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

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

QM/MM study of the C-C coupling reaction mechanism of CYP121, an essential cytochrome p450 of Mycobacterium tuberculosis

Among 20 p450s of Mycobacterium tuberculosis (Mt), CYP121 has received an outstanding interest, not only due to its essentiality for bacterial viability but also because it catalyzes an unusual carbon-carbon coupling reaction. Based on the structure of the substrate bound enzyme, several reaction mechanisms were proposed involving first Tyr radical formation, second Tyr radical formation, and C-C coupling. Key and unknown features, being the nature of the species that generate the first and second radicals, and the role played by the protein scaffold each step. In the present work we have used classical and quantum based computer simulation methods to study in detail its reaction mechanism. Our results show that substrate binding promotes formation of the initial oxy complex, Compound I is the responsible for first Tyr radical formation, and that the second Tyr radical is formed subsequently, through a PCET reaction, promoted by the presence of key residue Arg386. The final C-C coupling reaction possibly occurs in bulk solution, thus yielding the product in one oxygen reduction cycle. Our results thus contribute to a better comprehension of MtCYP121 reaction mechanism, with direct implications for inhibitor design, and also contribute to our general understanding of these type of enzymes.

Heme-dependent dioxygenases in tryptophan oxidation

L-Tryptophan is an essential amino acid for mammals. It is utilized not only for protein synthesis but also for the biosynthesis of serotonin and melatonin by the serotonin pathway as well as nicotinamide adenine dinucleotide by the kynurenine pathway. Although the kynurenine pathway is responsible for the catabolism of over 90% of l-tryptophan in the mammalian intracellular and extracellular pools, the scientific field was dominated in the last century by studies of the serotonin pathway, due to the physiological significance of the latter's catabolic intermediates and products. However, in the past decade, the focus gradually reversed as the link between the kynurenine pathway and various neurodegenerative disorders and immune diseases is increasingly highlighted. Notably, the first step of this pathway, which is catalyzed by heme-dependent dioxygenases, has been proven to be a potential target for immune regulation and cancer treatment. A thorough understanding of the intriguing chemistry of the heme-dependent dioxygenases may yield insight for the drug discovery of these prevalent illnesses. In this review, we survey enzymatic and mechanistic studies, initially started by Kotake and Masayama over 70 years ago, at the molecular level on the heme-dependent tryptophan dioxygenation reactions.

A dominant homolytic O-Cl bond cleavage with low-spin triplet-state Fe(IV)=O formed is revealed in the mechanism of heme-dependent chlorite dismutase

Chlorite dismutase (Cld) is a heme-dependent enzyme that catalyzes the decomposition of toxic chlorite (ClO2(-)) into innocuous chloride and O2. In this paper, using the hybrid B3LYP density functional theory (DFT) method including dispersion interactions, the Cld reaction mechanism has been studied with a chemical model constructed on the X-ray crystal structure. The calculations indicate that the reaction proceeds along a stepwise pathway in the doublet state, i.e. a homolytic O-Cl bond cleavage of the substrate leading to an O-Fe(heme) species and a ClO˙ radical, followed by a rebinding O-O bond formation between them. The O-Fe(heme) species is demonstrated to be a low-spin triplet-state Fe(IV)=O diradicaloid. A low-spin singlet-state Fe(IV)=O is much less stable than the former, with an energy difference of 9.2 kcal mol(-1). The O-Cl bond cleavage is rate-limiting with a barrier of 10.6 kcal mol(-1), in good agreement with the experimental reaction rate of 2.0 × 10(5) s(-1). Furthermore, a heterolytic O-Cl bond dissociation in the initial step is shown to be unreachable, which ensures the high efficiency of the Cld enzyme by avoiding the generation of chlorate byproduct observed in the reactions of synthetic Fe porphyrins. Also, the pathways in the quartet and sextet states are unfavorable for the Cld reaction. The present results reveal a detailed mechanism III (defined in the text) including an interesting di-radical intermediate composed of a low-spin triplet-state Fe(IV)=O and a ClO˙ radical. Compared to a competitive heterolytic Cl-O cleavage in synthetic Fe porphyrins, the revelation of the domination of homolysis in Cld indicates not only the high efficiency of enzyme, but also the sensitivity of a heme and the significance of the enzymatic active-site surroundings (the His170 and Arg183 residues in the present case), which gives more insights into heme chemistry.

Hydrogen peroxide as oxidant in bio-mimetic catalysis by manganese porphyrin: Theoretical DFT studies

The aim of this study is to elucidate the geometry and electronic structure of various adducts that may be formed between manganese(III) (Mn(III)) porphyrin and hydrogen peroxide. Hydrogen peroxide may interact with Mn(III) porphyrin either as H2O2 or, after dissociation, as OOH–. In the former, it may decompose into two hydroxo groups, which acquire OH– character or an oxo group (=O) and a water molecule. Therefore, the following systems are considered: MnP(H2O2)+, MnP(H2O2)(OH), MnP(OH)3, [Formula: see text], MnPO+, MnPO(OH), MnP(OOH), MnP(OOH)(OH)–, and the possible transformations between them are taken into account. The reported studies are performed within the Density Functional Theory (DFT) method with the GGA-BP functional. The geometry and electronic structures of the structures found along the studied reaction pathways are discussed in terms of interatomic distances, valence angles, Mulliken charges, and spin densities. It was found that different active oxygen species may be formed in the reaction between Mn(III) porphyrin and hydrogen peroxide. As manganese is a transition metal, numerous possible spin states for each of the studied structures are found, where the relative energies of different multiplicities depend strongly on the ligands present in the complex. In view of the catalytic properties, all oxygen-containing ligands are negatively charged, which results in their behaviour as nucleophiles towards hydrocarbons. Finally, the analysis of charge and spin populations on different parts of the studied systems indicate the porphyrin ligand as active in charge transfer processes.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

The use of deuterated camphor as a substrate in (1)H ENDOR studies of hydroxylation by cryoreduced oxy P450cam provides new evidence of the involvement of compound I

Electron paramagnetic resonance and (1)H electron nuclear double resonance (ENDOR) spectroscopies have been used to analyze intermediate states formed during the hydroxylation of (1R)-camphor (H(2)-camphor) and (1R)-5,5-dideuterocamphor (D(2)-camphor) as induced by cryoreduction (77 K) and annealing of the ternary ferrous cytochrome P450cam-O(2)-substrate complex. Hydroxylation of H(2)-camphor produced a primary product state in which 5-exo-hydroxycamphor is coordinated with Fe(III). ENDOR spectra contained signals derived from two protons [Fe(III)-bound C5-OH(exo) and C5-H(endo)] from camphor. When D(2)-camphor was hydroxylated under the same condition in H(2)O or D(2)O buffer, both ENDOR H(exo) and H(endo) signals are absent. For D(2)-camphor in H(2)O buffer, H/D exchange causes the C5-OH(exo) signal to reappear during relaxation upon annealing to 230 K; for H(2)-camphor in D(2)O, the magnitude of the C5-OH(exo) signal decreases via H/D exchange. These observations clearly show that Compound I is the reactive species in the hydroxylation of camphor in P450cam.

Model iron-oxo species and the oxidation of imidazole: insights into the mechanism of OvoA and EgtB?

A density functional theory cluster and first-principles quantum and statistical mechanics approach have been used to investigate the ability of iron-oxygen intermediates to oxidize a histidine cosubstrate, which may then allow for the possible formation of 2- and 5-histidylcysteine sulfoxide, respectively. Namely, the ability of ferric superoxo (Fe(III)O(2)(•-)), Fe(IV)═O, and ferrous peroxysulfur (Fe(III)OOS) complexes to oxidize the imidazole of histidine via an electron transfer (ET) or a proton-coupled electron transfer (PCET) was considered. While the high-valent mononuclear Fe(IV)═O species is generally considered the ultimate biooxidant, the free energies for its reduction (via ET or PCET) suggest that it is unable to directly oxidize histidine's imidazole. Instead, only the ferrous peroxysulfur complexes are sufficiently powerful enough oxidants to generate a histidyl-derived radical via a PCET process. Furthermore, while this process preferably forms a HisN(δ)(-H)(•) radical, several such oxidants are also suggested to be capable of generating the higher-energy HisC(δ)(-H)(•) and HisC(ε)(-H)(•) radicals. Importantly, the present results suggest that formation of the sulfoxide-containing products (seen in both OvoA and EgtB) is a consequence of the reduction of a powerful Fe(III)OOS oxidant via a PCET.

Axial and equatorial ligand effects on biomimetic cysteine dioxygenase model complexes

Density functional theory (DFT) calculations are presented on biomimetic model complexes of cysteine dioxygenase and focus on the effect of axial and equatorial ligand placement. Recent studies by one of us [Y. M. Badiei, M. A. Siegler and D. P. Goldberg, J. Am. Chem. Soc. 2011, 133, 1274] gave evidence of a nonheme iron biomimetic model of cysteine dioxygenase using an i-propyl-bis(imino)pyridine, equatorial tridentate ligand. Addition of thiophenol, an anion - either chloride or triflate - and molecular oxygen, led to several possible stereoisomers of this cysteine dioxygenase biomimetic complex. Moreover, large differences in reactivity using chloride as compared to triflate as the binding anion were observed. Here we present a series of DFT calculations on the origin of these reactivity differences and show that it is caused by the preference of coordination site of anion versus thiophenol binding to the chemical system. Thus, stereochemical interactions of triflate and the bulky iso-propyl substituents of the ligand prevent binding of thiophenol in the trans position using triflate. By contrast, smaller anions, such as chloride, can bind in either cis or trans ligand positions and give isomers with similar stability. Our calculations help to explain the observance of thiophenol dioxygenation by this biomimetic system and gives details of the reactivity differences of ligated chloride versus triflate.

1,2,3-triazole: unique ligand in promoting iron-catalyzed propargyl alcohol dehydration

A 1,2,3-traizole-promoted iron(III)-catalyzed propargyl alcohol dehydration was developed for the synthesis of conjugated enynes. The desired conjugated enynes were prepared in good to excellent yields (up to 95%) with a large substrate scope and excellent stereoselectivity (only Z-isomers).

Mechanism of S-oxygenation by a cysteine dioxygenase model complex

In this work, we present the first computational study on a biomimetic cysteine dioxygenase model complex, [Fe(II)(LN(3)S)](+), in which LN(3)S is a tetradentate ligand with a bis(imino)pyridyl scaffold and a pendant arylthiolate group. The reaction mechanism of sulfur dioxygenation with O(2) was examined by density functional theory (DFT) methods and compared with results obtained for cysteine dioxygenase. The reaction proceeds via multistate reactivity patterns on competing singlet, triplet, and quintet spin state surfaces. The reaction mechanism is analogous to that found for cysteine dioxygenase enzymes (Kumar, D.; Thiel, W.; de Visser, S. P. J. Am. Chem. Soc. 2011, 133, 3869-3882); hence, the computations indicate that this complex can closely mimic the enzymatic process. The catalytic mechanism starts from an iron(III)-superoxo complex and the attack of the terminal oxygen atom of the superoxo group on the sulfur atom of the ligand. Subsequently, the dioxygen bond breaks to form an iron(IV)-oxo complex with a bound sulfenato group. After reorganization, the second oxygen atom is transferred to the substrate to give a sulfinic acid product. An alternative mechanism involving the direct attack of dioxygen on the sulfur, without involving any iron-oxygen intermediates, was also examined. Importantly, a significant energetic preference for dioxygen coordinating to the iron center prior to attack at sulfur was discovered and serves to elucidate the function of the metal ion in the reaction process. The computational results are in good agreement with experimental observations, and the differences and similarities of the biomimetic complex and the enzymatic cysteine dioxygenase center are highlighted.

Comparative reactivity of ferric-superoxo and ferryl-oxo species in heme and non-heme complexes

Ferryl-oxo species have been recognized as a key oxidant in many heme and non-heme enzymes. Recently, less-characterized ferric-superoxo species have been found or suggested to be another electrophilic oxidant. Reactivity of several vital ferryl-oxo and ferric-superoxo model complexes was examined by DFT calculations. Reactivity is found to correlate well with thermodynamic driving force and can increase with higher electrophilicity of the oxidant. Reactivity of the ferric-superoxo oxidants generally is not "superior" to the ferryl-oxo ones. Compared to the high-spin non-heme ferric-superoxo, the lower reactivity of low-spin heme ferric-superoxo, seldom utilized in nature, can be attributed to lower electrophilicity and more pronounced quenching of anti-ferromagnetic coupling between the ferric and superoxo parts. The present comparison should shed some light on mechanistic strategies in heme and non-heme enzymes and provide clues to rational design of ferric-superoxo oxidants.