No Arabic abstract
We consider a mathematical model for the evolutionary dynamics of tumour cells in vascularised tumours under chemotherapy. The model comprises a system of coupled partial integro-differential equations for the phenotypic distribution of tumour cells, the concentration of oxygen and the concentration of a chemotherapeutic agent. In order to disentangle the impact of different evolutionary parameters on the emergence of intra-tumour phenotypic heterogeneity and the development of resistance to chemotherapy, we construct explicit solutions to the equation for the phenotypic distribution of tumour cells and provide a detailed quantitative characterisation of the long-time asymptotic behaviour of such solutions. Analytical results are integrated with numerical simulations of a calibrated version of the model based on biologically consistent parameter values. The results obtained provide a theoretical explanation for the observation that the phenotypic properties of tumour cells in vascularised tumours vary with the distance from the blood vessels. Moreover, we demonstrate that lower oxygen levels may correlate with higher levels of phenotypic variability, which suggests that the presence of hypoxic regions supports intra-tumour phenotypic heterogeneity. Finally, the results of our analysis put on a rigorous mathematical basis the idea, previously suggested by formal asymptotic results and numerical simulations, that hypoxia favours the selection for chemoresistant phenotypic variants prior to treatment. Consequently, this facilitates the development of resistance following chemotherapy.
We present a mathematical study of the emergence of phenotypic heterogeneity in vascularised tumours. Our study is based on formal asymptotic analysis and numerical simulations of a system of non-local parabolic equations that describes the phenotypic evolution of tumour cells and their nonlinear dynamic interactions with the oxygen, which is released from the intratumoural vascular network. Numerical simulations are carried out both in the case of arbitrary distributions of intratumour blood vessels and in the case where the intratumoural vascular network is reconstructed from clinical images obtained using dynamic optical coherence tomography. The results obtained support a more in-depth theoretical understanding of the eco-evolutionary process which underpins the emergence of phenotypic heterogeneity in vascularised tumours. In particular, our results offer a theoretical basis for empirical evidence indicating that the phenotypic properties of cancer cells in vascularised tumours vary with the distance from the blood vessels, and establish a relation between the degree of tumour tissue vascularisation and the level of intratumour phenotypic heterogeneity.
Tumour cells have to acquire a number of capabilities if a neoplasm is to become a cancer. One of these key capabilities is increased motility which is needed for invasion of other tissues and metastasis. This paper presents a qualitative mathematical model based on game theory and computer simulations using cellular automata. With this model we study the circumstances under which mutations that confer increased motility to cells can spread through a tumour made of rapidly proliferating cells. The analysis suggests therapies that could help prevent the progression towards malignancy and invasiveness of benign tumours.
In nature and human societies, the effects of homogeneous and heterogeneous characteristics on the evolution of collective behaviors are quite different from each other. It is of great importance to understand the underlying mechanisms of the occurrence of such differences. By incorporating pair pattern strategies and reference point strategies into an agent-based model, we have investigated the coupled effects of heterogeneous investment strategies and heterogeneous risk tolerance on price fluctuations. In the market flooded with the investors with homogeneous investment strategies or homogeneous risk tolerance, large price fluctuations are easy to occur. In the market flooded with the investors with heterogeneous investment strategies or heterogeneous risk tolerance, the price fluctuations are suppressed. For a heterogeneous population, the coexistence of investors with pair pattern strategies and reference point strategies causes the price to have a slow fluctuation around a typical equilibrium point and both a large price fluctuation and a no-trading state are avoided, in which the pair pattern strategies push the system far away from the equilibrium while the reference point strategies pull the system back to the equilibrium. A theoretical analysis indicates that the evolutionary dynamics in the present model is governed by the competition between different strategies. The strategy that causes large price fluctuations loses more while the strategy that pulls the system back to the equilibrium gains more. Overfrequent trading does harm to ones pursuit for more wealth.
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 endothelial 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.
Background: Analysing tumour architecture for metastatic potential usually focuses on phenotypic differences due to cellular morphology or specific genetic mutations, but often ignore the cells position within the heterogeneous substructure. Similar disregard for local neighborhood structure is common in mathematical models. Methods: We view the dynamics of disease progression as an evolutionary game between cellular phenotypes. A typical assumption in this modeling paradigm is that the probability of a given phenotypic strategy interacting with another depends exclusively on the abundance of those strategies without regard local heterogeneities. We address this limitation by using the Ohtsuki-Nowak transform to introduce spatial structure to the go vs. grow game. Results: We show that spatial structure can promote the invasive (go) strategy. By considering the change in neighbourhood size at a static boundary -- such as a blood-vessel, organ capsule, or basement membrane -- we show an edge effect that allows a tumour without invasive phenotypes in the bulk to have a polyclonal boundary with invasive cells. We present an example of this promotion of invasive (EMT positive) cells in a metastatic colony of prostate adenocarcinoma in bone marrow. Interpretation: Pathologic analyses that do not distinguish between cells in the bulk and cells at a static edge of a tumour can underestimate the number of invasive cells. We expect our approach to extend to other evolutionary game models where interaction neighborhoods change at fixed system boundaries.