Different Cytotoxic Effects of Vetiver Oil on Three Types of Cancer Cells, Mainly Targeting CNR2 on TNBC
Cancer Chemoprevention Research Center, Faculty of Pharmacy, Universitas Gadjah Mada, Sekip Utara II, 55281 Yogyakarta, Indonesia.
Herbal and Traditional Medicine Research Center, Ministry of Health, Republic of Indonesia.
Faculty of Pharmacy, Universitas Gadjah Mada, Sekip Utara II, 55281 Yogyakarta, Indonesia.
*For Correspondence: edy_meiyanto@ugm.ac.idAbstract
Objective:
To investigate vetiver oil (VO) selectivity effects on several cancer cell types and identify the β-caryophyllene role and mechanisms to prevent cancer development.
Methods:
Cytotoxic effects of VO on three types of cancer cells (WiDr, 4T1, T47D) were determined using MTT assay. VO’s effects on the cell cycle and apoptosis were analyzed using flow cytometry. Intracellular Reactive Oxygen Species (ROS) of cells after treatment with VO was observed with DCFDA staining. Bioinformatics study and molecular docking were used to determine the molecular targets of VO.
Results:
VO contained various essential oils in which β-caryophyllene was the most abundant. 4T1 cells performed the lowest IC50 value. WiDr and 4T1 cells showed an arrest in the G2/M phase, while T47D showed an increase of sub G1 population after VO treatment. On the other hand, apoptosis was only observed in WiDr and T47D cells. ROS levels were increased significantly in WiDr and T47D cells but not in 4T1 cells. Cannabinoids CB2 receptor (CNR2) was highly expressed in 4T1 cells and commonly exhibited a low survival rate on Triple Negative Breast Cancer (TNBC) patients. CNR2 was the notable target of β-caryophyllene and performed agonistic interaction, which might have contributed to its cytotoxic activity against 4T1 cells.
Conclusion:
The molecular interaction of VO cannabinoid agonists and the CNR2 receptor was the underlying cause of VO cytotoxicity, which is a VO distinction on TNBC. Therefore, VO is better suited for use as an anti-cancer agent in TNBC cells.
Introduction
Cannabinoids receptor (CNR) has been studied as a potential target for treating various diseases, including the progression of cancer and its deteriorating prognosis (Kisková et al., 2019). There are two types of CNR (CNR1 and CNR2), which belong to the seven trans-membranes superfamily receptors (Raup-Konsavage et al., 2018). CNR1 is highly expressed in brain areas, and it is found to be low in the peripheral nervous system, testes, prostate, uterus, and vascular endothelium. Therefore, this receptor is dominant in cognitive function, memory, anxiety, pain, motor regulation, and endocrine regulation (Dhopeshwarkar and Mackie, 2014). Meanwhile, high expression of CNR2 is found in the immune system cells, but it is also functionally expressed in the brain (Feng et al., 2015; Elbaz et al., 2017). The different tissue distributions of CNR1 and CNR2 allow different receptor effects with selective and specific activation pathways. Both types of CNR are highly expressed in various cancer tissues, but CNR2 plays a more critical role in carcinogenesis and cancer development (Raup-Konsavage et al., 2018).
CNR2 regulates diverse signaling pathways involving Gi/o-dependent inhibition of adenylyl cyclase, stimulation of mitogen-activated protein kinase (MAPK), phosphoinositide 3-kinase (PI3K), and cyclooxygenase-2 (COX-2) signaling, and activation of de novo ceramide synthesis (Raup-Konsavage et al., 2018). Several studies have found that when cancer cells are treated with cannabinoid agonists, the PI3K/Akt and MAPK/ERK signaling pathways are involved in controlling cell proliferation and survival (Khan et al., 2018). The inhibition of Akt leads to cell cycle arrest, which has been directly linked to cannabinoid agonists by controlling cell cycle checkpoints (Park et al., 2011). In leukemic cells, ceramide induces apoptosis by regulating the p38 MAPK signaling, while in lung cancer, the ceramide-dependent pro-apoptotic effect appears to be mediated by the up-regulation of COX-2 expression and the increased synthesis of the pro-apoptotic prostaglandin E-2 (PGE-2), which are triggered by cannabinoid agonists (Hinz and Ramer, 2019).
Cannabinoid compounds, both non-specific and non-selective agonists for CNR1 or CNR2, have been shown to modulate the tumor microenvironment through induction of apoptosis in different models of tumor cells (Hanlon et al., 2016). Interestingly, CNR2-selective cannabinoid ligands are able to treat malignancy of CNR2 expression without causing any psychotropic effects (Feng et al., 2015). β-caryophyllene belongs to the phytocannabinoids, which can activate the cannabinoid receptors (CNR1 and CNR2); however, β-caryophyllene has a solid affinity to CNR2, and it has no affinity to CNR1 (Fidyt et al., 2016). β-caryophyllene is one of the major active essential oils easily found in spices and plants, including in vetivers (Chrysopogon zizanioides L.) (Chahal et al., 2015). Vetiver cultivated from South India showed 0.37% of β-caryophyllene (Chahal et al., 2015). β-caryophyllene has also been widely suggested as a promising agent in cancer prevention and therapy. The anti-proliferative effects of β-caryophyllene on various cancer cell lines were reported by strong suppression of proliferation in two types of colon cancer cells (HCT-116 and HT-29) and of pancreatic cancer cells (PANC-1) (Dahham et al., 2015).
Our study aimed to determine the selectivity of VO effects against different types of cancer cells and identify the role and mechanism of β-caryophyllene to inhibit cancer progression. We used different types of cancer cells, which are 4T1 cells (considered as a triple-negative model), T47D (expresses estrogen receptors), and WiDr cells (with p53 mutations and high expression of COX-2) which also known to have different levels of CNR2 expression. (Martínez-Martínez et al., 2016; Elbaz et al., 2017; Xiang et al., 2018). We used an approach based on physical characteristics and molecular analysis to determine the effect of β-caryophyllene containing VO toward each cell line and the correlation to CNR2. This result will support further explorations of β-caryophyllene-containing essential oil as a chemopreventive agent.
Materials and Methods
Sample preparation and identification of the phytochemical profile
Vetiver (Chrysopogon zizanioides L.) obtained from Materia Medica Batu, Malang, Indonesia, has been identified. Vetiver rhizome weighing 500 grams was steam-distilled with 5 L of aquadest for 4 h in a Clevenger apparatus. The collected distillate was then added with anhydrous sodium sulfate to attract the water-in-oil footprint. Furthermore, VO was dissolved in methanol and then identified its phytochemical profile using Gas Chromatography-Mass Spectrometry (GC-MS) (Shimadzu) with Rtx 5 column, 30 meters long, helium carrier gas, and column temperature of 70-300°C.
Cell culture
The WiDr colon cancer cells (ATCC® CCL ™ 218) were obtained from the Cancer Chemoprevention Research Cancer collection, Faculty of Pharmacy, Universitas Gadjah Mada. TNBC 4T1 cells (ATCC® CRL-253) and luminal-type breast cancer T47D cells (ATCC® HTB-133) were obtained from Professor Masashi Kawaichi, Nara Institute of Science and Technology, Japan. WiDr cells were grown in RPMI medium (Gibco), while 4T1 and T47D cells were grown in DMEM medium (Sigma) supplemented with FBS (Sigma) 10% v/v and 1% penicillin-streptomycin (Sigma) at 37°C in 5% CO2 humidified atmosphere and 95% air in a CO2 incubator.
Cytotoxic assay
WiDr cells (8,000 cells/well), 4T1 cells (2,500 cells/well), and T47D cells (5,000 cells/well) were grown in 96 well plates and subjected to various concentrations of VO. After 24 h of incubation, the cell media was replaced and rinsed with PBS 1x, 5 mg/mL MTT reagent (Sigma) was added and incubated again for 4 h. When formazan was formed, an SDS stopper was added, then incubated overnight. The absorbance was measured using an ELISA microplate reader (BioRad) at 595 nm.
Cell cycle analysis
A density of 2 x 105 cells/well was planted on a 6-well plate and treated with various concentrations of VO, incubated for 24 hours. After cells were harvested, pellets were collected, washed with buffer solution, and stained with 1 mg/mL propidium iodide for cell cycle analysis (BD Cycletest™ Plus DNA Kit; BD Biosciences, San Jose, USA). The measurements were taken by an Accuri C6 (BD Bioscience) flow cytometer, and the data was analyzed with the BD C6 Software.
Apoptosis assay
A density of 2 x 105 cells/well was planted on a 6-well plate and treated with various concentrations of VO, incubated for 24 hours. Cells were then harvested, pellets were collected, washed with buffer solution, and stained with the Annexin-V-FLUOS staining kit (Roche) for 10 minutes for cell cycle analysis (BD Cycletest™ Plus DNA Kit; BD Biosciences, San Jose, USA). The measurements were taken with an Accuri C6 (BD Bioscience) flow cytometer, and the data was analyzed with the BD C6 Software.
Reactive Oxygen Species (ROS) intracellular assay
WiDr, 4T1, and T47D cells were planted on 24 well plates at a 5 × 104 cells/well density. Cells were rinsed using Phosphate Buffer Saline (PBS) and 200 µL of trypsin-EDTA 0.25% (Gibco). Followed by the addition of 500 µL Supplemented Buffer (PBS containing Fetal Bovine Serum (FBS) 10%) 1x, then the cell suspension was collected in a dark microtube and stained with 25 uM Dichlorodihydrofluorescein Diacetate (DCFDA), incubated at 37°C 5% CO2 for 30 minutes. Cells were treated with various concentrations of VO and Doxorubicin as a positive control for 4 hours. Using a BD Accuri C6 flow cytometer (BD Bioscience) with Ex 485 nm/Em 535 nm, measurements were taken by adjusting the number of 20,000 cells.
Gene Expression and Survival Rate
CNR2 and PPARA expression in breast cancer were analyzed through the UALCAN database with a BRCA dataset to obtain potential gene targets for further analysis. Expression was compared between normal versus tumor and normal versus major tumor subclasses in breast cancer. Correlation between survival rate and CNR2 and PPARA expression was analyzed through the OncoLnc database. Then, expression and survival rate data were re-visualized using GraphPad Prism 9 with boxplot and Kaplan-Meier plot, respectively.
Molecular Docking
Molecular docking was conducted using MOE 2010 (Licensed of Faculty of Pharmacy UGM) on CNR2 (PDB ID 5ZTY) as a potential target in breast cancer and GSK3B (PDB ID 3ZRM) as a potential target in colon cancer. The structure of small molecules was created using the Molecular Operating Environment (MOE) molecule builder (version 2010.10 licensed for Faculty of Pharmacy UGM) and then optimized using the MMFF94x force field and conformational search. Small molecules with the lowest docking score were collected, and protein-ligand interactions were analyzed.
Data analysis
All data collected in triplicate were analyzed through One Way ANOVA (IBM SPSS Statistics ver. 25) to identify the significant differences among the treatment.
Results
Distillation and Identification of VO’s Phytochemical Profile
The distilled VO produced a yield of 0.032% w/w, which was yellow in color and had a distinctive smell of vetiver. Phytochemical profile identification carried out by GC-MS showed 33 compounds, most of which were terpenoid compounds. The major active compound of VO was β-caryophyllene (peak 15) with an abundance of 15.43%, which was indicated by m/z = 204. In addition, nine other compounds had high abundance, namely biphenylene (peak 1), β-cedrene (peak 5), 2-carene (peak 8), cuminone (peak 14), ledane (peak 20), β-vetivenene (peak 21), α-humulene (peak 22 and 24), myrcenol (peak 23), propanedinitrile (peak 26) (Figure 1). These results confirmed the presence of β-caryophyllene, which was found to be very high in vetiver.
Interestingly, we realized that among the chemical compounds of VO, we found at least ten compounds possessing biological effects related to cancer. Based on the analysis, β-vetivenene was thought to be the marker of VO’s distinctive smell because it plays an essential role in the aroma characteristic of the oil. We then observed more deeply by literature studies among the ten compounds and their effects on cancer progression and evaluated the protein targets by Swiss-target prediction (Table 1). Most of the activities of VO were thought to come from the terpenoid group, which was known to have an extensive enough abundance that acts as a biomarker, namely β-caryophyllene (15.43%). The main compound of VO, which was β-caryophyllene, was an agonist of phytocannabinoid and a cannabinoids-2 receptor (CB2)-selective. β-caryophyllene had similar activity to Δ9-tetrahydrocannabinol (THC) through activation of CNR2 receptors, leading to cell cycle arrest-induced apoptosis (Fidyt et al., 2016; Chung et al., 2019; Francomano et al., 2019; Arul et al., 2020). In line with this mechanism, another terpenoid compound, namely α-humulene, which was anti-proliferative, pro-apoptotic, and had a strong anti-cancer effect when combined with 5-Fu and oxaliplatin (Ambrož et al., 2019; Chen et al., 2019). The cytotoxicity of VO was also due to the presence of biphenylene, β-cedrene, 2-carene, ledane, and propanedinitrile (Pojarová et al., 2007a; Jeong et al., 2014; Batra et al., 2016; Salim et al., 2016; Elshamy et al., 2017). The results of this study revealed various potential effects associated with the inhibition or killing of cancer cells.
Cytotoxic effect of VO
We then evaluated the cytotoxic effect of VO against cancer cell lines to figure out its selectivity on particular cancer cells. Cell proliferation profiles can be useful in understanding the mechanism of action of specific genes, proteins, and pathways involved in cell survival or death after exposure to toxic agents (Adan, Kiraz and Baran, 2016). We investigated this effect in three cell lines with different tissue origins and molecular markers, namely WiDr colon cancer cells expressing high levels of COX-2; 4T1 breast cancer cells categorized as TNBC; and T47D breast cancer cells that express estrogen receptors. The VO treatment resulted in different IC50 values, which were 302 µg/mL, 60 µg/mL, and 112 µg/mL, respectively (Figure 2). These results indicated that 4T1 was the most sensitive to VO treatment while WiDr was the most resistant. We may note that VO most likely preferred to be against breast cancer cells rather than colon cancer cells, and the TNBC sub-type was more suitable than the luminal subtype. This difference in response was probably caused by the distinctive expressions of protein markers that may affect the physiological process of each type of cell line.
Cell Cycle Inhibition of VO
The cell cycle is the main physiological process in cancer cells in which the cells undergo proliferation over the limit of cell division number. This mechanism is supported by protein kinases called cyclin-dependent kinases (CDKs), which control the activity of kinase and are only active when bound by a specific regulation component, “cycline,” with serine/threonine-specific catalytic core (Shah et al., 2020). Cytotoxic drugs can stop cell proliferation by abrogating cell cycle machinery resulting in cell cycle arrest or cell death. Interestingly VO treatment showed different phenomena among the three cell types. VO treatment caused cell cycle arrest in different phases, known from the experimental results that WiDr cells and 4T1 cells underwent arrest in the G2/M phase (p<0.01) (Figure 3). However, VO treatment on T47D increased sub G1 population (p<0.01), while 4T1 cells did not increase significantly compared to controls. These results indicated that the G2/M arrest was stronger in 4T1 cells. Meanwhile, treatment of VO on T47D cells gave no differences in G1 phase cells accumulation between the control and VO 30 µg/mL. Moreover, cell death induction was found after treatment of 30 and 60 µg/mL of VO.
Apoptosis Effect of VO
Apoptosis results in fragmentation of the cell nucleus genome. It has a vital role in eradicating cancer cells and is an essential target in discovering anti-cancer drugs (Meiyanto et al., 2019). To clarify whether the cell cycle arrest effect under VO treatments is correlated with apoptosis phenomena, we then carried out flow cytometry analysis with PI staining on VO treated cells. Flow cytograms confirmed that the treatment of VO at the concentration 30 and 60 µg/mL increased apoptosis in WiDr and T47D cells (p<0.01), while in 4T1 cells, VO treatment did not induce apoptosis (Figure 4). These findings lead us to give more attention to the molecular events that underlie the phenomena.
ROS scavenging system of VO
In general, cell cycle arrest and apoptosis are attributed to ROS levels in cells. Free radicals cause carcinogenesis and tumor cell progression, which are the characteristics of cancer cells. In cancer cells, low to medium intracellular ROS level plays a role in promoting cell proliferation. A high intracellular ROS level in normal cells causes oxidative DNA damage and leads to cancer formation and cell death. ROS plays an essential role at each stage of cancer cell development, namely in initiation, promotion, and progression. Increased intracellular ROS level leads to an activation of oncogene and inactivation of tumor suppressor genes, angiogenesis, and mitochondrial dysfunction (Kumari et al., 2018; Larasati et al., 2018). Our data revealed that ROS levels were increased significantly in WiDr and T47D cells but not in 4T1 cells (Figure 5). This phenomenon makes sense with the effect of VO treatment on cell cycle progression and apoptosis. VO’s apoptosis induction effect on WiDr and T47D is likely correlated with the increase of intracellular ROS. Hence, intracellular ROS seems to contribute to the molecular events of apoptosis. However, the phenomenon of cell cycle arrest of 4T1 cells might be triggered by other mechanisms.
Gene Expression and Survival Rate Analysis
We then analyzed the protein target that underlies these compounds to have anti-cancer activity. The results of our analysis with the Swiss Target Prediction showed that among the proteins that appeared, two proteins were likely to mediate the various effects seen, namely Cannabinoid CB2 receptor (CNR2) and peroxisome proliferator-activated receptor alpha (PPARA). Both proteins have a probability above 50%, so that it is very potential if used as targets. Interestingly, we found estrogen receptors (ESR1 and ESR2) as markers of T47D cells and GSK3B as markers of expressed WiDr cells. Therefore, we next analyzed the CNR2 and PPARA expression levels on luminal, TNBC, and colon cancer cells to determine precisely what causes the differences in activities between these cells.
TNBC cells expressed CNR2 and PPRA proteins significantly higher than normal cells (p<0.01) (Figure 6A). Luminal cells expressed both proteins in the same level as normal cells (p>0.01). Meanwhile, colon cancer cells showed a very low expression level of both proteins, even lower than controls (Figure 6B). Therefore, it is necessary to look further based on the survival rate of CNR2 and PPARA, which showed that CNR2 and PPARA had a poor survival probability with an endpoint of fewer than 8,000 days.
These findings may be related to the lower IC50 value of 4T1 cells and the sensitivity of the others. Meanwhile, the colon is not the target of CNR2, so it may tend to lead to GSK3B because it has high expression in the colon, but unfortunately, the probability is low, which is still below 50%. Likewise, T47D cells, which have ESR1 and ESR2, also have a high expression with low probability. Another confirmation was evidenced by molecular docking, where the docking scores between the CNR2 and β-caryophyllene ligands were equally low, supported by an RMSD value below 2. This result provides a bright insight into why the effect is not as strong as on 4T1 cells. Nevertheless, we still need to observe further to determine the pathways that elicit the activity by molecular docking.
Molecular Docking
Molecular docking was used to obtain the binding affinities described with docking scores (kcal/mol) between ligand and protein. This study examined the binding affinities of β-caryophyllene and α-humulene with cannabinoid receptor 2 (CNR2) as the potential target in breast cancer. We also subjected bicine and palmitoleic acid as agonists and antagonists to CNR2, respectively. Compared to each ligand, palmitoleic acid (antagonist) had a lower docking score, and bicine (agonist) had a higher docking score (Fig. 7). We could hypothesize that the lower the docking score, the more likely the ligand is to have an antagonist activity. α-humulene was more likely to have an antagonist activity compared to β-caryophyllene. The native ligand 9JU performed the lowest docking score of -12.056 kcal/mol, which had a direct bond with the amino acids Phenylalanine94, Trptophan194, Trptophan258, and Threonine114 (Figure 7). The VO metabolites, namely β-caryophyllene and α-humulene, both were bound to phenylalanine94 with not much different docking scores, namely -8.72 kcal/mol and -8.89 kcal/mol, respectively. On the other hand, Bicine (CNR2 agonist) was bound to Lysine109 and Phenylalanine183 while Palmitoleic acid (CNR2 antagonist) was bound to Leucine182. The amino acid interactions between native ligands, CNR2 agonists, respective CNR2 antagonists, and their respective VO metabolites indicated the potential of VO as a new candidate targeting CNR2. We also conducted molecular docking to glycogen synthase kinase 3 beta (GSK3B) as the potential target in colon cancer. Compared to the native ligand (ZRM), cuminone had a lower docking score than the native ligand. This possibility occured because cuminone did not show a direct interaction with the amino acid GSK3B while the native ligand interacted with Valine70 Lysine85 Valine135
| Peak | Compound | Abundance (%) | KI | Activity | Target protein based on Swiss Target Prediction |
|---|---|---|---|---|---|
| 1 | Biphenylene | 4.61 | 1357 | Cytotoxicity (Batra et al., 2016) | - |
| 5 | β-cedrene | 2.27 | 1424 | - | HSD11B1 |
| 14 | Cuminone | 8.85 | 1509 | - | GSK3B, TYR, IMPDH2, PTPN1, TDO2, HSD11B1 |
| 15 | β-caryophyllene | 15.43 | 1515 | Synergistic anti-cancer effects with doxorubicin, cytotoxicity, induction of apoptosis (Fidyt et al., 2016; Hanušová et al., 2017; Ramachandhiran et al., 2019) | PPARA*, CNR2*, FAAH |
| 19 | 2-carene | 4.3 | 1535 | Cytotoxicity (Liu et al., 2017) | PPARA, CNR2, FAAH, AR, ESR1, CHRM2 |
| 20 | Ledane | 3.6 | 1546 | - | PPARA, CNR2 |
| 21 | β-vetivenene | 3.62 | 1552 | Antioxidant (Kim et al., 2005) | ADORA1, ADORA2A, AR, ESR1, ESR2 |
| 22 | α-humulene | 8.12 | 1556 | Synergistic anti-cancer effect with 5-FU and oxaliplatin; antiproliferation and apoptosis induction (Ambrož et al., 2019; Chen et al., 2019) | PPARA, CNR2 |
| 23 | Myrcenol | 8.52 | 1560 | - | PTGS1, HSD11B1, SIGMAR1, DRD2, AR, HSD17B2 |
| 26 | Propanedinitrile | 3.54 | 1574 | Cytotoxicity (Seckl and Rozengurt, 1996) | - |
Discussion
The purpose of this study was to explore vetiver oil (VO) potential to inhibit the activity of cancer development based on in vitro and in silico studies. Bioinformatics analyses were also carried out to determine the specific target proteins thought to mediate these effects. We used two breast cancer cells, namely luminal type A that expressed ER+ (T47D cells), triple-negative breast cancer or TNBC cell line (4T1 cells), and colon cancer cells that are characterized with a high COX-2 expression (WiDr cells).
The different characteristics of the three cells that we used provide different responses for each cell (Figure 8). The first anti-cancer study evaluated was the degree of cytotoxicity. VO cytotoxicity showed that VO treatment had the highest IC50 value in WiDr cells and the lowest in 4T1 cells. A low IC50 value below 100 µg/mL generally indicates a strong and highly sensitive anti-cancer potential (Prayong, Barusrux, and Weerapreeyakul, 2008). Looking closer at the other breast cancer cells, T47D cells showed an IC50 value above 100 µg/mL, indicating that VO was weakly toxic to T47D cancer cells. This effect was also observed in WiDr cells. This result provides knowledge that VO has a good target for 4T1 breast cancer cells and can be projected as an anti-cancer agent.
Further confirmation on the anti-cancer effect was obtained through the cell death analysis. WiDr and 4T1 cells were in G2/M arrest, but interestingly, apoptosis was only observed in WiDr and T47D cells, whereas in 4T1 cells, there were none of the phenomena compared to controls. G2/M arrest can be caused by the failure of the cell mitotic phase, which in turn will lead to apoptosis and cell death (Pojarová et al., 2007b).
Cell death pathways can also be observed through intervening intracellular ROS metabolic rate. Intracellular ROS level that exceeds the threshold will eventually be toxic for cells (Larasati et al., 2018). Hence, to trigger cancer cell death, agents that can increase ROS levels are needed (Ahlina et al., 2020). Our study showed a significant increase in ROS levels in WiDr and T47D cells; this result corresponds with the previous analysis that showed how both cell lines experienced cell death due to changes in cellular metabolism after VO treatments.
Interestingly, 4T1 cells exhibited different results from the previous two cells. There was no increase in apoptotic cells and intracellular ROS levels compared to controls, making it hard to precisely predict how the mechanism mediates the previously described toxicity. In this study, we only carried out a new look at how the effects occurred during the 24 hours of testing. The effect of TNBC cell death may be obtained if the observation time is extended.
We did not stop there and tried to find other ways that caused this event. Recognizing that the major compounds of VO have various properties with diverse mechanisms of action provided us with a glimmer of insight (Table 1). Based on that fact, the metabolic pathways taken by cancer cells after VO administration were very interesting to be studied further, which would involve specific target proteins. CNR2 was found to have a higher expression than PPARA in all breast cancer cell lines (Figure 6). This prediction matched with a previous study, which revealed that CNR2 was responsible for breast cancer cells (Brown et al., 2020). In addition, CNR2 is also highly expressed in various cancer tissues.
We compared native ligands, bicine agonists, and palmitoleic acid antagonists from CNR2 to obtain compounds that can bind to CNR2 ligands selectively to determine the cause of the cancer-inhibiting effect. Compounds β-caryophyllene and α-humulene performed docking scores close to bicine, making them potential phytocannabinoids (Figure 7). Previous studies have reported that β-caryophyllene functions as phytocannabinoids, but there is no evidence for α-humulene. Our study provides an overview of the synergistic effect of the two phytocannabinoids on CNR2, which mediates the inhibition of cancer malignancy (Figure 8).
In colon cancer, CNR2 is not the main marker for the survival rate due to the low expression level while some are not expressed (Martínez-Martínez et al., 2016) (Figure 6). This finding on colon cancer turns the target to GSK3B protein because it is known for its high expression (Dupasquier et al., 2019; Wu et al., 2019; Evangelisti et al., 2020). GSK3B is one of the critical predictive targets that have a role in developing colorectal cancer through the intracellular signaling cascade and phosphorylation, so we pay attention to this finding (Wulandari et al., 2020). Based on this study, the prediction of GSK3B interaction with cuminone compounds from VO was stronger than native ligands (Figure 7). It can be predicted that cell death in WiDr colon cancer cells was caused by the presence of cuminone in VO (Figure 8).
These findings revealed a great potential of VO as an anti-cancer agent. VO was confirmed to be safe against normal cells by the result of the cytotoxicity test against NIH-3T3 cells, which are fibroblast cell lines representing normal cells (Jenie et al., 2021; Zulfin et al., 2021). Cytotoxicity to NIH-3T3 was considered to be low cytotoxic with an IC50 value of 266 µg/mL (data not shown). These interesting results are expected to open up opportunities for developing VO as an anti-cancer agent of the colon and co-chemotherapeutic agent to increase the effectiveness of chemotherapy agents that often experience resistance problems.
In conclusion, VO gave different cytotoxic effects on 4T1, T47D, and WiDr cells, and 4T1 cells were the most sensitive to VO. VO in WiDr and 4T1 cells caused cell cycle arrest in the G2/M phase, while in T47D cells, it caused cell accumulation in the subG1 phase. The most responsible compounds for these effects are β-caryophyllene and α-humulene, which specifically affect TNBC, targeted on CNR2.
Abbreviations
ROS: reactive oxygen species
TNBC: triple-negative breast cancer
DCFDA: 2′,7′-Dichlorodihydrofluorescein diacetate
CNR2: cannabinoids CB2 receptor
PPARA: peroxisome proliferator-activated receptor alpha
GSK3B: glycogen synthase kinase 3 beta
MAPK: mitogen-activated protein kinase
PI3K: phosphoinositide 3-kinase
COX-2: cyclooxygenase-2
PGE-2: prostaglandin E-2
Conflict of interest
We declare that we do not have a conflict of interest.
Acknowledgments
We thank CCRC, Faculty of Pharmacy, Universitas Gadjah Mada for the main support of this research.