Migration of individual microvessel endothelial cells: stochastic model and parameter measurement

Analysis of cell motility effects in physiological processes can be facilitated by a mathematical model capable of simulating individual cell movement paths. A quantitative description of motility of individual cells would be useful, for example, in the study of the formation of new blood vessel networks in angiogenesis by microvessel endothelial cell (MEC) migration. In this paper we propose a stochastic mathematical model for the random motility and chemotaxis of single cells, and evaluate migration paths of MEC in terms of this model. In our model, cell velocity under random motility conditions is described as a persistent random walk using the Ornstein-Uhlenbeck (O-U) process. Two parameters quantify this process: the magnitude of random movement accelerations, alpha, and a decay rate constant for movement velocity, beta. Two other quantities often used in measurements of individual cell random motility properties--cell speed, S, and persistence time in velocity, Pv--can be defined in terms of the fundamental stochastic parameters alpha and beta by: S =square root (alpha/beta) and Pv = 1/beta. We account for chemotactic cell movement in chemoattractant gradients by adding a directional bias term to the O-U process. The magnitude of the directional bias is characterized by the chemotactic responsiveness, kappa. A critical advantage of the proposed model is that it can generate, using experimentally measured values of alpha, beta and kappa, computer simulations of theoretical individual cell paths for use in evaluating the role of cell migration in specific physiological processes. We have used the model to assess MEC migration in the presence of absence of the angiogenic stimulus acidic fibroblast growth factor (aFGF). Time-lapse video was used to observe and track the paths of cells moving in various media, and the mean square displacement was measured from these paths. To test the validity of the model, we compared the mean square displacement measurements of each cell with model predictions of that displacement. The comparison indicates that the O-U process provides a satisfactory description of the random migration at this level of comparison. Using nonlinear regression in these comparisons, we measured the magnitude of random accelerations, alpha, and the velocity decay rate constant, beta, for each cell path. We consequently obtained values for the derived quantities, speed and persistence time. In control medium, we find that alpha = 250 +/- 100 microns 2h-3 and beta = 0.22 +/- 0.03h-1, while in stimulus medium (control plus unpurified aFGF) alpha = 1900 +/- 720 microns 2h-3 and beta = 0.99 +/- 0.37h-1.(ABSTRACT TRUNCATED AT 400 WORDS)