Theoretical Study on the Mechanism of the Oxygen Activation Process in Cysteine Dioxygenase Enzymes

Unveiling the mechanism of cysteamine dioxygenase: A combined HPLC-MS assay and metal-substitution approach

Mammalian cysteamine dioxygenase (ADO), a mononuclear non-heme Fe(II) enzyme with three histidine ligands, plays a key role in cysteamine catabolism and regulation of the N-degron signaling pathway. Despite its importance, the catalytic mechanism of ADO remains elusive. Here, we describe an HPLC-MS assay for characterizing thiol dioxygenase catalytic activities and a metal-substitution approach for mechanistic investigation using human ADO as a model. Two proposed mechanisms for ADO differ in oxygen activation: one involving a high-valent ferryl-oxo intermediate. We hypothesized that substituting iron with a metal that has a disfavored tendency to form high-valent states would discriminate between mechanisms. This chapter details the expression, purification, preparation, and characterization of cobalt-substituted ADO. The new HPLC-MS assay precisely measures enzymatic activity, revealing retained reactivity in the cobalt-substituted enzyme. The results obtained favor the concurrent dioxygen transfer mechanism in ADO. This combined approach provides a powerful tool for studying other non-heme iron thiol oxidizing enzymes.

New Frontiers in Nonheme Enzymatic Oxyferryl Species

Non-heme mononuclear iron dependent (NHM-Fe) enzymes exhibit exceedingly diverse catalytic reactivities. Despite their catalytic versatilities, the mononuclear iron centers in these enzymes show a relatively simple architecture, in which an iron atom is ligated with 2-4 amino acid residues, including histidine, aspartic or glutamic acid. In the past two decades, a common high-valent reactive iron intermediate, the S=2 oxyferryl (Fe(IV)-oxo or Fe(IV)=O) species, has been repeatedly discovered in NHM-Fe enzymes containing a 2-His-Fe or 2-His-1-carboxylate-Fe center. However, for 3-His/4-His-Fe enzymes, no common reactive intermediate has been identified. Recently, we have spectroscopically characterized the first S=1 Fe(IV) intermediate in a 3-His-Fe containing enzyme, OvoA, which catalyzes a novel oxidative carbon-sulfur bond formation. In this review, we summarize the broad reactivities demonstrated by S=2 Fe(IV)-oxo intermediates, the discovery of the first S=1 Fe(IV) intermediate in OvoA and the mechanistic implication of such a discovery, and the intrinsic reactivity differences of the S=2 and the S=1 Fe(IV)-oxo species. Finally, we postulate the possible reasons to utilize an S=1 Fe(IV) species in OvoA and their implications to other 3-His/4-His-Fe enzymes.

Thiol dioxygenases: from structures to functions

Thiol oxidation to dioxygenated sulfinic acid is catalyzed by an enzyme family characterized by a cupin fold. These proteins act on free thiol-containing molecules to generate central metabolism precursors and signaling compounds in bacteria, fungi, and animal cells. In plants and animals, they also oxidize exposed N-cysteinyl residues, directing proteins to proteolysis. Enzyme kinetics, X-ray crystallography, and spectroscopy studies prompted the formulation and testing of hypotheses about the mechanism of action and the different substrate specificity of these enzymes. Concomitantly, the physiological role of thiol dioxygenation in prokaryotes and eukaryotes has been studied through genetic and physiological approaches. Further structural characterization is necessary to enable precise and safe manipulation of thiol dioxygenases (TDOs) for therapeutic, industrial, and agricultural applications.

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.

Enhanced Reactivity through Equatorial Sulfur Coordination in Nonheme Iron(IV)-Oxo Complexes: Insights from Experiment and Theory

Sulfur ligation in metalloenzymes often gives the active site unique properties, whether it is the axial cysteinate ligand in the cytochrome P450s or the equatorial sulfur/thiol ligation in nonheme iron enzymes. To understand sulfur ligation to iron complexes and how it affects the structural, spectroscopic, and intrinsic properties of the active species and the catalysis of substrates, we pursued a systematic study and compared sulfur with amine-ligated iron(IV)-oxo complexes. We synthesized and characterized a biomimetic N4S-ligated iron(IV)-oxo complex and compared the obtained results with an analogous N5-ligated iron(IV)-oxo complex. Our work shows that the amine for sulfur replacement in the equatorial ligand framework leads to a rate enhancement for oxygen atom and hydrogen atom transfer reactions. Moreover, the sulfur-ligated iron(IV)-oxo complex reacts through a different reaction mechanism as compared to the N5-ligated iron(IV)-oxo complex, where the former reacts through hydride transfer with the latter reacting via radical pathways. We show that the reactivity differences are caused by a dramatic change in redox potential between the two complexes. Our studies highlight the importance of implementing a sulfur ligand into the equatorial ligand framework of nonheme iron(IV)-oxo complexes and how it affects the physicochemical properties of the oxidant and its reactivity.

Improved resolution of 3-mercaptopropionate dioxygenase active site provided by ENDOR spectroscopy offers insight into catalytic mechanism

3-mercaptopropionate (3MPA) dioxygenase (MDO) is a mononuclear nonheme iron enzyme that catalyzes the O2-dependent oxidation of thiol-bearing substrates to yield the corresponding sulfinic acid. MDO is a member of the cysteine dioxygenase family of small molecule thiol dioxygenases and thus shares a conserved sequence of active site residues (Serine-155, Histidine-157, and Tyrosine-159), collectively referred to as the SHY-motif. It has been demonstrated that these amino acids directly interact with the mononuclear Fe-site, influencing steady-state catalysis, catalytic efficiency, O2-binding, and substrate coordination. However, the underlying mechanism by which this is accomplished is poorly understood. Here, pulsed electron paramagnetic resonance spectroscopy [1H Mims electron nuclear double resonance spectroscopy] is applied to validate density functional theory computational models for the MDO Fe-site simultaneously coordinated by substrate and nitric oxide (NO), (3MPA/NO)-MDO. The enhanced resolution provided by electron nuclear double resonance spectroscopy allows for direct observation of Fe-bound substrate conformations and H-bond donation from Tyr159 to the Fe-bound NO ligand. Further inclusion of SHY-motif residues within the validated model reveals a distinct channel restricting movement of the Fe-bound NO-ligand. It has been argued that the iron-nitrosyl emulates the structure of potential Fe(III)-superoxide intermediates within the MDO catalytic cycle. While the merit of this assumption remains unconfirmed, the model reported here offers a framework to evaluate oxygen binding at the substrate-bound Fe-site and possible reaction mechanisms. It also underscores the significance of hydrogen bonding interactions within the enzymatic active site.

A Nonheme Iron(III) Superoxide Complex Leads to Sulfur Oxygenation

A new alkylthiolate-ligated nonheme iron complex, FeII(BNPAMe2S)Br (1), is reported. Reaction of 1 with O2 at -40 °C, or reaction of the ferric form with O2•- at -80 °C, gives a rare iron(III)-superoxide intermediate, [FeIII(O2)(BNPAMe2S)]+ (2), characterized by UV-vis, 57Fe Mössbauer, ATR-FTIR, EPR, and CSIMS. Metastable 2 then converts to an S-oxygenated FeII(sulfinate) product via a sequential O atom transfer mechanism involving an iron-sulfenate intermediate. These results provide evidence for the feasibility of proposed intermediates in thiol dioxygenases.

Mechanistic Insights into the N-Hydroxylations Catalyzed by the Binuclear Iron Domain of SznF Enzyme: Key Piece in the Synthesis of Streptozotocin

SznF, a member of the emerging family of heme-oxygenase-like (HO-like) di-iron oxidases and oxygenases, employs two distinct domains to catalyze the conversion of Nω-methyl-L-arginine (L-NMA) into N-nitroso-containing product, which can subsequently be transformed into streptozotocin. Using unrestricted density functional theory (UDFT) with the hybrid functional B3LYP, we have mechanistically investigated the two sequential hydroxylations of L-NMA catalyzed by SznF's binuclear iron central domain. Mechanism B primarily involves the O-O bond dissociation, forming Fe(IV)=O, induced by the H+/e- introduction to the FeA side of μ-1,2-peroxo-Fe2(III/III), the substrate hydrogen abstraction by Fe(IV)=O, and the hydroxyl rebound to the substrate N radical. The stochastic addition of H+/e- to the FeB side (mechanism C) can transition to mechanism B, thereby preventing enzyme deactivation. Two other competing mechanisms, involving the direct O-O bond dissociation (mechanism A) and the addition of H2O as a co-substrate (mechanism D), have been ruled out.

Enzyme-Catalyzed Oxidative Degradation of Ergothioneine

Ergothioneine is a sulfur-containing metabolite that is produced by bacteria and fungi, and is absorbed by plants and animals as a micronutrient. Ergothioneine reacts with harmful oxidants, including singlet oxygen and hydrogen peroxide, and may therefore protect cells against oxidative stress. Herein we describe two enzymes from actinobacteria that cooperate in the specific oxidative degradation of ergothioneine. The first enzyme is an iron-dependent thiol dioxygenase that produces ergothioneine sulfinic acid. A crystal structure of ergothioneine dioxygenase from Thermocatellispora tengchongensis reveals many similarities with cysteine dioxygenases, suggesting that the two enzymes share a common mechanism. The second enzyme is a metal-dependent ergothioneine sulfinic acid desulfinase that produces Nα-trimethylhistidine and SO2 . The discovery that certain actinobacteria contain the enzymatic machinery for O2 -dependent biosynthesis and O2 -dependent degradation of ergothioneine indicates that these organisms may actively manage their ergothioneine content.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.

Direct Electrophilic Attack of Compound I on the Indole Ring in the Peroxygenase Mechanism of Dehaloperoxidase DHP B in Degrading Haloindole: Electron Transfer Promotes the Reaction

The H2O2-dependent degradation of haloindole catalyzed by the dehaloperoxidase (DHP) from Amphitrite ornate has been reported to employ the peroxygenase mechanism, and the two oxidized products 5-halo-2-oxindole and 5-halo-3-oxindole have a similar amount. According to a previous experimental study, compound I (Cpd I) was suggested to be responsible for triggering the reaction, and the reaction may undergo three possible intermediates; however, the reaction details are still unclear. To clarify the reaction mechanism of DHP, the computational model was constructed on the basis of the high-resolution crystal structure, and a series of the quantum mechanical/molecular mechanical calculations were performed. Based on our calculation results, it is confirmed that the reaction starts from the direct electrophilic attack of Cpd I on the indole ring of the substrate, and the resulted intermediate contains both a carbocation and an oxygen anion, whereas the common hydrogen abstraction by Cpd I was calculated to correspond to a relatively higher barrier. In addition, a net electron transfer from the substrate to the iron center is observed during the attack of Cpd I on the indole ring; therefore, the carbocation/oxygen anion intermediate can easily undergo an intramolecular hydride transfer to form the product 5-halo-2-oxindole or isomerize to the epoxide intermediate which finally generates another product 5-halo-3-oxindole. It is the zwitterionic characteristic of the intermediate that makes the intermolecular hydride transfer quite easy, and it is the high electron affinity of the iron center that promotes the single-electron oxidation of the reaction intermediate. Our calculations well explain the formation of two oxidized products 5-halo-2-oxindole and 5-halo-3-oxindole.

Spectroscopic analysis of the mammalian enzyme cysteine dioxygenase

l-Cysteine (Cys) is an essential building block for the synthesis of new proteins and serves as a precursor for several biologically important sulfur-containing molecules, such as coenzyme A, taurine, glutathione, and inorganic sulfate. However, organisms must tightly regulate the concentration of free Cys, as elevated levels of this semi-essential amino acid can be extremely harmful. The non-heme iron enzyme cysteine dioxygenase (CDO) serves to maintain the proper levels of Cys by catalyzing its oxidation to cysteine sulfinic acid. Crystal structures of resting and substrate-bound mammalian CDO revealed two surprising structural motifs in the first and second coordination spheres of the Fe center. The first is the existence of a neutral three histidine (3-His) facial triad that coordinates the Fe ion, as opposed to an anionic 2-His-1-carboxylate facial triad that is typically observed in mononuclear non-heme Fe(II) dioxygenases. The second unusual structural feature exhibited by mammalian CDO is the presence of a covalent crosslink between the sulfur of a Cys residue and an ortho-carbon of a tyrosine residue. Spectroscopic studies of CDO have provided invaluable insights into the roles that these unusual features play with regards to substrate Cys and co-substrate O2 binding and activation. In this chapter, we summarize results obtained from electronic absorption, electron paramagnetic resonance, magnetic circular dichroism, resonance Raman, and Mössbauer spectroscopic studies of mammalian CDO carried out in the last two decades. Pertinent results obtained from complementary computational studies are also briefly summarized.

DFT insights into the mechanism of O2 activation catalyzed by a structural and functional model of cysteine dioxygenase with tris(2-pyridyl)methane-based ligand architecture

Cysteine dioxygenation is an important step in the metabolism of toxic L-cysteine (Cys) in the human body, carried out by cysteine dioxygenase enzyme (CDO). The disruption of this process is found to elicit neurological health issues. This work reports a computational investigation of mechanistic aspects of this reaction, using a recently reported tris(2-pyridyl)methane-based biomimetic model complex of CDO. The computed results indicate that, the initial SO2 bond formation process is the slowest step in the S-dioxygenation process, possessing an activation barrier of 12.7 kcal/mol. The remaining steps were found to be downhill requiring very small activation energies. The transition states were found to undergo spin crossover between triplet and quintet states, while the singlet surface remained unstable throughout the entire reaction. In essence, the mechanistic scheme and multistate reactivity pattern together with the relatively small computed rate-limiting activation barrier as well as the exothermic formation energy demonstrate that the model complex is an efficient biomimetic CDO model. In addition, the study also substantiates the involvement of Fe(IV)oxido intermediates in the mechanism of S-dioxygenation by the chosen model complex. The insights derived from the O₂ activation process might pave way for development of more accurate CDO model catalysts that might be capable of even more efficiently mimicking the geometric, spectroscopic and functional features of the CDO enzyme.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

Dioxygen Activation and Nδ,Nε-Dihydroxylation Mechanism Involved in the Formation of N-Nitrosourea Pharmacophore in Streptozotocin Catalyzed by Nonheme Diiron Enzyme SznF

SznF is a nonheme diiron-dependent enzyme that catalyzes the critical N-nitrosation involved in the formation of the N-nitrosourea moiety in the pancreatic cancer drug streptozotocin. The N-nitrosation contains two successive N-hydroxylation and N-nitrosation steps, which are carried out by two separate active sites, namely, the central domain and cupin domain. Recently, the crystal structure of SznF was obtained, and the central domain was proved to contain a diiron cofactor to catalyze the N-hydroxylation. In this work, to gain insights into the O₂ activation and the successive N-hydroxylation mechanism, on the basis of the high-resolution crystal structure, the enzyme-substrate complex models were constructed, and a series of combined QM/MM calculations were performed. Based on our calculations, the activation of O₂ starts from the diiron(II,III)-superoxo (S) to generate the diiron(IV)-oxo species (Q) via a diiron(III,III)-peroxo (P)-like transition state or unstable intermediate (P'), and species P' can be described as a hybridization of diiron(IV)-oxo species and diiron(III,III)-peroxo (P) owing to the long distances of Fe1-Fe2 (4.22 Å) and O1-O2 (1.89 Å), which is different from those of other nonheme diiron enzymes. In the following hydroxylation of N^δ and N^ε, the N^δ-hydroxylation was confirmed to occur first, agreeing with the experimental observations. Because the diiron(IV)-oxo species (Q) is responsible for hydroxylation, the reaction follows the H-abstraction/OH rebound mechanism, and the first abstraction occurs on the N^δ-H rather than N^ε-H, which may be attributed to the different orientation of Fe(IV)-oxo relative to N-H as well as the bond dissociation enthalpies of two N-H bonds. The hydroxylation of N-methyl-L-arginine does not employ the diiron(III,III)-hydroperoxo (P″) to trigger the electrophilic attack of the guanidine to directly form the N-O bond, as previously suggested. In addition, our calculations also revealed that the direct attack of the Fe(IV)═O unit to the N^δ of the substrate corresponds to a higher barrier than that in the H-abstraction/OH rebound mechanism. These results may provide useful information for understanding the formation of the di-hydroxylation intermediate involved in the biosynthesis of N-nitrosation.

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

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

Charge Maintenance during Catalysis in Nonheme Iron Oxygenases

Here, the choice of the first coordination shell of the metal center is analyzed from the perspective of charge maintenance in a binary enzyme-substrate complex and an O₂-bound ternary complex in the nonheme iron oxygenases. Comparing homogentisate 1,2-dioxygenase and gentisate dioxygenase highlights the significance of charge maintenance after substrate binding as an important factor that drives the reaction coordinate. We then extend the charge analysis to several common types of nonheme iron oxygenases containing either a 2-His-1-carboxylate facial triad or a 3-His or 4-His ligand motif, including extradiol and intradiol ring-cleavage dioxygenases, thiol dioxygenases, α-ketoglutarate-dependent oxygenases, and carotenoid cleavage oxygenases. After forming the productive enzyme-substrate complex, the overall charge of the iron complex at the 0, +1, or +2 state is maintained in the remaining catalytic steps. Hence, maintaining a constant charge is crucial to promote the reaction of the iron center beginning from the formation of the Michaelis or ternary complex. The charge compensation to the iron ion is tuned not only by protein-derived carboxylate ligands but also by substrates. Overall, these analyses indicate that charge maintenance at the iron center is significant when all the necessary components form a productive complex. This charge maintenance concept may apply to most oxygen-activating metalloenzymes systems that do not draw electrons and protons step-by-step from a separate reactant, such as NADH, via a reductase. The charge maintenance perception may also be useful in proposing catalytic pathways or designing prototypical reactions using artificial or engineered enzymes for biotechnological applications.

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) Å).

Computational Study of the C5-Hydroxylation Mechanism Catalyzed by the Diiron Monooxygenase PtmU3 as Part of the Platensimycin Biosynthesis

PtmU3 is a newly identified nonheme diiron monooxygenase, which installs a C-5 β-hydroxyl group into the C-19 CoA-ester intermediate involved in the biosynthesis of unique diterpene-derived scaffolds of platensimycin and platencin. PtmU3 possesses a noncanonical diiron active site architecture of a saturated six-coordinate iron center and lacks the μ-oxo bridge. Although the hydroxylation process is a simple reaction for nonheme mononuclear iron-dependent enzymes, how PtmU3 employs the diiron center to catalyze the H-abstraction and OH-rebound is still unknown. In particular, the electronic characteristic of diiron is also unclear. To understand the catalytic mechanism of PtmU3, we constructed two reactant models in which both the Fe1^II-Fe2^III-superoxo and Fe1^II-Fe2^IV═O are considered to trigger the H-abstraction and performed a series of quantum mechanics/molecular mechanics calculations. Our calculation results reveal that PtmU3 is a special monooxygenase, that is, both atoms of the dioxygen molecule can be incorporated into two molecules of the substrate by the successive reactions. In the first-round reaction, PtmU3 uses the Fe1^II-Fe2^III-superoxo to install a hydroxyl group into the substrate, generating the high-reactive Fe1^II-Fe2^IV═O complex. In the second-round reaction, the Fe1^II-Fe2^IV═O species is responsible for the hydroxylation of another molecule of the substrate. In the diiron center, Fe2 adopts the high spin state (S = 5/2) during the catalysis, whereas for Fe1, in addition to its structural role, it may also play an assistant role for Fe1 catalysis. In the two successive OH-installing steps, the H-abstraction is always the rate-liming step. E241 and D308 not only act as bridging ligands to connect two Fe ions but also take part in the electron reorganization. Owing to the high reactivity of Fe1^II-Fe2^IV═O compared to Fe1^II-Fe2^III-superoxo, besides the C5-hydroxylation, the C3- or C18-hydroxylation was also calculated to be feasible.

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.

Structure and Reactivity of a High-Spin, Nonheme Iron(III)- Superoxo Complex Supported by Phosphinimide Ligands

Nonheme iron oxygenases utilize dioxygen to accomplish challenging chemical oxidations. A further understanding of the Fe-O₂ intermediates implicated in these processes is challenged by their highly transient nature. To that end, we have developed a ligand platform featuring phosphinimide donors intended to stabilize oxidized, high-spin iron complexes. O₂ exposure of single crystals of a three-coordinate Fe(II) complex of this framework allowed for in crystallo trapping of a terminally bound Fe-O₂ complex suitable for XRD characterization. Spectroscopic and computational studies of this species support a high-spin Fe(III) center antiferromagnetically coupled to a superoxide ligand, similar to that proposed for numerous nonheme iron oxygenases. In addition to the apparent stability of this synthetic Fe-O₂ complex, its ability to engage in a range of stoichiometric and catalytic oxidation processes demonstrates that this iron-phosphinimide system is primed for development in modeling oxidizing bioinorganic intermediates and green oxidation chemistry.

Emerging roles for thiol dioxygenases as oxygen sensors

Cysteine dioxygenases, 3-mercaptopropionate dioxygenases and mercaptosuccinate dioxygenases are all thiol dioxygenases (TDOs) that catalyse oxidation of thiol molecules to sulphinates. They are Fe(II)-dependent dioxygenases with a cupin fold that supports a 3xHis metal-coordinating triad at the active site. They also have other, broadly common features including arginine residues involved in substrate carboxylate binding and a conserved trio of residues at the active site featuring a tyrosine important in substrate binding catalysis. Recently, N-terminal cysteinyl dioxygenase enzymes (NCOs) have been identified in plants (plant cysteine oxidases, PCOs), while human 2-aminoethanethiol dioxygenase (ADO) has been shown to act as both an NCO and a small molecule TDO. Although the cupin fold and 3xHis Fe(II)-binding triad seen in the small molecule TDOs are conserved in NCOs, other active site features and aspects of the overall protein architecture are quite different. Furthermore, the PCOs and ADO appear to act as biological O₂ sensors, as shown by kinetic analyses and hypoxic regulation of the stability of their biological targets (N-terminal cysteine oxidation triggers protein degradation via the N-degron pathway). Here, we discuss the emergence of these two subclasses of TDO including structural features that could dictate their ability to bind small molecule or polypeptide substrates. These structural features may also underpin the O₂ -sensing capability of the NCOs. Understanding how these enzymes interact with their substrates, including O₂ , could reveal strategies to manipulate their activity, relevant to hypoxic disease states and plant adaptive responses to flooding.

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.

Mechanistic Insights into the P450 TleB-Catalyzed Unusual Intramolecular C-N Bond Formation Involved in the Biosynthesis of Indolactam V

Indolactam V, a known biosynthetic precursor of indolactam alkaloids, is the main pharmacophore of this family and exhibits potential protein kinase C activation. A key step in the biosynthesis of indolactam V is the formation of an indole-fused nine-membered lactam core by intramolecular C-N bond formation. In this work, we report a computational study of the unique cytochrome P450 TleB enzyme-catalyzed direct and selective C-H bond amination reaction that can generate indolactam V from the dipeptide N-methylvalyl-tryptophanol. By performing molecular dynamics simulations and quantum-mechanical/molecular-mechanical calculations, we revealed that the C-H bond amination involves one step of proton transfer from N1-H of the indole ring to the Fe^IV═O unit, one step of hydrogen abstraction of N13-H in the side chain of the substrate by the Fe^IV-OH unit, and diradical coupling, in which two conformational changes of the side chain of the substrate are necessary. In the enzyme-substrate complex of TleB, the N-H bond of the indole ring of the substrate forms a strong hydrogen bond with the Fe^IV═O unit in compound I, and the porphyrin radical cation accepts an electron from the substrate to form the closed-shell electronic configuration. Thus, compound I in the enzyme-substrate complex cannot be described as Fe^IV═O coupled to a porphyrin radical cation, which is different from those of other P450 enzymes. Besides, two stages of conformational changes of the side chains of the substrate may increase the relative energies of reaction intermediates by 10-12 kcal/mol. From the structure point of view, it is the rotatable long side chain of the substrate and the large flexible active pocket of TleB that make the intramolecular diradical coupling feasible. Our findings may provide useful information to further understand the Tleb-catalyzed intramolecular C-H bond amination and the other bio-catalyzed intramolecular diradical coupling.

N-Nitrosation Mechanism Catalyzed by Non-heme Iron-Containing Enzyme SznF Involving Intramolecular Oxidative Rearrangement

The non-heme iron-dependent enzyme SznF catalyzes a critical N-nitrosation step during the N-nitrosourea pharmacophore biosynthesis in streptozotocin. The intramolecular oxidative rearrangement process is known to proceed at the Fe^II-containing active site in the cupin domain of SznF, but its mechanism has not been elucidated to date. In this study, based on the density functional theory calculations, a unique mechanism was proposed for the N-nitrosation reaction catalyzed by SznF in which a four-electron oxidation process is accomplished through a series of complicated electron transferring between the iron center and substrate to bypass the high-valent Fe^IV═O species. In the catalytic reaction pathway, the O₂ binds to the iron center and attacks on the substrate to form the peroxo bridge intermediate by obtaining two electrons from the substrate exclusively. Then, instead of cleaving the peroxo bridge, the C^ε-N^ω bond of the substrate is homolytically cleaved first to form a carbocation intermediate, which polarizes the peroxo bridge and promotes its heterolysis. After O-O bond cleavage, the following reaction steps proceed effortlessly so that the N-nitrosation is accomplished without NO exchange among reaction species.

Formation Mechanism of Cofactor Cys-Tyr in the Cysteine Dioxygenases (CDO and F2-CDO) and Its Influence on Catalysis: A QM/MM Study

Cysteine dioxygenase (CDO) is a nonheme mononuclear iron enzyme, which catalyzes the oxidation of cysteine to cysteine sulfinic acid. Crystal structure studies of mammalian CDO showed that there is a cross-linked cysteine-tyrosine (Cys-Tyr) cofactor in its active site. Moreover, the formation of the Cys-Tyr cofactor requires the metal cofactor (Fe²⁺) and O₂, and it was previously considered to substantially enhance the catalytic efficiency and half-life of CDO. Recently, a pure human CDO (F₂-CDO) without including the Cys-Tyr cofactor was crystalized by the site-directed mutagenesis approach in the anaerobic condition. In this work, to gain insights into the formation mechanism of the Cys-Tyr cofactor and whether it can really promote the catalytic reactivity of CDO, a series of computational models have been constructed, and quantum mechanical/molecular mechanical (QM/MM) calculations have been performed. Our calculation results reveal that WT-CDO and F₂-CDO follow different mechanisms for the formation of the Cys-Tyr cofactor. In F₂-CDO, the cofactor formation contains the H-abstraction, C-S bond formation, intramolecular F migration, and aromatization of the residue F₂Y157, in which the Fe-coordinate dioxygen can be recovered after the formation cofactor; however, in the WT-CDO, the cofactor formation shows some differences. During the reaction, hydrogen peroxide is generated, and the final aromatization requires the assistance of one water molecule. Furthermore, the overall barriers of cofactor formation are always higher than l-cysteine oxidation for both WT-CDO and F₂-CDO irrespective of the absence or presence of the cofactor. Thus, we can theoretically confirm that the Cys-Tyr cofactor is not essential for the oxidation activity of CDO, and cofactor formation is just an accompanying reaction but not a prerequisite for the oxidation reaction. These results may provide useful information for understanding the catalysis of CDO.

Superoxide Oxidation by a Thiolate-Ligated Iron Complex and Anion Inhibition

Superoxide (O₂^•-) is a toxic radical, generated via the adventitious reduction of dioxygen (O₂), which has been implicated in a number of human disease states. Nonheme iron enzymes, superoxide reductase (SOR) and superoxide dismutase (SOD), detoxify O₂^•- via reduction to afford H₂O₂ and disproportionation to afford O₂ and H₂O₂, respectively. The former contains a thiolate in the coordination sphere, which has been proposed to prevent O₂^•- oxidation to O₂. The work described herein shows that, in contrast to this, oxidized thiolate-ligated [Fe^III(S^Me2N₄(tren)(THF)]²⁺ (1^ox-THF) is capable of oxidizing O₂^•- to O₂. Coordinating anions, Cl^- and OAc^-, are shown to inhibit dioxygen evolution, implicating an inner-sphere mechanism. Previously we showed that the reduced thiolate-ligated [Fe^II(S^Me2N₄(tren))]⁺ (1) is capable of reducing O₂^•- via a proton-dependent inner-sphere mechanism involving a transient Fe(III)-OOH intermediate. A transient ferric-superoxo intermediate, [Fe^III(S^Me2N₄(tren))(O₂)]⁺ (3), is detected by electronic absorption spectroscopy at -130 °C in the reaction between 1^ox-THF and KO₂ and shown to evolve O₂ upon slight warming to -115 °C. The DFT calculated O-O (1.306 Å) and Fe-O (1.943 Å) bond lengths of 3 are typical of ferric-superoxo complexes, and the time-dependent DFT calculated electronic absorption spectrum of 3 reproduces the experimental spectrum. The electronic structure of 3 is shown to consist of two antiferromagnetically coupled (J^calc = -180 cm^-1) unpaired electrons, one in a superoxo π*(O-O) orbital and the other in an antibonding π*(Fe(d_yz)-S(p_y)) orbital.

A bioinspired oxoiron(IV) motif supported on a N2S2 macrocyclic ligand

A mononuclear oxoiron(iv) complex 1-trans bearing two equatorial sulfur ligations is synthesized and characterized as an active-site model of the elusive sulfur-ligated FeIV[double bond, length as m-dash]O intermediates in non-heme iron oxygenases. The introduction of sulfur ligands weakens the Fe[double bond, length as m-dash]O bond and enhances the oxidative reactivity of the FeIV[double bond, length as m-dash]O unit with a diminished deuterium kinetic isotope effect, thereby providing a compelling rationale for nature's use of the cis-thiolate ligated oxoiron(iv) motif in key metabolic transformations.

How external perturbations affect the chemoselectivity of substrate activation by cytochrome P450 OleTJE

Cytochrome P450 enzymes are versatile biocatalysts found in most forms of life. Generally, the cytochrome P450s react with dioxygen and hence are haem-based mono-oxygenases; however, in specific isozymes, H2O2 rather than O2 is used and these P450s act as peroxygenases. The P450 OleTJE is a peroxygenase that binds long to medium chain fatty acids and converts them to a range of products originating from Cα-hydroxylation, Cβ-hydroxylation, Cα-Cβ desaturation and decarboxylation of the substrate. There is still controversy regarding the details of the reaction mechanism of P450 OleTJE; how the products are formed and whether the product distributions can be influenced by external perturbations. To gain further insights into the structure and reactivity of P450 OleTJE, we set up a range of large active site model complexes as well as full enzymatic structures and did a combination of density functional theory studies and quantum mechanics/molecular mechanics calculations. In particular, the work focused on the mechanisms leading to these products under various reaction conditions. Thus, for a small cluster model, we find a highly selective Cα-hydroxylation pathway that is preferred over Cβ-H hydrogen atom abstraction by at least 10 kcal mol-1. Introduction of polar residues to the model, such as an active site protonated histidine residue or through external electric field effects, lowers the Cβ-H hydrogen atom abstraction barriers are lowered, while a full QM/MM model brings the Cα-H and Cβ-H hydrogen atom abstraction barriers within 1 kcal mol-1. Our studies; therefore, implicate that environmental effects in the second-coordination sphere can direct and guide selectivities in enzymatic reaction mechanisms.

Drastically inhibited nZVI-Fenton oxidation of organic pollutants by cysteine: Multiple roles in the nZVI/O2/hv system

The influence of l-cysteine, a common aliphatic amino acid, on the zero-valent iron (nZVI)/O₂ photo-Fenten degradation of rhodamine B (RhB) was investigated in this study. The oxidation rate of RhB in the nZVI/O₂/hv system was 91.2% after 40 min under the illumination and oxygen conditions and pH of 3, but when cysteine was introduced into the system, the oxidization process was inhibited. The removal of RhB was only about 50% after 40 min at a cysteine concentration ≥50 μM. It was shown experimentally that, under dark conditions, only 40.5% and 19.8% RhB was removed by the nZVI/O₂ and nZVI/O₂/cysteine systems, respectively. Electron paramagnetic resonance (EPR) and iron dissolving experiments revealed that the addition of cysteine clearly reduced the production of hydroxyl radicals (OH) and Fe²⁺ and Fe³⁺. In addition, Fourier transform infrared spectroscopy (FTIR) demonstrated that cysteine could form hydrogen bonds on the iron surface. These results indicated that the main inhibition mechanism of cysteine was the alleviation of the oxidation of nZVI to Fe²⁺ and Fe³⁺ through wrapping the nZVI particles. Moreover, cystine (the oxidized form of CYS) could partly react with OH to regenerate cysteine, which resulted in competition with RhB for OH. Another possible reason for the inhibitory effect of cysteine was the prevention of light utilization. These findings indicate a non-negligible inhibitory trait for heterogeneous Fenton process in wastewater treatment when amino acids are present.

The molecular mechanism of P450-catalyzed amination of the pyrrolidine derivative of lidocaine: insights from multiscale simulations

Nitrogen heterocycles are key and prevalent motifs in drugs. Evolved variants of cytochrome P450_BM3 (CYP102A1) from Bacillus megaterium employ high-valent oxo-iron(iv) species to catalyze the synthesis of imidazolidine-4-ones via an intramolecular C–H amination. Herein, we use multi-scale simulations, including classical molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical (QM/MM) calculations and QM calculations, to reveal the molecular mechanism of the intramolecular C–H amination of the pyrrolidine derivative of lidocaine bearing cyclic amino moieties catalyzed by the variant RP/FV/EV of P450_BM3, which bears five mutations compared to wild type. Our calculations show that overall catalysis includes both the enzymatic transformation in P450 and non-enzymatic transformation in water solution. The enzymatic transformation involves the exclusive hydroxylation of the C–H bond of the pyrrolidine derivative of lidocaine, leading to the hydroxylated intermediate, during which the substrate radical would be bypassed. The following dehydration and C–N coupling reactions are found to be much favored in aqueous situation compared to that in the non-polar protein environment. The present findings expand our understanding of the P450-catalyzed C(sp₃)–H amination reaction.

Nitrogen heterocycles are key and prevalent motifs in drugs.

Computational studies of DNA base repair mechanisms by nonheme iron dioxygenases: selective epoxidation and hydroxylation pathways

DNA base repair mechanisms of alkylated DNA bases is an important reaction in chemical biology and particularly in the human body. It is typically catalyzed by an α-ketoglutarate-dependent nonheme iron dioxygenase named the AlkB repair enzyme. In this work we report a detailed computational study into the structure and reactivity of AlkB repair enzymes with alkylated DNA bases. In particular, we investigate the aliphatic hydroxylation and C[double bond, length as m-dash]C epoxidation mechanisms of alkylated DNA bases by a high-valent iron(iv)-oxo intermediate. Our computational studies use quantum mechanics/molecular mechanics methods on full enzymatic structures as well as cluster models on active site systems. The work shows that the iron(iv)-oxo species is rapidly formed after dioxygen binding to an iron(ii) center and passes a bicyclic ring structure as intermediate. Subsequent cluster models explore the mechanism of substrate hydroxylation and epoxidation of alkylated DNA bases. The work shows low energy barriers for substrate activation and consequently energetically feasible pathways are predicted. Overall, the work shows that a high-valent iron(iv)-oxo species can efficiently dealkylate alkylated DNA bases and return them into their original form.

Mechanistic Insights into the Oxidative Rearrangement Catalyzed by the Unprecedented Dioxygenase ChaP Involved in Chartreusin Biosynthesis

ChaP is a non-heme iron-dependent dioxygenase belonging to the vicinal oxygen chelate (VOC) enzyme superfamily that catalyzes the final α-pyrone ring formation in the biosynthesis of chartreusin. In contrast to other common dioxygenases, for example, 2,3-catechol dioxygenase which uses the dioxygen molecule as the oxidant, ChaP requires the flavin-activated oxygen (O₂^2-) as the equivalent. Previous experiments showed that the ChaP-catalyzed ring rearrangement contains two successive C-C bond cleavages and one lactonization; however, the detailed reaction mechanism is unknown. In this work, on the basis of the recently obtained crystal structure of ChaP, the computational model was constructed and the catalytic mechanism of ChaP was explored by performing quantum mechanical/molecular mechanical (QM/MM) calculations. Our calculation results reveal that ChaP uses the proximal oxygen in iron-coordinated HOO^- to attack the carbonyl carbon of the substrate, whereas the previous proposal that Asp49 acts as a base to deprotonate the iron-coordinated HOO^- to generate O₂^2- is unlikely. In the first stage reaction, owing to the coordination of the substrate with iron, the substrate is activated by accepting an electron from iron and the resulting oxy-radical intermediate formed by O-O cleavage can easily trigger the ring rearrangement. In the final decarboxylation, the phenolic anion of the substrate cooperatively accepts the proton of iron-coordinated HOO^- to facilitate the attack of the distal oxygen, and the proton-coupled electron transfer (PCET) from the substrate to the Fe^IV═O plays a key role for the decarboxylation. These findings may provide useful information for understanding the ChaP-catalyzed oxidative rearrangement and other flavin-dependent non-heme dioxygenases.

Enhancing Tris(pyrazolyl)borate-Based Models of Cysteine/Cysteamine Dioxygenases through Steric Effects: Increased Reactivities, Full Product Characterization and Hints to Initial Superoxide Formation

The design of biomimetic model complexes for the cysteine dioxygenase (CDO) and cysteamine dioxygenase (ADO) is reported, where the 3-His coordination of the iron ion is simulated by three pyrazole donors of a trispyrazolyl borate ligand (Tp) and protected cysteine and cysteamine represent substrate ligands. It is found that the replacement of phenyl groups-attached at the 3-positions of the pyrazole units in a previous model-by mesityl residues has massive consequences, as the latter arrange to a more spacious reaction pocket. Thus, the reaction with O2 proceeds much faster and afterwards the first structural characterization of an iron(II) η2 -O,O-sulfinate product became possible. If one of the three Tp-mesityl groups is placed in the 5-position, an even larger reaction pocket results, which leads to yet faster rates and accumulation of a reaction intermediate at low temperatures, as shown by UV/Vis and Mössbauer spectroscopy. After comparison with the results of investigations on the cobalt analogues this intermediate is tentatively assigned to an iron(III) superoxide species.

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

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

Mechanism of Uncoupled Carbocyclization and Epimerization Catalyzed by Two Non-Heme Iron/α-Ketoglutarate Dependent Enzymes

The non-heme iron/α-ketoglutarate dependent enzymes SnoK and SnoN from Streptomyces nogalater are involved in the biosynthesis of anthracycline nogalamycin. Although they have similar active sites, SnoK is responsible for carbocyclization whereas SnoN solely catalyzes the hydroxyl epimerization. Herein, we performed docking, molecular simulations, and a series of combined quantum mechanics and molecular mechanics (QM/MM) calculations to illuminate the mechanisms of two enzymes. The catalytic reactions of two enzymes occur on the quintet state surface. For SnoK, the whole reaction includes two separated hydrogen-abstraction steps and one radical addition, and the latter step is calculated to be rate limiting with an energy barrier of 21.7 kcal/mol. Residue D106 is confirmed to participate in the construction of the hydrogen bond network, which plays a crucial role in positioning the bulky substrate in a specific orientation. Moreover, it is found that SnoN is only responsible for the hydrogen abstraction of the intermediate, and no residue was suggested to be suitable for donating a hydrogen atom to the substrate radical, which further confirms the suggestion based on experiments that either a cellular reductant or another enzyme protein could donate a hydrogen atom to the substrate. Our docking results coincide with the previous structural study that the different roles of two enzymes are achieved by minor changes in the alignment of the substrates in front of the reactive ferryl-oxo species. This work highlights the reaction mechanisms catalyzed by SnoK and SnoN, which is helpful for engineering the enzymes for the biosynthesis of anthracycline nogalamycin.

Mechanism of fatty acid decarboxylation catalyzed by a non-heme iron oxidase (UndA): a QM/MM study

UndA is a non-heme iron enzyme that was recognized to catalyze the decarboxylation of medium chain (C10-C14) fatty acids to produce trace amounts of 1-alkenes. Owing to the electron imbalance during the oxidative decarboxylation of the substrate and the reduction of O₂, only single turnover reactions were obtained in UndA in vitro assays. Unlike the general non-heme iron enzymes, the catalytic efficiency of UndA is quite low. According to the previous proposal, both Fe^III-OO˙^- and Fe^IV[double bond, length as m-dash]O complexes may abstract the β-H of fatty acids to trigger the oxidative decarboxylation reaction. Herein, on the basis of the crystal structures of UndA in complex with the substrate analogues, we constructed a series of computational models and performed quantum mechanics/molecular mechanics (QM/MM) calculations to explore the UndA-catalyzed decarboxylation using lauric acid as the substrate. Our calculation results reveal that only the Fe^III-OO˙^- complex can initiate the decarboxylation, and the substrate (lauric acid) should monodentately coordinate to the Fe center to facilitate the β-H abstraction. In addition, the monodentate coordination corresponds to higher relative energy than the bidentate mode, which may explain the low efficiency of UndA. It is also revealed that as long as the β-H is extracted by the Fe^III-OO˙^-, the decarboxylation of the substrate radical is quite easy, and an electron transfer from the substrate to the iron center is the prerequisite. For the Fe^IV[double bond, length as m-dash]O complex, since the β-H is far from the O_Fe atom and the angle of ∠Fe-O-H is 53.1°, the H-abstraction is calculated to be difficult.

Insights into the Mechanism and Enantioselectivity in the Biosynthesis of Ergot Alkaloid Cycloclavine Catalyzed by Aj_EasH from Aspergillus japonicus

Cycloclavine is a complex ergot alkaloid containing an unusual cyclopropyl moiety, which has a wide range of biological activities and pharmaceutical applications. The biosynthesis of cycloclavine requires a series of enzymes, one of which is a nonheme Fe^II/α-ketoglutarate-dependent (aKG) oxidase (Aj_EasH). According to the previous proposal, the cyclopropyl ring formation catalyzed by Aj_EasH follows an unprecedented oxidative mechanism; however, the reaction details are unknown. In this article, on the basis of the recently obtained crystal structure of Aj_EasH (EasH from Aspergillus japonicas), the reactant models were built, and the reaction details were investigated by performing QM-only and combined QM and MM calculations. Our calculation results reveal that the biosynthesis of cyclopropyl moiety involves a radical intermediate rather than a carbocationic or carbanionic intermediate as in the biosynthesis of terpenoid family. The iron(IV)-oxo first abstracts a hydrogen atom from the substrate to trigger the reaction, and then the generated radical intermediate undergoes ring rearrangement to form the fused 5-3 ring system of cycloclavine. On the basis of our calculations, the absolute configuration of the cycloclavine catalyzed by Aj_EasH from Aspergillus japonicus should be (5R,8R,10R), which is different from the product isolated from Ipomoea hildebrandtii (5R,8S,10S). Residues at the active site play an important role in substrate binding, ring rearrangement, and enantioselectivity.

Phenylalanine 4-monooxygenase: the "sulfoxidation polymorphism"

1. Consistent differences in the proportion of an orally administered dose of S-carboxymethyl-l-cysteine subsequently excreted in the urine as S-oxide metabolites were reported 40 years ago. This observation suggested the existence of inter-individual variation in the ability to undertake the enzymatic S-oxygenation of this compound. Pedigree studies and investigations employing twin pairs indicated a genetically controlled phenomenon overlaid with environmental influences. It was reproducible and not related to gender or age.2. Studies undertaken in several healthy volunteer cohorts always provided similar results that were not significantly different when statistically analysed. However, when compared to these healthy populations, a preponderance of subjects exhibiting the characteristic of poor sulfoxidation of S-carboxymethyl-l-cysteine was found within groups of patients suffering from various disease conditions. The most striking of these associations were witnessed amongst subjects diagnosed with neurodegenerative disorders; although, underlying mechanisms were unknown.3. Exhaustive investigation has identified the enzyme responsible for this S-oxygenation reaction as the tetrahydrobiopterin-dependent aromatic amino acid hydroxylase, phenylalanine 4-monooxygenase classically assigned the sole function of converting phenylalanine to tyrosine. The underlying principle is discussed that enzymes traditionally associated solely with intermediary metabolism may have as yet unrecognised alternative roles in protecting the organism from potential toxic assault.

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

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

Advances in Sustainable Catalysis: A Computational Perspective

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

Di- and Tetrairon(III) μ-Oxido Complexes of an N3S-Donor Ligand: Catalyst Precursors for Alkene Oxidations

The new di- and tetranuclear Fe(III) μ-oxido complexes [Fe₄(μ-O)₄(PTEBIA)₄](CF₃SO₃)₄(CH₃CN)₂] (1a), [Fe₂(μ-O)Cl₂(PTEBIA)₂](CF₃SO₃)₂ (1b), and [Fe₂(μ-O)(HCOO)₂(PTEBIA)₂](ClO₄)₂ (MeOH) (2) were prepared from the sulfur-containing ligand (2-((2,4-dimethylphenyl)thio)-N,N-bis ((1-methyl-benzimidazol-2-yl)methyl)ethanamine (PTEBIA). The tetrairon complex 1a features four μ-oxido bridges, while in dinuclear 1b, the sulfur moiety of the ligand occupies one of the six coordination sites of each Fe(III) ion with a long Fe-S distance of 2.814(6) Å. In 2, two Fe(III) centers are bridged by one oxido and two formate units, the latter likely formed by methanol oxidation. Complexes 1a and 1b show broad sulfur-to-iron charge transfer bands around 400–430 nm at room temperature, consistent with mononuclear structures featuring Fe-S interactions. In contrast, acetonitrile solutions of 2 display a sulfur-to-iron charge transfer band only at low temperature (228 K) upon addition of H₂O₂/CH₃COOH, with an absorption maximum at 410 nm. Homogeneous oxidative catalytic activity was observed for 1a and 1b using H₂O₂ as oxidant, but with low product selectivity. High valent iron-oxo intermediates could not be detected by UV-vis spectroscopy or ESI mass spectrometry. Rather, evidence suggest preferential ligand oxidation, in line with the relatively low selectivity and catalytic activity observed in the reactions.

An Alternative Active Site Architecture for O2 Activation in the Ergothioneine Biosynthetic EgtB from Chloracidobacterium thermophilum

Sulfoxide synthases are nonheme iron enzymes that catalyze oxidative carbon-sulfur bond formation between cysteine derivatives and N-α-trimethylhistidine as a key step in the biosynthesis of thiohistidines. The complex catalytic mechanism of this enzyme reaction has emerged as the controversial subject of several biochemical and computational studies. These studies all used the structure of the γ-glutamyl cysteine utilizing sulfoxide synthase, MthEgtB from Mycobacterium thermophilum (EC 1.14.99.50), as a structural basis. To provide an alternative model system, we have solved the crystal structure of CthEgtB from Chloracidobacterium thermophilum (EC 1.14.99.51) that utilizes cysteine as a sulfur donor. This structure reveals a completely different configuration of active site residues that are involved in oxygen binding and activation. Furthermore, comparison of the two EgtB structures enables a classification of all ergothioneine biosynthetic EgtBs into five subtypes, each characterized by unique active-site features. This active site diversity provides an excellent platform to examine the catalytic mechanism of sulfoxide synthases by comparative enzymology, but also raises the question as to why so many different solutions to the same biosynthetic problem have emerged.

A Nonheme Thiolate-Ligated Cobalt Superoxo Complex: Synthesis and Spectroscopic Characterization, Computational Studies, and Hydrogen Atom Abstraction Reactivity

The synthesis and characterization of a Co(II) dithiolato complex Co(Me₃TACN)(S₂SiMe₂) (1) are reported. Reaction of 1 with O₂ generates a rare thiolate-ligated cobalt-superoxo species Co(O₂)(Me₃TACN)(S₂SiMe₂) (2) that was characterized spectroscopically and structurally by resonance Raman, EPR, and X-ray absorption spectroscopies as well as density functional theory. Metal-superoxo species are proposed to S-oxygenate metal-bound thiolate donors in nonheme thiol dioxygenases, but 2 does not lead to S-oxygenation of the intramolecular thiolate donors and does not react with exogenous sulfur donors. However, complex 2 is capable of oxidizing the O-H bonds of 2,2,6,6-tetramethylpiperidin-1-ol derivatives via H atom abstraction. Complementary proton-coupled electron-transfer reactivity is seen for 2 with separated proton/reductant pairs. The reactivity studies indicate that 2 can abstract H atoms from weak X-H bonds with bond dissociation free energy (BDFE) ≤ 70 kcal mol^-1. DFT calculations predict that the putative Co(OOH) product has an O-H BDFE = 67 kcal mol^-1, which matches the observed pattern of reactivity seen for 2. These data provide new information regarding the selectivity of S-oxygenation versus H atom abstraction in thiolate-ligated nonheme metalloenzymes that react with O₂.