Virtual glioblastoma: growth, migration and treatment in a three-dimensional mathematical model

Biopsy location and tumor-associated macrophages in predicting malignant glioma recurrence using an in-silico model

Predicting the biological behavior and time to recurrence (TTR) of high-grade diffuse gliomas (HGG) after maximum safe neurosurgical resection and combined radiation and chemotherapy plays a pivotal role in planning clinical follow-up, selecting potentially necessary second-line treatment and improving the quality of life for patients diagnosed with a malignant brain tumor. The current standard-of-care (SoC) for HGG includes follow-up neuroradiological imaging to detect recurrence as early as possible and relies on several clinical, neuropathological, and radiological prognostic factors, which have limited accuracy in predicting TTR. In this study, using an in-silico analysis, we aim to improve predictive power for TTR by considering the role of (i) prognostically relevant information available through diagnostics used in the current SoC, (ii) advanced image-based information not currently part of the standard diagnostic workup, such as tumor-normal tissue interface (edge) features and quantitative data specific to biopsy positions within the tumor, and (iii) information on tumor-associated macrophages. In particular, we introduced a state-of-the-art spatio-temporal model of tumor-immune interactions, emphasizing the interplay between macrophages and glioma cells. This model serves as a synthetic reality for assessing the predictive value of various features. We generated a cohort of virtual patients based on our mathematical model. Each patient's dataset includes simulated T1Gd and Fluid-attenuated inversion recovery (FLAIR) MRI volumes. T1-weighted imaging highlights anatomical structures with high contrast, providing clear detail on brain morphology, whereas FLAIR suppresses fluid signals, improving the visualization of lesions near fluid-filled spaces, which is particularly helpful for identifying peritumoral edema. Additionally, we simulated results on macrophage density and proliferative activity, either in a specified part of the tumor, namely the tumor core or edge ("localized"), or unspecified ("non-localized"). To enhance the realism of our synthetic data, we imposed different levels of noise. Our findings reveal that macrophage density at the tumor edge contributed to a high predictive value of feature importance for the selected regression model. Moreover, there are lower MSE values for the "localized" biopsy in prediction accuracy toward recurrence post-resection compared with "non-localized" specimens in the noisy data. In conclusion, the results show that localized biopsies provided more information about tumor behavior, especially at the interface of tumor and normal tissue (Edge).

The self-organized structure of glioma oncostreams and the disruptive role of passive cells

Oncostreams are self-organized structures formed by spindle-like, elongated, self-propelled cells recently described in glioblastomas and especially in gliosarcomas. Cells within these structures either move as large clusters in one main direction, flocks, or as linear, intermingling collections of cells advancing in opposite directions, streams. Round, passive cells are also observed, either inside or segregated from the oncostreams. Here we generalize a recently formulated particle-field approach to investigate the genesis and evolution of these structures, first showing that, in systems consisting only of identical self-propelled cells, both flocks and streams emerge as self-organized dynamic configurations. Flocks are the more stable configurations, while streams are transient and usually originate in collisions between flocks. Stream degradation is easier at low self-propulsion speeds. In systems consisting of both motile and passive cells, the latter block stream formation and accelerate their degradation and flock stabilization. Since the flock appears to be the most effective invasive structure, we thus argue that a phenotype mixture (motile and passive cells) may favor glioblastoma invasion. hlBy relating cellular properties to the observed outcome, our model shows that oncostreams are self-organized structures that result from the interplay between speed, shape, and steric repulsion.

Spatial computational modelling illuminates the role of the tumour microenvironment for treating glioblastoma with immunotherapies

Glioblastoma is the most common and deadliest brain tumour in adults, with a median survival of 15 months under the current standard of care. Immunotherapies like immune checkpoint inhibitors and oncolytic viruses have been extensively studied to improve this endpoint. However, most thus far have failed. To improve the efficacy of immunotherapies to treat glioblastoma, new single-cell imaging modalities like imaging mass cytometry can be leveraged and integrated with computational models. This enables a better understanding of the tumour microenvironment and its role in treatment success or failure in this hard-to-treat tumour. Here, we implemented an agent-based model that allows for spatial predictions of combination chemotherapy, oncolytic virus, and immune checkpoint inhibitors against glioblastoma. We initialised our model with patient imaging mass cytometry data to predict patient-specific responses and found that oncolytic viruses drive combination treatment responses determined by intratumoral cell density. We found that tumours with higher tumour cell density responded better to treatment. When fixing the number of cancer cells, treatment efficacy was shown to be a function of CD4 + T cell and, to a lesser extent, of macrophage counts. Critically, our simulations show that care must be put into the integration of spatial data and agent-based models to effectively capture intratumoral dynamics. Together, this study emphasizes the use of predictive spatial modelling to better understand cancer immunotherapy treatment dynamics, while highlighting key factors to consider during model design and implementation.

Understanding glioblastoma invasion using physically-guided neural networks with internal variables

Microfluidic capacities for both recreating and monitoring cell cultures have opened the door to the use of Data Science and Machine Learning tools for understanding and simulating tumor evolution under controlled conditions. In this work, we show how these techniques could be applied to study Glioblastoma, the deadliest and most frequent primary brain tumor. In particular, we study Glioblastoma invasion using the recent concept of Physically-Guided Neural Networks with Internal Variables (PGNNIV), able to combine data obtained from microfluidic devices and some physical knowledge governing the tumor evolution. The physics is introduced in the network structure by means of a nonlinear advection-diffusion-reaction partial differential equation that models the Glioblastoma evolution. On the other hand, multilayer perceptrons combined with a nodal deconvolution technique are used for learning the go or grow metabolic behavior which characterises the Glioblastoma invasion. The PGNNIV is here trained using synthetic data obtained from in silico tests created under different oxygenation conditions, using a previously validated model. The unravelling capacity of PGNNIV enables discovering complex metabolic processes in a non-parametric way, thus giving explanatory capacity to the networks, and, as a consequence, surpassing the predictive power of any parametric approach and for any kind of stimulus. Besides, the possibility of working, for a particular tumor, with different boundary and initial conditions, permits the use of PGNNIV for defining virtual therapies and for drug design, thus making the first steps towards in silico personalised medicine.

In this work, we apply Physically-Guided Neural Networks with Internal Variables (PGNNIV) to the understanding of the Glioblastoma evolution process. We explain the metabolic changes between the proliferative and migrative activity of Glioblastoma cell cultures by using the go or grow activation functions as a pair of internal variables, whose dependence on the oxygen level is unravelled by some building blocks of the whole PGNNIV. Due to its model-free nature, our method is able to identify different classical mechanistic approaches and to outperform cell culture evolution predictions, as we demonstrate in the paper. Unlike Biologically-Informed Neural Networks we can assimilate data obtained from different boundary conditions and under different external stimuli to simulate the tumor progression under arbitrary conditions. We demonstrate this ability by comparing the predictions with different boundary conditions, resulting in different oxygenation conditions. This flexibility enables the use of our proposed method for personalised medical purposes, as the cell culture metabolic information, for a particular tumor, is encapsulated in a sub-network and may be used for arbitrary in silico tests.

Multiscale modeling of glioma pseudopalisades: contributions from the tumor microenvironment

Gliomas are primary brain tumors with a high invasive potential and infiltrative spread. Among them, glioblastoma multiforme (GBM) exhibits microvascular hyperplasia and pronounced necrosis triggered by hypoxia. Histological samples showing garland-like hypercellular structures (so-called pseudopalisades) centered around the occlusion site of a capillary are typical for GBM and hint on poor prognosis of patient survival. We propose a multiscale modeling approach in the kinetic theory of active particles framework and deduce by an upscaling process a reaction-diffusion model with repellent pH-taxis. We prove existence of a unique global bounded classical solution for a version of the obtained macroscopic system and investigate the asymptotic behavior of the solution. Moreover, we study two different types of scaling and compare the behavior of the obtained macroscopic PDEs by way of simulations. These show that patterns (not necessarily of Turing type), including pseudopalisades, can be formed for some parameter ranges, in accordance with the tumor grade. This is true when the PDEs are obtained via parabolic scaling (undirected tissue), while no such patterns are observed for the PDEs arising by a hyperbolic limit (directed tissue). This suggests that brain tissue might be undirected - at least as far as glioma migration is concerned. We also investigate two different ways of including cell level descriptions of response to hypoxia and the way they are related .

A neuro evolutionary algorithm for patient calibrated prediction of survival in Glioblastoma patients

BACKGROUND AND OBJECTIVES

Glioblastoma multiforme (GBM) is the most common and malignant type of primary brain tumors. Radiation therapy (RT) plus concomitant and adjuvant Temozolomide (TMZ) constitute standard treatment of GBM. Existing models for GBM growth do not consider the effect of different schedules on tumor growth and patient survival. However, clinical trials show that treatment schedule and drug dosage significantly affect patient survival. The goal is to provide a patient calibrated model for predicting survival according to the treatment schedule.

METHODS

We propose a top-down method based on artificial neural networks (ANN) and genetic algorithm (GA) to predict survival of GBM patients. A feed forward undercomplete Autoencoder network is integrated with the neuro-evolutionary (NE) algorithm in order to extract a compressed representation of input clinical data. The proposed NE algorithm uses GA to obtain optimal architecture of a multi-layer perceptron (MLP). Taguchi L₁₆ orthogonal design of experiments is used to tune parameters of the proposed NE algorithm. Finally, the optimal MLP is used to predict survival of GBM patients.

RESULTS

Data from 8 related clinical trials have been collected and integrated to train the model. From 847 evaluable cases, 719 were used for train and validation and the remaining 128 cases were used to test the model. Mean absolute error of the predictions on the test data is 0.087 months which shows excellent performance of the proposed model in predicting survival of the patients. Also, the results show that the proposed NE algorithm is superior to other existing models in both the mean and variability of the prediction error.

Improved detection of diffuse glioma infiltration with imaging combinations: a diagnostic accuracy study

Background

Surgical resection and irradiation of diffuse glioma are guided by standard MRI: T2/fluid attenuated inversion recovery (FLAIR)–weighted MRI for non-enhancing and T1-weighted gadolinium-enhanced (T1G) MRI for enhancing gliomas. Amino acid PET has been suggested as the new standard. Imaging combinations may improve standard MRI and amino acid PET. The aim of the study was to determine the accuracy of imaging combinations to detect glioma infiltration.

Methods

We included 20 consecutive adults with newly diagnosed non-enhancing glioma (7 diffuse astrocytomas, isocitrate dehydrogenase [IDH] mutant; 1 oligodendroglioma, IDH mutant and 1p/19q codeleted; 1 glioblastoma IDH wildtype) or enhancing glioma (glioblastoma, 9 IDH wildtype and 2 IDH mutant). Standardized preoperative imaging (T1-, T2-, FLAIR-weighted, and T1G MRI, perfusion and diffusion MRI, MR spectroscopy and O-(2-[18F]-fluoroethyl)-L-tyrosine ([¹⁸F]FET) PET) was co-localized with multiregion stereotactic biopsies preceding resection. Tumor presence in the biopsies was assessed by 2 neuropathologists. Diagnostic accuracy was determined using receiver operating characteristic analysis.

Results

A total of 174 biopsies were obtained (63 from 9 non-enhancing and 111 from 11 enhancing gliomas), of which 129 contained tumor (50 from non-enhancing and 79 from enhancing gliomas). In enhancing gliomas, the combination of apparent diffusion coefficient (ADC) with [¹⁸F]FET PET (area under the curve [AUC], 95% CI: 0.89, 0.79‒0.99) detected tumor better than T1G MRI (0.56, 0.39‒0.72; P < 0.001) and [¹⁸F]FET PET (0.76, 0.66‒0.86; P = 0.001). In non-enhancing gliomas, no imaging combination detected tumor significantly better than standard MRI. FLAIR-weighted MRI had an AUC of 0.81 (0.65–0.98) compared with 0.69 (0.56–0.81; P = 0.019) for [¹⁸F]FET PET.

Conclusion

Combining ADC and [¹⁸F]FET PET detects glioma infiltration better than standard MRI and [¹⁸F]FET PET in enhancing gliomas, potentially enabling better guidance of local therapy.

A Flux-Limited Model for Glioma Patterning with Hypoxia-Induced Angiogenesis

We propose a model for glioma patterns in a microlocal tumor environment under the influence of acidity, angiogenesis, and tissue anisotropy. The bottom-up model deduction eventually leads to a system of reaction–diffusion–taxis equations for glioma and endothelial cell population densities, of which the former infers flux limitation both in the self-diffusion and taxis terms. The model extends a recently introduced (Kumar, Li and Surulescu, 2020) description of glioma pseudopalisade formation with the aim of studying the effect of hypoxia-induced tumor vascularization on the establishment and maintenance of these histological patterns which are typical for high-grade brain cancer. Numerical simulations of the population level dynamics are performed to investigate several model scenarios containing this and further effects.

Radial basis function-generated finite difference scheme for simulating the brain cancer growth model under radiotherapy in various types of computational domains

BACKGROUND AND OBJECTIVES

We extend the original mathematical model, i.e., Swanson's reaction-diffusion equation to the surfaces with no boundary, and we find a new numerical method based on a meshless approach for solving numerically Swanson's reaction-diffusion model in the square and on the sphere.

METHODS

To solve numerically the Swanson's reaction-diffusion model and its extension version, a collocation meshless technique, namely radial basis function-generated finite difference (RBF-FD) scheme is employed for approximating the spatial variables in the square domain and on the sphere, respectively. Also, to approximate the time variable of the studied models, a first-order semi-implicit backward Euler scheme is used. The resulting fully discrete scheme is a linear system of algebraic equations per time step that is solved via the biconjugate gradient stabilized (BiCGSTAB) iterative algorithm with a zero-fill incomplete lower-upper (ILU) preconditioner.

RESULTS

The numerical simulations show the growth of untreated and treated brain tumors with radiotherapy using estimated and clinical data (given from magnetic resonance imaging (MRI) scans of patients). Moreover, the results reported here can be used for improving the treatment strategies of the invasive brain tumor.

CONCLUSIONS

Using the developed numerical scheme in this paper, we can simulate the behavior of the invasive form of brain tumor response to radiotherapy. Also, we can see the effects of radiation response on the brain tumor cell concentration of individual patients. The proposed meshless technique, which is applied for solving numerically the studied model, does not depend on any background mesh or triangulation for approximation in comparison with mesh-dependent methods. Moreover, we apply this technique to the sphere via any set of distributed points easily.

Reinforcement learning for optimal scheduling of Glioblastoma treatment with Temozolomide

BACKGROUND

Glioblastoma multiforme (GBM) is the most frequent primary brain tumor in adults and Temozolomide (TMZ) is an effective chemotherapeutic agent for its treatment. In Silico models of GBM growth provide an appropriate foundation for analysis and comparison of different regimens. We propose a mathematical frame for patient specific design of optimal chemotherapy regimens for GBM patients.

METHODS

The proposed frame includes online interaction of a virtual GBM with an optimizing agent. Spatiotemporal dynamics of GBM growth and its response to TMZ are simulated with a three dimensional hybrid cellular automaton. Q learning is tailored to the virtual GBM for treatment optimization aimed at minimizing tumor size at the end of treatment course. Q learning consists of a learning agent that interacts with the virtual GBM. System state is affected by the agent decisions and the obtained rewards guide Q learning to the optimal schedule.

RESULTS

Computational results confirm that the optimal chemotherapy schedule depends on some patient specific parameters including body weight, tumor size and its position in the brain. Furthermore, the algorithm is used for scheduling 2100 mg of TMZ on a virtual GBM and the obtained schedule is to administer150 mg of TMZ every other day. The obtained schedule is compared to the standard 7/14 regimen and the results show that it is superior to the 7/14 regimen in minimizing tumor size.

CONCLUSION

The proposed frame is an appropriate decision support system for patient specific design of TMZ administration regimens on GBM patients. Also, since the obtained optimal schedule outperforms the standard 7/14 regimen, it is worthy of further clinical testing.

Simulating glioblastoma growth consisting both visible and invisible parts of the tumor using a diffusion-reaction model followed by resection and radiotherapy

Glioblastoma is known to be among one of the deadliest brain tumors in the world today. There have been major improvements in the detection of cancerous cells in the twenty-first century. However, the threshold of detection of these cancerous cells varies in different scanning techniques such as magnetic resonance imaging (MRI) and computed tomography (CT). The growth of these tumors and different treatments have been modeled to assist medical experts in better predictions of the related tumor growth and in the selection of more accurate treatments. In clinical terms the tumor consisted of two parts known as the visible part, which is the part of the tumor that is above the threshold of the detecting device and the invisible part, which is below the detecting threshold. In this study, the common reaction-diffusion model of tumor growth is used to simulate the growth of the glioblastoma tumor. Also resection and radiotherapy have been modeled as methods to prevent the growth of the tumor. The results demonstrate that although the selected treatments were effective in reducing the number of cancerous cells to under the threshold of detection, they did not eliminate all cancerous cells and if no further treatments were applied, the cancerous cells would spread and become malignant again. Although previous studies have suggested that the ratio of proliferation to diffusion could describe the malignancy of the tumor, this study in addition shows the importance of each of the coefficients regarding the malignancy of the tumor.

Radiation protraction schedules for low-grade gliomas: a comparison between different mathematical models

We optimize radiotherapy (RT) administration strategies for treating low-grade gliomas. Specifically, we consider different tumour growth laws, both with and without spatial effects. In each scenario, we find the optimal treatment in the sense of maximizing the overall survival time of a virtual low-grade glioma patient, whose tumour progresses according to the examined growth laws. We discover that an extreme protraction therapeutic strategy, which amounts to substantially extending the time interval between RT sessions, may lead to better tumour control. The clinical implications of our results are also presented.

Patient-specific parameter estimates of glioblastoma multiforme growth dynamics from a model with explicit birth and death rates

Glioblastoma multiforme (GBM) is an aggressive primary brain cancer with a grim prog-nosis. Its morphology is heterogeneous, but prototypically consists of an inner, largely necrotic core surrounded by an outer, contrast-enhancing rim, and often extensive tumor-associated edema beyond. This structure is usually demonstrated by magnetic resonance imaging (MRI). To help relate the three highly idealized components of GBMs (i.e., necrotic core, enhancing rim, and maximum edema ex-tent) to the underlying growth "laws," a mathematical model of GBM growth with explicit motility, birth, and death processes is proposed. This model generates a traveling-wave solution that mimics tumor progression. We develop several novel methods to approximate key characteristics of the wave profile, which can be compared with MRI data. Several simplified forms of growth and death terms and their parameter identifiability are studied. We use several test cases of MRI data of GBM patients to yield personalized parameterizations of the model, and the biological and clinical implications are discussed.

Advances in Glioblastoma Multiforme Treatment: New Models for Nanoparticle Therapy

The most lethal form of brain cancer, glioblastoma multiforme, is characterized by rapid growth and invasion facilitated by cell migration and degradation of the extracellular matrix. Despite technological advances in surgery and radio-chemotherapy, glioblastoma remains largely resistant to treatment. New approaches to study glioblastoma and to design optimized therapies are greatly needed. One such approach harnesses computational modeling to support the design and delivery of glioblastoma treatment. In this paper, we critically summarize current glioblastoma therapy, with a focus on emerging nanomedicine and therapies that capitalize on cell-specific signaling in glioblastoma. We follow this summary by discussing computational modeling approaches focused on optimizing these emerging nanotherapeutics for brain cancer. We conclude by illustrating how mathematical analysis can be used to compare the delivery of a high potential anticancer molecule, delphinidin, in both free and nanoparticle loaded forms across the blood-brain barrier for glioblastoma.

A Brownian dynamics tumor progression simulator with application to glioblastoma

Tumor progression modeling offers the potential to predict tumor-spreading behavior to improve prognostic accuracy and guide therapy development. Common simulation methods include continuous reaction-diffusion (RD) approaches that capture mean spatio-temporal tumor spreading behavior and discrete agent-based (AB) approaches which capture individual cell events such as proliferation or migration. The brain cancer glioblastoma (GBM) is especially appropriate for such proliferation-migration modeling approaches because tumor cells seldom metastasize outside of the central nervous system and cells are both highly proliferative and migratory. In glioblastoma research, current RD estimates of proliferation and migration parameters are derived from computed tomography or magnetic resonance images. However, these estimates of glioblastoma cell migration rates, modeled as a diffusion coefficient, are approximately 1-2 orders of magnitude larger than single-cell measurements in animal models of this disease. To identify possible sources for this discrepancy, we evaluated the fundamental RD simulation assumptions that cells are point-like structures that can overlap. To give cells physical size (~10 μm), we used a Brownian dynamics approach that simulates individual single-cell diffusive migration, growth, and proliferation activity via a gridless, off-lattice, AB method where cells can be prohibited from overlapping each other. We found that for realistic single-cell parameter growth and migration rates, a non-overlapping model gives rise to a jammed configuration in the center of the tumor and a biased outward diffusion of cells in the tumor periphery, creating a quasi-ballistic advancing tumor front. The simulations demonstrate that a fast-progressing tumor can result from minimally diffusive cells, but at a rate that is still dependent on single-cell diffusive migration rates. Thus, modeling with the assumption of physically-grounded volume conservation can account for the apparent discrepancy between estimated and measured diffusion of GBM cells and provide a new theoretical framework that naturally links single-cell growth and migration dynamics to tumor-level progression.

Non-standard radiotherapy fractionations delay the time to malignant transformation of low-grade gliomas

Grade II gliomas are slowly growing primary brain tumors that affect mostly young patients. Cytotoxic therapies (radiotherapy and/or chemotherapy) are used initially only for patients having a bad prognosis. These therapies are planned following the “maximum dose in minimum time” principle, i. e. the same schedule used for high-grade brain tumors in spite of their very different behavior. These tumors transform after a variable time into high-grade gliomas, which significantly decreases the patient’s life expectancy. In this paper we study mathematical models describing the growth of grade II gliomas in response to radiotherapy. We find that protracted metronomic fractionations, i.e. therapeutical schedules enlarging the time interval between low-dose radiotherapy fractions, may lead to a better tumor control without an increase in toxicity. Other non-standard fractionations such as protracted or hypoprotracted schemes may also be beneficial. The potential survival improvement depends on the tumor’s proliferation rate and can be even of the order of years. A conservative metronomic scheme, still being a suboptimal treatment, delays the time to malignant progression by at least one year when compared to the standard scheme.

Mathematical Analysis of Glioma Growth in a Murine Model

Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm³ to 62 mm³, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.

A Feasibility Study of Personalized Prescription Schemes for Glioblastoma Patients Using a Proliferation and Invasion Glioma Model

Purpose: This study investigates the feasibility of personalizing radiotherapy prescription schemes (treatment margins and fractional doses) for glioblastoma (GBM) patients and their potential benefits using a proliferation and invasion (PI) glioma model on phantoms. Methods and Materials: We propose a strategy to personalize radiotherapy prescription schemes by simulating the proliferation and invasion of the tumor in 2D according to the PI glioma model. We demonstrate the strategy and its potential benefits by presenting virtual cases, where the standard and personalized prescriptions were applied to the tumor. Standard prescription was assumed to deliver 46 Gy in 23 fractions to the initial, gross tumor volume (GTV₁) plus a 2 cm margin and an additional 14 Gy in 7 fractions to the boost GTV₂ plus a 2 cm margin. The virtual cases include the tumors with a moving velocity of 0.029 (slow-move), 0.079 (average-move), and 0.13 (fast-move) mm/day for the gross tumor volume (GTV) with a radius of 1 (small) and 2 (large) cm. For each tumor size and velocity, the margin around GTV₁ and GTV₂ was varied between 0–6 cm and 1–3 cm, respectively. Equivalent uniform dose (EUD) to normal brain was constrained to the EUD value obtained by using the standard prescription. Various linear dose policies, where the fractional dose is linearly decreasing, constant, or increasing, were investigated to estimate the temporal effect of the radiation dose on tumor cell-kills. The goal was to find the combination of margins for GTV₁ and GTV₂ and a linear dose policy, which minimize the tumor cell-surviving fraction (SF) under a normal tissue constraint. The efficacy of a personalized prescription was evaluated by tumor EUD and the estimated survival time. Results: The personalized prescription for the slow-move tumors was to use 3.0–3.5 cm margins for GTV₁, and a 1.5 cm margin for GTV₂. For the average- and fast-move tumors, it was optimal to use a 6.0 cm margin for GTV₁ and then 1.5–3.0 cm margins for GTV₂, suggesting a course of whole brain therapy followed by a boost to a smaller volume. It was more effective to deliver the boost sequentially using a linearly decreasing fractional dose for all tumors. Personalized prescriptions led to surviving fractions of 0.001–0.465% compared to the standard prescription, and increased the tumor EUDs by 25.3–49.3% and estimated survival times by 7.6–22.2 months. Conclusions: Personalizing treatment margins based on the measured proliferative capacity of GBM tumor cells can potentially lead to significant improvements in tumor cell kill and related clinical outcomes.

A Numerical Handling of the Boundary Conditions Imposed by the Skull on an Inhomogeneous Diffusion-Reaction Model of Glioblastoma Invasion Into the Brain: Clinical Validation Aspects

A novel explicit triscale reaction-diffusion numerical model of glioblastoma multiforme tumor growth is presented. The model incorporates the handling of Neumann boundary conditions imposed by the cranium and takes into account both the inhomogeneous nature of human brain and the complexity of the skull geometry. The finite-difference time-domain method is adopted. To demonstrate the workflow of a possible clinical validation procedure, a clinical case/scenario is addressed. A good agreement of the in silico calculated value of the doubling time (ie, the time for tumor volume to double) with the value of the same quantity based on tomographic imaging data has been observed. A theoretical exploration suggests that a rough but still quite informative value of the doubling time may be calculated based on a homogeneous brain model. The model could serve as the main component of a continuous mathematics-based glioblastoma oncosimulator aiming at supporting the clinician in the optimal patient-individualized design of treatment using the patient's multiscale data and experimenting in silico (ie, on the computer).

Fast and high temperature hyperthermia coupled with radiotherapy as a possible new treatment for glioblastoma

Background

A new transcranial focused ultrasound device has been developed that can induce hyperthermia in a large tissue volume. The purpose of this work is to investigate theoretically how glioblastoma multiforme (GBM) can be effectively treated by combining the fast hyperthermia generated by this focused ultrasound device with external beam radiotherapy.

Methods/Design

To investigate the effect of tumor growth, we have developed a mathematical description of GBM proliferation and diffusion in the context of reaction–diffusion theory. In addition, we have formulated equations describing the impact of radiotherapy and heat on GBM in the reaction–diffusion equation, including tumor regrowth by stem cells. This formulation has been used to predict the effectiveness of the combination treatment for a realistic focused ultrasound heating scenario.

Our results show that patient survival could be significantly improved by this combined treatment modality.

Discussion

High priority should be given to experiments to validate the therapeutic benefit predicted by our model.

Electronic supplementary material

The online version of this article (doi:10.1186/s40349-016-0078-3) contains supplementary material, which is available to authorized users.

Why one-size-fits-all vaso-modulatory interventions fail to control glioma invasion: in silico insights

Gliomas are highly invasive brain tumours characterised by poor prognosis and limited response to therapy. There is an ongoing debate on the therapeutic potential of vaso-modulatory interventions against glioma invasion. Prominent vasculature-targeting therapies involve tumour blood vessel deterioration and normalisation. The former aims at tumour infarction and nutrient deprivation induced by blood vessel occlusion/collapse. In contrast, the therapeutic intention of normalising the abnormal tumour vasculature is to improve the efficacy of conventional treatment modalities. Although these strategies have shown therapeutic potential, it remains unclear why they both often fail to control glioma growth. To shed some light on this issue, we propose a mathematical model based on the migration/proliferation dichotomy of glioma cells in order to investigate why vaso-modulatory interventions have shown limited success in terms of tumour clearance. We found the existence of a critical cell proliferation/diffusion ratio that separates glioma responses to vaso-modulatory interventions into two distinct regimes. While for tumours, belonging to one regime, vascular modulations reduce the front speed and increase the infiltration width, for those in the other regime, the invasion speed increases and infiltration width decreases. We discuss how these in silico findings can be used to guide individualised vaso-modulatory approaches to improve treatment success rates.

Computational Trials: Unraveling Motility Phenotypes, Progression Patterns, and Treatment Options for Glioblastoma Multiforme

Glioblastoma multiforme is a malignant brain tumor with poor prognosis and high morbidity due to its invasiveness. Hypoxia-driven motility and concentration-driven motility are two mechanisms of glioblastoma multiforme invasion in the brain. The use of anti-angiogenic drugs has uncovered new progression patterns of glioblastoma multiforme associated with significant differences in overall survival. Here, we apply a mathematical model of glioblastoma multiforme growth and invasion in humans and design computational trials using agents that target angiogenesis, tumor replication rates, or motility. The findings link highly-dispersive, moderately-dispersive, and hypoxia-driven tumors to the patterns observed in glioblastoma multiforme treated by anti-angiogenesis, consisting of progression by Expanding FLAIR, Expanding FLAIR + Necrosis, and Expanding Necrosis, respectively. Furthermore, replication rate-reducing strategies (e.g. Tumor Treating Fields) appear to be effective in highly-dispersive and moderately-dispersive tumors but not in hypoxia-driven tumors. The latter may respond to motility-reducing agents. In a population computational trial, with all three phenotypes, a correlation was observed between the efficacy of the rate-reducing agent and the prolongation of overall survival times. This research highlights the potential applications of computational trials and supports new hypotheses on glioblastoma multiforme phenotypes and treatment options.

A mathematical approach to virus therapy of glioblastomas

BACKGROUND

It is widely believed that the treatment of glioblastomas (GBM) could benefit from oncolytic virus therapy. Clinical research has shown that Vesicular Stomatitis Virus (VSV) has strong oncolytic properties. In addition, mathematical models of virus treatment of tumors have been developed in recent years. Some experiments in vitro and in vivo have been done and shown promising results, but have been never compared quantitatively with mathematical models. We use in vitro data of this virus applied to glioblastoma.

RESULTS

We describe three increasingly realistic mathematical models for the VSV-GBM in vitro experiment with progressive incorporation of time-delay effects. For the virus dynamics, we obtain results consistent with the in vitro experimental speed data only when applying the more complex and comprehensive model, with time-delay effects both in the reactive and diffusive terms. The tumor speed is given by the minimum of a very simple function that nonetheless yields results within the experimental measured range.

CONCLUSIONS

We have improved a previous model with new ideas and carefully incorporated concepts from experimental results. We have shown that the delay time τ is the crucial parameter in this kind of models. We have demonstrated that our new model can satisfactorily predict the front speed for the lytic action of oncolytic VSV on glioblastoma observed in vitro. We provide a basis that can be applied in the near future to realistically simulate in vivo virus treatments of several cancers.

Mathematically modeling the biological properties of gliomas: A review

Although mathematical modeling is a mainstay for industrial and many scientific studies, such approaches have found little application in neurosurgery. However, the fusion of biological studies and applied mathematics is rapidly changing this environment, especially for cancer research. This review focuses on the exciting potential for mathematical models to provide new avenues for studying the growth of gliomas to practical use. In vitro studies are often used to simulate the effects of specific model parameters that would be difficult in a larger-scale model. With regard to glioma invasive properties, metabolic and vascular attributes can be modeled to gain insight into the infiltrative mechanisms that are attributable to the tumor's aggressive behavior. Morphologically, gliomas show different characteristics that may allow their growth stage and invasive properties to be predicted, and models continue to offer insight about how these attributes are manifested visually. Recent studies have attempted to predict the efficacy of certain treatment modalities and exactly how they should be administered relative to each other. Imaging is also a crucial component in simulating clinically relevant tumors and their influence on the surrounding anatomical structures in the brain.

Evaluation of current clinical target volume definitions for glioblastoma using cell-based dosimetry stochastic methods

OBJECTIVE

Determination of an optimal clinical target volume (CTV) is complex and remains uncertain. The aim of this study was to develop a glioblastoma multiforme (GBM) model to be used for evaluation of current CTV practices for external radiotherapy.

METHODS

The GBM model was structured as follows: (1) a Geant4 cellular model was developed to calculate the absorbed dose in individual cells represented by cubic voxels of 20 μm sides. The system was irradiated with opposing 6 MV X-ray beams. The beams encompassed planning target volumes corresponding to 2.0- and 2.5-cm CTV margins; (2) microscopic extension probability (MEP) models were developed using MATLAB(®) 2012a (MathWorks(®), Natick, MA), based on clinical studies reporting on GBM clonogenic spread; (3) the cellular dose distribution was convolved with the MEP models to evaluate cellular survival fractions (SFs) for both CTV margins.

RESULTS

A CTV margin of 2.5 cm, compared to a 2.0-cm CTV margin, resulted in a reduced total SF from 12.9% ± 0.9% to 3.6% ± 0.2%, 5.5% ± 0.4% to 1.2% ± 0.1% and 11.1% ± 0.7% to 3.0% ± 0.2% for circular, elliptical and irregular MEP distributions, respectively.

CONCLUSION

A Monte Carlo model was developed to quantitatively evaluate the impact of GBM CTV margins on total and penumbral SF. The results suggest that the reduction in total SF ranges from 3.5 to 5, when the CTV is extended by 0.5 cm.

ADVANCES IN KNOWLEDGE

The model provides a quantitative tool for evaluation of different CTV margins in terms of cell kill efficacy. Cellular platform of the tool allows future incorporation of cellular properties of GBM.

Extreme protraction for low-grade gliomas: theoretical proof of concept of a novel therapeutical strategy

Grade II gliomas are slowly growing primary brain tumours that affect mostly young patients and become fatal after a variable time period. Current clinical handling includes surgery as first-line treatment. Cytotoxic therapies (radiotherapy RT or chemotherapy QT) are used initially only for patients having a bad prognosis. Therapies are administered following the 'maximum dose in minimum time' principle, which is the same schedule used for high-grade brain tumours. Using mathematical models describing the growth of these tumours in response to radiotherapy, we find that an extreme protraction therapeutical strategy, i.e. enlarging substantially the time interval between RT fractions, may lead to better tumour control. Explicit formulas are found providing the optimal spacing between doses in a very good agreement with the simulations of the full 3D mathematical model approximating the tumour spatiotemporal dynamics. This idea, although breaking the well-established paradigm, has biological meaning since, in these slowly growing tumours, it may be more favourable to treat the tumour as the tumour cells leave the quiescent compartment and move into the cell cycle.

A mathematical model to elucidate brain tumor abrogation by immunotherapy with T11 target structure

T11 Target structure (T11TS), a membrane glycoprotein isolated from sheep erythrocytes, reverses the immune suppressed state of brain tumor induced animals by boosting the functional status of the immune cells. This study aims at aiding in the design of more efficacious brain tumor therapies with T11 target structure. We propose a mathematical model for brain tumor (glioma) and the immune system interactions, which aims in designing efficacious brain tumor therapy. The model encompasses considerations of the interactive dynamics of glioma cells, macrophages, cytotoxic T-lymphocytes (CD8+ T-cells), TGF-β, IFN-γ and the T11TS. The system undergoes sensitivity analysis, that determines which state variables are sensitive to the given parameters and the parameters are estimated from the published data. Computer simulations were used for model verification and validation, which highlight the importance of T11 target structure in brain tumor therapy.

Quiescence as an explanation of Gompertzian tumor growth revisited

Gompertz's empirical equation remains the most popular one in describing cancer cell population growth in a wide spectrum of bio-medical situations due to its good fit to data and simplicity. Many efforts were documented in the literature aimed at understanding the mechanisms that may support Gompertz's elegant model equation. One of the most convincing efforts was carried out by Gyllenberg and Webb. They divide the cancer cell population into the proliferative cells and the quiescent cells. In their two dimensional model, the dead cells are assumed to be removed from the tumor instantly. In this paper, we modify their model by keeping track of the dead cells remaining in the tumor. We perform mathematical and computational studies on this three dimensional model and compare the model dynamics to that of the model of Gyllenberg and Webb. Our mathematical findings suggest that if an avascular tumor grows according to our three-compartment model, then as the death rate of quiescent cells decreases to zero, the percentage of proliferative cells also approaches to zero. Moreover, a slow dying quiescent population will increase the size of the tumor. On the other hand, while the tumor size does not depend on the dead cell removal rate, its early and intermediate growth stages are very sensitive to it.

Tumor evolution in space: the effects of competition colonization tradeoffs on tumor invasion dynamics

We apply competition colonization tradeoff models to tumor growth and invasion dynamics to explore the hypothesis that varying selection forces will result in predictable phenotypic differences in cells at the tumor invasive front compared to those in the core. Spatially, ecologically, and evolutionarily explicit partial differential equation models of tumor growth confirm that spatial invasion produces selection pressure for motile phenotypes. The effects of the invasive phenotype on normal adjacent tissue determine the patterns of growth and phenotype distribution. If tumor cells do not destroy their environment, colonizer and competitive phenotypes coexist with the former localized at the invasion front and the latter, to the tumor interior. If tumors cells do destroy their environment, then cell motility is strongly selected resulting in accelerated invasion speed with time. Our results suggest that the widely observed genetic heterogeneity within cancers may not be the stochastic effect of random mutations. Rather, it may be the consequence of predictable variations in environmental selection forces and corresponding phenotypic adaptations.

Hypoxic cell waves around necrotic cores in glioblastoma: a biomathematical model and its therapeutic implications

Glioblastoma is a rapidly evolving high-grade astrocytoma that is distinguished pathologically from lower grade gliomas by the presence of necrosis and microvascular hyperplasia. Necrotic areas are typically surrounded by hypercellular regions known as "pseudopalisades" originated by local tumor vessel occlusions that induce collective cellular migration events. This leads to the formation of waves of tumor cells actively migrating away from central hypoxia. We present a mathematical model that incorporates the interplay among two tumor cell phenotypes, a necrotic core and the oxygen distribution. Our simulations reveal the formation of a traveling wave of tumor cells that reproduces the observed histologic patterns of pseudopalisades. Additional simulations of the model equations show that preventing the collapse of tumor microvessels leads to slower glioma invasion, a fact that might be exploited for therapeutic purposes.

In silico modelling of tumour margin diffusion and infiltration: review of current status

As a result of advanced treatment techniques, requiring precise target definitions, a need for more accurate delineation of the Clinical Target Volume (CTV) has arisen. Mathematical modelling is found to be a powerful tool to provide fairly accurate predictions for the Microscopic Extension (ME) of a tumour to be incorporated in a CTV. In general terms, biomathematical models based on a sequence of observations or development of a hypothesis assume some links between biological mechanisms involved in cancer development and progression to provide quantitative or qualitative measures of tumour behaviour as well as tumour response to treatment. Generally, two approaches are taken: deterministic and stochastic modelling. In this paper, recent mathematical models, including deterministic and stochastic methods, are reviewed and critically compared. It is concluded that stochastic models are more promising to provide a realistic description of cancer tumour behaviour due to being intrinsically probabilistic as well as discrete, which enables incorporation of patient-specific biomedical data such as tumour heterogeneity and anatomical boundaries.

Accurate state estimation from uncertain data and models: an application of data assimilation to mathematical models of human brain tumors

Background

Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.

Results

We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.

Conclusions

The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.

Reviewers

This article was reviewed by Anthony Almudevar, Tomas Radivoyevitch, and Kristin Swanson (nominated by Georg Luebeck).

Integrative physical oncology

Cancer is arguably the ultimate complex biological system. Solid tumors are microstructured soft matter that evolves as a consequence of spatio-temporal events at the intracellular (e.g., signaling pathways, macromolecular trafficking), intercellular (e.g., cell-cell adhesion/communication), and tissue (e.g., cell-extracellular matrix interactions, mechanical forces) scales. To gain insight, tumor and developmental biologists have gathered a wealth of molecular, cellular, and genetic data, including immunohistochemical measurements of cell type-specific division and death rates, lineage tracing, and gain-of-function/loss-of-function mutational analyses. These data are empirically extrapolated to a diagnosis/prognosis of tissue-scale behavior, e.g., for clinical decision. Integrative physical oncology (IPO) is the science that develops physically consistent mathematical approaches to address the significant challenge of bridging the nano (nm)-micro (µm) to macro (mm, cm) scales with respect to tumor development and progression. In the current literature, such approaches are referred to as multiscale modeling. In the present article, we attempt to assess recent modeling approaches on each separate scale and critically evaluate the current 'hybrid-multiscale' models used to investigate tumor growth in the context of brain and breast cancers. Finally, we provide our perspective on the further development and the impact of IPO.

Bright solitary waves in malignant gliomas

We put forward a nonlinear wave model describing the fundamental dynamical features of an aggressive type of brain tumors. Our model accounts for the invasion of normal tissue by a proliferating and propagating rim of active glioma cancer cells in the tumor boundary and the subsequent formation of a necrotic core. By resorting to numerical simulations, phase space analysis, and exact solutions we prove that bright solitary tumor waves develop in such systems. Possible implications of our model as a tool to extract relevant patient specific tumor parameters in combination with standard preoperative clinical imaging are also discussed.

Tumor-immune interaction, surgical treatment, and cancer recurrence in a mathematical model of melanoma

Malignant melanoma is a cancer of the skin arising in the melanocytes. We present a mathematical model of melanoma invasion into healthy tissue with an immune response. We use this model as a framework with which to investigate primary tumor invasion and treatment by surgical excision. We observe that the presence of immune cells can destroy tumors, hold them to minimal expansion, or, through the production of angiogenic factors, induce tumorigenic expansion. We also find that the tumor–immune system dynamic is critically important in determining the likelihood and extent of tumor regrowth following resection. We find that small metastatic lesions distal to the primary tumor mass can be held to a minimal size via the immune interaction with the larger primary tumor. Numerical experiments further suggest that metastatic disease is optimally suppressed by immune activation when the primary tumor is moderately, rather than minimally, metastatic. Furthermore, satellite lesions can become aggressively tumorigenic upon removal of the primary tumor and its associated immune tissue. This can lead to recurrence where total cancer mass increases more quickly than in primary tumor invasion, representing a clinically more dangerous disease state. These results are in line with clinical case studies involving resection of a primary melanoma followed by recurrence in local metastases.

Melanoma is a deadly skin cancer that invades into the dermis and metastasizes into the surrounding tissue. In clinical cases, surgical excision of the primary tumor has led to widespread and accelerated growth in metastases. We develop a mathematical model describing the basic process of melanoma invasion, metastatic spread, and the anti-tumor immune response. This model is formulated using partial differential equations that describe the spatial and temporal evolution of a number of different cellular populations, and it uses a realistic skin geometry. Using simulations, we examine the importance of the immune response when a primary tumor is spawning satellite metastases. We find that local metastases can be suppressed by the immune response directed against the primary tumor, but grow aggressively following surgical treatment. We also find that moderately metastatic tumors optimally activate the local immune response against disseminated disease, and in this case tumor excision may have profound effects on metastatic growth. We conclude that surgical perturbation of the immune response controlling local metastases is one mechanism by which cancer can recur. This could have implications as to the appropriate clinical management of melanomas and other solid tumors.