Abstract

In this paper we build on the mathematical model of Ward and King (1998) to study the effects of high molecular mass mitotic inhibitors released at cell death. The model assumes a continuum of living cells which, depending on the concentration of a generic nutrient, generate movement (described by a velocity field) due to the changes in volumes caused by cell birth and death. The necrotic material is assumed to consist of two diffusible materials: 1) basic cellular material which is used by living cells as raw material for mitosis; 2) a generic non-utilisable material which may inhibit mitosis. Numerical solutions of the resulting system of partial differential equations show all the main features of tumour growth and heterogeneity. Material 2) is found to act in an inhibitive fashion in two ways: i) directly, by reducing the mitotic rate and ii) indirectly, by occupying space, thereby reducing the availability of the basic cellular material. For large time the solutions to the model tend either to a steady-state, reflecting growth saturation, or to a travelling wave, indicating continual linear growth. The steady-state and travelling wave limits of the model are derived and studied, the regions of existence of these two types of long-time solution being explored in parameter space using numerical methods.