Modelling cell guidance and curvature control in evolving biological tissues.

Tissue geometry is an important influence on the evolution of many biological tissues. The local curvature of an evolving tissue induces tissue crowding or spreading, which leads to differential tissue growth rates, and to changes in cellular tension, which can influence cell behaviour. Here, we investigate how directed cell motion interacts with curvature control in evolving biological tissues. Directed cell motion is involved in the generation of angled tissue growth and anisotropic tissue material properties, such as tissue fibre orientation. We develop a new cell-based mathematical model of tissue growth that includes both curvature control and cell guidance mechanisms to investigate their interplay. The model is based on conservation principles applied to the density of tissue synthesising cells at or near the tissue's moving boundary. The resulting mathematical model is a partial differential equation for cell density on a moving boundary, which is solved numerically using a hybrid front-tracking method called the cell-based particle method. The inclusion of directed cell motion allows us to model new types of biological growth, where tangential cell motion is important for the evolution of the interface, or for the generation of anisotropic tissue properties. We illustrate such situations by applying the model to simulate both the resorption and infilling components of the bone remodelling process, and to simulate root hair growth. We also provide user-friendly MATLAB code to implement the algorithms.

Nonlinear studies of tumor morphological stability using a two-fluid flow model.

We consider the nonlinear dynamics of an avascular tumor at the tissue scale using a two-fluid flow Stokes model, where the viscosity of the tumor and host microenvironment may be different. The viscosities reflect the combined properties of cell and extracellular matrix mixtures. We perform a linear morphological stability analysis of the tumors, and we investigate the role of nonlinearity using boundary-integral simulations in two dimensions. The tumor is non-necrotic, although cell death may occur through apoptosis. We demonstrate that tumor evolution is regulated by a reduced set of nondimensional parameters that characterize apoptosis, cell-cell/cell-extracellular matrix adhesion, vascularization and the ratio of tumor and host viscosities. A novel reformulation of the equations enables the use of standard boundary integral techniques to solve the equations numerically. Nonlinear simulation results are consistent with linear predictions for nearly circular tumors. As perturbations develop and grow, the linear and nonlinear results deviate and linear theory tends to underpredict the growth of perturbations. Simulations reveal two basic types of tumor shapes, depending on the viscosities of the tumor and microenvironment. When the tumor is more viscous than its environment, the tumors tend to develop invasive fingers and a branched-like structure. As the relative ratio of the tumor and host viscosities decreases, the tumors tend to grow with a more compact shape and develop complex invaginations of healthy regions that may become encapsulated in the tumor interior. Although our model utilizes a simplified description of the tumor and host biomechanics, our results are consistent with experiments in a variety of tumor types that suggest that there is a positive correlation between tumor stiffness and tumor aggressiveness.

A level-set method for the evolution of cells and tissue during curvature-controlled growth.

Most biological tissues grow by the synthesis of new material close to the tissue's interface, where spatial interactions can exert strong geometric influences on the local rate of growth. These geometric influences may be mechanistic or cell behavioural in nature. The control of geometry on tissue growth has been evidenced in many in vivo and in vitro experiments, including bone remodelling, wound healing, and tissue engineering scaffolds. In this paper, we propose a generalisation of a mathematical model that captures the mechanistic influence of curvature on the joint evolution of cell density and tissue shape during tissue growth. This generalisation allows us to simulate abrupt topological changes such as tissue fragmentation and tissue fusion, as well as three dimensional cases, through a level-set-based method. The level-set method developed introduces another Eulerian field than the level-set function. This additional field represents the surface density of tissue-synthesising cells, anticipated at future locations of the interface. Numerical tests performed with this level-set-based method show that numerical conservation of cells is a good indicator of simulation accuracy, particularly when cusps develop in the tissue's interface. We apply this new model to several situations of curvature-controlled tissue evolutions that include fragmentation and fusion.