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

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

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