Understanding particle margination in blood flow – A step toward optimized drug delivery systems
Introduction
One of the major causes of human death in the world is cancer; for instance, about 8 million people died from cancer in 2012 [1]. Although significant progress has been achieved in developing treatments for cancer, the impact on patient survival remains rather moderate [2]. Therefore, the development of early detection and therapy strategies for cancer is under active research. One of the ideas with a great potential in this area is the targeted delivery of drugs and imaging agents through micro- and nano-particles [2]. A very challenging task is to design suitable carriers which would meet several demands including good adhesion properties to targeted sites in blood flow, efficient transport through the biological barriers (e.g., vessel walls, interstitial space, and cell membranes), and low clearance by various defense mechanisms of the body [3], [4], [5], [6], [7]. Proposed solutions are polymer conjugates, which are already in clinical use [8], [9], fabricated nano-particles [9], [10], and self-assembled structures from lipids or block copolymers forming liposomes, polymersomes, or worm-like micelles [9]. All these micro- and nano-carriers differ in shape, size, and deformability. Furthermore, the circulatory system of humans consists of blood vessels which vary from several centimeters to a few micrometers with a wide range of flow rates. Hence, to identify advantages and disadvantages of various particles for vascular drug delivery, their behavior in blood flow for different flow rates, hematocrits (i.e., volume fraction of red blood cells), vessel diameters, particle sizes, shapes, and deformability needs to be better understood.
Blood is a complex fluid which consists of red blood cells (RBCs), white blood cells (WBCs), platelets, and the blood plasma. The physiological systemic hematocrit is in the range 37%–54%, varying for males and females [11]. The local hematocrit in microcirculation can be even lower and depends on the vessel diameter [12], [13]. Experimental [14], [15], [16], [17], [18], [19], theoretical [20], [21], [22], and numerical [23], [24], [25], [26] studies have shown that vesicles and RBCs experience a lift force acting away from a wall due to the hydrodynamic interactions with the wall and the RBC shape and deformability. Consequently, RBCs in blood flow migrate toward the center of a vessel and close to the wall a region depleted of RBCs develops, called the RBC-free layer (RBC-FL). In contrast, WBCs and platelets migrate toward vessel walls, which is referred to as margination [27], [28], [29], [30]. The margination of WBCs and platelets is mediated by RBCs and appears to be important for these cells, since in order to perform their biological function they must have a possibility to adhere to the wall. Thus, margination can be thought of as a necessary precondition for the wall adhesion.
In similarity with WBCs and platelets, the margination process in blood flow is also expected to occur for micro- and nano-particles. Experiments [31] on margination of rigid platelet-like particles in channels with a characteristic size within 50 µm–200 µm have shown that a significant number of the particles gets accumulated near a wall if the hematocrit is larger than about 7% and the accumulation is enhanced with increasing hematocrit. Recent microfluidic experiments [32] have shown that rigid spheres with the size of 2 µm display a significantly higher adhesion density than particles with a size of 200 nm–500 nm. Thus, available experiments already show that there exists a strong dependence of particle margination on hematocrit, flow rate, and carrier size. Numerical simulations of blood flow are ideally suited to explore the margination trends for a wide range of conditions and to better understand the corresponding mechanisms.
Recent numerical simulations have considered segregation of a binary mixture of spherical elastic capsules in shear and channel flow [33], [34], [35]. Simulations for different capsule sizes have shown that in binary mixtures of particles, the large particles migrate toward the center and the small ones toward the wall if the volume fraction of the small ones is low. Furthermore, floppy particles tend to migrate to the center in a binary mixture together with rigid particles. Simulations of WBCs in two dimensions (2D) [36], [37] and three dimensions (3D) [38] have shown that WBCs are subject to strong margination for Ht ≈ 0.2–0.4 and the venular range of flow rates, characterized by the pseudo-shear rates in small vessels [39], [40]. Platelet margination has been studied in 2D for two shear rates, two hematocrit values ( and ), and two sizes (0.75 µm and 1.5 µm) of rigid particles [41] in a 50 µm wide channel. The results indicate that platelet margination increases with hematocrit, shear rate, and particle size. A similar trend is observed in another 2D and 3D simulation study of particles between 0.25 µm and 1 µm in a 20 µm wide channel [42], demonstrating that the margination is significantly worse for sub-micrometer particles in comparison to larger carriers in agreement with experimental findings [32]. This numerical investigation has also shown that an oblate discoidal shape has advantages for drug delivery systems compared to spheres. Other 3D simulations [43] have focused on the dependence of platelet margination on hematocrit, their aspect ratio, and different viscosities between the inner and outer fluids of RBCs. The simulation results have shown that margination of rigid particles increases with increasing hematocrit and that the RBC viscosity contrast influences margination as well. Additionally, the study indicated that discs marginate slower than spherical particles. Furthermore, platelet margination in blood flow has been studied numerically for different shear rates and hematocrit values of and [44], [45] for rigid discoidal and spherical particles in a channel of 34 µm. In these investigations, an increase of margination with hematocrit has been observed and related to the different RBC-FL thickness.
In this paper, we systematically investigate the dependence of margination on particle size for a wide range of hematocrit values and flow rates using a 3D model system. In this context, we concentrate on the margination into a thin layer near the wall, where particle adhesion to the wall would be possible. Thus, we assess the adhesion potential of particles in blood flow depending on various carrier and flow properties. In addition, we employ 2D simulations to study the effect of particle deformability on its margination properties. Finally, we consider blood flow and particle margination in channels of several widths. Our results from 3D simulations further support the observations that particles with a few micrometers in size marginate significantly better than their sub-micrometer counterparts. Deformable carriers are in general worse than rigid particles; however, they may marginate slightly better at high Ht values and low shear rates, where their deformability aids them to fit better within a relatively narrow RBC-FL at such conditions. Furthermore, the margination of particles is found to be most efficient in small channels with the sizes corresponding to capillaries in the human microvasculature. Finally, we will also discuss the physical mechanisms of margination and the relation between particle margination, its physical properties, and the thickness of the RBC-FL.
Section snippets
Simulation methods
Two particle-based hydrodynamic simulation approaches, the dissipative particle dynamics (DPD) [46], [47] method and the smoothed dissipative particle dynamics method with angular momentum conservation (SDPD+a) [48], have been employed. DPD has mainly been used for 2D simulations and SDPD+a for 3D simulations. An advantage of SDPD is the possibility to directly control the fluid compressibility and viscosity in comparison to DPD, where these quantities need to be pre-calculated in a separate
Results
The sampling of cross-sectional center-of-mass (COM) positions of the carriers in blood flow over time leads to a distribution, which reflects the probability of a carrier to be at a certain radial position r. Fig. 2(b) displays a COM distribution for RBCs and carriers with the diameter at and . In contrast to RBCs, the carriers migrate into the RBC-FL, a region depleted of RBCs next to the wall, and remain quasi-trapped there. To determine the RBC-FL thickness, the outer
Discussion
Our simulation results clearly support the idea that hematocrit, shear rate, particle size, and channel size strongly affect the efficiency of particle margination. The presence of RBCs is essential in this process, since almost all of the observed effects can be related to the change of the RBC-FL thickness. Without RBCs, which populate the channel center and result in the RBC-FL, suspended particles (e.g., platelets and drug carriers) would occupy the whole channel width. However, even
Conclusion
In conclusion, we have shown that the particle margination is mainly determined by the RBC-FL thickness in relation to the particle size. Accordingly, margination increases with increasing hematocrit, decreasing shear rate, increasing particle size, and decreasing channel width, as long as the particles are not too small in comparison to the RBC-FL. The margination results indicate that a diameter of about 2 µm–3 µm would be advantageous for drug carriers. Deformable particles seem to be
Conflict of interest
The authors declare no conflict of interest.
Ethic approval
No ethical approval is required for this work.
Acknowledgments
This work has been supported by the DFG Research Unit FOR 1543 SHENC – Shear Flow Regulation in Hemostasis. Dmitry A. Fedosov acknowledges funding by the Alexander von Humboldt Foundation. Kathrin Müller acknowledges support by the International Helmholtz Research School of Biophysics and Soft Matter (IHRS BioSoft). We also gratefully acknowledge a CPU time grant by the Jülich Supercomputing Center.
References (82)
- et al.
Control of endothelial targeting and intracellular delivery of therapeutic enzymes by modulating the size and shape of ICAM-1-targeted carriers
Mol Ther
(2008) - et al.
In vivo measurements of “apparent viscosity” and microvessel hematocrit in the mesentery of the cat
Microvasc Res
(1980) - et al.
Microvascular blood flow resistance: role of red blood cell migration and dispersion
Microvasc Res
(2015) - et al.
Maintenance of white blood cell margination at the passage through small venular junctions
Microvasc Res
(1980) - et al.
Margination of leukocytes in blood flow through small tubes
Microvasc Res
(1984) - et al.
The near-wall excess of platelet-sized particles in blood flow: its dependence on hematocrit and wall shear rate
Microvasc Res
(1987) - et al.
Smoothed dissipative particle dynamics with angular momentum conservation
J Comp Phys
(2015) - et al.
Velocity limit in DPD simulations of wall-bounded flows
J Comp Phys
(2008) - et al.
A multiscale red blood cell model with accurate mechanics, rheology, and dynamics
Biophys J
(2010) - et al.
The role of flexible tethers in multiple ligand-receptor bond formation between curved surfaces
Biophys J
(2006)
ADAMTS-13 rapidly cleaves newly secreted ultralarge von Willebrand factor multimers on the endothelial surface under flowing conditions
Blood
Acute hemodynamic effects and blood pool kinetics of polystyrene microspheres following intravenous administration
J Pharm Sci
The effect of shape on the margination dynamics of non-neutrally buoyant particles in two-dimensional shear flows
J Biomech
The margination propensity of ellipsoidal micro/nanoparticles to the endothelium in human blood flow
Biomaterials
The adhesive strength of non-spherical particles mediated by specific interactions
Biomaterials
Platelet adhesion from shear blood flow is controlled by near-wall rebounding collisions with erythrocytes
Biophys J
Finite platelet size could be responsible for platelet margination effect
Biophys J
Nanomedicine in cancer therapy: innovative trends and prospects
Cancer Sci
Cancer nanotechnology: opportunities and challenges
Nat Rev Cancer
Delivering nanomedicine to solid tumors
Nat Rev Clin Oncol
Seven challenges for nanomedicine
Nat Nanotechnol
Delivery of novel therapeutic agents in tumors: physiological barriers and strategies
J Natl Cancer Inst
Polymer conjugates as anticancer nanomedicines
Nat Rev Cancer
Advanced drug delivery systems that target the vascular endothelium
Mol Interv
Direct fabrication and harvesting of monodisperse, shape-specific nanobiomaterials
J Am Chem Soc
Human physiology
Generalization of the Fahraeus principle for microvessel networks
Am J Physiol
Red cell motions and wall interactions in tube flow
Fed Proc
Tank treading and unbinding of deformable vesicles in shear flow: determination of the lift force
Phys Rev Lett
Weakly adhering vesicles in shear flow: tanktreading and anomalous lift force
Europhys Lett
Hydrodynamic lift of vesicles under shear flow in microgravity
Europhys Lett
Separation of blood cells using hydrodynamic lift
Appl Phys Lett
Lift and down-gradient shear-induced diffusion in red blood cell suspensions
Phys Rev Lett
The role of tank-treading motions in the transverse migration of a spheroidal vesicle in a shear flow
J Phys A: Math Gen
Lift force and dynamical unbinding of adhering vesicles under shear flow
Phys Rev Lett
Hydrodynamic lift on bound vesicles
Phys Rev Lett
Influence of shear flow on vesicles near a wall: a numerical study
Phys Rev E
Dynamical regimes and hydrodynamic lift of viscous vesicles under shear
Phys Rev E
The dynamics of a vesicle in a wall-bound shear flow
Phys Fluids
Distribution of blood platelets flowing in arterioles
Am J Physiol
Cited by (71)
Particles in coronary circulation: A review on modelling for drug carrier design
2022, Materials and DesignReal-time visualization of morphology-dependent self-motion of hyaluronic acid nanomaterials in water
2021, International Journal of PharmaceuticsGeneral sites of nanoparticle biodistribution as a novel opportunity for nanomedicine
2021, European Journal of Pharmaceutics and BiopharmaceuticsBlood-nanomaterials interactions
2021, Nanotechnology for Hematology, Blood Transfusion, and Artificial BloodA correction scheme for wall-bounded two-way coupled point-particle simulations
2020, Journal of Computational Physics