DOCSLIB.ORG

Development and Dissemination of Computational Methods for Genome-Scale Modeling

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Development and Dissemination of Computational Methods for Genome-Scale Modeling UC San Diego UC San Diego Electronic Theses and Dissertations Title Development and Dissemination of Computational Methods for Genome-scale Modeling Permalink https://escholarship.org/uc/item/1w98620m Author Ebrahim, Ali Publication Date 2016 Peer reviewed|Thesis/dissertation eScholarship.org Powered by the California Digital Library University of California UNIVERSITY OF CALIFORNIA, SAN DIEGO Development and Dissemination of Computational Methods for Genome-scale Modeling A Dissertation submitted in partial satisfaction of the requirements for the degree Doctor of Philosophy in Bioengineering by Ali Ebrahim Committee in charge: Professor Bernhard Palsson, Chair Professor Nuno Bandeira Professor Christian Metallo Professor Glenn Tesler Professor Kun Zhang 2016 Copyright Ali Ebrahim, 2016 All rights reserved. The Dissertation of Ali Ebrahim is approved, and it is ac- ceptable in quality and form for publication on microfilm and electronically: Chair University of California, San Diego 2016 iii DEDICATION My heart breaks for the refugees who have been driven from their homes in Syria. Though it must mean very little, I dedicate this thesis to them, and I will try to apply some of what I have learned to help give them some of the opportunities I have had. iv EPIGRAPH What I cannot create, I do not understand. |Richard P. Feynman I like to call it low-input, high-throughput, no-output biology. |Sydney Brenner v TABLE OF CONTENTS Signature Page................................... iii Dedication...................................... iv Epigraph......................................v Table of Contents.................................. vi List of Figures................................... xi List of Tables.................................... xiii Acknowledgements................................. xiv Vita......................................... xix Abstract of the Dissertation............................ xxi Chapter 1 COBRApy: COnstraints-Based Reconstruction and Analysis for Python................................1 1.1 Abstract............................1 1.1.1 Background......................1 1.1.2 Results.........................2 1.1.3 Conclusion.......................2 1.2 Background..........................2 1.3 Implementation........................4 1.4 Results and discussion....................4 1.4.1 Core classes: model, metabolite, reaction, & gene.6 1.4.2 Key capabilities....................8 1.4.3 Advanced capabilities.................9 1.5 Conclusions.......................... 10 1.6 Acknowledgements...................... 10 Chapter 2 Distribution and Portability of Validated Genome-Scale Constraint-Based Models...................... 12 2.1 Introduction and Problem Statement............ 12 2.2 A Potential Solution..................... 13 2.3 Acknowledgments....................... 15 Chapter 3 Model-driven elucidation of transcriptional regulatory network in bacteria............................... 16 3.1 Abstract............................ 16 3.2 Introduction.......................... 17 vi 3.3 Results............................. 19 3.3.1 Model-driven prediction of activation conditions for TFs........................... 19 3.3.2 Experimental confirmation of predicted conditions. 21 3.3.3 Confirmation of ChIP-exo binding sites with motif analysis and comparison to known sites....... 23 3.3.4 Reconstruction of NtrC and Nac regulons...... 25 3.3.5 Association of TF with σ-factors........... 26 3.3.6 Contrasting functions of NtrC and Nac regulons.. 28 3.3.7 Primary response of the NtrC regulon........ 29 3.3.8 Carbon flux rebalancing by the Nac regulon.... 32 3.3.9 Lower activation of Nac on cytidine......... 33 3.3.10 Coupling of nitrogen and carbon metabolism.... 35 3.4 Discussion........................... 37 3.4.1 Selection of optimal experimental conditions by model-driven experimental design.......... 37 3.4.2 Model-guided regulon elucidation and characterization 39 3.4.3 Application of ChIP-exo method in bacterial studies 40 3.4.4 TF binding sites in non-regulatory regions..... 41 3.4.5 Nitrogen starvation is coupled to other responses.. 42 3.5 Materials and Methods.................... 43 3.6 Author Contributions..................... 44 3.7 Acknowledgements...................... 44 Chapter 4 Multi-omic data integration enables discovery of hidden biological regularities.............................. 52 4.1 Abstract............................ 52 4.2 Results and Discusssion.................... 53 4.3 Acknowledgements...................... 60 Chapter 5 Construction of Genome-scale models of Metabolism and Gene Expression for the Escherichia coli family............ 64 5.1 Abstract............................ 64 5.2 Introduction.......................... 65 5.3 Results and Discussion.................... 67 5.3.1 The \miniME" Reformulation of ME models.... 67 5.3.2 Global Parameter Characterization......... 69 5.3.3 Extending from E. coli K12 MG1655 to 40 related strains......................... 70 5.3.4 Conclusions...................... 72 5.4 Materials and Methods.................... 72 5.4.1 Optimization Procedure............... 72 5.4.2 Multiple ME Model Construction.......... 74 5.5 Acknowledgements...................... 75 vii Appendix A Documentation for COBRApy................... 79 A.1 Getting Started........................ 79 A.1.1 Reactions....................... 80 A.1.2 Metabolites...................... 82 A.1.3 Genes......................... 83 A.2 Building a Model....................... 85 A.3 Reading and Writing Models................. 89 A.3.1 SBML......................... 89 A.3.2 JSON......................... 90 A.3.3 MATLAB....................... 91 A.3.4 Pickle......................... 91 A.4 Simulating with FBA..................... 92 A.4.1 Running FBA..................... 92 A.4.2 Changing the Objectives............... 94 A.4.3 Running FVA..................... 95 A.4.4 Running pFBA.................... 98 A.5 Simulating Deletions..................... 98 A.5.1 Single Deletions.................... 98 A.5.2 Double Deletions................... 100 A.6 Phenotype Phase Plane.................... 101 A.7 Mixed-Integer Linear Programming............. 105 A.7.1 Ice Cream....................... 105 A.7.2 Restaurant Order................... 106 A.7.3 Boolean Indicators.................. 107 A.8 Quadratic Programming................... 110 A.9 Loopless FBA......................... 113 A.10 Gapfillling........................... 117 A.10.1 GrowMatch...................... 117 A.10.2 SMILEY........................ 119 A.11 Solver Interface........................ 119 A.11.1 Attributes and functions............... 120 A.11.2 Example with FVA.................. 124 A.12 Using the COBRA toolbox with cobrapy.......... 125 Appendix B Distributing Validated Constraint-Based Models in SBML... 128 B.1 Systems Biology Ontology Terms and COBRA models... 128 B.2 Strict Models......................... 130 B.3 Defining Flux Bounds..................... 130 B.4 Reaction Reversibility..................... 132 B.5 Gene Reaction Rules..................... 133 B.6 Biomass Objective Reactions................. 134 B.7 Conventions for Identifiers.................. 135 viii Appendix C Model-driven elucidation of transcriptional regulatory network in bacteria - Supplementary Information............... 136 C.1 Advantages of ChIP-exo over previous methods...... 136 C.2 Locally acting TFs regulated by Nac............. 137 C.3 Supplementary Experimental Procedures.......... 138 C.3.1 Prediction of activation conditions for transcription factor.......................... 138 C.3.2 Bacterial strains, media, and growth conditions.. 139 C.3.3 ChIP-exo experiment................. 140 C.3.4 RNA-seq expression profiling............. 141 C.3.5 Peak calling for ChIP-exo dataset.......... 143 C.3.6 Classification of regulatory and non-regulatory bind- ing sites........................ 143 C.3.7 Motif search from ChIP-exo peaks.......... 144 C.3.8 Calculation of differentially expressed gene..... 144 C.3.9 COG functional enrichment............. 145 C.3.10 Measuring growth rate and glucose uptake rates on different nitrogen sources............... 145 C.3.11 FBA analysis with glucose uptake rates....... 146 C.3.12 Comparison of in silico growth with experimental evidence........................ 146 C.3.13 Conservation analysis of nitrogen-related genes... 147 C.4 Supplementary Figures.................... 148 Appendix D Multi-omic data integration enables discovery of hidden biological regularities - Supplementary Information............. 163 D.1 Methods............................ 163 D.1.1 mRNA seq methods.................. 163 D.1.2 Ribosome profiling.................. 164 D.1.3 Genome-wide Secondary Structure Annotations.. 165 D.1.4 Tertiary Structure/ Protein Domain annotations.. 167 D.1.5 Identification of Shine-Dalgarno-like codons.... 167 D.1.6 Ribosome density/pause site accounting (Figures D.2- D.3).......................... 168 D.1.7 Computational Method for Predicting keff parameters 169 D.1.8 Predictions of mRNA expression under identical con- ditions to proteomic data............... 172 D.1.9 Predicting Differential Gene Expression with iOL1650-ME...................... 173 D.1.10 Sampling of M-model flux states in iJO1366.... 174 D.2 Structure of ME models................... 175 D.3 ME model coupling parameters............... 178 D.4 Simulation of Batch Growth with ME............ 180 ix D.5 Simulating ME with estimated parameters......... 180 D.6 Ribosome profiling pause site analysis............ 182 D.7 Supplementary Figures.................... 184
Load more

Recommended publications
  • Meeting Review: Bioinformatics and Medicine – from Molecules To
    Comparative and Functional Genomics Comp Funct Genom 2002; 3: 270–276. Published online 9 May 2002 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cfg.178 Feature Meeting Review: Bioinformatics And Medicine – From molecules to humans, virtual and real Hinxton Hall Conference Centre, Genome Campus, Hinxton, Cambridge, UK – April 5th–7th Roslin Russell* MRC UK HGMP Resource Centre, Genome Campus, Hinxton, Cambridge CB10 1SB, UK *Correspondence to: Abstract MRC UK HGMP Resource Centre, Genome Campus, The Industrialization Workshop Series aims to promote and discuss integration, automa- Hinxton, Cambridge CB10 1SB, tion, simulation, quality, availability and standards in the high-throughput life sciences. UK. The main issues addressed being the transformation of bioinformatics and bioinformatics- based drug design into a robust discipline in industry, the government, research institutes and academia. The latest workshop emphasized the influence of the post-genomic era on medicine and healthcare with reference to advanced biological systems modeling and simulation, protein structure research, protein-protein interactions, metabolism and physiology. Speakers included Michael Ashburner, Kenneth Buetow, Francois Cambien, Cyrus Chothia, Jean Garnier, Francois Iris, Matthias Mann, Maya Natarajan, Peter Murray-Rust, Richard Mushlin, Barry Robson, David Rubin, Kosta Steliou, John Todd, Janet Thornton, Pim van der Eijk, Michael Vieth and Richard Ward. Copyright # 2002 John Wiley & Sons, Ltd. Received: 22 April 2002 Keywords: bioinformatics;
    [Show full text]
  • Assessing Sequence Comparison Methods with Reliable Structurally Identified Distant Evolutionary Relationships
    Proc. Natl. Acad. Sci. USA Vol. 95, pp. 6073–6078, May 1998 Biochemistry Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships STEVEN E. BRENNER*†‡,CYRUS CHOTHIA*, AND TIM J. P. HUBBARD§ *MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 2QH, United Kingdom; and §Sanger Centre, Wellcome Trust Genome Campus, Hinxton, Cambs CB10 1SA, United Kingdom Communicated by David R. Davies, National Institute of Diabetes, Bethesda, MD, March 16, 1998 (received for review November 12, 1997) ABSTRACT Pairwise sequence comparison methods have Sequence comparison methodologies have evolved rapidly, been assessed using proteins whose relationships are known so no previously published tests has evaluated modern versions reliably from their structures and functions, as described in of programs commonly used. For example, parameters in the SCOP database [Murzin, A. G., Brenner, S. E., Hubbard, T. BLAST (1) have changed, and WU-BLAST2 (2)—which produces & Chothia C. (1995) J. Mol. Biol. 247, 536–540]. The evalua- gapped alignments—has become available. The latest version tion tested the programs BLAST [Altschul, S. F., Gish, W., of FASTA (3) previously tested was 1.6, but the current release Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. (version 3.0) provides fundamentally different results in the 215, 403–410], WU-BLAST2 [Altschul, S. F. & Gish, W. (1996) form of statistical scoring. Methods Enzymol. 266, 460–480], FASTA [Pearson, W. R. & The previous reports also have left gaps in our knowledge. Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444–2448], For example, there has been no published assessment of and SSEARCH [Smith, T.
    [Show full text]
  • A Novel and Effective Scoring Scheme for Structure Classification And
    A novel and effective scoring scheme for structure classification and pairwise similarity measurement Rezaul Karim∗, Md. Momin Al Aziz†,Swakkhar Shatabda‡ and M. Sohel Rahman∗ ∗Department of Computer Science and Engineering, Bangladesh University of Engineering and Technology †Department of Computer Science, University of Manitoba, Canada ‡Department of Computer Science and Engineering, United International University Abstract—Protein tertiary structure defines its functions, clas- However, a major drawback of these approaches with re- sification and binding sites. Similar structural characteristics gards to the global structure similarity comparison is that these between two proteins often lead to the similar characteristics are much sensitive to local structure dissimilarity. Furthermore, thereof. Determining structural similarity accurately in real these approaches need to find an optimal alignment before time is a crucial research issue. In this paper, we present computing the score despite that often there are no remarkable a novel and effective scoring scheme that is dependent on alignment possible. Also finding proper alignment between two novel features extracted from protein alpha carbon distance matrices. Our scoring scheme is inspired from pattern recog- protein structures is computationally expensive if the size of nition and computer vision. Our method is significantly better the protein is large [9]. In spite of these inherent limitations than the current state of the art methods in terms of fam- and drawbacks, there exist a number of such methods in the ily match of pairs of protein structures and other statistical literature and the most notable ones among these are DALI measurements. The effectiveness of our method is tested on [10], [11], CE [12], TM Align [13], [14] and SP Align [15].
    [Show full text]
  • Relative Orientation of Close-Packed ,8-Pleated Sheets in Proteins
    Proc. Nati Acad. Sci. USA Vol. 78, No. 7, pp. 4146-4150, July 1981 Biochemistry Relative orientation of close-packed ,8-pleated sheets in proteins (twisted ,3sheets/protein secondary and tertiary structure) CYRUS CHOTHIA*t AND JOEL JANINt *William Ramsay, Ralph Foster and Christopher Ingold Laboratories, University College London, Gower Street, London WC1E 6BT, England; tMedical Research Council Laboratory of Molecular Biology, Hills Road, Cambridge CB2 2QH, England; and *Unit6 de Biochimie Cellulaire, Department de Biochimie et Gen~tique Moleculaire, Institut Pasteur, 28 rue du Docteur Roux, 75724 Paris Cedex 15, France Communicated by Max F. Perutz, March 30, 1981 ABSTRACT When (3-pleated sheets pack face to face in pro- curve from left to right, while those that point down curve from teins, the angle between the strand directions ofthe two (3-sheets right to left. is observed to be near -30°. We propose a simple model for f3- Fig. 2 b and c illustrates some of the effects of twisting two sheet-to-(-sheet packing in concanavalin A, plastocyanin, y-crys- ,(3pleated sheets that are packed together. We show, on the left tallin, superoxide dismutase, prealbumin, and the immunoglobin ofthe figures, two hypothetical untwisted (3-sheets packed with fragment VREI. This model shows how the observed relative ori- their rows of side chains aligned: those in the top (-sheet are entation of two packed (3-sheets is a consequence of (i) the rows parallel to those in the bottom. On the right of the figures we of side chains at the interface being approximately aligned and show what happens when these packed (3-sheets are twisted (ii) the (3-sheet having a right-handed twist.
    [Show full text]
  • Transformer Neural Networks for Protein Prediction Tasks
    bioRxiv preprint doi: https://doi.org/10.1101/2020.06.15.153643; this version posted June 16, 2020. 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 4.0 International license. TRANSFORMING THE LANGUAGE OF LIFE:TRANSFORMER NEURAL NETWORKS FOR PROTEIN PREDICTION TASKS Ananthan Nambiar ∗ Maeve Heflin∗ Department of Bioengineering Department of Computer Science Carl R. Woese Inst. for Genomic Biol. Carl R. Woese Inst. for Genomic Biol. University of Illinois at Urbana-Champaign University of Illinois at Urbana-Champaign Urbana, IL 61801 Urbana, IL 61801 [email protected] Simon Liu∗ Sergei Maslov Department of Computer Science Department Bioengineering Carl R. Woese Inst. for Genomic Biol. Department of Physics University of Illinois at Urbana-Champaign Carl R. Woese Inst. for Genomic Biol. Urbana, IL 61801 University of Illinois at Urbana-Champaign Urbana, IL 61801 Mark Hopkinsy Anna Ritzy Department of Computer Science Department of Biology Reed College Reed College Portland, OR 97202 Portland, OR 97202 June 16, 2020 ABSTRACT The scientific community is rapidly generating protein sequence information, but only a fraction of these proteins can be experimentally characterized. While promising deep learning approaches for protein prediction tasks have emerged, they have computational limitations or are designed to solve a specific task. We present a Transformer neural network that pre-trains task-agnostic sequence representations. This model is fine-tuned to solve two different protein prediction tasks: protein family classification and protein interaction prediction.
    [Show full text]
  • Predicting RNA Binding Proteins Using Discriminative Methods
    Predicting RNA binding proteins using discriminative methods by Shweta Bhandare M.S. Department of Interdisciplinary Telecommunications, University of Colorado at Boulder 2003 B.E., College of Engineering, University of Goa, 1995 A thesis submitted to the Faculty of the Graduate School of the University of Colorado in partial fulfillment of the requirements for the degree of Doctor of Philosophy Department of Computer Science 2017 This thesis entitled: Predicting RNA binding proteins using discriminative methods written by Shweta Bhandare has been approved for the Department of Computer Science Prof. Robin Dowell Dr. Debra S. Goldberg Prof. Larry E. Hunter Dr. Daniel Weaver Date The final copy of this thesis has been examined by the signatories, and we find that both the content and the form meet acceptable presentation standards of scholarly work in the above mentioned discipline. iii Bhandare, Shweta (Ph.D., Computer Science) Predicting RNA binding proteins using discriminative methods Thesis directed by Prof. Robin Dowell and Dr. Debra S. Goldberg This thesis examines the role of computational methods in identifying the motifs utilized by RNA-binding proteins (RBPs). RBPs play an important role in post-transcriptional regulation and identify their targets in a highly specific fashion through recognition of primary sequence and/or secondary structure hence making the prediction a complex problem. I applied the existing k-spectrum kernel method to a support vector machine and verified the published binding sites of two RBPs: Human antigen R (HuR) and Tristetraprolin (TTP). These RBPs exhibit opposing effects to the bound messenger RNA (mRNA) transcript but have simi- lar binding preferences.
    [Show full text]
  • Wolfson Review Wolfson The
    2012 – 2013 2013 No.37 – 2012 The Wolfson Review Wolfson The THE Wolfson Review 2012 – 2013 2013 No.37 – 2012 Wolfson College Barton Road Cambridge CB3 9BB www.wolfson.cam.ac.uk Upon 50 Years by John McClenahen (1986), Press Fellow The College is stone and mortar, and wood and glass. The College is ideas, great and small. The College is books and the Internet. Published in 2013 by Wolfson College, Cambridge The College is gates, gardens, paths, Barton Road, Cambridge CB3 9BB courts, and plaques, and a sundial. © Wolfson College, 2013 The College is Lee Library, the Dining Hall, and Bredon House. The College is students, and tutors, and Fellows. The College is fellowship, principles, and ritual. And the College, this College, Wolfson College, Cambridge, is much more. For this remarkable College is a diverse universe, ever expanding. From this College, in their diversity, those who study, guide, and reside here seek knowledge and truth in myriad ways. From this College, this special place, those who study, guide, and reside here seek to create, to find, to explore, to challenge, and to validate. Now and forever may their efforts – all our efforts wherever we are – Cover photograph Coloured primary hypothalamic neuronal culture, labelled ring true to the diversity and distinguishing humanness of this College, with MAP2, GFAP and Dapi under microscope, part of our young College in this ancient University. Wolfson Fellow Giles Yeo’s research into the brain control of food intake. Image created by Dr Brian Lam and Mr Joseph Polex-Wolf from the Yeo laboratory. The paper used for the Review contains material sourced from responsibly managed forests, certified in accordance with the Forestry Stewardship Council, and is printed using vegetable based inks.
    [Show full text]
  • New Tools for Old Organelles
    New tools for old organelles Bioinformatics for comparative cell biology Marc Gouw Dissertation presented to obtain the Ph.D degree in Bioinformatics Instituto de Tecnologia Química e Biológica António Xavier | Universidade Nova de Lisboa Oeiras, July, 2015 New tools for old organelles Bioinformatics for comparative cell biology Marc Gouw Dissertation presented to obtain the Ph.D degree in Evolutionary Biology Instituto de Tecnologia Química e Biológica António Xavier | Universidade Nova de Lisboa Research work coordinated by: Oeiras July, 2015 Acknowledgments This thesis and the years of work and inspiration put into it would not have been possible without the help of many people. First of all, I just generally want to thank everyone for being so great. Especially those who in any way contributed to making the past few years in Portugal and at the IGC such a wonderful experience: without all of you, my life and this thesis would have turned out very different indeed. A special thanks to all of the past and present members of the Computational Genomics Laboratory, who have always been kind enough to be unforgivingly critical of my work, and equally enthusiastic to go enjoy a beer on a Friday afternoon. And especially Beatriz Ferreira Gomes and Yoan Diekmann, who were essential to my PhD project: You have been truly inspiring collogues and close friends. A lot of this work would not have been possible without you. The epic mtoc-explorer.org, and indeed this thesis would also not have been possible without Renato Alves. I blame you, your programming genius and incredible patience for teaching me so much.
    [Show full text]
  • 1975.Txt 1. Computer Simulation of Protein Folding Michael Levitt
    1975.txt 1. Computer simulation of protein folding Michael Levitt* & Arieh Warshel*† †Department of Chemical Physics, Weizmann Institute of Science, Rehovoth, Israel *Present Address: MRC Laboratory of Molecular Biology, Hills Road, Cambridge, UK. A new and very simple representation of protein conformations has been used together with energy minimisation and thermalisation to simulate protein folding. Under certain conditions, the method succeeds in 'renaturing' bovine pancreatic trypsin inhibitor from an open-chain conformation into a folded conformation close to that of the native molecule. References 1.Anfinsen, C. B., Science, 181, 223–230 (1973). 2.Schulz, G. E., Barry, C. D., Friedman, J., Chou, P. Y., Fasman, G. D., Finkelstein, A. V., Lim, V. I., Ptitsyn, O. B., Kabat, E. A., Wu, T. T., Levitt, M., Robson, B., Nagano, K., Nature, 250, 140–142 (1974). 3.Scheraga, H. A., in Current topics in biochemistry (edit. by Anfinsen, C. B., and Schechter, A. N.), 1–42 (Academic, New York, 1974). 4.Ptitsyn, O. B., Vestnik Akad. Nauk S.S.S.R., 5, 57–68 (1973). 5.Flory, P. J. in Statistical Mechanics of Chain Molecules, 248–306 (Wiley, New York, 1969). 6.Hill, T. L., in Statistical Mechanics, 15–17 (McGraw-Hill, New York, 1956). 7.Nozaki, Y., and Tanford, C., J. biol. Chem., 246, 2211–2217 (1971). 8.Simon, E. M., Biopolymers, 10, 973–989 (1971). 9.Huber, R., Kukla, D., Ruhlmann, A., and Steigemann, W., Cold Spring Harbor Symp. quant. Biol., 36, 141–148 (1971). 10.Creighton, T. E., J. molec. Biol., 87, 603–624 (1974). 11.Phillips, D. C., in British Biochemistry, Past and Present (edit.
    [Show full text]
  • Weighted Graphlets and Deep Neural Networks for Protein Structure
    Weighted graphlets and deep neural networks for protein structure classification Hongyu Guo1, Khalique Newaz2,3,4, Scott Emrich5, Tijana Milenković2,3,4, and Jun Li1,* 1Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN 46556, USA 2Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN 46556, USA 3Eck Institute for Global Health, University of Notre Dame, Notre Dame, IN 46556, USA 4Interdisciplinary Center for Network Science & Applications, University of Notre Dame, Notre Dame, IN 46556, USA 5Department of Electrical Engineering & Computer Science, University of Tennessee at Knoxville, TN 37996, USA *To whom correspondence should be addressed. Tel: +1 574 631 3429; Fax: +1 574 631 4822; Email: [email protected] Abstract As proteins with similar structures often have similar functions, analysis of protein structures can help predict protein functions and is thus important. We consider the problem of protein structure classification, which computationally classifies the struc- tures of proteins into pre-defined groups. We develop a weighted network that depicts the protein structures, and more importantly, we propose the first graphlet-based mea- sure that applies to weighted networks. Further, we develop a deep neural network (DNN) composed of both convolutional and recurrent layers to use this measure for classification. Put together, our approach shows dramatic improvements in performance arXiv:1910.02594v1 [stat.ML] 7 Oct 2019 over existing graphlet-based approaches on 36 real datasets. Even comparing with the state-of-the-art approach, it almost halves the classification error. In addition to protein structure networks, our weighted-graphlet measure and DNN classifier can potentially be applied to classification of other weighted networks in computational biology as well as in other domains.
    [Show full text]
  • Multi-Domain Protein Families and Domain Pairs: Comparison with Known Structures and a Random Model of Domain Recombination
    Journal of Structural and Functional Genomics 4: 67–78, 2003. 67 © 2003 Kluwer Academic Publishers. Printed in the Netherlands. Multi-domain protein families and domain pairs: comparison with known structures and a random model of domain recombination Gordana Apic1, Wolfgang Huber2 & Sarah A. Teichmann1* 1MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 2QH, UK and 2DKFZ (German Cancer Research Center), Im Neuenheimer Feld 280, 69120 Heidelberg, Germany; * To whom correspondence should be addressed at [email protected] Received 26 November 2002; Accepted in final form 11 March 2003 Key words: structural genomics, multi-domain proteins, recombination, domain combinations Abstract There is a limited repertoire of domain families in nature that are duplicated and combined in different ways to form the set of proteins in a genome. Most proteins in both prokaryote and eukaryote genomes consist of two or more domains, and we show that the family size distribution of multi-domain protein families follows a power law like that of individual families. Most domain pairs occur in four to six different domain architectures: in isolation and in combinations with different partners. We showed previously that within the set of all pairwise domain combinations, most small and medium-sized families are observed in combination with one or two other families, while a few large families are very versatile and combine with many different partners. Though this may appear to be a stochastic pattern, in which large families have more combination partners by virtue of their size, we establish here that all the domain families with more than three members in genomes are duplicated more frequently than would be expected by chance considering their number of neighbouring domains.
    [Show full text]
  • Signature of Author
    An Analysis of Representations for Protein Structure Prediction by Barbara K. Moore Bryant S.B., S.M., Massachusetts Institute of Technology (1986) Submitted to the Department of Electrical Engineering and Computer Science in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY September 1994 ( Massachusetts Institute of Technology 1994 SignatureofAuthor ............. ,...... ........-.. .. Department of Electrical Engineering and Computer Science September 2, 1994 Certified by... ..... -- Patrick. H. Winston Professor, Electrical Engineering and Computer Science This Supervisor %-"ULfrtifioAvL.Lv--k &y ho, . Tomas o) no-Perez and Computer Science L! Thp.i. .,nPrvi.qnr Accepted by .. ... , a. ,- · F' drY................... / Frederic~) R. Morgenthaler Departmental /Committee on Graduate Students Iv:r tfF:r _ ;: ! --.. t.8r'(itr i-7etFr An Analysis of Representations for Protein Structure Prediction by Barbara K. Moore Bryant Submitted to the Department of Electrical Engineering and Computer Science on September 2, 1994, in partial fulfillment of the requirements for the degree of Doctor of Philosophy Abstract Protein structure prediction is a grand challenge in' the fields of biology and computer science. Being able to quickly determine the structure of a protein from its amino acid sequence would be extremely useful to biologists interested in elucidating the mechanisms of life and finding ways to cure disease. In spite of a wealth of knowl- edge about proteins and their structure, the structure prediction problem has gone unsolved in the nearly forty years since the first determination of a protein structure by X-ray crystallography. In this thesis, I discuss issues in the representation of protein structure and se- quence for algorithms which perform structure prediction.
    [Show full text]