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

Investigating the competitive effects of cholesterol and melatonin in model lipid membranes

We have studied the impact of cholesterol and/or melatonin on the static and dynamical properties of bilayers made of DPPC or DOPC utilizing neutron scattering techniques, Raman spectroscopy and molecular dynamics simulations. While differing in the amplitude of the effect due to cholesterol or melatonin when comparing their interactions with the two lipids, their addition ensued recognizable changes to both types of bilayers. As expected, based on the two-component systems of lipid/cholesterol or lipid/melatonin studied previously, we show the impact of cholesterol and melatonin being opposite and competitive in the case of three-component systems of lipid/cholesterol/melatonin. The effect of cholesterol appears to prevail over that of melatonin in the case of structural properties of DPPC-based bilayers, which can be explained by its interactions targeting primarily the saturated lipid chains. The dynamics of hydrocarbon chains represented by the ratio of trans/gauche conformers reveals the competitive effect of cholesterol and melatonin being somewhat more balanced. The additive yet opposing effects of cholesterol and melatonin have been observed also in the case of structural properties of DOPC-based bilayers. We report that cholesterol induced an increase in bilayer thickness, while melatonin induced a decrease in bilayer thickness in the three-component systems of DOPC/cholesterol/melatonin. Commensurately, by evaluating the projected area of DOPC, we demonstrate a lipid area decrease with an increasing concentration of cholesterol, and a lipid area increase with an increasing concentration of melatonin. The demonstrated condensing effect of cholesterol and the fluidizing effect of melatonin appear in an additive manner upon their mutual presence.

A quantitative determination of lipid bilayer deposition efficiency using AFM

The efficacy of a number of different methods for depositing a dimyristoylphosphatidylcholine (DMPC) lipid bilayer or DMPC-cholesterol (3 : 1) mixed bilayer onto a silicon substrate has been investigated in a quantitative manner using atomic force microscopy (AFM) image analysis to extract surface coverage. Complementary AFM-IR measurements were used to confirm the presence of the lipids. For the Langmuir-Blodgett/Schaefer deposition method at temperatures below the chain-melting transition temperature (T m), a large number of bilayer defects resulted when DMPC was deposited from a water subphase. Addition of calcium ions to the trough led to smaller, more frequent defects, whereas addition of cholesterol to the lipid mixture led to a vast improvement in bilayer coverage. Poor coverage was achieved for deposition at temperatures above T m. Formation of the deposited bilayer from vesicle fusion proved a more reliable method for all systems, with formation of near-complete bilayers within 60 seconds at temperatures above T m, although this method led to a higher probability of multilayer formation and rougher bilayer surfaces.

Investigation of the Membrane Fluidity Regulation of Fatty Acid Intracellular Distribution by Fluorescence Lifetime Imaging of Novel Polarity Sensitive Fluorescent Derivatives

Free fatty acids are essential structural components of the cell, and their intracellular distribution and effects on membrane organelles have crucial roles in regulating the metabolism, development, and cell cycle of most cell types. Here we engineered novel fluorescent, polarity-sensitive fatty acid derivatives, with the fatty acid aliphatic chain of increasing length (from 12 to 18 carbons). As in the laurdan probe, the lipophilic acyl tail is connected to the environmentally sensitive dimethylaminonaphthalene moiety. The fluorescence lifetime imaging analysis allowed us to monitor the intracellular distribution of the free fatty acids within the cell, and to simultaneously examine how the fluidity and the microviscosity of the membrane environment influence their localization. Each of these probes can thus be used to investigate the membrane fluidity regulation of the correspondent fatty acid intracellular distribution. We observed that, in PC-12 cells, fluorescent sensitive fatty acid derivatives with increased chain length compartmentalize more preferentially in the fluid regions, characterized by a low microviscosity. Moreover, fatty acid derivatives with the longest chain compartmentalize in lipid droplets and lysosomes with characteristic lifetimes, thus making these probes a promising tool for monitoring lipophagy and related events.

Surface Activities of a Lipid Analogue Room-Temperature Ionic Liquid and Its Effects on Phospholipid Membrane

There are great efforts of synthesizing imidazolium-based ionic liquids (ILs) for developing new antibiotics as these molecules have shown strong antibacterial activities. Compared to a single-hydrocarbon-chained IL, the lipid analogues (LAs) with two chains are more effective. In the present study, the LA molecule ^MeIm(COOH)^Me(Oleylamine)Iodide has been synthesized and its surface activities along with the effectiveness in restructuring of a model cellular membrane have been quantified. The molecule is found to be highly surface active as estimated from the area-pressure isotherm of a monolayer of the molecules formed at the air-water interface. The X-ray reflectivity (XRR) studies of a monolayer dip-coated on a hydrophilic substrate have shown the structural properties of the layer which resembles to those of unsaturated phospholipids. The LA molecules are observed to fluidize a phospholipid bilayer formed by the saturated lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). At a lower surface pressure, the lipid monolayer of DPPC has exhibited a thickening effect at a low concentration of added LA and a thinning effect at higher concentration. However, at a high surface pressure of the monolayer, the thickness is found to decrease monotonically. The in-plane pressure-dependent interaction of LA molecules with model cellular membrane and the corresponding perturbation in the structure and physical properties of the membrane may be linked to the strong lysing effect of these types of molecules.

Influence of Cholesterol and Bilayer Curvature on the Interaction of PPO-PEO Block Copolymers with Liposomes

Interactions of nonionic poly(ethylene oxide)- b-poly(propylene oxide) (PEO-PPO) block copolymers, known as Pluronics or poloxamers, with cell membranes have been widely studied for a host of biomedical applications. Herein, we report how cholesterol within phosphatidylcholine (POPC) lipid bilayer liposomes and bilayer curvature affects the binding of several PPO-PEO-PPO triblocks with varying PPO content and a tPPO-PEO diblock, where t refers to a tert-butyl end group. Pulsed-field-gradient NMR was employed to quantify the extent of copolymer associated with liposomes prepared with cholesterol concentrations ranging from 0 to 30 mol % relative to the total content of POPC and cholesterol and vesicle extrusion radii of 25, 50, or 100 nm. The fraction of polymer bound to the liposomes was extracted from NMR data on the basis of the very different mobilities of the bound and free polymers in aqueous solution. Cholesterol concentration was manipulated by varying the molar percentage of this sterol in the POPC bilayer preparation. The membrane curvature was varied by adjusting the liposome size through a conventional pore extrusion technique. Although the PPO content significantly influences the overall amount of block copolymer adsorbed to the liposome, we found that polymer binding decreases with increasing cholesterol concentration in a universal fashion, with the fraction of bound polymer dropping 10-fold between 0 and 30 mol % cholesterol relative to the total content of POPC and cholesterol. Increasing the bilayer curvature (decreasing the radius of the liposome) in the absence of cholesterol increases polymer binding between 2- and 4-fold over the range of liposome sizes studied. These results demonstrate that cholesterol plays a dominant role, and bilayer curvature has a less significant impact as the curvature decreases, on polymer-membrane association.

Spatially Resolving the Condensing Effect of Cholesterol in Lipid Bilayers

We study the effect of cholesterol on the structure of dipalmitoylphosphatidylcholine phospholipid bilayers. Using extensive molecular dynamics computer simulations at atomistic resolution, we observe and quantify several structural changes upon increasing cholesterol content that are collectively known as the condensing effect: a thickening of the bilayer, an increase in lipid tail order, and a decrease in lateral area. We also observe a change in leaflet interdigitation and a lack thereof in the distributions of dipalmitoylphosphatidylcholine headgroup orientations. These results, obtained over a wide range of cholesterol mole fractions, are then used to calibrate the analysis of phospholipid properties in bilayers containing a single cholesterol molecule per leaflet, which we perform in a spatially resolved way. We find that a single cholesterol molecule affects phospholipids in its first and second solvation shells, which puts the range of this effective interaction to be on the order of 1-2 nm. We also observe a tendency of phospholipids to orient their polar headgroups toward the cholesterol, which provides additional support for the umbrella model of bilayer organization.

Direct Measurement of the Effect of Cholesterol and 6-Ketocholestanol on the Membrane Dipole Electric Field Using Vibrational Stark Effect Spectroscopy Coupled with Molecular Dynamics Simulations

Biological membranes are heterogeneous structures with complex electrostatic profiles arising from lipids, sterols, membrane proteins, and water molecules. We investigated the effect of cholesterol and its derivative 6-ketocholestanol (6-kc) on membrane electrostatics by directly measuring the dipole electric field (F⃗d) within lipid bilayers containing cholesterol or 6-kc at concentrations of 0-40 mol% through the vibrational Stark effect (VSE). We found that adding low concentrations of cholesterol, up to ∼10 mol %, increases F⃗d, while adding more cholesterol up to 40 mol% lowers F⃗d. In contrast, we measured a monotonic increase in F⃗d as 6-kc concentration increased. We propose that this membrane electric field is affected by multiple factors: the polarity of the sterol molecules, the reorientation of the phospholipid dipole due to sterol, and the impact of the sterol on hydrogen bonding with surface water. We used molecular dynamics simulations to examine the distribution of phospholipids, sterol, and helix in bilayers containing these sterols. At low concentrations, we observed clustering of sterols near the vibrational probe whereas at high concentrations, we observed spatial correlation between the positions of the sterol molecules. This work demonstrates how a one-atom difference in a sterol changes the physicochemical and electric field properties of the bilayer.

High Cholesterol/Low Cholesterol: Effects in Biological Membranes: A Review

Lipid composition determines membrane properties, and cholesterol plays a major role in this determination as it regulates membrane fluidity and permeability, as well as induces the formation of coexisting phases and domains in the membrane. Biological membranes display a very diverse lipid composition, the lateral organization of which plays a crucial role in regulating a variety of membrane functions. We hypothesize that, during biological evolution, membranes with a particular cholesterol content were selected to perform certain functions in the cells of eukaryotic organisms. In this review, we discuss the major membrane properties induced by cholesterol, and their relationship to certain membrane functions.

A computer model of a polyunsaturated monogalactolipid bilayer

1,2-di-O-acyl-3-O-β-D-galactopyranosyl-sn-glycerol (MGDG) is the main lipid component of thylakoid membranes of higher plants and algae. This monogalactolipid is thought of as a non-bilayer lipid but actually it can form both lamellar and nonlamellar phases. In this study, molecular dynamics (MD) simulations of the fully hydrated di-18:3 MGDG bilayer in the lamellar phase were carried out at 310 and 295 K for 200 and 450 ns, respectively, using the GROMACS 4 software package and OPLS-AA force field. At both temperatures, the lamellar phase of the systems was stable. The pure di-18:3 MGDG bilayer is the first step towards creating a computer model of the lipid matrix of the thylakoid membrane and the main aim of this study was to validate the computer model of di-18:3 MGDG in the bilayer and also to assess the properties of the bilayer. However, only a few of the predicted properties could be compared with those derived experimentally and in other MD simulations because of insufficient amount of such data. Thus, direct validation of the MGDG bilayer proved difficult. Therefore, in the validation process also an indirect approach was used, in which a computer model of the 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayer simulated at the same temperatures using the same force field as the MGDG bilayer was assessed. Successful validation of the DOPC bilayer parameterized in the OPLS-AA force field and similar properties of the MGDG molecules in the pure 18:3 MGDG and in binary 18:3 MGDG-PC bilayers indicate that the computer model of the MDGD molecule is faithful and the MGDG bilayer is representative on the time scales covered in these MD simulations.

Building Synthetic Sterols Computationally - Unlocking the Secrets of Evolution?

Cholesterol is vital in regulating the physical properties of animal cell membranes. While it remains unclear what renders cholesterol so unique, it is known that other sterols are less capable in modulating membrane properties, and there are membrane proteins whose function is dependent on cholesterol. Practical applications of cholesterol include its use in liposomes in drug delivery and cosmetics, cholesterol-based detergents in membrane protein crystallography, its fluorescent analogs in studies of cholesterol transport in cells and tissues, etc. Clearly, in spite of their difficult synthesis, producing the synthetic analogs of cholesterol is of great commercial and scientific interest. In this article, we discuss how synthetic sterols non-existent in nature can be used to elucidate the roles of cholesterol’s structural elements. To this end, we discuss recent atomistic molecular dynamics simulation studies that have predicted new synthetic sterols with properties comparable to those of cholesterol. We also discuss more recent experimental studies that have vindicated these predictions. The paper highlights the strength of computational simulations in making predictions for synthetic biology, thereby guiding experiments.

Evaluation of a possible role of Stigmatella aurantiaca ACE in Aβ peptide degradation: a molecular modeling approach

Amyloid-β (Aβ)-degrading enzymes are known to degrade Aβ peptides, a causative agent of Alzheimer's disease. These enzymes are responsible for maintaining Aβ concentration. However, loss of such enzymes or their Aβ-degrading activity because of certain genetic as well as nongenetic reasons initiates the accumulation of Aβ peptides in the human brain. Considering the limitations of the human enzymes in clearing Aβ peptide, the search for microbial enzymes that could cleave Aβ is necessary. Hence, we built a three-dimensional model of angiotensin-converting enzyme (ACE) from Stigmatella aurantiaca using homology modeling technique. Molecular docking and molecular dynamics simulation techniques were used to outline the possible cleavage mechanism of Aβ peptide. These findings suggest that catalytic residue Glu 434 of the model could play a crucial role to degrade Aβ peptide between Asp 7 and Ser 8. Thus, ACE from S. aurantiaca might cleave Aβ peptides similar to human ACE and could be used to design new therapeutic strategies against Alzheimer's disease.

Insight into molecular interactions of Aβ peptide and gelatinase from Enterococcus faecalis: a molecular modeling approach

Alzheimer's disease is characterized by the presence of extracellular deposition of amyloid beta (Aβ) peptides.

Sterol biosynthesis is required for heat resistance but not extracellular survival in leishmania

Sterol biosynthesis is a crucial pathway in eukaryotes leading to the production of cholesterol in animals and various C24-alkyl sterols (ergostane-based sterols) in fungi, plants, and trypanosomatid protozoa. Sterols are important membrane components and precursors for the synthesis of powerful bioactive molecules, including steroid hormones in mammals. Their functions in pathogenic protozoa are not well characterized, which limits the development of sterol synthesis inhibitors as drugs. Here we investigated the role of sterol C14α-demethylase (C14DM) in Leishmania parasites. C14DM is a cytochrome P450 enzyme and the primary target of azole drugs. In Leishmania, genetic or chemical inactivation of C14DM led to a complete loss of ergostane-based sterols and accumulation of 14-methylated sterols. Despite the drastic change in lipid composition, C14DM-null mutants (c14dm(-)) were surprisingly viable and replicative in culture. They did exhibit remarkable defects including increased membrane fluidity, failure to maintain detergent resistant membrane fraction, and hypersensitivity to heat stress. These c14dm(-) mutants showed severely reduced virulence in mice but were highly resistant to itraconazole and amphotericin B, two drugs targeting sterol synthesis. Our findings suggest that the accumulation of toxic sterol intermediates in c14dm(-) causes strong membrane perturbation and significant vulnerability to stress. The new knowledge may help improve the efficacy of current drugs against pathogenic protozoa by exploiting the fitness loss associated with drug resistance.

Refined OPLS all-atom force field for saturated phosphatidylcholine bilayers at full hydration

We report parametrization of dipalmitoyl-phosphatidylcholine (DPPC) in the framework of the Optimized Parameters for Liquid Simulations all-atom (OPLS-AA) force field. We chose DPPC as it is one of the most studied phospholipid species and thus has plenty of experimental data necessary for model validation, and it is also one of the highly important and abundant lipid types, e.g., in lung surfactant. Overall, PCs have not been previously parametrized in the OPLS-AA force field; thus, there is a need to derive its bonding and nonbonding parameters for both the polar and nonpolar parts of the molecule. In the present study, we determined the parameters for torsion angles in the phosphatidylcholine and glycerol moieties and in the acyl chains, as well the partial atomic charges. In these calculations, we used three methods: (1) Hartree-Fock (HF), (2) second order Møller-Plesset perturbation theory (MP2), and (3) density functional theory (DFT). We also tested the effect of the polar environment by using the polarizable continuum model (PCM), and for acyl chains the van der Waals parameters were also adjusted. In effect, six parameter sets were generated and tested on a DPPC bilayer. Out of these six sets, only one was found to be able to satisfactorily reproduce experimental data for the lipid bilayer. The successful DPPC model was obtained from MP2 calculations in an implicit polar environment (PCM).

Cholesterol-induced condensing and disordering effects on a rigid catanionic bilayer: a molecular dynamics study

Molecular dynamics simulation is applied to investigate the bilayer properties of a novel catanionic vesicle composed of an ion pair amphiphile, hexadecyltrimethylammonium-dodecylsulfate (HTMA-DS), with cholesterol. Structural properties, such as molecular organization, orientation, and conformation, were analyzed from the resulting trajectory. Simulation results showed that cholesterol could induce both condensing and disordering effects on the rigid HTMA-DS bilayer. The condensing effect of cholesterol was ascribed to the maximizing contact between cholesterol ring and the neighboring hydrocarbon chains. Thus, the inserted cholesterol ring restrained the neighboring hydrocarbon chain segments from motion and increased the order of the neighboring hydrocarbon chains. However, the presence of cholesterol would increase the distance between head groups of HTMA-DS and induce a shift of DS(-) head groups toward the inside of the bilayer. This led to the protrusion of the HTMA(+) head groups and conformational disorder in the front segments of HTMA(+) hydrocarbon chains. In addition, the cholesterol-induced void in the hydrophobic core of the HTMA-DS bilayer increased the motion freedom of the terminal segments of the hydrocarbon chains. The cholesterol-induced space in the polar region and void in the nonpolar region of the bilayer led to a conformational disorder. With high cholesterol contents, the conformational disorder effect would overwhelm the condensing effect, resulting in the apparent disordering effect on the rigid HTMA-DS bilayer.

The molecular structure of the liquid-ordered phase of lipid bilayers

Molecular dynamics simulations reveal substructures within the liquid-ordered phase of lipid bilayers. These substructures, identified in a 10 μs all-atom trajectory of liquid-ordered/liquid-disordered coexistence (L(o)/L(d)) are composed of saturated hydrocarbon chains packed with local hexagonal order and separated by interstitial regions enriched in cholesterol and unsaturated chains. Lipid hydrocarbon chain order parameters calculated from the L(o) phase are in excellent agreement with (2)H NMR measurements; the local hexagonal packing is also consistent with (1)H-MAS NMR spectra of the L(o) phase, NMR diffusion experiments, and small-angle X-ray and neutron scattering. The balance of cholesterol-rich to local hexagonal order is proposed to control the partitioning of membrane components into the L(o) regions. The latter have been frequently associated with formation of so-called rafts, platforms in the plasma membranes of cells that facilitate interaction between components of signaling pathways.

Carotenoids, versatile components of oxygenic photosynthesis

Carotenoids (CARs) are a group of pigments that perform several important physiological functions in all kingdoms of living organisms. CARs serve as protective agents, which are essential structural components of photosynthetic complexes and membranes, and they play an important role in the light harvesting mechanism of photosynthesizing plants and cyanobacteria. The protection against reactive oxygen species, realized by quenching of singlet oxygen and the excited states of photosensitizing molecules, as well as by the scavenging of free radicals, is one of the main biological functions of CARs. X-ray crystallographic localization of CARs revealed that they are present at functionally and structurally important sites of both the PSI and PSII reaction centers. Characterization of a CAR-less cyanobacterial mutant revealed that while the absence of CARs prevents the formation of PSII complexes, it does not abolish the assembly and function of PSI. CAR molecules assist in the formation of protein subunits of the photosynthetic complexes by gluing together their protein components. In addition to their aforementioned indispensable functions, CARs have a substantial role in the formation and maintenance of proper cellular architecture, and potentially also in the protection of the translational machinery under stress conditions.

Comparative computer simulation study of cholesterol in hydrated unary and binary lipid bilayers and in an anhydrous crystal

Models created with molecular dynamics simulations are used to compare the organization and dynamics of cholesterol (Chol) molecules in three different environments: (1) a hydrated pure Chol bilayer that models the Chol bilayer domain, which is a pure Chol domain embedded in the bulk membrane; (2) a 2-palmitoyl-3-oleoyl-d-glycerol-1-phosphorylcholine bilayer saturated with cholesterol (POPC-Chol50) that models the bulk membrane; (3) a Chol crystal. The computer model of the hydrated pure Chol bilayer is stable on the microsecond time scale. Some structural characteristics of Chol molecules in the Chol bilayer are similar to those in the POPC-Chol50 bilayer (e.g., tilt of Chol rings and chains), while others are similar to those in Chol crystals (e.g., surface area per Chol, bilayer thickness). The key result of this study is that the Chol bilayer has, unexpectedly, a dynamic structure, with Chol mobility similar to that in the POPC-Chol50 bilayer though slower. This is the major difference compared to Chol crystals, where Chol molecules are immobile. Also, water accessibility to Chol-OH groups in the Chol bilayer is not limited. On average, each Chol molecule makes 2.3 hydrogen bonds with water in the Chol bilayer, compared with 1.7 hydrogen bonds in the POPC-Col50 bilayer.

Using spin-label W-band EPR to study membrane fluidity profiles in samples of small volume

Conventional and saturation-recovery (SR) EPR at W-band (94GHz) using phosphatidylcholine spin labels (labeled at the alkyl chain [n-PC] and headgroup [T-PC]) to obtain profiles of membrane fluidity has been demonstrated. Dimyristoylphosphatidylcholine (DMPC) membranes with and without 50 mol% cholesterol have been studied, and the results have been compared with similar studies at X-band (9.4 GHz) (L. Mainali, J.B. Feix, J.S. Hyde, W.K. Subczynski, J. Magn. Reson. 212 (2011) 418-425). Profiles of the spin-lattice relaxation rate (T(1)(-1)) obtained from SR EPR measurements for n-PCs and T-PC were used as a convenient quantitative measure of membrane fluidity. Additionally, spectral analysis using Freed's MOMD (microscopic-order macroscopic-disorder) model (E. Meirovitch, J.H. Freed J. Phys. Chem. 88 (1984) 4995-5004) provided rotational diffusion coefficients (R(perpendicular) and R(||)) and order parameters (S(0)). Spectral analysis at X-band provided one rotational diffusion coefficient, R(perpendicular). T(1)(-1), R(perpendicular), and R(||) profiles reflect local membrane dynamics of the lipid alkyl chain, while the order parameter shows only the amplitude of the wobbling motion of the lipid alkyl chain. Using these dynamic parameters, namely T(1)(-1), R(perpendicular), and R(||), one can discriminate the different effects of cholesterol at different depths, showing that cholesterol has a rigidifying effect on alkyl chains to the depth occupied by the rigid steroid ring structure and a fluidizing effect at deeper locations. The nondynamic parameter, S(0), shows that cholesterol has an ordering effect on alkyl chains at all depths. Conventional and SR EPR measurements with T-PC indicate that cholesterol has a fluidizing effect on phospholipid headgroups. EPR at W-band provides more detailed information about the depth-dependent dynamic organization of the membrane compared with information obtained at X-band. EPR at W-band has the potential to be a powerful tool for studying membrane fluidity in samples of small volume, ~30 nL, compared with a representative sample volume of ~3 μL at X-band.

Impact of sterol tilt on membrane bending rigidity in cholesterol and 7DHC-containing DMPC membranes

Cholesterol is so essential to the proper function of mammalian cell membranes that even strikingly small inborn errors in cholesterol synthesis can be devastating. Here we combine molecular dynamics simulations with small angle x-ray diffraction experiments to compare mixed sterol/DMPC membranes over a wide range of sterol compositions for two types of sterols: cholesterol and its immediate metabolic precursor 7DHC, that differs from cholesterol by one double bond. We find that while most membrane properties are only slightly affected by the replacement of one sterol by the other, the tilt degree of freedom, as gauged by the tilt modulus, is significantly larger for cholesterol than for 7DHC over a large range of concentrations. In silico mutations of one sterol into the other further support these findings. Moreover, bending rigidities calculated from simulations and estimated in experiments show that cholesterol stiffens membranes to a larger extent than 7DHC. We discuss the possible mechanistic link between sterol tilt and the way it impacts the membrane mechanical properties, and comment on how this link may shed light on the way replacement of cholesterol by 7DHC leads to disease.

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.

Membrane fluidity profiles as deduced by saturation-recovery EPR measurements of spin-lattice relaxation times of spin labels

There are no easily obtainable EPR spectral parameters for lipid spin labels that describe profiles of membrane fluidity. The order parameter, which is most often used as a measure of membrane fluidity, describes the amplitude of wobbling motion of alkyl chains relative to the membrane normal and does not contain explicitly time or velocity. Thus, this parameter can be considered as nondynamic. The spin-lattice relaxation rate (T(1)(-1)) obtained from saturation-recovery EPR measurements of lipid spin labels in deoxygenated samples depends primarily on the rotational correlation time of the nitroxide moiety within the lipid bilayer. Thus, T(1)(-1) can be used as a convenient quantitative measure of membrane fluidity that reflects local membrane dynamics. T(1)(-1) profiles obtained for 1-palmitoyl-2-(n-doxylstearoyl)phosphatidylcholine (n-PC) spin labels in dimyristoylphosphatidylcholine (DMPC) membranes with and without 50 mol% cholesterol are presented in parallel with profiles of the rotational diffusion coefficient, R(⊥), obtained from simulation of EPR spectra using Freed's model. These profiles are compared with profiles of the order parameter obtained directly from EPR spectra and with profiles of the order parameter obtained from simulation of EPR spectra. It is shown that T(1)(-1) and R(⊥) profiles reveal changes in membrane fluidity that depend on the motional properties of the lipid alkyl chain. We find that cholesterol has a rigidifying effect only to the depth occupied by the rigid steroid ring structure and a fluidizing effect at deeper locations. These effects cannot be differentiated by profiles of the order parameter. All profiles in this study were obtained at X-band (9.5 GHz).

Structure of mixed phosphatidylethanolamine and cholesterol monolayers in a supported hybrid bilayer membrane studied by sum frequency generation vibrational spectroscopy

The structure of hybrid bilayer membranes (HBMs) containing either a pure cholesterol or mixed cholesterol/dipalmitoylphosphatidylethanolamine (DPPE) proximal layer adsorbed onto an octadecanethiol (ODT) self-assembled monolayer (SAM) on a gold substrate have been investigated by sum frequency generation (SFG) spectroscopy. The HBMs were formed by the adsorption of either a pure cholesterol or mixed DPPE/cholesterol monolayer from the air/water interface of a Langmuir-Blodgett trough at surface pressures of 1, 20, or 40 mN·m(-1). SFG spectra were also recorded of HBMs where cholesterol was replaced by cholesterol-d(7), in which the terminal isopropyl group of the alkyl chain of cholesterol was isotopically labeled. In order to isolate the contribution to the SFG spectra from the cholesterol in the mixed cholesterol/phospholipid films, DPPE-d was used, in which the alkyl chains of the phospholipid were deuterated. The infrared spectra of solvent-cast cholesterol and cholesterol-d(7) films were recorded to aid with assignment of the SFG spectra of the HBMs. Features corresponding to methyl, methylene, and methine stretches of cholesterol were identified in the SFG spectra. Information on the polar orientation of SFG-active groups was obtained from the phases of the spectral features. The structure of the HBMs showed little dependence on the surface pressure at which they were formed. SFG spectra of HBMs with a mixed cholesterol/DPPE proximal layer were very similar to the spectra of HBMs with a pure cholesterol proximal layer, although the features in the spectra were more intense than anticipated for a film with half the number of cholesterol molecules, indicating that DPPE did have some effect on the orientation of cholesterol molecules in the film.

Transbilayer organization of membrane cholesterol at low concentrations: Implications in health and disease

Cholesterol is an essential and representative lipid in higher eukaryotic cellular membranes and is often found distributed nonrandomly in domains in biological membranes. A large body of literature exists on the organization of cholesterol in plasma membranes or membranes with high cholesterol content. However, very little is known about organization of cholesterol in membranes containing low amounts of cholesterol such as the endoplasmic reticulum or inner mitochondrial membranes. In this review, we have traced the discovery and subsequent development of the concept of transbilayer cholesterol dimers (domains) in membranes at low concentrations. We have further discussed the role of membrane curvature and thickness on the transbilayer organization of cholesterol. Interestingly, this type of cholesterol organization could be relevant in cellular sorting and trafficking, and in pathological conditions.

Cholesterol orientation and tilt modulus in DMPC bilayers

We performed molecular dynamics (MD) simulations of hydrated bilayers containing mixtures of dimyristoylphosphatidylcholine (DMPC) and cholesterol at various ratios, to study the effect of cholesterol concentration on its orientation, and to characterize the link between cholesterol tilt and overall phospholipid membrane organization. The simulations show a substantial probability for cholesterol molecules to transiently orient perpendicular to the bilayer normal, and suggest that cholesterol tilt may be an important factor for inducing membrane ordering. In particular, we find that as cholesterol concentration increases (1-40% cholesterol) the average cholesterol orientation changes in a manner strongly (anti)correlated with the variation in membrane thickness. Furthermore, cholesterol orientation is found to be determined by the aligning force exerted by other cholesterol molecules. To quantify this aligning field, we analyzed cholesterol orientation using, to our knowledge, the first estimates of the cholesterol tilt modulus chi from MD simulations. Our calculations suggest that the aligning field that determines chi is indeed strongly linked to sterol composition. This empirical parameter (chi) should therefore become a useful quantitative measure to describe cholesterol interaction with other lipids in bilayers, particularly in various coarse-grained force fields.

Cholesterol induces specific spatial and orientational order in cholesterol/phospholipid membranes

BACKGROUND

In lipid bilayers, cholesterol facilitates the formation of the liquid-ordered phase and enables the formation of laterally ordered structures such as lipid rafts. While these domains have an important role in a variety of cellular processes, the precise atomic-level mechanisms responsible for cholesterol's specific ordering and packing capability have remained unresolved.

METHODOLOGY/PRINCIPAL FINDINGS

Our atomic-scale molecular dynamics simulations reveal that this ordering and the associated packing effects in membranes largely result from cholesterol's molecular structure, which differentiates cholesterol from other sterols. We find that cholesterol molecules prefer to be located in the second coordination shell, avoiding direct cholesterol-cholesterol contacts, and form a three-fold symmetric arrangement with proximal cholesterol molecules. At larger distances, the lateral three-fold organization is broken by thermal fluctuations. For other sterols having less structural asymmetry, the three-fold arrangement is considerably lost.

CONCLUSIONS/SIGNIFICANCE

We conclude that cholesterol molecules act collectively in lipid membranes. This is the main reason why the liquid-ordered phase only emerges for Chol concentrations well above 10 mol% where the collective self-organization of Chol molecules emerges spontaneously. The collective ordering process requires specific molecular-scale features that explain why different sterols have very different membrane ordering properties: the three-fold symmetry in the Chol-Chol organization arises from the cholesterol off-plane methyl groups allowing the identification of raft-promoting sterols from those that do not promote rafts.

Influence of cholesterol on the collective dynamics of the phospholipid acyl chains in model membranes

We have studied the packing and collective dynamics of the phospholipid acyl chains in a model membrane composed of 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) and cholesterol in varied phase state. After a structural characterization of this two-component model bilayer using X-ray reflectivity, we have carried out coherent inelastic neutron scattering to investigate the chain dynamics. Both DMPC/cholesterol membranes exhibited much sharper and more pronounced low-energy inelastic excitations than a pure DMPC membrane. In the high-energy regime above 10 meV, the insertion of cholesterol into the membrane was found to shift the position of the inelastic excitation towards values otherwise found in the pure lipids gel phase. Thus, the dissipative collective short-range dynamics of the acyl chains is strongly influenced by the presence of cholesterol.

Membrane cholesterol in the function and organization of G-protein coupled receptors

Cholesterol is an essential component of higher eukaryotic membranes and plays a crucial role in membrane organization, dynamics and function. The G-protein coupled receptors (GPCRs) are the largest class of molecules involved in signal transduction across membranes, and represent major targets in the development of novel drug candidates in all clinical areas. Membrane cholesterol has been reported to have a modulatory role in the function of a number of GPCRs. Two possible mechanisms have been previously suggested by which membrane cholesterol could influence the structure and function of GPCRs (i) through a direct/specific interaction with GPCRs, or (ii) through an indirect way by altering membrane physical properties in which the receptor is embedded, or due to a combination of both. Recently reported crystal structures of GPCRs have shown structural evidence of cholesterol binding sites. Against this backdrop, we recently proposed a novel mechanism by which membrane cholesterol could affect structure and function of GPCRs. According to our hypothesis, cholesterol binding sites in GPCRs could represent 'nonannular' binding sites. Interestingly, previous work from our laboratory has demonstrated that membrane cholesterol is required for the function of the serotonin(1A) receptor (a representative GPCR), which could be due to specific interaction of the receptor with cholesterol. Based on these results, we envisage that there could be specific/nonannular cholesterol binding site(s) in the serotonin(1A) receptor. We have analyzed putative cholesterol binding sites from protein databases in the serotonin(1A) receptor. Our analysis shows that cholesterol binding sites are inherent characteristic features of serotonin(1A) receptors and are conserved through natural evolution. Progress in deciphering molecular details of the GPCR-cholesterol interaction in the membrane would lead to better insight into our overall understanding of GPCR function in health and disease, thereby enhancing our ability to design better therapeutic strategies to combat diseases related to malfunctioning of GPCRs.

Molecular dynamics study of prolyl oligopeptidase with inhibitor in binding cavity

We used the crystal structure of prolyl oligopeptidase (POP) with bound Z-pro-prolinal (ZPP) inhibitor (Protein Data Bank (PDB) structure 1QFS) to perform an intensive molecular dynamics study of the POP-ZPP complex. We performed 100 ns of simulation with the hemiacetal bond, through which the ZPP is bound to the POP, removed in order to better investigate the binding cavity environment. From basic analysis, measuring the radius of gyration, root mean square deviation, solvent accessible surface area and definition of the secondary structure of protein, we determined that the protein structure is highly stable and maintains its structure over the entire simulation time. This demonstrates that such long time simulations can be performed without the protein structure losing stability. We found that water bridges and hydrogen bonds play a negligible role in binding the ZPP thus indicating the importance of the hemiacetal bond. The two domains of the protein are bound by a set of approximately 12 hydrogen bonds, specific to the particular POP protein.

Volumetric Properties, Compressibilities and Volume Fluctuations in Phospholipid-Cholesterol Bilayers

We conducted detailed measurements of the apparent specific volume of dipalmitoylphosphatidylcholine (DPPC)–cholesterol mixtures in excess water as a function of pressure up to 70MPa (700bar) at 20, 38 and 50°C. The volumetric properties and the isothermal compressibility κ T of the lipid vesicles were determined at cholesterol concentrations, χchol, ranging up to 50 mol. The thermodynamic data are compared with other physico-chemical properties of phospholipid–cholesterol mixtures. Furthermore, the thermodynamic properties of the system are discussed in the light of the various T, χchol–phase diagrams and computer simulation studies published in the literature.

A comparative study of the effect of cholesterol on bicelle model membranes using X-band and Q-band EPR spectroscopy

X-band and Q-band electron paramagnetic resonance (EPR) spectroscopic techniques were used to investigate the structure and dynamics of cholesterol containing phospholipid bicelles based upon molecular order parameters (S(mol)), orientational dependent hyperfine splittings and line shape analysis of the corresponding EPR spectra. The nitroxide spin-label 3-beta-doxyl-5-alpha-cholestane (cholestane) was incorporated into DMPC/DHPC bicelles to report the alignment of bicelles in the static magnetic field. The influence of cholesterol on aligned phospholipid bicelles in terms of ordering, the ease of alignment, phase transition temperature have been studied comparatively at X-band and Q-band. At a magnetic field of 1.25 T (Q-band), bicelles with 20 mol% cholesterol aligned at a much lower temperature (313 K), when compared to 318 K at a 0.35 T field strength for X-band, showed better hyperfine splitting values (18.29 G at X-band vs. 18.55 G at Q-band for perpendicular alignment and 8.25 G at X-band vs. 7.83 G at Q-band for the parallel alignment at 318 K) and have greater molecular order parameters (0.76 at X-band vs. 0.86 at Q-band at 318 K). Increasing cholesterol content increased the bicelle ordering, the bicelle-alignment temperature and the gel to liquid crystalline phase transition temperature. We observed that Q-band is more effective than X-band for studying aligned bicelles, because it yielded a higher ordered bicelle system for EPR spectroscopic studies.

Statistical thermodynamics of biomembranes

An overview of the major issues involved in the statistical thermodynamic treatment of phospholipid membranes at the atomistic level is summarized: thermodynamic ensembles, initial configuration (or the physical system being modeled), force field representation as well as the representation of long-range interactions. This is followed by a description of the various ways that the simulated ensembles can be analyzed: area of the lipid, mass density profiles, radial distribution functions (RDFs), water orientation profile, deuterium order parameter, free energy profiles and void (pore) formation; with particular focus on the results obtained from our recent molecular dynamic (MD) simulations of phospholipids interacting with dimethylsulfoxide (Me(2)SO), a commonly used cryoprotective agent (CPA).

Comparison of cholesterol and its direct precursors along the biosynthetic pathway: effects of cholesterol, desmosterol and 7-dehydrocholesterol on saturated and unsaturated lipid bilayers

Despite extensive studies, the remarkable structure-function relationship of cholesterol in cellular membranes has remained rather elusive. This is exemplified by the fact that the membrane properties of cholesterol are distinctly different from those of many other sterols. Here we elucidate this issue through atomic-scale simulations of desmosterol and 7-dehydrocholesterol (7DHC), which are immediate precursors of cholesterol in its two distinct biosynthetic pathways. While desmosterol and 7DHC differ from cholesterol only by one additional double bond, we find that their influence on saturated lipid bilayers is substantially different from cholesterol. The capability to form ordered regions in a saturated (dipalmitoyl-phosphatidylcholine) membrane is given by cholesterol > 7DHC > desmosterol, indicating the important role of cholesterol in saturated lipid environments. For comparison, in an unsaturated (dioleoyl-phosphatidylcholine) bilayer, the membrane properties of all sterols were found to be essentially identical. Our studies indicate that the different membrane ordering properties of sterols can be characterized by a single experimentally accessible parameter, the sterol tilt. The smaller the tilt, the more ordered are the lipids around a given sterol. The molecular level mechanisms responsible for tilt modulation are found to be related to changes in local packing around the additional double bonds.