Platelet Biorheology and Mechanobiology in Thrombosis and Hemostasis: Perspectives from Multiscale Computation.

Thrombosis is the pathological clot formation under abnormal hemodynamic conditions, which can result in vascular obstruction, causing ischemic strokes and myocardial infarction. Thrombus growth under moderate to low shear (<1000 s-1) relies on platelet activation and coagulation. Thrombosis at elevated high shear rates (>10,000 s-1) is predominantly driven by unactivated platelet binding and aggregating mediated by von Willebrand factor (VWF), while platelet activation and coagulation are secondary in supporting and reinforcing the thrombus. Given the molecular and cellular level information it can access, multiscale computational modeling informed by biology can provide new pathophysiological mechanisms that are otherwise not accessible experimentally, holding promise for novel first-principle-based therapeutics. In this review, we summarize the key aspects of platelet biorheology and mechanobiology, focusing on the molecular and cellular scale events and how they build up to thrombosis through platelet adhesion and aggregation in the presence or absence of platelet activation. In particular, we highlight recent advancements in multiscale modeling of platelet biorheology and mechanobiology and how they can lead to the better prediction and quantification of thrombus formation, exemplifying the exciting paradigm of digital medicine.

SIPA in 10 milliseconds: VWF tentacles agglomerate and capture platelets under high shear.

Shear-induced platelet aggregation (SIPA) occurs under elevated shear rates (10 000 s-1) found in stenotic coronary and carotid arteries. The pathologically high shear environment can lead to occlusive thrombosis by SIPA from the interaction of nonactivated platelets and von Willebrand factor (VWF) via glycoprotein Ib-A1 binding. This process under high shear rates is difficult to visualize experimentally with concurrent molecular- and cellular-resolutions. To understand this fast bonding, we employ a validated multiscale in silico model incorporating measured molecular kinetics and a thrombosis-on-a-chip device to delineate the flow-mediated biophysics of VWF and platelets assembly into mural microthrombi. We show that SIPA begins with VWF elongation, followed by agglomeration of platelets in the flow by soluble VWF entanglement before mural capture of the agglomerate by immobilized VWF. The entire SIPA process occurs on the order of 10 milliseconds with the agglomerate traveling a lag distance of a few hundred microns before capture, matching in vitro results. Increasing soluble VWF concentration by ∼20 times in silico leads to a ∼2 to 3 times increase in SIPA rates, matching the increase in occlusion rates found in vitro. The morphology of mural aggregates is primarily controlled by VWF molecular weight (length), where normal-length VWF leads to cluster or elongated aggregates and ultra-long VWF leads to loose aggregates seen by others' experiments. Finally, we present phase diagrams of SIPA, which provides biomechanistic rationales for a variety of thrombotic and hemostatic events in terms of platelet agglomeration and capture.

How the spleen reshapes and retains young and old red blood cells: A computational investigation.

The spleen, the largest secondary lymphoid organ in humans, not only fulfils a broad range of immune functions, but also plays an important role in red blood cell's (RBC) life cycle. Although much progress has been made to elucidate the critical biological processes involved in the maturation of young RBCs (reticulocytes) as well as removal of senescent RBCs in the spleen, the underlying mechanisms driving these processes are still obscure. Herein, we perform a computational study to simulate the passage of RBCs through interendothelial slits (IES) in the spleen at different stages of their lifespan and investigate the role of the spleen in facilitating the maturation of reticulocytes and in clearing the senescent RBCs. Our simulations reveal that at the beginning of the RBC life cycle, intracellular non-deformable particles in reticulocytes can be biomechanically expelled from the cell upon passage through IES, an insightful explanation of why this peculiar "pitting" process is spleen-specific. Our results also show that immature RBCs shed surface area by releasing vesicles after crossing IES and progressively acquire the biconcave shape of mature RBCs. These findings likely explain why RBCs from splenectomized patients are significantly larger than those from nonsplenectomized subjects. Finally, we show that at the end of their life span, senescent RBCs are not only retained by IES due to reduced deformability but also become susceptible to mechanical lysis under shear stress. This finding supports the recent hypothesis that transformation into a hemolyzed ghost is a prerequisite for phagocytosis of senescent RBCs. Altogether, our computational investigation illustrates critical biological processes in the spleen that cannot be observed in vivo or in vitro and offer insights into the role of the spleen in the RBC physiology.

Computational modeling of biomechanics and biorheology of heated red blood cells.

Because of their compromised deformability, heat denatured erythrocytes have been used as labeled probes to visualize spleen tissue or to assess the ability of the spleen to retain stiff red blood cells (RBCs) for over three decades, e.g., see Looareesuwan et al. N. Engl. J. Med. (1987). Despite their good accessibility, it is still an open question how heated RBCs compare to certain diseased RBCs in terms of their biomechanical and biorheological responses, which may undermine their effective usage and even lead to misleading experimental observations. To help answering this question, we perform a systematic computational study of the hemorheological properties of heated RBCs with several physiologically relevant static and hemodynamic settings, including optical-tweezers test, relaxation of prestretched RBCs, RBC traversal through a capillary-like channel and a spleen-like slit, and a viscometric rheology test. We show that our in silico RBC models agree well with existing experiments. Moreover, under static tests, heated RBCs exhibit deformability deterioration comparable to certain disease-impaired RBCs such as those in malaria. For RBC traversal under confinement (through microchannel or slit), heated RBCs show prolonged transit time or retention depending on the level of confinement and heating procedure, suggesting that carefully heat-treated RBCs may be useful for studying splenic- or vaso-occlusion in vascular pathologies. For the rheology test, we expand the existing bulk viscosity data of heated RBCs to a wider range of shear rates (1-1000 s-1) to represent most pathophysiological conditions in macro- or microcirculation. Although heated RBC suspension shows elevated viscosity comparable to certain diseased RBC suspensions under relatively high shear rates (100-1000 s-1), they underestimate the elevated viscosity (e.g., in sickle cell anemia) at low shear rates (<10 s-1). Our work provides mechanistic rationale for selective usage of heated RBC as a potentially useful model for studying the abnormal traversal dynamics and hemorheology in certain blood disorders.