Bending elasticity of vesicle membranes studied by Monte Carlo simulations of vesicle thermal shape fluctuations

There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning

We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.

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.

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.

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.

Polymer-Lipid Hybrid Materials

Hierarchic self-assembly underpins much of the form and function seen in synthetic or biological soft materials. Lipids are paramount examples, building themselves in nature or synthetically in a variety of meso/nanostructures. Synthetic block copolymers capture many of lipid's structural and functional properties. Lipids are typically biocompatible and high molecular weight polymers are mechanically robust and chemically versatile. The development of new materials for applications like controlled drug/gene/protein delivery, biosensors, and artificial cells often requires the combination of lipids and polymers. The emergent composite material, a "polymer-lipid hybrid membrane", displays synergistic properties not seen in pure components. Specific examples include the observation that hybrid membranes undergo lateral phase separation that can correlate in registry across multiple layers into a three-dimensional phase-separated system with enhanced permeability of encapsulated drugs. It is timely to underpin these emergent properties in several categories of hybrid systems ranging from colloidal suspensions to supported hybrid films. In this review, we discuss the form and function of a vast number of polymer-lipid hybrid systems published to date. We rationalize the results to raise new fundamental understanding of hybrid self-assembling soft materials as well as to enable the design of new supramolecular systems and applications.

The Elucidation of the Molecular Mechanism of the Extrusion Process

Extrusion is a popular method for producing homogenous population of unilamellar liposomes. The technique relies on forcing a lipid suspension through cylindrical pores in a polycarbonate membrane. The quantification of the extrusion and/or recalibration processes make possible the acquisition of experimental data, which can be correlated with the mechanical properties of the lipid bilayer. In this work, the force needed for the extrusion process was correlated with the mechanical properties of a lipid bilayer derived from other experiments. Measurements were performed using a home-made dedicated device capable of maintaining a stable volumetric flux of a liposome suspension through well-defined pores and to continuously measure the extrusion force. Based on the obtained results, the correlation between the lipid bilayer bending rigidity and extrusion force was derived. Specifically, it was found that the bending rigidity of liposomes formed from well-defined lipid mixtures agrees with data obtained by others using flicker-noise spectroscopy or micromanipulation. The other issue addressed in the presented studies was the identification of molecular mechanisms leading to the formation of unilamellar vesicles in the extrusion process. Finally, it was demonstrated that during the extrusion, lipids are not exchanged between vesicles, i.e., vesicles can divide but no membrane fusion or lipid exchange between bilayers was detected.

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.

Domain formation in bicomponent vesicles induced by composition-curvature coupling

Lipid vesicles composed of a mixture of two types of lipids are studied by intensive Monte Carlo numerical simulations. The coupling between the local composition and the membrane shape is induced by two different spontaneous curvatures of the components. We explore the various morphologies of these biphasic vesicles coupled to the observed patterns such as nano-domains or labyrinthine mesophases. The effect of the difference in curvatures, the surface tension, and the interaction parameter between components is thoroughly explored. Our numerical results quantitatively agree with the previous analytical results obtained by Gueguen et al. [Eur. Phys. J. E 37, 76 (2014)] in the disordered (high temperature) phase. Numerical simulations allow us to explore the full parameter space, especially close to and below the critical temperature, where analytical results are not accessible. Phase diagrams are constructed and domain morphologies are quantitatively studied by computing the structure factor and the domain size distribution. This mechanism likely explains the existence of nano-domains in cell membranes as observed by super-resolution fluorescence microscopy.

Magneto-mechanical actuation of barium-hexaferrite nanoplatelets for the disruption of phospholipid membranes

HYPOTHESIS

The magneto-mechanical actuation (MMA) of magnetic nanoparticles with a low-frequency alternating magnetic field (AMF) can be used to destroy cancer cells. So far, MMA was tested on different cells using different nanoparticles and different field characteristics, which makes comparisons and any generalizations about the results of MMA difficult. In this paper we propose the use of giant unilamellar vesicles (GUVs) as a simple model system to study the effect of MMA on a closed lipid bilayer membrane, i.e., a basic building block of any cell.

EXPERIMENTS

The GUVs were exposed to barium-hexaferrite nanoplatelets (NPLs, ~50 nm wide and 3 nm thick) with unique magnetic properties dominated by a permanent magnetic moment that is perpendicular to the platelet, at different concentrations (1-50 µg/mL) and pH values (4.2-7.4) of the aqueous suspension. The GUVs were observed with an optical microscope while being exposed to a uniaxial AMF (3-100 Hz, 2.2-10.6 mT).

FINDINGS

When the NPLs were electrostatically attached to the GUV membranes, the MMA induced cyclic fluctuations of the GUVs' shape corresponding to the AMF frequency at the low NPL concentration (1 µm/mL), whereas the GUVs were bursting at the higher concentration (10 µg/mL). Theoretical considerations suggested that the bursting of the GUVs is a consequence of the local action of an assembly of several NPLs, rather than a collective effect of all the absorbed NPLs.

Normal red blood cells' shape stabilized by membrane's in-plane ordering

Red blood cells (RBCs) are present in almost all vertebrates and their main function is to transport oxygen to the body tissues. RBCs' shape plays a significant role in their functionality. In almost all mammals in normal conditions, RBCs adopt a disk-like (discocyte) shape, which optimizes their flow properties in vessels and capillaries. Experimentally measured values of the reduced volume (v) of stable discocyte shapes range in a relatively broad window between v ~ 0.58 and 0.8. However, these observations are not supported by existing theoretical membrane-shape models, which predict that discocytic RBC shape is stable only in a very narrow interval of v values, ranging between v ~ 0.59 and 0.65. In this study, we demonstrate that this interval is broadened if a membrane's in-plane ordering is taken into account. We model RBC structures by using a hybrid Helfrich-Landau mesoscopic approach. We show that an extrinsic (deviatoric) curvature free energy term stabilizes the RBC discocyte shapes. In particular, we show on symmetry grounds that the role of extrinsic curvature is anomalously increased just below the nematic in-plane order-disorder phase transition temperature.

Membrane stretching elasticity and thermal shape fluctuations of nearly spherical lipid vesicles

One of the most widely used methods for determination of the bending elasticity modulus of model lipid membranes is the analysis of the shape fluctuations of nearly spherical lipid vesicles. The theoretical basis of this analysis is given by Milner and Safran [Phys. Rev. A 36, 4371 (1987)0556-279110.1103/PhysRevA.36.4371]. In their theory the stretching effects are not considered. In the present study we generalized their approach including the stretching effects deduced after application of the statistical mechanics to vesicles.

Fluctuation tension and shape transition of vesicles: renormalisation calculations and Monte Carlo simulations

It has been known for long that the fluctuation surface tension of membranes r, computed from the height fluctuation spectrum, is not equal to the bare surface tension σ, which is introduced in the theory either as a Lagrange multiplier to conserve the total membrane area or as an external constraint. In this work we relate these two surface tensions both analytically and numerically. They are also compared to the Laplace tension γ, and the mechanical frame tension τ. Using the Helfrich model and one-loop renormalisation calculations, we obtain, in addition to the effective bending modulus κ_eff, a new expression for the effective surface tension σ_eff = σ - εk_BT/(2a_p) where k_BT is the thermal energy, a_p the projected cut-off area, and ε = 3 or 1 according to the allowed configurations that keep either the projected area or the total area constant. Moreover we show that the crumpling transition for an infinite planar membrane occurs for σ_eff = 0, and also that it coincides with vanishing Laplace and frame tensions. Using extensive Monte Carlo (MC) simulations, triangulated membranes of vesicles made of N = 100-2500 vertices are simulated within the Helfrich theory. As compared to alternative numerical models, no local constraint is applied and the shape is only controlled by the constant volume, the spontaneous curvature and σ. It is shown that the numerical fluctuation surface tension r is equal to σ_eff both with radial MC moves (ε = 3) and with corrected MC moves locally normal to the fluctuating membrane (ε = 1). For finite vesicles of typical size R, two different regimes are defined: a tension regime for [small sigma, Greek, circumflex]_eff = σ_effR²/κ_eff > 0 and a bending one for -1 < [small sigma, Greek, circumflex]_eff < 0. A shape transition from a quasi-spherical shape imposed by the large surface energy, to more deformed shapes only controlled by the bending energy, is observed numerically at [small sigma, Greek, circumflex]_eff ≃ 0. We propose that the buckling transition, observed for planar supported membranes in the literature, occurs for [small sigma, Greek, circumflex]_eff ≃ -1, the associated negative frame tension playing the role of a compressive force. Hence, a precise control of the value of σ_eff in simulations cannot but enhance our understanding of shape transitions of vesicles and cells.

On the role of external force of actin filaments in the formation of tubular protrusions of closed membrane shapes with anisotropic membrane components

Biological membranes are composed of different components and there is no a priori reason to assume that all components are isotropic. It was previously shown that the anisotropic properties of membrane components may explain the stability of membrane tubular protrusions even without the application of external force. Our theoretical study focuses on the role of anisotropic membrane components in the stability of membrane tubular structures generated or stabilized by actin filaments. We show that the growth of the actin cytoskeleton inside the vesicle can induce the partial lateral segregation of different membrane components. The entropy of mixing of membrane components hinders the total lateral segregation of the anisotropic and isotropic membrane components. Self-assembled aggregates formed by anisotropic membrane components facilitate the growth of long membrane tubular protrusions. Protrusive force generated by actin filaments favors strong segregation of membrane components by diminishing the opposing effect of mixing entropy.