A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited

Neutron spin echo shows pHLIP is capable of retarding membrane thickness fluctuations

Cell membranes are responsible for a range of biological processes that require interactions between lipids and proteins. While the effects of lipids on proteins are becoming better understood, our knowledge of how protein conformational changes influence membrane dynamics remains rudimentary. Here, we performed experiments and computer simulations to study the dynamic response of a lipid membrane to changes in the conformational state of pH-low insertion peptide (pHLIP), which transitions from a surface-associated (SA) state at neutral or basic pH to a transmembrane (TM) α-helix under acidic conditions. Our results show that TM-pHLIP significantly slows down membrane thickness fluctuations due to an increase in effective membrane viscosity. Our findings suggest a possible membrane regulatory mechanism, where the TM helix affects lipid chain conformations, and subsequently alters membrane fluctuations and viscosity.

Effects of Normal and Lateral Electric Fields on Membrane Mechanical Properties

As a core component of biological and synthetic membranes, lipid bilayers are key to compartmentalizing chemical processes. Bilayer morphology and mechanical properties are heavily influenced by electric fields, such as those caused by biological ion concentration gradients. We present atomistic simulations exploring the effects of electric fields applied normally and laterally to lipid bilayers. We find that normal fields decrease membrane tension, while lateral fields increase it. Free energy perturbation calculations indicate the importance of dipole-dipole interactions to these tension changes, especially for lateral fields. We additionally show that membrane area compressibilities can be related to their cohesive energies, allowing us to estimate changes in membrane bending rigidity under applied fields. We find that normal and lateral fields decrease and increase bending rigidity, respectively. These results point to the use of directed electric fields to locally control membrane stiffness, thereby modulating associated cellular processes.

Exploration of lipid bilayer mechanical properties using molecular dynamics simulation

Important biological structures known for their exceptional mechanical qualities, lipid bilayers are essential to many cellular functions. Fluidity, elasticity, permeability, stiffness, tensile strength, compressibility, shear viscosity, line tension, and curvature elasticity are some of the fundamental characteristics affecting their behavior. The purpose of this review is to examine these characteristics in more detail by molecular dynamics simulation, elucidating their importance and the elements that lead to their appearance in lipid bilayers. Comprehending the mechanical characteristics of lipid bilayers is critical for creating medications, drug delivery systems, and biomaterials that interact with biological membranes because it allows one to understand how these materials respond to different stresses and deformations. The influence of mechanical characteristics on important lipid bilayer properties is examined in this review. The mechanical properties of lipid bilayers were clarified through the use of molecular dynamics simulation analysis techniques, including bilayer thickness, stress-strain analysis, lipid bilayer area compressibility, membrane bending rigidity, and time- or ensemble-averaged the area per lipid evaluation. We explain the significance of molecular dynamics simulation analysis methods, providing important new information about the stability and dynamic behavior of the bilayer. In the end, we hope to use molecular dynamics simulation to provide a comprehensive understanding of the mechanical properties and behavior of lipid bilayers, laying the groundwork for further studies and applications. Taken together, careful investigation of these mechanical aspects deepens our understanding of the adaptive capacities and functional roles of lipid bilayers in biological environments.

Cryo-EM images of phase-separated lipid bilayer vesicles analyzed with a machine-learning approach

Lateral lipid heterogeneity (i.e., raft formation) in biomembranes plays a functional role in living cells. Three-component mixtures of low- and high-melting lipids plus cholesterol offer a simplified experimental model for raft domains in which a liquid-disordered (Ld) phase coexists with a liquid-ordered (Lo) phase. Using such models, we recently showed that cryogenic electron microscopy (cryo-EM) can detect phase separation in lipid vesicles based on differences in bilayer thickness. However, the considerable noise within cryo-EM data poses a significant challenge for accurately determining the membrane phase state at high spatial resolution. To this end, we have developed an image-processing pipeline that utilizes machine learning (ML) to predict the bilayer phase in projection images of lipid vesicles. Importantly, the ML method exploits differences in both the thickness and molecular density of Lo compared to Ld, which leads to improved phase identification. To assess accuracy, we used artificial images of phase-separated lipid vesicles generated from all-atom molecular dynamics simulations of Lo and Ld phases. Synthetic ground-truth data sets mimicking a series of compositions along a tieline of Ld + Lo coexistence were created and then analyzed with various ML models. For all tieline compositions, we find that the ML approach can correctly identify the bilayer phase with >90% accuracy, thus providing a means to isolate the intensity profiles of coexisting Ld and Lo phases, as well as accurately determine domain-size distributions, number of domains, and phase-area fractions. The method described here provides a framework for characterizing nanoscopic lateral heterogeneities in membranes and paves the way for a more detailed understanding of raft properties in biological contexts.

Effect of leaflet asymmetry on the stretching elasticity of lipid bilayers with phosphatidic acid

The asymmetry of membranes has a significant impact on their biophysical characteristics and behavior. This study investigates the composition and mechanical properties of symmetric and asymmetric membranes in giant unilamellar vesicles (GUVs) made of palmitoyloleoyl phosphatidylcholine (POPC) and palmitoyloleoyl phosphatidic acid (POPA). A combination of fluorescence quantification, zeta potential measurements, micropipette aspiration, and bilayer molecular dynamics simulations are used to characterize these membranes. The outer leaflet composition in vesicles is found consistent across the two preparation methods we employed, namely electroformation and inverted emulsion transfer. However, characterizing the inner leaflet poses challenges. Micropipette aspiration of GUVs show that oil residues do not substantially alter membrane elasticity, but simulations reveal increased membrane thickness and decreased interleaflet coupling in the presence of oil. Asymmetric membranes with a POPC:POPA mixture in the outer leaflet and POPC in the inner leaflet display similar stretching elasticity values to symmetric POPC:POPA membranes, suggesting potential POPA insertion into the inner leaflet during vesicle formation and suppressed asymmetry. The inverse compositional asymmetry, with POPC in the outer leaflet and POPC:POPA in the inner one yield less stretchable membranes with higher compressibility modulus compared with their symmetric counterparts. Challenges in achieving and predicting compositional correspondence highlight the limitations of phase-transfer-based methods. In addition, caution is advised when using fluorescently labeled lipids (even at low fractions of 0.5 mol %), as unexpected gel-like domains in symmetric POPC:POPA membranes were observed only with a specific type of labeled DOPE (dioleoylphosphatidylethanolamine) and the same fraction of unlabeled DOPE. The latter suggest that such domain formation may result from interactions between lipids and membrane fluorescent probes. Overall, this study underscores the complexity of factors influencing GUV membrane asymmetry, emphasizing the need for further research and improvement of characterization techniques.

Bioprinting of Synthetic Cell-like Lipid Vesicles to Augment the Functionality of Tissues after Manufacturing

Bioprinting is an automated bioassembly method that enables the formation of human tissue-like constructs to restore or replace damaged tissues. Regardless of the employed bioprinting method, cells undergo mechanical stress that can impact their survival and function postprinting. In this study, we investigate the use of a synthetic cell-like unit, giant unilamellar vesicles (GUVs), as adjuvants of the cellular function of human cells postprinting, or in future as the complete replacement of human cells. We analyzed the impact of two nozzle-based bioprinting methods (drop-on-demand and extrusion bioprinting) on the structure, stability, and function of GUVs. We showed that over 65% of the GUVs remain intact when printing at 0.5 bar, demonstrating the potential of using GUVs as a synthetic cell source. We further increased the stability of GUVs in a cell culture medium by introducing polyethylene glycol (PEG) into the GUV lipid membrane. The presence of PEG, however, diminished the structural properties of GUVs postprinting, and reduced the interaction of GUVs with human cells. Although the design of PEG-GUVs can still be modified in future studies for better cell-GUV interactions, we demonstrated that GUVs are functional postprinting. Chlorin e6-PEG-GUVs loaded with a fluorescent dye were bioprinted, and they released the dye postprinting only upon illumination. This is a new strategy to deliver carriers, such as growth factors, drugs, nutrients, or gases, inside large bioprinted specimens on a millimeter to centimeter scale. Overall, we showed that printed GUVs can augment the functionality of manufactured human tissues.

Multifaceted Activity of Fabimycin: Insights from Molecular Dynamics Studies on Bacterial Membrane Models

Membranes─cells' essential scaffolds─are valid molecular targets for substances with an antimicrobial effect. While certain substances, such as octenidine, have been developed to target membranes for antimicrobial purposes, the recently reported molecule, fabimycin (F2B)─a novel agent targeting drug-resistant Gram-negative bacteria─has not received adequate attention regarding its activity on membranes in the literature. The following study aims to investigate the effects of F2B on different bacterial membrane models, including simple planar bilayers and more complex bilayer systems that mimic the Escherichia coli shell equipped with double inner and outer bilayers. Our results show that F2B exhibited more pronounced interactions with bacterial membrane systems compared to the control PC system. Furthermore, we observed significant changes in local membrane property homeostasis in both the inner and outer membrane models, specifically in the case of lateral diffusion, membrane thickness, and membrane resilience (compressibility, tilt). Finally, our results showed that the effect of F2B differed in a complex system and a single membrane system. Our study provides new insights into the multifaceted activity of F2B, demonstrating its potential to disrupt bacterial membrane homeostasis, indicating that its activity extends the currently known mechanism of FabI enzyme inhibition. This disruption, coupled with the ability of F2B to penetrate the outer membrane layers, sheds new light on the behavior of this antimicrobial molecule. This highlights the importance of the interaction with the membrane, crucial in combating bacterial infections, particularly those caused by drug-resistant strains.

Light-Activated Synthetic Rotary Motors in Lipid Membranes Induce Shape Changes Through Membrane Expansion

Membranes are the key structures to separate and spatially organize cellular systems. Their rich dynamics and transformations during the cell cycle are orchestrated by specific membrane-targeted molecular machineries, many of which operate through energy dissipation. Likewise, man-made light-activated molecular rotary motors have previously shown drastic effects on cellular systems, but their physical roles on and within lipid membranes remain largely unexplored. Here, the impact of rotary motors on well-defined biological membranes is systematically investigated. Notably, dramatic mechanical transformations are observed in these systems upon motor irradiation, indicative of motor-induced membrane expansion. The influence of several factors on this phenomenon is systematically explored, such as motor concentration and membrane composition., Membrane fluidity is found to play a crucial role in motor-induced deformations, while only minor contributions from local heating and singlet oxygen generation are observed. Most remarkably, the membrane area expansion under the influence of the motors continues as long as irradiation is maintained, and the system stays out-of-equilibrium. Overall, this research contributes to a comprehensive understanding of molecular motors interacting with biological membranes, elucidating the multifaceted factors that govern membrane responses and shape transitions in the presence of these remarkable molecular machines, thereby supporting their future applications in chemical biology.

Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations

Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.

Integration of an LPAR1 Antagonist into Liposomes Enhances Their Internalization and Tumor Accumulation in an Animal Model of Human Metastatic Breast Cancer

Lysophosphatidic acid receptor 1 (LPAR1) is elevated in breast cancer. The deregulation of LPAR1, including the function and level of expression, is linked to cancer initiation, progression, and metastasis. LPAR1 antagonists, AM095 or Ki16425, may be effective therapeutic molecules, yet their limited water solubility hinders in vivo delivery. In this study, we report on the synthesis of two liposomal formulations incorporating AM095 or Ki16425, embedded within the lipid bilayer, as targeted nanocarriers for metastatic breast cancer (MBC). The data show that the Ki16425 liposomal formulation exhibited a 50% increase in internalization by MBC mouse epithelial cells (4T1) and a 100% increase in tumor accumulation in a mouse model of MBC compared with that of a blank liposomal formulation (control). At the same time, normal mouse epithelial cells (EpH-4Ev) internalized the Ki16425 liposomal formulation 25% lesser than the control formulation. Molecular dynamics simulations show that the integration of AM095 or Ki16425 modified the physical and mechanical properties of the lipid bilayer, making it more flexible in these liposomal formulations compared with liposomes without drug. The incorporation of an LPAR1 antagonist within a liposomal drug delivery system represents a viable therapeutic approach for targeting the LPA-LPAR1 axis, which may hinder the progression of MBC.

Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.

Effect of the Lipid Landscape on the Efficacy of Cell-Penetrating Peptides

Every cell biological textbook teaches us that the main role of the plasma membrane is to separate cells from their neighborhood to allow for a controlled composition of the intracellular space. The mostly hydrophobic nature of the cell membrane presents an impenetrable barrier for most hydrophilic molecules larger than 1 kDa. On the other hand, cell-penetrating peptides (CPPs) are capable of traversing this barrier without compromising membrane integrity, and they can do so on their own or coupled to cargos. Coupling biologically and medically relevant cargos to CPPs holds great promise of delivering membrane-impermeable drugs into cells. If the cargo is able to interact with certain cell types, uptake of the CPP-drug complex can be tailored to be cell-type-specific. Besides outlining the major membrane penetration pathways of CPPs, this review is aimed at deciphering how properties of the membrane influence the uptake mechanisms of CPPs. By summarizing an extensive body of experimental evidence, we argue that a more ordered, less flexible membrane structure, often present in the very diseases planned to be treated with CPPs, decreases their cellular uptake. These correlations are not only relevant for understanding the cellular biology of CPPs, but also for rationally improving their value in translational or clinical applications.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

The Effect of Cholesterol in SOPC Lipid Bilayers at Low Temperatures

We study the behavior of lipid bilayers composed of SOPC (1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine) with different concentrations of cholesterol, ranging from 10 mol% to 50 mol% at 273 K. To this end, we carry out extensive atomistic molecular dynamic simulations with the aid of the Slipid force field aiming at computing basic bilayer parameters, as well as thermodynamic properties and structural characteristics. The obtained results are compared to available relevant experimental data and the outcome of atomistic simulations performed on bilayers composed of analogous phospholipids. Our results show a good quantitative, as well as qualitative, agreement with the main trends associated with the concentration increase in cholesterol. Moreover, it comes out that a change in the behavior of the bilayer is brought about at a concentration of about 30 mol% cholesterol. At this very concentration, some of the bilayer properties are found to exhibit a saturation and a significant long-range ordering of the lipid molecules in the membrane shows up.

Biophysics of Membrane Stiffening by Cholesterol and Phosphatidylinositol 4,5-bisphosphate (PIP2)

Cell membranes regulate a wide range of phenomena that are implicated in key cellular functions. Cholesterol, a critical component of eukaryotic cell membranes, is responsible for cellular organization, membrane elasticity, and other critical physicochemical parameters. Besides cholesterol, other lipid components such as phosphatidylinositol 4,5-bisphosphate (PIP2) are found in minor concentrations in cell membranes yet can also play a major regulatory role in various cell functions. In this chapter, we describe how solid-state deuterium nuclear magnetic resonance (²H NMR) spectroscopy together with neutron spin-echo (NSE) spectroscopy can inform synergetic changes to lipid molecular packing due to cholesterol and PIP2 that modulate the bending rigidity of lipid membranes. Fundamental structure-property relations of molecular self-assembly are illuminated and point toward a length and time-scale dependence of cell membrane mechanics, with significant implications for biological activity and membrane lipid-protein interactions.

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.

Cholesterol Stiffening of Lipid Membranes

Biomembrane order, dynamics, and other essential physicochemical parameters are controlled by cholesterol, a major component of mammalian cell membranes. Although cholesterol is well known to exhibit a condensing effect on fluid lipid membranes, the extent of stiffening that occurs with different degrees of lipid acyl chain unsaturation remains an enigma. In this review, we show that cholesterol locally increases the bending rigidity of both unsaturated and saturated lipid membranes, suggesting there may be a length-scale dependence of the bending modulus. We review our published data that address the origin of the mechanical effects of cholesterol on unsaturated and polyunsaturated lipid membranes and their role in biomembrane functions. Through a combination of solid-state deuterium NMR spectroscopy and neutron spin-echo spectroscopy, we show that changes in molecular packing cause the universal effects of cholesterol on the membrane bending rigidity. Our findings have broad implications for the role of cholesterol in lipid–protein interactions as well as raft-like mixtures, drug delivery applications, and the effects of antimicrobial peptides on lipid membranes.

Partial Volumes of Phosphatidylcholines and Vitamin E: α-Tocopherol Prefers Disordered Membranes

Despite its discovery over 95 years ago, the biological and nutritional roles of vitamin E remain subjects of much controversy. Though it is known to possess antioxidant properties, recent assertions have implied that vitamin E may not be limited to this function in living systems. Through densitometry measurements and small-angle X-ray scattering we observe favorable interactions between α-tocopherol and unsaturated phospholipids, with more favorable interactions correlating to an increase in lipid chain unsaturation. Our data provide evidence that vitamin E may preferentially associate with oxygen sensitive lipids─an association that is considered innate for a viable membrane antioxidant.

Potentiation of plant defense by bacterial outer membrane vesicles is mediated by membrane nanodomains

Outer membrane vesicles (OMVs) are released from the outer membranes of Gram-negative bacteria during infection and modulate host immunity during host-pathogen interactions. The mechanisms by which OMVs are perceived by plants and affect host immunity are unclear. Here, we used the pathogen Xanthomonas campestris pv. campestris to demonstrate that OMV-plant interactions at the Arabidopsis thaliana plasma membrane (PM) modulate various host processes, including endocytosis, innate immune responses, and suppression of pathogenesis by phytobacteria. The lipid phase of OMVs is highly ordered and OMVs directly insert into the Arabidopsis PM, thereby enhancing the plant PM's lipid order; this also resulted in strengthened plant defenses. Strikingly, the integration of OMVs into the plant PM is host nanodomain- and remorin-dependent. Using coarse-grained simulations of molecular dynamics, we demonstrated that OMV integration into the plant PM depends on the membrane lipid order. Our computational simulations further showed that the saturation level of the OMV lipids could fine-tune the enhancement of host lipid order. Our work unraveled the mechanisms underlying the ability of OMVs produced by a plant pathogen to insert into the host PM, alter host membrane properties, and modulate plant immune responses.

Full scale structural, mechanical and dynamical properties of HIV-1 liposomes

Enveloped viruses are enclosed by a lipid membrane inside of which are all of the components necessary for the virus life cycle; viral proteins, the viral genome and metabolites. Viral envelopes are lipid bilayers that adopt morphologies ranging from spheres to tubes. The envelope is derived from the host cell during viral replication. Thus, the composition of the bilayer depends on the complex constitution of lipids from the host-cell's organelle(s) where assembly and/or budding of the viral particle occurs. Here, molecular dynamics (MD) simulations of authentic, asymmetric HIV-1 liposomes are used to derive a unique level of resolution of its full-scale structure, mechanics and dynamics. Analysis of the structural properties reveal the distribution of thicknesses of the bilayers over the entire liposome as well as its global fluctuations. Moreover, full-scale mechanical analyses are employed to derive the global bending rigidity of HIV-1 liposomes. Finally, dynamical properties of the lipid molecules reveal important relationships between their 3D diffusion, the location of lipid-rafts and the asymmetrical composition of the envelope. Overall, our simulations reveal complex relationships between the rich lipid composition of the HIV-1 liposome and its structural, mechanical and dynamical properties with critical consequences to different stages of HIV-1's life cycle.

Measuring the bending modulus of lipid bilayers with cholesterol

The effect of cholesterol on the bending modulus of DOPC lipid bilayers has become a controversial topic that has implications for methods of measuring the bending modulus. Recent results using neutron spin echo and nuclear magnetic resonance relaxation methods that involve linear transport properties have conflicted with earlier results from purely equilibrium experiments that do not involve linear transport properties. A general discussion indicates how one can be misled by data obtained by methods that involve linear transport properties. It is then shown specifically how the recent neutron spin echo results can be interpreted to agree with the earlier purely equilibrium experimental results, thereby resolving that conflict. Regarding the nuclear magnetic resonance relaxation method, it is noted that current interpretation of the data is unclear regarding the identity of the modulus that is involved, and an alternative interpretation is explored that does not disagree with the results of the equilibrium experiments.

Vectorial insertion of a β-helical peptide into membrane: a theoretical study on polytheonamide B

Spontaneous unidirectional, or vectorial, insertion of transmembrane peptides is a fundamental biophysical process for toxin and viral actions. Polytheonamide B (pTB) is a potent cytotoxic peptide with a β^6.3-helical structure. Previous experimental studies revealed that the pTB inserts into the membrane in a vectorial fashion and forms a channel with its single molecular length long enough to span the membrane. Also, molecular dynamics simulation studies demonstrated that the pTB is prefolded in aqueous solution. These are unique features of pTB because most of the peptide toxins form channels through oligomerization of transmembrane helices. Here, we performed all-atom molecular dynamics simulations to examine the dynamic mechanism of the vectorial insertion of pTB, providing underlying elementary processes of the membrane insertion of a prefolded single transmembrane peptide. We find that the insertion of pTB proceeds with only the local lateral compression of the membrane in three successive phases: "landing," "penetration," and "equilibration" phases. The free energy calculations using the replica-exchange umbrella sampling simulations present an energy cost of 4.3 kcal/mol at the membrane surface for the membrane insertion of pTB from bulk water. The trajectories of membrane insertion revealed that the insertion process can occur in two possible pathways, namely "trapped" and "untrapped" insertions; in some cases, pTB is trapped in the upper leaflet during the penetration phase. Our simulations demonstrated the importance of membrane anchoring by the hydrophobic N-terminal blocking group in the landing phase, leading to subsequent vectorial insertion.

Simple Does Not Mean Trivial: Behavior of Phosphatidic Acid in Lipid Mono- and Bilayers

Phosphatidic acid (PA) is one of the simplest membrane phospholipids, yet it plays a crucial role in various biologically relevant processes that take place in cells. Since PA generation may be triggered by a variety of factors, very often of antagonistic character, the specific nature of physiological responses driven by PA is not clear. In order to shed more light on these issues, we carried out a systematic characterization of membranes containing one of the three biologically significant PA molecular species. The effect of these molecules on the properties of membranes composed of phosphatidylcholine and/or cholesterol was assessed in a multidisciplinary approach, including molecular dynamic simulations, flicker noise spectroscopy, and Langmuir monolayer isotherms. The first enables the determination of various macroscopic and microscopic parameters such as lateral diffusion, membrane thickness, and defect analysis. The obtained data revealed a strong interaction between unsaturated PA species and phosphatidylcholine. On the other hand, the behavior of saturated PA was greatly influenced by cholesterol. Additionally, a strong effect on mechanical properties was observed in the case of three-component systems, which could not be explained by the simple extrapolation of parameters of the corresponding two-component systems. Our data show that various PA species are not equivalent in terms of their influence on lipid mono- and bilayers and that membrane composition/properties, particularly those related to the presence of cholesterol, may strongly modulate PA behavior.

A Systematic Approach: Molecular Dynamics Study and Parametrisation of Gemini Type Cationic Surfactants

The spreading of antibiotic-resistant bacteria strains is one of the most serious problem in medicine to struggle nowadays. This triggered the development of alternative antimicrobial agents in recent years. One of such group is Gemini surfactants which are massively synthesised in various structural configurations to obtain the most effective antibacterial properties. Unfortunately, the comparison of antimicrobial effectiveness among different types of Gemini agents is unfeasible since various protocols for the determination of Minimum Inhibitory Concentration are used. In this work, we proposed alternative, computational, approach for such comparison. We designed a comprehensive database of 250 Gemini surfactants. Description of structure parameters, for instance spacer type and length, are included in the database. We parametrised modelled molecules to obtain force fields for the entire Gemini database. This was used to conduct in silico studies using the molecular dynamics to investigate the incorporation of these agents into model E. coli inner membrane system. We evaluated the effect of Gemini surfactants on structural, stress and mechanical parameters of the membrane after the agent incorporation. This enabled us to select four most likely membrane properties that could correspond to Gemini’s antimicrobial effect. Based on our results we selected several types of Gemini spacers which could demonstrate a particularly strong effect on the bacterial membranes.

Cholesterol Hinders the Passive Uptake of Amphiphilic Nanoparticles into Fluid Lipid Membranes

Plasma membranes represent pharmacokinetic barriers for the passive transport of site-specific drugs within cells. When engineered nanoparticles (NPs) are considered as transmembrane drug carriers, the plasma membrane composition can affect passive NP internalization in many ways. Among these, cholesterol-regulated membrane fluidity is probably one of the most biologically relevant. Herein, we consider small (2-5 nm in core diameter) amphiphilic gold NPs capable of spontaneously and nondisruptively entering the lipid bilayer of plasma membranes. We study their incorporation into model 1,2-dioleoyl-sn-glycero-3-phosphocholine membranes with increasing cholesterol content. We combine dissipative quartz crystal microbalance experiments, atomic force microscopy, and molecular dynamics simulations to show that membrane cholesterol, at biologically relevant concentrations, hinders the molecular mechanism for passive NP penetration within fluid bilayers, resulting in a dramatic reduction in the amount of NP incorporated.

The dynamic face of lipid membranes

Cell membranes - primarily composed of lipids, sterols, and proteins - form a dynamic interface between living cells and their environment. They act as a mechanical barrier around the cell while selectively facilitating material transport, signal transduction, and various other functions necessary for the cell viability. The complex functionality of cell membranes and the hierarchical motions and responses they exhibit demand a thorough understanding of the origin of different membrane dynamics and how they are influenced by molecular additives and environmental cues. These dynamic modes include single-molecule diffusion, thermal fluctuations, and large-scale membrane deformations, to name a few. This review highlights advances in investigating structure-driven dynamics associated with model cell membranes, with a particular focus on insights gained from neutron scattering and spectroscopy experiments. We discuss the uniqueness of neutron contrast variation and its remarkable potential in probing selective membrane structure and dynamics on spatial and temporal scales over which key biological functions occur. We also present a summary of current and future opportunities in synergistic combinations of neutron scattering with molecular dynamics (MD) simulations to gain further understanding of the molecular mechanisms underlying complex membrane functions.

Phenothiazines alter plasma membrane properties and sensitize cancer cells to injury by inhibiting annexin-mediated repair

Repair of damaged plasma membrane in eukaryotic cells is largely dependent on the binding of annexin repair proteins to phospholipids. Changing the biophysical properties of the plasma membrane may provide means to compromise annexin-mediated repair and sensitize cells to injury. Since, cancer cells experience heightened membrane stress and are more dependent on efficient plasma membrane repair, inhibiting repair may provide approaches to sensitize cancer cells to plasma membrane damage and cell death. Here, we show that derivatives of phenothiazines, which have widespread use in the fields of psychiatry and allergy treatment, strongly sensitize cancer cells to mechanical-, chemical-, and heat-induced injury by inhibiting annexin-mediated plasma membrane repair. Using a combination of cell biology, biophysics, and computer simulations, we show that trifluoperazine acts by thinning the membrane bilayer, making it more fragile and prone to ruptures. Secondly, it decreases annexin binding by compromising the lateral diffusion of phosphatidylserine, inhibiting the ability of annexins to curve and shape membranes, which is essential for their function in plasma membrane repair. Our results reveal a novel avenue to target cancer cells by compromising plasma membrane repair in combination with noninvasive approaches that induce membrane injuries.

Mechanical properties of anionic asymmetric bilayers from atomistic simulations

Mechanotransduction, the biological response to mechanical stress, is often initiated by activation of mechanosensitive (MS) proteins upon mechanically induced deformations of the cell membrane. A current challenge in fully understanding this process is in predicting how lipid bilayers deform upon the application of mechanical stress. In this context, it is now well established that anionic lipids influence the function of many proteins. Here, we test the hypothesis that anionic lipids could indirectly modulate MS proteins by alteration of the lipid bilayer mechanical properties. Using all-atom molecular dynamics simulations, we computed the bilayer bending rigidity (KC), the area compressibility (KA), and the surface shear viscosity (ηm) of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (PC) lipid bilayers with and without phosphatidylserine (PS) or phosphatidylinositol bisphosphate (PIP2) at physiological concentrations in the lower leaflet. Tensionless leaflets were first checked for each asymmetric bilayer model, and a formula for embedding an asymmetric channel in an asymmetric bilayer is proposed. Results from two different sized bilayers show consistently that the addition of 20% surface charge in the lower leaflet of the PC bilayer with PIP2 has minimal impact on its mechanical properties, while PS reduced the bilayer bending rigidity by 22%. As a comparison, supplementing the PIP2-enriched PC membrane with 30% cholesterol, a known rigidifying steroid lipid, produces a significant increase in all three mechanical constants. Analysis of pairwise splay moduli suggests that the effect of anionic lipids on bilayer bending rigidity largely depends on the number of anionic lipid pairs formed during simulations. The potential implication of bilayer bending rigidity is discussed in the framework of MS piezo channels.

Biomembrane Structure and Material Properties Studied With Neutron Scattering

Cell membranes and their associated structures are dynamical supramolecular structures where different physiological processes take place. Detailed knowledge of their static and dynamic structures is therefore needed, to better understand membrane biology. The structure–function relationship is a basic tenet in biology and has been pursued using a range of different experimental approaches. In this review, we will discuss one approach, namely the use of neutron scattering techniques as applied, primarily, to model membrane systems composed of lipid bilayers. An advantage of neutron scattering, compared to other scattering techniques, is the differential sensitivity of neutrons to isotopes of hydrogen and, as a result, the relative ease of altering sample contrast by substituting protium for deuterium. This property makes neutrons an ideal probe for the study of hydrogen-rich materials, such as biomembranes. In this review article, we describe isotopic labeling studies of model and viable membranes, and discuss novel applications of neutron contrast variation in order to gain unique insights into the structure, dynamics, and molecular interactions of biological membranes. We specifically focus on how small-angle neutron scattering data is modeled using different contrast data and molecular dynamics simulations. We also briefly discuss neutron reflectometry and present a few recent advances that have taken place in neutron spin echo spectroscopy studies and the unique membrane mechanical data that can be derived from them, primarily due to new models used to fit the data.

Curvature-Regulated Lipid Membrane Softening of Nano-Vesicles

The physico-mechanical properties of nanoscale lipid vesicles (e.g., natural nano-vesicles and artificial nano-liposomes) dictate their interaction with biological systems. Understanding the interplay between vesicle size and stiffness is critical to both the understanding of the biological functions of natural nano-vesicles and the optimization of nano-vesicle-based diagnostics and therapeutics. It has been predicted that, when vesicle size is comparable to its membrane thickness, the effective bending stiffness of the vesicle increases dramatically due to both the entropic effect as a result of reduced thermal undulation and the nonlinear curvature elasticity effect. Through systematic molecular dynamics simulations, we show that the vesicle membrane thins and softens with the decrease in vesicle size, which effectively counteracts the stiffening effects as already mentioned. Our simulations indicate that the softening of nano-vesicles results from a change in the bilayer's interior structure - a decrease in lipid packing order - as the membrane curvature increases. Our work thus leads to a more complete physical framework to understand the physico-mechanical properties of nanoscale lipid vesicles, paving the way to further advances in the biophysics of nano-vesicles and their biomedical applications.

Effect of quercetin on lipid membrane rigidity: assessment by atomic force microscopy and molecular dynamics simulations

Quercetin (3,3',4',5,7-pentahydroxyl-flavone) is a natural flavonoid with many valuable biological effects, but its solubility in water is low, posing major limitations in applications. Quercetin encapsulation in liposomes increases its bioavailability; the drug effect on liposome elastic properties is required for formulation development. Here, we quantify the effect of quercetin molecules on the rigidity of lipoid E80 liposomes using atomic force microscopy (AFM) and molecular dynamics (MD) simulations. AFM images show no effect of quercetin molecules on liposomes morphology and structure. However, AFM force curves suggest that quercetin softens lipid membranes; the Young modulus measured for liposomes encapsulating quercetin is smaller than that determined for blank liposomes. We then used MD simulations to interpret the effect of quercetin on membrane rigidity in terms of molecular interactions. The decrease in membrane rigidity was confirmed by the simulations, which also revealed that quercetin affects structural and dynamic properties: membrane thickness is decreased, acyl chains disorder is increased, and diffusion coefficients of lipid molecules are also increased. Such changes appear to be related to the preferential localization of quercetin within the membrane, near the interface between the hydrophobic core and polar head groups of the lipids.

Guardians of the Cell: State-of-the-Art of Membrane Proteins from a Computational Point-of-View

Membrane proteins (MPs) encompass a large family of proteins with distinct cellular functions, and although representing over 50% of existing pharmaceutical drug targets, their structural and functional information is still very scarce. Over the last years, in silico analysis and algorithm development were essential to characterize MPs and overcome some limitations of experimental approaches. The optimization and improvement of these methods remain an ongoing process, with key advances in MPs' structure, folding, and interface prediction being continuously tackled. Herein, we discuss the latest trends in computational methods toward a deeper understanding of the atomistic and mechanistic details of MPs.

Site-Specific Peroxidation Modulates Lipid Bilayer Mechanics

Peroxidation of plasma membranes, characterized by oxidative attack of lipidic carbon-carbon double-bonds in unsaturated fatty acids, has been identified as an important biochemical event in multiple pathological conditions, including neurodegenerative diseases, atherosclerosis, diabetes, preeclampsia, aging, cancer, etc. Changes to the lipid bilayer structure as a result of lipid peroxidation may lead to lipid membrane malfunction, and consequently initiate further downstream biochemical cascades. However, how lipid peroxidation modulates the mechanical properties of lipid membranes remains largely controversial. In this study, we investigate the peroxidation of lipids with polyunsaturated fatty acid tails using molecular dynamics simulations. By systematically varying the oxidation site, we find that lipid peroxidation alters the biophysical properties of bilayer membrane in a peroxidation site-specific manner. Specifically, our results suggest that peroxidation at sites in the bilayer interior disturbs and softens the membrane, whereas peroxidation at sites near the membrane-water interface results in a more ordered and stiffer membrane. Such a peroxidation site-specific modulation of lipid membrane mechanics provides an explanation for the contradictory results obtained in previous experiments. Our study paves the way for an improved understanding of the initiation of the downstream cellular dysfunction caused by lipid peroxidation.

Stabilizing Leaflet Asymmetry under Differential Stress in a Highly Coarse-Grained Lipid Membrane Model

We present a version of the coarse-grained Cooke lipid model, modified to simulate asymmetric lipid membranes. It is inspired by a method employed by Wang et al. [ Commun. Comput. Phys. 2013, 13, 1093-1106] for artificially penalizing lipid flip-flop but copes more robustly with differential stress, at the cost of one additional bead per lipid and the concomitant increase in computational overhead. Bilayer asymmetry ultimately breaks down beyond a system size dependent critical differential stress, which can be predicted from a simple analytical model. We remeasure many important material parameters for the new model and find it to be consistent with typical fluid lipid membranes. Maintaining a stable stress asymmetry has many applications, and we give two examples: (i) connecting monolayer stress to lipid number asymmetry in order to directly measure the monolayer area modulus and (ii) finding its strain-dependent higher-order correction by monitoring the equilibrium bilayer area.

Scaling relationships for the elastic moduli and viscosity of mixed lipid membranes

The elastic and viscous properties of biological membranes play a vital role in controlling cell functions that require local reorganization of the membrane components as well as dramatic shape changes such as endocytosis, vesicular trafficking, and cell division. These properties are widely acknowledged to depend on the unique composition of lipids within the membrane, yet the effects of lipid mixing on the membrane biophysical properties remain poorly understood. Here, we present a comprehensive characterization of the structural, elastic, and viscous properties of fluid membranes composed of binary mixtures of lipids with different tail lengths. We show that the mixed lipid membrane properties are not simply additive quantities of the single-component analogs. Instead, the mixed membranes are more dynamic than either of their constituents, quantified as a decrease in their bending modulus, area compressibility modulus, and viscosity. While the enhanced dynamics are seemingly unexpected, we show that the measured moduli and viscosity for both the mixed and single-component bilayers all scale with the area per lipid and collapse onto respective master curves. This scaling links the increase in dynamics to mixing-induced changes in the lipid packing and membrane structure. More importantly, the results show that the membrane properties can be manipulated through lipid composition the same way bimodal blends of surfactants, liquid crystals, and polymers are used to engineer the mechanical properties of soft materials, with broad implications for understanding how lipid diversity relates to biomembrane function.

How cholesterol stiffens unsaturated lipid membranes

Cholesterol is an integral component of eukaryotic cell membranes and a key molecule in controlling membrane fluidity, organization, and other physicochemical parameters. It also plays a regulatory function in antibiotic drug resistance and the immune response of cells against viruses, by stabilizing the membrane against structural damage. While it is well understood that, structurally, cholesterol exhibits a densification effect on fluid lipid membranes, its effects on membrane bending rigidity are assumed to be nonuniversal; i.e., cholesterol stiffens saturated lipid membranes, but has no stiffening effect on membranes populated by unsaturated lipids, such as 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). This observation presents a clear challenge to structure-property relationships and to our understanding of cholesterol-mediated biological functions. Here, using a comprehensive approach-combining neutron spin-echo (NSE) spectroscopy, solid-state deuterium NMR (²H NMR) spectroscopy, and molecular dynamics (MD) simulations-we report that cholesterol locally increases the bending rigidity of DOPC membranes, similar to saturated membranes, by increasing the bilayer's packing density. All three techniques, inherently sensitive to mesoscale bending fluctuations, show up to a threefold increase in effective bending rigidity with increasing cholesterol content approaching a mole fraction of 50%. Our observations are in good agreement with the known effects of cholesterol on the area-compressibility modulus and membrane structure, reaffirming membrane structure-property relationships. The current findings point to a scale-dependent manifestation of membrane properties, highlighting the need to reassess cholesterol's role in controlling membrane bending rigidity over mesoscopic length and time scales of important biological functions, such as viral budding and lipid-protein interactions.

Direct label-free imaging of nanodomains in biomimetic and biological membranes by cryogenic electron microscopy

The nanoscale organization of biological membranes into structurally and compositionally distinct lateral domains is believed to be central to membrane function. The nature of this organization has remained elusive due to a lack of methods to directly probe nanoscopic membrane features. We show here that cryogenic electron microscopy (cryo-EM) can be used to directly image coexisting nanoscopic domains in synthetic and bioderived membranes without extrinsic probes. Analyzing a series of single-component liposomes composed of synthetic lipids of varying chain lengths, we demonstrate that cryo-EM can distinguish bilayer thickness differences as small as 0.5 Å, comparable to the resolution of small-angle scattering methods. Simulated images from computational models reveal that features in cryo-EM images result from a complex interplay between the atomic distribution normal to the plane of the bilayer and imaging parameters. Simulations of phase-separated bilayers were used to predict two sources of contrast between coexisting ordered and disordered phases within a single liposome, namely differences in membrane thickness and molecular density. We observe both sources of contrast in biomimetic membranes composed of saturated lipids, unsaturated lipids, and cholesterol. When extended to isolated mammalian plasma membranes, cryo-EM reveals similar nanoscale lateral heterogeneities. The methods reported here for direct, probe-free imaging of nanodomains in unperturbed membranes open new avenues for investigation of nanoscopic membrane organization.

Membrane lipids are both the substrates and a mechanistically responsive environment of TMEM16 scramblase proteins

Recent discoveries about functional mechanisms of proteins in the TMEM16 family of phospholipid scramblases have illuminated the dual role of the membrane as both the substrate and a mechanistically responsive environment in the wide range of physiological processes and genetic disorders in which they are implicated. This is highlighted in the review of recent findings from our collaborative investigations of molecular mechanisms of TMEM16 scramblases that emerged from iterative functional, structural, and computational experimentation. In the context of this review, we present new MD simulations and trajectory analyses motivated by the fact that new structural information about the TMEM16 scramblases is emerging from cryo-EM determinations in lipid nanodiscs. Because the functional environment of these proteins in in vivo and in in vitro is closer to flat membranes, we studied comparatively the responses of the membrane to the TMEM16 proteins in flat membranes and nanodiscs. We find that bilayer shapes in the nanodiscs are very different from those observed in the flat membrane systems, but the function-related slanting of the membrane observed at the nhTMEM16 boundary with the protein is similar in the nanodiscs and in the flat bilayers. This changes, however, in the bilayer composed of longer-tail lipids, which is thicker near the phospholipid translocation pathway, which may reflect an enhanced tendency of the long tails to penetrate the pathway and create, as shown previously, a nonconductive environment. These findings support the correspondence between the mechanistic involvement of the lipid environment in the flat membranes, and the nanodiscs. © 2019 Wiley Periodicals, Inc.

Regulation of lipid saturation without sensing membrane fluidity

Cells maintain membrane fluidity by regulating lipid saturation, but the molecular mechanisms of this homeoviscous adaptation remain poorly understood. We have reconstituted the core machinery for regulating lipid saturation in baker’s yeast to study its molecular mechanism. By combining molecular dynamics simulations with experiments, we uncover a remarkable sensitivity of the transcriptional regulator Mga2 to the abundance, position, and configuration of double bonds in lipid acyl chains, and provide insights into the molecular rules of membrane adaptation. Our data challenge the prevailing hypothesis that membrane fluidity serves as the measured variable for regulating lipid saturation. Rather, we show that Mga2 senses the molecular lipid-packing density in a defined region of the membrane. Our findings suggest that membrane property sensors have evolved remarkable sensitivities to highly specific aspects of membrane structure and dynamics, thus paving the way toward the development of genetically encoded reporters for such properties in the future.

Cells maintain membrane fluidity by regulating lipid saturation, but the molecular mechanisms of this homeoviscous adaptation remain poorly understood. Here authors reconstituted the core machinery for regulating lipid saturation in baker’s yeast to directly characterize its response to defined membrane environments and uncover its mode-of-action.

Area Compressibility Moduli of the Monolayer Leaflets of Asymmetric Bilayers from Simulations

Extraction from simulations of the area compressibility moduli of the monolayers in a bilayer is considered theoretically. A statistical mechanical derivation shows that the bilayer modulus is the sum of the two monolayer moduli, as is often supposed but contrary to a recent study. Seemingly plausible assumptions regarding fluctuations are tested rigorously. Prospects for future research are discussed.

Modulation of Antimicrobial Peptide Potency in Stressed Lipid Bilayers

It is shown that the tendency of an archetypal antimicrobial peptide to insert into and perforate a simple lipid bilayer is strongly modulated by tensile stress in the membrane. The results, obtained through molecular dynamics simulations, have been demonstrated with several lipid compositions and appear to be general, although quantitative details differ. The findings imply that the potency of antimicrobial peptides may not be a purely intrinsic chemical property and, instead, depends on the mechanical state of the target membrane.