Cholesterol effects on the phosphatidylcholine bilayer polar region: a molecular simulation study

A molecular dynamics (MD) simulation of a fully hydrated, liquid-crystalline dimyristoylphosphatidylcholine (DMPC)-Chol bilayer membrane containing approximately 22 mol% Chol was carried out for 4.3 ns. The bilayer reached thermal equilibrium after 2.3 ns of MD simulation. A 2.0-ns trajectory generated during 2.3-4.3 ns of MD simulation was used for analyses to determine the effects of Chol on the membrane/water interfacial region. In this region, 70% of Chol molecules are linked to DMPC molecules via short-distance interactions, where the Chol hydroxyl group (OH-Chol) is 1) charge paired to methyl groups of the DMPC choline moiety ( approximately 34%), via the hydroxyl oxygen atom (Och); 2) water bridged to carbonyl ( approximately 19%) and nonester phosphate ( approximately 14%) oxygen atoms, via both Och and the hydroxyl hydrogen atom (Hch); and 3) directly hydrogen (H) bonded to carbonyl ( approximately 11%) and nonester phosphate ( approximately 5%) oxygen atoms, via Hch ( approximately 17% of DMPC-Chol links are multiple). DMPC's gamma-chain carbonyl oxygen atom is involved in 44% of water bridges and 51% of direct H bonds formed between DMPC and Chol. On average, a Chol molecule forms 0.9 links with DMPC molecules, while a DMPC molecule forms 2.2 and 0.3 links with DMPC and Chol molecules, respectively. OH-Chol makes hydrogen bonds with 1.1 water molecules, preferentially via Hch. The average number of water molecules H bonded to the DMPC headgroup is increased by 7% in the presence of Chol. These results indicate that inclusion of Chol decreases interlipid links and increases hydration in the polar region of the membrane.     

Cholesterol effects on a mixed-chain phosphatidylcholine bilayer: a molecular dynamics simulation study

A molecular dynamics simulation of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer containing approximately 22 mol% of cholesterol (POPC-Chol) was carried out for 15 ns. An 8-ns trajectory was analysed to determine the effects of Chol on the membrane properties and compare it with that on the fully saturated 1,2-dimyristoyl-phosphatidylcholine bilayer containing approximately 22 mol% of Chol (DMPC-Chol). The study suggests that the experimentally observed weaker effect of Chol on the POPC than DMPC bilayer might result from a different vertical localisation of the Chol hydroxyl group (OH-Chol) in both bilayers: in the POPC-Chol bilayer, OH-Chol is placed approximately 3 A higher in the bilayer interface than in the DMPC-Chol bilayer. Because of the rigid cis double bond in the beta-chain of POPC, Chol fits worse to the POPC-Chol membrane environment and is pushed up, in effect all Chol ring atoms are, on average, located above the double bond. Both in mono-cis-unsaturated and fully saturated PC bilayers, Chol induces stronger van der Waals interactions among the chains, whereas its interactions with the chains are weak. In contrast to DMPC, the smooth alpha-face of the Chol ring lowers the order of POPC chains, whereas the rough beta-face increases the order.     

Effects of phospholipid unsaturation on the membrane/water interface: a molecular simulation study.

Molecular dynamics (MD) simulations of fully hydrated bilayers in the liquid-crystalline state made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) or 1-palmitoyl-2-elaidoyl-phosphatidylcholine (PEPC) were carried out to investigate the effect of the incorporation of a double bond in the phosphatidylcholine (PC) beta-chain (cis or trans) on the membrane/water interface. The bilayers reached thermal equilibrium after 3 and 1 ns of MD simulations, respectively, and productive runs were carried out for 3 ns for each bilayer. As reference systems, the 1,2-dimyristoyl-phosphatidylcholine (DMPC) bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) and DMPC-cholesterol (Chol) bilayer containing 22 mol % Chol (M. Pasenkiewicz-Gierula, T. Róg, K. Kitamura, A. and Kusumi, 2000, Biophys. J. 78:1376-1389) were used. The study shows that at the interface of POPC, PEPC, and DMPC-Chol bilayers, average numbers of PC-water and PC-PC interactions are similar and, respectively, greater and smaller than in the DMPC bilayer. The average area/PC in mono-unsaturated bilayers is approximately 4 A(2) larger than in the DMPC bilayer; nevertheless, a strong correlation was found between a single molecular area (SMA) of a PC and the number of interactions this PC makes; i.e., PCs (either saturated or unsaturated) with the same SMA form similar numbers of intermolecular links. The numbers and corresponding SMAs are distributed about averages pertinent to each bilayer. No significant difference between cis and trans bonds was found.     

Molecular dynamics simulation studies of lipid bilayer systems

The main structural element of biological membranes is a liquid-crystalline lipid bilayer. Other constituents, i.e. proteins, sterols and peptides, either intercalate into or loosely attach to the bilayer. We applied a molecular dynamics simulation method to study membrane systems at various levels of compositional complexity. The studies were started from simple lipid bilayers containing a single type phosphatidylcholine (PC) and water molecules (PC bilayers). As a next step, cholesterol (Chol) molecules were introduced to the PC bilayers (PC-Chol bilayers). These studies provided detailed information about the structure and dynamics of the membrane/water interface and the hydrocarbon chain region in bilayers built of various types of PCs and Chol. This enabled studies of membrane systems of higher complexity. They included the investigation of an integral membrane protein in its natural environment of a PC bilayer, and the antibacterial activity of magainin-2. The latter study required the construction of a model bacterial membrane which consisted of two types of phospholipids and counter ions. Whenever published experimental data were available, the results of the simulations were compared with them.     

Cholesterol effects on the phosphatidylcholine bilayer nonpolar region: a molecular simulation study

A 15-ns molecular dynamics (MD) simulation of the fully hydrated dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer in the liquid-crystalline state was carried out to investigate the effect of Chol on the hydrocarbon chain region of the bilayer. The last 8-ns fragment of the generated trajectory was used for analyses. As a reference system, a pure DMPC bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) simulated for 14 ns was used. The study shows that a Chol-induced increase of the bulk molecular order parameter along both beta- and gamma-chain is mainly caused by a decrease of the average tilt of the chains, because the bulk average number of gauche rotamers/myristoyl chain is not significantly changed by Chol. Nevertheless, for DMPCs located near Chol molecules both the number of gauche rotamers/chain and the chain tilt are decreased. The magnitude of the Chol effect on the PC alkyl chains depends, in addition to the PC-Chol distance, on the side of the Chol molecule (alpha- or beta-face) that the chains are in contact with. This study provides some new insight into the properties of the coexistence region of the partial phase diagram for DMPC-Chol bilayers.     

Cholesterol effects on the phospholipid condensation and packing in the bilayer: a molecular simulation study

A 15-ns molecular dynamics simulation of the fully hydrated liquid-crystalline dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. The generated trajectory was analysed to investigate the mechanism of the Chol condensing effect on DMPC hydrocarbon chains and the influence of Chol on the chain packing in the membrane. Chol was found to induce stronger van der Waals interactions among the chains, whereas its interactions with the chains were weak. In the DMPC-Chol bilayer, as in the DMPC bilayer, DMPC chains were regularly packed around a chosen chain but around a Chol molecule they were not. DMPC gamma chains made closer contacts with Chol than the beta chains.     

Nonpolar interactions between trans-membrane helical EGF peptide and phosphatidylcholines, sphingomyelins and cholesterol. Molecular dynamics simulation studies.

A molecular dynamics simulation study of four lipid bilayers with inserted trans-membrane helical fragment of epithelial growth factor (EGF) receptor (EGF peptide) was performed. The lipid bilayers differ in their lipid composition and consist of (i) unsaturated phosphatidylcholine (palmitoyloleoylphosphatidylcholine, POPC), (ii) POPC and 20 mol% of cholesterol (Chol), (iii) sphingomyelin (SM) and 20 mol% of Chol, and (iv) SM and 50 mol% of Chol. Only 1 out of 26 residues in the EGF-peptide sequence is polar (Thr). The hydrophobic thickness of each bilayer is different but shorter than the length of the peptide and so, due to hydrophobic mismatch, the inserted peptide is tilted in each bilayer. Additionally, in the POPC bilayer, which is the thinnest, the peptide loses its helical structure in a short three-amino acid fragment. This facilitates bending of the peptide and burying all hydrophobic amino acids inside the membrane core (Figure 1(b)). Bilayer lipid composition affects interactions between the peptide and lipids in the membrane core. Chol increases packing of atoms relative to the peptide side chains, and thus increases van der Waals interactions. On average, the packing around the peptide is higher in SM-based bilayers than POPC-based bilayers but for certain amino acids, packing depends on their position relative to the bilayer center. In the bilayer center, packing is higher in POPC-based bilayers, while in regions closer to the interface packing is higher in SM-based bilayers. In general, amino acids with larger side chains interact strongly with lipids, and thus the peptide sequence is important for the pattern of interactions at different membrane depths. This pattern closely resembles the shape of recently published lateral pressure profiles [Ollila et alJ. Struct. Biol. DOI:10.1016/j.jsb.2007.01.012].     

Role of phosphatidylglycerols in the stability of bacterial membranes

An extensive 100-ns molecular dynamics simulation of lipid bilayer composed of mixture of phosphatidylethanolamine (PE) and phosphatidylglycerol (PG) was performed to elucidate the role of PGs to the stability of bacterial membranes. In addition, a control simulation of pure PE over 150 ns was performed. We observed that PGs decrease both the PE headgroup protrusions into the water phase, and the PE headgroup motion along bilayer normal. The above effects are caused by stronger inter-lipid interactions in the mixed bilayer: the number of hydrogen bonds created by PEs is 34% higher in the mixed than in the pure bilayer. Another contribution is due to the numerous ion-mediated inter-lipid links, which strongly enhance interface stability. That provides a plausible mechanism for preventing lipid desorption from the membrane, for example, under the influence of an organic solvent. A more compact and less dynamic interface structure also decreases membrane permeability. That provides a possible mechanism for stabilizing, e.g., bacterial membranes.     

Molecular dynamics simulations of pentapeptides at interfaces: salt bridge and cation-pi interactions

Peptide-membrane interactions are important for understanding the binding, partitioning, and folding of membrane proteins; the activity of antimicrobial and fusion peptides; and a number of other processes. We describe molecular dynamics simulations (10-25 ns) of two pentapeptides Ace-WLXLL (with X = Arg or Lys side chain) (White, S. H., and Wimley, W.C. (1996) Nat. Struct. Biol. 3, 842-848) in water and three different membrane mimetic systems: (i) a water/cyclohexane interface, (ii) water-saturated octanol, and (iii) a solvated dioleoylphosphatidylcholine bilayer. A salt bridge is found between the protonated Arg or Lys side chains with the carboxyl terminus at the three interfaces. In water/cyclohexane, the salt bridge is most exposed to the water phase and least stable. In water/octanol and the lipid bilayer systems, the salt bridge once formed persists throughout the simulations. In the lipid bilayer, the salt bridge is more stable when the peptide penetrates deeper into the bilayer. In one of two peptides, a cation-pi interaction between the Arg and the Trp side chains is stable in the lipid bilayer for about 15 ns before breaking. In all cases, the conformations of the peptides are restricted by their presence at the interface and can be assigned to a few major conformational clusters. Side chains facing the water phase are most mobile. In the lipid bilayer, the peptides remain in the interface area, where they overlap with the carbonyl area of the lipid bilayer and perturb the local density profile of the bilayer. The tryptophan side chain remains in the water-lipid interface, where it interacts with the lipid choline group and forms hydrogen bonds with the ester carbonyl of the lipid and with water in the interface.     

Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer

The glycophorin helix dimer is a paradigm for the exploration of helix-helix interactions in integral membrane proteins. Two NMR structures of the dimer are known, one in a detergent micelle and one in a lipid bilayer. Multiple (4 x 50 ns) molecular dynamics simulations starting from each of the two NMR structures, with each structure in either a dodecyl phosphocholine (DPC) micelle or a dimyristoyl phosphatidylcholine (DMPC) bilayer, have been used to explore the conformational dynamics of the helix dimer. Analysis of the helix-helix interaction, mediated by the GxxxG sequence motif, suggests convergence of the simulations to a common model. This is closer to the NMR structure determined in a bilayer than to micelle structure. The stable dimer interface in the final simulation model is characterized by (i) Gly/Gly packing and (ii) Thr/Thr interhelix H-bonds. These results demonstrate the ability of extended molecular dynamics simulations in a lipid bilayer environment to refine membrane protein structures or models derived from experimental data obtained in protein/detergent micelles.     

Non-polar interactions between cholesterol and phospholipids: a molecular dynamics simulation study

A 15-ns molecular dynamics simulation of the fully hydrated dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. An 8-ns trajectory was analysed to investigate the effect of Chol on the chain packing in the bilayer core. While the packing of DMPC chains on the smooth alpha-face side of the Chol ring is similar to that in the pure DMPC bilayer, the packing on the rough beta-face side is less regular and less tight. Two methyl groups located on the Chol beta-face disturb the packing; in effect, van der Waals (vdW) interactions between Chol rings and DMPC chains are weaker than the ones between sole DMPC chains. VdW interactions between an alkyl chain of DMPC and an isooctyl tail of Chol are similarly strong as those between two DMPC chains.     

Vstx1, a modifier of Kv channel gating, localizes to the interfacial region of lipid bilayers

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.     

Cholesterol-sphingomyelin interactions: a molecular dynamics simulation study

Stearoylsphingomyelin (SSM) bilayers containing 0, 22, and 50 mol % cholesterol (Chol) and a pentadecanoyl-stearoylphosphatidylcholine (15SPC) bilayer containing 22 mol % Chol were molecular dynamics simulated at two temperatures (37 degrees C and 60 degrees C). 15SPC is the best PC equivalent of SSM. The Chol effect on the SSM bilayer differs significantly from that on the 15SPC bilayer. At the same temperature and Chol content, H-bonding of Chol with SSM is more extensive than with 15SPC. SSM-Chol H-bonding anchors the OH group of Chol in the lower regions of the SSM-Chol bilayer interface. Such a location strengthens the influence of Chol on the SSM chains. In effect, the phase of the SSM-Chol bilayer containing 22 mol % Chol at 37 degrees C is shifted from the gel to the liquid-ordered phase, and the bilayer displays similar properties below and above the main phase-transition temperature for a pure SSM bilayer of approximately 45 degrees C. In contrast, due to a higher location, Chol is not able to change the phase of the 15SPC-Chol bilayer, which at 37 degrees C remains in the gel phase. Chol affects both the core and interface of the SSM bilayer. With increasing Chol content, the order of SSM chains and hydration of SSM headgroups increase, whereas polar interactions between lipids decrease.     

Differential dynamic and structural behavior of lipid-cholesterol domains in model membranes

Changes in the cholesterol (Chol) content of biological membranes are known to alter the physicochemical properties of the lipid lamella and consequently the function of membrane-associated enzymes. To characterize these changes, we used steady-state and time resolved fluorescence spectroscopy and two photon-excitation microscopy techniques. The membrane systems were chosen according to the techniques that were used: large unilamellar vesicles (LUVs) for cuvette and giant unilamellar vesicles (GUVs) for microscopy measurements; they were prepared from dipalmitoyl phosphatidylcholine (DPPC) and dioctadecyl phosphatidylcholine (DOPC) in mixtures that are well known to form lipid domains. Two fluorescent probes, which insert into different regions of the bilayer, were selected: 1,6-diphenyl-1,3,5-hexatriene (DPH) was located at the deep hydrophobic core of the acyl chain regions and 2-dimethylamino-6-lauroylnaphthalene (Laurdan) at the hydrophilic-hydrophobic membrane interface. Our spectroscopy results show that (i) the changes induced by cholesterol in the deep hydrophobic phospholipid acyl chain domain are different from the ones observed in the superficial region of the hydrophilic-hydrophobic interface, and these changes depend on the state of the lamella and (ii) the incorporation of cholesterol into the lamella induces an increase in the orientation dynamics in the deep region of the phospholipid acyl chains with a corresponding decrease in the orientation at the region close to the polar lipid headgroups. The microscopy data from DOPC/DPPC/Chol GUVs using Laurdan generalized polarization (Laurdan GP) suggest that a high cholesterol content in the bilayer weakens the stability of the water hydrogen bond network and hence the stability of the liquid-ordered phase (Lo).     

Saturation with cholesterol increases vertical order and smoothes the surface of the phosphatidylcholine bilayer: a molecular simulation study

Molecular dynamics (MD) simulations of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer and a POPC bilayer containing 50mol% cholesterol (POPC-Chol50) were carried out for 200ns to compare the spatial organizations of the pure POPC bilayer and the POPC bilayer saturated with Chol. The results presented here indicate that saturation with Chol significantly narrows the distribution of vertical positions of the center-of-mass of POPC molecules and POPC atoms in the bilayer. In the POPC-Chol50 bilayer, the same moieties of the lipid molecules are better aligned at a given bilayer depth, forming the following clearly separated membrane regions: the polar headgroup, the rigid core consisting of steroid rings and upper fragments of the acyl chains, and the fluid hydrocarbon core consisting of Chol chains and the lower fragments of POPC chains. The membrane surface of the POPC-Chol50 bilayer is smooth. The results have biological significance because the POPC-Chol50 bilayer models the bulk phospholipid portion of the fiber-cell membrane in the eye lens. It is hypothesized that in the eye lens cholesterol-induced smoothing of the membrane surface decreases light-scattering and helps to maintain lens transparency.     

Numerical studies of the membrane fluorescent dyes dynamics in ground and excited states

Fluorescence methods are widely used in studies of biological and model membranes. The dynamics of membrane fluorescent markers in their ground and excited electronic states and correlations with their molecular surrounding within the fully hydrated phospholipid bilayer are still not well understood. In the present work, Quantum Mechanical (QM) calculations and Molecular Dynamics (MD) simulations are used to characterize location and interactions of two membrane polarity probes (Prodan; 6-propionyl-2-dimethylaminonaphthalene and its derivative Laurdan; 2-dimethylamino-6-lauroylnaphthalene) with the dioleoylphosphatidylcholine (DOPC) lipid bilayer model. MD simulations with fluorophores in ground and excited states are found to be a useful tool to analyze the fluorescent dye dynamics and their immediate vicinity. The results of QM calculations and MD simulations are in excellent agreement with available experimental data. The calculation shows that the two amphiphilic dyes initially placed in bulk water diffuse within 10 ns towards their final location in the lipid bilayer. Analysis of solvent relaxation process in the aqueous phase occurs on the picoseconds timescale whereas it takes nanoseconds at the lipid/water interface. Four different relaxation time constants, corresponding to different relaxation processes, where observed when the dyes were embedded into the membrane.     

Lipid bilayer deformation and the free energy of interaction of a Kv channel gating-modifier toxin

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein.     

Sterol chemical configuration influences the thermotropic phase behaviour of dipalmitoylphosphatidylcholine bilayers containing 5α-cholestan-3β- and 3α-ol

It is commonly believed that all membrane sterols are rigid all-trans ring systems with a fully extended alkyl side-chain and that they similarly influence phospholipid bilayer physical properties. Here, we report the sterol concentration-dependent, thermotropic phase behaviour of binary dipalmitoylphosphatidylcholine (DPPC)/sterol mixtures containing two similar 5α-H sterols with different functional group orientations (3α-OH or 3β-OH), which adopt an ideal all-trans planar ring conformation but lack the deformed ring B conformation of cholesterol (Chol) and epicholesterol (Echol), using differential scanning calorimetry (DSC). Our deconvolution of the DSC main phase transition endotherms show differences in the proportions of sterol-poor (sharp) and sterol-rich (broad) domains in the DPPC bilayer with increasing sterol concentration, which delineate gel/liquid-crystalline (P_β′/L_α) and disordered gel (L_β)/liquid-ordered (l_o) phase regions. There are similarities in the DPPC main phase transition temperature, cooperativity and enthalpy for each 3β-ol and 3α-ol pair with increasing sterol concentration and differences in the parameters obtained for both the sterol-poor and sterol-rich regions. The sterol-poor domain persists over a greater concentration range in both 3α-ol/DPPC mixtures, suggesting that either those domains are more stable in the 3α-ols or that those sterols are less miscible in the sterol-rich domain. Corresponding parameters for the sterol-rich domain show that at sterol concentrations up to 20 mol%, the 5α-H,3β-ol is more effective at reducing the phase transition enthalpy of the broad component () than Chol, but is less effective at higher concentrations. Although mixtures containing Echol and 5α-cholestan-3α-ol have similar positive slopes below 7 mol% sterol, suggesting that they abolish the L_β/l_o phase transition equally effectively at low concentrations, Echol is more effective than the saturated 3α-ol at higher sterol concentrations. A comparison of obtained for the saturated and unsaturated pairs suggests that the latter sterols stabilize the l_o phase and broaden and abolish the DPPC main phase transition more effectively than the saturated sterols at physiologically relevant concentrations, supporting the idea that the double bond of Chol and Echol promotes greater sterol miscibility and the formation of l_o phase lipid bilayers relative to corresponding saturated sterols in biological membranes.     

Structural Effects and Translocation of Doxorubicin in a DPPC/Chol Bilayer: The Role of Cholesterol

We use molecular dynamics simulations to characterize the influence of cholesterol (Chol) on the interaction between the anticancer drug doxorubicin (DOX) and a dipalmitoyl phosphatidylcholine/Chol lipid bilayer. We calculate the potential of mean force, which gives us an estimate of the free energy barrier for DOX translocation across the membrane. We find free energy barriers of 23.1 ± 3.1 k_BT, 36.8 ± 5.1 k_BT, and 54.5 ± 4.7 k_BT for systems composed of 0%, 15%, and 30% Chol, respectively. Our predictions agree with Arrhenius activation energies from experiments using phospholipid membranes, including 20 k_BT for 0% Chol and 37.2 k_BT for 20% Chol. The location of the free energy barrier for translocation across the bilayer is dependent on composition. As Chol concentration increases, this barrier changes from the release of DOX into the water to flip-flop over the membrane center. The drug greatly affects local membrane structure by attracting dipalmitoyl phosphatidylcholine headgroups, curving the membrane, and allowing water penetration. Despite its hydrophobicity, DOX facilitates water transport via its polar groups.     

Constant-pressure molecular dynamics investigation of cholesterol effects in a dipalmitoylphosphatidylcholine bilayer.

We report a 1.4-ns constant-pressure molecular dynamics simulation of cholesterol at 12.5 mol% in a dipalmitoylphosphatidylcholine (DPPC) bilayer at 50 degrees C and compare the results to our previous simulation of a pure DPPC bilayer. The interlamellar spacing was increased by 2.5 A in the cholesterol-containing bilayer, consistent with x-ray diffraction results, whereas the bilayer thickness was increased by only 1 A. The bilayer/water interface was more abrupt because the lipid headgroups lie flatter to fill spaces left by the cholesterol molecules. This leads to less compensation by the lipid headgroups of the oriented water contribution to the membrane dipole potential and could explain the experimentally observed increase in the magnitude of the dipole potential by cholesterol. Our calculations suggested that 12.5 mol% cholesterol does not significantly affect the conformations and packing of the hydrocarbon chains and produces only a slight reduction in the empty free volume. However, cholesterol has a significant influence on the subnanosecond time scale lipid dynamics: the diffusion constant for the center-of-mass "rattling" motion was reduced by a factor of 3, and the reorientational motion of the methylene groups was slowed along the entire length of the hydrocarbon chains.