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.

We explored the dysregulation of G-protein-coupled receptor (GPCR) ligand systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes. Multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes and are associated to specific transcriptional programs and to patient survival patterns. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including, e.g., muscarinic, adenosine, 5-hydroxytryptamine, and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens, which we further validated. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

Immune checkpoint blockade (ICB) targeting PD-1 and CTLA-4 has revolutionized cancer treatment. However, many cancers do not respond to ICB, prompting the search for additional strategies to achieve durable responses. G-protein-coupled receptors (GPCRs) are the most intensively studied drug targets but are underexplored in immuno-oncology. Here, we cross-integrated large singe-cell RNA-sequencing datasets from CD8⁺ T cells covering 19 distinct cancer types and identified an enrichment of Gα_s-coupled GPCRs on exhausted CD8⁺ T cells. These include EP₂, EP₄, A_2AR, β₁AR and β₂AR, all of which promote T cell dysfunction. We also developed transgenic mice expressing a chemogenetic CD8-restricted Gα_s-DREADD to activate CD8-restricted Gα_s signaling and show that a Gα_s-PKA signaling axis promotes CD8⁺ T cell dysfunction and immunotherapy failure. These data indicate that Gα_s-GPCRs are druggable immune checkpoints that might be targeted to enhance the response to ICB immunotherapies.