Veronika
Schubertová
a,
Francisco J.
Martinez-Veracoechea
b and
Robert
Vácha
*c
aFaculty of Science, Masaryk University, Kotlárská 2, 611 37 Brno, Czech Republic
bInfochem-KBC, 23 Queen Elizabeth Street, London SE1 2LP, UK
cFaculty of Science and CEITEC – Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic. E-mail: robert.vacha@mail.muni.cz; robertvacha@gmail.com
First published on 16th February 2015
Cellular uptake is a crucial process in nanomedicine and drug-delivery; however, the factors that affect its efficiency/speed are not well understood. We report computer simulations on passive uptake via receptor-mediated endocytosis of nanoparticle coated with ligands. In particular, we study how the distribution of ligands on the nanoparticle surface influences the uptake rate. The speed of membrane wrapping and uptake was found to be the fastest for nanoparticles with homogeneous ligand distributions, where ligands are spread evenly on the surface. We show that the diffusion of the ligands on the nanoparticle can hinder its uptake, since upon the interaction with the membrane the ligand distribution becomes extremely inhomogeneous, with a large ligand-free patch. When the ligand-free-area was more than 20% of the surface, we did not observe uptake within the scale of our simulations.
Many nanoparticles are inspired by virus capsids, where ligand distribution corresponds to the surface distribution of receptor binding sites, also known as virus attachment protein positions. Such distribution follows the symmetry of the virus capsid but is generally not known. A few examples have been recently revealed by cryo-electron microscopy. The virus family Picornaviridiae, which includes the Rhinovirus and Poliovirus species, has 60 binding sites located between the binate and the quintet symmetry axis on an icosahedral capsid.16 On the Canine Parvovirus capsid the receptor binds to the center of the face on the capsid.17
As we are interested in the generic aspects of how the ligand distribution influences endocytosis, we use a simplified coarse grained model. We employed icosahedral shape particles, because it is the most common shape of virus capsids and is close to spherical nanoparticles. The effect of nanoparticle shape has been investigated before.4,18–20 We performed Molecular Dynamics (MD) simulations of nanoparticles in contact with a tensionless phospholipid membrane and studied the uptake process with focus on the influence on the uptake time by the concentration and distribution of binding sites on the nanoparticle surface.
Nanoparticles were made of beads of the same size as phospholipids and are mostly hydrophilic (i.e., purely repulsive). Specific numbers of beads (ligands) were made attractive to the membrane receptors (short range attraction with 8 kT minimum) and the positions of these beads were varied to represent different ligand distributions. The nanoparticle shell was made of 792 beads resulting in a diameter of about 13 nm. The choice of the model and its parameters is based on our previous studies where we investigated the uptake of nanoparticles with random homogeneous ligand distributions affected by the nanoparticle size, shape, total number of ligands, and strength of the ligand–receptor interaction. More details about the model and the uptake results can be found in ref. 4 and 9. Note that a tension-less membrane is a valid assumption to study nanoparticles of sizes up to a few hundred nanometers, beyond which the cellular tension (in order 3 × 103 to 1 dyn cm−1)23 will start to play a role.24 The simulation time unit τ can be related to 100 ns based on the lipid diffusion; however, this should be considered as a very rough estimate since hydrodynamic effects are missing in our implicit solvent model.
In Fig. 2 we plot the lateral dimension of the simulation box versus time for the different nanoparticle patterns. The lateral dimension of the simulation box (either x or y) is a good indicator of the amount of wrapping. This is because as the nanoparticle is wrapped with the membrane, the cross-section of the simulation box must diminish in order to satisfy the zero-tension constrain.
During our simulations we observed three types of uptake behaviors: (a) spontaneous direct uptake without metastable states; (b) uptake with metastable states; (c) no uptake (see Fig. 2).
Nanoparticles with a small (up to 72) number of ligands that cover up to 8% of its surface (V6F0E0, V0F0E2, V2F0E0, and V0F1E0) did not pass through the membrane and remained attached to the membrane within the timescale of our simulation. This means that nanoparticles were attached for more than 50000τ and ended with different levels of wrapping. Namely, the distributions V6F0E0, V0F0E2, V2F0E0, and V0F1E0 lead to about 50%, 50%, 25%, and 10% respectively nanoparticle wrapping at the end of the simulation (degree of wrapping was estimated visually from the surface area in contact with the membrane and from the ligand–receptor interaction energy).
Nanoparticles with more than 100 ligands, which is roughly one seventh of the nanoparticle surface, typically were uptaken fast without any metastable states (green points in Fig. 1B). This holds for distributions named V6F0E7, V6F7E0, V0F10E0, V16F0E0, V1F7E0, V0F0E4, and V4F3E0. The time needed for the uptake ranged from 1000τ to 5500τ with distributions V6F0E7 and V6F7E0 being the fastest. The estimated uptake time error is about 500–2000τ (see ESI†). The larger error (2000τ) is for simulations with metastable states due to their stochastic behavior. The error could be decreased by the repetition of simulations, which has not been carried out due to computational expenses.
The most interesting case was those nanoparticles with a medium number of ligands, for which the effect of ligand distribution was the most significant (red points in Fig. 1B). The effect was so dramatic that even though all the nanoparticles had roughly 80 ligands, the uptake behaviour could vary from fast uptake (below 10000τ) to no uptake at all within the whole range of our simulation (50
000τ).
The ligand distributions were more efficient when located at edges than at faces and the least efficient were distributions at vertices. This correlates well with the increasing number of position symmetry (2, 3, and 5) or the decreasing number of places (icosahedron has 30 edges, 20 faces, and 12 vertices). This clearly suggests that the most homogeneous distributions with ligands evenly spread are the most efficient. Indeed, the distributions V1F7E0 and V4F3E0 with 152 and 108 ligands had smaller uptake time (by 2000τ) than distributions V0F10E0 and V16F0E0 with 200 and 192 ligands, respectively (see the inset in Fig. 1B). Thus, it is demonstrated that on certain circumstances ligand distribution can be more critical than the number of ligands in determining the rate of uptake. The effect of ligand distribution is expected to be enhanced for larger particles and/or lower receptor concentrations.
To rationalize why a homogenous ligand distribution has a faster rate of uptake, we need to remember that the uptake process is a compromise between the energy increase due to bending of the membrane and the energy decrease due to ligand–receptor interactions.4 Once the first ligand–receptor contacts are formed, further wrapping is driven by the energy gain associated with further ligand–receptor contacts. However, if ligands are not distributed homogenously on the nanoparticle surface (i.e., they are highly concentrated on a few localized points which are far apart from each other) the wrapping process must incur a high penalty in bending energy before it can be rewarded by achieving new ligand–receptor contacts. Therefore, in this case, the wrapping process has a high activation energy, which consequently hinders the uptake process. However, wrapping can still occur via thermal fluctuations. We have observed two types of such fluctuations, particularly visible in longed lived metastable states: (1) fluctuations of the membrane shape and (2) nanoparticle rotations (see Fig. S2†). Naturally, larger fluctuations have higher activation energy and are therefore less likely to happen. As a result, heterogeneous distributions have a slower uptake rate via receptor-mediated endocytosis. As an example of homogeneous and heterogeneous distributions with about the same number of ligands see distributions V0F1E2 and V6F0E0, respectively in Fig. 1C.
In many applications, ligands may not be fixed on the nanoparticle surface as employed above, but instead they are mobile. Diffusion of ligands on the nanoparticle surface can dynamically change ligand distributions. We tested the effect of ligand mobility by replacing the nanoparticle with a small vesicle, where half of the lipid headgroups were modified as ligands, which remained free to diffuse as lipids. It has been shown that such vesicles can interact with a membrane in various ways.5 For our purpose we focused on a path where the membrane wraps around the vesicle as it was a hard nanoparticle (see Fig. 3). Most of the ligands diffused into contact with the membrane (i.e., the vesicle became polarized), creating a part of the vesicle (a patch) that was almost ligand free, seemingly hindering the uptake.
In order to avoid the effects of vesicle flexibility5,25 and further investigate the influence of ligand free patches on uptake efficiency of nanoparticles, we studied particles with three different diameters (5, 8, and 12 nm), where a ligand-free-area with 5-fold symmetry was constructed. Within the timescale of our simulations nanoparticles with a ligand-free area larger than 20% of the surface did not undergo the uptake. In the case of the nanoparticle with a diameter of only about 5 nm, already 10% of the ligand-free-area prevented the uptake (see Fig. 4). Note that the nanoparticles always rotated in such a way that the patch was wrapped last, if at all. This rotation is in line with the reported change of nanoparticle orientation during the membrane wrapping, which was shown to finish with the least energetically favorable part.4,18,20,26,27 This part can be calculated from the difference between the local membrane bending and receptor–ligand binding energy.
Therefore, nanoparticles where ligands can diffuse on the particle surface might be inefficient for the uptake unless almost the whole (more than 80%) of the surface is covered with ligands. As shown above the ligand-free patch on the nanoparticle can inhibit or significantly slow down the uptake. Note that in practice the polarization, i.e. accumulation of ligands on one side, can be hindered via slower diffusion (kinetically) or repulsion between the ligands (thermodynamically).
Our results suggest that nanoparticles where the ligands are all distributed on one side can be inefficient for the uptake unless almost the whole (more than 80%) of the surface is covered with ligands. This is in sharp contrast with only 10% ligands needed when the distribution is homogeneous.
We also found that ligand-free patches slow down the uptake, which can be completely prevented if the ligand-free area is large enough. Within the nanometer size of our nanoparticles and the timescale of our simulations, the threshold area for preventing uptake was only weakly dependent on the radius, and was about 20% of the surface. It means that nanoparticles where the ligands are all distributed on one side can be inefficient for the uptake unless more than 80% of the surface is covered with ligands. This is in sharp contrast with only 10% ligands needed when the distribution is homogeneous demonstrating the importance of ligand distribution. This has important consequences also for nanoparticles with ligands that can diffuse on the surface, as ligand-free patches can spontaneously occur upon interactions with membrane receptors.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4sm02815e |
This journal is © The Royal Society of Chemistry 2015 |