An advanced discrete state–discrete event multiscale simulation model of the response of a solid tumor to chemotherapy: Mimicking a clinical study

A comprehensive review of computational cell cycle models in guiding cancer treatment strategies

This article reviews the current knowledge and recent advancements in computational modeling of the cell cycle. It offers a comparative analysis of various modeling paradigms, highlighting their unique strengths, limitations, and applications. Specifically, the article compares deterministic and stochastic models, single-cell versus population models, and mechanistic versus abstract models. This detailed analysis helps determine the most suitable modeling framework for various research needs. Additionally, the discussion extends to the utilization of these computational models to illuminate cell cycle dynamics, with a particular focus on cell cycle viability, crosstalk with signaling pathways, tumor microenvironment, DNA replication, and repair mechanisms, underscoring their critical roles in tumor progression and the optimization of cancer therapies. By applying these models to crucial aspects of cancer therapy planning for better outcomes, including drug efficacy quantification, drug discovery, drug resistance analysis, and dose optimization, the review highlights the significant potential of computational insights in enhancing the precision and effectiveness of cancer treatments. This emphasis on the intricate relationship between computational modeling and therapeutic strategy development underscores the pivotal role of advanced modeling techniques in navigating the complexities of cell cycle dynamics and their implications for cancer therapy.

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.

Analysis of the effectiveness of the treatment of solid tumors in two cases of drug administration

A complete stability analysis of the equilibrium solutions of a system modeling tumor chemotherapy is performed in two cases of administration of the treatment, by continuous infusion and by periodic infusion. Several numerical simulations illustrate and complement the theory.

From tumour perfusion to drug delivery and clinical translation of in silico cancer models

In silico cancer models have demonstrated great potential as a tool to improve drug design, optimise the delivery of drugs to target sites in the host tissue and, hence, improve therapeutic efficacy and patient outcome. However, there are significant barriers to the successful translation of in silico technology from bench to bedside. More precisely, the specification of unknown model parameters, the necessity for models to adequately reflect in vivo conditions, and the limited amount of pertinent validation data to evaluate models' accuracy and assess their reliability, pose major obstacles in the path towards their clinical translation. This review aims to capture the state-of-the-art in in silico cancer modelling of vascularised solid tumour growth, and identify the important advances and barriers to success of these models in clinical oncology. Particular emphasis has been put on continuum-based models of cancer since they - amongst the class of mechanistic spatio-temporal modelling approaches - are well-established in simulating transport phenomena and the biomechanics of tissues, and have demonstrated potential for clinical translation. Three important avenues in in silico modelling are considered in this contribution: first, since systemic therapy is a major cancer treatment approach, we start with an overview of the tumour perfusion and angiogenesis in silico models. Next, we present the state-of-the-art in silico work encompassing the delivery of chemotherapeutic agents to cancer nanomedicines through the bloodstream, and then review continuum-based modelling approaches that demonstrate great promise for successful clinical translation. We conclude with a discussion of what we view to be the key challenges and opportunities for in silico modelling in personalised and precision medicine.

Approaches to combat hypoxia in cancer therapy and the potential for in silico models in their evaluation

AIM

The negative impact of tumour hypoxia on cancer treatment outcome has been long-known, yet there has been little success combating it. This paper investigates the potential role of in silico modelling to help test emerging hypoxia-targeting treatments in cancer therapy.

METHODS

A Medline search was undertaken on the current landscape of in silico models that simulate cancer therapy and evaluate their ability to test hypoxia-targeting treatments. Techniques and treatments to combat tumour hypoxia and their current challenges are also presented.

RESULTS

Hypoxia-targeting treatments include tumour reoxygenation, hypoxic cell radiosensitization with nitroimidazoles, hypoxia-activated prodrugs and molecular targeting. Their main challenges are toxicity and not achieving adequate delivery to hypoxic regions of the tumour. There is promising research toward combining two or more of these techniques. Different types of in silico therapy models have been developed ranging from temporal to spatial and from stochastic to deterministic models. Numerous models have compared the effectiveness of different radiotherapy fractionation schedules for controlling hypoxic tumours. Similarly, models could help identify and optimize new treatments for overcoming hypoxia that utilize novel hypoxia-targeting technology.

CONCLUSION

Current therapy models should attempt to incorporate more sophisticated modelling of tumour angiogenesis/vasculature and vessel perfusion in order to become more useful for testing hypoxia-targeting treatments, which typically rely upon the tumour vasculature for delivery of additional oxygen, (pro)drugs and nanoparticles.

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.

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

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.

Multiscale Modeling in the Clinic: Drug Design and Development

A wide range of length and time scales are relevant to pharmacology, especially in drug development, drug design and drug delivery. Therefore, multiscale computational modeling and simulation methods and paradigms that advance the linkage of phenomena occurring at these multiple scales have become increasingly important. Multiscale approaches present in silico opportunities to advance laboratory research to bedside clinical applications in pharmaceuticals research. This is achievable through the capability of modeling to reveal phenomena occurring across multiple spatial and temporal scales, which are not otherwise readily accessible to experimentation. The resultant models, when validated, are capable of making testable predictions to guide drug design and delivery. In this review we describe the goals, methods, and opportunities of multiscale modeling in drug design and development. We demonstrate the impact of multiple scales of modeling in this field. We indicate the common mathematical and computational techniques employed for multiscale modeling approaches used in pharmacometric and systems pharmacology models in drug development and present several examples illustrating the current state-of-the-art models for (1) excitable systems and applications in cardiac disease; (2) stem cell driven complex biosystems; (3) nanoparticle delivery, with applications to angiogenesis and cancer therapy; (4) host-pathogen interactions and their use in metabolic disorders, inflammation and sepsis; and (5) computer-aided design of nanomedical systems. We conclude with a focus on barriers to successful clinical translation of drug development, drug design and drug delivery multiscale models.

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.

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.

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.

A survey of MRI-based medical image analysis for brain tumor studies

MRI-based medical image analysis for brain tumor studies is gaining attention in recent times due to an increased need for efficient and objective evaluation of large amounts of data. While the pioneering approaches applying automated methods for the analysis of brain tumor images date back almost two decades, the current methods are becoming more mature and coming closer to routine clinical application. This review aims to provide a comprehensive overview by giving a brief introduction to brain tumors and imaging of brain tumors first. Then, we review the state of the art in segmentation, registration and modeling related to tumor-bearing brain images with a focus on gliomas. The objective in the segmentation is outlining the tumor including its sub-compartments and surrounding tissues, while the main challenge in registration and modeling is the handling of morphological changes caused by the tumor. The qualities of different approaches are discussed with a focus on methods that can be applied on standard clinical imaging protocols. Finally, a critical assessment of the current state is performed and future developments and trends are addressed, giving special attention to recent developments in radiological tumor assessment guidelines.

Altered fractionation outcomes for hypoxic head and neck cancer using the HYP-RT Monte Carlo model

OBJECTIVE

Altered fractionation radiotherapy is simulated on a set of virtual tumours to assess the total doses required for tumour control compared with clinical head and neck data and the doses required to control hypoxic vs well-oxygenated tumours with different radiobiological properties.

METHODS

The HYP-RT model is utilised to explore the impact of tumour oxygenation and the onset times of accelerated repopulation (AR) and reoxygenation (ROx) during radiotherapy. A biological effective dose analysis is used to rank the schedules based on their relative normal tissue toxicities.

RESULTS

Altering the onset times of AR and ROx has a large impact on the doses required to achieve tumour control. Immediate onset of ROx and 2-week onset time of AR produce results closely predicting average human outcomes in terms of the total prescription doses in clinical trials. Modifying oxygen enhancement ratio curves based on dose/fraction significantly reduces the dose (5-10 Gy) required for tumour control for hyperfractionated schedules. HYP-RT predicts 10×1.1 Gy per week to be most beneficial, whereas the conventional schedule is predicted as beneficial for early toxicity but has average-poor late toxicity.

CONCLUSION

HYP-RT predicts that altered radiotherapy schedules increase the therapeutic ratio and may be used to make predictions about the prescription doses required to achieve tumour control for tumours with different oxygenation levels and treatment responses.

ADVANCES IN KNOWLEDGE

Oxic and hypoxic tumours have large differences in total radiation dose requirements, affected by AR and ROx onset times by up to 15-25 Gy for the same fractionation schedule.

A multiscale road map of cancer spheroids – incorporating experimental and mathematical modelling to understand cancer progression

Computational models represent a highly suitable framework, not only for testing biological hypotheses and generating new ones but also for optimising experimental strategies. As one surveys the literature devoted to cancer modelling, it is obvious that immense progress has been made in applying simulation techniques to the study of cancer biology, although the full impact has yet to be realised. For example, there are excellent models to describe cancer incidence rates or factors for early disease detection, but these predictions are unable to explain the functional and molecular changes that are associated with tumour progression. In addition, it is crucial that interactions between mechanical effects, and intracellular and intercellular signalling are incorporated in order to understand cancer growth, its interaction with the extracellular microenvironment and invasion of secondary sites. There is a compelling need to tailor new, physiologically relevant in silico models that are specialised for particular types of cancer, such as ovarian cancer owing to its unique route of metastasis, which are capable of investigating anti-cancer therapies, and generating both qualitative and quantitative predictions. This Commentary will focus on how computational simulation approaches can advance our understanding of ovarian cancer progression and treatment, in particular, with the help of multicellular cancer spheroids, and thus, can inform biological hypothesis and experimental design.

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.

Simulating radiotherapy effect in high-grade glioma by using diffusive modeling and brain atlases

Applying diffusive models for simulating the spatiotemporal change of concentration of tumour cells is a modern application of predictive oncology. Diffusive models are used for modelling glioblastoma, the most aggressive type of glioma. This paper presents the results of applying a linear quadratic model for simulating the effects of radiotherapy on an advanced diffusive glioma model. This diffusive model takes into consideration the heterogeneous velocity of glioma in gray and white matter and the anisotropic migration of tumor cells, which is facilitated along white fibers. This work uses normal brain atlases for extracting the proportions of white and gray matter and the diffusion tensors used for anisotropy. The paper also presents the results of applying this glioma model on real clinical datasets.

A novel method for simulating the extracellular matrix in models of tumour growth

A novel hybrid continuum-discrete model to simulate tumour growth on a cellular scale is proposed. The lattice-based spatiotemporal model consists of reaction-diffusion equations that describe interactions between cancer cells and their microenvironment. The fundamental ingredients that are typically considered are the nutrient concentration, the extracellular matrix (ECM), and matrix degrading enzymes (MDEs). The in vivo processes are very complex and occur on different levels. This in turn leads to huge computational costs. The main contribution of the present work is therefore to describe the processes on the basis of simplified mathematical approaches, which, at the same time, depict realistic results to understand the biological processes. In this work, we discuss if we have to simulate the MDE or if the degraded matrix can be estimated directly with respect to the cancer cell distribution. Additionally, we compare the results for modelling tumour growth using the common and our simplified approach, thereby demonstrating the advantages of the proposed method. Therefore, we introduce variations of the positioning of the nutrient delivering blood vessels and use different initializations of the ECM. We conclude that the novel method, which does not explicitly model the matrix degrading enzymes, provides means for a straightforward and fast implementation for modelling tumour growth.

In silico modelling of treatment-induced tumour cell kill: developments and advances

Mathematical and stochastic computer (in silico) models of tumour growth and treatment response of the past and current eras are presented, outlining the aims of the models, model methodology, the key parameters used to describe the tumour system, and treatment modality applied, as well as reported outcomes from simulations. Fractionated radiotherapy, chemotherapy, and combined therapies are reviewed, providing a comprehensive overview of the modelling literature for current modellers and radiobiologists to ignite the interest of other computational scientists and health professionals of the ever evolving and clinically relevant field of tumour modelling.

Integrating intracellular dynamics using CompuCell3D and Bionetsolver: applications to multiscale modelling of cancer cell growth and invasion

In this paper we present a multiscale, individual-based simulation environment that integrates CompuCell3D for lattice-based modelling on the cellular level and Bionetsolver for intracellular modelling. CompuCell3D or CC3D provides an implementation of the lattice-based Cellular Potts Model or CPM (also known as the Glazier-Graner-Hogeweg or GGH model) and a Monte Carlo method based on the metropolis algorithm for system evolution. The integration of CC3D for cellular systems with Bionetsolver for subcellular systems enables us to develop a multiscale mathematical model and to study the evolution of cell behaviour due to the dynamics inside of the cells, capturing aspects of cell behaviour and interaction that is not possible using continuum approaches. We then apply this multiscale modelling technique to a model of cancer growth and invasion, based on a previously published model of Ramis-Conde et al. (2008) where individual cell behaviour is driven by a molecular network describing the dynamics of E-cadherin and β-catenin. In this model, which we refer to as the centre-based model, an alternative individual-based modelling technique was used, namely, a lattice-free approach. In many respects, the GGH or CPM methodology and the approach of the centre-based model have the same overall goal, that is to mimic behaviours and interactions of biological cells. Although the mathematical foundations and computational implementations of the two approaches are very different, the results of the presented simulations are compatible with each other, suggesting that by using individual-based approaches we can formulate a natural way of describing complex multi-cell, multiscale models. The ability to easily reproduce results of one modelling approach using an alternative approach is also essential from a model cross-validation standpoint and also helps to identify any modelling artefacts specific to a given computational approach.

Mathematical modeling of solid cancer growth with angiogenesis

Background

Cancer arises when within a single cell multiple malfunctions of control systems occur, which are, broadly, the system that promote cell growth and the system that protect against erratic growth. Additional systems within the cell must be corrupted so that a cancer cell, to form a mass of any real size, produces substances that promote the growth of new blood vessels. Multiple mutations are required before a normal cell can become a cancer cell by corruption of multiple growth-promoting systems.

Methods

We develop a simple mathematical model to describe the solid cancer growth dynamics inducing angiogenesis in the absence of cancer controlling mechanisms.

Results

The initial conditions supplied to the dynamical system consist of a perturbation in form of pulse: The origin of cancer cells from normal cells of an organ of human body. Thresholds of interacting parameters were obtained from the steady states analysis. The existence of two equilibrium points determine the strong dependency of dynamical trajectories on the initial conditions. The thresholds can be used to control cancer.

Conclusions

Cancer can be settled in an organ if the following combination matches: better fitness of cancer cells, decrease in the efficiency of the repairing systems, increase in the capacity of sprouting from existing vascularization, and higher capacity of mounting up new vascularization. However, we show that cancer is rarely induced in organs (or tissues) displaying an efficient (numerically and functionally) reparative or regenerative mechanism.

Dasatinib synergizes with both cytotoxic and signal transduction inhibitors in heterogeneous breast cancer cell lines--lessons for design of combination targeted therapy

Molecularly targeted therapies have emerged as the leading theme in cancer therapeutics. Multi-cytotoxic drug regimens have been highly successful, yet many studies in targeted therapeutics have centered on a single agent. We investigated whether the Src/Abl kinase inhibitor dasatinib displays synergy with other agents in molecularly heterogeneous breast cancer cell lines. MCF-7, SKBR-3, and MDA-MB-231 display different signaling and gene signatures profiles due to expression of the estrogen receptor, ErbB2, or neither. Cell proliferation was measured following treatment with dasatinib±cytotoxic (paclitaxel, ixabepilone) or molecularly targeted agents (tamoxifen, rapamycin, sorafenib, pan PI3K inhibitor LY294002, and MEK/ERK inhibitor U0126). Dose-responses for single or combination drugs were calculated and analyzed by the Chou-Talalay method. The drugs with the greatest level of synergy with dasatinib were rapamycin, ixabepilone, and sorafenib, for the MDA-MB-231, MCF-7, and SK-BR-3 cell lines respectively. However, dasatinib synergized with both cytotoxic and molecularly targeted agents in all three molecularly heterogeneous breast cancer cell lines. These results suggest that effectiveness of rationally designed therapies may not entirely rest on precise identification of gene signatures or molecular profiling. Since a systems analysis that reveals emergent properties cannot be easily performed for each cancer case, multi-drug regimens in the near future will still involve empirical design.

High-grade glioma diffusive modeling using statistical tissue information and diffusion tensors extracted from atlases

Glioma, especially glioblastoma, is a leading cause of brain cancer fatality involving highly invasive and neoplastic growth. Diffusive models of glioma growth use variations of the diffusion-reaction equation in order to simulate the invasive patterns of glioma cells by approximating the spatiotemporal change of glioma cell concentration. The most advanced diffusive models take into consideration the heterogeneous velocity of glioma in gray and white matter, by using two different discrete diffusion coefficients in these areas. Moreover, by using diffusion tensor imaging (DTI), they simulate the anisotropic migration of glioma cells, which is facilitated along white fibers, assuming diffusion tensors with different diffusion coefficients along each candidate direction of growth. Our study extends this concept by fully exploiting the proportions of white and gray matter extracted by normal brain atlases, rather than discretizing diffusion coefficients. Moreover, the proportions of white and gray matter, as well as the diffusion tensors, are extracted by the respective atlases; thus, no DTI processing is needed. Finally, we applied this novel glioma growth model on real data and the results indicate that prognostication rates can be improved.

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