A cellular automaton model for tumour growth in inhomogeneous environment

Most of the existing mathematical models for tumour growth and tumour-induced angiogenesis neglect blood flow. This is an important factor on which both nutrient and metabolite supply depend. In this paper we aim to address this shortcoming by developing a mathematical model which shows how blood flow and red blood cell heterogeneity influence the growth of systems of normal and cancerous cells. The model is developed in two stages. First we determine the distribution of oxygen in a native vascular network, incorporating into our model features of blood flow and vascular dynamics such as structural adaptation, complex rheology and red blood cell circulation. Once we have calculated the oxygen distribution, we then study the dynamics of a colony of normal and cancerous cells, placed in such a heterogeneous environment. During this second stage, we assume that the vascular network does not evolve and is independent of the dynamics of the surrounding tissue. The cells are considered as elements of a cellular automaton, whose evolution rules are inspired by the different behaviour of normal and cancer cells. Our aim is to show that blood flow and red blood cell heterogeneity play major roles in the development of such colonies, even when the red blood cells are flowing through the vasculature of normal, healthy tissue

Lattice and Continuum Modelling of a Bioactive Porous Tissue Scaffold

A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport, and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold, and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models, and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.Comment: 38 pages, 16 figures. This version includes a much-expanded introduction, and a new section on nonlinear diffusion in addition to polish throughou

Using Cellular Automata and Lattice Boltzmann Methods to Model Cancer Growth: Analysis of Combination Treatment Outcomes

In Canada it is estimated that 76,600 people will die of cancer in 2014. Cancer, a collection of over 200 diseases, has differences existing between globally, between individuals and overtime in one individual. Treatment options are similarly varied. These differences make selecting the best possible treatment for every type of cancer very challenging. In addition, with no single cure for cancer, treatments are often combined in different ways to form the best overall option. In an attempt to synthesize the properties of these diseases into a collection of common cellular changes, Hanahan and Weinberg proposed ``the hallmarks of cancer -- 10 differences between healthy cells and cancer cells, present in almost every cancer. There exists the potential for treatments that are broadly applicable if they reverse these general properties. This work seeks to simulate early cancer growth, specifically looking at these hallmarks, and detect the best combinations of hallmarks to remove in order to stop cancer growth. This hybrid simulation combines a discrete model of cancer cells using cellular automata, with a continuous model of blood flow using lattice Boltzmann methods. Hallmarks relevant during the early growth stages of solid tumour development are simulated using rules in the cellular automata. Hallmarks were removed in pairs, triplets and quadruplets in order to model combination therapy, abstracting drugs that target these properties as the removal of the hallmark from the system. Overall growth of the tumours with ``treatments applied were compared to tumours where all hallmarks were present. It was found that many combinations had no effect on tumour growth. In some cases combinations even increased growth, selecting for the most aggressive hallmarks since weaker hallmarks were unavailable. However, in general, as more treatments were applied, cancer growth decreased. This work is the first to simulate removing hallmarks in pairs, triplets and quadruplets from a model with biologically relevant oxygen flow. It provides a proof of concept that not all combinations are equally effective, even if the individual treatments are effective. This work suggests some combinations should be avoided while others could potentially be beneficial in a variety of diseases

A Three Species Model to Simulate Application of Hyperbaric Oxygen Therapy to Chronic Wounds

Chronic wounds are a significant socioeconomic problem for governments worldwide. Approximately 15% of people who suffer from diabetes will experience a lower-limb ulcer at some stage of their lives, and 24% of these wounds will ultimately result in amputation of the lower limb. Hyperbaric Oxygen Therapy (HBOT) has been shown to aid the healing of chronic wounds; however, the causal reasons for the improved healing remain unclear and hence current HBOT protocols remain empirical. Here we develop a three-species mathematical model of wound healing that is used to simulate the application of hyperbaric oxygen therapy in the treatment of wounds. Based on our modelling, we predict that intermittent HBOT will assist chronic wound healing while normobaric oxygen is ineffective in treating such wounds. Furthermore, treatment should continue until healing is complete, and HBOT will not stimulate healing under all circumstances, leading us to conclude that finding the right protocol for an individual patient is crucial if HBOT is to be effective. We provide constraints that depend on the model parameters for the range of HBOT protocols that will stimulate healing. More specifically, we predict that patients with a poor arterial supply of oxygen, high consumption of oxygen by the wound tissue, chronically hypoxic wounds, and/or a dysfunctional endothelial cell response to oxygen are at risk of nonresponsiveness to HBOT. The work of this paper can, in some way, highlight which patients are most likely to respond well to HBOT (for example, those with a good arterial supply), and thus has the potential to assist in improving both the success rate and hence the cost-effectiveness of this therapy

Optimisation of cancer drug treatments using cell population dynamics

International audienceCancer is primarily a disease of the physiological control on cell population proliferation. Tissue proliferation relies on the cell division cycle: one cell becomes two after a sequence of molecular events that are physiologically controlled at each step of the cycle at so-called checkpoints, in particular at transitions between phases of the cycle [105]. Tissue proliferation is the main physiological process occurring in development and later in maintaining the permanence of the organism in adults, at that late stage mainly in fast renewing tissues such as bone marrow, gut and skin

Cancer modeling: From mechanistic to data-driven approaches, and from fundamental insights to clinical applications

peer reviewe

The role of acidity in tumour development

Acidic pH is a common characteristic of human tumours. It has a significant impact on tumour progression and response to therapies. In this thesis, we utilise mathematical modelling to examine the role of acidosis in the interaction between normal and tumour cell populations. In the first section we investigate the cell–microenvironmental interactions that mediate somatic evolution of cancer cells. The model predicts that selective forces in premalignant lesions act to favour cells whose metabolism is best suited to respond to local changes in oxygen, glucose and pH levels. In particular the emergent cellular phenotype, displaying increased acid production and resistance to acid-induced toxicity, has a significant proliferative advantage because it will consistently acidify the local environment in a way that is toxic to its competitors but harmless to itself. In the second section we analyse the role of acidity in tumour growth. Both vascular and avascular tumour dynamics are investigated, and a number of different behaviours are observed. Whilst an avascular tumour always proceeds to a benign steady state, a vascular tumour may display either benign or invasive dynamics, depending on the value of a critical parameter. Extensions of the model show that cellular quiescence, or non-proliferation, may provide an explanation for experimentally observed cycles of acidity within tumour tissue. Analysis of both models allows assessment of novel therapies directed towards changing the level of acidity within the tumour. Finally we undertake a comparison between experimental tumour pH images and the models of acid dynamics set out in previous chapters. This analysis will allow us to assess and verify the previous modelling work, giving the mathematics a firm biological foundation. Moreover, it provides a methodology of calculating important diagnostic parameters from pH images

Exploring Complex Networks with Graph Investigator Research Application

This paper describes Graph Investigator, the application intended for analysis of complex networks. A rich set of application functions is briefly described including graph feature generation, comparison, visualization and edition. The program enables to analyze global and local structural properties of networks with the use of various descriptors derived from graph theory. Furthermore, it allows to quantify inter-graph similarity by embedding graph patterns into low-dimensional space or distance measurement based on feature vectors. The set of available graph descriptors includes over eighty statistical and algebraic measures. We present two examples of real-world networks analysis performed with Graph Investigator: comparison of brain vasculature with structurally similar artificial networks and analysis of vertices importance in a macaque cortical connectivity network. The third example describes tracking parameters of artificial vascular network evolving in the process of angiogenesis, modelled with the use of cellular automata