1. Introduction
Glioblastoma multiforme (GBM) is one of the most aggressive malignancies and the most common malignant primary tumor of the central nervous system. GBM patients have a median survival of less than 15 months from initial diagnosis. GBM originates from astrocytic glial cells and it is a high-grade malignant glioma (grade IV) [
1]. Standard therapy for GBM consists of surgical resection followed by radiotherapy and chemotherapy with temozolomide (TMZ). Despite maximal surgical resection and adjuvant therapy, almost all GBM tumors locally recur after treatment. Current challenges for the treatment of GBM include its incomplete resection, a high degree of genetic heterogeneity, the blood-brain barrier, and an immunosuppressive microenvironment. The highly infiltrative nature of GBM makes complete resection at the cellular level almost impossible. Although tumor margins can be roughly determined by imaging, GBM grows with microscopic finger-like projections, imperceptible to presurgical or even intraoperative imaging techniques [
2]
Fluorescence-guided surgery allows for better visualization of tumor tissue and prolongs patient survival. 5-aminolevulinic acid (ALA) is a precursor metabolized in the heme biosynthesis pathway to protoporphyrin IX (PpIX), an endogenous photosensitizer. Finally, the enzyme ferrochelatase incorporates iron into PpIX to form heme. PpIX accumulates in cancer cells after ALA administration due to abnormalities in the heme pathway, such as the poor ferrochelatase activity. ALA administration results in the PpIX fluorescence of high-grade gliomas [
3].
Photodynamic therapy (PDT) is based on light-excited photosensitizers to generate lethal reactive oxygen species (ROS), such as singlet oxygen (
1O
2), superoxide radicals (·O2-) and hydroxyl radicals (·OH). The resultant ROS can cause irreversible damage to cancer cells and lead to cell death, immune response, and microvascular damage. ALA and its derivatives (such as methyl aminolevulinate and hexyl aminolevulinate) are not photosensitive, but they are converted endogenously into PpIX photosensitizer by the heme biosynthesis pathway. PpIX tends to accumulate in malignant cells and can be activated by light to generate ROS, causing cancer cell damage [
4]. In recent years, PDT studies have increased for several cancers, including GBM. The United States Food and Drug Administration (FDA) has approved ALA for fluorescence-guided resection of high-grade gliomas. Consequently, the interest in using this PpIX precursor to administer PDT in these tumors has been renewed [
5].
The first well-established molecular alterations in GBM were relative to isocitrate dehydrogenase IDH1/2 mutations, which became a factor for categorization and prognostication of these tumors. Other alterations in primary GBMs were recognized comprising the proto-oncogenes EGFR, PDGFR, MET, CDK4, CDK6, CCND1, CCDN3, MDM2, MDM4, and MYCC and the tumor suppressor genes TP53, RB1, CDKN2A (p16INKA4/p14ARF), and PTEN [
6]. Epidermal growth factor receptor (EGFR) is a transmembrane tyrosine kinase encoded by a gene located on chromosome band 7p12. EGFR is a membrane of the ErbB receptors family. The amplification and over-expression of EGFR have shown to be associated with tumor cell proliferation, angiogenesis, tumor invasion and resistance to therapeutic approaches such as chemotherapy and radiotherapy [
7]. EGFR has a kinase domain composed of the N and C lobes and an adenosine triphosphate (ATP) binding cleft between the two lobes. When a ligand binds to its EGFR receptor, an asymmetric dimer is formed and the phosphate of ATP is transferred to the tyrosine residue of the regulatory domain. Various proteins are bound to this phosphorylated tyrosine and signals are transmitted downstream through the rat sarcoma (RAS)-rapidly accelerated fibrosarcoma (RAF)-mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3-kinase (PI3K)-protein kinase B (PKB/AKT) pathways. EGFR tyrosine kinase inhibitors (TKIs) competitively inhibit ATP with cleft of the kinase domain [
8]. For GBM, clinical trials have been developed inhibiting EGFR using drugs, peptide vaccines and antibodies [
9].
Populations of GBM cells were obtained from the T98G line resistant to PDT with a methylated derivative of ALA (Me-ALA) for characterization compared with the original parental population. We had previously determined by RT-qPCR that the resistant cells have increased EGFR mRNA levels compared to the parental populations [
10]. The objective of this work was to evaluate the general response of T98G GBM cells resistant to PDT when EGFR activity is inhibited with the drug erlotinib.
2. Materials and Methods
2.1. Bioinformatic Analysis of EGFR Expression in Human Samples
For a bioinformatic analysis of EGFR expression in human samples GEPIA 2 (Gene Expression Profiling Interactive Analysis 2) platform was used (
http://gepia2.cancer-pku.cn/). GEPIA uses the databases The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) [
11]. EGFR gene expression in human tumor samples and normal tissues and the GBM patients’ survival based on the expression status of EGFR gene were analyzed.
2.2. Cell Culture
T98 G (ATCC) cells were cultured in high glucose Dulbecco’s Modified Eagle medium (DMEM; Gibco), supplemented with 10 % v/v fetal bovine serum (Internegocios S.A.) and 1% v/v antibiotic-antimicotic (penicillin 10,000 units/mL, streptomycin 10,000 μg/mL, amphotericin B 25 μg/mL; Gibco). Cells were incubated at 37°C in an atmosphere containing 5% CO2.
2.3. Photosensitizer
The compound used in this study was the methyl-5-aminolevulinic acid (Me-ALA; Sigma-Aldrich), as a precursor of the photosensitizer PpIX. A stock solution of 100 mM Me-ALA was made in sterile phosphate buffered saline (PBS), from which 1 mM work solution was made employing DMEM without serum.
2.4. Light Source
For PDT treatment, cells were irradiated employing a monochromatic light source (635 nm ± 17 nm) with a multi-LED system (coherent light) at an irradiation intensity of 16.9 mW/cm2 (as measured by Coherent Lasermate power meter).
2.5. Obtainment of PDT-Me-ALA Resistant Cells
T98 G resistant to PDT/Me-ALA cells were obtained as described in Vilchez
et al., 2021 [
10]. The irradiation dose that caused cellular death rate of 70–90 % in parental cells (8.60 J/cm
2) was employed with the aim of obtaining resistant cells. T98G cells were cultured in 35 mm plates and incubated with Me-ALA 1 mM for 4 h. Thereafter, cells were exposed to red light. The surviving cells were harvested and replated at 24-48 h after PDT. After proliferation, cells were submitted to a new PDT treatment. As the number of rounds of PDT increased, less cell death was observed. The final population received a total of 8 cycles of PDT, at which time no death was observed. The initial population, not subjected to PDT, was called parental population; the cellular population submitted to one PDT treatment was called first resistant generation and so the following generations were named consecutively. Once resistant cells were obtained, they were kept in frozen stocks. To check the resistance abilities of 8th resistant generation respect to parental population, cells were defrosted and cell viability assay (3-[4,5-dimethylthiazol-2-yl]2, 5-diphenyltetrazolium bromide; MTT) was made after PDT treatment with several irradiation doses as indicated in the next section.
2.6. Viability Assay by MTT
2.6.1. Resistance Determination
The grade of resistance to PDT was measured in TG98 8th generation cells (resistant cells), with respect to parental cells. 12 × 104 cells/ml of each population were seeded in multiwell-96 plates (P96). After 24 h, cells were incubated with Me-ALA 1 mM for 4 h at 37°C and then irradiated at several light doses: 5.07, 6.08 and 8.11 J/cm2 (n = 8). Cells with drug, but without light were employed as drug controls and cells with light (at 8.60 J/cm2), but without drug were employed as light controls. Cells without drug and without light were employed as controls and they were expressed as 100 % of viability. After treatments, culture medium was replaced by complete medium. After 24 h of PDT, the cell viability was analyzed by MTT assay. 10 μL of MTT (Sigma-Aldrich) solution (5 mg/ml PBS) was added to the cells in each well and it was incubated for 3− 4 h. The culture medium was removed and formazan crystals were suspended in DMSO (Cicarelli). Measures of the absorbance were taken at 540 nm with a spectrophotometer (Thermo Scientific, Multiskan FC). Results are reported as the mean ± standard error of mean (SEM). The experiment was made 3 times. Results of a representative experiment are shown.
2.6.2. Cell Viability After Combined Treatment of PDT and Erlotinib
Parental and resistant T98G cells were seeded on P96 plates (12 × 10
4 cells/ml). The next day the cells were incubated with erlotinib at doses of 9.20, 25 and 50 µM (n=8) [
12]. Cells without drugs were used as controls. DMSO vehicle control groups were included. At 24 h the cell viability was analyzed by MTT assay, as described in the previous section. Also, cell viability was determined after combined treatment of PDT with erlotinib (50 µM). Cells were incubated with Me-ALA (1 mM) and erlotinib (50 µM) for 4 hours and then irradiated at 6.08 and 8.11 J/cm
2. After irradiation, the medium was replaced by a complete medium. Erlotinib was added in erlotinib and erlotinib+PDT groups. After 24 h, cell viability was measured by MTT. Results are reported as the mean ± standard error of mean (SEM). The experiment was made 3 times. Results of a representative experiment are shown.
2.7. Flow Cytometry
The fluorescence of the photosensitizer PpIX was measured by flow cytometry in parental and resistant T98G populations after incubation with Me-ALA and erlotinib. Cells (20 × 104 cells/ml) were seeded in Petri dishes (35 mm plates) (n = 6). After 24 h, cells were incubated with Me-ALA (1 mM) and erlotinib (50 µM) for 4 h, washed with PBS three times, trypsinized, resuspended in complete medium (to stop the effect of trypsin) and centrifuged. The pellet was resuspended in PBS and the cellular PpIX fluorescence was measured by flow cytometry (Guava® easyCyte; Cytek® Biosciences, USA). Excitation wavelength of 488 nm (blue laser) was used and the emission at 670 nm was determined (Red B). Flow cytometry data were analyzed employing FlowJo 7.6 and PpIX fluorescence histograms were compared. The percentage of positive cells for the photosensitizer fluorescence and the media fluorescence intensity were determined. Results are reported as the mean ± SEM. The experiment was made 3 times. Results of a representative experiment are shown.
2.8.3. D-Cultures Proliferation
T98G parental and resistant cells were seeded in P96 U-bottom wells coated with agarose (1% in deionized H
2O) [
13]. A total of 2,000 cells/well were seeded in DMEM medium 10% serum with erlotinib (50 µM; n=12). Controls without drug and vehicle controls with DMSO were included. The spheroids were photographed (Nikon DS-Qi1MC) 3 days after seeded and the number of cells by spheroids was counted employing a Neubauer chamber. Results are reported as the mean ± SEM. The experiments were made 3 times. Results of a representative experiment are shown.
2.9. Wound Healing Assay
Parental and resistant T98G cells were seeded (24 x 104 cells/ml) in P24. The next day, with a confluence of 100%, wounds were made in the cell monolayer using a micropipette tip (2 wounds per well). The cells were rinsed 3 times with a complete medium to remove lifted cells. Then, the cells were cultured with erlotinib at doses of 25 and 50 µM using culture medium with 1 and 10% serum (n=6 wells). Cells without erlotinib and without DMSO were used as controls. Vehicle-control groups with DMSO were included. When the control group closed the wounds, the cells were fixed with cold methanol and stained with toluidine blue (0.05%) for a better image visualization. Photographs were taken (Nikon DS-Qi1MC) at day 0 (t0) and final time (tf). The wound area was measured using ImageJ (n=8 wound areas). The percentage of wound closure was calculated using the formula: % migration= [100 - (wound area tf / wound area t0)] x 100. Results are reported as the mean ± SEM. The experiments were made 3 times. Results of a representative experiment are shown.
2.10. Statistical Analysis
The values in the figures are expressed as mean ± SEM. The R and R Studio software was used to carry out the statistical analyses. The statistical significance of the differences between the means was determined with the analysis of variance (ANOVA). When the means were found to be significantly different (with a p-value less than 0.05) multiple pairwise comparisons were performed using the Tukey HSD test.
4. Discussion
GBM is a grade IV brain tumor which represents one of the most lethal human cancers. Standard therapy for GBM encompasses surgical resection followed by chemotherapy and radiotherapy [
2]. A meta-analysis reported by Ren
et al., including 1,294 patients from 31 articles shows that the fluorescence guided resection (FGR) and PDT employing ALA are feasible for glioma treatment allowing a significant increase in the maximum resection rate and prolonging the survival after PDT [
14].
Advancement in precision oncology has renewed attention to the study of the receptor tyrosine kinase (RTK) EGFR as an important therapeutic target for GBM [
15]. Amplifications in the EGFR gene are detected in 57.4% of primary GBM patients, leading to high levels of EGFR protein and it has been observed that amplifications of EGFR are retained in recurrent gliomas [
16]. Repeated cycles of oncological treatments and cell growth in vitro can be performed in order to amplify the biochemical changes associated with cell resistance and to identify a selective target on surviving cells [
17]. We have found that PDT-resistant GBM cells from the T98G line have high levels of mRNA of EGFR, form spheroids with a greater number of cells in vitro and are more tumorigenic when injected into mice, compared to the parental population [
10].
Through the bioinformatics analysis by GEPIA that was shown in this work, it can be observed that high expression of EGFR in GBM patients can be related to poorer patient survival. It is known that EGFR is activated by ligands (EGF, EREG and EPGN) which induce dimerization structures and activates downstream effector pathways, including MAPK, PI3K, and STAT3, by auto-phosphorylation of tyrosine residues in its cytoplasmic domain [
15]. In GBM the activation of these pathways leads to tumor growth, invasion, angiogenesis and drug resistance to radio- and chemotherapy [
16].
GEPIA analysis also shows that compared to other cancer types, the highest expression of EGFR is observed in GBM and a high expression of this gene is found in GBM tumors compared to the low expression in its corresponding normal tissues. It has been described that the expression of EGFR in normal cells is about 4 × 104–10 × 104 receptors/cell, whereas, in cancer cells, more than 106 receptors/cell are observed. Except for hematopoietic cells, the majority of cell types express ErbB family members [
16]. It is important to consider the expression of tumor markers in normal cells for planning selective therapeutic strategies, which constitutes a fundamental challenge in oncology. Currently, there is a great impulse in the development of methods of selective drug delivery to tumors taking advantage of the differences in the microenvironment and morphology of tumor and healthy cells: effects on ion channels and the use of pH-sensitive intelligent formulations that allow selectively accumulating drug molecules in tumor cells, including the pH, membrane permeability, and ion channels [
18]. Using this approach, nanoparticles containing drugs can reduce side effects caused by drugs by encapsulating these drugs and releasing them at the target site. Researchers have examined different nanoparticle carrier categories such as polymers, lipids, and inorganic materials [
19].
In this work we evaluated the general response of T98G GBM cells resistant to PDT when EGFR activity is inhibited with the drug erlotinib. PDT in combination with erlotinib reduced viability in resistant populations compared to PDT alone, reaching viability levels similar to the parental cells. These results were related to an increase in the content of the photosensitizer PpIX in the resistant population when the cells were incubated with Me-ALA in combination with erlotinib. Fontana et al. found that EGFR inhibition increases the PpIX in GBM cells incubated with ALA because the expression of the enzyme heme oxygenase-1 (HO-1) is reduced. They showed that increased EGFR activity results in upregulation of HO-1-mediated heme clearance and ultimately reduced PpIX fluorescence. HO-1 accelerated heme depletion and thus a shift in enzymatic activity occurs in favor of increased PpIX metabolism by ferrochelatase. It is known that the expression of the HO-1 is primarily regulated at the transcriptional level by activating transcription factors such as NF-kβ, AP-2, and the heat shock-responsive element (HSE). EGF induces NF-kB activation through multiple EGFR dependent signaling molecules, including PI3K, protein kinase C (PKC), and IKK signaling pathways. So, EGF/EGFR signaling promotes HO-1 expression and activity in GBM cells through activation of the PI3K/AKT/NF-κB cascade [
20]. We could hypothesize that this axis which leads to a lower accumulation of the photosensitizer PpIX has been favored in the selective obtaining of resistant populations, among other resistance mechanisms. Then, the inhibition of EGFR activity by erlotinib leads to the increase of PpIX in resistant populations and to the consequent decrease in cell viability after PDT.
As we can see in the T98G 3D cultures, erlotinib inhibited the proliferation of both the parental T98G populations and the resistant populations. Also, erlotinib reduced wound closure in parental and resistant T98G using culture medium with low and high content of fetal bovine serum, suggesting that the decrease in EGFR activity affects cell migration and proliferation. A large body of work shows that EGFR inhibition reduces GBM cell migration, invasion, and proliferation. Half of GBMs that have EGFR amplifications overexpress the variant EGFRvIII, and expression of both receptors, EGFR and EGFRvIII, confers a worse prognosis. EGFR and EGFRvIII cooperate to promote tumor progression and invasion, in part, through activation of the Stat signaling pathway. Roos
et al. showed that EGFRvIII activates Stat5 and GBM invasion by inducing the expression of a mediator of glioma cell invasion and survival: fibroblast growth factor-inducible 14 (Fn14). EGFRvIII-mediated induction of Fn14 expression is Stat5 dependent and requires activation of Src, whereas EGFR regulation of Fn14 is dependent upon Src–MEK/ERK–Stat3 activation. Treatment of EGFRvIII-expressing GBM cells with the FDA-approved Stat5 inhibitor pimozide blocked Stat5 phosphorylation, Fn14 expression, cell migration and survival. Because EGFR inhibitors display limited therapeutic efficacy in GBM patients, the EGFRvIII–Stat5–Fn14 signaling pathway represents a node of vulnerability in the invasive GBM cell populations [
21]. Ji
et al. found that miR-615 expression was downregulated in GBM tissues and cells, and was inversely correlated with overall survival, relapse-free survival, tumor malignancy grade, and EGFR expression. They identified that miR-615 plays a tumor suppressor role in GBM cell proliferation, migration and invasion by targeting EGFR expression, and miR-615 may act as a novel biomarker for early diagnosis or therapeutic targets of GBM [
22]. Liu et al. evaluated the anti-tumor activity and effectiveness of the third-generation EGFR-targeted drug AZD9291 in vitro and in an orthotopic GBM model. AZD9291 inhibited significantly colony formation, migration, and invasion of GBM cells and reduced tumor growth and prolonged animal survival [
23].
Similar to T98G cells, in our work on PDT-resistant squamous cell carcinoma populations we found that PDT-resistant populations of the SCC-13 line accumulated less PpIX and were more tumorigenic than the parental cells [
24]. Furthermore, we determined that resistant SCC-13 populations have amplifications in the gene that encodes the mitogen-activated kinase MAPK/ERK, while these amplifications are not found in the parental SCC-13 population. These DNA changes were also observed at the mRNA and protein expression level and correlated with an increase in the expression and activation of EGFR. These results agree with those obtained by immunohistochemistry of human biopsies of squamous cell carcinomas resistant to two sessions of PDT with Me-ALA. Thus, our studies in squamous carcinoma propose that genomic imbalances related to EGFR and MAP3K1 may be an important factor in resistance to PDT [
25].
There are three generations of tyrosine kinase inhibitors (TKIs) approved for clinical use. First-generation TKIs (erlotinib, gefitinib, lapatinib and vandetanib) inhibit the receptor by reversible competitive binding with ATP in the tyrosine kinase domain. Subsequent generations of TKIs were created to overcome drug resistance. Second-generation TKIs (Afatinib, dacomitinib, tesevatinib) irreversibly inhibit all four ERBB receptors [
16]. The third-generation TKI (Almonertinib, Lazertinib, Furmonertinib, Osimertinib) is specifically designed to target the T790M resistance mutation [
26]. Efforts to target EGFR in glioblastoma continue to be investigated. Clinical trials currently active and in development for EGFR in GBM include the study of several novel agents, including the TKI epitinib, bi-armed activated T cells against EGFR/CD3, EGFR(V)-EDV-Dox, the monoclonal antibody GC1118, doxorubicin-loaded anti-EGFR immunoliposomes, and the anti-EGFRvIII/CD3 bispecific T-cell engager (BiTE) AMG 596, among others. Studies on the EGFR receptor in GBM are expected to lead to improved outcomes in this highly vulnerable patient population [
8].