Why Is the Oxidation State of Iron Crucial for the Activity of Heme-Dependent Aldoxime Dehydratase? A QM/MM Study

Bioinspired Iron-Catalyzed Dehydration of Aldoximes to Nitriles: A General N-O Redox-Cleavage Method

Inspired by OxdA that operates biocatalytic aldoxime dehydration, we have developed an efficient iron catalyst, Cp*Fe(1,2-Cy₂PC₆H₄O) (1), which rapidly converts various aliphatic and aromatic aldoximes to nitriles with release of H₂O at room temperature. The catalysis involves redox activation of the N-O bond by a 1e^- transfer from the iron catalyst to the oxime. Such redox-mediated N-O cleavage was demonstrated by the isolation of a ferrous iminato intermediate from the reaction of the ketoxime substrate. This iron-catalyzed acceptorless dehydration approach represents a general method for the preparation of nitriles, and it also delivers salicylonitriles by catalyzing the Kemp elimination reaction.

Systematic QM Region Construction in QM/MM Calculations Based on Uncertainty Quantification

While QM/MM studies of enzymatic reactions are widely used in computational chemistry, the results of such studies are subject to numerous sources of uncertainty, and the effect of different choices by the simulation scientist that are required when setting up QM/MM calculations is often unclear. In particular, the selection of the QM region is crucial for obtaining accurate and reliable results. Simply including amino acids by their distance to the active site is mostly not sufficient as necessary residues are missing or unimportant residues are included without evidence. Here, we take a first step toward quantifying uncertainties in QM/MM calculations by assessing the sensitivity of QM/MM reaction energies with respect to variations of the MM point charges. We show that such a point charge variation analysis (PCVA) can be employed to judge the accuracy of QM/MM reaction energies obtained with a selected QM region and devise a protocol to systematically construct QM regions that minimize this uncertainty. We apply such a PCVA to the example of catechol O-methyltransferase and demonstrate that it provides a simple and reliable approach for the construction of the QM region. Our PCVA-based scheme is computationally efficient and requires only calculations for a system with a minimal QM region. Our work highlights the promise of applying methods of uncertainty quantification in computational chemistry.

Bifunctional Naphthoquinone Aromatic Amide-Oxime Derivatives Exert Combined Immunotherapeutic and Antitumor Effects through Simultaneous Targeting of Indoleamine-2,3-dioxygenase and Signal Transducer and Activator of Transcription 3

Indoleamine-2,3-dioxygenase 1 (IDO1) and signal transducer and activator of transcription 3 (STAT3) are important targets in the tumor microenvironment for cancer therapy. In the present study, a set of naphthoquinone aromatic amide-oxime derivatives were designed, which stimulated the immune response via IDO1 inhibition and simultaneously displayed powerful antitumor activity against three selected cancer cell lines through suppressing STAT3 signaling. The representative compound 8u bound effectively to IDO1, with greater inhibitory activity relative to the commercial IDO1 inhibitor 4-amino-N-(3-chloro-4-fluorophenyl)-N'-hydroxy-1,2,5-oxadiazole-3-carboximidamide (IDO5L) in addition to the efficient suppression of nuclear translocation of STAT3. Consistently, in vivo assays demonstrated a higher antiproliferative activity of compound 8u in both wild-type B16-F10 isograft tumors and an athymic HepG2 xenograft model relative to 1-methyl-l-tryptophan (1-MT) and doxorubicin (DOX). This bifunctional compound with dual immunotherapeutic and anticancer efficacy may represent a new generation of highly efficacious drug candidates for cancer therapy.

Biocatalytic production of adiponitrile and related aliphatic linear α,ω-dinitriles

Linear α,ω-dinitriles are important precursors for the polymer industry. Most prominently, adiponitrile is produced on an annual scale of ca. 1 million tons. However, a drawback of today’s dominating process is the need for large amounts of highly toxic hydrogen cyanide. In this contribution, an alternative approach towards such linear dinitriles is presented based on dehydration of readily available α,ω-dialdoximes at ambient conditions by means of aldoxime dehydratases. In contrast to existing production routes this biocatalytic route enables a highly regio- and chemoselective approach towards dinitriles without the use of hydrogen cyanide or harsh reaction conditions. In addition, a selective synthesis of adiponitrile with substrate loadings of up to 100 g/L and high yields of up to 80% was achieved. Furthermore, a lab scale process on liter scale leading to > 99% conversion at 50 g/L underlines the potential and robustness of this method for technical applicability.

Typically, preparation of the polymer precursors α,ω-dinitriles requires hydrogen cyanide. Here, the authors use aldoxime hydratase to produce adiponitrile and related aliphatic linear dinitriles under ambient conditions starting from readily available substrates without needing hydrogen cyanide.

Biocatalytic Synthesis of Nitriles through Dehydration of Aldoximes: The Substrate Scope of Aldoxime Dehydratases

Nitriles, which are mostly needed and produced by the chemical industry, play a major role in various industry segments, ranging from high-volume, low-price sectors, such as polymers, to low-volume, high-price sectors, such as chiral pharma drugs. A common industrial technology for nitrile production is ammoxidation as a gas-phase reaction at high temperature. Further popular approaches are substitution or addition reactions with hydrogen cyanide or derivatives thereof. A major drawback, however, is the very high toxicity of cyanide. Recently, as a synthetic alternative, a novel enzymatic approach towards nitriles has been developed with aldoxime dehydratases, which are capable of converting an aldoxime in one step through dehydration into nitriles. Because the aldoxime substrates are easily accessible, this route is of high interest for synthetic purposes. However, whenever a novel method is developed for organic synthesis, it raises the question of substrate scope as one of the key criteria for application as a "synthetic platform technology". Thus, the scope of this review is to give an overview of the current state of the substrate scope of this enzymatic method for synthesizing nitriles with aldoxime dehydratases. As a recently emerging enzyme class, a range of substrates has already been studied so far, comprising nonchiral and chiral aldoximes. This enzyme class of aldoxime dehydratases shows a broad substrate tolerance and accepts aliphatic and aromatic aldoximes, as well as arylaliphatic aldoximes. Furthermore, aldoximes with a stereogenic center are also recognized and high enantioselectivities are found for 2-arylpropylaldoximes, in particular. It is further noteworthy that the enantiopreference depends on the E and Z isomers. Thus, opposite enantiomers are accessible from the same racemic aldehyde and the same enzyme.

New function of aldoxime dehydratase: Redox catalysis and the formation of an unexpected product

In general, hemoproteins are capable of catalyzing redox reactions. Aldoxime dehydratase (OxdA), which is a unique heme-containing enzyme, catalyzes the dehydration of aldoximes to the corresponding nitriles. Its reaction is a rare example of heme directly activating an organic substrate, unlike the utilization of H₂O₂ or O₂ as a mediator of catalysis by other heme-containing enzymes. While it is unknown whether OxdA catalyzes redox reactions or not, we here for the first time detected catalase activity (which is one of the redox activities) of wild-type OxdA, OxdA(WT). Furthermore, we constructed a His320 → Asp mutant of OxdA [OxdA(H320D)], and found it exhibits catalase activity. Determination of the kinetic parameters of OxdA(WT) and OxdA(H320D) revealed that their K_m values for H₂O₂ were similar to each other, but the k_cat value of OxdA(H320D) was 30 times higher than that of OxdA(WT). Next, we examined another redox activity and found it was the peroxidase activity of OxdAs. While both OxdA(WT) and OxdA(H320D) showed the activity, the activity of OxdA(H320D) was dozens of times higher than that of OxdA(WT). These findings demonstrated that the H320D mutation enhances the peroxidase activity of OxdA. OxdAs (WT and H320D) were found to catalyze another redox reaction, a peroxygenase reaction. During this reaction of OxdA(H320D) with 1-methoxynaphthalene as a substrate, surprisingly, the reaction mixture changed to a color different from that with OxdA(WT), which was due to the known product, Russig’s blue. We purified and identified the new product as 1-methoxy-2-naphthalenol, which has never been reported as a product of the peroxygenase reaction, to the best of our knowledge. These findings indicated that the H320D mutation not only enhanced redox activities, but also significantly altered the hydroxylation site of the substrate.

Kemp Elimination Catalyzed by Naturally Occurring Aldoxime Dehydratases

Recently, the Kemp elimination reaction has been extensively studied in computational enzyme design of new catalysts, as no natural enzyme has evolved to catalyze this reaction. In contrast to in silico enzyme design, we were interested in searching for Kemp eliminase activity in natural enzymes with catalytic promiscuity. Based on similarities of substrate structures and reaction mechanisms, we assumed that the active sites of naturally abundant aldoxime dehydratases have the potential to catalyze the non-natural Kemp elimination reaction. We found several aldoxime dehydratases that are efficient catalysts of this reaction. Although a few natural enzymes have been identified with promiscuous Kemp eliminase activity, to the best of our knowledge, this is a rare example of Kemp elimination catalyzed by naturally occurring enzymes with high catalytic efficiency.

Cleavage mechanism of the aliphatic C–C bond catalyzed by 2,4′-dihydroxyacetophenone dioxygenase from Alcaligenes sp. 4HAP: a QM/MM study

Calculations suggest that the reactant complex may firstly undergo a triplet–quintet crossing to initiate the reaction and then the subsequent chemistry occurs on the multiple-states surfaces. The key C–C bond cleavage is accompanied by an insertion reaction of oxygen radical.

Toluene Dioxygenase-Catalysed Oxidation of Benzyl Azide to Benzonitrile: Mechanistic Insights for an Unprecedented Enzymatic Transformation

Enzymatic dioxygenation of benzyl azide by toluene dioxygenase (TDO) produces significant amounts of the cis-cyclohexadienediol derived from benzonitrile, along with the expected azido diols. We demonstrate that TDO catalyses the oxidation of benzyl azide to benzonitrile, which is further dioxygenated to produce the observed cis-diol. A proposed mechanism for this transformation involves initial benzylic monooxygenation followed by a nitrene-mediated rearrangement to form an oxime, which is further dehydrated to afford the nitrile. To the best of our knowledge, this is the first report of enzymatic oxidation of an alkyl azide to a nitrile. In addition, the described oxime-dehydration activity has not been reported for Rieske dioxygenases.

Quantum mechanics/molecular mechanics study of oxygen binding in hemocyanin

We report a combined quantum mechanics/molecular mechanics (QM/MM) study on the mechanism of reversible dioxygen binding in the active site of hemocyanin (Hc). The QM region is treated by broken-symmetry density functional theory (DFT) with spin projection corrections. The X-ray structures of deoxygenated (deoxyHc) and oxygenated (oxyHc) hemocyanin are well reproduced by QM/MM geometry optimizations. The computed relative energies strongly depend on the chosen density functional. They are consistent with the available thermodynamic data for oxygen binding in hemocyanin and in synthetic model complexes when the BH&HLYP hybrid functional with 50% Hartree-Fock exchange is used. According to the QM(BH&HLYP)/MM results, the reaction proceeds stepwise with two sequential electron transfer (ET) processes in the triplet state followed by an intersystem crossing to the singlet product. The first ET step leads to a nonbridged superoxo CuB(II)-O2(•-) intermediate via a low-barrier transition state. The second ET step is even more facile and yields a side-on oxyHc complex with the characteristic Cu2O2 butterfly core, accompanied by triplet-singlet intersystem crossing. The computed barriers are very small so that the two ET processes are expected to very rapid and nearly simultaneous.

Convergence in the QM-only and QM/MM modeling of enzymatic reactions: A case study for acetylene hydratase

We report systematic quantum mechanics-only (QM-only) and QM/molecular mechanics (MM) calculations on an enzyme-catalyzed reaction to assess the convergence behavior of QM-only and QM/MM energies with respect to the size of the chosen QM region. The QM and MM parts are described by density functional theory (typically B3LYP/def2-SVP) and the CHARMM force field, respectively. Extending our previous work on acetylene hydratase with QM regions up to 157 atoms (Liao and Thiel, J. Chem. Theory Comput. 2012, 8, 3793), we performed QM/MM geometry optimizations with a QM region M4 composed of 408 atoms, as well as further QM/MM single-point calculations with even larger QM regions up to 657 atoms. A charge deletion analysis was conducted for the previously used QM/MM model (M3a, with a QM region of 157 atoms) to identify all MM residues with strong electrostatic contributions to the reaction energetics (typically more than 2 kcal/mol), which were then included in M4. QM/MM calculations with this large QM region M4 lead to the same overall mechanism as the previous QM/MM calculations with M3a, but there are some variations in the relative energies of the stationary points, with a mean absolute deviation (MAD) of 2.7 kcal/mol. The energies of the two relevant transition states are close to each other at all levels applied (typically within 2 kcal/mol), with the first (second) one being rate-limiting in the QM/MM calculations with M3a (M4). QM-only gas-phase calculations give a very similar energy profile for QM region M4 (MAD of 1.7 kcal/mol), contrary to the situation for M3a where we had previously found significant discrepancies between the QM-only and QM/MM results (MAD of 7.9 kcal/mol). Extension of the QM region beyond M4 up to M7 (657 atoms) leads to only rather small variations in the relative energies from single-point QM-only and QM/MM calculations (MAD typically about 1-2 kcal/mol). In the case of acetylene hydratase, a model with 408 QM atoms thus seems sufficient to achieve convergence in the computed relative energies to within 1-2 kcal/mol.