Comprehensive Omics Analysis of a Novel Small-Molecule Inhibitor of Chemoresistant Oncogenic Signatures in Colorectal Cancer Cell with Antitumor Effects

Tumor recurrence from cancer stem cells (CSCs) and metastasis often occur post-treatment in colorectal cancer (CRC), leading to chemoresistance and resistance to targeted therapy. MYC is a transcription factor in the nuclei that modulates cell growth and development, and regulates immune response in an antitumor direction by mediating programmed death ligand 1 (PD-L1) and promoting CRC tumor recurrence after adjuvant chemotherapy. However, the molecular mechanism through which c-MYC maintains stemness and confers treatment resistance still remains elusive in CRC. In addition, recent reports demonstrated that CRC solid colon tumors expresses C-X-C motif chemokine ligand 8 (CXCL8). Expression of CXCL8 in CRC was reported to activate the expression of PD-L1 immune checkpoint through c-MYC, this ultimately induces chemoresistance in CRC. Accumulating studies have also demonstrated increased expression of CXCL8, matrix metalloproteinase 7 (MMP7), tissue inhibitor of metalloproteinase 1 (TIMP1), and epithelial-to-mesenchymal transition (EMT) components, in CRC tumors suggesting their potential collaboration to promote EMT and CSCs. TIMP1 is MMP-independent and regulates cell development and apoptosis in various cancer cell types, including CRC. Recent studies showed that TIMP1 cleaves CXCL8 on its chemoattractant, thereby influencing its mechanistic response to therapy. This therefore suggests crosstalk among the c-MYC/CXCL8/TIMP1 oncogenic signatures. In this study, we explored computer simulations through bioinformatics to identify and validate that the MYC/CXCL8/TIMP1 oncogenic signatures are overexpressed in CRC, Moreover, our docking results exhibited putative binding affinities of the above-mentioned oncogenes, with our novel small molecule, RV59, Finally, we demonstrated the anticancer activities of RV59 against NCI human CRC cancer cell lines both as single-dose and dose-dependent treatments, and also demonstrated the MYC/CXCL8/TIMP1 signaling pathway as a potential RV59 drug target.

Furthermore, study evidence from a microarray analysis also outlined significant differences in several oncogenes between tumor tissues and normal tissues. Interestingly, CXCL8, MMP7, tissue inhibitor of metalloproteinase 1 (TIMP1), and components of the epithelial-to-mesenchymal transition (EMT) were found to be highly expressed in CRC tumors [33]. This suggests that CXCL8/MMP7/TIMP7 oncogenes may collaborate and promote the EMT and colon stem cells [34,35]. TIMP1 is matrix metalloproteinase (MMP)independent and regulates cell development and apoptosis in various cancer cell types, including CRC [36]. Recent studies showed that TIMP1 cleaves to CXCL8 on its chemoattractant, thereby influencing its mechanisms and responses to therapy [37]. TIMP1 was also reported to participate in various cell functions including proliferation and survival, leading to reduced sensitivity to chemotherapy in colon cancer [37][38][39]. Increased expression of TIPM1 was shown to be associated with poor clinical outcomes in CRC patients compared to normal samples [40]. This therefore suggests crosstalk among c-MYC/CXCL8/TIMP1 oncogenic signatures. In the current study, we used computation studies through a bioinformatics analysis to identify and validate expressions of the c-MYC/CXCL8/TIMP1 signaling pathway in CRC, and used in silico molecular docking to evaluate potential interactions of the RV59 with MYC/CXCL8/TIMP1 signaling. RV59 is a small molecule, a derivative of EGFR inhibitor osimertinib derivative [41,42], which was recently synthesized in our laboratory. The anticancer effects of the RV59 was evaluated using the US National Cancer Institute (NCI)-60 colon cancer cells to single-dose and dose-dependent treatments with RV59 [43].

Identifying Molecular Targets and Therapeutic Classes of RV59
To predict potential drug targets of RV59, we used computer-based Prediction of Biological Activity Spectra (PASS) (http://way2drug.com/PassOnline/) web resources to predict the spectrum of interactions for known protein kinase inhibitors [44]. Additionally, we applied the Developmental Therapeutics Program (DTP)-COMPARE, a public datamining website, to predict molecular targets and detect if RV59 retained activity similar to NCI synthetic compounds and standard agents. Herein, we used 50% growth inhibition (GI 50 ), which is the IC 50 as an endpoint, and National Safety Code (NSC) number (763967) as a delimiter [45].

Gene Expression Microarray Data Extraction
In total, three gene expression profiles (GSE41328, GSE44861, and GSE74602) were extracted from the Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/). The GEO database, a free public genomic database, stores array data and sequence data. GSE74602 included 10 CRC cancer samples and 10 matched adjacent normal samples [46]; GSE44861 included 56 colon cancer tissues and 55 adjacent noncancerous tissues [47]; and GSE74602 contained 30 CRC carcinoma tissues and 30 adjacent normal cancer tissues [48].

Validation of MYC/CXCL8/TIMP1 Expression Levels in CRC
Expression levels of the MYC/CXCL8/TIMP1 oncogenes were analyzed with UAL-CAN (http://ualcan.path.uab.edu/), an open-access public online tool for analysis of The Cancer Genome Atlas (TCGA) [49]. Expression levels of MYC/CXCL8/TIMP1 in CRC samples (red) were compared to adjacent normal samples (blue), with p < 0.05 indicating statistical significance. Furthermore, we used GECO, a gene expression correlation analytical tool, which distinguishes two expression datasets into positive and negative correlations [50], with positive Pearson correlation coefficients and p < 0.05 as statistically significant.

Protein-Protein Interaction (PPI) Network, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Analyses
Protein interactions were analyzed using the STRING tool (https://string-db.org/) to construct a PPI clustering network [51]. A confidence score >0.7 was considered most significant. The GeneMenia tool (https://genemania.org/) and cytoscape software (version 3.8.2) were used to build gene interactions and PPI networks respectively, and from analysis of the results, the interactive networks were based on gene co-expression, co-localization, genetic interactions, and various pathways involved within the network. The database for annotation, visualization, and integrated discovery (DAVID), (https://david.ncifcrf.gov/.jsp), was used to analyze enriched GO including biological processes and molecular functions involved, with the criterion set to p < 0.05.

Interpretation of Gene Co-Expression in MYC/CXCL8/TIMP1 Genes Network
Interpretation of gene expression network was analyzed using the Network Analyst 3.0 tool (https://www.networkanalyst.ca/) [52], a comprehensive analytical visual platform which integrates PPI networks and gene co-occurrence networks and interprets gene expression networks. Herein, we used Enrichment Map, a sub-tool of network analyst, with a Bipartite view, to determine the enrichment of coexpressed genes in KEGG pathways [53].

Correlation Analysis of MYC/CXCL8/TIMP1 Expressions and Tumor Infiltration Levels
Correlations between MYC/CXCL8/TIMP1 expressions and tumor infiltration levels were analyzed with the Tumor Immune Estimation Resource (TIMER) (https://cistrome. shinyapps.io/timer/), an online computational tool used to analyze the nature of tumor immune interactions across a variety of cancers [54]. Herein, we determined correlations of MYC, CXCL8, and TIMP1 with a set of gene markers of immune infiltration cells including CD8+ T cells and macrophages (with p < 0.05). Furthermore, we applied the SCNA statistical module, a sub-tool of TIMER to identify the abundances of tumor infiltrates including CD8+ T cells, macrophages, and dendritic cells (DCs) in colon cancer. The infiltration level was compared to the normal level using a two-sided Wilcoxon rank-sum test.

Drug Sensitivity Analysis of MYC/CXCL8/TIMP1 Oncogenes
To determine the correlation between MYC/CXCL8/TIMP1 oncogenes and drug sensitivity of the genomics of drug sensitivity in cancer (GDSC) top 30 drugs in pancancer, we used the Gene Set Cancer Analysis (GSCA), a web-based tool used to analyze differentially expressed genes (DEGs) and correlation to drug sensitivity [55]. All the drugs approved by the Food and Drug Administration (FDA) were displayed in this analysis.

Molecular Docking of Protein-Ligand Interactions
In silico molecular docking is a technique used to predict protein-ligand interactions, and this includes detection of the position and orientation of the ligand [56]. In order to evaluate the strength of interactions of RV59 with predicted and selected target genes from PASS and the DTP-compare algorithm, we performed molecular docking of RV59 with the MYC/CXCL8/TIMP1 oncogenes. Crystal structures of MYC (PDB: 6G6K), CXCL8 (PDB:5D14), and TIMP1 (1UEA), were downloaded from the Protein Data Bank (PDB). The 3D structure of RV59 was constructed with the Avogadro molecular visualization tool [57], while the standard inhibitor 3D structures for vorinostat (CID:5311), with a molecular weight (MW: 264.32 g/mol) with molecular formula (C 14 H 2 ON 2 O 3 ), repaxicin (CID:9838712), with a Mw of (283.39 g/mol and molecular formula (C 14 H 21 NO 3 ) and ilomastat (CID: 132519), with a Mw of (388.5 g/mol and molecular formula (C 20 H 28 N 4 O 4 ), were all retrieved from PubChem as SDF files. For further processing, we used PyMol software (https://pymol.org/2/) to visualize the ligands and convert them into PDB file format, and these files were subsequently converted into PDBQT format using the computer-based software, autodock (http://autodock.scripps.edu/resources/adt). For visualization and interpretation, docking results were analyzed using Discovery Studio [58].

Statistical Analysis
Pearson's correlations were used to assess correlations of MYC/CXCL8/TIMP1 expressions in CRC cancer types. The statistical significance of differentially expressed genes (DEGs) was evaluated using the Wilcoxon test. * p < 0.05 was accepted as being statistically significant.

MYC/CXCL8/TIMP1 Oncogenes Are Potential Drug Targets for RV59
We explored computer-based PASS and DTP-COMPARE, a public data-mining website to predict molecular targets and detect if RV59 retained activities similar to NCI synthetic compounds and standard agents. Similarities between our test compound and correlated drugs are presented as "Pearson's Correlation Coefficient" values, where a value of +1 indicates perfect positive correlation and a value of -1 indicates negative similar correlation. Herein, we used GI 50 as an endpoint and the NSC number as a delimiter [45], and identified several RV59 druggable genes. Interestingly, MYC, CXCL8, and TIMP1 were among the genes displayed by the prediction tools (Table 1). In addition to these finding, the activities of RV59 on the identified genes, shown by the PASS tool, were classified according to "probability to be active" (Pa) as opposed to "probability to be inactive" (Pi) values, based on the analyzed results (all Pa values were greater than Pi values), indicating the inhibitory activities of the compound on MYC and CXCL8, among other identified activities (Table 2). r, Pearson's correlation coefficient values range from −1 to +1 (the higher the positive number, the more positive correlation); CCLC, common cell lines count.

RV59 Passed the Required Drug-Likeness Criteria
RV59 small molecule, is a derivative of EGFR inhibitor osimertinib (CID: 71496458), which was recently synthesized in our laboratory. Using swissADME, a free web tool to evaluate pharmacokinetics, drug-likeness, and medicinal chemical friendliness of small molecules [47], we assessed the absorption, distribution, metabolism, and excretion (ADME) of RV59. Bioavailability radar displaying all six physical properties of RV59 showed that the compound successfully met the minimum requirements of drug likeness ( Figure 1). The criteria are based on the Mw of a compound (Mw recommended value ≤ 500 g/mol), flexibility (number of rotations: recommended value ≤ 10), solubility (log S (ESOL) recommended value 0-6), saturation (fraction Csp3 recommended value ≤ 1), polarity (TPSA recommended value ≤ 140 Å 2 ), and lipophilicity (XLOGP3 recommended value −0.7-5). Moreover, pharmacokinetics (PKs), drug-likeness, and medicinal chemical properties of RV59 (Table 3) indicated that RV59 had good synthetic accessibility of 2.95; this was evaluated according to the range from 1 (easy to make) to 10 (difficult to make). The compound also passed the criteria for Ghose, Veber (GSK), Egan (Pharmacia), and muegge and the Lipinski (Pfizer) rule-of-five of drug likeness and drug discovery. The bioavailability of the compound based on gastrointestinal absorptivity (GIA) indicated a score of 0.55 (55%) which indicates acceptable PK properties.   Pharmacokinetics displayed high gastrointestinal absorptivity (GIA) and low blood-brain barrier (BBB) permeability; drug-likeness, the compound passed all the criteria for Ghose, Veber (GSK), Egan (Pharmacia), and muegge, and the Lipinski (Pfizer) rule-of-five of drug likeness and drug discovery. The bioavailability of the compound based on GIA indicated a score of 0.55 (55%) and  Table 3 = 0.22), polarity (TPSA = 103.43 Å 2 ), and lipophilicity (XLOGP3 = 2.76). Table 3. Pharmacokinetics, drug-likeness, and medicinal chemistry of RV59.

Pharmacokinetics GI Absorption (High) BBB (Low)
Drug-likeness (Yes to all) Lipinski, Ghose, Veber, Egan, Muegge Bioavailability score 55% Medical Chemistry Synthetic accessibility: 2.95 Pharmacokinetics displayed high gastrointestinal absorptivity (GIA) and low blood-brain barrier (BBB) permeability; drug-likeness, the compound passed all the criteria for Ghose, Veber (GSK), Egan (Pharmacia), and muegge, and the Lipinski (Pfizer) rule-of-five of drug likeness and drug discovery. The bioavailability of the compound based on GIA indicated a score of 0.55 (55%) and synthetic accessibility of 2.95 evaluated according to the range from 1 (easy to make) to 10 (difficult to make).

Identification of Differentially Expressed Genes (DEGs) in CRC
DEGs from the microarray dataset were shown by identifying expressed genes between colon cancer samples and normal samples tallied from separate studies. Results showed that 536, 144, and 607 DEGs were respectively obtained from the GSE41328, GSE44861, and GSE74602 datasets. Eighty-five genes overlapped in these datasets as demonstrated in a Venn diagram ( Figure 2A). Red and blue dots in the heatmap and volcano plots respectively represent upregulated and downregulated genes ( Figure 2B-D). The expression density curve demonstrated that the expression scales of the three databases were between 0 and 16. Expression values of gene data in GSE41328 and GSE74602 were concentrated around 5 following standardization.

Validation of MYC/CXCL8/TIMP1 Expression Levels in CRC
Expression levels of the MYC/CXCL8/TIMP1 oncogenes were analyzed using UAL-CAN (http://ualcan.path.uab.edu/), a comprehensive open-access public online tool for analyzing cancer data [49] Expression levels of the MYC/CXCL8/TIMP1 oncogenes in CRC primary tumor samples (red) were compared to adjacent normal samples (blue), with p < 0.05 considered statistically significant ( Figure 3A-C). Interestingly, the analytical results showed high expressions of the MYC/CXCL8/TIMP1 oncogenes in tumor samples com-

Validation of MYC/CXCL8/TIMP1 Expression Levels in CRC
Expression levels of the MYC/CXCL8/TIMP1 oncogenes were analyzed using UAL-CAN (http://ualcan.path.uab.edu/), a comprehensive open-access public online tool for analyzing cancer data [49] Expression levels of the MYC/CXCL8/TIMP1 oncogenes in CRC primary tumor samples (red) were compared to adjacent normal samples (blue), with p < 0.05 considered statistically significant ( Figure 3A-C). Interestingly, the analytical results showed high expressions of the MYC/CXCL8/TIMP1 oncogenes in tumor samples compare d to normal samples. Furthermore, we used GECO, a gene expression correlation analytical tool, which displayed positive correlations of MYC with CXCL8, MYC with TIMP1, and TIMP1 with CXCL8. Positive Pearson correlation coefficients and p < 0.05 indicated statistical significance ( Figure 3D-F).

PPI Network and GO and KEGG Pathway Analysis
A functional interaction analysis was performed with STRING (https://stringdb.org/) to construct a PPI clustering network [51]. A confidence score of >0.7 was considered most significant. The GeneMenia tool (https://genemania.org/) and cytoscape software were used to build the PPI network. From analysis of the results, interactive networks were based on gene co-expression, co-localization, genetic interactions, and various pathways involved within the network. In short, we identified interactions of CXCL8 with MYC, TIMP1, MMP9 and ACKR; of MYC with CXCL8, TIMP1, MMP9, CDK4, CCND1, RIOX2, BCAT and DDX18; and of TIMP1 with MYC, CXCL8, MMP9, CCNA2, JUND, and STAT3 ( Figure 4A,B). The DAVID database (https://david.ncifcrf.gov/.jsp) was used to analyze the enriched GO including involved biological processes and biological pathways, with the criterion set to p < 0.05 ( Figure 4C,D).

PPI Network and GO and KEGG Pathway Analysis
A functional interaction analysis was performed with STRING (https://string-db. org/) to construct a PPI clustering network [51]. A confidence score of >0.7 was considered most significant. The GeneMenia tool (https://genemania.org/) and cytoscape software were used to build the PPI network. From analysis of the results, interactive networks were based on gene co-expression, co-localization, genetic interactions, and various pathways involved within the network. In short, we identified interactions of CXCL8 with MYC, TIMP1, MMP9 and ACKR; of MYC with CXCL8, TIMP1, MMP9, CDK4, CCND1, RIOX2, BCAT and DDX18; and of TIMP1 with MYC, CXCL8, MMP9, CCNA2, JUND, and STAT3 ( Figure 4A,B). The DAVID database (https://david.ncifcrf.gov/.jsp) was used to analyze the enriched GO including involved biological processes and biological pathways, with the criterion set to p < 0.05 ( Figure 4C,D).

MYC/CXCL8/TIMP1 Gene Co-Expression and Functional Enrichment Analysis
Gene co-expressions displayed the most factors which contributed to functional interactions. Herein, we determined the signaling network and KEGG pathway functional enrichment analysis, and found that the MYC, IL10, TP53, CXCL8, TIMP1, and CDK2 genes were co-expressed and were the most enriched ( Figure 5A). Moreover, the KEGG pathway enrichment analysis showed that co-expressed genes exhibited enrichment in the thyroid hormone signaling pathway, bladder cancer, cellular senescence, cell cycle, and HTLC-1 infections ( Figure 5B).
Gene co-expressions displayed the most factors which contributed to functional interactions. Herein, we determined the signaling network and KEGG pathway functional enrichment analysis, and found that the MYC, IL10, TP53, CXCL8, TIMP1, and CDK2 genes were co-expressed and were the most enriched ( Figure 5A). Moreover, the KEGG pathway enrichment analysis showed that co-expressed genes exhibited enrichment in the thyroid hormone signaling pathway, bladder cancer, cellular senescence, cell cycle, and HTLC-1 infections ( Figure 5B).

MYC/CXCL8/TIMP1 Expressions Were Correlated with Immune Cell Infiltration in Both Cancer and Normal Tissues
To identify associations of MYC/CXCL8/TIMP1 expressions with selected immune cells, we applied a correlation analysis between the above-mentioned oncogenes and immune infiltration cells (CD8+ and macrophages), where markers were adjusted by purity. As expected, results showed correlations of immune cell markers in colorectal adenocarcinoma (COAD), specifically CD8+ T cells and M2 macrophages ( Figure 6A-C), with p < 0.05 considered significant. Expressions of MYC/CXCL8/TIMP1 were also found to be correlated with infiltrating levels of CD8+ T cells, macrophages, and DCs, and red represents the most significant positive correlations with high amplification, while blue represents negative correlations ( Figure 6E-G). The infiltration level was compared to the normal level using a two-sided Wilcoxon rank-sum test.

MYC/CXCL8/TIMP1 Expressions Were Correlated with Immune Cell Infiltration in Both Cancer and Normal Tissues
To identify associations of MYC/CXCL8/TIMP1 expressions with selected immune cells, we applied a correlation analysis between the above-mentioned oncogenes and immune infiltration cells (CD8+ and macrophages), where markers were adjusted by purity. As expected, results showed correlations of immune cell markers in colorectal adenocarcinoma (COAD), specifically CD8+ T cells and M2 macrophages ( Figure 6A-C), with p < 0.05 considered significant. Expressions of MYC/CXCL8/TIMP1 were also found to be correlated with infiltrating levels of CD8+ T cells, macrophages, and DCs, and red represents the most significant positive correlations with high amplification, while blue represents negative correlations ( Figure 6E-G). The infiltration level was compared to the normal level using a two-sided Wilcoxon rank-sum test.

Drug Sensitivity Analysis of MYC/CXCL8/TIMP1 Oncogenes
To determine the drug sensitivity of MYC, CXCL8 and TIMP1, we used the GSCA tool to analyze the drug response (Figure 7). The correlation coefficients analysis, shows that upregulated gene expression is associated with drug resistance. From our analysis of results, we identified increased mRNA expression levels of MYC, CXCL8 and TIMP1 (indicated in orange bubbles), to be less sensitive to the drugs. Interestingly, high expression levels of MYC, TIMP1 and CXCL8 gene signatures, were shown to be resistance to Bx-912 (PDK-1 inhibitor) [59], navitoclax (Bcl-2 inhibitor) [60], vorinostat (HDAC and MYC inhibitor) [61,62] and tubastatin A (HDAC inhibitor) [63] among other FDA approved drugs. expressions were correlated with abundances of tumor infiltrates including CD8+ T cells, macrophages, and dendritic cells in colon cancer. Red represents significant positive correlations with high amplification, while blue represents negative correlations. The infiltration level was compared to the normal level using a two-sided Wilcoxon rank-sum test. p-value Significant Codes: 0 ≤ *** < 0.001 ≤ ** < 0.01 ≤ * < 0.05 ≤ . < 0.1.

Drug Sensitivity Analysis of MYC/CXCL8/TIMP1 Oncogenes
To determine the drug sensitivity of MYC, CXCL8 and TIMP1, we used the GSCA tool to analyze the drug response (Figure 7). The correlation coefficients analysis, shows that upregulated gene expression is associated with drug resistance. From our analysis of results, we identified increased mRNA expression levels of MYC, CXCL8 and TIMP1 (indicated in orange bubbles), to be less sensitive to the drugs. Interestingly, high expression levels of MYC, TIMP1 and CXCL8 gene signatures, were shown to be resistance to Bx-912 (PDK-1 inhibitor) [59], navitoclax (Bcl-2 inhibitor) [

Docking Results Displayed Strong Binding Energies between RV59 and the MYC/CXCL8/TIMP1 Oncogenes
Results from the in silico molecular docking analysis revealed unique binding affinities of the RV59 compound with selected target genes obtained from the PASS online prediction tool and DPT−COMPARE algorithm. Gibbs free binding energy results were obtained as follows: −7.6 kcal/mol for MYC, −7.7 kcal/mol for CXCL8, and −6.9 kcal/mol for TIMP1. In addition, a visualization analysis showed that docking of the small molecule (ligand) displayed shorter binding distances with the proteins (receptors) of 2.49 Å for MYC, of 2.03 Å for CXCL8, and of 2.51, 3.14, 3.4, and 3.29 Å for TIMP1 ( Figure 8A,B). Moreover, we compared results obtained from RV59 docking with standard inhibitors of MYC (vorinostat), CXCL8 (reparixin), and TIMP1 (ilomastat). Interestingly, the standard inhibitors exhibited lower binding energies for MYC (−6.3 kcal/mol) and CXCL8 (−6.3 kcal/mol), with the exception of TIM1 (7.4 kcal/mol), which showed a slightly higher binding energy compared to the RV59 compound ( Figure 9A-C). In addition, several interactions were identified for the ligand and protein complex, and these interactions included amino residues, a high number of conventional hydrogen bonds, van der Waals forces, carbon hydrogen bonds, Pi anions, Pi-sigma, Pi-Pi stacked, and amide Pi-stacked as shown in Table 4.

Docking Results Displayed Strong Binding Energies between RV59 and the MYC/CXCL8/TIMP1 Oncogenes
Results from the in silico molecular docking analysis revealed unique binding affinities of the RV59 compound with selected target genes obtained from the PASS online prediction tool and DPT−COMPARE algorithm. Gibbs free binding energy results were obtained as follows: −7.6 kcal/mol for MYC, −7.7 kcal/mol for CXCL8, and −6.9 kcal/mol for TIMP1. In addition, a visualization analysis showed that docking of the small molecule (ligand) displayed shorter binding distances with the proteins (receptors) of 2.49 Å for MYC, of 2.03 Å for CXCL8, and of 2.51, 3.14, 3.4, and 3.29 Å for TIMP1 ( Figure 8A,B). Moreover, we compared results obtained from RV59 docking with standard inhibitors of MYC (vorinostat), CXCL8 (reparixin), and TIMP1 (ilomastat). Interestingly, the standard inhibitors exhibited lower binding energies for MYC (−6.3 kcal/mol) and CXCL8 (−6.3 kcal/mol), with the exception of TIM1 (7.4 kcal/mol), which showed a slightly higher binding energy compared to the RV59 compound ( Figure 9A-C). In addition, several interactions were identified for the ligand and protein complex, and these interactions included amino residues, a high number of conventional hydrogen bonds, van der Waals forces, carbon hydrogen bonds, Pi anions, Pi-sigma, Pi-Pi stacked, and amide Pi-stacked as shown in Table 4.   Table 4.   Table 4.

Expressions of MYC/CXCL8/TIMP1 Oncogenes across Colon Cancer Cell Lines
To identify expression levels of the MYC/CXCL8/TIMP1 oncogenes in different CRC cell lines, we explored the expression database web tool, https://www.ebi.ac.uk/gxa/ home [64]. We used the RNA-sequence data of 675 commonly used human cancer cell lines for each gene. Results showed increased expression levels of the MYC/CXCL8/TIMP1 oncogenes across colon cancer cell lines ( Figure 10).

Expressions of MYC/CXCL8/TIMP1 Oncogenes across Colon Cancer Cell Lines
To identify expression levels of the MYC/CXCL8/TIMP1 oncogenes in different CRC cell lines, we explored the expression database web tool, https://www.ebi.ac.uk/gxa/home [64]. We used the RNA-sequence data of 675 commonly used human cancer cell lines for each gene. Results showed increased expression levels of the MYC/CXCL8/TIMP1 oncogenes across colon cancer cell lines ( Figure 10).

RV59 Displayed Anti-Proliferative and Cytotoxic Effects in NCI60 Human Colon Cancer Cell Lines
RV59 showed anticancer activities against NCI human colon cancer cell lines ( Figure  10). A single-dose treatment was administered at an initial dose of 10 μM, growth percentages revealed RV59 cytotoxic effects on the COLO 205, HCC-2998, and HCT-15 cell lines, and anti-proliferative activities on the HCT-116, HT29, KM12, and SW-620 cell lines ( Figure 11A). Since the compound showed potential anticancer effects at an initial dose of 10 μM, further dose-dependent investigations of the compound were performed to measure GI50, tumor growth inhibition (TGI), and 50% lethal concentration (LC50) values. Results showed potential anti-proliferative effects in a dose-dependent manner ( Figure  11B,C). Moreover, the in vitro IC50 results ranged 0.18-1.85 μM in the colon cancer cell lines, with HCT-15 being the more responsive at 0.

RV59 Displayed Anti-Proliferative and Cytotoxic Effects in NCI60 Human Colon Cancer Cell Lines
RV59 showed anticancer activities against NCI human colon cancer cell lines ( Figure 10). A single-dose treatment was administered at an initial dose of 10 µM, growth percentages revealed RV59 cytotoxic effects on the COLO 205, HCC-2998, and HCT-15 cell lines, and anti-proliferative activities on the HCT-116, HT29, KM12, and SW-620 cell lines ( Figure 11A). Since the compound showed potential anticancer effects at an initial dose of 10 µM, further dose-dependent investigations of the compound were performed to measure GI 50 , tumor growth inhibition (TGI), and 50% lethal concentration (LC 50 ) values. Results showed potential anti-proliferative effects in a dose-dependent manner ( Figure 11B,C). Moreover, the in vitro IC 50

Discussion
Advanced chemotherapy and targeted therapies still offer limited prolonged overall survival in CRC patients. One of the main causes of poor prognoses in patients is resistance to these therapeutic interventions [65]. Moreover, the molecular mechanisms through which cancer escapes chemotherapy and targeted therapy still remain elusive, mainly due to colon cancer's heterogenic properties. CRC is often diagnosed at a later, advanced stage, with distant metastasis present in most cases [66]. Therefore, there is a need for novel and effective targeted therapies, to improve patient's clinical outcomes and resistance in CRC patients. Protein kinases dysregulations have been reported to drive cancer, due to its association with genetic alterations, such as overexpression and mutations. Therefore, they have become pharmaceutical targets over the years [67]. RV59 is a novel small molecule derived in our lab from EGFR inhibitor osimertinib [68]. Our previous studies evaluated the anticancer activities of 20 nitrogen-substituted anthra [1,2c] [1,2,5]thiadiazole-6,11-dione derivatives, on cytoplasmic nuclear location sequence (NLS)-mutated Nrf2-transfecte, which promotes CRC tumor invasion and resistance to 5flurouracil (5-FU) chemotherapy. Among those 20 compounds, RV59 was more effective, overcame cNrf2-mediated chemoresistance, and suppressed tumor growth in colon cancer cells [69]. However, PK and toxicity analyses were not performed in that study. Further, in another study from our lab, we showed that RV59, exhibited a broad-spectrum of

Discussion
Advanced chemotherapy and targeted therapies still offer limited prolonged overall survival in CRC patients. One of the main causes of poor prognoses in patients is resistance to these therapeutic interventions [65]. Moreover, the molecular mechanisms through which cancer escapes chemotherapy and targeted therapy still remain elusive, mainly due to colon cancer's heterogenic properties. CRC is often diagnosed at a later, advanced stage, with distant metastasis present in most cases [66]. Therefore, there is a need for novel and effective targeted therapies, to improve patient's clinical outcomes and resistance in CRC patients. Protein kinases dysregulations have been reported to drive cancer, due to its association with genetic alterations, such as overexpression and mutations. Therefore, they have become pharmaceutical targets over the years [67]. RV59 is a novel small molecule derived in our lab from EGFR inhibitor osimertinib [68]. Our previous studies evaluated the anticancer activities of 20 nitrogen-substituted anthra[1,2-c][1,2,5]thiadiazole-6,11-dione derivatives, on cytoplasmic nuclear location sequence (NLS)-mutated Nrf2-transfecte, which promotes CRC tumor invasion and resistance to 5-flurouracil (5-FU) chemotherapy. Among those 20 compounds, RV59 was more effective, overcame cNrf2-mediated chemoresistance, and suppressed tumor growth in colon cancer cells [69]. However, PK and toxicity analyses were not performed in that study. Further, in another study from our lab, we showed that RV59, exhibited a broad-spectrum of cytotoxicity against various cancer cells, and interestingly, the compound displayed less cytotoxic effects as compared to chemotherapeutic agent, doxorubicin in normal tissues [70].
Herein, we used a bioinformatics simulation study to further expand on that previous study; however, the present study focused more on identifying and validating oncogenes associated with chemoresistance and alteration of oxaliplatin or doxorubicin treatment responses in advanced CRC. Studies showed that accumulation of tumor recurrence from cancer stem cells (CSCs), and metastasis often occur post-treatment in CRC, which leads to therapeutic resistance [65,71]. Using computer-based PASS and DTP-COMPARE drug target predictive tools, we identified MYC/CXCL8/TIMP1 as target genes for RV59; in addition, an analysis of the results also predicted that all probability active (Pa) values were greater that all predicted probability inactive (Pi) values, indicating the inhibitory or antagonistic activities of RV59 on MYC and CXCL8, among other identified activities. In the early stages of drug discovery and development, the disposition of a compound is assessed in terms of its absorption, distribution, metabolism, and excretion (ADME), with the final goal of identifying potential effectiveness of a medicine for patients [72]. After assessing the ADME of RV59, bioavailability radar displaying all six physical properties of the compound, showed that the compound successfully met the minimum requirements of drug-likeness, with an Mw of 352.41 g/mol, flexibility (rotations = 4), solubility (log S (ESOL) = −3.74), saturation (fraction Csp3 = 0.22), polarity (TPSA = 103.43 Å 2 ), and lipophilicity (XLOGP3 = 2.76), and all the values were within recommended values as described in Figure 1. The compound also passed the criteria for Ghose, Veber (GSK), Egan (Pharmacia), and muegge and the Lipinski (Pfizer) rule-of-five of drug likeness and drug discovery. The bioavailability of the compound based on GIA indicated a score of 0.55 (55%), which indicates acceptable PK properties. Moreover, we explored the UALCAN bioinformatics tool, and validated increased expressions of MYC/CXCL8/TIMP1 oncogenic signatures in CRC primary tumors and compared them to adjacent normal samples. The genes also exhibited positive correlations among each other with p < 0.05 indicating statistical significance.
In a further analysis, we determined interactions of the MYC/CXCL8/TIMP1 oncogenes at the protein level. Interestingly, the PPI network analysis from two independent databases showed interactions of all of the proteins with each other, as well as enrichment in GO terms, including biological processes and biological pathways involved in CRC, with p < 0.05 indicating statistical significance ( Figure 4). MYC plays significant roles in tumorigenesis and therapeutic resistance [15]. In 2019, Han et al. demonstrated that MYC induces immunogenic cell death of tumor cells, which leads to increased T cell infiltration and upregulation of the PD-L1 immune checkpoint protein in the TME [73].
Additional studies also showed that inhibition of programmed cell death protein 1 (PD1) by pembrolizumab enhanced regulation of cytotoxic T-cell tumoricidal activities, which ultimately leads to increased expression levels of tumor-enhanced CXCL8, which subsequently induces infiltration of tumor-associated macrophages (M2) in the immune system [74,75]. Others showed that the TIMP1 gene derived from tumor cells creates a metastatic niche, to which circulating tumor cells cling and promote CRC metastasis [76]. Hence, these finding validated the potential contribution of the MYC/CXCL8/TIMP1 oncogenic signatures to chemoresistance and resistance to targeted treatment. To identify associations of MYC/CXCL8/TIMP1 gene expressions with selected immune cells in CRC, we applied a correlation analysis between the above-mentioned oncogenes with immune infiltration cells (CD8+ and macrophages), where markers were adjusted by purity. As expected, MYC/CXCL8/TIMP1 gene expression levels displayed associations with tumor purity and were positively correlated with infiltrating levels of CD8+ T cells and macrophages (p < 0.05), and were correlated with the abundances of tumor infiltrates including CD8+ T cells, macrophages and DCs in CRC ( Figure 6).
We further predicted protein-ligand interactions using molecular docking between RV59 and the MYC/CXCL8/TIMP1 gene complex, and discovered that RV59 displayed strong binding energies to MYC (−7.6 kcal/mol), CXCL8 (−7.7 kcal/mol), and TIMP1 (−6.9 kcal/mol), and short binding distances with these protein (receptors): MYC (2.49 Å), CXCL8 (2.03 Å), and TIMP1 (2.51, 3.14, 3.4, and 3.29 Å). These results were compared to standard inhibitors of MYC (vorinostat), CXCL8 (reparixin), and TIMP1 (ilomastat). Interestingly, the inhibitors exhibited lower binding energies for MYC (−6.3 kcal/mol) and CXCL8 (−6.3 kcal/mol), with the exception of TIMP1 (7.4 kcal/mol) compared to RV59. Moreover, RV59 showed anticancer activities against NCI human colon cancer cell lines with single-dose treatment of 10 µM, toward the COLO 205, HCC-2998, and HCT-15 cell lines, and anti-proliferative activities toward the HCT-116, HT29, KM12, and SW-620 cell lines, in dose-dependent manners. measured Results of measuring the GI 50 , TGI, and LC 50 levels showed potential anti-proliferative effects in dose-dependent manners ( Figure 11B,C). Moreover, the in vitro IC 50 results ranged 0.18-1.85 µM on colon cancer cell lines, with HCT-15 cells being the most responsive at 0.18 µM, followed by HCT116 at 0.28 µM, SW-620 at 0.29 µM, HT-29 at 0.47 µM, KW12 at 0.84 µM, and COLO 205 at 1.15 µM, with HCC-2998 cells at 1.85 µM showing the least responsiveness compared to the aforementioned cell lines. This suggests that RV59 exhibits drug-like characteristics, and is a potential oral drug candidate. In summary, the findings from this study revealed the anticancer activities of RV59 in CRC, and highlight additional avenues for RV59 as a potential treatment for CRC chemotherapy and resistance to targeted therapy, particularly by targeting the MYC/CXCL8/TIMP1 signaling pathway. The compound is currently being assessed for its therapeutic potential effectiveness in CRC and breast cancer both in vitro and in vivo in our laboratory.

Conclusions
In conclusion, we revealed the anticancer activities of RV59 against NCI human colon cancer cell lines both as a single dose and dose-dependent treatment, and also demonstrated the MYC/CXCL8/TIMP1 signaling pathway, which is responsible for resistance to both chemotherapies and targeted therapies, to be a potential RV59 drug target. Moreover, the in silico molecular docking study exhibited putative binding affinities of RV59 with the above-mentioned oncogenes, which were even higher than the standard inhibitors of these genes. Currently in our lab, there are ongoing cell and animal experiments to evaluate the therapeutic effectiveness of RV59 in CRC and breast cancer.

Data Availability Statement:
The datasets generated and/or analyzed in this study are available on reasonable request.