DOCSLIB.ORG

Download Parameters Files Containing All Essential Information Regarding Docking Procedures

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Download Parameters Files Containing All Essential Information Regarding Docking Procedures Modelling the structure and interactions of leukocyte integrins Kyle-Richard Dawson Submitted in fulfilment of the requirements for the degree of Magister Scientiae: Biochemistry (Msc) in the Faculty of Science at the Nelson Mandela University 9 April 2019 Prof Vaughan Oosthuizen i Declaration I know that plagiarism is wrong. Plagiarism is to use another’s work and pretend that it is one’s own. I have used the Bioinformatics convention for citation and referencing. Each contribution to, and quotation in, this dissertation from the work(s) of other people has been attributed, and has been cited and referenced. This dissertation is my own work. This work has not been submitted to any institution other than Nelson Mandela University. I have not allowed, and will not allow, anyone to copy my work with the intention of passing it off as his or her own work. Signature ______________________________ Date __________________________________ ii Acknowledgements I would like to express my sincere gratitude and appreciation to: My supervisor, Prof V Oosthuizen, for his positive attitude and guidance and the National Research Foundation (NRF) for financial assistance. I would also like to thank Dr R Hatherley, Dr Vuyani Moses and Prof O Bishop for their guidance in the field of Bioinformatics. I would also like to thank the following; Abigail Sephton, Liza Findt, Eugen Schnautz, Jan Batelka, Travis Dugmore, Martin Dorfling and Blake Callahan for the encouragement, occasional idea and helpful assistance in maintaining my composure. iii Contents Abbreviations and acronyms ................................................................................................................ vii List of Figures ......................................................................................................................................... ix List of Tables .......................................................................................................................................... xi Abstract ................................................................................................................................................... 1 1. Integrin proteins ............................................................................................................................. 2 1.1. Structure ................................................................................................................................. 2 1.2. Activation and regulation ........................................................................................................ 5 1.3. Functions ................................................................................................................................. 6 1.3.1. Cell signalling ................................................................................................................... 6 1.3.2. Cell survival ........................................................................................................................... 7 1.3.3. Cell proliferation ................................................................................................................... 7 1.3.4. Cell differentiation ................................................................................................................ 9 1.3.5. Leukocyte recruitment, activation and adhesion ................................................................. 9 1.3.6. Integrin ligands .................................................................................................................... 11 2. Homology modelling ..................................................................................................................... 13 2.1. Comparative homology modelling overview ........................................................................ 14 2.1.1. Template searching ....................................................................................................... 14 2.1.2. Selecting templates ....................................................................................................... 16 2.1.3. Sequence alignment ...................................................................................................... 17 2.1.4. Model construction ....................................................................................................... 20 2.2. Modelling programs .............................................................................................................. 23 2.2.1. MODELLER ..................................................................................................................... 23 2.2.2. MEDELLER ..................................................................................................................... 23 2.2.3. PRIMO ........................................................................................................................... 25 2.2.4. Phyre2 ........................................................................................................................... 25 2.2.5. I-TASSER ........................................................................................................................ 27 2.3. Errors in homology modelling ............................................................................................... 28 2.4. Determining model accuracy ................................................................................................ 28 2.4.1. PROSA............................................................................................................................ 28 2.4.2. Verify 3D ........................................................................................................................ 29 2.4.3. PROSESS ........................................................................................................................ 29 2.5. Docking programs ................................................................................................................. 31 2.5.1. HADDOCK2.2 ................................................................................................................. 31 2.5.2. CLUSPRO ....................................................................................................................... 33 2.5.3. AutoDock Vina ............................................................................................................... 35 iv 3. Problem statement, approach and objectives for this study ....................................................... 37 4. Methods ........................................................................................................................................ 41 4.1. Obtaining target sequences .................................................................................................. 41 4.2. Modelling using online servers ............................................................................................. 41 4.2.1. PRIMO ........................................................................................................................... 42 4.2.2. I-TASSER and Phyre2 ..................................................................................................... 42 4.2.3. MEDELLER ..................................................................................................................... 42 4.3. MODELLER ............................................................................................................................. 42 4.3.1. Generating alignment files, PDB files and python scripts ............................................. 43 4.3.2. Modelling the complete protein ................................................................................... 43 4.3.3. Determining optimal model iterations ......................................................................... 43 4.3.4. Determining optimal refinement iterations ................................................................. 44 4.3.5. Determining optimal template numbers ...................................................................... 44 4.3.6. Very_slow against slow_large refinement .................................................................... 44 4.4. Modelling monomeric subunits ............................................................................................ 44 4.4.1. Initial modelling ............................................................................................................ 44 4.4.2. Separating “closed” and “open” templates .................................................................. 45 4.4.3. Modelling with secondary structure arguments .......................................................... 45 4.4.4. Modelling with forced transmembrane regions ........................................................... 45 4.4.5. Fragmented integrin modelling .................................................................................... 45 4.5. Final modelling ...................................................................................................................... 45 4.5.1. BLASTp........................................................................................................................... 46 4.5.2. Separation of “closed” and “open” templates ............................................................. 46 4.5.3. Preparing multiple runs under varying template numbers .......................................... 46 4.5.4. Generating alignment files ...........................................................................................
Load more

Recommended publications
  • Accurate and Automated Classification of Protein Secondary Structure with Psicsi
    Accurate and automated classification of protein secondary structure with PsiCSI LING-HONG HUNG AND RAM SAMUDRALA Computational Genomics, Department of Microbiology, University of Washington, Seattle, Washington 98109, USA (RECEIVED July 2, 2002; FINAL REVISION October 18, 2002; ACCEPTED October 31, 2002) Abstract PsiCSI is a highly accurate and automated method of assigning secondary structure from NMR data, which is a useful intermediate step in the determination of tertiary structures. The method combines information from chemical shifts and protein sequence using three layers of neural networks. Training and testing was performed on a suite of 92 proteins (9437 residues) with known secondary and tertiary structure. Using a stringent cross-validation procedure in which the target and homologous proteins were removed from the databases used for training the neural networks, an average 89% Q3 accuracy (per residue) was observed. This is an increase of 6.2% and 5.5% (representing 36% and 33% fewer errors) over methods that use chemical shifts (CSI) or sequence information (Psipred) alone. In addition, PsiCSI improves upon the and is able to use sequence (87.4% ס translation of chemical shift information to secondary structure (Q3 without 13C shifts and 86.9% ס information as an effective substitute for sparse NMR data (Q3 with only H␣ shifts available). Finally, errors made by PsiCSI almost exclusively involve the 86.8% ס Q3 interchange of helix or strand with coil and not helix with strand (<2.5 occurrences per 10000 residues). The automation, increased accuracy, absence of gross errors, and robustness with regards to sparse data make PsiCSI ideal for high-throughput applications, and should improve the effectiveness of hybrid NMR/de novo structure determination methods.
    [Show full text]
  • A Web Server for Rapid NMR-Based Protein Structure Determination Berjanskii, M.; Tang, P.; Liang, J.; Cruz, J
    NRC Publications Archive Archives des publications du CNRC GeNMR: a web server for rapid NMR-based protein structure determination Berjanskii, M.; Tang, P.; Liang, J.; Cruz, J. A.; Zhou, J.; Zhou, Y.; Bassett, E.; Macdonell, C.; Lu, P.; Lin, G.; Wishart, D. S. This publication could be one of several versions: author’s original, accepted manuscript or the publisher’s version. / La version de cette publication peut être l’une des suivantes : la version prépublication de l’auteur, la version acceptée du manuscrit ou la version de l’éditeur. For the publisher’s version, please access the DOI link below./ Pour consulter la version de l’éditeur, utilisez le lien DOI ci-dessous. Publisher’s version / Version de l'éditeur: https://doi.org/10.1093/nar/gkp280 Nucleic Acids Research, 37, Suppl. 2, pp. W670-W677, 2009-07-01 NRC Publications Record / Notice d'Archives des publications de CNRC: https://nrc-publications.canada.ca/eng/view/object/?id=8c24f92f-5a63-4bc9-9510-868e720b942e https://publications-cnrc.canada.ca/fra/voir/objet/?id=8c24f92f-5a63-4bc9-9510-868e720b942e Access and use of this website and the material on it are subject to the Terms and Conditions set forth at https://nrc-publications.canada.ca/eng/copyright READ THESE TERMS AND CONDITIONS CAREFULLY BEFORE USING THIS WEBSITE. L’accès à ce site Web et l’utilisation de son contenu sont assujettis aux conditions présentées dans le site https://publications-cnrc.canada.ca/fra/droits LISEZ CES CONDITIONS ATTENTIVEMENT AVANT D’UTILISER CE SITE WEB. Questions? Contact the NRC Publications Archive team at [email protected].
    [Show full text]
  • Chemical Shift Secondary Structure Population Inference
    bioRxiv preprint doi: https://doi.org/10.1101/2021.02.20.432095; this version posted February 20, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. 1 CheSPI: Chemical shift Secondary structure Population Inference Jakob Toudahl Nielsen1* and Frans A.A. Mulder1* 1Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry, Aarhus University, Gustav Wieds VeJ 14, 8000 Aarhus C, Denmark *Correspondence to: [email protected], [email protected] ORCID: Jtn http://orcid.org/0000-0002-9554-4662 faam http://orcid.org/0000-0001-9427-8406 Keywords: NMR, protein, order, disorder, chemical shifts Abstract NMR chemical shifts (CSs) are delicate reporters of local protein structure, and recent advances in random coil CS (RCCS) prediction and interpretation now offer the compelling prospect of inferring small populations of structure from small deviations from RCCSs. Here, we present CheSPI, a simple and efficient method that provides unbiased and sensitive aggregate measures of local structure and disorder. It is demonstrated that CheSPI can predict even very small amounts of residual structure and robustly delineate subtle differences into four structural classes for intrinsically disordered proteins. For structured regions and proteins, CheSPI can assign up to eight structural classes, which coincide with the well-known DSSP classification. The program is freely available, and can either be invoked from URL www.protein-nmr.org as a web implementation, or run locally from command line as a python program.
    [Show full text]
  • CSI 3.0: a Web Server for Identifying Secondary and Super-Secondary Structure in Proteins Using NMR Chemical Shifts Noor E
    W370–W377 Nucleic Acids Research, 2015, Vol. 43, Web Server issue Published online 15 May 2015 doi: 10.1093/nar/gkv494 CSI 3.0: a web server for identifying secondary and super-secondary structure in proteins using NMR chemical shifts Noor E. Hafsa1, David Arndt1 and David S. Wishart1,2,* 1Department of Computing Science, University of Alberta, Edmonton, AB T6G 2E8, Canada and 2Department of Biological Sciences, University of Alberta, Edmonton, AB T6G 2E8, Canada Received March 10, 2015; Revised April 23, 2015; Accepted May 02, 2015 ABSTRACT ture elements has long been an integral part of the pro- tein structure determination process. This is particularly The Chemical Shift Index or CSI 3.0 (http://csi3. true for nuclear magnetic resonance (NMR)-based protein wishartlab.com) is a web server designed to accu- structure determination where secondary structure is used rately identify the location of secondary and super- to help in structure generation and refinement (2,3). In pro- secondary structures in protein chains using only tein NMR, secondary structures are traditionally identi- nuclear magnetic resonance (NMR) backbone chem- fied and assigned using NOE-based (Nuclear Overhauser ical shifts and their corresponding protein sequence Effect) methods. By manually analyzing the positions and data. Unlike earlier versions of CSI, which only iden- patterns of weak, medium or strong NOEs it is possible to tified three types of secondary structure (helix, ␤- identify helices, ␤-turns and ␤-stands with reasonably good strand and coil), CSI 3.0 now identifies total of 11 accuracy. Even today NOE pattern measurements continue types of secondary and super-secondary structures, to be the most commonly used method for identifying sec- ondary structures in peptides and proteins (2).
    [Show full text]
  • BMC Structural Biology Biomed Central
    BMC Structural Biology BioMed Central Research article Open Access Relationship between chemical shift value and accessible surface area for all amino acid atoms Wim F Vranken*1 and Wolfgang Rieping2 Address: 1Protein Databank Europe, European Bioinformatics Institute, Cambridge, UK and 2Department of Biochemistry, University of Cambridge, Cambridge, UK Email: Wim F Vranken* - [email protected]; Wolfgang Rieping - [email protected] * Corresponding author Published: 2 April 2009 Received: 12 November 2008 Accepted: 2 April 2009 BMC Structural Biology 2009, 9:20 doi:10.1186/1472-6807-9-20 This article is available from: http://www.biomedcentral.com/1472-6807/9/20 © 2009 Vranken and Rieping; 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: Chemical shifts obtained from NMR experiments are an important tool in determining secondary, even tertiary, protein structure. The main repository for chemical shift data is the BioMagResBank, which provides NMR-STAR files with this type of information. However, it is not trivial to link this information to available coordinate data from the PDB for non-backbone atoms due to atom and chain naming differences, as well as sequence numbering changes. Results: We here describe the analysis of a consistent set of chemical shift and coordinate data, in which we focus on the relationship between the per-atom solvent accessible surface area (ASA) in the reported coordinates and their reported chemical shift value.
    [Show full text]
  • Noor E Hafsa
    Protein Structure Characterization From NMR Chemical Shifts by Noor E Hafsa A thesis submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy Department of Computing Science University of Alberta © Noor E Hafsa, 2016 Abstract In order to understand the complex biological functions of proteins, highly detailed, atomic resolution protein structures are needed. Experimental methods such as X-ray crystallography and NMR spectroscopy provide standard platforms for determining the atomic-resolution structures of proteins. However, a continuing bottleneck in conventional NOE-based NMR structure determination lies in the difficulty of measuring NOEs for medium-to-large proteins and the resulting time-costs and the corresponding reduction in structure accuracy and precision. This has led to an increased interest in using other easily identifiable NMR parameters, such as chemical shifts, to facilitate protein structure determination by NMR. Chemical shifts, often considered as mileposts of NMR, have long been used to decipher the structures of small molecules. However, chemical shifts are much less frequently utilized for structural interpretation of larger macromolecules such as peptides and proteins. Most existing macromolecular methods use chemical shifts and various heuristic, rule-based algorithms to identify and determine a small number of structural parameters (such as secondary structure). Other methods, such as CS-Rosetta and CS23D, which attempt to determine 3D structures from chemical shifts alone, are only modestly successful (~50% success). So while good progress has been made, I believe that there is still substantial room for improvement and that the “Shift-to- Structure” problem has not yet been fully solved. My PhD project involves investigating innovative computational and machine- learning approaches to develop chemical-shift based prediction models to determine protein structures with high efficiency and high accuracy (>90%).
    [Show full text]
  • An Efficient Backbone Assignment Tool for NMR Studies of Proteins
    Reduced Dimensionality (4,3)D-hnCOCANH Experiment: An Efficient Backbone Assignment tool for NMR studies of Proteins Dinesh Kumar* Centre of Biomedical Magnetic Resonance, SGPGIMS Campus, Raibareli Road, Lucknow-226014, India; *Email: [email protected] KEYWORDS: Multidimensional NMR; Backbone Assignment; hncoCANH, hnCOcaNH; Sequential Correlations. ABBREVIATIONS: NMR, Nuclear Magnetic Resonance; HSQC, Heteronuclear Single Quantum Correlation; CARA, Computer Aided Resonance Assignment; RD: Reduced dimensionality. 1 Abstract: Sequence specific resonance assignment and secondary structure determination of proteins form the basis for variety of structural and functional proteomics studies by NMR. In this context, an efficient standalone method for rapid assignment of backbone (1H, 15N, 13C and 13C') resonances and secondary structure determination of proteins has been presented here. Compared to currently available strategies used for the purpose, the method employs only a single reduced dimensionality (RD) experiment -(4,3)D-hnCOCANH and exploits the linear combinations of backbone (13C and 13C') chemical shifts to achieve a dispersion relatively better compared to those of individual chemical shifts (see the text) for efficient and rapid data analysis. Further, the experiment leads to the spectrum with direct distinction of self (intra-residue) and sequential (inter-residue) carbon correlation peaks; these appear opposite in signs and therefore can easily be discriminated without using an additional complementary experiment. On top of
    [Show full text]
  • Bioinformatics Methods for NMR Chemical Shift Data
    Bioinformatics Methods for NMR Chemical Shift Data Dissertation an der Fakult¨at fur¨ Mathematik, Informatik und Statistik der Ludwig-Maximilians-Universit¨at Munchen¨ Simon W. Ginzinger vorgelegt am 24.10.2007 Erster Gutachter: Prof. Dr. Volker Heun, Ludwig-Maximilians-Universit¨at Munchen¨ Zweiter Gutachter: Prof. Dr. Robert Konrat, Universit¨at Wien Rigorosum: 8. Februar 2008 Kurzfassung Die nukleare Magnetresonanz-Spektroskopie (NMR) ist eine der wichtigsten Meth- oden, um die drei-dimensionale Struktur von Biomolekulen¨ zu bestimmen. Trotz großer Fortschritte in der Methodik der NMR ist die Aufl¨osung einer Protein- struktur immer noch eine komplizierte und zeitraubende Aufgabe. Das Ziel dieser Doktorarbeit ist es, Bioinformatik-Methoden zu entwickeln, die den Prozess der Strukturaufkl¨arung durch NMR erheblich beschleunigen k¨onnen. Zu diesem Zweck konzentriert sich diese Arbeit auf bestimmte Messdaten aus der NMR, die so genannten chemischen Verschiebungen. Chemische Verschiebungen werden standardm¨aßig zu Beginn einer Struktur- aufl¨osung bestimmt. Wie alle Labordaten k¨onnen chemische Verschiebungen Fehler enthalten, die die Analyse erschweren, wenn nicht sogar unm¨oglich machen. Als erstes Resultat dieser Arbeit wird darum CheckShift pr¨asentiert, eine Meth- ode, die es erm¨oglich einen weit verbreiteten Fehler in chemischen Verschiebungs- daten automatisch zu korrigieren. Das Hauptziel dieser Doktorarbeit ist es jedoch, strukturelle Informationen aus chemischen Verschiebungen zu extrahieren. Als erster Schritt in diese Richtung wurde SimShift entwickelt. SimShift erm¨oglicht es zum ersten Mal, strukturelle Ahnlichkeiten¨ zwischen Proteinen basierend auf chemischen Verschiebungen zu identifizieren. Der Vergleich zu Methoden, die nur auf der Aminos¨aurensequenz basieren, zeigt die Uberlegenheit¨ des verschiebungsbasierten Ansatzes. Als eine naturliche¨ Erweiterung des paarweisen Vergleichs von Proteinen wird darauf fol- gend SimShiftDB vorgestellt.
    [Show full text]
  • Anna Katharina Dehof Novel Approaches for Bond Order Assignment and NMR Shift Prediction
    Anna Katharina Dehof Novel Approaches for Bond Order Assignment and NMR Shift Prediction Novel Approaches for Bond Order Assignment and NMR Shift Prediction Dissertation zur Erlangung des Grades eines Doktors der Naturwissenschaften (Dr. rer. nat.) der Naturwissenschaftlich–Technischen Fakult¨at I – Mathematik und Informatik – der Universit¨at des Saarlandes vorgelegt von Anna Katharina Dehof, M.Sc. Saarbr¨ucken 2011 c 2011 Anna Katharina Dehof Covergestaltung und Titelbild: c 2011 Anna Katharina Dehof, BALL developer team, http://www.ball-project.org Herstellung und Verlag: Books on Demand GmbH, Norderstedt ISBN 978-3-8482-0773-2 Bibliografische Informationen der Deutschen Nationalbibliothek Die Deutsche Nationalbibliothek verzeichnet diese Publikation in der Deutschen Nationalbibliographie; detaillierte bibliografische Daten sind im Internetuber ¨ http://dnb.d-ndb.de abrufbar. Tag des Kolloquiums 15.03.2012 Dekan Prof. Dr. Holger Hermanns Berichterstatter Prof. Dr. Hans-Peter Lenhof Prof. Dr. Sebastian B¨ocker Vorsitz Prof. Dr. Philipp Slusallek Akad. Mitarbeiter Dr. Oliver M¨uller Eidesstattliche Versicherung Hiermit versichere ich an Eides statt, dass ich die vorliegende Arbeit selbstst¨andig und ohne Benutzung anderer als der angegebenen Hilfsmittel angefertigt habe. Die aus anderen Quellen oder indirekt ubernommenen¨ Daten und Konzepte sind unter Angabe der Quelle gekennzeichnet. Die Arbeit wurde bisher weder im In- noch im Ausland in gleicher oder ¨ahnlicher Form in einem Verfahren zur Erlangung eines akademischen Grades vorgelegt. Saarbr¨ucken, den 17. November 2011 (Anna Katharina Dehof) Abstract Molecular modelling is one of the cornerstones of modern biological and pharmaceutical research. Accurate modelling approaches easily become computationally overwhelming and thus, different levels of approximations are typically employed. In this work, we develop such approximation approaches for problems arising in structural bioinformatics.
    [Show full text]
  • Equilibrium Simulations of Proteins Using Molecular Fragment Replacement and NMR Chemical Shifts
    Equilibrium simulations of proteins using molecular fragment replacement and NMR chemical shifts Wouter Boomsmaa,1,2, Pengfei Tianb,1, Jes Frellsenc, Jesper Ferkinghoff-Borgd, Thomas Hamelrycka, Kresten Lindorff-Larsena, and Michele Vendruscoloe,2 aDepartment of Biology, University of Copenhagen, 2200 Copenhagen, Denmark; bNiels Bohr Institute, University of Copenhagen, 2100 Copenhagen, Denmark; cDepartment of Engineering, University of Cambridge, Cambridge, CB2 1PZ, United Kingdom; dBiotech Research and Innovation Center, University of Copenhagen, Technical University of Denmark Campus, 2800 Kongens Lyngby, Denmark; and eDepartment of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom Edited by José N. Onuchic, Rice University, Houston, TX, and approved August 7, 2014 (received for review March 18, 2014) Methods of protein structure determination based on NMR chemical data provide information primarily about local structural prop- shifts are becoming increasingly common. The most widely used erties of a molecule. In such cases, without modifying the force approaches adopt the molecular fragment replacement strat- field, the data can be integrated directly in the sampling pro- egy, in which structural fragments are repeatedly reassembled cedure, which can dramatically increase the efficiency of the into different complete conformations in molecular simulations. simulations. For instance, local structural information can be Although these approaches are effective in generating individual integrated into the simulations
    [Show full text]
  • An Overview of Tools for the Validation of Protein NMR Structures
    J Biomol NMR (2014) 58:259–285 DOI 10.1007/s10858-013-9750-x ARTICLE An overview of tools for the validation of protein NMR structures Geerten W. Vuister • Rasmus H. Fogh • Pieter M. S. Hendrickx • Jurgen F. Doreleijers • Aleksandras Gutmanas Received: 27 March 2013 / Accepted: 4 June 2013 / Published online: 23 July 2013 Ó Springer Science+Business Media Dordrecht 2013 Abstract Biomolecular structures at atomic resolution examples from the PDB: a small, globular monomeric present a valuable resource for the understanding of biol- protein (Staphylococcal nuclease from S. aureus, PDB ogy. NMR spectroscopy accounts for 11 % of all structures entry 2kq3) and a small, symmetric homodimeric protein (a in the PDB repository. In response to serious problems with region of human myosin-X, PDB entry 2lw9). the accuracy of some of the NMR-derived structures and in order to facilitate proper analysis of the experimental Keywords NMR Á Structure Á Validation Á Restraints Á models, a number of program suites are available. We Quality Á Program Á Review Á Chemical shifts discuss nine of these tools in this review: PROCHECK- NMR, PSVS, GLM-RMSD, CING, Molprobity, Vivaldi, ResProx, NMR constraints analyzer and QMEAN. We Introduction evaluate these programs for their ability to assess the structural quality, restraints and their violations, chemical Biomolecular structures at atomic resolution are crucial for shifts, peaks and the handling of multi-model NMR interpreting cellular processes in a molecular context. In ensembles. We document both the input required by the addition, they serve important roles in drug discovery and programs and output they generate.
    [Show full text]
  • Model-Based Assignment and Inference of Protein Backbone Nuclear Magnetic Resonances
    Model-based assignment and inference of protein backbone nuclear magnetic resonances. Olga Vitek∗ Jan Vitek† Bruce Craig∗ Chris Bailey-Kellogg†‡ Abstract Nuclear Magnetic Resonance (NMR) spectroscopy is a key experimental technique used to study protein structure, dynamics, and interactions. NMR methods face the bottleneck of spectral analysis, in particular determining the resonance assignment, which helps define the mapping between atoms in the protein and peaks in the spectra. A substantial amount of noise in spectral data, along with ambiguities in interpretation, make this analysis a daunting task, and there exists no generally accepted measure of uncertainty associated with the resulting solutions. This paper develops a model- based inference approach that addresses the problem of characterizing uncertainty in backbone resonance assignment. We argue that NMR spectra are subject to random variation, and ignoring this stochasticity can lead to false optimism and erroneous conclusions. We propose a Bayesian statistical model that accounts for various sources of uncertainty and provides an automatable framework for inference. While assignment has previously been viewed as a deterministic optimization problem, we demonstrate the importance of considering all solutions consistent with the data, and develop an algorithm to search this space within our statistical framework. Our approach is able to characterize the uncertainty associated with backbone resonance assignment in several ways: 1) it quantifies of uncertainty in the individually assigned resonances in terms of their posterior standard deviations; 2) it assesses the information content in the data with a posterior distribution of plausible assignments; and 3) it provides a measure of the overall plausibility of assignments. We demonstrate the value of our approach in a study of experimental data from two proteins, Human Ubiquitin and Cold-shock protein A from E.
    [Show full text]