Modelling the radiotherapy effect in the reaction-diffusion equation

Applications of Machine Learning (ML) and Mathematical Modeling (MM) in Healthcare with Special Focus on Cancer Prognosis and Anticancer Therapy: Current Status and Challenges

The use of data-driven high-throughput analytical techniques, which has given rise to computational oncology, is undisputed. The widespread use of machine learning (ML) and mathematical modeling (MM)-based techniques is widely acknowledged. These two approaches have fueled the advancement in cancer research and eventually led to the uptake of telemedicine in cancer care. For diagnostic, prognostic, and treatment purposes concerning different types of cancer research, vast databases of varied information with manifold dimensions are required, and indeed, all this information can only be managed by an automated system developed utilizing ML and MM. In addition, MM is being used to probe the relationship between the pharmacokinetics and pharmacodynamics (PK/PD interactions) of anti-cancer substances to improve cancer treatment, and also to refine the quality of existing treatment models by being incorporated at all steps of research and development related to cancer and in routine patient care. This review will serve as a consolidation of the advancement and benefits of ML and MM techniques with a special focus on the area of cancer prognosis and anticancer therapy, leading to the identification of challenges (data quantity, ethical consideration, and data privacy) which are yet to be fully addressed in current studies.

WAM to SeeSaw model for cancer therapy - overcoming LQM difficulties

PURPOSE

The assessment of biological effects caused by radiation exposure has been currently carried out with the linear-quadratic (LQ) model as an extension of the linear non-threshold (LNT) model. In this study, we suggest a new mathematical model named as SeaSaw (SS) model, which describes proliferation and cell death effects by taking account of Bergonie-Tribondeau's law in terms of a differential equation in time. We show how this model overcomes the long-standing difficulties of the LQ model.

MATERIALS AND METHODS

We construct the SS model as an extended Wack-A-Mole (WAM) model by using a differential equation with respect to time in order to express the dynamics of the proliferation effect. A large number of accumulated data of such parameters as α and β in the LQ based models provide us with valuable pieces of information on the corresponding parameter b 1 and the maximum volume V m of the SS model. The dose rate b 1 and the notion of active cell can explain the present data without introduction of β, which is obtained by comparing the SS model with not only the cancer therapy data but also with in vitro experimental data. Numerical calculations are presented to grasp the global features of the SS model.

RESULTS

The SS model predicts the time dependence of the number of active- and inactive-cells. The SS model clarifies how the effect of radiation depends on the cancer stage at the starting time in the treatment. Further, the time dependence of the tumor volume is calculated by changing individual dose strength, which results in the change of the irradiation duration for the same effect. We can consider continuous irradiation in the SS model with interesting outcome on the time dependence of the tumor volume for various dose rates. Especially by choosing the value of the dose rate to be balanced with the total growth rate, the tumor volume is kept constant.

CONCLUSIONS

The SS model gives a simple equation to study the situation of clinical radiation therapy and risk estimation of radiation. The radiation parameter extracted from the cancer therapy is close to the value obtained from animal experiment in vitro and in vivo. We expect the SS model leads us to a unified description of radiation therapy and protection and provides a great development in cancer-therapy clinical-planning.

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.

The biology and mathematical modelling of glioma invasion: a review

Adult gliomas are aggressive brain tumours associated with low patient survival rates and limited life expectancy. The most important hallmark of this type of tumour is its invasive behaviour, characterized by a markedly phenotypic plasticity, infiltrative tumour morphologies and the ability of malignant progression from low- to high-grade tumour types. Indeed, the widespread infiltration of healthy brain tissue by glioma cells is largely responsible for poor prognosis and the difficulty of finding curative therapies. Meanwhile, mathematical models have been established to analyse potential mechanisms of glioma invasion. In this review, we start with a brief introduction to current biological knowledge about glioma invasion, and then critically review and highlight future challenges for mathematical models of glioma invasion.

Selection and Validation of Predictive Models of Radiation Effects on Tumor Growth Based on Noninvasive Imaging Data

The use of mathematical and computational models for reliable predictions of tumor growth and decline in living organisms is one of the foremost challenges in modern predictive science, as it must cope with uncertainties in observational data, model selection, model parameters, and model inadequacy, all for very complex physical and biological systems. In this paper, large classes of parametric models of tumor growth in vascular tissue are discussed including models for radiation therapy. Observational data is obtained from MRI of a murine model of glioma and observed over a period of about three weeks, with X-ray radiation administered 14.5 days into the experimental program. Parametric models of tumor proliferation and decline are presented based on the balance laws of continuum mixture theory, particularly mass balance, and from accepted biological hypotheses on tumor growth. Among these are new model classes that include characterizations of effects of radiation and simple models of mechanical deformation of tumors. The Occam Plausibility Algorithm (OPAL) is implemented to provide a Bayesian statistical calibration of the model classes, 39 models in all, as well as the determination of the most plausible models in these classes relative to the observational data, and to assess model inadequacy through statistical validation processes. Discussions of the numerical analysis of finite element approximations of the system of stochastic, nonlinear partial differential equations characterizing the model classes, as well as the sampling algorithms for Monte Carlo and Markov chain Monte Carlo (MCMC) methods employed in solving the forward stochastic problem, and in computing posterior distributions of parameters and model plausibilities are provided. The results of the analyses described suggest that the general framework developed can provide a useful approach for predicting tumor growth and the effects of radiation.

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.