Platelet accumulation at the site of a vascular injury is regulated by soluble platelet agonists, which induce various types of platelet responses, including integrin activation and granule secretion. The interplay between local biochemical cues, mechanical interactions between platelets and macroscopic thrombus dynamics is poorly understood. Here we describe a novel computational model of microvascular clot formation for the detailed analysis of thrombus mechanics. We adopt a previously developed two-dimensional particle-based model focused on the thrombus shell formation and revise it to introduce the platelet agonists. Blood flow is simulated via a computational fluid dynamics approach. In order to model soluble platelet activators, we apply Langevin dynamics to a large number of non-dimensional virtual particles. Taking advantage of the available data on platelet dense granule secretion kinetics, we model platelet degranulation as a stochastic agonist-dependent process. The new model qualitatively reproduces the enhanced thrombus formation due to dense granule secretion, in line with in vivo findings, and provides a mechanism for the thrombin confinement at the early stages of clot formation. Our calculations also predict that the release of platelet dense granules results in the additional mechanical stabilization of the inner layers of thrombus. Distribution of the inter-platelet forces throughout the aggregate reveals multiple weak spots in the outer regions of a thrombus, which are expected to result in the mechanical disruptions at the later stages of clot formation. © 2021 Elsevier Ltd