Elastic Membrane Deformations Govern Interleaflet Coupling of Lipid-Ordered Domains

Antivesiculation and Complete Unbinding of Tail-Tethered Lipids

We report the effect of tail-tethering on vesiculation and complete unbinding of bilayered membranes. Amphiphilic molecules of a bolalipid, resembling the tail-tethered molecular structure of archaeal lipids, with two identical zwitterionic phosphatidylcholine headgroups self-assemble into a large flat lamellar membrane, in contrast to the multilamellar vesicles (MLVs) observed in its counterpart, monopolar nontethered zwitterionic lipids. The antivesiculation is confirmed by small-angle X-ray scattering (SAXS) and cryogenic transmission electron microscopy (cyro-TEM). With the net charge of zero and higher bending rigidity of the membrane (confirmed by neutron spin echo (NSE) spectroscopy), the current membrane theory would predict that membranes should stack with each other (aka "bind") due to dominant van der Waals attraction, while the outcome of the nonstacking ("unbinding") membrane suggests that the theory needs to include entropic contribution for the nonvesicular structures. This report pioneers an understanding of how the tail-tethering of amphiphiles affects the structure, enabling better control over the final nanoscale morphology.

The Membrane-Mediated Interaction of Liquid-Ordered Lipid Domains in the Presence of Amphipathic Peptides

The lipid membranes of living cells are composed of a large number of lipid types and can undergo phase separation with the formation of nanometer-scale liquid-ordered lipid domains, also called rafts. Raft coalescence, i.e., the fusion of lipid domains, is involved in important cell processes, such as signaling and trafficking. In this work, within the framework of the theory of elasticity of lipid membranes, we explore how amphipathic peptides adsorbed on lipid membranes may affect the domain-domain fusion processes. We show that the elastic deformations of lipid membranes drive amphipathic peptides to the boundary of lipid domains, which leads to an increase in the average energy barrier of the domain-domain fusion, even if the surface concentration of amphipathic peptides is low and the domain boundaries are only partially occupied by the peptides. This inhibition of the fusion of lipid domains may lead to negative side effects of using amphipathic peptides as antimicrobial agents.

Phosphatidylinositol 4,5-Bisphosphate Sensing Lipid Raft via Inter-Leaflet Coupling Regulated by Acyl Chain Length of Sphingomyelin

Phosphatidylinositol 4,5-bisphosphate (PIP₂) is an important molecule located at the inner leaflet of cell membrane, where it serves as anchoring sites for a cohort of membrane-associated molecules and as a broad-reaching signaling intermediate. The lipid raft is thought as the major platform recruiting proteins for signal transduction and also known to mediate PIP₂ accumulation across the membrane. While the significance of this cross-membrane coupling is increasingly appreciated, it remains unclear whether and how PIP₂ senses the dynamic change of the ordered lipid domains over the packed hydrophobic core of the bilayer. Herein, by means of molecular dynamic simulation, we reveal that inner PIP₂ molecules can sense the outer lipid domain via inter-leaflet coupling, and the coupling manner is dictated by the acyl chain length of sphingomyelin (SM) partitioned to the lipid domain. Shorter SM promotes membrane domain registration, whereby PIP₂ accumulates beneath the domain across the membrane. In contrast, the anti-registration is thermodynamically preferred if the lipid domain has longer SM due to the hydrophobic mismatch between the corresponding acyl chains in SM and PIP₂. In this case, PIP₂ is expelled by the domain with a higher diffusivity. These results provide molecular insights into the regulatory mechanism of correlation between the outer lipid domain and inner PIP₂, both of which are critical components for cell signal transduction.

Determination of elastic parameters of lipid membranes from simulation under varied external pressure

Many cellular processes such as endocytosis, exocytosis, and vesicle trafficking involve membrane deformations, which can be analyzed in the framework of the elastic theories of lipid membranes. These models operate with phenomenological elastic parameters. A connection between these parameters and the internal structure of lipid membranes can be provided by three-dimensional (3D) elastic theories. Considering a membrane as a 3D layer, Campelo et al. [F. Campelo et al., Adv. Colloid Interface Sci. 208, 25 (2014)10.1016/j.cis.2014.01.018] developed a theoretical basis for the calculation of elastic parameters. In this work we generalize and improve this approach by considering a more general condition of global incompressibility instead of local incompressibility. Crucially, we find an important correction to the theory of Campelo et al., which if not taken into account leads to a significant miscalculation of elastic parameters. With the total volume conservation taken into account, we derive an expression for the local Poisson's ratio, which determines how the local volume changes upon stretching and permits a more precise determination of elastic parameters. Also, we substantially simplify the procedure by calculating the derivatives of the moments of the local tension with respect to stretching instead of calculating the local stretching modulus. We obtain a relation between the Gaussian curvature modulus as a function of stretching and the bending modulus, showing that these two elastic parameters are not independent, as was previously assumed. The proposed algorithm is applied to membranes composed of pure dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), and their mixture. The following elastic parameters of these systems are obtained: the monolayer bending and stretching moduli, spontaneous curvature, neutral surface position, and local Poisson's ratio. It is shown that the bending modulus of the DPPC/DOPC mixture follows a more complex trend than predicted by the classical Reuss averaging, which is often employed in theoretical frameworks.

A Mechanism of Double-Membrane Vesicle Formation from Liquid-Ordered/Liquid-Disordered Phase Separated Spherical Membrane

Genome replication of coronaviruses takes place in specific cellular compartments, in so-called double-membrane vesicles (DMVs), formed from the endoplasmic reticulum (ER). An intensive production of DMVs is induced by non-structural viral proteins. Here, we proposed a possible mechanism of the DMV formation from ER-derived spherical vesicles where liquid-ordered and liquid-disordered lipid phases coexist. These vesicles are supposed to divide into two homogeneous liquid-ordered and liquid-disordered vesicles. The formation of two spherical vesicles constituting DMV requires a mechanical work to be performed. We considered the excess energy of the boundary between the coexisting lipid phases as the main driving force behind the division of the initial vesicle. Explicitly accounting for the energy of elastic deformations and the interphase boundary energy, we analyzed a range of physical parameters where the DMV formation is possible. We concluded that this process can principally take place in a very narrow range of system parameters. The most probable diameter of DMVs formed according to the proposed mechanism appeared to be approximately 220 nm, in an agreement with the average diameter of DMVs observed in vivo. Our consideration predicts the DMV size to be strongly limited from above. The developed analysis can be utilized for the production of DMVs in model systems.

Determination of Elastic Parameters of Lipid Membranes with Molecular Dynamics: A Review of Approaches and Theoretical Aspects

Lipid membranes are abundant in living organisms, where they constitute a surrounding shell for cells and their organelles. There are many circumstances in which the deformations of lipid membranes are involved in living cells: fusion and fission, membrane-mediated interaction between membrane inclusions, lipid-protein interaction, formation of pores, etc. In all of these cases, elastic parameters of lipid membranes are important for the description of membrane deformations, as these parameters determine energy barriers and characteristic times of membrane-involved phenomena. Since the development of molecular dynamics (MD), a variety of in silico methods have been proposed for the determination of elastic parameters of simulated lipid membranes. These MD methods allow for the consideration of details unattainable in experimental techniques and represent a distinct scientific field, which is rapidly developing. This work provides a review of these MD approaches with a focus on theoretical aspects. Two main challenges are identified: (i) the ambiguity in the transition from the continuum description of elastic theories to the discrete representation of MD simulations, and (ii) the determination of intrinsic elastic parameters of lipid mixtures, which is complicated due to the composition-curvature coupling effect.

A Model of Lipid Monolayer-Bilayer Fusion of Lipid Droplets and Peroxisomes

Lipid droplets are unique organelles that store neutral lipids encapsulated by the lipid monolayer. In some processes of cellular metabolism, lipid droplets interact with peroxisomes resulting in the fusion of their envelopes and the formation of protrusions of the peroxisome monolayer, called pexopodia. The formation of pexopodia is facilitated by free fatty acids generated during lipolysis within lipid droplets. In this work, we studied the fusion of monolayer and bilayer membranes during the interaction between lipid droplets and peroxisomes. To this end, we built the energy trajectory of this process using the continuum elasticity theory and investigated the molecular details of the fusion structures utilizing molecular dynamics. We divided the fusion process into two stages: formation of a stalk and its consequent expansion into pexopodia. We found that in the considered system, the stalk was energetically more stable and had a lower energy barrier of formation compared to the case of bilayer fusion. The further evolution of the stalk depended on the value of the spontaneous curvature of the membrane in a threshold manner. We attributed the possible expansion of the stalk to the incorporation of free fatty acids into the stalk region. The developed model allowed describing quantitatively the process of monolayer-bilayer fusion.

Lysolipids regulate raft size distribution

The lipid matrix of cellular membranes, directly and indirectly, regulates many vital functions of the cell. The diversity of lipids in membranes leads to the formation of ordered domains called rafts, which play a crucial role in signal transduction, protein sorting and other cellular processes. Rafts are believed to impact the development of different neurodegenerative diseases, such as Alzheimer’s, Parkinson’s, Huntington’s ones, amyotrophic lateral sclerosis, some types of cancer, etc. These diseases correlate with the change in the membrane lipid composition resulting from an oxidative stress, age-related processes, dysfunction of proteins, and many others. In particular, a lot of studies report a significant rise in the level of lysolipids. Physicochemical properties of rafts are determined by membrane composition, in particular, by the content of lysolipids. Lysolipids may thus regulate raft-involving processes. However, the exact mechanism of such regulation is unknown. Although studying rafts in vivo still seems to be rather complicated, liquid-ordered domains are well observed in model systems. In the present study, we used atomic force microscopy (AFM) to examine how lysophospholipids influence the liquid-ordered domains in model ternary membranes. We demonstrated that even a small amount of lysolipids in a membrane significantly impacts domain size depending on the saturation of the lysolipid hydrocarbon tails and the amount of cholesterol. The mixture with the bigger relative fraction of cholesterol was more susceptible to the action of lysolipids. This data helped us to generalize our previous theoretical model of the domain size regulation by lipids with particular molecular shape expanding it to the case of lysolipids and dioleoylglycerol.

Fluid-gel coexistence in lipid membranes under differential stress

A widely conserved property of many biological lipid bilayers is their asymmetry. In addition to having distinct compositions on its two sides, a membrane can also exhibit different tensions in its two leaflets, a state known as differential stress. Here, we examine how this stress can influence the phase behavior of the constituent lipid monolayers of a single-component membrane. For temperatures sufficiently close to, but still above, the main transition, molecular dynamics simulations show the emergence of finite gel domains within the compressed leaflet. We describe the thermodynamics of this phenomenon by adding two empirical single-leaflet free energies for the fluid-gel transition, each evaluated at its respective asymmetry-dependent lipid density. Finite size effects arising in simulation are included in the theory through a geometry-dependent interfacial term. Our model reproduces the phase coexistence observed in simulation. It could therefore be used to connect the "hidden variable" of differential stress to experimentally observable properties of the main phase transition. These ideas could be generalized to any first-order bilayer phase transition in the presence of asymmetry, including liquid-ordered/liquid-disordered phase separation.

Enthalpic and entropic contributions to interleaflet coupling drive domain registration and antiregistration in biological membrane

Biological membrane is a complex self-assembly of lipids, sterols, and proteins organized as a fluid bilayer of two closely stacked lipid leaflets. Differential molecular interactions among its diverse constituents give rise to heterogeneities in the membrane lateral organization. Under certain conditions, heterogeneities in the two leaflets can be spatially synchronized and exist as registered domains across the bilayer. Several contrasting theories behind mechanisms that induce registration of nanoscale domains have been suggested. Following a recent study showing the effect of position of lipid tail unsaturation on domain registration behavior, we decided to develop an analytical theory to elucidate the driving forces that create and maintain domain registry across leaflets. Towards this, we formulated a Hamiltonian for a stacked lattice system where site variables capture the lipid molecular properties such as the position of unsaturation and various other interactions that could drive phase separation and interleaflet coupling. We solve the Hamiltonian using Monte Carlo simulations and create a complete phase diagram that reports the presence or absence of registered domains as a function of various Hamiltonian parameters. We find that the interleaflet coupling should be described as a competing enthalpic contribution due to interaction of lipid tail termini, primarily due to saturated-saturated interactions, and an interleaflet entropic contribution from overlap of unsaturated tail termini. A higher position of unsaturation is seen to provide weaker interleaflet coupling. Thermodynamically stable nanodomains could also be observed for certain points in the parameter space in our bilayer model, which were further verified by carrying out extended Monte Carlo simulations. These persistent noncoalescing registered nanodomains close to the lower end of the accepted nanodomain size range also point towards a possible "nanoscale" emulsion description of lateral heterogeneities in biological membrane leaflets.

Determinants of Lipid Domain Size

Lipid domains less than 200 nm in size may form a scaffold, enabling the concerted function of plasma membrane proteins. The size-regulating mechanism is under debate. We tested the hypotheses that large values of spontaneous monolayer curvature are incompatible with micrometer-sized domains. Here, we used the transition of photoswitchable lipids from their cylindrical conformation to a conical conformation to increase the negative curvature of a bilayer-forming lipid mixture. In contrast to the hypothesis, pre-existing micrometer-sized domains did not dissipate in our planar bilayers, as indicated by fluorescence images and domain mobility measurements. Elasticity theory supports the observation by predicting the zero free energy gain for splitting large domains into smaller ones. It also indicates an alternative size-determining mechanism: The cone-shaped photolipids reduce the line tension associated with lipid deformations at the phase boundary and thus slow down the kinetics of domain fusion. The competing influence of two approaching domains on the deformation of the intervening lipids is responsible for the kinetic fusion trap. Our experiments indicate that the resulting local energy barrier may restrict the domain size in a dynamic system.

Biophysical Reviews' "Meet the Councilor Series"-a profile of Peter Pohl

It is my pleasure to write a few words to introduce myself to the readers of Biophysical Reviews as part of the "Meet the Councilor Series." Currently, I am serving the second period as IUPAB councilor after having been elected first in 2017. Initially, I studied Biophysics in Moscow (Russia) and later Medicine in Halle (Germany). My scientific carrier took me from the Medical School of the Martin Luther University of Halle-Wittenberg, via the Leibniz Institute for Molecular Pharmacology (Berlin) and the Institute for Biology at the Humboldt University (Berlin) to the Physics Department of the Johannes Kepler University in Linz (Austria). My key research interests lie in the molecular mechanisms of transport phenomena occurring at the lipid membrane, including (i) spontaneous and facilitated transport of water and other small molecules across membranes in reconstituted systems, (ii) proton migration along the membrane surface, (iii) protein translocation, and (iv) bilayer mechanics. Training of undergraduate, graduate, and postdoctoral researchers from diverse academic disciplines has been-and shall remain-a consistent part of my work.

Characteristic lengths of transmembrane peptides controlling their tilt and lateral distribution between membrane domains

Lipids and proteins of plasma membranes of eukaryotic cells are supposed to form protein-lipid domains, characterized by a different molecular order, bilayer thickness, and elastic parameters. Several mechanisms of preferable distribution of transmembrane proteins to the ordered or disordered membrane domains have been revealed. The mismatch between the length of the protein transmembrane domain and hydrophobic thickness of the lipid bilayer is considered to be an important driving force of protein lateral sorting. Utilizing the continuum theory of elasticity, we analyzed optimal configurations and preferable membrane domains for single-pass transmembrane peptides of various hydrophobic lengths and effective molecular shapes. We obtained that short transmembrane peptides stand perpendicularly to the membrane plane. The exceedance of a certain characteristic length leads to the tilt of the peptide. This length depends on the bilayer thickness. Thus, in the membrane with coexisting ordered (thicker) and disordered (thinner) phases tilting of the peptide in each phase is governed by its individual characteristic length. The lateral distribution of the peptides between ordered and disordered membrane domains is shown to be described by two additional characteristic lengths. The exceedance of the smaller one drives the peptide towards a more ordered and thicker membrane, while the exceedance of the larger characteristic length switches the preferable membrane domain from ordered and thicker to disordered and thinner. Thus, membrane proteins with long enough transmembrane domains are predicted to accumulate in the thinner disordered membrane as compared to the thicker ordered bilayer. For hourglass-like and barrel-like shaped transmembrane peptides the specific regime of sorting was obtained: the peptides distributed almost equally between the phases in a wide range of peptide lengths. This finding allowed explaining the experimental data on lateral distribution of transmembrane peptide tLAT.

Amphipathic Peptides Impede Lipid Domain Fusion in Phase-Separated Membranes

Cell membranes are heterogeneous in lipid composition which leads to the phase separation with the formation of nanoscopic liquid-ordered domains, also called rafts. There are multiple cell processes whereby the clustering of these domains into a larger one might be involved, which is responsible for such important processes as signal transduction, polarized sorting, or immune response. Currently, antimicrobial amphipathic peptides are considered promising antimicrobial, antiviral, and anticancer therapeutic agents. Here, within the framework of the classical theory of elasticity adapted for lipid membranes, we investigate how the presence of the peptides in a phase-separated membrane influences the fusion of the domains. We show that the peptides tend to occupy the boundaries of liquid-ordered domains and significantly increase the energy barrier of the domain-domain fusion, which might lead to misregulation of raft clustering and adverse consequences for normal cell processes.

Additional contributions to elastic energy of lipid membranes: Tilt-curvature coupling and curvature gradient

Lipid bilayer membranes under biologically relevant conditions are flexible thin laterally fluid films consisting of two unimolecular layers (monolayers) each about 2 nm thick. On spatial scales much larger than the bilayer thickness, the membrane elasticity is well determined by its shape. The classical Helfrich theory considers the membrane as an elastic two-dimensional (2D) film, which has no particular internal structure. However, various local membrane heterogeneities can result in a lipids tilt relative to the membrane surface normal. On the basis of the classical elasticity theory of 3D bodies, Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)10.1007/s101890070003] derived the most general energy functional, taking into account the tilt and lipid monolayer curvature. Recently, Terzi and Deserno [J. Chem. Phys. 147, 084702 (2017)10.1063/1.4990404] showed that Hamm and Kozlov's derivation was incomplete because the tilt-curvature coupling term had been missed. However, the energy functional derived by Terzi and Deserno appeared to be unstable, thereby being invalid for applications that require minimizations of the overall energy of deformations. Here, we derive a stable elastic energy functional, showing that the squared gradient of the curvature was missed in both of these works. This change in the energy functional arises from a more accurate consideration of the transverse shear deformation terms and their influence on the membrane stability. We also consider the influence of the prestress terms on the stability of the energy functional, and we show that it should be considered small and the effective Gaussian curvature should be neglected because of the stability requirements. We further generalize the theory, including the stretching-compressing deformation modes, and we provide the geometrical interpretation of the terms that were previously missed by Hamm and Kozlov. The physical consequences of the new terms are analyzed in the case of a membrane-mediated interaction of two amphipathic peptides located in the same monolayer. We also provide the expression for director fluctuations, comparing it with that obtained by Terzi and Deserno.

Interleaflet Coupling of Lipid Nanodomains - Insights From in vitro Systems

The plasma membrane is a complex system, consisting of two layers of lipids and proteins compartmentalized into small structures called nanodomains. Despite the asymmetric composition of both leaflets, coupling between the layers is surprisingly strong. This can be evidenced, for example, by recent experimental studies performed on phospholipid giant unilamellar vesicles showing that nanodomains formed in the outer layer are perfectly registered with those in the inner leaflet. Similarly, microscopic phase separation in one leaflet can induce phase separation in the opposing leaflet that would otherwise be homogeneous. In this review, we summarize the current theoretical and experimental knowledge that led to the current view that domains are – irrespective of their size – commonly registered across the bilayer. Mechanisms inducing registration of nanodomains suggested by theory and calculations are discussed. Furthermore, domain coupling is evidenced by experimental studies based on the sparse number of methods that can resolve registered from independent nanodomains. Finally, implications that those findings using model membrane studies might have for cellular membranes are discussed.

Ordered Lipid Domains Assemble via Concerted Recruitment of Constituents from Both Membrane Leaflets

Lipid rafts serve as anchoring platforms for membrane proteins. Thus far they escaped direct observation by light microscopy due to their small size. Here we used differently colored dyes as reporters for the registration of both ordered and disordered lipids from the two leaves of a freestanding bilayer. Photoswitchable lipids dissolved or reformed the domains. Measurements of domain mobility indicated the presence of 120 nm wide ordered and 40 nm wide disordered domains. These sizes are in line with the predicted roles of line tension and membrane undulation as driving forces for alignment.

Elastic deformations mediate interaction of the raft boundary with membrane inclusions leading to their effective lateral sorting

Liquid-ordered lipid domains represent a lateral inhomogeneity in cellular membranes. These domains have elastic and physicochemical properties different from those of the surrounding membrane. In particular, their thickness exceeds that of the disordered membrane. Thus, elastic deformations arise at the domain boundary in order to compensate for the thickness mismatch. In equilibrium, the deformations lead to an incomplete register of monolayer ordered domains: the elastic energy is minimal if domains in opposing monolayers lie on the top of each other, and their boundaries are laterally shifted by about 3 nm. This configuration introduces a region, composed of one ordered and one disordered monolayers, with an intermediate bilayer thickness. Besides, a jump in a local monolayer curvature takes place in this intermediate region, concentrating here most of the elastic stress. This region can participate in a lateral sorting of membrane inclusions by offering them an optimal bilayer thickness and local curvature conditions. In the present study, we consider the sorting of deformable lipid inclusions, undeformable peripheral and deeply incorporated peptide inclusions, and undeformable transmembrane inclusions of different molecular geometry. With rare exceptions, all types of inclusions have an affinity to the ordered domain boundary as compared to the bulk phases. The optimal lateral distribution of inclusions allows relaxing the elastic stress at the boundary of domains.

The Effect of Transmembrane Protein Shape on Surrounding Lipid Domain Formation by Wetting

Signal transduction through cellular membranes requires the highly specific and coordinated work of specialized proteins. Proper functioning of these proteins is provided by an interplay between them and the lipid environment. Liquid-ordered lipid domains are believed to be important players here, however, it is still unclear whether conditions for a phase separation required for lipid domain formation exist in cellular membranes. Moreover, membrane leaflets are compositionally asymmetric, that could be an obstacle for the formation of symmetric domains spanning the lipid bilayer. We theoretically show that the presence of protein in the membrane leads to the formation of a stable liquid-ordered lipid phase around it by the mechanism of protein wetting by lipids, even in the absence of conditions necessary for the global phase separation in the membrane. Moreover, we show that protein shape plays a crucial role in this process, and protein conformational rearrangement can lead to changes in the size and characteristics of surrounding lipid domains.

Spontaneous and Stress-Induced Pore Formation in Membranes: Theory, Experiments and Simulations

The large plasticity, dynamics and adaptability of biological membranes allow different modes of intrinsic and inducible permeability. These phenomena are of physiological importance for a number of natural functions related to cell death and can also be manipulated artificially for practical purposes like gene transfer, drug delivery, prevention of infections or anticancer therapy. For these advances to develop in a controllable and specific way, we need a sufficient understanding of the membrane permeability phenomena. Since the formulation of early concepts of pore formation, there has been an enormous effort to describe membrane permeability by using theory, simulations and experiments. A major breakthrough has come recently through theoretical developments that allow building continuous trajectories of pore formation both in the absence and presence of stress conditions. The new model provides a coherent quantitative view of membrane permeabilization, useful to test the impact of known lipid properties, make predictions and postulate specific pore intermediates that can be studied by simulations. For example, this theory predicts unprecedented dependencies of the line tension on the pore radius and on applied lateral tension which explain previous puzzling results. In parallel, important concepts have also come from molecular dynamics simulations, of which the role of water for membrane permeabilization is of special interest. These advances open new challenges and perspectives for future progress in the study of membrane permeability, as experiments and simulations will need to test the theoretical predictions, while theory achieves new refinements that provide a physical ground for observations.

Experimental Evidence of the Existence of Interleaflet Coupled Nanodomains: An MC-FRET Study

Plasma membranes of living cells are compartmentalized into small submicroscopic structures (nanodomains) having potentially relevant biological functions. Despite this, structural features of these nanodomains remain elusive, primarily due to the difficulties in characterizing such small dynamic entities. It is unclear whether nanodomains found in the upper bilayer leaflet are transversally registered with those found in the lower leaflet. Experiments performed on larger microscopic domains indicate that the coupling between the leaflets is strong, forcing the domains to be in perfect registration, but can the same thing be said about the biologically more relevant nanodomains? This work provides experimental evidence that even small nanodomains of variable sizes between 10 and 160 nm are interleaflet coupled. Importantly, the alternative scenarios of partially registered, independent, or antiregistered nanodomains could be excluded.

Membrane Elastic Deformations Modulate Gramicidin A Transbilayer Dimerization and Lateral Clustering

Gramicidin A (gA) is a short β-helical peptide known to form conducting channels in lipid membranes because of transbilayer dimerization. The gA conducting dimer, being shorter than the lipid bilayer thickness, deforms the membrane in its vicinity, and the bilayer elastic energy contributes to the gA dimer formation energy. Likewise, membrane incorporation of a gA monomer, which is shorter than the lipid monolayer thickness, creates a void, thereby forcing surrounding lipid molecules to tilt to fill it. The energy of membrane deformation was calculated in the framework of the continuum elasticity theory, taking into account splay, tilt, lateral stretching/compression, Gaussian splay deformations, and external membrane tension. We obtained the interaction energy profiles for two gA monomers located either in the same or in the opposite monolayers. The profiles demonstrated the long-range attraction and short-range repulsion behavior of the monomers resulting from the membrane deformation. Analysis of the profile features revealed conditions under which clusters of gA monomers would not dissipate because of diffusion. The calculated dependence of the dimer formation and decay energy barriers on the membrane elastic properties was in good agreement with the available experimental data and suggested an explanation for a hitherto contentious phenomenon.

Lateral Membrane Heterogeneity Regulates Viral-Induced Membrane Fusion during HIV Entry

Sphingomyelin- and cholesterol- enriched membrane domains, commonly referred to as “rafts” play a crucial role in a large number of intra- and intercellular processes. Recent experiments suggest that not only the volumetric inhomogeneity of lipid distribution in rafts, but also the arrangement of the 1D boundary between the raft and the surrounding membrane is important for the membrane-associated processes. The reason is that the boundary preferentially recruits different peptides, such as HIV (human immunodeficiency virus) fusion peptide. In the present work, we report a theoretical investigation of mechanisms of influence of the raft boundary arrangement upon virus-induced membrane fusion. We theoretically predict that the raft boundary can act as an attractor for viral fusion peptides, which preferentially distribute into the vicinity of the boundary, playing the role of ‘line active components’ of the membrane (‘linactants’). We have calculated the height of the fusion energy barrier and demonstrated that, in the case of fusion between HIV membrane and the target cell, presence of the raft boundary in the vicinity of the fusion site facilitates fusion. The results we obtained can be further generalized to be applicable to other enveloped viruses.

Composition Fluctuations in Lipid Bilayers

Cell membranes contain multiple lipid and protein components having heterogeneous in-plane (lateral) distribution. Nanoscale rafts are believed to play an important functional role, but their phase state—domains of coexisting phases or composition fluctuations—is unknown. As a step toward understanding lateral organization of cell membranes, we investigate the difference between nanoscale domains of coexisting phases and composition fluctuations in lipid bilayers. We simulate model lipid bilayers with the MARTINI coarse-grained force field on length scales of tens of nanometers and timescales of tens of microseconds. We use a binary and a ternary mixture: a saturated and an unsaturated lipid, or a saturated lipid, an unsaturated lipid, and cholesterol, respectively. In these mixtures, the phase behavior can be tuned from a mixed state to a coexistence of a liquid-crystalline and a gel, or a liquid-ordered and a liquid-disordered phase. Transition from a two-phase to a one-phase state is achieved by raising the temperature and adding a hybrid lipid (with a saturated and an unsaturated chain). We analyze the evolution of bilayer properties along this transition: domains of two phases transform to fluctuations with local ordering and compositional demixing. Nanoscale domains and fluctuations differ in several properties, including interleaflet overlap and boundary length. Hybrid lipids show no enrichment at the boundary, but decrease the difference between the coexisting phases by ordering the disordered phase, which could explain their role in cell membranes.

Switching between Successful and Dead-End Intermediates in Membrane Fusion

Fusion of cellular membranes during normal biological processes, including proliferation, or synaptic transmission, is mediated and controlled by sophisticated protein machinery ensuring the preservation of the vital barrier function of the membrane throughout the process. Fusion of virus particles with host cell membranes is more sparingly arranged and often mediated by a single fusion protein, and the virus can afford to be less discriminative towards the possible different outcomes of fusion attempts. Formation of leaky intermediates was recently observed in some fusion processes, and an alternative trajectory of the process involving formation of π-shaped structures was suggested. In this study, we apply the methods of elasticity theory and Lagrangian formalism augmented by phenomenological and molecular geometry constraints and boundary conditions to investigate the traits of this trajectory and the drivers behind the choice of one of the possible scenarios depending on the properties of the system. The alternative pathway proved to be a dead end, and, depending on the parameters of the participating membranes and fusion proteins, the system can either reversibly enter the corresponding “leaky” configuration or be trapped in it. A parametric study in the biologically relevant range of variables emphasized the fusion protein properties crucial for the choice of the fusion scenario.

Pore formation in lipid membrane II: Energy landscape under external stress

Lipid membranes are extremely stable envelopes allowing cells to survive in various environments and to maintain desired internal composition. Membrane permeation through formation of transversal pores requires substantial external stress. Practically, pores are usually formed by application of lateral tension or transmembrane voltage. Using the same approach as was used for obtaining continuous trajectory of pore formation in the stress-less membrane in the previous article, we now consider the process of pore formation under the external stress. The waiting time to pore formation proved a non-monotonous function of the lateral tension, dropping from infinity at zero tension to a minimum at the tension of several millinewtons per meter. Transmembrane voltage, on the contrary, caused the waiting time to decrease monotonously. Analysis of pore formation trajectories for several lipid species with different spontaneous curvatures and elastic moduli under various external conditions provided instrumental insights into the mechanisms underlying some experimentally observed phenomena.

Membrane phase transition during heating and cooling: molecular insight into reversible melting

With increasing temperature, lipid bilayers undergo a gel-fluid phase transition, which plays an essential role in many physiological phenomena. In the present work, this first-order phase transition was investigated for variable heating and cooling rates for a dipalmitoylphosphatidylcholine (DPPC) lipid bilayer by means of atomistic molecular dynamics simulations. Alternative methods to track the melting temperature [Formula: see text] are compared. The resulting [Formula: see text] is shown to be independent of the scan rate for small heating rates (0.05-0.3 K/ns) implying reversible melting, and increases for larger heating (0.3-4 K/ns) or cooling rates (2-0.1 K/ns). The reported dependency of the melting temperature on the heating rate is in perfect agreement with a two-state kinetic rate model as suggested previously. Expansion and shrinkage, as well as the dynamics of melting seeds is described. The simulations further exhibit a relative shift between melting seeds in opposing membrane leaflets as predicted from continuum elastic theory.

Lipid Domain Co-localization Induced by Membrane Undulations

Multicomponent lipid bilayer membranes display rich phase transition and associated compositional lipid domain formation behavior. When both leaflets of the bilayer contain domains, they are often found co-localized across the leaflets, implying the presence of a thermodynamic interleaflet coupling. In this work, it is demonstrated that fluctuation-induced interactions between domains embedded within opposing membrane leaflets provide a robust means to co-localize the domains. In particular, it is shown via a combination of a mode-counting argument, a perturbative calculation, and a non-perturbative treatment of a special case, that spatial variations in membrane bending rigidity associated with lipid domains embedded within the background phase always lead to an attractive interleaflet coupling with a magnitude of ∼0.01k_BT/nm² in simple model membrane systems. Finally, it is demonstrated that the fluctuation-induced coupling is very robust against membrane tension and substrate interactions.

Undulations Drive Domain Registration from the Two Membrane Leaflets

Phase separation in biological membranes plays an important role in protein targeting and transmembrane signaling. Its occurrence in both membrane leaflets commonly gives rise to matching liquid or liquid-ordered domains in the opposing monolayers. The underlying mechanism of such co-localization is not fully understood. The decrease of the line tension around the thicker ordered domain constitutes an important driving force. Yet, robust domain coupling requires an additional energy source, which we have now identified as thermal undulations. Our theoretical analysis of elastic deformations in a lipid bilayer shows that stiffer lipid domains tend to distribute into areas with lower fluctuations of monolayer curvature. These areas naturally align in the opposing monolayers. Thus, coupling requires both membrane leafs to display a heterogeneity in splay rigidities. The heterogeneity may either originate from intrinsic lipid properties or be acquired by adsorption of peripheral molecules. Undulations and line tension act synergistically: the gain in energy due a minimized line tension is proportional to domain radius and thus primarily fuels the registration of smaller domains; whereas the energetic contribution of undulations increases with membrane area and thus primarily acts to coalesce larger domains.

Line Activity of Ganglioside GM1 Regulates the Raft Size Distribution in a Cholesterol-Dependent Manner

Liquid-ordered lipid domains, also called rafts, are assumed to be important players in different cellular processes, mainly signal transduction and membrane trafficking. They are thicker than the disordered part of the membrane and are thought to form to compensate for the hydrophobic mismatch between transmembrane proteins and the lipid environment. Despite the existence of such structures in vivo still being an open question, they are observed in model systems of multicomponent lipid bilayers. Moreover, the predictions obtained from model experiments allow the explanation of different physiological processes possibly involving rafts. Here we present the results of the study of the regulation of raft size distribution by ganglioside GM1. Combining atomic force microscopy with theoretical considerations based on the theory of membrane elasticity, we predict that this glycolipid should change the line tension of raft boundaries in two different ways, mainly depending on the cholesterol content. These results explain the shedding of gangliosides from the surface of tumor cells and the following ganglioside-induced apoptosis of T-lymphocytes in a raft-dependent manner. Moreover, the generality of the model allows the prediction of the line activity of different membrane components based on their molecular geometry.

Molecular Insight into the Line Tension of Bilayer Membranes Containing Hybrid Polyunsaturated Lipids

Line tension (γ) is a key parameter for the structure and dynamics of membrane domains. It was proposed that hybrid lipids, with mixed saturated and unsaturated acyl chains, participate in the relaxation of γ through different mechanisms. In this work, we used molecular dynamics simulations of the coarse-grained MARTINI model to measure γ in liquid-ordered-liquid-disordered (Lo-Ld) membranes, with increasingly larger relative proportion of the hybrid polyunsaturated lipid PAPC (4:0-5:4PC) to DAPC (di5:4PC) (i.e., X_H). We also calculated an elastic contribution to γ by the Lo-Ld thickness mismatch, tilt moduli, and bending moduli, as predicted by theory. We found that an increase in X_H decreased the overall γ value and the elastic contribution to line tension. The effect on the elastic line tension is driven by a reduced hydrophobic mismatch. Changes in the elastic constants of the phases due to an increase in X_H produced a slightly larger elastic γ term. In addition to this elastic energy, other major contributions to γ are found in these model membranes. Increasing X_H decreases both elastic and nonelastic contributions to γ. Finally, PAPC also behaves as a linactant, relaxing γ through an interfacial effect, as predicted by theoretical results. This study gives insight into the actual contribution of distinct energy terms to γ in bilayers containing polyunsaturated hybrid lipids.

Interleaflet Coupling, Pinning, and Leaflet Asymmetry-Major Players in Plasma Membrane Nanodomain Formation

The plasma membrane has a highly asymmetric distribution of lipids and contains dynamic nanodomains many of which are liquid entities surrounded by a second, slightly different, liquid environment. Contributing to the dynamics is a continuous repartitioning of components between the two types of liquids and transient links between lipids and proteins, both to extracellular matrix and cytoplasmic components, that temporarily pin membrane constituents. This make plasma membrane nanodomains exceptionally challenging to study and much of what is known about membrane domains has been deduced from studies on model membranes at equilibrium. However, living cells are by definition not at equilibrium and lipids are distributed asymmetrically with inositol phospholipids, phosphatidylethanolamines and phosphatidylserines confined mostly to the inner leaflet and glyco- and sphingolipids to the outer leaflet. Moreover, each phospholipid group encompasses a wealth of species with different acyl chain combinations whose lateral distribution is heterogeneous. It is becoming increasingly clear that asymmetry and pinning play important roles in plasma membrane nanodomain formation and coupling between the two lipid monolayers. How asymmetry, pinning, and interdigitation contribute to the plasma membrane organization is only beginning to be unraveled and here we discuss their roles and interdependence.

Helix-helix interactions in membrane domains of bitopic proteins: Specificity and role of lipid environment

Interaction between transmembrane helices often determines biological activity of membrane proteins. Bitopic proteins, a broad subclass of membrane proteins, form dimers containing two membrane-spanning helices. Some aspects of their structure-function relationship cannot be fully understood without considering the protein-lipid interaction, which can determine the protein conformational ensemble. Experimental and computer modeling data concerning transmembrane parts of bitopic proteins are reviewed in the present paper. They highlight the importance of lipid-protein interactions and resolve certain paradoxes in the behavior of such proteins. Besides, some properties of membrane organization provided a clue to understanding of allosteric interactions between distant parts of proteins. Interactions of these kinds appear to underlie a signaling mechanism, which could be widely employed in the functioning of many membrane proteins. Treatment of membrane proteins as parts of integrated fine-tuned proteolipid system promises new insights into biological function mechanisms and approaches to drug design. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Analytical calculation of the lipid bilayer bending modulus

Bending and Gaussian moduli of a homogenious single-component lipid bilayer are calculated analytically using microscopic model of the lipid hydrocarbon chains. The approach allows for thermodynamic averaging over different chains conformations. Each chain is modeled as a flexible string with finite bending rigidity and an incompressible cross-section area. The interchain steric repulsion is accounted for self-consistently determined single-chain confining parabolic potential. The model provides a simple analytical expression for the membrane bending modulus, which falls within a range of experimental values. An observed dependence of the modulus on hydrocarbon chain length is also reproduced. Correspondence between our microscopic model and the membrane theory of elasticity is established.

The aliphatic chain of cholesterol modulates bilayer interleaflet coupling and domain registration

Cholesterol is a necessary component and critical regulator of liquid-ordered membrane domains. However, the structural features that determine its unique physicochemical behaviors are not fully understood. In particular, very little is known about the specific functions of the terminal aliphatic chain of cholesterol, as previous studies have focused mainly on the rigid sterol ring structure and its hydroxyl head. In the current work, we used coarse-grained molecular dynamics simulations to investigate the effect of cholesterol aliphatic chain length on the dynamics and structure of coexisting lipid domains. We found that the aliphatic chain has no appreciable effect on phase separation per se, but it significantly affects the rate of cholesterol flip-flop and intermonolayer interaction. These effects are accompanied by changes in domain dynamics, lateral pressure, and interleaflet coupling. Our study provides useful insight into how biological sterols modulate communication between the outer and inner surfaces of the plasma membrane and, therefore, cellular signaling.

Transbilayer Colocalization of Lipid Domains Explained via Measurement of Strong Coupling Parameters

When micron-scale compositional heterogeneity develops in membranes, the distribution of lipids on one face of the membrane strongly affects the distribution on the other. Specifically, when lipid membranes phase separate into coexisting liquid phases, domains in each monolayer leaflet of the membrane are colocalized with domains in the opposite leaflet. Colocalized domains have never been observed to spontaneously move out of registry. This result indicates that the lipid compositions in one leaflet are strongly coupled to compositions in the opposing leaflet. Predictions of the interleaflet coupling parameter, Λ, vary by a factor of 50. We measure the value of Λ by applying high shear forces to supported lipid bilayers. This causes the upper leaflet to slide over the lower leaflet, moving domains out of registry. We find that the threshold shear stress required to deregister domains in the upper and lower leaflets increases with the inverse length of domains. We derive a simple, closed-form expression relating the threshold shear to Λ, and find Λ = 0.016 ± 0.004 kBT/nm2.

Roles of Interleaflet Coupling and Hydrophobic Mismatch in Lipid Membrane Phase-Separation Kinetics

Characterizing the nanoscale dynamic organization within lipid bilayer membranes is central to our understanding of cell membranes at a molecular level. We investigate phase separation and communication across leaflets in ternary lipid bilayers, including saturated lipids with between 12 and 20 carbons per tail. Coarse-grained molecular dynamics simulations reveal a novel two-step kinetics due to hydrophobic mismatch, in which the initial response of the apposed leaflets upon quenching is to increase local asymmetry (antiregistration), followed by dominance of symmetry (registration) as the bilayer equilibrates. Antiregistration can become thermodynamically preferred if domain size is restricted below ∼20 nm, with implications for the symmetry of rafts and nanoclusters in cell membranes, which have similar reported sizes. We relate our findings to theory derived from a semimicroscopic model in which the leaflets experience a “direct” area-dependent coupling, and an “indirect” coupling that arises from hydrophobic mismatch and is most important at domain boundaries. Registered phases differ in composition from antiregistered phases, consistent with a direct coupling between the leaflets. Increased hydrophobic mismatch purifies the phases, suggesting that it contributes to the molecule-level lipid immiscibility. Our results demonstrate an interplay of competing interleaflet couplings that affect phase compositions and kinetics, and lead to a length scale that can influence lateral and transverse bilayer organization within cells.

Elastic deformations of bolalipid membranes

Archaeal membranes have unique mechanical properties that enable these organisms to survive under extremely aggressive environmental conditions. The so-called bolalipids contribute to this exceptional stability. They have two polar heads joined by two hydrocarbon chains. The two headgroups can face different sides of the membrane (O-shape conformation) or the same side (U-shape conformation). We have developed an elasticity theory for bolalipid membranes and show that the energetic contributions of (i) tilt deformations, (ii) area compression/stretching deformations, (iii) as well as those of Gaussian splay from the two membrane surfaces are additive, while splay deformations yield a cross-term. The presence of a small fraction of U-shaped molecules resulted in spontaneous membrane curvature. We estimated the tilt modulus to be approximately equal to that of membranes in eukaryotic cells. In contrast to conventional lipids, the bolalipid membrane possesses two splay moduli, one of which is estimated to be an order of magnitude larger than that of conventional lipids. The projected values of elastic moduli act to hamper pore formation and to decelerate membrane fusion and fission.

Metabolic Precursor of Cholesterol Causes Formation of Chained Aggregates of Liquid-Ordered Domains

7-Dehydrocholesterol, an immediate metabolic predecessor of cholesterol, can accumulate in tissues due to some metabolic abnormalities, causing an array of symptoms known as Smith-Lemli-Opitz syndrome. Enrichment of cellular membranes with 7-dehydrocholesterol interferes with normal cell-signaling processes, which involve interaction between rafts and formation of the so-called signaling platforms. In model membranes, cholesterol-based ordered domains usually merge upon contact. According to our experimental data, ordered domains in the model systems where cholesterol is substituted for 7-dehydrocholesterol never merge on the time scale of the experiment, but clusterize into necklace-like aggregates. We attribute such different dynamical behavior to altered properties of the domain boundary. In the framework of thickness mismatch model, we analyzed changes of interaction energy profiles of two approaching domains caused by substitution of cholesterol by 7-dehydrocholesterol. The energy barrier for domain merger is shown to increase notably, with simultaneous appearance of another distinct local energy minimum. Such energy profile is in perfect qualitative agreement with the experimental observations. The observed change of domain dynamics can impair proper interaction between cellular rafts underlying pathologies associated with deviations in cholesterol metabolism.