The review of contemporary ideas about the influence of flow rate on blood clotting

Normal blood clotting is vitally important for mammals. The diffusion-convection transfer of clotting factors plays a key role in blood clot formation. Since the shear rates of blood flow are very high (up to 7000 s⁻¹), clot formation critically depends on the flow rate. The methods of study of the flow effect on clotting are indirect and the processes are rather complex; therefore, mathematical models of this process are significant for interpretation of results and understanding of the mechanisms. The review expounds the main experimental data on the effect of flow on the blood clotting cascade, some hypotheses and mathematical models explaining different regimes of the functioning of this system. The review is focused on specific problems encountered by researchers in this field. Some of the experimental works have shown that flow decreases the size of the formed fibrin clot and that the dependence of clot formation period on the flow shear rate is a threshold function. However, there are also experimental evidence that the flow can increase production of clotting factors (factor Xa), which must be expressed in the procoagulant action of the flow. Mathematical models of different aspects of clotting give no unified predictions either. Nevertheless, the combined analysis of results of detailed modeling and experiments, in our opinion, may result in understanding of the mechanisms, which determine the behavior of clotting in a flow.     

The Influence of Hindered Transport on the Development of Platelet Thrombi Under Flow

Vascular injury triggers two intertwined processes, platelet deposition and coagulation, and can lead to the formation of an intravascular clot (thrombus) that may grow to occlude the vessel. Formation of the thrombus involves complex biochemical, biophysical, and biomechanical interactions that are also dynamic and spatially-distributed, and occur on multiple spatial and temporal scales. We previously developed a spatial-temporal mathematical model of these interactions and looked at the interplay between physical factors (flow, transport to the clot, platelet distribution within the blood) and biochemical ones in determining the growth of the clot. Here, we extend this model to include reduction of the advection and diffusion of the coagulation proteins in regions of the clot with high platelet number density. The effect of this reduction, in conjunction with limitations on fluid and platelet transport through dense regions of the clot can be profound. We found that hindered transport leads to the formation of smaller and denser clots compared to the case with no protein hindrance. The limitation on protein transport confines the important activating complexes to small regions in the interior of the thrombus and greatly reduces the supply of substrates to these complexes. Ultimately, this decreases the rate and amount of thrombin production and leads to greatly slowed growth and smaller thrombus size. Our results suggest a possible physical mechanism for limiting thrombus growth.     

A computational thrombus formation model: application to an idealized two-dimensional aneurysm treated with bare metal coils

Cardiovascular implantable devices alter the biofluid dynamics and biochemistry of the blood in which they are placed. These perturbations can lead to thrombus formation which may or may not be desired, depending on the application. In this work, a computational model is developed that couples biofluid dynamics and biochemistry to predict the clotting response of blood to such devices. The model consists of 28 advection–diffusion–reaction partial differential equations to track proteins in the blood involved in clotting and utilizes boundary flux terms to model the initiation of the intrinsic clotting pathway at thrombogenic device surfaces. We use this model to simulate the transient clot growth within a 2D idealized bifurcation aneurysm filled with various distributions of bare metal coils with similar packing densities. The clot model predicts initial clot formation to occur in areas along coil surfaces where flow is minimal and where time-averaged shear rates are the smallest. Among the six coil-filled aneurysm cases simulated, maximum thrombus occlusion ranged between 80.8 and 92.2% of the post-treatment aneurysm volume and was achieved 325–450 s after treatment. With further refinement and validation, the computational clotting model will be a valuable engineering tool for evaluating and comparing the relative performance of cardiovascular implantable devices.     

Impact of Tissue Factor Localization on Blood Clot Structure and Resistance under Venous Shear

The structure and growth of a blood clot depend on the localization of tissue factor (TF), which can trigger clotting during the hemostatic process or promote thrombosis when exposed to blood under pathological conditions. We sought to understand how the growth, structure, and mechanical properties of clots under flow are shaped by the simultaneously varying TF surface density and its exposure area. We used an eight-channel microfluidic device equipped with a 20- or 100-μm-long collagen surface patterned with lipidated TF of surface densities ～0.1 and ～2 molecules/μm². Human whole blood was perfused at venous shear, and clot growth was continually measured. Using our recently developed computational model of clot formation, we performed simulations to gain insights into the clot’s structure and its resistance to blood flow. An increase in TF exposure area resulted not only in accelerated bulk platelet, thrombin, and fibrin accumulation, but also in increased height of the platelet mass and increased clot resistance to flow. Moreover, increasing the TF surface density or exposure area enhanced platelet deposition by approximately twofold, and thrombin and fibrin generation by greater than threefold, thereby increasing both clot size and its viscous resistance. Finally, TF effects on blood flow occlusion were more pronounced for the longer thrombogenic surface than for the shorter one. Our results suggest that TF surface density and its exposure area can independently enhance both the clot’s occlusivity and its resistance to blood flow. These findings provide, to our knowledge, new insights into how TF affects thrombus growth in time and space under flow.     

The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions

In the abdominal segment of the human aorta under a patient's average resting conditions, pulsatile blood flow exhibits complex laminar patterns with secondary flows induced by adjacent branches and irregular vessel geometries. The flow dynamics becomes more complex when there is a pathological condition that causes changes in the normal structural composition of the vessel wall, for example, in the presence of an aneurysm. This work examines the hemodynamics of pulsatile blood flow in hypothetical three-dimensional models of abdominal aortic aneurysms (AAAs). Numerical predictions of blood flow patterns and hemodynamic stresses in AAAs are performed in single-aneurysm, asymmetric, rigid wall models using the finite element method. We characterize pulsatile flow dynamics in AAAs for average resting conditions by means of identifying regions of disturbed flow and quantifying the disturbance by evaluating flow-induced stresses at the aneurysm wall, specifically wall pressure and wall shear stress. Physiologically realistic abdominal aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 < or = Rem < or = 300, corresponding to a range of peak Reynolds numbers 262.5 < or = Repeak < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is depicted by a sequence of four different flow phases in one period of the cardiac pulse. Peak wall shear stress and peak wall pressure are reported as a function of the time-average Reynolds number and aneurysm asymmetry. The effect of asymmetry in hypothetically shaped AAAs is to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole.     

Biophysical mechanisms of contact activation of blood-plasma clotting

This review considers the biochemical and biophysical mechanisms that trigger blood clotting upon contact of blood with an alien surface and leads via a cascade of enzymatic reactions to fibrin polymerization and the formation of a blood plasma clot, which permeates a primary platelet aggregate to produce a dense hemostatic clot. In spite of the substantial number of experimental and theoretical studies on the subject, there is still no consistent opinion as to what processes occur as the blood plasma contacts a surface. This review discusses the role that plasma protein factor XII and various surfaces play in triggering the contact pathway in vivo and in vitro. Current views of the molecular events that underlie the process are described.     

Flow-induced wall shear stress in abdominal aortic aneurysms: Part II--pulsatile flow hemodynamics

In continuing the investigation of AAA hemodynamics, unsteady flow-induced stresses are presented for pulsatile blood flow through the double-aneurysm model described in Part I. Physiologically realistic aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50< or =Re(m) < or =300. Hemodynamic disturbance is evaluated for a modified set of indicator functions which include wall pressure (p(w)), wall shear stress (tau(w)), Wall Shear Stress Gradient (WSSG), time-average wall shear stress (tau(w)*), and time-average Wall Shear Stress Gradient WSSG*. At peak flow, the highest shear stress and WSSG levels are obtained at the distal end of both aneurysms, in a pattern similar to that of steady flow. The maximum values of wall shear stresses and wall shear stress gradients are evaluated as a function of the time-average Reynolds number resulting in a fourth order polynomial correlation. A comparison between numerical predictions for steady and pulsatile flow is presented, illustrating the importance of considering time-dependent flow for the evaluation of hemodynamic indicators.     

Flow-induced Wall Shear Stress in Abdominal Aortic Aneurysms: Part II - Pulsatile Flow Hemodynamics

In continuing the investigation of AAA hemodynamics, unsteady flow-induced stresses are presented for pulsatile blood flow through the double-aneurysm model described in Part I. Physiologically realistic aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 h Re m h 300. Hemodynamic disturbance is evaluated for a modified set of indicator functions which include wall pressure ( p w ), wall shear stress ( w ), Wall Shear Stress Gradient (WSSG), time-average wall shear stress ( w *), and time-average Wall Shear Stress Gradient WSSG *. At peak flow, the highest shear stress and WSSG levels are obtained at the distal end of both aneurysms, in a pattern similar to that of steady flow. The maximum values of wall shear stresses and wall shear stress gradients are evaluated as a function of the time-average Reynolds number resulting in a fourth order polynomial correlation. A comparison between numerical predictions for steady and pulsatile flow is presented, illustrating the importance of considering time-dependent flow for the evaluation of hemodynamic indicators.     

Numerical modeling of pulsatile turbulent flow in stenotic vessels

Pulsatile turbulent flow in stenotic vessels has been numerically modeled using the Reynolds-averaged Navier-Stokes equation approach. The commercially available computational fluid dynamics code (CFD), FLUENT, has been used for these studies. Two different experiments were modeled involving pulsatile flow through axisymmetric stenoses. Four different turbulence models were employed to study their influence on the results. It was found that the low Reynolds number k-omega turbulence model was in much better agreement with previous experimental measurements than both the low and high Reynolds number versions of the RNG (renormalization-group theory) k-epsilon turbulence model and the standard k-epsilon model, with regard to predicting the mean flow distal to the stenosis including aspects of the vortex shedding process and the turbulent flow field. All models predicted a wall shear stress peak at the throat of the stenosis with minimum values observed distal to the stenosis where flow separation occurred.     

Modeling pulsatile flow in aortic aneurysms: effect of non-Newtonian properties of blood

Pulsatile flow in an axisymmetric rigid-walled model of an abdominal aorta aneurysm was analyzed numerically for various aneurysm dilations using physiologically realistic resting waveform at time-averaged Reynolds number of 300 and peak Reynolds number of 1607. Discretization of the governing equations was achieved using a finite element scheme based on the Galerkin method of weighted residuals. Comparisons with previously published work on the basis of special cases were performed and found to be in excellent agreement. Our findings indicate that the velocity fields are significantly affected by non-Newtonian properties in pathologically altered configurations. Non-Newtonian fluid shear stress is found to be greater than Newtonian fluid shear stress during peak systole. Further, the maximum shear stress is found to occur near the distal end of AAA during peak systole. The impact of non-Newtonian blood flow characteristics on pressure compared to Newtonian model is found insignificant under resting conditions. Viscous and inertial forces associated with blood flow are responsible for the changes in the wall that result in thrombus deposition and dilation while rupture of AAA is more likely determined by much larger mechanical stresses imposed by pulsatile pressure on the wall of AAA.     

Viscous flow simulation in a stenosis model using discrete particle dynamics: a comparison between DPD and CFD

Flow and stresses induced by blood flow acting on the blood cellular constituents can be represented to a certain extent by a continuum mechanics approach down to the order of the?μm level. However, the molecular effects of, e.g., adhesion/aggregation bonds of blood clotting can be on the order of nm. The coupling of the disparate length and timescales between such molecular levels and macroscopic transport represents a major computational challenge. To address this challenge, a multiscale numerical approach based on discrete particle dynamics (DPD) methodology derived from molecular dynamics (MD) principles is proposed. The feasibility of the approach was firstly tested for its ability to simulate viscous flow conditions. Simulations were conducted in low Reynolds numbers flows (Re = 25–33) through constricted tubes representing blood vessels with various degrees of stenosis. Multiple discrete particles interacting with each other were simulated, with 1.24–1.36 million particles representing the flow domain and 0.4 million particles representing the vessel wall. The computation was carried out on the massive parallel supercomputer NY BlueGene/L employing NAMD-a parallel MD package for high performance computing (HPC). Typical recirculation zones were formed distal to the stenoses. The velocity profiles and recirculation zones were in excellent agreement with computational fluid dynamics (CFD) 3D Navier–Stokes viscous fluid flow simulations and with classic numerical and experimental results by YC Fung in constricted tubes. This feasibility analysis demonstrates the potential of a methodology that widely departs from a continuum approach to simulate multiscale phenomena such as flow induced blood clotting.     

A mathematical model for the dissolution of non-occlusive blood clots in fast tangential blood flow

Our aim was to study the effect of an axially directed blood plasma flow on the dissolution rate of cylindrical non-occlusive blood clots in an in vitro flow system and to derive a mathematical model for the process. The model was based on the hypothesis that clot dissolution dynamics is proportional not only to the biochemical proteolysis of fibrin but also to the power of the flowing blood plasma dissipated along the clot. The predicted rate of thrombolysis is then proportional to the square of the average blood plasma velocity for laminar flow and to the third power of the average velocity for turbulent flow. To verify the model, the time dependence of the clot cross-sectional area was measured by dynamic magnetic resonance microscopy during fast (turbulent) and slow (laminar) flow of plasma through an axially directed channel along the clot. The flowing plasma contained a magnetic resonance imaging contrast agent (Gd-DTPA) and a thrombolytic agent (recombinant tissue-type plasminogen activator). The experimental data fitted well to the model, and confirmed the predicted increase in the dissolution rate when blood flow changed from a laminar to a turbulent flow regime.     

Modelling the impact of clot fragmentation on the microcirculation after thrombectomy

Many ischaemic stroke patients who have a mechanical removal of their clot (thrombectomy) do not get reperfusion of tissue despite the thrombus being removed. One hypothesis for this ‘no-reperfusion’ phenomenon is micro-emboli fragmenting off the large clot during thrombectomy and occluding smaller blood vessels downstream of the clot location. This is impossible to observe in-vivo and so we here develop an in-silico model based on in-vitro experiments to model the effect of micro-emboli on brain tissue. Through in-vitro experiments we obtain, under a variety of clot consistencies and thrombectomy techniques, micro-emboli distributions post-thrombectomy. Blood flow through the microcirculation is modelled for statistically accurate voxels of brain microvasculature including penetrating arterioles and capillary beds. A novel micro-emboli algorithm, informed by the experimental data, is used to simulate the impact of micro-emboli successively entering the penetrating arterioles and the capillary bed. Scaled-up blood flow parameters–permeability and coupling coefficients–are calculated under various conditions. We find that capillary beds are more susceptible to occlusions than the penetrating arterioles with a 4x greater drop in permeability per volume of vessel occluded. Individual microvascular geometries determine robustness to micro-emboli. Hard clot fragmentation leads to larger micro-emboli and larger drops in blood flow for a given number of micro-emboli. Thrombectomy technique has a large impact on clot fragmentation and hence occlusions in the microvasculature. As such, in-silico modelling of mechanical thrombectomy predicts that clot specific factors, interventional technique, and microvascular geometry strongly influence reperfusion of the brain. Micro-emboli are likely contributory to the phenomenon of no-reperfusion following successful removal of a major clot.     

A Constitutive Model for a Maturing Fibrin Network

Blood clot formation is crucial to maintain normal physiological conditions but at the same time involved in many diseases. The mechanical properties of the blood clot are important for its functioning but complicated due to the many processes involved. The main structural component of the blood clot is fibrin, a fibrous network that forms within the blood clot, thereby increasing its mechanical rigidity. A constitutive model for the maturing fibrin network is developed that captures the evolving mechanical properties. The model describes the fibrin network as a network of fibers that become thicker in time. Model parameters are related to the structural properties of the network, being the fiber length, bending stiffness, and mass-length ratio. Results are compared with rheometry experiments in which the network maturation is followed in time for various loading frequencies and fibrinogen concentrations. Three parameters are used to capture the mechanical behavior including the mass-length ratio. This parameter agrees with values determined using turbidimetry experiments and is subsequently used to derive the number of protofibrils and fiber radius. The strength of the model is that it describes the mechanical properties of the maturing fibrin network based on it structural quantities. At the same time the model is relatively simple, which makes it suitable for advanced numerical simulations of blood clot formation during flow in blood vessels.     

A Constitutive Model for a Maturing Fibrin Network

Blood clot formation is crucial to maintain normal physiological conditions but at the same time involved in many diseases. The mechanical properties of the blood clot are important for its functioning but complicated due to the many processes involved. The main structural component of the blood clot is fibrin, a fibrous network that forms within the blood clot, thereby increasing its mechanical rigidity. A constitutive model for the maturing fibrin network is developed that captures the evolving mechanical properties. The model describes the fibrin network as a network of fibers that become thicker in time. Model parameters are related to the structural properties of the network, being the fiber length, bending stiffness, and mass-length ratio. Results are compared with rheometry experiments in which the network maturation is followed in time for various loading frequencies and fibrinogen concentrations. Three parameters are used to capture the mechanical behavior including the mass-length ratio. This parameter agrees with values determined using turbidimetry experiments and is subsequently used to derive the number of protofibrils and fiber radius. The strength of the model is that it describes the mechanical properties of the maturing fibrin network based on it structural quantities. At the same time the model is relatively simple, which makes it suitable for advanced numerical simulations of blood clot formation during flow in blood vessels.     

Flow-Induced β-Hairpin Folding of the Glycoprotein Ibα β-Switch

Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces β-hairpin folding of the glycoprotein Ibα β-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-μs molecular dynamics simulations, we investigate the kinetics of flow-induced β-hairpin folding. Simulations sampling different flow velocities reveal that under flow, β-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded β-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the β-switch has two stable conformations imprinted on it, namely, loop and hairpin—with flow inducing a transition between the two.     

Microfluidic and computational study of structural properties and resistance to flow of blood clots under arterial shear

The ability of a blood clot to modulate blood flow is determined by the clot’s resistance, which depends on its structural features. For a flow with arterial shear, we investigated the characteristic patterns relating to clot shape, size, and composition on the one hand, and its viscous resistance, intraclot axial flow velocity, and shear distributions on the other. We used microfluidic technology to measure the kinetics of platelet, thrombin, and fibrin accumulation at a thrombogenic surface coated with collagen and tissue factor (TF), the key clot-formation trigger. We subsequently utilized the obtained data to perform additional calibration and validation of a detailed computational fluid dynamics model of spatial clot growth under flow. We then ran model simulations to gain insights into the resistance of clots formed under our experimental conditions. We found that increased thrombogenic surface length and TF surface density enhanced the bulk thrombin and fibrin generation in a nonadditive, synergistic way. The height of the platelet deposition domain—and, therefore, clot occlusivity—was rather robust to thrombogenic surface length and TF density variations, but consistently increased with time. Clot viscous resistance was non-uniform and tended to be higher in the fibrin-rich, inner “core” region of the clot. Interestingly, despite intraclot structure and viscous resistance variations, intraclot flow velocity variations were minor compared to the abrupt decrease in flow velocity around the platelet deposition region. Our results shed new light on the connection between the structure of clots under arterial shear and spatiotemporal variations in their resistance to flow.

     

Maintenance of silastic–teflon shunts for intermittent haemodialysis

The occurrence of infection in the tissues surrounding external arteriovenous shunts was studied and die important relationship of pyogenic infection to clotting was confirmed. The local application of fusidic add tulle and lanolin greatly reduced the occurrence of both infection and clotting and the need for cannula replacement.Urokinase used for declotting shunts when standard procedures had failed, restored blood flow whether dotting was related to infection or to local vascular factors. This treatment is not advised when clotting is associated with a local abscess, as it may make cannula replacement necessary. Severe local vascular factors, such as metastatic calcification, Raynaud''s phenomenon, and venous stenosis, may lead to poor blood flow, so that despite clot lysis elective cannula replacement or the creation of a subcutaneous arteriovenous fistula is required.     

Hemostasis as an optimal system

Zymogen and procofactor concentrations in physiological biochemical systems (PBS) have not yet been explained. The problem in question is to determine optimal plasma clotting factor (factors II, VII, IX, and X and cofactors V and VIII) concentrations for coagulation system (CS) as a whole. Constrained optimization technique is used to solve this problem. The constraint is determined by the ability of the CS to perform its physiological function--thrombin generation (and hence clot formation)--under vessel injury conditions. The constraint statement is based on the CS dynamics equations. In solving the problem the Lagrange multiplier is used. A hypothesis is advanced that this problem can be solved based on principle of minimum protein consumption subject to an above constraint. The results obtained indicate that the optimal clotting factor concentrations are in good agreement with those measured by biochemical techniques. A comparison between the theoretical results and experimental data lends support for our hypothesis that zymogen and procofactor concentrations in the CS (and, probably, for other biochemical systems) are determined by the principle of minimum protein consumption.