Blood flow and cell-free layer in microvessels

Imaging single-cell blood flow in the smallest to largest vessels in the living retina

Tissue light scatter limits the visualization of the microvascular network deep inside the living mammal. The transparency of the mammalian eye provides a noninvasive view of the microvessels of the retina, a part of the central nervous system. Despite its clarity, imperfections in the optics of the eye blur microscopic retinal capillaries, and single blood cells flowing within. This limits early evaluation of microvascular diseases that originate in capillaries. To break this barrier, we use 15 kHz adaptive optics imaging to noninvasively measure single-cell blood flow, in one of the most widely used research animals: the C57BL/6J mouse. Measured flow ranged four orders of magnitude (0.0002-1.55 µL min^-1) across the full spectrum of retinal vessel diameters (3.2-45.8 µm), without requiring surgery or contrast dye. Here, we describe the ultrafast imaging, analysis pipeline and automated measurement of millions of blood cell speeds.

In vitro and ex vivo measurement of the biophysical properties of blood using microfluidic platforms and animal models

Haemorheologically impaired microcirculation, such as blood clotting or abnormal blood flow, causes interrupted blood flows in vascular networks. The biophysical properties of blood, including blood viscosity, blood viscoelasticity, haematocrit, red blood bell (RBC) aggregation, erythrocyte sedimentation rate and RBC deformability, have been used to monitor haematological diseases. In this review, we summarise several techniques for measuring haemorheological properties, such as blood viscosity, RBC deformability and RBC aggregation, using in vitro microfluidic platforms. Several methodologies for the measurement of haemorheological properties with the assistance of an extracorporeal rat bypass loop are also presented. We briefly discuss several emerging technologies for continuous, long-term, multiple measurements of haemorheological properties under in vitro or ex vivo conditions.

Migration velocity of red blood cells in microchannels

The lateral migration of red blood cells (RBCs) in confined channel flows is an important ingredient of microcirculatory hydrodynamics and is involved in the development of a cell free layer near vessel walls and influences the distribution of RBCs in networks. It is also relevant to a number of lab-on-chip applications. This migration is a consequence of their deformability and is due to the combined effects of hydrodynamic wall repulsion and the curvature of the fluid velocity profile. We performed microfluidic experiments with dilute suspensions of RBCs in which the trajectories and migration away from the channel wall are analyzed to extract the mean behavior, from which we propose a generic scaling law for the transverse migration velocity valid in a whole range of parameters relevant to microcirculatory and practical situations. Experiments with RBCs of different mechanical properties (separated by density gradient sedimentation or fixed with glutaraldehyde) show the influence of this parameter which can induce significant dispersion of the trajectories.

A computational model for microcirculation including Fahraeus-Lindqvist effect, plasma skimming and fluid exchange with the tissue interstitium

We present a two-phase model for microcirculation that describes the interaction of plasma with red blood cells. The model takes into account of typical effects characterizing the microcirculation, such as the Fahraeus-Lindqvist effect and plasma skimming. Besides these features, the model describes the interaction of capillaries with the surrounding tissue. More precisely, the model accounts for the interaction of capillary transmural flow with the surrounding interstitial pressure. Furthermore, the capillaries are represented as one-dimensional channels with arbitrary, possibly curved configuration. The latter two features rely on the unique ability of the model to account for variations of flow rate and pressure along the axis of the capillary, according to a local differential formulation of mass and momentum conservation. Indeed, the model stands on a solid mathematical foundation, which is also addressed in this work. In particular, we present the model derivation, the variational formulation, and its approximation using the finite element method. Finally, we conclude the work with a comparative computational study of the importance of the Fahraeus-Lindqvist, plasma skimming, and capillary leakage effects on the distribution of flow in a microvascular network.

Concurrent coupling of atomistic simulation and mesoscopic hydrodynamics for flows over soft multi-functional surfaces

We develop an efficient parallel multiscale method that bridges the atomistic and mesoscale regimes, from nanometers to microns and beyond, via concurrent coupling of atomistic simulation and mesoscopic dynamics. In particular, we combine an all-atom molecular dynamics (MD) description for specific atomistic details in the vicinity of the functional surface with a dissipative particle dynamics (DPD) approach that captures mesoscopic hydrodynamics in the domain away from the functional surface. In order to achieve a seamless transition in dynamic properties we endow the MD simulation with a DPD thermostat, which is validated against experimental results by modeling water at different temperatures. We then validate the MD-DPD coupling method for transient Couette and Poiseuille flows, demonstrating that the concurrent MD-DPD coupling can resolve accurately the continuum-based analytical solutions. Subsequently, we simulate shear flows over grafted polydimethylsiloxane (PDMS) surfaces (polymer brushes) for various grafting densities, and investigate the slip flow as a function of the shear stress. We verify that a "universal" power law exists for the slip length, in agreement with published results. Having validated the MD-DPD coupling method, we simulate time-dependent flows past an endothelial glycocalyx layer (EGL) in a microchannel. Coupled simulation results elucidate the dynamics of the EGL changing from an equilibrium state to a compressed state under shear by aligning the molecular structures along the shear direction. MD-DPD simulation results agree well with results of a single MD simulation, but with the former more than two orders of magnitude faster than the latter for system sizes above one micron.

In-line whole blood fractionation for Raman analysis of blood plasma

Raman studies of dynamically expanded cell-free plasma domains in microfluidic blood flow.

A model of oxygen transport in the rat renal medulla

The renal medulla is prone to hypoxia. Medullary hypoxia is postulated to be a leading cause of acute kidney injury, so there is considerable interest in predicting the oxygen tension in the medulla. Therefore we have developed a computational model for blood and oxygen transport within a physiologically normal rat renal medulla, using a multilevel modeling approach. For the top-level model we use the theory of porous media and advection-dispersion transport through a realistic three-dimensional representation of the medulla’s gross anatomy to describe blood flow and oxygen transport throughout the renal medulla. For the lower-level models, we employ two-dimensional reaction-diffusion models describing the distribution of oxygen through tissue surrounding the vasculature. Steady-state model predictions at the two levels are satisfied simultaneously, through iteration between the levels. The computational model was validated by simulating eight sets of experimental data regarding renal oxygenation in rats (using 4 sets of control groups and 4 sets of treatment groups, described in 4 independent publications). Predicted medullary tissue oxygen tension or microvascular oxygen tension for control groups and for treatment groups that underwent moderate perturbation in hemodynamic and renal functions is within ±2 SE values observed experimentally. Diffusive shunting between descending and ascending vasa recta is predicted to be only 3% of the oxygen delivered. The validation tests confirm that the computational model is robust and capable of capturing the behavior of renal medullary oxygenation in both normal and early-stage pathological states in the rat.

Modeling Cell Adhesion and Extravasation in Microvascular System

The blood flow behaviors in the microvessels determine the transport modes and further affect the metastasis of circulating tumor cells (CTCs). Much biochemical and biological efforts have been made on CTC metastasis; however, precise experimental measurement and accurate theoretical prediction on its mechanical mechanism are limited. To complement these, numerical modeling of a CTC extravasation from the blood circulation, including the steps of adhesion and transmigration, is discussed in this chapter. The results demonstrate that CTCs prefer to adhere at the positive curvature of curved microvessels, which is attributed to the positive wall shear stress/gradient. Then, the effects of particulate nature of blood on CTC adhesion are investigated and are found to be significant in the microvessels. Furthermore, the presence of red blood cell (RBC) aggregates is also found to promote the CTC adhesion by providing an additional wall-directed force. Finally, a single cell passing through a narrow slit, mimicking CTC transmigration, was examined under the effects of cell deformability. It showed that the cell shape and surface area increase play a more important role than the cell elasticity in cell transit across the narrow slit.

Quantifying Platelet Margination in Diabetic Blood Flow

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

Shear rate dependent margination of sphere-like, oblate-like and prolate-like micro-particles within blood flow

This study investigates the shear rate dependent margination of micro-particles (MPs) with different shapes in blood flow through numerical simulations. We develop a multiscale computational model to handle the fluid-structure interactions involved in the blood flow simulations. The lattice Boltzmann method (LBM) is used to solve the plasma dynamics and a coarse-grained model is employed to capture the dynamics of red blood cells (RBCs) and MPs. These two solvers are coupled together by the immersed boundary method (IBM). The shear rate dependent margination of sphere MPs is firstly investigated. We find that margination of sphere MPs dramatically increases with the increment of wall shear rate [small gamma, Greek, dot above]ω under 800 s-1, induced by the breaking of rouleaux in blood flow. However, the margination probability only slowly grows when [small gamma, Greek, dot above]ω > 800 s-1. Furthermore, the shape effect of MPs is examined by comparing the margination behaviors of sphere-like, oblate-like and prolate-like MPs under different wall shear rates. We find that the margination of MPs is governed by the interplay of two factors: hydrodynamic collisions with RBCs including the collision frequency and collision displacement of MPs, and near wall dynamics. MPs that demonstrate poor performance in one process such as collision frequency may stand out in the other process like near wall dynamics. Specifically, the ellipsoidal MPs (oblate and prolate) with small aspect ratio (AR) outperform those with large AR regardless of the wall shear rate, due to their better performance in both the collision with RBCs and near wall dynamics. Additionally, we find there exists a transition shear rate region 700 s-1 < [small gamma, Greek, dot above]ω < 900 s-1 for all of these MPs: the margination probability dramatically increases with the shear rate below this region and slowly grows above this region, similar to sphere MPs. We further use the surface area to volume ratio (SVR) to distinguish different shaped MPs and illustrate their shear rate dependent margination in a contour in the shear rate-SVR plane. It is of significance that we can approximately predict the margination of MPs with a specific SVR. All these simulation results can be potentially applied to guide the design of micro-drug carriers for biomedical applications.

Dynamics of blood flow: modeling of Fåhraeus and Fåhraeus-Lindqvist effects using a shear-induced red blood cell migration model

Blood flow in micro capillaries of diameter approximately 15-500 μm is accompanied with a lower tube hematocrit level and lower apparent viscosity as the diameter decreases. These effects are termed the Fåhraeus and Fåhraeus-Lindqvist effects, respectively. Both effects are linked to axial accumulation of red blood cells. In the present investigation, we extend previous works using a shear-induced model for the migration of red blood cells and adopt a model for blood viscosity that accounts for the suspending medium viscosity and local hematocrit level. For fully developed hematocrit profiles (i.e., independent of axial location), the diffusion fluxes due to particle collision frequency and viscosity gradients are of equal magnitude and opposite directions. The ratio of the diffusion coefficients for the two fluxes affects both the Fåhraeus and Fåhraeus-Lindqvist effects and is found related to the capillary diameter and discharge hematocrit using a well-known data-fit correlation for apparent blood viscosity. The velocity and hematocrit profiles were determined numerically as functions of radial coordinate, tube diameter, and discharge hematocrit. The velocity profile determined numerically is consistent with the derived analytical expression and the results are in good agreement with published numerical results and experimental data for hematocrit ratio and hematocrit and velocity profiles.

The dynamics of a healthy and infected red blood cell in flow through constricted channels: A DPD simulation

Understanding the dynamics of red blood cell (RBC) motion under in silico conditions is central to the development of cost-effective diagnostic tools. Specifically, unraveling the relationship between the rheological properties and the nature of shape change in the RBC (healthy or infected) can be extremely useful. In case of malarial infection, RBC progressively loses its deformability and tends to occlude the microvessel. In the present study, detailed mesoscopic simulations are performed to investigate the deformation dynamics of an RBC in flow through a constricted channel. Specifically, the manifestation of viscous forces (through flow rates) on the passage and blockage characteristics of a healthy red blood cell (hRBC) vis-á-vis an infected red blood cell (iRBC) are investigated. A finite-sized dissipative particle dynamics framework is used to model plasma in conjunction with a discrete model for the RBC. Instantaneous wall boundary method was used to model no-slip wall boundary conditions with a good control on the near-wall density fluctuations and compressibility effects. To investigate the microvascular occlusion, the RBC motion through 2 types of constricted channels, viz, (1) a tapered microchannel and (2) a stenosed-type microchannel, were simulated. It was observed that the deformation of an infected cell was much less compared with a healthy cell, with an attendant increase in the passage time. Apart from the qualitative features, deformation indices were obtained. The deformation of hRBC was sudden, while the iRBC deformed slowly as it traversed through the constriction. For higher flow rates, both hRBC and iRBC were found to undergo severe deformation. Even under low flow rates, hRBC could easily traverse past the constricted channel. However, for sufficiently slow flow rates (eg, capillary flows), the microchannel was found to be completely blocked by the iRBC.

Design and Utility of a Point-of-Care Microfluidic Platform to Assess Hematocrit and Blood Coagulation

Purpose—

To develop a small volume whole blood analyzer capable of measuring the hematocrit and coagulation kinetics of whole blood.

Methods and Results—

A co-planar microfluidic chamber designed to facilitate self-driven capillary action across an internal electrical chip was developed and used to measure the electric parameters of whole human blood that had been anticoagulated or allowed to clot. To promote blood clotting, select chip surfaces were coated with a prothrombin time (PT) reagent containing lipidated tissue factor (TF), which activates the extrinsic pathway of coagulation to promote thrombin generation and fibrin formation. Whole human blood was added to the microfluidic device, and voltage changes within the platform were measured and interpreted using basic resistor-capacitor (RC) circuit and fluid dynamics theory. Upon wetting of the sensing zone, a circuit between two co-planar electrodes within the sensing zone was closed to generate a rapid voltage drop from baseline. The voltage then rose due to sedimentation of red blood cells (RBC) in the sensing zone. For anticoagulated blood samples, the time for the voltage to return to baseline was dependent on hematocrit. In the presence of coagulation, the initiation of fibrin formation in the presence of the PT reagent prevented the return of voltage to baseline due to the reduced packing of RBCs in the sensing zone.

Conclusions—

The technology presented in this study has potential for monitoring the hematocrit and coagulation parameters of patient samples using a small volume of whole blood, suggesting it may hold clinical utility as a point-of-care test.

Local Hematocrit Fluctuation Induced by Malaria-Infected Red Blood Cells and Its Effect on Microflow

We investigate numerically the microscale blood flow in which red blood cells (RBCs) are partially infected by Plasmodium falciparum, the malaria parasite. The infected RBCs are modeled as more rigid cells with less deformability than healthy ones. Our study illustrates that, in a 10 μm microvessel in low-hematocrit conditions (18% and 27%), the Plasmodium falciparum-infected red blood cells (Pf-IRBCs) and healthy ones first form a train of cells. Because of the slow moving of the Pf-IRBCs, the local hematocrit (H_ct) near the Pf-IRBCs is then increased, to approximately 40% or even higher values. This increase of the local hematocrit is temporary and is kept for a longer length of time because of the long RBC train formed in 27%-H_ct condition. Similar hematocrit elevation at the downstream region with 45%-H_ct in the same 10 μm microvessel is also observed with the cells randomly located. In 20 μm microvessels with 45%-H_ct, the Pf-IRBCs slow down the velocity of the healthy red blood cells (HRBCs) around them and then locally elevate the volume fraction and result in the accumulation of the RBCs at the center of the vessels, thus leaving a thicker cell free layer (CFL) near the vessel wall than normal. Variation of wall shear stress (WSS) is caused by the fluctuation of local H_ct and the distance between the wall and the RBCs. Moreover, in high-hematocrit condition (45%), malaria-infected cells have a tendency to migrate to the edge of the aggregates which is due to the uninterrupted hydrodynamic interaction between the HRBCs and Pf-IRBC. Our results suggest that the existence of Pf-IRBCs is a nonnegligible factor for the fluctuation of hematocrit and WSS and also contributes to the increase of CFL of pathological blood flow in microvessels. The numerical approach presented has the potential to be utilized to RBC disorders and other hematologic diseases.

Elastohydrodynamic Lift at a Soft Wall

We study experimentally the motion of nondeformable microbeads in a linear shear flow close to a wall bearing a thin and soft polymer layer. Combining microfluidics and 3D optical tracking, we demonstrate that the steady-state bead-to-surface distance increases with the flow strength. Moreover, such lift is shown to result from flow-induced deformations of the layer, in quantitative agreement with theoretical predictions from elastohydrodynamics. This study thus provides the first experimental evidence of "soft lubrication" at play at small scale, in a system relevant, for example, to the physics of blood microcirculation.

Manipulating nanoparticle transport within blood flow through external forces: an exemplar of mechanics in nanomedicine

A large number of nanoparticles (NPs) have been raised for diverse biomedical applications and some of them have shown great potential in treatment and imaging of diseases. Design of NPs is essential for delivery efficacy due to a number of biophysical barriers, which prevents the circulation of NPs in vascular flow and their accumulation at tumour sites. The physiochemical properties of NPs, so-called '4S' parameters, such as size, shape, stiffness and surface functionalization, play crucial roles in their life journey to be delivered to tumour sites. NPs can be modified in various ways to extend their blood circulation time and avoid their clearance by phagocytosis, and efficiently diffuse into tumour cells. However, it is difficult to overcome these barriers simultaneously by a simple combination of '4S' parameters for NPs. At this moment, external triggerings are necessary to guide the movement of NPs, which include light, ultrasound, magnetic field, electrical field and chemical interaction. The delivery system can be constructed to be sensitive to these external stimuli which can reduce the non-specific toxicity and improve the efficacy of the drug-delivery system. From a mechanics point of view, we discuss how different forces play their roles in the margination of NPs in blood flow and tumour microvasculature.

A parallel fluid-solid coupling model using LAMMPS and Palabos based on the immersed boundary method

The study of viscous fluid flow coupled with rigid or deformable solids has many applications in biological and engineering problems, e.g., blood cell transport, drug delivery, and particulate flow. We developed a partitioned approach to solve this coupled Multiphysics problem. The fluid motion was solved by Palabos (Parallel Lattice Boltzmann Solver), while the solid displacement and deformation was simulated by LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). The coupling was achieved through the immersed boundary method (IBM). The code modeled both rigid and deformable solids exposed to flow. The code was validated with the Jeffery orbits of an ellipsoid particle in shear flow, red blood cell stretching test, and effective blood viscosity flowing in tubes. It demonstrated essentially linear scaling from 512 to 8192 cores for both strong and weak scaling cases. The computing time for the coupling increased with the solid fraction. An example of the fluid-solid coupling was given for flexible filaments (drug carriers) transport in a flowing blood cell suspensions, highlighting the advantages and capabilities of the developed code.

Periodic and simultaneous quantification of blood viscosity and red blood cell aggregation using a microfluidic platform under in-vitro closed-loop circulation

To evaluate variations of blood circulating in closed loops, hemorheological properties including blood viscosity and red blood cells (RBCs) are quantitatively measured with independent in-vitro instruments after collecting blood from a closed loop. But, most previous methods require periodic blood collections which induce several problems such as geometric differences between the fluidic channel and the in-vitro method, hemodilution, storage time, and unspecific blood flow rates. To resolve these issues, in this study, blood viscosity and RBC aggregation of blood circulating within a closed loop are measured with a microfluidic platform periodically and simultaneously. To demonstrate the proposed method, in-vitro closed-loop circulation is established by connecting several components (peristaltic pump, air compliance unit, fluid divider, and reservoir) in series. In addition, to measure blood viscosity and RBC aggregation, a microfluidic platform composed of a microfluidic device, pinch valve, and syringe pump is created. During each period, blood viscosity and RBC aggregation are measured by monitoring blood flow at constant blood flow, and image intensity at stationary blood flow. The proposed method is first employed to evaluate the effect of hematocrits and dextran concentrations on the RBC aggregation and blood viscosity by using a syringe pump (i.e., specific blood flow-rate). The method is then applied to detect the blood viscosity and RBC aggregation under in-vitro closed-loop circulation (i.e., unspecific blood flow-rate). From these experimental demonstrations, it is found that the suggested method can be effectively used to monitor the RBC aggregation and blood viscosity under in-vitro closed-loop circulation. Since this method does not require periodic collection from closed-loop circulation or an additional procedure for estimating blood flow-rate with a syringe pump, it will be effectively used to monitor variations of blood circulating in extracorporeal bypass loops.

The quaternary state of polymerized human hemoglobin regulates oxygenation of breast cancer solid tumors: A theoretical and experimental study

A major constraint in the treatment of cancer is inadequate oxygenation of the tumor mass, which can reduce chemotherapeutic efficacy. We hypothesize that polymerized human hemoglobin (PolyhHb) can be transfused into the systemic circulation to increase solid tumor oxygenation, and improve chemotherapeutic outcomes. By locking PolyhHb in the relaxed (R) quaternary state, oxygen (O2) offloading at low O2 tensions (<20 mm Hg) may be increased, while O2 offloading at high O2 tensions (>20 mm Hg) is facilitated with tense (T) state PolyhHb. Therefore, R-state PolyhHb may deliver significantly more O2 to hypoxic tissues. Biophysical parameters of T and R-state PolyhHb were used to populate a modified Krogh tissue cylinder model to assess O2 transport in a tumor. In general, we found that increasing the volume of transfused PolyhHb decreased the apparent viscosity of blood in the arteriole. In addition, we found that PolyhHb transfusion decreased the wall shear stress at large arteriole diameters (>20 μm), but increased wall shear stress for small arteriole diameters (<10 μm). Therefore, transfusion of PolyhHb may lead to elevated O2 delivery at low pO2. In addition, transfusion of R-state PolyhHb may be more effective than T-state PolyhHb for O2 delivery at similar transfusion volumes. Reduction in the apparent viscosity resulting from PolyhHb transfusion may result in significant changes in flow distributions throughout the tumor microcirculatory network. The difference in wall shear stress implies that PolyhHb may have a more significant effect in capillary beds through mechano-transduction. Periodic top-load transfusions of PolyhHb into mice bearing breast tumors confirmed the oxygenation potential of both PolyhHbs via reduced hypoxic volume, vascular density, tumor growth, and increased expression of hypoxia inducible genes. Tissue section analysis demonstrated primary PolyhHb clearance occurred in the liver and spleen indicating a minimal risk for renal damage.

Hydrodynamic instabilities in miscible fluids

Hydrodynamic instabilities in miscible fluids are ubiquitous, from natural phenomena up to geological scales, to industrial and technological applications, where they represent the only way to control and promote mixing at low Reynolds numbers, well below the transition from laminar to turbulent flow. As for immiscible fluids, the onset of hydrodynamic instabilities in miscible fluids is directly related to the physics of their interfaces. The focus of this review is therefore on the general mechanisms driving the growth of disturbances at the boundary between miscible fluids, under a variety of forcing conditions. In the absence of a regularizing mechanism, these disturbances would grow indefinitely. For immiscible fluids, interfacial tension provides such a regularizing mechanism, because of the energy cost associated to the creation of new interface by a growing disturbance. For miscible fluids, however, the very existence of interfacial stresses that mimic an effective surface tension is debated. Other mechanisms, however, may also be relevant, such as viscous dissipation. We shall review the stabilizing mechanisms that control the most common hydrodynamic instabilities, highlighting those cases for which the lack of an effective interfacial tension poses deep conceptual problems in the mathematical formulation of a linear stability analysis. Finally, we provide a short overview on the ongoing research on the effective, out of equilibrium interfacial tension between miscible fluids.

Visible light optical coherence microscopy of the brain with isotropic femtoliter resolution in vivo

Most flying-spot optical coherence tomography and optical coherence microscopy (OCM) systems use a symmetric confocal geometry, where the detection path retraces the illumination path starting from and ending with the spatial mode of a single-mode optical fiber. Here we describe a visible light OCM instrument that breaks this symmetry to improve transverse resolution without sacrificing collection efficiency in scattering tissue. This was achieved by overfilling a water immersion objective on the illumination path while maintaining a conventional Gaussian mode detection path (1/e2 intensity diameter ∼0.82 Airy disks), enabling ∼1.1  μm full width at half-maximum (FWHM) transverse resolution. At the same time, a ∼0.9  μm FWHM axial resolution in tissue, achieved by a broadband visible light source, enabled femtoliter volume resolution. We characterized this instrument according to paraxial coherent microscopy theory and, finally, used it to image the meningeal layers, intravascular red blood cell-free layer, and myelinated axons in the mouse neocortex in vivo through the thinned skull.

Image-based modelling of skeletal muscle oxygenation

The supply of oxygen in sufficient quantity is vital for the correct functioning of all organs in the human body, in particular for skeletal muscle during exercise. Disease is often associated with both an inhibition of the microvascular supply capability and is thought to relate to changes in the structure of blood vessel networks. Different methods exist to investigate the influence of the microvascular structure on tissue oxygenation, varying over a range of application areas, i.e. biological in vivo and in vitro experiments, imaging and mathematical modelling. Ideally, all of these methods should be combined within the same framework in order to fully understand the processes involved. This review discusses the mathematical models of skeletal muscle oxygenation currently available that are based upon images taken of the muscle microvasculature in vivo and ex vivo Imaging systems suitable for capturing the blood vessel networks are discussed and respective contrasting methods presented. The review further informs the association between anatomical characteristics in health and disease. With this review we give the reader a tool to understand and establish the workflow of developing an image-based model of skeletal muscle oxygenation. Finally, we give an outlook for improvements needed for measurements and imaging techniques to adequately investigate the microvascular capability for oxygen exchange.

Variation in wall shear stress in channel networks of zebrafish models

Physiological functions of vascular endothelial cells (ECs) vary depending on wall shear stress (WSS) magnitude, and the functional change affects the pathologies of various cardiovascular systems. Several in vitro and in vivo models have been used to investigate the functions of ECs under different WSS conditions. However, these models have technical limitations in precisely mimicking the physiological environments of ECs and monitoring temporal variations of ECs in detail. Although zebrafish (Danio rerio) has several strategies to overcome these technical limitations, zebrafish cannot be used as a perfect animal model because applying various WSS conditions on blood vessels of zebrafish is difficult. This study proposes a new zebrafish model in which various WSS can be applied to the caudal vein. The WSS magnitude is controlled by blocking some parts of blood-vessel networks. The accuracy and reproducibility of the proposed method are validated using an equivalent circuit model of blood vessels in zebrafish. The proposed method is applied to lipopolysaccharide (LPS)-stimulated zebrafish as a typical application. The proposed zebrafish model can be used as an in vivo animal model to investigate the relationship between WSS and EC physiology or WSS-induced cardiovascular diseases.

Margination of Fluorescent Polylactic Acid-Polyaspartamide based Nanoparticles in Microcapillaries In Vitro: the Effect of Hematocrit and Pressure

The last decade has seen the emergence of vascular-targeted drug delivery systems as a promising approach for the treatment of many diseases, such as cardiovascular diseases and cancer. In this field, one of the major challenges is carrier margination propensity (i.e., particle migration from blood flow to vessel walls); indeed, binding of these particles to targeted cells and tissues is only possible if there is direct carrier–wall interaction. Here, a microfluidic system mimicking the hydrodynamic conditions of human microcirculation in vitro is used to investigate the effect of red blood cells (RBCs) on a carrier margination in relation to RBC concentration (hematocrit) and pressure drop. As model drug carriers, fluorescent polymeric nanoparticles (FNPs) were chosen, which were obtained by using as starting material a pegylated polylactic acid–polyaspartamide copolymer. The latter was synthesized by derivatization of α,β-poly(N-2-hydroxyethyl)-d,l-aspartamide (PHEA) with Rhodamine (RhB), polylactic acid (PLA) and then poly(ethyleneglycol) (PEG) chains. It was found that the carrier concentration near the wall increases with increasing pressure drop, independently of RBC concentration, and that the tendency for FNP margination decreases with increasing hematocrit. This work highlights the importance of taking into account RBC–drug carrier interactions and physiological conditions in microcirculation when planning a drug delivery strategy based on systemically administered carriers.

Margination and stretching of von Willebrand factor in the blood stream enable adhesion

The protein von Willebrand factor (VWF) is essential in primary hemostasis, as it mediates platelet adhesion to vessel walls. VWF retains its compact (globule-like) shape in equilibrium due to internal molecular associations, but is able to stretch when a high enough shear stress is applied. Even though the shear-flow sensitivity of VWF conformation is well accepted, the behavior of VWF under realistic blood flow conditions remains poorly understood. We perform mesoscopic numerical simulations together with microfluidic experiments in order to characterize VWF behavior in blood flow for a wide range of flow-rate and hematocrit conditions. In particular, our results demonstrate that the compact shape of VWF is important for its migration (or margination) toward vessel walls and that VWF stretches primarily in a near-wall region in blood flow making its adhesion possible. Our results show that VWF is a highly optimized protein in terms of its size and internal associations which are necessary to achieve its vital function. A better understanding of the relevant mechanisms for VWF behavior in microcirculation provides a further step toward the elucidation of the role of mutations in various VWF-related diseases.

Cellular Level In-silico Modeling of Blood Rheology with An Improved Material Model for Red Blood Cells

Many of the intriguing properties of blood originate from its cellular nature. Therefore, accurate modeling of blood flow related phenomena requires a description of the dynamics at the level of individual cells. This, however, presents several computational challenges that can only be addressed by high performance computing. We present Hemocell, a parallel computing framework which implements validated mechanical models for red blood cells and is capable of reproducing the emergent transport characteristics of such a complex cellular system. It is computationally capable of handling large domain sizes, thus it is able to bridge the cell-based micro-scale and macroscopic domains. We introduce a new material model for resolving the mechanical responses of red blood cell membranes under various flow conditions and compare it with a well established model. Our new constitutive model has similar accuracy under relaxed flow conditions, however, it performs better for shear rates over 1,500 s^-1. We also introduce a new method to generate randomized initial conditions for dense mixtures of different cell types free of initial positioning artifacts.

An Approach to Identifying Phenomena Accompanying Micro and Nanoparticles in Contact With Irregular Vessel Walls

The objective of this paper is to present the method for determining the nature and values of the forces needed to set micro and nanoparticles sitting immobile at the blood vessel wall in motion. The problem was tackled in two ways. Microparticles were examined as objects coming into contact with the wall with the actual large arteriole-type vessel structure. The forces acting on microparticles 10, 30, and [Formula: see text] in diameter were determined: drag force F_D , lift force F_L , electrostatic force F_E , and gravity force F_G . Fluid-structure interaction analysis was used to research the problem. However, nanoparticles were examined as objects coming into contact with the endothelial surface layer (ESL). Resistance forces during the movement of nanoparticles 20, 50, and 100 nm in diameter in the ESL were determined. The same was done for aggregates of nanoparticles 50 nm in diameter. Local irregularities in wall surface are important for microparticles. Small irregularities with the small values of electrostatic force F_E can effectively stop the particle. In the case of nanoparticles, the key is the interaction of the particle with ESL. The research methodology presented can be used to better understand the particle-blood vessel wall interaction phenomena, leading to a more informed particle movement control. The new application of known calculation methods presented in this paper can be successfully used as an additional tool that simplifies planning and design of strategies for drug delivery by means of micro and nanoparticles.

Drag reducing polymers decrease hepatic injury and metastases after liver ischemia-reperfusion

Introduction

Surgery, a crucial therapeutic modality in the treatment of solid tumors, can induce sterile inflammatory processes which can result in metastatic progression. Liver ischemia and reperfusion (I/R) injury, an inevitable consequence of hepatic resection of metastases, has been shown to foster hepatic capture of circulating cancer cells and accelerate metastatic growth. Efforts to reduce these negative consequences have not been thoroughly investigated. Drag reducing polymers (DRPs) are blood-soluble macromolecules that can, in nanomolar concentrations, increase tissue perfusion, decrease vascular resistance and decrease near-wall microvascular concentration of neutrophils and platelets thereby possibly reducing the inflammatory microenvironment. We hypothesize that DRP can potentially be used to ameliorate metastatic capture of tumor cells and tumor growth within the I/R liver.

Methods

Experiments were performed utilizing a segmental ischemia model of mice livers. Five days prior or immediately prior to ischemia, murine colon adenocarcinoma cells (MC38) were injected into the spleen. DRP (polyethylene oxide) or a control of low-molecular-weight polyethylene glycol without drag reducing properties were administered intraperitoneally at the onset of reperfusion.

Results

After three weeks from I/R, we observed that liver I/R resulted in an increased ability to capture and foster growth of circulating tumor cells; in addition, the growth of pre-existing micrometastases was accelerated three weeks later. These effects were significantly curtailed when mice were treated with DRPs at the time of I/R. Mechanistic investigations in vivo indicated that DRPs protected the livers from I/R injury as evidenced by significant decreases in hepatocellular damage, neutrophil recruitment into the liver, formation of neutrophil extracellular traps, deposition of platelets, formation of microthrombi within the liver sinusoids and release of inflammatory cytokines.

Conclusions

DRPs significantly attenuated metastatic tumor development and growth. DRPs warrant further investigation as a potential treatment for liver I/R injury in the clinical setting to improve cancer-specific outcomes.

Theory and algorithms to compute Helfrich bending forces: a review

Cell membranes are vital to shield a cell's interior from the environment. At the same time they determine to a large extent the cell's mechanical resistance to external forces. In recent years there has been considerable interest in the accurate computational modeling of such membranes, driven mainly by the amazing variety of shapes that red blood cells and model systems such as vesicles can assume in external flows. Given that the typical height of a membrane is only a few nanometers while the surface of the cell extends over many micrometers, physical modeling approaches mostly consider the interface as a two-dimensional elastic continuum. Here we review recent modeling efforts focusing on one of the computationally most intricate components, namely the membrane's bending resistance. We start with a short background on the most widely used bending model due to Helfrich. While the Helfrich bending energy by itself is an extremely simple model equation, the computation of the resulting forces is far from trivial. At the heart of these difficulties lies the fact that the forces involve second order derivatives of the local surface curvature which by itself is the second derivative of the membrane geometry. We systematically derive and compare the different routes to obtain bending forces from the Helfrich energy, namely the variational approach and the thin-shell theory. While both routes lead to mathematically identical expressions, so-called linear bending models are shown to reproduce only the leading order term while higher orders differ. The main part of the review contains a description of various computational strategies which we classify into three categories: the force, the strong and the weak formulation. We finally give some examples for the application of these strategies in actual simulations.

Artificial Red Blood Cells as Potential Photosensitizers in Dye Laser Treatment Against Port-Wine Stains

We suggest a novel method that uses artificial blood cells (hemoglobin vesicles, Hb-Vs) as photosensitizers in dye laser treatment (at 595-nm wavelength) for port-wine stains (i.e., capillary malformations presenting as red birthmarks) based on the results of animal experiments. As compared with human red blood cells, Hb-Vs have the same absorbance of 595 nm wavelength light and produce the same level of heat following dye laser irradiation. Small sized Hb-Vs (250 nm) distribute in the plasma phase in blood and tend to flow in the marginal zone of microvessels. Intravenous injections of Hb-Vs caused the dilatation of microvessels, and dye laser treatment with Hb-Vs destroyed the vessel wall effectively. Following the intravenous injection of Hb-Vs, the microvessels contained more Hb that absorbed laser photons and produced heat. This extra Hb tended to flow near the endothelial cells, which were the target of the laser treatment. These attributes of Hb-Vs will potentially contribute to enhancing the efficacy of dye laser treatment for port-wine stains. Hemoglobin is a type of porphyrin. Thus, our proposed treatment may have aspects of photodynamic therapy using porphyrin that leads to a cytotoxicity effect by active oxygen.

Microvasculature on a chip: study of the Endothelial Surface Layer and the flow structure of Red Blood Cells

Microvasculatures-on-a-chip, i.e. in vitro models that mimic important features of microvessel networks, have gained increasing interest in recent years. Such devices have allowed investigating pathophysiological situations involving abnormal biophysical interactions between blood cells and vessel walls. Still, a central question remains regarding the presence, in such biomimetic systems, of the endothelial glycocalyx. The latter is a glycosaminoglycans-rich surface layer exposed to blood flow, which plays a crucial role in regulating the interactions between circulating cells and the endothelium. Here, we use confocal microscopy to characterize the layer expressed by endothelial cells cultured in microfluidic channels. We show that, under our culture conditions, endothelial cells form a confluent layer on all the walls of the circuit and display a glycocalyx that fully lines the lumen of the microchannels. Moreover, the thickness of this surface layer is found to be on the order of 600 nm, which compares well with measurements performed ex or in vivo on microcapillaries. Furthermore, we investigate how the presence of endothelial cells in the microchannels affects their hydrodynamic resistance and the near-wall motion of red blood cells. Our study thus provides an important insight into the physiological relevance of in vitro microvasculatures.

Blood viscosity monitoring during cardiopulmonary bypass based on pressure-flow characteristics of a Newtonian fluid

We proposed a blood viscosity estimation method based on pressure-flow characteristics of oxygenators used during cardiopulmonary bypass (CPB) in a previous study that showed the estimated viscosity to correlate well with the measured viscosity. However, the determination of the parameters included in the method required the use of blood, thereby leading to high cost of calibration. Therefore, in this study we propose a new method to monitor blood viscosity, which approximates the pressure-flow characteristics of blood considered as a non-Newtonian fluid with characteristics of a Newtonian fluid by using the parameters derived from glycerin solution to enable ease of acquisition. Because parameters used in the estimation method are based on fluid types, bovine blood parameters were used to calculate estimated viscosity (ηe), and glycerin parameters were used to estimate deemed viscosity (ηdeem). Three samples of whole bovine blood with different hematocrit levels (21.8%, 31.0%, and 39.8%) were prepared and perfused into the oxygenator. As the temperature changed from 37 °C to 27 °C, the oxygenator mean inlet pressure and outlet pressure were recorded for flows of 2 L/min and 4 L/min, and the viscosity was estimated. The value of deemed viscosity calculated with the glycerin parameters was lower than estimated viscosity calculated with bovine blood parameters by 20-33% at 21.8% hematocrit, 12-27% at 31.0% hematocrit, and 10-15% at 39.8% hematocrit. Furthermore, deemed viscosity was lower than estimated viscosity by 10-30% at 2 L/min and 30-40% at 4 L/min. Nevertheless, estimated and deemed viscosities varied with a similar slope. Therefore, this shows that deemed viscosity achieved using glycerin parameters may be capable of successfully monitoring relative viscosity changes of blood in a perfusing oxygenator.

Predicting different adhesive regimens of circulating particles at blood capillary walls

A fundamental step in the rational design of vascular targeted particles is the firm adhesion at the blood vessel walls. Here, a combined lattice Boltzmann-immersed boundary model is presented for predicting the near-wall dynamics of circulating particles. A moving least squares algorithm is used to reconstruct the forcing term accounting for the immersed particle, whereas ligand-receptor binding at the particle-wall interface is described via forward and reverse probability distributions. First, it is demonstrated that the model predicts with good accuracy the rolling velocity of tumor cells over an endothelial layer in a microfluidic channel. Then, particle-wall interactions are systematically analyzed in terms of particle geometries (circular, elliptical with aspect ratios 2 and 3), surface ligand densities (0.3, 0.5, 0.7 and 0.9), ligand-receptor bond strengths (1 and 2) and Reynolds numbers (Re = 0.01, 0.1 and 1.0). Depending on these conditions, four different particle-wall interaction regimens are identified, namely not adhering, rolling, sliding and firmly adhering particles. The proposed computational strategy can be efficiently used for predicting the near-wall dynamics of particles with arbitrary geometries and surface properties and represents a fundamental tool in the rational design of particles for the specific delivery of therapeutic and imaging agents.

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

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

Facilitating the Clinical Integration of Nanomedicines: The Roles of Theoretical and Computational Scientists

Since the launch of multiple research initiatives on nanotechnology applied to medicine in the early 2000s, a plethora of nanomedicines have been developed that exhibit great therapeutic efficacy in preclinical models but yet minimal impact in daily clinical practice. The successful and complete clinical fruition of nanomedicines requires addressing three major technical challenges: improving loading efficacy and on-command release, modulating recognition and sequestration by immune cells, and maximizing accumulation at biological targets. In this Perspective, I describe how theoretical and computational models can help address each of these challenges. This armamentarium represents an ideal tool for maximizing the therapeutic efficacy of nanomedicines, thus facilitating their integration into daily clinical operations.

Flow of Red Blood Cells in Stenosed Microvessels

A computational study is presented on the flow of deformable red blood cells in stenosed microvessels. It is observed that the Fahraeus-Lindqvist effect is significantly enhanced due to the presence of a stenosis. The apparent viscosity of blood is observed to increase by several folds when compared to non-stenosed vessels. An asymmetric distribution of the red blood cells, caused by geometric focusing in stenosed vessels, is observed to play a major role in the enhancement. The asymmetry in cell distribution also results in an asymmetry in average velocity and wall shear stress along the length of the stenosis. The discrete motion of the cells causes large time-dependent fluctuations in flow properties. The root-mean-square of flow rate fluctuations could be an order of magnitude higher than that in non-stenosed vessels. Several folds increase in Eulerian velocity fluctuation is also observed in the vicinity of the stenosis. Surprisingly, a transient flow reversal is observed upstream a stenosis but not downstream. The asymmetry and fluctuations in flow quantities and the flow reversal would not occur in absence of the cells. It is concluded that the flow physics and its physiological consequences are significantly different in micro- versus macrovascular stenosis.

Microdevice for plasma separation from whole human blood using bio-physical and geometrical effects

In this research work, we present a simple and efficient passive microfluidic device for plasma separation from pure blood. The microdevice has been fabricated using conventional photolithography technique on a single layer of polydimethylsiloxane, and has been extensively tested on whole blood and enhanced (upto 62%) hematocrit levels of human blood. The microdevice employs elevated dimensions of about 100 μm; such elevated dimensions ensure clog-free operation of the microdevice and is relatively easy to fabricate. We show that our microdevice achieves almost 100% separation efficiency on undiluted blood in the flow rate range of 0.3 to 0.5 ml/min. Detailed biological characterization of the plasma obtained from the microdevice is carried out by testing: proteins by ultra-violet spectrophotometric method, hCG (human chorionic gonadotropin) hormone, and conducting random blood glucose test. Additionally, flow cytometry study has also been carried on the separated plasma. These tests attest to the high quality of plasma recovered. The microdevice developed in this work is an outcome of extensive experimental research on understanding the flow behavior and separation phenomenon of blood in microchannels. The microdevice is compact, economical and effective, and is particularly suited in continuous flow operations.

Hemodynamics

A review is presented of the physical principles governing the distribution of blood flow and blood pressure in the vascular system. The main factors involved are the pulsatile driving pressure generated by the heart, the flow characteristics of blood, and the geometric structure and mechanical properties of the vessels. The relationship between driving pressure and flow in a given vessel can be understood by considering the viscous and inertial forces acting on the blood. Depending on the vessel diameter and other physical parameters, a wide variety of flow phenomena can occur. In large arteries, the propagation of the pressure pulse depends on the elastic properties of the artery walls. In the microcirculation, the fact that blood is a suspension of cells strongly influences its flow properties and leads to a nonuniform distribution of hematocrit among microvessels. The forces acting on vessel walls include shear stress resulting from blood flow and circumferential stress resulting from blood pressure. Biological responses to these forces are important in the control of blood flow and the structural remodeling of vessels, and also play a role in major disease processes including hypertension and atherosclerosis. Consideration of hemodynamics is essential for a comprehensive understanding of the functioning of the circulatory system.

Modeling microcirculatory blood flow: current state and future perspectives

Microvascular blood flow determines a number of important physiological processes of an organism in health and disease. Therefore, a detailed understanding of microvascular blood flow would significantly advance biophysical and biomedical research and its applications. Current developments in modeling of microcirculatory blood flow already allow to go beyond available experimental measurements and have a large potential to elucidate blood flow behavior in normal and diseased microvascular networks. There exist detailed models of blood flow on a single cell level as well as simplified models of the flow through microcirculatory networks, which are reviewed and discussed here. The combination of these models provides promising prospects for better understanding of blood flow behavior and transport properties locally as well as globally within large microvascular networks.

Association of Early Atherosclerosis with Vascular Wall Shear Stress in Hypercholesterolemic Zebrafish

Although atherosclerosis is a multifactorial disease, the role of hemodynamic information has become more important. Low and oscillating wall shear stress (WSS) that changes its direction is associated with the early stage of atherosclerosis. Several in vitro and in vivo models were proposed to reveal the relation between the WSS and the early atherosclerosis. However, these models possess technical limitations in mimicking real physiological conditions and monitoring the developmental course of the early atherosclerosis. In this study, a hypercholesterolaemic zebrafish model is proposed as a novel experimental model to resolve these limitations. Zebrafish larvae are optically transparent, which enables temporal observation of pathological variations under in vivo condition. WSS in blood vessels of 15 days post-fertilisation zebrafish was measured using a micro particle image velocimetry (PIV) technique, and spatial distribution of lipid deposition inside the model was quantitatively investigated after feeding high cholesterol diet for 10 days. Lipids were mainly deposited in blood vessel of low WSS. The oscillating WSS was not induced by the blood flows in zebrafish models. The present hypercholesterolaemic zebrafish would be used as a potentially useful model for in vivo study about the effects of low WSS in the early atherosclerosis.

Hematocrit distribution and tissue oxygenation in large microcirculatory networks

OBJECTIVE

Oxygen tension in the brain is controlled by the microcirculatory supply of RBC, but the effect of non-Newtonian blood flow rheology on tissue oxygenation is not well characterized. This study assesses different biphasic blood flow models for predicting tissue oxygen tension as a function of microcirculatory hemodynamics.

METHODS

Two existing plasma-skimming laws are compared against measured RBC distributions in rat and hamster microcirculatory networks. A novel biphasic blood flow model is introduced. The computational models predict tissue oxygenation in the mesentery, cremaster muscle, and the human secondary cortex.

RESULTS

This investigation shows deficiencies in prior models, including inconsistent plasma-skimming trends and insufficient oxygen perfusion due to the high prevalence (33%) of RBC-free microvessels. Our novel method yields physiologically sound RBC distributions and tissue oxygen tensions within one standard deviation of experimental measurements.

CONCLUSIONS

A simple, novel biphasic blood flow model is introduced with equal or better predictive power when applied to historic raw data sets. It can overcome limitations of prior models pertaining to trifurcations, anastomoses, and loops. This new plasma-skimming law eases the computations of bulk blood flow and hematocrit fields in large microcirculatory networks and converges faster than prior procedures.

Red blood cells affect the margination of microparticles in synthetic microcapillaries and intravital microcirculation as a function of their size and shape

A key step in particle-based drug delivery throughmicrocirculation is particlemigration from blood flow to vesselwalls, also known as “margination”,which promotes particle contact and adhesion to the vesselwall. Margination and adhesion should be independently addressed as two distinct phenomena, considering that the former is a fundamental prerequisite to achieve particle adhesion and subsequent extravasation. Although margination has beenmodeled by numerical simulations and investigated inmodel systems in vitro, experimental studies including red blood cells (RBCs) are lacking. Here, we evaluate the effect of RBCs on margination through microfluidic studies in vitro and by intravital microscopy in vivo.We showthatmargination,which is almost absent when particles are suspended in a cell-free medium, is drastically enhanced by RBCs. This effect is size- and shape-dependent, larger spherical/discoid particles being more effectively marginated both in vitro and in vivo. Our findings can be explained by the collision of particles with RBCs that induces the drifting of the particles towards the vessel walls where they become trapped in the cell-free layer. These results are relevant for the design of drug delivery strategies based on systemically administered carriers.

Inflow/Outflow Boundary Conditions for Particle-Based Blood Flow Simulations: Application to Arterial Bifurcations and Trees

When blood flows through a bifurcation, red blood cells (RBCs) travel into side branches at different hematocrit levels, and it is even possible that all RBCs enter into one branch only, leading to a complete separation of plasma and RBCs. To quantify this phenomenon via particle-based mesoscopic simulations, we developed a general framework for open boundary conditions in multiphase flows that is effective even for high hematocrit levels. The inflow at the inlet is duplicated from a fully developed flow generated in a pilot simulation with periodic boundary conditions. The outflow is controlled by adaptive forces to maintain the flow rate and velocity gradient at fixed values, while the particles leaving the arteriole at the outlet are removed from the system. Upon validation of this approach, we performed systematic 3D simulations to study plasma skimming in arterioles of diameters 20 to 32 microns. For a flow rate ratio 6:1 at the branches, we observed the “all-or-nothing” phenomenon with plasma only entering the low flow rate branch. We then simulated blood-plasma separation in arteriolar bifurcations with different bifurcation angles and same diameter of the daughter branches. Our simulations predict a significant increase in RBC flux through the main daughter branch as the bifurcation angle is increased. Finally, we demonstrated the effectiveness of the new methodology in simulations of blood flow in vessels with multiple inlets and outlets, constructed using an angiogenesis model.

Blood tests, which provide a wealth of information on the state of human health, are often performed on cell-free samples. Therefore, blood-plasma separation needs to be achieved. A simple but effective solution for isolating plasma from blood utilizes capillary bifurcations. In a particle-based simulation study of plasma skimming in capillary bifurcations, the blood flow properties such as velocity and pressure fields differ drastically at the inlet and outlet regions. Therefore, a new open (non-periodic) boundary is required. In this paper, we have developed and validated a general parallel framework for open boundary conditions. This is a non-trivial enabling technology that could be used in all open boundary systems and all particle-based Lagrangian simulations. We performed systematic 3D simulations of blood flow in arteriolar bifurcations and elucidated the biophysical mechanism of blood-plasma separation as well as quantified the effects of branch size and bifurcation angle on cell separation efficiency, which have not been addressed before. We also demonstrated the applicability of the methodology in arterial trees with multiple inlets and outlets.

Effect of variation in hemorheology between human and animal blood on the binding efficacy of vascular-targeted carriers

Animal models are extensively used to evaluate the in vivo functionality of novel drug delivery systems (DDS). However, many variations likely exist in vivo between the animals and human physiological environment that significantly alter results obtained with animal models relative to human system. To date, it is not clear if the variation in hemorheology and hemodynamics between common animal and human models affect the functionality of DDS. This study investigates the role of hemorheology of humans and various animal models in dictating the binding efficiency of model vascular-targeted carriers (VTCs) to the wall in physiological blood flows. Specifically, the adhesion of sLe^A-coated nano- and micro-spheres to inflamed endothelial cells monolayers were conducted via a parallel plate flow chamber assay with steady and disturbed red blood cells (RBCs)-in-buffer and whole blood flows of common animal models. Our results suggest that the ratio of carrier size to RBC size dictate particle binding in blood flow. Additionally, the presence of white blood cells affects the trend of particle adhesion depending on the animal species. Overall, this work sheds light on some deviation in VTC vascular wall interaction results obtained with in vivo animal experimentation from expected outcome and efficiency in vivo in human.