Next Article in Journal
Capsaicin Targets tNOX (ENOX2) to Inhibit G1 Cyclin/CDK Complex, as Assessed by the Cellular Thermal Shift Assay (CETSA)
Previous Article in Journal
Is NO the Answer? The Nitric Oxide Pathway Can Support Bone Morphogenetic Protein 2 Mediated Signaling
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 8
Issue 10
10.3390/cells8101274
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Oncogenic Signaling Pathways in BRAF-Mutant Melanoma Cells Are Modulated by Naphthalene Diimide-Like G-Quadruplex Ligands

by
Marta Recagni
1,†,
Martina Tassinari
2,†,
Filippo Doria
3,
Graziella Cimino-Reale
1,
Nadia Zaffaroni
1,
Mauro Freccero
3,
Marco Folini
1,* and
Sara N. Richter
2,*
1
Department of Applied Research and Technological Development, Fondazione IRCCS Istituto Nazionale dei Tumori di Milano, Via G.A. Amadeo 42, 20133 Milan, Italy
2
Department of Molecular Medicine, University of Padua, via A. Gabelli 63, 35121 Padua, Italy
3
Department of Chemistry, University of Pavia, v. le Taramelli 10, 27100 Pavia, Italy
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this paper.
Cells 2019, 8(10), 1274; https://doi.org/10.3390/cells8101274
Submission received: 25 July 2019 / Revised: 27 September 2019 / Accepted: 16 October 2019 / Published: 18 October 2019
Browse Figures
Review Reports Versions Notes

Abstract

:
Melanoma is the most aggressive and deadly type of skin cancer. Despite the advent of targeted therapies directed against specific oncogene mutations, melanoma remains a tumor that is very difficult to treat, and ultimately remains incurable. In the past two decades, stabilization of the non-canonical nucleic acid G-quadruplex structures within oncogene promoters has stood out as a promising approach to interfere with oncogenic signaling pathways in cancer cells, paving the way toward the development of G-quadruplex ligands as antitumor drugs. Here, we present the synthesis and screening of a library of differently functionalized core-extended naphthalene diimides for their activity against the BRAFV600E-mutant melanoma cell line. The most promising compound was able to stabilize G-quadruplexes that formed in the promoter regions of two target genes relevant to melanoma, KIT and BCL-2. This activity led to the suppression of protein expression and thus to interference with oncogenic signaling pathways involved in BRAF-mutant melanoma cell survival, apoptosis, and resistance to drugs. This G-quadruplex ligand thus represents a suitable candidate for the development of melanoma treatment options based on a new mechanism of action and could reveal particular significance in the context of resistance to targeted therapies of BRAF-mutant melanoma cells.

Graphical Abstract

1. Introduction

Melanoma, a malignancy that develops from melanin-producing cells, is the most aggressive and deadly type of skin cancer. Though it is a rare tumor, accounting for only 4% of all skin cancers, it is responsible for the highest number of skin cancer-related deaths worldwide and it is the most difficult cancer to treat [1,2]. Chemotherapy is the standard curative option for advanced stage melanoma patients, even though it has an overall low efficacy and limited response rates [1].
The characterization of the melanoma mutational landscape has led to the identification of specific gene/pathway alterations that have had remarkable implications on the management of this disease [3]. In this context, a major breakthrough in the understanding of melanoma biology was the identification of the BRAFV600E mutation [3], which contributed to the development of anti-oncogenic compounds targeting activated mutant BRAF. Indeed, up to 50% of melanomas derived from the skin without chronic sun damage bear mutations in the BRAF gene, with BRAFV600E representing up to 90% of the activating mutations [4]. This scenario led to the development of the ATP-competitive BRAFV600E selective inhibitor, vemurafenib [5]. This drug has shown dramatic success in clinical trials, as it promotes tumor regression and increases the overall survival of metastatic melanoma patients [6]. Unfortunately, resistance mechanisms that render melanoma unresponsive to other therapeutic strategies develop within a few months from the beginning of treatment [6]. The second most common mutation in melanoma affects the NRAS gene in codon 61, producing such mutations as Q61R or Q61K [3]. These occur in 20–30% of melanoma patients and are mutually exclusive with BRAF mutations [1], except in resistant melanomas after targeted therapy, which may harbor co-occurring BRAF and NRAS mutations [3].
Recent evidence has indicated that the transmembrane receptor tyrosine kinase c-KIT may also be an attractive therapeutic target in melanoma [7]. Genetic alterations of KIT in melanoma include somatic gain-of-function mutations and copy number increases of wild-type KIT [7], whereas mutant receptors were found only in 2% of all cutaneous melanomas, thus representing a rare event for targeted treatment, and in up to 20% of mucosal, acral, and chronic sun-damaged skin melanomas [8]. A series of different KIT mutations, among which L576P was detected in one-third of all cases, was found, although many of them are not suitable targets [4].
The identification of druggable mutation-specific oncogene targets significantly contributed to the expansion of the arsenal of available therapies for patients with advanced melanoma over the past few years. The introduction of targeted therapies, such as BRAF (vemurafenib and dabrafenib) and MEK (trametinib and cobimetinib) inhibitors, as single agents or in combination [1,2], led to both improved response rates and mean overall survival of metastatic melanoma patients bearing the BRAFV600E mutation or mutant NRAS [3,8]. On the other hand, mutant c-KIT may be able to be targeted by tyrosine kinase receptor inhibitors (e.g., imatinib, sunitinib, and dasatinib), although, at present, clinical benefits have been reported only for imatinib in melanoma patients with KIT point mutations in exon 11 or 13, and not in those harboring gene amplification [8].
As per any oncogene-targeted therapy, treatment failure is associated with mechanisms of acquired drug resistance, which may rely on the reactivation of MAPK signaling, the activation of substitutive oncogenic pathways, such as that mediated by PI3K/AKT, as well as on the over-activation of growth factor receptors and the capability to evade apoptosis [1,8,9]. In this context, the deregulation of the BCL-2 family of proteins plays a major role in the evasion of melanoma cell apoptosis in response to treatment [9]. Notably, many BCL-2 proteins are downstream factors of the RAS/BRAF/MAPK and PI3K/AKT signaling pathways, the activation of which contributes to the relapse of melanoma from treatment with targeted therapies [9]. Multiple mechanisms have been reported to be responsible for the deregulation of BCL-2 protein family [9]. The development of strategies to target these pro-survival factors in melanoma has been a central theme for years [10], and may represent an alternative option to defeat melanoma as well as to overcome resistance to current targeted therapies [9].
This scenario supports the rationale for drug combination approaches [2] or, alternatively, for the use of single multi-targeting drug molecules, which are arising as valuable alternative tools to therapeutic regimens based on drug combinations [11], in order to overcome drug resistance and hopefully obtain long-term responses.
Nucleic acids can fold into several structural motifs to assemble the functional structural conformation for their precise biological roles in specific cellular environments. In particular, guanine (G)-rich sequences can self-associate into stacks of G-quartets using Hoogsten-type hydrogen bonds to form complex secondary structures knows as G-quadruplexes (G4s) [12], which are stabilized by K+ cations under physiological conditions [13]. In recent years, G4s have attracted great attention, largely due to both their peculiar polymorphisms [14] and critical regulatory roles in biological processes [15], such as modulation of gene expression [16], regulation of epigenetic modifications [17], telomerase dysfunction [18], transcription [19], genomic instability [20], and histone modifications [21]. Their implication in the pathogenesis of cancer [22] and neurodegenerative diseases [23,24] was extensively described, providing new possible targets in a number of different pathologies. In vivo formation of G4s was consolidated by the discovery of cellular proteins that specifically process G4s [25,26] and the development of G4-specific antibodies [27,28] that allowed both direct visualization and chromatin immunoprecipitation (ChIP) sequencing-based approaches [29,30,31]. Besides telomeres [32], G4s are mainly present in gene promoters and in regions that are close to transcription start sites (TSS) [33]. Based on the observation that G4s are located in the core or proximal promoter of genes related to the six hallmarks of cancer [34], a large number of G4 ligands were designed to target oncogene promoters. The hypothesis that G4 induction and stabilization by small molecules could be exploited to repress gene transcription was first demonstrated with the porphyrin ligand TMPyP4, which was able to down-regulate MYC expression [35], paving the way for the quest of specific G4-stabilizing molecules. Since then, many G4 ligands have been developed [36,37], with some showing interesting biological activity and proceeding to clinical trials [38,39]. These agents belong to different chemical families. Among them, naphthalene diimides (NDIs) are a class of small molecules that recognize, induce, and stabilize G4s with high affinity. Thanks to its electron-deficient aromatic systems, the NDI core is prone to interact via π–π stacking with an electron-rich partner, like the G-quartet [40]. Moreover, the extension of the NDI core offers an opportunity to enhance both G4 binding potency and selectivity [41,42]. In the present investigation, we explored the possibility to interfere with oncogenic signaling pathways in BRAF-mutant melanoma cells through NDI-based G4 ligands. A small library of core-extended NDIs (c-exNDIs 1–14, Supporting Information Figure S1) was screened in terms of cytotoxicity on A-375 melanoma cells, and the activity of the most active compound (1) was further investigated. Combining biophysical and molecular methods, we proved that 1 greatly stabilized G4s within human oncogene promoters. At the cellular level, the exposure of BRAF-mutant melanoma cells to the compound resulted in a marked down-modulation of KIT, which was paralleled by the shutdown of the MAPK and PI3K/AKT signaling pathways, both at the transcriptomic and post-translational levels. Moreover, 1-treated cells were characterized by an expression profile of BCL-2 family members indicating a trend toward a pro-apoptotic environment. We thus present a suitable candidate for the development of novel treatment options based on G4 targeting, which may be particularly relevant in the context of resistance to targeted therapies of BRAF-mutant melanoma cells.

2. Materials and Methods

2.1. Materials and General Procedures

Reagents, solvents, and chemicals were purchased from TCI (Milan, Italy) or Sigma-Aldrich (Milan, Italy) and used as supplied, without further purification. All oligonucleotides (Supporting Information Table S1) were purchased from Sigma-Aldrich (Milan, Italy). The tested compounds were dissolved in dimethylsulfoxide (DMSO), stored at −20 °C, and diluted in culture medium at the appropriate working concentrations immediately before use. The eluent for all of the HPLC analyses and purifications was 0.1% trifluoroacetic acid in water and acetonitrile. HPLC analysis was performed using Agilent system SERIES 1260. The column was XSelect HSS C18 (2.5 μm) (50 × 4.6 mm) (Waters, Milan, Italy). We used two different methods for analytical and preparative purposes. Analytical method: flow—1.4 mL/min; gradient—95% aqueous for 2 min, gradually to 100% of acetonitrile over 14 min and then isocratic flow for 4 min. Preparative HPLC purifications were attained using Agilent system SERIES 1260 and equipped with a preparative PUMP, a diode array system, and an automatic fraction collector. The preparative column was XSelect CSH Prep Phenyl-Hexyl 5 μm (150 × 30 mm) (Waters, Milan, Italy) using a 30 mL/min flow. The eluent was water + 0.1% trifluoroacetic acid, gradually decreased from 95 to 85%, over 4 min, then gradually to 65% water in 16 min. 1H- and 13C-NMR spectra were recorded on a Bruker ADVANCE 300 MHz.

2.2. Cell Lines

Human tumor cell lines and primary dermal fibroblasts from adult skin were obtained as indicated in Supporting Information Table S2. Cells were maintained as a monolayer in the logarithmic growth phase in the appropriate growth medium in the presence of 10% fetal calf serum at 37 °C in a humidified incubator at 5% CO2. Cells were periodically monitored for DNA profiling by short tandem repeats analysis using the AmpFISTR Identifiler PCR amplification kit (Thermo Fisher Scientific, Monza, Italy).

2.3. Cell Viability Assay

The MTT (3-(4,5-dimethylthiazol 2-yl)-2,5-diphenyltetrazolium bromide, Sigma-Aldrich, Milan, Italy) assay was performed to assess the cytotoxicity of c-exNDIs on a human malignant melanoma A375 cell line. A375 cells (2.5 × 103) were plated in 96-well plates and incubated for 24 h. Cells were treated with increasing concentrations (1.95–500.00 nM) of tested compounds and incubated for an additional 48 h. Cell viability was evaluated by MTT assay by adding 10 μL of freshly dissolved solution of 5 mg/mL MTT in phosphate buffered saline (PBS) to each well. After 4 h of incubation, MTT crystals were solubilized in solubilization solution (10% sodium dodecyl sulfate (SDS) and 0.01 M HCl). After overnight incubation at 37 °C, absorbance was read at 570 nm. The percentage of cell viability was calculated as follows: Cell survival = (Awell − Ablank)/(Acontrol − Ablank) × 100, where blank denotes the medium without cells. The IC50 was defined as the concentration of compound required to inhibit cell grown by 50%. Data were expressed as mean values ± standard deviation of three independent experiments.

2.4. Gene Expression Analysis

Total RNA was isolated using Qiagen RNeasy Mini kit (Qiagen S.r.l., Milan, Italy) and digested with 20 U RNase-free DNase I (Qiagen S.r.l., Milan, Italy), according to the manufacturer’s instructions. Five hundred nanograms of total RNA was randomly primed and reverse transcribed using the High Capacity cDNA Reverse Transcription kit (Applied Biosystems, Carlsbad, CA, USA). The expression levels of 92 genes associated with the molecular mechanisms of cancer were assessed by real-time RT-PCR using TaqMan Arrays (PN4418806; Applied Biosystems), according to the manufacturer’s instructions. Amplifications were run on 7900HT fast real-time PCR system. Data were analyzed by SDS 2.2.2 software and, if not otherwise specified, reported as relative quantity (RQ) with respect to untreated cells (calibrator sample), according to the 2−ΔΔCt method [43], where Ct represents the threshold cycle. The 18S housekeeping gene present in each array was used as a normalizer. Differentially expressed genes in treated vs. untreated cells were sorted based on p < 0.05 (Student’s t–test) and considered up- or down-regulated by setting a fold-change of 1.5 as the cut-off.

2.5. Circular Dichroism

For circular dichroism (CD) analysis, oligonucleotides (Supporting Information Table S1) were diluted to a final concentration (4 μM) in lithium cacodylate buffer (10 mM, pH 7.4) and KCl 10 mM. Samples were annealed by heating at 95 °C for 5 min and gradually cooled to room temperature and, where indicated, 1 was added at the final concentration of 16 μM. CD spectra were recorded using Chirascan-Plus (Applied Photophysisics, Leatherhead, United Kingdom) equipped with a Peltier temperature controller using a quartz cell of 5 mm optical path length, over a wavelength range of 230–320 nm. The reported spectrum of each sample represented the average of 2 scans at 20 °C and was baseline corrected for signal contributions due to the buffer. Observed ellipticities were converted to mean residue ellipticity (θ) = deg × cm2 × dmol−1 (molar ellipicity). For the determination of melting temperature (Tm), spectra were recorded over a temperature range of 20–90 °C, with temperature increase of 5 °C. Tm values were calculated according to Greenfield [44].

2.6. Fluorescence Resonance Energy Transfer

For fluorescence resonance energy transfer (FRET), 6-carboxyfluorescein (FAM) 5′-end-labeled and 6-carboxy-tetramethylrhodamine (TAMRA) 3′-end-labeled oligonucleotides (0.25 μM) (Supporting Information Table S1) were folded in lithium cacodylate buffer (10 mM, pH 7.4) and KCl 100 mM or 10 mM by heating at 95 °C for 5 min and gradually cooling to room temperature. Where indicated, 1 was added at a final concentration of 0.25 μM and, after stabilization at 4 °C, samples were processed in a Light Cycler (Roche, Milan, Italy). Oligonucleotide melting was monitored by measuring FAM emission in the 30–95 °C temperature range with a gradient of 1 °C/min. Melting profiles were normalized as previously described [45] and Tm was defined as the temperature corresponding to the 0.5 fraction of the normalized fluorescence.

2.7. Taq Polymerase Stop Assay

The DNA primer (Supporting Information Table S1) was 5′-end labeled with [γ-32P]ATP using T4 polynucleotide kinase (Thermo Scientific, Milan, Italy) at 37 °C for 30 min. The labeled primer (final concentration 72 nM) was annealed to the templates (final concentration 36 nM) (Supporting Information Table S1) in lithium cacodylate buffer (10 mM, pH 7.4) in the presence or absence of KCl 10 mM by heating at 95 °C for 5 min and gradually cooling to room temperature to allow both primer annealing and G4 folding. Where specified, samples were incubated with 1 (37.5–150.0 nM) and primer extension was performed with 2 U/reaction of AmpliTaq Gold DNA polymerase (Applied Biosystem, Carlsbad, CA, USA) at 56 °C for 30 min. Reactions were stopped by ethanol precipitation and primer extension products were separated on a 14% denaturing gel, and visualized with phosphorimaging (Typhoon FLA 9000, GE Healthcare, Milan, Italy). Markers were prepared based on Maxam and Gilbert sequencing by PCR reaction with 32P-labeled primer. PCR products were treated with formic acid for 5 min at 25 °C and then with piperidin for 30 min at 90 °C.

2.8. Western Immunoblotting

Total protein extracts were prepared according to standard methods. Forty micrograms of protein extract was fractioned by SDS-PAGE and transferred onto Hybond nitrocellulose filters (RPN 303D, GE Healthcare, Milan, Italy). Filters were blocked in PBS-Tween 20 and 5% skim milk and incubated overnight with the following primary antibodies: mouse monoclonal anti-BCL-2 (sc-509, Santa Cruz Biotechnology, Dallas, TX, USA), anti-KIT (#3308; Cell Signaling Technology, Danvers, MT, USA), and anti-MYC (ab32; Abcam, Cambridge, United Kingdom); rabbit polyclonal anti-AKT(S473) (sc-9271, Santa Cruz Biotechnology), anti-BRAF (ab33899, Abcam,), anti-ERK1/2 (T202/Y204) (#9101, Cell Signaling Technology), anti-γ-H2AX (ab11174, Abcam), anti-p21waf1 (ab7960, Abcam), and anti PARP-1 (#9542, Cell Signaling Technology). Mouse monoclonal anti-Vinculin (VCL, V9131, Sigma-Aldrich, Milan, Italy) or anti-β-Actin (Ab8226, Abcam) antibodies were used to ensure equal protein loading. The filters were then probed with secondary peroxidase-linked whole antibodies (GE Healthcare, Milan, Italy) and subjected to autoradiography using the Novex® Enhanced Chemoluminescent Horseradish Peroxidase detection system (Thermo Fisher Scientific, Monza, Italy). Films were scanned (ImageScanner III, GE Healthcare, Milan, Italy) and images were processed by Photoshop7.0.1 or analyzed using ImageJ 1.46r.

2.9. Analysis of Cell Cycle Phase Distribution

For the cell cycle phase distribution analysis, both adherent and floating cells were fixed in 70% EtOH and incubated at 4 °C for 30 min in staining solution containing 50 mg/mL of propidium iodide, 50 mg/mL of RNase, and 0.05% Nonidet-P40 in PBS. Samples were analyzed with a Fluorescence Activated Cell Sorting (FACS)-Calibur cytofluorimeter (Becton Dickinson, Milan, Italy). At least 30,000 events were read, and the histograms were analyzed using the CellQuest software according to the Modfit model (Becton Dickinson, Milan, Italy).

2.10. Statistical Analysis

If not otherwise specified, the data were reported as mean values ± s.d. from at least three independent experiments. Two-sided Student’s t-test was used to analyze any differences between samples. The association between the IC50 values and the cell doubling times in the cancer cell lines was assessed using the Spearman rank order correlation coefficient (GraphPad Software Inc., San Diego, CA, USA). p values of <0.05 were considered statistically significant.

3. Results

3.1. Chemistry

A family of fourteen c-exNDIs (Supporting Information Figure S1) characterized by a lateral extension of the NDI core was selected to be screened against melanoma cells. The compounds were designed to expand the π-stacking interaction with the G-quartet moiety in comparison to NDIs, thereby investigating whether binding and affinity to G4s were improved. Compounds 1–10, 13, and 14 were previously synthesized and characterized [41,42], while 11 and 12 were novel ligands and were prepared according to the synthetic protocol outlined in Supplementary Information Scheme S1. The screening and biophysical assays depicted below were performed using the hydrochloride amine salts (1–7 × 2 HCl, and 11–14 × 2 HCl) and the iodide quaternary ammonium salts 9 and 10.

3.2. Screening of the Anticancer Activity of c-exNDIs in Melanoma Cells

We first assessed the cytotoxic activity of the c-exNDI series of ligands via MTT assay in a representative in vitro model of mutant BRAFV600E melanoma. Specifically, 48 h of exposure of A375 cells to increasing concentrations of each compound resulted in a remarkable and dose-dependent inhibition of cell growth, with IC50 values in the nanomolar range for most compounds (Table 1). Among the tested c-exNDIs, the unsubstituted ligands at C4 (Y = H, see Supporting Information Figure S1 for numbering) were generally more active compared to the bromine derivatives (Y = Br). 1, 2, 3, and 11 displayed the highest cytotoxicity levels, with IC50 values in the low nM range. On the other hand, 9 and 10, which contained two quaternary ammonium moieties, did not display any cytotoxic effects up to the highest tested concentration, likely due to their inability to enter the cells, as already shown for NDIs with permanent positive charges [46]. Low activity was also detected with 5 and 6, which exhibited IC50 values higher than 500 nM, probably due to the anionic carboxylate moiety. In general, substituents at C11 had detrimental effects on cytotoxicity. The other c-exNDIs showed intermediate values of cytotoxicity.
Compound 1, which was the most active one in the series (IC50 of 8 ± 1 nM, Table 1), was selected and assessed against a panel of various cell lines representative of tumors of different histological origins. As reported in Table S2, the compound exhibited a pronounced cytotoxic activity on the tested tumor cell lines, whereas it did not remarkably cause cytotoxic effects on normal human primary skin fibroblasts at concentrations within the nanomolar range.
This evidence was in keeping with our previous findings showing that the compound exhibited a pronounced cytotoxic activity in a panel of metastatic castration-resistant prostate cancer cells, while sparing normal prostate epithelial cells [47] and thus contributing to the delineation of a therapeutic window for the compound.

3.3. C-exNDI 1 Induces a Global Down-Regulation of Gene Expression in BRAFV600E-Mutant Melanoma Cells

It was hypothesized that the effects at the transcriptomic level of small molecule-mediated stabilization of G4 located within gene transcriptional regulatory elements could be more pronounced in rapidly dividing cells than in slowly dividing cells. In this context, Spearman rank correlation analysis revealed the occurrence of a possible relationship between the cytotoxic activity of compound 1 and the doubling times for the 9 cancer cell lines listed in Supporting Information Table S2. Specifically, a statistically significant trend toward a strong, positive monotonic correlation between IC50 values and the cell doubling times was revealed (rs = 0.8201, n = 9, p = 0.0108). This observation indicated that G4 formation most likely occurs during DNA replication [15], as shown by the evidence that the number of G4 structures increased in human cells during the S-phase of the cell cycle when the duplex strands became separated at the replication forks, potentially causing single-stranded DNA to fold more easily into secondary structures [28], whereas it decreased following the exposure to aphidicolin, an inhibitor of DNA replication [28].
A comparative analysis of gene expression was carried out on a pair of malignant melanoma cell lines (A375 and SKMEL-2), which was characterized by both a distinct rate of growth in vitro and a pronouncedly different sensitivity to 1 (Figure 1A,B). Results showed that 48 h of exposure to an equitoxic amount (IC50) of the compound triggered G4 stabilization (Figure 1C) and induced a global down-regulation of gene expression in A375 compared to SKMEL-2 cells (Supporting Information Table S3). In particular, real-time RT-PCR assessment of the expression levels of 92 genes commonly associated with the molecular mechanisms of cancer revealed that ~61% of genes were differently expressed (fold-change > |1.5|) in 1-treated compared to untreated A375 cells (e.g., 53 and 3 genes were down- and up-modulated, respectively, Supporting Information Table S3). Conversely, only ~22% of genes were found to be differently expressed (10 down- and 10 up-regulated) in 1-treated SKMEL-2 (Supporting Information Table S3). Notably, an over-representation enrichment analysis carried out using the GEne SeT AnaLysis Toolkit (WebGestalt, http://webgestalt.org) revealed that the most relevant over-represented pathways (p < 0.01) defined by the list of differently expressed genes in 1-treated vs. untreated A375 cells with a fold-change >|2.0| comprised terms such as signaling by BRAF and RAF fusions, oncogenic MAPK signaling, and signal attenuation (Supporting Information Table S4), which are all relevant, targetable pathways in BRAF-mutant melanoma cells [3]. In more detail, real-time RT-PCR data showed a marked decrease in the expression levels of factors associated with RAS/RAF/MAPK and PI3K/AKT signaling pathways, as well as significant modulations in the expression levels of BCL-2 family members, indicating a trend toward a pro-apoptotic milieu, with down-regulation of the anti-apoptotic BCL2 and BCL2L1 transcripts and the up-regulation of pro-apoptotic factors, such as BAX and BCL2L11 (Figure 1D). On the contrary, although the paucity of differently expressed genes showing a fold-change >|2.0| did not allow us to retrieve significantly over-represented pathways (Supporting Information Table S4), it was clear that the exposure of SKMEL-2 cells to 1 resulted in a completely different gene expression landscape, without any remarkable modulation in the expression levels of factors belonging to RAS/RAF/MAPK except for the down-modulation of MAP3K5 (0.65 ± 0.03; Supporting Information Table S3) and/or PI3K/AKT pathways, and no significant perturbations in the levels of BCL-2 family members, except for a marked down-modulation of the BID transcript (0.47 ± 0.005; Supporting Information Table S3).
The different pattern of gene expression modulation observed in melanoma cell lines exposed to compound 1 was paralleled by a distinct response in terms of DNA damage induction. In particular, similarly to cells exposed to pyridostatin, which was used as a reference G4 ligand, the exposure of A375 cells to 1 resulted in a marked induction of DNA damage, as suggested by the prominent occurrence of DNA damage foci observed by immunofluorescence in cells probed with anti-γ-H2AX antibody (Supplementary Figure S2A). On the other hand, and differently from the pyridostatin treatment, 1-treated SKMEL-2 cells did not show remarkable accumulation of γ-H2AX foci (Supplementary Figure S2A).
Notably, it could not be excluded that the different observed scenarios between A375 and SKMEL-2 cells in terms of 1-dependent gene expression modulation represented a distinct cell response to the compound that reflected their opposite genetic background (TP53WT, BRAFV600E, and NRASWT in A375 cells, and TP53G245S, BRAFWT, and NRASQ61R in SKMEL-2 cells; see COSMIC, the Catalogue Of Somatic Mutations In Cancer, https://cancer.sanger.ac.uk/cell_lines). Nonetheless, it is worth nothing that exposure of BRAFWT/NRASMUT cells to 1 resulted in a peculiar gene expression pattern with over-represented pathways, including platelet degranulation, senescence-associated secretory phenotypes, and oxidative stress-induced senescence (Supporting Information Table S4). This scenario suggested that 1-treated SKMEL-2 cells underwent changes at the cellular/molecular level, likely resembling phenotype switching [1]. This phenomenon of cell plasticity, which includes, among others, epithelial-to-mesenchymal-like transition, switching from differentiated to de-differentiated, and slowing proliferating phenotypes, as well as metabolic rewiring (including increased Reactive Oxygen Species production), represented an adaptive process in response to treatment-induced insults by which melanoma cells may adapt to therapies [1].

3.4. Melanoma-Relevant Genes as Possible G4 Targets of c-exNDI 1

The finding that the exposure of BRAF-mutant melanoma cells to our G4 ligand resulted in modulation of the expression levels of genes belonging to signaling pathways that are relevant to melanoma cell survival and drug-resistance prompted us to investigate whether the compound interacted with specific G4 targets at the interplay of RAS/RAF/MAPK and PI3K/AKT, as well as the intrinsic apoptosis pathways. In particular, we focused on KIT, BRAF, and BCL-2, three relevant target genes in melanoma known to bear G4-forming sequences within their promoters [48]. Using biophysical and biomolecular assays, we comparatively assessed the ability of 1 to stabilize the G4-forming sequences within the promoter regions of the selected genes in comparison to MYC, which is known for having the highest G4 density among oncogenes [40].
Five different oligonucleotides corresponding to the G4-forming sequences within the promoter regions of MYC, BRAF, BCL-2, and KIT (KIT-1 and KIT-2) (Supporting Information Table S1) were assayed by circular dichroism (CD) to analyze both their conformation and their stability in the presence and absence of the compound. In-keeping with data reported in the literature and concerning the biophysical characterization of G4 motifs of the selected genes [49,50,51,52,53], we found that the CD spectra of the tested oligonucleotides were all distinctive of G4 topologies. In particular, the c-myc and c-kit2 sequences formed a parallel-type G4 structure showing a major positive peak at around 264 nm and a negative one at 242 nm. A prevalent parallel conformation was also observed in the bcl-2, c-kit1, and b-raf templates, the last two displaying an additional positive shoulder at around 290 nm, indicating a mixed or hybrid G4 conformation (Figure 2A). Upon addition of 1, c-myc, bcl-2, c-kit1, and c-kit2 maintained a prevalent parallel conformation, with the molar ellipticity at 264 nm of c-myc, c-kit1, and c-kit2 decreasing to different extents (Figure 2B). Additionally, a slight variation in the molar ellipticity signal was observed in bcl-2 and c-kit1 at 290 nm (an increase and a decrease, respectively). In contrast, 1 induced a structural transition from a parallel-type to a mixed parallel–antiparallel G4 on the b-raf template, as indicated by the presence of two positive peaks at 264 and 290 nm (Figure 2B).
Stability of the G4 structures was next evaluated by CD thermal unfolding (Supporting Information Figure S3); the melting temperatures (Tm) are reported in Table 2. For c-myc, bcl-2, c-kit1, and c-kit2, a single transition between 20 °C and 90 °C was appreciable at 264 nm, leading to discrete Tm values. A remarkable increase in stability with significant variation of the Tm (>90 °C) was recorded in the presence of the compound compared to the absence of the compound. Given that 1 induced a mixed G4 on the b-raf templates, Tm was calculated at both 264 nm and 290 nm; a notable increase in the Tm was recorded (Table 2) upon addition of 1, indicating effective stabilization of b-raf G4.
We next employed a dual-labeled system. The stabilizing activity of 1 was assessed by Fluorescence Resonance Energy Transfer (FRET) melting assay on c-kit2 (Supporting Information Table S1), which was chosen as the representative oncogene G4-forming sequence. C-kit2 displayed a Tm of 76.4 ± 0.3 °C and 65.7 ± 0.7 °C at 100 mM and 10 mM of K+, respectively. In both conditions, c-kit2 was greatly stabilized by 1 with Tm > 90 °C (Supporting Information Figure S4A,B). A double stranded (ds) oligonucleotide (Supporting Information Table S1) was also used as a non-G4-forming sequence control to check the selectivity of 1 for G4 vs. the duplex. This showed a Tm of 65.7 ± 0.7 °C and 61.2 ± 01 °C at 100 mM and 10 mM of K+, respectively, with the compound showing that it was basically devoid of a stabilizing effect on dsDNA with variations in Tm values of <2 °C (Supporting Information Figure S4C,D).
To corroborate the ability of 1 to stabilize the target G4s, a Taq polymerase stop assay was performed using bcl-2, c-myc, b-raf, and c-kit1+2 templates, the latter containing both c-kit1 and c-kit2 G4-forming motifs linked by a 16 nucleotide-long sequence (Supporting Information Table S1). All tested oligonucleotides were designed to contain additional flanking bases at both the 5′- and 3′-ends; in particular, an additional sequence at the 3′-end was used as a primer annealing region and 5-T linker regions were added to both the 3′-end, to separate the primer annealing region and the first G of the G4 portion, and the 5′-end. All G4-forming sequences, except for b-raf, stopped the polymerase at the first G-rich region encountered by the polymerase, even in the absence of K+ (Figure 2C, lane 1: bcl-2, c-myc, c-kit1+2). In the presence of K+ (Figure 2C, lane 2: b-raf, bcl-2, c-myc, c-kit1+2), stop sites corresponding to the first 3′G-tract were visible in all templates and were more intense than in the absence of K+, thus indicating that K+ stimulated G4 folding. Upon addition of increasing amounts (37.5–150.0 nM) of 1, the intensity of the stop bands greatly rose in all templates in a concentration-dependent manner (Figure 2C, lanes 3–5), accompanied by a considerable reduction in full-length amplicons, thus corroborating the effective stabilization of the G4s by the compound at nM concentrations. Furthermore, additional stop sites corresponding to other G-rich regions were observed in bcl-2, c-myc, and c-kit1+2 templates, possibly indicating the presence of multiple G4 structures. Quantification of the stop sites confirmed the concentration dependence of polymerase stalling of the compound (Figure 2D). No stop sites were observed with a control sequence devoid of G-tracts (Figure 2C, lanes 1–4: non-G4 cnt), indicating that the observed polymerase inhibition was G4-dependent. Overall, these data were in line with those obtained by CD analysis and confirmed the ability of the chosen sequences to be stabilized by 1.

3.5. The Exposure of BRAFV600E Melanoma Cells to c-exNDI 1 Results in KIT and BCL-2 Protein Down-Regulation, Inhibition of MAPK and PI3K/AKT Signaling Pathways, and Apoptosis Induction

Gene expression analyses revealed a reduction in the expression levels of BRAF, BCL-2, and MYC mRNAs in A375 cells upon 48 h of exposure to the ligand (Figure 1D and Figure 3A), whereas the levels of KIT mRNA remained undetermined both in treated and untreated cells, being undetectable even after 40 cycles of PCR amplification (Supporting Information Table S3). Nonetheless, the validation of qRT-PCR data by Western immunoblotting revealed a time-dependent reduction of both KIT and Bcl-2 protein amounts in 1-treated A375 cells with respect to untreated A375 cells (Figure 3B), whereas MYC and BRAF protein levels were not remarkably affected by the treatment (Figure 3B). On the other hand, KIT and MYC mRNA levels were mildly reduced (RQ = 0.67 ± 0.05 and 0.74 ± 0.03, respectively) in SKMEL-2 cells exposed for 48 h to an equitoxic amount of the compound (Figure 3A and Supporting Information Table S3).
The observed reduction in KIT protein levels was accompanied by perturbations of downstream signaling pathways. In particular, 72 h of exposure of A375 cells to the compound resulted in a mild reduction (p < 0.05) of phosphorylated MAPK3 (p44/ERK1T202)—although with no modification in the phosphorylation status of MAPK1 (p42/ERK2Y204)—as well as in a pronounced decrease (p < 0.01) in the phosphorylation of AKTS473 (Figure 3C). As already observed, 1-treated cells were characterized by the prominent accumulation of γ-H2AX (Figure 3C) and by a marked increase in the amounts of the cyclin-dependent kinase inhibitor p21waf1 (Figure 3C). Consistent with gene expression data (Supporting Information Table S3), a moderate increase in p21waf1 protein levels was also appreciable in SKMEL-2 cells exposed to an equitoxic amount of the compound (Supporting Information Figure S2B). Conversely, 1-treated SKMEL-2 cells were characterized by an increase in the phosphorylated forms of ERK1/2 and AKT, without showing a remarkable induction of DNA damage (Supporting Information Figure S2).
In addition, exposure to 1 led to significant perturbations in the progression of A375 cells throughout the cell cycle, as suggested by the significant increase over time in the percentage of cells residing in the G1 phase in treated vs. untreated cells (Figure 3D). Finally, in keeping with BCL-2 protein reduction (Figure 3B) and a trend in the expression of BCL-2 family transcripts suggesting a pro-apoptotic framework (Figure 1D), the exposure of A375 cells to 1 resulted in apoptosis induction, as revealed by the significant increase in the amounts of the cleaved form of PARP1 after 48 and 72 h of drug exposure compared to untreated cells (Figure 3E).

4. Discussion

Melanoma comprises a spectrum of malignancies characterized by targetable mutations in distinct oncogenic signaling pathways [3,9]. The development of BRAF and MEK inhibitors as treatment options for patients with BRAF-mutant tumors has contributed to improve the median overall survival of patients with metastatic melanoma [1,3]. Nevertheless, although these treatments lead to a quick initial response, most patients eventually relapse due to primary or secondary resistance mechanisms [1]. To date, distinct resistance mechanisms to targeted therapies have been described. These comprise the reactivation of MAPK signaling through NRAS mutations or amplification or alternative splicing of BRAF and MEK1/2 mutations [8]. Moreover, resistance may occur in a MAPK-independent fashion via the activation of alternative pathways (e.g., PI3K/AKT), the over-activation of tyrosine kinase receptors, and the evasion of apoptosis [1,8,9]. Hence, despite the advent of targeted therapies, melanoma still remains incurable and represents a tumor that is very difficult to treat [1,2].
In the last decades, G4 structures have emerged as attractive candidates for cancer therapy [22,48]. The growing interest in these oncogene promoters as therapeutic targets has fueled the rational design and development of small molecules that are able to recognize and stabilize such secondary structures [54,55]. A variety of G4-stabilizing small molecules identified in recent years are currently considered as fascinating weapons to therapeutically operate at the genomic level and represent a novel, though still challenging, tool to target oncogene promoters and downstream events in cancer cells.
In the present study, we showed that 1 was able to bind and stabilize G4 structures in vitro that form within the promoters of human KIT, BRAF, and BCL-2, which are three genes that may be relevant therapeutic targets in melanoma. Specifically, the biophysical data obtained by CD indicated that c-kit1, c-kit2, b-raf, and bcl-2 G4s were greatly stabilized by the compound to similar extents, with Tm values of >90 °C. Moreover, the Taq polymerase stop assay corroborated the ability of these sequences to be stabilized by 1 and also stalled polymerase progression.
At the cellular level, we observed that the exposure of BRAF-mutant melanoma cells to 1 resulted in a pronounced down-modulation of KIT, which was paralleled by the shutdown of the RAS/RAF/MAPK and PI3K/AKT signaling pathways, both at the transcriptomic and post-translational levels. Moreover, 1-treated cells were characterized by an expression profile of BCL-2 family members indicating a trend toward a pro-apoptotic milieu, which was also confirmed at the biochemical level by the remarkable reduction of BCL-2 protein amounts and the prominent accumulation of γ-H2AX and p21waf1, as well as the cleaved form of PARP-1. Conversely, our data revealed that NRAS-mutant melanoma cells were not particularly affected in oncogene-mediated signaling pathways, but rather showed a gene expression pattern in response to 1 exposure that was likely consistent with a phenotypic switch [1].
Altogether, our results supported the notion that the transcriptional activation of oncogenes in BRAF-mutant melanoma cells could be suppressed through the stabilization of G4s within gene promoters, and indicated that G4-stabilizing compounds could be used to interfere with oncogene expression, especially KIT and BCL-2, in a “double-hit” manner. However, we could not exclude that additional direct targets may have been amenable to 1-mediated G4 stabilization and, in keeping with previous evidence [56], our data suggested that the cell response to G4 ligands was largely dependent on the cellular context, as well as on the presence of specific genetic alterations.
Nonetheless, given the multifactoriality of cancer and the involvement of G4s in all its hallmarks [34], the low grade of specificity of the ligand toward one single G4 may be viewed as an advantage [22], and could be exploited to suppress the proliferation of tumor cells through activity at multiple targets. In fact, as a large body of preclinical and clinical evidence demonstrates, highly aggressive tumors such as melanomas may not be able to be defeated by monotherapies; rather, combinations of several signaling effectors should be envisioned [57].
While KIT may not represent a reliable primary therapeutic target in melanoma due to it being mutated in a small fraction of lesions [1,3], its relevance as a target may be envisaged in the context of acquired drug resistance of BRAF-mutant melanoma treated with targeted therapies. In fact, KIT activates a series of downstream effectors, such as the MAPK pathway with the activation of MAPK3/1 (ERK1/2) and the PI3K/AKT pathway [7]. Of note, the aberrant activation of these pathways plays a major role in conferring resistance to cell death (i.e., drug-induced apoptosis) in melanoma [7,9]. Notably, pro-death proteins of the BCL-2 family, such as BAD, BAX, BID, and BCL2L11, are inhibited by both AKT and MEK, whereas pro-survival factors belonging to the same family, such as BCL2L1, MCL1, and BCL-2, are activated by the MAPK pathway [7,9]. In this context, activation of MAPK or PI3K/AKT signaling pathways contributes to the relapse of melanoma patients treated with BRAF inhibitors [9]. Moreover, it was reported that BRAFV600E suppressed apoptosis by up-regulating the anti-apoptotic factor MCL1, whereas constitutively activated RAS/BRAF/MAPK signals inactivated the pro-apoptotic proteins BAD and BCL2L11. In addition, the activated PI3K/AKT pathway inhibited cell death by preventing the binding of BAD to BCL-2 or BCL2L1 and impeded the induction of apoptosis upon BRAF inhibition. Consequently, taking into account that BCL-2 family members may undergo regulation by MAPK or PI3K/AKT signaling pathways and that the activation of these pathways plays a key role in BRAF-inhibitor resistance, targeting downstream apoptotic proteins may represent a suitable option to overcome relapses [9]. In this context, different strategies to target Bcl-2 family members with promising pre-clinical evidence were reported [9]. For instance, it was shown that inhibiting BCL-2 in combination with MCL-1 inhibition or NOXA activation represented an effective approach to induce melanoma cell death, as well as that targeting BCL-2 with small molecule inhibitors (e.g., the BH3-mimetic ABT-737) sensitizes melanoma cells to BRAF inhibitors [9].

5. Conclusions

Overall, by affecting the expression levels of both KIT and BCL-2 proteins and, consequently, by impinging on MAPK and PI3K/AKT pathways, as well as on the balance between pro-death and pro-survival factors, our compound may be used to interfere with oncogene expression in melanoma cells. We thus present a starting point for the future development of therapeutic interventions based on G4 targeting that may be of particular relevance in the context of resistance to targeted therapies of BRAF-mutant melanoma cells.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/8/10/1274/s1, Scheme S1: Synthetic protocol for the preparation of core-extended naphthalene diimides; Figure S1: Structure and numbering of c-exNDIs; Figure S2: Immunofluorescence analysis of DNA-damage induction and Western immunoblotting of untreated (solvent) and 1-treated SKMEL-2 cells; Figure S3: CD analysis of c-myc, bcl-2, c-kit1, c-kit2, and b-raf G4 oligonucleotides; Figure S4: FRET analysis of c-kit2 and dsDNA oligonucleotides; Table S1: Oligonucleotides used in this study; Table S2: Cytotoxic activity of 1 on primary skin fibroblasts and cell lines established from human solid tumors of different histological origins; Table S3: Analysis of gene expression levels in melanoma cells upon 48 h of exposure to 1; Table S4: Over-representation analysis of the most relevant enriched pathways found in the list of differently expressed genes in treated vs. untreated melanoma cells; HPLC purity data; 1H- and 13C-NMR spectra of c-exNDIs 11 and 12; characterization of the c-exNDIs 11 and 12; supplementary references.

Author Contributions

Conceptualization, S.N.R. and M.F. (Marco Folini); data curation, M.F. (Marco Folini); formal analysis, M.R., M.F. (Marco Folini), and S.N.R.; funding acquisition, S.N.R. and M.F. (Marco Folini); investigation, M.R., M.T., G.C.-R., and F.D.; project administration, M.F. (Marco Folini) and S.N.R.; resources, M.F. (Marco Folini), N.Z., and S.N.R.; supervision, M.F. (Marco Folini) and S.N.R.; visualization, M.R., M.T., G.C.-R., M.F., (Marco Folini) and S.N.R.; writing—original draft preparation, G.C.-R., M.T., F.D., and M.F. (Marco Folini); writing—review and editing, S.N.R., M.F. (Marco Folini), N.Z., and M.F. (Mauro Freccero). M.F. (Marco Folini) and S.N.R. share senior authorship.

Funding

This research was funded by Fondazione AIRC per la Ricerca sul Cancro (Investigator grants #15194 to M.F. (Marco Folini) and #21850 to S.N.R.) and by the European Research Council (ERC-2013-CoG grant #615879 to S.N.R.).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kozar, I.; Margue, C.; Rothengatter, S.; Haan, C.; Kreis, S. Many ways to resistance: How melanoma cells evade targeted therapies. Biochim. Biophys. Acta. Rev. Cancer 2019, 1871, 313–322. [Google Scholar] [CrossRef] [PubMed]
  2. Najem, A.; Krayem, M.; Perdrix, A.; Kerger, J.; Awada, A.; Journe, F.; Ghanem, G. New drug combination strategies in melanoma: Current status and future directions. Anticancer Res. 2017, 37, 5941–5953. [Google Scholar] [PubMed]
  3. Kunz, M.; Hölzel, M. The impact of melanoma genetics on treatment response and resistance in clinical and experimental studies. Cancer Metastasis Rev. 2017, 36, 53–75. [Google Scholar] [CrossRef] [PubMed]
  4. Kunz, M. The genetic basis of new treatment modalities in melanoma. Curr. Drug Targets 2015, 16, 233–248. [Google Scholar] [CrossRef] [PubMed]
  5. Bollag, G.; Tsai, J.; Zhang, J.; Zhang, C.; Ibrahim, P.; Nolop, K.; Hirth, P. Vemurafenib: The first drug approved for BRAF-mutant cancer. Nat. Rev. Drug Discov. 2012, 11, 873–886. [Google Scholar] [CrossRef]
  6. Torres-Collado, A.X.; Knott, J.; Jazirehi, A.R. Reversal of resistance in targeted therapy of metastatic melanoma: Lessons learned from Vemurafenib (BRAFV600E-specific inhibitor). Cancers 2018, 10, 157. [Google Scholar] [CrossRef]
  7. Bastian, B.C.; Esteve-Puig, R. Targeting activated KIT signaling for melanoma therapy. J. Clin. Oncol. 2013, 31, 3288–3290. [Google Scholar] [CrossRef]
  8. Iams, W.T.; Sosman, J.A.; Chandra, S. Novel targeted therapies for metastatic melanoma. Cancer J. 2017, 23, 54–58. [Google Scholar] [CrossRef]
  9. Mukherjee, N.; Schwan, J.V.; Fujita, M.; Norris, D.A.; Shellman, Y.G. Alternative treatments for melanoma: Targeting BCL-2 family members to de-bulk and kill cancer stem cells. J. Invest. Dermatol. 2015, 135, 2155–2161. [Google Scholar] [CrossRef]
  10. Aird, D.; Teng, T.; Huang, C.-L.; Pazolli, E.; Banka, D.; Cheung-Ong, K.; Eifert, C.; Furman, C.; Wu, Z.J.; Seiler, M.; et al. Sensitivity to splicing modulation of BCL2 family genes defines cancer therapeutic strategies for splicing modulators. Nat. Commun. 2019, 10, 137. [Google Scholar] [CrossRef]
  11. Raghavendra, N.M.; Pingili, D.; Kadasi, S.; Mettu, A.; Prasad, S.V.U.M. Dual or multi-targeting inhibitors: The next generation anticancer agents. Eur. J. Med. Chem. 2018, 143, 1277–1300. [Google Scholar] [CrossRef] [PubMed]
  12. Huppert, J.L. Structure, location and interactions of G-quadruplexes. FEBS J. 2010, 277, 3452–3458. [Google Scholar] [CrossRef] [PubMed]
  13. Campbell, N.H.; Neidle, S. G-quadruplexes and metal ions. Met. Ions Life Sci. 2012, 10, 119–134. [Google Scholar] [PubMed]
  14. Burge, S.; Parkinson, G.N.; Hazel, P.; Todd, A.K.; Neidle, S. Quadruplex DNA: Sequence, topology and structure. Nucleic Acids Res. 2006, 34, 5402–5415. [Google Scholar] [CrossRef] [PubMed]
  15. Rhodes, D.; Lipps, H.J. G-quadruplexes and their regulatory roles in biology. Nucleic Acids Res. 2015, 43, 8627–8637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Shin, Y.-J.; Kumarasamy, V.; Camacho, D.; Sun, D. Involvement of G-quadruplex structures in regulation of human RET gene expression by small molecules in human medullary thyroid carcinoma TT cells. Oncogene 2015, 34, 1292–1299. [Google Scholar] [CrossRef]
  17. Zhou, J.; Fleming, A.M.; Averill, A.M.; Burrows, C.J.; Wallace, S.S. The NEIL glycosylases remove oxidized guanine lesions from telomeric and promoter quadruplex DNA structures. Nucleic Acids Res. 2015, 43, 4039–4054. [Google Scholar] [CrossRef] [Green Version]
  18. Berardinelli, F.; Siteni, S.; Tanzarella, C.; Stevens, M.F.; Sgura, A.; Antoccia, A. The G-quadruplex-stabilising agent RHPS4 induces telomeric dysfunction and enhances radiosensitivity in glioblastoma cells. DNA Repair 2015, 25, 104–115. [Google Scholar] [CrossRef]
  19. Perrone, R.; Nadai, M.; Frasson, I.; Poe, J.A.; Butovskaya, E.; Smithgall, T.E.; Palumbo, M.; Palù, G.; Richter, S.N. A dynamic G-quadruplex region regulates the HIV-1 long terminal repeat promoter. J. Med. Chem. 2013, 56, 6521–6530. [Google Scholar] [CrossRef]
  20. Paeschke, K.; Bochman, M.L.; Garcia, P.D.; Cejka, P.; Friedman, K.L.; Kowalczykowski, S.C.; Zakian, V.A. Pif1 family helicases suppress genome instability at G-quadruplex motifs. Nature 2013, 497, 458–462. [Google Scholar] [CrossRef] [Green Version]
  21. Takahama, K.; Takada, A.; Tada, S.; Shimizu, M.; Sayama, K.; Kurokawa, R.; Oyoshi, T. Regulation of telomere length by G-quadruplex telomere DNA- and TERRA-binding protein TLS/FUS. Chem. Biol. 2013, 20, 341–350. [Google Scholar] [CrossRef] [PubMed]
  22. Neidle, S. Quadruplex nucleic acids as targets for anticancer therapeutics. Nat. Rev. Chem. 2017, 1, 1–10. [Google Scholar] [CrossRef]
  23. Simone, R.; Balendra, R.; Moens, T.G.; Preza, E.; Wilson, K.M.; Heslegrave, A.; Woodling, N.S.; Niccoli, T.; Gilbert-Jaramillo, J.; Abdelkarim, S.; et al. G-quadruplex-binding small molecules ameliorate C9orf72 FTD/ALS pathology in vitro and in vivo. EMBO Mol. Med. 2018, 10, 22–31. [Google Scholar] [CrossRef] [PubMed]
  24. Zhou, B.; Geng, Y.; Liu, C.; Miao, H.; Ren, Y.; Xu, N.; Shi, X.; You, Y.; Lee, T.; Zhu, G. Characterizations of distinct parallel and antiparallel G-quadruplexes formed by two-repeat ALS and FTD related GGGGCC sequence. Sci. Rep. 2018, 8, 2366. [Google Scholar] [CrossRef]
  25. Mendoza, O.; Bourdoncle, A.; Boulé, J.-B.; Brosh, R.M.; Mergny, J.-L. G-quadruplexes and helicases. Nucleic Acids Res. 2016, 44, 1989–2006. [Google Scholar] [CrossRef] [Green Version]
  26. Scalabrin, M.; Frasson, I.; Ruggiero, E.; Perrone, R.; Tosoni, E.; Lago, S.; Tassinari, M.; Palù, G.; Richter, S.N. The cellular protein hnRNP A2/B1 enhances HIV-1 transcription by unfolding LTR promoter G-quadruplexes. Sci. Rep. 2017, 7, 45244. [Google Scholar] [CrossRef]
  27. Henderson, A.; Wu, Y.; Huang, Y.C.; Chavez, E.A.; Platt, J.; Johnson, F.B.; Brosh, R.M.; Sen, D.; Lansdorp, P.M. Detection of G-quadruplex DNA in mammalian cells. Nucleic Acids Res. 2014, 42, 860–869. [Google Scholar] [CrossRef]
  28. Biffi, G.; Tannahill, D.; McCafferty, J.; Balasubramanian, S. Quantitative visualization of DNA G-quadruplex structures in human cells. Nat. Chem. 2013, 5, 182–186. [Google Scholar] [CrossRef]
  29. Hänsel-Hertsch, R.; Beraldi, D.; Lensing, S.V.; Marsico, G.; Zyner, K.; Parry, A.; Di Antonio, M.; Pike, J.; Kimura, H.; Narita, M.; et al. G-quadruplex structures mark human regulatory chromatin. Nat. Genet. 2016, 48, 1267–1272. [Google Scholar] [CrossRef] [Green Version]
  30. Yang, S.Y.; Lejault, P.; Chevrier, S.; Boidot, R.; Robertson, A.G.; Wong, J.M.Y.; Monchaud, D. Transcriptome-wide identification of transient RNA G-quadruplexes in human cells. Nat. Commun. 2018, 9, 4730. [Google Scholar] [CrossRef]
  31. Renard, I.; Grandmougin, M.; Roux, A.; Yang, S.Y.; Lejault, P.; Pirrotta, M.; Wong, J.M.Y.; Monchaud, D. Small-molecule affinity capture of DNA/RNA quadruplexes and their identification in vitro and in vivo through the G4RP protocol. Nucleic Acids Res. 2019. [Google Scholar] [CrossRef] [PubMed]
  32. Neidle, S. Human telomeric G-quadruplex: The current status of telomeric G-quadruplexes as therapeutic targets in human cancer. FEBS J. 2010, 277, 1118–1125. [Google Scholar] [CrossRef] [PubMed]
  33. Huppert, J.L.; Balasubramanian, S. G-quadruplexes in promoters throughout the human genome. Nucleic Acids Res. 2007, 35, 406–413. [Google Scholar] [CrossRef] [PubMed]
  34. Brooks, T.A.; Kendrick, S.; Hurley, L. Making sense of G-quadruplex and i-motif functions in oncogene promoters. FEBS J. 2010, 277, 3459–3469. [Google Scholar] [CrossRef] [Green Version]
  35. Siddiqui-Jain, A.; Grand, C.L.; Bearss, D.J.; Hurley, L.H. Direct evidence for a G-quadruplex in a promoter region and its targeting with a small molecule to repress c-MYC transcription. Proc. Natl. Acad. Sci. USA 2002, 99, 11593–11598. [Google Scholar] [CrossRef]
  36. Musso, L.; Mazzini, S.; Rossini, A.; Castagnoli, L.; Scaglioni, L.; Artali, R.; Di Nicola, M.; Zunino, F.; Dallavalle, S. c-MYC G-quadruplex binding by the RNA polymerase I inhibitor BMH-21 and analogues revealed by a combined NMR and biochemical Approach. Biochim. Biophys. Acta Gen. Subj. 2018, 1862, 615–629. [Google Scholar] [CrossRef]
  37. O’Hagan, M.P.; Morales, J.C.; Galan, M.C. Binding and beyond: What else can G-quadruplex ligands do? Eur. J. Org. Chem. 2019, 2019, 4995–5017. [Google Scholar] [CrossRef]
  38. Xu, H.; Di Antonio, M.; McKinney, S.; Mathew, V.; Ho, B.; O’Neil, N.J.; Dos Santos, N.; Silvester, J.; Wei, V.; Garcia, J.; et al. CX-5461 is a DNA G-quadruplex stabilizer with selective lethality in BRCA1/2 deficient tumours. Nat. Commun. 2017, 8, 14432. [Google Scholar] [CrossRef]
  39. Balasubramanian, S.; Hurley, L.H.; Neidle, S. Targeting G-quadruplexes in gene promoters: A novel anticancer strategy? Nat. Rev. Drug Discov. 2011, 10, 261–275. [Google Scholar] [CrossRef]
  40. Pirota, V.; Nadai, M.; Doria, F.; Richter, S. Naphthalene diimides as multimodal G-quadruplex-selective ligands. Molecules 2019, 24, 426. [Google Scholar] [CrossRef]
  41. Perrone, R.; Doria, F.; Butovskaya, E.; Frasson, I.; Botti, S.; Scalabrin, M.; Lago, S.; Grande, V.; Nadai, M.; Freccero, M.; et al. Synthesis, binding and antiviral properties of potent core-extended naphthalene diimides targeting the HIV-1 long terminal repeat promoter G-quadruplexes. J. Med. Chem. 2015, 58, 9639–9652. [Google Scholar] [CrossRef] [PubMed]
  42. Zuffo, M.; Guédin, A.; Leriche, E.-D.; Doria, F.; Pirota, V.; Gabelica, V.; Mergny, J.-L.; Freccero, M. More is not always better: Finding the right trade-off between affinity and selectivity of a G-quadruplex ligand. Nucleic Acids Res. 2018, 46, e115. [Google Scholar] [CrossRef] [PubMed]
  43. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  44. Greenfield, N.J. Using circular dichroism collected as a function of temperature to determine the thermodynamics of protein unfolding and binding interactions. Nat. Protoc. 2006, 1, 2527–2535. [Google Scholar] [CrossRef] [PubMed]
  45. Rachwal, P.A.; Fox, K.R. Quadruplex melting. Methods 2007, 43, 291–301. [Google Scholar] [CrossRef] [PubMed]
  46. Doria, F.; Nadai, M.; Folini, M.; Di Antonio, M.; Germani, L.; Percivalle, C.; Sissi, C.; Zaffaroni, N.; Alcaro, S.; Artese, A.; et al. Hybrid ligand-alkylating agents targeting telomeric G-quadruplex structures. Org. Biomol. Chem. 2012, 10, 2798–2806. [Google Scholar] [CrossRef]
  47. Tassinari, M.; Cimino-Reale, G.; Nadai, M.; Doria, F.; Butovskaya, E.; Recagni, M.; Freccero, M.; Zaffaroni, N.; Richter, S.N.; Folini, M. Down-regulation of the androgen receptor by G-quadruplex ligands sensitizes castration-resistant prostate cancer cells to Enzalutamide. J. Med. Chem. 2018, 61, 8625–8638. [Google Scholar] [CrossRef]
  48. Cimino-Reale, G.; Zaffaroni, N.; Folini, M. Emerging role of G-quadruplex DNA as target in anticancer therapy. Curr. Pharm. Des. 2016, 22, 6612–6624. [Google Scholar] [CrossRef]
  49. Wei, D.; Todd, A.K.; Zloh, M.; Gunaratnam, M.; Parkinson, G.N.; Neidle, S. Crystal structure of a promoter sequence in the B-raf gene reveals an intertwined dimer quadruplex. J. Am. Chem. Soc. 2013, 135, 19319–19329. [Google Scholar] [CrossRef]
  50. You, H.; Wu, J.; Shao, F.; Yan, J. Stability and kinetics of c- MYC promoter G-quadruplexes studied by single-molecule manipulation. J. Am. Chem. Soc. 2015, 137, 2424–2427. [Google Scholar] [CrossRef]
  51. Agrawal, P.; Lin, C.; Mathad, R.I.; Carver, M.; Yang, D. The major G-quadruplex formed in the human BCL-2 proximal promoter adopts a parallel structure with a 13-nt loop in K + solution. J. Am. Chem. Soc. 2014, 136, 1750–1753. [Google Scholar] [CrossRef] [PubMed]
  52. Fernando, H.; Reszka, A.P.; Huppert, J.; Ladame, S.; Rankin, S.; Venkitaraman, A.R.; Neidle, S.; Balasubramanian, S. A conserved quadruplex motif located in a transcription activation site of the human c-kit oncogene. Biochemistry 2006, 45, 7854–7860. [Google Scholar] [CrossRef] [PubMed]
  53. Rankin, S.; Reszka, A.P.; Huppert, J.; Zloh, M.; Parkinson, G.N.; Todd, A.K.; Ladame, S.; Balasubramanian, S.; Neidle, S. Putative DNA quadruplex formation within the human c-kit oncogene. J. Am. Chem. Soc. 2005, 127, 10584–10589. [Google Scholar] [CrossRef] [PubMed]
  54. Neidle, S. The structures of quadruplex nucleic acids and their drug complexes. Curr. Opin. Struct. Biol. 2009, 19, 239–250. [Google Scholar] [CrossRef] [PubMed]
  55. Monchaud, D.; Teulade-Fichou, M.-P. A hitchhiker’s guide to G-quadruplex ligands. Org. Biomol. Chem. 2008, 6, 627–636. [Google Scholar] [CrossRef]
  56. Nadai, M.; Cimino-Reale, G.; Sattin, G.; Doria, F.; Butovskaya, E.; Zaffaroni, N.; Freccero, M.; Palumbo, M.; Richter, S.N.; Folini, M. Assessment of gene promoter G-quadruplex binding and modulation by a naphthalene diimide derivative in tumor cells. Int. J. Oncol. 2015, 46, 369–380. [Google Scholar] [CrossRef]
  57. Eberle, J.; Hossini, A.M. Expression and function of bcl-2 proteins in melanoma. Curr. Genom. 2008, 9, 409–419. [Google Scholar] [CrossRef]
Cells 08 01274 g001 550
Figure 1. (A) Time-dependent assessment of cell growth in untreated A375 and SKMEL-2 cells. Data are reported as the recorded absorbance (Abs) as a function of time and represent mean values ± s.d. from at least three independent experiments. The curves represent the exponential growth equation obtained by nonlinear fit of data points using GraphPad Prism 5.01 (GraphPad Software Inc., San Diego, California, United States). (B) Dose-response curves of melanoma cells exposed to increasing concentrations (0.1–10,000 nM) of 1. Data are reported as the percentage of viable cells with respect to untreated cells as a function of the Log10 of compound concentrations (for details see Supporting Information Table S2) and represent mean values ± s.d. from at least three independent experiments. (C) Representative photomicrographs showing the uptake of the G4 ligand (red) and the occurrence of G4 structures (green) in A375 and SKMEL-2 cells exposed to c-exNDI 1 (IC50) for 48 h and analyzed by fluorescence microscopy. The panels on the bottom show untreated cells. Nuclei were counterstained with DAPI. Magnification: ×100; bars: 10 μm. (D) Analysis of gene expression levels in A375 cells upon 48 h of exposure to 1 (IC50). Data are reported as Relative Quantity (mean RQ values ± s.d.) in treated vs. untreated cells, according to the 2−ΔΔCt method. Dotted line: calibrator sample. * p < 0.05; ** p < 0.01; *** p < 0.001 (Student’s t-test).
Figure 1. (A) Time-dependent assessment of cell growth in untreated A375 and SKMEL-2 cells. Data are reported as the recorded absorbance (Abs) as a function of time and represent mean values ± s.d. from at least three independent experiments. The curves represent the exponential growth equation obtained by nonlinear fit of data points using GraphPad Prism 5.01 (GraphPad Software Inc., San Diego, California, United States). (B) Dose-response curves of melanoma cells exposed to increasing concentrations (0.1–10,000 nM) of 1. Data are reported as the percentage of viable cells with respect to untreated cells as a function of the Log10 of compound concentrations (for details see Supporting Information Table S2) and represent mean values ± s.d. from at least three independent experiments. (C) Representative photomicrographs showing the uptake of the G4 ligand (red) and the occurrence of G4 structures (green) in A375 and SKMEL-2 cells exposed to c-exNDI 1 (IC50) for 48 h and analyzed by fluorescence microscopy. The panels on the bottom show untreated cells. Nuclei were counterstained with DAPI. Magnification: ×100; bars: 10 μm. (D) Analysis of gene expression levels in A375 cells upon 48 h of exposure to 1 (IC50). Data are reported as Relative Quantity (mean RQ values ± s.d.) in treated vs. untreated cells, according to the 2−ΔΔCt method. Dotted line: calibrator sample. * p < 0.05; ** p < 0.01; *** p < 0.001 (Student’s t-test).
Cells 08 01274 g001
Cells 08 01274 g002 550
Figure 2. Representative circular dichroism (CD) spectra of G4s formed by c-myc, bcl-2, c-kit1, c-kit2, and b-raf oligonucleotides (4 μM) recorded in the absence (A) or presence (B) of a 4-fold excess of 1 (16 μM). (C) Image of a typical Taq polymerase stop assay. The b-raf, bcl-2, c-myc, and c-kit1+2 templates were amplified by Taq polymerase in the absence (lanes 1) and presence of 10 mM K+, alone (lanes 2) or combined, with increasing amounts (37.5, 75.0, and 150.0 nM) of 1 (lanes 3–5). A template (non-G4 cnt) made of a scramble sequence unable to fold into G4 was also used as an internal control. Lane M: Ladder of markers obtained by the Maxam and Gilbert sequencing carried out on the amplified strand complementary to the template strand. Lane P: Unreacted labelled primer. Vertical bars indicate G4-specific Taq polymerase stop sites. (D) Quantification of the intensity of G4 stop bands obtained in the Taq polymerase stop assay. G4 stop band intensities are expressed as percentages with respect to the total elongated products.
Figure 2. Representative circular dichroism (CD) spectra of G4s formed by c-myc, bcl-2, c-kit1, c-kit2, and b-raf oligonucleotides (4 μM) recorded in the absence (A) or presence (B) of a 4-fold excess of 1 (16 μM). (C) Image of a typical Taq polymerase stop assay. The b-raf, bcl-2, c-myc, and c-kit1+2 templates were amplified by Taq polymerase in the absence (lanes 1) and presence of 10 mM K+, alone (lanes 2) or combined, with increasing amounts (37.5, 75.0, and 150.0 nM) of 1 (lanes 3–5). A template (non-G4 cnt) made of a scramble sequence unable to fold into G4 was also used as an internal control. Lane M: Ladder of markers obtained by the Maxam and Gilbert sequencing carried out on the amplified strand complementary to the template strand. Lane P: Unreacted labelled primer. Vertical bars indicate G4-specific Taq polymerase stop sites. (D) Quantification of the intensity of G4 stop bands obtained in the Taq polymerase stop assay. G4 stop band intensities are expressed as percentages with respect to the total elongated products.
Cells 08 01274 g002
Cells 08 01274 g003 550
Figure 3. (A) Analysis of KIT and MYC mRNA expression levels in A375 and SKMEL-2 cells upon 48 h of exposure to 1 (IC50). Data are reported as Relative Quantity (mean RQ values ± s.d.) in treated vs. untreated cells, according to the 2−ΔΔCt method. Dotted line: calibrator sample. (B) Representative Western immunoblotting showing the effects of 48 and 72 h of exposure of A375 cells to 1 (IC50) on the amounts of the indicated proteins. Cropped images of selected proteins are shown. Beta-actin (ACTB) was used to ensure equal protein loading. The graph on the right shows the extent of KIT (black bars) and BCL2 (white bars) protein inhibition in A375 cells exposed to 1 (IC50) for 48 and 72 h. Data are reported as the inhibition of protein expression in treated vs. untreated cells, and represent the mean values ± s.d. from at least three independent experiments. (C) Representative Western immunoblotting showing the amounts of the indicated proteins in untreated (−) A375 cells and upon 72 h of exposure (+) to 1 (IC50). Vinculin (VCL) was used to ensure equal protein loading. Cropped images of selected proteins are shown. The graph on the bottom reports the quantification of protein levels in untreated (white bars) and 1-treated (black bars) A375 cells. Data are reported as relative protein amounts following the normalization of densitometric signals toward Vinculin, and represent the mean values ± s.d. (D) Time-dependent assessment of the cell cycle in untreated and 1-treated A375 cells. Data are reported as the percentages of cells in G0/G1 (red), S (green), and G2/M (black) phases of the cell cycle, and represent the mean values ± s.d. from at least three independent experiments. (E) Representative Western immunoblotting showing the amounts of full length and cleaved PARP1 in untreated (−) A375 cells and upon 48 and 72 h of exposure (+) to 1 (IC50). Beta-actin (ACTB) was used to ensure equal protein loading. Cropped images of selected proteins are shown. The graph on the right reports the amounts of cleaved PARP1 in untreated (white bars) and 1-treated (black bars) cells. Data are reported as relative protein amounts vs. full length PARP1, and represent the mean values ± s.d. * p < 0.05; ** p < 0.01; *** p < 0.001 (Student’s t-test).
Figure 3. (A) Analysis of KIT and MYC mRNA expression levels in A375 and SKMEL-2 cells upon 48 h of exposure to 1 (IC50). Data are reported as Relative Quantity (mean RQ values ± s.d.) in treated vs. untreated cells, according to the 2−ΔΔCt method. Dotted line: calibrator sample. (B) Representative Western immunoblotting showing the effects of 48 and 72 h of exposure of A375 cells to 1 (IC50) on the amounts of the indicated proteins. Cropped images of selected proteins are shown. Beta-actin (ACTB) was used to ensure equal protein loading. The graph on the right shows the extent of KIT (black bars) and BCL2 (white bars) protein inhibition in A375 cells exposed to 1 (IC50) for 48 and 72 h. Data are reported as the inhibition of protein expression in treated vs. untreated cells, and represent the mean values ± s.d. from at least three independent experiments. (C) Representative Western immunoblotting showing the amounts of the indicated proteins in untreated (−) A375 cells and upon 72 h of exposure (+) to 1 (IC50). Vinculin (VCL) was used to ensure equal protein loading. Cropped images of selected proteins are shown. The graph on the bottom reports the quantification of protein levels in untreated (white bars) and 1-treated (black bars) A375 cells. Data are reported as relative protein amounts following the normalization of densitometric signals toward Vinculin, and represent the mean values ± s.d. (D) Time-dependent assessment of the cell cycle in untreated and 1-treated A375 cells. Data are reported as the percentages of cells in G0/G1 (red), S (green), and G2/M (black) phases of the cell cycle, and represent the mean values ± s.d. from at least three independent experiments. (E) Representative Western immunoblotting showing the amounts of full length and cleaved PARP1 in untreated (−) A375 cells and upon 48 and 72 h of exposure (+) to 1 (IC50). Beta-actin (ACTB) was used to ensure equal protein loading. Cropped images of selected proteins are shown. The graph on the right reports the amounts of cleaved PARP1 in untreated (white bars) and 1-treated (black bars) cells. Data are reported as relative protein amounts vs. full length PARP1, and represent the mean values ± s.d. * p < 0.05; ** p < 0.01; *** p < 0.001 (Student’s t-test).
Cells 08 01274 g003
Table
Table 1. Cytotoxic activity of core-extended naphthalene diimides (c-exNDIs) in A375 melanoma cells.
Table 1. Cytotoxic activity of core-extended naphthalene diimides (c-exNDIs) in A375 melanoma cells.
CompoundIC50 1 (nM)
18 ± 1
220 ± 2
350 ± 3
474 ± 6
5>500
6>500
7101 ± 11
894 ± 6
9>500
10>500
1154 ± 2
12104 ± 3
13415 ± 36
14195 ± 19
1 IC50 values (concentration of the compound that inhibited 50% of cell growth) were calculated from the dose-response curves obtained upon 48 h of exposure of A375 cells to increasing concentrations of the compounds and reported as mean values ± s.d. from at least three independent experiments.
Table
Table 2. Stabilization of oncogene promoter G4-folded sequences (4 μM) by compound 1 (16 μM), measured by CD thermal unfolding analysis.
Table 2. Stabilization of oncogene promoter G4-folded sequences (4 μM) by compound 1 (16 μM), measured by CD thermal unfolding analysis.
OligonucleotideTm Values 1 (°C)
No Compoundc-exNDI 1
λ = 264 nmλ = 290 nmλ = 264 nmλ = 290 nm
blc-267.5 ± 1.3/>90.0/
c-myc58.8 ± 0.7/>90.0/
b-raf>90.050.3 ± 0.3>90.068.7 ± 0.8
c-kit 148.0 ± 0.3/>90.0/
c-kit 257.0 ± 0.3/>90.0/
1 Melting temperature of oncogene promoter G4-folded sequences in the absence or presence of a 4-fold excess of 1. Data are reported as mean values ± sd from two independent experiments.

Share and Cite

MDPI and ACS Style

Recagni, M.; Tassinari, M.; Doria, F.; Cimino-Reale, G.; Zaffaroni, N.; Freccero, M.; Folini, M.; Richter, S.N. The Oncogenic Signaling Pathways in BRAF-Mutant Melanoma Cells Are Modulated by Naphthalene Diimide-Like G-Quadruplex Ligands. Cells 2019, 8, 1274. https://doi.org/10.3390/cells8101274

AMA Style

Recagni M, Tassinari M, Doria F, Cimino-Reale G, Zaffaroni N, Freccero M, Folini M, Richter SN. The Oncogenic Signaling Pathways in BRAF-Mutant Melanoma Cells Are Modulated by Naphthalene Diimide-Like G-Quadruplex Ligands. Cells. 2019; 8(10):1274. https://doi.org/10.3390/cells8101274

Chicago/Turabian Style

Recagni, Marta, Martina Tassinari, Filippo Doria, Graziella Cimino-Reale, Nadia Zaffaroni, Mauro Freccero, Marco Folini, and Sara N. Richter. 2019. "The Oncogenic Signaling Pathways in BRAF-Mutant Melanoma Cells Are Modulated by Naphthalene Diimide-Like G-Quadruplex Ligands" Cells 8, no. 10: 1274. https://doi.org/10.3390/cells8101274

APA Style

Recagni, M., Tassinari, M., Doria, F., Cimino-Reale, G., Zaffaroni, N., Freccero, M., Folini, M., & Richter, S. N. (2019). The Oncogenic Signaling Pathways in BRAF-Mutant Melanoma Cells Are Modulated by Naphthalene Diimide-Like G-Quadruplex Ligands. Cells, 8(10), 1274. https://doi.org/10.3390/cells8101274

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop