In this review we summarize our recent efforts in trying to understand the role of heterogeneity in cancer progression by using neural networks to characterise different aspects of the mapping from a cancer cells genotype and environment to its pheno
type. Our central premise is that cancer is an evolving system subject to mutation and selection, and the primary conduit for these processes to occur is the cancer cell whose behaviour is regulated on multiple biological scales. The selection pressure is mainly driven by the microenvironment that the tumour is growing in and this acts directly upon the cell phenotype. In turn, the phenotype is driven by the intracellular pathways that are regulated by the genotype. Integrating all of these processes is a massive undertaking and requires bridging many biological scales (i.e. genotype, pathway, phenotype and environment) that we will only scratch the surface of in this review. We will focus on models that use neural networks as a means of connecting these different biological scales, since they allow us to easily create heterogeneity for selection to act upon and importantly this heterogeneity can be implemented at different biological scales. More specifically, we consider three different neural networks that bridge different aspects of these scales and the dialogue with the micro-environment, (i) the impact of the micro-environment on evolutionary dynamics, (ii) the mapping from genotype to phenotype under drug-induced perturbations and (iii) pathway activity in both normal and cancer cells under different micro-environmental conditions.
Cancer cells are known to modify their micro-environment such that it can sustain a larger population, or, in ecological terms, they construct a niche which increases the carrying capacity of the population. It has however been argued that niche cons
truction, which benefits all cells in the tumour, would be selected against since cheaters could reap the benefits without paying the cost. We have investigated the impact of niche specificity on tumour evolution using an individual based model of breast tumour growth, in which the carrying capacity of each cell consists of two components: an intrinsic, subclone-specific part and a contribution from all neighbouring cells. Analysis of the model shows that the ability of a mutant to invade a resident population depends strongly on the specificity. When specificity is low selection is mostly on growth rate, while high specificity shifts selection towards increased carrying capacity. Further, we show that the long-term evolution of the system can be predicted using adaptive dynamics. By comparing the results from a spatially structured vs. well-mixed population we show that spatial structure restores selection for carrying capacity even at zero specificity, which a poses solution to the niche construction dilemma. Lastly, we show that an expanding population exhibits spatially variable selection pressure, where cells at the leading edge exhibit higher growth rate and lower carrying capacity than those at the centre of the tumour.
In this paper we analyse a previously proposed cell-based model of glioblastoma (brain tumour) growth, which is based on the assumption that the cancer cells switch phenotypes between a proliferative and motile state (Gerlee and Nelander, PLoS Comp.
Bio., 8(6) 2012). The dynamics of this model can be described by a system of partial differential equations, which exhibits travelling wave solutions whose wave speed depends crucially on the rates of phenotypic switching. We show that under certain conditions on the model parameters, a closed form expression of the wave speed can be obtained, and using singular perturbation methods we also derive an approximate expression of the wave front shape. These new analytical results agree with simulations of the cell-based model, and importantly show that the inverse relationship between wave front steepness and speed observed for the Fisher equation no longer holds when phenotypic switching is considered.
Mathematical modeling in cancer has been growing in popularity and impact since its inception in 1932. The first theoretical mathematical modeling in cancer research was focused on understanding tumor growth laws and has grown to include the competit
ion between healthy and normal tissue, carcinogenesis, therapy and metastasis. It is the latter topic, metastasis, on which we will focus this short review, specifically discussing various computational and mathematical models of different portions of the metastatic process, including: the emergence of the metastatic phenotype, the timing and size distribution of metastases, the factors that influence the dormancy of micrometastases and patterns of spread from a given primary tumor.
In this paper we consider chemotherapy in a spatial model of tumor growth. The model, which is of reaction-diffusion type, takes into account the complex interactions between the tumor and surrounding stromal cells by including densities of endotheli
al cells and the extra-cellular matrix. When no treatment is applied the model reproduces the typical dynamics of early tumor growth. The initially avascular tumor reaches a diffusion limited size of the order of millimeters and initiates angiogenesis through the release of vascular endothelial growth factor (VEGF) secreted by hypoxic cells in the core of the tumor. This stimulates endothelial cells to migrate towards the tumor and establishes a nutrient supply sufficient for sustained invasion. To this model we apply cytostatic treatment in the form of a VEGF-inhibitor, which reduces the proliferation and chemotaxis of endothelial cells. This treatment has the capability to reduce tumor mass, but more importantly, we were able to determine that inhibition of endothelial cell proliferation is the more important of the two cellular functions targeted by the drug. Further, we considered the application of a cytotoxic drug that targets proliferating tumor cells. The drug was treated as a diffusible substance entering the tissue from the blood vessels. Our results show that depending on the characteristics of the drug it can either reduce the tumor mass significantly or in fact accelerate the growth rate of the tumor. This result seems to be due to complicated interplay between the stromal and tumor cell types and highlights the importance of considering chemotherapy in a spatial context.
