Molecular dynamics simulation of structural changes of lipid bilayers induced by shock waves: Effects of incident angles

Unsteady and nonequilibrium molecular dynamics simulations of the response of dipalmitoylphosphatidylcholine (DPPC) bilayers to the shock waves of various incident angles are presented. The action of an incident shock wave is modeled by adding a momentum in an oblique direction to water molecules adjacent to a bilayer. We thereby elucidate the effects of incident shock angles on (i) collapse and rebound of the bilayer, (ii) lateral displacement of headgroups, (iii) tilts of lipid molecules, (iv) water penetration into the hydrophobic region of the bilayer, and (v) momentum transfer across the bilayer. The number of water molecules delivered into the hydrophobic region is found to be insensitive to incident shock angles. The most important structural changes are the lateral displacement of headgroups and tilts of lipid molecules, which are observed only in the half of the bilayer directly exposed to a shock wave for all incident shock angles studied here. As a result, only the normal component of the added oblique momentum is substantially transferred across the bilayer. This also suggests that the irradiation by shock waves may induce a jet-like streaming of the cytoplasm toward the nucleus.     

Molecular dynamics simulation of structural changes of lipid bilayers induced by shock waves: Effects of incident angles

Unsteady and nonequilibrium molecular dynamics simulations of the response of dipalmitoylphosphatidylcholine (DPPC) bilayers to the shock waves of various incident angles are presented. The action of an incident shock wave is modeled by adding a momentum in an oblique direction to water molecules adjacent to a bilayer. We thereby elucidate the effects of incident shock angles on (i) collapse and rebound of the bilayer, (ii) lateral displacement of headgroups, (iii) tilts of lipid molecules, (iv) water penetration into the hydrophobic region of the bilayer, and (v) momentum transfer across the bilayer. The number of water molecules delivered into the hydrophobic region is found to be insensitive to incident shock angles. The most important structural changes are the lateral displacement of headgroups and tilts of lipid molecules, which are observed only in the half of the bilayer directly exposed to a shock wave for all incident shock angles studied here. As a result, only the normal component of the added oblique momentum is substantially transferred across the bilayer. This also suggests that the irradiation by shock waves may induce a jet-like streaming of the cytoplasm toward the nucleus.     

爆炸波对生物膜微观创伤的分子动力学分析

爆炸冲击波作用到人体胸部时，肺部会出现肺出血及肺水肿等症状，这是人体爆炸创伤的主要原因，深入研究很有必要。为了更好地理解爆炸创伤的机理，应研究冲击波与微观组织作用的力学过程，但具有一定的难度，本文先从基本的生物膜做起。本文运用分子动力学研究了冲击波对DPPC膜造成的损伤，通过停止活塞来控制冲击波的冲量，观察了冲击过程中膜的恢复情况。观察了不同冲量下冲击波经过膜后磷脂分子及其周围水分子分布，发现随着冲量增大，膜越来越无序混乱，褶皱更严重，疏水区水分子越来越多。将膜冲击过程划为三个阶段，分别为冲击阶段、恢复阶段和后效阶段。发现当冲量大于153 mPa s时，在冲击过程中没有观察到膜的损伤恢复。     

Detachment and sonoporation of adherent HeLa-cells by shock wave-induced cavitation

The interaction of lithotripter-generated shock waves with adherent cells is investigated using high-speed optical techniques. We show that shock waves permeabilize adherent cells in vitro through the action of cavitation bubbles. The bubbles are formed in the trailing tensile pulse of a lithotripter-generated shock wave where the pressure drops below the vapor pressure. Upon collapse of cavitation bubbles, a strong flow field is generated which accounts for two effects: first, detachment of cells from the substrate; and second, the temporary opening of cell membranes followed by molecular uptake, a process called sonoporation. Comparison of observed cell detachment with results from a theoretical model considering peeling cell detachment by a wall jet-induced shear stress shows reasonable agreement.     

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

The cooperative behaviour of antimicrobial peptides in model membranes

A systematic analysis of the hypothesis of the antimicrobial peptides' (AMPs) cooperative action is performed by means of full atomistic molecular dynamics simulations accompanied by circular dichroism experiments. Several AMPs from the aurein family (2.5,2.6, 3.1), have a similar sequence in the first ten amino acids, are investigated in different environments including aqueous solution, trifluoroethanol (TFE), palmitoyloleoylphosphatidylethanolamine (POPE), and palmitoyloleoylphosphatidylglycerol (POPG) lipid bilayers. It is found that the cooperative effect is stronger in aqueous solution and weaker in TFE. Moreover, in the presence of membranes, the cooperative effect plays an important role in the peptide/lipid bilayer interaction. The action of AMPs is a competition of the hydrophobic interactions between the side chains of the peptides and the hydrophobic region of lipid molecules, as well as the intra peptide interaction. The aureins 2.5-COOH and 2.6-COOH form a hydrophobic aggregate to minimize the interaction between the hydrophobic group and the water. Once that the peptides reach the water/lipid interface the hydrophobic aggregate becomes smaller and the peptides start to penetrate into the membrane. In contrast, aurein 3.1-COOH forms only a transient aggregate which disintegrates once the peptides reached the membrane, and it shows no cooperativity in membrane penetration.     

Environments and conformations of tryptophan side chains of gramicidin A in phospholipid bilayers studied by Raman spectroscopy.

Raman spectroscopy has been used to investigate the hydrophobic interaction of the indole ring with the environments, the water accessibility to the N1H site, and the conformation about the C beta-C3 bond for the four tryptophan side chains of gramicidin A incorporated into phospholipid bilayers. Most of the tryptophan side chains of the head-to-head helical dimer transmembrane channel are strongly interacting with the lipid hydrocarbon chains, and the hydrophobic interactions for the rest increase with increasing hydrocarbon chain length of the lipid. One tryptophan side chain (probably Trp-15) is accessible to water molecules, another (Trp-9) is deeply buried in the bilayer and inaccessible, and the accessibilities of the remaining two (Trp-11 and Trp-13) depend on the bilayer thickness. The torsional angle about the C beta-C3 bond is found to be +/- 90 degrees for all the tryptophans irrespective of the membrane thickness. Binding of the sodium cation to the channel does not change the torsional angles but decreases the water accessibilities of two tryptophans (Trp-11 and Trp-13) considerably. In conjunction with a slight spectral change in the amide III region, it is suggested that the sodium binding causes a partial change in the main-chain conformation around Trp-11 and Trp-13, which results in the movements of these side chains toward the bilayer center. Two models consistent with the present Raman data are proposed for the tryptophan orientation in the dominant channel structure.     

Hydrophobic Barrier in Liposomes Studied by FT-ESEEM

An improvement of the electron spin echo envelope modulation technique (ESEEM) for studying the hydrophobic barrier in lipid membranes was proposed. Water penetration depth into the lipid bilayer composed of egg-phosphatidylcholine and dicetylphosphate was studied. For this purpose the stearic acid spin probes with nitroxide moiety at different positions on the acyl chain were dissolved in the bilayer of liposomes prepared in buffer solution with D2O. The ESEEM method is based on the deuterium effect to the 14N-hyperfine splitting constant of the spin probe. For the analysis of the spectra we suggested the frequency domain approach (Fourier transform-ESEEM) instead of the time domain. In order to test this approach, the influence of cholesterol and octanol on the water penetration depth into liposomes was measured. FT-ESEEM appeared to be a more sensitive method for the quantitative measurements of the water penetration depth than the time domain. This approach can become one of the most effective methods for studying the influence of surfactants and other biologically active molecules on water penetration depth into membranes.     

Identification of the hydrophobic thickness of a membrane protein using fluorescence spectroscopy: studies with the mechanosensitive channel MscL

The hydrophobic thickness of a membrane protein is an important parameter, defining how the protein sits within the hydrocarbon core of the lipid bilayer that surrounds it in a membrane. Here we show that Trp scanning mutagenesis combined with fluorescence spectroscopy can be used to define the hydrophobic thickness of a membrane protein. The mechanosensitive channel of large conductance (MscL) contains two transmembrane alpha-helices, of which the second (TM2) is lipid-exposed. The region of TM2 that spans the hydrocarbon core of the bilayer when MscL is reconstituted into bilayers of dioleoylphosphatidylcholine runs from Leu-69 to Leu-92, giving a hydrophobic thickness of ca. 25 A. The results obtained using Trp scanning mutagenesis were confirmed using Cys residues labeled with the N-methyl-amino-7-nitroben-2-oxa-1,3-diazole [NBD] group; both fluorescence emission maxima and fluorescence lifetimes for the NBD group are sensitive to solvent dielectric constant over the range (2-40) thought to span the lipid headgroup region of a lipid bilayer. Changing phospholipid fatty acyl chain lengths from C14 and C24 results in no significant change for the fluorescence of the interfacial residues, suggesting very efficient hydrophobic matching between the protein and the surrounding lipid bilayer.     

Molecular dynamics simulation of melittin in a dimyristoylphosphatidylcholine bilayer membrane.

Molecular dynamics trajectories of melittin in an explicit dimyristoyl phosphatidylcholine (DMPC) bilayer are generated to study the details of lipid-protein interactions at the microscopic level. Melittin, a small amphipathic peptide found in bee venom, is known to have a pronounced effect on the lysis of membranes. The peptide is initially set parallel to the membrane-solution interfacial region in an alpha-helical conformation with unprotonated N-terminus. Solid-state nuclear magnetic resonance (NMR) and polarized attenuated total internal reflectance Fourier transform infrared (PATIR-FTIR) properties of melittin are calculated from the trajectory to characterize the orientation of the peptide relative to the bilayer. The residue Lys7 located in the hydrophobic moiety of the helix and residues Lys23, Arg24, Gln25, and Gln26 at the C-terminus hydrophilic form hydrogen bonds with water molecules and with the ester carbonyl groups of the lipids, suggesting their important contribution to the stability of the helix in the bilayer. Lipid acyl chains are closely packed around melittin, contributing to the stable association with the membrane. Calculated density profiles and order parameters of the lipid acyl chains averaged over the molecular dynamics trajectory indicate that melittin has effects on both layers of the membrane. The presence of melittin in the upper layer causes a local thinning of the bilayer that favors the penetration of water through the lower layer. The energetic factors involved in the association of melittin at the membrane surface are characterized using an implicit mean-field model in which the membrane and the surrounding solvent are represented as structureless continuum dielectric material. The results obtained by solving the Poisson-Bolztmann equation numerically are in qualitative agreement with the detailed dynamics. The influence of the protonation state of the N-terminus of melittin is examined. After 600 ps, the N-terminus of melittin is protonated and the trajectory is continued for 400 ps, which leads to an important penetration of water molecules into the bilayer. These observations provide insights into how melittin interacts with membranes and the mechanism by which it enhances their lysis.     

Profile of changes in lipid bilayer structure caused by beta-amyloid peptide.

beta-Amyloid peptide (A beta) is the primary constituent of senile plaques, a defining feature of Alzheimer's disease. Aggregated A beta is toxic to neurons, but the mechanism of toxicity is uncertain. One hypothesis is that interactions between A beta aggregates and cell membranes mediate A beta toxicity. Previously, we described a positive correlation between the A beta aggregation state and surface hydrophobicity, and the ability of the peptide to decrease fluidity in the center of the membrane bilayer [Kremer, J. J., et al. (2000) Biochemistry 39, 10309--10318]. In this work, we report that A beta aggregates increased the steady-state anisotropy of 1,6-diphenyl-1,3,5-hexatriene (DPH) embedded in the hydrophobic center of the membrane in phospholipids with anionic, cationic, and zwitterionic headgroups, suggesting that specific charge--charge interactions are not required for A beta--membrane interactions. A beta did not affect the fluorescence lifetime of DPH, indicating that the increase in anisotropy is due to increased ordering of the phospholipid acyl chains rather than changes in water penetration into the bilayer interior. A beta aggregates affected membrane fluidity above, but not below, the lipid phase-transition temperature and did not alter the temperature or enthalpy of the phospholipid phase transition. A beta induced little to no change in membrane structure or water penetration near the bilayer surface. Overall, these results suggest that exposed hydrophobic patches on the A beta aggregates interact with the hydrophobic core of the lipid bilayer, leading to a reduction in membrane fluidity. Decreases in membrane fluidity could hamper functioning of cell surface receptors and ion channel proteins; such decreases have been associated with cellular toxicity.     

Perturbation of the hydrophobic core of lipid bilayers by the human antimicrobial peptide LL-37

LL-37 is a cationic, amphipathic alpha-helical antimicrobial peptide found in humans that kills cells by disrupting the cell membrane. To disrupt membranes, antimicrobial peptides such as LL-37 must alter the hydrophobic core of the bilayer. Differential scanning calorimetry and deuterium ((2)H) NMR experiments on acyl chain perdeuterated lipids demonstrate that LL-37 inserts into the hydrophobic region of the bilayer and alters the chain packing and cooperativity. The results show that hydrophobic interactions between LL-37 and the hydrophobic acyl chains are as important for the ability of this peptide to disrupt lipid bilayers as its electrostatic interactions with the polar headgroups. The (2)H NMR data are consistent with the previously determined surface orientation of LL-37 (Henzler Wildman, K. A., et al. (2003) Biochemistry 42, 6545) with an estimated 5-6 A depth of penetration of the hydrophobic face of the amphipathic helix into the hydrophobic interior of the bilayer. LL-37 also alters the material properties of lipid bilayers, including the area per lipid, hydrophobic thickness, and coefficient of thermal expansion in a manner that varies with lipid type and temperature. Comparison of the effect of LL-37 on 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC-d(31)) and 1,2-dimyristoyl-phosphatidylcholine (DMPC-d(54)) at different temperatures demonstrates the importance of bilayer order in determining the type and extent of disordering and disruption of the hydrophobic core by LL-37. One possible explanation, which accounts for both the (2)H NMR data presented here and the known surface orientation of LL-37 under identical conditions, is that bilayer order influences the depth of insertion of LL-37 into the hydrophobic/hydrophilic interface of the bilayer, altering the balance of electrostatic and hydrophobic interactions between the peptide and the lipids.     

Simulating induced interdigitation in membranes

In this study we introduce a mesoscopic lipid-water-alcohol model. Dissipative particle dynamics (DPD) simulations have been used to investigate the induced interdigitation of bilayers consisting of double-tail lipids by adding alcohol molecules to the bilayer. Our simulations nicely reproduce the experimental phase diagrams. We find that alcohol can induce an interdigitated structure where the common bilayer structure changes into monolayer in which the alcohol molecules screen the hydrophobic tails from the water phase. At low concentrations of alcohol the membrane has domains of the interdigitated phase that are in coexistence with the common membrane phase. We compute the effect of the chain length of the alcohol on the phase behavior of the membrane and show that the stability of the interdigitated phase depends on the length of the alcohol. We show that we can reproduce the experimental hydrophobic thickness of the bilayer for various combinations of lipids and alcohols. We use our model to clarify some of the experimental questions related to the structure of the interdigitated phase and put forward a simple model that explains the alcohol chain length dependence of the stability of this interdigitated phase.     

Interaction of substance P with phospholipid bilayers: A neutron diffraction study.

Neutron diffraction has been used to study the membrane-bound structure of substance P (SP), a member of the tachykinin family of neuropeptides. The depth of penetration of its C-terminus in zwitterionic and anionic phospholipid bilayers was probed by specific deuteration of leucine 10, the penultimate amino acid residue. The results show that the interaction of SP with bilayers, composed of either dioleoylphosphatidylcholine (DOPC), or a 50:50 mixture of DOPC and the anionic phospholipid dioleoylphosphatidylglycerol (DOPG), takes place at two locations. One requires insertion of the peptide into the hydrophobic region of the bilayer, the other is much more peripheral. The penetration of the peptide into the hydrophobic region of the bilayer is reflected in a marked difference in the water distribution profiles. SP is seen to insert into DOPC bilayers, but a larger proportion of the peptide is found at the surface when compared to the anionic bilayers. The positions of the two label populations show only minor differences between the two types of bilayer.     

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells’ membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane’s mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.     

Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study

All-atom molecular dynamics simulations have been performed on cimetidine in the presence of a palmitoyloleoylphosphatidylcholine (POPC) bilayer. The free energy profile of a single cimetidine molecule passing across POPC bilayer displays a minimum at the interface of bilayer and water. Ten cimetidine molecules were inserted into POPC bilayer to obtain an 8 mol % drug model, and molecular dynamics results showed that cimetidine molecules reside at the polar region of POPC bilayer with sulphur atoms directing to the hydrophobic region. By comparing the one drug model with 8 mol % drug model, one can see that the central barrier to cross the membrane increases while the free energy in bulk water decreases, indicating that the ability of cimetidine passing across the POPC bilayer weakens at increased concentration. In addition, the free energy minimum shifts closer to the hydrophobic core. Our results indicate that with the increased drug concentration, it is more difficult for cimetidine to enter and pass across POPC bilayer.     

Predicting doxorubicin drug delivery by single-walled carbon nanotube through cell membrane in the absence and presence of nicotine molecules: a molecular dynamics simulation study

Abstract

In order to study the interaction of the anticancer agent Doxorubicin with the single-walled carbon nanotubes with different diameters as drug delivery systems, the molecular dynamics (MD) simulations have been used. Also, for design and development of intracellular Doxorubicin drug delivery systems, a series of steered MD simulations are applied to explore the possibility of encapsulated Doxorubicin–carbon nanotube penetration through a lipid bilayer in presence and absence of Nicotine molecules at different pulling rates. Our simulation results showed that in spite of the adsorption of drug molecules on the outer sidewall of the nanotubes, the spontaneous localization of one Doxorubicin molecule into the cavity of the nanovectors with larger diameters is observed. It is found that the presence of Nicotine molecules in extracellular medium increases the required force for pulling nanotube-encapsulated drug as well as the required time for penetration process, especially at higher velocity. Also, the entering process of the Nicotine molecules into the carbon nanotube causes that the encapsulated drug molecule is fully released in the hydrophobic phase of the lipid bilayer.

Communicated by Ramaswamy H. Sarma     

Phosphatidylinositol induces fluid phase formation and packing defects in phosphatidylcholine model membranes

Liposomes consisted of phosphatidylinositol (PI) and phosphatidylcholine (PC) have been utilized as delivery vehicle for drugs and proteins. In the present work, we studied the effect of soy PI on physical properties of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) liposomes such as phase state of lipid bilayer, lipid packing and phase properties using multiple orthogonal biophysical techniques. The 6-dodecanoyl-2-dimethylamino naphthalene (Laurdan) fluorescence studies showed that presence of PI induces the formation of fluid phases in DMPC. Differential scanning calorimetry (DSC), temperature dependent fluorescence anisotropy measurements, and generalized polarization values for Laurdan showed that the presence of as low as 10mol% of PI induces substantial broadening and shift to lower temperature of phase transition of DMPC. The fluorescence emission intensity of DPH labeled, PI containing DMPC lipid bilayer decreased possibly due to deeper penetration of water molecules in lipid bilayer. In order to further delineate the effect of PI on the physico chemical properties of DMPC is due to either significant hydrophobic mismatch between the acyl chains of the DMPC and that of soy PI or due to the inositol head group, we systematically replaced soy PI with PC species of similar acyl chain composition (DPPC and 18:2 (Cis) PC) or with diacylglycerol (DAG), respectively. The anisotropy of PC membrane containing soy PI showed largest fluidity change compared to other compositions. The data suggests that addition of PI alters structure and dynamics of DMPC bilayer in that it promotes deeper water penetration in the bilayer, induces fluid phase characteristics and causes lipid packing defects that involve its inositol head group.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

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.