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

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

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

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

Funding: Cancer Research UK. Research at The Institute of Cancer Research is supported by Cancer Research UK under Programme C33589/A19727. Peter Ziegenhein is supported by Cancer Research UK under Programme C33589/A19908.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.Publisher PDFPeer reviewe

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

This work was supported by the Medical Research Council [grant number MR/P014704/1] and the PreDiCT-TB consortium (IMI Joint undertaking grant agreement number 115337, resources of which are composed of financial contribution from the European Union’s Seventh Framework Programme (FP7/2007-2013) and EFPIA companies’ in kind contribution.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.Publisher PDFPeer reviewe

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

Strategies of Eradicating Glioma Cells: A Multi-Scale Mathematical Model with MiR-451-AMPK-mTOR Control

The cellular dispersion and therapeutic control of glioblastoma, the most aggressive type of primary brain cancer, depends critically on the migration patterns after surgery and intracellular responses of the individual cancer cells in response to external biochemical and biomechanical cues in the microenvironment. Recent studies have shown that a particular microRNA, miR-451, regulates downstream molecules including AMPK and mTOR to determine the balance between rapid proliferation and invasion in response to metabolic stress in the harsh tumor microenvironment. Surgical removal of main tumor is inevitably followed by recurrence of the tumor due to inaccessibility of dispersed tumor cells in normal brain tissue. In order to address this multi-scale nature of glioblastoma proliferation and invasion and its response to conventional treatment, we propose a hybrid model of glioblastoma that analyses spatio-temporal dynamics at the cellular level, linking individual tumor cells with the macroscopic behaviour of cell organization and the microenvironment, and with the intracellular dynamics of miR-451-AMPK-mTOR signaling within a tumour cell. The model identifies a key mechanism underlying the molecular switches between proliferative phase and migratory phase in response to metabolic stress and biophysical interaction between cells in response to fluctuating glucose levels in the presence of blood vessels (BVs). The model predicts that cell migration, therefore efficacy of the treatment, not only depends on oxygen and glucose availability but also on the relative balance between random motility and strength of chemoattractants. Effective control of growing cells near BV sites in addition to relocalization of invisible migratory cells back to the resection site was suggested as a way of eradicating these migratory cells.Publisher PDFPeer reviewe