Computational modelling of perivascular-niche dynamics for the optimization of treatment schedules for glioblastoma

Agent-based modeling in cancer biomedicine: applications and tools for calibration and validation

Computational models are not just appealing because they can simulate and predict the development of biological phenomena across multiple spatial and temporal scales, but also because they can integrate information from well-established in vitro and in vivo models and test new hypotheses in cancer biomedicine. Agent-based models and simulations are especially interesting candidates among computational modeling procedures in cancer research due to the capability to, for instance, recapitulate the dynamics of neoplasia and tumor - host interactions. Yet, the absence of methods to validate the consistency of the results across scales can hinder adoption by turning fine-tuned models into black boxes. This review compiles relevant literature that explores strategies to leverage high-fidelity simulations of multi-scale, or multi-level, cancer models with a focus on verification approached as simulation calibration. We consolidate our review with an outline of modern approaches for agent-based models' validation and provide an ambitious outlook toward rigorous and reliable calibration.

A Multidisciplinary Hyper-Modeling Scheme in Personalized In Silico Oncology: Coupling Cell Kinetics with Metabolism, Signaling Networks, and Biomechanics as Plug-In Component Models of a Cancer Digital Twin

The massive amount of human biological, imaging, and clinical data produced by multiple and diverse sources necessitates integrative modeling approaches able to summarize all this information into answers to specific clinical questions. In this paper, we present a hypermodeling scheme able to combine models of diverse cancer aspects regardless of their underlying method or scale. Describing tissue-scale cancer cell proliferation, biomechanical tumor growth, nutrient transport, genomic-scale aberrant cancer cell metabolism, and cell-signaling pathways that regulate the cellular response to therapy, the hypermodel integrates mutation, miRNA expression, imaging, and clinical data. The constituting hypomodels, as well as their orchestration and links, are described. Two specific cancer types, Wilms tumor (nephroblastoma) and non-small cell lung cancer, are addressed as proof-of-concept study cases. Personalized simulations of the actual anatomy of a patient have been conducted. The hypermodel has also been applied to predict tumor control after radiotherapy and the relationship between tumor proliferative activity and response to neoadjuvant chemotherapy. Our innovative hypermodel holds promise as a digital twin-based clinical decision support system and as the core of future in silico trial platforms, although additional retrospective adaptation and validation are necessary.

Metabolic symbiosis between oxygenated and hypoxic tumour cells: An agent-based modelling study

Deregulated metabolism is one of the hallmarks of cancer. It is well-known that tumour cells tend to metabolize glucose via glycolysis even when oxygen is available and mitochondrial respiration is functional. However, the lower energy efficiency of aerobic glycolysis with respect to mitochondrial respiration makes this behaviour, namely the Warburg effect, counter-intuitive, although it has now been recognized as source of anabolic precursors. On the other hand, there is evidence that oxygenated tumour cells could be fuelled by exogenous lactate produced from glycolysis. We employed a multi-scale approach that integrates multi-agent modelling, diffusion-reaction, stoichiometric equations, and Boolean networks to study metabolic cooperation between hypoxic and oxygenated cells exposed to varying oxygen, nutrient, and inhibitor concentrations. The results show that the cooperation reduces the depletion of environmental glucose, resulting in an overall advantage of using aerobic glycolysis. In addition, the oxygen level was found to be decreased by symbiosis, promoting a further shift towards anaerobic glycolysis. However, the oxygenated and hypoxic populations may gradually reach quasi-equilibrium. A sensitivity analysis using Latin hypercube sampling and partial rank correlation shows that the symbiotic dynamics depends on properties of the specific cell such as the minimum glucose level needed for glycolysis. Our results suggest that strategies that block glucose transporters may be more effective to reduce tumour growth than those blocking lactate intake transporters.

Designing clinical trials for patients who are not average

The heterogeneity inherent in cancer means that even a successful clinical trial merely results in a therapeutic regimen that achieves, on average, a positive result only in a subset of patients. The only way to optimize an intervention for an individual patient is to reframe their treatment as their own, personalized trial. Toward this goal, we formulate a computational framework for performing personalized trials that rely on four mathematical techniques. First, mathematical models that can be calibrated with patient-specific data to make accurate predictions of response. Second, digital twins built on these models capable of simulating the effects of interventions. Third, optimal control theory applied to the digital twins to optimize outcomes. Fourth, data assimilation to continually update and refine predictions in response to therapeutic interventions. In this perspective, we describe each of these techniques, quantify their "state of readiness", and identify use cases for personalized clinical trials.

Predictive digital twin for optimizing patient-specific radiotherapy regimens under uncertainty in high-grade gliomas

We develop a methodology to create data-driven predictive digital twins for optimal risk-aware clinical decision-making. We illustrate the methodology as an enabler for an anticipatory personalized treatment that accounts for uncertainties in the underlying tumor biology in high-grade gliomas, where heterogeneity in the response to standard-of-care (SOC) radiotherapy contributes to sub-optimal patient outcomes. The digital twin is initialized through prior distributions derived from population-level clinical data in the literature for a mechanistic model's parameters. Then the digital twin is personalized using Bayesian model calibration for assimilating patient-specific magnetic resonance imaging data. The calibrated digital twin is used to propose optimal radiotherapy treatment regimens by solving a multi-objective risk-based optimization under uncertainty problem. The solution leads to a suite of patient-specific optimal radiotherapy treatment regimens exhibiting varying levels of trade-off between the two competing clinical objectives: (i) maximizing tumor control (characterized by minimizing the risk of tumor volume growth) and (ii) minimizing the toxicity from radiotherapy. The proposed digital twin framework is illustrated by generating an in silico cohort of 100 patients with high-grade glioma growth and response properties typically observed in the literature. For the same total radiation dose as the SOC, the personalized treatment regimens lead to median increase in tumor time to progression of around six days. Alternatively, for the same level of tumor control as the SOC, the digital twin provides optimal treatment options that lead to a median reduction in radiation dose by 16.7% (10 Gy) compared to SOC total dose of 60 Gy. The range of optimal solutions also provide options with increased doses for patients with aggressive cancer, where SOC does not lead to sufficient tumor control.

Understanding Glioblastoma Signaling, Heterogeneity, Invasiveness, and Drug Delivery Barriers

The most prevalent and aggressive type of brain cancer, namely, glioblastoma (GBM), is characterized by intra- and inter-tumor heterogeneity and strong spreading capacity, which makes treatment ineffective. A true therapeutic answer is still in its infancy despite various studies that have made significant progress toward understanding the mechanisms behind GBM recurrence and its resistance. The primary causes of GBM recurrence are attributed to the heterogeneity and diffusive nature; therefore, monitoring the tumor's heterogeneity and spreading may offer a set of therapeutic targets that could improve the clinical management of GBM and prevent tumor relapse. Additionally, the blood-brain barrier (BBB)-related poor drug delivery that prevents effective drug concentrations within the tumor is discussed. With a primary emphasis on signaling heterogeneity, tumor infiltration, and computational modeling of GBM, this review covers typical therapeutic difficulties and factors contributing to drug resistance development and discusses potential therapeutic approaches.

Phase I study of a novel glioblastoma radiation therapy schedule exploiting cell-state plasticity

BACKGROUND

Glioblastomas comprise heterogeneous cell populations with dynamic, bidirectional plasticity between treatment-resistant stem-like and treatment-sensitive differentiated states, with treatment influencing this process. However, current treatment protocols do not account for this plasticity. Previously, we generated a mathematical model based on preclinical experiments to describe this process and optimize a radiation therapy fractionation schedule that substantially increased survival relative to standard fractionation in a murine glioblastoma model.

METHODS

We developed statistical models to predict the survival benefit of interventions to glioblastoma patients based on the corresponding survival benefit in the mouse model used in our preclinical study. We applied our mathematical model of glioblastoma radiation response to optimize a radiation therapy fractionation schedule for patients undergoing re-irradiation for glioblastoma and developed a first-in-human trial (NCT03557372) to assess the feasibility and safety of administering our schedule.

RESULTS

Our statistical modeling predicted that the hazard ratio when comparing our novel radiation schedule with a standard schedule would be 0.74. Our mathematical modeling suggested that a practical, near-optimal schedule for re-irradiation of recurrent glioblastoma patients was 3.96 Gy × 7 (1 fraction/day) followed by 1.0 Gy × 9 (3 fractions/day). Our optimized schedule was successfully administered to 14/14 (100%) patients.

CONCLUSIONS

A novel radiation therapy schedule based on mathematical modeling of cell-state plasticity is feasible and safe to administer to glioblastoma patients.

HSC Niche Dynamics in Regeneration, Pre-malignancy, and Cancer: Insights From Mathematical Modeling

The hematopoietic stem cell (HSC) niche is a crucial driver of regeneration and malignancy. Its interaction with hematopoietic and malignant stem cells is highly complex and direct experimental observations are challenging. We here develop a mathematical model which helps relate processes in the niche to measurable changes of stem and non-stem cell counts. HSC attached to the niche are assumed to be quiescent. After detachment HSC become activated and divide or differentiate. To maintain their stemness, the progeny originating from division must reattach to the niche. We use mouse data from literature to parametrize the model. By combining mathematical analysis and computer simulations, we systematically investigate the impact of stem cell proliferation, differentiation, niche attachment, and detachment on clinically relevant scenarios. These include bone marrow transplantation, clonal competition, and eradication of malignant cells. According to our model, sampling of blood or bulk marrow provides only limited information about cellular interactions in the niche and the clonal composition of the stem cell population. Furthermore, we investigate how interference with processes in the stem cell niche could help to increase the effect of low-dose chemotherapy or to improve the homing of genetically engineered cells.

An update on the molecular biology of glioblastoma, with clinical implications and progress in its treatment

Glioblastoma multiforme (GBM) is the most aggressive and common malignant primary brain tumor. Patients with GBM often have poor prognoses, with a median survival of ∼15 months. Enhanced understanding of the molecular biology of central nervous system tumors has led to modifications in their classifications, the most recent of which classified these tumors into new categories and made some changes in their nomenclature and grading system. This review aims to give a panoramic view of the last 3 years' findings in glioblastoma characterization, its heterogeneity, and current advances in its treatment. Several molecular parameters have been used to achieve an accurate and personalized characterization of glioblastoma in patients, including epigenetic, genetic, transcriptomic and metabolic features, as well as age- and sex-related patterns and the involvement of several noncoding RNAs in glioblastoma progression. Astrocyte-like neural stem cells and outer radial glial-like cells from the subventricular zone have been proposed as agents involved in GBM of IDH-wildtype origin, but this remains controversial. Glioblastoma metabolism is characterized by upregulation of the PI3K/Akt/mTOR signaling pathway, promotion of the glycolytic flux, maintenance of lipid storage, and other features. This metabolism also contributes to glioblastoma's resistance to conventional therapies. Tumor heterogeneity, a hallmark of GBM, has been shown to affect the genetic expression, modulation of metabolic pathways, and immune system evasion. GBM's aggressive invasion potential is modulated by cell-to-cell crosstalk within the tumor microenvironment and altered expressions of specific genes, such as ANXA2, GBP2, FN1, PHIP, and GLUT3. Nevertheless, the rising number of active clinical trials illustrates the efforts to identify new targets and drugs to treat this malignancy. Immunotherapy is still relevant for research purposes, given the amount of ongoing clinical trials based on this strategy to treat GBM, and neoantigen and nucleic acid-based vaccines are gaining importance due to their antitumoral activity by inducing the immune response. Furthermore, there are clinical trials focused on the PI3K/Akt/mTOR axis, angiogenesis, and tumor heterogeneity for developing molecular-targeted therapies against GBM. Other strategies, such as nanodelivery and computational models, may improve the drug pharmacokinetics and the prognosis of patients with GBM.

Opportunities for improving brain cancer treatment outcomes through imaging-based mathematical modeling of the delivery of radiotherapy and immunotherapy

Immunotherapy has become a fourth pillar in the treatment of brain tumors and, when combined with radiation therapy, may improve patient outcomes and reduce the neurotoxicity. As with other combination therapies, the identification of a treatment schedule that maximizes the synergistic effect of radiation- and immune-therapy is a fundamental challenge. Mechanism-based mathematical modeling is one promising approach to systematically investigate therapeutic combinations to maximize positive outcomes within a rigorous framework. However, successful clinical translation of model-generated combinations of treatment requires patient-specific data to allow the models to be meaningfully initialized and parameterized. Quantitative imaging techniques have emerged as a promising source of high quality, spatially and temporally resolved data for the development and validation of mathematical models. In this review, we will present approaches to personalize mechanism-based modeling frameworks with patient data, and then discuss how these techniques could be leveraged to improve brain cancer outcomes through patient-specific modeling and optimization of treatment strategies.

Agent-based computational modeling of glioblastoma predicts that stromal density is central to oncolytic virus efficacy

Oncolytic viruses (OVs) are emerging cancer immunotherapy. Despite notable successes in the treatment of some tumors, OV therapy for central nervous system cancers has failed to show efficacy. We used an ex vivo tumor model developed from human glioblastoma tissue to evaluate the infiltration of herpes simplex OV rQNestin (oHSV-1) into glioblastoma tumors. We next leveraged our data to develop a computational, model of glioblastoma dynamics that accounts for cellular interactions within the tumor. Using our computational model, we found that low stromal density was highly predictive of oHSV-1 therapeutic success, suggesting that the efficacy of oHSV-1 in glioblastoma may be determined by stromal-to-tumor cell regional density. We validated these findings in heterogenous patient samples from brain metastatic adenocarcinoma. Our integrated modeling strategy can be applied to suggest mechanisms of therapeutic responses for central nervous system cancers and to facilitate the successful translation of OVs into the clinic.

Multiscale Computational Modeling of Vascular Adaptation: A Systems Biology Approach Using Agent-Based Models

The widespread incidence of cardiovascular diseases and associated mortality and morbidity, along with the advent of powerful computational resources, have fostered an extensive research in computational modeling of vascular pathophysiology field and promoted in-silico models as a support for biomedical research. Given the multiscale nature of biological systems, the integration of phenomena at different spatial and temporal scales has emerged to be essential in capturing mechanobiological mechanisms underlying vascular adaptation processes. In this regard, agent-based models have demonstrated to successfully embed the systems biology principles and capture the emergent behavior of cellular systems under different pathophysiological conditions. Furthermore, through their modular structure, agent-based models are suitable to be integrated with continuum-based models within a multiscale framework that can link the molecular pathways to the cell and tissue levels. This can allow improving existing therapies and/or developing new therapeutic strategies. The present review examines the multiscale computational frameworks of vascular adaptation with an emphasis on the integration of agent-based approaches with continuum models to describe vascular pathophysiology in a systems biology perspective. The state-of-the-art highlights the current gaps and limitations in the field, thus shedding light on new areas to be explored that may become the future research focus. The inclusion of molecular intracellular pathways (e.g., genomics or proteomics) within the multiscale agent-based modeling frameworks will certainly provide a great contribution to the promising personalized medicine. Efforts will be also needed to address the challenges encountered for the verification, uncertainty quantification, calibration and validation of these multiscale frameworks.