A predictive pattern of computed barriers for C-h hydroxylation by compound I of cytochrome p450

Mechanism of the Oxidative Ring-Closure Reaction during Gliotoxin Biosynthesis by Cytochrome P450 GliF

During gliotoxin biosynthesis in fungi, the cytochrome P450 GliF enzyme catalyzes an unusual C-N ring-closure step while also an aromatic ring is hydroxylated in the same reaction cycle, which may have relevance to drug synthesis reactions in biotechnology. However, as the details of the reaction mechanism are still controversial, no applications have been developed yet. To resolve the mechanism of gliotoxin biosynthesis and gain insight into the steps leading to ring-closure, we ran a combination of molecular dynamics and density functional theory calculations on the structure and reactivity of P450 GliF and tested a range of possible reaction mechanisms, pathways and models. The calculations show that, rather than hydrogen atom transfer from the substrate to Compound I, an initial proton transfer transition state is followed by a fast electron transfer en route to the radical intermediate, and hence a non-synchronous hydrogen atom abstraction takes place. The radical intermediate then reacts by OH rebound to the aromatic ring to form a biradical in the substrate that, through ring-closure between the radical centers, gives gliotoxin products. Interestingly, the structure and energetics of the reaction mechanisms appear little affected by the addition of polar groups to the model and hence we predict that the reaction can be catalyzed by other P450 isozymes that also bind the same substrate. Alternative pathways, such as a pathway starting with an electrophilic attack on the arene to form an epoxide, are high in energy and are ruled out.

Unraveling the Valence State and Reactivity of Copper Centers in Membrane-Bound Particulate Methane Monooxygenase

Particulate methane monooxygenase (pMMO) plays a critical role in catalyzing the conversion of methane to methanol, constituting the initial step in the C1 metabolic pathway within methanotrophic bacteria. However, the membrane-bound pMMO's structure and catalytic mechanism, notably the copper's valence state and genuine active site for methane oxidation, have remained elusive. Based on the recently characterized structure of membrane-bound pMMO, extensive computational studies were conducted to address these long-standing issues. A comprehensive analysis comparing the quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulated structures with cryo-EM data indicates that both the CuC and CuD sites tend to stay in the Cu(I) valence state within the membrane environment. Additionally, the concurrent presence of Cu(I) at both CuC and CuD sites leads to the significant reduction of the ligand-binding cavity situated between them, making it less likely to accommodate a reductant molecule such as durohydroquinone (DQH2). Subsequent QM/MM calculations reveal that the CuD(I) site is more reactive than the CuC(I) site in oxygen activation, en route to H2O2 formation and the generation of Cu(II)-O•- species. Finally, our simulations demonstrate that the natural reductant ubiquinol (CoQH2) assumes a productive binding conformation at the CuD(I) site but not at the CuC(I) site. This provides evidence that the true active site of membrane-bound pMMOs may be CuD rather than CuC. These findings clarify pMMO's catalytic mechanism and emphasize the membrane environment's pivotal role in modulating the coordination structure and the activity of copper centers within pMMO.

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.

Computational Mechanistic Investigations of Biocatalytic Nitrenoid C-H Functionalizations via Engineered Heme Proteins

Engineered heme proteins were developed to possess numerous excellent biocatalytic nitrenoid C-H functionalizations. Computational approaches such as density functional theory (DFT), hybrid quantum mechanics/molecular mechanics (QM/MM), and molecular dynamics (MD) calculations were employed to help understand some important mechanistic aspects of these heme nitrene transfer reactions. This review summarizes advances of computational reaction pathway results of these biocatalytic intramolecular and intermolecular C-H aminations/amidations, focusing on mechanistic origins of reactivity, regioselectivity, enantioselectivity, diastereoselectivity as well as effects of substrate substituent, axial ligand, metal center, and protein environment. Some important common and distinctive mechanistic features of these reactions were also described with brief outlook of future development.

Melatonin Activation by Cytochrome P450 Isozymes: How Does CYP1A2 Compare to CYP1A1?

Cytochrome P450 enzymes are versatile enzymes found in most biosystems that catalyze mono-oxygenation reactions as a means of biosynthesis and biodegradation steps. In the liver, they metabolize xenobiotics, but there are a range of isozymes with differences in three-dimensional structure and protein chain. Consequently, the various P450 isozymes react with substrates differently and give varying product distributions. To understand how melatonin is activated by the P450s in the liver, we did a thorough molecular dynamics and quantum mechanics study on cytochrome P450 1A2 activation of melatonin forming 6-hydroxymelatonin and N-acetylserotonin products through aromatic hydroxylation and O-demethylation pathways, respectively. We started from crystal structure coordinates and docked substrate into the model, and obtained ten strong binding conformations with the substrate in the active site. Subsequently, for each of the ten substrate orientations, long (up to 1 μs) molecular dynamics simulations were run. We then analyzed the orientations of the substrate with respect to the heme for all snapshots. Interestingly, the shortest distance does not correspond to the group that is expected to be activated. However, the substrate positioning gives insight into the protein residues it interacts with. Thereafter, quantum chemical cluster models were created and the substrate hydroxylation pathways calculated with density functional theory. These relative barrier heights confirm the experimental product distributions and highlight why certain products are obtained. We make a detailed comparison with previous results on CYP1A1 and identify their reactivity differences with melatonin.

Spotlight on CYP4B1

The mammalian cytochrome P450 monooxygenase CYP4B1 can bioactivate a wide range of xenobiotics, such as its defining/hallmark substrate 4-ipomeanol leading to tissue-specific toxicities. Similar to other members of the CYP4 family, CYP4B1 has the ability to hydroxylate fatty acids and fatty alcohols. Structural insights into the enigmatic role of CYP4B1 with functions in both, xenobiotic and endobiotic metabolism, as well as its unusual heme-binding characteristics are now possible by the recently solved crystal structures of native rabbit CYP4B1 and the p.E310A variant. Importantly, CYP4B1 does not play a major role in hepatic P450-catalyzed phase I drug metabolism due to its predominant extra-hepatic expression, mainly in the lung. In addition, no catalytic activity of human CYP4B1 has been observed owing to a unique substitution of an evolutionary strongly conserved proline 427 to serine. Nevertheless, association of CYP4B1 expression patterns with various cancers and potential roles in cancer development have been reported for the human enzyme. This review will summarize the current status of CYP4B1 research with a spotlight on its roles in the metabolism of endogenous and exogenous compounds, structural properties, and cancer association, as well as its potential application in suicide gene approaches for targeted cancer therapy.

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

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

Comparative oxidative ability of mononuclear and dinuclear high-valent iron-oxo species towards the activation of methane: does the axial/bridge atom modulate the reactivity?

Over the years, mononuclear Fe^IVO species have been extensively studied, but the presence of dinuclear Fe^IVO species in soluble methane monooxygenase (sMMO) has inspired the development of biomimic models that could activate inert substrates such as methane. There are some successful attempts; particularly the [(Por)(m-CBA) Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- species has been reported to activate methane and yield decent catalytic turnover numbers and therefore regarded as the closest to the sMMO enzyme functional model, as no mononuclear Fe^IVO analogues could achieve this feat. In this work, we have studied a series of mono and dinuclear models using DFT and ab initio DLPNO-CCSD(T) calculations to probe the importance of nuclearity in enhancing the reactivity. We have probed the catalytic activities of four complexes: [(HO)Fe^IV(O)(Por)]^- (1), [(HO)Fe^IV(O)(Por˙⁺)] (2), μ-oxo dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-O)Fe^IV(O) (Por˙⁺)]^- (3) and N-bridged dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- (4) towards the activation of methane. Additionally, calculations were performed on the mononuclear models [(X)Fe^IV(O)(Por˙⁺)]ⁿ {X = N 4a (n = -2), NH 4b (n = -1) and NH₂4c (n = 0)} to understand the role of nuclearity in the reactivity. DFT calculations performed on species 1-4 suggest an interesting variation among them, with species 1-3 possessing an intermediate spin (S = 1) as a ground state and species 4 possessing a high-spin (S = 2) as a ground state. Furthermore, the two Fe^IV centres in species 3 and 4 are antiferromagnetically coupled, yielding a singlet state with a distinct difference in their electronic structure. On the other hand, species 2 exhibits a ferromagnetic coupling between the Fe^IV and the Por˙⁺ moiety. Our calculations suggest that the higher barriers for the C-H bond activation of methane and the rebound step for species 1 and 3 are very high in energy, rendering them unreactive towards methane, while species 2 and 4 have lower barriers, suggesting their reactivity towards methane. Studies on the system reveal that model 4a has multiple FeN bonds facilitating greater reactivity, whereas the other two models have longer Fe-N bonds and less radical character with steeper barriers. Strong electronic cooperativity is found to be facilitated by the bridging nitride atom, and this cooperativity is suppressed by substituents such as oxygen, rendering them inactive. Thus, our study unravels that apart from enhancing the nuclearity, bridging atoms that facilitate strong cooperation between the metals are required to activate very inert substrates such as methane, and our results are broadly in agreement with earlier experimental findings.

Exploring the Factors which Result in Cytochrome P450 Catalyzed Desaturation Versus Hydroxylation

The cytochrome P450 family of monooxygenase enzymes have essential biological roles involving the selective oxidation of carbon-hydrogen bonds. They can also catalyze other important metabolic reactions including desaturation to form alkenes. Currently the factors that control the partition between P450 hydroxylation and desaturation pathways are poorly defined. The CYP199A4 enzyme from the bacterium Rhodopseudomonas palustris HaA2 catalyzes the oxidation of 4-ethyl- and 4-isopropyl- benzoic acids with hydroxylation and desaturation occurring in significant quantities. Here we demonstrate that 4-cyclopropylbenzoic acid is regioselectively hydroxylated by CYP199A4 at the benzylic carbon. In contrast, the oxidation of 4-n-propylbenzoic acid by CYP199A4 results in three major metabolites: an alkene from desaturation and two hydroxylation products at the benzylic (Cα) and Cβ carbons in similar quantities. Extending the length of the alkyl substituent resulted in 4-n-butylbenzoic acid being oxidized at the benzylic position (45%) and desaturated (55%). In contrast, 4-isobutylbenzoic generated very little alkene (5%) but was hydroxylated at the benzylic position (54%) and at the tertiary Cβ position (41%). The oxidation of 4-n-propylbenzoic acid by the F298 V mutant of CYP199A4 occurred with no hydroxylation at Cβ and a significant increase in metabolites arising from desaturation (73%). The X-ray crystal structures of CYP199A4 with each substrate revealed that they bind in the active site with the alkyl substituent positioned over the heme. However, the longer alkylbenzoic acids were bound in a different conformation as was 4-n-propylbenzoic acid in the F298 V mutant. Overall, the changes in metabolite distribution could be ascribed to bond strength differences and the position of the alkyl group relative to the heme.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

Rapid and Accurate Estimation of Activation Free Energy in Hydrogen Atom Transfer-Based C-H Activation Reactions: From Empirical Model to Artificial Neural Networks

A well-performing machine learning (ML) model is obtained by using proper descriptors and artificial neural network (ANN) algorithms, which can quickly and accurately predict activation free energy in hydrogen atom transfer (HAT)-based sp³ C-H activation. Density functional theory calculations (UωB97X-D) are used to establish the reaction system data sets of methoxyl (CH₃O·), trifluoroethoxyl (CF₃CH₂O·), tert-butoxyl (tBuO·), and cumyloxyl (CumO·) radicals. The simplified Roberts' equation proposed in our recent study works here [R ² = 0.84, mean absolute error (MAE) = 0.85 kcal/mol]. Its performance is comparable with univariate Mulliken-type electronegativity (χ) with the ANN model. The ANN model with bond dissociation free energy, χ, α-unsaturation, and Nolan buried volume (%V _buried) successively improves R ² and MAE to 0.93 and 0.54 kcal/mol, respectively. It reproduces the test sets of trichloroethoxyl (CCl₃CH₂O·) with R ² = 0.87 and MAE = 0.89 kcal/mol and accurately predicts the relative experimental barrier of the HAT reactions with CumO· and the site selectivity of CH₃O·.

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.

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

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

Using Physical Organic Chemistry Knowledge to Predict Unusual Metabolites of Synthetic Phenolic Antioxidants by Cytochrome P450

Biotransformation, especially by human CYP450 enzymes, plays a crucial role in regulating the toxicity of organic compounds in organisms, but is poorly understood for most emerging pollutants, as their numerous "unusual" biotransformation reactions cannot retrieve examples from the textbooks. Therefore, in order to predict the unknown metabolites with altering toxicological profiles, there is a realistic need to develop efficient methods to reveal the "unusual" metabolic mechanism of emerging pollutants. Combining experimental work with computational predictions has been widely accepted as an effective approach in studying complex metabolic reactions; however, the full quantum chemical computations may not be easily accessible for most environmentalists. Alternatively, this work practiced using the concepts from physical organic chemistry for studying the interrelationships between structure and reactivity of organic molecules, to reveal the "unusual" metabolic mechanism of synthetic phenolic antioxidants catalyzed by CYP450, for which the simple pencil-and-paper and property-computation methods based on physical organic chemistry were performed. The phenol-coupling product of butylated hydroxyanisole (BHA) (based on spin aromatic delocalization) and ipso-addition quinol metabolite of butylated hydroxytoluene (BHT) (based on hyperconjugative effect) were predicted as two "unusual" metabolites, which were further confirmed by our in vitro analysis. We hope this easily handled approach will promote environmentalists to attach importance to physical organic chemistry, with an eye to being able to use the knowledge gained to efficiently predict the fates of substantial unknown synthesized organic compounds in the future.

Hydrogen Abstraction by Alkoxyl Radicals: Computational Studies of Thermodynamic and Polarity Effects on Reactivities and Selectivities

Density functional theory calculations (ωB97X-D) are reported for the reactions of methoxy, tert-butoxy, trichloroethoxy, and trifluoroethoxy radicals with a series of 26 C-H bonds in different environments characteristic of a variety of hydrocarbons and substituted derivatives. The variations in activation barriers are analyzed with modified Evans-Polanyi treatments to account for polarity and unsaturation effects. The treatments by Roberts and Steel and by Mayer have inspired the development of a simple treatment involving the thermodynamics of reactions, the difference between the reactant radical and product radical electronegativities, and the absence or presence of α-unsaturation. The three-parameter equation (ΔH^⧧ = 0.52ΔH_rxn(1 - d) - 0.35Δχ_AB² + 10.0, where d = 0.44 when there is α-unsaturation to the reacting C-H bond), correlates well with quantum mechanically computed barriers and shows the quantitative importance of the thermodynamics of reactions (dictated by the reactant and the product bond dissociation energies) and polar effects.

Determining the inherent selectivity for carbon radical hydroxylation versus halogenation with high-spin oxoiron(iv)-halide complexes: a concerted rebound step

A synthetic iron model can process both halogenation and hydroxylation with vague selectivity, which is different from halogenase even though these structures are used for the simulation of halogenase. The key factor of the synthetic oxo-iron model mediated hydroxylation or the halogenation is still under debate. Herein density functional theory calculation is used to investigate the hydroxylation versus halogenation of propylene by the complex [Fe^IV(O)(TQA)(X)]⁺ (X = F, Cl, Br). Our results suggest that a concerted rebound mechanism (between the -X and the hydroxyl ligands after the hydrogen abstraction) leads to the formation of two different kinds of products.

DFT calculation for the hydroxylation versus halogenation of propylene by [Fe^IV(O)(TQA)X]⁺ (X = F, Cl and Br) reveals that after hydrogen abstraction, halogen and oxygen rebound reactions are a synergistic process.

Perspectives on Ligand Properties of N-Heterocyclic Carbenes in Iron Porphyrin Complexes

There has been considerable research interest in the ligand nature of N-heterocyclic carbenes (NHCs). In this work, two six-coordinate NHC iron porphyrin complexes [Fe^II(TTP)(1,3-Me₂Imd)₂] (TTP = tetratolylporphyrin, 1,3-Me₂Imd = 1,3-dimethylimidazol-2-ylidene) and [Fe^III(TDCPP)(1,3-Me₂Imd)₂]ClO₄ (TDCPP = 5,10,15,20-tetrakis(2,6-dichlorophenyl)porphyrin) are reported. Single-crystal X-ray characterizations demonstrate that both complexes have strongly ruffled conformations and relatively perpendicular ligand orientations which are forced by the sterically bulky 1,3-Me₂Imd NHC ligands. Multitemperature (4.2-300 K) and high magnetic field (0-9 T) Mössbauer and low-temperature (4.0 K) EPR spectroscopies definitely confirmed the low-spin states of [Fe^II(TTP)(1,3-Me₂Imd)₂] (S = 0) and [Fe^III(TDCPP)(1,3-Me₂Imd)₂]ClO₄ (S = 1/2). The similarity of 1,3-Me₂Imd and imidazole, as well as the well-established correlations between the ligand nature and spectroscopic characteristics of [Fe^II,III(Porph)(L)₂]^0,+ (Porph: porphyrin; L: planar base ligand) species, allowed direct comparisons between the pair of ligands which revealed for the first time that NHC has a stronger π-acceptor ability than imidazoles, in addition to its very strong σ-donation.

Water biocatalytic effect attenuates cytochrome P450-mediated carcinogenicity of diethylnitrosamine: A computational insight

The mechanism-based mutagenicity and carcinogenicity of diethylnitrosamine (DEN) are believed to act through interactions with cytochrome P450 (P450) enzymes. DFT calculations to explore the conceivable mechanisms underlying the reaction of P450 with DEN with and without water as a biocatalyst were performed. The results shed light on the biocatalytic role of water in lowering the H-abstraction energy barriers because of the electrostatic effect driven by hydrogen bonding. Our DFT analysis revealed how metabolites are formed in the dealkylation (toxification) and denitrosation (detoxification) pathways. Also, our findings uncovered the active position of DEN vulnerable to P450 interactions. Two factors control the toxification and detoxification rates: the stability of denitrosation products and the HS rebound barrier of the α-pathway. Thus, water biocatalytic attenuation of DEN carcinogenicity was attained by stabilizing denitrosation products and slowing the α-HS rebound process. Docking and MD simulations were performed to assess the binding modes of DEN to P450's active site and to inspect the denitrosation and dealkylation processes, respectively.

Bimodal Evans-Polanyi Relationships in Hydrogen Atom Transfer from C(sp3)-H Bonds to the Cumyloxyl Radical. A Combined Time-Resolved Kinetic and Computational Study

The applicability of the Evans-Polanyi (EP) relationship to HAT reactions from C(sp3)-H bonds to the cumyloxyl radical (CumO•) has been investigated. A consistent set of rate constants, kH, for HAT from the C-H bonds of 56 substrates to CumO•, spanning a range of more than 4 orders of magnitude, has been measured under identical experimental conditions. A corresponding set of consistent gas-phase C-H bond dissociation enthalpies (BDEs) spanning 27 kcal mol-1 has been calculated using the (RO)CBS-QB3 method. The log kH' vs C-H BDE plot shows two distinct EP relationships, one for substrates bearing benzylic and allylic C-H bonds (unsaturated group) and the other one, with a steeper slope, for saturated hydrocarbons, alcohols, ethers, diols, amines, and carbamates (saturated group), in line with the bimodal behavior observed previously in theoretical studies of reactions promoted by other HAT reagents. The parallel use of BDFEs instead of BDEs allows the transformation of this correlation into a linear free energy relationship, analyzed within the framework of the Marcus theory. The ΔG⧧HAT vs ΔG°HAT plot shows again distinct behaviors for the two groups. A good fit to the Marcus equation is observed only for the saturated group, with λ = 58 kcal mol-1, indicating that with the unsaturated group λ must increase with increasing driving force. Taken together these results provide a qualitative connection between Bernasconi's principle of nonperfect synchronization and Marcus theory and suggest that the observed bimodal behavior is a general feature in the reactions of oxygen-based HAT reagents with C(sp3)-H donors.

Product Distributions of Cytochrome P450 OleTJE with Phenyl-Substituted Fatty Acids: A Computational Study

There are two types of cytochrome P450 enzymes in nature, namely, the monooxygenases and the peroxygenases. Both enzyme classes participate in substrate biodegradation or biosynthesis reactions in nature, but the P450 monooxygenases use dioxygen, while the peroxygenases take H2O2 in their catalytic cycle instead. By contrast to the P450 monooxygenases, the P450 peroxygenases do not require an external redox partner to deliver electrons during the catalytic cycle, and also no external proton source is needed. Therefore, they are fully self-sufficient, which affords them opportunities in biotechnological applications. One specific P450 peroxygenase, namely, P450 OleTJE, reacts with long-chain linear fatty acids through oxidative decarboxylation to form hydrocarbons and, as such, has been implicated as a suitable source for the biosynthesis of biofuels. Unfortunately, the reactions were shown to produce a considerable amount of side products originating from Cα and Cβ hydroxylation and desaturation. These product distributions were found to be strongly dependent on whether the substrate had substituents on the Cα and/or Cβ atoms. To understand the bifurcation pathways of substrate activation by P450 OleTJE leading to decarboxylation, Cα hydroxylation, Cβ hydroxylation and Cα-Cβ desaturation, we performed a computational study using 3-phenylpropionate and 2-phenylbutyrate as substrates. We set up large cluster models containing the heme, the substrate and the key features of the substrate binding pocket and calculated (using density functional theory) the pathways leading to the four possible products. This work predicts that the two substrates will react with different reaction rates due to accessibility differences of the substrates to the active oxidant, and, as a consequence, these two substrates will also generate different products. This work explains how the substrate binding pocket of P450 OleTJE guides a reaction to a chemoselectivity.

Density Functional Theory Study into the Reaction Mechanism of Isonitrile Biosynthesis by the Nonheme Iron Enzyme ScoE

The nonheme iron enzyme ScoE catalyzes the biosynthesis of an isonitrile substituent in a peptide chain. To understand details of the reaction mechanism we created a large active site cluster model of 212 atoms that contains substrate, the active oxidant and the first- and second-coordination sphere of the protein and solvent. Several possible reaction mechanisms were tested and it is shown that isonitrile can only be formed through two consecutive catalytic cycles that both use one molecule of dioxygen and α-ketoglutarate. In both cycles the active species is an iron(IV)-oxo species that in the first reaction cycle reacts through two consecutive hydrogen atom abstraction steps: first from the N–H group and thereafter from the C–H group to desaturate the NH-CH2 bond. The alternative ordering of hydrogen atom abstraction steps was also tested but found to be higher in energy. Moreover, the electronic configurations along that pathway implicate an initial hydride transfer followed by proton transfer. We highlight an active site Lys residue that is shown to donate charge in the transition states and influences the relative barrier heights and bifurcation pathways. A second catalytic cycle of the reaction of iron(IV)-oxo with desaturated substrate starts with hydrogen atom abstraction followed by decarboxylation to give isonitrile directly. The catalytic cycle is completed with a proton transfer to iron(II)-hydroxo to generate the iron(II)-water resting state. The work is compared with experimental observation and previous computational studies on this system and put in a larger perspective of nonheme iron chemistry.

Modeling the hydroxylation of estragole via human liver cytochrome P450

Natural compounds derived from plants are generally regarded safe and devoid of adverse effects. However, there are individual ingredients that possess toxic, genotoxic, and carcinogenic activities. These compounds when exposed at specific level become hazardous to health. Estragole (1-allyl-4-methoxybenzene) is a common component of spice plants. Its toxicity gets activated with the hydroxylation at benzylic carbon (C1') position by P450 enzymes present in the human liver. The present study grounds to explore the reaction mechanism of conversion of estragole to hydroxylated metabolite using computational methodology. Density functional theory (DFT)-based calculations were employed to explore the cytochrome P450-catalyzed mechanism at C1 position aliphatic hydroxylation of estragole. Overall reaction energy profile, electronic configuration, and 3D structure of all intermediates, transition states, and product complexes formed during the reaction along with their free energies were tried to be investigated.

Interactive Regulation between Aliphatic Hydroxylation and Aromatic Hydroxylation of Thaxtomin D in TxtC: A Theoretical Investigation

TxtC is an unusual bifunctional cytochrome P450 that is able to perform sequential aliphatic and aromatic hydroxylation of the diketopiperazine substrate thaxtomin D in two distinct sites to produce thaxtomin A. Though the X-ray structure of TxtC complexed with thaxtomin D revealed a binding mode for its aromatic hydroxylation, the preferential hydroxylation site is aliphatic C₁₄. It is thus intriguing to unravel how TxtC accomplishes such two-step catalytic hydroxylation on distinct aliphatic and aromatic carbons and why the aliphatic site is preferred in the hydroxylation step. In this work, by employing molecular docking and molecular dynamics (MD) simulation, we revealed that thaxtomin D could adopt two different conformations in the TxtC active site, which were equal in energy with either the aromatic C₂₀-H or aliphatic C₁₄-H pointing toward the active Cpd I oxyferryl moiety. Further ONIOM calculations indicated that the energy barrier for the rate-limiting hydroxylation step on the aliphatic C₁₄ site was 9.6 kcal/mol more favorable than that on the aromatic C₂₀ site. The hydroxyl group on the monohydroxylated intermediate thaxtomin B C₁₄ site formed hydrogen bonds with Ser280 and Thr385, which induced the l-Phe moiety to rotate around the C_β-C_γ bond of the 4-nitrotryptophan moiety. Thus, it adopted an energetically favorable conformation with aromatic C₂₀ adjacent to the oxyferryl moiety. In addition, the hydroxyl group induced solvent water molecules to enter the active site, which propelled thaxtomin B toward the heme plane and resulted in heme distortion. Based on this geometrical layout, the rate-limiting aromatic hydroxylation energy barrier decreased to 15.4 kcal/mol, which was comparable to that of the thaxtomin D aliphatic hydroxylation process. Our calculations indicated that heme distortion lowered the energy level of the lowest Cpd I α-vacant orbital, which promoted electron transfer in the rate-limiting thaxtomin B aromatic hydroxylation step in TxtC.

Mechanism-Based Inactivation of Cytochrome P450 Enzymes: Computational Insights

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

Kinetics and mechanism of OH-mediated degradation of three pentanols in the atmosphere

Pentanols as potential biofuels have attracted considerable interest, and thus it is of great importance to gain insights into their combustion and atmospheric chemistry.

Deciphering Reaction Determinants of Altered-Activity CYP2D6 Variants by Well-Tempered Metadynamics Simulation and QM/MM Calculations

The xenobiotic metabolizing enzyme CYP2D6 is the P450 cytochrome family member with the highest rate of polymorphism. This causes changes in the enzyme activity and specificity, which can ultimately lead to adverse reactions during drug treatment. To avoid or lower CYP-related toxicity risks, prediction of the most likely positions within a molecule where a metabolic reaction might occur is paramount. In order to obtain accurate predictions, it is crucial to understand all phenomena within the active site of the enzyme that contribute to an efficient substrate recognition and the subsequent catalytic reaction together with their relative weight within the overall thermodynamic context. This study aims to define the weight of the driving forces upon the C-H bond activation within CYP2D6 wild-type and a clinically relevant allelic variant with increased activity (CYP2D6*53) featuring two amino acid mutations in close vicinity of the heme. First, we investigated the steric and electrostatic complementarity of the substrate bufuralol using well-tempered metadynamics simulations with the aim to obtain the free energy profiles for each site of metabolism (SoM) within the different active sites. Second, the stereoelectronic complementarity was determined for each SoM within the two different active-site environments. Relying on the well-tempered metadynamics simulation energy profiles of each SoM, we identified the binding mode that was closest to the preferred transition-state geometry for efficient C-H bond activation. The binding modes were then used as starting structures for the quantum mechanics/molecular mechanics calculations performed to quantify the corresponding activation barriers. Our results show the relevance of the steric component in orienting the SoM in an energetically accessible position toward the heme. However, the corresponding intrinsic reactivity and electronic complementarity within the active site must be accurately evaluated in order to obtain a meaningful reaction prediction, from which the predominant SoM can be determined. The F120I mutation lowered the activation barrier for the major site and one of the minor SoMs. However, it had an impact neither on the CYP2D6 enantioselectivity preference of the oxidation reaction nor on the stereoselectivity from the substrate point of view.

Deciphering the origin of million-fold reactivity observed for the open core diiron [HO-FeIII-O-FeIV[double bond, length as m-dash]O]2+ species towards C-H bond activation: role of spin-states, spin-coupling, and spin-cooperation

High-valent metal-oxo species have been characterised as key intermediates in both heme and non-heme enzymes that are found to perform efficient aliphatic hydroxylation, epoxidation, halogenation, and dehydrogenation reactions. Several biomimetic model complexes have been synthesised over the years to mimic both the structure and function of metalloenzymes. The diamond-core [Fe2(μ-O)2] is one of the celebrated models in this context as this has been proposed as the catalytically active species in soluble methane monooxygenase enzymes (sMMO), which perform the challenging chemical conversion of methane to methanol at ease. In this context, a report of open core [HO(L)FeIII-O-FeIV(O)(L)]2+ (1) gains attention as this activates C-H bonds a million-fold faster compared to the diamond-core structure and has the dual catalytic ability to perform hydroxylation as well as desaturation with organic substrates. In this study, we have employed density functional methods to probe the origin of the very high reactivity observed for this complex and also to shed light on how this complex performs efficient hydroxylation and desaturation of alkanes. By modelling fifteen possible spin-states for 1 that could potentially participate in the reaction mechanism, our calculations reveal a doublet ground state for 1 arising from antiferromagnetic coupling between the quartet FeIV centre and the sextet FeIII centre, which regulates the reactivity of this species. The unusual stabilisation of the high-spin ground state for FeIV[double bond, length as m-dash]O is due to the strong overlap of with the orbital, reducing the antibonding interactions via spin-cooperation. The electronic structure features computed for 1 are consistent with experiments offering confidence in the methodology chosen. Further, we have probed various mechanistic pathways for the C-H bond activation as well as -OH rebound/desaturation of alkanes. An extremely small barrier height computed for the first hydrogen atom abstraction by the terminal FeIV[double bond, length as m-dash]O unit was found to be responsible for the million-fold activation observed in the experiments. The barrier height computed for -OH rebound by the FeIII-OH unit is also smaller suggesting a facile hydroxylation of organic substrates by 1. A strong spin-cooperation between the two iron centres also reduces the barrier for second hydrogen atom abstraction, thus making the desaturation pathway competitive. Both the spin-state as well as spin-coupling between the two metal centres play a crucial role in dictating the reactivity for species 1. By exploring various mechanistic pathways, our study unveils the fact that the bridged μ-oxo group is a poor electrophile for both C-H activation as well for -OH rebound. As more and more evidence is gathered in recent years for the open core geometry of sMMO enzymes, the idea of enhancing the reactivity via an open-core motif has far-reaching consequences.

Strategies to Mitigate the Bioactivation of Aryl Amines

The bioactivation of xenobiotics to yield reactive metabolites can lead to tolerability and toxicity concerns within a drug discovery program. Development of strategies for mitigating the metabolic liability of commonly encountered toxicophores, such as anilines, relies on an understanding of the relative tendency of these functionalities to undergo bioactivation. In this report, we present the first systematic study of the structure-activity relationships of the bioactivation of aryl amine fragments (molecular weight < 250 Da) using a glutathione (GSH) trapping assay in the presence of human liver microsomes and the reduced form of nicotinamide adenine dinucleotide phosphate. This study demonstrates that conversion of anilines to nitrogen-containing heteroarylamines results in a lower abundance of GSH conjugates in the order phenyl > pyrimidine ≈ pyridine > pyridazine. Introduction of electron-withdrawing functionality on the aromatic ring had a less pronounced effect on the extent of GSH conjugation. Examination of more drug-like compounds sourced from in-house drug discovery programs revealed similar trends in bioactivation between matched pairs containing (hetero)aryl amines. This study provides medicinal chemists with insights and qualitative guidance for the minimization of risks related to aryl amine metabolism.

Bioengineering of Cytochrome P450 OleTJE: How Does Substrate Positioning Affect the Product Distributions?

The cytochromes P450 are versatile enzymes found in all forms of life. Most P450s use dioxygen on a heme center to activate substrates, but one class of P450s utilizes hydrogen peroxide instead. Within the class of P450 peroxygenases, the P450 OleT_JE isozyme binds fatty acid substrates and converts them into a range of products through the α-hydroxylation, β-hydroxylation and decarboxylation of the substrate. The latter produces hydrocarbon products and hence can be used as biofuels. The origin of these product distributions is unclear, and, as such, we decided to investigate substrate positioning in the active site and find out what the effect is on the chemoselectivity of the reaction. In this work we present a detailed computational study on the wild-type and engineered structures of P450 OleT_JE using a combination of density functional theory and quantum mechanics/molecular mechanics methods. We initially explore the wild-type structure with a variety of methods and models and show that various substrate activation transition states are close in energy and hence small perturbations as through the protein may affect product distributions. We then engineered the protein by generating an in silico model of the double mutant Asn242Arg/Arg245Asn that moves the position of an active site Arg residue in the substrate-binding pocket that is known to form a salt-bridge with the substrate. The substrate activation by the iron(IV)-oxo heme cation radical species (Compound I) was again studied using quantum mechanics/molecular mechanics (QM/MM) methods. Dramatic differences in reactivity patterns, barrier heights and structure are seen, which shows the importance of correct substrate positioning in the protein and the effect of the second-coordination sphere on the selectivity and activity of enzymes.

Cross-linking of aromatic phenolate groups by cytochrome P450 enzymes: a model for the biosynthesis of vancomycin by OxyB

Aromatic cross-linking by cytochrome P450 enzymes was studied computationally. P450 Compound I rapidly abstracts two weak phenolic H-atoms that link up via a rate-determining C–O bond formation.

Sluggish reactivity by a nonheme iron(iv)-tosylimido complex as compared to its oxo analogue

A comparative spectroscopic and computational study of reactivity between ferryl-tosylimido and ferryl-oxo complexes of two biomimetic model systems. The Fe(iv)-tosylimido complex was found to be sluggish in comparison to its fellow oxo counterpart.

Selective Hydrogen Atom Abstraction from Dihydroflavonol by a Nonheme Iron Center Is the Key Step in the Enzymatic Flavonol Synthesis and Avoids Byproducts

The plant non-heme iron dioxygenase flavonol synthase performs a regioselective desaturation reaction as part of the biosynthesis of the signaling molecule flavonol that triggers the growing of leaves and flowers. These compounds also have health benefits for humans. Desaturation of aliphatic compounds generally proceeds through two consecutive hydrogen atom abstraction steps from two adjacent carbon atoms and in nature often is performed by a high-valent iron(IV)-oxo species. We show that the order of the hydrogen atom abstraction steps, however, is opposite of those expected from the C-H bond strengths in the substrate and determines the product distributions. As such, flavonol synthase follows a negative catalysis mechanism. Using density functional theory methods on large active-site model complexes, we investigated pathways for desaturation and hydroxylation by an iron(IV)-oxo active-site model. Contrary to thermochemical predictions, we find that the oxidant abstracts the hydrogen atom from the strong C²-H bond rather than the weaker C³-H bond of the substrate first. We analyze the origin of this unexpected selective hydrogen atom abstraction pathway and find that the alternative C³-H hydrogen atom abstraction would be followed by a low-energy and competitive substrate hydroxylation mechanism hence, should give considerable amount of byproducts. Our computational modeling studies show that substrate positioning in flavonol synthase is essential, as it guides the reactivity to a chemo- and regioselective substrate desaturation from the C²-H group, leading to desaturation products efficiently.

Properties and reactivity of μ-nitrido-bridged dimetal porphyrinoid complexes: how does ruthenium compare to iron?

Methane hydroxylation by metal-oxo oxidants is one of the Holy Grails in biomimetic and biotechnological chemistry. The only enzymes known to perform this reaction in Nature are iron-containing soluble methane monooxygenase and copper-containing particulate methane monooxygenase. Furthermore, few biomimetic iron-containing oxidants have been designed that can hydroxylate methane efficiently. Recent studies reported that μ-nitrido-bridged diiron(IV)-oxo porphyrin and phthalocyanine complexes hydroxylate methane to methanol efficiently. To find out whether the reaction rates are enhanced by replacing iron by ruthenium, we performed a detailed computational study. Our work shows that the μ-nitrido-bridged diruthenium(IV)-oxo reacts with methane via hydrogen atom abstraction barriers that are considerably lower in energy (by about 5 kcal mol^‒1) as compared to the analogous diiron(IV)-oxo complex. An analysis of the electronic structure implicates similar spin and charge distributions for the diiron(IV)-oxo and diruthenium(IV)-oxo complexes, but the strength of the O‒H bond formed during the reaction is much stronger for the latter. As such a larger hydrogen atom abstraction driving force for the Ru complex than for the Fe complex is found, which should result in higher reactivity in the oxidation of methane.

Conformational turn triggers regio-selectivity in the bioactivation of thiophene-contained compounds mediated by cytochrome P450

In the present work, we performed Density Functional Theory calculations to explore the bioactivation mechanism of thiophene-containing molecules mediated by P450s. For this purpose, relatively large size compounds, 2,5-diaminothiophene derivatives were selected particularly for this investigation. Here we found the successive regio-selectivity triggered by conformational turn played a significant role in the occurrence of bioactivation. 2,5-Diaminothiophene was oxidized to a 2,5-diimine thiophene-reactive intermediate by Compound I (Cpd I) through successive activations of two N-H bonds (H3-N11 and H1-N6). This reaction exhibited three special characteristics: (1) self-controlled regio-selectivity during the oxidation process. There was a large scale of conformational turn in the abstraction of the first H atom which triggers the selection of the second H for abstraction. (2) Proton-shuttle mechanism. In high spin (HS) state, proton-shuttle mechanism was observed for the abstraction of the second H atom. (3) Spin-selective manner. In protein environment, the energy barrier in HS state was much lower than that in low spin state. The novel proposed bioactivation mechanism of 2,5-diaminothiophene compounds can help us in rational design of thiophene-contained drugs avoiding the occurrence of bioactivation.

Flavonol biosynthesis by nonheme iron dioxygenases: A computational study into the structure and mechanism

Plants produce flavonol compounds for vital functions regarding plant growth, fruit and flower colouring as well as fruit ripening processes. Several of these biosynthesis steps are stereo- and regioselective and are being carried out by nonheme iron enzymes. Using density functional theory calculations on a large active site model complex of flavanone-3β-hydroxylase (FHT), we established the mechanism for conversion of naringenin to its dihydroflavonol, which is a key step in the mechanism of flavonol biosynthesis. The reaction starts with dioxygen binding to the iron(II) centre and a reaction with α-ketoglutarate co-substrate gives succinate, an iron(IV)-oxo species and CO₂ with large exothermicity and small reaction barriers. The rate-determining reaction step in the mechanism; however, is hydrogen atom abstraction of an aliphatic CH bond by the iron(IV)-oxo species. We identify a large kinetic isotope effect for the replacement of the transferring hydrogen atom by deuterium. In a final step the OH and substrate radicals combine to form the alcohol product with a barrier of several kcal mol^-1. We show that the latter is the result of geometric constraints in the active site pocket. Furthermore, the calculations show that a weak tertiary CH bond is shielded from the iron(IV)-oxo species in the substrate binding position and therefore the enzyme is able to activate a stronger CH bond. As such, the flavanone-3β-hydroxylase enzyme reacts regioselectively with one specific CH bond of naringenin by avoiding activation of weaker bonds through tight substrate and oxidant positioning.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

Computational analysis of non-heme iron-oxo formation by direct NO release in nitrite reduction

A direct NO-releasing reaction of nitrite catalyzed by [N(afa^Cy)₃Fe(OTf)]⁺ (afa (azafulvene-amine); OTf (trifluoromethanesulfonate); Cy (cyclohexyl)) was investigated using density functional theory (DFT) with D3 dispersion correction. The complex featured a secondary coordination sphere that facilitated the formation of the iron-oxo product [N(afa^Cy)₃FeO]⁺ with three (Fe)OH-N hydrogen bonds. As a high-spin iron(ii), the O-binding initial intermediate Fe(O)-nitrito was thermodynamically favorable in the S = 2 state. The cleavage of the (Fe)O-NO bond was performed by a β-electron shift to produce Fe(iii)-O by electron rearrangement in the S = 5/2 state. The different electron configurations are responsible for the structural properties, the valence of iron in the complexes, and the pathways of the reactions. Moreover, the two different H-bonds, (Fe)OH-N and (Fe)O-HN (by O-protonation), in the product complexes played a role in determining the reaction channels by impacting the N-H bond rotation. Thus, an exothermic sequence of conversions Fe(ii) → Fe(iii)-O → Fe(iii)-OH → Fe(iii)-O was established for the targeted product formation. This process provided a clue to build two key intermediates, iron-oxo and iron-hydroxo, in a variety of biological and synthetic systems. The results of this study are in agreement with experimental observations and describe the roles of H-bonding in nitrite reduction catalyzed by the non-heme iron complex.

Interplay Between Steric and Electronic Effects: A Joint Spectroscopy and Computational Study of Nonheme Iron(IV)-Oxo Complexes

Iron is an essential element in nonheme enzymes that plays a crucial role in many vital oxidative transformations and metabolic reactions in the human body. Many of those reactions are regio- and stereospecific and it is believed that the selectivity is guided by second-coordination sphere effects in the protein. Here, results are shown of a few engineered biomimetic ligand frameworks based on the N4Py (N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine) scaffold and the second-coordination sphere effects are studied. For the first time, selective substitutions in the ligand framework have been shown to tune the catalytic properties of the iron(IV)-oxo complexes by regulating the steric and electronic factors. In particular, a better positioning of the oxidant and substrate in the rate-determining transition state lowers the reaction barriers. Therefore, an optimum balance between steric and electronic factors mediates the ideal positioning of oxidant and substrate in the rate-determining transition state that affects the reactivity of high-valent reaction intermediates.