On how the binding cavity of AsqJ dioxygenase controls the desaturation reaction regioselectivity: a QM/MM study

The Unique Role of the Second Coordination Sphere to Unlock and Control Catalysis in Nonheme Fe(II)/2-Oxoglutarate Histone Demethylase KDM2A

Nonheme Fe(II) and 2-oxoglutarate (2OG)-dependent histone lysine demethylases 2A (KDM2A) catalyze the demethylation of the mono- or dimethylated lysine 36 residue in the histone H3 peptide (H3K36me1/me2), which plays a crucial role in epigenetic regulation and can be involved in many cancers. Although the overall catalytic mechanism of KDMs has been studied, how KDM2 catalysis takes place in contrast to other KDMs remains unknown. Understanding such differences is vital for enzyme redesign and can help in enzyme-selective drug design. Herein, we employed molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) to explore the complete catalytic mechanism of KDM2A, including dioxygen diffusion and binding, dioxygen activation, and substrate oxidation. Our study demonstrates that the catalysis of KDM2A is controlled by the conformational change of the second coordination sphere (SCS), specifically by a change in the orientation of Y222, which unlocks the 2OG rearrangement from off-line to in-line mode. The study demonstrates that the variant Y222A makes the 2OG rearrangement more favorable. Furthermore, the study reveals that it is the size of H3K36me3 that prevents the 2OG rearrangement, thus rendering the enzyme inactivity with trimethylated lysine. Calculations show that the SCS and long-range interacting residues that stabilize the HAT transition state in KDM2A differ from those in KDM4A, KDM7B, and KDM6A, thus providing the basics for the enzyme-selective redesign and modulation of KDM2A without influencing other KDMs.

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

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

Elucidation of divergent desaturation pathways in the formation of vinyl isonitrile and isocyanoacrylate

Two different types of desaturations are employed by iron- and 2-oxoglutarate-dependent (Fe/2OG) enzymes to construct vinyl isonitrile and isocyanoacrylate moieties found in isonitrile-containing natural products. A substrate-bound protein structure reveals a plausible strategy to affect desaturation and hints at substrate promiscuity of these enzymes. Analogs are synthesized and used as mechanistic probes to validate structural observations. Instead of proceeding through hydroxylated intermediate as previously proposed, a plausible carbocation species is utilized to trigger C=C bond installation. These Fe/2OG enzymes can also accommodate analogs with opposite chirality and different functional groups including isonitrile-(D)-tyrosine, N-formyl tyrosine, and phloretic acid, while maintaining the reaction selectivity.

While canonical oxygen-rebound pathway has been utilized by several Fe/2OG enzymes, insights into the desaturation process remain to be elucidated. Here, the authors suggest that a carbocation might be deployed to enable chemically divergent desaturations through decarboxylation and deprotonation.

Peroxy Intermediate Drives Carbon Bond Activation in the Dioxygenase AsqJ

Dioxygenases catalyze stereoselective oxygen atom transfer in metabolic pathways of biological, industrial, and pharmaceutical importance, but their precise chemical principles remain controversial. The α-ketoglutarate (αKG)-dependent dioxygenase AsqJ synthesizes biomedically active quinolone alkaloids via desaturation and subsequent epoxidation of a carbon-carbon bond in the cyclopeptin substrate. Here, we combine high-resolution X-ray crystallography with enzyme engineering, quantum-classical (QM/MM) simulations, and biochemical assays to describe a peroxidic intermediate that bridges the substrate and active site metal ion in AsqJ. Homolytic cleavage of this moiety during substrate epoxidation generates an activated high-valent ferryl (Fe^IV = O) species that mediates the next catalytic cycle, possibly without the consumption of the metabolically valuable αKG cosubstrate. Our combined findings provide an important understanding of chemical bond activation principles in complex enzymatic reaction networks and molecular mechanisms of dioxygenases.

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.

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

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

Electrostatic Perturbations from the Protein Affect C-H Bond Strengths of the Substrate and Enable Negative Catalysis in the TmpA Biosynthesis Enzyme

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C¹ -position of the substrate as part of the biosynthesis of glycine betaine in bacteria and carnitine in humans. To understand how the enzyme avoids breaking the weak C² -H bond in favor of C¹ -hydroxylation, we set up a cluster model of 242 atoms representing the first and second coordination sphere of the metal center and substrate binding pocket, and investigated possible reaction mechanisms of substrate activation by an iron(IV)-oxo species by density functional theory methods. In agreement with experimental product distributions, the calculations predict a favorable C¹ -hydroxylation pathway. The calculations show that the selectivity is guided through electrostatic perturbations inside the protein from charged residues, external electric fields and electric dipole moments. In particular, charged residues influence and perturb the homolytic bond strength of the C¹ -H and C² -H bonds of the substrate, and strongly strengthens the C² -H bond in the substrate-bound orientation.

Fungal Dioxygenase AsqJ Is Promiscuous and Bimodal: Substrate-Directed Formation of Quinolones versus Quinazolinones

Previous studies showed that the Fe^II/α‐ketoglutarate dependent dioxygenase AsqJ induces a skeletal rearrangement in viridicatin biosynthesis in Aspergillus nidulans, generating a quinolone scaffold from benzo[1,4]diazepine‐2,5‐dione substrates. We report that AsqJ catalyzes an additional, entirely different reaction, simply by a change in substituent in the benzodiazepinedione substrate. This new mechanism is established by substrate screening, application of functional probes, and computational analysis. AsqJ excises H₂CO from the heterocyclic ring structure of suitable benzo[1,4]diazepine‐2,5‐dione substrates to generate quinazolinones. This novel AsqJ catalysis pathway is governed by a single substituent within the complex substrate. This unique substrate‐directed reactivity of AsqJ enables the targeted biocatalytic generation of either quinolones or quinazolinones, two alkaloid frameworks of exceptional biomedical relevance.

Regioselectivity of hyoscyamine 6β-hydroxylase-catalysed hydroxylation as revealed by high-resolution structural information and QM/MM calculations

Hyoscyamine 6β-hydroxylase (H6H) is a bifunctional non-heme 2-oxoglutarate/Fe2+-dependent dioxygenase that catalyzes the two final steps in the biosynthesis of scopolamine. Based on high resolution crystal structures of H6H from Datura metel, detailed information on substrate binding was obtained that provided insights into the onset of the enzymatic process. In particular, the role of two prominent residues was revealed - Glu-116 that interacts with the tertiary amine located on the hyoscyamine tropane moiety and Tyr-326 that forms CH-π hydrogen bonds with the hyoscyamine phenyl ring. The structures were used as the basis for QM/MM calculations that provided an explanation for the regioselectivity of the hydroxylation reaction on the hyoscyamine tropane moiety (C6 vs. C7) and quantified contributions of active site residues to respective barrier heights.

Selective Hydrogen Atom Abstraction from Dihydroflavonol by a Nonheme Iron Center Is the Key Step in the Enzymatic Flavonol Synthesis and Avoids Byproducts

The plant non-heme iron dioxygenase flavonol synthase performs a regioselective desaturation reaction as part of the biosynthesis of the signaling molecule flavonol that triggers the growing of leaves and flowers. These compounds also have health benefits for humans. Desaturation of aliphatic compounds generally proceeds through two consecutive hydrogen atom abstraction steps from two adjacent carbon atoms and in nature often is performed by a high-valent iron(IV)-oxo species. We show that the order of the hydrogen atom abstraction steps, however, is opposite of those expected from the C-H bond strengths in the substrate and determines the product distributions. As such, flavonol synthase follows a negative catalysis mechanism. Using density functional theory methods on large active-site model complexes, we investigated pathways for desaturation and hydroxylation by an iron(IV)-oxo active-site model. Contrary to thermochemical predictions, we find that the oxidant abstracts the hydrogen atom from the strong C²-H bond rather than the weaker C³-H bond of the substrate first. We analyze the origin of this unexpected selective hydrogen atom abstraction pathway and find that the alternative C³-H hydrogen atom abstraction would be followed by a low-energy and competitive substrate hydroxylation mechanism hence, should give considerable amount of byproducts. Our computational modeling studies show that substrate positioning in flavonol synthase is essential, as it guides the reactivity to a chemo- and regioselective substrate desaturation from the C²-H group, leading to desaturation products efficiently.

On the reaction mechanism of an endoperoxide ring formation by fumitremorgin B endoperoxidase. The right arrangement makes a difference

Fumitremorgin B endoperoxidase (FtmF) belongs to 2-oxoglutarate dependent dioxygenases and catalyzes an unusual oxidative reaction of endoperoxide formation at the final stage of biosynthesis of verruculogen - a mycotoxin produced by Aspergillus and Penicillinum strains. The published crystal structure of FtmF (PDB: ), which is of overall good quality, contains a model of the substrate bound in the active site, which, however, has very low occupancy and its conformation does not comply with the small molecule crystal structure. Moreover, a previous computational study that employed a model based on this crystal structure revealed a substantial reaction barrier, which might indicate that the model of FtmF/substrate complex can have serious errors. The purpose of this work was to model with computational methods the structure of the enzyme-substrate complex and to investigate the mechanisms of the enzymatic reaction. Docking, molecular dynamics simulation and DFT results, all indicate the substrate most likely binds in the active site in a configuration very different from that originally suggested. Moreover, for this newly proposed structure of the enzyme-substrate complex, the reaction energy profile is characterised exclusively by low barriers and it successfully explains the observed regiospecificity of the enzymatic process. Finally, a plausible binding site for ascorbate was found and it is suggested that ascorbate is involved in the final step of the FtmF reaction.

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.

Mechanistic Insights into a Stibene Cleavage Oxygenase NOV1 from Quantum Mechanical/Molecular Mechanical Calculations

NOV1, a stilbene cleavage oxygenase, catalyzes the cleavage of the central double bond of stilbenes to two phenolic aldehydes, using a 4-His Fe(II) center and dioxygen. Herein, we use in-protein quantum mechanical/molecular mechanical (QM/MM) calculations to elucidate the reaction mechanism of the central double bond cleavage of phytoalexin resveratrol by NOV1. Our results showed that the oxygen molecule prefers to bind to the iron center in a side-on fashion, as suggested from the experiment. The quintet Fe-O2 complex with the side-on superoxo antiferromagnetic coupled to the resveratrol radical is identified as the reactive oxygen species. The QM/MM results support the dioxygenase mechanism involving a dioxetane intermediate with a rate-limiting barrier of 10.0 kcal mol-1. The alternative pathway through an epoxide intermediate is ruled out due to a larger rate-limiting barrier (26.8 kcal mol-1). These findings provide important insight into the catalytic mechanism of carotenoid cleavage oxygenases and also the dioxygen activation of non-heme enzymes.

Mechanistic insights into ring cleavage of hydroquinone by PnpCD from quantum mechanical/molecular mechanical calculations

QM/MM calculations for ring cleavage of hydroquinone by PnpCD show that Asn258 loses coordination to the iron when the reaction begins. The first-sphere Glu262 can act as an acid–base catalyst to lower the rate-limiting barrier.

Computational Understanding of the Selectivities in Metalloenzymes

Metalloenzymes catalyze many different types of biological reactions with high efficiency and remarkable selectivity. The quantum chemical cluster approach and the combined quantum mechanics/molecular mechanics methods have proven very successful in the elucidation of the reaction mechanism and rationalization of selectivities in enzymes. In this review, recent progress in the computational understanding of various selectivities including chemoselectivity, regioselectivity, and stereoselectivity, in metalloenzymes, is discussed.

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.