Next Article in Journal
Metabolic Profiling of Human Plasma and Urine, Targeting Tryptophan, Tyrosine and Branched Chain Amino Acid Pathways
Previous Article in Journal
Lipid Profile Characterization and Lipoprotein Comparison of Extracellular Vesicles from Human Plasma and Serum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Screening of Molecular Targets of Action of Atractylodin in Cholangiocarcinoma by Applying Proteomic and Metabolomic Approaches

by
Kanawut Kotawong
1,
Wanna Chaijaroenkul
1,2,
Sittiruk Roytrakul
3,
Narumon Phaonakrop
3 and
Kesara Na-Bangchang
1,2,*
1
Graduate Program in Bioclinical Sciences, Chulabhorn International College of Medicine, Thammasat University, Paholyothin Road, Klonglung, Pathumthani 12120, Thailand
2
Center of Excellence in Pharmacology and Molecular Biology of Malaria and Cholangiocarcinoma, Thammasat University, Paholyothin Road, Klonglung, Pathumthani 12120, Thailand
3
Center for Genetic Engineering and Biotechnology (BIOTEC), National Science and Technology Development Agency, Pathumthani 12120, Thailand
*
Author to whom correspondence should be addressed.
Metabolites 2019, 9(11), 260; https://doi.org/10.3390/metabo9110260
Submission received: 2 October 2019 / Revised: 29 October 2019 / Accepted: 31 October 2019 / Published: 1 November 2019

Abstract

:
Cholangiocarcinoma (CCA) is cancer of the bile duct and the highest incidence of CCA in the world is reported in Thailand. Our previous in vitro and in vivo studies identified Atractylodes lancea (Thunb) D.C. as a promising candidate for CCA treatment. The present study aimed to examine the molecular targets of action of atractylodin, the bioactive compound isolated from A. lancea, in CCA cell line by applying proteomic and metabolomic approaches. Intra- and extracellular proteins and metabolites were identified by LC-MS/MS following exposure of CL-6, the CCA cell line, to atractylodin for 24 and 48 h. Analysis of the protein functions and pathways involved was performed using a Venn diagram, PANTHER, and STITCH software. Analysis of the metabolite functions and pathways involved, including the correlation between proteins and metabolites identified was performed using MetaboAnalyst software. Results suggested the involvement of atractylodin in various cell biology processes. These include the cell cycle, apoptosis, DNA repair, immune response regulation, wound healing, blood vessel development, pyrimidine metabolism, the citrate cycle, purine metabolism, arginine and proline metabolism, glyoxylate and dicarboxylate metabolism, the pentose phosphate pathway, and fatty acid biosynthesis. Therefore, it was proposed that the action of atractylodin may involve the destruction of the DNA of cancer cells, leading to cell cycle arrest and cell apoptosis.

1. Introduction

Cancer is uncontrolled cell proliferation, and it is the leading cause of mortality worldwide. Liver and colon cancers, including cholangiocarcinoma (CCA) are cancer types with high mortality and incidental rates. CCA is cancer of the bile duct and the highest incidence in the world is reported in Thailand. It is classified into intrahepatic (iCCA), perihilar (pCCA), and distal (dCCA) CCA based on anatomical position. iCCA develops in liver parenchymal cells, and accounts for 20–25% all CCA cases [1]. The 3 and 5-years survival rates are 30% and 18%, respectively [2]. pCCA originates from the epithelial cells of the hepatic duct and spreads to common hepatic and cystic ducts. The five-year survival rates of pCCA are reported at 11%. dCCA is located in the area from the cystic duct to the ampulla of Vater. The incidence of this cancer type is about 20–30% of all CCA cases. The five-year survival rate of pCCA is about 11%. Gemcitabine, cisplatin, and 5-fluorouracil (5-FU) are currently used to treat CCA, but their clinical efficacy is unsatisfactory. Research and development of alternative medicines, particularly those from natural sources are challenging. Our previous in vitro and in vivo studies identified Atractylodes lancea (Thunb) D.C. (AL) as a promising candidate for CCA treatment [3]. AL has long been used as a traditional medicine in China, Japan, and Thailand. The two main bioactive compounds isolated from the AL rhizome are atractylodin (14%), and β-eudesmol (6%) [4]. It has been shown that the inhibitory action of atractylodin on cancer cell growth is through suppression of dendritic cell maturation [5], the delay of gastric emptying [6], and regulation of interleukin-6 in HMC-1 by blocking NPM-ALK activation [7]. The present study aimed to further examine other molecular targets of action of atractylodin in the CCA cell line by applying proteomic and metabolomic approaches.

2. Results

2.1. Proteomics

A total of 4366 and 4383 proteins were identified in the intra- and extracellular components of the treated (atractylodin) and untreated CL-6 cells after 24 h of incubation. The corresponding numbers of proteins after 48 h incubation were 4391 and 4348 proteins, respectively. The analysis of Venn diagrams suggested that 659 and 507 proteins from the intracellular components were only found in atracylodin-treated and untreated cells after 24 h incubation, respectively. The corresponding number of proteins in the extracellular components after 24 h incubation was 674 and 532 proteins, respectively (Figure 1a,b). These proteins were selected for gene ontology analysis by PANTHER software. Results showed that 585 proteins were involved in biological processes. Most (192 proteins) were involved in cellular processes (GO:0009987), followed by metabolic processes (GO:0008152) (140 proteins), and biological regulation (GO:0065007) (97 proteins) (Figure 2a). Also, 25 proteins were involved in transporter activity, i.e., ATP-binding cassette sub-family G member 8, CSC1-like protein 1, gamma-aminobutyric acid receptor subunit epsilon, and melanoma inhibitory activity protein 2 (Table 1). Analysis of protein-protein interactions of the 585 proteins by STITCH software suggested the involvement of these proteins in various biological processes including the cell cycle, DNA repair, apoptosis, regulation of immune response, wound healing, and blood vessel development (Figure 3a).
For the 48 h incubation, the analysis of the Venn diagrams suggested that 606 and 564 proteins from the intracellular components were only found in atracylodin-treated and untreated cells, respectively. The corresponding number of proteins in the extracellular components after 48 h incubation was 655 and 538 proteins, respectively (Figure 1c,d). Gene ontology analysis by PANTHER showed that 645 proteins were involved in biological processes. Most (192 proteins) were involved in the cellular process (GO:0009987), followed by the metabolic process (GO:0008152) (148 proteins), and biological regulation (GO:0065007) (102 proteins) (Figure 2b). Also, 34 proteins were involved in transporter activity, i.e., adenosine 3’-phospho 5’-phosphosulfate transporter 2, amiloride-sensitive sodium channel subunit gamma, and ATP-binding cassette sub-family A member 12 (Table 1). Analysis of protein-protein interactions of the 645 proteins by STITCH software (Figure 3b), suggested the involvement of these proteins in cell morphogenesis, cell adhesion, cell cycle, DNA repairing, apoptosis, regulation of immune response, and wound healing.

2.2. Metabolomics

Following exposure of CL-6 cell to atractylodin for 24 and 48 h, 289 metabolites from both ionization modes (positive and negative) were identified in the intra- and extracellular components by LC-MS/MS. The area under the peak of each metabolite was calculated and compared with the control (untreated cell). Metabolites with a change in the peak area of at least 2-fold were selected for further analysis. For the 24 h exposure, 21, 19, 11, and 16 metabolites were up-regulated intracellular, down-regulated intracellular, up-regulated extracellular, and down-regulated extracellular metabolites, respectively. The metabolite that showed the highest fold-change was the up-regulated intracellular metabolite methyl nicotinamide (8.57-fold), followed by the up-regulated intracellular metabolite inosine (8.05-fold), the down-regulated intracellular metabolite glutathione disulfide (6.84-fold), and the up-regulated intracellular metabolite phenylpyruvate (6.25-fold). Metabolites that showed the lowest fold-change were the up-regulated intracellular metabolite pyridoxine (2.01-fold) and the down-regulated intracellular metabolite UTP (2.01-fold) (Table 2). The top five metabolic pathways identified by MetaboAnalyst software included aspartate metabolism, pyrimidine metabolism, lactose synthesis, ammonia recycling, and glutamate metabolism (Figure 4a). For the 48 h exposure, 58, 11, 8, and 57 metabolites were up-regulated intracellular, down-regulated intracellular, up-regulated extracellular, and down-regulated extracellular metabolites, respectively. The metabolite that showed the highest fold-change was the down-regulated extracellular metabolite glutathione disulfide (65.53-fold), followed by the down-regulated intracellular metabolite glutathione disulfide (61.87-fold), the up-regulated intracellular metabolite taurodeoxycholic acid (22.41-fold), and the up-regulated intracellular metabolite 2-aminooctanoic acid (21.35-fold). The metabolite that showed the lowest fold-change was the up-regulated intracellular metabolite pantothenate (2.00-fold) (Table 3). The top five metabolic pathways identified by MetaboAnalyst software included the citric acid cycle, pyrimidine metabolism, Warburg effect, glycine and serine metabolism, and aspartate metabolism (Figure 4b).

2.3. Co-Analysis of Proteomics and Metabolomics

The proteomics and metabolomics information were co-analyzed for a 24 h and 48 h exposure period. Proteins from the intra- and extracellular components that were found only in atractylodin-treated and untreated CL-6 cells, together with the up-regulated and down-regulated metabolites from both intra- and extracellular components following 24 and 48 h incubation were selected for analysis. For the 24 h exposure, results from MetaboAnalyst analysis revealed four proteins that are associated with six metabolites (Figure 5a). The pathways associated with these interactions are pyrimidine metabolism, purine metabolism, and the pentose phosphate pathway. The biological processes that are linked to these pathways are blood coagulation and wound healing (Figure 5a). For the 48 h exposure, results from MetaboAnalyst analysis revealed 29 proteins that are associated with 20 metabolites (Figure 5b). The pathways associated with these interactions are pyrimidine metabolism, the citrate cycle (TCA cycle), taurine and hypotaurine metabolism, biotin metabolism, purine metabolism, arginine and proline metabolism, glyoxylate and dicarboxylate metabolism, the pentose phosphate pathway, fatty acid biosynthesis, and primary bile acid biosynthesis. The biological processes that are linked to these pathways are blood coagulation, inflammatory response, regulation of the immune system process, regulation of cell proliferation, Ras protein signal transduction, and DNA repairing (Figure 5b).

3. Discussion

In this study, proteomic and metabolomic approaches were applied for the identification of the possible molecular targets of action of atractylodin in cholangiocarcinoma. The proteomics results showed that most proteins from the intra- and extracellular components after CL-6 exposure to atractylodin for 24 and 48 h were involved in cellular (GO:0009987) and metabolic (GO:0008152) processes. The proteins with transporter activity were low in number. Transporter proteins play an essential role in cancer chemotherapy and drug resistance. Modulation of drug transport through decreasing influx and increasing efflux has been shown to be one major mechanism of anticancer drug resistance. Transporter proteins are classified into two types, i.e., the ATP binding cassette (ABC) family and the solute carrier (SLC) family [8]. These transporters function in various pharmacokinetic processes including absorption and excretion of xenobiotics and endogenous molecules. Two transporter proteins identified in this study belong to the ATP-binding cassette sub-family A member 12 (ABCA12), and the ATP-binding cassette sub-family G member 8(ABCG8). The ABCA12 protein is a member of the ATP-binding cassette (ABC) family and it functions as lipid transporters. Up-regulation of ABCA12 was found to be associated with anticancer drug resistance in several types of cancer such as breast [9], ovarian [10], and colorectal cancer [11]. The ATP-binding cassette sub-family G member 8 (ABCG8) belongs to the ATP-binding cassette (ABC) family. This protein is located in enterocyte membranes and plays an essential role in controlling the transportation of dietary cholesterols into and out of the cells. Down-regulation of ABCG8 was observed in breast carcinoma after chemotherapy [9]. Besides, the suppression of ABCG8 induced the apoptosis of liver cancer by inhibiting doxorubicin efflux [12]. Other transporter proteins such as the calcium-binding mitochondrial carrier protein SCaMC-3 (CAD55563) are commonly found in the brain, liver, and small intestine at high density [13]. The function of this protein includes the transportation of adenine nucleotides into mitochondria. Increasing the expression of SCaMC-3 activates the influx of calcium to mitochondria, and thus, induction of cell apoptosis [14]. The transient receptor potential cation channel subfamily V member 6 protein (TRPV6), is a calcium selective cation channel that plays a critical role in calcium uptake in epithelial tissues in epididymis, placenta, and prostate cancer [15,16]. Free intracellular Ca2+ activates muscle contraction and cell cycle regulation. In gastric cancer, the increase in free Ca2+ in mitochondria via stimulation of TRPV6 protein was shown to result in cancer cell apoptosis [17].
Based on the results of the protein-protein interaction study, atractylodin action was associated with DNA repairing, the cell cycle, and apoptosis. DNA repairing involves several proteins, including breast cancer 2 (BRCA2), BRCA1 associated RING domain 1(BARD1), and proliferating cell nuclear antigen (PCNA). The BRCA2 protein belongs to the tumor suppressor gene family, which is located on chromosome 13. It functions in repairing DNA damage induced by chemical exposure. The PCNA protein is found in yeast, plant, and animal cells. The protein regulates cell division, DNA replication, and cell death through MyD118 and p21 [18]. The induction of cell apoptosis was shown to be associated with a decrease in PCNA expression [19]. BARD1 protein is synthesized in the cell cycle process at the S phase. The function of BARD1 involves DNA damage and repairing, and transcriptional regulation, including apoptosis induction. This protein activates cell apoptosis via p53-phosphorylation [20]. Altogether, results suggest that the mechanism of cytotoxicity of atractylodin is through induction of DNA damage similar to that observed with cisplatin, gemcitabine, and oxaliplatin [21]. This DNA damage results in the activation of cell cycle arrest and apoptosis. The cell cycle process is involved in cell division and is controlled by cell cycle checkpoints at G1, G2, and M checkpoints. Major proteins that are involved in the cell cycle process include cyclin A, cyclin B, cyclin D, cyclin E, cdk2, cdk4, and cdk6. Abnormal functioning of these proteins lead to cell cycle arrest. Atractylodin has previously been shown to induce cholangiocarcinoma cell cycle arrest at the G1 phase [22]. The apoptotic action of atractylodin has also been shown to be both concentration- and time-dependent [22]. Apoptosis, the cell process associated with cell death, involves two main mechanisms, i.e., intrinsic and extrinsic pathways. In the intrinsic pathway, activation of cell apoptosis is through DNA damage. In the extrinsic pathway, on the other hand, activation of cell apoptosis is through the death cell receptor.
The metabolomic study revealed the link between atractylodin action and several metabolic pathways, including aspartate metabolism, pyrimidine metabolism, the citric acid cycle, and ammonia recycling. The primary metabolites identified were glutathione disulfide (GSSG: 65.53-fold change), and taurodeoxycholic acid (TDCA: 22.41-fold change). GSSG consists of two molecules of glutathione linked by a disulfide bond. It plays a role in maintaining redox homeostasis [23]. NOV-002, the glutathione disulfide mimetic, was shown to suppress tumor cell invasion and metastasis via ErbB2, PI3K, Akt, and RhoA [24]. TDCA consists of the bile acid taurine (found in the liver), which is conjugated with deoxycholic acid. This metabolite activates IL-8 gene expression via RelA phosphorylation [25] and induction of HepG2 cell apoptosis [26]. For a better understanding of the atractylodin action, data from both proteomics together with metabolomics were further analyzed by MetaboAnalyst software. Results suggested that the key action of atractylodin is wound healing and blood coagulation. Wound healing is the repair process of skin and tissue following injury. It is an important process that leads to chronic fibrosis, and cancer development/progression. Prolonged and exacerbated healing could trigger chronic fibrosis and tumorigenesis [27]. The blood coagulation process is associated with tumor progression and metastasis [28]. Suppression of this process was previously demonstrated with atractylodin and the colorectal cancer drug avastin [29]. Atractylodin was shown to inhibit wound healing of the cholangiocarcinoma cell CL-6 after 24 and 48 h of exposure [30]. Besides, it also induced the suppression of HO-1 production, STAT1/3 phosphorylation, and NF-κB proteins [30]. The association between atractylodin action and blood coagulation remains to be clarified. The process is induced by two pathways, i.e., intrinsic and extrinsic pathways. The intrinsic pathway is activated by damaging the vessel wall. The factors that control this process are the Hageman factor (XII), plasma thromboplastin (XI), Christmas factor (IX), and Stuart factor (X) [31]. Abnormality of these factors causes hemophilia B and C [32]. The extrinsic pathway is activated by tissue damage outside the blood vessel, which is mainly controlled by the stable factor (VII) [31]. Abnormality of factor VII causes hemophilia A [32].

4. Materials and Methods

4.1. Cell Culture

The human CCA cell line (CL-6) is the intrahepatic cholangiocarcinoma that was kindly provided by Associate Professor Adisak Wongkajornsilp, Department of Pharmacology, Faculty of Medicine (Siriraj Hospital), Mahidol University. The cell was cultured in RPMI 1640 medium supplemented with 10% (v/v) heated fetal bovine serum (FBS), and 100 IU/mL anti-biotic-antimycotic solution, and maintained at 37 °C, 5% CO2 atmosphere, and 95% humidity (HERA CELL 150i, Thermo Scientific, USA).

4.2. Sample Preparation

The CL-6 (106 cells) was seeded into a 25 cm2 culture flask and incubated overnight at 37 °C, 5% CO2 atmosphere. The cell was incubated with atractylodin (40 µg/mL) for 24 h and 48 h. After incubation, proteins and metabolites from intracellular (treated cell) and extracellular components (culture media) were extracted and analyzed by LC-MS/MS.

4.3. Proteomics

4.3.1. Protein Extraction

Following exposure of CL-6 cells to atractylodin for 24 h and 48 h, the intra- and extracellular proteins were extracted. For intracellular proteins, the cell was washed with cold PBS and broken using SDS (0.5%). Cell supernatant was separated through centrifugation at 800× g (4 °C) for 10 min and stored at −80 °C until use. For extracellular proteins, proteins in the culture media were precipitated by mixing with pre-cold acetone (−20 °C) and incubated overnight at −20 °C. After incubation, the precipitated proteins were separated through centrifugation at 13,000× g for 15 min. Protein pellets were dissolved with 0.5% SDS and stored at −80 °C until use. Concentrations of both intra- and extracellular proteins were measured by the Lowry method [33].

4.3.2. SDS-PAGE Analysis and In-Gel Digestion by Trypsin

The proteins (30 µg/sample) extracted from both intra- and extracellular were separated using 12% SDS-PAGE (20 mA constant current for 80 min) and stained using Coomassie Brilliant Blue R. The stained proteins in the gel were cut into 19 small fragments per each sample. The gel fragments were incubated with dithiothreitol (DTT) and iodoacetamide (IAA) before being digested with trypsin and extracted with acetonitrile [34].

4.3.3. LC-MS/MS Analysis

The digested proteins extracted from the gel fragments were dissolved with 0.1% formic acid and injected into LC-MS/MS for protein identification. Briefly, the digested proteins were injected into a μ-pre-column (Monolithic Trap Column, 200 μm i.d. × 5 mm) coupled with the Ultimate 3000 LC system (Dionex), and ESI-Ion Trap MS (HCT ultra PTM Discovery System, BrukerDaltonik) with electrospray at a flow rate of 20 μL/min. Proteins were separated on a nano-column (Monolithic Nano Column, 100 μm i.d.× 5 cm) (Thermo Fisher Scientific, Amsterdam, Netherland)with a solvent gradient mobile phase consisting of solvent A (0.1% formic acid) and solvent B (50% water, 50% acetonitrile, 0.1% formic acid) running at a flow rate of 1 µL/min for 20 min.

4.3.4. Protein Identification and Analysis

The LC-MS/MS data were analyzed using DeCyderTM (Amersham Bioscience AB, Uppsala, Sweden) and MASCOT (http://www.matrixscience.com) programs for protein identification. The PANTHER software (http://www.pantherdb.org) was selected for studying the gene ontology, and the STITCH software was used to study the protein-protein interactions (http://stitch.embl.de).

4.4. Metabolomics

4.4.1. Metabolite Extraction

CL-6 cell was cultured in a 25 cm2 flask (106 cells/flask) and exposed to atractylodin (40 µg/mL) for 24 h and 48 h. Quenching reagent (cooled methanol) was immediately added (80% v/v of methanol solution) to halt any biochemical reactions. For the metabolite extraction, the samples were divided into two parts, i.e., intracellular metabolites, and extracellular metabolites. For intracellular metabolites, complete media was removed from CL-6 cells after the addition of a quenching reagent. The cell was immediately washed with cold PBS, and transferred to a new tube (1.5 mL) using a cell scraper. Following centrifugation (14,000× g, at 4–8 °C for 5 min), the cell supernatant was separated and dried using SpeedVac/lyophilize. These intracellular metabolites were stored at −80 °C until analysis. For extracellular metabolites, the metabolites in media were collected by centrifugation (14,000× g for 10 min, 4 °C) and transferred to a new tube (1.5 mL). The extracellular metabolites were dried by SpeedVac/lyophilize and stored at −80 °C until use.

4.4.2. Metabolite Identification

The dried metabolites from both the intra- and extracellular components were reconstituted with LC/MS grade water (500 µL), transferred to a 1.5 mL glass tube, and filtered through the phenix-RC filter (0.2 µm pore size). The metabolites were separated using HPLC (Agilent Technologies, 1200 Infinity Series) and identified by QTRAP® 5500 (AB SCIEX) [35]. The system consisted of security universal HPLC guard cartridge (Phenomenex, CA, USA), amide XBridge HPLC column (3.5 μm × 4.6 mm i.d. × 100 mm length (Waters, MA, USA). The mobile phase consisted of 50% buffer A (95% water, 5% acetonitrile, 20 mM ammonium hydroxide, and 20 mM ammonium acetate adjusted to pH 9.0) and 50% of buffer B (100% acetonitrile). The flow rate was 400 µL/minute with a run time of 23 min. For MS/MS, the spray voltage was set at 3200 V. Nitrogen and argon were used as auxiliary and collision gas, respectively. The scan time for each single-reaction monitoring (SRM) event transition was 0.1 s with a scan width of 1 m/z.

4.4.3. Metabolomics Analysis

The instrument control, data acquisition, and data analysis were achieved using XcalibarTM software (Thermo Scientific, MA, USA). The up- and down-regulated metabolites were selected for the identification of possible signally pathways using the free online web-based MetaboAnalyst 4.0 (http://old.metaboanalyst.ca/MetaboAnalyst/faces/ModuleView.xhtmL). Moreover, MetaboAnalyst 4.0 was used to determine the correlation between metabolites and proteins in this study.

5. Conclusions

In conclusion, by applying proteomic and metabolomic approaches, possible molecular targets of atractylodin in cholangiocarcinoma were shown to be associated with various cell biological processes including the cell cycle, apoptosis, DNA repair, immune response regulation, wound healing, blood vessel development, cell morphogenesis, cell proliferation, Ras protein signal transduction, pyrimidine metabolism, the citrate cycle (TCA cycle), purine metabolism, arginine and proline metabolism, glyoxylate and dicarboxylate metabolism, pentose phosphate pathway, fatty acid biosynthesis, and primary bile acid biosynthesis. Therefore, it was proposed that the action of atractylodin may involve the destruction of DNA of cancer cells, leading to cell cycle arrest and cell apoptosis. The exact molecular targets and their signaling pathways remain to be clarified.

Author Contributions

Writing—Original draft, K.K.; Formal analysis, K.K., W.C., N.P., and S.R.; Supervision, K.N.-B.; Writing—review & editing, K.N.-B.

Funding

This research was funded by Center of Excellence in Pharmacology and Molecular Biology of Malaria and Cholangiocarcinoma of Thammasat University (1/2556), National Research Council of Thailand (NRCT) (37/2562), and Thammasat University Research Fund under the TU Research Scholar, Contract No. 02/2561 and Thammasat University Ph.D. program scholarship.

Acknowledgments

The author would like to thank the Drug Discovery and Development Center, Thammasat University and Center for Genetic Engineering and Biotechnology (BIOTEC), National Science and Technology Development Agency for facilities support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khan, S.A.; Davidson, B.R.; Goldin, R.; Pereira, S.P.; Rosenberg, W.M.; Taylor-Robinson, S.D.; Thillainayagam, A.V.; Thomas, H.C.; Thursz, M.R.; Wasan, H.; et al. Guidelines for the diagnosis and treatment of cholangiocarcinoma: Consensus document. Gut 2002, 51, VI1–VI9. [Google Scholar] [CrossRef] [PubMed]
  2. Nathan, H.; Aloia, T.A.; Vauthey, J.N.; Abdalla, E.K.; Zhu, A.X.; Schulick, R.D.; Choti, M.A.; Pawlik, T.M. A proposed staging system for intrahepatic cholangiocarcinoma. Ann. Surg. Oncol. 2009, 16, 14–22. [Google Scholar] [CrossRef] [PubMed]
  3. Mahavorasirikul, W.; Viyanant, V.; Chaijaroenkul, W.; Itharat, A.; Na-Bangchang, K. Cytotoxic activity of Thai medicinal plants against human cholangiocarcinoma, laryngeal and hepatocarcinoma cells in vitro. BMC Complement. Altern. Med. 2010, 10, 55. [Google Scholar] [CrossRef] [PubMed]
  4. Na-Bangchang, K.; Karbwang, J. Traditional herbal medicine for the control of tropical diseases. Trop. Med. Health 2014, 42, 3–13. [Google Scholar] [CrossRef] [PubMed]
  5. Chuang, C.H.; Cheng, Y.C.; Lin, S.C.; Lehman, C.W.; Wang, S.P.; Chen, D.Y.; Tsai, S.W.; Lin, C.C. Atractylodin Suppresses Dendritic Cell Maturation and Ameliorates Collagen-Induced Arthritis in a Mouse Model. J. Agric. Food Chem. 2019, 67, 6773–6784. [Google Scholar] [CrossRef] [PubMed]
  6. Nakai, Y.; Kido, T.; Hashimoto, K.; Kase, Y.; Sakakibara, I.; Higuchi, M.; Sasaki, H. Effect of the rhizomes of Atractylodes lancea and its constituents on the delay of gastric emptying. J. Ethnopharmacol. 2003, 84, 51–55. [Google Scholar] [CrossRef]
  7. Chae, H.S.; Kim, Y.M.; Chin, Y.W. Atractylodin Inhibits Interleukin-6 by Blocking NPM-ALK Activation and MAPKs in HMC-1. Molecules 2016, 21, 1169. [Google Scholar] [CrossRef]
  8. Liang, Y.; Li, S.; Chen, L. The physiological role of drug transporters. Protein Cell 2015, 6, 334–350. [Google Scholar] [CrossRef] [Green Version]
  9. Hlavac, V.; Brynychova, V.; Vaclavikova, R.; Ehrlichova, M.; Vrana, D.; Pecha, V.; Kozevnikovova, R.; Trnkova, M.; Gatek, J.; Kopperova, D.; et al. The expression profile of ATP-binding cassette transporter genes in breast carcinoma. Pharmacogenomics 2013, 14, 515–529. [Google Scholar] [CrossRef]
  10. Elsnerova, K.; Mohelnikova-Duchonova, B.; Cerovska, E.; Ehrlichova, M.; Gut, I.; Rob, L.; Skapa, P.; Hruda, M.; Bartakova, A.; Bouda, J.; et al. Gene expression of membrane transporters: Importance for prognosis and progression of ovarian carcinoma. Oncol. Rep. 2016, 35, 2159–2170. [Google Scholar] [CrossRef]
  11. Hlavata, I.; Mohelnikova-Duchonova, B.; Vaclavikova, R.; Liska, V.; Pitule, P.; Novak, P.; Bruha, J.; Vycital, O.; Holubec, L.; Treska, V.; et al. The role of ABC transporters in progression and clinical outcome of colorectal cancer. Mutagenesis 2012, 27, 187–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Stokes, E.; Shuang, T.; Zhang, Y.; Pei, Y.; Fu, M.; Guo, B.; Parissenti, A.; Wu, L.; Wang, R.; Yang, G. Efflux inhibition by H2S confers sensitivity to doxorubicin-induced cell death in liver cancer cells. Life Sci. 2018, 213, 116–125. [Google Scholar] [CrossRef] [PubMed]
  13. Bassi, M.T.; Manzoni, M.; Bresciani, R.; Pizzo, M.T.; Della Monica, A.; Barlati, S.; Monti, E.; Borsani, G. Cellular expression and alternative splicing of SLC25A23, a member of the mitochondrial Ca2+-dependent solute carrier gene family. Gene 2005, 345, 173–182. [Google Scholar] [CrossRef] [PubMed]
  14. Hoffman, N.E.; Chandramoorthy, H.C.; Shanmughapriya, S.; Zhang, X.Q.; Vallem, S.; Doonan, P.J.; Malliankaraman, K.; Guo, S.; Rajan, S.; Elrod, J.W.; et al. SLC25A23 augments mitochondrial Ca(2)(+) uptake, interacts with MCU, and induces oxidative stress-mediated cell death. Mol. Biol. Cell 2014, 25, 936–947. [Google Scholar] [CrossRef]
  15. Wissenbach, U.; Niemeyer, B.A.; Fixemer, T.; Schneidewind, A.; Trost, C.; Cavalie, A.; Reus, K.; Meese, E.; Bonkhoff, H.; Flockerzi, V. Expression of CaT-like, a novel calcium-selective channel, correlates with the malignancy of prostate cancer. J. Biol. Chem. 2001, 276, 19461–19468. [Google Scholar] [CrossRef]
  16. Suzuki, Y.; Chitayat, D.; Sawada, H.; Deardorff, M.A.; McLaughlin, H.M.; Begtrup, A.; Millar, K.; Harrington, J.; Chong, K.; Roifman, M.; et al. TRPV6 Variants Interfere with Maternal-Fetal Calcium Transport through the Placenta and Cause Transient Neonatal Hyperparathyroidism. Am. J. Hum. Genet. 2018, 102, 1104–1114. [Google Scholar] [CrossRef] [Green Version]
  17. Chow, J.; Norng, M.; Zhang, J.; Chai, J. TRPV6 mediates capsaicin-induced apoptosis in gastric cancer cells—Mechanisms behind a possible new “hot” cancer treatment. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2007, 1773, 565–576. [Google Scholar] [CrossRef]
  18. Strzalka, W.; Ziemienowicz, A. Proliferating cell nuclear antigen (PCNA): A key factor in DNA replication and cell cycle regulation. Ann. Bot. 2011, 107, 1127–1140. [Google Scholar] [CrossRef]
  19. Paunesku, T.; Mittal, S.; Protic, M.; Oryhon, J.; Korolev, S.V.; Joachimiak, A.; Woloschak, G.E. Proliferating cell nuclear antigen (PCNA): Ringmaster of the genome. Int. J. Radiat. Biol. 2001, 77, 1007–1021. [Google Scholar] [CrossRef]
  20. Feki, A.; Jefford, C.E.; Berardi, P.; Wu, J.-Y.; Cartier, L.; Krause, K.-H.; Irminger-Finger, I. BARD1 induces apoptosis by catalysing phosphorylation of p53 by DNA-damage response kinase. Oncogene 2005, 24, 3726–3736. [Google Scholar] [CrossRef]
  21. Woods, D.; Turchi, J.J. Chemotherapy induced DNA damage response: Convergence of drugs and pathways. Cancer Biol. Ther. 2013, 14, 379–389. [Google Scholar] [CrossRef] [PubMed]
  22. Kotawong, K.; Chaijaroenkul, W.; Muhamad, P.; Na-Bangchang, K. Cytotoxic activities and effects of atractylodin and beta-eudesmol on the cell cycle arrest and apoptosis on cholangiocarcinoma cell line. J. Pharmacol. Sci. 2018, 136, 51–56. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, Y.; Hyde, A.S.; Simpson, M.A.; Barycki, J.J. Emerging regulatory paradigms in glutathione metabolism. Adv. Cancer Res. 2014, 122, 69–101. [Google Scholar] [PubMed]
  24. Gumireddy, K.; Li, A.; Cao, L.; Yan, J.; Liu, L.; Xu, X.; Pazoles, C.; Huang, Q. NOV-002, A Glutathione Disulfide Mimetic, Suppresses Tumor Cell Invasion and Metastasis. J. Carcinog. Mutagen. 2013, 2013, s7-002. [Google Scholar] [PubMed]
  25. Muhlbauer, M.; Allard, B.; Bosserhoff, A.K.; Kiessling, S.; Herfarth, H.; Rogler, G.; Scholmerich, J.; Jobin, C.; Hellerbrand, C. Differential effects of deoxycholic acid and taurodeoxycholic acid on NF-kappa B signal transduction and IL-8 gene expression in colonic epithelial cells. Am. J. Physiol. Gastrointest. Liver Physiol. 2004, 286, G1000–G1008. [Google Scholar] [CrossRef]
  26. Xie, Q.; Li, G.M.; Zhou, X.Q.; Liao, D.; Yu, H.; Guo, Q. Effect of Tauroursodeoxycholic acid on cytochrome C-mediated apoptosis in HepG2 cells. Zhonghua Gan Zang Bing Za Zhi 2003, 11, 298–301. [Google Scholar]
  27. Rybinski, B.; Franco-Barraza, J.; Cukierman, E. The wound healing, chronic fibrosis, and cancer progression triad. Physiol. Genom. 2014, 46, 223–244. [Google Scholar] [CrossRef]
  28. Lima, L.G.; Monteiro, R.Q. Activation of blood coagulation in cancer: Implications for tumour progression. Biosci. Rep. 2013, 33, e00064. [Google Scholar] [CrossRef]
  29. Scappaticci, F.A.; Fehrenbacher, L.; Cartwright, T.; Hainsworth, J.D.; Heim, W.; Berlin, J.; Kabbinavar, F.; Novotny, W.; Sarkar, S.; Hurwitz, H. Surgical wound healing complications in metastatic colorectal cancer patients treated with bevacizumab. J. Surg. Oncol. 2005, 91, 173–180. [Google Scholar] [CrossRef]
  30. Mathema, V.B.; Chaijaroenkul, W.; Na-Bangchang, K. Cytotoxic activity and molecular targets of atractylodin in cholangiocarcinoma cells. J. Pharm. Pharmacol. 2019, 71, 185–195. [Google Scholar] [CrossRef]
  31. Palta, S.; Saroa, R.; Palta, A. Overview of the coagulation system. Indian J. Anaesth. 2014, 58, 515–523. [Google Scholar] [CrossRef] [PubMed]
  32. Castaman, G.; Matino, D. Hemophilia A and B: Molecular and clinical similarities and differences. Haematologica 2019, 104, 1702–1709. [Google Scholar] [CrossRef] [PubMed]
  33. Lowry, O.H.; Rosebrough, N.J.; Farr, A.L.; Randall, R.J. Protein measurement with the Folin phenol reagent. J. Biol. Chem. 1951, 193, 265–275. [Google Scholar] [PubMed]
  34. Kotawong, K.; Thitapakorn, V.; Roytrakul, S.; Phaonakrop, N.; Viyanant, V.; Na-Bangchang, K. Plasma phosphoproteome and differential plasma phosphoproteins with opisthorchis viverrini-related cholangiocarcinoma. Asian Pac. J. Cancer Prev. 2015, 16, 1011–1018. [Google Scholar] [CrossRef] [PubMed]
  35. Yuan, M.; Breitkopf, S.B.; Yang, X.; Asara, J.M. A positive/negative ion-switching, targeted mass spectrometry-based metabolomics platform for bodily fluids, cells, and fresh and fixed tissue. Nat. Protoc. 2012, 7, 872–881. [Google Scholar] [CrossRef]
Figure 1. Venn diagrams of intracellular proteins at 24 h (a), extracellular proteins at 24 h (b), intracellular proteins at 48 h (c), and extracellular proteins at 48 h of CL-6 cell incubation (d).
Figure 1. Venn diagrams of intracellular proteins at 24 h (a), extracellular proteins at 24 h (b), intracellular proteins at 48 h (c), and extracellular proteins at 48 h of CL-6 cell incubation (d).
Metabolites 09 00260 g001
Figure 2. Protein classification based on the biological process of the identified proteins at 24 h (a) and 48 h (b) of CL-6 cell incubation.
Figure 2. Protein classification based on the biological process of the identified proteins at 24 h (a) and 48 h (b) of CL-6 cell incubation.
Metabolites 09 00260 g002
Figure 3. Protein-protein interactions of the identified proteins at 24 h (a) and 48 h (b) of CL-6 cell incubation by STITCH analysis.
Figure 3. Protein-protein interactions of the identified proteins at 24 h (a) and 48 h (b) of CL-6 cell incubation by STITCH analysis.
Metabolites 09 00260 g003
Figure 4. The analysis of metabolites for 24 h (a) and 48 h (b) of atractylodin exposure by MetaboAnalyst software.
Figure 4. The analysis of metabolites for 24 h (a) and 48 h (b) of atractylodin exposure by MetaboAnalyst software.
Metabolites 09 00260 g004
Figure 5. The co-analysis by MetaboAnalyst software of the proteins and metabolites detected at 24 h (a) and 48 h (b) following atractylodin exposure.
Figure 5. The co-analysis by MetaboAnalyst software of the proteins and metabolites detected at 24 h (a) and 48 h (b) following atractylodin exposure.
Metabolites 09 00260 g005
Table 1. Transporter proteins identified at 24 h and 48 h of CL-6 cell incubation.
Table 1. Transporter proteins identified at 24 h and 48 h of CL-6 cell incubation.
Transporter Proteins Identified at 24 h Incubation
Protein NamesProtein NamesProtein Names
ATP-binding cassette sub-family A member 12Polycystic kidney disease protein 1-like 2 Retinal-specific ATP-binding cassette transporter
ATP-binding cassette sub-family G member 8Polycystic kidney disease protein 1-like 3Sodium/hydrogen exchanger 9B2
CSC1-like protein 1Potassium voltage-gated channel subfamily A member 4Transient receptor potential cation channel subfamily V member 4
Gamma-aminobutyric acid receptor subunit epsilonPotassium voltage-gated channel subfamily H member 4Transient receptor potential cation channel subfamily V member 6
Melanoma inhibitory activity protein 2Potassium voltage-gated channel subfamily H member 8Voltage-dependent L-type calcium channel subunit beta-2
Mitochondrial inner membrane protein OXA1LPotassium voltage-gated channel subfamily KQT member 3V-type proton ATPase 116 kDa subunit α isoform 2
Mitochondrial pyruvate carrier 1Potassium voltage-gated channel subfamily S member 1Protein shisa-7
Neuronal acetylcholine receptor subunit alpha-10Potassium-transporting ATPase subunit beta
NTF2-related export protein 1Probable phospholipid-transporting ATPase VD
Transporter Proteins Identified at 48 h Incubation
Protein NamesProtein NamesProtein Names
Amiloride-sensitive sodium channel subunit gammaMitochondrial pyruvate carrier 1Protein shisa-7
ATP-binding cassette sub-family A member 12Neuronal acetylcholine receptor subunit alpha-10Sodium/bile acid cotransporter
ATP-binding cassette sub-family G member 8Neuronal acetylcholine receptor subunit alpha-7Sodium/hydrogen exchanger 9B2
Calcium permeable stress-gated cation channel 1Neuronal acetylcholine receptor subunit beta-3Synaptic vesicular amine transporter
Calcium-activated chloride channel regulator 4NTF2-related export protein 1Transient receptor potential cation channel subfamily V member 6
Calcium-binding mitochondrial carrier protein SCaMC-3Phospholipid-transporting ATPase IBTransmembrane channel-like protein 2
Ceramide-1-phosphate transfer proteinPolycystic kidney disease protein 1-like 2Vacuolar protein sorting-associated protein 26B
Excitatory amino acid transporter 5Polycystic kidney disease protein 1-like 3Vesicular acetylcholine transporter
Importin subunit alpha-4Potassium voltage-gated channel subfamily H member 4Zinc transporter ZIP6
Melanoma inhibitory activity protein 2Potassium voltage-gated channel subfamily S member 1
Mitochondrial glutamate carrier 2Potassium-transporting ATPase subunit beta
Table 2. The fold-change of metabolites identified in intra- and extracellular components at 24 h of atractylodin exposure.
Table 2. The fold-change of metabolites identified in intra- and extracellular components at 24 h of atractylodin exposure.
Up-Regulation of Intracellular Metabolites at 24 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
Inosine8.05Threonine3.37Taurodeoxycholic acid2.24
Phenylpyruvate6.25S-ribosyl-L-homocysteine3.31Phenylalanine2.22
6-phospho-D-glucono-1,5-lactone5.86Indole2.81Kynurenine2.16
Cytidine5.66Aconitate2.73Folate2.14
Choline5.60Creatine2.70Tyrosine2.13
UDP-D-glucose5.32L-arginino-succinate2.522-oxo-4-methylthiobutanoate2.12
Serine3.85D-gluconate2.24Pyridoxine2.01
Down-Regulation of Intracellular Metabolites at 24 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
Glutathione disulfide6.84AMP2.71Cholesteryl sulfate2.16
Lactate5.11Thiamine2.616-phospho-D-gluconate2.08
Methionine sulfoxide4.43Cystine2.53Glucosamine2.06
Glutamine4.21Aspartate2.43Glucose-1-phosphate2.04
Citrulline3.68Cellobiose2.41UTP2.01
4-aminobutyrate2.86Asparagine2.34
N-acetyl-glutamate2.83dAMP2.17
Up-Regulation of Extracellular Metabolites at 24 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
Methylnicotinamide8.57Glycerate3.00Orotidine-5-phosphate2.12
Adenine5.812,3-dihydroxybenzoic acid2.73Serine2.07
Pyroglutamic acid3.87Hydroxyisocaproic acid2.25Adenosine2.04
Cytidine3.61Dihydroorotate2.16
Down-Regulation of Extracellular Metabolites at 24 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
CDP4.83N-carbamoyl-L-aspartate2.79Acetylphosphate2.44
Glutathione disulfide3.96Orotate2.71Citraconic acid2.20
N6-Acetyl-L-lysine3.392-ketohaxanoic acid2.66Glucose-1-phosphate2.12
Citrate3.112-dehydro-D-gluconate2.61Taurodeoxycholic acid2.11
Phenylpyruvate3.05Shikimate-3-phosphate2.60
2-oxo-4-methylthiobutanoate2.84D-gluconate2.48
Table 3. The fold-change of metabolites identified in intra- and extracellular components at 48 h of atractylodin exposure.
Table 3. The fold-change of metabolites identified in intra- and extracellular components at 48 h of atractylodin exposure.
Up-Regulation of Intracellular Metabolites at 48 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
Taurodeoxycholic acid22.41Glucosamine3.09Taurine2.25
2-Aminooctanoic acid21.352-oxo-4-methylthiobutanoate2.97Citrate2.23
Cystine17.21Allantoate2.96Methylnicotinamide2.20
Orotidine-5-phosphate6.65Ribose-phosphate2.94N-acetyl-glutamate2.19
2-dehydro-D-gluconate5.77Asparagine2.89Ornithine2.19
Shikimate5.29Methylcysteine2.82Hexose-phosphate2.17
Acetylphosphate4.87Glycerate2.672,3-dihydroxybenzoic acid2.14
4-Pyridoxic acid4.76N-Acetyl-L-ornithine2.66Citrate-isocitrate2.11
2-ketohaxanoic acid4.57Pantothenate2.63Acetyllysine2.11
Shikimate-3-phosphate4.43Kynurenic acid2.56Thymidine2.08
Cytidine4.14Folate2.56Allantoin2.06
Threonine4.10Succinate2.54Tryptophan2.06
Aconitate4.066-phospho-D-gluconate2.54Histidine2.04
N6-Acetyl-L-lysine3.92Nicotinamide2.49Cholesteryl sulfate2.04
Cellobiose3.71Hydroxyisocaproic acid2.47Thymine2.03
Citrulline3.45Glucono-D-lactone2.41Orotate2.02
N-carbamoyl-L-aspartate3.27N-acetyl-glutamine2.40Deoxyribose-phosphate2.01
Serine3.19CDP2.31L-arginino-succinate2.01
Xanthine3.14Lysine2.30
Indole3.10Phenylpropiolic acid2.26
Down-Regulation of Intracellular Metabolites at 48 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
Glutathione disulfide61.87ATP4.17IMP2.78
AMP6.01NAD+4.071-Methyladenosine2.35
Glycerophosphocholine5.30GMP3.62Kynurenine2.27
UDP-D-glucose4.51UMP3.25
Up-Regulation of Extracellular Metabolites at 48 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
GMP2.24Serine2.94Adenine5.51
Glycerate2.73UMP4.12Cytidine6.18
S-ribosyl-L-homocysteine2.77Methylnicotinamide4.79
Down-Regulation of Extracellular Metabolites at 48 h of Exposure
NameFold ChangeNameFold ChangeNameFold Change
Glutathione disulfide65.53UDP-D-glucose3.63Allantoin2.47
CDP17.44Phenylpropiolic acid3.58Betaine2.42
6-phospho-D-glucono-1,5-lactone10.16Cholesteryl sulfate3.57AMP2.40
2-ketohaxanoic acid7.30Uric acid3.51Thymine2.37
Phenylpyruvate6.62Methylmalonic acid3.42Ng,NG-dimethyl-L-arginine2.36
Shikimate-3-phosphate6.57Deoxyribose-phosphate3.28p-hydroxybenzoate2.33
Citrate6.33Kynurenic acid3.21Ribose-phosphate2.29
Citraconic acid5.77Hydroxyphenylacetic acid3.17Acadesine2.26
Citrate-isocitrate5.23Acetylphosphate3.14Hydroxyisocaproic acid2.26
N-carbamoyl-L-aspartate5.082-oxo-4-methylthiobutanoate3.01Biotin2.24
Octulose-1,8-bisphosphate (OBP)5.07Nicotinamide2.95SN-glycerol-3-phosphate2.24
2-dehydro-D-gluconate5.01NAD+2.91Sarcosine2.21
D-gluconate4.79Shikimate2.84Imidazoleacetic acid2.18
Xanthine4.11Purine2.83Acetylcarnitine DL2.11
UTP4.11Lactate2.662-Hydroxy-2-methylthiobutanoate2.10
Allantoate4.00Thiamine2.64Malate2.05
L-arginino-succinate3.96Acetyllysine2.58Pantothenate2.00
A-ketoglutarate3.891,3-diphopshateglycerate2.51
4-Pyridoxic acid3.86Thymidine2.50
N-Acetyl-L-alanine3.83Hexose-phosphate2.48

Share and Cite

MDPI and ACS Style

Kotawong, K.; Chaijaroenkul, W.; Roytrakul, S.; Phaonakrop, N.; Na-Bangchang, K. Screening of Molecular Targets of Action of Atractylodin in Cholangiocarcinoma by Applying Proteomic and Metabolomic Approaches. Metabolites 2019, 9, 260. https://doi.org/10.3390/metabo9110260

AMA Style

Kotawong K, Chaijaroenkul W, Roytrakul S, Phaonakrop N, Na-Bangchang K. Screening of Molecular Targets of Action of Atractylodin in Cholangiocarcinoma by Applying Proteomic and Metabolomic Approaches. Metabolites. 2019; 9(11):260. https://doi.org/10.3390/metabo9110260

Chicago/Turabian Style

Kotawong, Kanawut, Wanna Chaijaroenkul, Sittiruk Roytrakul, Narumon Phaonakrop, and Kesara Na-Bangchang. 2019. "Screening of Molecular Targets of Action of Atractylodin in Cholangiocarcinoma by Applying Proteomic and Metabolomic Approaches" Metabolites 9, no. 11: 260. https://doi.org/10.3390/metabo9110260

APA Style

Kotawong, K., Chaijaroenkul, W., Roytrakul, S., Phaonakrop, N., & Na-Bangchang, K. (2019). Screening of Molecular Targets of Action of Atractylodin in Cholangiocarcinoma by Applying Proteomic and Metabolomic Approaches. Metabolites, 9(11), 260. https://doi.org/10.3390/metabo9110260

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop