Metabolic modelling-based in silico drug target prediction identifies six novel repurposable drugs for melanoma.

Despite high initial response rates to targeted kinase inhibitors, the majority of patients suffering from metastatic melanoma present with high relapse rates, demanding for alternative therapeutic options. We have previously developed a drug repurposing workflow to identify metabolic drug targets that, if depleted, inhibit the growth of cancer cells without harming healthy tissues. In the current study, we have applied a refined version of the workflow to specifically predict both, common essential genes across various cancer types, and melanoma-specific essential genes that could potentially be used as drug targets for melanoma treatment. The in silico single gene deletion step was adapted to simulate the knock-out of all targets of a drug on an objective function such as growth or energy balance. Based on publicly available, and in-house, large-scale transcriptomic data metabolic models for melanoma were reconstructed enabling the prediction of 28 candidate drugs and estimating their respective efficacy. Twelve highly efficacious drugs with low half-maximal inhibitory concentration values for the treatment of other cancers, which are not yet approved for melanoma treatment, were used for in vitro validation using melanoma cell lines. Combination of the top 4 out of 6 promising candidate drugs with BRAF or MEK inhibitors, partially showed synergistic growth inhibition compared to individual BRAF/MEK inhibition. Hence, the repurposing of drugs may enable an increase in therapeutic options e.g., for non-responders or upon acquired resistance to conventional melanoma treatments.

A multi-omics integrative approach unravels novel genes and pathways associated with senescence escape after targeted therapy in NRAS mutant melanoma.

Therapy Induced Senescence (TIS) leads to sustained growth arrest of cancer cells. The associated cytostasis has been shown to be reversible and cells escaping senescence further enhance the aggressiveness of cancers. Chemicals specifically targeting senescent cells, so-called senolytics, constitute a promising avenue for improved cancer treatment in combination with targeted therapies. Understanding how cancer cells evade senescence is needed to optimise the clinical benefits of this therapeutic approach. Here we characterised the response of three different NRAS mutant melanoma cell lines to a combination of CDK4/6 and MEK inhibitors over 33 days. Transcriptomic data show that all cell lines trigger a senescence programme coupled with strong induction of interferons. Kinome profiling revealed the activation of Receptor Tyrosine Kinases (RTKs) and enriched downstream signaling of neurotrophin, ErbB and insulin pathways. Characterisation of the miRNA interactome associates miR-211-5p with resistant phenotypes. Finally, iCell-based integration of bulk and single-cell RNA-seq data identifies biological processes perturbed during senescence and predicts 90 new genes involved in its escape. Overall, our data associate insulin signaling with persistence of a senescent phenotype and suggest a new role for interferon gamma in senescence escape through the induction of EMT and the activation of ERK5 signaling.

Single-cell transcriptomics of NRAS-mutated melanoma transitioning to drug resistance reveals P2RX7 as an indicator of early drug response.

Treatment options for patients with NRAS-mutant melanoma are limited and lack an efficient targeted drug combination that significantly increases overall and progression-free survival. In addition, targeted therapy success is hampered by the inevitable emergence of drug resistance. A thorough understanding of the molecular processes driving cancer cells' escape mechanisms is crucial to tailor more efficient follow-up therapies. We performed single-cell RNA sequencing of NRAS-mutant melanoma treated with MEK1/2 plus CDK4/6 inhibitors to decipher transcriptional transitions during the development of drug resistance. Cell lines resuming full proliferation (FACs [fast-adapting cells]) and cells that became senescent (SACs [slow-adapting cells]) over prolonged treatment were identified. The early drug response was characterized by transitional states involving increased ion signaling, driven by upregulation of the ATP-gated ion channel P2RX7. P2RX7 activation was associated with improved therapy responses and, in combination with targeted drugs, could contribute to the delayed onset of acquired resistance in NRAS-mutant melanoma.

Cross-Linking Ligation and Sequencing of Hybrids (qCLASH) Reveals an Unpredicted miRNA Targetome in Melanoma Cells.

MicroRNAs are key post-transcriptional gene regulators often displaying aberrant expression patterns in cancer. As microRNAs are promising disease-associated biomarkers and modulators of responsiveness to anti-cancer therapies, a solid understanding of their targetome is crucial. Despite enormous research efforts, the success rates of available tools to reliably predict microRNAs (miRNA)-target interactions remains limited. To investigate the disease-associated miRNA targetome, we have applied modified cross-linking ligation and sequencing of hybrids (qCLASH) to BRAF-mutant melanoma cells. The resulting RNA-RNA hybrid molecules provide a comprehensive and unbiased snapshot of direct miRNA-target interactions. The regulatory effects on selected miRNA target genes in predicted vs. non-predicted binding regions was validated by miRNA mimic experiments. Most miRNA-target interactions deviate from the central dogma of miRNA targeting up to 60% interactions occur via non-canonical seed pairing with a strong contribution of the 3' miRNA sequence, and over 50% display a clear bias towards the coding sequence of mRNAs. miRNAs targeting the coding sequence can directly reduce gene expression (miR-34a/CD68), while the majority of non-canonical miRNA interactions appear to have roles beyond target gene suppression (miR-100/AXL). Additionally, non-mRNA targets of miRNAs (lncRNAs) whose interactions mainly occur via non-canonical binding were identified in melanoma. This first application of CLASH sequencing to cancer cells identified over 8 K distinct miRNA-target interactions in melanoma cells. Our data highlight the importance non-canonical interactions, revealing further layers of complexity of post-transcriptional gene regulation in melanoma, thus expanding the pool of miRNA-target interactions, which have so far been omitted in the cancer field.

Systematic Transcriptional Profiling of Responses to STAT1- and STAT3-Activating Cytokines in Different Cancer Types.

Cytokines orchestrate responses to pathogens and in inflammatory processes, but they also play an important role in cancer by shaping the expression levels of cytokine response genes. Here, we conducted a large profiling study comparing miRNome and mRNA transcriptome data generated following different cytokine stimulations. Transcriptomic responses to STAT1- (IFNγ, IL-27) and STAT3-activating cytokines (IL6, OSM) were systematically compared in nine cancerous and non-neoplastic cell lines of different tissue origins (skin, liver and colon). The largest variation in our datasets was seen between cell lines of the three different tissues rather than stimuli. Notably, the variability in miRNome datasets was a lot more pronounced than in mRNA data. Our data also revealed that cells of skin, liver and colon tissues respond very differently to cytokines and that the cell signaling networks activated or silenced in response to STAT1- or STAT3-activating cytokines are specific to the tissue and the type of cytokine. However, globally, STAT1-activating cytokines had stronger effects than STAT3-inducing cytokines with most significant responses in liver cells, showing more genes upregulated and with higher fold change. A more detailed analysis of gene regulations upon cytokine stimulation in these cells provided insights into STAT1- versus STAT3-driven processes in hepatocarcinogenesis. Finally, independent component analysis revealed interconnected transcriptional networks distinct between cancer cells and their healthy counterparts.

Hypoxia-Induced Adaptations of miRNomes and Proteomes in Melanoma Cells and Their Secreted Extracellular Vesicles.

Reduced levels of intratumoural oxygen are associated with hypoxia-induced pro-oncogenic events such as invasion, metabolic reprogramming, epithelial-mesenchymal transition, metastasis and resistance to therapy, all favouring cancer progression. Small extracellular vesicles (EV) shuttle various cargos (proteins, miRNAs, DNA and others). Tumour-derived EVs can be taken up by neighbouring or distant cells in the tumour microenvironment, thus facilitating intercellular communication. The quantity of extracellular vesicle secretion and their composition can vary with changing microenvironmental conditions and disease states. Here, we investigated in melanoma cells the influence of hypoxia on the content and number of secreted EVs. Whole miRNome and proteome profiling revealed distinct expression patterns in normoxic or hypoxic growth conditions. Apart from the well-known miR-210, we identified miR-1290 as a novel hypoxia-associated microRNA, which was highly abundant in hypoxic EVs. On the other hand, miR-23a-5p and -23b-5p were consistently downregulated in hypoxic conditions, while the protein levels of the miR-23a/b-5p-predicted target IPO11 were concomitantly upregulated. Furthermore, hypoxic melanoma EVs exhibit a signature consisting of six proteins (AKR7A2, DDX39B, EIF3C, FARSA, PRMT5, VARS), which were significantly associated with a poor prognosis for melanoma patients, indicating that proteins and/or miRNAs secreted by cancer cells may be exploited as biomarkers.

Kinase inhibitor library screening identifies synergistic drug combinations effective in sensitive and resistant melanoma cells.

BACKGROUND: Melanoma is the most aggressive and deadly form of skin cancer with increasing case numbers worldwide. The development of inhibitors targeting mutated BRAF (found in around 60% of melanoma patients) has markedly improved overall survival of patients with late-stage tumors, even more so when combined with MEK inhibitors targeting the same signaling pathway. However, invariably patients become resistant to this targeted therapy resulting in rapid progression with treatment-refractory disease. The purpose of this study was the identification of new kinase inhibitors that do not lead to the development of resistance in combination with BRAF inhibitors (BRAFi), or that could be of clinical benefit as a 2nd line treatment for late-stage melanoma patients that have already developed resistance. METHODS: We have screened a 274-compound kinase inhibitor library in 3 BRAF mutant melanoma cell lines (each one sensitive or made resistant to 2 distinct BRAFi). The screening results were validated by dose-response studies and confirmed the killing efficacies of many kinase inhibitors. Two different tools were applied to investigate and quantify potential synergistic effects of drug combinations: the Chou-Talalay method and the Synergyfinder application. In order to exclude that resistance to the new treatments might occur at later time points, synergistic combinations were administered to fluorescently labelled parental and resistant cells over a period of > 10 weeks. RESULTS: Eight inhibitors targeting Wee1, Checkpoint kinase 1/2, Aurora kinase, MEK, Polo-like kinase, PI3K and Focal adhesion kinase killed melanoma cells synergistically when combined with a BRAFi. Additionally, combination of a Wee1 and Chk inhibitor showed synergistic killing effects not only on sensitive cell lines, but also on intrinsically BRAFi- and treatment induced-resistant melanoma cells. First in vivo studies confirmed these observations. Interestingly, continuous treatment with several of these drugs, alone or in combination, did not lead to emergence of resistance. CONCLUSIONS: Here, we have identified new, previously unexplored (in the framework of BRAFi resistance) inhibitors that have an effect not only on sensitive but also on BRAFi-resistant cells. These promising combinations together with the new immunotherapies could be an important step towards improved 1st and 2nd line treatments for late-stage melanoma patients.

A new ALK isoform transported by extracellular vesicles confers drug resistance to melanoma cells.

BACKGROUND: Drug resistance remains an unsolved clinical issue in oncology. Despite promising initial responses obtained with BRAF and MEK kinase inhibitors, resistance to treatment develops within months in virtually all melanoma patients. METHODS: Microarray analyses were performed in BRAF inhibitor-sensitive and resistant cell lines to identify changes in the transcriptome that might play a role in resistance. siRNA approaches and kinase inhibitors were used to assess the involvement of the identified Anaplastic Lymphoma Kinase (ALK) in drug resistance. The capability of extracellular vesicles (EVs) to transfer drug resistant properties was investigated in co-culture assays. RESULTS: Here, we report a new mechanism of acquired drug resistance involving the activation of a novel truncated form of ALK. Knock down or inhibition of ALK re-sensitised resistant cells to BRAF inhibition and induced apoptosis. Interestingly, truncated ALK was also secreted into EVs and we show that EVs were the vehicle for transferring drug resistance. CONCLUSIONS: To our knowledge, this is the first report demonstrating the functional involvement of EVs in melanoma drug resistance by transporting a truncated but functional form of ALK, able to activate the MAPK signalling pathway in target cells. Combined inhibition of ALK and BRAF dramatically reduced tumour growth in vivo. These findings make ALK a promising clinical target in melanoma patients.

Impact of BRAF kinase inhibitors on the miRNomes and transcriptomes of melanoma cells.

BACKGROUND: Melanoma is an aggressive skin cancer with increasing incidence worldwide. The development of BRAF kinase inhibitors as targeted treatments for patients with BRAF-mutant tumours contributed profoundly to an improved overall survival of patients with metastatic melanoma. Despite these promising results, the emergence of rapid resistance to targeted therapy remains a serious clinical issue. METHODS: To investigate the impact of BRAF inhibitors on miRNomes and transcriptomes, we used in vitro melanoma models consisting of BRAF inhibitor-sensitive and -resistant cell lines generated in our laboratory. Subsequently, microarray analyses were performed followed by RT-qPCR validations. RESULTS: Regarding miRNome and transcriptome changes, the long-term effects of BRAF inhibition differed in a cell line-specific manner with the two different BRAF inhibitors inducing comparable responses in three melanoma cell lines. Despite this heterogeneity, several miRNAs (e.g. miR-92a-1-5p, miR-708-5p) and genes (e.g. DOK5, PCSK2) were distinctly differentially expressed in drug-resistant versus -sensitive cell lines. Analyses of coexpressed miRNAs, as well as inversely correlated miRNA-mRNA pairs, revealed a low MITF/AXL ratio in two drug-resistant cell lines that might be regulated by miRNAs. CONCLUSION: Several genes and miRNAs were differentially regulated in the drug-resistant and -sensitive cell lines and might be considered as prognostic and/or diagnostic resistance biomarkers in melanoma drug resistance. GENERAL SIGNIFICANCE: Thus far, only little information is available on the significance and role of miRNAs with respect to kinase inhibitor treatments and emergence of drug resistance. In this study, promising miRNAs and genes were identified and associated to BRAF inhibitor-mediated resistance in melanoma. This article is part of a Special Issue entitled "Biochemistry of Synthetic Biology - Recent Developments" Guest Editor: Dr. Ilka Heinemann and Dr. Patrick O'Donoghue.

Comparison of a healthy miRNome with melanoma patient miRNomes: are microRNAs suitable serum biomarkers for cancer?

MiRNAs are increasingly recognized as biomarkers for the diagnosis of cancers where they are profiled from tumor tissue (intracellular miRNAs) or serum/plasma samples (extracellular miRNAs). To improve detection of reliable biomarkers from blood samples, we first compiled a healthy reference miRNome and established a well-controlled analysis pipeline allowing for standardized quantification of circulating miRNAs. Using whole miRNome and custom qPCR arrays, miRNA expression profiles were analyzed in 126 serum, whole blood and tissue samples of healthy volunteers and melanoma patients and in primary melanocyte and keratinocyte cell lines. We found characteristic signatures with excellent prognostic scores only in late stage but not in early stage melanoma patients. Upon comparison of melanoma tissue miRNomes with matching serum samples, several miRNAs were identified to be exclusively tissue-derived (miR-30b-5p, miR-374a-5p and others) while others had higher expression levels in serum (miR-3201 and miR-122-5p). Here we have compiled a healthy and widely applicable miRNome from serum samples and we provide strong evidence that levels of cell-free miRNAs only change significantly at later stages of melanoma progression, which has serious implications for miRNA biomarker studies in cancer.

New target genes of MITF-induced microRNA-211 contribute to melanoma cell invasion.

The non-coding microRNAs (miRNA) have tissue- and disease-specific expression patterns. They down-regulate target mRNAs, which likely impacts on most fundamental cellular processes. Differential expression patterns of miRNAs are currently being exploited for identification of biomarkers for early disease diagnosis, prediction of progression for melanoma and other cancers and as promising drug targets, since they can easily be inhibited or replaced in a given cellular context. Before successfully manipulating miRNAs in clinical settings, their precise expression levels, endogenous functions and thus their target genes have to be determined. MiR-211, a melanocyte lineage-specific small non-coding miRNA, is located in an intron of TRPM1, a target gene of the microphtalmia-associated transcription factor (MITF). By transcriptionally up-regulating TRPM1, MITF, which is critical for both melanocyte differentiation and survival and for melanoma progression, indirectly drives the expression of miR-211. Expression of this miRNA is often reduced in melanoma samples. Here, we investigated functional roles of miR-211 by identifying and studying new target genes. We show that MITF-correlated miR-211 expression levels are mostly but not always reduced in a panel of 11 melanoma cell lines and in primary and metastatic melanoma compared to normal melanocytes and nevi, respectively. MiR-211 itself only marginally impacted on cell invasion and migration, while perturbation of some new miR-211 target genes, such as AP1S2, SOX11, IGFBP5, and SERINC3 significantly increased invasion. These results and the variable expression levels of miR-211 raise serious doubts on the value of miR-211 as a melanoma tumor-suppressing miRNA and/or as a biomarker for melanoma.

Interplay of microRNAs, transcription factors and target genes: linking dynamic expression changes to function.

MicroRNAs (miRNAs) are ubiquitously expressed small non-coding RNAs that, in most cases, negatively regulate gene expression at the post-transcriptional level. miRNAs are involved in fine-tuning fundamental cellular processes such as proliferation, cell death and cell cycle control and are believed to confer robustness to biological responses. Here, we investigated simultaneously the transcriptional changes of miRNA and mRNA expression levels over time after activation of the Janus kinase/Signal transducer and activator of transcription (Jak/STAT) pathway by interferon-γ stimulation of melanoma cells. To examine global miRNA and mRNA expression patterns, time-series microarray data were analysed. We observed delayed responses of miRNAs (after 24-48 h) with respect to mRNAs (12-24 h) and identified biological functions involved at each step of the cellular response. Inference of the upstream regulators allowed for identification of transcriptional regulators involved in cellular reactions to interferon-γ stimulation. Linking expression profiles of transcriptional regulators and miRNAs with their annotated functions, we demonstrate the dynamic interplay of miRNAs and upstream regulators with biological functions. Finally, our data revealed network motifs in the form of feed-forward loops involving transcriptional regulators, mRNAs and miRNAs. Additional information obtained from integrating time-series mRNA and miRNA data may represent an important step towards understanding the regulatory principles of gene expression.

Interferon-γ-induced activation of Signal Transducer and Activator of Transcription 1 (STAT1) up-regulates the tumor suppressing microRNA-29 family in melanoma cells.

BACKGROUND: The type-II-cytokine IFN-γ is a pivotal player in innate immune responses but also assumes functions in controlling tumor cell growth by orchestrating cellular responses against neoplastic cells. The role of IFN-γ in melanoma is not fully understood: it is a well-known growth inhibitor of melanoma cells in vitro. On the other hand, IFN-γ may also facilitate melanoma progression. While interferon-regulated genes encoding proteins have been intensively studied since decades, the contribution of miRNAs to effects mediated by interferons is an emerging area of research.We recently described a distinct and dynamic regulation of a whole panel of microRNAs (miRNAs) after IFN-γ-stimulation. The aim of this study was to analyze the transcriptional regulation of miR-29 family members in detail, identify potential interesting target genes and thus further elucidate a potential signaling pathway IFN-γ → Jak→ P-STAT1 → miR-29 → miR-29 target genes and its implication for melanoma growth. RESULTS: Here we show that IFN-γ induces STAT1-dependently a profound up-regulation of the miR-29 primary cluster pri-29a~b-1 in melanoma cell lines. Furthermore, expression levels of pri-29a~b-1 and mature miR-29a and miR-29b were elevated while the pri-29b-2~c cluster was almost undetectable. We observed an inverse correlation between miR-29a/b expression and the proliferation rate of various melanoma cell lines. This finding could be corroborated in cells transfected with either miR-29 mimics or inhibitors. The IFN-γ-induced G1-arrest of melanoma cells involves down-regulation of CDK6, which we proved to be a direct target of miR-29 in these cells. Compared to nevi and normal skin, and metastatic melanoma samples, miR-29a and miR-29b levels were found strikingly elevated in certain patient samples derived from primary melanoma. CONCLUSIONS: Our findings reveal that the miR-29a/b1 cluster is to be included in the group of IFN- and STAT-regulated genes. The up-regulated miR-29 family members may act as effectors of cytokine signalling in melanoma and other cancer cells as well as in the immune system.

Dynamic regulation of microRNA expression following interferon-γ-induced gene transcription.

MicroRNAs are major players in post-transcriptional gene regulation. Even small changes in miRNA levels may have profound consequences for the expression levels of target genes. Hence, miRNAs themselves need to be tightly, albeit dynamically, regulated. Here, we investigated the dynamic behavior of miRNAs over a wide time range following stimulation of melanoma cells with interferon-γ (IFN-γ), which activates the transcription factor STAT1. By applying several bioinformatic and statistical software tools for visualization and identification of differentially expressed miRNAs derived from time-series microarray experiments, 8.9% of 1105 miRNAs appeared to be directly or indirectly regulated by STAT1. Focusing on distinct dynamic expression patterns, we found that the majority of robust miRNA expression changes occurred in the intermediate time range (24-48 h). Three miRNAs (miR-27a, miR-30a, miR-34a) had a delayed regulation occurring at 72 h while none showed significant expression changes at early time points between 30 min and 6 h. Expression patterns of individual miRNAs were altered gradually over time or abruptly increased or decreased between two time points. Furthermore, we observed coordinated dynamic transcription of most miRNA clusters while few were found to be regulated independently of their genetic cluster. Most interestingly, several "star" or passenger strand sequences were specifically regulated over time while their "guide" strands were not.

Signatures of microRNAs and selected microRNA target genes in human melanoma.

Small noncoding microRNAs (miRNA) regulate the expression of target mRNAs by repressing their translation or orchestrating their sequence-specific degradation. In this study, we investigated miRNA and miRNA target gene expression patterns in melanoma to identify candidate biomarkers for early and progressive disease. Because data presently available on miRNA expression in melanoma are inconsistent thus far, we applied several different miRNA detection and profiling techniques on a panel of 10 cell lines and 20 patient samples representing nevi and primary or metastatic melanoma. Expression of selected miRNAs was inconsistent when comparing cell line-derived and patient-derived data. Moreover, as expected, some discrepancies were also detected when miRNA microarray data were correlated with qPCR-measured expression levels. Nevertheless, we identified miRNA-200c to be consistently downregulated in melanocytes, melanoma cell lines, and patient samples, whereas miRNA-205 and miRNA-23b were markedly reduced only in patient samples. In contrast, miR-146a and miR-155 were upregulated in all analyzed patients but none of the cell lines. Whole-genome microarrays were performed for analysis of selected melanoma cell lines to identify potential transcriptionally regulated miRNA target genes. Using Ingenuity pathway analysis, we identified a deregulated gene network centered around microphthalmia-associated transcription factor, a transcription factor known to play a key role in melanoma development. Our findings define miRNAs and miRNA target genes that offer candidate biomarkers in human melanoma.

Recombinant interleukin-24 lacks apoptosis-inducing properties in melanoma cells.

IL-24, also known as melanoma differentiation antigen 7 (mda-7), is a member of the IL-10 family of cytokines and is mainly produced by Th(2) cells as well as by activated monocytes. Binding of IL-24 to either of its two possible heterodimeric receptors IL-20R1/IL-20R2 and IL-22R/IL-20R2 activates STAT3 and/or STAT1 in target tissues such as lung, testis, ovary, keratinocytes and skin. To date, the physiological properties of IL-24 are still not well understood but available data suggest that IL-24 affects epidermal functions by increasing proliferation of dermal cells. In stark contrast to its "normal" and physiological behaviour, IL-24 has been reported to selectively and efficiently kill a vast variety of cancer cells, especially melanoma cells, independent of receptor expression and Jak-STAT signalling. These intriguing properties have led to the development of adenovirally-expressed IL-24, which is currently being evaluated in clinical trials. Using three different methods, we have analysed a large panel of melanoma cell lines with respect to IL-24 and IL-24 receptor expression and found that none of the investigated cell lines expressed sufficient amounts of functional receptor pairs and therefore did not react to IL-24 stimulation with Jak/STAT activation. Results for three cell lines contrasted with previous studies, which reported presence of IL-24 receptors and activation of STAT3 following IL-24 stimulation. Furthermore, evaluating four different sources and modes of IL-24 administration (commercial recombinant IL-24, bacterially expressed GST-IL-24 fusion protein, IL-24 produced from transfected Hek cells, transiently over-expressed IL-24) no induction or increase in cell death was detected when compared to appropriate control treatments. Thus, we conclude that the cytokine IL-24 itself has no cancer-specific apoptosis-inducing properties in melanoma cells.