Detergent solubilization of phosphatidylcholine bilayers in the fluid state: Influence of the acyl chain structure

Seventeen different, chemically defined phosphatidylcholines, dispersed in aqueous medium in the form of large unilamellar vesicles, have been tested for solubilization by the non-ionic detergent Triton X-100. The temperatures (either 20 °C or 45 °C) were such that the bilayers were always in the liquid-disordered state. For each case, the solubilization parameters, D_on (total detergent: lipid mole ratio producing the onset of solubilization) and D₅₀ (total detergent: lipid mole ratio producing 50% solubilization), were determined under equilibrium conditions. Both parameters varied generally in parallel. When double bonds were introduced to the acyl chains, other factors remaining constant, solubilization became more difficult, i.e., more detergent was required. Cis-unsaturated phospholipids required more detergent than the corresponding trans-isomers. Increasing chain length in saturated phospholipids between C12 and C16 decreased moderately the detergent/lipid ratios causing solubilization. Acyl and alkyl phospholipids were equally susceptible to Triton X-100 solubilization. Lipid chain order, as measured by DPH fluorescence polarization, seemed to facilitate solubilization, perhaps because more ordered bilayers have a smaller capacity to accommodate detergent monomers without breaking down into lipid-detergent mixed micelles.     

Interaction of small peptides with lipid bilayers.

Molecular dynamics simulations of the tripeptide Ala-Phe-Ala-O-tert-butyl interacting with dimyristoylphosphatidylcholine lipid bilayers have been carried out. The lipid and aqueous environments of the peptide, the alkyl chain order, and the lipid and peptide dynamics have been investigated with use of density profiles, radial distribution functions, alkyl chain order parameter profiles, and time correlation functions. It appears that the alkyl chain region accommodates the peptides in the bilayer with minimal perturbation to this region. The peptide dynamics in the bilayer bound form has been compared with that of the free peptide in water. The peptide structure does not vary on the simulation time scale (of the order of hundreds of picoseconds) compared with the solution structure in which a random structure is observed.     

Solid-state NMR and MD simulations of the antiviral drug amantadine solubilized in DMPC bilayers

The interactions of ¹⁵N-labeled amantadine, an antiinfluenza A drug, with DMPC bilayers were investigated by solid-state NMR and by a 12.6-ns molecular dynamics (MD) simulation. The drug was found to assume a single preferred orientation and location when incorporated in these bilayers. The experimental and MD computational results demonstrate that the long axis of amantadine is on average parallel to the bilayer normal, and the amine group is oriented toward the headgroups of the lipid bilayers. The localization of amantadine was determined by paramagnetic relaxation and by the MD simulation showing that amantadine is within the interfacial region and that the amine interacts with the lipid headgroup and glycerol backbone, while the hydrocarbon portion of amantadine interacts with the glycerol backbone and much of the fatty acyl chain as it wraps underneath the drug. The lipid headgroup orientation changes on drug binding as characterized by the anisotropy of ³¹P chemical shielding and ¹⁴N quadrupolar interactions and by the MD simulation.     

Distribution of Pentachlorophenol in Phospholipid Bilayers: A Molecular Dynamics Study

Molecular dynamics computer simulations of pentachlorophenol (PCP) in palmitoyl-oleoyl-phosphatidylethanolamine and palmitoyl-oleoyl-phosphatidylcholine lipid bilayers were carried out to investigate the distribution of PCP and the effects of PCP on the phospholipid bilayer structure. Starting from two extreme starting structures, including PCP molecules outside the lipid bilayer, the PCP distribution converges in simulations of up to 50 ns. PCP preferentially occupies the region between the carbonyl groups and the double bonds in the acyl chains of the lipid molecules in the bilayer. In the presence of PCP, the lipid chain order increases somewhat in both chains, and the average tilt angle of the lipid chains decreases. The increase in the lipid chain order in the presence of PCP was more pronounced in the palmitoyl-oleoyl-phosphatidylcholine bilayer compared to the palmitoyl-oleoyl-phosphatidylethanolamine bilayer. The number of trans conformations of lipid chain dihedrals does not change significantly. PCP aligns parallel to the alkyl chains of the lipid to optimize the packing in the dense ordered chain region of the bilayer. The hydroxyl group of PCP forms hydrogen bonds with both water and lipid oxygen atoms in the water/lipid interface region.     

Lipid and Peptide Dynamics in Membranes upon Insertion of n-alkyl-β-D-Glucopyranosides

The effect of nonionic detergents of the n-alkyl-β-D-glucopyranoside class on the ordering of lipid bilayers and the dynamics of membrane-embedded peptides were investigated with ²H- and ³¹P-NMR. 1,2-dipalmitoyl-sn-glycero-3-phosphocholine was selectively deuterated at methylene segments C-2, C-7, and C-16 of the two fatty acyl chains. Two trans-membrane helices, WALP-19 and glycophorin A_71-98, were synthesized with Ala-d₃ in the central region of the α-helix. n-Alkyl-β-D-glucopyranosides with alkyl chains with 6, 7, 8, and 10 carbon atoms were added at increasing concentrations to the lipid membrane. The bilayer structure is retained up to a detergent/lipid molar ratio of 1:1. The insertion of the detergents leads to a selective disordering of the lipids. The headgroup region remains largely unaffected; the fatty acyl chain segments parallel to the detergent alkyl chain are only modestly disordered (10-20%), whereas lipid segments beyond the methyl terminus of the detergent show a decrease of up to 50%. The change in the bilayer order profile corresponds to an increase in bilayer entropy. Insertion of detergents into the lipid bilayers is completely entropy-driven. The entropy change accompanying lipid disorder is equivalent in magnitude to the hydrophobic effect. Ala-d₃ deuterated WALP-19 and GlycA_71-97 were incorporated into bilayers of 1,2-dimyristoyl-sn-glycero-3-phosphocholine at a peptide/lipid molar ratio of 1:100 and measured above the 1,2-dimyristoyl-sn-glycero-3-phosphocholine gel/liquid-crystal phase transition. Well-resolved ²H-NMR quadrupole splittings were observed for the two trans-membrane helices, revealing a rapid rotation of the CD₃ methyl rotor superimposed on an additional rotation of the whole peptide around the bilayer normal. The presence of detergent fluidizes the membrane and produces magnetic alignment of bilayer domains but does not produce essential changes in the peptide conformation or dynamics.     

Molecular dynamics simulations of charged and neutral lipid bilayers: treatment of electrostatic interactions

Molecular dynamics (MD) simulations complement experimental methods in studies of the structure and dynamics of lipid bilayers. The choice of algorithms employed in this computational method represents a trade-off between the accuracy and real calculation time. The largest portion of the simulation time is devoted to calculation of long-range electrostatic interactions. To speed-up evaluation of these interactions, various approximations have been used. The most common ones are the truncation of long-range interactions with the use of cut-offs, and the particle-mesh Ewald (PME) method. In this study, several multi-nanosecond cut-off and PME simulations were performed to establish the influence of the simulation protocol on the bilayer properties. Two bilayers were used. One consisted of neutral phosphatidylcholine molecules. The other was a mixed lipid bilayer consisting of neutral phosphatidylethanolamine and negatively charged phosphatidylglycerol molecules. The study shows that the cut-off simulation of a bilayer containing charge molecules generates artefacts; in particular the mobility and order of the charged molecules are vastly different from those determined experimentally. In the PME simulation, the bilayer properties are in general agreement with experimental data. The cut-off simulation of bilayers containing only uncharged molecules does not generate artefacts, nevertheless, the PME simulation gives generally better agreement with experimental data.     

Lipid tails modulate antimicrobial peptide membrane incorporation and activity

Membrane disrupting antimicrobial peptides (AMPs) are often amphipathic peptides that interact directly with lipid bilayers. AMPs are generally thought to interact mostly with lipid head groups, but it is less clear how the lipid alkyl chain length and saturation modulate interactions with membranes. Here, we used native mass spectrometry to measure the stoichiometry of three different AMPs—LL-37, indolicidin, and magainin-2—in lipid nanodiscs. We also measured the activity of these AMPs in unilamellar vesicle leakage assays. We found that LL-37 formed specific hexamer complexes but with different intermediates and affinities that depended on the bilayer thickness. LL-37 was also most active in lipid bilayers containing longer, unsaturated lipids. In contrast, indolicidin incorporated to a higher degree into more fluid lipid bilayers but was more active with bilayers with thinner, less fluid lipids. Finally, magainin-2 incorporated to a higher degree into bilayers with longer, unsaturated alkyl chains and showed more activity in these same conditions. Together, these data show that higher amounts of peptide incorporation generally led to higher activity and that AMPs tend to incorporate more into longer unsaturated lipid bilayers. However, the activity of AMPs was not always directly related to amount of peptide incorporated.     

Structure refinement of membrane proteins via molecular dynamics simulations

A refinement protocol based on physics‐based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge‐based or implicit membrane‐based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid‐facing residues. Scoring with knowledge‐based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane‐based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.     

Dynorphin induced magnetic ordering in lipid bilayers as studied by P NMR spectroscopy

Lipid bilayers of dimyristoyl phosphatidylcholine (DMPC) containing opioid peptide dynorphin A(1-17) are found to be spontaneously aligned to the applied magnetic field near at the phase transition temperature between the gel and liquid crystalline states (T_m=24°C), as examined by ³¹P NMR spectroscopy. The specific interaction between the peptide and lipid bilayer leading to this property was also examined by optical microscopy, light scattering, and potassium ion-selective electrode, together with a comparative study on dynorphin A(1-13). A substantial change in the light scattering intensity was noted for DMPC containing dynorphin A(1-17) near at T_m but not for the system containing A(1-13). Besides, reversible change in morphology of bilayer, from small lipid particles to large vesicles, was observed by optical microscope at T_m. These results indicate that lysis and fusion of the lipid bilayers are induced by the presence of dynorphin A(1-17). It turned out that the bilayers are spontaneously aligned to the magnetic field above T_m in parallel with the bilayer surface, because a single ³¹P NMR signal appeared at the perpendicular position of the ³¹P chemical shift tensor. In contrast, no such magnetic ordering was noted for DMPC bilayers containing dynorphin A(1-13). It was proved that DMPC bilayer in the presence of dynorphin A(1-17) forms vesicles above T_m, because leakage of potassium ion from the lipid bilayers was observed by potassium ion-selective electrode after adding Triton X-100. It is concluded that DMPC bilayer consists of elongated vesicles with the long axis parallel to the magnetic field, together with the data of microscopic observation of cylindrical shape of the vesicles. Further, the long axis is found to be at least five times longer than the short axis of the elongated vesicles in view of simulated ³¹P NMR lineshape.     

Shaped Apertures in Photoresist Films Enhance the Lifetime and Mechanical Stability of Suspended Lipid Bilayers

Planar lipid bilayers suspended in apertures provide a controlled environment for ion channel studies. However, short lifetimes and poor mechanical stability of suspended bilayers limit the experimental throughput of bilayer electrophysiology experiments. Although bilayers are more stable in smaller apertures, ion channel incorporation through vesicle fusion with the suspended bilayer becomes increasingly difficult. In an alternative bilayer stabilization approach, we have developed shaped apertures in SU8 photoresist that have tapered sidewalls and a minimum diameter between 60 and 100 μm. Bilayers formed at the thin tip of these shaped apertures, either with the painting or the folding method, display drastically increased lifetimes, typically >20 h, and mechanical stability, being able to withstand extensive perturbation of the buffer solution. Single-channel electrical recordings of the peptide alamethicin and of the proteoliposome-delivered potassium channel KcsA demonstrate channel conductance with low noise, made possible by the small capacitance of the 50 μm thick SU8 septum, which is only thinned around the aperture, and unimpeded proteoliposome fusion, enabled by the large aperture diameter. We anticipate that these shaped apertures with micrometer edge thickness can substantially enhance the throughput of channel characterization by bilayer lipid membrane electrophysiology, especially in combination with automated parallel bilayer platforms.     

Lipids, curvature stress, and the action of lipid prodrugs: free fatty acids and lysolipid enhancement of drug transport across liposomal membranes

Molecular shape and its impact on bilayer curvature stress are powerful concepts for describing the effects of lipids and fatty acids on fundamental membrane properties, such as passive permeability and derived properties like drug transport across liposomal membranes. We illustrate these relationships by studying the effects of fatty acids and lysolipids on the permeation of a potent anti-cancer drug, doxorubicin, across the bilayer of a liposome in which the drug is encapsulated. Using a simple fluorescence assay, we have systematically studied the passive permeation of doxorubicin across liposomal membranes in different lipid phases: the solid-ordered phase (DPPC bilayers), the liquid-disordered phase (POPC lipid bilayers), and the liquid-ordered phase induced by high levels of cholesterol (DOPC + cholesterol lipid bilayers). The effect of different free fatty acids (FA) and lysolipids (LL), separately and in combination, on permeability was assessed to elucidate the possible mechanism of phospholipase A₂-triggered release in cancer tissue of liposomal doxorubicin formulations. In all cases, FAs applied separately lead to significant enhancement of permeability, most pronounced in liquid-disordered bilayers and less pronounced in solid and solid-ordered bilayers. LLs applied separately had only a marginal effect on permeability. FA and LL applied in combination lead to a synergistic enhancement of permeability in solid bilayers, whereas in liquid-disordered bilayers, the combined effect suppressed the otherwise strong permeability enhancement due to the FAs.     

Comparative study of molecular dynamics,diffusion, and permeability for ligands in biomembranes of different lipid composition

A comparative study of several model lipid bilayers of different composition, which included analysis of kinetic parameters of model lipid bilayers and permeability of bilayer membranes for small molecules, has been carried out. The conformity of results of numeric experiments to experimental data (structure of membrane lipid bilayers, lateral diffusion coefficients, and relative permeability of biomembranes for ligands) is discussed in the framework of a standard molecular dynamics protocol.     

Surface areas of lipid membranes

Upon photolysis, alkyl pentacyanocobaltate complexes generate alkyl radicals which react rapidly and specifically with nitroxide radicals, and which do not penetrate phospholipid bilayers. By measuring the loss of paramagnetic resonance signal intensity when multilamellar liposomes containing a small amount of spin-labeled lipid are exposed to these radicals, we have measured the proportion of lipid on the external surface of liposomes. We have shown that liposomes prepared under specified conditions from dimyristoylphosphatidylcholine, dipalmitoylphosphatidylcholine, and binary mixtures of dipalmitoylphosphatidylcholine and cholesterol all have the same proportion of external lipid.     

Structure of Sphingomyelin Bilayers: A Simulation Study

We have carried out a molecular dynamics simulation of a hydrated 18:0 sphingomyelin lipid bilayer. The bilayer contained 1600 sphingomyelin (SM) molecules, and 50,592 water molecules. After construction and initial equilibration, the simulation was run for 3.8 ns at a constant temperature of 50°C and a constant pressure of 1 atm. We present properties of the bilayer calculated from the simulation, and compare with experimental data and with properties of dipalmitoyl phosphatidylcholine (DPPC) bilayers. The SM bilayers are significantly more ordered and compact than DPPC bilayers at the same temperature. SM bilayers also exhibit significant intramolecular hydrogen bonding between phosphate ester oxygen and hydroxyl hydrogen atoms. This results in a decreased hydration in the polar region of the SM bilayer compared with DPPC. Since our simulation system is very large we have calculated the power spectrum of bilayer undulation and peristaltic modes, and we compare these data with similar calculations for DPPC bilayers. We find that the SM bilayer has significantly larger bending modulus and area compressibility compared to DPPC.     

Lycopene,but not zeaxanthin,serves as a skeleton for the formation of an orthorhombic organization of intercellular lipids within the lamellae in the stratum corneum: Molecular dynamics simulations of the hydrated ceramide NS bilayer model

Carotenoids play an important role in the protection of biomembranes against oxidative damage. Their function depends on the surroundings and the organization of the lipid membrane they are embedded in. Carotenoids are located parallel or perpendicular to the surface of the lipid bilayer. The influence of carotenoids on the organization of the lipid bilayer in the stratum corneum has not been thoroughly considered. Here, the orientation of the exemplary cutaneous carotenoids lycopene and zeaxanthin in a hydrated ceramide NS24 bilayer model and the influence of carotenoids on the lateral organization of the lipid bilayer model were studied by means of molecular dynamics simulations for 32 °C and 37 °C. The results confirm that lycopene is located parallel and zeaxanthin perpendicular to the surface of the lipid bilayer. The lycopene-loaded lipid bilayer appeared to have a strong orthorhombic organization, while zeaxanthin-loaded and pure lipid bilayers were organized in a disordered hexagonal-like and liquid-like state, respectively. The effect is stronger at 32 °C compared to 37 °C based on p-values. Therefore, it was assumed that carotenoids without hydroxyl polar groups in their structure facilitate the formation of the orthorhombic organization of lipids, which provides the skin barrier function. It was shown that the distance between carotenoid atoms matched the distance between atoms in the lipids, indicating that parallel located carotenoids without hydroxyl groups serve as a skeleton for lipid membranes inside the lamellae. The obtained results provide reasonable prediction of the overall qualitative properties of lipid model systems and show the importance of parallel-oriented carotenoids in the development and maintenance of the skin barrier function.     

Disturb or Stabilize? A Molecular Dynamics Study of the Effects of Resorcinolic Lipids on Phospholipid Bilayers

Resorcinolic lipids, or resorcinols, are commonly found in plant membranes. They consist of a substituted benzene ring forming the hydrophilic lipid head, attached to an alkyl chain forming the hydrophobic tail. Experimental results show alternative effects of resorcinols on lipid membranes. Depending on whether they are added to lipid solutions before or after the formation of the liposomes, they either stabilize or destabilize these liposomes. Here we use atomistic molecular dynamics simulations to elucidate the molecular nature of this dual effect. Systems composed of either one of three resorcinol homologs, differing in the alkyl tail length, interacting with dimyristoylphosphatidylcholine lipid bilayers were studied. It is shown that resorcinols preincorporated into bilayers induce order within the lipid acyl chains, decrease the hydration of the lipid headgroups, and make the bilayers less permeable to water. In contrast, simulations in which the resorcinols are incorporated from the aqueous solution into a preformed phospholipid bilayer induce local disruption, leading to either transient pore formation or even complete rupture of the membrane. In line with the experimental data, our simulations thus demonstrate that resorcinols can either disturb or stabilize the membrane structure, and offer a detailed view of the underlying molecular mechanism.     

Evidence of pores and thinned lipid bilayers induced in oriented lipid membranes interacting with the antimicrobial peptides, magainin-2 and aurein-3.3

Dynamic structures of supramolecular lipid assemblies, such as toroidal pores and thinned bilayers induced in oriented lipid membranes, which are interacting with membrane-acting antimicrobial peptides (AMPs), magainin-2 and aurein-3.3, were explored by ³¹P and ²H solid-state NMR (ssNMR) spectroscopy. Various types of phospholipid systems, such as POPC-d₃₁, POPC-d₃₁/POPG, and POPC-d₃₁/cholesterol, were investigated to understand the membrane disruption mechanisms of magainin-2 and aurein-3.3 peptides at various peptide-to-lipid (P:L) ratios. The experimental lineshapes of anisotropic ³¹P and ²H ssNMR spectra measured on these peptide-lipid systems were simulated reasonably well by assuming the presence of supramolecular lipid assemblies, such as toroidal pores and thinned bilayers, in membranes. Furthermore, the observed decrease in the anisotropic frequency span of either ³¹P or ²H ssNMR spectra of oriented lipid bilayers, particularly when anionic POPG lipids are interacting with AMPs at high P:L ratios, can directly be explained by a thinned membrane surface model with fast lateral diffusive motions of lipids. The spectral analysis protocol we developed enables extraction of the lateral diffusion coefficients of lipids distributed on the curved surfaces of pores and thinned bilayers on a few nanometers scale.     

Areas of Monounsaturated Diacylphosphatidylcholines

We have studied the structural properties of monounsaturated diacylphosphatidylcholine lipid bilayers (i.e., diCn:1PC, where n = 14, 16, 18, 20, 22, and 24 is the number of acyl chain carbons). High-resolution x-ray scattering data were analyzed in conjunction with contrast-varied neutron scattering data using a technique we recently developed. Analyses of the data show that the manner by which bilayer thickness increases with increasing n in monounsaturated diacylphosphatidylcholines is dependent on the double bond's position. For commonly available monounsaturated diacylphosphatidylcholines, this results in the nonlinear behavior of both bilayer thickness and lipid area, whereas for diC18:1PC bilayers, lipid area assumes a maximum value. It is worthwhile to note that compared to previous data, our results indicate that lipid areas are smaller by ∼10%. This observation highlights the need to revisit lipid areas, as they are often used in comparisons with molecular dynamics simulations. Moreover, simulators are encouraged to compare their results not only to x-ray scattering data, but to neutron data as well.     

Development of a CP P NMR Broadline Simulation Methodology for Studying the Interactions of Antihypertensive AT1 Antagonist Losartan with Phospholipid Bilayers

A cross-polarization (CP) ³¹P NMR broadline simulation methodology was developed for studying the effects of drugs in phospholipids bilayers. Based on seven-parameter fittings, this methodology provided information concerning the conformational changes and dynamics effects of losartan in the polar region of the dipalmitoylphosphatidylcholine bilayers. The test molecule for this study was losartan, an antihypertensive drug known to exert its effect on AT₁ transmembrane receptors. The results were complemented and compared with those of differential scanning calorimetry, solid-state ¹³C NMR spectroscopy, Raman spectroscopy, and electron spin resonance. More specifically, these physical chemical methodologies indicated that the amphipathic losartan molecule interacts with the hydrophilic-head zone of the lipid bilayers. The CP ³¹P NMR broadline simulations showed that the lipid molecules in the bilayers containing losartan displayed greater collective tilt compared to the tilt displayed by the load-free bilayers, indicating improved packing. The Raman results displayed a decrease in the trans/gauche ratio and increased intermolecular interactions of the acyl chains in the liquid crystalline phase. Additional evidence, suggesting that losartan possibly anchors in the realm of the headgroup, was derived from upfield shift of the average chemical shift σ^iso of the ³¹P signal in the presence of losartan and from shift of the observed peak at 715 cm⁻¹ attributed to C-N stretching in the Raman spectra.     

Study of the topography of cannabinoids in model membranes using X-ray diffraction

Small-angle X-ray diffraction was used to determine the topography of (-)-delta 8-tetrahydrocannabinol in partially hydrated dimyristoylphosphatidylcholine bilayers. Electron density profiles of lipid bilayers in the presence and absence of the cannabinoid were calculated using Fourier transform. Step-function equivalent profiles were then constructed to obtain the absolute electron density scale. We have compared the electron density profiles of the above preparations to determine the location of the drug molecule in the bilayer. By using (-)-5'-iodo-delta 8-tetrahydrocannabinol in parallel experiments, we were also able to locate the iodine atom in the bilayer and deduce the conformation of the cannabinoid side alkyl chain. All comparisons were made between different preparations having the same mesomorphic form and total period repeat distance. To achieve this, we have carried out X-ray diffraction experiments at various temperatures to cover the different mesomorphic phases and combined our data with the corresponding results from differential scanning calorimetry. Based on the results of this work and previous data on the orientation of the cannabinoid in model membranes, we concluded that the phenolic hydroxy group of the drug molecule exists near the carbonyl groups of DMPC and that the average position of the iodine atom is approx. 5.5 A from the center (terminal methyl region) of the DMPC bilayer. This requires the cannabinoid side-chain to assume an orientation parallel to the bilayer chains.