We report a numerical study addressing the dynamics of compound vesicles confined in a channel under shear flow. The system comprises a smaller vesicle embedded within a larger one and can be used to mimic, for example, leukocytes or nucleate cells. A two-dimensional model, which combines molecular dynamics and mesoscopic hydrodynamics including thermal fluctuations, is adopted to perform an extended investigation. We are able to vary independently the swelling degree and the relative size of vesicles, the viscosities of fluids internal and external to vesicles, and the Capillary number, so to observe a rich dynamical phenomenology which goes well beyond what observed for single vesicles, matching quantitatively with experimental findings. Tank-treading, tumbling, and trembling motions are enriched by dynamical states where inner and outer vesicles can perform different motions. We show that thermal fluctuations are crucial during trembling and swinging dynamics, as observed in experiments. Undulating motion of the external vesicle, characterized by periodic oscillation of the inclination and buckling of the membrane, is observed at high filling fractions. This latter state exhibit
Motivated by recent studies of two-phase lipid vesicles possessing 2D solid domains integrated within a fluid bilayer phase, we study the shape equilibria of closed vesicles possessing a single planar, circular inclusion. While 2D solid elasticity tends to expel Gaussian curvature, topology requires closed vesicles to maintain an average, non-zero Gaussian curvature leading to an elementary mechanism of shape frustration that increases with inclusion size. We study elastic ground states of the Helfrich model of the planar-fluid composite vesicles, analytically and computationally, as a function of planar fraction and reduced volume. Notably, we show that incorporation of a planar inclusion of only a few percent dramatically shifts the ground state shapes of vesicles from predominantly {\it prolate} to {\it oblate}, and moreover, shifts the optimal surface to volume ratio far from spherical shapes. We show that for sufficiently small planar inclusions, the elastic ground states break symmetry via a complex variety of asymmetric oblate, prolate, and triaxial shapes, while inclusion sizes above about $8\%$ drive composite vesicles to adopt axisymmetric oblate shapes. These predictions
The wrapping of particles and vesicles by lipid bilayer membranes is a fundamental process in cellular transport and targeted drug delivery. Here, we investigate the wrapping behavior of nonspherical vesicles, such as ellipsoidal, prolate, oblate, and stomatocytes, by systematically varying the bending rigidity of the vesicle membrane and the tension of the planar membrane. Using the Helfrich Hamiltonian, triangulated membrane models, and energy minimization techniques, we predict multiple stable wrapping states and identify the conditions for their coexistence. Our results demonstrate that softer vesicles bind more easily to planar membranes; however, achieving complete wrapping requires significantly higher adhesion strengths compared to rigid particles. As membrane tension increases, deep-wrapped states disappear at a triple point where shallow-wrapped, deep-wrapped, and complete-wrapped states coexist. The coordinates of the triple point are highly sensitive to the vesicle shape and stiffness. For stomatocytes, increasing stiffness shifts the triple point to higher adhesion strengths and membrane tensions, while for oblates it shifts to lower values, influenced by shape changes
Uncovering the rules governing the nonequilibrium dynamics of the membranes that define biological cells is of central importance to understanding the physics of living systems. We theoretically and computationally investigate the behavior of flexible quasispherical vesicles that exchange membrane constituents, internal volume, and heat with an external reservoir. The excess chemical potential and osmotic pressure difference imposed by the reservoir act as generalized thermodynamic driving forces that modulate vesicle morphology. We show that the renormalization of membrane mechanical properties by nonequilibrium driving gives rise to a morphological transition between a weakly driven regime, in which growing vesicles remain quasispherical, and a strongly driven regime, in which vesicles accommodate rapid membrane uptake by developing surface wrinkles. Additionally, we propose a minimal vesicle growth-shape law, derived using insights from stochastic thermodynamics, that robustly describes vesicle growth dynamics even in strongly driven, far-from-equilibrium regimes.
Vesicles are important surrogate structures made up of multiple phospholipids and cholesterol distributed in the form of a lipid bilayer. Tubular vesicles can undergo pearling i.e., formation of beads on the liquid thread akin to the Rayleigh-Plateau instability. Previous studies have inspected the effects of surface tension on the pearling instabilities of single-component vesicles. In this study, we perform a linear stability analysis on a multicomponent cylindrical vesicle. We solve the Stokes equations along with the Cahn-Hilliard equations to develop the linearized dynamic equations governing the vesicle shape and surface concentration fields. This helps us show that multicomponent vesicles can undergo pearling, buckling, and wrinkling even in the absence of surface tension, which is a significantly different result from studies on single-component vesicles. This behaviour arises due to the competition between the free energies of phase separation, line tension, and bending for this multi-phospholipid system. We determine the conditions under which axisymmetric and non-axisymmetric modes are dominant, and supplement our results with an energy analysis that shows the sources fo
We examine the deformation of homogeneous spherical fluid vesicles along their equator by a circular rigid ring. We consider deformations preserving the axial and equatorial mirror symmetries of the vesicles. The configurations of the vesicle are determined employing the spontaneous curvature model subject to the constraints imposed by the ring as well as of having constant area or volume. We determine two expressions of the force exerted by the ring, one involving a discontinuity in the derivative of the curvature of the membrane across the ring, and another one in terms of the global quantities of the vesicle. For small enough values of the spontaneous curvature there is only one sequence of configurations for either fixed area or volume. The behavior of constricted vesicles is similar for both constraints, they follow a transition from prolate to dumbbell shapes, which culminates in two quasispherical vesicles connected by a small catenoid-like neck. We analyze the geometry and the force of the small neck employing a perturbative analysis about the catenoid. A stretched vesicle initially adopts an oblate shape for either constraint. If the area is fixed the vesicle increasingly
Polymer vesicles are stable robust vesicles made from block copolymer amphiphiles. Recent progress in the chemical design of block copolymers opens up the exciting possibility of creating a wide variety of polymer vesicles with varying fine structure, functionality and geometry. Polymer vesicles not only constitute useful systems for drug delivery and micro/nano-reactors but also provide an invaluable arena for exploring the ordering of matter on curved surfaces embedded in three dimensions. By choosing suitable liquid-crystalline polymers for one of the copolymer components one can create vesicles with smectic stripes. Smectic order on shapes of spherical topology inevitably possesses topological defects (disclinations) that are themselves distinguished regions for potential chemical functionalization and nucleators of vesicle budding. Here we report on glassy striped polymer vesicles formed from amphiphilic block copolymers in which the hydrophobic block is a smectic liquid crystal polymer containing cholesteryl-based mesogens. The vesicles exhibit two-dimensional smectic order and are ellipsoidal in shape with defects, or possible additional budding into isotropic vesicles, at t
The interaction of surfaces in relative motion in wet environments is dominated by lubrication forces, which play a pivotal role in the dynamics of microscopic systems. Here, we develop motile vesicles that exploit lubrication forces to roll on substrates. The activity of the vesicle comes from the confined rotational flow generated by a driven rotating particle encapsulated within the vesicle by droplet-microfluidics. Lubrication forces driving vesicle rolling are controlled by membrane mechanics and its tribological properties. This provides the design principles for motile vesicles that exploit frictional forces to efficiently navigate through complex environments.
Recent experiments on vesicles formed from block copolymers with liquid-crystalline side-chains reveal a rich variety of vesicle morphologies. The additional internal order ("structure") developed by these self-assembled block copolymer vesicles can lead to significantly deformed vesicles as a result of the delicate interplay between two-dimensional ordering and vesicle shape. The inevitable topological defects in structured vesicles of spherical topology also play an essential role in controlling the final vesicle morphology. Here we develop a minimal theoretical model for the morphology of the membrane structure with internal nematic/smectic order. Using both analytic and numerical approaches, we show that the possible low free energy morphologies include nano-size cylindrical micelles (nano-fibers), faceted tetrahedral vesicles, and ellipsoidal vesicles, as well as cylindrical vesicles. The tetrahedral vesicle is a particularly fascinating example of a faceted liquid-crystalline membrane. Faceted liquid vesicles may lead to the design of supra-molecular structures with tetrahedral symmetry and new classes of nano-carriers.
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 arr
Translocation of a daughter vesicle from a mother vesicle through a pore is experimentally studied by many groups using a model system of self-reproducing vesicles. However, the theoretical formulation of the problem is not fully understood. In the present study, we present a theoretical formulation of the process based on our previous work [P. Khunpetch et al., Phys. Fluids 33, 077103 (2021)]. In our previous work, we considered the daughter vesicle as a rigid body. In the present work, however, we allow the daughter vesicle to deform during the expulsion process. We thus derive the free energy constituting of the elastic moduli of both the mother and daughter vesicles, and of pressure-driven contribution. The minimum energy path of the translocation is searched by using the string method. Our improved model successfully suggests the disappearance of the energy barrier where all the elastic moduli are in agreement with the experimental reports, while the previous work is unsuccessful to do so. The equations of motion of the daughter vesicle have been derived within the framework of the Onsager principle. We found that the translocation time of the daughter vesicle can be reduced w
The dynamics of fluid vesicles is studied under flow in microchannels, in which the width varies periodically along the channel. Three types of flow instabilities of prolate vesicles are found. For small quasi-spherical vesicles -- compared to the average channel width -- perturbation theory predicts a transition from a state with orientational oscillations of a fixed prolate shape to a state with shape oscillations of symmetrical ellipsoidal or bullet-like shapes with increasing flow velocity. Experimentally, such orientational oscillations are observed during the slow migration of a vesicle towards the centerline of the channel. For larger vesicles, mesoscale hydrodynamics simulations and experiments show similar symmetric shape oscillation at reduced volumes $V^* \gtrsim 0.9$. However, for non-spherical vesicles with $V^* \lesssim 0.9$, shapes are found with two symmetric or a single asymmetric tail.
We investigated the thermodynamic relaxation of giant unilamellar vesicles (GUVs) which contained small vesicles within their interior. Quenching these vesicles from their fluid phase (T>Tm) through the phase transition in the gel state (T<Tm) drives the inner vesicles to be expelled from the larger mother vesicle via the accompanying decrease in vesicle area by ~25% which forces a pore to open in the mother vesicle. We demonstrate that the proceeding time evolution of the resulting efflux follows the relaxation of the membrane area and describe the entire relaxation process using an Onsager-like nonequilibrium thermodynamics ansatz. As a consequence of the volume efflux internal vesicles are expelled from the mother vesicle. Although complete sealing of the pore may occur during the expulsion, the global relaxation dynamics is conserved. Finally, comparison of these results to morphological relaxation phenomena found in earlier studies reveals a universal relaxation behaviour in GUVs. When quenched from the fluid to gel phase the typical time scale of relaxation shows little variation and ranges between 4-5 seconds.
Erythrocytes and vesicles shed by erythrocytes from human and rat blood were collected and analyzed after temperature was elevated by physical exercise or by exposure to external heat. The images of erythrocytes and vesicles were analyzed by the light microscopy system with spatial resolution of better than 90 nm. The samples were observed in an aqueous environment and required no freezing, dehydration, staining, shadowing, marking or any other manipulation. Temperature elevation, whether passive or through exercise, resulted in significant concentration increase of structurally transformed erythrocytes (echinocytes) and vesicles in blood. At temperature of 37 oC, mean vesicle concentrations and diameters in human and rat blood were (1.50+-0.35)x10^6 and (1.4+-0.2)x10^6 vesicles/μL, and 0.365+-0.065 and 0.436+-0.03 μm, respectively. It was estimated that 80% of all vesicles found in human blood are smaller than 0.4 μm. Thermodynamic analysis of experimental and literature data showed that erythrocyte transformation, vesicle release and other associated processes are driven by entropy with enthalpy-entropy compensation. It is suggested that physical state of hydrated cell membrane i
Piezo ion channels underlie many forms of mechanosensation in vertebrates, and have been found to bend the membrane into strongly curved dome shapes. We develop here a methodology describing the self-assembly of lipids and Piezo into polyhedral bilayer vesicles. We validate this methodology for bilayer vesicles formed from bacterial mechanosensitive channels of small conductance, for which experiments found a polyhedral arrangement of proteins with snub cube symmetry and a well-defined characteristic vesicle size. On this basis, we calculate the self-assembly diagram for polyhedral bilayer vesicles formed from Piezo. We find that the radius of curvature of the Piezo dome provides a critical control parameter for the self-assembly of Piezo vesicles, with high abundances of Piezo vesicles with octahedral, icosahedral, and snub cube symmetry with increasing Piezo dome radius of curvature.
In recent experiments, small, thick-walled vesicles with a narrow size distribution were formed from copolymers where the degree of polymerisation of the hydrophobic block, N_B, was significantly greater than that of the hydrophilic block, N_A. Using a mean-field theory, we reproduce several aspects of the behaviour of these vesicles. Firstly, we find a minimum in the free energy of the system of vesicles as a function of their radius, corresponding to a preferred size for the vesicles, when N_B is several times larger than N_A. Furthermore, the vesicle radius diverges as N_B is increased towards a critical value, consistent with the instability of the vesicles with respect to further aggregation seen in the experimental work. Finally, we find that this instability can also be triggered in our model by changing the interaction strength of the copolymers with the solvent.
Biological membranes constantly change their shape in response to external stimuli, and understanding the remodeling and stability of vesicles in heterogeneous environments is therefore of fundamental importance for a range of cellular processes. One crucial question is how vesicles respond to external osmotic stresses, imposed by differences in solute concentrations between the vesicle interior and exterior. Previous analyses of the membrane bending energy have predicted that micron-sized giant unilamellar vesicles (GUVs) should become globally deformed already for nanomolar concentration differences, in contrast to experimental findings that find deformations at much higher osmotic stresses. In this article, we analyze the mechanical stability of a spherical vesicle exposed to an external osmotic pressure in a statistical-mechanical model, including the effect of thermally excited membrane bending modes. We find that the inclusion of thermal fluctuations of the vesicle shape changes renders the vesicle deformation continuous, in contrast to the abrupt transition in the athermal picture. Crucially, however, the predicted critical pressure associated with global vesicle deformation
We investigate the influence of a non-ionic surfactant like polymer on phospholipid vesicles. Our results from cryogenic transmission electron microscopy (cryo-TEM), dynamic light scattering (DLS), small angle neutron and X-ray scattering (SANS/SAXS), identifies the existence of multilayer vesicles and an increase in size of the vesicles in presence of the polymers. We present a generalized model to obtain the bending rigidity from neutron spin echo spectroscopy (NSE) data for multilayer vesicles. We demonstrated that polymers are trapped in the lipid bilayer, causing a partial disruption in the vesicle, which is attributed to the reduction in bending rigidity per unit bilayer. We also observed substantial dampening of the trapped lipid tail motion in presence of the polymer. Our results highlighted the possibilities of using specialized polymers that can disrupt membrane and control their dynamics with possible application in topical drug or nutraceutical formulations.
The intracellular environment is a dynamic space filled with various organelles moving in all directions. Included in this diverse group of organelles are vesicles, which are involved in transport of molecular cargo throughout the cell. Vesicles move in either a directed or non-directed fashion, often depending on interactions with cytoskeletal proteins such as microtubules, actin filaments, and molecular motors. How these proteins affect the local fluctuations of vesicles in the cytoplasm is not clear since they have the potential to both facilitate and impede movement. Here we show that vesicle mobility is significantly affected by myosin-II, even though it is not a cargo transport motor. We find that myosin-II activity increases the effective diffusivity of vesicles and its inhibition facilitates longer states of non-directed motion. Our study suggests that altering myosin-II activity in the cytoplasm of cells can modulate the mobility of vesicles, providing a possible mechanism for cells to dynamically tune the cytoplasmic environment in space and time.
In this work we demonstrate how the first order phase transition in giant unilamellar vesicles (GUVs) can function as a trigger for membrane fission. When driven through their gel-fluid phase transition GUVs exhibit budding or pearl formation. These buds remain connected to the mother vesicle presumably by a small neck. Cooling these vesicles from the fluid phase (T>Tm) through the phase transition into the gel state (T<Tm), leads to complete rupture and fission of the neck, while the mother vesicle remains intact. Pearling tubes which formed upon heating break-up and decay into multiple individual vesicles which then diffuse freely. Finally we demonstrate that mimicking the intracellular bulk viscosity by increasing the bulk viscosity to 40cP does not affect the overall fission process, but leads to a significant decrease in size of the released vesicles.