Diffusion and reaction in the cell glycocalyx and the extracellular matrix

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.

Simulation of intrathrombus fluid and solute transport using in vivo clot structures with single platelet resolution

The mouse laser injury thrombosis model provides up to 0.22 μm-resolved voxel information about the pore architecture of the dense inner core and loose outer shell regions of an in vivo arterial thrombus. Computational studies were conducted on this 3D structure to quantify transport within and around the clot: Lattice Boltzmann method defined vessel hemodynamics, while passive Lagrangian Scalar Tracking with Brownian motion contribution simulated diffusive-convective transport of various inert solutes (released from lumen or the injured wall). For an input average lumen blood velocity of 0.478 cm/s (measured by Doppler velocimetry), a 0.2 mm/s mean flow rate was obtained within the thrombus structure, most of which occurred in the 100-fold more permeable outer shell region (calculated permeability of the inner core was 10(-11) cm(2)). Average wall shear stresses were 80-100 dyne/cm(2) (peak values >200 dyne/cm(2)) on the outer rough surface of the thrombus. Within the thrombus, small molecule tracers (0.1 kDa) experienced ~70,000 collisions/s and penetrated/exited it in about 1 s, whereas proteins (~50 kDa) had ~9000 collisions/s and required about 10 s (tortuosity ~2-2.5). These simulations help define physical processes during thrombosis and constraints for drug delivery to the thrombus.

Simulation of NO and O2 transport facilitated by polymerized hemoglobin solutions in an arteriole that takes into account wall shear stress-induced NO production

A mathematical model was developed to study nitric oxide (NO) and oxygen (O(2)) transport in an arteriole and surrounding tissues exposed to a mixture of red blood cells (RBCs) and hemoglobin (Hb)-based O(2) carriers (HBOCs). A unique feature of this model is the inclusion of blood vessel wall shear stress-induced production of endothelial-derived NO, which is very sensitive to the viscosity of the RBC and HBOC mixture traversing the blood vessel lumen. Therefore in this study, a series of polymerized bovine Hb (PolyHb) solutions with high viscosity, varying O(2) affinities, NO dioxygenation rate constants and O(2) dissociation rate constants that were previously synthesized and characterized by our group was evaluated via mathematical modeling, in order to investigate the effect of these biophysical properties on the transport of NO and O(2) in an arteriole and its surrounding tissues subjected to anemia with the commercial HBOC Oxyglobin® and cell-free bovine Hb (bHb) serving as appropriate controls. The computer simulation results indicated that transfusion of high viscosity PolyHb solutions promoted blood vessel wall shear stress dependent generation of the vasodilator NO, especially in the blood vessel wall and should transport enough NO inside the smooth muscle layer to activate vasodilation compared to the commercial HBOC Oxyglobin® and cell-free bHb. However, NO scavenging in the arteriole lumen was unavoidable due to the intrinsic high NO dioxygenation rate constant of the HBOCs being studied. This study also observed that all PolyHbs could potentially improve tissue oxygenation under hypoxic conditions, while low O(2) affinity PolyHbs were more effective in oxygenating tissues under normoxic conditions compared with high O(2) affinity PolyHbs. In addition, all ultrahigh molecular weight PolyHbs displayed higher O(2) transfer rates than the commercial HBOC Oxyglobin® and cell-free bHb. Therefore, these results suggest that ultrahigh molecular weight PolyHb solutions could be used as safe and efficacious O(2) carriers for use in transfusion medicine. It also suggests that future generations of PolyHb solutions should possess lower NO dioxygenation reaction rate constants in order to reduce NO scavenging, while maintaining high solution viscosity to take advantage of wall shear stress-induced NO production. Taken together, we suggest that this mathematical model can be used to predict the vasoactivity of HBOCs and help guide the design and optimization of the next generation of HBOCs for use in transfusion medicine.