Skip to main content

Insights into the mechanism of Arnebia euchroma on leukemia via network pharmacology approach

Abstract

Background

Arnebia euchroma (A. euchroma) is a traditional Chinese medicine (TCM) used for the treatment of blood diseases including leukemia. In recent years, many studies have been conducted on the anti-tumor effect of shikonin and its derivatives, the major active components of A. euchroma. However, the underlying mechanism of action (MoA) for all the components of A. euchroma on leukemia has not been explored systematically.

Methods

In this study, we analyzed the MoA of A. euchroma on leukemia via network pharmacology approach. Firstly, the chemical components and their concentrations in A. euchroma as well as leukemia-related targets were collected. Next, we predicted compound-target interactions (CTIs) with our balanced substructure-drug-target network-based inference (bSDTNBI) method. The known and predicted targets of A. euchroma and leukemia-related targets were merged together to construct A. euchroma-leukemia protein-protein interactions (PPIs) network. Then, weighted compound-target bipartite network was constructed according to combination of eight central attributes with concentration information through Cytoscape. Additionally, molecular docking simulation was performed to calculate whether the components and predicted targets have interactions or not.

Results

A total of 65 components of A. euchroma were obtained and 27 of them with concentration information, which were involved in 157 targets and 779 compound-target interactions (CTIs). Following the calculation of eight central attributes of targets in A. euchroma-leukemia PPI network, 37 targets with all central attributes greater than the median values were selected to construct the weighted compound-target bipartite network and do the KEGG pathway analysis. We found that A. euchroma candidate targets were significantly associated with several apoptosis and inflammation-related biological pathways, such as MAPK signaling, PI3K-Akt signaling, IL-17 signaling, and T cell receptor signaling pathways. Moreover, molecular docking simulation demonstrated that there were eight pairs of predicted CTIs had the strong binding free energy.

Conclusions

This study deciphered that the efficacy of A. euchroma in the treatment of leukemia might be attributed to 10 targets and 14 components, which were associated with inhibiting leukemia cell survival and inducing apoptosis, relieving inflammatory environment and inhibiting angiogenesis.

Peer Review reports

Background

Leukemia is a common malignancy in children and adults originating from pluripotent hematopoietic stem cells. Due to various reasons, leukemia cells undergo enhanced self-renewal, over proliferation, and blocked differentiation and apoptosis. Abnormal immature cells proliferate and accumulate in bone marrow and other hematopoietic tissues, leading to the suppression of normal hematopoiesis and infiltration of other organs and tissues. According to the types of cells affected and the developmental stage of the originating cells, leukemia is classified into four major categories: acute myeloid leukemia (AML), acute lymphoblastic leukemia (ALL), chronic myeloid leukemia (CML), and chronic lymphoblastic leukemia (CLL).

The usual approach to treat leukemia is the use of toxic compounds to kill cancer cells, which will eventually destroy the immune system, and make patients susceptible to fatal bacterial and fungal infections. For that reason, it was suggested that approaches such as inducing leukemia cell differentiation and apoptosis other than killing cells may be more effective in the treatment of leukemia [1]. Traditional Chinese medicine (TCM) has shown potentials for the treatment of leukemia, for example, arsenic trioxide (ATO), the main component of TCM white arsenic, has been applied in the treatment of leukemia [2].

Arnebia euchroma (Royle) Jonst., also called Xinjiang Zicao, is a TCM used in the treatment of blood diseases including leukemia [3]. It has been reported that among the 295 empirical prescriptions for treating haematological diseases in Chinese medicine, the frequency of A. euchroma is 10.3%, and the frequency of A. euchroma powder, a prescription that takes A. euchroma as the monarch herb is 53.1% [4]. The components of A. euchroma consist of two major categories: one is hydrophilic components, mainly a mixture of polysaccharides and glycoproteins; the other is lipophilic substances with a variety of biological activities, including naphthoquinones, alkaloids, monoterpene phenol and benzoquinones, organic acid ester, and so on. Among these, typical naphthoquinones include shikonin and a series of derivatives [5], which have a common scaffold of 5,8-dihydroxy-2-isohexene-1,4-naphthoquinone. According to the optical configuration of the chiral center in the side chain, naphthoquinones can be divided into two subtypes: R-configuration (named shikonin) and S-configuration (named alkannin) (see Supplementary Figure S1) [6]. Studies have shown that naphthoquinones are the major active ingredients of A. euchroma that exert its pharmacological effects, whose concentrations are no more than 7% [7, 8]. Especially, some studies demonstrated that shikonin has therapeutic effects on leukemia cells, mainly including HL-60 cells and K562 cells [9,10,11]. However, the underlying mechanism of action (MoA) for all the components of A. euchroma on leukemia is barely explored systematically.

In 2007, network pharmacology was proposed to describe the multiple interactions among drugs, targets and diseases [12, 13], which coincides with the “multi-component, multi-target” characteristics of TCM and has been widely used in the research and development of TCM. In this field, our group has developed a series of computational methods to predict compound-target interactions (CTIs), including network-based inference (NBI), substructure-drug-target network-based inference (SDTNBI), and balanced SDTNBI (bSDTNBI) [14,15,16]. In previous studies we have demonstrated that these methods could predict potential targets for TCM components and help to understand MoA of TCM reasonably [17,18,19].

Since NBI only can predict targets for old compounds within a known CTI network, while bSDTNBI has demonstrated better performance than SDTNBI in prediction of targets for new compounds outside the known CTI network in our previous publications. Therefore, in this study, we used bSDTNBI to predict targets for components of A. euchroma, and then investigated the material basis and MoA of A. euchroma on leukemia via network pharmacology approach, which may provide some basis and enlightenment to deeply explore the chemical and pharmacological basis of TCM.

Methods

Data collection and preparation

The components and their concentrations in A. euchroma were collected from several sources including TCM-MESH (http://mesh.tcm.microbioinformatics.org/, entering at Apr 2019) [20], TCMID (http://119.3.41.228/tcmid/, entering at Apr 2019) [21], Chinese Natural Product Chemical Composition Library (http://pharmdata.ncmi.cn/cnpc/, entering at Apr 2019), and Phytochemicals database (http://chemdb.sgst.cn/scdb/main/plant_introduce.asp, entering at Apr 2019), as well as related literature from PubMed (https://www.ncbi.nlm.nih.gov/pubmed/, entering at Apr 2019), Web of Science (http://apps.webofknowledge.com, entering at Apr 2019) and CNKI (http://www.cnki.net/, entering at Apr 2019). The structures of these compounds were obtained from NCBI PubChem Compound Database (https://www.ncbi.nlm.nih.gov/pccompound, entering at Apr 2019). Duplicates were removed.

Known compound-target interactions (KCTIs) for each component of A. euchroma were collected from four databases, including BindingDB (http://www.bindingdb.org/bind/index.jsp, entering at Apr 2019) [22], IUPHAR/BPS Guide to PHARMACOLOGY (https://www.guidetopharmacology.org/, entering at Apr 2019) [23], ChEMBL (https://www.ebi.ac.uk/chembl/, entering at Apr 2019) [24], and PubChem (https://pubchem.ncbi.nlm.nih.gov/, entering at Apr 2019) [25]. Target protein names of components were converted into gene names with the species limited into “Homo sapiens” by UniProt (https://www.uniprot.org/, entering at Apr 2019). Duplicated KCTIs were removed.

Leukemia related targets were obtained from five databases: OMIM (https://www.omim.org/, updated on Jan. 3, 2018) [26], PharmGkb (https://www.pharmgkb.org/, updated on Dec. 28, 2017) [27], TTD (http://bidd.nus.edu.sg/group/cjttd/, updated on Sep. 15, 2017) [28], DisGeNET (https://www.disgenet.org/search, entering at Apr 2019) [29], and MalaCards (https://www.malacards.org/, entering at Apr 2019) [30].

Prediction of compound-target interactions

The bSDTNBI method was used to predict potential compound-target interactions (PCTIs) for each component of A. euchroma. In bSDTNBI method, three parameters, namely α, β and γ, should be determined first. Among them, α controls the relative importance of substructure vs. target nodes, and β defines the importance of drug-substructure vs. drug-target edges, while γ determines the importance of hub nodes. The three parameters were optimized by grid search and 10-fold cross validation. First of all, α and β were optimized under the condition of ignoring the influence of hub nodes (γ = 0). Then, under the optimal values of α and β, the optimal value of parameter γ was searched. Therefore, with the optimal parameters α, β and γ, bSDTNBI method was used to predict 10 potential targets for each component of A. euchroma. The predicted targets were further standardized as official gene symbols. Eventually, the collected known and predicted CTIs were integrated.

Construction of PPI network and selection of leukemia-related targets

Targets of components of A. euchroma (TA) and targets related to leukemia (TL) were imported into Cytoscape 3.6.0, respectively. Protein-protein interactions (PPIs) of TA and TL were obtained by stringAPP, a Cytoscape plugin. High-confidence PPIs with scores above 0.7 were selected to construct A. euchroma PPI network (PPIA) and leukemia PPI network (PPIL). Then PPIA and PPIL were merged together to construct A. euchroma-leukemia PPI network.

A. euchroma-leukemia PPI network was analyzed by CytoNCA, another Cytoscape plugin. Each node in the network has eight central attributes, including subgraph centrality (SC), degree centrality (DC), eigenvector centrality (EC), information centrality (IC), local average connectivity-based method (LAC), betweenness centrality (BC), closeness centrality (CC), and network centrality (NC). For all target nodes, the targets of TA that meet the screening criteria “SC > median SC & DC > median DC & EC > median EC & IC > median IC & LAC > median LAC & BC > median BC & CC > median CC & NC > median NC” were regarded as leukemia-related targets of A. euchroma (TAL) and were retained [31].

Construction of weighted compound-target bipartite network

Compound-target interactions for TAL were obtained from all CTIs (KCTIs + PCTIs). Here the concentrations of components in A. euchroma were used to weight the importance of each component. Since the concentrations of components were determined under different experimental conditions and influenced by many factors, such as the origin and growth time of A. euchroma, we processed the original concentration data as follows. Firstly, all the concentrations of components were summarized and standardized to a uniform percentage (mass percentage). Then, the average value of percentage of each component was calculated after the maximum and minimum values were removed. Finally, with “percent concentration value = 0.1%” as the threshold, components were divided into three classes, and the formula was as following:

$$ Components\left\{\begin{array}{c} Class\ I: weight=3,\kern0.5em concentration>0.1\%\\ {} Class\ II:\kern0.5em weight=2,\kern0.5em 0< concentration<0.1\%\\ {} Class\ II I: weight=1,\kern0.75em without\ concentration\ data\end{array}\right\} $$

After all components were weighted, a weighted compound-target bipartite network was constructed through Cytoscape.

Enrichment of KEGG pathways

KEGG pathway enrichment analysis is usually used to describe the characteristics of query targets. Here, STRING 10.5 (https://string-db.org/cgi/input.pl, entering at Apr 2019) [32] was used to perform KEGG pathway enrichment analysis. The KEGG pathways with FDR (False Discovery Rate) < 0.05 were regarded as significant and useful. Targets in weighted compound-target bipartite network were inputted to STRING, and KEGG pathways with FDR < 0.05 were obtained. A scoring function was designed to screen pathways, and the equation was as following:

$$ Norm(Y)=\left(1-0.1\right)\ast \frac{Y_i-\mathit{\operatorname{MIN}}(Y)}{\mathit{\operatorname{MAX}}(Y)-\mathit{\operatorname{MIN}}(Y)}+0.1 $$
$$ {PS}_m= Norm\left(\sum \limits_{j=1}^{j={N}_T} Norm\left(\sum \limits_{i=1}^{i={N}_C}{X}_i\right)\right)+ Norm\left({LFDR}_m\right) $$

in which, PSm refers to the score of the m-th KEGG pathway; NT refers to the number of targets involved in the m-th pathway; NC refers to the number of components that interact with the j-th target in the compound-target bipartite network; Xi refers to the weight of the i-th component; LFDRm refers to the negative logarithmic value of FDR of the m-th pathway; Norm refers to the min-max normalization method with normalize data from 0.1 to 1.

After all pathways were scored, pathways that meet the role “score > 1.5 median score” were selected. Referring to the pathway classification standard of KEGG Pathway Database, here, these selected pathways were divided into four categories, including cell survival and death, immune system, endocrine system and specific human diseases. In order to facilitate analysis, we removed specific human disease pathways, which were not directly related to the treatment mechanism of leukemia, and the other pathways were remained.

Selection of hub genes and molecular docking simulation

Targets that were involved in the above qualified pathways were selected out and scored further by the following scoring function:

$$ Norm(X)=\left(1-0.1\right)\ast \frac{X_i-\mathit{\operatorname{MIN}}(X)}{\mathit{\operatorname{MAX}}(X)-\mathit{\operatorname{MIN}}(X)}+0.1 $$
$$ {TS}_{\boldsymbol{j}}={N}_p\ast Norm\left(\sum \limits_{i=1}^{\boldsymbol{i}={N}_{\boldsymbol{c}}}{X}_i\right) $$

in which, TSj refers to the score of the j-th target; Np refers to the number of pathways that the j-th target involved in; Nc refers to the number of components that have interactions with the j-th target; Xi refers to the weight of the i-th component; Norm refers to the min-max normalization method with normalize data from 0.1 to 1.

After all targets were scored, targets that meet “score > 1.5 median score” were regarded as hub genes, and CTIs that were associated with these targets were selected out. The molecular docking simulation was further performed to detect whether the components and predicted targets form reasonable interactions and estimate their binding affinities.

First, we collected the crystal structures of the selected protein from the RCSB Protein Data Bank (PDB, http://www.pdb.org/, updated on 2020-4-10) and selected the relatively higher resolution crystal structures with the ligands. Second, the structures of chemical components contained in A. euchroma were downloaded from the NCBI PubChem Compound Database. Then, protein preparation module of Schrödinger’s Maestro Molecular modeling suite (Schrödinger Release 2015–2) was utilized for preparation of the protein crystallographic structures. Water molecules were subsequently deleted from the structures, and the amide moieties in the side chain were adjusted to optimize their interactions with surrounding residues and groups of atoms. Force field OPLS_2005 was also added. A ligand grid generation was based on the ligand in the co-crystallographic structure. The docking region was centered on the ligand, and after the ligand grid was generated, the compounds were imported into Mastero. LigPrep module of the Maestro molecular modeling package was used to obtain the 3D structures and energy minimization of the identified compounds. Compounds were docked in the generated grid using the standard Glide docking mode.

Results

In this study, we took the concentrations of components in A. euchroma into account, and analyzed the MoA of A. euchroma on leukemia by network pharmacology approach, which involved three steps in the whole workflow (Fig. 1): (1) data collection. The components of A. euchroma along with their concentrations, corresponding targets and leukemia related targets were identified by various databases and literature; the predicted targets were obtained through bSDTNBI method; (2) network establishment. A. euchroma-leukemia PPI network and weighted compound-target bipartite network were constructed through Cytoscape; (3) network analysis. Hub genes were selected out according to KEGG pathways analysis, and molecular docking was used to calculate whether the components from A. euchroma and predicted targets have interactions or not and to construct the critical compound-target bipartite network. Figure 1 depicts a workflow of the technical strategy used in this study.

Fig. 1
figure1

Workflow for A. euchroma on leukemia

Components and concentrations of A. euchroma and leukemia related targets

During the process of data collection, all the components were collected from several databases and publications based on two criteria. On the one hand, the component should be found in at least two sources, such as two databases, or a database and a publication, or two publications; on the other hand, the component should have a clear structure. Finally, a total of 65 components of A. euchroma were obtained, and 27 of them with concentration information (see Supplementary Table S1). Considering that our bSDTNBI method could not distinguish stereoisomers, components that were mutually stereoisomers were treated as the same molecule among the 65 components.

According to the collected concentration information of each component, after the maximum and minimum values were removed, the average of percentage (mass percentage) was calculated. Finally, 11 components with concentration greater than 0.1% had a weight of 3 (Concentration Score, X = 3), 16 components with concentration less than 0.1% but greater than 0 had a weight of 2 (X = 2), and the remaining 38 components without collected concentration data had a weight of 1 (X = 1).

From OMIM, PharmGkb, TTD, DisGeNET, and MalaCards, a total of 332 targets were ultimately reserved after removing the duplicates.

Known and predicted compound-target interactions

A total of 129 KCTIs for 13 components (MOL16, MOL18, MOL36, MOL38, MOL42, MOL43, MOL46, MOL47, MOL49, MOL53, MOL55, MOL56, and MOL57) were collected from databases, which were involved in 92 targets (see Supplementary Table S2). Only 13 of the 65 components of A. euchroma have known targets. Therefore, it is especially necessary to predict potential targets for components.

Herein, our bSDTNBI method was used to predict 10 targets for each of the 65 components. After grid search and 10-fold cross validation, the optimal values of the three parameters were set as α = 0.3, β = 0.1, γ = − 0.4. Then a total of 79 potential targets were obtained, resulting in 650 new CTIs. After the known and predicted CTIs were merged, duplicated CTIs were then removed. Finally, the number of retained CTIs was 779, involved in 65 components and 157 targets (see Supplementary Table S2).

A. euchroma-leukemia PPI network

157 targets (TA = 157) of 65 components in A. euchroma and 332 leukemia related targets (TL = 332) were imported to Cytoscape, respectively. When the confidence (score) cutoff was set to 0.7, 129 targets out of TA had 364 PPIs (PPIA = 364), while 270 targets out of TL had 1335 PPIs (PPIL = 1335). When they were merged, 1658 PPIs were obtained. Among these, there were 709 PPIs for the 129 targets of A. euchroma components. These 709 PPIs were selected out and re-submitted to Cytoscape to construct A. euchroma-leukemia PPI network. In the PPI network, eight central attributes were calculated by CytoNCA (see Supplementary Table S3). According to the screening criteria, there were 37 targets of A. euchroma components (TAL) were qualified (Fig. 2), which corresponding to 55 components and 114 CTIs.

Fig. 2
figure2

Construction of A. euchroma-leukemia PPI network. a A. euchroma PPI network (PPIA); b Leukemia PPI network (PPIL); c A. euchroma-leukemia PPI network; d Targets of A. euchroma that meet the screening criteria

Weighted compound-target bipartite network

114 CTIs of the 37 TAL were imported into Cytoscape to construct the weighted compound-target bipartite network, as shown in Fig. 3 (see Supplementary Table S4). Three weight values were given to each of the components according to their concentrations in A. euchroma. Figure 3 illustrated the network as a bipartite graph for the components and their potential targets with color-coded nodes. In total, this network consisted of 92 nodes and 114 edges, with 55 components as hexagon nodes and 37 targets as circle nodes.

Fig. 3
figure3

Weighted compound-target bipartite network. Pink circle nodes represent potential targets, hexagon nodes remark components and each edge represent the interaction between them. Orange hexagon nodes represent Class I components, which have the highest concentration; green hexagon nodes remark Class II components, which have a lower concentration; gray hexagon nodes remark Class III components, which have without concentration data. Silver edges represent predicted CTIs and black edges remark known CTIs

To further clarify the relationships between the active compounds and their targets, we divided the 55 compounds into three groups according to their concentrations, namely Class I components, which have the highest concentration (orange node); Class II components, which have a lower concentration (green node); and Class III components, which do not have concentration information (gray node). There were 10 Class I components, 13 Class II components and 32 Class III components. Among Class I components, only MOL16 (acetylalkannin/acetylshikonin) and MOL18 (alkannin/shikonin) had interactions with more than one target; the remaining ones were associated with only one target. As for Class II components, eight components had interactions with one target, the others were associated with more than one target. NFE2L2 as the target with the most components, it could be hit by 42 components.

Analysis of KEGG pathways

37 TAL were enriched into 115 KEGG pathways. According to the scoring function, 34 pathways were qualified. We divided the 34 pathways into four categories, including four cell survival and death pathways, two immune system pathways, three endocrine system pathways and 25 specific human diseases pathways (see Supplementary Table S5). For ease of analysis, we remove 25 specific human diseases pathways that are not directly related to possible leukemia treatment mechanisms, and retain other nine pathways. We constructed a pathway-gene network based on the nine pathway-related genes by Cytoscape to identify hub genes in the pathway. The pathway-gene network contained 30 nodes and 53 edges, including 21 genes and nine pathways (Fig. 4).

Fig. 4
figure4

KEGG pathway analysis. a the score of 9 selected KEGG pathways. 9 pathways were divided into three categories, including four cell survival and death pathways (orange), two immune system pathways (cyan), three endocrine system pathways (green); b the selected 9 KEGG pathway-target interaction

There were four cell survival and death pathways include MAPK signaling pathway (hsa04010), PI3K-Akt signaling pathway (hsa04151), Sphingolipid signaling pathway (hsa04071), Apoptosis (hsa04210). These four pathways were related to cell growth, differentiation and apoptosis. MAPK signal transduction pathway plays a crucial role in cell proliferation, differentiation and other processes. As a caspase activator, MAPK also plays an important role on apoptosis. AKT1, CASP3, EGFR, HSPB1, MAPT, TNF, and TP53 participated in the MAPK signaling pathway, and 22 components could act on them. PI3K/Akt signaling pathway is extensively present in cells. It is believed to be one of the pivotal signaling pathways regulating cell growth, proliferation, differentiation and apoptosis. Akt (gene symbol as AKT1) is a direct downstream target of PI3K/Akt pathway. Activated Akt can inhibit the release of cytochrome c from the mitochondria, which blockade the apoptosis triggered by caspases. The PI3K/Akt signaling pathway enriched eight genes, including AKT1, EGFR, GSK3B, IL2, MCL1, PTEN, TLR2, TP53, and a total of 11 components could act on them. Apoptosis is also known as programmed cell death, which plays a critical role in maintaining development, homeostasis and defending against pathogens. In transmitting apoptotic signals, CASP3 is the major effect factor. The activation of CASP3 triggers apoptosis. MCL1 is a member of the Bcl-2 family, which is an apoptotic inhibitor. Inactivation of p53 functions also causes cancer cells to escape apoptosis. The p53 plays a pivotal role in leukemic hematopoiesis. A total of five genes were enriched in Apoptosis signaling pathway, they were AKT1, CASP3, MCL1, TNF, TP53, and 15 components were associated with them. Analysis of the above results indicates that induction of cell apoptosis may be a key factor in the treatment of leukemia by A. euchroma.

There were two immune system pathways, include IL-17 signaling pathway (hsa04657) and T cell receptor signaling pathway (hsa04660). The IL-17 signaling pathway were enriched by five genes, include CASP3, GSK3B, MMP1, MMP9 and TNF. As for T cell receptor signaling pathway (hsa04660), AKT1, FYN, GSK3B, IL2 and TNF were enriched. It’s well known that cytokines play a critical role in the context of inflammation, and extensive research has also proved that cytokines exert profound effects on the progression of hematopoietic malignancies including leukemia. It can be inferred that A. euchroma relieves inflammatory environment may be related to its treatment of leukemia.

The other three pathways were endocrine system-related pathways, include Estrogen signaling pathway (hsa04915), Thyroid hormone signaling pathway (hsa04919) and Relaxin signaling pathway (hsa04926). A total of 12 genes enriched by these three pathways, they were AKT1, EGFR, ESR1, GSK3B, HDAC3, HIF1A, MMP1, MMP2, MMP9, RARA, SP1 and TP53. Of which, EGFR, one of the transmembrane tyrosine kinase receptors, serves as a stimulus for cancer growth, which is aberrantly activated in most common tumors, including leukemia.

Molecular docking analysis and critical compound-target bipartite network

According to the scoring function, 10 hub targets (AKT1, CASP3, EGFR, GSK3B, HIF1A, MMP2, MMP9, PTEN, TNF and TP53) were qualified, and 18 components formed interactions with them in all CTIs. These components include arnebinone (MOL6), shikonofuran b (MOL8), acetylalkannin/acetylshikonin (MOL16), alkannin/shikonin (MOL18), arnebiabinone (MOL21), dehydroalkannin (MOL24), oleanolic acid (MOL47), shinjulactone k (MOL65), ursolic acid (MOL49), o9-angeloylretronecine (MOL41), caffeic acid (MOL42), linoleic acid (MOL53), daucosterol (MOL38), beta-sitosterol (MOL36), echimidine (Mol39), alkannin angelate (Mol19), 3-hydroxy-3-methyl butyric acid (Mol51), and tormentic acid (Mol48).

There were 27 CTIs between the 18 components and 10 targets, and 11 of them were KCTIs, 16 were PCTIs. All of the 16 PCTIs were composed of CTIs related to HIF1A and TNF while all of the CTIs related to HIF1A and TNF were PCTIs. In this study, we used the molecular docking stimulation to identify the binding ability of the 16 pairs of PCTIs that between components of A. euchroma and the obtained hub genes.

As the above described, only HIF1A and TNF-related CTIs were predicted, the other hub genes and their associated CTIs were KCTIs. Therefore, the molecular docking simulation was further performed to determine the binding level between the two predicted targets (HIF1A and TNF) and their corresponding components containing in A. euchroma. The docking scores of 16 pairs of PCTIs were listed in Table S6. There were eight pairs of PCTIs had strong binding free energy, since their docking scores were higher than the median value of all pairs. The eight pairs of PCTIs consisted of six pairs of TNF-related PCTIs and two of HIF1A. The action modes between TNF and its corresponding compounds were shown in Fig.5a-f, the interplay between HIF1A and its corresponding compounds were shown in Fig.5g-h. It can be found that dehydroalkannin (MOL24) could form a hydrogen bond with Tyr151, oleanolic acid (MOL47) could form a hydrogen bond with Ser60. In addition, dehydroalkannin (MOL24) could form a π-π interaction with Tyr59. Besides, arnebinone (MOL6) had the best binding ability with TNF (docking score = − 6.222), and followed by arnebiabinone (MOL21) (docking score = − 5.658), dehydroalkannin (MOL24) (docking score = − 5.405).

Fig. 5
figure5

The interactive modes of 2 targets and the corresponding 8 compounds. a arnebinone to TNF (docking score = − 6.222); b arnebiabinone to TNF (docking score = − 5.658); c dehydroalkannin to TNF (docking score = − 5.405); d oleanolic acid to TNF (docking score = − 4.875); e shinjulactone k to TNF (docking score = − 4.875); f ursolic acid to TNF (docking score = − 4.789); g shikonofuran b to HIF1A (docking score = − 4.802); h o9-angeloylretronecine to HIF1A (docking score = − 4.488)

The eight pairs of PCTIs and 11 pairs of KCTIs were used to critical compound-target bipartite network, which involved in 10 targets and 14 components, as shown in Fig. 6. Here, two (MOL6 and MOL8) components were monoterpene phenol and benzoquinones (green nodes), four (MOL24, MOL21, MOL18 and MOL16) were naphthoquinones (gold node), three (MOL65, MOL47 and MOL49) were triterpenoids (blue nodes), one (MOL41) was alkaloids (cyan node), two (MOL38 and MOL36) were steroids (purple nodes), and the other two (MOL42 and MOL53) were organic acids (brown nodes). Structures of the 14 components were shown in Supplementary Figure S2.

Fig. 6
figure6

Critical compound-target bipartite network. Pink circle nodes represent critical targets, hexagon nodes remark components and each edge represent the interaction between them. Orange nodes represent naphthoquinones; green nodes remark monoterpene phenol and benzoquinones; cyan nodes remark alkaloids; purple nodes represent steroids; brown nodes remark organic acids; blue nodes represent triterpenoids. Silver edges represent predicted CTIs and black edges remark known CTIs

Discussion

A. euchroma has been used to treat blood diseases, including leukemia in the clinical of Chinese medicine [3]. In recent years, more and more studies have indicated that shikonin and some of its derivatives can inhibit leukemia cell proliferation, induce leukemia cell differentiation and promote apoptosis [9,10,11]. Most of the studies focus on individual or several ingredients, and the effects of all components of A. euchroma have never been explored systematically. Network pharmacology focuses on “multi-constituents, multi-targets to treat diseases”, which coincides with the holistic and systematic concepts of TCM [33]. Moreover, although the importance of TCM dosage has received widespread attention, few TCM-related network pharmacology studies have considered the influence of ingredient concentration on the effectiveness of TCM. In this study, we employed network pharmacology approach to explore the MoA of A. euchroma on leukemia. We retrieved 65 components and 157 targets of A. euchroma from known databases and literature. Combined with network analysis and molecular docking simulation, we found that 14 components and 10 targets were potential critical components and targets, which might contribute to the effect of A. euchroma on leukemia treatment. The 14 components include alkannin/shikonin (MOL18), arnebinone (MOL6), shikonofuran b (MOL8), arnebiabinone (MOL21), acetylalkannin/acetylshikonin (MOL16), dehydroalkannin (MOL24), oleanolic acid (MOL47), shinjulactone k (MOL65), ursolic acid (MOL49), o9-angeloylretronecine (MOL41), caffeic acid (MOL42), linoleic acid (MOL53), daucosterol (MOL38), beta-sitosterol (MOL36). The 10 targets include AKT1, CASP3, EGFR, GSK3B, HIF1A, MMP2, MMP9, PTEN, TNF and TP53. Our results indicated that the possible MoA might be involved in three aspects: 1. A. euchroma inhibits leukemia cell survival and induces apoptosis; 2. A. euchroma relieves inflammatory environment; 3. A. euchroma inhibits angiogenesis.

Apoptosis, also known as programmed cell death, is a physiological process or some pathological condition, in which the cells take part in the death process after they are triggered by a certain signal. Inhibiting leukemia cell survival and inducing leukemia cell apoptosis is a promising therapeutic approach for leukemia. There were four pathways (MAPK signaling pathway (hsa04010), PI3K-Akt signaling pathway (hsa04151), Sphingolipid signaling pathway (hsa04071), Apoptosis (hsa04210)) were related to cell survival and apoptosis, and 13 targets, 26 components in weighted compound-target bipartite network were mapped onto the four pathways. When transmitting apoptotic signals, CASP3 is the major effect factor. The activation of CASP3 triggers apoptosis. As one of the caspase activators, MAPK is well known to play a crucial role in apoptosis. Extensive studies have shown that targeting to MAPK signaling cascades, alone or in combination with other drugs, results in enhanced anti-leukemic responses in AML [34, 35]. It has been shown that activation of p38 MAPK promoted BCL2 degradation, which partly induced K562 cell apoptosis [36]. Huang et al. reported that shikonin treatment activated p38 MAPK and JNK phosphorylation in human lens epithelial cells [37]. The PI3K/Akt signaling pathway regulates diverse cellular processes, including growth, proliferation, differentiation and apoptosis by phosphorylating its downstream target, including BCL2L11, BAD, MDM2, XIAP, CASP9, GSK3B and so on. TP53 is a tumor suppressor protein and regulator of the cell cycle, which plays a pivotal role in leukemic hematopoiesis, but its functions are frequently suppressed by MDM2. Therefore, inhibition of MDM2 is beneficial to the stability of TP53. However, Akt enhances MDM2 activity, which is harmful for the stability of TP53. PTEN is a main negative regulator of the PI3K/Akt pathway, and plays an important role in maintaining hematopoietic stem cells and preventing leukemia [38]. It has been demonstrated that impaired PTEN function and aberrantly activated PI3K/Akt were present in many leukemia cases [9, 39]. Pan et al. suggested that combined BCL2 inhibition and TP53 activation may be a promising therapeutic approach for AML [40].

It is well known that inflammation is one of the major barriers to cancer therapy. In recent years, it has become evident that inflammation plays a key role in leukemia [41]. There were two pathways, including IL-17 signaling pathway (hsa04657) and T cell receptor signaling pathway (hsa04660) were related to inflammation. A total of eight targets were enriched into inflammation-related pathway, they were TNF, IL2, CASP3, GSK3B, MMP1, MMP9, AKT1 and FYN. And these eight targets had interactions with 13 components in the weighted compound-target bipartite network. As a pro-inflammatory cytokine, TNF is recognized as a key mediator of inflammatory reactions in tumor tissues, and also responsible for increased NF-κB activity in many tumors. TNF promotes the progression of several hematopoietic tumors by participating in several signaling pathways, including NF-κB, and PI3K/Akt [42, 43]. It has been reported that TNF can be produced by various leukemia cells, including AML [34], ALL [44], CML [45], CLL [46], and so on. There was also evidence showing that TNF induced MMP9 expression or secretion in leukemia cells [47]. TNF has been used as an anti-cancer drugs in various cancer cells. Staniforth et al. found that shikonin (MOL18), isobutyrylshikonin (MOL26), acetylshikonin (MOL16), β, β-dimethylacrylshikonin (MOL31) and isovalerylshikonin (MOL27) showed significant dose-dependent inhibition of TNF in vivo [48].

Angiogenesis is a critical element to cancer cell survival, and increased vascularity was found in patients with AML [49]. HIF1A, a key regulator of the cellular response to hypoxia, controls a vast array of gene products involved in energy metabolism, glycolysis, angiogenesis, apoptosis and cell cycle [50]. HIF1A has been recognized as a strong promoter of tumor growth and it is responsible for VEGF gene expression [51]. VEGF is one of the major angiogenesis-activating protein, it has been implicated in leukemia-associated angiogenesis [52]. There were some data showed that HIF1A over-expressed in some leukemia cells [53, 54]. Frolova et al. reported that blockade of HIF1A-mediated signaling may enhance the efficacy of the therapeutic regiments in ALL [55]. As one of the important factors to promote angiogenesis, matrix metalloproteinase (MMPs) were found highly expressed in B-CLL cells, and promoted the migration and invasion of leukemia cells [56]. MMP9, secreted by leukemic cells, mediates opening of the blood–brain barrier by disrupting tight junction proteins in CNS leukemia [57]. Sustained or enhanced MMP9 secretion plays an important role in the pathophysiology of tumor progression.

On the basis of the targets of A. euchroma discovered by network pharmacology, molecular docking simulation was used to explore the binding ability of between components and proteins. This technique is structure-based and could help researchers discover the interactions between the components of TCM and network targets. In this study, HIF1A and TNF were selected for molecular docking studies because they were predicted targets for components. TNF was involved in five pathways, include IL-17 signaling pathway (hsa04657), MAPK signaling pathway (hsa04010), T cell receptor signaling pathway (hsa04660), Sphingolipid signaling pathway (hsa04071) and Apoptosis (hsa04210), which were related to both apoptosis and inflammation. The molecular docking results indicated that six components (MOL6, MOL21, MOL24, MOL47, MOL65, and MOL49) had strong binding free energy with TNF, and two components (MOL8 and MOL41) had good binding affinity with HIF1A, since their docking scores were higher than the median value of all pairs. Among the eight components, MOL6 and MOL8 were monoterpene phenol and benzoquinones, MOL24 and MOL21 were naphthoquinones, which were the two main types of active components in A. euchroma.

According to the above analysis, we selected 10 targets and 14 components as critical targets and key components. Among these, shikonin (MOL18) was the most studied components of A. euchroma, it was demonstrated to have anti-inflammatory and pro-apoptotic effects [10, 58]. Acetylshikonin (MOL16) was found to induce tumor cell apoptosis through activating the pro-apoptotic bcl-2 family and caspase-3 [59].

As one of the most famous Chinese medicine for treating leukemia, ATO is an effective and relatively safe drug in the treatment of APL (acute promyelocytic leukemia, a subtype of AML) [2]. Combining compound-target network of A. euchroma and related literature, we compared ATO and A. euchroma in terms of leukemia treatment, and we found some commonalities: (1) MAPK signaling pathway. ATO has been demonstrated to activate all three MAPKs pathways, ERKs, JNKs, and p38 kinases to exert its anti-leukemia effect [60]. The MAPK signaling pathway was also included in pathways that enriched by targets of A. euchroma. (2) Angiogenesis. ATO was reported to inhibit angiogenesis by inhibiting the production of VEGF in leukemic cell [61]. As for A. euchroma, seven components could act on HIF1A, which is responsible for increased expression of VEGF. Moreover, several studies have shown that some components of A. euchroma exhibited antiangiogenic activity by inhibiting VEGF [62, 63]. (3) Caspases. ATO was shown to increase caspases activity to induce apoptosis of leukemia cells [64]. Extensive studies have demonstrated that components of A. euchroma increased the expression of caspases to induce apoptosis [10, 65]. As a result, we have reasons to speculate that A. euchroma may be a potential anti-leukemia agent with less toxicity than ATO.

In this study, we found that the MoA of A. euchroma on leukemia involves the inhibition of leukemia cell survival and induction of leukemia cell apoptosis, the relief of inflammatory environment, and the inhibition of angiogenesis. Concretely speaking, it includes three aspects: (1) A. euchroma impacts the PI3K/AKT and MAPK signaling pathway, which inhibit leukemia cell survival and induce leukemia cell apoptosis. (2) A. euchroma inhibits the expression of inflammatory medium, such as TNF, IL6, IL2 and so on to relieve inflammatory environment. (3) A. euchroma inhibits angiogenesis by regulating HIF1A, MMP2, MMP9 and other factors related to angiogenesis. The results preliminarily validated and explained the therapeutic material basis and mechanism of A. euchroma on leukemia, which provided a preliminary information and basis for further in-depth exploration of its MoA, and a reference for the study of the more complex components of TCM.

Conclusions

Considering the concentration information of components in A. euchroma and combined with the methods of network analysis and molecular docking simulation, we found that 55 active components and 37 targets might be related to the anti-leukemia effects of A. euchroma, of which 14 components (MOL6, MOL8, MOL16, MOL18, MOL21, MOL24, MOL47, MOL65, MOL49, MOL41, MOL42, MOL53, MOL38, and MOL36) and 10 targets (AKT1, CASP3, EGFR, GSK3B, HIF1A, MMP2, MMP9, PTEN, TNF, and TP53) were considered to be critical. It was demonstrated that A. euchroma has potential anti-leukemia activities because of its effects on inhibiting leukemia cell survival and inducing apoptosis, relieving inflammatory environment and inhibiting angiogenesis. We further realized that, not only shikonin, the most fully studied ingredient of A. euchroma, one of compound of naphthoquinones, other components among the 14 critical components of A. euchroma, including monoterpene phenol and benzoquinones, triterpenoids, alkaloids, steroids and organic acids would also have therapeutic effects on leukemia.

Availability of data and materials

The data can be requested from the author.

Abbreviations

A. euchroma :

Arnebia euchroma

ALL:

Acute lymphoblastic leukemia

AML:

Acute myeloid leukemia

APL:

Acute promyelocytic leukemia

ATO:

Arsenic trioxide

BC:

Betweenness centrality

bSDTNBI:

balanced substructure-drug-target network-based inference

CC:

Closeness centrality

CLL:

Chronic lymphoblastic leukemia

CML:

Chronic myeloid leukemia

CTIs:

Compound-target interactions

DC:

Degree centrality

EC:

Eigenvector centrality

FDR:

False Discovery Rate

IC:

Information centrality

KCTIs:

Known compound-target interactions

LAC:

Local average connectivity-based method

MMPs:

Matrix metalloproteinase

MoA:

Mechanism of action

NBI:

Network-based inference

NC:

Network centrality

PCTIs:

Predicted compound-target interactions

PDB:

RCSB Protein Data Bank

PPIA :

A. euchroma PPI network

PPIL :

Leukemia PPI network

PPIs:

Protein-protein interactions

SC:

Subgraph centrality

SDTNBI:

Substructure-drug-target network-based inference

TA :

Targets of components of A. euchroma

TAL :

Leukemia-related targets of A. euchroma

TCM:

Traditional Chinese medicine

TL :

Targets related to leukemia

References

  1. 1.

    Bruserud Ø, Gjertsen BT, Huang T. Induction of differentiation and apoptosis— a possible strategy in the treatment of adult acute Myelogenous leukemia. Oncologist. 2000;5:454–62.

    CAS  Google Scholar 

  2. 2.

    Lengfelder E, Hofmann WK, Nowak D. Impact of arsenic trioxide in the treatment of acute promyelocytic leukemia. Leukemia. 2012;26:433–42. https://doi.org/10.1038/leu.2011.245.

    CAS  Article  Google Scholar 

  3. 3.

    Ma S, Geng Y, Ma L, Zhu L. Advances in Studies on Medicinal Radix Arnebia. Mod Chinese Med. 2020. https://doi.org/10.13313/j.issn.1673-4890.20190719004.

  4. 4.

    Shi L. Data analysis traditional Chinese medicine prescription combination rules of blood diseases treatment. J Basic Chin Med. 2016;22:853–5.

    Google Scholar 

  5. 5.

    Wang RB, Yin RT, Zhou W, Xu DF, Li SS. Shikonin and its derivatives: a patent review. Expert Opin Ther Pat. 2012;22:977–97.

    CAS  Google Scholar 

  6. 6.

    Papageorgiou VP, Assimopoulou AN, Couladouros EA, Hepworth D, Nicolaou KC. The chemistry and biology of alkannin, shikonin, and related naphthazarin natural products. Angew Chemie Int Ed. 1999;38:270–300.

    Google Scholar 

  7. 7.

    Hu Y, Jiang Z, Leung KSY, Zhao Z. Simultaneous determination of naphthoquinone derivatives in boraginaceous herbs by high-performance liquid chromatography. Anal Chim Acta. 2006;577:26–31.

    CAS  Google Scholar 

  8. 8.

    Zan K, Zheng H, Liu J, Guo L, Zheng J, Ma S. Simultaneous determination of eight hydroxyl naphthoquinones in different parts of Arnebiae euchroma by HPLC. China J Chinese Mater Medica. 2017;42:2532–7. https://doi.org/10.19540/j.cnki.cjcmm.20170419.001.

    Article  Google Scholar 

  9. 9.

    Chen Y, Wang T, Du J, Li Y, Wang X, Zhou Y, et al. The critical role of PTEN/PI3K/AKT signaling pathway in shikonin-induced apoptosis and proliferation inhibition of chronic myeloid leukemia. Cell Physiol Biochem. 2018;47:981–93.

    CAS  Google Scholar 

  10. 10.

    Mao X, Yu C, Li W, Li W. Induction of apoptosis by shikonin through a ROS/JNK-mediated process in Bcr/Abl-positive chronic myelogenous leukemia (CML) cells. Cell Res. 2008;18:879–88.

    CAS  Google Scholar 

  11. 11.

    Zhang B, Chen N, Chen H, Wang Z, Zheng Q. The critical role of redox homeostasis in Shikonin-induced HL-60 cell differentiation via unique modulation of the Nrf2/ARE pathway. Oxidative Med Cell Longev. 2012;2012:781516.

    Google Scholar 

  12. 12.

    Hopkins AL. Network pharmacology. Nat Biotechnol. 2007;25:1110–1.

    CAS  Google Scholar 

  13. 13.

    Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008;4:682–90.

    CAS  Google Scholar 

  14. 14.

    Cheng F, Liu C, Jiang J, Lu W, Li W, Liu G, et al. Prediction of drug-target interactions and drug repositioning via network-based inference. PLoS Comput Biol. 2012;8:1–12.

    Google Scholar 

  15. 15.

    Wu Z, Lu W, Wu D, Luo A, Bian H, Li J, et al. In silico prediction of chemical mechanism of action via an improved network-based inference method. Br J Pharmacol. 2016;173:3372–85.

    CAS  Google Scholar 

  16. 16.

    Wu Z, Cheng F, Li J, Li W, Liu G, Tang Y. SDTNBI: an integrated network and chemoinformatics tool for systematic prediction of drug–target interactions and drug repositioning. Brief Bioinform. 2017;18:333–47.

    CAS  Google Scholar 

  17. 17.

    Xiao X, Chen Z, Wu Z, Wang T, Li W, Liu G, et al. Insights into the antineoplastic mechanism of Chelidonium majus via systems pharmacology approach. Quant Biol. 2019;7:42–53.

    CAS  Google Scholar 

  18. 18.

    Wang T, Wu Z, Sun L, Li W, Liu G, Tang Y. A computational systems pharmacology approach to investigate molecular mechanisms of herbal formula Tian-Ma-Gou-Teng-Yin for treatment of Alzheimer’s disease. Front Pharmacol. 2018;9:1–16.

    Google Scholar 

  19. 19.

    Wang Y, Li J, Wu Z, Zhang B, Yang H, Wang Q, et al. Insights into the molecular mechanisms of Polygonum multiflorum Thunb-induced liver injury: a computational systems toxicology approach. Acta Pharmacol Sin. 2017;38:719–32. https://doi.org/10.1038/aps.2016.147.

    CAS  Article  Google Scholar 

  20. 20.

    Zhang R, Yu S, Bai H, Ning K. TCM-Mesh: the database and analytical system for network pharmacology analysis for TCM preparations. Sci Rep. 2017;7:1–14.

    Google Scholar 

  21. 21.

    Xue R, Fang Z, Zhang M, Yi Z, Wen C, Shi T. TCMID: traditional Chinese medicine integrative database for herb molecular mechanism analysis. Nucleic Acids Res. 2013;41:1089–95.

    Google Scholar 

  22. 22.

    Gilson MK, Liu T, Baitaluk M, Nicola G, Hwang L, Chong J. BindingDB in 2015: a public database for medicinal chemistry, computational chemistry and systems pharmacology. Nucleic Acids Res. 2016;44:D1045–53.

    CAS  Google Scholar 

  23. 23.

    Harding SD, Sharman JL, Faccenda E, Southan C, Pawson AJ, Ireland S, et al. The IUPHAR/BPS guide to PHARMACOLOGY in 2018: updates and expansion to encompass the new guide to IMMUNOPHARMACOLOGY. Nucleic Acids Res. 2018;46:D1091–106.

    CAS  Google Scholar 

  24. 24.

    Gaulton A, Hersey A, Nowotka ML, Patricia Bento A, Chambers J, Mendez D, et al. The ChEMBL database in 2017. Nucleic Acids Res. 2017;45:D945–54.

    CAS  Google Scholar 

  25. 25.

    Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Bryant SH. PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009;37(SUPPL. 2):623–33.

    Google Scholar 

  26. 26.

    Amberger J, Bocchini CA, Scott AF, Hamosh A. McKusick’s Online Mendelian Inheritance in Man (OMIM). Nucleic Acids Res. 2009;37 Database:D793–6. https://doi.org/10.1093/nar/gkn665.

    CAS  Article  Google Scholar 

  27. 27.

    Whirl-Carrillo M, McDonagh EM, Hebert JM, Gong L, Sangkuhl K, Thorn CF, et al. Pharmacogenomics knowledge for personalized medicine. Clin Pharmacol Ther. 2012;92:414–7.

    CAS  Google Scholar 

  28. 28.

    Li Y, Yu C, Li X, Zhang P, Tang J, Yang Q, et al. Therapeutic target database update 2018: enriched resource for facilitating bench-to-clinic research of targeted therapeutics. Nucleic Acids Res. 2018;46:D1121–7.

    CAS  Google Scholar 

  29. 29.

    Piñero J, Bravo Á, Queralt-Rosinach N, Gutiérrez-Sacristán A, Deu-Pons J, Centeno E, et al. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017;45:D833–9.

    Google Scholar 

  30. 30.

    Rappaport N, Twik M, Plaschkes I, Nudel R, Stein TI, Levitt J, et al. MalaCards: an amalgamated human disease compendium with diverse clinical and genetic annotation and structured search. Nucleic Acids Res. 2017;45:D877–87.

    CAS  Google Scholar 

  31. 31.

    Song W, Ni S, Fu Y, Wang Y. Uncovering the mechanism of maxing Ganshi decoction on asthma from a systematic perspective: a network pharmacology study. Sci Rep. 2018;8:1–11.

    Google Scholar 

  32. 32.

    Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.

    CAS  Google Scholar 

  33. 33.

    Li S, Fan TP, Jia W, Lu A, Zhang W. Network pharmacology in traditional chinese medicine. Evid Based Complement Altern Med. 2014;2014:10–2.

    Google Scholar 

  34. 34.

    Volk A, Li J, Xin J, You D, Zhang J, Liu X, et al. Co-inhibition of NF-κB and JNK is synergistic in TNF-expressing human AML. J Exp Med. 2014;211:1093–108. https://doi.org/10.1084/jem.20130990.

    CAS  Article  Google Scholar 

  35. 35.

    Geest CR, Coffer PJ. MAPK signaling pathways in the regulation of hematopoiesis. J Leukoc Biol. 2009;86:237–50.

    CAS  Google Scholar 

  36. 36.

    Changchien JJ, Chen YJ, Huang CH, Cheng TL, Lin SR, Chang L. Sen. Quinacrine induces apoptosis in human leukemia K562 cells via p38 MAPK-elicited BCL2 down-regulation and suppression of ERK/c-Jun-mediated BCL2L1 expression. Toxicol Appl Pharmacol. 2015;284:33–41. https://doi.org/10.1016/j.taap.2015.02.005.

    CAS  Article  Google Scholar 

  37. 37.

    Huang W, Zhang Y, Tang X. Shikonin inhibits the proliferation of human lens epithelial cells by inducing apoptosis through ROS and caspase-dependent pathway. Molecules. 2014;19:7785–97. https://doi.org/10.3390/molecules19067785.

    CAS  Article  Google Scholar 

  38. 38.

    Zhang J, Grindley JC, Yin T, Jayasinghe S, He X, Ross JT, et al. PTEN maintains haematopoietic stem cells and acts in lineage choice and leukaemia prevention. Nature. 2006;441:518–22.

    CAS  Google Scholar 

  39. 39.

    Lamy T, Moignet A, Loughran TP. LGL leukemia: from pathogenesis to treatment. Blood. 2017;129:1082–94.

    CAS  Google Scholar 

  40. 40.

    Pan R, Ruvolo V, Mu H, Leverson JD, Nichols G, Reed JC, et al. Synthetic Lethality of Combined Bcl-2 Inhibition and p53 Activation in AML: Mechanisms and Superior Antileukemic Efficacy. Cancer Cell. 2017;32:748–760.e6. https://doi.org/10.1016/j.ccell.2017.11.003.

    CAS  Article  Google Scholar 

  41. 41.

    Binder S, Luciano M, Horejs-hoeck J. The cytokine network in acute myeloid leukemia ( AML ): A focus on pro- and anti-inflammatory mediators. Cytokine Growth Factor Rev. 2018;43:8–15. https://doi.org/10.1016/j.cytogfr.2018.08.004.

    CAS  Article  Google Scholar 

  42. 42.

    Gu L, Findley HW, Zhu N, Zhou M. Endogenous TNFα mediates cell survival and chemotherapy resistance by activating the PI3K/Akt pathway in acute lymphoblastic leukemia cells. Leukemia. 2006;20:900–4.

    CAS  Google Scholar 

  43. 43.

    Tsai HJ, Kobayashi S, Izawa K, Ishida T, Watanabe T, Umezawa K, et al. Bioimaging analysis of nuclear factor-κB activity in Philadelphia chromosome-positive acute lymphoblastic leukemia cells reveals its synergistic upregulation by tumor necrosis factor-α-stimulated changes to the microenvironment. Cancer Sci. 2011;102:2014–21.

    CAS  Google Scholar 

  44. 44.

    Potapnev MP, Petyovka NV, Belevtsev MV, Savitskiy VP, Migal NV. Plasma level of tumor necrosis factor-alpha (TNF-α) correlates with leukocytosis and biological features of leukemic cells, but not treatment response of children with acute lymphoblastic leukemia. Leuk Lymphoma. 2003;44:1077–9.

    CAS  Google Scholar 

  45. 45.

    Gallipoli P, Pellicano F, Morrison H, Laidlaw K, Allan EK, Bhatia R, et al. Autocrine TNF-α production supports CML stem and progenitor cell survival and enhances their proliferation. Blood. 2013;122:3335–9.

    CAS  Google Scholar 

  46. 46.

    Lech-Maranda E, Grzybowska-Izydorczyk O, Wyka K, Mlynarski W, Borowiec M, Antosik K, et al. Serum tumor necrosis factor-α and interleukin-10 levels as markers to predict outcome of patients with chronic lymphocytic leukemia in different risk groups defined by the IGHV mutation status. Arch Immunol Ther Exp (Warsz). 2012;60:477–86.

    CAS  Google Scholar 

  47. 47.

    Heidinger M, Kolb H, Krell HW, Jochum M, Ries C. Modulation of autocrine TNF-α-stimulated matrix metalloproteinase 9 (MMP-9) expression by mitogen-activated protein kinases in THP-1 monocytic cells. Biol Chem. 2006;387:69–78.

    CAS  Google Scholar 

  48. 48.

    Staniforth V, Wang S, Shyur LF, Yang N. Shikonins, Phytocompounds from Lithospermum erythrorhizon, inhibit the transcriptional activation of human tumor necrosis factor α promoter in vivo. J Biol Chem. 2004;279:5877–85.

    CAS  Google Scholar 

  49. 49.

    Padró T, Ruiz S, Bieker R, Bürger H, Steins M, Kienast J, et al. Increased angiogenesis in the bone marrow of patients with acute myeloid leukemia. Blood. 2000;95:2637–44.

    Google Scholar 

  50. 50.

    Semenza GL. Hypoxia-inducible factor 1 and cardiovascular disease. Annu Rev Physiol. 2014;76:39–56.

    CAS  Google Scholar 

  51. 51.

    Forsythe JA, Jiang BH, Iyer NV, Agani F, Leung SW, Koos RD, et al. Activation of vascular endothelial growth factor gene transcription by hypoxia-inducible factor 1. Mol Cell Biol. 1996;16:4604–13.

    CAS  Google Scholar 

  52. 52.

    Bellamy WT, Richter L, Sirjani D, Roxas C, Glinsmann-Gibson B, Frutiger Y, et al. Vascular endothelial cell growth factor is an autocrine promoter of abnormal localized immature myeloid precursors and leukemia progenitor formation in myelodysplastic syndromes. Blood. 2001;97:1427–34.

    CAS  Google Scholar 

  53. 53.

    Mayerhofer M, Valent P, Sperr WR, Griffin JD, Sillaber C. BCR/ABL induces expression of vascular endothelial growth factor and its transcriptional activator, hypoxia inducible factor-1alpha, through a pathway involving phosphoinositide 3-kinase and the mammalian target of rapamycin. Blood. 2002;100:3767–75.

    CAS  Google Scholar 

  54. 54.

    Wang Y, Liu Y, Malek SN, Zheng P, Liu Y. Targeting HIF1α eliminates cancer stem cells in hematological malignancies. Cell Stem Cell. 2011;8:399–411.

    CAS  Google Scholar 

  55. 55.

    Frolova O, Samudio I, Benito J, Jacamo R, Kornblau SM, Markovic A, et al. Regulation of HIF-1α signaling and chemoresistance in acute lymphocytic leukemia under hypoxic conditions of the bone marrow microenvironment. Cancer Biol Ther. 2012;13:858–70.

    CAS  Google Scholar 

  56. 56.

    Redondo-muñoz J, Escobar-díaz E, Samaniego R, Terol MJ, José A, Dc W, et al. MMP-9 in B-cell chronic lymphocytic leukemia is up-regulated by α4β1 integrin or CXCR4 engagement via distinct signaling pathways, localizes to podosomes, and is involved in cell invasion and migration. Blood. 2006;108:3143–51.

    Google Scholar 

  57. 57.

    Feng S, Cen J, Huang Y, Shen H, Yao L, Wang Y, et al. Matrix metalloproteinase-2 and -9 secreted by leukemic cells increase the permeability of blood-brain barrier by disrupting tight junction proteins. PLoS One. 2011;6:e20599.

    CAS  Google Scholar 

  58. 58.

    Yang C, Liu P, Wang S, Zhao G, Zhang T, Guo S, et al. Shikonin exerts anti-inflammatory effects in LPS-induced mastitis by inhibiting NF-κB signaling pathway. Biochem Biophys Res Commun. 2018;505:1–6. https://doi.org/10.1016/j.bbrc.2018.08.198.

    CAS  Article  Google Scholar 

  59. 59.

    Xiong W, Luo G, Zhou L, Zeng Y, Yang W. In vitro and in vivo antitumor effects of acetylshikonin isolated from Arnebia euchroma (Royle) Johnst (Ruanzicao) cell suspension cultures. Chin Med. 2009;4:1–7.

    Google Scholar 

  60. 60.

    Bode AM, Dong Z. The paradox of arsenic: molecular mechanisms of cell transformation and chemotherapeutic effects. Crit Rev Oncol Hematol. 2002;42:5–24.

    Google Scholar 

  61. 61.

    Roboz GJ, Dias S, Lam G, Lane WJ, Soignet SL, Warrell J, et al. Arsenic trioxide induces dose- and time-dependent apoptosis of endothelium and may exert an antileukemic effect via inhibition of angiogenesis. Blood. 2000;96:1525–30.

    CAS  Google Scholar 

  62. 62.

    Lee H-J, Lee H-J, Magesh V, Nam D, Lee E-O, Ahn SK, et al. Shikonin, acetylshikonin, and isobutyroylshikonin inhibit VEGF-induced angiogenesis and suppress tumor growth in Lewis lung carcinoma-bearing mice. Yakugaku Zasshi. 2008;128:1681–8. https://doi.org/10.1248/yakushi.128.1681.

    CAS  Article  Google Scholar 

  63. 63.

    Xu Y, Xu X, Gao X, Chen H, Geng L. Shikonin suppresses IL-17-induced VEGF expression via blockage of JAK2/STAT3 pathway. Int Immunopharmacol. 2014;19:327–33. https://doi.org/10.1016/j.intimp.2014.01.027.

    CAS  Article  Google Scholar 

  64. 64.

    Perkins C, Kim CN, Fang G, Bhalla KN. Arsenic induces apoptosis of multidrug-resistant human myeloid leukemia cells thin express Bcr-Abl or overexpress MDR, MRP, Bcl-2, or Bcl-x(L). Blood. 2000;95:1014–22.

    CAS  Google Scholar 

  65. 65.

    Xiong Y, Ma X, Zhang Z, Shao Z, Zhang Y, Zhou L. Apoptosis induced by beta,beta-dimethylacrylshikonin is associated with Bcl-2 and NF-kappaB in human breast carcinoma MCF-7 cells. Oncol Lett. 2013;6:1789–93. https://doi.org/10.3892/ol.2013.1613.

    CAS  Article  Google Scholar 

Download references

Acknowledgments

n/a

Funding

This work was supported by the National Natural Science Foundation of China (Grants U1603122 and 81673356) and the 111 Project (Grant BP0719034).

Author information

Affiliations

Authors

Contributions

YT, JW, BZ and BW contributed to conception and design of the study. BW performed the experiments and wrote the manuscript. ZW provided the method of target prediction. WL and GL contributed to the writing of this article. YT implemented the study and modified the manuscript. All authors have discussed the results and approved the final manuscript.

Corresponding authors

Correspondence to Zengrui Wu or Yun Tang.

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.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1.

Figure S1. Two configurations of naphthoquinones.

Additional file 2.

Figure S2. Structures of the 14 critical components.

Additional file 3

Table S1. Components of A. euchroma and content score.

Additional file 4.

Table S2. Known and predicted component-target interactions.

Additional file 5.

Table S3. Eight central attributes of targets.

Additional file 6.

Table S4. Weighted component-important target interactions.

Additional file 7.

Table S5. 34 KEGG pathways.

Additional file 8.

Table S6. Docking scores of 16 pairs of PCTIs.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, B., Wu, Z., Wang, J. et al. Insights into the mechanism of Arnebia euchroma on leukemia via network pharmacology approach. BMC Complement Med Ther 20, 322 (2020). https://doi.org/10.1186/s12906-020-03106-z

Download citation

Keywords

  • Arnebia euchroma
  • Concentration
  • Leukemia
  • Traditional Chinese medicine
  • Network pharmacology
  • Weighted network