CPM Model Results

1 Cell

Resting Area of 4 Pixels

The simulation ran for 100 different instances each for a maximum of 60 Monte Carlo (MC) time steps.

From the simulations we were able to get statistics on the mean square displacement (MSD) and the displacement of the cell.

For 2-dimensional diffusion the MSD should obey the following relationship

Where is the diffusion coefficient and should be equal to 1.

Results

1cell_msd

So fitting our whole data we get the following values for and (in simulation units).

This is in relatively good agreement with expectations and with the result from the Szabo paper which states that displacement grows with time as

We can get closer agreement if we ignore the outlying, first three data points from our fit.

1cell_d_2

Which yields…

Resting Area of 50 Pixels

Results

Here is a movie of a typical cell.

1cell_msd_1

The results look good! The single cell with a relaxed area of 50 pixels does exhibit normal diffusion. We can also compare our results to that of the Szabo et al paper to see how our simulation compares. To do this one MC time step becomes one minute and the distance between one lattice site becomes 5.

We can translate our data into these units and plot the mean displacement.

1cell_msd_2

This appears to be in good agreement with the results from the paper. At 10 hours the mean displacement is around which is inline with the results from the paper.

I attempted to overlay this plot onto that from the paper in order to see how well the data matched up.

szabo_comparison

As you can see the data from the simulations (bottom blue) is very similar to the black curve from the paper. The slopes seem to be in agreement but there is an off set which corresponds to a difference in diffusion coefficients. I believe that some of the off set is due my plot not exactly lining up with how the plot from the paper scales but this still that does not explain away the whole off set.

Multiple Cells

In the case of multiple cells the relative values of the energy terms associated with free-cell and cell-cell boundaries are important.

Please be aware of the difference between this use of and that in the context of diffusion.

In the Szabo paper and . However, in the paper’s simulations of cell monolayers there are no free-cell boundaries so the parameter becomes irrelevant. The paper reports that increasing causes movement of cells within the monolayer to decrease. This makes sense since increasing causes the energy due to changing the cell boundaries to increase making movement less probable.

I believe that mimicking the values of and from the paper is not productive since we do want to consider a group of cells with free-cell boundaries. If then cell-cell boundaries are not favorable which will cause the group of cells to break apart. Therefore, we should restrict the parameters to in order to maintain cell-cell contact.

Here is a video for the case of and .

So I have chosen the following values for the parameters.

This satisfies the criteria for favorable cell-cell contact, and also keeps the same value of from single cell simulations.

Below is a a video for the case of and .

Results - Updated ( 6/3 )

Here are updated results with tracking of the whole group of cells center of mass properly implemented.

As before, and and a relaxed area of 4 pixels.

The results shown will be for the case of and and a relaxed area of 4 pixels. Below are the results for the MSD for groups of various number of cells.

small_mlt_com_msd_1

We can get a better idea for what is going on if we plot the fits for each respective data set.

small_mlt_com_msd_2

Notice that the groups move slower as the number of cells in a group increases. The relationship between MSD and time remains linear as the number of cells varies.

Lets see what happens to the diffusion coefficient as the group size changes.

small_mlt_com_diff

The overall behavior is that as the group size increases the diffusion coefficient decreases and it seems to be approaching some lower bound.

Old Results ( 5/ 22 )

small_mltcells_msd_1

We can get a better for what is going on if we plot the fits for each respective data set.

small_mltcells_msd_2

small_mltcells_msd_3

As you can see as the number of cells in a group increases so does . So as the number of cells increases diffusion of the group becomes more anomalous.

Lets see what happens to the diffusion coefficient as the group size changes.

small_mltcells_diff

The overall behavior is that as the group size increases the diffusion coefficient decreases and it seems to be approaching some lower bound. So even though increases this does not lead to faster diffusion for larger groups since the diffusion coefficient decreases.

Notes on Calculating Diffusion Coefficient and

Assume that the mean square displacement (MSD) follows the relationship

The simulations give us data for as a function of the number of time steps taken . Taking the logarithm of the data obtained we should be able to fit it to a linear function.

Julien Varennes

Written on May 20, 2015