On the surviving fraction in irradiated multicellular tumour spheroids: calculation of overall radiosensitivity parameters, influence of hypoxia and volume effects

We model the heterogeneous response to radiation of multicellular tumour spheroids assuming position- and volume-dependent radiosensitivity. We propose a method to calculate the overall radiosensitivity parameters to obtain the surviving fraction of tumours. A mathematical model of a spherical tumour with a hypoxic core and a viable rim which is a caricature of a real tumour is constructed. The model is embedded in a two-compartment linear-quadratic (LQ) model, assuming a mixed bivariated Gaussian distribution to attain the radiosensitivity parameters. Ergodicity, i.e., the equivalence between ensemble and volumetric averages is used to obtain the overall radiosensitivities for the two compartments. We obtain expressions for the overall radiosensitivity parameters resulting from the use of both a linear and a nonlinear dependence of the local radiosensitivity with position. The model's results are compared with experimental data of surviving fraction (SF) for multicellular spheroids of different sizes. We make one fit using only the smallest spheroid data and we are able to predict the SF for the larger spheroids. These predictions are acceptable particularly using bounded sensitivities. We conclude with the importance of taking into account the contribution of clonogenic hypoxic cells to radiosensitivity and with the convenience of using bounded local sensitivities to predict overall radiosensitivity parameters