Theoretical study of vesicle shapes driven by coupling curved proteins and active cytoskeletal forces

Stability of a biomembrane tube covered with proteins

Membrane tubes are essential structural features in cells that facilitate biomaterial transport and inter- and intracellular signaling. The shape of these tubes can be regulated by the proteins that surround and adhere to them. We study the stability of a biomembrane tube coated with proteins by combining linear stability analysis, out-of-equilibrium hydrodynamic calculations, and numerical solutions of a Helfrich-like membrane model. Our analysis demonstrates that both long- and short-wavelength perturbations can destabilize the tubes. Numerical simulations confirm the derived linear stability criteria and yield the nonlinearly perturbed vesicle shapes. Our study highlights the interplay between membrane shape and protein density, where the shape instability concurs with a redistribution of proteins into a banded pattern.

A minimal physical model for curvotaxis driven by curved protein complexes at the cell's leading edge

Cells often migrate on curved surfaces inside the body, such as curved tissues, blood vessels, or highly curved protrusions of other cells. Recent in vitro experiments provide clear evidence that motile cells are affected by the curvature of the substrate on which they migrate, preferring certain curvatures to others, termed "curvotaxis." The origin and underlying mechanism that gives rise to this curvature sensitivity are not well understood. Here, we employ a "minimal cell" model which is composed of a vesicle that contains curved membrane protein complexes, that exert protrusive forces on the membrane (representing the pressure due to actin polymerization). This minimal-cell model gives rise to spontaneous emergence of a motile phenotype, driven by a lamellipodia-like leading edge. By systematically screening the behavior of this model on different types of curved substrates (sinusoidal, cylinder, and tube), we show that minimal ingredients and energy terms capture the experimental data. The model recovers the observed migration on the sinusoidal substrate, where cells move along the grooves (minima), while avoiding motion along the ridges. In addition, the model predicts the tendency of cells to migrate circumferentially on convex substrates and axially on concave ones. Both of these predictions are verified experimentally, on several cell types. Altogether, our results identify the minimization of membrane-substrate adhesion energy and binding energy between the membrane protein complexes as key players of curvotaxis in cell migration.

Surfactin molecules with a cone-like structure promote the formation of membrane domains with negative spontaneous curvature and induce membrane invaginations

Surfactin uniquely influences lipid bilayer structure by initially inducing membrane invaginations before solubilization. In this study, we exposed DOPC giant vesicles to various surfactin concentrations at different temperatures and observed surfactin-induced membrane invaginations by using differential interference contrast and confocal laser fluorescence microscopy. These invaginations were stable at room temperature but not at higher temperatures. Surfactin molecules induce membrane nanodomains with negative spontaneous curvature and membrane invaginations despite their intrinsic conical shape and intrinsic positive curvature. Considering the experimentally observed capacity of surfactin to fluidize lipid acyl chains and induce partial dehydration of lipid headgroups, we propose that the resulting surfactin-lipid complexes exhibit a net negative spontaneous curvature. We further conducted 3D numerical Monte Carlo (MC) simulations to investigate the behaviour of vesicles containing negative curvature nanodomains within their membrane at varying temperatures. MC simulations demonstrated strong agreement with experimental results, revealing that invaginations are preferentially formed at low temperatures, while being less pronounced at elevated temperatures. Our findings go beyond the expectations of the Israelachvili molecular shape and packing concepts analysis. These concepts do not take into account the influence of specific interactions between neighboring molecules on the inherent shapes of molecules and their arrangement within curved membrane nanodomains. Our work contributes to a more comprehensive understanding of the complex factors governing vesicle morphology and membrane organization and provides insight into the role of detergent-lipid interactions in modulating vesicle morphology.

Investigation of nano- and microdomains formed by ceramide 1 phosphate in lipid bilayers

Biological membranes are renowned for their intricate complexity, with the formation of membrane domains being pivotal to the successful execution of numerous cellular processes. However, due to their nanoscale characteristics, these domains are often understudied, as the experimental techniques required for quantitative investigation present significant challenges. In this study we employ spot-variation z-scan fluorescence correlation spectroscopy (svzFCS) tailored for artificial lipid vesicles of varying composition and combine this approach with high-resolution imaging. This method has been harnessed to examine the lipid-segregation behavior of distinct types of ceramide-1-phosphate (C1P), a crucial class of signaling molecules, within these membranes. Moreover, we provide a quantitative portrayal of the lipid membranes studied and the domains induced by C1P at both nano and microscales. Given the lack of definitive conclusions from the experimental data obtained, it was supplemented with comprehensive in silico studies-including the analysis of diffusion coefficient via molecular dynamics and domain populations via Monte Carlo simulations. This approach enhanced our insight into the dynamic behavior of these molecules within model lipid membranes, confirming that nano- and microdomains can co-exist in lipid vesicles.

Shape transitions in a network model of active elastic shells

Morphogenesis involves the transformation of initially simple shapes, such as multicellular spheroids, into more complex 3D shapes. These shape changes are governed by mechanical forces including molecular motor-generated forces as well as hydrostatic fluid pressure, both of which are actively regulated in living matter through mechano-chemical feedback. Inspired by autonomous, biophysical shape change, such as occurring in the model organism hydra, we introduce a minimal, active, elastic model featuring a network of springs in a globe-like spherical shell geometry. In this model there is coupling between activity and the shape of the shell: if the local curvature of a filament represented by a spring falls below a critical value, its elastic constant is actively changed. This results in deformation of the springs that changes the shape of the shell. By combining excitation of springs and pressure regulation, we show that the shell undergoes a transition from spheroidal to either elongated ellipsoidal or a different spheroidal shape, depending on pressure. There exists a critical pressure at which there is an abrupt change from ellipsoids to spheroids, showing that pressure is potentially a sensitive switch for material shape. We thus offer biologically inspired design principles for autonomous shape transitions in active elastic shells.

Experimental and theoretical model for the origin of coiling of cellular protrusions around fibers

Protrusions at the leading-edge of a cell play an important role in sensing the extracellular cues during cellular spreading and motility. Recent studies provided indications that these protrusions wrap (coil) around the extracellular fibers. However, the physics of this coiling process, and the mechanisms that drive it, are not well understood. We present a combined theoretical and experimental study of the coiling of cellular protrusions on fibers of different geometry. Our theoretical model describes membrane protrusions that are produced by curved membrane proteins that recruit the protrusive forces of actin polymerization, and identifies the role of bending and adhesion energies in orienting the leading-edges of the protrusions along the azimuthal (coiling) direction. Our model predicts that the cell's leading-edge coils on fibers with circular cross-section (above some critical radius), but the coiling ceases for flattened fibers of highly elliptical cross-section. These predictions are verified by 3D visualization and quantitation of coiling on suspended fibers using Dual-View light-sheet microscopy (diSPIM). Overall, we provide a theoretical framework, supported by experiments, which explains the physical origin of the coiling phenomenon.

Molecular Relay Stations in Membrane Nanotubes: IRSp53 Involved in Actin-Based Force Generation

Membrane nanotubes are cell protrusions that grow to tens of micrometres and functionally connect cells. Actin filaments are semi-flexible polymers, and their polymerisation provides force for the formation and growth of membrane nanotubes. The molecular bases for the provision of appropriate force through such long distances are not yet clear. Actin filament bundles are likely involved in these processes; however, even actin bundles weaken when growing over long distances, and there must be a mechanism for their regeneration along the nanotubes. We investigated the possibility of the formation of periodic molecular relay stations along membrane nanotubes by describing the interactions of actin with full-length IRSp53 protein and its N-terminal I-BAR domain. We concluded that I-BAR is involved in the early phase of the formation of cell projections, while IRSp53 is also important for the elongation of protrusions. Considering that IRSp53 binds to the membrane along the nanotubes and nucleates actin polymerisation, we propose that, in membrane nanotubes, IRSp53 establishes molecular relay stations for actin polymerisation and, as a result, supports the generation of force required for the growth of nanotubes.

A minimal cell model for lamellipodia-based cellular dynamics and migration

One ubiquitous cellular structure for performing various tasks, such as spreading and migration over external surfaces, is the sheet-like protrusion called a lamellipodium, which propels the leading edge of the cell. Despite the detailed knowledge about the many components of this cellular structure, it is not yet fully understood how these components self-organize spatiotemporally to form lamellipodia. We review here recent theoretical works where we have demonstrated that membrane-bound protein complexes that have intrinsic curvature and recruit the protrusive forces of the cytoskeleton result in a simple, yet highly robust, organizing feedback mechanism that organizes the cytoskeleton and the membrane. This self-organization mechanism accounts for the formation of flat lamellipodia at the leading edge of cells spreading over adhesive substrates, allowing for the emergence of a polarized, motile 'minimal cell' model. The same mechanism describes how lamellipodia organize to drive robust engulfment of particles during phagocytosis and explains in simple physical terms the spreading and migration of cells over fibers and other curved surfaces. This Review highlights that despite the complexity of cellular composition, there might be simple general physical principles that are utilized by the cell to drive cellular shape dynamics.

Geometry-induced patterns through mechanochemical coupling

Intracellular protein patterns regulate a variety of vital cellular processes such as cell division and motility, which often involve dynamic cell-shape changes. These changes in cell shape may in turn affect the dynamics of pattern-forming proteins, hence leading to an intricate feedback loop between cell shape and chemical dynamics. While several computational studies have examined the rich resulting dynamics, the underlying mechanisms are not yet fully understood. To elucidate some of these mechanisms, we explore a conceptual model for cell polarity on a dynamic one-dimensional manifold. Using concepts from differential geometry, we derive the equations governing mass-conserving reaction-diffusion systems on time-evolving manifolds. Analyzing these equations mathematically, we show that dynamic shape changes of the membrane can induce pattern-forming instabilities in parts of the membrane, which we refer to as regional instabilities. Deformations of the local membrane geometry can also (regionally) suppress pattern formation and spatially shift already existing patterns. We explain our findings by applying and generalizing the local equilibria theory of mass-conserving reaction-diffusion systems. This allows us to determine a simple onset criterion for geometry-induced pattern-forming instabilities, which is linked to the phase-space structure of the reaction-diffusion system. The feedback loop between membrane shape deformations and reaction-diffusion dynamics then leads to a surprisingly rich phenomenology of patterns, including oscillations, traveling waves, and standing waves, even if these patterns do not occur in systems with a fixed membrane shape. Our paper reveals that the local conformation of the membrane geometry acts as an important dynamical control parameter for pattern formation in mass-conserving reaction-diffusion systems.

Modelling how curved active proteins and shear flow pattern cellular shape and motility

Cell spreading and motility on an adhesive substrate are driven by the active physical forces generated by the actin cytoskeleton. We have recently shown that coupling curved membrane complexes to protrusive forces, exerted by the actin polymerization that they recruit, provides a mechanism that can give rise to spontaneous membrane shapes and patterns. In the presence of an adhesive substrate, this model was shown to give rise to an emergent motile phenotype, resembling a motile cell. Here, we utilize this "minimal-cell" model to explore the impact of external shear flow on the cell shape and migration on a uniform adhesive flat substrate. We find that in the presence of shear the motile cell reorients such that its leading edge, where the curved active proteins aggregate, faces the shear flow. The flow-facing configuration is found to minimize the adhesion energy by allowing the cell to spread more efficiently over the substrate. For the non-motile vesicle shapes, we find that they mostly slide and roll with the shear flow. We compare these theoretical results with experimental observations, and suggest that the tendency of many cell types to move against the flow may arise from the very general, and non-cell-type-specific mechanism predicted by our model.

Theoretical model of membrane protrusions driven by curved active proteins

Eukaryotic cells intrinsically change their shape, by changing the composition of their membrane and by restructuring their underlying cytoskeleton. We present here further studies and extensions of a minimal physical model, describing a closed vesicle with mobile curved membrane protein complexes. The cytoskeletal forces describe the protrusive force due to actin polymerization which is recruited to the membrane by the curved protein complexes. We characterize the phase diagrams of this model, as function of the magnitude of the active forces, nearest-neighbor protein interactions and the proteins' spontaneous curvature. It was previously shown that this model can explain the formation of lamellipodia-like flat protrusions, and here we explore the regimes where the model can also give rise to filopodia-like tubular protrusions. We extend the simulation with curved components of both convex and concave species, where we find the formation of complex ruffled clusters, as well as internalized invaginations that resemble the process of endocytosis and macropinocytosis. We alter the force model representing the cytoskeleton to simulate the effects of bundled instead of branched structure, resulting in shapes which resemble filopodia.

Actin Filaments Couple the Protrusive Tips to the Nucleus through the I-BAR Domain Protein IRSp53 during the Migration of Cells on 1D Fibers

The cell migration cycle, well-established in 2D, proceeds with forming new protrusive structures at the cell membrane and subsequent redistribution of contractile machinery. Three-dimensional (3D) environments are complex and composed of 1D fibers, and 1D fibers are shown to recapitulate essential features of 3D migration. However, the establishment of protrusive activity at the cell membrane and contractility in 1D fibrous environments remains partially understood. Here the role of membrane curvature regulator IRSp53 is examined as a coupler between actin filaments and plasma membrane during cell migration on single, suspended 1D fibers. IRSp53 depletion reduced cell-length spanning actin stress fibers that originate from the cell periphery, protrusive activity, and contractility, leading to uncoupling of the nucleus from cellular movements. A theoretical model capable of predicting the observed transition of IRSp53-depleted cells from rapid stick-slip migration to smooth and slower migration due to reduced actin polymerization at the cell edges is developed, which is verified by direct measurements of retrograde actin flow using speckle microscopy. Overall, it is found that IRSp53 mediates actin recruitment at the cellular tips leading to the establishment of cell-length spanning fibers, thus demonstrating a unique role of IRSp53 in controlling cell migration in 3D.

Modulation of self-organizing circuits at deforming membranes by intracellular and extracellular factors

Mechanical forces exerted to the plasma membrane induce cell shape changes. These transient shape changes trigger, among others, enrichment of curvature-sensitive molecules at deforming membrane sites. Strikingly, some curvature-sensing molecules not only detect membrane deformation but can also alter the amplitude of forces that caused to shape changes in the first place. This dual ability of sensing and inducing membrane deformation leads to the formation of curvature-dependent self-organizing signaling circuits. How these cell-autonomous circuits are affected by auxiliary parameters from inside and outside of the cell has remained largely elusive. Here, we explore how such factors modulate self-organization at the micro-scale and its emerging properties at the macroscale.

A theoretical model of efficient phagocytosis driven by curved membrane proteins and active cytoskeleton forces

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.

Remeshing flexible membranes under the control of free energy

Cell membranes are flexible and often undergo large-scale morphological changes during processes like mitosis, protrusion and retraction, or vesicle fusion. Mathematical modeling of cell membranes depends on a representation of the free-form surface by discrete meshes. During morphological changes, these meshes must be adjusted under the minimization of the total free energy. Current methodology for meshing is limited in one of two ways: 1) Free energy-dependent methods have no restriction on the mesh geometry. The resulting irregular meshes cause artifacts in follow-up models of morphodynamics. 2) Geometry-dependent methods maintain mesh quality but violate the physics of free energy minimization. To fill this gap, we regulate mesh geometries via a free-energy-determined remeshing process: adding and removing mesh elements upon morphological changes based on barrier crossings in a double-barrier potential between neighboring vertices in the meshes. We test the method's robustness by reproducing the morphodynamics of red blood cells and vesicle fusions; and we demonstrate the method's adaptability by simulating the formation of filopodia, lamellipodia and invaginations. Finally, we use the method to study a mechanical decoupling effect of two connected membrane tethers that has been recently observed experimentally, but has not been mechanistically explained in the context of a complete membrane surface. We propose a biophysical model that strengthens the decoupling effect and broadens the original interpretation of the experiment. The method is developed in C/Matlab and distributed via https://github.com/DanuserLab/biophysicsModels.

Universal properties of active membranes

We put forward a general field theory for nearly flat fluid membranes with embedded activators and analyze their critical properties using renormalization group techniques. Depending on the membrane-activator coupling, we find a crossover between acoustic and diffusive scaling regimes, with mean-field dynamical critical exponents z=1 and 2, respectively. We argue that the acoustic scaling, which is exact in all spatial dimensions, leads to an early-time behavior, which is representative of the spatiotemporal patterns observed at the leading edge of motile cells, such as oscillations superposed on the growth of the membrane width. In the case of mean-field diffusive scaling, one-loop corrections to the mean-field exponents reveal universal behavior distinct from the Kardar-Parisi-Zhang scaling of passive interfaces and signs of strong-coupling behavior.

Renormalization group study of the dynamics of active membranes: Universality classes and scaling laws

Motivated by experimental observations of patterning at the leading edge of motile eukaryotic cells, we introduce a general model for the dynamics of nearly-flat fluid membranes driven from within by an ensemble of activators. We include, in particular, a kinematic coupling between activator density and membrane slope which generically arises whenever the membrane has a nonvanishing normal speed. We unveil the phase diagram of the model by means of a perturbative field-theoretical renormalization group analysis. Due to the aforementioned kinematic coupling the natural early-time dynamical scaling is acoustic, that is the dynamical critical exponent is 1. However, as soon as the the normal velocity of the membrane is tuned to zero, the system crosses over to diffusive dynamic scaling in mean field. Distinct critical points can be reached depending on how the limit of vanishing velocity is realized: in each of them corrections to scaling due to nonlinear coupling terms must be taken into account. The detailed analysis of these critical points reveals novel scaling regimes which can be accessed with perturbative methods, together with signs of strong coupling behavior, which establishes a promising ground for further nonperturbative calculations. Our results unify several previous studies on the dynamics of active membrane, while also identifying nontrivial scaling regimes which cannot be captured by passive theories of fluctuating interfaces and are relevant for the physics of living membranes.

Vesicle shape transformations driven by confined active filaments

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.

Dendrite tapering actuates a self-organizing signaling circuit for stochastic filopodia initiation in neurons

How signaling units spontaneously arise from a noisy cellular background is not well understood. Here, we show that stochastic membrane deformations can nucleate exploratory dendritic filopodia, dynamic actin-rich structures used by neurons to sample its surroundings for compatible transcellular contacts. A theoretical analysis demonstrates that corecruitment of positive and negative curvature-sensitive proteins to deformed membranes minimizes the free energy of the system, allowing the formation of long-lived curved membrane sections from stochastic membrane fluctuations. Quantitative experiments show that once recruited, curvature-sensitive proteins form a signaling circuit composed of interlinked positive and negative actin-regulatory feedback loops. As the positive but not the negative feedback loop can sense the dendrite diameter, this self-organizing circuit determines filopodia initiation frequency along tapering dendrites. Together, our findings identify a receptor-independent signaling circuit that employs random membrane deformations to simultaneously elicit and limit formation of exploratory filopodia to distal dendritic sites of developing neurons.

On the Role of Electrostatic Repulsion in Topological Defect-Driven Membrane Fission

Within a modified Langevin Poisson–Boltzmann model of electric double layers, we derived an analytical expression for osmotic pressure between two charged surfaces. The orientational ordering of the water dipoles as well as the space dependencies of electric potentials, electric fields, and osmotic pressure between two charged spheres were taken into account in the model. Thus, we were able to capture the interaction between the parent cell and connected daughter vesicle or the interactions between neighbouring beads in necklace-like membrane protrusions. The predicted repulsion between them can facilitate the topological antidefect-driven fission of membrane daughter vesicles and the fission of beads of undulated membrane protrusions.

A Monte Carlo study of giant vesicle morphologies in nonequilibrium environments

It is known that giant vesicles undergo dynamic morphological changes when exposed to a detergent. The solubilization process may take multiple pathways. In this work, we identify lipid vesicle shape dynamics before the solubilization of 1,2-dioleoyl-sn-glycero-3-phosphocholine giant vesicles with Triton X-100 (TR) detergent. The violent lipid vesicle dynamics was observed with laser confocal scanning microscopy and was qualitatively explained via a numerical simulation. A three-dimensional Monte Carlo scheme was constructed that emulated the nonequilibrium conditions at the beginning stages of solubilization, accounting for a gradual addition of TR detergent molecules into the lipid bilayers. We suggest that the main driving factor for morphology change in lipid vesicles is the associative tendency of the TR molecules, which induces spontaneous curvature of the detergent inclusions, an intrinsic consequence of their molecular shape. The majority of the observed lipid vesicle shapes in the experiments were found to correspond very well to the numerically calculated shapes in the phase space of possible solutions. The results give an insight into the early stages of lipid vesicle solubilization by amphiphilic molecules, which is nonequilibrium in nature and very difficult to study.

Membrane bending by protein phase separation

Membrane bending is a ubiquitous cellular process that is required for membrane traffic, cell motility, organelle biogenesis, and cell division. Proteins that bind to membranes using specific structural features, such as wedge-like amphipathic helices and crescent-shaped scaffolds, are thought to be the primary drivers of membrane bending. However, many membrane-binding proteins have substantial regions of intrinsic disorder which lack a stable three-dimensional structure. Interestingly, many of these disordered domains have recently been found to form networks stabilized by weak, multivalent contacts, leading to assembly of protein liquid phases on membrane surfaces. Here we ask how membrane-associated protein liquids impact membrane curvature. We find that protein phase separation on the surfaces of synthetic and cell-derived membrane vesicles creates a substantial compressive stress in the plane of the membrane. This stress drives the membrane to bend inward, creating protein-lined membrane tubules. A simple mechanical model of this process accurately predicts the experimentally measured relationship between the rigidity of the membrane and the diameter of the membrane tubules. Discovery of this mechanism, which may be relevant to a broad range of cellular protrusions, illustrates that membrane remodeling is not exclusive to structured scaffolds but can also be driven by the rapidly emerging class of liquid-like protein networks that assemble at membranes.

Mechanics of cup-shaped caveolae

Caveolae are cell membrane invaginations of defined lipid and protein composition that flatten with increasing membrane tension. Super-resolution light microscopy and electron microscopy have revealed that caveolae can take a variety of cuplike shapes. We show here that, for the range in membrane tension relevant for cell membranes, the competition between membrane tension and membrane bending yields caveolae with cuplike shapes similar to those observed experimentally. We find that the caveola shape and its sensitivity to changes in membrane tension can depend strongly on the caveola spontaneous curvature and on the size of caveola domains. Our results suggest that heterogeneity in caveola shape produces a staggered response of caveolae to mechanical perturbations of the cell membrane, which may facilitate regulation of membrane tension over the wide range of scales thought to be relevant for cell membranes.

Mechanical and Electrical Interaction of Biological Membranes with Nanoparticles and Nanostructured Surfaces

In this review paper, we theoretically explain the origin of electrostatic interactions between lipid bilayers and charged solid surfaces using a statistical mechanics approach, where the orientational degree of freedom of lipid head groups and the orientational ordering of the water dipoles are considered. Within the modified Langevin Poisson–Boltzmann model of an electric double layer, we derived an analytical expression for the osmotic pressure between the planar zwitterionic lipid bilayer and charged solid planar surface. We also show that the electrostatic interaction between the zwitterionic lipid head groups of the proximal leaflet and the negatively charged solid surface is accompanied with a more perpendicular average orientation of the lipid head-groups. We further highlight the important role of the surfaces’ nanostructured topography in their interactions with biological material. As an example of nanostructured surfaces, we describe the synthesis of TiO₂ nanotubular and octahedral surfaces by using the electrochemical anodization method and hydrothermal method, respectively. The physical and chemical properties of these nanostructured surfaces are described in order to elucidate the influence of the surface topography and other physical properties on the behavior of human cells adhered to TiO₂ nanostructured surfaces. In the last part of the paper, we theoretically explain the interplay of elastic and adhesive contributions to the adsorption of lipid vesicles on the solid surfaces. We show the numerically predicted shapes of adhered lipid vesicles corresponding to the minimum of the membrane free energy to describe the influence of the vesicle size, bending modulus, and adhesion strength on the adhesion of lipid vesicles on solid charged surfaces.

On the Role of Curved Membrane Nanodomains, and Passive and Active Skeleton Forces in the Determination of Cell Shape and Membrane Budding

Biological membranes are composed of isotropic and anisotropic curved nanodomains. Anisotropic membrane components, such as Bin/Amphiphysin/Rvs (BAR) superfamily protein domains, could trigger/facilitate the growth of membrane tubular protrusions, while isotropic curved nanodomains may induce undulated (necklace-like) membrane protrusions. We review the role of isotropic and anisotropic membrane nanodomains in stability of tubular and undulated membrane structures generated or stabilized by cyto- or membrane-skeleton. We also describe the theory of spontaneous self-assembly of isotropic curved membrane nanodomains and derive the critical concentration above which the spontaneous necklace-like membrane protrusion growth is favorable. We show that the actin cytoskeleton growth inside the vesicle or cell can change its equilibrium shape, induce higher degree of segregation of membrane nanodomains or even alter the average orientation angle of anisotropic nanodomains such as BAR domains. These effects may indicate whether the actin cytoskeleton role is only to stabilize membrane protrusions or to generate them by stretching the vesicle membrane. Furthermore, we demonstrate that by taking into account the in-plane orientational ordering of anisotropic membrane nanodomains, direct interactions between them and the extrinsic (deviatoric) curvature elasticity, it is possible to explain the experimentally observed stability of oblate (discocyte) shapes of red blood cells in a broad interval of cell reduced volume. Finally, we present results of numerical calculations and Monte-Carlo simulations which indicate that the active forces of membrane skeleton and cytoskeleton applied to plasma membrane may considerably influence cell shape and membrane budding.

Electric Double Layer and Orientational Ordering of Water Dipoles in Narrow Channels within a Modified Langevin Poisson-Boltzmann Model

The electric double layer (EDL) is an important phenomenon that arises in systems where a charged surface comes into contact with an electrolyte solution. In this work we describe the generalization of classic Poisson-Boltzmann (PB) theory for point-like ions by taking into account orientational ordering of water molecules. The modified Langevin Poisson-Boltzmann (LPB) model of EDL is derived by minimizing the corresponding Helmholtz free energy functional, which includes also orientational entropy contribution of water dipoles. The formation of EDL is important in many artificial and biological systems bound by a cylindrical geometry. We therefore numerically solve the modified LPB equation in cylindrical coordinates, determining the spatial dependencies of electric potential, relative permittivity and average orientations of water dipoles within charged tubes of different radii. Results show that for tubes of a large radius, macroscopic (net) volume charge density of coions and counterions is zero at the geometrical axis. This is attributed to effective electrolyte charge screening in the vicinity of the inner charged surface of the tube. For tubes of small radii, the screening region extends into the whole inner space of the tube, leading to non-zero net volume charge density and non-zero orientational ordering of water dipoles near the axis.

Non-uniform distribution of myosin-mediated forces governs red blood cell membrane curvature through tension modulation

The biconcave disk shape of the mammalian red blood cell (RBC) is unique to the RBC and is vital for its circulatory function. Due to the absence of a transcellular cytoskeleton, RBC shape is determined by the membrane skeleton, a network of actin filaments cross-linked by spectrin and attached to membrane proteins. While the physical properties of a uniformly distributed actin network interacting with the lipid bilayer membrane have been assumed to control RBC shape, recent experiments reveal that RBC biconcave shape also depends on the contractile activity of nonmuscle myosin IIA (NMIIA) motor proteins. Here, we use the classical Helfrich-Canham model for the RBC membrane to test the role of heterogeneous force distributions along the membrane and mimic the contractile activity of sparsely distributed NMIIA filaments. By incorporating this additional contribution to the Helfrich-Canham energy, we find that the RBC biconcave shape depends on the ratio of forces per unit volume in the dimple and rim regions of the RBC. Experimental measurements of NMIIA densities at the dimple and rim validate our prediction that (a) membrane forces must be non-uniform along the RBC membrane and (b) the force density must be larger in the dimple than the rim to produce the observed membrane curvatures. Furthermore, we predict that RBC membrane tension and the orientation of the applied forces play important roles in regulating this force-shape landscape. Our findings of heterogeneous force distributions on the plasma membrane for RBC shape maintenance may also have implications for shape maintenance in different cell types.

The spectrin-actin network of the membrane skeleton plays an important role in controlling specialized cell membrane morphology. In the paradigmatic red blood cell (RBC), where actin filaments are present exclusively in the membrane skeleton, recent experiments reveal that nonmuscle myosin IIA (NMIIA) motor contractility maintains the RBC biconcave disk shape. In this study, we have identified criteria for micron-scale distributions of NMIIA forces at the membrane required to maintain the biconcave disk shape of an RBC in the resting condition. Supported by experimental measurements of RBC NMIIA distribution, we showed that a heterogeneous force distribution with a larger force density at the dimple is able to capture the experimentally observed biconcave morphology of an RBC with better accuracy compared to previous models that did not consider the heterogeneity in the force distribution. Furthermore, we showed that the biconcave geometry of the RBC is closely regulated by the effective membrane tension and the direction of applied forces on the membrane. These findings can be generalized to any force-mediated membrane shape, providing insight into the role of actomyosin forces in prescribing and maintaining the morphology of different cell types.

Active Contact Forces Drive Nonequilibrium Fluctuations in Membrane Vesicles

We analyze the nonequilibrium shape fluctuations of giant unilamellar vesicles encapsulating motile bacteria. Owing to bacteria-membrane collisions, we experimentally observe a significant increase in the magnitude of membrane fluctuations at low wave numbers, compared to the well-known thermal fluctuation spectrum. We interrogate these results by numerically simulating membrane height fluctuations via a modified Langevin equation, which includes bacteria-membrane contact forces. Taking advantage of the lengthscale and timescale separation of these contact forces and thermal noise, we further corroborate our results with an approximate theoretical solution to the dynamical membrane equations. Our theory and simulations demonstrate excellent agreement with nonequilibrium fluctuations observed in experiments. Moreover, our theory reveals that the fluctuation-dissipation theorem is not broken by the bacteria; rather, membrane fluctuations can be decomposed into thermal and active components.

Cell-Substrate Patterns Driven by Curvature-Sensitive Actin Polymerization: Waves and Podosomes

Cells adhered to an external solid substrate are observed to exhibit rich dynamics of actin structures on the basal membrane, which are distinct from those observed on the dorsal (free) membrane. Here we explore the dynamics of curved membrane proteins, or protein complexes, that recruit actin polymerization when the membrane is confined by the solid substrate. Such curved proteins can induce the spontaneous formation of membrane protrusions on the dorsal side of cells. However, on the basal side of the cells, such protrusions can only extend as far as the solid substrate and this constraint can convert such protrusions into propagating wave-like structures. We also demonstrate that adhesion molecules can stabilize localized protrusions that resemble some features of podosomes. This coupling of curvature and actin forces may underlie the differences in the observed actin-membrane dynamics between the basal and dorsal sides of adhered cells.

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.

Numerical Model for the Determination of Erythrocyte Mechanical Properties and Wall Shear Stress in vivo From Intravital Microscopy

The mechanical properties and deformability of Red Blood Cells (RBCs) are important determinants of blood rheology and microvascular hemodynamics. The objective of this study is to quantify the mechanical properties and wall shear stress experienced by the RBC membrane during capillary plug flow in vivo utilizing high speed video recording from intravital microscopy, biomechanical modeling, and computational methods. Capillaries were imaged in the rat cremaster muscle pre- and post-RBC transfusion of stored RBCs for 2-weeks. RBC membrane contours were extracted utilizing image processing and parametrized. RBC parameterizations were used to determine updated deformation gradient and Lagrangian Green strain tensors for each point along the parametrization and for each frame during plug flow. The updated Lagrangian Green strain and Displacement Gradient tensors were numerically fit to the Navier-Lame equations along the parameterized boundary to determined Lame's constants. Mechanical properties and wall shear stress were determined before and transfusion, were grouped in three populations of erythrocytes: native cells (NC) or circulating cells before transfusion, and two distinct population of cells after transfusion with stored cells (SC1 and SC2). The distinction, between the heterogeneous populations of cells present after the transfusion, SC1 and SC2, was obtained through principle component analysis (PCA) of the mechanical properties along the membrane. Cells with the first two principle components within 3 standard deviations of the mean, were labeled as SC1, and those with the first two principle components greater than 3 standard deviations from the mean were labeled as SC2. The calculated shear modulus average was 1.1±0.2, 0.90±0.15, and 12 ± 8 MPa for NC, SC1, and SC2, respectively. The calculated young's modulus average was 3.3±0.6, 2.6±0.4, and 32±20 MPa for NC, SC1, and SC2, respectively. o our knowledge, the methods presented here are the first estimation of the erythrocyte mechanical properties and shear stress in vivo during capillary plug flow. In summary, the methods introduced in this study may provide a new avenue of investigation of erythrocyte mechanics in the context of hematologic conditions that adversely affect erythrocyte mechanical properties.

Cell confinement reveals a branched-actin independent circuit for neutrophil polarity

Migratory cells use distinct motility modes to navigate different microenvironments, but it is unclear whether these modes rely on the same core set of polarity components. To investigate this, we disrupted actin-related protein 2/3 (Arp2/3) and the WASP-family verprolin homologous protein (WAVE) complex, which assemble branched actin networks that are essential for neutrophil polarity and motility in standard adherent conditions. Surprisingly, confinement rescues polarity and movement of neutrophils lacking these components, revealing a processive bleb-based protrusion program that is mechanistically distinct from the branched actin-based protrusion program but shares some of the same core components and underlying molecular logic. We further find that the restriction of protrusion growth to one site does not always respond to membrane tension directly, as previously thought, but may rely on closely linked properties such as local membrane curvature. Our work reveals a hidden circuit for neutrophil polarity and indicates that cells have distinct molecular mechanisms for polarization that dominate in different microenvironments.

Cells display a high degree of plasticity in migration, but how polarity is organized in different microenvironments has remained unclear. This study uses mechanical perturbations to reveal that migration using actin-rich or bleb-based protrusions are both organized around Rac GTPase.

Inception Mechanisms of Tunneling Nanotubes

Tunneling nanotubes (TNTs) are thin membranous tubes that interconnect cells, representing a novel route of cell-to-cell communication and spreading of pathogens. TNTs form between many cell types, yet their inception mechanisms remain elusive. We review in this study general concepts related to the formation and stability of membranous tubular structures with a focus on a deviatoric elasticity model of membrane nanodomains. We review experimental evidence that tubular structures initiate from local membrane bending facilitated by laterally distributed proteins or anisotropic membrane nanodomains. We further discuss the numerical results of several theoretical and simulation models of nanodomain segregation suggesting the mechanisms of TNT inception and stability. We discuss the coupling of nanodomain segregation with the action of protruding cytoskeletal forces, which are mostly provided in eukaryotic cells by the polymerization of f-actin, and review recent inception mechanisms of TNTs in relation to motor proteins.