The effect of chain length on protein solubilization in polymer-based vesicles (polymersomes)

Using a mean-field analysis we derive a consistent model for the perturbation of a symmetric polymeric bilayer due to the incorporation of transmembrane proteins, as a function of the polymer molecular weight and the protein dimensions. We find that the mechanism for the inhibition of protein incorporation in polymeric bilayers differs from that of their inclusion in polymer-carrying lipid vesicles; in polymersomes, the equilibrium concentration of transmembrane proteins decreases as a function of the thickness mismatch between the protein and the bilayer core, whereas in liposomes the presence of polymer chains affects the protein adsorption kinetics. Despite the increased stiffness of polymer bilayers (when compared to lipid ones), their perturbation decay length and range of protein-protein interaction is found to be relatively long. The energetic penalty due to protein adsorption increases relatively slowly as a function of the polymer chain length due to the self-assembled nature of the polymer bilayer. As a result, we predict that transmembrane proteins may be incorporated in significant numbers even in bilayers where the thickness mismatch is large.     

Probing the lipid-protein interface using model transmembrane peptides with a covalently linked acyl chain

The aim of this study was to gain insight into how interactions between proteins and lipids in membranes are sensed at the protein-lipid interface. As a probe to analyze this interface, we used deuterium-labeled acyl chains that were covalently linked to a model transmembrane peptide. First, a perdeuterated palmitoyl chain was coupled to the Trp-flanked peptide WALP23 (Ac-CGWW(LA)₈LWWA-NH₂), and the deuterium NMR spectrum was analyzed in di-C18:1-phosphatidylcholine (PC) bilayers. We found that the chain order of this peptide-linked chain is rather similar to that of a noncovalently coupled perdeuterated palmitoyl chain, except that it exhibits a slightly lower order. Similar results were obtained when site-specific deuterium labels were used and when the palmitoyl chain was attached to the more-hydrophobic model peptide WLP23 (Ac-CGWWL₁₇WWA-NH₂) or to the Lys-flanked peptide KALP23 (Ac-CGKK(LA)₈LKKA-NH₂). The experiments showed that the order of both the peptide-linked chains and the noncovalently coupled palmitoyl chains in the phospholipid bilayer increases in the order KALP23 < WALP23 < WLP23. Furthermore, changes in the bulk lipid bilayer thickness caused by varying the lipid composition from di-C14:1-PC to di-C18:1-PC or by including cholesterol were sensed rather similarly by the covalently coupled chain and the noncovalently coupled palmitoyl chains. The results indicate that the properties of lipids adjacent to transmembrane peptides mostly reflect the properties of the surrounding lipid bilayer, and hence that (at least for the single-span model peptides used in this study) annular lipids do not play a highly specific role in protein-lipid interactions.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Effect of hydrostatic pressure on the bilayer phase behavior of symmetric and asymmetric phospholipids with the same total chain length

The bilayer phase transitions of palmitoylstearoyl-phosphatidylcholine (PSPC), diheptadecanoyl-PC (C17PC) and stearoylpalmitoyl-PC (SPPC) which have the same total carbon numbers in the two acyl chains were observed by differential scanning calorimetry and high-pressure optical method. As the temperature increased, these bilayers exhibited four phases of the subgel (L_c), lamellar gel (L_β′), ripple gel (P_β′) and liquid crystal (L_α), in turn. The L_c phase was observed only in the first heating scan after cold storage. The temperatures of the phase transitions were almost linearly elevated by applying pressure. The temperature-pressure phase diagrams and the thermodynamic quantities associated with the phase transitions were compared among the lipid bilayers. For all the bilayers studied, the pressure-induced interdigitated gel (L_βI) phase appeared above the critical interdigitation pressure (CIP) between the L_β′ and P_β′ phases. The CIPs for the PSPC, C17PC and SPPC bilayers were found to be 50.6, 79.1 and 93.0 MPa, respectively. Contribution of two acyl chains to thermodynamic properties for the phase transitions of asymmetric PSPC and SPPC bilayers was not even. The sn-2 acyl chain lengths of asymmetric PCs governed primarily the bilayer properties. The fluorescence spectra of Prodan in lipid bilayers showed the emission maxima characteristic of bilayer phases, which were dependent on the location of Prodan in the bilayers. Second derivative of fluorescent spectrum exhibited the original emission spectrum of Prodan to be composed of the distribution of Prodan into multiple locations in the lipid bilayer. The F″₄₉₇/F″₄₃₀ value, a ratio of second derivative of fluorescence intensity at 497 nm to that at 430 nm, is decisive evidence whether bilayer interdigitation will occur. With respect to the L_β′/L_βI phase transition in the SPPC bilayer, the emission maximum of Prodan exhibited the narrow-range red-shift from 441 to 449 nm, indicating that the L_βI phase in the SPPC bilayer has a less polar “pocket” formed by a space between uneven terminal methyl ends of the sn-1 and sn-2 chains, in which the Prodan molecule remains stably.     

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.     

New structural model for mixed-chain phosphatidylcholine bilayers

Multilamellar suspensions of a mixed-chain saturated phosphatidylcholine with 18 carbon atoms in the sn-1 chain and 10 carbon atoms in the sn-2 chain have been analyzed by X-ray diffraction techniques. The structural parameters for this lipid in the gel state are quite different than usual phosphatidylcholine bilayer phases. A symmetric and sharp wide-angle reflection at 4.11 A indicates that the hydrocarbon chains in hydrated C(18):C(10)PC bilayers are more tightly packed than in usual gel-state phosphatidylcholine bilayers and that there is no hydrocarbon chain tilt. The lipid thickness is about 12 A smaller than would be expected in a normal bilayer phase, and the area per molecule is 3 times the area per hydrocarbon chain. In addition, the bilayer thickness increases upon melting to the liquid-crystalline state, whereas normal bilayer phases decrease in thickness upon melting. On the basis of these data, we propose a new lipid packing model for gel-state C(18):C(10)PC bilayers in which the long C(18) chain spans the entire width of the hydrocarbon region of the bilayer and the short C(10) chain aligns or abuts with the C(10) chain from the apposing molecule. This model is novel in that there are three hydrocarbon chains per head group at the lipid-water interface. Calculations show that this phase is energetically favorable for mixed-chain lipids provided the long acyl chain is nearly twice the length of the shorter chain. In the liquid-crystalline state C(18):C(10)PC forms a normal fluid bilayer, with two chains per head group.(ABSTRACT TRUNCATED AT 250 WORDS)     

The 2004 Biophysical Society-Avanti Award in Lipids address: roles of bilayer structure and elastic properties in peptide localization in membranes

This review details how bilayer structural/elastic properties impact three distinct areas of biological significance. First, the partitioning of melittin into bilayers and melittin-induced bilayer leakage depended strongly on bilayer composition. The incorporation of cholesterol into phosphatidylcholine bilayers decreased melittin-induced leakage from 73 to 3%, and bilayers composed of lipopolysaccharide (LPS), the main lipid on the surface of Gram-negative bacteria, also had low (3%) melittin-induced leakage. Second, transbilayer peptides of different hydrophobic lengths were largely excluded from bilayer microdomains (“rafts”) enriched in sphingomyelin (SM) and cholesterol, even when the length of the transbilayer peptide domain matched the hydrocarbon thickness of the raft bilayer. This is likely due to the large area compressibility modulus of SM:cholesterol bilayers. Third, the major water barrier of skin, the extracellular lamellae of the stratum corneum, was found to contain tightly packed asymmetric lipid bilayers with cholesterol located preferentially on one side of the bilayer and a unique skin ceramide containing an unsaturated acyl chain on the opposite side. We argue that, in each of these three areas, key factors are differences in lipid hydrocarbon chain packing for different lipids, particularly the tight hydrocarbon chain packing caused by cholesterol’s strong interaction with saturated chains.     

Lipid lateral diffusion and membrane heterogeneity

The pulsed field gradient (pfg)-NMR method for measurements of translational diffusion of molecules in macroscopically aligned lipid bilayers is described. This technique is proposed to have an appreciable potential for investigations in the field of lipid and membrane biology. Transport of molecules in the plane of the bilayer can be successfully studied, as well as lateral phase separation of lipids and their dynamics within the bilayer organizations. Lateral diffusion coefficients depend on lipid packing and acyl chain ordering and investigations of order parameters of perdeuterated acyl chains, using ²H NMR quadrupole splittings, are useful complements. In this review we summarize some of our recent achievements obtained on lipid membranes. In particular, bilayers exhibiting two-phase coexistence of liquid disordered (l_d) and liquid ordered (l_o) phases are considered in detail. Methods for obtaining good oriented lipid bilayers, necessary for the pfg-NMR method to be efficiently used, are also briefly described. Among our major results, besides determinations of l_d and l_o phases, belongs the finding that the lateral diffusion is the same for all components, independent of the molecular structure (including cholesterol (CHOL)), if they reside in the same domain or phase in the membrane. Furthermore, quite unexpectedly CHOL seems to partition into the l_dand l_o phases to roughly the same extent, indicating that CHOL has no strong preference for any of these phases, i.e. CHOL seems to have similar interactions with all of the lipids. We propose that the lateral phase separation in bilayers containing one high-T_m and one low-T_m lipid together with CHOL is driven by the increasing difficulty of incorporating an unsaturated or prenyl lipid into the highly ordered bilayer formed by a saturated lipid and CHOL, i.e. the phase transition is entropy driven to keep the disorder of the hydrocarbon chains of the unsaturated lipid.     

Rhodopsin/Lipid Hydrophobic Matching—Rhodopsin Oligomerization and Function

Lipid composition of the membrane and rhodopsin packing density strongly modulate the early steps of the visual response of photoreceptor membranes. In this study, lipid-order and bovine rhodopsin function in proteoliposomes composed of the sn-1 chain perdeuterated lipids 14:0_d27-14:1-PC, 16:0_d31-16:1-PC, 18:0_d35-18:1-PC, or 20:0_d39-20:1-PC at rhodopsin/lipid molar ratios from 1:70 to 1:1000 (mol/mol) were investigated. Clear evidence for matching of hydrophobic regions on rhodopsin transmembrane helices and hydrophobic thickness of lipid bilayers was observed from ²H nuclear magnetic resonance order parameter measurements at low rhodopsin concentrations. Thin bilayers stretched to match the length of transmembrane helices observed as increase of sn-1 chain order, while thicker bilayers were compressed near the protein. A quantitative analysis of lipid-order parameter changes suggested that the protein adjusts its conformation to bilayer hydrophobic thickness as well, which confirmed our earlier circular-dichroism measurements. Changes in lipid order parameters upon rhodopsin incorporation vanished for bilayers with a hydrophobic thickness of 27 ± 1 Å, suggesting that this is the bilayer thickness at which rhodopsin packs in bilayers at the lowest membrane perturbation. The lipid-order parameter studies also indicated that a hydrophobic mismatch between rhodopsin and lipids triggers rhodopsin oligomerization with increasing rhodopsin concentrations. Both hydrophobic mismatch and rhodopsin oligomerization result in substantial shifts of the equilibrium between the photointermediates metarhodopsin I and metarhodopsin II; increasing bilayer thickness favors formation of metarhodopsin II while oligomerization favors metarhodopsin I. The results highlight the importance of hydrophobic matching for rhodopsin structure, oligomerization, and function.     

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.     

Hydrophobic thickness, lipid surface area and polar region hydration in monounsaturated diacylphosphatidylcholine bilayers: SANS study of effects of cholesterol and β-sitosterol in unilamellar vesicles

The influence of a mammalian sterol cholesterol and a plant sterol β-sitosterol on the structural parameters and hydration of bilayers in unilamellar vesicles made of monounsaturated diacylphosphatidylcholines (diCn:1PC, n = 14-22 is the even number of acyl chain carbons) was studied at 30 °C using small-angle neutron scattering (SANS). Recently published advanced model of lipid bilayer as a three-strip structure was used with a triangular shape of polar head group probability distribution (Ku?erka et al., Models to analyze small-angle neutron scattering from unilamellar lipid vesicles, Physical Review E 69 (2004) Art. No. 051903). It was found that 33 mol% of both sterols increased the thickness of diCn:1PC bilayers with n = 18-22 similarly. β-sitosterol increased the thickness of diC14:1PC and diC16:1PC bilayers a little more than cholesterol. Both sterols increased the surface area per unit cell by cca 12 Å² and the number of water molecules located in the head group region by cca 4 molecules, irrespective to the acyl chain length of diCn:1PC. The structural difference in the side chain between cholesterol and β-sitosterol plays a negligible role in influencing the structural parameters of bilayers studied.     

Properties of unsaturated phospholipid bilayers: Effect of cholesterol

Properties of hydrated unsaturated phosphatidylcholine (PC) lipid bilayers containing 40 mol % cholesterol and of pure PC bilayers have been studied. Various methods were applied, including molecular dynamics simulations, self-consistent field calculations, and the pulsed field gradient nuclear magnetic resonance technique. Lipid bilayers were composed of 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules. Lateral self-diffusion coefficients of the lipids in all these bilayers, mass density distributions of atoms and atom groups with respect to the bilayer normal, the C-H and C-C bond order parameter profiles of each phospholipid hydrocarbon chain with respect to the bilayer normal were calculated. It was shown that the lateral self-diffusion coefficient of PC molecules of the lipid bilayer containing 40 mol % cholesterol is smaller than that for a corresponding pure PC bilayer; the diffusion coefficients increase with increasing the degree of unsaturation of one of the PC chains in bilayers of both types (i.e., in pure bilayers or in bilayers with cholesterol). The presence of cholesterol in a bilayer promoted the extension of saturated and polyunsaturated lipid chains. The condensing effect of cholesterol on the order parameters was more pronounced for the double C=C bonds of polyunsaturated chains than for single C-C bonds of saturated chains.     

Atomic force microscopy of supported lipid membranes and their complexes with polycations

The method of atomic force microscopy has been used to investigate the morphology of mica-supported bilayer lipid membranes and stability of their complexes with a cationic polymer, poly-(N-ethyl-4-vinylpyridinium bromide). Lipid bilayers with a minimum of defects were obtained by the fusion of monolamellar neutral or mixed anionic bilayer vesicles (liposomes) on the mica surface, followed by excessive solvent removal by means of rapid rotation of a plate in horizontal plane (spin-coating). It has been shown that the cationic polymer does not interact with the bilayers, where the outer leaflet (i.e., the monolayer exposed to the surrounding aqueous solution) is made of an electroneutral phosphatidylcholine (PC). At the same time, the polymer irreversibly binds to the bilayer containing an anionic lipid.     

Properties of Hexadecaprenyl Monophosphate/Dioleoylphosphatidylcholine Vesicular Lipid Bilayers

In our study we investigated hemispherical phospholipid bilayer membranes and phospholipid vesicles made from hexadecaprenyl monophosphate (C₈₀-P), dioleoylphosphatidylocholine (DOPC) and their mixtures by voltammetric and transmission electron microscopy (TEM) techniques. The current-voltage characteristics, the membrane conductance-temperature relationships and the membrane breakdown voltage have been measured for different mixtures of C₈₀-P/DOPC. The membrane hydrophobic thickness and the activation energy of ion migration across the membrane have been determined. Hexadecaprenyl monophosphate decreased in comparison with DOPC bilayers, the membrane conductance, increased the activation energy and the membrane breakdown voltage for the various value of C₈₀-P/DOPC mole ratio, respectively. The TEM micrographs of C₈₀-P, DOPC and C₈₀-P/DOPC lipid vesicles showed several characteristic structures, which have been described. The data indicate that hexadecaprenyl monophosphate modulates the surface curvature of the membranes by the formation of aggregates in liquid-crystalline phospholipid membranes. We suggest that the dynamics and conformation of hexadecaprenyl monophosphate in membranes depend on the transmembrane electrical potential. The electron micrographs indicate that polyprenyl monophosphates with single isoprenyl chains form lipid vesicular bilayers. The thickness of the bilayer, evaluated from the micrographs, was 11 ± 1 nm. This property creates possibility of forming primitive bilayer lipid membranes by long single-chain polyprenyl phosphates in abiotic conditions. It can be the next step in understanding the origin of protocells. Received: 10 January 2000/Revised: 7 June 2000     

Polymerization in black lipid membranes

A variety of different lipids containing dienoyl groups in the side chains were tested for membrane formation using the planar lipid bilayer approach. One of these lipids formed stable bilayers which could be polymerized using UV-illumination. The influence of the polymerization was studied in monolayers, lipid vesicles and planar bilayers. The stability of the lipid bilayer membranes was increased by polymerization. Thus, the lifetime of the membranes increased from about 1 h to 4–5 h or longer. Furthermore, the specific conductance of unmodified membranes and of carrier-mediated transport is reduced. The transport of lipophilic ions was investigated as a function of polymerization using the charge-pulse method. The absorption of dipicrylamine (DPA^-) is not affected. The translocation of this compound and of tetraphenylborate (B(Ph) ₄ ^- ) showed a strong decrease with polymerization time. The influence of polymerization on the membrane structure may be explained on the basis of a strong viscosity increase in the lipid bilayer membrane.     

Influence of lipid membrane rigidity on properties of supporting polymer

Temperature-sensitive hydrogel polymers are utilized as responsive layers in various applications. Although the polymer's native characteristics have been studied extensively, details concerning its properties during interaction with biorelated structures are lacking. This work investigates the interaction between a thermoresponsive polymer cushion and different lipid membrane capping layers probed by neutron reflectometry. N-isopropylacrylamide copolymerized with methacroylbenzophenone first supported a lipid bilayer composed of 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE) and subsequently 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). The polymer-membrane systems were investigated above and below the polymer transition temperature (37 and 25°C). Although the same cushion supported each lipid membrane, the polymer hydration profile and thickness were markedly different for DPPE and DPPC systems. Because DPPE and DPPC have different bending rigidities, these results establish that the polymer-membrane interaction is critically mediated by the mechanics of the membrane, providing better insight into cell-hydrogel interactions.     

Influence of N-dodecyl-N,N-dimethylamine N-oxide on the activity of sarcoplasmic reticulum Ca(2+)-transporting ATPase reconstituted into diacylphosphatidylcholine vesicles: efects of bilayer physical parameters

Sarcoplasmic reticulum Ca-transporting ATPase (EC 3.6.1.38) was isolated from rabbit white muscle, purified and reconstituted into vesicles of synthetic diacylphosphatidylcholines with monounsaturated acyl chains using the cholate dilution method. In fluid bilayers at 37 degrees C, the specific activity of ATPase displays a maximum (31.5+/-0.8 IU/mg) for dioleoylphosphatidylcholine (diC18:1PC) and decreases progressively for both shorter and longer acyl chain lengths. Besides the hydrophobic mismatch between protein and lipid bilayer, changes in the bilayer hydration and lateral interactions detected by small angle neutron scattering (SANS) can contribute to this acyl chain length dependence. When reconstituted into dierucoylphosphatidylcholine (diC22:1PC), the zwitterionic surfactant N-dodecyl-N,N-dimethylamine N-oxide (C12NO) stimulates the ATPase activity from 14.2+/-0.6 to 32.5+/-0.8 IU/mg in the range of molar ratios C12NO:diC22:1PC=0/1.2. In dilauroylphosphatidylcholines (diC12:0PC) and diC18:1PC, the effect of C12NO is twofold-the ATPase activity is stimulated at low and inhibited at high C12NO concentrations. In diC18:1PC, it is observed an increase of activity induced by C12NO in the range of molar ratios C12NO:diC18:1PC< or =1.3 in bilayers, where the bilayer thickness estimated by SANS decreases by 0.4+/-0.1 nm. In this range, the 31P-NMR chemical shift anisotropy increases indicating an effect of C12NO on the orientation of the phosphatidylcholine dipole N(+)-P- accompanied by a variation of the local membrane dipole potential. A decrease of the ATPase activity is observed in the range of molar ratios C12NO:diC18:1PC=1.3/2.5, where mixed tubular micelles are detected by SANS in C12NO+diC18:1PC mixtures. It is concluded that besides hydrophobic thickness changes, the changes in dipole potential and curvature frustration of the bilayer could contribute as well to C12NO effects on Ca(2+)-ATPase activity.     

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 (D_B), hydrophobic thickness (D_C), and lipid volume (V_L). The only minor discrepancy between simulation and experiment is a measure of (D_B − D_HH) / 2 where D_HH 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.     

Phase behavior of hydrated lipid bilayer composed of binary mixture of phospholipids with different head groups

Phase behavior of hydrated lipid bilayer was investigated for the mixtures of two phospholipid species chosen from phosphatidic acid (PA), phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylglycerol (PG) with the same acyl chains. The pseudo-binary phase diagrams constructed by a differential scanning calorimetry (DSC) were analyzed based on a thermodynamic model applying the Bragg–Williams approximation for non-ideality of mixing. The interchange energy parameters, ρ₀, derived from this approach were positive for all mixture systems in both gel and liquid–crystalline phase bilayers, and increased in the order PG/PE<PC/PA<PC/PE<PG/PA with a few exception. This suggests that the energetical disadvantage for the mixed-pair formation relative to the like-pair formation in the hydrated bilayer increases in this order. In addition, the ρ₀ values increased with the increase in the acyl chain length of the phospholipids. These experimental results were discussed in terms of an intermolecular interaction of the phospholipid species in hydrated bilayer.     

Mean-field calculations of chain packing and conformational statistics in lipid bilayers: comparison with experiments and molecular dynamics studies.

A molecular, mean-field theory of chain packing statistics in aggregates of amphiphilic molecules is applied to calculate the conformational properties of the lipid chains comprising the hydrophobic cores of dipalmitoyl-phosphatidylcholine (DPPC), dioleoyl-phosphatidylcholine (DOPC), and palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers in their fluid state. The central quantity in this theory, the probability distribution of chain conformations, is evaluated by minimizing the free energy of the bilayer assuming only that the segment density within the hydrophobic region is uniform (liquidlike). Using this distribution we calculate chain conformational properties such as bond orientational order parameters and spatial distributions of the various chain segments. The lipid chains, both the saturated palmitoyl (-(CH2)14-CH3) and the unsaturated oleoyl (-(CH2)7-CH = CH-(CH2)7-CH3) chains are modeled using rotational isomeric state schemes. All possible chain conformations are enumerated and their statistical weights are determined by the self-consistency equations expressing the condition of uniform density. The hydrophobic core of the DPPC bilayer is treated as composed of single (palmitoyl) chain amphiphiles, i.e., the interactions between chains originating from the same lipid headgroup are assumed to be the same as those between chains belonging to different molecules. Similarly, the DOPC system is treated as a bilayer of oleoyl chains. The POPC bilayer is modeled as an equimolar mixture of palmitoyl and oleoyl chains. Bond orientational order parameter profiles, and segment spatial distributions are calculated for the three systems above, for several values of the bilayer thickness (or, equivalently, average area/headgroup) chosen, where possible, so as to allow for comparisons with available experimental data and/or molecular dynamics simulations. In most cases the agreement between the mean-field calculations, which are relatively easy to perform, and the experimental and simulation data is very good, supporting their use as an efficient tool for analyzing a variety of systems subject to varying conditions (e.g., bilayers of different compositions or thicknesses at different temperatures).