Tumor heterogeneity is widely considered to be a determinant factor in tumor progression and in particular in its recurrence after therapy. Unfortunately, current medical techniques are unable to deduce clinically relevant information about tumor heterogeneity by means of non-invasive methods. As a consequence, when radiotherapy is used as a treatment of choice, radiation dosimetries are prescribed under the assumption that the malignancy targeted is of a homogeneous nature. In this work we discuss the effects of different radiation dose distributions on heterogeneous tumors by means of an individual cell-based model. To that end, a case is considered where two tumor cell phenotypes are present, which we assume to strongly differ in their respective cell cycle duration and radiosensitivity properties. We show herein that, as a result of such differences, the spatial distribution of the corresponding phenotypes, whence the resulting tumor heterogeneity can be predicted as growth proceeds. In particular, we show that if we start from a situation where a majority of ordinary cancer cells (CCs) and a minority of cancer stem cells (CSCs) are randomly distributed, and we assume that the length of CSC cycle is significantly longer than that of CCs, then CSCs become concentrated at an inner region as tumor grows. As a consequence we obtain that if CSCs are assumed to be more resistant to radiation than CCs, heterogeneous dosimetries can be selected to enhance tumor control by boosting radiation in the region occupied by the more radioresistant tumor cell phenotype. It is also shown that, when compared with homogeneous dose distributions as those being currently delivered in clinical practice, such heterogeneous radiation dosimetries fare always better than their homogeneous counterparts. Finally, limitations to our assumptions and their resulting clinical implications will be discussed.

The choice of any radiotherapy treatment plan is usually made after the evaluation of a few preliminary isodose distributions obtained from different beam configurations. Despite considerable advances in planning techniques, such final decision remains a challenging task that would greatly benefit from efficient and reliable assessment tools. For any dosimetric plan considered, data on dose-volume histograms supplied by treatment planning systems are used to provide estimates on planning target coverage as well as on sparing of organs at risk and the remaining healthy tissue. These partial metrics are then combined into a dose distribution index (DDI), which provides a unified, easy-to-read score for each competing radiotherapy plan. To assess the performance of the proposed scoring system, DDI figures for fifty brain cancer patients were retrospectively evaluated. Patients were divided in three groups depending on tumor location and malignancy. For each patient, three tentative plans were designed and recorded during planning, one of which was eventually selected for treatment. We thus were able to compare the plans with better DDI scores and those actually delivered. When planning target coverage and organs at risk sparing are considered as equally important, the tentative plan with the highest DDI score is shown to coincide with that actually delivered in 32 of the 50 patients considered. In 15 (respectively 3) of the remaining 18 cases, the plan with highest DDI value still coincides with that actually selected, provided that organs at risk sparing is given higher priority (respectively, lower priority) than target coverage. DDI provides a straightforward and non-subjective tool for dosimetric comparison of tentative radiotherapy plans. In particular, DDI readily quantifies differences among competing plans with similar-looking dose-volume histograms and can be easily implemented for any tumor type and localization, irrespective of the planning system and irradiation technique considered. Moreover, DDI permits to estimate the dosimetry impact of different priorities being assigned to sparing of organs at risk or to better target coverage.

In this work a mathematical model for blood coagulation induced by an activator source is presented. Blood coagulation is viewed as a process resulting in fibrin polymerization, which is considered as the first step towards thrombi formation. We derive and study a system for the first moments of the polymer concentrations and the activating variables. Analysis of this last model allows us to identify parameter regions which could lead to thrombi formation, both in homeostatic and pathological situations.