Thymoquinone, piperine, and sorafenib combinations attenuate liver and breast cancers progression: epigenetic and molecular docking approaches
BMC Complementary Medicine and Therapies volume 23, Article number: 69 (2023)
Traditional herbal medicine has been used for centuries to cure many pathological disorders, including cancer. Thymoquinone (TQ) and piperine (PIP) are major bioactive constituents of the black seed (Nigella sativa) and black pepper (Piper nigrum), respectively. The current study aimed to explore the potential chemo-modulatory effects, mechanisms of action, molecular targets, and binding interactions after TQ and PIP treatments and their combination with sorafenib (SOR) against human triple-negative breast cancer (MDA-MB-231) and liver cancer (HepG2) cells.
We determined drug cytotoxicity by MTT assay, cell cycle, and death mechanism by flow cytometry. Besides, the potential effect of TQ, PIP, and SOR treatment on genome methylation and acetylation by determination of DNA methyltransferase (DNMT3B), histone deacetylase (HDAC3) and miRNA-29c expression levels. Finally, a molecular docking study was performed to propose potential mechanisms of action and binding affinity of TQ, PIP, and SOR with DNMT3B and HDAC3.
Collectively, our data show that combinations of TQ and/or PIP with SOR have significantly enhanced the SOR anti-proliferative and cytotoxic effects depending on the dose and cell line by enhancing G2/M phase arrest, inducing apoptosis, downregulation of DNMT3B and HDAC3 expression and upregulation of the tumor suppressor, miRNA-29c. Finally, the molecular docking study has identified strong interactions between SOR, PIP, and TQ with DNMT3B and HDAC3, inhibiting their normal oncogenic activities and leading to growth arrest and cell death.
This study reported TQ and PIP as enhancers of the antiproliferative and cytotoxic effects of SOR and addressed the mechanisms, and identified molecular targets involved in their action.
Cancer is a global burden and is ranked as one of the major causes of death especially in low-income nations . Globally, liver cancer is the 6th most common and the third leading cause of cancer mortality (8.3%) after lung (18.0%) and colorectal (9.4%) cancer, with the highest incidence rates found in Asia and Africa . On the other hand, worldwide, breast cancer has now exceeded lung cancer in incidence, with approximately 2.3 million new cases in 2020 and 685,000 deaths .
Chemotherapy is one of the major treatment patterns for cancer, used alone or in combination with other treatment modalities. Chemo-drugs kill cancer cells by interfering with the cell cycle regulating genes, therefore, inhibiting cell growth and proliferation, leading to cell death mainly by apoptosis [4,5,6]. Cancer patients receiving chemotherapy suffer from mild side effects, such as alopecia, constipation, and fatigue, to serious ones, such as sterility, cardiotoxicity, and pulmonary fibrosis. Therefore, reducing these side effects would contribute to a better lifestyle for these patients .
Sorafenib (SOR), doxorubicin, and cisplatin are the most used chemotherapeutic agents suitable for patients with liver and breast cancer. Sorafenib is a protein kinase inhibitor that is active against various protein kinases. In 2008, the US FDA approved SOR (NEXAVAR®) for treating patients with liver and advanced kidney cancers . SOR has well-known anticancer potential (Supplementary File 1), as stated in the CTD database (http://ctdbase.org/). Also, it is used as a therapy for about 890 cancer clinical trials (Supplementary File 2), as recognized in the clinical trials database (https://clinicaltrials.gov/).
Natural products derived from plants have always been used to treat various diseases and most of the present anticancer drugs use active ingredients extracted from plant sources [9,10,11,12,13]. These compounds usually exert their effect through several mechanisms such as activation of apoptotic cell death, cell cycle arrest, and inhibition of angiogenesis .
Thymoquinone (TQ) is a major bioactive constituent of the black seed (Nigella sativa) , it has shown different anticancer activities through cell proliferation inhibition, and apoptosis induction in cancer cells, the possible mechanisms of TQ anticancer activity against various proliferative cancer cells were summarized in the review article by El-Far . Interestingly, TQ was reported to have a significant selectivity against various malignant cells [17,18,19,20,21]. The anticancer potential and clinical trials of TQ were stated in Supplementary File 3 and Supplementary File 4, respectively. The taxonomic hierarchy of Nigella sativa L. is shown in Table 1 .
Piperine (PIP) is the major bioactive alkaloid found in black (Piper nigrum), white and long pepper (Piper longum), it has several actions, including anti-inflammatory and anticancer properties . The anticancer potential and clinical trials of TQ were stated in Supplementary File 5 and Supplementary File 6, respectively. Also, the taxonomic hierarchy of Piper nigrum  and Piper longum  is shown in Table 1.
This study examined the potential inhibitory effect of TQ, PIP, and SOR against human triple-negative breast cancer and hepatocellular carcinoma cells. Moreover, molecular targets and mechanisms involved in such activities were also investigated.
Materials and methodsChemicals and reagents
TQ and PIP were obtained from Sigma-Aldrich Chemical Co. (St. Louis, Missouri, USA), and SOR was purchased from Cipla Ltd, India. All cell culture materials were obtained from Gibco (New York, New York, USA).
Human hepatocellular carcinoma HepG2 and breast cancer MDA-MB-231 cells were supplied from American Type Culture Collection (ATCC). Cells were cultured in a complete DMEM medium and incubated at 37 °C in an atmosphere containing 5% CO2.
HepG-2 and MDA-MB-321 Cells were cultured at 15 × 103 per well in a 96-well plate with 100 µl of complete fresh medium for 24 h before treatment with different concentrations of PIP (12.5–200 µM), TQ (25–400 µM) and SOR (6.25–100 µM) for 48 hs. Cell viability was measured by MTT as previously described , and the IC50 was calculated by nonlinear regression analysis of the dose–response curve in each cell line.
For the determination of IC50 values in the combination treatments of TQ and/or PIP with SOR, HepG2, and MDA-MB-231were treated with (IC10-IC50) doses of TQ and PIP, together with SOR (1.0 – 40.0), then incubated for 48 hs before performing the MTT assay as mentioned above.
Cell line treatment
Cells were treated with half of the predetermined calculated IC50 values for all cellular and molecular analyses for each compound. Both cell lines were treated as follows: culture media or 0.1% dimethyl sulfoxide (DMSO), controls; single treatment with either TQ or PIP or SOR; double treatment with TQ + PIP or TQ + SOR or PIP + SOR and finally triple treatment with TQ + PIP + SOR. Treatment was performed 48 hs before the respective analysis, and experiments were repeated at least three times.
Cell cycle analysis
Cell cycle distribution analysis was carried out by cell cycle assay kit, Elabscience Biotechnology Co., Ltd (Houston, Texas, USA). Following trypsinization, cells were centrifuged at 300 × g for 5 min, resuspended using PBS, and 1.2 ml ethanol was added, and the tube was stored at -20 °C for 1 h then, cells were pelleted by centrifugation, the cell pellet was washed with PBS. 100 μl RNase A reagent was added to resuspend the cells and incubated at 37 °C water bath for 30 min then 400 μl propidium iodide (PI) staining solution was added, mixed, and incubated at 4 °C for 30 min. Finally, cells were analyzed using proper machine settings.
Assessment of apoptosis and necrosis by Annexin V-FITC/PI staining
The influence of TQ, PIP, and SOR on apoptosis in HepG2 and MDA-MB-321 cells were quantified by flow cytometry. In brief, cells were collected, washed with PBS, resuspended in 500 μl of annexin V binding buffer, and added 5 μl of annexin V-FITC/PI solution. Cells were resuspended and darkly incubated at 22º C for 20 min before FACS analysis.
RNA extraction, cDNA synthesis, and quantitative real-time PCR
The mRNA levels of DNA methyltransferase (DNMT3B), histone deacetylase (HDAC3) genes, and miRNA-29c were assessed by qRT-PCR. Total RNA was first isolated using the miRNeasy Mini Kit (Qiagen, Germany) and reverse‑transcribed to cDNA using the QuantiTect Reverse Transcription kit (Qiagen, Germany). Second, qRT-PCR was performed using the qPCR Master Mix kit (Enzynomics, Korea). The qRT-PCR cycles consisted of 10 min at 95 °C, 40 cycles of denaturation at 95 °C for 10 s, annealing at 60 °C for 15 s, and extension at 72 °C for 15 s. The primers for DNMT3B, HDAC3, miRNA-29c, β-actin, and U6 genes are listed in Table 2. The relative expression of DNMT3B and HDAC3 were calculated by the comparative 2−ΔΔCt method  using the endogenous β-actin as a housekeeping gene, while miRNA-29c was calculated using the U6 gene as an endogenous control.
To perform molecular docking of TQ, PIP, and SOR against DNMT3B (target site PDB ID: 6KDL), HDAC3 (target site PDB ID: 4A69), and vascular endothelial growth factor receptor-2 (VEGFR-2) (target site PDB ID:3V2A), we first downloaded from RCSB PDB database (https://www.rcsb.org/) and prepared by BIOVIA Discovery Studio (Vélizy-Villacoublay, France) [28,29,30]. The 6KDL retrieved from PDB is the human DNMT3B-DNMT3L complex, where the A and D chains represent DNMT3B. Also, HDAC3 was represented as A and B chains. Therefore, we selected the A chain of 6KDL, 4A69, and 3V2A for protein preparation by removal of water molecules and all ligands in addition to energy minimization and refinement processes.
In addition, the 3D structures of TQ, PIP, and SOR were obtained from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/). The binding free energy, binding affinity (pKi), and the ligand efficiency of TQ, PIP, and SOR against prepared DNMT3B (6KDL-A) and HDAC3 (4A69-A) were determined using InstaDock software . Finally, BIOVIA Discovery Studio Visualizer software did the visualization of target-ligand interaction.
Statistical analysis was performed using GraphPad prism 8.4.2 (https://www.graphpad.com/). Data were represented as mean ± SEM of three independent experiments. One-way ANOVA followed by Tukey’s multiple comparison tests were used to compare group differences. p < 0.05 was deemed to show statistical significance.
TQ and PIP enhanced the cytotoxicity of SOR in liver and breast cancer cells
The antiproliferative effect of TQ, PIP, and SOR was investigated in HepG2 and MDA-MB-231 cells after 48 hs treatments using the MTT assay. Dose–response curves were used to determine the IC50 values for TQ, PIP, and SOR. A dose-dependent growth inhibition was observed in TQ, PIP, and SOR-treated cells compared to control cells. In HepG2 cells, the IC50 values for TQ, PIP, and SOR were 31.57, 65.62, and 10.83 µM, respectively. On the other hand, the IC50 for the same compounds in MDA-MB-231 cells were 29.92, 102.6, and 23.69 µM, respectively (Figs. 1A-F).
The IC50 values for the combination treatments were investigated and presented in Table 3. In HepG2 cells, treatment with different concentrations (IC10—IC50) of TQ and PIP alone and in combination resulted in a significant decrease in the IC50 value for SOR, with the maximum reduction in IC50 (85.78%) detected in the combined treatment with the predetermined IC50 for TQ and PIP. Similarly, the same combined treatment made the SOR IC50 in MDA-MB-231 cells drop by 75.13%.
Microscopic examination of cells after exposure to TQ, PIP, and SOR for 48 hs showed various degrees of cytotoxicity depending on the treatment and cell line (Figs. 2A-B). HepG2 and MDA-MB-231 showed noticeable lower cell numbers, reduced viability, and morphological alterations such as losing normal architecture, cell shrinkage appearing as round, wrinkled in shape, condensation of cytoplasm, and increase in membrane roughness and losing the capacity to attach to the culture-plate surface. The highest cytotoxicity was observed after treatment of HepG2 with TQ + PIP + SOR and MDA-MB-231 with TQ + SOR. In contrast, control cells displayed no observable morphological change or death signs and grew efficiently, reaching 100% confluency while adhering to the cell culture plate.
Determination of cell death mechanism
The percentages of viable and dead cells were determined after 48 hs of treatment by double annexin V-FITC/PI staining (Figs. 3A and 4A). Flow cytometry analysis indicated that the highest cytotoxic effect showing the lowest number of viable cells was observed after treatment of HepG2 with TQ + PIP + SOR (8.10% ± 2.37) and MDA-MB-231 with TQ + SOR (15.33% ± 3.59) (Figs. 3B and 4B).
The cell death mechanism (apoptosis vs. necrosis) was investigated using flow cytometry by analyzing the cells in (early + late) apoptosis vs. necrotic cells. The highest percentage of apoptotic cells were detected after treatment of HepG2 with PIP (31.77% ± 0.29) and MDA-MB-231 cells with TQ + PIP + SOR (55.0% ± 0.46) (Figs. 3D and 4D).
Effect of TQ, PIP, and SOR on cell cycle phases
The effects of TQ, PIP, and SOR and their combination on the HepG2 and MDA-MB-231 cell cycle progression were evaluated by assessment of DNA content using flow cytometry. Representative histogram data files and cell cycle phases plotted as the percentage of total gated events are shown in Figs. 5A and 6A. In HepG2, all TQ, PIP, SOR, and their combinations treatments caused a significant G2/M cell cycle arrest, with SOR showing the highest (44.16% ± 0.1, p < 0.0001) significant arrest compared to control cells (27.01% ± 0.5) (Fig. 5B). Simultaneously, all treatments significantly reduced the cell percentage undergoing DNA synthesis compared to control cells. The highest significant subG1 peak representing cells with low fragmented DNA was detected in TQ-treated cells (19.61% ± 0.7, p < 0.0001) (Fig. 5B). On the other hand, in MDA-MB-231, TQ followed by SOR single treatments caused the highest significant arrest and cell accumulation in the G2/M (33.93% ± 0.10) and (33.45% ± 0.03), respectively (Fig. 6B). Finally, the highest significant subG1 cell population was detected in PIP + SOR treated cells (30.81% ± 0.5, p < 0.0001), (Fig. 6B).
TQ, PIP, and SOR reduced DNMT3B/HDAC3 and increased miRNA-29c expression
The effect of TQ, PIP, and SOR treatment on the expression of DNMT3B/HDAC3 genes and the tumor suppressor miRNA-29c was examined by RT- PCR.
DNMT3B expression analysis represented in Fig. 7A showed that in HepG2 cells, SOR treatment caused the highest reduction in expression (approximately 5 folds) compared to control. On the other hand, PIP and TQ + PIP + SOR treatments suppressed DNMT3B expression almost entirely in MDA-MB-231 cells (Fig. 7D).
The highest decrease in HDAC3 expression was detected after combined TQ + PIP + SOR (3.8-fold) in HepG2 cells and PIP single treatment (11-fold) in MDA-MB-231, (Figs. 7B and 7E).
Finally, PIP alone and combined with TQ caused the most significant increase in miRNA-29c expression in HepG2 (53-fold) and MDA-MB-231 (58 folds), respectively, compared to control cells (Figs. 7C and 7F).
Molecular Docking Study
Molecular docking of TQ, PIP, and SOR against DNMT3B (6KDL-A), HDAC3 (4A69-A), and VEGFR-2 (3V2A-A) are represented in Table 4. SOR showed the highest binding affinity (pKi: 6.45), (pKi: 5.79), and (pKi: 5.43) toward DNMT3B, HDAC3, and VEGFR-2, respectively. Also, PIP exhibited pKi of 5.21, 4.77, and 4.69 toward DNMT3B, HDAC3, and VEGFR-2, respectively, while TQ exhibited the lowest affinities (4.33, 3.96, and 3.74, respectively). Co-crystallized ligands redocked in DNMT3B, HDAC3, and VEGFR-2 are shown in Figs. 8, 9, and 10, respectively.
TQ, PIP, and SOR interacted with the amino acid residues in the binding site of DNMT3B by 6, 3, and 11 bonds, respectively (Table 5). Hydrophobic interactions were recognized between TQ (1 with PHE:581, 1 with VAL606, 1 with PRO650, and 2 with ARG832), PIP (1 with PRO650, 1 with THY776, and 1 with ARG832), and SOR (1 with VAL606, 2 with PRO650, 2 with VAL657, and 1 with ARG832) and DNMT3B binding site. Also, hydrogen bonds were reported between TQ and SER610 and three bonds between SOR and THY776. Furthermore, SOR interacted with three halogen bonds with THY776 and one charge bond with ARG832.
Eight hydrophobic interactions were recognized between TQ and PHE144 (2), PHE199 (2), PHE200 (2), and LEU266 (2) in HDAC3’s binding site (Table 6). Also, SOR made nine hydrophobic interactions with PRO095 (1), PHE144 (1), HIS172 (2), PHE199 (1), PHE200 (1), and LEU266 (3). In addition, SOR is bound with LEU233 by three halogen bonding, while bound with GLY91 (1), ASP93 (2), and VAL96 (1) with hydrogen bonds. PIP exhibited only three hydrogen bonds with GLU52, VAL73, and LEU100.
The molecular interaction between TQ, PIP, and SOR with VEGFR-2 is represented in Table 7, where they are bound by 4, 3, and 10 bonds, respectively with VEGFR-2’s binding site. TQ bound with PHE258 (2) and PHE288 (1) with three hydrophobic interactions, while PIP bound with VAL171 by one hydrophobic bond. SOR made five hydrophobic interactions with CYC150 (1), VAL159 (1), GLN210 (1), ILE212 (1), and TYR214 (1), while bound with two halogen bindings with GLY134 and TYR214.
SOR was the first systemic compound prescribed by US FDA that significantly increased the survival rate for liver cancer patients . However, the one major disadvantage of SOR is its high toxicity, especially at high doses . Therefore, testing new natural agents that might enhance its effects and allow lower doses, would potentially minimize this toxicity . TQ and PIP are two natural compounds commonly used for medicinal purposes [35, 36]. To the best of our knowledge, combining TQ or PIP with SOR has not been investigated. One study on breast carcinoma xenograft reported a combination between TQ and PIP .
In the current study, we showed that combinations of TQ and/or PIP with SOR have significantly enhanced the latter anti-proliferative and cytotoxic effects in both HepG2 and MDA-MB231 cells with variable potency. Overall, HepG2 cells were more responsive to PIP or SOR since both compounds showed lower IC50 values compared with MDA-MB231 cells. This finding agrees with the fact that SOR is prescribed mainly in the treatment of hepatocellular carcinomas . Interestingly, the IC50 values for SOR were significantly diminished in both HepG2 and MDA-MB-231 cells (up to 85% and 75%, respectively) after combined treatment with TQ and/or PIP at different doses. As a consequence, to reduce clinical doses and reduce the common side effects associated with high doses of chemotherapeutic drugs, it is crucial to decrease the 50% inhibitory cytotoxic dose of the chemotherapeutic drug SOR while maintaining its overall cytotoxicity. This could be done through the actions of the natural compounds TQ and PIP.
TQ and PIP alone treatments showed anti-proliferative effects, which agrees with previous studies that reported a decrease in viability of lymphocyte leukaemia cells when treated with TQ . Another study showed that TQ reduced the cell viability of the Huh-7 hepatocellular carcinoma cell line in a dose-dependent manner . On the other hand, Greenshields et al. observed that PIP reduced the proliferation and invasion of MDA-MB-231 cells .
Many recent reports, including the current study, indicated promising chemo-modulatory actions of TQ when combined with other chemotherapeutic agents used against various cancers [42, 43]. In 2020, it was reported that TQ enhanced docetaxel efficiency in MDA-MB-231 and MCF-7 cells by reducing its effective dose . Moreover, TQ synergistically improved the anticancer activity of doxorubicin and cisplatin in hepatocellular carcinoma HepG2 cells . Similarly, PIP was also reported to exert chemo-modulatory effects on chemotherapeutic drugs in different cancers [46, 47].
We have investigated the potential mechanisms causing this anti-proliferative effect, such as modulation of cell cycle interphase, cell death (apoptosis vs. necrosis), and epigenetic genome modifications. The cell cycle is vital in maintaining cell growth and tissue homeostasis, abnormalities in cell cycle progression results in serious diseases, including cancer . We showed that SOR single treatment has resulted in significant G2/M cell cycle arrest in HepG2 and MDA-MB231 cells, which agrees with a previous study that was performed on Hep3B, HepG2, PLC-PRF-5, and SK-Hep1, human hepatocellular carcinoma cell lines . However, other studies reported SOR to arrest cells at G0/G1 phase [50, 51]. Similar to SOR and TQ caused a significant G2/M arrest, TQ was previously reported to cause a G2/M phase arrest in HepG2 cells . However, others have reported that exposure of HepG2 cells to TQ arrested cells at the G0/G1 phase . Depending on the cell line under investigation and tumor behaviour, PIP and TQ were previously suggested to cause cell cycle arrest at G0/G1 and G2/M phases by binding to different cell cycle regulatory proteins [54, 55]. This observation explains the different cell cycle arrest phases that might be observed using the same compound.
Impaired apoptosis is reported as one of the mechanisms leading to SOR resistance in cancer cells and the antiapoptotic protein, B-cell lymphoma 2, is suggested to modulate this phenomenon . Many reports have attributed the anticancer effects of TQ and PIP to its selective proapoptotic actions through the regulation of p53, Bcl-2-associated X protein, and B-cell lymphoma 2 equilibrium and activation of the caspase enzymes [46, 56]. We have examined the apoptotic and necrotic cell death after individual and combined TQ, PIP, and SOR treatment. Notably, necrosis-related cell death was the most observed cause of cell death in HepG2 cells regardless of the treatment. In contrast, in MDA-MB231 cells, apoptotic cell death was prevalent in TQ and SOR single treatments and triple TQ + PIP + SOR treatments.
Interestingly, a combination of TQ or PIP with SOR increased the population of necrotic cells by 6.5 and 8 folds, respectively, compared to SOR alone treatment. Such data might be important in future combination therapy studies aiming to bypass apoptosis as a major cell death mechanism in SOR-resistant cancer cells. Bypassing cancer drug resistance by induction of necroptosis is under clinical investigation as a promising therapeutic mechanism against apoptosis-resistant cancer cells .
The current study showed that both TQ and PIP potentiated the anti-HDAC3 and anti-DNMT3B effect of SOR, causing a significant drop in HDAC3 and DNMT3B expression levels. HDACs and DNMT3B act as transcriptional repressors for tumor suppressor genes, high mRNA levels of HDAC3 and DNMT3B are associated with poor prognosis in patients with different types of cancer [58, 59]. Combination therapy with histone deacetylases (HDACi) or (DNMT3B) inhibitors are currently under investigation for achieving a full cancer therapeutic potential [60, 61]. Several studies have already reported the synergistic effects of treatment with HDACi and various chemotherapeutic drugs [62, 63]. SOR, TQ, and PIP were reported to act as HDACi and DNMT3B inhibitors [64,65,66]. Furthermore, molecular docking revealed that TQ, PIP, and SOR formed hydrophobic and hydrogen bonds with DNMT3B and HDACi, indicating their affinity for inhibiting DNMT3B and HDACi.
In the same way, TQ, PIP, and SOR can bind effectively to VEGFR-2 and inhibit angiogenesis. More research is needed to figure out how TQ, PIP, SOR, and their combinations affect cancer angiogenesis.
The miR-29C has been identified as a crucial miRNA in several cancers; it regulates several oncogenic processes, including epigenetics. The miR-29 family was reported to function as tumor suppressor microRNA through sponging DNMT3A and DNMT3B. However, the biological activity of miR-29C in cancer development and progression is still disputed, most studies reported miR-29C as a tumor suppressor, and others suggested it acts as an oncogene . miR-29b levels were oppositely related to HDAC levels in vitro and human tissue . To the best of our knowledge, the effect of SOR, TQ, and PIP on miR-29C expression levels was not investigated before; the current study confirmed the tumor suppressor role of miR-29C since all our treatments have caused a significant increase in miR-29C expression levels with the highest increase detected in PIP alone (53-fold in HepG2) and PIP + TQ (58-fold in MDA-MB-231).
The molecular docking study was based on assessing the binding energy and binding affinity of ligand-receptor interactions as reflected by the docking score in kcal/mol. A lower docking score indicates a higher binding affinity. The molecular docking results showed a strong affinity of SOR to DNMT3B, and HDAC3, followed by PIP and then TQ.
DNA methylation occurs almost exclusively in CpG islands in mammals catalyzed by DNA methyltransferases (DNMTs). DNMT1 has been shown to maintain methylation in somatic cells, and DNMT3a and DNMT3B are thought to be involved in de novo DNA methylation in embryonic stem cells and early embryos. It was recently found that DNMT1, DNMT3a, and DNMT3B are overexpressed in several human tumors, compared to levels in corresponding normal tissues [69,70,71]. Preclinical and clinical studies demonstrated that the anticancer properties of bioactive components (i.e., parthenolide, folate, retinoids, etc.) may be attributed to its influence on epigenetic processes through binding to DNMT1 enzymatic centre or/and disrupting DNMT transcription . Inhibiting DNMT3B with RNA interference or a selective DNMT3B inhibitor effectively suppressed the DNMT3B activity in vitro and in vivo. Our docking study showed strong binding between SOR, PIP, and TQ to DNMT3B, which inhibited DNMT3B activity. This finding agrees with another study that showed that DNMT3B inhibitor showed a synergistic effect with SOR in the SOR-resistant Hep3B cells .
Genome acetylation is one of the most important epigenetic modifications in promoters of genes that regulate the cell cycle, differentiation, cell growth, and survival in cancer . The acetylation of active genes is controlled by the expression and activity of HDAC . Recently, several reports have shown some dietary phytochemicals as HDAC inhibitors [76,77,78]. HDAC inhibitors lead to increased histone acetylation and transcriptional upregulation of gene expression, inducing cancer cell cycle arrest, apoptosis, and necrosis in a variety of transformed cell lines by several mechanisms depending on the cancer type; HDAC inhibitors, and doses . We showed strong binding between SOR, PIP, and TQ to HDAC, which would probably inhibit HDAC biological activity through disrupting protein–protein interactions important for HDAC activity, as recently suggested . Induction of apoptosis seems to be the predominant route of HDACi-induced cell death ; this agrees with our data showing apoptosis as the main cell death mechanism in PIP-treated HepG2 cells.
The present study showed that TQ and PIP enhance SOR's anti-tumor and cytotoxic effects against hepatocellular carcinoma and triple-negative breast cancer cells. In our study, we studied cell cycle arrest, cell toxicity, apoptosis, and necrosis induction, as well as molecular epigenetic mechanisms that reduce DNMT3B and HDAC3 expression and increase miR-29C expression. Additionally, the results of the molecular docking study validated and supported our in vitro findings.
Availability of data and materials
All data generated or analyzed during this study are included in this published article (and its supplementary information files).
Organization, W.H., World health statistics 2018: monitoring health for the SDGs, sustainable development goals. 2018: World Health Organization.
Petrick JL, et al. International trends in hepatocellular carcinoma incidence, 1978–2012. Int J Cancer. 2020;147(2):317–30.
Sung H, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Brighton, D. and M. Wood, The Royal Marsden Hospital handbook of cancer chemotherapy: a guide for the multidisciplinary team. 2005: Churchill Livingstone.
Kelland, L.R., Cancer cell biology, drug action and resistance. The royal marsden hospital handbook of cancer chemotherapy, 2005: p. 1–15.
Scurr, M., I. Judson, and T. Root, Combination chemotherapy and chemotherapy principles. Cancer chemotherapy’, Eds Brighton D, Wood M, London: Churchill Livingstone, 2005: p. 17–30.
Thackery, E., The Gale Encyclopedia of Cancer: LZ. Vol. 2. 2002: Gale Cengage.
Siegel AB, et al. Sorafenib: where do we go from here? Hepatology. 2010;52(1):360–9.
Cragg GM, Newman DJ. Plants as a source of anti-cancer agents. J Ethnopharmacol. 2005;100(1–2):72–9.
Khan H. Medicinal plants in light of history: recognized therapeutic modality. J Evid Based Complementary Altern Med. 2014;19(3):216–9.
Cortez-Trejo M. C. Potential Anticancer Activity of Pomegranate (Punica granatum L.) Fruits of Different Color: In Vitro and In Silico Evidence. Biomolecules. 2022;12:1649.
Malaluan I. N, et al. Antituberculosis and Antiproliferative Activities of the Extracts and Tetrahydrobisbenzylisoquinoline Alkaloids from Phaeanthus ophthalmicus: In Vitro and In Silico Investigations. Philipp J Sci. 2022;151:371–81.
Ma L, et al. Plant natural products: promising resources for cancer chemoprevention. Molecules. 2021;26:933.
Rayan A, Raiyn J, Falah M. Nature is the best source of anticancer drugs: Indexing natural products for their anticancer bioactivity. PLoS ONE. 2017;12(11):e0187925.
Imran M, et al. Thymoquinone: A novel strategy to combat cancer: A review. Biomed Pharmacother. 2018;106:390–402.
H El-Far, A., Thymoquinone anticancer discovery: possible mechanisms. Current drug discovery technologies, 2015. 12(2): p. 80–89.
El-Far AH, Darwish NH, Mousa SA. Senescent colon and breast cancer cells induced by doxorubicin exhibit enhanced sensitivity to curcumin, caffeine, and thymoquinone. Integr Cancer Ther. 2020;19:1534735419901160.
El-Far AH, et al. Thymoquinone and costunolide induce apoptosis of both proliferative and doxorubicin-induced-senescent colon and breast cancer cells. Integr Cancer Ther. 2021;20:15347354211035450.
El-Far AH, et al. Curcumin and Thymoquinone Combination Attenuates Breast Cancer Cell Lines’ Progression. Integr Cancer Ther. 2022;21:15347354221099536.
El-Far AH, et al. Thymoquinone and its nanoformulation attenuate colorectal and breast cancers and alleviate doxorubicin-induced cardiotoxicity. Nanomedicine. 2021;16(17):1457–69.
Saddiq, A.A., et al., Curcumin, thymoquinone, and 3, 3′-diindolylmethane combinations attenuate lung and liver cancers progression. Frontiers in Pharmacology, 2022: p. 2563.
(USDA Plants Database (2022). Nigella sativa L. black cumin. Available at: https://plants.sc.egov.usda.gov/home/plantProfile?symbol=NISA2. Accessed July 28.
Selvendiran K, Sakthisekaran D. Chemopreventive effect of piperine on modulating lipid peroxidation and membrane bound enzymes in benzo (a) pyrene induced lung carcinogenesis. Biomed Pharmacother. 2004;58(4):264–7.
(USDA Plants Database (2022). Piper nigrum. Available at: https://plants.sc.egov.usda.gov/home/plantProfile?symbol=PINI3. Accessed July 28.
(USDA Plants Database (2022). Piper longum. Available at: https://plants.sc.egov.usda.gov/home/plantProfile?symbol=PILO4. Accessed July 28.
Ismail N, et al. Novel combination of thymoquinone and resveratrol enhances anticancer effect on hepatocellular carcinoma cell line. Future J Pharm Sci. 2018;4(1):41–6.
Livak K.J, Schmittgen T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. methods. 2001;25(4):402–40.
Magpantay HD, et al. Antibacterial and COX-2 inhibitory tetrahydrobisbenzylisoquinoline alkaloids from the Philippine medicinal plant Phaeanthus ophthalmicus. Plants. 2021;10:462.
in silico study. Jayanti, D. A. P. I. S., Abimanyu, I. G. A. M. & Azzamudin, H. Spirulina platensis’s phycocyanobilin as an antiangiogenesis by inhibiting VEGFR2-VEGFA pathway in breast cancer. JSMARTech J Smart Bioprospecting Technol. 2021;2:87–91.
Teixeira IS, et al. Computer Modeling Explains the Structural Reasons for the Difference in Reactivity of Amine Transaminases Regarding Prochiral Methylketones. Int J Mol Sci. 2022;23:777.
Hassan NM, et al. Protein-ligand blind docking using QuickVina-W with inter-process spatio-temporal integration. Sci Rep. 2017;7(1):1–13.
Llovet JM, et al. Sorafenib in advanced hepatocellular carcinoma. N Engl J Med. 2008;359(4):378–90.
Almhanna K, Philip PA. Safety and efficacy of sorafenib in the treatment of hepatocellular carcinoma. Onco Targets Ther. 2009;2:261.
Morisaki T, et al. Combining celecoxib with sorafenib synergistically inhibits hepatocellular carcinoma cells in vitro. Anticancer Res. 2013;33(4):1387–95.
Khader M, Eckl PM. Thymoquinone: an emerging natural drug with a wide range of medical applications. Iran J Basic Med Sci. 2014;17(12):950.
Shityakov S, et al. Phytochemical and pharmacological attributes of piperine: A bioactive ingredient of black pepper. Eur J Med Chem. 2019;176:149–61.
Talib WH. Regressions of breast carcinoma syngraft following treatment with piperine in combination with thymoquinone. Sci Pharm. 2017;85(3):27.
Raoul J-L, et al. Sorafenib: experience and better management of side effects improve overall survival in hepatocellular carcinoma patients: a real-life retrospective analysis. Liver Cancer. 2019;8(6):457–67.
Salim LZA, et al. Thymoquinone induces mitochondria-mediated apoptosis in acute lymphoblastic leukaemia in vitro. Molecules. 2013;18(9):11219–40.
Ashour AE, et al. Thymoquinone suppression of the human hepatocellular carcinoma cell growth involves inhibition of IL-8 expression, elevated levels of TRAIL receptors, oxidative stress and apoptosis. Mol Cell Biochem. 2014;389(1):85–98.
Greenshields AL, et al. Piperine inhibits the growth and motility of triple-negative breast cancer cells. Cancer Lett. 2015;357(1):129–40.
Bashmail HA, et al. Thymoquinone enhances paclitaxel anti-breast cancer activity via inhibiting tumor-associated stem cells despite apparent mathematical antagonism. Molecules. 2020;25(2):426.
Fatfat Z, Fatfat M, Gali-Muhtasib H. Therapeutic potential of thymoquinone in combination therapy against cancer and cancer stem cells. World J Clin Oncol. 2021;12(7):522.
Alkhatib MH, Bawadud RS, Gashlan HM. Incorporation of docetaxel and thymoquinone in borage nanoemulsion potentiates their antineoplastic activity in breast cancer cells. Sci Rep. 2020;10(1):1–12.
Jehan S, et al. Thymoquinone selectively induces hepatocellular carcinoma cell apoptosis in synergism with clinical therapeutics and dependence of p53 status. Front Pharmacol. 2020;11:1453.
Fattah A, et al. The Synergistic Combination of Cisplatin and Piperine Induces Apoptosis in MCF-7 Cell Line. Iran J Public Health. 2021;50(5):1037–47.
Jeong S, et al. Piperine synergistically enhances the effect of temozolomide against temozolomide-resistant human glioma cell lines. Bioengineered. 2020;11(1):791–800.
Feitelson, M.A., et al. Sustained proliferation in cancer: Mechanisms and novel therapeutic targets. in Seminars in cancer biology. 2015. Elsevier.
Fernando J, et al. Sorafenib sensitizes hepatocellular carcinoma cells to physiological apoptotic stimuli. J Cell Physiol. 2012;227(4):1319–25.
Dattachoudhury S, et al. Sorafenib inhibits proliferation, migration and invasion of breast cancer cells. Oncology. 2020;98(7):478–86.
Youssef MM, et al. Novel combination of sorafenib and biochanin-A synergistically enhances the anti-proliferative and pro-apoptotic effects on hepatocellular carcinoma cells. Sci Rep. 2016;6(1):1–12.
Sutton KM, Greenshields AL, Hoskin DW. Thymoquinone, a bioactive component of black caraway seeds, causes G1 phase cell cycle arrest and apoptosis in triple-negative breast cancer cells with mutant p53. Nutr Cancer. 2014;66(3):408–18.
Hassan S, et al. In vitro challenge using thymoquinone on hepatocellular carcinoma (HepG2) cell line. Iran J Pharm Res. 2010;4:283–90.
Gali-Muhtasib HU, et al. Molecular pathway for thymoquinone-induced cell-cycle arrest and apoptosis in neoplastic keratinocytes. Anticancer Drugs. 2004;15(4):389–99.
Rather RA, Bhagat M. Cancer chemoprevention and piperine: molecular mechanisms and therapeutic opportunities. Front Cell Dev Biol. 2018;6:10.
Tutusaus A, et al. Antiapoptotic BCL-2 proteins determine sorafenib/regorafenib resistance and BH3-mimetic efficacy in hepatocellular carcinoma. Oncotarget. 2018;9(24):16701.
Jafri A, et al. Induction of apoptosis by piperine in human cervical adenocarcinoma via ROS mediated mitochondrial pathway and caspase-3 activation. EXCLI J. 2019;18:154.
Hu X, Xuan Y. Bypassing cancer drug resistance by activating multiple death pathways–a proposal from the study of circumventing cancer drug resistance by induction of necroptosis. Cancer Lett. 2008;259(2):127–37.
Niederwieser C, et al. Prognostic and biologic significance of DNMT3B expression in older patients with cytogenetically normal primary acute myeloid leukemia. Leukemia. 2015;29(3):567–75.
Zhou L, et al. Prognosis analysis of histone deacetylases mRNA expression in ovarian cancer patients. J Cancer. 2018;9(23):4547.
Cheray M, et al. DNMT inhibitors in cancer, current treatments and future promising approach: inhibition of specific DNMT-including complexes. Epigenetic Diagn Ther. 2015;1(1):37–48.
Suraweera A, O’Byrne KJ, Richard DJ. Combination therapy with histone deacetylase inhibitors (HDACi) for the treatment of cancer: achieving the full therapeutic potential of HDACi. Front Oncol. 2018;8:92.
Liang B-Y, et al. Synergistic suppressive effect of PARP-1 inhibitor PJ34 and HDAC inhibitor SAHA on proliferation of liver cancer cells. J Huazhong Univ Sci Technolog Med Sci. 2015;35(4):535–40.
Venturelli S, et al. Epigenetic combination therapy as a tumor‐selective treatment approach for hepatocellular carcinoma. Cancer. 2007;109(10):2132–41.
Abdullah O, et al. Thymoquinone Is a Multitarget Single Epidrug That Inhibits the UHRF1 Protein Complex. Genes. 2021;12(5):622.
Bayat S, et al. HDACis (class I), cancer stem cell, and phytochemicals: Cancer therapy and prevention implications. Biomed Pharmacother. 2018;97:1445–53.
Garmpis N, et al. Histone Deacetylase Inhibitors in the Treatment of Hepatocellular Carcinoma: Current Evidence and Future Opportunities. J Pers Med 2021;11(3):223.
Kwon JJ, et al. A systematic review of miR-29 in cancer. Mol Ther Oncolytics. 2019;12:173–94.
Amodio N, et al. Therapeutic targeting of miR-29b/HDAC4 epigenetic loop in multiple myeloma. Mol Cancer Ther. 2016;15(6):1364–75.
Nagai M, et al. Expression of DNA (5-cytosin)-methyltransferases (DNMTs) in hepatocellular carcinomas. Hepatol Res. 2003;26(3):186–91.
Oh B-K, et al. DNA methyltransferase expression and DNA methylation in human hepatocellular carcinoma and their clinicopathological correlation. Int J Mol Med. 2007;20(1):65–73.
Robertson KD, et al. The human DNA methyltransferases (DNMTs) 1, 3a and 3b: coordinate mRNA expression in normal tissues and overexpression in tumors. Nucleic Acids Res. 1999;27(11):2291–8.
Liu Z, et al. Modulation of DNA methylation by a sesquiterpene lactone parthenolide. J Pharmacol Exp Ther. 2009;329(2):505–14.
Lai S-C, et al. DNMT3b/OCT4 expression confers sorafenib resistance and poor prognosis of hepatocellular carcinoma through IL-6/STAT3 regulation. J Exp Clin Cancer Res. 2019;38(1):1–18.
Audia JE, Campbell RM. Histone modifications and cancer. Cold Spring Harb Perspect Biol. 2016;8(4):a019521.
Parbin S, et al. Histone deacetylases: a saga of perturbed acetylation homeostasis in cancer. J Histochem Cytochem. 2014;62(1):11–33.
Evans LW, Ferguson BS. Food bioactive HDAC inhibitors in the epigenetic regulation of heart failure. Nutrients. 2018;10(8):1120.
Rajendran P, et al. Dietary phytochemicals, HDAC inhibition, and DNA damage/repair defects in cancer cells. Clin Epigenetics. 2011;3(1):1–23.
Eckschlager T, et al. Histone deacetylase inhibitors as anticancer drugs. Int J Mol Sci. 2017;18(7):1414.
Maolanon AR, Madsen AS, Olsen CA. Innovative strategies for selective inhibition of histone deacetylases. Cell Chem Biol. 2016;23(7):759–68.
Frew AJ, Johnstone RW, Bolden JE. Enhancing the apoptotic and therapeutic effects of HDAC inhibitors. Cancer Lett. 2009;280(2):125–33.
Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB). NA.
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sorafenib anticancer potentials.
Cancer clinical trials of sorafenib.
The anticancer potential of Thymoquinone.
Thymoquinone anticancer clinical trials.
Piperine anticancer potentials.
Piperine anticancer clinical trials.
About this article
Cite this article
El-Shehawy, A.A., Elmetwalli, A., El-Far, A.H. et al. Thymoquinone, piperine, and sorafenib combinations attenuate liver and breast cancers progression: epigenetic and molecular docking approaches. BMC Complement Med Ther 23, 69 (2023). https://doi.org/10.1186/s12906-023-03872-6