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

Following publication of the original article [1], the authors noticed that the following errors were introduced by pdf/html formatting issues.

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.

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.

Advancing Colorectal Cancer Diagnosis with AI-Powered Breathomics: Navigating Challenges and Future Directions.

Early detection of colorectal cancer is crucial for improving outcomes and reducing mortality. While there is strong evidence of effectiveness, currently adopted screening methods present several shortcomings which negatively impact the detection of early stage carcinogenesis, including low uptake due to patient discomfort. As a result, developing novel, non-invasive alternatives is an important research priority. Recent advancements in the field of breathomics, the study of breath composition and analysis, have paved the way for new avenues for non-invasive cancer detection and effective monitoring. Harnessing the utility of Volatile Organic Compounds in exhaled breath, breathomics has the potential to disrupt colorectal cancer screening practices. Our goal is to outline key research efforts in this area focusing on machine learning methods used for the analysis of breathomics data, highlight challenges involved in artificial intelligence application in this context, and suggest possible future directions which are currently considered within the framework of the European project ONCOSCREEN.

Immune Phenotype Correlates With Survival in Patients With GBM Treated With Standard Temozolomide-based Therapy and Immunotherapy.

BACKGROUND/AIM: The need for more effective treatment modalities that can improve the clinical outcome of patients with glioblastoma multiforme remains imperative. Dendritic cell vaccination is a fast-developing treatment modality, currently under exploration. Functional immune cell subpopulations may play a role in the final outcome. MATERIALS AND METHODS: Data from 101 patients drawn from the HGG-2010 trial, including baseline patient characteristics and fluorescence-activated cell sorting of immune cell subpopulations, were analyzed by statistical and machine-learning methods. RESULTS: The analysis revealed strong correlations between immune profiles and overall survival, when the extent of resection and the vaccination schedule were used as stratification variables. CONCLUSION: A systematic, in silico workflow detecting strong and statistically significant correlations between overall survival and immune profile-derived quantities obtained at the start of dendritic cell vaccination was devised. The derived correlations could serve as a basis for the identification of prognostic markers discriminating between potential long- and short-term survivors of patients with glioblastoma multiforme.

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

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.

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.

A spatio-temporal simulation model of the response of solid tumours to radiotherapy in vivo: parametric validation concerning oxygen enhancement ratio and cell cycle duration.

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

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 silico predictive 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.

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.

Computational Horizons In Cancer (CHIC): Developing Meta- and Hyper-Multiscale Models and Repositories for In Silico Oncology - a Brief Technical Outline of the Project.

This paper briefly outlines the aim, the objectives, the architecture and the main building blocks of the ongoing large scale integrating transatlantic research project CHIC (http://chic-vph.eu/).

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.

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.

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.

Multiscale modeling for image analysis of brain tumor studies.

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multiscale, multiphysics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlas-based segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

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.

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.