Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. I. Scaling of neutron data and the distributions of double bonds and water.

We described in two previous papers a method for the joint refinement of the structure of fluid bilayers using neutron and x-ray diffraction data (Wiener, M. C., and S. H. White 1991a, b. Biophys. J. 59: 162-173 and 174-185). An essential part of the method is the appropriate scaling of the diffraction data. Here we describe the scaling of the neutron data and the determination of the transbilayer distribution of double bonds in liquid-crystalline (L alpha phase) phospholipid bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The distribution was determined by neutron diffraction of oriented multilayers (66% RH) of DOPC specifically deuterated at the 9- and 10-position of both acyl chains. The double-bond distribution is described accurately by a pair of Gaussian functions each located at a position Zcc = 7.88 +/- 0.09 A from the bilayer center with 1/e-halfwidths of Acc = 4.29 +/- 0.16 A. Previously, we determined the transbilayer distribution of bromine atoms in a specifically halogenated lipid, 1-oleoyl-2-9,10-dibromostearoyl-sn-glycero-3-phosphocholine (OBPC), and showed it to be an isomorphous replacement for DOPC (Wiener, M. C., and S. H. White, 1991c. Biochemistry. In press). A comparison of the double-bond and bromine profiles indicates that the positions of the centers of the deuterated double bond and the brominated methylene Gaussian distributions are equal within experimental error and that each label undergoes similar average thermal motions with respect to the bilayer normal.(ABSTRACT TRUNCATED AT 250 WORDS)     

Transbilayer distribution of bromine in fluid bilayers containing a specifically brominated analogue of dioleoylphosphatidylcholine

We describe in this paper the transbilayer distribution of the bromines of the specifically halogenated phospholipid 1-oleoyl-2-(9,10-dibromostearoyl)-sn-glycero-3- phosphocholine (OBPC). The distribution was determined by X-ray diffraction of oriented multilayers of mixtures of OBPC and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) at 66% relative humidity by the general approach of Franks et al. (1978) [Nature 276, 530-532]. The bromine distribution of OBPC in the fluid L alpha phase is described accurately by a pair of Gaussian functions located 7.97 +/- 0.27 A from the center of the bilayer with l/e half-widths of 4.96 +/- 0.62 A. We find that OBPC bilayers are accurately described as DOPC bilayers with an additional bromine distribution centered at the position of the double bond of DOPC and conclude that OBPC is an excellent structural isomorph for DOPC under the conditions of these experiments. The distribution obtained is the complete and fully resolved transbilayer image of the halogen label because the broad distribution of the bromines is due entirely to thermal disorder and not to experimental limitations [Wiener, M. C., & White, S. H. (1991a) Biophys. J. 59, 162-173]. The observed width of the bromine distribution indicates that virtually all of the hydrocarbon interior is accessible to the bromines. The distance between the bromine/double-bond position and the headgroup phosphate position was determined from one-dimensional Patterson maps and found to be approximately 12 A. The application of accurately determined bromine distributions to the quantitative interpretation of fluorescence quenching experiments is discussed. A method for the self-consistent global analysis of diffraction data from mixtures that permits the use of data sets with different instrumental scale factors is developed in an Appendix.     

Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. II. Distribution and packing of terminal methyl groups.

We continue in this paper the presentation of theoretical and experimental methods for the joint refinement of neutron and x-ray lamellar diffraction data for the analysis of fluid (L alpha phase) bilayer structure (Wiener, M. C., and S. H. White. 1991 a, b, c. Biophys. J. 59:162-173 and 174-185; Biochemistry. 30:6997-7008; Wiener, M. C., G. I. King, and S. H. White. Biophys. J. 60: 568-576). We show how to obtain the distribution and packing of the terminal methyls in the interior of a fluid dioleoylphosphatidylcholine bilayer (66% RH) by combining x-ray and neutron scattering-length transbilayer profiles with no a priori assumptions about the functional form of the distribution. We find that the methyls can be represented by a Gaussian function with 1/e-halfwidth of 2.95 +/- 0.28 A situated at the bilayer center. There is substantial mixing of the methyls and methylenes in the bilayer center. The Gaussian representation of the methyl distribution is narrower and has a different shape than predicted by several simulations of fluid bilayers (Gruen, D. W. R., and E. H. B. de Lacey. 1984. Surfactants in Solution, Vol. 1. Plenum Publishing Corp., New York. 279-306; de Loof, H., et al. 1991. Biochemistry. 30:2099-2133) but this may be due to the smaller area/lipid of our experiments and the presence of the double-bonds. Determination of the absolute specific volume of DOPC and an analysis of bulk alkane volumetric data over a range of hydrostatic pressures lead to estimates of methylene and methyl volumes at the bilayer center of 27 +/- 1 A3 and 57.2 +/- 3.6 A3, respectively. This result provides direct confirmation of the common assumption that the molecular packing of methyl and methylene groups in bilayers is the same as in bulk liquid alkanes.     

Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. III. Complete structure.

We present in this paper the complete structure of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in the L alpha phase (66% RH, 23 degrees C) obtained by the joint refinement of neutron and x-ray lamellar diffraction data. The structural details obtained have previously required a large number of neutron diffraction experiments, using numerous specifically-deuterated phospholipid isomorphs (Büldt et al., 1978. Nature (Lond.). 271:182-184). The joint-refinement approach minimizes specific deuteration by utilizing independent neutron and x-ray data sets. The method yields a quasimolecular structure consisting of a series of multiatomic fragments that are each represented by one or several Gaussian distributions whose positions and widths can be determined to within 0.06 to 0.52 A exclusive of the methylene region. The image of DOPC at 66% RH (5.36 +/- 0.08 waters per lipid) is consistent with many aspects of bilayer structure previously determined by structural and spectroscopic studies. The most striking feature of the structure is the large amount of transbilayer thermal motion suggested by the widths and overlaps of the Gaussian envelopes of the quasimolecular fragments. We discuss the "dynamic bilayer thickness" which describes the minimum effective thickness of the hydrocarbon permeability barrier in terms of the thermal motion of the water. A gradient of thermal motion exists that increases in either direction away from the glycerol backbone which is the most constrained portion of the bilayer. The steric interactions between headgroups of apposed bilayers, expected at the hydration level of our experiments, are clearly revealed. A useful consequence of the quasimolecular structure is that average boundaries within bilayers calculated using composition and volumetric data and ad hoc assumptions can be related to the positions of the principal structural groups. Several measures of "bilayer thickness" in common use can be identified as the positions of the cholines for Luzzati's d1 (Luzzati and Husson. 1962. J. Cell Biol. 12:207-219) and the glycerols for Small's dL (Small. 1967. J. Lipid Res. 8:551-556). We do not know if these relations will be true at other hydrations or for other lipids. Of particular interest is the fact that the position of the carbonyl groups marks the average hydrocarbon/headgroup boundary. It must be emphasized, however, that this region of the bilayer must be generally characterized as one of tumultuous chemical heterogeneity because of the thermal motion of the bilayer.     

Diffraction-based density restraints for membrane and membrane-peptide molecular dynamics simulations

We have recently shown that current molecular dynamics (MD) atomic force fields are not yet able to produce lipid bilayer structures that agree with experimentally-determined structures within experimental errors. Because of the many advantages offered by experimentally validated simulations, we have developed a novel restraint method for membrane MD simulations that uses experimental diffraction data. The restraints, introduced into the MD force field, act upon specified groups of atoms to restrain their mean positions and widths to values determined experimentally. The method was first tested using a simple liquid argon system, and then applied to a neat dioleoylphosphatidylcholine (DOPC) bilayer at 66% relative humidity and to the same bilayer containing the peptide melittin. Application of experiment-based restraints to the transbilayer double-bond and water distributions of neat DOPC bilayers led to distributions that agreed with the experimental values. Based upon the experimental structure, the restraints improved the simulated structure in some regions while introducing larger differences in others, as might be expected from imperfect force fields. For the DOPC-melittin system, the experimental transbilayer distribution of melittin was used as a restraint. The addition of the peptide caused perturbations of the simulated bilayer structure, but which were larger than observed experimentally. The melittin distribution of the simulation could be fit accurately to a Gaussian with parameters close to the observed ones, indicating that the restraints can be used to produce an ensemble of membrane-bound peptide conformations that are consistent with experiments. Such ensembles pave the way for understanding peptide-bilayer interactions at the atomic level.     

Area per molecule and distribution of water in fully hydrated dilauroylphosphatidylethanolamine bilayers

The area per lipid molecule for fully hydrated dilauroylphosphatidylethanolamine (DLPE) has been obtained in both the gel and liquid-crystalline states by combining wide-angle X-ray diffraction, electron density profiles, and previously published dilatometry results [Wilkinson, D. A., & Nagle, J. F. (1981) Biochemistry 20, 187-192]. The molecular area increases from 41.0 +/- 0.2 to 49.1 +/- 1.2 A2 upon melting from the gel to liquid-crystalline phase. The thickness of the bilayer, as measured from the electron density profiles, decreases about 4 A upon melting, from 45.2 +/- 0.3 to 41.0 +/- 0.6 A. A somewhat unexpected result is that the fluid layer between fully hydrated bilayers is the same in both gel and liquid-crystalline phases and is only about 5 A thick. From these data, plus the volume of the anhydrous DLPE molecule, it is possible to determine the number of water molecules per lipid and their approximate distribution relative to the lipid molecule. Our analysis shows that there are about 7 and 9 waters per DLPE molecule in the gel and liquid-crystalline phases, respectively. About half of the water is located in the fluid space between adjacent bilayers, and the remaining waters are intercalated into the bilayer, presumably in the head group region. There are significantly fewer water molecules in the fluid spaces between DLPE bilayers than in the fluid spaces in gel- or liquid-crystalline-phase phosphatidylcholine bilayers. This small fluid space in PE bilayers could arise from interbilayer hydrogen bond formation through the water molecules or electrostatic interactions between the amine and phosphate groups on apposing bilayers.     

Acyl-chain methyl distributions of liquid-ordered and -disordered membranes

A central feature of the lipid raft concept is the formation of cholesterol-rich lipid domains. The introduction of relatively rigid cholesterol molecules into fluid liquid-disordered (L_d) phospholipid bilayers can produce liquid-ordered (L_o) mixtures in which the rigidity of cholesterol causes partial ordering of the flexible hydrocarbon acyl chains of the phospholipids. Several lines of evidence support this concept, but direct structural information about L_o membranes is lacking. Here we present the structure of L_o membranes formed from cholesterol and dioleoylphosphatidylcholine (DOPC). Specific deuteration of the DOPC acyl-chain methyl groups and neutron diffraction measurements reveal an extraordinary disorder of the acyl chains of neat L_d DOPC bilayers. The disorder is so great that >20% of the methyl groups are in intimate contact with water in the bilayer interface. The ordering of the DOPC acyl chains by cholesterol leads to retraction of the methyl groups away from the interface. Molecular dynamics simulations based on experimental systems reveal asymmetric transbilayer distributions of the methyl groups associated with each bilayer leaflet.     

Decrease of elastic moduli of DOPC bilayers induced by a macrolide antibiotic, azithromycin

The elastic properties of membrane bilayers are key parameters that control its deformation and can be affected by pharmacological agents. Our previous atomic force microscopy studies revealed that the macrolide antibiotic, azithromycin, leads to erosion of DPPC domains in a fluid DOPC matrix [A. Berquand, M. P. Mingeot-Leclercq, Y. F. Dufrene, Real-time imaging of drug-membrane interactions by atomic force microscopy, Biochim. Biophys. Acta 1664 (2004) 198-205.]. Since this observation could be due to an effect on DOPC cohesion, we investigated the effect of azithromycin on elastic properties of DOPC giant unilamellar vesicles (GUVs). Microcinematographic and morphometric analyses revealed that azithromycin addition enhanced lipid membranes fluctuations, leading to eventual disruption of the largest GUVs. These effects were related to change of elastic moduli of DOPC, quantified by the micropipette aspiration technique. Azithromycin decreased both the bending modulus (k(c), from 23.1+/-3.5 to 10.6+/-4.5 k(B)T) and the apparent area compressibility modulus (K(app), from 176+/-35 to 113+/-25 mN/m). These data suggested that insertion of azithromycin into the DOPC bilayer reduced the requirement level of both the energy for thermal fluctuations and the stress to stretch the bilayer. Computer modeling of azithromycin interaction with DOPC bilayer, based on minimal energy, independently predicted that azithromycin (i) inserts at the interface of phospholipid bilayers, (ii) decreases the energy of interaction between DOPC molecules, and (iii) increases the mean surface occupied by each phospholipid molecule. We conclude that azithromycin inserts into the DOPC lipid bilayer, so as to decrease its cohesion and to facilitate the merging of DPPC into the DOPC fluid matrix, as observed by atomic force microscopy. These investigations, based on three complementary approaches, provide the first biophysical evidence for the ability of an amphiphilic antibiotic to alter lipid elastic moduli. This may be an important determinant for drug: lipid interactions and cellular pharmacology.     

Formation of supported phospholipid bilayers via co-adsorption with beta-D-dodecyl maltoside

We have investigated the formation of supported model membranes via the adsorption of phospholipid-surfactant mixtures at the Si-water interface by specular neutron reflection. The adsorption of mixed micelles of the nonionic surfactant beta-D-dodecyl maltoside and DOPC or POPC was determined as a function of bulk concentration, and using d25-beta-D-dodecyl maltoside, the composition of DOPC and POPC bilayers was determined. Bilayer thicknesses of 39+/-3 A for DOPC and 41+/-3 A for POPC agree well with data from bulk lamellar phases for both lipids, and the average area per lipid molecule can be varied from 62 to 115 A2 by varying the bulk concentrations used. The amount of surfactant in the bilayer is very sensitive to the bulk volume-to-surface area ratio, but it can be fully eliminated by ensuring a sufficiently large dilution/rinsing volume of the solution.     

Abnormal transbilayer mobility of phosphatidylcholine in hereditary pyropoikilocytosis reflects the increased heat sensitivity of the membrane skeleton

We determined whether the membrane defect in hereditary pyropoikilocytosis (HPP) is associated with thermally induced changes in the lipid bilayer, the stability of which was probed by the rate of translocation of phosphatidylcholine (PC) over the two leaflets. [14C]PC was incorporated into the outer leaflet of the lipid bilayer of the intact erythrocytes using a PC-specific phospholipid exchange protein. The transbilayer equilibration of this PC was determined by measuring the time-dependent changes in its accessibility to exogenous phospholipase A2. The rate of transbilayer equilibration of PC was increased in HPP cells at 37 degrees C when compared to normal erythrocytes (rate constants, 0.07 +/- 0.02 and 0.03 +/- 0.01 h-1, respectively). A further dramatic increase in PC transbilayer equilibration was noted in HPP cells incubated at 44 degrees C (rate constant, 0.15 +/- 0.02 h-1). A similar marked acceleration in transbilayer movement of PC was also seen in normal erythrocytes when incubated at 46 degrees C (rate constant, 0.13 +/- 0.03 h-1). Despite the enhanced transbilayer mobility of PC in HPP cells when compared to normal erythrocytes, no major alteration in the asymmetric distribution could be observed when probed with phospholipase A2. Since changes in transbilayer mobility of PC and cell morphology occur in HPP cells at lower temperature than in normal red cells, it may be concluded that the enhanced thermal sensitivity of spectrin is the major factor responsible for these changes. Our results therefore support the view that the structural integrity of the skeletal network is essential for stabilization of the lipid bilayer of the red cell membrane.     

Formation of "solvent-free" black lipid bilayer membranes from glyceryl monooleate dispersed in squalene.

A simple technique for forming "black" lipid bilayer membranes containing negligible amounts of alkyl solvent is described. The membranes are formed by the method of Mueller et al (Circulation. 1962. 26:1167.) from glyceryl monooleate (GMO) dispersed in squalene. The squalene forms an annulus to satisfy the boundary conditions of the planar bilayer but does not appear to dissolve noticeably in the bilayer itself. The specific geometric capacitance (Cg) of the membranes at 20 degrees C formed by this technique is 0.7771 +/- 0.0048 muF/cm2. Theoretical estimates of Cg for solvent-free bilayers range from 0.75 to 0.81 muF/cm2. Alkane-free GMO bilayers formed from n-octadecane by the solvent freeze-out method of White (Biochim. Biophys. Acta. 1974. 356:8) have values of Cg = 0.7903 +/- 0.0013 muF/cm2 at 20.5 degrees C. The agreement between the various values of Cg strongly suggests that the bilayers are free of squalene. DC potentials applied to the bilayers have no detectable effect on the value of Cg, as expected for solvent-free films. The ability to form bilayers essentially free of the solvent used in the forming solution makes it possible to determine the area per molecule of the surface active lipid in the bilayer. The area per molecule of GMO at 20 degrees C is estimated to be 37.9 +/- 0.2 A2.     

The nature of the hydrophobic binding of small peptides at the bilayer interface: implications for the insertion of transbilayer helices

One method of obtaining useful information about the physical chemistry of peptide/bilayer interactions is to relate thermodynamic parameters of the interactions to structural parameters obtained by diffraction methods. We report here the results of the application of this approach to interactions of hydrophobic tripeptides of the form Ala-X-Ala-O-tert-butyl with lipid bilayers. The thermodynamic constants (delta Gt, delta Ht, and delta St) for the transfer of the tripeptides from water into DMPC vesicles were determined for X = Leu, Phe, and Trp and found to be consistent with those expected for hydrophobic interactions above the phase transition of DMPC. Combining these results with the earlier ones of Jacobs and White [(1986) Biochemistry 25, 2605-2612], the favorable free energies of transfer with different amino acids in the -X- position increase in the order Gly less than Ala less than Leu less than Phe less than Trp in agreement with the Nozaki and Tanford [(1971) J. Biol. Chem. 246, 2211-2217] hydrophobicity scale. Determination of the location of Ala-[2H5]Trp-Ala-O-tert-butyl in oriented DOPC bilayers by neutron diffraction shows that the most hydrophobic peptide of the series is confined to the bilayer headgroup/water region. Refinement of the diffraction measurements shows that only 13% of the tryptophan is associated with the hydrocarbon core. The distribution of the water tends to mirror that of the peptide. Unlike peptide-free bilayers, 5% of the water penetrates the hydrocarbon, which is about 100-fold greater than expected. A quantitative thermodynamic analysis of the interfacial binding of the peptides suggests that (1) the hydrophobic interactions are 60-70% complete upon binding at the bilayer interface, (2) the interface is likely to play an important role in helix formation and insertion, (3) the hydrogen bond status of amino acid side chains is crucial to insertion, and (4) an a priori lack of knowledge of the status of such bonds could limit the precision of hydrophobicity plots. We introduce an interfacial hydrophobicity scale, IFH(h), with a variable hydrogen bond parameter (h) that permits one to consider explicitly hydrogen bonding in transbilayer helix searches.     

Differing modes of interaction between monomeric Aβ(1-40) peptides and model lipid membranes: an AFM study

Membrane interactions with β-amyloid peptides are implicated in the pathology of Alzheimer's disease and cholesterol has been shown to be key modulator of this interaction, yet little is known about the mechanism of this interaction. Using atomic force microscopy, we investigated the interaction of monomeric Aβ(1-40) peptides with planar mica-supported bilayers composed of DOPC and DPPC containing varying concentrations of cholesterol. We show that below the bilayer melting temperature, Aβ monomers adsorb to, and assemble on, the surface of DPPC bilayers to form layers that grow laterally and normal to the bilayer plane. Above the bilayer melting temperature, we observe protofibril formation. In contrast, in DOPC bilayers, Aβ monomers exhibit a detergent-like action, forming defects in the bilayer structure. The kinetics of both modes of interaction significantly increases with increasing membrane cholesterol content. We conclude that the mode and rate of the interaction of Aβ monomers with lipid bilayers are strongly dependent on lipid composition, phase state and cholesterol content.     

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.     

Membrane-thinning effect of curcumin

Interaction of curcumin with lipid bilayers is not well understood. A recent experiment showed that curcumin significantly affected the single-channel lifetime of gramicidin in a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayer without affecting its single-channel conductance. We performed two experiments to understand this result. By isothermal titration calorimetry, we measured the partition coefficient of curcumin binding to DOPC bilayers. By x-ray lamellar diffraction, we measured the thickness change of DOPC bilayers as a function of the curcumin/lipid ratio. A nonlinear membrane-thinning effect by curcumin was discovered. The gramicidin data were qualitatively interpreted by the combination of isothermal titration calorimetry and x-ray results. We show that not only does curcumin thin the lipid bilayer, it might also weaken its elasticity moduli. The result implies that curcumin may affect the function of membrane proteins by modifying the properties of the host membrane.     

Effect of amphiphilic surfactant LDAO on the solubilization of DOPC vesicles and on the activity of Ca(2+)-ATPase reconstituted in DOPC vesicles

Solubilization of large unilamellar 1,2-dioleoylphosphatidylcholine (DOPC) vesicles by N-dodecyl-N,N-dimethylamine-N-oxide (LDAO) was studied using turbidimetry. From turbidity data, the LDAO partition coefficient between the aqueous phase and DOPC bilayers was obtained. Using this partition coefficient, the LDAO:DOPC molar ratio in the bilayer was calculated and effects of LDAO on the bilayer stability, bilayer thickness and on the phosphohydrolase activity of sarcoplasmic reticulum Ca(2+) transporting ATPase (SERCA) reconstituted into DOPC were compared at the same LDAO:DOPC molar ratios in the bilayer. The sequence "bilayers in vesicles - bilayer fragments (flat mixed micelles) - tubular mixed micelles - globular mixed micelles" was suggested for the solubilization mechanism of DOPC vesicles from the combined turbidimetric and small-angle neutron scattering (SANS) results. The effective molecular packing parameter delta = 0.5, corresponding to the mixed bilayer - mixed tubular micelle transition, was calculated from fragmental DOPC and LDAO volumes at the molar ratio LDAO:DOPC = 2.00 in bilayers, in the middle of transition region observed earlier experimentally by small-angle neutron scattering (SANS). The bilayer thickness decrease induced by LDAO in DOPC observed by SANS did not result in the SERCA phosphohydrolase activity decrease and this indicates that some other factors compensated this bilayer effect of LDAO. The ATPase activity decrease at higher LDAO concentrations was caused by the bilayer deformation. This deformation resulted in the formation of non-bilayer aggregates in LDAO+DOPC system.     

Structural determinants of water permeability through the lipid membrane

Despite intense study over many years, the mechanisms by which water and small nonelectrolytes cross lipid bilayers remain unclear. While prior studies of permeability through membranes have focused on solute characteristics, such as size, polarity, and partition coefficient in hydrophobic solvent, we focus here on water permeability in seven single component bilayers composed of different lipids, five with phosphatidylcholine headgroups and different chain lengths and unsaturation, one with a phosphatidylserine headgroup, and one with a phosphatidylethanolamine headgroup. We find that water permeability correlates most strongly with the area/lipid and is poorly correlated with bilayer thickness and other previously determined structural and mechanical properties of these single component bilayers. These results suggest a new model for permeability that is developed in the accompanying theoretical paper in which the area occupied by the lipid is the major determinant and the hydrocarbon thickness is a secondary determinant. Cholesterol was also incorporated into DOPC bilayers and X-ray diffuse scattering was used to determine quantitative structure with the result that the area occupied by DOPC in the membrane decreases while bilayer thickness increases in a correlated way because lipid volume does not change. The water permeability decreases with added cholesterol and it correlates in a different way from pure lipids with area per lipid, bilayer thickness, and also with area compressibility.     

Helix tilt of the M2 transmembrane peptide from influenza A virus: an intrinsic property

Solid-state NMR has been used to study the influence of lipid bilayer hydrophobic thickness on the tilt of a peptide (M2-TMP) representing the transmembrane portion of the M2 protein from influenza A. Using anisotropic (15)N chemical shifts as orientational constraints, single-site isotopically labeled M2-TMPs were studied in hydrated dioleoylphosphatidylcholine (DOPC) and dimyristoylphosphatidylcholine (DMPC) lipid bilayers oriented between thin glass plates. These chemical shifts provide orientational information for the molecular frame with respect to the magnetic field in the laboratory frame. When modeled as a uniform ideal alpha-helix, M2-TMP has a tilt of 37(+/-3) degrees in DMPC and 33(+/-3) degrees in DOPC with respect to the bilayer normal in these lipid environments. The difference in helix tilt between the two environments appears to be small. This lack of a substantial change in tilt further suggests that significant interactions occur between the helices, as in an oligomeric state, to prevent a change in tilt in thicker lipid bilayers.     

Effects of cholesterol concentration on the interaction of cytarabine with lipid membranes: a molecular dynamics simulation study

Liposomal cytarabine, DepoCyt, is a chemotherapy agent which is used in cancer treatment. This form of cytarabine has more efficacy and fewer side effects relative to the other forms. Since DepoCyt contains the cytarabine encapsulated within phosphatidylcholine and the sterol molecules, we modeled dioleoylphosphatidylcholine (DOPC)/cholesterol bilayer membrane as a carrier for cytarabine to study drug–bilayer interactions. For this purpose, we performed a series of united-atom molecular dynamics (MD) simulations for 25?ns to investigate the interactions between cytarabine and cholesterol-containing DOPC lipid bilayers. Only the uncharged form of cytarabine molecule was investigated. In this study, different levels of the cholesterol content (0, 20, and 40%) were used. MD simulations allowed us to determine dynamical and structural properties of the bilayer membrane and to estimate the preferred location and orientation of the cytarabine molecule inside the bilayer membrane. Properties such as membrane thickness, area per lipid, diffusion coefficient, mass density, bilayer packing, order parameters, and intermolecular interactions were examined. The results show that by increasing the cholesterol concentration in the lipid bilayers, the bilayer thickness increases and area per lipid decreases. Moreover, in accordance with the experiments, our calculations show that cholesterol molecules have ordering effect on the hydrocarbon acyl chains. Furthermore, the cytarabine molecule preferentially occupies the polar region of the lipid head groups to form specific interactions (hydrogen bonds). Our results fully support the experimental data. Our finding about drug–bilayer interaction is crucial for the liposomal drug design.     

Fluid bilayer structure determination by the combined use of x-ray and neutron diffraction. I. Fluid bilayer models and the limits of resolution.

This is the first in a series of papers concerned with methods for the determination of the structures of fluid phospholipid bilayers in the liquid-crystalline (L alpha) phase. The basic approach is the joint refinement of quasimolecular models (King and White, 1986. Biophys. J. 49:1047-1054) using x-ray and neutron diffraction data. We present here (a) the rationale for quasimolecular models, (b) the nature of the resolution problem for thermally disordered bilayers, and (c) an analysis of the resolution of experiments in which Gaussian functions are used to describe the distribution of submolecular components. We show that multilamellar liquid-crystalline bilayers are best described by the convolution of a perfect lattice function with a thermally disordered bilayer unit cell. Lamellar diffraction measurements on such a system generally yield only 5-10 orders of diffraction data from which transbilayer profiles of the unit cell can be constructed. The canonical resolution of these transbilayer profiles, defined as the Bragg spacing divided by the index of the highest recorded diffraction order, is typically 5-10 A. Using simple model calculations, we show that the canonical resolution is a measure of the widths of the distributions of constituents of the unit cell rather than a measure of the spatial separation of the distributions. The widths provide a measure of the thermal motion of the bilayer constituents which can be described by Gaussian functions. The equilibrium positions of the centers of the distributions can be determined with a precision of 0.1-0.5 A based upon typical experimental errors.