In-silico oncology: an approximate model of brain tumor mass effect based on directly manipulated free form deformation

Brain shift in neuronavigation of brain tumors: A review

PURPOSE

Neuronavigation based on preoperative imaging data is a ubiquitous tool for image guidance in neurosurgery. However, it is rendered unreliable when brain shift invalidates the patient-to-image registration. Many investigators have tried to explain, quantify, and compensate for this phenomenon to allow extended use of neuronavigation systems for the duration of surgery. The purpose of this paper is to present an overview of the work that has been done investigating brain shift.

METHODS

A review of the literature dealing with the explanation, quantification and compensation of brain shift is presented. The review is based on a systematic search using relevant keywords and phrases in PubMed. The review is organized based on a developed taxonomy that classifies brain shift as occurring due to physical, surgical or biological factors.

RESULTS

This paper gives an overview of the work investigating, quantifying, and compensating for brain shift in neuronavigation while describing the successes, setbacks, and additional needs in the field. An analysis of the literature demonstrates a high variability in the methods used to quantify brain shift as well as a wide range in the measured magnitude of the brain shift, depending on the specifics of the intervention. The analysis indicates the need for additional research to be done in quantifying independent effects of brain shift in order for some of the state of the art compensation methods to become useful.

CONCLUSION

This review allows for a thorough understanding of the work investigating brain shift and introduces the needs for future avenues of investigation of the phenomenon.

Identification of crucial parameters in a mathematical multiscale model of glioblastoma growth

Glioblastomas are highly malignant brain tumours. Mathematical models and their analysis provide a tool to support the understanding of the development of these tumours as well as the design of more effective treatment strategies. We have previously developed a multiscale model of glioblastoma progression that covers processes on the cellular and molecular scale. Here, we present a novel nutrient-dependent multiscale sensitivity analysis of this model that helps to identify those reaction parameters of the molecular interaction network that influence the tumour progression on the cellular scale the most. In particular, those parameters are identified that essentially determine tumour expansion and could be therefore used as potential therapy targets. As indicators for the success of a potential therapy target, a deceleration of the tumour expansion and a reduction of the tumour volume are employed. From the results, it can be concluded that no single parameter variation results in a less aggressive tumour. However, it can be shown that a few combined perturbations of two systematically selected parameters cause a slow-down of the tumour expansion velocity accompanied with a decrease of the tumour volume. Those parameters are primarily linked to the reactions that involve the microRNA-451 and the thereof regulated protein MO25.

A generic framework for modeling brain deformation as a constrained parametric optimization problem to aid non-diffeomorphic image registration in brain tumor imaging

OBJECTIVES

In the present paper a novel computational framework for modeling tumor induced brain deformation as a biophysical prior for non-rigid image registration is described. More precisely, we aim at providing a generic building block for non-rigid image registration that can be used to resolve inherent irregularities in non-diffeomorphic registration problems that naturally arise in serial and cross-population brain tumor imaging studies due to the presence (or progression) of pathology.

METHODS

The model for the description of brain cancer dynamics on a tissue level is based on an initial boundary value problem (IBVP). The IBVP follows the accepted assumption that the progression of primary brain tumors on a tissue level is governed by proliferation and migration of cancerous cells into surrounding healthy tissue. The model of tumor induced brain deformation is phrased as a parametric, constrained optimization problem. As a basis of comparison and to demonstrate generalizability additional soft constraints (penalties) are considered. A back-tracking line search is implemented in conjunction with a limited memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) method in order to handle the numerically delicate log-barrier strategy for confining volume change.

RESULTS

Numerical experiments are performed to test the flexible control of the computed deformation patterns in terms of varying model parameters. The results are qualitatively and quantitatively related to patterns in patient individual magnetic resonance imaging data.

CONCLUSIONS

Numerical experiments demonstrate the flexible control of the computed deformation patterns. This in turn strongly suggests that the model can be adapted to patient individual imaging patterns of brain tumors. Qualitative and quantitative comparison of the computed cancer profiles to patterns in medical imaging data of an exemplary patient demonstrates plausibility. The designed optimization problem is based on computational tools widely used in non-rigid image registration, which in turn makes the model generally applicable for integration into non-rigid image registration algorithms.

Biophysical modeling of brain tumor progression: from unconditionally stable explicit time integration to an inverse problem with parabolic PDE constraints for model calibration

PURPOSE

A novel unconditionally stable, explicit numerical method is introduced to the field of modeling brain cancer progression on a tissue level together with an inverse problem (IP) based on optimal control theory that allows for automated model calibration with respect to observations in clinical imaging data.

METHODS

Biophysical models of cancer progression on a tissue level are in general based on the assumption that the spatiotemporal spread of cancerous cells is determined by cell division and net migration. These processes are typically described in terms of a parabolic partial differential equation (PDE). In the present work a parallelized implementation of an unconditionally stable, explicit Euler (EE(⋆)) time integration method for the solution of this PDE is detailed. The key idea of the discussed EE(⋆) method is to relax the strong stability requirement on the spectral radius of the coefficient matrix by introducing a subdivision regime for a given outer time step. The performance is related to common implicit numerical methods. To quantify the numerical error, a simplified model that has a closed form solution is considered. To allow for a systematic, phenomenological validation a novel approach for automated model calibration on the basis of observations in medical imaging data is developed. The resulting IP is based on optimal control theory and manifests as a large scale, PDE constrained optimization problem.

RESULTS

The numerical error of the EE(⋆) method is at the order of standard implicit numerical methods. The computing times are well below those obtained for implicit methods and by that demonstrate efficiency. Qualitative and quantitative analysis in 12 patients demonstrates that the obtained results are in strong agreement with observations in medical imaging data. Rating simulation success in terms of the mean overlap between model predictions and manual expert segmentations yields a success rate of 75% (9 out of 12 patients).

CONCLUSIONS

The discussed EE(⋆) method provides desirable features for image-based model calibration or hybrid image registration algorithms in which the model serves as a biophysical prior. This is due to (i) ease of implementation, (ii) low memory requirements, (iii) efficiency, (iv) a straightforward interface for parameter updates, and (v) the fact that the method is inherently matrix-free. The explicit time integration method is confirmed via experiments for automated model calibration. Qualitative and quantitative analysis demonstrates that the proposed framework allows for recovering observations in medical imaging data and by that phenomenological model validity.

A novel method for simulating the extracellular matrix in models of tumour growth

A novel hybrid continuum-discrete model to simulate tumour growth on a cellular scale is proposed. The lattice-based spatiotemporal model consists of reaction-diffusion equations that describe interactions between cancer cells and their microenvironment. The fundamental ingredients that are typically considered are the nutrient concentration, the extracellular matrix (ECM), and matrix degrading enzymes (MDEs). The in vivo processes are very complex and occur on different levels. This in turn leads to huge computational costs. The main contribution of the present work is therefore to describe the processes on the basis of simplified mathematical approaches, which, at the same time, depict realistic results to understand the biological processes. In this work, we discuss if we have to simulate the MDE or if the degraded matrix can be estimated directly with respect to the cancer cell distribution. Additionally, we compare the results for modelling tumour growth using the common and our simplified approach, thereby demonstrating the advantages of the proposed method. Therefore, we introduce variations of the positioning of the nutrient delivering blood vessels and use different initializations of the ECM. We conclude that the novel method, which does not explicitly model the matrix degrading enzymes, provides means for a straightforward and fast implementation for modelling tumour growth.