SPH simulations of WBC adhesion to the endothelium: the role of haemodynamics and endothelial binding kinetics

Investigation of the hemodynamic flow conditions and blood-induced stresses inside an abdominal aortic aneurysm by means of a SPH numerical model

The estimation of blood flow-induced loads occurring on the artery wall is affected by uncertainties hidden in the complex interaction of the pulsatile flow, the mechanical parameters of the artery, and the external support conditions. To circumvent these difficulties, a specific tool is developed by combining the aorta displacements measured by an electrocardiogram-gated-computed tomography angiography, with the blood velocity field computed by a smoothed particle hydrodynamics (SPH) numerical model. In the present work, the SPH model has been specifically adapted to the solution of the 3D Navier-Stokes equations inside a domain with boundaries of prescribed motion. Images of the abdominal aorta aneurysm (AAA) of a 44-year-old female patient were acquired during a stabilized cardiac cycle by electrocardiogram-gated-computed tomography angiography. The in vivo kinematic field inside the pulsating arterial wall was estimated by using recent technology, which makes it possible to follow the shape of the arterial wall during a cardiac cycle. We compare the flow conditions and the blood-induced loads, computed by the numerical model under the assumption of a moving arterial wall, with the corresponding results obtained assuming three rigid wall geometries of the vessel during the cardiac cycle. Significant differences were found for the wall shear stress distribution.

The key events of thrombus formation: platelet adhesion and aggregation

Thrombus formation is a complex, dynamic and multistep process, involving biochemical reactions, mechanical stimulation, hemodynamics, and so on. In this study, we concentrate on its two crucial steps: (i) platelets adhered to a vessel wall, or simply platelet adhesion, and (ii) platelets clumping and arrested to the adherent platelets, named platelet aggregation. We report the first direct simulation of three modes of platelet adhesion, detachment, rolling adhesion and firm adhesion, as well as the formation, disintegration, arrestment and consolidation of platelet plugs. The results show that the bond dissociation in the detachment mode is mainly attributed to a high probability of rupturing bonds, such that any existing bond can be quickly ruptured and all bonds would be completely broken. In the rolling adhesion, however, it is mainly attributed to the strong traction from the shear flow or erythrocytes, causing that the bonds are ruptured at the trailing edge of the platelet. The erythrocytes play an important role in platelet activities, such as the formation, disintegration, arrestment and consolidation of platelet plugs. They exert an aggregate force on platelets, a repulsion at a near distance but an attraction at a far distance to the platelets. This aggregate force can promote platelets to form a plug and/or bring along a part of a platelet plug causing its disintegration. It also greatly influences the arrestment and consolidation of platelet plugs, together with the adhesive force from the thrombus.

Modeling thrombosis in silico: Frontiers, challenges, unresolved problems and milestones

Hemostasis is a complex physiological mechanism that functions to maintain vascular integrity under any conditions. Its primary components are blood platelets and a coagulation network that interact to form the hemostatic plug, a combination of cell aggregate and gelatinous fibrin clot that stops bleeding upon vascular injury. Disorders of hemostasis result in bleeding or thrombosis, and are the major immediate cause of mortality and morbidity in the world. Regulation of hemostasis and thrombosis is immensely complex, as it depends on blood cell adhesion and mechanics, hydrodynamics and mass transport of various species, huge signal transduction networks in platelets, as well as spatiotemporal regulation of the blood coagulation network. Mathematical and computational modeling has been increasingly used to gain insight into this complexity over the last 30 years, but the limitations of the existing models remain profound. Here we review state-of-the-art-methods for computational modeling of thrombosis with the specific focus on the analysis of unresolved challenges. They include: a) fundamental issues related to physics of platelet aggregates and fibrin gels; b) computational challenges and limitations for solution of the models that combine cell adhesion, hydrodynamics and chemistry; c) biological mysteries and unknown parameters of processes; d) biophysical complexities of the spatiotemporal networks' regulation. Both relatively classical approaches and innovative computational techniques for their solution are considered; the subjects discussed with relation to thrombosis modeling include coarse-graining, continuum versus particle-based modeling, multiscale models, hybrid models, parameter estimation and others. Fundamental understanding gained from theoretical models are highlighted and a description of future prospects in the field and the nearest possible aims are given.

Effects of membrane deformability and bond formation/dissociation rates on adhesion dynamics of a spherical capsule in shear flow

Cellular adhesion plays a critical role in biological systems and biomedical applications. Cell deformation and biophysical properties of adhesion molecules are of significance for the adhesion behavior. In the present work, dynamic adhesion of a deformable capsule to a planar substrate, in a linear shear flow, is numerically simulated to investigate the combined influence of membrane deformability (quantified by the capillary number) and bond formation/dissociation rates on the adhesion behavior. The computational model is based on the immersed boundary-lattice Boltzmann method for the capsule-fluid interaction and a probabilistic adhesion model for the capsule-substrate interaction. Three distinct adhesion states, detachment, rolling adhesion and firm adhesion, are identified and presented in a state diagram as a function of capillary number and bond dissociation rate. The impact of bond formation rate on the state diagram is further investigated. Results show that the critical bond dissociation rate for the transition of rolling or firm adhesion to detachment is strongly related to the capsule deformability. At the rolling-adhesion state, smaller off rates are needed for larger capillary number to increase the rolling velocity and detach the capsule. In contrast, the critical off rate for firm-to-detach transition slightly increases with the capillary number. With smaller on rate, the effect of capsule deformability on the critical off rates is more pronounced and capsules with moderate deformability are prone to detach by the shear flow. Further increasing of on rate leads to large expansion of both rolling-adhesion and firm-adhesion regions. Even capsules with relatively large deformability can maintain stable rolling adhesion at certain off rate.