Modes of Action of a Novel c-MYC Inhibiting 1,2,4-Oxadiazole Derivative in Leukemia and Breast Cancer Cells

The c-MYC oncogene regulates multiple cellular activities and is a potent driver of many highly aggressive human cancers, such as leukemia and triple-negative breast cancer. The oxadiazole class of compounds has gained increasing interest for its anticancer activities. The aim of this study was to investigate the molecular modes of action of a 1,2,4-oxadiazole derivative (ZINC15675948) as a c-MYC inhibitor. ZINC15675948 displayed profound cytotoxicity at the nanomolar range in CCRF-CEM leukemia and MDA-MB-231-pcDNA3 breast cancer cells. Multidrug-resistant sublines thereof (i.e., CEM/ADR5000 and MDA-MB-231-BCRP) were moderately cross-resistant to this compound (<10-fold). Molecular docking and microscale thermophoresis revealed a strong binding of ZINC15675948 to c-MYC by interacting close to the c-MYC/MAX interface. A c-MYC reporter assay demonstrated that ZINC15675948 inhibited c-MYC activity. Western blotting and qRT-PCR showed that c-MYC expression was downregulated by ZINC15675948. Applying microarray hybridization and signaling pathway analyses, ZINC15675948 affected signaling routes downstream of c-MYC in both leukemia and breast cancer cells as demonstrated by the induction of DNA damage using single cell gel electrophoresis (alkaline comet assay) and induction of apoptosis using flow cytometry. ZINC15675948 also caused G2/M phase and S phase arrest in CCRF-CEM cells and MDA-MB-231-pcDNA3 cells, respectively, accompanied by the downregulation of CDK1 and p-CDK2 expression using western blotting. Autophagy induction was observed in CCRF-CEM cells but not MDA-MB-231-pcDNA3 cells. Furthermore, microarray-based mRNA expression profiling indicated that ZINC15675948 may target c-MYC-regulated ubiquitination, since the novel ubiquitin ligase (ELL2) was upregulated in the absence of c-MYC expression. We propose that ZINC15675948 is a promising natural product-derived compound targeting c-MYC in c-MYC-driven cancers through DNA damage, cell cycle arrest, and apoptosis.


Introduction
Cancer is mainly caused by genomic alterations that result from the loss of tumor suppressor genes and the activation of oncogenes [1]. The MYC gene is one of the oncogenes that has been described as a master regulator of gene expression in multiple biological processes [2,3]. This gene encodes a family of basic helix-loop-helix zipper (bHLHZip) proteins consisting of c-MYC, N-MYC, and L-MYC [4]. MYC acts as a transcription factor that binds to its obligatory partner, MYC-associated factor X (MAX). The MYC-MAX heterodimer activates a large number of genes by binding to E box sequences (5 -CACGTG-3 ) within gene promoters and enhancers [5].
MYC-regulated transcription is tightly controlled in non-transformed cells. However, MYC is estimated to be overexpressed up to 70% in various human cancer types [6].
MYC contributes to several hallmarks of cancer, including the escape from programmed cell death, promoting sustainable proliferation, genome instability, escape from immunosurveillance, and change of cellular metabolism [7,8]. The major mechanisms of MYC deregulation are gene alternation and the activation of upstream MYC-related signaling pathways (e.g., NOTCH, WNT, and EGFR) [9]. Moreover, MYC is an unstable protein with a short half-life [10], and its stability is also regulated by post-translational modifications [11]. For example, MYC gene amplification largely drives breast carcinogenesis [12], and the chromosomal translocation of the MYC gene contributes to the development of T-cell acute leukemia [13]. In mice with MYC-induced hematological cancers, an inactivation of the MYC transgene led to tumor regression [14]. To date, there is no approved MYC inhibitor, despite several promising MYC inhibitors being investigated in preclinical or clinical studies [8]. Therefore, there is an ongoing quest for novel effective MYC inhibitors for an improvement of cancer therapy with targeted drugs.
Natural products are a promising resource for drug discovery [15]. Oxadiazole represents a five-membered heterocyclic scaffold with one oxygen and two nitrogen atoms. Oxadiazole-based compounds are a rapidly growing field in drug development. They are fundamental pharmacophores due to their stability in aqueous medium and are commonly employed as bioisosteric substitutes [16,17]. Based on the position of oxygen and nitrogen in the ring, oxadiazoles are classified into four regioisomeric structures [18]. In addition, 1,2,4-Oxadiazoles received considerable attention as witnessed by increasing numbers of published studies [19]. The first natural products from the class of 1,2,4-oxadiazoles were phidianidine A and B that were isolated from marine mollusk Phidiana militaris by Alder and Hancock [20]. Interestingly, phidianidine A was identified as an antifoulant agent which is a chemical defense of slow-moving marine organisms to deter predators, other colonizers, or competitors [21]. In the past two decades, a great number of 1,2,4oxadiazole derivates were synthesized and studied for their pharmacological properties, which exhibit numerous bioactivities, including anti-cancer [22,23], anti-inflammatory [24], anti-bacterial [25], anti-viral [26], anti-malarial [27], anti-diabetic [28], and anti-Alzheimer activities [29]. Zibotentan represents an example of an oxadiazole that reached a clinical phase III trial for the treatment of hormone-resistant prostate carcinoma [30]. Studies revealed that 1,2,4-oxadiazoles derivatives exhibited more powerful anticancer activities than established drugs (e.g., doxorubicin and etoposide) by inhibiting cell proliferation which make them appealing as prospective therapeutic candidates [16,18]. However, the modes of action of 1,2,4-oxadiazoles derivatives in cancer have not been fully understood yet.
Breast cancer is the most common malignancy in females worldwide. Triple-negative breast cancer (TNBC) is characterized by the absence of estrogen (ER)/progesterone (PR) expression and human epidermal growth factor receptor-2 (HER2) amplification [31]. MYC is remarkably elevated in TNBC compared with other breast cancer subtypes. TNBC is highly invasive leading to a poor five-year survival rate and distant recurrence rates [12,32]. On the other hand, acute lymphoblastic leukemia (ALL) is an aggressive hematological cancer frequently appearing in children and adolescents. A multitude of evidence supports the important role of MYC in the initiation and progression of ALL [33].
The aim of this study was to investigate the 1,2,4-oxadiazoles derivative ZINC15675948 as a novel c-MYC inhibitor. We determined the anticancer activity of ZINC15675948 against drug-sensitive and -resistant leukemia (CCRF-CEM and CEM/ADR5000) and triple-negative breast cancer (MDA-MB-231-pcDNA3 and MDA-MB-BCRP) cell lines. The effects of ZINC15675948 were studied by using western blot, qRT-PCR, molecular docking, microscale thermophoresis, and c-MYC reporter assay. By means of microarray hybridization and Ingenuity Pathway Analysis, we explored the underlying modes of action of ZINC15675948 regarding cell cycle, apoptosis, autophagy, and DNA damage. In addition, the interactions of ZINC15675948 with the multidrug-resistance-mediating ATP-binding cassette transporters' P-glycoprotein and breast cancer resistance protein were investigated using molecular docking.

Growth Inhibition Assay
ZINC15675948 exhibited profound cytotoxicity toward both leukemia and breast cancer cell lines using resazurin reduction assays. Doxorubicin was used as a positive control [34]. CCRF-CEM and MDA-MB-231-pcDNA3 cells were extremely sensitive to ZINC15675948 with IC 50 values of 0.008 ± 0.001 µM, and 0.08 ± 0.004 µM, respectively. It is interesting that CEM/ADR5000 cells as well as MDA-MB-BCRP cells displayed crossresistance toward ZINC15675948 (degrees of resistance: 8.37 and 9.0, respectively). The dose-response curves are shown in Figure 1C,D, and the IC 50 values and degrees of resistance are presented in Table 1.
resistance protein were investigated using molecular docking.

Growth Inhibition Assay
ZINC15675948 exhibited profound cytotoxicity toward both leukemia and breast cancer cell lines using resazurin reduction assays. Doxorubicin was used as a positive control [34]. CCRF-CEM and MDA-MB-231-pcDNA3 cells were extremely sensitive to ZINC15675948 with IC50 values of 0.008 ± 0.001 µM, and 0.08 ± 0.004 µM, respectively. It is interesting that CEM/ADR5000 cells as well as MDA-MB-BCRP cells displayed crossresistance toward ZINC15675948 (degrees of resistance: 8.37 and 9.0, respectively). The dose-response curves are shown in Figure 1C,D, and the IC50 values and degrees of resistance are presented in Table 1.
As multidrug-resistant cell lines were cross-resistant to ZINC15675948, further investigations using ZINC15675948 were performed only in CCRF-CEM and MDA-MB-231-pcDNA3 cells.   As multidrug-resistant cell lines were cross-resistant to ZINC15675948, further investigations using ZINC15675948 were performed only in CCRF-CEM and MDA-MB-231-pcDNA3 cells.

Molecular Docking
The oncogene MYC is deregulated in various human cancers and drives several cancerrelated hallmarks. Despite its unquestionable contribution to cancer development, MYC has been regarded as undruggable, and there are only a few MYC inhibitors so far [35,36]. The aim of this study was targeting c-MYC by ZINC15675948, as demonstrated below with different in vitro verifications and pathway analysis of microarray data. To investigate the possible interaction and binding affinity of ZINC15675948 to c-MYC, molecular docking was performed using AutoDock 4.2.6. As shown in Figure 2A, the binding site of ZINC15675948 bound to c-MYC was almost the same as for the known MYC-inhibiting control drugs. ZINC15675948 was particularly near to 10058-F4 and 10074-A4, and Arg925 and Gln927 were two common amino acid residues shown in the interactions, and Gln 927 displayed hydrogen-bonding ( Figure 2B-D). Table 2 illustrates that ZINC15675948 displayed a strong binding to c-MYC, revealed by a lowest binding energy (LBE) value of −9.91 kcal/mol and a predicted inhibition constant (pKi) of 0.05 µM. In comparison, the binding affinities of known c-MYC inhibitors were weak (LBE value of 10058-F4: −4.92 kcal/mol; LBE value of 10074-A4: −6.42 ± 0.01 kcal/mol; LBE value of 10074-G5: −6.96 ± 0.01 kcal/mol). To better understand the cross-resistance of CCRF-CEM and MDA-MB-231-pcDNA3 cells toward ZINC15675948, molecular docking of ZINC15675948 towards P-gp and BCRP was carried out. Table 3 revealed that ZINC15675948 showed high binding affinities to P-gp and BCRP (LBE values: −10.55 ± 0.24 kcal/mol and −11.49 kcal/mol, respectively). The binding site of ZINC15675948 to P-gp was similar to doxorubicin ( Figure 3A). They shared the same amino acid with Trp232 and Gln990 on P-gp ( Figure 3B,C). ZINC15675948, likewise doxorubicin, bound to the substrate-binding site on P-gp. However, the binding site of ZINC15675948 was similar to both doxorubicin and Ko143 ( Figure 3E-H). We further performed a doxorubicin uptake assay (Section 2.11) to pinpoint whether ZINC15675948 is an inhibitor or substrate of P-gp and BCRP.

Microscale Thermophoresis
To confirm the in silico binding of ZINC15675948 to c-MYC with an in vitro assay, we applied microscale thermophoresis (MST). As shown in Figure 2F, the measurement of the concentration-dependent fluorescence signals revealed an interaction between the fluorescently labeled c-MYC protein and ZINC15675948. ZINC15675948 bound to c-MYC with a K d of 1.08 ± 0.1 µM.    determined by a c-MYC reporter assay. Statistical significance (* p ≤ 0.05) was compared to DMSO (negative control). All experiments were performed three times independently.

c-MYC Reporter Assay
To determine whether the c-MYC activity could be diminished by ZINC15675948 binding, we performed c-MYC reporter assays in HEK293 cells with a transfected c-MYCluciferase reporter construct. Notably, the c-MYC activity was suppressed by ZINC15675948 in a concentration-dependent manner. A significant inhibition was observed at a concentration of 320 nM ( Figure 2G). Surprisingly, 10058-F4 as a positive control only showed a slight inhibition. Therefore, ZINC15675948 indeed inhibited c-MYC activity as consistently also demonstrated with molecular docking, microscale thermophoresis, and as illustrated below by western blotting and qRT-PCR.

Gene Expression Profile of Cell Lines Using Microarray Analyses
The gene expression measured using mRNA microarray hybridization was filtered by Chipster software (version 3.16.3). A total of 329 and 314 genes were significantly deregulated in CCRF-CEM cells and MDA-MB-231-pcDNA3 cells, respectively, compared with their untreated samples (Supplementary Tables S1 and S2). The deregulated genes were further analyzed by Ingenuity Pathway Analysis (IPA) to predict canonical pathways, networks, and cellular functions and diseases affected by ZIN15675948. Here, we did not observe any potentially impacted canonical pathways. Interestingly, several affected cellular functions were commonly revealed in the two cell lines, including "cell death and survival", "cell cycle", "cellular growth and proliferation", as well as "DNA replication, recombination and repair" (Figure 4). "Cancer" and "hematological disease" were affected correspondingly by ZINC15675948 in CCRF-CEM cells. MDA-MB-231-pcDNA3 cells also showed "cancer" as an important affected disease. Therefore, we further investigated the roles of the cell cycle, apoptosis, autophagy, and DNA damage to verify these microarraybased results by independent other methods.
We further accessed the networks under different cellular functions to unravel the genes that were involved. The IPA-based comparison of untreated and ZINC15675948treated CCRF-CEM cells indeed also revealed a downregulation of the c-MYC gene as illustrated in Figure 5A,B. These networks are related to cell cycle and cell death, implying that they were downstream and affected by the c-MYC gene. However, the c-MYC gene did not appear in the gene expression profiles of MDA-MB-231-pcDNA3 cells. The cell death network revealed that MCL-1 and BAD were downregulated accompanied by an upregulation of SQSTM1 (p62) ( Figure 5C). Furthermore, Figure 5D shows that the ELL2 gene was upregulated, which is involved in the proteasomal degradation of c-MYC [37,38]. The tumor suppressor TP53 was upregulated, which may be a consequence of c-MYC downregulation as we will discuss below.

Quantitative Reverse Transcription PCR (qRT-PCR)
The technical verifications of the results obtained from microarray hybridization were conducted with the top deregulated genes. Two upregulated and downregulated genes (RAD21, HMGCS1, PGK1, and ATP5MF) in CCRF-CEM cells, and two upregulated and downregulated genes (HSPD1, CITED2, H4C3, and DHFR) in MDA-MB-231-pcDNA3 cells were subjected to perform qRT-PCR ( Figure 6A). The Pearson correlation coefficients were calculated between the determined fold change of microarray hybridization and qRT-PCR data. As shown in Figure 6B,C, the r value was 0.98 in CCRF-CEM cells and 0.97 in MDA-MB-231-pcDNA3 cells, confirming a high degree of concordance between these two different methods. Furthermore, the expression level of c-MYC was determined by qRT-PCR ( Figure 6D). Treatment of ZINC15675948 (IC 50 ) significantly downregulated c-MYC (p ≤ 0.05) in both leukemia and breast cancer cell lines, indicating that ZINC15675948 inhibited c-MYC at the gene expression level.
Molecules 2023, 28, x FOR PEER REVIEW 11 of 30 2.6. Quantitative Reverse Transcription PCR (qRT-PCR) The technical verifications of the results obtained from microarray hybridization were conducted with the top deregulated genes. Two upregulated and downregulated genes (RAD21, HMGCS1, PGK1, and ATP5MF) in CCRF-CEM cells, and two upregulated and downregulated genes (HSPD1, CITED2, H4C3, and DHFR) in MDA-MB-231-pcDNA3 cells were subjected to perform qRT-PCR ( Figure 6A). The Pearson correlation coefficients were calculated between the determined fold change of microarray hybridization and qRT-PCR data. As shown in Figure 6B,C, the r value was 0.98 in CCRF-CEM cells and 0.97 in MDA-MB-231-pcDNA3 cells, confirming a high degree of concordance between these two different methods. Furthermore, the expression level of c-MYC was determined by qRT-PCR ( Figure 6D). Treatment of ZINC15675948 (IC50) significantly downregulated c-MYC (p ≤ 0.05) in both leukemia and breast cancer cell lines, indicating that ZINC15675948 inhibited c-MYC at the gene expression level.

Single Cell Gel Electrophoresis (Alkaline Comet Assay)
As "DNA replication, recombination and repair" appeared in our IPA analysis as "cellular functions" affected by ZINC15675948, we performed alkaline comet assays to detect DNA damage at the level of single cells. Representative images are shown in Figure  7. Compared with non-damaged control cells (DMSO), there was an increase in ZINC15675948-induced comet tails in both CCRF-CEM and MDA-MB-231-pcDNA3 cells,

Single Cell Gel Electrophoresis (Alkaline Comet Assay)
As "DNA replication, recombination and repair" appeared in our IPA analysis as "cellular functions" affected by ZINC15675948, we performed alkaline comet assays to detect DNA damage at the level of single cells. Representative images are shown in Figure 7. Compared with non-damaged control cells (DMSO), there was an increase in ZINC15675948-induced comet tails in both CCRF-CEM and MDA-MB-231-pcDNA3 cells, suggesting that DNA was indeed damaged. H 2 O 2 as a positive control also led to clearly visible comet tails. The analysis of tails of each 50 cells revealed that ZINC15675948 induced DNA damage in both cell lines in a concentration-dependent manner.
Molecules 2023, 28, x FOR PEER REVIEW 12 suggesting that DNA was indeed damaged. H2O2 as a positive control also led to cl visible comet tails. The analysis of tails of each 50 cells revealed that ZINC1567594 duced DNA damage in both cell lines in a concentration-dependent manner.

Cell Cycle Arrest
The IPA analysis of the microarray data also revealed that cell cycle progression disturbed by ZINC15675948 in both CCRF-CEM and MDA-MB-231-pcDNA3 cells. T fore, we investigated the cell cycle distribution by flow cytometry. Figure 8 shows ZINC15675948 significantly increased the fraction of CCRF-CEM cells in the G2/M p after treatment for 72 h, which was in a range of 28.25-33.3% at different concentra compared with DMSO (15.0%). Whereas the S phase fraction of MDA-MB-231-pcD increased after 24 h. At the highest concentration (2 × IC50), S phase arrest reached 2 (p < 0.01) compared with DMSO as mock control, which was close to the positive con cisplatin (36.2%) in MDA-MB-231-pcDNA3 cells. Hence, ZINC15675948 arrested the kemia cells in the G2/M phase and the breast cancer cells in the S phase of the cell cy

Cell Cycle Arrest
The IPA analysis of the microarray data also revealed that cell cycle progression was disturbed by ZINC15675948 in both CCRF-CEM and MDA-MB-231-pcDNA3 cells. Therefore, we investigated the cell cycle distribution by flow cytometry. The results were represented as mean ± SD from three independent measurements. Statistical significance was analyzed using Student's t-test, * p < 0.05, ** p < 0.01 compared with DMSO.

Detection of Apoptosis
As "cell death and survival" was revealed as cellular function in the microarray analysis, we measured apoptosis induction by ZINC15675948 using flow cytometry and annexin V and PI double staining in CCRF-CEM cells and with F2N12S and SYTOX ® AADvancedTM dead cell staining in MDA-MB-231-pcDNA3 cells. As shown in Figure  9A, ZINC15675948 induced late apoptosis in CCRF-CEM cells in a time-and concentration-dependent manner. The late apoptotic cells remarkably increased to 41.4% at a concentration of 4 × IC50 after 72 h, which was similar to the effect of 5 µM vincristine (66%). In parallel, ZINC15675948 also induced apoptosis in MDA-MB-231-pcDNA3 cells after 48 h ( Figure 9B). Especially, treatment with 4 × IC50 resulted in 14.7% apoptotic cells compared to DMSO (5.9%). Therefore, ZINC15675948 significantly induced apoptosis in both cell lines. The results were represented as mean ± SD from three independent measurements. Statistical significance was analyzed using Student's t-test, * p < 0.05, ** p < 0.01 compared with DMSO.

Detection of Apoptosis
As "cell death and survival" was revealed as cellular function in the microarray analysis, we measured apoptosis induction by ZINC15675948 using flow cytometry and annexin V and PI double staining in CCRF-CEM cells and with F2N12S and SYTOX ® AADvancedTM dead cell staining in MDA-MB-231-pcDNA3 cells. As shown in Figure 9A, ZINC15675948 induced late apoptosis in CCRF-CEM cells in a time-and concentrationdependent manner. The late apoptotic cells remarkably increased to 41.4% at a concentration of 4 × IC 50 after 72 h, which was similar to the effect of 5 µM vincristine (66%). In parallel, ZINC15675948 also induced apoptosis in MDA-MB-231-pcDNA3 cells after 48 h ( Figure 9B). Especially, treatment with 4 × IC 50 resulted in 14.7% apoptotic cells compared to DMSO (5.9%). Therefore, ZINC15675948 significantly induced apoptosis in both cell lines.  ( Figure 10B). Therefore, ZINC15675948 induced autophagy in CCRF-CEM cells but not in MDA-MB-231-pcDNA3 cells.
Furthermore, we verified biomarkers of the G2/M and S phases of the cell cycle that were arrested by ZINC15675948 in CCRF-CEM cells and MDA-MB-231-pcDNA3 cells, respectively. Indeed, the expression of CDK1 (indicative for G2/M phase arrest) was downregulated in CCRF-CEM cells, and the expression of phosphorylated CDK2 (Thr160) (indicative for S phase arrest) decreased in MDA-MB-231-pcDNA3 cells. To determine whether autophagy contributes to the "cell death and survival" pathway predicted by IPA, we examined two autophagy biomarkers, Beclin and p62. As shown in Figure 10A, the expression of p62 was significantly downregulated by ZINC15675948 at a concentration of 4 × IC 50 in CCRF-CEM cells (p = 0.02) with a slightly raised expression of Beclin, while the expression of p62 was significantly upregulated in MDA-MB-231-pcDNA3 cells, which was consistent with p62 gene upregulation found with microarray hybridization. The expression of Beclin was accordingly downregulated ( Figure 10B). Therefore, ZINC15675948 induced autophagy in CCRF-CEM cells but not in MDA-MB-231-pcDNA3 cells.
Furthermore, we verified biomarkers of the G2/M and S phases of the cell cycle that were arrested by ZINC15675948 in CCRF-CEM cells and MDA-MB-231-pcDNA3 cells, respectively. Indeed, the expression of CDK1 (indicative for G2/M phase arrest) was downregulated in CCRF-CEM cells, and the expression of phosphorylated CDK2 (Thr160) (indicative for S phase arrest) decreased in MDA-MB-231-pcDNA3 cells.

Doxorubicin Uptake Assay
Since multidrug-resistant cell lines were cross-resistant to ZINC15675948 in the growth inhibition assay, the uptake assay of doxorubicin using flow cytometry was further carried out in P-gp-overexpressing (CEM/ADR5000) and BCRP-overexpressing (MDA-MB-BCRP) cells to verify the functional role of ZINC15675948 with P-gp and BCRP. As shown in Figure 11A, treatment with positive control verapamil significantly enhanced the ability of CEM/ADR5000 cells to uptake doxorubicin by 3.07-fold (p = 0.02). In comparison, ZINC15675948-treated ADR/CEM5000 cells with three different concentrations did not observe doxorubicin accumulation. In MDA-MB-BCRP cells ( Figure 11B), Ko143 as a BCRP inhibitor also remarkably increased doxorubicin uptake (p = 0.003), while there was no evidence of ZINC15675948 treatments increasing doxorubicin accumulation in MDA-MB-BCRP cells. Therefore, these results confirmed that ZINC15675948 is a substrate of P-gp and BCRP.

Discussion
Over a time period from 1946 to 2019, natural products and their derivatives accounted for a majority of all anticancer drugs [39], suggesting the eminent importance of natural products in drug development. Seeking novel c-MYC inhibitors represents an urgent demand as c-MYC broadly initiates cancer development. A c-MYC deregulation is frequently involved in the carcinogenesis of lymphoblastic leukemia and triple-negative breast cancer [12,33]. Recently, compounds harboring an oxadiazole scaffold efficiently blocked c-MYC activity. The well-known c-MYC inhibitors, 10074-G5, and 10074-A4, are excellent examples [40,41]. The aim of this study was to investigate a natural product derivate of 1,2,4-oxadiazole, ZINC15675948, as a potential c-MYC inhibitor and its modes of action in leukemia and breast cancer cells.
We first demonstrated that ZINC15675948 exhibited profound growth inhibitory activity against leukemia and breast cancer cell lines at nanomolar ranges. Another 1,2,4oxadiazole derivative with reported strong cytotoxic activities toward leukemia and breast cancer cell lines supported our results [18]. Since CEM/ADR5000 and MDA-MB-BCRP cell lines were cross-resistant to ZINC15675948, we selected the drug-sensitive cell lines CCRF-CEM and MDA-MB-231-pcDNA3 to analyze the molecular mechanisms of this compound.
Next, we addressed the question of whether ZINC15675948 can inhibit c-MYC activity or c-MYC expression. MYC was previously considered as undruggable [42]. Recently, strategies to directly or indirectly inhibit MYC expression with small molecules showed promising results, such as targeting the MYC-MAX interface or disrupting MYC-MAX DNA binding, inducing MYC degradation by ubiquitination protease system, and targeting upstream regulators of MYC [43][44][45]. In our study, molecular docking revealed that ZINC15675948 displayed even better LBE to c-MYC than the three positive drugs. The binding of ZINC15675948 to c-MYC was close to the MYC-MAX interaction site, similar to 10058-F4, 10074-A4, and 10074-G5, which have been demonstrated to bind to the MYC monomer and inhibit its ability to interact with MAX [46]. It is known that one of the most important properties of oxadiazoles is that they act as hydrogen bond acceptors and thus are stable in water-based medium [19]. Our molecular docking definitely showed that ZINC15675948 interacted with Gln927 as a hydrogen bond. Microscale thermophoresis was correlated well with in silico binding, which indicated that ZINC15675948 bound to c-MYC with a K d value of 1.08 ± 0.1 µM. In comparison, 10074-G5 showed the best lowest binding energy of all known c-MYC inhibitors. It bound to MYC with a K d value of 4.4 µM as reported [47]. Meanwhile, the c-MYC reporter assay confirmed that ZINC15675948 showed a better inhibition activity to c-MYC than 10058-F4. Furthermore, the expression of c-MYC was downregulated by ZINC15675948 determined using western blotting and qRT-PCR. Therefore, ZINC15675948 can inhibit c-MYC expression. ZINC15675948 may directly target c-MYC by interfering with the MYC-MAX dimer and thereby affecting c-MYC activity.
In our microarray analysis, c-MYC expression was only downregulated in CCRF-CEM cells but not in MDA-MB-231-pcDNA3 cells by ZINC15675948. However, ZINC15675948 upregulated ELL2 in MDA-MB-231-pcNDA3 cells. Besides transcriptional and post-transcriptional mechanisms, the stability and activity of c-MYC are regulated by a series of post-translational modifications, including phosphorylation, methylation, acetylation, glycosylation, SUMOylation, proline isomerization, and ubiquitination. The ubiquitination-proteasome pathway is one of the prime mechanisms to degrade c-MYC [10,48]. Several E3 ubiquitin ligases were identified to positively or negatively regulate c-MYC inhibition through degradation [49]. The product of the ELL (eleven-nineteen lysine-rich leukemia) gene is a transcription elongation factor that accelerates the transcription elongation process of RNA polymerase II. ELL1 (ELL) is the first identified potential c-MYC suppressor in the ELL family. ELL1 functions as an E3 ubiquitin ligase through a CEYLH region in the C-terminus of the protein and targets c-MYC for proteasomal degradation, which thereby inhibits c-MYCdependent transcriptional activity and cell proliferation [38]. This CEYLH region shares homology within the ELL family suggesting that ELL2 could also have a similar role as ELL1 in the ubiquitin-proteasomal degradation of c-MYC [37]. Therefore, we assume that ELL2 ubiquitinated c-MYC protein in MDA-MB-231-pcDNA3 cells, explaining why c-MYC mRNA was not downregulated in the microarray analysis. It is worth pointing out that the c-MYC activity has been indeed reported to be regulated by ubiquitination [49]. Hence, ZINC15675948 may inhibit c-MYC by ubiquitination in MDA-MB-231-pcDNA3 cells.
c-MYC plays a critical role in regulating cell growth, differentiation, cell cycle, apoptosis, metabolism, angiogenesis, DNA repair, immune response, protein translocation, and stem cell formation [6]. However, c-MYC activation is not sufficient alone to cause carcinogenesis. Key mediators of c-MYC-driven cellular processes such as cell cycle checkpoints, apoptosis regulators, as well as tumor suppressor proteins are also required for tumorigenesis. Therefore, synthetic lethal interactions may provide additional opportunities to indirectly inhibit c-MYC [9]. Our microarray hybridization data showed that the "cell cycle", "cell death and survival" and "DNA damage" functions were affected both in leukemia and breast cancer cells by ZINC15675948.
DNA damage is an important mechanism of cancer chemotherapeutics that is associated with cell cycle arrest and apoptosis to limit cancer progression [50]. Using single cell gel electrophoresis (alkaline comet assay), we observed that ZINC15675948 remarkably caused DNA damage as seen in increased comet tails upon treatment. In line with this result, ZINC15675948 arrested CCRF-CEM cells in the G2/M phase and MDA-MB-231-pcDNA3 cells in the S phase of the cell cycle as determined by flow cytometry. These results were further confirmed by a downregulated expression of CDK1 or p-CDK2 (Thr160) using western blotting, respectively. C-MYC overexpression promotes cancer cells to reenter the cell cycle by activating cyclin-dependent kinases (CDKs) and cyclins [9]. CDK1 represents a core component in mammalian cell division. Inhibition of CDK1 induced c-MYC-dependent apoptosis in lymphoblastoid cell lines [51]. CDK2 is activated in the late G1 phase and continues into the S phase. Activation of CDK2 requires phosphorylation of Thr160 [52]. Hence, we conclude that ZINC15675948 treatment resulted in DNA damage and cell cycle arrest as downstream events of c-MYC inhibition.
ZINC15675948 significantly induced apoptosis in both leukemia and breast cancer cell lines. Induction of cell death in CCRF-CEM cells was also indicated as a c-MYC-dependent event by our IPA analysis due to the downregulated c-MYC expression. Although c-MYC did not appear in the microarray-based dataset of MDA-MB-231-pcDNA3 cells, TP53 was upregulated. Evidence for crosstalk of c-MYC and p53 networks showed that c-MYCinduced apoptosis is dependent on p53 stabilization and activation [53]. In addition, the induction of apoptosis in MDA-MB-231-pcDNA3 cells correlated well with the microarraybased IPA prediction that MCL-1 and BAD (BCL-xL) were downregulated. MCL-1 and BCL-xL are members of the anti-apoptotic B-cell lymphoma 2 (BCL-2) family. Overexpression of anti-apoptotic genes implies an inhibition of apoptosis and leads to cancer development. A recent study revealed that c-MYC-driven human Burkitt lymphoma cells were eradicated by targeting MCL-1 [54]. Therefore, these data favor the view that ZINC15675948-driven inhibition of c-MYC and related signals caused apoptosis.
In addition to apoptosis, autophagy is a type II programmed cell death that might also be relevant for the modes of action of ZINC15675948. Autophagy is a mechanism that transports damaged proteins or organelles to lysosomes for degradation, playing important roles in maintaining energy homeostasis and ensuring cellular materials' quality control [55]. If autophagy is induced, Beclin is recruited to organize the pre-autophagosomal structure that is necessary for autophagosome formation [56]. P62 (SQSTM1) transports ubiquitinated cargoes for autophagic degradation. The activation of autophagy declines the expression of p62 [57]. In this study, we demonstrated that the induction of autophagy by ZINC15675948 was cell line-dependent. ZINC15675948 induced autophagy in CCRF-CEM cells as evidenced by a significantly downregulated expression of p62 and an upregulation of Beclin. Induction of autophagy could result from the suppression of c-MYC. These results are supported by our recent report from our group that artemisinin-induced autophagy followed upon c-MYC inhibition in CCRF-CEM cells [58]. Autophagy induction was not observed in MDA-MB-231-pcDNA3 cells, while the upregulated expression of p62 was correlated with that of microarray data.
ATP-binding cassette (ABC) transporters are widely expressed in various tissues protecting them from endogenous molecules and xenobiotics. ABC transporters are involved in the efflux of cytotoxic drugs from cancer cells leading to the multidrug resistance (MDR) phenotype [59] with broad cross-resistance profiles to diverse chemotherapeutics without structural or pharmacological commonalities, including Vinca alkaloids, taxanes, anthracyclines, epipodophyllotoxins, and others [60]. Human P-glycoprotein and BCRP are well-known ABC transporters mediating MDR. P-gp is encoded by the ABCB1/MDR1 gene and is expressed in the normal intestine, kidneys, liver, placenta, brain, and adrenal glands. BCRP is encoded by the ABCG2 gene and is also expressed in many tissues, including breast, colon, liver, bile, and in brain capillaries [61][62][63]. P-gp and BCRP confer anticancer drug resistance and were clinically found to correlate with poor diagnosis and survival of cancer patients [64,65]. Cytotoxic agents displaying degrees of cross-resistance above two can be considered as moderate substrates of ABC transporters, while high degrees of resistance (>1000) can also be observed in MDR cells [66,67]. P-gp overexpressing CEM/ADR5000 cells were 8.37 times more cross-resistant to ZINC15675948, while BCRP overexpressing cells were nine times more resistant to ZINC15675948 compared to their sensitive cell lines, indicating that ZINC15675948 was moderately cross-resistant to P-g and BCRP. Using molecular docking, we observed that ZINC15675948 bound at the substrate binding site on P-gp, but its binding site on BCRP is ambiguous concerning the substrate or inhibitor pockets. In the further doxorubicin uptake assay, our data confirmed that ZINC15675948 failed to block the efflux. Therefore, there was no intracellular accumulation of doxorubicin, which can be taken as a further hint that ZINC15675948 is a substrate of P-gp and BCRP. This finding is consistent with our recent report [68]. Future structure-activity relationship (SAR) investigations of derivatives of ZINC15675948 can be carried out to search for more therapeutically promising P-gp and BCRP inhibitors.

Growth Inhibition Assay
The resazurin reduction assay was used to access the effect of growth inhibition of ZINC15675948 [74]. The non-fluorescent dye resazurin is metabolically converted to the fluorescent dye resorufin by living cells [75]. Briefly, the CCRF-CEM and CEM/ADR5000 suspension cells were seeded into 96-well plates (1 × 10 4 cells/well) and then directly treated with 10 concentrations of ZINC15675948 in a range of 0.3-100 µM, respectively, in a total volume of 200 µL. MDA-MB-231-pcDNA3 and MDA-MB-231-BCRP clone 23 cells were seeded into 96 well plates (5 × 10 3 cells/well) overnight and treated in the same series of concentrations as in leukemia cells with ZINC15675948 on the following day. After 72 h incubation, 20 µL of 0.01% resazurin (Promega, Germany) were added to each well and incubated for 4 h at 37 • C. The resazurin fluorescence was detected using an Infinite M2000 Pro TM plate reader (Tecan, Crailsheim, Germany) at Ex/Em = 550 nm/590 nm wavelength. Relative cell viability was calculated in comparison to the DMSO-treated control. The final concentration of DMSO was 0.5%. The growth inhibition was accessed according to the effectiveness in inhibiting cell proliferation by half and was expressed as half-maximal inhibitory concentration (IC 50 ) values. All IC 50 values were expressed as mean ± standard deviation (SD). This experiment was repeated three times independently with six wells for each concentration. The figures were analyzed using GraphPad Prism Software (version 9.0.2) (GraphPad Software Inc., San Diego, CA, USA).

Molecular Docking
In silico molecular docking was performed using the AutoDock 4.2.6 software (The Scripps Research Institute, CA, USA). The protocol was recently described by us [69]. The SDF format of ZINC15675948 was downloaded from ZINC 15 database (https://zinc15. docking.org/ accessed on 20 October 2019). ZINC15675948 was considered to perform molecular docking to P-glycoprotein (P-gp, MDR1, ABCB1) and the breast cancer resistance protein (BCRP, ABCG2), since growth inhibition of ZINC15675948 was determined in P-gpor BCRP-overexpressing multidrug-resistant cell lines. The 3D structures of c-MYC (PDB code: 1NKP), P-gp (PDB code: 6QEX) and BCRP (PDB code: 6FFC) were downloaded from the RCSB Protein Data Bank (http://www.rcsb.org/, accessed on 20 October 2019) as PDB files. The known c-MYC inhibitors 10058-F4, 10074-A4 and 10074-G5 were used as positive control drugs to compare their affinity and binding mode with ZINC15675948 [76]. C-MYC was covered with grid box for the defined docking mode. The center of the grid box was set at x = 67.083, y = 64.481, and z = 42.943 with a spacing of 0.397 and with a number of grid points of 68 in x, 64 in y, and 72 in z). Doxorubicin as a known substrate of P-gp and BCRP was used as a control drug. Verapamil and Ko143 are inhibitors of P-gp or BCRP, respectively. Both were also used as control drugs. The grid boxes of P-gp and BCRP were placed around the drug-binding site [74]. Hydrogens were added to each protein structure, missing atoms were checked. The Lamarckian GA (4.2) was applied as an algorithm for 250 runs and 25,000,000 energy evaluations for each cycle. Docking was performed three times independently and the predicted inhibition constants were obtained from the docking log files (dlg). VMD (Visual Molecular Dynamics) software (version 1.

Microscale Thermophoresis
Microscale thermophoresis (MST) was performed with human recombinant c-MYC protein and ZINC15675948. The protein was purchased from Abcam (Cambridge, UK) with a concentration of 0.5 mg/mL (ab169901). C-MYC was labeled with Monolith Protein Labeling Kit RED-NHS 2nd Generation (Nano Temper Technologies GmbH, Munich, Germany) according to the manufacturer's instructions. The final concentration of c-MYC after labeling was 595 nM. ZINC15675948 was diluted from 400 µM to a series of concentrations in assay buffer (50 mM Tris buffer (pH 7.4) containing 10 mM MgCl 2 , 150 mM NaCl, and 0.05% Tween-20). The labeled protein and diluted compounds were mixed (1:1). After 30 min incubation in the dark at room temperature, the fluorescence signal was measured on a Monolith NT.115 instrument (Nano Temper Technologies) with Monolith NT.115 standard capillaries. The MST with ZINC15675948 was performed with 70% LED power and 10% MST power. Fitting curved and dissociation constant (K d ) values were calculated with MO. Affinity Analysis software (version 2.2.4) (Nano Temper Technologies). The measurements were repeated three times independently.

c-MYC Reporter Assay
A signal MYC reporter assay (Qiagen, Germantown, MD, USA) was used to determine the impact of ZINC15675948 on c-MYC activity as we described recently [58]. Briefly, a c-MYC-luciferase reporter construct was transfected into human embryonic kidney HEK293 cells and incubated according to the manufacturer's instructions. Subsequently, cells were treated with two concentrations (4 × IC 50 in CCRF-CEM or MDA-MB-231-pcDNA3) of ZINC15675948 (32 nM and 320 nM) and DMSO (negative control) or of the known c-MYC inhibitor 10058-F4 (positive control) for 48 h. A Dual-glo ® Luciferase Reporter Assay System (Promega, Madison, WI, USA) was applied for the measurement of c-MYC promoter activity. Renilla and firefly luciferase luminescence were measured using an Infinite M2000 Pro TM plate reader (Tecan, Crailsheim, Germany).

Gene Expression Profiles
Total mRNA was isolated using the InviTrap ® Spin Universal RNA Mini Kit (Invitek Molecular, Berlin, Germany). CCRF-CEM (1 × 10 6 cells/well) and MDA-MB-231-pcDNA3 cells (5 × 10 5 cells/well) were treated with DMSO as a solvent control and ZINC15675948 for 24 h in duplicate. ZINC15675948 was applied with a concentration according to the IC 50 value (CCRF-CEM cells: 0.008 µM, MDA-MB-231-pcDNA3 cells: 0.08 µM). The RNA concentrations were determined using NanoDrop1000 (PEQLAB, Erlangen, Germany). Afterwards, quality control of total RNA, probe labeling, hybridization, scanning and data analysis of the samples were performed at the Genomics and Proteomics Core Facility of the German Cancer Research Center (DKFZ, Heidelberg, Germany). Affymetrix GeneChips ® with human Clariom TM S assays (Affymetrix, Santa Clara, CA, USA) were applied for microarray hybridizations as previously described in detail [34].

Pathway Analysis of Microarray Data
The Chipster software (http://chipster.csc.fi/, accessed on 20 October 2019) (The Finnish IT Center for Science CSC, Espoo, Finland) was used to filter a set of differentially expressed genes acquired from microarray hybridization [77]. The Empirical Bayes t-test (p < 0.05) was applied to access the deregulated genes between DMSO and ZINC15675948treated groups (accessed in July 2021). The filtered genes were analyzed with the Ingenuity Pathway Analysis (IPA) software (Qiagen, Redwood City, CA, USA) (content version 51963813) by the core analysis tool to determine the cellular functions and networks affected by drug treatment (accessed in August 2021).

Quantitative Real-Time Reverse Transcription PCR
The same total RNA samples (DMSO control and IC 50 ) used for the microarray analyses were also used for qRT-PCR experiments [78]. One microgram RNA was converted to cDNA using the LunaScript ® RT SuperMix Kit cDNA Synthesis Kit (New England Bio Labs, Darmstadt, Germany) according to the manufacturer's instructions. All PCR primers were designed using the NCBI Primer-BLAST (https://www.ncbi.nlm. nih.gov/tools/primer-blast/) website and purchased from Eurofins genomics (Ebersberg, Germany) (https://eurofinsgenomics.eu/en/dna-rna-oligonucleotides/optimisedapplication-oligos/pcr-primers/) (accessed in September 2021). The primer sequences are shown in Table 4. The GAPDH gene served as an internal control. The reaction mix-ture contained 4 µL master mix (5 × Hot Start Taq EveGreen ® qPCR Mix (no Rox) (Axon Labortechnik, Kaiserslautern, Germany), 1 µL forward or reversed primer (250 nM final concentration), 13 µL nuclease-free water (Thermo Fisher), and 1 µL cDNA converted from 300 ng RNA. The qRT-PCR was performed on CFX384 Touch Real-Time PCR Detection System (Bio-Rad, Munich, Germany) using the 384-well plate. The initial denaturation of qRT-PCR was at 95 • C for 10 min followed by 40 cycles including strand separation at 95 • C for 15 s, annealing at 57.5 • C for 40 s and extension at 72 • C for 1 min. CFX Manager Software (version 3.1) was used to generate the C q values. The fold-change of gene expression was calculated using comparative 2 −∆∆CT method as reported in the literature [79]. Specifically, for the DMSO control sample, the data of the target gene are presented as fold change in gene expression normalized to GAPDH gene. For the sample treat with the IC 50 of the compound, the evaluation of 2 −∆∆CT after normalization to GAPDH gene reveals the fold change of the target gene in gene expression relative to the untreated control.

Gene Symbol
Forward Primer Reverse Primer

Single Cell Gel Electrophoresis (Alkaline Comet Assay)
DNA damage was detected by alkaline comet assay using the OxiSelect™ Comet Assay Kit (Cell Biolabs/Biocat, Heidelberg, Germany) as described in [78]. The alkaline comet assay is a sensitive method for monitoring the migration of DNA fragments from nuclei under alkaline conditions. It detects DNA single-and double-strand breaks at a cellular level [50]. Briefly, CCRF-CEM cells (1 × 10 6 cells/well) were seeded into a 6-well plate and treated with different concentrations of ZINC15675948 (IC 50 , 2 × IC 50 and 4 × IC 50 ) and DMSO as a negative control for 24 h. MDA-MB-231-pcDNA3 cells (5 × 10 5 cells/well) were first seeded for 24 h to allow adhesion and then incubated with the same treatments for 24 h. Both cell lines were treated with H 2 O 2 (50 µM) as a positive control for 30 min [80]. Cells were harvested and centrifuged at 3000× g for 10 min and were suspended in 1 mL cold PBS. Next, 1 × 10 5 cells were counted and mixed with agarose at 37 • C at a ratio of 1:6 and then spread on a comet slide. The following steps were conducted in the dark. Slides were left at 4 • C for 30 min to solidify and then immersed in pre-chilled lysis solution (NaCl 14.6 g, EDTA solution 20 mL, 10× lysis solution, pH 10.0, fulfill to 100 mL with distilled water, stored at 4 • C) for 1 h at 4 • C. Afterwards, slides were immersed in pre-chilled alkaline electrophoresis solution buffer (NaOH 12 g, EDTA solution 2 mL, fulfill to 100 mL with distilled water, stored at 4 • C) for 40 min. Next, electrophoresis was performed for 20 min at 20 V with alkaline electrophoresis solution buffer. Subsequently, slides were transferred into pre-chilled distilled water for 2 × 5 min for washing, followed by immersion in 70% ethanol for 5 min. After slides were dry, Vista Green DNA dye was diluted in TE buffer (Tris 121.14 mg, EDTA 200 µL, pH 7.5, fulfill to 100 mL distilled water) at a ratio of 1:1000. Then 100 µL diluted Vista Green DNA dye was added to each sample. DNA damage was captured by EVOS digital inverted microscope (Life Technologies GmbH, Darmstadt, Germany). At least 50 comets of each treatment were analyzed by Imaged J software (version 1.53q) using a plugin OpenComet (National University of Singapore, Singapore). The tail DNA percentage was used as a parameter of DNA damage [81,82].

Cell Cycle Arrest
CCRF-CEM cells (1 × 10 6 cells/well) were seeded into a 6-well-plate and treated with ZINC15675948 at concentrations of IC 50 , 2 × IC 50 or 4 × IC 50 , DMSO (negative control), or vincristine (positive control, 0.3 µM) (University Hospital Pharmacy, Mainz, Germany) for 72 h. MDA-MB-231-pcDNA3 cells (3 × 10 5 cells/well) were seeded and on the second day treated with ZINC15675948 at concentrations of 0.5 × IC 50 , IC 50 or 2 × IC 50 , and DMSO (negative control) and cisplatin (positive control, 0.5 µM) (University Hospital Pharmacy) for 24 h. The cells were harvested and centrifuged with cold PBS (4 • C) twice (1500 rpm for 5 min). Ice-cold ethanol (80%) was used for fixation. Samples were kept at −20 • C at least for 24 h. Afterwards, cells were spun down from the ethanol by centrifugation at 4000 rpm for 10 min and washed twice with cold PBS. Before measurement, the cells were resuspended with 500 µL cold PBS containing 20 g/mL RNase (Roche Diagnostics, Mannheim, Germany) and incubated at room temperature for 30 min, followed by staining with 50 µg/mL propidium iodide (PI) (Sigma-Aldrich). After 15 min incubation in the dark at 4 • C, CCRF-CEM cells were analyzed on a BD Accuri TM C6 Flow Cytometer (Becton-Dickinson, Heidelberg, Germany). MDA-MB-231-pcDNA3 cells were analyzed on a BD LSRFortessa SORP (Becton Dickinson, Heidelberg, Germany). The cells were gated firstly using FSC-A/SSC-A gate in linear scale, 10 4 cells were recorded. Then doublets were removed using FL2-A/FL2-H gate also in linear scale. The DNA histogram was generated using FL2-A/histogram properties. All the experiments were repeated three times independently. The cell cycle distributions were analyzed by the FlowJo software (version 10.8.1) (Celeza, Olten, Switzerland) [34].

Detection of Apoptosis in Adherent Cells
MDA-MB-231-pcDNA3 adherent cells generally need trypsinization for cell-harvesting, which may result in false positive results by applying the Annexin V-FITC apoptosis kit. Therefore, the Violet Ratiometric Membrane Asymmetry Prob/Dead Cell Apoptosis Kit (Thermo Fisher Scientific, Darmstadt, Germany) was alternatively utilized to detect apoptosis in adherent cells [84]. Furthermore, 4'-N,N-diethylamino-6-(N,N,N dodecylmethylamino-sulfopropyl)-methyl-3-hydroxyflavone (F2N12S) is a novel fluorescent probe to monitor the lipid composition on the plasma membrane in the early stage of apoptosis. SYTOX ® AADvanced TM can pass through the cell membrane only in necrosis or late apoptosis. Briefly, 1 × 10 5 MDA-MB-231-pcDNA3 cells were seeded in a 6-well plate overnight and then treated with different concentrations (0.5 × IC 50, IC 50 , 2 × IC 50 or 4 × IC 50 ) of ZINC1565948, DMSO (negative control), or vincristine (positive control, 1 µM). After the incubation for 48 h, cells were detached using 500 µL Accutase (Thermo Fisher Scientific, Darmstadt, Germany) at room temperature and then washed with 1 mL cold Hank's balanced salt solution (HBSS, 4 • C) twice. Cells were suspended with 1 mL HBSS. Subsequently, 1 µL F2N12S solution was added to each sample at a final concentration of 200 nM. Another 1 µL SYTOX ® AADvanced TM dead cell stain solution was added at a final concentration of 1 µM. Samples were incubated at room temperature for 5 min and then analyzed using a BD LSRFortessa SORP (Becton Dickinson). The F2N12S was excited at 405 nm, and the emissions were collected with orange fluorescence (585/15 bandpass filter) and green fluorescence (530/30 bandpass filter). A ratio of the orange fluorescence channel to the green fluorescence channel was set as a derive parameter. The SYTOX ® AADvanced TM dead cell stain was excited at 488 nm, and the emission was collected with a 670/30 bandpass filter. The detections were repeated three times independently. The data were analyzed by the FlowJo software (version 10.8.1) (Celeza).

Protein Analyses by SDS-PAGE and Immunoblotting
CCRF-CEM cells (6,000,000 cells/flask) and MDA-MB-231-pcDNA3 cells (500,000 cells/well) were treated with varying concentrations of ZINC15675948 (IC 50 , 2 × IC 50 , or 4 × IC 50 ) or DMSO as a negative control. After 24 h incubation, the cells were harvested and washed with PBS. Then 100 µL of ice-cold M-PER Mammalian Protein Extraction Reagent (Thermo Fisher Scientific, Darmstadt, Germany) containing 1% Halt Protease Inhibitor Cocktail and phosphatase inhibitor (Thermo Fisher Scientific) were added to each sample. The cell lysis solutions were shaken for 30 min at 4 • C and centrifuged at 14,000× g for 15 min at 4 • C. The supernatants were harvested, and the protein concentrations were measured with NanoDrop1000 (PEQLAB, Erlangen, Germany).

Doxorubicin Uptake Assay
Doxorubicin is a substrate of P-gp. P-gp-overexpressing CEM/ADR5000 cells were seeded in 12-well plates (10 4 cells/well). Doxorubicin was obtained from the University Hospital Pharmacy (Mainz, Germany) and used at a concentration of 10 µM in all samples. In parallel, different concentrations of ZINC15675948 (IC 50 , 2 × IC 50 or 4 × IC 50 ) were combined with doxorubicin treatment. Verapamil (20 µM; Sigma-Aldrich) was applied as a positive control for P-gp inhibition. Unstained CEM/ADR5000 cells were used as a negative control. In comparison, non-P-gp-expressing CCRF-CEM cells were used as a positive control for doxorubicin uptake. After 3 h incubation at 37 • C, cells were centrifuged to discard the old medium and resuspended in 1 mL PBS. The doxorubicin fluorescence assay was carried out on a BD Accuri TM C6 Flow Cytometer (Becton-Dickinson) by blue laser at excitation 488 nm and emission 530 nm. In each sample, 3000 cells were counted. Dead cells and debris were removed by gating the cells in forward vs. side scatter, and forward-area vs. forward-height scatter, respectively. All experiments were performed in triplicate. The protocol has been described by us in [74].

Conclusions
In conclusion, ZINC15675948 displayed remarkable cytotoxicity in CCRF-CEM and MDA-MB-231-pcDNA3 cells. ZINC15675948 bound to the c-MYC/MAX interface and inhibited c-MYC activity and expression. The inhibition of c-MYC in MDA-MB-231-pcDNA3 cells involved ubiquitination. By means of microarray-based mRNA expression profiling, we verified that ZINC15675948 induced DNA damage and apoptosis in both cell lines as downstream effects of c-MYC inhibition. ZINC15675948 caused G2/M phase cell cycle arrest in CCRF-CEM cells and S phase arrest in MDA-MB-231-pcDNA3 cells. ZINC15675948 also induced autophagy in CCRF-CEM cells. Furthermore, ZINC15675948 was a substrate of two ABC transporters, P-gp and BCRP. These results illustrate that ZINC15675948 is a promising inhibitor of c-MYC worth further development as a novel anticancer drug.
Supplementary Materials: The following supporting information can be downloaded at: https://www. mdpi.com/article/10.3390/molecules28155658/s1, Table S1: Deregulated gene expression upon treatment of CCRF-CEM leukemia cells with the IC 50 concentration of ZINC15675948 for 24 h; Table S2

Data Availability Statement:
The authors declare that the data supporting the findings of this study are available within the paper.