Mechanistic modelling of prostate-specific antigen dynamics shows potential for personalized prediction of radiation therapy outcome

External beam radiation therapy is a widespread treatment for prostate cancer. The ensuing patient follow-up is based on the evolution of the prostate-specific antigen (PSA). Serum levels of PSA decay due to the radiation-induced death of tumour cells and cancer recurrence usually manifest as a rising PSA. The current definition of biochemical relapse requires that PSA reaches nadir and starts increasing, which delays the use of further treatments. Also, these methods do not account for the post-radiation tumour dynamics that may contain early information on cancer recurrence. Here, we develop three mechanistic models of post-radiation PSA evolution. Our models render superior fits of PSA data in a patient cohort and provide a biological justification for the most common empirical formulation of PSA dynamics. We also found three model-based prognostic variables: the proliferation rate of the survival fraction, the ratio of radiation-induced cell death rate to the survival proliferation rate, and the time to PSA nadir since treatment termination. We argue that these markers may enable the early identification of biochemical relapse, which would permit physicians to subsequently adapt patient monitoring to optimize the detection and treatment of cancer recurrence.

[1]  Alicia Martínez-González,et al.  Delay effects in the response of low-grade gliomas to radiotherapy: a mathematical model and its therapeutical implications. , 2014, Mathematical medicine and biology : a journal of the IMA.

[2]  J. Murray,et al.  A quantitative model for the dynamics of serum prostate-specific antigen as a marker for cancerous growth: an explanation for a medical anomaly. , 2001, The American journal of pathology.

[3]  M. Nowak,et al.  Tumour and immune cell dynamics explain the PSA bounce after prostate cancer brachytherapy , 2016, British Journal of Cancer.

[4]  D. Cui,et al.  An integrative model of prostate cancer interaction with the bone microenvironment. , 2017, Mathematical biosciences.

[5]  Trachette L. Jackson,et al.  A mathematical investigation of the multiple pathways to recurrent prostate cancer: comparison with experimental data. , 2004, Neoplasia.

[6]  R. Vollmer,et al.  The dynamics of prostate-specific antigen after definitive radiation therapy for prostate cancer. , 1999, Clinical cancer research : an official journal of the American Association for Cancer Research.

[7]  A. Pollack,et al.  Kinetics of serum prostate-specific antigen after external beam radiation for clinically localized prostate cancer. , 1997, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[8]  B. Ross,et al.  Mathematical Modeling of PDGF-Driven Glioblastoma Reveals Optimized Radiation Dosing Schedules , 2014, Cell.

[9]  Y. Yamada,et al.  The natural history and predictors of outcome following biochemical relapse in the dose escalation era for prostate cancer patients undergoing definitive external beam radiotherapy. , 2015, European urology.

[10]  E. Venkatraman,et al.  Outcome predictors for the increasing PSA state after definitive external-beam radiotherapy for prostate cancer. , 2005, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[11]  Johan Pallud,et al.  Dynamic imaging response following radiation therapy predicts long-term outcomes for diffuse low-grade gliomas. , 2012, Neuro-oncology.

[12]  C. Pinnock,et al.  Prostate-specific antigen (PSA) rate of decline post external beam radiotherapy predicts prostate cancer death. , 2013, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[13]  Andrew J Vickers,et al.  Systematic review of pretreatment prostate-specific antigen velocity and doubling time as predictors for prostate cancer. , 2009, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[14]  K Hendrickson,et al.  Predicting the efficacy of radiotherapy in individual glioblastoma patients in vivo: a mathematical modeling approach , 2010, Physics in medicine and biology.

[15]  Yang Kuang,et al.  Mechanisms of resistance to intermittent androgen deprivation in patients with prostate cancer identified by a novel computational method. , 2014, Cancer research.

[16]  Menggang Yu,et al.  Individualized Predictions of Disease Progression Following Radiation Therapy for Prostate Cancer , 2005 .

[17]  Guillermo Lorenzo,et al.  Tissue-scale, personalized modeling and simulation of prostate cancer growth , 2016, Proceedings of the National Academy of Sciences.

[18]  Kazuyuki Aihara,et al.  Development of a mathematical model that predicts the outcome of hormone therapy for prostate cancer. , 2010, Journal of theoretical biology.

[19]  Jacob G. Scott,et al.  A genome-based model for adjusting radiotherapy dose (GARD): a retrospective, cohort-based study. , 2017, The Lancet. Oncology.

[20]  D. Kuban,et al.  PSA nadir predicts biochemical and distant failures after external beam radiotherapy for prostate cancer: a multi-institutional analysis. , 2004, International journal of radiation oncology, biology, physics.

[21]  M. Cooperberg,et al.  Temporal Trends and the Impact of Race, Insurance, and Socioeconomic Status in the Management of Localized Prostate Cancer. , 2017, European urology.

[22]  Xavier Filella,et al.  Emerging biomarkers in the diagnosis of prostate cancer , 2018, Pharmacogenomics and personalized medicine.

[23]  Robin T Vollmer,et al.  Dissecting the dynamics of serum prostate-specific antigen. , 2010, American journal of clinical pathology.

[24]  J. S. St. Sauver,et al.  Natural history of benign prostatic enlargement: long‐term longitudinal population‐based study of prostate volume doubling times , 2010, BJU international.

[25]  X. Allen Li,et al.  HOW LOW IS THE / RATIO FOR PROSTATE CANCER? , 2003 .

[26]  M. Eble,et al.  Prostate-specific antigen kinetics following external-beam radiotherapy and temporary (Ir-192) or permanent (I-125) brachytherapy for prostate cancer. , 2010, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[27]  Renaud de Crevoisier,et al.  Towards an integrative computational model for simulating tumor growth and response to radiation therapy , 2017, Symposium on Medical Information Processing and Analysis.

[28]  A. Hanlon,et al.  Modeling postradiation prostate specific antigen level kinetics , 1998, Cancer.

[29]  C. Fuller,et al.  Early prostate‐specific antigen (PSA) kinetics following prostate carcinoma radiotherapy , 2004, Cancer.

[30]  R W Veltri,et al.  Implication of cell kinetic changes during the progression of human prostatic cancer. , 1995, Clinical cancer research : an official journal of the American Association for Cancer Research.

[31]  P. Carroll,et al.  Prognostic value of Ki67 in localized prostate carcinoma: a multi-institutional study of >1000 prostatectomies , 2016, Prostate Cancer and Prostatic Disease.

[32]  David A Hormuth,et al.  Biophysical Modeling of In Vivo Glioma Response After Whole-Brain Radiation Therapy in a Murine Model of Brain Cancer. , 2017, International journal of radiation oncology, biology, physics.

[33]  Jian Z. Wang,et al.  Impact of tumor repopulation on radiotherapy planning. , 2005, International journal of radiation oncology, biology, physics.

[34]  Marcial García-Rojo,et al.  Non-standard radiotherapy fractionations delay the time to malignant transformation of low-grade gliomas , 2017, PloS one.

[35]  S. F. C. O’Rourke,et al.  Linear quadratic and tumour control probability modelling in external beam radiotherapy , 2009, Journal of mathematical biology.

[36]  Gouhei Tanaka,et al.  A Mathematical Model of Intermittent Androgen Suppression for Prostate Cancer , 2008, J. Nonlinear Sci..

[37]  J. Oden,et al.  Selection and Validation of Predictive Models of Radiation Effects on Tumor Growth Based on Noninvasive Imaging Data. , 2017, Computer methods in applied mechanics and engineering.

[38]  Mark A. J. Chaplain,et al.  Towards Predicting the Response of a Solid Tumour to Chemotherapy and Radiotherapy Treatments: Clinical Insights from a Computational Model , 2013, PLoS Comput. Biol..

[39]  Kristin R. Swanson,et al.  Toward Patient-Specific, Biologically Optimized Radiation Therapy Plans for the Treatment of Glioblastoma , 2013, PloS one.

[40]  T. Pickles,et al.  PSA doubling time kinetics during prostate cancer biochemical relapse after external beam radiation therapy. , 2005, International journal of radiation oncology, biology, physics.

[41]  H. Kal,et al.  HOW LOW IS THE / RATIO FOR PROSTATE CANCER? , 2003 .

[42]  Avner Friedman,et al.  Mathematical modeling of prostate cancer progression in response to androgen ablation therapy , 2011, Proceedings of the National Academy of Sciences.

[43]  J. Davies,et al.  Molecular Biology of the Cell , 1983, Bristol Medico-Chirurgical Journal.

[44]  J. Grasman,et al.  Comparison to Experimental Data , 2019, Relativistic Fluid Dynamics in and out of Equilibrium.

[45]  Martin H Kroll,et al.  Kinetics of tumor growth of prostate carcinoma estimated using prostate-specific antigen. , 2005, Urology.

[46]  Albert Lai,et al.  Prognostic significance of growth kinetics in newly diagnosed glioblastomas revealed by combining serial imaging with a novel biomathematical model. , 2009, Cancer research.

[47]  D. Kuban,et al.  First-year PSA kinetics and minima after prostate cancer radiotherapy are predictive of overall survival. , 2005, International journal of radiation oncology, biology, physics.

[48]  Laurent Pujo-Menjouet,et al.  Mathematical models of radiation action on living cells: From the target theory to the modern approaches. A historical and critical review. , 2016, Journal of theoretical biology.

[49]  Eduardo G Moros,et al.  The Evolution of Tumour Composition During Fractionated Radiotherapy: Implications for Outcome , 2017, Bulletin of Mathematical Biology.

[50]  I. Kaplan,et al.  Prostate-specific antigen kinetics after external beam irradiation for carcinoma of the prostate. , 1994, International journal of radiation oncology, biology, physics.

[51]  Tumor volume in prostate cancer and serum prostate-specific antigen. Analysis from a kinetic viewpoint. , 2003 .

[52]  Thomas E Yankeelov,et al.  Predicting the Response of Breast Cancer to Neoadjuvant Therapy Using a Mechanically Coupled Reaction-Diffusion Model. , 2015, Cancer research.

[53]  M. Eble,et al.  Long-term prognostic significance of rising PSA levels following radiotherapy for localized prostate cancer – focus on overall survival , 2017, Radiation oncology.