Computer simulations of cell sorting due to differential adhesion

The actions of cell adhesion molecules, in particular, cadherins during embryonic development and morphogenesis more generally, regulate many aspects of cellular interactions, regulation and signaling. Often, a gradient of cadherin expression levels drives collective and relative cell motions generating macroscopic cell sorting. Computer simulations of cell sorting have focused on the interactions of cells with only a few discrete adhesion levels between cells, ignoring biologically observed continuous variations in expression levels and possible nonlinearities in molecular binding. In this paper, we present three models relating the surface density of cadherins to the net intercellular adhesion and interfacial tension for both discrete and continuous levels of cadherin expression. We then use then the Glazier-Graner-Hogeweg (GGH) model to investigate how variations in the distribution of the number of cadherins per cell and in the choice of binding model affect cell sorting. We find that an aggregate with a continuous variation in the level of a single type of cadherin molecule sorts more slowly than one with two levels. The rate of sorting increases strongly with the interfacial tension, which depends both on the maximum difference in number of cadherins per cell and on the binding model. Our approach helps connect signaling at the molecular level to tissue-level morphogenesis.