Active particles induce large shape deformations in giant lipid vesicles.

Biological cells generate intricate structures by sculpting their membrane from within to actively sense and respond to external stimuli or to explore their environment1-4. Several pathogenic bacteria also provide examples of how localized forces strongly deform cell membranes from inside, leading to the invasion of neighbouring healthy mammalian cells5. Giant unilamellar vesicles have been successfully used as a minimal model system with which to mimic biological cells6-11, but the realization of a minimal system with localized active internal forces that can strongly deform lipid membranes from within and lead to dramatic shape changes remains challenging. Here we present a combined experimental and simulation study that demonstrates how self-propelled particles enclosed in giant unilamellar vesicles can induce a plethora of non-equilibrium shapes and active membrane fluctuations. Using confocal microscopy, in the experiments we explore the membrane response to local forces exerted by self-phoretic Janus microswimmers. To quantify dynamic membrane changes, we perform Langevin dynamics simulations of active Brownian particles enclosed in thin membrane shells modelled by dynamically triangulated surfaces. The most pronounced shape changes are observed at low and moderate particle loadings, with the formation of tether-like protrusions and highly branched, dendritic structures, whereas at high volume fractions globally deformed vesicle shapes are observed. The resulting state diagram predicts the conditions under which local internal forces generate various membrane shapes. A controlled realization of such distorted vesicle morphologies could improve the design of artificial systems such as small-scale soft robots and synthetic cells.

Wrapping anisotropic microgel particles in lipid membranes: Effects of particle shape and membrane rigidity.

Cellular engulfment and uptake of macromolecular assemblies or nanoparticles via endocytosis can be associated to both healthy and disease-related biological processes as well as delivery of drug nanoparticles and potential nanotoxicity of pollutants. Depending on the physical and chemical properties of the system, the adsorbed particles may remain at the membrane surface, become wrapped by the membrane, or translocate across the membrane through an endocytosis-like process. In this paper, we address the question of how the wrapping of colloidal particles by lipid membranes can be controlled by the shape of the particles, the particle-membrane adhesion energy, the membrane phase behavior, and the membrane-bending rigidity. We use a model system composed of soft core-shell microgel particles with spherical and ellipsoidal shapes, together with phospholipid membranes with varying composition. Confocal microscopy data clearly demonstrate how tuning of these basic properties of particles and membranes can be used to direct wrapping and membrane deformation and the organization of the particles at the membrane. The deep-wrapped states are more favorable for ellipsoidal than for spherical microgel particles of similar volume. Theoretical calculations for fixed adhesion strength predict the opposite behavior-wrapping becomes more difficult with increasing aspect ratio. The comparison with the experiments implies that the microgel adhesion strength must increase with increasing particle stretching. Considering the versatility offered by microgels systems to be synthesized with different shapes, functionalizations, and mechanical properties, the present findings further inspire future studies involving nanoparticle-membrane interactions relevant for the design of novel biomaterials and therapeutic applications.

Collective behavior of squirmers in thin films.

Bacteria in biofilms form complex structures and can collectively migrate within mobile aggregates, which is referred to as swarming. This behavior is influenced by a combination of various factors, including morphological characteristics and propulsive forces of swimmers, their volume fraction within a confined environment, and hydrodynamic and steric interactions between them. In our study, we employ the squirmer model for microswimmers and the dissipative particle dynamics method for fluid modeling to investigate the collective motion of swimmers in thin films. The film thickness permits a free orientation of non-spherical squirmers, but constraints them to form a two-layered structure at maximum. Structural and dynamic properties of squirmer suspensions confined within the slit are analyzed for different volume fractions of swimmers, motility types (e.g., pusher, neutral squirmer, puller), and the presence of a rotlet dipolar flow field, which mimics the counter-rotating flow generated by flagellated bacteria. Different states are characterized, including a gas-like phase, swarming, and motility-induced phase separation, as a function of increasing volume fraction. Our study highlights the importance of an anisotropic swimmer shape, hydrodynamic interactions between squirmers, and their interaction with the walls for the emergence of different collective behaviors. Interestingly, the formation of collective structures may not be symmetric with respect to the two walls. Furthermore, the presence of a rotlet dipole significantly mitigates differences in the collective behavior between various swimmer types. These results contribute to a better understanding of the formation of bacterial biofilms and the emergence of collective states in confined active matter.

Motion of microswimmers in cylindrical microchannels.

Biological and artificial microswimmers often have to propel through a variety of environments, ranging from heterogeneous suspending media to strong geometrical confinement. Under confinement, local flow fields generated by microswimmers, and steric and hydrodynamic interactions with their environment determine the locomotion. We propose a squirmer-like model to describe the motion of microswimmers in cylindrical microchannels, where propulsion is generated by a fixed surface slip velocity. The model is studied using an approximate analytical solution for cylindrical swimmer shapes, and by numerical hydrodynamics simulations for spherical and spheroidal shapes. For the numerical simulations, we employ the dissipative particle dynamics method for modelling fluid flow. Both the analytical model and simulations show that the propulsion force increases with increasing confinement. However, the swimming velocity under confinement remains lower than the swimmer speed without confinement for all investigated conditions. In simulations, different swimming modes (i.e. pusher, neutral, puller) are investigated, and found to play a significant role in the generation of propulsion force when a swimmer approaches a dead end of a capillary tube. Propulsion generation in confined systems is local, such that the generated flow field generally vanishes beyond the characteristic size of the swimmer. These results contribute to a better understanding of microswimmer force generation and propulsion under strong confinement, including the motion in porous media and in narrow channels.

Dynamic shapes of floppy vesicles enclosing active Brownian particles with membrane adhesion.

Recent advances in micro- and nano-technologies allow the construction of complex active systems from biological and synthetic materials. An interesting example is active vesicles, which consist of a membrane enclosing self-propelled particles, and exhibit several features resembling biological cells. We investigate numerically the behavior of active vesicles, where the enclosed self-propelled particles can adhere to the membrane. A vesicle is represented by a dynamically triangulated membrane, while the adhesive active particles are modelled as active Brownian particles (ABPs) that interact with the membrane via the Lennard-Jones potential. Phase diagrams of dynamic vesicle shapes as a function of ABP activity and particle volume fraction inside the vesicle are constructed for different strengths of adhesive interactions. At low ABP activity, adhesive interactions dominate over the propulsion forces, such that the vesicle attains near static configurations, with protrusions of membrane-wrapped ABPs having ring-like and sheet-like structures. At moderate particle densities and strong enough activities, active vesicles show dynamic highly-branched tethers filled with string-like arrangements of ABPs, which do not occur in the absence of particle adhesion to the membrane. At large volume fractions of ABPs, vesicles fluctuate for moderate particle activities, and elongate and finally split into two vesicles for large ABP propulsion strengths. We also analyze membrane tension, active fluctuations, and ABP characteristics (e.g., mobility, clustering), and compare them to the case of active vesicles with non-adhesive ABPs. The adhesion of ABPs to the membrane significantly alters the behavior of active vesicles, and provides an additional parameter for controlling their behavior.

Erythrocyte Sedimentation: Collapse of a High-Volume-Fraction Soft-Particle Gel.

The erythrocyte sedimentation rate is one of the oldest medical diagnostic methods whose physical mechanisms remain debatable today. Using both light microscopy and mesoscale cell-level simulations, we show that erythrocytes form a soft-particle gel. Furthermore, the high volume fraction of erythrocytes, their deformability, and weak attraction lead to unusual properties of this gel. A theoretical model for the gravitational collapse is developed, whose predictions are in agreement with detailed macroscopic measurements of the interface velocity.

Correction: Non-equilibrium shapes and dynamics of active vesicles.

Correction for 'Non-equilibrium shapes and dynamics of active vesicles' by Priyanka Iyer et al., Soft Matter, 2022, 18, 6868-6881, https://doi.org/10.1039/D2SM00622G.

Non-equilibrium shapes and dynamics of active vesicles.

Active vesicles, constructed through the confinement of self-propelled particles (SPPs) inside a lipid membrane shell, exhibit a large variety of non-equilibrium shapes, ranging from the formation of local tethers and dendritic conformations, to prolate and bola-like structures. To better understand the behavior of active vesicles, we perform simulations of membranes modelled as dynamically triangulated surfaces enclosing active Brownian particles. A systematic analysis of membrane deformations and SPP clustering, as a function of SPP activity and volume fraction inside the vesicle is carried out. Distributions of membrane local curvature, and the clustering and mobility of SPPs obtained from simulations of active vesicles are analysed. There exists a feedback mechanism between the enhancement of membrane curvature, the formation of clusters of active particles, and local or global changes in vesicle shape. The emergence of active tension due to the activity of SPPs can well be captured by the Young-Laplace equation. Furthermore, a simple numerical method for tether detection is presented and used to determine correlations between the number of tethers, their length, and local curvature. We also provide several geometrical arguments to explain different tether characteristics for various conditions. These results contribute to the future development of steerable active vesicles or soft micro-robots whose behaviour can be controlled and used for potential applications.

Emergent collective behavior of active Brownian particles with visual perception.

Systems comprised of self-steering active Brownian particles are studied via simulations for a minimal cognitive flocking model. The dynamics of the active Brownian particles is extended by an orientational response with limited maneuverability to an instantaneous visual input of the positions of neighbors within a vision cone and a cut-off radius. The system exhibits large-scale self-organized structures, which depend on selected parameter values, and, in particular, the presence of excluded-volume interactions. The emergent structures in two dimensions, such as worms, worm-aggregate coexistence, and hexagonally close-packed structures, are analysed and phase diagrams are constructed. The analysis of the particle's mean-square displacement shows ABP-like dynamics for dilute systems and the worm phase. In the limit of densely packed structures, the active diffusion coefficient is significantly smaller and depends on the number of particles in the cluster. Our analysis of the cluster-growth dynamics shows distinct differences to processes in systems of short-range attractive colloids in equilibrium. Specifically, the characteristic time for the growth and decay of clusters of a particular size is longer than that of isotropically attractive colloids, which we attribute to the non-reciprocal nature of the directed visual perception. Our simulations reveal a strong interplay between ABP-characteristic interactions, such as volume exclusion and rotational diffusion, and cognitive-based interactions and navigation.

Alignment and propulsion of squirmer pusher-puller dumbbells.

The properties of microswimmer dumbbells composed of pusher-puller pairs are investigated by mesoscale hydrodynamic simulations employing the multiparticle collision dynamics approach for the fluid. An individual microswimmer is represented by a squirmer, and various active-stress combinations in a dumbbell are considered. The squirmers are connected by a bond, which does not impose any geometrical restriction on the individual rotational motion. Our simulations reveal a strong influence of the squirmers' flow fields on the orientation of their propulsion directions, their fluctuations, and the swimming behavior of a dumbbell. The properties of pusher-puller pairs with an equal magnitude of the active stresses depend only weakly on the stress magnitude. This is similar to dumbbells of microswimmers without hydrodynamic interactions. However, for non-equal stress magnitudes, the active stress implies strong orientational correlations of the swimmers' propulsion directions with respect to each other, as well as the bond vector. The orientational coupling is most pronounced for pairs with large differences in the active-stress magnitude. The alignment of the squirmers' propulsion directions with respect to each other is preferentially orthogonal in dumbbells with a strong pusher and weak puller, and antiparallel in the opposite case when the puller dominates. These strong correlations affect the active motion of dumbbells, which is faster for strong pushers and slower for strong pullers.

Tangentially driven active polar linear polymers-An analytical study.

The conformational and dynamical properties of isolated flexible active polar linear polymers (APLPs) are studied analytically. The APLPs are modeled as Gaussian bead-spring linear chains augmented by tangential active forces, both in a discrete and continuous representation. The polar forces lead to linear non-Hermitian equations of motion, which are solved by an eigenfunction expansion in terms of a biorthogonal basis set. Our calculations show that the polymer conformations are independent of activity. However, tangential propulsion strongly impacts the polymer dynamics and yields an active ballistic regime as well as an activity-enhanced long-time diffusive regime, which are both absent in passive systems. The polar forces imply a coupling of modes in the eigenfunction representation, in particular with the translational mode, with a respective strong influence on the polymer dynamics. The total polymer mean-square displacement on scales smaller than the radius of gyration is determined by the active internal dynamics rather than the collective center-of-mass motion, in contrast to active Brownian polymers at large Péclet numbers, reflecting the distinct difference in the propulsion mechanism.

The role of thickness inhomogeneities in hierarchical cortical folding.

The mammalian brain cortex is highly folded, with several developmental disorders affecting folding. On the extremes, lissencephaly, a lack of folds in humans, and polymicrogyria, an overly folded brain, can lead to severe mental retardation, short life expectancy, epileptic seizures, and tetraplegia. Not only a specific degree of folding, but also stereotyped patterns, are required for normal brain function. A quantitative model on how and why these folds appear during the development of the brain is the first step in understanding the cause of these conditions. In recent years, there have been various attempts to understand and model the mechanisms of brain folding. Previous works have shown that mechanical instabilities play a crucial role in the formation of brain folds, and that the geometry of the fetal brain is one of the main factors in dictating its folding characteristics. However, modeling higher-order folding, one of the main characteristics of the highly gyrencephalic brain, has not been fully tackled. The simulations presented in this work are used to study the effects of thickness inhomogeneity in the gyrogenesis of the mammalian brain in silico. Finite-element simulations of rectangular slabs are performed. These slabs are divided into two distinct regions, where the outer region mimicks the gray matter, and the inner region the underlying white matter. Differential growth is introduced by growing the top region tangentially, while keeping the underlying region untouched. The brain tissue is modeled as a neo-Hookean hyperelastic material. Simulations are performed with both homogeneous and inhomogeneous cortical thicknesses. Our results show that the homogeneous cortex folds into a single wavelength, as is common for bilayered materials, while the inhomogeneous cortex folds into more complex conformations. In the early stages of development of the inhomogeneous cortex, structures reminiscent of the deep sulci in the brain are obtained. As the cortex continues to develop, secondary undulations, which are shallower and more variable than the structures obtained in earlier gyrification stage emerge, reproducing well-known characteristics of higher-order folding in the mammalian, and particularly the human, brain.

Effect of cytosol viscosity on the flow behavior of red blood cell suspensions in microvessels.

OBJECTIVE: The flow behavior of blood is strongly affected by red blood cell (RBC) properties, such as the viscosity ratio C between cytosol and suspending medium, which can significantly be altered in several pathologies (e.g. sickle-cell disease, malaria). The main objective of this study is to understand the effect of C on macroscopic blood flow properties such as flow resistance in microvessels, and to link it to the deformation and dynamics of single RBCs. METHODS: We employ mesoscopic hydrodynamic simulations to investigate flow properties of RBC suspensions with different cytosol viscosities for various flow conditions in cylindrical microchannels. RESULTS: Starting from a dispersed cell configuration which approximates RBC dispersion at vessel bifurcations in the microvasculature, we find that the flow convergence and development of RBC-free layer (RBC-FL) depend only weakly on C, and require a convergence length in the range of 25D-50D, where D is channel diameter. In vessels with D≤20µm , the final resistance of developed flow is nearly the same for C = 5 and C = 1, while for D=40µm , the flow resistance for C = 5 is about 10% larger than for C = 1. The similarities and differences in flow resistance can be explained by viscosity-dependent RBC-FL thicknesses, which are associated with the viscosity-dependent dynamics of single RBCs. CONCLUSIONS: The weak effect on the flow resistance and RBC-FL explains why RBCs can contain a high concentration of hemoglobin for efficient oxygen delivery, without a pronounced increase in the flow resistance. Furthermore, our results suggest that significant alterations in microvascular flow in various pathologies are likely not due to mere changes in cytosolic viscosity.

Correction: Hydrodynamic interactions in squirmer dumbbells: active stress-induced alignment and locomotion.

Correction for 'Hydrodynamic interactions in squirmer dumbbells: active stress-induced alignment and locomotion' by Judit Clopés et al., Soft Matter, 2020, 16, 10676-10687, DOI: .

Multi-ciliated microswimmers-metachronal coordination and helical swimming.

The dynamics and motion of multi-ciliated microswimmers with a spherical body and a small number N (with [Formula: see text]) of cilia with length comparable to the body radius, is investigated by mesoscale hydrodynamics simulations. A metachronal wave is imposed for the cilia beat, for which the wave vector has both a longitudinal and a latitudinal component. The dynamics and motion is characterized by the swimming velocity, its variation over the beat cycle, the spinning velocity around the main body axis, as well as the parameters of the helical trajectory. Our simulation results show that the microswimmer motion strongly depends on the latitudinal wave number and the longitudinal phase lag. The microswimmers are found to swim smoothly and usually spin around their own axis. Chirality of the metachronal beat pattern generically generates helical trajectories. In most cases, the helices are thin and stretched, i.e., the helix radius is about an order of magnitude smaller than the pitch. The rotational diffusion of the microswimmer is significantly smaller than the passive rotational diffusion of the body alone, which indicates that the extended cilia contribute strongly to the hydrodynamic radius. The swimming velocity is found to increase with the cilia number N with a slightly sublinear power law, consistent with the behavior expected from the dependence of the transport velocity of planar cilia arrays on the cilia separation.

Wall-anchored semiflexible polymer under large amplitude oscillatory shear flow.

The properties of semiflexible polymers tethered by one end to an impenetrable wall and exposed to oscillatory shear flow are investigated by mesoscale simulations. A polymer, confined in two dimensions, is described by a linear bead-spring chain, and fluid interactions are incorporated by the Brownian multiparticle collision dynamics approach. At small strain, the polymers follow the applied flow field. However, at high strain, we find a strongly nonlinear response with major conformational changes. Polymers are stretched along the flow direction and exhibit U-shaped conformations while following the flow. As a consequence of confinement in the half-space, frequency doubling in the time-dependent polymer properties appears along the direction normal to the wall.

Active bath-induced localization and collapse of passive semiflexible polymers.

The conformational and dynamical properties of a passive polymer embedded in a bath of active Brownian particles (ABPs) are studied by Langevin dynamics simulations. Various activities and ABP concentrations below and above the critical values for motility-induced phase separation (MIPS) are considered. In a homogeneous ABP fluid, the embedded polymer swells with increasing bath activity, with stronger swelling for larger densities. The polymer dynamics is enhanced, with the diffusion coefficient increasing by a power-law with increasing activity, where the exponent depends on the ABP concentration. For ABP concentrations in the MIPS regime, we observe a localization of the polymer in the low-density ABP phase associated with polymer collapse for moderate activities and a reswelling for high activities accompanied by a preferred localization in the high-density ABP phase. Localization and reswelling are independent of the polymer stiffness, with stiff polymers behaving similarly to flexible polymers. The polymer collapse is associated with a slowdown of its dynamics and a significantly smaller center-of-mass diffusion coefficient. In general, the polymer dynamics can only partially be described by an effective (bath) temperature. Moreover, the properties of a polymer embedded in a homogeneous active bath deviate quantitatively from those of a polymer composed of active monomers, i.e., linear chains of ABPs; however, such a polymer exhibits qualitatively similar activity-dependent features.

Osmotic Concentration-Controlled Particle Uptake and Wrapping-Induced Lysis of Cells and Vesicles.

In vivo, high protein and ion concentrations determine the preferred volumes of cells, organelles, and vesicles. Deformations of their lipid-bilayer membranes by nanoparticle wrapping reduce the interior volumes available to solutes and thus induce large osmotic pressure differences. Osmotic concentration can therefore be an important control parameter for wrapping of nanoparticles. We employ a curvature-elasticity model of the membrane and contact interaction with spherical particles to study their wrapping at initially spherical vesicles. Although the continuous particle-binding transition is independent of the presence of solutes, the discontinuous envelopment transition shifts to higher adhesion strengths and the corresponding energy barrier increases with increasing osmotic concentration. High osmotic concentrations stabilize partial-wrapped, membrane-bound states for both, particle attachment to the inside and the outside. In this regime, wrapping of particles controls membrane tension, with power-law dependencies on osmotic concentration and adhesion strength. For high adhesion strengths, particle wrapping can lead to the opening of mechanosensitive channels in cell membranes and to lysis. Membrane tension-induced stabilization of partial-wrapped states as well as wrapping-induced lysis play important roles not only for desired mechano-bacteriocidal effects of engineered nanomaterials but may also determine viral burst sizes of bacteria and control endocytosis for mammalian cells.

Enhanced Rotational Motion of Spherical Squirmer in Polymer Solutions.

The rotational diffusive motion of a self-propelled, attractive spherical colloid immersed in a solution of self-avoiding polymers is studied by mesoscale hydrodynamic simulations. A drastic enhancement of the rotational diffusion by more than an order of magnitude in the presence of activity is obtained. The amplification is a consequence of two effects, a decrease of the amount of adsorbed polymers by active motion and an asymmetric encounter with polymers on the squirmer surface, which yields an additional torque and random noise for the rotational motion. Our simulations suggest a way to control the rotational dynamics of squirmer-type microswimmers by the degree of polymer adsorption and system heterogeneity.

Chiral-filament self-assembly on curved manifolds.

Rod-like and banana-shaped proteins, like BAR-domain proteins and MreB proteins, adsorb on membranes and regulate the membrane curvature. The formation of large filamentous complexes of these proteins plays an important role in cellular processes like membrane trafficking, cytokinesis and cell motion. We propose a simplified model to investigate such curvature-dependent self-assembly processes. Anisotropic building blocks, modeled as trimer molecules, which have a preferred binding site, interact via pair-wise Lennard-Jones potentials. When several trimers assemble, they form an elastic ribbon with an intrinsic curvature and twist, controlled by bending and torsional rigidity. For trimer self-assembly on the curved surface of a cylindrical membrane, this leads to a preferred spatial orientation of the ribbon. We show that these interactions can lead to the formation of helices with several windings around the cylinder. The emerging helix angle and pitch depend on the rigidities and the intrinsic curvature and twist values. In particular, a well-defined and controllable helix angle emerges in the case of equal bending and torsional rigidity. The dynamics of filament growth is characterized by three regimes, in which filament length increases with the power laws tz in time, with z≃ 3/4, z = 1/2, and z = 0 for short, intermediate, and long times, respectively. A comparison with the solutions of the Smoluchowski aggregation equation allows the identification of the underlying mechanism in the short-time regime as a crossover from size-independent to diffusion-limited aggregation. Thus, helical structures, as often observed in biology, can arise by self-assembly of anisotropic and chiral proteins.