# Target shape dependence in a simple model of receptor-mediated endocytosis and phagocytosis

See allHide authors and affiliations

Edited by William Bialek, Princeton University, Princeton, NJ, and approved April 19, 2016 (received for review November 10, 2015)

## Significance

Most previous work on endocytosis (the way cells ingest particles) has focused on engulfment of spherical targets. However, this is rarely relevant outside the laboratory, where cells must engulf a huge variety of target shapes, from capped cylinders (like many bacteria) to hourglasses (such as dividing cells). Here, we map this problem to the well-studied physics problem of the freezing of water. This allows us to study a wide range of different, biologically relevant shapes. In particular, for the first time, to our knowledge, it allows us to examine the orientation dependence of phagocytosis, to explain the shape dependence of phagocytosis and receptor-meditated endocytosis, and to explain why nonspherical targets often engulf faster when the most highly-curved tip is engulfed first.

## Abstract

Phagocytosis and receptor-mediated endocytosis are vitally important particle uptake mechanisms in many cell types, ranging from single-cell organisms to immune cells. In both processes, engulfment by the cell depends critically on both particle shape and orientation. However, most previous theoretical work has focused only on spherical particles and hence disregards the wide-ranging particle shapes occurring in nature, such as those of bacteria. Here, by implementing a simple model in one and two dimensions, we compare and contrast receptor-mediated endocytosis and phagocytosis for a range of biologically relevant shapes, including spheres, ellipsoids, capped cylinders, and hourglasses. We find a whole range of different engulfment behaviors with some ellipsoids engulfing faster than spheres, and that phagocytosis is able to engulf a greater range of target shapes than other types of endocytosis. Further, the 2D model can explain why some nonspherical particles engulf fastest (not at all) when presented to the membrane tip-first (lying flat). Our work reveals how some bacteria may avoid being internalized simply because of their shape, and suggests shapes for optimal drug delivery.

Cells are capable of ingesting a huge range of particle shapes, from rod-shaped *Escherichia coli* to doubly lobed budding yeast, from helical *Borrelia* to filamentous *Legionella* (1⇓–3). Further, the shape is critical in determining whether engulfment is possible and, if so, how long it takes. In fact, it has been argued that the target shape plays an even more critical role than size (4, 5). In addition, target orientation strongly affects internalization, with ellipsoidal particles engulfed far more easily when the highly curved tip is presented first (4, 6). However, despite this, there is relatively little work studying shape and orientation dependence, with the vast majority of previous papers, both experimental and theoretical, involving only spherical targets.

Endocytosis encompasses a number of different mechanisms of cellular internalization, including clathrin-mediated endocytosis, caveolar-type endocytosis, macropinocytosis, and phagocytosis (7). Although sharing some similarities, these processes often proceed in strikingly different manners. For example, phagocytosis is typically highly active and involves the membrane extending outward, with finger-like protrusions surrounding the target in an actin-dependent process (8, 9). Conversely, clathrin-mediated endocytosis is a more passive process, with targets appearing to sink into the cell (10). Here we focus on types of endocytosis that involve target recognition via receptors. We assume that receptors bind irreversibly to ligands on the target so that engulfment proceeds monotonically (11, 12). In particular, we distinguish phagocytosis from other less active forms of receptor-driven internalization. We refer to all these latter processes as receptor-mediated endocytosis, including clathrin-mediated endocytosis.

Various mathematical models have attempted to understand the mechanism and dynamics of engulfment. All forms of endocytosis (in particular, phagocytosis) are extremely complex, involving hundreds of different protein species in signaling cascades and cytoskeleton remodeling (13⇓⇓–16). As a result, overly complicated models, that try to include every component, are unlikely to be useful. Instead, more progress can be made by considering simpler, intuitive models that capture the essential mechanisms. For example, viral engulfment has been modeled by examining the free energy of membrane bending and adhesion (17, 18). Similarly, van Effenterre et al. (19) used a thermodynamic approach that involved considering an ensemble of target particles. Recently, the role of actin during endocytosis in yeast was addressed by using a variational approach and arguing that the final pinch-off stage is due to a pearling-like instability (20).

There are also models that focus on phagocytosis. For example, van Zon et al. (6) included simple dynamics for both actin and receptors, which they used to understand why engulfment normally either stalls before halfway or reaches completion. Herant et al. (21, 22) focused on the forces required to explain the cup shape and showed, using a continuum mechanics approach, the need for both repulsion at the cup edge and flattening within the cup. Various other approaches focus only on energetic requirements, such as those due to membrane bending and receptor−ligand binding, equating the phagocytic cup with some minimum energy state. For example, Dasgupta et al. (23, 24) used this approach to argue that ellipsoids are harder to engulf than spheres, whereas Tollis et al. (25) found that an actin-driven ratchet mechanism can lead to robust engulfment.

One particularly elegant approach by Gao et al. (26) models endocytosis by considering only the motion of receptors within the membrane. They argued that the essence of receptor-mediated endocytosis is related to the dynamics of the receptors themselves, which can be mapped to the supercooled Stefan problem, a simple physical model of how the boundary between ice and water moves during freezing. Although they were able to understand how particle radius affects the rate of engulfment (predicting an optimum radius corresponding to the quickest possible engulfment time), they only considered spherical particles. Cylindrical ellipsoids were considered by Decuzzi et al. (27), but, as we will show, not in a consistent manner. In ref. 28, we extended the Gao et al. model to phagocytosis, arguing that a similar mathematical model can explain receptor motion in both cases. This involved including receptor drift and signaling, and showed that the two distinct stages of engulfment (an initial slow stage followed by a much quicker second stage) can be explained by different receptor dynamics during different stages. However, again, we almost exclusively focused on only spherical particles.

Here, to address the role of particle shape during engulfment, we extend these models in three ways. First, by allowing nonconstant curvature, we are able to study nonspherical particle shapes, such as ellipsoids, capped cylinders, and hourglasses (Fig. 1). All of the shapes we consider are relevant for real biological systems, such as rod-shaped bacteria and dividing budding yeast. Second, we include signaling and a role for actin to examine the different target shape dependencies between phagocytosis and other types of receptor-mediated endocytosis. Third, we extend the formalism from a 1D to a 2D model, which allows, for the first time to our knowledge, lower-symmetry shapes, such as lying ellipsoids, to be studied. With this model we are able to address how orientation affects engulfment, finding an explanation for why prolate spheroids engulf fastest when presented to the cell tip-first. Finally, we compare our model with known experimental results.

## Results

We first study endocytosis in general and consider particles both with and without circular symmetry. Then we extend the model to apply to phagocytosis by including signaling and actin. By particle symmetry, we mean the symmetry of the particle looking from directly above the membrane. For example, a prolate spheroid will have circular symmetry only when presented tip-first (as in the shape in Fig. 1*B*, *Left*). As we now show, our model for circularly symmetric particles becomes effectively 1D, allowing easier analytic and numerical treatment. Conversely, nonsymmetric particles require a full 2D model, introducing extra physical and numerical issues. We first consider circularly symmetric particles—spheres, tip-first spheroids, capped cylinders, hourglasses—and next examine nonsymmetric particles, for example, ellipsoids that lie flat (as in the shape in Fig. 1*B*, *Left*, lying on its side).

### A 1D Model for Endocytosis of Circularly Symmetric Particles.

We first extend the model in Gao et al. (26) to include nonspherical particles. This is a model of receptor motion with the inclusion of membrane bending and ligand−receptor binding, which focuses on the cup size as a function of time. Consider a particle whose cross-section parallel to the membrane is circularly symmetric (such as all of the shapes in Fig. 1). Then receptor motion within the membrane and the fraction of engulfment can be parameterized by just one spatial parameter, *r*, the distance from the center of the cup. Thus, at time *t*, the receptors are described by a density *A*).

Initially, *ρ* is assumed uniform and given by *B*). We assume, for simplicity, that receptors are never destroyed or created, so that the evolution of the receptor density, *ρ*, in the nonengulfed region (*D*. At the cup boundary, receptors flow into the bound region, increasing the engulfment *a*. By considering the conservation of receptors, the rate of increase of *a* can be related to the flux of receptors. Full details are given in *SI Appendix*, where we show that our 1D model is given in the nonengulfed region

To find a unique solution, we impose one extra condition, the value of *ρ* at the edge of the cup, *SI Appendix*, we consider three contributions to the free energy (receptor−ligand binding, membrane curvature, and receptor entropy), which leads to the condition*H* is the mean curvature of the particle at the cup edge. We ignore any contribution from the Gaussian curvature because, as we explain in *SI Appendix*, its effect is normally similar to that from the mean curvature for biologically relevant target shapes.

The original model in ref. 26 only considered spherical particles, so that the curvature *H* is constant. To consider other particle shapes, we allow *H* to depend on the engulfment *a*. This, in turn, means that *SI Appendix*). To make progress, we numerically solve the system (see *Materials and Methods*) and consider, in turn, spheroids, capped cylinders, and hourglasses.

### Spheroids.

A spheroid is an ellipsoid where two of the principal axes have the same length. Later, in our 2D model, we will be able to consider spheroids presented to the cell both end-on and lying down. However, in the 1D model, which requires circularly symmetric particles, only standing spheroids can be examined. A spheroid is characterized by two radii: *A*). We use a standard parameterization of the surface of a spheroid given by *z* axis represents the perpendicular distance from the membrane, so that

We must relate the mean curvature *H* to the arc length *s* (measured from *v* as*SI Appendix*),*v* cannot be directly eliminated between Eqs. **3** and **4**, we proceed as follows: For a given engulfment (arc length) *a*, we numerically solve Eq. **4** to find *v* and then use Eq. **3** to find the curvature.

We find that, for *B*). The speed of engulfment now depends on the local curvature at the cup edge. Higher curvatures decrease the right-hand side of Eq. **2** and increase *B*). In our mathematical model, this occurs when **1** becomes negative, which we interpret as

To compare different spheroids, we vary *SI Appendix*). As shown in Fig. 3 *B* and *C*, we find five classes of behavior, depending on the value of

To better understand the effect of changing *D*). Spheroids with *B*, *Left*) always take longer than the equivalent sphere to reach both half and full engulfment. Conversely, spheroids with

To study the orientation dependence of engulfment, we compare ellipsoids with {*R*_{1}, *R*_{2}} = {*R*, 1 μm} to those with {*R*_{1}, *R*_{2}} = {1 μm, *E*). For *R* > 1 μm, we find that the oblate spheroid (with {*R*_{1}, *R*_{2}} = {*R*, 1 μm}) engulfs fastest. This is simply because the total engulfment length is smaller. For the same reason, for *R* slightly smaller than 1 μm, engulfment is also quickest for oblate spheroids (which are now those with {*R*_{1}, *R*_{2}} = {1 μm, *R*, the required value of *R*_{1}, *R*_{2}} = {1 μm,

### Capped Cylinders and Hourglasses.

We now briefly consider two other circularly symmetric shapes; full details can be found in *SI Appendix*. Capped cylinders, consisting of a cylinder with a hemisphere attached at either end, are good models for various rod-shaped bacteria, such as *E. coli* and *Bacillus subtilis*. These display a slow−fast−stall−slow engulfment behavior, related to the lower curvature found in the neck region. The stalling is caused by a jump in

Hourglasses are more complicated shapes, but are useful for modeling cases of dividing cells, such as budding yeast during division (1, 30). The engulfment behavior and the total engulfment time depend on exactly how narrow the neck is compared with the spherical caps, with the neck region engulfed sometimes faster and sometimes slower than the caps. Interestingly, the quickest engulfed hourglass does not correspond to the shape with the smallest neck curvature. See *SI Appendix* for full details.

### Extension to Phagocytosis.

Although actin plays a role in many types of endocytosis, it seems to be far more important for successful engulfment in phagocytosis (25, 31, 32). To this end, we extend our model to apply to phagocytosis by adding a simple role for actin. Motivated by our previous work (28), we achieve this by including a signaling molecule *S* which recruits actin to push at the cup edge. We assume that *S* is produced only within the cup (with rate *τ*), and diffuses (with constant *S* is described by (28)

Even though the actin network provides a pushing force at the cup edge, this will not necessarily increase the engulfed target length, *a*, unless sufficient receptors are also present: simply pushing the membrane further around the target without a suitable zipper mechanism risks subsequent unwrapping (25). Thus, the real effect of actin (at least as far as our model is concerned) is in facilitating receptor motion to the cup edge. One way this could be achieved is by actin helping to bend the membrane around the target, thus effectively reducing the free energy needed for membrane bending. For example, an actin network structure that matches the target curvature would naturally encourage the membrane to bend in the same manner via scaffolding (33). We implement this mechanism by allowing the bending modulus, *S* at the cup edge,

We first examine the difference between receptor-mediated endocytosis and phagocytosis for oblate spheroids (*D*). However, with the introduction of signaling and a role for actin, such ellipsoids can often be ingested by phagocytosis (Fig. 4*A*). In these cases, the initial rate of engulfment, when *S* is small, is similar in both models. However, when engulfment slows down around halfway, the recruitment of actin to the cup edge provides an extra push. If the bending-actin constant is large enough, this then leads to complete engulfment (Fig. 4*A*). The situation is markedly different for prolate spheroids (

For oblate spheroids, we can examine the effect of the bending-actin constant on engulfment time in cases where nonphagocytic endocytosis stalls around halfway (Fig. 4*A*). Initially, as ^{2} with the parameters in Fig. 4*A*), complete engulfment occurs. As

Finally, we examine how the range of ingestible spheroids changes as the bending-actin constant increases (Fig. 4*B*). For any fixed value of

### A 2D Model for Noncircularly Symmetric Particles.

A limitation of the 1D model is that only circularly symmetric particles can be studied. In particular, this excludes general ellipsoids (with all of the axes different lengths) and even spheroids that are not presented tip-first. To study such shapes, we extend the above model to two dimensions, which involves studying a 2D Stefan problem. We parameterize the surface by polar coordinates *r* and *θ*, with *a* can be derived from the conservation of receptors. As we show in *SI Appendix*, the second part of Eq. **1** is replaced by**2**. (see *SI Appendix*). To find solutions, we developed a numerical approach that involves a 2D lattice that must be continually updated as the cup enlarges (see *Materials and Methods*). To test our numerics, we checked that spherical particles (with constant curvature

We focus on studying ellipsoids, which are parameterized by their three semiprincipal axes, *A*). Because the cup does not proceed in unison around the particle (different angles engulf at different rates), there is now no unique definition of half engulfment. We consider two definitions: half-area engulfment when the engulfed area is one-half the ellipsoid surface area, and half-circumference engulfment when all points of the cup are over halfway to the top of the particle. Half-area engulfment always occurs before (or at the same time as) half-circumference engulfment.

As an example of engulfment, consider an ellipsoid with *B*, *Left*, if it were lying on its side (see Movie S1). The lack of circular symmetry (when viewed from above) means that engulfment now proceeds at different rates at different angles. In particular, consider

As seen in Fig. 5*B* and Movie S1, the *θ* direction) causes even more receptors to move away from *SI Appendix*, Fig. S6). This amplification explains why engulfment takes so much longer at

For the first time, to our knowledge, we are now able to compare spheroids presented to the cell with different orientations. To do this, we compare spheroids with radii *R*, we always choose *C* shows the results for the total engulfment time. Both symmetrically and asymmetrically presented spheroids demonstrate the expected characteristic behavior, with quickest engulfment corresponding to some intermediate spheroid: More squashed or more pointed ellipsoids take longer to engulf, with sufficiently eccentric shapes never reaching complete engulfment. Further, the symmetrically presented spheroid always engulfs first, both for prolate and oblate spheroids. This is due to the amplification of angular diffusion described above, where the

For symmetrically presented spheroids, we find, as we should, similar results to those for spheroids in the 1D model. For fixed surface area, the sphere is not the optimal shape. Rather, slightly oblate spheroids (*SI Appendix*, Fig. S5*B*). However, the situation is different for the half-area engulfment time, where prolate spheroids (*SI Appendix*, Fig. S5*A*).

Thus, our 2D model leads to three important conclusions for spheroids. First, prolate spheroids engulf fastest when presented to the membrane tip-first, which agrees with previous phagocytic measurements (4) and is a prediction for other types of endocytosis. Sufficiently spherical targets will completely engulf in either orientation, whereas targets that are too eccentric are not engulfed in any orientation. This could well explain the observation in ref. 4, where the membrane sometimes only spreads along the target with no internalization. However, there also exists an intermediate regime where prolate spheroids are engulfed only when one of the highly curved tips is engulfed first. Although, in phagocytosis, this orientation dependence can be explained by adding active processes and a second engulfment stage to the 1D model (28), our 2D model shows that this (for all types of endocytosis) can also be understood simply by receptor diffusion. Second, we predict that oblate spheroids engulf fastest when presented flat to the membrane. This finding is perhaps surprising, because it is then the low curvature region that is engulfed first (and last) rather than the highly curved tips as for prolate spheroids. Hence, oblate and prolate spheroids have completely different engulfment behaviors even for superficially similar shapes (such as the two standing spheroids in Fig. 5*C*). Further experiments with spheroidal targets should help clarify this prediction. Third, the 2D model predicts that the absolute quickest engulfment for a spheroid (with fixed surface area) occurs for a slightly oblate spheroid that lies flat with respect to the membrane. The last finding is confirmed by Sharma et al. (35) for macrophages.

## Discussion

In analogy with the freezing of water, as described by the heat equation with a moving boundary, we used a minimal model of receptor diffusion and capture to study shape and orientation dependence during endocytosis. The process of engulfment is then equivalent to a first-order phase transition characterized by a transition from mobile, unbound receptors to immobile, ligand-bound receptors. By including signaling and a role of actin, we were then able to adapt our model to phagocytosis and examine the different role shape plays for receptor-mediated endocytosis and phagocytosis.

We both extended a previous 1D model to include nonspherical targets and developed a more realistic 2D version, allowing lower-symmetry particles to be studied. The 1D model shows that a whole host of engulfment behaviors emerge from different particle types. This supports the idea that the local curvature at the cup edge is a crucial factor in determining whether (and how quickly) targets are engulfed (4). Highly curved regions may prevent complete engulfment and so are likely to be factors that, at least partially, drive the evolution of pathogen shapes. We found that, for some shapes (such as capped cylinders), the more spherical the target the quicker the engulfment, whereas, for other shapes (such as spheroids), nonspherical particles can be engulfed faster. Thus, long, narrow filamentous shapes could, in some situations, be more beneficial for pathogens that try to avoid ingestion. Further, our hourglass-shaped particles also explain the observation that phagocytosis of budded yeast can stall at half engulfment (1, 30).

The 2D model allowed, for the first time to our knowledge, the orientation dependence of engulfment to be probed. We found that, for nonspherical targets, engulfment can sometimes proceed at completely different rates in different directions, with higher-curved directions engulfing much slower. This should be directly observable with, for example, scanning electron microscopy. Further, we discovered that this effect can explain the experimental observation that prolate spheroids phagocytose fastest when presented to the cell tip-first (4). In addition, the 2D model predicts both that oblate spheroids engulf fastest when presented flat to the cell and that slightly oblate spheroids can engulf even faster than the equivalent sphere. These perhaps counterintuitive results could readily be tested in the near future.

Although phagocytosis is distinct to other forms of receptor-mediated endocytosis, the role of receptors and their motion toward the cup is equally important (28). By adding signaling and a simple role for actin, our models predict that all types of endocytosis, including phagocytosis, have similar shape dependence. In addition, we found that phagocytosis not only fulfils its traditional role of allowing larger particles to be ingested but also allows more highly curved targets to be engulfed.

It is important to point out the limitations of our models. We have focused only on passive diffusion of receptors, neglecting more active receptor motion (15, 16). Further, in our phagocytosis model, we have only included one signaling molecule and have added actin in a highly simplified manner without any attempt to model the underlying structure of the cytoskeleton. Although these simplifications will not capture many aspects of endocytosis, it is likely that the interplay between target (and hence membrane) curvature and the motion of receptors is correctly represented even in our minimal models. Our phagocytosis model also assumes that signaling reduces the membrane bending modulus, effectively making bending around highly curved objects easier. An alternative model representation could have focused on forces instead of receptor dynamics so that, for example, signaling and actin polymerization would lead to a membrane pushing force. Taking spatial derivatives of the cup energy would translate our energy-based model to a force-based model (6). Finally, we have assumed that targets do not rotate during engulfment. It would be interesting to also include this in our models, which would then allow further questions, such as the stability of different orientations, to be addressed.

Shape dependence is a fascinating aspect of endocytosis, which is probably even more important for successful engulfment than target size (4). Understanding the role of shape may also have direct applications to, for example, how some pathogens invade cells (36), how other pathogens avoid uptake, the optimal design for drug delivery carriers (37, 38), and even cell behavior (39). We have shown that there is great variety in the uptake of different target shapes, which is also heavily orientation-dependent. Future work, using differently coated particles, will need to address how geometric, mechanical, and chemical signals are integrated to explain the wide range of observed engulfment outcomes.

## Materials and Methods

### The 1D Numerical Simulations.

For nonspherical particles, even the 1D model must be solved numerically. We used a grid with spacing *L* = 50 μm. Eq. **1** was solved using the Euler method with time step **2**. We checked convergence by showing that smaller values of *SI Appendix* for full details.

### The 2D Numerical Simulations.

Our 2D numerical solutions were calculated on a 2D grid with *r*, a fact that had to be taken into account. We solved Eqs. **6** and **7** with the usual Euler method using *t* and angle *θ*, the normal to the cup surface was estimated from *a* away from the discrete *θ*-lattice points (*SI Appendix*.

## Acknowledgments

D.M.R. and R.G.E. were supported by BBSRC Grant BB/I019987/1. D.M.R. was also supported by the Wellcome Trust Institutional Strategic Support Award (WT105618MA). R.G.E. also acknowledges funding from ERC Starting Grant 280492-PPHPI.

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. Email: david.richards{at}exeter.ac.uk.

Author contributions: D.M.R. and R.G.E. designed research; D.M.R. performed research; and D.M.R. and R.G.E. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1521974113/-/DCSupplemental.

Freely available online through the PNAS open access option.

## References

- ↵
- ↵.
- Rittig MG, et al.

*Borrelia burgdorferi*. Infect Immun 60(10):4205–4212. - ↵
- ↵.
- Champion JA,
- Mitragotri S

- ↵
- ↵
- ↵
- ↵.
- Diakonova M,
- Bokoch G,
- Swanson JA

- ↵.
- Dart AE,
- Tollis S,
- Bright MD,
- Frankel G,
- Endres RG

- ↵
- ↵.
- Griffin FM Jr,
- Griffin JA,
- Leider JE,
- Silverstein SC

- ↵.
- Griffin FM Jr,
- Griffin JA,
- Silverstein SC

- ↵.
- Šamaj J, et al.

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Herant M,
- Heinrich V,
- Dembo M

- ↵
- ↵
- ↵
- ↵
- ↵.
- Gao H,
- Shi W,
- Freund LB

- ↵
- ↵
- ↵
- ↵.
- Dieckmann R, et al.

- ↵
- ↵
- ↵.
- McMahon HT,
- Boucrot E

- ↵.
- Hoppe AD,
- Swanson JA

- ↵
- ↵.
- Hagedorn M,
- Rohde KH,
- Russell DG,
- Soldati T

- ↵
- ↵
- ↵

## Citation Manager Formats

## Article Classifications

- Physical Sciences
- Applied Physical Sciences

- Biological Sciences
- Biophysics and Computational Biology