L M Dillon1, J R Bean1, W Yang1, K Shee1, L K Symonds1, J M Balko2, W H McDonald3, S Liu4, A M Gonzalez-Angulo4, G B Mills5, C L Arteaga6, T W Miller7. 1. Department of Pharmacology and Toxicology, Norris Cotton Cancer Center, Geisel School of Medicine at Dartmouth, Dartmouth-Hitchcock Medical Center, Lebanon, NH, USA. 2. 1] Department of Medicine, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, TN, USA [2] Breast Cancer Research Program, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, TN, USA. 3. Proteomics Laboratory, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, TN, USA. 4. 1] Department of Breast Medical Oncology, MD Anderson Cancer Center, University of Texas, Houston, TX, USA [2] Department of Systems Biology, MD Anderson Cancer Center, University of Texas, Houston, TX, USA. 5. Department of Systems Biology, MD Anderson Cancer Center, University of Texas, Houston, TX, USA. 6. 1] Department of Medicine, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, TN, USA [2] Breast Cancer Research Program, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, TN, USA [3] Department of Cancer Biology, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, TN, USA. 7. 1] Department of Pharmacology and Toxicology, Norris Cotton Cancer Center, Geisel School of Medicine at Dartmouth, Dartmouth-Hitchcock Medical Center, Lebanon, NH, USA [2] Comprehensive Breast Cancer Program, Norris Cotton Cancer Center, Geisel School of Medicine at Dartmouth, Lebanon, NH, USA.
Abstract
Phosphatidylinositol 3-kinase (PI3K) promotes cancer cell survival, migration, growth and proliferation by generating phosphatidylinositol 3,4,5-trisphosphate (PIP3) in the inner leaflet of the plasma membrane. PIP3 recruits pleckstrin homology domain-containing proteins to the membrane to activate oncogenic signaling cascades. Anticancer therapeutics targeting the PI3K/AKT/mTOR (mammalian target of rapamycin) pathway are in clinical development. In a mass spectrometric screen to identify PIP3-regulated proteins in breast cancer cells, levels of the Rac activator PIP3-dependent Rac exchange factor-1 (P-REX1) increased in response to PI3K inhibition, and decreased upon loss of the PI3K antagonist phosphatase and tensin homolog (PTEN). P-REX1 mRNA and protein levels were positively correlated with ER expression, and inversely correlated with PI3K pathway activation in breast tumors as assessed by gene expression and phosphoproteomic analyses. P-REX1 increased activation of Rac1, PI3K/AKT and MEK/ERK signaling in a PTEN-independent manner, and promoted cell and tumor viability. Loss of P-REX1 or inhibition of Rac suppressed PI3K/AKT and MEK/ERK, and decreased viability. P-REX1 also promoted insulin-like growth factor-1 receptor activation, suggesting that P-REX1 provides positive feedback to activators upstream of PI3K. In support of a model where PIP3-driven P-REX1 promotes both PI3K/AKT and MEK/ERK signaling, high levels of P-REX1 mRNA (but not phospho-AKT or a transcriptomic signature of PI3K activation) were predictive of sensitivity to PI3K inhibitors among breast cancer cell lines. P-REX1 expression was highest in estrogen receptor-positive breast tumors compared with many other cancer subtypes, suggesting that neutralizing the P-REX1/Rac axis may provide a novel therapeutic approach to selectively abrogate oncogenic signaling in breast cancer cells.
Phosphatidylinositol 3-kinase (PI3K) promotes cancer cell survival, migration, growth and proliferation by generating phosphatidylinositol 3,4,5-trisphosphate (PIP3) in the inner leaflet of the plasma membrane. PIP3 recruits pleckstrin homology domain-containing proteins to the membrane to activate oncogenic signaling cascades. Anticancer therapeutics targeting the PI3K/AKT/mTOR (mammalian target of rapamycin) pathway are in clinical development. In a mass spectrometric screen to identify PIP3-regulated proteins in breast cancer cells, levels of the Rac activator PIP3-dependent Rac exchange factor-1 (P-REX1) increased in response to PI3K inhibition, and decreased upon loss of the PI3K antagonist phosphatase and tensin homolog (PTEN). P-REX1 mRNA and protein levels were positively correlated with ER expression, and inversely correlated with PI3K pathway activation in breast tumors as assessed by gene expression and phosphoproteomic analyses. P-REX1 increased activation of Rac1, PI3K/AKT and MEK/ERK signaling in a PTEN-independent manner, and promoted cell and tumor viability. Loss of P-REX1 or inhibition of Rac suppressed PI3K/AKT and MEK/ERK, and decreased viability. P-REX1 also promoted insulin-like growth factor-1 receptor activation, suggesting that P-REX1 provides positive feedback to activators upstream of PI3K. In support of a model where PIP3-driven P-REX1 promotes both PI3K/AKT and MEK/ERK signaling, high levels of P-REX1 mRNA (but not phospho-AKT or a transcriptomic signature of PI3K activation) were predictive of sensitivity to PI3K inhibitors among breast cancer cell lines. P-REX1 expression was highest in estrogen receptor-positive breast tumors compared with many other cancer subtypes, suggesting that neutralizing the P-REX1/Rac axis may provide a novel therapeutic approach to selectively abrogate oncogenic signaling in breast cancer cells.
The phosphatidylinositol 3-kinase (PI3K)/AKT/mTOR pathway is the most commonly aberrantly activated pathway in humancancer, promoting cell growth, migration, survival, and proliferation. Genes encoding proteins in this pathway are mutated in >70% of breast cancers, and pathway activation has been linked with resistance to anti-estrogen and anti-HER2 therapies in estrogen receptor α (ER)-positive and HER2-positive breast cancers, respectively [1]. PI3K pathway-targeted therapeutics are in clinical development. Class IA PI3Ks are activated by receptor tyrosine kinases (RTKs) such as members of the epidermal growth factor receptor family (EGFR, HER3) and the insulin receptor family (InsR, IGF-1R), and by Gβγ subunits released from activated G-protein coupled receptors (GPCRs) [2]. PI3K phosphorylates phosphatidylinositol (4,5)-bisphosphate (PIP2) in the inner leaflet of the plasma membrane to yield phosphatidylinositol (3,4,5)-trisphosphate (PIP3) [3]. Lipid phosphatase activity of the tumor suppressor PTEN antagonizes PI3K by dephosphorylating PIP3 [4]. PIP3 recruits AKT and other pleckstrin homology (PH) domain-containing proteins to the membrane to promote their interaction and initiate signaling cascades.In a proteomic screen to identify PIP3-regulated proteins as effectors and markers of PI3K activation in breast cancer cells, we found that levels of PIP3-dependent Rac exchange factor 1 (P-REX1) were inversely correlated with PI3K activation. This association was confirmed in humanbreast tumors. P-REX1 was most highly expressed in ER+ breast cancer compared with various other cancer subtypes, and P-REX1 mRNA levels were predictive of sensitivity to PI3K inhibition. We also found that P-REX1 activates IGF-1R/InsR, Rac1, PI3K/AKT, and MEK/ERK signaling, and promotes cell and tumor viability. These findings place P-REX1 both downstream and upstream of PI3K, suggesting that P-REX1 is a key component of a positive feedback loop that drives oncogenic signaling.
Results
PI3K activity regulates P-REX1 expression
We performed a mass spectrometry (MS)-based proteomic screen to identify proteins modulated by PI3K pathway activity in breast cancer cells. MCF-7 cells stably expressing shRNA targeting PTEN or mismatch control (shMM) were treated ± the pan-PI3K inhibitor BKM-120 [5] for 16 h. Lysates were analyzed by multi-dimensional protein identification technology (MudPIT) coupled with tandem MS/MS. Proteins exhibiting altered abundance in response to PI3K inhibition and/or PTEN knockdown are listed in Table S1. Comparing the lists of proteins altered ≥1.5 fold in abundance, we identified P-REX1 as a protein downregulated by PTEN knockdown (which activates PI3K), and upregulated by PI3K inhibition (Fig. 1A-B). P-REX1 is a guanine exchange factor (GEF) that integrates inputs from GPCRs via Gβγ, and from RTKs and GPCRs via PIP3 to activate the Rho family GTPases Rac1/2/3 [6]. Given that P-REX1 was previously implicated in prostate cancer metastasis [7] and shown to promote oncogenic signaling [8, 9], we further investigated its role in PI3K signaling in breast cancer.
Figure 1
P-REX1 levels decrease upon PTEN loss and increase with PI3K inhibition
A) Summary of results from mass spectrometry-based proteomic screen in MCF-7 cells (red asterisk indicates P-REX1). B) P-REX1 spectral counts from mass spectrometry shown for each of duplicate runs (run B for MCF-7/shMM + BKM-120 sample failed and was excluded from analysis). C-E) Immunoblots of lysates from cells treated with inhibitors were probed with the indicated antibodies.
We validated MS results by immunoblotting of lysates from breast cancer cells treated with inhibitors of PI3K or AKT. MCF-7/shMM and T47D/shMM cells treated with 1 μM BKM-120 showed reduced p-AKT and increased P-REX1 protein levels over time (Fig. 1C-D). The reason that P-REX1 levels did not increase as markedly in shPTEN cells compared to shMM cells upon PI3K inhibition is unclear but may be due in part to residual PIP3 in shPTEN cells, as evidenced by residual p-AKT in BKM-120-treated T47D/shPTEN cells (Fig. 1D). Since the P-REX1 homologue P-REX2a was found to bind PTEN [10], we considered that PTEN may bind and stabilize P-REX1. However, we detected only very low levels of P-REX1/PTEN interaction in cells (data not shown). Treatment with the pan-PI3K inhibitor GDC-0941 [11] similarly increased P-REX1 levels in parental MCF-7 and T47D cells, while AKT inhibition with MK-2206 did not (Fig. 1E). Thus, P-REX1 levels are modulated downstream of PI3K, but independently of AKT.We also explored whether PI3K activity modulates P-REX1 gene expression. PI3K inhibition for 4 or 24 h increased P-REX1 mRNA levels in MCF-7 cells (Fig. S1A-B). Mining of prior gene expression microarray data from isogenic MCF-7/shPTEN, T47D/shPTEN, and MDA-MB-361/shPTEN cells compared to respective shMM controls [12] revealed that P-REX1 mRNA levels were significantly downregulated in all shPTEN lines. In addition, PI3K inhibition increased P-REX1 protein levels in cells pretreated with the translation inhibitor cycloheximide (Fig. S1C). These data suggest that P-REX1 levels are modulated by PI3K signaling at the transcriptional and post-translational levels: PIP3-mediated activation of P-REX1 may induce P-REX1 protein turnover as a form of negative feedback, and PI3K inhibition blocks negative feedback to induce transcription of genes encoding proteins that activate and transduce PI3K signaling [e.g., ERBB3, IGF1R [13], TSC1 [14], and PREX1].
Tumor levels of P-REX1 are inversely correlated with PI3K/AKT/mTOR pathway activation, and positively correlated with ER expression in breast tumors
Since PI3K activity modulates P-REX1 levels in breast cancer cells, we explored this relationship in humantumors. Using reverse-phase protein array (RPPA) analysis of lysates from 712 primary breast tumors, we measured levels of P-REX1, phosphoprotein markers of PI3K pathway activation (p-AKTT308, p-AKTS473, p-4E-BP1T37/46, p-4E-BP1S65, p-p70S6KT389, p-mTORS2448, p-TSC2T1462, p-GSK3α/βS21/9), ER, and HER2. P-REX1 expression was inversely correlated with all PI3K pathway markers, positively correlated with ER levels, and modestly negatively correlated with HER2 levels (Fig. S2). Hierarchical clustering of P-REX1 and phosphoprotein data revealed 1) two (phospho)protein clusters, where P-REX1 clustered independently from PI3K pathway markers, and 2) several tumor clusters, where the two most distinct clusters exhibited high P-REX1 with low PI3K activation, and low P-REX1 with high PI3K activation (Fig. 2A).
Figure 2
P-REX1 levels are inversely correlated with PI3K activation, and positively correlated with ER expression in human tumors
A) Hierarchical clustering of signal values of P-REX1 and markers of PI3K pathway activation quantified by RPPA of 712 breast tumors. Levels of ER and HER2 are shown below, but were excluded from clustering. Red = high; blue = low. B) Correlation of P-REX1 mRNA levels and PI3K activation score in all tumor types (black) and breast tumors (red) from TCGA. Pearson correlation coefficients (r) and p-values are indicated. C) Tukey box plot comparing P-REX1 mRNA levels in breast tumors stratified based on expression of an mRNA signature of PTEN-loss vs. -intact (determined in ref. [17]) by Student’s t-test. D) Comparison of P-REX1 protein levels [from tumors in (A)] in ER+ vs. ER− breast tumors as in (C) (clinical ER status determined by IHC). E) Comparison of P-REX1 and ER mRNA levels in breast tumors from TCGA as in (B).
We then evaluated patterns of P-REX1 gene expression in humantumors. P-REX1 mRNA and protein levels were well-correlated (r2=0.6, p<0.0001)) in a set of 667 breast tumors (Fig. S3A). P-REX1 mRNA and protein levels were inversely correlated with a PI3K activation score derived from an expression signature of 1,793 genes (Table S2, derived from ref. [15]; signature does not include P-REX1) across 1,950 solid tumors and 919 breast tumors from The Cancer Genome Atlas (TCGA; Figs. 2B, S3B), and within two other datasets of 420 and 553 breast tumors [16] (Fig. S3C-D). In a fourth dataset of 100 breast tumors [17], P-REX1 mRNA levels were lower in tumors exhibiting a gene expression signature of PTEN loss compared to “PTEN-intact” tumors (Fig. 2C). These findings validate our observations in breast cancer cells (Figs. 1, S1), and confirm that P-REX1 levels are inversely correlated with PI3K activity.Stratifying breast tumors by clinical ER status or molecular subtype, we found that P-REX1 protein and mRNA levels were significantly higher in ER+/luminal tumors compared with ER−/basal tumors, and compared with a diversity of other tumor subtypes (Figs. 2D, S3E-J). P-REX1 and ER mRNA levels were strongly correlated in breast tumors (Fig. 2E). These observations suggest that P-REX1 may an ER+ breast tumor-specific protein.
Genetic lesions in PREX1 are enriched in cancers with mutations in the PI3K pathway
Co-existent mutations in genes encoding proteins that lie in the same signaling cascade (e.g., HER2 and PIK3CA [18]) are thought to provide robustness to promote oncogenic phenotypes and fitness [19]. We analyzed genomic datasets to determine whether PREX1 is genetically altered in humantumors, and whether PREX1 lesions co-exist with other PI3K pathway alterations. PREX1 is amplified or mutated in 3.65% (163/4,462) of cancers, and in 3.65% (25/685) of primary breast tumors (Table S3). PREX1 gene copy number significantly correlates with P-REX1 mRNA and protein levels in breast tumors (Fig. S3K-L), suggesting that PREX1 amplification confers increased expression. We then compared the genetic status of PREX1 and 79 PI3K pathway-related genes across 1,523 solid tumors (482 breast, 212 colorectal, 143 glioblastoma, 179 lung, 207 ovarian, 93 prostate, 207 sarcoma). PREX1 lesions significantly co-occurred with lesions in 51 of 79 PI3K pathway-related genes (Fisher’s exact test p<0.05; Table S4); however, PREX1 lesions were not significantly enriched in tumors with lesions in PIK3CA or PTEN. These results suggest that PREX1 alterations may be enriched in PI3K pathway-driven tumors to enhance PI3K pathway robustness, and/or vice versa.
P-REX1 activates IGF-1R/InsR, PI3K/AKT, and MEK/ERK signaling
Since PIP3 activates P-REX1 [6] and PI3K inhibition increases P-REX1 levels (Figs. 1, S1), we hypothesized that P-REX1 levels are regulated by negative feedback signaling from PI3K. To determine whether P-REX1 provides feedback to PI3K to form a complete circuit, we overexpressed exogenous myc-PREX1-HA or knocked-down endogenous P-REX1 using RNA interference in MCF-7 and T47D cells. P-REX1 overexpression increased AKT phosphorylation compared to control under IGF-1-stimulated and heregulin (HER3 ligand)-stimulated conditions (Fig. 3A-B). Conversely, P-REX1 knockdown decreased growth factor-induced and steady-state P-AKT in MCF-7 cells (Fig. 3A). These effects were confirmed using a second siRNA against the 3’ UTR of P-REX1, and restoration of P-REX1 expression using the myc-PREX1-HA cDNA construct (not targeted by siPREX1-2) rescued P-AKT levels (Fig. 3C).
Figure 3
P-REX1 activates IGF-1R/InsR, PI3K/AKT, and MEK/ERK signaling in breast cancer cells
Immunoblot analysis of lysates from cells transfected with vectors encoding myc-PREX1-HA, vector (vec) control, siRNA against P-REX1, or siControl. Panels A-C and E-G describe transiently transfected cells. Panel D describes stably transfected cells. Cells were maintained in 10% FBS (A,C-F), or serum-starved for 24 hours then stimulated ± 10-100 ng/ml IGF-1 for 10 min. (A,G), or 10 ng/ml heregulin (HRG) for 5 min. (B). In panel G, cells were pretreated ± 1 μM GDC-0941 for 2 h, then stimulated ± IGF-1.
Since P-REX2a directly inhibits PTENlipid phosphatase activity to increase PIP3 levels [10], we also evaluated P-REX1 effects in intrinsically PTEN-deficient breast cancer cells. P-REX1 overexpression in ZR75-1 and MDA-MB-415 cells increased p-AKT, while P-REX1 knockdown reduced p-AKT in ZR75-1 cells (Fig. 3D-F; note- ZR75-1 cells express high levels of endogenous P-REX1, so levels of myc-PREX1-HA overexpression are modest). Thus, the effects of P-REX1 on AKT activation are PTEN-independent. In contrast to reported effects of P-REX2a [10], we did not detect an effect of P-REX1 on PTENlipid phosphatase activity (data not shown).Since P-REX1 activates the PI3K/AKT pathway (Fig. 3A-F) that can crosstalk with Rac/Pak/Raf/MEK/ERK pathways [20], and P-REX1GEF function activates Rac1/2/3 [6, 9] that promote Pak/Raf/MEK/ERK activation [21-24], we tested whether P-REX1 also modulates MEK/ERK signaling in breast cancer. Transfected MCF-7 and T47D cells were serum-starved, then stimulated with IGF-1 or heregulin. P-REX1 overexpression increased and knockdown decreased basal and growth factor-induced p-MEK (Fig. 3A-B). Similar effects on p-MEK and p-ERK were observed under steady-state conditions in MDA-MB-415 and ZR75-1 cells (Fig. 3D-F).We observed that P-REX1 knockdown in MCF-7 and T47D cells reduced levels of p-IGF-1Rβ/p-InsRβ (Fig. 3A). To determine whether P-REX1 effects on IGF-1R/InsR activation were PI3K-dependent, we knocked-down P-REX1 expression in MCF-7 cells, pre-treated ± GDC-0941 (to inhibit PI3K), then stimulated ± IGF-1. While P-REX1 knockdown reduced p-IGF-1Rβ/p-InsRβ and increased total IGF-1Rβ, PI3K inhibition had no effect on receptor phosphorylation but decreased p-AKT (Fig. 3G). This suggests that P-REX1 promotes IGF-1Rβ/InsRβ activation independently of PI3K/AKT activity.
P-REX1 promotes PI3K/AKT and MEK/ERK signaling via Rac activation
Since P-REX1 activates Rac1/2/3 [6, 9], which can promote Pak/Raf/MEK/ERK activation [21-24], and Rac1 was recently found to activate the p110β isoform of Class IA PI3K [25, 26], we postulated that P-REX1 acts through Rac to drive activation of both PI3K/AKT and MEK/ERK signaling in ER+ breast cancer cells. We first assessed the role of P-REX1 in Rac1 activation in T47D cells with or without RNAi-induced P-REX1 knock-down. In an activated Rac1 pull-down assay, cell lysates treated with GDP (to inactivate Rac1) or GTPγS (non-hydrolyzable GTP to activate Rac1) were used as negative and positive controls, respectively. Immunoblot analysis of pull-down eluates showed that levels of activated Rac1 were lower in GDP-treated than GTPγS-treated lysates, confirming assay specificity for activated (GTP-bound) Rac1 (Fig. S4A). Surprisingly, synthetic Rac1 activation with GTPγS increased P-REX1/Rac1 interaction compared to synthetic Rac1 inactivation with GDP, suggesting that P-REX1 preferentially interacts with activated Rac1. P-REX1 knock-down decreased activated Rac1, which was found in complex with p-ERK. Rac1 has been found in a triple-complex with Pak1 and MEK1/2 [21], so the MEK1/2 substrates ERK1/2 could also be part of this complex. These findings confirm that P-REX1 activates Rac1 in ER+ breast cancer cells, and loss of P-REX1-induced Rac1 activation decreases p-ERK1/2.We then tested whether Rac signaling is required for activation of PI3K/AKT and MEK/ERK. Treatment with the Rac1/2/3 small molecule inhibitor EHT1864 [27] reduced steady-state, basal, IGF-1-induced, and heregulin-induced p-AKT and p-ERK1/2 in PTEN-deficient (ZR75-1, MDA-MB-415), PIK3CA-mutant (T47D, BT-474, MCF-7), and HER2+ (BT-474) breast cancer cells, all of which are ER+ (Figs. 4A-D, S4B). These results suggest that Rac GTPases drive activation of the PI3K/AKT and MEK/ERK pathways in ER+ breast cancer cells harboring genetic lesions in the PI3K pathway. Since P-REX1 activates Rac, these findings place P-REX1/Rac upstream of these signaling cascades.
Figure 4
Rac promotes activation of PI3K/AKT and MEK/ERK and is required for P-REX1-induced activation of these pathways in breast cancer
A-D) Immunoblot analysis of lysates from cells treated as indicated ± serum-starvation for 24 h, then pretreated ± 50 μM EHT1864 × 2 h, and stimulated ± 20 ng/mL HRG or 100 ng/mL IGF-1. E-H) Immunoblot analysis of lysates from cells transiently (E,G,H) or stably (F) transfected with myc-PREX1-HA or vector (vec) control, maintained in 10% FBS (E-F) or serum-starved for 24 hours (G-H), pretreated ± 50 μM EHT1864 × 2 h (E-H), then stimulated ± IGF-1 for 10 min. (G-H).
To confirm the involvement of Rac activation in P-REX1-induced oncogenic signaling, we tested the effects of Rac inhibition with EHT1864 on cells overexpressing myc-PREX1-HA or control. While P-REX1 overexpression increased p-MEK, p-ERK, and p-AKT, such effects were abrogated by treatment with EHT1864 (Fig. 4E-H). Interestingly, Rac inhibition increased levels of endogenous P-REX1 (Fig. 4E), suggesting that Rac signaling provides negative feedback to control P-REX1 levels.
P-REX1 promotes viability of breast cancer cells and tumors
Since P-REX1 increases oncogenic signaling, we tested whether P-REX1 also promotes cell viability. Cells transfected with siRNA against P-REX1, vectors encoding myc-PREX1-HA, or controls were assayed using XTT assay, which measures mitochondrial enzymatic activity. PREX1 overexpression in T47D, MCF-7, and MDA-MB-415 cells increased viable cell number, while PREX1 knockdown decreased viability in T47D cells (Fig. 5A-D). Similarly, inhibition of Rac with EHT1864 decreased viability in all lines tested (Fig. 5A-E). In agreement with prior findings [8], knockdown of P-REX1 also decreased MCF-7 xenograft growth in immunodeficientmice (Figs. 5F, S5).
Figure 5
P-REX1 increases viability of breast cancer cells
A-E) Cells transiently (A-B) or stably (C) transfected with vectors encoding myc-PREX1-HA or vector (vec) control, transiently transfected with siRNA (D), or untransfected (E), were evaluated by XTT assay. In (B-E), cells were treated ± 50 μM EHT1864 for 3 days prior to assay. Mean of 3-4 replicates ± SEM is shown. *p<0.05, **p<0.01 by Student’s t-test vs. control. F) Mice were subcutaneously implanted with MCF-7/shPREX1 or MCF-7/shNSC cells. Once tumor volumes reached ~150 mm3, tumors were measured twice weekly. Data are presented as % baseline volume (mean + SEM). Data were analyzed by two-way ANOVA followed by Bonferonni post-hoc comparison between cell lines at each time point (*p<0.05; **p<0.01).G) Forty-three breast cancer cell lines were dichotomized based on P-REX1 mRNA levels, and IC50 values were compared by Mann-Whitney U-test. Legend indicates molecular subtypes. Horizontal bars indicate median values.
P-REX1 contains amino-terminal tandem DH/PH domains characteristic of GEFs (Fig. S6A). The DH domain is required for binding Gβγ and GEF activity. The PH domain is needed for PIP3 binding and intracellular membrane localization [28]. Prior evidence indicates that PH domain binding to PIP3 may increase GTPase (i.e., Rac) binding activity of the DH domain in GEFs, although this remains controversial and may differ between GEFs (reviewed in ref. [29]). Adding to this understanding, we observed that deletion of the DH domain (ΔDH) increases P-REX1 binding to PIP3, and facilitates interaction with other phosphatidylinositol species (Figs. S6B-C), suggesting that the DH domain regulates interactions of the PH domain with lipid membranes. Adjacent to the DH/PH domains are two Dishevelled, Egl-10, and Pleckstrin (DEP) domains reported to be required for interaction with mTOR [30]; however, we did not detect an interaction of P-REX1 with mTOR in immunoprecipitation experiments (data not shown). Overexpression in T47D cells revealed that full-length P-REX1 increased viability, while the ΔDH and ΔDH/PH mutants were less effective (Fig. S6D). These findings suggest that P-REX1GEF activity and membrane localization are critical for P-REX1-mediated promotion of cell viability.Since P-REX1 requires PIP3 for activation [6], P-REX1 promotes PI3K/AKT and MEK/ERK signaling (Fig. 3), and combined inhibition of PI3K and MEK has strong anti-tumor effects (reviewed in ref. [20]), we hypothesized that P-REX1-expressing cells would be sensitized to pharmacologic inhibition of PI3K. Breast cancer cell lines (n=43) were dichotomized based on P-REX1 mRNA levels (extracted from the Cancer Cell Line Encyclopedia [31]), and IC50 values in response to GDC-0941 and the PI3K/mTOR dual inhibitor PI-103 (IC50 values extracted from [32]) were compared between P-REX1-high vs. -low cell lines. P-REX1-high cells were significantly more sensitive to GDC-0941 (p=0.01), and showed a trend towards sensitivity to PI-103 (p=0.07) compared to P-REX1-low cells (Fig. 5G). In agreement with P-REX1 levels being highest in ER+ breast cancers (Figs. 2D-E, S3), P-REX1-high cells classified mainly as luminal or HER2-enriched by molecular subtyping, while P-REX1-low cells were mostly basal or normal-like. In a multivariate model including molecular subtype, P-REX1 status (high/low) was not independently predictive of sensitivity to PI3K inhibition, suggesting that P-REX1 expression is tightly linked with molecular subtype. Transcriptomic PI3K activation score and P-AKTT308 levels were not predictive of sensitivity to PI3K inhibition (Fig. S7).
Discussion
Herein, we globally identified proteins quantitatively modulated by PI3K pathway inhibition and activation in breast cancer cells. We found that P-REX1 levels are inversely correlated with PI3K activation and positively correlated with ER expression in humanbreast tumors. P-REX1 activates Rac GTPases to promote activation of IGF-1R/InsR, PI3K/AKT, and MEK/ERK to increase breast cancer cell viability in a PI3K/PIP3-dependent, PTEN-independent manner. Therefore, P-REX1 lies both downstream and upstream of PI3K, and creates a positive feedback loop to drive oncogenic signaling (Fig. 6). In support of a model where P-REX1/Rac drives both PI3K/AKT and MEK/ERK signaling, breast cancer cells expressing high levels of P-REX1 are more sensitive to PI3K inhibition than P-REX1-low cells.
Figure 6
Model of PI3K-dependent P-REX1 positive feedback to activate Rac, IGF-1R/InsR, PI3K/AKT, and MEK/ERK signaling
Growth factor receptor tyrosine kinases (RTKs) and G-protein coupled receptors (GPCRs) activate PI3K, which produces PIP3. P-REX1 is recruited to the membrane via PH domain binding to PIP3. Activation of GPCRs induces dissociation of Gβγ subunits, which bind the P-REX1 DH domain to promote its GEF activity and Rac activation. Rac activates PI3K/p110β to increase PIP3 production, stimulates Raf/MEK signaling via Pak, interacts with p-ERK, and promotes actin cytoskeleton remodeling. P-REX1 also promotes PI3K-independent activation of IGF-1R/InsR, which can lead to further activation of PI3K and MEK signaling.
We and others found that P-REX1 knockdown reduces levels of activated Rac in breast cancer cells (Fig. S4A and refs. [9, 33]). In addition, we demonstrated that P-REX1 activates IGF-1R/InsR, PI3K/AKT, and MEK/ERK signaling in breast cancer cells in a Rac activation-dependent manner (Figs. 3-5, S4). Recent studies have implicated Ras and Rho family GTPases in the regulation of PI3K activity. Fritsch et al. showed that activated Rac1 binds and activates PI3K/p110β to increase PIP3 production [25, 26]. Therefore, P-REX1 may activate PI3K/p110β via Rac in breast cancer cells. P-REX1 also promotes activation of IGF-1R/InsR in a PI3K/PIP3-independent, Rac-dependent manner (Figs. 3G, 4G-H), which in turn can activate PI3K/p110α [34]. We did not consistently detect P-REX1/IGF-1Rβ or P-REX1/IRS-1 interaction in immunoprecipitation trials (data not shown), and the mechanism underlying P-REX1-mediated receptor activation requires further study.Cross-talk between the PI3K/AKT and Pak/Raf/MEK/ERK pathways exists at several levels [20], and we demonstrated that P-REX1 modulates both pathways in ER+ breast cancer cells. P-REX1 activates MEK/ERK signaling via activation of Rac (Figs. 4, S4). The Rac effector Pak1 can serve as an intermediary to activate Raf/MEK/ERK signaling (Fig. 6 & refs. [21-24]). In support of our findings, a recent study showed that P-REX1 knock-down reduces Rac1/Pak/Raf/MEK/ERK signaling [33]. Furthermore, Rac/Pak1-mediated activation of ERK and AKT is required for HER2-induced transformation of breast cancer cells [35]. We thus infer that P-REX1GEF activity promotes Rac/Pak-induced activation of MEK/ERK and PI3K/AKT in ER+ breast cancer cells. Additionally, we detected binding of activated Rac1 to p-ERK1/2, and synthetic Rac1 activation with GTPγS increased P-REX1/Rac1 interaction (Fig. S4A). A Rac1/Pak1/MEK1/2 complex has been described [21]; we speculate that this complex may also bind ERK1/2. Whether P-REX1 interacts with Rac1 when complexed with Pak1/MEK1/2/ERK1/2 is unclear. Given that GEFs typically interact with GDP-bound GTPases, GTPase-activating proteins (GAPs) interact with GTP-bound GTPases, and other GEF/GAP/Rho GTPase complexes have been described [36], it is conceivable that P-REX1 interacts with a complex containing activated Rac1 and a GAP.In support of a role for P-REX1 driving both the PI3K/AKT and MEK/ERK pathways, we observed that high P-REX1 mRNA levels, but not PI3K activation score (by gene expression) or P-AKTT308 level, are predictive of increased sensitivity to PI3K inhibitors (Figs. 5G, S7). Drug combinations targeting PI3K and MEK elicit dramatic anti-tumor effects in preclinical studies, but such combinations are accompanied by considerable toxicity in patients [20]. Thus, P-REX1 may be a marker of cells in which PI3K inhibitors alone will block both PI3K and MEK. Indeed, PI3K inhibition decreases P-ERK levels which, in turn, induces upregulation of the pro-apoptotic protein Bim in P-REX1-expressing cells [33]. Interestingly, the two luminal/ER+ breast cancer cell lines most resistant to GDC-0941 (HCC-1500 and MDA-MB-134VI; Fig. 5G) harbor PAK1 amplification, and Pak1 knock-down or treatment with the Pak inhibitor PF-3758309 was shown to induce apoptosis and decrease MEK/ERK phosphorylation in these cells [37]. PAK1 is amplified in 8.4%, 13.8%, and 3.7% of luminal, HER2-enriched, and basal breast cancers, respectively (extracted from ref. [38] using cBio [39]). Thus, P-REX1 expression and PAK1 amplification may serve as respective positive and negative biomarkers of sensitivity to PI3K inhibitors in breast cancer.The P-REX1 homolog P-REX2a binds PTEN and inhibits its lipid phosphatase activity [10]. While we detected a low level of P-REX1/PTEN binding in situ, in vitro and in situ analyses did not reveal an effect of P-REX1 on PTEN phosphatase activity (data not shown). Furthermore, we observed that P-REX1 activates PI3K/AKT and MEK/ERK signaling in PTEN-null cells (MDA-MB-415 and ZR75-1; Figs. 3D-F, 4F). Thus, P-REX1 and P-REX2a may promote PI3K pathway activation via different mechanisms.In summary, we report that P-REX1 mediates a PI3K/PIP3-dependent positive feedback loop to promote activation of Rac, IGF-1R/InsR, PI3K/AKT, and MEK/ERK signaling in ER+ breast cancer. Combined targeting of two key signaling nodes in cancer, PI3K and MEK, is being investigated in several ongoing clinical studies [20]. Thus, targeting the P-REX1/Rac axis may inhibit both PI3K and MEK signaling to elicit anti-cancer effects. Our mechanistic findings on P-REX1 are from ER+ breast cancer models harboring genetic lesions in the PI3K pathway; P-REX1 effects in ERBB2/PIK3CA/PTEN-wild-type ER+ models, and in ER− models remain to be determined. Furthermore, P-REX1 mRNA expression is predictive of sensitivity to PI3K inhibitors in breast cancer cells, offering P-REX1 as a biomarker to identify tumors likely to respond to these agents. P-REX1 is most highly expressed in ER+ breast cancers, suggesting that P-REX1 is an oncoprotein mainly in this subtype of cancer. Since P-REX1 knock-out mice are generally healthy with mild neutrophilia [40], targeting of P-REX1 may provide a high therapeutic index in patients.
Material and Methods
Cell culture
Cell lines (ATCC) were grown in DMEM (Life Technologies) with 10% FBS (Hyclone). MCF-7 and T47D cells stably expressing shRNA targeting PTEN or mismatch control (shMM) were described previously [41]. Generation of other stable cell lines is described in Supplementary Methods. In transient transfection experiments, cells were transfected with pcDNA plasmid encoding myc-PREX1-HA (from Atanasio Pandiella, CSIC-Universidad de Salamanca) or FLAG-HA control (Addgene) plus Lipofectamine 2000. Cells were transfected with siRNA targeting P-REX1 or AllStars Negative Control (Qiagen, Dharmacon) plus Lipofectamine RNAiMAX (Life Technologies). Cells were treated with IGF-1 (PeproTech), heregulin (HRG), BKM-120, GDC-0941, MK-2206 (Selleckchem), or EHT1864 (Tocris Bioscience) as indicated in figures.
Mass spectrometry (MS)
Lysates from MCF-7/shPTEN and MCF-7/shMM cells treated with 1 μM BKM-120 or DMSO (control) for 16 h were analyzed by MudPIT coupled with tandem MS/MS [42, 43]. Proteins were identified using the Sequest algorithm [44] and assembled into groups using IDPicker [45] (details in Supplementary Methods).
Immunoblotting
Cells were lysed in RIPA buffer containing protease and phosphatase inhibitors. Lysates were sonicated for 15 sec., centrifuged at 18k × g for 10 min. at 4°C, and protein in supernatants was quantified using BCA assay (Pierce). Reduced/denatured lysates were separated by SDS-PAGE, and proteins were transferred to nitrocellulose. Even protein loading across lanes was confirmed with Ponceau S staining. Blots were probed with antibodies against P-REX1 (Abcam), HA, HER3 (Santa Cruz), Rac1 (Cytoskeleton), p-AKTS473, p-AKTT308, AKT, p-S6S240/4, p-ERK1/2T202/Y204, p-MEK1/2S217/221, p-IGF-1RβY1135/6/p-InsRβY1150/1, IGF-1R, IRS-1, p-HER3Y1289, or Actin (Cell Signaling). HRP-labeled secondary antibodies (GE Healthcare) and ECL substrate (Pierce) were used for signal detection.
Cell viability assay
XTT cell viability assay (Biotium) was performed per manufacturer’s instructions. Cells were transfected with siRNA or plasmid, then reseeded 48 h later at 1-1.5 × 104/well in 96-well plates. Assay was performed 48-72 h later by adding 25 μL activated XTT solution, and reading absorbance at 450 nm (reflective of live cells) and 655 nm (background) after 2-6 hours. Background-subtracted values (Abs450 - Abs655) of 3-4 replicate wells were compared by Student’s t-test. Similarly, cells were seeded, then treated the next day ± 50 μM EHT1864, and XTT assay was performed 48 h later.
Breast cancer xenografts
Animal studies were approved by the Dartmouth College IACUC. Female NOD-scid IL2Rγ−/− (NSG; NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ) mice (6-7 wks old; obtained from the Norris Cotton Cancer Center Transgenic & Genetic Construct Shared Resource) were implanted subcutaneously with a 60-day-release 17β-estradiol pellet (0.72 mg; Innovative Research of America) and 5 × 106 cells (MCF-7/shPREX1 or /shNSC; described in Supplementary Methods) suspended in 50% matrigel (BD Biosciences). Tumor diameters were measured twice per week using calipers (volume = length2 × width/2), and volumes were analyzed as percent baseline by two-way ANOVA followed by Bonferroni post-hoc comparison.
RPPA analysis
Tumors analyzed by RPPA were part of TCGA project. P-REX1 was quantified specifically for this study. Other (phospho)proteins were quantified as part of TCGA [46]. All tumor proteomic data has been deposited at The Cancer Proteome Atlas (TCPA; http://app1.bioinformatics.mdanderson.org/tcpa/_design/basic/index.html). Tumor lysates were analyzed as described in Supplementary Methods.
Comparing P-REX1 gene expression with sensitivity to PI3K inhibition
RMA-normalized P-REX1 mRNA expression values from cancer cell lines were obtained from the Cancer Cell Line Encyclopedia [31]. Molecular subtypes of 43 breast cancer cell lines were determined by PAM50 subtyping using R software with package ‘genefu,’ median-centered data, and the 306 intrinsic subtype genes and centroids identified in ref. [47]. Cell lines were dichotomized based on high vs. low P-REX1 mRNA levels, and IC50 values in response to treatment with GDC-0941 and PI-103 (obtained from ref. [32]) were compared between groups by Mann-Whitney U-test.
Authors: Heidi C E Welch; W John Coadwell; Christian D Ellson; G John Ferguson; Simon R Andrews; Hediye Erdjument-Bromage; Paul Tempst; Phillip T Hawkins; Len R Stephens Journal: Cell Date: 2002-03-22 Impact factor: 41.582
Authors: Michael J MacCoss; W Hayes McDonald; Anita Saraf; Rovshan Sadygov; Judy M Clark; Joseph J Tasto; Kathleen L Gould; Dirk Wolters; Michael Washburn; Avery Weiss; John I Clark; John R Yates Journal: Proc Natl Acad Sci U S A Date: 2002-06-11 Impact factor: 11.205
Authors: Maria Soledad Sosa; Cynthia Lopez-Haber; Chengfeng Yang; Hongbin Wang; Mark A Lemmon; John M Busillo; Jiansong Luo; Jeffrey L Benovic; Andres Klein-Szanto; Hiroshi Yagi; J Silvio Gutkind; Ramon E Parsons; Marcelo G Kazanietz Journal: Mol Cell Date: 2010-12-22 Impact factor: 17.970
Authors: Lao H Saal; Peter Johansson; Karolina Holm; Sofia K Gruvberger-Saal; Qing-Bai She; Matthew Maurer; Susan Koujak; Adolfo A Ferrando; Per Malmström; Lorenzo Memeo; Jorma Isola; Pär-Ola Bendahl; Neal Rosen; Hanina Hibshoosh; Markus Ringnér; Ake Borg; Ramon Parsons Journal: Proc Natl Acad Sci U S A Date: 2007-04-23 Impact factor: 11.205
Authors: Ariella B Hanker; Adam D Pfefferle; Justin M Balko; María Gabriela Kuba; Christian D Young; Violeta Sánchez; Cammie R Sutton; Hailing Cheng; Charles M Perou; Jean J Zhao; Rebecca S Cook; Carlos L Arteaga Journal: Proc Natl Acad Sci U S A Date: 2013-08-12 Impact factor: 11.205
Authors: Heng-Jia Liu; Lisa M Ooms; Nuthasuda Srijakotre; Joey Man; Jessica Vieusseux; JoAnne E Waters; Yue Feng; Charles G Bailey; John E J Rasko; John T Price; Christina A Mitchell Journal: J Biol Chem Date: 2016-06-29 Impact factor: 5.157
Authors: Meagan B Ryan; Alexander J Finn; Katherine H Pedone; Nancy E Thomas; Channing J Der; Adrienne D Cox Journal: Mol Cancer Res Date: 2016-07-14 Impact factor: 5.852
Authors: Wei Yang; Sarah R Hosford; Nicole A Traphagen; Kevin Shee; Eugene Demidenko; Stephanie Liu; Todd W Miller Journal: FASEB J Date: 2018-01-03 Impact factor: 5.191
Authors: Wei Yang; Sarah R Hosford; Lloye M Dillon; Kevin Shee; Stephanie C Liu; Jennifer R Bean; Laurent Salphati; Jodie Pang; Xiaolin Zhang; Michelle A Nannini; Eugene Demidenko; Darcy Bates; Lionel D Lewis; Jonathan D Marotti; Alan R Eastman; Todd W Miller Journal: Clin Cancer Res Date: 2016-01-05 Impact factor: 12.531
Authors: Nuthasuda Srijakotre; Heng-Jia Liu; Max Nobis; Joey Man; Hon Yan Kelvin Yip; Antonella Papa; Helen E Abud; Kurt I Anderson; Heidi C E Welch; Tony Tiganis; Paul Timpson; Catriona A McLean; Lisa M Ooms; Christina A Mitchell Journal: Proc Natl Acad Sci U S A Date: 2020-10-23 Impact factor: 11.205