A constitutive model for developing blood clots with various compositions and their nonlinear viscoelastic behavior

Mechanical properties of clot made from human and bovine whole blood differ significantly

Thromboembolism - that is, clot formation and the subsequent fragmentation of clot - is a leading cause of death worldwide. Clots' mechanical properties are critical determinants of both the embolization process and the pathophysiological consequences thereof. Thus, understanding and quantifying the mechanical properties of clots is important to our ability to treat and prevent thromboembolic disease. However, assessing these properties from in vivo clots is experimentally challenging. Therefore, we and others have turned to studying in vitro clot mimics instead. Unfortunately, there are significant discrepancies in the reported properties of these clot mimics, which have been hypothesized to arise from differences in experimental techniques and blood sources. The goal of our current work is therefore to compare the mechanical behavior of clots made from the two most common sources, human and bovine blood, using the same experimental techniques. To this end, we tested clots under pure shear with and without initial cracks, under cyclic loading, and under stress relaxation. Based on these data, we computed and compared stiffness, strength, work-to-rupture, fracture toughness, relaxation time constants, and prestrain. While clots from both sources behaved qualitatively similarly, they differed quantitatively in almost every metric. We also correlated each mechanical metric to measures of blood composition. Thereby, we traced this inter-species variability in clot mechanics back to significant differences in hematocrit, but not platelet count. Thus, our work suggests that the results of past studies that have used bovine blood to make in vitro mimics - without adjusting blood composition - should be interpreted carefully. Future studies about the mechanical properties of blood clots should focus on human blood alone.

Clots reveal anomalous elastic behavior of fiber networks

The adaptive mechanical properties of soft and fibrous biological materials are relevant to their functionality. The emergence of the macroscopic response of these materials to external stress and intrinsic cell traction from local deformations of their structural components is not well understood. Here, we investigate the nonlinear elastic behavior of blood clots by combining microscopy, rheology, and an elastic network model that incorporates the stretching, bending, and buckling of constituent fibrin fibers. By inhibiting fibrin cross-linking in blood clots, we observe an anomalous softening regime in the macroscopic shear response as well as a reduction in platelet-induced clot contractility. Our model explains these observations from two independent macroscopic measurements in a unified manner, through a single mechanical parameter, the bending stiffness of individual fibers. Supported by experimental evidence, our mechanics-based model provides a framework for predicting and comprehending the nonlinear elastic behavior of blood clots and other active biopolymer networks in general.

Quick model-based viscoelastic clot strength predictions from blood protein concentrations for cybermedical coagulation control

Cybermedical systems that regulate patient clotting in real time with personalized blood product delivery will improve treatment outcomes. These systems will harness popular viscoelastic assays of clot strength such as thromboelastography (TEG), which help evaluate coagulation status in numerous conditions: major surgery (e.g., heart, vascular, hip fracture, and trauma); liver cirrhosis and transplants; COVID-19; ICU stays; sepsis; obstetrics; diabetes; and coagulopathies like hemophilia. But these measurements are time-consuming, and thus impractical for urgent care and automated coagulation control. Because protein concentrations in a blood sample can be measured in about five minutes, we develop personalized, phenomenological, quick, control-oriented models that predict TEG curve outputs from input blood protein concentrations, to facilitate treatment decisions based on TEG curves. Here, we accurately predict, experimentally validate, and mechanistically justify curves and parameters for common TEG assays (Functional Fibrinogen, Citrated Native, Platelet Mapping, and Rapid TEG), and verify results with trauma patient clotting data.

Probing interactions of red blood cells and contracting fibrin platelet clots

Contraction of blood clots plays an important role in blood clotting, a natural process that restores hemostasis and regulates thrombosis in the body. Upon injury, a chain of events culminate in the formation of a soft plug of cells and fibrin fibers attaching to wound edges. Platelets become activated and apply contractile forces to shrink the overall clot size, modify clot structure, and mechanically stabilize the clot. Impaired blood clot contraction results in unhealthy volumetric, mechanical, and structural properties of blood clots associated with a range of severe medical conditions for patients with bleeding and thrombotic disorders. Due to the inherent mechanical complexity of blood clots and a confluence of multiple interdependent factors governing clot contraction, the mechanics and dynamics of clot contraction and the interactions with red blood cells (RBCs) remain elusive. Using an experimentally informed, physics-based mesoscale computational model, we probe the dynamic interactions among platelets, fibrin polymers, and RBCs, and examine the properties of contracted blood clots. Our simulations confirm that RBCs strongly affect clot contraction. We find that RBC retention and compaction in thrombi can be solely a result of mechanistic contraction of fibrin mesh due to platelet activity. Retention of RBCs hinders clot contraction and reduces clot contractility. Expulsion of RBCs located closer to clot outer surface results in the development of a dense fibrin shell in thrombus clots commonly observed in experiments. Our simulations identify the essential parameters and interactions that control blood clot contraction process, highlighting its dependence on platelet concentration and the initial clot size. Furthermore, our computational model can serve as a useful tool in clinically relevant studies of hemostasis and thrombosis disorders, and post thrombotic clot lysis, deformation, and breaking.

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.

Elasticity of whole blood clots measured via Volume Controlled Cavity Expansion

Measuring and understanding the mechanical properties of blood clots can provide insights into disease progression and the effectiveness of potential treatments. However, several limitations hinder the use of standard mechanical testing methods to measure the response of soft biological tissues, like blood clots. These tissues can be difficult to mount, and are inhomogeneous, irregular in shape, scarce, and valuable. To remedy this, we employ in this work Volume Controlled Cavity Expansion (VCCE), a technique that was recently developed, to measure local mechanical properties of soft materials in their natural environment. Through highly controlled volume expansion of a water bubble at the tip of an injection needle, paired with simultaneous measurement of the resisting pressure, we obtain a local signature of whole blood clot mechanical response. Comparing this data with predictive theoretical models, we find that a 1-term Ogden model is sufficient to capture the nonlinear elastic response observed in our experiments and produces shear modulus values that are comparable to values reported in the literature. Moreover, we find that bovine whole blood stored at 4 °C for greater than 2 days exhibits a statistically significant shift in the shear modulus from 2.53 ± 0.44 kPa on day 2 (N = 13) to 1.23 ± 0.18 kPa on day 3 (N = 14). In contrast to previously reported results, our samples did not exhibit viscoelastic rate sensitivity within strain rates ranging from 0.22 - 21.1 s-1. By surveying existing data on whole blood clots for comparison, we show that this technique provides highly repeatable and reliable results, hence we propose the more widespread adoption of VCCE as a path forward to building a better understanding of the mechanics of soft biological materials.

Computational analysis of effects of clot length on Acute ischemic stroke recanalization under different cyclic aspiration loading conditions

Acute ischemic stroke, the second leading cause of death worldwide, results from occlusion of a cerebral artery by a blood clot. Application of cyclic aspiration using an aspiration catheter is a current therapy for the removal of lodged clots. In this study, we perform finite element simulations to analyze deformation of long clots, having length to radius ratio of 2-10, which corresponds to clot-length of 2.85-14.25 mm, under peak-to-peak cyclic aspiration pressures of 10-50 mmHg, and frequencies of 0.5, 1, and 2 Hz. Our computational system comprises of a nonlinear viscoelastic solid clot, a hyperelastic artery, and a nonlinear viscoelastic cohesive zone, the latter modeling the clot-artery interface. We observe that clots having length-to-radius ratio approximately greater than two separate from the inner arterial surface somewhere between the axial and distal ends, irrespective of the cyclic aspiration loading conditions. The stress distribution within the clot shows large tensile stresses in the clot interior, indicating the possibility of simultaneous fragmentation of the clot. Thus, this study shows us the various failure mechanisms simultaneously present in the clot during cyclic aspiration. Similarly, the stress distribution within the artery implies a possibility of endothelial damage to the arterial wall near the end where the aspiration pressure is applied. This framework provides a foundation for further investigation to clot fracture and adhesion characterization.

Passive high-frequency microrheology of blood

Indicative of various pathologies, blood properties are under intense scrutiny. The hemorheological characteristics are traditionally gauged by bulk, low-frequency indicators that average out critical information about the complex, multi-scale, and multi-component structure. In particular, one cannot discriminate between the erythrocytes contribution to global rheology and the impact of plasma. Nevertheless, in their fast stochastic movement, before they encounter each other, the erythrocytes probe the subtle viscoelasticity of their protein-rich environment. Thus, if these short time scales can be resolved experimentally, the plasma properties could be determined without having to separate the blood components; the blood is practically testing itself. This microrheological description of blood plasma provides a direct link between the composition of whole blood and its coagulability status. We present a parametric model for the viscoelasticity of plasma, which is probed by the erythrocytes over frequency ranges of kilohertz in a picoliter-sized volume. The model is validated both in vitro, using artificial hemo-systems where the composition is controlled, as well as on whole blood where continuous measurements provide real-time information. We also discuss the possibility of using this passive microrheology as an in vivo assay for clinically relevant situations where the blood clotting condition must be observed and managed continuously for diagnosis or during therapeutic procedures at different stages of hemostatic and thrombotic processes.

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.

Preclinical modeling of mechanical thrombectomy

Mechanical thrombectomy to treat large vessel occlusions (LVO) causing a stroke is one of the most effective treatments in medicine, with a number needed to treat to improve clinical outcomes as low as 2.6. As the name implies, it is a mechanical solution to a blocked artery and modeling these mechanics preclinically for device design, regulatory clearance and high-fidelity physician training made clinical applications possible. In vitro simulation of LVO is extensively used to characterize device performance in representative vascular anatomies with physiologically accurate hemodynamics. Embolus analogues, validated against clots extracted from patients, provide a realistic simulated use experience. In vitro experimentation produces quantitative results such as particle analysis of distal emboli generated during the procedure, as well as pressure and flow throughout the experiment. Animal modeling, used mostly for regulatory review, allows estimation of device safety. Other than one recent development, nearly all animal modeling does not incorporate the desired target organ, the brain, but rather is performed in the extracranial circulation. Computational modeling of the procedure remains at the earliest stages but represents an enormous opportunity to rapidly characterize and iterate new thrombectomy concepts as well as optimize procedure workflow. No preclinical model is a perfect surrogate; however, models available can answer important questions during device development and have to date been successful in delivering efficacious and safe devices producing excellent clinical outcomes. This review reflects on the developments of preclinical modeling of mechanical thrombectomy with particular focus on clinical translation, as well as articulate existing gaps requiring additional research.

A review on the association of thrombus composition with mechanical and radiological imaging characteristics in acute ischemic stroke

Thrombus composition and mechanical properties significantly impact the ease and outcomes of thrombectomy procedures in patients with acute ischemic stroke. A wide variation exists in the composition of thrombi between patients. If a relationship can be determined between the composition of a thrombus and its mechanical behaviour, as well as between the composition of a thrombus and its radiological imaging characteristics, then there is the potential to personalise thrombectomy treatment based on each individual thrombus. This review aims to give an overview of the current literature addressing this issue. Here, we present a scoping review detailing associations between thrombus composition, mechanical behaviour and radiological imaging characteristics. We conducted two searches 1) on the association between thrombus composition and the mechanical behaviour of the tissue and 2) on the association between radiological imaging characteristics and thrombus composition in the acute stroke setting. The review suggests that higher fibrin and lower red blood cell (RBC) content contribute to stiffer thrombi independent of the loading mode. Further, platelet-contracted thrombi are stiffer than non-contracted compositional counterparts. Fibrin content contributes to the elastic portion of viscoelastic behaviour while RBC content contributes to the viscous portion. It is possible to identify fibrin-rich or RBC-rich thrombi with computed tomography and magnetic resonance imaging vessel signs. Standardisation is required to quantify the association between thrombus density on non-contrast computed tomography and the RBC content. The characterisation of the thrombus fibrin network has not been addressed so far in radiological imaging but may be essential for the prediction of device-tissue interactions and distal thrombus embolization. The association between platelet-driven clot contraction and radiological imaging characteristics has not been explicitly investigated. However, evidence suggests that perviousness may be a marker of clot contraction.

Modeling acute ischemic stroke recanalization through cyclic aspiration

We model the deformation of a thromboembolus lodged in a cerebral artery under the application of aspiration pressure as it would be provided by an aspiration catheter during a mechanical thrombectomy procedure. The system considered consists of (i) a clot modeled as a viscoelastic solid; (ii) an artery modeled as a hyperelastic solid; and (iii) a viscoelastic cohesive interface between the clot and the artery. For the chosen system and geometry, we show that the application of aspiration pressure results in the impingement of the thrombus against the inner arterial wall near the aspiration location. Conditions leading to interfacial failure are nucleated at the distal end of the clot and, depending on the details of the loading conditions, propagate toward the proximal end. The results provide useful information in identifying the circumstances that play a decisive role for clot removal by aspiration alone.

Finite element analysis of blood clots based on the nonlinear visco-hyperelastic model

Mechanical thrombectomy has become the standard treatment for patients with an acute ischemic stroke. In this approach, to remove blood clots, mechanical force is applied using thrombectomy devices, in which the interaction between the clot and the device could significantly affect the clot retrieval performance. It is expected that the finite element method (FEM) could visualize the mechanical interaction by the visualization of the stress transmission from the device to the clot. This research was aimed at verifying the constitutive theory by implementing FEM based on the visco-hyperelastic theory, using a three-dimensional clot model. We used the visco-hyperelastic FEM to reproduce the mechanical behavior of blood clots, as observed in experiments. This study is focused on the mechanical responses of clots under tensile loading and unloading because in mechanical thrombectomy, elongation is assumed to occur locally on the clots during the retrieval process. Several types of cylindrical clots were created by changing the fibrinogen dose. Tensile testing revealed that the stiffness (E_0.45-value) of clots with fibrinogen could be more than three times higher than that of clots without fibrinogen. It was also found that the stiffness was not proportional to the fibrinogen dose. By fitting to the theoretical curve, it was revealed that the Mooney-Rivlin model could reproduce the hyperelastic characteristics of clots well. From the stress-relaxation data, the three-chain Maxwell model could accurately fit the experimental viscoelastic data. FEM, taking the theoretical models into account, was then carried out, and the results matched well with the experimental visco-hyperelastic characteristics of clots under tensile load, reproducing the mechanical hysteresis during unloading, the stress dependence on the strain rate, and the time-dependent stress decrease in the stress-relaxation test.

Synthetic hydrogels as blood clot mimicking wound healing materials

Excessive bleeding-or hemorrhage-causes millions of civilian and non-civilian casualties every year. Additionally, wound sequelae, such as infections, are a significant source of chronic morbidity, even if the initial bleeding is successfully stopped. To treat acute and chronic wounds, numerous wound healing materials have been identified, tested, and adopted. Among them are topical dressings, such as gauzes, as well as natural and biomimetic materials. However, none of these materials successfully mimic the complex and dynamic properties of the body's own wound healing material: the blood clot. Specifically, blood clots exhibit complex mechanical and biochemical properties that vary across spatial and temporal scales to guide the wound healing response, which make them the ideal wound healing material. In this manuscript, we review blood clots' complex mechanical and biochemical properties, review current wound healing materials, and identify opportunities where new materials can provide additional functionality, with a specific focus on hydrogels. We highlight recent developments in synthetic hydrogels that make them capable of mimicking a larger subset of blood clot features: as plugs and as stimuli for tissue repair. We conclude that future hydrogel materials designed to mimic blood clot biochemistry, mechanics, and architecture can be combined with exciting platelet-like particles to serve as hemostats that also promote the biological wound healing response. Thus, we believe synthetic hydrogels are ideal candidates to address the clear need for better wound healing materials.

Hyper-viscoelastic damage modeling of whole blood clot under large deformation

Blood clots play a diametric role in our bodies as they are both vital as a wound sealant, as well as the source for many devastating diseases. In blood clots' physiological and pathological roles, their mechanics play a critical part. These mechanics are non-trivial owing to blood clots' complex nonlinear, viscoelastic behavior. Casting this behavior into mathematical form is a fundamental step toward a better basic scientific understanding of blood clots, as well as toward diagnostic and prognostic computational models. Here, we identify a hyper-viscoelastic damage model that we fit to original data on the nonlinear, viscoelastic behavior of blood clots. Our model combines the classic Ogden hyperelastic constitutive law, a finite viscoelastic model for large deformations, and a non-local, gradient-enhanced damage formulation. By fitting our model to cyclic tensile test data and extension-to-failure data, we inform the model's nine unknown material parameters. We demonstrate the predictability of our model by validating it against unseen cyclic tensile test and stress-relaxation data. Our original data, model formulation, and the identified constitutive parameters of this model are openly available for others to use, which will aid in developing accurate, quantitative simulations of blood clot mechanics.

Analysis of human emboli and thrombectomy forces in large-vessel occlusion stroke

OBJECTIVE

This study's purpose was to improve understanding of the forces driving the complex mechanical interaction between embolic material and current stroke thrombectomy devices by analyzing the histological composition and strength of emboli retrieved from patients and by evaluating the mechanical forces necessary for retrieval of such emboli in a middle cerebral artery (MCA) bifurcation model.

METHODS

Embolus analogs (EAs) were generated and embolized under physiological pressure and flow conditions in a glass tube model of the MCA. The forces involved in EA removal using conventional endovascular techniques were described, analyzed, and categorized. Then, 16 embolic specimens were retrieved from 11 stroke patients with large-vessel occlusions, and the tensile strength and response to stress were measured with a quasi-static uniaxial tensile test using a custom-made platform. Embolus compositions were analyzed and quantified by histology.

RESULTS

Uniaxial tension on the EAs led to deformation, elongation, thinning, fracture, and embolization. Uniaxial tensile testing of patients' emboli revealed similar soft-material behavior, including elongation under tension and differential fracture patterns. At the final fracture of the embolus (or dissociation), the amount of elongation, quantified as strain, ranged from 1.05 to 4.89 (2.41 ± 1.04 [mean ± SD]) and the embolus-generated force, quantified as stress, ranged from 63 to 2396 kPa (569 ± 695 kPa). The ultimate tensile strain of the emboli increased with a higher platelet percentage, and the ultimate tensile stress increased with a higher fibrin percentage and decreased with a higher red blood cell percentage.

CONCLUSIONS

Current thrombectomy devices remove emboli mostly by applying linear tensile forces, under which emboli elongate until dissociation. Embolus resistance to dissociation is determined by embolus strength, which significantly correlates with composition and varies within and among patients and within the same thrombus. The dynamic intravascular weakening of emboli during removal may lead to iatrogenic embolization.

Blood Clot Phenotyping by Rheometry: Platelets and Fibrinogen Chemistry Affect Stress-Softening and -Stiffening at Large Oscillation Amplitude

(1) Background: Together with treatment protocols, viscoelastic tests are widely used for patient care. Measuring at broader ranges of deformation than currently done will add information on a clot’s mechanical phenotype because fibrin networks follow different stretching regimes, and blood flow compels clots into a dynamic non-linear response. (2) Methods: To characterize the influence of platelets on the network level, a stress amplitude sweep test (LAOStress) was applied to clots from native plasma with five platelet concentrations. Five species were used to validate the protocol (human, cow, pig, rat, horse). By Lissajous plots the oscillation cycle for each stress level was analyzed. (3) Results: Cyclic stress loading generates a characteristic strain response that scales with the platelet quantity at low stress, and that is independent from the platelet count at high shear stress. This general behavior is valid in the animal models except cow. Here, the specific fibrinogen chemistry induces a stiffer network and a variant high stress response. (4) Conclusions: The protocol provides several thresholds to connect the softening and stiffening behavior of clots with the applied shear stress. This points to the reversible part of deformation, and thus opens a new route to describe a blood clot’s phenotype.

The effect of hematocrit, fibrinogen concentration and temperature on the kinetics of clot formation of whole blood

BACKGROUND

Dynamic mechanical analysis of blood clots can be used to detect the coagulability of blood.

OBJECTIVE

We investigated the kinetics of clot formation by changing several blood components, and we looked into the clot "signature" at its equilibrium state by using viscoelastic and dielectric protocols.

METHODS

Oscillating shear rheometry, ROTEM, and a dielectro-rheological device was used.

RESULTS

In fibrinogen- spiked samples we found the classical high clotting ability: shortened onset, faster rate of clotting, and higher plateau stiffness. Electron microscopy explained the gain of stiffness. Incorporated RBCs weakened the clots. Reduction of temperature during the clotting process supported the development of high moduli by providing more time for fiber assembly. But at low HCT, clot firmness could be increased by elevating the temperature from 32 to 37°C. In contrast, when the fibrinogen concentration was modified, acceleration of clotting via temperature always reduced clot stiffness, whatever the initial fibrinogen concentration. Electrical resistance increased continuously during clotting; loss tangent (D) (relaxation frequency 249 kHz) decreased when clots became denser: fewer dipoles contributed to the relaxation process. The relaxation peak (Dmax) shifted to lower frequencies at higher platelet count.

CONCLUSION

Increasing temperature accelerates clot formation but weakens clots. Rheometry and ROTEM correlate well.

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.

Development of a computational model for acute ischemic stroke recanalization through cyclic aspiration

Acute ischemic stroke (AIS), the result of embolic occlusion of a cerebral artery, is responsible for 87% of the 6.5 million stroke-related deaths each year. Despite improvements from first-generation thrombectomy devices for treating AIS, 80% of eligible stroke patients will either die or suffer a major disability. In order to maximize the number of patients with good outcomes, new AIS therapies need to be developed to achieve complete reperfusion on the first pass. One such therapy that has shown promise experimentally is the application of cyclic aspiration pressure, which led to higher recanalization rates at lower pressure magnitudes. In order to investigate AIS and cyclic aspiration recanalization, an improved computational modeling framework was developed, combining a viscoelastic thromboembolus model with a cohesive zone (CZ) model for the thromboembolus-artery interface. The model was first validated against experimental displacement data of a cyclically aspirated thromboembolus analog. The CZ model parameters, including the addition of a damage accumulation model, were then investigated computationally to determine their individual effects on the thromboembolus and CZ behavior. The relaxation time and the damage model critical opening length were shown to have the greatest effect on the CZ opening and led to increased displacement that accumulated with repeated loading. Additional simulations were performed with parameters relevant to AIS including internal carotid artery dimensions and thromboemboli mechanical properties. In these AIS cases, more upstream CZ opening was observed compared to the thromboembolus analog cases and greater displacement was achieved with the lower-frequency aspiration (0.5 vs 1 Hz).

Monitoring of blood coagulation with non-contact drop oscillation rheometry

BACKGROUND

Thromboelastography is widely used as a tool to assess the coagulation status of critical-care patients. It allows observation of changes in the material properties of whole blood brought about by clot formation and clot lysis. However, contact activation of the coagulation cascade at surfaces of thromboelastographic systems leads to inherent variability and unreliability in predicting bleeding or thrombosis risks, while also requiring large sample volumes.

OBJECTIVES

To develop a non-contact drop oscillation rheometry (DOR) method to measure the viscoelastic properties of blood clots and to compare the results with current laboratory standard measurements.

METHODS

Drops of human blood and plasma (5-10 μL) were acoustically levitated. Acoustic field modulation induced drop shape oscillations, and the viscoelastic properties of the sample were calculated by measuring the resonance frequency and damping ratio.

RESULTS

DOR showed sensitivity to coagulation parameters. An increase in platelet count resulted in an increase in the maximum clot stiffness. An increase in the calcium ion level enhanced the coagulation rate prior to saturation. An increase in hematocrit resulted in a higher rate of clot formation and increased clot stiffness. Comparison of the results with those obtained with thromboelastography showed that coagulation started sooner with DOR, but with a lower rate and lower maximum stiffness.

CONCLUSIONS

DOR can be used as a monitoring tool to assess blood coagulation status. The advantages of small sample size, the lack of contact and small strain (linear viscoelasticity) makes this technique unique for real-time monitoring of blood coagulation.

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.

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.

Modeling Ramp-hold Indentation Measurements based on Kelvin-Voigt Fractional Derivative Model

Interpretation of experimental data from micro- and nano-scale indentation testing is highly dependent on the constitutive model selected to relate measurements to mechanical properties. The Kelvin-Voigt Fractional Derivative model (KVFD) offers a compact set of viscoelastic features appropriate for characterizing soft biological materials. This paper provides a set of KVFD solutions for converting indentation testing data acquired for different geometries and scales into viscoelastic properties of soft materials. These solutions, which are mostly in closed-form, apply to ramp-hold relaxation, load-unload and ramp-load creep-testing protocols. We report on applications of these model solutions to macro- and nano-indentation testing of hydrogels, gastric cancer cells and ex vivo breast tissue samples using an Atomic Force Microscope (AFM). We also applied KVFD models to clinical ultrasonic breast data using a compression plate as required for elasticity imaging. Together the results show that KVFD models fit a broad range of experimental data with a correlation coefficient typically R2 > 0.99. For hydrogel samples, estimation of KVFD model parameters from test data using spherical indentation versus plate compression as well as ramp relaxation versus load-unload compression all agree within one standard deviation. Results from measurements made using macro- and nano-scale indentation agree in trend. For gastric cell and ex vivo breast tissue measurements, KVFD moduli are, respectively, 1/3 - 1/2 and 1/6 of the elasticity modulus found from the Sneddon model. In vivo breast tissue measurements yield model parameters consistent with literature results. The consistency of results found for a broad range of experimental parameters suggest the KVFD model is a reliable tool for exploring intrinsic features of the cell/tissue microenvironments.

Buffers Strongly Modulate Fibrin Self-Assembly into Fibrous Networks

Fibrin is a plasma protein with a central role in blood clotting and wound repair. Upon vascular injury, fibrin forms resilient fibrillar networks (clots) via a multistep self-assembly process, from monomers, to double-stranded protofibrils, to a branched network of thick fibers. In vitro, fibrin self-assembly is sensitive to physicochemical conditions like the solution pH and ionic strength, which tune the strength of the noncovalent driving forces. Here we report a surprising finding that the buffer-which is necessary to control the pH and is typically considered to be inert-also significantly influences fibrin self-assembly. We show by confocal microscopy and quantitative light scattering that various common buffering agents have no effect on the initial assembly of fibrin monomers into protofibrils but strongly hamper the subsequent lateral association of protofibrils into thicker fibers. We further find that the structural changes are independent of the molecular structure of the buffering agents as well as of the activation mechanism and even occur in fibrin networks formed from platelet-poor plasma. This buffer-mediated decrease in protofibril bundling results in a marked reduction in the permeability of fibrin networks but only weakly influences the elastic modulus of fibrin networks, providing a useful tuning parameter to independently control the elastic properties and the permeability of fibrin networks. Our work raises the possibility that fibrin assembly in vivo may be regulated by variations in the acute-phase levels of bicarbonate and phosphate, which act as physiological buffering agents of blood pH. Moreover, our findings add a new example of buffer-induced effects on biomolecular self-assembly to recent findings for a range of proteins and lipids.

Assessing composition and structure of soft biphasic media from Kelvin-Voigt fractional derivative model parameters

Kelvin-Voigt fractional derivative (KVFD) model parameters have been used to describe viscoelastic properties of soft tissues. However, translating model parameters into a concise set of intrinsic mechanical properties related to tissue composition and structure remains challenging. This paper begins by exploring these relationships using a biphasic emulsion materials with known composition. Mechanical properties are measured by analyzing data from two indentation techniques - ramp-stress relaxation and load-unload hysteresis tests. Material composition is predictably correlated with viscoelastic model parameters. Model parameters estimated from the tests reveal that elastic modulus E₀ closely approximates the shear modulus for pure gelatin. Fractional-order parameter α and time constant τ vary monotonically with the volume fraction of the material's fluid component. α characterizes medium fluidity and the rate of energy dissipation, and τ is a viscous time constant. Numerical simulations suggest that the viscous coefficient η is proportional to the energy lost during quasi-static force-displacement cycles, E_A . The slope of E_A versus η is determined by α and the applied indentation ramp time T_r. Experimental measurements from phantom and ex vivo liver data show close agreement with theoretical predictions of the η - E_A relation. The relative error is less than 20% for emulsions 22% for liver. We find that KVFD model parameters form a concise features space for biphasic medium characterization that described time-varying mechanical properties.

Fibrin mechanical properties and their structural origins

Fibrin is a protein polymer that is essential for hemostasis and thrombosis, wound healing, and several other biological functions and pathological conditions that involve extracellular matrix. In addition to molecular and cellular interactions, fibrin mechanics has been recently shown to underlie clot behavior in the highly dynamic intra- and extravascular environments. Fibrin has both elastic and viscous properties. Perhaps the most remarkable rheological feature of the fibrin network is an extremely high elasticity and stability despite very low protein content. Another important mechanical property that is common to many filamentous protein polymers but not other polymers is stiffening occurring in response to shear, tension, or compression. New data has begun to provide a structural basis for the unique mechanical behavior of fibrin that originates from its complex multi-scale hierarchical structure. The mechanical behavior of the whole fibrin gel is governed largely by the properties of single fibers and their ensembles, including changes in fiber orientation, stretching, bending, and buckling. The properties of individual fibrin fibers are determined by the number and packing arrangements of double-stranded half-staggered protofibrils, which still remain poorly understood. It has also been proposed that forced unfolding of sub-molecular structures, including elongation of flexible and relatively unstructured portions of fibrin molecules, can contribute to fibrin deformations. In spite of a great increase in our knowledge of the structural mechanics of fibrin, much about the mechanisms of fibrin's biological functions remains unknown. Fibrin deformability is not only an essential part of the biomechanics of hemostasis and thrombosis, but also a rapidly developing field of bioengineering that uses fibrin as a versatile biomaterial with exceptional and tunable biochemical and mechanical properties.