Section 4
RPMI 1640 medium, fetal bovine serum (FBS), L-glutamine, amphotericin B (25 μg/mL), penicillin, streptomycin, trypsin EDTA solution, cell proliferation Kit I (MTT)—Thiazolyl Blue Tetrazolium Bromide, Calcein-AM (206700-1MG), and Rhodamine 123 (Rho 123) were obtained from Sigma (St. Louis, Missouri, USA). DPBS was bought from Corning (Corning, NY, USA). Cytotoxic drugs: PAC, CIS, DOX, and the inhibitor elacridar were sourced from Selleckchem (Houston, TX, USA). The anti-P-glycoprotein antibody (C219) (mouse monoclonal, Alx-801-002-c100) came from Enzo, (Farmingdale, NY, USA). The for GAPDH, β-actin, and goat anti-mouse horseradish peroxidase (HRP)-conjugated antibody were purchased from Proteintech, (Rosemont, IL, USA) (SA00001-1-A), and the fluorescent Alexa Fluor ® 488 AffiniPure™ Donkey anti-mouse IgG (H + L) (715-545-150, 1:200) was acquired from Jackson ImmunoResearch Laboratories (Cambridgeshire, UK). The RnaseOUT (10777019) and M-MLV reverse transcriptase kit (28025013) were obtained from Invitrogen by Thermo Fisher (Waltham, MA, USA). The materials used in the Western blot procedure (gels, protein marker, and membranes) were purchased from Bio-Rad (Bio-Rad Laboratories Ltd., Watford, Hertfordshire, UK).
The W1 human primary ovarian cancer cell line was established in our laboratory using ovarian cancer tissue obtained from an untreated patient as described previously [ 91 ]. Paclitaxel-resistant sublines (W1PR1, W1PR1-C7, and W1PR2) were created by gradually exposing W1 cells to increasing concentrations of PAC, eventually reaching a final concentration of 1100 ng/mL. These cells were grown as a monolayer in RPMI 1640 medium supplemented with 10% ( v / v ) FBS, 2 mM L-glutamine, 0.1 mg/mL penicillin, 100 U/mL streptomycin, and 25 µg/mL amphotericin B, incubated at 37 °C in a 5% CO 2 atmosphere.
RNA was extracted from all cell lines using the Gene Matrix Universal RNA Purification Kit (EURx, Ltd., Gdańsk, Poland) according to the manufacturer’s guidelines. The RNA concentration was determined by measuring absorbance at 260 and 280 nm. To synthesize cDNA, the reverse transcription was performed with 1.5 µg of RNA using the CFX Opus 96 Real-Time PCR System (Bio-Rad, Hercules, CA, USA). 1 μL oligodT18A (IBB PAN, Warsaw, Poland) and 1 μL dNTP Mix (ThermoFisher, Waltham, MA, USA, R0192) were added to each RNA sample, and denaturation followed at 65 °C for 5 min. Subsequently, 4 µL 5x First Strand Buffer, 2 µL DTT, 0.5 μL RnaseOUT, and 0.5 μL M-MLV RT were added. The reaction was carried out at 37 °C for 60 min, followed by 75 °C for 15 min.
Real-time PCR was conducted with the CFX Opus 96 Real-Time PCR System (Bio-Rad, USA). Each reaction sample contained: 1 μL cDNA solution, 12.5 μL Takyon™ ROX SYBR ® MasterMix blue dTTP (Eurogentec, Kaneka Eurogentec, Liège, Belgium), 1 μL of each sequence-specific primer (7.5 µM) from Oligo.pl (Warsaw, Poland) ( Table 6 ) and 9.5 μL water free of nucleases. A negative control was also included, using water free of nucleases in place of cDNA.
The real-time PCR was carried out under the following conditions: an initial denaturation step at 95 °C for 15 min, followed by 45 cycles of denaturation at 95 °C for 15 s, primer annealing at 60 °C for 30 s, and extension at 72 °C for 30 s. A final extension step was performed at 72 °C for 30 s. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) served as the housekeeping gene (internal reference control).
The results were analyzed with the CFX Maestro software (Bio-Rad v2.3 (5.3.022.1030)). Gene expression was assessed through the relative quantification (RQ) method. Gene expression levels for each sample were assessed relative to the 2D culture condition, which served as the calibrator (RQ = 1). The relative quantification was determined using the formula: RQ = [expression in sample (3D condition)/expression in calibrator (2D condition)]. Data visualization was performed with SigmaPlot software, version 15.0 (Systat Software GmbH, Erkrath, Germany).
All W1 cell lines at a concentration of 3000 cells per well were seeded into a 96-well plate in 200 µL of cell culture medium. After 48 h, the medium was replaced, and the cells were treated with elacridar (0 µM, 0.1 µM, or 1 µM) along with varying concentrations of PAC, DOX, or CIS drugs. Following a 72 h incubation period, an MTT assay was performed. The medium was replaced with 170 µL of fresh medium, and 100 µg of MTT reagent was added to each well, followed by a 1 h incubation. The resulting formazan crystals were dissolved in 200 µL of DMSO, and the absorbance was measured at wavelengths of 570 nm and 720 nm using the Synergy LX Multi-Mode Reader from BioTek Instruments, Inc. (Winooski, VT, USA) The IC 50 values were then determined from the results.
All W1 cell lines (15k cells per well) were seeded into non-adherent, 96-wells plates (BRAND plates inter Grade, F-bottom, 781902 (Merck, Darmstadt, Germany)) and maintained for 48 h. Following this initial incubation, 100 µL of the medium in each well was replaced with 100 µL of medium containing elacridar (final concentration: 0 µM, 0.1 µM, or 1 µM) together with the designated doses of PAC or DOX. Next, the spheroids were incubated for an additional 72 h. Cell viability was assessed using the MTT-based Cell Proliferation Kit I (Roche 11465007001). For the assay, 100 μL of medium was removed from each well, and replaced with 10 μL of the MTT labeling reagent. After a 4 h incubation, 100 µL of solubilization solution was added, and the plates were left overnight in the incubator to allow complete dissolution of the formazan crystals. Absorbance was measured (at 570 nm and 720 nm) using Synergy LX Multi-Mode Reader (BioTek Instruments, Inc.) IC 50 values were calculated based on the absorbance data.
Glass coverslips were placed into 24-well plates (one coverslip per well). Next, the cells were cultured in the wells at a density of 100,000 cells per well. Once cultures reached approximately 70% confluence, the cells were rinsed with phosphate-buffered saline (PBS), followed by fixation and permeabilization using a cold acetone/methanol mixture (1:1) for 15 min at room temperature (RT). After an additional PBS wash, nonspecific binding sites were blocked with 3% BSA for 30 min at RT. The samples were subsequently incubated with the primary antibody against P-gp (C219, 1:50) for 2 h at RT. After 3 PBS washes, the coverslips the Alexa Fluor ® 488-conjugated donkey anti-mouse secondary antibody for 1 h at RT. Following another series of PBS washes, the coverslips were mounted with a DAPI-containing medium (Sigma, F6057) for nuclei staining. P-gp protein expression was examined using a Zeiss Axio-Imager.Z1 fluorescence microscope (Oberkochen, Germany) and analyzed with Zen Blue v3.3 software.
All cell lines were seeded into 24-well plates at a density of 1 × 10 4 cells per well. After 72 h of growth, the cells were treated with the appropriate dosage of elacridar (final concentration 0.1 and 1 µM). After 1 h of incubation, the medium was changed to the medium with the addition of elacridar in proper concentration and Rh123 (1 µg/mL) and incubated 1 h. Next, the cells were washed three times with a cold solution of 50 µM Verapamil in PBS. Microscopic observation followed, and images were captured at 20× magnification for FITC and VIS channels (Leica DMi8 microscope, Leica, Germany).
The cells were seeded into non-adherent 96-well plates at a density of 5 × 10 4 cells per well and cultured for 5 days. The medium was then replaced with medium containing the appropriate amount of elacridar, and the cells were incubated for 1 h. Next, the spheroids were treated again with the same concentration of elacridar, with 0.25 μM Calcein-AM added for 1 h. The spheroids were washed three times with a cold solution of 50 µM Verapamil in PBS and observed under a microscope (Leica DMi8 microscope, Leica, Germany) taken at 40× magnification for FITC and VIS channels.
Cells were washed three times with PBS containing Ca/Mg ions and lysed using RIPA buffer (300 µL for 1 × 10 6 cells) with a protease inhibitor cocktail (Roche Diagnostics GmbH, Mannheim, Germany) for 80 min on ice at 4 °C. The lysates were then centrifuged at 13.4 × 10 3 rpm for 30 min at 4 °C. Protein concentrations were measured using the Bradford protein assay system (Bio-Rad Laboratories, Hemel Hempstead, UK). Samples for electrophoresis were prepared from 15 μg of protein mixed with 4× loading buffer (Bio-Rad Laboratories, Hemel Hempstead, UK) containing 10% β-mercaptoethanol. Proteins were separated by SDS-PAGE on a 4–20% mini-PROTEAN ® TGX™ precast gel and transferred to a nitrocellulose membrane using the Trans-Blot ® ™ transfer system (Bio-Rad Laboratories, Hemel Hempstead, UK). The membrane was then blocked with 5% milk in TBS/Tween (0.1 M Tris-HCl, 0.15 M NaCl, 0.1% Tween 20) and incubated overnight with primary antibodies against P-gp (Enzo, Alx-801-002-c100, 1:400), followed by a 3 h incubation with HPR-labeled mouse secondary antibody (Proteintech, SA00001-1-A, 1:10,000). The same membrane was subsequently washed with TBS/Tween (0.1 M Tris-HCl, 0.15 M NaCl, 0.1% Tween 20) for 1 h and incubated with the primary antibody against the reference protein, either β-actin (1:7500) or GADPH (1:40,000), followed by a 3 h incubation with HRP-conjugated secondary antibody. Signals were visualized using a chemiluminescence detection system (ECL, Femto Super Signal Reagent) and ChemiDoc™ (Bio-Rad Laboratories, Hemel Hempstead, UK).
P-gp efflux function after elacridar exposure was evaluated using Rhodamine 123 (Rho123) fluorescence. Cells were suspended at 1 × 106 cells/mL in culture medium and incubated with 0.1 or 1 µM elacridar for 1 h at 37 °C and mixed at 800 rpm. Rho123 was then added to a concentration of 1 µg/mL, and the samples were incubated for an additional hour under the same conditions. After incubation, the cells were placed on ice, centrifuged at 200× g for 5 min at 4 °C, and washed twice with ice-cold PBS supplemented with 50 μM verapamil, and immediately analyzed on a MACSQuant X flow cytometer (Miltenyi Biotec GmbH, Cologne, Germany). A total of 10,000 events were collected per sample, and fluorescence was detected at 488 nm. Data analysis was performed with FlowJo software (v10.9.0).
Intro
Ovarian cancer is the eighth most common cancer diagnosed in women. In 2020, it caused 4.7% of cancer deaths. Moreover, it is one of the most common gynecological cancers, and ranks third in terms of mortality, right after cervical and uterine cancer [ 1 ]. The reason for the high mortality rate is the late detection of cancer, which is a result of the initial stages being asymptomatic. Most cases are diagnosed in advanced stages, specifically stages III and IV, according to the FIGO classification, which accounts for 85% of diagnosed cases. Another cause of high mortality is low treatment effectiveness caused by primary or acquired resistance to chemotherapy [ 1 , 2 , 3 ].
Among ovarian cancers, epithelial ovarian cancer (EOC) is responsible for 90% of ovarian cancer cases. This type of cancer is characterized by a high degree of heterogeneity, which contributes to the histological and genetic complexity of this tumor. EOC is divided into two types that differ in tumor progression and genetic mutations. Type I tumors are low-grade tumors arising from endometriosis or borderline tumors, which include endometrioid, clear cell, mucinous, and low-grade serous carcinoma. Type II tumors are high-grade tumors with unidentified precursors, which include high-grade serous carcinoma, carcinosarcoma, and undifferentiated carcinoma [ 3 , 4 ].
The primary treatment for EOC involves cytoreductive surgery combined with chemotherapy [ 5 ]. The first line treatment comprises chemotherapy based on platinum drugs and taxane compounds (Paclitaxel—PAC). Depending on the response to platinum compounds, patients are divided into four main groups: platinum sensitive, partially platinum-sensitive, platinum-resistant, or platinum refractory [ 6 ]. Unfortunately, as many as 70% of patients experience the development of resistance and disease relapse. Patients with EOC resistant to platinum compounds are treated in the second line of chemotherapy with liposomal doxorubicin (DOX), gemcitabine, or topotecan (TOP) [ 3 , 7 ].
Low response rate to chemotherapy results from the presence of different mechanisms of drug resistance, which can be intrinsic or can develop during chemotherapy. Functionally, the mechanisms of drug resistance can be divided into those specific to cancer tissue and cancer cells specific ones [ 8 ].
The cancer-tissue-specific resistance is related to the tumor microenvironment and organization of the tumor tissue: high cell density in the tumor, the presence of blood vessels, the presence of other cells, like Cancer-Associated Fibroblasts (CAF), and non-cellular components (extracellular matrix (ECM), ECM-remodeling enzymes, and growth factors) [ 9 , 10 ]. Increased expression of ECM components is observed in chemotherapy-resistant cancer cells, which suggests an important role of ECM in drug resistance [ 11 , 12 ]. High tumor tissue density and the presence of ECM components may reduce drug availability by limiting drug diffusion or binding [ 13 , 14 ]. Moreover, ECM components can interact with cellular receptors (mainly integrins) and activate pro-survival signaling pathways, which is defined as cell adhesion-dependent drug resistance (CAM-DR). This resistance may also affect the upregulation of ABC drug transporters [ 15 , 16 ].
The active removal of cytotoxic agents from cancer cells is the most important mechanism of resistance at the cellular level. This activity results from the overexpression of ATP-binding cassette (ABC) transport proteins. ABC transporters have low substrate specificity and are thus able to remove different drugs, which are not related in chemical structure, from the cells. This phenomenon is designated as multidrug resistance (MDR) [ 17 ]. One of the most important ABC transporters is the MRP2 protein (ABCC2, encoded by the ABCC2 gene), which is responsible for resistance to several drugs, including cisplatin (CIS). BCRP (ABCG2) is another MDR protein that plays an important role in resistance to various drugs, such as anthracyclines, mitoxantrone, and TOP. The most important ABC transporter is P-glycoprotein (P-gp/ABCB1/MDR1, encoded by the MDR1 gene), which is responsible for the transport of various compounds [ 18 , 19 ].
P-gp belongs to the eukaryotic ABC superfamily of membrane transporters and is classified as a surface ATPase [ 20 ]. It is encoded by the MDR1 gene, which is located on chromosome 7q2, and consists of 29 exons [ 21 ]. The protein consists of 1280 amino acids and has a molecular mass between 140 and 150 kDa. The molecular mass of P-gp increases to 170 kDa after N-glycosylation of amino acids 91, 94 and 99 [ 22 , 23 ]. P-gp, as a transmembrane protein, has 12 transmembrane domains (TMD), two nucleotide-binding domains (NBD) and two ATP-binding domains located in the intracellular part of the protein. This protein does not have a well-defined ligand-binding pocket [ 24 ].
In physiological conditions, due to its transport capabilities, P-gp is expressed in epithelial cells, mainly with secretory and excretory functions, e.g., epithelium lining the trachea, mammary glands, prostate, biliary ducts, kidneys, intestine, and brain–blood barrier [ 25 ]. This protein plays an important role in liver cells, where it prevents xenobiotics or harmful chemicals from entering the cells [ 24 ]. P-gp is also responsible for the secretion of various endogenous and exogenous compounds into the bile [ 26 ]. Expression of P-gp protein in the intestine may reduce the bioavailability of some drugs [ 27 ].
Unfortunately, P-gp protein is also expressed in cancer cells. In some cancers, the expression can be primary. However, in most cases, expression of P-gp protein increases after chemotherapy [ 28 , 29 ]. High P-gp expression is observed in cancers such as breast, colon, lung, prostate, and ovarian cancer, leukemia, or glioblastoma [ 28 ]. Expression of P-gp and other ABC proteins in cancer stem cells (CSCs) also appears to influence the transport of signaling molecules contributing to tumorigenesis [ 30 ].
Overexpression of ABC transporters, including P-gp, in cancer cells causes the development of MDR and is associated with poorer patient prognosis [ 24 , 25 ]. The P-gp protein is responsible for removing many chemotherapeutic agents used in cancer treatment from the cells, such as anthracyclines (DOX, daunorubicin, epirubicin), actinomycin D, etoposide, mitomycin C, mitoxantrone, vinblastine, vincristine, taxanes (e.g., PAC), and TOP used in ovarian cancer treatment [ 24 ]. Among P-gp substrates used in chemotherapy of ovarian cancer, PAC stabilizes microtubules, preventing cell mitosis, in contrast, DOX intercalates between DNA’s double helix base pairs, resulting in the inhibition of replication and transcription [ 31 , 32 ].
The acquisition of MDR by tumors significantly complicates the treatment and, consequently, patient prognosis. Therefore, attempts are being made to overcome resistance to drugs. Among many different targets, ABC transporters are very promising. Compounds that inhibit ABC transporters or compete with chemotherapeutics for transport through ABC proteins are particularly sought after [ 20 ].
The inhibition of ABC protein activity causes increased accumulation of therapeutic agents inside cancer cells, which results in their concentration reaching therapeutic values [ 20 ]. One of the P-gp inhibitors is elacridar, which is a third-generation inhibitor of ABC transporters. The third-generation inhibitors are compounds synthesized specifically for the optimal inhibition of a given protein [ 24 ]. This inhibitor is a derivative of acridone carboxamide, a tricyclic acridine-based drug used in chemotherapy. It is also known as GG918, GF120918, or N-{4-[2-(1,2,3,4-tetrahydro-6,7-dimethoxy-2-isoquinolinyl)-ethyl]-phenyl}-9,10-dihydro-5-methoxy-9-oxo-4-acridine carboxamide. The molar mass of this compound is 563.64 g/mol [ 33 ]. Elacridar binding to P-gp most likely occurs in the transmembrane domains and competes for the binding site of [3H]azidopine [ 28 ]. This inhibitor is not specific to P-gp. It can also inhibit the activity of BCRP, another ABC transporter [ 34 , 35 ].
The effect of elacridar on drug sensitivity has been tested in many cancer cell lines and resulted in sensitization of cancer cells to drug treatment [ 36 ]. Drug-resistant prostate cancer cells with concomitant positive P-gp expression have been shown to regain sensitivity to docetaxel (DTX) when treated with elacridar [ 37 ]. Zhou et al. showed that elacridar sensitized highly resistant leukemia cells to daunorubicin and mitoxantrone [ 38 ]. Also, hepatoblastoma cells with increased P-gp expression increased their sensitivity to DOX after treatment with this inhibitor [ 39 ]. Elacridar was used to treat PAC-resistant human gastric cancer cells, resulting in sensitization of these cell lines to PAC [ 40 ]. This inhibitor also enhanced the efficacy of TOP in an ovarian cancer cell line and in cultured canine kidney cells overexpressing BCRP [ 35 , 41 ].
In preclinical studies, elacridar was studied in mice, rats, dogs, and monkeys [ 42 , 43 ]. The aim of these studies was to evaluate pharmacokinetic parameters, which showed good absorption and no effect on P450 enzymes [ 43 ]. Additionally, oral administration of taxanes with elacridar increased the concentrations of PAC and DTX in mouse plasma [ 43 ]. Jonker et al. showed that the bioavailability of TOP was increased by co-administration of elacridar in mice [ 44 ].
The use of elacridar in DOX co-treatment resulted in the reversal of resistance to DOX in P388/DOX tumors in a mouse xenograft model [ 45 ]. Moreover, administration of elacridar with levatinib to mice with hepatocellular carcinoma resulted in a better treatment response than administration of levatinib alone [ 46 ]. Omori et al. showed that combination therapy with elacridar and a topoisomerase inhibitor, Etoposide inhibits the growth of small cell lung cancer tumors in vivo [ 47 ]. Co-administration of elacridar and sunitinib inhibited tumor growth in sunitinib-resistant metastatic clear cell renal cancer in chicken chorioallantoic membrane and mouse models [ 48 ].
In clinical trials, the oncology patients with solid tumors were treated with a combination of DOX and elacridar. Increased exposure to doxorubicinol (a metabolite of DOX) was observed in patients treated with co-therapy compared to patients treated with DOX alone [ 49 ]. In other clinical studies, elacridar has been shown to significantly increase the systemic exposure to oral PAC in oncology patients [ 27 ]. Elacridar has also been shown to increase the bioavailability of TOP in cancer patients [ 44 ].
In the presented studies, PAC-sensitive and -resistant ovarian cancer cell lines were used. Previously, PAC-resistant cell lines W1PR1, W1PR1-C7, and W1PR2 were derived from the drug-sensitive W1 cell line. All PAC-resistant lines used in this study are characterized by high levels of P-gp expression [ 50 , 51 ] and are capable of growing in 2D and 3D cell culture conditions. The 2D model is a common model used in research. However, it does not reflect the conditions prevailing in the tumor, such as tumor structure, tumor cell packing, or overexpression of ECM molecules, which affect the treatment effect by lower drug diffusion. The 3D model allows for a better, although still not perfect, representation of the conditions inside the tumor [ 49 ]. Previous studies show that in a 3D model, cells exhibit greater resistance to chemotherapeutic drugs than those cultured in a monolayer [ 52 , 53 , 54 ].
The aim of the following studies was to investigate whether inhibition of P-gp activity by elacridar affects the resistance of cancer cells to cytotoxic drugs. To better understand the elacridar effect on drug resistance, we elucidated in 2D cultures a 3D model of cancer cell growth. This work provides a comprehensive analysis of P-gp activity, expression, and inhibition in PAC-sensitive and-resistant ovarian cancer cells.
Results
To study the potential of elacridar to reverse chemoresistance, we selected a model of the W1 ovarian cancer cell line and its PAC-resistant derivative. As the MDR1/P-gp overexpression plays a main role in PAC-resistance, we compared the expression of the MDR1 gene and P-gp protein in our model. MDR1 gene expression was determined using Q-PCR analysis, which showed a statistically significant increase in expression of this gene in PAC-resistant cell lines compared to sensitive cell lines ( Figure 1 A). We observed over 500-fold higher expression of the MDR1 gene in the W1PR1 cell line ( p < 0.001) and WIPR1-C7 cell line ( p < 0.01) and about 2000-fold higher expression in the W1PR2 cell line ( p < 0.001) than in sensitive cells.
The expression of P-gp protein in all cell lines was checked using the Western blot method ( Figure 1 B). We observed bands corresponding to P-gp protein in PAC-resistant cell lines, but they were not present in the W1 drug-sensitive cell line. The highest amount of P-gp protein was detected in the W1PR2 cell line.
Immunofluorescence staining was performed to confirm the presence of P-gp protein ( Figure 1 C). This assay confirmed the presence of P-gp protein in resistant cell lines with a clear localization of P-gp protein in the cell membrane. However, P-gp protein was not detected in the sensitive W1 cell line.
P-gp protein transport activity in W1 and PAC-resistant cell lines was examined by determining the accumulation of Rhodamine 123 (Rho123) using flow cytometry. The fluorescence intensity of accumulated Rho123 in PAC-resistant cell lines was lower than that in the sensitive cell line ( Figure 2 A–D). The lowest accumulation was observed in the W1PR2 and W1PR1-C7 lines, which is due to greater transport of Rho123 out of the cells. A slightly higher accumulation of Rho123 was observed in the W1PR1 cell line, but it is lower than in the PAC-sensitive line.
The accumulation of Rh123 in cells was also investigated using an inverted fluorescence microscope ( Figure 2 E). The intravital imaging of PAC-resistant cell lines did not show Rho123 accumulation. In contrast, in the W1 cell line, we observed a clear Rho123 fluorescence.
As we determined increased expression and activity of P-gp drug transporter in PAC-resistant cell lines, we were interested in whether P-gp expression/activity is reflected by the level of drug resistance. Thus, we used drugs that are the substrates of P-gp PAC and DOX, as well as CIS, a drug that is not a P-gp substrate. For all the cell lines, we determined the drug concentration that caused a 50% reduction in cell viability (IC 50 ). All PAC-resistant cell lines showed a significant increase in PAC resistance. In the W1PR1 cell line, we observed a 616-fold increase in IC 50 value, compared to the W1 cell line (IC 50 = 1960 ng/mL vs. IC 50 = 3.18 ng/mL). Lower resistance was demonstrated in the W1PR1-C7 line, which showed a 377-fold increase in resistance (IC 50 = 1199 ng/mL vs. IC 50 = 3.18 ng/mL), compared to the sensitive cell line. The highest PAC-resistance was observed in the W1PR2 cell line, with a 1064-fold increase in resistance in comparison to the W1 cell line (IC 50 = 3384 ng/mL vs. IC 50 = 3.18 ng/mL) ( Table 1 ).
DOX was the next cytotoxic drug that we analyzed. The IC 50 value for the sensitive cell line was 32.07 ng/mL. In all resistant cell lines, we observed an increased resistance to DOX compared to the W1 cell line. DOX resistance increased 215-fold for the W1PR1 line (IC 50 = 6904 ng/mL) and 86-fold for the W1PR1-C7 cell line (IC 50 = 2769 ng/mL). The highest IC 50 was determined in the W1PR2 cell line, which was 7836 ng/mL, which resulted in 244-fold higher resistance compared to the sensitive cell line ( Table 1 ).
The last cytotoxic agent tested was CIS, the drug used in the first line of ovarian cancer chemotherapy. We did not observe significant differences in resistance between the investigated cell lines. For the W1 cell line, the IC 50 was determined to be 19,741 ng/mL. PAC-resistant lines showed IC 50 of 21,566 ng/mL (W1PR1), 14,118 ng/mL (W1PR1-C7), and 16,779 ng/mL (W1PR2) ( Table 1 ).
To determine whether inhibiting P-gp activity would increase the sensitivity of cells to cytotoxic drugs, we used a selective P-gp inhibitor, elacridar, at concentrations of 0.1 µM and 1 µM. The first drug analyzed was PAC ( Figure 3 , Table 2 ). In the W1 drug-sensitive cell line, we did not observe any effect of elacridar on PAC-resistance ( Figure 3 A): IC 50 = 3.18 ng/mL without elacridar vs. IC 50 = 2.7 ng/mL with 0.1 µM elacridar and IC 50 = 2.65 ng/mL with 1 µM elacridar. In all PAC-resistant cell lines, we observed a strong effect of elacridar treatment on PAC resistance. We observed the effect of elacridar in both concentrations used on PAC resistance, starting from 2 ng/mL in the W1PR1 cell line ( Figure 3 B) and from 1 ng/mL in the W1PR1-C7 cell line, respectively ( Figure 3 C). This resulted in a decrease in IC 50 value. In the W1PR1 line, we observed a 404-fold decrease in IC 50 value in the presence of 0.1 μM elacridar (IC 50 = 4.85 ng/mL vs. IC 50 = 1960 ng/mL) and a 607-fold decrease in IC 50 value in the presence of 1 µM elacridar (IC 50 = 3.23 ng/mL vs. IC 50 = 1960 ng/mL). In the W1PR1-C7 cell line, we observed a 538-fold decrease in IC 50 value (IC 50 = 2.29 ng/mL vs. IC 50 = 1199 ng/mL in the presence of 0.1 µM elacridar) and a 620-fold decrease in IC 50 value in the presence of 1 µM elacridar (IC 50 = 1.93 ng/mL vs. IC 50 = 1199 ng/mL). In the W1PR2 cell line, a clear effect of elacridar in the concentration of 0.1 µM was observed in PAC concentration of 5 ng/mL and higher, and 1 µM of elacridar in the PAC concentration of 2 ng/mL and higher ( Figure 3 D). This was reflected by the 650-fold decrease in IC 50 value in elacridar concentration of 0.1 µM (IC 50 = 5.20 ng/mL vs. IC 50 = 3384 ng/mL) and a 945-fold decrease in IC 50 value for elacridar in concentration of 1 µM (3.58 ng/mL vs. IC 50 = 3384 ng/mL).
Next, we investigated the effect of elacridar on resistance to another P-gp substrate—DOX ( Figure 4 , Table 3 ). In the W1 cell line, we observed a slight, though not statistically significant, decrease in IC 50 value after elacridar treatment (IC 50 = 32.07 ng/mL vs. IC 50 = 25.0 ng/mL (0.1 µM E) vs. IC 50 = 21.6 ng/mL (1 µM E). For PAC-resistant cell lines, the elacridar treatment resulted in a reduction in resistance to DOX, and the effect was clearly dependent on the concentration used. Namely, in the W1PR1 cell line, elacridar in concentrations of 0.1 µM or 1 µM resulted in 82-fold (IC 50 = 6904 ng/mL vs. IC 50 = 84.6 ng/mL) or 152-fold (IC 50 = 6904 ng/mL vs. IC 50 = 45.5 ng/mL) decrease in IC 50 value, respectively, compared to DOX alone. The W1PR1-C7 line showed a 77-fold and 98-fold decrease in IC 50 value in the presence of 0.1 µM or 1 µM elacridar, respectively (IC 50 = 2769 ng/mL vs. IC 50 = 36.1 ng/mL vs. IC 50 = 28.4 ng/mL). The most significant effect of elacridar on DOX resistance was observed in the W1PR2 cell line, where we observed 123-fold and 436-fold reduction in IC 50 value in comparison to cells treated with DOX only (IC 50 = 7836 ng/mL vs. IC 50 = 63.8 ng/mL vs. IC 50 = 18.0 ng/mL).
The effect of elacridar and CIS co-treatment was also analyzed ( Figure 5 ). In both the sensitive and PAC-resistant cell lines, no statistically significant changes in cell viability were observed during co-treatment. For the W1 cell line, the IC 50 was determined to be 19,741 ng/mL (CIS alone), 16,355 ng/mL (CIS + 0.1 µM E), and 12,752 ng/mL (CIS + 1 µM E). The IC 50 determined for resistant cell lines was: W1PR1—IC 50 = 21,566 ng/mL vs. IC 50 = 20,099 ng/mL vs. IC 50 = 22,470 ng/mL; W1PR1-C7—IC 50 = 1418 ng/mL vs. IC 50 = 14,906 ng/mL vs. IC 50 = 14,998 ng/mL; W1PR2—IC 50 = 16,779 ng/mL vs. IC 50 = 20,541 ng/mL vs. IC 50 = 21,291 ng/mL ( Table 4 ).
Next, we performed an analysis of P-gp protein expression in PAC-resistant cells lines following elacridar treatment ( Figure 6 ). Cancer cells were treated with elacridar at different concentrations (0, 1, 2 or 5 µM) for 72 h. After this time, P-gp protein levels were analyzed by Western blot. The level of the tested protein did not change in any PAC-resistant line after elacridar treatment.
To check the activity of P-gp protein after treatment of cancer cells with elacridar, Rho123 fluorescence analysis using flow cytometry was performed. Rho123 fluorescence intensity assay in cells showed that elacridar does not affect the accumulation of the dye in PAC-sensitive cells ( Figure 7 A). In PAC-resistant cell lines, treatment with elacridar in concentration of 0.1 µM and 1 µM resulted in increased accumulation of the Rho123 dye with very similar effect in both concentration used ( Figure 7 B–D).
To further confirm the effect of elacridar on P-gp activity we analyzed the Rho123 accumulation in live cells using inverted fluorescence microscope. In W1 drug sensitive cell line we observed a similar Rho123 fluorescence in the presence and absence of elacridar in both concentration used. In contrast in all PAC-resistant cell lines Rho123 fluorescence was observed only after elacridar treatment and was independent of elacridar concentration used ( Figure 8 ).
To further investigate the effect of elacridar on drug resistance, a 3D cell culture model was used. It can be observed that the W1 cell line forms compact spheroids, while the PAC-resistant cell lines form a less compact spheroid structure ( Figure 9 B). First, the expression and activity of P-gp protein were checked in the studied cell lines. P-gp protein levels were compared in 2D and 3D cell cultures ( Figure 9 A). The sensitive cell line did not express the tested protein under any conditions. All PAC-resistant cell lines (W1PR1, W1PR1-C7, W1PR2) showed P-gp expression in both culture models used. The W1PR1-C7 cell line showed higher P-gp expression for 3D culture. However, the remaining PAC-resistant cell lines did not show any difference in the level of the tested protein depending on the cell culture model used.
The ability of cancer cells to accumulate Calcein-AM in a 3D model was also tested, which translates into the activity of P-gp proteins ( Figure 9 B). The dye was accumulated only in the PAC-sensitive cell line sphere. In contrast, in spheroids formed from PAC-resistant cell lines we did not observe CA accumulation.
The ability of 3D spheroids to accumulate Calcein-AM after elacridar treatment was tested ( Figure 10 ). In the PAC-sensitive cell line, we did not observe a change in the accumulation of this compound. At both elacridar concentrations, there is no visible change in Calcein-AM accumulation in the spheres formed from the W1 line in comparison to the control. In spheroids generated from PAC-resistant cell lines, no dye accumulation was observed in the control. In the presence of 0.1 µM elacridar, we observed Calcein-AM accumulations in spheroids from all three Pac-resistant cell lines. Increasing the inhibitor concentration to 1 µM did not significantly increase Calcein-AM accumulation.
The effect of elacridar on the response of spheroids to PAC treatment was assessed. The cytotoxicity of PAC was checked using the MTT method and determining the IC 50 for given conditions. Spheroids obtained from 1.5 × 10 4 cells were treated with 1 µM elacridar and increasing concentrations of PAC for 72 h. In the W1 cell line, no significant changes in IC 50 values were observed between cells treated with PAC alone (IC 50 = 4539 ng/mL) and cotreatment with elacridar (IC 50 = 4597 ng/mL) ( Figure 11 A, Table 5 ). However, for PAC-resistant cell lines, the presence of elacridar resulted in a reduction in the IC 50 value for PAC ( Figure 11 B,C, Table 5 ). Namely, for the W1PR1 cell line, the IC 50 value reduction was the greatest: 129-fold (decrease from IC 50 = 8405 ng/mL to IC 50 = 65.3 ng/mL). For the W1PR1-C7 cell line, a decrease in IC 50 was noted (from IC 50 = 4187 ng/mL to IC 50 = 53 ng/mL, which is a 79-fold reduction in value). The IC 50 value for the co-treated W1PR2 cell line decreased 73.5-fold (decrease from IC 50 = 5894 ng/mL to IC 50 = 78.3 ng/mL) ( Figure 11 D).
Discussion
Despite the enormous progress of cancer research that has occurred in over past 20 years, EOC statistics and prognoses for patients remain poor [ 55 ]. The nonspecific symptoms, late diagnosis, and eventual development of drug resistance continue to be obstacles to successful treatment [ 56 ]. As overexpression of ABC transporters plays a significant role in drug resistance, the use of ABC transporter inhibitors to tackle the issue of drug resistance is one of the approaches to improve effective treatment [ 57 ].
In our work, we have tested whether elacridar, a third-generation ABC transporter inhibitor, could be considered of use to combat PAC and DOX resistance. To study the impact of elacridar, we used a model of MDR development composed of the drug-sensitive W1 ovarian cancer cell line and its sublines—W1PR1, W1PR1-C7, and W1PR2—developed by the treatment with PAC [ 51 , 52 ].
As PAC and DOX resistance is most often associated with P-gp overexpression [ 17 ], we started by analyzing the expression of MDR1 transcript and P-gp protein expression in our model. We concluded that the expression of MDR1 /P-gp in PAC-resistant cell lines is much higher compared to the drug-sensitive cell line, and P-gp in drug-resistant cell lines is mainly localized in the plasma membrane, which is consistent with our previous studies [ 51 , 58 ].
To further investigate whether P-gp expression corresponds to increased activity in the researched cells, we have checked the cellular ability to retain the common fluorescent substrate of this protein, Rho123. This compound is a popular fluorone dye used to examine MDR1-related membrane transport in vitro [ 59 ]. Both the tests conducted on the Flow Cytometer and under the microscope concluded that Rho123 is accumulated only in the sensitive cell line, whereas it is actively removed by P-gp in all of the PAC-resistant cell lines. This is consistent with our previous works concerning different PAC-resistant ovarian cancer cell lines [ 60 ], as well as with the works of other authors [ 61 ].
As P-gp transports many cytotoxic drugs, our next step was to assess the resistance of the investigated cell lines to PAC and DOX, a commonly used P-gp substrate and a cytotoxic drug that is not a substrate of this protein—CIS [ 18 ]. We noted a significant increase in PAC and DOX resistance among all of the PAC-resistant cell lines when compared to the sensitive W1 cell line, with the lower resistance in the W1PR1-C7 cell line and the higher resistance in the W1PR2 cell line. The level of resistance seems to correlate with P-gp protein expression. A similar correlation was observed by us previously [ 60 ]. On the other hand, no significant change in CIS resistance was observed in the tested cell lines. This further confirms that drug resistance in this model is facilitated by P-gp. For further study, we used this model to investigate the effect of elacridar on P-gp activity and drug resistance.
Thus, using the MTT assay, we tested how elacridar in concentrations of 0.1 μM and 1 μM affects the PAC, DOX, and CIS resistance of the investigated cell lines. In case of PAC and DOX, we observed that the addition of elacridar in both concentrations caused a powerful decrease in PAC IC 50 value to a level similar to those in W1 drug-sensitive cell line, but it did not have a significant effect on PAC or DOX resistance in the W1 cell line. The effect of elacridar was concentration-dependent and increased the higher concentration was used. Similar results were observed by our team previously, in PAC-resistant sublines of A2780 cell culture [ 60 ]. In case of CIS resistance, elacridar had no relevant effect, regardless of the concentration used, and this was also observed by our team previously [ 60 ]. This was to be expected, as CIS is not a P-gp substrate [ 55 ].
In studies on chronic myeloid leukemia (CML) cell lines K562 and LAMA-84, as well as imatinib-resistant cultures, elacridar at a concentration of 0.25 μM effectively inhibited ABC transporter activity and led to a 5- to 10-fold increase in sensitivity to imatinib, as confirmed by flow cytometry [ 62 ]. Similarly, in the daunorubicin-resistant human promyelocytic leukemia cell line HL60/DNR, treatment with 0.1 μM elacridar restored drug sensitivity, resulting in a 40-fold reduction in IC 50 value for daunorubicin and a 57-fold decrease in IC 50 value for mitoxantrone [ 38 ]. In DTX-resistant human prostate cancer cell lines DU-145R and 22RV1, which show positive P-gp expression, the addition of 0.25 μM elacridar reversed resistance to DTX [ 37 ]. In non-small cell lung cancer (NSCLC) cell lines resistant to DTX, the combination of 0.25 μM elacridar with DTX significantly enhanced drug sensitivity, with IC 50 reductions ranging from 29-fold to over 3000-fold, depending on the cell line [ 63 ]. For the CIS-resistant ovarian cancer IGROVCDDP cell line characterized by P-gp overexpression, treatment with 0.25 μM elacridar increased sensitivity to multiple P-gp substrates—including DTX, PAC, epirubicin, and vinblastine—while it had no effect on drugs not transported by P-gp such as SN-38, 5-FU, and methotrexate [ 64 ]. In renal carcinoma 786-O cells, elacridar used at higher concentrations (2.5 μM and 5 μM) in combination with sunitinib enhanced cytotoxicity by inhibiting P-gp activity [ 65 ]. Moreover, elacridar at 0.5 µM concentration improved sunitinib efficacy in P-gp and/or BCRP overexpressing colorectal cancer cell lines (HT29, COLO320, DLD1, and CACO2) as well as patient-derived organoids [ 66 ].
Thus, we and others observed a strong impact of elacridar at a concentration of 1 µM or lower (in most cases) on resistance to cytotoxic drugs that are P-gp substrates. In our study, we observed a stronger effect in the W1PR2 cell line with the highest MDR1 /P-gp expression level, which suggests high effectiveness of elacridar in combating the resistance in cell lines with high P-gp expression level. A very similar effect was observed by our team previously [ 60 ]. This led us to investigate what is the mechanism of elacridar action in our model.
As a next step, we checked P-gp expression and activity in elacridar-treated cell lines. Although the protein expression following inhibitor treatment remained constant, there was a notable change in P-gp activity, which was confirmed both by flow cytometry and live-cell fluorescent microscopy, with Rho123 used as a fluorescent substrate for the protein. Similarly, we previously observed that elacridar inhibits the activity of another ABC protein—BCRP, without any changes in protein expression [ 35 ]. Thus, our results are consistent with the results of other investigators and suggest that the main mechanism of elacridar action is the inhibition of ABC transporters [ 27 ].
This indicates that elacridar successfully inhibits P-gp activity in our 2D cell culture model and is a good candidate for a drug that can increase the effectiveness of chemotherapy. However, a tumor is not composed of single cancer cells but is a complicated three-dimensional (3D) structure composed of cellular and non-cellular components that effectively limit drug diffusion, bind some anticancer drugs, and induce drug resistance in the CAM-DR mechanism [ 9 , 10 , 13 , 14 ].
Thus, we conducted similar experiments for the cells grown in the 3D model. First, we compared the expression of P-gp protein in the same cell lines growing as 2D monolayer and 3D spheroids, and we observed a very similar P-gp expression level in both cell culture conditions in W1PR1 and W1PR2 cell lines. In contrast, increased P-gp expression in 3D conditions was observed in W1PR1-C7 cell line. A similar level of MDR1 transcript between cells growing in 2D and 3D cell culture conditions and expression of P-gp protein in 3D spheroids were observed by our team previously [ 54 , 67 ]. These results indicate that overexpression of P-gp protein is a universal feature of PAC-resistant cell lines, regardless of cell culture conditions.
Next, we performed a fluorescent microscopy analysis in 3D spheroids with Calcein-AM as a P-gp substrate. Calcein-AM is a non-fluorescent compound that undergoes hydrolysis in cells, turning into a fluorescent calcein [ 68 ]. Unlike Rho123, this compound becomes fluorescent once it is inside the cell, resulting in a significantly reduced background signal [ 69 ]. This is particularly important when working with spheroids, as they are fragile and often become damaged during washing off the excess substrate. The effect of elacridar on Calcein-AM accumulation in 3D spheres was similar to the effect on Rho123 accumulation in the 2D model. The W1 cell line with no P-gp expression retained the substrate. In contrast, the PAC-resistant cell lines did not retain the substrate, suggesting high P-gp activity in these spheroids. Very similar observations were made by us previously [ 60 ]. This experiment also indicates that Calcein-AM can easily diffuse into dense spheroids like spheroids formed by the W1 cell line, making it a good fluorochrome to study P-gp activity in a 3D cell culture model.
Next, we investigated the effect of elacridar on P-gp activity in 3D spheroids. No further increase in elacridar accumulation was observed in the drug-sensitive cell line. However, following elacridar treatment, every drug-resistant cell line showed increased substrate retention. The increase in Calcein-AM accumulation in 3D spheroids indicates that elacridar easily diffuses in the spheroid structure, and makes it a good candidate for increasing the effectiveness of chemotherapy based on P-gp substrates.
MTT tests carried out in the 3D model showed a few-fold increase in resistance to PAC in 3D cell culture conditions in comparison to 2D cell culture conditions in all three PAC-resistant cell lines. However, over a thousand-fold increase in PAC resistance was observed in the W1 cell line in 3D spheroids in comparison to a 2D monolayer. Elacridar treatment did not change the IC 50 value in W1 3D spheroids. In contrast, elacridar was able to reverse PAC-resistance dramatically in PAC-resistant cell lines. Although the method chosen to determine the IC 50 value is not perfect and allows for determining only approximate IC 50 values due to the a very strong reduction in cell viability at a PAC concentration of 100 ng/mL in the presence of elacridar at a concentration of 1 µM, the obtained results indicate an effective decrease in the IC 50 value in the presence of elacridar under 3D culture conditions in the resistant cell lines. However, the IC 50 value in 3D spheroids was still over 20-fold higher than the IC 50 value in the same cell line in a 2D monolayer following elacridar treatment. In summary, after elacridar treatment, the spheroids from drug-resistant cell lines were over 20-fold more resistant than the cells grown in a monolayer, and spheroids from W1 cell line were over a thousand-fold more resistant than the same cells grown in a monolayer. This indicates that other factors can be responsible for drug resistance in 3D spheroids. Looking at the structure of spheroids created from investigated cell lines, it is easy to distinguish two different kinds of spheroids. The regular, very dense spheroids that comprise the outer proliferation zone, inner quiescent zone and necrotic zone present in the biggest spheroids [ 70 ] were formed by W1 cells. In contrast, all PAC-resistant cell lines formed low-density spheroids. The very high increase in resistance of the W1 cell line between 2D and 3D cell culture conditions seems to result from the very high density of W1 spheroids, which effectively limits the diffusion of PAC, which is a big molecule with very limited diffusion properties [ 71 , 72 , 73 , 74 ]. This makes the high-density spheroids much more resistant to PAC action, as was observed by us previously [ 54 , 67 ] and is supported by studies of other investigators in breast cancer spheroids [ 75 , 76 ].
Another factor that can increase resistance to cytotoxic drugs is overexpression of extracellular matrix (ECM) molecules and expression changes in other drug-resistant genes, as it was observed by others [ 77 ] and by our team [ 74 , 78 , 79 , 80 , 81 ] in drug-resistant cell lines. It is widely recognized that the high expression of the ECM components in tumors restricts the penetration of therapeutic agents [ 82 , 83 ]. This effect is further intensified when certain cytotoxic drugs bind to cellular macromolecules, ultimately reducing the amount of drug available to target cancer cells [ 82 ]. Among the investigated drug-resistant cell lines, we previously observed a strong expression of COL3A1 gene and its corresponding protein in the W1PR2 cell line, which seems to be related to PAC resistance [ 54 , 71 ].
Beyond acting as a physical barrier, the ECM also contributes to drug resistance through CAM-DR [ 82 , 83 , 84 ]. In this process, interactions between cancer cells and ECM proteins trigger intracellular signaling cascades and induce the expression of pro-survival genes, leading to increased resistance to drug-induced apoptosis [ 82 , 83 ]. This mechanism seems to be especially important in 3D spheroids, where cells are in close contact.
Increased intracellular signaling leading to increased expression of drug-resistant genes can also be related to decreased expression of some genes in drug-resistant cell lines. Previously, we observed decreased expression of some genes in PAC-resistant cell lines, among them: SEMA3A (Semaphorin 3 A) [ 79 ] and PTPRK (Receptor-type tyrosine-protein phosphatase kappa) [ 85 ]. SEMA3A is a tumor suppressor gene, and its decreased expression was observed in gastric and ovarian cancer and correlated with disease progression and poor prognosis [ 86 , 87 ]. In contrast, overexpression of SEMA3A correlated with inhibition of Src and FAK (focal adhesion kinase) kinases phosphorylation in tongue squamous cell carcinoma cell line SSC-9 [ 88 ]. PTPRK is a phosphatase, and its decreased expression was observed in breast cancer and correlated with the disease progression [ 89 ]. Mutation in the PTPRK gene correlated with chemoresistance in glioma [ 90 ]. Previously, we observed downregulation of PTPRK in fifteen ovarian cancer cell lines [ 85 ] and restoration of its expression by piperine treatment correlated with increased sensitivity to PAC [ 58 ]. Thus, the decreased expression of SEMA3A and PTPRK can be another factor that increases resistance in PAC-resistant spheroids.
Another aspect that should be discussed is the difference in the PAC response curves between spheroids formed from W1 and drug-resistant cell lines. In spheroids formed from PAC-resistant cell lines, we observed a two-step response curve in the absence of elacridar. These spheroids remained resistant even at very high PAC concentrations. Only when the PAC concentration exceeded the transport capacity of P-gp did cell death occur, which was reflected by a sudden collapse of the response curve. In contrast, after blocking of P-gp activity with elacridar, a sudden collapse of the response curve was observed at low PAC concentrations. However, some cells still survived, which was reflected by a flat response curve in the concentration range of 100 ng/mL to approximately 5000 ng/mL. Thus, about 50-fold increase in PAC concentration did not change spheroid viability. Furthermore, the level of viability seems to be related to spheroid density, being lowest in the W1PR1-C7 cell line (5–10% of control) and the highest in W1PR2 cell line (approximately 35% of control).
In contrast, in W1 spheroids, we observed a three-step response curve. At PAC concentration of 100 ng/mL, cell viability remained above 60% of the control, and did not change up to 2000 ng/mL (over a 20-fold increase). Only higher PAC concentrations were able to kill all cells within the spheroids.
These responsive curves appear to be associated with distinct cellular zones within the spheroids. In the outer zone proliferating cells are present, and this region is easily distinguishable in W1 spheroids. Since PAC is an antimitotic agent, it primarily kills proliferating cells. As most cells in drug-resistant spheroids are highly proliferative, blocking P-gp activity with elacridar results in PAC-induced cell death in these spheroids. In contrast, the inner zone contains quiescent cells that are much more resistant to PAC. This zone is considerably larger in W1 spheroids, which is reflected by their higher viability at increased PAC concentrations. However, very high PAC concentrations can also kill quiescent cells, likely through a different mechanism.
In conclusion, the conducted experiments suggest that elacridar, a dual P-gp and BCRP inhibitor, is a promising target for further research on combating MDR in ovarian cancer. The conducted research demonstrates its effectiveness in 2D and 3D in vitro conditions, by making the cells more susceptible to PAC and DOX. However, further research in more complicated models is necessary to further assess the effectiveness of elacridar as a drug in chemotherapy.
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.