Molecular structures of fluid phase phosphatidylglycerol bilayers as determined by small angle neutron and X-ray scattering

Microfluidic Formulation of DNA-Loaded Multicomponent Lipid Nanoparticles for Gene Delivery

In recent years, lipid nanoparticles (LNPs) have gained considerable attention in numerous research fields ranging from gene therapy to cancer immunotherapy and DNA vaccination. While some RNA-encapsulating LNP formulations passed clinical trials, DNA-loaded LNPs have been only marginally explored so far. To fulfil this gap, herein we investigated the effect of several factors influencing the microfluidic formulation and transfection behavior of DNA-loaded LNPs such as PEGylation, total flow rate (TFR), concentration and particle density at the cell surface. We show that PEGylation and post-synthesis sample concentration facilitated formulation of homogeneous and small size LNPs with high transfection efficiency and minor, if any, cytotoxicity on human Embryonic Kidney293 (HEK-293), spontaneously immortalized human keratinocytes (HaCaT), immortalized keratinocytes (N/TERT) generated from the transduction of human primary keratinocytes, and epidermoid cervical cancer (CaSki) cell lines. On the other side, increasing TFR had a detrimental effect both on the physicochemical properties and transfection properties of LNPs. Lastly, the effect of particle concentration at the cell surface on the transfection efficiency (TE) and cell viability was largely dependent on the cell line, suggesting that its case-by-case optimization would be necessary. Overall, we demonstrate that fine tuning formulation and microfluidic parameters is a vital step for the generation of highly efficient DNA-loaded LNPs.

Stable DOPG/Glycyrrhizin Vesicles with a Wide Range of Mixing Ratios: Structure and Stability as Seen by Scattering Experiments and Cryo-TEM

Phosphatidylglycerols represent a large share of the lipids in the plasmamembrane of procaryotes. Therefore, this study investigates the role of charged lipids in the plasma membrane with respect to the interaction of the antiviral saponin glycyrrhizin with such membranes. Glycyrrhizin is a natural triterpenic-based surfactant found in licorice. Vesicles made of 1,2-dioleoyl-sn-glycero-3-phospho-rac-(1’-glycerol) (DOPG)/glycyrrhizin are characterized by small-angle scattering with neutrons and X-rays (SANS and SAXS). Small-angle scattering data are first evaluated by the model-independent modified Kratky–Porod method and afterwards fitted by a model describing the shape of small unilamellar vesicles (SUV) with an internal head-tail contrast. Complete miscibility of DOPG and glycyrrhizin was revealed even at a ratio of lipid:saponin of 1:1. Additional information about the chain-chain correlation distance of the lipid/saponin mixtures in the SUV structures is obtained from wide-angle X-ray scattering (WAXS).

Ionic environment, thickness and line tension as determinants of phase separation in whole Purified Myelin Membranes monolayers

Purified myelin membranes (PMM) were spread as monomolecular films at the air/aqueous solution interface, and visualized by Brewster Angle Microscopy (BAM) at different lateral pressures (π) on three specific aqueous solutions: absence of salts, physiological conditions and presence of calcium. Coexistence of Liquid-Expanded (LE) and Liquid Ordered (LO) phases persisted up to collapse in the presence of salts, whereas monolayers became homogeneous at π ≥ 35-40 mN/m when salts are absent. This PMM phase-mixing behavior in monolayers is similar to the previously reported behavior of PMM multilamellar vesicles. Reflectivities (R_p) of p-polarized light from both phases were assessed throughout the whole π -range, and film thicknesses (t) were calculated from the R_p values and measured film refractive indices (n). The LO phase was found to be more reflective and thicker than the LE phase at π ≤ 15 mN/m, but less reflective and thinner at higher π. We also determined the line tension (λ) of PMM monolayers at the domain boundaries from the rate of domain shape relaxation, which turned out to be of the order of picoNewtons (pN) and decreased as π increased. A correlation between λ and thickness differences (Δt) was found, suggesting that Δt is a molecular determinant for λ in PMM monolayers. Both λ and Δt were found to increase markedly when calcium was present in the subphase. This result corroborates the concept of divalent cations as a stabilizing factor for phase separation, in line with earlier studies on this mixture forming multilamellar membrane arrangements.

Structure and Interdigitation of Chain-Asymmetric Phosphatidylcholines and Milk Sphingomyelin in the Fluid Phase

We addressed the frequent occurrence of mixed-chain lipids in biological membranes and their impact on membrane structure by studying several chain-asymmetric phosphatidylcholines and the highly asymmetric milk sphingomyelin. Specifically, we report trans-membrane structures of the corresponding fluid lamellar phases using small-angle X-ray and neutron scattering, which were jointly analyzed in terms of a membrane composition-specific model, including a headgroup hydration shell. Focusing on terminal methyl groups at the bilayer center, we found a linear relation between hydrocarbon chain length mismatch and the methyl-overlap for phosphatidylcholines, and a non-negligible impact of the glycerol backbone-tilting, letting the sn1-chain penetrate deeper into the opposing leaflet by half a CH2 group. That is, penetration-depth differences due to the ester-linked hydrocarbons at the glycerol backbone, previously reported for gel phase structures, also extend to the more relevant physiological fluid phase, but are significantly reduced. Moreover, milk sphingomyelin was found to follow the same linear relationship suggesting a similar tilt of the sphingosine backbone. Complementarily performed molecular dynamics simulations revealed that there is always a part of the lipid tails bending back, even if there is a high interdigitation with the opposing chains. The extent of this back-bending was similar to that in chain symmetric bilayers. For both cases of adaptation to chain length mismatch, chain-asymmetry has a large impact on hydrocarbon chain ordering, inducing disorder in the longer of the two hydrocarbons.

Calcium mediated DNA binding in non-lamellar structures formed by DOPG/glycerol monooleate

In order to test an encapsulation method of short fragmented DNA (∼ 20-300 bp), we study the solubilisation in 150 mM solution of NaCl of a cubic phase formed by glycerol monooleate (GMO) with negatively charged dioleoylphosphatidylglycerol (DOPG) up to the level of unilamellar vesicles and, subsequently, the restoration of the cubic phase using Ca²⁺ cations. We performed small angle X-ray and neutron scattering (SAXS and SANS) to follow structural changes in DOPG/GMO mixtures induced by increasing DOPG content. The cubic phase (Pn3m space group) is preserved up to ∼ 11 mol% of DOPG in DOPG/GMO. Above 20 mol%, the SANS curves are typical of unilamellar vesicles. The thickness of the DOPG/GMO lipid bilayer (d_L) decreases slightly with increasing fraction of DOPG. The addition of 15 mM of CaCl₂ solution shields the electrostatic repulsions of DOPG molecules, increases slightly d_L and restores the cubic structures in the mixtures up to ∼ 37 mol% of DOPG. Zeta potential shows negative surface charge. The analysis of the data provides the radius of the water nano-channels of the formed non-lamellar structures. We discuss their dimensions with respect to DNA binding. In addition, Ca²⁺ mediates DNA - DOPG/GMO binding. The formed hexagonal phase, H_II, binds less of DNA in comparison with cubic phases (∼ 6 wt% and ∼ 20 wt% of the total amount, respectively). The studied system can be utilized as anionic Q_II delivery vector for genetic material.

To Close or to Collapse: The Role of Charges on Membrane Stability upon Pore Formation

Resealing of membrane pores is crucial for cell survival. Membrane surface charge and medium composition are studied as defining regulators of membrane stability. Pores are generated by electric field or detergents. Giant vesicles composed of zwitterionic and negatively charged lipids mixed at varying ratios are subjected to a strong electric pulse. Interestingly, charged vesicles appear prone to catastrophic collapse transforming them into tubular structures. The spectrum of destabilization responses includes the generation of long-living submicroscopic pores and partial vesicle bursting. The origin of these phenomena is related to the membrane edge tension, which governs pore closure. This edge tension significantly decreases as a function of the fraction of charged lipids. Destabilization of charged vesicles upon pore formation is universal-it is also observed with other poration stimuli. Disruption propensity is enhanced for membranes made of lipids with higher degree of unsaturation. It can be reversed by screening membrane charge in the presence of calcium ions. The observed findings in light of theories of stability and curvature generation are interpreted and mechanisms acting in cells to prevent total membrane collapse upon poration are discussed. Enhanced membrane stability is crucial for the success of electroporation-based technologies for cancer treatment and gene transfer.

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.

Membrane Interactions of Virus-like Mesoporous Silica Nanoparticles

In the present study, we investigated lipid membrane interactions of silica nanoparticles as carriers for the antimicrobial peptide LL-37 (LLGDFFRKSKEKIGKEFKRIVQRIKDFLRNLVPRTES). In doing so, smooth mesoporous nanoparticles were compared to virus-like mesoporous nanoparticles, characterized by a "spiky" external surface, as well as to nonporous silica nanoparticles. For this, we employed a combination of neutron reflectometry, ellipsometry, dynamic light scattering, and ζ-potential measurements for studies of bacteria-mimicking bilayers formed by palmitoyloleoylphosphatidylcholine/palmitoyloleoylphosphatidylglycerol. The results show that nanoparticle topography strongly influences membrane binding and destabilization. We found that virus-like particles are able to destabilize such lipid membranes, whereas the corresponding smooth silica nanoparticles are not. This effect of particle spikes becomes further accentuated after loading of such particles with LL-37. Thus, peptide-loaded virus-like nanoparticles displayed more pronounced membrane disruption than either peptide-loaded smooth nanoparticles or free LL-37. The structural basis of this was clarified by neutron reflectometry, demonstrating that the virus-like nanoparticles induce trans-membrane defects and promote incorporation of LL-37 throughout both bilayer leaflets. The relevance of such effects of particle spikes for bacterial membrane rupture was further demonstrated by confocal microscopy and live/dead assays on Escherichia coli bacteria. Taken together, these findings demonstrate that topography influences the interaction of nanoparticles with bacteria-mimicking lipid bilayers, both in the absence and presence of antimicrobial peptides, as well as with bacteria. The results also identify virus-like mesoporous nanoparticles as being of interest in the design of nanoparticles as delivery systems for antimicrobial peptides.

Overlapping Properties of the Short Membrane-Active Peptide BP100 With (i) Polycationic TAT and (ii) α-helical Magainin Family Peptides

BP100 is a short, designer-made membrane-active peptide with multiple functionalities: antimicrobial, cell-penetrating, and fusogenic. Consisting of five lysines and 6 hydrophobic residues, BP100 was shown to bind to lipid bilayers as an amphipathic α-helix, but its mechanism of action remains unclear. With these features, BP100 embodies the characteristics of two distinctly different classes of membrane-active peptides, which have been studied in detail and where the mechanism of action is better understood. On the one hand, its amphiphilic helical structure is similar to the pore forming magainin family of antimicrobial peptides, though BP100 is much too short to span the membrane. On the other hand, its length and high charge density are reminiscent of the HIV-TAT family of cell penetrating peptides, for which inverted micelles have been postulated as translocation intermediates, amongst other mechanisms. Assays were performed to test the antimicrobial and hemolytic activity, the induced leakage and fusion of lipid vesicles, and cell uptake. From these results the functional profiles of BP100, HIV-TAT, and the magainin-like peptides magainin 2, PGLa, MSI-103, and MAP were determined and compared. It is observed that the activity of BP100 resembles most closely the much longer amphipathic α-helical magainin-like peptides, with high antimicrobial activity along with considerable fusogenic and hemolytic effects. In contrast, HIV-TAT shows almost no antimicrobial, fusogenic, or hemolytic effects. We conclude that the amphipathic helix of BP100 has a similar membrane-based activity as magainin-like peptides and may have a similar mechanism of action.

CHARMM36 Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Phosphatidylethanolamine, Phosphatidylglycerol, and Ether Lipids

Long-range Lennard-Jones (LJ) interactions have been incorporated into the CHARMM36 (C36) lipid force field (FF) using the LJ particle-mesh Ewald (LJ-PME) method in order to remove the inconsistency of bilayer and monolayer properties arising from the exclusion of long-range dispersion [Yu, Y.; Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion. J. Chem. Theory Comput. 2021, 10.1021/acs.jctc.0c01326. (preceding article in this issue)]. The new FF is denoted C36/LJ-PME. While the first optimization was based on three phosphatidylcholines (PCs), this work extends the validation and parametrization to more lipids including PC, phosphatidylethanolamine (PE), phosphatidylglycerol (PG), and ether lipids. The agreement with experimental structure data is excellent for PC, PE, and ether lipids. C36/LJ-PME also compares favorably with scattering data of PG bilayers but less so with NMR deuterium order parameters of 1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DMPG) at 303.15 K, indicating a need for future optimization regarding PG-specific parameters. Frequency dependence of NMR T1 spin-lattice relaxation times is well-described by C36/LJ-PME, and the overall agreement with experiment is comparable to C36. Lipid diffusion is slower than C36 due to the added long-range dispersion causing a higher viscosity, although it is still too fast compared to experiment after correction for periodic boundary conditions. When using a 10 Å real-space cutoff, the simulation speed of C36/LJ-PME is roughly equal to C36. While more lipids will be incorporated into the FF in the future, C36/LJ-PME can be readily used for common lipids and extends the capability of the CHARMM FF by supporting monolayers and eliminating the cutoff dependence.

Composition effects on photooxidative membrane destabilization by TiO2 nanoparticles

Membrane interactions and photooxidative membrane destabilization of titanium dioxide (TiO₂) nanoparticles were investigated, focusing on the effects of membrane composition, notably phospholipid headgroup charge and presence of cholesterol. For this, we employed a battery of state-of-the-art methods for studies of bilayers formed by zwitterionic palmitoyloleoylphosphatidylcholine (POPC) containing also polyunsaturated palmitoylarachidonoylphosphocholine (PAPC), as well as its mixtures with anionic palmitoyloleoylphosphatidylglycerol (POPG) and cholesterol. It was found that the TiO₂ nanoparticles display close to zero charge at pH 7.4, resulting in aggregation. At pH 3.4, in contrast, the 6 nm TiO₂ nanoparticles are well dispersed due to a strongly positive ζ-potential. Mirroring this pH dependence, TiO₂ nanoparticles were observed to bind to negatively charged lipid bilayers at pH 3.4, but much less so at pH 7.4. While nanoparticle binding has some destabilizing effect alone, illumination with ultraviolet (UV) light accentuates membrane destabilization, a result of oxidative stress caused by generated reactive oxygen species (ROS). Neutron reflectivity (NR), quartz crystal microbalance (QCM), and small-angle X-ray scattering (SAXS) results all demonstrate that membrane composition strongly influences membrane interactions and photooxidative destabilization of lipid bilayers. In particular, the presence of anionic POPG makes the bilayers more sensitive to oxidative destabilization, whereas a stabilizing effect was observed in the presence of cholesterol. Also, structural aspects of peroxidation were found to depend strongly on membrane composition, notably the presence of anionic phospholipids. The results show that membrane interactions and UV-induced ROS generation act in concert and need to be considered together to understand effects of lipid membrane composition on UV-triggered oxidative destabilization by TiO₂ nanoparticles, e.g., in the context of oxidative damage of bacteria and cells.

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.

Understanding the Structure and Dynamics of Complex Biomembrane Interactions by Neutron Scattering Techniques

The membrane is one of the key structural materials of biology at the cellular level. Composed predominantly of a bilayer of lipids with embedded and bound proteins, it defines the boundaries of the cell and many organelles essential to life and therefore is involved in almost all biological processes. Membrane-specific interactions, such as drug binding to a membrane receptor or the interactions of an antimicrobial compound with the lipid matrix of a pathogen membrane, are of interest across the scientific disciplines. Herein we present a review, aimed at nonexperts, of the major neutron scattering techniques used in membrane studies: small-angle neutron scattering, neutron membrane diffraction, neutron reflectometry, quasielastic neutron scattering, and neutron spin echo. Neutron scattering techniques are well suited to studying biological membranes. The nondestructive nature of cold neutrons means that samples can be measured for long periods without fear of beam damage from ultraviolet, electron, or X-ray radiation, and neutron beams are highly penetrating, thus offering flexibility in samples and sample environments. Most important is the strong difference in neutron scattering lengths between the two most abundant forms of hydrogen, protium and deuterium. Changing the relative amounts of protium/deuterium in a sample allows the production of a series of neutron scattering data sets, enabling the observation of differing components within complex membrane architectures. This approach can be as simple as using the naturally occurring neutron contrast between different biomolecules to study components in a complex by changing the solution H₂O/D₂O ratio or as complex as selectively labeling individual components with hydrogen isotopes. This review presents an overview of each experimental technique with the neutron instrument configuration, related sample preparation and sample environment, and data analysis, highlighted by a special emphasis on using prominent neutron contrast to understand structure and dynamics. This review gives researchers a practical introduction to the often enigmatic suite of neutron beamlines, thereby lowering the barrier to taking advantage of these large-facility techniques to achieve new understandings of membranes and their interactions with other molecules.

Constant-pH Brownian Dynamics Simulations of a Protein near a Charged Surface

We have developed a rigid-body Brownian dynamics algorithm that allows for simulations of a globular protein suspended in an ionic solution confined by a charged planar boundary, with an explicit treatment of pH-dependent protein protonation equilibria and their couplings to the electrostatic potential of the plane. Electrostatic interactions are described within a framework of the continuum Poisson-Boltzmann model, whereas protein-plane hydrodynamic interactions are evaluated based on analytical expressions for the position- and orientation-dependent near-wall friction tensor of a spheroid. The algorithm was applied to simulate near-surface diffusion of lysozyme in solutions having pH in the range 4-10 and ionic strengths of 10 and 150 mM. As a reference, we performed Brownian dynamics simulations in which the protein is assigned a fixed, most probable protonation state, appropriate for given solution conditions and unaffected by the presence of the charged plane, and Brownian dynamics simulations in which the protein probes possible protonation states with the pH-dependent probability, but these variations are not coupled to the electric field generated by the boundary. We show that electrostatic interactions with the negatively charged plane substantially modify probabilities of different protonation states of lysozyme and shift protonation equilibria of both acidic and basic amino acid side chains toward higher pH values. Consequently, equilibrium energy distributions, equilibrium position-orientation distributions, and functions that characterize rotational dynamics, which for a protein with multiple ionization sites, such as lysozyme, in the presence of a charged obstacle are pH-dependent, are significantly affected by the approach taken to incorporate the solution pH into simulations.

Fluctuation spectroscopy of giant unilamellar vesicles using confocal and phase contrast microscopy

A widely used method to measure the bending rigidity of bilayer membranes is fluctuation spectroscopy, which analyses the thermally-driven membrane undulations of giant unilamellar vesicles recorded with either phase-contrast or confocal microscopy. Here, we analyze the fluctuations of the same vesicle using both techniques and obtain consistent values for the bending modulus. We discuss the factors that may lead to discrepancies.

Lipid membrane interactions of self-assembling antimicrobial nanofibers: effect of PEGylation

Supramolecular assembly and PEGylation (attachment of a polyethylene glycol polymer chain) of peptides can be an effective strategy to develop antimicrobial peptides with increased stability, antimicrobial efficacy and hemocompatibility. However, how the self-assembly properties and PEGylation affect their lipid membrane interaction is still an unanswered question. In this work, we use state-of-the-art small angle X-ray and neutron scattering (SAXS/SANS) together with neutron reflectometry (NR) to study the membrane interaction of a series of multidomain peptides, with and without PEGylation, known to self-assemble into nanofibers. Our approach allows us to study both how the structure of the peptide and the membrane are affected by the peptide-lipid interactions. When comparing self-assembled peptides with monomeric peptides that are not able to undergo assembly due to shorter chain length, we found that the nanofibers interact more strongly with the membrane. They were found to insert into the core of the membrane as well as to absorb as intact fibres on the surface. Based on the presented results, PEGylation of the multidomain peptides leads to a slight net decrease in the membrane interaction, while the distribution of the peptide at the interface is similar to the non-PEGylated peptides. Based on the structural information, we showed that nanofibers were partially disrupted upon interaction with phospholipid membranes. This is in contrast with the considerable physical stability of the peptide in solution, which is desirable for an extended in vivo circulation time.

Update of the CHARMM36 United Atom Chain Model for Hydrocarbons and Phospholipids

Accurate lipid force field (FF) parameters used in molecular dynamics (MD) simulations are crucial for understanding the properties of lipid-containing systems and biological processes related to lipids. The last update of the CHARMM36 united atom chain model (C36UA) was in 2013 [Lee, S. J. Phys. Chem. B 2014, 118, 547 556]; it utilized CHARMM36 (C36) lipid FF parameters for headgroups and OPLS-UA Lennard-Jones (LJ) parameters for tails. Simulations with the FF were able to reproduce many experimental observables of lipid bilayers accurately, but to be more applicable for a wide range of lipids, additional FF parameter optimization was needed. In this work, we present an update of the model, named C36UAr. The parameterization included the LJ parameters for hydrocarbons and related dihedrals. Bulk liquid properties (density, heat of vaporization, isothermal compressibility, and diffusion constant) of model compounds were used as targets for the LJ parameter fitting, and dihedrals were fit to either quantum mechanical (QM) or potential of mean force (PMF) calculations using C36. Thermodynamic reweighting was used to further improve the parameters. Bilayer simulations of various lipid headgroups (phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol) and tails (saturated, monounsaturated, and polyunsaturated) were performed to validate the model, and significant improvements were seen in bilayer properties, including surface area, membrane thicknesses, NMR deuterium order parameters, and density profiles. C36UAr was also compared to the hydrogen mass repartitioning (HMR) method. The high accuracy and competitive efficiency shown in this study make C36UAr one of the best choices for studies of membrane structure and membrane-associated proteins.

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 structures of polyunsaturated lipid bilayers by joint refinement of neutron and X-ray scattering data

We present the detailed structural analysis of polyunsaturated fatty acid-containing phospholipids namely, 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PDPC) and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC). A newly developed molecular dynamics (MD) simulation parsing scheme for lipids containing fatty acids with multiple double bonds was implemented into the scattering density profile (SDP) model to simultaneously refine differently contrasted neutron and X-ray scattering data. SDP analyses of scattering data at 30 °C yielded lipid areas of 71.1 Ų and 70.4 Ų for PDPC and SDPC bilayers, respectively, and a model free analysis of PDPC at 30 °C resulted in a lipid area of 72 Ų. In addition to bilayer structural parameters, using area-constrained MD simulations we determined the area compressibility modulus, K_A, to be 246.4 mN/m, a value similar to other neutral phospholipids.

Delayed nucleation in lipid particles

Metastable states in first-order phase-transitions have been traditionally described by classical nucleation theory (CNT). However, recently an increasing number of systems displaying such a transition have not been successfully modelled by CNT. The delayed crystallization of phospholipids upon super-cooling is an interesting case, since the extended timescales allow access into the dynamics. Herein, we demonstrate the controllable behavior of the long-lived metastable liquid-crystalline phase of dilauroyl-phosphatidylethanolamine (DLPE), arranged in multi-lamellar vesicles, and the ensuing cooperative transition to the crystalline state. Experimentally, we find that the delay in crystallization is a bulk phenomenon, which is tunable and can be manipulated to span two orders of magnitude in time by changing the quenching temperature, solution salinity, or adding a secondary phospholipid. Our results reveal the robust persistence of the metastability, and showcase the apparent deviation from CNT. This distinctive suppression of the transition may be explained by the resistance of the multi-lamellar vesicle to deformations caused by nucleated crystalline domains. Since phospholipids are used as a platform for drug-delivery, a programmable design of cargo hold and release can be of great benefit.

How Well Can You Tailor the Charge of Lipid Vesicles?

Knowledge and control of surface charge or potential is important for tailoring colloidal interactions. In this work, we compare widely used zeta potential (ζ) measurements of charged lipid vesicle surface potential to direct measurements using the surface force apparatus (SFA). Our measurements show good agreement between the two techniques. On varying the fraction of anionic lipids dimyristoylphosphatidylserine (DMPS) or dimyristoylphosphatidylglycerol (DMPG) mixed with zwitterionic dimyristoylphosphatidylcholine (DMPC) from 0 to 100 mol % we observed a near-linear increase in membrane surface charge or potential up to 20-30 mol % charged lipids beyond which charge saturation occurred in physiological (high) salt conditions. Similarly, in low salt concentrations, a linear increase in charge/potential was found but only up to ∼5-10 mol % charged lipids beyond which the surface charge or potential leveled off. While a lower degree of ionization is expected due to the lower dielectric constant (ε ∼ 4) of the lipid acyl chain environment, increasing intramembrane electrostatic repulsion between neighboring charged lipid head groups at higher charge loading contributes to charge suppression. Measured potentials in physiological salt solutions were consistent with predictions using the Gouy-Chapman-Stern-Grahame (GCSG) model of the electrical double layer with Langmuir binding of counterions, but in low salt conditions, the model significantly overestimated the surface charge/potential. The much lower ionization in low salt (maximum ∼1-2% of total lipids ionized) instead was consistent with counterion condensation at the bilayer surface which limited the charge that could be obtained. The strong interplay between membrane composition, lipid headgroup ionization, electrolyte concentration, and solution pH complicates exact prediction and tuning of membrane surface charge for applications. However, the theoretical frameworks used here can provide guidelines to understand this interplay and establish a range of achievable potentials for a system and predict the response to triggers like pH and salt concentration changes.

On the Mechanism of Bilayer Separation by Extrusion, or Why Your LUVs Are Not Really Unilamellar

Extrusion through porous filters is a widely used method for preparing biomimetic model membranes. Of primary importance in this approach is the efficient production of single bilayer (unilamellar) vesicles that eliminate the influence of interlamellar interactions and strictly define the bilayer surface area available to external reagents such as proteins. Submicroscopic vesicles produced using extrusion are widely assumed to be unilamellar, and large deviations from this assumption would impact interpretations from many model membrane experiments. Using three probe-free methods-small angle X-ray and neutron scattering and cryogenic electron microscopy-we report unambiguous evidence of extensive multilamellarity in extruded vesicles composed of neutral phosphatidylcholine lipids, including for the common case of neutral lipids dispersed in physiological buffer and extruded through 100-nm diameter pores. In such preparations, only ∼35% of lipids are externally accessible and this fraction is highly dependent on preparation conditions. Charged lipids promote unilamellarity as does decreasing solvent ionic strength, indicating the importance of electrostatic interactions in determining the lamellarity of extruded vesicles. Smaller extrusion pore sizes also robustly increase the fraction of unilamellar vesicles, suggesting a role for membrane bending. Taken together, these observations suggest a mechanistic model for extrusion, wherein the formation of unilamellar vesicles involves competition between bilayer bending and adhesion energies. The findings presented here have wide-ranging implications for the design and interpretation of model membrane studies, especially ensemble-averaged observations relying on the assumption of unilamellarity.

Lipid Head Group Parameterization for GROMOS 54A8: A Consistent Approach with Protein Force Field Description

Membranes are a crucial component of both bacterial and mammalian cells, being involved in signaling, transport, and compartmentalization. This versatility requires a variety of lipid species to tailor the membrane’s behavior as needed, increasing the complexity of the system. Molecular dynamics simulations have been successfully applied to study model membranes and their interactions with proteins, elucidating some crucial mechanisms at the atomistic detail and thus complementing experimental techniques. An accurate description of the functional interplay of the diverse membrane components crucially depends on the selected parameters that define the adopted force field. A coherent parameterization for lipids and proteins is therefore needed. In this work, we propose and validate new lipid head group parameters for the GROMOS 54A8 force field, making use of recently published parametrizations for key chemical moieties present in lipids. We make use additionally of a new canonical set of partial charges for lipids, chosen to be consistent with the parameterization of soluble molecules such as proteins. We test the derived parameters on five phosphocholine model bilayers, composed of lipid patches four times larger than the ones used in previous studies, and run 500 ns long simulations of each system. Reproduction of experimental data like area per lipid and deuterium order parameters is good and comparable with previous parameterizations, as well as the description of liquid crystal to gel-phase transition. On the other hand, the orientational behavior of the head groups is more realistic for this new parameter set, and this can be crucial in the description of interactions with other polar molecules. For that reason, we tested the interaction of the antimicrobial peptide lactoferricin with two model membranes showing that the new parameters lead to a weaker peptide–membrane binding and give a more realistic outcome in comparing binding to antimicrobial versus mammal membranes.

Bending rigidity of charged lipid bilayer membranes

We experimentally investigate the effect of lipid charge on the stiffness of bilayer membranes. The bending rigidity of membranes with composition 0-100 mol% of charged lipids, in the absence and presence of salt at different concentrations, is measured with the flicker spectroscopy method, using the shape fluctuations of giant unilamellar vesicles. The analysis considers both the mean squared amplitudes and the time autocorrelations of the shape modes. Our results show that membrane charge increases the bending rigidity relative to the charge-free membrane. The effect is diminished by the addition of monovalent salt to the suspending solutions. The trend shown by the membrane bending rigidity correlates with zeta potential measurements, confirming charge screening at different salt concentrations. The experimental results in the presence of salt are in good agreement with existing theories of membrane stiffening by surface charge.

PACKMOL-Memgen: A Simple-To-Use, Generalized Workflow for Membrane-Protein-Lipid-Bilayer System Building

We present PACKMOL-Memgen, a simple-to-use, generalized workflow for automated building of membrane-protein-lipid-bilayer systems based on open-source tools including Packmol, memembed, pdbremix, and AmberTools. Compared with web-interface-based related tools, PACKMOL-Memgen allows setup of multiple configurations of a system in a user-friendly and efficient manner within minutes. The generated systems are well-packed and thus well-suited as starting configurations in MD simulations under periodic boundary conditions, requiring only moderate equilibration times. PACKMOL-Memgen is distributed with AmberTools and runs on most computing platforms, and its output can also be used for CHARMM or adapted to other molecular-simulation packages.

An introduction to classical molecular dynamics simulation for experimental scattering users

Classical molecular dynamics simulations are a common component of multi-modal analyses of scattering measurements, such as small-angle scattering and diffraction. Users of these experimental techniques often have no formal training in the theory and practice of molecular dynamics simulation, leading to the possibility of these simulations being treated as a 'black box' analysis technique. This article describes an open educational resource (OER) designed to introduce classical molecular dynamics to users of scattering methods. This resource is available as a series of interactive web pages, which can be easily accessed by students, and as an open-source software repository, which can be freely copied, modified and redistributed by educators. The topics covered in this OER include classical atomistic modelling, parameterizing interatomic potentials, molecular dynamics simulations, typical sources of error and some of the approaches to using simulations in the analysis of scattering data.

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.

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.

Assessing the Interactions of Auristatin Derivatives with Mixed Phospholipid-Sodium Dodecyl Sulfate Aggregate Dispersions

The aim of this study was to assess what properties of the pseudostationary phases in electrokinetic capillary chromatography affect the interactions between monomethyl auristatin E (MMAE) and hydrophilically modified structural analogues thereof with various lipophilic phases. MMAE is a widely used cytotoxic agent in antibody-drug conjugates (ADC), which are used as selective biopharmaceutical drugs in the treatment of cancers. MMAE and its derivatives are highly lipophilic, yet they fail to interact with biomimicking phosphatidylcholine-phosphatidylserine liposomes. To reveal what properties affect the interaction of the auristatin derivatives with cell plasma membrane-mimicking vesicles, capillary electrokinetic chromatography was used with four different types of micellar and vesicular pseudostationary phases: pure vesicles, mixed vesicles, mixed micelles, and pure micelles. Vesicular phases were composed of pure phospholipids [dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC)] and phospholipid-surfactant mixtures [sodium dodecyl sulfate, (SDS) with DMPC and DLPC] while the micellar phases comprised pure surfactant (SDS) and surfactant-phospholipid mixtures (SDS-DMPC and SDS-DLPC). In addition, differential scanning calorimetry and dynamic light scattering were used to monitor the aggregate composition. Our data shows that the interaction between hydrophobic auristatin derivatives and hydrophobic pseudostationary phases critically depends on the type, size, and hydrogen bonding capability of the pseudostationary phases.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Membrane-Protein-Hydration Interaction of α-Synuclein with Anionic Vesicles Probed via Angle-Resolved Second-Harmonic Scattering

Amyloid formation of the protein α-synuclein promotes neurodegeneration in Parkinson's disease. The normal function of α-synuclein includes synaptic vesicle transport and fusion, and the protein binds strongly to negatively charged vesicles in vitro. Here, we demonstrate that nonresonant angle-resolved second-harmonic scattering detects α-synuclein binding to liposomes through changes in water orientational correlations and can thus be used as a high-accuracy and high-throughput label-free probe of protein-liposome interactions. The obtained results support a binding model in which the N-terminus of α-synuclein adopts an α-helical conformation that lies flat on the vesicle surface while the negatively charged C-terminus remains in solution.

Nanoparticles at Biomimetic Interfaces: Combined Experimental and Simulation Study on Charged Gold Nanoparticles/Lipid Bilayer Interfaces

The poor understanding of the interaction of nanomaterials with biologically relevant interfaces is recognized as one of the major issues currently limiting the development of nanomedicine. The central purpose of this study is to compare experimental (confocal microscopy, fluorescence correlation spectroscopy, X-ray reflectivity) and computational (molecular dynamics simulations) results to thoroughly describe the interaction of cationic gold nanoparticles (AuNPs) with mixed zwitterionic/anionic lipid membranes. The adhesion of AuNPs to the lipid membrane is investigated on different length scales from a structural and dynamical point of view; with this approach, a series of complex phenomena, spanning from lipid extraction, localized membrane disruption, lateral phase separation, and slaved diffusion, are characterized and interpreted from a molecular level to macroscopic observations.

Bayesian determination of the effect of a deep eutectic solvent on the structure of lipid monolayers

A novel reflectometry analysis method reveals the structure of lipid monolayers at the air-DES interface.

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.

Probing the mechanism of bupivacaine drug release from multivesicular liposomes

The mechanism of drug release from complex dosage forms, such as multivesicular liposomes (MVLs), is complex and oftentimes sensitive to the release environment. This challenges the design and development of an appropriate in vitro release test (IVRT) method. In this study, a commercial bupivacaine MVL product was selected as a model product and an IVRT method was developed using a modified USP 2 apparatus in conjunction with reverse-dialysis membranes. This setup allowed the use of in situ UV-Vis probes to continuously monitor the drug concentration during release. In comparison to the traditional sample-and-separate methods, the new method allowed for better control of the release conditions allowing for study of the drug release mechanism. Bupivacaine (BPV) MVLs exhibited distinct tri-phasic release characteristics comprising of an initial burst release, lag phase and a secondary release. Temperature, pH, agitation speed and release media composition were observed to impact the mechanism and rate of BPV release from MVLs. The size and morphology of the MVLs as well as their inner vesicle compartments were analyzed using cryogenic-scanning electron microscopy (cryo-SEM), confocal laser scanning microscopy and laser diffraction, where the mean diameters of the MVLs and their inner "polyhedral" vesicles were found to be 23.6 ± 11.5 μm and 1.52 ± 0.44 μm, respectively. Cryo-SEM results further showed a decrease in particle size and loss of internal "polyhedral" structure of the MVLs over the duration of release, indicating erosion and rearrangement of the lipid layers. Based on these results a potential MVL drug release mechanism was proposed, which may assist with the future development of more biorelevant IVRT method for similar formulations.

Lipid composition in fungal membrane models: effect of lipid fluidity

The creation of effective fungal membrane models for neutron and X-ray reflectometry experiments is a key step in the development of new antifungal pharmaceuticals and agrochemicals to allow in vitro investigation of their mode of interaction with target cells. The structure of the obtained models depends on the properties of the lipids used and the final composition of the leaflets, and can be subject to the spontaneous translocation of phospholipids across the bilayer. The effect of phospholipid acyl-chain unsaturation and the presence of steroids in the membrane on the bilayer asymmetry were examined by means of neutron reflectometry. The measurements showed that membrane stability was higher if a zwitterionic, saturated acyl-chain phospholipid is present as the inner leaflet. Furthermore, membrane asymmetry was higher in the case of fully saturated lipid systems. As a result, membrane models consisting of fully saturated acyl chains within the inner leaflet are recommended as the starting point for subsequent studies of antifungal interactions owing to the simplicity of the models and their relative stability, thus allowing better control over the exact lipid composition facing the tested antifungal.

Resolving the structural interactions between antimicrobial peptides and lipid membranes using small-angle scattering methods: the case of indolicidin

Using small angle X-ray and neutron scattering (SAXS/SANS) and detailed theoretical modelling we have elucidated the structure of the antimicrobial peptide, indolicidin, and the interaction with model lipid membranes of different anionic lipid compositions mimicking typical charge densities found in the cytoplasmic membrane of bacteria. First, we show that indolicidin displays a predominantly disordered, random chain conformation in solution with a small fraction (≈1%) of fiber-like nanostructures that are not dissolved at higher temperatures. The peptide is shown to strongly interact with the membranes at all charge densities without significantly perturbing the lipid bilayer structure. Instead, the results show that indolicidin inserts into the outer leaflet of the lipid vesicles causing a reduced local order of the lipid packing. This result is supported by an observed change in the melting point of the lipids upon addition of the peptide, as seen by differential scanning calorimetry experiments. The peptide does not to our observation affect the thickness of the membrane or form distinct structural pores in the membrane at physiologically relevant concentrations as has been previously suggested as an important mode of action. Finally, using sophisticated contrast variation SANS, we show that the peptide does not affect the random lateral distribution of anionic lipids in the membrane. Together, these results demonstrate that the structural aspects of the mode of action of antimicrobial peptides can be elucidated in detail using SAS techniques with liposomes as model systems.

Molecular dynamics modeling of Pseudomonas aeruginosa outer membranes

Pseudomonas aeruginosa is a common Gram-negative bacterium and opportunistic human pathogen. The distinctive structure of its outer membrane (OM) and outer membrane vesicles (OMVs) plays a fundamental role in bacterial virulence, colonization ability, and antibiotic resistance. To provide critical insights into OM and OMV functionality, we conducted an all-atom molecular dynamics study of asymmetric membranes that are biologically relevant to P. aeruginosa. We hybridized a GLYCAM06-based lipopolysaccharides force field with the Stockholm lipids force field (Slipids) to model bilayer membranes with Lipid A molecules in one leaflet and physiologically relevant phospholipid molecules in the other, including 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE), 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPG), and 1,2-dioleoyl-sn-glycero-3-phosphoglycerol (DOPG). In particular, a membrane with phospholipid composition representing the P. aeruginosa OM was constructed and modeled by mixing the physiologically dominant components. The detailed structure of membranes was characterized by area per lipid, transmembrane mass and charge densities, radial distribution function (RDF), deuterium order parameter (SCD) of acyl chains, and inclination angles of phosphates and disaccharide in Lipid A. The membrane fluidity in equilibrium and the hydration of functional groups were probed and characterized quantitatively. The consistent properties of the Lipid A leaflets in different membranes demonstrate its compatibility with various phospholipids present in the P. aeruginosa OM. The more ordered acyl chains of Lipid A compared to the cytoplasmic cell membrane contribute to the low permeability of bacterial outer membrane. The findings of this computational investigation of P. aeruginosa OM will further the understanding of microbial pathogenesis and enable future study of OMV biogenesis.

Characterization of self-assembled hybrid siloxane-phosphocholine bilayers

We have synthesized six new hybrid siloxane phosphocholines (SiPCs) and examined their self-assembly behaviour in aqueous dispersions. Employing small angle X-ray scattering we have characterized SiPC bilayers. SiPCs exhibit differential self-assembly behaviour that results from the interplay between the siloxane fatty acid in the sn-2 position and the differing chain length fatty acids in the sn-1 position. SiPCs that possess a fatty acid chain of a C₈-C₁₄ chain length in the sn-1 position form unilamellar vesicles. Extending the fatty acid chain length to C₁₆ and C₁₈ allows for the formation of both unilamellar and multilamellar vesicles. We propose that the preferential formation of unilamellar vesicles is the result of an enhanced hydrophobic effect imparted by siloxane chains at the termini of lipid tails.

Charged giant unilamellar vesicles prepared by electroformation exhibit nanotubes and transbilayer lipid asymmetry

Giant unilamellar vesicles (GUVs) are increasingly used as a versatile research tool to investigate membrane structure, morphology and phase state. In these studies, GUV preparation is typically enhanced by an externally applied electric field, a process called electroformation. We find that upon osmotic deflation, GUVs electroformed from charged and neutral lipids exhibit inward pointing lipid nanotubes, suggesting negative spontaneous curvature of the membrane. By quenching a fluorescent analog of the charged lipid, zeta potential measurements and experiments with the lipid marker annexin A5, we show that electroformed GUVs exhibit an asymmetric lipid distribution across the bilayer leaflets. The asymmetry is lost either after storing electroformed GUVs at room temperature for one day or by applying higher voltages and temperatures during electroformation. GUVs having the same lipid composition but grown via gel-assisted swelling do not show asymmetric lipid distribution. We discuss possible mechanisms for the generation and relaxation of lipid asymmetry, as well as implications for studies using electroformed vesicles. The observed effects allow to control the molecular assembly of lipid bilayer leaflets. Vesicle tubulation as reported here is an example of protein-free reshaping of membranes and is caused by compositional lipid asymmetry between leaflets.