No abstract

There is growing clinical interest in the use of 125I (half-life 59.4 days) and 103Pd (half-life 16.97 days) for permanent brachytherapy implants. These radionuclides pose interesting radiobiological challenges because, even with slowly growing tumours, significant tumour cell repopulation may occur during the long period taken to deliver the full radiation dose. This results in a considerable amount of the prescribed dose being wasted. There may also be changes in the tumour volume during treatment (due to oedema and/or shrinkage), thus altering the relative geometry of the implanted seeds and causing additional dose rate variations. This assessment examines the interaction between the above effects and additionally includes allowance for the influence of the relative biological effectiveness (RBE) of the radiations emitted by the two radionuclides. The results are presented in terms of the biologically effective doses (BEDs) and likely tumour control probabilities (TCPs) associated with the various parameter combinations. The overall BED enhancement due to the RBE effect is shown always to be greater than the RBE itself and is greatest in tumours which are radio-resistive and/or fast growing. The biological dose uncertainties are found to be less with 103Pd and the TCPs associated with this radionuclide are expected to be significantly higher in the treatment of some 'difficult' tumours. Using typically prescribed doses 125i appears to be better for treating radiosensitive tumours with long doubling times and which shrink fairly rapidly. However, unless 125I doses are reduced, this advantage may well be offset by the greatly enhanced biological doses delivered to adjacent normal structures.

Tumours behave as complex, self-organizing, opportunistic dynamic systems. In an attempt to better understand and describe the highly complicated tumour behaviour, a novel four-dimensional simulation model of in vivo tumour growth and response to radiotherapy has been developed. This paper presents the latest improvements to the model as well as a parametric validation of it. Improvements include an advanced algorithm leading to conformal tumour shrinkage, a quantitative consideration of the influence of oxygenation on radiosensitivity and a more realistic, imaging based description of the neovasculature distribution. The tumours selected for the validation of the model are a wild type and a mutated p53 gene glioblastomas multiforme. According to the model predictions, a whole tumour with larger cell cycle duration tends to repopulate more slowly. A lower oxygen enhancement ratio value leads to a more radiosensitive whole tumour. Higher clonogenic cell density (CCD) produces a higher number of proliferating tumour cells and, therefore, a more difficult tumour to treat. Simulation predictions agree at least semi-quantitatively with clinical experience, and particularly with the outcome of the Radiation Therapy Oncology Group (RTOG) Study 83-02. It is stressed that the model allows a quantitative study of the interrelationship between the competing influences in a complex, dynamic tumour environment. Therefore, the model can already be useful as an educational tool with which to study, understand and demonstrate the role of various parameters in tumour growth and response to irradiation. A long term quantitative clinical adaptation and validation of the model aiming at its integration into the treatment planning procedure is in progress.

Advanced bio-simulation methods are expected to substantially improve radiotherapy treatment planning. To this end a novel spatio-temporal patient-specific simulation model of the in vivo response of malignant tumours to radiotherapy schemes has been recently developed by our group. This paper discusses recent improvements to the model: an optimized algorithm leading to conformal shrinkage of the tumour as a response to radiotherapy, the introduction of the oxygen enhancement ratio (OER), a realistic initial cell phase distribution and finally an advanced imaging-based algorithm simulating the neovascularization field. A parametric study of the influence of the cell cycle duration Tc, OER, OERbeta for the beta LQ parameter on tumour growth. shrinkage and response to irradiation under two different fractionation schemes has been made. The model has been applied to two glioblastoma multiforme (GBM) cases, one with wild type (wt) and another one with mutated (mt) p53 gene. Furthermore, the model has been applied to a hypothetical GBM tumour with alpha and beta values corresponding to those of generic radiosensitive tumours. According to the model predictions, a whole tumour with shorter Tc tends to repopulate faster, as is to be expected. Furthermore, a higher OER value for the dormant cells leads to a more radioresistant whole tumour. A small variation of the OERbeta value does not seem to play a major role in the tumour response. Accelerated fractionation proved to be superior to the standard scheme for the whole range of the OER values considered. Finally, the tumour with mt p53 was shown to be more radioresistant compared to the tumour with wt p53. Although all simulation predictions agree at least qualitatively with the clinical experience and literature, a long-term clinical adaptation and quantitative validation procedure is in progress.

A novel four-dimensional, patient-specific Monte Carlo simulation model of solid tumor response to chemotherapeutic treatment in vivo is presented. The special case of glioblastoma multiforme treated by temozolomide is addressed as a simulation paradigm. Nevertheless, a considerable number of the involved algorithms are generally applicable. The model is based on the patient's imaging, histopathologic and genetic data. For a given drug administration schedule lying within acceptable toxicity boundaries, the concentration of the prodrug and its metabolites within the tumor is calculated as a function of time based on the drug pharamacokinetics. A discretization mesh is superimposed upon the anatomical region of interest and within each geometrical cell of the mesh the most prominent biological "laws" (cell cycling, necrosis, apoptosis, mechanical restictions, etc.) are applied. The biological cell fates are predicted based on the drug pharmacodynamics. The outcome of the simulation is a prediction of the spatiotemporal activity of the entire tumor and is virtual reality visualized. A good qualitative agreement of the model's predictions with clinical experience supports the applicability of the approach. The proposed model primarily aims at providing a platform for performing patient individualized in silico experiments as a means of chemotherapeutic treatment optimization.