1. Introduction
Plants have been an inestimable source of anticancer compounds, providing approximately 50% of the approved chemotherapeutics [
1,
2,
3]. In fact, natural compounds (NPs) show promising features, such as scaffold diversity and structural complexity, that make them ideal for drug discovery [
4,
5].
Arbutus unedo L. (
A. unedo) is a Mediterranean plant used in traditional medicine to treat various illnesses, such as gastrointestinal, dermatological, cardio-vascular, and urological disorders, kidney diseases, or also used as a diuretic and antidiabetic, suggesting the presence of numerous NPs with promising pharmacological activities [
6]. Regarding potential uses in cancer, studies have reported that NPs extracted from the entire plant, fruits, honey, and leaves could promote cytotoxicity in cellular tumor models. In particular, it has been recently reported that the leaf extract of
A. unedo reduced the viability of U2OS, a cellular model of osteosarcoma, without toxicity in human umbilical vein endothelial (HUVEC) cells, a healthy cellular model [
7]. However, further investigations are required to identify the NPs responsible for the observed anticancer activity.
GISTs are rare mesenchymal neoplasms caused by a gain of function mutation in KIT or PDGFRα, two tyrosine kinase receptors (TKRs). For this reason, they are treated with imatinib, a tyrosine kinase inhibitor (TKI) [
8]. Although imatinib has significantly improved patient prognosis, providing remarkable amelioration of their life expectations, the therapy is not conclusive. Indeed, imatinib mainly promotes partial response or stable disease in most patients, while complete response is only observed in 5% of patients [
9]. Moreover, imatinib treatment is frequently associated with drug resistance, mostly due to acquired secondary mutations in TKRs (50–60% of patients) [
10,
11], fostering the development and approvals of novel multi-target TKIs as further treatment lines [
12]. Unfortunately, resistance has not been effectively treated with a strategy exclusively based on TKIs, which have provided a modest increase in progressive free survival (PFS). Preclinical studies have reported that imatinib resistance could be more multifaceted than initially hypothesized, supporting the research of novel compounds that could broadly target imatinib-resistant GIST cells [
13]. Traditional chemotherapy could potentially represent an alternative strategy for treating GISTs. Indeed, the lack of cell specificity, a limitation of traditional chemotherapy, could be an advantage, potentially targeting both TKI-sensitive and -resistant subclones through a non-specific mechanism of action. However, the available chemotherapeutics, including those effective in other soft tissue sarcomas, such as doxorubicin, had already been tested in GISTs before imatinib approval, showing no effect [
14,
15]. Thus, no standard chemotherapy has been approved in GIST so far [
16].
Therefore, considering that identifying novel chemotherapeutics could be promising for treating GISTs and the interesting anticancer properties of A. unedo already reported in the literature, we hypothesized that A. unedo could harbor NPs with chemotherapeutic properties. Therefore, we observed that an extract obtained from A. unedo leaves affected the viability of GIST cell lines characterized by different mutations and degrees of imatinib response. Furthermore, a bio-guided assay fractionation procedure was employed to investigate the different compounds’ contribution to the crude extract’s activity.
3. Discussion and Conclusions
Arbutus unedo L. (Ericaceae family), known as the strawberry tree, is an evergreen shrub growing in circum-Mediterranean regions. Leaves are used in the traditional medicine of the Iberian Peninsula and Sardinia (Italy) to treat many illnesses [
6], demonstrating a number of health-promoting properties associated with this plant. In this work, we showed that an extract from
Arbutus unedo L. leaves (indicated as AUN) induces early apoptosis in GIST cells, exerting chemotherapeutic properties. This result is of great importance, considering chemotherapy is generally ineffective in advanced GISTs, and no options are available.
Patients’ treatment first relies on imatinib as the first-line therapy. Indeed, primary GISTs are commonly associated with gain-of-function mutations in KIT or PDGFRα (about 90% of cases), the disease drivers which, once mutated, show a ligand-independent dysregulated activity [
19,
20]. Moreover, given that secondary mutations in both KIT or PDGFRα have been found in imatinib-resistant GISTs, multi-target TKIs, such as sunitinib, regorafenib, and ripretinib, are administered as further treatment lines after the onset of resistance [
12]. However, multi-target TKIs have not kept their promises, only modestly improving the outcomes of imatinib-resistant GIST patients. Indeed, imatinib-resistant subclones escape through alternative pathways that are not targetable from multi-target TKIs, suggesting that a TKI-based approach alone may not represent the only or best option for GIST patients [
13]. Therefore, GIST patients urgently require the identification of novel therapies other than TKIs. In this context, we hypothesized that plant extracts could represent an important source of novel NPs to be studied in GIST. Among these, AUN showed antiproliferative activity in osteosarcoma (U2OS) cells [
7], representing a promising unexplored therapeutic strategy in GIST. In view of this consideration, we tested the anticancer activity of AUN in GIST cellular models. AUN significantly impaired GIST cell viability, triggering a significant Annexin (+)/7-AAD (−) cell population, a marker of early apoptosis. AUN was further explored by bio-guided fractionation to deepen what could be the most active part of AUN. Hence, AUN was first fractionated by a liquid/liquid partition, separating roughly AUN phytochemicals on the basis of their polarity. This first step led to the identification of the FR2 fraction as the most active, able to impair cell viability and promote early apoptosis in different models of imatinib-sensitive GISTs (GIST-882 and GIST-T1), i.e., characterized by different primary kit mutations. In view of this consideration, we performed a second step of fractionation, subjecting the FR2 fraction to a reverse-phase MPLC, which led to the identification of the FR2-A fraction as the most active one in imatinib-sensitive cellular models (GIST-882 and GIST-T1), regardless of primary KIT mutations. Moreover, it resulted in being even more potent than FR2, as it significantly impaired cell viability at lower concentrations (12.5 μg/mL vs. 50 μg/mL). Afterward, FR2-A was tested in imatinib-resistant cell lines (GIST48 and GIST-48b), once again proving effective in impairing cell viability, providing new insight into the FR2-A mechanism of action. Indeed, the cell line GIST-48b is characterized by the absence of KIT expression, supporting the hypothesis that the mechanism is independent from KIT, possibly paving the way for developing therapeutic strategies other than TKIs in imatinib resistance. Additionally, FR2-A was more efficient with respect to imatinib in sensitive cells, inducing a higher percentage of Annexin V (+)/7-AAD (−) cells. This result confirmed the capacity of AUN to efficiently induce early apoptosis, prompting us to deepen the mechanism. Nevertheless, we did not observe the cleavage of caspase 3, caspase 9, and PARP-1, suggesting a cell death mechanism other than classical apoptosis. Interestingly, although not cleaved, a remarkable and rapid PARP-1 downregulation after FR2-A treatment was observed, possibly suggesting an epigenetic mechanism. Aware of the limitation, mainly due to the limited number of GIST cell models analyzed—though it depends on the availability of commercial cell lines—we are conscious that our results are speculative. However, recent findings in the literature seem to support our hypothesis of an epigenetic mechanism of action. Indeed, a study has identified PARP-1 as a target of miR-7-5p in a model of lung cancer cells, demonstrating that sponging miR-7-5p promotes the homologous repair (HR) path through upregulating PARP1 expression [
21]. In another study, Garmutin-A (GA), an NP isolated from
Garcinia multiflora, induced apoptosis via the upregulation of miR-17-5p, triggering PARP-1 downregulation, in leukemic CB3 cells [
22]. Overall, this finding is in line with the evidence that many NPs could exert their anticancer activity by affecting miRNAs expression [
23]. Therefore, a possible speculation is that FR2-A could upregulate different miRNAs, thereby altering the PARP-1-mediated repair capability. Another thought offered by PARP-1 downregulation is the possibility of re-thinking imatinib treatment in GIST, considering synergistic/combinatory approaches. Indeed, most GIST patients commonly achieve partial response or stable disease, suggesting that GIST cells could survive the treatment by activating adaptive responses. In this context, our study corroborates that imatinib mainly stabilizes GIST cells in a non-proliferative state rather than inducing cell death. Therefore, a combination of imatinib and PARP-1 inhibitors could represent a promising strategy that has not yet been investigated in GIST.
Regarding the possible bio-active compound characterizing FR2-A, we initially focused on β-arbutin. Indeed, this phytochemical, already identified as one of the primary metabolites in AUN [
24], was identified as the most representative compound in FR2-A. Strengthened by this result, together with the knowledge that β-arbutin was already associated with proapoptotic activity in melanoma cells [
25], we tested pure β-arbutin, unfortunately with negative results. This prompted us to perform an additional fractionation step, this time with a size-exclusion approach, to increase the concentrations of active phytochemicals. This approach led to sixty-one different subfractions, of which nine, forming four different clusters, exhibited a significant biological activity in GIST cells. Interestingly, all clusters contained phytochemicals characterized by the pyrogallol moiety, a polyphenolic structure already associated with proapoptotic activity. A study conducted on As4.1 juxtaglomerular cells, a model of benign kidney tumors, demonstrated that pyrogallol treatment could efficiently induce apoptosis [
26]. Interestingly, the authors observed caspase 3 and PARP-1 cleavage; however, treatment with caspase 3 inhibitors did not prevent apoptosis, representing an intriguing result. Indeed, the authors did not evaluate the PARP-1 status in the presence of a caspase 3 inhibitor. Therefore, other mechanisms besides caspase activation could not be excluded. Another study demonstrated that the pyrogallol moiety in polyphenols had an important role to play in apoptosis induction in human embryonic kidney cells (HEK293T) and the chronic myelogenous leukemia cell line (K562), even though the mechanism was not investigated [
27]. In the same study, the authors highlighted that the pyrogallol moiety was important for cytotoxic activity. In agreement with this, all the FR2-A active subfractions showed cytotoxicity in PBMCs, similar to those observed with doxorubicin, a traditional chemotherapeutic. These findings suggest a non-targeted but chemotherapy-like mechanism of action, which should not be underestimated considering the lack of effective chemotherapeutic agents for GIST patients. Indeed, doxorubicin, approved in several solid tumor treatments, has been tested in clinical trials with poor response in GISTs, results that led to its non-approval [
14]. In confirmation of this, unlike FR2-A clusters, doxorubicin does not impair the viability of imatinib-resistant cells [
28]. Therefore, deepening FR2-A composition could be of great interest in GIST, leading to the identification of novel chemotherapeutic agents to be used as alternatives or in combination with the approved TKI. In view of these considerations, the most active FR2-A subfractions (2A-35 and 2A-36) were mainly characterized by the presence of gallic acid, a phenolic acid which has been identified as a promising health-promoting agent in many conditions, including cancer [
29]. Despite the expectations, our finding ascribed an activity to gallic acid, which, however, was weaker than 2A-35. Therefore, gallic acid cannot be considered a solely bioactive compound but rather one that acts in concert with other components, highlighting the importance of the phytocomplex.
In conclusion, although bio-guided fractionation did not lead to the identification of a single compound active in GIST cellular models, it was confirmed as a feasible tool for circumscribing phytochemicals of interest. In particular, in the present study, we identified pyrogallol-bearing compounds as the main players in GIST cell cytotoxicity, harboring chemotherapeutic properties in imatinib-sensitive and resistant cells. These data may fuel further studies searching for novel treatments in GIST, including additional strategies that augment the existing treatments, resulting in better patient outcomes. Therefore, the chemotherapeutic-like activity should not discourage further preclinical studies on pyrogallol-bearing compounds in GISTs to conclude whether they could represent an alternative strategy for treating GISTs.
4. Material and Methods
4.1. GIST Cell Lines and Cell Culture Conditions
GIST-882 and GIST-T1 are defined as primary mutated and imatinib-sensitive cellular models. They harbor exon 9 (K642E) homozygous mutations and exon 13 heterozygous (V560-Y579) in the KIT gene, respectively. GIST-882 and GIST-T1 were grown in RPMI-1640 supplemented with 15% FBS. GIST-48 is instead reported as an imatinib- and sunitinib-resistant cell line harboring a primary homozygous mutation on KIT exon 11 (V560D) and an additional secondary heterozygous mutation in exon 17 (D820A). GIST-48b was established in vitro, starting from GIST-48 after HSP-90 inhibitor (17-AAG) drug pressure selection, resulting in a subline characterized by nearly undetectable KIT transcript and protein. GIST-48 and GIST-48b were grown in Iscove’s Modified Dulbecco’s Medium (IMDM) supplemented with 15% FBS. All the indicated cell lines were routinely tested to avoid mycoplasma contamination with MycoBlue Mycoplasma Detector (Vazyme, Nanjing, China). GIST-T1 was purchased from CosmoBio (Tokyo, Japan), while GIST-48 and GIST-882 cells were kindly provided by Fletcher JA, MD (Harvard Medical School).
4.2. Plant Material
The leaves of
Arbutus unedo L. were harvested in Sardinia in 2018 (voucher specimen CAG 878/v, deposited at the General Herbarium of the Department of Life and Environmental Sciences, University of Cagliari). To ensure reproducibility, the plant material was analyzed by
1H NMR profiling [
25].
4.3. Cell Viability and Apoptotic Profile
105 GIST cells were seeded in a 24-well cell culture plate the day before treatment. Cells were treated with AUN or AUN subfractions at the indicated final concentrations in the culture medium for the indicated time points. The analysis of cell viability and the presence of apoptotic cells were evaluated on the entire cellular population, including cells that were in suspension due to the treatment. For this reason, the culture medium and the phosphate-buffered saline (PBS) used to rinse cells before trypsinization were also collected. Adherent cells were harvested by trypsinization and combined with cells in the culture medium and PBS. The entire sample was then centrifugated and resuspended in 500 µL of fresh medium. Samples were stained with Guava® ViaCount™ or Guava® Nexin Reagent according to the manufacturer’s instructions. The Guava® ViaCount™ distinguishes viable and non-viable cells based on the differential permeabilities of two DNA-binding dyes. The nuclear dye (in the graphs indicated as “Nucleated Cells” on the y-axis) differentiates nucleated cells from cellular debris, while the viability dye brightly stains dying cells (in the graphs indicated as “Viability” on the x-axis). The number of viable cells in the sample was defined starting from the cell concentration (cells/μL) calculated by the flow cytometer instrument. The monitoring of apoptosis activation was performed using Guava® Nexin Reagent. The reagent combines fluorescently labeled Annexin V (Annexin V-PE) and 7-Aminoactinomycin D (7-AAD). Annexin V-PE, a calcium-dependent phospholipid-binding protein with a high affinity for phosphatidylserine (PS), is a membrane component early exposed on the external cell surface during apoptotic pathway stimulation (indicated in the graphs as “Annexin-V” on the x-axis). Instead, late apoptotic/dead cells are recognized by 7-AAD. This fluorescent DNA intercalator is membrane-impermeant and excluded from live, healthy, and early apoptotic cells (indicated in the graphs as “Viability” on the y-axis).
4.4. AUN Extraction and Fractionation Steps of A. unedo
Five hundred grams of dried and grounded plant material was extracted using 2 L of CH3OH 80% v/v. The extract was filtered and dried in a rotary evaporator; the extraction procedure was repeated four times (24 h each) on the same plant material, obtaining 67.7 g of extract. The extract (AUN) was suspended in 700 mL of water and subjected to liquid/liquid partition using in-series chloroform (CHCl3) or ethyl acetate (EtOAc). Phytochemical extractions were performed using each solvent three times. Three AUN-derived fractions were obtained: CHCl3 (FR1; yield: 0.33% w/w), EtOAc (FR2; yield: 8.8% w/w), and H2O (FR3; yield 70%). The stock solutions for the bioassays of FR3 were prepared by solubilizing dried extracts in distilled H2O at an initial concentration of 20 mg/mL, while FR2 and FR3 were instead solubilized in DMSO 10% at the same initial concentration. FR1, FR2, and FR3 fractions were then diluted (1:1000) during the biological assays, leading both H2O and DMSO to a final concentration of 0.1% in the medium supplemented with FBS 15%. Being the bioactivity localized in the FR2, 4 g of the fraction was suspended in 5 mL of water and injected in a medium-pressure liquid chromatography (MPLC) instrument (Reveleris®, Büchi, Switzerland) using a reverse-phase stationary phase (40 g of C18 column). A gradient of water (solvent A) and methanol (solvent B) was used as an eluent. The gradient was composed of an isocratic phase of 10 min (90% A and 10% B), a gradient from 90% A to 80% A in 1.1 min, an isocratic phase of 20 min (80% A and 20% B), a gradient from 80% A to 70% A in 1.1 min, an isocratic phase of 10 min (70% A and 30% B), a gradient from 70% A to 50% A in 1.1 min, an isocratic phase of 10 min (50% A and 50% B) a gradient from 50% A to 30% A in 1.1 min, an isocratic phase of 5 min (30% A and 70% B), a gradient from 30% A to 0% A in 1.1 min, an isocratic phase of 5 min (0% A and 100% B). The flow rate was 20 mL/min, and the run length was 70 min. Based on the chromatogram and UV–Vis absorbance signals at three wavelengths (UV1 = 254 nm, UV2 = 270 nm, and UV3 = 340 nm), FR2-derived fractions were collected and dried, obtaining a total of four FR2 subfractions (FR2-A, FR2-B, FR2-C, and FR2-D). An analogous preparation of FR2 was followed for its subfractions in the biological assays. Being the bioactivity localized in the FR2-A, 879.8 mg of the fraction was suspended in the minimum amount of methanol and then subjected to size exclusion chromatography using a chromatography column (1800 mm × 25 mm) filled with 220 g of Sephadex (LH-20) and, as an eluent, methanol. The flow rate was 0.4 mL/min. The eluate in each tube was concentrated in a rotary evaporator, and a small quantity was used to perform thin-layer chromatography (TLC). The stationary phase of the TLC used a silica gel matrix with a 254 nm florescent indicator (Sigma-Aldrich), while EtOH:MeOH:H2O (10:1.35:1) was used as the mobile phase. This latter was employed to acquire a first overview of the chemical composition of the fractions by UV–Vis light exposure, enabling us to group the most similar ones, obtaining 84 final fractions. Excluding the fraction yielding less than 1 mg of material, 61 were tested in the bioassays. The most active ones were phytochemically investigated compared to the inactive or weakly active ones, which were eluted immediately before or after.
4.5. NMR Spectra Measurement
1H NMR spectra were recorded at 25 °C on a Varian Inova instrument (equipped with a reverse triple resonance probe) operating at 600.13 MHz. Each 1H NMR spectrum consisted of 256 scans (corresponding to 16 min) with a relaxation delay (RD) of 2 s, acquisition time of 0.707 s, and spectral width of 9595.8 Hz (corresponding to δ 16.0). All the 1H NMR spectra were uploaded onto Zenodo.
4.6. UHPLC-MS Analysis
The UHPLC-MS analysis was run on a Waters ACQUITY ARC UHPLC/MS system consisting of a QDa mass spectrometer equipped with an electrospray ionization interface and a 2489 UV–Vis detector. The detected wavelengths (λ) were 254 nm and 365 nm. The analyses were performed on an XBridge BEH C18 column (10 × 2.1 mm i.d., particle size 2.5 µm) with an XBridge BEH C18 VanGuard Cartridge precolumn (5 mm × 2.1 mm i.d., particle size 1.8 µm). The mobile phases were H2O (0.1% formic acid) (A) and MeCN (0.1% formic acid) (B). Electrospray ionization in positive and negative modes was applied in the 50–1200 Da mass scan range. The FR2-A subfractions were diluted to 100 µg/mL, and a volume of 3 µL was injected. The samples were eluted with the following method: 20% B for one minute, followed by a gradient reaching 95% B in 3 min; 95% B was kept for 1 min; then the gradient reached 20% B in 0.2 min again; and 20% B was kept for 2 min. The flow rate was 0.8 mL/min.
4.7. Software and Statistical Analysis
Statistical analysis was performed using GraphPad Prism software 8.0.2, applying the appropriate statistical test. Details are indicated below each figure. IC50 was instead calculated with the free “aatbioquest” tool (
https://www.aatbio.com/tools/ic50-calculator, accessed on 11 April 2024).
4.8. Western Blot
Whole-cell protein lysates were prepared using NP40 buffer containing protease inhibitors (halt protease and phosphatase inhibitor cocktail; Thermo Fisher Scientific) and 1 mM of phenylmethylsulfonyl fluoride (Sigma-Aldrich, St. Louis, MI, USA). Proteins were separated in SDS-PAGE (12%) and transferred onto nitrocellulose membranes. Membranes were blocked by 5% skimmed milk, followed by incubation at 4 °C overnight with primary antibodies. Primary antibodies against KIT (A4502; Dako, Glostrup, Denmark), phospo-KIT (3391; Cell Signaling, Danvers, MA, USA), PARP-1 (9542, Cell-Signaling), and actin (A1978; Sigma-Aldrich) were used. After rinsing, membranes were incubated with horseradish peroxidase-conjugated secondary antibody (Thermo Fisher Scientific) at room temperature for 2 h. After further rinsing, immunoreactive bands were visualized by enhanced chemiluminescence (BioRad, Hercules, CA, USA), and signals were captured and quantified using ChemiDoc (BioRad).
4.9. Chemical Compounds
Sigma-Aldrich supplied β-arbutin and gallic acid for the biological assays. Selleck Chemicals supplied imatinib.
4.10. Cell Viability by MTT
2.5 × 104 cells were seeded in a 96-well culture plate the day before the treatment. Cells were treated with FR2-A-derived fractions, and viability was analyzed 24 h after treatment. Treatment was removed, and cells were incubated with MTT reagent (0.5 mg/mL) in a medium without serum for 2 h. At the end of the incubation, the MTT solution was removed carefully, and formazan crystals were dissolved in DMSO. Absorbance was read at 492 nm using a TECAN spectrophotometer.
4.11. PBMC Isolation
Peripheral blood lymphocyte (PBMC) isolation was performed as described [
30]. Briefly, PBMCs were isolated from blood samples of healthy donors (Buffy coat) using a density gradient centrifugation with Histopaque-1077. Donors were healthy, non-smokers, and with no known exposure to genotoxic chemicals or radiation. Authorization for the use of human blood samples for research purposes was received from AUSL Bologna IT, S. Orsola-Malpighi Hospital.