The Essential Oil from Conyza bonariensis (L.) Cronquist (Asteraceae) Exerts an In Vitro Antimelanoma Effect by Inducing Apoptosis and Modulating the MAPKs, NF-κB, and PKB/AKT Signaling Pathways

The characterization and cytotoxicity of the essential oil from Conyza bonariensis (L.) aerial parts (CBEO) were previously conducted. The major compound was (Z)-2-lachnophyllum ester (EZ), and CBEO exhibited significant ROS-dependent cytotoxicity in the melanoma cell line SK-MEL-28. Herein, we employed the Molegro Virtual Docker v.6.0.1 software to investigate the interactions between the EZ and Mitogen-Activated Protein Kinases (MAPKs), the Nuclear Factor kappa B (NF-κB), and the Protein Kinase B (PKB/AKT). Additionally, in vitro assays were performed in SK-MEL-28 cells to assess the effect of CBEO on the cell cycle, apoptosis, and these signaling pathways by flow cytometry and the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay using MAPKs inhibitors. CBEO induced a significant increase in the sub-G1 peak, as well as biochemical and morphological changes characteristic of apoptosis. The in-silico results indicated that EZ interacts with Extracellular Signal-Regulated Kinase 1 (ERK1), c-Jun N-terminal Kinase 1 (JNK1), p38α MAPK, NF-κB, and PKB/AKT. Moreover, CBEO modulated the ERK1/2, JNK, p38 MAPK, NF-κB, and PKB/AKT activities in SK-MEL-28 cells. Furthermore, CBEO’s cytotoxicity against SK-MEL-28 cells was significantly altered in the presence of MAPKs inhibitors. These findings support the in vitro antimelanoma effect of CBEO through apoptosis induction, and the modulation of ERK, JNK, p38 MAPK, NF-κB, and PKB/AKT activities.


Introduction
Cutaneous melanoma is the most aggressive type of skin cancer [1], originating from the malignant transformation of melanocytes [2].This cancer represents a global public health problem, with 57,000 recorded deaths in 2020 [3].
Chemotherapy remains a significant pharmacological modality for the treatment of many cancer types [4], including melanoma [5].However, issues are associated with antineoplastic chemotherapy, such as chemoresistance [6][7][8].Therefore, research should be conducted to obtain new therapeutic agents [9].In this context, natural products standout as a valuable source of molecules with potential applications for cancer treatment [10].
be conducted to obtain new therapeutic agents [9].In this context, natural products standout as a valuable source of molecules with potential applications for cancer treatment [10].
Apoptosis, or programmed cell death type I, involves a series of tightly controlled events.Cells undergoing apoptosis exhibit biochemical and morphological characteristics, such as DNA condensation and fragmentation, and the exposure of phosphatidylserine on the outer surface of the plasma membrane [24].During carcinogenesis, cells acquire the ability to evade apoptosis, resulting in proliferative and survival advantages [25].
Mitogen-activated protein kinases (MAPKs) include Extracellular Signal-Regulated Kinases 1 and 2 (ERK1 and ERK2); c-Jun N-terminal Kinases 1, 2, and 3 (JNK1, JNK2, and JNK 3); and p38 Mitogen-Activated Protein Kinases α, β, γ, and δ (p38α, p38β, p38γ, and p38δ) [26].Regarding the involvement of the MAPKs' signaling pathways in cancer, it is known that approximately 50% of melanomas exhibit mutations in the BRAF oncogene, leading to constitutive activation of the Extracellular Signal-Regulated Kinases 1 and 2 (ERK1/2) pathway and uncontrolled proliferation [27].Zhou et al. (2021) [28] observed that the expression of utrophin, a protein encoded by the tumor suppressor gene UTRN, is reduced in the human melanoma cell line A375, and its antiproliferative activity in these cells involves the inhibition of the c-Jun N-terminal Protein Kinase (JNK) pathway.In BRAF mutant melanoma cells (A375 and SK-MEL-28), the activation of p38 MAPK is associated with the negative regulation of subtype 4b plasma membrane Ca 2+ pumps The antimelanoma effect of various EOs and some of their components occurs through multiple mechanisms, including cell cycle arrest and the induction of apoptosis [20].Additionally, it has been demonstrated that the increase in ROS is involved in the antitumor effect by modulating intracellular signaling pathways, including the Mitogen-Activated Protein Kinases (MAPKs) [21], Nuclear Factor kappa B (NF-κB) [22], and Protein Kinase B (PKB/AKT) [23].
Apoptosis, or programmed cell death type I, involves a series of tightly controlled events.Cells undergoing apoptosis exhibit biochemical and morphological characteristics, such as DNA condensation and fragmentation, and the exposure of phosphatidylserine on the outer surface of the plasma membrane [24].During carcinogenesis, cells acquire the ability to evade apoptosis, resulting in proliferative and survival advantages [25].
Mitogen-activated protein kinases (MAPKs) include Extracellular Signal-Regulated Kinases 1 and 2 (ERK1 and ERK2); c-Jun N-terminal Kinases 1, 2, and 3 (JNK1, JNK2, and JNK 3); and p38 Mitogen-Activated Protein Kinases α, β, γ, and δ (p38α, p38β, p38γ, and p38δ) [26].Regarding the involvement of the MAPKs' signaling pathways in cancer, it is known that approximately 50% of melanomas exhibit mutations in the BRAF oncogene, leading to constitutive activation of the Extracellular Signal-Regulated Kinases 1 and 2 (ERK1/2) pathway and uncontrolled proliferation [27].Zhou et al. (2021) [28] observed that the expression of utrophin, a protein encoded by the tumor suppressor gene UTRN, is reduced in the human melanoma cell line A375, and its antiproliferative activity in these cells involves the inhibition of the c-Jun N-terminal Protein Kinase (JNK) pathway.In BRAF mutant melanoma cells (A375 and SK-MEL-28), the activation of p38 MAPK is associated with the negative regulation of subtype 4b plasma membrane Ca 2+ pumps (PMCA4b), a metastasis suppressor [29].However, it has been reported that these intracellular signaling pathways have pleiotropic roles.In fact, several clinically used antineoplastic drugs, such as doxorubicin, induce the activation of ERK1/2 and p38 MAPK [30][31][32].Additionally, in ovarian carcinoma cells, the activation of autophagy-dependent cell death can be induced by the positive modulation of JNK [33].
The PKB/AKT protein signaling pathway plays essential roles in cellular physiology, including growth, survival, proliferation, and metabolism [40].Various diseases are associated with the dysregulation of this signaling pathway, such as diabetes and cardiovascular and neurological disorders [41].In cancer, the overexpression and/or aberrant activation of proteins in this pathway are linked to the development of various malignant tumors.Constitutive activation of PKB/AKT is present in approximately 70% of melanomas [42,43].Conversely, Sun et al. (2019) [44] demonstrated the involvement of the activation of the PKB/AKT pathway in the antitumor effect through the induction of apoptosis by aloin and metformin, alone or in combination, in liver cancer models.
Therefore, considering that cell cycle arrest, apoptosis induction, and the modulation of the MAPK, NF-κB, and PKB/AKT signaling pathways are promising targets for cancer treatment, we present here the in vitro effects of CBEO on the cell cycle and the induction of apoptosis.Additionally, we investigated the molecular docking prediction between the major compound of CBEO, EZ, and the proteins ERK1, JNK1, p38α MAPK, NF-κB (p50/p65), and PKB/AKT.Finally, our study provides insight into the in vitro effect of CBEO on the modulation of MAPK, NF-κB, and PKB/AKT proteins in the human melanoma cell line SK-MEL-28.
plastic drugs, such as doxorubicin, induce the activation of ERK1/2 and p38 MAPK [30][31][32].Additionally, in ovarian carcinoma cells, the activation of autophagy-dependent cell death can be induced by the positive modulation of JNK [33].
The NF-κB family consists of five transcription factors that form homo-or heterodimers: c-Rel, RelA (p65), RelB, NF-κB1 (p105/p50), and NF-κB2 (p100/p52) [34], all possessing a Rel homology domain (RHD) for DNA binding [35].In cancer, the NF-κB signaling pathway contributes to tumor progression by activating the transcription of multiple genes, such as those involved in apoptosis inhibition and angiogenesis induction [36].On the other hand, it has also been reported that the cytotoxicity of some antitumor drugs depends on the activation of this pathway [37][38][39].
The PKB/AKT protein signaling pathway plays essential roles in cellular physiology, including growth, survival, proliferation, and metabolism [40].Various diseases are associated with the dysregulation of this signaling pathway, such as diabetes and cardiovascular and neurological disorders [41].In cancer, the overexpression and/or aberrant activation of proteins in this pathway are linked to the development of various malignant tumors.Constitutive activation of PKB/AKT is present in approximately 70% of melanomas [42,43].Conversely, Sun et al. (2019) [44] demonstrated the involvement of the activation of the PKB/AKT pathway in the antitumor effect through the induction of apoptosis by aloin and metformin, alone or in combination, in liver cancer models.
Therefore, considering that cell cycle arrest, apoptosis induction, and the modulation of the MAPK, NF-κB, and PKB/AKT signaling pathways are promising targets for cancer treatment, we present here the in vitro effects of CBEO on the cell cycle and the induction of apoptosis.Additionally, we investigated the molecular docking prediction between the major compound of CBEO, EZ, and the proteins ERK1, JNK1, p38α MAPK, NF-κB (p50/p65), and PKB/AKT.Finally, our study provides insight into the in vitro effect of CBEO on the modulation of MAPK, NF-κB, and PKB/AKT proteins in the human melanoma cell line SK-MEL-28.
Before carrying out the molecular docking simulations, the validation of the compounds under study was carried out through redocking between the ligands and cocrystallized proteins.In redocking, the RMSD (Root Mean Square Deviation) value is determined.This value is calculated by comparing the coordinates of the heaviest atoms in the experimentally determined crystallographic structure with their counterparts in the docked position.An RMSD value of up to 2 Å is considered acceptable for this metric.
The redocking calculation was performed only for enzymes that present a cocrystallized ligand and corresponded to 1.2477 for the ERK1 in the complex with the ligand carbamate (PDB: 5LCJ) and 0.1135 for the JNK1 in the complex with the ligand 6-chloro-9-hydroxy-1,3-dimethyl-1,9-dihydro-4H-pyrazolo (3,4-b)quinolin-4-one (PDB: 2G01).During the redocking analysis, it was observed that the RMSD value was less than 2.0 Å, indicating that the generated pose correctly positioned the ligand within the active site.This indicates that the program provided values considered satisfactory for the docking validation.The docking results for the targets under study, according to the energy score of the MolDock Score algorithm, can be viewed in Table 1.
As observed in Figure 8, CBEO treatment induced an increase in the percentage of cells labeled with anti-JNK antibodies (40 µg/mL: 61.20 ± 1.90%, p < 0.05), compared to the control (2.82 ± 0.27%).The standard drug DXR also induced a significant increase in the percentage of cells labeled with anti-JNK antibodies (16.19 ± 1.33%, p < 0.05), compared to the control.

Discussion
Natural products represent many compounds currently used in antineoplastic chemotherapy.[45].In this context, research has been conducted to investigate the antitumor effects of natural products, such as essential oils (EOs) [46].The EO from Conyza bonariensis (L.) Cronquist aerial parts (CBEO) was previously obtained and characterized by Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences among conditions; p < 0.05.

Discussion
Natural products represent many compounds currently used in antineoplastic chemotherapy.[45].In this context, research has been conducted to investigate the antitumor effects of natural products, such as essential oils (EOs) [46].The EO from Conyza bonariensis (L.) Cronquist aerial parts (CBEO) was previously obtained and characterized by Ferreira et al. (2023) [19].The major compound identified was (Z)-2-lachnophyllum ester (EZ) (57.24%).CBEO exhibited toxicity in the zebrafish model, cytotoxicity in human peripheral blood mononuclear cells (PBMCs), and in vitro antitumor effects in the human melanoma cell line SK-MEL-28 by inducing oxidative stress.To better characterize the CBEO antimelanoma effect, we investigated its effect on cell cycle arrest, the induction of apoptosis, and the modulation of the Mitogen-Activated Protein Kinases (MAPKs), Nuclear Factor kappa B (NF-κB), and Protein Kinase B (PKB/AKT) signaling pathways through in silico and in vitro assays using the SK-MEL-28 cell line.
In our study, doxorubicin (DXR) was used as the standard drug.DXR is an anthracycline isolated from the species Streptomyces peucetius var.caesius, and it is still widely employed in antineoplastic chemotherapy, including the treatment of leukemias and breast cancer [47].Its antitumor mechanism of action includes the inhibition of the enzyme topoisomerase II, the induction of cell cycle arrest, the induction of oxidative stress [48], and the modulation of the activity of Mitogen-Activated Protein Kinases (MAPKs) [49][50][51], nuclear factor kappa B (NF-κB) [52], and protein kinase B (PKB/AKT) [53].
The cell cycle is a tightly regulated process that leads to cell division.In cancer, this event becomes dysregulated, resulting in uncontrolled cell proliferation.Therefore, agents that interfere with the cell cycle are considered promising for antitumor chemotherapy [54].The significant increase in the sub-G1 peak observed after CBEO treatment is characteristic of apoptotic cells [55].This effect has also been observed after the exposure of tumor cells to other EOs [56][57][58].Therefore, considering that the increase in the sub-G1 peak is related to apoptosis, the possible induction of apoptosis by CBEO was also investigated.
Apoptosis is a regulated process of cell death [59], characterized by cellular changes such as protein degradation, membrane blebbing, and DNA cleavage [60].In cancer, cells exhibit resistance to apoptosis [61], and thus, the induction of this type of cell death has been a major goal in antitumor therapy [62].CBEO induced an increase in the percentage of cells undergoing apoptosis.In addition, the analysis of the images obtained by laser confocal microscopy showed morphological changes characteristic of the apoptotic process, such as chromatin condensation, membrane blebs, and DNA fragmentation.Apoptosis is a shared mechanism of action for various EOs of the Asteraceae family.The EO from Tridax procumbens leaves induced apoptosis in murine melanoma B16-F10 cells and in the in vivo melanoma model [63].Eos from Artemisia arborescens flowers and leaves were able to induce the apoptotic process [64].Carpesium abrotanoides' EO induces apoptosis via the intrinsic pathway in human hepatocellular carcinoma cells (HepG2) [65].Finally, the Conyza canadensis EO induced apoptosis in the human cervical cancer HeLa cell line [66].Based on our results, CBEO induces an antimelanoma effect by apoptosis induction.
In melanoma, tumoral development is associated with oxidative stress, and Reactive Oxygen Species (ROS) are involved in the modulation of intracellular signaling pathways related to cell proliferation and apoptosis [21].The ROS increase is involved in the antitumor effect of several drugs by modulating intracellular signaling pathways, including the MAPKs pathway [67], NF-κB pathway [22], and PKB/AKT pathway [23].Ferreira et al. (2023) [19] observed that the CBEO antimelanoma effect is ROS-dependent.Therefore, as a preliminary analysis of the possible involvement of these intracellular signaling pathways in the antitumor effect of this EO, molecular docking was performed on crystallographic structures of these proteins to obtain protein-ligand complexes with the major compound, EZ, of CBEO.
MAPKs' signaling pathways are involved in various cellular processes, including proliferation and survival [68].In cancer, the dual role of these proteins has been reported.Therefore, stimulating the activity of these kinases can be involved in either a pro-tumoral or anti-tumoral effect [69,70].A molecular docking analysis showed favorable interaction between EZ and the ERK1 protein through amino acid residues that also contribute to the interaction of this target with the standard drug DXR (Ala35) or the ligand from the Protein Data Bank (PDB) (Ala35 and Leu170).We did not find molecular docking data between ERK1 compounds structurally similar to EZ.Nevertheless, the literature reports such data for compounds from different chemical classes.Nagalakshmamma et al. (2021) [71] designed and synthesized compounds derived from 1,4-diisocyanatobenzene bisurea, namely 3a-3j.These compounds exhibited a significant in vitro antitumor effect against the HeLa cell line.In the molecular docking study, it was observed that the derivatives 3d and 3e, which were the most cytotoxic against the cell line used, had binding with the ERK protein through amino acid residues, including Ala35 and Leu170.EZ was also able to interact with the JNK1 protein for amino acid residues that are involved in the binding of DXR (Ala53, Ile32, and Val158) or the PDB ligand (Ala53, Ile32, Val158, and Met111) to this target.In addition, EZ showed a higher binding probability value compared to the PDB ligand.We did not find reports in the literature of in silico molecular interaction between JNK1 and compounds structurally similar to EZ. Sugara et al. (2021) [72] showed new haloxanthones, coded as 3CX, 10CX, 5BX, and 4BX, designed for potential antitumor activity and bound to JNK, involved amino acids that participated in the interaction between EZ and this protein (Ala53, Ile32, Val158, Leu110, Met111, and Glu109).Favorable binding was also observed between EZ and the p38α MAPK in Gln120 amino acid residue that is involved in the interaction of DXR with this protein.Similarly, Leu216 and Val117 involved in the binding of EZ are also involved in the interaction between p38α MAPK and (4-(4-fluorophenyl)-2-(4-methylsulfinylphenyl)-5-(4-pyridyl)-imidazole) (SB203580), a specific inhibitor of this protein [73].We did not find molecular docking data between p38α MAPK and compounds structurally similar to EZ. Shankar et al. (2021) [74] investigated the interaction between agasthisflavone, a phytochemical compound found in Anacardium occidentale (Anacardiaceae), and p38α MAPK.These researchers observed that there is binding between agasthisflavone and p38α MAPK through the interaction with Pro153, which also participates in the binding between EZ and this protein.
NF-κB is a transcription factor that plays important roles in cellular functions such as proliferation, differentiation, and apoptosis [75].Currently, there is substantial evidence pointing to the dual role of NF-κB in cancer [76].There was a favorable interaction between EZ and NF-κB, with amino acid residues involved in this binding also participating in the interaction between DXR and NF-κB (Cys120, Arg187, and Tyr36).EZ showed a higher binding affinity to NF-κB when compared to the binding affinity of the ligand 4-methyl-N1-(3-phenylpropyl)benzene-1,2-diamine to this protein.We did not find molecular docking data between NF-κB and compounds structurally similar to EZ.Nevertheless, the literature reports such data for compounds from different chemical classes.Shiroma et al. (2020) [77], after a virtual screening, identified compound A55 (2-(3-carbamoyl-6-hydroxy-4-methyl-2oxopyridin-1(2H)-yl)acetic acid) as a selective inhibitor of RelA(p65)-DNA binding.The molecular docking result showed that this compound interacted with the Tyr36 residue of NF-κB, which is also involved in the interaction between EZ and this protein target.Zazeri et al. (2020) [78], in the search for structural analogs of piperine and piperlongumine (PPL), which are alkaloids found in Piper longum [79] and known for their various biological activities including antitumor effects [80,81], obtained a series of new substances.In the in silico study, in addition to the natural compounds piperine and PPL, the synthesized compounds exhibited an interaction with NF-κB.The interaction between PPL and NF-κB involved all the residues that are also related to the interaction of EZ and this protein (Cys120, Arg187, Tyr36, Leu154, Phe34, and Asp185).It is known that the antitumor activity of PPL involves the modulation of NF-κB [82].On the other hand, among the amino acids involved in the interaction between EZ and NF-κB, the interaction of piperine with this protein involves only the Arg187 residue.Among the analogs, compound 4a stood out, and its interaction with the molecular target involved residues Cys120, Arg187, Tyr36, Leu154, and Asp185.
The PKB/AKT protein regulates various cellular events, including cell proliferation and angiogenesis [83].In cancer, its effect has been reported to be either pro-tumoral or antitumoral [44,84].For this target, EZ showed significant interaction energy values, indicating a favorable binding.This interaction involved amino acid residues that also participate in the interaction between the ligand 8-[4-(1-aminocyclobutyl)phenyl]-9-phenyl-2H- [1,2,4]triazolo [3,4-f][1,6]naphthyridin-3-one (MK-2206) and this protein (Tyr316, Leu348, Leu317, Val331, Ala318, Lys277, Val338, and Glu279).As with the other molecular targets in this study, we did not find reports in the literature of an in silico molecular interaction between PKB/AKT and compounds structurally similar to EZ. Brevilin A, a molecule isolated from Centipeda minima (Asteraceae) with significant antitumor activity, interacts with PKB/AKT in amino acid residues such as Lys277 [85], which also belongs to the group of amino acids involved in the interaction between EZ and this target.This amino acid residue is also involved in the interaction between epirubicin, an anthracycline used in the treatment of cancers such as breast cancer [86], and PKB/AKT [87].Another compound with significant antitumor activity, pelargonidin, found in red or pink fruits, interacts with PKB/AKT in residues that include Leu317 and Lys277 [88].
Therefore, our results demonstrate that the proteins ERK, JNK, p38 MAPK, NF-κB, and PKB/AKT are potential targets for investigating the antitumor mechanisms of CBEO.Thus, in vitro tests were conducted to assess the effect of this EO on the signaling pathways of these proteins.
SK-MEL-28 cell treatment with CBEO induced ERK and JNK activation, as evidenced by the increased percentage of fluorescent cells marked with anti-ERK and anti-JNK antibodies.Furthermore, corroborating that the CBEO antitumor effect involves the activation of these intracellular signaling pathways, the ERK inhibitor (U0126) and JNK inhibitor (SP600125) treatments partially prevented the cytotoxicity of this EO.This antitumor mechanism is also shared with other EOs.For instance, after treating human oral squamous carcinoma cells (KB) with the EO from Artemisia lavandulaefolia (Asteraceae) aerial parts, the induction of apoptosis was observed, related to the phosphorylation of ERK and the JNK increased [89].Similarly, the EO from Pinus densiflora (Pinaceae) leaves induced apoptosis in human oral carcinoma cells (YD-8), involving ERK and JNK activation, ROS generation, and caspase activation [90].Additionally, the EO from Teucrium alopecurus (Poaceae) leaves induced cytotoxicity in human colon cancer cells (HCT-116) through ERK and JNK activation [91], which were also observed in the in vitro antitumor effect of the EO from Thymus hirtus sp.algeriensis (Lamiaceae) aerial parts on this cell line [92].Hence, when we consider that ROS have the potential to trigger ERK1/2 [32] and JNK activation [93], we can propose a connection between the apoptosis observed in SK-MEL-28 cells following CBEO treatment and the stimulation of ERK1/2 and JNK proteins.Furthermore, it appears that these proteins are influenced positively by the ROS generated as a result of EO treatment.
Regarding the p38 MAPK, CBEO reduced the percentage of cells marked with anti-p38 MAPK antibodies, and in the MTT assay, the p38 MAPK inhibitor (PD 169316) treatment significantly increased the CBEO cytotoxicity.The inhibition of p38 MAPK activity, which is involved in the antitumor effect of Eos, has also been reported in the literature.The antitumor effect of bergamot EOs, obtained from the Citrus bergamia (Rutaceae) fruit peel, on human neuroblastoma cell line SH-SY5Y, involves p38 MAPK inhibition [94].Similarly, the EO from Saussurea lappa (Asteraceae) roots induced apoptosis in human hepatocellular carcinoma cell lines (SMMC-7721 and Hep3B) by concentration-dependent p38 MAPK inhibition [95].Therefore, taking into account the dual role of p38 MAPK in cancer, we propose that the in vitro antimelanoma effect of CBEO is linked to its inactivation.Additionally, we suggest that the activation of JNK exerts a negative regulatory impact on p38 MAPK, as described in the literature [96].
The labeling of SK-MEL-28 cells with anti-NF-κB (p65) antibodies showed that this EO induces the activation of this protein.A similar effect was observed in KB cells after treatment with the EO from Artemisia capillaris (Asteraceae).In this study, cytosolic IκBα degradation increased nuclear translocation of the RelA (p65) protein as well as NF-κB binding to DNA [97].Furthermore, isolated molecules from EOs have also demonstrated antitumor effects by inducing NF-κB activation.The limonene antitumor effect, one of the monoterpenes found in CBEO, in the murine lymphoma cell line BW5147, involved the increased of NF-κB nuclear translocation [98].Additionally, patchouli alcohol isolated from the EO of Pogostemon cablin (Lamiaceae) exhibited significant cytotoxicity against HCT-116 and SW480 cells (human colon adenocarcinoma) by increasing NF-κB transcriptional activity.Moreover, a significant increase in p65 levels in the nucleus of HCT-116 cells was observed, confirming the involvement of the NF-κB pathway activation in the patchouli alcohol antitumor effect [99].Therefore, it is conceivable that the activation of ROS-dependent apoptosis following the treatment of SK-MEL-28 cells with CBEO is contingent upon the activation of NF-κB.Furthermore, ROS may be associated with the activation of this protein, as described in the literature [100].
Finally, PKB/AKT activation was evident after CBEO treatment.This antitumor action mechanism was also observed for different fractions of the EO from Boswellia sacra (Burseraceae) gum resin in Panc-28 cells (human pancreatic adenocarcinoma).In this study, it was observed that the expression of phosphorylated PKB/AKT increased rapidly, peaking at 15 min after the EO from Boswellia sacra (Burseraceae) gum resin treatment.Subsequently, gradual decreases in this phosphorylated PKB/AKT expression were observed [101].It has been demonstrated that the activation of PKB/AKT can be mediated by ROS [102].Therefore, we propose that the activation of PKB/AKT is associated with the induction of oxidative stress induced by CBEO.
The drugs and reagent solutions were prepared immediately before use.

Human Tumor Cell Line
The SK-MEL-28 (human melanoma) cell line was obtained from Rio de Janeiro Cell Bank (BCRJ), Brazil, and cultured in Dulbecco's Modified Eagle's Medium (DMEM) supplemented with 10% Fetal Bovine Serum and 1% penicillin-streptomycin at 37 • C with 5% CO 2 .

Cell Cycle Analysis
To assess the effect of CBEO on cell cycle progression, SK-MEL-28 cells were stained with PI (propidium iodide) for the DNA content analysis and analyzed by flow cytometry.The fluorescence intensity emitted by PI bound to DNA allows for the quantification of DNA in cells at each phase of the cell cycle [103].For this purpose, SK-MEL-28 cells were plated (2 × 10 5 cells/mL) in 24-well plates and incubated with CBEO at concentrations of 20 or 40 µg/mL (corresponding to IC 50 and double the CI 50 , respectively) or DXR (4 µM) for 48 h.After this period, the cells were collected, centrifuged (500× g, 20 • C, 5 min), gently resuspended in PBS, fixed with pre-chilled 70% ethanol while being agitated, and frozen (−20 • C) until the analysis.At the time of the analysis, the cells were recovered by centrifugation (400× g, 10 min, 4 • C), washed with PBS, and incubated with RNase (0.1 mg/mL) and propidium iodide (PI; 0.05 mg/mL) in the dark (25 • C, 30 min).Subsequently, the analysis was conducted using a flow cytometer, acquiring 10,000 events/sample.The percentages of cells in different phases of the cell cycle, as well as sub-diploid (sub-G1) cells, were determined.Data were obtained using the DIVA 6.0 software.Three independent experiments were performed in duplicate [104].

Apoptosis Analysis by Flow Cytometry
For the evaluation of the type of cell death induced by CBEO in SK-MEL-28 cells, an Annexin V conjugated with a fluorescein isothiocyanate (FITC) and propidium iodide (PI) staining assay was performed.In this assay, SK-MEL-28 cells were plated (2 × 10 5 cells/mL) in 24-well plates and incubated with CBEO (20 or 40 µg/mL) or DXR (4 µM) for 48 h.After, the cells were collected and centrifuged (500× g, 20 • C, 5 min), they were resuspended in a binding buffer.Subsequently, the cells were labeled with Annexin V-FITC and incubated at room temperature for 10 min in the dark.After this period, the cells were washed with the binding buffer, centrifuged (500× g, 20 • C, 5 min), resuspended in the binding buffer, and labeled with PI (20 µg/mL).Then, the analysis was conducted using a flow cytometer, acquiring 10,000 events/sample.Data were obtained using the DIVA 6.0 software.Three independent experiments were performed in duplicate [104].

Apoptosis Analysis by Laser Confocal Microscopy
In this assay, a morphological analysis of SK-MEL-28 cells after CBEO treatment was conducted.The cells were stained with acridine orange (AO) and propidium iodide (PI).SK-MEL-28 cells were plated (2 × 10 5 cells/mL) in 24-well plates and incubated with CBEO (20 or 40 µg/mL) or DXR (4 µM) for 48 h.Subsequently, the cells were collected and centrifuged (500× g, 20 • C, 5 min), resuspended in PBS, and stained with AO (1 mg/mL) and the PI solution (10 µg/mL).The stained cells were observed under a laser scanning confocal microscope.The established criteria were as follows: (a) viable cells exhibited a clear green nucleus with an intact structure; (b) early apoptotic cells displayed a bright green nucleus, indicating chromatin condensation; (c) late apoptotic cells exhibited dense orange areas (stained in both green and red) of chromatin condensation and the formation of apoptotic blebs on the membrane; (d) necrotic or dead cells were stained only in red [105][106][107].Three independent experiments were performed in duplicate.

Docking Prediction
For the in silico assays, the Molegro Virtual Docker (MVD) v.6.0.1 software was employed (molexus.io/molegro-virtual-docker/).Molecular docking was carried out with enzymes that already had their co-crystallized ligands in their structures: Extracellular Signal-Regulated Kinase 1 (ERK1) (Protein Data Bank (PDB) ID: 5LCJ) and c-Jun N-terminal Kinase 1 (JNK1) (PDB ID: 2G01).Additionally, structures of p38 Mitogen-Activated Protein Kinase α (p38α MAPK) (PDB ID: 1R39), Nuclear Factor kappa B (NF-κB; p50/p65) (PDB ID: 1VKX), and Protein Kinase B (PKB/AKT) (PDB ID: 1GZN) were used.In the absence of their ligands in the PDB, catalytic site coordinates were obtained from the literature [73,[108][109][110].The tested molecules were designed using Marvin Sketch v. 19.18 software.Subsequently, the compounds were standardized using ChemAxon's Standardizer v. 21.2.0 software, which added hydrogen atoms, standardized the aromatic ring, removed salts, converted the structure to 3D, and then subjected the compounds to molecular docking.All water molecules were removed from the crystalline structure.A template was created between the enzyme and the co-crystallized ligand to delineate the active site of the macromolecule.Subsequently, the molecules were added, and the docking prediction was performed.Prior to the docking simulation, redocking was performed, which corresponds to the root mean square deviation (RMSD) calculated from the poses, indicating the degree of confidence in

Figure 2 .
Figure 2. Percentage of cells in the different phases of the cell cycle and at the sub-G1 peak after 48 h of treatment with essential oil from Conyza bonariensis (L.) aerial parts (CBEO) (20 or 40 µg/mL) or doxorubicin (DXR) (4 µM).Data are obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences between experimental groups for each condition; p < 0.05.

Figure 3 .
Figure 3.Effect of the essential oil from Conyza bonariensis (L.) aerial parts (CBEO) and doxorubicin (DXR) on human melanoma SK-MEL-28 cells labeled with annexin V-fluorescein isothiocyanate (FITC) and/or propidium iodide (PI) after 48 h of treatment.(A) Representative dotplots were obtained through a flow cytometry analysis.The dotplots were divided into four quadrants representing different cellular populations: viable cells (annexin V-FITC-/PI−, lower left quadrant), cells in early apoptosis (annexin V-FITC+/PI−, lower right quadrant), cells in late apoptosis/necrosis (annexin V-FITC+/PI+, upper right quadrant), and dead cells (annexin V-FITC−/PI+, upper left quadrant).A total of 10,000 events/samples were acquired using red fluorescence detectors (PI, 325-488 nm) and green fluorescence detectors (FITC, 493-525 nm).(B) Percentage of cells labeled with annexin V-FITC and/or PI.Data were obtained from three independent experiments performed in duplicate and analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters indicate significant differences between experimental groups for each condition.p < 0.05.

Figure 4 .
Figure 4. Apoptosis induction by the essential oil from Conyza bonariensis (L.) aerial parts (CBEO) or doxorubicin (DXR) analyzed through acridine orange (AO) and propidium iodide (PI) staining in human melanoma SK-MEL-28 cells after 48 h of treatment.Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences between experimental groups for each condition; p < 0.05.In addition, representative images of cells treated with CBEO (20 µg/mL) or DXR, which exhibit characteristics of early and late apoptosis, are shown in Figure 5.

Figure 4 . 26 Figure 5 Figure 5 .
Figure 4. Apoptosis induction by the essential oil from Conyza bonariensis (L.) aerial parts (CBEO) or doxorubicin (DXR) analyzed through acridine orange (AO) and propidium iodide (PI) staining in human melanoma SK-MEL-28 cells after 48 h of treatment.Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences between experimental groups for each condition; p < 0.05.In addition, representative images of cells treated with CBEO (20 µg/mL) or DXR, which exhibit characteristics of early and late apoptosis, are shown in Figure 5.Pharmaceuticals 2023, 16, x FOR PEER REVIEW 6 of 26

Figure 5 .
Figure 5. Representative images of SK-MEL-28 cells after 48 h of treatment with essential oil from Conyza bonariensis (L.) aerial parts (CBEO) (20 µg/mL) or DXR (4 µM), stained with acridine orange (AO) and/or propidium iodide (PI).(A) Viable cells exhibit a light green nucleus and an intact structure.Early apoptotic cells display a light green nucleus showing chromatin condensation.Late apoptotic cells show dense orange areas (green/red) of chromatin condensation and membrane blebs, while necrotic cells have a red nucleus.(B) Representative enlarged images of SK-MEL-28 cells after 48 h of incubation with CBEO (20 µg/mL) or DXR (4 µM).The yellow arrow indicates chromatin condensation, the blue arrow indicates membrane blebs, and the red arrow indicates DNA fragmentation.

Figure 7 .
Figure 7. Extracellular Signal-Regulated Kinase 1 and 2 (ERK1/2) modulation in SK-MEL-28 after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/m DXR (4 µM).(A) The cell population was determined from the analysis of cell size dotplots ( Forward Scatter, X-axis) and cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by cytometry.For the analysis of ERK1/2 fluorescence histograms (PerCP-Cy5.5), a region corresp ing to the autofluorescence of cells not labeled with the antibody was initially delimited.Onl percentage of fluorescent cells after this demarcated region was considered positive.(B) Overl representative histograms showing ERK1/2 fluorescence (PerCP-Cy5.5,X-axis) and the numb cells/events (Y-axis) in different experimental groups.(C) Graphical representation of the re obtained by flow cytometry.Data were obtained from three independent experiments carrie in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analy Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differe among conditions; p < 0.05.

Figure 7 .
Figure 7. Extracellular Signal-Regulated Kinase 1 and 2 (ERK1/2) modulation in SK-MEL-28 cells after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) The cell population was determined from the analysis of cell size dotplots (FSC-Forward Scatter, X-axis) and cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by flow cytometry.For the analysis of ERK1/2 fluorescence histograms (PerCP-Cy5.5), a region corresponding to the autofluorescence of cells not labeled with the antibody was initially delimited.Only the percentage of fluorescent cells after this demarcated region was considered positive.(B) Overlay of representative histograms showing ERK1/2 fluorescence (PerCP-Cy5.5,X-axis) and the number of cells/events (Y-axis) in different experimental groups.(C) Graphical representation of the results obtained by flow cytometry.Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences among conditions; p < 0.05.

Figure 8 .
Figure 8. c-Jun N-terminal Kinase (JNK) modulation in SK-MEL-28 cells after treatment with tial oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) T population was determined from the analysis of cell size dotplots (FSC-Forward Scatter, X-ax cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by flow cytometry.For the ana JNK fluorescence histograms (APC), a region corresponding to the autofluorescence of cells beled with the antibody was initially delimited.Only the percentage of fluorescent cells af demarcated region was considered positive.(B) Overlay of representative histograms showi fluorescence (APC, X-axis) and the number of cells/events (Y-axis) in different experimental g (C) Graphical representation of the results obtained by flow cytometry.Data were obtaine three independent experiments carried out in duplicate and expressed as mean ± standard e the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's te ferent letters denote significant differences among conditions; p < 0.05.

Figure 8 . 26 Figure 9 .
Figure 8. c-Jun N-terminal Kinase (JNK) modulation in SK-MEL-28 cells after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) The cell population was determined from the analysis of cell size dotplots (FSC-Forward Scatter, X-axis) and cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by flow cytometry.For the analysis of JNK fluorescence histograms (APC), a region corresponding to the autofluorescence of cells not labeled with the antibody was initially delimited.Only the percentage of fluorescent cells after this demarcated region was considered positive.(B) Overlay of representative histograms showing JNK fluorescence (APC, X-axis) and the number of cells/events (Y-axis) in different experimental groups.(C) Graphical representation of the results obtained by flow cytometry.Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences among conditions; p < 0.05.

Figure 9 .
Figure 9. p38 Mitogen-Activated Protein Kinase (p38 MAPK) modulation in SK-MEL-28 cells after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) The cell population was determined from the analysis of cell size dotplots (FSC-Forward Scatter, X-axis) and cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by flow cytometry.For the analysis of p38 MAPK fluorescence histograms (PE-Cy7), a region corresponding to the autofluorescence of cells not labeled with the antibody was initially delimited.Only the percentage of fluorescent cells after this demarcated region was considered positive.(B) Overlay of representative histograms showing p38 MAPK fluorescence (PE-Cy7, X-axis) and the number of cells/events (Y-axis) in different experimental groups.(C) Graphical representation of the results obtained by flow cytometry.Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences among conditions; p < 0.05.

Pharmaceuticals 2023 , 26 Figure 11 .Figure 11 .
Figure 11.Nuclear Factor kappa B (NF-κB) modulation in SK-MEL-28 cells after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) The cell population was determined from the analysis of cell size dotplots (FSC-Forward Scatter, X-axis) and cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by flow cytometry.For the analysis of NF-κB fluorescence histograms (PE), a region corresponding to the autofluorescence of cells not labeled with the antibody was initially delimited.Only the percentage of fluorescent cells after this demarcated region was considered positive.(B) Overlay of representative histograms showing NF-κB fluorescence (PE, X-axis) and the number of cells/events (Y-axis) in different experimentalFigure 11.Nuclear Factor kappa B (NF-κB) modulation in SK-MEL-28 cells after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) The cell

26 Figure 12 .
Figure 12.Protein Kinase B (PKB/AKT) modulation in SK-MEL-28 cells after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) The cell population was determined from the analysis of cell size dotplots (FSC-Forward Scatter, X-axis) and cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by flow cytometry.For the analysis of PKB/AKT fluorescence histograms (PE-Cy7), a region corresponding to the autofluorescence of cells not labeled with the antibody was initially delimited.Only the percentage of fluorescent cells after this demarcated region was considered positive.(B) Overlay of representative histograms showing PKB/AKT fluorescence (PE-Cy7, X-axis) and the number of cells/events (Y-axis) in different experimental groups.(C) Graphical representation of the results obtained by flow cytometry.Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences among conditions; p < 0.05.

Figure 12 .
Figure 12.Protein Kinase B (PKB/AKT) modulation in SK-MEL-28 cells after treatment with essential oil from Conyza bonariensis aerial parts (CBEO) (20 or 40 µg/mL) or DXR (4 µM).(A) The cell population was determined from the analysis of cell size dotplots (FSC-Forward Scatter, X-axis) and cytoplasmic granularity (SSC-Side Scatter, Y-axis) obtained by flow cytometry.For the analysis of PKB/AKT fluorescence histograms (PE-Cy7), a region corresponding to the autofluorescence of cells not labeled with the antibody was initially delimited.Only the percentage of fluorescent cells after this demarcated region was considered positive.(B) Overlay of representative histograms showing PKB/AKT fluorescence (PE-Cy7, X-axis) and the number of cells/events (Y-axis) in different experimental groups.(C) Graphical representation of the results obtained by flow cytometry.Data were obtained from three independent experiments carried out in duplicate and expressed as mean ± standard error of the mean (SEM) analyzed by Analysis of Variance (ANOVA) one-way followed by Tukey's test.Different letters denote significant differences among conditions; p < 0.05.