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.

A study on the predictability of acute lymphoblastic leukaemia response to treatment using a hybrid oncosimulator.

Efficient use of Virtual Physiological Human (VPH)-type models for personalized treatment response prediction purposes requires a precise model parameterization. In the case where the available personalized data are not sufficient to fully determine the parameter values, an appropriate prediction task may be followed. This study, a hybrid combination of computational optimization and machine learning methods with an already developed mechanistic model called the acute lymphoblastic leukaemia (ALL) Oncosimulator which simulates ALL progression and treatment response is presented. These methods are used in order for the parameters of the model to be estimated for retrospective cases and to be predicted for prospective ones. The parameter value prediction is based on a regression model trained on retrospective cases. The proposed Hybrid ALL Oncosimulator system has been evaluated when predicting the pre-phase treatment outcome in ALL. This has been correctly achieved for a significant percentage of patient cases tested (approx. 70% of patients). Moreover, the system is capable of denying the classification of cases for which the results are not trustworthy enough. In that case, potentially misleading predictions for a number of patients are avoided, while the classification accuracy for the remaining patient cases further increases. The results obtained are particularly encouraging regarding the soundness of the proposed methodologies and their relevance to the process of achieving clinical applicability of the proposed Hybrid ALL Oncosimulator system and VPH models in general.

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.

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).

A Modular Repository-based Infrastructure for Simulation Model Storage and Execution Support in the Context of In Silico Oncology and In Silico Medicine.

The plethora of available disease prediction models and the ongoing process of their application into clinical practice - following their clinical validation - have created new needs regarding their efficient handling and exploitation. Consolidation of software implementations, descriptive information, and supportive tools in a single place, offering persistent storage as well as proper management of execution results, is a priority, especially with respect to the needs of large healthcare providers. At the same time, modelers should be able to access these storage facilities under special rights, in order to upgrade and maintain their work. In addition, the end users should be provided with all the necessary interfaces for model execution and effortless result retrieval. We therefore propose a software infrastructure, based on a tool, model and data repository that handles the storage of models and pertinent execution-related data, along with functionalities for execution management, communication with third-party applications, user-friendly interfaces to access and use the infrastructure with minimal effort and basic security features.

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.

Differentiation resistance through altered retinoblastoma protein function in acute lymphoblastic leukemia: in silico modeling of the deregulations in the G1/S restriction point pathway.

BACKGROUND: As in many cancer types, the G1/S restriction point (RP) is deregulated in Acute Lymphoblastic Leukemia (ALL). Hyper-phosphorylated retinoblastoma protein (hyper-pRb) is found in high levels in ALL cells. Nevertheless, the ALL lymphocyte proliferation rate for the average patient is surprisingly low compared to its normal counterpart of the same maturation level. Additionally, as stated in literature, ALL cells possibly reside at or beyond the RP which is located in the late-G1 phase. This state may favor their differentiation resistant phenotype. A major phenomenon contributing to this fact is thought to be the observed limited redundancy in the phosphorylation of retinoblastoma protein (pRb) by the various Cyclin Dependent Kinases (Cdks). The latter may result in partial loss of pRb functions despite hyper-phosphorylation. RESULTS: To test this hypothesis, an in silico model aiming at simulating the biochemical regulation of the RP in ALL is introduced. By exploiting experimental findings derived from leukemic cells and following a semi-quantitative calibration procedure, the model has been shown to satisfactorily reproduce such a behavior for the RP pathway. At the same time, the calibrated model has been proved to be in agreement with the observed variation in the ALL cell cycle duration. CONCLUSIONS: The proposed model aims to contribute to a better understanding of the complex phenomena governing the leukemic cell cycle. At the same time it constitutes a significant first step in the creation of a personalized proliferation rate predictor that can be used in the context of multiscale cancer modeling. Such an approach is expected to play an important role in the refinement and the advancement of mechanistic modeling of ALL in the context of the emergent and promising scientific domains of In Silico Oncology and more generally In Silico Medicine.

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.

Towards in silico oncology: adapting a four dimensional nephroblastoma treatment model to a clinical trial case based on multi-method sensitivity analysis.

In the past decades a great progress in cancer research has been made although medical treatment is still widely based on empirically established protocols which have many limitations. Computational models address such limitations by providing insight into the complex biological mechanisms of tumor progression. A set of clinically-oriented, multiscale models of solid tumor dynamics has been developed by the In Silico Oncology Group (ISOG), Institute of Communication and Computer Systems (ICCS)-National Technical University of Athens (NTUA) to study cancer growth and response to treatment. Within this context using certain representative parameter values, tumor growth and response have been modeled under a cancer preoperative chemotherapy protocol in the framework of the SIOP 2001/GPOH clinical trial. A thorough cross-method sensitivity analysis of the model has been performed. Based on the sensitivity analysis results, a reasonable adaptation of the values of the model parameters to a real clinical case of bilateral nephroblastomatosis has been achieved. The analysis presented supports the potential of the model for the study and eventually the future design of personalized treatment schemes and/or schedules using the data obtained from in vitro experiments and clinical studies.

Coupling biomechanics to a cellular level model: an approach to patient-specific image driven multi-scale and multi-physics tumor simulation.

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

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.

Introduction of hypermatrix and operator notation into a discrete mathematics simulation model of malignant tumour response to therapeutic schemes in vivo. Some operator properties.

The tremendous rate of accumulation of experimental and clinical knowledge pertaining to cancer dictates the development of a theoretical framework for the meaningful integration of such knowledge at all levels of biocomplexity. In this context our research group has developed and partly validated a number of spatiotemporal simulation models of in vivo tumour growth and in particular tumour response to several therapeutic schemes. Most of the modeling modules have been based on discrete mathematics and therefore have been formulated in terms of rather complex algorithms (e.g. in pseudocode and actual computer code). However, such lengthy algorithmic descriptions, although sufficient from the mathematical point of view, may render it difficult for an interested reader to readily identify the sequence of the very basic simulation operations that lie at the heart of the entire model. In order to both alleviate this problem and at the same time provide a bridge to symbolic mathematics, we propose the introduction of the notion of hypermatrix in conjunction with that of a discrete operator into the already developed models. Using a radiotherapy response simulation example we demonstrate how the entire model can be considered as the sequential application of a number of discrete operators to a hypermatrix corresponding to the dynamics of the anatomic area of interest. Subsequently, we investigate the operators' commutativity and outline the "summarize and jump" strategy aiming at efficiently and realistically address multilevel biological problems such as cancer. In order to clarify the actual effect of the composite discrete operator we present further simulation results which are in agreement with the outcome of the clinical study RTOG 83-02, thus strengthening the reliability of the model developed.

The "Oncosimulator": a multilevel, clinically oriented simulation system of tumor growth and organism response to therapeutic schemes. Towards the clinical evaluation of in silico oncology.

The "Oncosimulator" is at the same time a concept of multilevel integrative cancer and (treatment affected) normal tissue biology, an algorithmic construct and a software tool which aims at supporting the clinician in the process of optimizing cancer treatment on the patient individualized basis. Additionally it is a platform for better understanding and exploring the natural phenomenon of cancer as well as training doctors and interested patients alike. In order to achieve all of these goals it has to undergo a thorough clinical optimization and validation process. This is one of the goals of the European Commission funded integrated project "ACGT: Advancing Clinicogenomic Trials on Cancer". Nephroblastoma (Wilms' tumor) and breast cancer have been selected to serve as two paradigms to clinically specify and evaluate the "Oncosimulator" as well as the emerging domain of in silico oncology.

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 spatiotemporal, patient individualized simulation model of solid tumor response to chemotherapy in vivo: the paradigm of glioblastoma multiforme treated by temozolomide.

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.

A computer simulation of in vivo tumour growth and response to radiotherapy: new algorithms and parametric results.

The aim of this paper is to present the newest algorithms and simulation results of a computer model of in vivo tumour growth and response to radiotherapy. The new algorithms are analytically presented. A set of parametric simulations has been performed with special emphasis on the influence of the genetic profile of a tumour on its radiosensitivity. The results of the simulation procedure are three-dimensionally visualized and critically surveyed. The long-term goal of this work is twofold: the development of a computational tool for getting insight into cancer biology and the development of a patient-specific decision support system.

Simulating chemotherapeutic schemes in the individualized treatment context: the paradigm of glioblastoma multiforme treated by temozolomide in vivo.

A novel patient individualized, spatiotemporal Monte Carlo simulation model of tumor response to chemotherapeutic schemes in vivo is presented. Treatment of glioblastoma multiforme by temozolomide is considered as a paradigm. The model is based on the patient's imaging, histopathologic and genetic data. 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, apoptosis, etc.) in conjunction with pharmacokinetics and pharmacodynamics information are applied. A good qualitative agreement of the model's predictions with clinical experience supports the applicability of the approach to chemotherapy optimization.

Simulating growth dynamics and radiation response of avascular tumour spheroids-model validation in the case of an EMT6/Ro multicellular spheroid.

The goal of this paper is to provide both the basic scientist and the clinician with an advanced computational tool for performing in silico experiments aiming at supporting the process of biological optimisation of radiation therapy. Improved understanding and description of malignant tumour dynamics is an additional intermediate objective. To this end an advanced three-dimensional (3D) Monte-Carlo simulation model of both the avascular development of multicellular tumour spheroids and their response to radiation therapy is presented. The model is based upon a number of fundamental biological principles such as the transition between the cell cycle phases, the diffusion of oxygen and nutrients and the cell survival probabilities following irradiation. Efficient algorithms describing tumour expansion and shrinkage are proposed and applied. The output of the biosimulation model is introduced into the (3D) visualisation package AVS-Express, which performs the visualisation of both the external surface and the internal structure of the dynamically evolving tumour based on volume or surface rendering techniques. Both the numerical stability and the statistical behaviour of the simulation model have been studied and evaluated for the case of EMT6/Ro spheroids. Predicted histological structure and tumour growth rates have been shown to be in agreement with published experimental data. Furthermore, the underlying structure of the tumour spheroid as well as its response to irradiation satisfactorily agrees with laboratory experience.