Introduction: We developed a multiscale model to simulate the dynamics of platelet aggregation by recruitment of unactivated platelets flowing in viscous shear flows by an activated platelet deposited onto a blood vessel wall. This model uses coarse grained molecular dynamics (CGMD) for platelets at the microscale and dissipative particle dynamics (DPD) for the shear flow at the macroscale. Under conditions of relatively low shear, aggregation is mediated by fibrinogen via αIIbβ3 receptors.

Methods: The binding of αIIbβ3 and fibrinogen is modeled by a molecular-level hybrid force field consisting of Morse potential and Hooke law for the nonbonded and bonded interactions, respectively. The force field, parametrized in two different interaction scales, is calculated by correlating with the platelet contact area measured in vitro and the detaching force between αIIbβ3 and fibrinogen.

Results: Using our model, we derived, the relationship between recruitment force and distance between the centers of mass of two platelets, by integrating the molecular-scale inter-platelet interactions during recruitment aggregation in shear flows. Our model indicates that assuming a rigid-platelet model, underestimates the contact area by 89% and the detaching force by 93% as compared to a model that takes into account the platelet deformability leading to a prediction of a significantly lower attachment during recruitment.

Conclusions: The molecular-level predictive capability of our model sheds a light on differences observed between transient and permanent platelet aggregation patterns. The model and simulation framework can be further adapted to simulate initial thrombus formation involving multiple flowing platelets as well as deposition and adhesion onto blood vessels.