
==== Front
Iran J Basic Med Sci
Iran J Basic Med Sci
IJBMS
Iranian Journal of Basic Medical Sciences
2008-3866
2008-3874
Mashhad University of Medical Sciences Mashhad, Iran

10.22038/ijbms.2024.76691.16599
Original Article
Andrographolide demonstrates anti-proliferative activity in oral cancer by promoting apoptosis, the programmed cell death process
Kumbhar Gauri Mansinh 1
Jadhav Amol Dilip 2
Kheur Supriya 1
Vaibhav Sunil Ladke 3*
1 Dr. D. Y. Patil Dental College & Hospital Dr. D. Y. Patil Vidyapeeth, Sant Tukaram Nagar, Pimpri, Pune. Maharashtra. India Pune: 411018
2 Institute of Applied Biological Research and Development, a Division of Nirav BioSolutions Pvt Ltd, Aundh, Pune, India
3 Dr. D. Y. Patil Madical College, Hospital & Research Centre. Dr. D. Y. Patil Vidyapeeth, (Deemed to be University), Sant Tukaram Nagar, Pimpri, Pune. Maharashtra. India Pune: 411018
* Corresponding author: Ladke Vaibhav Sunil. Dr. D. Y. Patil Madical College, Hospital & Research Centre. Dr. D. Y. Patil Vidyapeeth, (Deemed to be University), Sant Tukaram Nagar, Pimpri, Pune. Maharashtra. India Pune: 411018. Tel/ Fax: +91-02027805217, Email: drvaibhavladke@gmail.com
2024
27 10 13001308
7 12 2023
26 2 2024
https://creativecommons.org/licenses/by/3.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License, (http://creativecommons.org/licenses/by/3.0/) which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Objective(s):

Andrographolide has been studied on different types of human cancer cells, but very few studies have been conducted on oral cancer. The study aimed to evaluate the anticancer potential of Andrographolide on an oral cancer cell line (KB) through in-silico network analysis and in vitro assays.

Materials and Methods:

The in-silico analysis involved the determination of drug-likeness prediction, prediction of common targets between oral cancer and andrographolide, Protein-Protein Interactions (PPI), hub genes, top 10 associated pathways by Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway, gene ontology (GO), and molecular docking experiments. In vitro assays comprised MTT assay, apoptosis assay, cell cycle analysis, intracellular reactive oxygen species (ROS) measurement, mitochondrial membrane potential (MMP), anti-migration activity, and gene expressions using polymerase chain reaction (PCR).

Results:

Fifteen common genes were obtained and were seen to be involved in cellular proliferation, regulation of apoptosis, migration of cells, regulation of MAPK cascade, and regulation of cell cycle. The most common genes involved in the top 10 pathways were MAPK1, MAPK8, MAPK14, and IL6 which were seen to be associated with the MAPK signaling pathway which may be the key pathway through which andrographolide may aid in treating oral cancer. In vitro assays showed anti-proliferative properties, late apoptosis, and anti-migratory properties.

Conclusion:

According to the results obtained, andrographolide has shown anticancer properties and has the potential to be used as a chemotherapeutic drug. The in-silico approach used in the present study can aid as a model for future research in developing efficient cancer treatments.

Key Words

Andrographolide Apoptosis
Cell cycle
Gene expression
In-silico analysis
Mitochondrial membrane-
potential
Oral cancer
==== Body
pmcIntroduction

Oral cancer (OC) has the unfortunate distinction of being the leading cancer in India, with an increasing incidence rate of 10.4% and a mortality rate of 9.3% for both men and women. On a global scale, it ranks as the eighth most prevalent cancer form (1). Despite the continued importance of surgical resection in cancer treatment, a major obstacle lies in the substantial risk of cancer recurrence observed in many patients, frequently occurring within a short period (2).

Andrographis paniculata, a herbal plant belonging to the Acanthaceae family, has been traditionally used in medicine for addressing a variety of health conditions. Widely cultivated in India, Thailand, and China (3), this plant’s leaves and stem are rich in the diterpene lactone known as andrographolide. A. paniculata showcases diverse bioactive properties, including anticancer effects (4), anti-inflammatory benefits (5), hepatoprotection (6), immunomodulation, and anticancer properties (7), as well as anti-infection properties (8). Through dichloromethane extraction, the primary bioactive compound andrographolide from A. paniculata has demonstrated the ability to hinder the proliferation of various human cancer cells, spanning a broad spectrum of cancer types (9). The primary bioactive compound found in the plant is andrographolide (C20H30O5), constituting 1.84% of the plant extract. This colorless, crystal-like compound possesses an intensely bitter taste and features a lactone structure, specifically a bicyclic diterpenoid lactone (10). Its molar mass is 350.455 g/mol and its melting point is 230-231 °C. Its solubility in water is limited. In contemporary contexts, andrographolide and its various derivatives have been documented to exhibit numerous pharmacological attributes. These include anti-inflammatory, hepatoprotective, anti-viral, neuroprotective, antioxidant, anti-fibrotic, anti-hyperglycemic, anti-tumor, anti-atherosclerotic, antimicrobial, and cardiovascular protective properties (11-16).

Andrographolide has the capability to traverse the blood-brain barrier (BBB) and exerts a potent anti-inflammatory influence on various types of leukocytes, such as T-cells, macrophages, and neutrophils (17), as well as endothelial cells (18). Beyond its anticancer attributes, numerous studies have suggested that andrographolide exhibits a comprehensive anti-inflammatory impact, including the inhibition of NFκB binding to DNA (19, 20).

Andrographolide has been observed to regulate various signaling pathways associated with cancer and angiogenesis, including PI3K/AKT/mTOR (21, 22), SRC/MAPKs/AP-1 (23), TLR4/NF-κB/MMP-9 (24), and VEGF/VEGFR2/AKT (25). Recent findings underscore that the administration of andrographolide in human cancer cells elevates apoptosis rates and impedes cell proliferation. While this effect has been extensively studied in numerous cancer cells, the specific mechanistic pathways through which andrographolide exerts its anti-tumorigenic actions in OC remain unclear and await comprehensive understanding.

To explore the potential anti-cancer attributes of Andrographolide, particularly concerning OC, the present study utilized a combination of in-silico network pharmacology and in vitro assays.

Materials and Methods

In-silico analysis

Druglikeness prediction

To assess the drug-likeness of andrographolide, Lipinski’s rule of five (RO5) was employed, for screening oral drugs in humans. Various parameters were evaluated. Andrographolide SMILES format, CC12CCC(C(C1CCC(=C)C2CC=C3C(COC3=O)O)(C)CO), was entered into the “SwissADME server (http://www.swissadme.ch)”, an online tool that calculates parameters like absorption, distribution, metabolism, excretion (ADME), oral bioavailability (OB), and drug-likeness (DL)(26).

Target proteins of andrographolide and oral cancer

To predict the targets of andrographolide, the “Swiss Target Prediction database (http://www.swisstargetprediction.ch/)”(27) was employed for predicting the corresponding genes associated with andrographolide. The DisGeNet database was utilized to obtain the genes related to OC. The common genes were subsequently identified and selected for further analysis.

Gene ontology and pathway enrichment

To perform GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis, we utilized two different tools: the Database for Annotation, Visualization, and Integrated Discovery (DAVID, https://david.ncifcrf.gov/, ver. 6.8) for GO analysis and the ShinyGO database (ShinyGO, http://bioinformatics.sdstate.edu/) for pathway enrichment analysis. DAVID is a versatile tool for annotating and interpreting gene lists, while ShinyGO specializes in GO and pathway enrichment analysis. KEGG is a comprehensive pathway database that provides graphical representations of biochemical pathways(28, 29). GO is a valuable resource for functional genomics, offering definitions and classifications of gene functions (30). To present and analyze the data, we generated bubble charts and histograms using the Bioinformatics cloud platform (http://www.bioinformatics.com.cn/), an online platform designed for data processing and visualization.

Protein-protein interaction analysis

Protein-protein interactions play a crucial role in biological processes and are essential for understanding the complex systems within living cells (31). To map the PPI network, the cluster of target genes was analyzed using the Search Tool for the Retrieval of Interacting Genes database (http://string-db.org/; version 11.5). The analysis focused on “Homo sapiens” as the species and a threshold of >0.9 was applied to ensure high-confidence information. Subsequently, the PPI network was constructed using Cytoscape (https://cytoscape.org/; version 3.9.1), a widely used bioinformatics software for data visualization and integration (32). To identify clusters or highly interconnected regions within the PPI network, the Cytoscape plugin cytoHubba (https://apps.cytoscape.org/apps/cytohubba; version 0.1) was employed. Proteins with the highest MNC (Maximum Neighborhood Component) level rankings were identified as hub targets within the network.

Molecular docking assessment between hub genes and andrographolide

The molecular docking simulations were executed using CB-Dock, a tool capable of automatically identifying active sites within a given protein, determining their centers and sizes, based on the query ligands (33). The Protein Data Bank (http://www.rcsb.org) was used to access the crystal structures of the target proteins. Similarly, the 3D structure of Andrographolide was obtained from the PubChem compound database (https://pubchem.ncbi.nlm.nih.gov/). These protein and ligand structures were used as inputs for CB-Dock, where the docking analysis explored the binding activities between the proteins and andrographolide. The Discovery Studio Visualizer software (Accelrys Software Inc.) was employed for the visualization and analysis of the docking results (34).

Top of form

Gene expression levels of hub genes

In this study, Gene Expression Profiling Interactive Analysis (GEPIA; http://gepia2.cancer-pku.cn/) was employed to verify the varied expressions of the hub genes in OC and normal oral tissues. GEPIA is an online server that offers interactive and customizable functionalities utilizing data from the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) database Furthermore, GEPIA facilitated the analysis of these genes based on pathological stages, providing valuable insights into their expression patterns in different disease stages (35).

Overall survival analysis of hub genes

To investigate the impact of the hub targets on the overall survival (OS) of patients with OC, the Kaplan-Meier [KM] Plotter (http://kmplot.com/analysis/index.php?p=service)(36), a cancer genomics dataset, was utilized. This dataset allows for the assessment of the prognostic significance of genes on survival outcomes. The patients with OC were categorized into two groups based on the presentation levels of the hub genes: low and high expression. A KM survival plot was generated to compare the survival outcomes between the two groups.

Anticancer activity of Andrographolide

Cell culture and maintenance

The KB (Oral Cancer) cell line was obtained from the National Centre for Cell Sciences (NCCS) in Pune. These cells were cultured in DMEM (GibcoTM) medium at a temperature of 37 °C. The culture medium was supplemented with 10% Fetal Bovine Serum (FBS, GibcoTM) and 1% antimycotic-antibiotic solution (GibcoTM). The cells were maintained in a CO2-enriched environment with a concentration of 5% to support their growth and viability.

Preparation of stock solutions

Andrographolide (Sigma-Aldrich) was prepared as a stock solution with a concentration of 10 mg/ml by dissolving 10 mg of andrographolide in 1 ml of Dimethyl Sulfoxide (DMSO). This stock solution was then stored at -80 °C. From the stock solution, a working solution with a concentration of 1 mg/ml was prepared by diluting the stock solution in a complete medium and passing it through a sterile 0.22 µm filter to ensure sterility. Dilutions of the compound were made in a complete culture medium to obtain concentrations of 110, 100, 90, 80, 70, and 60 µg/ml.

Cytotoxicity assay of andrographolide

The MTT conversion assay was employed to assess cell cytotoxicity. KB cells were seeded in a 96-well plate at a density of 5×104 cells/ml. Various concentrations of andrographolide, ranging from 60 µg/ml to 110 µg/ml, were applied to the cells. After treatment, 20 µl of 5 mg/ml MTT solution was added to each well and incubated for 4 hr. Subsequently, 100 µl of DMSO was added to each well to dissolve the formazan crystals. The absorbance at 570 nm was measured using the Multi-scanGo Thermo Fischer Scientific ELISA plate reader. The IC50 values were calculated using a Microsoft Office Excel worksheet (37).

Apoptosis analysis

The apoptosis assay was performed using the FITC Annexin V/Dead Cell Apoptosis Kit (Invitrogen-Molecular Probes®). The protocol prescribed by the manufacturer was followed. Flow cytometry was used to examine these stained cells at emission wavelengths of 530 and >575 nm. In this assay, Andrographolide and Paclitaxel (PTX) IC50 and IC25 values were used for performing the Apoptosis assay (38).

Cell cycle analysis

KB cells, both with and without treatment of andrographolide and PTX, were harvested 24 hr after washing with 1X PBS and trypsinization. To the collected cells, a mixture of 25 μl RNase A (20 mg/ml Invitrogen), 2 mM MgCl2 (Sigma), and 5-10 μl of 100 µg/ml propidium iodide (Invitrogen) was added. The cells were then incubated at room temperature for 10-15 min and subsequently analyzed using a FACS-caliber instrument from BD Bioscience (39).

Wound scratch assay

A density of 2×105 cells/ml was seeded in 12-well plates and allowed to reach over 90% confluency. A linear wound was created in the center of each well using a 200 μl plastic tip. The wounded cell monolayers were washed three times with 1X PBS to remove cell debris and then treated with andrographolide and PTX, with a control group, before being incubated for 24 hr. Photos of the scratches were recorded at 0, 6, and 24 hr using an OLYMPUS CKX53 microscope. Subsequently, the cells were allowed to migrate by incubating them at 37 °C in a medium containing 5% serum in the presence or absence of the drugs (40).

Estimation of generation of intracellular reactive oxygen species [ROS]

To assess ROS generation, flow cytometry was employed using DCFHDA(41). In this assay, KB cells were cultured in 6-well plates for 24 hr. Upon reaching 70-80% confluence, the cells were treated with andrographolide and PTX at concentrations corresponding to their respective IC50 and IC25. Flow cytometry [Beckman Coulter Cytomics FC 500 instrument] (495 nm and 520 nm) was employed to measure the fluorescence.

Mitochondrial membrane potential (ΔΨm):

To assess ΔΨm in KB cells, we employed the MitoProbe™ DiIC1 (5) Assay Kit as per the manufacturer’s instructions. This kit contains DiIC1 (5), a cyanine dye sensitive to changes in membrane potential, and CCCP, a disrupter of mitochondrial membrane potential used for research. DiIC1 (5) can readily enter cell cytoplasm and accumulate in mitochondria with active ΔΨm, yielding a bright far-red fluorescence.

RNA isolation and RT-PCR

The culture cells were treated with the desirable drug concentrations. After drug incubation cells were harvested by discarding the growth medium and total RNA was isolated from cells using the TRIzol method (TRIzol™ Reagent, Invitrogen cat no. 15596018), following the manufacturer’s instructions. The primers used in the present study are listed in Supplementary Table 1[ST: 1].

Statistical analysis

The experiments were replicated three times, and the results are presented as the mean±standard deviation (SD). Data analysis was done using GraphPad Prism 8. “Two-way ANOVA” followed by respective post hoc tests at ****P<0.0001, ***P<0.001, and **P<0.01 to determine the statistical significance between the groups as compared to control.

Table 1 Molecular docking scores for andrographolide with hub target proteins depicting binding energy levels indicating strong interaction of hub gene and andrographolide

Receptor	PDB ID	Binding energy (kcal/Mol)	
IL6	1P9M	-7.1	
MMP9	1L6J	-7.2	
MAPK1	1TVO	-9	
AR	1XQ3	-7.5	
CDK4	2W96	-8.1	
MAPK14	1A9U	-7	
MAPK8	3PZE	-7.8	
ADAM17	2DDF	-7.7	
WEE1	5VC5	-8.4	
CDK1	4YC6	-7.9	

Figure 1 Gene ontology (GO) and kyoto encyclopaedia of genes and genomes (KEGG) enrichment analysis

A): GO enrichment analysis was obtained by taking 10 biological processes (BP), Molecular functions (MF), and Cellular components (CC) from the DAVID Database for the targets in andrographolide treating oral cancer. B): A Bubble plot depicting the enrichment of the top 15 signaling pathways associated with Oral Cancer was generated. The X-axis represents the enrichment factor of the genes, while the Y-axis represents the different pathways. The circles in the plot are color-coded and sized based on the Log10 [P-value], where red signifies pathways with the highest number of genes and light green indicates pathways with fewer genes

Figure 2 Sketch and matching diagrams of molecular docking of Andrographolide and top 10 hub genes (target proteins)

Figure 3 mRNA expression levels of hub genes in the cancer genome atlas (TCGA) and genotype-tissue expression (GTEx) databases

(A) mRNA expression levels in the GEPIA database (Boxplot of hub genes). Red depicts Oral Cancer tissue and grey depicts normal oral mucosa. (B) mRNA expression level and pathological stage in the GEPIA database. (Stage plot of hub genes)

Figure 4 Kaplan–Meier overall survival analyses of patients with oral cancer based on expression of the ten hub genes. HR, hazard ratio (“http://kmplot.com/analysis/index.php?p=service&cancer=pancancer_rnaseq”)

Figure 5 Apoptosis (A-C)

A) Flow-cytometry assay of andrographolide-induced apoptosis with andrographolide at IC50 (106.2 µg/ml) and IC25 (53.1µg/ml) after 24 hr, (B) represents the cells treated with PTX at IC50 (92.27 µg/l) and IC25 (46.13µg/ml) after 24 hr. (C) The data was expressed as mean±standard deviation (SD), comparing apoptotic phases for the respective concentrations of andrographolide and PTX. 5: Cell Cycle (D-F). D) Assessment of the impact of andrographolide and PTX on cell cycle distribution with concentration of andrographolide IC50 (106.2 µg/ml) and IC25 (53.1 µg/ml). E) PTX IC50 (92.27 µg/l) and IC25 (46.13µg/ml) for 24 hr. F) The data was expressed as mean±standard deviation (SD); comparing different phases of the cell cycle for the respective concentrations of andrographolide and PTX

Figure 6 Representative Image for scratch assay

A) KB cells were incubated with andrographolide at IC50 (106.2 µg/ml) and PTX IC50 (92.21µg/ml) after 24 hr. The migration of KB cells was suppressed by andrographolide IC50 after 6 and 24 hr, as depicted in the second vertical panel. The first panel represented the control group. Additionally, the scratch area of KB cells treated with PTX IC50 was shown in the third panel. B) The data was expressed as Mean ± Standard Deviation (SD) of n=3 of independent experiments

Figure 7 mRNA Expression levels of KB cells (Real time qRT-PCR)

Comprehensive analysis of Bax, c-Myc, SLUG, and TGF-β genes in KB cells treated with andrographolide and PTX at IC50 and IC25 with medium only. The b-actin gene is used as a control. The results were subjected to Two-way ANOVA followed by Tukey post hoc tests at ****P<0.0001, ***P<0.001, and **P<0.01 to determine the statistical significance between the groups as compared to the control of respective time points

Results

In-silico analysis

Molecular properties of andrographolide

Our findings demonstrate that andrographolide adheres to Lipinski’s Rule of Five (RO5). The molecular properties of andrographolide, followed the RO5 criteria, indicating that it possesses favorable drug-like characteristics.

Target identification and analysis

Screening of OC-related targets using the search term “Lip and oral Cavity carcinoma” resulted in the identification of a total of 734 targets. Additionally, the Swiss Target Prediction database was employed to search for targets of andrographolide, leading to the identification of 100 targets. A comparison between the targets of andrographolide and the OC-related targets revealed a common set of 15 genes. (Supplementary Figure 1. SF:1).

Development of protein–protein interaction network (PPI) and determination of key targets

The study used the STRING Database to analyze protein-protein interactions among 15 identified targets. Cluster analysis in Cytoscape 3.9.1 revealed two distinct clusters: Cluster 1 (P-value-0.00297)(SF:2) associated with cell proliferation (red), apoptosis (yellow), cell migration (blue), and the MAPK cascade (green), and Cluster 2 (P-value-0.00889)(SF:2) linked to cell division, G2/M transition regulation (red), and G1/S transition regulation in the mitotic cell cycle (blue). These clusters suggest close associations and shared biological processes among the target proteins.

Top 10 hub genes analysis

Through the application of various algorithms, the study identified the top 10 hub genes: IL6, MMP9, MAPK1, AR, CDK4, MAPK14, MAPK8, ADAM17, WEE1, and CDK1. Among these hub genes, IL6 emerged as the most prominently active gene (ST: 2, SF:3).

GO and KEGG enrichment analysis

GO enrichment analysis of the 10 hub genes revealed around 69 GO terms. In terms of biological processes (BP), these targets are involved in cellular responses to ROS, regulation of the G2/M phase, and apoptotic processes. The cellular component (CC) results included the cyclin-dependent protein kinase, mitotic spindles, and the endoplasmic reticulum lumen. In relation to molecular function (MF), the targets predominantly play roles in IL6 receptor binding and MAP kinase activity (Figure 1). Through KEGG pathway analysis, a total of 40 pathways were identified. Among these pathways, the top 15 pathways were selected for further examination. It was found that the top 10 hub genes had close associations with various pathways, predominantly the TNF signaling pathway, IL-17 signaling pathway, and cellular senescence (Figure 1). The top 10 pathways consistently involved genes such as MAPK 1, 8, 14, and IL-6. These genes were notably associated with the MAPK signaling pathway (ST:3).

Confirmation of hub target by molecular docking

In a recent evaluation of drug-target interactions, ten hub genes were chosen as targets for molecular docking analysis. Using CB-DOCK, the structure of andrographolide was assessed for its docking potential with IL6, MMP9, MAPK1, AR, CDK4, MAPK14, MAPK8, ADAM17, WEE1, and CDK1. The results revealed binding energies lower than -5.0 for all core target proteins, indicating a robust binding activity between andrographolide and the core targets. The specific binding energies are presented in Table 1, and the docking sketch maps illustrating the interactions between the target proteins and andrographolide are depicted in Figure 2.

mRNA expression levels of hub genes

Using the GEPIA database, we found that CDK1, CDK4, and MMP9 mRNA levels were significantly higher in OC tissues compared to normal oral mucosa samples (P<0.01)(Figure 3A). Furthermore, results showed that CDK1, CDK4, IL6, MAPK8, and MMP9 exhibited significant changes across different pathological stages. CDK1, CDK4, IL6, and MAPK8 showed substantial increases in stage IV, while MMP9 increased in stage II (Figure 3B).

Survival analysis of the hub genes

Survival analysis was conducted on the 10 hub genes. The analysis was performed on a cohort of 500 OC patients from the TCGA database. The results demonstrated that all of the hub genes exhibited a significant association with poor prognosis (P<0.05)(Figure 4).

Anticancer activity of andrographolide on KB cell line:

Cytotoxicity result

Using different concentrations of Andrographolide (60, 70, 80, 90, 100, and 110 µg/ml) and PTX (40, 50, 75, 100, 150, and 200 µg/ml) the experiment was performed and the findings demonstrate that andrographolide exhibited a potent cytotoxic effect on the KB cell line in a dose-dependent manner, with an IC50 value of 106±1 µg/ml and an IC25 value of 53±1 µg/ml. Similarly, PTX displayed an IC50 value of 92±4.43 µg/ml and an IC25 value of 46.1±2.21 µg/ml (SF:4).

Effect of andrographolide on apoptosis regulation

Apoptosis analysis of KB cells was conducted with two concentrations of andrographolide: IC50 and IC25 for 24 hr. Results showed that exposure to the IC50 concentration led to 22.22±0.1 cells in the early stage of apoptosis and 37.48±3.1 cells in the late stage of apoptosis. Similar observations were made with the IC25 concentration of Andrographolide (Figure 5A). Whereas PTX, at its IC50, resulted predominantly in late stage of apoptosis. Similar observations were made at the IC25 value of PTX (Figure 5B). Additionally, when comparing the effects of andrographolide and PTX, it was noted that both compounds had a more pronounced impact on late apoptosis compared to early apoptosis (Figure 5C).

Effect of Andrographolide on cell cycle regulation of KB cells

Andrographolide at concentrations corresponding to IC50 and IC25 for 24 hr showed that approximately 46.7±0.8% of cells were arrested in the G1 phase, with 36.15±0.9% in the S phase at IC25 value. Similarly, at the IC50 concentration, around 7642±1.0% of cells were arrested in the G1 phase, and 16.8±0.6% were in the S phase (Figure 5D, 5E). After 24 hr of treatment with PTX at the IC50 concentration, there was an accumulation of 75.85±0.21 cells in the G1 phase and 16.75±0.92 cells in the S phase. Similarly, at the IC25 concentration, 46.35±0.07 cells were observed in the G1 phase and 21.4±4.38 cells in the S phase (Figure 5F). This finding underscores the specific impact of andrographolide on halting cell cycle progression during the G1 phase, shedding light on its potential implications for therapeutic interventions.

Effect of andrographolide on intracellular ROS measurement and mitochondrial membrane potential (ΔΨm) measurement

While Andrographolide at IC25 showed a slight ROS decrease in cancer cells, it wasn’t statistically significant (SF:5). We used potentiometric dyes in flow cytometry to detect early apoptosis stages marked by ΔΨm loss, but Andrographolide at IC50 and IC25, along with PTX, did not affect mitochondrial membrane potential in KB cells (SF:6). Thus, these treatments did not impact the mitochondrial membrane potential of KB cells.

Wound scratch assay results

After seeding the cells, a 24-hour incubation period was provided before treating them with IC50 concentrations of Andrographolide and PTX, alongside the control group. The wells were examined at 0, 6, and 24 hr to assess the effects of the treatments. The results indicated that the control group exhibited cell proliferation and migration after 6 hr, which was not observed in the group treated with Andrographolide and PTX (Figure 6A). The group treated with andrographolide and PTX showed significant inhibition of cell proliferation and migration at both 6 and 24 hr compared to the control group (Figure 6B).

Andrographolide-induced apoptosis in KB cells by targeting genes

Andrographolide caused an up-regulation of Bax gene expression predominantly by IC50 value, whereas IC25 did not show a significant effect on Bax expression. Whereas IC50 did not show any changes in SLUG expression. However, IC25 showed a slight increase in SLUG expression. indicating the induction of apoptosis in cancer cells and not affecting Epithelial-Mesenchymal Transition. Andrographolide induced apoptosis of OC cells by the TGF-β/c-myc pathway. At both IC-25 and IC-50 concentrations, andrographolide significantly reduced the expression of TGF-β and c-myc. IC50 value showed a significant reduction in c-myc expression compared to IC25 value. This observation suggests that andrographolide may inhibit c-myc expression by suppressing TGF-β, thereby regulating cellular proliferation and apoptosis (Figure 7).

Overall Andrographolide IC50 value exhibited a significant effect on the above-mentioned gene expression contributing to the anti-cancer activity of andrographolide against OC predominantly as an apoptosis promoter and regulation of cellular proliferation.

Discussion

Despite ongoing research and treatment advancements, clinical outcomes and overall survival rates for HNSCC have seen limited improvement in recent decades, with a discouraging 5-year survival rate as low as 50% (42, 43). Given the unsatisfactory results and significant toxicity associated with current treatment approaches for HNSCC, current research is concentrated on exploring alternative therapies with reduced toxicity. Complementary and alternative medicine has gained increasing attention as a promising area for cancer management, leading to a greater focus on exploring these options in recent years. In the recent past, over 3000 anti-cancer products derived from plants have been introduced, and what’s intriguing is that they tend to have significantly fewer side effects compared to conventional chemotherapy drugs (44).

An in-silico approach identified anti-cancer genes targeted by andrographolide. We assessed its cytotoxicity using the MTT assay and examined apoptotic potential through multiple assays in an OC cell line.

We identified 15 common genes associated with andrographolide and OC, impacting processes like cell proliferation, apoptosis regulation, cell migration, and the MAPK cascade. Notably, IL6, MMP9, MAPK1, MAPK14, and ADAM17 were prominent among the top 10 hub genes. Andrographolide demonstrated cytotoxic effects with an IC50 value of 106.2 µg/ml, inducing late apoptosis. This aligns with previous research where A. paniculata leaves’ methanol extract yielded four cytotoxic compounds, with andrographolide exhibiting the highest cytotoxicity and caspase-3 activation in HSC-2 cells (45).

Andrographolide was found to induce apoptosis in HT-29 cells, linked to increased intracellular ROS levels and disruption of ΔΨm. Interestingly, it caused G2/M phase cell cycle arrest at lower doses and G0/G1 phase arrest at higher doses (46). In contrast to some previous studies, we didn’t observe significant changes in ROS and MMP activity in our study. Andrographolide led to cell cycle arrest in the G0/G1 phase in our experiments.

Cancer cells often evade apoptosis by boosting anti-apoptotic BCL-2 proteins. Conversely, more resistant cancer types may down-regulate or inactivate pro-apoptotic proteins like Bax to suppress apoptosis. In our study, we noticed an up-regulation of Bax, suggesting that andrographolide exerts its anticancer effects by targeting the apoptotic pathway. Top of Form

Elevated c-Myc levels have been associated with reduced expression of immune checkpoints, which can suppress the immune response. In another study, inhibiting c-Myc expression was found to enhance apoptosis in CAL-27 cells (47). Another study by Marconi et al. emphasized the critical role of c-Myc in cell survival, proliferation, and tumor growth in Cal-27 cells (48). In accordance with these study results, our study also suggested that andrographolide may act as a negative regulator of c-Myc activity.

Conclusion

Our network pharmacology analysis demonstrates its impact on various targets, pathways, and biological processes, effectively regulating cell proliferation and apoptosis to combat OC. This indicates it inhibits cell viability, induces apoptosis, and suppresses TGF-β and downstream gene c-myc, highlighting its anti-OC activity. By inhibiting c-myc, it stabilizes the cell cycle, proliferation, DNA synthesis, and genomic stability, potentially enhancing chemotherapy efficacy against OC. It is a promising candidate for tumor therapies and a chemopreventive agent in human OC treatment. Nonetheless, further research is needed to validate its clinical effectiveness and comprehensively understand its mechanisms in the context of OC.

Acknowledgment

This work was financially supported by Dr D. Y. Patil Vidyapeeth, Pimpri, Pune [DPU/557 (4)/2020)]. The authors would like to thank Dr. Rojatkar for their assistance in the study.

Authors’ Contributions

G K provided methodology, data curation, software and formal analysis, and writing the original draft. A J helped with data curation, formal analysis, software, and methodology. S K supervised and helped review and edit. V L conceived the study and performed formal analysis, investigation, validation, and helped write the original daft. All authors approved the submitted manuscript.

Conflicts of Interest

All authors declare that they have no known competing financial interests.
==== Refs
References

1 Freddie Bray Jacques Ferlay Isabelle Soerjomataram Rebecca L Siegel Lindsey A Torre Ahmedin Jemal Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries CA Cancer J Clin 2018 68 394 424 30207593
2 Chen Zh Zhang P Xu Y Yan J Liu Z Lau WB Surgical stress and cancer progression: The twisted tango Mol Cancer 2019 18 132 142 31477121
3 Arvind S Negi J K Kumar Suaib Luqman Karuna Shanker M M Gupta S P S Khanuja Recent advances in plant hepatoprotectives: a chemical and biological profile of some important leads Med Res Rev 2008 28 746 772 17979145
4 L-J Wang X Zhou W Wang F Tang C-L Qi X Yang Andrographolide inhibits oral squamous cell carcinogenesis through NF-κB inactivation J Dent Res 2011 90 1246 1252 21841043
5 Haw-Wen Chen Ai-Hsuan Lin Hsing-Chin Chu Chien-Chun Li Chia-Wen Tsai Che-Yi Chao Inhibition of TNF-α-induced inflammation by andrographolide via down-regulation of the PI3K/Akt signaling pathway J Nat Prod 2011 74 2408 2413 22026410
6 Neha P Trivedi Upendra M Rawal Beena P Patel Potency of andrographolide as an antitumor compound in BHC-induced liver damage Integr Cancer Ther 2009 8 177 189 19679627
7 Wei Wang Jing Wang Sheng-fu Dong Chun-hong Liu Paola Italiani Shu-hui Sun Immunomodulatory activity of andrographolide on macrophage activation and specific antibody response Acta Pharmacol Sin 2010 31 191 201 20139902
8 Xiaojian Jiang Pei Yu Jie Jiang Zaijun Zhang Zhongli Wang Zhaoqi Yang Synthesis and evaluation of antibacterial activities of andrographolide analogues Eur J Med Chem 2009 44 2936 2943 19152987
9 R Ajaya Kumar K Sridevi N Vijaya Kumar S Nanduri S Rajagopal Anticancer and immunostimulatory compounds from Andrographis paniculata J Ethnopharmacol 2004 92 291 295 15138014
10 Bera R Ahmed S K Sarkar L Sen T Karmakar S Pharmacokinetic analysis and tissue distribution of andrographolide in rat by a validated LC-MS/MS method Pharm Biol 2014 52 321 329 24171780
11 Zhang H Li S Si Y Xu H Andrographolide and its derivatives: Current achievements and future perspectives Eur J Med Chem 2021 224 113710 34315039
12 Kumar G Singh D Tali JA Dheer D Shankar R Andrographolide: Chemical modification and its effect on biological activities Bioorg Chem 2020 95 103511 31884143
13 Zeng B Wei A Zhou Q Yuan M Lei K Liu Y Andrographolide: A review of its pharmacology, pharmacokinetics, toxicity and clinical trials and pharmaceutical researches Phytother Res 2022 36 336 364 34818697
14 Jayakumar T Hsieh CY Lee JJ Sheu JR Experimental and clinical pharmacology of Andrographis paniculata and its major bioactive phytoconstituent andrographolide Evid Based Complement Alternat Med 2013 2013 846740 846755 23634174
15 Islam MT Ali ES Uddin SJ Islam MA Shaw S Khan IN Andrographolide, a diterpene lactone from andrographis paniculata and its therapeutic promises in cancer Cancer Lett 2018 420 129 145 29408515
16 Kandanur SGS Tamang N Golakoti NR Nanduri S Andrographolide: A natural product template for the generation of structurally and biologically diverse diterpenes Eur J Med Chem 2019 176 513 533 31151068
17 Chiou WF Chen CF Lin JJ Mechanisms of suppression of inducible nitric oxide synthase (iNOS) expression in RAW 264 7 cells by andrographolide Br J Pharmacol 2000 129 1553 1560 10780958
18 Chang Ch Duann YF Yen TL Chen YY Jayakumar Th Ong ETh Andrographolide, a novel NF-κB inhibitor, inhibits vascular smooth muscle cell proliferation and cerebral endothelial cell inflammation Acta Cardiol Sin 2014 4 308 315
19 Xia YF Ye BQ Li YD Wang JG He XJ Lin X Andrographolide attenuates inflammation by inhibition of NF-kappa B activation through covalent modification of reduced cysteine 62 of p50 J Immunol 2004 173 4207 4217 15356172
20 Hidalgo MA Romero A Figueroa J Cortés P Concha II Hancke JL Andrographolide interferes with binding of nuclear factor-kappaB to DNA in HL-60-derived neutrophilic cells Br J Pharmacol 2005 144 680 686 15678086
21 Kumar D Das B Sen R Kundu P Manna A Sarkar A Andrographolide analogue induces apoptosis and autophagy mediated cell death in U937 cells by inhibition of PI3K/Akt/mTOR pathway PLoS One 2015 10 e0139657 26436418
22 Li Y Zhang P Qiu F Chen L Miao C Li J Inactivation of PI3K/Akt signaling mediates proliferation inhibition and G2/M phase arrest induced by andrographolide in human glioblastoma cells Life Sci 2012 90 962 967 22634579
23 Li Y Zhang P Qiu F Chen L Miao C Li J Andrographolide antagonizes TNF-α-induced IL-8 via inhibition of NADPH oxidase/ROS/NF-κB and Src/MAPKs/AP-1 axis in human colorectal cancer HCT116 cells J Agric Food Chem 2018 66 5139 5148 29672044
24 Zhang R Zhao J Xu J Jiao DX Wang J Gong ZhQ Andrographolide suppresses proliferation of human colon cancer SW620 cells through the TLR4/NF-κB/MMP-9 signaling pathway Oncol Lett 2017 14 4305 4310 28943944
25 Yang W Zhao J Wang Y Xu H Wu Zh Hu Y In vivo inhibitory activity of andrographolide derivative ADN-9 against liver cancer and its mechanisms involved in inhibition of tumor angiogenesis Toxicol Appl Pharmacol 2017 327 1 12 28438631
26 Daina A Michielin O Zoete V SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules Sci Rep 2017 7 42717 42729 28256516
27 Daina A Michielin O Zoete V SwissTargetPrediction: Updated data and new features for efficient prediction of protein targets of small molecules Nucleic Acids Res 2019 47 W357 W3664 31106366
28 Ladke VS Kumbhar GM Joshi K Kheur S Systemic explanation of Glycyrrhiza glabra’s analyzed compounds and anti-cancer mechanism based on network pharmacology in oral cancer J Oral Biosci 2022 64 452 460 36113760
29 Kanehisa M Furumichi M Sato Y Ishiguro-Watanabe M Tanabe M KEGG: Integrating viruses and cellular organisms Nucleic Acids Res 2021 49 D545 D551 33125081
30 M Ashburner C A Ball J A Blake D Botstein H Butler J M Cherry Gene ontology: Tool for the unification of biology The Gene Ontology Consortium. Nat Genet 2000 25 25 29 10802651
31 Lin JSh Lai EM Protein-protein interactions: Co-immunoprecipitation Methods Mol Biol 2017 1615 211 219 28667615
32 Shannon P Markiel A Ozier O Baliga NS Wang JT Ramage D Cytoscape: A software environment for integrated models of biomolecular interaction networks Genome Res 2003 13 2498 2504 14597658
33 Liu Y Grimm M Dai WT Hou MCH Xiao ZHX Cao Y CB-Dock: A web server for cavity detection-guided protein-ligand blind docking Acta Pharmacol Sin 2020 41 138 144 31263275
34 Yadav RP Ibrahim KS Gurusubramanian G Kumar NS In silico docking studies of non-azadirachtin limonoids against ecdysone receptor of Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) Med Chem Res 2015 24 2621 2631
35 Tang Z Li Ch Kang B Gao G Li Ch Zhang Z GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses Nucleic Acids Res 2017 45 W98 W102 28407145
36 Nagy Á Munkácsy G Győrffy B Pancancer survival analysis of cancer hallmark genes Sci Rep 2021 11 6047 6056 33723286
37 Javir G Joshi K Evaluation of the combinatorial effect of Tinospora cordifolia and Zingiber officinale on human breast cancer cells 3 Biotech 2019 9 428 439
38 Wyrębska A Gach K Lewandowska U Szewczyk K Hrabec E Modranka J Anticancer activity of new synthetic α-Methylene-δ-Lactones on two breast cancer cell lines Basic Clin Pharmacol Toxicol 2013 113 391 400 24034139
39 Sanap A Chandravanshi B Shah T Tillu G Dhanushkodi A Bhonde R Herbal pre-conditioning induces proliferation and delays senescence in Wharton’s jelly mesenchymal stem cells Biomed Pharmacother 2017 93 772 778 28724259
40 Liang Ch Park AY Guan JL In vitro scratch assay: A convenient and inexpensive method for analysis of cell migration in vitro Nat Protoc 2007 2 329 333 17406593
41 Brandt R Keston AS Synthesis of diacetyldichlorofluorescin: a stable reagent for fluorometric analysis Anal Biochem 1965 11 6 9 14328648
42 Argiris A Brockstein BE Haraf DJ Stenson KM Mittal BB Kies MS Competing causes of death and second primary tumors in patients with locoregionally advanced head and neck cancer treated with chemoradiotherapy Clin Cancer Res 2004 10 1956 1962 15041712
43 Vokes EE Weichselbaum RR Lippman SM Hong WK Head and neck cancer N Engl J Med 1993 328 184 194 8417385
44 Koehn FE Carter GT The evolving role of natural products in drug discovery Nat Rev Drug Discov 2005 4 206 220 15729362
45 Suzuki R Matsushima Y Okudaira N Sakagami H Shirataki Y Cytotoxic components against human oral squamous cell carcinoma isolated from andrographis paniculata Anticancer Res 2016 36 5931 5935 27793918
46 Khan I Khan F Farooqui A Ansari IA Andrographolide exhibits anticancer potential against human colon cancer cells by inducing cell cycle arrest and programmed cell death via augmentation of intracellular reactive oxygen species level Nutr Cancer 2018 70 787 803 29781715
47 Zhao S An L Yang X Wei Zh Zhang H Wang Y Identification and validation of the role of c-Myc in head and neck squamous cell carcinoma Front Oncol 2022 12 820587 820600 36119473
48 Marconi GD Rocca YD Fonticoli L Melfi F Rajan ThS C-Myc expression in oral squamous cell carcinoma: Molecular mechanisms in cell survival and cancer progression Pharmaceuticals (Basel) 2022 15 890 903 35890188
