Proton-Electron Transfer to the Active Site Is Essential for the Reaction Mechanism of Soluble Δ9-Desaturase

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.

DFT mechanistic study of biomimetic diiron complex catalyzed dehydrogenation: Unexpected Fe(III)Fe(III)-1,1-μ-hydroperoxy active species for hydride abstraction

The diiron active site is pivotal in catalyzing transformations in both biological and chemical systems. Recently, a range of biomimetic diiron catalysts have been synthesized, drawing inspiration from the active architecture of soluble methane monooxygenase (sMMO). These catalysts have been successfully deployed for the dehydrogenation of indolines, marking a significant advancement in the field. Using density functional theory (DFT) calculations, we have identified a novel mechanistic pathway that governs the dehydrogenation of indolines catalyzed by a biomimetic diiron complex. Specifically, this reaction is facilitated by the transfer of a hybrid atom from the C1 position of the substrate to the distal oxygen atom of the Fe(III)Fe(III)-1,1-μ-hydroperoxy active species. This transfer serves as the rate-limiting step for the heterolytic cleavage of the OO bond, ultimately generating the substrate cation. The mechanism we propose aligns well with mechanistic investigations incorporating both kinetic isotope effect (KIE) measurements and evaluations of stereochemical selectivity. This research contributes to the broader scientific understanding of catalysis involving biomimetic diiron complexes and offers valuable insights into the catalytic behaviors of non-heme diiron metalloenzymes.

Direct and remote control of electronic structures and redox potentials in μ-oxo diferric complexes

Non-heme diiron enzymes activate O2 for the oxidation of substrates in the form of peroxo FeIII2 or high-valent FeIV2 intermediates. We have developed a dinucleating bis(tetradentate) ligand system that stabilizes peroxo and hydroperoxo FeIII2 complexes with terminal 6-methylpyridine donors, while the peroxo FeIII2 intermediate is reactive with terminal pyridine donors presumably via conversion to a fluent high-valent FeIV2 intermediate. We present here a derivative with electron-donating methoxy substituents at the pyridine donors and its diferric complexes with an {FeIIIX(μ-O)FeIIIX} (X- = Cl-, OAc-, and OH-) or an {FeIII(μ-O)(μ-OAc)FeIII} core. The complex-induced oxidation of EtOH with H2O2 provides μ-OAc-, and in acetone, the complex with mixed OH-/OAc- exogenous donors is obtained. Both reactivities indicate a reactive fluent peroxo FeIII2 intermediate. The coupling constant J and the LMCT transitions are insensitive to the nature of the directly bound ligands X- and reflect mainly the electronic structure of the central {FeIII(μ-O)FeIII} core, while Mössbauer spectroscopy and d-d transitions probe the local FeIII sites. The remote methoxy substituents decrease the potential for the oxidation to FeIV by ∼100 mV, while directly bound OH- in {FeIII(OH)(μ-O)FeIII(OH)} with a short 1.91 Å FeIII-OOH bond decreases the potential by 590 mV compared to {FeIII(OAc)(μ-O)FeIII(OAc)} with a 2.01 Å FeIII-OOAc bond. Interestingly, this FeIII-OH bond is even shorter (1.87 Å) in the mixed OH-/OAc- complex but the potential is the mean value of the potentials of the OH-/OH- and OAc-/OAc- complexes, thus reflecting the electron density of the central {FeIII(μ-O)FeIII} core and not of the local FeIII-OH unit.

Discovery of a Druggable, Cryptic Pocket in SARS-CoV-2 nsp16 Using Allosteric Inhibitors

A collaborative, open-science team undertook discovery of novel small molecule inhibitors of the SARS-CoV-2 nsp16-nsp10 2'-O-methyltransferase using a high throughput screening approach with the potential to reveal new inhibition strategies. This screen yielded compound 5a, a ligand possessing an electron-deficient double bond, as an inhibitor of SARS-CoV-2 nsp16 activity. Surprisingly, X-ray crystal structures revealed that 5a covalently binds within a previously unrecognized cryptic pocket near the S-adenosylmethionine binding cleft in a manner that prevents occupation by S-adenosylmethionine. Using a multidisciplinary approach, we examined the mechanism of binding of compound 5a to the nsp16 cryptic pocket and developed 5a derivatives that inhibited nsp16 activity and murine hepatitis virus replication in rat lung epithelial cells but proved cytotoxic to cell lines canonically used to examine SARS-CoV-2 infection. Our study reveals the druggability of this newly discovered SARS-CoV-2 nsp16 cryptic pocket, provides novel tool compounds to explore the site, and suggests a new approach for discovery of nsp16 inhibition-based pan-coronavirus therapeutics through structure-guided drug design.

Comprehensive Mechanistic View of the Hydrolysis of Oxadiazole-Based Inhibitors by Histone Deacetylase 6 (HDAC6)

Histone deacetylase (HDAC) inhibitors used in the clinic typically contain a hydroxamate zinc-binding group (ZBG). However, more recent work has shown that the use of alternative ZBGs, and, in particular, the heterocyclic oxadiazoles, can confer higher isoenzyme selectivity and more favorable ADMET profiles. Herein, we report on the synthesis and biochemical, crystallographic, and computational characterization of a series of oxadiazole-based inhibitors selectively targeting the HDAC6 isoform. Surprisingly, but in line with a very recent finding reported in the literature, a crystal structure of the HDAC6/inhibitor complex revealed that hydrolysis of the oxadiazole ring transforms the parent oxadiazole into an acylhydrazide through a sequence of two hydrolytic steps. An identical cleavage pattern was also observed both in vitro using the purified HDAC6 enzyme as well as in cellular systems. By employing advanced quantum and molecular mechanics (QM/MM) and QM calculations, we elucidated the mechanistic details of the two hydrolytic steps to obtain a comprehensive mechanistic view of the double hydrolysis of the oxadiazole ring. This was achieved by fully characterizing the reaction coordinate, including identification of the structures of all intermediates and transition states, together with calculations of their respective activation (free) energies. In addition, we ruled out several (intuitively) competing pathways. The computed data (ΔG^‡ ≈ 21 kcal·mol^-1 for the rate-determining step of the overall dual hydrolysis) are in very good agreement with the experimentally determined rate constants, which a posteriori supports the proposed reaction mechanism. We also clearly (and quantitatively) explain the role of the -CF₃ or -CHF₂ substituent on the oxadiazole ring, which is a prerequisite for hydrolysis to occur. Overall, our data provide compelling evidence that the oxadiazole warheads can be efficiently transformed within the active sites of target metallohydrolases to afford reaction products possessing distinct selectivity and inhibition profiles.

A local pair-natural orbital-based complete-active space perturbation theory using orthogonal localized virtual molecular orbitals

The multireference second-order perturbation theory (CASPT2) is known to deliver a quantitative description of various complex electronic states. Despite its near-size-consistent nature, the applicability of the CASPT2 method to large, real-life systems is mostly hindered by large computational and storage costs for the two-external tensors, such as two-electron integrals, amplitudes, and residuum. To this end, Menezes and co-workers developed a reduced-scaling CASPT2 scheme by incorporating the local pair-natural orbital (PNO) representation of the many-body wave functions using non-orthonormal projected atomic orbitals (PAOs) into the CASPT theory [F. Menezes et al., J. Chem. Phys. 145, 124115 (2016)]. Alternatively, in this paper, we develop a new PNO-based CASPT2 scheme using the orthonormal localized virtual molecular orbitals (LVMOs) and assess its performance and accuracy in comparison with the conventional PAO-based counterpart. Albeit the compactness, the LVMOs were considered to perform somewhat poorly compared to PAOs in the local correlation framework because they caused enormously large orbital domains. In this work, we show that the size of LVMO domains can be rendered comparable to or even smaller than that of PAOs by the use of the differential overlap integrals for domain construction. Optimality of the MOs from the CASSCF treatment is a key to reducing the LVMO domain size for the multireference case. Due to the augmented Hessian-based localization algorithm, an additional computational cost for obtaining the LVMOs is relatively minor. We demonstrate that the LVMO-based PNO-CASPT2 method is routinely applicable to large, real-life molecules such as Menshutkin SN2 reaction in a single-walled carbon nanotube reaction field.

Properties of the Reactants and Their Interactions within and with the Enzyme Binding Cavity Determine Reaction Selectivities. The Case of Fe(II)/2-Oxoglutarate Dependent Enzymes

Fe(II)/2-oxoglutarate dependent dioxygenases (ODDs) share a double stranded beta helix (DSBH) fold and utilise a common reactive intermediate, ferryl species, to catalyse oxidative transformations of substrates. Despite the structural similarities, ODDs accept a variety of substrates and facilitate a wide range of reactions, that is hydroxylations, desaturations, (oxa)cyclisations and ring rearrangements. In this review we present and discuss the factors contributing to the observed (regio)selectivities of ODDs. They span from inherent properties of the reactants, that is, substrate molecule and iron cofactor, to the interactions between the substrate and the enzyme's binding cavity; the latter can counterbalance the effect of the former. Based on results of both experimental and computational studies dedicated to ODDs, we also line out the properties of the reactants which promote reaction outcomes other than the "default" hydroxylation. It turns out that the reaction selectivity depends on a delicate balance of interactions between the components of the investigated system.

Generation of a μ-1,2-hydroperoxo FeIIIFeIII and a μ-1,2-peroxo FeIVFeIII Complex

μ-1,2-Peroxo-diferric intermediates (P) of non-heme diiron enzymes are proposed to convert upon protonation either to high-valent active species or to activated P′ intermediates via hydroperoxo-diferric intermediates. Protonation of synthetic μ-1,2-peroxo model complexes occurred at the μ-oxo and not at the μ-1,2-peroxo bridge. Here we report a stable μ-1,2-peroxo complex {Fe^III(μ-O)(μ-1,2-O₂)Fe^III} using a dinucleating ligand and study its reactivity. The reversible oxidation and protonation of the μ-1,2-peroxo-diferric complex provide μ-1,2-peroxo Fe^IVFe^III and μ-1,2-hydroperoxo-diferric species, respectively. Neither the oxidation nor the protonation induces a strong electrophilic reactivity. Hence, the observed intramolecular C-H hydroxylation of preorganized methyl groups of the parent μ-1,2-peroxo-diferric complex should occur via conversion to a more electrophilic high-valent species. The thorough characterization of these species provides structure-spectroscopy correlations allowing insights into the formation and reactivities of hydroperoxo intermediates in diiron enzymes and their conversion to activated P′ or high-valent intermediates.

Iron coordination complexes can be used to gain insight on biologically relevant iron-oxygen compounds generated in iron metalloenzymes. Here, the authors characterise a μ-1,2-hydroperoxo Fe^IIIFe^III and a μ-1,2-peroxo Fe^IVFe^III, and study their reactivity in C-H activation.

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.

Understanding desaturation/hydroxylation activity of castor stearoyl Δ9-Desaturase through rational mutagenesis

•

Rationally designed mutations in the Δ⁹ desaturase promoted hydroxylation activity.

•

Proton and electron transfer to the active site is crucial for the Δ⁹D to desaturate

•

Detailed analysis of all enzymatic products of the Δ⁹D was carried out

•

Insight into the chemo-, and stereoselectivity of non-heme diiron enzymes was obtained

A recently proposed reaction mechanism of soluble Δ⁹ desaturase (Δ⁹D) allowed us to identify auxiliary residues His203, Asp101, Thr206 and Cys222 localized near the di-iron active site that are supposedly involved in the proton transfer (PT) to and from the active site. The PT, along with the electron transfer (ET), seems to be crucial for efficient desaturation. Thus, perturbing the major PT chains is expected to impair the native reaction and (potentially) amplify minor reaction channels, such as the substrate hydroxylation. To verify this hypothesis, we mutated the four residues mentioned above into their counterparts present in a soluble methane monooxygenase (sMMO), and determined the reaction products of mutants. We found that the mutations significantly promote residual monohydroxylation activities on stearoyl-CoA, often at the expense of native desaturation activity. The favored hydroxylation positions are C₉, followed by C₁₀ and C₁₁. Reactions with unsaturated substrate, oleoyl-CoA, yield erythro-9,10-diol, cis-9,10-epoxide and a mixture of allylic alcohols. Additionally, using 9- and 11-hydroxystearoyl-CoA, we showed that the desaturation reaction can proceed only with the hydroxyl group at position C₁₁, whereas the hydroxylation reaction is possible in both cases, i.e. with hydroxyl at position C₉ or C₁₁. Despite the fact that the overall outcome of hydroxylation is rather modest and that it is mostly the desaturation/hydroxylation ratio that is affected, our results broaden understanding of the origin of chemo- and stereoselectivity of the Δ⁹D and provide further insight into the catalytic action of the NHFe₂ enzymes.

A Desaturase-Like Enzyme Catalyzes Oxazole Formation in Pseudomonas Indolyloxazole Alkaloids

Indolyloxazole alkaloids occur in diverse micro- and macroorganisms and exhibit a wide range of pharmacological activities. Despite their ubiquitous occurrence and simple structures, the biosynthetic pathway remained unknown. Here, we used transposon mutagenesis in the labradorin producer Pseudomonas entomophila to identify a cryptic biosynthetic locus encoding an N-acyltransferase and a non-heme diiron desaturase-like enzyme. Heterologous expression in E. coli demonstrates that both enzymes are sufficient to produce indolyloxazoles. Probing their function in stable-isotope feeding experiments, we provide evidence for an unusual desaturase mechanism that generates the oxazole by decarboxylative cyclization.