Monte Carlo radiotherapy simulations of accelerated repopulation and reoxygenation for hypoxic head and neck cancer

Abstract

A temporal Monte Carlo tumour growth and radiotherapy effect model (HYP-RT) simulating hypoxia in head and neck cancer has been developed and used to analyse parameters influencing cell kill during conventionally fractionated radiotherapy. The model was designed to simulate individual cell division up to 10(8) cells, while incorporating radiobiological effects, including accelerated repopulation and reoxygenation during treatment. Reoxygenation of hypoxic tumours has been modelled using randomised increments of oxygen to tumour cells after each treatment fraction. The process of accelerated repopulation has been modelled by increasing the symmetrical stem cell division probability. Both phenomena were onset immediately or after a number of weeks of simulated treatment. The extra dose required to control (total cell kill) hypoxic vs oxic tumours was 15-25% (8-20 Gy for 5 × 2 Gy per week) depending on the timing of accelerated repopulation onset. Reoxygenation of hypoxic tumours resulted in resensitisation and reduction in total dose required by approximately 10%, depending on the time of onset. When modelled simultaneously, accelerated repopulation and reoxygenation affected cell kill in hypoxic tumours in a similar manner to when the phenomena were modelled individually; however, the degree was altered, with non-additive results. Simulation results were in good agreement with standard linear quadratic theory; however, differed for more complex comparisons where hypoxia, reoxygenation as well as accelerated repopulation effects were considered. Simulations have quantitatively confirmed the need for patient individualisation in radiotherapy for hypoxic head and neck tumours, and have shown the benefits of modelling complex and dynamic processes using Monte Carlo methods.