CH bond activation in heme proteins: the role of thiolate ligation in cytochrome P450

An Oxoiron(IV) Complex Supported by an N-Alkylated Cyclam Ligand System Containing a Pendant Alcohol Moiety

The effect of a pendant neutral alcohol moiety in the N-alkylated cyclam (1,4,8,11-tetraazacyclotetradecane) ligand backbone is examined for the non-heme mononuclear oxoiron(IV) unit in [FeIV(Osyn)(TMC-HOR)(NCCH3)]2+ (1-syn) (TMC-HOR=2-(4,8,11-trimethyl-1,4,8,11-tetraazacyclotetradecan1-yl)ethan-1-ol). Unlike in the related [FeIV(Oanti)(TMC-SR)]+ (3-anti) (TMC-SR=1-mercaptoethyl-4,8,11-trimethyl-1,4,8,11-tetraazacyclotetradecane) complex, bearing an axial mono-anionic thiolate ligand trans to the oxo unit, the alcohol moiety in 1-syn stays protonated and does not axially coordinate to iron. The protonation of the alcohol moiety is a prerequisite for the stabilization of the oxoiron(IV) core; it presumably serves as a hydrogen bonding donor to the oxoiron(IV) unit, which is positioned syn to the three methyl groups. Comparative reactivity studies reveal 1-syn to be a stronger hydrogen atom abstraction but weaker oxygen atom transfer agent relative to the [FeIV(Osyn)(TMC)(NCCH3)]2+ (2-syn) complex, bearing the N-tetramethylated cyclam (TMC) ligand.

Insights into Active Site Cysteine Residues in Mycobacterium tuberculosis Enzymes: Potential Targets for Anti-Tuberculosis Intervention

Cysteine, a semi-essential amino acid, is found in the active site of a number of vital enzymes of the bacterium Mycobacterium tuberculosis (Mtb) and in particular those that relate to its survival, adaptability and pathogenicity. Mtb is the causative agent of tuberculosis, an infectious disease that affects millions of people globally. Common anti-tuberculosis targets are focused on immobilizing a vital cysteine amino acid residue in enzymes that plays critical roles in redox and non-redox catalysis, the modulation of the protein, enzyme activity, protein structure and folding, metal coordination, and posttranslational modifications of newly synthesized proteins. This review examines five Mtb enzymes that contain an active site cysteine residue and are considered as key targets for anti-tuberculosis drugs, namely alkyl hydroperoxide reductase (AhpC), dihydrolipoamide dehydrogenase (Lpd), aldehyde dehydrogenase (ALDH), methionine aminopeptidase (MetAP) and cytochromes P450. AhpC and Lpd protect Mtb against oxidative and nitrosative stress, whereas AhpC neutralizes peroxide/peroxynitrite substrates with two active site cysteine residues. Mtb ALDH detoxifies aldehydes, using a nucleophilic active site cysteine to form an oxyanion thiohemiacetal intermediate, whereas MtMetAP's active site cysteine is essential for substrate recognition. The P450s metabolize various endogenous and exogenous compounds. Targeting these critical active site cysteine residues could disrupt enzyme functions, presenting a promising avenue for developing anti-mycobacterial agents.

Computational Insights into Hydrogen Atom Transfer Mediators in C-H Activation Catalysis of Nonheme Fe(IV)O Complexes

This study presents a detailed density functional theory (DFT) investigation into the mechanism and energetics of C-H activations catalyzed by bioinspired Fe(IV)O complexes, particularly in the presence of N-hydroxy mediators. The findings show that these mediators significantly enhance the reactivity of the iron-oxo complex. The study examines three substrates with varying bond dissociation energies─ethylbenzene, cyclohexane, and cyclohexadiene─alongside the [Fe(IV)O(N4Py)]2+ complex. Mediators N-hydroxyphthalimide (NHPI) and N-hydroxyquinolinimide (NHQI) were chosen for their strong oxidative abilities. The results reveal that NO-H bond cleavage in N-hydroxy compounds occurs more readily than C-H bond cleavage in hydrocarbons, as supported by the Marcus cross-relation applied to H-abstraction. This leads to the rapid formation of aminoxyl radicals, which are more reactive than Fe(IV)O species, lowering the activation energy and enhancing the reaction rate. The C-H bond activation aligns with the Bell-Evans-Polanyi principle, correlating the activation energy with the substrate bond dissociation energy. The investigation reveals that the mediator pathway is favored both thermodynamically and kinetically. Additionally, distortion energy provides a compelling explanation for the observed reactivity trends, further highlighting NHQI's superior efficiency compared to NHPI. Additionally, quantum mechanical tunneling plays a significant role, as evidenced by the computed kinetic isotope effect, which matches experimental data.

Tutorial Review on the Set-Up and Running of Quantum Mechanical Cluster Models for Enzymatic Reaction Mechanisms

Enzymes turnover substrates into products with amazing efficiency and selectivity and as such have great potential for use in biotechnology and pharmaceutical applications. However, details of their catalytic cycles and the origins surrounding the regio- and chemoselectivity of enzymatic reaction processes remain unknown, which makes the engineering of enzymes and their use in biotechnology challenging. Computational modelling can assist experimental work in the field and establish the factors that influence the reaction rates and the product distributions. A popular approach in modelling is the use of quantum mechanical cluster models of enzymes that take the first- and second coordination sphere of the enzyme active site into consideration. These QM cluster models are widely applied but often the results obtained are dependent on model choice and model selection. Herein, we show that QM cluster models can give highly accurate results that reproduce experimental product distributions and free energies of activation within several kcal mol-1, regarded that large cluster models with >300 atoms are used that include key hydrogen bonding interactions and charged residues. In this tutorial review, we give general guidelines on the set-up and applications of the QM cluster method and discuss its accuracy and reproducibility. Finally, several representative QM cluster model examples on metal-containing enzymes are presented, which highlight the strength of the approach.

Exploring the influence of H-bonding and ligand constraints on thiolate ligated non-heme iron mediated dioxygen activation

Converting triplet dioxygen into a powerful oxidant is fundamentally important to life. The study reported herein quantitatively examines the formation of a well-characterized, reactive, O2-derived thiolate ligated FeIII-superoxo using low-temperature stopped-flow kinetics. Comparison of the kinetic barriers to the formation of this species via two routes, involving either the addition of (a) O2 to [FeII(S2 Me2N3(Pr,Pr))] (1) or (b) superoxide to [FeIII(S2 Me2N3(Pr,Pr))]+ (3) is shown to provide insight into the mechanism of O2 activation. Route (b) was shown to be significantly slower, and the kinetic barrier 14.9 kJ mol-1 higher than route (a), implying that dioxygen activation involves inner-sphere, as opposed to outer sphere, electron transfer from Fe(ii). H-bond donors and ligand constraints are shown to dramatically influence O2 binding kinetics and reversibility. Dioxygen binds irreversibly to [FeII(S2 Me2N3(Pr,Pr))] (1) in tetrahydrofuran, but reversibly in methanol. Hydrogen bonding decreases the ability of the thiolate sulfur to stabilize the transition state and the FeIII-superoxo, as shown by the 10 kJ mol-1 increase in the kinetic barrier to O2 binding in methanol vs. tetrahydrofuran. Dioxygen release from [FeIII(S2 Me2N3(Pr,Pr))O2] (2) is shown to be 24 kJ mol-1 higher relative to previously reported [FeIII(SMe2N4(tren))(O2)]+ (5), the latter of which contains a more flexible ligand. These kinetic results afford an experimentally determined reaction coordinate that illustrates the influence of H-bonding and ligand constraints on the kinetic barrier to dioxygen activation an essential step in biosynthetic pathways critical to life.

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.

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

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

Correlating Structure with Spectroscopy in Ascorbate Peroxidase Compound II

Structural and spectroscopic investigations of compound II in ascorbate peroxidase (APX) have yielded conflicting conclusions regarding the protonation state of the crucial Fe(IV) intermediate. Neutron diffraction and crystallographic data support an iron(IV)-hydroxo formulation, whereas Mössbauer, X-ray absorption (XAS), and nuclear resonance vibrational spectroscopy (NRVS) studies appear consistent with an iron(IV)-oxo species. Here we examine APX with spectroscopy-oriented QM/MM calculations and extensive exploration of the conformational space for both possible formulations of compound II. We establish that irrespective of variations in the orientation of a vicinal arginine residue and potential reorganization of proximal water molecules and hydrogen bonding, the Fe-O distances for the oxo and hydroxo forms consistently fall within distinct, narrow, and nonoverlapping ranges. The accuracy of geometric parameters is validated by coupled-cluster calculations with the domain-based local pair natural orbital approach, DLPNO-CCSD(T). QM/MM calculations of spectroscopic properties are conducted for all structural variants, encompassing Mössbauer, optical, X-ray absorption, and X-ray emission spectroscopies and NRVS. All spectroscopic observations can be assigned uniquely to an Fe(IV)═O form. A terminal hydroxy group cannot be reconciled with the spectroscopic data. Under no conditions can the Fe(IV)═O distance be sufficiently elongated to approach the crystallographically reported Fe-O distance. The latter is consistent only with a hydroxo species, either Fe(IV) or Fe(III). Our findings strongly support the Fe(IV)═O formulation of APX-II and highlight unresolved discrepancies in the nature of samples used across different experimental studies.

High-valent nonheme Fe(IV)O/Ru(IV)O complexes catalyze C-H activation reactivity and hydrogen tunneling: a comparative DFT investigation

A comprehensive density functional theory investigation has been presented towards the comparison of the C-H activation reactivity between high-valent iron-oxo and ruthenium-oxo complexes. A total of four compounds, e.g., [Ru(IV)O(tpy-dcbpy)] (1), [Fe(IV)O(tpy-dcbpy)] (1'), [Ru(IV)O(TMCS)] (2), and [Fe(IV)O(TMCS)] (2'), have been considered for this investigation. The macrocyclic ligand framework tpy(dcbpy) implies tpy = 2,2':6',2''-terpyridine, dcbpy = 5,5'-dicarboxy-2,2'-bipyridine, and TMCS is TMC with an axially tethered -SCH2CH2 group. Compounds 1 and 2' are experimentally synthesized standard complexes with Ru and Fe, whereas compounds 1' and 2 were considered to keep the macrocycle intact when switching the central metal atom. Three reactants including benzyl alcohol, ethyl benzene, and dihydroanthracene were selected as substrates for C-H activation. It is noteworthy to mention that Fe(IV)O complexes exhibit higher reactivity than those of their Ru(IV)O counterparts. Furthermore, regardless of the central metal, the complex featuring a tpy-dcbpy macrocycle demonstrates higher reactivity than that of TMCS. Here, a thorough analysis of the reactivity-controlling characteristics-such as spin state, steric factor, distortion energy, energy of the electron acceptor orbital, and quantum mechanical tunneling-was conducted. Fe(IV)O exhibits the exchanged enhanced two-state-reactivity with the quintet reactive state, whereas Ru(IV)O has only a triplet reactive state. Both the distortion energy and acceptor orbital energy are low in the case of Fe(IV)O supporting its higher reactivity. All the investigated C-H activation processes involve a significant contribution from hydrogen tunneling, which is more pronounced in the case of Ru, although it cannot alter the reactivity pattern. Furthermore, it has also been found that, independent of the central metal, aliphatic hydroxylation is always preferable to aromatic hydroxylation. Overall, this work is successful in establishing and investigating the cause of enzymes' natural preference for Fe over Ru as a cofactor for C-H activation enzymes.

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.

Influence of Equatorial Co-Ligands on the Reactivity of LFeIIIOIPh

Previous biomimetic studies clearly proved that equatorial ligands significantly influence the redox potential and thus the stability/reactivity of biologically important oxoiron intermediates; however, no such studies were performed on FeIIIOIPh species. In this study, the influence of substituted pyridine co-ligands on the reactivity of iron(III)-iodosylbenzene adduct has been investigated in sulfoxidation and epoxidation reactions. Selective oxidation of thioanisole, cis-cyclooctene, and cis- and trans-stilbene in the presence of a catalytic amount of [FeII(PBI)3](OTf)2 with PhI(OAc)2 provide products in good to excellent yields through an FeIIIOIPh intermediate depending on the co-ligand (4R-Py) used. Several mechanistic studies were performed to gain more insight into the mechanism of oxygen atom transfer (OAT) reactions to support the reactive intermediate and investigate the effect of the equatorial co-ligands. Based on competitive experiments, including a linear free-energy relationship between the relative reaction rates (logkrel) and the σp (4R-Py) parameters, strong evidence has been observed for the electrophilic character of the reactive species. The presence of the [(PBI)2(4R-Py)FeIIIOIPh]3+ intermediates and the effect of the co-ligands was also supported by UV-visible measurements, including the color change from red to green and the hypsochromic shifts in the presence of co-ligands. This is another indication that the title iron(III)-iodosylbenzene adduct is able to oxygenate sulfides and alkenes before it is transformed into the oxoiron form by cleavage of the O-I bond.

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.

Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides

Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations.

Improved Bioproduction of the Nylon 12 Monomer by Combining the Directed Evolution of P450 and Enhancing Heme Synthesis

The nylon 12 (PA12) monomer ω-aminododecanoic acid (ω-AmDDA) could be synthesized from lauric acid (DDA) through multi-enzyme cascade transformation using engineered E. coli, with the P450 catalyzing terminal hydroxylation of DDA as a rate-limiting enzyme. Its activity is jointly determined by the heme domain and the reductase domain. To obtain a P450 mutant with higher activity, directed evolution was conducted using a colorimetric high-throughput screening (HTS) system with DDA as the real substrate. After two rounds of directed evolution, a positive double-site mutant (R14R/D629G) with 90.3% higher activity was obtained. Molecular docking analysis, kinetic parameter determination and protein electrophoresis suggested the improved soluble expression of P450 resulting from the synonymous mutation near the N-terminus and the shortened distance of the electron transfer between FMN and FAD caused by D629G mutation as the major reasons for activity improvement. The significantly increased k_cat and unchanged K_m provided further evidence for the increase in electron transfer efficiency. Considering the important role of heme in P450, its supply was strengthened by the metabolic engineering of the heme synthesis pathway. By combining P450-directed evolution and enhancing heme synthesis, 2.02 ± 0.03 g/L of ω-AmDDA was produced from 10 mM DDA, with a yield of 93.6%.

Construction of Biocatalysts Using the P450 Scaffold for the Synthesis of Indigo from Indole

With the increasing demand for blue dyes, it is of vital importance to develop a green and efficient biocatalyst to produce indigo. This study constructed a hydrogen peroxide-dependent catalytic system for the direct conversion of indole to indigo using P450BM3 with the assistance of dual-functional small molecules (DFSM). The arrangements of amino acids at 78, 87, and 268 positions influenced the catalytic activity. F87G/T268V mutant gave the highest catalytic activity with k_cat of 1402 min^-1 and with a yield of 73%. F87A/T268V mutant was found to produce the indigo product with chemoselectivity as high as 80%. Moreover, F87G/T268A mutant was found to efficiently catalyze indole oxidation with higher activity (k_cat/K_m = 1388 mM^-1 min^-1) than other enzymes, such as the NADPH-dependent P450BM3 (2.4-fold), the Ngb (32-fold) and the Mb (117-fold). Computer simulation results indicate that the arrangements of amino acid residues in the active site can significantly affect the catalytic activity of the protein. The DFSM-facilitated P450BM3 peroxygenase system provides an alternative, simple approach for a key step in the bioproduction of indigo.

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

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

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

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

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.

Electronic Structure and Reactivity of Dioxygen-Derived Aliphatic Thiolate-Ligated Fe-Peroxo and Fe(IV) Oxo Compounds

Herein, we examine the electronic and geometric structural properties of O₂-derived aliphatic thiolate-ligated Fe-peroxo, Fe-hydroxo, and Fe(IV) oxo compounds. The latter cleaves strong C-H bonds (96 kcal mol^-1) on par with the valine C-H bond cleaved by isopencillin N synthase (IPNS). Stopped-flow kinetics studies indicate that the barrier to O₂ binding to [Fe^II(S^Me2N₄(tren))]⁺ (3) is extremely low (E_a = 36(2) kJ mol^-1), as theoretically predicted for IPNS. Dioxygen binding to 3 is shown to be reversible, and a superoxo intermediate, [Fe^III(S^Me2N₄(tren))(O₂)]⁺ (6), forms in the first 25 ms of the reaction at -40 °C prior to the rate-determining (E_a = 46(2) kJ mol^-1) formation of peroxo-bridged [(S^Me2N₄(tren))Fe(III)]₂(μ-O₂)²⁺ (7). A log(k_obs) vs log([Fe]) plot for the formation of 7 is consistent with the second-order dependence on iron, and H₂O₂ assays are consistent with a 2:1 ratio of Fe/H₂O₂. Peroxo 7 is shown to convert to ferric-hydroxo [Fe^III(S^Me2N(tren))(OH)]⁺ (9, g_⊥ = 2.24, g_∥ = 1.96), the identity of which was determined via its independent synthesis. Rates of the conversion 7 → 9 are shown to be dependent on the X-H bond strength of the H-atom donor, with a k_H/k_D = 4 when CD₃OD is used in place of CH₃OH as a solvent. A crystallographically characterized cis thiolate-ligated high-valent iron oxo, [Fe^IV(O)(S^Me2N₄(tren))]⁺ (11), is shown to form en route to hydroxo 9. Electronic structure calculations were shown to be consistent with 11 being an S = 1 Fe(IV)═O with an unusually high ν_Fe-O stretching frequency at 918 cm^-1 in line with the extremely short Fe-O bond (1.603(7) Å).

The Cytochrome P450 OxyA from the Kistamicin Biosynthesis Cyclization Cascade is Highly Sensitive to Oxidative Damage

Cytochrome P450 enzymes (P450s) are a superfamily of monooxygenases that utilize a cysteine thiolate–ligated heme moiety to perform a wide range of demanding oxidative transformations. Given the oxidative power of the active intermediate formed within P450s during their active cycle, it is remarkable that these enzymes can avoid auto-oxidation and retain the axial cysteine ligand in the deprotonated—and thus highly acidic—thiolate form. While little is known about the process of heme incorporation during P450 folding, there is an overwhelming preference for one heme orientation within the P450 active site. Indeed, very few structures to date contain an alternate heme orientation, of which two are OxyA homologs from glycopeptide antibiotic (GPA) biosynthesis. Given the apparent preference for the unusual heme orientation shown by OxyA enzymes, we investigated the OxyA homolog from kistamicin biosynthesis (OxyA_kis), which is an atypical GPA. We determined that OxyA_kis is highly sensitive to oxidative damage by peroxide, with both UV and EPR measurements showing rapid bleaching of the heme signal. We determined the structure of OxyA_kis and found a mixed population of heme orientations present in this enzyme. Our analysis further revealed the possible modification of the heme moiety, which was only present in samples where the alternate heme orientation was present in the protein. These results suggest that the typical heme orientation in cytochrome P450s can help prevent potential damage to the heme—and hence deactivation of the enzyme—during P450 catalysis. It also suggests that some P450 enzymes involved in GPA biosynthesis may be especially prone to oxidative damage due to the heme orientation found in their active sites.

Local Charge Distributions, Electric Dipole Moments, and Local Electric Fields Influence Reactivity Patterns and Guide Regioselectivities in α-Ketoglutarate-Dependent Non-heme Iron Dioxygenases

Non-heme iron dioxygenases catalyze vital processes for human health related to the biosynthesis of essential products and the biodegradation of toxic metabolites. Often the natural product biosyntheses by these non-heme iron dioxygenases is highly regio- and chemoselective, which are commonly assigned to tight substrate-binding and positioning. However, recent high-level computational modeling has shown that substrate-binding and positioning is only part of the story and long-range electrostatic interactions can play a major additional role.In this Account, we review and summarize computational viewpoints on the high regio- and chemoselectivity of α-ketoglutarate-dependent non-heme iron dioxygenases and how external perturbations affect the catalysis. In particular, studies from our groups have shown that often a regioselectivity in enzymes can be accomplished by stabilization of the rate-determining transition state for the reaction through external charges, electric dipole moments, or local electric field effects. Furthermore, bond dissociation energies in molecules are shown to be influenced by an electric field effect, and through targeting a specific bond in an electric field, this can lead to an unusually specific reaction. For instance, in the carbon-induced starvation protein, we studied two substrate-bound conformations and showed that regardless of what C-H bond of the substrate is closest to the iron(IV)-oxo oxidant, the lowest hydrogen atom abstraction barrier is always for the pro-S C²-H abstraction due to an induced dipole moment of the protein that weakens this bond. In another example of the hygromycin biosynthesis enzyme, an oxidative ring-closure reaction in the substrate forms an ortho-δ-ester ring. Calculations on this enzyme show that the selectivity is guided by a protonated lysine residue in the active site that, through its positive charge, triggers a low energy hydrogen atom abstraction barrier. A final set of examples in this Account discuss the viomycin biosynthesis enzyme and the 2-(trimethylammonio)ethylphosphonate dioxygenase (TmpA) enzyme. Both of these enzymes are shown to possess a significant local dipole moment and local electric field effect due to charged residues surrounding the substrate and oxidant binding pockets. The protein dipole moment and local electric field strength changes the C-H bond strengths of the substrate as compared to the gas-phase triggers the regioselectivity of substrate activation. In particular, we show that in the gas phase and in a protein environment C-H bond strengths are different due to local electric dipole moments and electric field strengths. These examples show that enzymes have an intricately designed structure that enables a chemical reaction under ambient conditions through the positioning of positively and negatively charged residues that influence and enhance reaction mechanisms. These computational insights create huge possibilities in bioengineering to apply local electric field and dipole moments in proteins to achieve an unusual selectivity and specificity and trigger a fit-for-purpose biocatalyst for unique biotransformations.

NRVS investigation of ascorbate peroxidase compound II: Observation of Iron(IV)oxo stretching

The protonation state of ascorbate peroxidase compound II (APX-II) has been a subject of debate. A combined X-ray/neutron crystallographic study reported that APX-II is best described as an iron(IV)hydroxide species with an FeO distance of 1.88 Å (Kwon, et al. Nat Commun2016, 7, 13,445), while X-ray absorption spectroscopy (XAS) experiments (utilizing extended X-ray absorption fine structure (EXAFS) and pre-edge analyses) indicate APX-II is an authentic iron(IV)oxo species with an FeO distance 1.68 Å (Ledray, et al. Journal of the American Chemical Society2020,142, 20,419). Previous debates concerning ferryl protonation states have been resolved through the application of Badger's rule, which correlates FeO bond distances with FeO vibrational frequencies. To obtain the required vibrational data, we have collected Nuclear Resonance Vibrational Spectroscopy (NRVS) data for APX-II. We observe a broad vibrational feature in the range associated with iron(IV)oxo stretching (700-800 cm-1). This feature appears to have two peaks at 732 cm-1 and 770 cm-1, corresponding to FeO distances of 1.69 and 1.67 Å, respectively. The broad vibrational envelope and the presence of multiple resonances could reflect a distribution of hydrogen bonding interactions within the active-site pocket.

Efficient Base-Free Aqueous Reforming of Methanol Homogeneously Catalyzed by Ruthenium Exhibiting a Remarkable Acceleration by Added Catalytic Thiol

Production of H2 by methanol reforming is of particular interest due the low cost, ready availability, and high hydrogen content of methanol. However, most current methods either require very high temperatures and pressures or strongly rely on the utilization of large amounts of base. Here we report an efficient, base-free aqueous-phase reforming of methanol homogeneously catalyzed by an acridine-based ruthenium pincer complex, the activity of which was unexpectedly improved by a catalytic amount of a thiol additive. The reactivity of this system is enhanced by nearly 2 orders of magnitude upon addition of the thiol, and it can maintain activity for over 3 weeks, achieving a total H2 turnover number of over 130 000. On the basis of both experimental and computational studies, a mechanism is proposed which involves outer-sphere dehydrogenations promoted by a unique ruthenium complex with thiolate as an assisting ligand. The current system overcomes the need for added base in homogeneous methanol reforming and also highlights the unprecedented acceleration of catalytic activity of metal complexes achieved by the addition of a catalytic amount of thiol.

Nitrene transfers mediated by natural and artificial iron enzymes

Amines are ubiquitous in biology and pharmacy. As a consequence, introducing N functionalities in organic molecules is attracting strong continuous interest. The past decade has witnessed the emergence of very efficient and selective catalytic systems achieving this goal thanks to engineered hemoproteins. In this review, we examine how these enzymes have been engineered focusing rather on the rationale behind it than the methodology employed. These studies are put in perspective with respect to in vitro and in vivo nitrene transfer processes performed by cytochromes P450. An emphasis is put on mechanistic aspects which are confronted to current molecular knowledge of these reactions. Forthcoming developments are delineated.

Rejigging Electron and Proton Transfer to Transition between Dioxygenase, Monooxygenase, Peroxygenase, and Oxygen Reduction Activity: Insights from Bioinspired Constructs of Heme Enzymes

Nature has employed heme proteins to execute a diverse set of vital life processes. Years of research have been devoted to understanding the factors which bias these heme enzymes, with all having a heme cofactor, toward distinct catalytic activity. Among them, axial ligation, distal super structure, and substrate binding pockets are few very vividly recognized ones. Detailed mechanistic investigation of these heme enzymes suggested that several of these enzymes, while functionally divergent, use similar intermediates. Furthermore, the formation and decay of these intermediates depend on proton and electron transfer processes in the enzyme active site. Over the past decade, work in this group, using in situ surface enhanced resonance Raman spectroscopy of synthetic and biosynthetic analogues of heme enzymes, a general idea of how proton and electron transfer rates relate to the lifetime of different O2 derived intermediates has been developed. These findings suggest that the enzymatic activities of all these heme enzymes can be integrated into one general cycle which can be branched out to different catalytic pathways by regulating the lifetime and population of each of these intermediates. This regulation can further be achieved by tuning the electron and proton transfer steps. By strategically populating one of these intermediates during oxygen reduction, one can navigate through different catalytic processes to a desired direction by altering proton and electron transfer steps.

Increasing reactivity by incorporating π-acceptor ligands into coordinatively unsaturated thiolate-ligated iron(II) complexes

Reported herein is the structural, spectroscopic, redox, and reactivity properties of a series of iron complexes containing both a π-donating thiolate, and π-accepting N-heterocycles in the coordination sphere, in which we systematically vary the substituents on the N-heterocycle, the size of the N-heterocycle, and the linker between the imine nitrogen and tertiary amine nitrogen. In contrast to our primary amine/thiolate-ligated Fe(II) complex, [Fe^II(S^Me2N₄(tren))]⁺ (1), the Fe(II) complexes reported herein are intensely colored, allowing us to visually monitor reactivity. Ferrous complexes with R = H substituents in the 6-position of the pyridines, [Fe^II(S^Me2N₄(6-H-DPPN)]⁺ (6) and [Fe^II(S^Me2N₄(6-H-DPEN))(MeOH)]⁺ (8-MeOH) are shown to readily bind neutral ligands, and all of the Fe(II) complexes are shown to bind anionic ligands regardless of steric congestion. This reactivity is in contrast to 1 and is attributed to an increased metal ion Lewis acidity assessed via aniodic redox potentials, E_p,a, caused by the π-acid ligands. Thermodynamic parameters (ΔH, ΔS) for neutral ligand binding were obtained from T-dependent equilibrium constants. All but the most sterically congested complex, [Fe^II(S^Me2N₄(6-Me-DPPN)]⁺ (5), react with O₂. In contrast to our Mn(II)-analogues, dioxygen intermediates are not observed. Rates of formation of the final mono oxo-bridged products were assessed via kinetics and shown to be inversely dependent on redox potentials, E_p,a, consistent with a mechanism involving electron transfer.

Calculation of CYP450 protein-ligand binding and dissociation free energy paths

The function of an enzyme depends on its dynamic structure, and the catalytic mechanism has long been an active focus of research. The principle for interpreting protein selectivity and fidelity stems from optimization of the active site upon protein-substrate complexation, i.e., a lock-and-key configuration, on which most protein-substrate molecule binding recognition, and hence drug discovery, relies. Yet another thought has been to incorporate the protein folding interior tunnels for stereo- and regio-selectivity along the protein-substrate or protein-ligand/inhibitor binding process. Free energy calculations provide valuable information for molecular recognition and protein-ligand binding dynamics and kinetics. In this study, we focused on the kinetics of cytochrome P450 proteins (CYP450s) and the protein interior tunnel structure-dynamics relationship in terms of the substrate binding and leaving mechanism. A case in point is given by the prostaglandin H₂ (PGH₂) homologous isomerase of prostacyclin synthase. To calculate the reactant and product traversing the tunnels to and from the heme site, the free energy paths and tunnel potentials of mean force are constructed from steered molecular dynamics simulations and adaptive basing force umbrella sampling simulations. We explore the binding tunnels and critical residue lining characteristics for the ligand traverse and the underlying mechanism of CYP450 activity. Our theoretical analysis provides insights into the decisive role of the substrate tunnel binding process of the CYP450 mechanism and may be useful in drug design and protein engineering contexts.

pH Changes That Induce an Axial Ligand Effect on Nonheme Iron(IV) Oxo Complexes with an Appended Aminopropyl Functionality

Nonheme iron enzymes often utilize a high-valent iron(IV) oxo species for the biosynthesis of natural products, but their high reactivity often precludes structural and functional studies of these complexes. In this work, a combined experimental and computational study is presented on a biomimetic nonheme iron(IV) oxo complex bearing an aminopyridine macrocyclic ligand and its reactivity toward olefin epoxidation upon changes in the identity and coordination ability of the axial ligand. Herein, we show a dramatic effect of the pH on the oxygen-atom-transfer (OAT) reaction with substrates. In particular, these changes have occurred because of protonation of the axial-bound pendant amine group, where its coordination to iron is replaced by a solvent molecule or anionic ligand. This axial ligand effect influences the catalysis, and we observe enhanced cyclooctene epoxidation yields and turnover numbers in the presence of the unbound protonated pendant amine group. Density functional theory studies were performed to support the experiments and highlight that replacement of the pendant amine with a neutral or anionic ligand dramatically lowers the rate-determining barriers of cyclooctene epoxidation. The computational work further establishes that the change in OAT is due to electrostatic interactions of the pendant amine cation that favorably affect the barrier heights.

Electrostatic Perturbations from the Protein Affect C-H Bond Strengths of the Substrate and Enable Negative Catalysis in the TmpA Biosynthesis Enzyme

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C¹ -position of the substrate as part of the biosynthesis of glycine betaine in bacteria and carnitine in humans. To understand how the enzyme avoids breaking the weak C² -H bond in favor of C¹ -hydroxylation, we set up a cluster model of 242 atoms representing the first and second coordination sphere of the metal center and substrate binding pocket, and investigated possible reaction mechanisms of substrate activation by an iron(IV)-oxo species by density functional theory methods. In agreement with experimental product distributions, the calculations predict a favorable C¹ -hydroxylation pathway. The calculations show that the selectivity is guided through electrostatic perturbations inside the protein from charged residues, external electric fields and electric dipole moments. In particular, charged residues influence and perturb the homolytic bond strength of the C¹ -H and C² -H bonds of the substrate, and strongly strengthens the C² -H bond in the substrate-bound orientation.

A comprehensive insight into aldehyde deformylation: mechanistic implications from biology and chemistry

Aldehyde deformylation is an important reaction in biology, organic chemistry and inorganic chemistry and the process has been widely applied and utilized. For instance, in biology, the aldehyde deformylation reaction has wide differences in biological function, whereby cyanobacteria convert aldehydes into alkanes or alkenes, which are used as natural products for, e.g., defense mechanisms. By contrast, the cytochromes P450 catalyse the biosynthesis of hormones, such as estrogen, through an aldehyde deformylation reaction step. In organic chemistry, the aldehyde deformylation reaction is a common process for replacing functional groups on a molecule, and as such, many different synthetic methods and procedures have been reported that involve an aldehyde deformylation step. In bioinorganic chemistry, a variety of metal(iii)-peroxo complexes have been synthesized as biomimetic models and shown to react efficiently with aldehydes through deformylation reactions. This review paper provides an overview of the various aldehyde deformylation reactions in organic chemistry, biology and biomimetic model systems, and shows a broad range of different chemical reaction mechanisms for this process. Although a nucleophilic attack at the carbonyl centre is the consensus reaction mechanism, several examples of an alternative electrophilic reaction mechanism starting with hydrogen atom abstraction have been reported as well. There is still much to learn and to discover on aldehyde deformylation reactions, as deciphered in this review paper.

Histidine versus Cysteine-Bearing Heme-Dependent Halogen Peroxidases: Parallels and Differences for Cl- Oxidation

The homodimeric myeloperoxidase (MPO) features a histidine as a proximal ligand and a sulfonium linkage covalently attaching the heme porphyrin ring to the protein. MPO is able to catalyze Cl^- oxidation with about the same efficiency as chloroperoxidase at pH 7.0. In this study, we seek to explore the parallels and differences between the histidine and cysteine heme-dependent halogen peroxidases. Transition states, reaction barriers, and relevant thermodynamic properties are calculated on protein models. Together with electronic structure calculations, it gives an overview of the reaction mechanisms and of the factors that determine the selectivity between one- and two-electron paths. Conclusions point to the innate oxidizing nature of MPO with the ester and sulfonium linkages hiking up the reactivity to enable chloride oxidation. The installation of a deprotonated imidazolate as a proximal ligand does not shift the equilibrium from one- to two-electron events without influencing the chemistry of the oxidation reaction.

EPR spectroscopy elucidates the electronic structure of [FeV(O)(TAML)] complexes

The complete hyperfine tensor of ¹⁷O of the Fe^V-oxo moeity was probed by ENDOR spectroscopy. The EPR spectroscopic results reported here provide a conclusive experimental basis for elucidating the electronic structure of the Fe^V-oxo complex.

Ascorbate Peroxidase Compound II Is an Iron(IV) Oxo Species

The protonation state of the iron(IV) oxo (or ferryl) form of ascorbate peroxidase compound II (APX-II) is a subject of debate. It has been reported that this intermediate is best described as an iron(IV) hydroxide species. Neutron diffraction data obtained from putative APX-II crystals indicate a protonated oxygenic ligand at 1.88 Å from the heme iron. This finding, if correct, would be unprecedented. A basic iron(IV) oxo species has yet to be spectroscopically observed in a histidine-ligated heme enzyme. The importance of ferryl basicity lies in its connection to our fundamental understanding of C-H bond activation. Basic ferryl species have been proposed to facilitate the oxidation of inert C-H bonds, reactions that are unknown for histidine-ligated hemes enzymes. To provide further insight into the protonation status of APX-II, we examined the intermediate using a combination of Mössbauer and X-ray absorption spectroscopies. Our data indicate that APX-II is an iron(IV) oxo species with an Fe-O bond distance of 1.68 Å, a K-edge pre-edge absorption of 18 units, and Mössbauer parameters of ΔEq = 1.65 mm/s and δ = 0.03 mm/s.

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.

Nature's Machinery, Repurposed: Expanding the Repertoire of Iron-Dependent Oxygenases

Iron is an especially important redox-active cofactor in biology because of its ability to mediate reactions with atmospheric O2. Iron-dependent oxygenases exploit this earth-abundant transition metal for the insertion of oxygen atoms into organic compounds. Throughout the astounding diversity of transformations catalyzed by these enzymes, the protein framework directs reactive intermediates toward the precise formation of products, which, in many cases, necessitates the cleavage of strong C-H bonds. In recent years, members of several iron-dependent oxygenase families have been engineered for new-to-nature transformations that offer advantages over conventional synthetic methods. In this Perspective, we first explore what is known about the reactivity of heme-dependent cytochrome P450 oxygenases and nonheme iron-dependent oxygenases bearing the 2-His-1-carboxylate facial triad by reviewing mechanistic studies with an emphasis on how the protein scaffold maximizes the catalytic potential of the iron-heme and iron cofactors. We then review how these cofactors have been repurposed for abiological transformations by engineering the protein frameworks of these enzymes. Finally, we discuss contemporary challenges associated with engineering these platforms and comment on their roles in biocatalysis moving forward.

Lignin Biodegradation by a Cytochrome P450 Enzyme: A Computational Study into Syringol Activation by GcoA

A recently characterized cytochrome P450 isozyme GcoA activates lignin components through a selective O-demethylation or alternatively an acetal formation reaction. These are important reactions in biotechnology and, because lignin is readily available; it being the main component in plant cell walls. In this work we present a density functional theory study on a large active site model of GcoA to investigate syringol activation by an iron(IV)-oxo heme cation radical oxidant (Compound I) leading to hemiacetal and acetal products. Several substrate-binding positions were tested and full energy landscapes calculated. The study shows that substrate positioning determines the product distributions. Thus, with the phenol group pointing away from the heme, an O-demethylation is predicted, whereas an initial hydrogen-atom abstraction of the weak phenolic O-H group would trigger a pathway leading to ring-closure to form acetal products. Predictions on how to engineer P450 GcoA to get more selective product distributions are given.

Catalytic Mechanism of Aromatic Nitration by Cytochrome P450 TxtE: Involvement of a Ferric-Peroxynitrite Intermediate

The cytochromes P450 are heme-dependent enzymes that catalyze many vital reaction processes in the human body related to biodegradation and biosynthesis. They typically act as mono-oxygenases; however, the recently discovered P450 subfamily TxtE utilizes O₂ and NO to nitrate aromatic substrates such as L-tryptophan. A direct and selective aromatic nitration reaction may be useful in biotechnology for the synthesis of drugs or small molecules. Details of the catalytic mechanism are unknown, and it has been suggested that the reaction should proceed through either an iron(III)-superoxo or an iron(II)-nitrosyl intermediate. To resolve this controversy, we used stopped-flow kinetics to provide evidence for a catalytic cycle where dioxygen binds prior to NO to generate an active iron(III)-peroxynitrite species that is able to nitrate l-Trp efficiently. We show that the rate of binding of O₂ is faster than that of NO and also leads to l-Trp nitration, while little evidence of product formation is observed from the iron(II)-nitrosyl complex. To support the experimental studies, we performed density functional theory studies on large active site cluster models. The studies suggest a mechanism involving an iron(III)-peroxynitrite that splits homolytically to form an iron(IV)-oxo heme (Compound II) and a free NO₂ radical via a small free energy of activation. The latter activates the substrate on the aromatic ring, while compound II picks up the ipso-hydrogen to form the product. The calculations give small reaction barriers for most steps in the catalytic cycle and, therefore, predict fast product formation from the iron(III)-peroxynitrite complex. These findings provide the first detailed insight into the mechanism of nitration by a member of the TxtE subfamily and highlight how the enzyme facilitates this novel reaction chemistry.

Engineering Plant Cytochrome P450s for Enhanced Synthesis of Natural Products: Past Achievements and Future Perspectives

Cytochrome P450s (P450s) are the most versatile catalysts and are widely used by plants to synthesize a vast array of structurally diverse specialized metabolites that not only play essential ecological roles but also constitute a valuable resource for the development of new drugs. To accelerate the metabolic engineering of these high-value metabolites, it is imperative to identify and characterize pathway P450s, and to further improve their activities through protein engineering. In this review, we focus on P450 engineering and summarize the major strategies for enhancing the stability and activity of P450s and successful cases of P450 engineering. Studies in which the functions of P450s were altered to create de novo metabolic pathways or novel compounds are discussed as well. We also overview emerging tools, specifically DNA synthesis, machine learning, and de novo protein design, as well as the evolutionary patterns of P450s unveiled from a massive number of DNA sequences that could be integrated to accelerate the engineering of these enzymes. These approaches would greatly aid in the exploitation of plant-specialized metabolites or derivatives for various uses including medical applications.

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.

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

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

How Metal Ion Lewis Acidity and Steric Properties Influence the Barrier to Dioxygen Binding, Peroxo O-O Bond Cleavage, and Reactivity

Herein we quantitatively investigate how metal ion Lewis acidity and steric properties influence the kinetics and thermodynamics of dioxygen binding versus release from structurally analogous Mn-O2 complexes, as well as the barrier to Mn peroxo O-O bond cleavage, and the reactivity of Mn oxo intermediates. Previously we demonstrated that the steric and electronic properties of MnIII-OOR complexes containing N-heterocyclic (NAr) ligand scaffolds can have a dramatic influence on alkylperoxo O-O bond lengths and the barrier to alkylperoxo O-O bond cleavage. Herein, we examine the dioxygen reactivity of a new MnII complex containing a more electron-rich, less sterically demanding NAr ligand scaffold, and compare it with previously reported MnII complexes. Dioxygen binding is shown to be reversible with complexes containing the more electron-rich metal ions. The kinetic barrier to O2 binding and peroxo O-O bond cleavage is shown to correlate with redox potentials, as well as the steric properties of the supporting NAr ligands. The reaction landscape for the dioxygen chemistry of the more electron-rich complexes is shown to be relatively flat. A total of four intermediates, including a superoxo and peroxo species, are observed with the most electron-rich complex. Two new intermediates are shown to form following the peroxo, which are capable of cleaving strong X-H bonds. In the absence of a sacrificial H atom donor, solvent, or ligand, serves as a source of H atoms. With TEMPOH as sacrificial H atom donor, a deuterium isotope effect is observed (kH/kD = 3.5), implicating a hydrogen atom transfer (HAT) mechanism. With 1,4-cyclohexadiene, 0.5 equiv of benzene is produced prior to the formation of an EPR detected MnIIIMnIV bimetallic species, and 0.5 equiv after its formation.

Regulating the Basicity of Metal-Oxido Complexes with a Single Hydrogen Bond and Its Effect on C-H Bond Cleavage

The functionalization of C-H bonds is an essential reaction in biology and chemistry. Metalloenzymes that often exhibit this type of reactivity contain metal-oxido intermediates that are directly involved in the initial cleavage of the C-H bonds. Regulation of the cleavage process is achieved, in part, by hydrogen bonds that are proximal to the metal-oxido units, yet our understanding of their exact role(s) is still emerging. To gain further information into the role of H-bonds on C-H bond activation, a hybrid set of urea-containing tripodal ligands has been developed in which a single H-bond can be adjusted through changes in the properties of one ureayl N-H bond. This modularity is achieved by appending a phenyl ring with different para-substituents from one ureayl NH group. The ligands have been used to prepare a series of Mn^III-oxido complexes, and a Hammett correlation was found between the pK_a values of the complexes and the substituents on the phenyl ring that was explained within the context of changes to the H-bonds involving the Mn^III-oxido unit. The complexes were tested for their reactivity toward 9,10-dihydroanthracene (DHA), and a Hammett correlation was found between the second-order rate constants for the reactions and the pK_a values. Studies to determine activation parameters and the kinetic isotope effects are consistent with a mechanism in which rate-limiting proton transfer is an important contributor. However, additional reactivity studies with xanthene found a significant increase in the rate constant compared to DHA, even though the substrates have the same pK_a(C-H) values. These results do not support a discrete proton-transfer/electron-transfer process, but rather an asynchronous mechanism in which the proton and electron are transferred unequally at the transition state.

Site-Selective Supramolecular Complexation Activates Catalytic Ethane Oxidation by a Nitrido-Bridged Iron Porphyrinoid Dimer

Development of supramolecular methods to further activate a highly reactive intermediate is a fascinating strategy to create novel potent catalysts for activation of inert chemicals. Herein, a supramolecular approach to enhance the oxidizing ability of a high-valent oxo species of a nitrido-bridged iron porphyrinoid dimer that is a known potent molecular catalyst for light alkane oxidation is reported. For this purpose, a nitrido-bridged dinuclear iron complex of porphyrin-phthalocyanine heterodimer 3⁵⁺ , which is connected through a fourfold rotaxane, was prepared. Heterodimer 3⁵⁺ catalyzed ethane oxidation in the presence of H₂ O₂ at a relatively low temperature. The site-selective complexation of 3⁵⁺ with an additional anionic porphyrin (TPPS^4- ) through π-π stacking and electrostatic interactions afforded a stable 1:1 complex. It was demonstrated that the supramolecular post-synthetic modification of 3⁵⁺ enhances its catalytic activity efficiently. Moreover, supramolecular conjugates achieved higher catalytic ethane oxidation activity than nitrido-bridged iron phthalocyanine dimer, which is the most potent iron-oxo-based molecular catalyst for light-alkane oxidation reported so far. Electrochemical measurements proved that the electronic perturbation from TPPS^4- to 3⁵⁺ enhanced the catalytic activity.