Simulation of platelets suspension flowing through a stenosis model using a dissipative particle dynamics approach

Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications

Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.

Coarse Grained Modeling of Multiphase Flows with Surfactants

Coarse-grained modeling methods allow simulations at larger scales than molecular dynamics, making it feasible to simulate multifluid systems. It is, however, critical to use model parameters that represent the fluid properties with fidelity under both equilibrium and dynamic conditions. In this work, dissipative particle dynamics (DPD) methods were used to simulate the flow of oil and water in a narrow slit under Poiseuille and Couette flow conditions. Large surfactant molecules were also included in the computations. A systematic methodology is presented to determine the DPD parameters necessary for ensuring that the boundary conditions were obeyed, that the oil and water viscosities were represented correctly, and that the velocity profile for the multifluid system agreed with the theoretical expectations. Surfactant molecules were introduced at the oil–water interface (sodium dodecylsulfate and octaethylene glycol monododecyl ether) to determine the effects of surface-active molecules on the two-phase flow. A critical shear rate was found for Poiseuille flow, beyond which the surfactants desorbed to form the interface forming micelles and destabilize the interface, and the surfactant-covered interface remained stable under Couette flow even at high shear rates.

Shear-Mediated Platelet Activation is Accompanied by Unique Alterations in Platelet Release of Lipids

INTRODUCTION

Platelet activation by mechanical means such as shear stress exposure, is a vital driver of thrombotic risk in implantable blood-contacting devices used in the treatment of heart failure. Lipids are essential in platelets activation and have been studied following biochemical activation. However, little is known regarding lipid alterations occurring with mechanical shear-mediated platelet activation.

METHODS

Here, we determined if shear-activation of platelets induced lipidome changes that differ from those associated with biochemically-mediated platelet activation. We performed high-resolution lipidomic analysis on purified platelets from four healthy human donors. For each donor, we compared the lipidome of platelets that were non-activated or activated by shear, ADP, or thrombin treatment.

RESULTS

We found that shear activation altered cell-associated lipids and led to the release of lipids into the extracellular environment. Shear-activated platelets released 21 phospholipids and sphingomyelins at levels statistically higher than platelets activated by biochemical stimulation.

CONCLUSIONS

We conclude that shear-mediated activation of platelets alters the basal platelet lipidome. Further, these alterations differ and are unique in comparison to the lipidome of biochemically activated platelets. Many of the released phospholipids contained an arachidonic acid tail or were phosphatidylserine lipids, which have known procoagulant properties. Our findings suggest that lipids released by shear-activated platelets may contribute to altered thrombosis in patients with implanted cardiovascular therapeutic devices.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12195-021-00692-x.

A multiscale model for multiple platelet aggregation in shear flow

We developed a multiscale model for simulating aggregation of multiple, free-flowing platelets in low-intermediate shear viscous flow, in which aggregation is mediated by the interaction of αIIbβ3 receptors on the platelet membrane and fibrinogen (Fg). This multiscale model uses coarse grained molecular dynamics (CGMD) for platelets at the microscales and dissipative particle dynamics (DPD) for the shear flow at the macroscales, employing our hybrid aggregation force field for modeling molecular level receptor ligand bonds. We define an aggregation tensor and use it to quantify the molecular level contact characteristics between platelets in an aggregate. We perform numerical studies under different flow conditions for platelet doublets and triplets and evaluate the contact area, detaching force and minimum distance between different pairs of platelets in an aggregate. We also present the dynamics of applied stress and velocity magnitude distributions on the platelet membrane during aggregation and quantify the increase in stress in the contact region under different flow conditions. Integrating the knowledge from our previously validated models, together with new aggregation scenarios, our model can dynamically quantify aggregation characteristics and map stress and velocity distribution on the platelet membrane which are difficult to measure in vitro, thus providing an insight into mechanotransduction bond formation response of platelets to flow-induced shear stresses. This modeling framework, together with the tensor method for quantifying inter-platelet contact, can be extended to simulate and analyze larger aggregates and their adhesive properties.

Shear-mediated platelet activation in the free flow II: Evolving mechanobiological mechanisms reveal an identifiable signature of activation and a bi-directional platelet dyscrasia with thrombotic and bleeding features

Shear-mediated platelet activation (SMPA) in the "free flow" is the net result of a range of cell mechanobiological mechanisms. Previously, we outlined three main groups of mechanisms including: 1) mechano-destruction - i.e. additive platelet (membrane) damage; 2) mechano-activation - i.e. activation of shear-sensitive ion channels and pores; and 3) mechano-transduction - i.e. "outside-in" signaling via a range of transducers. Here, we report on recent advances since our original report which describes additional features of SMPA. A clear "signature" of SMPA has been defined, allowing differentiation from biochemically-mediated activation. Notably, SMPA is characterized by mitochondrial dysfunction, platelet membrane eversion, externalization of anionic phospholipids, and increased thrombin generation on the platelet surface. However, SMPA does not lead to integrin αIIbβ3 activation or P-selectin exposure due to platelet degranulation, as is commonly observed in biochemical activation. Rather, downregulation of GPIb, αIIbβ3, and P-selectin surface expression is evident. Furthermore, SMPA is accompanied by a decrease in overall platelet size coupled with a concomitant, progressive increase in microparticle generation. Shear-ejected microparticles are highly enriched in GPIb and αIIbβ3. These observations indicate the enhanced diffusion, migration, or otherwise dispersion of platelet adhesion receptors to membrane zones, which are ultimately shed as receptor-rich PDMPs. The pathophysiological consequence of this progressive shear accumulation phenomenon is an associated dyscrasia of remaining platelets - being both reduced in size and less activatable via biochemical means - a tendency to favor bleeding, while concomitantly shed microparticles are highly prothrombotic and increase the tendency for thrombosis in both local and systemic milieu. These mechanisms and observations offer direct clinical utility in allowing measurement and guidance of the net balance of platelet driven events in patients with implanted cardiovascular therapeutic devices.

Device-Induced Hemostatic Disorders in Mechanically Assisted Circulation

Mechanically assisted circulation (MAC) sustains the blood circulation in the body of a patients undergoing cardiac surgery with cardiopulmonary bypass (CPB) or on ventricular assistance with a ventricular assist device (VAD) or on extracorporeal membrane oxygenation (ECMO) with a pump-oxygenator system. While MAC provides short-term (days to weeks) support and long-term (months to years) for the heart and/or lungs, the blood is inevitably exposed to non-physiological shear stress (NPSS) due to mechanical pumping action and in contact with artificial surfaces. NPSS is well known to cause blood damage and functional alterations of blood cells. In this review, we discussed shear-induced platelet adhesion, platelet aggregation, platelet receptor shedding, and platelet apoptosis, shear-induced acquired von Willebrand syndrome (AVWS), shear-induced hemolysis and microparticle formation during MAC. These alterations are associated with perioperative bleeding and thrombotic events, morbidity and mortality, and quality of life in MCS patients. Understanding the mechanism of shear-induce hemostatic disorders will help us develop low-shear-stress devices and select more effective treatments for better clinical outcomes.

Artificial Intelligence for Accelerating Time Integrations in Multiscale Modeling

We developed a novel data-driven Artificial Intelligence-enhanced Adaptive Time Stepping algorithm (AI-ATS) that can adapt timestep sizes to underlying biophysical dynamics. We demonstrated its values in solving a complex biophysical problem, at multiple spatiotemporal scales, that describes platelet dynamics in shear blood flow. In order to achieve a significant speedup of this computationally demanding problem, we integrated a framework of novel AI algorithms into the solution of the platelet dynamics equations. Our framework involves recurrent neural network-based autoencoders by the Long Short-Term Memory and the Gated Recurrent Units as the first step for memorizing the dynamic states in long-term dependencies for the input time series, followed by two fully-connected neural networks to optimize timestep sizes and step jumps. The computational efficiency of our AI-ATS is underscored by assessing the accuracy and speed of a multiscale simulation of the platelet with the standard time stepping algorithm (STS). By adapting the timestep size, our AI-ATS guides the omission of multiple redundant time steps without sacrificing significant accuracy of the dynamics. Compared to the STS, our AI-ATS achieved a reduction of 40% unnecessary calculations while bounding the errors of mechanical and thermodynamic properties to 3%.

A predictive multiscale model for simulating flow-induced platelet activation: Correlating in silico results with in vitro results

Flow-induced platelet activation prompts complex filopodial formation. Continuum methods fail to capture such molecular-scale mechanisms. A multiscale numerical model was developed to simulate this activation process, where a Dissipative Particle Dynamics (DPD) model of viscous blood flow is interfaced with a Coarse Grained Molecular Dynamics (CGMD) platelet model. Embedded in DPD blood flow, the macroscopic dynamic stresses are interactively transferred to the CGMD model, inducing intra-platelet associated events. The platelets activate by a biomechanical transductive linkage chain and dynamically change their shape in response. The models are fully coupled via a hybrid-potential interface and multiple time-stepping (MTS) schemes for handling the disparity between the spatiotemporal scales. Cumulative hemodynamic stresses that may lead to platelet activation are mapped on the surface membrane and simultaneously transmitted to the cytoplasm and cytoskeleton. Upon activation, the flowing platelets lose their quiescent discoid shape and evolve by forming filopodia. The model predictions were validated by a set of in vitro experiments, Platelets were exposed to various combinations of shear stresses and durations in our programmable hemodynamic shearing device (HSD). Their shape change was measured at multiple time points using scanning electron microscopy (SEM). The CGMD model parameters were fine-tuned by interrogating a parameter space established in these experiments. Segmentation of the SEM imaging streams was conducted by a deep machine learning system. This model can be further employed to simulate shear mediated platelet activation thrombosis initiation and to study the effects of modulating platelet properties to enhance their shear resistance via mechanotransduction pathways.

Numerical study on the adhesion of a circulating tumor cell in a curved microvessel

The adhesion of a circulating tumor cell (CTC) in a three-dimensional curved microvessel was numerically investigated. Simulations were first performed to characterize the differences in the dynamics and adhesion of a CTC in the straight and curved vessels. After that, a parametric study was performed to investigate the effects of the applied driven force density f (or the flow Reynolds number Re) and the CTC membrane bending modulus K_b on the CTC adhesion. Our simulation results show that the CTC prefers to adhere to the curved vessel as more bonds are formed around the transition region of the curved part due to the increased cell-wall contact by the centrifugal force. The parametric study also indicates that when the flow driven force f (or Re) increases or when the CTC becomes softer (K_b decreases), the bond formation probability increases and the bonds will be formed at more sites of a curved vessel. The increased f (or Re) brings a larger centrifugal force, while the decreased K_b generates more contact areas at the cell-wall interface, both of which are beneficial to the bond formation. In the curved vessel, it is found that the site where bonds are formed the most (hotspot) varies with the applied f and the K_b. For our vessel geometry, when f is small, the hotspot tends to be within the first bend of the vessel, while as f increases or K_b decreases, the hotspot may shift to the second bend of the vessel.

Neutrophils promote venular thrombosis by shaping the rheological environment for platelet aggregation

In advanced inflammatory disease, microvascular thrombosis leads to the interruption of blood supply and provokes ischemic tissue injury. Recently, intravascularly adherent leukocytes have been reported to shape the blood flow in their immediate vascular environment. Whether these rheological effects are relevant for microvascular thrombogenesis remains elusive. Employing multi-channel in vivo microscopy, analyses in microfluidic devices, and computational modeling, we identified a previously unanticipated role of leukocytes for microvascular clot formation in inflamed tissue. For this purpose, neutrophils adhere at distinct sites in the microvasculature where these immune cells effectively promote thrombosis by shaping the rheological environment for platelet aggregation. In contrast to larger (lower-shear) vessels, this process in high-shear microvessels does not require fibrin generation or extracellular trap formation, but involves GPIbα-vWF and CD40-CD40L-dependent platelet interactions. Conversely, interference with these cellular interactions substantially compromises microvascular clotting. Thus, leukocytes shape the rheological environment in the inflamed venular microvasculature for platelet aggregation thereby effectively promoting the formation of blood clots. Targeting this specific crosstalk between the immune system and the hemostatic system might be instrumental for the prevention and treatment of microvascular thromboembolic pathologies, which are inaccessible to invasive revascularization strategies.

HEMODYNAMIC INVESTIGATION OF SUBCLAVIAN-CORONARY STEAL SYNDROME IN DIALYSIS PATIENTS WITH CORONARY ARTERY OCCLUSION AND DIFFERENT STENOSIS PERCENTAGES IN SUBCLAVIAN ARTERY

Purpose: Diseases of the coronary arteries represent critical cardiovascular bleakness and mortality around the world. The main focus of this study is to simulate hemodialysis patients undergoing Coronary Artery Bypass Graft (CABG) surgery by the left internal mammary artery. The survey is to investigate whether subclavian-coronary steal syndrome will occur in patients with stenosis in their left subclavian.

Methods: A three-dimensional model of Brachial—Radial—Ulnar tree arterial system and dialysis graft are constructed and then simulated by the implementation of fluid–structure interaction (FSI) and non-FSI models. After consent was obtained, data attributed to dialysis patients who had undergone coronary bypass surgery at the Tehran Heart Center (THD) was collected.

Results: Stenosis in the subclavian artery causes the left internal mammary artery (LIMA) graft to lack flow that, in CABG cases, is supplying the heart muscle. With the increase of stenosis from zero to 54%, the flow results show a negative flow indicating reversed flow in the artery. Meanwhile, the comparison between flow quantity of a normal Left Anterior Descending (LAD) and LIMA in case of 38% stenosis in subclavian shows that the amount of flow reaching the heart is less than standard. The threshold stenosis in the subclavian artery for which subclavian-coronary steal phenomena occurs is simulated to be around 54%. Furthermore, investigation of the effect of flow quantity in cases of hemodialysis versus nonhemodialysis shows that the flow is lessened in case of hemodialysis.

Conclusions: In hemodialysis patients with more than 38% stenosis in left subclavian, it is well advised that LIMA bypass graft is not the preferred graft choice, and at 54% of subclavian stenosis, the steal phenomenon occurs.

A Multiscale Model for Recruitment Aggregation of Platelets by Correlating with In Vitro Results

Introduction

We developed a multiscale model to simulate the dynamics of platelet aggregation by recruitment of unactivated platelets flowing in viscous shear flows by an activated platelet deposited onto a blood vessel wall. This model uses coarse grained molecular dynamics (CGMD) for platelets at the microscale and dissipative particle dynamics (DPD) for the shear flow at the macroscale. Under conditions of relatively low shear, aggregation is mediated by fibrinogen via αIIbβ3 receptors.

Methods

The binding of αIIbβ3 and fibrinogen is modeled by a molecular-level hybrid force field consisting of Morse potential and Hooke law for the nonbonded and bonded interactions, respectively. The force field, parametrized in two different interaction scales, is calculated by correlating with the platelet contact area measured in vitro and the detaching force between αIIbβ3 and fibrinogen.

Results

Using our model, we derived, the relationship between recruitment force and distance between the centers of mass of two platelets, by integrating the molecular-scale inter-platelet interactions during recruitment aggregation in shear flows. Our model indicates that assuming a rigid-platelet model, underestimates the contact area by 89% and the detaching force by 93% as compared to a model that takes into account the platelet deformability leading to a prediction of a significantly lower attachment during recruitment.

Conclusions

The molecular-level predictive capability of our model sheds a light on differences observed between transient and permanent platelet aggregation patterns. The model and simulation framework can be further adapted to simulate initial thrombus formation involving multiple flowing platelets as well as deposition and adhesion onto blood vessels.

Quantifying Platelet Margination in Diabetic Blood Flow

Patients with type 2 diabetes mellitus (T2DM) develop thrombotic abnormalities strongly associated with cardiovascular diseases. In addition to the changes of numerous coagulation factors such as elevated levels of thrombin and fibrinogen, the abnormal rheological effects of red blood cells (RBCs) and platelets flowing in blood are crucial in platelet adhesion and thrombus formation in T2DM. An important process contributing to the latter is the platelet margination. We employ the dissipative particle dynamics method to seamlessly model cells, plasma, and vessel walls. We perform a systematic study on RBC and platelet transport in cylindrical vessels by considering different cell shapes, sizes, and RBC deformabilities in healthy and T2DM blood, as well as variable flowrates and hematocrit. In particular, we use cellular-level RBC and platelet models with parameters derived from patient-specific data and present a sensitivity study. We find T2DM RBCs, which are less deformable compared to normal RBCs, lower the transport of platelets toward the vessel walls, whereas platelets with higher mean volume (often observed in T2DM) lead to enhanced margination. Furthermore, increasing the flowrate or hematocrit enhances platelet margination. We also investigated the effect of platelet shape and observed a nonmonotonic variation with the highest near-wall concentration corresponding to platelets with a moderate aspect ratio of 0.38. We examine the role of white blood cells (WBCs), whose count is increased notably in T2DM patients. We find that WBC rolling or WBC adhesion tends to decrease platelet margination due to hydrodynamic effects. To the best of our knowledge, such simulations of blood including all blood cells have not been performed before, and our quantitative findings can help separate the effects of hydrodynamic interactions from adhesive interactions and potentially shed light on the associated pathological processes in T2DM such as increased inflammatory response, platelet activation and adhesion, and ultimately thrombus formation.

A multiscale biomechanical model of platelets: Correlating with in-vitro results

Using dissipative particle dynamics (DPD) combined with coarse grained molecular dynamics (CGMD) approaches, we developed a multiscale deformable platelet model to accurately describe the molecular-scale intra-platelet constituents and biomechanical properties of platelets in blood flow. Our model includes the platelet bilayer membrane, cytoplasm and an elaborate elastic cytoskeleton. Correlating numerical simulations with published in-vitro experiments, we validated the biorheology of the cytoplasm, the elastic response of membrane to external stresses, and the stiffness of the cytoskeleton actin filaments, resulting in an accurate representation of the molecular-level biomechanical microstructures of platelets. This enabled us to study the mechanotransduction process of the hemodynamic stresses acting onto the platelet membrane and transmitted to these intracellular constituents. The platelets constituents continuously deform in response to the flow induced stresses. To the best of our knowledge, this is the first molecular-scale platelet model that can be used to accurately predict platelets activation mechanism leading to thrombus formation in prosthetic cardiovascular devices and in vascular disease processes. This model can be further employed to study the effects of novel therapeutic approaches of modulating platelet properties to enhance their shear resistance via mechanotransduction pathways.

Shear-mediated platelet activation in the free flow: Perspectives on the emerging spectrum of cell mechanobiological mechanisms mediating cardiovascular implant thrombosis

Shear-mediated platelet activation (SMPA) is central in thrombosis of implantable cardiovascular therapeutic devices. Despite the morbidity and mortality associated with thrombosis of these devices, our understanding of mechanisms operative in SMPA, particularly in free flowing blood, remains limited. Herein we present and discuss a range of emerging mechanisms for consideration for "free flow" activation under supraphysiologic shear. Further definition and manipulation of these mechanisms will afford opportunities for novel pharmacologic and mechanical strategies to limit SMPA and enhance overall implant device safety.

Sub-cellular modeling of platelet transport in blood flow through microchannels with constriction

Platelet transport through arterial constrictions is one of the controlling processes influencing their adhesive functions and the formation of thrombi. We perform high-fidelity mesoscopic simulations of blood flow in microchannels with constriction, resembling arterial stenoses. The wall shear rates inside the constrictions reach levels as high as ≈8000 s(-1), similar to those encountered in moderate atherosclerotic plaques. Both red blood cells and platelets are resolved at sub-cellular resolution using the Dissipative Particle Dynamics (DPD) method. We perform a systematic study on the red blood cell and platelet transport by considering different levels of constriction, blood hematocrit and flow rates. We find that higher levels of constriction and wall shear rates lead to significantly enhanced margination of platelets, which may explain the experimental observations of enhanced post-stenosis platelet aggregation. We also observe similar margination effects for stiff particles of spherical shapes such as leukocytes. To our knowledge, such numerical simulations of dense blood through complex geometries have not been performed before, and our quantitative findings could shed new light on the associated physiological processes such as ATP release, plasma skimming, and thrombus formation.

Numerical simulation of a single cell passing through a narrow slit

The narrow slit between endothelial cells that line the microvessel wall is the principal pathway for tumor cell extravasation to the surrounding tissue. To understand this crucial step for tumor hematogenous metastasis, we used dissipative particle dynamics method to investigate an individual cell passing through a narrow slit numerically. The cell membrane was simulated by a spring-based network model which can separate the internal cytoplasm and surrounding fluid. The effects of the cell elasticity, cell shape, nucleus and slit size on the cell transmigration through the slit were investigated. Under a fixed driving force, the cell with higher elasticity can be elongated more and pass faster through the slit. When the slit width decreases to 2/3 of the cell diameter, the spherical cell becomes jammed despite reducing its elasticity modulus by 10 times. However, transforming the cell from a spherical to ellipsoidal shape and increasing the cell surface area by merely 9.3 % can enable the cell to pass through the narrow slit. Therefore, the cell shape and surface area increase play a more important role than the cell elasticity in cell passing through the narrow slit. In addition, the simulation results indicate that the cell migration velocity decreases during entrance but increases during exit of the slit, which is qualitatively in agreement with the experimental observation.

Damage in spherical cellular membrane generated by the shock waves: coarse-grained molecular dynamics simulation of lipid vesicle

Traumatic Brain Injury is a major health issue that is hard to diagnose since it often occurs without signs of external injuries. While it is well known that exposure of biological cells to shock waves causes damage to the cell membrane, it is currently unknown by which mechanisms damage is caused, and how it depends on physical parameters such as shock wave velocity, shock pulse duration, or shock pulse shape. In this computational study, we use a coarse-grained model of the lipid vesicle as a generic model of a cell membrane to elucidate the general principles of the cellular damage induced by the shock wave direct passage through the cranium. Results indicate that the extent of the liposome compression does not strongly depend on the pressure pulse and that liposome extension is very sensitive to the change in the negative pressure phase. The structural integrity of the vesicle is altered as pores form in the lipid membrane at overall pressure impulses generated by supersonic shock waves, which are greater than 5 Pa·s at single or repetitive exposure. Consequently, these permeability changes may lead to changes in the influx of sodium, potassium, and calcium ions.

A Multiple Time Stepping Algorithm for Efficient Multiscale Modeling of Platelets Flowing in Blood Plasma

We developed a multiple time-stepping (MTS) algorithm for multiscale modeling of the dynamics of platelets flowing in viscous blood plasma. This MTS algorithm improves considerably the computational efficiency without significant loss of accuracy. This study of the dynamic properties of flowing platelets employs a combination of the dissipative particle dynamics (DPD) and the coarse-grained molecular dynamics (CGMD) methods to describe the dynamic microstructures of deformable platelets in response to extracellular flow-induced stresses. The disparate spatial scales between the two methods are handled by a hybrid force field interface. However, the disparity in temporal scales between the DPD and CGMD that requires time stepping at microseconds and nanoseconds respectively, represents a computational challenge that may become prohibitive. Classical MTS algorithms manage to improve computing efficiency by multi-stepping within DPD or CGMD for up to one order of magnitude of scale differential. In order to handle 3-4 orders of magnitude disparity in the temporal scales between DPD and CGMD, we introduce a new MTS scheme hybridizing DPD and CGMD by utilizing four different time stepping sizes. We advance the fluid system at the largest time step, the fluid-platelet interface at a middle timestep size, and the nonbonded and bonded potentials of the platelet structural system at two smallest timestep sizes. Additionally, we introduce parameters to study the relationship of accuracy versus computational complexities. The numerical experiments demonstrated 3000x reduction in computing time over standard MTS methods for solving the multiscale model. This MTS algorithm establishes a computationally feasible approach for solving a particle-based system at multiple scales for performing efficient multiscale simulations.

Mechanical platelet activation potential in abdominal aortic aneurysms

Intraluminal thrombus (ILT) in abdominal aortic aneurysms (AAA) has potential implications to aneurysm growth and rupture risk; yet, the mechanisms underlying its development remain poorly understood. Some researchers have proposed that ILT development may be driven by biomechanical platelet activation within the AAA, followed by adhesion in regions of low wall shear stress. Studies have investigated wall shear stress levels within AAA, but platelet activation potential (AP) has not been quantified. In this study, patient-specific computational fluid dynamic (CFD) models were used to analyze stress-induced AP within AAA under rest and exercise flow conditions. The analysis was conducted using Lagrangian particle-based and Eulerian continuum-based approaches, and the results were compared. Results indicated that biomechanical platelet activation is unlikely to play a significant role for the conditions considered. No consistent trend was observed in comparing rest and exercise conditions, but the functional dependence of AP on stress magnitude and exposure time can have a large impact on absolute levels of anticipated platelet AP. The Lagrangian method obtained higher peak AP values, although this difference was limited to a small percentage of particles that falls below reported levels of physiologic background platelet activation.

Multiscale Particle-Based Modeling of Flowing Platelets in Blood Plasma Using Dissipative Particle Dynamics and Coarse Grained Molecular Dynamics

We developed a multiscale particle-based model of platelets, to study the transport dynamics of shear stresses between the surrounding fluid and the platelet membrane. This model facilitates a more accurate prediction of the activation potential of platelets by viscous shear stresses - one of the major mechanisms leading to thrombus formation in cardiovascular diseases and in prosthetic cardiovascular devices. The interface of the model couples coarse-grained molecular dynamics (CGMD) with dissipative particle dynamics (DPD). The CGMD handles individual platelets while the DPD models the macroscopic transport of blood plasma in vessels. A hybrid force field is formulated for establishing a functional interface between the platelet membrane and the surrounding fluid, in which the microstructural changes of platelets may respond to the extracellular viscous shear stresses transferred to them. The interaction between the two systems preserves dynamic properties of the flowing platelets, such as the flipping motion. Using this multiscale particle-based approach, we have further studied the effects of the platelet elastic modulus by comparing the action of the flow-induced shear stresses on rigid and deformable platelet models. The results indicate that neglecting the platelet deformability may overestimate the stress on the platelet membrane, which in turn may lead to erroneous predictions of the platelet activation under viscous shear flow conditions. This particle-based fluid-structure interaction multiscale model offers for the first time a computationally feasible approach for simulating deformable platelets interacting with viscous blood flow, aimed at predicting flow induced platelet activation by using a highly resolved mapping of the stress distribution on the platelet membrane under dynamic flow conditions.