Bioinformatics and network-based screening and discovery of potential molecular targets and small molecular drugs for breast cancer

Abstract

Accurate identification of molecular targets of disease plays an important role in diagnosis, prognosis, and therapies. Breast cancer (BC) is one of the most common malignant cancers in women worldwide. Thus, the objective of this study was to accurately identify a set of molecular targets and small molecular drugs that might be effective for BC diagnosis, prognosis, and therapies, by using existing bioinformatics and network-based approaches. Nine gene expression profiles (GSE54002, GSE29431, GSE124646, GSE42568, GSE45827, GSE10810, GSE65216, GSE36295, and GSE109169) collected from the Gene Expression Omnibus (GEO) database were used for bioinformatics analysis in this study. Two packages, LIMMA and clusterProfiler, in R were used to identify overlapping differential expressed genes (oDEGs) and significant GO and KEGG enrichment terms. We constructed a PPI (protein–protein interaction) network through the STRING database and identified eight key genes (KGs) EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5 by using six topological measures, betweenness, closeness, eccentricity, degree, MCC, and MNC, in the Analyze Network tool in Cytoscape. Three online databases GSCALite, Network Analyst, and GEPIA were used to analyze drug enrichment, regulatory interaction networks, and gene expression levels of KGs. We checked the prognostic power of KGs through the prediction model using the popular machine learning algorithm support vector machine (SVM). We suggested four TFs (TP63, MYC, SOX2, and KDM5B) and four miRNAs (hsa-mir-16-5p, hsa-mir-34a-5p, hsa-mir-1-3p, and hsa-mir-23b-3p) as key transcriptional and posttranscriptional regulators of KGs. Finally, we proposed 16 candidate repurposing drugs YM201636, masitinib, SB590885, GSK1070916, GSK2126458, ZSTK474, dasatinib, fedratinib, dabrafenib, methotrexate, trametinib, tubastatin A, BIX02189, CP466722, afatinib, and belinostat for BC through molecular docking analysis. Using BC cell lines, we validated that masitinib inhibits the mTOR signaling pathway and induces apoptotic cell death. Therefore, the proposed results might play an effective role in the treatment of BC patients.

Introduction

Breast cancer (BC) is the most common cancer diagnosis and the leading cause of cancer-related deaths in women worldwide. It was calculated that 2,261,419 (11.7% of all cancers) new cases and 684,996 (6.9% of all cancers) deaths occurred in BC in 2020 (Sung et al., 2021). In China, 416,371 new cases and 117,414 BC-related deaths occurred in 2020, accounting for 18.4% of diagnoses and 17.1% of deaths (Cao et al., 2021). It was also estimated that there were 1,700,000 new cases and 521,900 deaths in 2012 BC worldwide (Torre et al., 2015). We noticed that BC has a high prevalence and mortality rate and is steadily increasing every year, although it is claimed that there have been significant advances in systematic treatment over the decades. The 5-year overall survival rate for BC patients is still low, although it depends on various factors (Sarveazad et al., 2018). Therefore, the progress of the existing treatment has not yet reached a satisfactory level. Therefore, it is urgent to improve research to discover potential molecular targets and effective candidate drugs for the development of innovative therapies for BC.

Novel drug discovery is challenging, time consuming, and expensive due to the very low rate of approval through clinical trials. Discovering a new drug using the de novo technique takes approximately 14 years and costs 800 million dollars, yet many pharmaceutical companies are working on it (Song et al., 2009; Lavecchia and Di Giovanni, 2013). Drug repurposing (DR) is a promising strategy to identify new indications for a specific disease by using approved drugs (existing drugs) (Simsek et al., 2018). The DR strategy is safer, cheaper, and less time consuming than the de novo strategy due to the knowledge of in vivo screening, chemical optimization, and toxicology of existing drugs (Ko, 2020). Thus, the DR technique in computers widely used by pharmaceutical companies and researchers over the past few decades has achieved significant success (Shi et al., 2020; Alam et al., 2022b). Molecular docking analysis is a momentous strategy for validating drug-target structural binding performance in computational DR processes. Bioinformatics analysis plays a significant role in accurately identifying the key genes (KGs)/targets with candidate drugs of the disease to inform diagnosis, prognosis, and therapies. The gene expression profile analysis is one of the most popular platforms for disease-guided KG identification. Thus, we used multiple datasets collected from different environments to identify more common and stable BC-guided KG through gene expression profiles (Alam et al., 2022a).

In the present study, we identified eight stable BC-guided KGs (EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5) and highlighted their role as molecular targets. We also proposed 16 KG-guided candidate drugs (YM201636, masitinib, SB590885, GSK1070916, GSK2126458, ZSTK474, dasatinib, fedratinib, dabrafenib, methotrexate, trametinib, tubastatin A, BIX02189, CP466722, Afatinib, and Belinostat) for BC treatment and validated them in silico through molecular docking analysis. Using BC cell lines, we found that masitinib inhibits the mTOR signaling pathway and induces apoptotic cell death.

Materials and methods

Microarray data

The gene expression profiles analyzed in our study were downloaded from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) (Edgar and Barrett, 2006). The workflow of the present study is presented in Supplementary Figure S1. Eight sets of gene expression profiles with accession numbers GSE54002 (Tan et al., 2014), GSE29431 (Wang et al., 2021), GSE124646 (Sinn et al., 2019), GSE42568 (Clarke et al., 2013), GSE45827 (Gruosso et al., 2016), GSE10810 (Pedraza et al., 2010), GSE65216 (Maubant et al., 2015), and GSE36295 (Karim et al., 2016) were used to identify differentially expressed genes (DEGs) between BC and normal samples. A total of 864 BC samples and 123 normal samples were included in these eight sets of profiles. An independent gene expression profile (GSE109169) (Yu et al., 2020) including 25 BC and 25 normal samples was collected to investigate the prognostic performance of KGs through the cancer prediction model. Furthermore, the six profiles (GSE54002, GSE29431, GSE42568, GSE45827, GSE10810, and GSE65216) based on the GPL570 platform, profile GSE124646 based on the GPL96 platform, profile GSE36295 based on the GPL6244 platform, and profile GSE109169 based on the GPL96 platform were identified. Details of the datasets are described in Table 1.

Data processing and differentially expressed gene identification

We normalized all datasets by using log2-transformation, after which we used the normalized data for further analysis. Then, the Limma (version: 3.14) package (Ritchie et al., 2015) in RStudio (version: 1.2.5019) was used to identify DEGs between BC samples and normal samples for eight datasets individually. The moderated T-statistics and the Benjamini and Hochberg false discovery rate method were used to calculate the p value and adjusted p value (adj.P.Val) and logFC (fold change), accordingly. The cutoff criteria adj.P.Val < 0.01 and | logFC | > 1.5 were considered to select significant DEGs. Overlapping upregulated and downregulated DEGs were screened using the Venn diagram web tool (https://bioinformatics.psb.ugent.be/webtools/Venn/).

Protein–protein interaction network analysis of overlapping differentially expressed genes to identify key genes

The PPI network of oDEGs was constructed using the online database STRING (Search Tool for the Retrieval of Interacting Genes) (v 11.5) and was visualized through the Cytoscape software (Shannon et al., 2003; Szklarczyk et al., 2019). Then, we screened the top-ranked common KGs by using the six topological measures betweenness, closeness, eccentricity, degree, maximal clique centrality (MCC), and maximum neighborhood component (MNC) through the Analyze Network tool in the Cytoscape.

Functional and pathway enrichment analysis of differentially expressed genes

GO (Gene Ontology) function in three categories BPs (Biological Processes), CC (Cellular Component), and MF (MF) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyzed through the R/Bioconductor package “clusterProfiler” (version: 3.14). p-value cutoff = 0.05 and q-value cutoff = 0.05 were selected as the cutoff criteria (Yu et al., 2012).

Regulatory interaction network analysis

The regulatory interaction network TFs (Transcription factors)-KGs-miRNAs (microRNAs) to detect transcriptional and posttranscriptional regulatory factors of KGs. The Network Analyst (Zhou et al., 2019) web-based tool was used to construct the regulatory interaction networks and was visualized by using the Cytoscape software (Shannon et al., 2003).

Validation of prognostic power and expression pattern of key genes

We validated the expression pattern of KGs through TCGA (The Cancer Genome Atlas) RNA-seq data (independent data) by using the GEPIA database (Tang et al., 2017). Then, we performed the cancer prediction/classification model using the popular machine learning algorithm SVM to test the prognostic power of our proposed KGs. We used the expression profiles of KGs from eight sets of data (GSE54002, GSE29431, GSE124646, GSE42568, GSE45827, GSE10810, GSE65216, and GSE36295) as the prediction set and another independent validation set GSE26964. The ROCR package in R was used to generate the ROC curves (Sing et al., 2005).

Drug repurposing by molecular docking analysis

To explore KG-guided candidate drugs for BC by molecular docking analysis, we considered eight KGs as drug target proteins (receptors). Subsequently, the 77 KG-associated meta-drug agents (ligands) were selected from the GSCALite database on the basis of a positive Spearman rank correlation with at least one KG (Liu et al., 2018). We downloaded 3-D structures of target proteins and drug agents for molecular docking analysis. The 3-D structures of the eight KGs EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5 were downloaded from the Protein Data Bank (PDB) database with the PDB IDs 3g5z, 2haz, 4mi0, 6hyg, 6gu6, 3dj5, 1zxm, and 1xox, respectively (Berman et al., 2002). The 3D structures of the meta-drug agents were downloaded from the PubChem database (Kim et al., 2019). The 3-D structure of receptors and ligands was preprocessed for molecular docking analysis using the PyMol software. Molecular docking analysis was performed using the AutoDock Vina software in PyRx to check the structural binding performance between receptors and ligands and computed the binding affinity scores (BAS) (kcal/mol) (Trott and Olson, 2010; Dallakyan and Olson, 2015). The “plot.matrix” package in R was used to visualize the molecular docking results.

Statistical analysis

All gene expression profiles were statistically analyzed by the R software (Version 4.0.4). Moderated t-statistics were utilized to testH₀ (equally expressed gene (EEG) in both the case and control groups) versus H₁ (differentially expressed gene (DEG) between the case and control groups). Statistically significant differences were considered when p < 0.05.

Cell culture and drug treatment

The BC cell lines MCF-7 and MDA-MB-231 cells were cultured in Dulbecco’s modified Eagle’s medium (Gibco, Los Angeles, CA) containing 10% fetal bovine serum (FBS, Gibco) supplemented with penicillin (100  µ/ml) and streptomycin (100 μg/ml, Gibco). Masitinib was purchased from MedChemExpress (Monmouth Junction, NJ, United States) and dissolved in dimethyl sulfoxide (DMSO). The cells were treated with different doses of masitinib for 24 or 48 h. After treatment, the cells were harvested with lysis buffer for immunoblot analyses or fixed with 4% paraformaldehyde for immunofluorescent staining.

Cell viability assays

The cells were seeded in 96-well plates and cultured overnight. After treatment with masitinib at the indicated dose, a 10 µl solution of cell counting kit-8 (CCK-8) reagent was added to each well and incubated for 1 h in the dark. The absorbance was measured at 450 nm to calculate cell viability.

Propidium iodide staining assay

The cells were treated with masitinib at different concentrations. After treatment, the cells were incubated with Hoechst 33342 (Sigma, St. Louis, MO, United States) and PI (Sigma, St. Luis, MO, United States ) for 10 min. The cells were then imaged with an inverted IX71 microscope system (Olympus, Tokyo, Japan).

Immunoblot analyses

Immunoblot analyses were performed as described previously (Zhang et al., 2022). In brief, the cells were lysed in 1 × SDS lysis buffer supplemented with a protease inhibitor cocktail (Roche, Basel, Switzerland). The cell lysates were subjected to SDS–PAGE and transferred onto a PVDF membrane (Millipore, Billerica, MA, United States). The membranes were incubated with the following primary antibodies: anti-β-actin, anti-cleaved caspase-3, anti-cyclin D, anti-mTOR, anti-PARP, and anti-phospho-mTOR. The secondary antibodies, sheep anti-rabbit or anti-mouse IgG-HRP, were obtained from Thermo Fisher (Waltham, MA, United States). The proteins were visualized using an ECL detection kit (Thermo Fisher, Waltham, MA, United States).

Results

Identification of differentially expressed genes

We used the LIMMA statistical approach to identify DEGs and considered the threshold adj.P.Val < 0.01 and |log2 fold change (FC) | >1.5 for selecting significantly upregulated and downregulated DEGs. The DEGs in each dataset are presented using volcano plots (Figure 1A), where red dots indicate downregulated DEGs and green dots indicate upregulated DEGs. Then, we generated a Venn diagram of the overlapping upregulated and downregulated DEGs in (Figure 1B). Finally, we obtained a total of 68 oDEGs (36 upregulated and 32 downregulated oDEGs) in Table 2, which are referred to as DEGs in this study.

FIGURE 1

Protein–protein interaction network of overlapping differentially expressed genes to identify key genes

We visualized the PPI network of oDEGs in Figure 2, where green indicates upregulated oDEGs, pink indicates downregulated oDEGs, and large size indicates KGs. After that, we used six network scoring measures, including betweenness, closeness, eccentricity, degree, MCC, and MNC in the Analyze Network tool in the Cytoscape, and selected the top 25 genes for each measure. We extracted eight common genes (EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5) among the six lists and considered them as KGs in this study (Table 3). Further analyzed-based information of the KGs is presented in Table 4.

FIGURE 2

Functional and pathway enrichment analysis of overlapping differentially expressed genes

The GO functional and KEGG pathway enrichment analysis showed that 148 GO-BP terms, 27 GO-CC terms, 12 GO-MF terms, and 13 KEGG terms were enriched by oDEGs, and the top 10 terms of each category are represented in Figures 3A–C. Among them four significantly enriched BP terms, six CC terms, three MF terms, and five KEGG terms were directly associated with at least two KGs. The four BP terms were nuclear division, organelle fission, mitotic nuclear division, and chromosome segregation. The six CC terms were spindle, condensed chromosome, “chromosome, centromeric region”, chromosomal region, midbody, and spindle microtubule. The three MF terms were histone kinase activity, heparin binding, and protein serine/threonine kinase activity. Finally, the five KEGG terms were identified as progesterone-mediated oocyte maturation, focal adhesion, oocyte meiosis, adherens junction, and bacterial invasion of epithelial cells.

FIGURE 3

Regulatory interaction network analysis

Regulatory interaction networks (TFs-KGs-miRNAs) are visualized in Figure 4, where green indicates miRNAs, pink indicates KGs, blue indicates TFs, and large size indicates key factors. We observed that the four TFs, such as TP63, MYC, SOX2, and KDM5B are associated with all KGs; therefore, we considered these four TFs to be key transcriptional regulators of KGs. Similarly, the four miRNAs, including hsa-mir-16-5p, hsa-mir-34a-5p, hsa-mir-1-3p, and hsa-mir-23b-3p were also associated with all KGs; thus, we considered these four miRNAs to be key posttranscriptional regulators of KGs.

FIGURE 4

Validation of the expression pattern of key genes

The expression patterns of KGs are displayed by box plots for independent data (TCGA RNA-seq data) in Figure 5. We observed that the expression patterns of all KGs were highly differentiated; among them, six KGs (FN1, EZH2, CDK1, AURKA, TOP2A, and BIRC5) were upregulated, and the remaining two KGs (EGFR and MET) were downregulated, which supported our original results.

FIGURE 5

Prognostic power analysis of key genes

The ROC curves of the cancer prediction/classification model to test the prognostic power of our proposed KGs are presented in Figure 6. The AUC ranged from 94.7 to 99.7 for the prediction set (red color) and 85.5 to 93.5 for the validation set (blue color). Overall, EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5 achieved the best performance (AUC > 85.5) for each of the training and independent test datasets, which indicates that there is a strong prognostic power of the identified KGs for discriminating between the tumor and the normal samples in BC patients.

FIGURE 6

Exploring repurposing candidate drugs by molecular docking analysis

To explore the candidate drugs for BC, we considered eight KG-based proteins EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5 as drug targets. We collected eight KG-associated drugs from the GSCALite database and considered them meta-drug agents. Then, we performed molecular docking analysis between our proposed receptors and the meta-drug agents. The binding affinity score matrix between the ordered receptors and the ordered drug agents is displayed in Figure 7. We observed that the top-order eleven leading compounds/drugs (YM201636, masitinib, SB590885, GSK1070916, GSK2126458, ZSTK474, dasatinib, TG101348, dabrafenib, methotrexate, and trametinib) produced highly significant binding affinity scores (BAS) < −7 with all target proteins. Next, six drugs (tubastatin_A, lapatinib, BIX02189, CP466722, afatinib, and belinostat) had highly significant BAS < −7 with seven target proteins. The drug laptinib was approved by the FDA in 2007 for BC. Therefore, we suggest in this computational study that 16 repurposing candidate drugs may be effective against BC patients (Table 5).

FIGURE 7

Induction of breast cancer cell death by the candidate drug masitinib

To further confirm the potential roles of the candidate drugs on BC, we selected one of the candidates masitinib to examine the effects on two BC cell lines. In both MCF-7 (Figure 8A) and MDA-MB-231 (Figure 8B) cells, masitinib decreased cell viability in a dose-dependent manner. Moreover, masitinib inhibited the phosphorylation of mTOR (Figure 8C). Furthermore, masitinib decreased the levels of cyclin D, an essential regulator of the G1 to S phase transition and increased the cleavage of PARP1 (Figure 8D), suggesting that masitinib influences the cell cycle and induces apoptosis. The induction of apoptotic cell death by masitinib was further confirmed with PI staining (Figure 8E) and immunofluorescent staining using anti-cleaved caspase-3 antibodies (Figure 8F), showing that with higher doses of masitinib, more cells were labeled with PI and cleaved caspase-3 in both MCF-7 and MDA-MB-231 cells.

FIGURE 8

Discussion

BC is the cause of a high prevalence and leading mortality rate for women worldwide and is constantly increasing every year. The 5-year overall survival rate for BC patients is still very low, although there has been significant progress in systematic treatment over the past few decades. Therefore, research needs to be improved to discover potential molecular targets and effective candidate drugs that could play a significant role in improving survival rates and reducing mortality in BC patients.

In this study, we identified the potential biomarkers and the candidate drugs by highlighting their pathogenetic processes, transcriptional and posttranscriptional regulatory factors, expression levels, prognostic power, and drug molecules through integrated bioinformatics and network-based techniques. First, we screened 68 oDEGs between BC and normal samples using gene expression profiles. Then, we selected BC-causing two downregulated KGs (EGFR and MET) and six upregulated KGs (FN1, EZH2, CDK1, AURKA, TOP2A, and BIRC5) through the PPI network analysis (Figure 2; Table 3). Some studies have also proposed our identified KGs as BC-causing genes (Navolanic et al., 2003; Gastaldi et al., 2010; Li et al., 2017; Komoto et al., 2018; Deng et al., 2019). In particular, EGFR and its downstream signaling pathways are related to the progression of BC and play an important role against BC treatment using cytotoxic drugs (Navolanic et al., 2003). Analysis of gene expression profiles revealed that MET oncogene immunoreactivity is significantly higher in the progression of basal-like BC in humans than in other types of cancer (Gastaldi et al., 2010). Decreased stability of EZH2 is responsible for the progression and metastasis of BC, but ANCR plays an important role in controlling the stability of EZH2 which inhibits the progression of BC (Li et al., 2017). Bioinformatic analysis has proposed CDK1 as an upregulated gene related to BC progression and tumorigenesis (Deng et al., 2019). Chalcones inhibited the expression of AURKA protein and affected the anti-resistance and anti-metastatic properties of the BC cell lines MCF-7 and BT-20 (Komoto et al., 2018). Based on in silico analysis proposed that TOP2A is related to BC tumorigenesis and progression and has been highlighted as a pathogenetic process (Deng et al., 2019). BIRC5 was proposed as an upregulated gene and could be a significant marker for the detection and prognosis of BC at an early age (Ghaffari et al., 2016; Wang et al., 2018).

To investigate the pathogenetic processes of KGs, we selected four BP terms, six CC terms, three MF terms, and five KEGG terms based on directly associated with at least two KGs from the top ten terms of each category (Figure 3). Several studies also claimed that our identified KG-associated functional and pathway terms are responsible for BC progression (Emmert-Streib et al., 2014; Liu B. et al., 2018; Xiu et al., 2019; Han et al., 2020; Hermawan et al., 2020; Liu et al., 2020; Peng et al., 2021; Xiao et al., 2021; Zeng et al., 2021). GO functional analysis revealed that the BC-causing genes were significantly enriched in nuclear division (associated with KGs: TOP2A, AURKA, and BIRC5), organelle fission (associated with KGs: EGFR, AURKA, and BIRC5), and mitotic nuclear division (associated with KGs: MET and BIRC5) (Emmert-Streib et al., 2014; Han et al., 2020; Xiao et al., 2021). The three CC terms are condensed chromosome (associated with KGs: TOP2A, AURKA, and BIRC5), chromosomal region (associated with KGs: CDK1, AURKA, BIRC5, and EZH2) and midbody (associated with KGs: CDK1, AURKA, and BIRC5), which are associated with BC progression (Xiu et al., 2019; Liu et al., 2020). The MF term protein serine/threonine kinase activity (associated with KGs: CDK1, AURKA, and EGFR) is responsible for BC progression (Hermawan et al., 2020). Finally, the four KEGG pathway terms focal adhesion (associated with KGs: FN1, EGFR, and MET), progesterone-mediated oocyte maturation (associated with KGs: CDK1 and AURKA), oocyte meiosis (associated with KGs: CDK1 and AURKA), and adherens junction (associated with KGs: EGFR and MET), were also related to BC progression (Liu B. et al., 2018; Han et al., 2020; Peng et al., 2021; Zeng et al., 2021).

Four TFs (TP63, MYC, SOX2, and KDM5B) and four miRNAs (hsa-mir-16-5p, hsa-mir-34a-5p, hsa-mir-1-3p, and hsa-mir-23b-3p) were considered transcriptional and posttranscriptional and were connected to all KGs in the regulatory interaction network (Figure 4). Furthermore, we also checked the prognostic power of KGs for BC patients through the popular machine learning algorithm SVM. We developed a prediction model through SVM and showed good performance for both the training (GSE54002, GSE29431, GSE124646, GSE42568, GSE45827, GSE10810, GSE65216, and GSE36295) and test (GSE109169) (independent data) datasets. We found AUC values for eight KGs through the training dataset: EGFR (AUC: 96.1), FN1 (AUC: 96.8), EZH2 (AUC: 97.8), MET (AUC: 94.7), CDK1 (AUC: 98.6), AURKA (AUC: 99.7), TOP2A (AUC: 98.2), and BIRC5 (AUC: 98.7). Similarly, we found AUC values for eight KGs through the test dataset: EGFR (AUC: 93.4), FN1 (AUC: 91.5), EZH2 (AUC: 84.9), MET (AUC: 88.3), CDK1 (AUC: 93.5), AURKA (AUC: 92.7), TOP2A (AUC: 90.3), and BIRC5 (AUC: 85.5) (Figure 6). We noticed that the predictor model showed good performance with independent test datasets. These results indicate that our identified KGs have strong prognostic power for BC patients.

To discover potential candidate repurposing drugs for the treatment against BC patients, 77 KG-related drug agents were selected from the GSCALite database. Among them 16 candidate drugs (YM201636, masitinib, SB590885, GSK1070916, GSK2126458, ZSTK474, dasatinib, fedratinib, dabrafenib, methotrexate, trametinib, tubastatin A, BIX02189, CP466722, afatinib, and belinostat) were proposed based on BAS < −7 through molecular docking analysis (Figure 7; Table 5). Some of our proposed 16 drugs have already been approved, some are under clinical trial, and some are under in silico, in vivo, and in vitro analyses (Folch et al., 2015; Lee et al., 2015; Agency, 2017; FDA, 2017; European Commission, 2018; Odogwu et al., 2018; FDA, 2019; Hou et al., 2019; Drayman et al., 2021). Specifically, the drug YM201636 has been recommended which is effective against the progression of liver cancer through PIKfyve target inhibition (Hou et al., 2019). The drug masitinib was approved for mast cell disease and amyotrophic lateral sclerosis but later denied approval by the EU in 2017 and 2018 (Agency, 2017; European Commission, 2018). This drug has been studied for various human diseases and is now under clinical trial such as COVID-19 and Alzheimer’s disease (Folch et al., 2015; Drayman et al., 2021). Dasatinib was approved by the FDA in 2010 and 2017 for the treatment of adults with CP-CML and children with Ph + -CML, respectively (FDA, 2017). Fedratinib was approved by the FDA on 16 August 2019 for the treatment of MPN disease (FDA, 2019). Dabrafenib was approved by the FDA on 29 May 2013 for the treatment of BRAF V600E mutation-positive advanced melanoma disease (Odogwu et al., 2018). Tramatinib was approved by the FDA in May 2013 for the treatment of V600E mutated metastatic melanoma disease (Odogwu et al., 2018). Belinostat was approved by the FDA in 2013 for the treatment of peripheral T-cell lymphoma disease (Lee et al., 2015). However, the current study emphasizes experimental-lab verification for proposed candidate drugs with target proteins for treatment against BC.

Masitinib is an antitumor drug used primarily for mast cell tumors, but also used in solid tumors (Hahn et al., 2008). It inhibits the tyrosine kinase c-kit of mast cells, leading to a decrease in cell proliferation (Hahn et al., 2008). In BC, the cell surface receptors, including EGFR, HER-2, c-MET, and Trop2A, are cell surface markers that are potential targets for drugs (Butti et al., 2018). EGFR activation can induce the PI3K/AKT/mTOR signaling pathway to promote cell growth and proliferation (Miyamoto et al., 2017). Hyperactivation of the PI3K/AKT/mTOR pathway is commonly observed in BC, including triple-negative BC (Hussain et al., 2022). Inhibition of PI3K/AKT/mTOR by a dual PI3K/mTOR inhibitor apitolisib induces apoptotic cell death and decreases cell proliferation (Omeljaniuk et al., 2021). In the present study, we identified that masitinib as a candidate drug that potentially targets multiple BC-related products, including the cell surface marker EGFR. Using the breast cancer cell lines MCF-7 and MDA-MB-231, we observed that masitinib inhibits the activity of mTOR, a downstream effector of EGFR. It is well-documented that mTOR signaling is involved in cell proliferation, growth, metabolism, and survival (Zoungrana et al., 2022). In our observations, masitinib decreased cell migration and induced apoptotic cell death in MCF-7 and MDA-MB-231 cells. Thus, our data suggest that masitinib possibly acts on the EGFR/mTOR pathway to mediate antitumor activity.

Conclusion

In this study, we identified eight stable BC-guided KGs EGFR, FN1, EZH2, MET, CDK1, AURKA, TOP2A, and BIRC5 using well-established bioinformatics and network-based tools and highlighted their pathogenetic processes, regulatory factors, prognostic power, and drug molecules. Then, we proposed 16 potential candidate repurposing drugs YM201636, masitinib, SB590885, GSK1070916, GSK2126458, ZSTK474, dasatinib, fedratinib, dabrafenib, methotrexate, trametinib, tubastatin A, BIX02189, CP466722, afatinib, and belinostat through molecular docking analysis. Furthermore, we verified the proposed results through literature review and different databases. Finally, we examined one of the candidates masitinib to identify the anti-breast cancer effects. We found that masitinib inhibits the mTOR signaling pathway and induces apoptotic cell death. Therefore, the proposed molecular biomarkers and repurposing candidate drugs may play an important role in the diagnosis and therapy of BC.

References

• 1

  AgencyE. M. (2017). Refusal of the marketing authorisation for Masipro (masitinib). Available at: https://www.ema.europa.eu/en/documents/smop-initial/questions-answers-refusal-marketing-authorisation-masipro-masitinib_en.pdf.

  • Google Scholar

• 2

  AlamM. S.RahamanM. M.SultanaA.WangG.MollahM. N. H. (2022a). Statistics and network-based approaches to identify molecular mechanisms that drive the progression of breast cancer. Comput. Biol. Med.145, 105508. 10.1016/j.compbiomed.2022.105508

  • CrossRef
  • Google Scholar

• 3

  AlamM. S.SultanaA.RezaM. S.AmanullahM.KabirS. R.MollahM. N. H. (2022b). Integrated bioinformatics and statistical approaches to explore molecular biomarkers for breast cancer diagnosis, prognosis and therapies. PLoS One17, e0268967. 10.1371/journal.pone.0268967

  • CrossRef
  • Google Scholar

• 4

  BermanH. M.BattistuzT.BhatT. N.BluhmW. F.BourneP. E.BurkhardtK.et al (2002). The protein Data Bank. Acta Crystallogr. D. Biol. Crystallogr.58, 899–907. 10.1107/s0907444902003451

  • CrossRef
  • Google Scholar

• 5

  ButtiR.DasS.GunasekaranV. P.YadavA. S.KumarD.KunduG. C. (2018). Receptor tyrosine kinases (RTKs) in breast cancer: Signaling, therapeutic implications and challenges. Mol. Cancer17, 34. 10.1186/s12943-018-0797-x

  • CrossRef
  • Google Scholar

• 6

  CaoW.ChenH. D.YuY. W.LiN.ChenW. Q. (2021). Changing profiles of cancer burden worldwide and in China: A secondary analysis of the global cancer statistics 2020. Chin. Med. J.134, 783–791. 10.1097/CM9.0000000000001474

  • CrossRef
  • Google Scholar

• 7

  ClarkeC.MaddenS. F.DoolanP.AherneS. T.JoyceH.O'DriscollL.et al (2013). Correlating transcriptional networks to breast cancer survival: A large-scale coexpression analysis. Carcinogenesis34, 2300–2308. 10.1093/carcin/bgt208

  • CrossRef
  • Google Scholar

• 8

  DallakyanS.OlsonA. J. (2015). Small-molecule library screening by docking with PyRx. Methods Mol. Biol.1263, 243–250. 10.1007/978-1-4939-2269-7_19

  • CrossRef
  • Google Scholar

• 9

  DengJ. L.XuY. H.WangG. (2019). Identification of potential crucial genes and key pathways in breast cancer using bioinformatic analysis. Front. Genet.10, 695. 10.3389/fgene.2019.00695

  • CrossRef
  • Google Scholar

• 10

  DraymanN.DeMarcoJ. K.JonesK. A.AziziS. A.FroggattH. M.TanK.et al (2021). Masitinib is a broad coronavirus 3CL inhibitor that blocks replication of SARS-CoV-2. Science373, 931–936. 10.1126/science.abg5827

  • CrossRef
  • Google Scholar

• 11

  EdgarR.BarrettT. (2006). NCBI GEO standards and services for microarray data. Nat. Biotechnol.24, 1471–1472. 10.1038/nbt1206-1471

  • CrossRef
  • Google Scholar

• 12

  Emmert-StreibF.de Matos SimoesR.MullanP.Haibe-KainsB.DehmerM. (2014). The gene regulatory network for breast cancer: Integrated regulatory landscape of cancer hallmarks. Front. Genet.5, 15. 10.3389/fgene.2014.00015

  • CrossRef
  • Google Scholar

• 13

  European Commission (2018). Union Register of refused medicinal products for human use [Online]. Available at: https://ec.europa.eu/health/documents/community-register/html/ho26580.htm.

  • Google Scholar

• 14

  FDA (2017). FDA approves dasatinib for pediatric patients with CML [Online]. Available at: https://www.fda.gov/drugs/resources-information-approved-drugs/fda-approves-dasatinib-pediatric-patients-cml.

  • Google Scholar

• 15

  FDA (2019). FDA approves treatment for patients with rare bone marrow disorder [Online]. Available at: https://www.fda.gov/news-events/press-announcements/fda-approves-treatment-patients-rare-bone-marrow-disorder.

  • Google Scholar

• 16

  FolchJ.PetrovD.EttchetoM.PedrosI.AbadS.Beas-ZarateC.et al (2015). Masitinib for the treatment of mild to moderate Alzheimer's disease. Expert Rev. Neurother.15, 587–596. 10.1586/14737175.2015.1045419

  • CrossRef
  • Google Scholar

• 17

  GastaldiS.ComoglioP. M.TrusolinoL. (2010). The met oncogene and basal-like breast cancer: Another culprit to watch out for?Breast Cancer Res.12, 208. 10.1186/bcr2617

  • CrossRef
  • Google Scholar

• 18

  GhaffariK.HashemiM.EbrahimiE.ShirkoohiR. (2016). BIRC5 genomic copy number variation in early-onset breast cancer. Iran. Biomed. J.20, 241–245. 10.7508/ibj.2016.04.009

  • CrossRef
  • Google Scholar

• 19

  GruossoT.MieuletV.CardonM.BourachotB.KiefferY.DevunF.et al (2016). Chronic oxidative stress promotes H2AX protein degradation and enhances chemosensitivity in breast cancer patients. EMBO Mol. Med.8, 527–549. 10.15252/emmm.201505891

  • CrossRef
  • Google Scholar

• 20

  HahnK. A.OgilvieG.RuskT.DevauchelleP.LeblancA.LegendreA.et al (2008). Masitinib is safe and effective for the treatment of canine mast cell tumors. J. Vet. Intern. Med.22, 1301–1309. 10.1111/j.1939-1676.2008.0190.x

  • CrossRef
  • Google Scholar

• 21

  HanB.YangX.ZhangP.ZhangY.TuY.HeZ.et al (2020). DNA methylation biomarkers for nasopharyngeal carcinoma. PLoS One15, e0230524. 10.1371/journal.pone.0230524

  • CrossRef
  • Google Scholar

• 22

  HermawanA.PutriH.UtomoR. Y. (2020). Comprehensive bioinformatics study reveals targets and molecular mechanism of hesperetin in overcoming breast cancer chemoresistance. Mol. Divers.24, 933–947. 10.1007/s11030-019-10003-2

  • CrossRef
  • Google Scholar

• 23

  HouJ. Z.XiZ. Q.NiuJ.LiW.WangX.LiangC.et al (2019). Inhibition of PIKfyve using YM201636 suppresses the growth of liver cancer via the induction of autophagy. Oncol. Rep.41, 1971–1979. 10.3892/or.2018.6928

  • CrossRef
  • Google Scholar

• 24

  HussainY.KhanH.AlamW.AschnerM.AbdullahAlsharifK. F.et al (2022). Flavonoids targeting the mTOR signaling cascades in cancer: A potential crosstalk in anti-breast cancer therapy. Oxid. Med. Cell. Longev.2022, 4831833. 10.1155/2022/4831833

  • CrossRef
  • Google Scholar

• 25

  KarimS.MerdadA.SchultenH. J.JayapalM.DallolA.BuhmeidaA.et al (2016). Low expression of leptin and its association with breast cancer: A transcriptomic study. Oncol. Rep.36, 43–48. 10.3892/or.2016.4806

  • CrossRef
  • Google Scholar

• 26

  KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al (2019). PubChem 2019 update: Improved access to chemical data. Nucleic Acids Res.47, D1102–D1109. 10.1093/nar/gky1033

  • CrossRef
  • Google Scholar

• 27

  KoY. (2020). Computational drug repositioning: Current progress and challenges. Appl. Sci.10, 5076. 10.3390/app10155076

  • CrossRef
  • Google Scholar

• 28

  KomotoT. T.BernardesT. M.MesquitaT. B.BortolottoL. F. B.SilvaG.BitencourtT. A.et al (2018). Chalcones repressed the AURKA and MDR proteins involved in metastasis and multiple drug resistance in breast cancer cell lines. Molecules23, 2018. 10.3390/molecules23082018

  • CrossRef
  • Google Scholar

• 29

  LavecchiaA.Di GiovanniC. (2013). Virtual screening strategies in drug discovery: A critical review. Curr. Med. Chem.20, 2839–2860. 10.2174/09298673113209990001

  • CrossRef
  • Google Scholar

• 30

  LeeH. Z.KwitkowskiV. E.Del ValleP. L.RicciM. S.SaberH.HabtemariamB. A.et al (2015). FDA approval: Belinostat for the treatment of patients with relapsed or refractory peripheral T-cell lymphoma. Clin. Cancer Res.21, 2666–2670. 10.1158/1078-0432.CCR-14-3119

  • CrossRef
  • Google Scholar

• 31

  LiZ.HouP.FanD.DongM.MaM.LiH.et al (2017). The degradation of EZH2 mediated by lncRNA ANCR attenuated the invasion and metastasis of breast cancer. Cell Death Differ.24, 59–71. 10.1038/cdd.2016.95

  • CrossRef
  • Google Scholar

• 32

  LiuB.TanJ.WangX.LiuX. (2018a). Identification of recurrent risk-related genes and establishment of support vector machine prediction model for gastric cancer. Neoplasma65, 360–366. 10.4149/neo_2018_170507N326

  • CrossRef
  • Google Scholar

• 33

  LiuC. J.HuF. F.XiaM. X.HanL.ZhangQ.GuoA. Y. (2018b). GSCALite: A web server for gene set cancer analysis. Bioinformatics34, 3771–3772. 10.1093/bioinformatics/bty411

  • CrossRef
  • Google Scholar

• 34

  LiuD.ZhangJ.LiL.WangQ.LanY. (2020). Detection of critical genes associated with poor prognosis in breast cancer via integrated bioinformatics analyses. J. BUON25, 2537–2545.

  • Google Scholar

• 35

  MaubantS.TessonB.MaireV.YeM.RigaillG.GentienD.et al (2015). Transcriptome analysis of Wnt3a-treated triple-negative breast cancer cells. PLoS One10, e0122333. 10.1371/journal.pone.0122333

  • CrossRef
  • Google Scholar

• 36

  MiyamotoY.SuyamaK.BabaH. (2017). Recent advances in targeting the EGFR signaling pathway for the treatment of metastatic colorectal cancer. Int. J. Mol. Sci.18, 752. 10.3390/ijms18040752

  • CrossRef
  • Google Scholar

• 37

  NavolanicP. M.SteelmanL. S.McCubreyJ. A. (2003). EGFR family signaling and its association with breast cancer development and resistance to chemotherapy (Review). Int. J. Oncol.22, 237–252. 10.3892/ijo.22.2.237

  • CrossRef
  • Google Scholar

• 38

  OdogwuL.MathieuL.BlumenthalG.LarkinsE.GoldbergK. B.GriffinN.et al (2018). FDA approval summary: Dabrafenib and trametinib for the treatment of metastatic non-small cell lung cancers harboring BRAF V600E mutations. Oncologist23, 740–745. 10.1634/theoncologist.2017-0642

  • CrossRef
  • Google Scholar

• 39

  OmeljaniukW. J.KretowskiR.Ratajczak-WronaW.JablonskaE.Cechowska-PaskoM. (2021). Novel dual PI3K/mTOR inhibitor, apitolisib (GDC-0980), inhibits growth and induces apoptosis in human glioblastoma cells. Int. J. Mol. Sci.22, 11511. 10.3390/ijms222111511

  • CrossRef
  • Google Scholar

• 40

  PedrazaV.Gomez-CapillaJ. A.EscaramisG.GomezC.TorneP.RiveraJ. M.et al (2010). Gene expression signatures in breast cancer distinguish phenotype characteristics, histologic subtypes, and tumor invasiveness. Cancer116, 486–496. 10.1002/cncr.24805

  • CrossRef
  • Google Scholar

• 41

  PengJ.WuA.YuX.ZhongQ.DengX.ZhuY. (2021). Combined network Pharmacology and cytology experiments to identify potential anti-breast cancer targets and mechanisms of delphinidin. Nutr. Cancer74, 2591–2606. 10.1080/01635581.2021.2012582

  • CrossRef
  • Google Scholar

• 42

  RitchieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res.43, e47. 10.1093/nar/gkv007

  • CrossRef
  • Google Scholar

• 43

  SarveazadA.BabahajianA.ShamseddinJ.BahardoustM. (2018). 5-Year survival rates and prognostic factors in patients with synchronus and metachronus breast cancer from 2010 to 2015. Asian pac. J. Cancer Prev.19, 3489–3493. 10.31557/APJCP.2018.19.12.3489

  • CrossRef
  • Google Scholar

• 44

  ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res.13, 2498–2504. 10.1101/gr.1239303

  • CrossRef
  • Google Scholar

• 45

  ShiW.ChenX.DengL. (2020). A review of recent developments and progress in computational drug repositioning. Curr. Pharm. Des.26, 3059–3068. 10.2174/1381612826666200116145559

  • CrossRef
  • Google Scholar

• 46

  SimsekM.MeijerB.van BodegravenA. A.de BoerN. K. H.MulderC. J. J. (2018). Finding hidden treasures in old drugs: The challenges and importance of licensing generics. Drug Discov. Today23, 17–21. 10.1016/j.drudis.2017.08.008

  • CrossRef
  • Google Scholar

• 47

  SingT.SanderO.BeerenwinkelN.LengauerT. (2005). Rocr: Visualizing classifier performance in R. Bioinformatics21, 3940–3941. 10.1093/bioinformatics/bti623

  • CrossRef
  • Google Scholar

• 48

  SinnB. V.FuC.LauR.LittonJ.TsaiT. H.MurthyR.et al (2019). SETER/PR: A robust 18-gene predictor for sensitivity to endocrine therapy for metastatic breast cancer. NPJ Breast Cancer5, 16. 10.1038/s41523-019-0111-0

  • CrossRef
  • Google Scholar

• 49

  SongC. M.LimS. J.TongJ. C. (2009). Recent advances in computer-aided drug design. Brief. Bioinform.10, 579–591. 10.1093/bib/bbp023

  • CrossRef
  • Google Scholar

• 50

  SungH.FerlayJ.SiegelR. L.LaversanneM.SoerjomataramI.JemalA.et al (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca. Cancer J. Clin.71, 209–249. 10.3322/caac.21660

  • CrossRef
  • Google Scholar

• 51

  SzklarczykD.GableA. L.LyonD.JungeA.WyderS.Huerta-CepasJ.et al (2019). STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res.47, D607–D613. 10.1093/nar/gky1131

  • CrossRef
  • Google Scholar

• 52

  TanT. Z.MiowQ. H.MikiY.NodaT.MoriS.HuangR. Y.et al (2014). Epithelial-mesenchymal transition spectrum quantification and its efficacy in deciphering survival and drug responses of cancer patients. EMBO Mol. Med.6, 1279–1293. 10.15252/emmm.201404208

  • CrossRef
  • Google Scholar

• 53

  TangZ.LiC.KangB.GaoG.LiC.ZhangZ. (2017). GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res.45, W98-W102–W102. 10.1093/nar/gkx247

  • CrossRef
  • Google Scholar

• 54

  TorreL. A.BrayF.SiegelR. L.FerlayJ.Lortet-TieulentJ.JemalA. (2015). Global cancer statistics, 2012. Ca. Cancer J. Clin.65, 87–108. 10.3322/caac.21262

  • CrossRef
  • Google Scholar

• 55

  TrottO.OlsonA. J. (2010). AutoDock vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem.31, 455–461. 10.1002/jcc.21334

  • CrossRef
  • Google Scholar

• 56

  WangY.ZhangY.HuangQ.LiC. (2018). Integrated bioinformatics analysis reveals key candidate genes and pathways in breast cancer. Mol. Med. Rep.17, 8091–8100. 10.3892/mmr.2018.8895

  • CrossRef
  • Google Scholar

• 57

  WangZ.HuS.LiX.LiuZ.HanD.WangY.et al (2021). MiR-16-5p suppresses breast cancer proliferation by targeting ANLN. BMC Cancer21, 1188. 10.1186/s12885-021-08914-1

  • CrossRef
  • Google Scholar

• 58

  XiaoX.ZhangZ.LuoR.PengR.SunY.WangJ.et al (2021). Identification of potential oncogenes in triple-negative breast cancer based on bioinformatics analyses. Oncol. Lett.21, 363. 10.3892/ol.2021.12624

  • CrossRef
  • Google Scholar

• 59

  XiuY.LiuW.WangT.LiuY.HaM. (2019). Overexpression of ECT2 is a strong poor prognostic factor in ER(+) breast cancer. Mol. Clin. Oncol.10, 497–505. 10.3892/mco.2019.1832

  • CrossRef
  • Google Scholar

• 60

  YuG.WangL. G.HanY.HeQ. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS16, 284–287. 10.1089/omi.2011.0118

  • CrossRef
  • Google Scholar

• 61

  YuZ.HeQ.XuG. (2020). Screening of prognostic factors in early-onset breast cancer. Technol. Cancer Res. Treat.19, 1533033819893670. 10.1177/1533033819893670

  • CrossRef
  • Google Scholar

• 62

  ZengX.ShiG.HeQ.ZhuP. (2021). Screening and predicted value of potential biomarkers for breast cancer using bioinformatics analysis. Sci. Rep.11, 20799. 10.1038/s41598-021-00268-9

  • CrossRef
  • Google Scholar

• 63

  ZhangS.WangY.CaoY.WuJ.ZhangZ.RenH.et al (2022). Inhibition of the PINK1-parkin pathway enhances the lethality of sorafenib and regorafenib in hepatocellular carcinoma. Front. Pharmacol.13, 851832. 10.3389/fphar.2022.851832

  • CrossRef
  • Google Scholar

• 64

  ZhouG.SoufanO.EwaldJ.HancockR. E. W.BasuN.XiaJ. (2019). NetworkAnalyst 3.0: A visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res.47, W234–W241. 10.1093/nar/gkz240

  • CrossRef
  • Google Scholar

• 65

  ZoungranaL. I.Krause-HauchM.WangH.FatmiM. K.BatesL.LiZ.et al (2022). The interaction of mTOR and Nrf2 in neurogenesis and its implication in neurodegenerative diseases. Cells11, 2048. 10.3390/cells11132048

  • CrossRef
  • Google Scholar