System Prediction and Validation of TCM for Chronic Myeloid Leukemia Treatment from the Perspective of Low-Toxicity Chemotherapy: A Stilbene α-Viniferin Has a Proapoptotic Effect on K562 Cells via the Mitochondrial Pathway

Objective In traditional Chinese medicine (TCM), chronic myeloid leukemia (CML) has been attributed to “poisoned bone marrow,” which is viewed as a loss of Qi or blood, a deficiency in Yin or Yang that causes a complex imbalance between cell growth and death. Malignant myeloid progenitor cells display excessive growth that is difficult to control without toxicity. More than 60 herbs in TCM have shown efficacy against CML. However, the key molecules and mechanisms involved in the holistic-level characterization, as well as the effective target associations, are still unknown. Methods The present study employed a computational approach with filtering potential compounds via admetSAR, systems biology-based functional data prediction, and biochemical and molecular biological validation. Results We generated 118 bioactive compounds from 11 herbs within four dialectical therapy groups that are most commonly used to treat CML and predicted 141 potential targets. The stilbene resveratrol and its derivatives were found to be highly related to these targets. Among them, α-viniferin was predicted to target Bcl-2, caspase-3, 8, and 9, MAPK14, CDK2, HSP90AA1, and others, reflecting CML therapeutic strategies. In vitro, experimental data showed a nonnecrotic growth limitation of K562 cells caused by α-viniferin was predicted to target Bcl-2, caspase-3, 8, and 9, MAPK14, CDK2, HSP90AA1, and others, reflecting CML therapeutic strategies. μg·mL−1 at 24 h. Finally, we validated the chemotherapeutic effect of α-viniferin was predicted to target Bcl-2, caspase-3, 8, and 9, MAPK14, CDK2, HSP90AA1, and others, reflecting CML therapeutic strategies. Conclusions Our work sheds light on the mechanism of the efficacy of the stilbene α-viniferin in TCM for the prevention of CML. This work also predicts and validates targets in the mitochondrial signaling pathway, providing a novel strategy for CML treatment.α-viniferin was predicted to target Bcl-2, caspase-3, 8, and 9, MAPK14, CDK2, HSP90AA1, and others, reflecting CML therapeutic strategies.


Introduction
Chronic myelogenous leukemia (CML) is a hematologic malignancy of stem and progenitor cells that accounts for 15% of adult leukemia cases, with a global annual incidence of 1.6-2.0/10 million [1]. In the Chinese population, CML accounts for 20% of all cases of leukemia and 95% of cases of chronic leukemia, with an annual incidence of 0.36/10 million [2].
Currently, the BCR-ABL gene mutation is considered the main molecular mechanism of CML. e translocation (9;22) q (34;11) entails the movement of the Abl protooncogene to the break-point cluster on the chromosome 22 long arm area (BCR) to form the BCR-ABL fusion gene [3]. e fusion gene encodes a protein with tyrosine kinase activity, leading to the expansion of the pool of myeloid progenitor cells. Myeloid progenitor cells display excessive growth as a hematologic disorder, which is difficult to treat without toxicity. Chemotherapy is the main treatment of CML, and the tyrosine kinase BCR-ABL is the main target. BCR-ABL has become the focus of drug research and development, of which significant advances have been made. Drugs such as imatinib, nilotinib, and dasatinib inhibit the phosphorylation activity of BCR-ABL and activation of downstream molecules, thus limiting the growth of CML cells [4]. However, the high costs make many medium-and low-income people unable to afford it. During chemotherapy, the targets of these drugs are associated with serious gastrointestinal toxicity, severe diarrhea, and blood toxicity.
In order to find new CML treatments, natural products in traditional Chinese medicine (TCM) have become attractive clinical candidates because of their long-proven clinical use with therapeutic effects in China [5]. In TCM theory, CML is viewed as an imbalance of many integral body processes that coordinate with each other, complement each other functionally, and influence each other pathologically. In long-term TCM clinical practice, a unique theory has been established for the treatment of hematological diseases, which is focused on the rebalance of the Qi and blood and of Yin and Yang. us, CML is believed to result from "poisoned bone marrow," which is observed as a loss of Qi or blood or a deficiency in Yin or Yang, resulting in a complex imbalance [6]. According to the principles of disease and syndrome classification, the general treatments and prescriptions for CML are divided into the following four dialectical therapy groups: "clearing heat and cooling blood," "relieving poisoning and detoxification," "eliminating wind and dispersing wet," and "supplementing Qi and nourishing Yin" [7]. e long history and multiple target effects of Chinese herbal medicine make the efficacy of TCM more durable and milder, with fewer by-effects. e use of Chinese medicine in the treatment of leukemia began as early as the 1960s. Paeonia anomala L., Rheum palmatum L., indigo, and others have been found to have efficacy for the treatment of leukemia [8][9][10].
e clinical application of TCM has been accepted for CML treatment in western China because of its fewer associated adverse reactions and low cost. It is widely used in local hospitals and has achieved good clinical outcomes [11,12]. TCM with active indigo ingredients has been successfully applied as an antileukemia drug. Although nearly hundreds of prescriptions for CML treatment have been recorded, the mechanisms of most TCM drugs remain unclear.
Recently, network pharmacology, genomics, proteomics, and metabolomics provide holistic approaches to study the essence of TCM and the functions of herbal compounds [13]. Meanwhile, systems biology-based analysis of TCM and drug-target interaction network approaches has made a significant contribution to study TCM [14]. is approach provides a possibility for a multitarget model study of TCM treatment for leukemia. Systems pharmacology studies of some TCM formulas have shown good results, such as that for the application of the Chinese herbal Radix Curcumae to treat cardiovascular disease, and these studies have provided insights into cardiovascular drug discovery and therapy [15].
e key molecules and mechanisms underlying CML treatment must be elucidated on a holistic level. us, we aimed at characterizing the active ingredients and target associations via both systems pharmacological and molecular biological methods, using an in silico approach combined with an in vitro approach to uncover the mechanism of TCM therapy in CML.

In Silico Screening for Potential Active Molecules via admetSAR Prediction.
To find the potential active molecules in CML-related TCM, all collected chemicals in the library were screened according to the predicted in silico oral bioavailability records via admetSAR (http://lmmd.ecust. edu.cn/admetsar1) [16]. e compounds with negative parameters of human intestinal absorption (HIA) and Caco-2 permeability (Caco-2) and negative metabolism parameters of inhibiting the five major CYP isoforms (CYP450 1A2, CYP450 2C9, CYP450 2C19, CYP450 2D6, and CYP450 3A4) were discarded, thereby remarkably reducing a great quantity of primitive chemicals to a smaller group, which was used to develop a TCM formula [17]. e inhibitory parameters of the five major CYP isoforms which contribute to the metabolism of drugs were calibrated comprehensively and calculated with the following equation [18]: Equation (1) [20]. To elucidate the effects of drugs, we collect compound targets to study compound-target interactions. e corresponding targets are obtained from TCMSP database searching by compound name and PubChem database on Document Mining. In addition, the predicted targets are further screened using the sysDT drug target prediction model [21]. To explain the relationship between candidate compounds and predicted targets and disease, we established "candidate compound-candidate target network (cC-cT network)," and "target-disease network" (T-D network), and all target networks were created by Cytoscape (http://www. cytoscape.org/). In the network, compounds and targets are represented by nodes, and interactions between two nodes are represented by edges [22]. e degree of a node is the number of edges connected to the node. e greater the degree is, the more nodes in the network that are directly related to the node, indicating that the node is more important in the network.

Cell
Culture. e human CML K562 cell line was obtained from Cancer Cell Repository (Shanghai Cell Bank, Shanghai, China). K562 cells were cultured in RPMI-1640 containing 10% fetal bovine serum and 100 U·mL − 1 each of penicillin and streptomycin in a humidified atmosphere maintained at 37°C and 5% CO 2 .

Analysis of Cell Viability.
e cells were inoculated with 4 × 10 4 /mL per well in a 96-well culture plate, and the volume of the culture medium was 200 μL. Different concentrations of α-viniferin were added and cells were cultured at 37°C for 12, 24, and 48 h. 20 μL of MTT solution (5 g/L) was added to each well, the cells were continued to culture for 4 h and centrifuged at 1000 r/min for 10 min, and the supernatant was discarded and 150 μL of DMSO was added. After shaking for 10 min, a plate reader (Varioskan Flash 3001, USA) detects the A value of each well at 490 nm and calculates the cell proliferation inhibition rate according to the following formula:

Analysis of Cell Apoptosis.
To detect morphological changes in nuclear chromatin of apoptotic cells using Hoechst 33258 and AO/EB staining, K562 cells treated with drugs or untreated were collected and washed with PBS. en, the cells were stained with either 10 mg·mL − 1 Hoechst 33258 or 100 mg·mL − 1 AO/EB double staining in PBS. Fluorescence microscopy (Axio Observer A1; Zeiss, Germany) was used to record apoptotic and necrotic cells.
To analyze apoptosis, K562 cells were stained with FITC Annexin V and PI and a flow cytometer was used (Calibur; BD Biosciences, Franklin Lakes, NJ, USA) to measure fluorescence.

Determination of Mitochondrial Membrane Potentials.
For JC-1 analysis, the K562 cells were treated with or without different α-viniferin concentrations for 24 h. We analyzed mitochondrial membrane potentials by JC-1. We analyzed samples by measuring the fluorescence value of flow cytometry (FACS Calibur; BD Biosciences).
Detection of Bcl-2 family gene expression was carried out using quantitative real-time RT-PCR (qRT-PCR), which was performed using a single-tube SYBR Green Kit (QIA-GEN, Valencia, CA, USA), a Rotor-Gene Q Real-Time PCR System (Rotor-Gene Q, QIAGEN), and specific primer sets. Bcl-2 family mRNA expression was calculated by the 2 − △△Ct method, and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an endogenous control.

Evaluation of Bcl-2 Family Protein Expression.
K562 cells were treated with or without different α-viniferin concentrations for 24 h, and then RIPA lysate was used to extract protein. e levels of Bcl-2 and Bax were analyzed using western blotting with the primary antibody (1 : 400) and its corresponding secondary antibody (1 : 10000) according to the manufacturer's instructions [25].

Statistical Analyses.
All data were used as the mean ± SD of at least three independent experiments and assessed by ANOVA. Student's t-test for multiple comparisons was used to identify differences between groups. P < 0.05 was considered to be statistically significant value.

Compounds Collecting and Sorting in CML Herbs.
In the Chinese ethnopharmacological system, CML treatments can be divided into four dialectical therapy groups according to syndrome differentiation, and they are highly connected to imbalances in the following four elements: Qi, blood, Yin, and Yang [26]. According to the four dialectical therapy groups, we collected prescriptions for treating CML (Supplemental Table 1), from the database of traditional Chinese medicine prescriptions (https://db.yaozh.com/fangji) and Chinese prescription database in National Scientific Data Sharing Platform for population and health (http:// cowork.cintcm.com/engine/outline?page�2&channelid�37595& ispage�yes). e databases contain the information of prescriptions contained in Treatise on Typhoid Fever, Synopsis of the Golden Chamber, Criteria of Syndrome and Treatment, Orthodox Surgery, and other books. e 11 most frequently used herbs in CML treatments are classified as follows: clearing heat and cooling blood ("Rheum palmatum L., Reynoutria japonica Houtt., Moutan officinalis (L.) Lindl. & Paxton, cortex, Paeonia anomala L."); relieving poisoning and detoxifying ("Smilax China L., Smilax glabra Roxb."); eliminating wind and dispersing wet ("Notopterygium incisum K.C.Ting ex H.T.Chang, Morus alba L., twig, dry and Eucommia ulmoides Oliv."); and nourishing Yin and invigorating Yang ("Stemona sessilifolia (Miq.) Miq. and Ginkgo biloba L., leaf, dry"). Recently, there has been increasing attention on validating and explaining the combined principles of TCM by modern approach. In those studies, a systems pharmacology method has been introduced and applied to investigate TCM, which may allow for the understanding of the combined effects of herbal formulas from the ADMET in silico, active chemical, target protein, and disease network perspectives [5]. We predicted the absorption and metabolism-associated properties of all compounds via admetSAR.  Table 2).
irty-two candidate compounds also met the positive metabolic properties and pharmacokinetic data ( Table 1).

Candidate Compounds-Candidate Targets (cC-cT)
Networks. Systems pharmacology can also provide new approaches for drug discovery for complex diseases. By considering drug actions in the context of the regulatory networks within which the drug targets and disease gene products function, network analysis has the potential to greatly increase our knowledge of the mechanisms underlying the multiple actions of drugs [27]. Here, we used a systems pharmacology method to determine the active components of TCM in the treatment of CML. e cC-cT network has 259 nodes including 118 candidate compounds and 141 candidate targets (Figure 1). e targets in the outer circle exhibited considerably fewer interactions with the candidate compounds than those in the inner circles (the green circles indicate a high degree, whereas the blue circles indicate a low degree). In this result, there were many candidate targets connected to one or two candidate compounds; however, some candidate targets could be affected by multiple compounds rather than a single compound. In the inner circles, the number of targets (degree ≥ 10) was 46 (32% of the total number of targets). e other targets (degree < 10, 68% of the total number of targets) were likely less related to leukemia. e results indicate that, under the network pharmacology research paradigm, candidate compounds and their candidate targets exhibit a relationship in which multiple compounds are acting on multiple targets [5].

T-D Network Analysis.
In this section, target-related diseases were searched from GeneCards (http://www. genecards.org) and TTD websites (http://bidd.nus.edu.sg/ group/TTD/ttd.asp) and used to construct a T-D network ( Figure 2). According to ICD-10 records in Medical Subject Headings (MeSH, http://www.nlm.nih.gov), 139 diseases were classified into 24 groups. For instance, acute myelogenous leukemia and chronic myelogenous leukemia belong to both Neoplasm Diseases (C04) and Hemic and Lymphatic Diseases (C15). To identify the CML diseaseassociated targets, 1526 targets were collected from TCMSP, DisGeNET, and OMIM. After mapping the disease targetscandidate targets network, we identified 50 targets related to CML (Supplemental Table 3), and twelve of them belong to 46 potential targets in the inner circles in Figure 1.
We selected 121 chemicals with positive absorption as a prerequisite for analysis. Among them, proapoptotic activities on K562 cells were reported previously. For example, diosmetin induces OCI-AML2 and K562 cell apoptosis via the caspase pathway [28]. Resveratrol induces apoptosis of K562 through p38 and JNK-regulated H2AX phosphorylation [29]. Kaempferol also showed a potent cytotoxic effect on K562 cells and U937 cells through Bcl-2 signaling pathway-induced apoptosis [30]. en, we constructed a cC-cT network and found that there were 46 targets with a close relationship to leukemia. From the chemical structure perspective, flavonoid (11 compounds) and stilbene (7 compounds) family were displayed frequently (Supplemental Table 4). Additionally, literature mining displayed inhibition data of those flavonoids on leukemic cell proliferation, but there have been few reports on stilbenes besides resveratrol. 11 candidate stilbenes from 121 compounds in the whole network of Figure 1 were predicted to have candidate protein targets ( Figure 3). We found that α-viniferin, a basic resveratrol trimer, was predicted to be closely related to the main targets of leukemia. is finding agrees with the report that resveratrol has antileukemic activities in vitro [31]. Next, experiments were performed to clarify the potential antileukemic activity of α-viniferin.

cC-cT Networks for Stilbene Compounds.
is cC-cT network included 11 stilbenes and their candidate targets, which consists of 152 nodes (11 candidate compounds and 141 candidate targets). . ese results indicate that resveratrol and α-viniferin were common stilbene compounds. In limited studies, stilbenes were proved to be highly related to targets of leukemia treatments, for example, pterostilbene, a 3′, 5′-dimethoxy-resveratrol interacted with 25 targets, including MAPK14, CDK2, and HSP90AA1, which are highly related to leukemia. Hsiao et al. have reported that pterostilbene induces apoptosis in HL-60 cells by MAPK-mediated regulation of the mitochondrial apoptotic pathway while also causing cell cycle arrest in the G0/G1 phase and inhibiting CDK2 expression [32]. In the present study, resveratrol and α-viniferin were found to be associated with antitumor targets, including Bcl-2, caspase-3, 8, and 9, SIRT1, and others. Taken together, these results show that α-viniferin, a resveratrol dimer, is predicted to be highly related to targets of leukemia, which are closely associated with its structure. us, the mechanism of leukemia treatment using α-viniferin needs to be clarified.

α-Viniferin Inhibits Cell Proliferation in K562 Cells.
Stilbene compounds were found to be capable of having high sensitivity towards leukemia cells by network analysis.   [35], 50 μg·mL − 1 (A549) [33], and 37 to 45 μg·mL − 1 (WRL-68) [36]. Why the low concentration of the stilbene compounds had a high sensitivity to leukemia cells? Previously, it has been reported that resveratrol can cause cell apoptosis in various kinds of cancer cell lines, such as leukemia cells [37], breast carcinoma cells [38], and hepatoma cells [39]. According to the National Cancer Institute Developmental erapeutics Program records (NSC 655524) (https://dtp.cancer.gov/ dtpstandard/dwindex/index.jsp), α-viniferin is more sensitive to K562 cells in 32 candidate compounds by comparing GI 50 and one dose growth percent (Supplemental Table 5). Resveratrol had remarkable cytotoxic effects and induced apoptosis in K562 cells, and IC 50 value is 9 μg·mL − 1 . ere were few records about α-viniferin inhibiting K562 cell growth. To verify the underlying antitumor activity of α-viniferin, MTT assay was used to detect the cytotoxic activity of different concentrations of α-viniferin ranging from 0.5 to 128 μg·mL − 1 on K562 cells. e results revealed that α-viniferin significantly inhibited the proliferation of K562 cells in both dose-and time-dependent manners ( Figure 4). e inhibition rates sharply increased at the concentration of α-viniferin from 8 to 128 μg·mL − 1 at 12, 24, and 48 h. We observed α-viniferin inhibits K562 cell growth IC 50 values of 13.61 μg·mL − 1 at 24 h. Compared with resveratrol, there were few records about α-viniferin, a resveratrol trimer, inhibiting K562 cell growth. But the effect of α-viniferin is not simply triple of resveratrol; it is probably due to another way of inhibiting the growth of K562 cells. We observed high-dose α-viniferin (32 to 128 μg·mL − 1 ) caused serious cell death, cell fragmentation, and nuclei lysis.  Evidence-Based Complementary and Alternative Medicine us, we chose a comparatively lower concentration range from 0 to 32 μg·mL − 1 at 24 h for the investigation of mechanisms relating to the inhibitory effects of α-viniferin on K562 cell proliferation. In order to verify the selectivity of α-viniferin on leukemia cells, peripheral blood mononuclear cells (PBMCs) from healthy C57BL/6 mice were isolated and incubated with Con A or different concentrations of α-viniferin for 24 h and 48 h to detect cell viability. In total, α-viniferin promoted PBMC proliferation, but the proliferation effect of PBMC was observed only in 32 μg·mL − 1 with significance (Supplemental Figure 1). We also noticed that α-viniferin inhibits the proliferation of K562 cells by reducing the expression of the BCR-ABL protein (Supplemental Figure 2).

α-Viniferin Induces Proapoptotic Morphology in K562
Cells. As shown in the microscopy images in Figure 5(a), following exposure to 0 to 128 μg·mL − 1 α-viniferin for 24 h, K562 cells exhibited apoptotic morphological features, including distortion and disruption, all of which showed no difference with the cell proliferation assay results. In Figures 5(b)-5(d), Hoechst-stained control cells showed a circular nucleus containing diffuse Hoechst chromatin staining. After α-viniferin treatment, the cell nuclei became granular, which is consistent with the results of apoptosis-positive drug H 2 O 2 . In addition, cells stained with AO/EB appeared as bright green arcs and had condensed yellow/orange nuclei during the early stage and the late stage of apoptosis; however, no necrotic cells were observed (Figures 5(e)-5(g)). Besides, the ultrastructural analysis further proved the absence of apoptotic cells in the control group ( Figure 5(h)). But apoptotic bodies were spotted in K562 cells treated with 32 μg·mL − 1 α-viniferin for 24 h (Figures 5(i) and 5(j)).  Evidence-Based Complementary and Alternative Medicine 7 staining analysis to evaluate apoptotic cell death. Despite the pronounced concentration-dependent cell death that was observed in the cytotoxicity assay, apoptotic cell death was pronounced only at a high concentration of α-viniferin ( Figure 6). To gain insights into the antileukemic effect of α-viniferin, we performed FITC Annexin V and PI staining of K562 cells. Cells treated for 24 h with α-viniferin revealed marked flow cytometry findings, including a few cells appearing in the second quadrant (FITC Annexin V-negative and PI-positive), indicating that necrosis rarely occurred in the drug-treated cells. e finding of drug-induced cell apoptosis, indicated by cells moving from the third (FITC Annexin V-and PI-negative) to the first quadrant (FITC Annexin V-and PI-positive) and fourth quadrants (FITC Annexin V-positive and PI-negative), suggested that apoptosis was induced by the treatment.

Mechanism of α-Viniferin-Induced Apoptosis.
To verify the potential mechanism by which α-viniferin induces apoptosis, we detected changes in the mitochondrial membrane potential involved in apoptosis. To explore the effects of α-viniferin on mitochondrial membrane potential, JC-1stained K562 cells were analyzed by flow cytometry, which were cultured in the presence of different drug concentrations for 24 h (Figure 7). Compared with untreated cells, mitochondrial membrane potential decreases in a dosedependent manner. is result showed that α-viniferin disrupts the mitochondrial membrane potential leading to cytosolic accumulation of monomeric JC-1, activating the intrinsic pathway to indicate apoptosis. e mRNA levels of Bad, Bax, Bcl-2, and Bcl-xL were measured in α-viniferin-induced and control K562 cells by quantitative real-time RT-PCR (qRT-PCR). As shown in Figure 8(a), when α-viniferin was added, the mRNA levels of apoptosis-promoting genes were upregulated with increasing drug concentrations. In contrast, antiapoptosis mRNA  Evidence-Based Complementary and Alternative Medicine expression levels were significantly downregulated. However, the results of Bcl-2 and Bax are not very significant compared to Bcl-xL and Bad. To further verify the observed changes in Bcl-2 family underlying α-viniferin-induced apoptosis, the protein levels of Bcl-2 and Bax were measured by western blotting. Compared with the control group, Bcl-2 expression decreased, and Bax expression increased after treated with α-viniferin for 24 h (Figures 8(b)-8(d)). And the mRNA levels of caspase-9 were upregulated in a dose-dependent manner measured by RT-PCR (Supplemental Figure 3).

Discussion
We found changes in the apoptosis pathway after drug treatment, and the expression of Bcl-2 family changed significantly.
ese findings suggest that stilbene compounds do not rely on toxicity to kill cancer cells; instead, cancer cells are killed through mitochondrial pathway interference of their growth and induction of apoptosis. Scientists working in the field of apoptosis agree that programmed cell death occurs via mitochondrial damage. For example, a recent study has found that resveratrol is a tumor suppressor compound in grapes that induces apoptosis through a mitochondrial pathway regulated by Bcl-2 [40]. Studies have shown that resveratrol induces apoptosis in different acute lymphoblastic leukemia cells by depolarizing mitochondrial membranes and activating caspase-9 [41]. Wieder et al. reported that a hydroxylation analogue of resveratrol, named piceatannol, could induce apoptosis in the lymphoma cell line BJAB and primary leukemic lymphoblasts [42]. Resveratrol-mediated stimulation of C7H2 cells has been shown to lead to the production of reactive oxygen species (ROS), which is remarkably reduced by Bcl-2 [32]. Another study has shown that resveratrol induces apoptosis in K562 cells by reducing mitochondrial membrane potential [43].
is study observed the same phenomenon for 32 μg·mL − 1 α-viniferin treated with K562 cells for 24 h.
Lower cytotoxicity and higher protective potency appear to be common characteristics of natural stilbenes. However, α-viniferin may play a role in plant disease resistance, human health, and other diseases. A large number of Chinese medicines analyzed by systems pharmacology have been found to contain α-viniferin and to exhibit an effect on leukemia cells. According to the National Cancer Institute Developmental erapeutics Program records (NSC 655524) (https://dtp.cancer.gov/dtpstandard/dwindex/index.jsp), leukemia and central nervous system cell lines are sensitive to α-viniferin treatments in vitro. We used the Comparative Toxicogenomics Database to evaluate the specificity of α-viniferin against different diseases according to inferred scores [44]. Among the top 100 inferred diseases related to α-viniferin, cancer and cardiovascular disease were highly correlated with α-viniferin. ese results indicate the possible proapoptotic activity of α-viniferin for potential application in treating cancer, whose effects were observed in our work. However, the antitumor mechanism of α-viniferin is likely greater than this single effect. Its role in cell cycle arrest, autophagy, and other functions must be investigated further.

Conclusions
In summary, our study is the first systems pharmacology investigation of the treatment of CML by TCM. We generated 118 bioactive compounds from 11 herbs within four dialectical therapy groups and 141 predicted targets. Our data showed that α-viniferin, a natural stilbene, was found to be highly related to these targets. And it has lower cytotoxicity and high sensitivity to K562 cells. We found α-viniferin could induce early apoptosis on K562 cells, especially considering the Bcl-2 family of proteins in regulating K562 cell apoptosis through unique modulation of the mitochondrial pathway. α-Viniferin caused mitochondrial damage and increased Bcl-2 family proapoptotic protein activity, leading to apoptosis. us, we hypothesize that α-viniferin merits additional experimental and clinical research in the treatment of CML.

Supplementary Materials
Supplemental Table 1: TCM prescriptions for CML. Supplemental Table 2: full ADME properties in silico of 121 candidate compounds via admetSAR Prediction. Supplemental Table 3: intersection of CML targets and prediction Targets. Supplemental Table 4: literature mining of 32 active compounds. Supplemental Table 5: NCI Cancer Screen Data of 32 candidate compounds on K562 cells. Supplemental Figure 1: effects of α-viniferin on peripheral blood mononuclear cells (PBMC) proliferation. Supplemental Figure 2: evaluation of α-viniferin on the protein expression level of BCR-ABL. Supplemental Figure 3: evaluation of α-viniferin on the mRNA level of caspase-9. (Supplementary Materials)