# Nonmotile Single-Cell Migration as a Random Walk in Nonuniformity: The "Extreme Dumping Limit" for Cell-to-Cell Communications.

1. IntroductionTo understand many physiological processes in living organisms, such as embryogenesis and wound healing, among others, as well as their malfunctions, e.g., inflammatory diseases, tumor growth, and metastasis, it is fundamental and of great interest to comprehend the process of relocalization of cells, commonly known as cell migration. This term is used to refer to different processes that involve the movement of cells from one location to another. In living animals, embryonic development provides a clear example of importance of accuracy in cell migration, as errors in this process can result in birth defects. It is also known that proper cell migration is necessary for functional immune response and tissue repair in adults. Conversely, failure in cell migration or inappropriate migratory movements may result in life-threatening scenarios, such as autoimmune diseases, defective wound repair, inflammatory diseases, and tumor dissemination, promoting metastatic cancer progression [1-3].

The process of cell migration is very specific and depends on the cell type and the context of the migration process, thus several modes of cell migration have been described [4]. There are migrating cells that are self-propelled (self-driven, with autonomous mobility) and others that are nonmobile. In the case of bacteria, flagella-associated self-propulsion is an important virulence factor for some strains such as Escherichia coli and plays a crucial role in attachment to biomaterial surfaces and infection [5, 6]. Additionally, cells can move either as separate entities or by exhibiting a collective behavior. A mathematical model of cell collective movement has been reported [7] as well as a model that explains how swimming velocity of self-driven cells can increase in viscosity [8]. However, in this work we are interested in single-cell movement, which plays a crucial role in maintaining the homeostasis of the body (i.e., leukocyte migration through blood vessels), as well as in tumor progression and metastasis [9]. During these processes, a migrating cell travels through the body by a motion called random walk, yet this process does not fulfill the necessary understanding of the migration process of different cell types. For example, recent evidence has shown that leukocytes exhibit types of migratory behavior which differs from the previously described random walk [10] and that cells can undergo directed migration by the influence of chemical or mechanical signals originating from the environment [11-14] (Figure 1).

At first sight, the migration of nonmotile cells is essentially a random walk, very similar to thermally driven Brownian particles. This is also true for bacteria in suspension. The observation that, when suspended in water, small pollen grains are found to be in a very animated and irregular state of motion was first systematically investigated by Scottish botanist Robert Brown in 1827, and the observed phenomenon took the name of Brownian motion. Albert Einstein in 1905 [15] and sometime later Paul Langevin in 1908 [16] explained Brownian motion using different but equally successful mathematical approaches. Einstein's analysis was based on the diffusion equation:

[mathematical expression not reproducible], (1)

where D is diffusion coefficient, with initial condition f(x, t = 0) = [delta](x), where [delta](x) is Dirac's delta function. Solution of diffusion equation is given by [17]

[mathematical expression not reproducible], (2)

and the mean of the square of displacement is given by <[x.sup.2]> = 2Dt. On the other hand, Langevin started from Newton's equation of motion assuming a Stokes's drag force and a random thermal force due to continuous bombardment from molecules of the liquid. Although he did not exploit all richness of his model, Langevin obtained in the long-time regime Einstein's result, namely, by <[x.sup.2]> = 2([k.sub.B]T/6[pi]a[mu])t, where a is the radius of particle, [mu] is fluid viscosity, T is temperature, and [k.sub.B] is the Boltzmann's constant, hence bridging Brownian motion, random walk, and diffusion, a view soon quantified experimentally by Perrin [18]. Therefore, diffusion coefficient can be computed in terms of properties of the fluid in Brownian particles, which is the Einstein-Stokes formula D = ([k.sub.B]T)/ (6[pi]a[mu]).

In modern times, Langevin's approach is still used. In fact, during the last decades, physicist and mathematician modelers have viewed cell movement as a persistent random walk, which can be modeled using well-known stochastic differential equations. The most widely used model is the Ornstein-Uhlenbeck model (OU) [19], with certain predictions that until some years ago were in good agreement with observations and experimental results ([20]). However, recent discoveries seem to question and challenge the archetypical OU model. Specifically, Wua et al., reported that in three dimensions, the path of cells is more directional than random [21]. In addition, recent experiments show a scaling at the long-range regime [22], <[x.sup.2]> ~ [t.sup.[beta]], where the power [beta] is just a number in the interval 1 < [beta] < 2. It turns out that it is useful to introduce two functions of interest, namely, (i) logarithmic derivative of mean square of displacement (MS) [23], defined previously as

[beta](t) = [dMSD(t)/dt][t/MSD(t)], (3)

and (ii) kurtosis [24]:

[mathematical expression not reproducible], (4)

where [bar.x] is the mean value of the distribution function. In the special case, where the mean value is zero, kurtosis can be expressed simply, as

C(t) = <[x.sup.4]>/[<[x.sup.2]>.sup.2], (5)

since they allow us to make contact between experimental data and predictions of models. Both functions are dimensionless, and in addition, [beta](t) is defined such as in the power law cases, <[x.sup.2]> ~ [t.sup.[beta]], and it is a constant and precisely coincides with the power [beta]. In the framework of the OU model, kurtosis is a monotonic function and asymptotically reaches value 3 from below [19], while in [25], data show that kurtosis reaches value 2.3, while in [26] reaches asymptotic value from above. Therefore, in both cases, a departure from standard OU behavior is observed. In light of these findings, the persistent random motion of cell migration has been questioned, and other more complicated theoretical models have been introduced [24, 27] to confront with experimental results.

Therefore, the aim of the present work is to revisit the idea that single-cell motility can be described as a random walk. We point out that previous statements are only valid in case of free Brownian particles. However, if we introduce an external applied force and work in extreme dumping limit, the dynamics change completely, and predictions of the model depend each time on the form of assumed external potential. We define the concept of extreme dumping limit as the extreme nonuniform behavior observed for cell responses and cell-to-cell communications in vivo, within a given biological context. Contrary to previous studies, where authors usually solve the Newton-Langevin equation, here we work with corresponding Fokker-Planck equation and explain why it admits an exact solution for cases we have considered, and we show in plots the kurtosis as well as the logarithmic derivative of mean square of displacement versus time for three different simple models. We show that the nonstandard behavior seen experimentally can be reproduced in the framework of random motion with an applied external force within the extreme dumping limit scenario. Therefore, the random motion paradigm minimally extended can still be used to describe cell motility successfully without the need of more complicated and sophisticated models introduced previously in the literature.

2. Results and Discussion

2.1. Cell Movement as a Stochastic Process

2.1.1. Newton-Langevin and Fokker-Planck Equations. Let us consider a Brownian particle in one dimension in an external potential V (x) with a drag force from the medium proportional to the velocity-[lambda]v with A given by Stokes' law [lambda] = 6[pi]a[mu] and a random thermal force (stochastic variable) [xi](t) that considers the random collisions of the Brownian particle with the molecules of the liquid. This generalizes the OU model where there is no external applied force, and thus describes a free Brownian particle. Newton's equation of motion takes the form

ma = f(x) - [lambda]v + [xi](t), (6)

where m is the mass of the Brownian particle, f (x) = -V' (x) is the applied external force, and the random force [xi](t) is assumed to be a Gaussian white noise:

[mathematical expression not reproducible], (7)

where g = [lambda]kBT required by the so-called fluctuation-dissipation theorem [28]. In the case of a free particle, there is no external force, f(x) = 0, and we obtain the standard OU model [24]:

ma = -[lambda]v + [xi] (t), (8)

and assuming as initial condition x0 = x(0) = 0, we obtain for the mean squared of displacement (MSD) (in one dimension), the Furth formula [29]:

[mathematical expression not reproducible], (9)

where [tau] = m/[lambda] is the so-called persistent time. In the short-time regime and in the long-time regime, we obtain [30]

[mathematical expression not reproducible]. (10)

The system eventually exhibits diffusive behavior at late times, while at early times the dynamics are dominated by the inertia of the particle, and the behavior is ballistic, which has been observed in [31].

As it is known that cells are able to feel and sense certain environmental cues such as the stiffness of their environment [32] or surface nanoscale patterning [33], it would be more suitable to extend the standard OU model by considering an applied force that models anything from the external environment that perturbs the movement of the cell, such as signals, stimuli, etc. Furthermore, our focus on extreme dumping limit in which coefficient [lambda] in drag force from the fluid is so large that acceleration term m a can be neglected. In other words, we can write Newton's equation of motion in the following form:

ma + [lambda]v = f(x) + [xi](t), (11)

and assume that [lambda]v [much greater than] ma. In this case the Langevin-Newton equation takes the simpler form:

[mathematical expression not reproducible], (12)

where we have defined A(x) = f (x)/[lambda] and [eta](t) = [xi](t)/[lambda]. The latter is the Langevin equation for the process x( t) where the noise [eta](t) satisfies <[eta]([t.sub.1])[eta]([t.sub.2])> = 2D[delta]([t.sub.1] - [t.sub.2]) with D = [k.sub.B]T/[lambda]. Since the random force [eta](t) is not known, we can only compute mean values of powers of the position, or the moments, <[x.sup.n]>, once the density probability function is known. The density probability function u(t, x) satisfies the corresponding Fokker-Planck (FP) equation [34]:

[mathematical expression not reproducible], (13)

where the first term is the diffusion term while the second term is due to the external force with a constant diffusion coefficient D. If we ignore the external applied force, the FP equation reduces to standard diffusion equation. That explains why Einstein's approach and Langevin's approach were equally successful. Solving the FP equation, we then can compute the moments performing the integrals:

[mathematical expression not reproducible], (14)

and it is a function of time. It is known that this type of Fokker-Planck equation can sometimes be recast in the usual diffusion equation. This happens when the following condition is satisfied [35]:

2DA' (x) + A[(x).sup.2] = [c.sub.0] + [c.sub.1]x + [c.sub.2][x.sup.2], (15)

and the reason why this happens is that the diffusion equation and the FP equation at hand have the same number of symmetries.

2.1.2. Applications: Three Concrete Simple Models. Next, we shall consider three cases in which we can find exact analytical solution of the FP equation.

(i) Constant force or linear potential V (x) = -[v.sub.drift][lambda]x. A constant force could be for example the gravitational force. In this case the FP takes the form

[mathematical expression not reproducible], (16)

where we have put [DELTA] = 2D, and it is trivial to check that the condition above is satisfied. So, the FK equation can be recast in the diffusion equation [[omega].sub.[tau]] = [[omega].sub.yy], and the solution is given by u(t,x) = f(t, x)[omega]([tau](t, x), y(t, x)) where f(t, x), y(t, x) and f(t, x) are given by [35]

[mathematical expression not reproducible]. (17)

Therefore, the solution finally is given by [36]

[mathematical expression not reproducible], (18)

and the mean values <[x.sup.2]> and <[x.sup.4]> are given by

[mathematical expression not reproducible], (19)

where we have made use of the Gaussian integrals:

[mathematical expression not reproducible]. (20)

It is easy to check that when [v.sub.drift] = 0, we recover Einstein's results for the pure diffusion case. Furthermore, it is easy to check that one can obtain the same result by solving the Langevin equation for x(t) with the initial condition x(0) = 0, namely,

[mathematical expression not reproducible], (21)

and then by squaring this expression and using the properties of the Gaussian white noise <[eta](t)> = 0 and <[eta]([t.sub.1])[eta]([t.sub.2])> = 2D[delta]([t.sub.1] - [t.sub.2]), one obtains the previous expression for MSD.

Looking at the expression for MSD obtained, the first term is the contribution of the diffusion term, while the second term is the contribution from the applied external force.

Therefore, the system exhibits the diffusive behavior only in the beginning of the evolution, contrary to the OU model, and eventually the deterministic force takes over.

This can explain the observation that sometimes cell motion is more directed than random, such as in the case of cancer cell migration [37] or during zebrafish development [38].

Now the kurtosis C(t) as well as the logarithmic derivative of MSD [beta](t) can be computed. We see that for the linear potential case, the kurtosis C(t) [right arrow] 3 when t [right arrow] 0 and C(t) [right arrow] 1 when t [right arrow] [infinity], while the mean squared displacement <[x.sup.2]> ~ 2Dt at early times (or [beta](t) = 1) and <[x.sup.2]> ~ [t.sup.2] at late times (or [beta](t) [right arrow] 2).

(ii) Harmonic oscillator or parabolic potential V (x) = [bar.x][lambda][x.sup.2]/2 that corresponds to a force f(x) = -[theta][lambda]x. It is quite common to model attractive forces with springs ([39, 40]), and thus a harmonic trap is a reasonable potential to consider. In this case, the FP takes the form

[mathematical expression not reproducible], (22)

and the condition above is again satisfied. So, the FK equation can be recast in the diffusion equation [[omega].sub.[tau]] = [[omega].sub.yy], and the solution is given by u(t,x) = f(t, x)[omega]([tau](t, x), y(t, x)), where f(t, x), y(t, x) and [tau](t, x) are given by [36]

[mathematical expression not reproducible]. (23)

Therefore, the solution finally is given by [34]

[mathematical expression not reproducible], (24)

and using the same Gaussian integrals as before, the mean values <[x.sup.2]> and <[x.sup.4]> are given by

[mathematical expression not reproducible], (25)

and one can check that when [theta] = 0, we recover Einstein's results for pure diffusion, and we make use of the fact that

[mathematical expression not reproducible]. (26)

Therefore, the kurtosis in this case is always a constant in time C(t) = 3. On the other hand, at early times t [right arrow] 0 the mean squared <[x.sup.2]> ~ 2Dt (diffusive behavior or [beta](0) = 1), while at late times t [right arrow] [infinity]<[x.sup.2]> [right arrow] D/[theta](or [beta](t) [right arrow] 0).

(iii) Constant force and within a harmonic trap V (x) = [theta][lambda][x.sup.2]/2 - [v.sub.drift][lambda]x. This model combines the two previous cases, and by redefinition z = x - ([v.sub.drift]/[theta]), we recover the FK equation of the harmonic trap. Therefore, the solution reads

[mathematical expression not reproducible]. (27)

Finally, with help of Gaussian integrals and defining

[mathematical expression not reproducible], (28)

the moments are computed to be

[mathematical expression not reproducible]. (29)

The MSD in short-time regime is diffusive <[x.sup.2]> ~ 2Dt, while in long-time regime approaches a constant value <[x.sup.2]> [right arrow] (D/[theta]) + ([v.sup.2.sub.drift]/[[theta].sup.2]). Kurtosis at early times starts from the value 3, while eventually approaches the value

[C.sub.[infinity]] = (1 + 6b + 3[b.sup.2])/(1 + 2b + [b.sup.2]), (30)

where we have defined dimensionless quantity b = [theta]D/[v.sup.2.sub.drift]. This value depends on the interplay between [theta] and [v.sub.drift]. The drift velocity (linear term in potential) dominates b [much less than] 1 and [C.sub.[infinity]] [congruent to] 1, while [theta] (quadratic term in potential) dominates b [much greater than] 1 and [C.sub.[infinity]] [congruent to] 3. Therefore, it is possible to reproduce observation of [25] that kurtosis approaches 2.3 value. This can be achieved for b [congruent to] 0.33.

In both the second and third models, MSD asymptotically in time goes to a constant value due to the harmonic trap, and their logarithmic derivatives of MSD exhibit similar behavior, namely, they both are a monotonically decreasing function of time in the interval 0 < [beta](t) < 1.

In the corresponding Figures, we show functions C(t) and [beta](t) versus time for all three models. Figure 2 shows the kurtosis, while Figure 3 shows [beta](t) for the models considered here.

We see that each model exhibits its own dynamics, and they behave differently at late times, although in the short-range regime (t [right arrow] 0), they all exhibit diffusive behavior.

The simplest model with a constant applied force can explain (i) the scaling behavior with power [beta] = 2, (ii) the fact that cell movement can be more directed than random, and (iii) the monotonic decrease of kurtosis as seen in [26].

In addition, a more complicated model with a harmonic trap and constant force can explain asymptotic value of 2.3 seen previously [25]. Therefore, our results show that departures from OU model seen in recent experiments are also present here, and thus our findings suggest that the random motion paradigm minimally extended is still compatible with biological observations, at least qualitatively.

It would be interesting to obtain more data that could verify (or falsify) the predictions of the models considered here in a quantitative manner.

2.1.3. Biological Significance. Although it is widely accepted that cell migration is complex and multifactorial, our results show that this process can be described with a modification of the random walk model that also considers the application of external forces within a complex biological environment. Therefore, it remains possible to model cell migration in a minimalistic fashion without the need of using more complex calculations. Furthermore, this model is in line with several experimental observations in the laboratory. For example, despite the fact that fibroblasts can display random migration patterns in 2D tissue culture [20], many groups have observed that fibroblasts are also able to alter their migratory behavior when changed from a 2D to a 3D culture [4, 41]. This change in migration patterns is believed to be due to the ability of cells to detect and respond to stimuli in their environment, such as an increased number of adhesion points or mechanical differences of the matrix. Thus, given the 3D nature of the in vivo setting, it is possible to determine that cell migration within a living organism is more complex than a simple random walk, being also influenced by a number of external cues.

Furthermore, it has also been discussed that macrophages and neutrophils migrate using various modes of random walks (i.e., biased random walks) in response to acute injury, most of which include external factors that are guiding the cell towards an area of interest [42-44]. The presence of these factors (i.e., chemokines) are crucial in the biological setting, as they are responsible for generating effective inflammatory and wound healing responses [45]. Also, chemokine-based migration is imperative for osteoclast recruitment into bone tissue, which plays a big role in bone resorption observed in many chronic inflammatory diseases such as rheumatoid arthritis, periodontitis, and peri-implantitis [46, 47]. In all these scenarios, either constant or complex external forces are at play guiding cell migration towards the area of interest. However, our results strengthen the idea that cell migration is in fact a combination of stochastic processes and directed motion by external stimuli, which is supported by the calculations discussed in this research.

3. Conclusion

Studying properties of cell migration is of fundamental interest to understand many physiological processes in living organisms, as well as some pathological processes such as tumor metastasis or bacterial infection. Single-cell motility of nonmotile cells can be viewed as a random walk assuming that cells are thermally driven Brownian particles and can be modeled using well-known stochastic differential equations such as Langevin and Fokker-Planck equations. Recent experimental results have questioned the archetypical Ornstein-Uhlenbeck model, as they have shown some departures from standard predictions of persistent random motion paradigm based on a free Brownian particle driven by a drag Stokes's force, as well as by a random force due to molecule thermal motion. In the present work, we have revisited the issue of cell migration viewed as random walk by adding an applied external force and working in extreme dumping limit. We have studied three concrete cases for which the Fokker-Planck equation can be solved exactly, and we have provided analytical expressions for MSD and for certain quantities of interest that can be used to make contact observations. Our results show that predictions of the model and behavior of the system depend on form of the applied force (although all models exhibit diffusive behavior in short-time regime), and they all differ compared to the standard OU model. Our work shows that random motion paradigm minimally extended can still be used to describe cell motility successfully without introduction of sophisticated models. Overall, this model could be potentially beneficial to understand the migration behavior of cells during relevant biological processes such as wound healing, inflammation, and embryonic development from a minimalistic approach.

https://doi.org/10.1155/2018/9680713

Data Availability

Mathematical calculations utilized to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge grants provided to BioMAT'X Research Group, (Laboratorio de Biomateriales, Farmaceuticos y Bioingenieria de Tejidos CThneo Maxilo-Facial), member of CIIB (Centro de Investigation e Innovation Biomedica), Faculty of Dentistry, Universidad de los Andes, Santiago de Chile. The corresponding author acknowledges supplementary operating funding provided from CONICYT-FONDEF, Chile, under awarded project/grant (national) no. ID16I10366 (2016-2019) and Fondo de Ayuda a la Investigacion (FAI) Universidad de los Andes no. INV-IN-2015-101 (2015-2019).

References

[1] A. J. Ridley, M. A. Schwartz, K. Burridge et al., "Cell migration: integrating signals from front to back," Science, vol. 302, no. 5651, pp. 1704-1709, 2003.

[2] N. A. Wright, W. Hoffmann, W. R. Otto, M. C. Rio, and L. Thim, "Rolling in the clover: trefoil factor family (TFF)-domain peptides, cell migration and cancer," FEBS Letters, vol. 408, no. 2, pp. 121-123, 1997.

[3] A. J. Ridley, "Rho GTPases and cell migration," Journal of Cell Biology, vol. 150, no. 4, pp. 807-822, 2001.

[4] R. J. Petrie, A. D. Doyle, and K. M. Yamada, "Random versus directionally persistent cell migration," Nature Reviews Molecular Cell Biology, vol. 10, no. 8, pp. 538-549, 2009.

[5] B. Rodenborn, C.-H. Chen, H. L. Swinney, B. Liu, and H. P. Zhang, "Propulsion of microorganisms by a helical flagellum," Proceedings of the National Academy of Sciences, vol. 110, no. 5, pp. E338-E347, 2013.

[6] R. S. Friedlander, H. Vlamakis, P. Kim, M. Khan, R. Kolter, and J. Aizenberg, "Bacterial flagella explore microscale hummocks and hollows to increase adhesion," Proceedings of the National Academy of Sciences, vol. 110, no. 14, pp. 5624-5629, 2013.

[7] D. P. Stonko, L. Manning, M. Starz-Gaiano, and B. E. Peercy, "A mathematical model of collective cell migration in a three-dimensional, heterogeneous environment," PLoS One, vol. 10, no. 4, Article ID e0122799, 2015.

[8] Y. Magariyama and S. Kudo, "A mathematical explanation of an increase in bacterial swimming speed with viscosity in linear-polymer solutions," Biophysical Journal, vol. 83, no. 2, pp. 733-739, 2002.

[9] C. De Pascalis and S. Etienne-Manneville, "Single and collective cell migration: the mechanics of adhesions," Molecular Biology of the Cell, vol. 28, no. 14, pp. 1833-1846, 2017.

[10] H. B. Taylor, J. Liepe, C. Barthen et al., "P38 and JNK have opposing effects on persistence of in vivo leukocyte migration in zebrafish," Immunology and Cell Biology, vol. 91, no. 1, pp. 60-69, 2013.

[11] S. Wong, W.-H. Guo, and Y.-L. Wang, "Fibroblasts probe substrate rigidity with filopodia extensions before occupying an area," Proceedings of the National Academy of Sciences, vol. 11, no. 48, pp. 17176-17181, 2014.

[12] J. Ballester-Beltran, M. J. P. Biggs, M. J. Dalby, M. Salmeran-Sanchez, and A. Leal-Egana, "Sensing the difference: the influence of anisotropic cues on cell behavior," Frontiers in Materials, vol. 2, 2015.

[13] H. R. Bourne and O. Weiner, "A chemical compass," Nature, vol. 419, no. 6902, p. 21, 2002.

[14] T. Ishii, J. Kikuta, A. Kubo, and M. Ishii, "Control of osteoclast precursor migration: a novel point of control for osteoclastogenesis and bone homeostasis," IBMS BoneKEy, vol. 7, no. 8, pp. 279-286, 2010.

[15] a. Einstein, "Uber d4e von der rnolekuZarh%net4schem Theorie der WUrme geforderte Bewegumg vow in ruhenden P Z ikss4ykeitem," Annalen der Physik, vol. 322, no. 8, pp. 549-560, 1905.

[16] D. S. Lemons and A. Gythiel, "Paul Langevin ' s 1908 paper '" on the theory of brownian motion "' and on the theory of brownian motion, a note from M. P. Langevin, presented by M. Mascart," American Journal of Physics, vol. 65, no. 11, pp. 1079-1081, 1997.

[17] G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physicists, 2013.

[18] J. Baudrillard, Extrait de la publication, Wanda Sacher Masoch, Confession Ma Vie, 1972.

[19] G. E. Uhlenbeck and L. S. Ornstein, "On the theory of the Brownian motion," Physical Review, vol. 36, no. 5, pp. 823841, 1930.

[20] M. H. Gail and C. W. Boone, "The locomotion of mouse fibroblasts in tissue culture," Biophysical Journal, vol. 10, no. 10, pp. 980-993, 1970.

[21] P.-H. Wu, A. Giri, S. X. Sun, and D. Wirtz, "Three-dimensional cell migration does not follow a random walk," Proceedings of the National Academy of Sciences, vol. 111, no. 11, pp. 3949-3954, 2014.

[22] P. Bursac, G. Lenormand, B. Fabry et al., "Cytoskeletal remodelling and slow dynamics in the living cell," Nature Materials, vol. 4, no. 7, pp. 557-561, 2005.

[23] C. Metzner, C. Raupach, D. Paranhos Zitterbart, and B. Fabry, "Simple model of cytoskeletal fluctuations," Physical Review E, vol. 76, no. 2, 2007.

[24] D. Campos, V. Mendez, and I. Llopis, "Persistent random motion: uncovering cell migration dynamics," Journal of Theoretical Biology, vol. 267, no. 4, pp. 526-534, 2010.

[25] P. Dieterich, R. Klages, R. Preuss, and A. Schwab, "Anomalous dynamics of cell migration," Proceedings of the National Academy of Sciences, vol. 105, no. 2, pp. 459-463, 2008.

[26] H. Takagi, M. J. Sato, T. Yanagida, and M. Ueda, "Functional analysis of spontaneous cell movement under different physiological conditions," PLoS One, vol. 3, no. 7, Article ID e2648, 2008.

[27] D. Campos and V. Mendez, "Superdiffusive-like motion of colloidal nanorods," Journal of Chemical Physics, vol. 130, no. 13, article 134711, 2009.

[28] H. Y. Wang, S. Taylor, J. Simkin, J. M. Oakley, C. Emson, and M. J. Lancaster, "Simulation of microstrip small antennas," in Proceedings of 11th International Conference on Antennas and Propagation (ICAP 2001), pp. 17-20, Manchester, UK, April 2001.

[29] R. Furth, "Die Brownsche Bewegung bei Beriicksichtigung einer Persistenz der Bewegungsrichtung. Mit Anwendungen auf die Bewegung lebender Infusorien," Zeitschrift Fur Phys. vol. 2, no. 3, pp. 244-256, 1920.

[30] K. R. Naqvi, "The origin of the Langevin equation and the calculation of the mean squared displacement: let's set the record straight," Science and Technology, pp. 1-10, 2005.

[31] T. Li, S. Kheifets, D. Medellin, and M. G. Raizen, "Measurement of the instantaneous velocity of a brownian particle," Science, vol. 80, no. 328, pp. 1673-1675, 2010.

[32] M. Escude, M. K. Rigozzi, and E. M. Terentjev, "How cells feel: stochastic model for a molecular mechanosensor," Biophysical Journal, vol. 106, no. 1, pp. 124-133, 2014.

[33] J. H. Slater, P. J. Boyce, M. P. Jancaitis et al., "Modulation of endothelial cell migration via manipulation of adhesion site growth using nanopatterned surfaces," ACS Applied Materials and Interfaces, vol. 7, no. 7, pp. 4390-4400, 2015.

[34] H. Risken, The Fokker-Planck Equation, Springer, Berlin, Heidelberg, 1989.

[35] V. Stohny, "Symmetry properties and exact solutions of the fokker-planck equation," Journal of Nonlinear Mathematical Physics, vol. 4, no. 1-2, pp. 132-136, 1997.

[36] M. Brics, J. Kaupuzs, and R. Mahnke, "How to solve fokkerplanck equation treating mixed eigenvalue spectrum?," Condensed Matter Physics, vol. 16, no. 1, article 13002, 2013.

[37] E. T. Roussos, J. S. Condeelis, and A. Patsialou, "Chemotaxis in cancer," Nature Reviews Cancer, vol. 11, no. 8, pp. 573-587, 2011.

[38] J. Bussmann and E. Raz, "Chemokine-guided cell migration and motility in zebrafish development," EMBO Journal, vol. 34, no. 10, pp. 1309-1318, 2015.

[39] U. S. Schwarz, T. Erdmann, and I. B. Bischofs, "Focal adhesions as mechanosensors: the two-spring model," BioSystems, vol. 83, no. 2-3, pp. 225-232, 2006.

[40] O. K. Dudko, G. Hummer, and A. Szabo, "Intrinsic rates and activation free energies from single-molecule pulling experiments," Physical Review Letters, vol. 96, no. 10, 2006.

[41] B. Garber, "Quantitative studies on the dependence of cell morphology and motility upon the fine structure of the medium in tissue culture," Experimental Cell Research, vol. 5, no. 1, pp. 132-146, 1953.

[42] J. Liepe, H. Taylor, C. P. Barnes et al., "Calibrating spatiotemporal models of leukocyte dynamics against in vivo liveimaging data using approximate Bayesian computation," Integrative Biology, vol. 4, no. 3, p. 335, 2012.

[43] W. Alt, "Biased random walk models for chemotaxis and related diffusion approximations," Journal of Mathematical Biology, vol. 9, no. 2, pp. 147-177, 1980.

[44] P. J. M. Jones, A. Sim, H. B. Taylor et al., "Inference of random walk models to describe leukocyte migration," Physical Biology, vol. 12, no. 6, article 066001, 2015.

[45] L. B. Ford, V. Cerovic, S. W. F. Milling, G. J. Graham, C. A. H. Hansell, and R. J. B. Nibbs, "Characterization of conventional and atypical receptors for the chemokine CCL2 on mouse leukocytes," Journal of Immunology, vol. 193, no. 1, pp. 400-411, 2014.

[46] X. Yu, Y. Huang, P. Collin-Osdoby, and P. Osdoby, "CCR1 chemokines promote the chemotactic recruitment, RANKL development, and motility of osteoclasts and are induced by inflammatory cytokines in osteoblasts," Journal of Bone and Mineral Research, vol. 19, no. 12, pp. 2065-2077, 2004.

[47] S. E. Sahingur and W. A. Yeudall, "Chemokine function in periodontal disease and oral cavity cancer," Frontiers in Immunology, vol. 6, 2015.

Grigorios P. Panotopoulos, (1) Sebastian Aguayo (iD), (2) and Ziyad S. Haidar (iD) (1,3)

(1) BioMAT'X, Facultad de Odontologia, Universidad de Los Andes, Santiago, Chile

(2) School of Dentistry, Pontificia Universidad Catolica de Chile, Santiago, Chile

(3) CIIB, Facultad de Medicina, Universidad de Los Andes, Santiago, Chile

Correspondence should be addressed to Ziyad S. Haidar; zhaidar@uandes.cl

Received 10 July 2018; Revised 24 October 2018; Accepted 6 November 2018; Published 25 November 2018

Academic Editor: Panagiotis Kosmas

Caption: Figure 1: Cell migration is an important factor in physiological and pathological processes. (a) Inflammatory cells can migrate towards a site of interest via the sensing of chemokines and other inflammatory molecules. However, within the blood vessels, they are also subject to other relevant forces such as blood flow. (b) Osteoclast precursor cells are recruited into the tissues, where they can become activated by factors such as receptor activator of nuclear factor kappa-B ligand (RANKL) and cause bone resorption in health and disease. (c) Cell migration is also an important process that allows nonmotile bacteria to attach to surfaces, initiating biofilm formation. Migration of floating bacteria is also determined by important factors such as gravitational forces and flow. In all the above situations, migration is generated by a combination of stochastic (i.e., Brownian motion and random walk) and external forces (i.e., chemokines and flow).

Caption: Figure 2: Kurtosis versus time for three models considered in the present work. The constant function corresponds to the parabolic potential (second model, blue color), the one that eventually goes to 1 corresponds to the constant force (first model, red color), and the last one corresponds to the third model that combines the two (constant force in a harmonic trap, black color) for b = 0.5. Its asymptotic value depends on the interplay between [theta] and Vdrift and lies between 1 and 3.

Caption: Figure 3: Logarithmic derivative of MSD [beta](t) versus time for the three models considered in the present work. The decreasing function with values in the interval (0, 1] corresponds to the harmonic traps (second model in blue and third model in black), while the increasing function with values in the interval [1, 2) corresponds to the constant force case (first model red color).

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Panotopoulos, Grigorios P.; Aguayo, Sebastian; Haidar, Ziyad S. |

Publication: | Journal of Healthcare Engineering |

Date: | Jan 1, 2018 |

Words: | 5545 |

Previous Article: | Medical Mechatronics for Healthcare. |

Next Article: | Mobile Aid to Assist with Care Decisions in Children with Autism Spectrum Disorder (ASD). |

Topics: |