## Cooperative SIR dynamics as a model for spontaneous blood clot initiation

## Philip Greulich a

a Mathematical Sciences, University of Southampton, Southampton, SO17 1BJ, United Kingdom b Institute for Life Sciences, University of Southampton, Southampton, SO17 1BJ, United Kingdom

## Abstract

Blood clotting is an important physiological process to suppress bleeding upon injury, but when it occurs inadvertently, it can cause thrombosis, which can lead to life threatening conditions. Hence, understanding the microscopic mechanistic factors for inadvertent, spontaneous blood clotting, in absence of a vessel breach, can help in predicting and adverting such conditions. Here, we present a minimal model - reminiscent of the SIR model - for the initiating stage of spontaneous blood clotting, the collective activation of blood platelets. This model predicts that in the presence of very small initial activation signals, macroscopic activation of the platelet population requires a su ffi cient degree of heterogeneity of platelet sensitivity. To propagate the activation signal and achieve collective activation of the bulk platelet population, it requires the presence of, possibly only few, hyper-sensitive platelets, but also a su ffi cient proportion of platelets with intermediate, yet higher-than-average sensitivity. A comparison with experimental results demonstrates a qualitative agreement for high platelet signalling activity.

Keywords: blood clotting, heterogeneity, SIR model, collective phenomena

## 1. Introduction

Blood clotting is the formation of macroscopic aggregates of blood platelets which are stabilised by a fibrin polymer network (coagulation) (Gale, 2011). The physiological purpose of such an aggregate, called a blood clot , or thrombus , is to seal breached vessels and to protect an individual from blood loss. However, when blood clots emerge in the blood stream, away from vessel breaches, they can occlude vessels ( thrombosis ) and obstruct blood flow, which is associated with pathologies such as stroke and heart attacks (Raskob et al., 2014), leading to more than 15 million deaths per annum world-wide (British Heart Foundation, 2021).

The formation of a platelet plug through platelet aggregation is the first stage of platelet-mediated blood clotting; this is then followed by the coagulation cascade that involves the development of a fibrin network to stabilise the platelet plug and form a mature blood clot (Gale, 2011). Platelets can only aggregate once they have been activated, either directly by agonists in the bloodstream and from vessel walls (usually arising from damaged vasculature) or via paracrine signalling from other activated platelets (Stalker et al., 2012; van der Meijden and Heemskerk, 2019).

It has been shown that platelets are highly diverse, and their sensitivity -the threshold and propensity to activate when exposed to stimulating agonists - di ff ers greatly between platelets within a person (van der Meijden and Heemskerk, 2019; Heemskerk and West, 2022; Jongen et al., 2020). Some platelets are extremely sensitive (Baaten et al., 2017), and it has been hypothesised that activation of such hyper-sensitive platelets -which might be rare, but can activate at very low stimulant concentrations - can propagate activation via paracrine signalling to activate also less sensitive platelets, thereby generating a large population of activated platelets. This could initiate platelet aggregation in situations where the abundance of stimulants is much lower than the threshold for the bulk activation of platelets (Baaten et al., 2017; Jongen et al., 2020; Lesyk and Jurasz, 2019). Thus, the diversity of platelet sensitivity, and not bulk platelet sensitivity alone, may drive platelet aggregation.

Physiological blood clotting occurs at vessel breaches, but blood clots can sometimes emerge without apparent vascular damage and with stimulant levels far below threshold levels expected to activate platelets ( hypercoagulation ) (LaPelusa and Dave, 2023). Although genetic and lifestyle risk factors for this are well established, the microscopic mechanistic origin of the emergence of such spontaneous blood clotting has so far eluded our full understanding, and quantitative modelling approaches could not reproduce or predict spontaneous blood clotting.

Whether activated clusters remain small or the activation signal percolates through the platelet population, leading to collective activation aggregation, is a non-trivial critical emergent / phenomenon which has not been understood yet. Although models for blood clot growth have been widely used, those established models usually model blood clot growth in the presence of a vessel breach, which secrets strong activation stimulants and where shear flow can trigger activation (Fogelson, 1993; Wang and Fogelson, 1999; Fogelson and Neeves, 2015; Fogelson, 2016; Schoeman et al., 2017; Link et al., 2020). On the other hand, no model has so far been able to explain how tiny stimulant signals can lead to macroscopic activation.

In this work, we present a minimal mechanistic model for paracrine (platelet-platelet) activation via signalling, which captures the main qualitative features of the first phase of spon-

taneous blood clotting - the macroscopic, collective activation of platelet population, upon a weak stimulant signal, which usually would not be able to activate a macroscopic platelet population alone. We thus propose a reaction-rate model which can be translated into a dynamical system for the concentration of activated and naive platelets. This system has similarities to the SIR model, a model for the spread of epidemic diseases, with the di ff erence that multiple contagion carrying agents (here: stimulant secreting activated platelets) need to meet close to each other to spread the contagion (here: activation signal). We will show that in order to achieve macroscopic activation of the platelet population, following minor stimulation - which from a medical point of view can lead to a dangerous thrombus that occludes otherwise undamaged vessels - can only propagate to macroscopic proportions if there is a substantial degree of heterogeneity of platelet sensitivity.

## 2. Model

Our aim is to model the transition from a naive, non-activated population of platelets, which are not able to aggregate, to a situation where a substantial proportion of the platelet population is activated. Here, we only wish to model the first, initiating phase of thrombus formation, namely the collective activation of platelets, which allows platelets to adhere to each other and is a pre-requisite for aggregation of platelets to form a platelet plug.

We make the following simplifying, yet biologically motivated assumptions:

- · The platelet population is well mixed.
- · N aive, non-activated platelets ( N ) can be activated either by external stimuli (e.g. collagen, convulxin), or by paracrine signalling (e.g. ADP)) from other platelets in their immediate vicinity. They require a threshold amount of signalling molecules at their surface (van der Meijden and Heemskerk, 2019; Jongen et al., 2020).
- · Freshly a ctivated platelets ( As ) s ecret signalling molecules which can activate other platelets ('degranulation') (van der Meijden and Heemskerk, 2019).
- · Signalling molecules disperse quickly through di ff usion and advection in sometimes turbulent blood flow (Saqr et al., 2020). Hence, the range of the signalling interaction is very short-ranged, and can be characterised by a stimulant concentration c 0 added by each activated secreting platelet, As , to its immediate environment.
- · The threshold amount of stimulus required to activate a naive platelet, c , might be higher than c 0 , the concentration provided by a single neighbouring secreting platelet.
- · Degranulation occurs for a limited amount of time (Polasek, 2006), after which platelets turn into an idle activated state ( Ai ), that is, they remain activated (and thus able to aggregate and initiate coagulation), but stop secreting activation signals.
- · We consider the scenario in flow or in droplets in which shear is small. Hence, we neglect shear-flow induced activation of platelets, that plays a role for clotting at vessel walls (Rana et al., 2019; Hellmuth et al., 2016).

Figure 1: Illustration of the model for the case c = 3 c 0, meaning m = 3. (Top right:) The propensity to activate as function of the provided stimulant concentration c ′ for c = 3 c 0. Only if c ′ ≥ c , platelets can activate. (Bottom:) Illustration of the cooperative activation of platelets. A single secreting activated platelet provides only a stimulant level c 0, not su ffi cient to activate other naive platelets. However, if 3 secreting platelets are close to a naive platelet, their stimulants add linearly to provide c ′ = 3 c 0 = c , which is su ffi cient to activate the naive platelet.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-64334d73-a639-48fe-8516-e914a626e971.jpg)

From these assumptions, we can formalise the activation dynamics of platelets. Let us denote the concentration of externally provided stimulant as c ex . If m ′ secreting platelets are immediately next to a naive platelet (within the distance over which the stimulant level, c 0 , secreted by a platelet is maintained), they generate a maximal concentration c ′ = c ex + m c ′ 0 in their immediate vicinity. To activate a naive platelet with threshold concentration c , it requires that c ′ &gt; c , and thus,

$$m = \left \lfloor \frac { c - c _ { \text{ex} } } { c _ { 0 } } \right \rfloor \quad \text{(1)}$$

is the minimal number of secreting platelets required to activate that platelet. Hence, we can formally write the events of this model as:

$$N _ { m } + m A _ { s } \to ( m + 1 ) A _ { s }, \ \ A _ { s } \to A _ { i } \quad \quad ( 2 )$$

The first term describes how m activated secreting platelets, As , turn a naive platelet, Nm , into another activated secreting platelet. The second term denotes that a secreting activated platelet stops secreting after some time and becomes an i dle activated platelet Ai . These dynamics are also illustrated in Fig. 1. Wenote that in principle also higher numbers of meeting secreting platelets m ′ &gt; m can activate a platelet, but the probability for this to happen, which decreases exponentially with the number of platelets to meet, is much lower and thus neglected.

As the population is well mixed, we can use the law of mass action to derive a dynamical model for the concentrations of naive platelets, nm , secreting activated platelets, as , and idle ac-

Figure 2: Time course of proportion of activated platelets, a = ae + ai as predicted from model (6) for r = 300 and platelet sensitivities being distributed normally with mean ¯ m = 4 and di ff erent values of standard deviation σ . We have, from left to right, σ = 1 1 2 1 5 2. , . , . , (Top row:) Activated platelets separated for their initial sensitivity, am = ae + ai , with m = 1 (black), m = 2 (red), m = 3 (orange), m = 4 (yellow), m = 5 (green), m = 6 (cyan), m = 7 (blue). (Bottom row:) Total proportion of activated platelets a = ae + ai .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-b960d141-bc0f-4217-b674-4e049ac794e4.jpg)

tivated platelets, ai , from the 'reactions' (2):

the non-dimensional form of Eq. (3):

$$\dot { n } _ { m } & = - \omega n _ { m } a _ { s } ^ { m } & ( 3 ) \\ \dot { a } _ { s } & = \omega n _ { m } a _ { s } ^ { m } - \gamma a _ { s } \\ \dot { a } _ { i } & = \gamma a _ { s } \, & \quad \quad \quad$$

where ˙ : x = dx t ( ) dt for x = nm , as , ai , while ω is the specific activation rate. This rate depends on the range a secreting platelet travels during the time it keeps secreting ('degranulation time' t s (Polasek, 2006)), the signalling range, and the the propensity to activate ˜ ω (as in Fig. 1). The parameter γ = 1 t s is the rate at which a previously activated platelet ceases secreting stimulants.

Notably, for m = 1, this model is equivalent to the SIR model, the paradigmatic model for the spread of a contagion (Kermack, W. O. and McKendrick, 1927). In the context of infectious diseases, the case m &gt; 1 can be seen as a variant of the SIR model in which several infected individuals need to be present at the same time to infect others, i.e. they 'cooperate' to spread the contagion. We therefore refer to this model, in its general form, as a cooperative SIR model . The following analysis will study this model, and its heterogeneous version, for its general properties. Due to its relationship to the SIR model, this may have implications for understanding of infectious diseases, yet here we wish to focus on its implications for collective platelet activation and thus, eventually, for blood clotting.

To simplify the analysis, we will use a non-dimensionalised version of Eqs. (3), in which we use the degranulation time as time unit. Thus, we use rescaled time ˜ t = γ t , and the reproductive number r = ωρ γ , where ρ = n + as + ai is the total concentration of platelets. Furthermore, we express the equations in terms of the proportions of subpopulations, ˜ n = n ρ , ˜ as = as ρ , ˜ ai = ai ρ , so that ˜ n + ˜ as + ˜ ai = 1. We can then eliminate ˜ ai = 1 -˜ nm -˜ as , and express the system by two equations. For convenience, we rename the non-dimensionalised quantities to remove the tilde: ˜ n → n as , ˜ → nm , ˜ ai → ai , ˜ t → t , to arrive at

$$\dot { n } _ { m } = - r n _ { m } a _ { s } ^ { m }$$

$$\dot { a } _ { s } = r n _ { m } a _ { s } ^ { m } - a _ { s } \.$$

We will also study a heterogeneous version of the model. First, we note that heterogeneity in r is not qualitatively different to a homogeneous system. Assuming a distribution of naive platelets with di ff erent r , n r ( ), we have, instead of Eq. (5), ˙ as = R ∞ 0 r n r a ( ) m dr s -as . Since a m s R ∞ 0 r n r dr ( ) = ¯ rnas , the corresponding equation is the same as (5) when replacing r with the mean value ¯. Thus, we do not consider this type of heteror geneity explicitly. Another type of heterogeneity is when the threshold concentration c varies between platelets, according to a probability density distribution n c ( ). Associated with this is are sub-populations nm with di ff erent threshold numbers m = 1 2 , , ... , which are related to n c ( ) by nm = R c ex + ( m + 1) c 0 c ex + mc 0 n c ( ) dc . Then, we get the non-dimensionalised time evolution of the system:

$$\dot { n } _ { m } & = - r n _ { m } a _ { s } ^ { m } \quad \text{ for all } m \in \mathbb { N } \quad \text{ (6)} \\ \dot { a } _ { s } & = r \sum _ { m } n _ { m } a _ { s } ^ { m } - a _ { s }$$

While the sensitivity of a platelet is not relevant for the dynamics after it is activated, we will for our analysis also distinguish activated platelets by sensitivity, defining the proportions of activated platelets with activation threshold m as am : = as ( m ) + a ( m ) i : = nm (0) -nm for m = 1 2 , , ... .

## 3. Results

## 3.1. Collective activation

Our main goal is to study under which circumstances a macroscopic population of platelets activates. First, we consider the situation with diminishing amounts of external stimulant c ex ≈ 0. For a convenient terminology, we define the terms

Figure 3: Equilibrium proportions of activated platelets, separated for sensitivities, am = nm (0) -nm , as function of the width of sensitivity distribution, σ , for time t = 10, as predicted by model (6). Plots are shown for di ff erent values of the distribution's mean ¯ m (rows) and r (columns). Rows have, from top to bottom, ¯ m = 2 3 4 5, columns have, from left to right, , , , r = 5 25 125 625. Colours represent , , , am for di ff erent m as in Fig. 2.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-502ee905-66f0-4f1e-934a-2e0c70995a99.jpg)

'microscopic' 'macroscopic', as follows: / We consider a total of X platelets (or, alternatively, stimulant molecules), which is large, that is, X → ∞ , and a sub-population thereof with Xs platelets. The sub-population Xs is microscopic , if Xs &gt; 0, but its proportion on the total population X , x : = Xs X → 0 for X → ∞ . On the other hand, a population is macroscopic if Xs diverges for X →∞ , in a way that 0 &lt; x ≤ 1.

We assess under which circumstances a macroscopic population of activated platelets, a : = as + ai = 1 -P m nm , emerges, if initially only a microscopic population of the size as = ϵ → 0 of secreting activated platelets is present, if external stimulant is absent or only microscopic, c ex = 0. This is equivalent to asking under which circumstances an epidemic breaks out in a contagion model such as the SIR model. As for m = 1, the model (4),(5) is equivalent to the SIR model, it is well known that the condition for an epidemic to break out is for r &gt; 1, if initially all individuals, except for a microscopic proportion, are susceptible - which in our case corresponds to all but a few platelets being naive.

We now wish to study this question for general m &gt; 1 in Eqs. (4),(5) and then for the heterogeneous model, Eq. (6). We first note that the condition for a macroscopic platelet population to become activated upon exposure to a microscopic population of secreting activated platelets, as ( t = 0) = ϵ → 0, is equivalent to the fixed point x ∗ = ( n = 1 , as = 0) being unstable. Thus, to assess this, we have a look at the Jacobian matrix of Eqs. (4),(5)

at that fixed point for m &gt; 1,

$$\mathfrak { h } _ { 0 }, \quad J | _ { x ^ { * } } = \begin{pmatrix} - r a _ { s } ^ { m } & - r m n a _ { s } ^ { m - 1 } \\ r a _ { s } ^ { m } & r m n a _ { s } ^ { m - 1 } - 1 \end{pmatrix} | _ { ( n = 1, a _ { s } = 0 ) } = \begin{pmatrix} 0 & 0 \\ 0 & - 1 \end{pmatrix} \. \quad ( 7 )$$

This matrix has eigenvalues 0 and -1, therefore the fixed point is Lyapunov stable. This means that for m &gt; 1, no macroscopic activation of the platelet population can occur when seeded with an infinitesimally small population of activated platelets, as ( t = 0) = ϵ . This has some advantages from the biological point of view: since macroscopic activation of platelets can lead to blood clotting, and inadvertent clotting can lead to dangerous pathologies through thrombosis, it could be a dangerous situation if a tiny population of such platelets were su ffi cient to trigger this. Hence, a situation with m &gt; 1 is protecting against inadvertent clotting.

However, it has been shown that platelet sensitivity is highly heterogeneous (Jongen et al., 2020), therefore, we study under which circumstances the heterogeneous model, Eq. (6), with a distribution of platelet thresholds nm ( t = 0) = nm (0) , can lead to macroscopic activation. Taking the Jacobian matrix of the heterogeneous model, Eq. (6), gives:

$$\intertext { o } \intertext { t o } J = \begin{pmatrix} - r a _ { s } & 0 & 0 & \cdots & - r n _ { 1 } \\ 0 & - r a _ { s } ^ { 2 } & 0 & \cdots & - 2 r n _ { 2 } a _ { s } \\ \vdots & \vdots & \cdots & \vdots & \vdots \\ r a _ { s } & r a _ { s } ^ { 2 } & \cdots & \cdots & r \sum _ { m } m n _ { m } a _ { s } ^ { m - 1 } - 1 \end{pmatrix} \,, \quad ( 8 )$$

Figure 4: Equilibrium values of total proportion of activated platelets, a = ae + ai , as function of the width of (Normal) sensitivity distribution, σ , for nondimensionalised time t = 10, as predicted by model (6). Plots are shown for di ff erent values of the distribution's mean ¯ m (rows) and signalling range r (columns). Rows have, from top to bottom, ¯ m = 2 3 4 5, columns have, from left to right, , , , r = 5 25 125 625. , , ,

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-34cf2bfa-6d6d-4f18-99a1-23e0aa95c7e8.jpg)

where the size of the matrix is m max + 1 with m max being the largest m for which nm (0) &gt; 0. Taking the Jacobian at the fixed point nm = nm (0) for m = 1 2 , , ... and as (0) = 0, we get,

$$J = \begin{pmatrix} 0 & 0 & 0 & \cdots & - r n _ { 1 } ^ { ( 0 ) } \\ 0 & 0 & 0 & \cdots & - 2 r n _ { 2 a _ { s } } \\ \vdots & \vdots & \cdots & \vdots \\ 0 & 0 & \cdots & \cdots & r n _ { 1 } ^ { ( 0 ) } - 1 \end{pmatrix} \,. \quad \text{ (9)} \quad \text{tion} \\ \quad \text{time} \\ m, w$$

This is a triangular matrix, which thus has the eigenvalues 0 (with multiplicity m max) and r n (0) 1 -1. As the latter eigenvalue is positive for r n (0) 1 &gt; 1, we get that for those values the fixed point is unstable and thus a macroscopic proportion of the platelet population becomes activated, while for r n (0) 1 &lt; 1, no macroscopic activation occurs, i.e. a ≈ 0. For notational convenience, we thus call n ∗ 1 : = 1 r and r ∗ : = 1 n (0) 1 , respectively, the critical values which, when exceeded, lead to macroscopic activation. This is not surprising, since the sub-population of hyper-sensitive platelets, with m = 1, follows the dynamics of the SIR model. In that case (when rn (0) 1 &gt; 1), we have that ˙ as &gt; rn 1 -1 &gt; 0 as long as n 1 &gt; 1 r . Since the activation wave ceases only once ˙ a &lt; 0, this means that eventually n 1 &lt; 1 r , and thus the final proportion of activated platelets, for t = t f →∞ , has a lower bound: a = n (0) 1 -n 1 &gt; n (0) 1 -1 r .

whether platelets with less sensitivity, that is, higher thresholds of m , will be activated to macroscopic proportions. To assess this, we consider numerical solutions of Eqs. (6) for di ff erent distributions of platelet sensitivities and values of r . We first consider the scenario where platelet sensitivities are distributed according to a Normal distribution with mean activation threshold ¯ m and standard deviation σ . In Fig. 2 we see time courses of activation, separated for platelets with di ff erent m , with am = nm (0) -nm for fixed ¯ m = 4 and di ff erent values of σ . Notably, we see that for low values of σ , no activation occurs, for larger values, activation occurs, but only for population n 1 , while for even larger values activation of populations nm with higher m occurs as well. Thus, we observe that in a heterogeneous population, the population with the highest sensitivity m = 1, but also populations with the lower sensitivity (larger m &gt; 1) are activated if signalling strength r and heterogeneity, characterised by σ , are su ffi ciently high.

At first glance, it is only obvious that the sub-population of platelets with m = 1 is activated, while it is not clear

To further assess how the mean sensitivity and the heterogeneity of the sensitivity a ff ect the collective activation of platelets, we show in Figs. 3 and 4 the end points of the time courses of Figure 2 (as the curves reach a plateau) as a function of the standard deviation of the distribution, σ . Total activated proportions a (Fig. 4) and individual proportions am (Fig. 3, as curves of di ff erent colours) are shown for various values of ¯ m and r . Wesee that for fixed values of ¯ , but small heterogeneity m σ , there is no macroscopic activation, but as the heterogeneity increases, macrosopic activation with a &gt; 0 occurs. For small

Figure 5: Phase portrait projected on as , n 1 (collection of trajectories as ( n 1)), rescaled by n ∗ 1 = 1 r for small as when O a s ( 3 ) can be neglected and nm ≈ const for m &gt; 1. Parameter r = 100. The emphasised blue trajectory shows the trajectory converging to the fixed point ( n 1 = 0 , as = 1 n 2 r ). This trajectory denotes the stable manifold which separates basins of attraction.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-b5f79a90-48f2-47ad-b0fd-77bc6c523d4e.jpg)

values of r , we see that when increasing σ , first only the population with m = 1 is activated and then, for higher σ populations with lower sensitivity ( m &gt; 1) are activated. Notably, for higher r , the activated proportions of the populations with m &gt; 1 exceed substantially the activated population with m = 1. While we know that lower sensitive populations with m &gt; 1 cannot be activated macroscopically by microscopic as , it appears that the population with m = 1 is initially activated, generating a nonzero proportion as , which serves as a 'seed' to activate lower sensitive populations. This activation may then sustain itself if as is large enough.

In order to understand this behaviour, we study analytically under which circumstances the population of activated platelets remains small, i.e. when a ≪ 1 and thus as ≪ 1. We wish to study this for arbitrary distributions of platelet sensitivities nm ( t = 0) = nm (0) , m = 1 2 , , ... , but focus on situations where the population of hyper-sensitive platelets with m = 1 is very small, while the other populations are much larger: n (0) 1 ≪ 1 and nm (0) ≫ n (0) 1 are such that n (0) 1 nm (0) ∼ O as ( ) ≪ 1 for all m &gt; 1. Furthermore, since macroscopic activation can only occur for n (0) 1 &gt; 1 r we additionally assume that n ∗ 1 = 1 r ≪ 1, that is, r is large. This is also a biologically realistic regime: as the degranulation time is around 60s (Polasek, 2006), which is roughly the time blood circulates once through the entire circulatory system, any secreting platelet can potentially get in contact with a very large number of naive platelets and activate them.

Since, according to these assumptions, am ≪ 1, we can assume nm ≈ nm (0) being approximately constant for m &gt; 1, and it is therefore su ffi cient to consider only n 1 , as as dynamical variables. We note that the first two terms of Eq. (6), rn as 1 and rn a s 2 2 are both of order O a s ( 2 ) while the other terms are of order O a s ( 3 ), thus we have, in leading order of (small) as for Eq. (6):

$$\dot { n } _ { 1 } & = - r n _ { 1 } a _ { s } & ( 1 0 ) \\ \dot { a } _ { s } & = r n _ { 1 } a _ { s } + r n _ { 2 } a _ { s } ^ { 2 } - a _ { s } + O ( a _ { s } ^ { 3 } ) \approx r n _ { 1 } a _ { s } + r n _ { 2 } a _ { s } ^ { 2 } - a _ { s }$$

We can now study the fixed points of this system when neglecting O a s ( 3 ). We see that there are fixed points for any as = 0, and there is one fixed point for as &gt; 0, namely for n 1 = 0 and as = 1 n r 2 = : a s ∗ . A linear stability analysis (see Appendix A) shows that this is a saddle point, which means that there is a stable manifold, which, since this is a two-dimensional system, consists of the trajectories that converge to the point x ∗ = ( n 1 = 0 , as = 1 n r 2 ). This stable manifold separates the phase space into two basins of attraction, one which converges to as = 0 and one where trajectories diverge. To illustrate this, we show the phase portrait of system (10) for r = 100 in Fig. 6. Here, we see that the highlighted trajectory separates the other trajectories according to their fate. If n (0) 1 lies below this curve, as cannot exceed a s ∗ = 1 n r 2 . On the other hand, if it lies above that curve, as will diverge, meaning that our approximation will break down and thus as ≫ 1 n r 2 .

To find the stable manifold which separates basins of attraction, we express trajectories as functions as ( n 1 ) in the n 1 -as -plane. They can be found as solutions to the di ff erential equation,

$$\frac { d a _ { s } } { d n _ { 1 } } = \frac { \frac { d a _ { s } } { d t } } { \frac { d n _ { 1 } } { d t } } \approx - 1 - \frac { n _ { 2 } } { n _ { 1 } } a _ { s } ( n _ { 1 } ) + \frac { 1 } { r n _ { 1 } } \. \quad \ \ ( 1 1 )$$

This is a linear ODE, whose solution can be found, for example, through the integrating factor method. As we are looking for a trajectory converging to ( n 1 = 0 , as = 1 n r 2 ), we can choose the initial condition as ( n 1 = 0) = 1 n r 2 , and with this we get the solution,

$$a _ { s } ( n _ { 1 } ) = \frac { 1 + n _ { 2 } - r n _ { 1 } n _ { 2 } } { r n _ { 2 } ( 1 + n _ { 2 } ) } \. \quad \quad ( 1 2 )$$

Furthermore, since ˙1 n &lt; 0, the direction of this trajectory is towards x ∗ , hence this is indeed the stable manifold. Now, we need to distinguish whether n (0) 1 is below or above that trajectory. To this end, we determine the intersection of this trajectory with the n 1 axis, which we find as the solution to 0 = 1 + n 2 -rn n 1 2 rn 2(1 + n 2) , which is n ∗∗ 1 = 1 + n 2 rn 2 . Hence, for n (0) 1 &lt; n ∗∗ 1 , as , and thus a , remain below a s ∗ = 1 n r 2 , while for n (0) 1 &gt; n ∗∗ 1 , as , and thus a , exceed this value substantially. Note that in the case r → ∞ , 1 n r 2 = a s ∗ must be microscopic, so the condition for a remaining microscopic is n (0) 1 &lt; 1 + n 2 rn 2 . In that case, we thus get that substantial proportions of am emerge once n (0) 1 exceeds n ∗∗ 1 = n ∗ 1 1 + n 2 n 2 . This is shown in Fig. 6 to be an accurate prediction: there final activation a is shown as function of n 1 for di ff erent values of n 2 and a mean value ¯ m = 3 1 .

1 Symmetry of the sensitivity distribution then implies n (0) 3 = 1 -2( n (0) 1 + n (0) 2 ) , n (0) 4 = n (0) 2 , n (0) 5 = n (0) 1

Figure 6: Equilibrium proportions of activated platelets, a = as + ai , as function of n (0) 1 , rescaled by n ∗ 1 = 1 / r . The plots show di ff erent values of n 2 in di ff erent columns, from left to right n 2 = 0 0 2 0 3 0 4. We assume a symmetric distribution with , . , . , . ¯ m = 3, which further implies n (0) 3 = 1 -2 n (0) 1 -2 n (0) 2 , n (0) 4 = n (0) 2 , n (0) 5 = n (0) 1 . The vertical dashed line indicates the critical value n ∗∗ 1 = n ∗ 1 1 + n 2 n 2 . (Top row:) Illustration of the distributions nm (0) for m = 1 2 , , ..., 5 used below. (Middle row:) Total activated proportion a . (Bottom row:) shows the activated populations separated by m am , = nm (0) -nm .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-5d6d3921-6a03-4144-b1c2-fcf3fe26fdbd.jpg)

This demonstrates that a microscopic proportion of hypersensitive platelets, nm (0) , is su ffi cient to trigger a macroscopic activation also of platelet populations with less sensitivity, m &gt; 1, including those with average sensitivity ¯ . m Nonetheless, for n (0) 2 = 0, no substantial activation is observed. Thus, we conclude that the presence of platelets with intermediate sensitivity, between hyper-sensitivity with m = 1, and the mean sensitivity ¯ , are essential for macroscopic activation. Although we do m not have analytical means to study further terms in Eq. (6), we expect that more populations with m &gt; 2 can further amplify the spread of the activation signal if nm (0) ≫ n (0) 1 , as has been indicated above in Figs. 2-4.

## 3.2. Comparison with experiments with external stimulants

We now consider the situation where an external stimulant is present which is able to directly activate platelets, c ex &gt; 0. According to Eq. (1), the threshold number of secreting platelets m required to activate platelets with threshold concentration c is m = ⌊ c -c ex c 0 ⌋ . Hence, given a probability density distribution n c ( ) of platelets with threshold concentrations c , we get a distribution of initially naive platelets of nm = n GLYPH&lt;22&gt; c + c ex c 0 GLYPH&lt;23&gt; . This means that with increasing c ex the distribution of platelets gets shifted towards lower m . Importantly, platelets with c &lt; c ex are activated from the beginning, when the stimulant is given. We thus have as initial condition as ( t = 0) = R c ′ &lt; c ex n c ( ′ ) dc ′ ≈ P m &lt; c ex c 0 nm . Hence, with increasing c ex , both the proportion of initially secreting activated platelets increases, and the distribution of platelets is shifted towards lower m , consequently enhancing activation.

We wish to compare our model results with experimental results from droplet microfluidics assays with platelets, performed by Jongen et al. (2020). In one of those experiments, platelets were enclosed in micro-droplets and exposed to the external stimulant convulxin, at varying concentration, and then the abundance of the activation marker P-selectin was measured via fluorescent cytometry, that is, the fluorescent intensity, I , was recorded (Fig. 2B in (Jongen et al., 2020)). Two scenarios are compared: (A) where each droplet contains a single platelet, (B) where droplets contain platelet collectives. In situation (A) no platelet can activate another platelet through paracrine signalling since they are separated by droplets, hence c = c ex and thus the proportion of activated platelets at a given concentration of convulxin is the cumulative distribution F c ( ) = as ( t = 0) = R c ′ &lt; c n c ( ′ ) dc ′ . Hence, the distribution of sensitivities can directly be inferred from this data. In scenario (B), platelets can activate each other, and thus this measures the collective activation of a platelet population. This can be compared with the predicted proportion of activated platelets, a , from our model.

However, our model has limitations modelling this data, as the experimental setup does not meet all assumptions of our model. In particular, since platelet collectives are encapsulated in a flow-free confined environment (droplets), the agonists secreted by platelets are not quickly dispersed, and thus can accumulate and reach platelets which are even further away. This means that we can expect an enhanced paracrine activation between platelets, resulting in a higher value of r than would be expected in an in vivo situation. Nonetheless, we wish to test whether the qualitative behaviour, namely the predicted

collective activation of the population through a small, hypersensitive population of platelets, prevails. To that end, rather than finding accurate precise numerical model fits, we explore the parameter space to find values for which the model qualitatively reproduces the data.

To determine an estimate for the distribution n c ( ) and thus nm , we fit a cumulative distribution function to the data of single-platelet droplets (A). First, we need to renormalise the experimental fluorescence measurements, I , from Jongen et al. (2020), as there is always some background signal even in absence of activation. Hence, we define the relative activation as a exp = I -I min I max -I min where I min and I max are minimum and maximum fluorescence values. As the data in (Jongen et al., 2020) shows a symmetric sigmoidal form on a logarithmic scale of c , we will attempt to fit a log-normal distribution, with mean ¯ and logac rithmic variance σ 2 c to the data. For an easier comparison with our model, which uses m = c c 0 to define platelet populations, we will try only values ¯ c = ¯ mc 0 and σ c = σ mc 0 with integer ¯ m , σ m = 1 2 , , ... . To that end, we first need to find an estimate for c 0 . As for c ex = c 0 the sub-population with m = 1 is being activated for the first time upon increasing values of c ex , we identify c 0 as the concentration for which, upon increasing concentrations of c ex , for the first time substantial collective activation emerges in data (B). While due to the lack of intermediate data points, this estimate can only be a rough one, we can estimate this from inspection of data (B) to be around c 0 = 0 3 . ng mL. While some activation is present already at / c ex = 0 1 . ng mL, this could be due to spontaneous activation that is ex-/ pected from our model according to our previous discussion.

We find that for ¯ m = 14 and σ m = 2, corresponding to ¯ c = 4 2 ng mL and . / σ m = 0 6 ng mL, a good match of the . / log-normal distribution with the single-platelet droplet data, (A), is achieved, as is shown in Fig. 7A. We will thus use -(ln m -ln ¯ ) m 2

nm (0) = N e 2 σ 2 m as the initial distribution of platelet activation thresholds in absence of external stimulant, where N is a normalisation factor that accounts for the discreteness of the values m . Then, we determine the initial distribution of platelet sensitivities with external stimulant c ex as nm ( t = 0) = n (0) m c + ex / c 0 and as ( t = 0) = P m &lt; c ex c 0 nm (0) , and we solve the system (6) numerically using Scipy's solve\_ivp function. The result is shown in Fig. 7B for di ff erent values of r , together with the experimental data. We note that for smaller values of r no reasonable match is achieved, but for very high values of r , in particular for r = 30000 (green curve in Fig. 7B), a good qualitative agreement is achieved. This is consistent with the expectation that due to the accumulation of paracrine signalling molecules within droplets, a single platelet can reach and activate a much higher population of platelets than would be in a (possibly turbulent) in vivo flow environment. Overall, this demonstrates that our model is able to qualitatively reproduce the experimental data on collective platelet activation, though at parameter values which are likely beyond the corresponding in vivo values.

Figure 7: Comparison with experimental results from (Jongen et al., 2020). Pluses show experimental results measuring overall P-selectin activity (a marker for platelet activation) as fluorescence levels via flow cytometry, rescaled for background fluorescence and as proportion of maximal fluorescence (see text). Error bars are smaller than the symbol sizes. Curves show model predictions for a = ae + ai . The horizontal axis shows varying concentrations of stimulant (convulxin) to which platelets were exposed in the experiments, the y-axis shows relative activation. (A) Shows a log-normal distribution with mean c = 4 2 ng mL and standard deviation . / σ = 0 6 ng mL, together with . / experimental results where single platelets are encapsulated in single droplets. This is a proxy for the cumulative distribution F c ( ) = R c ′ &lt; c n c ( ′ ) dc ′ . (B) Shows model output for a , compared with experimental results where many platelets are in single droplets, for di ff erent values of r : r = 1000 (black), r = 3000 (red), r = 10000 (orange), r = 30000 (green), r = 100000 (cyan).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-513a1573-7710-4121-92d5-1b6d6193727e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00039v1-figure-9876edfc-2626-4677-8773-8183b01cdf81.jpg)

## 4. Conclusions

Blood clotting is an important physiological process, but when it occurs inadvertently, it can cause thrombosis which can be live threatening. Thus, understanding and eventually predicting blood clotting can aid in preventing thrombosis, while retaining the physiological features of blood clotting. The initiating stage on the path towards blood clotting is the collective activation of platelets, hence understanding the platelet activation dynamics will aid in assessing blood clotting risk.

We have introduced a minimal model for the activation dynamics of blood platelets, which includes activation by external stimuli and through stimulation via paracrine, platelet-platelet signalling. In particular, we considered the scenario in blood flow, where due to fast dispersion of stimulants, one or more activated platelets that secret activation signals are required in the immediate vicinity of a non-activated, naive platelet, to activate the latter. From these and other biologically motivated assumptions, we formulated a minimal dynamical model, (3), for the time evolution of naive and activated platelets, with the aim to predict the proportion of activated platelets when starting with a predominantly naive platelet population. If the number of secreting platelets required to activate a naive one, m , is equal to one, this model is equivalent to the SIR model the paradigmatic contagion model - and macroscopic activation occurs whenever the 'reproductive number' r is above a critical value, even for microscopic proportions of initially activated secreting platelets. For m &gt; 1 the model, which we then call cooperative SIR model , still has the characteristic of a contagion model, yet, if the population is uniform, no macroscopic activation, of a non-zero proportion of a large platelet population, is possible, if initially only few platelets are secreting.

Nonetheless, it is known that platelet sensitivity varies considerably between platelets, therefore we considered a heterogeneous version of the model where m may di ff er between platelets. We showed that in this case, a macroscopic activation of the population is possible when only few activated platelets are initially present, if r is su ffi ciently high. In that case, the hyper-sensitive sub-population with m = 1 becomes activated, and this activated, secreting sub-population then serves as a 'seed' to activate sub-populations with higher m , in a staggered cascade of activating ever less sensitive sub-populations (with higher m &gt; 1). Notably, we showed analytically and numerically that a microscopic proportion of platelets with m = 1 is su ffi cient to achieve this, if r is su ffi ciently large. This supports the hypothesis that a rare population of hyper-sensitive platelets could be able to mediate macroscopic activation and thus trigger blood clotting (Baaten et al., 2017; Jongen et al., 2020; Lesyk and Jurasz, 2019). We further showed that in order to activate the bulk of the platelet population, also platelets with less, yet higher-than-average sensitivity, with thresholds m between that of hyper-sensitive and average platelets (thresholds 1 &lt; m &lt; ¯ ), m are required to achieve macroscopic activation of the bulk platelet population. This is biologically required to then achieve macroscopic platelet aggregation and eventually blood clotting. A 'gap' in sensitivities, for example if there are no, or only few platelets with m = 2, may break the activation wave, and no substantial activation of the platelet population as a whole emerges.

Finally, we compared our model results with experimental results where platelet activation and aggregation was measured with platelets encapsulated in micro droplets (Jongen et al., 2020). While this experimental setup did not meet all our model assumptions - in particular since secretory, stimulating signals can accumulate - we showed that our model is capable of qualitatively reproducing the observed dose-response curves, namely the activation intensity as function of provided external stimulant, if the reproductive number r is very large.

As most clinical platelet function tests for assessing thrombosis risk measure bulk platelet properties to assess thrombosis risk (Anghel et al., 2020; Pabinger and Ay, 2009), these findings suggests that clinical practice may need to be revised: if the distribution of (higher-than average) platelet sensitivities is determining the onset of blood clotting, rather than bulk platelet sensitivity, then measurements of the latter - for example via droplet microfluidics as done by Jongen et al. (2020) - would be required for appropriate assessment of the thrombosis risk.

To summarise, we showed that a simple model, reminiscent of the SIR model, is able to reproduce measured qualitative features of collective platelet activation. This model explains mechanistically how even very rare hyper-sensitive platelets can serve as an activation seed that may propagate the activation through the whole population to achieve a macroscopic activation, if su ffi cient amounts of less, but larger-than-average sensitivities are present. Once a substantial fraction of the platelet population have been activated, they may then aggregate and form the body of a blood clot. This model serves as a first step towards a development of a comprehensive model for whole blood clotting, that may allow quantitative predictions for blood clotting probability, and hence the risk of thrombosis.

## Acknowledgements

I thank Jonathan West for guidance in processing the experimental data of (Jongen et al., 2020).

## Appendix A.

We wish to analyse the stability and stable manifold of the fixed point of Eq. (10). This system has fixed points for any as = 0 and for as = 1 n r 2 and n 1 = 0. The Jacobian of the system at this fixed point is:

$$\underset { \bullet \cdot } { \text{if} } \underset { \bullet } { \text{if} } J = \begin{pmatrix} - r a _ { s } & - r n _ { 1 } \\ r a _ { s } & r n _ { 1 } + 2 r n _ { 2 } a _ { s } - 1 \end{pmatrix} | _ { n _ { 1 } = 0, a _ { s } = 1 / r } = \begin{pmatrix} - 1 & 0 \\ 1 & 1 \end{pmatrix} \quad ( \text{A.1} )$$

This matrix is in triangular form and thus has eigenvalues -1 and 1. Hence it is a saddle point, which possesses a stable trajectory, which separates the space into to basins of attraction, defined as the trajectory coverging to the fixed point. In the main text, we show that as ( n 1 ) = 1 + n 2 -rn n 1 2 rn 2(1 + n 2) is this trajectory.

## References

Anghel, L., Sasc˘ au, R., Radu, R., St˘ atescu, C., 2020. From classical laboratory parameters to novel biomarkers for the diagnosis of venous thrombosis. International Journal of Molecular Sciences 21. doi: 10.3390/ ijms21061920 .

Baaten, C.C., ten Cate, H., van der Meijden, P.E., Heemskerk, J.W., 2017. Platelet populations and priming in hematological diseases. Blood Reviews 31, 389-399. doi: 10.1016/j.blre.2017.07.004 .

British Heart Foundation, 2021. Global Heart and Circulatory Diseases Factsheet.

Fogelson, A.L., 1993. Continuum Models of Platelet Aggregation: Mechanical Properties and Chemically-Induced Phasen Transition. Contemporary Mathematics 141, 279.

Fogelson, A.L., 2016. Continuum Models of Platelet Aggregation : Formulation and Mechanical Properties. SIAM Journal on Applied Mathematics 52, 1089-1110.

Fogelson, A.L., Neeves, K.B., 2015. Fluid Mechanics of Blood Clot Formation. Annu. Rev. Fluid. Mech. 47, 377. doi: 10.1146/ annurev-fluid-010814-014513.Fluid .

Gale, A.J., 2011. Current Understanding of Hemostasis. Toxicologic Pathology 39, 273-280. doi: 10.1177/0192623310389474 .

Heemskerk, J.W., West, J., 2022. Emerging technologies for understanding platelet diversity,. Arteriosclerosis, Thrombosis, and Vascular Biology 42, 540-552.

Hellmuth, R., Bruzzi, M.S., Quinlan, N.J., 2016. Analysis of Shear-Induced Platelet Aggregation and Breakup. Annals of Biomedical Engineering 44, 914-928. doi: 10.1007/s10439-015-1409-1 .

Jongen, M.S., MacArthur, B.D., Englyst, N.A., West, J., 2020. Single platelet variability governs population sensitivity and initiates intrinsic heterotypic responses. Communications Biology 3, 281. URL: http://dx.doi.org/ 10.1038/s42003-020-1002-5 , doi: 10.1038/s42003-020-1002-5 .

Kermack, W. O., McKendrick, G.A., 1927. A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society A 115, 700. doi: 10. 21136/cpmf.1924.109357 .

LaPelusa, A., Dave, H.D., 2023. Physiology, Hemostasis.

Lesyk, G., Jurasz, P., 2019. Advances in Platelet Subpopulation Research. Frontiers in Cardiovascular Medicine 6, 138. doi: 10.3389/fcvm.2019. 00138 .

Link, K.G., Sorrells, M.G., Danes, N.A., Neeves, K.B., Leiderman, K., Fogelson, A.L., 2020. A mathematical model of platelet aggregation in an extravascular injury under flow. Multiscale Modeling and Simulation 18, 1489-1524. doi: 10.1137/20M1317785 arXiv:2003.02251 , .

van der Meijden, P.E.J., Heemskerk, J.W.M., 2019. Platelet biology and functions: new concepts and clinical perspectives. Nature Reviews Cardiology 16, 166-179. doi: 10.1038/s41569-018-0110-0 .

Pabinger, I., Ay, C., 2009. Biomarkers and venous thromboembolism. Arteriosclerosis, Thrombosis, and Vascular Biology 29, 332-336. doi: 10.1161/ ATVBAHA.108.182188 .

Polasek, J., 2006. Three Modes of Platelet Degranulation. Pathophysiology of Haemostasis and Thrombosis 35, 408-409. doi: 10.1159/000097698 .

Rana, A., Westein, E., Niego, B., Hagemeyer, C.E., 2019. Shear-Dependent Platelet Aggregation: Mechanisms and Therapeutic Opportunities. Frontiers in Cardiovascular Medicine 6, 1-21. doi: 10.3389/fcvm.2019.00141 .

Raskob, G., Angchaisuksiri, P., Blanco, A., Buller, H., Gallus, A., Hunt, B., Hylek, E., Kakkar, A., Konstantinides, S., McCumber, M., Ozaki, Y., Wendelboe, A., Weitz, J.I., 2014. Thrombosis - A Major Contributor to Global Disease Burden. Arterioscler. Thromb. Vasc. Biol. 34, 2363-2371. doi: 10.1161/ATVBAHA.114.304488 .

Saqr, K.M., Tupin, S., Rashad, S., Endo, T., Niizuma, K., Tominaga, T., Ohta, M., 2020. Physiologic blood flow is turbulent. Scientific Reports 10, 1-12. doi: 10.1038/s41598-020-72309-8 .

- Schoeman, R.M., Rana, K., Danes, N., Lehmann, M., Di Paola, J.A., Fogelson, A.L., Leiderman, K., Neeves, K.B., 2017. A Microfluidic Model of Hemostasis Sensitive to Platelet Function and Coagulation. Cellular and Molecular Bioengineering 10, 3-15. doi: 10.1007/s12195-016-0469-0 .

Stalker, T.J., Newman, D.K., Ma, P., Wannemacher, K.M., Brass, L.F., 2012. Handbook of Experimental Pharmacology 210, 59-85.

Platelet signalling. doi: 10.1007/978-3-642-29423-5 .

Wang, N.T., Fogelson, A.L., 1999. Computational Methods for Continuum Models of Platelet Aggregation. Journal of Computational Physics 151, 649-675. doi: 10.1006/jcph.1999.6212 .
Philip Greulich
Cooperative SIR dynamics as a model for spontaneous blood clot initiation
Blood clotting is an important physiological process to suppress bleeding
upon injury, but when it occurs inadvertently, it can cause thrombosis, which
can lead to life threatening conditions. Hence, understanding the microscopic
mechanistic factors for inadvertent, spontaneous blood clotting, in absence of
a vessel breach, can help in predicting and adverting such conditions. Here, we
present a minimal model -- reminiscent of the SIR model -- for the initiating
stage of spontaneous blood clotting, the collective activation of blood
platelets. This model predicts that in the presence of very small initial
activation signals, macroscopic activation of the platelet population requires
a sufficient degree of heterogeneity of platelet sensitivity. To propagate the
activation signal and achieve collective activation of the bulk platelet
population, it requires the presence of, possibly only few, hyper-sensitive
platelets, but also a sufficient proportion of platelets with intermediate, yet
higher-than-average sensitivity. A comparison with experimental results
demonstrates a qualitative agreement for high platelet signalling activity.