University of Groningen A Model for a Correlated Random Walk Based on the Ordered Extension of Pseudopodia

Cell migration in the absence of external cues is well described by a correlated random walk. Most single cells move by extending protrusions called pseudopodia. To deduce how cells walk, we have analyzed the formation of pseudopodia by Dictyostelium cells. We have observed that the formation of pseudopodia is highly ordered with two types of pseudopodia: First, de novo formation of pseudopodia at random positions on the cell body, and therefore in random directions. Second, pseudopod splitting near the tip of the current pseudopod in alternating right/left directions, leading to a persistent zig-zag trajectory. Here we analyzed the probability frequency distributions of the angles between pseudopodia and used this information to design a stochastic model for cell movement. Monte Carlo simulations show that the critical elements are the ratio of persistent splitting pseudopodia relative to random de novo pseudopodia, the Left/Right alternation, the angle between pseudopodia and the variance of this angle. Experiments confirm predictions of the model, showing reduced persistence in mutants that are defective in pseudopod splitting and in mutants with an irregular cell surface. Citation: Van Haastert PJM (2010) A Model for a Correlated Random Walk Based on the Ordered Extension of Pseudopodia. PLoS Comput Biol 6(8): e1000874. doi:10.1371/journal.pcbi.1000874 Editor: Edward C. Cox, Princeton University, United States of America Received April 14, 2010; Accepted June 30, 2010; Published August 12, 2010 Copyright: 2010 Peter J. M. Van Haastert. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: The author received no specific funding for this work. Competing Interests: The author has declared that no competing interests exist. * E-mail: P.J.M.van.haastert@rug.nl


Introduction
Eukaryotic cells move by extending pseudopodia, which are actin-filled protrusions of the cell surface [1]. Pseudopod formation by Dictyostelium cells, like many other moving cells, shows a typical pseudopod cycle: upon their initiation, pseudopodia grow at a constant rate during their first ,15 s and then stop. The next pseudopod is typically formed a few seconds later, but sometimes commences while the present pseudopod is still growing, giving rise to a cell with two pseudopodia. The fate of the pseudopod after its initial growth phase determines its role in cell movement: the pseudopod is either retracted, or is maintained by flow of the cytoplasm into the pseudopod thereby moving the cell body. The frequency, position and directions of the maintained pseudopodia form the basis of cell movement, because they determine the speed and trajectory of the cell. An important aspect of cell motility is the ability of cells to respond to directional cues with oriented movement. Gradients of chemicals give rise to chemotaxis [2]. Other directional cues that can induce oriented movement are temperature gradients (thermotaxis) or electric fields (electrotaxis) [3,4]. These signals somehow modulate basal pseudopod extension such that, on average, cells move in the direction of the positional cues. In this respect, studies on cell movement are critical for understanding directional movement.
Cells in the absence of external cues do not move in random directions but exhibit a so-called correlated random walk [5][6][7][8][9]. This tendency to move in the same direction is called persistence, and the duration of the correlation is the persistence time. Cells with strong persistence make fewer turns, move for prolonged periods of time in the same direction, and thereby effectively penetrate into the surrounding space. Other search strategies for efficient exploration are local diffusive search and Levi walks [8,10]. Can we understand the cell trajectory by analyzing how cells extend pseudopodia?
To obtain large data sets of extending pseudopodia we developed a computer algorithm that identifies the cell contour and its protrusions. The extending pseudopod is characterized by a vector that connects the x,y,t coordinates of the pseudopod at the beginning and end of the growth phase, respectively [11]. A picture of ordered cell movement has emerged from the analysis of ,6000 pseudopodia that are extended by wild type and mutant cells in buffer [12]. Dictyostelium cells, as many other eukaryotic cells, may extend two types of pseudopodia: de novo at regions devoid of recent pseudopod activity, or by splitting of an existing pseudopod [12,13]. Pseudopod splitting occurs very frequently alternating to the right and left at a relatively small angle of ,55 degrees. Therefore, pseudopod splitting may lead to a persistent zig-zag trajectory [14]. In contrast, de novo pseudopodia are extended in all directions and do not exhibit a right/left bias, suggesting that de novo pseudopodia induce a random turn of the cells. We observed strong persistence for cells that extend many splitting pseudopodia. Conversely, mutants that extend mostly de novo pseudopodia have very short persistence time and exhibit a near Brownian random walk [12].
In this report we investigated the theory of correlated random walks in the context of the observed ordered extension of pseudopodia. The aim is to define the descriptive persistence time or average turn angle with primary experimentally-derived pseudopod properties. First we obtained detailed quantitative data on the probability frequency distributions of the size and direction of pseudopod activity. We then formulated a model that consists of five components: pseudopod size, fraction of splitting pseudopodia, alternating right/left bias, angle between pseudopodia and variance of this angle due to irregularity of cell shape. We measured the parameter values of these components for several Dictyostelium mutants with defects in signaling pathways or cytoskeleton functions. Subsequently, we used these observed parameters in Monte Carlo simulations of the model and compared the predicted trajectories with the observed trajectories of the mutants. The results demonstrate two critical components in these correlated random walks: the ratio of pseudopod splitting relative to de novo pseudopodia, and the shape of the cell.

Dictyostelium strains and cell recordings
The strains used are wild type AX3, pi3k-null strain GMP1 with a deletion of pi3k1 and pi3k2 genes [15], pla2-null with a deletion of the plaA gene [16], sgc/gca-null cells (abbreviated as gc-null cells) with a deletion of gca and sgc genes [17], sgc/pla2-null cells with a deletion of sgc and pla2A genes [18], and ddia2-null cells lacking the forH gene encoding the Dictyostelium homologue of formin [19]. Cells were grown in HG5 medium (contains per liter: 14.3 g oxoid peptone, 7.15 g bacto yeast extract, 1.36 g Na 2 HPO 4 ?12H 2 O, 0.49 g KH 2 PO 4 , 10.0 g glucose), harvested in PB (10 mM KH2PO4/Na2HPO4, pH 6.5), and allowed to develop in 1 ml PB in a well of a 6-wells plate (Nunc). Movies were recorded at a rate of 1 frame per second for at least 15 minutes with an inverted light microscope (Olympus Type CK40 with 206 objective) and images were captured with a JVC CCD camera. Cell trajectories were recorded as the movement of the centroid of the cell as described [20].

Pseudopod analysis
Images were analyzed with the fully automatic pseudopodtracking algorithm Quimp3, which is described in detail [11]. Briefly, the program uses an active contour analysis to represent the outline of the cell using ,150 nodes [21]. Extending pseudopodia that satisfied the user-defined minimum number of adjacent convex nodes and the minimum area change were identified. The direction of each extending pseudopod was identified by the x,y and time coordinates of the central convex node of the convex area at the start and end of growth, respectively. The tangent to the surface at the node where the pseudopod started was calculated using the position of the adjacent nodes. The automated algorithm annotates each pseudopod as de novo versus splitting using the criterion that the convex area of the new pseudopod exhibits overlap with the convex area of the current pseudopod or is within a user-defined distance. The output files containing the x,y-coordinates of the start and end position of the pseudopod, the tangent of the surface and the annotation of the pseudopod were imported in Excel to calculate pseudopod size, interval, direction to gradient, direction to tangent, etc for de novo and splitting pseudopodia (see Fig. 1), as well as fraction s of pseudopod splitting and alternating Right/Left bias a (RL +LR)/total splitting; Table 1).
The cell shape parameter Y was determined as follows: Using the outline of the cell with ,150 nodes, two ellipsoids were constructed, the largest ellipse inside the cell outline and the smallest ellipse outside the cell outline.  Movies at a rate of 1 frame per second were recorded for Dictyostelium cells moving on a solid support in buffer. The program Quimp3 represents the cell outline using a polygon of ,150 nodes, and then identifies extending convex protrusions as pseudopodia, described by the x,y,t coordinates of the start and end of their growing phase. The program also calculates the tangent to the surface at the position where the pseudopod started. These data were used to calculate for each pseudopod the size l, the angle a relative to a specific point in space, the angle b relative to the tangent, and the angle Q relative to the previous pseudopod. doi:10.1371/journal.pcbi.1000874.g001

Author Summary
Even in the absence of external information, many organisms do not move in purely random directions. Usually, the current direction is correlated with the direction of prior movement. This persistent random walk is the typical way that simple cells or complex organisms move. Cells with poor persistence exhibit Brownian motion with little displacement. In contrast, cells with strong persistence explore much larger areas. We have explored the principle of the persistent random walk by analyzing how Dictyostelium cells extend protrusions called pseudopodia. These cells can extend a new pseudopod in a random direction. However, usually cells use the current pseudopod for alternating right/left splittings, by which they move in a persistent zig-zag trajectory. A stochastic model was designed for the persistent random walk, which is based on the observed angular frequencies of pseudopod extensions. Critical elements for persistent movement are the ratio of de novo and splitting pseudopodia, and, unexpectedly, the shape of the cell. A relatively round cell moves with much more persistence than a cell with an irregular shape. These predictions of the model were confirmed by experiments that record the movement of mutant cells that are specifically defective in pseudopod splitting or have a very irregular shape.
surface area I; see Fig. S4). The intermediate ellipse was positioned in such a way that O~I (this also implies that the surface area of the cell (T) is identical to the surface area of the intermediate ellipse). The cell shape parameter is defined as Y~(OzI)=T; it holds that 0ƒYƒ2. For a cell with a regular shape that approaches a smooth ellipsoid, the surface areas O and I are very small and Y approaches zero. In contrast, O and I are larger for a cell with a very irregular shape; the largest value observed among ,600 cells was Y = 0.92.
With the exception of 5h starved cells, each database contains information from 200-300 pseudopodia, obtained from 6-10 cells, using two independent movies. For 5h starved cells, we collected a larger database containing 835 pseudopodia from 28 cells using 4 independent movies, and typical databases for each mutant. The data are presented as the means and standard deviation (SD) or standard error of the means (SEM), where n represents the number of pseudopodia or number of cells analyzed, as indicated in Table 1.

Statistical analysis of pseudopod angles
The probability density functions of angles can not be analyzed as the common distribution on a line. Angular distributions belong to the family of circular distributions, which are constructed by wrapping the usual distribution on the real line around a circle. The data were analyzed with two circular distributions, the von Mises distribution (vMD), which matches reasonably well with the wrapped normal distribution, and the wrapped Cauchy distribution (WCD), which has fatter tails [22]. The vMD is given by where I 0 (k) is the modified Bessel function of the first kind of order zero The WCD is given by , 0vrv1: ð3Þ

Monte Carlo simulations
Pseudopod extension is an ordered stochastic event [12]. The position of the tip of the formed pseudopodia depends on pseudopod size l p , splitting fraction s, Left/Right alternating ratio a, angle between split pseudopodia Q and variance of this angle s Q . A Monte Carlo simulation starts with a random angle a(1) of the first pseudopod. For the next and all subsequent pseudopodia the simulation uses four uniformly distributed random numbers R i,n (i = 1, .., 4) to calculate a(n), the angle of the n th pseudopod: R 1,n [ ½0,1 with the decision to split if R 1,n ,s; R 2,n [ ½0,1 with the decision for alternating splitting if R 2,n ,a; R 3,n [ ½0,1 for direction of split after de novo with decision right if R 3,n ,0.5; and R 4,n [ ½{180,180 for the direction of the de novo pseudopod. These probabilities result in a projected angle of extension in degrees. Finally, the actual pseudopod direction is drawn from a wrapped von Mises distribution with this projected angle as mean and s Q 2 as variance (k = 1/s Q

2
; variance converted to radians). The obtained a(n) and the pseudopod size l p are used to calculate the x,y coordinates of the tip of the pseudopod, followed by a next round of four random numbers to calculate a(n+1). In the simulations reported here we did not include stochastic variation in pseudopod length and pseudopod frequency, since we observed that they had only minor effects on the trajectories over several cell lengths.
Please note that in the simulations the direction of the simulated de novo pseudopodia is random; consequently, a small fraction of de novo pseudopodia are in the same direction of the previous pseudopod, which would be recognized in experiments as splitting pseudopodia. Conversely, a small fraction of the simulated splitting pseudopodia have angles The observed pseudopod data were derived from wild type (WT) cells starved for the indicated times, or from mutants starved for 5h. Indicated are the number of cells and pseudopodia analyzed. The data for l p are the means and SEM for pseudopodia; the data for Q and s Q are the means and SD for pseudopodia; the data for a and s are the means and SEM for cells. The data for the observed c obs were obtained by least square curve fitting to Eq. 4 of the dispersion of 20 cells during 15 minutes; the error indicates the 95% confidence level of the fit; h = cos 21 c. The data for the MC simulated c MC were obtained by simulating 100,000 trajectories using the observed pseudopod parameters, and the corresponding value of c was computed as described for the experimentally observed c. Finally, c step was calculated with Eq. 9, only using pseudopod parameters. doi:10.1371/journal.pcbi.1000874.t001 much larger than 55 degrees and would be recognized in experiments as de novo pseudopodia. From the geometry of the cell, we estimate that the number of simulated de novo in the current pseudopod and the number of splitting pseudopodia outside the current pseudopod are approximately the same, suggesting that the simulations represent the observed ratio of splitting and de novo pseudopodia.

Pseudopod extensions
The angles between pseudopodia were analyzed in detail and the results are presented in Fig. 2. For splitting pseudopodia, the angle between the current and next pseudopod (Q 1,2 ) has a clear bimodal distribution (Fig. 2A). A probability density function (PDF) of angles belongs to the family of circular or wrapped distributions. The data reported in this study were all fitted well by a von Mises distribution (vMD), which is the circular analog of the normal distribution. The wrapped Cauchy distribution has fatter tails and provided a poorer fit of the data (data not shown). The bimodal vMD presented in Fig. 2A is symmetric, yielding two means (Q 1,2 = +/255) that have the same variance k = 1/s Q 2 ; s Q1,2 = 28 degrees). Figure 2B shows the PDF of the angle between the current and next-next pseudopod (Q 1,3 ), which is best described by a single vMD with a mean of Q 1,3 = 2 degrees and s Q1,3 = 42 degrees. Figure 2C reveals that there is no significant correlation between the magnitude of angles between first/second pseudopod and the magnitude of the angles between second/ third pseudopod (thus e.g. splitting at a larger angle is not followed by a split at a smaller angle). The extension of splitting pseudopodia is summarized in Fig. 2D, and is based on the previous observation that a pseudopod split to the right is frequently followed by a split to the left and visa versa [12]. Thus the next pseudopod is extended at an angle of ,55 degrees to the right or left relative to the current pseudopod, and the next-next pseudopod is extended in roughly the same direction as the current pseudopod.
The angle between a de novo pseudopod and the previous pseudopod shows a very broad distribution (Fig. 2E). Nearly all angles between 2180 and +180 are well represented with a somewhat lower abundance of angles around 0 degrees. This suggests that a de novo pseudopod can be extended in any direction, but with slightly lower probability of the direction of the current pseudopod.

Trajectories of wild type and mutant cells
To investigate the consequence of the observed ordered extension of pseudopodia for cell movement on a coarse time scale for many pseudopodia we recorded the movement of Dictyostelium cells during 15 minutes; in this period about 30 pseudopodia are extended. Previously we have presented the cell trajectories for several strains and developmental stages [12] (see also Fig. S1). The mean square displacement as a function time, SD 2 (t)T, exhibits a slow approach to a linear function (Fig. 3A), which is typical for a transition of a correlated random walk at short times to a Brownian random walk after longer times [6,23]. Previously, the often used equation for a correlated random walk were fit to the data points to estimate persistence time and speed of the cells [12]. The aim of the present study is to analyze the mechanism of cell movement from the perspective of the extending pseudopodia, which have a specific length and direction. A correlated random walk in two dimensions can also be described with steps and turns [24,25]. With the replacement of the number of steps (n) in Eq. 7 in reference [25] for n = Ft we obtain where l is the step size in mm, F is the step frequency, and c is the correlation factor of dispersion (0,c,1), defined as the arithmetic mean of the cosine of the turn angle h between steps With three variables (F, l, c) the estimates of the parameters become uncertain. Fortunately, the step size can be deduced accurately from experimental data. As will be shown below in Eq. 10, the step size is given by l = l p cos(Q/2), where measurements for l p and Q are presented in Table 1. Using this value for l, the dispersion data were fitted to obtain the observed correlation factor of dispersion (c obs ) with the corresponding turn angle (h). In cells starved for 1 or 3 hours the correlation factor is only ,0.5 with turn angle of ,60 degrees. At 5 and 7 hours of starvation, cells move with much stronger persistence (correlation factor of 0.74 and 0.81 and a small turn angle of 42 and 36 degrees). Deletion of PLA2 or guanylyl cyclases prevents this increase of correlation factor, persistence is very low and cells disperse poorly.
How is pseudopod extension related to the observed correlation factor of dispersion c obs ? As previously stated (see Fig. 2), Dictyostelium cells may extend either de novo pseudopodia in nearly random directions, or splitting pseudopodia in a direction similar to the previous direction. Therefore, cells that extend exclusively de novo pseudopodia are expected to exhibit a random walk with c obs = 0 (turn angle h = 90 degrees), whereas cells extending exclusively splitting pseudopodia will exhibit strong persistence with large c and small turn angle h. As a consequence, c obs is expected to depend on the ratio s of splitting/de novo pseudopodia. Fig. 3B demonstrates that within experimental error this relationship is approximately linear; this holds true for the mutants as well as for wild type cells at different stages of development. The linear regression of all data yields c obs = 0.921s20.044. Thus, when all pseudopodia are de novo (s = 0) the correlation factor is small (c obs = 20.044) giving a turn angle h dn = 93 degrees, close to the expected value of 90 degrees for random turns. In contrast, when all pseudopodia are the result of splitting (s = 1) the correlation factor is large (c obs = 0.88) yielding a small turn angle (h s = 29 degrees). The implication of small turn angles for splitting pseudopodia will be discussed later.

Model for persistent movement
The alternating right/left extension of splitting pseudopodia can be used to simplify a description of the movement of Dictyostelium cells over longer distances. In this approach, the simplification may be valid for movement on a longer time scale only, as we study here, but may not be appropriate over shorter time scales of a few pseudopodia. Because pseudopodia are frequently extended alternating right/left, we consider movement by pairs of two pseudopodia. Figure 4 shows four possibilities of pairs of splitting pseudopodia, which are the RL, LR, RR and LL, each with corresponding probabilities and angles as indicated. In addition to these splitting pairs, three combinations with de novo pseudopodia are possible: split-de novo, de novo-split, and de novo-de novo. The correlation factor of dispersion yields for the seven pairs: De novo pseudopodia are extended in a random direction, i.e. Scos h ds T, Scos h sd T and Scos h dd T equal zero. The turn angles of the four splitting pairs are 0, Q and 2Q, as indicated in Fig. 4A, and the variance is approximately 2s Q 2 (see Fig. 2B). Consequently Eq. 6 reduces to:  where Scos (Q,k)T denotes the expected value of the cosines of the angles on a circle with weights given by the vMD with mean Q and variance given by k = 1/(2s Q

2
). Since all splitting pseudopodia show the same variance this can be further reduced to In this equation Scos (0,k)T is obtained by calculating the probabilities of all turn angles on a circle with the vMD using Eqs. 1 and 2 and then taking the weighted average of the cosines of these angles. Although this procedure is straightforward, Eq. 8 can be further simplified, because for s Q smaller than ,50 degrees a good approximation is Scos (0,k)T&½cos (2s Q = ffiffiffi p p ) 2 (see Fig. S2). Finally, on a longer time scale and averaged over many steps, the correlation factor of pairs is related to the correlation factor of its underlying two steps by c pair~c 2 step . With these replacements we obtain the analytical expression for the correlation factor Thus, the correlation factor c is the product of three terms: the splitting ratio s, a noise term with the variance s Q , and a term with right/left bias a and angle Q.
Finally, by considering movement in pairs of steps, Fig. 4 reveals that the step size of the displacement is given by Monte Carlo simulations of cell movement We used Monte Carlo simulations to investigate how l and c depend on the pseudopod parameters size l p , splitting fraction s, alternating ratio a, angle between split pseudopodia Q and  variance of this angle s Q 2 . These simulations are also useful to inspect whether step size l and correlation factor c are correctly described by Eqs. 8-10. The direction in which a pseudopod is extended appears to be an ordered stochastic event [12] that depends on multiple decisions according to the following scheme: The next pseudopod is a splitting or de novo according to the ratio s. A splitting pseudopod is extended with angle Q to the right or left relative to the previous splitting according to the alternating right/ left bias a. A de novo pseudopod is extended in a random direction. The splitting pseudopod that appears after a de novo pseudopod has an equal probability to be extended to the left or right. Finally, the direction of the emerging pseudopod has a variance s Q 2 . The Monte Carlo simulation starts with a random angle a(1) of the first pseudopod and then uses the probabilities for splitting fraction s, alternating ratio a, angle between split pseudopodia Q and variance of this angle s Q 2 to stochastically simulate the angle of the next pseudopod (see methods and Table 1 for pseudopod parameters of wild type and mutant cells). The simulated trajectories are qualitatively similar to the experimentally observed trajectories (see Fig. S1B): Fed wild type cells or mutants with abundant de novo pseudopodia make many turns and have small displacement, whilst the trajectories of starved wild type cells with abundant pseudopod splitting are more persistent with large displacements.
To investigate how the correlation factor c depends on pseudopod parameters, the displacement SD 2 (n)T was calculated from 100,000 trajectories obtained by MC simulation using a unit pseudopod size and different values of s, a, Q and s Q . The obtained displacement SD 2 (n)T was then fitted to Eq. 4 to obtain estimates for the step size l and Monte Carlo correlation factor, c MC . The symbols in Fig. 5 show the results of the MC simulation, whereas the curves are the result of Eqs. 8-10.
We first investigated the angle Q between splitting pseudopodia and the alternating right/left bias a. When all splitting pseudopodia are alternating (a = 1), the cells make a nearly perfect zig-zag trajectory, and therefore the angle Q has very little effect on the persistence factor c (Fig. 5A). When splitting pseudopodia are extended in a random fashion to the right or left (a = 0.5), the persistence factor c MC decreases sharply as Q becomes larger than ,30 degrees. At an intermediate right/left bias (a = 0.75) the persistence factor c MC remains relatively high as long as the angle between pseudopodia is below 60 degrees. The results of the MC simulation appear to be described very well by the simplified model (Eqs. [8][9][10]. Furthermore, at the observed angle of Q = 55 degrees and alternating factor of a = 0.77, the deduced persistence factor c MC is 0.88 (see asterisk in Fig. 5A).
The fraction of splitting pseudopodia has a major impact on the persistence factor c. In the MC simulations, the value of c declines approximately linearly with the value of s (Fig. 5B), as was also observed experimentally (Fig. 3B), and obtained in Eqs. 8 and 9. Finally, we investigated the contribution of the variance s Q 2 of the splitting angle to the persistence factor c. This reveals that the persistence decreases strongly with increasing variance (Fig. 5C), with c MC following an approximately linear relationship with cos(s Q ). The MC simulations are well described by Eq. 8, but deviate from Eq. 9 at s Q .30 degrees, as expected (see Fig. S2).
We also used these Monte Carlo simulations to obtain an estimate of the step size l. It appears that l does not to depend on s and a, but depends on Q according to l~l p cos (Q=2) (Inset Fig 5A), as was obtained in Eq. 10.
In summary, the obtained correlation factor from the MC simulation (c MC ) are nearly identical to the correlation factor calculated with Eq. 9 (c step ). This suggests that the movement of Dictyostelium cells is qualitatively and quantitatively described very well by the model of persistent steps and random turns with the observed pseudopod parameters l p , s, a, Q and s Q .

Movement of the centroid of the cell
How does the movement of pseudopodia relate to the movement of the centroid of the cell? The data presented in Table 1 reveal that the observed correlation factor c obs of the centroid for different cell types correspond well with the deduced correlation factors of the pseudopods (c MC and c step ), but is always larger by ,15% (Table 1). Apparently, the observed turn angle of the cell's centroid is smaller than the turn angle of the extending pseudopod. Inspection of movies of 5h starved AX3 cells confirm this notion: the average angle between splitting pseudopodia is 55628 degrees (Fig. 2A), while the centroid moves during period at an angle of only 31623 degrees (mean and SD). Equation 9 reveals that the correlation factor c step increases by 15% when Q = 55628 degrees for the pseudopod is replaced by Q = 31623 degrees for the cells centroid.
Probably two phenomena are responsible for the difference between pseudopod and centroid: extension of multiple pseudopodia and geometry of cells. When cells extend multiple pseudopodia it is likely that at any given instant of time, the front of the cell moves with a fixed fraction of the vector sum of velocities possessed by the pseudopodia active at that instant in time. The temporal overlap of two pseudopodia was deduced from the measured probability distributions of pseudopod extensions (Fig 2F in [26]), which reveal that ,25% of the pseudopodia overlap with another pseudopod during on average ,40% of their extension time. This suggest that the tip of the cell moves at an angle that is ,6 degrees smaller than 55 degrees. Secondly, geometry predicts that the rear of the cell makes smaller changes of direction than the tip of the cell, comparable to the differences in curvature made by the front and rear wheels of a car. Figure S3 indicates that for a stereotypic pseudopod at 55 degrees the directional change of the centroid is ,40 degrees (see Fig. S3). Together, multiple pseudopodia and cell geometry can explain observed difference between pseudopod and centroid changes of direction, leading to the small 15% difference between deduced pseudopod correlation factor (c MC and c step ) and observed centroid correlation factor (c obs ).

Directional displacement
The directional displacement SD D(0,n)T is the displacement after n steps in the direction of the first step. An expression for the directional displacement is especially relevant when the organism is exposed to positional cues leading to a drift in one direction, such as during chemotaxis. The directional displacement of a cell after extending one pseudopod at an angle h is l cos h, and for a population of cells SD D(0,1)T~lScos hT. By Eqs. 3 and 10, the displacement at the first step may be written as vD D(0,1)wc l p cos (Q=2), and at the i th step SD D(i{1,i)T~c i l p cos (Q=2), see Eq. 6 of reference [25]. The cumulative displacement after n steps is which at n?! is given by In essence, this equation describes the displacement of a cell population in which all cells extend the first pseudopod in the same direction. Subsequent pseudopodia are extended with a bias, which reduces geometrically with each step; the correlation factor c indicates how many pseudopodia have correlated direction and therefore how far the population will disperse in the direction of the first pseudopod. Figure 6 presents the directional displacement as observed experimentally in wild type cells. The displacement in the direction of the first pseudopod slowly decreases at each subsequent pseudopod, approaching random movement after ,10 pseudopodia. On average a cell moves ,15 mm in the direction of the first pseudopod, which is the equivalent of about 3 pseudopodia (given a pseudopod size of ,5 mm). This figure also presents the directional displacement as modeled by Eq. 11a with observed data for l p , Q and c, which is in very close agreement with experimental data, again suggesting that the movement of a cell is satisfactory described by the model with five pseudopod parameters.
Cell shape, the variance in the direction of pseudopod extension, and dispersal The variation in pseudopod direction s Q 2 plays an important role in Eqs. 8-11 describing cell dispersal. Previously [12] we have shown that the next pseudopod emerges at a specific distance d from the tip of the current pseudopod, and is then extended perpendicular to the cell surface (i.e. perpendicular to the tangent to the surface curvature at the position where the pseudopod emerges). The pseudopod direction is expected to have high confidence for cells with a smooth ellipsoid shape, because the local bending is very predictable. However, this confidence is much smaller for cells with a very irregular shape. We investigated the role of cell shape using three experiments. First we demonstrate that the variance s Q 2 indeed depends on the variance of the tangent and the normal to the tangent. Second, we show that wild type or mutant cells with irregular shape exhibit increased variance s Q 2 . Finally we show that, due to the increased variance, the mutant exhibits poor dispersal.
Quimp3 was used to construct the tangent to the surface curvature at the position where the pseudopod emerges. We first determined for wild-type cells the angle a t of this tangent relative to the previous pseudopod (a t = 34.5624.9 degrees), and the angle b of the new pseudopod relative to this tangent (b = 89.1613.3 degrees). As mentioned above, the observed angle of the new pseudopod relative to the previous pseudopod is Q = 55.2627.8 degrees. We expect that the angle of the tangent relative to the previous pseudopod is independent from the angle of the pseudopod relative to the tangent; therefore we expect s s 2 Q~s 2 t zs 2 b~2 4:9 2 z13:3 2~2 8:2 2 . Indeed, the observed standard deviation of 27.8 degrees is close to this expected value of Figure 6. Directional displacement. Pseudopod formation and trajectories were recorded for 5h starved Dictyostelium cells; 8 cells were followed during 15 min. The directional displacement is the distance moved after n pseudopodia in the direction of the first pseudopod. Data points are the means of ,120 measurements, the line is the outcome of Eq. 11 with pseudopod parameters as indicated in Table 1 28.2 degrees. Importantly, the largest contribution to s Q 2 is derived from the variance of the tangent s t 2 , which is related to the local shape of the cell.
In the collection of Dictyostelium mutants, we selected a strain with an irregular shape. Mutant ddia2-null with a deletion of the forH gene encoding the formin dDia2 has a star-like shape (Fig. 7C). In this mutant, new pseudopodia are extended at about the same frequency and distance from the present pseudopodia as in wild type cells, pseudopodia also grow perpendicular to the surface, and are extended roughly in the same direction of Q = 55 degrees as wild type cells (Fig. 7A). However pseudopodia exhibit much more variation in direction (s Q = 47 degrees compared to s Q = 28 degrees for wild type cells). Finally, we determined a shape parameter Y that indicates how much the cell outline deviates from an ellipse (see method section and Fig. S4). Figure 7B reveals that cells with increased irregular shape, either being wild-type or mutant, exhibit strongly increased variance s Q 2 . Importantly, the distance d and angle Q of the pseudopodia does not change with cell shape (Fig. 7A).
Using the observed values for s, a, Q, and s Q for ddia2-null cells we expect from Eq. 9 to obtain c step = 0.43, significantly lower compared to c step = 0.69. for wild type cells. Fig. 7D shows that the dispersion of ddia2-null cells is strongly reduced. The observed mean square dispersion was fitted to Eq. 4 yielding a correlation factor of c obs = 0.53 (Table 1), close to the value that was predicted from the extension of pseudopodia from an irregular surface.
In summary, these and previous results [12] suggest that a splitting pseudopod is induced at some distance d from the tip of the current pseudopod, and then grows perpendicular to the surface. In a cell with a regular shape, the tangent and therefore pseudopod direction can be approximated using the distance d; alternating R/L extensions lead to a relative straight zig-zag trajectory, providing strong persistence of movement. In a cell with a very irregular shape, the local curvature of the membrane at distance d is unpredictable. Consequently, alternating R/L splitting occur with large variation of directions, leading to frequent turns and poor persistence.

Discussion
The movement of many organisms in the absence of external cues is not purely random, but shows properties of a correlated random walk. The direction of future movement is correlated with the direction of prior movement. For organisms moving in two dimensions, such as most land-living organisms, this implies that movement to the right is balanced on a short term by movement to the left to assure a long-term persistence of the direction. In bipedal locomotion, the alternating steps with the left and right foot will yield a persistent trajectory. Amoeboid cells in the absence of external cues show ordered extension of pseudopodia: a new pseudopod emerges preferentially just after the previous pseudopod has stopped growth [12]. Importantly, the position at the cell Figure 7. Movement of a mutant with irregular shape. Pseudopod formation and trajectories were recorded for 5h starved wild type and mutant ddia2-null cells. For the next pseudopod we measured the angle Q and distance d relative to the current pseudopod, the angle a t of the tangent relative to the current pseudopod and the angle b of the pseudopod relative to the tangent. A. The data show the means and SD measured for 312 and 289 pseudopodia from wild type and ddia2-null cells, respectively. The means are not significantly different between wild type and ddia2null cells, but the SD s Q and s t are significantly larger for ddia2-null compared to wild type (Chi-square, P,0. surface where this new pseudopod emerges is highly biased. When the current pseudopod has been extended to the left (relative to the previous pseudopod), the next pseudopod emerges preferentially nearby the tip at the right side of the current pseudopod. Since pseudopodia are extended perpendicular to the cell surface, this next pseudopod is extended at a small angle relative to the current pseudopod [12]. Therefore, this (imperfect) alternating right/left pseudopod splitting resembles bipedal locomotion. Cells may also extend a de novo pseudopod somewhere at the cell body, which is extended in a random direction. In starved Dictyostelium cells, the probability of extending a de novo pseudopod is ,10-fold lower than of pseudopod splitting (probability calculated per mm circumference of the cell [12]).
The model for pseudopod-based cell dispersion depends on five parameters, the pseudopod size (l), the fraction of split pseudopodia (s), the alternating left/right bias (a), the angle between pseudopodia (Q) and the variance of this angle (s Q 2 ). With these parameters the experimental data on mean square displacement and directional displacement are well-explained using Eqs. 9 and 11, respectively. Pseudopodia are the fundamental instruments for amoeboid movement. The notion that the trajectories are described well by the five pseudopod parameters probably implies that we have identified the basic concept of the amoeboid correlated random walk: persistent alternating pseudopod splitting and formation of de novo pseudopodia in random directions.
The cells may modify one or more of these five pseudopod parameters in order to modulate the trajectories (see Table 1). Nearly all mutants, as well as wild type cells at different stages of starvation and development, have approximately the same average pseudopod size l p . In addition, the alternating right/left bias (a) fluctuates between 0.67 and 0.82, and the angle between splitting pseudopodia (Q) between 50 and 62 degrees. [pla2-null cells are the only exception [12]; emerging pseudopodia in pla2-null cells exhibit longer growth periods (,27 s) than wild type cells (,13 s), and are thus longer]. This suggests that all strains use the same mechanism for pseudopod splitting. In contrast to these constant properties of split pseudopodia, the fraction of split pseudopodia (s) changes dramatically upon starvation, and appears to be regulated by cGMP and PLA2 signaling. Well-fed cells extend pseudopodia that are predominantly de novo in random directions, leading to a nearly Brownian random walk [20]. Upon starvation, the appearance of cGMP and PLA2 signaling enhances splitting and suppresses de novo pseudopod extensions, which leads to more persistent movement. The important role of the fraction of splitting pseudopodia for cell movement is also depicted by the linear dependence of the correlation factor c on the fraction s of splitting pseudopodia (Fig. 3B and Eq. 9).
The variance of the angle of pseudopod extension (s Q 2 ) plays an important role in movement. In wild type cells, as well as in many mutant strains, s Q is about 28 degrees. The primary source of the variance of pseudopod angles lies in the variation of cell shape, by which the normal to the cell surface at a specific position on this surface will have significant variation. Since the direction of pseudopodia is given by this normal, it is predicted that a cell with irregular shape should have more variation in pseudopod direction, and consequently shows poor dispersion. The experiments with mutant ddia2-null cells strongly support this interpretation. Wild-type cells have a relatively regular spherical shape by which two nearby pseudopodia are extended in nearly the same direction (small s Q ). In contrast, mutant ddia2-null cells have an irregular star-like shape; therefore, two nearby pseudopodia are often extended in very different directions (large s Q ). The variance s Q 2 can be regarded as the noise of the system. It indicates how fast a cell that extends only alternating splitting pseudopodia (a = 1 and s = 1) will lose correlation of directionality. With s Q = 28 degrees for wild type cells it follows from Eq. 10 that after ten pseudopodia the correlation of direction is still ,0.5. In contrast, for ddia2-null cells we obtained s Q = 46.5 degrees, which implies that already after four pseudopodia the correlation of direction has declined to ,0.5. Supported by Monte Carlo simulations using the parameters of the mutant, we conclude that poor dispersion of ddia2-null cells is due to the increased variance of pseudopod angles, which is caused by its irregular shape. The correlation factor c is the product of three terms (see Fig. 5 and Eq. 9), namely: splitting fraction (s), alternating pseudopod angles (a and Q), and the SD of the pseudopod angle (s Q ). Strong persistence of cell movement is attained when all three terms are large and about equal in magnitude. Starved wild type cells follow this strategy: each term is ,0.9, resulting in the observed correlation factor of 0.74. Mutants in which one of these terms is compromised, such as reduced splitting in sgc/pla2-null cells or enhanced noise of ddia2-null cells, have poor dispersion.
In summary, the correlated random walk of amoeboid cells is well described by the balanced bipedal movement, mediated by the alternating right/left extension of splitting pseudopodia. Cells deviate from movement in a straight line due to noise and because cells occasionally hop or make random turns. The turns in particular are used by the cells to modulate the persistence time, thereby shifting between nearly Brownian motion during growth and strong persistent movement during starvation. Figure S1 Trajectories. Movies were recorded during 15 minutes and the trajectories of the centroid of ten cells were determined. A. Wild type Dictyostelium cells at different times after removing of food. The frequency and size of pseudopod extension is not very different, but starved cells extend predominantly splitting pseudopodia (see (16)). B. Trajectories of 5 hour starved mutant cells with deletions of genes encoding guanylyl cyclases (sGC and GCA) and PLA2. C. Monte Carlo simulations of the trajectories calculated with the pseudopod parameters that were obtained experimentally for the mutants as indicated in Table 1. In this equation, the notation Scos (0,k)T is the average of the cosines of the angles on a circle with weights given by the von Mises Probability Distribution (vMD) with mean of 0 degrees and variance given by k = 1/s Q 2 . The figure reveals that Scos (0,k)T deviates less than 2% from the simple expression ½cos ( The parameter of cell shape is defined as Y = (O+I)/T, where T is the surface area of the cell (grey + blue). The minimal value is Y = 0 when the cell is an ellipse, and the maximal value is Y = 2 for a cell with extreme long extensions; the maximal value observed among ,600 Dictyostelium cells is Y = 0.92. The neuron cell is shown for comparison. Found at: doi:10.1371.journal.pcbi.1000874.s004 (0.13 MB PDF) Author Contributions