Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature

Molecular packing of amphipathic peptides on the surface of lipid membranes

When polypeptides bind to the membrane surface, they become confined to a restricted quasi-two-dimensional space where peptide-peptide interactions become highly relevant, and the concept of a crowded medium is appropriate. Within this crowded environment interesting effects like clustering, separation of phases, cooperative alignment, and common movements occur. Here we investigated such effects by measuring distances between fluorophore-labeled peptides in the range ≤1 nm by fluorescence self-quenching. For helical peptides with dimensions of approximately 1 × 3 nm such a small "ruler" is sensitive to the packing of the labeled peptides and thereby to their molecular arrangement. A novel approach to characterize peptide-peptide interactions within membranes is presented using the designer peptide LAH4. This sequence changes membrane topology in a controlled manner being transmembrane at neutral conditions but oriented parallel to the surface at low pH. Experimental measurements of the fluorescence self-quenching of close-by chromophores and the changes that occur upon dilution with unlabeled peptides are used to analyze the peptide distribution within the membrane surface. The data show a strong effect of electrostatic interactions and under some experimental conditions clustering of the peptides. Furthermore, the results suggest that at pH 4 the peptides arrange along the membrane surface in an ordered mesophase-like arrangement.

Phase segregation in bio-inspired multi-component vesicles encompassing double tail phospholipid species

Our aim is to investigate the phase segregation and the structure of multi-component bio-inspired phospholipid vesicles via dissipative particle dynamics. The chemical distinction in the phospholipid species arises due to different head and tail group moieties, and molecular stiffness of the hydrocarbon tails. The individual amphiphilic phospholipid molecular species are represented by a hydrophilic head group and two hydrophobic tails. The distinct chemical nature of the moieties is modeled effectively via soft repulsive interaction parameters, and the molecular rigidity is tuned via suitable three-body potential constants. We demonstrate the formation of a stable hybrid vesicle through the self-assembly of the amphiphilic phospholipid molecules in the presence of a hydrophilic solvent. We investigate and characterize the phase segregation and the structure of the binary vesicles for different phospholipid mixtures. Our results demonstrate macroscopic phase separation for phospholipid mixtures composed of species with different hydrocarbon tail groups. We also investigate the relationship between the phase segregation and thermodynamic variables such as interfacial line tension and surface tension, and obtain correspondence between existing theory and experiments, and our simulation results. We report variations in the molecular chain stiffness to have negligible contributions to the phase segregation in the mixed bilayer, and to demonstrate shape transformations of the hybrid vesicle. Our results can be used to design novel bio-inspired hybrid vehicles for potential applications in biomedicine, sensing, imaging and sustainability.

Highly hydrated deformable polyethylene glycol-tethered lipid bilayers

The realization of a solid-supported lipid bilayer acting as a workbench for the study of membrane processes is a difficult task. For robustness, the bilayer has to be tethered to the substrate. At the same time, diffusion of the lipids and plastic deformations of the membrane should not be obstructed. Furthermore, a highly hydrated surrounding is mandatory. Here, we show that grafting of a polyethylene glycol-lipid construct (PEG2000-DSPE) to a silicon oxide surface via multiple-step silane chemistry and subsequent deposition of lipids by spin-coating result in a cushioned membrane that has the desired properties. Neutron and X-ray reflectometry measurements are combined to access thickness, density, and hydration of the bilayer and the PEG cushion. We observe a spacer of 55 Å thickness between lipid bilayer and silicon-oxide surface with a rather high hydration of up to 90 ± 3% water. While 11.5 ± 3% of the lipids are grafted to the surface, as determined from the neutron data, the diffusion constant of the lipids, as probed by diffusion of 0.5% Texas Red labeled lipids, remains rather large (D = 2.1 ± 0.1 μm(2)/s), which is a reduction of only 12% compared to a supported lipid bilayer reference without immobilized lipids. Finally, AFM indentation confirms the plastic behavior of the membrane against deformation. We show that rupture of the bilayer does not occur before the deformation exceeds 40 Å. Altogether, the presented PEG-tethered lipid bilayer mimics the deformability of natural cell membranes much better than standard solid-supported lipid bilayers.

Molecular simulation of ibuprofen passing across POPC membrane

Permeability assessment is an important procedure in the drug development process, and drug partitioning in membrane bilayer is related to permeability. To investigate the pH dependence on drug partitioning, the process of different ionization state of ibuprofen passing across POPC bilayer was studied using molecular dynamics simulation. The results show that both atomic charge scheme and ionization state of the drug affect the value and shape of energy profile when passing across the POPC bilayer. The neutral ibuprofen (ibuprofen under acidic condition) has a much lower energy barrier as compared with the anionic ibuprofen (ibuprofen under basic condition). Meantime, hydrogen bond analysis also certifies that it is easy for neutral ibuprofen to pass from bulk water to bilayer center. Our calculation suggests that the ionization state of ibuprofen may be changed between neutral and anionic state when passing across membrane: it may be ionized outside the membrane and neutralized inside the membrane.

The Bio-SANS instrument at the High Flux Isotope Reactor of Oak Ridge National Laboratory

Small-angle neutron scattering (SANS) is a powerful tool for characterizing complex disordered materials, including biological materials. The Bio-SANS instrument of the High Flux Isotope Reactor of Oak Ridge National Laboratory (ORNL) is a high-flux low-background SANS instrument that is, uniquely among SANS instruments, dedicated to serving the needs of the structural biology and biomaterials communities as an open-access user facility. Here, the technical specifications and performance of the Bio-SANS are presented. Sample environments developed to address the needs of the user program of the instrument are also presented. Further, the isotopic labeling and sample preparation capabilities available in the Bio-Deuteration Laboratory for users of the Bio-SANS and other neutron scattering instruments at ORNL are described. Finally, a brief survey of research performed using the Bio-SANS is presented, which demonstrates the breadth of the research that the instrument's user community engages in.

A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties

Although lipid force fields (FFs) used in molecular dynamics (MD) simulations have proved to be accurate, there has not been a systematic study on their accuracy over a range of temperatures. Motivated by the X-ray and neutron scattering measurements of common phosphatidylcholine (PC) bilayers (Kučerka et al. BBA. 1808: 2761, 2011), the CHARMM36 (C36) FF accuracy is tested in this work with MD simulations of six common PC lipid bilayers over a wide range of temperatures. The calculated scattering form factors and deuterium order parameters from the C36 MD simulations agree well with the X-ray, neutron, and NMR experimental data. There is excellent agreement between MD simulations and experimental estimates for the surface area per lipid, bilayer thickness (DB), hydrophobic thickness (DC), and lipid volume (VL). The only minor discrepancy between simulation and experiment is a measure of (DB-DHH)/2 where DHH is the distance between the maxima in the electron density profile along the bilayer normal. Additional MD simulations with pure water and heptane over a range of temperatures provide explanations of possible reasons causing the minor deviation. Overall, the C36 FF is accurate for use with liquid crystalline PC bilayers of varying chain types and over biologically relevant temperatures.

The molecular structure of a phosphatidylserine bilayer determined by scattering and molecular dynamics simulations

Phosphatidylserine (PS) lipids play essential roles in biological processes, including enzyme activation and apoptosis. We report on the molecular structure and atomic scale interactions of a fluid bilayer composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS). A scattering density profile model, aided by molecular dynamics (MD) simulations, was developed to jointly refine different contrast small-angle neutron and X-ray scattering data, which yielded a lipid area of 62.7 Å(2) at 25 °C. MD simulations with POPS lipid area constrained at different values were also performed using all-atom and aliphatic united-atom models. The optimal simulated bilayer was obtained using a model-free comparison approach. Examination of the simulated bilayer, which agrees best with the experimental scattering data, reveals a preferential interaction between Na(+) ions and the terminal serine and phosphate moieties. Long-range inter-lipid interactions were identified, primarily between the positively charged ammonium, and the negatively charged carboxylic and phosphate oxygens. The area compressibility modulus KA of the POPS bilayer was derived by quantifying lipid area as a function of surface tension from area-constrained MD simulations. It was found that POPS bilayers possess a much larger KA than that of neutral phosphatidylcholine lipid bilayers. We propose that the unique molecular features of POPS bilayers may play an important role in certain physiological functions.

Tris-thiourea tripodal-based molecules as chloride transmembrane transporters: insights from molecular dynamics simulations

The interaction of six tripodal synthetic chloride transmembrane transporters with a POPC bilayer was investigated by means of molecular dynamics simulations using the general Amber force field (GAFF) for the transporters and the LIPID11 force field for phospholipids. These transporters are structurally simple molecules, based on the tris(2-aminoethyl)amine scaffold, containing three thiourea binding units coupled with three n-butyl (1), phenyl (2), fluorophenyl (3), pentafluorophenyl (4), trifluoromethylphenyl (5), or bis(trifluoromethyl)phenyl (6) substituents. The passive diffusion of 1-6⊃ Cl(-) was evaluated with the complexes initially positioned either in the water phase or inside the bilayer. In the first scenario the chloride is released in the water solution before the synthetic molecules achieve the water-lipid interface and permeate the membrane. In the latter one, only when the chloride complex reaches the interface is the anion released to the water phase, with the transporter losing the initial ggg tripodal shape. Independently of the transporter used in the membrane system, the bilayer structure is preserved and the synthetic molecules interact with the POPC molecules at the phosphate headgroup level, via N-H···O hydrogen bonds. Overall, the molecular dynamics simulations' results indicate that the small tripodal molecules in this series have a low impact on the bilayer and are able to diffuse with chloride inside the lipid environment. Indeed, these are essential conditions for these molecules to promote the transmembrane transport as anion carriers, in agreement with experimental efflux data.

Molecular dynamics simulations of homo-oligomeric bundles embedded within a lipid bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)3-NH2 (LS2) or Ac-(LSSLLSL)3-NH2 (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine-serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a "dimer of dimers." In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.

Do lipids show state-dependent affinity to the voltage-gated potassium channel KvAP?

As all integral membrane proteins, voltage-gated ion channels are embedded in a lipid matrix that regulates their channel behavior either by physicochemical properties or by direct binding. Because manipulation of the lipid composition in cells is difficult, we investigated the influence of different lipids on purified KvAP channels reconstituted in planar lipid bilayers of known composition. Lipids developed two distinct and independent effects on the KvAP channels; lipids interacting with the pore lowered the energy barriers for the final transitions, whereas voltage sensor-bound lipids shifted the midpoint of activation dependent on their electrostatic charge. Above all, the midpoint of activation was determined only by those lipids the channels came in contact with first after purification and can seemingly only be exchanged if the channel resides in the open state. The high affinity of the bound lipids to the binding site has implications not only on our understanding of the gating mechanism but also on the general experimental design of any lipid dependence study.

The mechanism of membrane disruption by cytotoxic amyloid oligomers formed by prion protein(106-126) is dependent on bilayer composition

The formation of fibrillar aggregates has long been associated with neurodegenerative disorders such as Alzheimer and Parkinson diseases. Although fibrils are still considered important to the pathology of these disorders, it is now widely understood that smaller amyloid oligomers are the toxic entities along the misfolding pathway. One characteristic shared by the majority of amyloid oligomers is the ability to disrupt membranes, a commonality proposed to be responsible for their toxicity, although the mechanisms linking this to cell death are poorly understood. Here, we describe the physical basis for the cytotoxicity of oligomers formed by the prion protein (PrP)-derived amyloid peptide PrP(106-126). We show that oligomers of this peptide kill several mammalian cells lines, as well as mouse cerebellar organotypic cultures, and we also show that they exhibit antimicrobial activity. Physical perturbation of model membranes mimicking bacterial or mammalian cells was investigated using atomic force microscopy, polarized total internal reflection fluorescence microscopy, and NMR spectroscopy. Disruption of anionic membranes proceeds through a carpet or detergent model as proposed for other antimicrobial peptides. By contrast, when added to zwitterionic membranes containing cholesterol-rich ordered domains, PrP(106-126) oligomers induce a loss of domain separation and decreased membrane disorder. Loss of raft-like domains may lead to activation of apoptotic pathways, resulting in cell death. This work sheds new light on the physical mechanisms of amyloid cytotoxicity and is the first to clearly show membrane type-specific modes of action for a cytotoxic peptide.

Self-assembling peptides form nanodiscs that stabilize membrane proteins

New methods to handle membrane bound proteins, e.g. G-protein coupled receptors (GPCRs), are highly desirable. Recently, apoliprotein A1 (ApoA1) based lipoprotein particles have emerged as a new platform for studying membrane proteins, and it has been shown that they can self-assemble in combination with phospholipids to form discoidal shaped particles that can stabilize membrane proteins. In the present study, we have investigated an ApoA1 mimetic peptide with respect to its solution structure when in complex with phospholipids. This was achieved using a powerful combination of small-angle X-ray scattering (SAXS) and small-angle neutron scattering (SANS) supported by coarse-grained molecular dynamics simulations. The detailed structure of the discs was determined in unprecedented detail and it was found that they adopt a discoidal structure very similar to the ApoA1 based nanodiscs. We furthermore show that, like the ApoA1 and derived nanodiscs, these peptide discs can accommodate and stabilize a membrane protein. Finally, we exploit their dynamic properties and show that the 18A discs may be used for transferring membrane proteins and associated phospholipids directly and gently into phospholipid nanodiscs.

Global small-angle X-ray scattering data analysis for multilamellar vesicles: the evolution of the scattering density profile model

The highly successful scattering density profile (SDP) model, used to jointly analyze small-angle X-ray and neutron scattering data from unilamellar vesicles, has been adapted for use with data from fully hydrated, liquid crystalline multilamellar vesicles (MLVs). Using a genetic algorithm, this new method is capable of providing high-resolution structural information, as well as determining bilayer elastic bending fluctuations from standalone X-ray data. Structural parameters such as bilayer thickness and area per lipid were determined for a series of saturated and unsaturated lipids, as well as binary mixtures with cholesterol. The results are in good agreement with previously reported SDP data, which used both neutron and X-ray data. The inclusion of deuterated and non-deuterated MLV neutron data in the analysis improved the lipid backbone information but did not improve, within experimental error, the structural data regarding bilayer thickness and area per lipid.

Conformation and lipid interaction of the fusion peptide of the paramyxovirus PIV5 in anionic and negative-curvature membranes from solid-state NMR

Viral fusion proteins catalyze the merger of the virus envelope and the target cell membrane through multiple steps of protein conformational changes. The fusion peptide domain of these proteins is important for membrane fusion, but how it causes membrane curvature and dehydration is still poorly understood. We now use solid-state NMR spectroscopy to investigate the conformation, topology, and lipid and water interactions of the fusion peptide of the PIV5 virus F protein in three lipid membranes, POPC/POPG, DOPC/DOPG, and DOPE. These membranes allow us to investigate the effects of lipid chain disorder, membrane surface charge, and intrinsic negative curvature on the fusion peptide structure. Chemical shifts and spin diffusion data indicate that the PIV5 fusion peptide is inserted into all three membranes but adopts distinct conformations: it is fully α-helical in the POPC/POPG membrane, adopts a mixed strand/helix conformation in the DOPC/DOPG membrane, and is primarily a β-strand in the DOPE membrane. (31)P NMR spectra show that the peptide retains the lamellar structure and hydration of the two anionic membranes. However, it dehydrates the DOPE membrane, destabilizes its inverted hexagonal phase, and creates an isotropic phase that is most likely a cubic phase. The ability of the β-strand conformation of the fusion peptide to generate negative Gaussian curvature and to dehydrate the membrane may be important for the formation of hemifusion intermediates in the membrane fusion pathway.

Lipid14: The Amber Lipid Force Field

The AMBER lipid force field has been updated to create Lipid14, allowing tensionless simulation of a number of lipid types with the AMBER MD package. The modular nature of this force field allows numerous combinations of head and tail groups to create different lipid types, enabling the easy insertion of new lipid species. The Lennard-Jones and torsion parameters of both the head and tail groups have been revised and updated partial charges calculated. The force field has been validated by simulating bilayers of six different lipid types for a total of 0.5 μs each without applying a surface tension; with favorable comparison to experiment for properties such as area per lipid, volume per lipid, bilayer thickness, NMR order parameters, scattering data, and lipid lateral diffusion. As the derivation of this force field is consistent with the AMBER development philosophy, Lipid14 is compatible with the AMBER protein, nucleic acid, carbohydrate, and small molecule force fields.

Structural properties of POPC monolayers under lateral compression: computer simulations analysis

1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a lipid comprising a saturated and an unsaturated acyl chain, belongs to the class of glycerophosphatidylcholines, major lipids in eukaryotic cell membranes. To get insight into the structural properties of this lipid within monolayers as membrane models, we performed molecular dynamics (MD) simulations of POPC monolayers under compression at the air/water interface. MD simulations were carried out at 300 K and at different surface pressures using the all-atom general Amber force field (GAFF). A good agreement was found between the simulated data and experimental isotherms. At surface pressures greater than 15 mN/m, two orientations of the head groups clearly appear: one nearly parallel to the monolayer interface and another one pointing toward the water. On the basis of the analysis of headgroup dihedral angles, we propose that the conformational variations around the bonds connecting the phosphorus atom to the adjacent oxygens are involved in these two orientations of the headgroup. The glycerol group orientation is characterized by a large distribution centered around 50° with respect to the monolayer normal. The acyl chains are predominantly in trans configuration from 7.5 to 43 mN/m surface pressures. Moreover, the calculated order parameter profiles of both chains suggest an independent behavior of the saturated and unsaturated chains that could be correlated with the formation of chain-type clusters observed along the simulated trajectories.

Size-dependent ultrafast structural dynamics inside phospholipid vesicle bilayers measured with 2D IR vibrational echoes

The ultrafast structural dynamics inside the bilayers of dilauroylphosphatidylcholine (DLPC) and dipalmitoylphosphatidylcholine vesicles with 70, 90, and 125 nm diameters were directly measured with 2D IR vibrational echo spectroscopy. The antisymmetric CO stretch of tungsten hexacarbonyl was used as a vibrational probe and provided information on spectral diffusion (structural dynamics) in the alkyl region of the bilayers. Although the CO stretch absorption spectra remain the same, the interior structural dynamics become faster as the size of the vesicles decrease, with the size dependence greater for dipalmitoylphosphatidylcholine than for DLPC. As DLPC vesicles become larger, the interior dynamics approach those of the planar bilayer.

Dimyristoyl phosphatidylcholine: a remarkable exception to α-tocopherol's membrane presence

Using data obtained from different physical techniques (i.e., neutron diffraction, NMR and UV spectroscopy), we present evidence which explains some of the conflicting and inexplicable data found in the literature regarding α-tocopherol's (aToc's) behavior in dimyristoyl phosphatidylcholine (di-14:0PC) bilayers. Without exception, the data point to aToc's active chromanol moiety residing deep in the hydrophobic core of di-14:0PC bilayers, a location that is in stark contrast to aToc's location in other PC bilayers. Our result is a clear example of the importance of lipid species diversity in biological membranes and importantly, it suggests that measurements of aToc's oxidation kinetics, and its associated byproducts observed in di-14:0PC bilayers, should be reexamined, this time taking into account its noncanonical location in this bilayer.

A microscopic multiphase diffusion model of viable epidermis permeability

A microscopic model of passive transverse mass transport of small solutes in the viable epidermal layer of human skin is formulated on the basis of a hexagonal array of cells (i.e., keratinocytes) bounded by 4-nm-thick, anisotropic lipid bilayers and separated by 1-μm layers of extracellular fluid. Gap junctions and tight junctions with adjustable permeabilities are included to modulate the transport of solutes with low membrane permeabilities. Two keratinocyte aspect ratios are considered to represent basal and spinous cells (longer) and granular cells (more flattened). The diffusion problem is solved in a unit cell using a coordinate system conforming to the hexagonal cross section, and an efficient two-dimensional treatment is applied to describe transport in both the cell membranes and intercellular spaces, given their thinness. Results are presented in terms of an effective diffusion coefficient, D¯(epi), and partition coefficient, K¯(epi/w), for a homogenized representation of the microtransport problem. Representative calculations are carried out for three small solutes-water, L-glucose, and hydrocortisone-covering a wide range of membrane permeability. The effective transport parameters and their microscopic interpretation can be employed within the context of existing three-layer models of skin transport to provide more realistic estimates of the epidermal concentrations of topically applied solutes.

Monolayer spontaneous curvature of raft-forming membrane lipids

Monolayer spontaneous curvatures for cholesterol, DOPE, POPE, DOPC, DPPC, DSPC, POPC, SOPC, and egg sphingomyelin were obtained using small-angle X-ray scattering (SAXS) on inverted hexagonal phases (HII). Spontaneous curvatures of bilayer forming lipids were estimated by adding controlled amounts to a HII forming template following previously established protocols. Spontaneous curvatures of both phosphatidylethanolamines and cholesterol were found to be at least a factor of two more negative than those of phosphatidylcholines, whose J0 values are closer to zero. Interestingly, a significant positive J0 value was retrieved for DPPC. We further determined the temperature dependence of the spontaneous curvatures J0(T) in the range from 15 to 55 °C, resulting in a quite narrow distribution of -1 to -3 × 10-3 (nm °C)-1 for most investigated lipids. The data allowed us to estimate the monolayer spontaneous curvatures of ternary lipid mixtures showing liquid ordered/liquid disordered phase coexistence. We report spontaneous curvature phase diagrams for DSPC/DOPC/Chol, DPPC/DOPC/Chol and SM/POPC/Chol and discuss effects on protein insertion and line tension.

Xenon and other volatile anesthetics change domain structure in model lipid raft membranes

Inhalation anesthetics have been in clinical use for over 160 years, but the molecular mechanisms of action continue to be investigated. Direct interactions with ion channels received much attention after it was found that anesthetics do not change the structure of homogeneous model membranes. However, it was recently found that halothane, a prototypical anesthetic, changes domain structure of a binary lipid membrane. The noble gas xenon is an excellent anesthetic and provides a pivotal test of the generality of this finding, extended to ternary lipid raft mixtures. We report that xenon and conventional anesthetics change the domain equilibrium in two canonical ternary lipid raft mixtures. These findings demonstrate a membrane-mediated mechanism whereby inhalation anesthetics can affect the lipid environment of transmembrane proteins.

Solid supported membranes doped with PIP2: influence of ionic strength and pH on bilayer formation and membrane organization

Phosphoinositides and in particular L-α-phosphatidylinositol-4,5-bisphosphate (PIP2) are key lipids controlling many cellular events and serve as receptors for a large number of intracellular proteins. To quantitatively analyze protein-PIP2 interactions in vitro in a time-resolved manner, planar membranes on solid substrates are highly desirable. Here, we describe an optimized protocol to form PIP2 containing planar solid supported membranes on silicon surfaces by vesicle spreading. Supported lipid bilayers (SLBs) were obtained by spreading POPC/PIP2 (92:8) small unilamellar vesicles onto hydrophilic silicon substrates at a low pH of 4.8. These membranes were capable of binding ezrin, resulting in large protein coverage as concluded from reflectometric interference spectroscopy and fluorescence microscopy. As deduced from fluorescence microscopy, only under low pH conditions, a homogeneously appearing distribution of fluorescently labeled PIP2 molecules in the membrane was achieved. Fluorescence recovery after photobleaching experiments revealed that PIP2 is not mobile in the bottom layer of the SLBs, while PIP2 is fully mobile in the top layer with diffusion coefficients of about 3 μm(2)/s. This diffusion coefficient was considerably reduced by a factor of about 3 if ezrin has been bound to PIP2 in the membrane.

Iterative layer-by-layer assembly of polymer-tethered multi-bilayers using maleimide–thiol coupling chemistry

The current study reports on the layer-by-layer assembly of a polymer-tethered lipid multi-bilayer stack using the iterative addition and roll out of giant unilamellar vesicles (GUVs) containing constituents with thiol and maleimide functional groups, respectively. Confocal microscopy and photobleaching experiments confirm stack integrity and stability over time, as well as the lateral fluidity of individual bilayers within the stacks. Complementary wide-field single molecule fluorescence microscopy and atomic force microscopy experiments show that increasing bilayer-substrate distances are associated with changes in lipid lateral mobility and bilayer morphology. Importantly, the described iterative approach can be employed to assemble multi-bilayer stacks with more than two bilayers, thus further reducing the influence of the underlying solid substrate on membrane behavior. Furthermore, the presence of lipopolymers within the multi-bilayer stacks results in fascinating membrane dynamics and organization properties, with interesting parallels to those found in plasma membranes. In that sense, the described multi-bilayer architecture represents an attractive model membrane platform for a variety of different biophysical studies.

Simulations of biased agonists in the β(2) adrenergic receptor with accelerated molecular dynamics

The biased agonism of the G protein-coupled receptors (GPCRs), where in addition to a traditional G protein-signaling pathway a GPCR promotes intracellular signals though β-arrestin, is a novel paradigm in pharmacology. Biochemical and biophysical studies have suggested that a GPCR forms a distinct ensemble of conformations signaling through the G protein and β-arrestin. Here we report on the dynamics of the β2 adrenergic receptor bound to the β-arrestin and G protein-biased agonists and the empty receptor to further characterize the receptor conformational changes caused by biased agonists. We use conventional and accelerated molecular dynamics (aMD) simulations to explore the conformational transitions of the GPCR from the active state to the inactive state. We found that aMD simulations enable monitoring of the transition within the nanosecond time scale while capturing the known microscopic characteristics of the inactive states, such as the ionic lock, the inward position of F6.44, and water clusters. Distinct conformational states are shown to be stabilized by each biased agonist. In particular, in simulations of the receptor with the β-arrestin-biased agonist N-cyclopentylbutanepherine, we observe a different pattern of motions in helix 7 when compared to simulations with the G protein-biased agonist salbutamol that involves perturbations of the network of interactions within the NPxxY motif. Understanding the network of interactions induced by biased ligands and the subsequent receptor conformational shifts will lead to development of more efficient drugs.

Environmental conditions that influence the ability of humic acids to induce permeability in model biomembranes

The interaction of humic acids (HAs) with 1-palmitoyl-2-oleoyl-Sn-glycero-3-phosphocholine (POPC) large unilamellar vesicle (LUV) model biomembrane system was studied by fluorescence spectroscopy. HAs from aquatic and terrestrial (including coal) sources were studied. The effects of HA concentration and temperature over environmentally relevant ranges of 0 to 20 mg C/L and 10 to 30 °C, respectively, were investigated. The dosage studies revealed that the aquatic Suwannee River humic acid (SRHA) causes an increased biomembrane perturbation (percent leakage of the fluorescent dye, Sulforhodamine B) over the entire studied concentration range. The two terrestrial HAs, namely Leonardite humic acid (LAHA) and Florida peat humic acid (FPHA), at concentrations above 5 mg C/L, show a decrease or a plateau effect attributable to the competition within the HA mixture and/or the formation of "partial aggregates". The temperature studies revealed that biomembrane perturbation increases with decreasing temperature for all three HAs. Kinetic studies showed that the membrane perturbation process is complex with both fast and slow absorption (sorption into the bilayer) components and that the slow component could be fitted by first order kinetics. A mechanism based on "lattice errors" within the POPC LUVs is put forward to explain the fast and slow components. A rationale behind the concentration and temperature findings is provided, and the environmental implications are discussed.

The Observation of Highly Ordered Domains in Membranes with Cholesterol

Rafts, or functional domains, are transient nano- or mesoscopic structures in the exoplasmic leaflet of the plasma membrane, and are thought to be essential for many cellular processes. Using neutron diffraction and computer modelling, we present evidence for the existence of highly ordered lipid domains in the cholesterol-rich (32.5 mol%) liquid-ordered ([Formula: see text]) phase of dipalmitoylphosphatidylcholine membranes. The liquid ordered phase in one-component lipid membranes has previously been thought to be a homogeneous phase. The presence of highly ordered lipid domains embedded in a disordered lipid matrix implies non-uniform distribution of cholesterol between the two phases. The experimental results are in excellent agreement with recent computer simulations of DPPC/cholesterol complexes [Meinhardt, Vink and Schmid (2013). Proc Natl Acad Sci USA 110(12): 4476-4481], which reported the existence of nanometer size [Formula: see text] domains in a liquid disordered lipid environment.

Tocopherol activity correlates with its location in a membrane: a new perspective on the antioxidant vitamin E

We show evidence of an antioxidant mechanism for vitamin E which correlates strongly with its physical location in a model lipid bilayer. These data address the overlooked problem of the physical distance between the vitamin's reducing hydrogen and lipid acyl chain radicals. Our combined data from neutron diffraction, NMR, and UV spectroscopy experiments all suggest that reduction of reactive oxygen species and lipid radicals occurs specifically at the membrane's hydrophobic-hydrophilic interface. The latter is possible when the acyl chain "snorkels" to the interface from the hydrocarbon matrix. Moreover, not all model lipids are equal in this regard, as indicated by the small differences in vitamin's location. The present result is a clear example of the importance of lipid diversity in controlling the dynamic structural properties of biological membranes. Importantly, our results suggest that measurements of aToc oxidation kinetics, and its products, should be revisited by taking into consideration the physical properties of the membrane in which the vitamin resides.

Permeability of fluid-phase phospholipid bilayers: assessment and useful correlations for permeability screening and other applications

Permeability data (P(lip/w) ) for liquid crystalline phospholipid bilayers composed of egg lecithin and dimyristoylphosphatidylcholine (DMPC) are analyzed in terms of a mathematical model that accounts for free surface area and chain-ordering effects in the bilayer as well as size and lipophilicity of the permeating species. Free surface area and chain ordering are largely determined by temperature and cholesterol content of the membrane, molecular size is represented by molecular weight, and lipophilicity of the barrier region is represented by the 1,9-decadiene/water partition coefficient, following earlier work by Xiang, Anderson, and coworkers. A correlating variable χ = MW(n) σ/(1 -σ) is used to link the results from different membrane systems, where different values of n are tried, and σ denotes a reduced phospholipid density. The group (1 -σ)/σ is a measure of free surface area, but can also be interpreted in terms of free volume. A single exponential function of χ is developed that is able to correlate 39 observations of P(lip/w) for different compounds in egg lecithin at low density, and 22 observations for acetic acid in DMPC at higher densities, spanning nine orders of magnitude to within an rms error for log 10 P(lip/w) of 0.20. The best fit found for n = 0.87 ultimately makes χ much closer to the ratio of molecular to free volumes than surface areas. The results serve as a starting point for estimating passive permeability of cell membranes to nonionized solutes as a function of temperature and cholesterol content of the membrane.

Materials Science and Engineering of the Low Temperature Sensitive Liposome (LTSL): Composition-Structure-Property Relationships That Underlie its Design and Performance

This chapter presents the material science and materials engineering concepts that went into the design and testing of the Low Temperature-Sensitive Liposome (LTSL), including: the roles of each of the components that make up the composite membrane; how the molecular and nanostructures that they form might influence the already anomalous permeability at the phase transition of the bilayer; and how this thermally sensitive “Smart Drug Delivery System” leads to ultrafast release of a loaded doxorubicin drug, triggered and controlled in the micro-vasculature of tumors by applied mild hyperthermia. This formulation approach, as ThermoDox®, has been used in a completed 700-patient Phase III human clinical trial in liver cancer (HEAT study), is in a Phase II trial in chest wall recurrence of cancer (DIGNITY study) and has been used in a Phase I trial of patients with colorectal liver metastases (ABLATE study). With additional research and preclinical studies underway, and a range of other drugs, imaging agents and biological modifiers poised for encapsulation, the LTSL could provide a new paradigm for drug and agent delivery for the treatment of localized tumors: rapid triggered drug release in the tumor bloodstream and deep penetration of drug into the tumor tissue.

Using small-angle neutron scattering to detect nanoscopic lipid domains

The cell plasma membrane is a complex system, which is thought to be capable of exhibiting non-random lateral organization. Studies of live cells and model membranes have yielded mechanisms responsible for the formation, growth, and maintenance of nanoscopic heterogeneities, although the existence and mechanisms that give rise to these heterogeneities remain controversial. Small-angle neutron scattering (SANS) is a tool ideally suited to interrogate lateral heterogeneity in model membranes, primarily due to its unique spatial resolution (i.e., ~5-100nm) and its ability to resolve structure with minimal perturbation to the membrane. In this review we examine several methods used to analyze the SANS signal arising from freely suspended unilamellar vesicles containing lateral heterogeneity. Specifically, we discuss an analytical model for a single, round domain on a spherical vesicle. We then discuss a numerical method that uses Monte Carlo simulation to describe systems with multiple domains and/or more complicated morphologies. Also discussed are several model-independent approaches that are sensitive to membrane heterogeneity. The review concludes with several recent applications of SANS to the study of membrane raft mixtures.

The backbone dynamics of the amyloid precursor protein transmembrane helix provides a rationale for the sequential cleavage mechanism of γ-secretase

The etiology of Alzheimer's disease depends on the relative abundance of different amyloid-β (Aβ) peptide species. These peptides are produced by sequential proteolytic cleavage within the transmembrane helix of the 99 residue C-terminal fragment of the amyloid precursor protein (C99) by the intramembrane protease γ-secretase. Intramembrane proteolysis is thought to require local unfolding of the substrate helix, which has been proposed to be cleaved as a homodimer. Here, we investigated the backbone dynamics of the substrate helix. Amide exchange experiments of monomeric recombinant C99 and of synthetic transmembrane domain peptides reveal that the N-terminal Gly-rich homodimerization domain exchanges much faster than the C-terminal cleavage region. MD simulations corroborate the differential backbone dynamics, indicate a bending motion at a diglycine motif connecting dimerization and cleavage regions, and detect significantly different H-bond stabilities at the initial cleavage sites. Our results are consistent with the following hypotheses about cleavage of the substrate: First, the GlyGly hinge may precisely position the substrate within γ-secretase such that its catalytic center must start proteolysis at the known initial cleavage sites. Second, the ratio of cleavage products formed by subsequent sequential proteolysis could be influenced by differential extents of solvation and by the stabilities of H-bonds at alternate initial sites. Third, the flexibility of the Gly-rich domain may facilitate substrate movement within the enzyme during sequential proteolysis. Fourth, dimerization may affect substrate processing by decreasing the dynamics of the dimerization region and by increasing that of the C-terminal part of the cleavage region.

A Coarse-Grained Molecular Dynamics Study of DLPC, DMPC, DPPC, and DSPC Mixtures in Aqueous Solution

The structural and dynamics properties of the bilayer comprising 128 molecules of dipalmitoylphosphatidylcholine (DPPC), dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), and distearoylphosphatidylcholine (DSPC) in water were investigated using a coarse-grained molecular dynamics (CG-MD) simulation technique. The model mixture system was simulated at 298 K under semi-isotropic pressure conditions. The aggregation was initiated from the random configurations followed by the formation of a bilayer over a period of 500 ns. The calculated values of the area per lipid, thickness, and lateral diffusion for the mixed model were different from when a single lipid was used. Our results confirmed that the chain length of the lipid molecules strongly affects the phospholipid bilayer’s physical properties.

Introductory lecture: basic quantities in model biomembranes

One of the many aspects of membrane biophysics dealt with in this Faraday Discussion regards the material moduli that describe energies at a supramolecular level. This introductory lecture first critically reviews differences in reported numerical values of the bending modulus K(C), which is a central property for the biologically important flexibility of membranes. It is speculated that there may be a reason that the shape analysis method tends to give larger values of K(C) than the micromechanical manipulation method or the more recent X-ray method that agree very well with each other. Another theme of membrane biophysics is the use of simulations to provide exquisite detail of structures and processes. This lecture critically reviews the application of atomic level simulations to the quantitative structure of simple single component lipid bilayers and diagnostics are introduced to evaluate simulations. Another theme of this Faraday Discussion was lateral heterogeneity in biomembranes with many different lipids. Coarse grained simulations and analytical theories promise to synergistically enhance experimental studies when their interaction parameters are tuned to agree with experimental data, such as the slopes of experimental tie lines in ternary phase diagrams. Finally, attention is called to contributions that add relevant biological molecules to bilayers and to contributions that study the exciting shape changes and different non-bilayer structures with different lipids.

Interactions between ether phospholipids and cholesterol as determined by scattering and molecular dynamics simulations

Cholesterol and ether lipids are ubiquitous in mammalian cell membranes, and their interactions are crucial in ether lipid mediated cholesterol trafficking. We report on cholesterol's molecular interactions with ether lipids as determined using a combination of small-angle neutron and X-ray scattering, and all-atom molecular dynamics (MD) simulations. A scattering density profile model for an ether lipid bilayer was developed using MD simulations, which was then used to simultaneously fit the different experimental scattering data. From analysis of the data the various bilayer structural parameters were obtained. Surface area constrained MD simulations were also performed to reproduce the experimental data. This iterative analysis approach resulted in good agreement between the experimental and simulated form factors. The molecular interactions taking place between cholesterol and ether lipids were then determined from the validated MD simulations. We found that in ether membranes cholesterol primarily hydrogen bonds with the lipid headgroup phosphate oxygen, while in their ester membrane counterparts cholesterol hydrogen bonds with the backbone ester carbonyls. This different mode of interaction between ether lipids and cholesterol induces cholesterol to reside closer to the bilayer surface, dehydrating the headgroup's phosphate moiety. Moreover, the three-dimensional lipid chain spatial density distribution around cholesterol indicates anisotropic chain packing, causing cholesterol to tilt. These insights lend a better understanding of ether lipid-mediated cholesterol trafficking and the roles that the different lipid species have in determining the structural and dynamical properties of membrane associated biomolecules.