×

Chemotaxis and random motility in unsteady chemoattractant fields: a computational study. (English) Zbl 1445.92050

Summary: We discuss a generic computational model which captures the effects of transient chemoattractant concentration on the chemotactic motility of individual cells. The model solves the appropriate unsteady chemoattractant transport equation using finite differences, while simultaneously executing biased random walks representing individual cells. The simulations were implemented for a 2D homogeneous domain, and two case studies were considered. In the first case study, we consider a single-point source at the origin of the domain which produces chemoattractant, while other cells execute biased random walks toward this point source. We observe that for continuous chemoattractant production, chemoattractant diffusivity has no effect on cell motility, as measured by the mean of time to reach the source. However, in the case of pulsed random production with a specific average duty cycle, the mean time-to-contact is generally minimal with respect to chemoattractant diffusivity over a moderate range of diffusivities. In the second case study, two mobile cells which simultaneously secrete chemoattractant are initially placed a certain distance apart and are then allowed to execute biased random walks. Our model shows that a pulsed random protocol for chemoattractant production facilitates the two cells “finding” one another compared to continuous production. From this case study we also learn that there exists a range of moderate chemoattractant diffusivities for which the mean time-to-contact is minimal when cells both produce/detect chemoattractant and chemotactically migrate. Using these case studies, we discuss how transience in chemoattractant concentration becomes important in characterizing the effectiveness of chemotaxis.

MSC:

92C17 Cell movement (chemotaxis, etc.)
Full Text: DOI

References:

[1] Alt, M., Biased random walk models for chemotaxis and related diffusion approximations, J. Math. Biol., 9, 147-177 (1980) · Zbl 0434.92001
[2] Anderson, A. R.A.; Chaplain, M. A.J., A mathematical model for capillary network formation in the absence of endothelial cell proliferation, Appl. Math. Lett., 11, 3, 109-114 (1998) · Zbl 0935.92024
[3] Asthagiri, A. R.; Lauffenburger, D. A., Bioengineering models of cell signaling, Annu. Rev. Biomed. Eng., 2, 31-53 (2000)
[4] Chaplain, M. A.J., Mathematical modeling of angiogenesis, J. Neuro-Oncol., 50, 37-51 (2000)
[5] Dallon, J. C.; Othmer, H. G., A continuum analysis of the chemotatic signal seen by Dictyostelium discoideum, J. Theor. Biol., 194, 461-483 (1998)
[6] Devreotes, P. N., Dictyostelium discoideuma model system for cell-cell interactions in development, Science, 245, 1054-1058 (1989)
[7] Dickinson, R. B.; Tranquillo, R. T., Transport equations and indices for random and biased cell migration based on single cell properties, SIAM J. Appl. Math., 55, 1419-1454 (1995) · Zbl 0843.60073
[8] Dunn, G. A.; Brown, A. F., A unified approach to analyzing cell motility, J. Cell Sci. Suppl., 8, 81-102 (1987)
[9] Francis, K.; Palsson, B. O., Effective intercellular communication distances are determined by the relative time constant for cyto/chemokine secretion and diffusion, Proc. Natl Acad. Sci. USA, 94, 12258-12262 (1997)
[10] Gordon, P., Nonsymmetric difference equations, J. Soc. Ind. Appl. Math., 13, 667-673 (1965) · Zbl 0144.18104
[11] Gourlay, A. R., Hopscotha fast second-order partial differential equation solver, J. Inst. Math. Appl., 6, 375-390 (1970) · Zbl 0218.65029
[12] Iijima, M.; Huang, Y. E.; Devreotes, P., Temporal and spatial regulation of chemotaxis, Dev. Cell, 3, 469-478 (2002)
[13] Keller, E. F.; Segel, L. A., Model for chemotaxis, J. Theor. Biol., 30, 225-235 (1971) · Zbl 1170.92307
[14] Lapidus, L.; Pinder, G. F., Numerical Solution of Partial Differential Equations in Science and Engineering (1982), Wiley: Wiley New York · Zbl 0584.65056
[15] Maheshwari, G.; Lauffenburger, D. A., Deconstructing (and reconstructing) cell migration, Micro. Res. Tech., 43, 358-368 (1998)
[16] Markus, M.; Böhm, D.; Schmick, M., Simulation of vessel morphogenesis using cellular automata, Math. Biosci., 156, 191-206 (1999) · Zbl 0933.92009
[17] Matsumoto, M.; Nishimura, T., Mersenne twisterA 623-dimensionally equidistributed uniform pseudorandom number generator, ACM Trans. Mod. Comp. Sim., 8, 1, 3-30 (1998) · Zbl 0917.65005
[18] Muthukkaruppan, V. R.; Kubai, L.; Auerbach, R., Tumor-induced neovascularization in the mouse eye, J. Natl Cancer Inst., 69, 699-708 (1982)
[19] Othmer, H. G.; Stevens, A., Aggregation, blowup, and collapsethe ABC’s of taxis in reinforced random walks, SIAM J. Appl. Math., 57, 1044-1081 (1997) · Zbl 0990.35128
[20] Palsson, E.; Othmer, H. G., A model for individual and collective cell movement in Dictyostelium discoideum, Proc. Natl Acad. Sci., 97, 10448-10453 (2000)
[21] Rice, R. G.; Do, D. D., Applied Mathematical and Modeling for Chemical Engineering (1995), Wiley: Wiley New York
[22] Rivero, M. A.; Tranquillo, R. T.; Buettner, H. M.; Lauffenburger, D. A., Transport models for chemotactic cell populations based on individual cell behavior, Chem. Eng. Sci., 44, 2881-2897 (1989)
[23] Scalerandi, M., Capogrosso Sansone, B., Condat, C.A., 2001. Diffusion with evolving sources and competing sinks: development of angiogenesis. Phys. Rev. E 65, Art. No. 011902.
[24] Serini, G.; Ambrosi, D.; Giraudo, E.; Gamba, A.; Preziosi, L.; Bussolino, F., Modeling the early stages of vascular network assembly, EMBO J., 22, 8, 1771-1779 (2003)
[25] Shenderov, A. D.; Sheets, M. P., Inversely correlated cycles in speed and turning in an amoebaan oscillatory model of cell locomotion, Biophys. J., 72, 2382-2389 (1997)
[26] Sleeman, B. D.; Wallace, I. P., Tumour induced angiogenesis as a reinforced random walkmodelling capillary network formation without endothelial cell proliferation, Math. Comput. Modelling, 36, 339-358 (2002) · Zbl 1021.92016
[27] Stokes, C. L.; Lauffenburger, D. A., Analysis of the roles of microvessel endothelial cell random motility and chemotaxis in angiogenesis, J. Theor. Biol., 3, 377-403 (1991)
[28] Tong, S., 2003. Private communication.
[29] Tong, S.; Yuan, F., Numerical simulation of angiogenesis in the cornea, Microvas. Res., 61, 14-27 (2001)
[30] Uber, J. G.; Murdock, L. C., Evaluation of the Hopscotch method for transient ground water flow, J. Hydraul Eng., 126, 615-626 (2000)
[31] Zigmond, S. H., Ability of polymorphonuclear leukocytes to orient in gradients of chemotactic factors, J. Cell Biol., 75, 606-616 (1977)
[32] Zygourakis, K., Quantification and regulation of cell migration, Tissue Eng., 2, 1, 1-26 (1996)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.