A three-dimensional model of cell movement in multicellular systems☆
Introduction
Cell and tissue movement plays a vital role throughout the lifespan of many organisms. Bacteria and other single-cell organisms find food and avoid repellents by chemotaxis. Coordinated movements of cells and tissues occur throughout early embryonic development, and special terminology such as gastrulation and invagination is used to describe them.
In multicellular systems the relative movement of cells can lead to sorting out of different cell types and formation of specific structures and patterns, such as during gastrulation [1], [10], cancer cell invasion into tissues [39], the development of Dictyostelium discoideum [3], limb regeneration [43], in wound healing [1], and white blood cell movements through blood vessels [23]. In these multicellular systems the collective motion can be quite different than the motion of isolated individuals. Since the combined effect of cell–cell interactions, and production and propagation of chemotactic signals are often very complex, it is not always obvious what underlying mechanism is responsible for the observed patterns. Thus, there is a need to develop models of cell movements in multicellular systems that will give insight into the mechanism governing this complex behavior. There have been several different approaches used to model the collective motion of cells in multicellular systems:
- 1.
A cellular automata approach [34], in which formal rules for cell–cell interactions and for the reactions of the cells to their chemical environment are postulated. Simulations of such models show patterns that resemble experimental observations, but limited insight is gained from such models because they do not incorporate the physics of cell–cell interactions and cell–cell signaling.
- 2.
Models designed to simulate cell sorting [7], [20], [24], [25], [38], [41]. One class are fluid type models where the cell aggregate is modeled as a mixture of two fluids. Here the cells have zero size, so the sorting time only depends on the “temperature” (randomness) of the fluid mixture. The problem with these models is that the size and stiffness of the cells are ignored, so they really verify only the theoretical predictions regarding mixing of two fluids. Another class of models use some kind of Potts model [17], [18], [34] to represent the energy between the cells, and cells exchange position with higher probabilities if that exchange is energetically favorable. The problem with Potts type models is that they fail to explain how the cells exchange positions and in some sense they assume that the cell knows beforehand which positions are energetically favorable.
- 3.
Continuum approaches such as those used by Odell and Bonner [27] and by Vasiev et al. [45] to model the slug movement in Dictyostelium discoideum. In the Odell–Bonner model, cells respond to cAMP chemotactically, as they do during the aggregation phase, and the active component of the propulsive force enters as a contribution to the stress tensor. Active stresses are derived from an active relative velocity field set up at the cell surface. Each amoeba is thought of as a “cortical tractor” pulling the surrounding viscous fluid past itself by rotating its membrane at different rates on different parts of the surface. While this is an important attempt at modeling the motion of the slug, there are limitations to the model. Firstly, it is known that cytoplasm is not a simple Newtonian fluid, as they assume [26], [33]. Secondly, the properties of the extracellular fluid, which are critical for the model, are not known. Thirdly, traction forces play no role in the model. Finally, no molecular or cell-based mechanism is proposed that gives rise to the constitutive relations used in the model, and thus one cannot relate them to properties of the cells. Vasiev et al. incorporate the cAMP dynamics, via a FitzHugh–Nagumo model, into a continuum model of mound formation in Dictyostelium. These authors begin with the Navier–Stokes equations for the velocity field and simply add forces intended to model those arising from chemotaxis. While this approach can lead to solutions that apparently reflect aggregation, it is completely phenomenological and no procedure is given to connect the forces postulated with experimentally measurable quantities such as the force exerted by a single cell.
The model we present here is a three-dimensional mathematical model of cell movement and cell–cell interactions which allows us to analyze and gain insight into the various processes that govern cell movement inside multicellular systems. To make it feasible to model the motion of tissues comprising thousands of cells, we stipulate that all cells are ellipsoids, and thus restrict the admissible deformations to those which only change the relative lengths of the semi-axes but conserve the cell volume. These computational cells have characteristics and responses that correspond to those observed in real cells. These characteristics include the stiffness of the cell, cell adhesion, locomotive force generation and response to environmental cues. Most cell types have common characteristics like, how they move and interact, deform, exert forces onto other cells and move in response to chemical or mechanical cues, but their exact responses differ. For instances the stiffness of a cell can vary by up to two orders of magnitude; a fibroblast has stiffness 0.6 mdyn/μm whereas a red blood cell has a stiffness of 10−3 mdyn/μm [5]. The adhesion can also differ for different type of cells [30], the surface tension for a limb bud aggregate is 20 dyn/cm but it is 1.6 dyn/cm for Neural retina aggregate [15]. The locomotive force a cell applies also varies more than an order of magnitude [19], [21], [22], [28], [42]. In our model, these values can be changed for each cell individually depending on what cell types we are trying to model, thus making the model very flexible, and applicable to a large number of different multicellular systems. Since our model is based on individual cells and the movement and deformation of the cells is calculated directly from all the forces acting on the cells from the equations of motion and deformation, it does not have many of the limitations of the other models mentioned above. This feature allows us to study the effect that finite cell size, the stiffness of the cells and random movement, has on the rate and completeness of sorting.
Before designing the model, it is important to know how an individual cell moves. When an ameoboid cell moves, either randomly or in response to a chemotactic signal, it sends out pseudopods, one of which eventually dominates and attaches it to the surroundings [46], [48]. This determines the direction of motion, as the cell realigns its axis towards the pseudopod, and the rest of the cell body is “pulled” towards the attached pseudopod. The realignment is not achieved by rotation but by disassembly and reassembly of several proteins [43], [44]. The establishment of a new pseudopod and realignment of the axis does take some time [46], [48], but once established the cell becomes polarized and moves in the new direction for up to several minutes; for Dictyostelium new pseudopods form about once every 3 min [44]. When the cell senses a chemical gradient, most of the new pseudopods are formed towards the gradient, but in the absence of a gradient the probability of lateral pseudopods is increased [13], [44]. Once the cell moves in a gradient, it becomes polarized and is even less likely to form lateral pseudopods.
The extensions of pseudopods and the retraction of the rest of the cell body require active force generation, and this force is applied at the site of pseudopod attachment (Fig. 1a). When the cell is moving on a surface, the pseudopod attachment will be onto the surface and the applied force is transmitted directly to the ground. However, when the cell is inside a multicellular aggregate, it must attach the pseudopod to another cell (Fig. 1b). It is important to realize that when a cell attaches its pseudopod to another cell, and pulls itself forward, the other cell experiences an equal force in the opposite direction. So if that cell is not firmly attached to the ground it is pulled backwards; for a good discussion see [27]. Only forces that are transmitted down to the ground give rise to any net movement of the multicellular system and therefore a complete model of cell movements must balance all the forces. Neglecting to balance the forces has been the main fault of many previous models. In 2D it is not as important since one can assume that each cell is in direct contact with the surface where it applies the active force. However this does not hold true in 3D where the active force most often is applied on a neighbor cell, and this force also affects the motion of that cell. These interactions can have a significant effect on the motion of the system as a whole.
Section snippets
Design of model
The basic units in the model are individual cells, each of which is characterized by its location and orientation within the aggregate, its state of stress, and the active forces it can exert in response to the local micro-environment. Knowing this for each cell, the movement of all cells, and hence of an aggregate of cells can be calculated. We describe below how the present state of a cell is determined, how the forces it can exert are computed, and the algorithm for time-stepping the
Results
Malcolm Steinberg’s group performed a number of experiments [9], [14], [15], where they compressed cell aggregates of embryonic tissues between two parallel plates and measured the force required to keep the plates at a given separation. The purpose of these experiments was to show that the cells behave as viscoelastic liquids and to measure the surface tension of these embryonic tissues, and thus get a measure of cell–cell adhesion strength. Foty [15] used the following approximation for the
Discussion
We have developed a model that allows us to simulate and visualize in 3D the movement of cells in multicellular systems. This model allows us to explore how specific characteristics of the cells, such as adhesion and cell stiffness affect the cell motion. We can also use the model to explore how chemotactic signals direct cell movement and how in turn the cell movement affects the signal propagation. Being able to visualize the movement and trajectories of individual cells can help us gain
Supplementary data
Acknowledgements
Supported in part by NSF Grant No. 9805494.
Eirikur Palsson received his MA in Plasma Physics in 1991 and MA in Applied Mathematics in 1993 and his PhD in Applied Mathematics/Molecular Biology in 1996, all from Princeton University. His PhD thesis focused on pattern formation and on the initiation and evolution of spiral and target patterns due to the cAMP signaling in Dictyostelium discoideum. He had a postdoctoral position from 1996 to 1999 in the Department of Mathematics at the University of Utah, where he worked on models of cell
References (50)
- et al.
Surface tensions and tissue affinities in Rana Pipiens gastrulae: quantitative measurements
Dev. Biol.
(1997) - et al.
Cellular mechanism underlying neural convergent extension in Xenopus Laevis embryos
Dev. Biol.
(1997) Detailed mechanics of membrane–membrane adhesion and separation. I. Continuum of molecular cross-bridges
Biophys. J. Biophys. Soc.
(1985)Detailed mechanics of membrane–membrane adhesion and separation. II. Discrete kinetically trapped molecular cross-bridges
Biophys. J. Biophys. Soc.
(1985)- et al.
Viscoelastic properties of living embryonic tissues: a quantitative study
Biophys. J.
(1998) Geometrical models for tissues
Int. Rev. Cytol.
(1983)Leukocite–endothelial cell interactions
Curr. Opin. Cell Biol.
(1992)- et al.
A stochastic model for cell sorting and measuring cell–cell adhesion
J. Theoret. Biol.
(1996) On the elasticity of cytoskeletal networks
Biophys. J.
(1988)- et al.
Analysis of optical density wave propagation and cell movement during mound formation in Dictyostelium discoideum
Dev. Biol.
(1996)
Modeling morphogenesis: from single cells to crawling slugs
J. Theoret. Biol.
Analysis of optical density wave propagation and cell movement in the cellular slime mould Dictyostelium discoideum
Physica D
Adhesion-guided multicellular assembly: a commentary upon the postulates, real and imagined, of the differential adhesion hypothesis, with special attention to computer simulations of cell sorting
J. Theoret. Biol.
A model of cell sorting
J. Theoret. Biol.
Cadherins in cancer: implications for invasion and metastasis
Curr. Opin. Cell Biol.
Theoretical model for morphogenesis and cell sorting in Dictyostelium discoideum
Physica D
Locomotion forces generated by a polymorphonuclear leukocyte
Biophys. J.
A hydrodynamic model for Dictyostelium discoideum mound formation
J. Theoret. Biol.
Models for the specific adhesion of cells to cells
Science
A way of following individual cells in the migration slugs of Dictyostelium discoideum
Proc. Natl. Acad. Sci. USA
Cited by (0)
Eirikur Palsson received his MA in Plasma Physics in 1991 and MA in Applied Mathematics in 1993 and his PhD in Applied Mathematics/Molecular Biology in 1996, all from Princeton University. His PhD thesis focused on pattern formation and on the initiation and evolution of spiral and target patterns due to the cAMP signaling in Dictyostelium discoideum. He had a postdoctoral position from 1996 to 1999 in the Department of Mathematics at the University of Utah, where he worked on models of cell movements in multicellular tissues. Currently, he is an Assistant Professor in the Department of Biology at the College of Staten Island, where he is continuing to develop his model on cell movements. He has been modeling in detail the movement and cell–cell signaling in Dictyostelium discoideum. His research interests pertain to mathematical biology. In particular, designing and programming of mathematical models in molecular and cellular biology and using these models to get an understanding of the underlying biological mechanism. Simulations of his models can be viewed at http://www.math.utah.edu/ ∼ epalsson/research.html.
- ☆
Electronic annexes available. See http://www.elsevier.nl/locate/future.