EXPANSION: a webserver to explore the functional consequences of protein-coding alternative splice variants in cancer genomics

Summary

EXPANSION (https://expansion.bioinfolab.sns.it/) is an integrated web-server to explore the functional consequences of protein-coding alternative splice variants. We combined information from Differentially Expressed (DE) protein-coding transcripts from cancer genomics, together with domain architecture, protein interaction network, and gene enrichment analysis to provide an easy-to-interpret view of the effects of protein-coding splice variants. We retrieved all the protein-coding Ensembl transcripts and mapped Interpro domains and post-translational modifications on canonical sequences to identify functionally relevant splicing events. We also retrieved isoform-specific protein-protein interactions and binding regions from IntAct to uncover isoform-specific functions via gene-set over-representation analysis. Through EXPANSION, users can analyze precalculated or user-inputted DE transcript datasets, to easily gain functional insights on any protein spliceform of interest.

Availability and Implementation

EXPANSION is freely available at http://expansion.bioinfolab.sns.it/. The code of the scripts used for EXPASION is available at: https://github.com/raimondilab/expansion. Datasets associated to this resource are available at the following URL: https://doi.org/10.5281/zenodo.8229120. The web-server was developed using Apache2 (https://https.apache.org/) and Flask (v2.0.2) (http://flask.pocoo.org/) for the web frontend and for the internal pipeline to handle back-end processes. We additionally used the following Python and JavaScript libraries at both back- and front-ends: D3 (v4), jQuery (v3.2.1), DataTables (v2.3.2), biopython (v1.79), gprofiler-officia l(v1.0.0), Mysql-connector-python (v8.0.31). To construct the API, Fast API library (v0.95.1) was used.

GPCRome-wide analysis of G-protein-coupling diversity using a computational biology approach

GPCRs are master regulators of cell signaling by transducing extracellular stimuli into the cell via selective coupling to intracellular G-proteins. Here we present a computational analysis of the structural determinants of G-protein-coupling repertoire of experimental and predicted 3D GPCR-G-protein complexes. Interface contact analysis recapitulates structural hallmarks associated with G-protein-coupling specificity, including TM5, TM6 and ICLs. We employ interface contacts as fingerprints to cluster G_s vs G_i complexes in an unsupervised fashion, suggesting that interface residues contribute to selective coupling. We experimentally confirm on a promiscuous receptor (CCKAR) that mutations of some of these specificity-determining positions bias the coupling selectivity. Interestingly, G_s-GPCR complexes have more conserved interfaces, while G_i/o proteins adopt a wider number of alternative docking poses, as assessed via structural alignments of representative 3D complexes. Binding energy calculations demonstrate that distinct structural properties of the complexes are associated to higher stability of G_s than G_i/o complexes. AlphaFold2 predictions of experimental binary complexes confirm several of these structural features and allow us to augment the structural coverage of poorly characterized complexes such as G_12/13.

PRECOGx: exploring GPCR signaling mechanisms with deep protein representations

In this study we show that protein language models can encode structural and functional information of GPCR sequences that can be used to predict their signaling and functional repertoire. We used the ESM1b protein embeddings as features and the binding information known from publicly available studies to develop PRECOGx, a machine learning predictor to explore GPCR interactions with G protein and β-arrestin, which we made available through a new webserver (https://precogx.bioinfolab.sns.it/). PRECOGx outperformed its predecessor (e.g. PRECOG) in predicting GPCR-transducer couplings, being also able to consider all GPCR classes. The webserver also provides new functionalities, such as the projection of input sequences on a low-dimensional space describing essential features of the human GPCRome, which is used as a reference to track GPCR variants. Additionally, it allows inspection of the sequence and structural determinants responsible for coupling via the analysis of the most important attention maps used by the models as well as through predicted intramolecular contacts. We demonstrate applications of PRECOGx by predicting the impact of disease variants (ClinVar) and alternative splice forms from healthy tissues (GTEX) of human GPCRs, revealing the power to dissect system biasing mechanisms in both health and disease.