Differential Interactome Based Drug Repositioning Unraveled Abacavir, Exemestane, Nortriptyline Hydrochloride, and Tolcapone as Potential Therapeutics for Colorectal Cancers

There is a critical requirement for alternative strategies to provide the better treatment in colorectal cancer (CRC). Hence, our goal was to propose novel biomarkers as well as drug candidates for its treatment through differential interactome based drug repositioning. Differentially interacting proteins and their modules were identified, and their prognostic power were estimated through survival analyses. Drug repositioning was carried out for significant target proteins, and candidate drugs were analyzed via in silico molecular docking prior to in vitro cell viability assays in CRC cell lines. Six modules (mAPEX1, mCCT7, mHSD17B10, mMYC, mPSMB5, mRAN) were highlighted considering their prognostic performance. Drug repositioning resulted in eight drugs (abacavir, ribociclib, exemestane, voriconazole, nortriptyline hydrochloride, theophylline, bromocriptine mesylate, and tolcapone). Moreover, significant in vitro inhibition profiles were obtained in abacavir, nortriptyline hydrochloride, exemestane, tolcapone, and theophylline (positive control). Our findings may provide new and complementary strategies for the treatment of CRC.

Introduction

Colorectal cancer (CRC) is the most commonly diagnosed cancer in women and men worldwide. It occurs in the colon or rectum and affects the large intestine or large bowel. It is also known as bowel cancer, rectal cancer and colon cancer. The predicted new cases are over 1, 9 million and the number of deaths is 9,35,000 in 2020 (Sung et al., 2021). It is assumed to increase to 2.2 million new cases and 1.1 million deaths by 2030 (Douaiher et al., 2017). According to American Cancer Society, only 4 out of 10 CRC patients are detected at an early stage. If detected at an early stage, the 5-years survival rate can be as high as 90%. Otherwise, if metastases occur in all parts of the body, the 5-years survival rate drops to 14% (Rahib et al., 2014).

For many years, various drugs have been improved for the treatment of CRC, 5-fluorouracil, bevacizumab, cetuximab (Mauricio, 2019) also in combination with 5-fluorouracil/leucovorin with either oxaliplatin (FOLFOX) or irinotecan (FOLFIRI) (Neugut et al., 2019). However, patients with advanced CRC are resistant to 5-fluorouracil (Douillard et al., 2000). Therefore, new therapeutic agents may be needed for treatment.

In the last decades, many systems biology approaches have been used to overcome problems in the diagnosis, prognosis, or therapy of many cancers including CRC (Turanli et al., 2017a). For instance, five potential target genes (DUSP4, ETV5, GNB5, NT5E, and PHLDA1 to overcome cetuximab resistance in CRC utilizing gene expression profiling and mathematical modelling, and the results were subsequently validated through wet-lab studies (Park et al., 2019). Furthermore, several source genes and target proteins were identified through the employment of differential gene expression analysis, protein-protein interactions (PPIs), and genome-scale metabolic modelling (Zhang et al., 2019; Ilyas et al., 2020).

Physical interactions between proteins in all living organisms are the basis for cellular signaling pathways that mediate essential biological processes such as gene expression and metabolism (Sevimoglu and Arga, 2014; Calimlioglu et al., 2015). By integrating high-throughput genomic data with PPI networks under the concept of network science, disease mechanisms and genes associated with diseases or proteins can thus be discovered. Clarification of disease mechanisms associated with biological processes and identification of disease candidate genes or proteins through the use of PPI networks may suggest potential biomarkers and drug targets for these diseases (Safari-Alighiarloo et al., 2014).

Biomarkers play a crucial role in the concept of personalized medicine in identifying subtype phenotypes, identifying the convenient therapeutic approach and estimating clinical course and prognosis. Due to the limited efficiency of molecular biomarkers in the diagnosis and treatment of complex diseases such as cancer, researchers have recently focused on the detection of biological clusters of molecules (mostly gene and/or protein clusters) that have functional relationships with each other and are referred to as systems biomarkers (Turanli et al., 2017a). Accordingly, studies have shown that the identified systems biomarkers play an active role in the diagnosis and prognosis of diseases, as well as in the development of treatment strategies (Gov and Arga, 2017; Kori and Arga, 2018; Sevimoglu et al., 2018).

The differential interactome approach relies on important alterations that appeared in protein-protein interactions (PPIs) between phenotypes. The differential interactome algorithm allows to predict the probability distributions for each possible co-expression profile of gene pairs (encoding proteins that interact with each other) across phenotypes and to determine the uncertainty of whether a PPI matches the corresponding phenotype (Ayyildiz et al., 2017). This approach has shown success by being utilized effectively in various cancers and subtypes. (Ayyildiz et al., 2017; Turanli et al., 2019a; Turanli et al., 2019b; Gulfidan et al., 2020; Caliskan et al., 2021).

Drug repositioning (DR), in which existing drugs are repurposed for a new therapeutic indication is a promising approach because it reduces time and cost in drug development and circumvents problems due to safety and efficacy issues (Shim and Liu, 2014; Würth et al., 2016; Xue et al., 2018). Recently, a review paper pointed out repositioning in drug-resistant CRC (Nowak-Sliwinska et al., 2019). Citalopram, amantadine, and captopril are repurposed drugs for the prevention or treatment of disease (Koh et al., 2014; Van Noort et al., 2014; Diaz-Carballo et al., 2015). Also, there are computerized efforts to reuse drugs such as GW-8510, etacrynic acid, ginkgolide A, and 6-azathymin with the use of Functional Module Connectivity Map for CRC (Chung et al., 2014). Currently, the identification of candidate biological targets and new potential drugs could be done by using in silico methods for DR by collecting clinical data at different omics levels and analyzing them in systematic and integrative pipelines (Turanli et al., 2019a). As an instance, the well-established drug metformin is used to treat Type 2 diabetes and was found to be a preventative agent in CRC via in silico methods (Higurashi and Nakajima, 2018). In addition, signature-based DR is another strategy in the DR approach to identify existing drugs for potential treatment or to fulfill new indications. This method considers gene expression signatures and compares drug-gene and disease-gene expression profiles. (Yella et al., 2018). One of the earliest examples of signature-based approaches is “Connectivity Map” and a public resource used to find small molecules and mechanisms of their action, chemicals or physiological processes, diseases, and drugs (Lamb et al., 2006). Another category of DR is network-based DR, which utilize to identify molecular mechanisms and diagnostic/prognostic biomarkers in many diseases including cancer. Besides, network modelling is an important approach for computational drug repositioning, forming a triangle of disease, genes, and drug (Yella et al., 2018).

Considering the urgent need for the development of new therapeutic strategies in CRC, in this study, we applied a network-based DR approach to propose novel drug candidates for CRC treatment. For this purpose, taking into account the fact that colorectal carcinomas arise predominantly from adenomas, we evaluated comprehensive datasets for human colorectal adenomas and carcinomas together, employed the differential interactome algorithm and evaluated the potential of differentially interacting proteins as drug targets. Then, we repurposed candidate drugs which are later analyzed in silico by molecular docking simulations and in vitro by using viability assays to determine their potential in CRC treatment.

Methods

Gene Expression Datasets

Two comprehensive transcriptome datasets were employed: 1) the microarray dataset (GSE8671) (Sabates-Bellver et al., 2007) which compiled on Affymetrix Human Genome U133 Plus 2.0 Array platform (Affymetrix Inc.,Santa Clara, CA, United States) with 32 colorectal adenoma samples and 32 matched normal tissue samples obtained from NCBI-Gene Expression Omnibus (GEO) (Barrett et al., 2013), and 2) the RNA-seq COAD and READ datasets (normalized as FPKM) from The Cancer Genome Atlas (TCGA) consisting of 644 tumor tissue samples, and 51 normal tissue samples (Tomczak et al., 2015). TCGA-COAD consisted of 478 primary tissue samples and 41 normal colon tissue samples. TCGA-READ consisted of 166 primary tissue samples and 10 normal rectal tissue samples.

Protein-Protein Interactions Data

BioGRID database (v.3.5.167) (Chatr-aryamontri et al., 2017) was used, which contains 35,688 physically and experimentally detected PPIs among 8,570 human proteins for the human protein interactome. Filtering the interactome dataset with proteins encoded by the genes represented in the transcriptome datasets resulted in a network of 34,603 PPIs among 8,322 proteins considering the TCGA dataset and 32,259 PPIs among 7,951 proteins considering the GSE8671 dataset.

Identification of Differential Interactome and Differentially Interacting Proteins

The differential interactome algorithm was recruited to identify the differential PPIs (dPPIs) between the tumor (or adenoma) phenotype and the normal phenotype, considering the relative observation frequencies (q-value) of each PPI as previously described (Ayyildiz et al., 2017; Gulfidan et al., 2020). Briefly, genes were grouped into three different levels as −1, 0, 1 according to their expression levels within each sample. When the expression level of a gene was lower than the average expression value, it was labeled as “−1’; when the expression level of a gene was higher than the average, it was labeled as “1”; otherwise, it was labeled as “0”. The false discovery rate was 0.05. According to the three-level expression categorization, there were nine possible gene expression states (i.e., [0 0], [0 1], [0 −1], [1 0], [1 1], [1 −1], [−1 0], [−1 1], [−1 −1]) for each interacting protein pair. The number of times the conditions occurred in the normal group (N₀) and the number of times they occurred in the tumor group (N₁) were calculated. Taking into account the imbalance between the sample sizes of the groups, the count parameters were normalized considering the total sizes of the normal (N_N) and tumor (N_T) groups (the maximum possible numbers of N₀ and N₁, respectively), and the normalized observation frequencies in each group were obtained. The q value (the estimate of the probability of that state occurring in the tumor state) was calculated as follows:

$q= \frac{\frac{N_1}{N_T}}{\frac{N_0}{N_N}+\frac{N_1}{N_T}}$

PPIs with a normalized observation frequency in either the normal or tumor phenotype greater than 20% and q-values of less than 0.10 (significantly repressed in the tumor phenotype) or greater than 0.90 (significantly activated in the tumor phenotype) were considered dPPIs. Differential interacting proteins (DIPs) were defined as proteins that exhibit significant alterations in their interaction patterns during the transition from normal phenotype to tumor state (Gulfidan et al., 2020), and categorized into two groups depending on their interaction patterns: 1) DIPs with repressed interactions under tumor conditions, and 2) DIPs with activated interactions under tumor conditions. Modules were constructed around DIPs using their dPPIs and visualized using Cytoscape 3.7.2 (Shannon et al., 2003).

Gene Set Over-Representation Analysis

Gene set over-representation analyses were carried out using the ConsensusPathDB (ver.34) (Kamburov et al., 2013). Preferred data sources for metabolic and signalling pathways were KEGG (ver.88.0) (Kanehisa et al., 2017), Reactome (ver.64) (Fabregat et al., 2017), and Biocarta (ver.2009_05_12) (Nishimura, 2001), and Gene Ontology annotations (Ashburner et al., 2000) were employed to represent associated biological processes. Statistical significance was defined by a p-value cut-off < 0.05 for all functional enrichment analyses. Each p-value was subsequently converted to a z-score by using the inverse normal cumulative distribution.

Prognostic Power Analysis

The prognostic analyses were performed using two datasets, GSE17536 (Smith et al., 2010) and TCGA-COADREAD containing patient survival data. The dataset of GSE17536 consists of 177 samples (24 patients with stage I, 57 patients with stage II, 57 patients with stage III, and 39 patients with stage IV) obtained from a patient cohort having the average age of 65.5 ± 13.1, whereas the dataset of TCGA COADREAD consists of 466 samples (89 patients with stage I, 170 patients with stage II, 130 patients with stage III, 62 patients with stage IV, and 15 patients with no stage information) acquired from a patient cohort having average age of 66.9 ± 12.5.

Survival analyses were performed by stratifying patients into high- and low-risk groups based on their prognostic index (PI), which is the linear component of the Cox model (PI = β₁x₁ + β₂x₂ + … + β_px_p, where β_p is the coefficient obtained from the Cox fit, x_p is the expression value of each gene) to investigate the prognostic performance of each DIP module. Analyses were performed with the web server SurvExpress (Aguirre-Gamboa et al., 2013) using datasets with clinical data. Signatures of survival in each risk group were estimated by Kaplan–Meier curves and Hazard Ratios (HR). The statistical significance of each curve was assessed by the cut-off for the log-rank p-value < 0.05. The hazard ratio (HR = (O1/E1)/(O2/E2)) was calculated to discover the significance of the survival curves based on the ratio between the relative death rate in group 1 (O₁/E₁) and the relative death rate in group 2 (O₂/E₂), where O denotes the observed number of deaths, and E denotes the expected number of deaths.

Drug Repositioning

GeneXpharma (Turanli et al., 2017b), which is a publicly available platform presenting 50,304 gene-drug interactions among 4,344 genes and 11,939 drugs and employing statistical tests for the disease-gene-drug triad, was used for the network-based DR considering protein targets CDKN2A, GSK3B, HDAC2, and PML. A hypergeometric probability test was used to associate drugs with target proteins, and simulation results with p < 10⁻³ were accepted as statistically significant.

Molecular Docking

The 3-D structures of target proteins were obtained from Protein Data Bank (PDB) (Berman et al., 2002). The PDB identifiers were 1DC2 (CDKN2A) (Yuan et al., 2000), 6Y9R (GSK3B) (Buonfiglio et al., 2020), 6WBZ (HDAC2) (Yu et al., 2020), and 5YUF (PML) (Wang et al., 2018). Ligand binding sites of proteins were determined by PDBe (Velankar et al., 2016) and BIOVIA Discovery Studio (ver.20.1) (BIOVIA, 2016). The structures of candidate drugs were taken from PubChem (Kim et al., 2019). Autodock Vina software (ver.1.1.2) (Trott and Olson, 2019) was used for molecular docking analyses. All water molecules were deleted, and all polar hydrogens were added to the structure of a protein in the preparation of macromolecules. Molecular docking was performed five times for each ligand and each simulation yielded twenty poses. Exhaustiveness was set to eight for all docking calculations. Binding affinities were estimated to determine the importance of binding between protein targets and drug candidates. The top-scoring pose (with the lowest binding free energy) was selected for further analysis.

Drugs

Abacavir (Selleckchem, S5215), Ribociclib (Selleckchem, S7440), Exemestane (Selleckchem, S1196), Voriconazole (Selleckchem, S1442), Bromocriptine mesylate (R&D systems, 0427/50), Tolcapone (Selleckchem, S4021) were dissolved in dimethyl sulfoxide (DMSO) and then diluted in 1x PBS. For each treatment step, 100 mM stock solution was diluted to a final concentration. Total DMSO concentration did not exceed 0.1% and control groups (CTRL) were also treated with the same concentration of DMSO. Nortriptyline hydrochloride (Selleckchem, S3698) and Theophylline (Sigma, T1633) were dissolved in 1x PBS.

In Vitro Cell Viability Assay

CCD-841-Con cells (ATCC, CRL-1790, human healthy colon epithelial cell line) and HCT-116 cells (ATCC-CCL-247, human colorectal carcinoma cell line) were cultured in Dulbecco’s modified Eagle’s medium and Eagle’s Minimum Essential Medium, both supplemented with penicillin (100 units/ml), streptomycin (100 g/ml), and 10% fetal calf serum in a humidified atmosphere of 5% CO₂ and 95% air at 37°C. Later, cells were treated with varying abacavir, voriconazole, nortriptyline hydrochloride, tolcapone, and theophylline concentrations (10, 50, 100, 150, and 200 µM) in the first phase. Another wide exposure dose was used as 10, 50, 100, and 200 nM and 200–500 µM for bromocriptine mesylate, budesonide, exemestane, voriconazole, and theophylline. Also, nortriptyline hydrochloride concentration in the range of 10–50 µM was used to identify IC₅₀. Cell viability was determined using MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2H-tetrazolium bromide) reduction by viable cells following the exposure durations of 24, 48, and 72 h 2 µl of MTT solution was mixed with 100 μL medium and incubated for 3 h at 37°C and 5% CO₂. At the end of this incubation period, the medium with MTT was discarded and 100 μL DMSO was added. The formazan crystals were dissolved by shaking at 150 rpm for 10 min. The intensity (OD) of the color formed was measured with a spectrophotometer with a microplate reader at a wavelength of 590 nm (Reference filter; 660 nm). Calculations were performed related to the absorbance of control samples which was equaled to 100%.

Results

Differentially Interacting Proteins, Their Modules, and Biological Interpretation

Two comprehensive transcriptome datasets associated with colorectal adenomas and tumors were recruited to apply differential interactome methodology for predicting high probability PPIs in tumor states and identifying differential PPIs. Two datasets were analyzed independently, and common dPPI signatures were considered in further analyzes. In the current formulation, upper and lower bounds (0.90 and 0.10, respectively) were used for q-values which represent the probability estimates. PPIs that had a normalized frequency of observation in either the normal or tumor phenotype greater than 20% were considered as significant. As a result, a total of 2,214 differential PPIs (dPPIs) were identified in the GEO dataset and 1,625 dPPIs were identified in the TCGA dataset. dPPIs were further classified as “significantly repressed in tumor phenotype” (if q < 0.10) or “significantly activated in tumor phenotype” (if q > 0.90). To this end, among the dPPIs, 718 interactions were repressed and 1707 interactions were activated in the GEO dataset, while 81 interactions were repressed and 1,557 interactions were activated in the TCGA dataset (Figure 1A).

FIGURE 1

FIGURE 1. (A) Interaction networks of both datasets and common interactions. GEO-specific edges were shown in blue, TCGA-specific edges were shown in red and common interactions were shown in green. (B) Common modules in both datasets. Node size increases with the size of each module, edge thickness increases with the number of shared interactions between modules. The color scale of nodes changes according to the presence of each hub protein in other modules.

The scale-free topology of the differential interactome network highlights the presence of hubs termed differentially interacting proteins (DIPs), indicating major changes in their interaction patterns during tumorigenesis (i.e., transition from “normal” to “tumor” phenotypes) (Ayyildiz et al., 2017). In this study, we identified DIPs and constructed a module around each DIP with its interaction partners, and the corresponding module was named adding the letter “m” to the beginning of the central DIPs (e.g., mMYC). The number of modules with more than 5 PPIs was 86 in TCGA and 177 in the GEO dataset. In both datasets, modules with at least 5 PPIs were considered and filtered based on their common interactions only. The number of common modules in two datasets (GEO and TCGA) were 24 (Figure 1B, Supplementary Table S1).

Among these modules, mMYC was the largest module with 116 members. Besides, the MYC protein was the most common protein observed in other modules. Moreover, mPSMB5, mCCT7, mRAN, mAPEX1, mHSD17B10, and mHDGF were also modules with more than 10 interactions. mPSMB5 had 33 interactions, as well as 30 interactions for the mCCT7, and mRAN had 29 interactions. For others, the number of interactions was 15 and 11, respectively.

Considering the members of the modules, gene set over-representation analyses were performed to identify biological processes, molecular pathways, and cancer hallmarks associated with the DIPs (Figure 2, Supplementary Table S1). The members of the modules were significantly enriched with several biological processes such as deubiquitination, protein modifications, cell cycle, neddylation, and transcription; signaling pathways such as PI3K-Akt, NOTCH, and Wnt signaling; cancer pathways such as transcriptional misregulation of cancer and microRNAs in cancer; and cancer pathways associated with oncogenic viruses such as Hepatitis B, Human papillomavirus (HPV), Epstein-Barr virus, and Human T-cell leukemia virus 1 (Figure 2A, Supplementary Tables S2, S3). Furthermore, the members of the modules were significantly associated with all cancer hallmarks (Senga and Grose, 2021), i.e., sustaining proliferative signaling, inducing angiogenesis, resisting cell death, deregulating cellular energetics, evading growth suppressors, activating invasion and metastasis, enabling replicative immortality, avoiding immune destruction, and genome instability and mutation (Figure 2B).

FIGURE 2

FIGURE 2. (A) Heatmap indicating top-scored 30 pathways enriched with the modules. Cells were colored depending on statistical significance (i.e., z-score). (B) Circos plot linking cancer hallmarks to modules via the number of module members associated with the cancer hallmarks.

Evaluation of DIP Centered Modules as Potential Prognostic Systems Biomarkers

The evaluation of prognostic performances of DIP-centered modules was performed via Kaplan-Meier survival analyses through partitioning low- and high-risk groups regarding the expression levels of module genes. To cross-validate the results, the survival analyses were independently evaluated using two different datasets. As a result, prognostic performances of six modules (namely, mAPEX1, mCCT7, mHSD17B10, mMYC, mPSMB5, and mRAN) were found significant in both CRC datasets considering log-rank test p-values (p < 0.01) and hazard ratios (HR > 2) (Figure 3).

FIGURE 3

FIGURE 3. Kaplan–Meier plots estimating survival probabilities of patients considering GEO and TCGA datasets. Patients were stratified considering their prognostic index determined through gene expression profiles. The low-risk group was shown in green, and the high-risk group was shown in red.

Identification of Candidate Drugs for CRC Treatment Through Drug Repositioning

To identify potential drug targets, modules were further evaluated considering their sizes, members, prognostic performances, and enriched processes, pathways, and cancer hallmarks. Considering its relevance, we focused on the interactions of MYC in the module and filtered the proteins according to available 3-D structure information in PDB (Berman et al., 2002), and drug-gene interaction information in geneXpharma (Turanli et al., 2017b). The structures of 60 proteins out of 116 proteins were found in PDB, and drug interactions were identified for 39 proteins with available structural information. Durg repositioning simulation through geneXpharma resulted in 651 drugs associated with 39 target proteins. Then, these drugs were filtered according to FDA-approval, association with CRC in literature, and having the 3-D structures on PubChem (Kim et al., 2019). Finally, we identified eight drugs (abacavir, voriconazole, exemestane, nortriptyline hydrochloride, tolcapone, bromocriptine mesylate, ribociclib, and theophylline) targeting 4 proteins (PML, GSK3B, CDKN2A, HDAC2) (Figure 4).

FIGURE 4

FIGURE 4. The network between drug targets and drugs obtained from geneXpharma. Edge thickness varies according to the hypergeometric probability (p) describing the significance of the drug-target association.

In silico Investigation of Potential Drug-Target Interactions by Molecular Docking Simulations

In silico validation studies were carried out to evaluate the potential candidate drugs for the therapeutic strategies in CRC prior to in vitro viability assays. For this purpose, the 3-D structures of target proteins, i.e., 1DC2 (CDKN2A), 6Y9R (GSK3B), 6WBZ (HDAC2), and 5YUF (PML), and candidate drugs were retrieved from PDB (Berman et al., 2002) and PubChem (Kim et al., 2019), respectively, and the binding affinities of the drugs to their targets were estimated via protein-ligand molecular docking simulations using Autodock Vina (Trott and Olson, 2019) (Figure 5).

FIGURE 5

FIGURE 5. The structures of the bound protein-ligand complexes for each drug target and their binding affinities (mean ± standard deviation (S.D).

The binding affinity of the positive control case (theophylline) was predicted to be −5.1 kcal/mol. The most significant binding affinities were obtained for bromocriptine mesylate, nortriptyline hydrochloride, tolcapone, and exemestane, which were predicted as −9.7, −8.2, −7.3, and −7 kcal/mol, respectively. Similarly, the binding affinities of ribociclib, abacavir, and voriconazole were −6.9, −5.9, −5.7 kcal/mol, respectively. Besides, the binding affinity of nortriptyline hydrochloride which targets PML protein was -6.1 kcal/mol. According to the molecular docking simulations, all potential drugs (bromocriptine mesylate, tolcapone, nortriptyline hydrochloride, ribociclib, exemestane, voriconazole, and abacavir) display significantly higher binding affinities when docked to their target proteins, compared to the positive control case (Figure 5). The details of the top-scoring pose of five docking simulations for each ligand are displayed in Supplementary Figure S1.

In vitro Cell Viability Assay of Repurposed Drugs

To have an insight into the antitumor potential of drugs for CRC treatment, in vitro cell viability assay was performed in HCT-116 (carcinoma) cell line. Concentration ranges as 10–200 uM for abacavir and tolcapone, 5–100 uM for ribociclib, 200–500 uM for exemestane, voriconazole, and bromocriptine mesylate, 10–50 µM for nortriptyline hydrochloride, and 20–500 µM for theophylline were tested to find out the IC₅₀ values in HCT-116 cells.

After 24, 48, and 72 h drug treatments, IC₅₀ values were determined as 100 μM at 72 h for abacavir, 100 μM at 72 h for ribociclib, 500 μM at 48 h for exemestane, 40 μM at 48 h for nortriptyline hydrochloride, 500 μM at 72 h for theophylline and 10 μM at 72 h for tolcapone (Figure 6). Bromocriptine mesylate and voriconazole did not decrease viability at tested concentrations and time points.

FIGURE 6

FIGURE 6. Effects on viability in HCT116 cells following 24, 48 and 72 h of drug treatments. Data denote mean ± S.D. *p < 0.05 vs. CTRL (n = 3); **p < 0.01 vs. CTRL (n = 3); ***p < 0.001 vs. CTRL (n = 3). Viability following of Abacavir (A), Ribociclib (B), Exemestane (C), Nortriptyline (D), Theophylline (E), Tolcapone (F) treatment in HCT116 cells.

In addition, the toxicity of the drugs was tested in and CCD-841-Con healthy epithelial cells in the obtained IC₅₀ values. Ribociclib was the only drug that caused death in healthy cells when compared to control (Figure 7).

FIGURE 7

FIGURE 7. Effects on viability in CCD 841 CoN cells following drug treatments in IC50 conditions. Data denote mean ± S.D. ***p < 0.001 vs. CTRL (n = 3). Abacavir, 100 uM and 72 h; Ribociclib, 100 uM and 24 h; Exemestane, 500 uM and 48 h; Nortriptyline, 40 uM and 72 h; Theophylline 500 uM and 72 h; Tolcapone, 10 uM, and 72 h treatment in CCD 841 CON cells.

Discussion

In the present study, we aimed to identify potential systems biomarkers and drug candidates for the treatment of CRC. To this end, we considered the transcriptome profiles of colorectal adenomas and carcinomas together, as colorectal carcinomas predominantly arise from adenomas, and applied the differential interactome approach for the first time in CRC to identify prognostic system biomarkers and drug repositioning candidates. The potential of drug candidates was demonstrated by in silico molecular docking simulations and in vitro cytotoxicity assays on CRC and healthy cell lines. In summary, 24 modules were found in both GEO and TCGA datasets, six of which showed tumor prognosis. As a result of the detailed evaluation of modules within the relevant literature, MYC, a transcription factor that coordinates many biological processes, came to the fore in all aspects. The MYC-centered module (mMYC) showed high prognostic performance in both datasets; MYC was significantly overexpressed in human cancers (Dang et al., 2009), and it promotes cell proliferation, sensitizes to apoptosis, and induces cells to undergo apoptosis (Massó-Vallés and Soucek, 2020). Therefore, it was promising to target MYC and its interacting partners as therapeutic targets in cancer therapy. On the other hand, previous studies reported that targeting MYC was problematic (Wolf and Eilers, 2020) and direct targeting of MYC for cancer therapy has not been proposed as a rational strategy (Chen et al., 2018); therefore, we hypothesized that breaking the bond between MYC and its interacting partners would be a promising strategy for the treatment of CRC. For this reason, we focused on four MYC-interacting proteins, i.e., Promyelocytic leukemia protein (PML), Glycogen synthase kinase-3 beta (GSK3B), Cyclin-Dependent Kinase Inhibitor 2A (CDKN2A), and Histone deacetylase 2 (HDAC2), which were reported to be mis-regulated in cancer in previous studies (Guan and Kao, 2015; Domoto et al., 2016; Zhao et al., 2016; Shan et al., 2017).

According to the GeneCards database (Stelzer et al., 2016), the PML protein is a member of the triple motif family (TRIM). This phosphoprotein localizes in nuclear bodies where it functions as a transcription factor and tumor suppressor by preventing cells from growing and dividing in an uncontrolled manner. The tumor suppressor PML regulates cell cycle, apoptosis, senescence, migration, angiogenesis, and DNA repair pathways as well as the p53 response to oncogenic signals (Guan and Kao, 2015). GSK3B protein is a serine-threonine kinase that belongs to the glycogen synthase kinase subfamily. It is a regulator of glucose homeostasis and plays tumor promoting roles in cell survival, evasion of apoptosis and proliferation (Domoto et al., 2016). CDKN2A has two different promoters, which are involved in the retinoblastoma protein (Rb) and p53 pathways, and acts as an inhibitor of CDK4/6 kinase. CDKN2A contributes to the regulation of cell cycle progression by inhibiting the S phase. It binds to CDK4/6, inhibiting cyclin D-CDK4/6 complex formation and phosphorylation of Rb family members so preventing exit from G1 phase of the cell cycle (Romagosa et al., 2011). HDAC2 protein plays role in gene transcription, DNA repair, immune stability, and related signaling pathways (Shan et al., 2017). Moreover, the PML protein showed a high expression pattern in endothelial cells and peripheral nerves/ganglia, and similarly, CDKN2A protein showed a high expression pattern in glandular cells and peripheral nerve/ganglion. When HDAC2 protein was examined in Human Protein Atlas, it showed a high expression pattern in glandular cells and an intermediate expression pattern in endothelial cells and peripheral nerve/ganglion. GSK3B protein also showed an intermediate expression pattern in glandular cells and peripheral nerve/ganglion and was detected in plasma (Uhlen et al., 2015).

Among the proposed drugs, abacavir, voriconazole, exemestane, and nortriptyline hydrochloride targeted the PML protein; whereas tolcapone, bromocriptine mesylate, and nortriptyline hydrochloride targeted GSK3B. CDKN2A was the target protein for ribociclib. Additionally, theophylline targeted the HDAC2 protein. Recently, researchers investigated the antitumor effects of the reverse transcription inhibitor abacavir in several studies. Carlini and coworkers (Carlini et al., 2010) investigated the effect of abacavir in prostate cancer cell lines and observed that abacavir reduced cell growth, migration, and invasion processes. Another study investigated the effect of abacavir and its combination with other reverse transcriptase inhibitors in breast cancer cells and found an increase in apoptosis and a decrease in migration in treated cells (Sekeroglu et al., 2021). On the other hand, Exemestane is already used in breast cancer, but Koutras and colleagues (Koutras et al., 2009) observed the antiplatelet effect of Exemestane also in lung cancer cells. Moreover, nortriptyline, a tricyclic antidepressant used to treat depression, however, has antineoplastic activity in various cancers such as bladder, prostate, myeloma (Pan et al., 2010; Mao et al., 2011; Yuan et al., 2015). In addition, recent studies have shown that tolcapone, a catechol-O-methyltransferase inhibitor used in Parkinson’s disease, decreases cell viability in lung cancer and neuroblastoma cell lines (Forester and Lambert, 2014; Maser et al., 2017).

According to DrugBank (Wishart et al., 2006), abacavir was a reverse transcriptase inhibitor and was used against Human Immunodeficiency Virus Type 1 (HIV-1). It targets HIV reverse transcriptase by forming the pharmacologically active compound carbovir 5′-triphosphate, which is an analogue of guanosine. Voriconazole is used to treat severe fungal infections. It binds to 14-alpha-sterol demethylase, which is known as CYP51, and inhibits the demethylation of lanosterol as part of the ergosterol synthesis pathway in yeast and other fungi. Exemestane is an irreversible steroidal aromatase inactivator structurally related to the natural substrate androstenedione that has been used to treat estrogen receptor-positive breast cancer. It acts as a false substrate for the aromatase enzyme and is processed into an intermediate that irreversibly binds to the active site of the enzyme and causes its inactivation. Nortriptyline hydrochloride is a tricyclic antidepressant and is used to treat depression by inhibiting serotonin and norepinephrine reuptake in neuronal cell membranes. It also exerts antimuscarinic effects through its action on the acetylcholine receptor. Tolcapone is a selective and reversible inhibitor of catechol-O-methyltransferase (COMT). It is used for the symptomatic treatment of Parkinson’s Disease besides levodopa/carbidopa therapy. Bromocriptine mesylate is a dopamine D2 receptor that stimulates centrally located dopaminergic receptors, resulting in a number of pharmacological effects. It is used for signs and symptoms of Parkinsonian syndrome. Ribociclib is a cyclin-dependent kinase inhibitor and helps slow the growth of cancer cells by inhibiting CDK4/6 by arresting cells at the G1 checkpoint, preventing tumor cells from proliferating. Theophylline is used to treat lung diseases such as asthma. It has two distinct effects in the airways of patients with reversible obstruction; smooth muscle relaxation and suppression of airway response to stimuli. Besides, it was investigated whether these drugs had been used previously in CRC. With the exception of theophylline, the drugs were novel in CRC. Theophylline was used as a positive control in this study.

Drugs that are candidates for the targets (PML, CDKN2A, GSK3B, and HDAC2) were further evaluated using the current literature and screened for association with CRC. Theophylline was used in SW480 cell lines and the IC₅₀ value of theophylline was found to be 10⁻⁴ M at 48 h in rectal cancer (Peng et al., 2018). In our study, theophylline showed its IC₅₀ effect at 500 µM concentration at 72 h after viability assay in the carcinoma cell line. The difference in cell lines might have caused the difference in inhibitory concentrations. It was also observed that ribociclib, a cyclin-dependent kinase inhibitor that slows cancer cell growth by inhibiting CDK4/6 proteins, significantly decreased viability in the healthy cell line in addition to its effect in the carcinoma cell line. We speculate that the possible reason for the selectivity between healthy and CRC cells may be that CRC cells are less dependent on the CDK4/6 axis for proliferation and survival than healthy cells.

Conclusion

This study aimed to find new therapeutic targets and drugs for the treatment of CRC. There were 24 common protein modules, whose Kaplan-Meier plots showed that six modules were prognostically important. In addition, functional enrichment analysis of the modules revealed that they were involved in signaling pathways related to cancer and metabolism. Molecular docking results showed that each potential drug and respective protein has significantly higher binding affinities compared to the positive control. Four of the drug candidates (abacavir, exemestane, nortriptyline hydrochloride, tolcapone) showed statistically significant inhibition profiles on the CRC cell line. The efficacy of four of these drugs on CRC cell lines was demonstrated in this study for the first time in the literature. Therefore, they can be used for further studies. The positive control (theophylline) showed that our methods can be used as new therapeutics in cancer therapy. Our results may provide new and complementary strategies for the treatment of CRC. On the other hand, these results should be supported by further experiments to elucidate the mechanisms of drug action in colorectal cancer cells. Moreover, our results highlight the value of studying DIPs to propose potential therapeutics.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: GSE8671-https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE8671; TCGA-COAD- https://portal.gdc.cancer.gov/projects/TCGA-COAD; TCGA-READ- https://portal.gdc.cancer.gov/projects/TCGA-READ.

Author Contributions

HB, BT, PO, and KA designed the data analysis framework. HB and GG performed the data analyses. SA and BK designed in vitro experiments. HB, GG, BK, PO, BT, and KA interpreted the results. KA conceived and directed the study. HB drafted the article. HB, BT, PO, BK and KA revised the article. All authors read and approved the final article.

Funding

This research was financially supported by TUBITAK through project 119S174, and Marmara University Research Fund (BAPKO) through project FEN-C-YLP-120619-0199.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The scholarships under the TUBITAK 2211-C and YOK 100/2000 Doctoral Fellowship Programs provided to GG, and the fellowship under the TUBITAK project 119S174 provided to HB are greatly acknowledged.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbinf.2021.710591/full#supplementary-material

References

Aguirre-Gamboa, R., Gomez-Rueda, H., Martínez-Ledesma, E., Martínez-Torteya, A., Chacolla-Huaringa, R., Rodriguez-Barrientos, A., et al. (2013). SurvExpress: An Online Biomarker Validation Tool and Database for Cancer Gene Expression Data Using Survival Analysis. PLoS One 8, e74250. doi:10.1371/journal.pone.0074250

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene Ontology: Tool for the Unification of Biology. The Gene Ontology Consortium. Nat. Genet. 25, 25–29. doi:10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayyildiz, D., Gov, E., Sinha, R., and Arga, K. Y. (2017). Ovarian Cancer Differential Interactome and Network Entropy Analysis Reveal New Candidate Biomarkers. OMICS 21, 285–294. doi:10.1089/omi.2017.0010

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2013). NCBI GEO: Archive for Functional Genomics Data Sets-Update. Nucleic Acids Res. 41, D991–D995. doi:10.1093/nar/gks1193

PubMed Abstract | CrossRef Full Text | Google Scholar

Berman, H. M., Battistuz, T., Bhat, T. N., Bluhm, W. F., Bourne, P. E., Burkhardt, K., et al. (2002). The Protein Data Bank. Acta Crystallogr. D Biol. Crystallogr. 58, 899–907. doi:10.1107/S0907444902003451

PubMed Abstract | CrossRef Full Text | Google Scholar

Biovia, D. S. (2016). Discovery Studio Modeling Environment. San Diego: Dassault Systèmes. Release 2017.

Buonfiglio, R., Prati, F., Bischetti, M., Cavarischia, C., Furlotti, G., and Ombrato, R. (2020). Discovery of Novel Imidazopyridine GSK-3β Inhibitors Supported by Computational Approaches. Molecules 25, 2163. doi:10.3390/molecules25092163

PubMed Abstract | CrossRef Full Text | Google Scholar

Calimlioglu, B., Karagoz, K., Sevimoglu, T., Kilic, E., Gov, E., and Arga, K. Y. (2015). Tissue-Specific Molecular Biomarker Signatures of Type 2 Diabetes: An Integrative Analysis of Transcriptomics and Protein-Protein Interaction Data. OMICS 19, 563–573. doi:10.1089/omi.2015.0088

PubMed Abstract | CrossRef Full Text | Google Scholar

Caliskan, A., Gulfidan, G., Sinha, R., and Arga, K. Y. (2021). Differential Interactome Proposes Subtype-specific Biomarkers and Potential Therapeutics in Renal Cell Carcinomas. Jpm 11, 158. doi:10.3390/jpm11020158

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlini, F., Ridolfi, B., Molinari, A., Parisi, C., Bozzuto, G., Toccacieli, L., et al. (2010). The Reverse Transcription Inhibitor Abacavir Shows Anticancer Activity in Prostate Cancer Cell Lines. PLoS One 5, e14221. doi:10.1371/journal.pone.0014221

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatr-aryamontri, A., Oughtred, R., Boucher, L., Rust, J., Chang, C., Kolas, N. K., et al. (2017). The BioGRID Interaction Database: 2017 Update. Nucleic Acids Res. 45, D369–D379. doi:10.1093/nar/gkw1102

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Liu, H., and Qing, G. (2018). Targeting Oncogenic Myc as a Strategy for Cancer Treatment. Signal. Transduct. Target. Ther. 3, 5–7. doi:10.1038/s41392-018-0008-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Chung, F. H., Chiang, Y. R., Tseng, A. L., Sung, Y. C., Lu, J., Huang, M. C., et al. (2014). Functional Module Connectivity Map (FMCM): A Framework for Searching Repurposed Drug Compounds for Systems Treatment of Cancer and an Application to Colorectal Adenocarcinoma. PLoS One 9, e86299. doi:10.1371/journal.pone.0086299

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang, C. V., Le, A., and Gao, P. (2009). MYC-induced Cancer Cell Energy Metabolism and Therapeutic Opportunities. Clin. Cancer Res. 15, 6479–6483. doi:10.1158/1078-0432.CCR-09-0889

PubMed Abstract | CrossRef Full Text | Google Scholar

Díaz-Carballo, D., Acikelli, A. H., Klein, J., Jastrow, H., Dammann, P., Wyganowski, T., et al. (2015). Therapeutic Potential of Antiviral Drugs Targeting Chemorefractory Colorectal Adenocarcinoma Cells Overexpressing Endogenous Retroviral Elements. J. Exp. Clin. Cancer Res. 34, 81–13. doi:10.1186/s13046-015-0199-5

CrossRef Full Text | Google Scholar

Domoto, T., Pyko, I. V., Furuta, T., Miyashita, K., Uehara, M., Shimasaki, T., et al. (2016). Glycogen Synthase Kinase-3β Is a Pivotal Mediator of Cancer Invasion and Resistance to Therapy. Cancer Sci. 107, 1363–1372. doi:10.1111/cas.13028

PubMed Abstract | CrossRef Full Text | Google Scholar

Douaiher, J., Ravipati, A., Grams, B., Chowdhury, S., Alatise, O., and Are, C. (2017). Colorectal Cancer-Global burden, Trends, and Geographical Variations. J. Surg. Oncol. 115, 619–630. doi:10.1002/jso.24578

CrossRef Full Text | Google Scholar

Douillard, J. Y., Cunningham, D., Roth, A. D., Navarro, M., James, R. D., Karasek, P., et al. (2000). Irinotecan Combined with Fluorouracil Compared with Fluorouracil Alone as First-Line Treatment for Metastatic Colorectal Cancer: a Multicentre Randomised Trial. Lancet 355, 1041–1047. doi:10.1016/S0140-6736(00)02034-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabregat, A., Sidiropoulos, K., Viteri, G., Forner, O., Marin-Garcia, P., Arnau, V., et al. (2017). Reactome Pathway Analysis: A High-Performance In-Memory Approach. BMC Bioinformatics 18, 142. doi:10.1186/s12859-017-1559-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Forester, S. C., and Lambert, J. D. (2014). Synergistic Inhibition of Lung Cancer Cell Lines by (-)-Epigallocatechin-3-Gallate in Combination with Clinically Used Nitrocatechol Inhibitors of Catechol-O-Methyltransferase. Carcinogenesis 35, 365–372. doi:10.1093/carcin/bgt347

PubMed Abstract | CrossRef Full Text | Google Scholar

Gov, E., and Arga, K. Y. (2017). Differential Co-expression Analysis Reveals a Novel Prognostic Gene Module in Ovarian Cancer. Sci. Rep. 7, 4996–5010. doi:10.1038/s41598-017-05298-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Guan, D., and Kao, H. Y. (2015). The Function, Regulation and Therapeutic Implications of the Tumor Suppressor Protein, PML. Cell Biosci 5, 60–14. doi:10.1186/s13578-015-0051-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulfidan, G., Turanli, B., Beklen, H., Sinha, R., and Arga, K. Y. (2020). Pan-cancer Mapping of Differential Protein-Protein Interactions. Sci. Rep. 10, 3272. doi:10.1038/s41598-020-60127-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Higurashi, T., and Nakajima, A. (2018). Metformin and Colorectal Cancer. Front. Endocrinol. 9, 1–7. doi:10.3389/fendo.2018.00622

PubMed Abstract | CrossRef Full Text | Google Scholar

Ilyas, U., Zaman, S. U., Altaf, R., Nadeem, H., and Muhammad, S. A. (2020). Genome Wide Meta-Analysis of cDNA Datasets Reveals New Target Gene Signatures of Colorectal Cancer Based on Systems Biology Approach. J. Biol. Res. (Thessalon) 27, 8–13. doi:10.1186/s40709-020-00118-1

CrossRef Full Text | Google Scholar

Kamburov, A., Stelzl, U., Lehrach, H., and Herwig, R. (2013). The ConsensusPathDB Interaction Database: 2013 Update. Nucleic Acids Res. 41, D793–D800. doi:10.1093/nar/gks1055

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y., and Morishima, K. (2017). KEGG: New Perspectives on Genomes, Pathways, Diseases and Drugs. Nucleic Acids Res. 45, D353–D361. doi:10.1093/nar/gkw1092

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., et al. (2019). PubChem 2019 Update: Improved Access to Chemical Data. Nucleic Acids Res. 47, D1102–D1109. doi:10.1093/nar/gky1033

PubMed Abstract | CrossRef Full Text | Google Scholar

Koh, S. L., Ager, E. I., Costa, P. L., Malcontenti-Wilson, C., Muralidharan, V., and Christophi, C. (2014). Blockade of the Renin-Angiotensin System Inhibits Growth of Colorectal Cancer Liver Metastases in the Regenerating Liver. Clin. Exp. Metastasis 31, 395–405. doi:10.1007/s10585-014-9635-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kori, M., and Yalcin Arga, K. (2018). Potential Biomarkers and Therapeutic Targets in Cervical Cancer: Insights from the Meta-Analysis of Transcriptomics Data within Network Biomedicine Perspective. PLoS One 13, e0200717. doi:10.1371/journal.pone.0200717

PubMed Abstract | CrossRef Full Text | Google Scholar

Koutras, A., Giannopoulou, E., Kritikou, I., Antonacopoulou, A., Evans, T. R., Papavassiliou, A. G., et al. (2009). Antiproliferative Effect of Exemestane in Lung Cancer Cells. Mol. Cancer 8, 109. doi:10.1186/1476-4598-8-109

PubMed Abstract | CrossRef Full Text | Google Scholar

Lamb, J., Crawford, E. D., Peck, D., Modell, J. W., Blat, I. C., Wrobel, M. J., et al. (20061929–1935). The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease. Science 313, 1929–1935. doi:10.1126/science.1132939

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao, X., Hou, T., Cao, B., Wang, W., Li, Z., Chen, S., et al. (2011). The Tricyclic Antidepressant Amitriptyline Inhibits D-Cyclin Transactivation and Induces Myeloma Cell Apoptosis by Inhibiting Histone Deacetylases: In Vitro and In Silico Evidence. Mol. Pharmacol. 79, 672–680. doi:10.1124/mol.110.068122

PubMed Abstract | CrossRef Full Text | Google Scholar

Maser, T., Rich, M., Hayes, D., Zhao, P., Nagulapally, A. B., Bond, J., et al. (2017). Tolcapone Induces Oxidative Stress Leading to Apoptosis and Inhibition of Tumor Growth in Neuroblastoma. Cancer Med. 6, 1341–1352. doi:10.1002/cam4.1065

PubMed Abstract | CrossRef Full Text | Google Scholar

Massó-Vallés, D., and Soucek, L. (2020). Blocking Myc to Treat Cancer: Reflecting on Two Decades of Omomyc. Cells 9, 883. doi:10.3390/cells9040883

CrossRef Full Text | Google Scholar

Mauricio, M. M. (2019). Modern Strategies in Cancer Study: Drug Repositioning in Colorectal Cancer Treatment. Nacs 2, 212–216. doi:10.31031/nacs.2019.02.000544

CrossRef Full Text | Google Scholar

Neugut, A. I., Lin, A., Raab, G. T., Hillyer, G. C., Keller, D., O'Neil, D. S., et al. (2019). FOLFOX and FOLFIRI Use in Stage IV Colon Cancer: Analysis of SEER-Medicare Data. Clin. Colorectal Cancer 18, 133–140. doi:10.1016/j.clcc.2019.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishimura, D. (2001). BioCarta. Biotech. Softw. Internet Rep. 2, 117–120. doi:10.1089/152791601750294344

CrossRef Full Text | Google Scholar

Nowak-Sliwinska, P., Scapozza, L., and Ruiz i Altaba, A. (2019). Drug Repurposing in Oncology: Compounds, Pathways, Phenotypes and Computational Approaches for Colorectal Cancer. Biochim. Biophys. Acta Rev. Cancer 1871, 434–454. doi:10.1016/j.bbcan.2019.04.0051871

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, C.-C., Shaw, C.-F., Huang, J.-K., Kuo, C.-C., Kuo, D.-H., Shieh, P., et al. (2010). Effect of Nortriptyline on Cytosolic Ca2+ Regulation and Viability in PC3 Human Prostate Cancer Cells. Drug Dev. Res. 71, 323–330. doi:10.1002/ddr.20377

CrossRef Full Text | Google Scholar

Park, S. M., Hwang, C. Y., Cho, S. H., Lee, D., Gong, J. R., Lee, S., et al. (2019). Systems Analysis Identifies Potential Target Genes to Overcome Cetuximab Resistance in Colorectal Cancer Cells. FEBS J. 286, 1305–1318. doi:10.1111/febs.14773

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, H., Su, Q., Lin, Z. C., Zhu, X. H., Peng, M. S., and Lv, Z. B. (2018). Potential Suppressive Effects of Theophylline on Human Rectal Cancer SW480 Cells In Vitro by Inhibiting YKL-40 Expression. Oncol. Lett. 15, 7403–7408. doi:10.3892/ol.2018.8220

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahib, L., Smith, B. D., Aizenberg, R., Rosenzweig, A. B., Fleshman, J. M., and Matrisian, L. M. (2014). Projecting Cancer Incidence and Deaths to 2030: The Unexpected burden of Thyroid, Liver, and Pancreas Cancers in the united states. Cancer Res. 74, 2913–2921. doi:10.1158/0008-5472.CAN-14-0155

PubMed Abstract | CrossRef Full Text | Google Scholar

Romagosa, C., Simonetti, S., López-Vicente, L., Mazo, A., Lleonart, M. E., Castellvi, J., et al. (2011). p16(Ink4a) Overexpression in Cancer: a Tumor Suppressor Gene Associated with Senescence and High-Grade Tumors. Oncogene 30, 2087–2097. doi:10.1038/onc.2010.614

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabates-Bellver, J., Van Der Flier, L. G., De Palo, M., Cattaneo, E., Maake, C., Rehrauer, H., et al. (2007). Transcriptome Profile of Human Colorectal Adenomas. Mol. Cancer Res. 5, 1263–1275. doi:10.1158/1541-7786.MCR-07-0267

PubMed Abstract | CrossRef Full Text | Google Scholar

Safari-Alighiarloo, N., Taghizadeh, M., Rezaei-Tavirani, M., Goliaei, B., and Peyvandi, A. A. (2014). Protein-protein Interaction Networks (PPI) and Complex Diseases. Gastroenterol. Hepatol. Bed Bench 7, 17–31. doi:10.22037/ghfbb.v7i1.511

PubMed Abstract | CrossRef Full Text | Google Scholar

Şekeroğlu, Z. A., Şekeroğlu, V., and Küçük, N. (2021). Effects of Reverse Transcriptase Inhibitors on Proliferation, Apoptosis, and Migration in Breast Carcinoma Cells. Int. J. Toxicol. 40, 52–61. doi:10.1177/1091581820961498

CrossRef Full Text | Google Scholar

Senga, S. S., and Grose, R. P. (2021). Hallmarks of Cancer-The New Testament. Open Biol. 11, 200358. doi:10.1098/rsob.200358

PubMed Abstract | CrossRef Full Text | Google Scholar

Sevimoglu, T., and Arga, K. Y. (2014). The Role of Protein Interaction Networks in Systems Biomedicine. Comput. Struct. Biotechnol. J. 11, 22–27. doi:10.1016/j.csbj.2014.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Sevimoglu, T., Turanli, B., Bereketoglu, C., Arga, K. Y., and Karadag, A. S. (2018). Systems Biomarkers in Psoriasis: Integrative Evaluation of Computational and Experimental Data at Transcript and Protein Levels. Gene 647, 157–163. doi:10.1016/j.gene.2018.01.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Shan, W., Jiang, Y., Yu, H., Huang, Q., Liu, L., Guo, X., et al. (2017). HDAC2 Overexpression Correlates with Aggressive Clinicopathological Features and DNA-Damage Response Pathway of Breast Cancer. Am. J. Cancer Res. 7, 1213–1226.

Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Shim, J. S., and Liu, J. O. (2014). Recent Advances in Drug Repositioning for the Discovery of New Anticancer Drugs. Int. J. Biol. Sci. 10, 654–663. doi:10.7150/ijbs.9224

CrossRef Full Text | Google Scholar

Smith, J. J., Deane, N. G., Wu, F., Merchant, N. B., Zhang, B., Jiang, A., et al. (2010). Experimentally Derived Metastasis Gene Expression Profile Predicts Recurrence and Death in Patients with Colon Cancer. Gastroenterology 138, 958–968. doi:10.1053/j.gastro.2009.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinformatics 54. doi:10.1002/cpbi.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A. Cancer J. Clin. 71, 209–249. doi:10.3322/caac.21660

CrossRef Full Text | Google Scholar

Tomczak, K., Czerwińska, P., and Wiznerowicz, M. (2015). The Cancer Genome Atlas (TCGA): An Immeasurable Source of Knowledge. Contemp. Oncol. (Pozn) 19, A68–A77. doi:10.5114/wo.2014.47136

PubMed Abstract | CrossRef Full Text | Google Scholar

Trott, O., and Olson, A. J. (2019). AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem. 31, 455–461. doi:10.1002/jcc.21334.AutoDock

CrossRef Full Text | Google Scholar

Turanli, B., Gulfidan, G., and Arga, K. Y. (2017b). Transcriptomic-Guided Drug Repositioning Supported by a New Bioinformatics Search Tool: GeneXpharma. OMICS 21, 584–591. doi:10.1089/omi.2017.0127

PubMed Abstract | CrossRef Full Text | Google Scholar

Turanli, B., Karagoz, K., Bidkhori, G., Sinha, R., Gatza, M. L., Uhlen, M., et al. (2019b). Multi-omic Data Interpretation to Repurpose Subtype Specific Drug Candidates for Breast Cancer. Front. Genet. 10, 420. doi:10.3389/fgene.2019.00420

PubMed Abstract | CrossRef Full Text | Google Scholar

Turanli, B., Altay, O., Borén, J., Turkez, H., Nielsen, J., Uhlen, M., et al. (2021a). Systems Biology Based Drug Repositioning for Development of Cancer Therapy. Semin. Cancer Biol. 68, 47–58. doi:10.1016/j.semcancer.2019.09.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Turanli, B., Gulfidan, G., and Arga, K. Y. (2017a). “Systems Biomedicine : From Genome Scale Molecular Networks to Efficacious Biomarkers and Therapeutic Target,” in In Researches on Science and Art in 21st Century Turkey. Editors H. Arapligirlioglu, A. Atik, R. Elliott, and E. Turgeon (Ankara: Gece Kitaplığı), 2380.

Google Scholar

Uhlén, M., Fagerberg, L., Hallström, B. M., Lindskog, C., Oksvold, P., Mardinoglu, A., et al. (2015). Proteomics. Tissue-Based Map of the Human Proteome. Science 347, 1260419. doi:10.1126/science.1260419

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Noort, V., Schölch, S., Iskar, M., Zeller, G., Ostertag, K., Schweitzer, C., et al. (2014). Novel Drug Candidates for the Treatment of Metastatic Colorectal Cancer through Global Inverse Gene-Expression Profiling. Cancer Res. 74, 5690–5699. doi:10.1158/0008-5472.CAN-13-3540

PubMed Abstract | CrossRef Full Text | Google Scholar

Velankar, S., van Ginkel, G., Alhroub, Y., Battle, G. M., Berrisford, J. M., Conroy, M. J., et al. (2016). PDBe: Improved Accessibility of Macromolecular Structure Data from PDB and EMDB. Nucleic Acids Res. 44, D385–D395. doi:10.1093/nar/gkv1047

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P., Benhenda, S., Wu, H., Lallemand-Breitenbach, V., Zhen, T., Jollivet, F., et al. (2018). RING Tetramerization Is Required for Nuclear Body Biogenesis and PML Sumoylation. Nat. Commun. 9, 1277. doi:10.1038/s41467-018-03498-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wishart, D. S., Knox, C., Guo, A. C., Shrivastava, S., Hassanali, M., Stothard, P., et al. (2006). DrugBank: a Comprehensive Resource for In Silico Drug Discovery and Exploration. Nucleic Acids Res. 34, D668–D672. doi:10.1093/nar/gkj067

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolf, E., and Eilers, M. (2020). Targeting MYC Proteins for Tumor Therapy. Annu. Rev. Cancer Biol. 4, 61–75. doi:10.1146/annurev-cancerbio-030518-055826

CrossRef Full Text | Google Scholar

Würth, R., Thellung, S., Bajetto, A., Mazzanti, M., Florio, T., and Barbieri, F. (2016). Drug-repositioning Opportunities for Cancer Therapy: Novel Molecular Targets for Known Compounds. Drug Discov. Today 21, 190–199. doi:10.1016/j.drudis.2015.09.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, H., Li, J., Xie, H., and Wang, Y. (2018). Review of Drug Repositioning Approaches and Resources. Int. J. Biol. Sci. 14, 1232–1244. doi:10.7150/ijbs.24612

PubMed Abstract | CrossRef Full Text | Google Scholar

Yella, J. K., Yaddanapudi, S., Wang, Y., and Jegga, A. G. (2018). Changing Trends in Computational Drug Repositioning. Pharmaceuticals (Basel) 11, 57. doi:10.3390/ph11020057

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, W., Liu, J., Yu, Y., Zhang, V., Clausen, D., Kelly, J., et al. (2020). Discovery of Ethyl Ketone-Based HDACs 1, 2, and 3 Selective Inhibitors for HIV Latency Reactivation. Bioorg. Med. Chem. Lett. 30, 127197. doi:10.1016/j.bmcl.2020.127197

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, C., Selby, T. L., Li, J., Byeon, I. J., and Tsai, M. D. (2000). Tumor Suppressor INK4: Refinement of p16INK4A Structure and Determination of p15INK4B Structure by Comparative Modeling and NMR Data. Protein Sci. 9, 1120–1128. doi:10.1110/ps.9.6.1120

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, S. Y., Cheng, C. L., Ho, H. C., Wang, S. S., Chiu, K. Y., Su, C. K., et al. (2015). Nortriptyline Induces Mitochondria and Death Receptor-Mediated Apoptosis in Bladder Cancer Cells and Inhibits Bladder Tumor Growth In Vivo. Eur. J. Pharmacol. 761, 309–320. doi:10.1016/j.ejphar.2015.06.007

CrossRef Full Text | Google Scholar

Zhang, Q., Wang, J. Y., Zhou, S. Y., Yang, S. J., and Zhong, S. L. (2019). Circular RNA Expression in Pancreatic Ductal Adenocarcinoma. Oncol. Lett. 18, 2923–2930. doi:10.3892/ol.2019.10624

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, R., Choi, B. Y., Lee, M. H., Bode, A. M., and Dong, Z. (2016). Implications of Genetic and Epigenetic Alterations of CDKN2A (p16(INK4a)) in Cancer. EBioMedicine 8, 30–39. doi:10.1016/j.ebiom.2016.04.017

PubMed Abstract | CrossRef Full Text | Google Scholar