Catalytic Mechanism of Cofactor-Free Dioxygenases and How They Circumvent Spin-Forbidden Oxygenation of Their Substrates

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

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

An Enzymatic Oxidation Cascade Converts δ-Thiolactone Anthracene to Anthraquinone in the Biosynthesis of Anthraquinone-Fused Enediynes

Anthraquinone-fused enediynes are anticancer natural products featuring a DNA-intercalating anthraquinone moiety. Despite recent insights into anthraquinone-fused enediyne (AQE) biosynthesis, the enzymatic steps involved in anthraquinone biogenesis remain to be elucidated. Through a combination of in vitro and in vivo studies, we demonstrated that a two-enzyme system, composed of a flavin adenine dinucleotide (FAD)-dependent monooxygenase (DynE13) and a cofactor-free enzyme (DynA1), catalyzes the final steps of anthraquinone formation by converting δ-thiolactone anthracene to hydroxyanthraquinone. We showed that the three oxygen atoms in the hydroxyanthraquinone originate from molecular oxygen (O2), with the sulfur atom eliminated as H2S. We further identified the key catalytic residues of DynE13 and A1 by structural and site-directed mutagenesis studies. Our data support a catalytic mechanism wherein DynE13 installs two oxygen atoms with concurrent desulfurization and decarboxylation, whereas DynA1 acts as a cofactor-free monooxygenase, installing the final oxygen atom in the hydroxyanthraquinone. These findings establish the indispensable roles of DynE13 and DynA1 in AQE biosynthesis and unveil novel enzymatic strategies for anthraquinone formation.

Impact of Cytochrome P450 Enzyme on Fruit Quality

Cytochrome P450 enzymes are monooxygenases widely diffused in nature ranging from viruses to man. They can catalyze a very wide range of reactions, including the ketonization of C-H bonds, N/O/S-dealkylation, C-C bond cleavage, N/S-oxidation, hydroxylation, and the epoxidation of C=C bonds. Their versatility makes them valuable across various fields such as medicine, chemistry, and food processing. In this review, we aim to highlight the significant contribution of P450 enzymes to fruit quality, with a specific focus on the ripening process, particularly in grapevines. Grapevines are of particular interest due to their economic importance in the fruit industry and their significance in winemaking. Understanding the role of P450 enzymes in grapevine fruit ripening can provide insights into enhancing grape quality, flavor, and aroma, which are critical factors in determining the market value of grapes and derived products like wine. Moreover, the potential of P450 enzymes extends beyond fruit ripening. They represent promising candidates for engineering crop species that are resilient to both biotic and abiotic stresses. Their involvement in metabolic engineering offers opportunities for enhancing fruit quality attributes, such as taste, nutritional content, and shelf life. Harnessing the capabilities of P450 enzymes in crop improvement holds immense promise for sustainable agriculture and food security.

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

Cofactor-independent C-C bond cleavage reactions catalyzed by the AlpJ family of oxygenases in atypical angucycline biosynthesis

Biosynthesis of atypical angucyclines involves unique oxidative B-ring cleavage and rearrangement reactions, which are catalyzed by AlpJ-family oxygenases, including AlpJ, JadG, and GilOII. Prior investigations established the essential requirement for FADH2/FMNH2 as cofactors when utilizing the quinone intermediate dehydrorabelomycin as a substrate. In this study, we unveil a previously unrecognized facet of these enzymes as cofactor-independent oxygenases when employing the hydroquinone intermediate CR1 as a substrate. The enzymes autonomously drive oxidative ring cleavage and rearrangement reactions of CR1, yielding products identical to those observed in cofactor-dependent reactions of AlpJ-family oxygenases. Furthermore, the AlpJ- and JadG-catalyzed reactions of CR1 could be quenched by superoxide dismutase, supporting a catalytic mechanism wherein the substrate CR1 reductively activates molecular oxygen, generating a substrate radical and the superoxide anion O2 •-. Our findings illuminate a substrate-controlled catalytic mechanism of AlpJ-family oxygenases, expanding the realm of cofactor-independent oxygenases. Notably, AlpJ-family oxygenases stand as a pioneering example of enzymes capable of catalyzing oxidative reactions in either an FADH2/FMNH2-dependent or cofactor-independent manner.

Cofactorless oxygenases guide anthraquinone-fused enediyne biosynthesis

The anthraquinone-fused enediynes (AFEs) combine an anthraquinone moiety and a ten-membered enediyne core capable of generating a cytotoxic diradical species. AFE cyclization is triggered by opening the F-ring epoxide, which is also the site of the most structural diversity. Previous studies of tiancimycin A, a heavily modified AFE, have revealed a cryptic aldehyde blocking installation of the epoxide, and no unassigned oxidases could be predicted within the tnm biosynthetic gene cluster. Here we identify two consecutively acting cofactorless oxygenases derived from methyltransferase and α/β-hydrolase protein folds, TnmJ and TnmK2, respectively, that are responsible for F-ring tailoring in tiancimycin biosynthesis by comparative genomics. Further biochemical and structural characterizations reveal that the electron-rich AFE anthraquinone moiety assists in catalyzing deformylation, epoxidation and oxidative ring cleavage without exogenous cofactors. These enzymes therefore fill important knowledge gaps for the biosynthesis of this class of molecules and the underappreciated family of cofactorless oxygenases.

Evolutionary adaptation from hydrolytic to oxygenolytic catalysis at the α/β-hydrolase fold

Protein fold adaptation to novel enzymatic reactions is a fundamental evolutionary process. Cofactor-independent oxygenases degrading N-heteroaromatic substrates belong to the α/β-hydrolase (ABH) fold superfamily that typically does not catalyze oxygenation reactions. Here, we have integrated crystallographic analyses under normoxic and hyperoxic conditions with molecular dynamics and quantum mechanical calculations to investigate its prototypic 1-H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase (HOD) member. O2 localization to the "oxyanion hole", where catalysis occurs, is an unfavorable event and the direct competition between dioxygen and water for this site is modulated by the "nucleophilic elbow" residue. A hydrophobic pocket that overlaps with the organic substrate binding site can act as a proximal dioxygen reservoir. Freeze-trap pressurization allowed the structure of the ternary complex with a substrate analogue and O2 bound at the oxyanion hole to be determined. Theoretical calculations reveal that O2 orientation is coupled to the charge of the bound organic ligand. When 1-H-3-hydroxy-4-oxoquinaldine is uncharged, O2 binds with its molecular axis along the ligand's C2-C4 direction in full agreement with the crystal structure. Substrate activation triggered by deprotonation of its 3-OH group by the His-Asp dyad, rotates O2 by approximately 60°. This geometry maximizes the charge transfer between the substrate and O2, thus weakening the double bond of the latter. Electron density transfer to the O2(π*) orbital promotes the formation of the peroxide intermediate via intersystem crossing that is rate-determining. Our work provides a detailed picture of how evolution has repurposed the ABH-fold architecture and its simple catalytic machinery to accomplish metal-independent oxygenation.

The case for an oxidopyrylium intermediate in the mechanism of quercetin dioxygenases

The quercetin dioxygenases (QDOs) are unusual metalloenzymes in that they display ring-opening dioxygenase activity with several different first-row transition metal ions which do not undergo redox changes during turnover. The QDOs are also unique in that the substrate binds as an η1-flavonolate rather than the η2 -bidentate mode seen in all reported model complexes. The flavonol substrates were early examples of excited state intramolecular proton transfer (ESIPT) phenomena, in which photoexcitation causes an H-atom exchange between the adjacent hydroxyl and ketone, generating an oxidopyrylium emissive state. These oxidopyryliums undergo ring-opening dioxygenations analogous to the enzymatic reactions. Our hypothesis is that lability of the divalent metal ion may allow access to a reactive oxidopyrylium intermediate via coordination switching from the oxy to ketone position, which allows reaction with O2. In this report, we use a straight-forward methylation strategy to generate a panel of flavonol and thioflavonol derivatives modeling several η1- and η2-coordination modes. Methylation of 3-hydroxythioflavone generates an air stable η1 hydroxopyrylium salt, which undergoes rapid ring-opening dioxygenation by deprotonation or photoexcitation. By comparison, the η1-methoxyflavonol does not react with O2 under any condition. We find that any of the studied flavonol derivatives, η1 or η2, which demonstrates ESIPT-like oxidopyrylium emissions undergo QDO-like ring-opening reactions with dioxygen. The implications of these results concerning the mechanism of QDOs and related dioxygenases is discussed.

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.

Spin-Forbidden Addition of Molecular Oxygen to Stable Enol Intermediates-Decarboxylation of 2-Methyl-1-tetralone-2-carboxylic Acid

The deprotonation of an organic substrate is a common preactivation step for the enzymatic cofactorless addition of O2 to this substrate, as it promotes charge-transfer between the two partners, inducing intersystem crossing between the triplet and singlet states involved in the process. Nevertheless, the spin-forbidden addition of O2 to uncharged ligands has also been observed in the laboratory, and the detailed mechanism of how the system circumvents the spin-forbiddenness of the reaction is still unknown. One of these examples is the cofactorless peroxidation of 2-methyl-3,4-dihydro-1-naphthol, which will be studied computationally using single and multi-reference electronic structure calculations. Our results show that the preferred mechanism is that in which O2 picks a proton from the substrate in the triplet state, and subsequently hops to the singlet state in which the product is stable. For this reaction, the formation of the radical pair is associated with a higher barrier than that associated with the intersystem crossing, even though the absence of the negative charge leads to relatively small values of the spin-orbit coupling.

Cofactor-Free Dioxygenases-Catalyzed Reaction Pathway via Proton-Coupled Electron Transfer

Understanding the general mechanism of the metal-free and cofactor-free oxidases and oxygenases catalyzed activation of triplet O₂ is one of the most challenging questions in the field of enzymatic catalysis. Herein, we have performed Quantum Mechanics/Molecular Mechanics (QM/MM) multiscale simulations to reveal the detailed mechanism of the HOD catalyzed (i.e., 1-H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase from Arthrobacter nitroguajacolicus Rü61a) decomposition of N-heteroaromatic compounds. The complete catalytic mechanism includes four steps: (1) proton transfer from 1-H-3-hydroxy-4-oxoquinaldine (QND) substrate to His251 residue coupled with an electron transfer from QND to triplet O₂ (i.e., PCET), (2) formation of C-O bond via an open-shell singlet diradical recombination pathway, (3) ring-closure to form a bicyclic ring, and (4) dissociation of CO. The dissociation of CO is determined as the rate-limiting step, and its calculated energy barrier of 14.9 kcal/mol is consistent with the 15.5 kcal/mol barrier derived from experimental kinetic data. The mechanistic profile is not only valuable for understanding the fundamental pathway of cofactor-free oxidases and oxygenases-catalyzed reactions involving the triplet O₂ activation but also discloses a new pathway that undergoes the processes of PCET and open-shell singlet transition state.

NHC-Catalyzed Transformation Reactions of Imines: Electrophilic versus Nucleophilic Attack

It is generally accepted that N-heterocyclic carbene (NHC)-catalyzed imine transformations are initiated by a nucleophilic attack (NA) by NHC. However, due to significant nucleophilicity of the iminyl nitrogen atom in imines, the electrophilic attack (EA) by electrophiles onto imine would also be a possible mechanism of these kinds of reactions. Therefore, we use the quantum mechanical approach to disclose that both the NA and EA modes could be switchable for a wide range of NHC-catalyzed transformations of imines.

Cryptic Oxidative Transamination of Hydroxynaphthoquinone in Natural Product Biosynthesis

Pyridoxal 5'-phosphate (PLP)-dependent enzymes are a group of versatile enzymes that catalyze various reactions, but only a small number of them react with O₂. Here, we report an unprecedented PLP-dependent enzyme, NphE, that catalyzes both transamination and two-electron oxidation using O₂ as an oxidant. Our intensive analysis reveals that NphE transfers the l-glutamate-derived amine to 1,3,6,8-tetrahydroxynaphthalene-derived mompain to form 8-amino-flaviolin (8-AF) via a highly conjugated quinonoid intermediate that is reactive with O₂. During the NphE reaction, O₂ is reduced to yield H₂O₂. An integrated technique involving NphE structure prediction by AlphaFold v2.0 and molecular dynamics simulation suggested the O₂-accessible cavity. Our in vivo results demonstrated that 8-AF is a genuine biosynthetic intermediate for the 1,3,6,8-tetrahydroxynaphthalene-derived meroterpenoid naphterpin without an amino group, which was supported by site-directed mutagenesis. This study clearly establishes the NphE reaction product 8-AF as a common intermediate with a cryptic amino group for the biosynthesis of terpenoid-polyketide hybrid natural products.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol^-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol^-1 can be inferred for the known k_cat ≈ 1 s^-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (k_cat = 8.6 × 10^-5 s^-1 to 2.7 s^-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH^⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.

Bienzyme-Catalytic and Dioxygenation-Mediated Anthraquinone Ring Opening

The C-10-C-4a bond cleavage of anthraquinone is believed to be a crucial step in fungal seco-anthraquinone biosynthesis and has long been proposed as a classic Baeyer-Villiger oxidation. Nonetheless, genetic, enzymatic, and chemical information on ring opening remains elusive. Here, a revised questin ring-opening mechanism was elucidated by in vivo gene disruption, in vitro enzymatic analysis, and 18O chasing experiments. It has been confirmed that the reductase GedF is responsible for the reduction of the keto group at C-10 in questin to a hydroxyl group with the aid of NADPH. The C-10-C-4a bond of the resultant questin hydroquinone is subsequently cleaved by the atypical cofactor-free dioxygenase GedK, giving rise to desmethylsulochrin. This proposed bienzyme-catalytic and dioxygenation-mediated anthraquinone ring-opening reaction shows universality.

Dioxygen Activation and Pyrrole α-Cleavage with Calix[4]pyrrolato Aluminates: Enzyme Model by Structural Constraint

The present work describes the reaction of triplet dioxygen with the porphyrinogenic calix[4]pyrrolato aluminates to alkylperoxido aluminates in high selectivity. Multiconfigurational quantum chemical computations disclose the mechanism for this spin‐forbidden process. Despite a negligible spin–orbit coupling constant, the intersystem crossing (ISC) is facilitated by singlet and triplet state degeneracy and spin–vibronic coupling. The formed peroxides are stable toward external substrates but undergo an unprecedented oxidative pyrrole α‐cleavage by ligand aromatization/dearomatization‐initiated O−O σ‐bond scission. A detailed comparison of the calix[4]pyrrolato aluminates with dioxygen‐related enzymology provides insights into the ISC of metal‐ or cofactor‐free enzymes. It substantiates the importance of structural constraint and element–ligand cooperativity for the functions of aerobic life.

Theoretical Study on the Mechanism of the Acylate Reaction of β-Lactamase

Using density functional theory and a cluster approach, we study the reaction potential surface and compute Gibbs free energies for the acylate reaction of β-lactamase with penicillin G, where the solvent effect is important and taken into consideration. Two reaction paths are investigated: one is a multi-step process with a rate-limit energy barrier of 19.1 kcal/mol, which is relatively small, and the reaction can easily occur; the other is a one-step process with a barrier of 45.0 kcal/mol, which is large and thus makes the reaction hard to occur. The reason why the two paths have different barriers is explained.

Cofactor-free oxidase-mimetic nanomaterials from self-assembled histidine-rich peptides

Natural oxidases mainly rely on cofactors and well-arranged amino acid residues for catalysing electron-transfer reactions but suffer from non-recovery of their activity upon externally induced protein unfolding. However, it remains unknown whether residues at the active site can catalyse similar reactions in the absence of the cofactor. Here, we describe a series of self-assembling, histidine-rich peptides, as short as a dipeptide, with catalytic function similar to that of haem-dependent peroxidases. The histidine residues of the peptide chains form periodic arrays that are able to catalyse H₂O₂ reduction reactions efficiently through the formation of reactive ternary complex intermediates. The supramolecular catalyst exhibiting the highest activity could be switched between inactive and active states without loss of activity for ten cycles of heating/cooling or acidification/neutralization treatments, demonstrating the reversible assembly/disassembly of the active residues. These findings may aid the design of advanced biomimetic catalytic materials and provide a model for primitive cofactor-free enzymes.

DpgC-Catalyzed Peroxidation of 3,5-Dihydroxyphenylacetyl-CoA (DPA-CoA): Insights into the Spin-Forbidden Transition and Charge Transfer Mechanisms*

Despite being a very strong oxidizing agent, most organic molecules are not oxidized in the presence of O₂ at room temperature because O₂ is a diradical whereas most organic molecules are closed-shell. Oxidation then requires a change in the spin state of the system, which is forbidden according to non-relativistic quantum theory. To overcome this limitation, oxygenases usually rely on metal or redox cofactors to catalyze the incorporation of, at least, one oxygen atom into an organic substrate. However, some oxygenases do not require any cofactor, and the detailed mechanism followed by these enzymes remains elusive. To fill this gap, here the mechanism for the enzymatic cofactor-independent oxidation of 3,5-dihydroxyphenylacetyl-CoA (DPA-CoA) is studied by combining multireference calculations on a model system with QM/MM calculations. Our results reveal that intersystem crossing takes place without requiring the previous protonation of molecular oxygen. The characterization of the electronic states reveals that electron transfer is concomitant with the triplet-singlet transition. The enzyme plays a passive role in promoting the intersystem crossing, although spontaneous reorganization of the water wire connecting the active site with the bulk presets the substrate for subsequent chemical transformations. The results show that the stabilization of the singlet radical-pair between dioxygen and enolate is enough to promote spin-forbidden reaction without the need for neither metal cofactors nor basic residues in the active site.

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

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

Transcriptomic analysis of polyaromatic hydrocarbon degradation by the halophilic fungus Aspergillus sydowii at hypersaline conditions

Polycyclic aromatic hydrocarbons (PAHs) are among the most persistent xenobiotic compounds, with high toxicity effects. Mycoremediation with halophilic Aspergillus sydowii was used for their removal from a hypersaline medium (1 M NaCl). A. sydowii metabolized PAHs as sole carbon sources, resulting in the removal of up to 90% for both PAHs [benzo [a] pyrene (BaP) and phenanthrene (Phe)] after 10 days. Elimination of Phe and BaP was almost exclusively due to biotransformation and not adsorption by dead mycelium and did not correlate with the activity of lignin modifying enzymes (LME). Transcriptomes of A. sydowii grown on PAHs, or on glucose as control, both at hypersaline conditions, revealed 170 upregulated and 76 downregulated genes. Upregulated genes were related to starvation, cell wall remodelling, degradation and metabolism of xenobiotics, DNA/RNA metabolism, energy generation, signalling and general stress responses. Changes of LME expression levels were not detected, while the chloroperoxidase gene, possibly related to detoxification processes in fungi, was strongly upregulated. We propose that two parallel metabolic pathways (mitochondrial and cytosolic) are involved in degradation and detoxification of PAHs in A. sydowii resulting in intracellular oxidation of PAHs. To the best of our knowledge, this is the most comprehensive transcriptomic analysis on fungal degradation of PAHs.

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.

An α/β-Hydrolase Fold Subfamily Comprising Pseudomonas Quinolone Signal-Cleaving Dioxygenases

The quinolone ring is a common core structure of natural products exhibiting antimicrobial, cytotoxic, and signaling activities. A prominent example is the Pseudomonas quinolone signal (PQS), a quorum-sensing signal molecule involved in the regulation of virulence of Pseudomonas aeruginosa The key reaction to quinolone inactivation and biodegradation is the cleavage of the 3-hydroxy-4(1H)-quinolone ring, catalyzed by dioxygenases (HQDs), which are members of the α/β-hydrolase fold superfamily. The α/β-hydrolase fold core domain consists of a β-sheet surrounded by α-helices, with an active site usually containing a catalytic triad comprising a nucleophilic residue, an acidic residue, and a histidine. The nucleophile is located at the tip of a sharp turn, called the "nucleophilic elbow." In this work, we developed a search workflow for the identification of HQD proteins from databases. Search and validation criteria include an [H-x(2)-W] motif at the nucleophilic elbow, an [HFP-x(4)-P] motif comprising the catalytic histidine, the presence of a helical cap domain, the positioning of the triad's acidic residue at the end of β-strand 6, and a set of conserved hydrophobic residues contributing to the substrate cavity. The 161 candidate proteins identified from the UniProtKB database originate from environmental and plant-associated microorganisms from all domains of life. Verification and characterization of HQD activity of 9 new candidate proteins confirmed the reliability of the search strategy and suggested residues correlating with distinct substrate preferences. Among the new HQDs, PQS dioxygenases from Nocardia farcinica, N. cyriacigeorgica, and Streptomyces bingchenggensis likely are part of a catabolic pathway for alkylquinolone utilization.IMPORTANCE Functional annotation of protein sequences is a major requirement for the investigation of metabolic pathways and the identification of sought-after biocatalysts. To identify heterocyclic ring-cleaving dioxygenases within the huge superfamily of α/β-hydrolase fold proteins, we defined search and validation criteria for the primarily motif-based identification of 3-hydroxy-4(1H)-quinolone 2,4-dioxygenases (HQD). HQDs are key enzymes for the inactivation of metabolites, which can have signaling, antimicrobial, or cytotoxic functions. The HQD candidates detected in this study occur particularly in environmental and plant-associated microorganisms. Because HQDs active toward the Pseudomonas quinolone signal (PQS) likely contribute to interactions within microbial communities and modulate the virulence of Pseudomonas aeruginosa, we analyzed the catalytic properties of a PQS-cleaving subset of HQDs and specified characteristics to identify PQS-cleaving dioxygenases within the HQD family.

Constructing Spin-Adiabatic States for the Modeling of Spin-Crossing Reactions. I. A Shared-Orbital Implementation

In the modeling of spin-crossing reactions, it has become popular to directly explore the spin-adiabatic surfaces. Specifically, through constructing spin-adiabatic states from a two-state Hamiltonian (with spin-orbit coupling matrix elements) at each geometry, one can readily employ advanced geometry optimization algorithms to acquire a "transition state" structure, where the spin crossing occurs. In this work, we report the implementation of a fully-variational spin-adiabatic approach based on Kohn-Sham density functional theory spin states (sharing the same set of molecular orbitals) and the Breit-Pauli one-electron spin-orbit operator. For three model spin-crossing reactions [predissociation of N₂O, singlet-triplet conversion in CH₂, and CO addition to Fe(CO)₄], the spin-crossing points were obtained. Our results also indicated the Breit-Pauli one-electron spin-orbit coupling can vary significantly along the reaction pathway on the spin-adiabatic energy surface. On the other hand, due to the restriction that low-spin and high-spin states share the same set of molecular orbitals, the acquired spin-adiabatic energy surface shows a cusp (i.e. a first-order discontinuity) at the crossing point, which prevents the use of standard geometry optimization algorithms to pinpoint the crossing point. An extension with this restriction removed is being developed to achieve the smoothness of spin-adiabatic surfaces.

Monooxygenation of aromatic compounds by flavin-dependent monooxygenases

Many flavoenzymes catalyze hydroxylation of aromatic compounds especially phenolic compounds have been isolated and characterized. These enzymes can be classified as either single-component or two-component flavin-dependent hydroxylases (monooxygenases). The hydroxylation reactions catalyzed by the enzymes in this group are useful for modifying the biological properties of phenolic compounds. This review aims to provide an in-depth discussion of the current mechanistic understanding of representative flavin-dependent monooxygenases including 3-hydroxy-benzoate 4-hydroxylase (PHBH, a single-component hydroxylase), 3-hydroxyphenylacetate 4-hydroxylase (HPAH, a two-component hydroxylase), and other monooxygenases which catalyze reactions in addition to hydroxylation, including 2-methyl-3-hydroxypyridine-5-carboxylate oxygenase (MHPCO, a single-component enzyme that catalyzes aromatic-ring cleavage), and HadA monooxygenase (a two-component enzyme that catalyzes additional group elimination reaction). These enzymes have different unique structural features which dictate their reactivity toward various substrates and influence their ability to stabilize flavin intermediates such as C4a-hydroperoxyflavin. Understanding the key catalytic residues and the active site environments important for governing enzyme reactivity will undoubtedly facilitate future work in enzyme engineering or enzyme redesign for the development of biocatalytic methods for the synthesis of valuable compounds.

How Chemical Environment Activates Anthralin and Molecular Oxygen for Direct Reaction

The role of the chemical environment in promoting anthralin/O₂ reactions was discovered using neat solvents to model the amino acids of a cofactor-independent oxygenase. Experimental and computational results highlight the importance of the substrate-enolate, which is accessed via energetically small, escalating steps in which the ground-state keto-isomer is tautomerized to an enol and then ionized by solvent. The resulting ion-pair is poised for spontaneous electron transfer to O₂. Similar activation may be exploited in biological/nonbiological oxidations involving O₂.

Hydroxyl Radical-Coupled Electron-Transfer Mechanism of Flavin-Dependent Hydroxylases

Class A flavin-dependent hydroxylases (FdHs) catalyze the hydroxylation of organic compounds in a site- and stereoselective manner. In stark contrast, conventional synthetic routes require environmentally hazardous reagents and give modest yields. Thus, understanding the detailed mechanism of this class of enzymes is essential to their rational manipulation for applications in green chemistry and pharmaceutical production. Both electrophilic substitution and radical intermediate mechanisms have been proposed as interpretations of FdH hydroxylation rates and optical spectra. While radical mechanistic steps are often difficult to examine directly, modern quantum chemistry calculations combined with statistical mechanical approaches can yield detailed mechanistic models providing insights that can be used to differentiate reaction pathways. In the current work, we report quantum mechanical/molecular mechanical (QM/MM) calculations on the fungal TropB enzyme that shows an alternative reaction pathway in which hydroxylation through a hydroxyl radical-coupled electron-transfer mechanism is significantly favored over electrophilic substitution. Furthermore, QM/MM calculations on several modified flavins provide a more consistent interpretation of the experimental trends in the reaction rates seen experimentally for a related enzyme, para-hydroxybenzoate hydroxylase. These calculations should guide future enzyme and substrate design strategies and broaden the scope of biological spin chemistry.

Bridging the Homogeneous-Heterogeneous Divide: Modeling Spin for Reactivity in Single Atom Catalysis

Single atom catalysts (SACs) are emergent catalytic materials that have the promise of merging the scalability of heterogeneous catalysts with the high activity and atom economy of homogeneous catalysts. Computational, first-principles modeling can provide essential insight into SAC mechanism and active site configuration, where the sub-nm-scale environment can challenge even the highest-resolution experimental spectroscopic techniques. Nevertheless, the very properties that make SACs attractive in catalysis, such as localized d electrons of the isolated transition metal center, make them challenging to study with conventional computational modeling using density functional theory (DFT). For example, Fe/N-doped graphitic SACs have exhibited spin-state dependent reactivity that remains poorly understood. However, spin-state ordering in DFT is very sensitive to the nature of the functional approximation chosen. In this work, we develop accurate benchmarks from correlated wavefunction theory (WFT) for relevant octahedral complexes. We use those benchmarks to evaluate optimal DFT functional choice for predicting spin state ordering in small octahedral complexes as well as models of pyridinic and pyrrolic nitrogen environments expected in larger SACs. Using these guidelines, we determine Fe/N-doped graphene SAC model properties and reactivity as well as their sensitivities to DFT functional choice. Finally, we conclude with broad recommendations for computational modeling of open-shell transition metal single-atom catalysts.

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.

Mechanism of Oxygenase-Pathway Reactions Catalyzed by Rubisco from Large-Scale Kohn-Sham Density Functional Calculations

Ribulose 1,5-bisphosphate carboxylase-oxygenase (Rubisco) is the primary carbon-fixing enzyme in photosynthesis, fixing CO₂ to a 5-carbon sugar, RuBP, in a series of five reactions. However, it also catalyzes an oxygenase reaction by O₂ addition to the same enolized RuBP substrate in an analogous reaction series in the same active site, producing a waste product and loss of photosynthetic efficiency. Starting from RuBP, the reactions are enolization to the enediolate form, addition of CO₂ or O₂ to form the carboxy or peroxo adduct, hydration to form a gemdiolate, scission of the C2-C3 bond of the original RuBP, and stereospecific or nonstereospecific protonation to form two molecules of the 3-carbon PGA product, or one molecule of PGA, one of 2-carbon PG (waste product), and one water molecule. Reducing the loss of efficiency from the oxygenase reaction is an attractive means to increase crop productivity. However, lack of understanding of key aspects of the catalytic mechanisms for both the carboxylase and oxygenase reactions, particularly those involving proton exchanges and roles of water molecules, has stymied efforts at re-engineering Rubisco to reduce losses from the oxygenation reaction. As the stable form of molecular oxygen is the triplet biradical state (³O₂), its reaction with near-universal singlet-state molecules is formally spin forbidden. Although in oxygenase enzymes, ³O₂ activation is usually achieved by one-electron transfers using transition-metal ions or organic cofactors, recently, cofactor-less oxygenases in which the substrate itself is the source of the electron for ³O₂ activation have been identified, but in all such cases an aromatic ring stabilizes the substrate's negative charge. Here we present the first large-scale Kohn-Sham density functional theory study of the reaction mechanism of the Rubisco oxygenase pathway. First, we show that the enediolate substrate complexed to Mg²⁺ and its ligands extends the region for charge delocalization and stabilization of its negative charge to allow formation of a caged biradical enediolate-O₂ complex. Thus, Rubisco is a unique type of oxygenase without precedent in the literature. Second, for the O₂ addition to proceed to the singlet peroxo-adduct intermediate, the system must undergo an intersystem crossing. We found that the presence of protonated LYS334 is required to stabilize this intermediate and that both factors (strongly stabilized anion and protonated LYS334) facilitate a barrier-less activation of ³O₂. This finding supports our recent proposal that deoxygenation, that is, reversal of gas binding, is possible. Third, as neither CO₂ nor O₂ binds to the enzyme, our findings support the proposal from our recent carboxylase study that the observed K_C or K_O (Michaelis-Menten constants) in the steady-state kinetics reflect the respective adducts, carboxy or peroxo. Fourth, after computing hydration pathways with water addition both syn and anti to C3, we found, in contrast to the results of our carboxylation study indicating anti addition, that in the oxygenation reaction only syn-hydration is capable of producing a stable gemdiolate that facilitates the rate-limiting C2-C3 bond scission to final products. Fifth, we propose that an excess proton we previously found was required in the carboxylation reaction for activating the C2-C3 bond scission is utilized in the oxygenation reaction for the required elimination of a water molecule. In summary, despite its oxygenase handicap, Rubisco's success in directing 75% of its substrate through the carboxylation pathway can be considered impressively effective. Although native C3 Rubiscos are in a fix with unwanted activity of ³O₂ hampering its primary carboxylase function, mechanistic differences presented here with findings in our recent carboxylase study for both the gas-addition and subsequent reactions provide some clues as to how creative Rubisco re-engineering may offer a solution to reducing the oxygenase activity.

3-Hydroxyflavones and 3-Hydroxy-4-oxoquinolines as Carbon Monoxide-Releasing Molecules

Carbon monoxide-releasing molecules (CORMs) that enable the delivery of controlled amounts of CO are of strong current interest for applications in biological systems. In this review, we examine the various conditions under which CO is released from 3-hydroxyflavones and 3-hydroxy-4-oxoquinolines to advance the understanding of how these molecules, or derivatives thereof, may be developed as CORMs. Enzymatic pathways from quercetin dioxygenases and 3-hydroxy-4-oxoquinoline dioxygenases leading to CO release are examined, along with model systems for these enzymes. Base-catalyzed and non-redox-metal promoted CO release, as well as UV and visible light-driven CO release from 3-hydroxyflavones and 3-hydroxy-4-oxoquinolines, are summarized. The visible light-induced CO release reactivity of recently developed extended 3-hydroxyflavones and a 3-hydroxybenzo[g]quinolone, and their uses as intracellular CORMs, are discussed. Overall, this review provides insight into the chemical factors that affect the thermal and photochemical dioxygenase-type CO release reactions of these heterocyclic compounds.

Mechanisms of Cytochrome P450-Catalyzed Oxidations

Enzymes are complex biological catalysts and are critical to life. Most oxidations of chemicals are catalyzed by cytochrome P450 (P450, CYP) enzymes, which generally utilize mixed-function oxidase stoichiometry, utilizing pyridine nucleotides as electron donors: NAD(P)H + O₂ + R → NAD(P)⁺ + RO + H₂O (where R is a carbon substrate and RO is an oxidized product). The catalysis of oxidations is largely understood in the context of the heme iron-oxygen complex generally referred to as Compound I, formally FeO³⁺, whose basis was in peroxidase chemistry. Many X-ray crystal structures of P450s are now available (≥ 822 structures from ≥146 different P450s) and have helped in understanding catalytic specificity. In addition to hydroxylations, P450s catalyze more complex oxidations, including C-C bond formation and cleavage. Enzymes derived from P450s by directed evolution can even catalyze more unusual reactions, e.g. cyclopropanation. Current P450 questions under investigation include the potential role of the intermediate Compound 0 (formally Fe^III-O₂ ^-) in catalysis of some reactions, the roles of high- and low-spin forms of Compound I, the mechanism of desaturation, the roles of open and closed structures of P450s in catalysis, the extent of processivity in multi-step oxidations, and the role of the accessory protein cytochrome b ₅. More global questions include exactly how structure drives function, prediction of catalysis, and roles of multiple protein conformations.

Does Substrate Positioning Affect the Selectivity and Reactivity in the Hectochlorin Biosynthesis Halogenase?

In this work we present the first computational study on the hectochlorin biosynthesis enzyme HctB, which is a unique three-domain halogenase that activates non-amino acid moieties tethered to an acyl-carrier, and as such may have biotechnological relevance beyond other halogenases. We use a combination of small cluster models and full enzyme structures calculated with quantum mechanics/molecular mechanics methods. Our work reveals that the reaction is initiated with a rate-determining hydrogen atom abstraction from substrate by an iron (IV)-oxo species, which creates an iron (III)-hydroxo intermediate. In a subsequent step the reaction can bifurcate to either halogenation or hydroxylation of substrate, but substrate binding and positioning drives the reaction to optimal substrate halogenation. Furthermore, several key residues in the protein have been identified for their involvement in charge-dipole interactions and induced electric field effects. In particular, two charged second coordination sphere amino acid residues (Glu223 and Arg245) appear to influence the charge density on the Cl ligand and push the mechanism toward halogenation. Our studies, therefore, conclude that nonheme iron halogenases have a chemical structure that induces an electric field on the active site that affects the halide and iron charge distributions and enable efficient halogenation. As such, HctB is intricately designed for a substrate halogenation and operates distinctly different from other nonheme iron halogenases.

A Comparative Review on the Catalytic Mechanism of Nonheme Iron Hydroxylases and Halogenases

Enzymatic halogenation and haloperoxidation are unusual processes in biology; however, a range of halogenases and haloperoxidases exist that are able to transfer an aliphatic or aromatic C–H bond into C–Cl/C–Br. Haloperoxidases utilize hydrogen peroxide, and in a reaction with halides (Cl−/Br−), they react to form hypohalides (OCl−/OBr−) that subsequently react with substrate by halide transfer. There are three types of haloperoxidases, namely the iron-heme, nonheme vanadium, and flavin-dependent haloperoxidases that are reviewed here. In addition, there are the nonheme iron halogenases that show structural and functional similarity to the nonheme iron hydroxylases and form an iron(IV)-oxo active species from a reaction of molecular oxygen with α-ketoglutarate on an iron(II) center. They subsequently transfer a halide (Cl−/Br−) to an aliphatic C–H bond. We review the mechanism and function of nonheme iron halogenases and hydroxylases and show recent computational modelling studies of our group on the hectochlorin biosynthesis enzyme and prolyl-4-hydroxylase as examples of nonheme iron halogenases and hydroxylases. These studies have established the catalytic mechanism of these enzymes and show the importance of substrate and oxidant positioning on the stereo-, chemo- and regioselectivity of the reaction that takes place.

Quantum Mechanics/Molecular Mechanics Studies on the Relative Reactivities of Compound I and II in Cytochrome P450 Enzymes

The cytochromes P450 are drug metabolizing enzymes in the body that typically react with substrates through a monoxygenation reaction. During the catalytic cycle two reduction and protonation steps generate a high-valent iron (IV)-oxo heme cation radical species called Compound I. However, with sufficient reduction equivalents present, the catalytic cycle should be able to continue to the reduced species of Compound I, called Compound II, rather than a reaction of Compound I with substrate. In particular, since electron transfer is usually on faster timescales than atom transfer, we considered this process feasible and decided to investigate the reaction computationally. In this work we present a computational study using density functional theory methods on active site model complexes alongside quantum mechanics/molecular mechanics calculations on full enzyme structures of cytochrome P450 enzymes. Specifically, we focus on the relative reactivity of Compound I and II with a model substrate for O⁻H bond activation. We show that generally the barrier heights for hydrogen atom abstraction are higher in energy for Compound II than Compound I for O⁻H bond activation. Nevertheless, for the activation of such bonds, Compound II should still be an active oxidant under enzymatic conditions. As such, our computational modelling predicts that under high-reduction environments the cytochromes P450 can react with substrates via Compound II but the rates will be much slower.

Mechanistic Insight on the Activity and Substrate Selectivity of Nonheme Iron Dioxygenases

Nonheme iron dioxygenases catalyze vital reactions for human health particularly related to aging processes. They are involved in the biosynthesis of amino acids, but also the biodegradation of toxic compounds. Typically they react with their substrate(s) through oxygen atom transfer, although often with the assistance of a co-substrate like α-ketoglutarate that is converted to succinate and CO₂ . Many reaction processes catalyzed by the nonheme iron dioxygenases are stereoselective or regiospecific and hence understanding the mechanism and protein involvement in the selectivity is important for the design of biotechnological applications of these enzymes. To this end, I will review recent work of our group on nonheme iron dioxygenases and include background information on their general structure and catalytic cycle. Examples of stereoselective and regiospecific reaction mechanisms we elucidated are for the AlkB repair enzyme, prolyl-4-hydroxylase and the ergothioneine biosynthesis enzyme. Finally, I cover an example where we bioengineered S-p-hydroxymandelate synthase into the R-p-hydroxymandelate synthase.

Rubisco is not really so bad

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) is the most widespread carboxylating enzyme in autotrophic organisms. Its kinetic and structural properties have been intensively studied for more than half a century. Yet important aspects of the catalytic mechanism remain poorly understood, especially the oxygenase reaction. Because of its relatively modest turnover rate (a few catalytic events per second) and the competitive inhibition by oxygen, Rubisco is often viewed as an inefficient catalyst for CO₂ fixation. Considerable efforts have been devoted to improving its catalytic efficiency, so far without success. In this review, we re-examine Rubisco's catalytic performance by comparison with other chemically related enzymes. We find that Rubisco is not especially slow. Furthermore, considering both the nature and the complexity of the chemical reaction, its kinetic properties are unremarkable. Although not unique to Rubisco, oxygenation is not systematically observed in enolate and enamine forming enzymes and cannot be considered as an inevitable consequence of the mechanism. It is more likely the result of a compromise between chemical and metabolic imperatives. We argue that a better description of Rubisco mechanism is still required to better understand the link between CO₂ and O₂ reactivity and the rationale of Rubisco diversification and evolution.

Keto-Enol Tautomerization Triggers an Electrophilic Aldehyde Deformylation Reaction by a Nonheme Manganese(III)-Peroxo Complex

Oxygen atom transfer by high-valent enzymatic intermediates remains an enigma in chemical catalysis. In particular, manganese is an important first-row metal involved in key biochemical processes, including the biosynthesis of molecular oxygen (through the photosystem II complex) and biodegradation of toxic superoxide to hydrogen peroxide by superoxide dismutase. Biomimetic models of these biological systems have been developed to gain understanding on the structure and properties of short-lived intermediates but also with the aim to create environmentally benign oxidants. In this work, we report a combined spectroscopy, kinetics and computational study on aldehyde deformylation by two side-on manganese(III)-peroxo complexes with bispidine ligands. Both manganese(III)-peroxo complexes are characterized by UV-vis and mass spectrometry techniques, and their reactivity patterns with aldehydes was investigated. We find a novel mechanism for the reaction that is initiated by a hydrogen atom abstraction reaction, which enables a keto-enol tautomerization in the substrate. This is an essential step in the mechanism that makes an electrophilic attack on the olefin bond possible as the attack on the aldehyde carbonyl is too high in energy. Kinetics studies determine a large kinetic isotope effect for the replacement of the transferring hydrogen atom by deuterium, while replacing the transferring hydrogen atom by a methyl group makes the substrate inactive and hence confirm the hypothesized mechanism. Our new mechanism is confirmed with density functional theory modeling on the full mechanism and rationalized through valence bond and thermochemical cycles. Our unprecedented new mechanism may have relevance to biological and biomimetic chemistry processes in general and gives insight into the reactivity patterns of metal-peroxo and metal-hydroperoxo intermediates in general.

How Are Substrate Binding and Catalysis Affected by Mutating Glu127 and Arg161 in Prolyl-4-hydroxylase? A QM/MM and MD Study

Prolyl-4-hydroxylase is a vital enzyme for human physiology involved in the biosynthesis of 4-hydroxyproline, an essential component for collagen formation. The enzyme performs a unique stereo- and regioselective hydroxylation at the C4 position of proline despite the fact that the C5 hydrogen atoms should be thermodynamically easier to abstract. To gain insight into the mechanism and find the origin of this regioselectivity, we have done a quantum mechanics/molecular mechanics (QM/MM) study on wildtype and mutant structures. In a previous study (Timmins et al., 2017) we identified several active site residues critical for substrate binding and positioning. In particular, the Glu127 and Arg161 were shown to form multiple hydrogen bonding and ion-dipole interactions with substrate and could thereby affect the regio- and stereoselectivity of the reaction. In this work, we decided to test that hypothesis and report a QM/MM and molecular dynamics (MD) study on prolyl-4-hydroxylase and several active site mutants where Glu127 or Arg161 are mutated for Asp, Gln, or Lys. Thus, the R161D and R161Q mutants give very high barriers for hydrogen atom abstraction from any proline C-H bond and therefore will be inactive. The R161K mutant, by contrast, sees the regio- and stereoselectivity of the reaction change but still is expected to hydroxylate proline at room temperature. By contrast, the Glu127 mutants E127D and E127Q show possible changes in regioselectivity with the former being more probable to react compared to the latter.

Indole Biodegradation in Acinetobacter sp. Strain O153: Genetic and Biochemical Characterization

Indole is a molecule of considerable biochemical significance, acting as both an interspecies signal molecule and a building block of biological elements. Bacterial indole degradation has been demonstrated for a number of cases; however, very little is known about genes and proteins involved in this process. This study reports the cloning and initial functional characterization of genes (iif and ant cluster) responsible for indole biodegradation in Acinetobacter sp. strain O153. The catabolic cascade was reconstituted in vitro with recombinant proteins, and each protein was assigned an enzymatic function. Degradation starts with oxidation, mediated by the IifC and IifD flavin-dependent two-component oxygenase system. Formation of indigo is prevented by IifB, and the final product, anthranilic acid, is formed by IifA, an enzyme which is both structurally and functionally comparable to cofactor-independent oxygenases. Moreover, the iif cluster was identified in the genomes of a wide range of bacteria, suggesting the potential of widespread Iif-mediated indole degradation. This work provides novel insights into the genetic background of microbial indole biodegradation.IMPORTANCE The key finding of this research is identification of the genes responsible for microbial biodegradation of indole, a toxic N-heterocyclic compound. A large amount of indole is present in urban wastewater and sewage sludge, creating a demand for an efficient and eco-friendly means to eliminate this pollutant. A common strategy of oxidizing indole to indigo has the major drawback of producing insoluble material. Genes and proteins of Acinetobacter sp. strain O153 (DSM 103907) reported here pave the way for effective and indigo-free indole removal. In addition, this work suggests possible novel means of indole-mediated bacterial interactions and provides the basis for future research on indole metabolism.

Understanding How Prolyl-4-hydroxylase Structure Steers a Ferryl Oxidant toward Scission of a Strong C-H Bond

Prolyl-4-hydroxylase (P4H) is a non-heme iron hydroxylase that regio- and stereospecifically hydroxylates proline residues in a peptide chain into R-4-hydroxyproline, which is essential for collagen cross-linking purposes in the human body. Surprisingly, in P4H, a strong aliphatic C-H bond is activated, while thermodynamically much weaker aliphatic C-H groups, that is, at the C³ and C⁵ positions, are untouched. Little is known on the origins of the high regio- and stereoselectivity of P4H and many non-heme and heme enzymes in general, and insight into this matter may be relevant to Biotechnology as well as Drug Development. The active site of the protein contains two aromatic residues (Tyr₁₄₀ and Trp₂₄₃) that we expected to be crucial for guiding the regioselectivity of the reaction. We performed a detailed quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) study on wild-type and mutant structures. The work shows that Trp₂₄₃ is involved in key protein loop-loop interactions that affect the shape and size of the substrate binding pocket and its mutation has major long-range effects. By contrast, the Tyr₁₄₀ residue is shown to guide the regio- and stereoselectivity by holding the substrate and ferryl oxidant in a specific orientation through hydrogen bonding and π-stacking interactions. Compelling evidence is found that the Tyr₁₄₀ residue is involved in expelling the product from the binding pocket after the reaction is complete. It is shown that mutations where the hydrogen bonding network that involves the Tyr₁₄₀ and Trp₂₄₃ residues is disrupted lead to major changes in folding of the protein and the size and shape of the substrate binding pocket. Specifically, the Trp₂₄₃ residue positions the amino acid side chains of Arg₁₆₁ and Glu₁₂₇ in specific orientations with substrate. As such, the P4H enzyme is a carefully designed protein with a subtle and rigid secondary structure that enables the binding of substrate, guides the regioselectivity, and expels product efficiently.