**Myron Larrick ^{1}, Nasim Nosoudi^{1}, Surinder P Singh^{2} and Jaime E Ramirez-Vick^{1*}**

^{1}Department of Biomedical, Industrial and Human Factors Engineering, Wright State University, Dayton, OH, USA

^{2}CSIR-National Physical Laboratory, Dr. K S Krishnan Marg, New Delhi, India

**Received:** February 09, 2019

**Accepted:** March 18, 2019

**Version of Record Online: **March 29, 2019

Larrick M, Nosoudi N, Singh SP, Ramirez-Vick JE (2019) Physicochemical Properties of Nanoparticles and Their Effect on Transport Across the Microvasculature. J Ann Bioeng 2019(1): 09-25.

Correspondence should be addressed to

Jaime E Ramirez-Vick

**E-mail:** jaime.ramirez-vick@wright.edu

Copyright © 2019 Jaime E Ramirez-Vick et al. This is an open access article distributed under the Creative Commons Attribution License which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and work is properly cited.

In the past couple of decades nanomedicine has become a well-defined area of interdisciplinary science, in which advances in nanotechnology are allowing the incorporation of multiple therapeutic, sensing and targeting agents into nanostructured materials. This is reflected by the large number of promising nanomedical formulations being currently evaluated by the US Food and Drug Administration. A problem with many, if not all of these, is that mathematical modeling, as a tool to narrow the design parameter space to obtain an optimal formulation, is not part of the design process. Adoption of these tool will bring nanomedicine closer to pharmaceutics, a field in which mathematical modeling has been an essential part of the design process for many years.

Nanoparticle delivery occurs mainly through the microvasculature, comprised of arterioles, capillaries, and venules ranging in diameter from 100 to 300 m, 5 to 10 m, and 7 to 50 m, respectively. The study of how nanoparticle design parameters affect their transport properties as it pertains to targeting and drug delivery has become a growing field in the development of nanomedical formulations. The physicochemical properties of nanoparticles, such as size, shape and surface charge can be optimized to improve their performance in these applications. For instance, it allows modulation of interactions with the immune system, better control over blood clearance and interactions with target tissue, permitting effective delivery of payload within cells or tissues. Their ability to target and enter tissues from the microvasculature is highly dependent on their behavior under blood flow. Here we present a review of the current approaches to mathematical modeling of nanoparticle transport across the vasculature and its effects on targeting efficacy.

Nanomedicine; Nanoparticle Transport; Physicochemical Properties; Vascular Targeting

Cancer has now been noted as the number one cause of death leading those from cardiovascular heart disease for anyone under 85 years old [1]. Nanoparticles (NPs) in the size range of 1 to 100 nm have been effectively used to target tumors like Kapok’s sarcoma, breast cancer, Non-Hodgkin’s lymphoma, ovarian cancer, acute lymphoblastic leukemia, and myeloma cancer [2,3]. Since 1995, there have been 50 nanotherapeutic drugs approved by the US Food and Drug Administration with the global nanomedicine market expected to reach values of $177 billion by 2019, with markets involving *in vitro *diagnostics, *in vivo* imaging, drug therapy, and biomaterials. Liposomal drugs, which are spherical in shape, constitute 80% of this market [4,5], due to their biocompatibility, capacity for self-assembly and ability to carry large drug payloads. Some of these include non-PEGylated and PEGylated liposomes for the drugs DaunoXome, Myocet, Onco TCS, Doxil/Caelyx, Oncaspar, etc., among others [3,6]. It has been shown that NPs that are discoidal tend to marginate toward the vessel wall as compared to the quasi-hemispherical and spherical NPs [7]. It has been seen that the biodistribution of spherical silica particles 700 nm to 3 m in size do not accumulate in non-reticulo-endothelial organs, with increasing NP diameter, when compared to discoidal particles [8]. It has been seen in the pulmonary vasculature of mice that anti-ICAM spheres 0.1, 1, 5 and 10 m in diameter were endocytosed quicker than with delivery using anti-ICAM disks 0.1×3 m^{2} in size, which remained longer in circulation and were more successful targeting the lung versus the liver [9]. Size is a critical parameter in determining NP fate. It has been also shown that delivering the drug methotrexate in 5 nm polyamidoamine-based G5 dendrimers is 10 times more effective targeting tumors than the free drug, with lower toxicity effects [10]. It should also be noted that NPs with diameters smaller than 100 nm have longer circulation time than those larger than 100 nm [11]. NP pharmacokinetics are greatly affected by geometrical factors, for instance high aspect ratios in rods/tubes have been noted to show lower macrophage uptake than spheres due to the large curvature angles, thus allowing longer circulation times [12].

There are many factors that have been investigated that affect the physicochemical properties of NP transport across the microvasculature, which include the shape, size, surface density, surface charge, bulk concentration, length of hydrodynamic forces, blood circulation half-life, and diffusivity [4,7,13]. Nanoparticles with different shapes, including spherical, rod, and disk are discussed in this review. Special consideration must be taken when designing drugs for particular applications. This review presents findings from mathematical modeling completed over the last 10 years for discussions on NPs traversing the microvasculature (Table 1). It outlines findings focusing on results showing shape and size in relation to the NP’s binding affinity. It also compares size and shape to currently approved FDA drugs made in a NP formulation compared to the results from the models.

Model Assumptions |
Governing Equations |
Observations/Results |
References |

Navier-Stokes, laminar shear flow | Small NPs adhesion limited to number of ligand-receptor bonds; larger NPs adhesion related to the force to interaction, Ac. Ellipsoids with aspect ratio of 2 adhere more effectively than spherical NPs | [14] | |

Brownian motion, equation upon total magnetic force on the NPs | β_{vf,100}_{ }is the volume fraction of magnetic NPs required to be in the carrier NP to be captured by the center of the magnet; β_{vf,100}_{ }is approximately 1/d showing that with increased carrier particles fewer magnetic particles need to be contained in the volume fraction; model revealed not NPs can be captured |
[15] | |

Adhesion Model developed for Ligand-receptor pair to follow a Hookean elastic spring | Model showed that for a 50% predicted probability of adhesion with increased NP diameter there is a need for an increased amount of ligands in order to bind to the vessel wall; Design maps were created to show where adhesion occurs and can provide how to design NPs based upon size, surface characteristics of ligand-receptor density and affinity and how to prevent NP endocytosed. | [16] | |

Newton’s Law with Stokes Law applied for NP drag force; magnetic force applied to magnetic NPs | NP capture occurs quicker for increased NP size and increased magnetic flux | [17] | |

laminar flow with Brownian, van der Waals, centrifugal, hydrodynamic, electrostatic forces; NPs follows Newton’s law of motion | Reynolds number less than 0.1 due to size of capillaries. Showed as Stokes number increases it is dependent upon size, density, rotation of NPs where Stokes number of 20 was determined to have maximum margination velocity; disc NPs with smallest aspect ratio of 0.20 marginated toward vessel better than ellipsoidal, hemispherical, and spherical NPs | [18] | |

Navier-Stokes equations in COMSOL Multiphysics with laminar flow; fluid as non-Newtonian Carreau-Model | Larger NPs have a greater chance to adhere to the vessel wall with smaller NPs being swept away; model does not account for pressure changes in the vessel | [19] | |

Non-Newtonian fluid model for magnetic NPs following Herschel-Bulkley fluid | Larger radius of NPs has a greater tendency to be captured by the tumor site; trajectory dependent upon NP radius, volume fraction, and magnetic force location. | [20] | |

Navier-Stokes, laminar shear flow | Mid-sized Disks showed the best adhesion by 2× compared to rods and spheres showed lest ability to adhere. Hydrodynamic forces on the spheres and rods was shown to be greater than on the discs. | [21] | |

2D finite element computational model; Brownian motion; Einstein-Stokes equation for diffusion; apply Fick’s law | Convection to diffusion model developed. Probability of adherence is higher for non-spherical NPs (oblate, rod and disc) compared to spheres where discs displayed the highest adhesion by 300×, rods by 20×, and oblate by 10× | [22] | |

Brownian motion; Naiver-Stokes; incompressible fluid | Model captured effect of ligand-receptor binding. Model used aspect ratio of 3 and 5 where using 5 showed highest binding compared to spheres. Also, smaller NPs had stronger binding. Shear rates going past 2,000 only the rods with aspect ratio of 5 could bind to the vessel wall | [23] | |

Agent based model, Navier-Stokes, laminar flow, fluid having non-Newtonian properties | Model noted polydispersity is important in designing NPs for targeting tumors showing smaller NPs had better uptake. Brownian motion negligible on NPs. Center of vessel has high laminar flow but near the wall lower flow velocity allowing greater NPs to adhere to the wall. Cellular trafficking is incorporated into the model with fenestrations as large as 240 nm for those with tumors | [24] | |

Course-grained model incorporating dissipative particle dynamics | NPs 2, 4, or 6 nm in diameter were not affected by shear rates 0-2000 s^{-1}. Model demonstrated larger NPs that were rod shaped compared to spherical in shape have increased binding abilities. Th model showed the bigger the difference between the ligand-receptor and receptor-solvent the stronger the NP adhesion. |
[25] | |

Navier-Stokes, nonlinear flow, Brownian motion | Flow is modeled as nonlinear partial differentials including magnetic force effects of magnetic NPs; as magnetic field is increased from 1 to 5 kOe the NP velocity decreases from 0.05 to 0.027 cm/s and the fluid velocity decreases from 0.07 to 0.036 cm/s along the vessel wall | [26] |

**Table 1:** Comparison of Computational Models.

Decuzzi, et al., modeled the effect of shape, size and adhesion strength within a laminar shear flow [14]. This model developed the relationship between the optimal volume *V _{opt}* at which the NP had the maximum capacity to be adhere to the surface of receptors as given as:

where the *V _{opt}*

where *a* is the characteristic NP size, *r _{o}* is radius of the NP, the aspect ratio is noted as

For small NPs the adhesion is limited by the number of ligand-receptor bonds, while for large NPs adhesion is related to the force to contact area, *A _{c}*. The model assumes a wall shear stress in the range of 1 to 10 Pa, with surface receptor density of 10

Decuzzi et al., derived the probability of adhesion, *P _{a}(f)*,

where the *β* is the ratio of the number of ligands, *N _{l}*, to receptors,

is a dimensionless number regarding the NP surface binding, is the dimensionless NP radius, is the dimensionless radius of circular area of interaction in the vessel with the NP, is the dimensionless area of adhesion, and

The model showed that for a 50% predicted probability that with increased NP diameter there is a need for an increased number of ligands in order to bind to the surface as hydrodynamic forces increase. The design maps created show three regions where adhesion and endocytosis occur, adhesion and no endocytosis occur, and where no adhesion occurs. These design maps can be used to provide information on design of NPs based upon size, surface characteristics of ligand to receptor density and ligand to receptor binding affinity and how to prevent NPs from being endocytosed. It was revealed that adhesion increases with the surface density of the receptors and ligands.

Gentile et al., developed a mathematical model the flow of NPs here with radius *a =* 200 nm in capillary vessels with treating blood as a Casson fluid where the blood flows with a plug zone of no radial velocity in the center and a parabolic velocity outer zone [27,28]. Immune and red blood cells flowing within these vessels remain within the central plug region, with the remaining annular volume in contact with the vessel wall being known as the cell-free layer. The governing equations for this model are as follows:

represents the fraction of the vessel flow area occupied by the central plug region. Here the model considers effective longitudinal diffusion where *C* is the solute concentration, *u(r)* is the fluid velocity with mean velocity defined as *U*, *Q* is the flow rate, and *D _{m}* is the solute Brownian diffusion coefficient in a quiescent fluid.

Here the model establishes that *D _{eff}* is affected by the molecular diffusion coefficient,

The model showed laminar flow in the capillaries, venules, and arterioles reaching *P _{e}* numbers of 833-41667, (1.66-4.16)×10

Lee et al., formulated a mathematical model considering NPs as spheres, ellipses, discs, and hemispheres within laminar flow conditions with external forces of Hydrodynamic, centrifugal, electrostatic, van der Waals, and Brownian acting on the NPs [18]. Here hydrodynamic forces showed the strongest effect while centrifugal forces were shown to be almost negligible. The governing model considered the flow of particles moving in orientation along the x-axis with the conservation of mass (19) and momentum (20) equations expressed respectively as:

Here the u is the fluid velocity, *ρ _{f}*

The NP was described to adhere to Newton’s law of motion, (*m _{p}* mass of particle and

Here the boundary condition of ** u** = 0 is set at the wall,

This model assumes a low Reynolds number of Re ≤ 0.1, consistent with what is found in capillaries, venules and alveoli vessels. An important parameter for describing NP behavior in blood vessels, is the Stokes number (*St*), which is a dimensionless parameter characterizing the behavior of particles suspended in a fluid flow. As *St* increases its dependent upon size, density and rotation, which in turn results in the NPs having increased lateral drift velocity towards the wall where a maximum margination velocity at *St* ~ 20 was determined in the study. Whereas it was seen that the NPs has a minimal margination occur at *St* equal to 0.1. It was noticed that disc-shaped NPs with the smallest aspect ratio of 0.20 tended to marginate towards the vessel better, when compared to the other shapes, such as ellipsoidal, hemispherical, and spherical. Here the NPs have been shown to increase margination speed, as buoyancy forces are increased through the vessel. Here external hydrodynamic forces need to be applied within the microvasculature for margination to occur as the shear rate is below 100 s^{-1} and non-spherical NPs tend to oscillate as they approach the desired targeted site.

In nanomedicine, the development of magnetic NPs and techniques for their transport and concentration at specific sites in the human body using magnetic forces has shown great interest for many years [29,30]. Furlani et al., presented a mathematical model to predict magnetic trajectory in the targeting of NPs noninvasively with dominant magnetic, *F _{m}*, and fluidic,

where *H** _{a}* is the magnetic field applied to the center of the carrier NP,

and where the *υ _{f}* within the microvasculature is expressed as:

where is the average blood velocity, and *R _{bv}* is vessel radius (i.e., 50 m). While the equation of motion of the carrier NP was given as:

The volume fraction of magnetic particles was derived under the assumption that *-R _{bv}* <

It should be noted that the *β _{vf,100}*

Sharma et al., proposed a mathematical model where a magnetic field is applied to understand the effects of magnetic particles, to reduce reaching healthy cells during delivery [26]. Navier-Stokes equations are established as nonlinear partial differential equations, which include effects from magnetic forces, viscous drag, particle and blood-cell, inertia, buoyancy, gravity, Brownian motion, particle-fluid, and magnetic dipole effects. It demonstrates that NPs have the maximum velocity at the center of the vessel. It was also shown that as the magnetic field is increased from 1 to 5 kOe, the NP velocity decreases from 0.05 to 0.027 cm/s and fluid velocity decreases from 0.07 to 0.036 cm/s along the vessel, respectively.

Babincova et al., proposed a finite element method to solve the model for the trajectory of magnetic NPs with a 50 nm radius in a magnetic field consisting of neodymium rare earth magnets with magnetic energy of 37 MG.Oe for a viscous fluid following Newton’s law [17]:

where *m _{p}* is mass and

where *v _{f }(velocity)* is considered 0 m/s for an ambient quiescent fluid. The model showed the NPs could be captured faster when the particle size was increased. It was also found that by increasing the magnetic flux will decrease the time it takes for the magnetic NPs to be captured.

Heidsieck et al., developed a model for tracking magnetic NPs through blood flow using hydrodynamic forces where the particle trajectories are modeled with external magnetic forces (*F** _{mag}*) consisting of NdFeB N50 Neodymium magnets applied to accelerate the magnetic NPs traveling in the direction of the magnetic force until it adheres to the vessel wall [19]. Buoyancy and inertia forces are present but considered negligible. In addition, Brownian motion is present but considered negligible in this model as well. The model assumed a vessel diameter of 1.2 mm with length 30 mm. The magnetic forces acting on the NPs are expressed as:

where the is the magnetic dipole moment and ** B** is the magnetic flux density field. The fluid flow was modeled using Navier-Stokes equations in COMSOL Multiphysics modeling software with laminar flow occurring at 0.2 m/s. Flow around the NPs was defined using the Stokes drag force on spherical NPs given as:

where *r _{np}* is radius of the NP, γ is the shear rate,

where ** v** is NP velocity and

Shaw et al., developed a non-Newtonian fluid model for magnetic NPs (Fe_{3}O_{4}) in various sizes with the radius of the carrier NP *R _{cp}* = 300 to 1000 nm traveling in blood flow, following the Herschel-Bulkley fluid in micro-vessels with 50 m radius [20]. The applied magnet for this model had a diameter of 6 cm with magnetization

where *m _{cp}* and

Following the Herschel-Bulkley for the fluid the fluid force on a spherical NP and the blood flow laminar the blood velocity can be written as:

where (*1-**ξ _{c}*

Adriani et al., modeled rods and disks using both an experimental flow chamber and computational modeling [21]. The rods were made using an aspect ratio of 4.5 to keep them with a similar loading volume to the disks. The disks ranged in sizes (diameter × height) from 600 × 200 nm, 1000 × 400 nm, and 1800 × 600 nm. The rods ranged in size from 1500 × 200 nm to 1800 × 400 nm. The experimental method conducted utilized mesoporous silicon particles passed through a flow chamber with fluid dynamic trajectories captured using a fluorescence. The results showed that the disks adhered more than the rods by a factor of 2 with significance of shear to 100 s^{-1}. The computational model assumed an incompressible fluid and unsteady flow and applying laminar flow conditions with the Navier-Stokes equation reduced to:

where the ** u** is the fluid velocity,

The hydrodynamic forces (e.g., shear rate) on NPs was observed to be within the range of 10 to 200 s^{-1}. It was found that as the shear rate at the wall increases, the NP adhesion capacity decreases. It was shown that the hydrodynamic forces on the rods was somewhat greater than on the disks. The disks at 1000 × 400 nm showed the most adhesion, when compared to the experimental results. The computational model demonstrated that the hydrodynamic forces affect spheres more when compared to the rods and disks, as can be seen by spheres showing the lowest adhesion, followed by rods and disks, respectively.

Liu et al. presented a 2D finite element computational model incorporating NPs in Brownian motion with flow velocity, *U *(0-25 dyne/cm^{2}), towards adhesion with ligand to receptor binding in the microvasculature 5 m long and 2 m in diameter [22]. Here the convection to diffusion, *D* (10^{-9} m^{2}/s), governing equation for concentration, *c *(*c _{0 }*= 1000 mol/m

where the Boltzmann constant *k _{B}*,

The NPs are modeled following Brownian motion given by velocity (i.e., *V _{s}* for solid and

where, the correction factors are defined as

The model incorporates the fluid flow as an incompressible fluid with the following Navier-Stokes equation where *v _{f}* is the fluid velocity:

The model shows a depletion layer when the drug concentration is reduced as the NP adheres to the receptor surface. Here when the flow rate goes from 0.1 mm/s to 1 mm/s the depletion layer becomes smaller. The model incorporated the probability of NP adhesion, *P _{a}*, given as:

where *m _{r}* and

Tan et al., modeled NPs traversing in Brownian motion through bifurcation vessels with parent diameter of 2 m and daughter diameter of 1 m [23]. The model examined the effects of NP spheres between 100 and 200 nm and rods 63 × 189 nm and 126 × 378 nm, with an aspect ratio of 3, and 52 × 261 nm and 104 × 522 nm, with an aspect ratio of 5. The fluid was treated as an incompressible viscous fluid following Navier-Stokes with the NPs having ligand to receptor, *N _{b}*, binding along the vessel surface, moving with velocity,

The velocity of a NP moving under a deterministic force is given by

where the friction coefficients derived from Stokes’ law are

With the receptor-ligand binding process can be described by

where *k _{f}* is forward binding rate,

The model focused on how the NPs will target a tumor site by having ligands coated on the NPs and they affected the binding on the receptors of the tumor site with a probability of adhesion, *P _{a}*, equivalent to equation (2) from [14] as:

where is the association constant at zero load of the ligand to receptor pair, *λ* is the bond length of the ligand to receptor, *m _{r}* and

It was demonstrated that rod-shaped NPs will bind based on orientation to the vessel wall and inversely proportional to shear stress, with a decrease in binding with shear rates due to increased drag force. Rod-shaped NPs with aspect ratio of 5 showed the highest binding capacity compared to spheres or lower aspect ratio rods. It was shown that smaller NPs (i.e., 100 nm versus 200 nm), tend to show stronger binding capacity due to a lower drag force. In bifurcation areas it was seen that NPs tend to accumulate at higher levels, when compared with straight vessels. It was noticed that spheres cannot bind at shear rate above 1,200 s^{-1} and that shear rates of 2,000 only the NP rods with aspect ratio of 5 were able to bind. Thus, adhesion is affected by flow rate, vessel size, NP geometry, and NP velocity through the fluid.

Hence, the Peclet number (*P _{e}*) seen in Equation (16) was noticed to affect the adhesion ability with length as shown as:

where *L* is the vessel length (0.1 m), *D* is the NP diffusivity coefficient (4.4×10^{-12 }m^{2}/s and 2.2×10^{-12} m^{2}/s) for 100 and 200 nm NPs, respectively, and *U* is the fluid velocity.

Fullstone et al., developed an agent-based mathematical model of NPs being transported in capillaries at 8 m in laminar blood flow to provide method showing how NP size is possible to predict reaching the targeted site [24].

The governing equations established for momentum, mass, and energy using Navier-Stokes for the fluid flow are shown as:

where *ρ* is the density, *δt* is the step time, is the del operator, *v* is the velocity, *p* is the pressure, *T* is the total stress tensor, *F* is the force, *C _{p}* is the specific heat capacity at constant pressure,

Peng et al., developed a coarse-grained mathematical model incorporating Dissipative Particle Dynamics (DPD) to study the effects of binding spherical NPs to vessel walls, while considering shear flow [25]. Here DPD is used to model the hydrodynamic interaction of the vessel and NPs regarding adhesion under shear flow situations where various aspect ratios are compared. Here NPs had ligands attached to the NPs and receptors attached to the vessel surface. The model observed that shear rates that ranged between 0 - 2000 s^{-1} were not affecting the ability of NPs 2, 4, or 6 nm in size to adhere to capillaries. The model compared rods with aspect ratios of 5, 10, and 15 and NPs with diameters of 2, 4, and 6 nm. The model demonstrated that larger spherical NPs and rod-shaped NPs with larger aspect ratios, have increased binding capacity in relation to the time it takes for adhesion to take place, confirming the importance of Brownian motion in this process. The model showed that the bigger the difference between the ligand to receptor and receptor to solvent repulsion, the stronger NP adhesion. It was shown that a 2 nm NPs adhere strongly with a difference of 20 at 1,000 s^{-1}. In DPD simulations, the binding factor Δ*a* is defined as the difference between the ligand-receptor repulsion factor and the receptor-solvent repulsion factor. The model showed that for a weak Δ*a* of 5 to 11, the NPs would remain within the middle of the flow while at larger Δ*a* of 13 to 25, binding occurred faster with a linear shift in improvement going from a 10 to 30% linear increase to almost 100% linear when the difference is at 11, representing the point where NP vessel wall adhesion becomes stable.

Physicochemical properties, such as size, shape, and surface charge, are key design parameters, which optimized to improve NP targeting efficacy. This have been shown by *in vivo* studies in which variation in these parameters influence underlying processes which affect the performance of these nanotherapeutics [32-34]. In this section, we show the effect of size, shape, and surface charge on NP targeting efficacy, based on *in vivo* results.

**Size**

The properties of materials change as their size approaches the nanoscale, this includes melting point, fluorescence, electrical conductivity, magnetic permeability, and chemical reactivity. The variation in these properties as a function of size can be put into practical use in nanomedicine. For instance, at this size scale something called quantum confinement occurs, in which the bandgap varies inversely with size, resulting in electrons and holes being confined into a dimension that approaches a critical quantum measurement, called the exciton Bohr radius [35]. These quantum effects allow the tuning of the NP properties, such a changing the color it fluoresces when excited by changing its size, thus allowing the development of multiple unique labels. In addition, as the size of these NPs decreases, the surface area per volume increases, affecting the amount of material it can interact with, in the physiological environment it is immersed, thus making them more chemically reactive, also affecting their strength and electrical properties [35]. This increase in surface area per unit volume increases the NP contact with the biological environment, affecting its capacity to interact with cells and tissues, circulation time, accumulation and penetration in the tumor, the cell entry mechanism used, and the adsorption of proteins forming the corona [36].

Size significantly influences the blood circulation time and the biodistribution of the NPs, with smaller size having longer circulation times and reduced accumulation in the liver and spleen [37]. In the case of NP accumulation and penetration in the tumor mass, these leak out leak into tumor tissue through permeable tumor vessels and are then retained in the tumor bed due to reduced lymphatic drainage as a result of the enhanced permeability and retention effect [38]. Their accumulation is significantly increased with size, while their penetration decreases with size [39]. When NPs are exposed to a biological environment, proteins start to adsorb on the surface in a competitive manner forming the protein corona. This corona changes over time due to the Vroman effect, where proteins in the environment adsorb on to the NP surface and are replaced over time with other proteins with a higher affinity for the NP surface [40]. It alters the size and interfacial composition of a NP, affecting how these interact with their environment, and specifically cells. This change in size and interfacial composition affects the physiological response including agglomeration, cellular uptake, circulation lifetime, signaling, kinetics, transport, accumulation, and toxicity. The composition of the protein corona is unique to each nanomaterial and depends on many NP parameters (i.e., material, surface, size, charge, and shape) or environmental parameters (i.e., composition, exposure time, pH, temperature, and shear stress) [41,42].

**Shape**

Along with size, NP shape is another key parameter influencing NP behavior in circulation (i.e., fluid dynamics, margination towards the vessel wall) and cell-NP interactions (i.e*.*, binding, internalization, clearance) [42,43]. Although various NP shapes have been synthesized (e.g., spherical, triangular, cubic, rods, platelets, etc.), the NPs under preclinical or clinical studies, are spherical in shape [44]. The aspect ratio is a parameter used to describe the shape of non-spherical NPs, and which has been shown to influence their margination dynamics via lateral drifting toward the blood vessel wall [34]. This hydrodynamic behavior deviating from spherical NPs could be the reason why non-spherical NPs of many shapes have shown extended circulation times, when compared to spherical NPs [45].

Studies involving macrophage uptake or phagocytosis as a function of NP shape suggested that it plays an important role. Specifically, the point of contact between the NP and macrophage, determining if these are capable of forming the necessary actin structures to initiate phagocytosis, or if the NPs will simply spread over the cell membrane [46]. In addition, mammalian epithelial and immune cells preferentially internalize disc-shaped NPs with high aspect ratios compared with nanorods and lower aspect-ratio disc-shaped NPs. While endothelial cells prefer disc-shaped NPs of an intermediate aspect ratio [47]. With epithelial cells internalizing the NPs using the caveolae-mediated pathway, the human umbilical vein endothelial cells, used a clathrin-mediated mechanism.

**Surface charge**

The effective electric charge on the NP surface is measured by the zeta potential, which describes the electrostatic interactions of cells and NPs in a fluid environment, with values of ±30 mV representating stable NPs [48,49]. NP internalization into the cancer cells is promoted with a positive surface, since the surface charge of tumor cells is highly negative compared to healthy cells [39]. The NP surface charge also affects the surface binding of serum proteins during corona formation, their selective adsorption on cellular membranes, and transmembrane permeability. As it was mentioned before, the protein composition of the NP corona is dynamic, changing with time, making the NP surface charge also vary with time [41,42]. Current attempts to control the effect of the protein corona on NP on targeting includes surface modifications to control protein adsorption on to the NP surface (i.e., surface charge, hydrophobicity, and smoothness) [40-42].

FDA approved Nano therapeutics benefits have been shown to reduce drug toxicity, provide longer half-life in circulation by reducing immunogenicity and releasing at projected delivery rates in turn reducing the number of drugs taken. Determining the correct physicochemical properties needed is crucial to the targeting and enhancement of drug efficacy. It is therefore important is design the NPs used for medical drug treatment with optimization to further enhance the efficacy already being witnessed by NPs. A major limitation of current modeling approaches is the lack of consideration of the protein corona. Future mathematical models will greatly improve their capacity to predicts NP targeting and drug delivery efficacy, once the dynamic nature of the protein corona is considered.* *This review focused on the progress development of mathematical models used to predict the ability of a NPs trajectory, adherence in relationship to its physicochemical properties, such as size, shape, and surface charge. Providing valuable design insight regarding the importance of tailoring and optimizing NP physicochemical properties to enhance their margination and wall-adhesion during vascular targeted drug delivery applications.

All authors declare no conflicts of interest in this article.

- Siegel RL, Miller KD, Jemal A (2019) Cancer statistics, 2019. CA Cancer J Clin 69: 7-34.
- Shi J, Kantoff PW, Wooster R, Farokhzad OC (2017) Cancer nanomedicine: progress, challenges and opportunities. Nat Rev Cancer 17: 20-37.
- Sanna V, Pala N, Sechi M (2014) Targeted therapy using nanotechnology: focus on cancer. Int J Nanomedicine 9: 467-483.
- Zhang L, Gu FX, Chan JM, Wang AZ, Langer RS, et al. (2008) Nanoparticles in medicine: therapeutic applications and developments. Clin Pharmacol Ther 83: 761-769.
- Ventola CL (2017) Progress in nanomedicine: approved and investigational nanodrugs. P T 42: 742-755.
- Jones AD, Mi G, Webster TJ (2019) A status report on FDA approval of medical devices containing nanostructured materials. Trends Biotechnol 37: 117-120.
- Decuzzi P, Pasqualini R, Arap W, Ferrari M (2009) Intravascular delivery of particulate systems: does geometry really matter? Pharm Res 26: 235-243.
- Decuzzi P, Godin B, Tanaka T, Lee SY, Chiappini C, et al. (2010) Size and shape effects in the biodistribution of intravascularly injected particles. J Control Release 141: 320-327.
- Muro S, Garnacho C, Champion JA, Leferovich J, Gajewski C, et al. (2008) Control of endothelial targeting and intracellular delivery of therapeutic enzymes by modulating the size and shape of ICAM-1-targeted carriers. Mol Ther 16: 1450-1458.
- Kukowska-Latallo JF, Candido KA, Cao Z, Nigavekar SS, Majoros IJ, et al. (2005) Nanoparticle targeting of anticancer drug improves therapeutic response in animal model of human epithelial cancer. Cancer Res 65: 5317-5324.
- Palombo M, Deshmukh M, Myers D, Gao J, Szekely Z, et al. (2014) Pharmaceutical and toxicological properties of engineered nanomaterials for drug delivery. Annu Rev Pharmacol Toxicol 54: 581-598.
- Champion J, Mitragotri S (2009) Shape induced inhibition of phagocytosis of polymer particles. Pharma Res 26: 244-249.
- Muzykantov VR, Radhakrishnan R, Eckmann DM (2012) Dynamic factors controlling targeting nanocarriers to vascular endothelium. Curr Drug Metab 13: 70-81.
- Decuzzi P, Ferrari M (2006) The adhesive strength of non-spherical particles mediated by specific interactions. Biomaterials 27: 5307-5314.
- Furlani EP, Ng KC (2006) Analytical model of magnetic nanoparticle transport and capture in the microvasculature. Phys Rev E 73: 1-10.
- Decuzzi P, Ferrari M (2008) Design maps for nanoparticles targeting the diseased microvasculature. Biomaterials 29: 377-384.
- Babincova M, Babinec P (2009) Magnetic drug delivery and targeting principles and applications. Biomed Pap Med Fac 153: 243-250.
- Lee SY, Ferrari M, Decuzzi P (2008) Shaping nano-/micro-particles for enhanced vascular interaction in laminar flows. Nanotechnology 20: 495101.
- Heidsieck A, Gleich B (2010) Analysis of particle trajectories for magnetic drug targeting. Proceedings of the 6th COMSOL Conference, Paris, France.
- Shaw S, Murthy PVSN (2010) Magnetic targeting in the impermeable microvessel with two-phase fluid model-non-newtonian characteristics of blood. J Microvasc Res 80: 209-220.
- Adriani G, de Tullio MD, Ferrari M, Hussain F, Pascazio G, et al. (2012) The preferential targeting of the diseased microvasculature by disk-like particles. Biomaterials 33: 5504-5513.
- Liu Y, Shah S, Tan J (2012) Computational modelling of nanoparticle targeted drug delivery. Rev Nanosci Nanotechnol 1: 66-83.
- Tan J, Shah S, Thomas A, Ou-Yang HD, Liu Y (2013) The influence of size, shape and vessel geometry on nanoparticle distribution. Microfluid Nanofluid 14: 77-87.
- Fullstone G, Wood J, Holcombe M, Battaglia G (2015) Modeling the transport of nanoparticles under blood flow using an agent-based approach. Sci Rep 5: 10649.
- Peng B, Liu Y, Zhou Y, Yang L, Zhang G, et al. (2015) Modeling nanoparticle targeting to a vascular surface in shear flow through diffusive particle dynamics. Nanoscale Res Let 10: 235.
- Sharma S, Katiar VK, Singh U (2015) Mathematical modeling of magnetic nano particles in a blood vessel in a magnetic field. J Magnet Mag Mater 379: 102-107.
- Gentile F, Ferrari M, Decuzzi P (2008) The transport of nanoparticles in blood vessels: the effect of vessel permeability and blood rheology. Ann Biomed Eng 36: 254-261.
- Gentile F, Chiappini C, Fine D, Bhavane RC, Peluccio MS, et al. (2008) The effect of shape on the margination dynamics of non-neutrally buoyant particles in two-dimensional shear flows. J. Biomech 41: 2312-2318.
- Kumar A, Jena PK, Behera S, Lockey RF, Mohapatra S, et al. (2010) Multifunctional magnetic nanoparticles for targeted delivery. Nanomedicine 6: 64-69.
- Amirfazli A (2007) Nanomedicine: magnetic nanoparticles hit the target. Nat Nanotechnol 2: 467-468.
- Shibeshi SS, Collins WE (2005) The rheology of blood flow in a branched arterial system. App Rheol 6: 398-405.
- Duan X, Li Y (2013) Physicochemical characteristics of nanoparticles affect circulation, biodistribution, cellular internalization, and trafficking. Small 9: 1521-1532.
- Dufort S, Sancey L, Coll JL (2012) Physico-chemical parameters that govern nanoparticles fate also dictate rules for their molecular evolution. Adv Drug Deliv Rev 64: 179-189.
- Kinnear C, Moore TL, Rodriguez-Lorenzo L, Rothen-Rutishauser B, Petri-Fink A (2017) Form follows function: nanoparticle shape and its implications for nanomedicine. Chem Rev 117: 11476-11521.
- Atluri R, Jensen KA (2017) Engineered nanomaterials: their physicochemical characteristics and how to measure them. Adv Exp Med Biol 947: 3-23.
- Ta HT, Truong NP, Whittaker AK, Davis TP, Peter K (2018) The effects of particle size, shape, density and flow characteristics on particle margination to vascular walls in cardiovascular diseases. Expert Opin Drug Deliv 15: 33-45.
- Hoshyar N, Gray S, Han H, Bao G (2016) The effect of nanoparticle size on
*in vivo*pharmacokinetics and cellular interaction. Nanomedicine 11: 673-692. - Nakamura Y, Mochida A, Choyke PL, Kobayashi H (2016) Nanodrug Delivery: Is the Enhanced Permeability and Retention Effect Sufficient for Curing Cancer? Bioconjug Chem 27: 2225-2238.
- Chakraborty S, Dhakshinamurthy GS, Misra SK (2017) Tailoring of physicochemical properties of nanocarriers for effective anti-cancer applications. J Biomed Mater Res A 105: 2906-2928.
- Mortensen NP, Hurst GB, Wang W, Foster CM, Nallathamby PD, et al. (2013) Dynamic development of the protein corona on silica nanoparticles: composition and role in toxicity. Nanoscale 5: 6372-6380.
- Nguyen VH, Lee BJ (2017) Protein corona: a new approach for nanomedicine design. Int J Nanomedicine 12: 3137-3151.
- Xiao W, Gao H (2018) The impact of protein corona on the behavior and targeting capability of nanoparticle-based delivery system. Int J Pharm 552: 328-339.
- Zhao Z, Ukidve A, Krishnan V, Mitragotri S (2019) Effect of physicochemical and surface properties on
*in vivo*fate of drug nanocarriers. Adv Drug Deliv Rev In Press. - Bobo D, Robinson KJ, Islam J, Thurecht KJ, Corrie SR (2016) Nanoparticle-Based Medicines: A Review of FDA-Approved Materials and Clinical Trials to Date. Pharm Res 33: 2373-87.
- Shi J, Kantoff PW, Wooster R, Farokhzad OC (2017) Cancer nanomedicine: progress, challenges and opportunities. Nat Rev Cancer 17: 20-37.
- Champion JA, Mitragotri S (2006) Role of target geometry in phagocytosis. Proc Natl Acad Sci USA 103: 4930-4934.
- Agarwal R, Singh V, Jurney P, Shi L, Sreenivasan SV, et al. (2013) Mammalian cells preferentially internalize hydrogel nanodiscs over nanorods and use shape-specific uptake mechanisms. Proc Natl Acad Sci USA 110: 17247-17252.
- Honary S, Zahir F (2013) Effect of zeta potential on the properties of nano-drug delivery systems - A review (Part 1). Trop J Pharm Res 12: 255-264.
- Honary S, Zahir F (2013) Effect of zeta potential on the properties of nano-drug delivery systems - A review (Part 2). Trop J Pharm Res 12: 265-273.

This work is licensed under a Creative Commons Attribution 4.0 International License.

Copyright © 2019 - All Rights Reserved - ocimumpublishers.com