X‐ray structural studies of freely suspended ordered hydrated DMPC multimembrane films

Cryo-TEM Reveals the Influence of Multivalent Charge and PEGylation on Shape Transitions in Fluid Lipid Assemblies: From Vesicles to Discs, Rods, and Spheres

Lipids, and cationic lipids in particular are of interest as delivery vectors for hydrophobic drugs such as the cancer therapeutic paclitaxel, and the structures of lipid assemblies affect their efficacy. We investigated the effect of incorporating the multivalent cationic lipid MVL5 (+5e) and poly(ethylene glycol)-lipids (PEG-lipids), alone and in combination, on the structure of fluid-phase lipid assemblies of the charge-neutral lipid 1,2-dioleoyl-sn-glycero-phosphocholine (DOPC). This allowed us to elucidate lipid-assembly structure correlations in sonicated formulations with high charge density, which are not accessible with univalent lipids such as the well-studied DOTAP (+1e). Cryogenic transmission electron microscopy (cryo-TEM) allowed us to determine the structure of the lipid assemblies, revealing diverse combinations of vesicles and disc-shaped, worm-like, and spherical micelles. Remarkably, MVL5 forms an essentially pure phase of disc micelles at 50 mol % MVL5. At a higher (75 mol %) content of MVL5, short- and intermediate-length worm-like micellar rods were observed, and in ternary mixtures with PEG-lipid, longer and highly flexible worm-like micelles formed. Independent of their length, the worm-like micelles coexisted with spherical micelles. In stark contrast, DOTAP forms mixtures of vesicles, disc micelles, and spherical micelles at all studied compositions, even when combined with PEG-lipids. The observed similarities and differences in the effects of charge (multivalent versus univalent) and high curvature (multivalent charge versus PEG-lipid) on the assembly structure provide insights into parameters that control the size of fluid lipid nanodiscs, relevant for future applications.

Insight into structural biophysics from solution X-ray scattering

The current challenges of structural biophysics include determining the structure of large self-assembled complexes, resolving the structure of ensembles of complex structures and their mass fraction, and unraveling the dynamic pathways and mechanisms leading to the formation of complex structures from their subunits. Modern synchrotron solution X-ray scattering data enable simultaneous high-spatial and high-temporal structural data required to address the current challenges of structural biophysics. These data are complementary to crystallography, NMR, and cryo-TEM data. However, the analysis of solution scattering data is challenging; hence many different analysis tools, listed in the SAS Portal (http://smallangle.org/), were developed. In this review, we start by briefly summarizing classical X-ray scattering analyses providing insight into fundamental structural and interaction parameters. We then describe recent developments, integrating simulations, theory, and advanced X-ray scattering modeling, providing unique insights into the structure, energetics, and dynamics of self-assembled complexes. The structural information is essential for understanding the underlying physical chemistry principles leading to self-assembled supramolecular architectures and computational structural refinement.

Paclitaxel-Loaded Cationic Fluid Lipid Nanodiscs and Liposomes with Brush-Conformation PEG Chains Penetrate Breast Tumors and Trigger Caspase-3 Activation

Novel approaches are required to address the urgent need to develop lipid-based carriers of paclitaxel (PTX) and other hydrophobic drugs for cancer chemotherapy. Carriers based on cationic liposomes (CLs) with fluid (i.e., chain-melted) membranes (e.g., EndoTAG-1) have shown promise in preclinical and late-stage clinical studies. Recent work found that the addition of a cone-shaped poly(ethylene glycol)-lipid (PEG-lipid) to PTX-loaded CLs (CLsPTX) promotes a transition to sterically stabilized, higher-curvature (smaller) nanoparticles consisting of a mixture of PEGylated CLsPTX and PTX-containing fluid lipid nanodiscs (nanodiscsPTX). These CLsPTX and nanodiscsPTX show significantly improved uptake and cytotoxicity in cultured human cancer cells at PEG coverage in the brush regime (10 mol % PEG-lipid). Here, we studied the PTX loading, in vivo circulation half-life, and biodistribution of systemically administered CLsPTX and nanodiscsPTX and assessed their ability to induce apoptosis in triple-negative breast-cancer-bearing immunocompetent mice. We focused on fluid rather than solid lipid nanodiscs because of the significantly higher solubility of PTX in fluid membranes. At 5 and 10 mol % of a PEG-lipid (PEG5K-lipid, molecular weight of PEG 5000 g/mol), the mixture of PEGylated CLsPTX and nanodiscsPTX was able to incorporate up to 2.5 mol % PTX without crystallization for at least 20 h. Remarkably, compared to preparations containing 2 and 5 mol % PEG5K-lipid (with the PEG chains in the mushroom regime), the particles at 10 mol % (with PEG chains in the brush regime) showed significantly higher blood half-life, tumor penetration, and proapoptotic activity. Our study suggests that increasing the PEG coverage of CL-based drug nanoformulations can improve their pharmacokinetics and therapeutic efficacy.

Polymer-Lipid Hybrid Materials

Hierarchic self-assembly underpins much of the form and function seen in synthetic or biological soft materials. Lipids are paramount examples, building themselves in nature or synthetically in a variety of meso/nanostructures. Synthetic block copolymers capture many of lipid's structural and functional properties. Lipids are typically biocompatible and high molecular weight polymers are mechanically robust and chemically versatile. The development of new materials for applications like controlled drug/gene/protein delivery, biosensors, and artificial cells often requires the combination of lipids and polymers. The emergent composite material, a "polymer-lipid hybrid membrane", displays synergistic properties not seen in pure components. Specific examples include the observation that hybrid membranes undergo lateral phase separation that can correlate in registry across multiple layers into a three-dimensional phase-separated system with enhanced permeability of encapsulated drugs. It is timely to underpin these emergent properties in several categories of hybrid systems ranging from colloidal suspensions to supported hybrid films. In this review, we discuss the form and function of a vast number of polymer-lipid hybrid systems published to date. We rationalize the results to raise new fundamental understanding of hybrid self-assembling soft materials as well as to enable the design of new supramolecular systems and applications.

Effect of pH on the phase behavior of DMPC bilayers

We have studied the effect of acidic pH on the phase behavior of the zwitterionic lipid 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) using differential scanning calorimetry and x-ray scattering. Dispersions of DMPC in HCl solutions of pH = 4 and 3 behave identical to dispersions in water. The main transition temperature increases sharply and the pre-transition disappears at lower pH. An untilted gel phase is observed at pH = 2 and 1, in contrast to the tilted gel phase found at higher pH. The relatively large periodicity of the untilted gel phase, in comparison to that of the tilted gel phase occurring near neutral pH, clearly demonstrates the simultaneous charging and dehydration of the headgroups as the pH approaches the pK of the phosphate group. Headgroup dehydration at low pH also leads to the formation of DMPC crystallites and the inverted hexagonal phase at low and high temperatures, respectively, after a few days of incubation. These results show the significant effect of acidic pH on the phase behavior of zwitterionic lipids.

Two-fluid model for the breakdown of flow alignment in nematic liquid crystals

We present a macroscopic two-fluid model to explain the breakdown of flow alignment in nematic liquid crystals under shear flow due to smectic clusters. We find that the velocity difference of the two fluids plays a key role to mediate the time-dependent behavior as soon as a large enough amount of smectic order is induced by flow. For the minimal model it is sufficient to keep the nematic degrees of freedom, the mass density of the smectic clusters and the degree of smectic order, the density, and two velocities as macroscopic variables. While frequently a smectic A or C phase arises at lower temperatures, this is not required for the applicability of the present model. Indeed, as pointed out before by Gähwiller, there are compounds showing a breakdown of flow alignment over a large temperature range and no smectic phase, but a solid phase at lower temperatures. We also demonstrate that, using a one velocity model, there is no coupling under shear flow between induced smectic order and the director orientation in stationary situations thus rendering such a model to be unsuitable to describe the breakdown of flow alignment. In a two-fluid description, flow alignment breaks down and becomes unstable with regard to a space- and time-dependent state due to an induced finite velocity difference. In an Appendix we outline a mesoscopic model to account for the sign change in the anisotropy of the electric conductivity observed in nematics with smectic clusters.

Fluid-fluid coexistence in phospholipid membranes induced by decanol

We have observed fluid-fluid coexistence in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membrane containing 1-decanol, using different experimental techniques and membrane morphologies. This phase behavior is reversible and occurs over a temperature range just above the chain melting transition temperature of the membrane. Although earlier experimental studies and computer simulations have shown the ability of decanol to enhance lipid chain ordering, its potential to induce fluid-fluid coexistence in membranes has not been hitherto recognized. Being the only binary membrane system known so far to exhibit fluid-fluid coexistence, the present system can serve as a simple model to gain a better understanding of mechanisms that drive this unusual phase behavior, which is believed to play an important role in the functioning of cell membranes.

Behavior of Hydrated Lipid Bilayers at Cryogenic Temperatures

Neutron diffraction was used to study the behavior of water present in phospholipid multilamellar stacks from 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) at cryogenic temperatures. Evidence was found for the existence of a highly viscous phase of water that exists between 180 and 220 K based on the observation that water can leave the intermembrane space at these low temperatures. Similar measurements are described in the literature for purple membrane (PM) samples. From a comparison with results from this natural membrane by using the same flash-cooling protocol, it is found that in the case of pure lipid samples, less water is trapped and the water flows out at lower temperatures. This suggests that the water is less hindered in its movements than in the PM case. It is shown that at least the Lβ′-phase of DMPC can be trapped likely by flash cooling; upon heating to about 260 K, it transforms to another phase that was not fully characterized.

Examining Tail and Headgroup Effects on Binary and Ternary Gel-Phase Lipid Bilayer Structure

The structural properties of two- and three-component gel-phase bilayers were studied using molecular dynamics simulations. The bilayers contain distearoylphosphatidylcholine (DSPC) phospholipids mixed with alcohols and/or fatty acids of varying tail lengths, with carbon chain lengths of 12, 16, and 24 studied. Changes in both headgroup chemistry and tail length are found to affect the balance between steric repulsion and van der Waals attraction within the bilayers, manifesting in different bilayer structural properties. Lipid components are found to be located at different depths within the bilayer depending on both chain length and headgroup chemistry. The highest bilayer ordering and lowest area per tail are found in systems with medium-length tails. While longer tails can enhance van der Waals attractions, the increased tail-length asymmetry is found to induce disorder and reduce tail packing. Bulkier headgroups further increase steric repulsion, as reflected in increased component offsets and reduced tail packing. These findings help explain how bilayer composition affects the structure of gel-phase bilayers.

X-ray scattering reveals molecular tilt is an order parameter for the main phase transition in a model biomembrane

Synchrotron diffuse x-ray scattering data reveal a dramatic softening of the molecular tilt modulus K_{θ} of the model biomembrane composed of DMPC lipids as the temperature is lowered towards the main phase transition temperature at T_{M}=24^{∘}C. Spontaneous tilt occurs below T_{M}, suggesting that tilt is a symmetry breaking order parameter. Consistent with this hypothesis, it is also found that a different lipid POPS has no spontaneous tilt below its T_{M} at 14^{∘}C and correspondingly its tilt modulus did not soften as T_{M} was approached from above. As previously known, the bending modulus K_{C} of DMPC also softens close to T_{M}, but unlike the tilt modulus, K_{C} has a maximum 3^{∘} above T_{M}, which also marks the limit of the well-known anomalous swelling regime. Tilt adds a different perspective to our previous understanding of the main phase transition in lipid bilayers.

Dynamical Behavior of Hydration Water Molecules between Phospholipid Membranes

The dynamical behavior of hydration water sandwiched between 1,2-dimyristyl-sn-glycero-3-phosphocholine (DMPC) bilayers was investigated by quasi-elastic neutron scattering (QENS) in the range between 275 and 316 K, where the main transition temperature of DMPC is interposed. The results revealed that the hydration water could be categorized into three types of water: (1) free water, whose dynamical behavior is slightly different from that of bulk water; (2) loosely bound water, whose dynamical behavior is 1 order of magnitude slower than that of the free water; and (3) tightly bound water, whose dynamical behavior is comparable with that of DMPC molecules. The number of loosely bound and tightly bound water molecules per DMPC molecule monotonically decreased and increased with decreasing temperature, respectively, and the sum of these water molecules remained constant. The number of free water molecules per DMPC molecule was constant in the measured temperature range.

Hydration forces between aligned DNA helices undergoing B to A conformational change: In-situ X-ray fiber diffraction studies in a humidity and temperature controlled environment

Hydration forces between DNA molecules in the A- and B-Form were studied using a newly developed technique enabling simultaneous in situ control of temperature and relative humidity. X-ray diffraction data were collected from oriented calf-thymus DNA fibers in the relative humidity range of 98%-70%, during which DNA undergoes the B- to A-form transition. Coexistence of both forms was observed over a finite humidity range at the transition. The change in DNA separation in response to variation in humidity, i.e. change of chemical potential, led to the derivation of a force-distance curve with a characteristic exponential decay constant of∼2Å for both A- and B-DNA. While previous osmotic stress measurements had yielded similar force-decay constants, they were limited to B-DNA with a surface separation (wall-to-wall distance) typically>5Å. The current investigation confirms that the hydration force remains dominant even in the dry A-DNA state and at surface separation down to∼1.5Å, within the first hydration shell. It is shown that the observed chemical potential difference between the A and B states could be attributed to the water layer inside the major and minor grooves of the A-DNA double helices, which can partially interpenetrate each other in the tightly packed A phase. The humidity-controlled X-ray diffraction method described here can be employed to perform direct force measurements on a broad range of biological structures such as membranes and filamentous protein networks.

Comparative Study of the Condensing Effects of Ergosterol and Cholesterol

Cholesterol, due to its condensing effect, is considered an important regulator of membrane thickness. Other sterols, due to their structural similarities to cholesterol, are often assumed to have a universal effect on membrane properties similar to the condensing effect of cholesterol, albeit possibly to different degrees. We used x-ray diffraction to investigate this assumption. By the combination of lamellar diffraction and grazing-angle scattering, we measured the membrane thickness and the tilt-angle distribution of the lipid's hydrocarbon chains. This method is sensitive to phase separation, which is important for examining the miscibility of sterols and phospholipids. Mixtures of ergosterol or cholesterol with dimyristoylphosphatidylcholine, palmitoyloleoylphosphatidylcholine, and dioleoylphosphatidylcholine were systematically studied. We found that mixing ergosterol with phospholipids into a single phase became increasingly difficult with higher sterol concentrations and also with higher concentrations of unsaturated lipid chains. The only condensing effect of ergosterol was found in dimyristoylphosphatidylcholine, although the effect was less than one-third of the effect of cholesterol. Unlike cholesterol, ergosterol could not maintain a fixed electron density profile of the surrounding lipids independent of hydration. In dioleoylphosphatidylcholine and palmitoyloleoylphosphatidylcholine, ergosterol made the membranes thinner, opposite to the effect of cholesterol. In all cases, the tilt-angle variation of the chain diffraction was consistent with the membrane thickness changes measured by lamellar diffraction, i.e., a thickening was always associated with a reduction of chain tilt angles. Our findings do not support the notion that different sterols have a universal behavior that differs only in degree.

Effect of cholesterol on the molecular structure and transitions in a clinical-grade lung surfactant extract

The lipid-protein film covering the interface of the lung alveolar in mammals is vital for proper lung function and its deficiency is related to a range of diseases. Here we present a molecular-level characterization of a clinical-grade porcine lung surfactant extract using a multitechnique approach consisting of [Formula: see text]-[Formula: see text] solid-state nuclear magnetic spectroscopy, small- and wide-angle X-ray scattering, and mass spectrometry. The detailed characterization presented for reconstituted membranes of a lung extract demonstrates that the molecular structure of lung surfactant strongly depends on the concentration of cholesterol. If cholesterol makes up about 11% of the total dry weight of lung surfactant, the surfactant extract adopts a single liquid-ordered lamellar phase, [Formula: see text], at physiological temperatures. This [Formula: see text] phase gradually changes into a liquid-disordered lamellar phase, [Formula: see text], when the temperature is increased by a few degrees. In the absence of cholesterol the system segregates into one lamellar gel phase and one [Formula: see text] phase. Remarkably, it was possible to measure a large set of order parameter magnitudes [Formula: see text] from the liquid-disordered and -ordered lamellar phases and assign them to specific C-H bonds of the phospholipids in the biological extract with no use of isotopic labeling. These findings with molecular details on lung surfactant mixtures together with the presented NMR methodology may guide further development of pulmonary surfactant pharmaceuticals that better mimic the physiological self-assembly compositions for treatment of pathological states such as respiratory distress syndrome.

Localization and Ordering of Lipids Around Aquaporin-0: Protein and Lipid Mobility Effects

Hydrophobic matching, lipid sorting, and protein oligomerization are key principles by which lipids and proteins organize in biological membranes. The Aquaporin-0 channel (AQP0), solved by electron crystallography (EC) at cryogenic temperatures, is one of the few protein-lipid complexes of which the structure is available in atomic detail. EC and room-temperature molecular dynamics (MD) of dimyristoylglycerophosphocholine (DMPC) annular lipids around AQP0 show similarities, however, crystal-packing and temperature might affect the protein surface or the lipids distribution. To understand the role of temperature, lipid phase, and protein mobility in the localization and ordering of AQP0-lipids, we used MD simulations of an AQP0-DMPC bilayer system. Simulations were performed at physiological and at DMPC gel-phase temperatures. To decouple the protein and lipid mobility effects, we induced gel-phase in the lipids or restrained the protein. We monitored the lipid ordering effects around the protein. Reducing the system temperature or inducing lipid gel-phase had a marginal effect on the annular lipid localization. However, restraining the protein mobility increased the annular lipid localization around the whole AQP0 surface, resembling EC. The distribution of the inter-phosphate and hydrophobic thicknesses showed that stretching of the DMPC annular layer around AQP0 surface is the mechanism that compensates the hydrophobic mismatch in this system. The distribution of the local area-per-lipid and the acyl-chain order parameters showed particular fluid- and gel-like areas that involved several lipid layers. These areas were in contact with the surfaces of higher and lower protein mobility, respectively. We conclude that the AQP0 surfaces induce specific fluid- and gel-phase prone areas. The presence of these areas might guide the AQP0 lipid sorting interactions with other membrane components, and is compatible with the squared array oligomerization of AQP0 tetramers separated by a layer of annular lipids.

Pre-transition effects mediate forces of assembly between transmembrane proteins

We present a mechanism for a generic, powerful force of assembly and mobility for transmembrane proteins in lipid bilayers. This force is a pre-transition (or pre-melting) effect for the first-order transition between ordered and disordered phases in the membrane. Using large-scale molecular simulation, we show that a protein with hydrophobic thickness equal to that of the disordered phase embedded in an ordered bilayer stabilizes a microscopic order–disorder interface. The stiffness of that interface is finite. When two such proteins approach each other, they assemble because assembly reduces the net interfacial energy. Analogous to the hydrophobic effect, we refer to this phenomenon as the 'orderphobic effect'. The effect is mediated by proximity to the order–disorder phase transition and the size and hydrophobic mismatch of the protein. The strength and range of forces arising from this effect are significantly larger than those that could arise from membrane elasticity for the membranes considered.

DOI:http://dx.doi.org/10.7554/eLife.13150.001

The membrane that surrounds cells provides a selective barrier that allows some molecules through, but blocks the path of others. A cell’s membrane is made up of two layers of molecules with oily tails, and is therefore known as a bilayer. Many proteins are dotted within and on the inner and outer surfaces of the bilayer: some act as channels that control what goes in and out of the cell, while others protrude outside the cell so that they can sense changes in the environment.

Membrane proteins can move and interact within the bilayer, and various models have emerged to try to explain this dynamic system. These models are based on the membrane having some fluidity but also having regions where there is more structure, and typically describe the proteins as ordered clusters floating in an otherwise disordered fluid membrane. However, many researchers now think some proteins that pass through both layers of the bilayer (i.e., transmembrane proteins) make membranes more ordered, with a possibly gel-like state. However, it is not clear how transmembrane proteins can move and assemble together within such a relatively rigid membrane.

To investigate this, Katira, Mandadapu, Vaikuntanathan et al. carried out computer simulations using a model of a simple bilayer membrane. This membrane can exist in an ordered state, where the oily tails are neatly aligned, or a disordered state, where they are irregularly packed. Virtual ‘heating’ of the membrane caused it to shift from an ordered to a disordered state. When a simple transmembrane protein favoring the disordered state was inserted into the ordered state of the modeled membrane, disordered regions formed locally around the protein and the protein was able to diffuse within the membrane.

Modeling what would happen if two transmembrane proteins approached each other revealed that a consequence of the order–disorder transition is a strong attractive force that assembles the proteins together. Katira, Mandadapu, Vaikuntanathan et al. named this new phenomenon the 'orderphobic effect'. The forces arising from this effect were much greater than those currently believed to contribute to the assembly of membrane protein complexes, such as those generated by the elasticity of the membrane. This means that the orderphobic effect may be responsible for generating the protein clusters commonly seen in cell membranes.

Future work should next explore the opposite effect, where proteins favoring the ordered state are inserted into the disordered state of a membrane. This is expected to cause clustering of such proteins and thus large ordered regions in an otherwise disordered membrane.

DOI:http://dx.doi.org/10.7554/eLife.13150.002

Reparameterized United Atom Model for Molecular Dynamics Simulations of Gel and Fluid Phosphatidylcholine Bilayers

A new united atom parametrization of diacyl lipids like dipalmitoylphosphatidylcholine (DPPC) and the dimyristoylphosphatidylcholine (DMPC) has been constructed based on ab initio calculations to obtain fractional charges and the dihedral potential of the hydrocarbon chains, while the Lennard-Jones parameters of the acyl chains were fitted to reproduce the properties of liquid hydrocarbons. The results have been validated against published experimental X-ray and neutron scattering data for fluid and gel phase DPPC. The derived charges of the lipid phosphatidylcholine (PC) headgroup are shown to yield dipole components in the range suggested by experiments. The aim has been to construct a new force field that retains and improves the good agreement for the fluid phase and at the same time produces a gel phase at low temperatures, with properties coherent with experimental findings. The gel phase of diacyl-PC lipids forms a regular triangular lattice in the hydrocarbon region. The global bilayer tilt obtains an azimuthal value of 31° and is aligned between lattice vectors in the bilayer plane. We also show that the model yields a correct heat of melting as well as decent heat capacities in the fluid and gel phase of DPPC.

Cyclic and Linear Monoterpenes in Phospholipid Membranes: Phase Behavior, Bilayer Structure, and Molecular Dynamics

Monoterpenes are abundant in essential oils extracted from plants. These relatively small and hydrophobic molecules have shown important biological functions, including antimicrobial activity and membrane penetration enhancement. The interaction between the monoterpenes and lipid bilayers is considered important to the understanding of the biological functions of monoterpenes. In this study, we investigated the effect of cyclic and linear monoterpenes on the structure and dynamics of lipids in model membranes. We have studied the ternary system 1,2-dimyristoyl-sn-glycero-3-phosphocholine-monoterpene-water as a model with a focus on dehydrated conditions. By combining complementary techniques, including differential scanning calorimetry, solid-state nuclear magnetic resonance, and small- and wide-angle X-ray scattering, bilayer structure, phase transitions, and lipid molecular dynamics were investigated at different water contents. Monoterpenes cause pronounced melting point depression and phase segregation in lipid bilayers, and the extent of these effects depends on the hydration conditions. The addition of a small amount of thymol to the fluid bilayer (volume fraction of 0.03 in the bilayer) leads to an increased order in the acyl chain close to the bilayer interface. The findings are discussed in relation to biological systems and lipid formulations.

Control of membrane permeability in air-stable droplet interface bilayers

Air-stable droplet interface bilayers (airDIBs) on oil-infused surfaces are versatile model membranes for synthetic biology applications, including biosensing of airborne species. However, airDIBs are subject to evaporation, which can, over time, destabilize them and reduce their useful lifetime compared to traditional DIBs that are fully submerged in oil. Here, we show that the lifetimes of airDIBs can be extended by as much as an order of magnitude by maintaining the temperature just above the dew point. We find that raising the temperature from near the dew point (which was 7 °C at 38.5% relative humidity and 22 °C air temperature) to 20 °C results in the loss of hydrated water molecules from the polar headgroups of the lipid bilayer membrane due to evaporation, resulting in a phase transition with increased disorder. This dehydration transition primarily affects the bilayer electrical resistance by increasing the permeability through an increasingly disordered polar headgroup region of the bilayer. Temperature and relative humidity are conveniently tunable parameters for controlling the stability and composition of airDIB membranes while still allowing for operation in ambient environments.

Structure of the DMPC lipid bilayer ripple phase

High resolution structure is presented for the ripple (Pβ') phase of the phospholipid dimyristoylphosphatidylcholine. Low angle X-ray scattering from oriented samples yielded 57 orders, more than twice as many as recorded previously. The determined electron density map has a sawtooth profile similar to the result from lower resolution data, but the features are sharper allowing better estimates for the modulated bilayer profile and the distribution of headgroups along the aqueous interface. Analysis of high resolution wide angle X-ray data shows that the hydrocarbon chains in the longer, major side of the asymmetric sawtooth are packed similarly to the LβF gel phase, with chains in both monolayers coupled and tilted by 18° in the same direction. The absence of Bragg rods that could be associated with the minor side is consistent with disordered chains, as often suggested in the literature. However, the new high resolution bilayer profile strongly suggests that the chains in the two monolayers in the minor side and the curved region are not in registry. This staggered monolayer modulated melting suggests a direction for improving theories of the ripple phase.

Mesoscale computational studies of membrane bilayer remodeling by curvature-inducing proteins

Biological membranes constitute boundaries of cells and cell organelles. These membranes are soft fluid interfaces whose thermodynamic states are dictated by bending moduli, induced curvature fields, and thermal fluctuations. Recently, there has been a flood of experimental evidence highlighting active roles for these structures in many cellular processes ranging from trafficking of cargo to cell motility. It is believed that the local membrane curvature, which is continuously altered due to its interactions with myriad proteins and other macromolecules attached to its surface, holds the key to the emergent functionality in these cellular processes. Mechanisms at the atomic scale are dictated by protein-lipid interaction strength, lipid composition, lipid distribution in the vicinity of the protein, shape and amino acid composition of the protein, and its amino acid contents. The specificity of molecular interactions together with the cooperativity of multiple proteins induce and stabilize complex membrane shapes at the mesoscale. These shapes span a wide spectrum ranging from the spherical plasma membrane to the complex cisternae of the Golgi apparatus. Mapping the relation between the protein-induced deformations at the molecular scale and the resulting mesoscale morphologies is key to bridging cellular experiments across the various length scales. In this review, we focus on the theoretical and computational methods used to understand the phenomenology underlying protein-driven membrane remodeling. Interactions at the molecular scale can be computationally probed by all atom and coarse grained molecular dynamics (MD, CGMD), as well as dissipative particle dynamics (DPD) simulations, which we only describe in passing. We choose to focus on several continuum approaches extending the Canham - Helfrich elastic energy model for membranes to include the effect of curvature-inducing proteins and explore the conformational phase space of such systems. In this description, the protein is expressed in the form of a spontaneous curvature field. The approaches include field theoretical methods limited to the small deformation regime, triangulated surfaces and particle-based computational models to investigate the large-deformation regimes observed in the natural state of many biological membranes. Applications of these methods to understand the properties of biological membranes in homogeneous and inhomogeneous environments of proteins, whose underlying curvature fields are either isotropic or anisotropic, are discussed. The diversity in the curvature fields elicits a rich variety of morphological states, including tubes, discs, branched tubes, and caveola. Mapping the thermodynamic stability of these states as a function of tuning parameters such as concentration and strength of curvature induction of the proteins is discussed. The relative stabilities of these self-organized shapes are examined through free-energy calculations. The suite of methods discussed here can be tailored to applications in specific cellular settings such as endocytosis during cargo trafficking and tubulation of filopodial structures in migrating cells, which makes these methods a powerful complement to experimental studies.

Stratum corneum molecular mobility in the presence of natural moisturizers

The outermost layer of the skin, the stratum corneum (SC), is a lipid-protein membrane that experiences considerable osmotic stress from a dry and cold climate. The natural moisturizing factor (NMF) comprises small and polar substances, which like osmolytes can protect living systems from osmotic stress. NMF is commonly claimed to increase the water content in the SC and thereby protect the skin from dryness. In this work we challenge this proposed mechanism, and explore the influence of NMF on the lipid and protein components in the SC. We employ natural-abundance (13)C solid-state NMR methods to investigate how the SC molecular components are influenced by urea, glycerol, pyrrolidone carboxylic acid (PCA), and urocanic acid (UCA), all of which are naturally present in the SC as NMF compounds. Experiments are performed with intact SC, isolated corneocytes and model lipids. The combination of NMR experiments provides molecularly resolved qualitative information on the dynamics of different SC lipid and protein components. We obtain completely novel molecular information on the interaction of these NMF compounds with the SC lipids and proteins. We show that urea and glycerol, which are also common ingredients in skin care products, increase the molecular mobility of both SC lipids and proteins at moderate relative humidity where the SC components are considerably more rigid in the absence of these compounds. This effect cannot be attributed to increased SC water content. PCA has no detectable effect on SC molecular mobility under the conditions investigated. It is finally shown that the more apolar compound, UCA, specifically influences the mobility of the SC lipid regions. The present results show that the NMF components act to retain the fluidity of the SC molecular components under dehydrating conditions in such a way that the SC properties remain largely unchanged as compared to more hydrated SC. These findings provide a new molecular insight into how small polar molecules in NMF and skin care products act to protect the human skin from drying.

Physics of Free-Standing Lyotropic Films

We explore the structures and properties of stable, free-standing films of lyotropic mesophases drawn on apertures of various shapes in an atmosphere of controlled humidity. New phenomena are uncovered and interpreted.

Comparing lipid membranes in different environments

When engineering lipid membranes for applications, it is essential to characterize them to avoid artifacts introduced by manipulation and the experimental environment. Wide-angle X-ray scattering is a powerful structural characterization tool for well-ordered lipid systems. It reveals remarkable differences in rotational order parameters for samples prepared in different ways. New data and perspectives are presented here for multilamellar systems that support and extend the characterization work on unilamellar systems that is reported by Watkins et al. in this issue of ACS Nano.

Competition between anisometric and aliphatic entities: an unusual phase sequence with the induction of a phase in an n-alkane-liquid crystal binary system

In this work, we demonstrate two important features that arise out of introducing a liquid-crystalline (LC) compound into the rotator phase matrix and the consequent competition between the anisometric segments of the LC moieties and the aliphatic units. First, we show that the change in the structural character of the mixed medium depends on which of the entities forms the minority concentration: in the case of this being the alkane, the two components of the binary system are nanophase segregated, whereas if the LC molecules are present in a small concentration, then the layered structure merely gets roughened without any segregation. The second and more significant result of the calorimetric and X-ray experiments, at low LC concentrations, is the induction of a rotator phase that leads to unusual phase sequence not reported hitherto. Possible scenarios for the molecular arrangement are discussed. A Landau model is also presented that explains some of the observed features.

Structure and conformational dynamics of DMPC/dicationic surfactant and DMPC/dicationic surfactant/DNA systems

Amphiphilic dicationic surfactants, known as gemini surfactants, are currently studied for gene delivery purposes. The gemini surfactant molecule is composed of two hydrophilic “head” groups attached to hydrophobic chains and connected via molecular linker between them. The influence of different concentrations of 1,5-bis (1-imidazolilo-3- decyloxymethyl) pentane chloride (gemini surfactant) on the thermotropic phase behaviour of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers with and without the presence of DNA was investigated using Fourier transformed infrared (FTIR) and circular dichroism (CD) spectroscopies, small angle scattering of synchrotron radiation and differential scanning calorimetry. With increasing concentration of surfactant in DMPC/DNA systems, a disappearance of pretransition and a decrease in the main phase transition enthalpy and temperature were observed. The increasing intensity of diffraction peaks as a function of surfactant concentration also clearly shows the ability of the surfactant to promote the organisation of lipid bilayers in the multilayer lamellar phase.

Structural evolution of environmentally responsive cationic liposome-DNA complexes with a reducible lipid linker

Environmentally responsive materials (i.e., materials that respond to changes in their environment with a change in their properties or structure) are attracting increasing amounts of interest. We recently designed and synthesized a series of cleavable multivalent lipids (CMVLn, with n = 2-5 being the number of positive headgroup charges at full protonation) with a disulfide bond in the linker between their cationic headgroup and hydrophobic tails. The self-assembled complexes of the CMVLs and DNA are a prototypical environmentally responsive material, undergoing extensive structural rearrangement when exposed to reducing agents. We investigated the structural evolution of CMVL-DNA complexes at varied complex composition, temperature, and incubation time using small-angle X-ray scattering (SAXS) and wide-angle X-ray scattering (WAXS). A related lipid with a stable linker, TMVL4, was used as a control. In a nonreducing environment, CMVL-DNA complexes form the lamellar (L(α)(C)) phase, with DNA rods sandwiched between lipid bilayers. However, new self-assembled phases form when the disulfide linker is cleaved by dithiothreitol or the biologically relevant reducing agent glutathione. The released DNA and cleaved CMVL headgroups form a loosely organized phase, giving rise to a characteristic broad SAXS correlation profile. CMVLs with high headgroup charge also form condensed DNA bundles. Intriguingly, the cleaved hydrophobic tails of the CMVLs reassemble into tilted chain-ordered L(β') phases upon incubation at physiological temperature (37 °C), as indicated by characteristic WAXS peaks. X-ray scattering further reveals that two of the three phases (L(βF), L(βL), and L(βI)) constituting the L(β') phase coexist in these samples. The described system may have applications in lipid-based nanotechnologies.

Modulated phases of phospholipid bilayers induced by tocopherols

The influence of α-, γ- and δ-tocopherols on the structure and phase behavior of dipalmitoyl phosphatidylcholine (DPPC) bilayers has been determined from X-ray diffraction studies on oriented multilayers. In all the three cases the main-transition temperature (T(m)) of DPPC was found to decrease with increasing tocopherol concentration up to around 25 mol%. Beyond this the main transition is suppressed in the case of γ-tocopherol, whereas T(m) becomes insensitive to composition in the other two cases. The pre-transition is found to be suppressed over a narrow tocopherol concentration range between 7.5 and 10 mol% in DPPC-γ-tocopherol and DPPC-δ-tocopherol bilayers, and the ripple phase occurs down to the lowest temperature studied. In all the three cases a modulated phase is observed above a tocopherol concentration of about 10 mol%, which is similar to the P(β) phase reported in DPPC-cholesterol bilayers. This phase is found to occur even in excess water conditions at lower tocopherol concentrations, and consists of bilayers with periodic height modulation. These results indicate the ability of tocopherols to induce local curvature in membranes, which could be important for some of their biological functions.

Phase behavior of two-component lipid membranes: theory and experiments

The structure of the ripple phase of phospholipid membranes remains poorly understood in spite of a large number of theoretical studies, with many experimentally established structural features of this phase unaccounted for. In this article we present a phenomenological theory of phase transitions in single- and two-component achiral lipid membranes in terms of two coupled order parameters: a scalar order parameter describing lipid chain melting, and a vector order parameter describing the tilt of the hydrocarbon chains below the chain-melting transition. This model reproduces all the salient structural features of the ripple phase, providing a unified description of the phase diagram and microstructure. In addition, it predicts a variant of this phase that does not seem to have been experimentally observed. Using this model we have calculated generic phase diagrams of two-component membranes. We have also determined the phase diagram of a two-component lipid membrane from x-ray diffraction studies on aligned multilayers. This phase diagram is found to be in good agreement with that calculated from the model.

Environmental chamber for in situ dynamic control of temperature and relative humidity during x-ray scattering

We have designed, constructed, and evaluated an environmental chamber that has in situ dynamic control of temperature (25 to 90 °C) and relative humidity (0% to 95%). The compact specimen chamber is designed for x-ray scattering in transmission with an escape angle of 2θ = ±30°. The specimen chamber is compatible with a completely evacuated system such as the Rigaku PSAXS system, in which the specimen chamber is placed inside a larger evacuated chamber (flight path). It is also compatible with x-ray systems consisting of evacuated flight tubes separated by small air gaps for sample placement. When attached to a linear motor (vertical displacement), the environmental chamber can access multiple sample positions. The temperature and relative humidity inside the specimen chamber are controlled by passing a mixture of dry and saturated gas through the chamber and by heating the chamber walls. Alternatively, the chamber can be used to control the gaseous environment without humidity. To illustrate the value of this apparatus, we have probed morphology transformations in Nafion(®) membranes and a polymerized ionic liquid as a function of relative humidity in nitrogen.

Hydration effects on membrane structure probed by single molecule orientations

Single molecule fluorescence measurements are used to probe the structural changes in glass-supported DPPC bilayers as a function of relative humidity (RH). Defocused polarized total internal reflection fluorescence microscopy is employed to determine the three-dimensional orientation of the fluorescent lipid analogue BODIPY-PC, doped into DPPC membranes in trace amounts. Supported DPPC bilayers formed using vesicle fusion and Langmuir-Blodgett/Langmuir-Schäfer (LB/LS) transfer are compared and show similar trends as a function of relative humidity. Population histograms of the emission dipole tilt angle reveal bimodal distributions as observed previously for BODIPY-PC in DPPC. These distributions are dominated by large populations of BODIPY-PC molecules with emission dipoles oriented parallel (≥81°) and normal (≤10°) to the membrane plane, with less than 25% oriented at intermediate tilts. As the relative humidity is increased from 13% to 95%, the population of molecules oriented normal to the surface decreases with a concomitant increase in those oriented parallel to the surface. The close agreement in trends observed for bilayers formed from vesicle fusion and LB/LS transfer supports the assignment of an equivalent surface pressure of 23 mN/m for bilayers formed from vesicle fusion. At each RH condition, a small population of BODIPY-PC dye molecules are laterally mobile in both bilayer preparations. This population exponentially increases with RH but never exceeds 6% of the total population. Interestingly, even under conditions where there is little lateral diffusion, fluctuations in the single molecule orientations can be observed which suggests there is appreciable freedom in the acyl chain region. Dynamic measurements of single molecule orientation changes, therefore, provide a new view into membrane properties at the single molecule level.

X-Ray Diffraction Studies of Tubules Formed from a Diacetylenic Phosphocholine Lipid

The diacetylenic phosphocholine lipid 1,2-bis(10, 12-tricosadiyonol)-sn-glycerol- 3-phosphocholine forms, upon cooling below the chain melting temperature, hollow cylindrical structures known as ‘tubules’. These tubules are approximately one micron in diameter and may range from tens to hundreds of microns in length. We have carried out x-ray diffraction studies of tubules (1) in 75/25 ethanol/water by volume and (2) in a dessicated state. We present results which show marked differences in the nature of both the in-plane and inter-layer correlations associated with the wrapped multilayer structure of the tubules in the two states.

Dynamic structure factors from lipid membrane molecular dynamics simulations

Dynamic structure factors for a lipid bilayer have been calculated from molecular dynamics simulations. From trajectories of a system containing 1024 lipids we obtain wave vectors down to 0.34 nm(-1), which enables us to directly resolve the Rayleigh and Brillouin lines of the spectrum. The results confirm the validity of a model based on generalized hydrodynamics, but also improves the line widths and the position of the Brillouin lines. The improved resolution shows that the Rayleigh line is narrower than in earlier studies, which corresponds to a smaller thermal diffusivity. From a detailed analysis of the power spectrum, we can, in fact, distinguish two dispersive contributions to the elastic scattering. These translate to two exponential relaxation processes in separate time domains. Further, by including a first correction to the wave-vector-dependent position of the Brillouin lines, the results agree favorably to generalized hydrodynamics even up to intermediate wave vectors, and also yields a 20% higher adiabatic sound velocity. The width of the Brillouin lines shows a linear, not quadratic, dependence to low wave vectors.

Effects of ether vs. ester linkage on lipid bilayer structure and water permeability

The structure and water permeability of bilayers composed of the ether-linked lipid, dihexadecylphosphatidylcholine (DHPC), were studied and compared with the ester-linked lipid, dipalmitoylphosphaditdylcholine (DPPC). Wide angle X-ray scattering on oriented bilayers in the fluid phase indicate that the area per lipid A is slightly larger for DHPC than for DPPC. Low angle X-ray scattering yields A=65.1A(2) for DHPC at 48 degrees C. LAXS data provide the bending modulus, K(C)=4.2x10(-13)erg, and the Hamaker parameter H=7.2x10(-14)erg for the van der Waals attractive interaction between neighboring bilayers. For the low temperature phases with ordered hydrocarbon chains, we confirm the transition from a tilted L(beta') gel phase to an untilted, interdigitated L(beta)I phase as the sample hydrates at 20 degrees C. Our measurement of water permeability, P(f)=0.022cm/s at 48 degrees C for fluid phase DHPC is slightly smaller than that of DPPC (P(f)=0.027cm/s) at 50 degrees C, consistent with our triple slab theory of permeability.

Atomistic simulations of mixed-lipid bilayers in gel and fluid phases

The slow rate of diffusive mixing poses a challenge for molecular dynamics (MD) simulation studies of mixed-lipid bilayers. A mixed Monte Carlo-molecular dynamics (MC-MD) approach, which uses mutation moves to swap lipid types throughout the system within the semi-grand canonical ensemble, is here applied to a comparison of binary mixtures in the gel and liquid crystalline phases. The two lipid components modeled, distearoylphosphatidylcholine (DSPC) and dimyristoylphosphatidylcholine (DMPC), differ by four carbons in the lengths of their acyl tails and are investigated here at full hydration at a temperature (313 K) between their transition temperatures, where coexistence between a DSPC-rich gel phase and a DMPC-rich liquid crystalline phase is expected. An analysis of DSPC-DMPC mixtures in the gel phase indicates strong deviation from ideality in the thermodynamics of mixing, accompanied by a tendency of the shorter-tailed component DMPC to associate laterally and for DMPC headgroups to be displaced toward the bilayer midplane. The liquid crystal phase mixtures, in contrast, show more mild deviation from thermodynamically ideal mixing with no apparent tendency for similar lipids to cluster laterally and no difference in headgroup normal distribution profiles.

Defect-mediated emulsification in two dimensions

We consider two-dimensional dispersions of droplets of isotropic phase in a liquid with an XY -like order parameter, tilt, nematic, and hexatic symmetries being included. Strong anchoring boundary conditions are assumed. Textures for a single droplet and a pair of droplets are calculated and a universal droplet-droplet pair potential is obtained. The interaction of dispersed droplets via the ordered phase is attractive at large distances and repulsive at short distances, which results in a well defined preferred separation for two droplets and topological stabilization of the emulsion. This interaction also drives self-assembly into chains. Preferred separations and energy barriers to coalescence are calculated, and the effects of thermal fluctuations and film thickness are discussed.

Electrostatic interactions in phospholipid membranes revealed by coherent 2D IR spectroscopy

The inter- and intramolecular interactions of the carbonyl moieties at the polar interface of a phospholipid membrane are probed by using nonlinear femtosecond infrared spectroscopy. Two-dimensional IR correlation spectra separate homogeneous and inhomogeneous broadenings and show a distinct cross-peak pattern controlled by electrostatic interactions. The inter- and intramolecular electrostatic interactions determine the inhomogeneous character of the optical response. Using molecular dynamics simulation and the nonlinear exciton equations approach, we extract from the spectra short-range structural correlations between carbonyls at the interface.

Vesicle shape, molecular tilt, and the suppression of necks

Can the presence of molecular-tilt order significantly affect the shapes of lipid bilayer membranes, particularly membrane shapes with narrow necks? Motivated by the propensity for tilt order and the common occurrence of narrow necks in the intermediate stages of biological processes such as endocytosis and vesicle trafficking, we examine how tilt order inhibits the formation of necks in the equilibrium shapes of vesicles. For vesicles with a spherical topology, point defects in the molecular order with a total strength of +2 are required. We study axisymmetric shapes and suppose that there is a unit-strength defect at each pole of the vesicle. The model is further simplified by the assumption of tilt isotropy: invariance of the energy with respect to rotations of the molecules about the local membrane normal. This isotropy condition leads to a minimal coupling of tilt order and curvature, giving a high energetic cost to regions with Gaussian curvature and tilt order. Minimizing the elastic free energy with constraints of fixed area and fixed enclosed volume determines the allowed shapes. Using numerical calculations, we find several branches of solutions and identify them with the branches previously known for fluid membranes. We find that tilt order changes the relative energy of the branches, suppressing thin necks by making them costly, leading to elongated prolate vesicles as a generic family of tilt-ordered membrane shapes.

Structure and Phase Behavior of Self-Assembled DPPC−DNA−Metal Cation Complexes

Multilamellar liposomes of dipalmitoylphosphatidylcholine (DPPC) in solution with DNA and bivalent metal cations (Ca2+, Mn2+, Mg2+) self-assemble into a ternary DPPC-DNA-Me2+ complex. The supramolecular structure of the complex consists of an ordered multilamellar assembly where hydrated DNA helices are sandwiched between the lipid bilayers and the metal cations bind the phosphate groups of DNA to the lipid polar heads. In the range of explored incubation times, the complex coexists with the uncomplexed DPPC over the whole temperature range investigated (20-55 degrees C). Accordingly, two distinct coexisting lamellar phases are observed, one corresponding to the ternary complex and the other to the uncomplexed lipid. The structure and thermotropic phase behavior of both of these have been investigated by means of synchrotron X-ray diffraction, and the relevant structural data are deduced from experimental electron density profiles. While the uncomplexed lipid exhibits the same phase behavior as pure DPPC, that is, L beta'-P beta'-L alpha, the thermotropic behavior of the bound lipid in the complex is partially altered. This is manifested as an increase in the main transition temperature and the disappearance of the ripple phase leading to the single -phase transition. The role of the different metal cations in promoting and stabilizing the DNA condensation into the ternary complex is also discussed.

Solvent-free simulations of fluid membrane bilayers

A molecular level model for lipid bilayers is presented. Lipids are represented by rigid, asymmetric, soft spherocylinders in implicit solvent. A simple three parameter potential between pairs of lipids gives rise to a rich assortment of phases including (but not limited to) micelles, fluid bilayers, and gel-like bilayers. Monte Carlo simulations have been carried out to verify self-assembly, characterize the phases corresponding to different potential parametrizations, and to quantify the physical properties associated with those parameter sets corresponding to fluid bilayer behavior. The studied fluid bilayers have compressibility moduli in agreement with experimental systems, but display bending moduli at least three times larger than typical biological membranes without cholesterol.

Implicit solvent simulation models for biomembranes

Fully atomic simulation strategies are infeasible for the study of many processes of interest to membrane biology, biophysics and biochemistry. We review various coarse-grained simulation methodologies with special emphasis on methods and models that do not require the explicit simulation of water. Examples from our own research demonstrate that such models have potential for simulating a variety of biologically relevant phenomena at the membrane surface.