Skip to main content

Deciphering the mechanism of Indirubin and its derivatives in the inhibition of Imatinib resistance using a “drug target prediction-gene microarray analysis-protein network construction” strategy



The introduction of imatinib revolutionized the treatment of chronic myeloid leukaemia (CML), substantially extending patient survival. However, imatinib resistance is currently a clinical problem for CML. It is very importantto find a strategy to inhibit imatinib resistance.


(1) We Identified indirubin and its derivatives and predicted its putative targets; (2) We downloaded data of the gene chip GSE2810 from the Gene Expression Omnibus (GEO) database and performed GEO2R analysis to obtain differentially expressed genes (DEGs); and (3) we constructed a P-P network of putative targets and DEGs to explore the mechanisms of action and to verify the results of molecular docking.


We Identified a total of 42 small-molecule compounds, of which 15 affected 11 putative targets, indicating the potential to inhibit imatinib resistance; the results of molecular docking verified these results. Six biomarkers of imatinib resistance were characterised by analysing DEGs.


The 15 small molecule compounds inhibited imatinib resistance through the cytokine-cytokine receptor signalling pathway, the JAK-stat pathway, and the NF-KB signalling pathway. Indirubin and its derivatives may be new drugsthat can combat imatinib resistance.

Peer Review reports


Chronic myeloid leukaemia (CML) is a clonal haematopoietic stem cell proliferation-induced myeloproliferative disease [1]. Because of its high heterogeneity and distinct molecular genetic features, it has attracted extensive attention from researchers. The unique cytogenetic features of CML include the Philadelphia chromosome t (9; 22) (q34; q11), forming a BCR-ABL fusion gene; this gene complex encodes a constitutively active form of the BCR–ABL fusion tyrosine kinase protein. The active site of the tyrosine kinase has a binding site for ATP [2]. Most signalling pathways activated by BCR-ABL are involved in promoting the development of cancer in bone marrow cells, including the Ras-MAPK pathway, the Src-Pax-Fak-Rac pathway, the phosphoinositide-3 kinase (PI3K)–Akt pathway, and the JAK-STAT pathway [3,4,5,6].

The development of the tyrosine kinase inhibitor (TKI) imatinib represents a milestone in CML treatment. Imatinib binds specifically to the ATP-binding site of BCR-ABL to form a fusion protein complex, locking in the active site [7]. This blocks CML cells whose active sites limit repeated cell growth and cell proliferation, killing the cancer cells. However, TKI treatment is long-term and induces resistance to TKI, often leading to poor clinical outcomes in CML patients. Drug resistance to TKIs is currently a clinical problem for CML. It isvery importantto find a strategy to inhibit imatinib resistance.

Classical traditional Chinese medicine (TCM) in China has been used for thousands of years. Especially in recent years, Chinese medicine has made some progress in the treatment of cancer. For example, Bu-Zhong-Yi-Qi-Decoction (BZYQD) has been reported to induce gastric cancer cell death by nonapoptotic mechanisms and to induce human ovarian cancer cell death by apoptotic mechanisms [8, 9]. Yu Ning, et al., through the combination BZYQD with cisplatin in cisplatin-resistant A549/DDP cells, showed that BZYQD exhibited direct cytotoxic and chemosensitising effects, suggesting that cotreatment with BZYQD and cisplatin might reverse cisplatin resistance by inducing ROS accumulation, activating apoptosis and autophagy by oxidative stress [10]. It was reported that Qingdai acted on a variety of pathways for the treatment of chronic myeloid leukaemia, including cytokine-cytokine receptor interaction, cell cycle, p53 signalling pathway, MAPK signalling pathway, and immune system-related pathways [11]. Indirubin is the most important and valuable compound in Qingdai; it has been determined to be the quality marker of Qingdai in the Chinese Pharmacopoeia (the State Pharmacopoeia Commission of China, 2015). Studies showed that indirubin and its derivatives inhibited imatinib resistance. For example, the AGM130 compound, derived from indirubin, known as a cyclin-dependent kinase inhibitor, was a strong candidate for treating imatinib-resistant CML [12]. Therefore, in this study, we will use the strategy of ‘Drug Target Prediction-Gene Microarray Analysis-Protein Network Construction’ to explore the mechanism of indirubin and its derivatives in inhibiting imatinib resistance.


To decipher the mechanisms by which indirubin and its derivatives reverse imatinib resistance, we adopted the following strategies: (1) we Identified the 2D structure of indirubin and its derivatives through data mining; (2) we downloaded GSE2810 from the GEO database and Identified imatinib-resistant DEGs; (3) we predicted targets of indirubin and its derivatives using related databases; (4) we analysed the possible molecular mechanisms of indirubin and its derivatives reversing imatinib resistance; and (5) we verified the results through computer network molecular docking technology.

Data preparation

Identify indirubin and its derivatives

We identified indirubin and its derivatives from two sources: first by searching the PubChem database and then by manually searching PubMed to augment the data. PubChem ( is a public repository for information on chemical substances and their biological activities. As of September 2015, it contained more than 157 million depositor-provided chemical substance descriptions, 60 million unique chemical structures and 1 million biological assay descriptions, covering approximately 10 thousand unique protein target sequences [13]. We searched the PubChem database with “indirubin” as the key word to identify indirubin and its derivatives, downloaded the compound 2D structures and finally downloaded the “smile” format. In order to increase the comprehensiveness of the data, we manually searched the relevant literature in the PubMed database for titles dealing with indirubin derivatives.

Identify the putative target of indirubin and its derivatives

It requires much manpower, material and financial resources to Identify targets of indirubin and its derivatives through experimentation. Therefore, we used a computerized virtual platform to screen for targets and then validated the targets by molecular docking or experimental verification. Swiss Target Prediction (, a web server to accurately predict the targets of bioactive molecules based on a combination of 2D and 3D similarity measures with known ligands, was used to predict the putative targets of the indirubin and its derivatives. Predictions can be carried out in five different organisms, and mapping predictions by homology within and between different species is enabled for close paralogs and orthologs [14]. The “smiles” formats of indirubin and its derivatives were imported into Swiss Target Prediction to predict their putative targets of action. It is noteworthy that the predicted putative target was limited to Homo sapiens. In order to improve the reliability of predictions goal, only high-probability targets were selected. All putative targets Identified were sent to the Therapeutic Target Database (TTD) (, 2015-09-10), the Comparative Toxicogenomics Database (CTD) (, 2017-12-05) and the PharmGKB ( to verify whether these putative targets had some connection to CML.

Identify imatinib resistance related genes

Gene expression profiling analysis is a useful method with broad clinical application in the identification of tumour-related genes in various types of cancer, from molecular diagnosis to pathological classification, from therapeutic evaluation to prognosis prediction, and from drug sensitivity to neoplasm recurrence [15]. Gene expression profile GSE2810 was downloaded from the Gene Expression Omnibus (GEO) database, GSE2810 data is based on the GPL2531 (Novusgene type 3 Hematology/Oncology TMU 667 array) platform,including 4 samples (2 imatinib-resistant samples and 2 imatinib-sensitive samples). It was submitted by Ohyashiki JH [16]. Quality control of gene expression data was performed using gene-specific probes. The analysis was carried out by using GEO2R, an online analysis tool for the GEO database, based on R language. We applied the analysis to classify the sample into two groups that had similar expression patterns in imatinib-sensitive and imatinib-resistant. We defined genes as differentially expressed (DEGs) when logFC > 1 or logFC < − 1(FC:Fold Change,the difference in the amount of gene expression in the sample). A p value < 0.05 was considered statistically significant. To further study the characteristics of DEGs and their functions, we analysed the DEGs with Gene Ontology and KEGG Pathway. Gene Ontology annotates and classifies genes by Molecular Function (MF), biological process (BP) and cellular component (CC). The pvalue of the GO term of the DEGs was calculated, and the most likely related GO term of the differential gene was located [17].KEGG is an online biochemical energy database that contains a set of genomic and enzymatic methods and is an information resource for the systematic analysis of gene functions and associated high-level genomic functions [18]. ClueGo, a plugin for Cytoscape 3.5.1 software, provides systematic and comprehensive biologically functional annotation of high-throughput gene expression [19]. Therefore, ClueGo online tools were employed for GO and KEGG pathway analysis. P < 0.05 was considered significant.

Network construction

Protein-protein network (P-P network). P-P network was built using the relationship between the putative targets of Indirubin and its derivatives and Imatinib resistance related DEGs.

Cytoscape 3.5.1 ( is an open software application for visualizing, integrating, modeling and analyzing interactive networks. All networks are built by it.

Analysis the protein-protein network

If the degree of a node is more than 2 fold of the median degree of all nodes in a network, such gene hub is believed to play a critical role in the network, and we treat it as major hub. The topological features of the target-target network are analysed by several important topological properties such as degree (the number of links to node) [20], betweenness (the number of shortest paths between pairs of nodes which run through node) [20], closeness(the sum of the distances of node to all other nodes) [20], and K-coreness (a measure of the centrality of node) [21]. The larger a protein’s degree/node betweenness/closeness centrality, the more important that protein is in the PPI network [22]. Subsequently, the targets were screened for topological importance. Then, the major hubs were screened. DAVIDwebserver ( was used to perform KEGG pathway enrichment analysis of the main targets.

Molecular docking simulation

Using computer molecular docking simulation techniques to verify the credibility of the study. SystemsDOCK ( were performed to Molecule docking [23]. SystemsDock, a web server for network pharmacology-based prediction and analysis, which permits docking simulation and molecular pathway map for comprehensive characterization of ligand selectivity and interpretation of ligand action on a complex molecular network, the score reported by docK-IN is a negative logarithm of the experimental dissociation/inhibition constant, usually ranging from 0 to 10 (i.e. from weak to strong binding). We conducted molecular docking between the small molecule compounds and their putative targets that are included in the major hubs selected by the P-P network map to evaluate whether indirubin and its derivatives inhibited imatinib resistance.


Data preparation

Indirubin and 41 derivatives and putative targets

We Identified indirubin and 41 derivatives from the database and downloaded “smiles” format and 2D structures. The putative targets of indirubin and its derivatives were predicted through structural similarities. Indirubin and 41 derivatives and putative targets are shown in Table 1.

Table 1 Indirubin and 41 derivatives and putative targets

Imatinib resistance related genes

After gene chip data analysis, we obtained a heat map of the differentially expressed genes of the gene chip G2810 (Additional file 1: Fig. S1), we Identified a total of 125 DEGs with imatinib resistance (Fig. 1), of which 66 were up-regulated and 59 were down-regulated. According to FC,the top 10 significantly up-regulated DEGs and down-regulated DEGs are shown in Table 2. Go analysis and KEGG analysis of DEGs, we found that DEGs of imatinib resistance were closely related to biological processes including immune responses, regulation of protein modification process, regulation of phosphorylation, and regulation of cellular protein metabolic processes. DEGs were mainly involved in cytokine-cytokine receptor interaction pathways.

Fig. 1

Based on GEO2R analysis, differentially expressed genes of imatinib resistance in chronic myeloid leukemia were Identified from GEO2810(logFC > 1 or logFC < − 1;P < 0.05), and a P-P network about DEGswas constructed. The red nodes represent up-regulated differentially expressed genes, and the blue nodes represent down-regulated differentially expressed genes

Table 2 The top 10 significantly up-regulated DEGs and down-regulated DEGs

CCL13, the first significantly up-regulated chemokine, is a chemotactic factor that attracts monocytes, lymphocytes, basophils and eosinophils [24]. MAPK11, the second significantly up-regulated chemokine, plays an important role in the cascades of cellular responses evoked by extracellular stimuli, including proinflammatory cytokines and physical stress leading to direct activation of transcription factors. The study of Huang J et al. showed that the ERK signalling pathway was more activated in epirubicin treated triple-negative breast cancer (TNBC), possibly contributing to epirubicin resistance, suggesting that the ERK pathway could be used as a novel candidate for targeting therapy in refractory and relapse TNBC [25]. MLH1, the first significantly down-regulated DEG, has been shown to play an important role in haematologic malignancies. The novel mutation was also revealed to be a somatic aberration occurring prior to the initiation of the blast phase in a chronic myelogenous leukaemia (CML) patient. Among the possible MLH1 partners involved in signalling MMR or apoptosis is the proto-oncogene c-MYC, closely associated with cellular proliferation [26]. BCL10, the second significantly down-regulated chemokine, was involved in adaptive immune responses. Proliferation of NIK and IKK cells is promoted by pro-caspase-9 maturation and NF-κB activation.

To further explain the function of differentially expressed genes, we performed functional enrichment analysis of all differential genes based on GO analysis, and performed passway enrichment analysis of all differential genes based on KEGG analysis. we chose significantly up-regulated and down-regulated GO categories based on functional enrichment, The analysis results are shown in Figs. 2 and 3. Through GO analysis, we reached the following conclusions: up-regulated differentially expressed genes were primarily involved in the regulation of cell apoptosis, including immune responses, regulation of apoptosis, regulation of programmed cell death, regulation of cell death, regulation of transcription, cell death, death and DNA binding. The down-regulated DEGs were primarily related to cellular structures, such as cytoplasm, nucleus, extracellular space, positive regulation of transcription from the RNA polymerase II promoter, transcription factor activity and sequence-specific DNA binding growth factor activity. We performed pathway enrichment analysis of differentially expressed genes to Identify the biological pathways. Up-regulated differentially expressed genes were primarily involved in cytokine-cytokine receptor interaction, chemokine signalling pathways, the Toll-like receptor signalling pathway, the neurotrophin signalling pathway, leukocyte transendothelial migration, the MAPK signalling pathway, haematopoietic cell lineage, apoptosis, the T cell receptor signalling pathway and the JAK-STAT signalling pathway. Pathways dramatically altered among down-regulated genes were the cytokine-cytokine receptor interaction, Toll-like receptor signalling pathway, Jak-STAT signalling pathway, pathways in cancer, the NOD-like receptor signalling pathway, apoptosis, cell cycle and the p53 signalling pathway.

Fig. 2

The significantly up-regulated GO categories and enrichment pathways of DEGs(P < 0.05)

Fig. 3

The significantly down-regulated GO categories and enrichment pathways of DEGs(P < 0.05)

To identify the relationship between the putative targets of indirubin and its derivatives and DEGs of imatinib resistance, we constructed a P-P network of putative targets and DEGs (Fig. 4). The T-T network consisted of 171 nodes and 1082 edges. The major hubs in the hub interaction network were determined by calculating four features: degree,betweenness,closeness and K-coreness. We showed the major hubs in Fig. 3. After screening, we identified a total of 62 major hubs (Table 3), including 11 (EGFR, JAK2, ERBB2, CHUK, CDK5, KIF11, DRD2, CDK3, HTR1A, JAK3 and TYK2) indirubin and derivative targets and 51 DEGs for imatinib resistance. These 11 major hubs were closely related to DEGs that were resistant to imatinib. Indirubin and its derivatives may inhibit imatinib resistance through the regulation of these genes.

Fig. 4

a P-P network, a co-expression network of the predicted target of indirubin and its derivatives and imatinib-resistant differentially expressed genes,the size of the node increases as the degree increases; b a network of 62 key nodes of the P-P network,the 11 nodes of yellow are not only the predicted targets of indirubin and its derivatives, but also the differentially expressed genes related to imatinib resistance

Table 3 The 62 major targets information of P-P network

We manually screened out small molecule compounds that affected 11 major hubs in the putative target. After screening, a total of 15 small molecule compounds affected these putative targets, including 1, 3, 4, 5, 6, 8, 11, 14, 21, 24, 26, 33,36, 40, 41. These derivatives may all inhibit imatinib resistance. To further verify this conclusion, we evaluated docking of small molecule compounds and their putative targets that were included in the major hubs. The docking results are shown in Table 4.

Table 4 The docking results of molecule compounds and their putative targets. ‘4 + EGFR’ represents the molecular docking of the indirubin derivative numbered 4 with EGFR, and Score represents the score Identified by molecular docking


Qingdai is a traditional Chinese medicine used to treat CML; it is the major active TCM of Qing-Huang-San [27], a Chinese traditional medicine used for the treatment of CML symptoms. It has been widely used in China and has achieved good clinical results. Indirubin is the major active component of Qingdai. Numerous studies have shown that indirubin and its derivatives not only promote apoptosis of CML cells but also inhibit imatinib resistance, including indirubin, indirubin derivative E804, and indirubin-3-acetoxime [28,29,30]. The exact mechanism of action remains unclear. Therefore, We used the Drug Target Prediction-Gene Microarray Analysis-Protein Network Construction model to investigate the mechanism by which indirubin and its derivatives inhibit imatinib resistance. Various methods, including indirubin derivative screening, drug target search screening, gene chip analysis, network construction, network target analysis, and molecular docking were combined to perform this study. A total of 42 small-molecule compounds were collected and predicted for putative targets. A total of 125 DEGs were selected for imatinib resistance. A total of 15 small-molecule compounds were found to inhibit imatinib resistance by 11 related genes. In our research, data mining of existing databases allows for the objective and rapid discovery of associations and identification of potential drug targets to facilitate the discovery of drugs that inhibit imatinib resistance.

CML is a major haematological malignancy. Imatinib is one of the primary drugs for the treatment of chronic myelogenous leukaemia; however, due to the resistance to imatinib, we were forced to study new drugs to inhibit the resistance to imatinib [31]. Drug resistance involves multiple steps and multiple genes. Therefore, various studies have analysed the differences in gene expression in imatinib-resistant and non-resistant genes by genomic microarrays. In the present study, we performed Go analysis and KEGG analysis on 125 differentially expressed genes and found that the resistance to imatinib was closely related to the following signalling pathways: (1) cell cycle, cell transcription, proliferation, apoptosis, and angiogenesis-related pathways; (2) cytokine-cytokine receptor interaction and chemokine signalling pathways; (3) cancer system related pathways, including pathways in cancer, the p53 signalling pathway and Jak-STAT signalling; (4) the immune system signalling pathway, the T cell receptor signalling pathway, the Toll-like receptor signalling pathway and the NOD-like receptor signalling pathway.

By analysing DEGs, we found that individual genes can serve as biomarkers for imatinib resistance. In up-regulated DEGs, CCL-13, the most significant up-regulated DEGs, is a chemokine that induces eosinophilic chemicals [32]; it can be involved in the interaction between haematopoietic stem cells and the bone marrow microenvironment [33]. In addition, the cytokine-cytokine receptor and chemokine signalling pathways involved in CCL-13 are important pathways involved in imatinib resistance. MAPK11 is the second most prominently expressed gene in the up-regulated differentially expressed genes for imatinib resistance, and MAPK11 is an important constituent gene of the MAPK signalling pathway and is involved in the regulation of various angiogenesis-related diseases [34]. The MAPK signalling pathway is significantlyaugmentedafter imatinib resistance and may be closely related to imatinib resistance. MAPK11 is also involved in up-regulating multiple regulatory pathways for DEGs, including the Toll-like receptor signalling pathway and leukocyte transendothelial migration. PIK3CD is involved in almost all pathways involved in the up-regulation of differentially expressed genes and is significantly augmentedin the course of imatinib resistance. Mesenchymal stem cells (MSC) from BM of chronic myeloid leukaemia (CML) patients on interaction with CML cells or its secreted factors, secreted high levels of IL6, providing a survival advantage to CML cells from imatinib-induced apoptosis [35]; Thus, IL6 may contribute to CML immune escape. Moreover, IL6 is involved in the cytokine-cytokine receptor interaction, the Jak-STAT signalling pathway, and pathways in cancer; therefore, it is closely related to imatinib resistance.

In the down-regulated DEGs, CASP8, an apoptosis-related factor, is an important apoptosis-related gene. Investigators used quantitative PCR to study apoptotic gene expression profile before and after imatinib treatment; they suggested that apoptosis-related gene expression profiles were associated with primary resistance to imatinib [36]. IL12A enhances cellular immunity in the treatment of CML. Studies have shown that immunotherapy enhanced the efficacy of imatinib, and low expression of IL12A led to immune escape of CML cells [37]. Therefore, CCL13, MAPK11, PIK3CD, IL6, CASP8, and IL12A play an important role in the process of imatinib resistance and can be used as biomarkers for imatinib resistance.

To elucidate the relationship between indirubin and its derivatives and imatinib resistance, we constructed a P-P network [38]. By analysing the P-P network, we found that there was a close relationship between the putative target of indirubin and its derivatives and DEGs of imatinib resistance. Through screening, we characterised a total of 11 putative targets [39]. Indirubin and its derivatives may inhibit imatinib resistance through these 11 putative targets. Based on 11 putative targets, we screened 15 small molecule compounds.

Among the 11 putative targets, gefitinib, an EGFR inhibitor, was tested in combination with imatinib in K562 CML cell line using MTT cell proliferation assay and was found to have a synergistic antiproliferative activity; EGFR inhibits or reverses imatinib resistance by enhancing the ability of imatinib to bind at the ATP-binding site of Bcr-Abl kinase [40]. The study found that JAK2 and JAK3 had antiproliferative effects on imatinib-resistant BCR-ABL(+) cells [41], and the administration of imatinib plus a JAK inhibitor reduced expression of stem cells markers, enhancing the antitumour effects of imatinib in CML cells [42]. Human ERBB2 is a proto-oncogene that codes for the erbB-2 epithelial growth factor receptor [43]. CHUK plays an important role in the NF-κB signalling pathway; indirubin and its derivatives inhibited CML cell proliferation by inhibiting CHUK activation of the NF-κB signalling pathway [44]. A study showed that NF-κB represents a potential target for molecular therapies in CML [45]. KIF11 inhibited cell proliferation by blocking the cycle of CML cells. The data showed that KIF11 was overexpressed in BCR-ABL+ CML cells and may become a novel treatment agent for patients with CML [46]. Administration of the imatinib plus JAK inhibitor reduces the expression of stem cell markers, such as ABCG2 and ALDH1A1. Blocking JAK3 with imatinib and JAK3 inhibitors may represent a new therapeutic strategy for eradicating LSCs and preventing CML recurrence [47].

We Identified a total of 15 small-molecule compounds that showed potential inhibition or reversal of resistance to imatinib. Active indirubins might inhibit T315I Abl kinase through unprecedented binding to both active and Src-like inactive conformations [30]. The AGM130 compound is derived from indirubin; data showed that the AGM130 compound efficiently decreased the viability of CML-derived K562 cells. Moreover, this compound also efficiently decreased the viability of imatinib-resistant CML cells in in vitro and in vivo systems [5]. E804, the most potent in indirubin derivative, blocked Stat5 signalling in human K562 CML cells, inhibiting the SFK/Stat5 signalling pathway downstream of Bcr-Abl, leading to apoptosis of K562, KCL-22 M and primary CML cells [48]. In the present study, we Identified small-molecule compounds of indirubin and its derivatives that could potentially inhibit imatinib resistance through drug target prediction, gene microarray analysis, and network construction, accelerating the discovery of new drugs for the treatment of imatinib resistance.

Finally, we used computer simulation techniques to dock selected small-molecule compounds to putative targets, and docking scores showed meaningful results, indicating that our series of strategies can achieve the desired results.


Definition of a potential drug target is an important first step in the process of drug discovery and drug design. Gene microarray analysis and protein network mapping can be key tools for identification of the factors that play a role in disease progression and thus are the potential drug targets. Subsequently, molecular docking experiments in silico can be used to predict putative interaction of small molecule compounds with the identified targets. In this study, based on the above methods, the mechanism of action of indirubin and its derivatives in inhibiting or reversing the resistance to imatinib was explored, and biomarkers and novel therapeutic targets that inhibited the resistance to imatinib were discovered. We validated experimental results by computerized molecular docking techniques. A limitation of this study was that the results were initially verified by computer simulation, and further verification can be achieved through experimental research.



Biological process




Cellular component


Chronic myeloid leukaemia


Comparative toxicogenomics database


The Database for Annotation, Visualization and Integrated Discovery


Differentially expressed genes


Gene expression omnibus


Gene ontology analysis


Kyoto encyclopedia of genes and genomes


Molecular function

P-P network:

Protein-protein network


Traditional Chinese medicine


Tyrosine kinase inhibitor


Therapeutic target database


  1. 1.

    Osorio S, Escudero-Vilaplana V, Gómez-Centurión I, González-Arias E, García-González X, Díez JL. Inadequate response to imatinib treatment in chronic myeloid leukemia due to a drug interaction with phenytoin. J Oncol Pharm Pract. 2017;1:1078155217743565.

    Google Scholar 

  2. 2.

    Kidan N, Khamaisie H, Ruimi N, Roitman S, Eshel E, Dally N. Ectopic Expression of Snail and Twist in Ph+ Leukemia Cells Upregulates CD44 Expression and Alters Their Differentiation Potential. J Cancer. 2017;8:3952–68.

    Article  Google Scholar 

  3. 3.

    Hilger RA, Scheulen ME, Strumberg D. The Ras-Raf-MEK-ERK pathway in the treatment of cancer. Onkologie. 2002;25:511–8.

    CAS  PubMed  Google Scholar 

  4. 4.

    Burridge K, Wennerberg K. Rho and Rac take center stage. Cell. 2004;116:167–79.

    CAS  Article  Google Scholar 

  5. 5.

    Bratton MR, Duong BN, Elliott S, Weldon CB, Beckman BS, McLachlan JA. Regulation of ERalpha-mediated transcription of Bcl-2 by PI3K–Akt crosstalk: implications for breast cancer cell survival. Int J Oncol. 2010;37:541–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Constantinescu SN, Girardot M, Pecquet C. Mining for JAK-STAT mutations in cancer. Trends Biochem Sci. 2008;33:122–31.

    CAS  Article  Google Scholar 

  7. 7.

    Eck MJ, Manley PW. The interplay of structural information and functional studies in kinase drug design, insights from BCR–ABL. Curr Opin Cell Biol. 2009;21:288–95.

    CAS  Article  Google Scholar 

  8. 8.

    Chao T, Fu P, Chang C, Chang S, Mao FC, Lin C. Prescription patterns of Chinese herbal products for post-surgery colon cancer patients in Taiwan. J Ethnopharmacol. 2014;155:702–8.

    Article  Google Scholar 

  9. 9.

    Zhu K, Fukasawa I, Furuno M, Inaba F, Yamazaki T, Kamemori T. Inhibitory effects of herbal drugs on the growth of human ovarian cancer cell lines through the induction of apoptosis. Gynecol Oncol. 2004;97:405–9.

    Article  Google Scholar 

  10. 10.

    Yu N, Xiong Y, Wang C. Bu-Zhong-Yi-Qi Decoction, the water extract of chinese traditional herbal medicine, enhances cisplatin cytotoxicity in A549/DDP cells through induction of apoptosis and autophagy. Biomed Res Int. 2017.

  11. 11.

    Li HY, Liu LJ, Liu C, et al. Deciphering key pharmacological pathways of Qingdai acting on chronic myeloid leukemia using a network pharmacology-based strategy. Med Sci Monit. 2018;24:5668–88.

    Article  Google Scholar 

  12. 12.

    Kim WS, Lee MJ, Kim DH, Lee JE, Kim JI, Kim YC. 5-OH-5-nitro-Indirubin oxime (AGM130), an Indirubin derivative, induces apoptosis of Imatinib-resistant chronic myeloid leukemia cells. Leuk Res. 2013;37:427–33.

    CAS  Article  Google Scholar 

  13. 13.

    Kim S, Thiessen PA, Bolton EE, Chen J, Fu G, Gindulyte A, Han L, He J, He S, Shoemaker BA, Wang J, Yu B, Zhang J, Bryant SH. PubChem substance and compound databases. Nucleic Acids Res. 2016;44:D1202-13.

  14. 14.

    David G, Aurelien G, Matthias W, Antoine D, Olivier M, Vincent Z. SwissTargetPrediction: a web server for target prediction of bioactive small molecules. Nucleic Acids Res. 2014;42:32–8.

    Google Scholar 

  15. 15.

    Freije WA, Castro-Vargas FE, Fang Z, Horvath S, Cloughesy T, Liau LM. Gene expression profiling of gliomas strongly predicts survival. Cancer Res. 2004;64:6503–10.

    CAS  Article  Google Scholar 

  16. 16.

    Nunoda K, Tauchi T, Takaku T, Okabe S, et al. Identification and functional signature of genes regulated by structurally different ABL kinase inhibitors. Oncogene. 2007;26:4179–88.

    CAS  Article  Google Scholar 

  17. 17.

    Gene Ontology Consortium. The gene ontology (GO) project in 2006. Nucleic Acids Res. 2006;34:D322–6.

    Article  Google Scholar 

  18. 18.

    Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    CAS  Article  Google Scholar 

  19. 19.

    Paul S, Andrew M, Owen O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    Article  Google Scholar 

  20. 20.

    Zhang YQ, Guo QY, Li QY, et al. Main active constituent identification in Guanxinjing capsule, a traditional Chinese medicine, for the treatment of coronary heart disease complicated with depression. Acta Pharmacol Sin. 2018;39:975–87.

    CAS  Article  Google Scholar 

  21. 21.

    Stefan W, Eivind A. Evolutionary cores of domain co-occurrence networks. BMC Evol Biol. 2005;5:24.

    Article  Google Scholar 

  22. 22.

    Wang Y, Liu Z, Li C. Drug target prediction based on the herbs components: the study on the multitargets pharmacological mechanism of Qishenkeli acting on the coronary heart disease. Evid Based Complement Alternat Med. 2012.

  23. 23.

    Hsin KY, Matsuoka Y, Asai Y, Kamiyoshi K, Watanabe T, Kawaoka Y. SystemsDock: a web server for network pharmacology-based prediction and analysis. Nucleic Acids Res. 2016;44:507–13.

    Article  Google Scholar 

  24. 24.

    Cossio-Ayala M, Domínguez-López M, Mendez-Enriquez E, Portillo-Téllez MDC, García-Hernández E. In vitro and in vivo antimicrobial activity of a synthetic peptide derived from the C-terminal region of human chemokine CCL13 against Pseudomonas aeruginosa. Peptides. 2017;94:49–55.

    CAS  Article  Google Scholar 

  25. 25.

    Huang J, Luo Q, Xiao Y, Li H, Kong L, Ren G. The implication from RAS/RAF/ERK signaling pathway increased activation in epirubicin treated triple negative breast cancer. Oncotarget. 2017;8:108249–60.

    PubMed  PubMed Central  Google Scholar 

  26. 26.

    Amikam D, Leshanski L, Sagi M. A novel MLH1 mutation harbored as a germ line aberration by a young woman of an HNPCC-like family and exhibited by a CML patient when occurring prior to the initiation of the blast phase concomitant with a c-MYC amplification. Int J Mol Med. 2006;17:1023–6.

    CAS  PubMed  Google Scholar 

  27. 27.

    Liu C, LIU LJ, ZHOU C, ZHUANG J, LIU H, SUN CG. Analysis of mechanism of indigo Naturalis in treating chronic Myelocytic leukemia based on three-dimentional model of protein-protein interaction network-molecular docking technique-in vitro experiment. Chin J Experiment Trad Med Form. 2017;23:206–11.

    Google Scholar 

  28. 28.

    Lee MY, Liu YW, Chen MH, Wu JY, Ho HY, Wang QF. Indirubin-3′-monoxime promotes autophagic and apoptotic death in JM1 human acute lymphoblastic leukemia cells and K562 human chronic myelogenous leukemia cells. Oncol Rep. 2013;29:2072–8.

    CAS  Article  Google Scholar 

  29. 29.

    Heshmati N, Wagner B, Cheng X, Scholz T, Kansy M, Eisenbrand G. Physicochemical characterization and in vitro permeation of an indirubin derivative. Eur J Pharm Sci. 2013;50:467–75.

    CAS  Article  Google Scholar 

  30. 30.

    Gaboriaud-Kolar N, Myrianthopoulos V, Vougogiannopoulou K, Gerolymatos P, Horne DA, Jove R. Natural-based Indirubins display potent cytotoxicity toward wild-type and T315I-resistant leukemia cell lines. J Nat Prod. 2016;79:2464–71.

  31. 31.

    Larocque EA, Naganna N, Opoku-Temeng C, Lambrecht AM, Sintim HO. Alkynylnicotinamide-Based Compounds as ABL1 Inhibitors with Potent Activities against Drug-Resistant CML Harboring ABL1(T315I) Mutant Kinase. ChemMedChem. 2018; p. 1172-80.

  32. 32.

    Yang T, Li Y, Lyu Z, Huang K, Corrigan CJ, Ying S. Characteristics of Proinflammatory Cytokines and Chemokines in Airways of Asthmatics: Relationships with Disease Severity and Infiltration of Inflammatory Cells. Chin Med J (Engl). 2017;130:2033–40.

    Article  Google Scholar 

  33. 33.

    Mukaida N, Tanabe Y, Baba T. Chemokines as a conductor of bone marrow microenvironment in chronic myeloid leukemia. Int J Mol Sci. 2017;18:1824.

    Article  Google Scholar 

  34. 34.

    Chun YX, Hui ZD, Bond LW, Qiang LK, Yu TJ, Chao CZ, Liang MX, Hua LJ, Fan Q. lncRNA ENSMUST00000134285 Increases MAPK11 Activity, Regulating Aging-Related Myocardial Apoptosis. J Gerontol A Biol Sci Med Sci. 2018;73:1010-17.

  35. 35.

    Kumar A, Anand T, Bhattacharyya J, Sharma A, Jaganathan BG. K562 chronic myeloid leukemia cells modify osteogenic differentiation and gene expression of bone marrow stromal cells. J Cell Commun Signal. 2018;12:441-50.

  36. 36.

    Ferreira AF, Oliveira GL, Tognon R, Collassanti MD, Zanichelli MA, Hamerschlak N. Apoptosis-related gene expression profile in chronic myeloid leukemia patients after imatinib mesylate and dasatinib therapy. Acta Haematol. 2015;133:354–64.

    CAS  Article  Google Scholar 

  37. 37.

    Tao K, Li YJ, Wang D, Qi JY, Deng YP, Wang HX. Enhancement of specific cellular immune response induced by glycosyl-phosphatidylinositol-anchored BCR/ABL and mIL-12. Cancer Biol Ther. 2011;12:881–7.

    CAS  Article  Google Scholar 

  38. 38.

    Yu GH, ZhangYQ RWQ, Dong L, Li JF, Geng Y. Network pharmacology-based identification okey pharmacological pathways of Yin–Huang–Qing–Fei capsule acting on chronic bronchitis. Int J Chron Obstruct Pulmon Dis. 2016;12:85–94.

    Article  Google Scholar 

  39. 39.

    Yue SJ, Liu J, Feng WW, Zhang FL, Chen JX, Xin LT. System pharmacology-BaseDissection of the synergistic mechanism of Huangqi and Huanglian for diabetes mellitus. Front Pharmacol. 2017;8:694.

    Article  Google Scholar 

  40. 40.

    Singh VK, Chang HH, Kuo CC, Shiao HY, Hsieh HP, Coumar MS. Drug repurposing for chronic myeloid leukemia: in silico and in vitro investigation of DrugBank database for allosteric Bcr-Abl inhibitors. J Biomol Struct Dyn. 2017;35:1833–48.

    CAS  Article  Google Scholar 

  41. 41.

    Yagi K, Shimada A, Sendo T. Pharmacological inhibition of JAK3 enhances the antitumor activity of imatinib in human chronic myeloid leukemia. Eur J Pharmacol. 2018;825:28–33.

    CAS  Article  Google Scholar 

  42. 42.

    Tanaka R, Squires MS, Kimura S, Yokota A, Nagao R, Yamauchi T. Activity of the multitargeted kinase inhibitor, AT9283, in imatinib-resistant BCR-ABL-positive leukemic cells. Blood. 2010;116:2089–95.

    CAS  Article  Google Scholar 

  43. 43.

    Santos S, Baptista CS, Abreu RM, Bastos E, Amorim I, Gut IG. ERBB2 in cat mammary neoplasias disclosed a positive correlation between RNA and protein low expression levels: a model for erbB-2 negative human breast cancer. PLoS One. 2013;8:e83673.

    Article  Google Scholar 

  44. 44.

    Kirchner D, Duyster J, Ottmann O, Schmid RM, Bergmann L, Munzert G. Mechanisms of Bcr-Abl-mediated NF-kappaB/Rel activation. Exp Hematol. 2003;31:504–11.

    CAS  Article  Google Scholar 

  45. 45.

    Lu Y, Jamieson L, Brasier AR, Fields AP. NF-kappaB/RelA transactivation is required for atypical protein kinase C iota-mediated cell survival. Oncogene. 2001;20:4777–92.

    CAS  Article  Google Scholar 

  46. 46.

    Yin Y, Sun H, Xu J. Kinesin spindle protein inhibitor SB743921 induces mitotic arrest and apoptosis and overcomes imatinib resistance of chronic myeloid leukemia cells. Leuk Lymphoma. 2015;56:1813–20.

    CAS  Article  Google Scholar 

  47. 47.

    Kenta Y, Akira S, Toshiaki S. Pharmacological inhibition of JAK3 enhances the antitumor activity of imatinib in human chronic myeloid leukemia. Eur J Pharmacol. 2018;825:28–33.

    Article  Google Scholar 

  48. 48.

    Nam S, Scuto A, Yang F, Chen W, Park S, Yoo HS. Indirubin derivatives induce apoptosis of chronic myelogenous leukemia cells involving inhibition of Stat5 signaling. Mol Oncol. 2012;6:276–83.

    CAS  Article  Google Scholar 

Download references


Not applicable.


This work is supported by the grants from National Natural Science Foundation of China (No.81673799) and National Natural Science Foundation of China Youth Fund (No.81703915).

Availability of data and materials

All data generated or analysed during this study are included in this published article.

Author information




SCG and LHY conceived and designed the study; LLJ, LC and ZJ performed the study; LLJ, LHY and ZC analyzed the data; YJ, GCD, LGX, and LQL contributed analysis tools; LHY and LLJ wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Changgang Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional information

Huayao Li and Lijuan Liu are co-first author.

Additional file

Additional file 1:

Figure S1. Heat maps of differentially expressed genes associated with imatinib resistance (we selected 100 genes with the most significant differential expression) (P < 0.05). The color from blue to red shows a trend from low to high expression. (JPG 298 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, H., Liu, L., Zhuang, J. et al. Deciphering the mechanism of Indirubin and its derivatives in the inhibition of Imatinib resistance using a “drug target prediction-gene microarray analysis-protein network construction” strategy. BMC Complement Altern Med 19, 75 (2019).

Download citation


  • Indirubin
  • Derivatives
  • Imatinib resistance
  • Drug target prediction
  • Gene microarray analysis
  • Protein network construction