Quantitative analysis of castration resistant prostate cancer progression through phosphoproteome signaling
© Lescarbeau and Kaplan; licensee BioMed Central Ltd. 2014
Received: 16 September 2013
Accepted: 21 April 2014
Published: 8 May 2014
Although recent progress has been made in treating castration resistant prostate cancer, the interplay of signaling pathways which enable castration resistant growth is incompletely understood. A data driven, multivariate approach, was used in this study to predict prostate cancer cell survival based on the phosphorylation levels of key proteins in PC3, LNCaP, and MDA-PCa-2b cell lines in response to EGF, IGF1, IL6, TNFα, dihydrotestosterone, and docetaxel treatment.
The prostate cancer cell lines were treated with ligands or inhibitors, cell lyates were collected, and the amount of phosphoprotein quantified using 384 well ELISA assays. In separate experiments, relative cell viability was determined using an MTT assay. Normalized data was imported into Matlab where regression analysis was performed.
Based on a linear model developed using partial least squares regression, p-Erk1/2 was found to correlate with castration resistant survival along with p-RPS6, and this model was determined to have a leave-one-out cross validated R2 value of 0.61. The effect of androgen on the phosphoproteome was examined, and increases in PI3K related phosphoproteins (p-Akt, p-RPS6, and p-GSK3) were observed which accounted for the majority of the significant increase in androgen-mediated cell survival. Simultaneous inhibition of the PI3K pathway and treatment with androgen resulted in a non-significant increase in survival. Given the strong effect of PI3K related signaling in enabling castration resistant survival, the specific effect of mTor versus complete inhibition was examined using targeted inhibitors. It was determine that mTor inhibition accounts for 52% of the effect of complete PI3K inhibition on cell survival. The differences in signaling between the cell lines were explored it was observed that MDA-PCa-2b exhibited far less activation of p-Erk in response to varying treatments, explaining one of the reasons for the lack of castration resistance.
In this work, regression analysis to the phosphoproteome was used to illustrate the sources of castration resistance between the cell lines including reduced p-Erk signaling in MDA-PCa-2b and variations in p-JNK across the cell lines, as well as studying the signaling pathways which androgen acts through, and determining the response to treatment with targeted inhibitors.
KeywordsProstate cancer Phosphoproteome Castration resistance Regression analysis Cell signaling
Every year 223,000 men will be diagnosed with prostate cancer in the United States with most patients having androgen dependent disease at the initial stages . Although there have been recent advances in treating castration resistant prostate cancer, prognoses are still poor once the disease progresses to the castration resistant, metastatic state [2, 3]. There have been numerous mechanisms reported which can enable castration resistant growth including intracrine synthesis of androgen, upregulation of the androgen receptor (AR), co-activation of the AR by other pathways, or complete bypass of androgen signaling through the activation of other pathways [4–6]. These mechanisms can include activation of oncogenes, mutation of tumor suppressors, epigenetic alterations, or activation of a pathway through extracellular matrix or ligand cues contained in the microenvironment.
The signaling mechanisms which enable castration resistant growth have been studied using various cell line models, including PC3, LNCaP, and MDA-PCa-2b cells lines. These cell lines display a range of phenotypes, including aggressive castration resistant growth in PC3 cells and androgen-dependent growth in LNCaP and MDA-PCa-2b cells. These cell lines additionally display various mutations in their genome with LNCaP and PC3 cells having inactivated PTEN (Phosphatase and tensin homolog) and MDA-PCa-2b cells having intact PTEN [7, 8].These differences are used to model the variation present in patients with differing stages of disease progression. Depending on the cell line, certain growth factor treatments such as EGF (Epidermal growth factor) or IGF1 (Insulin-like growth factor 1), or targeted kinase inhibitor treatments, can enhance castration resistant growth or treat castration resistant cancer through modulating signal transduction pathways.
The analysis of prostate cancer signaling often involves the examination of numerous pathways through genomics, transcriptomics, or proteomics. The relationship of these data sets to cell phenotype is often multivariate and non-intuitive. To investigate these relationships, multivariate linear regression techniques have been utilized over the last decade, and have been successful in correlating the signaling of multiple pathways using phosphoproteomic data to phenotypic outcomes including apoptosis, proliferation, invasion, and migration [9–11]. Partial least squares (PLS) regression is a multiple linear regression algorithm which correlates variation in the Y matrix (cell survival) to the X matrix (phosphoprotein levels) by identifying vectors which simultaneously describe variation in both data sets. These latent variables are able to account for the multicollinearity of similarly regulated phosphoproteins (i.e. phosphosites which may be part of the same pathway or crosstalk between pathways).
In the present work, the objective was to correlate castration resistant growth to pathway activation via phosphoproteomic signaling using regression analysis. The use of the PC3, LNCaP, and MDA-PCa-2b cell lines allowed us to capture diversity in different prostate cancer genotypes, and make comparisons across cell lines. The epigenetic and genetic variations are assumed to be abstracted into the levels of phosphoprotein activation (i.e., PTEN inactivating mutations causing higher levels of p-Akt) with differences in unmeasured pathways across cell lines being sources of error in the model. This approach enables the exploration of a range of hypotheses to understand how cell signaling drives castration resistance, the importance of various signaling proteins in enabling castration resistant growth, the correlation between these signaling proteins, and the specific effect of various targeted kinase inhibitors in modulating the effect of these signaling proteins. This work will aid in the long term goal of optimizing the inhibition of signaling pathways to prevent castration resistant prostate cancer progression.
Cell culture and reagents
LNCaP, MDA-PCa-2b, and PC3 cell lines were acquired from ATCC (Manassas, VA, USA). PC3 and LNCaP cells lines were cultured in 10% fetal bovine serum (FBS), RPMI 1640, and 1% antibiotic-antimycotic. The MDA-PCa-2b cell line was cultured in BRFF-HPC1 media purchased from AthenaES (Baltimore, MD, USA) supplemented with 20% FBS. Dihydrotestosterone was acquired from Sigma-Aldrich (St. Louis, MO, USA). Androgen depleted media consisted of 10% charcoal stripped FBS with phenol red free RPMI 1640. Docetaxel was acquired from Sigma-Aldrich. Temsirolimus and SB202190 were purchased from Selleckchem (Houston, TX, USA). All other inhibitors were purchased from EMD Millipore (Billerica, MA, USA). Unless otherwise stated all other cell culture reagents were acquired from Invitrogen (Grand Island, NY, USA).
Cell survival assay
Relative cell viability was assessed using an MTT ((3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay acquired from Invitrogen. As previously determined by our lab, MTT correlates to relative cell number as confirmed via DNA quantification and manual cell counting . All three cell lines were plated at a concentration of 5,000 cells/cm2 in a 24 well plate in their respective growth media. The cells were allowed to adhere for 24 hours. The media was then changed to androgen depleted media which the cells were cultured in for an additional 72 hours. Finally, relative cell viability was determined using an MTT assay according to the manufacturer’s instructions. Targeted kinase inhibitors were used at the following concentrations: LY294002 at 7 μm, U0126 at 325 nm, Wedelolactone at 10 μm, Temsirolimus at 50 nm, and SB202190 at 500 nm. Additionally, the total protein amount of biological replicates from each cell type was measured using a Bicinchoninic assay purchased from Thermo Scientific (Rockford, IL, USA). After measuring cell survival with an MTT assay the results were normalized to total protein measured to account for variations in cell size between the cell lines.
Measuring phosphoprotein levels
Each prostate cancer cell line was plated to six well plates at a density of 7,500 cells/cm2 in their respective growth media and allowed to adhere for 24 hours. After 24 hours cells were treated with androgen depleted media supplemented with the appropriate treatment. For studies involving the use of inhibitors on LNCaP cells, the cells were first pretreated for 30 minutes with the inhibitor before additional treatments were added to ensure complete inhibition. Following the appropriate amount of time (30 minutes, 4 hours, or 24 hours) the media was removed and the cells were lysed. R&D Systems (Minneapolis, MN, USA) Duoset ELISA kits were used to quantify the amount of phosphoprotein present in each sample. Lysates were processed and the assays performed according to manufacturer’s instructions. A Bicinchoninic acid assay was performed on each lysate and the lysates were diluted such that 20 ug of protein lysate was used in each ELISA assay. The lysis buffer was made by combining 20 ml of PBS, 1nM of EDTA, 5 mM NaF, 6 M Urea, 0.1 ml of Triton ×100, and 2 packets on Halt protease/phosphatase inhibitors from Thermo Scientific (Waltham, MA, USA). Briefly, the antibody pairs for the ELISA assays were optimized on 384 well ELISA plates from Santa Cruz Biotechnology (Dallas, Texas, USA) using the accompanied positive control samples. An eight point standard curve was generated and fitted using a second order polynomial. The amount of phosphoprotein in ng per 20 ug of total protein lysate was then determined by comparing the measured absorbance of the sample to the standard curve.
Following data acquisition, calibration to the ELISA standard curve, and normalization to total protein content, the data was imported into Matlab (The Mathworks, Natick, MA) where both protein (X matrix) and survival data (Y matrix) were mean centered and unit variance scaled. The data was arranged such that each column of the X matrix represented a phosphoprotein at a specific time (8 phosphoproteins × 3 time points = 24 columns). The rows represent the cell treatments with the values in the X matrix corresponding to phosphorylation levels and the rows of the Y matrix corresponding to relative cell survival in response to that treatment. The X and Y matrices were then inputted into a function which utilizes the native plsregress function packaged with Matlab to employ the SIMPLS algorithm and calculate the regression coefficients. This was repeated with each row (treatment condition) left out. The calculated model was applied to the left out data to determine a predicted Y value. The R2 value was then calculated using the measured (Y vector) and predicted survival data. Partial least squares regression is a multiple regression algorithm which attempts to explain the Y matrix by finding a multidimensional direction in the X space which explains the maximum variation in both matrices . This algorithm is especially suited to applications where the X matrix contains many more variables than observations, or when many of the X variables are multicollinear, as is often the case in cell signaling data.
An approach for calculating significance in PLS regression models was employed which randomizes the X matrix (phosphoprotein data) as compared to the Y matrix (survival data) and performs regression analysis. From this randomized regression a R2 is calculated and saved. We repeated this procedure 3,000 times and determined a mean R2 and standard deviation for these calculated random models. The randomized R2 values were assumed to follow a normal distribution. Using the mean and standard deviation from the R2 values calculated for randomized regression, and the R2 of the correctly calculated model, the number of standard deviations away from the random mean was determined (z-score), and from this a p-value determined.
The level of phosphoprotein activation in response to ligand treatment was calculated as a percent increase over untreated controls. This data was imported into Cytoscape and used as relative measures of edge thickness between ligand and the resulting phosphoproteins . Decreases in phosphoprotein levels in response to treatment were depicted as a red edge.
To model the correlation between the phosphosites in the three different cell lines, the Pearson correlation between all possible unique pairs of phosphosites within the same cell line were assessed and a P-value calculated which represents the statistical significance of the correlation. This was completed on observations for untreated cells, EGF, IGF1, IL6, TNFα, DHT, and docetaxel treated cells using all three time points (30 minutes, 4 hours, and 24 hours). The Q-value (P-value equivalent adjusted for multiple hypothesis testing) was also determined using the Q-value software downloaded from the Storey lab website to adjust for multiple hypothesis testing .
Measuring phosphorylation and castration resistant survival in response to treatment
Erk1: T202/Y204 
RPS6: S235/S236 
GSK3α: S21 
p38δ: T180/Y182 
JNK1 and JNK2: T183/Y185 
HSP27: S78/S82 
Stat3: Y705 
After the phosphoprotein data was collected and normalized (see methods), hierarchical clustering analysis was applied across the phosphosites at the three time points as well as the treatment groups. This analysis measures the similarity between each observation using a Euclidean distance metric (Figure 1B). Across the y dimension of the X matrix, the treatments were found to cluster first by cell line and then by inhibitor treatment (for LNCaP cells only), with little clustering in the ligand treatment groups (Figure 1B and Additional file 2: Table S1). In the x dimension the phosphoprotein activation was generally found to cluster the three time points of each phosphoprotein together (Figure 1B and Additional file 3: Table S2). This clustering indicated that the cell line, and then inhibitor, and finally the ligand treatment imparted the most substantial changes in the cells in the y dimension (treatment conditions). In the × dimension (phosphoproteins), the data indicated that the change by time point tended to cause the most substantial response in phosphoprotein levels.
For each treatment, biological duplicates were measured and the absolute percentage difference between the two replicates was determined (Additional file 4: Table S3). A mean difference of 20.4% was observed across all cell lines which when compared to the finding that the phosphosites varied by approximately 670% on average over untreated controls, was considered an acceptable amount of error.
Regression analysis correlating phosphoprotein measurements to cell survival in androgen depleted conditions
Upon determining that this model was significantly more accurate than a randomized model, we examined the regression coefficients to determine weights calculated on the different phosphoproteins. Consistently positive coefficients for p-Erk (Extracellular signal-regulated kinases) were noted, as well as consistently increased p-RPS6 (Ribosomal protein S6) across all time points (Figure 2B). p-JNK regression coefficients were negative at all time points along with p-Akt and p-Stat3. p-GSK3 (Glycogen synthase kinase 3) additionally had minimal early and late time point regression coefficients, however had a substantially increased 4 hour regression coefficient.
In order to better assess the contribution of the regression coefficients to the model outcome the absolute value of the coefficients was taken for each time point and the mean plotted for each phosphoprotein in descending order (Figure 2C). From this, p-Erk was determined to most strongly contribute to the model, followed by p-RPS6 and p-JNK. We used this data to plot the R2 value of models built on increasing amounts of data, starting with p-Erk and adding phosphoproteins in order of their mean absolute value of regression coefficients. It can be seen that a model built solely on p-Erk, p-RPS6, and p-JNK resulted in R2 values of 0.4655 as compared to the complete model which gave us a R2 value of 0.616 (Figure 2D). Beyond these phosphoproteins, only the Akt phosphoprotein added substantial further information to the model, increasing the R2 from 0.484 to 0.570, indicating this data added substantial accuracy to the model without having a large regression coefficient. From these results it was concluded that the phosphorylation levels of Erk, RPS6, JNK, and Akt were able to explain the majority of variation in castration resistant survival across these three cell lines.
The amount of error between the predicted values from the model and the measured values were also grouped by treatment, cell line, and inhibitor (for LNCaP cells treated in combination with targeted inhibitors) (Additional file 5: Figure S1A, B, and C). The only significant difference that was observed between any conditions was a much higher docetaxel error (Additional file 5: Figure S1A). This is likely due to the fact that docetaxel is a chemotherapeutic which causes cell death, however little variation in the phosphoproteome as compared to controls was seen. Therefore a model of phosphoproteomic signaling was unable to predict docetaxel’s apoptotic effect.
The effect of androgen treatment on phosphoprotein signaling
The survival of LNCaP cells in response to DHT treatment was examined and an increase of 38% was observed as compared to the control condition (Figure 3B). This survival advantage was completely abrogated when treated in combination with LY294002 (PI3Ki) which reduced p-Akt, p-GSk3, and p-RPS6 to below baseline levels at all time points. The combination of DHT plus LY294002 caused a non-significant increase in survival of 25% over the treatment of LY294002. There was little difference in phosphoprotein levels from LY294002 treatment alone, indicating direct activation of the cell cycle by AR or activation of other non-measured pathways by AR other than PI3K.
Based on these observations we propose a modification of the model originally proposed by Gosh et al. (Figure 3C) . Here, the PI3K pathway can activate the AR which can activate the cell cycle. However, activation of the AR can also activate the PI3K pathway. Additionally, activation of the PI3K pathway can activate cell cycle through bypassing the AR via mTOR/RPS6.
Comparison of phosphoprotein alterations between LNCaP, MDA-PCa-2b, and PC3 cell lines
A final observation made regarding the mean phosphoprotein levels across all treatments was the decreasing levels of phosphorylation in JNK from MDA-PCa-2b cells to LNCaPs and then PC3 cells (Figure 4B, C, and D). Initially, this was a counterintuitive observation due to the fact that this phosphosite has previously been described as an oncogene, and we have measured castration resistance in the cell lines inverse to the amount of p-JNK (Additional file 6: Figure S2) . However, this observation corroborates recent work indicating that JNK acts as an oncogene in tumor development and a tumor suppressor in regards to castration resistant growth .
Modeling the effect of treatments and targeted inhibitors
The effect of complete PI3K inhibition with LY294002 versus mTor inhibition alone with temsirolimus was also examined. Based on the relative survival levels of LNCaP cells treated with LY294002 versus temsirolimus it was determined that the temsirolimus treated group had 31% increased cell survival over cells treated with LY294002. However, both treatments reduced the p-RPS6 to similar levels which were near complete inhibition from basal levels, while LY294002 also strongly reduced measured p-Akt and p-GSK3 levels (Figure 6D and 6E). Based on this observation it was concluded that signaling upstream of mTor (such as p-GSK3 which was observed to be highly correlated to p-Akt) accounted for the difference in survival between complete PI3K inhibition and inhibition of mTor alone.
Modeling the correlation between phosphosites’ activation
The goal of this work was to examine how variation in disparate signaling pathways altered castration resistant growth of three different prostate cancer cell lines in response to activating treatments and targeted inhibitors. In future work, an understanding of how multiple signaling pathways enable castration resistance in patients will be critical to optimizing patient specific treatments using targeted therapies. Differences in the basal level of castration resistant growth across the three cell lines were observed, as was their response to the treatments. A regression model was developed for predicting castration resistant growth and survival, using an MTT assay, which far exceeded randomized data sets (P-value < 0.0001), and was able to account for over half of the variation in cell survival (R2 = 0.616). The MTT assay acted as an approximate metric of cell survival and abstracted the proliferation and apoptosis balance as well as other cellular processes such as neuroendocrine differentiation into one value representing total cell survival in androgen depleted conditions in response to treatment. There are numerous other pathways which are perturbed in prostate cancer by the treatments used here, as well as epigenetic and genetic variability which likely account for the remaining unexplained variance in cell survival, however a majority of cell survival can be explained by these 8 phosphoproteins’ activation level at three time points.
When the effect of androgen treatment on phosphoproteomic signaling was examined we observed an increase in PI3K related phosphoprotein activation (p-Akt, p-RPS6, p-GSK3) at later time points. This is consistent with the observation that AR activation can cause activation of the PI3K pathway, at least in part, through induction of IGF1 secretion . Previous work has indicated that activation of the PI3K pathway can coactivate the AR, causing reciprocal feedback . Additionally, the AR can cause the transcription of cell cycle-related genes directly through binding to the promoter elements and transcribing genes such as c-Myc .
Phosphoprotein levels across cell lines were also examined and there was a clear inverse trend between innate castration resistance and p-JNK levels which did not substantially vary in response to treatment. As previously discussed, this effect may play a role in castration resistance . This variation between cell lines was also seen in the lack of consistent correlation between phosphosites indicating that the genetic and epigenetic differences between the cell lines significantly alters how cell signaling networks respond to treatment. PI3K-related signaling was the only exception to this which had somewhat conserved correlation values across cell lines.
To make additional comparisons PLS regression was performed on the individual cell line data yielding models of cell survival with high R-squared values (LNCaP: R2 = 0.986, PC3: R2 = 0.998, and MDA-PCa-2b: R2 = 0.969). Upon examining the regression coefficients from these models PC3 cells generally weighted positively p-Erk, p-Stat3, p-RPS6, and p-GSK3 as compared to LNCaP which generally weighted p-Erk, p-Stat3, and p-GSK3 positively. Finally, MDA-PCa-2b weighted positively p-Akt, p-RPS6, and p-GSK3 in determining cell survival (Additional file 8: Figure S4). From this data it can be seen that survival appears to be largely mediated through PI3K-related signaling in MDA-PCa-2b cells with an increasing role of p-Erk and p-Stat3 in LNCaP and PC3 cells. Additionally, given the few preserved correlations observed across all cell lines, the data indicates that variations between cell lines cause substantial changes in signaling crosstalk.
The use of targeted kinase inhibitors allowed the elucidation of the role of particular phosphoproteins. Specifically, we identified the role of phosphoproteins upstream of mTor in the PI3K in enabling survival. In a recent phase II clinical trial Temsirolimus as a single agent had an effect on 32% of patients, and numerous PI3K inhibitors are being investigated for use in prostate cancer . Additionally, an increase in p-Erk was noted in response to treatment with the p38 inhibitor SB202190 which is consistent with the observation that p38 inhibition can increase survival [29, 30].
In this work, regression analysis was used to determine how cell signaling correlates with castration resistant growth across three cell lines. Based on the data presented, in vitro prostate cancer cell androgen independent growth could be largely described via MAPK (Mitogen-activated protein kinases) and PI3K signaling. Androgen mediated signaling also largely acted through PI3K signaling. p-JNK appeared to potentially play a role in the fundamental castration resistance of a cell line, and MDA-PCa-2b cells did not utilize p-Erk to enable androgen-independent growth. Given the myriad of targeted inhibitors currently in development, approaches similar to this work which determine drivers of disease progression may enable optimizing treatment given the unique signaling of each patient.
Epidermal growth factor
Extracellular signal-regulated kinases
Glycogen synthase kinase 3
Half maximal inhibitory concentration
Insulin-like growth factor 1
c-Jun N-terminal kinases
Mitogen-activated protein kinases
Mitogen-activated protein kinase kinase
Mammalian target of rapamycin
Partial least squares
Phosphatase and tensin homolog
Ribosomal protein S6
Tumor necrosis factor alpha.
The authors would like to thank Rebecca Lescarbeau for assistance performing regression analysis, and Andrew Moon for assistance taking cell lysates. This work was supported by National Institutes of Health grant P41 EB002520-05 (Tissue Engineering Resource Center) (DL Kaplan).
- Howlader N, Noone AM, Krapcho M, Garshell J, Miller D, Altekruse SF, Kosary CL, Yu M, Ruhl J, Tatalovich Z, Mariotto A, Lewis DR, Chen HS, Feuer EJ, Cronin KA: SEER cancer statistics review. http://seer.cancer.gov/csr/1975_2011/,
- De Bono JS, Logothetis CJ, Molina A, Fizazi K, North S, Chu L, Chi KM, Jones RL, Goodman OB, Saad F, Staffurth JN, Mainwaring P, Harland S, Flaig TW, Hutson TE, Cheng T, Patterson H, Hainsworth JD, Ryan CJ, Sternberg CN, Ellard SL, Flechon A, Saleh M, Scolz M, Kheoh T, Haqq CM, Scher HI: Abiraterone and increased survival in metastatic prostate cancer. N Engl J Med. 2011, 364: 1995-2005. 10.1056/NEJMoa1014618.View ArticlePubMedPubMed CentralGoogle Scholar
- Scher HI, Beer TM, Higano CS, Anand A, Taplin ME, Efstathiou E, Rathkopf D, Shelkey J, Yu EY, Alumkal J, Hung D, Hirmand M, Seely L, Morris MJ, Danila DC, Humm J, Larson S, Fleisher M, Sawyers CL: Antitumor activity of MDV3100 in castration-resistant prostate cancer: a phase 1–2 study. Lancet. 2010, 375: 1437-1446. 10.1016/S0140-6736(10)60172-9.View ArticlePubMedPubMed CentralGoogle Scholar
- Mostaghel EA, Nelson PS: Intracrine androgen metabolism in prostate cancer progression: mechanisms of castration resistance and therapeutic implications. J Clin Endocrinol Metab. 2008, 22: 243-258.Google Scholar
- Wang G, Wang J, Sadar MD: Crosstalk between the androgen receptor and β-Catenin in castrate resistant prostate cancer. Cancer Res. 2008, 68: 9918-9927. 10.1158/0008-5472.CAN-08-1718.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu M, Kyprianou N: Androgen receptor and growth factor signaling cross-talk in prostate cancer cells. Endocr Relat Cancer. 2008, 15: 841-849. 10.1677/ERC-08-0084.View ArticlePubMedPubMed CentralGoogle Scholar
- Alimonti A, Nardella C, Chen Z, Clohessy JG, Carracedo A, Trotman LC, Cheng K, Varmeh S, Kozma SC, Thomas G, Rosivatz E, Woscholski R, Cognetti F, Scher HI, Pandolfi PP: A novel type of cellular senescence that can be enhanced in mouse models and human tumor xenografts to suppress prostate tumorigenesis. J Clin Invest. 2010, 120: 681-693. 10.1172/JCI40535.View ArticlePubMedPubMed CentralGoogle Scholar
- Chetram MA, Odero-Marah V, Hinton CV: Loss of PTEN permits CXCR4-mediated tumorigenesis through ERK1/2 in prostate cancer cells. Mol Cancer Res. 2010, 9: 90-102. 10.1186/1476-4598-9-90.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller-Jensen K, Janes KA, Brugge JS, Lauffenburger DA: Common effector processing mediates cell-specific responses to stimuli. Nature. 2007, 448: 604-608. 10.1038/nature06001.View ArticlePubMedGoogle Scholar
- Kumar N, Wolf-Yadlin A, White F, Lauffenburger DA: Modeling HER2 effects on cell behavior from mass spectrometry phosphotyrosine data. PLoS Comput Biol. 2007, 3: e4-10.1371/journal.pcbi.0030004.View ArticlePubMedPubMed CentralGoogle Scholar
- Saez-Rodriguez J, Alexopoulos LG, Epperlein J, Samaga R, Lauffenburger DA, Klamt S, Sorger PK: Discrete logic modeling as a means to link protein signaling networks with functional analysis of mammalian signal transduction. Mol Syst Biol. 2009, 2: 331-Google Scholar
- Lescarbeau RM, Seib FP, Prewitz M, Werner C, Kaplan DL: In vitro model of metastasis to bone marrow mediates prostate cancer castration resistant growth through paracrine and extracellular matrix factors. PLoS One. 2012, 7: e40372-10.1371/journal.pone.0040372.View ArticlePubMedPubMed CentralGoogle Scholar
- Lopes CT, Franz M, Kazi F, Donaldson SL, Morris Q, Bader GD: Cytoscape web: an interactive web-based network browser. Bioinformatics. 2010, 26: 2347-2348. 10.1093/bioinformatics/btq430.View ArticlePubMedPubMed CentralGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.View ArticlePubMedPubMed CentralGoogle Scholar
- Kinkade CW, Castillo-Martin M, Puzio-Kuter A, Yan J, Foster TH, Gao H, Sun Y, Ouyang X, Gerald WL, Cordon-Cardo C, Abate-Shen C: Targeting AKT/mTOR and ERK MAPK signaling inhibits hormone-refractory prostate cancer in a preclinical mouse model. J Clin Invest. 2008, 118: 3051-3064.PubMedPubMed CentralGoogle Scholar
- Floc’h N, Kinkade CW, Kobayashi T, Aytes A, Lefebvre C, Mitrofanova A, Cardiff RD, Califano A, Shen MM, Abate-Shen C: Dual targeting of the Akt/mTOR signaling pathway inhibits castration-resistant prostate cancer in a genetically engineered mouse model. Cancer Res. 2012, 72: 4483-4493. 10.1158/0008-5472.CAN-12-0283.View ArticlePubMedPubMed CentralGoogle Scholar
- Schutz S, Schrader AJ, Zengerling F, Genze F, Cronauer MV, Schrader M: Inhibition of glycogen synthase kinase-3β counteracts ligand-independent activity of the androgen receptor in castration resistant prostate cancer. PLoS One. 2011, 6: e25341-10.1371/journal.pone.0025341.View ArticlePubMedPubMed CentralGoogle Scholar
- Gan L, Wang J, Xu H, Yang X: Resistance to docetaxel-induced apoptosis in prostate cancer cells by p38/p53/p21 signaling. Prostate. 2011, 71: 1158-1166. 10.1002/pros.21331.View ArticlePubMedGoogle Scholar
- Hubner A, Mulholland DJ, Standen CL, Karasarides M, Cavanagh-Kyros J, Barrett T, Chi H, Greiner DL, Tournier C, Sawyers CL, Flavell RA, Wu H, Davis RJ: JNK and PTEN cooperatively control the development of invasive adenocarcinoma of the prostate. Proc Natl Acad Sci USA. 2012, 109: 12046-12051. 10.1073/pnas.1209660109.View ArticlePubMedPubMed CentralGoogle Scholar
- Zoubeidi A, Zardan A, Wiedmann RM, Locke J, Beraldi E, Fazli L, Gleave ME: Hsp27 promotes insulin-like growth factor-I survival signaling in prostate cancer via p90Rsk-dependent phosphorylation and inactivation of BAD. Cancer Res. 2010, 70: 2307-2317. 10.1158/0008-5472.CAN-09-3252.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee SO, Lou W, Johnson CS, Trump DL, Gao AC: Interlukin-6 protects LNCaP cells from apoptosis induced by androgen deprivation through the Stat3 pathway. Prostate. 2004, 60: 178-186. 10.1002/pros.20045.View ArticlePubMedGoogle Scholar
- Vellaichamy A, Dezso Z, JeBailey L, Chinnaiyan AM, Sreekumar A, Nesvizhskii A, Omenn GS, Bugrim A: Topological significance: analysis of gene expression and proteomic profiles from prostate cancer cells reveals key mechanisms of androgen response. PLoS One. 2009, 5: e10936-View ArticleGoogle Scholar
- Bonaccorsi L, Carloni V, Muratori M, Salvadori A, Giannini A, Carini M, Serio M, Forti G, Baldi E: Androgen receptor expression in prostate carcinoma cells suppresses alpha6beta4 integrin-mediated invasive phenotype. Endocrinol. 2000, 141: 3172-3182.Google Scholar
- Gosh PM, Malik SN, Bedolla RG, Wang Y, MIkhailova M, Prihoda TJ, Troyer DA, Kreisberg JI: Signal transduction pathways in androgen-dependent and independent prostate cancer cell proliferation. Endocr Relat Cancer. 2005, 12: 119-134. 10.1677/erc.1.00835.View ArticleGoogle Scholar
- Leppa S, Bohmann D: Diverse functions of JNK and c-Jun in stress response and apoptosis. Oncogene. 1999, 18: 6158-6162. 10.1038/sj.onc.1203173.View ArticlePubMedGoogle Scholar
- Carver BS, Chapinski C, Wongvipat J, Hieronymus H, Chen Y, Chanadarlapaty S, Arora VK, Le C, Koutcher J, Scher H, Scardino PT, Rosen N, Sawyers CL: Reciprocal feedback regulation of PI3K and androgen receptor signaling in PTEN-deficient prostate cancer. Cancer Cell. 2011, 19: 575-586. 10.1016/j.ccr.2011.04.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Shiota M, Yokomizo A, Naito S: Increased androgen receptor transcription: a cause of castration-resistant prostate cancer and a possible therapeutic target. J Mol Endocrinol. 2011, 47: R25-R41. 10.1530/JME-11-0018.View ArticlePubMedGoogle Scholar
- Nabhan C, Kruczek K, Tolzien K, Galvez AG, Lestingi TM: A phase II study evaluating the toxicity and efficacy of single-agent temsorilmus (TEM) in chemotherapy-naive castration-resistant prostate cancer (CRPC): first report of suggested activity. J Clin Oncol. 2012, 5: 165-Google Scholar
- Shimo T, Matsumura S, Ibaragi S, Isowa S, Kishimoto K, Mese H, Nishiyama A, Sasaki A: Specific inhibitor of MEK-mediated cross-talk between ERK and p38 MAPK during differentiation of human osteosarcoma cells. J Cell Commun Signal. 2007, 1: 103-111. 10.1007/s12079-007-0010-2.View ArticlePubMedPubMed CentralGoogle Scholar
- Hotokezaka H, Sakai E, Kanaoka K, Saito K, Matsuo K, Kitaura H, Yoshida N, Nakayama K: U0126 and PD98059, specific inhibitors of Mek, accelerate differentiation of RAW264.7 cells into osteoclast-like cells. J Biol Chem. 2002, 277: 47366-47372. 10.1074/jbc.M208284200.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/14/325/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.