Prediction of Distal Residue Participation in Enzyme Catalysis

Total Page:16

File Type:pdf, Size:1020Kb

Prediction of Distal Residue Participation in Enzyme Catalysis Article Protein Science DOI 10.1002/pro.2648 Prediction of Distal Residue Participation in Enzyme Catalysis Heather R. Brodkina,b,1 Nicholas A. DeLateura,2 Srinivas Somarowthua,3 Caitlyn L. Millsa Walter R. Novakb,4 Penny J. Beuninga Dagmar Ringeb Mary Jo Ondrechena* aDepartment of Chemistry and Chemical Biology, Northeastern University, Boston, MA, 02115, USA, and bDepartments of Biochemistry and Chemistry, Rosenstiel Basic Medical Sciences Research Center, Brandeis University, Waltham, MA 02454-9110, USA1Present address: Alkermes, Inc., 852 Winter Street, Waltham, MA 02451, USA 2Present address: Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139 3Present address: Molecular, Cellular and Developmental Biology, Yale University, New Haven, CT 06520 4Present address: Department of Chemistry, Wabash College, Crawfordsville, IN 47933 *To whom correspondence should be sent. Corresponding author information: Professor Mary Jo Ondrechen Department of Chemistry and Chemical Biology Northeastern University 360 Huntington Avenue Boston, MA 02115 USA Telephone: +1-617-373-2856 Fax: +1-617-373-8795 E-mail: [email protected] This article has been accepted for publication and undergone full peer review but has not been through the copyediting, typesetting, pagination and proofreading process which may lead to differences between this version and the Version of Record. Please cite this article as doi: 10.1002/pro.2648 © 2014 The Protein Society Received: Nov 22, 2014; Revised: Jan 10, 2015; Accepted: Jan 26, 2015 This article is protected by copyright. All rights reserved. Protein Science Page 2 of 42 Abstract A scoring method for the prediction of catalytically important residues in enzyme structures is presented and used to examine the participation of distal residues in enzyme catalysis. Scores are based on the Partial Order Optimum Likelihood (POOL) machine learning method, using computed electrostatic properties, surface geometric features, and information obtained from the phylogenetic tree as input features. Predictions of distal residue participation in catalysis are compared with experimental kinetics data from the literature on variants of the featured enzymes; some additional kinetics measurements are reported for variants of Pseudomonas putida nitrile hydratase (ppNH) and for E. coli alkaline phosphatase (AP). The multi-layer active sites of Pseudomonas putida nitrile hydratase and of human phosphoglucose isomerase are predicted by the POOL logZP scores, as is the single-layer active site of Pseudomonas putida ketosteroid isomerase. The logZP score cutoff utilized here results in over- prediction of distal residue involvement in E. coli alkaline phosphatase. While fewer experimental data points are available for Pseudomonas putida mandelate racemase and for human carbonic anhydrase II, the POOL logZP scores properly predict the previously reported participation of distal residues. Keywords Remote residues; catalytic efficiency; Partial Order Optimum Likelihood; nitrile hydratase; phosphoglucose isomerase; ketosteroid isomerase; alkaline phosphatase Abbreviations CSA – Catalytic Site Atlas INTREPID Information Theoretic Tree Traversal for Protein Functional Site Identification KSI – Ketosteroid isomerase PGI – Phosphoglucose isomerase POOL – Partial Order Optimum Likelihood THEMATICS – Theoretical Microscopic Anomalous Titration Curve Shapes 2 John Wiley & Sons This article is protected by copyright. All rights reserved. Page 3 of 42 Protein Science Introduction Many chemical reactions that require high temperature and/or extreme pH in the laboratory are catalyzed by enzymes at physiological temperature and neutral pH. One of the central challenges in biochemistry is to understand how enzymes achieve this catalytic power under mild conditions. Structural and biochemical studies have now characterized the active sites of hundreds of enzymes; nearly all of these studies have focused on the amino acid residues in direct contact with the reacting substrate molecule(s). A typical enzyme might consist of hundreds of residues, but only a handful of them have been identified in the literature as catalytically important. In the manually curated portion of the Catalytic Site Atlas (CSA),1 a database of enzyme active site composition compiled from the experimental literature, an average of only 1.2% of the residues in the biological units of these 170 enzymes is listed as catalytically important. Virtually all of these reported catalytic residues are, in the active conformation of the enzyme, in direct contact with the reacting substrate molecule(s). A recent review surveys examples from the literature of proteins where distal residues have been shown to participate directly in the biochemical function.2 The present paper presents a model, based on computed chemical properties, that successfully predicts both first-shell and more distant residues that are important for function. Model predictions are compared with available experimental data and, for a few cases, new experimental results are presented. The ability to predict, with a simple calculation, the involvement of distal residues in catalytic and other processes is a very important capability for protein engineering and for understanding how nature builds catalytic sites. Here new computational and experimental results, together with experimental data from the literature, are presented to show that residues that are not in direct contact 3 John Wiley & Sons This article is protected by copyright. All rights reserved. Protein Science Page 4 of 42 with the reacting substrate molecule often play significant roles in catalysis, but more importantly, that their participation in catalysis is predictable with a simple calculation. Previously we have shown3,4 that all of the residues in a protein structure may be assigned scores based on computed chemical properties and that these scores may be used to rank-order all of the residues according to their probability of functional importance. In the present work, these scores are transformed so that they are more comparable across proteins of different sizes. We now show that the pattern of these scores predicts whether an enzyme has a highly localized, compact active site, or a spatially extended, multi-layer active site with extensive involvement of distal residues in catalysis and ligand binding. For purposes of the present study, residues interacting directly with the substrate molecule are referred to as first-shell residues; other residues interacting directly with one or more first-shell residues are called second-shell residues; and finally, other residues interacting with one or more second-shell residues are called third-shell residues. All residues defined as first-, second- or third-shell thus may be thought of as members of respective spheres surrounding the reacting substrate molecule. In the present work, enzymes with available experimental data on the participation of distal residues in catalysis are the focus. THEMATICS (THEoretical Microscopic Anomalous TItration Curve Shapes) is a computational method for the prediction of catalytic and binding residues in proteins solely from the 3D structure of the query protein.5-7 THEMATICS takes advantage of the unique chemical and electrostatic properties of catalytic and binding residues. In particular, the ionizable residues in active sites tend to have anomalously shaped theoretical titration curves. We have established metrics6 to quantify the degree of deviation from normal titration behavior and have shown how these metrics may be used 4 John Wiley & Sons This article is protected by copyright. All rights reserved. Page 5 of 42 Protein Science to predict, with high sensitivity and specificity,7 the active site residues in a protein 3D structure, just from a simple set of calculations. POOL (Partial Order Optimum Likelihood)3,4,8 is a machine learning method for the prediction of functionally important residues in protein structures. POOL, a multidimensional, monotonicity-constrained maximum likelihood technique, utilizes input for which the output variable (in this case the probability of functional importance of a residue) is a monotonic function of each of the input features. POOL starts with the hypothesis that the larger the THEMATICS metrics for a given residue, the higher the probability that that residue is important for the biochemical function. Electrostatics features from THEMATICS are combined with multidimensional isotonic regression to form maximum likelihood estimates of probabilities that specific residues are active in biochemical function. All residues are then rank-ordered according to their probability of functional importance. Through the introduction of environment variables in POOL, THEMATICS metrics are used to characterize the chemical environments of all residues, not just the ionizable ones. Furthermore, the input to POOL is flexible, in that any property upon which the probability of the functional importance of a residue depends monotonically can be a POOL input feature. These include features based solely on the structure of the query protein, such as THEMATICS metrics6,7,9 and surface geometric properties,10 but may also include other features if they are available, such as phylogenetic tree-based scores.11 It has been demonstrated3 that POOL with purely structure-based features performs very nearly as well as the best methods that utilize both sequence alignments and structure. POOL with input features that include
Recommended publications
  • Computational Protein Function Predictions
    Methods 93 (2016) 1–2 Contents lists available at ScienceDirect Methods journal homepage: www.elsevier.com/locate/ymeth Guest Editor’s Introduction Computational protein function predictions What is the role of this protein? Where does this protein local- in CATH-Gene3D using a hidden Markov model and thereby ize in a cell? Are there any ligands that bind to this protein? If so, predict function of the query. what are they? Which residues constitute a functional site of the The next three methods [4–6] take a protein tertiary structure protein? These questions, which in a broader sense seek the biolog- as input and predict binding ligands for the query protein. Binding ical function of a protein, are fundamental and central in modern ligand prediction is not only for function prediction of proteins but biology. Ultimately, the biological function of a protein needs to it is also useful for drug design. The method developed by be determined by experiments. However, a hypothesis is needed Nakamura and Tomii represents a binding pocket by an ensemble to design an assay because it determines whether a target protein of all triangles consisting of three Ca atoms in the pocket [4]. The has a particular function or not. Biologists come up with hypothe- triangles are classified into 171,700 types considering amino acid ses of protein function from circumstantial evidence, and computa- types of vertices and edge distances. Thus, a pocket is represented tional function prediction can play an important part. by a vector that shows the frequency of each triangle type in the Computational function prediction methods are also useful for pocket, which is further reduced to a vector of 11 elements by analyzing protein function in a proteomic scale since methods multi-dimensional scaling.
    [Show full text]
  • Bioorganic & Medicinal Chemistry, Volume 11, Dec, 2003
    MMCC Results Volume 14 No. 3 1 April 2005 Vol. 14, No. 3 1 April 2005 Editorial and News A groundbreaker in membrane biophysics was Coverage period: 1 Mar. 2005 through 1 April 2005 published by Gunther Peters and colleagues this month in About 253 papers from 13 journals are cited. J. Chem. Physics (see Methodology, Molecular Dynamics). It deals with a problem in membrane 1. APPLICATIONS (105) . page 2 biophysics that has achieved nearly universal interest in the 1.1 Small Molecules – 29 past year or so, the computation of the lateral pressure profile in the lipid bilayer. The reason this is so interesting Water and Solvation – 6 QSAR — 12 Med Chem & Drug Design - 9 Carbon Nanoparticles - 2 is the heterogeneous and variable shapes of membrane proteins. Although we still don’t know much about them, 1.2 Biopolymers – 73 membrane NMR and spectroscopic studies of membrane Bioinformatics – 6 Protein Dynamics – 5 proteins are starting to reveal small but significant Protein Seq Anal and Align – 3 Ligand Binding – 4 variations in countour at their interfaces with lipid Protein Structure Prediction – 2 Enzyme Catalysis – 8 Comp and Homol Modeling – 3 Protein-Protein Inter. – 2 membranes and large conformational changes during Peptide Conform Anal – 1 Membrane Proteins - 9 function. Although most membrane proteins are essentially Protein Structure Anal – 8 Protein Nuc Acid Inter - 1 Protein Folding – 6 Nucleic Acids – 4 right cylinders (either β-barrels or helix bundles) with axes Protein Design& Engineering – 2 Lipids and Surfactants - 4 parallel to the membrane normal, border helices or strands Protein Hydration – 3 Carbohydrates - 1 can sometimes have tilts, kinks, or extrusions that require Protein Electr and Titr – 1 accommodation by bilayer lipids.
    [Show full text]
  • Application of Molecular Modeling to Drug Discovery and Functional Genomics
    APPLICATION OF MOLECULAR MODELING TO DRUG DISCOVERY AND FUNCTIONAL GENOMICS A Dissertation Presented By Zhouxi (Josie) Wang To The Department of Chemistry and Chemical Biology in partial fulfillment of the requirements For the degree of Doctor of Philosophy in the field of Chemistry Northeastern University Boston, Massachusetts May, 2012 1 © 2012 Zhouxi Wang ALL RIGHTS RESERVED 2 APPLICATION OF MOLECULAR MODELING TO DRUG DISCOVERY AND FUNCTIONAL GENOMICS By Zhouxi (Josie) Wang ABSTRACT OF DISSERTATION Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Chemistry in the Department of Chemistry and Chemical Biology in the Graduate School of Northeastern University, Boston, Massachusetts May 2012 3 ABSTRACT Molecular modeling can accelerate and guide drug design and contribute to the understanding of the biochemical functions of gene products. This thesis applies molecular modeling to facilitate the drug design for human African trypanosomiasis and develops a new modeling technique for protein biochemical function annotations. A special technique employed in this work is the prediction of the individual amino acids in protein structures that are involved in ligand interactions; these predicted local interaction sites are used for drug discovery and for the prediction of the biochemical function of proteins of unknown function. Chapter 2 applies molecular modeling to the structure based drug design for human African trypanosomiasis (HAT) at the Aurora kinase -1 target. HAT is a vector borne disease caused by several species of trypanosomes, affecting thousands of people every year. This disease is fatal if untreated. Current therapeutic interventions are unsatisfactory, all with limited efficacy or life-threatening side effects.
    [Show full text]
  • Specificity of DNA Damage Inducible DNA Polymerase IV from Escherichia Coli
    Specificity of DNA Damage Inducible DNA Polymerase IV from Escherichia coli A thesis presented by Jason M. Walsh to The Department of Chemistry and Chemical Biology in partial fulfillment of the requirements for the degree of Master of Science in the field of Chemistry Northeastern University Boston, Massachusetts December 2009 1 © 2009, Jason M Walsh All rights reserved. 2 Specificity of DNA Damage Inducible DNA Polymerase IV from Escherichia coli by Jason M. Walsh ABSTRACT OF THESIS Submitted in partial fulfillment of the requirements for the degree of Master of Science in Chemistry and Chemical Biology in the Graduate School of Arts and Sciences of Northeastern University, December 2009 3 ABSTRACT DNA polymerases are responsible for DNA replication during cell division. There are multiple families of polymerases (A, B, C, D, X) responsible for copying DNA during replication and repair. There is also a class of polymerases conserved throughout evolution, known as the Y family polymerases, that have reduced replication fidelity on undamaged DNA (Tang et al. 2000). However Y family DNA polymerases have the specialized property of replicating DNA by copying damaged DNA, a process known as translesion synthesis (TLS). Structural differences between Y family and replicative polymerases may account for the difference in enzymatic activity. However we demonstrate that the Klenow fragment (A family) can bypass a fluorescent cytosine analog known as 1, 3-diaza-2-oxophenothiazine (tC), that DinB, a Y family polymerase, cannot bypass. We show that DinB inserts dGTP faithfully, but cannot extend the DNA primer beyond that. Verifying which amino acid residues are responsible for both function and specificity of the Y family polymerases is accomplished by assessing the kinetic data of nucleotide incorporation events of DinB variants as compared to wild-type DinB.
    [Show full text]
  • Identification of Functional Subclasses in the DJ-1 Superfamily Proteins
    University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Biochemistry -- Faculty Publications Biochemistry, Department of 1-2007 Identification of unctionalF Subclasses in the DJ-1 Superfamily Proteins Ying Wei Northeastern University, Dagmar Ringe University of Nebraska - Lincoln Mark A. Wilson Brandeis University, [email protected] Mary Jo Ondrechen Northeastern University, Follow this and additional works at: https://digitalcommons.unl.edu/biochemfacpub Part of the Biochemistry, Biophysics, and Structural Biology Commons Wei, Ying; Ringe, Dagmar; Wilson, Mark A.; and Ondrechen, Mary Jo, "Identification of unctionalF Subclasses in the DJ-1 Superfamily Proteins" (2007). Biochemistry -- Faculty Publications. 25. https://digitalcommons.unl.edu/biochemfacpub/25 This Article is brought to you for free and open access by the Biochemistry, Department of at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Biochemistry -- Faculty Publications by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. Identification of Functional Subclasses in the DJ-1 Superfamily Proteins Ying Wei1,2, Dagmar Ringe3,4,5*, Mark A. Wilson3,4,5¤, Mary Jo Ondrechen1,2 1 Department of Chemistry and Chemical Biology, Northeastern University, Boston, Massachusetts, United States of America, 2 Institute for Complex Scientific Software, Northeastern University, Boston, Massachusetts, United States of America, 3 Department of Biochemistry, Brandeis University, Waltham, Massachusetts, United States of America, 4 Department of Chemistry, Brandeis University, Waltham, Massachusetts, United States of America, 5 Rosenstiel Basic Medical Sciences Research Center, Brandeis University, Waltham, Massachusetts, United States of America Genomics has posed the challenge of determination of protein function from sequence and/or 3-D structure. Functional assignment from sequence relationships can be misleading, and structural similarity does not necessarily imply functional similarity.
    [Show full text]
  • Che Newsletter
    SUMMER 2020 DIVISION OF CHEMISTRY CHE NEWSLETTER INSIDE CHE FUNDED RESEARCH 2 Update from the Division Director HIGHLIGHT Making single molecule electronic components, like rectifiers or switches, is 3 Chemistry Staffing News and Recruiting an ultimate goal for fabricating electronics a. Special Thank-You to Carol Bessel in the nano realm. Traditionally, single b. Congratulations to Dr. Evelyn Goldfield This work was published molecule junctions are fabricated by on her Retirement in Angewandte Chemie anchoring the molecule between two metal (https://doi.org/10.1002/ c. Welcome the new Rotators and Farewell electrodes and charge transport is measured to the departing Rotators anie.201903898, CHE-1508567). only along a single axis. When studying 8 Recruiting Rotators: Would you like to join the charge conduction through benzene the NSF team as a rotator? derivatives, Eric Borguet’s group at Temple University discovered that the geometry of a molecule between the electrodes has a significant impact 8 CHE COVID-19 Effort on the charge transport through that molecule. Harnessing that anisotropic effect in trimesic acid, Borguet’s group developed a single‐molecule 10 CHE Virtual Office Hours electromechanical switch. Using an applied voltage at the electrode, they controlled the geometry of the molecule and found that a planar orientation 11 Funding Opportunities gave rise to a conductance that is more than 400 times higher than that of a. Division of Chemistry’s updated Disciplinary Research Programs Solicitation (CHE-DRP) an upright molecule. b. MPS Alliances for Graduate Education and the Professoriate (AGEP) – Graduate Research Supplements (GRS) Dear Colleague Letter 12 Chemistry Outreach in 2020 POLARITON CHEMISTRY: 12 2021 NSF Chemistry Early Career Investigator Workshop EXPLOITING CONFINED LIGHT 12 Chemistry Highlights FOR CHEMICAL Figure from J.B.
    [Show full text]
  • Selective Prediction of Interaction Sites in Protein Structures with THEMATICS Ying Wei1, Jaeju Ko1,2, Leonel F Murga1,3 and Mary Jo Ondrechen*1
    BMC Bioinformatics BioMed Central Research article Open Access Selective prediction of interaction sites in protein structures with THEMATICS Ying Wei1, Jaeju Ko1,2, Leonel F Murga1,3 and Mary Jo Ondrechen*1 Address: 1Department of Chemistry and Chemical Biology and Institute for Complex Scientific Software, Northeastern University, Boston, Massachusetts 02115 USA, 2NSF-ROA awardee. Department of Chemistry, Indiana University of Pennsylvania, 975 Oakland Avenue, Indiana, Pennsylvania 15705 USA and 3Rosenstiel Basic Medical Sciences Center, Brandeis University, Waltham, MA 02454 USA Email: Ying Wei - [email protected]; Jaeju Ko - [email protected]; Leonel F Murga - [email protected]; Mary Jo Ondrechen* - [email protected] * Corresponding author Published: 9 April 2007 Received: 1 November 2006 Accepted: 9 April 2007 BMC Bioinformatics 2007, 8:119 doi:10.1186/1471-2105-8-119 This article is available from: http://www.biomedcentral.com/1471-2105/8/119 © 2007 Wei et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Background: Methods are now available for the prediction of interaction sites in protein 3D structures. While many of these methods report high success rates for site prediction, often these predictions are not very selective and have low precision. Precision in site prediction is addressed using Theoretical Microscopic Titration Curves (THEMATICS), a simple computational method for the identification of active sites in enzymes. Recall and precision are measured and compared with other methods for the prediction of catalytic sites.
    [Show full text]
  • Using Computation to Identify Catalytic Features of Enzymes And
    Graduate Category: Physical and Life Sciences Degree Seeking: PhD Abstract ID# 2277 Using Computation to Identify Catalytic Features of Enzymes and to POOL Server: www.pool.neu.edu Engineer Enzymes with Enhanced or Altered Function DNA Lab Site: www.dna.neu.edu 205-R2-CX-0011 Timothy A. Coulther, Israel Adam, Penny J. Beuning, Mary Jo Ondrechen (TAC) NSF MCB-1158176 ORG Site: nuweb15.neu.edu/org NSF-MCB-1517290 Retro-aldolases: de novo Engineering a replicative polymerase for accurate damage bypass The catalytic properties of enzymes make them useful for wide application, THEoretical Microscopic Anomalous designed and evolved in Identify mutations that promote the accurate lesion conformation while such as in medicine, biofuel production, or sample analysis. However, we are TItration Curve Shapes1 laboratories3,4,5,6,7 retaining catalytic properties currently restricted to the enzymes that nature has to offer, limiting their Typical H-H Accurate Mutagenic applicability. Unfortunately, most initial enzyme design attempts Curve Identification of catalytic residues G stays G G becomes T have resulted in very low activities or no activity at all. Long-term directed evolution computationally for understanding activity and further engineering experiments have increased enzyme activities, but these require Approach years, high cost, and sometimes sophisticated equipment, that limit their development and use. A major issue is the sheer size of mutational space, as there Perturbed H-H 48 Charge are over 1.0x10 possibilities for a protein with 200 residues. The goal of our 4 Curves research is to use computation to reduce the number of possibilities, identifying a smaller number of potential designs that are more likely to be successful.
    [Show full text]
  • Ronald E. Hatcher Science on Saturday Lecture Series 24 January 2015
    Ronald E. Hatcher Science on Saturday Lecture Series 24 January 2015 Using Physics to Understand the Genome Mary Jo Ondrechen Department of Chemistry & Chemical Biology Northeastern University, Boston, MA 02115 USA ABSTRACT: The sequencing of the genome, the complete DNA, for humans and for thousands of other species at the dawn of the new millennium brought promises of untold advances in medicine and of many other novel technologies of benefit to the world. However, big questions remain about the roles played by these DNA sequences, including many of the protein-coding genes. Structural Genomics (SG) projects have since determined the 3D structures of over 13,000 gene products, but most of them are of unknown or uncertain function. In order to transform the promises of genome science into useful and practical advances, we must understand the function of these protein gene products. Therefore the functional classification of protein 3D structures is an important current problem. A new approach to this problem is described based primarily on chemical and electrostatic properties computed using techniques from chemical physics. Our function prediction method utilizes electrical potential functions that we calculate from the protein 3D structure. These potential functions are in turn used to predict the amino acids that are involved directly in the protein’s function. For proteins of known function, these sets of amino acids are used to establish “chemical signatures” for each functional type. Then, we predict the sets of active amino acids for protein structures of unknown function and, using matching with the chemical signatures, seek a functional assignment for the proteins of unknown function.
    [Show full text]
  • Machine Learning Differentiates Enzymatic and Non- Enzymatic Metals in Proteins ✉ Ryan Feehan 1,3, Meghan W
    ARTICLE https://doi.org/10.1038/s41467-021-24070-3 OPEN Machine learning differentiates enzymatic and non- enzymatic metals in proteins ✉ Ryan Feehan 1,3, Meghan W. Franklin1,3 & Joanna S. G. Slusky 1,2 Metalloenzymes are 40% of all enzymes and can perform all seven classes of enzyme reactions. Because of the physicochemical similarities between the active sites of metal- loenzymes and inactive metal binding sites, it is challenging to differentiate between them. fi 1234567890():,; Yet distinguishing these two classes is critical for the identi cation of both native and designed enzymes. Because of similarities between catalytic and non-catalytic metal binding sites, finding physicochemical features that distinguish these two types of metal sites can indicate aspects that are critical to enzyme function. In this work, we develop the largest structural dataset of enzymatic and non-enzymatic metalloprotein sites to date. We then use a decision-tree ensemble machine learning model to classify metals bound to proteins as enzymatic or non-enzymatic with 92.2% precision and 90.1% recall. Our model scores electrostatic and pocket lining features as more important than pocket volume, despite the fact that volume is the most quantitatively different feature between enzyme and non- enzymatic sites. Finally, we find our model has overall better performance in a side-to-side comparison against other methods that differentiate enzymatic from non-enzymatic sequences. We anticipate that our model’s ability to correctly identify which metal sites are responsible for enzymatic activity could enable identification of new enzymatic mechanisms and de novo enzyme design. 1 Center for Computational Biology, The University of Kansas, Lawrence, KS, USA.
    [Show full text]
  • THEMATICS: a Simple Computational Predictor of Enzyme Function from Structure
    THEMATICS: A simple computational predictor of enzyme function from structure Mary Jo Ondrechen*†‡, James G. Clifton†, and Dagmar Ringe†§ *Department of Chemistry, Northeastern University, Boston, MA 02115; and †Rosenstiel Basic Medical Sciences Research Center and §Departments of Biochemistry and Chemistry, Brandeis University, Waltham, MA 02454 Communicated by Gregory A. Petsko, Brandeis University, Waltham, MA, August 18, 2001 (received for review June 19, 2001) We show that theoretical microscopic titration curves (THEMATICS) sometimes be verified by spectroscopic measurements. Also, the can be used to identify active-site residues in proteins of known pKa values of the catalytic residues can be determined indirectly structure. Results are featured for three enzymes: triosephosphate by the measured reaction rate as a function of pH. This measured isomerase (TIM), aldose reductase (AR), and phosphomannose pKa can then be compared with the computed pKa’s of the isomerase (PMI). We note that TIM and AR have similar structures residues of the protein to identify the active site. We now show but catalyze different kinds of reactions, whereas TIM and PMI that the shapes of the theoretical titration functions carry have different structures but catalyze similar reactions. Analysis of information about active-site location, even in the absence of the theoretical microscopic titration curves for all of the ionizable experimental kinetic data. residues of these proteins shows that a small fraction (3–7%) of the The titration curves of most of the ionizable residues in a curves possess a flat region where the residue is partially proton- protein have a typical shape: as pH is increased, there is a ated over a wide pH range.
    [Show full text]
  • Form and Function in Enzyme Activity 6 April 2012, by Angela Herring
    Form and function in enzyme activity 6 April 2012, By Angela Herring are structurally important while others are required for the enzyme's catalytic properties. "There are cavities on the surface of a protein where a molecule can come in and sit down," Ondrechen said. "The enzyme does a reaction on it and the product goes away." The current body of research on enzyme activity mostly focuses on the amino acids in that cavity, which come into direct contact with the reactive molecule. But over the years, some research has Mary Jo Ondrechen and Penny Beuning recently suggested that amino acids far away from the received an NSF grant to study enzyme activity in the active site also play a role in catalysis. lab and computationally. Credit: Dominick Reuter Ondrechen's team, using a method she developed 10 years ago, will be able to predict which remote amino acids will impact reactivity. Beuning's team Many industrial chemistry applications, such as will test these predictions experimentally. drug or biofuel synthesis, require large energy inputs and often produce toxic pollutants. But "My lab is really interested in specificity of chemistry and chemical biology professor Mary Jo enzymes," Beuning said. "We look enzymes and Ondrechen said enzymes - proteins that increase figure out how they recognize their substrates." the rate of chemical reactions in the body - could be used to effectively replace standard industrial To do this, her team takes a protein engineering processes. approach in which they manipulate the enzyme's composition and observe how it affects its function. "Enzymatic reactions are cleaner, produce fewer byproducts and use less energy," she explained.
    [Show full text]