Structural and functional studies reveal the molecular basis of substrate promiscuity of a glycosyltransferase originating from a major agricultural pest.

The two-spotted spider mite, Tetranychus urticae, is a major cosmopolitan pest that feeds on more than 1100 plant species. Its genome contains an unprecedentedly large number of genes involved in detoxifying and transporting xenobiotics, including 80 genes that code for UDP glycosyltransferases (UGTs). These enzymes were acquired via horizontal gene transfer from bacteria after loss in the Chelicerata lineage. UGTs are well-known for their role in phase II metabolism; however, their contribution to host adaptation and acaricide resistance in arthropods, such as T. urticae, is not yet resolved. TuUGT202A2 (Tetur22g00270) has been linked to the ability of this pest to adapt to tomato plants. Moreover, it was shown that this enzyme can glycosylate a wide range of flavonoids. To understand this relationship at the molecular level, structural, functional, and computational studies were performed. Structural studies provided specific snapshots of the enzyme in different catalytically relevant stages. The crystal structure of TuUGT202A2 in complex with UDP-glucose was obtained and site-directed mutagenesis paired with molecular dynamic simulations revealed a novel lid-like mechanism involved in the binding of the activated sugar donor. Two additional TuUGT202A2 crystal complexes, UDP-(S)-naringenin and UDP-naringin, demonstrated that this enzyme has a highly plastic and open-ended acceptor-binding site. Overall, this work reveals the molecular basis of substrate promiscuity of TuUGT202A2 and provides novel insights into the structural mechanism of UGTs catalysis.

Bacteria at Work - Experimental and Theoretical Studies Reveal the Catalytic Mechanism of Ectoine Synthase.

Ectoine synthase (EctC) catalyses the ultimate step of ectoine biosynthesis, a kosmotropic compound produced as compatible solute by many bacteria and some archaea or eukaryotes. EctC is an Fe2+-dependent homodimeric cytoplasmic protein. Using Mössbauer spectroscopy, molecular dynamics simulations and QM/MM calculations, we determined the most likely coordination number and geometry of the Fe2+ ion and proposed a mechanism of the EctC-catalysed reaction. Most notably, we show that apart from the three amino acids binding to the iron ion (Glu57, Tyr84 and His92), one water molecule and one hydroxide ion are required as additional ligands for the reaction to occur. They fill the first coordination sphere of the Fe2+-cofactor and act as critical proton donors and acceptors during the cyclization reaction.

Reaction Mechanism Catalyzed by the Dissimilatory Sulfite Reductase - The Role of the Siroheme-[4FeS4] Cofactor.

The present work is another part of our investigation on the pathway of dissimilatory sulfate reduction and covers a theoretical study on the reaction catalyzed by dissimilatory sulfite reductase (dSIR). dSIR is the terminal enzyme involved in this metabolic pathway, which uses the siroheme-[4Fe4S] cofactor for six-electron reduction of sulfite to sulfide. In this study we use a large cluster model containing siroheme-[4Fe4S] cofactor and protein residues involved in the direct interactions with the substrate, to get insight into the most feasible reaction mechanism and to understand the role of each considered active site component. In combination with earlier studies reported in the literature, our results lead to several interesting insights. One of the most important conclusions is that the reaction mechanism consists of three steps of two-electron reduction of sulfur and the probable role of the siroheme-[4Fe4S] cofactor is to ensure the delivery of packages of two electrons to the reactant.

Development of fluorophore labeled or biotinylated anticancer small molecule NSC243928.

Small molecule NSC243928 binds with LY6K, a potential target for the treatment of triple-negative breast cancer, and induces cancer cell death with an unclear mechanism. We have developed chemical tools to identify the molecular mechanisms of NSC243928-LY6K interaction. Herein, we report on the development and synthesis of biotinylated and fluorophore-tethered derivatives of NSC243928 guided by docking studies and molecular dynamics. Surface plasmon resonance assay indicates that these derivatives retained a direct binding with LY6K protein. Confocal analysis revealed that nitrobenzoxadiazole (NBD) fluorophore tagged NSC243928 is retained in LY6K expressing cancer cells. These novel modified compounds will be employed in future in vitro and in vivo studies to understand the molecular mechanisms of NSC243928 mediated cancer cell death. These studies will pave the path for developing novel targeted therapeutics and understanding any potential side-effects of these treatments for hard-to-treat cancers such as triple-negative breast cancer or other cancers with high expression of LY6K.

Experimental and Computational Studies Reveal Novel Interaction of Lymphocytes Antigen 6K to TGF-ß Receptor Complex.

TGF-ß signaling promotes migration, invasion, and distant colonization of cancer cells in advanced metastatic cancers. TGF-ß signaling suppresses the anti-tumor immune response in a tumor microenvironment, allowing sustained tumor growth. TGF-ß plays an important role in normal physiology; thus it is no surprise that the clinical development of effective and safe TGF-ß inhibitors has been hampered due to their high toxicity. We discovered that increased expression of LY6K in cancer cells led to increased TGF-ß signaling and that inhibition of LY6K could lead to reduced TGF-ß signaling and reduced in vivo tumor growth. LY6K is a highly cancer-specific protein, and it is not expressed in normal organs except in the testes. Thus, LY6K is a valid target for developing therapeutic strategies to inhibit TGF-ß signaling in cancer cells. We employed in vitro pull-down assays and molecular dynamics simulations to understand the structural determinants of the TGF-ß receptor complex with LY6K. This combined approach allowed us to identify the critical residues and dynamics of the LY6K interaction with the TGF-ß receptor complex. These data are critical in designing novel drugs for the inhibition of TGF-ß in LY6K expressing cancer, induction of anti-tumor immune response, and inhibition of tumor growth and metastatic spread.

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.

Enzyme Multifunctionality by Control of Substrate Positioning Within the Catalytic Cycle-A QM/MM Study of Clavaminic Acid Synthase.

Clavaminic acid synthase from Streptomyces clavuligerus is an FeII /2-oxoglutarate-dependent dioxygenase, crucial for the biosynthesis of the ß-lactamase inhibitor clavulanic acid. It catalyses three consecutive oxidative reactions, that is, hydroxylation, cyclisation and desaturation, in a single binding cavity. As follows from the results of this QM/MM study, CAS versatility and selectivity depends on the binding cavity, which interplays differently with the substrate for each reaction. The enzyme-substrate interactions affect the substrate's ability to re-position during the reaction, either constraining it in its primary position, which impedes processes other than oxygen rebound, or allowing change, which facilitates desaturation. This differential effect originates from two aspartate residues, which strongly interact with the guanidine group of the hydroxylation substrate and stabilise the orientation of the molecule. These residues interact less effectively with the smaller amine group of the desaturation substrate(s), aiding their re-positioning and the subsequent formation of a double bond.

Substitution of histidine 30 by asparagine in manganese superoxide dismutase alters biophysical properties and supports proliferation in a K562 leukemia cell line.

We have generated a mutant of C. elegans manganese superoxide dismutase at histidine 30 by site-directed mutagenesis. The structure was solved at a resolution of 1.52 Å by X-ray crystallography (pdb: 6S0D). His30 was targeted, as it forms as a gateway residue at the top of the solvent access funnel to the active site, together with Tyr34. In the wild-type protein, these gateway residues are involved in the hydrogen-bonding network providing the protons necessary for the catalytic reaction at the metal center. However, biophysical characterization and cell viability experiments reveal that a mutation from histidine to asparagine in the H30N mutant modifies metal selectivity in the protein, favoring the uptake of iron over manganese in minimal media conditions, alters active-site coordination from the characteristic trigonal bipyramidal to octahedral geometry, and encourages cellular proliferation in K562 cells, when added exogenously to the cells.

Crystal structure of thebaine 6-O-demethylase from the morphine biosynthesis pathway.

Thebaine 6-O-demethylase (T6ODM) from Papaver somniferum (opium poppy), which belongs to the non-heme 2-oxoglutarate/Fe(II)-dependent dioxygenases (ODD) family, is a key enzyme in the morphine biosynthesis pathway. Initially, T6ODM was characterized as an enzyme catalyzing O-demethylation of thebaine to neopinone and oripavine to morphinone. However, the substrate range of T6ODM was recently expanded to a number of various benzylisoquinoline alkaloids. Here, we present crystal structures of T6ODM in complexes with 2-oxoglutarate (T6ODM:2OG, PDB: 5O9W) and succinate (T6ODM:SIN, PDB: 5O7Y). Both metal and 2OG binding sites display similarity to other proteins from the ODD family, but T6ODM is characterized by an exceptionally large substrate binding cavity, whose volume can partially explain the promiscuity of this enzyme. Moreover, the size of the cavity allows for binding of multiple molecules at once, posing a question about the substrate-driven specificity of the enzyme.

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

The Fe(II)/2-oxoglutarate-dependent dioxygenase AsqJ from Aspergillus nidulans catalyses two pivotal steps in the synthesis of quinolone antibiotic 4'-methoxyviridicatin, i.e., desaturation and epoxidation of a benzodiazepinedione. The previous experimental results signal that, during the desaturation reaction, hydrogen atom transfer (HAT) from the benzylic carbon atom (C10) is a more likely step to initiate the reaction than the alternative HAT from the ring moiety (C3 atom). To unravel the origins of this regioselectivity and to explain why the observed reaction is desaturation and not the "default" hydroxylation, we performed a QM/MM study on the reaction catalysed by AsqJ. Herein, we report results that complement the experimental findings and suggest that HAT at the C10 position is the preferred reaction due to favourable interactions between the substrate and the binding cavity that compensate for the relatively high intrinsic barrier associated with the process. For the resultant radical intermediate, the desaturation/hydroxylation selectivity is governed by electronic properties of the reactants, i.e., the energy gap between the orbital that hosts the unpaired electron and the sigma orbital for the C-H bond as well as the gap between the orbitals mixing in transition state structures for each elementary step. Regiospecificity of the AsqJ dehydrogenation reaction is dictated by substrate-protein interactions. 82 × 44 mm (300 × 300 dpi).

On the Structure and Reaction Mechanism of Human Acireductone Dioxygenase.

Acireductone dioxygenase (ARD) is an intriguing enzyme from the methionine salvage pathway that is capable of catalysing two different oxidation reactions with the same substrate depending on the type of the metal ion in the active site. To date, the structural information regarding the ARD-acireductone complex is limited and possible reaction mechanisms are still under debate. The results of joint experimental and computational studies undertaken to advance knowledge about ARD are reported. The crystal structure of an ARD from Homo sapiens was determined with selenomethionine. EPR spectroscopy suggested that binding acireductone triggers one protein residue to dissociate from Fe2+ , which allows NO (and presumably O2 ) to bind directly to the metal. Mössbauer spectroscopic data (interpreted with the aid of DFT calculations) was consistent with bidentate binding of acireductone to Fe2+ and concomitant dissociation of His88 from the metal. Major features of Fe vibrational spectra obtained for the native enzyme and upon addition of acireductone were reproduced by QM/MM calculations for the proposed models. A computational (QM/MM) study of the reaction mechanisms suggests that Fe2+ promotes O-O bond homolysis, which elicits cleavage of the C1-C2 bond of the substrate. Higher M3+ /M2+ redox potentials of other divalent metals do not support this pathway, and instead the reaction proceeds similarly to the key reaction step in the quercetin 2,3-dioxygenase mechanism.

A Single Mutation is Sufficient to Modify the Metal Selectivity and Specificity of a Eukaryotic Manganese Superoxide Dismutase to Encompass Iron.

We have generated a site-directed mutant of the manganese superoxide dismutase SOD-3 of C.elegans (MnSOD-3) which modifies the metal specificity of the enzyme. While wild-type MnSOD-3 functions with manganese in the active site (3600 U mg-1 of protein) it has little or no activity when iron is incorporated. However, when histidine replaces glutamine 142 in the active site, the enzyme retains 50 % of its activity and becomes cambialistic for its metal cofactor exhibiting very similar specific activity with either manganese or iron.

Bacterial steroid hydroxylases: enzyme classes, their functions and comparison of their catalytic mechanisms.

The steroid superfamily includes a wide range of compounds that are essential for living organisms of the animal and plant kingdoms. Structural modifications of steroids highly affect their biological activity. In this review, we focus on hydroxylation of steroids by bacterial hydroxylases, which take part in steroid catabolic pathways and play an important role in steroid degradation. We compare three distinct classes of metalloenzymes responsible for aerobic or anaerobic hydroxylation of steroids, namely: cytochrome P450, Rieske-type monooxygenase 3-ketosteroid 9α-hydroxylase, and molybdenum-containing steroid C25 dehydrogenases. We analyze the available literature data on reactivity, regioselectivity, and potential application of these enzymes in organic synthesis of hydroxysteroids. Moreover, we describe mechanistic hypotheses proposed for all three classes of enzymes along with experimental and theoretical evidences, which have provided grounds for their formulation. In case of the 3-ketosteroid 9α-hydroxylase, such a mechanistic hypothesis is formulated for the first time in the literature based on studies conducted for other Rieske monooxygenases. Finally, we provide comparative analysis of similarities and differences in the reaction mechanisms utilized by bacterial steroid hydroxylases.

Correction to: Bacterial steroid hydroxylases: enzyme classes, their functions and comparison of their catalytic mechanisms.

The published online version contains mistake in the author list. The correct presentation should have been "Rita Bernhardt" instead of "Rita Bernhard". There was a missing "t" on the original publication.

DFT study of the mechanism of manganese quercetin 2,3-dioxygenase: quest for origins of enzyme unique nitroxygenase activity and regioselectivity.

Quercetin 2,3-dioxygenase (QDO) is an enzyme which accepts various transition metal ions as cofactors, and cleaves the heterocyclic ring of quercetin with consumption of dioxygen and release of carbon monoxide. QDO from B. subtilis that binds Mn(II) displays an unprecedented nitroxygenase activity, whereby nitroxyl (HNO) is incorporated into quercetin cleavage products instead of dioxygen. Interestingly, the reaction proceeds with high regiospecificity, i.e., nitrogen and oxygen atoms of HNO are incorporated into specific fragments of the cleavage product. A nonenzymatic base-catalyzed reaction, which occurs in pH above 7.5, yields the same reaction products. Herein, we report results of quantum chemical studies on the mechanisms of the nitroxygenase reaction of Mn-QDO. Density functional method with dispersion correction (B3LYP-D3) was applied to the Mn-QDO active site model and the reactants of the nonenzymatic reaction. Co(II)- and Fe(II)-variants of the active site were also considered. Analysis of reaction energy profiles suggests that the regiospecificity of the reaction is an inherent property of the reactants, whereas the unique reactivity of Mn-QDO, as opposed to Co- or Fe-QDO that do not catalyze nitroxygenation, stems from weak HNO binding and lack of strong preference for coordination of HNO through the nitrogen atom. Moreover, the enzyme activates quercetin through deprotonation and the proton acceptor-Glu69 needs to reorient for the reaction to proceed.

Quantum Mechanics/Molecular Mechanics Modeling of Enzymatic Processes: Caveats and Breakthroughs.

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.

Anion Effects in Oxidative Aliphatic Carbon-Carbon Bond Cleavage Reactions of Cu(II) Chlorodiketonate Complexes.

Aliphatic oxidative carbon-carbon bond cleavage reactions involving Cu(II) catalysts and O2 as the terminal oxidant are of significant current interest. However, little is currently known regarding how the nature of the Cu(II) catalyst, including the anions present, influence the reaction with O2. In previous work, we found that exposure of the Cu(II) chlorodiketonate complex [(6-Ph2TPA)Cu(PhC(O)CClC(O)Ph)]ClO4 (1) to O2 results in oxidative aliphatic carbon-carbon bond cleavage within the diketonate unit, leading to the formation of benzoic acid, benzoic anhydride, benzil, and 1,3-diphenylpropanedione as organic products. Kinetic studies of this reaction revealed a slow induction phase followed by a rapid decay of the absorption features of 1. Notably, the induction phase is not present when the reaction is performed in the presence of a catalytic amount of chloride anion. In the studies presented herein, a combination of spectroscopic (UV-vis, EPR) and density functional theory (DFT) methods have been used to examine the chloride and benzoate ion binding properties of 1 under anaerobic conditions. These studies provide evidence that each anion coordinates in an axial position of the Cu(II) center. DFT studies reveal that the presence of the anion in the Cu(II) coordination sphere decreases the barrier for O2 activation and the formation of a Cu(II)-peroxo species. Notably, the chloride anion more effectively lowers the barrier associated with O-O bond cleavage. Thus, the nature of the anion plays an important role in determining the rate of reaction of the diketonate complex with O2. The same type of anion effects were observed in the O2 reactivity of the simple Cu(II)-bipyridine complex [(bpy)Cu(PhC(O)C(Cl)C(O)Ph)ClO4] (3).

Mechanism of O2 Activation by α-Ketoglutarate Dependent Oxygenases Revisited. A Quantum Chemical Study.

Four mechanisms previously proposed for dioxygen activation catalyzed by α-keto acid dependent oxygenases (α-KAO) were studied with dispersion-corrected DFT methods employing B3LYP and TPSSh functionals in combination with triple-ζ basis set (cc-pVTZ). The aim of this study was to revisit mechanisms suggested in the past decade and resolve remaining issues related to dioxygen activation. Mechanism A, which runs on the quintet potential energy surface (PES) and includes formation of an Fe(III)-superoxide radical anion complex, subsequent oxidative decarboxylation, and O-O bond cleavage, was found to be most likely. However, mechanism B taking place on the septet PES involves a rate limiting barrier comparable to the one found for mechanism A, and thus it cannot be excluded, though two other mechanisms (C and D) were ruled out. Mechanism C is a minor variation of mechanism A, whereas mechanism D proceeds through formation of a triplet Fe(IV)-alkyl peroxo bridged intermediate. The study covered also full optimization of relevant minimum energy crossing points (MECPs). The relative energy of critical intermediates was also studied with the CCSD(T) method in order to benchmark TPSSh and B3LYP functionals with respect to their credibility in predicting relative energies of septet and triplet spin states of the ternary enzyme-Fe-α-keto glutarate (α-KG)-O2 complex.

Role of substrate positioning in the catalytic reaction of 4-hydroxyphenylpyruvate dioxygenase-A QM/MM Study.

Ring hydroxylation and coupled rearrangement reactions catalyzed by 4-hydroxyphenylpyruvate dioxygenase were studied with the QM/MM method ONIOM(B3LYP:AMBER). For electrophilic attack of the ferryl species on the aromatic ring, five channels were considered: attacks on the three ring atoms closest to the oxo ligand (C1, C2, C6) and insertion of oxygen across two bonds formed by them (C1-C2, C1-C6). For the subsequent migration of the carboxymethyl substituent, two possible directions were tested (C1→C2, C1→C6), and two different mechanisms were sought (stepwise radical, single-step heterolytic). In addition, formation of an epoxide (side)product and benzylic hydroxylation, as catalyzed by the closely related hydroxymandelate synthase, were investigated. From the computed reaction free energy profiles it follows that the most likely mechanism of 4-hydroxyphenylpyruvate dioxygenase involves electrophilic attack on the C1 carbon of the ring and subsequent single-step heterolytic migration of the substituent. Computed values of the kinetic isotope effect for this step are inverse, consistent with available experimental data. Electronic structure arguments for the preferred mechanism of attack on the ring are also presented.

Halide-promoted dioxygenolysis of a carbon-carbon bond by a copper(II) diketonate complex.

A mononuclear Cu(II) chlorodiketonate complex was prepared, characterized, and found to undergo oxidative aliphatic carbon-carbon bond cleavage within the diketonate unit upon exposure to O2 at ambient temperature. Mechanistic studies provide evidence for a dioxygenase-type C-C bond cleavage reaction pathway involving trione and hypochlorite intermediates. Significantly, the presence of a catalytic amount of chloride ion accelerates the oxygen activation step via the formation of a Cu-Cl species, which facilitates monodentate diketonate formation and lowers the barrier for O2 activation. The observed reactivity and chloride catalysis is relevant to Cu(II) halide-catalyzed reactions in which diketonates are oxidatively cleaved using O2 as the terminal oxidant. The results of this study suggest that anion coordination can play a significant role in influencing copper-mediated oxygen activation in such systems.