Global analysis of a cancer model with drug resistance due to Lamarckian induction and microvesicle transfer

Abstract

Development of resistance to chemotherapy in cancer patients strongly effects the outcome of the treatment. Due to chemotherapeutic agents, resistance can emerge by Darwinian evolution. Besides this, acquired drug resistance may arise via changes in gene expression. A recent discovery in cancer research uncovered a third possibility, indicating that this phenotype conversion can occur through the transfer of microvesicles from resistant to sensitive cells, a mechanism resembling the spread of an infectious agent. We present a model describing the evolution of sensitive and resistant tumour cells considering Darwinian selection, Lamarckian induction and microvesicle transfer. We identify three threshold parameters which determine the existence and stability of the three possible equilibria. Using a simple Dulac function, we give a complete description of the dynamics of the model depending on the three threshold parameters.We also establish an agent based model as a spatial version of the ODE model and compare the outputs of the two models. We find that although the ODE model does not provide spatial information about the structure of the tumour, it is capable to determine the outcome in terms of tumour size and distribution of cell types.We demonstrate the possible effects of increasing drug concentration, and characterize the possible bifurcation sequences. Our results show that the presence of microvesicle transfer cannot ruin a therapy that otherwise leads to extinction, however it may doom a partially successful therapy to failure.