Background: In metastatic prostate cancer (PCa) patients, androgen-deprivation therapy (ADT) combined with chemotherapy or next-generation androgen receptor targeted agents is a new standard treatment. The objective of the present study is to assess longitudinal PSA kinetics during treatment using mathematical modeling, to identify the modeled PSA kinetic parameters able to exhibit early prognostic/predictive values. Methods: Phase III clinical trial dataset (NCT00764166) comparing ADT +/− docetaxel in 250 locally treated patients for PCa with rising PSA levels, who were at high risk of metastatic disease was assessed. A kinetic-pharmacodynamic (K-PD) model was used to fit PSA kinetics during the first 100 treatment days, to estimate the modeled PSA production rate K (KPROD) and elimination constant rate K (KELIM). The prognostic value of these parameters, considered as categorized (favorable vs. unfavorable) covariates regarding PSA progression-free survival (PSA-PFS) and overall survival (OS), was assessed using univariate/multivariate analyses. Results: Data from 177/250 patients was assessed. KELIM exhibited a significant prognostic value regarding PSA-PFS and KPROD regarding OS (univariate analysis). In the PSA-PFS final multivariate model, KELIM and the primary therapy type were significant. The OS multivariate model integrated both KPROD and baseline PSA doubling-time. Conclusion: In this first study assessing the modeled PSA kinetics prognostic value in PCa patients treated with systemic treatments, KELIM and KPROD exhibited respective prognostic values regarding PSA-PFS and OS.