A four-dimensional simulation model of tumour response to radiotherapy in vivo: parametric validation considering radiosensitivity, genetic profile and fractionation

A Multidisciplinary Hyper-Modeling Scheme in Personalized In Silico Oncology: Coupling Cell Kinetics with Metabolism, Signaling Networks, and Biomechanics as Plug-In Component Models of a Cancer Digital Twin

The massive amount of human biological, imaging, and clinical data produced by multiple and diverse sources necessitates integrative modeling approaches able to summarize all this information into answers to specific clinical questions. In this paper, we present a hypermodeling scheme able to combine models of diverse cancer aspects regardless of their underlying method or scale. Describing tissue-scale cancer cell proliferation, biomechanical tumor growth, nutrient transport, genomic-scale aberrant cancer cell metabolism, and cell-signaling pathways that regulate the cellular response to therapy, the hypermodel integrates mutation, miRNA expression, imaging, and clinical data. The constituting hypomodels, as well as their orchestration and links, are described. Two specific cancer types, Wilms tumor (nephroblastoma) and non-small cell lung cancer, are addressed as proof-of-concept study cases. Personalized simulations of the actual anatomy of a patient have been conducted. The hypermodel has also been applied to predict tumor control after radiotherapy and the relationship between tumor proliferative activity and response to neoadjuvant chemotherapy. Our innovative hypermodel holds promise as a digital twin-based clinical decision support system and as the core of future in silico trial platforms, although additional retrospective adaptation and validation are necessary.

Intraoperative navigation system with a multi-modality fusion of 3D virtual model and laparoscopic real-time images in laparoscopic pancreatic surgery: a preclinical study

BACKGROUND

Laparoscopy is widely used in pancreatic surgeries nowadays. The efficient and correct judgment of the location of the anatomical structures is crucial for a safe laparoscopic pancreatic surgery. The technologies of 3-dimensional(3D) virtual model and image fusion are widely used for preoperative planning and intraoperative navigation in the medical field, but not in laparoscopic pancreatic surgery up to now. We aimed to develop an intraoperative navigation system with an accurate multi-modality fusion of 3D virtual model and laparoscopic real-time images for laparoscopic pancreatic surgery.

METHODS

The software for the navigation system was developed ad hoc. The preclinical study included tests with the laparoscopic simulator and pilot cases. The 3D virtual models were built using preoperative Computed Tomography (CT) Digital Imaging and Communications in Medicine (DICOM) data. Manual and automatic real-time image fusions were tested. The practicality of the navigation system was evaluated by the operators using the National Aeronautics and Space Administration-Task Load Index (NASA-TLX) method.

RESULTS

The 3D virtual models were successfully built using the navigation system. The 3D model was correctly fused with the real-time laparoscopic images both manually and automatically optical orientation in the preclinical tests. The statistical comparative tests showed no statistically significant differences between the scores of the rigid model and those of the phantom model(P > 0.05). There was statistically significant difference between the total scores of automatic fusion function and those of manual fusion function (P = 0.026). In pilot cases, the 3D model was correctly fused with the real-time laparoscopic images manually. The Intraoperative navigation system was easy to use. The automatic fusion function brought more convenience to the user.

CONCLUSIONS

The intraoperative navigation system applied in laparoscopic pancreatic surgery clearly and correctly showed the covered anatomical structures. It has the potentiality of helping achieve a more safe and efficient laparoscopic pancreatic surgery.

A Fractional Analysis of Hyperthermia Therapy on Breast Cancer in a Porous Medium along with Radiative Microwave Heating

Cancer is a prominent source of mortality and morbidity globally, but little is known about how it develops and spreads. Tumor cells are unable to thrive in high-temperature environments, according to recent research. Hyperthermia is the name for this therapy method. This study provides insights into hyperthermia therapy on breast cancer in the presence of a porous material with fractional derivative access when using radiative microwave heating. The mathematical model is formulated by PDE, while the time-fractional Caputo derivative is applied to make our equation more general as compared to the classical model. To produce a more efficient analysis of blood temperature distributions inside the tissues of the breast, the unsteady state is calculated by using the Laplace transform technique. The Laplace inversion is found by Durbin’s and Zakian’s algorithms. The treatment involves mild temperature hyperthermia, which causes cell death by enhancing cell sensitivity to radiation therapy and blood flow in the tumor. The variations of different parameters to control the temperate profile during therapy are discussed; we can also see how a fractional parameter makes our study more realistic for further experimental study.

In Silico Mathematical Modelling for Glioblastoma: A Critical Review and a Patient-Specific Case

Glioblastoma extensively infiltrates the brain; despite surgery and aggressive therapies, the prognosis is poor. A multidisciplinary approach combining mathematical, clinical and radiological data has the potential to foster our understanding of glioblastoma evolution in every single patient, with the aim of tailoring therapeutic weapons. In particular, the ultimate goal of biomathematics for cancer is the identification of the most suitable theoretical models and simulation tools, both to describe the biological complexity of carcinogenesis and to predict tumor evolution. In this report, we describe the results of a critical review about different mathematical models in neuro-oncology with their clinical implications. A comprehensive literature search and review for English-language articles concerning mathematical modelling in glioblastoma has been conducted. The review explored the different proposed models, classifying them and indicating the significative advances of each one. Furthermore, we present a specific case of a glioblastoma patient in which our recently proposed innovative mechanical model has been applied. The results of the mathematical models have the potential to provide a relevant benefit for clinicians and, more importantly, they might drive progress towards improving tumor control and patient's prognosis. Further prospective comparative trials, however, are still necessary to prove the impact of mathematical neuro-oncology in clinical practice.

Three-dimensional tumor growth in time-varying chemical fields: a modeling framework and theoretical study

Background

Contemporary biological observations have revealed a large variety of mechanisms acting during the expansion of a tumor. However, there are still many qualitative and quantitative aspects of the phenomenon that remain largely unknown. In this context, mathematical and computational modeling appears as an invaluable tool providing the means for conducting in silico experiments, which are cheaper and less tedious than real laboratory experiments.

Results

This paper aims at developing an extensible and computationally efficient framework for in silico modeling of tumor growth in a 3-dimensional, inhomogeneous and time-varying chemical environment. The resulting model consists of a set of mathematically derived and algorithmically defined operators, each one addressing the effects of a particular biological mechanism on the state of the system. These operators may be extended or re-adjusted, in case a different set of starting assumptions or a different simulation scenario needs to be considered.

Conclusion

In silico modeling provides an alternative means for testing hypotheses and simulating scenarios for which exact biological knowledge remains elusive. However, finer tuning of pertinent methods presupposes qualitative and quantitative enrichment of available biological evidence. Validation in a strict sense would further require comprehensive, case-specific simulations and detailed comparisons with biomedical observations.

Studying the regression profiles of cervical tumours during radiotherapy treatment using a patient-specific multiscale model

Apart from offering insight into the biomechanisms involved in cancer, many recent mathematical modeling efforts aspire to the ultimate goal of clinical translation, wherein models are designed to be used in the future as clinical decision support systems in the patient-individualized context. Most significant challenges are the integration of multiscale biodata and the patient-specific model parameterization. A central aim of this study was the design of a clinically-relevant parameterization methodology for a patient-specific computational model of cervical cancer response to radiotherapy treatment with concomitant cisplatin, built around a tumour features-based search of the parameter space. Additionally, a methodological framework for the predictive use of the model was designed, including a scoring method to quantitatively reflect the similarity and bilateral predictive ability of any two tumours in terms of their regression profile. The methodology was applied to the datasets of eight patients. Tumour scenarios in accordance with the available longitudinal data have been determined. Predictive investigations identified three patient cases, anyone of which can be used to predict the volumetric evolution throughout therapy of the tumours of the other two with very good results. Our observations show that the presented approach is promising in quantifiably differentiating tumours with distinct regression profiles.

Mathematical Modeling of the Effects of Tumor Heterogeneity on the Efficiency of Radiation Treatment Schedule

Radiotherapy uses high doses of energy to eradicate cancer cells and control tumors. Various treatment schedules have been developed and tried in clinical trials, yet significant obstacles remain to improving the radiotherapy fractionation. Genetic and non-genetic cellular diversity within tumors can lead to different radiosensitivity among cancer cells that can affect radiation treatment outcome. We propose a minimal mathematical model to study the effect of tumor heterogeneity and repair in different radiation treatment schedules. We perform stochastic and deterministic simulations to estimate model parameters using available experimental data. Our results suggest that gross tumor volume reduction is insufficient to control the disease if a fraction of radioresistant cells survives therapy. If cure cannot be achieved, protocols should balance volume reduction with minimal selection for radioresistant cells. We show that the most efficient treatment schedule is dependent on biology and model parameter values and, therefore, emphasize the need for careful tumor-specific model calibration before clinically actionable conclusions can be drawn.

Clinical implications of in silico mathematical modeling for glioblastoma: a critical review

Glioblastoma remains a clinical challenge in spite of years of extensive research. Novel approaches are needed in order to integrate the existing knowledge. This is the potential role of mathematical oncology. This paper reviews mathematical models on glioblastoma from the clinical doctor's point of view, with focus on 3D modeling approaches of radiation response of in vivo glioblastomas based on contemporary imaging techniques. As these models aim to provide a clinically useful tool in the era of personalized medicine, the integration of the latest advances in molecular and imaging science and in clinical practice by the in silico models is crucial for their clinical relevance. Our aim is to indicate areas of GBM research that have not yet been addressed by in silico models and to point out evidence that has come up from in silico experiments, which may be worth considering in the clinic. This review examines how close these models have come in predicting the outcome of treatment protocols and in shaping the future of radiotherapy treatments.

Applying a 4D Multiscale In Vivo Tumor Growth Model to the Exploration of Radiotherapy Scheduling: The Effects of Weekend Treatment Gaps and P53 Gene Status on the Response of Fast Growing Solid Tumors

The present paper aims at demonstrating clinically oriented applications of the multiscale four dimensional in vivo tumor growth simulation model previously developed by our research group. To this end the effect of weekend radiotherapy treatment gaps and p53 gene status on two virtual glioblastoma tumors differing only in p53 gene status is investigated in silico. Tumor response predictions concerning two rather extreme dose fractionation schedules (daily dose of 4.5 Gy administered in 3 equal fractions) namely HART (Hyperfractionated Accelerated Radiotherapy weekend less) 54 Gy and CHART (Continuous HART) 54 Gy are presented and compared. The model predictions suggest that, for the same p53 status, HART 54 Gy and CHART 54 Gy have almost the same long term effects on locoregional tumor control. However, no data have been located in the literature concerning a comparison of HART and CHART radiotherapy schedules for glioblastoma. As non small cell lung carcinoma (NSCLC) may also be a fast growing and radiosensitive tumor, a comparison of the model predictions with the outcome of clinical studies concerning the response of NSCLC to HART 54 Gy and CHART 54 Gy is made. The model predictions are in accordance with corresponding clinical observations, thus strengthening the potential of the model.

A Numerical Handling of the Boundary Conditions Imposed by the Skull on an Inhomogeneous Diffusion-Reaction Model of Glioblastoma Invasion Into the Brain: Clinical Validation Aspects

A novel explicit triscale reaction-diffusion numerical model of glioblastoma multiforme tumor growth is presented. The model incorporates the handling of Neumann boundary conditions imposed by the cranium and takes into account both the inhomogeneous nature of human brain and the complexity of the skull geometry. The finite-difference time-domain method is adopted. To demonstrate the workflow of a possible clinical validation procedure, a clinical case/scenario is addressed. A good agreement of the in silico calculated value of the doubling time (ie, the time for tumor volume to double) with the value of the same quantity based on tomographic imaging data has been observed. A theoretical exploration suggests that a rough but still quite informative value of the doubling time may be calculated based on a homogeneous brain model. The model could serve as the main component of a continuous mathematics-based glioblastoma oncosimulator aiming at supporting the clinician in the optimal patient-individualized design of treatment using the patient's multiscale data and experimenting in silico (ie, on the computer).

Comparison between model-predicted tumor oxygenation dynamics and vascular-/flow-related Doppler indices

PURPOSE

Mathematical modeling is a powerful and flexible method to investigate complex phenomena. It discloses the possibility of reproducing expensive as well as invasive experiments in a safe environment with limited costs. This makes it suitable to mimic tumor evolution and response to radiotherapy although the reliability of the results remains an issue. Complexity reduction is therefore a critical aspect in order to be able to compare model outcomes to clinical data. Among the factors affecting treatment efficacy, tumor oxygenation is known to play a key role in radiotherapy response. In this work, we aim at relating the oxygenation dynamics, predicted by a macroscale model trained on tumor volumetric data of uterine cervical cancer patients, to vascularization and blood flux indices assessed on Ultrasound Doppler images.

METHODS

We propose a macroscale model of tumor evolution based on three dynamics, namely active portion, necrotic portion, and oxygenation. The model parameters were assessed on the volume size of seven cervical cancer patients administered with 28 fractions of intensity modulated radiation therapy (IMRT) (1.8 Gy/fraction). For each patient, five Doppler ultrasound tests were acquired before, during, and after the treatment. The lesion was manually contoured by an expert physician using 4D View^® (General Electric Company - Fairfield, Connecticut, United States), which automatically provided the overall tumor volume size along with three vascularization and/or blood flow indices. Volume data only were fed to the model for training purpose, while the predicted oxygenation was compared a posteriori to the measured Doppler indices.

RESULTS

The model was able to fit the tumor volume evolution within 8% error (range: 3-8%). A strong correlation between the intrapatient longitudinal indices from Doppler measurements and oxygen predicted by the model (about 90% or above) was found in three cases. Two patients showed an average correlation value (50-70%) and the remaining two presented poor correlations. The latter patients were the ones featuring the smallest tumor reduction throughout the treatment, typical of hypoxic conditions. Moreover, the average oxygenation value predicted by the model was close to the average vascularization-flow index (average difference: 7%).

CONCLUSIONS

The results suggest that the modeled relation between tumor evolution and oxygen dynamics was reasonable enough to provide realistic oxygenation curves in five cases (correlation greater than 50%) out of seven. In case of nonresponsive tumors, the model failed in predicting the oxygenation trend while succeeded in reproducing the average oxygenation value according to the mean vascularization-flow index. Despite the need for deeper investigations, the outcomes of the present work support the hypothesis that a simple macroscale model of tumor response to radiotherapy is able to predict the tumor oxygenation. The possibility of an objective and quantitative validation on imaging data discloses the possibility to translate them as decision support tools in clinical practice and to move a step forward in the treatment personalization.

A radiobiological model of radiotherapy response and its correlation with prognostic imaging variables

Radiobiological models of tumour control probability (TCP) can be personalized using imaging data. We propose an extension to a voxel-level radiobiological TCP model in order to describe patient-specific differences and intra-tumour heterogeneity. In the proposed model, tumour shrinkage is described by means of a novel kinetic Monte Carlo method for inter-voxel cell migration and tumour deformation. The model captures the spatiotemporal evolution of the tumour at the voxel level, and is designed to take imaging data as input. To test the performance of the model, three image-derived variables found to be predictive of outcome in the literature have been identified and calculated using the model's own parameters. Simulating multiple tumours with different initial conditions makes it possible to perform an in silico study of the correlation of these variables with the dose for 50% tumour control ([Formula: see text]) calculated by the model. We find that the three simulated variables correlate with the calculated [Formula: see text]. In addition, we find that different variables have different levels of sensitivity to the spatial distribution of hypoxia within the tumour, as well as to the dynamics of the migration mechanism. Finally, based on our results, we observe that an adequate combination of the variables may potentially result in higher predictive power.

In Silico Oncology: Quantification of the In Vivo Antitumor Efficacy of Cisplatin-Based Doublet Therapy in Non-Small Cell Lung Cancer (NSCLC) through a Multiscale Mechanistic Model

The 5-year survival of non-small cell lung cancer patients can be as low as 1% in advanced stages. For patients with resectable disease, the successful choice of preoperative chemotherapy is critical to eliminate micrometastasis and improve operability. In silico experimentations can suggest the optimal treatment protocol for each patient based on their own multiscale data. A determinant for reliable predictions is the a priori estimation of the drugs’ cytotoxic efficacy on cancer cells for a given treatment. In the present work a mechanistic model of cancer response to treatment is applied for the estimation of a plausible value range of the cell killing efficacy of various cisplatin-based doublet regimens. Among others, the model incorporates the cancer related mechanism of uncontrolled proliferation, population heterogeneity, hypoxia and treatment resistance. The methodology is based on the provision of tumor volumetric data at two time points, before and after or during treatment. It takes into account the effect of tumor microenvironment and cell repopulation on treatment outcome. A thorough sensitivity analysis based on one-factor-at-a-time and latin hypercube sampling/partial rank correlation coefficient approaches has established the volume growth rate and the growth fraction at diagnosis as key features for more accurate estimates. The methodology is applied on the retrospective data of thirteen patients with non-small cell lung cancer who received cisplatin in combination with gemcitabine, vinorelbine or docetaxel in the neoadjuvant context. The selection of model input values has been guided by a comprehensive literature survey on cancer-specific proliferation kinetics. The latin hypercube sampling has been recruited to compensate for patient-specific uncertainties. Concluding, the present work provides a quantitative framework for the estimation of the in-vivo cell-killing ability of various chemotherapies. Correlation studies of such estimates with the molecular profile of patients could serve as a basis for reliable personalized predictions.

Less than 14% of medically treated patients with locally advanced and metastatic non-small cell lung cancer are expected to be alive 5 years after diagnosis. Standard therapeutic strategies include the administration of two drugs in combination, aiming at shrinking the tumor before surgery and improving overall survival. Knowing the sensitivity profile of each patient to different treatment strategies at diagnosis may help choose the most appropriate ones. We develop a methodology for the quantitative estimation of the cytotoxic efficacy of cisplatin-based doublets on cancer cells by applying a simulation model of cancer progression and response. The model incorporates the proliferation cycle, quiescence, differentiation and loss of tumor cells. We evaluate the effect of in vivo microenvironment of real tumors, as expressed by measurable tumor proliferation kinetics, such as how fast the tumor grows, the percentage of cells that are actively dividing, the resistance of stem cells, etc. on treatment outcome so as to derive more accurate estimates. A literature survey guides the selection of values. The methodology is applied to a real clinical dataset of patients. Correlation studies between the derived cytotoxicities and the patients’ molecular profile could lead to predictions of treatment response at the time of diagnosis.

Optimal treatment and stochastic modeling of heterogeneous tumors

In this work we review past articles that have mathematically studied cancer heterogeneity and the impact of this heterogeneity on the structure of optimal therapy. We look at past works on modeling how heterogeneous tumors respond to radiotherapy, and take a particularly close look at how the optimal radiotherapy schedule is modified by the presence of heterogeneity. In addition, we review past works on the study of optimal chemotherapy when dealing with heterogeneous tumors.

REVIEWERS

This article was reviewed by Thomas McDonald, David Axelrod, and Leonid Hanin.

Modeling the efficacy of the extent of surgical resection in the setting of radiation therapy for glioblastoma

Standard therapy for glioblastoma (GBM) includes maximal surgical resection and radiation therapy. While it is established that radiation therapy provides the greatest survival benefit of standard treatment modalities, the impact of the extent of surgical resection (EOR) on patient outcome remains highly controversial. While some studies describe no correlation between EOR and patient survival even up to total resection, others propose either qualitative (partial versus subtotal versus complete resection) or quantitative EOR thresholds, below which there is no correlation with survival. This work uses a mathematical model in the form of a reaction–diffusion partial differential equation to simulate tumor growth and treatment with radiation therapy and surgical resection based on tumor‐specific rates of diffusion and proliferation. Simulation of 36 tumors across a wide spectrum of diffusion and proliferation rates suggests that while partial or subtotal resections generally do not provide a survival advantage, complete resection significantly improves patient outcomes. Furthermore, our model predicts a tumor‐specific quantitative threshold below which EOR has no effect on patient survival and demonstrates that this threshold increases with tumor aggressiveness, particularly with the rate of proliferation. Thus, this model may serve as an aid for determining both when surgical resection is indicated as well as the surgical margins necessary to provide clinically significant improvements in patient survival. In addition, by assigning relative benefits to radiation and surgical resection based on tumor invasiveness and proliferation, this model confirms that (with the exception of the least aggressive tumors) the survival benefit of radiation therapy exceeds that of surgical resection.

Feasibility of four-dimensional preoperative simulation for elbow debridement arthroplasty

BACKGROUND

Recent advances in imaging modalities have enabled three-dimensional preoperative simulation. A four-dimensional preoperative simulation system would be useful for debridement arthroplasty of primary degenerative elbow osteoarthritis because it would be able to detect the impingement lesions.

METHODS

We developed a four-dimensional simulation system by adding the anatomical axis to the three-dimensional computed tomography scan data of the affected arm in one position. Eleven patients with primary degenerative elbow osteoarthritis were included. A "two rings" method was used to calculate the flexion-extension axis of the elbow by converting the surface of the trochlea and capitellum into two rings. A four-dimensional simulation movie was created and showed the optimal range of motion and the impingement area requiring excision. To evaluate the reliability of the flexion-extension axis, interobserver and intraobserver reliabilities regarding the assessment of bony overlap volumes were calculated twice for each patient by two authors. Patients were treated by open or arthroscopic debridement arthroplasties. Pre- and postoperative examinations included elbow range of motion measurement, and completion of the patient-rated questionnaire Hand20, Japanese Orthopaedic Association-Japan Elbow Society Elbow Function Score, and the Mayo Elbow Performance Score.

RESULTS

Measurement of the bony overlap volume showed an intraobserver intraclass correlation coefficient of 0.93 and 0.90, and an interobserver intraclass correlation coefficient of 0.94. The mean elbow flexion-extension arc significantly improved from 101° to 125°. The mean Hand20 score significantly improved from 52 to 22. The mean Japanese Orthopaedic Association-Japan Elbow Society Elbow Function Score significantly improved from 67 to 88. The mean Mayo Elbow Performance Score significantly improved from 71 to 91 at the final follow-up evaluation.

CONCLUSION

We showed that four-dimensional, preoperative simulation can be generated by adding the rotation axis to the one-position, three-dimensional computed tomography image of the affected arm. This method is feasible for elbow debridement arthroplasty.

Numerical simulation of vascular tumour growth under antiangiogenic treatment: addressing the paradigm of single-agent bevacizumab therapy with the use of experimental data

Background

Antiangiogenic agents have been recently added to the oncological armamentarium with bevacizumab probably being the most popular representative in current clinical practice. The elucidation of the mode of action of these agents is a prerequisite for personalized prediction of antiangiogenic treatment response and selection of patients who may benefit from this kind of therapy. To this end, having used as a basis a preexisting continuous vascular tumour growth model which addresses the targeted nature of antiangiogenic treatment, we present a paper characterized by the following three features. First, the integration of a two-compartmental bevacizumab specific pharmacokinetic module into the core of the aforementioned preexisting model. Second, its mathematical modification in order to reproduce the asymptotic behaviour of tumour volume in the theoretical case of a total destruction of tumour neovasculature. Third, the exploitation of a range of published animal datasets pertaining to antitumour efficacy of bevacizumab on various tumour types (breast, lung, head and neck, colon).

Results

Results for both the unperturbed growth and the treatment module reveal qualitative similarities with experimental observations establishing the biologically acceptable behaviour of the model. The dynamics of the untreated tumour has been studied via a parameter analysis, revealing the role of each relevant input parameter to tumour evolution. The combined effect of endogenous proangiogenic and antiangiogenic factors on the angiogenic potential of a tumour is also studied, in order to capture the dynamics of molecular competition between the two key-players of tumoural angiogenesis. The adopted methodology also allows accounting for the newly recognized direct antitumour effect of the specific agent.

Conclusions

Interesting observations have been made, suggesting a potential size-dependent tumour response to different treatment modalities and determining the relative timing of cytotoxic versus antiangiogenic agents administration. Insight into the comparative effectiveness of different antiangiogenic treatment strategies is revealed. The results of a series of in vivo experiments in mice bearing diverse types of tumours (breast, lung, head and neck, colon) and treated with bevacizumab are successfully reproduced, supporting thus the validity of the underlying model.

Reviewers

This article was reviewed by L. Hanin, T. Radivoyevitch and L. Edler.

Integrating Multiscale Modeling with Drug Effects for Cancer Treatment

In this paper, we review multiscale modeling for cancer treatment with the incorporation of drug effects from an applied system's pharmacology perspective. Both the classical pharmacology and systems biology are inherently quantitative; however, systems biology focuses more on networks and multi factorial controls over biological processes rather than on drugs and targets in isolation, whereas systems pharmacology has a strong focus on studying drugs with regard to the pharmacokinetic (PK) and pharmacodynamic (PD) relations accompanying drug interactions with multiscale physiology as well as the prediction of dosage-exposure responses and economic potentials of drugs. Thus, it requires multiscale methods to address the need for integrating models from the molecular levels to the cellular, tissue, and organism levels. It is a common belief that tumorigenesis and tumor growth can be best understood and tackled by employing and integrating a multifaceted approach that includes in vivo and in vitro experiments, in silico models, multiscale tumor modeling, continuous/discrete modeling, agent-based modeling, and multiscale modeling with PK/PD drug effect inputs. We provide an example application of multiscale modeling employing stochastic hybrid system for a colon cancer cell line HCT-116 with the application of Lapatinib drug. It is observed that the simulation results are similar to those observed from the setup of the wet-lab experiments at the Translational Genomics Research Institute.

Spatio-Temporal Dynamics of Hypoxia during Radiotherapy

Tumour hypoxia plays a pivotal role in cancer therapy for most therapeutic approaches from radiotherapy to immunotherapy. The detailed and accurate knowledge of the oxygen distribution in a tumour is necessary in order to determine the right treatment strategy. Still, due to the limited spatial and temporal resolution of imaging methods as well as lacking fundamental understanding of internal oxygenation dynamics in tumours, the precise oxygen distribution map is rarely available for treatment planing. We employ an agent-based in silico tumour spheroid model in order to study the complex, localized and fast oxygen dynamics in tumour micro-regions which are induced by radiotherapy. A lattice-free, 3D, agent-based approach for cell representation is coupled with a high-resolution diffusion solver that includes a tissue density-dependent diffusion coefficient. This allows us to assess the space- and time-resolved reoxygenation response of a small subvolume of tumour tissue in response to radiotherapy. In response to irradiation the tumour nodule exhibits characteristic reoxygenation and re-depletion dynamics which we resolve with high spatio-temporal resolution. The reoxygenation follows specific timings, which should be respected in treatment in order to maximise the use of the oxygen enhancement effects. Oxygen dynamics within the tumour create windows of opportunity for the use of adjuvant chemotherapeutica and hypoxia-activated drugs. Overall, we show that by using modelling it is possible to follow the oxygenation dynamics beyond common resolution limits and predict beneficial strategies for therapy and in vitro verification. Models of cell cycle and oxygen dynamics in tumours should in the future be combined with imaging techniques, to allow for a systematic experimental study of possible improved schedules and to ultimately extend the reach of oxygenation monitoring available in clinical treatment.

A voxel-based multiscale model to simulate the radiation response of hypoxic tumors

PURPOSE

In radiotherapy, it is important to predict the response of tumors to irradiation prior to the treatment. This is especially important for hypoxic tumors, which are known to be highly radioresistant. Mathematical modeling based on the dose distribution, biological parameters, and medical images may help to improve this prediction and to optimize the treatment plan.

METHODS

A voxel-based multiscale tumor response model for simulating the radiation response of hypoxic tumors was developed. It considers viable and dead tumor cells, capillary and normal cells, as well as the most relevant biological processes such as (i) proliferation of tumor cells, (ii) hypoxia-induced angiogenesis, (iii) spatial exchange of cells leading to tumor growth, (iv) oxygen-dependent cell survival after irradiation, (v) resorption of dead cells, and (vi) spatial exchange of cells leading to tumor shrinkage. Oxygenation is described on a microscopic scale using a previously published tumor oxygenation model, which calculates the oxygen distribution for each voxel using the vascular fraction as the most important input parameter. To demonstrate the capabilities of the model, the dependence of the oxygen distribution on tumor growth and radiation-induced shrinkage is investigated. In addition, the impact of three different reoxygenation processes is compared and tumor control probability (TCP) curves for a squamous cells carcinoma of the head and neck (HNSSC) are simulated under normoxic and hypoxic conditions.

RESULTS

The model describes the spatiotemporal behavior of the tumor on three different scales: (i) on the macroscopic scale, it describes tumor growth and shrinkage during radiation treatment, (ii) on a mesoscopic scale, it provides the cell density and vascular fraction for each voxel, and (iii) on the microscopic scale, the oxygen distribution may be obtained in terms of oxygen histograms. With increasing tumor size, the simulated tumors develop a hypoxic core. Within the model, tumor shrinkage was found to be significantly more important for reoxygenation than angiogenesis or decreased oxygen consumption due to an increased fraction of dead cells. In the studied HNSSC-case, the TCD50 values (dose at 50% TCP) decreased from 71.0 Gy under hypoxic to 53.6 Gy under the oxic condition.

CONCLUSIONS

The results obtained with the developed multiscale model are in accordance with expectations based on radiobiological principles and clinical experience. As the model is voxel-based, radiological imaging methods may help to provide the required 3D-characterization of the tumor prior to irradiation. For clinical application, the model has to be further validated with experimental and clinical data. If this is achieved, the model may be used to optimize fractionation schedules and dose distributions for the treatment of hypoxic tumors.

In silico oncology: exploiting clinical studies to clinically adapt and validate multiscale oncosimulators

This paper presents a brief outline of the notion and the system of oncosimulator in conjunction with a high level description of the basics of its core multiscale model simulating clinical tumor response to treatment. The exemplary case of lung cancer preoperatively treated with a combination of chemotherapeutic agents is considered. The core oncosimulator model is based on a primarily top-down, discrete entity - discrete event multiscale simulation approach. The critical process of clinical adaptation of the model by exploiting sets of multiscale data originating from clinical studies/trials is also outlined. Concrete clinical adaptation results are presented. The adaptation process also conveys important aspects of the planned clinical validation procedure since the same type of multiscale data - although not the same data itself- is to be used for clinical validation. By having exploited actual clinical data in conjunction with plausible literature-based values of certain model parameters, a realistic tumor dynamics behavior has been demonstrated. The latter supports the potential of the specific oncosimulator to serve as a personalized treatment optimizer following an eventually successful completion of the clinical adaptation and validation process.

Optimization of radiation dosing schedules for proneural glioblastoma

Glioblastomas are the most aggressive primary brain tumor. Despite treatment with surgery, radiation and chemotherapy, these tumors remain uncurable and few significant increases in survival have been observed over the last half-century. We recently employed a combined theoretical and experimental approach to predict the effectiveness of radiation administration schedules, identifying two schedules that led to superior survival in a mouse model of the disease (Leder et al., Cell 156(3):603-616, 2014). Here we extended this approach to consider fractionated schedules to best minimize toxicity arising in early- and late-responding tissues. To this end, we decomposed the problem into two separate solvable optimization tasks: (i) optimization of the amount of radiation per dose, and (ii) optimization of the amount of time that passes between radiation doses. To ensure clinical applicability, we then considered the impact of clinical operating hours by incorporating time constraints consistent with operational schedules of the radiology clinic. We found that there was no significant loss incurred by restricting dosage to an 8:00 a.m. to 5:00 p.m. window. Our flexible approach is also applicable to other tumor types treated with radiotherapy.

Kinetic Models for Predicting Cervical Cancer Response to Radiation Therapy on Individual Basis Using Tumor Regression Measured In Vivo With Volumetric Imaging

This article describes a macroscopic mathematical modeling approach to capture the interplay between solid tumor evolution and cell damage during radiotherapy. Volume regression profiles of 15 patients with uterine cervical cancer were reconstructed from serial cone-beam computed tomography data sets, acquired for image-guided radiotherapy, and used for model parameter learning by means of a genetic-based optimization. Patients, diagnosed with either squamous cell carcinoma or adenocarcinoma, underwent different treatment modalities (image-guided radiotherapy and image-guided chemo-radiotherapy). The mean volume at the beginning of radiotherapy and the end of radiotherapy was on average 23.7 cm(3) (range: 12.7-44.4 cm(3)) and 8.6 cm(3) (range: 3.6-17.1 cm(3)), respectively. Two different tumor dynamics were taken into account in the model: the viable (active) and the necrotic cancer cells. However, according to the results of a preliminary volume regression analysis, we assumed a short dead cell resolving time and the model was simplified to the active tumor volume. Model learning was performed both on the complete patient cohort (cohort-based model learning) and on each single patient (patient-specific model learning). The fitting results (mean error: ∼ 16% and ∼ 6% for the cohort-based model and patient-specific model, respectively) highlighted the model ability to quantitatively reproduce tumor regression. Volume prediction errors of about 18% on average were obtained using cohort-based model computed on all but 1 patient at a time (leave-one-out technique). Finally, a sensitivity analysis was performed and the data uncertainty effects evaluated by simulating an average volume perturbation of about 1.5 cm(3) obtaining an error increase within 0.2%. In conclusion, we showed that simple time-continuous models can represent tumor regression curves both on a patient cohort and patient-specific basis; this discloses the opportunity in the future to exploit such models to predict how changes in the treatment schedule (number of fractions, doses, intervals among fractions) might affect the tumor regression on an individual basis.

Optimal weekly scheduling in fractionated radiotherapy: effect of an upper bound on the dose fraction size

This work concerns the optimization of the dose fractionation for cancer radiotherapy schedules of the kind one fraction/day, five fractions/week, assuming a fixed overall treatment time. Constraints are set to limit the radiation damages to surrounding normal tissues, as well as the daily fraction size. The response to radiation of tumour and normal tissues is represented by the classical LQ model, including the exponential repopulation term. We provide a framework to analytically determine the optimal weekly scheme of radiation doses as a function of the tumour type, the fraction upper bound and the normal tissue parameters. For a comparison with the literature, we present some numerical examples of optimal treatment schedules for specific tumour types.

Identification of crucial parameters in a mathematical multiscale model of glioblastoma growth

Glioblastomas are highly malignant brain tumours. Mathematical models and their analysis provide a tool to support the understanding of the development of these tumours as well as the design of more effective treatment strategies. We have previously developed a multiscale model of glioblastoma progression that covers processes on the cellular and molecular scale. Here, we present a novel nutrient-dependent multiscale sensitivity analysis of this model that helps to identify those reaction parameters of the molecular interaction network that influence the tumour progression on the cellular scale the most. In particular, those parameters are identified that essentially determine tumour expansion and could be therefore used as potential therapy targets. As indicators for the success of a potential therapy target, a deceleration of the tumour expansion and a reduction of the tumour volume are employed. From the results, it can be concluded that no single parameter variation results in a less aggressive tumour. However, it can be shown that a few combined perturbations of two systematically selected parameters cause a slow-down of the tumour expansion velocity accompanied with a decrease of the tumour volume. Those parameters are primarily linked to the reactions that involve the microRNA-451 and the thereof regulated protein MO25.

Mathematical modeling of PDGF-driven glioblastoma reveals optimized radiation dosing schedules

Glioblastomas (GBMs) are the most common and malignant primary brain tumors and are aggressively treated with surgery, chemotherapy, and radiotherapy. Despite this treatment, recurrence is inevitable and survival has improved minimally over the last 50 years. Recent studies have suggested that GBMs exhibit both heterogeneity and instability of differentiation states and varying sensitivities of these states to radiation. Here, we employed an iterative combined theoretical and experimental strategy that takes into account tumor cellular heterogeneity and dynamically acquired radioresistance to predict the effectiveness of different radiation schedules. Using this model, we identified two delivery schedules predicted to significantly improve efficacy by taking advantage of the dynamic instability of radioresistance. These schedules led to superior survival in mice. Our interdisciplinary approach may also be applicable to other human cancer types treated with radiotherapy and, hence, may lay the foundation for significantly increasing the effectiveness of a mainstay of oncologic therapy. PAPERCLIP:

A mathematical model of tumor volume changes during radiotherapy

PURPOSE

To develop a clinically viable mathematical model that quantitatively predicts tumor volume change during radiotherapy in order to provide treatment response assessment for prognosis, treatment plan optimization, and adaptation.

METHOD AND MATERIALS

The correction factors containing hypoxia, DNA single strand breaks, potentially lethal damage, and other factors were used to develop an improved cell survival model based on the popular linear-quadratic model of cell survival in radiotherapy. The four-level cell population model proposed by Chvetsov et al. was further simplified by removing the initial hypoxic fraction and reoxygenation parameter, which are hard to obtain in routine clinics, such that an easy-to-use model can be developed for clinical applications. The new model was validated with data of nine lung and cervical cancer patients.

RESULTS

Out of the nine cases, the new model can predict tumor volume change in six cases with a correlation index R (2) greater than 0.9 and the rest of three with R (2) greater than 0.85.

CONCLUSION

Based on a four-level cell population model, a more practical and simplified cell survival curve was proposed to model the tumor volume changes during radiotherapy. Validation study with patient data demonstrated feasibility and clinical usefulness of the new model in predicting tumor volume change in radiotherapy.

An architecture for integrating cancer model repositories

The TUMOR project aims at developing a European clinically oriented semantic-layered cancer digital model repository from existing EU projects that will be interoperable with the US grid-enabled semantic-layered digital model repository platform at CViT.org (Center for the Development of a Virtual Tumor, Massachusetts General Hospital (MGH), Boston, USA) which is NIH/NCI-caGRID compatible. In this paper we describe the modular and federated architecture of TUMOR that effectively addresses model integration, interoperability, and security related issues.

Modelling the interplay between hypoxia and proliferation in radiotherapy tumour response

A tumour control probability computational model for fractionated radiotherapy was developed, with the goal of incorporating the fundamental interplay between hypoxia and proliferation, including reoxygenation over a course of radiotherapy. The fundamental idea is that the local delivery of oxygen and glucose limits the amount of proliferation and metabolically-supported cell survival a tumour sub-volume can support. The model has three compartments: a proliferating compartment of cells receiving oxygen and glucose; an intermediate, metabolically-active compartment receiving glucose; and a highly hypoxic compartment of starving cells. Following the post-mitotic cell death of proliferating cells, intermediate cells move into the proliferative compartment and hypoxic cells move into the intermediate compartment. A key advantage of the proposed model is that the initial compartmental cell distribution is uniquely determined from the assumed local growth fraction (GF) and volume doubling time (TD) values. Varying initial cell state distributions, based on the local (voxel) GF and TD, were simulated. Tumour response was simulated for head and neck squamous cell carcinoma using relevant parameter values based on published sources. The tumour dose required to achieve a 50% local control rate (TCD50) was found for various GFs and TD's, and the effect of fraction size on TCD50 was also evaluated. Due to the advantage of reoxygenation over a course of radiotherapy, conventional fraction sizes (2-2.4 Gy fx(-1)) were predicted to result in smaller TCD50's than larger fraction sizes (4-5 Gy fx(-1)) for a 10 cc tumour with GFs of around 0.15. The time to eliminate hypoxic cells (the reoxygenation time) was estimated for a given GF and decreased as GF increased. The extra dose required to overcome accelerated stem cell accumulation in longer treatment schedules was estimated to be 0.68 Gy/day (in EQD26.6), similar to published values derived from clinical data. The model predicts, for a 2 Gy/weekday fractionation, that increased initial proliferation (high GF) should, surprisingly, lead to moderately higher local control values. Tumour hypoxia is predicted to increase the required dose for local control by approximately 30%. Predicted tumour regression patterns are consistent with clinical observations. This simple yet flexible model shows how the local competition for chemical resources might impact local control rates under varying fractionation conditions.

Dealing with diversity in computational cancer modeling

This paper discusses the need for interconnecting computational cancer models from different sources and scales within clinically relevant scenarios to increase the accuracy of the models and speed up their clinical adaptation, validation, and eventual translation. We briefly review current interoperability efforts drawing upon our experiences with the development of in silico models for predictive oncology within a number of European Commission Virtual Physiological Human initiative projects on cancer. A clinically relevant scenario, addressing brain tumor modeling that illustrates the need for coupling models from different sources and levels of complexity, is described. General approaches to enabling interoperability using XML-based markup languages for biological modeling are reviewed, concluding with a discussion on efforts towards developing cancer-specific XML markup to couple multiple component models for predictive in silico oncology.

In silico modelling of tumour margin diffusion and infiltration: review of current status

As a result of advanced treatment techniques, requiring precise target definitions, a need for more accurate delineation of the Clinical Target Volume (CTV) has arisen. Mathematical modelling is found to be a powerful tool to provide fairly accurate predictions for the Microscopic Extension (ME) of a tumour to be incorporated in a CTV. In general terms, biomathematical models based on a sequence of observations or development of a hypothesis assume some links between biological mechanisms involved in cancer development and progression to provide quantitative or qualitative measures of tumour behaviour as well as tumour response to treatment. Generally, two approaches are taken: deterministic and stochastic modelling. In this paper, recent mathematical models, including deterministic and stochastic methods, are reviewed and critically compared. It is concluded that stochastic models are more promising to provide a realistic description of cancer tumour behaviour due to being intrinsically probabilistic as well as discrete, which enables incorporation of patient-specific biomedical data such as tumour heterogeneity and anatomical boundaries.

The HYP-RT hypoxic tumour radiotherapy algorithm and accelerated repopulation dose per fraction study

The HYP-RT model simulates hypoxic tumour growth for head and neck cancer as well as radiotherapy and the effects of accelerated repopulation and reoxygenation. This report outlines algorithm design, parameterisation and the impact of accelerated repopulation on the increase in dose/fraction needed to control the extra cell propagation during accelerated repopulation. Cell kill probabilities are based on Linear Quadratic theory, with oxygenation levels and proliferative capacity influencing cell death. Hypoxia is modelled through oxygen level allocation based on pO₂ histograms. Accelerated repopulation is modelled by increasing the stem cell symmetrical division probability, while the process of reoxygenation utilises randomised pO₂ increments to the cell population after each treatment fraction. Propagation of 10⁸ tumour cells requires 5–30 minutes. Controlling the extra cell growth induced by accelerated repopulation requires a dose/fraction increase of 0.5–1.0 Gy, in agreement with published reports. The average reoxygenation pO₂ increment of 3 mmHg per fraction results in full tumour reoxygenation after shrinkage to approximately 1 mm. HYP-RT is a computationally efficient model simulating tumour growth and radiotherapy, incorporating accelerated repopulation and reoxygenation. It may be used to explore cell kill outcomes during radiotherapy while varying key radiobiological and tumour specific parameters, such as the degree of hypoxia.

In-depth analysis and evaluation of diffusive glioma models

Glioma is one of the most aggressive types of brain tumor. Several mathematical models have been developed during the past two decades, toward simulating the mechanisms that govern the development of glioma. The most common models use the diffusion-reaction equation (DRE) for simulating the spatiotemporal variation of tumor cell concentration. Nevertheless, despite the applications presented, there has been little work on studying the details of the mathematical solution and implementation of the 3-D diffusion model and presenting a qualitative analysis of the algorithmic results. This paper presents a complete mathematical framework on the solution of the DRE using different numerical schemes. This framework takes into account all characteristics of the latest models, such as brain tissue heterogeneity, anisotropic tumor cell migration, chemotherapy, and resection modeling. The different numerical schemes presented have been evaluated based upon the degree to which the DRE exact solution is approximated. Experiments have been conducted both on real datasets and a test case for which there is a known algebraic expression of the solution. Thus, it is possible to calculate the accuracy of the different models.

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

OBJECTIVE

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.

METHOD

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.

RESULTS

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.

CONCLUSION

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.

Clinically driven design of multi-scale cancer models: the ContraCancrum project paradigm

The challenge of modelling cancer presents a major opportunity to improve our ability to reduce mortality from malignant neoplasms, improve treatments and meet the demands associated with the individualization of care needs. This is the central motivation behind the ContraCancrum project. By developing integrated multi-scale cancer models, ContraCancrum is expected to contribute to the advancement of in silico oncology through the optimization of cancer treatment in the patient-individualized context by simulating the response to various therapeutic regimens. The aim of the present paper is to describe a novel paradigm for designing clinically driven multi-scale cancer modelling by bringing together basic science and information technology modules. In addition, the integration of the multi-scale tumour modelling components has led to novel concepts of personalized clinical decision support in the context of predictive oncology, as is also discussed in the paper. Since clinical adaptation is an inelastic prerequisite, a long-term clinical adaptation procedure of the models has been initiated for two tumour types, namely non-small cell lung cancer and glioblastoma multiforme; its current status is briefly summarized.

Evaluation framework for the multilevel macroscopic models of solid tumor growth in the glioma case

This paper investigates the applicability of multilevel macroscopic models for simulating solid tumor growth in the invasive glioblastoma multiforme (GBM) case. The continuum case approach tumor model based on the diffusion reaction equation is evaluated on a pre-segmented tomographic atlas where all tissue properties are known a priori. The atlas is further registered on a real clinical case where the tumor invasion status is gauged in two successive points in time. Based on the latter, the model attempts to fully replicate tumor growth taking into account tissue based properties as identified from the atlas template. The whole process is performed on a clinical platform specially designed to facilitate precise identification and delineation of tumors of large number of 3D tomographic datasets by an expert clinician. The promising results presented encourage the potential clinical applicability of the proposed model in the glioma case and identify crucial points and direction of further model refinement and research.

Exploiting clinical trial data drastically narrows the window of possible solutions to the problem of clinical adaptation of a multiscale cancer model

The development of computational models for simulating tumor growth and response to treatment has gained significant momentum during the last few decades. At the dawn of the era of personalized medicine, providing insight into complex mechanisms involved in cancer and contributing to patient-specific therapy optimization constitute particularly inspiring pursuits. The in silico oncology community is facing the great challenge of effectively translating simulation models into clinical practice, which presupposes a thorough sensitivity analysis, adaptation and validation process based on real clinical data. In this paper, the behavior of a clinically-oriented, multiscale model of solid tumor response to chemotherapy is investigated, using the paradigm of nephroblastoma response to preoperative chemotherapy in the context of the SIOP/GPOH clinical trial. A sorting of the model's parameters according to the magnitude of their effect on the output has unveiled the relative importance of the corresponding biological mechanisms; major impact on the result of therapy is credited to the oxygenation and nutrient availability status of the tumor and the balance between the symmetric and asymmetric modes of stem cell division. The effect of a number of parameter combinations on the extent of chemotherapy-induced tumor shrinkage and on the tumor's growth rate are discussed. A real clinical case of nephroblastoma has served as a proof of principle study case, demonstrating the basics of an ongoing clinical adaptation and validation process. By using clinical data in conjunction with plausible values of model parameters, an excellent fit of the model to the available medical data of the selected nephroblastoma case has been achieved, in terms of both volume reduction and histological constitution of the tumor. In this context, the exploitation of multiscale clinical data drastically narrows the window of possible solutions to the clinical adaptation problem.

Comparing finite elements and finite differences for developing diffusive models of glioma growth

Glioma is the most aggressive type of brain tumor. Several mathematical models have been developed during the last two decades, towards simulating the mechanisms that govern the development of glioma. The most common models use the diffusion-reaction equation (DRE) for simulating the spatiotemporal variation of tumor cell concentration. The proposed diffusive models have mainly used finite differences (FDs) or finite elements (FEs) for the approximation of the solution of the partial differential DRE. This paper presents experimental results on the comparison of the FEs and FDs, especially focused on the glioma model case. It is studied how the different meshes of brain can affect computational consistency, simulation time and efficiency of the model. The experiments have been studied on a test case, for which there is a known algebraic expression of the solution. Thus, it is possible to calculate the error that the different models yield.

Computer simulation of tumour control probabilities after irradiation for varying intrinsic radio-sensitivity using a single cell based model

BACKGROUND

Currently, optimisation of the dose distribution and clinical acceptance are almost entirely based on the physical dose distribution and tumour control probability modelling is far from being routinely used as objective in treatment planning. For future individualised radiotherapeutic strategies, a reliable patient specific simulation model, taking into account customised tumour features, is needed to predict and improve treatment outcome.

MATERIALS AND METHODS

To approach these demands, a single cell and Monte-Carlo based model was developed, which enables three-dimensional tumour growth and radiation response simulation. Tumour cells were characterised by cell-associated features such as age, intrinsic radio-sensitivity, proliferation ability, and oxygenation status, while capillary cells were considered as sources of a radial-dependent oxygen profile. Response to radiation was simulated by the linear-quadratic model, taking into account the lower radio-sensitivity of poorly oxygenated tumour cells.

RESULTS

The present study shows the influence of the model components and demonstrates the impact of the intra- and inter-tumoural radio-sensitivity heterogeneity on the treatment response.

CONCLUSION

The simulation model adequately delineates the importance of the above described selected parameters on tumour control probability, providing an insight into the interplay of different physical and biological parameters, and its relevance for an individual tumour response.