Budding and vesiculation induced by conical membrane inclusions

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.

Membrane-Mediated Interactions Between Nonspherical Elastic Particles

The transport of particles across lipid-bilayer membranes is important for biological cells to exchange information and material with their environment. Large particles often get wrapped by membranes, a process which has been intensively investigated in the case of hard particles. However, many particles in vivo and in vitro are deformable, e.g., vesicles, filamentous viruses, macromolecular condensates, polymer-grafted nanoparticles, and microgels. Vesicles may serve as a generic model system for deformable particles. Here, we study nonspherical vesicles with various sizes, shapes, and elastic properties at initially planar lipid-bilayer membranes. Using the Helfrich Hamiltonian, triangulated membranes, and energy minimization, we predict the interplay of vesicle shapes and wrapping states. Increasing particle softness enhances the stability of shallow-wrapped and deep-wrapped states over nonwrapped and complete-wrapped states. The free membrane mediates an interaction between partial-wrapped vesicles. For the pair interaction between deep-wrapped vesicles, we predict repulsion. For shallow-wrapped vesicles, we predict attraction for tip-to-tip orientation and repulsion for side-by-side orientation. Our predictions may guide the design and fabrication of deformable particles for efficient use in medical applications, such as targeted drug delivery.

Mechanics of diffusion-mediated budding and implications for virus replication and infection

Budding allows virus replication and macromolecular secretion in cells through the formation of a membrane protrusion (bud) that evolves into an envelope. The largest energetic barrier to bud formation is membrane deflection and is trespassed primarily thanks to nucleocapsid-membrane adhesion. Transmembrane proteins (TPs), which later form the virus ligands, are the main promotors of adhesion and can accommodate membrane bending thanks to an induced spontaneous curvature. Adhesive TPs must diffuse across the membrane from remote regions to gather on the bud surface, thus, diffusivity controls the kinetics. This paper proposes a simple model to describe diffusion-mediated budding unravelling important size limitations and size-dependent kinetics. The predicted optimal virion radius, giving the fastest budding, is validated against experiments for coronavirus, HIV, flu and hepatitis. Assuming exponential replication of virions and hereditary size, the model can predict the size distribution of a virus population. This is verified against experiments for SARS-CoV-2. All the above comparisons rely on the premise that budding poses the tightest size constraint. This is true in most cases, as demonstrated in this paper, where the proposed model is extended to describe virus infection via receptor- and clathrin-mediated endocytosis, and via membrane fusion.

Amphiphilic Gold Nanoparticles: A Biomimetic Tool to Gain Mechanistic Insights into Peptide-Lipid Interactions

Functional peptides are now widely used in a myriad of biomedical and clinical contexts, from cancer therapy and tumor targeting to the treatment of bacterial and viral infections. Underlying this diverse range of applications are the non-specific interactions that can occur between peptides and cell membranes, which, in many contexts, result in spontaneous internalization of the peptide within cells by avoiding energy-driven endocytosis. For this to occur, the amphipathicity and surface structural flexibility of the peptides play a crucial role and can be regulated by the presence of specific molecular residues that give rise to precise molecular events. Nevertheless, most of the mechanistic details regulating the encounter between peptides and the membranes of bacterial or animal cells are still poorly understood, thus greatly limiting the biomimetic potential of these therapeutic molecules. In this arena, finely engineered nanomaterials—such as small amphiphilic gold nanoparticles (AuNPs) protected by a mixed thiol monolayer—can provide a powerful tool for mimicking and investigating the physicochemical processes underlying peptide-lipid interactions. Within this perspective, we present here a critical review of membrane effects induced by both amphiphilic AuNPs and well-known amphiphilic peptide families, such as cell-penetrating peptides and antimicrobial peptides. Our discussion is focused particularly on the effects provoked on widely studied model cell membranes, such as supported lipid bilayers and lipid vesicles. Remarkable similarities in the peptide or nanoparticle membrane behavior are critically analyzed. Overall, our work provides an overview of the use of amphiphilic AuNPs as a highly promising tailor-made model to decipher the molecular events behind non-specific peptide-lipid interactions and highlights the main affinities observed both theoretically and experimentally. The knowledge resulting from this biomimetic approach could pave the way for the design of synthetic peptides with tailored functionalities for next-generation biomedical applications, such as highly efficient intracellular delivery systems.

The morphological role of ligand inhibitors in blocking receptor- and clathrin-mediated endocytosis

Cells often internalize particles through endocytic pathways that involve the binding between cell receptors and particle ligands, which drives the cell membrane to wrap the particle into a delivery vesicle. Previous findings showed that receptor-mediated endocytosis is impossible for spherical particles smaller than a minimum size because of the energy barrier created by membrane bending. In this study, we investigate the morphological role of ligand inhibitors in blocking endocytosis, inspired by antibodies that inhibit virus ligands to prevent infection. While ligand inhibitors have the obvious effect of reducing the driving force due to adhesion, they also have a nontrivial (morphological) impact on the entropic and elastic energy of the system. We determine the necessary conditions for endocytosis by considering the additional energy barrier due to the membrane bending to wrap the inhibiting protrusions. We find that inhibitors increase the minimum radius previously reported, depending on their density and size. In addition, we extend this result to the case of clathrin-mediated endocytosis, which is the most common pathway for virus entry. The assembly of a clathrin coat with a spontaneous curvature increases the energy barrier and sets a maximum particle size (in agreement with experimental observations on spherical particles). Our investigation suggests that morphological considerations can inform the optimal design of neutralizing viral antibodies and new strategies for targeted nanomedicine.

Mechanochemical coupling of lipid organization and protein function through membrane thickness deformations

Cell membranes are composed of a great variety of protein and lipid species with distinct unperturbed hydrophobic thicknesses. To achieve hydrophobic matching, the lipid bilayer tends to deform around membrane proteins so as to match the protein hydrophobic thickness at bilayer-protein interfaces. Such protein-induced distortions of the lipid bilayer hydrophobic thickness incur a substantial energy cost that depends critically on the bilayer-protein hydrophobic mismatch, while distinct conformational states of membrane proteins often show distinct hydrophobic thicknesses. As a result, hydrophobic interactions between membrane proteins and lipids can yield a rich interplay of lipid-protein organization and transitions in protein conformational state. We combine here the membrane elasticity theory of protein-induced lipid bilayer thickness deformations with the Landau-Ginzburg theory of lipid domain formation to systematically explore the coupling between local lipid organization, lipid and protein hydrophobic thickness, and protein-induced lipid bilayer thickness deformations in membranes with heterogeneous lipid composition. We allow for a purely mechanical coupling of lipid and protein composition through the energetics of protein-induced lipid bilayer thickness deformations as well as a chemical coupling driven by preferential interactions between particular lipid and protein species. We find that the resulting lipid-protein organization can endow membrane proteins with diverse and controlled mechanical environments that, via protein-induced lipid bilayer thickness deformations, can strongly influence protein function. The theoretical approach employed here provides a general framework for the quantitative prediction of how membrane thickness deformations influence the joint organization and function of lipids and proteins in cell membranes.

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.

Self-assembly of polyhedral bilayer vesicles from Piezo ion channels

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 a methodology describing the self-assembly of lipids and Piezo proteins 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 proteins. 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.

Vesicle budding induced by binding of curvature-inducing proteins

Vesicle budding induced by protein binding that generates an isotropic spontaneous curvature is studied using a mean-field theory. Many spherical buds are formed via protein binding. As the binding chemical potential increases, the proteins first bind to the buds and then to the remainder of the vesicle. For a high spontaneous curvature and/or high bending rigidity of the bound membrane, it is found that a first-order transition occurs between a small number of large buds and a large number of small buds. These two states coexist around the transition point. The proposed scheme is simple and easily applicable to many interaction types, so we investigate the effects of interprotein interactions, the protein-insertion-induced changes in area, the variation of the saddle-splay modulus, and the area-difference-elasticity energy. The differences in the preferred curvatures for curvature sensing and generation are also clarified.

Unraveling topology-induced shape transformations in dendrimersomes

The vital functions of cell membranes require their ability to quickly change shape to perform complex tasks such as motion, division, endocytosis, and apoptosis. Membrane curvature in cells is modulated by very complex processes such as changes in lipid composition, the oligomerization of curvature-scaffolding proteins, and the reversible insertion of protein regions that act like wedges in the membrane. But, could much simpler mechanisms support membrane shape transformation? In this work, we demonstrate how the change of amphiphile topology in the bilayer can drive shape transformations of cell membrane models. To tackle this, we have designed and synthesized new types of amphiphiles-Janus dendrimers-that self-assemble into uni-, multilamellar, or smectic-ordered vesicles, named dendrimersomes. We synthesized Janus dendrimers containing a photo-labile bond that upon UV-Vis irradiation cleavage lose a part of the hydrophilic dendron. This leads to a change from a cylindrically to a wedge-shaped amphiphile. The high mobility of these dendrimers allows for the concentration of the wedge-shaped amphiphiles and the generation of transmembrane asymmetries. The concentration of the wedges and their rate of segregation allowed control of the budding and generation of structures such as tubules and high genus vesicles.

Regulation of membrane proteins through local heterogeneity in lipid bilayer thickness

Cell membranes show an intricate organization of lipids and membrane proteins into domains with distinct composition and hydrophobic thickness. Using mechanosensitive ion channels as a model system, we employ the membrane elasticity theory of lipid-protein interactions together with the Landau-Ginzburg theory of lipid domain formation to quantify protein-induced lipid bilayer thickness deformations in lipid bilayers with heterogeneous hydrophobic thickness. We show that protein-induced lipid bilayer thickness deformations yield, without any assumptions about preferential interactions between particular lipid and protein species, organization of lipids and membrane proteins according to their preferred hydrophobic thickness, and couple the conformational states of membrane proteins to the local membrane composition. Our calculations suggest that protein-induced lipid bilayer thickness deformations endow proteins in cell membranes with diverse and controlled mechanical environments that, in turn, allow targeted regulation of membrane proteins.

Diffusion on Membrane Domes, Tubes, and Pearling Structures

Diffusion is a fundamental mechanism for protein distribution in cell membranes. These membranes often exhibit complex shapes, which range from shallow domes to elongated tubular or pearl-like structures. Shape complexity of the membrane influences the diffusive spreading of proteins and molecules. Despite the importance membrane geometry plays in these diffusive processes, it is challenging to establish the dependence between diffusion and membrane morphology. We solve the diffusion equation numerically on various static curved shapes representative for experimentally observed membrane shapes. Our results show that membrane necks become diffusion barriers. We determine the diffusive half-time, i.e., the time that is required to reduce the amount of protein in the budded region by one half, and find a quadratic relation between the diffusive half-time and the averaged mean curvature of the membrane shape, which we rationalize by a scaling law. Our findings thus help estimate the characteristic diffusive timescale based on the simple measure of membrane mean curvature.

Effect of protein steric constraints on the symmetry of membrane protein polyhedra

Experiments have shown that, in an aqueous environment, lipids and membrane proteins can self-assemble into membrane protein polyhedral nanoparticles (MPPNs). MPPNs are closed, spherical vesicles composed of a lipid bilayer membrane and membrane proteins, with a polyhedral arrangement of membrane proteins. The observed symmetry and size of MPPNs can be understood from the interplay of protein-induced lipid bilayer deformations in MPPNs, topological defects in protein packing necessitated by the spherical shape of MPPNs, and thermal fluctuations in MPPN self-assembly. We explore here the effect of protein steric constraints on MPPN shape. The protein steric constraints considered here may arise from a well-defined shape of protein domains outside the membrane, entropic repulsion between membrane proteins with flexible domains outside the membrane, or binding of other molecules to membrane proteins. Calculating MPPN self-assembly diagrams under protein steric constraints we find that protein steric constraints can strongly affect MPPN self-assembly. Depending on the specific scenario considered, protein steric constraints can leave large portions of the MPPN self-assembly diagrams with no clearly defined MPPN symmetry or substantially expand the regions of MPPN self-assembly diagrams dominated by highly symmetric MPPN states, such as MPPNs with icosahedral or snub cube symmetry. Our results suggest that modification of protein steric constraints may allow the directed self-assembly of MPPNs with specified symmetry, size, and protein composition and may thus facilitate the further utilization of MPPNs for membrane protein structural analysis or targeted drug delivery.

Gaussian curvature and the budding kinetics of enveloped viruses

The formation of a membrane-enveloped virus starts with the assembly of a curved layer of capsid proteins lining the interior of the plasma membrane (PM) of the host cell. This layer develops into a spherical shell (capsid) enveloped by a lipid-rich membrane. In many cases, the budding process stalls prior to the release of the virus. Recently, Brownian dynamics simulations of a coarse-grained model system reproduced protracted pausing and stalling, which suggests that the origin of pausing/stalling is to be found in the physics of the budding process. Here, we propose that the pausing/stalling observed in the simulations can be understood as a purely kinetic phenomenon associated with the neck geometry. A geometrical potential energy barrier develops during the budding that must be overcome by capsid proteins diffusing along the membrane prior to incorporation into the capsid. The barrier is generated by a conflict between the positive Gauss curvature of the assembling capsid and the negative Gauss curvature of the neck region. A continuum theory description is proposed and is compared with the Brownian simulations of the budding of enveloped viruses.

Despite intense study, the life-cycle of the HIV-1 virus continues to pose mysteries. One of these is the fact that the assembly of an HIV-1 virus along the plasma membrane (PM) of the host cell—the budding process—stalls prior to release of the virus. Many other important viral pathogens with a surrounding lipid membrane envelope display similar stalling. Combining numerical and analytical methods, we demonstrate that the neck-like shape of the membrane that forms prior to release of the virus creates a barrier that blocks the proteins required for the assembly process from reaching the budding virus. An improved understanding of the physics of the blocking process could enable new strategies to combat enveloped viruses.

Supramolecular organization of membrane proteins with anisotropic hydrophobic thickness

Experiments have revealed that membrane proteins often self-assemble into locally ordered clusters. Such membrane protein lattices can play key roles in the functional organization of cell membranes. Membrane protein organization can be driven, at least in part, by bilayer-mediated elastic interactions between membrane proteins. For membrane proteins with anisotropic hydrophobic thickness, bilayer-mediated protein interactions are inherently directional. Here we establish general relations between anisotropy in membrane protein hydrophobic thickness and supramolecular membrane protein organization. We show that protein symmetry is distinctively reflected in the energy landscape of bilayer-mediated protein interactions, favoring characteristic lattice architectures of membrane protein clusters. We find that, in the presence of thermal fluctuations, anisotropy in protein hydrophobic thickness can induce membrane proteins to form mesh-like structures dividing the membrane into compartments. Our results help to elucidate the physical principles and mechanisms underlying the functional organization of cell membranes.

Piezo's membrane footprint and its contribution to mechanosensitivity

Piezo1 is an ion channel that gates open when mechanical force is applied to a cell membrane, thus allowing cells to detect and respond to mechanical stimulation. Molecular structures of Piezo1 reveal a large ion channel with an unusually curved shape. This study analyzes how such a curved ion channel interacts energetically with the cell membrane. Through membrane mechanical calculations, we show that Piezo1 deforms the membrane shape outside the perimeter of the channel into a curved 'membrane footprint'. This membrane footprint amplifies the sensitivity of Piezo1 to changes in membrane tension, rendering it exquisitely responsive. We assert that the shape of the Piezo channel is an elegant example of molecular form evolved to optimize a specific function, in this case tension sensitivity. Furthermore, the predicted influence of the membrane footprint on Piezo gating is consistent with the demonstrated importance of membrane-cytoskeletal attachments to Piezo gating.

Modeling Membrane Curvature Generation due to Membrane⁻Protein Interactions

To alter and adjust the shape of the plasma membrane, cells harness various mechanisms of curvature generation. Many of these curvature generation mechanisms rely on the interactions between peripheral membrane proteins, integral membrane proteins, and lipids in the bilayer membrane. Mathematical and computational modeling of membrane curvature generation has provided great insights into the physics underlying these processes. However, one of the challenges in modeling these processes is identifying the suitable constitutive relationships that describe the membrane free energy including protein distribution and curvature generation capability. Here, we review some of the commonly used continuum elastic membrane models that have been developed for this purpose and discuss their applications. Finally, we address some fundamental challenges that future theoretical methods need to overcome to push the boundaries of current model applications.

Lipid self-assembly and lectin-induced reorganization of the plasma membrane

The plasma membrane represents an outstanding example of self-organization in biology. It plays a vital role in protecting the integrity of the cell interior and regulates meticulously the import and export of diverse substances. Its major building blocks are proteins and lipids, which self-assemble to a fluid lipid bilayer driven mainly by hydrophobic forces. Even if the plasma membrane appears-globally speaking-homogeneous at physiological temperatures, the existence of specialized nano- to micrometre-sized domains of raft-type character within cellular and synthetic membrane systems has been reported. It is hypothesized that these domains are the origin of a plethora of cellular processes, such as signalling or vesicular trafficking. This review intends to highlight the driving forces of lipid self-assembly into a bilayer membrane and the formation of small, transient domains within the plasma membrane. The mechanisms of self-assembly depend on several factors, such as the lipid composition of the membrane and the geometry of lipids. Moreover, the dynamics and organization of glycosphingolipids into nanometre-sized clusters will be discussed, also in the context of multivalent lectins, which cluster several glycosphingolipid receptor molecules and thus create an asymmetric stress between the two membrane leaflets, leading to tubular plasma membrane invaginations.This article is part of the theme issue 'Self-organization in cell biology'.

Nanoparticle wrapping at small non-spherical vesicles: curvatures at play

Nanoparticles in biological systems encounter lipid-bilayer membranes as barriers. They interact with plasma membranes, membranous organelles, such as the endoplasmic reticulum and the Golgi apparatus, the nucleus, and intracellular and extracellular vesicles, such as autophagosomes, lysosomes, and exosomes. Extracellular vesicles have recently attracted particular attention, as they are involved in the transmission of biological signals and as regulators for biological processes. For example, exosomes, small vesicles containing proteins, mRNA, and miRNA, that are released by cells into the extracellular environment, have been suggested to participate in tumor metastasis. Furthermore, vesicles can be applied as targeted-drug-delivery systems. We systematically characterize wrapping of spherical nanoparticles that enter and exit vesicles, depending on particle size, vesicle size, vesicle reduced volume, and membrane spontaneous curvature. We predict the complex wrapping behavior, in particular for large particle-to-vesicle size ratios, where the shape changes of the free membrane contribute significantly to the deformation energy and where nanoparticle wrapping transitions and vesicle shape transitions are coupled. Partial-wrapped membrane-bound particles impose boundary conditions on the membrane that stabilise oblates and stomatocytes for particle entry, and prolates and stomatocytes for particle exit. Our results suggest that nanoparticles may stimulate autophagocytic engulfment, which would facilitate transport of the nanoparticles into lysosomes and would lead to subsequent degradation of nanoparticle-attached proteins.

Clustering on Membranes: Fluctuations and More

Clustering of extracellular ligands and proteins on the plasma membrane is required to perform specific cellular functions, such as signaling and endocytosis. Attractive forces that originate in perturbations of the membrane's physical properties contribute to this clustering, in addition to direct protein-protein interactions. However, these membrane-mediated forces have not all been equally considered, despite their importance. In this review, we describe how line tension, lipid depletion, and membrane curvature contribute to membrane-mediated clustering. Additional attractive forces that arise from protein-induced perturbation of a membrane's fluctuations are also described. This review aims to provide a survey of the current understanding of membrane-mediated clustering and how this supports precise biological functions.

Membrane-Mediated Cooperativity of Proteins

Besides direct protein-protein interactions, indirect interactions mediated by membranes play an important role for the assembly and cooperative function of proteins in membrane shaping and adhesion. The intricate shapes of biological membranes are generated by proteins that locally induce membrane curvature. Indirect curvature-mediated interactions between these proteins arise because the proteins jointly affect the bending energy of the membranes. These curvature-mediated interactions are attractive for crescent-shaped proteins and are a driving force in the assembly of the proteins during membrane tubulation. Membrane adhesion results from the binding of receptor and ligand proteins that are anchored in the apposing membranes. The binding of these proteins strongly depends on nanoscale shape fluctuations of the membranes, leading to a fluctuation-mediated binding cooperativity. A length mismatch between receptor-ligand complexes in membrane adhesion zones causes repulsive curvature-mediated interactions that are a driving force for the length-based segregation of proteins during membrane adhesion.

Curvature-Mediated Assembly of Janus Nanoparticles on Membrane Vesicles

Besides direct particle-particle interactions, nanoparticles adsorbed to biomembranes experience indirect interactions that are mediated by the membrane curvature arising from particle adsorption. In this Letter, we show that the curvature-mediated interactions of adsorbed Janus particles depend on the initial curvature of the membrane prior to adsorption, that is, on whether the membrane initially bulges toward or away from the particles in our simulations. The curvature-mediated interaction can be strongly attractive for Janus particles adsorbed to the outside of a membrane vesicle, which initially bulges away from the particles. For Janus particles adsorbed to the vesicle inside, in contrast, the curvature-mediated interactions are repulsive. We find that the area fraction of the adhesive Janus particle surface is an important control parameter for the curvature-mediated interaction and assembly of the particles, besides the initial membrane curvature.

Acceleration and suppression of banana-shaped-protein-induced tubulation by addition of small membrane inclusions of isotropic spontaneous curvatures

Membrane tubulation induced by banana-shaped protein rods is investigated by using coarse-grained meshless membrane simulations. It is found that tubulation is promoted by laterally isotropic membrane inclusions that generate the same sign of spontaneous curvature as the adsorbed protein rods. The inclusions are concentrated in the tubules and reduce the bending energy of the tip of the tubules. On the other hand, inclusions with an opposite curvature suppress tubulation by percolated-network formation at a high protein-rod density while they induce the formation of a spherical membrane bud at a low rod density. When equal amounts of the two types of inclusions (with positive and negative curvatures) are added, their effects cancel each other for the first short period but later the tubulation is slowly accelerated. Positive surface tension suppresses tubulation. Our results suggest that the cooperation of scaffolding of BAR (Bin/Amphiphysin/Rvs) domains and isotropic membrane inclusions is important for tubulation.

Conformations of a charged vesicle interacting with an oppositely charged particle

Endocytotic and exocytotic processes are usually studied using particle-vesicle systems in theory, but most of them are electroneutral. Nevertheless, charged particle-vesicle systems are much closer to real biological systems. Therefore, wrapping behaviors of a negatively charged vesicle wrapping a positively charged particle are systematically investigated by a series of 2D dynamical simulations in this article. The competition between the elastic bending energy and the electrostatic energy dictates the vesicle configuration and charge distribution. It is found that only for intermediate charge concentrations and small particle sizes a vesicle can completely engulf the particle. When the charge density is high, the interaction between vesicle and particle is unexpectedly weakened by both the hardening effect of the charged membrane and the effective-transportation-frozen effect of the charged components. When the particle is strongly charged, multi-layer folding conformations are observed. These studies may provide important insights into mechanism of endocytotic and exocytotic processes in biological systems.

Conformations and membrane-driven self-organization of rodlike fd virus particles on freestanding lipid membranes

Membrane-mediated interactions and aggregation of colloidal particles adsorbed to responsive elastic membranes are challenging problems relevant for understanding the microscopic organization and dynamics of biological membranes. We experimentally study the behavior of rodlike semiflexible fd virus particles electrostatically adsorbed to freestanding cationic lipid membranes and find that their behavior can be controlled by tuning the membrane charge and ionic strength of the surrounding medium. Three distinct interaction regimes of rodlike virus particles with responsive elastic membranes can be observed. (i) A weakly charged freestanding cationic lipid bilayer in a low ionic strength medium represents a gentle quasi-2D substrate preserving the integrity, structure, and mechanical properties of the membrane-bound semiflexible fd virus, which under these conditions is characterized by a monomer length of 884 ± 4 nm and a persistence length of 2.5 ± 0.2 μm, in perfect agreement with its properties in bulk media. (ii) An increase in the membrane charge leads to the membrane-driven collapse of fd virus particles on freestanding lipid bilayers and lipid nanotubes into compact globules. (iii) When the membrane charge is low, and the mutual electrostatic repulsion of membrane-bound virus particles is screened to a considerable degree, membrane-driven self-organization of membrane-bound fd virus particles into long linear tip-to-tip aggregates showing dynamic self-assembly/disassembly and quasi-semiflexible behavior takes place. These observations are in perfect agreement with the results of recent theoretical and simulation studies predicting that membrane-mediated interactions can control the behavior of colloidal particles adsorbed on responsive elastic membranes.

Formation and Stability of Lipid Membrane Nanotubes

Lipid membrane nanotubes are abundant in living cells, even though tubules are energetically less stable than sheet-like structures. According to membrane elastic theory, the tubular endoplasmic reticulum (ER), with its high area-to-volume ratio, appears to be particularly unstable. We explore how tubular membrane structures can nevertheless be induced and why they persist. In Monte Carlo simulations of a fluid-elastic membrane model subject to thermal fluctuations and without constraints on symmetry, we find that a steady increase in the area-to-volume ratio readily induces tubular structures. In simulations mimicking the ER wrapped around the cell nucleus, tubules emerge naturally as the membrane area increases. Once formed, a high energy barrier separates tubules from the thermodynamically favored sheet-like membrane structures. Remarkably, this barrier persists even at large area-to-volume ratios, protecting tubules against shape transformations despite enormous driving forces toward sheet-like structures. Molecular dynamics simulations of a molecular membrane model confirm the metastability of tubular structures. Volume reduction by osmotic regulation and membrane area growth by lipid production and by fusion of small vesicles emerge as powerful factors in the induction and stabilization of tubular membrane structures.

Nano- and microparticles at fluid and biological interfaces

Systems with interfaces are abundant in both technological applications and biology. While a fluid interface separates two fluids, membranes separate the inside of vesicles from the outside, the interior of biological cells from the environment, and compartmentalize cells into organelles. The physical properties of interfaces are characterized by interface tension, those of membranes are characterized by bending and stretching elasticity. Amphiphilic molecules like surfactants that are added to a system with two immiscible fluids decrease the interface tension and induce a bending rigidity. Lipid bilayer membranes of vesicles can be stretched or compressed by osmotic pressure; in biological cells, also the presence of a cytoskeleton can induce membrane tension. If the thickness of the interface or the membrane is small compared with its lateral extension, both can be described using two-dimensional mathematical surfaces embedded in three-dimensional space. We review recent work on the interaction of particles with interfaces and membranes. This can be micrometer-sized particles at interfaces that stabilise emulsions or form colloidosomes, as well as typically nanometer-sized particles at membranes, such as viruses, parasites, and engineered drug delivery systems. In both cases, we first discuss the interaction of single particles with interfaces and membranes, e.g. particles in external fields, non-spherical particles, and particles at curved interfaces, followed by interface-mediated interaction between two particles, many-particle interactions, interface and membrane curvature-induced phenomena, and applications.

Membrane structure formation induced by two types of banana-shaped proteins

The assembly of banana-shaped rodlike proteins on membranes and the associated membrane shape transformations are investigated by analytical theory and coarse-grained simulations. The membrane-mediated interactions between two banana-shaped inclusions are derived theoretically using a point-like formalism based on fixed anisotropic curvatures, both for zero surface tension and for finite surface tension. On a larger scale, the interactions between the assemblies of such rodlike inclusions are determined analytically. Meshless membrane simulations are performed in the presence of a large number of inclusions of two types, corresponding to the curved rods of opposite curvatures, both for flat membranes and vesicles. Rods of the same type aggregate into linear assemblies perpendicular to the rod axis, leading to membrane tubulation. However, rods of the other type, those of opposite curvature, are attracted to the lateral sides of these assemblies, and stabilize a straight bump structure that prevents tubulation. When the two types of rods have almost opposite curvatures, the bumps attract one another, forming a striped structure. Positive surface tension is found to stabilize stripe formation. The simulation results agree well with the theoretical predictions provided the point-like curvatures of the model are scaled-down to account for the effective flexibility of the simulated rods.

Curvature-Induced Spatial Ordering of Composition in Lipid Membranes

Phase segregation of membranal components, such as proteins, lipids, and cholesterols, leads to the formation of aggregates or domains that are rich in specific constituents. This process is important in the interaction of the cell with its surroundings and in determining the cell's behavior and fate. Motivated by published experiments on curvature-modulated phase separation in lipid membranes, we formulate a mathematical model aiming at studying the spatial ordering of composition in a two-component biomembrane that is subjected to a prescribed (imposed) geometry. Based on this model, we identified key nondimensional quantities that govern the biomembrane response and performed numerical simulations to quantitatively explore their influence. We reproduce published experimental observations and extend them to surfaces with geometric features (imposed geometry) and lipid phases beyond those used in the experiments. In addition, we demonstrate the possibility for curvature-modulated phase separation above the critical temperature and propose a systematic procedure to determine which mechanism, the difference in bending stiffness or difference in spontaneous curvatures of the two phases, dominates the coupling between shape and composition.

Pointlike Inclusion Interactions in Tubular Membranes

Membrane tubes and tubular networks are ubiquitous in living cells. Inclusions like proteins are vital for both the stability and the dynamics of such networks. These inclusions interact via the curvature deformations they impose on the membrane. We analytically study the resulting membrane mediated interactions in strongly curved tubular membranes. We model inclusions as constraints coupled to the curvature tensor of the membrane tube. First, as special test cases, we analyze the interaction between ring- and rod-shaped inclusions. Using Monte Carlo simulations, we further show how pointlike inclusions interact to form linear aggregates. To minimize the curvature energy of the membrane, inclusions self-assemble into either line- or ringlike patterns. Our results show that the global curvature of the membrane strongly affects the interactions between proteins embedded in it, and can lead to the spontaneous formation of biologically relevant structures.

Symmetry and Size of Membrane Protein Polyhedral Nanoparticles

In recent experiments [T. Basta et al., Proc. Natl. Acad. Sci. U.S.A. 111, 670 (2014)] lipids and membrane proteins were observed to self-assemble into membrane protein polyhedral nanoparticles (MPPNs) with a well-defined polyhedral protein arrangement and characteristic size. We develop a model of MPPN self-assembly in which the preferred symmetry and size of MPPNs emerge from the interplay of protein-induced lipid bilayer deformations, topological defects in protein packing, and thermal effects. With all model parameters determined directly from experiments, our model correctly predicts the observed symmetry and size of MPPNs. Our model suggests how key lipid and protein properties can be modified to produce a range of MPPN symmetries and sizes in experiments.

Shape deformation of lipid membranes by banana-shaped protein rods: Comparison with isotropic inclusions and membrane rupture

The assembly of curved protein rods on fluid membranes is studied using implicit-solvent meshless membrane simulations. As the rod curvature increases, the rods on a membrane tube assemble along the azimuthal direction first and subsequently along the longitudinal direction. Here, we show that both transition curvatures decrease with increasing rod stiffness. For comparison, curvature-inducing isotropic inclusions are also simulated. When the isotropic inclusions have the same bending rigidity as the other membrane regions, the inclusions are uniformly distributed on the membrane tubes and vesicles even for large spontaneous curvature of the inclusions. However, the isotropic inclusions with much larger bending rigidity induce shape deformation and are concentrated on the region of a preferred curvature. For high rod density, high rod stiffness, and/or low line tension of the membrane edge, the rod assembly induces vesicle rupture, resulting in the formation of a high-genus vesicle. A gradual change in the curvature suppresses this rupture. Hence, large stress, compared to the edge tension, induced by the rod assembly is the key factor determining rupture. For rod curvature with the opposite sign to the vesicle curvature, membrane rupture induces inversion of the membrane, leading to division into multiple vesicles as well as formation of a high-genus vesicle.

Bilayer-thickness-mediated interactions between integral membrane proteins

Hydrophobic thickness mismatch between integral membrane proteins and the surrounding lipid bilayer can produce lipid bilayer thickness deformations. Experiment and theory have shown that protein-induced lipid bilayer thickness deformations can yield energetically favorable bilayer-mediated interactions between integral membrane proteins, and large-scale organization of integral membrane proteins into protein clusters in cell membranes. Within the continuum elasticity theory of membranes, the energy cost of protein-induced bilayer thickness deformations can be captured by considering compression and expansion of the bilayer hydrophobic core, membrane tension, and bilayer bending, resulting in biharmonic equilibrium equations describing the shape of lipid bilayers for a given set of bilayer-protein boundary conditions. Here we develop a combined analytic and numerical methodology for the solution of the equilibrium elastic equations associated with protein-induced lipid bilayer deformations. Our methodology allows accurate prediction of thickness-mediated protein interactions for arbitrary protein symmetries at arbitrary protein separations and relative orientations. We provide exact analytic solutions for cylindrical integral membrane proteins with constant and varying hydrophobic thickness, and develop perturbative analytic solutions for noncylindrical protein shapes. We complement these analytic solutions, and assess their accuracy, by developing both finite element and finite difference numerical solution schemes. We provide error estimates of our numerical solution schemes and systematically assess their convergence properties. Taken together, the work presented here puts into place an analytic and numerical framework which allows calculation of bilayer-mediated elastic interactions between integral membrane proteins for the complicated protein shapes suggested by structural biology and at the small protein separations most relevant for the crowded membrane environments provided by living cells.

Membrane tubule formation by banana-shaped proteins with or without transient network structure

In living cells, membrane morphology is regulated by various proteins. Many membrane reshaping proteins contain a Bin/Amphiphysin/Rvs (BAR) domain, which consists of a banana-shaped rod. The BAR domain bends the biomembrane along the rod axis and the features of this anisotropic bending have recently been studied. Here, we report on the role of the BAR protein rods in inducing membrane tubulation, using large-scale coarse-grained simulations. We reveal that a small spontaneous side curvature perpendicular to the rod can drastically alter the tubulation dynamics at high protein density, whereas no significant difference is obtained at low density. A percolated network is intermediately formed depending on the side curvature. This network suppresses tubule protrusion, leading to the slow formation of fewer tubules. Thus, the side curvature, which is generated by protein–protein and membrane–protein interactions, plays a significant role in tubulation dynamics. We also find that positive surface tensions and the vesicle membrane curvature can stabilize this network structure by suppressing the tubulation.

Non-Brownian diffusion in lipid membranes: Experiments and simulations

The dynamics of constituents and the surface response of cellular membranes-also in connection to the binding of various particles and macromolecules to the membrane-are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder-through the addition of cholesterol or proteins-and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments-the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane-binding particles. We discuss how membrane compartmentalisation and the particle-membrane binding energy may impact the dynamics and response of lipid membranes. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Architecture and Function of Mechanosensitive Membrane Protein Lattices

Experiments have revealed that membrane proteins can form two-dimensional clusters with regular translational and orientational protein arrangements, which may allow cells to modulate protein function. However, the physical mechanisms yielding supramolecular organization and collective function of membrane proteins remain largely unknown. Here we show that bilayer-mediated elastic interactions between membrane proteins can yield regular and distinctive lattice architectures of protein clusters, and may provide a link between lattice architecture and lattice function. Using the mechanosensitive channel of large conductance (MscL) as a model system, we obtain relations between the shape of MscL and the supramolecular architecture of MscL lattices. We predict that the tetrameric and pentameric MscL symmetries observed in previous structural studies yield distinct lattice architectures of MscL clusters and that, in turn, these distinct MscL lattice architectures yield distinct lattice activation barriers. Our results suggest general physical mechanisms linking protein symmetry, the lattice architecture of membrane protein clusters, and the collective function of membrane protein lattices.

Modulation of apical constriction by Wnt signaling is required for lung epithelial shape transition

In lung development the apically constricted columnar epithelium forms numerous buds during the pseudoglandular stage and subsequently changes the shape into flat or cuboidal pneumocytes that compose the air sacs during the canalicular and saccular (canalicular-saccular) stages, yet the impact of cell shapes on tissue morphogenesis remains unclear. The expression of Wnt components were decreased in the canalicular-saccular stages, and genetically constitutive activation of Wnt signaling impaired air sac formation by inducing apical constriction in the epithelium as seen in the pseudoglandular stage. Organ culture models also demonstrated that Wnt signaling induces apical constriction through the apical actomyosin cytoskeletal organization. Mathematical modeling revealed that apical constriction induces bud formation and loss of apical constriction is required for the formation of an air sac-like structure. MAP/Microtubule affinity-regulating kinase (MARK1) was identified as a downstream molecule of Wnt signaling and required for the apical cytoskeletal organization and bud formation. These results suggest that Wnt signaling is required for bud formation by inducing apical constriction during the pseudoglandular stage, while loss of Wnt signaling is for air sac formation in the canalicular-saccular stages.

Application of a free-energy-landscape approach to study tension-dependent bilayer tubulation mediated by curvature-inducing proteins

We investigate the phenomenon of protein-induced tubulation of lipid bilayer membranes within a continuum framework using Monte Carlo simulations coupled with the Widom insertion technique to compute excess chemical potentials. Tubular morphologies are spontaneously formed when the density and the curvature-field strength of the membrane-bound proteins exceed their respective thresholds and this transition is marked by a sharp drop in the excess chemical potential. We find that the planar to tubular transition can be described by a micellar model and that the corresponding free-energy barrier increases with an increase in the curvature-field strength (i.e., of protein-membrane interactions) and also with an increase in membrane tension.

Membrane-wrapping contributions to malaria parasite invasion of the human erythrocyte

The blood stage malaria parasite, the merozoite, has a small window of opportunity during which it must successfully target and invade a human erythrocyte. The process of invasion is nonetheless remarkably rapid. To date, mechanistic models of invasion have focused predominantly on the parasite actomyosin motor contribution to the energetics of entry. Here, we have conducted a numerical analysis using dimensions for an archetypal merozoite to predict the respective contributions of the host-parasite interactions to invasion, in particular the role of membrane wrapping. Our theoretical modeling demonstrates that erythrocyte membrane wrapping alone, as a function of merozoite adhesive and shape properties, is sufficient to entirely account for the first key step of the invasion process, that of merozoite reorientation to its apex and tight adhesive linkage between the two cells. Next, parasite-induced reorganization of the erythrocyte cytoskeleton and release of parasite-derived membrane can also account for a considerable energetic portion of actual invasion itself, through membrane wrapping. Thus, contrary to the prevailing dogma, wrapping by the erythrocyte combined with parasite-derived membrane release can markedly reduce the expected contributions of the merozoite actomyosin motor to invasion. We therefore propose that invasion is a balance between parasite and host cell contributions, evolved toward maximal efficient use of biophysical forces between the two cells.

How hydrophobic nanoparticles aggregate in the interior of membranes: A computer simulation

Lipid-based dispersion of hydrophobic nanoparticles (NPs) not only gives fundamental insight into how nanomaterials distribute in live cells and organisms, but also provides a quite general route to designing nanocarrier agents in triggered drug delivery and medical imaging. It is not clearly understood how hydrophobic NPs arrange in the interior of a membrane. In this paper, with computer simulation techniques, we demonstrate that hydrophobic NPs having a diameter compared to the hydrophobic thickness of the membrane are capable of clustering in the hydrophobic interior of a cell membrane. Except from the isotropic aggregation, an unexpected linear arrangement of spherical NPs, which is still not found from experiments, is identified here. The free-energy costs associated with linear and isotropic aggregations are computed explicitly to interpret aggregation behavior and the obtained phase diagrams give us a comprehensive understanding of where linear aggregation is expected. In this work we also shows that NP size and membrane tension play key roles in determining the NP aggregate, while the effects of NP concentration and membrane curvature seem to be relatively weak.

Capillary assembly of microscale ellipsoidal, cuboidal, and spherical particles at interfaces

Micron-sized anisotropic particles with homogeneous surface properties at a fluid interface can deform the interface due to their shape. The particles thereby create excess interfacial area and interact in order to minimize this area, which lowers the total interfacial energy. We present a systematic investigation of the interface deformations around single ellipsoidal particles and cuboidal particles with rounded edges in the near field for various contact angles and particle aspect ratios. The correlation of these deformations with capillary bond energies-the interaction energies of two particles at contact-quantifies the relation between the interactions and the near-field deformations. We characterize the interactions using effective power laws and investigate how anisotropic particles self-assemble by capillary forces. Interface deformations and particle interactions for cuboidal particles are weaker compared with those for ellipsoidal particles with the same aspect ratios. For both particle shapes, the bound state in side-by-side orientation is most stable, while the interaction in tip-to-side orientation is repulsive. Furthermore, we find capillary attraction between spherical and ellipsoidal particles. Our calculations therefore suggest cluster formation of spherical and ellipsoidal particles, which elucidates the role of spherical particles as stoppers for the growth of worm-like chains of ellipsoidal particles. The interaction between spherical and ellipsoidal particles might also explain the suppression of the "coffee-ring effect" that has been observed for evaporating droplets with mixtures of spherical and ellipsoidal particles. In general, our calculations of the near-field interactions complement previous calculations in the far field and help to predict colloidal assembly and rheological properties of particle-laden interfaces.

Shape and orientation matter for the cellular uptake of nonspherical particles

Recent advances in nanotechnology have made a whole zoo of particles of different shapes available for applications, but their interaction with biological cells and their toxicity is often not well understood. Experiments have shown that particle uptake by cells is determined by an intricate interplay between physicochemical particle properties like shape, size, and surface functionalization, but also by membrane properties and particle orientation. Our work provides systematic understanding, based on a mechanical description, for membrane wrapping of nanoparticles, viruses, and bacterial forms. For rod-like particles, we find stable endocytotic states with small and high wrapping fraction; an increased aspect ratio is unfavorable for complete wrapping. For high aspect ratios and round tips, the particles enter via a submarine mode, side-first with their long edge parallel to the membrane. For small aspect ratios and flat tips, the particles enter tip-first via a rocket mode.

Membrane-elasticity model of Coatless vesicle budding induced by ESCRT complexes

The formation of vesicles is essential for many biological processes, in particular for the trafficking of membrane proteins within cells. The Endosomal Sorting Complex Required for Transport (ESCRT) directs membrane budding away from the cytosol. Unlike other vesicle formation pathways, the ESCRT-mediated budding occurs without a protein coat. Here, we propose a minimal model of ESCRT-induced vesicle budding. Our model is based on recent experimental observations from direct fluorescence microscopy imaging that show ESCRT proteins colocalized only in the neck region of membrane buds. The model, cast in the framework of membrane elasticity theory, reproduces the experimentally observed vesicle morphologies with physically meaningful parameters. In this parameter range, the minimum energy configurations of the membrane are coatless buds with ESCRTs localized in the bud neck, consistent with experiment. The minimum energy configurations agree with those seen in the fluorescence images, with respect to both bud shapes and ESCRT protein localization. On the basis of our model, we identify distinct mechanistic pathways for the ESCRT-mediated budding process. The bud size is determined by membrane material parameters, explaining the narrow yet different bud size distributions in vitro and in vivo. Our membrane elasticity model thus sheds light on the energetics and possible mechanisms of ESCRT-induced membrane budding.

Curvature-driven lipid sorting in biomembranes

It has often been suggested that the high curvature of transport intermediates in cells may be a sufficient means to segregate different lipid populations based on the relative energy costs of forming bent membranes. In this review, we present in vitro experiments that highlight the essential physics of lipid sorting at thermal equilibrium: It is driven by a trade-off between bending energy, mixing entropy, and interactions between species. We collect evidence that lipid sorting depends strongly on lipid-lipid and protein-lipid interactions, and hence on the underlying composition of the membrane and on the presence of bound proteins.

Curling and local shape changes of red blood cell membranes driven by cytoskeletal reorganization

Human red blood cells (RBCs) lack the actin-myosin-microtubule cytoskeleton that is responsible for shape changes in other cells. Nevertheless, they can display highly dynamic local deformations in response to external perturbations, such as those that occur during the process of apical alignment preceding merozoite invasion in malaria. Moreover, after lysis in divalent cation-free media, the isolated membranes of ruptured ghosts show spontaneous inside-out curling motions at the free edges of the lytic hole, leading to inside-out vesiculation. The molecular mechanisms that drive these rapid shape changes are unknown. Here, we propose a molecular model in which the spectrin filaments of the RBC cortical cytoskeleton control the sign and dynamics of membrane curvature depending on whether the ends of the filaments are free or anchored to the bilayer. Computer simulations of the model reveal that curling, as experimentally observed, can be obtained either by an overall excess of weakly-bound filaments throughout the cell, or by the flux of such filaments toward the curling edges. Divalent cations have been shown to arrest the curling process, and Ca2+ ions have also been implicated in local membrane deformations during merozoite invasion. These effects can be replicated in our model by attributing the divalent cation effects to increased filament-membrane binding. This process converts the curl-inducing loose filaments into fully bound filaments that arrest curling. The same basic mechanism can be shown to account for Ca2+-induced local and dynamic membrane deformations in intact RBCs. The implications of these results in terms of RBC membrane dynamics under physiological, pathological, and experimental conditions is discussed.

Phase diagrams for multi-component membrane vesicles: a coarse-grained modeling study

We develop a multicomponent membrane coarse-grained model in which the effects of spontaneous curvature and fluidity are included. This model is used to perform computer simulations to study the phase segregation, domain coarsening, budding and budding off of multicomponent membrane vesicles. Three types of phase diagram are presented with variation of composition, spontaneous curvature, and line tension (the unlike bond energy). Various phases have been observed, including, sphere, biconcave, starfish, capsule, budding, and budding off. Our simulations show that budding off occurs only when the size of a domain is larger than a critical value, and the critical condition is closely related to the spontaneous curvature and line tension. A continuum model is used to predict the critical condition for budding off. Quantitative comparisons are made between the present simulation results and the continuum model predictions, and good agreements have been achieved.