key: cord-003539-tazd6dvm authors: Yang, Kun; Jin, Ming-Ji; Quan, Zhe-Shan; Piao, Hu-Ri title: Design and Synthesis of Novel Anti-Proliferative Emodin Derivatives and Studies on their Cell Cycle Arrest, Apoptosis Pathway and Migration date: 2019-03-02 journal: Molecules DOI: 10.3390/molecules24050884 sha: doc_id: 3539 cord_uid: tazd6dvm Emodin is a cell arrest and apoptosis-inducing compound that is widely distributed in different plants (rhubarb, aloe), lichens and terrestrial fungi, and also isolated from marine-derived fungi and marine sponge-associated fungi. In this study, we designed and synthesized a novel series of emodin derivatives by binding emodin to an amino acid using linkers of varying lengths and composition, and evaluated their anti-proliferative activities using HepG2 cells (human hepatic carcinoma), MCF-7 cells (human breast cancer) and human normal liver L02 cells. Most of these derivatives showed moderate to potent anti-proliferative activities. Notably, compound 7a exhibited potent anti-proliferative activity against HepG2 cells with the half maximal inhibitory concentration (IC(50)) value of 4.95 µM, which was enhanced 8.8-fold compared to the parent compound emodin (IC(50) = 43.87 µM), and it also exhibited better selective anti-proliferative activity and specificity than emodin. Moreover, further experiments demonstrated that compound 7a displayed a significant efficacy of inducing apoptosis through mitochondrial pathway via release of cytochrome c from mitochondria and subsequent activation of caspase-9 and caspase-3, inducing cell arrest at G0/G1 phase, as well as suppression of cell migration of tumor cells. The preliminary results suggested that compound 7a could be a promising lead compound for the discovery of novel anti-tumor drugs and has the potential for further investigations as an anti-cancer drug. Emodin (3-methyl-1,6,8-trihydroxyanthraquinone, Figure 1 ) is a naturally occurring anthraquinone derivative that is widely distributed in different plants (rhubarb, aloe), lichens and terristrial fungi, that has also been isolated from marine-derived fungi and marine sponge-associated fungi [1] [2] [3] [4] . It is also an active ingredient in Chinese medicinal herbs used in the treatment of constipation, jaundice, gastro-intestinal hemorrhage, and ulcers. Emodin has the same tricyclic planar chromophore skeleton as certain anti-tumor antibiotics such as daunorubicin, doxorubicin, epirubicin, and several synthetic analogs distributed such as mitoxantrone and pixantrone ( Figure 1) ; all of which are in used for clinical treatment of various cancers. Many reports have demonstrated that emodin possesses a wide spectrum of pharmacological effects such as anti-tumour [5] [6] [7] , anti-inflammatory [8, 9] , antiviral [10] , antibacterial [11] , anti-allergic [12] , anti-osteoporotic [13] , anti-diabetic [14, 15] , immunosuppressive [16] , neuroprotective [17] , and hepatoprotective [18, 19] activities, etc. Among these effects, anti-tumor anti-osteoporotic [13] , anti-diabetic [14, 15] , immunosuppressive [16] , neuroprotective [17] , and hepatoprotective [18, 19] activities, etc. Among these effects, anti-tumor activity is the most widely reported. In recent decades, pharmacological studies have shown that emodin is capable of inhibiting cellular proliferation [20] , inducing of cell differentiation [21] and apoptosis [22] , and activating of caspase cascade pathway [23, 24] and mitochondrial death pathway [25, 26] in different cancer cells. Chemical transformation of bioactive compounds of natural products is one of the most common approaches in drug discovery to improve therapeutic properties. To date, the anti-tumor activity of emodin has been improved by structural modification, mainly to the side chain including methyl, hydroxyl, and/or aryl ring groups. For example, Wang et al. reported that emodin quaternary ammonium salt derivatives showed significant anti-cancer activities against hepatoma cells [27] [28] [29] . Xing et al. reported the effects of emodin rhamnoside derivatives against human cancer cells [30] , and Tan et al. reported DNA-binding pyrazole emodin derivatives [31] . However, emodin is not an ideal chemotherapeutic agent for cancer due to its poor bioavailability, low solubility, and high toxicity in vivo. The bioavailability of a drug is positively correlated with its solubility; thus amino acids, which possess the carboxylic and amino functionality, are ideal as a moiety for the structural modification of bioactive natural products. Amino acids can interact with other biomolecules via secondary interactions such as hydrogen bonding to improve their pharmacological profiles in both potency and bioavailability [32] . For example, conjugation of an amino acid to betulinic acid is able to improve its water-solubility as well as its anti-melanoma activity [33] . These findings prompted us to design and synthesize a series of novel emodin derivatives linked with various amino acids as TFA salt ( Figure 2 ). Among these derivatives, compound 3b derived from (R)-2-aminopropanoic acid exhibited the most anti-proliferative activity against both HepG2 cells and MCF-7 cells. This prompted us to screen the different length of diol linker between 2-aminopropanoic acid and emodin. Our results also suggest that compared to compound 3a, compound 7a exhibited a slightly more potent Chemical transformation of bioactive compounds of natural products is one of the most common approaches in drug discovery to improve therapeutic properties. To date, the anti-tumor activity of emodin has been improved by structural modification, mainly to the side chain including methyl, hydroxyl, and/or aryl ring groups. For example, Wang et al. reported that emodin quaternary ammonium salt derivatives showed significant anti-cancer activities against hepatoma cells [27] [28] [29] . Xing et al. reported the effects of emodin rhamnoside derivatives against human cancer cells [30] , and Tan et al. reported DNA-binding pyrazole emodin derivatives [31] . However, emodin is not an ideal chemotherapeutic agent for cancer due to its poor bioavailability, low solubility, and high toxicity in vivo. The bioavailability of a drug is positively correlated with its solubility; thus amino acids, which possess the carboxylic and amino functionality, are ideal as a moiety for the structural modification of bioactive natural products. Amino acids can interact with other biomolecules via secondary interactions such as hydrogen bonding to improve their pharmacological profiles in both potency and bioavailability [32] . For example, conjugation of an amino acid to betulinic acid is able to improve its water-solubility as well as its anti-melanoma activity [33] . These findings prompted us to design and synthesize a series of novel emodin derivatives linked with various amino acids as TFA salt ( Figure 2 ). anti-osteoporotic [13] , anti-diabetic [14, 15] , immunosuppressive [16] , neuroprotective [17] , and hepatoprotective [18, 19] activities, etc. Among these effects, anti-tumor activity is the most widely reported. In recent decades, pharmacological studies have shown that emodin is capable of inhibiting cellular proliferation [20] , inducing of cell differentiation [21] and apoptosis [22] , and activating of caspase cascade pathway [23, 24] and mitochondrial death pathway [25, 26] in different cancer cells. Chemical transformation of bioactive compounds of natural products is one of the most common approaches in drug discovery to improve therapeutic properties. To date, the anti-tumor activity of emodin has been improved by structural modification, mainly to the side chain including methyl, hydroxyl, and/or aryl ring groups. For example, Wang et al. reported that emodin quaternary ammonium salt derivatives showed significant anti-cancer activities against hepatoma cells [27] [28] [29] . Xing et al. reported the effects of emodin rhamnoside derivatives against human cancer cells [30] , and Tan et al. reported DNA-binding pyrazole emodin derivatives [31] . However, emodin is not an ideal chemotherapeutic agent for cancer due to its poor bioavailability, low solubility, and high toxicity in vivo. The bioavailability of a drug is positively correlated with its solubility; thus amino acids, which possess the carboxylic and amino functionality, are ideal as a moiety for the structural modification of bioactive natural products. Amino acids can interact with other biomolecules via secondary interactions such as hydrogen bonding to improve their pharmacological profiles in both potency and bioavailability [32] . For example, conjugation of an amino acid to betulinic acid is able to improve its water-solubility as well as its anti-melanoma activity [33] . These findings prompted us to design and synthesize a series of novel emodin derivatives linked with various amino acids as TFA salt ( Figure 2 ). Among these derivatives, compound 3b derived from (R)-2-aminopropanoic acid exhibited the most anti-proliferative activity against both HepG2 cells and MCF-7 cells. This prompted us to screen the different length of diol linker between 2-aminopropanoic acid and emodin. Our results also suggest that compared to compound 3a, compound 7a exhibited a slightly more potent Among these derivatives, compound 3b derived from (R)-2-aminopropanoic acid exhibited the most anti-proliferative activity against both HepG2 cells and MCF-7 cells. This prompted us to screen the different length of diol linker between 2-aminopropanoic acid and emodin. Our results also suggest that compared to compound 3a, compound 7a exhibited a slightly more potent anti-proliferative activity against HepG2 cells and MCF-7 cells. Following this, the anti-cancer mechanism of compound 7a was evaluated against HepG2 cells. A novel series of emodin derivatives linked with amino acids were designed and synthesized. The synthesis strategies for these emodin derivatives are outlined in Schemes 1-3. For the synthesis of compounds 3a-3z, compound 1 was subjected to an alkylation reaction with 2-iodoethanol and Cs 2 CO 3 in DMF that resulted in compound 2. Then coupling reactions under the conventional coupling condition (DCC/DMAP) between compound 2 and various commercially available N-Boc protected amino acids resulted in the formation of condensation products. After removal of the protecting group usinng 20% TFA in DCM, the target compounds 3a-3z were obtained as TFA salts. Molecules 2019, 24, x FOR PEER REVIEW 3 of 23 anti-proliferative activity against HepG2 cells and MCF-7 cells. Following this, the anti-cancer mechanism of compound 7a was evaluated against HepG2 cells. A novel series of emodin derivatives linked with amino acids were designed and synthesized. The synthesis strategies for these emodin derivatives are outlined in Scheme 1-3. For the synthesis of compounds 3a-3z, compound 1 was subjected to an alkylation reaction with 2-iodoethanol and Cs2CO3 in DMF that resulted in compound 2. Then coupling reactions under the conventional coupling condition (DCC/DMAP) between compound 2 and various commercially available N-Boc protected amino acids resulted in the formation of condensation products. After removal of the protecting group usinng 20% TFA in DCM, the target compounds 3a-3z were obtained as TFA salts. Following the DCC/DMAP coupling condition and deprotection under 20% TFA in DCM for the synthesis of compounds 3a-3z, the target compounds 4a and 4b were obtained as TFA salts (Scheme 2). Scheme 1. Synthesis of compounds 3a-3z. Reagents and conditions: (a) 2-iodoethanol, Cs 2 CO 3 , DMF, 60 • C, 65%; (b) (i) various N-Boc amino acids, DCC, DCM, 0 • C; (ii) 20% TFA in DCM, r.t., 20%-50% over two steps. Following the DCC/DMAP coupling condition and deprotection under 20% TFA in DCM for the synthesis of compounds 3a-3z, the target compounds 4a and 4b were obtained as TFA salts (Scheme 2). These products were then coupled with formaldehyde (37% aqueous solution) and reduced utilizing NaBH3CN to provide 5a and 5b as TFA salt in 85% and 89% yields, respectively. As shown in Scheme 3, by employing the route for the synthesis of compounds 3a-3z, compounds 7a-7l were obtained as TFA salts. Scheme 3. Synthesis of compounds 7a-7l. Reagents and conditions: (a) various hydroxyalkyl bromides or iodides, Cs2CO3, DMF, 60 °C, 30%-55%; (b) (i) R/S-N-Boc-Ala-OH, DCC, DCM, 0 °C; (ii) 20% TFA in DCM, r.t., 65%-70% over two steps. The in vitro anti-proliferative activities of all the novel synthesized compounds were evaluated against HepG2 cells and MCF-7 cells by CellTiter-Glo ® Luminescent Cell Viability Assay, using paclitaxel (a clinically used drug) as a positive control. Emodin (1) and compound 2 were also included for comparison. As shown in Table 1 , emodin (1) showed weak inhibitory activity against HepG2 cells and MCF-7 cells with IC50 values of 43.87 ± 1.28 µM and 52.72 ± 2.22 µM, respectively. Compared to emodin (1), compound 2 showed reduced activity due to the introduction of the hydroxyethyl group at the 3-position of emodin. Among compounds 3a-3z containing amino acid groups, it is worth noting that most of the compounds displayed more potent activity than emodin (1) . Compound 3a with the Gly group displayed near double the activity of emodin (1) against HepG2 cells, but had slightly decreased activity against MCF-7 cells. Compound 3b with the D-Ala group displayed the most potent anti-proliferative activity against HepG2 cells and MCF-7 cells with IC50 values of 8.22 ± 0.13 µM and 8.99 ± 0.64 µM, respectively. Compound 3c with the L-Ala group displayed slightly weaker activity than compound 3b, but still displayed stronger activity than emodin (1). Compounds 3d-3j containing alkyl chain amino acid groups also displayed stronger inhibitory activities than emodin (1) . Overall, with the exception of compound 3n, compounds 3k-3m with cyclic alkyl amino acid groups displayed worse inhibitory activities than emodin (1); especially compound 3m containing a heterocycle. Interestingly, compound 3n which contains a cyclopentane amino acid group displayed higher anti-proliferative activity against HepG2 cells and MCF-7 cells with IC50 values of 12.48 ± 0.59 µM Scheme 2. Synthesis of compounds 4a, 4b, 5a and 5b. Reagents and conditions: (a) (i) Boc-N-Me-R/ S-Ala-OH, DCC, DCM, 0 • C; (ii) 20% TFA in DCM, r.t., 40%-50% over two steps; (b) formaldehyde (37% aqueous solution), NaBH 3 CN, MeOH, r.t., 85%-89%. These products were then coupled with formaldehyde (37% aqueous solution) and reduced utilizing NaBH 3 CN to provide 5a and 5b as TFA salt in 85% and 89% yields, respectively. As shown in Scheme 3, by employing the route for the synthesis of compounds 3a-3z, compounds 7a-7l were obtained as TFA salts. These products were then coupled with formaldehyde (37% aqueous solution) and reduced utilizing NaBH3CN to provide 5a and 5b as TFA salt in 85% and 89% yields, respectively. As shown in Scheme 3, by employing the route for the synthesis of compounds 3a-3z, compounds 7a-7l were obtained as TFA salts. The in vitro anti-proliferative activities of all the novel synthesized compounds were evaluated against HepG2 cells and MCF-7 cells by CellTiter-Glo ® Luminescent Cell Viability Assay, using paclitaxel (a clinically used drug) as a positive control. Emodin (1) and compound 2 were also included for comparison. As shown in Table 1 , emodin (1) showed weak inhibitory activity against HepG2 cells and MCF-7 cells with IC50 values of 43.87 ± 1.28 µM and 52.72 ± 2.22 µM, respectively. Compared to emodin (1), compound 2 showed reduced activity due to the introduction of the hydroxyethyl group at the 3-position of emodin. Among compounds 3a-3z containing amino acid groups, it is worth noting that most of the compounds displayed more potent activity than emodin (1) . Compound 3a with the Gly group displayed near double the activity of emodin (1) against HepG2 cells, but had slightly decreased activity against MCF-7 cells. Compound 3b with the D-Ala group displayed the most potent anti-proliferative activity against HepG2 cells and MCF-7 cells with IC50 values of 8.22 ± 0.13 µM and 8.99 ± 0.64 µM, respectively. Compound 3c with the L-Ala group displayed slightly weaker activity than compound 3b, but still displayed stronger activity than emodin (1). Compounds 3d-3j containing alkyl chain amino acid groups also displayed stronger inhibitory activities than emodin (1) . Overall, with the exception of compound 3n, compounds 3k-3m with cyclic alkyl amino acid groups displayed worse inhibitory activities than emodin (1); especially compound 3m containing a heterocycle. Interestingly, compound 3n which contains a cyclopentane amino acid group displayed higher anti-proliferative activity against HepG2 cells and MCF-7 cells with IC50 values of 12.48 ± 0.59 µM The in vitro anti-proliferative activities of all the novel synthesized compounds were evaluated against HepG2 cells and MCF-7 cells by CellTiter-Glo ® Luminescent Cell Viability Assay, using paclitaxel (a clinically used drug) as a positive control. Emodin (1) and compound 2 were also included for comparison. As shown in Table 1 , emodin (1) showed weak inhibitory activity against HepG2 cells and MCF-7 cells with IC 50 values of 43.87 ± 1.28 µM and 52.72 ± 2.22 µM, respectively. Compared to emodin (1), compound 2 showed reduced activity due to the introduction of the hydroxyethyl group at the 3-position of emodin. Among compounds 3a-3z containing amino acid groups, it is worth noting that most of the compounds displayed more potent activity than emodin (1) . Compound 3a with the Gly group displayed near double the activity of emodin (1) against HepG2 cells, but had slightly decreased activity against MCF-7 cells. Compound 3b with the D-Ala group displayed the most potent anti-proliferative activity against HepG2 cells and MCF-7 cells with IC 50 values of 8.22 ± 0.13 µM and 8.99 ± 0.64 µM, respectively. Compound 3c with the L-Ala group displayed slightly weaker activity than compound 3b, but still displayed stronger activity than emodin (1). Compounds 3d-3j containing alkyl chain amino acid groups also displayed stronger inhibitory activities than emodin (1) . Overall, with the exception of compound 3n, compounds 3k-3m with cyclic alkyl amino acid groups displayed worse inhibitory activities than emodin (1); especially compound 3m containing a heterocycle. Interestingly, compound 3n which contains a cyclopentane amino acid group displayed higher anti-proliferative activity against HepG2 cells and MCF-7 cells with IC 50 values of 12.48 ± 0.59 µM and 23.51 ± 1.93 µM, respectively. Compounds 3o and 3p with D-and L-Ala groups exhibited similar inhibitory activity to compound 3n. Compared to emodin (1), compounds 3q-3z with aryl (phenyl, substituted phenyl, benzyl, substituted benzyl, etc.) amino acid groups displayed enhanced activity in varying degrees when tested against two cancer cell lines, with IC 50 values ranging from 11.41 ± 0.77 to 42.44 ± 1.16 µM. Following our initial results, we performed further modification of compounds 3b and 3c to improve the anti-proliferative activity, thus mono-and dimethylation of the amino moiety in the amino acid afforded compounds 4a, 4b, 5a and 5b. Unfortunately, this modification led to decreased activity to varying degrees against both the HepG2 and MCF-7 cell lines, with IC 50 values ranging from 13.69 ± 0.89 to 68.81 ± 2.25 µM, as shown in Table 2 . Notably, compound 4b displayed a 6-fold activity reduction against MCF-7 cell lines compared to emodin (1). The results revealed that "NH 2 " was the moiety that favors improving the anti-proliferative activity. The preliminary explanation was that amino group can interact with other biomolecules to improve the anti-proliferative activity via formation of hydrogen bond. Mono-or dimethylation of amino group might weaken the trend of hydrogen bonding and decrease the anti-proliferative activity. Next, as shown in Table 3 , the linking group of compounds 3b and 3c was further optimized by increasing the linker length in compounds 3b and 3c by one additional methylene group. The resulting compounds, 7a and 7b, displayed further enhanced potency at inhibiting cell growth against both HepG2 cells and MCF-7 cells. However, compounds 7c-7f initially failed to increase the potency of cell growth inhibition against both tested cancer cell lines; this was successfully addressed by increasing the linker length using additional methylene bridge. On the other hand, elongation of the linker by two to four ethanediol groups resulted in compounds 7g-7l, which showed decreased activity compared to compounds 3b and 3c. Three pairs of compounds possessing stronger anti-proliferative activities were screened out and their cytotoxic activities were preformed against human normal liver cells (L02) in vitro. As shown in Table 4 , the parent compound emodin exhibited stronger cytotoxic activity against L02 cells with an IC 50 value of 22.52 ± 0.18 µM than against cancer cell lines (IC 50 = 43.87 ± 1.28 µM against HepG2 cells and IC 50 = 52.72 ± 2.22 µM against MCF-7 cells). However, all the six compounds exhibited weaker cytotoxic activities against L02 cells compared with their corresponding cancer cell lines. Above all, we finally chose to investigate the mechanism of action for compound 7a because it had the strongest anti-proliferative activity against HepG2 cells. Lack of selective cytotoxicity is the main factor that hinders conventional chemotherapeutic agents. Thus, to evaluate the selective anti-proliferative activity of the compound 7a, the selectivity index (SI) between cancer and normal cells was calculated and the results are summarized in Table 5 . The SI was calculated by dividing the IC 50 values in normal cells by the IC 50 values in cancer cells. Emodin displayed stronger cytotoxic activity, with SI values of 0.51 (HepG2) and 0.43 (MCF-7), respectively, indicating that both normal cells and cancer cells would be killed. However, compound 7a exhibited weaker cytotoxic activity with SI value of 2.82 (HepG2) and 2.36 (MCF-7), respectively. It means that 7a exhibited a better selective anti-proliferative activity and specificity than emodin. In order to verify whether compound 7a is able to induce apoptosis in HepG2 cells, we utilized FITC-Annexin V/PI staining and estimated the percentage of apoptotic cells by flow cytometry. We noted a concentration-dependent increase in the percentage of apoptotic cells when the cells were treated with compound 7a for 48 h at concentrations 2.5, 5, and 10 µM. As shown in Figure 3A , few (5.5%) apoptotic cells were present in the control panel. In contrast, the percentage of apoptotic cells increased to 22.2% after treatment with compound 7a at 5 µM for 48 h and further increased to 50.7% after treatment with 7a at the concentration of 10 µM. As illustrated in Figure 3B , the quantitative analysis of apoptosis strongly suggests that treatment with compound 7a effectively induced apoptosis in HepG2 cells in a concentration-dependent manner in comparison to the control. In order to verify whether compound 7a is able to induce apoptosis in HepG2 cells, we utilized FITC-Annexin V/PI staining and estimated the percentage of apoptotic cells by flow cytometry. We noted a concentration-dependent increase in the percentage of apoptotic cells when the cells were treated with compound 7a for 48 h at concentrations 2.5, 5, and 10 µM. As shown in Figure 3A , few (5.5%) apoptotic cells were present in the control panel. In contrast, the percentage of apoptotic cells increased to 22.2% after treatment with compound 7a at 5 µM for 48 h and further increased to 50.7% after treatment with 7a at the concentration of 10 µM. As illustrated in Figure 3B , the quantitative analysis of apoptosis strongly suggests that treatment with compound 7a effectively induced apoptosis in HepG2 cells in a concentration-dependent manner in comparison to the control. Statistical significance is determined by two-tailed Student t-test: "***" denote p < 0.001, "**" denote p < 0.01, respectively (Supplementary Table S1 ). (C, E) Western blot analysis effect of compound 7a on the levels of Bax, Bcl-2, cytochrome c, procaspase-3, caspase-3 and procaspase-9 expression in HepG2 cells. (D, F) An equal amount of protein was loaded on SDS-PAGE gel for western blot analysis. Data are expressed as means ± SD of the percentages of apoptotic cells from three independent experiments. Statistical significance is determined by two-tailed Student t-test: "***" denote p < 0.001, "**" denote p < 0.01, "*" denote p < 0.05, respectively (Supplementary Table S2 ). To verify the molecular mechanisms of apoptosis induction of compound 7a, we performed a western blot assay. It is well known that the Bcl-2 family of pro-apoptotic and anti-apoptotic The quantitative analysis of apoptosis. Data are expressed as means ± SD of the percentages of apoptotic cells from three independent experiments. Statistical significance is determined by two-tailed Student t-test: "***" denote p < 0.001, "**" denote p < 0.01, respectively (Supplementary Table S1 ). (C, E) Western blot analysis effect of compound 7a on the levels of Bax, Bcl-2, cytochrome c, procaspase-3, caspase-3 and procaspase-9 expression in HepG2 cells. (D, F) An equal amount of protein was loaded on SDS-PAGE gel for western blot analysis. Data are expressed as means ± SD of the percentages of apoptotic cells from three independent experiments. Statistical significance is determined by two-tailed Student t-test: "***" denote p < 0.001, "**" denote p < 0.01, "*" denote p < 0.05, respectively (Supplementary Table S2 ). To verify the molecular mechanisms of apoptosis induction of compound 7a, we performed a western blot assay. It is well known that the Bcl-2 family of pro-apoptotic and anti-apoptotic proteins regulates the mitochondrial pathway of apoptosis. These Bcl-2 family proteins stimulate the permeabilization of the mitochondrial outer membrane, which results in the release of cytochrome c into the cytosol and in turn promotes the activation of the caspase cascade. The activation of the caspase cascade ultimately leads to the induction of apoptotic cell death. As shown in Figure 3C and 3D, in comparison with the control cells, compound 7a induced an increase in the levels of Bax and a decrease in the expression of Bcl-2 in a concentration-dependent manner. Meanwhile, the release of cytochrome c from mitochondria increased after the treatment of compound 7a, while procaspase 9 and procaspase 3 decreased after treatment with 7a, indicating that the caspase 9 and caspase 3 were activated. As shown in Figure 3E ,F, the increased expression of cleaved caspase-3 after treatment with 7a provided a further evidence that compound 7a induced cell apoptosis through mitochondrial pathway in a concentration-dependent manner. The apoptosis process can be summarized as follows: The mitochondrial apoptosis-induced channel (MAC) of HepG2 cells was formed by pro-apoptotic protein Bax after the treatment of compound 7a. The formation of MAC led to the releasing of cytochrome c from mitochondria. Once cytochrome c was released, it binded with apoptotic protease activating factor-1 (Apaf-1) and ATP, which then binded to procaspase-9 to create a protein complex known as apoptosome. The apoptosome cleaved the pro-caspase-9 to its active form of initiator caspase-9, which in turn activated procaspase-3 and then the effector caspase-3 and finally resulted in cell apoptosis. To further examine how compound 7a suppressed the growth of HepG2 cells, the effect of compound 7a on cell cycle distribution with different concentrations was investigated by flow cytometric analysis following staining the DNA with propidium iodide (PI). The results of a typical experiment are shown in Figure 4A . As determined by flow cytometry, the exposure of HepG2 cells to compound 7a for 48 h resulted in an obvious increase in the percentage of cells in G0/G1 phase in comparison with the control. Treatment with compound 7a resulted in an increase of G1 phase cells (63.37%) compared to the control (33.16%). Inversely, S phase cell population decreased to 14.60% compared to the control (48.13%). Therefore, compound 7a resulted in a significant G0/G1 phase arrest in a concentration-dependent manner with a concomitant decrease in the number of cells in the S phase of the cycle. proteins regulates the mitochondrial pathway of apoptosis. These Bcl-2 family proteins stimulate the permeabilization of the mitochondrial outer membrane, which results in the release of cytochrome c into the cytosol and in turn promotes the activation of the caspase cascade. The activation of the caspase cascade ultimately leads to the induction of apoptotic cell death. As shown in Figure 3C and 3D, in comparison with the control cells, compound 7a induced an increase in the levels of Bax and a decrease in the expression of Bcl-2 in a concentration-dependent manner. Meanwhile, the release of cytochrome c from mitochondria increased after the treatment of compound 7a, while procaspase 9 and procaspase 3 decreased after treatment with 7a, indicating that the caspase 9 and caspase 3 were activated. As shown in Figure 3E ,F, the increased expression of cleaved caspase-3 after treatment with 7a provided a further evidence that compound 7a induced cell apoptosis through mitochondrial pathway in a concentration-dependent manner. The apoptosis process can be summarized as follows: The mitochondrial apoptosis-induced channel (MAC) of HepG2 cells was formed by pro-apoptotic protein Bax after the treatment of compound 7a. The formation of MAC led to the releasing of cytochrome c from mitochondria. Once cytochrome c was released, it binded with apoptotic protease activating factor-1 (Apaf-1) and ATP, which then binded to procaspase-9 to create a protein complex known as apoptosome. The apoptosome cleaved the pro-caspase-9 to its active form of initiator caspase-9, which in turn activated procaspase-3 and then the effector caspase-3 and finally resulted in cell apoptosis. To further examine how compound 7a suppressed the growth of HepG2 cells, the effect of compound 7a on cell cycle distribution with different concentrations was investigated by flow cytometric analysis following staining the DNA with propidium iodide (PI). The results of a typical experiment are shown in Figure 4A . As determined by flow cytometry, the exposure of HepG2 cells to compound 7a for 48 h resulted in an obvious increase in the percentage of cells in G0/G1 phase in comparison with the control. Treatment with compound 7a resulted in an increase of G1 phase cells (63.37%) compared to the control (33.16%). Inversely, S phase cell population decreased to 14.60% compared to the control (48.13%). Therefore, compound 7a resulted in a significant G0/G1 phase arrest in a concentration-dependent manner with a concomitant decrease in the number of cells in the S phase of the cycle. Table S3 ). To evaluate the effect of compound 7a on cancer migration, a wound healing assay was conducted to determine whether compound 7a could prevent HepG2 cell migration. After culturing HepG2 cells for 48 h in the presence and absence of compound 7a at 0, 0.25, 0.5, and 1 µM, a pipette tip was streaked through the cell culture, resulting in a cell deficient space. The territory recovered by the HepG2 cells was used to observe the migration inhibition capability of compound 7a. The empty space reduced significantly in size in the absence of compound 7a because of cell proliferation and migration and results indicate compound 7a suppressed the mobility of tumor cells in a concentration-dependent manner ( Figure 5 ). cell cycle distribution after treatment of compound 7a. Data are expressed as means ± SD of the percentages of apoptotic cells from three independent experiments. Statistical significance is determined by two-tailed Student t-test: "***" denote p < 0.001, "*" denote p < 0.05 and ns means no significance, respectively (Supplementary Table S3 ). To evaluate the effect of compound 7a on cancer migration, a wound healing assay was conducted to determine whether compound 7a could prevent HepG2 cell migration. After culturing HepG2 cells for 48 h in the presence and absence of compound 7a at 0, 0.25, 0.5, and 1 µM, a pipette tip was streaked through the cell culture, resulting in a cell deficient space. The territory recovered by the HepG2 cells was used to observe the migration inhibition capability of compound 7a. The empty space reduced significantly in size in the absence of compound 7a because of cell proliferation and migration and results indicate compound 7a suppressed the mobility of tumor cells in a concentration-dependent manner ( Figure 5 ). Table S4 ). The starting materials and reagents, purchased from commercial suppliers, were used without further purification. Emodin extracted from Polygonum cuspidatum was purchased from China Xi'an Sino-Herb Bio-technology Co., Ltd. (Xi'an, China). All reactions were monitored by thin-layer chromatography (TLC) on aluminum sheets (Silica gel 60-F254, E. Merck, Darmstadt, Germany). Compounds were visualized by UV light. Column chromatography was carried out using silica gel (200-300 mesh). All reaction solvents were dried prior to use according to standard procedures. All primary reagents were commercially available. Silica gel chromatography solvents were of analytical grade. NMR spectra were recorded in DMSO-d6 on a Bruker-250 spectrometer (Bruker Biospin, Fällanden, Switzerlahd), at 400 MHz for 1 H-NMR, 101 MHz for 13 C-NMR and 376 MHz for 19 F-NMR with TMS as the internal standard. Chemical shifts were expressed in δ (ppm) and coupling constants (J) in Hz. Multiplicity was indicated as follows: s (singlet), d (doublet), t (triplet), p (quintet), dd (doublet of doublets), brs (broad singlet), etc. Mass spectra were obtained on an Agilent 1100 Series LC/MSD Trap mass spectrometer (ESI-MS, Agilent, Santa Clara, CA, USA). Table S4 ). The starting materials and reagents, purchased from commercial suppliers, were used without further purification. Emodin extracted from Polygonum cuspidatum was purchased from China Xi'an Sino-Herb Bio-technology Co., Ltd. (Xi'an, China). All reactions were monitored by thin-layer chromatography (TLC) on aluminum sheets (Silica gel 60-F254, E. Merck, Darmstadt, Germany). Compounds were visualized by UV light. Column chromatography was carried out using silica gel (200-300 mesh). All reaction solvents were dried prior to use according to standard procedures. All primary reagents were commercially available. Silica gel chromatography solvents were of analytical grade. NMR spectra were recorded in DMSO-d 6 on a Bruker-250 spectrometer (Bruker Biospin, Fällanden, Switzerlahd), at 400 MHz for 1 H-NMR, 101 MHz for 13 C-NMR and 376 MHz for 19 F-NMR with TMS as the internal standard. Chemical shifts were expressed in δ (ppm) and coupling constants (J) in Hz. Multiplicity was indicated as follows: s (singlet), d (doublet), t (triplet), p (quintet), dd (doublet of doublets), brs (broad singlet), etc. Mass spectra were obtained on an Agilent 1100 Series LC/MSD Trap mass spectrometer (ESI-MS, Agilent, Santa Clara, CA, USA). 1,8-Dihydroxy-3-(2-hydroxyethoxy)-6-methylanthracene-9,10-dione (2) To a mixture of emodin (10.0 g, 37.0 mmol) in dry DMF (150 mL) were added Cs 2 CO 3 (13.2 g, 40.5 mmol) and 2-iodoethanol (19.1 g, 111 mmol) at room temperature. After stirring for 36 h at 60 • C, the resulting mixture was evaporated under reduced pressure and then mixed with water (500 mL). The pH value of aqueous phase was adjusted to around 5 with 10% hydrochloric acid solution. The yellow precipitate was collected and washed with water to give the crude product, which was in further purification by triturating twice with ethyl acetate (100 mL) and filtered to afford compound 2 (7.6 g, 65%) as a brown solid; 1 General procedure A for preparation of compounds 3a-3z, 4a-4b, and 7a-7l To a mixture of compound 2 (1.00 mmol) in dry dichloromethane (20 mL) were added various N-Boc amino acids (1.20 mmol), dicyclohexyl carbodiimide (DCC) (4.50 mmol) and 4-(N,N-dimethlyamino) pyridine (DMAP) (1.00 mmol) at 0 • C. After stirring about 30 min to 5 h at 0 • C, TLC analysis showed the complete consumption of compound 2, and then the resulting mixture was added dropwise TFA (5 mL) at the same temperature and kept stirring for anther about 2 h. The insoluble side product was filtered out and the filtrate was evaporated to give the residue, which was purified by reverse phase flash chromatography with the following conditions: Column: Spherical C18, 20-40 µm, 330 g; Mobile Phase A: Water (plus 5 mM TFA); Mobile Phase B: ACN; Flow rate: 80 mL/min; Gradient: 5% B gradient in 10 min, 25% B-45% B gradient in 25 min; Detector: 254 nm. The fractions containing the desired product were collected at around 40% B and concentrated under reduced pressure to afford compounds 3a-3z in 20% to 50% yields. 2-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-2-oxoethanaminium 2,2,2-trifluoroacetate (3a). According to the general procedure A, compound 2 was treated with N-Boc-glycine and then purified by reverse phase flash chromatography to give compound 3a: Yellow solid; yield, 45%; 1 (S)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-1-oxopropan-2-aminium 2,2,2-trifluoroacetate (3c). According to the general procedure A, compound 2 was treated with N-Boc-L-alanine and then purified by reverse phase flash chromatography to give compound 3c: Yellow solid; yield, 48%; 1 H-NMR (S)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-3-hydroxy-1-oxopropan-2-aminium 2,2,2-trifluoroacetate (3d). According to the general procedure A, compound 2 was treated with N-Boc-O-tert-butyl-L-serine and then purified by reverse phase flash chromatography to give compound 3d: Yellow solid; yield, 20%; 1 (R)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-3-methyl-1-oxobutan-2-aminium 2,2,2-trifluoroacetate (3e). According to the general procedure A, compound 2 was treated with N-Boc-D-valine and then purified by reverse phase flash chromatography to give compound 3e: Yellow solid; yield, 46%; 1 189.89, 181.04, 168.97, 164.68, 164.26, 161.46, 148.57, 134.85, 132.73, 124.23, 120.57, 113.37, 110.09, 107.74, 107.05, 66.65, 63.45, 57. 161.48, 148.62, 134.89, 132.78, 124.26, 120.61, 113.42, 110.13, 107.79, 107.11, 66.63, 63.71 1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-2-methyl-1-oxopropan-2aminium 2,2,2-trifluoroacetate (3j). According to the general procedure A, compound 2 was treated with N-Boc-Aib-OH and then purified by reverse phase flash chromatography to give compound 3j: Yellow solid; yield, 38%; 1 H-NMR δ 11.96 (br, 2H), 8.51 (br, 3H), 7.54 (s, 1H), 7.34-7.14 (m, 2H), 6 1-((2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)carbonyl)cyclopropan-aminium 2,2,2-trifluoroacetate (3k). According to the general procedure A, compound 2 was treated with 1-(Boc-amino)cyclopropanecarboxylic acid and then purified by reverse phase flash chromatography to give compound 3k: Yellow solid; yield, 28%; 1 H-NMR δ 11.96 (brs, 1H), 8 189.28, 180.10, 169.83, 164.69, 164.22, 161.39, 148.39, 134.15, 132.07, 123.98, 120.37, 112.73, 109.54, 107.67, 106.74, 66.57, 63.78 1-((2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)carbonyl)cyclobutan-aminium 2,2,2-trifluoroacetate (3l). According to the general procedure A, compound 2 was treated with 1-(Boc-amino)cyclobutanecarboxylic acid and then purified by reverse phase flash chromatography to give compound 3l: Yellow solid; yield, 28%; 1 H-NMR δ 11.95 (brs, 1H), 8 3-((2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)carbonyl)oxetan-3-aminium 2,2,2-trifluoroacetate (3m). According to the general procedure A, compound 2 was treated with 3-Boc-amino-3-oxetanecarboxylic acid and then purified by reverse phase flash chromatography to give compound 3m: Yellow solid; yield, 37%; 1 H-NMR δ 11.97 (brs, 2H), 8 1-((2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)carbonyl)cyclopentan-aminium 2,2,2-trifluoroacetate (3n). According to the general procedure A, compound 2 was treated with N-Boc-aminocyclopentanecarboxylic acid and then purified by reverse phase flash chromatography to give compound 3n: Yellow solid; yield, 31%; 1 H-NMR δ 11.98 (brs, 2H), 8 (R)-2-((2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)carbonyl)-pyrrolidinium 2,2,2-trifluoroacetate (3o). According to the general procedure A, compound 2 was treated with N-Boc-D-proline and then purified by reverse phase flash chromatography to give compound 3o: Yellow solid; yield, 44%; 1 H-NMR δ 7.52 (d, J = 1. 7 (S)-2-((2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)carbonyl)-pyrrolidinium 2,2,2-trifluoroacetate (3p). According to the general procedure A, compound 2 was treated with N-Boc-L-proline and then purified by reverse phase flash chromatography to give compound 3p: (R)-2-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-2-oxo-1-phenylethanaminium 2,2,2-trifluoroacetate (3q). According to the general procedure A, compound 2 was treated with N-Boc-D-phenylglycine and then purified by reverse phase flash chromatography to give compound 3q: Yellow solid; yield, 39%; 1 H-NMR δ 11.97 (brs, 2H), 8 (R)-2-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-1-(4-hydroxyphenyl)-2-oxoethanaminium 2,2,2-trifluoroacetate (3s). According to the general procedure A, compound 2 was treated with N-Boc-D-4-hydroxyphenylglycine and then purified by reverse phase flash chromatography to give compound 3s: Yellow solid; yield, 42%; 1 H-NMR δ 11.99 (brs, 2H), 9. (S)-2-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-1-(4-hydroxyphenyl)-2-oxoethanaminium 2,2,2-trifluoroacetate (3t). According to the general procedure A, compound 2 was treated with N-Boc-L-4-hydroxyphenylglycine and then purified by reverse phase flash chromatography to give compound 3t: Yellow solid; yield, 40%; 1 H-NMR δ 11.96 (brs, 1H), 9. (R)-2-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-1-(4-fluorophenyl)-2-oxoethanaminium 2,2,2-trifluoroacetate (3u). According to the general procedure A, compound 2 was treated with (R)-N-Boc-4-fluorophenylglycine and then purified by reverse phase flash chromatography to give compound 3u: Yellow solid; yield, 45%; 1 H-NMR δ 11.97 (brs, 2H), 8 (S)-2-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-1-(4-fluorophenyl)-2oxoethanaminium 2,2,2-trifluoroacetate (3v). According to the general procedure A, compound 2 was treated with (S)-N-Boc-4-fluorophenylglycine and then purified by reverse phase flash chromatography to give compound 3v: Yellow solid; yield, 43%; 1 H-NMR (400 MHz, DMSO-d 6 ) δ 11.95 (brs, 2H), 8 (R)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-1-oxo-3-phenyl-propan-2-aminium 2,2,2-trifluoroacetate (3w). According to the general procedure A, compound 2 was treated with N-Boc-D-phenylalanine and then purified by reverse phase flash chromatography to give compound 3w: Yellow solid; yield, 48%; 1 H-NMR δ 11.90 brs, 1H), 8 (S)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-1-oxo-3-phenyl-propan-2-aminium 2,2,2-trifluoroacetate (3x). According to the general procedure A, compound 2 was treated with N-Boc-L-phenylalanine and then purified by reverse phase flash chromatography to give compound 3x: Yellow solid; yield, 45%; 1 H-NMR δ 11.90 (brs, 1H), 8 (R)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-3-(4-hydroxyphenyl)-1-oxopropan-2-aminium 2,2,2-trifluoroacetate (3y). According to the general procedure A, compound 2 was treated with N-Boc-O-tert-butyl-L-tyrosine and then purified by reverse phase flash chromatography to give compound 3y: Yellow solid; yield, 26%; 1 H-NMR δ 9.33 (s, 1H), 7.54 (d, J = 1. 6 (S)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-3-(4-hydroxyphenyl)-1-oxopropan-2-aminium 2,2,2-trifluoroacetate (3z). According to the general procedure A, compound 2 was treated with N-Boc-L-tyrosine and then purified by reverse phase flash chromatography to give compound 3z: Yellow solid; yield, 23%; 1 H-NMR δ 11.04 (s, 1H), 7.50 (d, J = 5. 7 (R)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-N-methyl-1-oxo-propan-2-aminium 2,2,2-trifluoroacetate (4a). According to the general procedure A, compound 2 was treated with Boc-N-methyl-D-alanine and then purified by reverse phase flash chromatography to give compound 4a: Yellow solid; yield, 48%; 1 H-NMR δ 11.93 (brs, 2H), 9.18 (S)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-N-methyl-1-oxo-propan-2-aminium 2,2,2-trifluoroacetate (4b). According to the general procedure A, compound 2 was treated with Boc-N-methyl-L-alanine and then purified by reverse phase flash chromatography to give compound 4b: Yellow solid; yield, 45%; 1 H-NMR δ 12.15 (brs, 1H), 11.96 (brs, 1H), 9.11 General procedure B for preparation of compounds 5a and 5b To a solution of compound 4a or 4b (100 mg, 0.25 mmol) and paraformaldehyde (37% aqueous solution, 62 mg, 0.76 mmol) in MeOH (15 mL) was added NaBH 3 CN (24 mg, 0.38 mmol) at 0 • C. After stirring 2 h at room temperature, the reaction was quenched by TFA (0.1 mL). The resulting reaction solution was used directly in purification by reverse phase flash chromatography with the following conditions: Column: Spherical C18, 20-40 µm, 330 g; Mobile Phase A: Water (plus 5 mM TFA); Mobile Phase B: ACN; Flow rate: 80 mL/min; Gradient: 5% B gradient in 10 min, 25% B-45% B gradient in 25 min; Detector: 254 nm. The fractions containing the desired product were collected at around 40% B and concentrated under reduced pressure to afford compounds 5a or 5b in 85% or 89% yield, respectively. (R)-1-(2-(4,5-Dihydroxy-7-methyl-9,10-dioxo-9,10-dihydroanthracen-2-yloxy)ethoxy)-N,N-dimethyl-1-oxopropan-2-aminium 2,2,2-trifluoroacetate (5a). Yellow solid; yield, 85%; 1 General procedure C for preparation of compounds 6a-6f To a mixture of emodin (10 mmol) in dry DMF (50 mL) were added Cs 2 CO 3 (12 mmol) and hydroxybromides or iodides (30 mmol) at room temperature. After stirring for 36 h at 60 • C, the resulting mixture was evaporated under reduced pressure and then mixed with water (100 mL). The pH value of aqueous phase was adjusted to around 5 with 10% hydrochloric acid solution, extracted with dichloromethane (2 × 100 mL). The combined organic layer was washed with brine (200 mL), dried over anhydrous sodium sulfate and evaporated to dryness. The crude product was purified by silica gel column chromatography with 1%-10% ethyl acetate in petroleum to afford compounds 6a-6f. Total cell lysates from cultured HepG2 cells treated with different concentrations of compound 7a for 48 h were obtained by lysing the cells in ice-cold RIPA buffer (1 PBS, 1% NP-40, 0.5% sodium deoxycholate and 0.1% SDS) containing 100 mg/mL PMSF, 5 mg/mL aprotinin, 5 mg/mL leupeptin, 5 mg/mL pepstatin and 100 mg/mL NaF. After centrifugation at 12,000 rpm for 10 min, the protein in the supernatant was quantified by the Bradford method (BIO-RAD, Hercules, CA, USA) using a Multimode Varioscan instrument (Thermo Fischer Scientific, Waltham, MA, USA). Twenty micrograms of protein per lane was applied in 12% SDS polyacrylamide gel. After electrophoresis, the protein was transferred to a polyvinylidine difluoride membrane (Amersham Biosciences, Marlborough, MA, USA). The membrane was blocked at room temperature for 2 h in TBST containing 5% blocking powder (Santa Cruz, Dallas, TX, USA). The membrane was washed with TBST for 5 min, and the primary antibody was added and incubated at 4 • C overnight (O/N). Bax (5023, CST, Danvers, MA, USA), Bcl-2 (15071, CST), cytochrome c (4280, CST), procaspase-9 (AB138412, Abcam, Cambridge, MA, USA), procaspase-3 (AB32150, Abcam), caspase-3 (66470-2-lg, Proteintech Group, Rosemont, IL, USA) and GAPDH (AB8245, Abcam) antibodies were employed. After three TBST washes, the membrane was incubated with the corresponding horseradish peroxidase-labeled secondary antibody (1:5000) (Santa Cruz) at room temperature for 2 h. Membranes were washed with TBST for 15 min five times and the protein blots were visualized with chemiluminescence reagent (Thermo Fischer Scientific Ltd.). The X-ray films were developed with a developer and fixed with fixer solution. The grey levels were analyzed using ImageQuant LAS 4000 system (GE, Marlborough, MA, USA). The HepG2 cells were treated with indicated concentrations of compound 7a. After incubation for 48 h, cells were washed twice with ice-cold PBS, fixed and permeabilized with ice cold 70% ethanol at −20 • C overnight. The cells were treated with 100 µg/mL RNase A at 37 • C for 30min, then washed with ice-cold PBS and finally stained with 1mg/mL PI in the dark at 4 • C for 30 min. The cellular DNA content for the cell cycle distribution analysis was performed with the system software (CellQuest; BD Biosciences), plotting at least 30,000 events per sample. The percentage of cells in the G1, S and G2 phases of the cell cycle were determined using the ModFit LT version 5.0 software package (Verity Software, Topsham, ME, USA). HepG2 cells were grown in DMEM medium containing growth factors at a cell density of 1 × 10 5 cells/mL for 24 h. A disposable 200 mL plastic pipette tip was used to scratch the monolayer of cells in a streaking motion. Compounds were added to the streaked cell culture at the indicated concentrations. The streaked cells were then cultured in serum-free medium for an additional 48 h and photographed. To quantify the experimental results, the % cell inhibitory rate was calculated by the equation: cell inhibitory rate (%) = (1 − D drug /D control ) × 100%, where D drug is the mean distance of cell migration in drug group and D control is the mean distance of cell migration in control group. Pictures of the initial wounded monolayers were compared with the corresponding pictures of cells at the end of the incubation, and data were presented as mean ± SD from three independent experiments. In conclusion, a novel series of emodin derivatives via the introduction of an amino acid were designed and synthesized. Their in vitro anti-proliferation tests revealed that these derivatives exhibited moderate to potent anti-proliferative activity against HepG2 cells and MCF-7 cells. Among these compounds, the most potent compound, 7a, exhibited better selective anti-proliferative activity and specificity than emodin, displayed a significant effect in inducing cell cycle arrest at G0/G1 phase and inducing cell apoptosis in HepG2 cells via release of cytochrome c and subsequent activation of caspase-9 and caspase-3, which revealed that the possible molecular mechanism of 7a apoptosis induction may mainly through the mitochondrial death pathway. Moreover, Compound 7a also resulted in inhibition of HepG2 cells migration in the wound healing assay. These preliminary molecular mechanism results suggest that compound 7a could be a promising lead compound for the development of novel antitumor drugs and has the potential for further investigations as an anti-cancer drug. Supplementary Materials: The following are available online at http://www.mdpi.com/1420-3049/24/5/884/s1, Tables S1-S4: Biochemistry analytical data from three independent experiments; Figures S1-S140: 1 H-NMR, 13 C-NMR and 19 F-NMR spectra of these compounds. According to the general procedure C, emodin was treated with 3-iodopropan-1-ol and then purified by silica gel column chromatography to give compound 6a: Brown solid; yield, 55%; 1 H-NMR δ 12.14 (s, 1H), 11.96 (s, 1H), 7.51 (s, 1H According to the general procedure C, emodin was treated with 4-bromobutan-1-ol and then purified by silica gel column chromatography to give compound 6b: Brown solid; yield, 47%; 1 H-NMR δ 12.16 (s, 1H), 11.98 (s, 1H), 7.53 (s, 1H) According to the general procedure C, emodin was treated with 5-bromopentan-1-ol and then purified by silica gel column chromatography to give compound 6c: Brown solid; yield, 52%; 1 H-NMR δ 12.15 (s, 1H), 11.96 (s, 1H) According to the general procedure C, emodin was treated with 2-(2-bromoethoxy)ethanol and then purified by silica gel column chromatography to give compound 6d: Brown solid; yield, 38%; 1 H-NMR δ 12.14 (s, 1H), 11.95 (s, 1H), 7.50 (d, J = 1 -bromoethoxy)ethoxy)ethanol and then purified by silica gel column chromatography to give compound 6e: Brown solid; yield, 42%; 1 H-NMR δ 12.01 (s, 1H) -bromoethoxy)ethoxy)ethoxy) ethanol and then purified by silica gel column chromatography to give compound 6f: Brown solid H-NMR δ 12.13 (s, 1H), 11.94 (s, 1H), 7.49 (d, J = 1.6 Hz, 1H), 7.18 (t, J = 1.3 Hz, 1H), 7.15 (d, J = 2.5 Hz, 1H) According to the general procedure A, compound 6a was treated with N-Boc-D-alanine and then purified by reverse phase flash chromatography to give compound 7a: Yellow solid; yield, 70%; 1 H-NMR δ 7.52 (d, J = 1 According to the general procedure A, compound 6a was treated with N-Boc-L-alanine and then purified by reverse phase flash chromatography to give compound 7b: Yellow solid According to the general procedure A, compound 6b was treated with N-Boc-D-alanine and then purified by reverse phase flash chromatography to give compound 7c: Yellow solid According to the general procedure A, compound 6b was treated with N-Boc-L-alanine and then purified by reverse phase flash chromatography to give compound 7d: Yellow solid Hz, 1H), 2.44 (s, 3H), 1.92-1.76 (m, 4H), 1.41 (d, J = 7.2 Hz, 3H) According to the general procedure A, compound 6c was treated with N-Boc-D-alanine and then purified by reverse phase flash chromatography to give compound 7e: Yellow solid; yield, 70%; 1 H-NMR δ 12.00 (brs, 2H), 8.32 (brs, 3H), 7.54 (d, J = 1 According to the general procedure A, compound 6c was treated with N-Boc-L-alanine and then purified by reverse phase flash chromatography to give compound 7f: Yellow solid 43 (s, 3H), 1.83-1.76 (m, 2H), 1.75-1.68 (m, 2H), 1.56-1.46 (m, 2H), 1.40 (d, J = 7.1 Hz, 3H) According to the general procedure A, compound 6d was treated with N-Boc-D-alanine and then purified by reverse phase flash chromatography to give compound 7g: Yellow solid; yield, 65%; 1 H-NMR δ 9 According to the general procedure A, compound 6d was treated with N-Boc-L-alanine and then purified by reverse phase flash chromatography to give compound 7h: Yellow solid According to the general procedure A, compound 6e was treated with N-Boc-D-alanine and then purified by reverse phase flash chromatography to give compound 7i: Yellow solid; yield, 66%; 1 H-NMR δ 9 According to the general procedure A, compound 6e was treated with N-Boc-L-alanine and then purified by reverse phase flash chromatography to give compound 7j: Yellow solid; yield, 65%; 1 H-NMR δ 9 According to the general procedure A, compound 6f was treated with N-Boc-D-alanine and then purified by reverse phase flash chromatography to give compound 7k: Yellow solid; yield, 66%; 1 H-NMR δ 9.76 (brs, 3H), 7.46 (d, J = 1.6 Hz, 1H) According to the general procedure A, compound 6f was treated with N-Boc-L-alanine and then purified by reverse phase flash chromatography to give compound 7l: Yellow solid; yield, 65%; 1 H-NMR δ 9 ) containing 10% Fetal bovine serum (FBS) and 1% penicillin (100 units/mL)-streptomycin (100 µg/mL) in a humidified incubator at 37 • C under 5% CO 2 and 95% relative humidity (RH) atmosphere. The cells were harvested for metabolomics experiments using cell scrapers Cell Viability Assay HepG2 cells and MCF-7 cells in 40 µL DMEM supplemented with 10% FBS were seeded at a density of 500 cells/well in 384-well cell culture plate (3570, Corning), respectively. After 24 h, the medium was substituted with fresh medium containing various sample solutions in DMSO USA) was added to each well containing cells and the contents were mixed for 2 min on an orbital shaker to induce cell lysis. The plates were then incubated at room temperature for 20 min to stabilize the luminescent signal. The luminescence was recorded on an Envision Cells were seeded at 2 × 10 5 /well in 10% FBS-DMEM into 6-well plates and treated with compounds 7a for 48 h. The cells were washed twice with cold phosphate buffered saline (PBS) and then Analysis is as follows: lower left quadrant, viable cells (Annexin V−/PI−); lower right quadrant, early apoptotic cells (Annexin V+/PI−); upper right quadrant, late apoptotic cells (annexin V+/PI+); upper left quadrant, necrotic cells A new diphenyl ether from marine-derived fungus Aspergillus sp B-F-2 A: A new antimicrobial anthraquinone from a sea urchin-derived fungus Monodictys sp A New Ergosterol Analog, a New Bis-Anthraquinone and Anti-Obesity Activity of Anthraquinones from the Marine Sponge-Associated Fungus Talaromyces stipitatus KUFA 0207. Mar. Drugs Bis-Indolyl Benzenoids, Hydroxypyrrolidine Derivatives and Other Constituents from Cultures of the Marine Sponge-Associated Fungus Aspergillus candidus KUFA0062 Emodin and Aloe-Emodin Suppress Breast Cancer Cell Proliferation through ER alpha Inhibition. Evid.-Based Complement Emodin suppresses wnt signaling in human colorectal cancer cells SW480 and SW620 Emodin Suppresses Hyperglycemia-Induced Proliferation and Fibronectin Expression in Mesangial Cells via Inhibiting cFLIP Emodin suppresses inflammatory responses and joint destruction in collagen-induced arthritic mice Emodin blocks the SARS coronavirus spike protein and angiotensin-converting enzyme 2 interaction Effect of emodin on the cariogenic properties of Streptococcus mutans and the development of caries in rats Emodin, a naturally occurring anthraquinone derivative, suppresses IgE-mediated anaphylactic reaction and mast cell activation CT Imaging Biomarker for Evaluation of Emodin as a Potential Drug on LPS-mediated Osteoporosis Mice Formula Optimization of the Jiashitang Scar Removal Ointment and Antiinflammatory Compounds Screening by NF-kappa B Bioactivity-guided Dual-luciferase Reporter Assay System Emodin Protects against Diabetic Cardiomyopathy by Regulating the AKT/GSK-3 beta Signaling Pathway in the Rat Model Emodin Prolongs Recipient Survival Time After Orthotopic Liver Transplantation in Rats by Polarizing the Th1/Th2 Paradigm to Th2 Emodin induces neurite outgrowth through PI3K/Akt/GSK-3 beta-mediated signaling pathways in Neuro2a cells Emodin Ameliorates Ethanol-Induced Fatty Liver Injury in Mice Emodin Protects Against Concanavalin A-Induced Hepatitis in Mice through Inhibiting Activation of the p38 MAPK-NF-kappa B Signaling Pathway Emodin suppresses cell proliferation and fibronectin expression via p38MAPK pathway in rat mesangial cells cultured under high glucose Emodin accelerates osteoblast differentiation through phosphatidylinositol 3-kinase activation and bone morphogenetic protein-2 gene expression Emodin-induced apoptosis through p53-dependent pathway in human hepatoma cells Emodin azide methyl anthraquinone derivative triggers mitochondrial-dependent cell apoptosis involving in caspase-8-mediated Bid cleavage Aloe-emodin induces apoptosis of human nasopharyngeal carcinoma cells via caspase-8-mediated activation of the mitochondrial death pathway Emodin inhibits TNF-alpha-induced human aortic smooth-muscle cell proliferation via caspase-and mitochondrial-dependent apoptosis Emodin induces apoptosis in human lung adenocarcinoma cells through a reactive oxygen species-dependent mitochondrial signaling pathway Synthesis and antitumor activity of emodin quaternary ammonium salt derivatives Synthesis and biological activity, evaluation of emodin quaternary ammonium salt derivatives as potential anticancer agents Synthesis, SAR and pharmacological characterization of novel anthraquinone cation compounds as potential anticancer agents Antitumor Effects and Mechanism of Novel Emodin Rhamnoside Derivatives against Human Cancer Cells In Vitro Synthesis, DNA binding and cytotoxicity of new pyrazole emodin derivatives Synthesis and cytotoxic activities of beta-carboline amino acid ester conjugates Preparation of amino acid conjugates of betulinic acid with activity against human melanoma This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license The authors declare no conflict of interest.