PDF) Tributions Package from CPAN Was Used for Obtaining P-Values

Total Page:16

File Type:pdf, Size:1020Kb

PDF) Tributions Package from CPAN Was Used for Obtaining P-Values A Measure of the Promiscuity of Proteins and Characteristics of Residues in the Vicinity of the Catalytic Site That Regulate Promiscuity Sandeep Chakraborty, Basuthkar J. Rao* Department of Biological Sciences, Tata Institute of Fundamental Research, Mumbai, India, Abstract Promiscuity, the basis for the evolution of new functions through ‘tinkering’ of residues in the vicinity of the catalytic site, is yet to be quantitatively defined. We present a computational method Promiscuity IndicesEstimator (PROMISE) - based on signatures derived from the spatial and electrostatic properties of the catalytic residues, to estimate the promiscuity (PromIndex) of proteins with known active site residues and 3D structure. PromIndex reflects the number of different active site signatures that have congruent matches in close proximity of its native catalytic site, the quality of the matches and difference in the enzymatic activity. Promiscuity in proteins is observed to follow a lognormal distribution (m = 0.28, s = 1.1 reduced chi-square = 3.0E-5). The PROMISE predicted promiscuous functions in any protein can serve as the starting point for directed evolution experiments. PROMISE ranks carboxypeptidase A and ribonuclease A amongst the more promiscuous proteins. We have also investigated the properties of the residues in the vicinity of the catalytic site that regulates its promiscuity. Linear regression establishes a weak correlation (R2,0.1) between certain properties of the residues (charge, polar, etc) in the neighborhood of the catalytic residues and PromIndex. A stronger relationship states that most proteins with high promiscuity have high percentages of charged and polar residues within a radius of 3 A˚ of the catalytic site, which is validated using one-tailed hypothesis tests (P-values,0.05). Since it is known that these characteristics are key factors in catalysis, their relationship with the promiscuity index cross validates the methodology of PROMISE. Citation: Chakraborty S, Rao BJ (2012) A Measure of the Promiscuity of Proteins and Characteristics of Residues in the Vicinity of the Catalytic Site That Regulate Promiscuity. PLoS ONE 7(2): e32011. doi:10.1371/journal.pone.0032011 Editor: Bin Xue, University of South Florida, United States of America Received October 19, 2011; Accepted January 19, 2012; Published February 16, 2012 Copyright: ß 2012 Chakraborty, Rao. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This work was funded by the Tata Institute of Fundamental Research (Department of Atomic Energy). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * E-mail: [email protected] Introduction analysis of multispecific enzymes. ‘Rigorous and quantitative measures of promiscuity’ which will measure the ‘magnitude and Jensen first proposed that promiscuity shaped the evolution of degree of promiscuity in a wide range of proteins’ is thus required primitive cells which presumably had minimal gene content [1–3]. [8]. Biochemical ‘leakiness’/‘messiness’ supplemented by gene dupli- We present a computational method for assigning a relative cation ensured increased gene content and specialization [4,5]. promiscuous index (PromIndex) to proteins with known active site The catalysis of reactions distinct from the one the protein has and 3D structure - Promiscuity IndicesEstimator (PROMISE). evolved to perform, but using the same domain, is a definition of PROMISE is based on the previously described method for active promiscuity adopted by many researchers including the current site detection which relies on the spatial and electrostatic work [6–8]. Promiscuity is distinguished from moonlighting properties of the catalytic residues (CLASP) [19]. Modular functions which are typically catalyzed using a domain of the approaches similar to CLASP using only spatial congruence have protein different from the active site scaffold [9,10]. been previously adopted [20–24]. Pruning based on electrostatic Protein engineers tailor innovative proteins by honing an existing moonlighting or promiscuous activity [11–14], often properties done by CLASP reduces considerable false positives without compromising its native function [15]. Inspite of the compared to purely 3D matching method, as potential difference mushrooming number of promiscuous proteins being discovered, congruence implies an appropriate milieu in the catalytic site [19]. till date there is no quantitative measure of promiscuity of proteins. Another innovative method [25] tries to dock high-energy A recent study traces the evolution of domain promiscuity along intermediates of various metabolites listed in a database [26] various evolutionary pathways [16], while another method and successfully predicts the function of an unknown protein. proposed a numerical analysis of promiscuous domains from CLASP was applied on alkaline phosphatases (APs), one of the genomic sequences [17]. An attempt at quantifying promiscuity widely studied promiscuous enzymes [27], which led to the provides a measure of the catalytic efficiencies of an enzyme discovery of a promiscuous protease activity in shrimp AP [19], toward a pre-defined set of substrates, but is limited in its scope and a promiscuous metallo-beta-lactamase activity in E. Coli AP and scalability [18]. Since this method assumes a uniform (S. Chakraborty, R. Minda, L. Salaye, J.M. Frere, Basuthkar J. chemical transformation on all substrates, it is more apt for the Rao, Manuscript in preparation). PLoS ONE | www.plosone.org 1 February 2012 | Volume 7 | Issue 2 | e32011 Quantitative Measure of Promiscuity of Proteins A set of non-homologous proteins with known active sites and the PROMISE computed promiscuity. The PromIndex of each structures from the Catalytic Atlas Site (CSA) [28], ‘‘an excellent protein has been calculated for at least 30 values, i.e. each protein starting point for characterizing mechanistically diverse superfam- figures at least 30 times in the randomly selected sets used for ilies’’ [29], was supplemented to include other promiscuous calculating PromIndex. Figure 1B shows the PromIndex comput- proteins described in a recent review [8] (Sproteins: |Spro- ed on the full population (which is skewed with respect to the teins| = 305).. PromIndex is now computed for each protein in the representation of each EC number) as compared with the mean set reflecting the number of matching active site signatures from and median obtained from balanced sets used in Figure 1A. other proteins in the vicinity of the native catalytic site, the quality Figure 1A and Figure 1B have the same ordering of proteins - but of the matches and how much the EC numbers differ. Since the set Figure 1A is based on (many) balanced sets and Figure 1B is based of proteins is unbalanced with respect to the EC numbers, random on the full set. As can be seen, the indices follow the same trend - sample sets such that all Enzyme Commission (EC) numbers are and is not randomly scattered. At first glance it appears that this equally represented are created from Sproteins. Computation of imbalance biases the promiscuity - since there are more proteins PromIndex was repeated multiple times by such random sampling from EC 3 (hydrolases), we expect hydrolases to be more to obtain statistics like mean, median and standard deviation, and promiscuous in this computation. However, this bias can be PromIndex was fitted for a lognormal distribution. PROMISE partially negated by assigning more weightage to matches across ranks a carboxypeptidase A, a catechol 2,3-dioxygenase, a ECs. Thus when we consider a protein from EC 6, even a single phosphoenolpyruvate carboxykinase, a kynureninase and a extra match with a protein from EC 3 (which are more ribonuclease-A as the most promiscuous, and a ribonucleotide predominant in this set), will add more to the promiscuity index reductase, a lipase and a haloalkane dehalogenase as the least neutralizing the effect of the higher number of EC 3 proteins. It is promiscuous proteins in Sproteins. We also establish a weak linear seen that the value of the weights chosen has virtually no effect on relationship (R2,0.1) between PromIndex and certain character- PromIndex. This is expected since the PromIndex is normalized. istics (% of charged and polar residues) in the immediate Figure 1C shows the variation in PromIndex when the radius neighborhood of the catalytic residues. Since these properties of around the catalytic residues which is to be considered as part of residues are known to be intrinsic to catalysis, this relationship the active site is varied from 1–6 A˚ . Note that the proteins are cross validates the methodology to quantify promiscuity. While the sorted according to the PromIndex as computed for a radius of forward relationship is weak (i.e. high percentage of polar residues 5A˚ . Figure 1D shows the distribution of PromIndex for the whole does not necessarily imply high promiscuity), the reverse population. This was fitted for a lognormal distribution relationship is quite strong (i.e. if the percentage of polar residues . is low, we can say with high confidence that the protein is not very ÀÁ ~ -y ~ 2 2 promiscuous). We establish low P-values (,0.05) based on fxðÞx; m,s e xsHðÞ2p ,y lnðÞ x -mÞ 2s hypothesis testing for proportion to establish that proteins with high promiscuity mostly have high % of charged and polar and has m = 0.28, s = 1.1 and reduced chi-square of 3.0E-5. residues within a radius of 3 A˚ from the catalytic residues, and that Pearson’s chi-squared tests the goodness of fit of an observed proteins which have low % of charged and polar residues (again frequency distribution with a theoretical distribution, and a within a 3 A˚ radius) are rarely promiscuous. reduced chi-square statistic is the Pearson’s chi-squared divided PROMISE is an automated computational method for by the number of degrees of freedom.
Recommended publications
  • Ribonuclease L Mediates the Cell-Lethal Phenotype of the Double-Stranded RNA Editing
    1 2 Ribonuclease L mediates the cell-lethal phenotype of the double-stranded RNA editing 3 enzyme ADAR1 in a human cell line 4 5 Yize Lia,#, Shuvojit Banerjeeb,#, Stephen A. Goldsteina, Beihua Dongb, Christina Gaughanb, Sneha 6 Rathc, Jesse Donovanc, Alexei Korennykhc, Robert H. Silvermanb,* and Susan R Weissa,* 7 aDepartment of Microbiology, Perelman School of Medicine, University of Pennsylvania, 8 Philadelphia, PA, USA, 19104; b Department of Cancer Biology, Lerner Research Institute, 9 Cleveland Clinic, Cleveland, OH, USA 44195; c Department of Molecular Biology, Princeton 10 University, Princeton, NJ 08544 11 12 13 14 15 16 17 18 # These authors contributed equally to this work 19 * Corresponding authors 20 21 22 23 24 25 26 27 28 29 30 Abstract 31 ADAR1 isoforms are adenosine deaminases that edit and destabilize double-stranded RNA 32 reducing its immunostimulatory activities. Mutation of ADAR1 leads to a severe neurodevelopmental 33 and inflammatory disease of children, Aicardi-Goutiéres syndrome. In mice, Adar1 mutations are 34 embryonic lethal but are rescued by mutation of the Mda5 or Mavs genes, which function in IFN 35 induction. However, the specific IFN regulated proteins responsible for the pathogenic effects of 36 ADAR1 mutation are unknown. We show that the cell-lethal phenotype of ADAR1 deletion in human 37 lung adenocarcinoma A549 cells is rescued by CRISPR/Cas9 mutagenesis of the RNASEL gene or 38 by expression of the RNase L antagonist, murine coronavirus NS2 accessory protein. Our result 39 demonstrate that ablation of RNase L activity promotes survival of ADAR1 deficient cells even in the 40 presence of MDA5 and MAVS, suggesting that the RNase L system is the primary sensor pathway 41 for endogenous dsRNA that leads to cell death.
    [Show full text]
  • Chapter 1 Introduction
    The Conditional Protein Splicing of Alpha-Sarcin: A model for inducible assembly of protein toxins in vivo. by Spencer C. Alford B.Sc., University of Victoria, 2004 A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of MASTER OF SCIENCE in the Department of Biochemistry and Microbiology Spencer C. Alford, 2007 University of Victoria All rights reserved. This thesis may not be reproduced in whole or in part, by photocopy or other means, without the permission of the author. ii Supervisory Committee The Conditional Protein Splicing of Alpha-Sarcin: A model for inducible assembly of protein toxins in vivo. by Spencer C. Alford B.Sc, University of Victoria, 2004 Supervisory Committee Dr. Perry Howard, Supervisor (Department of Biochemistry and Microbiology) Dr. Juan Ausio, Departmental Member (Department of Biochemistry and Microbiology) Dr. Robert Chow, Outside Member (Department of Biology) iii Supervisory Committee Dr. Perry Howard, Supervisor (Department of Biochemistry and Microbiology) Dr. Juan Ausio, Departmental Member (Department of Biochemistry and Microbiology) Dr. Robert Chow, Outside Member (Department of Biology) Abstract Conditional protein splicing (CPS) is an intein-mediated post-translational modification. Inteins are intervening protein elements that autocatalytically excise themselves from precursor proteins to ligate flanking protein sequences, called exteins, with a native peptide bond. Artificially split inteins can mediate the same process by splicing proteins in trans, when intermolecular reconstitution of split intein fragments occurs. An established CPS model utilizes an artificially split Saccharomyces cerevisiae intein, called VMA. In this model, VMA intein fragments are fused to the heterodimerization domains, FKBP and FRB, which selectively form a complex with the immunosuppressive drug, rapamycin.
    [Show full text]
  • Yeast Genome Gazetteer P35-65
    gazetteer Metabolism 35 tRNA modification mitochondrial transport amino-acid metabolism other tRNA-transcription activities vesicular transport (Golgi network, etc.) nitrogen and sulphur metabolism mRNA synthesis peroxisomal transport nucleotide metabolism mRNA processing (splicing) vacuolar transport phosphate metabolism mRNA processing (5’-end, 3’-end processing extracellular transport carbohydrate metabolism and mRNA degradation) cellular import lipid, fatty-acid and sterol metabolism other mRNA-transcription activities other intracellular-transport activities biosynthesis of vitamins, cofactors and RNA transport prosthetic groups other transcription activities Cellular organization and biogenesis 54 ionic homeostasis organization and biogenesis of cell wall and Protein synthesis 48 plasma membrane Energy 40 ribosomal proteins organization and biogenesis of glycolysis translation (initiation,elongation and cytoskeleton gluconeogenesis termination) organization and biogenesis of endoplasmic pentose-phosphate pathway translational control reticulum and Golgi tricarboxylic-acid pathway tRNA synthetases organization and biogenesis of chromosome respiration other protein-synthesis activities structure fermentation mitochondrial organization and biogenesis metabolism of energy reserves (glycogen Protein destination 49 peroxisomal organization and biogenesis and trehalose) protein folding and stabilization endosomal organization and biogenesis other energy-generation activities protein targeting, sorting and translocation vacuolar and lysosomal
    [Show full text]
  • Anti-Inflammatory Role of Curcumin in LPS Treated A549 Cells at Global Proteome Level and on Mycobacterial Infection
    Anti-inflammatory Role of Curcumin in LPS Treated A549 cells at Global Proteome level and on Mycobacterial infection. Suchita Singh1,+, Rakesh Arya2,3,+, Rhishikesh R Bargaje1, Mrinal Kumar Das2,4, Subia Akram2, Hossain Md. Faruquee2,5, Rajendra Kumar Behera3, Ranjan Kumar Nanda2,*, Anurag Agrawal1 1Center of Excellence for Translational Research in Asthma and Lung Disease, CSIR- Institute of Genomics and Integrative Biology, New Delhi, 110025, India. 2Translational Health Group, International Centre for Genetic Engineering and Biotechnology, New Delhi, 110067, India. 3School of Life Sciences, Sambalpur University, Jyoti Vihar, Sambalpur, Orissa, 768019, India. 4Department of Respiratory Sciences, #211, Maurice Shock Building, University of Leicester, LE1 9HN 5Department of Biotechnology and Genetic Engineering, Islamic University, Kushtia- 7003, Bangladesh. +Contributed equally for this work. S-1 70 G1 S 60 G2/M 50 40 30 % of cells 20 10 0 CURI LPSI LPSCUR Figure S1: Effect of curcumin and/or LPS treatment on A549 cell viability A549 cells were treated with curcumin (10 µM) and/or LPS or 1 µg/ml for the indicated times and after fixation were stained with propidium iodide and Annexin V-FITC. The DNA contents were determined by flow cytometry to calculate percentage of cells present in each phase of the cell cycle (G1, S and G2/M) using Flowing analysis software. S-2 Figure S2: Total proteins identified in all the three experiments and their distribution betwee curcumin and/or LPS treated conditions. The proteins showing differential expressions (log2 fold change≥2) in these experiments were presented in the venn diagram and certain number of proteins are common in all three experiments.
    [Show full text]
  • Supplementary Table S4. FGA Co-Expressed Gene List in LUAD
    Supplementary Table S4. FGA co-expressed gene list in LUAD tumors Symbol R Locus Description FGG 0.919 4q28 fibrinogen gamma chain FGL1 0.635 8p22 fibrinogen-like 1 SLC7A2 0.536 8p22 solute carrier family 7 (cationic amino acid transporter, y+ system), member 2 DUSP4 0.521 8p12-p11 dual specificity phosphatase 4 HAL 0.51 12q22-q24.1histidine ammonia-lyase PDE4D 0.499 5q12 phosphodiesterase 4D, cAMP-specific FURIN 0.497 15q26.1 furin (paired basic amino acid cleaving enzyme) CPS1 0.49 2q35 carbamoyl-phosphate synthase 1, mitochondrial TESC 0.478 12q24.22 tescalcin INHA 0.465 2q35 inhibin, alpha S100P 0.461 4p16 S100 calcium binding protein P VPS37A 0.447 8p22 vacuolar protein sorting 37 homolog A (S. cerevisiae) SLC16A14 0.447 2q36.3 solute carrier family 16, member 14 PPARGC1A 0.443 4p15.1 peroxisome proliferator-activated receptor gamma, coactivator 1 alpha SIK1 0.435 21q22.3 salt-inducible kinase 1 IRS2 0.434 13q34 insulin receptor substrate 2 RND1 0.433 12q12 Rho family GTPase 1 HGD 0.433 3q13.33 homogentisate 1,2-dioxygenase PTP4A1 0.432 6q12 protein tyrosine phosphatase type IVA, member 1 C8orf4 0.428 8p11.2 chromosome 8 open reading frame 4 DDC 0.427 7p12.2 dopa decarboxylase (aromatic L-amino acid decarboxylase) TACC2 0.427 10q26 transforming, acidic coiled-coil containing protein 2 MUC13 0.422 3q21.2 mucin 13, cell surface associated C5 0.412 9q33-q34 complement component 5 NR4A2 0.412 2q22-q23 nuclear receptor subfamily 4, group A, member 2 EYS 0.411 6q12 eyes shut homolog (Drosophila) GPX2 0.406 14q24.1 glutathione peroxidase
    [Show full text]
  • Inclusion of a Furin Cleavage Site Enhances Antitumor Efficacy
    toxins Article Inclusion of a Furin Cleavage Site Enhances Antitumor Efficacy against Colorectal Cancer Cells of Ribotoxin α-Sarcin- or RNase T1-Based Immunotoxins Javier Ruiz-de-la-Herrán 1, Jaime Tomé-Amat 1,2 , Rodrigo Lázaro-Gorines 1, José G. Gavilanes 1 and Javier Lacadena 1,* 1 Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid 28040, Spain; [email protected] (J.R.-d.-l.-H.); [email protected] (J.T.-A.); [email protected] (R.L.-G.); [email protected] (J.G.G.) 2 Centre for Plant Biotechnology and Genomics (UPM-INIA), Universidad Politécnica de Madrid, Pozuelo de Alarcón, Madrid 28223, Spain * Correspondence: [email protected]; Tel.: +34-91-394-4266 Received: 3 September 2019; Accepted: 10 October 2019; Published: 12 October 2019 Abstract: Immunotoxins are chimeric molecules that combine the specificity of an antibody to recognize and bind tumor antigens with the potency of the enzymatic activity of a toxin, thus, promoting the death of target cells. Among them, RNases-based immunotoxins have arisen as promising antitumor therapeutic agents. In this work, we describe the production and purification of two new immunoconjugates, based on RNase T1 and the fungal ribotoxin α-sarcin, with optimized properties for tumor treatment due to the inclusion of a furin cleavage site. Circular dichroism spectroscopy, ribonucleolytic activity studies, flow cytometry, fluorescence microscopy, and cell viability assays were carried out for structural and in vitro functional characterization. Our results confirm the enhanced antitumor efficiency showed by these furin-immunotoxin variants as a result of an improved release of their toxic domain to the cytosol, favoring the accessibility of both ribonucleases to their substrates.
    [Show full text]
  • A Novel Mouse Model for Inhibition of DOHH-Mediated Hypusine Modification Reveals a Crucial Function in Embryonic Development, P
    © 2014. Published by The Company of Biologists Ltd | Disease Models & Mechanisms (2014) 7, 963-976 doi:10.1242/dmm.014449 RESEARCH ARTICLE A novel mouse model for inhibition of DOHH-mediated hypusine modification reveals a crucial function in embryonic development, proliferation and oncogenic transformation Henning Sievert1, Nora Pällmann1,2,*, Katharine K. Miller3,*, Irm Hermans-Borgmeyer3, Simone Venz4, Ataman Sendoel5,6, Michael Preukschas1, Michaela Schweizer3, Steffen Boettcher6, P. Christoph Janiesch3, Thomas Streichert7, Reinhard Walther4, Michael O. Hengartner5, Markus G. Manz6, Tim H. Brümmendorf8, Carsten Bokemeyer1, Melanie Braig1, Joachim Hauber2, Kent E. Duncan3 and Stefan Balabanov1,6,‡ ABSTRACT factor 5A (eIF5A), represents an essential mechanism in the control The central importance of translational control by post-translational of proliferation of eukaryotic cells (Cooper et al., 1982). This modification has spurred major interest in regulatory pathways that modification leads to the activation of eIF5A and is mediated by control translation. One such pathway uniquely adds hypusine to deoxyhypusine synthase (DHS), which catalyses the transfer of a 4- eukaryotic initiation factor 5A (eIF5A), and thereby affects protein aminobutyl moiety of spermidine to the ε-amino group of Lys50 to synthesis and, subsequently, cellular proliferation through an form an intermediate residue, deoxyhypusine (Dhp50) (Park et al., unknown mechanism. Using a novel conditional knockout mouse 1981). Subsequently, deoxyhypusine hydroxylase (DOHH)
    [Show full text]
  • XXI Fungal Genetics Conference Abstracts
    Fungal Genetics Reports Volume 48 Article 17 XXI Fungal Genetics Conference Abstracts Fungal Genetics Conference Follow this and additional works at: https://newprairiepress.org/fgr This work is licensed under a Creative Commons Attribution-Share Alike 4.0 License. Recommended Citation Fungal Genetics Conference. (2001) "XXI Fungal Genetics Conference Abstracts," Fungal Genetics Reports: Vol. 48, Article 17. https://doi.org/10.4148/1941-4765.1182 This Supplementary Material is brought to you for free and open access by New Prairie Press. It has been accepted for inclusion in Fungal Genetics Reports by an authorized administrator of New Prairie Press. For more information, please contact [email protected]. XXI Fungal Genetics Conference Abstracts Abstract XXI Fungal Genetics Conference Abstracts This supplementary material is available in Fungal Genetics Reports: https://newprairiepress.org/fgr/vol48/iss1/17 : XXI Fungal Genetics Conference Abstracts XXI Fungal Genetics Conference Abstracts Plenary sessions Cell Biology (1-87) Population and Evolutionary Biology (88-124) Genomics and Proteomics (125-179) Industrial Biology and Biotechnology (180-214) Host-Parasite Interactions (215-295) Gene Regulation (296-385) Developmental Biology (386-457) Biochemistry and Secondary Metabolism(458-492) Unclassified(493-502) Index to Abstracts Abstracts may be cited as "Fungal Genetics Newsletter 48S:abstract number" Plenary Abstracts COMPARATIVE AND FUNCTIONAL GENOMICS FUNGAL-HOST INTERACTIONS CELL BIOLOGY GENOME STRUCTURE AND MAINTENANCE COMPARATIVE AND FUNCTIONAL GENOMICS Genome reconstruction and gene expression for the rice blast fungus, Magnaporthe grisea. Ralph A. Dean. Fungal Genomics Laboratory, NC State University, Raleigh NC 27695 Rice blast disease, caused by Magnaporthe grisea, is one of the most devastating threats to food security worldwide.
    [Show full text]
  • Very Fast Empirical Prediction and Rationalization of Protein Pka Values
    PROTEINS: Structure, Function, and Bioinformatics 61:704–721 (2005) Very Fast Empirical Prediction and Rationalization of Protein pKa Values Hui Li,1† Andrew D. Robertson,2 and Jan H. Jensen1* 1Department of Chemistry and Center for Biocatalysis and Bioprocessing, The University of Iowa, Iowa City, Iowa 2Department of Biochemistry, The University of Iowa, Iowa City, Iowa 8,14–17 ABSTRACT A very fast empirical method is added to a model pKa value. These models usually presented for structure-based protein pKa predic- have a root-mean-square deviation (RMSD) from experi- tion and rationalization. The desolvation effects ment of Ͻ1. However, the currently available LPBE and intra-protein interactions, which cause varia- models tend to overestimate the intraprotein charge– 18,19 tions in pKa values of protein ionizable groups, are charge interactions and underestimate the hydrogen- 20 empirically related to the positions and chemical bonding and desolvation effects in calculating pKa shifts. nature of the groups proximate to the pKa sites. A The LPBE methods typically require tens of minutes to computer program is written to automatically pre- hours of computer time. dict pKa values based on these empirical relation- The pKa values of small organic acids and bases can be 21 ships within a couple of seconds. Unusual pKa val- predicted with the Hammet-Taft method based on empiri- ues at buried active sites, which are among the most cal relationship between pKa shifts and substituents. The interesting protein pKa values, are predicted very Hammet-Taft method is accurate for molecules similar to well with the empirical method. A test on 233 car- the ones included in the parameterization and is very fast.
    [Show full text]
  • Method for Prediction of Protein Function from Sequence Using The
    Article No. mb981993 J. Mol. Biol. (1998) 281, 949±968 Method for Prediction of Protein Function from Sequence using the Sequence-to-Structure-to- Function Paradigm with Application to Glutaredoxins/ Thioredoxins and T1 Ribonucleases Jacquelyn S. Fetrow1 and Jeffrey Skolnick2* 1Department of Biological The practical exploitation of the vast numbers of sequences in the gen- Sciences, Center for ome sequence databases is crucially dependent on the ability to identify Biochemistry and Biophysics the function of each sequence. Unfortunately, current methods, including University at Albany, SUNY global sequence alignment and local sequence motif identi®cation, are 1400 Washington Avenue limited by the extent of sequence similarity between sequences of Albany, NY 12222, USA unknown and known function; these methods increasingly fail as the sequence identity diverges into and beyond the twilight zone of sequence 2Department of Molecular identity. To address this problem, a novel method for identi®cation of Biology, The Scripps Institute protein function based directly on the sequence-to-structure-to-function 10550 North Torrey Pines paradigm is described. Descriptors of protein active sites, termed ``fuzzy Road, La Jolla, CA 92037, USA functional forms'' or FFFs, are created based on the geometry and confor- mation of the active site. By way of illustration, the active sites respon- sible for the disul®de oxidoreductase activity of the glutaredoxin/ thioredoxin family and the RNA hydrolytic activity of the T1 ribonuclease family are presented. First, the FFFs are shown to correctly identify their corresponding active sites in a library of exact protein models produced by crystallography or NMR spectroscopy, most of which lack the speci- ®ed activity.
    [Show full text]
  • Mimicking the Action of Ribonucleases: Studies on Rnase a and Design of PNA Based Artificial Enzymes
    From DEPARTMENT OF BIOSCIENCES AND NUTRITION Karolinska Institutet, Stockholm, Sweden Mimicking the action of ribonucleases: studies on RNase A and design of PNA based artificial enzymes Alice Ghidini Stockholm 2015 All previously published papers were reproduced with permission from the publisher. Published by Karolinska Institutet. Printed by Eprint AB 2015 © Alice Ghidini, 2015 ISBN 978-91-7676-039-0 Mimicking the action of ribonucleases: studies on RNase A and design of PNA based artificial enzymes THESIS FOR DOCTORAL DEGREE (Ph.D.) By Alice Ghidini Principal Supervisor: Opponent: Professor Roger Strömberg Professor Michael J. Gait Karolinska Institutet Medical Research Council Department of Bioscience and nutrition Department of Laboratory of Molecular Biology Division of Bioorganic chemistry Examination Board: Co-supervisor(s): Lars Baltzer Merita Murtola Uppsala Universitet Karolinska Institutet Department of Institutionen för kemi Department of Bioscience and nutrition Division of BMC, Fysikalisk organisk kemi Division of Bioorganic chemistry Mikael Leijon Malgorzata Honcharenko National Veterinary Institute (SVA) Karolinska Institutet Division of Virology, Immunobiology and Department of Bioscience and nutrition Parasitology (VIP) Division of Bioorganic chemistry Marcus Wilhelmsson Chalmers University of Technology Department of Chemical and Chemical Engineering/Chemistry and Biochemistry With love to my family, lontani ma sempre vicini… ABSTRACT A 3’-deoxy-3’-C-methylenephosphonate modified diribonucleotide is highly resistant to degradation by spleen phosphodiesterase and not cleaved at all by snake venom phosphodiesterase. Despite that both the vicinal 2-hydroxy nucleophile and the 5’-oxyanion leaving group are intact, the 3’-methylenephosponate RNA modification is also highly resistant towards the action of RNase A. Several different approaches were explored for conjugation of oligoethers to PNA with internally or N-terminal placed diaminopropionic acid residues.
    [Show full text]
  • Bacterial Retrons Function in Anti-Phage Defense
    Article Bacterial Retrons Function In Anti-Phage Defense Graphical Abstract Authors Adi Millman, Aude Bernheim, Retrons appear in an operon Retrons generate an RNA-DNA hybrid via reverse transcription with additional “effector” genes Avigail Stokar-Avihail, ..., Azita Leavitt, ncRNA Reverse Transcriptase (RT) msDNA Ribosyltransferase Yaara Oppenheimer-Shaanan, (RNA-DNA hybrid) DNA-binding RT Rotem Sorek RNA Retron function 2 transmembrane 5’ G RT domains 2’-5’ 3’ was unknown 3’ RT Correspondence cDNA RT Cold-shock [email protected] 5’ G 2’ 3’ G G cDNA RT RT ATPase Nuclease In Brief reverse transcription RNase H Retrons are part of a large family of anti- Retrons protect bacteria from phage Inhibition of RecBCD by phages triggers retron Ec48 defense phage defense systems that are widespread in bacteria and confer Bacterial density during phage infection Growth resistance against a broad range of B Effector arrest RT D RT activation with retron C B phages, mediated by abortive infection. Effector no retron Ec48 “guards” the bacterial Effector activation leads RecBCD complex to abortive infection bacterial density Retron Ec48 RT 2TM time RecBCD inhibitor B B D D RT Bacteria without retron Bacteria with retron C C B B Phage proteins inhibit RecBCD Retron Ec48 senses RecBCD inhibition Highlights d Retrons are preferentially located in defense islands d Retrons, together with their effector genes, protect bacteria from phages d Protection from phage is mediated by abortive infection d Retron Ec48 guards RecBCD. Inhibition of RecBCD by phages triggers retron defense Millman et al., 2020, Cell 183, 1–11 December 10, 2020 ª 2020 Elsevier Inc.
    [Show full text]