Carina M. Edwards

Mathematical modeling of cell population dynamics in the colonic crypt and in colorectal cancer

Colorectal cancer is initiated in colonic crypts. A succession of genetic mutations or epigenetic changes can lead to homeostasis in the crypt being overcome, and subsequent unbounded growth. We consider the dynamics of a single colorectal crypt by using a compartmental approach [Tomlinson IPM, Bodmer WF (1995) Proc Natl Acad Sci USA 92:11130–11134], which accounts for populations of stem cells, differentiated cells, and transit cells. That original model made the simplifying assumptions that each cell population divides synchronously, but we relax these assumptions by adopting an age-structured approach that models asynchronous cell division, and by using a continuum model. We discuss two mechanisms that could regulate the growth of cell numbers and maintain the equilibrium that is normally observed in the crypt. The first will always maintain an equilibrium for all parameter values, whereas the second can allow unbounded proliferation if the net per capita growth rates are large enough. Results show that an increase in cell renewal, which is equivalent to a failure of programmed cell death or of differentiation, can lead to the growth of cancers. The second model can be used to explain the long lag phases in tumor growth, during which new, higher equilibria are reached, before unlimited growth in cell numbers ensues.

On the proportion of cancer stem cells in a tumour

It is now generally accepted that cancers contain a sub-population, the cancer stem cells (CSCs), which initiate and drive a tumours growth. At least until recently it has been widely assumed that only a small proportion of the cells in a tumour are CSCs. Here we use a mathematical model, supported by experimental evidence, to show that such an assumption is unwarranted. We show that CSCs may comprise any possible proportion of the tumour, and that the higher the proportion the more aggressive the tumour is likely to be.

Force localization in contracting cell layers.

Epithelial cell layers on soft elastic substrates or pillar arrays are commonly used as model systems for investigating the role of force in tissue growth, maintenance, and repair. Here we show analytically that the experimentally observed localization of traction forces to the periphery of the cell layers does not necessarily imply increased local cell activity, but follows naturally from the elastic problem of a finite-sized contractile layer coupled to an elastic foundation. For homogeneous contractility, the force localization is determined by one dimensionless parameter interpolating between linear and exponential force profiles for the extreme cases of very soft and very stiff substrates, respectively. If contractility is sufficiently increased at the periphery, outward directed displacements can occur at intermediate positions. We also show that anisotropic extracellular stiffness can lead to force localization in the stiffer direction, as observed experimentally.

Examples of Mathematical Modeling: Tales from the Crypt

Mathematical modeling is being increasingly recognized within the biomedical sciences as an important tool that can aid the understanding of biological systems. The heavily regulated cell renewal cycle in the colonic crypt provides a good example of how modeling can be used to find out key features of the system kinetics, and help to explain both the breakdown of homeostasis and the initiation of tumorigenesis.We use the cell population model by Johnston et al. (2007) Proc. Natl. Acad. Sci. USA 104, 4008-4013, to illustrate the power of mathematical modeling by considering two key questions about the cell population dynamics in the colonic crypt. We ask: how can a model describe both homeostasis and unregulated growth in tumorigenesis; and to which parameters in the system is the model most sensitive? In order to address these questions, we discuss what type of modeling approach is most appropriate in the crypt.We use the model to argue why tumorigenesis is observed to occur in stages with long lag phases between periods of rapid growth, and we identify the key parameters.

Comparing a discrete and continuum model of the intestinal crypt.

The integration of processes at different scales is a key problem in the modelling of cell populations. Owing to increased computational resources and the accumulation of data at the cellular and subcellular scales, the use of discrete, cell-level models, which are typically solved using numerical simulations, has become prominent. One of the merits of this approach is that important biological factors, such as cell heterogeneity and noise, can be easily incorporated. However, it can be difficult to efficiently draw generalizations from the simulation results, as, often, many simulation runs are required to investigate model behaviour in typically large parameter spaces. In some cases, discrete cell-level models can be coarse-grained, yielding continuum models whose analysis can lead to the development of insight into the underlying simulations. In this paper we apply such an approach to the case of a discrete model of cell dynamics in the intestinal crypt. An analysis of the resulting continuum model demonstrates that there is a limited region of parameter space within which steady-state (and hence biologically realistic) solutions exist. Continuum model predictions show good agreement with corresponding results from the underlying simulations and experimental data taken from murine intestinal crypts.