Eur. Phys. J. E (2013) 36: 65
https://doi.org/10.1140/epje/i2013-13065-4

Regular Article

Model on cell movement, growth, differentiation and de-differentiation: Reaction-diffusion equation and wave propagation

¹⁹⁸⁷⁸ Department of Physics, Graduate Institute of Biophysics, and Center for Complex Systems, National Central University, Chungli, Taiwan, 320, R.O.C.
²⁹⁸⁷⁸ Institute of Physics, Academia Sinica, Nankang, Taipei, Taiwan, 115, R.O.C.
³⁹⁸⁷⁸ School of Science, Nanjing University of Science and Technology, Nanjing, 210094, China

^* e-mail: pylai@phy.ncu.edu.tw
^** e-mail: ckchan@gate.sinica.edu.tw

Received: 11 September 2012
Revised: 17 April 2013
Accepted: 4 June 2013
Published online: 27 June 2013

Abstract

We construct a model for cell proliferation with differentiation into different cell types, allowing backward de-differentiation and cell movement. With different cell types labeled by state variables, the model can be formulated in terms of the associated transition probabilities between various states. The cell population densities can be described by coupled reaction-diffusion partial differential equations, allowing steady wavefront propagation solutions. The wavefront profile is calculated analytically for the simple pure growth case (2-states), and analytic expressions for the steady wavefront propagating speeds and population growth rates are obtained for the simpler cases of 2-, 3- and 4-states systems. These analytic results are verified by direct numerical solutions of the reaction-diffusion PDEs. Furthermore, in the absence of de-differentiation, it is found that, as the mobility and/or self-proliferation rate of the down-lineage descendant cells become sufficiently large, the propagation dynamics can switch from a steady propagating wavefront to the interesting situation of propagation of a faster wavefront with a slower waveback. For the case of a non-vanishing de-differentiation probability, the cell growth rate and wavefront propagation speed are both enhanced, and the wavefront speeds can be obtained analytically and confirmed by numerical solution of the reaction-diffusion equations.