Structural and mechanical properties of cardiolipin lipid bilayers determined using neutron spin echo, small angle neutron and X-ray scattering, and molecular dynamics simulations

Binding of yeast and human cytochrome c to cardiolipin nanodiscs at physiological ionic strength

Binding of cytochrome c (Cytc) to membranes containing cardiolipin (CL) is of considerable interest because of the importance of this interaction in the early stages of apoptosis. The molecular-level determinants of this interaction are still not well defined and there appear to be species-specific differences in Cytc affinity for CL-containing membranes. Many studies are carried out at low ionic strength far from the 100-150 mM ionic strength within mitochondria. Similarly, most binding studies are done at Cytc concentrations of 10 μM or less, much lower that the estimated range of 0.1 to 5 mM Cytc present in mitochondria. In this study, we evaluate binding of human and yeast Cytc to CL nanodiscs using size exclusion chromatography at 25 μM Cytc concentration and 100 mM ionic strength. We find that yeast Cytc affinity for CL nanodiscs is much stronger than that of human Cytc. Mutational analysis of the site A binding surface shows that lysines 86 and 87 are more important for yeast Cytc binding to CL nanodiscs than lysines 72 and 73, counter to results at lower ionic strength. Analysis of the electrostatic surface potential of human versus yeast Cytc shows that the positive potential due to lysines 86 and 87 and other nearby lysines (4, 5, 11, 89) is stronger than that due to lysines 72 and 73. In the case of human Cytc the positive potential around site A is less uniform and likely weakens electrostatic binding to CL membranes through site A.

Effect of cardiolipin on the lamellarity and elongation of liposomes hydrated in PBS

Lamellarity and shape are important factors in the formation of vesicles and determine their role in biological systems and pharmaceutical applications. Cardiolipin (CL) is a major lipid in many biological membranes and exerts a great influence on their structural organization due to its particular structure and physico-chemical properties. Here, we used small-angle X-ray and neutron scattering to study the effects of CL with different acyl chain lengths and saturations (CL14:0, CL18:1, CL18:2) on vesicle morphology and lamellarity in membrane models containing mixtures of phosphatidylcholine and phosphatidylethanolamine with different acyl chain lengths and saturations (C14:0 and C 18:1). Measurements were performed in the presence of Phosphate Buffer Saline (PBS), at 37°C, to better reflect physiological conditions, which resulted in strong effects on vesicle morphology, depending on the type and amount of CL used. The presence of small quantities of CL (from 2.5%) reduced inter-membrane correlations and increased perturbation of the membrane, an effect which is enhanced in the presence of matched shorter saturated acyl chains, and mainly unilamellar vesicles (ULV) are formed. In extruded vesicles, employed for SANS experiments, flattened vesicles are observed partly due to the hypertonic effect of PBS, but also influenced by the type of CL added. Our experimental data from SAXS and SANS revealed a strong dependence on CL content in shaping the membrane microstructure, with an apparent optimum in the PC:CL mixture in terms of promoting reduced correlations, preferred curvature and elongation. However, the use of PBS caused distinct differences from previously published studies in water in terms of vesicle shape, and highlights the need to investigate vesicle formation under physiological conditions in order to be able to draw conclusions about membrane formation in biological systems.

Small-angle X-ray and neutron scattering applied to lipid-based nanoparticles: Recent advancements across different length scales

Lipid-based nanoparticles (LNPs), ranging from nanovesicles to non-lamellar assemblies, have gained significant attention in recent years, as versatile carriers for delivering drugs, vaccines, and nutrients. Small-angle scattering methods, employing X-rays (SAXS) or neutrons (SANS), represent unique tools to unveil structure, dynamics, and interactions of such particles on different length scales, spanning from the nano to the molecular scale. This review explores the state-of-the-art on scattering methods applied to unveil the structure of lipid-based nanoparticles and their interactions with drugs and bioactive molecules, to inform their rational design and formulation for medical applications. We will focus on complementary information accessible with X-rays or neutrons, ranging from insights on the structure and colloidal processes at a nanoscale level (SAXS) to details on the lipid organization and molecular interactions of LNPs (SANS). In addition, we will review new opportunities offered by Time-resolved (TR)-SAXS and -SANS for the investigation of dynamic processes involving LNPs. These span from real-time monitoring of LNPs structural evolution in response to endogenous or external stimuli (TR-SANS), to the investigation of the kinetics of lipid diffusion and exchange upon interaction with biomolecules (TR-SANS). Finally, we will spotlight novel combinations of SAXS and SANS with complementary on-line techniques, recently enabled at Large Scale Facilities for X-rays and neutrons. This emerging technology enables synchronized multi-method investigation, offering exciting opportunities for the simultaneous characterization of the structure and chemical or mechanical properties of LNPs.

Curvature sensing lipid dynamics in a mitochondrial inner membrane model

Membrane curvature is essential for many cellular structures and processes, and factors such as leaflet asymmetry, lipid composition, and proteins all play important roles. Cardiolipin is the signature lipid of mitochondrial membranes and is essential for maintaining the highly curved shapes of the inner mitochondrial membrane (IMM) and the spatial arrangement of membrane proteins. In this study, we investigate the partitioning behavior of various lipids present in the IMM using coarse-grained molecular dynamics simulations. This study explores curved bilayer systems containing phosphatidylcholine (PC), phosphatidylethanolamine (PE), and cardiolipin (CDL) in binary and ternary component mixtures. Curvature properties such as mean and Gaussian curvatures, as well as the distribution of lipids into the various curved regions of the cristae models, are quantified. Overall, this work represents an advance beyond previous studies on lipid curvature sensing by simulating these systems in a geometry that has the morphological features and scales of curvature consistent with regions of the IMM. We find that CDL has a stronger preference for accumulating in regions of negative curvature than PE lipids, in agreement with previous results. Furthermore, we find lipid partitioning propensity is dominated by sensitivity to mean curvature, while there is a weaker correlation with Gaussian curvature.

Structure and Elasticity of Mitochondrial Membranes: A Molecular Dynamics Simulation Study

Mitochondria are known as the powerhouse of the cell because they produce energy in the form of adenosine triphosphate. They also have other crucial functions such as regulating apoptosis, calcium homeostasis, and reactive oxygen species production. To perform these diverse functions, mitochondria adopt specific structures and frequently undergo dynamic shape changes, indicating that their mechanical properties play an essential role in their functions. To gain a detailed understanding at the molecular level of the structure and mechanical properties of mitochondria, we carry out atomistic molecular dynamics simulations for three inner mitochondrial membranes and three outer mitochondrial membrane models. These models take into account variations in cardiolipin and cholesterol concentrations as well as the symmetry/asymmetry between the two leaflets. Our simulations allow us to calculate various structural quantities and the bending, twisting, and tilting elastic moduli of the membrane models. Our results indicate that the structures of the inner and outer mitochondrial membranes are quite similar and do not depend much on the variation in lipid compositions. However, the bending modulus of the membranes increases with increasing concentrations of cardiolipin or cholesterol but decreases with a membrane asymmetry. Notably, we found that the dipole potential of the membrane increases with an increasing cardiolipin concentration. Finally, possible roles of cardiolipin in regulating ion and proton currents and maintaining the cristate are discussed in some details.

A review on the measurement of the bending rigidity of lipid membranes

This review provides an overview of the latest developments in both experimental and simulation techniques used to assess the bending rigidity of lipid membranes. It places special emphasis on experimental methods that utilize model vesicles to manipulate lipid compositions and other experimental parameters to determine the bending rigidity of the membrane. It also describes two commonly used simulation methods for estimating bending rigidity. The impact of various factors on membrane bending rigidity is summarized, including cholesterol, lipids, salt concentration, surface charge, membrane phase state, peptides, proteins, and polyethylene glycol. These factors are shown to influence the bending rigidity, contributing to a better understanding of the biophysical properties of membranes and their role in biological processes. Furthermore, the review discusses future directions and potential advancements in this research field, highlighting areas where further investigation is required.

Mitochondrial membrane model: Lipids, elastic properties, and the changing curvature of cardiolipin

Mammalian and Drosophila melanogaster model mitochondrial membrane compositions are constructed from experimental data. Simplified compositions for inner and outer mitochondrial membranes are provided, including an asymmetric inner mitochondrial membrane. We performed atomistic molecular dynamics simulations of these membranes and computed their material properties. When comparing these properties to those obtained by extrapolation from their constituting lipids, we find good overall agreement. Finally, we analyzed the curvature effect of cardiolipin, considering ion concentration effects, oxidation, and pH. We draw the conclusion that cardiolipin-negative curvature is most likely due to counterion effects, such as cation adsorption, in particular of H3O+. This oft-neglected effect might account for the puzzling behavior of this lipid.

The influence of Leucidal - eco-preservative from radish - on model lipid membranes and selected pathogenic bacteria

In this work the effect of Leucidal - a natural preservative from radish dedicated to be used in cosmetics - on bacteria cells and model bacteria membranes was investigated. To get insight into the mechanism of action of this formulation the lipid Langmuir monolayers imitating Escherichia coli (E. coli) and Staphylococcus aureus (S. aureus) membranes were prepared. Then, the influence of Leucidal on model systems was investigated by means of the surface pressure/area measurements, penetration studies and Brewster Angle Microscopy (BAM) visualization. Similar experiments were done also for one component monolayers formed from the model membrane lipids. The in vitro tests were done on five different bacteria species (E. coli, Enterococcus faecalis, S. aureus, Salmonella enterica, Pseudomonas aeruginosa). Leucidal was found to decrease packing of the monolayers, however, it was excluded from the films at higher concentrations. Model membrane experiments evidenced also a stronger affinity of the components of this eco-preservative to E. coli vs S. aureus membrane. Among one component films, those formed from phosphatidylglycerols and cardiolipins were more sensitive to the presence of Leucidal. However, in vitro tests evidenced that Leucidal exerts stronger inhibitory effect against S. aureus bacteria as compared to E. coli strain. These findings were discussed from the point of view of the role of Leucidal components and the lipid membrane properties in the membrane - based mechanism of action of this preservative. The results allow one to suggest that the membrane may not be the main site of action of Leucidal on bacteria. Moreover, since high concentration of the tested preparation exerted antibacterial activity in relation to all tested bacteria, a low selectivity of Leucidal can be postulated, which may be problematic from the point of view of its effect on the skin microbiome.

Small-angle scattering from flat bilayers containing correlated scattering length density inhomogeneities

Model lipid bilayers have been widely employed as a minimal system to investigate the structural properties of biological membranes by small-angle X-ray (SAXS) and neutron scattering (SANS) techniques. These have nanometre resolution and can give information regarding membrane thickness and scattering length densities (SLDs) of polar and apolar regions. However, biological membranes are complex systems containing different lipids and protein species, in which lipid domains can be dynamically assembled and disassembled. Therefore, SLD variations can occur within the biomembrane. In this work, a novel method has been developed to simulate SAXS and SANS profiles obtained from large unilamellar vesicles containing SLD inhomogeneities that are spatially correlated over the membrane surface. Such inhomogeneities are represented by cylindrical entities with equivalent SLDs. Stacking of bilayers is also included in the model, with no correlation between horizontal and vertical order. The model is applied to a lipid bilayer containing SLD inhomogeneities representing pores, lipid domains, and transmembrane, partially immersed and anchored proteins. It is demonstrated that all the structural information from the host lipid bilayer and from the SLD inhomogeneity can be consistently retrieved by a combined analysis of experimental SAXS and SANS data through the methodology proposed here.

Neutron scattering studies on dynamics of lipid membranes

Neutron scattering methods are powerful tools for the study of the structure and dynamics of lipid bilayers in length scales from sub Å to tens to hundreds nm and the time scales from sub ps to μs. These techniques also are nondestructive and, perhaps most importantly, require no additives to label samples. Because the neutron scattering intensities are very different for hydrogen- and deuterium-containing molecules, one can replace the hydrogen atoms in a molecule with deuterium to prepare on demand neutron scattering contrast without significantly altering the physical properties of the samples. Moreover, recent advances in neutron scattering techniques, membrane dynamics theories, analysis tools, and sample preparation technologies allow researchers to study various aspects of lipid bilayer dynamics. In this review, we focus on the dynamics of individual lipids and collective membrane dynamics as well as the dynamics of hydration water.

The Transport of Charged Molecules across Three Lipid Membranes Investigated with Second Harmonic Generation

Subtle variations in the structure and composition of lipid membranes can have a profound impact on their transport of functional molecules and relevant cell functions. Here, we present a comparison of the permeability of bilayers composed of three lipids: cardiolipin, DOPG (1,2-dioleoyl-sn-glycero-3-phospho-(1'-rac-glycerol), and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol)). The adsorption and cross-membrane transport of a charged molecule, D289 (4-(4-diethylaminostyry)-1-methyl-pyridinium iodide), on vesicles composed of the three lipids were monitored by second harmonic generation (SHG) scattering from the vesicle surface. It is revealed that structural mismatching between the saturated and unsaturated alkane chains in POPG leads to relatively loose packing structure in the lipid bilayers, thus providing better permeability compared to unsaturated lipid bilayers (DOPG). This mismatching also weakens the efficiency of cholesterol in rigidifying the lipid bilayers. It is also revealed that the bilayer structure is somewhat disturbed by the surface curvature in small unilamellar vesicles (SUVs) composed of POPG and the conical structured cardiolipin. Such subtle information on the relationship between the lipid structure and the molecular transport capability of the bilayers may provide clues for drug development and other medical and biological studies.

Small-Angle Neutron Scattering Reveals the Nanostructure of Liposomes with Embedded OprF Porins of Pseudomonas aeruginosa

The use of liposomes as drug delivery systems emerged in the last decades in view of their capacity and versatility to deliver a variety of therapeutic agents. By means of small-angle neutron scattering (SANS), we performed a detailed characterization of liposomes containing outer membrane protein F (OprF), the main porin of the Pseudomonas aeruginosa bacterium outer membrane. These OprF-liposomes are the basis of a novel vaccine against this antibiotic-resistant bacterium, which is one of the main hospital-acquired pathogens and causes each year a significant number of deaths. SANS data were analyzed by a specific model we created to quantify the crucial information about the structure of the liposome containing OprF, including the lipid bilayer structure, the amount of protein in the lipid bilayer, the average protein localization, and the effect of the protein incorporation on the lipid bilayer. Quantification of such structural information is important to enhance the design of liposomal delivery systems for therapeutic applications.

Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes

Small-angle neutron scattering (SANS) is a powerful tool for studying biological membranes and model lipid bilayer membranes. The length scales probed by SANS, being from 1 nm to over 100 nm, are well-matched to the relevant length scales of the bilayer, particularly when it is in the form of a vesicle. However, it is the ability of SANS to differentiate between isotopes of hydrogen as well as the availability of deuterium labeled lipids that truly enable SANS to reveal details of membranes that are not accessible with the use of other techniques, such as small-angle X-ray scattering. In this work, an overview of the use of SANS for studying unilamellar lipid bilayer vesicles is presented. The technique is briefly presented, and the power of selective deuteration and contrast variation methods is discussed. Approaches to modeling SANS data from unilamellar lipid bilayer vesicles are presented. Finally, recent examples are discussed. While the emphasis is on studies of unilamellar vesicles, examples of the use of SANS to study intact cells are also presented.

Mitochondrial membrane models built from native lipid extracts: Interfacial and transport properties

The mitochondrion is an essential organelle enclosed by two membranes whose functionalities depend on their very specific protein and lipid compositions. Proteins from the outer mitochondrial membrane (OMM) are specialized in mitochondrial dynamics and mitophagy, whereas proteins of the inner mitochondrial membrane (IMM) have dedicated functions in cellular respiration and apoptosis. As for lipids, the OMM is enriched in glycerophosphatidyl choline but cardiolipin is exclusively found within the IMM. Though the lipid topology and distribution of the OMM and IMM are known since more than four decades, little is known about the interfacial and dynamic properties of the IMM and OMM lipid extracts. Here we build monolayers, supported bilayers and giant unilamellar vesicles (GUVs) of native OMM and IMM lipids extracts from porcine heart. Additionally, we perform a comparative analysis on the interfacial, phase immiscibility and mechanical properties of both types of extract. Our results show that IMM lipids form more expanded and softer membranes than OMM lipids, allowing a better understanding of the physicochemical and biophysical properties of mitochondrial membranes.

New Insights into the Interaction of Class II Dihydroorotate Dehydrogenases with Ubiquinone in Lipid Bilayers as a Function of Lipid Composition

The fourth enzymatic reaction in the de novo pyrimidine biosynthesis, the oxidation of dihydroorotate to orotate, is catalyzed by dihydroorotate dehydrogenase (DHODH). Enzymes belonging to the DHODH Class II are membrane-bound proteins that use ubiquinones as their electron acceptors. We have designed this study to understand the interaction of an N-terminally truncated human DHODH (HsΔ29DHODH) and the DHODH from Escherichia coli (EcDHODH) with ubiquinone (Q₁₀) in supported lipid membranes using neutron reflectometry (NR). NR has allowed us to determine in situ, under solution conditions, how the enzymes bind to lipid membranes and to unambiguously resolve the location of Q₁₀. Q₁₀ is exclusively located at the center of all of the lipid bilayers investigated, and upon binding, both of the DHODHs penetrate into the hydrophobic region of the outer lipid leaflet towards the Q₁₀. We therefore show that the interaction between the soluble enzymes and the membrane-embedded Q₁₀ is mediated by enzyme penetration. We can also show that EcDHODH binds more efficiently to the surface of simple bilayers consisting of 1-palmitoyl, 2-oleoyl phosphatidylcholine, and tetraoleoyl cardiolipin than HsΔ29DHODH, but does not penetrate into the lipids to the same degree. Our results also highlight the importance of Q₁₀, as well as lipid composition, on enzyme binding.

PoET: automated approach for measuring pore edge tension in giant unilamellar vesicles

Motivation

A reliable characterization of the membrane pore edge tension of single giant unilamellar vesicles (GUVs) requires the measurement of micrometer sized pores in hundreds to thousands of images. When manually performed, this procedure has shown to be extremely time-consuming and to generate inconsistent results among different users and imaging systems. A user-friendly software for such analysis allowing quick processing and generation of reproducible data had not yet been reported.

Results

We have developed a software (PoET) for automatic pore edge tension measurements on GUVs. The required image processing steps and the characterization of the pore dynamics are performed automatically within the software and its use allowed for a 30-fold reduction in the analysis time. We demonstrate the applicability of the software by comparing the pore edge tension of GUVs of different membrane compositions and surface charges. The approach was applied to electroporated GUVs but is applicable to other means of pore formation.

Availability and implementation

The complete software is implemented in Python and available for Windows at https://dx.doi.org/10.17617/3.7h.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Vesicle Viewer: Online visualization and analysis of small-angle scattering from lipid vesicles

Small-angle X-ray and neutron scattering are among the most powerful experimental techniques for investigating the structure of biological membranes. Much of the critical information contained in small-angle scattering (SAS) data is not easily accessible to researchers who have limited time to analyze results by hand or to nonexperts who may lack the necessary scientific background to process such data. Easy-to-use data visualization software can allow them to take full advantage of their SAS data and maximize the use of limited resources. To this end, we developed an internet-based application called Vesicle Viewer to visualize and analyze SAS data from unilamellar lipid bilayer vesicles. Vesicle Viewer utilizes a modified scattering density profile (SDP) analysis called EZ-SDP in which key bilayer structural parameters, such as area per lipid and bilayer thickness, are easily and robustly determined. Notably, we introduce a bilayer model that is able to describe an asymmetric bilayer, whether it be chemically or isotopically asymmetric. The application primarily uses Django, a Python package specialized for the development of robust web applications. In addition, several other libraries are used to support the more technical aspects of the project; notable examples are Matplotlib (for graphs) and NumPy (for calculations). By eliminating the barrier of downloading and installing software, this web-based application will allow scientists to analyze their own vesicle scattering data using their preferred operating system. The web-based application can be found at https://vesicleviewer.dmarquardt.ca/.

Cardiolipin prevents pore formation in phosphatidylglycerol bacterial membrane models

Several antimicrobial peptides, including magainin and the human cathelicidin LL-37, act by forming pores in bacterial membranes. Bacteria such as Staphylococcus aureus modify their membrane's cardiolipin composition to resist such types of perturbations that compromise their membrane stability. Here, we used molecular dynamic simulations to quantify the role of cardiolipin on the formation of pores in simple bacterial-like membrane models composed of phosphatidylglycerol and cardiolipin mixtures. Cardiolipin modified the structure and ordering of the lipid bilayer, making it less susceptible to mechanical changes. Accordingly, the free-energy barrier for the formation of a transmembrane pore and its kinetic instability augmented by increasing the cardiolipin concentration. This is attributed to the unfavorable positioning of cardiolipin near the formed pore, due to its small polar head and bulky hydrophobic body. Overall, our study demonstrates how cardiolipin prevents membrane-pore formation and this constitutes a plausible mechanism used by bacteria to act against stress perturbations and, thereby, gain resistance to antimicrobial agents.

Nanostructural Characterization of Cardiolipin-Containing Tethered Lipid Bilayers Adsorbed on Gold and Silicon Substrates for Protein Incorporation

A key to the development of lipid membrane-based devices is a fundamental understanding of how the molecular structure of the lipid bilayer membrane is influenced by the type of lipids used to build the membrane. This is particularly important when membrane proteins are included in these devices since the precise lipid environment affects the ability to incorporate membrane proteins and their functionality. Here, we used neutron reflectometry to investigate the structure of tethered bilayer lipid membranes and to characterize the incorporation of the NhaA sodium proton exchanger in the bilayer. The lipid membranes were composed of two lipids, dioleoyl phosphatidylcholine and cardiolipin, and were adsorbed on gold and silicon substrates using two different tethering architectures based on functionalized oligoethylene glycol molecules of different lengths. In all of the investigated samples, the addition of cardiolipin caused distinct structural rearrangement including crowding of ethylene glycol groups of the tethering molecules in the inner head region and a thinning of the lipid tail region. The incorporation of NhaA in the tethered bilayers following two different protocols is quantified, and the way protein incorporation modulates the structural properties of these membranes is detailed.

Mitochondrial Cristae Architecture and Functions: Lessons from Minimal Model Systems

Mitochondria are known as the powerhouse of eukaryotic cells. Energy production occurs in specific dynamic membrane invaginations in the inner mitochondrial membrane called cristae. Although the integrity of these structures is recognized as a key point for proper mitochondrial function, less is known about the mechanisms at the origin of their plasticity and organization, and how they can influence mitochondria function. Here, we review the studies which question the role of lipid membrane composition based mainly on minimal model systems.

Activation of Cytochrome C Peroxidase Function Through Coordinated Foldon Loop Dynamics upon Interaction with Anionic Lipids

Cardiolipin (CL) is a mitochondrial anionic lipid that plays important roles in the regulation and signaling of mitochondrial apoptosis. CL peroxidation catalyzed by the assembly of CL-cytochrome c (cyt c) complexes at the inner mitochondrial membrane is a critical checkpoint. The structural changes in the protein, associated with peroxidase activation by CL and different anionic lipids, are not known at a molecular level. To better understand these peripheral protein-lipid interactions, we compare how phosphatidylglycerol (PG) and CL lipids trigger cyt c peroxidase activation, and correlate functional differences to structural and motional changes in membrane-associated cyt c. Structural and motional studies of the bound protein are enabled by magic angle spinning solid state NMR spectroscopy, while lipid peroxidase activity is assayed by mass spectrometry. PG binding results in a surface-bound state that preserves a nativelike fold, which nonetheless allows for significant peroxidase activity, though at a lower level than binding its native substrate CL. Lipid-specific differences in peroxidase activation are found to correlate to corresponding differences in lipid-induced protein mobility, affecting specific protein segments. The dynamics of omega loops C and D are upregulated by CL binding, in a way that is remarkably controlled by the protein:lipid stoichiometry. In contrast to complete chemical denaturation, membrane-induced protein destabilization reflects a destabilization of select cyt c foldons, while the energetically most stable helices are preserved. Our studies illuminate the interplay of protein and lipid dynamics in the creation of lipid peroxidase-active proteolipid complexes implicated in early stages of mitochondrial apoptosis.

Native-like membrane models of E. coli polar lipid extract shed light on the importance of lipid composition complexity

BACKGROUND

Lipid-protein interactions stabilize protein oligomers, shape their structure, and modulate their function. Whereas in vitro experiments already account for the functional importance of lipids by using natural lipid extracts, in silico methods lack behind by embedding proteins in single component lipid bilayers. However, to accurately complement in vitro experiments with molecular details at very high spatio-temporal resolution, molecular dynamics simulations have to be performed in natural(-like) lipid environments.

RESULTS

To enable more accurate MD simulations, we have prepared four membrane models of E. coli polar lipid extract, a typical model organism, each at all-atom (CHARMM36) and coarse-grained (Martini3) representations. These models contain all main lipid headgroup types of the E. coli inner membrane, i.e., phosphatidylethanolamines, phosphatidylglycerols, and cardiolipins, symmetrically distributed between the membrane leaflets. The lipid tail (un)saturation and propanylation stereochemistry represent the bacterial lipid tail composition of E. coli grown at 37∘C until 3/4 of the log growth phase. The comparison of the Simple three lipid component models to the complex 14-lipid component model Avanti over a broad range of physiologically relevant temperatures revealed that the balance of lipid tail unsaturation and propanylation in different positions and inclusion of lipid tails of various length maintain realistic values for lipid mobility, membrane area compressibility, lipid ordering, lipid volume and area, and the bilayer thickness. The only Simple model that was able to satisfactory reproduce most of the structural properties of the complex Avanti model showed worse agreement of the activation energy of basal water permeation with the here performed measurements. The Martini3 models reflect extremely well both experimental and atomistic behavior of the E. coli polar lipid extract membranes. Aquaporin-1 embedded in our native(-like) membranes causes partial lipid ordering and membrane thinning in its vicinity. Moreover, aquaporin-1 attracts and temporarily binds negatively charged lipids, mainly cardiolipins, with a distinct cardiolipin binding site in the crevice at the contact site between two monomers, most probably stabilizing the tetrameric protein assembly.

CONCLUSIONS

The here prepared and validated membrane models of E. coli polar lipids extract revealed that lipid tail complexity, in terms of double bond and cyclopropane location and varying lipid tail length, is key to stabilize membrane properties over a broad temperature range. In addition, they build a solid basis for manifold future simulation studies on more realistic lipid membranes bridging the gap between simulations and experiments.

Curvature-Dependent Binding of Cytochrome c to Cardiolipin

Cytochrome c binds cardiolipin on the concave surface of the inner mitochondrial membrane, before oxidizing the lipid and initiating the apoptotic pathway. This interaction has been studied in vitro, where mimicking the membrane curvature of the binding environment is difficult. Here we report binding to concave, cardiolipin-containing, membrane surfaces and compare findings to convex binding under the same conditions. For binding to the convex outer surface of cardiolipin-containing vesicles, a two-step structural rearrangement is observed with a small rearrangement detectable by Soret circular dichroism (CD) occurring at an exposed lipid-to-protein ratio (LPR) near 10 and partial unfolding detectable by Trp59 fluorescence occurring at an exposed LPR near 23. On the concave inner surface of cardiolipin-containing vesicles, the structural transitions monitored by Soret CD and Trp59 fluorescence are coincident and occur at an exposed LPR near 58. On the concave inner surface of mitochondrial cristae, we estimate the LPR of cardiolipin to cytochrome c is between 50 and 100. Thus, cytochrome c may have adapted to its native environment so that it can undergo a conformational change that switches on its peroxidase activity when it binds to CL-containing membranes in the cristae early in apoptosis. Our results show that membrane curvature qualitatively affects peripheral protein-lipid interactions and also highlights the disparity between in vitro binding studies and their physiological counterparts where cone-shaped lipids, like cardiolipin, are involved.

Influence of External NaCl Salt on Membrane Rigidity of Neutral DOPC Vesicles

Sodium chloride (NaCl) is a very common molecule in biotic and abiotic aqueous environments. In both cases, variation of ionic strength is inevitable. In addition to the osmotic variation posed by such perturbations, the question of whether the interactions of monovalent ions Na⁺ and Cl^-, especially with the neutral head groups of phospholipid membranes are impactful enough to change the membrane rigidity, is still not entirely understood. We investigated the dynamics of 1,2-di-(octadecenoyl)-sn-glycero-3-phosphocholine (DOPC) vesicles with zwitterionic neutral head groups in the fluid phase with increasing external salt concentration. At higher salt concentrations, we observe an increase in bending rigidity from neutron spin echo (NSE) spectroscopy and an increase in bilayer thickness from small-angle X-ray scattering (SAXS). We compared different models to distinguish membrane undulations, lipid tail motions, and the translational diffusion of the vesicles. All of the models indicate an increase in bending rigidity by a factor of 1.3-3.6. We demonstrate that even down to t > 10 ns and for Q > 0.07 Å^-1, the observed NSE relaxation spectra are influenced by translational diffusion of the vesicles. For t < 5 ns, the lipid tail motion dominates the intermediate dynamic structure factor. As the salt concentration increases, this contribution diminishes. We introduced a time-dependent analysis for the bending rigidity that highlights only a limited Zilman-Granek time window in which the rigidity is physically meaningful.

Molecular Structure of Sphingomyelin in Fluid Phase Bilayers Determined by the Joint Analysis of Small-Angle Neutron and X-ray Scattering Data

We have determined the fluid bilayer structure of palmitoyl sphingomyelin (PSM) and stearoyl sphingomyelin (SSM) by simultaneously analyzing small-angle neutron and X-ray scattering data. Using a newly developed scattering density profile (SDP) model for sphingomyelin lipids, we report structural parameters including the area per lipid, total bilayer thickness, and hydrocarbon thickness, in addition to lipid volumes determined by densitometry. Unconstrained all-atom simulations of PSM bilayers at 55 °C using the C36 CHARMM force field produced a lipid area of 56 Ų, a value that is 10% lower than the one determined experimentally by SDP analysis (61.9 Ų). Furthermore, scattering form factors calculated from the unconstrained simulations were in poor agreement with experimental form factors, even though segmental order parameter (S_CD) profiles calculated from the simulations were in relatively good agreement with S_CD profiles obtained from NMR experiments. Conversely, constrained area simulations at 61.9 Ų resulted in good agreement between the simulation and experimental scattering form factors, but not with S_CD profiles from NMR. We discuss possible reasons for the discrepancies between these two types of data that are frequently used as validation metrics for molecular dynamics force fields.

The mitochondria-targeted peptide SS-31 binds lipid bilayers and modulates surface electrostatics as a key component of its mechanism of action

Mitochondrial dysfunction underlies many heritable diseases, acquired pathologies, and aging-related declines in health. Szeto-Schiller (SS) peptides comprise a class of amphipathic tetrapeptides that are efficacious toward a wide array of mitochondrial disorders and are believed to target mitochondrial membranes because they are enriched in the anionic phospholipid cardiolipin (CL). However, little is known regarding how SS peptides interact with or alter the physical properties of lipid bilayers. In this study, using biophysical and computational approaches, we have analyzed the interactions of the lead compound SS-31 (elamipretide) with model and mitochondrial membranes. Our results show that this polybasic peptide partitions into the membrane interfacial region with an affinity and a lipid binding density that are directly related to surface charge. We found that SS-31 binding does not destabilize lamellar bilayers even at the highest binding concentrations; however, it did cause saturable alterations in lipid packing. Most notably, SS-31 modulated the surface electrostatics of both model and mitochondrial membranes. We propose nonexclusive mechanisms by which the tuning of surface charge could underpin the mitoprotective properties of SS-31, including alteration of the distribution of ions and basic proteins at the interface, and/or modulation of bilayer physical properties. As a proof of concept, we show that SS-31 alters divalent cation (calcium) distribution within the interfacial region and reduces the energetic burden of calcium stress in mitochondria. The mechanistic details of SS-31 revealed in this study will help inform the development of future compound variants with enhanced efficacy and bioavailability.

Wheel and Deal in the Mitochondrial Inner Membranes: The Tale of Cytochrome c and Cardiolipin

Cardiolipin oxidation and degradation by different factors under severe cell stress serve as a trigger for genetically encoded cell death programs. In this context, the interplay between cardiolipin and another mitochondrial factor—cytochrome c—is a key process in the early stages of apoptosis, and it is a matter of intense research. Cytochrome c interacts with lipid membranes by electrostatic interactions, hydrogen bonds, and hydrophobic effects. Experimental conditions (including pH, lipid composition, and post-translational modifications) determine which specific amino acid residues are involved in the interaction and influence the heme iron coordination state. In fact, up to four binding sites (A, C, N, and L), driven by different interactions, have been reported. Nevertheless, key aspects of the mechanism for cardiolipin oxidation by the hemeprotein are well established. First, cytochrome c acts as a pseudoperoxidase, a process orchestrated by tyrosine residues which are crucial for peroxygenase activity and sensitivity towards oxidation caused by protein self-degradation. Second, flexibility of two weakest folding units of the hemeprotein correlates with its peroxidase activity and the stability of the iron coordination sphere. Third, the diversity of the mode of interaction parallels a broad diversity in the specific reaction pathway. Thus, current knowledge has already enabled the design of novel drugs designed to successfully inhibit cardiolipin oxidation.

A molecular view on the escape of lipoplexed DNA from the endosome

The use of non-viral vectors for in vivo gene therapy could drastically increase safety, whilst reducing the cost of preparing the vectors. A promising approach to non-viral vectors makes use of DNA/cationic liposome complexes (lipoplexes) to deliver the genetic material. Here we use coarse-grained molecular dynamics simulations to investigate the molecular mechanism underlying efficient DNA transfer from lipoplexes. Our computational fusion experiments of lipoplexes with endosomal membrane models show two distinct modes of transfection: parallel and perpendicular. In the parallel fusion pathway, DNA aligns with the membrane surface, showing very quick release of genetic material shortly after the initial fusion pore is formed. The perpendicular pathway also leads to transfection, but release is slower. We further show that the composition and size of the lipoplex, as well as the lipid composition of the endosomal membrane, have a significant impact on fusion efficiency in our models.

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF-Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born-Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. Here, using tetrahydrofuran (THF)-water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D₈THF to examine the molecular origins of the isotope-dependent loss of the THF-water miscibility gap.

Decoupling copolymer, lipid and carbon nanotube interactions in hybrid, biomimetic vesicles

Bilayer vesicles that mimic a real biological cell can be tailored to carry out a specific function by manipulating the molecular composition of the amphiphiles. These bio-inspired and bio-mimetic structures are increasingly being employed for a number of applications from drug delivery to water purification and beyond. Complex hybrid bilayers are the key building blocks for fully synthetic vesicles that can mimic biological cell membranes, which often contain a wide variety of molecular species. While the assembly and morpholgy of pure phospholid bilayer vesicles is well understood, the functionality and structure dramaticlly changes when copolymer and/or carbon nanotube porins (CNTP) are added. The aim of this study is to understand how the collective molecular interactions within hybrid vesicles affect their nanoscale structure and properties. In situ small and wide angle X-ray scattering (SAXS/WAXS) and molecular dynamics simulations (MD) are used to investigate the morphological effect of molecular interactions between polybutadiene polyethylene oxide, lipids and carbon nanotubes (CNT) within the hybrid vesicle bilayer. Within the lipid/copolymer system, the hybrid bilayer morphology transitions from phase separated lipid and compressed copolymer at low copolymer loadings to a mixed bilayer where opposing lipids are mostly separated from the inner region. This transition begins between 60 wt% and 70 wt%, with full homogenization observed by 80 wt% copolymer. The incorporation of CNT into the hybrid vesicles increases the bilayer thickness and enhances the bilayer symmetry. Analysis of the WAXS and MD indicate that the CNT-dioleoyl interactions are much stronger than the CNT-polybutadiene.

Antituberculosis Drug Interactions with Membranes: A Biophysical Approach Applied to Bedaquiline

This work focuses on the interaction of the novel and representative antituberculosis (anti-TB) drug bedaquiline (BDQ) with different membrane models of eukaryotic and prokaryotic cells. The effect of BDQ on eukaryotic cell membrane models was assessed using liposomes, namely, multilamellar vesicles (MLVs) made of 1,2-dimyristoyl-rac-glycero-3-phosphocholine (DMPC) and also a mixture of DMPC and cholesterol (CHOL) (8:2 molar ratio). To mimic the prokaryotic cell membrane, 1,2-dimyristoyl-sn-glycero-3-phospho-rac-(1-glycerol) (DMPG) and 1,1′2,2′-tetra-oleoyl-cardiolipin (TOCL) were chosen. Powerful biophysical techniques were employed, including small-angle X-ray scattering (SAXS) and wide-angle X-ray scattering (WAXS), to understand the effect of BDQ on the nanostructure of the membrane models. The results showed that BDQ demonstrated a pronounced disordering effect in the bacterial cell membrane models, especially in the membrane model with cardiolipin (CL), while the human cell membrane model with large fractions of neutral phospholipids remained less affected. The membrane models and techniques provide detailed information about different aspects of the drug–membrane interaction, thus offering valuable information to better understand the effect of BDQ on their target membrane-associated enzyme as well as its side effects on the cardiovascular system.

Cardiolipin-Dependent Properties of Model Mitochondrial Membranes from Molecular Simulations

Cardiolipin is an anionic lipid found in the mitochondrial membranes of eukaryotes ranging from unicellular microorganisms to metazoans. This unique lipid contributes to various mitochondrial functions, including metabolism, mitochondrial membrane fusion and/or fission dynamics, and apoptosis. However, differences in cardiolipin content between the two mitochondrial membranes, as well as dynamic fluctuations in cardiolipin content in response to stimuli and cellular signaling events, raise questions about how cardiolipin concentration affects mitochondrial membrane structure and dynamics. Although cardiolipin’s structural and dynamic roles have been extensively studied in binary mixtures with other phospholipids, the biophysical properties of cardiolipin in higher number lipid mixtures are still not well resolved. Here, we used molecular dynamics simulations to investigate the cardiolipin-dependent properties of ternary lipid bilayer systems that mimic the major components of mitochondrial membranes. We found that changes to cardiolipin concentration only resulted in minor changes to bilayer structural features but that the lipid diffusion was significantly affected by those alterations. We also found that cardiolipin position along the bilayer surfaces correlated to negative curvature deflections, consistent with the induction of negative curvature stress in the membrane monolayers. This work contributes to a foundational understanding of the role of cardiolipin in altering the properties in ternary lipid mixtures composed of the major mitochondrial phospholipids, providing much-needed insights to help understand how cardiolipin concentration modulates the biophysical properties of mitochondrial membranes.

Membrane curvature induces cardiolipin sorting

Cardiolipin is a cone-shaped lipid predominantly localized in curved membrane sites of bacteria and in the mitochondrial cristae. This specific localization has been argued to be geometry-driven, since the CL's conical shape relaxes curvature frustration. Although previous evidence suggests a coupling between CL concentration and membrane shape in vivo, no precise experimental data are available for curvature-based CL sorting in vitro. Here, we test this hypothesis in experiments that isolate the effects of membrane curvature in lipid-bilayer nanotubes. CL sorting is observed with increasing tube curvature, reaching a maximum at optimal CL concentrations, a fact compatible with self-associative clustering. Observations are compatible with a model of membrane elasticity including van der Waals entropy, from which a negative intrinsic curvature of -1.1 nm-1 is predicted for CL. The results contribute to understanding the physicochemical interplay between membrane curvature and composition, providing key insights into mitochondrial and bacterial membrane organization and dynamics.

Molecular Dynamics Analysis of Cardiolipin and Monolysocardiolipin on Bilayer Properties

The mitochondrial lipid cardiolipin (CL) contributes to the spatial protein organization and morphological character of the inner mitochondrial membrane. Monolysocardiolipin (MLCL), an intermediate species in the CL remodeling pathway, is enriched in the multisystem disease Barth syndrome. Despite the medical relevance of MLCL, a detailed molecular description that elucidates the structural and dynamic differences between CL and MLCL has not been conducted. To this end, we performed comparative atomistic molecular dynamics studies on bilayers consisting of pure CL or MLCL to elucidate similarities and differences in their molecular and bulk bilayer properties. We describe differential headgroup dynamics and hydrogen bonding patterns between the CL variants and show an increased cohesiveness of MLCL's solvent interfacial region, which may have implications for protein interactions. Finally, using the coarse-grained Martini model, we show that substitution of MLCL for CL in bilayers mimicking mitochondrial composition induces drastic differences in bilayer mechanical properties and curvature-dependent partitioning behavior. Together, the results of this work reveal differences between CL and MLCL at the molecular and mesoscopic levels that may underpin the pathomechanisms of defects in cardiolipin remodeling.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Lipid-modulation of membrane insertion and refolding of the apoptotic inhibitor Bcl-xL

Bcl-xL is a member of the Bcl-2 family of apoptotic regulators, responsible for inhibiting the permeabilization of the mitochondrial outer membrane, and a promising anti-cancer target. Bcl-xL exists in the following conformations, each believed to play a role in the inhibition of apoptosis: (a) a soluble folded conformation, (b) a membrane-anchored (by its C-terminal α8 helix) form, which retains the same fold as in solution and (c) refolded membrane-inserted conformations, for which no structural data are available. Previous studies established that in the cell Bcl-xL exists in a dynamic equilibrium between soluble and membranous states, however, no direct evidence exists in support of either anchored or inserted conformation of the membranous state in vivo. In this in vitro study, we employed a combination of fluorescence and EPR spectroscopy to characterize structural features of the bilayer-inserted conformation of Bcl-xL and the lipid modulation of its membrane insertion transition. Our results indicate that the core hydrophobic helix α6 inserts into the bilayer without adopting a transmembrane orientation. This insertion disrupts the packing of Bcl-xL and releases the regulatory N-terminal BH4 domain (α1) from the rest of the protein structure. Our data demonstrate that both insertion and refolding of Bcl-xL are modulated by lipid composition, which brings the apparent pK_a of insertion to the threshold of physiological pH. We hypothesize that conformational rearrangements associated with the bilayer insertion of Bcl-xL result in its switching to a so-called non-canonical mode of apoptotic inhibition. Presented results suggest that the alteration in lipid composition before and during apoptosis can serve as an additional factor regulating the permeabilization of the mitochondrial outer membrane.

Surface-Binding to Cardiolipin Nanodomains Triggers Cytochrome c Pro-apoptotic Peroxidase Activity via Localized Dynamics

The peroxidation of cardiolipins by reactive oxygen species, which is regulated and enhanced by cytochrome c (cyt c), is a critical signaling event in mitochondrial apoptosis. We probe the molecular underpinnings of this mitochondrial death signal through structural and functional studies of horse heart cyt c binding to mixed-lipid membranes containing cardiolipin with mono- and polyunsaturated acyl chains. Lipidomics reveal the selective oxidation of polyunsaturated fatty acid (PUFA) cardiolipin (CL), while multidimensional solid-state NMR probes the structure and dynamics of the membrane and the peripherally bound protein. The hydrophilic milieu at the membrane interface stabilizes a native-like fold, but also leads to localized flexibility at the membrane-interacting protein face. PUFA CL acts as both a preferred substrate and a dynamic regulator by affecting the dynamics of the cyt c N70-I85 Ω loop, which covers the heme cavity.

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

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity (κ_C) and area compressibility (K_A). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate K_A values, current approaches suffer significant technical limitations. Here, we present a novel, to our knowledge, computational framework that allows for a direct estimation of K_A values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives K_A from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model's predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.

Could Cardiolipin Protect Membranes against the Action of Certain Antimicrobial Peptides? Aurein 1.2, a Case Study

The activity of a host of antimicrobial peptides has been examined against a range of lipid bilayers mimicking bacterial and eukaryotic membranes. Despite this, the molecular mechanisms and the nature of the physicochemical properties underlying the peptide-lipid interactions that lead to membrane disruption are yet to be fully elucidated. In this study, the interaction of the short antimicrobial peptide aurein 1.2 was examined in the presence of an anionic cardiolipin-containing lipid bilayer using molecular dynamics simulations. Aurein 1.2 is known to interact strongly with anionic lipid membranes. In the simulations, the binding of aurein 1.2 was associated with buckling of the lipid bilayer, the degree of which varied with the peptide concentration. The simulations suggest that the intrinsic properties of cardiolipin, especially the fact that it promotes negative membrane curvature, may help protect membranes against the action of peptides such as aurein 1.2 by counteracting the tendency of the peptide to induce positive curvature in target membranes.

Phospholipid Bilayer Softening Due to Hydrophobic Gold Nanoparticle Inclusions

Liposome-nanoparticle assemblies (LNAs) are vital in the context of novel targeted drug-delivery systems, in addition to investigating nanoparticle-lipid bilayer interactions. Quantifying membrane structural properties and dynamics in presence of nanoparticle inclusions provides a simple model to elucidate nanoparticle effects on membrane biophysical properties. We present experimental evidences of bilayer softening due to small hydrophobic gold nanoparticle inclusions. LNA structure has been investigated by a combination of cryo-transmission electron microscopy, dynamic light scattering, and small-angle neutron scattering. Neutron spin echo spectroscopy demonstrated a remarkable ∼15% bending modulus decrease for LNAs relative to pure liposomes. Clear dependence of bending modulus on gold nanoparticle diameter and concentration was observed from our observations. Our findings point toward local bilayer fluidization by nanoparticle inclusions leading to an overall bilayer softening. These findings add valuable information to liposomal drug-delivery vehicle design and membrane biophysics research.

Determination of bending rigidity and tilt modulus of lipid membranes from real-space fluctuation analysis of molecular dynamics simulations

We have recently developed a novel computational methodology (termed RSF for Real-Space Fluctuations) to quantify the bending rigidity and tilt modulus of lipid membranes from real-space analysis of fluctuations in the tilt and splay degrees of freedom as sampled in molecular dynamics (MD) simulations. In this article, we present a comprehensive study that combines results from the application of the RSF method to a wide range of lipid bilayer systems that encompass membranes of different fluidities and sizes, including lipids with saturated and unsaturated lipid tails, single and multi-component lipid systems, as well as non-standard lipids such as the four-tailed cardiolipin. By comparing the material properties calculated with the RSF method to those obtained from experimental data and from other computational methodologies, we rigorously demonstrate the validity of our approach and show its robustness. This should allow for future applications of even more complex lipidic assemblies, whose material properties are not tractable by other computational techniques. In addition, we discuss the relationship between different definitions of the tilt modulus appearing in current literature to address some important unresolved discrepancies in the field.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.