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

Quantifying ERK activity in response to inhibition of the BRAFV600E-MEK-ERK cascade using mathematical modelling

Funding: UK Medical Research Council (Grant Number(s): MR/R017506/1).Background Simultaneous inhibition of multiple components of the BRAF-MEK-ERK cascade (vertical inhibition) has become a standard of care for treating BRAF-mutant melanoma. However, the molecular mechanism of how vertical inhibition synergistically suppresses intracellular ERK activity, and consequently cell proliferation, are yet to be fully elucidated. Methods We develop a mechanistic mathematical model that describes how the mutant BRAF inhibitor, dabrafenib, and the MEK inhibitor, trametinib, affect BRAFV600E-MEK-ERK signalling. The model is based on a system of chemical reactions that describes cascade signalling dynamics. Using mass action kinetics, the chemical reactions are re-expressed as ordinary differential equations that are parameterised by in vitro data and solved numerically to obtain the temporal evolution of cascade component concentrations. Results The model provides a quantitative method to compute how dabrafenib and trametinib can be used in combination to synergistically inhibit ERK activity in BRAFV600E-mutant melanoma cells. The model elucidates molecular mechanisms of vertical inhibition of the BRAFV600E-MEK-ERK cascade and delineates how elevated BRAF concentrations generate drug resistance to dabrafenib and trametinib. The computational simulations further suggest that elevated ATP levels could be a factor in drug resistance to dabrafenib. Conclusions The model can be used to systematically motivate which dabrafenib–trametinib dose combinations, for treating BRAFV600E-mutated melanoma, warrant experimental investigation.Publisher PDFPeer reviewe

Blackboard-to-bedside:a mathematical modeling bottom-up approach towards 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

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

SH was supported by the Medical Research Council [grant code MR/R017506/1] and Swansea University PhD Research Studentship. URLs: https://mrc.ukri.org/ and https://www.swansea.ac.uk/ MK acknowledges the financial support from the Canadian Institutes of Health Research (CIHR). URL: https://cihr-irsc.gc.ca LJD, AY and PL acknowledge financial support from European Research Council (ERC) advanced grant (ERC-ADG-2015, n 694812 – Hypoximmuno) and EUROSTARS (COMPACT 12053). URLs: https://erc.europa.eu/ and https://www.eurostars-eureka.eu The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.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.Publisher PDFPeer reviewe

Combining hypoxia-activated prodrugs and radiotherapy <i>in silico</i>: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

Quantifying ERK activity in response to inhibition of the BRAFV600E-MEK-ERK cascade using mathematical modelling

BackgroundSimultaneous inhibition of multiple components of the BRAF-MEK-ERK cascade (vertical inhibition) has become a standard of care for treating BRAF-mutant melanoma. However, the molecular mechanism of how vertical inhibition synergistically suppresses intracellular ERK activity, and consequently cell proliferation, are yet to be fully elucidated.MethodsWe develop a mechanistic mathematical model that describes how the mutant BRAF inhibitor, dabrafenib, and the MEK inhibitor, trametinib, affect BRAFV600E-MEK-ERK signalling. The model is based on a system of chemical reactions that describes cascade signalling dynamics. Using mass action kinetics, the chemical reactions are re-expressed as ordinary differential equations that are parameterised by in vitro data and solved numerically to obtain the temporal evolution of cascade component concentrations.ResultsThe model provides a quantitative method to compute how dabrafenib and trametinib can be used in combination to synergistically inhibit ERK activity in BRAFV600E-mutant melanoma cells. The model elucidates molecular mechanisms of vertical inhibition of the BRAFV600E-MEK-ERK cascade and delineates how elevated BRAF concentrations generate drug resistance to dabrafenib and trametinib. The computational simulations further suggest that elevated ATP levels could be a factor in drug resistance to dabrafenib.ConclusionsThe model can be used to systematically motivate which dabrafenib–trametinib dose combinations, for treating BRAFV600E-mutated melanoma, warrant experimental investigation

Spatial cumulant models enable spatially informed treatment strategies and analysis of local interactions in cancer systems

Theoretical and applied cancer studies that use individual-based models (IBMs) have been limited by the lack of a mathematical formulation that enables rigorous analysis of these models. However, spatial cumulant models (SCMs), which have arisen from theoretical ecology, describe population dynamics generated by a specific family of IBMs, namely spatio-temporal point processes (STPPs). SCMs are spatially resolved population models formulated by a system of differential equations that approximate the dynamics of two STPP-generated summary statistics: first-order spatial cumulants (densities), and second-order spatial cumulants (spatial covariances). We exemplify how SCMs can be used in mathematical oncology by modelling theoretical cancer cell populations comprising interacting growth factor-producing and non-producing cells. To formulate model equations, we use computational tools that enable the generation of STPPs, SCMs and mean-field population models (MFPMs) from user-defined model descriptions (Cornell et al. Nat Commun 10:4716, 2019). To calculate and compare STPP, SCM and MFPM-generated summary statistics, we develop an application-agnostic computational pipeline. Our results demonstrate that SCMs can capture STPP-generated population density dynamics, even when MFPMs fail to do so. From both MFPM and SCM equations, we derive treatment-induced death rates required to achieve non-growing cell populations. When testing these treatment strategies in STPP-generated cell populations, our results demonstrate that SCM-informed strategies outperform MFPM-informed strategies in terms of inhibiting population growths. We thus demonstrate that SCMs provide a new framework in which to study cell-cell interactions, and can be used to describe and perturb STPP-generated cell population dynamics. We, therefore, argue that SCMs can be used to increase IBMs' applicability in cancer research.Peer reviewe