Mathematical model for characterization of bacterial migration through sand cores

The migration of chemotactic bacteria in liquid media has previously been characterized in terms of two fundamental transport coefficients-the random motility coefficient and the chemotactic sensitivity coefficient. For modeling migration in porous media, we have shown that these coefficients which appear in macroscopic balance equations can be replaced by effective values that reflect the impact of the porous media on the swimming behavior of individual bacteria. Explicit relationships between values of the coefficients in porous and liquid media were derived. This type of quantitative analysis of bacterial migration is necessary for predicting bacterial population distributions in subsurface environments for applications such as in situ bioremediation in which bacteria respond chemotactically to the pollutants that they degrade.We analyzed bacterial penetration times through sand columns from two different experimental studies reported in the literature within the context of our mathematical model to evaluate the effective transport coefficients. Our results indicated that the presence of the porous medium reduced the random motility of the bacterial population by a factor comparable to the theoretical prediction. We were unable to determine the effect of the porous medium on the chemotactic sensitivity coefficient because no chemotactic response was observed in the experimental studies. However, the mathematical model was instrumental in developing a plausible explanation for why no chemotactic response was observed. The chemical gradients may have been too shallow over most of the sand core to elicit a measurable response. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 53: 487-496, 1997.     

A simple expression for quantifying bacterial chemotaxis using capillary assay data: application to the analysis of enhanced chemotactic responses from growth-limited cultures.

An individual cell-based mathematical model of Rivero et al. provides a framework for determining values of the chemotactic sensitivity coefficient chi 0, an intrinsic cell population parameter that characterizes the chemotactic response of bacterial populations. This coefficient can theoretically relate the swimming behavior of individual cells to the resulting migration of a bacterial population. When this model is applied to the commonly used capillary assay, an approximate solution can be obtained for a particular range of chemotactic strengths yielding a very simple analytical expression for estimating the value of chi 0, [formula: see text] from measurements of cell accumulation in the capillary, N, when attractant uptake is negligible. A0 and A infinity are the dimensionless attractant concentrations initially present at the mouth of the capillary and far into the capillary, respectively, which are scaled by Kd, the effective dissociation constant for receptor-attractant binding. D is the attractant diffusivity, and mu is the cell random motility coefficient. NRM is the cell accumulation in the capillary in the absence of an attractant gradient, from which mu can be determined independently as mu = (pi/4t)(NRM/pi r2bc)2, with r the capillary tube radius and bc the bacterial density initially in the chamber. When attractant uptake is significant, a slightly more involved procedure requiring a simple numerical integration becomes necessary. As an example, we apply this approach to quantitatively characterize, in terms of the chemotactic sensitivity coefficient chi 0, data from Terracciano indicating enhanced chemotactic responses of Escherichia coli to galactose when cultured under growth-limiting galactose levels in a chemostat.     

Quantification of random motility and chemotaxis bacterial transport coefficients using individual-cell and population-scale assays

A number of individual-cell and population-scale assays have been introduced to quantify bacterial motility and chemotaxis. The transport coefficients reported in the literature, however, span several orders of magnitude, making it difficult to ascertain to what degree variations in bacterial species/strain, growth medium, growth and experimental conditions, and experiment type contribute to the reported differences in coefficient values. We quantified the random motility of Escherichia coli AW405 using the capillary assay, stopped-flow diffusion chamber (SFDC), and tracking microscope. We obtained good agreement for the random motility coefficient between these assays when using the same bacterial strain and consistent growth and experimental conditions. Chemotaxis of E. coli toward the attractant alpha-methylaspartate was quantified using the SFDC and capillary assay. Good agreement for the chemotactic sensitivity coefficient between the SFDC and the capillary assay was obtained across a limited attractant concentration range. Three different mathematical models were considered for analyzing capillary assay data to obtain a chemotactic sensitivity coefficient. These models differed by their treatment of the bacterial concentration in the chamber and the attractant concentration at the mouth. Results from our study indicate that the capillary assay, the most commonly used bacterial random motility and chemotaxis assay, can be used to accurately quantify bacterial transport coefficients over a limited range of attractant concentrations, provided experiments are performed carefully and appropriate mathematical models are used to interpret the experimental data.     

Measurement of bacterial random motility and chemotaxis coefficients: I. Stopped-flow diffusion chamber assay

Bacterial chemotaxis, the directed movement of a cell population in response to a chemical gradient, plays a critical role in the distribution and dynamic interaction of bacterial populations in nonmixed systems. Therefore, in order to make reliable predictions about the migratory behavior of bacteria within the environment, a quantitative characterization of the chemotactic response in terms of intrinsic cell properties is needed.The design of the stopped-flow diffusion chamber (SFDC) provides a well-characterized chemical gradient and reliable method for measuring bacterial migration behavior. During flow through the chamber, a step change in chemical concentration is imposed on a uniform suspension of bacteria. Once flow is stopped, diffusion causes a transient chemical gradient to develop, and bacteria respond by forming a band of high cell density which travels toward higher concentrations of the attractant. Changes in bacterial spatial distributions observed through light scattering are recorded on photomicrographs during a 10-min period. Computer-aided image analysis converts absorbance of the photographic negatives to a digital representation of bacterial density profiles. A mathematical model (part II) is used to quantitatively characterize these observations in terms of intrinsic cell parameters: a chemotactic sensitivity coefficient, mu(0), from the aggregate cell density accumulated in the band and a random motility coefficient, mu, from population dispersion in the absence of a chemical gradient.Using the SFDC assay and an individual-cell-based mathematical model, we successfully determined values for both of these population parameters for Escherichia coli K12 responding to fucose. The values obtained were mu = 1.1 +/- 0. 4 x 10(-5) cm(2)/s and chi(o) = 8 +/- 3 +/- 10(-5) cm(2)/s. We have demonstrated a method capable of determining these parameter values from the now validated mathematical model which will be useful for predicting bacterial migration in application systems.     

Quasi-elastic light scattering from migrating chemotactic bands of Escherichia coli. III. Studies of band formation propagation and motility in oxygen and serine substrates.

A series of light scattering experiments have been performed to study both macroscopic aspects of band formation and propagation and microscopic motility parameters of Escherichia coli in the combined substrate gradients of oxygen and serine. From the band formation experiment the conclusion is drawn that a minimum threshold gradient of the substrate is required for bacteria to form a band. From the band propagation experiment in the serine substrate the motility coefficient mu and chemotactic coefficient delta are determined. A separate quasi-elastic scattering experiment has been made with a propagating band to obtain three microscopic motility parameters: mean twiddle time tau 1, mean run time tau 2, and mean run speed V2. Finally, a scaling argument is made to connect the macroscopic parameters mu and delta with the microscopic parameters tau 1, tau 2, and V2, thus achieving a unified understanding of macroscopic and microscopic aspect of chemotaxis.     

Measurement of the chemotaxis coefficient for human neutrophils in the under-agarose migration assay

Clinical and scientific investigations of leukocyte chemotaxis will be greatly aided by an ability to measure quantitative parameters characterizing the intrinsic random motility, chemokinetic, and chemotactic properties of cell populations responding to a given attractant. Quantities typically used at present, such as leading front distances, migrating cell numbers, etc., are unsatisfactory in this regard because their values are affected by many aspects of the assay system unrelated to cell behavioral properties. In this paper we demonstrate the measurement of cell migration parameters that do, in fact, characterize the intrinsic cell chemosensory movement responses using cell density profiles obtained in the linear under-agarose assay. These parameters are the random motility coefficient, mu, and the chemotaxis coefficient, chi, which appear in a theoretical expression for cell population migration. We propose a priori the dependence of chi on attractant concentration, based on an independent experimental correlation of individual cell orientation bias in an attractant gradient with a spatial difference in receptor occupancy. Our under-agarose population migration results are consistent with this proposition, allowing chemotaxis to be reliably characterized by a chemotactic sensitivity constant, chi 0, to which chi is directly proportional. Further, chi 0 has fundamental significance; it represents the reciprocal of the difference in number of bound receptors across cell dimensions required for directional orientation bias. In particular, for the system of human peripheral blood polymorphonuclear neutrophil leukocytes responding to FNLLP, we find that the chemotaxis coefficient is a function of attractant concentration, a following the expression: chi = chi 0NT0 f(a) S(a) Kd/(Kd + a)2 where Kd is the FNLLP-receptor equilibrium dissociation constant and NT0 is the total number of cell surface receptors for FNLLP. f(a) is the fraction of surface receptors remaining after down-regulation, and S(a) is the cell movement speed, both known functions of FNLLP concentration. We find that chi 0NT0 = 0.2 cm; according to a theoretical argument outlined in the Appendix this means that these cells exhibit 75% orientation toward higher attractant concentration when the absolute spatial difference in bound receptors is 0.0025NT0 over 10 micron. (For example, if NT0 = 50,000 this would correspond to a spatial difference of 125 bound receptors over 10 micron.) This result can be compared with estimates obtained from visual studies of individual neutrophils.     

Analysis of chemotactic bacterial distributions in population migration assays using a mathematical model applicable to steep or shallow attractant gradients

The mathematical model developed by Riveroet al. (1989,Chem. Engng Sci. 44, 2881–2897) is applied to literature data measuring chemotactic bacterial population distributions in response to steep as well as shallow attractant gradients. This model is based on a fundamental picture of the sensing and response mechanisms of individual bacterial cells, and thus relates individual cell properties such as swimming speed and tumbling frequency to population parameters such as the random motility coefficient and the chemotactic sensitivity coefficient. Numerical solution of the model equations generates predicted bacterial density and attractant concentration profiles for any given experimental assay. We have previously validated the mathematical model from experimental work involving a step-change in the attractant gradient (Fordet al., 1991Biotechnol. Bioengng.37, 647–660; For and Lauffenburger, 1991,Biotechnol. Bioengng,37, 661–672). Within the context of this experimental assay, effects of attractant diffusion and consumption, random motility, and chemotactic sensitivity on the shape of the profiles are explored to enhance our understanding of this complex phenomenon. We have applied this model to various other types of gradients with successful intepretation of data reported by Dalquistet al. (1972,Nature New Biol. 236, 120–123) forSalmonella typhimurum validating the mathematical model and supportin the involvement of high and low affinity receptors for serine chemotaxis by these cells.     

Growth rate effects on fundamental transport properties of bacterial populations

In many natural environments, bacterial populations experience suboptimal growth due to the competition with other microorganisms for limited resources. The chemotactic response provides a mechanism by which bacterial populations can improve their situation by migrating toward more favorable growth conditions. For bacteria cultured under suboptimal growth conditions, evidence for an enhanced chemotactic response has been observed previously. In this article, for the first time, we have quantitatively characterized this behavior in terms of two macroscopic transport coefficients, the random motility and chemotactic sensitivity coefficients, measured in the stopped-flow diffusion chamber assay. Escherichia coli cultured over a range of growth rates in a chemostat exhibits a dramatic increase in the chemotactic sensitivity coefficient for D-fucose at low growth rates, while the random motility coefficient remains relatively constant by comparison. The change in the chemotactic sensitivity coefficient is accounted for by an independently measured increase in the number of galactose-binding proteins which mediate the chemotactic signal. This result is consistent with the relationship between macroscopic and microscopic parameters for chemotaxis, which was proposed in the mathematical model of Rivero and co-workers. (c) 1993 John Wiley & Sons, Inc.     

Quantitative analysis of experiments on bacterial chemotaxis to naphthalene

A mathematical model was developed to quantify chemotaxis to naphthalene by Pseudomonas putida G7 (PpG7) and its influence on naphthalene degradation. The model was first used to estimate the three transport parameters (coefficients for naphthalene diffusion, random motility, and chemotactic sensitivity) by fitting it to experimental data on naphthalene removal from a discrete source in an aqueous system. The best-fit value of naphthalene diffusivity was close to the value estimated from molecular properties with the Wilke-Chang equation. Simulations applied to a non-chemotactic mutant strain only fit the experimental data well if random motility was negligible, suggesting that motility may be lost rapidly in the absence of substrate or that gravity may influence net random motion in a vertically oriented experimental system. For the chemotactic wild-type strain, random motility and gravity were predicted to have a negligible impact on naphthalene removal relative to the impact of chemotaxis. Based on simulations using the best-fit value of the chemotactic sensitivity coefficient, initial cell concentrations for a non-chemotactic strain would have to be several orders of magnitude higher than for a chemotactic strain to achieve similar rates of naphthalene removal under the experimental conditions we evaluated. The model was also applied to an experimental system representing an adaptation of the conventional capillary assay to evaluate chemotaxis in porous media. Our analysis suggests that it may be possible to quantify chemotaxis in porous media systems by simply adjusting the model's transport parameters to account for tortuosity, as has been suggested by others.     

Spatio-temporal structure of migrating chemotactic band of Escherichia coli. I. Traveling band profile.

We developed a rapid-scanning, light-scattering densitometer by which extensive measurements of band migration speeds and band profiles of chemotactic bands of Escherichia coli in motility buffer both with and without serine have been made. The purpose is to test the applicability of the phenomenological model proposed by Keller and Segel (J. Theor. Biol. 1971. 30:235) and to determine the motility (mu) and chemotactic (delta) coefficients of the bacteria. We extend the previous analytical solution of the simplified Keller-Segel model by taking into account the substrate diffusion which turns out to be significant in the case of oxygen. We demonstrate that unique sets of values of mu and delta can be obtained for various samples at different stages of migration by comparing the numerical solution of the model equation and the experimental data. The rapid-scanning technique also reveals a hitherto unobserved time-dependent fine structure in the bacterial band. We give a qualitative argument to show that the fine structure is an example of the dissipative structure that arises from a nonlinear coupling between the bacterial density and the oxygen concentration gradient. Implications for a further study of the dissipative structure in testing the Keller-Segel model of chemotaxis are briefly discussed.     

Quantitative relationships between single-cell and cell-population model parameters for chemosensory migration responses of alveolar macrophages to C5a

Phenomenological parameters from a mathematical model of cell motility are used to quantitatively characterize chemosensory migration responses of rat alveolar macrophages migrating to C5a in the linear under-agarose assay, simultaneously at the levels of both single cells and cell populations. This model provides theoretical relationships between single-cell and cell-population motility parameters. Our experiments offer a critical test of these theoretical linking relationships, by comparison of results obtained at the cell population level to results obtained at the single-cell level. Random motility of a cell population is characterized by the random motility coefficient, mu (analogous to a particle diffusion coefficient), whereas single-cell random motility is described by cell speed, s, and persistence time, P (related to the period of time that a cell moves in one direction before changing direction). Population chemotaxis is quantified by the chemotactic sensitivity, chi 0, which provides a measure of the minimum attractant gradient necessary to elicit a specified chemotactic response. Single-cell chemotaxis is characterized by the chemotactic index, CI, which ranges from 0 for purely random motility to 1 for perfectly directed motility. Measurements of cell number versus migration distance were analyzed in conjunction with the phenomenological model to determine the population parameters while paths of individual cells in the same experiment were analyzed in order to determine the single-cell parameters. The parameter mu shows a biphasic dependence on C5a concentration with a maximum of 1.9 x 10(-8) cm2/sec at 10(-11) M C5a and relative minima of 0.86 x 10(-8) cm2/sec at 10(-7) M C5a and 1.1 x 10(-8) cm2/sec in the absence of Ca; s and P remain fairly constant with C5a concentration, with s ranging from 2.1 to 2.5 microns/min and P varying from 22 to 32 min. chi 0 is equal to 1.0 x 10(-6) cm/receptor for all C5a concentrations tested, corresponding to 60% correct orientation for a difference of 500 bound C5a receptors across a 20 microns cell length. The maximum CI measured was 0.2. Values for the population parameters mu and chi 0 were calculated from single-cell parameter values using the aforementioned theoretical linking relationships. The values of mu and chi 0 calculated from single-cell parameters agreed with values of mu and chi 0 determined independently from population migrations, over the full range of C5a concentrations, confirming the validity of the linking equations. Experimental confirmation of such relationships between single-cell and cell-population parameters has not previously been reported.     

Stopped-flow chamber and image analysis system for quantitative characterization of bacterial population migration: Motility and chemotaxis ofEscherichia coli K12 to fucose

The directed movement of a bacterial population in response to a chemical gradient is known as bacterial chemotaxis and plays a critical role in the distribution and dynamic interaction of bacterial populations. A quantitative characterization of the chemotactic response in terms of intrinsic cell properties is necessary for making reliable predictions about the migratory behavior of bacterial populations within the environment. The design of the stopped-flow diffusion chamber (SFDC) provides a well-characterized chemical gradient and reliable method for measuring bacterial migration behavior. During flow through the chamber a step change in the chemical concentration is imposed on a uniform suspension of bacteria. Once flow is stopped a transient chemical gradient forms due to diffusion; bacteria respond by forming a band of high cell density that travels toward higher concentrations of the attractant. Sequential observations of bacterial spatial distributions over a period of about ten minutes are recorded on photomicrographs. Computer-aided image analysis of the photographic negatives converts light-scattering information to a digital representation of the bacterial density profiles. A mathematical model is used to quantitatively characterize these observations in terms of intrinsic cell parameters: a chemotactic sensitivity coefficient, χ₀, from the aggregate cell density accumulated in the band and a random motility coefficient, μ₀, from population dispersion in the absence of a chemical gradient. Using the SFDC assay and an individual cell-based mathematical model we successfully determined values for both of these population parameters forEscherichia coli K12 responding to fucose. The values we obtained were μ₀=1.1 ± 0.4 x 10^-5 cm²/sec and χ₀=8 ± 3 x 10^-5 cm²/sec. These parameters will be useful for predicting population behavior in application systems such as biofilm development, population dynamics of genetically-engineered bacteria released into the environment, and in situ bioremediation technologies.     

Studies of bacterial chemotaxis in defined concentration gradients. A model for chemotaxis toward L-serine

The details of the chemotactic response of Salmonella typhimurium to gradients of L-serine have been examined in some detail. Two relatively macroscopic techniques have been employed to measure the bacterial response. These include measurements of the average velocity as the bacterial population moves toward attractants, and measurement of the upward-to-downward flux ratio, R, in the stable preformed attractant gradients. The dependence of the average velocity on gradient appears to be hyperbolic in nature, while the flux ratio depends linearly on the gradient. These data suggest a microscopic model for the dependence of bacterial behavior on the serine gradient. The model involves a linear dependence of the mean lifetime of a bacterial trajectory on the gradient for those bacteria moving toward higher attractant concentration. Those moving toward low concentrations of attractant do not change the mean duration of their trajectories, or the speed at which a given bacterium swims through the solution. This model generates the observed dependences of the average velocity and flux ratio on gradient. Interpretation of the experimental data suggests that a gradient which increases serine concentration by a factor of 2 in 10 mm is sufficient to double the average duration of a trajectory for a bacterium moving directly up the gradient. The concentration dependence of the chemotactic response to serine is more complicated. It suggests that more than one receptor of serine may be involved in determining chemotactic behavior to this attractant.     

Quantification of bacterial chemotaxis by measurement of model parameters using the capillary assay

The capillary assay for quantitative characterization of bacterial motility and chemotaxis is analyzed in terms of a mathematical model for cell population migration, in order to determine values for the cell random motility coefficient, mu and the cell chemotaxis coefficient, chi. The analysis involves both analytical perturbation methods and numerical finite-difference techniques. Transient cell density profiles within the capillary tube are determined as they depend upon mu and chi, providing a means for estimating mu and chi from the common protocol measurements of cell accumulation in the tube at specified observation times. The effects of extraneous factors such as assay geometry, stimulus diffusivity, bacterial density, and observation time are thus separated from the intrinsic cell-stimulus interaction and response. This allows independent population measurements of cell chemosensory movement properties to be extrapolated to situations involving growth and competition of populations, for purposes of better understanding microbial population dynamics in systems of biotechnological and microbial ecological importance.     

Measurement of leukocyte motility and chemotaxis parameters with a linear under-agarose migration assay

The interpretation of quantitative assays for leukocyte chemotactic migration is usually made in terms of measurements such as leading front distance, total migrating cells, and leukotactic index. These quantities allow comparison of cellular migration behavior under specified conditions. They are not useful; however, for comparisons between systems or for correlation with in vivo performance, because they depend upon specific physical aspects of the assay system, such as the geometry, chemoattractant concentration and diffusivity, and observation time. It would be more helpful to measure intrinsic properties of cell movement that could be used for comparison between systems, for correlation with in vivo studies, and to increase our understanding of the cell physiology. In this paper we demonstrate a means of quantitating leukocyte random motility, chemokinesis, and chemotaxis in terms of parameters that do characterize intrinsic cell properties. These parameters are the random motility coefficient and the chemotaxis coefficient, which appear in theoretical models of cell migration. We examine how well such a model describes the leukocyte density profile data observed in a modified under-agarose assay having a linear geometry. Furthermore, we obtain values for the random motility coefficient (and its dependence upon the concentration of the attractant peptide FNLLP) and for the chemotaxis coefficient for leukocytes responding to FNLLP.     

Residence time calculation for chemotactic bacteria within porous media

Local chemical gradients can have a significant impact on bacterial population distributions within subsurface environments by evoking chemotactic responses. These local gradients may be created by consumption of a slowly diffusing nutrient, generation of a local food source from cell lysis, or dissolution of nonaqueous phase liquids trapped within the interstices of a soil matrix. We used a random walk simulation algorithm to study the effect of a local microscopic gradient on the swimming behavior of bacteria in a porous medium. The model porous medium was constructed using molecular dynamics simulations applied to a fluid of equal-sized spheres. The chemoattractant gradient was approximated with spherical symmetry, and the parameters for the swimming behavior of soil bacterium Pseudomonas putida were based on literature values. Two different mechanisms for bacterial chemotaxis, one in which the bacteria responded to both positive and negative gradients, and the other in which they responded only to positive gradients, were compared. The results of the computer simulations showed that chemotaxis can increase migration through a porous medium in response to microscopic-scale gradients. The simulation results also suggested that a more significant role of chemotaxis may be to increase the residence time of the bacteria in the vicinity of an attractant source.     

Chemotaxis increases vertical migration and apparent transverse dispersion of bacteria in a bench-scale microcosm

The success of in situ bioremediation is often limited by the inability to bring bacteria in contact with the pollutant, which they will degrade. A bench-scale model aquifer was used to evaluate the impact of chemotaxis on the migration of bacteria toward the source of a chemical pollutant. The model was packed with sand and aqueous media was pumped across horizontally, simulating groundwater flow in a homogenous aquifer. A vertical gradient in chemoattractant was created by either a continuous injection of sodium benzoate or a pulse injection of sodium acetate. A pulse of chemotactic Pseudomonas putida F1 or a non-chemotactic mutant of the same species was injected below the attractant. The eluent was sampled at the microcosm outlet to generate vertical concentration profiles of the bacteria and chemoattractant. Moment analysis was used to determine the center and variance of the bacterial profiles. The center of the chemotactic bacterial population was located at an average of 0.74 ± 0.07 cm closer to the level at which the chemoattractant was injected than its non-chemotactic mutant in benzoate experiments (P < 0.015) and 0.4 ± 0.2 cm closer in acetate experiments (P < 0.05). The transverse dispersivity of the chemotactic bacteria was 4 ± 1 × 10(-3) cm higher in benzoate experiments than the transverse dispersivity of the non-chemotactic mutant and 1 ± 2 × 10(-3) cm higher in acetate experiments. These results underscore the contribution of chemotaxis to improve transport of bacteria to contaminant sources, potentially enhancing the effectiveness of in situ bioremediation.     

Platelet-derived growth factor in chemotactic for fibroblasts

Chemotaxis assays in modified Boyden chambers were used to detect fibroblast chemoattractants in materials released from early-stage inflammatory cells, namely, mast cells, platelets, and neutrophils. Strong attractant activity was found in substances released from platelets. This activity was accounted for mainly by the platelet- derived growth factor (PDGF), which is released from the platelets and which was active as a chemoattractant at 0.5-1.0 mitogenic units/ml. The mitogenic activity of purified PDGF, measured by [3H]thymidine incorporation, occurs at a similar concentration range. By varying the gradient of PDGF, we demonstrated that PDGF stimulates chemotaxis rather than random motility. Preincubation of suspensions of fibroblasts in the presence of PDGF decreased the subsequent migration of cells to a gradient of PDGF as well as to a gradient of fibronectin, which is also in attractant for fibroblasts. The chemotactic response of fibroblasts to PDGF was not inhibited by hydroxyurea or azidocytidine but was inhibited by actinomycin D and cycloheximide, suggesting that synthesis of RNA and proteins but not of DNA is required for the chemotactic response to occur. Fibroblast growth factor, epidermal growth factor, nerve growth factor, and insulin were not chemotactic for human skin fibroblasts, suggesting that the chemoattractant activity of PDGF for fibroblasts is not a general property of growth factors and mitogens. These results suggest that PDGF could have two functions in wound healing: to attract fibroblasts to migrate into the clot and then to induce their proliferation.