Vladimir P.
Zhdanov
*ab and
Bengt
Kasemo
a
aDepartment of Applied Physics, Chalmers University of Technology, S-412 96 Göteborg, Sweden
bBoreskov Institute of Catalysis, Russian Academy of Sciences, Novosibirsk, 630090, Russia. E-mail: zhdanov@catalysis.nsk.su; zhdanov@fy.chalmers.se
First published on 1st December 2003
We present Monte Carlo lattice simulations illustrating how the growth of stem and differentiated cells on a uniform surface can be governed (guided differentiation) by tuning long-range cell–cell communication. Due to spatial constraints on cell division, the growth rate is shown to strongly depend on the interplay of cellular division, diffusion, and signaling. With decreasing diffusion rate, the kinetics become slower. To increase the selectivity of the growth with respect to “desirable” cells, one can use induced messenger-mediated generation of these cells by stem cells. The simplest schemes based on this idea may, however, result in appreciable slowdown of the growth. A more optimal strategy may be simultaneous promotion of generation of desirable cells and stem cells.
In addition to the experiments, it is of interest to form a basis of theoretical tissue engineering or, more specifically, to develop models describing the growth of ensembles of stem and differentiated cells. At present such models are in fact lacking (for a review, see ref. 10). In the nearest future, one can hardly expect the models constructed to be fully predictive. However, mathematical modelling in this area is still motivated, firstly because it may help to conceptually understand what may occur in real systems, and secondly it is a concise way of structuring existing knowledge. Thirdly, it is an invaluable tool to guide experiments.
The goal of our present Monte Carlo (MC) simulations is to illustrate how the 2D cellular growth on an uniform surface depend on and/or can be governed by varying the rate of cellular mobility and/or cell–cell communication. (For reviews of cellular mobility and communication, see refs. 11, 12 and 13, respectively.)
In principle, the number of cell types could be easily increased, but we believe that the 3-type case is generic because it makes it possible to illustrate various aspects of cellular growth with a minimal number of model parameters. In particular, a model with this number cell types can be used to show how the selectivity of the growth with respect to desirable cells (e.g. of type 2) depends on the details of diffusion, division, and cell–cell communication.
In the simulations, cells are located on a 2D L×
L square lattice. Each lattice site is either vacant or occupied by one cell. This means that the reconfiguration of a single cell is not described explicitly (for the multi-site models making it possible to mimic changes of the cell shape, see ref. 14). The cell size is assumed to be slightly larger than the site size but smaller than the distance between the next-nearest-neighbour (nnn) sites. With this prescription, the interactions between nearest-neighbour (nn) and nnn cells, εnnij and εnnnij, should be repulsive and attractive due to deformation and adhesion, respectively. In addition, the attractive nnn interactions between similar cell should be stronger than that between different cells. These interactions, introduced to take into account softness, deformation, adhesion and segregation of cells, are of course “effective”.
Cell diffusion is realized via jumps (spreading) to vacant nn sites. A jump is allowed provided that in the final state (after the jump) the cell has at least one contact with other cells, i.e., after a jump at least one of the nn or nnn sites should be occupied (this condition is introduced in order to suppress rare events of escape of single cells from the central area where they are located; its role in the simulations is minor, because the probability of such events is low anyway). The normalized jump probability is defined by the standard Metropolis rule. In the case under consideration, the Metropolis algorithm is appropriate, because it correctly describes the tendency of cells to increase the number of adhesive contacts and segregate. One should however bear in mind that in our context the thermal energy, kBT, used in the Metropolis rule, is effective as well as the cell–cell interactions (for discussion of this aspect of the lattice models, see ref. 15).
In our previous general simulations,10 we allowed division of cells of all the types. Specifically, the division rate of differentiated cells was assumed to be higher that that of the stem cells. From the point of view of stem-cell tissue engineering, it is of interest to analyze the opposite case, because if the rate of division of differentiated cells were high, there would be no need of stem cells. For this reason, we prohibit division of differentiated cells. The division of a stem cell is allowed provided that (i) all the nn sites are vacant (the cells located in nn sites are deformed and accordingly the division probability in such cases is assumed to be negligibly small) and (ii) the newly-born cell located (after division) in one of the nn sites has at least one nn vacant site.
In reality, the rate of division of stem cells and, what is more important, the differentiation probabilities depend on cell–cell signaling including transport of messenger molecules between cells, signal transmission via the cell membrane, gene regulation, etc. The signaling may occur via ligand diffusion between cells12 or via the so-called gap-junction communication.13 In the former case, the cells may be far from each other. In the latter case, the cells should be in close contact. Earlier,10 we discussed both these schemes, but the simulations were performed for the simplest case of communication between adjacent (nnn) cells. The present simulations are focused on a more general case of long-range ligand-diffusion-mediated signaling.
The cell–cell communication is assumed to influence primarily the “decision-taking” process in stem cells, i.e. it changes the distribution of the newly born cells but does not change the division rate. Thus, the effect of signaling on the rate of division of these cells is neglected. The differentiation probabilities of stem cells are considered to depend on the signals arriving from cells of type 2 or, more specifically, on the signal molecule concentration due to signaling of type-2 cells (the signaling of cells of types 1 and 3 is for simplicity neglected). The signal molecule concentration near a given stem cell is represented as10
![]() | (1) |
![]() | (2) |
The probabilities Pdivi specify the distribution of cells after division but do not characterize the rate of division. To run MC simulations, we actually need a ratio of the rates of cell division and diffusion. To characterize the relative rates of these processes, we use a dimensionless parameter pdiv. During the simulations, the division and diffusion trials are performed with probabilities pdiv and 1−
pdiv, respectively. In all the examples below, diffusion of cells is considered to be faster or much faster than division, i.e.pdiv
≪
1. In reality, diffusion of cells on surfaces is usually slow, but division is slow as well. Thus, the condition pdiv
≪
1 may be fulfilled.
(1) A site is chosen at random. If the site is vacant, the trial ends. Otherwise, a cell located in the site has chances to diffuse irrespective of its type [item (2)] if ρ>
pdiv, or to divide [item (3)] provided that it belongs to type 1 and ρ
<
pdiv, where ρ
(0
≤
ρ
≤
1) is a random number.
(2) For diffusion, an adjacent site is randomly selected, and if the latter site is vacant, the cell is transferred to it provided that (i) in the final state the cell has at least one contact with other cells and (ii)
ρ′<
Pdifi, where ρ′
(0
≤
ρ′
≤
1) is a new random number, and Pdifi is the standard Metropolis jump probability dependent on cell–cell lateral interactions.
(3) Division of a stem cell may occur if all the nn sites are vacant. In this case, one of the nn sites (site 2) is selected at random and if this site has at least one vacant nn site the division is realized with unit probability. After division of a stem cell, there are two cells. One of them is a stem cell and another one is a stem or differentiated cell (the type of this cell is defined by probabilities Pdivi). The former cell occupies the site which was occupied by the stem cell before division. The latter cell is situated in site 2.
All the MC runs were started from the clean lattice with a single stem cell located in the center. During simulations, we used the so-called “no-flux” boundary conditions (this means that a MC trial ends if the nn site chosen for diffusion jump or division is located outside the L×
L array of sites). The duration of the runs was usually selected so that in the end there were about 2000 cells.
Time was measured in MC steps (MCS). Usually, one MCS is defined as L×
L MC trials. In our case, this definition would result in a time scale connected primarily with cell migration, because this process is considered to be rapid compared to cell division. Under such circumstances, the conventionally defined MC time is physically somewhat artificial. The situation here is similar to MF calculations of reaction kinetics, where the time scale is usually independent of the rate of reactant diffusion. For this reason, we define one MC step as L
×
L MC trials of cell division on average. Practically, this means that the MC time is calculated by dividing the total number of trials by L
×
L and multiplying by pdiv. One of the advantages of this choice is that we can directly, with the same MC time, compare the kinetics corresponding to very different rates of cell mobility.
The simulations presented are executed on a lattice with L=
150. The energetic parameters are selected so that different cells tend to segregate and this process is not too slow. Practically, this means that the ratio of the effective cell–cell interactions and kBT should not be too large. In particular, we use εnnij/kBT
=
1, εnnn11/kBT
=
−1, εnnnii/kBT
=
−2 for i
≥
2, and εnnnij/kBT
=
−1 for i
≠
j. The stem-cell division probability, pdiv, is chosen to be a governing parameter and varied from 10−3 to 0.1 (this means that the ratio of the rates of cell diffusion and division is between one and three orders of magnitude). For single stem cells, the distribution of birth probabilities for the three cell types is considered to be P01
=
0.1, P02
=
0.5, and P03
=
0.4. To describe cell–cell communication, we measure distance in terms of the lattice spacing and employ α
=
0.1 (this means that the length scale of communication is about ten lattice spacings). The parameters Bi
[in eqn. (2)] are varied in order to illustrate typical situations.
![]() | ||
Fig. 1 Number of cells of types 1, 2, and 3 (thin solid, thick solid, and dashed lines) as a function of time for pdiv![]() ![]() |
![]() | ||
Fig. 2 [(a)–(c)] Patterns obtained in the end of the MC runs shown in Fig. 1 (a)–(c). Open circles, filled circles, and plus signs represent cells of types 1, 2, and 3, respectively. |
Typical results of the simulations for negligible cell–cell communication (Bi=
0) are shown in Fig. 1. For all the values of pdiv, the characteristics of the cellular growth kinetics are similar. The distribution of cell types is the same as well, because in the absence of cell signaling it is defined simply by probabilities P0i. The times, characterizing the ensemble growth, and cellular patterns are however different. If cell mobility is relatively low (pdiv
=
0.1), the growth is slow (the 2000 cell population is reached after about 700 MCS). With increasing diffusion rate compared to rate of division, the growth becomes much faster. For pdiv
=
10−3, e.g., the 2000 cell population is observed already after about 80 MCS. This difference in the growth rate is related to the spatial constraints on cell division and also related to cell segregation. If cell diffusion is relatively slow (or, in other words, division is rapid), the cellular patterns are rather compact [Fig. 2 (a) and (b)], segregation is poor, stem cells are surrounded by other cells, the conditions for division of stem cells inside the cellular patch are as a rule not fulfilled, and accordingly the growth rate is relatively low. With increasing diffusion rate, the cell arrangement is less compact [Fig. 2(c)], segregation is better, and accordingly the kinetics are faster.
The simplest cell–cell communication favourable for production of cells of type 2 includes signal-molecule-mediated promotion of birth of these cells by stem cells due to signals from type-2 cells. Mathematically, this situation is described by using B1=
B3
=
0 and B2
>
0. Typical results of the corresponding simulations, performed with B2
=
0.1 (this value is relatively small), are exhibited in Figs. 3 and 4. As expected, the self-promotion of type-2 cells is seen to result in an increase of their population (cf.Figs. 1 and 3). The price for this effect is that the kinetics become slower, especially for slow diffusion [cf.Figs. 1(a) and 3(a)], because type-2 cells suppress birth of stem cells. With increase of B2 up to 0.2 (Figs. 5 and 6), the population of type-2 cells further increases but not appreciably. The slowdown of the cellular growth for slow diffusion is in this case dramatic [cf.Figs. 1(a), 3(a), and 5(a)].
![]() | ||
Fig. 3 Number of cells of type 1, 2, and 3 (thin solid, thick solid, and dashed lines) as a function of time for pdiv![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
![]() | ||
Fig. 4 [(a)–(c)] Snapshots obtained in the end of the MC runs presented in Fig. 3(a)–(c). The designations are as in Fig. 2. |
![]() | ||
Fig. 5 As Fig. 3 for B2![]() ![]() |
![]() | ||
Fig. 6 [(a)–(c)] Patterns obtained in the end of the MC runs shown in Fig. 5(a)–(c). The designations are as in Fig. 2. |
To avoid the slowdown in growth, cell–cell communication can be tuned so that in addition to self-promotion the type-2 cells promote birth of stem cells. This situation is illustrated by employing B1=
0.05, B1
=
0.05, and B3
=
0 (Figs. 7 and 8). In this case, the growth kinetics are rapid even if diffusion is relatively slow [Fig. 7(a)]. The price for this improvement is that the selectivity of production of type-2 cells becomes somewhat lower [cf.Figs. 5(a) and 7(a)].
![]() | ||
Fig. 7 As Fig. 3 for B1![]() ![]() ![]() ![]() |
![]() | ||
Fig. 8 [(a)–(c)] Snapshots obtained in the end of the MC runs presented in Fig. 7(a)–(c). The designations are as in Fig. 2. |
This journal is © the Owner Societies 2004 |