First Passage Time Simulations

Phase Space

We vary the parameters and and observe the resulting mean first passage time. Other parameters are left fixed unless otherwise noted, all plots are of groups of 4 cells.

No Perimeter Energy Constraint

no_lambdaP

Above is the resulting phase space for simulations in which there is no constraint on the perimeter length of cells. Since there is no perimeter constraint, cells are allowed to stretch into long-stringy shapes which are not biologically relevant.

From the plot we see that FPT is at a minimum for a ratio of and equal to a little more than 1 . If then cell polarization strength is very strong and works to pull the group of cells apart causing FPT to increase. In the case of the energy cost of the cell moving into the ECM is too large resulting in very slow net movement of the group.

Perimeter Energy Constraint

Although we can optimize parameters to get short FPT, in all cases the cells tend to stretch-out into unphysical shapes. In order to prevent this from happening we add a perimeter energy contraint term to the Hamiltonian of the form

This will make it energetically favorable for cells to retain circle-like shapes.

Strong Perimeter Constraint

In this case we have set .

lambdaP_0.1_scaled

The FPT data has been scaled such that 100,000 time steps is the max value represented in the heat map.

In this case it seems that there are two separate regions of minimal FPT. One occurs at low values of for both parameters. The other is again for ratio a of and equal to a little more than 1, but only for sufficiently large values of both parameters.

Lets explore what the cell dynamics actually look like for different regions on this phase space. I hope you like pink!

For and

For and

For and

Weak Perimeter Constraint

In this case we set , ten times smaller than the previous case.

lambdaP_0.01

The resulting phase space looks much more like the original, no perimeter constraint result. Again, for a ratio of to of a little more than one the FPT is minimized. Like the other simulations, when is too large the FPT increases since movement into the ECM becomes energetically unfavorable. When is too large the group of cells is torn apart decreasing cell’s ability to sense its environment.

Lets look at the cell dynamics for different parameter values in the phase space.

For and

For and

For and

Julien Varennes

Written with StackEdit.

Written on September 3, 2015