Molecular mechanisms of spontaneous curvature and softening in complex lipid bilayer mixtures

Softening in Two-Component Lipid Mixtures by Spontaneous Curvature Variance

The bending modulus of a lipid bilayer quantifies its mechanical resistance to curvature. It is typically understood in terms of thickness; e.g., thicker bilayers are usually stiffer. Here, we describe an additional and powerful molecular determinant of stiffness─the variance in the distribution of curvature sensitivity of lipids and lipid conformations. Zwitterionic choline and ethanolamine headgroups of glycerophospholipids dynamically explore inter- and intraspecies interactions, leading to transient clustering. We demonstrate that these clusters couple strongly to negative curvature, exciting undulatory membrane modes and reducing the apparent bending modulus. Three force fields (Martini 2, Martini 3, and all-atom CHARMM C36) each show the effect to a different extent, with the coarse-grained Martini models showing the most clustering and thus the most softening. The theory is a guide to understanding the stiffness of biological membranes with their complex composition, as well as how choices of force field parameterization are translated into mechanical stiffness.

The Role of Lipid Intrinsic Curvature in the Droplet Interface Bilayer

Model bilayers are constructed from lipids having different intrinsic curvatures using the droplet interface bilayer (DIB) method, and their static physicochemical properties are determined. Geometrical and tensiometric measurements are used to derive the free energy of formation (ΔF) of a two-droplet DIB relative to a pair of isolated aqueous droplets, each decorated with a phospholipid monolayer. The lipid molecules employed have different headgroup sizes but identical hydrophobic tail structure, and each is characterized by an intrinsic curvature value (c0) that increases in absolute value with decreasing size of headgroup. Mixtures of lipids at different ratios were also investigated. The role of curvature stress on the values of ΔF of the respective lipid bilayers in these model membranes is discussed and is illuminated by the observation of a decrement in ΔF that scales as a near linear function of c02. Overall, the results reveal an association that should prove useful in studies of ion channels and other membrane proteins embedded in model droplet bilayer systems that will impact the understanding of protein function in cellular membranes composed of lipids of high and low curvature.

Computing the influence of lipids and lipid complexes on membrane mechanics

Methodology for extracting the spontaneous curvature, bending modulus, and neutral surface of a lipid bilayer is described. The "SPEX" method is a robust technique for computing the bilayer bending modulus while allowing for resolution of the spontaneous curvature of specific interacting lipids and complexes, and the dependence of spontaneous curvature on wavelength. The method is described referring to the publicly available MembraneAnalysis.jl software package.

Building complex membranes with Martini 3

The Martini model is a popular force field for coarse-grained simulations. Membranes have always been at the center of its development, with the latest version, Martini 3, showing great promise in capturing more and more realistic behavior. In this chapter we provide a step-by-step tutorial on how to construct starting configurations, run initial simulations and perform dedicated analysis for membrane-based systems of increasing complexity, including leaflet asymmetry, curvature gradients and embedding of membrane proteins.

Membranes get in shape: Biophysics of curving bilayers

Membrane curvature is ubiquitous and essential in cell biology. Curved membranes have several distinct features, including specific protein and lipid sorting, distinct lipid ordering, and changes in transbilayer stress. Curvature also interplays with membrane tension to generate forces that change membrane shape. This research highlight summarizes recent contributions to this topic published in Biophysical Journal.

LEGO-lipophosphonoxins: length of hydrophobic module affects permeabilizing activity in target membranes of different phospholipid composition

In the past few decades, society has faced rapid development and spreading of antimicrobial resistance due to antibiotic misuse and overuse and the immense adaptability of bacteria. Difficulties in obtaining effective antimicrobial molecules from natural sources challenged scientists to develop synthetic molecules with antimicrobial effect. We developed modular molecules named LEGO-Lipophosphonoxins (LEGO-LPPO) capable of inducing cytoplasmic membrane perforation. In this structure-activity relationship study we focused on the role of the LEGO-LPPO hydrophobic module directing the molecule insertion into the cytoplasmic membrane. We selected three LEGO-LPPO molecules named C9, C8 and C7 differing in the length of their hydrophobic chain and consisting of an alkenyl group containing one double bond. The molecule with the long hydrophobic chain (C9) was shown to be the most effective with the lowest MIC and highest perforation rate both in vivo and in vitro. We observed high antimicrobial activity against both G+ and G- bacteria with significant differences in LEGO-LPPOs mechanism of action on these two cell types. We observed a highly cooperative mechanism of LEGO-LPPO action on G- bacteria as well as on liposomes resembling G- bacteria. LEGO-LPPO action on G- bacteria was significantly slower compared to G+ bacteria suggesting the role of the outer membrane in affecting the LEGO-LPPOs perforation rate. This notion was supported by the higher sensitivity of the E. coli strain with a compromised outer membrane. Finally, we noted that the composition of the cytoplasmic membrane affects the activity of LEGO-LPPOs since the presence of phosphatidylethanolamine increases their membrane disrupting activity.

Origin of the nonlinear structural and mechanical properties in oppositely curved lipid mixtures

Structural and mechanical properties of membranes such as thickness, tail order, bending modulus and curvature energetics play crucial role in controlling various cellular functions that depend on the local lipid organization and membrane reshaping. While behavior of these biophysical properties are well understood in single component membranes, very little is known about how do they change in the mixed lipid membranes. Often various properties of the mixed lipid bilayers are assumed to change linearly with the mole fractions of the constituent lipids which, however, is true for "ideal" mixing only. In this study, using molecular dynamics simulations, we show that structural and mechanical properties of binary lipid mixture change nonlinearly with the lipid mole fractions, and the strength of the nonlinearity depends on two factors - spontaneous curvature difference and locally inhomogeneous interactions between the lipid components.

Membrane free-energy landscapes derived from atomistic dynamics explain nonuniversal cholesterol-induced stiffening

All lipid membranes have inherent morphological preferences and resist deformation. Yet adaptations in membrane shape can and do occur at multiple length scales. While this plasticity is crucial for cellular physiology, the factors controlling the morphological energetics of lipid bilayers and the dominant mechanisms of membrane remodeling remain to be fully understood. An ongoing debate regarding the universality of the stiffening effect of cholesterol underscores the challenges facing this field, both experimentally and theoretically, even for simple lipid mixtures. On the computational side, we have argued that enhanced-sampling all-atom molecular dynamics simulations are uniquely suited for the quantification of membrane conformational energetics, as they minimize a priori assumptions and permit analysis of bilayers in deformed states. To showcase this approach, we examine reported inconsistencies between alternative experimental measurements of bending moduli for cholesterol-enriched membranes. Specifically, we analyze lipid bilayers with different chain saturation and compute free-energy landscapes for curvature deformations distributed over areas from ∼5 to ∼60 nm2 . These enhanced simulations, totaling over 100 μs of sampling time, enable us to directly quantify both bending and tilt moduli and to dissect the contributing factors and molecular mechanisms of curvature generation at each length scale. Our results show that the effects of cholesterol on bending rigidity are lipid-specific and suggest that this specificity arises from differences in the torsional dynamics of the acyl chains. In summary, we demonstrate that quantitative relationships can now be established between lipid structure and bending energetics, paving the way for addressing open fundamental questions in cell membrane mechanics.

Simulated dynamic cholesterol redistribution favors membrane fusion pore constriction

Endo- and exocytosis proceed through a highly strained membrane fusion pore topology regardless of the aiding protein machinery. The membrane's lipid components bias fusion pores toward expansion or closure, modifying the necessary work done by proteins. Cholesterol, a key component of plasma membranes, promotes both inverted lipid phases with concave leaflets (i.e., negative total curvature, which thins the leaflet) and flat bilayer phases with thick, ordered hydrophobic interiors. We demonstrate by theory and simulation that both leaflets of nascent catenoidal fusion pores have negative total curvature. Furthermore, the hydrophobic core of bilayers with strong negative Gaussian curvature is thinned. Therefore, it is an open question whether cholesterol will be enriched in these regions because of the negative total curvature or depleted because of the membrane thinning. Here, we compare all-atom molecular dynamics simulations (built using a procedure to create specific fusion pore geometries) and theory to understand the underlying reasons for lipid redistribution on fusion pores. Our all-atom molecular dynamics simulations resolve this question by showing that cholesterol is strongly excluded from the thinned neck of fusion and fission pores, revealing that thickness (and/or lipid order) influences cholesterol distributions more than curvature. The results imply that cholesterol exclusion can drive fusion pore closure by creating a small, cholesterol-depleted zone in the neck. This model agrees with literature evidence that membrane reshaping is connected to cholesterol-dependent lateral phase separation.

Interaction of detergent with complex mimics of bacterial membranes

Detergents are valuable tools to extract membrane proteins for biophysical, biochemical, and structural scrutiny. The detergent-driven solubilization of bilayers made from a single lipid species is commonly described in terms of pseudo-phase diagrams and a three-stage model accounting for three ranges comprising (i) intact vesicles, (ii) vesicle/micelle co-existence, or (iii) mixed micelles. Moreover, the pseudo-phase boundaries thus determined can often be quantitatively rationalized in terms of the molecular shapes of the lipid and the detergent used. Yet, it has remained unclear to what extent this approach can be applied to multi-component lipid membranes that more closely mimic the compositional complexity of cellular membranes. Here, we studied how lipid mixtures composed of palmitoyl oleoyl phosphatidylethanolamine (POPE), palmitoyl oleoyl phosphatidylglycerol (POPG), and tetraoleoyl cardiolipin (TOCL) are solubilized by the commonly used zwitterionic detergent lauryldimethylamine N-oxide using isothermal titration calorimetry. While phase diagrams of the diverse lipid mixtures showed the typical ranges of the three-stage model, we found that POPG-rich POPE/POPG bilayers are more difficult to solubilize than POPG-poor POPE/POPG bilayers. In turn, POPE/POPG/TOCL bilayers became increasingly resistant to detergent with increasing TOCL content. Since POPG is nearly cylindrically shaped and TOCL adopts inverted cone-like shapes under current buffer conditions, our solubilization data do not align with shape-based arguments. Instead, additional electrostatic interactions between lipids and detergents lead to non-additive mixing behavior affecting the resilience of complex lipid bilayers against solubilization.

Kinetic relaxation of giant vesicles validates diffusional softening in a binary lipid mixture

The stiffness of biological membranes determines the work required by cellular machinery to form and dismantle vesicles and other lipidic shapes. Model membrane stiffness can be determined from the equilibrium distribution of giant unilamellar vesicle surface undulations observable by phase contrast microscopy. With two or more components, lateral fluctuations of composition will couple to surface undulations depending on the curvature sensitivity of the constituent lipids. The result is a broader distribution of undulations whose complete relaxation is partially determined by lipid diffusion. In this work, kinetic analysis of the undulations of giant unilamellar vesicles made of phosphatidylcholine-phosphatidylethanolamine mixtures validates the molecular mechanism by which the membrane is made 25% softer than a single-component one. The mechanism is relevant to biological membranes, which have diverse and curvature-sensitive lipids.

The full model of micropipette aspiration of cells: A mesoscopic simulation

Studies on the interaction between cells and micromanipulation tools are necessary to optimize the procedures and improve the developmental potential of cells. The molecular dynamics simulation is not possible for such a large-scale simulation, and the spring-damped viscoelastic models and the constitutive equations of the continuum are usually adopted to model the cells as a whole without consideration of the different properties presented by the heterogeneous subcellular components. In this study, we utilized coarse-grained modeling to develop a subcellular model of suspension cell dynamics and a model of a holding micropipette for the fixation of a suspension cell, and designed a large-scale, accurate mesoscopic simulation environment for specific cell micromanipulation. We established a triangular mesh cell membrane and a uniformly distributed, non-intersecting cytoskeleton network and added polymerization/depolymerization processes to connect the cytoskeleton chains with the membrane and cross-linking proteins. In the cell aspiration model, we adopted the profile of the reversed Poiseuille flow to calibrate the viscosity of the fluid and set the bounce-back condition and the appropriate solid-fluid force coefficient to realize non-slip flow at the boundary. The rheological properties of the cells during micropipette aspiration were further analyzed in the simulation by varying parameters such as the inner diameter of the micropipette, negative pressure, and maximum bond length. The model well reproduced the experimentally observed cell deformation phenomenon at low and high pressures. The dynamic response of the cell elongation observed from the simulation was consistent with that obtained from the analysis of the experimental data collected from a custom-designed micromanipulation system. STATEMENT OF SIGNIFICANCE: In this study, we extended the coarse-grained modeling of cells by developing a relatively large-scale micromanipulation environment consisting of a subcellular cell dynamics model and a fluid flow model for cell aspiration. We simulated cytoskeleton filaments that were uniformly distributed in space via applying Harmonic energy to model cytoskeleton with a high level of fidelity. The shortcoming of the soft repulsion in the solid-fluid interaction in the current simulation technique was solved by implementing the bounce-back boundary and the condition that the total force imposed by the wall particles on the fluid particles was equal to the pressure of the fluid. This work paved the way for understanding the mechanical properties of cells and improving the biological efficacy of micromanipulation.