Fluid Mechanics of Blood Clot Formation

Numerical analysis of blood flow through stenosed microvessels using a multi-phase model

Blood flow in arterioles have attracted considerable research attention due to their clinical implications. However, the fluid structure interaction between red blood cells and plasma in the blood poses formidable difficulty to the computational efforts. In this contribution, we seek to represent the red blood cells in the blood as a continuous non-Newtonian phase and construct a multi-phase model for the blood flow in microvessels. The methods are presented and validated using a channel with sudden expansion. And the resulting blood flow inside a stenosed microvessel is investigated at different inlet velocity amplitudes and hematocrits. It is show that the increase of both inlet velocity amplitude and inlet hematocrit leads to longer and thicker cell-rich layer downstream the stenosis. Besides, it is found that the maximum values of wall shear stress scales up with inlet velocity amplitudes and hematocrits. These results show the validity of the proposed computational model and provide helpful insights into blood flow behaviors inside stenosed vessels.

Numerical methods for hemolysis and thrombus evaluation in the percutaneous ventricular assist device

BACKGROUND

A percutaneous ventricular assist device (pVAD) is an effective method to treat heart failure, but its complications, mainly hemolysis and thrombus formation, cannot be ignored. Accurate evaluation of hemolysis and thrombus formation in pVAD is essential to guide the development of pVAD and reduce the incidence of complications.

METHODS

This study optimized the numerical model to predict hemolysis and thrombus formation in pVAD. The hemolysis model is based on the power law function, and the multi-component thrombus prediction model is improved by introducing the von Willebrand factor.

RESULTS

The error between the numerical simulation and the hydraulic performance experiment is within 5%. The numerical results of hemolysis are in good agreement with those of in vitro experiments. Meanwhile, the thrombus location predicted by the numerical model is the same as that found in the in vivo experiment.

CONCLUSION

The numerical model suggested in this study may therefore accurately assess the possible hemolytic and thrombotic dangers in pVAD, making it an effective tool to support the development of pVAD.

In situ polymer gelation in confined flow controls intermittent dynamics

Polymer flows through pores, nozzles and other small channels govern engineered and naturally occurring dynamics in many processes, from 3D printing to oil recovery in the earth's subsurface to a wide variety of biological flows. The crosslinking of polymers can change their material properties dramatically, and it is advantageous to know a priori whether or not crosslinking polymers will lead to clogged channels or cessation of flow. In this study, we investigate the flow of a common biopolymer, alginate, while it undergoes crosslinking by the addition of a crosslinker, calcium, driven through a microfluidic channel at constant flow rate. We map the boundaries defining complete clogging and flow as a function of flow rate, polymer concentration, and crosslinker concentration. Interestingly, the boundaries of the dynamic behavior qualitatively match the thermodynamic jamming phase diagram of attractive colloidal particles. That is, polymer clogging occurs in a region analogous to colloids in a jammed state, while the polymer flows in regions corresponding to colloids in a liquid phase. However, between the dynamic regimes of complete clogging and unrestricted flow, we observe a remarkable phenomenon in which the crosslinked polymer intermittently clogs the channel. This pattern of deposition and removal of a crosslinked gel is simultaneously highly reproducible, long-lasting, and controllable by system parameters. Higher concentrations of polymer and cross-linker result in more frequent ablation, while gels formed at lower component concentrations ablate less frequently. Upon ablation, the eluted gel maintains its shape, resulting in micro-rods several hundred microns long. Our results suggest both rich dynamics of intermittent flows in crosslinking polymers and the ability to control them.

A computational investigation of occlusive arterial thrombosis

The generation of occlusive thrombi in stenotic arteries involves the rapid deposition of millions of circulating platelets under high shear flow. The process is mediated by the formation of molecular bonds of several distinct types between platelets; the bonds capture the moving platelets and stabilize the growing thrombi under flow. We investigated the mechanisms behind occlusive thrombosis in arteries with a two-phase continuum model. The model explicitly tracks the formation and rupture of the two types of interplatelet bonds, the rates of which are coupled with the local flow conditions. The motion of platelets in the thrombi results from competition between the viscoelastic forces generated by the interplatelet bonds and the fluid drag. Our simulation results indicate that stable occlusive thrombi form only under specific combinations for the ranges of model parameters such as rates of bond formation and rupture, platelet activation time, and number of bonds required for platelet attachment.

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

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

Surface modification of polyvinyl chloride with sodium alginate/carboxymethyl chitosan and heparin for realizing the anticoagulation

Thrombosis of extracorporeal circuits causes significant morbidity and mortality worldwide. In this study, plasma treatment technology and chemical grafting method were used to construct heparinized surfaces on the PVC substrate, which could not only reduce thrombosis but also decrease the side effects of the direct injection of anticoagulants. The PVC substrate was modified by plasma treatment technology firstly to obtain the active surface with the hydroxyl groups used for grafting. Then, heparin was grafted onto the modified PVC surface using different grafting strategies to prepare different heparinized surfaces. The experimental results indicated that the sodium alginate (SA) and carboxymethyl chitosan (CCS) used as interlayers could significantly increase the graft density of heparin to improve the anticoagulant effects and hemocompatibility of heparinized surfaces. In addition, the modification of heparin can further improve the anticoagulant effects. The CCS/low-molecular-weight heparin (LWMH) surface has the best anticoagulant properties, which can prolong the activated partial thromboplastin time (APTT) values of human plasma for about 35 s, reduce the hemolysis rates to <0.3 %, and perform well in the in-vitro blood circulation test. The heparinized surfaces prepared in this work have great application potential in anticoagulant treatment for medical devices.

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.

Shear Stress and Endothelial Mechanotransduction in Trauma Patients with Hemorrhagic Shock: Hidden Coagulopathy Pathways and Novel Therapeutic Strategies

Massive trauma remains a leading cause of death and a global public health burden. Post-traumatic coagulopathy may be present even before the onset of resuscitation, and correlates with severity of trauma. Several mechanisms have been proposed to explain the development of abnormal coagulation processes, but the heterogeneity in injuries and patient profiles makes it difficult to define a dominant mechanism. Regardless of the pattern of death, a significant role in the pathophysiology and pathogenesis of coagulopathy may be attributed to the exposure of endothelial cells to abnormal physical forces and mechanical stimuli in their local environment. In these conditions, the cellular responses are translated into biochemical signals that induce/aggravate oxidative stress, inflammation, and coagulopathy. Microvascular shear stress-induced alterations could be treated or prevented by the development and use of innovative pharmacologic strategies that effectively target shear-mediated endothelial dysfunction, including shear-responsive drug delivery systems and novel antioxidants, and by targeting the venous side of the circulation to exploit the beneficial antithrombogenic profile of venous endothelial cells.

Computational mesoscale framework for biological clustering and fractal aggregation

Hierarchical clustering due to diffusion and reaction is a widespread occurrence in natural phenomena, displaying fractal behavior with non-integer size scaling. The study of this phenomenon has garnered interest in both biological systems such as morphogenesis and blood clotting, and synthetic systems such as colloids and polymers. The modeling of biological clustering can be difficult, as it can occur on a variety of scales and involve multiple mechanisms, necessitating the use of various methods to capture its behavior. Here, we propose a novel framework, the generalized-mesoscale-clustering (GMC), for the study of complex hierarchical clustering phenomena in biological systems. The GMC framework incorporates the effects of hydrodynamic interactions, bonding, and surface tension, and allows for the analysis of both static and dynamic states of cluster development. The framework is applied to a range of biological clustering mechanisms, with a focus on blood-related clustering from fibrin network formation to platelet aggregation. Our study highlights the importance of a comprehensive characterization of the structural properties of the cluster, including fractal dimension, pore-scale diffusion, initiation time, and consolidation time, in fully understanding the behavior of biological clustering systems. The GMC framework also provides the potential to investigate the temporal evolution and mechanical properties of the clusters by tracking bond density and including hydrodynamic interactions.

Advanced in vitro hemocompatibility assessment of biomaterials using a new flow incubation system

Physiologically relevant in vitro hemocompatibility assessment of biomaterials remains challenging. We present a new setup that enables standardized whole blood incubation of biomedical materials under flow. A blood volume of 2 mL is recirculated over test surfaces in a custom-made parallel plate incubation system to determine the activation of hemostasis and inflammation. Controlled physiological shear rates between 125 s-1 and 1250 s-1 and minimized contact to air are combined with a natural-like pumping process. A unique feature of this setup allows tracing adhesion of blood cells to test surfaces microscopically in situ. Validation testing was performed in comparison to previously applied whole blood incubation methodologies. Experiments with the newly developed setup showed that even small obstacles to blood flow activate blood (independent of materials-induced blood activation levels); that adhesion of blood cells to biomaterials equilibrates within 5 to 10 min; that high shear rates (1250 compared to 375 s-1) induce platelet activation; and that hemolysis, platelet factor 4 (PF4) release and platelet loss - but not thrombin formation - depend on shear rate (within the range investigated, 125 to 1250 s-1).

Critical evaluation of kinetic schemes for coagulation

Two well-established numerical representations of the coagulation cascade either initiated by the intrinsic system (Chatterjee et al., PLOS Computational Biology 2010) or the extrinsic system (Butenas et al., Journal of Biological Chemistry, 2004) were compared with thrombin generation assays under realistic pathological conditions. Biochemical modifications such as the omission of reactions not relevant to the case studied, the modification of reactions related to factor XI activation and auto-activation, the adaptation of initial conditions to the thrombin assay system, and the adjustment of some of the model parameters were necessary to align in vitro and in silico data. The modified models are able to reproduce thrombin generation for a range of factor XII, XI, and VIII deficiencies, with the coagulation cascade initiated either extrinsically or intrinsically. The results emphasize that when existing models are extrapolated to experimental parameters for which they have not been calibrated, careful adjustments are required.

Liver Clot after Periodontal Surgery: A Rare Case Report

Due to the primary closure of the fragile tissues, large postoperative drainage is unexpected following a periodontal medical procedure. This case study elaborately illustrates the formation of a "liver coagulation" or "currant jam clot" following a periodontal fold surgical operation. Contamination, internal damage, and the presence of foreign particles such as bone splinters or bits of dental restorative dressing can all increase the risk of an incident happening.

Personalized biomechanical insights in atrial fibrillation: opportunities & challenges

INTRODUCTION

Atrial fibrillation (AF) is an increasingly prevalent and significant worldwide health problem. Manifested as an irregular atrial electrophysiological activation, it is associated with many serious health complications. AF affects the biomechanical function of the heart as contraction follows the electrical activation, subsequently leading to reduced blood flow. The underlying mechanisms behind AF are not fully understood, but it is known that AF is highly correlated with the presence of atrial fibrosis, and with a manifold increase in risk of stroke.

AREAS COVERED

In this review, we focus on biomechanical aspects in atrial fibrillation, current and emerging use of clinical images, and personalized computational models. We also discuss how these can be used to provide patient-specific care.

EXPERT OPINION

Understanding the connection betweenatrial fibrillation and atrial remodeling might lead to valuable understanding of stroke and heart failure pathophysiology. Established and emerging imaging modalities can bring us closer to this understanding, especially with continued advancements in processing accuracy, reproducibility, and clinical relevance of the associated technologies. Computational models of cardiac electromechanics can be used to glean additional insights on the roles of AF and remodeling in heart function.

Near-wall hemodynamic changes in subclavian artery perfusion induced by retrograde inner branched thoracic endograft implantation

Objective

Left subclavian artery (LSA)-branched endografts with retrograde inner branch configuration (thoracic branch endoprosthesis [TBE]) offer a complete endovascular solution when LSA preservation is required during zone 2 thoracic endovascular aortic repair. However, the hemodynamic consequences of the TBE have not been well-investigated. We compared near-wall hemodynamic parameters before and after the TBE implantation using computational fluid dynamic simulations.

Methods

Eleven patients who had undergone TBE implantation were included. Three-dimensional aortic arch geometries were constructed from the pre- and post-TBE implantation computed tomography images. The resulting 22 three-dimensional aortic arch geometries were then discretized into finite element meshes for computational fluid dynamic simulations. Inflow boundary conditions were prescribed using normal physiological pulsatile circulation. Outlet boundary conditions consisted of Windkessel models with previously published values. Blood flow, modeled as Newtonian fluid, simulations were performed with rigid wall assumptions using SimVascular's incompressible Navier-Stokes solver. We compared well-established hemodynamic descriptors: pressure, flow rate, time-averaged wall shear stress (TAWSS), the oscillatory shear index (OSI), and percent area with an OSI of >0.2. Data were presented on the stented portion of the LSA.

Results

TBE implantation was associated with a small decrease in peak LSA pressure (153 mm Hg; interquartile range [IQR], 151-154 mm Hg vs 159 mm Hg; IQR, 158-160 mm Hg; P = .005). No difference was observed in peak LSA flow rates before and after implantation: 40.4 cm³/ (IQR, 39.5-41.6 cm³/s) vs 41.3 cm³/s (IQR, 37.2-44.8 cm³/s; P = .59). There was a significant postimplantation increase in TAWSS (15.2 dynes/cm² [IQR, 12.2-17.7 dynes/cm²] vs 6.2 dynes/cm² [IQR, 5.7-10.3 dynes/cm²]; P = .003), leading to decreases in both the OSI (0.088 [IQR, 0.063 to -0.099] vs 0.1 [IQR, 0.096-0.16]; P = .03) and percentage of area with an OSI of >0.2 (10.4 [IQR, 5.8-15.8] vs 15.7 [IQR, 10.7-31.9]; P = .13). Neither LSA side branch angulation (median, 81°, IQR, 77°-109°) nor moderate compression (16%-58%) seemed to have an impact on the pressure, flow rate, TAWSS, or percentage of area with an OSI of >0.2 in the stented LSA.

Conclusions

The implantation of TBE produces modest hemodynamic disturbances that are unlikely to result in clinically relevant changes.

A Computational Investigation of Occlusive Arterial Thrombosis

The generation of occlusive thrombi in stenotic arteries involves the rapid deposition of millions of circulating platelets under high shear flow. The process is mediated by the formation of molecular bonds of several distinct types between platelets; the bonds capture the moving platelets and stabilize the growing thrombi under flow. We investigated the mechanisms behind occlusive thrombosis in arteries with a two-phase continuum model. The model explicitly tracks the formation and rupture of the two types of interplatelet bonds, the rates of which are coupled with the local flow conditions. The motion of platelets in the thrombi results from competition between the viscoelastic forces generated by the interplatelet bonds and the fluid drag. Our simulation results indicate that stable occlusive thrombi form only under specific combinations for the ranges of model parameters such as rates of bond formation and rupture, platelet activation time, and number of bonds required for platelet attachment.

Wall shear gradient dependent thrombosis studied in blood-on-a-chip with stenotic, branched, and valvular constructions

Thrombosis is the leading cause of death, while the effect of the shear flow on the formation of thrombus in vascular constructions has not been thoroughly understood, and one of the challenges is to observe the origination of thrombus with a controlled flow field. In this work, we use blood-on-a-chip technology to mimic the flow conditions in coronary artery stenosis, neonatal aortic arch, and deep venous valve. The flow field is measured by the microparticle image velocimeter (μPIV). In the experiment, we find that the thrombus often originates at the constructions of stenosis, bifurcation, and the entrance of valve, where the flow stream lines change suddenly, and the maximum wall shear rate gradient appears. Using the blood-on-a-chip technology, the effect of the wall shear rate gradients on the formation of the thrombus has been illustrated, and the blood-on-a-chip is demonstrated to be a perspective tool for further studies on the flow-induced formation of thrombosis.

A model of fluid-structure and biochemical interactions for applications to subclinical leaflet thrombosis

Subclinical leaflet thrombosis (SLT) is a potentially serious complication of aortic valve replacement with a bioprosthetic valve in which blood clots form on the replacement valve. SLT is associated with increased risk of transient ischemic attacks and strokes and can progress to clinical leaflet thrombosis. SLT following aortic valve replacement also may be related to subsequent structural valve deterioration, which can impair the durability of the valve replacement. Because of the difficulty in clinical imaging of SLT, models are needed to determine the mechanisms of SLT and could eventually predict which patients will develop SLT. To this end, we develop methods to simulate leaflet thrombosis that combine fluid-structure interaction and a simplified thrombosis model that allows for deposition along the moving leaflets. Additionally, this model can be adapted to model deposition or absorption along other moving boundaries. We present convergence results and quantify the model's ability to realize changes in valve opening and pressures. These new approaches are an important advancement in our tools for modeling thrombosis because they incorporate both adhesion to the surface of the moving leaflets and feedback to the fluid-structure interaction.

Thrombosis in the pathogenesis of abdominal aortic aneurysm

Background

Abdominal aortic aneurysms (AAAs) are a relatively common vascular pathology of the elderly with high morbidity potential. Irreversible degeneration of the aortic wall leads to lethal rupture if left untreated. Nearly all AAAs contain intraluminal thrombus (ILT) to a varying degree, yet the mechanisms explaining how thrombosis is disturbed in AAA are relatively unknown. This review examined the thrombotic complications associated with AAA, the impact of thrombosis on AAA surgical outcomes and AAA pathogenesis, and the use of antithrombotic therapy in the management of this disease.

Methods

A literature search of the PubMed database was conducted using relevant keywords related to thrombosis and AAAs.

Results

Thrombotic complications are relatively infrequent in AAA yet carry significant morbidity risks. The ILT can impact endovascular aneurysm repair by limiting anatomic suitability and influence the risk of endoleaks. Many of the pathologic mechanisms involved in AAA development, including hemodynamics, inflammation, oxidative stress, and aortic wall remodeling, contain pathways that interact with thrombosis. Conversely, the ILT can also be a source of biochemical stress and exacerbate these aneurysmal processes. In animal AAA models, antithrombotic therapies have shown favorable results in preventing and stabilizing AAA. Antiplatelet agents may be beneficial for reducing risks of major adverse cardiovascular events in AAA patients; however, neither antiplatelet nor anticoagulation is currently used solely for the management of AAA.

Conclusions

Thrombosis and ILT may have detrimental effects on AAA growth, rupture risk, and patient outcomes, yet there is limited understanding of the pathologic thrombotic mechanisms in aneurysmal disease at the molecular level. Preventing ILT using platelet and coagulation inhibitors may be a reasonable theoretical target for aneurysm progression and stability; however, the practical benefits of current antithrombotic therapies in AAA are unclear. Further research is needed to demonstrate the extent to which thrombosis impacts AAA pathogenesis and to develop novel pharmacologic strategies for the medical management of this disease.

Building Blood Vessel Chips with Enhanced Physiological Relevance

Blood vessel chips are bioengineered microdevices, consisting of biomaterials, human cells, and microstructures, which recapitulate essential vascular structure and physiology and allow a well-controlled microenvironment and spatial-temporal readouts. Blood vessel chips afford promising opportunities to understand molecular and cellular mechanisms underlying a range of vascular diseases. The physiological relevance is key to these blood vessel chips that rely on bioinspired strategies and bioengineering approaches to translate vascular physiology into artificial units. Here, we discuss several critical aspects of vascular physiology, including morphology, material composition, mechanical properties, flow dynamics, and mass transport, which provide essential guidelines and a valuable source of bioinspiration for the rational design of blood vessel chips. We also review state-of-art blood vessel chips that exhibit important physiological features of the vessel and reveal crucial insights into the biological processes and disease pathogenesis, including rare diseases, with notable implications for drug screening and clinical trials. We envision that the advances in biomaterials, biofabrication, and stem cells improve the physiological relevance of blood vessel chips, which, along with the close collaborations between clinicians and bioengineers, enable their widespread utility.

Absence of hyperfibrinolysis may explain lack of efficacy of tranexamic acid in hypoproliferative thrombocytopenia

The American Trial Using Tranexamic Acid (TXA) in Thrombocytopenia (A-TREAT, NCT02578901) demonstrated no superiority of TXA over placebo in preventing World Health Organization (WHO) grade 2 or higher bleeding in patients with severe thrombocytopenia requiring supportive platelet transfusion following myeloablative therapy for hematologic disorders. In this ancillary study, we sought to determine whether this clinical outcome could be explained on the basis of correlative assays of fibrinolysis. Plasma was collected from A-TREAT participants (n = 115) before the initiation of study drug (baseline) and when TXA was at steady-state trough concentration (follow-up). Global fibrinolysis was measured by 3 assays: euglobulin clot lysis time (ECLT), plasmin generation (PG), and tissue-type plasminogen activator (tPA)-challenged clot lysis time (tPA-CLT). TXA was quantified in follow-up samples by tandem mass spectrometry. Baseline samples did not demonstrate fibrinolytic activation by ECLT or tPA-CLT. Furthermore, neither ECLT nor levels of plasminogen activator inhibitor-1, tPA, plasminogen, alpha2-antiplasmin, or plasmin-antiplasmin complexes were associated with a greater risk of WHO grade 2+ bleeding. TXA trough concentrations were highly variable (range, 0.7-10 μg/mL) and did not correlate with bleeding severity, despite the fact that plasma TXA levels correlated strongly with pharmacodynamic assessments by PG (Spearman r, -0.78) and tPA-CLT (r, 0.74). We conclude that (1) no evidence of fibrinolytic activation was observed in these patients with thrombocytopenia, (2) trough TXA concentrations varied significantly between patients receiving the same dosing schedule, and (3) tPA-CLT and PG correlated well with TXA drug levels.

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.

Superhydrophobic Biological Fluid-Repellent Surfaces: Mechanisms and Applications

Superhydrophobic biological fluid-repellent surfaces (SBFRSs) have attracted great attention in the treatment of blood and urine-related diseases because of their unique wettability and compatibility, which creates a new path for the development of medical apparatus and instruments, and are expected to create advances in various fields. Here, this review provides an up-to-date summary of research progress on the repellent mechanism and application of SBFRSs. The underlying physical and chemical principles for designing superhydrophobic surfaces are first introduced. Then, the dialectical influences of solid-liquid interactions between superhydrophobic surfaces and biological fluids on the wettability and compatibility are emphatically expounded. Subsequently, attention is drawn to the recent applications of SBFRSs in biomedical fields, such as surgical medical apparatus, implant materials, extracorporeal circulation devices, and biological fluid detection. Finally, the outlook and challenges in terms of employing SBFRSs are also discussed. This review is expected to provide a comprehensive guidance for the preparation of SBFRSs with compatibility and long-term superhydrophobic stability that is closely related to clinical applications.

The advances of blood clots used as biomaterials in regenerative medicine

The physiologic process of blood clot formation is well understood and occurs naturally in the setting of tissue injury to achieve hemostasis and begin the process of wound healing. While the investigation of blood clots as a biomaterial is still in the early stages, there has been some research with similar biomaterials made of the components of blood clots that support the innovative idea of using an autologous blood clot as a scaffold or delivery method for therapeutic agents. Here, we review the physiology of blood clots in wound healing and how using blood clots as a biomaterial and delivery system can potentially promote wound healing, provide targeted therapeutic agent delivery and use it as an innovative tool in regenerative medicine.

A Nanomedicine Structure-Activity Framework for Research, Development, and Regulation of Future Cancer Therapies

Despite their clinical success in drug delivery applications, the potential of theranostic nanomedicines is hampered by mechanistic uncertainty and a lack of science-informed regulatory guidance. Both the therapeutic efficacy and the toxicity of nanoformulations are tightly controlled by the complex interplay of the nanoparticle's physicochemical properties and the individual patient/tumor biology; however, it can be difficult to correlate such information with observed outcomes. Additionally, as nanomedicine research attempts to gradually move away from large-scale animal testing, the need for computer-assisted solutions for evaluation will increase. Such models will depend on a clear understanding of structure-activity relationships. This review provides a comprehensive overview of the field of cancer nanomedicine and provides a knowledge framework and foundational interaction maps that can facilitate future research, assessments, and regulation. By forming three complementary maps profiling nanobio interactions and pathways at different levels of biological complexity, a clear picture of a nanoparticle's journey through the body and the therapeutic and adverse consequences of each potential interaction are presented.

Platelet-Neutrophil Association in NETs-Rich Areas in the Retrieved AIS Patient Thrombi

Histological structure of thrombi is a strong determinant of the outcome of vascular recanalization therapy, the only treatment option for acute ischemic stroke (AIS) patients. A total of 21 AIS patients from this study after undergoing non-enhanced CT scan and multimodal MRI were treated with mechanical stent-based and manual aspiration thrombectomy, and thromboembolic retrieved from a cerebral artery. Complementary histopathological and imaging analyses were performed to understand their composition with a specific focus on fibrin, von Willebrand factor, and neutrophil extracellular traps (NETs). Though distinct RBC-rich and platelet-rich areas were found, AIS patient thrombi were overwhelmingly platelet-rich, with 90% of thrombi containing <40% total RBC-rich contents (1.5 to 37%). Structurally, RBC-rich areas were simple, consisting of tightly packed RBCs in thin fibrin meshwork with sparsely populated nucleated cells and lacked any substantial von Willebrand factor (VWF). Platelet-rich areas were structurally more complex with thick fibrin meshwork associated with VWF. Plenty of leukocytes populated the platelet-rich areas, particularly in the periphery and border areas between platelet-rich and RBC-rich areas. Platelet-rich areas showed abundant activated neutrophils (myeloperoxidase+ and neutrophil-elastase+) containing citrullinated histone-decorated DNA. Citrullinated histone-decorated DNA also accumulated extracellularly, pointing to NETosis by the activated neutrophils. Notably, NETs-containing areas showed strong reactivity to VWF, platelets, and high-mobility group box 1 (HMGB1), signifying a close interplay between these components.

Heparin-Peptide Nanogranules for Thrombosis-Actuated Anticoagulation

Despite nearly a century of clinical use as a blood thinner, heparin's rapid serum clearance and potential to induce severe bleeding events continue to urge the development of more effective controlled delivery strategies. Subcutaneous depots that steadily release the anticoagulant into circulation represent a promising approach to reducing overdose frequency, sustaining therapeutic concentrations of heparin in plasma, and prolonging anticoagulant activity in a safe and effective manner. Subcutaneously deliverable heparin-peptide nanogranules that allow for long-lasting heparin bioavailability in the circulatory system, while enabling on-demand activation of heparin's anticoagulant effects in the thrombus microenvironment, are reported. Biophysical studies demonstrate this responsive behavior is due to the sequestration of heparin within self-assembling peptide nanofibrils and its mechanically actuated decoupling to elicit antithrombotic effects at the clotting site. In vivo studies show these unique properties converge to allow subcutaneous nanogranule depots to extend heparin serum concentrations for an order of magnitude longer than standard dosing regimens while enabling prolonged and controlled anticoagulant activity. This biohybrid delivery system demonstrates a potentially scalable platform for the development of safer, easier to administer, and more effective antithrombotic nanotechnologies.

Model reduction of coagulation cascade based on genetic algorithm

Fibrin is an important product of the coagulation cascade, and plays an eminent role in platelet stabilization. Since coagulation cascade models typically involve the reaction kinetics of dozens of proteins, which will incur burdensome computational costs when coupled to blood flow in complex geometries, researchers often ignore this process when constructing thrombosis models. However, previous studies have shown that fundamental aspects of coagulation can be reproduced with simpler models, which motivated us to obtain a reduced-order model of fibrin generation through a systematic approach. Therefore, we introduced a semi-automatic framework to perform model-reduction of cascade reactions in this study, which consisted of two processes. Specifically, the retained protein species and cascade reactions were determined based on published studies and simulation results from the full cascade model, while the optimal reaction rates for the new cascade network were determined using a genetic algorithm. The framework has been applied to a 19-species coagulation model that triggers fibrin generation in internal fields via reactive boundaries, and a 10-species reduced-order model was obtained to reproduce the kinetics of fibrinogenesis in the full cascade model at different boundary tissue factor concentrations. This reduced-order model of fibrinogenesis would be valuable for thrombosis modeling that considers both the coagulation cascade and platelet activity. Furthermore, the framework proposed herein can also be applied to the reductions of other cascade reaction models.

Interaction of von Willebrand factor with blood cells in flow models: a systematic review

The presence of blood flow influences the interaction between von Willebrand factor (VWF) and blood cells, affecting characteristics of forming blood clots. The interactions between coagulation and inflammation have mainly been studied in thrombosis models, but it remains unclear whether these interactions might also play a role in reduced bleeding in patients with bleeding disorders. In this systematic review, we provide an overview of the literature investigating the interactions between VWF and blood cells in flow models. For article selection, a systematic search was performed in Embase, Medline-Ovid, Cochrane Library, Web of Science databases, and Google Scholar. After selection, 24 articles were included. These articles describe direct or platelet-dependent interactions between VWF and neutrophils, monocytes, erythrocytes, or lymphocytes under different flow conditions. Almost all the described interactions required the presence of activated platelets. Only erythrocytes, monocytes, and natural killer cells were capable of directly binding the VWF multimers. Overall, interactions between VWF and blood cells mainly occurred in the presence of platelets. Because of the large variation in study design and used flow rates, further research is necessary to compare the results between studies and draw firm conclusions on when and under what conditions these interactions can occur. After our findings, many questions remained unanswered. This review might provide a starting point for future research. Extended knowledge on the influence of blood flow on VWF and blood cell interactions can contribute to improved understanding of the variation in bleeding in patients with bleeding disorders.

Nuclear Progestin Receptor-mediated Linkage of Blood Coagulation and Ovulation

Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, ~3145-fold) and tissue factor (f3a, ~120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wild-type preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting-triggered ovulation in vertebrates.

Computational Analysis of the Related Factors of Deep Vein Thrombosis (DVT) Formation in Patients Undergoing Hip Fracture Surgery

A retrospective study was conducted on 51 patients undergoing hip fracture surgery to investigate the factors associated with the formation of deep venous thrombosis (DVT). The independent sample t-test and correlation analysis were used to sort out and analyze the data. The findings are as follows. (1) Different gender samples showed significant differences in the Caprini score and thrombus location. Most DVTs in females are located in the posterior tibial vein and intermuscular veins. The Caprini score of females was significantly higher than that of males. (2) Age displays a positive correlation with DVT, coronary heart disease, hypertension, and different surgical types, respectively. (3) There is a correlation between age and operation duration. (4) Hyperlipidemia and cerebrovascular disease show a positive correlation with DVT. (5) There was a significant negative correlation between the Caprini score and the quantification of D-dimer. This indicates that in this sample, the higher the patients' Caprini score is, the lower the quantitation of D-dimer will be. (6) Hyperlipidemia and cardiac insufficiency show a positive correlation with cerebrovascular disease. Patients with hyperlipidemia and cardiac insufficiency may also suffer from cerebrovascular diseases.

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.

Models of Shear-Induced Platelet Activation and Numerical Implementation With Computational Fluid Dynamics Approaches

Shear-induced platelet activation is one of the critical outcomes when blood is exposed to elevated shear stress. Excessively activated platelets in the circulation can lead to thrombus formation and platelet consumption, resulting in serious adverse events such as thromboembolism and bleeding. While experimental observations reveal that it is related to the shear stress level and exposure time, the underlying mechanism of shear-induced platelet activation is not fully understood. Various models have been proposed to relate shear stress levels to platelet activation, yet most are modified from the empirically calibrated power-law model. Newly developed multiscale platelet models are tested as a promising approach to capture a single platelet's dynamic shape during activation, but it would be computationally expensive to employ it for a large-scale analysis. This paper summarizes the current numerical models used to study the shear-induced platelet activation and their computational applications in the risk assessment of a particular flow pattern and clot formation prediction.

Ethnobotany, phytochemistry, pharmacology, and conservation of the genus Calanthe R. Br. (Orchidaceae)

ETHNOPHARMACOLOGICAL RELEVANCE

The Genus Calanthe (family Orchidaceae) consists of more than 207 species distributed in both tropical and subtropical regions. In traditional medicine, Calanthe species provide remedies against various conditions such as arthritis, rheumatism, traumatic injuries, snake-bites, abdominal discomfort, nose bleeding, common colds, ulcers, chronic coughs, and others. Some species are also used as aphrodisiacs, tonics, and as pain relievers on joints and toothaches.

AIM OF THE REVIEW

This review provides comprehensive information on the herbal uses, chemical components, pharmacological activities, and conservation of Calanthe, which might be useful in the future development of potent herbal medicines and facilitate the enactment of better conservation strategies.

MATERIALS AND METHODS

Relevant information was obtained from online databases including SCI-Finder, Google Scholar, Web of Science, Science Direct, PubMed, Springer, IOP Science, and other web sources such as PubChem, The Plant List, and World Flora Online. Books, Ph.D. and MSc dissertations were used for unpublished literature. Information from Chinese literature was obtained from the CNKI database.

RESULTS

In total, 19 species of the genus Calanthe have been reported to be used in traditional medicine in different countries of Asia. A total of 265 chemical compounds from different chemical classes including, alkaloids, terpenoids, phenolic compounds and phenolic derivatives, phenanthrenes, and others, have been identified from Calanthe species. Calanquinone A isolated from C. arisanensis has been reported to exhibit antitumor activity against six malignant cell lines. Other bioactive compounds from Calanthe with pharmacological activity include phenanthrenes, phenanthrenequinones, 6'-O-β-D-apiofuranosylindican, 4H-Pyran-4one, 2, 3-dihydro-3,5 dihydroxy-6-methyl, and calanthoside. These compounds exhibit valuable biological properties such as hair restoration, anticancer activity, anti-inflammatory and antiarthritic activity, antidiabetic and hepatoprotective potency, antiplatelet aggregation action, and antibacterial and antifungal activities. Some Calanthe species, including C. ecallosa and C. yuana, are endangered in the IUCN red list. The high risk of extinction is attributed to illegal trade and unsustainable harvesting and utilization.

CONCLUSIONS

This review summarizes the herbal uses, chemical components, biological activity, and conservation of Calanthe. The pharmacological studies on this genus are limited; thus, extensive research on the toxicology, pharmaceutical standardization, and mechanism of action of the isolated bioactive compounds are needed. Since some species of Calanthe are listed as endangered, stringent guidelines on trade, collection, and sustainable utilization of medicinal orchids should be set up to facilitate the conservation of these species.

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.

Improved Drainage Cannula Design to Reduce Thrombosis in Veno-Arterial Extracorporeal Membrane Oxygenation

Thrombosis is a potentially life-threatening complication in veno-arterial extracorporeal membrane oxygenation (ECMO) circuits, which may originate from the drainage cannula due to unfavorable blood flow dynamics. This study aims to numerically investigate the effect of cannula design parameters on local fluid dynamics, and thus thrombosis potential, within ECMO drainage cannulas. A control cannula based on the geometry of a 17 Fr Medtronic drainage cannula concentrically placed in an idealized, rigid-walled geometry of the right atrium and superior and inferior vena cava was numerically modeled. Simulated flow dynamics in the control cannula were systematically compared with 10 unique cannula designs which incorporated changes to side hole diameter, the spacing between side holes, and side hole angles. Local blood velocities, maximum wall shear stress (WSS), and blood residence time were used to predict the risk of thrombosis. Numerical results were experimentally validated using particle image velocimetry. The control cannula exhibited low blood velocities (59 mm/s) at the cannula tip, which may promote thrombosis. Through a reduction in the side hole diameter (2 mm), the spacing between the side holes (3 mm) and alteration in the side hole angle (30° relative to the flow direction), WSS was reduced by 52%, and cannula tip blood velocity was increased by 560% compared to the control cannula. This study suggests that simple geometrical changes can significantly alter the risk of thrombosis in ECMO drainage cannulas.

DEVELOPMENT OF FIBRIN BRANCH STRUCTURE BEFORE AND AFTER GELATION

In [Fogelson and Keener, Phys. Rev. E, 81 (2010), 051922], we introduced a kinetic model of fibrin polymerization during blood clotting that captured salient experimental observations about how the gel branching structure depends on the conditions under which the polymerization occurs. Our analysis there used a moment-based approach that is valid only before the finite time blow-up that indicates formation of a gel. Here, we extend our analyses of the model to include both pre-gel and post-gel dynamics using the PDE-based framework we introduced in [Fogelson and Keener, SIAM J. Appl. Math., 75 (2015), pp. 1346-1368]. We also extend the model to include spatial heterogeneity and spatial transport processes. Studies of the behavior of the model reveal different spatial-temporal dynamics as the time scales of the key processes of branch formation, monomer introduction, and diffusion are varied.

Cerebral Aneurysm Occlusion at 12-Month Follow-Up After Flow-Diverter Treatment: Statistical Modeling for V&V With Real-World Data

Background: Flow-Diverter (FD) porosity has been pointed as a critical factor in the occlusion of cerebral aneurysms after treatment.

Objective: Verification and Validation of computational models in terms of predictive capacity, relating FD porosity and occlusion after cerebral aneurysms treatment.

Methods: Sixty-four aneurysms, with pre-treatment and follow-up images, were considered. Patient demographics and aneurysm morphological information were collected. The computational simulation provided by ANKYRAS provided FD porosity, expansion, and mesh angle. FD occlusion was assessed and recorded from follow-up images. Multiple regression Logit and analysis of covariance (ANCOVA) models were used to model the data with both categorical and continuous models.

Results: Occlusion of the aneurysm after 12 months was affected by aneurysm morphology but not by FD mesh morphology. A Time-To-Occlusion (TTO) of 6.92 months on average was observed with an SE of 0.24 months in the aneurysm population surveyed. TTO was estimated with statistical significance from the resulting model for the data examined and was capable of explaining 92% of the data variation.

Conclusions: Porosity was found to have the most correction power when assessing TTO, proving its importance in the process of aneurysm occlusion. Still, further Verification and Validation (V&V) of treatment simulation in more extensive, multi-center, and randomized databases is required.

A Hybrid Semi-Lagrangian Cut Cell Method for Advection-Diffusion Problems with Robin Boundary Conditions in Moving Domains

We present a new discretization approach to advection-diffusion problems with Robin boundary conditions on complex, time-dependent domains. The method is based on second order cut cell finite volume methods introduced by Bochkov et al. [8] to discretize the Laplace operator and Robin boundary condition. To overcome the small cell problem, we use a splitting scheme along with a semi-Lagrangian method to treat advection. We demonstrate second order accuracy in the L ¹, L ², and L ^∞ norms for both analytic test problems and numerical convergence studies. We also demonstrate the ability of the scheme to convert one chemical species to another across a moving boundary.

Continuum modeling of thrombus formation and growth under different shear rates

Obstruction of blood flow due to thrombosis is a major cause of ischemic stroke, myocardial infarction, and in severe cases, mortality. In particular, in blood wetting medical devices, thrombosis is a common reason for failure. The prediction of thrombosis by understanding signaling pathways using computational models, lead to identify the risk of thrombus formation in blood-contacting devices and design improvements. In this study, a mathematical model of thrombus formation and growth is presented. A biochemical model of platelet activation and aggregation is developed to predict thrombus size and shape at the site of vascular injury. Computational fluid dynamics using the finite volume method is employed to compute the velocity and pressure fields which are influenced by the growing thrombi. The passive transport of platelets, agonists, the platelet activation kinetics, their adhesion to the growing thrombi and embolization of platelets are solved by a fully coupled set of convection-diffusion-reaction equations. The thrombogenic surface representing blood-contacting material or injured blood vessel was incorporated into the model as a surface flux boundary condition to initiate thrombus formation. The blood is considered as a Newtonian fluid, while the thrombus is treated as a porous medium. The results are compared with in vitro experiments of a microfluidic chamber at an initial inlet venous shear rate of 200s^-1 using a pressure-inlet boundary condition. The thrombus development due to agonist concentrations and change in the shear rate as well as thromboembolism for this benchmark problem is successfully computed. Furthermore, to extend the current model to a physiologically relevant configuration, thrombus formation in a blood tube is simulated. Two different heterogeneous reaction rates for platelet aggregation are used to simulate thrombus growth under a constant inlet flow rate. The findings show that the thrombus shape can be substantially altered by the hemodynamic conditions, increase in the shear rate and due to the combined effects of shear induced platelet activation and the heterogeneous reaction rates. It is also concluded that the model is able to predict thrombus formation in different physiological and pathological hemodynamics.

Computational investigation of blood cell transport in retinal microaneurysms

Microaneurysms (MAs) are one of the earliest clinically visible signs of diabetic retinopathy (DR). MA leakage or rupture may precipitate local pathology in the surrounding neural retina that impacts visual function. Thrombosis in MAs may affect their turnover time, an indicator associated with visual and anatomic outcomes in the diabetic eyes. In this work, we perform computational modeling of blood flow in microchannels containing various MAs to investigate the pathologies of MAs in DR. The particle-based model employed in this study can explicitly represent red blood cells (RBCs) and platelets as well as their interaction in the blood flow, a process that is very difficult to observe in vivo. Our simulations illustrate that while the main blood flow from the parent vessels can perfuse the entire lumen of MAs with small body-to-neck ratio (BNR), it can only perfuse part of the lumen in MAs with large BNR, particularly at a low hematocrit level, leading to possible hypoxic conditions inside MAs. We also quantify the impacts of the size of MAs, blood flow velocity, hematocrit and RBC stiffness and adhesion on the likelihood of platelets entering MAs as well as their residence time inside, two factors that are thought to be associated with thrombus formation in MAs. Our results show that enlarged MA size, increased blood velocity and hematocrit in the parent vessel of MAs as well as the RBC-RBC adhesion promote the migration of platelets into MAs and also prolong their residence time, thereby increasing the propensity of thrombosis within MAs. Overall, our work suggests that computational simulations using particle-based models can help to understand the microvascular pathology pertaining to MAs in DR and provide insights to stimulate and steer new experimental and computational studies in this area.

Recent advances in blood rheology: a review

Due to the potential impact on the diagnosis and treatment of various cardiovascular diseases, work on the rheology of blood has significantly expanded in the last decade, both experimentally and theoretically. Experimentally, blood has been confirmed to demonstrate a variety of non-Newtonian rheological characteristics, including pseudoplasticity, viscoelasticity, and thixotropy. New rheological experiments and the development of more controlled experimental protocols on more extensive, broadly physiologically characterized, human blood samples demonstrate the sensitivity of aspects of hemorheology to several physiological factors. For example, at high shear rates the red blood cells elastically deform, imparting viscoelasticity, while at low shear rates, they form "rouleaux" structures that impart additional, thixotropic behavior. In addition to the advances in experimental methods and validated data sets, significant advances have also been made in both microscopic simulations and macroscopic, continuum, modeling, as well as novel, multiscale approaches. We outline and evaluate the most promising of these recent developments. Although we primarily focus on human blood rheology, we also discuss recent observations on variations observed across some animal species that provide some indication on evolutionary effects.

Mathematical and Computational Modeling of Device-Induced Thrombosis

Given the extensive and routine use of cardiovascular devices, a major limiting factor to their success is the thrombotic rate that occurs. This both poses direct risk to the patient and requires counterbalancing with anticoagulation and other treatment strategies, contributing additional risks. Developing a better understanding of the mechanisms of device-induced thrombosis to aid in device design and medical management of patients is critical to advance the ubiquitous use and durability. Thus, mathematical and computational modelling of device-induced thrombosis has received significant attention recently, but challenges remain. Additional areas that need to be explored include microscopic/macroscopic approaches, reconciling physical and numerical timescales, immune/inflammatory responses, experimental validation, and incorporating pathologies and blood conditions. Addressing these areas will provide engineers and clinicians the tools to provide safe and effective cardiovascular devices.

Healing through the lens of immunothrombosis: Biology-inspired, evolution-tailored, and human-engineered biomimetic therapies

Evolution, from invertebrates to mammals, has yielded and shaped immunoclotting as a defense and repair response against trauma and infection. This mosaic of immediate and local wound-sealing and pathogen-killing mechanisms results in survival, restoration of homeostasis, and tissue repair. In mammals, immunoclotting has been complemented with the neuroendocrine system, platelets, and contact system among other embellishments, adding layers of complexity through interconnecting blood-born proteolytic cascades, blood cells, and the neuroendocrine system. In doing so, immunothrombosis endows humans with survival advantages, but entails vulnerabilities in the current unprecedented and increasingly challenging environment. Immunothrombosis and tissue repair appear to go hand in hand with common mechanisms mediating both processes, a fact that is underlined by recent advances that are deciphering the mechanisms of the repair process and of the biochemical pathways that underpins coagulation, hemostasis and thrombosis. This review is intended to frame both the universal aspects of tissue repair and the therapeutic use of autologous fibrin matrix as a biology-as-a-drug approach in the context of the evolutionary changes in coagulation and hemostasis. In addition, we will try to shed some light on the molecular mechanisms underlying the use of the autologous fibrin matrix as a biology-inspired, evolution-tailored, and human-engineered biomimetic therapy.

Semi-empirical anticoagulation model (SAM): INR monitoring during Warfarin therapy

The International Normalized Ratio (INR) monitoring is an essential component to manage thrombotic disease therapy. This study presents a semi-empirical model of INR as a function of time and assigned therapy (Warfarin, k-vitamin). With respect to other methodologies, this model is able to describe the INR using a limited number of parameters and is able to describe the time variation of INR described in the literature. The presented methodology showed great accuracy in model calibration [(trueness (precision)]: 0.2% (0.1%) to 1.2% (0.3%) for coagulation factors, from 5% (9%) to 9.7% (12%) for Warfarin-related parameters and 38% (40%) for K-vitamin-related parameters. The latter value was considered acceptable given the assumptions made in the model. It has two other important results: the first is that it was able to correctly estimate INR with respect to daily therapy doses taken from the literature. The second is that it introduces a single numeric semi-empirical parameter that is able to correlate INR/dose response to physiological and environmental condition of patients.

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.

Heterogeneity between proximal and distal aspects of occlusive thrombi on pretreatment imaging in acute ischemic stroke

BACKGROUND

The potential heterogeneity in occlusive thrombi caused by in situ propagation by secondary thrombosis after embolic occlusion could obscure the characteristics of original thrombi, preventing the clarification of a specific thrombus signature for the etiology of ischemic stroke. We aimed to investigate the heterogeneity of occlusive thrombi by pretreatment imaging.

METHODS

Among consecutive stroke patients with acute embolic anterior circulation large vessel occlusion treated with thrombectomy, we retrospectively reviewed 104 patients with visible occlusive thrombi on pretreatment non-contrast computed tomography admitted from January 2015 to December 2018. A region of interest was set on the whole thrombus on non-contrast computed tomography under the guidance of computed tomography angiography. The region of interest was divided equally into the proximal and distal segments and the difference in Hounsfield unit densities between the two segments was calculated.

RESULTS

Hounsfield unit density in the proximal segment was higher than that in the distal segment (mean difference 4.45; p < 0.001), regardless of stroke subtypes. On multivariate analysis, thrombus length was positively correlated (β = 0.25; p < 0.001) and time from last-known-well to imaging was inversely correlated (β = -0.0041; p = 0.002) with the difference in Hounsfield unit densities between the proximal and distal segments.

CONCLUSIONS

The difference in density between the proximal and distal segments increased as thrombi became longer and decreased as thrombi became older after embolic occlusion. This time/length-dependent thrombus heterogeneity between the two segments is suggestive of secondary thrombosis initially occurring on the proximal side of the occlusion.

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.