A Generalised Age- and Phase-Structured Model of Human Tumour Cell Populations Both Unperturbed and Exposed to a Range of Cancer Therapies

We develop a general mathematical model for a population of cells differentiated by their position within the cell division cycle. A system of partial differential equations governs the kinetics of cell densities in certain phases of the cell division cycle dependent on time t (hours) and an age-like variable τ (hours) describing the time since arrival in a particular phase of the cell division cycle. Transition rate functions control the transfer of cells between phases. We first obtain a theoretical solution on the infinite domain −∞ < t < ∞. We then assume that age distributions at time t=0 are known and write our solution in terms of these age distributions on t=0. In practice, of course, these age distributions are unknown. All is not lost, however, because a cell line before treatment usually lies in a state of asynchronous balanced growth where the proportion of cells in each phase of the cell cycle remain constant. We assume that an unperturbed cell line has four distinct phases and that the rate of transition between phases is constant within a short period of observation (‘short’ relative to the whole history of the tumour growth) and we show that under certain conditions, this is equivalent to exponential growth or decline. We can then gain expressions for the age distributions. So, in short, our approach is to assume that we have an unperturbed cell line on t ≤ 0, and then, at t=0 the cell line is exposed to cancer therapy. This corresponds to a change in the transition rate functions and perhaps incorporation of additional phases of the cell cycle. We discuss a number of these cancer therapies and applications of the model.