A scalable solver for a stochastic, hybrid cellular automaton model of personalized breast cancer therapy.

Mathematical modeling and simulation is a promising approach to personalized cancer medicine. Yet, the complexity, heterogeneity and multi-scale nature of cancer pose significant computational challenges. Coupling discrete cell-based models with continuous models using hybrid cellular automata (CA) is a powerful approach for mimicking biological complexity and describing the dynamical exchange of information across different scales. However, when clinically relevant cancer portions are taken into account, such models become computationally very expensive. While efficient parallelization techniques for continuous models exist, their coupling with discrete models, particularly CA, necessitates more elaborate solutions. Building upon FEniCS, a popular and powerful scientific computing platform for solving partial differential equations, we developed parallel algorithms to link stochastic CA with differential equations (https://bitbucket.org/HTasken/cansim). The algorithms minimize the communication between processes that share CA neighborhood values while also allowing for reproducibility during stochastic updates. We demonstrated the potential of our solution on a complex hybrid cellular automaton model of breast cancer treated with combination chemotherapy. On a single-core processor, we obtained nearly linear scaling with an increasing problem size, whereas weak parallel scaling showed moderate growth in solving time relative to increase in problem size. Finally, we applied the algorithm to a problem that is 500 times larger than previous work, allowing us to run personalized therapy simulations based on heterogeneous cell density and tumor perfusion conditions estimated from magnetic resonance imaging data on an unprecedented scale.

Apparent diffusion coefficient estimates based on 24 hours tracer movement support glymphatic transport in human cerebral cortex.

The recently proposed glymphatic system suggests that bulk flow is important for clearing waste from the brain, and as such may underlie the development of e.g. Alzheimer's disease. The glymphatic hypothesis is still controversial and several biomechanical modeling studies at the micro-level have questioned the system and its assumptions. In contrast, at the macro-level, there are many experimental findings in support of bulk flow. Here, we will investigate to what extent the CSF tracer distributions seen in novel magnetic resonance imaging (MRI) investigations over hours and days are suggestive of bulk flow as an additional component to diffusion. In order to include the complex geometry of the brain, the heterogeneous CSF flow around the brain, and the transport over the time-scale of days, we employed the methods of partial differential constrained optimization to identify the apparent diffusion coefficient (ADC) that would correspond best to the MRI findings. We found that the computed ADC in the cortical grey matter was 5-26% larger than the ADC estimated with DTI, which suggests that diffusion may not be the only mechanism governing transport.

Toward Personalized Computer Simulation of Breast Cancer Treatment: A Multiscale Pharmacokinetic and Pharmacodynamic Model Informed by Multitype Patient Data.

The usefulness of mechanistic models to disentangle complex multiscale cancer processes, such as treatment response, has been widely acknowledged. However, a major barrier for multiscale models to predict treatment outcomes in individual patients lies in their initialization and parametrization, which needs to reflect individual cancer characteristics accurately. In this study, we use multitype measurements acquired routinely on a single breast tumor, including histopathology, MRI, and molecular profiling, to personalize parts of a complex multiscale model of breast cancer treated with chemotherapeutic and antiangiogenic agents. The model accounts for drug pharmacokinetics and pharmacodynamics. We developed an open-source computer program that simulates cross-sections of tumors under 12-week therapy regimens and used it to individually reproduce and elucidate treatment outcomes of 4 patients. Two of the tumors did not respond to therapy, and model simulations were used to suggest alternative regimens with improved outcomes dependent on the tumor's individual characteristics. It was determined that more frequent and lower doses of chemotherapy reduce tumor burden in a low proliferative tumor while lower doses of antiangiogenic agents improve drug penetration in a poorly perfused tumor. Furthermore, using this model, we were able to correctly predict the outcome in another patient after 12 weeks of treatment. In summary, our model bridges multitype clinical data to shed light on individual treatment outcomes. SIGNIFICANCE: Mathematical modeling is used to validate possible mechanisms of tumor growth, resistance, and treatment outcome.

Variational data assimilation for transient blood flow simulations: Cerebral aneurysms as an illustrative example.

Several cardiovascular diseases are caused from localised abnormal blood flow such as in the case of stenosis or aneurysms. Prevailing theories propose that the development is caused by abnormal wall shear stress in focused areas. Computational fluid mechanics have arisen as a promising tool for a more precise and quantitative analysis, in particular because the anatomy is often readily available even by standard imaging techniques such as magnetic resonance and computed tomography angiography. However, computational fluid mechanics rely on accurate initial and boundary conditions, which are difficult to obtain. In this paper, we address the problem of recovering high-resolution information from noisy and low-resolution physical measurements of blood flow (for example, from phase-contrast magnetic resonance imaging [PC-MRI]) using variational data assimilation based on a transient Navier-Stokes model. Numerical experiments are performed in both 3D (2D space and time) and 4D (3D space and time) and with pulsatile flow relevant for physiological flow in cerebral aneurysms. The results demonstrate that, with suitable regularisation, the model accurately reconstructs flow, even in the presence of significant noise.

In vivo estimation of elastic heterogeneity in an infarcted human heart.

In myocardial infarction, muscle tissue of the heart is damaged as a result of ceased or severely impaired blood flow. Survivors have an increased risk of further complications, possibly leading to heart failure. Material properties play an important role in determining post-infarction outcome. Due to spatial variation in scarring, material properties can be expected to vary throughout the tissue of a heart after an infarction. In this study we propose a data assimilation technique that can efficiently estimate heterogeneous elastic material properties in a personalized model of cardiac mechanics. The proposed data assimilation is tested on a clinical dataset consisting of regional left ventricular strains and in vivo pressures during atrial systole from a human with a myocardial infarction. Good matches to regional strains are obtained, and simulated equi-biaxial tests are carried out to demonstrate regional heterogeneities in stress-strain relationships. A synthetic data test shows a good match of estimated versus ground truth material parameter fields in the presence of no to low levels of noise. This study is the first to apply adjoint-based data assimilation to the important problem of estimating cardiac elastic heterogeneities in 3-D from medical images.