Next Article in Journal
The Intriguing Sigma-1 and Sigma-2 Receptors and Their Potential Therapeutic Roles 2.0
Previous Article in Journal
Transfection of Vein Grafts with Early Growth Response Factor-1 Oligodeoxynucleotide Decoy: Effects on Stem-Cell Genes and Toll-like Receptor-Mediated Inflammation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 21
10.3390/ijms242115867
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Dual Functions of Andrographolide in the Epstein–Barr Virus-Positive Head-and-Neck Cancer Cells: The Inhibition of Lytic Reactivation of the Epstein–Barr Virus and the Induction of Cell Death

by
Chukkris Heawchaiyaphum
1,2,
Praphatson Malat
1,3,
Chamsai Pientong
1,4,*,
Sittiruk Roytrakul
5,
Yodying Yingchutrakul
5,
Sirinart Aromseree
3,4,
Supawadee Suebsasana
6,
Panupong Mahalapbutr
7 and
Tipaya Ekalaksananan
1,4,*
1
Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Biotechnology, Faculty of Science and Technology, Rangsit Center, Thammasart University, Pathum Thani 12120, Thailand
3
Faculty of Agriculture and Technology, Nakhon Phanom University, Nakhon Phanom 48000, Thailand
4
HPV & EBV and Carcinogenesis Research Group, Khon Kaen University, Khon Kaen 40002, Thailand
5
Functional Ingredients and Food Innovation Research Group, National Center for Genetic Engineering and Biotechnology, National Science and Technology Development Agency, Pathum Thani 12120, Thailand
6
Department of Pharmaceutical Sciences, Faculty of Pharmacy, Rangsit Center, Thammasat University, Pathum Thani 12120, Thailand
7
Department of Biochemistry, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(21), 15867; https://doi.org/10.3390/ijms242115867
Submission received: 13 September 2023 / Revised: 21 October 2023 / Accepted: 30 October 2023 / Published: 1 November 2023
(This article belongs to the Section Molecular Microbiology)
Browse Figures
Versions Notes

Abstract

:
Andrographolide, a medicinal compound, exhibits several pharmacological activities, including antiviral and anticancer properties. Previously, we reported that andrographolide inhibits Epstein–Barr virus (EBV) lytic reactivation, which is associated with viral transmission and oncogenesis in epithelial cancers, including head-and-neck cancer (HNC) cells. However, the underlying mechanism through which andrographolide inhibits EBV lytic reactivation and affects HNC cells is poorly understood. Therefore, we investigated these mechanisms using EBV-positive HNC cells and the molecular modeling and docking simulation of protein. Based on the results, the expression of EBV lytic genes and viral production were significantly inhibited in andrographolide-treated EBV-positive HNC cells. Concurrently, there was a reduction in transcription factors (TFs), myocyte enhancer factor-2D (MEF2D), specificity protein (SP) 1, and SP3, which was significantly associated with a combination of andrographolide and sodium butyrate (NaB) treatment. Surprisingly, andrographolide treatment also significantly induced the expression of DNA Methyltransferase (DNMT) 1, DNMT3B, and histone deacetylase (HDAC) 5 in EBV-positive cells. Molecular modeling and docking simulation suggested that HDAC5 could directly interact with MEF2D, SP1, and SP3. In our in vitro study, andrographolide exhibited a stronger cytotoxic effect on EBV-positive cells than EBV-negative cells by inducing cell death. Interestingly, the proteome analysis revealed that the expression of RIPK1, RIPK3, and MLKL, the key molecules for necroptosis, was significantly greater in andrographolide-treated cells. Taken together, it seems that andrographolide exhibits concurrent activities in HNC cells; it inhibits EBV lytic reactivation by interrupting the expression of TFs and induces cell death, probably via necroptosis.

1. Introduction

The Epstein–Barr virus (EBV), a member of the gamma-herpesvirus family, is an enveloped virus containing a linear, double-stranded DNA genome. EBV is ubiquitous, infecting more than 90% of the world’s adult population, with transmission mostly mediated via the oral route. The primary infection by EBV occurs in the oral cavity. EBV infection is closely associated with various types of human malignancies, including Burkitt’s lymphoma, Hodgkin’s lymphoma, NK/T-cell lymphoma, nasopharyngeal carcinoma (NPC), EBV-associated gastric carcinoma (EBVaGC), breast carcinoma, and oral squamous cell carcinoma (OSCC) [1,2,3].
EBV exhibits two alternative phases in the course of infection: latent and lytic. During the lytic phase, new virions are produced and released from the host cell to infect new target cells. EBV lytic reactivation is mediated by both host and virus factors. The molecular mechanism underlying the regulation of EBV lytic reactivation has yet to be fully elucidated. Accumulating evidence has demonstrated that the host master transcriptional regulators of cell differentiation, such as BLIMP1 and KLF4, and transcription factors (TFs), such as ATF, Sp1/3, MEF2D, XBPs, CREB family members, C/EBP family members, AP1, and HIF1-α, can activate the Zta promoter (Zp) and Rta promoter (Rp) of EBV by binding to the responsive elements located on these promoters. On the other hand, host TFs, such as YY1, ZEB1, and ZEB2, can suppress the activation of Zp and Rp [4,5,6,7,8,9,10].
Recent evidence suggests that both latent EBV infection and lytic EBV infection contribute to the genesis of EBV-associated malignancies. BZLF1 is a lytic gene of EBV that plays a critical role in transcriptional transactivation to regulate the switching of EBV from latency to lytic replication [7]. The Zta, a product of the BZLF1 gene, contributes to the oncogenesis of EBV-associated malignancies by inducing genome instability, inflammation, cell proliferation, and inhibiting cell death through the activation of several signaling pathways [11,12,13]. In addition, Rta, a product of the BRLF1 gene, can also promote the oncogenesis of EBV-associated malignancies by inducing cell migration via the activation of the IL6/JAK/STAT3 signaling pathway [14]. Moreover, BRLF1 can also activate the expression of LMP1 during EBV lytic replication in epithelial cells. The expression of LMP1 enhances the promotor activity of Zp and Rp [15,16,17,18,19,20]. It is well documented that LMP1 plays critical roles in EBV-mediated malignant transformation, including resistance to apoptosis and increased angiogenesis, invasion, and metastasis in EBV-harboring cancer cells through the activation of oncogenic signaling pathways, such as the NF-κB, PI3-K/Akt and JAK/STAT pathways [1,21,22,23,24].
As mentioned previously, the lytic infection of EBV is not only associated with the transmission of virus, but it also plays a role in the development of various cancers. Therefore, the inhibition of lytic reactivation by using natural compounds is valuable for the prevention and treatment of EBV-associated diseases. Andrographolide is a medicinal compound that is extracted from the herb Andrographis paniculata (Burm. f.) Nees. This compound exhibits antiviral activity, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and EBV [25,26]. Andrographolide inhibits EBV replication by inhibiting the expression of EBV lytic proteins, including Rta, Zta, and EA-D, during lytic activation in P3HR1 cells and thus inhibits the promoter activities of Zta and Rta [27]. Similarly, our previous study also demonstrated that andrographolide inhibits lytic gene expression and hinders the production of viral particles [28,29]. However, the underlying molecular mechanism through which andrographolide performs these functions in each type of epithelial malignancy is still unknown. Therefore, we aimed to fill this knowledge gap using EBV-positive head-and-neck cancer (HNC) cell lines. Here, we demonstrate for the first time that andrographolide exhibits two concurrent activities in such cells. Firstly, it limits EBV lytic reactivation through the inhibition of the expression of relevant EBV genes and EBV virion production via the inhibition of TFs, MEF2D, SP1, and SP3, possibly through epigenetic mechanisms. Secondly, andrographolide also induces cell death, probably via necroptosis in HNC cells.

2. Results

2.1. Andrographolide Modulates the Pattern of Protein Expression in the EBV-Infected NPC Cell Line

To identify the proteins that were markedly differentially expressed in the EBV-infected NPC cell line, HONE1-EBV, we analyzed the proteomic profiles of the treated cells in the four different treatments. A Venn diagram was constructed, which is shown in Figure 1A. In total, 1015 proteins were expressed only in the cells treated with andrographolide in combination with NaB (Figure 1A). In addition, gene ontology (GO) enrichment analysis was conducted on these 1015 proteins using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) as a tool for annotating gene functions. The 1015 proteins analyzed in this study were significantly enriched in biological processes, such as cellular processes (85.5%) and macromolecule metabolic processes (56.4%), among others (Figure 1B). In addition, proteins for binding (76.2%), organic cyclic-compound binding (35.7%), etc., were significantly enriched in the molecular function category (Figure 1C).

2.2. Andrographolide Inhibits EBV Lytic Reactivation in HNC Cell Lines

Our previous study demonstrated that andrographolide inhibits EBV lytic reactivation through the suppression of lytic protein expression and the inhibition of EBV virion production [28,29]. In this study, we further examined the inhibitory effect of andrographolide on the inhibition of EBV lytic reactivation in HNC cell lines via the treatment of EBV-positive HNC cell lines (HONE1-EBV, SCC25-EBV, and HSC1-EBV) with NaB alone, andrographolide alone, or andrographolide combined with NaB. We observed a significant decrease in the expression of EBV lytic genes in HNC cell lines treated with andrographolide and the combination of andrographolide and NaB when compared with NaB treatment (Figure 2A–D).
Additionally, the copy number of the EBV genome was estimated via qPCR. Consistently, the EBV genome’s copy number was significantly lower in EBV-positive HNC cell lines treated with andrographolide alone and the combination of andrographolide and NaB than in those treated with NaB only (Figure 2E). Therefore, these results indicate that andrographolide inhibits EBV lytic reactivation in HNC cells through the suppression of EBV lytic protein and virion production.

2.3. Andrographolide Inhibits EBV Lytic Reactivation via the Dysregulation of TFs

Accumulating evidence has demonstrated that the activation of EBV lytic replication can be initiated or inhibited by TFs, both transcriptional transactivators and transcriptional repressors. Hence, we hypothesized that andrographolide modulates the expression of these TFs to regulate EBV lytic reactivation in HNC cells. To test our hypothesis, liquid chromatography with tandem mass spectrometry (LC-MS/MS) was used to examine the expression levels of TFs. Interestingly, MEF2D, SP1, and SP3 proteins were significantly decreased in EBV-positive HNC cells treated with the combination of andrographolide and NaB when compared with NaB treatment alone (Figure 3A). Consistent with this finding, the expression of MEF2D, SP1, and SP3 genes was significantly decreased in EBV-positive HNC cell lines treated with the combination of andrographolide and NaB when compared with NaB treatment alone (Figure 3B–D). In addition, we also examined the expression of transcriptional repressors via LC-MS/MS. As shown in Figure 3E, the expression of transcriptional repressors ZEB1, ZEB2, PIAS1, and SAMHD1 was significantly increased in EBV-positive HNC cell lines treated with the combination of andrographolide and NaB when compared with NaB treatment alone. Surprisingly, in cells subjected to the combination treatment, the expression of key molecules involved in epigenetic mechanisms, including DNA methyltransferase (DNMT) 1, DNMT3B, histone deacetylase (HDAC) 5, and HDAC9, was significantly greater (Figure 3E–H). Therefore, these results highlight the mechanism used by andrographolide to inhibit the lytic reactivation of EBV by modulating the expression of TFs and epigenetic changes.

2.4. Prediction of the Interaction of Histone Deacetylases (HDACs and TFs) Using Molecular Modeling and Docking

It is well documented that HDACs regulate gene expression via deacetylation and that some HDACs can directly interact with target proteins [30,31]. Therefore, we further examined whether HDAC5 and HDAC9 can directly interact with MEF2D, SP1, and SP3 by using the PSOPIA web tool. The results showed that, based on protein sequences, both HDAC5 and HDAC9 could directly interact with MEF2D, SP1, and SP3, as evidenced by the Averaged One-Dependence Estimator (AODE) scores (Table 1). Consistent with this result, using the Search Tool for the Retrieval of Interacting Genes (STRING) database, we predicted that proteins related to the regulation of EBV lytic replication are able to interact with each other, except for SAMHD1 (Figure 4A).
To confirm whether HDACs can directly bind with TFs, molecular modeling and docking simulation were performed. As expected, HDACs could bind to TFs via several electrostatic interactions and hydrogen bonds (Figure 4B–G and Tables S1 and S2). The coefficient weights of the lowest energy of the interactions of HDAC5/MEF2D, HDAC5/SP1, HDAC5/SP3, HDAC9/MEF2D, HDAC9/SP1, and HDAC9/SP3 were −765.9, −1155.5, −940.7, −815.3, −1019.8, and −1065.3, respectively. These results suggest that HDAC5 and HDAC9 can directly interact with MEF2D, SP1, and SP3, and thus they can inhibit their transactivation activity.

2.5. Andrographolide Induces Cytotoxicity and Cell Death of HNC Cells

To determine whether andrographolide is cytotoxic for HNC cells and induces cell death, the Cell Counting Kit-8 (CCK-8) was used to determine cytotoxicity in both EBV-positive and EBV-negative HNC cell lines. As shown in Figure 5, andrographolide treatment significantly reduced the cell viability of EBV-positive cell lines but not EBV-negative cell lines, when compared with the NaB treatment. This result suggests that andrographolide is cytotoxic for HNC cells and induces cell death, particularly in EBV-positive HNC cells.
To gain more insight into the functional pathway that was influenced by andrographolide treatment, we used the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis to investigate differentially expressed proteins. The result revealed that, in addition to apoptosis, the necroptosis pathway was enriched in andrographolide-treated cells when compared with NaB-treated cells. Furthermore, we confirmed the proteomic results via the determination of cell death using flow cytometry. As expected, andrographolide dramatically induced cell death in HNC cell lines, especially EBV-positive cell lines, via both necrosis and apoptosis when compared with NaB treatment (Figure 6A,B).
Furthermore, we confirmed the expression of key molecules in the necroptosis pathway, namely RIPK1, RIPK3, and MLKL, via qRT-PCR. As expected, the expression of these genes was significantly greater in andrographolide-treated cells (Figure 6C–E). These results suggest that andrographolide induces cell death in HNC cells, probably via the necroptosis pathway.
Overall, we demonstrated, for the first time, that the mechanism through which andrographolide inhibits EBV lytic reactivation in HNC cells is the dysregulation of TFs through epigenetic mechanisms, namely DNA methylation and histone deacetylation. HDAC5 and HDAC9 may also directly interact with MEF2D, SP1, and SP3 proteins and thus inhibit their transactivator activities. At the same time, andrographolide also induces HNC cell death, probably via the necroptosis pathway (Figure 7).

3. Discussion

As stated earlier, EBV infection is closely associated with various types of tumors, including B-cell tumors and epithelial tumors. EBV has two infection cycles: lytic and latent. The lytic replication of EBV is associated with viral transmission and plays an important role in the oncogenesis of EBV-associated malignancies [3,11,32]. Therefore, the inhibition of EBV lytic replication could be beneficial for the treatment of EBV-associated diseases.
Our study is the first report pointing to the inhibition of EBV lytic reactivation and the induction of cell death in EBV-positive HNC cell lines via andrographolide treatment. We demonstrated this using proteomic analysis, PPI network prediction, molecular docking, and in vitro study. Andrographolide treatment inhibits EBV lytic replication in P3HR1 cells by inhibiting the expression of lytic proteins, Zta, Rta, and EA-D, and suppressing the promoter activity of Zp and Rp [27]. Consistent with previous reports, we also showed that andrographolide can inhibit EBV lytic reactivation in EBV-positive HNC cell lines via the inhibition of EBV gene expression and EBV particle production. Similarly, other plant-derived compounds, such as glycyrrhizic acid, (-)-epigallocatechin gallate (EGCG), moronic acid, resveratrol, and protoapigenone, can also inhibit EBV lytic replication via the inhibition of the expression of EBV lytic genes through the suppression of Zp and Rp [33,34,35,36,37].
As mentioned previously, host factors and TFs, among other factors, play an important role in the regulation of EBV lytic replication. In this study, we also clearly showed that andrographolide treatment inhibited the expression of the transcription factors, MEF2D, SP1, and SP3. Dehydroandrographolide, another major compound extracted from A. paniculata, significantly reduces the expression of TFs, c-Fos, c-Jun, and SP1 in SCC9 cells (an oral cancer cell line) [38]. In addition, luteolin, a natural compound extracted from plants, disrupts the binding of SP1 to the Zta and Rta promoters in NA cells during the induction of EBV lytic reactivation via TPA and NaB treatment, as evidenced by chromatin immunoprecipitation (CHIP) [39]. Thus, our findings are consistent with these previous studies in indicating that andrographolide inhibits EBV lytic reactivation by inhibiting the expression of TFs.
Several host factors are involved in the inhibition of EBV lytic replication. The zinc-finger E-box binding factors, ZEB1 and ZEB2, maintain EBV latency by binding to the ZV element on Zp of EBV in EBV-positive B-lymphocyte cell lines [40,41]. In addition, PIAS1 acts synergistically with SAMHD1 to inhibit EBV lytic replication through protein–protein interactions and SUMOylation [42,43]. Our proteome analysis also showed that PIAS1, SAMHD1, ZEB1, and ZEB2, repressors for EBV lytic replication, were significantly increased in the cells treated with andrographolide and were subsequently induced to become lytic with NaB. This result suggests that andrographolide induces the expression of cellular proteins, playing a role in the inhibition of EBV lytic reactivation.
In addition, the expression of the epigenetic factors DNMT1, DNMT3B, HDAC5, and HDAC9 also increased in the cells treated with the combination of andrographolide and NaB. DNA methyltransferase enzymes, especially DNMT1 and DNMT3B, but not DNMT3A, play a critical role in the maintenance of EBV latency via the restriction of EBNA and LMP expression through DNA methylation in B-lymphocyte cell lines [44,45]. On the other hand, the inhibition of DNMT activities by 5-azacytidine can activate the lytic cycle of EBV by inducing the expression of the lytic gene, BZLF1 [46]. In addition to DNMTs, HDACs can modulate the expression of multiple genes of both the host and EBV by modulating chromosome structure [46,47]. The switching of EBV infection phases from latent to lytic can be induced by stimuli, such as HDAC inhibitors [48,49,50,51,52]. In this study, the expression of HDAC5 and HDAC9 was greater in the cells treated with andrographolide and NaB, suggesting that andrographolide inhibits EBV lytic reactivation via the induction of epigenetic machinery.
Accumulating evidence has demonstrated that MEF2 family members can interact with HDAC family members, especially HDAC4 and HDAC5, to regulate the biological and pathophysiological functions of cells [53,54,55,56,57,58]. Consistent with this finding, by using the STRING database, we also found that HDAC5 and HDAC9 potentially interacted with MEF2D. In addition, HDAC4 and HDAC5 inhibited EBV lytic replication by directly interacting with MEF2 protein that was bound with Zp [59].
Previously, we demonstrated that andrographolide suppresses the lytic reactivation of EBV in gastric cancer cell lines by attenuating the expression of host transcription factors, SP1 and MEF2, via epigenetic modifications, especially through the function of HDAC6 and DNMT3A [27,28]. However, in HNC, we demonstrated the inhibitory effect of andrographolide on the EBV lytic replication by modulating the expression of transcription factors SP1, SP3, and MEF2, at least in part, via the epigenetic alterations, especially through the induction of HDAC5 and DNMT3B, suggesting the cell-type-specific mechanism through which andrographolide inhibits the lytic reactivation of EBV.
Andrographolide treatment can induce cell apoptosis via the activation of various signaling pathways in human cancers, such as skin cancer, colon cancer, cervical cancer, and gastric cancer [60,61,62,63]. In NPC, andrographolide treatment also inhibited cell proliferation and induced cell apoptosis via the activation of various signaling pathways [64,65]. In addition, our previous study also demonstrated a strong cytotoxic effect of andrographolide in gastric cancer cell lines by inducing cell death via apoptosis through the induction of proapoptotic proteins, including BCL2L1, EDOG, HRK, and PUMA [28]. By contrast, in this study, andrographolide exhibited strong cytotoxic effects in EBV-positive HNC cells by inducing cell death. By using bioinformatic tools, we found that the necroptosis-associated proteins were predominant in EBV-positive HNC cells treated with andrographolide and NaB. The key genes of necroptosis, namely RIPK1, RIPK3, and MLKL, were also upregulated. Therefore, this study reveals that a novel mechanism through which andrographolide induces cell death, is, at least in part, necroptosis.

4. Materials and Methods

4.1. Cell Lines and Culture Conditions

The NPC cell lines HONE1 and HONE1-EBV (kindly provided by Professor Hironori Yoshiyama, Shimane University, Shimane, Japan) were cultured in an RPMI-1640 medium (Sigma, St. Louis, MO, USA). EBV-negative OSCC cell lines SCC25 and HSC1 (kindly provided by Dr. Tohru Kiyono, National Cancer Center, Tokyo, Japan) and previously established EBV-positive OSCC cell lines SCC25-EBV and HSC1-EBV c [66] were maintained in Dulbecco’s modified Eagle medium/F12 (DMEM/F12; Sigma, St. Louis, MO, USA). All cell lines were supplemented with 10% fetal bovine serum (Gibco, Breda, The Netherlands) and a penicillin–streptomycin solution (Nacalai Tesque Inc., Kyoto, Japan). Cells were cultured at 37 °C in a 5% CO2 incubator.

4.2. Natural Compounds

The andrographolide compound was prepared as previously described [67] by the Department of Pharmaceutical Sciences, Faculty of Pharmacy, Thammasat University, Thailand. The concentrations of andrographolide used for the treatment of HONE1-EBV, HSC1-EBV, and SCC25-EBV cells were 36.3, 35.0, and 35.0 µM, respectively.

4.3. Proteomic Analysis via Liquid Chromatography with Tandem Mass Spectrometry (LC-MS/MS)

Cells were treated with 0.1% DMSO, or 3 mM NaB, or andrographolide at sub-toxic concentration and incubated for 48 h. In addition, cells were pretreated with andrographolide for 3 h, and the lytic cycle was subsequently induced with NaB and further incubated for 48 h. Protein was extracted from cells using TRIzol™ reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The concentration of total proteins was determined using a Bio-Rad protein-assay dye reagent concentration kit. To prepare protein samples for mass spectrometry, tryptic in-gel digestion was performed as described [68].
The trypsin-digested peptides were introduced into LC-MS/MS using an Ultimate3000 Nano/Capillary LC System (Thermo Scientific, Waltham, MA, USA) coupled with a hybrid quadrupole Q-Tof Impact II™ (Bruker, Billerica, MA, USA) equipped with a nano-captive spray ion source. The proteomic analysis via LC-MS/MS was performed as previously described [68].

4.4. Data Processing and Analysis

To identify and quantify peptides and proteins, MaxQuant version 1.6.5.0 was used. Our data were run against the UniProtKB Human reference database (accessed in October 2020) using the Andromeda search engine. Fragment ion tolerance was set to 20 ppm, and the option of matching between runs (0.4 min match time window) was enabled. Trypsin was specified as the proteolytic enzyme with up to two missed cleavage sites authorized. The N-terminal protein acetylation and oxidation of methionine were set as variable modifications, while carbamidomethylated cysteine was set as a fixed modification. Exploration outcomes were determined for a minimum length of seven amino acids (1% peptide and protein FDR). If a minimum of two peptides was compared between sample groups, the label-free protein quantification (LFQ) was calculated. The match-between-runs feature of MaxQuant was only used in experimental replicates.
The functional annotation (GO) of differentially expressed proteins was performed using the DAVID v6.8 web tool (https://david.ncifcrf.gov/, accessed on 3 March 2021). Differentially expressed proteins were analyzed in terms of two factors: biological process and molecular function. Categories with a p value > 0.05 were considered significant, and the top ten categories are presented in this paper. The KEGG pathway analysis was conducted to identify the biological function of differentially expressed proteins. Differentially expressed proteins that were found among the different treatments were identified using the jvenn web tool (http://jvenn.toulouse.inra.fr, accessed on 21 March 2021).
To construct a PPI network, STRING version 11.5 (https://string-db.org/, accessed on 9 September 2021) was used to explore the possible interactions of differentially expressed proteins at the protein level [69]. A PPI score of >0.4 was considered significant. The PPI networks were visualized using Cytoscape software (version 3.8.2) [70], and p < 0.05 was considered a statistically significant difference. The interactions of HDAC5 and HDAC9 with TFs were predicted by using the PSOPIA web tool [71].
To further confirm protein–protein interactions, molecular modeling and molecular docking simulation were performed. The sequences of candidate proteins, including HDAC5 (Uniprot ID: Q9UQL6), HDAC9 (Uniprot ID: Q9UKV0), MEF2D (Uniprot ID: Q14814), SP1 (Uniprot ID: P08047), and SP3 (Uniprot ID: Q02447), were submitted to SWISS-MODEL, https://swissmodel.expasy.org/ (accessed on 24 November 2021) [72], a protein prediction software, to build the putative 3D structures. The stereochemistry of proteins was evaluated using Ramachandran plot analysis (https://zlab.umassmed.edu/bu/rama/, accessed on 28 November 2021). All ionizable amino acid residues in the predicted models were protonated at pH 7.4 using H++ 1.0 (http://biophysics.cs.vt.edu/, accessed on 28 November 2021). Molecular docking simulation was performed using the web-based service ClusPro [73] to analyze the interactions among the candidate proteins. The PyMOL and Discovery Studio 2021 software version 21.1 were used for the visualization of 3D models and molecular docking. The coefficient weights of protein binding were calculated as follows: E = 0.40Erep + (−0.40)Eatt + 600Eelec + 1.00EDARS [74].

4.5. Quantification of Gene Expression via qRT-PCR

The total RNA was extracted from cells using TRIzol™ reagent (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s instructions, and 1 µg of RNA was used to synthesize cDNA using a RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. Gene expression was quantified with a qRT-PCR assay using SsoAdvancedTM SYBR® Green SuperMix (Bio-Rad, Hercules, CA, USA) in the QuantStudio 6 Flex Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as an internal control. The relative mRNA expression level was quantified using the 2−ΔΔCT method. The primers used in the present study are listed in Table 2.

4.6. Quantification of EBV Genome’s Copy Number via qPCR

The EBV genome’s copy number was estimated with qPCR targeting the EBNA1 gene. Briefly, genomic DNA was extracted from the cell-culture supernatant using the lysis buffer–proteinase K method [75]. Calibration curves were generated in parallel with each analysis using sequential dilutions of DNA from the pGEM-T-EBNA1 plasmid. The copy number of the EBV genome in cells was analyzed with a calibration curve using CT values that were obtained from the samples. The CT values were used to calculate the EBV genome’s copy number from a linear regression equation. The results are expressed as the mean (based on three independent experiments) of EBV copies/50 ng of host genomic DNA.

4.7. Determination of Cell Cytotoxicity Using CCK-8

Cell cytotoxicity was determined using CCK-8 (DOJINDO, Kumamoto, Japan). Cells were seeded into a 96-well plate at 5000 cells/well and preincubated at 37 °C under 5% CO2 for 24 h. The cells were treated with DMSO, NaB, or andrographolide. In addition, for the treatment using the combination of NaB and andrographolide, the cells were pretreated with andrographolide for 3h, followed by treatment with NaB, to induce the lytic cycle. Subsequently, the cells were further incubated for 24 and 48 h. Then, the cells were incubated with 10 µL/well of CCK-8 solution for 1–4 h, and absorbance was measured at 450 nm using a Varioskan LUX Multimode Microplate Reader (Thermo Scientific, Waltham, MA, USA).

4.8. Quantification of Cell Death via Flow Cytometry

The dead cells were quantified using the Dead Cell Apoptosis Kit with Annexin V Alexa Fluor™ 488 and Propidium Iodide (PI) Kit (Life Technologies, Carlsbad, CA, USA). Cells were plated in 6-well plates at a seeding density of 2.5 × 105 cells/well and preincubated at 37 °C under 5% CO2 until 90% confluent. The cells were then treated with DMSO, NaB, or andrographolide, or they were pretreated with andrographolide followed by NaB to induce the lytic cycle and further incubated for 18 h. The cells were stained with Annexin V and PI at room temperature for 15 min and further analyzed via flow cytometry using a BD FACSCanto™ II Cell Analyzer (BD bioscience, Franklin Lakes, NJ, USA). Data were analyzed with FlowJo software version 9 (Flowjo, Treestar Inc., Ashland, OR, USA). Staurosporine treatment was used as a positive control for cell apoptosis. The cells stained by both Annexin V and PI were considered late apoptotic cells, those only positive for Annexin V staining were considered early apoptotic cells, and those only positive for PI staining were considered necrotic cells.

4.9. Statistical Analysis

GraphPad Prism (https://www.graphpad.com/, accessed on 28 November 2021) (GraphPad Software Inc., San Diego, CA, USA) was used for all data analysis. Mann–Whitney tests were used to test whether there was a difference between two independent groups, expressed as mean ± standard deviation (SD). All experiments were repeated three times. A probability (p) value of <0.05 was considered statistically significant.

5. Conclusions

In the present study, we have demonstrated for the first time that andrographolide inhibits EBV lytic reactivation in HNC cell lines via the inhibition of TF expression. The expression of TFs is probably regulated through epigenetic mechanisms, such as DNA methylation and histone modifications. Concurrently, andrographolide also induces cell death via necroptosis through the upregulation of key mediators, including genes encoding RIPK1, RIPK3, and MLKL.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms242115867/s1.

Author Contributions

Conceptualization, C.H., T.E. and C.P.; methodology, C.H., P.M. (Praphatson Malat), Y.Y. and S.R.; software, C.H., P.M. (Panupong Mahalapbutr), and S.R.; validation, C.H.; formal analysis, C.H.; investigation, C.H. and P.M. (Praphatson Malat); resources, S.S.; data curation, C.H., T.E. and C.P.; writing—original draft preparation, C.H.; writing—review and editing, C.H., T.E., S.A. and C.P.; visualization, C.H.; supervision, T.E.; project administration, T.E. and C.P.; funding acquisition, S.A. and C.P. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the research program of the Research and Graduate Studies Department at Khon Kaen University and supported by Thammasat University Research Fund, Contract No. TUFT 63/2566.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The published article includes all the datasets generated or analyzed during this study.

Acknowledgments

We would like to acknowledge David Blair, for editing this manuscript via Publication Clinic KKU, Thailand.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dawson, C.W.; Port, R.J.; Young, L.S. The role of the EBV-encoded latent membrane proteins LMP1 and LMP2 in the pathogenesis of nasopharyngeal carcinoma (NPC). Semin. Cancer Biol. 2012, 22, 144–153. [Google Scholar] [CrossRef]
  2. Ko, Y.H. EBV and human cancer. Exp. Mol. Med. 2015, 47, e130. [Google Scholar] [CrossRef] [PubMed]
  3. Young, L.S.; Yap, L.F.; Murray, P.G. Epstein–Barr virus: More than 50 years old and still providing surprises. Nat. Rev. Cancer 2016, 16, 789–802. [Google Scholar] [CrossRef] [PubMed]
  4. Kenney, S.C.; Mertz, J.E. Regulation of the latent-lytic switch in Epstein–Barr virus. Semin. Cancer Biol. 2014, 26, 60–68. [Google Scholar] [CrossRef] [PubMed]
  5. Reusch, J.A.; Nawandar, D.M.; Wright, K.L.; Kenney, S.C.; Mertz, J.E. Cellular differentiation regulator BLIMP1 induces Epstein-Barr virus lytic reactivation in epithelial and B cells by activating transcription from both the R and Z promoters. J. Virol. 2015, 89, 1731–1743. [Google Scholar] [CrossRef]
  6. Nawandar, D.M.; Wang, A.; Makielski, K.; Lee, D.; Ma, S.; Barlow, E.; Reusch, J.; Jiang, R.; Wille, C.K.; Greenspan, D.; et al. Differentiation-dependent KLF4 expression promotes lytic Epstein-Barr virus infection in epithelial cells. PLoS Pathog. 2015, 11, e1005195. [Google Scholar] [CrossRef] [PubMed]
  7. McKenzie, J.; El-Guindy, A. Epstein-Barr virus lytic cycle reactivation. Curr. Top. Microbiol. Immunol. 2015, 391, 237–261. [Google Scholar] [PubMed]
  8. Murata, T.; Noda, C.; Saito, S.; Kawashima, D.; Sugimoto, A.; Isomura, H.; Kanda, T.; Yokoya, K.K.; Tsurumi, T. Involvement of Jun dimerization protein 2 (JDP2) in the maintenance of Epstein-Barr virus latency. J. Biol. Chem. 2011, 286, 22007–22016. [Google Scholar] [CrossRef]
  9. Yu, X.; Wang, Z.; Mertz, J.E. ZEB1 regulates the latent-lytic switch in infection by Epstein-Barr virus. PLoS Pathog. 2007, 3, e194. [Google Scholar] [CrossRef]
  10. Montalvo, E.A.; Cottam, M.; Hill, S.; Wang, Y. YY1 binds to and regulates cis-acting negative elements in the Epstein-Barr virus BZLF1 promoter. J. Virol. 1995, 69, 4158–4165. [Google Scholar] [CrossRef]
  11. Münz, C. Latency and lytic replication in Epstein–Barr virus-associated oncogenesis. Nat. Rev. Microbiol. 2019, 17, 691–700. [Google Scholar] [CrossRef]
  12. Germini, D.; Sall, F.B.; Shmakova, A.; Wiels, J.; Dokudovskaya, S.; Drouet, E.; Vassetzky, Y. Oncogenic properties of the EBV ZEBRA protein. Cancers 2020, 12, 1479. [Google Scholar] [CrossRef] [PubMed]
  13. Deng, Y.; Münz, C. Roles of Lytic Viral Replication and Co-Infections in the Oncogenesis and Immune Control of the Epstein–Barr Virus. Cancers 2021, 13, 2275. [Google Scholar] [CrossRef] [PubMed]
  14. Tung, K.L.; Wu, Y.T.; Liu, C.; Lin, S.C.; Wu, C.H.; Wu, S.Y.; Chang, Y.; Lan, Y.Y. EBV rta-induced IL-6 promotes migration of bystander tumor cells through IL-6R/JAK/STAT3 pathway in vitro. Anticancer Res. 2020, 40, 3255–3264. [Google Scholar] [CrossRef] [PubMed]
  15. Ahsan, N.; Kanda, T.; Nagashima, K.; Takada, K. Epstein-Barr virus transforming protein LMP1 plays a critical role in virus production. J. Virol. 2005, 79, 4415–4424. [Google Scholar] [CrossRef]
  16. Nawandar, D.M.; Ohashi, M.; Djavadian, R.; Barlow, E.; Makielski, K.; Ali, A.; Lee, D.; Lambert, P.F.; Johannsen, E.; Kenney, S.C. Differentiation-dependent LMP1 expression is required for efficient lytic Epstein-Barr virus reactivation in epithelial cells. J. Virol. 2017, 91, e02438-16. [Google Scholar] [CrossRef]
  17. Chang, Y.; Lee, H.H.; Chang, S.S.; Hsu, T.Y.; Wang, P.W.; Chang, Y.S.; Takada, K.; Tsai, C.H. Induction of Epstein-Barr virus latent membrane protein 1 by a lytic transactivator Rta. J. Virol. 2004, 78, 13028–13036. [Google Scholar] [CrossRef]
  18. Boos, H.; Berger, R.; Kuklik-Roos, C.; Iftner, T.; Mueller-Lantzsch, N. Enhancement of Epstein-Barr virus membrane protein (LMP) expression by serum, TPA, or n-butyrate in latently infected Raji cells. Virology 1987, 159, 161–165. [Google Scholar] [CrossRef]
  19. Lu, C.C.; Jeng, Y.Y.; Tsai, C.H.; Liu, M.Y.; Yeh, S.W.; Hsu, T.Y.; Chen, M.R. Genome-wide transcription program and expression of the Rta responsive gene of Epstein–Barr virus. Virology 2006, 345, 358–372. [Google Scholar] [CrossRef]
  20. Yuan, J.; Cahir-McFarland, E.; Zhao, B.; Kieff, E. Virus and cell RNAs expressed during Epstein-Barr virus replication. J. Virol. 2006, 80, 2548–2565. [Google Scholar] [CrossRef]
  21. Zheng, H.; Li, L.L.; Hu, D.S.; Deng, X.Y.; Cao, Y. Role of Epstein-Barr virus encoded latent membrane protein 1 in the carcinogenesis of nasopharyngeal carcinoma. Cell. Mol. Immunol. 2007, 4, 185–196. [Google Scholar] [PubMed]
  22. Ma, X.; Yang, L.; Xiao, L.; Tang, M.; Liu, L.; Li, Z.; Deng, M.; Sun, L.; Cao, Y. Down-regulation of EBV-LMP1 radio-sensitizes nasal pharyngeal carcinoma cells via NF-κB regulated ATM expression. PLoS ONE 2011, 6, e24647. [Google Scholar] [CrossRef] [PubMed]
  23. Lu, Z.X.; Ma, X.Q.; Yang, L.F.; Wang, Z.L.; Zeng, L.; Li, Z.J.; Li, X.N.; Tang, M.; Yi, W.; Gong, J.P.; et al. DNAzymes targeted to EBV-encoded latent membrane protein-1 induce apoptosis and enhance radiosensitivity in nasopharyngeal carcinoma. Cancer Lett. 2008, 265, 226–238. [Google Scholar] [CrossRef] [PubMed]
  24. Xiao, L.; Hu, Z.; Dong, X.; Tan, Z.; Li, W.; Tang, M.; Chen, L.; Yang, L.; Tao, Y.; Jiang, Y.; et al. Targeting Epstein–Barr virus oncoprotein LMP1-mediated glycolysis sensitizes nasopharyngeal carcinoma to radiation therapy. Oncogene 2014, 33, 4568–4578. [Google Scholar] [CrossRef]
  25. Islam, M.T.; Ali, E.S.; Uddin, S.J.; Islam, M.A.; Shaw, S.; Khan, I.N.; Saravi, S.S.S.; Ahmad, S.; Rehman, S.; Gupta, V.K.; et al. Andrographolide, a diterpene lactone from Andrographis paniculata and its therapeutic promises in cancer. Cancer Lett. 2018, 420, 129–145. [Google Scholar] [CrossRef]
  26. Li, X.; Song, Y. Structure and function of SARS-CoV and SARS-CoV-2 main proteases and their inhibition: A comprehensive review. Eur. J. Med. Chem. 2023, 260, 115772. [Google Scholar] [CrossRef]
  27. Lin, T.P.; Chen, S.Y.; Duh, P.D.; Chang, L.K.; Liu, Y.N. Inhibition of the Epstein-Barr virus lytic cycle by andrographolide. Biol. Pharm. Bull. 2008, 31, 2018–2023. [Google Scholar] [CrossRef]
  28. Malat, P.; Ekalaksananan, T.; Heawchaiyaphum, C.; Suebsasana, S.; Roytrakul, S.; Yingchutrakul, Y.; Pientong, C. Andrographolide Inhibits Lytic Reactivation of Epstein-Barr Virus by Modulating Transcription Factors in Gastric Cancer. Microorganisms 2021, 9, 2561. [Google Scholar] [CrossRef]
  29. Malat, P.; Ekalaksananan, T.; Heawchaiyaphum, C.; Suebsasana, S.; Roytrakul, S.; Yingchutrakul, Y.; Pientong, C. Andrographolide Inhibits Epstein–Barr Virus Lytic Reactivation in EBV-Positive Cancer Cell Lines through the Modulation of Epigenetic-Related Proteins. Molecules 2022, 27, 4666. [Google Scholar] [CrossRef]
  30. Park, S.Y.; Kim, J.S. A short guide to histone deacetylases including recent progress on class II enzymes. Exp. Mol. Med. 2020, 52, 204–212. [Google Scholar] [CrossRef]
  31. Gallinari, P.; Di Marco, S.; Jones, P.; Pallaoro, M.; Steinkühler, C. HDACs, histone deacetylation and gene transcription: From molecular biology to cancer therapeutics. Cell Res. 2007, 17, 195–211. [Google Scholar] [PubMed]
  32. Young, L.S.; Rickinson, A.B. Epstein–Barr virus: 40 years on. Nat. Rev. Cancer 2004, 4, 757–768. [Google Scholar] [CrossRef] [PubMed]
  33. Lin, J.C. Mechanism of action of glycyrrhizic acid in inhibition of Epstein-Barr virus replication in vitro. Antivir. Res. 2003, 59, 41–47. [Google Scholar] [CrossRef] [PubMed]
  34. Chang, L.K.; Wei, T.T.; Chiu, Y.F.; Tung, C.P.; Chuang, J.Y.; Hung, S.K.; Li, C.; Liu, S.T. Inhibition of Epstein-Barr virus lytic cycle by (-)-epigallocatechin gallate. Biochem. Biophys. Res. Commun. 2003, 301, 1062–1068. [Google Scholar] [CrossRef]
  35. Chang, F.R.; Hsieh, Y.C.; Chang, Y.F.; Lee, K.H.; Wu, Y.C.; Chang, L.K. Inhibition of the Epstein-Barr virus lytic cycle by moronic acid. Antivir. Res. 2010, 85, 490–495. [Google Scholar] [CrossRef] [PubMed]
  36. Yiu, C.Y.; Chen, S.Y.; Chang, L.K.; Chiu, Y.F.; Lin, T.P. Inhibitory effects of resveratrol on the Epstein-Barr virus lytic cycle. Molecules 2010, 15, 7115–7124. [Google Scholar] [CrossRef]
  37. Tung, C.P.; Chang, F.R.; Wu, Y.C.; Chuang, D.W.; Hunyadi, A.; Liu, S.T. Inhibition of the Epstein-Barr virus lytic cycle by protoapigenone. J. Gen. Virol. 2011, 92, 1760–1768. [Google Scholar] [CrossRef] [PubMed]
  38. Hsieh, M.J.; Chen, J.C.; Yang, W.E.; Chien, S.Y.; Chen, M.K.; Lo, Y.S.; His, Y.T.; Chuang, Y.C.; Lin, C.C.; Yang, S.F. Dehydroandrographolide inhibits oral cancer cell migration and invasion through NF-κB-, AP-1-, and SP-1-modulated matrix metalloproteinase-2 inhibition. Biochem. Pharmacol. 2017, 130, 10–20. [Google Scholar] [CrossRef]
  39. Wu, C.C.; Fang, C.Y.; Hsu, H.Y.; Chen, Y.J.; Chou, S.P.; Huang, S.Y.; Cheng, Y.J.; Lin, S.F.; Chang, Y.; Tsai, C.H.; et al. Luteolin inhibits Epstein-Barr virus lytic reactivation by repressing the promoter activities of immediate-early genes. Antivir. Res. 2016, 132, 99–110. [Google Scholar] [CrossRef]
  40. Kraus, R.J.; Perrigoue, J.G.; Mertz, J.E. ZEB negatively regulates the lytic-switch BZLF1 gene promoter of Epstein-Barr virus. J. Virol. 2003, 77, 199–207. [Google Scholar] [CrossRef]
  41. Ellis, A.L.; Wang, Z.; Yu, X.; Mertz, J.E. Either ZEB1 or ZEB2/SIP1 can play a central role in regulating the Epstein-Barr virus latent-lytic switch in a cell-type-specific manner. J. Virol. 2010, 84, 6139–6152. [Google Scholar] [CrossRef] [PubMed]
  42. Zhang, K.; Lv, D.W.; Li, R. Conserved herpesvirus protein kinases target SAMHD1 to facilitate virus replication. Cell Rep. 2019, 28, 449–459. [Google Scholar] [CrossRef] [PubMed]
  43. Saiada, F.; Zhang, K.; Li, R. PIAS1 potentiates the anti-EBV activity of SAMHD1 through SUMOylation. Cell Biosci. 2021, 11, 127. [Google Scholar] [CrossRef]
  44. Hughes, D.J.; Marendy, E.M.; Dickerson, C.A.; Yetming, K.D.; Sample, C.E.; Sample, J.T. Contributions of CTCF and DNA methyltransferases DNMT1 and DNMT3B to Epstein-Barr virus restricted latency. J. Virol. 2012, 86, 1034–1045. [Google Scholar] [CrossRef]
  45. Guo, R.; Zhang, Y.; Teng, M.; Jiang, C.; Schineller, M.; Zhao, B.; Doench, J.D.; O’Reilly, R.J.; Cesarman, E.; Giulino-Roth, L.; et al. DNA methylation enzymes and PRC1 restrict B-cell Epstein–Barr virus oncoprotein expression. Nat. Microbiol. 2020, 5, 1051–1063. [Google Scholar] [CrossRef] [PubMed]
  46. Saha, A.; Jha, H.C.; Upadhyay, S.K.; Robertson, E.S. Epigenetic silencing of tumor suppressor genes during in vitro Epstein–Barr virus infection. Proc. Natl. Acad. Sci. USA 2015, 112, E5199–E5207. [Google Scholar] [CrossRef]
  47. Ropero, S.; Esteller, M. The role of histone deacetylases (HDACs) in human cancer. Mol. Oncol. 2007, 1, 19–25. [Google Scholar] [CrossRef]
  48. Gorres, K.L.; Daigle, D.; Mohanram, S.; McInerney, G.E.; Lyons, D.E.; Miller, G. Valpromide inhibits lytic cycle reactivation of Epstein-Barr virus. mBio 2016, 7, e00113-16. [Google Scholar] [CrossRef]
  49. Hui, K.F.; Leung, Y.Y.; Yeung, P.L.; Middeldorp, J.M.; Chiang, A.K. Combination of SAHA and bortezomib up-regulates CDKN2A and CDKN1A and induces apoptosis of Epstein-Barr virus-positive Wp-restricted Burkitt lymphoma and lymphoblastoid cell lines. Br. J. Haematol. 2014, 167, 639–650. [Google Scholar] [CrossRef]
  50. Kim, J.H.; Kim, W.S.; Park, C. Sildenafil prevents HDACi-induced Epstein-Barr virus reactivation through the PKG pathway in NK/T cell lymphoma; potential implications for HDACi-mediated fatal complications. Antivir. Res. 2021, 189, 105063. [Google Scholar] [CrossRef]
  51. Feng, W.H.; Kenney, S.C. Valproic acid enhances the efficacy of chemotherapy in EBV-positive tumors by increasing lytic viral gene expression. Cancer Res. 2006, 66, 8762–8769. [Google Scholar] [CrossRef] [PubMed]
  52. Luka, J.; Kallin, B.; Klein, G. Induction of the Epstein-Barr virus (EBV) cycle in latently infected cells by n-butyrate. Virology 1979, 94, 228–231. [Google Scholar] [CrossRef] [PubMed]
  53. Miska, E.A.; Karlsson, C.; Langley, E.; Nielsen, S.J.; Pines, J.; Kouzarides, T. HDAC4 deacetylase associates with and represses the MEF2 transcription factor. EMBO J. 1999, 18, 5099–5107. [Google Scholar] [CrossRef]
  54. Lu, J.; McKinsey, T.A.; Nicol, R.L.; Olson, E.N. Signal-dependent activation of the MEF2 transcription factor by dissociation from histone deacetylases. Proc. Natl. Acad. Sci. USA 2000, 97, 4070–4075. [Google Scholar] [CrossRef] [PubMed]
  55. Lu, J.; McKinsey, T.A.; Zhang, C.L.; Olson, E.N. Regulation of skeletal myogenesis by association of the MEF2 transcription factor with class II histone deacetylases. Mol. Cell 2000, 6, 233–244. [Google Scholar] [CrossRef]
  56. Choi, S.J.; Park, S.Y.; Han, T.H. 14-3-3τ associates with and activates the MEF2D transcription factor during muscle cell differentiation. Nucleic Acids Res. 2001, 29, 2836–2842. [Google Scholar] [CrossRef] [PubMed]
  57. Belfield, J.L.; Whittaker, C.; Cader, M.Z.; Chawla, S. Differential effects of Ca2+ and cAMP on transcription mediated by MEF2D and cAMP-response element-binding protein in hippocampal neurons. J. Biol. Chem. 2006, 281, 27724–27732. [Google Scholar] [CrossRef]
  58. Aude-Garcia, C.; Collin-Faure, V.; Bausinger, H.; Hanau, D.; Rabilloud, T.; Lemercier, C. Dual roles for MEF2A and MEF2D during human macrophage terminal differentiation and c-Jun expression. Biochem. J. 2010, 430, 237–244. [Google Scholar] [CrossRef]
  59. Gruffat, H.; Manet, E.; Sergeant, A. MEF2-mediated recruitment of class II HDAC at the EBV immediate early gene BZLF1 links latency and chromatin remodeling. EMBO Rep. 2002, 3, 141–146. [Google Scholar] [CrossRef]
  60. Asasutjarit, R.; Sooksai, N.; Fristiohady, A.; Lairungruang, K.; Ng, S.F.; Fuongfuchat, A. Optimization of Production Parameters for Andrographolide-Loaded Nanoemulsion Preparation by Microfluidization and Evaluations of Its Bioactivities in Skin Cancer Cells and UVB Radiation-Exposed Skin. Pharmaceutics 2021, 13, 1290. [Google Scholar] [CrossRef]
  61. Sharda, N.; Ikuse, T.; Hill, E.; Garcia, S.; Czinn, S.J.; Bafford, A.; Blanchard, T.G.; Banerjee, A. Impact of Andrographolide and Melatonin Combinatorial Drug Therapy on Metastatic Colon Cancer Cells and Organoids. Clin. Med. Insights Oncol. 2021, 15, 11795549211012672. [Google Scholar] [CrossRef] [PubMed]
  62. Pasha, A.; Kumbhakar, D.V.; Doneti, R.; Kumar, K.; Dharmapuri, G.; Poleboyina, P.K.; SK, H.; Basavaraju, P.; Pasumarthi, D.; SD, A.; et al. Inhibition of Inducible Nitric Oxide Synthase (iNOS) by Andrographolide and In Vitro Evaluation of Its Antiproliferative and Proapoptotic Effects on Cervical Cancer. Oxidative Med. Cell. Longev. 2021, 16, 6692628. [Google Scholar] [CrossRef] [PubMed]
  63. Lim, S.C.; Jeon, H.J.; Kee, K.H.; Lee, M.J.; Hong, R.; Han, S.I. Andrographolide induces apoptotic and non-apoptotic death and enhances tumor necrosis factor-related apoptosis-inducing ligand-mediated apoptosis in gastric cancer cells. Oncol. Lett. 2017, 13, 3837–3844. [Google Scholar] [CrossRef] [PubMed]
  64. Peng, T.; Hu, M.; Wu, T.T.; Zhang, C.; Chen, Z.; Huang, S.; Zhou, X.H. Andrographolide suppresses proliferation of nasopharyngeal carcinoma cells via attenuating NF-B pathway. BioMed Res. Int. 2015, 2015, 735056. [Google Scholar] [CrossRef]
  65. Wu, B.; Chen, X.; Zhou, Y.; Hu, P.; Wu, D.; Zheng, G.; Cai, Y. Andrographolide inhibits proliferation and induces apoptosis of nasopharyngeal carcinoma cell line C666-1 through LKB1-AMPK-dependent signaling pathways. Die Pharm.-Int. J. Pharm. Sci. 2018, 73, 594–597. [Google Scholar]
  66. Heawchaiyaphum, C.; Iizasa, H.; Ekalaksananan, T.; Burassakarn, A.; Kiyono, T.; Kanehiro, Y.; Yoshiyama, H.; Pientong, C. Epstein–Barr Virus Infection of Oral Squamous Cells. Microorganisms 2020, 8, 419. [Google Scholar] [CrossRef]
  67. Aromdee, C.; Suebsasana, S.; Ekalaksananan, T.; Pientong, C.; Thongchai, S. Stage of action of naturally occurring andrographolides and their semisynthetic analogues against herpes simplex virus type 1 in vitro. Planta Med. 2011, 77, 915–921. [Google Scholar] [CrossRef]
  68. Kaewseekhao, B.; Naranbhai, V.; Roytrakul, S.; Namwat, W.; Paemanee, A.; Lulitanond, V.; Chaiprasert, A.; Faksri, K. Comparative proteomics of activated THP-1 cells infected with Mycobacterium tuberculosis identifies putative clearance biomarkers for tuberculosis treatment. PLoS ONE 2015, 10, e0134168. [Google Scholar] [CrossRef]
  69. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING v10: Protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef]
  70. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  71. Murakami, Y.; Mizuguchi, K. Homology-based prediction of interactions between proteins using Averaged One-Dependence Estimators. BMC Bioinform. 2014, 15, 213. [Google Scholar] [CrossRef] [PubMed]
  72. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, 296–303. [Google Scholar] [CrossRef] [PubMed]
  73. Comeau, S.R.; Gatchell, D.W.; Vajda, S.; Camacho, C.J. ClusPro: An automated docking and discrimination method for the prediction of protein complexes. Bioinformatics 2004, 20, 45–50. [Google Scholar] [CrossRef] [PubMed]
  74. Kozakov, D.; Hall, D.R.; Xia, B.; Porter, K.A.; Padhorny, D.; Yueh, C.; Vajda, S. The ClusPro web server for protein-protein docking. Nat. Protoc. 2017, 12, 255–278. [Google Scholar] [CrossRef] [PubMed]
  75. Sambrook, J.; Fritsch, E.F.; Maniatis, T. Molecular Cloning: A Laboratory Manual, 2nd ed.; Cold Spring Harbor Laboratory Press: Long Island, NY, USA, 1989. [Google Scholar]
Ijms 24 15867 g001
Figure 1. Andrographolide alters the proteomic profile of an EBV-infected NPC cell line (HONE1-EBV): (A) Venn diagram for proteins expressed in HNC cells treated with DMSO (control), NaB alone, andrographolide alone, or a combination of andrographolide and NaB. GO enrichment analysis of differentially expressed proteins was carried out using DAVID. The ten most significantly (p < 0.05) enriched GO terms in biological process (B) and molecular function (C) are presented as bubble diagrams. The Y-axis reveals the top ten functional enrichment results. The X-axis indicates the percentage of genes involved in each biological process or molecular function. The color represents the p value: a range from blue to red indicates a low to high p value, respectively. The size of each bubble indicates the gene numbers involved in each biological process or molecular function. NaB: sodium butyrate, Androg: andrographolide.
Figure 1. Andrographolide alters the proteomic profile of an EBV-infected NPC cell line (HONE1-EBV): (A) Venn diagram for proteins expressed in HNC cells treated with DMSO (control), NaB alone, andrographolide alone, or a combination of andrographolide and NaB. GO enrichment analysis of differentially expressed proteins was carried out using DAVID. The ten most significantly (p < 0.05) enriched GO terms in biological process (B) and molecular function (C) are presented as bubble diagrams. The Y-axis reveals the top ten functional enrichment results. The X-axis indicates the percentage of genes involved in each biological process or molecular function. The color represents the p value: a range from blue to red indicates a low to high p value, respectively. The size of each bubble indicates the gene numbers involved in each biological process or molecular function. NaB: sodium butyrate, Androg: andrographolide.
Ijms 24 15867 g001
Ijms 24 15867 g002
Figure 2. Andrographolide inhibits EBV lytic reactivation in HNC cell lines. Cells were treated with DMSO (control), NaB, andrographolide, or a combination of andrographolide and NaB for 48 h, and the expression of EBV lytic genes, including BZLF1 (A), BRLF1 (B), and BMRF1 (C), was examined using qRT-PCR. The expression of EBV lytic and latent proteins was examined via LC-MS/MS (D). The EBV genome’s copy number was quantified using qPCR (E). NaB: sodium butyrate, Androg: andrographolide. * p < 0.05, *** p < 0.001.
Figure 2. Andrographolide inhibits EBV lytic reactivation in HNC cell lines. Cells were treated with DMSO (control), NaB, andrographolide, or a combination of andrographolide and NaB for 48 h, and the expression of EBV lytic genes, including BZLF1 (A), BRLF1 (B), and BMRF1 (C), was examined using qRT-PCR. The expression of EBV lytic and latent proteins was examined via LC-MS/MS (D). The EBV genome’s copy number was quantified using qPCR (E). NaB: sodium butyrate, Androg: andrographolide. * p < 0.05, *** p < 0.001.
Ijms 24 15867 g002
Ijms 24 15867 g003
Figure 3. Andrographolide inhibits the expression of TFs related to EBV lytic reactivation. Cells were treated for 48 h with NaB, andrographolide, or a combination of andrographolide and NaB. The expression of TFs was examined via LC-MS/MS (A). The expression of MEF2D (B), SP1 (C), and SP3 (D) was determined using qRT-PCR. The expression of repressors was examined via LC-MS/MS (E), and the expression of DNMT3B (F), HDAC5 (G), and HDAC9 (H) was examined using qRT-PCR. NaB: sodium butyrate, Androg: andrographolide. * p < 0.05, *** p < 0.001.
Figure 3. Andrographolide inhibits the expression of TFs related to EBV lytic reactivation. Cells were treated for 48 h with NaB, andrographolide, or a combination of andrographolide and NaB. The expression of TFs was examined via LC-MS/MS (A). The expression of MEF2D (B), SP1 (C), and SP3 (D) was determined using qRT-PCR. The expression of repressors was examined via LC-MS/MS (E), and the expression of DNMT3B (F), HDAC5 (G), and HDAC9 (H) was examined using qRT-PCR. NaB: sodium butyrate, Androg: andrographolide. * p < 0.05, *** p < 0.001.
Ijms 24 15867 g003
Ijms 24 15867 g004
Figure 4. HDAC5 and HDAC9 can directly interact with MEF2D, SP1, and SP3. The PPI was constructed using the STRING database and Cytoscape software version 3.9.1 (A). The interactions of HDAC5/MEF2D (B), HDAC5/SP1 (C), HDAC5/SP3 (D), HDAC9/MEF2D (E), HDAC9/SP1 (F), and HDAC9/SP3 (G) were predicted using ClusPro and visualized with PyMOL and Discovery Studio 2021 software version 21.1.
Figure 4. HDAC5 and HDAC9 can directly interact with MEF2D, SP1, and SP3. The PPI was constructed using the STRING database and Cytoscape software version 3.9.1 (A). The interactions of HDAC5/MEF2D (B), HDAC5/SP1 (C), HDAC5/SP3 (D), HDAC9/MEF2D (E), HDAC9/SP1 (F), and HDAC9/SP3 (G) were predicted using ClusPro and visualized with PyMOL and Discovery Studio 2021 software version 21.1.
Ijms 24 15867 g004
Ijms 24 15867 g005
Figure 5. Andrographolide exhibits a strong cytotoxic effect on EBV-positive HNC cell lines. Cells were treated with DMSO (1% v/v), NaB (3 mM), andrographolide (HONE1 and HONE1-EBV: 36.3 μM; SCC25 and SCC25-EBV: 35.0 μM; and HSC1 and HSC1-EBV: 35.0 μM), or the combination of andrographolide and NaB for 48 h, and cytotoxicity in HONE1 (A), SCC25 (B), and HSC1 (C) was examined using CCK-8. Staurosporine (100 nM) was used as a positive control for cytotoxicity. *** p < 0.001.
Figure 5. Andrographolide exhibits a strong cytotoxic effect on EBV-positive HNC cell lines. Cells were treated with DMSO (1% v/v), NaB (3 mM), andrographolide (HONE1 and HONE1-EBV: 36.3 μM; SCC25 and SCC25-EBV: 35.0 μM; and HSC1 and HSC1-EBV: 35.0 μM), or the combination of andrographolide and NaB for 48 h, and cytotoxicity in HONE1 (A), SCC25 (B), and HSC1 (C) was examined using CCK-8. Staurosporine (100 nM) was used as a positive control for cytotoxicity. *** p < 0.001.
Ijms 24 15867 g005
Ijms 24 15867 g006
Figure 6. Andrographolide induces cell death via the necroptosis pathway. The effect of andrographolide on the induction of cell death in NPC (A) and OSCC (B) cells was assessed via flow cytometry. The expression levels of key mediators of necroptosis, namely RIPK1 (C), RIPK3 (D), and MLKL (E), were examined using qRT-PCR. Staurosporine was used as a positive control for cell apoptosis. *** p < 0.001.
Figure 6. Andrographolide induces cell death via the necroptosis pathway. The effect of andrographolide on the induction of cell death in NPC (A) and OSCC (B) cells was assessed via flow cytometry. The expression levels of key mediators of necroptosis, namely RIPK1 (C), RIPK3 (D), and MLKL (E), were examined using qRT-PCR. Staurosporine was used as a positive control for cell apoptosis. *** p < 0.001.
Ijms 24 15867 g006
Ijms 24 15867 g007
Figure 7. Possible mechanisms through which andrographolide inhibits EBV lytic reactivation and induces cell death in HNC cells. Andrographolide treatment induces the expression of the key mediators of necroptosis, namely RIPK1, RIPK3, and MLKL, resulting in cell death. In addition, andrographolide also induces the silencing of target genes, in particular TF genes (MEF2, SP1, and SP3), via the increased expression of DNMT1, DNMT3B, HDAC5, and HDAC9. HDAC5 and HDAC9 may directly interact with MEF2, SP1, and SP3, resulting in the inhibition of the transcriptional activity of these TFs.
Figure 7. Possible mechanisms through which andrographolide inhibits EBV lytic reactivation and induces cell death in HNC cells. Andrographolide treatment induces the expression of the key mediators of necroptosis, namely RIPK1, RIPK3, and MLKL, resulting in cell death. In addition, andrographolide also induces the silencing of target genes, in particular TF genes (MEF2, SP1, and SP3), via the increased expression of DNMT1, DNMT3B, HDAC5, and HDAC9. HDAC5 and HDAC9 may directly interact with MEF2, SP1, and SP3, resulting in the inhibition of the transcriptional activity of these TFs.
Ijms 24 15867 g007
Table 1. Prediction of protein–protein interactions (PPIs).
Table 1. Prediction of protein–protein interactions (PPIs).
ProteinProteinAODE Scores
Sseq *Sdom *Snet *All *
HDAC5MEF2D0.99620.99650.83510.9991
SP10.8130.51460.83510.9011
SP30.8130.51460.83510.9011
SAMHD10.35370.51460.00000.3197
HDAC9MEF2D0.99620.99650.83510.9991
SP10.8130.51460.83510.9011
SP30.8130.51460.83510.9011
SAMHD10.35370.51460.00000.3197
* Sseq is a score assigned using sequence similarities to a known interacting protein pair. Sdom is a score assigned using the statistical propensities of domain–domain interactions. Snet is a score assigned using a sum of edge weights along the shortest path between homologous proteins in a PPI network. ALL is a score assigned using all three features: Sseq, Sdom, and Snet.
Table 2. Primer sequences.
Table 2. Primer sequences.
GeneForward (5′-3′)Reverse (5′-3′)
BMRF1ACCTGCCGTTGGATCTTAGTGGGCGTTGTTGGAGTCCTGTG
BRLF1TGTTTCAACCGCTCCGACTGGGGTTATGTCGGAGACTGGG
BZLF1TGTTTCAACCGCTCCGACTGGGGTTATGTCGGAGACTGGG
EBNA1CCACAATGTCGTCTTACACCATAACAGACAATGGACTCCCT
HDAC5CCTCAACCATTCCCTCCCACGTTCAGAGGCTGTTTTGCGG
HDAC9CCCCTGCTGCCTCTGTTTTAGGAATTGCCACAAACGCACT
GAPDHTCATCAGCAATGCCTCCTGCATGGGTGGCAGTGATGGCA
MEF2DCATGCCCACTGCCTACAACATGACATTGCCTAGCGACAGC
MLKLCGGCCCTCTGTGGATGAAATGCCTCTCCCAGCTTCTTGTC
RIPK1CGACCTTCTGAGCAGCTTGATCTGAATGCTCTGAGGCAGC
RIPK3CATGGAGAACGGCTCCTTGTGGTTCTGGTCGTGCAGGTAA
SP1CTGTGATACGGATCAGAAACCGTCCACCAAACAATAAAGAGTGCT
SP3CAGAAAGGGTGGGCCTTGAAGCCATCTGTTAAGAGGGCGT
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Heawchaiyaphum, C.; Malat, P.; Pientong, C.; Roytrakul, S.; Yingchutrakul, Y.; Aromseree, S.; Suebsasana, S.; Mahalapbutr, P.; Ekalaksananan, T. The Dual Functions of Andrographolide in the Epstein–Barr Virus-Positive Head-and-Neck Cancer Cells: The Inhibition of Lytic Reactivation of the Epstein–Barr Virus and the Induction of Cell Death. Int. J. Mol. Sci. 2023, 24, 15867. https://doi.org/10.3390/ijms242115867

AMA Style

Heawchaiyaphum C, Malat P, Pientong C, Roytrakul S, Yingchutrakul Y, Aromseree S, Suebsasana S, Mahalapbutr P, Ekalaksananan T. The Dual Functions of Andrographolide in the Epstein–Barr Virus-Positive Head-and-Neck Cancer Cells: The Inhibition of Lytic Reactivation of the Epstein–Barr Virus and the Induction of Cell Death. International Journal of Molecular Sciences. 2023; 24(21):15867. https://doi.org/10.3390/ijms242115867

Chicago/Turabian Style

Heawchaiyaphum, Chukkris, Praphatson Malat, Chamsai Pientong, Sittiruk Roytrakul, Yodying Yingchutrakul, Sirinart Aromseree, Supawadee Suebsasana, Panupong Mahalapbutr, and Tipaya Ekalaksananan. 2023. "The Dual Functions of Andrographolide in the Epstein–Barr Virus-Positive Head-and-Neck Cancer Cells: The Inhibition of Lytic Reactivation of the Epstein–Barr Virus and the Induction of Cell Death" International Journal of Molecular Sciences 24, no. 21: 15867. https://doi.org/10.3390/ijms242115867

APA Style

Heawchaiyaphum, C., Malat, P., Pientong, C., Roytrakul, S., Yingchutrakul, Y., Aromseree, S., Suebsasana, S., Mahalapbutr, P., & Ekalaksananan, T. (2023). The Dual Functions of Andrographolide in the Epstein–Barr Virus-Positive Head-and-Neck Cancer Cells: The Inhibition of Lytic Reactivation of the Epstein–Barr Virus and the Induction of Cell Death. International Journal of Molecular Sciences, 24(21), 15867. https://doi.org/10.3390/ijms242115867

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