Monte Carlo simulation of the heterotypic aggregation kinetics of platelets and neutrophils.

The heterotypic aggregation of cell mixtures or colloidal particles such as proteins occurs in a variety of settings such as thrombosis, immunology, cell separations, and diagnostics. Using the set of population balance equations (PBEs) to predict dynamic aggregate size and composition distributions is not feasible. The stochastic algorithm of Gillespie for chemical reactions (. J. Comput. Phys. 22:403-434) was reformulated to simulate the kinetic behavior of aggregating systems. The resulting Monte Carlo (MC) algorithm permits exact calculation of the decay rates of monomers and the temporally evolving distribution of sizes and compositions of the aggregates. Moreover, it permits calculation of all moments of these distributions. Using this method, we explored the heterotypic aggregation of fully activated platelets and neutrophils in a linear shear flow of shear rate G = 335 s(-1). At plasma concentrations, the half-lives of homotypically aggregating platelet and neutrophil singlets were 8.5 and 2.4 s, respectively. However, for heterotypic aggregation, the half-lives for platelets and neutrophils decreased to 2.0 and 0.11 s, respectively, demonstrating that flowing neutrophils accelerate capture of platelets and growth of aggregates. The required number of calculations per time step of the MC algorithm was typically a small fraction of Omega(1/2), where Omega is the initial number of particles in the system, making this the fastest MC method available. The speed of the algorithm makes feasible the deconvolution of kernels for general biological heterotypic aggregation processes.

[1]  Themis Matsoukas,et al.  Constant-number Monte Carlo simulation of population balances , 1998 .

[2]  Daniel T. Gillespie,et al.  The Stochastic Coalescence Model for Cloud Droplet Growth. , 1972 .

[3]  S. Diamond,et al.  Kinetics of beta2-integrin and L-selectin bonding during neutrophil aggregation in shear flow. , 1998, Biophysical journal.

[4]  Daniel T. Gillespie,et al.  An Exact Method for Numerically Simulating the Stochastic Coalescence Process in a Cloud , 1975 .

[5]  J. Ramírez,et al.  Platelet activation In mice and human Helicobacter pylori infection. , 1997, The Journal of clinical investigation.

[6]  A. Perelson,et al.  Kinetics of rouleau formation. I. A mass action approach with geometric features. , 1982, Biophysical journal.

[7]  F. Neumann,et al.  Increased neutrophil-platelet adhesion in patients with unstable angina. , 1996, Circulation.

[8]  S Neelamegham,et al.  Modeling the reversible kinetics of neutrophil aggregation under hydrodynamic shear. , 1997, Biophysical journal.

[9]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[10]  S Neelamegham,et al.  Venous levels of shear support neutrophil-platelet adhesion and neutrophil aggregation in blood via P-selectin and beta2-integrin. , 1998, Circulation.

[11]  H. Goldsmith,et al.  Adenosine diphosphate-induced aggregation of human platelets in flow through tubes. II. Effect of shear rate, donor sex, and ADP concentration. , 1989, Biophysical journal.

[12]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[13]  G. Benedek,et al.  Experimental measurements of the temporal evolution of cluster size distributions for high-functionality antigens cross-linked by antibody , 1983 .

[14]  C. DeLisi.  Theory of clustering of cell surface receptors by ligands of arbitrary valence: dependence of dose response patterns on a coarse cluster characteristic , 1980 .

[15]  Doraiswami Ramkrishna,et al.  Simulation of particulate systems using the concept of the interval of quiescence , 1977 .

[16]  J. D. Hellums,et al.  Analysis of shear-induced platelet aggregation with population balance mathematics. , 1986, Biophysical journal.

[17]  R. McEver,et al.  GMP-140 mediates adhesion of stimulated platelets to neutrophils. , 1990, Blood.

[18]  J. O'brien,et al.  Shear-induced platelet aggregation , 1990, The Lancet.

[19]  Y. Matsuzawa,et al.  Endogenous adenosine inhibits P-selectin-dependent formation of coronary thromboemboli during hypoperfusion in dogs. , 1998, The Journal of clinical investigation.

[20]  Doraiswami Ramkrishna,et al.  A puristic analysis of population balance—II. , 1973 .

[21]  M. Hounslow,et al.  A discretized population balance for nucleation, growth, and aggregation , 1988 .

[22]  M Long,et al.  Probabilistic modeling of shear-induced formation and breakage of doublets cross-linked by receptor-ligand bonds. , 1999, Biophysical journal.

[23]  F. W. Wiegel,et al.  The equilibrium size distribution of rouleaux. , 1982, Biophysical journal.

[24]  M. Bednar,et al.  Neutrophil Depletion Suppresses 111In-Labeled Platelet Accumulation in Infarcted Myocardium , 1985, Journal of cardiovascular pharmacology.

[25]  D. Altieri,et al.  Activation of Mac-1 (CD11b/CD18)-bound factor X by released cathepsin G defines an alternative pathway of leucocyte initiation of coagulation. , 1996, The Biochemical journal.

[26]  AlbertSchömig,et al.  Induction of Cytokine Expression in Leukocytes by Binding of Thrombin-Stimulated Platelets , 1997 .

[27]  J. Mathew,et al.  Cardiopulmonary bypass induces leukocyte-platelet adhesion. , 1992, Blood.

[28]  H. Goldsmith,et al.  Adenosine diphosphate-induced aggregation of human platelets in flow through tubes. I. Measurement of concentration and size of single platelets and aggregates. , 1989, Biophysical journal.

[29]  W. T. Scott,et al.  Analytic Studies of Cloud Droplet Coalescence I , 1968 .

[30]  A S Perelson,et al.  Aggregation of cell surface receptors by multivalent ligands , 1982, Journal of mathematical biology.

[31]  C. Benjamin,et al.  Leukocyte accumulation promoting fibrin deposition is mediated in vivo by P-selectin on adherent platelets , 1992, Nature.

[32]  J D Hellums,et al.  Aggregation and disaggregation kinetics of human blood platelets: Part I. Development and validation of a population balance method. , 1993, Biophysical journal.

[33]  S Neelamegham,et al.  Molecular dynamics of the transition from L-selectin- to beta 2-integrin-dependent neutrophil adhesion under defined hydrodynamic shear. , 1996, Biophysical journal.

[34]  E. Ponder,et al.  ON SEDIMENTATION AND ROULEAUX FORMATION—II , 1925 .

[35]  C. Cerletti,et al.  Platelet/polymorphonuclear leukocyte interaction in dynamic conditions: evidence of adhesion cascade and cross talk between P-selectin and the beta 2 integrin CD11b/CD18. , 1996, Blood.

[36]  S L Diamond,et al.  Hydrodynamic effects and receptor interactions of platelets and their aggregates in linear shear flow. , 1997, Biophysical journal.

[37]  K. Brown,et al.  Neutrophil-platelet adhesion: relative roles of platelet P-selectin and neutrophil beta2 (DC18) integrins. , 1998, American journal of respiratory cell and molecular biology.

[38]  J D Hellums,et al.  Aggregation and disaggregation kinetics of human blood platelets: Part II. Shear-induced platelet aggregation. , 1993, Biophysical journal.

[39]  Alan S. Perelson,et al.  Receptor clustering on a cell surface. III. theory of receptor cross-linking by multivalent ligands: description by ligand states , 1981 .

[40]  F. Neumann,et al.  Induction of cytokine expression in leukocytes by binding of thrombin-stimulated platelets. , 1997, Circulation.