Correlated Ab Initio and Density Functional Studies on H2 Activation by FeO+

Evidence for the formation of 6PPD-quinone from antioxidant 6PPD by cytochrome P450

N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine (6PPD) as a rubber antioxidant has attracted global concern, since its ozone-oxidation product 6PPD-quinone (6PPDQ) was found to be the primary toxicant responsible for urban runoff mortality syndrome in coho salmon. However, the biotransformation fate and associated toxicological mechanism of 6PPD have not received much study yet. In this work, the in vitro assays showed 6PPD can be transformed into 6PPDQ by cytochromes P450 (CYP450) in human liver microsomes (HLMs) with 0.98 % production rate, and the adducts of 6PPDQ with calf thymus DNA and the N-N coupling product between 6PPD and 6PPDQ were further identified after 6PPD incubation in HLMs. Further evidence for the 6PPDQ formation can be obtained from the in vivo assays that the 6PPDQ-DNA adducts and 6PPD-N-N-6PPDQ dimer were detected in mice by oral gavage with 6PPD, and the latter dimer species was detected as well in 6PPD exposure to zebrafish larvae. Especially, the bioaccumulation property and high reactivity of 6PPDQ result in the continuous formation of the significant DNA adducts and 6PPD-N-N-6PPDQ dimer even in case of low production rate of biotransformation of 6PPD to 6PPDQ, which may provide potentially effective biomarkers for such process. DFT computations revealed the formation mechanism of 6PPDQ is the (N)H-abstraction of 6PPD by CYP450, followed by amino radical rebound at the nearby ortho-carbon, yielding a quinol intermediate due to spin delocalization, that might readily undergo further oxidation by CYP450 into 6PPDQ.

Computational simulations uncover enantioselective metabolism of chiral triazole fungicides by human CYP450 enzymes: A case study of tebuconazole

Tebuconazole (TEB), a prominent chiral triazole fungicide, has been extensively utilized for plant pathogen control globally. Despite experimental evidence of TEB metabolism in mammals, the enantioselectivity in the biotransformation of R- and S-TEB enantiomers by specific CYP450s remains elusive. In this work, integrated in silico simulations were employed to unveil the binding interactions and enantioselective metabolic fate of TEB enantiomers within human CYP1A2, 2B6, 2E1, and 3A4. Molecular dynamics (MD) simulations clearly delineated the binding specificity of R- and S-TEB to the four CYP450s, crucially determining their differences in metabolic activity and enantioselectivity. The primary driving force for robust ligand binding was identified as van der Waals interactions with CYP450s, particularly involving the hydrophobic residues. Mechanistic insights derived from quantum mechanics/molecular mechanics (QM/MM) calculations established C2-methyl hydroxylation as the predominant route of R-/S-TEB metabolism, while C6-hydroxylation and triazol epoxidation were deemed kinetically infeasible pathways. Specifically, the resulting hydroxy-R-TEB metabolite primarily originates from R-TEB biotransformation by 1A2, 2E1 and 3A4, whereas hydroxy-S-TEB is preferentially produced by 2B6. These findings significantly contribute to our comprehension of the binding specificity and enantioselective metabolic fate of chiral TEB by CYP450s, potentially informing further research on human health risk assessment associated with TEB exposure.

Impact of carboxylate ligation on the C-H activation reactivity of a non-heme Fe(IV)O complex: a computational investigation

A comprehensive DFT investigation has been presented to predict how a carboxylate-rich macrocycle would affect the reactivity of a non-heme Fe(IV)O complex towards C-H activation. The popular non-heme iron oxo complex [FeIV(O)(N4Py)]2+, (N4Py = N,N-(bis(2-pyridyl)methyl)N-bis(2-pyridylmethyl)amine) (1), has been selected here as the primary compound. It is transformed to the compound [FeIV(O)(nBu-P2DA)], where nBu-P2DA = N-(1',1'-bis(2-pyridyl)pentyl)iminodiacetate (2) after the replacement of two pyridine donors of N4Py with carboxylate groups. Two other complexes, namely 3 and 4, have been predicted sequentially substituting nitrogen with the carboxylate groups. Ethylbenzene and dihydrotoluene were chosen as substrates. In terms of C-H activation reactivity, an interesting pattern emerges: as the carboxylate group becomes more equatorially enriched, the reactivity increases, following the trend 1 < 2 < 3 < 4. This also aligns with available experimental reports related to complexes 1 and 2. Fe(IV)O complexes exhibit two-state reactivity (triplet and quintet), whereas the quintet state is more favourable due to the stabilization of the transition states through exchange interactions involving a greater number of unpaired electrons. A detailed analysis of the factors influencing reactivity has been performed, including distortion energy (which decreases for the transition state with the addition of carboxylate groups), the triplet-quintet oxidant energy gap (which consistently decreases as carboxylate group enrichment increases), steric factors, and quantum mechanical tunneling. This investigation provides a detailed explanation of the observed outcomes and predicts the higher reactivity of carboxylate-enriched Fe(IV)O complexes. After potential experimental verification, this could lead to the development of new, optimal catalysts for C-H activation.

Effects of Clinical Mutations in the Second Coordination Sphere and Remote Regions on the Catalytic Mechanism of Non-Heme Fe(II)/2-Oxoglutarate-Dependent Aspartyl Hydroxylase AspH

Aspartyl/asparaginyl hydroxylase (AspH) catalyzes the post-translational hydroxylations of vital human proteins, playing an essential role in maintaining their biological functions. Single-point mutations in the Second Coordination Sphere (SCS) and long-range (LR) residues of AspH have been linked to pathological conditions such as the ophthalmologic condition Traboulsi syndrome and chronic kidney disease (CKD). Although the clinical impacts of these mutations are established, there is a critical knowledge gap regarding their specific atomistic effects on the catalytic mechanism of AspH. In this study, we report integrated computational investigations on the potential mechanistic implications of four mutant forms of human AspH with clinical importance: R735W, R735Q, R688Q, and G434V. All the mutant forms exhibited altered binding interactions with the co-substrate 2-oxoglutarate (2OG) and the main substrate in the ferric-superoxo and ferryl complexes, which are critical for catalysis, compared to the wild-type (WT). Importantly, the mutations strongly influence the energetics of the frontier molecular orbitals (FMOs) and, thereby, the activation energies for the hydrogen atom transfer (HAT) step compared to the WT AspH. Insights from our study can contribute to enzyme engineering and the development of selective modulators for WT and mutants of AspH, ultimately aiding in treating cancers, Traboulsi syndrome and, CKD.

Hey ho, where'd the proton go? Final deprotonation of O6 within the S3 state of photosystem II

The deprotonation of O6 within the S3 state marks the final deprotonation event before the formation of oxygen‑oxygen bond interactions and eventual production and release of dioxygen. Gaining a thorough understanding of this event, from the proton acceptors involved, to the exfiltration pathways available, is key in determining the nature of the resulting oxygen species, influencing the mechanism through which the first oxygen‑oxygen bond forms. Computational analysis, using BS-DFT methodologies, showed that proton abstraction by the local Glu189 residue provides consistent evidence against this being a viable mechanistic pathway due to the lack of a stable product structure. In contrast, abstraction via W3 shows an increasingly stable oxo-oxo product state between r[O5O6] = 2.1 Å & 1.9 Å. The resulting oxo-oxo state is stabilised through donation of β electron character from O6 to Mn1 and α electron character from O6 to O5. This donation from the O6 lone pair is shown to be a key factor in stabilising the oxo-oxo state, in addition to showing the initiation of first O5-O6 bond.

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

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

A density functional theory analysis of the C-H activation reactivity of iron(IV)-oxo complexes with an 'O' substituted tetramethylcyclam macrocycle

In this article, we present a meticulous computational study to foresee the effect of an oxygen-rich macrocycle on the reactivity for C-H activation. For this study, a widely studied nonheme Fe(IV)O molecule with a TMC (1,4,8,11-tetramethyl 1,4,8,11-tetraazacyclotetradecane) macrocycle that is equatorially attached to four nitrogen atoms (designated as N4) and acetonitrile as an axial ligand has been taken into account. For the goal of hetero-substitution, step-by-step replacement of the N4 framework with O atoms, i.e., N4, N3O1, N2O2, N1O3, and O4 systems, has been considered, and dihydroanthracene (DHA) has been used as the substrate. In order to neutralise the system and prevent the self-interaction error in DFT, triflate counterions have also been included in the calculations. The study of the energetics of these C-H bond activation reactions and the potential energy surfaces mapped therefore reveal that the initial hydrogen abstraction, which is the rate-determining step, follows the two-state reactivity (TSR) pattern, which means that the originally excited quintet state falls lower in the transition state and the product. The reaction follows the hydrogen atom transfer (HAT) mechanism, as indicated by the spin density studies. The results revealed a fascinating reactivity order, in which the reactivity increases with the enrichment of the oxygen atom in the equatorial position, namely the order follows N4 < N3O1 < N2O2 < N1O3 < O4. The impacts of oxygen substitution on quantum mechanical tunneling and the H/D kinetic isotope effect have also been investigated. When analysing the causes of this reactivity pattern, a number of variables have been identified, including the reactant-like transition structure, spin density distribution, distortion energy, and energies of the electron acceptor orbital, i.e., the energy of the LUMO (σ*z2), which validate the obtained outcome. Our results also show very good agreement with earlier combined experimental and theoretical studies considering TMC and TMCO-type complexes. The DFT predictions reported here will undoubtedly encourage experimental research in this biomimetic field, as they provide an alternative with higher reactivity in which heteroatoms can be substituted for the traditional nitrogen atom.

Insights on cyclophosphamide metabolism and anticancer mechanism of action: A computational study

The oxazaphosphorine cyclophosphamide (CP) is a DNA-alkylating agent commonly used in cancer chemotherapy. This anticancer agent is administered as a prodrug activated by a liver cytochrome P450-catalyzed 4-hydroxylation reaction that yields the active, cytotoxic metabolite. The primary metabolite, 4-hydroxycyclophosphamide, equilibrates with the ring-open aldophosphamide that undergoes β-elimination to yield the therapeutically active DNA cross-linking phosphoramide mustard and the byproduct acrolein. The present paper presents a DFT investigation of the different metabolic phases and an insight into the mechanism by which CP exerts its cytotoxic action. A detailed computational analysis of the energy profiles describing all the involved transformations and the mechanism of DNA alkylation is given with the aim to contribute to an increase of knowledge that, after more than 60 years of unsuccessful attempts, can lead to the design and development of a new generation of oxazaphosphorines.

Dynamical Origin of Rebound versus Dissociation Selectivity during Fe-Oxo-Mediated C-H Functionalization Reactions

The mechanism of catalytic C-H functionalization of alkanes by Fe-oxo complexes is often suggested to involve a hydrogen atom transfer (HAT) step with the formation of a radical-pair intermediate followed by diverging pathways for radical rebound, dissociation, or desaturation. Recently, we showed that in some Fe-oxo reactions, the radical pair is a nonstatistical-type intermediate and dynamic effects control rebound versus dissociation pathway selectivity. However, the effect of the solvent cage on the stability and lifetime of the radical-pair intermediate has never been analyzed. Moreover, because of the extreme complexity of motion that occurs during dynamics trajectories, the underlying physical origin of pathway selectivity has not yet been determined. For the reaction between [(TQA_Cl)FeIVO]+ and cyclohexane, here, we report explicit solvent trajectories and machine learning analysis on transition-state sampled features (e.g., vibrational, velocity, and geometric) that identified the transferring hydrogen atom kinetic energy as the most important factor controlling rebound versus nonrebound dynamics trajectories, which provides an explanation for our previously proposed dynamic matching effect in fast rebound trajectories that bypass the radical-pair intermediate. Manual control of the reaction trajectories confirmed the importance of this feature and provides a mechanism to enhance or diminish selectivity for the rebound pathway. This led to a general catalyst design principle and proof-of-principle catalyst design that showcases how to control rebound versus dissociation reaction pathway selectivity.

Insights into PSII's S3YZ• State: An Electronic and Magnetic Analysis

Using BS-DFT (broken-symmetry density functional theory), the electronic and magnetic properties of the S3YZ• state of photosystem II were investigated and compared to those of the S3 state. While the O5 oxo-O6 hydroxo species presents little difference between the two states, a previously identified [O5O6]3- exhibits reduced stabilization of the O5-O6 shared spin. This species is shown to have some coupling with the YZ• center through Mn1 and O6. Similarly, a peroxo species is found to exhibit significant exchange couplings between the YZ• center and the Mn cluster through Mn1. Mechanistic changes in O-O bond formation in S3YZ• are highlighted by analysis of IBOs (intrinsic bonding orbitals) showing deviation for Mn1 and O6 centered IBOs. This change in coupling interactions throughout the complex as a result of S3YZ• formation presents implications for the determination of the mechanism spanning the end of the S3 and the start of the S4 states, affecting both electron movement and oxygen bond formation.

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.

M11pz: A Nonlocal Meta Functional with Zero Hartree-Fock Exchange and with Broad Accuracy for Chemical Energies and Structures

The accuracy of Kohn-Sham density functional theory depends strongly on the approximation to the exchange-correlation functional. In this work, we present a new exchange-correlation functional called M11pz (M11 plus rung-3.5 terms with zero Hartree-Fock exchange) that is built on the M11plus functional with the goal of using its rung-3.5 terms without a Hartree-Fock exchange term, especially to improve the accuracy for strongly correlated systems. The M11pz functional is optimized with the same local and rung-3.5 ingredients that are used in M11plus but without any percentage of Hartree-Fock exchange. The performance of M11pz is compared with eight local functionals, and M11pz is found to be in top three when the errors or ranks are averaged over eight grouped and partially overlapping databases: AME418/22, atomic and molecular energies; MGBE172, main-group bond energies; TMBE40, transition-metal bond energies; SR309, single-reference systems; MR54, multireference systems; BH192, barrier heights; NC579, noncovalent interaction energies; and MS20, molecular structures. For calculations of band gaps of solids, M11pz is the second best of the nine tested functionals that have zero Hartree-Fock exchange.

New insight into biodegradation mechanism of phenylurea herbicides by cytochrome P450 enzymes: Successive N-demethylation mechanism

Phenylurea herbicides (PUHs) present one of the most important herbicides, which have cause serious effects on ecological environment and humans. Nowadays enzyme strategy shows great advantages in degradation of PUHs. Here density functional theory (DFT), quantitative structure - activity relationship (QSAR) and quantum mechanics/molecular mechanics (QM/MM) approaches are used to investigate the degradation mechanism of PUHs catalyzed by P450 enzymes. Two successive N-demethylation pathways are identified and two hydrogen abstraction (H-abstraction) reaction pathways are identified as the rate-determining step through high-throughput DFT calculations. The Boltzmann-weighted average energy barrier of the second H-abstraction pathway (19.95 kcal/mol) is higher than that of the first H-abstraction pathway (16.80 kcal/mol). Two QSAR models are established to predict the energy barriers of the two H-abstraction pathways based on the quantum chemical descriptors and mordred molecular descriptors. The determination coefficient (R2) values of QSAR models are > 0.9, which reveal that the established QSAR models have great predictive capability. QM/MM calculations indicate that human P450 enzymes are more efficient in degradation of PUHs than crop and weed P450 enzymes. Correlations between energy barriers and key structural/charge parameters are revealed and key parameters that have influence on degradation efficiency of PUHs are identified. This study provides lateral insights into the biodegradation strategy and removal method of PUHs and valuable information for designing or engineering of highly efficient degradation enzymes and genetically modified crops.

Understanding the Multifaceted Mechanism of Compound I Formation in Unspecific Peroxygenases through Multiscale Simulations

Unspecific peroxygenases (UPOs) can selectively oxyfunctionalize unactivated hydrocarbons by using peroxides under mild conditions. They circumvent the oxygen dilemma faced by cytochrome P450s and exhibit greater stability than the latter. As such, they hold great potential for industrial applications. A thorough understanding of their catalysis is needed to improve their catalytic performance. However, it remains elusive how UPOs effectively convert peroxide to Compound I (CpdI), the principal oxidizing intermediate in the catalytic cycle. Previous computational studies of this process primarily focused on heme peroxidases and P450s, which have significant differences in the active site from UPOs. Additionally, the roles of peroxide unbinding in the kinetics of CpdI formation, which is essential for interpreting existing experiments, have been understudied. Moreover, there has been a lack of free energy characterizations with explicit sampling of protein and hydration dynamics, which is critical for understanding the thermodynamics of the proton transport (PT) events involved in CpdI formation. To bridge these gaps, we employed multiscale simulations to comprehensively characterize the CpdI formation in wild-type UPO from Agrocybe aegerita (AaeUPO). Extensive free energy and potential energy calculations were performed in a quantum mechanics/molecular mechanics setting. Our results indicate that substrate-binding dehydrates the active site, impeding the PT from H2O2 to a nearby catalytic base (Glu196). Furthermore, the PT is coupled with considerable hydrogen bond network rearrangements near the active site, facilitating subsequent O-O bond cleavage. Finally, large unbinding free energy barriers kinetically stabilize H2O2 at the active site. These findings reveal a delicate balance among PT, hydration dynamics, hydrogen bond rearrangement, and cosubstrate unbinding, which collectively enable efficient CpdI formation. Our simulation results are consistent with kinetic measurements and offer new insights into the CpdI formation mechanism at atomic-level details, which can potentially aid the design of next-generation biocatalysts for sustainable chemical transformations of feedstocks.

A larger basis set describes atomization energy core-valence correction better than a higher-order coupled-cluster method

The accuracy of coupled-cluster methods for the computation of core-valence correction to atomization energy was assessed. Truncation levels up to CCSDTQP were considered together with (aug-)cc-pwCVnZ (n = D, T, Q, 5) basis sets and three different extrapolation techniques (canonical and flexible Helgaker formula and Riemann zeta function extrapolation). With the exception of CCSD, a more accurate correction can be obtained from a larger basis set using a lower-level coupled-cluster method, and not vice versa. For the CCSD(T) level, it also implies faster computations with modern codes. We also discussed the importance of moving to higher-order or all-electron methods for geometry optimizations. The present study provides the general knowledge needed for the most accurate state-of-the-art computations.

Enantioselective metabolism of novel chiral insecticide Paichongding by human cytochrome P450 3A4: A computational insight

As a novel chiral neonicotinoid insecticide, Paichongding (IPP) has been widely applied in agriculture due to its excellent insecticidal activity. However, the enantioselective metabolism of IPP stereoisomers (5R7R-IPP, 5S7S-IPP, 5R7S-IPP, and 5S7R-IPP) mediated by enzymes in non-target organisms, especially the cytochrome P450s (CYPs), remains unknown. To address this knowledge gap, we developed an integrated computational framework to elucidate the binding interactions and enantioselective metabolism of IPP stereoisomers in human CYP3A4. The results reveal that 5R7R-IPP shows much stronger binding affinity to CYP3A4 than 5S7S-IPP, while enantiomers 5R7S-IPP and 5S7R-IPP have no essential difference in their binding potential, owing to their specific interactions with key CYP3A4 residues. Although enantiomers 5R7R-IPP and 5S7S-IPP feature distinct binding modes resulting from the chiral differences, their transformation activities are slightly different, with C5 and C13 being the primary metabolic sites, respectively. In contrast, CYP3A4 preferably metabolizes 5R7S-IPP over 5S7R-IPP. The metabolism of epimers 5R7R-IPP and 5R7S-IPP share C5-hydroxylation routes due to the conserved 5R-conformaitons, but differ with the transformation routes at C11/C13 and C3 sites. The 7R-chirality of 5S7R-IPP significantly reduces the metabolic potency compared to 5S7S-IPP. CYP3A4-catalyzed hydroxylation and desaturation of IPP stereoisomers generate various chiral metabolites, with C5- and C13-hydroxyIPPs further transforming into depropylated products. Furthermore, the toxicity assessment reveals that IPP, along with the majority of its hydroxylated, desaturated, and depropylated metabolites, can potentially induce adverse effects on human health, specifically hepatotoxicity, respiratory toxicity, and carcinogenicity. This study provides valuable insights into the enantioselective fate of chiral IPP metabolism by CYP3A4, and the identified metabolites can serve as potential biomarkers for monitoring IPP exposure and associated health risk in human body.

Emerging Metabolic Profiles of Sulfonamide Antibiotics by Cytochromes P450: A Computational-Experimental Synergy Study on Emerging Pollutants

Metabolism, especially by CYP450 enzymes, is the main reason for mediating the toxification and detoxification of xenobiotics in humans, while some uncommon metabolic pathways, especially for emerging pollutants, probably causing idiosyncratic toxicity are easily overlooked. The pollution of sulfonamide antibiotics in aqueous system has attracted increasing public attention. Hydroxylation of the central amine group can trigger a series of metabolic processes of sulfonamide antibiotics in humans; however, this work parallelly reported the coupling and fragmenting initiated by amino H-abstraction of sulfamethoxazole (SMX) catalyzed by human CYP450 enzymes. Elucidation of the emerging metabolic profiles was mapped via a multistep synergy between computations and experiments, involving preliminary DFT computations and in vitro and in vivo assays, profiling adverse effects, and rationalizing the fundamental factors via targeted computations. Especially, the confirmed SMX dimer was shown to potentially act as a metabolism disruptor in humans, while spin aromatic delocalization resulting in the low electron donor ability of amino radicals was revealed as the fundamental factor to enable coupling of sulfonamide antibiotics by CYP450 through the nonconventional nonrebound pathway. This work may further strengthen the synergistic use of computations prior to experiments to avoid wasteful experimental screening efforts in environmental chemistry and toxicology.

Direct Dynamics Trajectories Demonstrate Dynamic Matching and Nonstatistical Radical Pair Intermediates during Fe-Oxo-Mediated C-H Functionalization Reactions

The generally proposed mechanism for the reaction between non-heme Fe-oxo complexes and alkane C-H bonds involves a hydrogen atom transfer (HAT) reaction step with a radical pair intermediate that then has competitive radical rebound, dissociation, or desaturation pathways. Here, we report density functional theory-based quasiclassical direct dynamics trajectories that examine post-HAT reaction dynamics. Trajectories revealed that the radical pair intermediate can be a nonstatistical type intermediate without complete internal vibrational redistribution and post-HAT selectivity is generally determined by dynamic effects. Fast rebound trajectories occur through dynamic matching between the rotational motion of the newly formed Fe-OH bond and collision with the alkane radical, and all of this occurs through a nonsynchronous dynamically concerted process that circumvents the radical pair intermediate structure. For radical pair dissociation, trajectories proceeded to the radical pair intermediate for a very brief time, followed by complete dissociation. These trajectories provide a new viewpoint and model to understand the inherent reaction pathway selectivity for non-heme Fe-oxo-mediated C-H functionalization reactions.

How the Conformational Movement of the Substrate Drives the Regioselective C-N Bond Formation in P450 TleB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations

P450 TleB catalyzes the oxidative cyclization of the dipeptide N-methylvalyl-tryptophanol into indolactam V through selective intramolecular C-H bond amination at the indole C4 position. Understanding its catalytic mechanism is instrumental for the engineering or design of P450-catalyzed C-H amination reactions. Using multiscale computational methods, we show that the reaction proceeds through a diradical pathway, involving a hydrogen atom transfer (HAT) from N1-H to Cpd I, a conformational transformation of the substrate radical species, and a second HAT from N13-H to Cpd II. Intriguingly, the conformational transformation is found to be the key to enabling efficient and selective C-N coupling between N13 and C4 in the subsequent diradical coupling reaction. The underlined conformational transformation is triggered by the first HAT, which proceeds with an energy-demanding indole ring flip and is followed by the facile approach of the N13-H group to Cpd II. Detailed analysis shows that the internal electric field (IEF) from the protein environment plays key roles in the transformation process, which not only provides the driving force but also stabilizes the flipped conformation of the indole radical. Our simulations provide a clear picture of how the P450 enzyme can smartly modulate the selective C-N coupling reaction. The present findings are in line with all available experimental data, highlighting the crucial role of substrate dynamics in controlling this highly valuable reaction.

Exploring the Accuracy Limits of PNO-Based Local Coupled-Cluster Calculations for Transition-Metal Complexes

While the domain-based local pair natural orbital coupled-cluster method with singles, doubles, and perturbative triples (DLPNO-CCSD(T)) has proven instrumental for computing energies and properties of large and complex systems accurately, calculations on first-row transition metals with a complex electronic structure remain challenging. In this work, we identify and address the two main error sources that influence the DLPNO-CCSD(T) accuracy in this context, namely, (i) correlation effects from the 3s and 3p semicore orbitals and (ii) dynamic correlation-induced orbital relaxation (DCIOR) effects that are not described by the local MP2 guess. We present a computational strategy that allows us to completely eliminate the DLPNO error associated with semicore correlation effects, while increasing, at the same time, the efficiency of the method. As regards the DCIOR effects, we introduce a diagnostic for estimating the deviation between DLPNO-CCSD(T) and canonical CCSD(T) for systems with significant orbital relaxation.

Endowing homodimeric carbamoyltransferase GdmN with iterative functions through structural characterization and mechanistic studies

Iterative enzymes, which catalyze sequential reactions, have the potential to improve the atom economy and diversity of industrial enzymatic processes. Redesigning one-step enzymes to be iterative biocatalysts could further enhance these processes. Carbamoyltransferases (CTases) catalyze carbamoylation, an important modification for the bioactivity of many secondary metabolites with pharmaceutical applications. To generate an iterative CTase, we determine the X-ray structure of GdmN, a one-step CTase involved in ansamycin biosynthesis. GdmN forms a face-to-face homodimer through unusual C-terminal domains, a previously unknown functional form for CTases. Structural determination of GdmN complexed with multiple intermediates elucidates the carbamoylation process and identifies key binding residues within a spacious substrate-binding pocket. Further structural and computational analyses enable multi-site enzyme engineering, resulting in an iterative CTase with the capacity for successive 7-O and 3-O carbamoylations. Our findings reveal a subclade of the CTase family and exemplify the potential of protein engineering for generating iterative enzymes.

Role of "S" Substitution on C-H Activation Reactivity of Iron(IV)-Oxo Cyclam Complexes: a Computational Investigation

A comprehensive density functional theory (DFT) investigation has been presented in this article to address the role of equatorial sulfur ligation in C-H activation. A non-heme iron-oxo compound with four nitrogen atoms constituting the equatorially connected macrocyclic framework (represented as N₄) [Fe(IV)═O(THC)(CH₃CN)]²⁺(THC = 1,4,8,11-tetrahydro1,4,8,11-tetraazacyclotetradecane) has been considered as the base compound. Other complexes have been anticipated by the sequential replacement of this nitrogen by sulfur, that is, N₄, N₃S₁, N₂S₂, N₁S₃, and S₄. Counterions, as always, have been considered to avoid the self-interaction error in DFT. Generally, the anti-conformers (with respect to equatorial N-H and Fe═O) turned out to be the most stable. It was found that with the enrichment of the equatorial sulfur atom, reactivity increases successively, that is, we get the trend N₄ < N₃S₁ < N₂S₂ < N₁S₃ < S_4. Our investigations have also verified the available experimental results where it has been reported that N₂S₂ is more reactive than N₄ in their mixed conformation. In search of insights into this typical pattern of reactivity, the interplay of several factors has been recognized, such as the distortion energy which decreases for the transition states with the addition of sulfur; the spin density on the oxygen atom which increases implying that the radical character of abstractor increases on sulfur ligation; the energy of the electron acceptor orbital (the lowest unoccupied molecular orbital (σ_z²*)) which decreases continuously with the sulfur substitution; and the triplet-quintet oxidant energy gap which decreases consistently with S enrichment in the equatorial position. The computational predictions reported here, if further validated by experiments, will definitely encourage the synthesis of sulfur-ligated bio-inspired complexes instead of the ones constituting nitrogen exclusively.

Does an Enol Pathway Preclude High Stereoselectivity in Iron-Catalyzed Indole C-H Functionalization via Carbene Insertion?

C-H functionalization of indoles via Fe carbenoids presents an attractive strategy to obtain biologically important structural motifs. However, obtaining good stereoselectivity with Fe has been a significant challenge. It is unclear whether the low selectivity is due to a radical pathway or an ionic mechanism involving metal-free species. We therefore present a density functional theory (DFT) study of indole alkylation with diazoacetates catalyzed by Fe(ClO₄)TMEDA/spirobisoxazoline and myoglobin. We explore three mechanistic pathways: nucleophilic, radical, and oxocarbenium routes. The nucleophilic pathway is the most feasible with the formation of an enol species that tautomerizes to furnish the alkylated indole. While this mechanism is routinely proposed, the stereochemical model has been conspicuously absent until now. We show that the conventionally invoked enol pathway is not responsible for the low enantiomeric excess. The enol intermediate can stay coordinated to the catalyst via different binding sites placing the enol in proximity to the chiral environment and affecting the stereoselective proton transfer. Both the binding strength and the chiral environment are crucial for obtaining high selectivity. Our study provides the much needed insights for the modest-low selectivities of Fe systems and could help in expediting the discovery of an efficient catalytic system. These mechanistic underpinnings could also be applicable to other metal (Rh, Pd, Cu, etc.)-catalyzed X-H insertion reactions.

Mechanistic insights into the chemistry of compound I formation in heme peroxidases: quantum chemical investigations of cytochrome c peroxidase

Peroxidases are heme containing enzymes that catalyze peroxide-dependant oxidation of a variety of substrates through forming key ferryl intermediates, compounds I and II. Cytochrome c peroxidase (Ccp1) has served for decades as a chemical model toward understanding the chemical biology of this heme family of enzymes. It is known to feature a distinctive electronic behaviour for its compound I despite significant structural similarity to other peroxidases. A water-assisted mechanism has been proposed over a dry one for the formation of compound I in similar peroxidases. To better identify the viability of these mechanisms, we employed quantum chemistry calculations for the heme pocket of Ccp1 in three different spin states. We provided comparative energetic and structural results for the six possible pathways that suggest the preference of the dry mechanism energetically and structurally. The doublet state is found to be the most preferable spin state for the mechanism to proceed and for the formation of the Cpd I ferryl-intermediate irrespective of the considered dielectric constant used to represent the solvent environment. The nature of the spin state has negligible effects on the calculated structures but great impact on the energetics. Our analysis was also expanded to explain the major contribution of key residues to the peroxidase activity of Ccp1 through exploring the mechanism at various in silico generated Ccp1 variants. Overall, we provide valuable findings toward solving the current ambiguity of the exact mechanism in Ccp1, which could be applied to peroxidases with similar heme pockets.

Discerning the feasibility of a no-water peroxidase mechanism in the doublet spin state irrespective of the environment surrounding the heme pocket.

Modeling the Ligand Effect on the Structure of CYP 450 Within the Density Functional Theory

An improved understanding of the P450 structure is relevant to the development of biomimetic catalysts and inhibitors for controlled CH-bond activation, an outstanding challenge of synthetic chemistry. Motivated by the experimental findings of an unusually short Fe-S bond of 2.18 Å for the wild-type (WT) OleT P450 decarboxylase relative to a cysteine pocket mutant form (A369P), a computational model that captures the effect of the thiolate axial ligand on the iron-sulfur distance is presented. With the computational efficiency and streamlined analysis in mind, this model combines a cluster representation of the enzyme─40-110 atoms, depending on the heme and ligand truncation level─with a density functional theory (DFT) description of the electronic structure (ES) and is calibrated against the experimental data. The optimized Fe-S distances show a difference of 0.25 Å between the low and high spin states, in agreement with the crystallographic structures of the OleT WT and mutant forms. We speculate that this difference is attributable to the packing of the ligand; the mutant is bulkier due to an alanine-to-proline replacement, meaning that it is excluded from the energetically favored low-spin minimum because of steric constraints. The presence of pure spin-state pairs and the intersection of the low/high spin states for the enzyme model is indicative of the limitations of single-reference ES methods in such systems and emphasizes the significance of using the proper state when modeling the hydrogen atom transfer (HAT) reaction catalyzed by OleT. At the same time, the correct characterization of both the short and long Fe-S bonds within a small DFT-based model of 42 atoms paves the way for quantum dynamics modeling of the HAT step, which initiates the OleT decarboxylation reaction.

Catalytic properties of the ferryl ion in the solid state: a computational review

This review summarises the last findings in the emerging field of heterogeneous catalytic oxidation of light alkanes by ferryl species supported on solid-state systems such as the conversion of methane into methanol by FeO-MOF74.

Computational insight into biotransformation of halophenols by cytochrome P450: Mechanism and reactivity for epoxidation

Halophenols (XPs) have aroused great interests due to their high toxicity and low biodegradability. Previous experimental studies have shown that XPs can be catalytically transformed into epoxides and haloquinones by cytochrome P450 enzymes (CYPs). However, these metabolites have never been detected directly. Moreover, the effects of the reaction site and the type and number of halogen substituents on the biotransformation reactivity of halophenols still remain unknown. In this work, we performed density functional theory (DFT) calculations to simulate the CYP-mediated biotransformation of 36 XPs with mono-, di-, and tri-halogen (F, Cl, and Br) substitutions to unravel the mechanism and relevant kinetics of XPs epoxidation. The whole epoxidation process consists of initial rate-determining O-addition and subsequent ring-closure steps. The simulation results show that the epoxidation in low-spin (LS) state is kinetically preferred over that in high-spin (HS) state, and the formation of epoxide metabolite is strongly exothermic. For all XPs, the epoxidation reactivity follows the order of ortho/para O-addition > meta O-addition. Moreover, the O-addition with higher energy barriers roughly corresponds to chlorophenols and fluorophenols with more halogen atoms. Compared with dichlorophenols, the additional ortho-Cl substitution on trichlorophenols can slightly increase the energy barriers of meta O-addition. By contrast, the additional inclusion of an ortho-Cl to monochlorophenols enhances the meta O-addition reactivity of dichlorophenols. Overall, the present work clarifies the biotransformation routes of XPs to produce epoxides, and identifies the key factors affecting the epoxidation reactivity, which are beneficial in understanding comprehensively the metabolic fate and toxicity of XPs.

Precision Biotransformation of Emerging Pollutants by Human Cytochrome P450 Using Computational-Experimental Synergy: A Case Study of Tris(1,3-dichloro-2-propyl) Phosphate

Precision biotransformation is an envisioned strategy offering detailed insights into biotransformation pathways in real environmental settings using experimentally guided high-accuracy quantum chemistry. Emerging pollutants, whose metabolites are easily overlooked but may cause idiosyncratic toxicity, are important targets of such a strategy. We demonstrate here that complex metabolic reactions of tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) catalyzed by human CYP450 enzymes can be mapped via a three-step synergy strategy: (i) screening the possible metabolites via high-throughout (moderate-accuracy) computations; (ii) analyzing the proposed metabolites in vitro by human liver microsomes and recombinant human CYP450 enzymes; and (iii) rationalizing the experimental data via precise mechanisms using high-level targeted computations. Through the bilateral dialogues from qualitative to semi-quantitative to quantitative levels, we show how TDCIPP metabolism especially by CYP3A4 generates bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) as an O-dealkylation metabolite and bis(1,3-dichloro-2-propyl) 3-chloro-1-hydroxy-2-propyl phosphate (alcohol_β-dehalogen) as a dehalogenation/reduction metabolite via the initial rate-determining H-abstraction from αC- and βC-positions. The relative yield ratio [dehalogenation/reduction]/[O-dealkylation] is derived from the relative barriers of H-abstraction at the βC- and αC-positions by CYP3A4, estimated as 0.002 to 0.23, viz., an in vitro measured ratio of 0.04. Importantly, alcohol_β-dehalogen formation points to a new mechanism involving successive oxidation and reduction functions of CYP450, with its precursor aldehyde_β-dehalogen being a key intermediate detected by trapping assays and rationalized by computations. We conclude that the proposed three-step synergy strategy may meet the increasing challenge of elucidating biotransformation mechanisms of substantial synthesized organic compounds in the future.

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

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

Increased adverse effects during metabolic transformation of short-chain chlorinated paraffins by cytochrome P450: A theoretical insight into 1-chlorodecane

Short-chain chlorinated paraffins (SCCPs), frequently detected in human tissues or organs, can result in threat to human health by disturbing normal metabolism. However, their metabolism mechanisms and fates are largely unclear. Therefore, to better understand the impacts of SCCPs and their metabolites on the human health, the metabolic mechanism and kinetics of SCCPs by cytochrome P450 enzymes (CYPs) were explored using density functional theory employed 1-chlorodecane as a model SCCPs. The results show that 1-chlorodecane could be readily metabolized by CYPs, and the rate constant reaches up 42.3 s^-1 in human body. Dechlorination of 1-chlorodecane is unlikely to occur and hydroxylation is dominated via H-abstraction pathways, especially from the intermediate C atom of 1-chlorodecane. The toxicity assessments suggest that the two metabolites, 10-chloro-decan-5-ol and 1-chlorodecanol could exhibit higher bioaccumulation, carcinogenicity and more serious damage on cardiovascular system after the metabolism of 1-chlorodecane. To our knowledge, this is the first study from the viewpoint of theoretical analysis to explore the metabolism of typical SCCPs in human body. It may provide deep insight into the metabolic transformation mechanism of SCCPs and cause the concerns about the adverse effects of their metabolites in human body.

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe²⁺, Fe³⁺, and Co³⁺, including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. We find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method-in terms of computational cost and in terms of accuracy on monatomic systems-and therefore we chose this method for molecular benchmarks. Then Kohn-Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF₂, FeCl₂, and CoF₂. We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

Electronic Structure and Magnetic Properties of a Titanium(II) Coordination Complex

Stable coordination complexes of Ti^II (3d²) are relatively uncommon, but are of interest as synthons for low oxidation state titanium complexes for application as potential catalysts and reagents for organic synthesis. Specifically, high-spin Ti^II ions supported by redox-inactive ligands are still quite rare due to the reducing power of this soft ion. Among such Ti^II complexes is trans-[TiCl₂(tmeda)₂], where tmeda = N,N,N',N'-tetramethylethane-1,2-diamine. This complex was first reported by Gambarotta and co-workers almost 30 years ago, but it was not spectroscopically characterized and theoretical investigation by quantum chemical theory (QCT) was not feasible at that time. As part of our interest in low oxidation state early transition metal complexes, we have revisited this complex and report a modified synthesis and a low temperature (100 K) crystal structure that differs slightly from that originally reported at ambient temperature. We have used magnetometry, high-frequency and -field EPR (HFEPR), and variable-temperature variable-field magnetic circular dichroism (VTVH-MCD) spectroscopies to characterize trans-[TiCl₂(tmeda)₂]. These techniques yield the following S = 1 spin Hamiltonian parameters for the complex: D = -5.23(1) cm^-1, E = -0.88(1) cm^-1, (E/D = 0.17), g = [1.86(1), 1.94(2), 1.77(1)]. This information, in combination with electronic transitions from MCD, was used as input for both classical ligand-field theory (LFT) and detailed QCT studies, the latter including both density functional theory (DFT) and ab initio methods. These computational methods are seldom applied to paramagnetic early transition metal complexes, particularly those with S > 1/2. Our studies provide a complete picture of the electronic structure of this complex that can be put into context with the few other high-spin and mononuclear Ti^II species characterized to date.

Spin-Forbidden Reactivity of Transition Metal Oxo Species: Exploring the Potential Energy Surfaces

Spin-forbidden reactions are frequently encountered when transition metal oxo species are involved, particularly in oxygen transfer reactivity. The computational study of such reactions is challenging, because reactants and products are located on different spin potential energy surfaces (PESs). One possible approach to describe these reactions is the so-called minimum energy crossing point (MECP) between the diabatic reactants and products PESs. Alternatively, inclusion of spin-orbit coupling (SOC) effects allows to locate a saddle point on a single adiabatic PES (TS SOC). The TS SOC approach is rarely applied because of its high computational cost. Recently evidence for a TS SOC impact on significantly lowering the activation barrier in dioxygen addition to a carbene-gold(I)-hydride complex reaction (Chem. Sci. 2016, 7, 7034-7039) or even on predicting a qualitatively different reaction mechanism in mercury methylation by cobalt corrinoid (Angew. Chem. Int. Ed. 2016, 55, 11503-11506) has been put forward. Using MECP and TS SOC approaches a systematic analysis is provided here of three prototypical transition metal oxo spin-forbidden processes to investigate their implications on reactivity. Cycloaddition of ethylene to chromyl chloride (CrO₂ Cl₂ +C₂ H₄ ), iron oxide cation insertion into the hydrogen molecule (FeO⁺ +H₂ ) and H-abstraction from toluene by a Mn^V -oxo-porphyrin cation (MnOP(H₂ O)⁺ +C₆ H₅ CH₃ ) are case studies. For all these processes the MECP and TS SOC results are compared, which show that the spin-forbidden reactivity of transition metal oxo species can be safely described by a MECP approach, at least for the first-row transition metals investigated here, where the spin-orbit coupling is relatively weak. However, for the Mn-oxo reactivity, the MECP and TS SOC have been found to be crucial for a correct description of the reaction mechanism. In particular, the TS SOC approach allows to straightforwardly explore detailed features of the adiabatic potential energy surface which in principle could affect the overall reaction rate in cases where the involved diabatic PESs are tricky.

Computational Biotransformation Profile of Emerging Phenolic Pollutants by Cytochromes P450: Phenol-Coupling Mechanism

Phenols are ubiquitous environmental pollutants, whose biotransformation involving phenol coupling catalyzed by cytochromes P450 may produce more lipophilic and toxic metabolites. Density functional theory (DFT) computations were performed to explore the debated phenol-coupling mechanisms, taking triclosan as a model substrate. We find that a diradical pathway facilitated by compound I and protonated compound II of P450 is favored vs alternative radical addition or electron-transfer mechanisms. The identified diradical coupling resembles a "two-state reactivity" from compound I characterized by significantly high rebound barriers of the phenoxy radicals, which can be formulated into three equations for calculating the ratio [coupling]/[hydroxylation]. A higher barrier for rebound than for H-abstraction in high-spin triclosan can facilitate the phenoxy radical dissociation and thus enable phenol coupling, while H-abstraction/radical rebound causing phenol hydroxylation via minor rebound barriers mostly occurs via the low-spin state. Therefore, oxidation of triclosan by P450 fits the first equation with a ratio [coupling]/[hydroxylation] of 1:4, consistent with experimental data indicating different extents of triclosan coupling (6-40%). The high rebound barrier of phenoxy radicals, as a key for the mechanistic identification of phenol coupling vs hydroxylation, originates from their weak electron donor ability due to spin aromatic delocalization. We envision that the revealed mechanism can be extended to the cross-coupling reactions between different phenolic pollutants, and the coupling reactions of several other aromatic pollutants, to infer unknown metabolites.

Understanding and Predicting Post H-Atom Abstraction Selectivity through Reactive Mode Composition Factor Analysis

The selective functionalization of C-H bonds is one of the Grails of synthetic chemistry. In this work, we demonstrate that the selectivity toward fast hydroxylation or radical diffusion (known as the OH-rebound and dissociation mechanisms) following H-atom abstraction (HAA) from a substrate C-H bond by high-valent iron-oxo oxidants is already encoded in the HAA step when the post-HAA barriers are much lower than the preceding one. By applying the reactive mode composition factor (RMCF) analysis, which quantifies the kinetic energy distribution (KED) at the reactive mode (RM) of transition states, we show that reactions following the OH-rebound coordinate concentrate the RM kinetic energy on the motion of the reacting oxygen atom and the nascent substrate radical, whereas reactions following the dissociation channel localize most of their kinetic energy in H-atom motion. These motion signatures serve to predict the post-HAA selectivity, and since KED is affected by the free energy of reaction and asynchronicity (factor η) of HAA, we show that bimolecular HAA reactions in solution that are electron transfer-driven and highly exergonic have the lowest fraction of KED on the transferred H-atom and the highest chance to follow rebound hydroxylation. Finally, the RMCF analysis predicts that the H/D primary kinetic isotope effect can serve as a probe for these mechanisms, as confirmed in virtually all reported examples in the literature.

Molecular Basis for Metabolic Regioselectivity and Mechanism of Cytochrome P450s toward Carcinogenic 4-(Methylnitrosamino)-(3-pyridyl)-1-butanone

As an abundantly present tobacco component, carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) has also been detected in atmospheric particulate matter, suggesting the ineluctable exposure risk of this contaminant. NNK metabolic activation by cytochrome P450 enzymes (CYPs) is a prerequisite to exerting its genotoxicity, but the metabolic regioselectivity and mechanism are still unknown. Here the binding feature and regioselectivity of CYPs 1A1, 1A2, 2A6, 2A13, 2B6, and 3A4 toward NNK are unraveled through molecular docking and molecular dynamics (MD) simulations. Binding mode analyses reveal that 1A2 and 2B6 have definite preferences for NNK α-methyl hydroxylation, while the other four CYPs preferentially catalyze α-methylene hydroxylation. The binding affinities between NNK and CYPs evaluated by the binding free energies follow the order 2A13 > 2B6 > 1A2 > 2A6 > 1A1 > 3A4. Density functional theory (DFT) calculations are further performed to characterize the mechanism of NNK biotransformation. Results show that the α-hydroxyNNK generated from α-hydroxylation may undergo nonenzymatic decomposition to form genotoxic diazohydroxide and aldehyde, and further oxidation by P450 to yield nitrosamide, which mainly contributes to NNK toxification capacity. Meanwhile the pyridine N-oxidation and denitrosation of C_α-radical intermediate play an important role in detoxifying NNK. Overall, the present study provides the molecular basis for CYP-catalyzed regioselectivity and mechanism of NNK biotransformation, which can enable the identification of metabolites for assessing the health risk of individual NNK exposure.

Chlorination versus hydroxylation selectivity mediated by the non-heme iron halogenase WelO5

O₂ activation in WelO5 (an α-KG dependent halogenase) leads to a Fe(IV)O species with an equatorial conformation. After hydrogen abstraction, the hydroxyl ligand is far from the substrate radical which leads to the chlorination selectivity

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.

QM/MM MD simulations reveal an asynchronous PCET mechanism for nitrite reduction by copper nitrite reductase

The nitrite reduction in copper nitrite reductase is found to proceed through an asynchronous proton-coupled electron transfer (PCET) mechanism, with electron transfer from T1-Cu to T2-Cu preceding the proton transfer from Asp98 to nitrite.

Fenton-Derived OH Radicals Enable the MPnS Enzyme to Convert 2-Hydroxyethylphosphonate to Methylphosphonate: Insights from Ab Initio QM/MM MD Simulations

The mechanism for dioxygen activation represents one of the core issues in metalloenzymes. In most cases, the activation of the O₂ molecule requires additional electrons from an external reducant. However, nonheme hydroxyethylphosphonate dioxygenase (HEPD) and methylphosphonate synthase (MPnS) are exceptional C-H oxygenases. Both enzymes do not utilize reductants, rather they employ directly iron(III)-superoxide species to initiate H-abstraction reactions and lead thereby to catalysis of the C-C cleavage in 2-hydroxyethylphosphonate (2-HEP). Using the recently characterized MPnS structure and QM(B3LYP)/MM-based metadynamics simulations, we deciphered the chemical mechanism for MPnS. Our simulations demonstrate O₂ activation in MPnS is mediated by an adjacent Lysine residue (Lys28) in the active site, leading to an unusual H ₂ O ₂ intermediate in the reductant-independent nonheme MPnS enzyme. Furthermore, the so-generated H ₂ O ₂ intermediate is subsequently employed in a Fenton-type reaction, leading to a locked •OH radical that spontaneously attaches to the substrate carbonyl group. Meanwhile, the proton from the Fe(III)-OH is shuttled back to the deprotonated Lys28, affording the Fe(IV)-oxo species that is identified by experiment in HEPD. Thus, our calculations demonstrate an unusual proton-shuttle mechanism for O ₂ activation in metalloenzymes.