help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Biophys J. BioFAST on November 8, 2004.
doi:10.1529/biophysj.104.047001
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.104.047001v1
88/2/1432    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mody, N. A.
Right arrow Articles by King, M. R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Mody, N. A.
Right arrow Articles by King, M. R.
Biophysical Journal 88:1432-1443 (2005)
© 2005 The Biophysical Society

Mechanics of Transient Platelet Adhesion to von Willebrand Factor under Flow

Nipa A. Mody *, Oleg Lomakin *, Teresa A. Doggett {dagger}, Thomas G. Diacovo {dagger} and Michael R. King *

* Department of Biomedical Engineering, University of Rochester, New York; and {dagger} Departments of Pediatrics and Pathology, Washington University School of Medicine, St. Louis, Missouri

Correspondence: Address reprint requests to Michael R. King, Dept. of Biomedical Engineering, University of Rochester, Box 639, Medical Center, 601 Elmwood Ave., Rochester, NY 14642. Tel.: 585-275-3285; Fax: 585-273-4746; E-mail: mike_king{at}urmc.rochester.edu.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX 1
 APPENDIX 2
 ACKNOWLEDGEMENTS
 REFERENCES
 
A primary and critical step in platelet attachment to injured vascular endothelium is the formation of reversible tether bonds between the platelet glycoprotein receptor Ib{alpha} and the A1 domain of surface-bound von Willebrand factor (vWF). Due to the platelet's unique ellipsoidal shape, the force mechanics involved in its tether bond formation differs significantly from that of leukocytes and other spherical cells. We have investigated the mechanics of platelet tethering to surface-immobilized vWF-A1 under hydrodynamic shear flow. A computer algorithm was used to analyze digitized images recorded during flow-chamber experiments and track the microscale motions of platelets before, during, and after contact with the surface. An analytical two-dimensional model was developed to calculate the motion of a tethered platelet on a reactive surface in linear shear flow. Through comparison of the theoretical solution with experimental observations, we show that attachment of platelets occurs only in orientations that are predicted to result in compression along the length of the platelet and therefore on the bond being formed. These results suggest that hydrodynamic compressive forces may play an important role in initiating tether bond formation.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX 1
 APPENDIX 2
 ACKNOWLEDGEMENTS
 REFERENCES
 
Platelet adhesion dynamics is of central importance in hemostasis and thrombosis. During these biological processes, circulating platelets are recruited from the blood stream to sites of vascular injury, an event dependent on two key platelet glycoprotein (GP) receptors: GPIb{alpha} (subunit of the GPIb-IX-V complex) and GPIIb-IIIa ({alpha}IIbßIII integrin) (Lopez and Dong, 1997Go; Ruggeri, 1997Go; Du and Ginsberg, 1997Go). One key ligand that supports platelet adhesion by interacting with both of these receptors is von Willebrand factor (vWF), a multimeric plasma glycoprotein that is deposited on the exposed subendothelium region of the damaged vessel (Sakariassen et al., 1979Go; Turitto et al., 1980Go; Coller et al., 1983Go). In regard to GPIb{alpha}-mediated adhesion, it is the A1 region of this multidomain molecule that is critical for supporting the initial attachment (tethering) and subsequent translocation of platelets in flow (Savage et al., 1996Go; Cruz et al., 2000Go). This interaction also stimulates activation of the platelet GPIIb-IIIa receptor, which is necessary for firm irreversible adhesion and platelet aggregation (Savage et al., 1996Go; Parise, 1999Go; Kasirer-Friede et al., 2004Go).

Experiments have been conducted by a number of research groups to determine the effects of shear stress on platelet adhesion and the complementary roles of various platelet receptors (GPIb{alpha} and GPIIb-IIIa) and their ligands (vWF, fibrinogen, vitronectin, fibronectin) in promoting normal and pathological thrombus formation (Turitto et al., 1980Go; Weiss et al., 1989Go; Ikeda et al., 1991Go; Savage et al., 1992Go, 1996Go; Alevriadou et al., 1993Go; Kroll et al., 1996Go; Grunemeier et al., 2000Go). In shear flow, platelets were observed to transiently tether as well as translocate on the injured vessel surface via GPIb{alpha}-vWF tether bonds before firm adhesion, an interaction reminiscent of selectin-mediated rolling of leukocytes on inflamed endothelium (Springer, 1994Go; Savage et al., 1996Go; Cruz et al., 2000Go). Interestingly, Doggett et al. (2002Go, 2003Go) investigated the kinetics of the GPIb{alpha}-vWF-A1 tether bond and demonstrated that it exhibits similar biophysical attributes as the selectins; these include fast association and dissociation rates, force-dependent kinetics, and the requirement of a critical level of shear flow for adhesion to occur. This work has recently been confirmed by Kumar et al. (2003)Go. In contrast to leukocytes, however, tethered platelets do not roll, but instead exhibit a flipping motion in the direction of flow due to their flattened ellipsoidal shape. Thus, distinct mechanics imposed by the differences in the shapes of platelets versus leukocytes must govern tethering of these cell types in a flow environment even though they possess similar kinetic attributes.

In the current study, we sought to further explore and define the mechanics of platelet-surface contact by studying interactions between GPIb{alpha} and the vWF-A1 domain under varying hydrodynamic shear stress conditions. Unlike previous platelet experiments that have primarily involved bulk measurements of platelets interacting with a surface to determine the total number of adhering platelets (Savage et al., 1996Go; Cruz et al., 2000Go; Grunemeier et al., 2000Go), our work focused on examining the microscale motions of individual platelets during their transient interactions with vWF-A1. Images of platelet motion over vWF-A1-coated surfaces were captured during flow chamber experiments, for a range of fluid shear stresses (0.2–8.0 dyn/cm2). We used a platelet image processing algorithm developed in MATLAB to determine, from the digitized images, platelet three-dimensional orientation with respect to the surface as a function of time.

Since theoretical models have been successfully developed to aid in our understanding of primary and secondary leukocyte recruitment at sites of inflammation (Hammer and Apte, 1992Go; King and Hammer, 2001aGo,bGo), we anticipated that a suitable model for platelet mechanics would prove useful in providing physical insight into the platelet recruitment process. To this end, we developed an analytical model from first principles to characterize the flipping motion of a tethered platelet on an adhesive surface under linear shear flow and then compared predicted platelet trajectories with those obtained from flow chamber experiments. The hydrodynamic forces acting on the tethered platelet for various orientations during its flipping motion on the surface were also obtained from the model. Finally, the instantaneous orientation of platelets with respect to the surface at the time of surface capture and release were recorded and compared to the theory. We found a consistent correlation between the nature of the hydrodynamic force acting on the platelet and the angular orientation of the platelet at the time of attachment to or detachment from the surface. The results obtained may have important implications on the role of compressive hydrodynamic forces in enabling platelet adhesion.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS
 DISCUSSION
 APPENDIX 1
 APPENDIX 2
 ACKNOWLEDGEMENTS
 REFERENCES
 
Experiments
Platelet tethering experiments were carried out in a parallel-plate flow chamber apparatus at a range of fluid shear stresses as previously described (Cruz et al., 2000Go; Doggett et al., 2002Go). Briefly, platelets were purified from citrated whole blood (5 x 107 per ml) and perfused over plasma vWF (20 µg/ml) at fluid shear stresses ranging from 0.2 to 8 dyn/cm2. Images of platelets flowing over, attaching to, flipping on, and detaching from vWF-A1-coated surfaces were captured at 120 fps (235 fps in one set) using a Nikon microscope (Nikon, Tokyo, Japan) equipped with a 20x objective and high speed video camera (Speed Vision Technologies, San Diego, CA), and converted to digital movies. The total magnification corresponded to 7.57 pixels per micron.

Cell tracking methodology
For the purposes of cell tracking, the platelet was modeled as a rigid ellipsoid. The MATLAB image processing program first analyzed the entire sequence of digital grayscale images captured from phase contrast microscopic images of platelets to determine the dimensions of each cell. Each digital movie studied contained most of the range of possible platelet orientations; thus the major and minor platelet axes could be readily identified. To calculate the three-dimensional platelet orientation for a particular frame, an initial orientation was assumed based on the most recent platelet position, and its elliptical two-dimensional projection on the x,y plane was compared with the actual platelet projection. The simplex method for minimization of the sum of errors (distance between complementary points on the experimental platelet projection and the predicted projection) was used to determine, in proximity of the initial assumed position, the true platelet orientation. Note that the platelet can be tracked with high fidelity during the entire interaction using this technique. Fig. 1 illustrates graphically the algorithm's tracking results of one platelet flipping event. We have made this program available to others for use in tracking the motion of ellipsoidal particles (contact M.R.K.).



View larger version (51K):
[in this window]
[in a new window]
  		FIGURE 1  Graphical results of the platelet tracking algorithm developed to analyze three-dimensional position of tethering platelets from digital movies. (A) Digital movie frames of a tethered platelet flipping on a vWF-coated surface showing platelet orientation at five different times (i–v) during the flip. (B) Corresponding diagrams of ellipsoidal mesh models rotated approximately three axes and its calculated two-dimensional projected shadow. (C) Regression of calculated shadow outline (represented by thick solid line) to experimental image outline (represented by {circ}).

 
Analytical platelet flipping model
The platelet was modeled as a thin plate, hinged at one end and extending infinitely in the direction perpendicular to the flow and velocity gradient directions. We assumed in our model that the platelet-surface bond(s) provides no resistance to rotation and behaves as a frictionless hinge. The flow near the surface was approximated as linear shear flow. Because of the microscopic dimensions of this fluid system (Reynolds number Re = {rho}L2{gamma}/µ << 1, where {rho} is the fluid density, L is the characteristic length, {gamma} is the shear rate, and µ is the fluid viscosity) the flow near the wall is within the Stokes regime and inertial effects can be neglected. Taking advantage of the linearity of the Stokes flow equations, where v is the fluid velocity and p is the fluid pressure, the hydrodynamic problem of platelet flipping can be decomposed into the sum of two simpler problems (Fig. 2):
  1. Fence problem. Linear shear flow around a stationary fence inclined at angle {alpha} to an infinite rigid wall. The governing boundary conditions for this flow are v = 0 at {theta} = 0, {pi} and v = 0 at {theta} = {alpha} for 0 ≤ r ≤ 1, and v -> {gamma} zî as r -> {infty}, where {alpha} is the angle made by the fence with the plane and is measured in counterclockwise direction, {gamma} is the shear rate, and î is a unit vector in the x direction.
  2. Hinge problem. Flow due to rotation at angular velocity {omega} of a hinged plate oriented at angle {alpha} to a planar surface. The surrounding fluid is motionless at infinity. Edge effects due to the finite extent of the platelet were neglected. The analysis was done in two parts with flow fields for fluid regions on either side of the plate determined separately. The boundary conditions consist of the no-slip and no-penetration conditions at the solid surfaces of {theta} = 0, {pi}, {alpha}.



View larger version (19K):
[in this window]
[in a new window]
  		FIGURE 2  Schematic of how the Stokes flow problem of a platelet flipping on a surface may be modeled as a two-dimensional analytical problem by decomposing it into the sum of two simpler cases. Flipping is a result of the hydrodynamic effect of shear flow on the platelet as it interacts with surface-bound vWF. (A) Platelet attached to the surface by GPIb{alpha}-vWF bond(s) and subjected to linear shear flow. (B) A stationary and rigid inclined fence subjected to linear shear flow. (C) A rigid, flat, hinged plate rotating at angular velocity {omega} toward the surface in a quiescent fluid.

 
A solution to the fence problem has been previously determined by Jeong and Kim (1983)Go. They derived analytical expressions for force and torque exerted on the fence by the fluid, for all orientations of the fence with respect to the wall (Appendix 1). For the hinge problem, the necessary expressions for velocity profiles, force, and torque were derived, and are discussed in Appendix 2.


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
DISCUSSION
 APPENDIX 1
 APPENDIX 2
 ACKNOWLEDGEMENTS
 REFERENCES
 
Flipping of tethered platelets
A platelet tethering event was defined as the occurrence of a platelet binding to the surface such that the anchorage point did not move while the platelet flipped about this point. After completing one flip in the direction of flow, platelets were observed to either dissociate after some finite time, or undergo a second flip at a new anchorage point. The platelet flipping events that were analyzed lasted between 0.35 and 4 s. This is significantly longer than the largest value of characteristic time for force loading of the bond, which can be calculated as 1/(minimum shear rate) = 1/100 = 0.01 s. This is the average time taken for a (spherical) cell to move a distance equal to its radius in the direction of flow resulting in maximum stressing of the bond. In the experimental characterization, platelets that did not exhibit significant rotation about the x axis were selected for study, i.e., complex, fully three-dimensional trajectories were not analyzed. Seven out of 39 observed flipping events exhibited three-dimensional motion. These so-called three-dimensional trajectories involved either diagonal flipping with a velocity component normal to the flow direction or significant bending of the platelet. Platelets rolling on their edges were not observed. Nine representative flipping trajectories, spanning wall shear stresses ranging from 1.0 to 8.0 dyn/cm2, were chosen for comparison with the theoretical prediction. Fig. 3 compares the experimentally observed platelet flipping trajectories with that predicted by the two-dimensional analytical model. Importantly, the motion of the platelets was in good agreement with the predicted trajectories over a wide range of shear stresses (1–8 dyn/cm2). It is evident from Fig. 3 that when time is non-dimensionalized with respect to wall shear rate, the experimental platelet trajectories obtained for a range of shear stresses all collapse onto a single curve (within experimental error) as expected from the model.



View larger version (24K):
[in this window]
[in a new window]
  		FIGURE 3  Plot of analytically predicted and experimentally observed trajectories of bound platelets during platelet flips as a function of the dimensionless time. At time t = 0, the platelet is oriented at 90° to the surface. (Thick solid line represents the theoretical prediction; the symbols {diamond}, {bigtriangleup}, {square}, and {circ} represent four experimental flipping events at 8 dyn/cm2; the symbols {triangleright}, {diamond}, {bigtriangleup}, and {square} represent four experimental flipping events at 1.5 dyn/cm2; x represents a flipping event at 1 dyn/cm2.)

 
Radial bond force prediction
Our theoretical analysis also provides information on the magnitude and direction of force acting on the platelet, and thus yields an estimate of the hydrodynamic radial force experienced by the tethered platelet as a function of its orientation with respect to the surface. Making the assumption that the receptor-ligand bond(s) is oriented in the direction of the major axis of the platelet, one can thus calculate the radial force on the bond (see Discussion for motivation for this assumption). Fig. 4 B shows the magnitude and direction of non-dimensionalized radial force acting on the platelet-surface bond during the period of the flipping event. The theory predicts that the radial force on the bond changes sign during the platelet flip. The force on the bond is calculated to be compressive over the first half of the flip and tensile over the second half. Tees et al. (1993)Go demonstrated in their investigations of dissociation of sphere doublets under shear stress, that the doublet undergoes alternating tensile and compressive forces during its periodic rotation in the fluid. During one quarter of the complete orbit of the doublet (prolate spheroid) long axis, the doublet experiences compressive force which promote surface contact and bond formation (Neelamegham et al., 1997Go), and during the next quarter of the orbit, tensile forces dominate promoting bond dissociation. To estimate the magnitudes of compressive and tensile radial force on a platelet from the model, we assume a platelet size of 3 µm x 3 µm square and a fluid viscosity of 1.0 cP. The maximum absolute value of radial force on the bond is calculated as 0.93 pN and 10.2 pN at a shear rate of 73 s–1 and 800 s–1, respectively. Note that Doggett et al. (2002)Go showed a minimum shear threshold for tether bond formation at 73 s–1.



View larger version (22K):
[in this window]
[in a new window]
  		FIGURE 4  Plot of non-dimensionalized theoretical hydrodynamic radial bond force as a function of decreasing platelet orientation angle, compared with experimentally observed platelet attachment/detachment angles. (A) Schematic diagrams of three positions of a platelet during its flip, where {alpha} is the angle the platelet major axis makes with the surface, measured in counterclockwise direction. (B) Non-dimensionalized theoretical radial bond force, F/(2{gamma}µA), plotted as a function of platelet orientation angle where the angles vary from {pi} to 0 as the platelet flips from one side to the other; the experimentally observed attachment angles ({diamond}) and detachment angles ({circ}) obtained from experiments conducted at shear stresses of 1.0, 1.5, 2.0, 4.0, and 8.0 dyn/cm2 are plotted with the theoretical bond force estimate. Note that the occurrence of attachment/detachment events corresponds to the negative/positive sign of the predicted radial force, respectively.

 
Analysis of platelet attachment/detachment events
We analyzed experimental platelet tethering events using our platelet tracking/image processing algorithm to determine angles of platelet attachment to, and detachment from, the surface, and to correlate these angles with the predicted radial force acting on the ellipsoidal cell. Twenty-two attachment events and 30 detachment events were recorded for a total of 52. These events include all monitored platelet-surface contacts that resulted in attachment or detachment. We were unable to experimentally determine the number of contact events that result in attachment since we could not identify a non-adhesive contact event by visual observation; however, we could precisely determine an attachment event. The angles at which the platelets attached to the surface were all found to be >{pi}/2 whereas the detachment angles were found to be <{pi}/2, i.e., attachment occurred during the first-half of a platelet rotation ({alpha} = {pi} to {pi}/2) whereas detachment occurred during the second-half of a platelet flip ({alpha} = {pi}/2 to 0). Fig. 4 B compares the experimentally observed attachment and detachment angles with the theoretical estimate of the radial force acting on the platelet. From these data, we noted a consistent record of platelet attachment events occurring when the platelet is predicted to experience hydrodynamic compressive force, with no observations of attachment when the radial force is predicted to be tensile in nature.

We examined how the average platelet attachment and detachment angles depend on wall shear stress for the 52 attachment and detachment events analyzed. Fig. 5 shows that, within experimental error, the attachment angle is independent of shear stress over a range of 1.5–4 dyn/cm2. Interestingly, the detachment angle was found to be a decreasing function of shear stress. Pooling the low shear data (1.5–2 dyn/cm2) and similarly pooling the high shear data (4–8 dyn/cm2) and comparing the mean of the two samples, we find that this difference reaches statistical significance (p < 0.001, student's t-test). Both the independence of the attachment angle and the dependence of the detachment angle with wall shear stress agree with our intuitive expectations. Initial platelet attachment should depend primarily on the incidence of cell contact, and in Stokes flows the trajectories of suspended particles are independent of wall shear rate (or equivalently, shear stress divided by fluid viscosity; see Appendices). Conversely, during the cell detachment process both the rotation rate and the magnitude of the force exerted on the formed bond(s) each depend linearly on the wall shear stress. Since the GPIb{alpha}-vWf bond has been found to exhibit selectin-like bond kinetics, dissociation rate will depend strongly on force in an exponential manner (Doggett et al., 2002Go). Since the full force-loading history on the receptor-ligand bond during platelet rotation is highly nonlinear, it will require detailed adhesive dynamics simulations to fully analyze and predict how detachment angles should correlate with wall shear stress and other physical determinants (King and Hammer, 2001aGo,bGo).



View larger version (15K):
[in this window]
[in a new window]
  		FIGURE 5  Plot of the experimentally measured average attachment ({diamond}) and detachment angles ({circ}), as a function of wall shear stress.

 
Characterization of the trajectories of freely flowing platelets
We characterized the trajectories of freely flowing platelets as they rotated in the direction of flow near the vWF-A1 functionalized surface. According to Frojmovic et al. (1990)Go, the platelet can be viewed as an oblate spheroid with aspect ratio 0.25. Jeffery orbits (Jeffery, 1922Go) describe the periodic trajectories of oblate spheroids rotating in unbounded linear shear flow and were used to predict the rotational motion of free-stream platelets subjected to linear shear flow. The non-dimensionalized Jeffery equation yielding the periodic free-stream platelet trajectories is given by where {alpha} is defined as shown in Fig. 4 A, r is the aspect ratio of the platelet taken as 0.25, and is time non-dimensionalized with respect to fluid shear rate. Fig. 6 compares the rotational trajectories of free stream platelets flowing near the wall with the Jeffrey orbit prediction.



View larger version (25K):
[in this window]
[in a new window]
  		FIGURE 6  Plot of the trajectories of free-stream platelets as predicted by the Jeffery orbits theory, and platelet trajectories from experiments conducted at 0.2 dyn/cm2 shear stress. The experimental data ({circ}, *, {bigtriangledown}, and {bigtriangleup}) are compared to the Jeffery orbit model for spheroid rotation in unbounded fluid (dashed line), and the theoretical prediction corrected to include the wall effect on plate rotational motion (thick solid line).

 
The Jeffery orbit theory does not account for the retardation of the rotational velocity of spheroids due to the presence of a nearby wall. Kim et al. (2001)Go showed that the viscous torque resistance to the rotation of a circular disk when oriented parallel to the plane, at a dimensionless distance h/a = 1 from an infinite plane, is 1.53 times greater than that for a circular disk rotating in unbounded fluid. At a distance of h/a = 0.6, where a is the disk radius, the relative torque increases to 1.88. Considering a platelet with aspect ratio 0.25 and in near contact with a wall, the maximum separation between the platelet centroid and the wall will have value h/a = 1.0 (corresponding to {alpha} = {pi}/2), whereas the minimum separation will be h/a = 0.25 (corresponding to {alpha} = 0, {pi}). As a first approximation to inclusion of the influence of the wall we chose an intermediate value of h/a = 0.6 and scaled the freestream platelet rotational velocity by the corresponding factor of 1.88. Here, rescaling the rotational velocity is equivalent to stretching the time domain. Fig. 6 compares the calculated rotational motion with the experimentally measured orientation angles. Note that the corrected trajectory is based on the results for a disk oriented horizontally and therefore may not be valid for nonparallel platelet orientations. Considering these approximations, as shown in Fig. 6, the agreement with our experiments is quite good and suggests that the effect of the wall is of sufficient magnitude to account for the discrepancy with the unbounded fluid Jeffrey orbit theory. Thus, together with the results of Fig. 4, these data show that although unbound flowing platelets sample all possible orientations, they consistently attach only during the first-half of their rotation (orientation angle {alpha} = {pi} to {pi}/2).


   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
APPENDIX 1
 APPENDIX 2
 ACKNOWLEDGEMENTS
 REFERENCES
 
Over the past century, much theoretical and experimental work has elucidated the motion of spherical particles in Stokes flow for a variety of flow fields and bounding conditions. In the last four decades, significant progress has been made in the development of exact and approximate analytical solutions for hydrodynamic interactions in unbounded shear flow between two spheres (Lin et al., 1970Go; Batchelor and Green, 1972Go; Brenner and O'Neill, 1972Go) and between multiple spheres (Haber and Brenner, 1999Go). Dean and O'Neill (1963)Go, O'Neill (1964)Go, Brenner (1961)Go, and Goldman et al. (1967aGo,bGo) developed exact solutions using bi-spherical coordinates for the translational and rotational motion of a sphere in bounded Stokes flow near a plane wall, and Brenner and Happel (1958)Go developed approximations for determining flow characteristics of viscous flow past a sphere in a cylinder. Solutions for the flow of nonspherical particles are less commonly available mainly because of the difficulties in determining the hydrodynamic interactions between nonspherical particles or a nonspherical particle and a plane. Jeffery (1922)Go solved the hydrodynamic problem of the motion of an ellipsoidal particle in fluid that is motionless at infinity. Yoon and Kim (1990)Go used boundary collocation methods to solve the mobility problem of two spheroids in Stokes flow. However, we are aware of no studies that have considered the motion of an ellipsoid or oblate spheroid near a plane wall under linear shear flow, a scenario highly relevant to platelet flow and adhesion in blood vessels.

Computational modeling of blood cell adhesion is a highly useful tool to clarify the influence of physicochemical factors (hydrodynamic flows, hydrodynamic cell-cell interactions, chemical kinetics) on adhesion dynamics of cells under flow. A number of theoretical models have been developed to describe blood cell adhesion to the vascular wall and cell-cell aggregation at physiological flow rates (Hammer and Apte, 1992Go; Helmke et al., 1998Go; Tandon and Diamond, 1997Go, 1998Go; Long et al., 1999Go; King and Hammer, 2001aGo,bGo). These models approximate blood cells as spherical in shape. Computational models of platelet adhesion to a surface under flow are lacking due to the scarcity of theoretical studies on nonspherical particulate-wall interactions under flow. Since platelet-shaped cells behave differently from spherical cells under shear flow, the unique shape of the platelet brings new challenges in understanding the effect of hydrodynamic forces on its tether bond formation.

In this article, we have presented a clear two-dimensional analytical model to characterize the flipping of tethered platelets at various physiological shear rates. Transient tethering and translocation of platelets by means of flipping on the injured endothelial surface is a characteristic feature of GPIb{alpha}-vWF-A1 interactions (Doggett et al., 2002Go), and has been likened to the rolling of leukocytes. Our theoretical model compares well with the experimentally observed platelet motion on an adhesive surface. Additionally, this theoretical analysis estimates the hydrodynamic radial force acting on the platelet, and predicts the force to be compressive (directed toward the wall) during the first half of the tethered platelet's flip on the surface and tensile (directed away from the wall) during the second half. For the case of a platelet freely flowing over the surface, the rotational motion of the platelet caused by shearing fluid forces results in the platelet edges being periodically pushed toward and pulled away from the surface. Note that a rigid spherical body does not experience a lateral force when flowing parallel to a planar wall in shear Stokes flow (Goldman et al., 1967aGo), whereas deformable droplets, on the other hand, experience a lift force away from the wall even at zero Reynolds number (Leal, 1980Go); such flow behaviors may be expected to influence contacting of cells with the wall. Finally, we analyzed the orientation of platelets at the time of binding to, and release from, the surface, and showed that all of the experimental observations of bond formation and bond dissociation correlate with the predicted sign of the hydrodynamic radial forces acting on the platelet, i.e., tethers are only observed to form when the predicted hydrodynamic force is compressive, whereas tethers are only observed to dissociate when the force is tensile.

Our findings suggest that the pattern of platelet binding to, and release from, a vWF-coated surface is governed by radial forces experienced by the platelet. Since we are unable to obtain a direct, experimental measurement of the hydrodynamic force acting on the platelet at the time of cell-surface binding, we must rely on the theoretical model's estimate of this force, which is derived through solution of the creeping flow equations and shown to accurately predict the instantaneous velocity of the observed flipping platelets. Tether bond formation appears to occur only in those orientations where a hydrodynamic compressive force is predicted to act along the platelet length, pushing it against the surface, thereby initiating cell adhesion. Compression promotes bond formation by bringing the surfaces into close molecular contact and helping to overcome repulsive long-range interactions between the surfaces (Bell et al., 1984Go). After initial bond formation, this force aids in maintaining close contact between the cell surfaces and strengthens the initial tether by promoting further bond formation.

Platelet-endothelial collision and resulting contact is expected to occur only under conditions of hydrodynamic compressive force on the platelet. Figs. 4 and 5 clearly demonstrate that although shear flow promotes contact and subsequent attachment at angles of incidence >90°, the angle of attachment is not found to be a function of shear rate. As explained earlier, this result is expected because particle trajectories are independent of the magnitude of shear rate in Stokes flow. Compressive force on the platelet, on the other hand, is a function of both the incident angle and shear rate. If a critical magnitude of compressive force was required to promote bond formation, one might expect to see a trend in the angles of bond formation with shear rate. Here, it is important to note that once the platelet collides with the surface at an angle (>90°), it remains in contact until it assumes a 90° orientation. Although the GPIb{alpha}-vWf-A1 bond formation kinetics have yet to be fully characterized, we speculate that a finite time of application of this compressive force may induce bond formation at smaller angles (but >90°) associated with less favorable compressive forces. It is therefore difficult to ultimately conclude the essentiality or non-essentiality of a minimum compressive force for binding from these data alone.

Although it is not obvious from the experiments and theoretical model whether cell compression or cell contact is directly responsible for the good agreement between the experimental observations of tether formation and dissociation and the theoretical calculation of hydrodynamic forces, the conclusion that compressive forces play a role in initiating tether formation is nevertheless consistent with the experimental and theoretical results, and is indirectly supported by these results. It has been shown that platelet tethering on immobilized vWf increases with increase in shear rate (Doggett et al., 2002Go), which implies that tether bond formation is positively influenced by shear rate and hence compressive hydrodynamic forces. The role of compressive forces in initiating receptor-ligand binding is a fundamental question in blood cell adhesion. Our results suggest the importance of sub-picoNewton compressive forces in the initial capture of human platelets to vWF under physiological flow conditions. Extending our findings to other families of bonds, there is a possibility that leukocyte selectin-ligand bonds may associate only under suitable conditions of compressive force, which is consistent with recent studies of leukocyte bond formation using micromanipulation techniques (Chesla et al., 1998Go).

The GPIb{alpha}-vWF-A1 tether bond requires a threshold level of fluid shear stress (73 s–1) (Doggett et al., 2002Go) to promote adhesive interactions. Below this shear threshold, platelets cease to interact with the immobilized vWF-A1 on the surface and continue to flow at hydrodynamic velocity. This similar property is also shown by L-selectin interactions with its sialylated glycoprotein ligand (shear threshold ~40 s–1) (Finger et al., 1996Go). The mechanism by which shearing forces govern vWF-GPIb{alpha} binding is unknown. Interestingly, whereas GPIb{alpha} platelet receptors do not bind to circulating vWF multimeric plasma proteins under normal conditions, pathological high-shear stress conditions (above 60 dyn/cm2) induces binding of large circulating vWF to GPIb{alpha} (Konstantopoulos et al., 1997Go; Li et al., 2004Go) even in the absence of exogenous agonists or activating agents. Shear-induced GPIb{alpha}-vWF binding subsequently initiates platelet activation and aggregation. GPIb{alpha} cannot bind either surface-immobilized or circulating vWF under static conditions, which suggests that the mechanism of shear stress-induced binding is the same for both types of interactions. Our model (Eq. A31) shows that the radial compressive force acting on the platelet increases linearly with increasing shear rate, which implies that a critical shear threshold is equivalent to a critical normal force. According to our theory, the threshold shear rate of 73 s–1 for GPIb{alpha}-vWF-A1 adhesive interactions translates to a critical normal force of only 0.93 pN. By gaining a better understanding of the nature of hydrodynamic forces during cell adhesion, it will be possible to elucidate the link between the critical shear threshold and bond formation of GPIb{alpha} with vWF, a prime motivational factor for this study.

The hydrodynamic forces exerted on leukocyte bonds, during transient tethering or rolling on a reactive surface, can be determined by carrying out a simple force and torque balance on the model spherical leukocyte. This technique was first employed by Alon et al. (1995)Go for calculating the first-order kinetic parameters of the P-selectin-PSGL-1 bond, and again (1997) for studying L-selectin tethers. Shao et al. (1998)Go and Park et al. (2002)Go used this bond force calculation method for studying microvilli extension. Doggett et al. (2002)Go used this same force balance technique to determine the kinetic parameters of the GPIb{alpha}-vWF-A1 bond. They circumvented the problem of the platelet being nonspherical by coating microspheres with vWF-A1 and then allowing them to flow over surface-immobilized platelets. It is important to note that the majority of platelet bonds with the surface, unlike leukocyte tethers, will not experience a sizeable lever-arm effect due to distinct differences in the geometry of these two adhesive systems. In the case of a spherical leukocyte binding with the surface, a lever arm exerts force and torque on the cell at a location several microns upstream of the lowest point of the spherical cell body (Alon et al., 1995Go; Park et al., 2002Go). However, platelet bonds are usually located close to the point of platelet-surface contact (the so-called hinge region) while flipping. Additionally, due to the ellipsoidal shape of the platelet, normal forces cause the platelet to flip or rotate over the surface (which is not possible for tethered spherical cells) and so the lever effect is minimal here. Therefore, we assume that the force on the bond is of the same order of magnitude as the hydrodynamic radial force acting on the platelet.

A key element in these force and torque balances is the hydrodynamic shear force and torque that the cell experiences when stationary on a planar surface. These are provided by the solution for an immobilized sphere under linear shear flow contacting a plane wall, developed by Goldman et al. (1967b)Go. However, as mentioned earlier, there are no existing solutions for the force exerted on a nonspherical cell tethered to a surface. What is of prime importance in any physiological model for blood cell adhesion is an accurate estimation of the forces acting on the cell near, or on, a surface. Without this knowledge it becomes difficult to predict the flow behavior of the cell, the behavior of bonds when the cell is tethered, and how quickly bonds will form or dissociate based on the range of forces acting on the bonds and cell. Such analysis of platelet tethering behavior is possible only if the hydrodynamic effects on platelet-shaped cells are well defined. For the case of a platelet tethered to, and translocating on a surface, the hydrodynamic forces acting on it will depend on its shape, size, and orientation. In the present article, we estimated the platelet as a flat plate, and proceeded to calculate the forces and torques acting on the plate for all orientations. Therefore, our results can be viewed as first estimates of the hydrodynamic force and torque on a platelet tethered to a surface, which can be compared with the observed platelet adhesive behavior to gain a better understanding of hydrodynamic effects on platelet adhesive phenomena.

Although the model presented here merits simplicity and provides insight into the force mechanics acting on the tethered platelet during its flipping motion on the surface, its utility is nevertheless limited since it does not consider binding kinetics and simplifies the platelet shape as a thin two-dimensional planar surface. A logical next step would be to develop a numerical simulation that examines the full three-dimensional motion of platelets during approach, collision, adhesion, and release from the vessel wall in order to extend the range of validity of the analytical model, and further define the precise force history on the GPIb{alpha}-vWF-A1 bond during tethering and translocation.


   APPENDIX 1
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX 1
APPENDIX 2
 ACKNOWLEDGEMENTS
 REFERENCES
 
Jeong and Kim (1983)Go derived the analytical expressions for forces and moments due to linear shear flow acting on a stationary fence with unit length, mounted at an arbitrary angle {alpha} with respect to an infinite rigid wall (Fig. 2). The flow around the fence was assumed to be slow and viscous so that the inertial terms are negligible and flow may be approximated as Stokes flow. If are unit vectors in the x and z directions respectively, the governing equations describing the flow are

(A1)

(A2)

(A3)

The Stokes equation and the continuity equation along with the above boundary conditions were solved by describing pressure and velocity as vector harmonic functions, transforming the resulting equations by Mellin transformation and then subsequently solving a set of simultaneous Wiener-Hopf equations in four unknowns using the technique described by Noble (1958)Go. The normal force Fn, shear force Fr, and moment (or torque) Tf exerted on the fence by the fluid were calculated using the following relations:

(A4)

(A5)

(A6)
We found that the radial force Fr (Eq. A5) can be approximated (with maximum error of 3.6% and average error of 1.76%) by a polynomial fitted to the Jeong-Kim plot of analytically determined tangential force as a function of inclination angle {alpha}, as

(A7)
where A is the area of one face of the plate.

The torque Tf generated by shear flow can be similarly approximated (with maximum error of 3.5% and average error of 1.1%) by a polynomial fitted to the Jeong-Kim plot of analytically determined torque as a function of inclination angle {alpha}, as

(A8)
where is the height of the plate.

Fig. 7 demonstrates the good fit of the approximating polynomials determined by a least-squares regression fit to the analytically derived solution of force (Eq. A7) and torque (Eq. A8).



View larger version (15K):
[in this window]
[in a new window]
  		FIGURE 7  Plots of the analytically derived solution for force and torque acting on a stationary inclined fence subjected to linear shear flow (fence problem; Jeong and Kim, 1983Go), and the polynomial approximations fitted to the analytical solution by least-squares regression. (A) Plot of force values taken from the analytically derived solution (Eq. A5) ({diamondsuit}) and force calculated from the polynomial approximation (Eq. A7) (thick solid line). (B) Plot of torque values taken from the analytically derived solution (Eq. A6) ({diamondsuit}) and torque calculated from the polynomial approximation (Eq. A8) (thick solid line).

 

   APPENDIX 2
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX 1
 APPENDIX 2
ACKNOWLEDGEMENTS
 REFERENCES
 
We derived the analytical expressions for the force and torque acting on a two-dimensional plate that rotates at an angular velocity {omega}, and is hinged to a rigid infinite wall, with the surrounding fluid being motionless at infinity. Fluid motion in the vorticity direction is neglected. The Stokes approximation is valid in general for flow of a viscous fluid near a sharp corner between two planes. Because of the microscopic dimensions of our fluid system, Stokes flow prevails everywhere in the fluid. Fluid motion on either side of the rotating plate is considered separately in our derivations. The full solution is then equal to the sum of both flows.

The flow due to the motion of the plate is unsteady and therefore the acceleration term in the equation of motion for a viscous incompressible fluid, does not vanish. On non-dimensionalizing the time-dependent Stokes equation, we obtain where NRE (Reynolds Number) = L2{rho}{omega}/µ. Since for the experimental conditions studied NRE << 1, the flow is in the Stokes regime and we may assume quasi-steady-state flow.

Hinge flow in the region downstream of the rotating plate
The surface at {theta} = 0 is stationary and the surface at {theta} = {alpha} moves at an angular velocity of {omega}(= –d{alpha}/dt) toward the horizontal surface (Fig. 2). Both angle {alpha} and angular velocity {omega} are functions of time. The derived results are valid at an instantaneous time t. This problem was solved using stream functions. The boundary conditions for this problem include the no-slip and no-penetration conditions at the solid surfaces, {theta} = 0, {theta} = {alpha}, and are expressed as

(A9)

Solving the Stokes flow equation we obtain the non-dimensionalized velocity components as

(A10a)

(A10b)
where the variables are non-dimensionalized, as shown below, and where is the length of the hinged plate,

The pressure dependence on r,{theta} is determined by substituting Eq. A10 into the Navier-Stokes equations in cylindrical coordinates and integrating, yielding the result

(A11)
where C1 is an integration constant. Note that the pressure is independent of {theta}.

The r- and {theta}-components of the dimensionless drag force can be evaluated from the expressions

(A12a)

(A12b)
where is the viscous stress tensor, and is the pressure as determined above. The radial drag force was calculated as

(A13)
which, evaluated at the planar surface, gives Similarly, the {theta}-component of the force is obtained as

which, when evaluated at the planar surface, yields

(A14)

Note that the viscous stresses do not contribute to drag force on the hinged plate. Also, there is no radial force on the plate due to the hinge flow.

Hinge flow in the region upstream of the rotating plate
Next we considered the fluid motion induced on the other side of the rotating plate. The surface at {theta} = {pi} is stationary and the surface at {theta} = {alpha} moves at an angular velocity of {omega}(= –d{alpha}/dt) away from the horizontal surface (Fig. 2). The solution of this problem is solved using the same methodology as in Hinge Flow in the Region Downstream of the Rotating Plate, above. The non-dimensionalized pressure and drag forces are calculated as

(A15)
where C2 is an integration constant,

(A16)

and

(A17)

(A18)

Evaluation of torque on the rotating plate
The torque on the plate—generated by pressure forces—which tends to oppose motion of the hinged plate, is calculated by substituting Eqs. A11 and A15 in

(A19)
to yield

(A20)

From the constraints of symmetry in Stokes flow for the geometry under consideration, C = C1 = –C2, where C is a constant taken to be independent of {alpha}, and Eq. A19 becomes

(A21)

Determination of plate angle with respect to the surface as a function of time using combined fence and hinge flow solutions
The torques (calculated from the fence and hinge problems) acting on the plate are equated (Tf = Th) to obtain an expression for angular velocity {omega} = –d{alpha}/dt of the plate. The steady-state torque balance incorporating fence flow (Eq. A8) and hinge flow solutions (Eq. A21) yields an expression for {omega},

which is rearranged to give

(A22)

Numerical integration of Eq. A22 with respect to time yields the theoretical dependence of plate angle {alpha} (platelet orientation {alpha}) with time as the hinged plate (platelet) rotates (flips) with respect to the rigid wall under an imposed linear shear flow. Changing the single adjustable parameter C, representing the unknown pressure integration constant, results in a change in the slope of the curve. C takes on a value of 30 in our solution, which produces good agreement between the theoretical prediction and the experimental observations. To verify the Stokes approximation, NRE is calculated for the largest value of {omega}, which evaluated from Eq. A22 is 0.133{gamma}. Substituting into the dimensionless group L2{rho}{omega}/µ, one obtains 0.133 {gamma} L2{rho}/µ. Substituting pertinent values for these physical parameters, L = 3 µm, {rho} = 1 g/cm3, µ = 1 cp, and shear rate = 100–800 s–1, the range of values derived are 1.2 x 10–4 to 9.6 x 10–4 which are << 1. Thus the Stokes flow and quasi-steady-state assumption remain valid.


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX 1
 APPENDIX 2
 ACKNOWLEDGEMENTS
REFERENCES
 
This work has been supported by a grant from the Whitaker Foundation to M.R.K., and a National Institutes of Health grant HL063244 and an American Heart Association grant 02-40009N to T.G.D.

Submitted on June 2, 2004; accepted for publication November 1, 2004.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX 1
 APPENDIX 2
 ACKNOWLEDGEMENTS
 REFERENCES
 
Alevriadou, B. R., J. L. Moake, N. A. Turner, Z. M. Ruggeri, B. J. Folie, M. D. Phillips, A. B. Schreiber, M. E. Hrinda, and L. V. McIntire. 1993. Real-time analysis of shear-dependent thrombus formation and its blockade by inhibitors of von Willebrand factor binding to platelets. Blood. 81:1263–1276.[Abstract/Free Full Text]

Alon, R., D. A. Hammer, and T. A. Springer. 1995. Lifetime of the P-selectin carbohydrate bond and its response to tensile force in hydrodynamic flow. Nature. 374:539–542.[CrossRef][Medline]

Alon, R., S. Chen, K. D. Puri, E. B. Finger, and T. A. Springer. 1997. The kinetics of L-selectin tethers and the mechanics of selectin-mediated rolling. J. Cell Biol. 138:1169–1180.[Abstract/Free Full Text]

Batchelor, G. K., and J. T. Green. 1972. The hydrodynamic interaction of two small freely-moving spheres in a linear flow field. J. Fluid Mech. 56:375–400.[CrossRef]

Bell, G. I., M. Dembo, and P. Bongrand. 1984. Cell adhesion. Competition between non-specific repulsion and specific bonding. Biophys. J. 45:1051–1064.[Abstract/Free Full Text]

Brenner, H., and J. Happel. 1958. Slow viscous flow past a sphere in a cylindrical tube. J. Fluid Mech. 4:195–213.[CrossRef]

Brenner, H. 1961. The slow motion of a sphere through a viscous fluid towards a plane surface. Chem. Eng. Sci. 16:242–251.[CrossRef]

Brenner, H., and M. E. O'Neill. 1972. On the Stokes resistance of multiparticle systems in a linear shear field. Chem. Eng. Sci. 27:1421–1439.[CrossRef]

Chesla, S. E., P. Selvaraj, and C. Zhu. 1998. Measuring two-dimensional receptor-ligand binding kinetics by micropipette. Biophys. J. 75:1553–1572.[Abstract/Free Full Text]

Coller, B. S., E. I. Peerschke, L. E. Scudder, and C. A. Sullivan. 1983. Studies with a murine monoclonal antibody that abolishes ristocetin-induced binding of von Willebrand factor to platelets: additional evidence in support of GPIb as a platelet receptor for von Willebrand factor. Blood. 61:99–110.[Abstract/Free Full Text]

Cruz, M. A., T. G. Diacovo, J. Emsley, R. Liddington, and R. I. Handin. 2000. Mapping the glycoprotein Ib-binding site in the von Willebrand factor A1 domain. J. Biol. Chem. 275:252–255.

Dean, W. R., and M. E. O'Neill. 1963. A slow motion of viscous liquid caused by the rotation of a solid sphere. Mathematika. 10:13–24.

Doggett, T. A., G. Girdhar, A. Lawshé, D. W. Schmidtke, I. J. Laurenzi, S. L. Diamond, and T. G. Diacovo. 2002. Selectin-like kinetics and biomechanics promote rapid platelet adhesion in flow: the GPIb{alpha}-vWF tether bond. Biophys. J. 83:194–205.[Abstract/Free Full Text]

Doggett, T. A., G. Girdhar, A. Lawshe, J. L. Miller, I. J. Laurenzi, S. L. Diamond, and T. G. Diacovo. 2003. Alterations in the intrinsic properties of the GPIb{alpha}-VWF tether bond define the kinetics of the platelet-type von Willebrand disease mutation, Gly233Val. Blood. 102:152–160.[Abstract/Free Full Text]

Du, X., and M. H. Ginsberg. 1997. Integrin {alpha} IIb ß3 and platelet function. Thromb. Haemost. 78:96–100.[Medline]

Finger, E. B., K. D. Puri, R. Alon, M. B. Lawrence, U. H. von Andrian, and T. A. Springer. 1996. Adhesion through L-selectin requires a threshold hydrodynamic shear. Nature. 379:266–269.[CrossRef][Medline]

Frojmovic, M., K. Longmire, and T. G. M. Van de Ven. 1990. Long-range interactions in mammalian platelet aggregation. II. The role of platelet pseudopod number and length. Biophys. J. 58:309–318.[Abstract/Free Full Text]

Goldman, A. J., R. G. Cox, and H. Brenner. 1967a. Slow viscous motion of a sphere parallel to a plane wall. I. Motion through a quiescent fluid. Chem. Eng. Sci. 22:637–651.[CrossRef]

Goldman, A. J., R. G. Cox, and H. Brenner. 1967b. Slow viscous motion of a sphere parallel to a plane wall. II. Couette flow. Chem. Eng. Sci. 22:653–660.[CrossRef]

Grunemeier, J. M., W. B. Tsai, C. D. McFarland, and T. A. Horbett. 2000. The effect of adsorbed fibrinogen, fibronectin, von Willebrand factor and vitronectin on the procoagulant state of adherent platelets. Biomaterials. 21:2243–2252.[CrossRef][Medline]

Haber, S., and H. Brenner. 1999. Hydrodynamic interactions of spherical particles in quadratic Stokes flows. Intl. J. Multiphase Flow. 25:1009–1032.[CrossRef]

Hammer, D. A., and S. M. Apte. 1992. Simulation of cell rolling and adhesion on surfaces in shear flow: general results and analysis of selectin-mediated neutrophil adhesion. Biophys. J. 62:35–57.

Helmke, B. P., M. Sugihara-Seki, R. Skalak, and G. W. Schmid-Schönbein. 1998. A mechanism for erythrocyte-mediated elevation of apparent viscosity by leukocytes in vivo without adhesion to the endothelium. Biorheology. 35:437–448.[CrossRef][Medline]

Ikeda, Y., M. Handa, K. Kawano, T. Kamata, M. Murata, Y. Araki, H. Anbo, Y. Kawai, K. Watanabe, I. Itagaki, K. Sakai, and Z. M. Ruggeri. 1991. The role of von Willebrand factor and fibrinogen in platelet aggregation under varying shear stress. J. Clin. Invest. 87:1234–1240.[Medline]

Jeffery, G. B. 1922. The motion of ellipsoidal particles immersed in a viscous fluid. Proc. Roy. Soc. Phys. A. 102:161–179.

Jeong, J., and M. Kim. 1983. Slow viscous flow around an inclined fence on a plane. J. Phys. Soc. Jpn. 52:2356–2363.[CrossRef]

Kasirer-Friede, A., M. R. Cozzi, M. Mazzucato, L. De Marco, Z. M. Ruggeri, and S. J. Shattil. 2004. Signaling through GP Ib-IX-V activates {alpha} IIb ß3 independently of other receptors. Blood. 103:3403–3411.[Abstract/Free Full Text]

Kim, M.-U., K. W. Kim, Y.-H. Cho, and B. M. Kwak. 2001. Hydrodynamic force on a plate near the plane wall. I. Plate in sliding motion. Fluid Dyn. Res. 29:137–170.[CrossRef]

King, M. R., and D. A. Hammer. 2001a. Multiparticle adhesive dynamics: hydrodynamic recruitment of rolling leukocytes. Proc. Natl. Acad. Sci. USA. 98:14919–14924.[Abstract/Free Full Text]

King, M. R., and D. A. Hammer. 2001b. Multiparticle adhesive dynamics. Interactions between stably rolling cells. Biophys. J. 81:799–813.[Abstract/Free Full Text]

Konstantopoulos, K., T. W. Chow, N. A. Turner, J. D. Hellums, and J. L. Moake. 1997. Shear stress-induced binding of von Willebrand factor to platelets. Biorheology. 34:57–71.[CrossRef][Medline]

Kroll, M. H., J. D. Hellums, L. V. McIntire, A. I. Schafer, and J. L. Moake. 1996. Platelets and shear stress. Blood. 88:1525–1541.[Free Full Text]

Kumar, R. A., J. F. Dong, J. A. Thaggard, M. A. Cruz, J. A. López, and L. V. McIntire. 2003. Kinetics of GPIb{alpha}-vWF-A1 tether bond under flow: effect of GPIb{alpha} mutations on the association and dissociation rates. Biophys. J. 85:4099–4109.[Abstract/Free Full Text]

Leal, L. G. 1980. Particle motion in a viscous fluid. Ann. Rev. Fluid Mech. 12:435–476.[CrossRef]