Research Papers the SUPERFAMILY Database in Structural Genomics

Total Page:16

File Type:pdf, Size:1020Kb

Research Papers the SUPERFAMILY Database in Structural Genomics research papers Acta Crystallographica Section D Biological The SUPERFAMILY database in structural genomics Crystallography ISSN 0907-4449 Julian Gough The SUPERFAMILY hidden Markov model library repre- Received 31 January 2002 senting all proteins of known structure predicts the domain Accepted 21 August 2002 architecture of protein sequences and classi®es them at the MRC Laboratory of Molecular Biology, Hills Road, cambridge CB2 2QH, England SCOP superfamily level. This analysis has been carried out on all completely sequenced genomes. The ways in which the database can be useful to crystallographers is discussed, in Correspondence e-mail: [email protected] particular with a view to high-throughput structure determi- nation. The application of the SUPERFAMILY database to different target-selection strategies is suggested: novel folds, novel domain combinations and targeted attacks on genomes. Use of the database for more general inquiry in the context of structural studies is also explained. The database provides evolutionary relationships between target proteins and other proteins of known structure through the SCOP database, genome assignments and multiple sequence alignments. 1. Introduction As a result of the increase in the rate of experimental deter- mination of DNA sequences and the consequent success of the genome-sequencing projects, there are now (publicly avail- able) over 60 completely sequenced genomes spanning all kingdoms of life. More recently, there have arisen structural genomics projects (Smith, 2000), which are expected to begin reaching their production phases during 2002±2003. These projects will accelerate the production of new protein struc- tures and hence signi®cantly increase the available structural data. Because of the dif®culty and cost of solving three- dimensional structures, it is not possible to attempt to solve the structure of every protein. A more targeted approach than that used by the sequencing projects is needed to maximize the return on the projects' efforts. Using experimental and computational tools, targets can be selected which are expected to be in some way novel. Some of the ways in which the SUPERFAMILY (Gough & Chothia, 2002) database can contribute to targeting strategies are discussed here. Most targeting strategies aim to achieve as complete as possible coverage of something, for example a proteome or a functional pathway. There are a limited number of common structural folds (Chothia, 1992) and some targeting strategies of structural genomics projects will signi®cantly increase the proportion of this limited number for which we have a struc- tural representative. The improved completeness of fold coverage will change the current view of protein-structure space. 2. SUPERFAMILY # 2002 International Union of Crystallography The SUPERFAMILY database (Gough & Chothia, 2002; Printed in Denmark ± all rights reserved Gough et al., 2001) is a library of hidden Markov models Acta Cryst. (2002). D58, 1897±1900 Gough SUPERFAMILY database 1897 research papers (HMMs; Eddy, 1996; Hughey & Krogh, 1996; Krogh et al., consist of a single domain, whereas medium-sized proteins 1994) of domains of known structure created using the SAM may consist of one or more domains. Large proteins consist of (Karplus et al., 1998) software package. Services and data are multiple domains. A domain is de®ned as the minimum available at http://supfam.org. evolutionary unit, so a protein will only have parts classi®ed into separate domains if those parts are observed indepen- 2.1. What it does dently in nature either on their own or in combination with The purpose of SUPERFAMILY is to detect and classify in other domains. protein sequences evolutionary domains for which there is a Structural, functional and sequence information is used to known structural representative. Given a protein about which group together in the hierarchy at the superfamily level nothing is known other than the amino-acid sequence, the domains for which there is evidence for a common evolu- object is to assign known structural domains or more speci®- tionary ancestor. Domains belonging to the same superfamily cally domains at the SCOP (Murzin et al., 1995) superfamily have very similar structure and hence usually the same or a level. related function. 2.1.1. SCOP. The SCOP database classi®es all proteins in the PDB (Berman et al., 2000) into domains which are hier- 2.2. How it does it archically organized at levels of similarity. Small proteins SUPERFAMILY uses a library of HMMs representing all superfamilies in SCOP. HMMs are sequence pro®les very similar to PSI-BLAST (Altschul et al., 1997) pro®les. Pro®les are built from multiple sequence alignments and represent a group of sequences (in this case a superfamily) rather than a single sequence. The pro®le should embody the features that characterize the superfamily which it is supposed to represent. By comparing a sequence to a pro®le, far more distant rela- tionships can be detected than by comparing two sequences, which is what pairwise methods such as BLAST (Altschul et al., 1990) and FASTA (Pearson & Lipman, 1988) do. 2.3. Why use it? To compare the ability of different methods to detect and classify domains in SCOP (Murzin et al., 1995), a test was carried out, the results of which are shown in Fig. 1. The test Figure 1 The number of true and false positives found by four different methods comprises of an all-against-all search of sequences of struc- when scoring all sequences of structures in SCOP ®ltered to 40% tures in SCOP ®ltered to 40% sequence identity (Brenner et sequence identity. Figure 3 Figure 2 The SUPERFAMILY coverage of 58 complete genomes. The coverage of The number of domains assigned by SUPERFAMILY to the Escherichia genes is the percentage of genes with at least one domain assignment; the coli genome at different error rates (using different E-value thresholds). coverage of sequence is the percentage of amino acids in all genes In this case, the E values are calculated such that at 0.01 the expected covered by domain assignments. The two-letter codes for the genomes error rate is 1%. The curve shows that in the critical region around the shown in this graph are those de®ned by the SUPERFAMILY database. 1% error rate the number of domains assigned is not very sensitive to the The full names are available at http://supfam.org/SUPERFAMILY/ threshold. cgi-bin/genome_names.cgi. 1898 Gough SUPERFAMILY database Acta Cryst. (2002). D58, 1897±1900 research papers al., 2000). The SCOP classi®cation is then used to decide 3.1.2. Domain combinations. The HMM library has been whether the relationships detected are true or false. used to assign structural domains to sequences in all of the It is clear from Fig. 1 that pro®le methods (SAM-T99, PSI- completely sequenced genomes (see Fig. 3). Assignments BLAST) perform far better in this test of remote homology currently cover approximately 40% of the amino acids in detection than pairwise methods (WU-BLAST). SUPER- eukaryote genomes (50% of the sequences) and 45% of the FAMILY is based upon SAM-T99 (Karplus et al., 1998), but prokaryote genomes (55% of the sequences). Genome improves upon and adds to it, speci®cally aiming at SCOP sequences may have had some of their domains assigned but superfamily classi®cation. For more information, please refer not all, hence the discrepancy between the coverage of amino to Gough et al. (2001). acids and number of sequences. 2.3.1. Accuracy and reliability.TheE-value scores In the PDB there are structures of proteins with different provided by the SAM software which are used for the combinations of domains. By examining the assignments of assignments provide a theoretical value for the expected error domains to genome sequences, it is possible to observe protein rate. Large-scale tests show that these theoretical expectations sequences with combinations of domains for which a structure are very close to the observed error rates (Gough et al., 2001). has not yet been solved (Apic et al., 2001a,b). These novel Close examination of hundreds of assignments where the combinations provide targets for structural genomics projects results are checked (Bujnicki et al., 2001) indicate a 1% error and are available at http://www.mrc-lmb.cam.ac.uk/genomes/ rate. DomCombs. Although the novel aspect of these targets is not Examination of the assignments made by SUPERFAMILY the fold but the combination of folds, these targets offer at different error rates (E-value threshold) in Fig. 2 shows that certainty in this aspect. It is not possible (using any method) to in the critical region the number of assignments are not highly ®nd targets which are certain to have a novel fold. sensitive to the threshold chosen. 3.1.3. Targeted attacks. The aim of some projects is to achieve maximum structural coverage of a particular genome, e.g. Mycobacterium tuberculosis (http://www.doe-mbi.ucla.edu/ TB/pubs.php). As explained above, approximately half of all 3. Applications to structural genomics genome sequences are assigned to a superfamily for which The HMM library was designed for genome analysis leading to there is a known structural representative. The genome evolutionary studies (Apic et al., 2001a; Teichmann et al., analysis provides a good starting point for target selection on 2001), but also has applications in structural genomics. any project targeted at a particular genome. 3.1. Target selection 3.2. General inquiry The aim of structural genomics projects is to solve new As well as target selection, the model library can be used structures which give us a more complete view of the world of more generally to obtain information about a protein which protein structure. There are several different views of the could be relevant to structural studies. structural world, which lead to alternative approaches to the 3.2.1.
Recommended publications
  • Applied Category Theory for Genomics – an Initiative
    Applied Category Theory for Genomics { An Initiative Yanying Wu1,2 1Centre for Neural Circuits and Behaviour, University of Oxford, UK 2Department of Physiology, Anatomy and Genetics, University of Oxford, UK 06 Sept, 2020 Abstract The ultimate secret of all lives on earth is hidden in their genomes { a totality of DNA sequences. We currently know the whole genome sequence of many organisms, while our understanding of the genome architecture on a systematic level remains rudimentary. Applied category theory opens a promising way to integrate the humongous amount of heterogeneous informations in genomics, to advance our knowledge regarding genome organization, and to provide us with a deep and holistic view of our own genomes. In this work we explain why applied category theory carries such a hope, and we move on to show how it could actually do so, albeit in baby steps. The manuscript intends to be readable to both mathematicians and biologists, therefore no prior knowledge is required from either side. arXiv:2009.02822v1 [q-bio.GN] 6 Sep 2020 1 Introduction DNA, the genetic material of all living beings on this planet, holds the secret of life. The complete set of DNA sequences in an organism constitutes its genome { the blueprint and instruction manual of that organism, be it a human or fly [1]. Therefore, genomics, which studies the contents and meaning of genomes, has been standing in the central stage of scientific research since its birth. The twentieth century witnessed three milestones of genomics research [1]. It began with the discovery of Mendel's laws of inheritance [2], sparked a climax in the middle with the reveal of DNA double helix structure [3], and ended with the accomplishment of a first draft of complete human genome sequences [4].
    [Show full text]
  • Functional Effects Detailed Research Plan
    GeCIP Detailed Research Plan Form Background The Genomics England Clinical Interpretation Partnership (GeCIP) brings together researchers, clinicians and trainees from both academia and the NHS to analyse, refine and make new discoveries from the data from the 100,000 Genomes Project. The aims of the partnerships are: 1. To optimise: • clinical data and sample collection • clinical reporting • data validation and interpretation. 2. To improve understanding of the implications of genomic findings and improve the accuracy and reliability of information fed back to patients. To add to knowledge of the genetic basis of disease. 3. To provide a sustainable thriving training environment. The initial wave of GeCIP domains was announced in June 2015 following a first round of applications in January 2015. On the 18th June 2015 we invited the inaugurated GeCIP domains to develop more detailed research plans working closely with Genomics England. These will be used to ensure that the plans are complimentary and add real value across the GeCIP portfolio and address the aims and objectives of the 100,000 Genomes Project. They will be shared with the MRC, Wellcome Trust, NIHR and Cancer Research UK as existing members of the GeCIP Board to give advance warning and manage funding requests to maximise the funds available to each domain. However, formal applications will then be required to be submitted to individual funders. They will allow Genomics England to plan shared core analyses and the required research and computing infrastructure to support the proposed research. They will also form the basis of assessment by the Project’s Access Review Committee, to permit access to data.
    [Show full text]
  • Hmms Representing All Proteins of Known Structure. SCOP Sequence Searches, Alignments and Genome Assignments Julian Gough* and Cyrus Chothia
    268–272 Nucleic Acids Research, 2002, Vol. 30, No. 1 © 2002 Oxford University Press SUPERFAMILY: HMMs representing all proteins of known structure. SCOP sequence searches, alignments and genome assignments Julian Gough* and Cyrus Chothia MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 2QH, UK Received September 28, 2001; Revised and Accepted October 30, 2001 ABSTRACT (currently 59) covering approximately half of the The SUPERFAMILY database contains a library of soluble protein domains. The assignments, super- family breakdown and statistics on them are available hidden Markov models representing all proteins of from the server. The database is currently used by known structure. The database is based on the SCOP this group and others for genome annotation, structural ‘superfamily’ level of protein domain classification genomics, gene prediction and domain-based genomic which groups together the most distantly related studies. proteins which have a common evolutionary ancestor. There is a public server at http://supfam.org which provides three services: sequence searching, INTRODUCTION multiple alignments to sequences of known structure, The SUPERFAMILY database is based on the SCOP (1) and structural assignments to all complete genomes. classification of protein domains. SCOP is a structural domain- Given an amino acid or nucleotide query sequence based heirarchical classification with several levels including the server will return the domain architecture and the ‘superfamily’ level. Proteins grouped together at the super- SCOP classification. The server produces alignments family level are defined as having structural, functional and sequence evidence for a common evolutionary ancestor. It is at of the query sequences with sequences of known this level, as the name suggests, that SUPERFAMILY operates structure, and includes multiple alignments of because it is the level with the most distantly related protein genome and PDB sequences.
    [Show full text]
  • AUTOPHY by Deepika Prasad a Thesis
    ANALYZING MARKER GENE DIVERSITY USING AN AUTOMATED PHYLOGENETIC TOOL: AUTOPHY by Deepika Prasad A thesis submitted to the Faculty of the University of Delaware in partial fulfillment of the requirements for the degree of Master of Science in Bioinformatics and Computational Biology Spring 2017 © 2017 Deepika Prasad All Rights Reserved ANALYZING MARKER GENE DIVERSITY USING AN AUTOMATED PHYLOGENETIC TOOL: AUTOPHY by Deepika Prasad Approved: __________________________________________________________ Shawn Polson, Ph.D. Professor in charge of thesis on behalf of the Advisory Committee Approved: __________________________________________________________ Kathleen F.McCoy, Ph.D. Chair of the Department of Computer and Information Sciences Approved: __________________________________________________________ Babatunde A. Ogunnaike, Ph.D. Dean of the College of Engineering Approved: __________________________________________________________ Ann L. Ardis, Ph.D. Senior Vice Provost for Graduate and Professional Education ACKNOWLEDGMENTS I want to express my sincerest gratitude to Professor Shawn Polson for his patience, continuous support, and enthusiasm. I could not have imagined having a better mentor to guide me through the process of my Master’s thesis. I would like to thank my committee members, Professor Eric Wommack and Honzhan Huang for their insightful comments and guidance in the thesis. I would also like to express my gratitude to Barbra Ferrell for asking important questions throughout my thesis research and helping me with the writing process. I would like to thank my friend Sagar Doshi, and lab mates Daniel Nasko, and Prasanna Joglekar, for helping me out with data, presentations, and giving important advice, whenever necessary. Last but not the least, I would like to thank my family, for being my pillars of strength.
    [Show full text]
  • Protein Domain Coocurences Reveal Functional Changes of Regulatory Mechanisms During Evolution
    Protein Domain Coocurences Reveal Functional Changes of Regulatory Mechanisms During Evolution A.A.Parikesit*, S.J. Prohaska, P.F. Stadler Chair of Bioinformatics, University of Leipzig Correlation between transcription factors and chromatin related proteins Introduction Figure 3 : Correlation of the number of transcription factors and chromatin domains. The emergence of higher organisms was facilitated by a Blue box. Species with few chromatin-related domains but many dramatic increases in the complexity of gene regulatory transcription factors: human, kangaroo rat, mouse, opossum, fugu. Green box: species with many chromatin-related domains but few mechanisms. This is achieved not only by addition of novel transcription factors: seq squirt, medaka, dolphin, yeast but also by expansion of existing mechanisms. Such an There are no organisms (known) that have a lot of both. expansion is usually characterized by the proliferation of functionally paralogous proteins and the appearance of novel combinations of functional domains. Large scale phylogenetic analysis can shed light on the relative amounts of functional domains and their combinations and interactions involved in Results show massive problems with data quality: closely related species certain regulatory networks. (e.g, dolphin and human) show dramatically different distributions of transcription factors and chromatin domains. This is not reasonable within Methods mammals and contradicts biological knowledge. We performed comparative and functional analysis of three SUPERFAMILY thus cannot be used for large-scale quantitative regulatory mechanisms: (1) transcriptional regulation by comparisons across species due to several sources of bias: transcription factors, (2) post-transcriptional regulation by - different completeness of protein annotation for different genomes miRNAs, and (3) chromatin regulation across all domains of - differences in transcript coverage life.
    [Show full text]
  • Renaissance SUPERFAMILY in the Decade Since Its Launch, a Repository of Information About Proteins in Genomes Has Developed Into a Primary Reference
    Renaissance SUPERFAMILY In the decade since its launch, a repository of information about proteins in genomes has developed into a primary reference. Dr Julian Gough, its creator, describes current work enhancing the scope, content and functionality of this key service Could you begin by outlining the reasons for How do HMMs feed into the service? creating SUPERFAMILY? HMMs are profi les which represent multiple SUPERFAMILY was originally created to better sequence alignments of homologous proteins understand molecular evolution, initially by in a rigorous statistical framework. They enabling comparison of the repertoire of proteins can be used to classify sequences based on and domains across the genomes of different homology and to create sequence alignments. species. The starting point was the Structural We use sequences of domains of known Classifi cation of Proteins (SCOP), and the most structure via iterative search procedures on basic purpose of SUPERFAMILY is to detect and large background sequence databases to classify these domains with known structural build alignments and, subsequently, models representatives in the protein sequences of representing domains of known structure at genomes. There are tens of thousands of known the superfamily level. These models are then protein structures, each the result of a costly searched against genome sequences to detect three-dimensional atomic resolution structure and classify the structural domains. determination by experiment, usually X-ray crystallography or nuclear magnetic resonance.
    [Show full text]
  • An Expanded Evaluation of Protein Function Prediction Methods Shows an Improvement in Accuracy
    The University of Southern Mississippi The Aquila Digital Community Faculty Publications 9-7-2016 An Expanded Evaluation of Protein Function Prediction Methods Shows an Improvement In Accuracy Yuxiang Jiang Indiana University TFollowal Ronnen this and Or onadditional works at: https://aquila.usm.edu/fac_pubs Buck P arInstitutet of the forGenomics Research Commons On Agin Wyatt T. Clark RecommendedYale University Citation AsmaJiang, YR.., OrBankapuron, T. R., Clark, W. T., Bankapur, A. R., D'Andrea, D., Lepore, R., Funk, C. S., Kahanda, I., Verspoor, K. M., Ben-Hur, A., Koo, D. E., Penfold-Brown, D., Shasha, D., Youngs, N., Bonneau, R., Lin, A., Sahraeian, S. Miami University M., Martelli, P. L., Profiti, G., Casadio, R., Cao, R., Zhong, Z., Cheng, J., Altenhoff, A., Skunca, N., Dessimoz, DanielC., Dogan, D'Andr T., Hakala,ea K., Kaewphan, S., Mehryar, F., Salakoski, T., Ginter, F., Fang, H., Smithers, B., Oates, M., UnivGough,ersity J., ofTör Romeönen, P., Koskinen, P., Holm, L., Chen, C., Hsu, W., Bryson, K., Cozzetto, D., Minneci, F., Jones, D. T., Chapan, S., BKC, D., Khan, I. K., Kihara, D., Ofer, D., Rappoport, N., Stern, A., Cibrian-Uhalte, E., Denny, P., Foulger, R. E., Hieta, R., Legge, D., Lovering, R. C., Magrane, M., Melidoni, A. N., Mutowo-Meullenet, P., SeePichler next, K., page Shypitsyna, for additional A., Li, B.,authors Zakeri, P., ElShal, S., Tranchevent, L., Das, S., Dawson, N. L., Lee, D., Lees, J. G., Stilltoe, I., Bhat, P., Nepusz, T., Romero, A. E., Sasidharan, R., Yang, H., Paccanaro, A., Gillis, J., Sedeño- Cortés, A. E., Pavlidis, P., Feng, S., Cejuela, J.
    [Show full text]
  • Simulation and Graph Mining Tools for Improving Gene Mapping Efficiency
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Helsingin yliopiston digitaalinen arkisto DEPARTMENT OF COMPUTER SCIENCE SERIES OF PUBLICATIONS A REPORT A-2011-3 Simulation and graph mining tools for improving gene mapping efficiency Petteri Hintsanen To be presented, with the permission of the Faculty of Science of the University of Helsinki, for public criticism in Auditorium XIV (Univer- sity Main Building, Unioninkatu 34) on 30 September 2011 at twelve o’clock noon. UNIVERSITY OF HELSINKI FINLAND Supervisor Hannu Toivonen, University of Helsinki, Finland Pre-examiners Tapio Salakoski, University of Turku, Finland Janne Nikkila,¨ University of Helsinki, Finland Opponent Michael Berthold, Universitat¨ Konstanz, Germany Custos Hannu Toivonen, University of Helsinki, Finland Contact information Department of Computer Science P.O. Box 68 (Gustaf Hallstr¨ omin¨ katu 2b) FI-00014 University of Helsinki Finland Email address: [email protected].fi URL: http://www.cs.Helsinki.fi/ Telephone: +358 9 1911, telefax: +358 9 191 51120 Copyright c 2011 Petteri Hintsanen ISSN 1238-8645 ISBN 978-952-10-7139-3 (paperback) ISBN 978-952-10-7140-9 (PDF) Computing Reviews (1998) Classification: G.2.2, G.3, H.2.5, H.2.8, J.3 Helsinki 2011 Helsinki University Print Simulation and graph mining tools for improving gene mapping efficiency Petteri Hintsanen Department of Computer Science P.O. Box 68, FI-00014 University of Helsinki, Finland petteri.hintsanen@iki.fi http://iki.fi/petterih PhD Thesis, Series of Publications A, Report A-2011-3 Helsinki, September 2011, 136 pages ISSN 1238-8645 ISBN 978-952-10-7139-3 (paperback) ISBN 978-952-10-7140-9 (PDF) Abstract Gene mapping is a systematic search for genes that affect observable character- istics of an organism.
    [Show full text]
  • Specialized Hidden Markov Model Databases for Microbial Genomics
    Comparative and Functional Genomics Comp Funct Genom 2003; 4: 250–254. Published online 1 April 2003 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cfg.280 Conference Review Specialized hidden Markov model databases for microbial genomics Martin Gollery* University of Nevada, Reno, 1664 N. Virginia Street, Reno, NV 89557-0014, USA *Correspondence to: Abstract Martin Gollery, University of Nevada, Reno, 1664 N. Virginia As hidden Markov models (HMMs) become increasingly more important in the Street, Reno, NV analysis of biological sequences, so too have databases of HMMs expanded in size, 89557-0014, USA. number and importance. While the standard paradigm a short while ago was the E-mail: [email protected] analysis of one or a few sequences at a time, it has now become standard procedure to submit an entire microbial genome. In the future, it will be common to submit large groups of completed genomes to run simultaneously against a dozen public databases and any number of internally developed targets. This paper looks at some of the readily available HMM (or HMM-like) algorithms and several publicly available Received: 27 January 2003 HMM databases, and outlines methods by which the reader may develop custom Revised: 5 February 2003 HMM targets. Copyright 2003 John Wiley & Sons, Ltd. Accepted: 6 February 2003 Keywords: HMM; Pfam; InterPro; SuperFamily; TLfam; COG; TIGRfams Introduction will be a true homologue. As a result, HMMs have become very popular in the field of bioinformatics Over the last few years, hidden Markov models and a number of HMM databases have been (HMMs) have become one of the pre-eminent developed.
    [Show full text]
  • Download Flyer
    COLD SPRING HARBOR ASIA ORGANIZERS(Speaker, Affiliation, Country/Region) Steven E. Brenner Frontiers in University of California, Berkeley, USA A Keith Dunker Indiana University School of Medicine, USA Computational Biology Julian Gough MRC Laboratory of Molecular Biology, UK Suzhou, China September 3-7, 2018 Luhua Lai & Bioinformatics Peking University, China Yunlong Liu Abstract deadline: July 13, 2018 Indiana University School of Medicine, USA MAJOR TOPICS Precision medicine, human genome variation, disease & diagnosis Molecular evolution Pathways, networks & developmental biology Molecular structure, with pioneering techniques Molecular machines, their functions & dynamics Intrinsically disordered proteins & their functions RNA function, regulation & splicing 3D genomics & regulatory inferences Single cell analysis KEYNOTE SPEAKERS (Speaker, Affiliation, Country/Region) Nancy Cox, Vanderbilt University, USA Yoshihide Hayashizaki, RIKEN Research Cluster for Innovation, JAPAN INVITED SPEAKERS (Speaker, Affiliation, Country/Region) Russ Altman, Stanford University, USA Manolis Kellis, MIT Computer Science and Broad Institute, USA Lukasz Kurgan, Virginia Commonwealth University, USA Peer Bork, European Molecular Biology Laboratory, GERMANY Luhua Lai, Peking University, CHINA Steven Brenner, University of California, Berkeley, USA Michal Linial, The Hebrew University of Jerusalem, ISRAEL Angela Brooks, University of California, Santa Cruz, USA Yunlong Liu, Indiana University School of Medicine, USA Luonan Chen, Shanghai Institutes for
    [Show full text]
  • Assignment of Homology to Genome
    doi:10.1006/jmbi.2001.5080availableonlineathttp://www.idealibrary.comon J. Mol. Biol. (2001) 313, 903±919 AssignmentofHomologytoGenomeSequences using a Library of Hidden Markov Models that Represent all Proteins of Known Structure JulianGough1*,KevinKarplus2,RichardHughey2andCyrusChothia1 1MRC, Laboratory of Molecular Of the sequence comparison methods, pro®le-based methods perform Biology, Hills Road, Cambridge with greater selectively than those that use pairwise comparisons. Of the CB2 2QH, UK pro®le methods, hidden Markov models (HMMs) are apparently the best. The ®rst part of this paper describes calculations that (i) improve 2Department of Computer the performance of HMMs and (ii) determine a good procedure for creat- Engineering, Jack Baskin School ing HMMs for sequences of proteins of known structure. For a family of of Engineering, University of related proteins, more homologues are detected using multiple models California, Santa Cruz built from diverse single seed sequences than from one model built from CA 95064, USA a good alignment of those sequences. A new procedure is described for detecting and correcting those errors that arise at the model-building stage of the procedure. These two improvements greatly increase selectiv- ity and coverage. The second part of the paper describes the construction of a library of HMMs, called SUPERFAMILY, that represent essentially all proteins of known structure. The sequences of the domains in proteins of known structure, that have identities less than 95 %, are used as seeds to build the models. Using the current data, this gives a library with 4894 models. The third part of the paper describes the use of the SUPERFAMILY model library to annotate the sequences of over 50 genomes.
    [Show full text]
  • Sequences and Topology: Disorder, Modularity, and Post/Pre Translation Modification
    COSTBI-1118; NO. OF PAGES 3 Available online at www.sciencedirect.com Sequences and topology: disorder, modularity, and post/pre translation modification Editorial overview Julian Gough and A Keith Dunker Current Opinion in Structural Biology 2013, 23:xx– yy This review comes from a themed issue on Se- quences and topology 0959-440X/$ – see front matter, # 2013 Elsevier Ltd. All rights reserved. http://dx.doi.org/10.1016/j.sbi.2013.05.001 Julian Gough In this collection of reviews we bring together discussions of intrinsically Department of Computer Science, University disordered proteins and their functional mechanisms operating at different of Bristol, Merchant Venturers Building, loci across the ordered-to-disordered spectrum. Our previous editorial [1] for Woodland Road, Bristol BS8 1UB, UK sequences and topology in 2011 contrasts and compares structured and e-mail: [email protected] disordered proteins but viewed in light of their evolution. Here we move forward, focusing on function across the spectrum of disorder, but still within Julian Gough began his research career in structural bioinformatics, but soon combined the context of evolution. Disordered regions can also participate as modules this with genomics and sequence analysis to in multi-domain proteins; the reviews herein cover both structured [2] and study evolution, also producing the disordered modularity [3]. A growing cohort of researchers suggest that we SUPERFAMILY database of structural need to improve the connections between structure, disorder and their domain assignments to genomes. Gough functions, and several reviews present computational and experimental continues to examine molecular evolution methodologies for doing so. A running theme is the importance to function across all species and applies bioinformatics at multiple levels; most recently to intrinsic of modifications, both post and pre-translation.
    [Show full text]