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

An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models (A, B, and C) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C4-H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C, the dipole moment is positioned along the C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr233 phenol loses a proton to the nearby Asp179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

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.

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

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

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways.

Using BpyAla to generate copper artificial metalloenzymes: a catalytic and structural study

Artificial metalloenzymes (ArMs) have emerged as a promising avenue in the field of biocatalysis, offering new reactivity. However, their design remains challenging due to the limited understanding of their protein dynamics and how the introduced cofactors alter the protein scaffold structure. Here we present the structures and catalytic activity of novel copper ArMs capable of (R)- or (S)-stereoselective control, utilizing a steroid carrier protein (SCP) scaffold. To incorporate 2,2'-bipyridine (Bpy) into SCP, two distinct strategies were employed: either Bpy was introduced as an unnatural amino acid (2,2'-bipyridin-5-yl)alanine (BpyAla) using amber stop codon expression or via bioconjugation of bromomethyl-Bpy to cysteine residues. The resulting ArMs proved to be effective at catalysing an enantioselective Friedel-Crafts reaction with SCP_Q111BpyAla achieving the best selectivity with an enantioselectivity of 72% ee (S). Interestingly, despite using the same protein scaffold, different attachment strategies for Bpy at the same residue (Q111) led to a switch in the enantiopreference of the ArM. X-ray crystal structures of SCP_Q111CBpy and SCP_Q111BpyAla ArMs with bound Cu(ii) ions unveiled crucial differences in the orientation of the catalytic centre. Combining structural information, alanine scanning studies, and computational analysis shed light on the distinct active sites of the ArMs, clarifying that these active sites stabilise the nucleophilic substrate on different sides of the electrophile leading to the observed switch in enantioselectivity. This work underscores the importance of integrating structural studies with catalytic screening to unravel the intricacies of ArM behaviour and facilitate their development for targeted applications in biocatalysis.

CO2 Reduction by an Iron(I) Porphyrinate System: Effect of Hydrogen Bonding on the Second Coordination Sphere

Transforming CO2 into valuable materials is an important reaction in catalysis, especially because CO2 concentrations in the atmosphere have been growing steadily due to extensive fossil fuel usage. From an environmental perspective, reduction of CO2 to valuable materials should be catalyzed by an environmentally benign catalyst and avoid the use of heavy transition-metal ions. In this work, we present a computational study into a novel iron(I) porphyrin catalyst for CO2 reduction, namely, with a tetraphenylporphyrin ligand and analogues. In particular, we investigated iron(I) tetraphenylporphyrin with one of the meso-phenyl groups substituted with o-urea, p-urea, or o-2-amide groups. These substituents can provide hydrogen-bonding interactions in the second coordination sphere with bound ligands and assist with proton relay. Furthermore, our studies investigated bicarbonate and phenol as stabilizers and proton donors in the reaction mechanism. Potential energy landscapes for double protonation of iron(I) porphyrinate with bound CO2 are reported. The work shows that the bicarbonate bridges the urea/amide groups to the CO2 and iron center and provides a tight bonding pattern with strong hydrogen-bonding interactions that facilitates easy proton delivery and reduction of CO2. Specifically, bicarbonate provides a low-energy proton shuttle mechanism to form CO and water efficiently. Furthermore, the o-urea group locks bicarbonate and CO2 in a tight orientation and helps with ideal proton transfer, while there is more mobility and lesser stability with an o-amide group in that position instead. Our calculations show that the o-urea group leads to reduction in proton-transfer barriers, in line with experimental observation. We then applied electric-field-effect calculations to estimate the environmental effects on the two proton-transfer steps in the reaction. These calculations describe the perturbations that enhance the driving forces for the proton-transfer steps and have been used to make predictions about how the catalysts can be further engineered for more enhanced CO2 reduction processes.

Biocatalytic role of cytochrome P450s to produce antibiotics: A review

Cytochrome P450s belong to a family of heme-binding monooxygenases, which catalyze regio- and stereospecific functionalisation of C-H, C-C, and C-N bonds, including heteroatom oxidation, oxidative C-C bond cleavages, and nitrene transfer. P450s are considered useful biocatalysts for the production of pharmaceutical products, fine chemicals, and bioremediating agents. Despite having tremendous biotechnological potential, being heme-monooxygenases, P450s require either autologous or heterologous redox partner(s) to perform chemical transformations. Randomly distributed P450s throughout a bacterial genome and devoid of particular redox partners in natural products biosynthetic gene clusters (BGCs) showed an extra challenge to reveal their pharmaceutical potential. However, continuous efforts have been made to understand their involvement in antibiotic biosynthesis and their modification, and this review focused on such BGCs. Here, particularly, we have discussed the role of P450s involved in the production of macrolides and aminocoumarin antibiotics, nonribosomal peptide (NRPSs) antibiotics, ribosomally synthesized and post-translationally modified peptide (RiPPs) antibiotics, and others. Several reactions catalyzed by P450s, as well as the role of their redox partners involved in the BGCs of various antibiotics and their derivatives, have been primarily addressed in this review, which would be useful in further exploration of P450s for the biosynthesis of new therapeutics.

Mechanism of CO2 Reduction to Methanol with H2 on an Iron(II)-scorpionate Catalyst

CO2 utilization is an important process in the chemical industry with great environmental power. In this work we show how CO2 and H2 can be reacted to form methanol on an iron(II) center and highlight the bottlenecks for the reaction and what structural features of the catalyst are essential for efficient turnover. The calculations predict the reactions to proceed through three successive reaction cycles that start with heterolytic cleavage of H2 followed by sequential hydride and proton transfer processes. The H2 splitting process is an endergonic process and hence high pressures will be needed to overcome this step and trigger the hydrogenation reaction. Moreover, H2 cleavage into a hydride and proton requires a metal to bind hydride and a nearby source to bind the proton, such as an amide or pyrazolyl group, which the scorpionate ligand used here facilitates. As such the computations highlight the non-innocence of the ligand scaffold through proton shuttle from H2 to substrate as an important step in the reaction mechanism.

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.

Direct aromatic nitration by bacterial P450 enzymes

Nitro aromatics have broad applications in industry, agriculture, and pharmaceutics. However, their industrial production is faced with many challenges including poor selectivity, heavy pollution and safety concerns. Nature provides multiple strategies for aromatic nitration, which opens the door for the development of green and efficient biocatalysts. Our group's efforts focused on a unique bacterial cytochrome P450 TxtE that originates from the biosynthetic pathway of phytotoxin thaxtomins, which can install a nitro group at C4 of l-Trp indole ring. TxtE is a Class I P450 and its reaction relies on a pair of redox partners ferredoxin and ferredoxin reductase for essential electron transfer. To develop TxtE as an efficient nitration biocatalyst, we created artificial self-sufficient P450 chimeras by fusing TxtE with the reductase domain of the bacterial P450BM3 (BM3R). We evaluated the catalytic performance of the chimeras with different lengths of the linker connecting TxtE and BM3R domains and identified one with a 14-amino-acid linker (TB14) to give the best activity. In addition, we demonstrated the broad substrate scope of the engineered biocatalyst by screening diverse l-Trp analogs. In this chapter, we provide a detailed procedure for the development of aromatic nitration biocatalysts, including the construction of P450 fusion chimeras, biochemical characterization, determination of catalytic parameters, and testing of enzyme-substrate scope. These protocols can be followed to engineer other P450 enzymes and illustrate the processes of biocatalytic development for the synthesis of nitro chemicals.

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.

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

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

Insights into Substrate Recognition by the Unusual Nitrating Enzyme RufO

Nitration reactions are crucial for many industrial syntheses; however, current protocols lack site specificity and employ hazardous chemicals. The noncanonical cytochrome P450 enzymes RufO and TxtE catalyze the only known direct aromatic nitration reactions in nature, making them attractive model systems for the development of analogous biocatalytic and/or biomimetic reactions that proceed under mild conditions. While the associated mechanism has been well-characterized in TxtE, much less is known about RufO. Herein we present the first structure of RufO alongside a series of computational and biochemical studies investigating its unusual reactivity. We demonstrate that free l-tyrosine is not readily accepted as a substrate despite previous reports to the contrary. Instead, we propose that RufO natively modifies l-tyrosine tethered to the peptidyl carrier protein of a nonribosomal peptide synthetase encoded by the same biosynthetic gene cluster and present both docking and molecular dynamics simulations consistent with this hypothesis. Our results expand the scope of direct enzymatic nitration reactions and provide the first evidence for such a modification of a peptide synthetase-bound substrate. Both of these insights may aid in the downstream development of biocatalytic approaches to synthesize rufomycin analogues and related drug candidates.

Structural and spectroscopic characterization of RufO indicates a new biological role in rufomycin biosynthesis

Rufomycins constitute a class of cyclic heptapeptides isolated from actinomycetes. They are secondary metabolites that show promising treatment against Mycobacterium tuberculosis infections by inhibiting a novel drug target. Several nonproteinogenic amino acids are integrated into rufomycins, including a conserved 3-nitro-tyrosine. RufO, a cytochrome P450 (CYP)-like enzyme, was proposed to catalyze the formation of 3-nitro-tyrosine in the presence of O2 and NO. To define its biological function, the interaction between RufO and the proposed substrate tyrosine is investigated using various spectroscopic methods that are sensitive to the structural change of a heme center. However, a low- to high-spin state transition and a dramatic increase in the redox potential that are commonly found in CYPs upon ligand binding have not been observed. Furthermore, a 1.89-Å crystal structure of RufO shows that the enzyme has flexible surface regions, a wide-open substrate access tunnel, and the heme center is largely exposed to solvent. Comparison with a closely related nitrating CYP reveals a spacious and hydrophobic distal pocket in RufO, which is incapable of stabilizing a free amino acid. Molecular docking validates the experimental data and proposes a possible substrate. Collectively, our results disfavor tyrosine as the substrate of RufO and point to the possibility that the nitration occurs during or after the assembly of the peptides. This study indicates a new function of the unique nitrating enzyme and provides insights into the biosynthesis of nonribosomal peptides.

How Is Substrate Halogenation Triggered by the Vanadium Haloperoxidase from Curvularia inaequalis?

Vanadium haloperoxidases (VHPOs) are unique enzymes in biology that catalyze a challenging halogen transfer reaction and convert a strong aromatic C-H bond into C-X (X = Cl, Br, I) with the use of a vanadium cofactor and H₂O₂. The VHPO catalytic cycle starts with the conversion of hydrogen peroxide and halide (X = Cl, Br, I) into hypohalide on the vanadate cofactor, and the hypohalide subsequently reacts with a substrate. However, it is unclear whether the hypohalide is released from the enzyme or otherwise trapped within the enzyme structure for the halogenation of organic substrates. A substrate-binding pocket has never been identified for the VHPO enzyme, which questions the role of the protein in the overall reaction mechanism. Probing its role in the halogenation of small molecules will enable further engineering of the enzyme and expand its substrate scope and selectivity further for use in biotechnological applications as an environmentally benign alternative to current organic chemistry synthesis. Using a combined experimental and computational approach, we elucidate the role of the vanadium haloperoxidase protein in substrate halogenation. Activity studies show that binding of the substrate to the enzyme is essential for the reaction of the hypohalide with substrate. Stopped-flow measurements demonstrate that the rate-determining step is not dependent on substrate binding but partially on hypohalide formation. Using a combination of molecular mechanics (MM) and molecular dynamics (MD) simulations, the substrate binding area in the protein is identified and even though the selected substrates (methylphenylindole and 2-phenylindole) have limited hydrogen-bonding abilities, they are found to bind relatively strongly and remain stable in a binding tunnel. A subsequent analysis of the MD snapshots characterizes two small tunnels leading from the vanadate active site to the surface that could fit small molecules such as hypohalide, halide, and hydrogen peroxide. Density functional theory studies using electric field effects show that a polarized environment in a specific direction can substantially lower barriers for halogen transfer. A further analysis of the protein structure indeed shows a large dipole orientation in the substrate-binding pocket that could enable halogen transfer through an applied local electric field. These findings highlight the importance of the enzyme in catalyzing substrate halogenation by providing an optimal environment to lower the energy barrier for this challenging aromatic halide insertion reaction.

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

The Engagement of Cytochrome P450 Enzymes in Tryptophan Metabolism

Tryptophan is metabolized along three main metabolic pathways, namely the kynurenine, serotonin and indole pathways. The majority of tryptophan is transformed via the kynurenine pathway, catalyzed by tryptophan-2,3-dioxygenase or indoleamine-2,3-dioxygenase, leading to neuroprotective kynurenic acid or neurotoxic quinolinic acid. Serotonin synthesized by tryptophan hydroxylase, and aromatic L-amino acid decarboxylase enters the metabolic cycle: serotonin → N-acetylserotonin → melatonin → 5-methoxytryptamine→serotonin. Recent studies indicate that serotonin can also be synthesized by cytochrome P450 (CYP), via the CYP2D6-mediated 5-methoxytryptamine O-demethylation, while melatonin is catabolized by CYP1A2, CYP1A1 and CYP1B1 via aromatic 6-hydroxylation and by CYP2C19 and CYP1A2 via O-demethylation. In gut microbes, tryptophan is metabolized to indole and indole derivatives. Some of those metabolites act as activators or inhibitors of the aryl hydrocarbon receptor, thus regulating the expression of CYP1 family enzymes, xenobiotic metabolism and tumorigenesis. The indole formed in this way is further oxidized to indoxyl and indigoid pigments by CYP2A6, CYP2C19 and CYP2E1. The products of gut-microbial tryptophan metabolism can also inhibit the steroid-hormone-synthesizing CYP11A1. In plants, CYP79B2 and CYP79B3 were found to catalyze N-hydroxylation of tryptophan to form indole-3-acetaldoxime while CYP83B1 was reported to form indole-3-acetaldoxime N-oxide in the biosynthetic pathway of indole glucosinolates, considered to be defense compounds and intermediates in the biosynthesis of phytohormones. Thus, cytochrome P450 is engaged in the metabolism of tryptophan and its indole derivatives in humans, animals, plants and microbes, producing biologically active metabolites which exert positive or negative actions on living organisms. Some tryptophan-derived metabolites may influence cytochrome P450 expression, affecting cellular homeostasis and xenobiotic metabolism.

Engineering C-C Bond Cleavage Activity into a P450 Monooxygenase Enzyme

The cytochrome P450 (CYP) superfamily of heme monooxygenases has demonstrated ability to facilitate hydroxylation, desaturation, sulfoxidation, epoxidation, heteroatom dealkylation, and carbon-carbon bond formation and cleavage (lyase) reactions. Seeking to study the carbon-carbon cleavage reaction of α-hydroxy ketones in mechanistic detail using a microbial P450, we synthesized α-hydroxy ketone probes based on the physiological substrate for a well-characterized benzoic acid metabolizing P450, CYP199A4. After observing low activity with wild-type CYP199A4, subsequent assays with an F182L mutant demonstrated enzyme-dependent C-C bond cleavage toward one of the α-hydroxy ketones. This C-C cleavage reaction was subject to an inverse kinetic solvent isotope effect analogous to that observed in the lyase activity of the human P450 CYP17A1, suggesting the involvement of a species earlier than Compound I in the catalytic cycle. Co-crystallization of F182L-CYP199A4 with this α-hydroxy ketone showed that the substrate bound in the active site with a preference for the (S)-enantiomer in a position which could mimic the topology of the lyase reaction in CYP17A1. Molecular dynamics simulations with an oxy-ferrous model of CYP199A4 revealed a displacement of the substrate to allow for oxygen binding and the formation of the lyase transition state proposed for CYP17A1. This demonstration that a correctly positioned α-hydroxy ketone substrate can realize lyase activity with an unusual inverse solvent isotope effect in an engineered microbial system opens the door for further detailed biophysical and structural characterization of CYP catalytic intermediates.

Nitration of Pyrrolo[2,1-a]isoquinolines

We have successfully modified a series of pyrrolo[2,1-a]isoquinolines via direct nitration under mild reaction conditions. Easily accessible nitrates including CAN, Cu(NO₃)₂·H₂O, and Fe(NO₃)₃·9H₂O all can serve as effective nitrating reagents for functionalizing pyrrolo[2,1-a]isoquinolines. Various nitro-bearing pyrrolo[2,1-a]isoquinolines have been efficiently prepared in acceptable to good yields.

QM Calculations Revealed that Outer-Sphere Electron Transfer Boosted O-O Bond Cleavage in the Multiheme-Dependent Cytochrome bd Oxygen Reductase

The cytochrome bd oxygen reductase catalyzes the four-electron reduction of dioxygen to two water molecules. The structure of this enzyme reveals three heme molecules in the active site, which differs from that of heme-copper cytochrome c oxidase. The quantum chemical cluster approach was used to uncover the reaction mechanism of this intriguing metalloenzyme. The calculations suggested that a proton-coupled electron transfer reduction occurs first to generate a ferrous heme b₅₉₅. This is followed by the dioxygen binding at the heme d center coupled with an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the dioxygen moiety, affording a ferric ion superoxide intermediate. A second proton-coupled electron transfer produces a heme d ferric hydroperoxide, which undergoes efficient O-O bond cleavage facilitated by an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the O-O σ* orbital and an inner-sphere proton transfer from the heme d hydroxyl group to the leaving hydroxide. The synergistic benefits of the two types of hemes rationalize the highly efficient oxygen reduction repertoire for the multi-heme-dependent cytochrome bd oxygen reductase family.

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.

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.

An iron(ii)-based metalloradical system for intramolecular amination of C(sp2)-H and C(sp3)-H bonds: synthetic applications and mechanistic studies

A catalytic system for intramolecular C(sp2)-H and C(sp3)-H amination of substituted tetrazolopyridines has been successfully developed. The amination reactions are developed using an iron-porphyrin based catalytic system. It has been demonstrated that the same iron-porphyrin based catalytic system efficiently activates both the C(sp2)-H and C(sp3)-H bonds of the tetrazole as well as azide-featuring substrates with a high level of regioselectivity. The method exhibited an excellent functional group tolerance. The method affords three different classes of high-value N-heterocyclic scaffolds. A number of important late-stage C-H aminations have been performed to access important classes of molecules. Detailed studies (experimental and computational) showed that both the C(sp2)-H and C(sp3)-H amination reactions involve a metalloradical activation mechanism, which is different from the previously reported electro-cyclization mechanism. Collectively, this study reports the discovery of a new class of metalloradical activation modes using a base metal catalyst that should find wide application in the context of medicinal chemistry, drug discovery and industrial applications.

The Mechanism of Biochemical NO‐Sensing: Insights from Computational Chemistry

The binding of small gas molecules such as NO and CO plays a major role in the signaling routes of the human body. The sole NO‐receptor in humans is soluble guanylyl cyclase (sGC) – a histidine‐ligated heme protein, which, upon NO binding, activates a downstream signaling cascade. Impairment of NO‐signaling is linked, among others, to cardiovascular and inflammatory diseases. In the present work, we use a combination of theoretical tools such as MD simulations, high‐level quantum chemical calculations and hybrid QM/MM methods to address various aspects of NO binding and to elucidate the most likely reaction paths and the potential intermediates of the reaction. As a model system, the H‐NOX protein from Shewanella oneidensis (So H‐NOX) homologous to the NO‐binding domain of sGC is used. The signaling route is predicted to involve NO binding to form a six‐coordinate intermediate heme‐NO complex, followed by relatively facile His decoordination yielding a five‐coordinate adduct with NO on the distal side with possible isomerization to the proximal side through binding of a second NO and release of the first one. MD simulations show that the His sidechain can quite easily rotate outward into solvent, with this motion being accompanied in our simulations by shifts in helix positions that are consistent with this decoordination leading to significant conformational change in the protein.

Non-Ionic Inulin-Based Polymer Nanodiscs Enable Functional Reconstitution of a Redox Complex Composed of Oppositely Charged CYP450 and CPR in a Lipid Bilayer Membrane

Although polymer-based lipid nanodiscs are increasingly used in the structural studies of membrane proteins, the charge of the belt-forming polymer is a major limitation for functional reconstitution of membrane proteins possessing an opposite net charge to that of the polymer. This limitation also rules out the reconstitution of a protein-protein complex composed of oppositely charged membrane proteins. In this study, we report the first successful functional reconstitution of a membrane-bound redox complex constituting a cationic cytochrome P450 (CYP450) and an anionic cytochrome P450 reductase (CPR) in non-ionic inulin-based lipid nanodiscs. The gel-to-liquid-crystalline phase-transition temperature (Tm) of DMPC:DMPG (7:3 w/w) lipids in polymer nanodiscs was determined by differential scanning calorimetry (DSC) and 31P NMR experiments. The CYP450-CPR redox complex reconstitution in polymer nanodiscs was characterized by size-exclusion chromatography (SEC), and the electron transfer kinetics was carried out using the stopped-flow technique under anaerobic conditions. The Tm of DMPC:DMPG (7:3 w/w) in polymer nanodiscs measured from 31P NMR agrees with that obtained from DSC and was found to be higher than that for liposomes due to the decreased cooperativity of lipids present in the nanodiscs. The stopped-flow measurements revealed the CYP450-CPR redox complex reconstituted in nanodiscs to be functional, and the electron transfer kinetics was found to be temperature-dependent. Based on the successful demonstration of the use of non-ionic inulin-based polymer nanodiscs reported in this study, we expect them to be useful in studying the function and structures of a variety of membrane proteins/complexes irrespective of the charge of the molecular components. Since the polymer nanodiscs were shown to align in an externally applied magnetic field, they can also be used to measure residual dipolar couplings (RDCs) and residual quadrupolar couplings (RQCs) for various molecules ranging from small molecules to soluble proteins and nucleic acids.

Redesigning Enzymes for Biocatalysis: Exploiting Structural Understanding for Improved Selectivity

The discovery of new enzymes, alongside the push to make chemical processes more sustainable, has resulted in increased industrial interest in the use of biocatalytic processes to produce high-value and chiral precursor chemicals. Huge strides in protein engineering methodology and in silico tools have facilitated significant progress in the discovery and production of enzymes for biocatalytic processes. However, there are significant gaps in our knowledge of the relationship between enzyme structure and function. This has demonstrated the need for improved computational methods to model mechanisms and understand structure dynamics. Here, we explore efforts to rationally modify enzymes toward changing aspects of their catalyzed chemistry. We highlight examples of enzymes where links between enzyme function and structure have been made, thus enabling rational changes to the enzyme structure to give predictable chemical outcomes. We look at future directions the field could take and the technologies that will enable it.

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.

Electrostatic perturbations in the substrate-binding pocket of taurine/α-ketoglutarate dioxygenase determine its selectivity

Taurine/α‐ketoglutarate dioxygenase is an important enzyme that takes part in the cysteine catabolism process in the human body and selectively hydroxylates taurine at the C¹‐position. Recent computational studies showed that in the gas‐phase the C²−H bond of taurine is substantially weaker than the C¹−H bond, yet no evidence exists of 2‐hydroxytaurine products. To this end, a detailed computational study on the selectivity patterns in TauD was performed. The calculations show that the second‐coordination sphere and the protonation states of residues play a major role in guiding the enzyme to the right selectivity. Specifically, a single proton on an active site histidine residue can change the regioselectivity of the reaction through its electrostatic perturbations in the active site and effectively changes the C¹−H and C²−H bond strengths of taurine. This is further emphasized by many polar and hydrogen bonding interactions of the protein cage in TauD with the substrate and the oxidant that weaken the pro‐R C¹−H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently.

Biocatalysis making waves in organic chemistry

Biocatalysis has an enormous impact on chemical synthesis. The waves in which biocatalysis has developed, and in doing so changed our perception of what organic chemistry is, were reviewed 20 and 10 years ago. Here we review the consequences of these waves of development. Nowadays, hydrolases are widely used on an industrial scale for the benign synthesis of commodity and bulk chemicals and are fully developed. In addition, further enzyme classes are gaining ever increasing interest. Particularly, enzymes catalysing selective C-C-bond formation reactions and enzymes catalysing selective oxidation and reduction reactions are solving long-standing synthetic challenges in organic chemistry. Combined efforts from molecular biology, systems biology, organic chemistry and chemical engineering will establish a whole new toolbox for chemistry. Recent developments are critically reviewed.

Biodegradation of herbicides by a plant nonheme iron dioxygenase: mechanism and selectivity of substrate analogues

Aryloxyalkanoate dioxygenases are unique herbicide biodegrading nonheme iron enzymes found in plants and hence, from environmental and agricultural point of view they are important and valuable. However, they often are substrate specific and little is known on the details of the mechanism and the substrate scope. To this end, we created enzyme models and calculate the mechanism for 2,4-dichlorophenoxyacetic acid biodegradation and 2-methyl substituted analogs by density functional theory. The work shows that the substrate binding is tight and positions the aliphatic group close to the metal center to enable a chemoselective reaction mechanism to form the C 2 -hydroxy products, whereas the aromatic hydroxylation barriers are well higher in energy. Subsequently, we investigated the metabolism of R - and S -methyl substituted inhibitors and show that these do not react as efficiently as 2,4-dichlorophenoxyacetic acid substrate due to stereochemical clashes in the active site and particularly for the R -isomer give high rebound barriers.

Bio-inspired nitrogen oxide (NOx) interconversion reactivities of synthetic heme Compound-I and Compound-II intermediates

Dioxygen activating heme enzymes have long predicted to be powerhouses for nitrogen oxide interconversion, especially for nitric oxide (NO) oxidation which has far-reaching biological and/or environmental impacts. Lending credence, reactivity of NO with high-valent heme‑oxygen intermediates of globin proteins has recently been implicated in the regulation of a variety of pivotal physiological events such as modulating catalytic activities of various heme enzymes, enhancing antioxidant activity to inhibit oxidative damage, controlling inflammatory and infectious properties within the local heme environments, and NO scavenging. To reveal insights into such crucial biological processes, we have investigated low temperature NO reactivities of two classes of synthetic high-valent heme intermediates, Compound-II and Compound-I. In that, Compound-II rapidly reacts with NO yielding the six-coordinate (NO bound) heme ferric nitrite complex, which upon warming to room temperature converts into the five-coordinate heme ferric nitrite species. These ferric nitrite complexes mediate efficient substrate oxidation reactions liberating NO; i.e., shuttling NO₂^- back to NO. In contrast, Compound-I and NO proceed through an oxygen-atom transfer process generating the strong nitrating agent NO₂, along with the corresponding ferric nitrosyl species that converts to the naked heme ferric parent complex upon warmup. All reaction components have been fully characterized by UV-vis, ²H NMR and EPR spectroscopic methods, mass spectrometry, elemental analyses, and semi-quantitative determination of NO₂^- anions. The clean, efficient, potentially catalytic NO_x interconversions driven by high-valent heme species presented herein illustrate the strong prospects of a heme enzyme/O₂/NO_x dependent unexplored territory that is central to human physiology, pathology, and therapeutics.

Cometabolism of 17α-ethynylestradiol by nitrifying bacteria depends on reducing power availability and leads to elevated nitric oxide formation

17α-ethynylestradiol (EE2) is a priority emerging contaminant (EC) in diverse environments that can be cometabolized by ammonia oxidizing bacteria (AOB). However, its transformation kinetics and the underlying molecular mechanism are unclear. In this study, kinetic parameters, including maximum specific EE2 transformation rate, EE2 half-saturation coefficient, and EE2transformation capacity of AOBwere obtained by using the model AOB strain, Nitrosomonas europaea 19718. The relationship between EE2 cometabolism and ammonia oxidation was divided into three phases according to reducing power availability, namely "activation", "coupling", and "saturation". Specifically, there was a universal lag of EE2 transformation after ammonia oxidation was initiated, suggesting that sufficient reducing power (approximately 0.95 ± 0.06 mol NADH/L) was required to activate EE2 cometabolism. Interestingly, nitric oxide emission increased by 12 ± 2% during EE2 cometabolism, along with significantly upregulated nirK cluster genes. The findings are of importance to understanding the cometabolic behavior and mechanism of EE2 in natural and engineered environments. Maintaining relatively high and stable reducing power supply from ammonia oxidation can potentially improve the cometabolic removal of EE2 and other ECs during wastewater nitrification processes.

The Ferric-Superoxo Intermediate of the TxtE Nitration Pathway Resists Reduction, Facilitating Its Reaction with Nitric Oxide

TxtE is a cytochrome P450 (CYP) homologue that mediates the nitric oxide (NO)-dependent direct nitration of l-tryptophan (Trp) to form 4-nitro-l-tryptophan (4-NO₂-Trp). A recent report showed evidence that TxtE activity requires NO to react with a ferric-superoxo intermediate. Given this minimal mechanism, it is not clear how TxtE avoids Trp hydroxylation, a mechanism that also traverses the ferric-superoxo intermediate. To provide insight into canonical CYP intermediates that TxtE can access, electron coupling efficiencies to form 4-NO₂-Trp under single- or limited-turnover conditions were measured and compared to steady-state efficiencies. As previously reported, Trp nitration by TxtE is supported by the engineered self-sufficient variant, TB14, as well as by reduced putidaredoxin. Ferrous (Fe^II) TxtE exhibits excellent electron coupling (70%), which is 50-fold higher than that observed under turnover conditions. In addition, two- or four-electron reduced TB14 exhibits electron coupling (∼6%) that is 2-fold higher than that of one-electron reduced TB14 (3%). The combined results suggest (1) autoxidation is the sole TxtE uncoupling pathway and (2) the TxtE ferric-superoxo intermediate cannot be reduced by these electron transfer partners. The latter conclusion is further supported by ultraviolet-visible absorption spectral time courses showing neither spectral nor kinetic evidence for reduction of the ferric-superoxo intermediate. We conclude that resistance of the ferric-superoxo intermediate to reduction is a key feature of TxtE that increases the lifetime of the intermediate and enables its reaction with NO and efficient nitration activity.

Theoretical Investigation on the Elusive Reaction Mechanism of Spirooxindole Formation Mediated by Cytochrome P450s: A Nascent Feasible Charge-Shift C-O Bond Makes a Difference

Spirooxindoles are pivotal biofunctional groups widely distributed in natural products and clinic drugs. However, construction of such subtle chiral skeletons is a long-standing challenge to both organic and bioengineering scientists. The knowledge of enzymatic spirooxindole formation in nature may inspire rational design of new catalysts. To this end, we presented a theoretical investigation on the elusive mechanism of the spiro-ring formation at the 3-position of oxindole mediated by cytochrome P450 enzymes (P450). Our calculated results demonstrated that the electrophilic attack of CpdI, the active species of P450, to the substrate, shows regioselectivity, i.e., the attack at the C9 position forms a tetrahedral intermediate involving an unusual feasible charge-shift C9^δ+-O^δ- bond, while the attack at the C1 position forms an epoxide intermediate. The predominant route is the first route with the charge-shift bonding intermediate due to holding a relatively lower barrier by >5 kcal mol^-1 than the epoxide route, which fits the experimental observations. Such a delocalized charge-shift bond facilitates the formation of a spiro-ring mainly through elongation of the C1-C9 bond to eliminate the aromatization of the tricyclic beta-carboline. Our theoretical results shed profound mechanistic insights for the first time into the elusive spirooxindole formation mediated by P450s.

Structure and Functional Differences of Cysteine and 3-Mercaptopropionate Dioxygenases: A Computational Study

Thiol dioxygenases are important enzymes for human health; they are involved in the detoxification and catabolism of toxic thiol-containing natural products such as cysteine. As such, these enzymes have relevance to the development of Alzheimer's and Parkinson's diseases in the brain. Recent crystal structure coordinates of cysteine and 3-mercaptopropionate dioxygenase (CDO and MDO) showed major differences in the second-coordination spheres of the two enzymes. To understand the difference in activity between these two analogous enzymes, we created large, active-site cluster models. We show that CDO and MDO have different iron(III)-superoxo-bound structures due to differences in ligand coordination. Furthermore, our studies show that the differences in the second-coordination sphere and particularly the position of a positively charged Arg residue results in changes in substrate positioning, mobility and enzymatic turnover. Furthermore, the substrate scope of MDO is explored with cysteinate and 2-mercaptosuccinic acid and their reactivity is predicted.

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.

Expanding the Substrate Scope of Nitrating Cytochrome P450 TxtE by Active Site Engineering of a Reductase Fusion

Aromatic nitration reactions are a cornerstone of organic chemistry, but are challenging to scale due to corrosive reagents and elevated temperatures. The cytochrome P450 TxtE nitrates the indole 4-position of l-tryptophan at room temperature using NO, O2 and NADPH, and has potential to be developed into a useful aromatic nitration biocatalyst. However, its narrow substrate scope (requiring both the α-amino acid and indole functionalities) have hindered this. Screening of an R59 mutant library of a TxtE-reductase fusion protein identified a variant (R59C) that nitrates tryptamine, which is not accepted by native TxtE. This variant exhibits a broader substrate scope than the wild type enzyme and is able to nitrate a range of tryptamine analogues, with significant alterations to the aromatic and aminoethyl moieties.

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

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

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.

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.