Microstructure aware modeling of biochemical transport in arterial blood clots

Investigating clot-flow interactions by integrating intravital imaging with in silico modeling for analysis of flow, transport, and hemodynamic forces

As a blood clot forms, grows, deforms, and embolizes following a vascular injury, local clot-flow interactions lead to a highly dynamic flow environment. The local flow influences transport of biochemical species relevant for clotting, and determines the forces on the clot that in turn lead to clot deformation and embolization. Despite this central role, quantitative characterization of this dynamic clot-flow interaction and flow environment in the clot neighborhood remains a major challenge. Here, we propose an approach that integrates dynamic intravital imaging with computer geometric modeling and computational flow and transport modeling to develop a unified in silico framework to quantify the dynamic clot-flow interactions. We outline the development of the methodology referred to as Intravital Integrated In Silico Modeling or IVISim, and then demonstrate the method on a sample set of simulations comprising clot formation following laser injury in two mouse cremaster arteriole injury model data: one wild-type mouse case, and one diYF knockout mouse case. Simulation predictions are verified against experimental observations of transport of caged fluorescent Albumin (cAlb) in both models. Through these simulations, we illustrate how the IVISim methodology can provide insights into hemostatic processes, the role of flow and clot-flow interactions, and enable further investigations comparing and contrasting different biological model scenarios and parameter variations.

Decoding thrombosis through code: a review of computational models

From the molecular level up to a blood vessel, thrombosis and hemostasis involves many interconnected biochemical and biophysical processes over a wide range of length and time scales. Computational modeling has gained eminence in offering insights into these processes beyond what can be obtained from in vitro or in vivo experiments, or clinical measurements. The multiscale and multiphysics nature of thrombosis has inspired a wide range of modeling approaches that aim to address how a thrombus forms and dismantles. Here, we review recent advances in computational modeling with a focus on platelet-based thrombosis. We attempt to summarize the diverse range of modeling efforts straddling the wide-spectrum of physical phenomena, length scales, and time scales; highlighting key advancements and insights from existing studies. Potential information gleaned from models is discussed, ranging from identification of thrombus-prone regions in patient-specific vasculature to modeling thrombus deformation and embolization in response to fluid forces. Furthermore, we highlight several limitations of current models, future directions in the field, and opportunities for clinical translation, to illustrate the state-of-the-art. There are a plethora of opportunity areas for which models can be expanded, ranging from topics of thromboinflammation to platelet production and clearance. Through successes demonstrated in existing studies described here, as well as continued advancements in computational methodologies and computer processing speeds and memory, in silico investigations in thrombosis are poised to bring about significant knowledge growth in the years to come.

Image-based flow simulation of platelet aggregates under different shear rates

Hemodynamics is crucial for the activation and aggregation of platelets in response to flow-induced shear. In this paper, a novel image-based computational model simulating blood flow through and around platelet aggregates is presented. The microstructure of aggregates was captured by two different modalities of microscopy images of in vitro whole blood perfusion experiments in microfluidic chambers coated with collagen. One set of images captured the geometry of the aggregate outline, while the other employed platelet labelling to infer the internal density. The platelet aggregates were modelled as a porous medium, the permeability of which was calculated with the Kozeny-Carman equation. The computational model was subsequently applied to study hemodynamics inside and around the platelet aggregates. The blood flow velocity, shear stress and kinetic force exerted on the aggregates were investigated and compared under 800 s-1, 1600 s-1 and 4000 s-1 wall shear rates. The advection-diffusion balance of agonist transport inside the platelet aggregates was also evaluated by local Péclet number. The findings show that the transport of agonists is not only affected by the shear rate but also significantly influenced by the microstructure of the aggregates. Moreover, large kinetic forces were found at the transition zone from shell to core of the aggregates, which could contribute to identifying the boundary between the shell and the core. The shear rate and the rate of elongation flow were investigated as well. The results imply that the emerging shapes of aggregates are highly correlated to the shear rate and the rate of elongation. The framework provides a way to incorporate the internal microstructure of the aggregates into the computational model and yields a better understanding of the hemodynamics and physiology of platelet aggregates, hence laying the foundation for predicting aggregation and deformation under different flow conditions.