skip to main content

Title: Swinging and tumbling of multicomponent vesicles in flow
Biological membranes are host to proteins and molecules which may form domain-like structures resulting in spatially varying material properties. Vesicles with such heterogeneous membranes can exhibit intricate shapes at equilibrium and rich dynamics when placed into a flow. Under the assumption of small deformations and a two-dimensional system, we develop a reduced-order model to describe the fluid-structure interaction between a viscous background shear flow and an inextensible membrane with spatially varying bending stiffness and spontaneous curvature. Material property variations of a critical magnitude, relative to the flow rate and internal/external viscosity contrast, can set off a qualitative change in the vesicle dynamics. A membrane of nearly constant bending stiffness or spontaneous curvature undergoes a small amplitude swinging motion (which includes tangential tank-treading), while for large enough material variations the dynamics pass through a regime featuring tumbling and periodic phase-lagging of the membrane material, and ultimately for very large material variation to a rigid-body tumbling behaviour. Distinct differences are found for even and odd spatial modes of domain distribution. Full numerical simulations are used to probe the theoretical predictions, which appear valid even when studying substantially deformed membranes.
; ;
Award ID(s):
2003819 1661900
Publication Date:
Journal Name:
Journal of Fluid Mechanics
Sponsoring Org:
National Science Foundation
More Like this
  1. Protein domains, such as ENTH (epsin N-terminal homology) and BAR (bin/amphiphysin/rvs), contain amphipathic helices that drive preferential binding to curved membranes. However, predicting how the physical parameters of these domains control this ‘curvature sensing’ behavior is challenging due to the local membrane deformations generated by the nanoscopic helix on the surface of a large sphere. We here use a deformable continuum model that accounts for the physical properties of the membrane and the helix insertion to predict curvature sensing behavior, with direct validation against multiple experimental datasets. We show that the insertion can be modeled as a local change to the membrane's spontaneous curvature, c ins0, producing excellent agreement with the energetics extracted from experiments on ENTH binding to vesicles and cylinders, and of ArfGAP helices to vesicles. For small vesicles with high curvature, the insertion lowers the membrane energy by relieving strain on a membrane that is far from its preferred curvature of zero. For larger vesicles, however, the insertion has the inverse effect, de-stabilizing the membrane by introducing more strain. We formulate here an empirical expression that accurately captures numerically calculated membrane energies as a function of both basic membrane properties (bending modulus κ and radius R )more »as well as stresses applied by the inserted helix ( c ins0 and area A ins ). We therefore predict how these physical parameters will alter the energetics of helix binding to curved vesicles, which is an essential step in understanding their localization dynamics during membrane remodeling processes.« less
  2. We show in the companion paper that the free membrane shape of lipid bilayer vesicles containing the mechanosensitive ion channel Piezo can be predicted, with no free parameters, from membrane elasticity theory together with measurements of the protein geometry and vesicle size [C. A. Haselwandter, Y. R. Guo, Z. Fu, R. MacKinnon, Proc. Natl. Acad. Sci. U.S.A. , 10.1073/pnas.2208027119 (2022)]. Here we use these results to determine the force that the Piezo dome exerts on the free membrane and hence, that the free membrane exerts on the Piezo dome, for a range of vesicle sizes. From vesicle shape measurements alone, we thus obtain a force–distortion relationship for the Piezo dome, from which we deduce the Piezo dome’s intrinsic radius of curvature, 42 ± 12 nm, and bending stiffness, 18 ± 2.1   k B T , in freestanding lipid bilayer membranes mimicking cell membranes. Applying these estimates to a spherical cap model of Piezo embedded in a lipid bilayer, we suggest that Piezo’s intrinsic curvature, surrounding membrane footprint, small stiffness, and large area are the key properties of Piezo that give rise to low-threshold, high-sensitivity mechanical gating.
  3. Cellular membranes are elastic lipid bilayers that contain a variety of proteins, including ion channels, receptors and scaffolding proteins. These proteins are known to diffuse in the plane of the membrane and to influence the bending of the membrane. Experiments have shown that lipid flow in the plane of the membrane is closely coupled with the diffusion of proteins. Thus, there is a need for a comprehensive framework that accounts for the interplay between these processes. Here, we present a theory for the coupled in-plane viscous flow of lipids, diffusion of transmembrane proteins and elastic deformation of lipid bilayers. The proteins in the membrane are modelled such that they influence membrane bending by inducing a spontaneous curvature. We formulate the free energy of the membrane with a Helfrich-like curvature elastic energy density function modified to account for the chemical potential energy of proteins. We derive the conservation laws and equations of motion for this system. Finally, we present results from dimensional analysis and numerical simulations and demonstrate the effect of coupled transport processes in governing the dynamics of membrane bending and protein diffusion.
  4. Many cellular lipid bilayers consist of leaflets that differ in their lipid composition — a non-equilibrium state actively maintained by cellular sorting processes that counter passive lipid flip-flop. While this lipidomic aspect of membrane asymmetry has been known for half a century, its elastic and thermodynamic ramifications have garnered attention only fairly recently. Notably, the torque arising when lipids of different spontaneous curvature reside in the two leaflets can be counterbalanced by a difference in lateral mechanical stress between them. Such membranes can be essentially flat in their relaxed state, despite being compositionally strongly asymmetric, but they harbor a surprisingly large but macroscopically invisible differential stress. This hidden stress can affect a wide range of other membrane properties, such as the resistance to bending, the nature of phase transitions in its leaflets, and the distribution of flippable species, most notably sterols. In this short note we offer a concise overview of our recently proposed basic framework for capturing the interplay between curvature, lateral stress, leaflet phase behavior, and cholesterol distribution in generally asymmetric membranes, and how its implied signatures might be used to learn more about the hidden but physically consequential differential stress.
  5. Abstract

    Three-dimensional dynamics of flexible fibers in shear flow are studied numerically, with a qualitative comparison to experiments. Initially, the fibers are straight, with different orientations with respect to the flow. By changing the rotation speed of a shear rheometer, we change the ratioAof bending to shear forces. We observe fibers in the flow-vorticity plane, which gives insight into the motion out of the shear plane. The numerical simulations of moderately flexible fibers show that they rotate along effective Jeffery orbits, and therefore the fiber orientation rapidly becomes very close to the flow-vorticity plane, on average close to the flow direction, and the fiber remains in an almost straight configuration for a long time. This ‘ordering’ of fibers is temporary since they alternately bend and straighten while tumbling. We observe numerically and experimentally that if the fibers are initially in the compressional region of the shear flow, they can undergo compressional buckling, with a pronounced deformation of shape along their whole length during a short time, which is in contrast to the typical local bending that originates over a long time from the fiber ends. We identify differences between local and compressional bending and discuss their competition, which depends onmore »the initial orientation of the fiber and the bending stiffness ratioA. There are two main finding. First, the compressional buckling is limited to a certain small range of the initial orientations, excluding those from the flow-vorticity plane. Second, since fibers straighten in the flow-vorticity plane while tumbling, the compressional buckling is transient—it does not appear for times longer than 1/4 of the Jeffery period. For larger times, bending of fibers is always driven by their ends.

    « less