A synthetic lethal screen for Snail-induced enzalutamide resistance identifies JAK/STAT signaling as a therapeutic vulnerability in prostate cancer

Despite substantial improvements in the treatment landscape of prostate cancer, the evolution of hormone therapy-resistant and metastatic prostate cancer remains a major cause of cancer-related death globally. The mainstay of treatment for advanced prostate cancer is targeting of androgen receptor signaling, including androgen deprivation therapy plus second-generation androgen receptor blockade (e.g., enzalutamide, apalutamide, darolutamide), and/or androgen synthesis inhibition (abiraterone). While these agents have significantly prolonged the lives of patients with advanced prostate cancer, is nearly universal. This therapy resistance is mediated by diverse mechanisms, including both androgen receptor-dependent mechanisms, such as androgen receptor mutations, amplifications, alternative splicing, and amplification, as well as non-androgen receptor-mediated mechanisms, such as lineage plasticity toward neuroendocrine-like or epithelial-mesenchymal transition (EMT)-like lineages. Our prior work identified the EMT transcriptional regulator Snail as critical to hormonal therapy resistance and is commonly detected in human metastatic prostate cancer. In the current study, we sought to interrogate the actionable landscape of EMT-mediated hormone therapy resistant prostate cancer to identify synthetic lethality and collateral sensitivity approaches to treating this aggressive, therapy-resistant disease state. Using a combination of high-throughput drug screens and multi-parameter phenotyping by confluence imaging, ATP production, and phenotypic plasticity reporters of EMT, we identified candidate synthetic lethalities to Snail-mediated EMT in prostate cancer. These analyses identified multiple actionable targets, such as XPO1, PI3K/mTOR, aurora kinases, c-MET, polo-like kinases, and JAK/STAT as synthetic lethalities in Snail+ prostate cancer. We validated these targets in a subsequent validation screen in an LNCaP-derived model of resistance to sequential androgen deprivation and enzalutamide. This follow-up screen provided validation of inhibitors of JAK/STAT and PI3K/mTOR as therapeutic vulnerabilities for both Snail+ and enzalutamide-resistant prostate cancer.


Introduction
The treatment landscape of prostate cancer exemplifies the "two truths" of cancer treatment (Somarelli et al., 2022): While tremendous progress has been made to improve patient outcomes, there also remains an urgent need to overcome the significant challenges imposed by the evolution of treatment resistance and metastasis. From the groundbreaking studies of Huggins and Hodges (Huggins and Hodges, 1972) to the development of novel, second-generation androgen receptor inhibitors (Scher et al., 2012;Beer et al., 2014;Smith et al., 2018;Chi et al., 2019;Fizazi et al., 2019;Armstrong et al., 2022), and anti-androgens (Fizazi et al., 2017;James et al., 2017), much of the existing treatment options for prostate cancer are currently focused on targeting the androgen receptor (AR) signaling axis. These agents have demonstrated significant clinical benefit; however, progression of men treated with these agents in the metastatic, castration-resistant setting is nearly universal.
The evolution of resistance to AR signaling inhibitors is mediated by heterogeneous genetic and non-genetic pathways that include both AR-dependent and AR-independent mechanisms [reviewed in (Blatt and Raj, 2019)]. Among these heterogeneous mechanisms, phenotypic plasticity is a central hallmark of AR signaling inhibitor resistance (Somarelli et al., 2020a). This phenotypic plasticity occurs along multiple, interconnected cellular lineage axes, such as stemness (Alumkal et al., 2020;Rodriguez et al., 2022), epithelial/mesenchymal Miao et al., 2017;Bai et al., 2019;Ware et al., 2020), luminal/basal (Stoyanova et al., 2013;Li et al., 2018), and neuroendocrine-like lineages or cell states (Ku et al., 2017;Mu et al., 2017). Phenotypic plasticity along these axes often leads to a loss of AR expression/activity and dependency (Formaggio et al., 2021), as well as additional aggressive features that promote survival and metastasis (Schroeder et al., 2014;Labrecque et al., 2021). New approaches are needed to capitalize on these emerging phenotypic states for therapeutic benefit.
Targeted therapy alters the ecological fitness landscapes of cancer in multiple ways (Somarelli, 2021). The altered fitness landscape of the drugged environment can promote aggressive biology, but can also induce "collateral sensitivities" to novel agents (Acar et al., 2020). This concept, also known as negative cross resistance, has been applied to identify new strategies to treat the evolution of resistance in bacterial infections (Imamovic and Sommer, 2013), malaria (Kirkman et al., 2018), herbicides (Cutti et al., 2021), and pesticides (Wazir and Shad, 2022).
In the present study, we combined high-throughput screens with multiparameter endpoint measurements from transcription-based reporters, confluence, and cell viability assays to characterize the therapeutic landscapes of Snail-mediated EMT, enzalutamide resistance, and AR activity ( Figure 1A). Our analyses pinpoint histone deacetylases (HDAC), protein kinase A (PKA), PI3K/ mTOR, and Janus Kinase (JAK) as key collateral sensitivities to Snail-mediated enzalutamide resistance in prostate cancer cells. Follow-up screens in a model of progressive adaptation to ADT and enzalutamide resistance verified the relevance of these pathways as novel therapeutic vulnerabilities for enzalutamide-resistant prostate cancer ( Figure 1B). These analyses provide a deeper understanding of the therapeutic vulnerabilities induced by epithelial plasticity and enzalutamide resistance.

Materials and methods
Cell culture models LNCaP95-Snail (Hu et al., 2012) and CS2 enzalutamide-resistant cells Jindal et al., 2023) were cultured in RPMI containing 10% charcoal stripped Fetal Bovine Serum (Sigma, St. Louis, MO) and 1% penicillin/streptomycin (Life Technologies, Carlsbad, CA). CS2 enzalutamide-resistant cell populations were maintained in the presence of 50 μM enzalutamide (Selleckchem, Houston, TX). Cell lines were maintained in standard tissue culture-treated plasticware within a humidified incubator at 37°C and 5% CO 2 . LNCaP95 cells stably expressing inducible Snail were generated as previously described . Induction of Snail nuclear translocation was mediated by the addition of 4-hydroxy-tamoxifen (4OHT) at a concentration of 20 nM. Ethanol (EtOH) was used as a vehicle control. All cells were authenticated by the Duke DNA Analysis Facility using analysis of short tandem repeats and were verified to be mycoplasma-free.

Development and testing of MET and PSA reporter lines
We adapted the GIIIcI 2 MET reporter (Somarelli et al., 2013;Somarelli et al., 2016) for lentiviral transduction by cloning the previously-described vector into the lentiviral vector pLVX-puro using restriction enzymes EcoRI/SmaI (NEB, Ipswich, MA). The PSA reporter was synthesized in the lentiviral expression plasmid, pLV[Exp]-Puro by VectorBuilder (Chicago, IL) to include 2 Kb of the proximal PSA promoter upstream of the enhanced green fluorescent protein (EGFP) open reading frame (Zhang et al., 1996). Cells stably expressing inducible Snail (Addgene plasmid #18798) or indicated reporter plasmids were generated by transduction of LNCaP95 or CS2 cells as described: https://www.addgene.org/ protocols/generating-stable-cell-lines/. Confluence and fluorescence were measured with and without EMT induction using Snail activation as described above. For PSA-GFP expressing cells, confluence and fluorescence was quantified with and without AR activation using synthetic androgen R1881(Sigma, St. Louis, MO) at 1 nM.

High-throughput drug screening
High-throughput screens were performed in collaboration with the Duke Functional Genomics Shared Resource as previously described Somarelli et al., 2020b;Somarelli et al., 2020c;Rao et al., 2020). Briefly, compounds from the Bioactives library (SelleckChem, Houston, TX) were stamped in triplicate into 384-well plates at a final concentration of 1 μM using an Echo Acoustic Dispenser (Labcyte, Indianapolis, IN, United States). Cells and media were subsequently dispensed into plates using a WellMate (Thermo Fisher, Waltham, MA, United States) at a density of 2,000 cells/well for each cell line. Confluence was quantified using an IncuCyte S3 (Gottingen, Germany) live cell imaging system. GIIIcI 2 and PSA-GFP readouts were quantified by IncuCyte imaging at 24, 48, 72, and 96 h. CellTiter Glo (Promega, Madison, WI) was added at 96 h, and luminescence was read using a Clariostar plate reader (BMG, Berlin, Germany).

RNA-Seq analysis of EMT scores
Quantification of EMT status for each sample was done using three distinct methods, 76GS, KS, MLR, each of which uses a unique algorithm and gene set. The 76GS scores were calculated based on the expression of 76 genes (Byers et al., 2013). Higher scores correspond to more epithelial states. A 76GS score >0 typically indicates an epithelial phenotype and <0 indicates a mesenchymal phenotype. The score for each sample is computed as the weighted sum of expression values of 76 genes, with the weight factor being the correlation of expression values of that gene with that of CDH1 in the given dataset. KS score was determined based on a Kolmogorov-Smirnov two-samples test (Tan et al., 2014). Using a 218 gene signature, the cumulative distribution functions are estimated for mesenchymal and epithelial signatures, and the maximum difference in cumulative distribution functions is retained as the statistic for the two sample-KS test. KS score ranges from [−1, 1], with negative and positive scores representing mesenchymal and epithelial phenotypes, respectively. MLR scores are provided on a scale of [0, 2]; higher scores are associated with more mesenchymal samples (George et al., 2017). Using an ordinal multinomial logistic regression, the score encompasses an order structure, with a hybrid epithelial/ mesenchymal signature situated between the epithelial and mesenchymal phenotypes. Scores are calculated based on the probability assigned for each sample to belong to one of the three phenotypes.

Data analysis
The primary objective for the high-throughput screen of LNCaP95-Snail cells was to identify synthetic lethality for Snail+ cells. Snail−cells (EtOH-treated vehicle controls) were used as a reference control to calculate differential effects across all parameters. The primary objective for the high-

Frontiers in Molecular Biosciences
frontiersin.org throughput screen of CS2 enzalutamide-resistant cells was to identify collateral sensitivities for enzalutamide-resistant cells.
The central hypothesis for this work was that activation of key pathways in Snail+, enzalutamide-resistant prostate cancer can be exploited for therapeutic benefit through synthetic lethal and collateral sensitivity approaches. Experimental data were visualized and analyzed in GraphPad Prism 9 (Boston, MA). Analysis of cell viability by CellTiter Glo was performed by normalizing to the average of all empty (non-drug) wells. Quantification of immunofluorescence images was performed using ImageJ. Briefly, images were converted to 8-bit, adjusted for threshold, and analyzed for particle count and area. For nuclear proportion experiments, nuclear Snail was counted and normalized to total number of cells per field. Imaging of confluence and GFP was compared using repeated measures ANOVA. Linear regression was used to assess correlations between screen analysis parameters, and outliers were considered to fall outside the 95% confidence interval bands. p-values < 0.05 were considered statistically reliable.

Fluorescence-based reporters enable realtime monitoring of epithelial plasticity
Prior studies have pinpointed the epithelial plasticity regulator, Snail, as both upregulated during AR inhibition (Miao et al., 2017) and a mediator of enzalutamide resistance through sustained androgen receptor signaling . LNCaP95 cells are a model of castration-resistant prostate cancer that demonstrate sustained AR expression, cellular plasticity (Hu et al., 2012) and acquire enzalutamide resistance through activation of the transcription factor Snail . In the present work we sought to develop a system to identify novel synthetic lethality to Snail-induced resistance to enzalutamide. To do this we turned to a Snail inducible LNCaP95 cell line system in which Snail is fused to an estrogen receptor mutant (ER mut ) whereby 4-hydroxytamoxifen (4OHT) acts as an agonist ( Figure 2A). Addition of 4OHT induces estrogen receptor-Snail fusion nuclear localization and activation of Snail ( Figure 2A). Addition of 4OHT in the Snail-inducible LNCaP95 prostate cancer cell line leads to cell scattering and upregulation of the mesenchymal marker, vimentin ( Figure 2B; Supplementary Figure S1). To track dynamics of Snail-mediated epithelial plasticity we adapted the GIIIcI 2 fluorescence-based reporter (Somarelli et al., 2013) for lentiviral transduction. The GIIIcI 2 reporter utilizes the lineage-specific alternative splicing within the ligand binding domain of FGFR2 to control EGFP (Somarelli et al., 2013) expression based on epithelial or mesenchymal phenotype. The EGFP open reading frame is interrupted by the FGFR2-IIIc exon and flanking introns ( Figure 2C). Splicing of FGFR2-IIIc in epithelial cells leads to fusion of the EGFP reading frame and subsequent EGFP expression while inclusion of the IIIc exon interrupts the EGFP reading frame and prevents EGFP expression ( Figure 2C). Treatment of LNCaP95-Snail cells with 4OHT leads to a modest reduction in confluence, consistent with the known relationship between Snail and cell cycle arrest (Vega et al., 2004) (Figure 2D). Similarly, Snail induction also induces robust inhibition of EGFP expression ( Figure 2E) consistent with inclusion of the mesenchymal FGFR2-IIIc exon. A loss of EGFP signal in Snail+ cells is also evident by fluorescence imaging of Snail− (EtOH) and

Frontiers in Molecular Biosciences
frontiersin.org Snail+ (4OHT) cells ( Figure 2F). EGFP expression from the GIIIcI 2 reporter is also consistent with endogenous FGFR2 splicing, in which 4OHT induces a switch from the IIIb to IIIc isoforms, as observed by isoform-specific restriction digestion of FGFR2 RT-PCR products ( Figure 2G).

High-throughput screens identify synthetic lethality to Snail-induced epithelial plasticity
We applied this Snail-inducible plasticity reporter system to identify compounds with synthetic lethality for Snail+ prostate cancer that could be subsequently validated for activity in models of enzalutamide resistance given the association between Snail expression and enzalutamide resistance . To do this, we performed a high-throughput small molecule screen using the SelleckChem Bioactives compound library. The Bioactives library contains 2,100 small molecules annotated by target and pathway. The library was designed to include compounds that are structurally diverse, medicinally active, and cell permeable, including both FDA-approved and nonapproved compounds Somarelli et al., 2020b;Somarelli et al., 2020c). Screen results were analyzed for cell viability/ATP production by CellTiter Glo at the fourday endpoint, and for cell growth rate and epithelial plasticity status by daily IncuCyte imaging of confluence and GIIIcI 2 EGFP levels, respectively, for 4 days ( Figure 3A; Supplementary Table  S1). Analysis of CellTiter Glo values for empty wells revealed a significant reduction in growth for Snail+ cells (Supplementary Figure S2A), which is consistent with the known role of Snail as a mediator of cell cycle arrest. Across the entire compound library 3.8% of compounds inhibited CellTiter Glo signal for Snail− cells by 50% or more, while 22% of the library inhibited Snail+ cells 50% or more (Supplementary Table S1). To identify compounds with differential sensitivity based on Snail expression, we analyzed the differential sensitivity of Snail− and Snail+ cells to all compounds in the library, with a 1.0 representing no difference in sensitivity. Drugs with values <1.0 differentially inhibit CellTiter Glo output of Snail+ cells while drugs with values >1.0 differentially inhibit CellTiter Glo output in Snail− cells ( Figure 3B).
In parallel to CellTiter Glo, we also quantified differences in growth rate for all screen compounds with and without Snail induction. Cell confluence was moderately, but significantly, correlated with CellTiter Glo values when comparing all treatment conditions (Supplementary Figure S2B). To identify collateral sensitivities based on growth rate we first calculated differences in slope of the growth rates between Snail− (EtOH) and Snail+ (4OHT) cells. This analysis is shown for a subset of compounds in Supplementary Figure S2C, with compounds in gray having little to no effect on cell growth of Snail− (EtOH) cells and these same compounds inhibiting growth in Snail+ (4OHT) cells. Subsequent annotation by target enabled identification of targets for which >2 drugs hit the same target. Top hits were ranked by their differential slope when comparing Snail+ to Snail−cells. Among these hits were inhibitors targeting signaling molecules and Overlapping drugs with differential sensitivity in Snail+ cells for both confluence and CTG assays. (F) Candidate EMT/MET inducers ranked by GIIIcI 2 induction (higher GFP = more epithelial; lower GFP = more mesenchymal). (G) Top 10 candidate MET inducing compounds, as estimated by EGFP expression from the GIIIcI 2 reporter.

Frontiers in Molecular Biosciences
frontiersin.org pathways known to be involved in lineage plasticity and prostate cancer therapy resistance, such as aurora kinase, c-MET, and mTOR/PI3K ( Figure 3C). Other targets included inhibitors of CRM1 (XPO1), a nuclear shuttling protein, cyclin-dependent kinases, polo-like kinases, and protein kinase C ( Figure 3C). To identify synthetic lethality for Snail+ cells, we focused on agents with <50% killing in Snail− (EtOH) cells and >50% killing in Snail+ cells by CellTiter Glo. Among these compounds, comparison of drugs that inhibited both CellTiter Glo production and growth rate by greater than 2-fold in Snail+ cells as compared to Snail− cells revealed four candidate compounds ( Figure 3D), including ONX-0914 (immunoproteasome inhibitor), AZ-960 (JAK2 inhibitor), WHI-P154 (JAK3 and EGFR inhibitor), and CCT137690 (aurora kinase inhibitor) ( Figure 3E). We next attempted to identify compounds and pathways that inhibit Snail-induced EMT. To do this we first calculated the fold change in EGFP expression for each compound at day 4 as compared to day 1. The fold change in EGFP expression for 4OHT (Snail+) cells was divided by EtOH (Snail−) cells for each compound to identify drugs that were capable of overcoming Snailmediated EMT. To ensure the gain in EGFP expression was not simply a function of cell growth inhibition or cell death, we next compared the EGFP expression to the differential confluence in 4OHT-treated versus EtOH-treated cells. This analysis revealed a subset of compounds that led to differential re-activation of EGFP expression from the GIIIcI 2 EMT/MET reporter while maintaining at least 50% viability or greater ( Figure 3E). These included GSK2126458 (mTOR/PI3K), three microtubule associated agents, TAK-875 (GPR40 agonist), PIK-75 (DNA-PK, p110α), Sparfloxacin (antibiotic), LY2228820 (p38/MAPK), AUY922 (HSP90), and Edoxaban (Factor Xa) (Figures 3F, G).

The chemical landscape of collateral sensitivity to enzalutamide-resistant prostate cancer
Given the association between Snail-mediated EMT and enzalutamide resistance, we hypothesized that the evolution of enzalutamide resistance may also enrich for this EMT-like plasticity. To better understand these relationships between phenotypic plasticity and enzalutamide resistance we applied a series of EMT scoring metrics (Chakraborty et al., 2021;Subbalakshmi et al., 2021;Pillai et al., 2022) to analyze RNA-Seq data from four independent pairs of enzalutamide-sensitive and enzalutamide-resistant cell line models . Consistent with our hypothesis, enzalutamide-resistant cells exhibited a significant shift in scores toward a more mesenchymal-like gene expression signature ( Figure 4A). These overall trends were consistent across scoring metrics, with some exceptions for specific cell line pairs, depending on the scoring metric used (Supplementary Figures S3A, B). Also consistent with this, treatment of LNCaP95-Snail(−) cells with enzalutamide led to an increase in nuclear localization of Snail ( Figures 4B, D). The enzalutamide-treated LNCaP95-Snail cells mirrored induction of Snail nuclear localization with 4OHT treatment (Figures 4C, D). These analyses indicate that, compared to enzalutamide-sensitive cells, enzalutamide-resistant cells exhibit a more EMT-like phenotype.
To further extend the analysis of Snail-specific synthetic lethality, we next attempted to identify potential collateral sensitivities to this EMT-like enzalutamide-resistant phenotype. In order to accomplish this we performed a separate highthroughput compound screen on enzalutamide-resistant

Frontiers in Molecular Biosciences
frontiersin.org CS2 cells. The CS2 model is an LNCaP-derived subclone that was generated from long-term exposure to androgen deprivation through chronic culture in media containing charcoal-stripped fetal bovine serum . Subsequent exposure of enzalutamide-sensitive CS2 cells to increasing doses of enzalutamide over approximately 6 months led to the development of an enzalutamide-resistant CS2 cell line model . The CS2 enzalutamide-resistant model was transduced with a lentiviral PSA reporter in which the proximal promoter of PSA harboring androgen responsive elements is inserted upstream of the GFP reading frame ( Figure 5A). These CS2 PSA−GFP enzalutamideresistant cells were screened using the Bioactives library to interrogate AR signaling (GFP), ATP production (CellTiter Glo), and cell growth (IncuCyte imaging) ( Figure 5A). To ensure the PSA reporter is responsive to androgen receptor signaling, cells were treated with the anabolic-androgenic steroid derivative, R1881, or enzalutamide. R1881 treatment led to a significant increase in GFP signal while enzalutamide had no effect on GFP expression in the enzalutamide-resistant CS2 model ( Figure 5B). The increase in GFP during R1881 treatment was not due to a change in confluence, as these treatments did not increase cell confluence ( Figure 5C). Analysis of cell growth inhibition for the Bioactives screen at the pathway level in the CS2 enzalutamide-resistant cells pinpointed candidate collateral sensitivities of interest, including DNA-PK, cyclin-dependent kinases, histone deacetylases, PI3K, mTOR, CRM1, and PLK ( Figure 5D). Analysis of PSA reporter expression as a function of cell viability also revealed compounds targeting multiple receptors (androgen receptor, estrogen receptor, glucocorticoid receptor) as inducers of PSA reporter activity ( Figure 5E) and compounds that target epigenetic modifiers as repressors of PSA reporter activity ( Figure 5F).
To provide further validation of candidates, we plotted the relative cell viability by CellTiter Glo for compounds in the CS2 enzaR screen by cell viability (CellTiter Glo) in the LNCaP95-Snail screen ( Figure 6A). This analysis revealed a subset of drugs active in both screens. We ranked these top hits by a sum rank statistic that includes the rank of cell death by CellTiter for both screens as well as the differential confluence for Snail− vs. Snail+ cells ( Figure 6B). Top targets from this analysis including PI3K, mTOR, and the proteasome ( Figure 6B). Among this subset, AZ 960 (JAK2 inhibitor) and BGT226 (dual PI3K/ mTOR inhibitor) were the most effective at inhibiting Snail+ cell confluence ( Figures 6C, D). Differences in confluence are not likely due to changes in cell proliferation, as each value was normalized to the average confluence of untreated cells. These results were further validated in both LNCaP95-Snail and CS2 EnzaR cells with AZ 960 treatment using IC50 dose response curves (Supplementary Figure S4). Consistent with our observations of sensitivity to JAK2 inhibition, analysis of phospho-proteomics data from three previously-characterized pairs of enzalutamide-resistant lines , including CS2 enzalutamide-sensitive and -resistant lines demonstrates increased phosphorylation of STAT1, STAT2, JAK1, and JAK2 ( Figures 6E-H), pinpointing the JAK/STAT signaling axis as a potential therapeutic vulnerability for Snail+ and enzalutamide-resistant prostate cancer.

Discussion
In the present study we sought to characterize the therapeutic vulnerabilities for enzalutamide-resistant prostate cancer. To do this we combined high-throughput small molecule screens with real- Frontiers in Molecular Biosciences frontiersin.org time imaging and endpoint assays to reveal chemical landscapes of synthetic lethality for Snail-mediated EMT and collateral sensitivities for enzalutamide-resistant cells. These screens identified multiple therapeutic vulnerabilities of Snail+ prostate cancer cells, including several with known functions in prostate cancer and/or EMT, such as aurora kinases (Beltran et al., 2011;Kivinummi et al., 2017;Beltran et al., 2019), MET (Chu et al., 2014;Lucas et al., 2014), polo-like kinases (Weichert et al., 2004;Liu et al., 2011;Deeraksa et al., 2013), and CRM1/XPO1 (Gravina et al., 2017;Wei et al., 2018). The screen also pinpointed several inhibitors that differentially inhibited EMT as indicated by a change in the EMT reporter output without a change in confluence, including inhibitors of mTOR/PI3K, DNA-PK, and p38/MAPK ( Figure 3H). All of these pathways have been previously connected to EMT biology in prostate cancer (Mulholland et al., 2012;Thakur et al., 2014;Zhang et al., 2019;Li et al., 2021). We also identified the GPR40 agonist, TAK-875, and Factor Xa inhibitor, Edoxoban, as potential inducers of MET. Consistent with these observations, another GPR40 agonist, GW9508, has been shown to prevent cytokine-induced airway epithelial barriers disruption of claudin, occludin, and ZO-1 (Moonwiriyakit et al., 2019), and Factor Xa inhibition has been shown to reduce EMT in chronic kidney disease (Fang et al., 2022). These agents represent promising candidates for follow-up studies to inhibit EMT and prevent or delay invasive and metastatic phenotypes associated with hormone therapy resistance. Similar to the screen for Snail+ prostate cancer the follow-up screen for therapeutic vulnerabilities in enzalutamide-resistant CS2 cells pinpointed targets and pathways known to be involved in prostate cancer and hormone therapy resistance, including histone deacetylases, the PI3K/mTOR pathway, JAK-STAT signaling, DNA-PK, and Syk. For example, the identification of histone deacetylases and other epigenetic modifying agents is consistent with the known importance of epigenetic regulation of androgen receptor signaling (Rokhlin et al., 2006;Welsbie et al., 2009). Other targets, however, are linked to AR signaling bypass, as in the case of PTEN loss and subsequent constitutive activation of PI3K signaling (Carver et al., 2011), activation of JAK/STAT and FGFR signaling during the acquisition of AR independence and lineage plasticity (Chan et al., 2022;Deng et al., 2022), and the role of Syk as a potential mediator of invasive features and bone metastasis (Ghotra et al., 2015). While the relevance of these targets is well supported by preclinical evidence, the clinical utility of these targets is more varied. For example, our identification of mTOR/PI3K signaling inhibition as a key vulnerability may be the result of PTEN loss in these LNCaP-derived models (Lotan et al., 2011); however, while PTEN loss is also common among patients, these agents have been unsuccessful in clinical trials (Cham et al., 2021). Likewise, currently-available HDAC inhibitors have largely failed in clinical trials, mostly due to their toxicity (Bradley et al., 2009) or lack of efficacy (Molife et al., 2010;Eigl et al., 2015). Conversely, we also identified compounds targeting microtubules and microtubule dynamics in both screens. This is consistent with the use of taxane chemotherapy in the hormone therapy-resistant setting (Tannock et al., 2004). In terms of novel agents and pathways, there are a number of ongoing clinical trials for JAK inhibitors-particularly JAK2 inhibitors-in advanced prostate cancer, but thus far these (B) Top hits for both screens based on a sum rank statistic that includes (CS2 enzaR confluence, Snail+ differential confluence, and Snail+ differential slope of growth rate). have not demonstrated sufficient monotherapy activity in men with mCRPC [ (Beinhoff et al., 2021) and see NCT00638378; closed due to lack of efficacy]. Our data suggest that a number of critical and nonredundant pathways may be involved in enzalutamide resistance and lineage plasticity, suggesting the need for combination trial approaches.
Comparison across both screens identified drugs with distinct effects in a single model as well as drugs that were common hits in both screens. There are multiple possible reasons for the observed differences in hits targeting each cell line, including, but not limited to, differences in the genetic and gene expression features of each cell line . For example, LNCaP95-Snail cells express AR-V7 while enzalutamide-resistant CS2 cells lack AR-V7. Enzalutamide-resistant CS2 cells also harbor dual loss of BRCA2 and RB1 and have a greater number of mutations and copy number alterations than LNCaP95 cells. These unique features may explain, at least in part, some of the differences in the list of hits from each screen.
A major limitation of the present study is the lack of in vivo modeling to validate the impact of our identified in vitro hits. This work is ongoing and also requires an assessment of the immune consequences of drug effects in the tumor microenvironment. Given the expression of mTOR, p38, and JAK/STAT signaling, for example, in immune cells and the immune suppressive impact of these agents in patients, assessing the net benefits of any drugs identified in our in vitro screen requires in vivo validation in a range of immunocompetent models either as monotherapy, in selected combinations and ideally in patient correlative samples.
The current study provides a platform to quantify the effects of thousands of compounds across multiple parameters and phenotypes simultaneously to identify and prioritize candidates for follow up in a rapid and cost-effective manner. While this study is limited by the exclusive use of in vitro cell line models, the integration of data from phenotypic reporters, confluence imaging, and CellTiter Glo readouts across multiple models rapidly identified a prioritized list of top hits, including the dual mTOR/PI3K inhibitor, BGT-226 and the JAK2 inhibitor, AZ-960, as promising candidates for future in vivo studies.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1104505/ full#supplementary-material SUPPLEMENTARY FIGURE S1 (A) Immunofluorescence staining of LNCaP95-Snail cells. EtOH serves as a vehicle for Snail induction. 4OHT induces localization of Snail and concomitant downregulation of E-cadherin and upregulation of vimentin.