Systems oncology: Towards patient-specific treatment regimes informed by multiscale mathematical modelling

Advancing cancer drug development with mechanistic mathematical modeling: bridging the gap between theory and practice

Quantitative predictive modeling of cancer growth, progression, and individual response to therapy is a rapidly growing field. Researchers from mathematical modeling, systems biology, pharmaceutical industry, and regulatory bodies, are collaboratively working on predictive models that could be applied for drug development and, ultimately, the clinical management of cancer patients. A plethora of modeling paradigms and approaches have emerged, making it challenging to compile a comprehensive review across all subdisciplines. It is therefore critical to gauge fundamental design aspects against requirements, and weigh opportunities and limitations of the different model types. In this review, we discuss three fundamental types of cancer models: space-structured models, ecological models, and immune system focused models. For each type, it is our goal to illustrate which mechanisms contribute to variability and heterogeneity in cancer growth and response, so that the appropriate architecture and complexity of a new model becomes clearer. We present the main features addressed by each of the three exemplary modeling types through a subjective collection of literature and illustrative exercises to facilitate inspiration and exchange, with a focus on providing a didactic rather than exhaustive overview. We close by imagining a future multi-scale model design to impact critical decisions in oncology drug development.

Modelling the Tumour Microenvironment, but What Exactly Do We Mean by "Model"?

The Oxford English Dictionary includes 17 definitions for the word "model" as a noun and another 11 as a verb. Therefore, context is necessary to understand the meaning of the word model. For instance, "model railways" refer to replicas of railways and trains at a smaller scale and a "model student" refers to an exemplary individual. In some cases, a specific context, like cancer research, may not be sufficient to provide one specific meaning for model. Even if the context is narrowed, specifically, to research related to the tumour microenvironment, "model" can be understood in a wide variety of ways, from an animal model to a mathematical expression. This paper presents a review of different "models" of the tumour microenvironment, as grouped by different definitions of the word into four categories: model organisms, in vitro models, mathematical models and computational models. Then, the frequencies of different meanings of the word "model" related to the tumour microenvironment are measured from numbers of entries in the MEDLINE database of the United States National Library of Medicine at the National Institutes of Health. The frequencies of the main components of the microenvironment and the organ-related cancers modelled are also assessed quantitatively with specific keywords. Whilst animal models, particularly xenografts and mouse models, are the most commonly used "models", the number of these entries has been slowly decreasing. Mathematical models, as well as prognostic and risk models, follow in frequency, and these have been growing in use.

Modelling Radiation Cancer Treatment with a Death-Rate Term in Ordinary and Fractional Differential Equations

Fractional calculus has recently been applied to the mathematical modelling of tumour growth, but its use introduces complexities that may not be warranted. Mathematical modelling with differential equations is a standard approach to study and predict treatment outcomes for population-level and patient-specific responses. Here, we use patient data of radiation-treated tumours to discuss the benefits and limitations of introducing fractional derivatives into three standard models of tumour growth. The fractional derivative introduces a history-dependence into the growth function, which requires a continuous death-rate term for radiation treatment. This newly proposed radiation-induced death-rate term improves computational efficiency in both ordinary and fractional derivative models. This computational speed-up will benefit common simulation tasks such as model parameterization and the construction and running of virtual clinical trials.

Data-driven spatio-temporal modelling of glioblastoma

Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.

Biological modeling in thermoradiotherapy: present status and ongoing developments toward routine clinical use

Biological modeling for anti-cancer treatments using mathematical models can be very supportive in gaining more insight into dynamic processes responsible for cellular response to treatment, and predicting, evaluating and optimizing therapeutic effects of treatment. This review presents an overview of the current status of biological modeling for hyperthermia in combination with radiotherapy (thermoradiotherapy). Various distinct models have been proposed in the literature, with varying complexity; initially aiming to model the effect of hyperthermia alone, and later on to predict the effect of the combined thermoradiotherapy treatment. Most commonly used models are based on an extension of the linear-quadratic (LQ)-model enabling an easy translation to radiotherapy where the LQ model is widely used. Basic predictions of cell survival have further progressed toward 3 D equivalent dose predictions, i.e., the radiation dose that would be needed without hyperthermia to achieve the same biological effect as the combined thermoradiotherapy treatment. This approach, with the use of temperature-dependent model parameters, allows theoretical evaluation of the effectiveness of different treatment strategies in individual patients, as well as in patient cohorts. This review discusses the significant progress that has been made in biological modeling for hyperthermia combined with radiotherapy. In the future, when adequate temperature-dependent LQ-parameters will be available for a large number of tumor sites and normal tissues, biological modeling can be expected to be of great clinical importance to further optimize combined treatments, optimize clinical protocols and guide further clinical studies.

Cells and material-based strategies for regenerative endodontics

<p class = "Abstract" style = "margin: 0 cm; line-height: 32px; font-size: 12 pt; font-family: "Times New Roman", serif; color: rgb(0, 0, 0); "><span lang = "EN-US">The carious process leads to inflammation of pulp tissue. Current care options include root canal treatment or apexification. These procedures, however, result in the loss of tooth vitality, sensitivity, and healing. Pulp capping and dental pulp regeneration are continually evolving techniques to regenerate pulp tissue, avoiding necrosis and loss of vitality. Many studies have successfully employed stem/progenitor cell populations, revascularization approaches, scaffolds or material-based strategies for pulp regeneration. Here we outline advantages and disadvantages of different methods and techniques which are currently being used in the field of regenerative endodontics. We also summarize recent findings on efficacious peptide-based materials which target the dental niche.<o:p></o:p></span></p>

•

Pulp infection necessitates removal of necrotic, inflamed and infected tissue.

•

Materials used clinically are inert (such as gutta percha, mineral trioxide aggregate).

•

Recent developments in materials (angiogenic hydrogels, stem cell composites) have tuneable bioactivity.

•

Dental pulp regeneration may now be possible through the use of bioactive systems, that guide regeneration.

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

Targeting Cellular DNA Damage Responses in Cancer: An In Vitro-Calibrated Agent-Based Model Simulating Monolayer and Spheroid Treatment Responses to ATR-Inhibiting Drugs

We combine a systems pharmacology approach with an agent-based modelling approach to simulate LoVo cells subjected to AZD6738, an ATR (ataxia–telangiectasia-mutated and rad3-related kinase) inhibiting anti-cancer drug that can hinder tumour proliferation by targeting cellular DNA damage responses. The agent-based model used in this study is governed by a set of empirically observable rules. By adjusting only the rules when moving between monolayer and multi-cellular tumour spheroid simulations, whilst keeping the fundamental mathematical model and parameters intact, the agent-based model is first parameterised by monolayer in vitro data and is thereafter used to simulate treatment responses in in vitro tumour spheroids subjected to dynamic drug delivery. Spheroid simulations are subsequently compared to in vivo data from xenografts in mice. The spheroid simulations are able to capture the dynamics of in vivo tumour growth and regression for approximately 8 days post-tumour injection. Translating quantitative information between in vitro and in vivo research remains a scientifically and financially challenging step in preclinical drug development processes. However, well-developed in silico tools can be used to facilitate this in vitro to in vivo translation, and in this article, we exemplify how data-driven, agent-based models can be used to bridge the gap between in vitro and in vivo research. We further highlight how agent-based models, that are currently underutilised in pharmaceutical contexts, can be used in preclinical drug development.

Bridging cell-scale simulations and radiologic images to explain short-time intratumoral oxygen fluctuations

Radiologic images provide a way to monitor tumor development and its response to therapies in a longitudinal and minimally invasive fashion. However, they operate on a macroscopic scale (average value per voxel) and are not able to capture microscopic scale (cell-level) phenomena. Nevertheless, to examine the causes of frequent fast fluctuations in tissue oxygenation, models simulating individual cells’ behavior are needed. Here, we provide a link between the average data values recorded for radiologic images and the cellular and vascular architecture of the corresponding tissues. Using hybrid agent-based modeling, we generate a set of tissue morphologies capable of reproducing oxygenation levels observed in radiologic images. We then use these in silico tissues to investigate whether oxygen fluctuations can be explained by changes in vascular oxygen supply or by modulations in cellular oxygen absorption. Our studies show that intravascular changes in oxygen supply reproduce the observed fluctuations in tissue oxygenation in all considered regions of interest. However, larger-magnitude fluctuations cannot be recreated by modifications in cellular absorption of oxygen in a biologically feasible manner. Additionally, we develop a procedure to identify plausible tissue morphologies for a given temporal series of average data from radiology images. In future applications, this approach can be used to generate a set of tissues comparable with radiology images and to simulate tumor responses to various anti-cancer treatments at the tissue-scale level.

Low levels of oxygen, called hypoxia, are observable in many solid tumors. They are associated with more aggressive malignant cells that are resistant to chemo-, radio-, and immunotherapies. Recently developed imaging techniques provide a way to measure the magnitude of frequent short-term oxygen fluctuations, but they operate on a macro-scale voxel level. To examine the possible causes of rapid oxygen fluctuations at the cell level, we developed a hybrid agent-based mathematical model. We tested two different mechanisms that may be responsible for these cyclic effects on tissue oxygenation: temporal variations in vascular influx of oxygen and modulations in cellular oxygen absorption. Additionally, we developed a procedure to identify plausible tissue morphologies from data collected from radiological images. This can provide a bridge between the micro-scale simulations with individual cells and the longitudinal medical images containing average values. In future applications, this approach can be used to generate a set of tissues compatible with radiology images and to simulate tumor responses to various anticancer treatments at the cell-scale level.

Patient-Specific Organoid and Organ-on-a-Chip: 3D Cell-Culture Meets 3D Printing and Numerical Simulation

The last few decades have witnessed diversified in vitro models to recapitulate the architecture and function of living organs or tissues and contribute immensely to advances in life science. Two novel 3D cell culture models: 1) Organoid, promoted mainly by the developments of stem cell biology and 2) Organ-on-a-chip, enhanced primarily due to microfluidic technology, have emerged as two promising approaches to advance the understanding of basic biological principles and clinical treatments. This review describes the comparable distinct differences between these two models and provides more insights into their complementarity and integration to recognize their merits and limitations for applicable fields. The convergence of the two approaches to produce multi-organoid-on-a-chip or human organoid-on-a-chip is emerging as a new approach for building 3D models with higher physiological relevance. Furthermore, rapid advancements in 3D printing and numerical simulations, which facilitate the design, manufacture, and results-translation of 3D cell culture models, can also serve as novel tools to promote the development and propagation of organoid and organ-on-a-chip systems. Current technological challenges and limitations, as well as expert recommendations and future solutions to address the promising combinations by incorporating organoids, organ-on-a-chip, 3D printing, and numerical simulation, are also summarized.

Coupled mathematical modeling of cisplatin electroporation

The use of electrochemotherapy (ECT) is a well-established technique to increase the cellular uptake of cytotoxic agents within certain cancer treatment strategies. The study of the mechanisms that take part in this complex process is of high interest to gain a deeper knowledge of it, enabling the improvement of these strategies. In this work, we present a coupled multi-physics electroporation model based on a related previous one, to describe the effect of a set of electric pulses on cisplatin transport across the plasma membrane. The model applies a system of partial differential equations that includes Poisson's equation for the electric field, Nernst-Planck's equation for species transport, Maxwell's tensor and mechanical equilibrium equation for membrane deformation and Smoluchowski's equation for pore creation dynamics. Our numerical results were compared with previous numerical and experimental published data with good qualitative and quantitative agreement. These results indicate that pore aperture is favored at the cell poles by the electric field and mechanical stress forces, giving support to the dominant hypothesis of hydrophilic pore creation as the main mechanism of drug entry during an ECT treatment.

A neuro evolutionary algorithm for patient calibrated prediction of survival in Glioblastoma patients

BACKGROUND AND OBJECTIVES

Glioblastoma multiforme (GBM) is the most common and malignant type of primary brain tumors. Radiation therapy (RT) plus concomitant and adjuvant Temozolomide (TMZ) constitute standard treatment of GBM. Existing models for GBM growth do not consider the effect of different schedules on tumor growth and patient survival. However, clinical trials show that treatment schedule and drug dosage significantly affect patient survival. The goal is to provide a patient calibrated model for predicting survival according to the treatment schedule.

METHODS

We propose a top-down method based on artificial neural networks (ANN) and genetic algorithm (GA) to predict survival of GBM patients. A feed forward undercomplete Autoencoder network is integrated with the neuro-evolutionary (NE) algorithm in order to extract a compressed representation of input clinical data. The proposed NE algorithm uses GA to obtain optimal architecture of a multi-layer perceptron (MLP). Taguchi L₁₆ orthogonal design of experiments is used to tune parameters of the proposed NE algorithm. Finally, the optimal MLP is used to predict survival of GBM patients.

RESULTS

Data from 8 related clinical trials have been collected and integrated to train the model. From 847 evaluable cases, 719 were used for train and validation and the remaining 128 cases were used to test the model. Mean absolute error of the predictions on the test data is 0.087 months which shows excellent performance of the proposed model in predicting survival of the patients. Also, the results show that the proposed NE algorithm is superior to other existing models in both the mean and variability of the prediction error.

Pathway-Based Drug-Repurposing Schemes in Cancer: The Role of Translational Bioinformatics

Cancer is a set of complex pathologies that has been recognized as a major public health problem worldwide for decades. A myriad of therapeutic strategies is indeed available. However, the wide variability in tumor physiology, response to therapy, added to multi-drug resistance poses enormous challenges in clinical oncology. The last years have witnessed a fast-paced development of novel experimental and translational approaches to therapeutics, that supplemented with computational and theoretical advances are opening promising avenues to cope with cancer defiances. At the core of these advances, there is a strong conceptual shift from gene-centric emphasis on driver mutations in specific oncogenes and tumor suppressors-let us call that the silver bullet approach to cancer therapeutics-to a systemic, semi-mechanistic approach based on pathway perturbations and global molecular and physiological regulatory patterns-we will call this the shrapnel approach. The silver bullet approach is still the best one to follow when clonal mutations in driver genes are present in the patient, and when there are targeted therapies to tackle those. Unfortunately, due to the heterogeneous nature of tumors this is not the common case. The wide molecular variability in the mutational level often is reduced to a much smaller set of pathway-based dysfunctions as evidenced by the well-known hallmarks of cancer. In such cases "shrapnel gunshots" may become more effective than "silver bullets". Here, we will briefly present both approaches and will abound on the discussion on the state of the art of pathway-based therapeutic designs from a translational bioinformatics and computational oncology perspective. Further development of these approaches depends on building collaborative, multidisciplinary teams to resort to the expertise of clinical oncologists, oncological surgeons, and molecular oncologists, but also of cancer cell biologists and pharmacologists, as well as bioinformaticians, computational biologists and data scientists. These teams will be capable of engaging on a cycle of analyzing high-throughput experiments, mining databases, researching on clinical data, validating the findings, and improving clinical outcomes for the benefits of the oncological patients.

Combining hypoxia-activated prodrugs and radiotherapy in silico: Impact of treatment scheduling and the intra-tumoural oxygen landscape

Hypoxia-activated prodrugs (HAPs) present a conceptually elegant approach to not only overcome, but better yet, exploit intra-tumoural hypoxia. Despite being successful in vitro and in vivo, HAPs are yet to achieve successful results in clinical settings. It has been hypothesised that this lack of clinical success can, in part, be explained by the insufficiently stringent clinical screening selection of determining which tumours are suitable for HAP treatments. Taking a mathematical modelling approach, we investigate how tumour properties and HAP-radiation scheduling influence treatment outcomes in simulated tumours. The following key results are demonstrated in silico: (i) HAP and ionising radiation (IR) monotherapies may attack tumours in dissimilar, and complementary, ways. (ii) HAP-IR scheduling may impact treatment efficacy. (iii) HAPs may function as IR treatment intensifiers. (iv) The spatio-temporal intra-tumoural oxygen landscape may impact HAP efficacy. Our in silico framework is based on an on-lattice, hybrid, multiscale cellular automaton spanning three spatial dimensions. The mathematical model for tumour spheroid growth is parameterised by multicellular tumour spheroid (MCTS) data.

Blackboard to Bedside: A Mathematical Modeling Bottom-Up Approach Toward Personalized Cancer Treatments

Cancers present with high variability across patients and tumors; thus, cancer care, in terms of disease prevention, detection, and control, can highly benefit from a personalized approach. For a comprehensive personalized oncology practice, this personalization should ideally consider data gathered from various information levels, which range from the macroscale population level down to the microscale tumor level, without omission of the central patient level. Appropriate data mined from each of these levels can significantly contribute in devising personalized treatment plans tailored to the individual patient and tumor. Mathematical models of solid tumors, combined with patient-specific tumor profiles, present a unique opportunity to personalize cancer treatments after detection using a bottom-up approach. Here, we discuss how information harvested from mathematical models and from corresponding in silico experiments can be implemented in preclinical and clinical applications. To conceptually illustrate the power of these models, one such model is presented, and various pertinent tumor and treatment scenarios are demonstrated in silico. The presented model, specifically a multiscale, hybrid cellular automaton, has been fully validated in vitro using multiple cell-line-specific data. We discuss various insights provided by this model and other models like it and their role in designing predictive tools that are both patient, and tumor specific. After refinement and parametrization with appropriate data, such in silico tools have the potential to be used in a clinical setting to aid in treatment protocols and decision making.

Optimising Hydrogel Release Profiles for Viro-Immunotherapy Using Oncolytic Adenovirus Expressing IL-12 and GM-CSF with Immature Dendritic Cells

Sustained-release delivery systems, such as hydrogels, significantly improve cancer therapies by extending the treatment efficacy and avoiding excess wash-out. Combined virotherapy and immunotherapy (viro-immunotherapy) is naturally improved by these sustained-release systems, as it relies on the continual stimulation of the antitumour immune response. In this article, we consider a previously developed viro-immunotherapy treatment where oncolytic viruses that are genetically engineered to infect and lyse cancer cells are loaded onto hydrogels with immature dendritic cells (DCs). The time-dependent release of virus and immune cells results in a prolonged cancer cell killing from both the virus and activated immune cells. Although effective, a major challenge is optimising the release profile of the virus and immature DCs from the gel so as to obtain a minimum tumour size. Using a system of ordinary differential equations calibrated to experimental results, we undertake a novel numerical investigation of different gel-release profiles to determine the optimal release profile for this viro-immunotherapy. Using a data-calibrated mathematical model, we show that if the virus is released rapidly within the first few days and the DCs are released for two weeks, the tumour burden can be significantly decreased. We then find the true optimal gel-release kinetics using a genetic algorithm and suggest that complex profiles present unnecessary risk and that a simple linear-release model is optimal. In this work, insight is provided into a fundamental problem in the growing field of sustained-delivery systems using mathematical modelling and analysis.

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.

Quantitative Predictive Modelling Approaches to Understanding Rheumatoid Arthritis: A Brief Review

Rheumatoid arthritis is a chronic autoimmune disease that is a major public health challenge. The disease is characterised by inflammation of synovial joints and cartilage erosion, which lead to chronic pain, poor life quality and, in some cases, mortality. Understanding the biological mechanisms behind the progression of the disease, as well as developing new methods for quantitative predictions of disease progression in the presence/absence of various therapies is important for the success of therapeutic approaches. The aim of this study is to review various quantitative predictive modelling approaches for understanding rheumatoid arthritis. To this end, we start by briefly discussing the biology of this disease and some current treatment approaches, as well as emphasising some of the open problems in the field. Then, we review various mathematical mechanistic models derived to address some of these open problems. We discuss models that investigate the biological mechanisms behind the progression of the disease, as well as pharmacokinetic and pharmacodynamic models for various drug therapies. Furthermore, we highlight models aimed at optimising the costs of the treatments while taking into consideration the evolution of the disease and potential complications.

Optimization of Cancer Treatment in the Frequency Domain

Thorough exploration of alternative dosing frequencies is often not performed in conventional pharmacometrics approaches. Quantitative systems pharmacology (QSP) can provide novel insights into optimal dosing regimen and drug behaviors which could add a new dimension to the design of novel treatments. However, methods for such an approach are currently lacking. Recently, we illustrated the utility of frequency-domain response analysis (FdRA), an analytical method used in control engineering, using several generic pharmacokinetic-pharmacodynamic case studies. While FdRA is not applicable to models harboring ever increasing variables such as those describing tumor growth, studying such models in the frequency domain provides valuable insight into optimal dosing frequencies. Through the analysis of three distinct tumor growth models (cell cycle-specific, metronomic, and acquired resistance), we demonstrate the application of a simulation-based analysis in the frequency domain to optimize cancer treatments. We study the response of tumor growth to dosing frequencies while simultaneously examining treatment safety, and found for all three models that above a certain dosing frequency, tumor size is insensitive to an increase in dosing frequency, e.g., for the cell cycle-specific model, one dose per 3 days, and an hourly dose yield the same reduction of tumor size to 3% of the initial size after 1 year of treatment. Additionally, we explore the effect of drug elimination rate changes on the tumor growth response. In summary, we show that the frequency-domain view of three models of tumor growth dynamics can help in optimizing drug dosing regimen to improve treatment success.

Lattice and continuum modelling of a bioactive porous tissue scaffold

A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.

Hybrid modeling frameworks of tumor development and treatment

Tumors are complex multicellular heterogeneous systems comprised of components that interact with and modify one another. Tumor development depends on multiple factors: intrinsic, such as genetic mutations, altered signaling pathways, or variable receptor expression; and extrinsic, such as differences in nutrient supply, crosstalk with stromal or immune cells, or variable composition of the surrounding extracellular matrix. Tumors are also characterized by high cellular heterogeneity and dynamically changing tumor microenvironments. The complexity increases when this multiscale, multicomponent system is perturbed by anticancer treatments. Modeling such complex systems and predicting how tumors will respond to therapies require mathematical models that can handle various types of information and combine diverse theoretical methods on multiple temporal and spatial scales, that is, hybrid models. In this update, we discuss the progress that has been achieved during the last 10 years in the area of the hybrid modeling of tumors. The classical definition of hybrid models refers to the coupling of discrete descriptions of cells with continuous descriptions of microenvironmental factors. To reflect on the direction that the modeling field has taken, we propose extending the definition of hybrid models to include of coupling two or more different mathematical frameworks. Thus, in addition to discussing recent advances in discrete/continuous modeling, we also discuss how these two mathematical descriptions can be coupled with theoretical frameworks of optimal control, optimization, fluid dynamics, game theory, and machine learning. All these methods will be illustrated with applications to tumor development and various anticancer treatments.

This article is characterized under:

Analytical and Computational Methods > Computational Methods

Translational, Genomic, and Systems Medicine > Therapeutic Methods

Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models

Towards personalized computational oncology: from spatial models of tumour spheroids, to organoids, to tissues

A main goal of mathematical and computational oncology is to develop quantitative tools to determine the most effective therapies for each individual patient. This involves predicting the right drug to be administered at the right time and at the right dose. Such an approach is known as precision medicine. Mathematical modelling can play an invaluable role in the development of such therapeutic strategies, since it allows for relatively fast, efficient and inexpensive simulations of a large number of treatment schedules in order to find the most effective. This review is a survey of mathematical models that explicitly take into account the spatial architecture of three-dimensional tumours and address tumour development, progression and response to treatments. In particular, we discuss models of epithelial acini, multicellular spheroids, normal and tumour spheroids and organoids, and multi-component tissues. Our intent is to showcase how these in silico models can be applied to patient-specific data to assess which therapeutic strategies will be the most efficient. We also present the concept of virtual clinical trials that integrate standard-of-care patient data, medical imaging, organ-on-chip experiments and computational models to determine personalized medical treatment strategies.

pH as a potential therapeutic target to improve temozolomide antitumor efficacy : A mechanistic modeling study

Despite intensive treatments including temozolomide (TMZ) administration, glioblastoma patient prognosis remains dismal and innovative therapeutic strategies are urgently needed. A systems pharmacology approach was undertaken to investigate TMZ pharmacokinetics-pharmacodynamics (PK-PD) incorporating the effect of local pH, tumor spatial configuration and micro-environment. A hybrid mathematical framework was designed coupling ordinary differential equations describing the intracellular reactions, with a spatial cellular automaton to individualize the cells. A differential drug impact on tumor and healthy cells at constant extracellular pH was computationally demonstrated as TMZ-induced DNA damage was larger in tumor cells as compared to normal cells due to less acidic intracellular pH in cancer cells. Optimality of TMZ efficacy defined as maximum difference between damage in tumor and healthy cells was reached for extracellular pH between 6.8 and 7.5. Next, TMZ PK-PD in a solid tumor was demonstrated to highly depend on its spatial configuration as spread cancer cells or fragmented tumors presented higher TMZ-induced damage as compared to compact tumor spheroid. Simulations highlighted that smaller tumors were less acidic than bigger ones allowing for faster TMZ activation and their closer distance to blood capillaries allowed for better drug penetration. For model parameters corresponding to U87 glioma cells, inter-cell variability in TMZ uptake play no role regarding the mean drug-induced damage in the whole cell population whereas this quantity was increased by inter-cell variability in TMZ efflux which was thus a disadvantage in terms of drug resistance. Overall, this study revealed pH as a new potential target to significantly improve TMZ antitumor efficacy.

Representing dynamic biological networks with multi-scale probabilistic models

Dynamic models analyzing gene regulation and metabolism face challenges when adapted to modeling signal transduction networks. During signal transduction, molecular reactions and mechanisms occur in different spatial and temporal frames and involve feedbacks. This impedes the straight-forward use of methods based on Boolean networks, Bayesian approaches, and differential equations. We propose a new approach, ProbRules, that combines probabilities and logical rules to represent the dynamics of a system across multiple scales. We demonstrate that ProbRules models can represent various network motifs of biological systems. As an example of a comprehensive model of signal transduction, we provide a Wnt network that shows remarkable robustness under a range of phenotypical and pathological conditions. Its simulation allows the clarification of controversially discussed molecular mechanisms of Wnt signaling by predicting wet-lab measurements. ProbRules provides an avenue in current computational modeling by enabling systems biologists to integrate vast amounts of available data on different scales.

The many faces of mathematical modelling in oncology

The application of modelling to solve problems in biology and medicine, and specifically in oncology and radiation therapy, is increasingly established and holds big promise. We provide an overview of the basic concepts of the field and its current state, along with new tools available and future directions for research. We will outline radiobiology models, examples of other anticancer therapy models, multiscale modelling, and we will discuss mechanistic and phenomenological approaches to modelling.

An agent-based model for drug-radiation interactions in the tumour microenvironment: Hypoxia-activated prodrug SN30000 in multicellular tumour spheroids

Multicellular tumour spheroids capture many characteristics of human tumour microenvironments, including hypoxia, and represent an experimentally tractable in vitro model for studying interactions between radiotherapy and anticancer drugs. However, interpreting spheroid data is challenging because of limited ability to observe cell fate within spheroids dynamically. To overcome this limitation, we have developed a hybrid continuum/agent-based model (ABM) for HCT116 tumour spheroids, parameterised using experimental models (monolayers and multilayers) in which reaction and diffusion can be measured directly. In the ABM, cell fate is simulated as a function of local oxygen, glucose and drug concentrations, determined by solving diffusion equations and intracellular reactions. The model is lattice-based, with cells occupying discrete locations on a 3D grid embedded within a coarser grid that encompasses the culture medium; separate solvers are employed for each grid. The generated concentration fields account for depletion in the medium and specify concentration-time profiles within the spheroid. Cell growth and survival are determined by intracellular oxygen and glucose concentrations, the latter based on direct measurement of glucose diffusion/reaction (in multilayers) for the first time. The ABM reproduces known features of spheroids including overall growth rate, its oxygen and glucose dependence, peripheral cell proliferation, central hypoxia and necrosis. We extended the ABM to describe in detail the hypoxia-dependent interaction between ionising radiation and a hypoxia-activated prodrug (SN30000), again using experimentally determined parameters; the model accurately simulated clonogenic cell killing in spheroids, while inclusion of reversible cell cycle delay was required to account for the marked spheroid growth delay after combined radiation and SN30000. This ABM of spheroid growth and response exemplifies the utility of integrating computational and experimental tools for investigating radiation/drug interactions, and highlights the critical importance of understanding oxygen, glucose and drug concentration gradients in interpreting activity of therapeutic agents in spheroid models.

Studies in 3D cultures, notably multicellular tumour spheroids that mimic many features of solid tumours, have great potential for speeding up anticancer drug development. However the increased complexity of 3D cultures makes interpretation of experiments more difficult. We have developed a hybrid continuum/agent-based mathematical model, validated by experiments, to aid interpretation of spheroid experiments in developing drugs designed to eliminate radiation-resistant hypoxic cells. This model includes key features of the tumour microenvironment including oxygen and glucose transport and regions of hypoxia where the cells are resistant to radiation, but sensitive to hypoxia-activated prodrugs such as SN30000. This enables us to predict the growth and cell response in untreated spheroids and compare the results to spheroids treated with radiation and SN30000. We demonstrate good prediction of cellular responses in spheroids treated with radiation and SN30000 and good agreement with spheroid regrowth after treatment when additional effects of cellular growth delay are added. This demonstrates that the modelling approach has potential to improve interpretation of experimental investigations of drug and radiation combinations.

When a duck is not a duck; a new interdisciplinary synthesis for environmental radiation protection

This consensus paper presents the results of a workshop held in Essen, Germany in September 2017, called to examine critically the current approach to radiological environmental protection. The meeting brought together participants from the field of low dose radiobiology and those working in radioecology. Both groups have a common aim of identifying radiation exposures and protecting populations and individuals from harmful effects of ionising radiation exposure, but rarely work closely together. A key question in radiobiology is to understand mechanisms triggered by low doses or dose rates, leading to adverse outcomes of individuals while in radioecology a key objective is to recognise when harm is occurring at the level of the ecosystem. The discussion provided a total of six strategic recommendations which would help to address these questions.

Micropharmacology: An In Silico Approach for Assessing Drug Efficacy Within a Tumor Tissue

Systemic chemotherapy is one of the main anticancer treatments used for most kinds of clinically diagnosed tumors. However, the efficacy of these drugs can be hampered by the physical attributes of the tumor tissue, such as tortuous vasculature, dense and fibrous extracellular matrix, irregular cellular architecture, tumor metabolic gradients, and non-uniform expression of the cell membrane receptors. This can impede the transport of therapeutic agents to tumor cells in sufficient quantities. In addition, tumor microenvironments undergo dynamic spatio-temporal changes during tumor progression and treatment, which can also obstruct drug efficacy. To examine ways to improve drug delivery on a cell-to-tissue scale (single-cell pharmacology), we developed the microscale pharmacokinetics/pharmacodynamics (microPKPD) modeling framework. Our model is modular and can be adjusted to include only the mathematical equations that are crucial for a biological problem under consideration. This modularity makes the model applicable to a broad range of pharmacological cases. As an illustration, we present two specific applications of the microPKPD methodology that help to identify optimal drug properties. The hypoxia-activated drugs example uses continuous drug concentrations, diffusive-advective transport through the tumor interstitium, and passive transmembrane drug uptake. The targeted therapy example represents drug molecules as discrete particles that move by diffusion and actively bind to cell receptors. The proposed modeling approach takes into account the explicit tumor tissue morphology, its metabolic landscape and/or specific receptor distribution. All these tumor attributes can be assessed from patients' diagnostic biopsies; thus, the proposed methodology can be developed into a tool suitable for personalized medicine, such as neoadjuvant chemotherapy.

Combining radiation with hyperthermia: a multiscale model informed by in vitro experiments

Combined radiotherapy and hyperthermia offer great potential for the successful treatment of radio-resistant tumours through thermo-radiosensitization. Tumour response heterogeneity, due to intrinsic, or micro-environmentally induced factors, may greatly influence treatment outcome, but is difficult to account for using traditional treatment planning approaches. Systems oncology simulation, using mathematical models designed to predict tumour growth and treatment response, provides a powerful tool for analysis and optimization of combined treatments. We present a framework that simulates such combination treatments on a cellular level. This multiscale hybrid cellular automaton simulates large cell populations (up to 107 cells) in vitro, while allowing individual cell-cycle progression, and treatment response by modelling radiation-induced mitotic cell death, and immediate cell kill in response to heating. Based on a calibration using a number of experimental growth, cell cycle and survival datasets for HCT116 cells, model predictions agreed well (R2 > 0.95) with experimental data within the range of (thermal and radiation) doses tested (0-40 CEM43, 0-5 Gy). The proposed framework offers flexibility for modelling multimodality treatment combinations in different scenarios. It may therefore provide an important step towards the modelling of personalized therapies using a virtual patient tumour.

The proportion of cancer-related entries in PubMed has increased considerably; is cancer truly "The Emperor of All Maladies"?

In this work, the public database of biomedical literature PubMed was mined using queries with combinations of keywords and year restrictions. It was found that the proportion of Cancer-related entries per year in PubMed has risen from around 6% in 1950 to more than 16% in 2016. This increase is not shared by other conditions such as AIDS, Malaria, Tuberculosis, Diabetes, Cardiovascular, Stroke and Infection some of which have, on the contrary, decreased as a proportion of the total entries per year. Organ-related queries were performed to analyse the variation of some specific cancers. A series of queries related to incidence, funding, and relationship with DNA, Computing and Mathematics, were performed to test correlation between the keywords, with the hope of elucidating the cause behind the rise of Cancer in PubMed. Interestingly, the proportion of Cancer-related entries that contain “DNA”, “Computational” or “Mathematical” have increased, which suggests that the impact of these scientific advances on Cancer has been stronger than in other conditions. It is important to highlight that the results obtained with the data mining approach here presented are limited to the presence or absence of the keywords on a single, yet extensive, database. Therefore, results should be observed with caution. All the data used for this work is publicly available through PubMed and the UK’s Office for National Statistics. All queries and figures were generated with the software platform Matlab and the files are available as supplementary material.

Multiscale image analysis reveals structural heterogeneity of the cell microenvironment in homotypic spheroids

Three-dimensional multicellular aggregates such as spheroids provide reliable in vitro substitutes for tissues. Quantitative characterization of spheroids at the cellular level is fundamental. We present the first pipeline that provides three-dimensional, high-quality images of intact spheroids at cellular resolution and a comprehensive image analysis that completes traditional image segmentation by algorithms from other fields. The pipeline combines light sheet-based fluorescence microscopy of optically cleared spheroids with automated nuclei segmentation (F score: 0.88) and concepts from graph analysis and computational topology. Incorporating cell graphs and alpha shapes provided more than 30 features of individual nuclei, the cellular neighborhood and the spheroid morphology. The application of our pipeline to a set of breast carcinoma spheroids revealed two concentric layers of different cell density for more than 30,000 cells. The thickness of the outer cell layer depends on a spheroid’s size and varies between 50% and 75% of its radius. In differently-sized spheroids, we detected patches of different cell densities ranging from 5 × 10⁵ to 1 × 10⁶cells/mm³. Since cell density affects cell behavior in tissues, structural heterogeneities need to be incorporated into existing models. Our image analysis pipeline provides a multiscale approach to obtain the relevant data for a system-level understanding of tissue architecture.

A SPATIOTEMPORAL MODEL TO SIMULATE CHEMOTHERAPY REGIMENS FOR HETEROGENEOUS BLADDER CANCER METASTASES TO THE LUNG

Tumors are composed of heterogeneous populations of cells. Somatic genetic aberrations are one form of heterogeneity that allows clonal cells to adapt to chemotherapeutic stress, thus providing a path for resistance to arise. In silico modeling of tumors provides a platform for rapid, quantitative experiments to inexpensively study how compositional heterogeneity contributes to drug resistance. Accordingly, we have built a spatiotemporal model of a lung metastasis originating from a primary bladder tumor, incorporating in vivo drug concentrations of first-line chemotherapy, resistance data from bladder cancer cell lines, vascular density of lung metastases, and gains in resistance in cells that survive chemotherapy. In metastatic bladder cancer, a first-line drug regimen includes six cycles of gemcitabine plus cisplatin (GC) delivered simultaneously on day 1, and gemcitabine on day 8 in each 21-day cycle. The interaction between gemcitabine and cisplatin has been shown to be synergistic in vitro, and results in better outcomes in patients. Our model shows that during simulated treatment with this regimen, GC synergy does begin to kill cells that are more resistant to cisplatin, but repopulation by resistant cells occurs. Post-regimen populations are mixtures of the original, seeded resistant clones, and/or new clones that have gained resistance to cisplatin, gemcitabine, or both drugs. The emergence of a tumor with increased resistance is qualitatively consistent with the five-year survival of 6.8% for patients with metastatic transitional cell carcinoma of the urinary bladder treated with a GC regimen. The model can be further used to explore the parameter space for clinically relevant variables, including the timing of drug delivery to optimize cell death, and patient-specific data such as vascular density, rates of resistance gain, disease progression, and molecular profiles, and can be expanded for data on toxicity. The model is specific to bladder cancer, which has not previously been modeled in this context, but can be adapted to represent other cancers.

Tracking the evolution of cancer cell populations through the mathematical lens of phenotype-structured equations

Background

A thorough understanding of the ecological and evolutionary mechanisms that drive the phenotypic evolution of neoplastic cells is a timely and key challenge for the cancer research community. In this respect, mathematical modelling can complement experimental cancer research by offering alternative means of understanding the results of in vitro and in vivo experiments, and by allowing for a quick and easy exploration of a variety of biological scenarios through in silico studies.

Results

To elucidate the roles of phenotypic plasticity and selection pressures in tumour relapse, we present here a phenotype-structured model of evolutionary dynamics in a cancer cell population which is exposed to the action of a cytotoxic drug. The analytical tractability of our model allows us to investigate how the phenotype distribution, the level of phenotypic heterogeneity, and the size of the cell population are shaped by the strength of natural selection, the rate of random epimutations, the intensity of the competition for limited resources between cells, and the drug dose in use.

Conclusions

Our analytical results clarify the conditions for the successful adaptation of cancer cells faced with environmental changes. Furthermore, the results of our analyses demonstrate that the same cell population exposed to different concentrations of the same cytotoxic drug can take different evolutionary trajectories, which culminate in the selection of phenotypic variants characterised by different levels of drug tolerance. This suggests that the response of cancer cells to cytotoxic agents is more complex than a simple binary outcome, i.e., extinction of sensitive cells and selection of highly resistant cells. Also, our mathematical results formalise the idea that the use of cytotoxic agents at high doses can act as a double-edged sword by promoting the outgrowth of drug resistant cellular clones. Overall, our theoretical work offers a formal basis for the development of anti-cancer therapeutic protocols that go beyond the ‘maximum-tolerated-dose paradigm’, as they may be more effective than traditional protocols at keeping the size of cancer cell populations under control while avoiding the expansion of drug tolerant clones.

Reviewers

This article was reviewed by Angela Pisco, Sébastien Benzekry and Heiko Enderling.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0143-4) contains supplementary material, which is available to authorized users.

Mathematical modelling of microtumour infiltration based on in vitro experiments

The present mathematical models of microtumours consider, in general, volumetric growth and spherical tumour invasion shapes. Nevertheless in many cases, such as in gliomas, a need for more accurate delineation of tumour infiltration areas in a patient-specific manner has arisen. The objective of this study was to build a mathematical model able to describe in a case-specific way as well as to predict in a probabilistic way the growth and the real invasion pattern of multicellular tumour spheroids (in vitro model of an avascular microtumour) immersed in a collagen matrix. The two-dimensional theoretical model was represented by a reaction-convection-diffusion equation that considers logistic proliferation, volumetric growth, a rim with proliferative cells at the tumour surface and invasion with diffusive and convective components. Population parameter values of the model were extracted from the experimental dataset and a shape function that describes the invasion area was derived from each experimental case by image processing. New possible and aleatory shape functions were generated by data mining and Monte Carlo tools by means of a satellite EGARCH model, which were fed with all the shape functions of the dataset. Then the main model is used in two different ways: to reproduce the growth and invasion of a given experimental tumour in a case-specific manner when fed with the corresponding shape function (descriptive simulations) or to generate new possible tumour cases that respond to the general population pattern when fed with an aleatory-generated shape function (predictive simulations). Both types of simulations are in good agreement with empirical data, as it was revealed by area quantification and Bland-Altman analysis. This kind of experimental-numerical interaction has wide application potential in designing new strategies able to predict as much as possible the invasive behaviour of a tumour based on its particular characteristics and microenvironment.

Bystander effects and their implications for clinical radiation therapy: Insights from multiscale in silico experiments

Radiotherapy is a commonly used treatment for cancer and is usually given in varying doses. At low radiation doses relatively few cells die as a direct response to radiation but secondary radiation effects, such as DNA mutation or bystander phenomena, may affect many cells. Consequently it is at low radiation levels where an understanding of bystander effects is essential in designing novel therapies with superior clinical outcomes. In this paper, we use a hybrid multiscale mathematical model to study the direct effects of radiation as well as radiation-induced bystander effects on both tumour cells and normal cells. We show that bystander responses play a major role in mediating radiation damage to cells at low-doses of radiotherapy, doing more damage than that due to direct radiation. The survival curves derived from our computational simulations showed an area of hyper-radiosensitivity at low-doses that are not obtained using a traditional radiobiological model.

Microenvironmental Niches and Sanctuaries: A Route to Acquired Resistance

A tumor vasculature that is functionally abnormal results in irregular gradients of metabolites and drugs within the tumor tissue. Recently, significant efforts have been committed to experimentally examine how cellular response to anti-cancer treatments varies based on the environment in which the cells are grown. In vitro studies point to specific conditions in which tumor cells can remain dormant and survive the treatment. In vivo results suggest that cells can escape the effects of drug therapy in tissue regions that are poorly penetrated by the drugs. Better understanding how the tumor microenvironments influence the emergence of drug resistance in both primary and metastatic tumors may improve drug development and the design of more effective therapeutic protocols. This chapter presents a hybrid agent-based model of the growth of tumor micrometastases and explores how microenvironmental factors can contribute to the development of acquired resistance in response to a DNA damaging drug. The specific microenvironments of interest in this work are tumor hypoxic niches and tumor normoxic sanctuaries with poor drug penetration. We aim to quantify how spatial constraints of limited drug transport and quiescent cell survival contribute to the development of drug resistant tumors.

Model-Based Tumor Growth Dynamics and Therapy Response in a Mouse Model of De Novo Carcinogenesis

Tumorigenesis is a complex, multistep process that depends on numerous alterations within the cell and contribution from the surrounding stroma. The ability to model macroscopic tumor evolution with high fidelity may contribute to better predictive tools for designing tumor therapy in the clinic. However, attempts to model tumor growth have mainly been developed and validated using data from xenograft mouse models, which fail to capture important aspects of tumorigenesis including tumor-initiating events and interactions with the immune system. In the present study, we investigate tumor growth and therapy dynamics in a mouse model of de novo carcinogenesis that closely recapitulates tumor initiation, progression and maintenance in vivo. We show that the rate of tumor growth and the effects of therapy are highly variable and mouse specific using a Gompertz model to describe tumor growth and a two-compartment pharmacokinetic/ pharmacodynamic model to describe the effects of therapy in mice treated with 5-FU. We show that inter-mouse growth variability is considerably larger than intra-mouse variability and that there is a correlation between tumor growth and drug kill rates. Our results show that in vivo tumor growth and regression in a double transgenic mouse model are highly variable both within and between subjects and that mathematical models can be used to capture the overall characteristics of this variability. In order for these models to become useful tools in the design of optimal therapy strategies and ultimately in clinical practice, a subject-specific modelling strategy is necessary, rather than approaches that are based on the average behavior of a given subject population which could provide erroneous results.

Diagnostic assessment of osteosarcoma chemoresistance based on Virtual Clinical Trials

Osteosarcoma is the most common primary bone tumor in pediatric and young adult patients. Successful treatment of osteosarcomas requires a combination of surgical resection and systemic chemotherapy, both neoadjuvant (prior to surgery) and adjuvant (after surgery). The degree of necrosis following neoadjuvant chemotherapy correlates with the subsequent probability of disease-free survival. Tumors with less than 10% of viable cells after treatment represent patients with a more favorable prognosis. However, being able to predict early, such as at the time of the pre-treatment tumor biopsy, how the patient will respond to the standard chemotherapy would provide an opportunity for more personalized patient care. Patients with unfavorable predictions could be studied in a protocol, rather than a standard setting, towards improving therapeutic success. The onset of necrotic cells in osteosarcomas treated with chemotherapeutic agents is a measure of tumor sensitivity to the drugs. We hypothesize that the remaining viable cells, i.e., cells that have not responded to the treatment, are chemoresistant, and that the pathological characteristics of these chemoresistant tumor cells within the osteosarcoma pre-treatment biopsy can predict tumor response to the standard-of-care chemotherapeutic treatment. This hypothesis can be tested by comparing patient histopathology samples before, as well as after treatment to identify both morphological and immunochemical cellular features that are characteristic of chemoresistant cells, i.e., cells that survived treatment. Consequently, using computational simulations of dynamic changes in tumor pathology under the simulated standard of care chemotherapeutic treatment, one can couple the pre- and post-treatment morphological and spatial patterns of chemoresistant cells, and correlate them with patient clinical diagnoses. This procedure, that we named 'Virtual Clinical Trials', can serve as a potential predictive biomarker providing a novel value-added decision support tool for oncologists.