The Cancermuts software package for the prioritization of missense cancer variants: a case study of AMBRA1 in melanoma

AMBRA1 levels predict resistance to MAPK inhibitors in melanoma

Intrinsic and acquired resistance to mitogen-activated protein kinase inhibitors (MAPKi) in melanoma remains a major therapeutic challenge. Here, we show that the clinical development of resistance to MAPKi is associated with reduced tumor expression of the melanoma suppressor Autophagy and Beclin 1 Regulator 1 (AMBRA1) and that lower expression levels of AMBRA1 predict a poor response to MAPKi treatment. Functional analyses show that loss of AMBRA1 induces phenotype switching and orchestrates an extracellular signal-regulated kinase (ERK)-independent resistance mechanism by activating focal adhesion kinase 1 (FAK1). In both in vitro and in vivo settings, melanomas with low AMBRA1 expression exhibit intrinsic resistance to MAPKi therapy but higher sensitivity to FAK1 inhibition. Finally, we show that the rapid development of resistance in initially MAPKi-sensitive melanomas can be attributed to preexisting subclones characterized by low AMBRA1 expression and that cotreatment with MAPKi and FAK1 inhibitors (FAKi) effectively prevents the development of resistance in these tumors. In summary, our findings underscore the value of AMBRA1 expression for predicting melanoma response to MAPKi and supporting the therapeutic efficacy of FAKi to overcome MAPKi-induced resistance.

Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks

The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.

Structure of the DDB1-AMBRA1 E3 ligase receptor complex linked to cell cycle regulation

AMBRA1 is a tumor suppressor protein that functions as a substrate receptor of the ubiquitin conjugation system with roles in autophagy and the cell cycle regulatory network. The intrinsic disorder of AMBRA1 has thus far precluded its structural determination. To solve this problem, we analyzed the dynamics of AMBRA1 using hydrogen deuterium exchange mass spectrometry (HDX-MS). The HDX results indicated that AMBRA1 is a highly flexible protein and can be stabilized upon interaction with DDB1, the adaptor of the Cullin4A/B E3 ligase. Here, we present the cryo-EM structure of AMBRA1 in complex with DDB1 at 3.08 Å resolution. The structure shows that parts of the N- and C-terminal structural regions in AMBRA1 fold together into the highly dynamic WD40 domain and reveals how DDB1 engages with AMBRA1 to create a binding scaffold for substrate recruitment. The N-terminal helix-loop-helix motif and WD40 domain of AMBRA1 associate with the double-propeller fold of DDB1. We also demonstrate that DDB1 binding-defective AMBRA1 mutants prevent ubiquitination of the substrate Cyclin D1 in vitro and increase cell cycle progression. Together, these results provide structural insights into the AMBRA1-ubiquitin ligase complex and suggest a mechanism by which AMBRA1 acts as a hub involved in various physiological processes.

Computational analysis of five neurodegenerative diseases reveals shared and specific genetic loci

Neurodegenerative diseases (ND) are heterogeneous disorders of the central nervous system that share a chronic and selective process of neuronal cell death. A computational approach to investigate shared genetic and specific loci was applied to 5 different ND: Amyotrophic lateral sclerosis (ALS), Alzheimer's disease (AD), Parkinson's disease (PD), Multiple sclerosis (MS), and Lewy body dementia (LBD). The datasets were analyzed separately, and then we compared the obtained results. For this purpose, we applied a genetic correlation analysis to genome-wide association datasets and revealed different genetic correlations with several human traits and diseases. In addition, a clumping analysis was carried out to identify SNPs genetically associated with each disease. We found 27 SNPs in AD, 6 SNPs in ALS, 10 SNPs in PD, 17 SNPs in MS, and 3 SNPs in LBD. Most of them are located in non-coding regions, with the exception of 5 SNPs on which a protein structure and stability prediction was performed to verify their impact on disease. Furthermore, an analysis of the differentially expressed miRNAs of the 5 examined pathologies was performed to reveal regulatory mechanisms that could involve genes associated with selected SNPs. In conclusion, the results obtained constitute an important step toward the discovery of diagnostic biomarkers and a better understanding of the diseases.

TRAP1 S-nitrosylation as a model of population-shift mechanism to study the effects of nitric oxide on redox-sensitive oncoproteins

S-nitrosylation is a post-translational modification in which nitric oxide (NO) binds to the thiol group of cysteine, generating an S-nitrosothiol (SNO) adduct. S-nitrosylation has different physiological roles, and its alteration has also been linked to a growing list of pathologies, including cancer. SNO can affect the function and stability of different proteins, such as the mitochondrial chaperone TRAP1. Interestingly, the SNO site (C501) of TRAP1 is in the proximity of another cysteine (C527). This feature suggests that the S-nitrosylated C501 could engage in a disulfide bridge with C527 in TRAP1, resembling the well-known ability of S-nitrosylated cysteines to resolve in disulfide bridge with vicinal cysteines. We used enhanced sampling simulations and in-vitro biochemical assays to address the structural mechanisms induced by TRAP1 S-nitrosylation. We showed that the SNO site induces conformational changes in the proximal cysteine and favors conformations suitable for disulfide bridge formation. We explored 4172 known S-nitrosylated proteins using high-throughput structural analyses. Furthermore, we used a coarse-grained model for 44 protein targets to account for protein flexibility. This resulted in the identification of up to 1248 proximal cysteines, which could sense the redox state of the SNO site, opening new perspectives on the biological effects of redox switches. In addition, we devised two bioinformatic workflows ( https://github.com/ELELAB/SNO_investigation_pipelines ) to identify proximal or vicinal cysteines for a SNO site with accompanying structural annotations. Finally, we analyzed mutations in tumor suppressors or oncogenes in connection with the conformational switch induced by S-nitrosylation. We classified the variants as neutral, stabilizing, or destabilizing for the propensity to be S-nitrosylated and undergo the population-shift mechanism. The methods applied here provide a comprehensive toolkit for future high-throughput studies of new protein candidates, variant classification, and a rich data source for the research community in the NO field.

Ambra1 modulates the tumor immune microenvironment and response to PD-1 blockade in melanoma

BACKGROUND

Loss of Ambra1 (autophagy and beclin 1 regulator 1), a multifunctional scaffold protein, promotes the formation of nevi and contributes to several phases of melanoma development. The suppressive functions of Ambra1 in melanoma are mediated by negative regulation of cell proliferation and invasion; however, evidence suggests that loss of Ambra1 may also affect the melanoma microenvironment. Here, we investigate the possible impact of Ambra1 on antitumor immunity and response to immunotherapy.

METHODS

This study was performed using an Ambra1-depleted Braf^V600E /Pten^{-/ -} genetically engineered mouse (GEM) model of melanoma, as well as GEM-derived allografts of Braf^V600E /Pten^{-/ -} and Braf^V600E /Pten^{-/ -}/Cdkn2a^{-/ -} tumors with Ambra1 knockdown. The effects of Ambra1 loss on the tumor immune microenvironment (TIME) were analyzed using NanoString technology, multiplex immunohistochemistry, and flow cytometry. Transcriptome and CIBERSORT digital cytometry analyses of murine melanoma samples and human melanoma patients (The Cancer Genome Atlas) were applied to determine the immune cell populations in null or low-expressing AMBRA1 melanoma. The contribution of Ambra1 on T-cell migration was evaluated using a cytokine array and flow cytometry. Tumor growth kinetics and overall survival analysis in Braf^V600E /Pten^{-/ -}/Cdkn2a^{-/ -} mice with Ambra1 knockdown were evaluated prior to and after administration of a programmed cell death protein-1 (PD-1) inhibitor.

RESULTS

Loss of Ambra1 was associated with altered expression of a wide range of cytokines and chemokines as well as decreased infiltration of tumors by regulatory T cells, a subpopulation of T cells with potent immune-suppressive properties. These changes in TIME composition were associated with the autophagic function of Ambra1. In the Braf^V600E /Pten^{-/ -}/Cdkn2a^{-/ -} model inherently resistant to immune checkpoint blockade, knockdown of Ambra1 led to accelerated tumor growth and reduced overall survival, but at the same time conferred sensitivity to anti-PD-1 treatment.

CONCLUSIONS

This study shows that loss of Ambra1 affects the TIME and the antitumor immune response in melanoma, highlighting new functions of Ambra1 in the regulation of melanoma biology.