Nanoparticles are regarded as promising carriers for targeted drug delivery and imaging probes. A fundamental understanding of the dynamics of polymeric nanoparticle targeting to receptor-coated vascular surfaces is therefore of great importance to enhance the design of nanoparticles toward improving binding ability. Although the effects of particle size and shear flow on the binding of nanoparticles to a vessel wall have been studied at the particulate level, a computational model to investigate the details of the binding process at the molecular level has not been developed. In this research, dissipative particle dynamics simulations are used to study nanoparticles with diameters of several nanometers binding to receptors on vascular surfaces under shear flow. Interestingly, shear flow velocities ranging from 0 to 2000 s^{−1} had no effect on the attachment process of nanoparticles very close to the capillary wall. Increased binding energy between the ligands and wall caused a corresponding linear increase in bonding ability. Our simulations also indicated that larger nanoparticles and those of rod shape with a higher aspect ratio have better binding ability than those of smaller size or rounder shape.

## Background

Nanoparticulate systems have been widely used for drug and gene delivery, imaging, and photodynamic therapy [1–12]. A typical nanoparticulate system consists of a nanoplatform, such as liposomes, polymeric micelles, quantum dots, nanoshells, or dendrimers, coated with ligands like hydrophobic drugs, DNA, or imaging agent. Ligands direct the nanoplatforms to specific locations and help to improve their bioavailability during circulation in a biological system [2, 3, 7–10, 13, 14]. Two main methods are used to transport ligand-coated nanoparticles (NPs) to diseased sites: passive and active targeting. In passive targeting, the accumulation of NPs is achieved by the enhanced permeability and retention effect [3, 7, 10, 15, 16] because the leaky vasculature and low lymphatic drainage prolong the residence time of NPs in the tumor. Conversely, active targeting is mediated by specific interactions between ligands that are connected via flexible spring tethers and receptors that are overexpressed at the pathological site. The highly concentrated receptors around pathological sites are preferred for ligand interaction because they can enhance NP internalization and retention [3, 7, 10, 15, 16].

Understanding the effects of NP size, hydrodynamic force, and multivalent interactions with a targeted biosurface on the mechanisms of a targeted delivery process is essential to aid the design and fabrication of NP systems [15]. Experimental techniques, such as fluorescence spectroscopy combined with microfluidics [17] and surface plasmon resonance [18], have been developed to investigate the ligand–receptor binding kinetics in vivo. The acquired experimental data indicate that the process of NP binding to a targeted surface is a synergic result of many factors, including the shape and diffusion of NPs, the flow effects [17, 19], as well as binding and internalization kinetics [20]. However, exploring this phenomenon experimentally is a very time-consuming task because of the small size of NPs and the dynamic nature of the transportation–deposition process; moreover, many details are difficult to capture because the binding process is very fast.

Therefore, theoretical modeling and numerical simulation have been performed to study the margination and adhesion processes of NPs in a fluid. For instance, Liu et al. [21] investigated the shape-dependent adhesion kinetics of non-spherical NPs through theoretical modeling. The influences of NP shape, ligand density, and shear rate on bonding ability under Brownian dynamics were systematically studied. They also investigated the distribution of NPs with different shapes and sizes in a mimetic branched blood vessel and found that NPs with smaller size and rod shape have better bonding ability [19].

Dissipative particle dynamics (DPD) simulations can precisely model hydrodynamic interactions at a mesoscopic scale with acceptable time scales [22, 23], which can overcome the limitations of molecular dynamics simulations [24, 25] to predict complex hydrodynamics with much higher efficiency. Although DPD was first introduced to simulate the dynamics of fluids [26–28], it has been successfully used to reproduce hydrodynamic forces [27], explore the phase behavior of lipid molecules [29], and study the interactions of biomembranes and NPs [30–33]. For example, Filipovic et al. [34] used DPD to simulate the motions of circular and elliptical particles in 2D shear flow and compared their results with those obtained from finite element (FE) calculations to validate the ability of the DPD method to model the motion of micro/nanoparticles at the mesoscale. They also combined the multiscale mesoscopic FE bridging procedure with DPD and the lattice Boltzmann method to model the motion of circular and elliptical particles in 2D laminar flow [35]. This approach proved to be an effective way to model the motion of NPs in drug delivery systems. Meanwhile, Ding et al. [36] studied the effects of the coating ligands on the cellular uptake of NPs and found that the strength of the receptor–ligand interaction along with the density, length, and rigidity of the ligand can markedly affect the final equilibrium in receptor-mediated endocytosis.

Despite these exciting advances, theoretical modeling using approaches such as Brownian adhesive dynamics can provide some insights into adsorption kinetics and the dynamics of adsorbed NPs but lacks specific details about the binding process [37, 38]. This paper presents the details of dynamic transportation and adhesion of NPs to a vascular wall under shear flow determined using DPD simulations. Parameters such as bonding time and the mean-square displacement of NPs are considered. Results obtained for spherical NPs with different binding forces and diameters and for NPs with different shapes or aspect ratios but the same volume are compared to assess the influence of such parameters on the binding of NPs to a vascular wall.

## Methods

## Coarse-Grained (CG) Model: DPD Simulation

To achieve targeted drug delivery, NPs are usually coated with polymers that specifically bind to a particular type of receptor on the vessel cell surface [37, 38]. It is computationally expensive to model the transportation and adhesion processes using an atomistic molecular dynamics simulation. However, the coarse-grained (CG) method guarantees that the general trend of the simulation will be determined without entirely erasing the structural details [39]. In this work, CG models were used to represent the components within the simulation system. The ligands on each NP were modeled as a chain of ten coarse-grain beads, namely, three hydrophobic beads and seven hydrophilic beads connected linearly to represent the polar head groups. As shown in Fig. 1, chains were uniformly distributed on the surface of a spherical NP. A harmonic potential was used to model the diblock copolymer chain, and the spring constant was set to 100 (reduced DPD units).

The size of the simulation box in our work was 22*r*
_{c} × 22*r*
_{c} × 22*r*
_{c} (*r*
_{c} is interaction length) with periodic boundary conditions in the *x* and *y* directions. The vascular surface was simplified as a fixed wall and placed at the boundary of the system consisting of fixed CG beads during the simulation. The wall was impenetrable with a “no slip” boundary condition where both the normal and tangential components of the particle momentum were inverted [40]. During the whole process, the wall particles did not move, acting as the location of the receptor that could interact with the ligands on the NPs. NPs with a diameter of 2 nm were also modeled by rigid beads placed in the middle of the box and filling the rest of the space with 27,783 explicit fluid particles [41]. The number densities *ρ* of the vascular wall and fluids were set as 3, as suggested elsewhere [42].

## Interaction Forces and Units in DPD Simulations

The interaction forces between different beads in the DPD formulation include a conservative force *f*
^{C}, a bead-spring force of the bonded monomers *f*
^{S}, a dissipative force *f*
^{D}, and a random force *f*
^{R} [43]:$\begin{array}{l}{f}_{ij}={f}_{ij}^{\mathrm{C}}+{f}_{ij}^{\mathrm{S}}+{f}_{ij}^{\mathrm{D}}+{f}_{ij}^{\mathrm{R}}\hfill \\ \phantom{\rule{1.12em}{0ex}}=\left[-{a}_{ij}\left({r}_{\mathrm{c}}-{r}_{ij}\right)-K\left({r}_{\mathrm{s}}-{r}_{ij}\right)-y{w}^{\mathrm{D}}\left(\overrightarrow{r}/{r}_{ij},{v}_{ij}\right)+\sigma {w}^{\mathrm{R}}{\left({\varsigma}_{ij}\left(\Delta t\right)\right)}^{-0.5}\right]\overrightarrow{r}/{r}_{ij},{r}_{ij}<{r}_{\mathrm{c}\text{,}}\hfill \end{array}$1

where *a*
_{ij} is the repulsion factor, *r*
_{ij} and *v*
_{ij} are the respective distance and velocity vectors of particle *i* with respect to particle *j*, *r*
_{c,} and *r*
_{s} are the respective cutoff distances for conservative and bead-spring forces, *K*, *γ*, and *σ* denote the spring constant, friction coefficient, and noise amplitude, respectively, *ω*
^{D} and *ω*
^{R} are the weight functions (*ω*
^{D} = (*ω*
^{R})^{2} = (*r*
_{c} − *r*
_{ij})^{2}), *ɛ*
_{ij} is the Gaussian random number, and Δ*t* is the simulation time step.

The random noise strength is expressed as a function of the dissipation strength and temperature *T* via the fluctuation–dissipation relation [28],${{\sigma}_{ij}}^{2}=2{\gamma}_{ij}{k}_{\mathrm{b}}T$2

where *σ*
_{ij} and *γ*
_{ij} are the random noise strength and dissipation strength between beads *i* and *j*, respectively. We carried out the simulations using a frictional coefficient *γ* of 3.

The default interactions between elements in DPD were described by the repulsion factor *a*
_{ii} = 25*k*
_{b}
*T*/*r*
_{c} to guarantee the compressibility of water at room temperature. Interaction factors between hydrophobic and hydrophilic beads were set to 45, while the others were set to 25 [44]. The repulsion factors between different elements used in the DPD simulation are listed in Table 1, where FE, HL, TL, VS, and WM denote the functional end of the ligand, the head of the ligand, the tail of the ligand, the vascular surface, and a water molecule, respectively.

To implement DPD simulation, *r*
_{c}, the bead mass *m*, and the thermostat temperature *k*
_{b}
*T* were set as unit elements [43, 45, 46]. All simulations were performed in the *NVE* ensemble with constant particle number *N*, simulation box volume *V*, and energy *E*. The velocity Verlet algorithm was used to integrate with a relatively small time step of Δ*t* = 0.02*τ*, and each simulation was run for 4 × 10^{5} steps.

## Simplified Shear Rate

To apply shear flow in the flow region, we employed the SLLOD algorithm [47, 48] using the following equations:$\frac{d{r}_{i,v}}{dt}=\frac{{P}_{i,v}}{{m}_{v}}+y{z}_{v}{\delta}_{i}\text{,}$3 $\frac{d{p}_{i,v}}{dt}={F}_{i,v}-\stackrel{\u2022}{y}{{p}_{z,}}_{v}{\delta}_{i}\text{,}$4

where *r*
_{i,v}, *p*
_{i,v}, and *m*
_{v} are the position vector, peculiar momentum, and mass of the *v*th bead, respectively, *y*
^{•} is the shear rate, and *δ*
_{i} is the unit vector in the *x* direction. This approach allows us to impose a linear velocity profile in the *x* direction with a constant gradient in the *z* direction.

## Statistical Analysis

We examined the significance of the data presented in Figs. 2, 3, and 4 below. All *P* values were <0.05, so these data are significant at 0.05 level, indicating that it is highly unlikely that these results would be observed under the null hypothesis, and bonding times for different *x* values are significantly different. Therefore, even though the error bars in these figures look wide, the results are reasonable.

## Results and Discussion

The bonding processes of NPs under different shear flows were simulated with the developed CG model. A typical NP bonding process is shown in Fig. 5. As illustrated in Fig. 5a, the functional ends of an NP sense the attraction force from the vascular surface, and the ligands start to move toward the wall surface. Then, the NP is attracted to the wall until it is firmly attached to it, as shown in Fig. 5b, c.

## Effect of Binding Energy on Ligand–Receptor Binding Kinetics

Experimental results have shown that binding energies from 5 to 35 *k*
_{b}
*T* strongly influence ligand–receptor binding kinetics [49]. In DPD simulation, the binding factor Δ*a* is defined as the difference between the ligand–receptor repulsion factor and the receptor–solvent repulsion factor. To simulate different attractive forces between ligands and the vascular surface, we varied the ligand–receptor repulsion factor while keeping the receptor–solvent repulsion factor constant. Here, Δ*a* = 5 indicates weak binding and the ligand–receptor repulsion factor is very close to that of the ligand with solvent, while Δ*a* = 25 indicates strong binding and the ligand–receptor repulsion factor is close to zero [44].

For relatively weak binding strength (5 < Δ*a* < 11), NPs mostly lingered in the middle of the flow domain and only a fraction of them established a stable contact with the receptor surface. For relatively strong binding strength (13 < Δ*a* < 25), bonding readily occurred and the binding time (from the beginning of the simulation to stable attachment) was measured. Figure 6 reveals that the probability of attachment initially increased linearly from about 10 to 30 % when 5 < Δ*a* < 9, and then increased abruptly to nearly 100 % when Δ*a* = 11. Figure 2 depicts ten different simulations run using various binding strengths. The mean bonding time decreased almost linearly as Δ*a* increased. The standard deviation of bonding time for each Δ*a* also decreased as Δ*a* increased. Therefore, NPs with a large bonding force have a higher probability of bonding and take less time to bond than those with a small bonding force. These results agree well with a previous report [44], which stated that when Δ*a* ≈ 12 or larger, any initial contact between ligands and a vascular surface leads to stable attachment.

## Effect of Shear Flow on the Binding Process

The physiological range of shear rate in blood flow is approximately 40–2000 s^{−1}, including flow within postcapillary venules, large arteries, and arterioles/capillaries [50]. In this paper, simulations were carried out for shear rates ranging from 0 to 2000 s^{−1} to study the NP bonding process under different shear flow conditions. The mean-square displacement (MSD) of 2-nm NPs under different shear rates was determined, and the results are shown in Fig. 7.

To simplify the analysis, we set Δ*a* = 13, which meant that the NPs would readily attach to the wall and the effects of the adhesion force can be minimized. In Fig. 7, the trajectory of NPs can be traced from their MSD. Initially, an NP moves randomly in the middle of the box under the influence of shear force and Brownian motion, and thus, the curve rises and falls at the outset. When the NP moves near the wall, it is attracted by the receptors, so it progresses to the wall and is bound to it, at which point the curve reaches the maximum MSD. After binding, the NP can still move because of the drive velocity [44] originating from the Brownian movement. The chains of each NP are long enough to allow reasonable vibration of the NP. Therefore, the curve subsequently fluctuates around the ultimate MSD.

Higher shear rate is reported to result in lower bonding possibility for NPs when the NP diameter is larger than 100 nm [51]. However, we found that the bonding efficiency (time needed for the NPs to reach equilibrium) did not decrease with increasing shear rate, as seen in Fig. 7. The bonding situation may be different for NPs with a diameter of 2 nm in fluid at a position 20 nm above the capillary wall, so we compared the bonding time and ultimate MSD of NPs for shear rates ranging from 0 to 2000 s^{−1}, as shown in Fig. 8. With increasing shear rate, the bonding time and ultimate MSD do not decrease accordingly. Instead, the spots appear randomly in relation to shear rate, which means that the shear rate has no effect on the bonding process in a certain area. Even when the shear rate was set to 0 s^{−1}, the bonding time was still longer than that in most cases with shear rate. This is because Brownian force outweighs the drag force and is the dominant force for NPs larger than 100 nm [51], a phenomenon that is even more obvious for smaller NPs. This is the reason for the heterogeneous binding time and MSD in Fig. 8. As a result, for NPs with a diameter of 2 nm, the bonding probability is not influenced by shear rate. NPs will firmly bond to the wall once in contact with it, and thus, bonding condition is dominated by diffusive process and independent of shear flow. We also simulated the behavior of NPs with diameters of 4 and 6 nm under the same conditions and obtained equivalent results.

## Effect of NP Size on the Binding Process

The Brownian force $\left({F}_{\mathrm{B}}\propto {R}^{\frac{1}{2}}\phantom{\rule{0.12em}{0ex}}\mathrm{f}\mathrm{o}\mathrm{r}\phantom{\rule{0.12em}{0ex}}\mathrm{round}\phantom{\rule{0.12em}{0ex}}\mathrm{nanoparticle}\phantom{\rule{0.12em}{0ex}}\mathrm{in}\phantom{\rule{0.12em}{0ex}}\mathrm{fluid}\right)$ will increase with the size of the NPs. Therefore, NPs with diameters of 4 and 6 nm were also simulated to study how NP size affects its bonding process. The size limitation of the simulation box meant that 6 nm was the largest size of NP we could investigate here. The bonding times of NPs with diameters of 2, 4, and 6 nm are illustrated in Fig. 9. The results indicate that the bonding time is shorter for larger NPs.

To allow quantitative analysis, the mean and standard deviations of bonding time for NPs of different sizes are presented in Fig. 3. Figure 3 reveals that the average bonding times are shorter and the standard deviations smaller for larger NPs. In our simulation, Brownian force is the determinative force, indicating that the motion of small particles in fluids is controlled by random collisions with surrounding fluid molecules. The random collision of NPs of smaller size is less likely to be balanced than that of larger NPs. If the unbalanced force is not in the direction of the wall, which is more likely to happen than the force being in the direction of the wall, the NP is less likely to be attracted to the wall and its bonding time will be prolonged. The track of a smaller NP is more disordered than that of a larger NP because of the unbalanced random force, which results in a larger standard deviation. Thus, the simulation results demonstrate that a bigger NP has higher bonding ability than a smaller one.

## Effect of NP Shape on the Binding Process

Both simulation and experimental results show that shape strongly affects pharmacokinetics and pharmacodynamics [52]. Different shapes cause contact areas, random forces of Brownian motion, and drag forces induced by shear flow to vary. Therefore, rod-shaped NPs (nanorods) were investigated in addition to spherical NPs to study the effect of NP geometry on the particle bonding process [38]. Nanorods with aspect ratios (*γ*; ratio of long axis to short axis) [51] of 5, 10, and 15 were considered. Nanorods with *γ* = 3 is shown in Fig. 10. The volume of the nanorods used was the same as that of the 2-nm NPs to ensure the same drug load capacity.

Figure 11 illustrates the bonding process of a nanorod. In Fig. 11a, the ligands on the end of the nanorod start to sense the attraction force from the vascular surface. Then, the ligands on the side of the nanorod are gradually absorbed by the wall until they are firmly attached to it, as depicted in Fig. 11b, c.

Figure 12 compares bonding times for NPs of different shapes. Bonding efficiency was higher for nanorods than spherical NPs, and increased with *γ*. Figure 4 reveals that the average bonding time and standard deviation were smaller for nanorods than spherical NPs. This is because of the tumbling motion and the larger contact area of the nanorods compared with the spherical NPs [38]. The contact area of a spherical NP is irrelevant to its orientation, and the binding area remains constant within interacting distance. Conversely, for nanorods, the contact area depends on orientation, and there is a higher chance to initiate contact with the wall because of its longer length compared with spherical NPs. Both the mean and standard deviations of bonding time decrease with increasing *γ* because NPs with a larger *γ* are thinner and longer, so the ligands on the end of the nanorod have a higher chance of interacting with the wall. These results are consistent with the finding that the strength of adhesion increases with *γ* [53], and thus, the bonding time decreases.

## Conclusions

We used DPD simulations to study the dynamics of polymerized NP binding with a vascular surface. We described in detail how the shear rate, bonding energy, size, and shape of an NP affect its bonding ability. The results indicate that the bonding ability increases linearly with bonding energy. Interestingly, the shear rate does not influence the bonding process of NPs with a diameter of 2–6 nm in a liquid environment 20 nm above the capillary wall. Compared with small spherical NPs, those with larger diameter or rod shape will move in a more orderly manner and require less time to reach the surface of the capillary wall, which means that they have better bonding efficiency. Additionally, the bonding ability of nanorods increased with *γ*. Our results provide some useful theoretical bases for designing NPs, which may aid in the development of new types of NPs with advantageous functionalities for biomedicine applications.

## Abbreviations

- CG
- DPD
- FE
- NP

Coarse-grained

Diffusive particle dynamics

Finite element

Nanoparticle

## Footnotes

**Competing Interests**The authors declare that they have no competing interests.**Authors’ Contributions**YL finished the model, analyzed the results, and acquired the original data in this article. BP, YHZ, LXY, and YLL made substantial contributions to the conception and design of this article. All the authors read and approved the final manuscript.**Authors’ Information**BP is the Vice Dean of the School of Mechatronics Engineering, University of Electronic Science and Technology of China, 2006, Xiyuan Ave, West Hi-Tech Zone, Chengdu, Sichuan 611731, People’s Republic of China.

## Acknowledgements

The authors would like to acknowledge the partial supports provided by the National Natural Science Foundation of China (Nos. 91123023, 51205047, and 51305066), the Fundamental Research Funds for the Central Universities (No. ZYGX2014Z004), and the National Youth Top-Notch Talent Support Program.

## References

- 1.
Allen, TM; Moase, EH: Therapeutic opportunities for targeted liposomal drug delivery, Adv Drug Deliv Rev, 1996, 21, 2, S. 117–33Crossref

- 2.
Hubbell, JA: Enhancing drug function, Science, 2003, 300, 5619, S. 595–6CrossrefPubMed

- 3.
Torchilin, VP: Targeted pharmaceutical nanocarriers for cancer therapy and imaging, AAPS J, 2007, 9, 2, S. E128–47CrossrefPubMed

- 4.
Murphy M, Ting K, Zhang X, Soo C, Zheng Z. Current development of silver nanoparticle preparation, investigation, and application in the field of medicine. J Nanomater. 2015. doi.org/10.1155/2015/696918Google Scholar

- 5.
Zhang, L; Gu, FX; Chan, JM: Nanoparticles in medicine: therapeutic applications and developments, Clin Pharmacol Therap, 2008, 83, 5, S. 761–9CrossrefPubMed

- 6.
Panyam, J; Labhasetwar, V: Biodegradable nanoparticles for drug and gene delivery to cells and tissue, Adv Drug Deliv Rev, 2003, 55, 3, S. 329–47CrossrefPubMed

- 7.
Byrne, JD; Betancourt, T; Brannon-Peppas, L: Active targeting schemes for nanoparticle systems in cancer therapeutics, Adv Drug Deliv Rev, 2008, 60, 15, S. 1615–26CrossrefPubMed

- 8.
Davis, ME; Shin, DM: Nanoparticle therapeutics: an emerging treatment modality for cancer, Nat Rev Drug Discov, 2008, 7, 9, S. 771–82CrossrefPubMed

- 9.
Sutton, D; Nasongkla, N; Blanco, E; Gao, J: Functionalized micellar systems for cancer targeted drug delivery, Pharm Res, 2007, 24, 6, S. 1029–46CrossrefPubMed

- 10.
McNeil, SE: Nanotechnology for the biologist, J Leukoc Biol, 2005, 78, 3, S. 585–94CrossrefPubMed

- 11.
Shi, C; Zhu, N; Cao, Y; Wu, P: Biosynthesis of gold nanoparticles assisted by the intracellular protein extract of Pycnoporus sanguineus and its catalysis in degradation of 4-nitroaniline, Nanoscale Res Lett, 2015, 10, 1, S. 1–8PubMed

- 12.
Liu, M; Feng, B; Shi, Y; Su, C; Song, H; Cheng, W: Protamine nanoparticles for improving shRNA-mediated anti-cancer effects, Nanoscale Res Lett, 2015, 10, 1, S. 1–7PubMed

- 13.
Goel, V; Pietrasik, J; Dong, H; Sharma, J; Matyjaszewski, K; Krishnamoorti, R: Structure of polymer tethered highly grafted nanoparticles, Macromolecules, 2011, 44, 20, S. 8129–35Crossref

- 14.
Spenley, N: Scaling laws for polymers in dissipative particle dynamics, EPL (Europhysics Letters), 2000, 49, 4, S. 534Crossref

- 15.
Hirsjarvi, S; Passirani, C; Benoit, J-P: Passive and active tumour targeting with nanocarriers, Curr Drug Discov Technol, 2011, 8, 3, S. 188–96CrossrefPubMed

- 16.
Danhier, F; Ucakar, B; Magotteaux, N; Brewster, ME; Préat, V: Active and passive tumor targeting of a novel poorly soluble cyclin dependent kinase inhibitor, JNJ-7706621, Int J Pharm, 2010, 392, 1, S. 20–8CrossrefPubMed

- 17.
Haun, JB; Hammer, DA: Quantifying nanoparticle adhesion mediated by specific molecular interactions, Langmuir, 2008, 24, 16, S. 8821–32CrossrefPubMed

- 18.
Xia, J; Zhong, C: Dissipative particle dynamics study of the formation of multicompartment micelles from ABC star triblock copolymers in water, Macromol Rapid Commun, 2006, 27, 14, S. 1110–4Crossref

- 19.
Goldstein, B; Coombs, D; He, X; Pineda, AR; Wofsy, C: The influence of transport on the kinetics of binding to surface receptors: applications to cells and BIA core, J Mol Recognit, 1999, 12, 5, S. 293–9CrossrefPubMed

- 20.
Wilhelm, C; Gazeau, F; Roger, J; Pons, J; Bacri, J-C: Interaction of anionic superparamagnetic nanoparticles with cells: kinetic analyses of membrane adsorption and subsequent internalization, Langmuir, 2002, 18, 21, S. 8148–55Crossref

- 21.
Shah, S; Liu, Y; Hu, W; Gao, J: Modeling particle shape-dependent dynamics in nanomedicine, J Nanosci Nanotechnol, 2011, 11, 2, S. 919–28CrossrefPubMed

- 22.
Huang, J; Wang, Y: Control of aggregation of nanoparticles by double-hydrophilic block copolymers: a dissipative particle dynamics study, J Phys Chem B, 2007, 111, 27, S. 7735–41CrossrefPubMed

- 23.
Groot, RD; Madden, TJ: Dynamic simulation of diblock copolymer microphase separation, J Chem Phys, 1998, 108, S. 8713Crossref

- 24.
Peng, B; He, W; Hao, X: Interfacial thermal conductance and thermal accommodation coefficient of evaporating thin liquid films: A molecular dynamics study[J], Computational Materials Science, 2014, 87, S. 260–266Crossref

- 25.
Li, Y; Guo, ZY; Peng, B: Buckling Behaviors of Imperfect Single-Walled Carbon Nanotubes: A Molecular Dynamic Simulation[J], Applied Mechanics and Materials, 2012, 110, S. 3831–3837

- 26.
Groot, RD; Warren, PB: Dissipative particle dynamics: bridging the gap between atomistic and mesoscopic simulation, J Chem Phys, 1997, 107, 11, S. 4423Crossref

- 27.
Hoogerbrugge, P; Koelman, J: Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics, EPL (Europhysics Letters), 1992, 19, 3, S. 155Crossref

- 28.
Espanol, P; Warren, P: Statistical mechanics of dissipative particle dynamics, EPL (Europhysics Letters), 1995, 30, 4, S. 191Crossref

- 29.
Rodgers, JM; Sorensen, J; de Meyer, FJ; Schiott, B; Smit, B: Understanding the phase behavior of coarse-grained model lipid bilayers through computational calorimetry, J Phys Chem B, 2012, 116, 5, S. 1551–69CrossrefPubMed

- 30.
Smith, KA; Jasnow, D; Balazs, AC: Designing synthetic vesicles that engulf nanoscopic particles, J Chem Phys, 2007, 127, 8, S. 084703CrossrefPubMed

- 31.
Yue, T; Li, S; Zhang, X; Wang, W: The relationship between membrane curvature generation and clustering of anchored proteins: a computer simulation study, Soft Matter, 2010, 6, 24, S. 6109–18Crossref

- 32.
Yue, T; Zhang, X: Molecular understanding of receptor-mediated membrane responses to ligand-coated nanoparticles, Soft Matter, 2011, 7, 19, S. 9104–12Crossref

- 33.
Yang, K; Ma, Y-Q: Computer simulation of the translocation of nanoparticles with different shapes across a lipid bilayer, Nat Nanotechnol, 2010, 5, 8, S. 579–83CrossrefPubMed

- 34.
Filipovic, N; Kojic, M; Ferrari, M: Dissipative particle dynamics simulation of circular and elliptical particles motion in 2D laminar shear flow, Microfluid Nanofluid, 2011, 10, 5, S. 1127–34Crossref

- 35.
Filipovic, N; Isailovic, V; Đukić, T; Ferrari, M; Kojic, M: Multiscale modeling of circular and elliptical particles in laminar shear flow, IEEE Trans Biomed Eng, 2012, 59, 1, S. 50–3CrossrefPubMed

- 36.
H-m, D; Y-q, M: Role of physicochemical properties of coating ligands in receptor-mediated endocytosis of nanoparticles, Biomaterials, 2012, 33, 23, S. 5798–802CrossrefPubMed

- 37.
English, TJ; Hammer, DA: Brownian adhesive dynamics (BRAD) for simulating the receptor-mediated binding of viruses, Biophys J, 2004, 86, 6, S. 3359–72CrossrefPubMed

- 38.
Liu, Y; Tan, J; Thomas, A; Ou-Yang, D; Muzykantov, VR: The shape of things to come: importance of design in nanotechnology for drug delivery, Ther Deliv, 2012, 3, 2, S. 181–94CrossrefPubMed

- 39.
Spaeth, JR; Kevrekidis, IG; Panagiotopoulos, AZ: Dissipative particle dynamics simulations of polymer-protected nanoparticle self-assembly, J Chem Phys, 2011, 135, 18, S. 184903CrossrefPubMed

- 40.
Pivkin, IV; Karniadakis, GE: A new method to impose no-slip boundary conditions in dissipative particle dynamics, J Comput Phys, 2005, 207, 1, S. 114–28Crossref

- 41.
Huang, M-J; Kapral, R; Mikhailov, AS; Chen, H-Y: Coarse-grain model for lipid bilayer self-assembly and dynamics: multiparticle collision description of the solvent, J Chem Phys, 2012, 137, 5, S. 055101CrossrefPubMed

- 42.
Qian H-J, Chen L-J, Lu Z-Y, Li Z-S. Surface diffusion dynamics of a single polymer chain in dilute solution. Phys Rev Lett. 2007;99(6). doi:10.1103/PhysRevLett.99.068301.Google Scholar

- 43.
Duong-Hong, D; Phan-Thien, N; Fan, X: An implementation of no-slip boundary conditions in DPD, Comput Mech, 2004, 35, 1, S. 24–9Crossref

- 44.
Djohari, H; Dormidontova, EE: Kinetics of nanoparticle targeting by dissipative particle dynamics simulations, Biomacromolecules, 2009, 10, 11, S. 3089–97CrossrefPubMed

- 45.
Cui, Y; Zhong, C; Xia, J: Multicompartment micellar solutions in shear: a dissipative particle dynamics study, Macromol Rapid Commun, 2006, 27, 17, S. 1437–41Crossref

- 46.
Chen, C; Gao, C; Zhuang, L; Li, X; Wu, P; Dong, J: A many-body dissipative particle dynamics study of spontaneous capillary imbibition and drainage, Langmuir, 2010, 26, 12, S. 9533–8CrossrefPubMed

- 47.
Evans, DJ; Morriss, G: Statistical mechanics of nonequilibrium liquids, 2008, Cambridge, University Press

- 48.
Kalra, V; Escobedo, F; Joo, YL: Effect of shear on nanoparticle dispersion in polymer melts: a coarse-grained molecular dynamics study, J Chem Phys, 2010, 132, 2, S. 024901CrossrefPubMed

- 49.
Moore, NW; Kuhl, TL: The role of flexible tethers in multiple ligand-receptor bond formation between curved surfaces, Biophys J, 2006, 91, 5, S. 1675–87CrossrefPubMed

- 50.
Ghitescu, L; Bendayan, M: Immunolabeling efficiency of protein A-gold complexes, J Histochem Cytochem, 1990, 38, 11, S. 1523–30CrossrefPubMed

- 51.
Tan, J; Shah, S; Thomas, A; Ou-Yang, HD; Liu, Y: The influence of size, shape and vessel geometry on nanoparticle distribution, Microfluid Nanofluid, 2013, 14, 1–2, S. 77–87CrossrefPubMed

- 52.
Tao, L; Hu, W; Liu, Y; Huang, G; Sumer, BD; Gao, J: Shape-specific polymeric nanomedicine: emerging opportunities and challenges, Exp Biol Med, 2011, 236, 1, S. 20–9Crossref

- 53.
Decuzzi, P; Ferrari, M: The adhesive strength of non-spherical particles mediated by specific interactions, Biomaterials, 2006, 27, 30, S. 5307–14CrossrefPubMed