Understanding particle margination in blood flow – A step toward optimized drug delivery systems

https://doi.org/10.1016/j.medengphy.2015.08.009Get rights and content

Highlights

  • Particle margination has been studied for a wide range of flow, channel, and particle properties.

  • Efficient margination depends on the relation between particle size and the thickness of red blood cell-free layer.

  • Particle margination decreases as a particle gets more deformable

  • Particle margination is less efficient in larger channels.

Abstract

Targeted delivery of drugs and imaging agents is very promising to develop new strategies for the treatment of various diseases such as cancer. For an efficient targeted adhesion, the particles have to migrate toward the walls in blood flow – a process referred to as margination. Due to a huge diversity of available carriers, a good understanding of their margination properties in blood flow depending on various flow conditions and particle properties is required. We employ a particle-based mesoscopic hydrodynamic simulation approach to investigate the margination of different carriers for a wide range of hematocrits (volume fraction of red blood cells) and flow rates. Our results show that margination strongly depends on the thickness of the available free space close to the wall, the so-called red blood cell-free layer (RBC-FL), in comparison to the carrier size. The carriers with a few micrometers in size are comparable with the RBC-FL thickness and marginate better than their sub-micrometer counterparts. Deformable carriers, in general, show worse margination properties than rigid particles. Particle margination is also found to be most pronounced in small channels with a characteristic size comparable to blood capillaries. Finally, different margination mechanisms are discussed.

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 γ˙¯80s1 in small vessels [39], [40]. Platelet margination has been studied in 2D for two shear rates, two hematocrit values (Ht=0.2 and Ht=0.4), 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 Ht=0.1 and Ht=0.2 [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 Dc=0.25Dr at Ht=0.3 and γ˙*59. 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)

  • DongJ.-F. et al.

    ADAMTS-13 rapidly cleaves newly secreted ultralarge von Willebrand factor multimers on the endothelial surface under flowing conditions

    Blood

    (2002)
  • SlackJ.D. et al.

    Acute hemodynamic effects and blood pool kinetics of polystyrene microspheres following intravenous administration

    J Pharm Sci

    (1981)
  • GentileF. et al.

    The effect of shape on the margination dynamics of non-neutrally buoyant particles in two-dimensional shear flows

    J Biomech

    (2008)
  • ThompsonA.J. et al.

    The margination propensity of ellipsoidal micro/nanoparticles to the endothelium in human blood flow

    Biomaterials

    (2013)
  • DecuzziP. et al.

    The adhesive strength of non-spherical particles mediated by specific interactions

    Biomaterials

    (2006)
  • TokarevA.A. et al.

    Platelet adhesion from shear blood flow is controlled by near-wall rebounding collisions with erythrocytes

    Biophys J

    (2011)
  • TokarevA.A. et al.

    Finite platelet size could be responsible for platelet margination effect

    Biophys J

    (2011)
  • Cancer fact sheet. http://www.who.int/mediacentre/factsheets/fs297/en/ (February...
  • BlancoE. et al.

    Nanomedicine in cancer therapy: innovative trends and prospects

    Cancer Sci

    (2011)
  • FerrariM.

    Cancer nanotechnology: opportunities and challenges

    Nat Rev Cancer

    (2005)
  • JainR.K. et al.

    Delivering nanomedicine to solid tumors

    Nat Rev Clin Oncol

    (2010)
  • SanhaiW.R. et al.

    Seven challenges for nanomedicine

    Nat Nanotechnol

    (2008)
  • JainR.K.

    Delivery of novel therapeutic agents in tumors: physiological barriers and strategies

    J Natl Cancer Inst

    (1989)
  • DuncanR.

    Polymer conjugates as anticancer nanomedicines

    Nat Rev Cancer

    (2006)
  • DingB.-S. et al.

    Advanced drug delivery systems that target the vascular endothelium

    Mol Interv

    (2006)
  • RollandJ.P. et al.

    Direct fabrication and harvesting of monodisperse, shape-specific nanobiomaterials

    J Am Chem Soc

    (2005)
  • PocockG. et al.

    Human physiology

    (2006)
  • PriesA.R. et al.

    Generalization of the Fahraeus principle for microvessel networks

    Am J Physiol

    (1986)
  • GoldsmithH.L.

    Red cell motions and wall interactions in tube flow

    Fed Proc

    (1971)
  • AbkarianM. et al.

    Tank treading and unbinding of deformable vesicles in shear flow: determination of the lift force

    Phys Rev Lett

    (2002)
  • LorzB. et al.

    Weakly adhering vesicles in shear flow: tanktreading and anomalous lift force

    Europhys Lett

    (2000)
  • CallensN. et al.

    Hydrodynamic lift of vesicles under shear flow in microgravity

    Europhys Lett

    (2008)
  • GeislingerT.M. et al.

    Separation of blood cells using hydrodynamic lift

    Appl Phys Lett

    (2012)
  • GrandchampX. et al.

    Lift and down-gradient shear-induced diffusion in red blood cell suspensions

    Phys Rev Lett

    (2013)
  • OllaP.

    The role of tank-treading motions in the transverse migration of a spheroidal vesicle in a shear flow

    J Phys A: Math Gen

    (1997)
  • CantatI. et al.

    Lift force and dynamical unbinding of adhering vesicles under shear flow

    Phys Rev Lett

    (1999)
  • SeifertU.

    Hydrodynamic lift on bound vesicles

    Phys Rev Lett

    (1999)
  • SukumaranS. et al.

    Influence of shear flow on vesicles near a wall: a numerical study

    Phys Rev E

    (2001)
  • MesslingerS. et al.

    Dynamical regimes and hydrodynamic lift of viscous vesicles under shear

    Phys Rev E

    (2009)
  • ZhaoH. et al.

    The dynamics of a vesicle in a wall-bound shear flow

    Phys Fluids

    (2011)
  • TangelderG.J. et al.

    Distribution of blood platelets flowing in arterioles

    Am J Physiol

    (1985)
  • Cited by (71)

    • General sites of nanoparticle biodistribution as a novel opportunity for nanomedicine

      2021, European Journal of Pharmaceutics and Biopharmaceutics
    • Blood-nanomaterials interactions

      2021, Nanotechnology for Hematology, Blood Transfusion, and Artificial Blood
    View all citing articles on Scopus
    View full text