Antagonism of membrane compression effects by high pressure gas mixtures in a phospholipid bilayer system

Binary mixtures of helium with nitrogen, xenon or nitrous oxide were applied to suspensions of phosphatidylcholine-cholesterol vesicles to determine those mixtures of lipid soluble gases which would exactly antagonize the membrane rigidifying effect of 100 ATA compression. A previous study has shown that the initial application of 100 ATA compression by gas produces a significant reduction in the fluidity of the phospholipid bilayer. However, as the high pressure gas dissolves into the lipid region it creates disorder and increases fluidity. Fluidity of the bilayer at equilibrium represents the sum of the compression-ordering and dissolved-gas disordering effects and is dependent on the gas/lipid partition coefficient of the particular gas. The beneficial effect of a narcotic gas added to Trimix mixtures to ameliorate HPNS in deep divers may be due to a balance of compression-ordering and solubility-disordering effects achieved within the nerve membrane. It is therefore valuable to determine those gas mixtures which achieve balance of these two effects and result in zero net change in phospholipid bilayer fluidity at an established pressure of 100 ATA. Binary mixtures of helium with 88% nitrogen, 3.8% xenon or 2.8% nitrous oxide resulted in zero net change in bilayer fluidity with our model system at 100 ATA. A graph of the percent of narcotic gas needed to produce zero net effect as a function of pressure, however, was nonlinear. This would suggest the ratio of gases in Trimix must be varied as a function of pressure. While the phosphatidylcholine-cholesterol bilayer is a good model for certain components of the nerve membrane, it does not allow for study of protein-lipid or gas-protein interactions. The data presented thus aid in our understanding of HPNS but are yet incomplete for precise use in predicting diving mixtures.     

Effects of small nonpolar molecules on membrane compressibility and permeability: A theoretical study of the effects of anesthetic gases

We explore from a theoretical perspective the effects of small nonpolar molecules, such as anesthetic gases, on membrane compressibility and permeability. As a model system we expand a previously proposed generalization of Nagle's model for biomembrane phase transitions. In this model anesthetic gases alter membrane compressibility, causing profound changes in membrane permeability. Anesthetics either increase or decrease membrane permeability, depending on whether the membrane lipid is originally in the solid or melted state, or in a two-phase region. These changes are reversed by high pressure, in agreement with experimental results. Anesthetic-induced changes in compressibility are predicted to inhibit fusion of phospholipid vesicles to each other and to planar bilayers, and thus might be expected to inhibit the fusion of presynaptic vesicles with the presynaptic nerve membrane. This work provides a detailed molecular theory for many of the effects of anesthetic gases on both synapse and axon, and provides a coherent framework for understanding diverse experimental results.     

Pressure-induced exclusion of a local anesthetic from model and nerve membranes

Because it is well established that the anesthetic state can be reversed by pressure, a number of molecular theories that have been proposed for the mechanism of action of both local and general anesthetics can be tested by varying the pressure. Using Fourier transform infrared spectroscopy, we report here the first direct observation of the expulsion from lipid bilayers of a local anesthetic, tetracaine, by pressure. Moreover, we establish for the first time that this phenomenon is common to both model membranes and to myelinated and unmyelinated nerve membranes, vindicating the utility of model membrane systems. A distinctive feature of this behavior in model systems is that, in saturated phosphatidylcholines at high pH, expulsion only occurs in the presence of cholesterol, whose ordering effect on the acyl chains evidently assists pressure in squeezing the anesthetic out of the bilayer. This pressure-induced phenomenon may provide insight into the molecular mechanisms underlying the antagonistic effect of pressure against anesthesia.     

High-pressure infrared study of phosphatidylserine bilayers and their interactions with the local anesthetic tetracaine

High-pressure Fourier-transform infrared (FT-IR) spectroscopy was used to study the barotropic behavior of phosphatidylserine bilayers and their interactions with the local anesthetic tetracaine. The model membrane systems studied were multilamellar aqueous dispersions of 1,2-dimyristoyl-sn-glycero-3-phospho-L-serine (DMPS) and 1,2-dioleoyl-sn-glycero-3-phospho-L-serine (DOPS) in the absence and the presence of tetracaine at pH 5.5 and 9.5. The infrared spectra were measured at 28 degrees C in a diamond anvil cell as a function of pressure up to 25 kbar. The results show that the barotropic behavior of the negatively charged phosphatidylserine bilayers is very similar to that observed for zwitterionic phospholipids, such as phosphatidylcholine and phosphatidylethanolamine, with corresponding acyl chains. The results also indicate that the local anesthetic partitions into phosphatidylserine bilayers in an environment close to the membrane-water interface and interacts electrostatically with the lipid head group. Application of high hydrostatic pressure on the lipid-anesthetic systems results in the pressure-induced expulsion of the anesthetic from a membrane to an aqueous environment. The pressures required for expulsion of anesthetic from bilayers are much higher for the unsaturated lipid (DOPS) than for the saturated lipid (DMPS) (approximately 6 kbar vs approximately 2 kbar, respectively). Whereas incorporation of the anesthetic into DOPS bilayers does not affect significantly the structural and dynamic properties of the disordered acyl chains in the liquid-crystalline phase, it orders the DMPS acyl chains in the gel phase.     

Solute partitioning into lipid bilayer membranes

We have measured the membrane/water partition coefficients of benzene into lipid bilayers as a function of the surface density of the phospholipid chains. A simple 2H NMR method was used for the measurement of surface densities; it is shown to give results similar to those obtained from more demanding X-ray diffraction measurements. We observe that benzene partitioning into the bilayer is dependent not only on the partitioning chemistry, characterized by the oil/water partition coefficient, but also on the surface density of the bilayer chains. Increasing surface density leads to solute exclusion: benzene partitioning decreases by an order of magnitude as the surface density increases from 50% to 90% of its maximum value, a range readily accessible in bilayers and biomembranes under physiological conditions. This effect is independent of the nature of the agent used to alter surface density: temperature, cholesterol, and phospholipid chain length were tested here. These observations support the recent statistical thermodynamic theory of solute partitioning into chain molecule interphases, which predicts that the expulsion of solute is due to entropic effects of the orientational ordering among the phospholipid chains. We conclude that the partitioning of solutes into bilayer membranes, which are interfacial phases, is of a fundamentally different nature than partitioning into bulk oil and octanol phases.     

Multiple binding sites for local anesthetics in membranes: characterization of the sites and their equilibria by deuterium NMR of specifically deuterated procaine and tetracaine

Anesthetics bound to model membranes were observed directly by means of deuterium nuclear magnetic resonance (NMR). The specifically deuterated local anesthetics procaine and tetracaine were synthesized, and their partition coefficients (water:phosphatidylcholine) and pKa values determined. The interaction of these anesthetics with lamellar dispersions of egg phosphatidylcholine was studied by 2H nuclear magnetic resonance and by electron spin resonance (ESR) of a spin-labelled phospholipid at low (5.5) and high (9.5) pH. The ESR experiments suggest that tetracaine intercalates in the membrane and that it equilibrates between water and the phospholipid bilayers of the multilamellar system. The NMR results are consistent with a model where the anesthetic is (1) free in water, (2) weakly bound, and (3) strongly bound to the membrane. A fast exchange exists between the two first sites, but exchange is slow with the third site. Binding of type 3 is observed only at high pH for procaine, whereas it is found both at low and high pH for tetracaine. Calculations of the partition coefficients for the charged and uncharged forms of tetracaine indicate that both sites, 2 and 3, are occupied by the charged form at low pH and by the uncharged form at high pH. The partition coefficient for the weakly bound species was estimated from an analysis of the dependence of line width on the lipid to water ratio. The NMR data suggest that the binding sites for the strongly bound charged and uncharged species are different, the former probably being closer to the membrane-water interface. Estimates of molecular order parameters for the strongly bound species indicate that it is located with its long molecular axis approximately parallel to the director for ordering of the fatty acyl chains. A small increase in lipid ordering by tetracaine is observed at low pH, as evidenced by 2H NMR of the deuterated N-methyl groups of phosphatidylcholine; the reverse occurs at high pH.     

Interactions of the local anesthetic tetracaine with membranes containing phosphatidylcholine and cholesterol: a 2H NMR study

The interactions of the local anesthetic tetracaine with multilamellar dispersions of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and cholesterol have been investigated by deuterium nuclear magnetic resonance of specifically deuteriated tetracaines, DMPC and cholesterol. Experiments were performed at pH 5.5, when the anesthetic is primarily charged, and at pH 9.5, when it is primarily uncharged. The partition coefficients of the anesthetic in the membrane have been measured at both pH values for phosphatidylcholine bilayers with and without cholesterol. The higher partition coefficients obtained at pH 9.5 reflect the hydrophobic interactions between the uncharged form of the anesthetic and the hydrocarbon region of the bilayer. The lower partition coefficients for the DMPC/cholesterol system at both pH values suggest that cholesterol, which increases the order of the lipid chains, decreases the solubility of tetracaine into the bilayer. For phosphatidylcholine bilayers, it has been proposed [Boulanger, Y., Schreier, S., & Smith, I. C. P. (1981) Biochemistry 20, 6824-6830] that the charged tetracaine at low pH is located mostly at the phospholipid headgroup level while the uncharged tetracaine intercalates more deeply into the bilayer. The present study suggests that the location of tetracaine in the cholesterol-containing system is different from that in pure phosphatidylcholine bilayers: the anesthetic sits higher in the membrane. An increase in temperature results in a deeper penetration of the anesthetic into the bilayer. Moreover, the incorporation of the anesthetic into DMPC bilayers with or without cholesterol results in a reduction of the lipid order parameters both in the plateau and in the tail regions of the acyl chains, this effect being greater with the charged form of the anesthetic.     

HgCl2 disrupts the structure of the human erythrocyte membrane and model phospholipid bilayers

The structural effects of Hg(II) ions on the erythrocyte membrane were studied through the interactions of HgCl2 with human erythrocytes and their isolated resealed membranes. Studies were carried out by scanning electron microscopy and fluorescence spectroscopy, respectively. Hg(II) induced shape changes in erythrocytes, which took the form of echinocytes and stomatocytes. This finding means that Hg(II) locates in both the outer and inner monolayers of the erythrocyte membrane. Fluorescence spectroscopy results indicate strong interactions of Hg(II) ions with phospholipid amino groups, which also affected the packing of the lipid acyl chains at the deep hydrophobic core of the membrane. HgCl2 also interacted with bilayers of dimyristoylphosphatidylcholine and dimyristoylphosphatidylethanolamine, representative of phospholipid classes located in the outer and inner monolayers of the erythrocyte membrane, respectively. X-ray diffraction indicated that Hg(II) ions induced molecular disorder to both phospholipid bilayers, while fluorescence spectroscopy of dimyristoylphosphatidylcholine large unilamellar vesicles confirmed the interaction of Hg(II) ions with the lipid polar head groups. All these findings point to the important role of the phospholipid bilayers in the interaction of Hg(II) on cell membranes.     

Pressure effects on dipalmitoylphosphatidylcholine bilayers measured by 2H nuclear magnetic resonance

The effects of pressure, up to 5 kbar, on multilamellar vesicles of 1,2-dipalmitoyl-sn-phosphatidylcholine perdeuterated in the acyl chains (DPPC-d62) were examined by using high-pressure NMR techniques. A deuterium probe was built, and the quadrupole splitting was measured against pressure at various temperatures. The experiments were performed on pure lipid bilayers in the liquid-crystalline state and on bilayers in the liquid-crystalline state containing the local anesthetic tetracaine. The results show that the order parameter of all segments of the acyl chains increases with pressure in the liquid-crystalline state. The more highly ordered regions of the chains are affected slightly more than the regions near the methyl ends. The addition of tetracaine increases the disorder of the chains, and pressure reverses the effect of anesthetic on the lipid as seen by the reversal of the changes in line shape and the measured order parameter.     

Insight into the putative specific interactions between cholesterol, sphingomyelin, and palmitoyl-oleoyl phosphatidylcholine

The effects of cholesterol (Chol) on phospholipid bilayers include ordering of the fatty acyl chains, condensing of the lipids in the bilayer plane, and promotion of the liquid-ordered phase. These effects depend on the type of phospholipids in the bilayer and are determined by the nature of the underlying molecular interactions. As for Chol, it has been shown to interact more favorably with sphingomyelin than with most phosphatidylcholines, which in given circumstances leads to formation of lateral domains. However, the exact origin and nature of Chol-phospholipid interactions have recently been subjects of speculation. We examine interactions between Chol, palmitoylsphingomyelin (PSM) and palmitoyl-oleoyl-phosphatidylcholine (POPC) in hydrated lipid bilayers by extensive atom-scale molecular dynamics simulations. We employ a tailored lipid configuration: Individual PSM and Chol monomers, as well as PSM-Chol dimers, are embedded in a POPC lipid bilayer in the liquid crystalline phase. Such a setup allows direct comparison of dimeric and monomeric PSMs and Chol, which ultimately shows how the small differences in PSM and POPC structure can lead to profoundly different interactions with Chol. Our analysis shows that direct hydrogen bonding between PSM and Chol does not provide an adequate explanation for their putative specific interaction. Rather, a combination of charge-pairing, hydrophobic, and van der Waals interactions leads to a lower tilt in PSM neighboring Chol than in Chol with only POPC neighbors. This implies improved Chol-induced ordering of PSM's chains over POPC's chains. These findings are discussed in the context of the hydrophobic mismatch concept suggested recently.     

Effects of the local anesthetic tetracaine on the structural and dynamic properties of lipids in model membranes: a high-pressure Fourier transform infrared study

High-pressure Fourier transform infrared (FT-IR) spectroscopy was used to study the effects of a local anesthetic, tetracaine, on the structural and dynamic properties of lipids in model membranes. The model membrane systems studied were multilamellar aqueous dispersions of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-di-O-hexadecyl-sn-glycero-3-phosphocholine (DHPC) in the absence and presence of a physiological concentration of cholesterol (30 mol %). The infrared spectra were measured at 28 degrees C in a diamond anvil cell as a function of pressure up to 25 kbar. The results indicate that the effects of tetracaine on the structure of pure DMPC bilayers in the gel state are dependent on the state of charge of the anesthetic. The uncharged tetracaine disorders the lipid acyl chains while the charged form induces the formation of an interdigitated gel phase. The presence of cholesterol in the latter system prevents the formation of the interdigitated phase, whereas in the former system it disorders the lipid acyl chains in the gel state. Moreover, it is shown that the addition of uncharged tetracaine to interdigitated DHPC bilayers does not alter the interdigitated state of the hydrocarbon chains.     

Elastic response of bilayers with intrinsic proteins

Intrinsic membrane proteins affect the ordering of neighbouring lipid chains. We have used a model of protein-lipid interactions in bilayers proposed by Owicki et al. (Owicki, J.C., Springgate, M.W. and McConnell, H.M. (1978) Proc. Natl. Acad. Sci. U.S.A. 75, 1616–1619) to show that near the lipid phase transition this effect may significantly increase the magnitude of a membrane's lateral compressibility (or correspondingly, decrease the magnitude of the membrane's elastic moduli).     

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

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

Structural aspects of pressure effects on infrared spectra of mixed-chain phosphatidylcholine assemblies in D2O

The barotropic behavior of D2O dispersions of 1-stearoyl-2-caproyl-sn-glycero-3-phosphocholine, C(18):C(10)PC, a highly asymmetric phospholipid in which the length of the fully extended acyl chain at the sn-1 position of the glycerol backbone is twice as long as that at the sn-2 position, has been investigated by high-pressure Fourier transform infrared spectroscopy. This asymmetric phosphatidylcholine bilayer at room temperature displays a pressure-induced phase transition corresponding to the liquid-crystalline----gel phase transition at 1.4 kbar. A conformational ordering of the lipid acyl chains is observed to take place abruptly at the transition pressure of 1.4 kbar. However, the lamellar lipid molecules and their acyl chains remain to be orientationally disordered in the gel phase until the applied pressure reaches 5.5 kbar. In the gel phase of fully hydrated C(18):C(10)PC, the asymmetric lipid molecules assemble into mixed interdigitated bilayers with perpendicular orientation of the zigzag planes among neighboring acyl chains. The role of excess water played in the interchain structure and the behavior of excess water and bound water under high pressure are also discussed.     

The effect of isoflurane on erythrocyte membranes studied by ATR-FTIR

The effect of isoflurane on erythrocyte membranes has been investigated by means of attenuated total reflection infrared spectroscopy. Infrared spectra were measured on sonicated erythrocyte ghosts layered upon a ZnSe crystal covered with D(2)O saline solutions containing increasing amounts of isoflurane. At clinically relevant anesthetic concentrations and 37 degrees C, significant changes in the structural and dynamic properties of the membrane phospholipid bilayers are observed. Both the acyl chain methylene symmetric and asymmetric stretching modes and the carbonyl ester stretching band displayed frequency shifts interpreted as transitions toward disordered liquid-like structure accompanied by dehydration of the phospholipid polar heads. In turn, no secondary structure-linked changes are observed in the amide I region of membrane proteins. Higher anesthetic concentrations (500-900 microM), resulted in progressive detachment of the multilayers from the ATR crystal and irreversible formation of denatured protein. Polarization studies in correspondence of the acyl lipid methylene stretching bands indicated that isoflurane decreases the dichroic ratio thus inducing disorder in the orientation of the lipid acyl chains.     

Effects of hydrostatic pressure on the location of PRODAN in lipid bilayers: a FT-IR study

The effects of hydrostatic pressure on the location of 6-propionyl-2-(dimethylamino)naphthalene (PRODAN), an environmentally sensitive fluorescent probe, in phosphatidylcholine lipid bilayers have been studied by Fourier-transform infrared spectroscopy (FT-IR) over the pressure range of 0.001-25 kbar. The results derived from the PRODAN C = O stretching band, the correlation field splitting of the methylene scissoring mode, and the methylene symmetric stretching mode as well as the absorption of the naphthalene ring show that in the sample of 4% (w/w) PRODAN in dimyristoyl-L-alpha-phosphatidylcholine (DMPC) at pH 6.8, most of the PRODAN molecules are embedded in the bilayers. In contrast, at pH 3.0, PRODAN was found to reside either on the membrane surface or dispersed in water. Compared to DMPC, egg yolk phosphatidylcholine (egg PC), which contains a substantial amount of unsaturated fatty acyl chains, is more susceptible to PRODAN permeation. The present study shows that the pressure dependence of the location of PRODAN in lipid membranes is different from that of tetracaine, a local anesthetic, in lipid bilayers. The model regarding the PRODAN location in lipid bilayers derived from the present infrared data has been compared with that obtained with previous fluorescence studies.     

Factors affecting molecular packing in mixed lipid monolayers and bilayers

The nature of the molecular interactions and the factors determining molecular packing in mixed phospholipid/glyceride monolayers and bilayers were investigated by monolayer and nuclear magnetic resonance (NMR) techniques. Force-area curves were obtained at various temperatures for monolayers, at the air-water interface, of synthetic lecithins and a phosphatidyl ethanolamine mixed with di- and triglycerides in different molar ratios. The linewidths of peaks in the high resolution NMR spectra of lecithin/glyceride co-dispersions in excess water at different temperatures were used to obtain information about molecular mobilities.It was found that the molecular packing in mixed lipid monolayers and bilayers is determined by the following factors: (1) Whether lipid chains are above or below their melting point (T_C). (2) The difference between experimental temperature and T_C: the larger the difference, the smaller the effect of one component on the other. (3) The degree of similarity of the chains of the components; this influences the degree of cooperativity of chain motions and the degree of mixing of the components. (4) The nature, orientation, mutual interaction and degree of hydration of the polar groups.It is shown that mean molecular area does not always reflect the state of chain motions in lipid films, because of heterogeneity of motion and structure along the molecules. Cooperativity of motion may reduce steric requirements; other effects which are of particular importance for lecithins are interactions of zwitterions, and the influence of polar group hydration.     

Cholesterol dynamics in membranes.

Time-resolved fluorescence anisotropy of the sterol analogue, cholestatrienol, and 13C nuclear magnetic resonance (NMR) spin lattice relaxation time (T1c) measurements of [13C4] labeled cholesterol were exploited to determine the correlation times characterizing the major modes of motion of cholesterol in unsonicated phospholipid multilamellar liposomes. Two modes of motion were found to be important: (a) rotational diffusion and (b) time dependence of the orientation of the director for axial diffusion, or "wobble." From the time-resolved fluorescence anisotropy decays of cholestatrienol in egg phosphatidylcholine (PC) bilayers, a value for tau perpendicular, the correlation time for wobble, of 0.9 x 10(-9) s and a value for S perpendicular, the order parameter characterizing the same motion, of 0.45 s were calculated. Both tau perpendicular and S perpendicular were relatively insensitive to temperature and cholesterol content of the membranes. The T1c measurements of [13C4] labeled cholesterol did not provide a quantitative determination of tau parallel, the correlation time for axial diffusion. T1c from the lipid hydrocarbon chains suggested a value for tau perpendicular similar to that for cholesterol. Steady-state anisotropy measurements and time-resolved anisotropy measurements of cholestatrienol were used to probe sterol behavior in a variety of pure and mixed lipid multilamellar liposomes. Both the lipid headgroups and the lipid hydrocarbons chains contributed to the determination of the sterol environment in the membrane, as revealed by these fluorescence measurements. In particular, effects of the phosphatidylethanolamine (PE) headgroup and of multiple unsaturation in the lipid hydrocarbon chains were observed. However, while the steady-state anisotropy was sensitive to these factors, the time-resolved fluorescence analysis indicated that tau perpendicular was not strongly affected by the lipid composition of the membrane. S perpendicular may be increased by the presence of PE. Both steady-state anisotropy measurements and time-resolved anisotropy measurements of cholestatrienol were used to probe sterol behavior in three biological membranes: bovine rod outer segment (ROS) disk membranes, human erythrocyte plasma membranes, and light rabbit muscle sarcoplasmic reticulum membranes. In the ROS disk membranes the value for S perpendicular was marginally higher than in the PC membranes, perhaps reflecting the influence of PE. The dramatic difference noted was in the value for tau perpendicular. In both the ROS disk membranes and the erythrocyte membranes, tau perpendicular was one-third to one-fifth of tau perpendicular in the phospholipid bilayers. This result may reveal an influence of membrane proteins on sterol behavior.     

Exploring the effect of xenon on biomembranes

We report the initial findings of 100 ns molecular dynamics simulations of the role of cellular membranes in general anaesthesia. The effect of xenon on hydrated dipalmitoylphosphatidylcholine bilayers is described. The xenon atoms were found to prefer the interfacial and central regions of the bilayer. The presence of xenon was observed to lead to a small increase in the surface area, membrane thickness, and order of the acyl chains.     

The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.