Eleni A. Kolokotroni

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.

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 models 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 tumors 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.

The Technologically Integrated Oncosimulator: Combining Multiscale Cancer Modeling With Information Technology in the In Silico Oncology Context

This paper outlines the major components and function of the technologically integrated oncosimulator developed primarily within the Advancing Clinico Genomic Trials on Cancer (ACGT) project. The Oncosimulator is defined as an information technology system simulating in vivo tumor response to therapeutic modalities within the clinical trial context. Chemotherapy in the neoadjuvant setting, according to two real clinical trials concerning nephroblastoma and breast cancer, has been considered. The spatiotemporal simulation module embedded in the Oncosimulator is based on the multiscale, predominantly top-down, discrete entity-discrete event cancer simulation technique developed by the In Silico Oncology Group, National Technical University of Athens. The technology modules include multiscale data handling, image processing, invocation of code execution via a spreadsheet-inspired environment portal, execution of the code on the grid, and the visualization of the predictions. A refining scenario for the eventual coupling of the oncosimulator with immunological models is also presented. Parameter values have been adapted to multiscale clinical trial data in a consistent way, thus supporting the predictive potential of the oncosimulator. Indicative results demonstrating various aspects of the clinical adaptation and validation process are presented. Completion of these processes is expected to pave the way for the clinical translation of the system.

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.

Web-Based Workflow Planning Platform Supporting the Design and Execution of Complex Multiscale Cancer Models

Significant Virtual Physiological Human efforts and projects have been concerned with cancer modeling, especially in the European Commission Seventh Framework research program, with the ambitious goal to approach personalized cancer simulation based on patient-specific data and thereby optimize therapy decisions in the clinical setting. However, building realistic in silicopredictive models targeting the clinical practice requires interactive, synergetic approaches to integrate the currently fragmented efforts emanating from the systems biology and computational oncology communities all around the globe. To further this goal, we propose an intelligent graphical workflow planning system that exploits the multiscale and modular nature of cancer and allows building complex cancer models by intuitively linking/interchanging highly specialized models. The system adopts and extends current standardization efforts, key tools, and infrastructure in view of building a pool of reliable and reproducible models capable of improving current therapies and demonstrating the potential for clinical translation of these technologies.

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.

Studying the growth kinetics of untreated clinical tumors by using an advanced discrete simulation model

Abstract Prior to an eventual clinical adaptation and validation of any clinically oriented model, a thorough study of its dynamic behavior is a sine qua non. Such a study can also elucidate aspects of the interplay of the involved biological mechanisms. Toward this goal, the paper focuses on an in-depth investigation of the free growth behavior of a macroscopically homogeneous malignant tumor system, using a discrete model of tumor growth. We demonstrate that when a clinical tumor grows exponentially, the following preconditions must be fulfilled: (a) time- and space-independent tumor dynamics, in terms of the transition rates among the considered cell categories and the duration of the cell cycle phases, and (b) a tumor system in a state of population equilibrium. Moreover, constant tumor dynamics during the simulation are assumed. In order to create a growing tumor, a condition that the model parameters must fulfill has been derived based on an analytical treatment of the model’s assumptions. A detailed parametric analysis of the model has been performed, in order to determine the impact and the interdependences of its parameters with focus on the free growth rate and the composition of cell population. Constraining tumor cell kinetics, toward limiting the number of possible solutions (i.e., sets of parameters) to the problem of adaptation to the real macroscopic features of a tumor, is also discussed. After completing all parametric studies and after adapting and validating the model on clinical data, it is envisaged to end up with a reliable tool for supporting clinicians in selecting the most appropriate pattern, extracted from several candidate therapeutic schemes, by exploiting tumor- and patient-specific imaging, molecular and histological data.

Translating multiscale cancer models into clinical trials: Simulating breast cancer tumor dynamics within the framework of the “Trial of Principle” clinical trial and the ACGT project.

The potential of cancer multilevel modeling has been particularly emphasized over the past years. Integration of multiscale experimental and clinical information pertaining to cancer via advanced computer models seems to considerably accelerate optimization of cancer treatment in the patient individualized context. However, a sine qua non prerequisite for such models to reach clinical practice is to be thoroughly tested through clinical trials for validation and optimization purposes. This is one of the major goals of the European Commission funded ldquoadvancing clinico-genomic trials on cancerrdquo (ACGT) project. This paper presents a discrete state based, four dimensional, multiscale tumor dynamics model that has been specially developed by the in silico oncology group in order to mimick the trial of principle (TOP) clinical trial concerning breast cancer treated with epirubicin. The TOP trial constitutes one of the ACGT clinical trials. A substantial part of the model can address other tumor types as well. The actual pseudoanonymized imaging, histopathological, molecular and clinical data of the patient are exploited. Special emphasis is put on the effect of cancer stem/clonogenic, progenitor, differentiated and dead cells, the cell category transition rates and the cell category relative populations within the tumor from the treatment baseline onwards. The importance of adaptation of the cell category relative populations to the cell category transition rates for free tumor growth is revealed and the concept of a pertinent nomogram is introduced. A method which ensures adaptation of these two sets of entities at the beginning of the simulation execution is proposed and subsequently successfully applied. Convergence and code checking issues are addressed. Indicative parametric/sensitivity studies are presented along with specific numerical findings. The modelpsilas behavior substantiates its potential to serve as the basis of a treatment optimization system following an eventually succesful completion of the clinical validation and optimization process.

Multilevel cancer modeling in the clinical environment: Simulating the behavior of Wilms tumor in the context of the SIOP 2001/GPOH clinical trial and the ACGT project

Mathematical and computational tumor dynamics models can provide considerable insight into the relative importance and interdependence of related biological mechanisms. They may also suggest the existence of optimal treatment windows in the generic setting. Nevertheless, they cannot be translated into clinical practice unless they undergo a strict and thorough clinical validation and adaptation. In this context one of the major actions of the EC funded project ldquoAdvancing Clinico-Genomic Trials on Cancerrdquo (ACGT) is dedicated to the development of a patient specific four dimensional multiscale tumor model mimicking the nephroblastoma tumor response to chemotherapeutic agents according to the SIOP 2001/GPOH clinical trial. Combined administration of vincristine and dactinomycin is considered. The patient#x2019;s pseudoanonymized imaging, histopathological, molecular and clinical data are carefully exploited. The paper briefly outlines the basics of the model developed by the In Silico Oncology Group and particularly stresses the effect of stem/clonogenic, progenitor and differentiated tumor cells on the overall tumor dynamics. The need for matching the cell category transition rates to the cell category relative populations of free tumor growth for an already large solid tumor at the start of simulation has been clarified. A technique has been suggested and succesfully applied in order to ensure satisfaction of this condition. The concept of a nomogram matching the cell category transition rates to the cell category relative populations at the treatment baseline is introduced. Convergence issues are addressed and indicative numerical results are presented. Qualitative agreement of the modelpsilas behavior with the corresponding clinical trial experience supports its potential to constitute the basis for an optimization system within the clinical environment following completion of its clinical validation and optimization. In silico treatment experimentation in the patient individualized context is expected to constitute the primary application of the model.

In Silico oncology: a top-down multiscale simulator of cancer dynamics. Studying the effect of symmetric stem cell division on the cellular constitution of a tumour

The tremendous rate of accumulation of experimentally and clinically extracted knowledge concerning cancer at all levels of biocomplexity dictates the development of integrative in silico models of tumour dynamics in order to better understand and treat the disease. Since the eventual translation of biomodels into clinical practice presupposes successful clinical validation we have developed a number of multiscale cancer simulation models oriented towards patient individualized treatment optimization. A top-down modelling approach based primarily on discrete event/state simulation has been proposed, developed and implemented. The emerging simulators have been serving as the core simulation modules (oncosimulators) of both the EC funded research projects Contra- Cancrum and ACGT. In this paper a brief outline of the basics of the approach along with paradigmal results demonstrating the effect of symmetric division of cancer stem cells on the cellular constitution of a tumour are presented. The potential and extensibility of the models are discussed.