Two States and Two More in the Mechanisms of Hydroxylation and Epoxidation by Cytochrome P450

Computational Modeling of the Enzymatic Achmatowicz Rearrangement by Heme-Dependent Chloroperoxidase: Reaction Mechanism, Enantiopreference, Regioselectivity, and Substrate Specificity

The chloroperoxidase from Caldariomyces fumago (CfCPO) catalyzes the oxidative ring expansion of α-heterofunctionalized furans via the Achmatowicz rearrangement, providing an elegant tool to convert furan rings into complex-prefunctionalized scaffolds. However, the mechanism of this transformation remains unclear. Herein, the CfCPO-catalyzed reaction of rac-1-(2-furyl)ethanol (1a) is studied by quantum chemical calculations and molecular dynamics simulations. The calculations reveal that the conversion follows the general mechanism of the Achmatowicz reaction. Notably, the binding of 1a to the enzyme's active site influences the Compound I (Cpd I) formation, and the (R)-1a enantiomer binding results in a lower barrier compared to (S)-1a, explaining the observed (R)-enantiopreference toward a racemic substrate. Additionally, due to the weaker steric hindrance between the porphyrin ring and substrate, the nucleophilic attack of Cpd I on the furan core of 1a is preferred at the less-substituted C4=C5 bond, providing a rationale for the experimentally observed regioselectivity. Finally, the bottleneck residues in the substrate delivery channel and also the active site surroundings are proposed to be responsible for the substrate specificity of CfCPO. This study lays a theoretical foundation for the rational design of new CPOs that catalyze the Achmatowicz rearrangement with a broader substrate spectrum or specific stereopreference.

Hydrogen-Bond-Assisted Catalysis: Hydroxylation of Paclitaxel by Human CYP2C8

Paclitaxel (PTX, or Taxol), a chemotherapeutic agent widely employed in the treatment of various cancers, undergoes metabolic transformations through the cytochrome P450 enzymes CYP3A4 and CYP2C8. CYP3A4 catalyzes the aromatic hydroxylation reaction of PTX, whereas CYP2C8 demonstrates a distinct reactivity pattern, producing 6α-hydroxypaclitaxel via alkane hydroxylation. Despite the significant impact of PTX metabolism on its anticancer efficacy, the detailed mechanisms underlying these transformations have remained largely unclear. In this study, we employed hybrid quantum mechanics and molecular mechanics (QM/MM) calculations to elucidate the mechanism of PTX metabolism by human CYP2C8. Our QM/MM results reveal that the hydroxylation of PTX by CYP2C8 follows an atypical rebound mechanism. Either of the two hydrogen atoms at the C6 position of PTX can be abstracted, leading to a common radical intermediate. Although the subsequent rebound barrier is unusually high, stereochemical scrambling is unlikely, as the rebound barrier for the formation of the 6α-hydroxylated PTX─the actual product─is significantly lower than that for the 6β-hydroxylated metabolite. Thus, product selectivity is determined by the non-rate-determining rebound step. Furthermore, the hydroxyl group at the C7 position of PTX plays a catalytic role by facilitating the hydrogen abstraction and rebound steps. Our study also confirms a pronounced stability of the transition state in the high-spin sextet spin state, enabled by the enzyme's specific substrate positioning.

Enhancing the high-spin reactivity in C-H bond activation by Iron (IV)-Oxo species: insights from paclitaxel hydroxylation by CYP2C8

Previous theoretical studies have revealed that high-spin states possess flatter potential energy surfaces than low-spin states in reactions involving iron(IV)-oxo species of cytochrome P450 enzymes (P450s), nonheme enzymes, or biomimetic complexes. Therefore, actively utilizing high-spin states to enhance challenging chemical transformations, such as C-H bond activation, represents an intriguing research avenue. However, the inherent instability of high-spin states relative to low-spin states in pre-reaction complexes often hinders their accessibility around the transition state, especially in heme systems with strong ligand fields. Counterintuitively, our investigation of the metabolic hydroxylation of paclitaxel by human CYP2C8 using a hybrid quantum mechanics and molecular mechanics (QM/MM) approach showed that the high-spin sextet state exhibits unusually high stability, when the reaction follows a secondary reaction pathway leading to 6β-hydroxypaclitaxel. We thoroughly analyzed the factors contributing to the enhanced stabilization of the high-spin state, and the knowledge obtained could be instrumental in designing competent biomimetic catalysts and biocatalysts for C-H bond activation.

Exploring the Bonding Nature of Iron(IV)-Oxo Species through Valence Bond Calculations and Electron Density Analysis

Compound I (Cpd I) plays a pivotal role in substrate transformations within the catalytic cycle of cytochrome P450 enzymes (P450s). A key constituent of Cpd I is the iron(IV)-oxo unit, a structural motif also found in other heme enzymes and nonheme enzymes. In this study, we performed ab initio valence bond (VB) calculations, employing the valence bond self-consistent field (VBSCF) and breathing orbital valence bond (BOVB) methods, to unveil the bonding nature of this vital "Fe(IV)═O″ unit in bioinorganic chemistry. Comparisons were drawn with the triplet O2 molecule, which shares some electronic characteristics with iron(IV)-oxo. Additionally, Cpd I models of horseradish peroxidase (HRP) and catalase (CAT) were analyzed to assess the proximal ligand effect on the electronic structure of iron(IV)-oxo. Our VB analysis underscores the significant role of noncovalent resonance effects in shaping the iron(IV)-oxo bonding. The resonance stabilization within the π and σ frameworks occurs to comparable degrees, with additional stabilization resulting from resonance between VB structures from these frameworks. Furthermore, we elucidated the substantial influence of proximal and equatorial ligands in modulating the relative significance of different VB structures. Notably, in the presence of these ligands, iron(IV)-oxo is better described as iron(III)-oxyl or iron(II)-oxygen, displaying weak covalent character but enhanced by resonance effects. Although both species exhibit diradicaloid characters, resonance stabilization in iron(IV)-oxo is weaker than in O2. Further exploration using the Laplacian of electron density shows that, unlike O2, which exhibits a charge concentration region between its two oxygen atoms, iron(IV)-oxo species display a charge depletion region.

Unraveling the Geometry-Driven C═C Epoxidation and C-H Hydroxylation Reactivity of Tetra-Coordinated Nonheme Iron(IV)-Oxo Complexes

The electronic structure and reactivity of tetra-coordinated nonheme iron(IV)-oxo complexes have remained unexplored for years. The recent synthesis of a closed-shell iron(IV)-oxo complex [(quinisox)FeIV(O)]+ (1) has set up a platform to understand how such complexes compare with the celebrated open-shell iron-oxo chemistry. Herein, using density functional theory and ab initio calculations, we present an in-depth electronic structure investigation of the C═C epoxidation [oxygen atom transfer (OAT)] and C-H hydroxylation [hydrogen atom transfer (HAT)] reactivity of 1. Using a solvent-coordinated geometry of 1 (1') and other potential tetra-coordinated iron(IV)-oxo complexes bearing rigid ligands (2 and 3), we established the geometric origin of spin-state energetics and reactivity of 1. Complex 1 featuring a strong Fe-O bond exhibits OAT and HAT reactivity in its quintet state. The lowest quintet OAT pathway has a lower barrier by ∼4 kcal/mol than the quintet HAT pathway, corroborating the experimentally observed gas-phase OAT reactivity preference. A conventional HAT reactivity preference for 2 and a comparable OAT and HAT reactivity for 3 are observed. This further supports the geometry-driven reactivity preference for 1. Noncovalent interaction analyses reveal a pronounced π-π interaction between the substrate and ligand in the OAT transition state, rationalizing the origin of the observed reactivity preference for 1.

Reaction engineering blocks ether cleavage for synthesizing chiral cyclic hemiacetals catalyzed by unspecific peroxygenase

Hemiacetal compounds are valuable building blocks in synthetic chemistry, but their enzymatic synthesis is limited and often hindered by the instability of hemiacetals in aqueous environments. Here, we show that this challenge can be addressed through reaction engineering by using immobilized peroxygenase from Agrocybe aegerita (AaeUPO) under neat reaction conditions, which allows for the selective C-H bond oxyfunctionalization of environmentally significant cyclic ethers to cyclic hemiacetals. A wide range of chiral cyclic hemiacetal products are prepared in >99% enantiomeric excess and 95170 turnover numbers of AaeUPO. Furthermore, by changing the reaction medium from pure organic solvent to alkaline aqueous conditions, cyclic hemiacetals are in situ transformed into lactones. Lactams are obtained under the applied conditions, albeit with low enzyme activity. These findings showcase the synthetic potential of AaeUPO and offer a practical enzymatic approach to produce chiral cyclic hemiacetals through C-H oxyfunctionalization under mild conditions.

Mechanism of Selective Aromatic Hydroxylation in the Metabolic Transformation of Paclitaxel Catalyzed by Human CYP3A4

Paclitaxel (PTX) is heralded as one of the most successful natural-product drugs for the treatment of refractory cancers. In humans, the hepatic metabolic transformation of PTX is primarily mediated by two cytochrome P450 enzymes (P450s): CYP3A4 and CYP2C8. The impact of P450 metabolism on the anticancer effectiveness of PTX is significant. However, the precise mechanism underlying selective P450-catalyzed reactions in PTX metabolism remains elusive. To address this knowledge gap, we conducted molecular docking and molecular dynamics simulations using multiple crystal structures of CYP3A4, which originally contained other ligands. These methods enabled us to determine the most plausible binding structure of PTX within the enzyme. By further employing hybrid quantum mechanics and molecular mechanics calculations, we successfully identified two primary pathways for the reaction between compound I (Cpd I) of CYP3A4 and PTX. One of these pathways involves the formation of an epoxide, while the other proceeds through a ketone intermediate.

The Nature of the Chemical Bonds of High-Valent Transition-Metal Oxo (M=O) and Peroxo (MOO) Compounds: A Historical Perspective of the Metal Oxyl-Radical Character by the Classical to Quantum Computations

This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).

Bidirectional Charge Transfer at the Heme Iron in Reversible and Quasi-irreversible Cytochrome P450 Inhibition

The coordination bonding between inhibitor ligands and heme iron plays a critical role in disrupting the essential catalytic functions of cytochrome P450 enzymes (P450s). Despite its intrinsic importance and consequential implications for human health, our current understanding of coordination bonding in P450 inhibition remains limited. To address this knowledge gap, we conducted a systematic theoretical analysis of the complexes between a ferric or a ferrous heme model and representative inhibitor ligands. Specifically, we evaluated the charge-transfer (CT) effect within these complexes by employing a series of theoretical methods based on density functional theory (DFT). Through a comprehensive analysis, we unveiled the relative significance of ligand-to-heme forward CT in the ferric and ferrous complexes of reversible inhibitors. In contrast, backward CT dominates over forward CT in the ferrous heme complexes of quasi-irreversible inhibitors. Further analysis using the compact frontier orbital method underscores the elevated electron-accepting abilities of quasi-irreversible inhibitors for π backdonation, which greatly amplifies their binding affinity for the ferrous heme. This study sheds light on the intricate mechanisms underlying P450 inhibition and provides valuable insights for future inhibitor design and development.

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.

The role of agostic interaction in the mechanism of ethylene polymerisation using Cr(III) half-sandwich complexes: What dictates the reactivity?

Chromium-based catalysts play a significant role in the production of ultra-high molecular weight polyethylene, and half-sandwich functionalised-metallocene complexes were proven to be one of the most suitable candidates as catalysts for generating such large polymeric-length olefins. Earlier experimental studies on olefin polymerisation using a series of catalysts such as [L1-2CrCl2] (where L1 = 1-((pyridin-2-yl)methyl)indenyl (1) and L2 = 2-methyl-1-{[4-(yridinene-1-yl)yridine-2-yl]methyl}-1H-indenyl (2)) reveal significant variation where peripheral substitution on the ligand was found to influence the reactivity significantly. However, the specific ligand position that affects the reactivity has not been established. As these reactions are fast and robust, it is challenging to establish reactive intermediates via experiments, and therefore, mechanistic clues for such reactions are elusive. Here we have undertaken a detailed computational study by employing an array of DFT (uB3LYP-D3/def2-TZVP, CASSCF/NEVPT2, and DLPNO-CCSD(T) methods to explore the substituted and non-substituted pyridine-cyclo-pentadienyl chromium complexes and their influence on the catalytic activity in ethylene polymerisation. Our study not only unravels the catalytic pathway for olefin polymerisation for such Cr(III)-half-sandwich complexes but also reveals that the energetics of the formation of pseudo-three-coordinate alkyl intermediates is key to the variation in the reactivity observed. A detailed examination using MO and NBO analysis unveils the presence of a C-H⋯Cr agostic interaction that is found to significantly stabilise this intermediate when the pyridine ligand has strong electron-donating groups at its para position. The other substitutions, such as on the cyclopentadienyl ligand, neither yield the desired stability nor the desired interaction. Further studies on models support this proposal. In order to improve the efficiency and selectivity of catalytic systems in olefin polymerisation, we strongly advocate for the integration of agostic interactions as a crucial criterion in the design of future catalysts. Considering the prevalence of electron-deficient metal centres in successful olefin polymerisation catalysts, this research prompts a broader mechanistic inquiry to propose a unified approach for this industrially crucial reaction.

The role of basicity in selective C-H bond activation by transition metal-oxidos

The development of (bio)catalysts capable of selectively activating strong C-H bonds is a continuing challenge in modern chemistry. In both metalloenzymes and synthetic systems capable of activating C-H bonds, transition metal-oxido intermediates serve as the active species for reactivity whose thermodynamic properties influence the bond strengths they are capable of activating. In this Frontier article, we present current ideas of how the basicity of transition metal-oxidos impacts their reactivity with C-H bonds and present new opportunities within this field. We highlight recent insights into the role basicity plays in the activation process and its influence on mechanism, as well as the important role that secondary coordination sphere effects, such as hydrogen bonds, in tuning the basicity of the metal-oxido species is discussed.

C-N Bond Forming Radical Rebound Is the Enantioselectivity-Determining Step in P411-Catalyzed Enantioselective C(sp3)-H Amination: A Combined Computational and Experimental Investigation

Engineered metalloenzymes represent promising catalysts for stereoselective C-H functionalization reactions. Recently, P450 enzymes have been evolved to allow for new-to-nature intramolecular C(sp3)-H amination reactions via a nitrene transfer mechanism, giving rise to diamine derivatives with excellent enantiocontrol. To shed light on the origin of enantioselectivity, a combined computational and experimental study was carried out. Hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the activation energies and enantioselectivities of both the hydrogen atom transfer (HAT) and the subsequent C-N bond forming radical rebound steps. Contrary to previously hypothesized enantioinduction mechanisms, our calculations show that the radical rebound step is enantioselectivity-determining, whereas the preceding HAT step is only moderately stereoselective. Furthermore, the selectivity in the initial HAT is ablated by rapid conformational change of the radical intermediate prior to C-N bond formation. This finding is corroborated by our experimental study using a set of enantiomerically pure, monodeuterated substrates. Furthermore, classical and ab initio molecular dynamics simulations were carried out to investigate the conformational flexibility of the carbon-centered radical intermediate. This key radical species undergoes a facile conformational change in the enzyme active site from the pro-(R) to the pro-(S) configuration, whereas the radical rebound is slower due to the spin-state change and ring strain of the cyclization process, thereby allowing stereoablative C-N bond formation. Together, these studies revealed an underappreciated enantioinduction mechanism in biocatalytic C(sp3)-H functionalizations involving radical intermediates, opening up new avenues for the development of other challenging asymmetric C(sp3)-H functionalizations.

Assessing Alkene Reactivity toward Cytochrome P450-Mediated Epoxidation through Localized Descriptors and Regression Modeling

The prediction of sites of epoxidation by cytochrome P450s during metabolism is particularly important in drug design, as epoxides are capable of alkylating biological macromolecules. Reliable methods are needed to quantitatively predict P450-mediated epoxidation barriers for inclusion in high-throughput screening campaigns alongside protein-ligand docking. Utilizing the fractional occupation number weighted density (FOD) and orbital-weighted Fukui index (f_w⁺) as descriptors of local reactivity and a data set of 36 alkene epoxidation barriers computed with density functional theory (DFT), we developed and validated a multiple linear regression model for the reliable estimation of epoxidation barriers using only substrate structures as input. Using our recommended level of theory (GFN2-xTB//GFN-FF), mean absolute errors in the training and test sets were found to be 0.66 and 0.70 kcal/mol, respectively, with coefficients of determination of ca. 0.80. We demonstrate the utility of this approach on three known substrates of CYP101A1 and further show that this approach is inappropriate for particularly electron-rich alkenes. By employing a modern semiempirical method on force-field-generated geometries, the required descriptors can be calculated on the millisecond timescale per structure, making the approach well suited for incorporation into high-throughput methodologies alongside docking.

Gold nanoclusters as prospective carriers and detectors of pramipexole

Pramipexole (PPX) is known in the treatment of Parkinson's disease and restless legs syndrome. We carried out a theoretical investigation on pramipexole-Au cluster interactions for the applications of drug delivery and detection. Three Au N clusters with sizes N = 6, 8 and 20 were used as reactant models to simulate the metallic nanostructured surfaces. Quantum chemical computations were performed in both gas phase and aqueous environments using density functional theory (DFT) with the PBE functional and the cc-pVDZ-PP/cc-pVTZ basis set. The PPX drug is mainly adsorbed on gold clusters via its nitrogen atom of the thiazole ring with binding energies of ca. -22 to -28 kcal mol-1 in vacuum and ca. -18 to -24 kcal mol-1 in aqueous solution. In addition to such Au-N covalent bonding, the metal-drug interactions are further stabilized by electrostatic effects, namely hydrogen-bond NH⋯Au contributions. Surface-enhanced Raman scattering (SERS) of PPX adsorbed on the Au surfaces and its desorption process were also examined. In comparison to Au8, both Au6 and Au20 clusters undergo a shorter recovery time and a larger change of energy gap, being possibly conducive to electrical conversion, thus signaling for detection of the drug. A chemical enhancement mechanism for SERS procedure was again established in view of the formation of nonconventional hydrogen interactions Au⋯H-N. The binding of PPX to a gold cluster is expected to be reversible and triggered by the presence of cysteine residues in protein matrices or lower-shifted alteration of environment pH. These findings would encourage either further theoretical probes to reach more accurate views on the efficiency of pramipexole-Au interactions, or experimental attempts to build appropriate gold nanostructures for practical trials, harnessing their potentiality for applications.

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.

Copper catalysis with redox-active ligands

Copper catalysis finds applications in various synthetic fields by utilizing the ability of copper to sustain mono- and bielectronic elementary steps. Further to the development of well-defined copper complexes with classical ligands such as phosphines and N-heterocyclic carbenes, a new and fast-expanding area of research is exploring the possibility of a complementing metal-centered reactivity with electronic participation by the coordination sphere. To achieve this electronic flexibility, redox-active ligands can be used to engage in a fruitful “electronic dialogue” with the metal center, and provide additional venues for electron transfer. This review aims to present the latest results in the area of copper-based cooperative catalysis with redox-active ligands.

The oxidation of cyclo-olefin by the S = 2 ground-state complex [FeIV(O)(TQA)(NCMe)]2

Density functional theory calculation is used to investigate the oxidation of cyclo-olefin (cyclobutene, cyclopentene, cyclohexene, cycloheptene, and cyclo-octene) by the complex [Fe^IV(O)(TQA)(NCMe)]²⁺, which has S = 2 ground state, and the effect of electronic factors and steric hindrance on reaction barriers. Our results suggest that the oxo-iron(IV) complex can oxidise C-H and C = C bonds via a single-state mechanism, and two different ways of electron transport exist. The energy barriers initially decrease with increasing substrate size, and the trend then reverses. Comparison of the energy barrier in different systems reveals that except for the reaction between [Fe^IV(O)(TQA)(NCMe)]²⁺ and cycloheptene, oxo-iron(IV) complexes prefer epoxidation to hydroxylation. However, the hydroxylated product is more stable than the corresponding epoxidated product. This result indicates that the products of epoxidation tend to decompose first. The energy barrier of hydroxylation and epoxidation originates from the balance of orbital interaction and Pauli repulsion from the equatorial ligand and protons on the approaching substrate. In this regard, we calculate the weak interaction between two fragments (oxo-iron complex and substrates) using the independent gradient model and drawn the corresponding 3D isosurface representations of reactants.

Mechanistic insights into the allylic oxidation of aliphatic compounds by tetraamido iron(v) species: A C–H vs. O–H bond activation

This work is based on a deep insight into a comparative study of C–H vs. O–H bond activation of allylic compound by the high valent iron complex. Our theoretical findings can help to design catalysts with better efficiency for catalytic reactions.

Mechanistic dichotomies in redox reactions of mononuclear metal–oxygen intermediates

This review article focuses on various mechanistic dichotomies in redox reactions of metal–oxygen intermediates with the emphasis on understanding and controlling their redox reactivity from experimental and theoretical points of view.

Computational Exploration of Chiral Iron Porphyrin-Catalyzed Asymmetric Hydroxylation of Ethylbenzene Where Stereoselectivity Arises from π-π Stacking Interaction

The mechanism and origins of stereoselectivity of chiral iron porphyrin-catalyzed asymmetric hydroxylation of ethylbenzene were explored with density functional theory. The hydrogen atom abstraction is the rate- and stereoselectivity-determining step. In good agreement with experimental results, the formation of the (R)-1-phenylethanol product is found to be the most favorable pathway. The transition state of hydrogen atom abstraction which leads to the (S)-1-phenylethanol product is unfavorable by 1.7 kcal/mol compared to the corresponding transition state which leads to the (R)-1-phenylethanol product. Enantioselectivity arises from an attractive π-π stacking interaction between the phenyl group of ethylbenzene substrate and the naphthyl group of the porphyrin ligand.

Chemoselectivity in the Oxidation of Cycloalkenes with a Non-Heme Iron(IV)-Oxo-Chloride Complex: Epoxidation vs. Hydroxylation Selectivity

We report and analyze chemoselectivity in the gas phase reactions of cycloalkenes (cyclohexene, cycloheptene, cis-cyclooctene, 1,4-cyclohexadiene) with a non-heme iron(IV)-oxo complex [(PyTACN)Fe(O)(Cl)]⁺, which models the active species in iron-dependent halogenases. Unlike in the halogenases, we did not observe any chlorination of the substrate. However, we observed two other reaction pathways: allylic hydrogen atom transfer (HAT) and alkene epoxidation. The HAT is clearly preferred in the case of 1,4-cyclohexadiene, both pathways have comparable reaction rates in reaction with cyclohexene, and epoxidation is strongly favored in reactions with cycloheptene and cis-cyclooctene. This preference for epoxidation differs from the reactivity of iron(IV)-oxo complexes in the condensed phase, where HAT usually prevails. To understand the observed selectivity, we analyze effects of the substrate, spin state, and solvation. Our DFT and CASPT2 calculations suggest that all the reactions occur on the quintet potential energy surface. The DFT-calculated energies of the transition states for the epoxidation and hydroxylation pathways explain the observed chemoselectivity. The SMD implicit solvation model predicts the relative increase of the epoxidation barriers with solvent polarity, which explains the clear preference of HAT in the condensed phase.

Metabolic Activation and Carcinogenesis of Tobacco-Specific Nitrosamine N'-Nitrosonornicotine (NNN): A Density Function Theory and Molecular Docking Study

N'-nitrosonornicotine (NNN) is one of the tobacco-specific nitrosamines (TSNAs) that exists widely in smoke and smokeless tobacco products. NNN can induce tumors in various laboratory animal models and has been identified by International Agency for Research on Cancer (IARC) as a human carcinogen. Metabolic activation of NNN is primarily initiated by cytochrome P450 enzymes (CYP450s) via 2'-hydroxylation or 5'-hydroxylation. Subsequently, the hydroxylating intermediates undergo spontaneous decomposition to generate diazohydroxides, which can be further converted to alkyldiazonium ions, followed by attacking DNA to form various DNA damages, such as pyridyloxobutyl (POB)-DNA adducts and pyridyl-N-pyrrolidinyl (py-py)-DNA adducts. If not repaired correctly, these lesions would lead to tumor formation. In the present study, we performed density functional theory (DFT) computations and molecular docking studies to understand the mechanism of metabolic activation and carcinogenesis of NNN. DFT calculations were performed to explore the 2'- or 5'- hydroxylation reaction of (R)-NNN and (S)-NNN. The results indicated that NNN catalyzed by the ferric porphyrin (Compound I, Cpd I) at the active center of CYP450 included two steps, hydrogen abstraction and rebound reactions. The free energy barriers of the 2'- and 5'-hydroxylation of NNN are 9.82/8.44 kcal/mol (R/S) and 7.99/9.19 kcal/mol (R/S), respectively, suggesting that the 2'-(S) and 5'-(R) pathways have a slight advantage. The free energy barriers of the decomposition occurred at the 2'-position and 5'-position of NNN are 18.04/18.02 kcal/mol (R/S) and 18.33/19.53 kcal/mol (R/S), respectively. Moreover, we calculated the alkylation reactions occurred at ten DNA base sites induced by the 2'-hydroxylation product of NNN, generating the free energy barriers ranging from 0.86 to 4.72 kcal/mol, which indicated that these reactions occurred easily. The docking study showed that (S)-NNN had better affinity with CYP450s than that of (R)-NNN, which was consistent with the experimental results. Overall, the combined results of the DFT calculations and the docking obtained in this study provide an insight into the understanding of the carcinogenesis of NNN and other TSNAs.

Density functional theory studies on the conversion of hydroxyheme to iron-verdoheme in the presence of dioxygen

Detailed insight into the second step of heme degradation by heme oxygenase, oxophlorin to verdoheme and biliverdin, is presented. Density functional theory methods are reported for the conversion of oxophlorin to verdoheme. Since it is currently unclear whether dioxygen binding to iron oxophlorin is followed by a reduction or not, in this work we have focused on the difference in reactivity between [(Im)(O₂˙)Fe^III(PO˙)] (PO˙ is the oxophlorin dianion radical) and [(Im)(O₂˙)Fe^III(PO)]^- (PO is the oxophlorin trianion). Thus, we have shown that in [(Im)(O₂˙)Fe^III(PO˙)] and [(Im)(O₂˙)Fe^III(PO)]^-, the mechanisms are stepwise with an initial C-O bond activation to form a ring-structure where the oxophlorin is distorted from planarity. This is followed by homolytic dioxygen bond breaking that directly leads to iron-oxo verdoheme products. The [(Im)(O₂˙)Fe^III(PO˙)] mechanism proceeds via two-state-reactivity patterns on the adjacent doublet and quartet spin state surfaces, whereas the [(Im)(O₂˙)Fe^III(PO)]^- route shows single-state-reactivity on a triplet spin state surface. In both, the rate determining step is the C-O bond activation, with substantially lower barriers on the [(Im)(O₂˙)Fe^III(PO˙)] surface of 12.15 kcal mol^-1 in the gas phase compared to 22.55 kcal mol^-1 for the intermediate-spin of [(Im)(O₂˙)Fe^III(PO)]^-. The complete active space self-consistent-field wave functions with second-order multi-reference perturbation theory were also studied. Finally, the effects of the solvent and the medium on the reaction barriers were tested and shown to be considerable.

Copper-Catalyzed Aziridination with Redox-Active Ligands: Molecular Spin Catalysis

Small-molecule catalysts as mimics of biological systems illustrate the chemists' attempts at emulating the tantalizing abilities displayed by nature's metalloenzymes. Among these innate behaviors, spin multistate reactivity is used by biological systems as it offers thermodynamic leverage towards challenging chemical reactivity but this concept is difficult to translate into the realm of synthetic organometallic catalysis. Here, we report a rare example of molecular spin catalysis involving multistate reactivity in a small-molecule biomimetic copper catalyst applied to aziridination. This behavior is supported by spin state flexibility enabled by the redox-active ligand.

Do two oxidants (ferric-peroxo and ferryl-oxo species) act in the biosynthesis of estrogens? A DFT calculation

Density functional theory calculations were performed in order to reveal the mysterious catalytic step of the biosynthesis of estrogens. The results indicated two reactive oxidants, ferric-peroxo and ferryl-oxo (compound I) species, to participate in the conversion of androgens to estrogens. The ferric-peroxo species was determined, according to our derived mechanism, to act in the oxidation of 19-OH androgen to yield the 19,19-gem-diol intermediate and generate the ferryl-oxo (compound I) species. This species was then modeled to effect, in the final step, an abstraction of H from an O–H group of 19,19-gem-diol to give the experimentally observed products. We considered our new mechanistic scenario to reasonably explain the latest experimental observations and to provide deep insight complementing the newly accepted compound I (Cpd I) mechanism.

Density functional theory calculations were performed in order to reveal the mysterious catalytic step of the biosynthesis of estrogens.

VTST/MT studies of the catalytic mechanism of C-H activation by transition metal complexes with [Cu2(μ-O2)], [Fe2(μ-O2)] and Fe(IV)-O cores based on DFT potential energy surfaces

High-valent Cu and Fe species, which are generated from dioxygen activation in metalloenzymes, carry out the functionalization of strong C-H bonds. Understanding the atomic details of the catalytic mechanism has long been one of the main objectives of bioinorganic chemistry. Large H/D kinetic isotope effects (KIEs) were observed in the C-H activation by high-valent non-heme Cu or Fe complexes in enzymes and their synthetic models. The H/D KIE depends significantly on the transition state properties, such as structure, energies, frequencies, and shape of the potential energy surface, when the tunneling effect is large. Therefore, theoretical predictions of kinetic parameters such as rate constants and KIEs can provide a reliable link between atomic-level quantum mechanical mechanisms and experiments. The accurate prediction of the tunneling effect is essential to reproduce the kinetic parameters. The rate constants and HD/KIE have been calculated using the variational transition-state theory including multidimensional tunneling based on DFT potential energy surfaces along the reaction coordinate. Excellent agreement was observed between the predicted and experimental results, which assures the validity of the DFT potential energy surfaces and, therefore, the proposed atomic-level mechanisms. The [Cu₂(μ-O)₂], [Fe₂(μ-O)₂], and Fe(IV)-oxo species were employed for C-H activation, and their role as catalysts was discussed at an atomic level.

Assessment of electronic structure methods for the determination of the ground spin states of Fe(ii), Fe(iii) and Fe(iv) complexes

We studied spin states of Fe²⁺ ion, gaseous FeO, and 14 Fe(ii), Fe(iii) and Fe(iv) complexes using density functional theory.

Computational analysis of site differences in selective aliphatic C–H hydroxylation by nonheme iron–oxo complexes

The CH₃CN-solvent influences selective C–H hydroxylation by nonheme iron complexes due to interactions in the H-bonding formation in H-abstraction.

To rebound or dissociate? This is the mechanistic question in C-H hydroxylation by heme and nonheme metal-oxo complexes

Enzymatic reactions that involve C-H bond activation of alkanes by high-valent iron-oxo species can be explained by the rebound mechanism (RM). Hydroxylation reactions of alkane substrates effected by the reactive compound I (Cpd I) species of cytochrome P450 enzymes are good examples. There was initially little doubt that the rebound paradigm could be carried over in the same form to the arena of synthetic nonheme high-valent iron-oxo or other metal-oxo complexes. However, the active reaction centres of these synthetic systems are not well-caged, in contrast to the active sites of enzymes; therefore, the relative importance of the radical dissociation pathway can become prominent. Indeed, accumulating experimental and theoretical evidence shows that introduction of the non-rebound mechanism (non-RM) is necessary to rationalise the different reactivity patterns observed for synthetic nonheme complexes. In this tutorial review, we discuss several specific examples involving the non-RM while making frequent comparisons to the RM, mainly from the perspective of computational chemistry. We also provide a technical guide to DFT calculations of RM and non-RM and to the interpretations of computational outcomes.

Metabolic fate of 3α,5-cycloandrostanes in the endogenous lactonization pathway of Aspergillus tamarii KITA

A series of 3α,5-cycloandrostane analogues with a range of functionality (6α and 6β alcohols and ketone) at carbon 6 were tested in the endogenous lactonization pathway in Aspergillus tamarii KITA. This metabolic route converts progesterone to testololactone in high yield through a four step enzymatic pathway. To date, no studies have looked at the effect of steroids devoid of polar functionality at carbon 3 and their subsequent metabolic fate by fungi which contain Baeyer-Villiger monooxygenases. Incubation of all of the cycloandrostane analogues resulted in lactonization of ring-D irrespective of C-6 stereochemistry or absence of C-3 functionality. Presence of 6β-hydroxy group and the C-17 ketone was required in order for these analogues to undergo hydroxylation at C-15β position. All metabolites were isolated by column chromatography and were identified by (1)H, (13)C NMR, DEPT analysis and other spectroscopic data.

What factors influence the reactivity of C-H hydroxylation and C=C epoxidation by [Fe(IV)(L(ax))(1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane)(O)](n+)

Density functional theory is used to investigate geometric structures and mechanisms for hydroxylation and epoxidation from propene for four non-heme iron complexes, [Fe(IV)(L(ax))(TMC)(O)](n+), which are the inverted isomers of [Fe(IV)(O)(TMC)(Lax)](n+) (Lax = acetonitrile (AN), monoanionic trifluoroacetate (TF), azide (N3), thiolate (SR)). The Fe(IV)O unit is found to be sterically less hindered in [Fe(IV)(L(ax))(TMC)(O)](n+) than that in [Fe(IV)(O)(TMC)(L(ax))](n+). Becke, three-parameter, Lee-Yang-Parr (B3LYP) calculations show that hydroxylation and epoxidation proceed via a two-state-reactivity on competing triplet and quintet spin surfaces; and the reactions have been invariably mediated by the S = 2 state. The reaction pathways computed reveal that 2-AN is the most reactive in the four [Fe(IV)(L(ax))(TMC)(O)](n+) complexes; along the reaction pathway, the axial ligand moves away from the iron center, and thus, the energy of the LUMO decreases. The anionic axial ligand, which is more electron releasing than neutral AN, shows a strong overlap of iron orbitals. Thus, the anionic ligand does not move away from the iron center. The H-abstraction is affected by the tunneling contribution, the more electron donation power of the axial ligand, the more effect of the tunneling contribution. Adding the tunneling correction, the relative reactivity of the hydroxylation follows the trend: 2-AN > 2-SR ≈ 2-N3 > 2-TF. However, for the epoxidation, the reactivity is in the following order of 2-AN > 2-TF > 2-N3 > 2-SR. Except for 2-AN, 2-X (L(ax) = TF, N3, SR) complexes chemoselectively hydroxylate even in the presence of a C=C double bond.

Nonheme Fe(IV) Oxo Complexes of Two New Pentadentate Ligands and Their Hydrogen-Atom and Oxygen-Atom Transfer Reactions

Two new pentadentate {N5} donor ligands based on the N4Py (N4Py = N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine) framework have been synthesized, viz. [N-(1-methyl-2-benzimidazolyl)methyl-N-(2-pyridyl)methyl-N-(bis-2-pyridyl methyl)amine] (L(1)) and [N-bis(1-methyl-2-benzimidazolyl)methyl-N-(bis-2-pyridylmethyl)amine] (L(2)), where one or two pyridyl arms of N4Py have been replaced by corresponding (N-methyl)benzimidazolyl-containing arms. The complexes [Fe(II)(CH3CN)(L)](2+) (L = L(1) (1); L(2) (2)) were synthesized, and reaction of these ferrous complexes with iodosylbenzene led to the formation of the ferryl complexes [Fe(IV)(O)(L)](2+) (L = L(1) (3); L(2) (4)), which were characterized by UV-vis spectroscopy, high resolution mass spectrometry, and Mössbauer spectroscopy. Complexes 3 and 4 are relatively stable with half-lives at room temperature of 40 h (L = L(1)) and 2.5 h (L = L(2)). The redox potentials of 1 and 2, as well as the visible spectra of 3 and 4, indicate that the ligand field weakens as ligand pyridyl substituents are progressively substituted by (N-methyl)benzimidazolyl moieties. The reactivities of 3 and 4 in hydrogen-atom transfer (HAT) and oxygen-atom transfer (OAT) reactions show that both complexes exhibit enhanced reactivities when compared to the analogous N4Py complex ([Fe(IV)(O)(N4Py)](2+)), and that the normalized HAT rates increase by approximately 1 order of magnitude for each replacement of a pyridyl moiety; i.e., [Fe(IV)(O)(L(2))](2+) exhibits the highest rates. The second-order HAT rate constants can be directly related to the substrate C-H bond dissociation energies. Computational modeling of the HAT reactions indicates that the reaction proceeds via a high spin transition state.

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.

Ab initio dynamics of the cytochrome P450 hydroxylation reaction

The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.

Transition metal complexes and the activation of dioxygen

In order to address how diverse metalloprotein active sites, in particular those containing iron and copper, guide O₂binding and activation processes to perform diverse functions, studies of synthetic models of the active sites have been performed. These studies have led to deep, fundamental chemical insights into how O₂coordinates to mono- and multinuclear Fe and Cu centers and is reduced to superoxo, peroxo, hydroperoxo, and, after O-O bond scission, oxo species relevant to proposed intermediates in catalysis. Recent advances in understanding the various factors that influence the course of O₂activation by Fe and Cu complexes are surveyed, with an emphasis on evaluating the structure, bonding, and reactivity of intermediates involved. The discussion is guided by an overarching mechanistic paradigm, with differences in detail due to the involvement of disparate metal ions, nuclearities, geometries, and supporting ligands providing a rich tapestry of reaction pathways by which O₂is activated at Fe and Cu sites.

Alternative Pathway for the Reaction Catalyzed by DNA Dealkylase AlkB from Ab Initio QM/MM Calculations

AlkB is the title enzyme of a family of DNA dealkylases that catalyze the direct oxidative dealkylation of nucleobases. The conventional mechanism for the dealkylation of N1-methyl adenine (1-meA) catalyzed by AlkB after the formation of FeIV-oxo is comprised by a reorientation of the oxo moiety, hydrogen abstraction, OH rebound from the Fe atom to the methyl adduct, and the dissociation of the resulting methoxide to obtain the repaired adenine base and formaldehyde. An alternative pathway with hydroxide as a ligand bound to the iron atom is proposed and investigated by QM/MM simulations. The results show OH- has a small impact on the barriers for the hydrogen abstraction and OH rebound steps. The effects of the enzyme and the OH- ligand on the hydrogen abstraction by the FeIV-oxo moiety are discussed in detail. The new OH rebound step is coupled with a proton transfer to the OH- ligand and results in a novel zwitterion intermediate. This zwitterion structure can also be characterized as Fe-O-C complex and facilitates the formation of formaldehyde. In contrast, for the pathway with H2O bound to iron, the hydroxyl product of the OH rebound step first needs to unbind from the metal center before transferring a proton to Glu136 or other residue/substrate. The consistency between our theoretical results and experimental findings is discussed. This study provides new insights into the oxidative repair mechanism of DNA repair by nonheme FeII and α-ketoglutarate (α-KG) dependent dioxygenases and a possible explanation for the substrate preference of AlkB.