Abstract Purpose: Bone is the most predominant site of distant metastasis in prostate cancer, and patients have limited therapeutic options at this stage. Experimental Design: We performed a system-wide quantitative proteomic analysis of bone metastatic prostate tumors from 22 patients operated to relieve spinal cord compression. At the time of surgery, most patients had relapsed after androgen-deprivation therapy, while 5 were previously untreated. An extended cohort of prostate cancer bone metastases (n = 65) was used for immunohistochemical validation. Results: On average, 5,067 proteins were identified and quantified per tumor. Compared with primary tumors (n = 26), bone metastases were more heterogeneous and showed increased levels of proteins involved in cell-cycle progression, DNA damage response, RNA processing, and fatty acid β-oxidation; and reduced levels of proteins were related to cell adhesion and carbohydrate metabolism. Within bone metastases, we identified two phenotypic subgroups: BM1, expressing higher levels of AR canonical targets, and mitochondrial and Golgi apparatus resident proteins; and BM2, with increased expression of proliferation and DNA repair–related proteins. The two subgroups, validated by the inverse correlation between MCM3 and prostate specific antigen immunoreactivity, were related to disease prognosis, suggesting that this molecular heterogeneity should be considered when developing personalized therapies. Conclusions: This work is the first system-wide quantitative characterization of the proteome of prostate cancer bone metastases and a valuable resource for understanding the etiology of prostate cancer progression. Clin Cancer Res; 24(21); 5433–44. ©2018 AACR.