Illuminating G-Protein-Coupling Selectivity of GPCRs.

Heterotrimetic G proteins consist of four subfamilies (Gs, Gi/o, Gq/11, and G12/13) that mediate signaling via G-protein-coupled receptors (GPCRs), principally by receptors binding Gα C termini. G-protein-coupling profiles govern GPCR-induced cellular responses, yet receptor sequence selectivity determinants remain elusive. Here, we systematically quantified ligand-induced interactions between 148 GPCRs and all 11 unique Gα subunit C termini. For each receptor, we probed chimeric Gα subunit activation via a transforming growth factor-α (TGF-α) shedding response in HEK293 cells lacking endogenous Gq/11 and G12/13 proteins, and complemented G-protein-coupling profiles through a NanoBiT-G-protein dissociation assay. Interrogation of the dataset identified sequence-based coupling specificity features, inside and outside the transmembrane domain, which we used to develop a coupling predictor that outperforms previous methods. We used the predictor to engineer designer GPCRs selectively coupled to G12. This dataset of fine-tuned signaling mechanisms for diverse GPCRs is a valuable resource for research in GPCR signaling.

Biochemically validated structural model of the 15-subunit intraflagellar transport complex IFT-B.

Cilia are ubiquitous eukaryotic organelles impotant for cellular motility, signaling, and sensory reception. Cilium formation requires intraflagellar transport of structural and signaling components and involves 22 different proteins organized into intraflagellar transport (IFT) complexes IFT-A and IFT-B that are transported by molecular motors. The IFT-B complex constitutes the backbone of polymeric IFT trains carrying cargo between the cilium and the cell body. Currently, high-resolution structures are only available for smaller IFT-B subcomplexes leaving > 50% structurally uncharacterized. Here, we used Alphafold to structurally model the 15-subunit IFT-B complex. The model was validated using cross-linking/mass-spectrometry data on reconstituted IFT-B complexes, X-ray scattering in solution, diffraction from crystals as well as site-directed mutagenesis and protein-binding assays. The IFT-B structure reveals an elongated and highly flexible complex consistent with cryo-electron tomographic reconstructions of IFT trains. The IFT-B complex organizes into IFT-B1 and IFT-B2 parts with binding sites for ciliary cargo and the inactive IFT dynein motor, respectively. Interestingly, our results are consistent with two different binding sites for IFT81/74 on IFT88/70/52/46 suggesting the possibility of different structural architectures for the IFT-B1 complex. Our data present a structural framework to understand IFT-B complex assembly, function, and ciliopathy variants.

DLG1 functions upstream of SDCCAG3 and IFT20 to control ciliary targeting of polycystin-2.

Polarized vesicular trafficking directs specific receptors and ion channels to cilia, but the underlying mechanisms are poorly understood. Here we describe a role for DLG1, a core component of the Scribble polarity complex, in regulating ciliary protein trafficking in kidney epithelial cells. Conditional knockout of Dlg1 in mouse kidney causes ciliary elongation and cystogenesis, and cell-based proximity labeling proteomics and fluorescence microscopy show alterations in the ciliary proteome upon loss of DLG1. Specifically, the retromer-associated protein SDCCAG3, IFT20, and polycystin-2 (PC2) are reduced in the cilia of DLG1-deficient cells compared to control cells. This phenotype is recapitulated in vivo and rescuable by re-expression of wild-type DLG1, but not a Congenital Anomalies of the Kidney and Urinary Tract (CAKUT)-associated DLG1 variant, p.T489R. Finally, biochemical approaches and Alpha Fold modelling suggest that SDCCAG3 and IFT20 form a complex that associates, at least indirectly, with DLG1. Our work identifies a key role for DLG1 in regulating ciliary protein composition and suggests that ciliary dysfunction of the p.T489R DLG1 variant may contribute to CAKUT.

WeSA: a web server for improving analysis of affinity proteomics data.

Protein-protein interaction experiments still yield many false positive interactions. The socioaffinity metric can distinguish true protein-protein interactions from noise based on available data. Here, we present WeSA (Weighted SocioAffinity), which considers large datasets of interaction proteomics data (IntAct, BioGRID, the BioPlex) to score human protein interactions and, in a statistically robust way, flag those (even from a single experiment) that are likely to be false positives. ROC analysis (using CORUM-PDB positives and Negatome negatives) shows that WeSA improves over other measures of interaction confidence. WeSA shows consistently good results over all datasets (up to: AUC = 0.93 and at best threshold: TPR = 0.84, FPR = 0.11, Precision = 0.98). WeSA is freely available without login (wesa.russelllab.org). Users can submit their own data or look for organized information on human protein interactions using the web server. Users can either retrieve available information for a list of proteins of interest or calculate scores for new experiments. The server outputs either pre-computed or updated WeSA scores for the input enriched with information from databases. The summary is presented as a table and a network-based visualization allowing the user to remove those nodes/edges that the method considers spurious.

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.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

Genomic landscape of follicular lymphoma across a wide spectrum of clinical behaviors.

While some follicular lymphoma (FL) patients do not require treatment or experience prolonged responses, others relapse early, and little is known about genetic alterations specific to patients with a particular clinical behavior. We selected 56 grade 1-3A FL patients according to their need of treatment or timing of relapse: never treated (n = 7), non-relapsed (19), late relapse (14), early relapse or POD24 (11), and primary refractory (5). We analyzed 56 diagnostic and 12 paired relapse lymphoid tissue biopsies and performed copy number alteration (CNA) analysis and next generation sequencing (NGS). We identified six focal driver losses (1p36.32, 6p21.32, 6q14.1, 6q23.3, 9p21.3, 10q23.33) and 1p36.33 copy-neutral loss of heterozygosity (CN-LOH). By integrating CNA and NGS results, the most frequently altered genes/regions were KMT2D (79%), CREBBP (67%), TNFRSF14 (46%) and BCL2 (40%). Although we found that mutations in PIM1, FOXO1 and TMEM30A were associated with an adverse clinical behavior, definitive conclusions cannot be drawn, due to the small sample size. We identified common precursor cells harboring early oncogenic alterations of the KMT2D, CREBBP, TNFRSF14 and EP300 genes and 16p13.3-p13.2 CN-LOH. Finally, we established the functional consequences of mutations by means of protein modeling (CD79B, PLCG2, PIM1, MCL1 and IRF8). These data expand the knowledge on the genomics behind the heterogeneous FL population and, upon replication in larger cohorts, could contribute to risk stratification and the development of targeted therapies.

Mechnetor: a web server for exploring protein mechanism and the functional context of genetic variants.

Advances in DNA sequencing and proteomics mean that researchers must now regularly interrogate thousands of positional gene/protein changes in order to find those relevant for potential clinical application or biological insights. The abundance of already known information on protein interactions, mechanism, and tertiary structure provides the possible means to understand these changes rapidly, though a careful and systematic integration of these diverse datasets is first needed. For this purpose, we developed Mechnetor, a tool that allows users to quickly explore and visualize integrated mechanistic data for proteins or interactions of interest. Central to the system is a careful cataloguing of diverse sources of protein interaction mechanism, and an efficient means to visualize interactions between relevant and/or known protein regions. The result is a finer resolution interaction network that provides more immediate clues as to points of intervention or mechanistic understanding. Users can import protein, interactions, genetic variants or post-translational modifications and see these data in the best known mechanistic context. We demonstrate the tool with topical examples in human genetic diseases and cancer genomics. The tool is freely available at: mechnetor.russelllab.org.

The mutational landscape of Burkitt-like lymphoma with 11q aberration is distinct from that of Burkitt lymphoma.

The new recently described provisional lymphoma category Burkitt-like lymphoma with 11q aberration comprises cases similar to Burkitt lymphoma (BL) on morphological, immunophenotypic and gene-expression levels but lacking the IG-MYC translocation. They are characterized by a peculiar imbalance pattern on chromosome 11, but the landscape of mutations is not yet described. Thus, we investigated 15 MYC-negative Burkitt-like lymphoma with 11q aberration (mnBLL,11q,) cases by copy-number analysis and whole-exome sequencing. We refined the regions of 11q imbalance and identified the INO80 complex-associated gene NFRKB as a positional candidate in 11q24.3. Next to recurrent gains in 12q13.11-q24.32 and 7q34-qter as well as losses in 13q32.3-q34, we identified 47 genes recurrently affected by protein-changing mutations (each ≥3 of 15 cases). Strikingly, we did not detect recurrent mutations in genes of the ID3-TCF3 axis or the SWI/SNF complex that are frequently altered in BL, or in genes frequently mutated in germinal center-derived B-cell lymphomas like KMT2D or CREBBP An exception is GNA13, which was mutated in 7 of 15 cases. We conclude that the genomic landscape of mnBLL,11q, differs from that of BL both at the chromosomal and mutational levels. Our findings implicate that mnBLL,11q, is a lymphoma category distinct from BL at the molecular level.

PRECOG: PREdicting COupling probabilities of G-protein coupled receptors.

G-protein coupled receptors (GPCRs) control multiple physiological states by transducing a multitude of extracellular stimuli into the cell via coupling to intra-cellular heterotrimeric G-proteins. Deciphering which G-proteins couple to each of the hundreds of GPCRs present in a typical eukaryotic organism is therefore critical to understand signalling. Here, we present PRECOG (precog.russelllab.org): a web-server for predicting GPCR coupling, which allows users to: (i) predict coupling probabilities for GPCRs to individual G-proteins instead of subfamilies; (ii) visually inspect the protein sequence and structural features that are responsible for a particular coupling; (iii) suggest mutations to rationally design artificial GPCRs with new coupling properties based on predetermined coupling features.

Illuminating the Onco-GPCRome: Novel G protein-coupled receptor-driven oncocrine networks and targets for cancer immunotherapy.

G protein-coupled receptors (GPCRs) are the largest gene family of cell membrane-associated molecules mediating signal transmission, and their involvement in key physiological functions is well-established. The ability of GPCRs to regulate a vast array of fundamental biological processes, such as cardiovascular functions, immune responses, hormone and enzyme release from endocrine and exocrine glands, neurotransmission, and sensory perception (e.g. vision, odor, and taste), is largely due to the diversity of these receptors and the layers of their downstream signaling circuits. Dysregulated expression and aberrant functions of GPCRs have been linked to some of the most prevalent human diseases, which renders GPCRs one of the top targets for pharmaceutical drug development. However, the study of the role of GPCRs in tumor biology has only just begun to make headway. Recent studies have shown that GPCRs can contribute to the many facets of tumorigenesis, including proliferation, survival, angiogenesis, invasion, metastasis, therapy resistance, and immune evasion. Indeed, GPCRs are widely dysregulated in cancer and yet are underexploited in oncology. We present here a comprehensive analysis of GPCR gene expression, copy number variation, and mutational signatures in 33 cancer types. We also highlight the emerging role of GPCRs as part of oncocrine networks promoting tumor growth, dissemination, and immune evasion, and we stress the potential benefits of targeting GPCRs and their signaling circuits in the new era of precision medicine and cancer immunotherapies.

Defining clinical subgroups and genotype-phenotype correlations in NBAS-associated disease across 110 patients.

PURPOSE: Pathogenic variants in neuroblastoma-amplified sequence (NBAS) cause an autosomal recessive disorder with a wide range of symptoms affecting liver, skeletal system, and brain, among others. There is a continuously growing number of patients but a lack of systematic and quantitative analysis. METHODS: Individuals with biallelic variants in NBAS were recruited within an international, multicenter study, including novel and previously published patients. Clinical variables were analyzed with log-linear models and visualized by mosaic plots; facial profiles were investigated via DeepGestalt. The structure of the NBAS protein was predicted using computational methods. RESULTS: One hundred ten individuals from 97 families with biallelic pathogenic NBAS variants were identified, including 26 novel patients with 19 previously unreported variants, giving a total number of 86 variants. Protein modeling redefined the ß-propeller domain of NBAS. Based on the localization of missense variants and in-frame deletions, three clinical subgroups arise that differ significantly regarding main clinical features and are directly related to the affected region of the NBAS protein: ß-propeller (combined phenotype), Sec39 (infantile liver failure syndrome type 2/ILFS2), and C-terminal (short stature, optic atrophy, and Pelger-Huët anomaly/SOPH). CONCLUSION: We define clinical subgroups of NBAS-associated disease that can guide patient management and point to domain-specific functions of NBAS.

Landscape of nuclear transport receptor cargo specificity.

Nuclear transport receptors (NTRs) recognize localization signals of cargos to facilitate their passage across the central channel of nuclear pore complexes (NPCs). About 30 different NTRs constitute different transport pathways in humans and bind to a multitude of different cargos. The exact cargo spectrum of the majority of NTRs, their specificity and even the extent to which active nucleocytoplasmic transport contributes to protein localization remains understudied because of the transient nature of these interactions and the wide dynamic range of cargo concentrations. To systematically map cargo-NTR relationships in situ, we used proximity ligation coupled to mass spectrometry (BioID). We systematically fused the engineered biotin ligase BirA* to 16 NTRs. We estimate that a considerable fraction of the human proteome is subject to active nuclear transport. We quantified the specificity and redundancy in NTR interactions and identified transport pathways for cargos. We extended the BioID method by the direct identification of biotinylation sites. This approach enabled us to identify interaction interfaces and to discriminate direct versus piggyback transport mechanisms. Data are available via ProteomeXchange with identifier PXD007976.

Capturing protein communities by structural proteomics in a thermophilic eukaryote.

The arrangement of proteins into complexes is a key organizational principle for many cellular functions. Although the topology of many complexes has been systematically analyzed in isolation, their molecular sociology in situ remains elusive. Here, we show that crude cellular extracts of a eukaryotic thermophile, Chaetomium thermophilum, retain basic principles of cellular organization. Using a structural proteomics approach, we simultaneously characterized the abundance, interactions, and structure of a third of the C. thermophilum proteome within these extracts. We identified 27 distinct protein communities that include 108 interconnected complexes, which dynamically associate with each other and functionally benefit from being in close proximity in the cell. Furthermore, we investigated the structure of fatty acid synthase within these extracts by cryoEM and this revealed multiple, flexible states of the enzyme in adaptation to its association with other complexes, thus exemplifying the need for in situ studies. As the components of the captured protein communities are known-at both the protein and complex levels-this study constitutes another step forward toward a molecular understanding of subcellular organization.

Systematic identification of phosphorylation-mediated protein interaction switches.

Proteomics techniques can identify thousands of phosphorylation sites in a single experiment, the majority of which are new and lack precise information about function or molecular mechanism. Here we present a fast method to predict potential phosphorylation switches by mapping phosphorylation sites to protein-protein interactions of known structure and analysing the properties of the protein interface. We predict 1024 sites that could potentially enable or disable particular interactions. We tested a selection of these switches and showed that phosphomimetic mutations indeed affect interactions. We estimate that there are likely thousands of phosphorylation mediated switches yet to be uncovered, even among existing phosphorylation datasets. The results suggest that phosphorylation sites on globular, as distinct from disordered, parts of the proteome frequently function as switches, which might be one of the ancient roles for kinase phosphorylation.

Dissecting the genomic complexity underlying medulloblastoma.

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

Mechismo: predicting the mechanistic impact of mutations and modifications on molecular interactions.

Systematic interrogation of mutation or protein modification data is important to identify sites with functional consequences and to deduce global consequences from large data sets. Mechismo (mechismo.russellab.org) enables simultaneous consideration of thousands of 3D structures and biomolecular interactions to predict rapidly mechanistic consequences for mutations and modifications. As useful functional information often only comes from homologous proteins, we benchmarked the accuracy of predictions as a function of protein/structure sequence similarity, which permits the use of relatively weak sequence similarities with an appropriate confidence measure. For protein-protein, protein-nucleic acid and a subset of protein-chemical interactions, we also developed and benchmarked a measure of whether modifications are likely to enhance or diminish the interactions, which can assist the detection of modifications with specific effects. Analysis of high-throughput sequencing data shows that the approach can identify interesting differences between cancers, and application to proteomics data finds potential mechanistic insights for how post-translational modifications can alter biomolecular interactions.

Mutagenesis of N-terminal residues of feline foamy virus Gag reveals entirely distinct functions during capsid formation, particle assembly, Gag processing and budding.

BACKGROUND: Foamy viruses (FVs) of the Spumaretrovirinae subfamily are distinct retroviruses, with many features of their molecular biology and replication strategy clearly different from those of the Orthoretroviruses, such as human immunodeficiency, murine leukemia, and human T cell lymphotropic viruses. The FV Gag N-terminal region is responsible for capsid formation and particle budding via interaction with Env. However, the critical residues or motifs in this region and their functional interaction are currently ill-defined, especially in non-primate FVs. RESULTS: Mutagenesis of N-terminal Gag residues of feline FV (FFV) reveals key residues essential for either capsid assembly and/or viral budding via interaction with the FFV Env leader protein (Elp). In an in vitro Gag-Elp interaction screen, Gag mutations abolishing particle assembly also interfered with Elp binding, indicating that Gag assembly is a prerequisite for this highly specific interaction. Gradient sedimentation analyses of cytosolic proteins indicate that wild-type Gag is mostly assembled into virus capsids. Moreover, proteolytic processing of Gag correlates with capsid assembly and is mostly, if not completely, independent from particle budding. In addition, Gag processing correlates with the presence of packaging-competent FFV genomic RNA suggesting that Pol encapsidation via genomic RNA is a prerequisite for Gag processing. Though an appended heterogeneous myristoylation signal rescues Gag particle budding of mutants unable to form capsids or defective in interacting with Elp, it fails to generate infectious particles that co-package Pol, as evidenced by a lack of Gag processing. CONCLUSIONS: Changes in proteolytic Gag processing, intracellular capsid assembly, particle budding and infectivity of defined N-terminal Gag mutants highlight their essential, distinct and only partially overlapping roles during viral assembly and budding. Discussion of these findings will be based on a recent model developed for Gag-Elp interactions in prototype FV.

Genes encoding members of the JAK-STAT pathway or epigenetic regulators are recurrently mutated in T-cell prolymphocytic leukaemia.

T-cell prolymphocytic leukaemia (T-PLL) is an aggressive leukaemia. The primary genetic alteration in T-PLL are the inv(14)(q11q32)/t(14;14)(q11;q32) leading to TRD/TRA-TCL1A fusion, or the t(X;14)(q28;q11) associated with TRD/TRA-MTCP1 fusion. However, additional cooperating abnormalities are necessary for emergence of the full neoplastic phenotype. Though the pattern of secondary chromosomal aberrations is remarkably conserved, targets of the changes are largely unknown. We analysed a cohort of 43 well-characterized T-PLL for hotspot mutations in the genes JAK3, STAT5B and RHOA. Additionally, we selected a subset of 23 T-PLL cases for mutational screening of 54 genes known to be recurrently mutated in T-cell and other haematological neoplasms. Activating mutations in the investigated regions of the JAK3 and STAT5B genes were detected in 30% (13/43) and 21% (8/39) of the cases, respectively, and were mutually exclusive. Further, we identified mutations in the genes encoding the epigenetic regulators EZH2 in 13% (3/23), TET2 in 17% (4/23) and BCOR in 9% (2/23) of the cases. We confirmed that the JAK-STAT pathway is a major mutational target, and identified epigenetic regulators recurrently mutated in T-PLL. These findings complement the mutational spectrum of secondary aberrations in T-PLL and underscore the potential therapeutical relevance of epigenetic regulators in T-PLL.

A recurrent 11q aberration pattern characterizes a subset of MYC-negative high-grade B-cell lymphomas resembling Burkitt lymphoma.

The genetic hallmark of Burkitt lymphoma (BL) is the t(8;14)(q24;q32) and its variants leading to activation of the MYC oncogene. It is a matter of debate whether true BL without MYC translocation exists. Here, we identified 59 lymphomas concordantly called BL by 2 gene expression classifiers among 753 B-cell lymphomas. Only 2 (3%) of these 59 molecular BL lacked a MYC translocation, which both shared a peculiar pattern of chromosome 11q aberration characterized by interstitial gains including 11q23.2-q23.3 and telomeric losses of 11q24.1-qter. We extended our analysis to 17 MYC-negative high-grade B-cell lymphomas with a similar 11q aberration and showed this aberration to be recurrently associated with morphologic and clinical features of BL. The minimal region of gain was defined by high-level amplifications in 11q23.3 and associated with overexpression of genes including PAFAH1B2 on a transcriptional and protein level. The recurrent region of loss contained a focal homozygous deletion in 11q24.2-q24.3 including the ETS1 gene, which was shown to be mutated in 4 of 16 investigated cases. These findings indicate the existence of a molecularly distinct subset of B-cell lymphomas reminiscent of BL, which is characterized by deregulation of genes in 11q.