A multiscale model of thrombus development

[Progress in the analysis of hemolysis and coagulation models for interventional micro-axial flow blood pumps]

Interventional micro-axial flow blood pump is widely used as an effective treatment for patients with cardiogenic shock. Hemolysis and coagulation are vital concerns in the clinical application of interventional micro-axial flow pumps. This paper reviewed hemolysis and coagulation models for micro-axial flow blood pumps. Firstly, the structural characteristics of commercial interventional micro-axial flow blood pumps and issues related to clinical applications were introduced. Then the basic mechanisms of hemolysis and coagulation were used to study the factors affecting erythrocyte damage and platelet activation in interventional micro-axial flow blood pumps, focusing on the current models of hemolysis and coagulation on different scales (macroscopic, mesoscopic, and microscopic). Since models at different scales have different perspectives on the study of hemolysis and coagulation, a comprehensive analysis combined with multi-scale models is required to fully consider the influence of complex factors of interventional pumps on hemolysis and coagulation.

Balancing competing effects of tissue growth and cytoskeletal regulation during Drosophila wing disc development

How a developing organ robustly coordinates the cellular mechanics and growth to reach a final size and shape remains poorly understood. Through iterations between experiments and model simulations that include a mechanistic description of interkinetic nuclear migration, we show that the local curvature, height, and nuclear positioning of cells in the Drosophila wing imaginal disc are defined by the concurrent patterning of actomyosin contractility, cell-ECM adhesion, ECM stiffness, and interfacial membrane tension. We show that increasing cell proliferation via different growth-promoting pathways results in two distinct phenotypes. Triggering proliferation through insulin signaling increases basal curvature, but an increase in growth through Dpp signaling and Myc causes tissue flattening. These distinct phenotypic outcomes arise from differences in how each growth pathway regulates the cellular cytoskeleton, including contractility and cell-ECM adhesion. The coupled regulation of proliferation and cytoskeletal regulators is a general strategy to meet the multiple context-dependent criteria defining tissue morphogenesis.

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.

A Fluid-Solid-Growth Solver for Cardiovascular Modeling

We implement full, three-dimensional constrained mixture theory for vascular growth and remodeling into a finite element fluid-structure interaction (FSI) solver. The resulting "fluid-solid-growth" (FSG) solver allows long term, patient-specific predictions of changing hemodynamics, vessel wall morphology, tissue composition, and material properties. This extension from short term (FSI) to long term (FSG) simulations increases clinical relevance by enabling mechanobioloigcally-dependent studies of disease progression in complex domains.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

Computational modeling of hypercoagulability in COVID-19

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has infected more than 100 million people worldwide and claimed millions of lives. While the leading cause of mortality in COVID-19 patients is the hypoxic respiratory failure from acute respiratory distress syndrome, there is accumulating evidence that shows excessive coagulation also increases the fatalities in COVID-19. Thus, there is a pressing demand to understand the association between COVID-19-induced hypercoagulability and the extent of formation of undesired blood clots. Mathematical modeling of coagulation has been used as an important tool to identify novel reaction mechanisms and to identify targets for new drugs. Here, we employ the coagulation factor data of COVID-19 patients reported from published studies as inputs for two mathematical models of coagulation to identify how the concentrations of coagulation factors change in these patients. Our simulation results show that while the levels of many of the abnormal coagulation factors measured in COVID-19 patients promote the generation of thrombin and fibrin, two key components of blood clots, the increased level of fibrinogen and then the reduced level of antithrombin are the factors most responsible for boosting the level of fibrin and thrombin, respectively. Altogether, our study demonstrates the potential of mathematical modeling to identify coagulation factors responsible for the increased clot formation in COVID-19 patients where clinical data is scarce.

Precision Navigation of Venous Thrombosis Guided by Viscosity-Activatable Near-Infrared Fluorescence

Thrombus are blood clots formed by abnormal hemostasis in blood vessels and are closely associated with various diseases such as pulmonary embolism, myocardial infarction and stroke. Early diagnosis and treatment of thrombus is the key to reducing the high risk of thrombotic disease. Given that early thrombus is small in early size, free instability, wide regional distribution and fast formation, it is urgent to develop all-inclusive detection methods that combine high signal-to-noise ratio, in situ dynamic and rapid in-depth tissue imaging. Near-infrared (NIR) fluorescence imaging, with its excellent high spatiotemporal resolution and tissue penetration depth, is a powerful technique for direct visualization of thrombotic events in situ. Considering the fibrin highly expressed in the thrombus is a typical thrombotic target. Moreover, the viscosity of the thrombus is markedly higher than its surroundings. Therefore, we developed a fibrin-targeting and viscosity-activating thrombus NIR fluorescent probe (TIR-V) for high-resolution and high-sensitivity in situ lighten-up thrombus. TIR-V has the advantages of good thrombus targeting, significant "off-on" fluorescence specific response to viscosity, bright NIR fluorescence and good biocompatibility. The thrombus is clearly delineated by a high signal-to-noise ratio NIR fluorescence imaging, enabling imaging detection and precise navigation of thrombotic regions. This work demonstrates the potential of TIR-V as a bifunctional probe for definitive diagnostic imaging and direct navigation of thrombotic lesions in clinical applications.

Uncertainty quantification of a thrombosis model considering the clotting assay PFA-100®

Mathematical models of thrombosis are currently used to study clinical scenarios of pathological thrombus formation. As these models become more complex to predict thrombus formation dynamics high computational cost must be alleviated and inherent uncertainties must be assessed. Evaluating model uncertainties allows to increase the confidence in model predictions and identify avenues of improvement for both thrombosis modeling and anti-platelet therapies. In this work, an uncertainty quantification analysis of a multi-constituent thrombosis model is performed considering a common assay for platelet function (PFA-100®). The analysis is facilitated thanks to time-evolving polynomial chaos expansions used as a parametric surrogate for the full thrombosis model considering two quantities of interest; namely, thrombus volume and occlusion percentage. The surrogate is thoroughly validated and provides a straightforward access to a global sensitivity analysis via computation of Sobol' coefficients. Six out of 15 parameters linked to thrombus consitution, vWF activity, and platelet adhesion dynamics were found to be most influential in the simulation variability considering only individual effects; while parameter interactions are highlighted when considering the total Sobol' indices. The influential parameters are related to thrombus constitution, vWF activity, and platelet to platelet adhesion dynamics. The surrogate model allowed to predict realistic PFA-100® closure times of 300,000 virtual cases that followed the trends observed in clinical data. The current methodology could be used including common anti-platelet therapies to identify scenarios that preserve the hematological balance.

Multiphysics and multiscale modeling of microthrombosis in COVID-19

Emerging clinical evidence suggests that thrombosis in the microvasculature of patients with Coronavirus disease 2019 (COVID-19) plays an essential role in dictating the disease progression. Because of the infectious nature of SARS-CoV-2, patients’ fresh blood samples are limited to access for in vitro experimental investigations. Herein, we employ a novel multiscale and multiphysics computational framework to perform predictive modeling of the pathological thrombus formation in the microvasculature using data from patients with COVID-19. This framework seamlessly integrates the key components in the process of blood clotting, including hemodynamics, transport of coagulation factors and coagulation kinetics, blood cell mechanics and adhesive dynamics, and thus allows us to quantify the contributions of many prothrombotic factors reported in the literature, such as stasis, the derangement in blood coagulation factor levels and activities, inflammatory responses of endothelial cells and leukocytes to the microthrombus formation in COVID-19. Our simulation results show that among the coagulation factors considered, antithrombin and factor V play more prominent roles in promoting thrombosis. Our simulations also suggest that recruitment of WBCs to the endothelial cells exacerbates thrombogenesis and contributes to the blockage of the blood flow. Additionally, we show that the recent identification of flowing blood cell clusters could be a result of detachment of WBCs from thrombogenic sites, which may serve as a nidus for new clot formation. These findings point to potential targets that should be further evaluated, and prioritized in the anti-thrombotic treatment of patients with COVID-19. Altogether, our computational framework provides a powerful tool for quantitative understanding of the mechanism of pathological thrombus formation and offers insights into new therapeutic approaches for treating COVID-19 associated thrombosis.

Emerging clinical evidence suggests that thrombosis in the microvasculature of patients with Coronavirus disease 2019 (COVID-19) plays an essential role in dictating the disease progression. We employ a novel multiphysics and multiscale computational framework to investigate the underlying mechanism of the pathological formation of microthrombi and circulating cell clusters in COVID-19. We quantify the contributions of many prothrombotic factors reported in the literature, such as stasis, the derangement in blood coagulation factor levels and activities, inflammatory responses of endothelial cells and leukocytes to the microthrombus formation in COVID-19, through which we identify the potential targets that should be further evaluated, and prioritized in the anti-thrombotic treatment of patients with COVID-19.

An integrated fluid-structure interaction and thrombosis model for type B aortic dissection

False lumen thrombosis (FLT) in type B aortic dissection has been associated with the progression of dissection and treatment outcome. Existing computational models mostly assume rigid wall behavior which ignores the effect of flap motion on flow and thrombus formation within the FL. In this study, we have combined a fully coupled fluid–structure interaction (FSI) approach with a shear-driven thrombosis model described by a series of convection–diffusion reaction equations. The integrated FSI-thrombosis model has been applied to an idealized dissection geometry to investigate the interaction between vessel wall motion and growing thrombus. Our simulation results show that wall compliance and flap motion can influence the progression of FLT. The main difference between the rigid and FSI models is the continuous development of vortices near the tears caused by drastic flap motion up to 4.45 mm. Flap-induced high shear stress and shear rates around tears help to transport activated platelets further to the neighboring region, thus speeding up thrombus formation during the accelerated phase in the FSI models. Reducing flap mobility by increasing the Young’s modulus of the flap slows down the thrombus growth. Compared to the rigid model, the predicted thrombus volume is 25% larger using the FSI-thrombosis model with a relatively mobile flap. Furthermore, our FSI-thrombosis model can capture the gradual effect of thrombus growth on the flow field, leading to flow obstruction in the FL, increased blood viscosity and reduced flap motion. This model is a step closer toward simulating realistic thrombus growth in aortic dissection, by taking into account the effect of intimal flap and vessel wall motion.

Seven challenges in the multiscale modeling of multicellular tissues

The growth and dynamics of multicellular tissues involve tightly regulated and coordinated morphogenetic cell behaviors, such as shape changes, movement, and division, which are governed by subcellular machinery and involve coupling through short- and long-range signals. A key challenge in the fields of developmental biology, tissue engineering and regenerative medicine is to understand how relationships between scales produce emergent tissue-scale behaviors. Recent advances in molecular biology, live-imaging and ex vivo techniques have revolutionized our ability to study these processes experimentally. To fully leverage these techniques and obtain a more comprehensive understanding of the causal relationships underlying tissue dynamics, computational modeling approaches are increasingly spanning multiple spatial and temporal scales, and are coupling cell shape, growth, mechanics, and signaling. Yet such models remain challenging: modeling at each scale requires different areas of technical skills, while integration across scales necessitates the solution to novel mathematical and computational problems. This review aims to summarize recent progress in multiscale modeling of multicellular tissues and to highlight ongoing challenges associated with the construction, implementation, interrogation, and validation of such models. This article is categorized under: Reproductive System Diseases > Computational Models Metabolic Diseases > Computational Models Cancer > Computational Models.

Recent Advances in Computational Modeling of Biomechanics and Biorheology of Red Blood Cells in Diabetes

Diabetes mellitus, a metabolic disease characterized by chronically elevated blood glucose levels, affects about 29 million Americans and more than 422 million adults all over the world. Particularly, type 2 diabetes mellitus (T2DM) accounts for 90-95% of the cases of vascular disease and its prevalence is increasing due to the rising obesity rates in modern societies. Although multiple factors associated with diabetes, such as reduced red blood cell (RBC) deformability, enhanced RBC aggregation and adhesion to the endothelium, as well as elevated blood viscosity are thought to contribute to the hemodynamic impairment and vascular occlusion, clinical or experimental studies cannot directly quantify the contributions of these factors to the abnormal hematology in T2DM. Recently, computational modeling has been employed to dissect the impacts of the aberrant biomechanics of diabetic RBCs and their adverse effects on microcirculation. In this review, we summarize the recent advances in the developments and applications of computational models in investigating the abnormal properties of diabetic blood from the cellular level to the vascular level. We expect that this review will motivate and steer the development of new models in this area and shift the attention of the community from conventional laboratory studies to combined experimental and computational investigations, aiming to provide new inspirations for the development of advanced tools to improve our understanding of the pathogenesis and pathology of T2DM.

A three-dimensional multiscale model for the prediction of thrombus growth under flow with single-platelet resolution

Modeling thrombus growth in pathological flows allows evaluation of risk under patient-specific pharmacological, hematological, and hemodynamical conditions. We have developed a 3D multiscale framework for the prediction of thrombus growth under flow on a spatially resolved surface presenting collagen and tissue factor (TF). The multiscale framework is composed of four coupled modules: a Neural Network (NN) that accounts for platelet signaling, a Lattice Kinetic Monte Carlo (LKMC) simulation for tracking platelet positions, a Finite Volume Method (FVM) simulator for solving convection-diffusion-reaction equations describing agonist release and transport, and a Lattice Boltzmann (LB) flow solver for computing the blood flow field over the growing thrombus. A reduced model of the coagulation cascade was embedded into the framework to account for TF-driven thrombin production. The 3D model was first tested against in vitro microfluidics experiments of whole blood perfusion with various antiplatelet agents targeting COX-1, P₂Y₁, or the IP receptor. The model was able to accurately capture the evolution and morphology of the growing thrombus. Certain problems of 2D models for thrombus growth (artifactual dendritic growth) were naturally avoided with realistic trajectories of platelets in 3D flow. The generalizability of the 3D multiscale solver enabled simulations of important clinical situations, such as cylindrical blood vessels and acute flow narrowing (stenosis). Enhanced platelet-platelet bonding at pathologically high shear rates (e.g., von Willebrand factor unfolding) was required for accurately describing thrombus growth in stenotic flows. Overall, the approach allows consideration of patient-specific platelet signaling and vascular geometry for the prediction of thrombotic episodes.

The excessive formation of blood clots under flow within the circulatory system (thrombosis) is known to initiate heart attacks and strokes. Therefore, obtaining insights into the formation and progression of these clots will be useful in evaluating pharmacological options. To this end, we have developed a 3D computational model that tracks the growth of a blood clot under flow from initial platelet deposition to full vessel occlusion in the presence of soluble platelet agonists. We first validated the model against experimental predictions of blood clots formed in vitro. Due to the construction of the model in 3D, we were able to carry out simulations of clot formation under important clinical situations, namely cylindrical blood vessels and acute flow narrowings (stenoses). To our knowledge, our model is the first of its kind that can account for patient-specific platelet phenotypes to perform robust 3D simulations of thrombus growth in geometries of clinical relevance.

A two-fluid blood stasis model for false lumen thrombosis after type B dissection repair

The formation of thrombosis is a major concern in the false lumen (FL) for post-TEVAR (thoracic endovascular aortic repair) patients. Blood stasis is one of the key factors which lead to the formation of thrombosis in the arterial systems. This study proposed a computational model for blood stasis, using a two-fluid principle to track the locations of blood residual over time. The current study applied this novel model to evaluate blood stasis and thrombosis potential in four patient-specific post-TEVAR FLs of type B aortic dissection, with their follow-up in-vivo observations two years after TEVAR. The locations and topologies of residual blood in the FL predicted by the model agreed well with the in-vivo observations of thrombus. In addition, the results corresponded better with clinical observations in terms of interpatient comparison of degree of thrombosis, compared with conventional hemodynamic parameters. The blood stasis model serves as a valuable addition to conventional metrics to better predict thrombosis potential. Collectively, these metrics can provide an efficient non-invasive method for evaluating blood stasis and thrombosis potential in arterial system, and useful guidance for clinicians' operative planning and postoperative evaluation.

Sensitivity analysis of a reduced model of thrombosis under flow: Roles of Factor IX, Factor XI, and γ'-Fibrin

A highly reduced extrinsic pathway coagulation model (8 ODEs) under flow considered a thin 15-micron platelet layer where transport limitations were largely negligible (except for fibrinogen) and where cofactors (FVIIa, FV, FVIII) were not rate-limiting. By including thrombin feedback activation of FXI and the antithrombin-I activities of fibrin, the model accurately simulated measured fibrin formation and thrombin fluxes. Using this reduced model, we conducted 10,000 Monte Carlo (MC) simulations for ±50% variation of 5 plasma zymogens and 2 fibrin binding sites for thrombin. A sensitivity analysis of zymogen concentrations indicated that FIX activity most influenced thrombin generation, a result expected from hemophilia A and B. Averaging all MC simulations confirmed both the mean and standard deviation of measured fibrin generation on 1 tissue factor (TF) molecule per μm2. Across all simulations, free thrombin in the layer ranged from 20 to 300 nM (mean: 50 nM). The top 2% of simulations that produced maximal fibrin were dominated by conditions with low antithrombin-I activity (decreased weak and strong sites) and high FIX concentration. In contrast, the bottom 2% of simulations that produced minimal fibrin were dominated by low FIX and FX. The percent reduction of fibrin by an ideal FXIa inhibitor (FXI = 0) ranged from 71% fibrin reduction in the top 2% of MC simulations to only 34% fibrin reduction in the bottom 2% of MC simulations. Thus, the antithrombotic potency of FXIa inhibitors may vary depending on normal ranges of zymogen concentrations. This reduced model allowed efficient multivariable sensitivity analysis.

Microstructure aware modeling of biochemical transport in arterial blood clots

Flow-mediated transport of biochemical species is central to thrombotic phenomena. Comprehensive three-dimensional modeling of flow-mediated transport around realistic macroscale thrombi poses challenges owing to their arbitrary heterogeneous microstructure. Here, we develop a microstructure aware model for species transport within and around a macroscale thrombus by devising a custom preconditioned fictitious domain formulation for thrombus-hemodynamics interactions, and coupling it with a fictitious domain advection-diffusion formulation for transport. Microstructural heterogeneities are accounted through a hybrid discrete particle-continuum approach for the thrombus interior. We present systematic numerical investigations on unsteady arterial flow within and around a three-dimensional macroscale thrombus; demonstrate the formation of coherent flow structures around the thrombus which organize advective transport; illustrate the role of the permeation processes at the thrombus boundary and subsequent intra-thrombus transport; and characterize species transport from bulk flow to the thrombus boundary and vice versa.

A 1D-3D Hybrid Model of Patient-Specific Coronary Hemodynamics

PURPOSE

Coronary flow is affected by evolving events such as atherosclerotic plaque formation, rupture, and thrombosis, resulting in myocardial ischemia and infarction. Highly resolved 3D hemodynamic data at the stenosis is essential to model shear-sensitive thrombotic events in coronary artery disease.

METHODS

We developed a hybrid 1D-3D simulation framework to compute patient-specific coronary hemodynamics efficiently. A 1D model of the coronary flow is coupled to an image-based 3D model of the region of interest. This framework affords the advantages of reduced-order modeling, decreasing the global computational cost, without sacrificing the accuracy of the quantities of interest.

RESULTS

We validated our 1D-3D model against full 3D coronary simulations in healthy and diseased conditions. Our results showed good agreement between the 3D and the 1D-3D models while reducing the computational cost by 40-fold compared to the 3D simulation. The 1D-3D model predicted left/right coronary flow distribution within 3% and provided an accurate estimation of fractional flow reserve and wall shear stress distribution at the stenosis comparable to the 3D simulation.

CONCLUSION

Savings in computational cost may be significant in situations with changing geometry, such as growing thrombosis. Also, this approach would allow quantifying the time-dependent effect of thrombotic growth and occlusion on the global coronary circulation.

Computational modeling of blood component transport related to coronary artery thrombosis in Kawasaki disease

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case—mechanical or biochemical—could help define improved strategies for thrombosis prevention tailored for each patient.

Computational studies aiming to model thrombosis often rely on an arterial wall injury. Collagen and other extracellular matrix components are exposed to the bloodstream, which facilitates platelet adhesion to the wall and subsequent clot formation. However, these models are not adequate to explain thrombosis in other settings where even in the absence of a focal lesion, clots may still form under certain flow conditions. Coronary artery aneurysm thrombosis following KD is an example of the need to understand the mechanisms of thrombus initiation in the absence of an extrinsic factor. This study provides a new framework to investigate thrombus initiation in KD from a patient-specific perspective, which integrates fluid mechanics and biochemistry and which could help quantify the pro-coagulant environment induced by the aneurysm and become a predictive tool. The work presented here has broad relevance to other clinical situations where flow stagnation and transport are driving factors in thrombus formation.

An integrated fluid-chemical model toward modeling the thrombus formation in an idealized model of aortic dissection

Type B aortic dissection is a major aortic catastrophe that can be acutely complicated by rapid expansion, rupture, and malperfusion syndromes. The separation of the intima from aortic walls will form a second blood-filled lumen defined as "false lumen (FL)", where the thrombus is more likely to form due to the local stasis hemodynamic conditions. Complete thrombosis of FL is associated with a beneficial outcome while patency and partial thrombosis will lead to later complications. However, the thrombosis mechanism is still unclear and little is known about the impact of chemical species transported by blood flow on this process. The proteins involved in the coagulation cascade (CC) may play an important role in the process of thrombosis, especially in the activation and stabilization of platelets. Based on this hypothesis, a reduced-order fluid-chemical model was established to simulate CC in an aortic dissection phantom with two tears. A high level of fibrin is continuously observed at the top of the FL and some time-varying areas between two tears, indicating a high likelihood of thrombus formation there. This finding is consistent with the clinical observation. The time evolution of coagulation factors is greatly affected by local hemodynamics, especially in the high disturbance zone where the evolution has characteristics of periodic changes consistent with the flow field. The ability of the proposed model to reproduce the CC response provides a potential application to integrate with a model that can simulate platelet activities, forming a biochemical-based model which would help unveil the mechanisms of thrombosis in FL and the clinical decision of appropriate treatment.

Modeling Thrombus Shell: Linking Adhesion Receptor Properties and Macroscopic Dynamics

Damage to arterial vessel walls leads to the formation of platelet aggregate, which acts as a physical obstacle for bleeding. An arterial thrombus is heterogeneous; it has a dense inner part (core) and an unstable outer part (shell). The thrombus shell is very dynamic, being composed of loosely connected discoid platelets. The mechanisms underlying the observed mobility of the shell and its (patho)physiological implications are unclear. To investigate arterial thrombus mechanics, we developed a novel, to our knowledge, two-dimensional particle-based computational model of microvessel thrombosis. The model considers two types of interplatelet interactions: primary reversible (glycoprotein Ib (GPIb)-mediated) and stronger integrin-mediated interaction, which intensifies with platelet activation. At high shear rates, the former interaction leads to adhesion, and the latter is primarily responsible for stable platelet aggregation. Using a stochastic model of GPIb-mediated interaction, we initially reproduced experimental curves that characterize individual platelet interactions with a von Willebrand factor-coated surface. The addition of the second stabilizing interaction results in thrombus formation. The comparison of thrombus dynamics with experimental data allowed us to estimate the magnitude of critical interplatelet forces in the thrombus shell and the characteristic time of platelet activation. The model predicts moderate dependence of maximal thrombus height on the injury size in the absence of thrombin activity. We demonstrate that the developed stochastic model reproduces the observed highly dynamic behavior of the thrombus shell. The presence of primary stochastic interaction between platelets leads to the properties of thrombus consistent with in vivo findings; it does not grow upstream of the injury site and covers the whole injury from the first seconds of the formation. А simplified model, in which GPIb-mediated interaction is deterministic, does not reproduce these features. Thus, the stochasticity of platelet interactions is critical for thrombus plasticity, suggesting that interaction via a small number of bonds drives the dynamics of arterial thrombus shell.

Non-invasive Inference of Thrombus Material Properties with Physics-Informed Neural Networks

We employ physics-informed neural networks (PINNs) to infer properties of biological materials using synthetic data. In particular, we successfully apply PINNs on inferring permeability and viscoelastic modulus from thrombus deformation data, which can be described by the fourth-order Cahn-Hilliard and Navier-Stokes Equations. In PINNs, the partial differential equations are encoded into a loss function, where partial derivatives can be obtained through automatic differentiation (AD). In addition to tackling the challenge of calculating the fourth-order derivative in the Cahn-Hilliard equation with AD, we introduce an auxiliary network along with the main neural network to approximate the second-derivative of the energy potential term. Our model can simultaneously predict unknown material parameters and velocity, pressure, and deformation gradient fields by merely training with partial information among all data, i.e., phase field and pressure measurements, while remaining highly flexible in sampling within the spatio-temporal domain for data acquisition. We validate our model by numerical solutions from the spectral/hp element method (SEM) and demonstrate its robustness by training it with noisy measurements. Our results show that PINNs can infer the material properties from noisy synthetic data, and thus they have great potential for inferring these properties from experimental multi-modality and multi-fidelity data.

Experimental and mathematical modelling of magnetically labelled mesenchymal stromal cell delivery

A key challenge for stem cell therapies is the delivery of therapeutic cells to the repair site. Magnetic targeting has been proposed as a platform for defining clinical sites of delivery more effectively. In this paper, we use a combined in vitro experimental and mathematical modelling approach to explore the magnetic targeting of mesenchymal stromal cells (MSCs) labelled with magnetic nanoparticles using an external magnet. This study aims to (i) demonstrate the potential of magnetic tagging for MSC delivery, (ii) examine the effect of red blood cells (RBCs) on MSC capture efficacy and (iii) highlight how mathematical models can provide both insight into mechanics of therapy and predictions about cell targeting in vivo. In vitro MSCs are cultured with magnetic nanoparticles and circulated with RBCs over an external magnet. Cell capture efficacy is measured for varying magnetic field strengths and RBC percentages. We use a 2D continuum mathematical model to represent the flow of magnetically tagged MSCs with RBCs. Numerical simulations demonstrate qualitative agreement with experimental results showing better capture with stronger magnetic fields and lower levels of RBCs. We additionally exploit the mathematical model to make hypotheses about the role of extravasation and identify future in vitro experiments to quantify this effect.

Integrating blood cell mechanics, platelet adhesive dynamics and coagulation cascade for modelling thrombus formation in normal and diabetic blood

Normal haemostasis is an important physiological mechanism that prevents excessive bleeding during trauma, whereas the pathological thrombosis especially in diabetics leads to increased incidence of heart attacks and strokes as well as peripheral vascular events. In this work, we propose a new multiscale framework that integrates seamlessly four key components of blood clotting, namely transport of coagulation factors, coagulation kinetics, blood cell mechanics and platelet adhesive dynamics, to model the development of thrombi under physiological and pathological conditions. We implement this framework to simulate platelet adhesion due to the exposure of tissue factor in a three-dimensional microchannel. Our results show that our model can simulate thrombin-mediated platelet activation in the flowing blood, resulting in platelet adhesion to the injury site of the channel wall. Furthermore, we simulate platelet adhesion in diabetic blood, and our results show that both the pathological alterations in the biomechanics of blood cells and changes in the amount of coagulation factors contribute to the excessive platelet adhesion and aggregation in diabetic blood. Taken together, this new framework can be used to probe synergistic mechanisms of thrombus formation under physiological and pathological conditions, and open new directions in modelling complex biological problems that involve several multiscale processes.

Computational investigation of blood flow and flow-mediated transport in arterial thrombus neighborhood

A pathologically formed blood clot or thrombus is central to major cardiovascular diseases like heart attack and stroke. Detailed quantitative evaluation of flow and flow-mediated transport processes in the thrombus neighborhood within large artery hemodynamics is crucial for understanding disease progression and assessing treatment efficacy. This, however, remains a challenging task owing to the complexity of pulsatile viscous flow interactions with arbitrary shape and heterogeneous microstructure of realistic thrombi. Here, we address this challenge by conducting a systematic parametric simulation-based study on characterizing unsteady hemodynamics and flow-mediated transport in the neighborhood of an arterial thrombus. We use a hybrid particle-continuum-based finite element approach to handle arbitrary thrombus shape and microstructural variations. Results from a cohort of 50 different unsteady flow scenarios are presented, including unsteady vortical structures, pressure gradient across the thrombus boundary, finite time Lyapunov exponents, and dynamic coherent structures that organize advective transport. We clearly illustrate the combined influence of three key parameters-thrombus shape, microstructure, and extent of wall disease-in terms of: (a) determining hemodynamic features in the thrombus neighborhood and (b) governing the balance between advection, permeation, and diffusion to regulate transport processes in the thrombus neighborhood.

Modelling the linkage between influenza infection and cardiovascular events via thrombosis

There is a heavy burden associated with influenza including all-cause hospitalization as well as severe cardiovascular and cardiorespiratory events. Influenza associated cardiac events have been linked to multiple biological pathways in a human host. To study the contribution of influenza virus infection to cardiovascular thrombotic events, we develop a dynamic model which incorporates some key elements of the host immune response, inflammatory response, and blood coagulation. We formulate these biological systems and integrate them into a cohesive modelling framework to show how blood clotting may be connected to influenza virus infection. With blood clot formation inside an artery resulting from influenza virus infection as the primary outcome of this integrated model, we demonstrate how blood clot severity may depend on circulating prothrombin levels. We also utilize our model to leverage clinical data to inform the threshold level of the inflammatory cytokine TNFα which initiates tissue factor induction and subsequent blood clotting. Our model provides a tool to explore how individual biological components contribute to blood clotting events in the presence of influenza infection, to identify individuals at risk of clotting based on their circulating prothrombin levels, and to guide the development of future vaccines to optimally interact with the immune system.

Predictive modelling of thrombus formation in diabetic retinal microaneurysms

Microaneurysms (MAs) are one of the earliest clinically visible signs of diabetic retinopathy (DR). Vision can be reduced at any stage of DR by MAs, which may enlarge, rupture and leak fluid into the neural retina. Recent advances in ophthalmic imaging techniques enable reconstruction of the geometries of MAs and quantification of the corresponding haemodynamic metrics, such as shear rate and wall shear stress, but there is lack of computational models that can predict thrombus formation in individual MAs. In this study, we couple a particle model to a continuum model to simulate the platelet aggregation in MAs with different shapes. Our simulation results show that under a physiologically relevant blood flow rate, thrombosis is more pronounced in saccular-shaped MAs than fusiform-shaped MAs, in agreement with recent clinical findings. Our model predictions of the size and shape of the thrombi in MAs are consistent with experimental observations, suggesting that our model is capable of predicting the formation of thrombus for newly detected MAs. This is the first quantitative study of thrombosis in MAs through simulating platelet aggregation, and our results suggest that computational models can be used to predict initiation and development of intraluminal thrombus in MAs as well as provide insights into their role in the pathophysiology of DR.

Modeling Thrombin Generation in Plasma under Diffusion and Flow

We investigate the capacity of published numerical models of thrombin generation to reproduce experimentally observed threshold behavior under conditions in which diffusion and/or flow are important. Computational fluid dynamics simulations incorporating species diffusion, fluid flow, and biochemical reactions are compared with published data for thrombin generation in vitro in 1) quiescent plasma exposed to patches of tissue factor and 2) plasma perfused through a capillary coated with tissue factor. Clot time is correctly predicted in individual cases, and some models qualitatively replicate thrombin generation thresholds across a series of tissue factor patch sizes or wall shear rates. Numerical results suggest that there is not a genuine patch size threshold in quiescent plasma-clotting always occurs given enough time-whereas the shear rate threshold observed under flow is a genuine physical limit imposed by flow-mediated washout of active coagulation factors. Despite the encouraging qualitative results obtained with some models, no single model robustly reproduces all experiments, demonstrating that greater understanding of the underlying reaction network, and particularly of surface reactions, is required. In this direction, additional simulations provide evidence that 1) a surface-localized enzyme, speculatively identified as meizothrombin, is significantly active toward the fluorescent thrombin substrate used in the experiments or, less likely, 2) thrombin is irreversibly inhibited at a faster-than-expected rate, possibly explained by a stimulatory effect of plasma heparin on antithrombin. These results highlight the power of simulation to provide novel mechanistic insights that augment experimental studies and build our understanding of complex biophysicochemical processes. Further validation work is critical to unleashing the full potential of coagulation models as tools for drug development and personalized medicine.

A three-dimensional phase-field model for multiscale modeling of thrombus biomechanics in blood vessels

Mechanical interactions between flowing and coagulated blood (thrombus) are crucial in dictating the deformation and remodeling of a thrombus after its formation in hemostasis. We propose a fully-Eulerian, three-dimensional, phase-field model of thrombus that is calibrated with existing in vitro experimental data. This phase-field model considers spatial variations in permeability and material properties within a single unified mathematical framework derived from an energy perspective, thereby allowing us to study effects of thrombus microstructure and properties on its deformation and possible release of emboli under different hemodynamic conditions. Moreover, we combine this proposed thrombus model with a particle-based model which simulates the initiation of the thrombus. The volume fraction of a thrombus obtained from the particle simulation is mapped to an input variable in the proposed phase-field thrombus model. The present work is thus the first computational study to integrate the initiation of a thrombus through platelet aggregation with its subsequent viscoelastic responses to various shear flows. This framework can be informed by clinical data and potentially be used to predict the risk of diverse thromboembolic events under physiological and pathological conditions.

Simulation of thrombosis in a stenotic microchannel: The effects of vWF-enhanced shear activation of platelets

This study was undertaken to develop a numerical/computational simulation of von Willebrand Factor (vWF) - mediated platelet shear activation and deposition in an idealized stenosis. Blood is treated as a multi-constituent mixture comprised of a linear fluid component and a porous solid component (thrombus). Chemical and biological species involved in coagulation are modeled using a system of coupled convection-reaction-diffusion (CRD) equations. This study considers the cumulative effect of shear stress (history) on platelet activation. The vWF activity is modeled as an enhancement function for the shear stress accumulation and is related to the experimentally-observed unfolding rate of vWF. A series of simulations were performed in an idealized stenosis in which the predicted platelets deposition agreed well with previous experimental observations spatially and temporally, including the reduction of platelet deposition with decreasing expansion angle. Further simulation indicated a direct relationship between vWF-mediated platelet deposition and degree of stenosis. Based on the success with these benchmark simulations, it is hoped that the model presented here may provide additional insight into vWF-mediated thrombosis and prove useful for the development of more hemo-compatible blood-wetted devices in the future.

Recent advances in computational modeling of fibrin clot formation: A review

The field of thrombosis and hemostasis is crucial for understanding and developing new therapies for pathologies such as deep vein thrombosis, diabetes related strokes, pulmonary embolisms, and hemorrhaging related diseases. In the last two decades, an exponential growth in studies related to fibrin clot formation using computational tools has been observed. Despite this growth, the complete mechanism behind thrombus formation and hemostasis has been long and rife with obstacles; however, significant progress has been made in the present century. The computational models and methods used in this context are diversified into different spatiotemporal scales, yet there is no single model which can predict both physiological and mechanical properties of fibrin clots. In this review, we list the major strategies employed by researchers in modeling fibrin clot formation using recent and existing computational techniques. This review organizes the computational strategies into continuum level, system level, discrete particle (DPD), and multi-scale methods. We also discuss strengths and weaknesses of various methods and future directions in which computational modeling of fibrin clots can advance.

A density-dependent FEM-FCT algorithm with application to modeling platelet aggregation

Upon injury to a blood vessel, flowing platelets will aggregate at the injury site, forming a plug to prevent blood loss. Through a complex system of biochemical reactions, a stabilizing mesh forms around the platelet aggregate forming a blood clot that further seals the injury. Computational models of clot formation have been developed to a study intravascular thrombosis, where a vessel injury does not cause blood leakage outside the blood vessel but blocks blood flow. To model scenarios in which blood leaks from a main vessel out into the extravascular space, new computational tools need to be developed to handle the complex geometries that represent the injury. We have previously modeled intravascular clot formation under flow using a continuum approach wherein the transport of platelet densities into some spatial location is limited by the platelet fraction that already reside in that location, i.e., the densities satisfy a maximum packing constraint through the use of a hindered transport coefficient. To extend this notion to extravascular injury geometries, we have modified a finite element method flux-corrected transport (FEM-FCT) scheme by prelimiting antidiffusive nodal fluxes. We show that our modified scheme, under a variety of test problems, including mesh refinement, structured vs unstructured meshes, and for a range of reaction rates, produces numerical results that satisfy a maximum platelet-density packing constraint.

Reduced model to predict thrombin and fibrin during thrombosis on collagen/tissue factor under venous flow: Roles of γ'-fibrin and factor XIa

During thrombosis, thrombin generates fibrin, however fibrin reversibly binds thrombin with low affinity E-domain sites (K_D = 2.8 μM) and high affinity γ’-fibrin sites (K_D = 0.1 μM). For blood clotting on collagen/tissue factor (1 TF-molecule/μm²) at 200 s^-1 wall shear rate, high μM-levels of intraclot thrombin suggest robust prothrombin penetration into clots. Setting intraclot zymogen concentrations to plasma levels (and neglecting cofactor rate limitations) allowed the linearization of 7 Michaelis-Menton reactions between 6 species to simulate intraclot generation of: Factors FXa (via TF/VIIa or FIXa), FIXa (via TF/FVIIa or FXIa), thrombin, fibrin, and FXIa. This reduced model [7 rates, 2 K_D’s, enzyme half-lives~1 min] predicted the measured clot elution rate of thrombin-antithrombin (TAT) and fragment F1.2 in the presence and absence of the fibrin inhibitor Gly-Pro-Arg-Pro. To predict intraclot fibrin reaching 30 mg/mL by 15 min, the model required fibrinogen penetration into the clot to be strongly diffusion-limited (actual rate/ideal rate = 0.05). The model required free thrombin in the clot (~100 nM) to have an elution half-life of ~2 sec, consistent with measured albumin elution, with most thrombin (>99%) being fibrin-bound. Thrombin-feedback activation of FXIa became prominent and reached 5 pM FXIa at >500 sec in the simulation, consistent with anti-FXIa experiments. In predicting intrathrombus thrombin and fibrin during 15-min microfluidic experiments, the model revealed “cascade amplification” from 30 pM levels of intrinsic tenase to 15 nM prothrombinase to 15 μM thrombin to 90 μM fibrin. Especially useful for multiscale simulation, this reduced model predicts thrombin and fibrin co-regulation during thrombosis under flow.

During blood clotting events, a complex series of reaction are involved. Simulation gives insights to the concentration of different enzymes which are at too low of concentration to be detected. However, the models are often large and difficult to solve for clotting under flow conditions. With a thin film approximation, we were able to simplify clotting under flow with parameters from literature, with only 3 adjusted in order to fit the experimental data. This model gave insights into the dynamics of the species involved, and the roles of γ’-fibrin and thrombin feedback activation. This reduced model may be useful in further multiscale simulations.

Utilizing Plasma Composition Data to Help Determine Procoagulant Dynamics in Patients with Thermal Injury: A Computational Assessment

INTRODUCTION

The development of methods that generate individualized assessments of the procoagulant potential of burn patients could improve their treatment. Beyond its role as an essential intermediate in the formation of thrombin, factor (F)Xa has systemic effects as an agonist to inflammatory processes. In this study, we use a computational model to study the FXa dynamics underlying tissue factor-initiated thrombin generation in a small cohort of burn patients.

MATERIALS AND METHODS

Plasma samples were collected upon admission (Hour 0) from nine subjects (five non-survivors) with major burn injuries and then at 48 hours. Coagulation factor concentrations (II, V, VII, VIII, IX, X, TFPI, antithrombin (AT), protein C (PC)) were measured and used in a computational model to generate time course profiles for thrombin (IIa), FXa, extrinsic tenase, intrinsic tenase and prothrombinase complexes upon a 5 pM tissue factor stimulus in the presence of 1 nM thrombomodulin. Parameters were extracted from the thrombin and FXa profiles (including max rate (MaxRIIa and MaxRFXa) and peak level (MaxLIIa and MaxLFXa)). Procoagulant potential was also evaluated by determining the concentration of the complexes at select times. Parameter values were compared between survivors and non-survivors in the burn cohort and between the burn cohort and a simulation based on the mean physiological (100%) concentration for all factor levels.

RESULTS

Burn patients differed at Hour 0 (p < 0.05) from 100% mean physiological levels for all coagulation factor levels except FV and FVII. The concentration of FX, FII, TFPI, AT and PC was lower; FIX and FVIII were increased. The composition differences resulted in all nine burn patients at Hour 0 displaying a procoagulant phenotype relative to 100% mean physiological simulation (MaxLIIa (306 ± 90 nM vs. 52 nM), MaxRIIa (2.9 ± 1.1 nM/s vs. 0.3 nM/s), respectively p < 0.001); MaxRFXa and MaxLFXa were also an order of magnitude greater than 100% mean physiological simulation (p < 0.001). When grouped by survival status and compared at the time of admission, non-survivors had lower PC levels (56 ± 18% vs. 82 ± 9%, p < 0.05), and faster MaxRFXa (29 ± 6 pM/s vs. 18 ± 6 pM/s, p < 0.05) than those that survived; similar trends were observed for all other procoagulant parameters. At 48 hours when comparing non-survivors to survivors, TFPI levels were higher (108 ± 18% vs. 59 ± 18%, p < 0.05), and MaxRIIa (1.5 ± 1.4 nM/s vs. 3.6 ± 0.7 nM/s, p < 0.05) and MaxRFXa (13 ± 12 pM/s vs. 35 ± 4 pM/s, p < 0.05) were lower; similar trends were observed with all other procoagulant parameters. Overall, between admission and 48 hours, procoagulant potential, as represented by MaxR and MaxL parameters for thrombin and FXa, in non-survivors decreased while in survivors they increased (p < 0.05). In patients that survived, there was a positive correlation between FX levels and MaxLFXa (r = 0.96) and reversed in mortality (r= -0.91).

CONCLUSIONS

Thrombin and FXa generation are increased in burn patients at admission compared to mean physiological simulations. Over the first 48 hours, burn survivors became more procoagulant while non-survivors became less procoagulant. Differences between survivors and non-survivors appear to be present in the underlying dynamics that contribute to FXa dynamics. Understanding how the individual specific balance of procoagulant and anticoagulant proteins contributes to thrombin and FXa generation could ultimately guide therapy and potentially reduce burn injury-related morbidity and mortality.

Coagulopathy implications using a multiscale model of traumatic bleeding matching macro- and microcirculation

Quantifying the relationship between vascular injury and the dynamic bleeding rate requires a multiscale model that accounts for changing and coupled hemodynamics between the global and microvascular levels. A lumped, global hemodynamic model of the human cardiovascular system with baroreflex control was coupled to a local 24-level bifurcating vascular network that spanned diameters from the muscular artery scale (0.1-1.3 mm) to capillaries (5-10 μm) via conservation of momentum and conservation of mass boundary conditions. For defined injuries of severing all vessels at each nth-level, the changing pressures and flowrates were calculated using prescribed shear-dependent hemostatic clot growth rates (normal or coagulopathic). Key results were as follows: 1) the upstream vascular network rapidly depressurizes to reduce blood loss; 2) wall shear rates at the hemorrhaging wound exit are sufficiently high (~10,000 s^-1) to drive von Willebrand factor unfolding; 3) full coagulopathy results in >2-liter blood loss in 2 h for severing all vessels of 0.13- to 0.005-mm diameter within the bifurcating network, whereas full hemostasis limits blood loss to <100 ml within 2 min; and 4) hemodilution from transcapillary refill increases blood loss and could be implicated in trauma-induced coagulopathy. A sensitivity analysis on length-to-diameter ratio and branching exponent demonstrated that bleeding was strongly dependent on these tissue-dependent network parameters. This is the first bleeding model that prescribes the geometry of the injury to calculate the rate of pressure-driven blood loss and local wall shear rate in the presence or absence of coagulopathic blood. NEW & NOTEWORTHY We developed a multiscale model that couples a lumped, global hemodynamic model of a patient to resolved, single-vessel wounds ranging from the small artery to capillary scale. The model is able to quantify wall shear rates, seal rates, and blood loss rates in the presence and absence of baroreflex control and hemodilution.

In vivo measurement of blood clot mechanics from computational fluid dynamics based on intravital microscopy images

Ischemia which leads to heart attacks and strokes is one of the major causes of death in the world. Whether an occlusion occurs or not depends on the ability of a growing thrombus to resist flow forces exerted on its structure. This manuscript provides the first known in vivo measurement of how much stress a clot can withstand, before yielding to the surrounding blood flow. Namely, Lattice-Boltzmann Method flow simulations are performed based on 3D clot geometries, which are estimated from intravital microscopy images of laser-induced injuries in cremaster microvasculature of live mice. In addition to reporting the blood clot yield stresses, we also show that the thrombus "core" does not experience significant deformation, while its "shell" does. This indicates that the shell is more prone to embolization. Therefore, drugs should be designed to target the shell selectively, while leaving the core intact to minimize excessive bleeding. Finally, we laid down a foundation for a nondimensionalization procedure which unraveled a relationship between clot mechanics and biology. Hence, the proposed framework could ultimately lead to a unified theory of thrombogenesis, capable of explaining all clotting events. Thus, the findings presented herein will be beneficial to the understanding and treatment of heart attacks, strokes and hemophilia.

Multiscale simulation of thrombus growth and vessel occlusion triggered by collagen/tissue factor using a data-driven model of combinatorial platelet signalling

During clotting under flow, platelets bind and activate on collagen and release autocrinic factors such as ADP and thromboxane, while tissue factor (TF) on the damaged wall leads to localized thrombin generation. Towards patient-specific simulation of thrombosis, a multiscale approach was developed to account for: platelet signalling [neural network (NN) trained by pairwise agonist scanning (PAS), PAS-NN], platelet positions (lattice kinetic Monte Carlo, LKMC), wall-generated thrombin and platelet-released ADP/thromboxane convection–diffusion (partial differential equation, PDE) and flow over a growing clot (lattice Boltzmann). LKMC included shear-driven platelet aggregate restructuring. The PDEs for thrombin, ADP and thromboxane were solved by finite element method using cell activation-driven adaptive triangular meshing. At all times, intracellular calcium was known for each platelet by PAS-NN in response to its unique exposure to local collagen, ADP, thromboxane and thrombin. When compared with microfluidic experiments of human blood clotting on collagen/TF driven by constant pressure drop, the model accurately predicted clot morphology and growth with time. In experiments and simulations at TF at 0.1 and 10 molecule-TF/\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\mu$\end{document}m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$^{2}$\end{document} and initial wall shear rate of 200 s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$^{-1}$\end{document}, the occlusive blockade of flow for a 60-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\mu$\end{document}m channel occurred relatively abruptly at 600 and 400 s, respectively (with no occlusion at zero TF). Prior to occlusion, intrathrombus concentrations reached 50 nM thrombin, ~ 1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\mu$\end{document}M thromboxane and ~ 10 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\mu$\end{document}M ADP, while the wall shear rate on the rough clot peaked at ~ 1000–2000 s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$^{-1}$\end{document}. Additionally, clotting on TF/collagen was accurately simulated for modulators of platelet cyclooxygenase-1, P2Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{1}$\end{document} and IP-receptor. This multiscale approach facilitates patient-specific simulation of thrombosis under hemodynamic and pharmacological conditions.

Three-phase Model of Visco-elastic Incompressible Fluid Flow and its Computational Implementation

Energetic Variational Approach is used to derive a novel thermodynamically consistent three-phase model of a mixture of Newtonian and visco-elastic fluids. The model which automatically satisfies the energy dissipation law and is Galilean invariant, consists of coupled Navier-Stokes and Cahn-Hilliard equations. Modified General Navier Boundary Condition with fluid elasticity taken into account is also introduced for using the model to study moving contact line problems. Energy stable numerical scheme is developed to solve system of model equations efficiently. Convergence of the numerical scheme is verified by simulating a droplet sliding on an inclined plane under gravity. The model can be applied for studying various biological or biophysical problems. Predictive abilities of the model are demonstrated by simulating deformation of venous blood clots with different visco-elastic properties and experimentally observed internal structures under different biologically relevant shear blood flow conditions.

Model predictions of deformation, embolization and permeability of partially obstructive blood clots under variable shear flow

Thromboembolism, one of the leading causes of morbidity and mortality worldwide, is characterized by formation of obstructive intravascular clots (thrombi) and their mechanical breakage (embolization). A novel two-dimensional multi-phase computational model is introduced that describes active interactions between the main components of the clot, including platelets and fibrin, to study the impact of various physiologically relevant blood shear flow conditions on deformation and embolization of a partially obstructive clot with variable permeability. Simulations provide new insights into mechanisms underlying clot stability and embolization that cannot be studied experimentally at this time. In particular, model simulations, calibrated using experimental intravital imaging of an established arteriolar clot, show that flow-induced changes in size, shape and internal structure of the clot are largely determined by two shear-dependent mechanisms: reversible attachment of platelets to the exterior of the clot and removal of large clot pieces. Model simulations predict that blood clots with higher permeability are more prone to embolization with enhanced disintegration under increasing shear rate. In contrast, less permeable clots are more resistant to rupture due to shear rate-dependent clot stiffening originating from enhanced platelet adhesion and aggregation. These results can be used in future to predict risk of thromboembolism based on the data about composition, permeability and deformability of a clot under specific local haemodynamic conditions.

Structural modelling of the cardiovascular system

Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.

Computational Fluid Dynamics Assessment Associated with Transcatheter Heart Valve Prostheses: A Position Paper of the ISO Working Group

The governing international standard for the development of prosthetic heart valves is International Organization for Standardization (ISO) 5840. This standard requires the assessment of the thrombus potential of transcatheter heart valve substitutes using an integrated thrombus evaluation. Besides experimental flow field assessment and ex vivo flow testing, computational fluid dynamics is a critical component of this integrated approach. This position paper is intended to provide and discuss best practices for the setup of a computational model, numerical solving, post-processing, data evaluation and reporting, as it relates to transcatheter heart valve substitutes. This paper is not intended to be a review of current computational technology; instead, it represents the position of the ISO working group consisting of experts from academia and industry with regards to considerations for computational fluid dynamic assessment of transcatheter heart valve substitutes.

Thrombosis in Cerebral Aneurysms and the Computational Modeling Thereof: A Review

Thrombosis is a condition closely related to cerebral aneurysms and controlled thrombosis is the main purpose of endovascular embolization treatment. The mechanisms governing thrombus initiation and evolution in cerebral aneurysms have not been fully elucidated and this presents challenges for interventional planning. Significant effort has been directed towards developing computational methods aimed at streamlining the interventional planning process for unruptured cerebral aneurysm treatment. Included in these methods are computational models of thrombus development following endovascular device placement. The main challenge with developing computational models for thrombosis in disease cases is that there exists a wide body of literature that addresses various aspects of the clotting process, but it may not be obvious what information is of direct consequence for what modeling purpose (e.g., for understanding the effect of endovascular therapies). The aim of this review is to present the information so it will be of benefit to the community attempting to model cerebral aneurysm thrombosis for interventional planning purposes, in a simplified yet appropriate manner. The paper begins by explaining current understanding of physiological coagulation and highlights the documented distinctions between the physiological process and cerebral aneurysm thrombosis. Clinical observations of thrombosis following endovascular device placement are then presented. This is followed by a section detailing the demands placed on computational models developed for interventional planning. Finally, existing computational models of thrombosis are presented. This last section begins with description and discussion of physiological computational clotting models, as they are of immense value in understanding how to construct a general computational model of clotting. This is then followed by a review of computational models of clotting in cerebral aneurysms, specifically. Even though some progress has been made towards computational predictions of thrombosis following device placement in cerebral aneurysms, many gaps still remain. Answering the key questions will require the combined efforts of the clinical, experimental and computational communities.

Multiscale systems biology of trauma-induced coagulopathy

Trauma with hypovolemic shock is an extreme pathological state that challenges the body to maintain blood pressure and oxygenation in the face of hemorrhagic blood loss. In conjunction with surgical actions and transfusion therapy, survival requires the patient's blood to maintain hemostasis to stop bleeding. The physics of the problem are multiscale: (a) the systemic circulation sets the global blood pressure in response to blood loss and resuscitation therapy, (b) local tissue perfusion is altered by localized vasoregulatory mechanisms and bleeding, and (c) altered blood and vessel biology resulting from the trauma as well as local hemodynamics control the assembly of clotting components at the site of injury. Building upon ongoing modeling efforts to simulate arterial or venous thrombosis in a diseased vasculature, computer simulation of trauma-induced coagulopathy is an emerging approach to understand patient risk and predict response. Despite uncertainties in quantifying the patient's dynamic injury burden, multiscale systems biology may help link blood biochemistry at the molecular level to multiorgan responses in the bleeding patient. As an important goal of systems modeling, establishing early metrics of a patient's high-dimensional trajectory may help guide transfusion therapy or warn of subsequent later stage bleeding or thrombotic risks. This article is categorized under: Analytical and Computational Methods > Computational Methods Biological Mechanisms > Regulatory Biology Models of Systems Properties and Processes > Mechanistic Models.

Mathematical Techniques for Understanding Platelet Regulation and the Development of New Pharmacological Approaches

Mathematical and computational modeling is currently in the process of becoming an accepted tool in the arsenal of methods utilized for the investigation of complex biological systems. For some problems in the field, like cellular metabolic regulation, neural impulse propagation, or cell cycle, progress is already unthinkable without use of such methods. Mathematical models of platelet signaling, function, and metabolism during the last years have not only been steadily increasing in their number, but have also been providing more in-depth insights, generating hypotheses, and allowing predictions to be made leading to new experimental designs and data. Here we describe the basic approaches to platelet mathematical model development and validation, highlighting the challenges involved. We then review the current theoretical models in the literature and how these are being utilized to increase our understanding of these complex cells.

Dynamic cellular finite-element method for modelling large-scale cell migration and proliferation under the control of mechanical and biochemical cues: a study of re-epithelialization

Computational modelling of cells can reveal insight into the mechanisms of the important processes of tissue development. However, current cell models have limitations and are challenged to model detailed changes in cellular shapes and physical mechanics when thousands of migrating and interacting cells need to be modelled. Here we describe a novel dynamic cellular finite-element model (DyCelFEM), which accounts for changes in cellular shapes and mechanics. It also models the full range of cell motion, from movements of individual cells to collective cell migrations. The transmission of mechanical forces regulated by intercellular adhesions and their ruptures are also accounted for. Intra-cellular protein signalling networks controlling cell behaviours are embedded in individual cells. We employ DyCelFEM to examine specific effects of biochemical and mechanical cues in regulating cell migration and proliferation, and in controlling tissue patterning using a simplified re-epithelialization model of wound tissue. Our results suggest that biochemical cues are better at guiding cell migration with improved directionality and persistence, while mechanical cues are better at coordinating collective cell migration. Overall, DyCelFEM can be used to study developmental processes when a large population of migrating cells under mechanical and biochemical controls experience complex changes in cell shapes and mechanics.

Modelling of thrombin generation under flow in realistic left anterior descending geometries

Currently there are no available methods for prediction of thrombotic complications in Coronary Artery disease. Additionally, blood coagulation tests are mainly performed in a steady system while coagulation in vivo occurs under flow conditions. In this work, a phenomenological model for coagulation up-to thrombin generation is proposed; the model is mainly based on the results of thrombin generation assays and therefore it can account for the variation of the coagulability that is observed in different individuals. The model is applied on 3 cases of left anterior descending arteries (LAD) with 50% maximum stenosis placed at a different location and have been statistically assessed as of different complication risk. The simulations showed that parameters of thrombin generation assays obtain different values when they refer to thrombin generation under realistic coronary flow conditions. The flow conditions prevailing locally because of the geometric differences among the arterial trees can lead to different initiation times and thrombin production rates and it also alters the spatial distribution of the coagulation products. Similarly, small changes of the coagulation characteristics of blood under identical flow conditions can allow or prevent the initiation of coagulation. The results indicate that combined consideration of geometry and coagulation characteristics of blood can lead to entirely different conclusions compared to independent assessment of each factor.