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.

Development of a coupled simulation toolkit for computational radiation biology based on Geant4 and CompuCell3D.

Understanding and designing clinical radiation therapy is one of the most important areas of state-of-the-art oncological treatment regimens. Decades of research have gone into developing sophisticated treatment devices and optimization protocols for schedules and dosages. In this paper, we presented a comprehensive computational platform that facilitates building of the sophisticated multi-cell-based model of how radiation affects the biology of living tissue. We designed and implemented a coupled simulation method, including a radiation transport model, and a cell biology model, to simulate the tumor response after irradiation. The radiation transport simulation was implemented through Geant4 which is an open-source Monte Carlo simulation platform that provides many flexibilities for users, as well as low energy DNA damage simulation physics, Geant4-DNA. The cell biology simulation was implemented using CompuCell3D (CC3D) which is a cell biology simulation platform. In order to couple Geant4 solver with CC3D, we developed a 'bridging' module, RADCELL, that extracts tumor cellular geometry of the CC3D simulation (including specification of the individual cells) and ported it to the Geant4 for radiation transport simulation. The cell dose and cell DNA damage distribution in multicellular system were obtained using Geant4. The tumor response was simulated using cell-based tissue models based on CC3D, and the cell dose and cell DNA damage information were fed back through RADCELL to CC3D for updating the cell properties. By merging two powerful and widely used modeling platforms, CC3D and Geant4, we delivered a novel tool that can give us the ability to simulate the dynamics of biological tissue in the presence of ionizing radiation, which provides a framework for quantifying the biological consequences of radiation therapy. In this introductory methods paper, we described our modeling platform in detail and showed how it can be applied to study the application of radiotherapy to a vascularized tumor.

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.

Role of hypoxia-activated prodrugs in combination with radiation therapy: An in silico approach.

Tumour hypoxia has been associated with increased resistance to various cancer treatments, particularly radiation therapy. Conversely, tumour hypoxia is a validated and ideal target for guided cancer drug delivery. For this reason, hypoxia-activated prodrugs (HAPs) have been developed, which remain inactive in the body until in the presence of tissue hypoxia, allowing for an activation tendency in hypoxic regions. We present here an experimentally motivated mathematical model predicting the effectiveness of HAPs in a variety of clinical settings. We first examined HAP effectiveness as a function of the amount of tumour hypoxia and showed that the drugs have a larger impact on tumours with high levels of hypoxia. We then combined HAP treatment with radiation to examine the effects of combination therapies. Our results showed radiation-HAP combination therapies to be more effective against highly hypoxic tumours. The analysis of combination therapies was extended to consider schedule sequencing of the combination treatments. These results suggested that administering HAPs before radiation was most effective in reducing total cell number. Finally, a sensitivity analysis of the drug-related parameters was done to examine the effect of drug diffusivity and enzyme abundance on the overall effectiveness of the drug. Altogether, the results highlight the importance of the knowledge of tumour hypoxia levels before administration of HAPs in order to ensure positive results.

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.

What does not kill a tumour may make it stronger: In silico insights into chemotherapeutic drug resistance.

Tumour recurrence post chemotherapy is an established clinical problem and many cancer types are often observed to be increasingly drug resistant subsequent to chemotherapy treatments. Drug resistance in cancer is a multipart phenomenon which can be derived from several origins and in many cases it has been observed that cancer cells have the ability to possess, acquire and communicate drug resistant traits. Here, an in silico framework is developed in order to study drug resistance and drug response in cancer cell populations exhibiting various drug resistant features. The framework is based on an on-lattice hybrid multiscale mathematical model and is equipped to simulate multiple mechanisms on different scales that contribute towards chemotherapeutic drug resistance in cancer. This study demonstrates how drug resistant tumour features may depend on the interplay amongst intracellular, extracelluar and intercellular factors. On a cellular level, drug resistant cell phenotypes are here derived from inheritance or mutations that are spontaneous, drug-induced or communicated via exosomes. Furthermore intratumoural heterogeneity and spatio-temporal drug dynamics heavily influences drug delivery and the development of drug resistant cancer cell subpopulations. Chemotherapy treatment strategies are here optimised for various in silico tumour scenarios and treatment objectives. We demonstrate that optimal chemotherapy treatment strategies drastically depend on which drug resistant mechanisms are activated, and that furthermore suboptimal chemotherapy administration may promote drug resistance.

Modelling the effects of bacterial cell state and spatial location on tuberculosis treatment: Insights from a hybrid multiscale cellular automaton model.

If improvements are to be made in tuberculosis (TB) treatment, an increased understanding of disease in the lung is needed. Studies have shown that bacteria in a less metabolically active state, associated with the presence of lipid bodies, are less susceptible to antibiotics, and recent results have highlighted the disparity in concentration of different compounds into lesions. Treatment success therefore depends critically on the responses of the individual bacteria that constitute the infection. We propose a hybrid, individual-based approach that analyses spatio-temporal dynamics at the cellular level, linking the behaviour of individual bacteria and host cells with the macroscopic behaviour of the microenvironment. The individual elements (bacteria, macrophages and T cells) are modelled using cellular automaton (CA) rules, and the evolution of oxygen, drugs and chemokine dynamics are incorporated in order to study the effects of the microenvironment in the pathological lesion. We allow bacteria to switch states depending on oxygen concentration, which affects how they respond to treatment. This is the first multiscale model of its type to consider both oxygen-driven phenotypic switching of the Mycobacterium tuberculosis and antibiotic treatment. Using this model, we investigate the role of bacterial cell state and of initial bacterial location on treatment outcome. We demonstrate that when bacteria are located further away from blood vessels, less favourable outcomes are more likely, i.e. longer time before infection is contained/cleared, treatment failure or later relapse. We also show that in cases where bacteria remain at the end of simulations, the organisms tend to be slower-growing and are often located within granulomas, surrounded by caseous material.

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.

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

The multiscale complexity of cancer as a disease necessitates a corresponding multiscale modelling approach to produce truly predictive mathematical models capable of improving existing treatment protocols. To capture all the dynamics of solid tumour growth and its progression, mathematical modellers need to couple biological processes occurring at various spatial and temporal scales (from genes to tissues). Because effectiveness of cancer therapy is considerably affected by intracellular and extracellular heterogeneities as well as by the dynamical changes in the tissue microenvironment, any model attempt to optimise existing protocols must consider these factors ultimately leading to improved multimodal treatment regimes. By improving existing and building new mathematical models of cancer, modellers can play important role in preventing the use of potentially sub-optimal treatment combinations. In this paper, we analyse a multiscale computational mathematical model for cancer growth and spread, incorporating the multiple effects of radiation therapy and chemotherapy in the patient survival probability and implement the model using two different cell based modelling techniques. We show that the insights provided by such multiscale modelling approaches can ultimately help in designing optimal patient-specific multi-modality treatment protocols that may increase patients quality of life.

Towards predicting the response of a solid tumour to chemotherapy and radiotherapy treatments: clinical insights from a computational model.

In this paper we use a hybrid multiscale mathematical model that incorporates both individual cell behaviour through the cell-cycle and the effects of the changing microenvironment through oxygen dynamics to study the multiple effects of radiation therapy. The oxygenation status of the cells is considered as one of the important prognostic markers for determining radiation therapy, as hypoxic cells are less radiosensitive. Another factor that critically affects radiation sensitivity is cell-cycle regulation. The effects of radiation therapy are included in the model using a modified linear quadratic model for the radiation damage, incorporating the effects of hypoxia and cell-cycle in determining the cell-cycle phase-specific radiosensitivity. Furthermore, after irradiation, an individual cell's cell-cycle dynamics are intrinsically modified through the activation of pathways responsible for repair mechanisms, often resulting in a delay/arrest in the cell-cycle. The model is then used to study various combinations of multiple doses of cell-cycle dependent chemotherapies and radiation therapy, as radiation may work better by the partial synchronisation of cells in the most radiosensitive phase of the cell-cycle. Moreover, using this multi-scale model, we investigate the optimum sequencing and scheduling of these multi-modality treatments, and the impact of internal and external heterogeneity on the spatio-temporal patterning of the distribution of tumour cells and their response to different treatment schedules.

Modelling the effects of cell-cycle heterogeneity on the response of a solid tumour to chemotherapy: biological insights from a hybrid multiscale cellular automaton model.

The therapeutic control of a solid tumour depends critically on the responses of the individual cells that constitute the entire tumour mass. A particular cell's spatial location within the tumour and intracellular interactions, including the evolution of the cell-cycle within each cell, has an impact on their decision to grow and divide. They are also influenced by external signals from other cells as well as oxygen and nutrient concentrations. Hence, it is important to take these into account when modelling tumour growth and the response to various treatment regimes ('cell-kill therapies'), including chemotherapy. In order to address this multiscale nature of solid tumour growth and its response to treatment, we propose a hybrid, individual-based approach that analyses spatio-temporal dynamics at the level of cells, linking individual cell behaviour with the macroscopic behaviour of cell organisation and the microenvironment. The individual tumour cells are modelled by using a cellular automaton (CA) approach, where each cell has its own internal cell-cycle, modelled using a system of ODEs. The internal cell-cycle dynamics determine the growth strategy in the CA model, making it more predictive and biologically relevant. It also helps to classify the cells according to their cell-cycle states and to analyse the effect of various cell-cycle dependent cytotoxic drugs. Moreover, we have incorporated the evolution of oxygen dynamics within this hybrid model in order to study the effects of the microenvironment in cell-cycle regulation and tumour treatments. An important factor from the treatment point of view is that the low concentration of oxygen can result in a hypoxia-induced quiescence (G0/G1 arrest) of the cancer cells, making them resistant to key cytotoxic drugs. Using this multiscale model, we investigate the impact of oxygen heterogeneity on the spatio-temporal patterning of the cell distribution and their cell-cycle status. We demonstrate that oxygen transport limitations result in significant heterogeneity in HIF-1 α signalling and cell-cycle status, and when these are combined with drug transport limitations, the efficacy of the therapy is significantly impaired.