Skip to content

This code was written and used to perform the comparative interactome analysis of β-arrestin1 and 2.

License

Notifications You must be signed in to change notification settings

mo-yoda/Haider_Review

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

Comparative Analysis of the Interactome of both β-Arrestin Isoforms

Conformational flexibility of β-arrestins – how these scaffolding proteins guide and transform the functionality of GPCRs

Raphael Silvanus Haider*[1,2] , Mona Reichel*[1], Edda Sofie Fabienne Matthees[1], Carsten Hoffmann[1]

[1]: Institut für Molekulare Zellbiologie, CMB – Center for Molecular Biomedicine, Universitätsklinikum Jena, Friedrich-Schiller-Universität Jena, Hans-Knöll Straße 2, D-07745 Jena, Germany

[2]: Division of Physiology, Pharmacology and Neuroscience, School of Life Sciences, Queen's Medical Centre, University of Nottingham, Nottingham, UK

[*] contributed equally


This code was written and used to perform the comparative interactome analysis of β-arrestin1 and 2 found in DOI.


Resource data

Details about enriched GO term clustering


Short description

Interacting proteins were retrieved from STRINGdb 11.5 (Szklarczyk et al. 2018). GO enrichment analysis was performed using the online tool Database for Annotation, Visualization and Integrated Discovery (DAVID) 6.8 (Huang et al. 2009, Sherman et al. 2022). Enriched GO terms were clustered and visualized using the R/Bioconductor package simplifyEnrichment (Gu & Huebschmann 2022).


Workflow

1. Retrieving interacting proteins of β-arrestin1 and 2 from STRINGdb

  • retrieves all proteins interacting with β-arrestin1 or 2 from STRINGdb using STRINGdb API
  • entries with confidence score < 0.5 are removed
  • returns xlsx file with ENSP of interacting proteins, confidence score of this interaction and the "uniqueness" of the interaction (whether this interaction is unique to one of the β-arrestin isoforms)
  • splits xlsx file based on the "uniqueness" column
  • returns xlsx files for β-arrestin1, β-arrestin2 and interactors of both isoforms

(Note, that proteins interacting with both isoforms are found twice in these exports)


2. Translating ENSP IDs to uniprot accession number

  • translates ENSP IDs exported from STRINGdb to uniprot accession number using the uniprot ID mapping tool via the API
  • Note, that order of fetched IDs do not match order of requested IDs - to prevent mistranslation, the script creates a translation dictionary which is then applied to the imported ENSP IDs
  • returns input xlsx with additional column containing uniprot accession IDs

3. Extract information from uniprot entries

  • uniprot accession number duplicates in input table are removed (which are all proteins interacting with both β-arrestins as described above)
  • using the uniprot accession number, selected information from each uniprot entry is retrieved via APIs (to test how certain information can be accessed from xml entry, test.xml was created)
  • according to the uniprot Keyword 'G-protein coupled receptor', resulting table is separated into interactors which are GPCRs and non-GPCRs

is unused for analysis, but gave an initial overview

  • retrieves GO terms of all three ontologies from uniprot using APIs
  • lists each GO term separately with corresponding proteins

4. GO enrichment analysis via DAVID

GO enrichment analysis of the biological process ontology was performed using the online tool DAVID 6.8. Uniprot accession numbers of non-GPCR proteins obtained from extract_uniprot_info.py were used as input.

In detail, functional annotation analysis was performed using the uniprot accession IDs of non-GPCR proteins interacting with exclusively one β-arrestin isoform (gene list input) against all non-GPCR interacting proteins of this isoform (background).The functional annotation chart created by DAVID was downloaded and used in the last step visualizeEnrichment.R.


5. Clustering and visualization of enriched GO terms

  • uses functional annotation chart created by DAVID 6.8 as input

  • enriched GO terms for β-arrestin1-only or β-arrestin2-only interacting proteins are clustered according to their similarity and displayed as heatmap using simplifyGOFromMultipleLists() of the R/Bioconductor package simplifyEnrichment (Gu & Huebschmann 2022)

  • GO terms found in each cluster are exported as xlsx file

  • additionally, proteins which contributed to each cluster due to their GO term assignment are exported as xlsx file

    (Note, that one protein can contribute to several clusters as each protein has several GO terms assigned)


This project is licensed under the terms of the MIT license.