Raj Kumar
Sadhu
*a,
Sarah R.
Barger
b,
Samo
Penič
c,
Aleš
Iglič
c,
Mira
Krendel
d,
Nils C.
Gauthier
e and
Nir S.
Gov
*a
aDepartment of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 7610001, Israel. E-mail: raj-kumar.sadhu@curie.fr; nir.gov@weizmann.ac.il
bMolecular, Cellular, Developmental Biology, Yale University, New Haven, USA
cLaboratory of Physics, Faculty of Electrical Engineering, University of Ljubljana, Ljubljana, Slovenia
dDepartment of Cell and Developmental Biology, State University of New York Upstate Medical University, Syracuse, USA
eIFOM, FIRC Institute of Molecular Oncology, Milan, Italy
First published on 18th November 2022
Phagocytosis is the process of engulfment and internalization of comparatively large particles by cells, and plays a central role in the functioning of our immune system. We study the process of phagocytosis by considering a simplified coarse grained model of a three-dimensional vesicle, having a uniform adhesion interaction with a rigid particle, and containing curved membrane-bound protein complexes or curved membrane nano-domains, which in turn recruit active cytoskeletal forces. Complete engulfment is achieved when the bending energy cost of the vesicle is balanced by the gain in the adhesion energy. The presence of curved (convex) proteins reduces the bending energy cost by self-organizing with a higher density at the highly curved leading edge of the engulfing membrane, which forms the circular rim of the phagocytic cup that wraps around the particle. This allows the engulfment to occur at much smaller adhesion strength. When the curved membrane-bound protein complexes locally recruit actin polymerization machinery, which leads to outward forces being exerted on the membrane, we found that engulfment is achieved more quickly and at a lower protein density. We consider spherical and non-spherical particles and found that non-spherical particles are more difficult to engulf in comparison to the spherical particles of the same surface area. For non-spherical particles, the engulfment time crucially depends on the initial orientation of the particles with respect to the vesicle. Our model offers a mechanism for the spontaneous self-organization of the actin cytoskeleton at the phagocytic cup, in good agreement with recent high-resolution experimental observations.
Phagocytosis (and endocytosis) in general requires receptors on the cell membrane that bind to the surface of the target particles.10 Binding-driven engulfment, which does not involve the forces of the cytoskeleton, was modelled in ref. 10–15 with purely diffusive dynamics of the receptors10,12 and with an active drift.11,14 The engulfment of charged particles was previously modelled and shown to depend on the lateral segregation of mobile charged membrane components (lipids and proteins), leading to a discontinuous wrapping transition.16 The effect of signalling molecules was also considered in ref. 11, where using drift and diffusion the model successfully recovered the two stages of the process: an initial slow engulfment followed by a fast engulfment of the second half of the spherical object, compared to the experiments of the neutrophil phagocytosis of polystyrene beads.11 Phagocytosis is, however, an active process which is known to involve the forces of actin polymerization that push the engulfing membrane forward.17,18 In subsequent models, a detailed description of active forces18–21 and also the bending and elastic energies of the membrane and cortex are considered.22–26 The cytoskeletal components of phagocytosis were investigated during the engulfment of artificially produced particles14,27–29 or during bacterial engulfment by the immune cells.30–32 In ref. 14, it was found that the engulfment process either stalls before the half engulfment of a spherical object or leads to complete engulfment. In ref. 27 and 33, the experiment was performed using micro-pipettes to hold both the immune cell and the target particle, allowing the measurement of the cortical tension. The recent experiments provide high-resolution details of the organization of the actin cytoskeleton at the rim of the phagocytic cup,34–36 challenging the theoretical modelling of this process.
In reality, however, the shape of the engulfed particles, such as bacteria or viruses, may be highly non-spherical, which motivated several experiments with the artificial particles of various geometries.28,37–39 In ref. 28, it was found that the engulfment time for non-spherical shapes can be five times larger than that for spherical particles.37 In ref. 38, it was observed that an oblate spheroid was engulfed more easily compared to spheres, while spheres were easier to engulf compared to prolate spheroids, rod-shaped particles, or needles.40 Note that the surface area was not kept constant for the particles with different geometries. In theoretical studies, several non-spherical shapes were modelled, such as ellipsoids, rod-like particles, and capped cylinders;41–43 however, in these studies, the effect of active forces was not considered.
From these studies, we know that the phagocytosis process involves highly complex dynamics of cytoskeleton rearrangements, membrane shape deformations and protein aggregations.7 At present, there is no complete theoretical understanding of the dynamics of the self-organization of the membrane and the actin cytoskeleton, including the active forces it exerts, during the engulfment process. Here, we address this problem using a coarse-grained theoretical model.44–46 In this model, a three-dimensional vesicle that contains curved membrane complexes,47–49 which recruit the protrusive forces of the cytoskeleton,50 was shown to describe spreading over a flat substrate.46 The spreading of the vesicle in this model is driven by the self-organization of the curved proteins in the form of a leading edge aggregate, stabilized by the protrusive forces that describe the actin polymerization pressure. We therefore propose to use this model to describe the dynamics of the spontaneous recruitment of the actin cytoskeleton and the membrane spreading during phagocytosis. This modelling approach can be further motivated by the observations of curved membrane proteins that are also related to the actin recruitment at the leading edge of the cell membrane protrusions.51–54
Using this model, we aim to expose the role of curved membrane proteins, coupled with adhesion, in self-organizing and regulating the spatio-temporal dynamics of the active cytoskeletal forces, during phagocytosis. This minimal model allows us to explore the basic physical mechanisms during this process, which exist in addition to the more complex biochemical signaling, which plays a role in regulating phagocytosis.6,55
![]() | (1) |
![]() | (2) |
![]() | (3) |
The above equation indicates that when the proteins are distributed inhomogeneously, there will be a net force on the vesicle in a particular direction.46 However, in the present work, we only simulate vesicles that are adhered to particles that are fixed in their location, which effectively links our vesicle to the lab frame thereby restoring momentum conservation (fixing the engulfed particle acts as an infinite momentum reservoir).
Finally, the adhesion energy is given by
![]() | (4) |
We update the vesicle with mainly two moves: (1) vertex movement and (2) bond flip. In a vertex movement, a vertex is randomly chosen and attempted to move by a random length and direction, with a maximum possible distance restricted by 0.15 lmin. In a bond flip move, a single bond is chosen, which is a common side of two neighbouring triangles, and this bond is cut and reestablished between the other two unconnected vertices.45,46,59 While the vertex movement allows the vesicle to change its shape, a bond flip mechanism provides lateral fluidity to the membrane, by allowing vertices (and the vertices that are occupied by protein complexes) to diffuse on the membrane surface. This lateral motion allows the protein complexes to move on the membrane surface and aggregate with other protein complexes to form larger aggregates. The maximum bond length is restricted to lmax = 1.7lmin. The particle is assumed to be static and does not have any independent dynamics. We update the system using the Metropolis algorithm, where any movement that increases the energy of the system by an amount ΔW = ΔWb + ΔWd + ΔWF + ΔWA occurs with the rate exp(−ΔW/kBT); otherwise, it occurs with rate unity. Note that our model only captures the wrapping and engulfment events, while the other upcoming processes, such as the membrane fission,60 and the movement of the engulfed particle inside the vesicle are not captured in the present model.
The forces exerted at the locations of the curved proteins are “active” since they give an effective energy (work) term that is unbounded from below (3). This property defines them as out-of-equilibrium, as they can not arise in an energy conserving system. These active forces represent the pushing forces exerted on the membrane by polymerizing actin filaments,61–64 which we do not model explicitly to reduce the complexity of the model.65
Unless specified, we use a vesicle of the total number of vertices, N = 3127 (radius ∼ 20 lmin), where lmin is the unit of length in our model and defines the minimum bond length. The bending rigidity is κ = 20 kBT, the protein–protein attraction strength is w = 1 kBT, and ρ = Nc/N is the protein density, with Nc vertices occupied by curved (convex) membrane proteins having spontaneous curvature: C0 = 1.0 lmin−1. Note that we do not conserve the vesicle volume, although this could be maintained using an osmotic pressure term as the volume constraint.45 The membrane area is very well conserved (ΔA/A < 1%).
Note that our MC simulations do not describe the correct dynamics of the fluid motion around the membrane and the engulfed object. However, it gives the energetically most favorable trajectory that controls the dynamics, without the correct time-scales, which depend on hydrodynamic dissipation. In addition, the actin filaments that exert the pushing forces at the leading edge of the phagocytic cup also transmit traction forces to the engulfed particle, which can therefore move and rotate it. These traction forces are not explicitly described, as we maintain the fixed position of the engulfed particle, similar to some experimental set-ups.66
In this numerical study, we chose a fixed vesicle and spherical particle size, for which we have a reasonable balance between accuracy and simulation time. The vesicle radius was chosen to be Rvesicle ∼ 20 lmin, while the particle radius was taken to be R = 10 lmin (smaller and larger particles are explored in the ESI† Section S3 and Fig. S1). The mesh size for these dimensions allows us to describe the membrane shapes during the engulfment process with good accuracy, even along the sharp (highly curved) leading edge. The leading edge curvature is determined by the spontaneous curvature of the curved proteins, which we kept from our previous work on vesicle spreading.46 We quantified the accuracy by comparing the calculated mean curvature of the membrane with the discrete version and found that, even for spherical vesicles with a radius of ∼1 lmin, the error is only ∼10% (Fig. S21, ESI†). Clearly, membrane features that are on smaller lengthscales than lmin are not captured by this calculation.
Note that our continuum description of the membrane is valid for lengthscales that are larger than that of the membrane width or the single-protein size. Our simulations are therefore valid for particles that are larger than this lengthscale. Furthermore, we use the term “curved proteins” to describe any membrane-bound or membrane embedded protein complexes, or larger nano-domains, that have the properties of the spontaneous curvature and ability to recruit the nucleation of actin polymerization. Let us emphasize that our continuum model of the membrane is not a microscopic model, where each node of the mesh represents a single protein or a lipid.
The length scale lmin that we use in our simulation does not necessarily correspond to a particular length in reality. By assigning a different real length to the basic length-scale of the simulation (lmin), we do not change the resulting dynamics for the following reason: the dynamics of the engulfment process depends on the energy gain/loss per additional adhered area increment (vertex). The adhesion energy per adhered vertex rescales as 1/lmin2, as does the bending energy per vertex (because the total bending energy of a specific shape is scale-invariant, provided we also rescale the spontaneous curvature (C0) accordingly). The active work is also independent of the absolute scale of lmin, as we define the force exerted by each vertex to be in units of kBT/lmin, and therefore, its contribution to the total energy increment per unit area, also scales as 1/lmin2. The absolute scaling of lmin does not therefore change the qualitative nature of our simulation results: both the adhesion energy and the bending energy per unit area, scale as 1/lmin2.
There are important length scales in this system which we wish to discuss. The first is the lengthscale given by the ratio between the adhesion and bending energy of the adhered membrane. This ratio determines the minimal radius of the spherical particle: , below which the membrane does not lower its energy when adhering and therefore adhesion vanishes completely.67 Above the critical radius, the engulfment is inhibited by the interplay between the adhesion energy of the adhered membrane and the bending energy of the non-adhered membrane (see the ESI†, Fig. S2). The bending energy in this regime depends on the ratio of the vesicle to the particle size,67 such that large particles are not engulfed due to the finite size of the vesicle. We demonstrate this using our simulations for the protein-free vesicle in the ESI† (Fig. S4). This analysis and the dynamics of the engulfment in this regime are modified when the vesicle contains curved proteins, as well as active forces, and is the main focus of this work. We leave the detailed exploration of the effects of some of the parameters of our model for future studies, such as C0.
We chose the surface area of the vesicle to be much larger (at least four times larger) than the particles, so that the engulfment is not inhibited by the membrane surface area limitation. During phagocytosis, the cell is able to increase the local membrane area,6 to allow for the phagocytosis of larger particles. Note that our vesicle does not represent the entire shape of the cell and does not include the nucleus and other constraints. Furthermore, we do not explicitly take into account the lateral68 and surface69 membrane tension, but our calculation conserves the membrane area (total area fluctuations of <1%). In addition, we do not fix the vesicle volume, corresponding to the vanishing osmotic pressure gradient across the membrane, as maintained in cells.70 This allows our vesicle to undergo both large volume and shape changes during the engulfment process. Introducing a high osmotic pressure inside our simulated vesicle would act to conserve the volume,45 and in this way hinder the engulfment process by inhibiting large vesicle shape deformations. The effects of osmotic pressure on the phagocytosis process can be further investigated in the future.
![]() | ||
Fig. 2 Sphere engulfment by a vesicle with passive proteins. (a) Phase diagram in the Ead–ρ plane, where Ead is the adhesion strength and ρ is the density of curved membrane proteins. The background color shows the adhered area fraction. The red dashed line is approximately separating the fully engulfed state and the partially engulfed state. The snapshots are shown for Ead = 0.60, 0.80, 1.0, and 1.20 (in units of kBT) and ρ = 1.6%, 4.8%, 9.6%, 12.8%, 16%, and 19.2%. (b) Comparison of the adhesion energy gain (dashed horizontal line) and the bending energy cost (ΔWb, points) per adhered node (both having units of kBT) as a function of the engulfment fraction, for a fixed adhesion strength Ead = 1.0 kBT and various values of ρ. Blue symbols denote ρ = 0.8%, green circles denote ρ = 2.4% and red triangles denote ρ = 4.0%. The origin of the large error bars is due to the large fluctuations in the bending energy at the leading edge of the membrane and having a relatively small ensemble of simulations for each case. (c) Mean cluster size (the mean number of proteins in a cluster) as a function of the engulfed area fraction for both partial and complete engulfment cases. The color codes are same as shown in Fig. 2b. The inset shows the increase in the density of curved proteins at the rim of the phagocytic cup during the engulfment process for ρ = 4.0%. (d) Examples of the cross-sections of the vesicle membrane around the particle, at the end of the engulfment process, for different protein densities. We use Ead = 1.0 kBT here. the other parameters are the total number of vertices (nodes) in the vesicle is N = 3127, the radius of the spherical particle is R = 10 lmin, the bending rigidity is κ = 20 kBT, the protein–protein attraction strength is w = 1 kBT, and the spontaneous curvature of curved (convex) membrane proteins is C0 = 1.0 lmin−1. |
Next, we consider vesicles with passive (F = 0, where F is the active protrusive force exerted by the proteins) curved proteins. In Fig. 2a, we show the phase diagram in the ρ–Ead plane (where Ead is the adhesion energy per adhered node and ρ is the areal density of the proteins on the vesicle membrane). For small Ead, below ∼0.60 kBT, we only obtain partial engulfment, with the engulfment area increasing with ρ but not reaching the full engulfment. For intermediate Ead (0.80–1.0 kBT), there is an incomplete engulfment for small ρ (Movie S3, ESI†), but the engulfment is complete for larger ρ (Movie S4, ESI†). This is the interesting regime, where engulfment is facilitated by the presence of the curved proteins.
Note that we calculate the adhered fraction of the spherical particle (Aad/Amax), using the maximum adhered area for Amax. This may be larger than the true surface area of the particle, since we defined a finite width of the adhesive interactions, and within this width the adhered membrane may fluctuate, therefore containing a larger membrane area than the area of the rigid bead surface (see for example the cross-sections in Fig. 2d, and for more details see the ESI,† Section S2).
Increasing ρ reduces the bending energy cost of the engulfment by the spontaneous aggregation of the curved proteins at the sharp leading edge of the engulfing membrane and allows the vesicle to adhere more and engulf the particle. This is the same mechanism that allows the passive curved proteins to drive vesicle spreading on flat adhesive substrates.46
The engulfment process can be understood by considering the competition between the adhesion and the bending energy cost of the vesicle as the engulfment proceeds. In Fig. 2b, we compare the adhesion energy gain (Ead) and the bending energy cost (ΔWb), per additional adhered membrane node, as a function of the adhered fraction of the spherical particle. The bending energy cost increases with the adhered fraction, and for small ρ, the engulfment process stops where the bending energy cost becomes higher than the adhesion energy gain per node. In Fig. 2b, we show two cases of partial engulfment, for ρ = 0.8% and 2.4% (blue and green points), where the bending energy cost becomes larger than the adhesion energy gain, and the system reaches steady-state at Aad/Amax values of ∼0.2 and 0.35, respectively. At these low densities, the curved proteins fail to form large aggregates at the leading-edge of the engulfing membrane, as shown by the mean cluster size (Fig. 2c), and therefore do not lower the bending energy cost sufficiently (see the insets of Fig. 2c for the adhered fractions that are smaller than 0.5).
The complete engulfment is achieved for higher densities, where the adhesion energy gain is larger than the bending energy cost. For ρ = 4.0%, the bending energy cost is always smaller than the adhesion energy gain, and the system proceeds to complete engulfment (red points in Fig. 2b). This is driven by the formation of a large protein cluster along the leading-edge of the engulfing membrane (Fig. 2c and insets for adhered fractions that are larger than 0.5). However, the aggregation at the leading edge is transient, since as the radius of the hole gets too small, the mean curvature of the rim of the hole decreases and is less favorable for the curved proteins (ESI† Section S5 and Fig. S5), which leads to the dispersal of the curved proteins on the membrane (Fig. 2d, ρ = 4.8%). This is the origin for the appearance of a maximum value of the mean cluster size during complete engulfment (Fig. 2c).
Note that in addition to the adhesion energy gain, there is also the gain in the protein–protein binding energy, which further helps drive the membrane over the bending energy barrier, but this is a small contribution due to the small size of the ring-like clusters of proteins (Fig. 2c). Following the completion of the engulfment, this ring-like cluster of curved proteins disperses spontaneously. Since the topology of the membrane is fixed, we do not allow for membrane fusion, and there remains a tiny open membrane neck that connects the engulfed sphere with the space outside the vesicle (Fig. 2d).
For very large ρ (Fig. 2a), however, the wrapping is again incomplete. The large number of curved proteins stabilizes the cluster at the membrane rim and prevent the hole from closing up (Movie S5, ESI†). In Fig. 2d, we show an example of the cross-sectional view of the vesicle at the end of the engulfment process. For high densities (ρ = 9.6% or larger), the protein cluster remains stable at the rim of the phagocytic cup. This cluster maintains the curvature of the rim to fit better the spontaneous curvature of the proteins, and this prevents the rim from shrinking, resulting in incomplete engulfment (ESI† Section S6 and Fig. S6). Note that at large densities the curved proteins form pearled clusters, which are not directly related to the rim cluster at the edge of the phagocytic cup.
For large Ead(≥1.2 kBT), the adhesion interaction between the vesicle and the particle is strong enough to drive complete engulfment even in the absence of curved proteins. The red dashed line in Fig. 2a encloses the region within which we find complete engulfment. The shape of this complete engulfment transition line can be qualitatively reproduced by a simplified analytical model (ESI† Section S7 and Fig. S8, S9).
Finally, we also calculated the engulfment time as a function of ρ (see the ESI† Section S8 and Fig. S10), where we found that it is decreasing with increasing ρ. The engulfment time diverges for ρ ∼ 2.25%, below which there is no engulfment (for Ead = 1.0 kBT).
For the passive case, we analyze the engulfment process by plotting the adhesion energy gain and the bending energy cost per adhered node, as a function of Aad/Amax. In Fig. 3b, we compare both the case of the complete and incomplete engulfment, as the force is increased. As in the passive case, incomplete engulfment occurs when the bending energy cost per node becomes larger than Ead (blue boxes in Fig. 3b). For a higher force (F = 0.40 kBT/lmin), the active force enhances the aggregation of the proteins around the rim and thereby lowers the bending energy cost, resulting in complete engulfment (green circles in Fig. 3b).
Even higher forces (F = 2 kBT/lmin) produce large shape fluctuations of the vesicle, and the bending energy cost increases. However, a complete engulfment is achieved nevertheless (red triangles in Fig. 3b and top insets in Fig. 3c). In this case, the engulfment is enabled by the active forces that directly push the membrane over the spherical surface, and complete engulfment is achieved even when the bending energy cost per node goes above Ead. In order to extract the contribution due to active forces in this case, we define a region close to the particle (phagocytic cup or rim) and assume that only the forces from these proteins are helping the vesicle to engulf the particle. We then calculate the magnitude of active forces in the tangential direction to the surface of the particle (Fig. 3d). We then divide this total force Ftotal by the length of the leading rim of the membrane to extract the effective active work done per adhered node (see the ESI† Section S9 and Fig. S11 for details). We add this active work to the adhesion energy gain per adhered node (the red dashed line in Fig. 3b), which is shown to be sufficient to offset the bending energy cost per node. The magnitude of the total active force (without scaling by the rim length) as a function of the adhered fraction is shown in the ESI,† Fig. S11(c). It is clear from Fig. S11(c) (ESI†) that the total active force is very small during the early stages when the curved protein complexes are still dispersed over the membrane surface, so that their active forces are pointing in random directions, before they form the organized aggregate at the leading edge of the phagocytic cup.
Note that in the active case the complete engulfment can proceed at much lower densities of proteins, compared to the passive case, such that there are not necessarily enough of them to complete a ring-like cluster at the leading edge. An arc-like cluster of active proteins, which forms at an earlier stage of the engulfment process (red triangles in Fig. 3c), can be sufficient to drive a flat protrusion that spreads over the spherical particle to complete the engulfment (Fig. 3c and d). This is further shown in Fig. 3e for larger protein densities (Movie S9 and S10, ESI†). The spherical particle is transparent to ensure the visibility of the protein clusters along the rim of the phagocytic cup. Due to the higher density, the rim clusters form at a lower adhered fraction (compared to Fig. 3c), but they are still highly fragmented and do not merge to form a complete circular ring until very late stages of the engulfment (Fig. 3e).
Many of the qualitative spatio-temporal features that we found for the curved active proteins in our simulations match very well with the recent high-resolution experimental images of the actin organization along the leading edge of the engulfing membrane protrusion during phagocytosis (Fig. 4 and Movie S11, S12, ESI†).35 In both the experiments and the simulations, we found that the actin ring50 is initially discontinuous, and highly fragmented, in the form of dispersed “teeth” and “arcs” along the leading edge. The actin aggregate becomes more cohesive in the final stages of the engulfment, before it disperses after the engulfment is complete (also see the ESI† Section S10 and Movie S13). Note that in our model the curved membrane complexes are also the sites of the recruitment of actin polymerization, which therefore ends up exerting a protrusive force on the membrane at these sites. We therefore identify the membrane accumulation of the curved proteins in our simulations with the accumulation of cortical actin in the experimental observations.
When comparing to the experiments, we have to remember the limitations of our numerical simulation, which can only describe membrane features that are on lengthscales larger than lmin (for a particle of radius 3.5 μm, lmin = 350 nm). However, the main features of the actin organization that we highlight in Fig. 4a, namely the fragmented ring and its continuous closure, are on larger lengthscales and therefore reliably emerge in our simulations. We therefore found this choice of mesh resolution to be sufficient in order to explain the outstanding features of the membrane and actin organization during engulfment, as observed in the experiments, while maintaining a reasonable computation time. Future studies using a finer mesh may allow smaller scale features of this process to be studied.
For the protein-free and passive proteins systems, driven by adhesion alone, we found in our simulations the typical “two-stage” behavior: an initial slow growth up to an adhered area fraction of ∼0.5, and then a very fast stage to complete engulfment. A similar behavior is also observed for the active system in the regime of low protein densities. This dynamics was observed experimentally, and for the passive engulfment system it was also motivated by theoretical modelling.11,16
However, for the active system at high protein densities, where the rim clusters form at an early stage of the engulfment process (Fig. 3e), we found a much more uniform dynamics with steady growth of the adhered fraction (Fig. 5a). This agrees qualitatively with the observations in,66 as well as,34,35 which indeed correspond to this regime (Fig. 3e and 4).
In Fig. 5a, the engulfment dynamics are shown as a function of the normalized time, so the actual engulfment duration is not presented. The dependence of the mean engulfment time on the active force is found to be non-monotonic (Fig. 5b). Below a critical force, in the regime of low Ead and low ρ, where the passive system does not fully engulf (Fig. 3a), there is no complete engulfment (the green shaded area). Above this regime, the increasing force induces complete engulfment with faster and more robust engulfment as F increases. The engulfment time distribution in this regime is relatively narrow (Fig. 5c).
However, increasing F beyond some values, which are set by the relationship of the adhesion energy and the active work, leads to increasing both the mean value and the variability of the engulfment time. This occurs due to active forces stretching the membrane away from the adhered particle (top insets in Fig. 3c), which gives rise to long “waiting times” during which the engulfment does not progress. Following this “waiting time”, once the protein aggregate forms at the membrane rim (Fig. 3c), the engulfment proceeds rapidly as shown in Fig. 5a (yellow circles). This variability is manifested in the very broad distribution of engulfment times (Fig. 5d).
At even larger forces, the stretching of the membrane sideways by the active forces pulls the membrane away from the adhesive particle and prevents engulfment (the blue shaded area in Fig. 5b, and ESI† Section S11, Fig. S13(a), Movie S14).
In this regime of large F, if we further increase ρ, we reach the regime where a free vesicle will form a flat pancake-like shape, where all the proteins are aggregated at the rim.45 In this regime, the vesicle can engulf the particle if this process is faster than the timescale of the formation of the pancake-like shape (ESI† Section S11, Fig. S13(b) and Movie S15). For very large forces, we found again a non-engulfed regime (similar to the blue shaded area in Fig. 5b and ESI† Fig. S13(c), Movie S16).
We first consider spheroids. We keep the surface area constant, i.e. same as for the sphere of radius 10 lmin. We use prolate shapes (Rx = Ry < Rz) and oblate shapes (Rx = Ry > Rz), such that the aspect ratios are Rx/Rz < 1 for prolate and Rx/Rz > 1 for oblate.
In Fig. 6, we show the mean engulfment time for oblate and prolate spheroids, for both the cases when the vesicle is initially in contact from the top (the poles of the shape) or from the side (along the equator), as shown in the insets. For both passive and active curved proteins, we found that the engulfment times increase as the shape deviates from a sphere. This is in agreement with previous experimental observations28 and can be attributed to the higher local curvature on the non-spherical shapes, which increase the bending energy barrier that opposes the progression of the leading edge membrane during their engulfment. Furthermore, we found that the engulfment time increases more rapidly for the oblate shapes, compared to the prolate shapes, for both passive and active, as was observed in experiments.28
For the particles with a small aspect ratio (Rx/Rz close to unity), the engulfment time of all the cases is comparable. However, for the particles with large aspect ratios, there are significant differences in the engulfment times, between the passive and active systems. For the passive curved proteins (blue triangles and red boxes), we found that the engulfment times increase sharply with the deviation from the spherical shape. Furthermore, we found that initialising from the side (red boxes) results in relatively shorter engulfment times, in particular for the oblate shapes. For an oblate particle with a large aspect ratio (for Rx = Ry > 11 lmin), the particle is engulfed from side, while it is not engulfed from the top (Fig. 6). For prolate shapes with large aspect ratios, however, the engulfment times seem comparable for both initial conditions. The differences between the engulfment dynamics on the oblate and prolate shapes arise from the different local curvatures that the edge membrane encounters during engulfment: when adhering to the side (rim) of the oblate shape, the leading edge moves mostly on the flat surfaces, which allow the balance of adhesion versus bending to be favorable for the smooth progression of the forward motion. When spreading from the top of the prolate shape, the membrane edge encounters the highly bent surface all along the leading edge, and therefore faces a higher bending energy barrier for the forward movement of the membrane.72 When spreading from the side of the prolate shape, it is similar to the spreading of the vesicle on a cylinder, where one direction has low curvature, and therefore faster engulfment (see these dynamics in terms of the energy of adhesion and bending in Fig. S15, ESI†).
By comparison, for active curved proteins, the engulfment time is lower than that for the passive case, even with a lower protein density, and the active curved proteins enable complete engulfment over a larger range of aspect ratios. We also note that for the active curved proteins it is the prolate shape that exhibits a dependence of the engulfment time on the initial conditions, where starting from the side (black circles) leads to faster engulfment. For the oblate shapes, however, the engulfment time is largely independent of the initial conditions, opposite to the behavior of the passive curved proteins.
For a non-spherical particle with high aspect ratios, due to the highly curved poles (in the prolate) or equator (oblate) regions, the engulfment process is slower and the variability in engulfment times increases. The slower engulfment allows the vesicle to reorient on the surface of the particle, driven by both passive energy minimization and active work. We show the details of this reorientation process in the ESI† Section S12 (Movies S17–S22).
When the vesicle is initially at the top of a prolate shape with a high aspect ratio, the vesicle takes a longer time to form the initial adhesion since the particle's surface facing the vesicle has high curvature in all directions. Following this slow initial adhesion time, the vesicle reorients itself always to the side and engulfs the particle, along the easier direction where the long axis of the prolate shape presents a direction of the low curvature (Fig. S15e and f, ESI†). For passive curved proteins, this reorientation process is faster and thereby the engulfment time when started from the top or from the side are comparable (blue triangles and red boxes in Fig. 6 for Rx/Rz < 1). For active curved proteins, this reorientation process is slower because of large shape fluctuations of the vesicle, and thus the engulfment time when started from the side is smaller than that from the top (green inverted-triangles and black circles in Fig. 6 for Rx/Rz < 1).
For the oblate shapes, we found a similar reorientation process when the vesicle starts the adhesion from the highly curved sides. In this geometry, the passive vesicle can engulf from side without any reorientation, which is very fast, or it reorients to the top of the oblate shape and ends up with no engulfment. Therefore, when we calculate the engulfment time, we only average over those realizations where the particle was engulfed without any reorientation (red boxes for Rx > 11). Otherwise, the mean engulfment time of oblate shapes by passive vesicles diverges for both the top and side orientations around an aspect ratio of ∼1.25 (Fig. 6). This value can be explained by the bending energy barrier along the oblate rim and has the local radius of the curvature similar to the smallest spherical particle that can still be engulfed (ESI† Section S13 and Fig. S17).
For the active case, the vesicle reorients in every realization and engulfs from the top. This reorientation process is very fast and the engulfment time is similar for starting either from the side or from the top (green inverted-triangles and black circles in Fig. 6 for Rx/Rz > 1).
We further explore more extreme non-spherical shapes, such as elongated sphero-cylinders and dumb-bell shapes in the ESI† (Section S14 and Movies S23–S26). The elongated particles resemble the shapes of different bacteria or fungi such as E. coli, Bacillus subtilis and budding yeast.12,50,73
We found that the curved proteins spontaneously aggregate around the highly curved rim of the phagocytic cup and enable the engulfment to proceed efficiently by reducing the bending energy costs. The active forces in our model, which represent the protrusive forces due to the recruited actin polymerization, further facilitate the engulfment, due to several effects that we have identified: the protrusive force at the membrane leading edge promotes the aggregation of the curved proteins (reducing the bending energy cost) and directly pushing the membrane over the engulfed particle. However, our model predicts that if the active forces are too large, they induce membrane ruffles that are inefficient for smooth engulfment and may hinder the phagocytosis process.
In addition, the active forces can stretch the membrane to wrap around sharp corners, which enables the engulfment of non-spherical particles. Compared to the passive curved proteins, the curved active proteins drive more efficient engulfment at lower adhesion energy, at lower curved protein density, and of highly non-spherical particles. For oblate and prolate shapes, we recover the experimentally observed increase in engulfment time with the aspect ratio. We expose the mechanisms that allow the active forces to more robustly overcome the local curvature barriers during the engulfment of such non-speherical objects. These results of our model highlight the central role that the recruitment of actin plays during phagocytosis, which is usually abolished by the inhibition of actin polymerization (see for example74).
Our model predicts that the actin aggregation at the rim of the phagocytic cup is typically fragmented and does not form a complete circular ring. Nevertheless, it is still highly effective in driving the engulfment. Furthermore, the actin aggregate forms a more cohesive ring in the final stages of the engulfment, before it spontaneously disperses as the engulfment is complete. These features of our model are verified and observed in recent high-resolution imaging of the actin during phagocytosis35 (Fig. 4).
The model we presented currently lacks the contractile forces that arise at the phagocytic cup due to myosin-II activity.35 The contractile forces play an important role during the engulfment of soft objects, such as cells,75–77 but they are dispensable for the phagocytosis of more rigid particles.35,74 The contribution of myosin-II-mediated contractility may be more or less pronounced depending on the particle size and stiffness and, in addition, could influence mechanosensitive signaling within phagocytes rather than directly drive phagocytic cup extension and closure. Our results should therefore be important even without this feature, which we plan to add to our model in the future, in particular when modelling the engulfment of flexible particles.
To conclude, we demonstrate that curved membrane proteins can play a key role during phagocytosis. Convex-shaped proteins or a curved complex of proteins, which furthermore recruit actin polymerization, spontaneously aggregate at the leading edge of the membrane and drive robust engulfment. This is similar to the role of these proteins that was recently found at the leading edge of the lamellipodia, in experiments53,54 and theory.46 We also show that concave-shaped proteins can further enhance the engulfment process by reducing the bending energy cost of the membrane that adheres to the particle's surface78,79 (ESI† Section S15, Fig. S20 and Movie S27). Our results should motivate future experimental studies aimed at the identification and characterization of the different curved membrane components that participate in phagocytosis.3,52,80
Footnote |
† Electronic supplementary information (ESI) available: High resolution supplementary movies. See DOI: https://doi.org/10.1039/d2sm01152b |
This journal is © The Royal Society of Chemistry 2023 |