- Research article
- Open Access
Identification of the anti-breast cancer targets of triterpenoids in Liquidambaris Fructus and the hints for its traditional applications
BMC Complementary Medicine and Therapies volume 20, Article number: 369 (2020)
Liquidambaris Fructus is the infructescences of Liquidambar formosana Hance and it has been used to treat some breast disease in Traditional Chinese Medicine. In the previous study we found the anti-breast cancer effect of triterpenoid in Liquidambaris Fructus. This study is a further investigation of the triterpenoids in Liquidambaris Fructus and aims to identify their anti-breast cancer targets, meanwhile, to estimate the rationality of the traditional applications of Liquidambaris Fructus.
Triterpenoids in Liquidambaris Fructus were isolated and their structures were identified by NMR spectrums. Potential targets of these triterpenoids were predicted using a reverse pharmacophore mapping strategy. Associations between these targets and the therapeutic targets of breast cancer were analyzed by constructing protein-protein interaction network, and targets played important roles in the network were identified using Molecular Complex Detection method. Binding affinity between the targets and triterpenoids was studied using molecular docking method. Gene ontology enrichment analysis was conducted to reveal the biological process and signaling pathways that the identified targets were involved in.
Thirteen triterpenoids were identified and 6 of them were the first time isolated from Liquidambaris Fructus. Predicted ADME properties revealed a good druggability of these triterpenoids. We identified 18 protein targets which were closely related to breast cancer progression, especially triple-negative, basal-like or advanced stage breast cancers. The triterpenoids could bind with these targets as their inhibitors: hydrophobic skeleton is a favorable factor for them to stabilize at binding site and polar C17- or C3- substituent was necessary for binding. GO enrichment analysis indicated that inhibition of protein tyrosine kinases autophosphorylation might be the primary mechanism for the anti-breast cancer effect of the triterpenoids, and ErbB4 and EGFR were the most relevant targets.
The study revealed that triterpenoids from Liquidambaris Fructus might exert anti-breast cancer effect by directly inhibit multiple protein targets and signaling pathways, especially ErbB4 and EGFR and related pathways. This study also brings up another hint that the traditional applications of Liquidambaris Fructus on hypogalactia should be reassessed systematically because it might suppress rather than promote lactation by inhibiting the activity of ErbB4.
Liquidambaris Fructus (LF), also called Lu Lu Tong, is the dry infructescences of Liquidambar formosana Hance. LF is believed to unclog the meridians (Jingluo) in Traditional Chinese Medicine, and because it mainly affect the liver meridian which is considered closely related to mammary function, therefore LF has been used to treat some breast disease such as hypogalactia , mastitis , mammary duct ectasia , mammary gland hyperplasia , etc. Breast cancer is the severest breast disease threatened women’s health as its high incidence rates and mortality rates . Many natural products, including Chinese herbal medicines, have been used in breast cancer treatment, and pharmacological evidence has been accumulated for their application. In the previous research we found that betulonic acid, a triterpenoid come from LF, could inhibit human breast cancer MCF-7 cells proliferation, and induce the cell cycle arrest and apoptosis . Hence it is necessary to further investigate the chemical ingredients and the mechanism of anti-breast cancer effect of LF.
Therefore, in the present study, 13 triterpenoids were identified and 6 of them were the first time isolated from LF. Considering that vast majority of the compounds reported in LF were triterpenoids [7,8,9,10,11,12], investigating macromolecular targets of LF triterpenoids could provide valuable information. Accordingly we adopted PharmMapper server as the computational tool which utilize a large-scale reverse pharmacophore mapping strategy . On the other hand, we also collected the known and explored therapeutic protein targets of breast cancer. Because proteins do not function in isolation, and proteins which interact with the disease proteins may also play an important role in disease progression , thus we constructed protein-protein interaction (PPI) network to display associations between the predict targets of LF triterpenoids and the therapeutic targets of breast cancer. However, this guilt-by-association method generated a large PPI data set that needs further analysis to provide more effectively information . Because cellular processes are mainly carried out by protein complexes , we identified the protein complexes from the PPI network using Molecular Complex Detection (MCODE) method, a graph theoretic clustering algorithm that detects densely connected regions that may represent molecular complexes [15, 17]. Molecular docking method could be used to analyses the conformations and orientations of the ligand binding to protein receptor, and predict molecular interactions and calculate binding energies , thus it was adopted to identify the targets that had high binding affinity to LF triterpenoids. Gene ontology (GO) enrichment analysis was finally conducted to reveal the biological process and signaling pathways that the identified targets were involved in. By means of the described methodology, this research could help us to acquire a better understanding about the anti-breast cancer potential and the traditional applications of LF.
Chemical ingredients investigation
Open column chromatography (CC) separation was carried out using silica gel (200–300 mesh) and Sephadex LH-20 as stationary phase. NMR spectra were run on a Bruker AVANCE II 600 spectrometer (600 MHz for 1H NMR and 150 MHz for 13C NMR) with trimethylsilane as an internal standard. LF was collected in Jiangsu Province of China in February, 2014, and was identified by one of the authors, Prof. Dong-Fang Zhang. Voucher specimens have been deposited in the Department of pharmacognosy, School of Pharmacy, China Medical University (Liaoning, China).
LF (9 kg) was grinded and refluxed with 70% ethanol, and then concentrated under reduced pressure to afford a residue 430 g. The residue was dispersed in water and successively extracted with cyclohexane and EtOAc, and 13 compounds were isolated from the EtOAc extraction (separation process was shown in Additional file 1). These compounds were determined using 1H-NMR and 13C-NMR spectrum analysis, and were identified by comparing the results with literature [7,8,9,10,11,12, 19,20,21,22,23,24, 75, 76].
Absorption and metabolism properties of LF triterpenoids were predicted using admetSAR (http://lmmd.ecust.edu.cn/admetsar1) . Potential targets of LF triterpenoids were predicted using PharmMapper server (http://www.lilab-ecust.cn/pharmmapper/) with the target set Human Protein Targets Only (v2010, 2241) . The known and explored therapeutic targets of breast cancer (ICD10: C50) were collected from the Therapeutic Target Database (http://bidd.nus.edu.sg/group/cjttd/) . Target gene names were normalized using the Uniprot database (https://www.uniprot.org/) . PPI network was constructed using the STRING database (version 10.5, https://string-db.org/) , and the network data was processed using Cytoscape software (version 3.6.1)  and the isolated nodes were removed. Three parameters, Degree, Closeness Centrality and Betweenness Centrality were used to assess the topological properties of the nodes in the PPI network, and higher values means more topological importance. The nodes whose values were less than the medians were defined as non-critical nodes and were removed from the PPI network. Then the network was analyzed using MCODE (version 1.5.1) , and the values of Degree cutoff, Node Score Cutoff, K-Core, Max Depth were 2, 0.2, 2, 100, respectively. Four clusters were obtained and Cluster 4 was abandoned because it had only 3 nodes and contained no therapeutic targets. Cluster 3 was merged into cluster 2. Then cluster 1 and 2 were used to create the compounds−predicted targets−therapeutic targets interaction networks, respectively.
Crystal structures of protein targets were collected from the Protein Data Bank (PDB, http://www.rcsb.org/)  and inhibitors or agonists contained in the crystal structures were set as reference ligands. Besides, approved drugs or clinical trial drugs of the targets were collected from the Therapeutic Target Database and were also set as reference ligands (Additional file 2). Molecular docking experiments were performed at Molecular Operating Environment (MOE, 2018.01, Chemical Computing Group ULC). Pocket atoms were set as docking sites. A triangle matcher method was used as placement method. London dG scoring was used to select the top 30 poses (conformations and orientations of the ligand). Selected poses were refined in induced fit mode with the sidechains set as ‘free to move’ to achieve flexible receptor docking. GBVI/WSA dG scoring was finally used to select 5 poses with the lowest binding free energy, so a lower score means a more stable conformation.
GO enrichment analysis
GO enrichment analysis was performed using ClueGO (version 2.5.2)  and CluePedia (version 1.5.2) . GO biological process and GO molecular function were employed as the annotation sources and result show only enriched terms with p value < 0.01. Network was created with kappa statistics and reflected the similarity of the terms’ gene memberships. 0.4 kappa score was applied as threshold to create functional groups which were represented by their most significant term (leading term).
Isolation and identification of LF triterpenoids
Thirteen triterpenoids were isolated from LF and were identified as betulonic acid (LF01), 28-nor-β-amyrenone (LF02), ambronal (LF03), betulin (LF04), oleanonic acid (LF05), ursonic acid (LF06), liquidambaric lactone (LF07), betulinic acid (LF08), epibetulinic acid (LF09), 3-oxo-12β-hydroxy-oleanan-28,13β-olide (LF10), hydroxyoleanonic lactone (LF11), 3α-acetoxy-25-hydroxyolean-12-en-28-oic acid (LF12), Arjunolic acid (LF13), respectively. Six of them, LF03, LF04, LF05, LF08, LF09 and LF10 (NMR spectrums were shown in Additional file 3), were the first time isolated from LF. Other known triterpenoids contained in LF, including 6β-Hydroxy-3-oxo-lup-20(29)-en-28-oic acid (LFa), Oleanolic acid (LFb), Erythrodiol (LFc), 11α-Methoxy-28-nor-β-amyrenone (LFd), 2α,3β-dihydroxy-23-norolean-4(24), 12(13)-dien-28-oic acid (LFe), lantanolic acid (LFf) and ursolic acid (LFg), were collected from literature [7,8,9,10,11,12]. Among the 20 triterpenoids, LF01, LF04, LF08, LF09 and LFa belong to lupane-type triterpenoids, and LF06, LFg belong to ursane-type triterpenoids, and the others are oleanane-type triterpenoids (Fig. 1).
Absorption and metabolism properties of LF triterpenoids
Predicted ADME properties indicated that LF triterpenoids were appropriate for oral administration as they were classified as HIA+ (human intestinal absorption was no less than 30%) and Caco2+ (Caco-2 permeability value (Papp) > 8×10− 6 cm/s) except LF12 and LF13 (HIA+, Caco2-). Although LF triterpenoids were labeled as P-glycoprotein substrates, LF02, LF03, LF04, LF07, LF10, LF11 and LFd were also potential P-glycoprotein inhibitors, which might facilitate the transport of other triterpenoids. With respect to the predicted metabolism properties, all the 20 triterpenoids were the substrates of CYP450 3A4 and they did not inhibit the CYP450 metabolic activities.
Exploration of anti-breast cancer targets of LF triterpenoids
Potential targets of LF triterpenoids were ranked according to the normalized fit score, and the top 100 targets were selected for each compound, then these targets were merged and 174 targets were finally obtained. On the other hand, 111 therapeutic targets (including 27 subtypes) of breast cancer were collected (Additional file 4). Based on these targets, 2 optimized compounds−predicted targets−therapeutic targets interaction networks were constructed (Fig. 2). The network contained 36 predicted targets and 35 therapeutic targets, which had 11 targets in common. The result indicated that the 36 predicted targets remained in the network were more closely related to breast cancer progression than the others.
Molecular docking experiment was carried out on LF triterpenoids and the selected 36 targets except for calmodulin-1 (CALM1), eNOS (NOS3), serum albumin (ALB), Mn-SOD (SOD2), BMP2 and placenta growth factor (PGF) because no binding site for reference. Molecular docking results were summarized in Additional file 5 and Fig. 3. Eighteen of the 30 targets, including ErbB4 (ERBB4), EGFR, CDK2, CDK6, p38 MAPK (MAPK14), ERK2 (MAPK1), Lck (LCK), Src (SRC), IGF1R, PI3Kγ (PIK3CG), VEGFR2 (KDR), aromatase (CYP19A1), PTP1B (PTPN1), SHP2 (PTPN11), caspase-3 (CASP3), PARP-1, HSP90α (HSP90AA1) and HSC70 (HSPA8), were identified as the anti-breast cancer targets of LF triterpenoids. To acquire a better understanding about the results, binding characteristics and functions of these targets were analyzed as follows.
Targets that belong to protein kinases family
Half of the 30 targets belong to protein kinase superfamily, among which 14 targets were protein tyrosine kinases or CMGC serine/threonine protein kinases (Additional file 6). These two types of protein kinases play important roles in cancer progression and they are promising targets in breast cancer therapy [31, 32]. The catalytic domain of protein tyrosine kinases have high homology with that of serine/threonine protein kinase , which may explain why LF triterpenoids showed the features of ATP-competitive inhibitors targeting both the two types of protein kinases.
EGFR subfamily members (belonging to receptor tyrosine kinases) play critical role in the pathogenesis of breast cancer. Besides the well-known EGFR (ErbB1/HER1) and ErbB2 (HER2, ERBB2), there is growing concern about another family member ErbB4 (HER4, ERBB4). Clinical studies suggested that to the patients with advanced stage or triple-negative breast cancer, increased expression of ErBb4 is associated with poor prognosis [34, 35]. Molecular docking result showed that hydrophobic skeleton is a favorable factor for LF triterpenoids to stabilize at the ATP binding site of ErbB4 (Fig. 4a). For active protein kinases, K/E/D/D signature motif plays important structural and catalytic roles . It can be observed that Lys751 (the K of K/E/D/D in ErbB4)  formed a salt bridge with Asp861, therefore the -NH3+ group of Lys751 could be a stable hydrogen bond donor to LF triterpenoids (Fig. 4b). For most of the LF triterpenoids, hydrogen bond donor electrostatically preferred locations were at the opposite ends of skeleton: substituent located at C-17 and C-3 (Fig. 4c, Table 1). In the LF01-ErbB4 complex, it was the 17-carboxyl group bound to Lys751 (Fig. 4d); in the other LF triterpenoids-ErbB4 complexes, however, it was the C-3 substituent (Table 1). Among the three types of LF triterpenoids, lupane-type triterpenoids could form more stable hydrogen bonds with Lys751 because of the low binding energies, and the most stable hydrogen bond came from LF01. However, LF04 could not form strong interaction with Lys751 might because the 17-carboxyl group was displaced by hydroxymethyl group. LF triterpenoids bound to EGFR in a similar fashion, but Lys745 (the K of K/E/D/D in EGFR)  bound to the 17-carboxyl group more often than the C-3 substituent, and the most stable hydrogen bonds came from LF05 and LFg (Table 1). However, lupane-type triterpenoids did not show strong interactions with EGFR, which meant that lupane-type triterpenoids from LF might have higher affinity and inhibition to ErbB4 than EGFR.
CDKs (belonging to CMGC serine/threonine protein kinase family) are also important anti-breast cancer targets. Although selective CDK4/6 inhibition results in less toxicity, a huge challenge is that breast cancer (including basal-like breast cancer) adapt quickly to CDK4/6 inhibition, while inhibiting CDK2 simultaneously is an effective solution [38, 39]. Molecular docking result showed that most LF triterpenoids could bind with both CDK6 and CDK2. In LF triterpenoids-CDK6 complexes, binding sites were catalytic loop (especially was Lys147) and activation segment (especially was Asp163, the D of DFG and the second D of K/E/D/D). In LF triterpenoids-CDK2 complexes, however, binding sites mainly were hinge residues (Asp86, Gln85, His84 and Leu83) and the K of K/E/D/D (Lys33) (Table 2). Because hinge region played an important role in inhibitor binding , the result indicated that LF triterpenoids probably had a higher inhibitory activity to CKD2 than CDK6.
It was reported that phosphorylated MAPKs including ERK2 (MAPK1), JNK1 (MAPK8) and p38 MAPK (MAPK14) over-expressed in breast cancer patients and is associated with poor prognosis [41,42,43]. Molecular docking result showed that LF triterpenoids could bind to the hinge region (Thr106, His107, Met109 and Gly110) and the DFG motif (Asp168 and Phe169) of p38 MAPK, but in LF triterpenoids-ERK2 complexes the interactions were weaker, and LF triterpenoids could not form hydrogen bonds with the hinge region of JNK1 (Table 3), indicating LF triterpenoids had higher inhibitions to p38 MAPK than ERK2 and JNK1. Lck (LCK) plays a pathological role in breast cancer progression probably by activating MAPK signaling pathway [44,45,46]. Molecular docking result showed that LF triterpenoids could bind with key residues of Lck including the conserved ATP binding lysine (Lys273), the linker sequence between the SH2 domain and the catalytic domain (Ser323), the catalytic domain (Ala368 and Asn369) and the activation loop (Asp382, Glu390, Tyr394, Thr395, and Arg397), which meant that LF triterpenoids could compete with ATP for binding Lck and inhibit its activity.
In addition, Src (SRC), VEGFR2 (KDR), IGF1R and PI3Kγ (PIK3CG) have already been successful targets in breast cancer therapy. Molecular docking result showed that LF triterpenoids could bind with these protein kinases and act as ATP competitive inhibitors.
Targets that belong to nuclear hormone receptor family and cytochrome P450 family
For nuclear hormone receptors ERα (ESR1), AR and PPARγ (PPARG), both agonist and antagonist bind at the same site of LBD, therefore it’s hard to figure out whether the ligand were agonist or antagonist through molecular docking result. Besides, AR and PPARγ are still complicated targets for breast cancer treatment because studies gave contradictory conclusions [47, 48]. Hence effects of LF triterpenoids caused by binding with these 3 targets need further studies. Human cytochrome P450 aromatase (AROM, CYP19A1) use androgens as substrates with high specificity and catalyze them to estrogens , and its inhibitors, such as exemestane, were used to treat breast cancer for many years. Molecular docking result showed that LF triterpenoids could compete with testosterone for binding aromatase, forming hydrogen bonds with Met374 and Arg115 in the catalytic cleft, and it is remarkable that LF triterpenoids had more stable conformations (with scores from − 8.1861 to − 9.3923) than testosterone (with scores from − 6.5824 to − 7.8243) and exemestane (scores − 7.7675) in the binding pocket, indicating that LF triterpenoids might be promising inhibitors of aromatase.
Although PTP1B (PTPN1) was previously known could dephosphorylate some receptor protein tyrosine kinases, studies showed that inhibition of PTP1B could be a new therapeutic strategy for breast cancer, including triple-negative breast cancer [50, 51]. The catalytic domain of PTP1B can be divided into four sites and A site is catalytic pocket. However because of the high homology of the catalytic domain between PTP1B and TCPTP (a phosphatase related to regulating T-cell activation), inhibitors only binding to A site have low selectivity [52, 53]. Molecular docking result showed that LF triterpenoids could bind on A site (Ala217, Gly220, Arg221, Gln266), C site (Tyr46, Arg47 and Asp48) and D site (Glu115, Lys116, Asp181). Lys116 on D site was a preferred residue to half of LF triterpenoids, such as LF05, LF06 and LF12 (the lowest binding energy was-9.1, − 8.3 and − 9.0 kcal/mol, respectively), but could not form bond with LF triterpenoids without 17-carboxyl group. Nearly half of LF triterpenoids could bind with more than one site, such as LF12 (occupied AC sites) and LFa (occupied ACD sites) (Fig. 5a, b). These characters of LF triterpenoids probably could increase their selectivity and inhibitory activity to PTP1B.
SHP2 (PTPN11) is required for the full activation of MAPK signaling pathway and could promotes breast cancer (including triple-negative breast cancer) progression [54, 55]. Important binding sites for SHP2 inhibitors include P-loop (the PTP signature motif for recognition of phosphoryl moiety in substrate), WPD loop (lid of the active site pocket) and pTyr recognition loop . Molecular docking result of SHP2 showed that LF triterpenoids could bind with Cys459 and Arg465 in P-loop (such as LF03), and Gly427 in the flexible WPD loop, and form C-H···π interaction with Tyr279 in pTyr recognition loop (such as LF04) (Fig. 5c, d). Besides, about half of LF triterpenoids could form hydrogen bond with the side chain of Lys366, such as LF05, LF12 and LFg (the lowest binding energy was − 8.5, − 8.2 and − 9.9 kcal/mol, respectively).
Targets related to apoptosis
Although the proteolytic enzyme caspase-3 (CASP3) could initiates apoptosis by cleaving and inactivating PARP1, it also could stimulate survived breast cancer cell growth by promoting PGE2 release . Molecular docking result showed that LF triterpenoids could bind with caspase-3 like its inhibitors. However, LF triterpenoids could promote apoptosis by directly inhibiting PARP1. PARP1 is a base excision repair protein that could promote cancer cell survival, and is an effective target in triple-negative breast cancer [58, 59]. Molecular docking result showed that C-17 substituent of LF triterpenoids could form hydrogen bound with the backbone of Arg878 (or the sidechain of Asp770) on adenine ribose binding site (AD site) of PARP1, such as LF08 (Fig. 6a). And binding the AD site is considered to gain selectivity for the inhibitors . Besides, some LF triterpenoids could form C-H···π interaction with Tyr896 in nicotinamide binding pocket (NI site), or form hydrogen bound with Asp766 and Glu763 on phosphate binding site (PH site), such as LFc (Fig. 6b). The result indicated that LF triterpenoids could compete with NAD+ for binding PARP1 then inhibit its activity.
Molecular chaperones HSP90α (HSP90AA1) and HSC70 (HSPA8) could assist the stability and function of many oncogenic proteins then keep cancer cells surviving . Molecular docking result showed that LF triterpenoids could bind with HSP90α and HSC70 at ATP binding site, indicating they may act as HSPs’ inhibitors then promote cell death.
The p53 protein can initiate apoptosis while its negative regulator MDM2 could bind and inactive p53 thus acts as a tumor promoter . Molecular docking result showed that LF triterpenoids could occupy the hydrophobic pockets of MDM2 that bind with p53, and formed hydrogen bonds with residues Met62, Tyr67, Tyr100, or formed C-H···π/O-H···π interactions with His96, but these interactions were not strong enough (binding energy − 0.6 ~ − 2.9 kcal/mol) to compete with the macromolecule p53, indicating LF triterpenoids might be not high affinity MDM2 antagonists.
HPGDS catalyzes the conversion of PGH2 to PGD2, and closely related to the progression of inflammation. Although its role in breast cancer is not entirely clear, it is reported that HPGDS suppress gastrointestinal cancer progression [63, 64]. However, molecular docking result showed that LF triterpenoids (except for LF01 and LF04) could form hydrogen bonds and C-H···π/O-H···π interactions with Arg14 and Trp104 of HPGDS, which are important residues for binding affinity and binding pocket stability, indicating LF triterpenoids could compete with PGH2 for binding HPGDS and acting as an inhibitor.
For other targets including annexin A5 (ANXA5), prothrombin (F2) and renin (REN), there was no firm evidence to link their functions directly to breast cancer progression, and they were less significant in the PPI network, which meant that they might be not important targets for LF triterpenoids.
GO enrichment analysis
GO enrichment analysis was conducted on the 18 identified targets according to the molecular docking result. As a result, totaling 29 terms were acquired, and they could be clustered into 7 functional groups (Additional file 7). The biggest group that covered 44.19% of enriched terms was peptidyl-tyrosine autophosphorylation group, indicating that inhibition of protein tyrosine kinases autophosphorylation might be the primary mechanism for the anti-breast cancer effect of LF triterpenoids. Figure 7 revealed that the most relevant targets for peptidyl-tyrosine autophosphorylation functional group were ErbB4 and EGFR, because they had more connections (each had 9 edges) with terms belonging to this group than the other targets. Besides, targets of LF triterpenoids also involved in response to epidermal growth factor, regulation of cyclin-dependent protein serine/threonine kinase activity, insulin receptor binding, regulation of mitochondrial depolarization, regulation of telomerase activity, vascular endothelial growth factor receptor signaling pathway, which were closed related to breast cancer progression.
According to the above result, we then analyzed signaling pathways concerned with the identified targets based on related studies [45, 46, 50, 51, 55, 57] and KEGG pathway database. By binding with the identified targets, LF triterpenoids could inhibit ErbB signaling pathway, MAPK signaling pathway, PI3K/Akt signaling pathway and VEGF signaling pathway, then inhibit cancer cells proliferation, and induce the cell cycle arrest and apoptosis (Fig. 8), which could explain the result of our previous research . There have been studies showing that betulonic acid and other LF triterpenoids could induce cancer cell apoptosis and cell cycle arrest via inhibiting PI3K/Akt pathway [65,66,67], which could also support the predicted mechanism.
According to the concept of polypharmacology, for complex diseases like cancers, inhibition of single target often induces chemo-resistance and lead to therapy failure [68,69,70]. In previous study we found the anti-breast cancer effect of betulonic acid , and this study revealed that LF triterpenoids might exert anti-breast cancer effect by directly inhibit multiple protein targets and signaling pathways, especially ErbB4 and EGFR and related pathways. It is worth noting that many of these targets were investigated or potential therapeutic targets of triple-negative, basal-like or advanced stage breast cancers, and patients with these types of breast cancer have rapid progression, high mortality rates, and high risk of recurrence and metastasis [59, 71]. Therefore the multi-targeted activity of LF triterpenoids revealed in this study could bring some inspirations in breast cancer treatment (Additional file 8).
In addition, this study left an extra hint about the traditional applications of LF: as a high binding affinity target of LF triterpenoids, ErbB4 activity is essential for mammary gland development and the initiation of lactation , which meant that LF might aggravate rather than alleviate hypogalactia, contradicted with the description in Chinese Pharmacopoeia. Actually in early literature that involved LF, such as Bencao Gangmu Shiyi (1765 AD), no record showed that LF could promote lactation. Therefore this usage probably was an extension of its original efficacy because hypogalactia is considered a symptom caused by blocking of meridians. However, this extended usage need more evidence and this study gave a negative conclusion (Additional file 8).
Compared to synthetic drugs, Traditional Chinese Medicine is generally considered to have better safety. Although LF has been used for more than 250 years and no obvious adverse effects were reported, based on the result of the present study, there are still some precautions should be noted. Because ErbB4 play essential roles in the development of cardiac muscle and central nervous system during pregnancy , LF triterpenoids might have adverse effect on fetal development by inhibiting the activity of ErbB4. Besides, because Lck activation could initiate T cell antigen receptor signaling, and is a critical step in immune response , therefore LF triterpenoids probably could cause immune suppression by inhibiting Lck, and should be used very carefully on patients with immunodeficiency.
To further investigate the chemical ingredients of LF, 13 triterpenoids were identified and 6 of them were the first time isolated from LF. Using network pharmacology-based analysis in combination with molecular docking method, 18 protein targets, ErbB4, EGFR, CDK2, CDK6, p38 MAPK, ERK2, Lck, Src, IGF1R, PI3Kγ, VEGFR2, aromatase, PTP1B, SHP2, caspase-3, PARP-1, HSP90α and HSC70, were identified and by binding with these targets, LF triterpenoids could inhibit ErbB signaling pathway, MAPK signaling pathway, PI3K/Akt signaling pathway and VEGF signaling pathway, then inhibit breast cancer cells proliferation, and induce the cell cycle arrest and apoptosis. GO enrichment analysis indicated that ErbB4 and EGFR and related pathways might play more important roles during the process.
On the other hand, this study also revealed that LF triterpenoids could inhibit the activity of ErbB4 and accordingly suppressed rather than promoted lactation. Therefore the traditional applications of LF on hypogalactia should be reassessed systematically. Besides, this study also indicated that LF must be used very cautiously to treat women during pregnancy or patients with immunodeficiency, because LF triterpenoids might cause embryo toxicity and immune suppression by inhibiting the activity of ErbB4 and Lck.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
The combination of CDKs
GSKs and CLKs
Epidermal growth factor receptor
Estrogen receptor α
Extracellular signal-regulated kinase 2
Hematopoietic prostaglandin D synthase
Lymphocyte cell-specific protein-tyrosine kinase
Mitogen-activated protein kinase
Protein tyrosine phosphatases
Src-homology 2 domain-containing phosphatase 2
T-cell protein tyrosine phosphatase
Chinese Pharmacopoeia Commission. Pharmacopoeia of the People’s Republic of China, vol. I. Beijing: China Medical Science and Technology Press; 2015.
Li Y, Gong X, Ding X. Research on the law of Chinese medicine in the treatment of early lactation acute mastitis. Chin J Ethnomed Ethnopharm. 2019;28:8–10,16.
Xu J, Wang F, Li Z, Chen J. Randomized parallel control study of Chupixiaoyi granule combined drug lavage treatment on nipple discharge mammary duct ectasia. J Pract Tradit Chin Intern Med. 2017;31:70–3.
Du X. The flat by Sanjie powder in the treatment of breast cystic hyperplasia randomized parallel controlled study. J Pract Tradit Chin Intern Med. 2013;27:12–3.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Mu XT, Qian P, Jiang LL, Wang Y, Jin X, Zhang DF. Effect of betulonic acid on proliferation of human breast cancer MCF-7 cells and human cervical cancer C-33A cells. Pract Pharm Clin Remedies. 2017;20:254–7.
Lai Z, Dong Y. Studies on the chemical constituents of Lu Lu Tong (1) structure determination of a novel pentacyclic triterpene. Acta Sci Nat Univ Sunyatseni. 1996;35:64–9.
Lai Z, Dong Y. Studies on the chemical constituents of Lu Lu Tong. (2). Acta Sci Nat Univ Sunyatseni. 1996;35:40–4.
Sun Y, Sun Y. Structure elucidation of liquidambaric lactone. Acta Pharm Sin. 1996;31:437–40.
Li C, Sun Y, Sun Y. A new triterpene from the fruits of Liquidambar formosana Hance-Lulutong. J Chin Pharm Sci. 2002;11:1–4.
Dat NT, Lee IS, Cai XF, Shen G, Kim YH. Oleanane triterpenoids with inhibitory activity against NFAT transcription factor from Liquidambar formosana. Biol Pharm Bull. 2004;27:426–8.
Shang H, Wang W, Li D, Hua H, Li Z. A new triterpenoid isolated from fruits of Liquidambaris Fructus. Chin Tradit Herb Drugs. 2014;45:1207–10.
Wang X, Shen Y, Wang S, Li S, Zhang W, Liu X, et al. PharmMapper 2017 update: a web server for potential drug target identification with a comprehensive target pharmacophore database. Nucleic Acids Res. 2017;45:W356–60.
Ideker T, Sharan R. Protein networks in disease. Genome Res. 2008;18:644–52.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.
Gavin AC, Bösche M, Krause R, Grandi P, Marzioch M, Bauer A, et al. Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature. 2002;415:141–7.
Tong AHY, Drees B, Nardelli G, Bader GD, Brannetti B, Castagnoli L, et al. A combined experimental and computational strategy to define protein interaction networks for peptide recognition modules. Science (80- ). 2002;295:321–4.
Torres PHM, Sodero ACR, Jofily P, Silva-Jr FP. Key topics in molecular docking for drug design. Int J Mol Sci. 2019;20:1–29.
Irungu BN, Orwa JA, Gruhnojic A, Fitzpatrick PA, Landberg G, Kimani F, et al. Constituents of the roots and leaves of Ekebergia capensis and their potential antiplasmodial and cytotoxic activities. Molecules. 2014;19:14235–46.
Liu C, Xu JF, He QM, Yu QS. Chemical components of the Liquidambar formosana H. Resin Youji Huaxue. 1991;11:508–10.
Mahato SB, Kundu AP. 13C NMR spectra of pentacyclic triterpenoids-a compilation and some salient features. Phytochemistry. 1994;37:1517–75.
Peng C, Bodenhausen G, Qiu S, Fong HHS, Farnsworth NR, Yuan S, et al. Computer-assisted structure elucidation: application of CISOC-SES to the resonance assignment and structure generation of betulinic acid. Magn Reson Chem. 1998;36:267–78.
Cheng F, Li W, Zhou Y, Shen J, Wu Z, Liu G, et al. AdmetSAR: a comprehensive source and free tool for assessment of chemical ADMET properties. J Chem Inf Model. 2012;52:3099–105.
Li YH, Yu CY, Li XX, 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.
The UniProt Consortium. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017;45:D158–69. https://doi.org/10.1093/nar/gkw1099.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45:D362–8. https://doi.org/10.1093/nar/gkw937.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504. https://doi.org/10.1101/gr.1239303.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The protein data Bank. Nucleic Acids Res. 2000;28:235–42. https://doi.org/10.1093/nar/28.1.235.
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091–3. https://doi.org/10.1093/bioinformatics/btp101.
Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. 2013;29:661–3. https://doi.org/10.1093/bioinformatics/btt019.
Varjosalo M, Keskitalo S, VanDrogen A, Nurkkala H, Vichalkovski A, Aebersold R, et al. The protein interaction landscape of the human CMGC kinase group. Cell Rep. 2013;3:1306–20. https://doi.org/10.1016/j.celrep.2013.03.027.
Butti R, Das S, Gunasekaran VP, Yadav AS, Kumar D, Kundu GC. Receptor tyrosine kinases (RTKs) in breast cancer: signaling, therapeutic implications and challenges. Mol Cancer. 2018;17:1–18.
Hubbard SR, Till JH. Protein tyrosine kinase structure and function. Annu Rev Biochem. 2000;69:373–98. https://doi.org/10.1146/annurev.biochem.69.1.373.
Lodge AJ, Anderson JJ, Gullick WJ, Haugk B, Leonard RCF, Angus B. Type 1 growth factor receptor expression in node positive breast cancer: adverse prognostic significance of c-erbB-4. J Clin Pathol. 2003;56:300–4.
Kim JY, Jung HH, Do IG, Bae SY, Lee SK, Kim SW, et al. Prognostic value of ERBB4 expression in patients with triple negative breast cancer. BMC Cancer. 2016;16:1–12.
Roskoski R. A historical overview of protein kinases and their targeted small molecule inhibitors. Pharmacol Res. 2015;100:1–23. https://doi.org/10.1016/j.phrs.2015.07.010.
Roskoski R. Small molecule inhibitors targeting the EGFR/ErbB family of protein-tyrosine kinases in human cancers. Pharmacol Res. November 2018;2019(139):395–411. https://doi.org/10.1016/j.phrs.2018.11.014.
Herrera-Abreu MT, Palafox M, Asghar U, Rivas MA, Cutts RJ, Garcia-Murillas I, et al. Early adaptation and acquired resistance to CDK4/6 inhibition in estrogen receptor-positive breast cancer. Cancer Res. 2016;76:2301–13.
Asghar US, Barr AR, Cutts R, Beaney M, Babina I, Sampath D, et al. Single-cell dynamics determines response to CDK4/6 inhibition in triple-negative breast cancer. Clin Cancer Res. 2017;23:5561–72.
Roskoski R. Cyclin-dependent protein serine/threonine kinase inhibitors as anticancer drugs. Pharmacol Res. 2019;139:471–88. https://doi.org/10.1016/j.phrs.2018.11.035.
Gee JMW, Robertson JFR, Ellis IO, Nicholson RI. Phosphorylation of ERK1/2 mitogen-activated protein kinase is associated with poor response to anti-hormonal therapy and decreased patient survival in clinical breast cancer. Int J Cancer. 2001;95:247–54. https://doi.org/10.1002/1097-0215(20010720)95:4<247::AID-IJC1042>3.0.CO;2-S.
Yeh YT, Hou MF, Chung YF, Chen YJ, Yang SF, Chen DC, et al. Decreased expression of phosphorylated JNK in breast infiltrating ductal carcinoma is associated with a better overall survival. Int J Cancer. 2006;118:2678–84.
Esteva FJ, Sahin AA, Smith TL, Yang Y, Pusztai L, Nahta R, et al. Prognostic significance of phosphorylated P38 mitogen-activated protein kinase and HER-2 expression in lymph node-positive breast carcinoma. Cancer. 2004;100:499–506.
Wu T, Wang X, Li J, Song X, Wang Y, Wang Y, et al. Identification of personalized chemoresistance genes in subtypes of basal-like breast cancer based on functional differences using pathway analysis. PLoS One. 2015;10:1–18.
Chakraborty G, Rangaswami H, Jain S, Kundu GC. Hypoxia regulates cross-talk between Syk and Lck leading to breast cancer progression and angiogenesis. J Biol Chem. 2006;281:11322–31.
Sibilski C, Mueller T, Kollipara L, Zahedi RP, Rapp UR, Rudel T, et al. Tyr728 in the kinase domain of the murine kinase suppressor of RAS 1 regulates binding and activation of the mitogen-activated protein kinase kinase. J Biol Chem. 2013;288:35237–52.
Narayanan R, Dalton JT. Androgen receptor: a complex therapeutic target for breast cancer. Cancers (Basel). 2016;8:1–17.
Yousefnia S, Momenzadeh S, Seyed Forootan F, Ghaedi K, Nasr Esfahani MH. The influence of peroxisome proliferator-activated receptor γ (PPARγ) ligands on cancer cell tumorigenicity. Gene. 2018;649:14–22. https://doi.org/10.1016/j.gene.2018.01.018.
Ghosh D, Griswold J, Erman M, Pangborn W. Structural basis for androgen specificity and oestrogen synthesis in human aromatase. Nature. 2009;457:219–23. https://doi.org/10.1038/nature07614.
Julien SG, Dubé N, Read M, Penney J, Paquet M, Han Y, et al. Protein tyrosine phosphatase 1B deficiency or inhibition delays ErbB2-induced mammary tumorigenesis and protects from lung metastasis. Nat Genet. 2007;39:338–46.
Yu M, Liu Z, Liu Y, Zhou X, Sun F, Liu Y, et al. PTP1B markedly promotes breast cancer progression and is regulated by miR-193a-3p. FEBS J. 2019;286:1136–53.
Zhang Y, Du Y. The development of protein tyrosine phosphatase1B inhibitors defined by binding sites in crystalline complexes. Future Med Chem. 2018;10:2345–67.
Li XQ, Wang LJ, Shi DY. The design strategy of selective PTP1B inhibitors over TCPTP. Bioorg Med Chem. 2016;24:3343–52. https://doi.org/10.1016/j.bmc.2016.06.035.
Butterworth S, Overduin M, Barr AJ. Targeting protein tyrosine phosphatase SHP2 for therapeutic intervention. Future Med Chem. 2014;6:1423–37. https://doi.org/10.4155/fmc.14.88.
Aceto N, Sausgruber N, Brinkhaus H, Gaidatzis D, Martiny-Baron G, Mazzarol G, et al. Tyrosine phosphatase SHP2 promotes breast cancer progression and maintains tumor-initiating cells via activation of key transcription factors and a positive feedback signaling loop. Nat Med. 2012;18:529–37.
Zhang X, He Y, Liu S, Yu Z, Jiang Z-X, Yang Z, et al. Salicylic acid based small molecule inhibitor for the oncogenic Src homology-2 domain containing protein tyrosine phosphatase-2 (SHP2). J Med Chem. 2010;53:2482–93. https://doi.org/10.1021/jm901645u.
Huang Q, Li F, Liu X, Li W, Shi W, Liu F, et al. Caspase 3-mediated repopulation of cancer cels after radiotherapy. Nat Med. 2012;17:860–6.
Gibson BA, Kraus WL. New insights into the molecular and cellular functions of poly (ADP-ribose) and PARPs. Nat Rev Mol Cell Biol. 2012;13:411–24.
Lee A, Djamgoz MBA. Triple negative breast cancer: emerging therapeutic modalities and novel combination therapies. Cancer Treat Rev. 2018;62:110–22. https://doi.org/10.1016/j.ctrv.2017.11.003.
Steffen JD, Brody JR, Armen RS, Pascal JM. Structural implications for selective targeting of PARPs. Front Oncol. 2013;3(December):301.
Jego G, Hazoumé A, Seigneuric R, Garrido C. Targeting heat shock proteins in cancer. Cancer Lett. 2013;332:275–85. https://doi.org/10.1016/j.canlet.2010.10.014.
Vazquez A, Bond EE, Levine AJ, Bond GL. The genetics of the p53 pathway, apoptosis and cancer therapy. Nat Rev Drug Discov. 2008;7:979–87.
Park JM, Kanaoka Y, Eguchi N, Aritake K, Grujic S, Materi AM, et al. Hematopoietic prostaglandin D synthase suppresses intestinal adenomas in Apc min/+ mice. Cancer Res. 2007;67:881–9. https://doi.org/10.1158/0008-5472.CAN-05-3767.
Fukuoka T, Yashiro M, Kinoshita H, Morisaki T, Hasegawa T, Hirakawa T, et al. Prostaglandin d synthase is a potential novel therapeutic agent for the treatment of gastric carcinomas expressing PPARγ. Int J Cancer. 2015;137:1235–44.
Zhang XJ, Han L, Ji YB, Fang GZ. Studies of betulonic acid on cell cycle and related protein expression on mice of bearing H22 tumor cells. China J Chin Mater Medica. 2008;33:1739–43.
Leong KH, Mahdzir MA, Din MFM, Awang K, Tanaka Y, Kulkeaw K, et al. Induction of intrinsic apoptosis in leukaemia stem cells and in vivo zebrafish model by betulonic acid isolated from Walsura pinnata Hassk (Meliaceae). Phytomedicine. 2017;26:11–21. https://doi.org/10.1016/j.phymed.2016.12.018.
Zhang H, Kang Y, Li N, Wang H, Bao Y, Li Y, et al. Triterpenoids from Liquidambar Fructus induced cell apoptosis via a PI3K-AKT related signal pathway in SMMC7721 cancer cells. Phytochemistry. 2020;171:112228. https://doi.org/10.1016/j.phytochem.2019.112228.
Knight ZA, Lin H, Shokat KM. Targeting the cancer kinome through polypharmacology. Nat Rev Cancer. 2010;10:130–7.
Huang S, Kauffman S. How to escape the cancer attractor: rationale and limitations of multi-target drugs. Semin Cancer Biol. 2013;23:270–8. https://doi.org/10.1016/j.semcancer.2013.06.003.
Raghavendra NM, Pingili D, Kadasi S, Mettu A, Prasad SVUM. Dual or multi-targeting inhibitors: the next generation anticancer agents. Eur J Med Chem. 2018;143:1277–300. https://doi.org/10.1016/j.ejmech.2017.10.021.
Dent R, Trudeau M, Pritchard KI, Hanna WM, Kahn HK, Sawka CA, et al. Triple-negative breast cancer: clinical features and patterns of recurrence. Clin Cancer Res. 2007;13:4429–34.
Jones FE. HER4 intracellular domain (4ICD) activity in the developing mammary gland and breast cancer. J Mammary Gland Biol Neoplasia. 2008;13:247–58.
Gassmann M, Casagranda F, Orioli D, Simon H, Lai C, Klein R, et al. Aberrant neural and cardiac development in mice lacking the ErbB4 neuregulin receptor. Nature. 1995;378:390–4. https://doi.org/10.1038/378390a0.
Meyn M III, Smithgall T. Small molecule inhibitors of Lck: the search for specificity within a kinase family. Mini Rev Med Chem. 2008;8:628–37.
Poehland BL, Carté BK, Francis TA, Hyland LJ, Allaudeen HS, Troupe N. In vitro antiviral activity of dammar resin triterpenoids. J Nat Prod. 1987;50:706–13. https://doi.org/10.1021/np50052a022.
Shirane N, Hashimoto Y, Ueda K, Takenaka H, Katoh K. Ring-a cleavage of 3-oxo-olean-12-en-28-oic acid by the fungus Chaetomium longirostre. Phytochemistry. 1996;43:99–104.
This work was supported by the National Natural Science Foundation of China (grant number 81703756) and the Natural Science Foundation of Liaoning Province (grant number 20170541028).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Separation process of the chemical compounds in LF.
Protein crystal structures and reference ligands used in the molecular docking experiment.
The NMR spectrum of LF03, LF04, LF05, LF08, LF09 and LF10.
Predicted targets of LF triterpenoids and therapeutic targets of breast cancer.
Molecular docking results of LF triterpenoids.
Protein families that the molecular docking targets belonged to.
Functionally grouped network view of the GO enrichment analysis result.
The framework and research strategy of the article.
About this article
Cite this article
Qian, P., Mu, XT., Su, B. et al. Identification of the anti-breast cancer targets of triterpenoids in Liquidambaris Fructus and the hints for its traditional applications. BMC Complement Med Ther 20, 369 (2020). https://doi.org/10.1186/s12906-020-03143-8
- Liquidambaris Fructus
- Breast cancer
- Network pharmacology
- Molecular docking