Lipid exchange mechanism of the cholesteryl ester transfer protein clarified by atomistic and coarse-grained simulations.

Cholesteryl ester transfer protein (CETP) transports cholesteryl esters, triglycerides, and phospholipids between different lipoprotein fractions in blood plasma. The inhibition of CETP has been shown to be a sound strategy to prevent and treat the development of coronary heart disease. We employed molecular dynamics simulations to unravel the mechanisms associated with the CETP-mediated lipid exchange. To this end we used both atomistic and coarse-grained models whose results were consistent with each other. We found CETP to bind to the surface of high density lipoprotein (HDL) -like lipid droplets through its charged and tryptophan residues. Upon binding, CETP rapidly (in about 10 ns) induced the formation of a small hydrophobic patch to the phospholipid surface of the droplet, opening a route from the core of the lipid droplet to the binding pocket of CETP. This was followed by a conformational change of helix X of CETP to an open state, in which we found the accessibility of cholesteryl esters to the C-terminal tunnel opening of CETP to increase. Furthermore, in the absence of helix X, cholesteryl esters rapidly diffused into CETP through the C-terminal opening. The results provide compelling evidence that helix X acts as a lid which conducts lipid exchange by alternating the open and closed states. The findings have potential for the design of novel molecular agents to inhibit the activity of CETP.

Role of lipids in spheroidal high density lipoproteins.

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.

Membrane simulations mimicking acidic pH reveal increased thickness and negative curvature in a bilayer consisting of lysophosphatidylcholines and free fatty acids.

Phospholipids are key components of biological membranes and their lipolysis with phospholipase A(2) (PLA(2)) enzymes occurs in different cellular pH environments. Since no studies are available on the effect of pH on PLA(2)-modified phospholipid membranes, we performed 50-ns atomistic molecular dynamics simulations at three different pH conditions (pH 9.0, 7.5, and 5.5) using a fully PLA(2)-hydrolyzed phosphatidylcholine (PC) bilayer which consists solely of lysophosphatidylcholine and free fatty acid molecules. We found that a decrease in pH results in lateral squeezing of the membrane, i.e. in decreased surface area per headgroup. Thus, at the decreased pH, the lipid hydrocarbon chains had larger S(CD) order parameter values, and also enhanced membrane thickness, as seen in the electron density profiles across the membrane. From the lateral pressure profiles, we found that the values of spontaneous curvature of the two opposing monolayers became negative when the pH was decreased. At low pH, protonation of the free fatty acid headgroups reduces their mutual repulsion and accounts for the pH dependence of all the above-mentioned properties. The altered structural characteristics may significantly affect the overall surface properties of biomembranes in cellular vesicles, lipid droplets, and plasma lipoproteins, play an important role in membrane fission and fusion, and modify interactions between membrane lipids and the proteins embedded within them.

Atomistic simulations of phosphatidylcholines and cholesteryl esters in high-density lipoprotein-sized lipid droplet and trilayer: clues to cholesteryl ester transport and storage.

Cholesteryl esters (CEs) are the water-insoluble transport and storage form of cholesterol. For both transport and storage, phospholipids and proteins embrace the CEs to form an amphipathic monolayer that surrounds the CEs. CEs are transported extracellularly in lipoproteins and are stored intracellularly as cytoplasmic lipid droplets. To clarify the molecular phenomena related to the above structures, we conducted atomistic molecular-dynamics simulations for a spherical, approximately high density lipoprotein sized lipid droplet comprised of palmitoyl-oleoyl-phosphatidylcholine (POPC) and cholesteryl oleate (CO) molecules. An additional simulation was conducted for a lamellar lipid trilayer consisting of the same lipid constituents. The density profiles showed that COs were located in the core of the spherical droplet. In trilayer simulations, CO molecules were also in the core and formed two denser strata. This is remarkable because the intra- and intermolecular behaviors of the COs were similar to previous findings from bulk COs in the fluid phase. In accordance with previous experimental studies, the solubility of COs in the POPC monolayers was found to be low. The orientation distribution of the sterol moiety with respect to the normal of the system was found to be broad, with mainly isotropic or slightly parallel orientations observed deep in the core of the lipid droplet or the trilayer, respectively. In both systems, the orientation of the sterol moiety changed to perpendicular with respect to the normal close to the phopsholipid monolayers. Of interest, within the POPC monolayers, the intramolecular conformation of the COs varied from the previously proposed horseshoe-like conformation to a more extended one. From a metabolic point of view, the observed solubilization of CEs into the phospholipid monolayers, and the conformation of CEs in the phospholipid monolayers are likely to be important regulatory factors of CE transport and hydrolysis.

Atom-scale molecular interactions in lipid raft mixtures.

We review the relationship between molecular interactions and the properties of lipid environments. A specific focus is given on bilayers which contain sphingomyelin (SM) and sterols due to their essential role for the formation of lipid rafts. The discussion is based on recent atom-scale molecular dynamics simulations, complemented by extensive comparison to experimental data. The discussion is divided into four sections. The first part investigates the properties of one-component SM bilayers and compares them to bilayers with phosphatidylcholine (PC), the focus being on a detailed analysis of the hydrogen bonding network in the two bilayers. The second part deals with binary mixtures of sterols with either SM or PC. The results show how the membrane properties may vary substantially depending on the sterol and SM type available, the membrane order and interdigitation being just two of the many examples of this issue. The third part concentrates on the specificity of intermolecular interactions in three-component mixtures of SM, PC and cholesterol (CHOL) under conditions where the concentrations of SM and CHOL are dilute with respect to that of PC. The results show how SM and CHOL favor one another, thus acting as nucleation sites for the formation of highly ordered nanosized domains. Finally, the fourth part discusses the large-scale properties of raft-like membrane environments and compares them to the properties of non-raft membranes. The differences turn out to be substantial. As a particularly intriguing example of this, the lateral pressure profiles of raft-like and non-raft systems indicate that the lipid composition of membrane domains may have a major impact on membrane protein activation.

Molecular dynamics simulations of a lipovitellin-derived amphiphilic beta-sheet homologous to apoB-100 beta-sheets at a hydrophobic decane-water interface.

Lipovitellin, an egg-yolk lipoprotein, transports lipids in a pocket surrounded by amphiphilic beta-sheets. Its X-ray structure provides possibilities to study interactions between lipophilic beta-sheets and lipids at the atomic level. Here, we studied a 67-residue-long amphiphilic beta-sheet of lipovitellin previously suggested a suitable working model for studies of the lipid-binding behaviour of amphiphilic beta-sheet regions in apolipoprotein B-100 (apoB-100). We performed four molecular dynamics simulations with different starting configurations to define characteristics of the amphiphilic beta-sheet model at a decane-water interface. In each simulation the model beta-sheet bound keenly to the decane layer via its hydrophobic surface. The structural profiles showed unchanged secondary structure of the beta-sheet during the attachment. Also, aromatic side chains, especially tryptophans and tyrosines, mediated the attachment to the hydrophobic layer and influenced the orientation of the decane molecules that are in contact with the beta-sheet. In conclusion, the present simulations reveal high affinity of a lipovitellin-derived amphiphilic beta-sheet to a hydrophobic decane layer. They lay thereby the basis for further studies of the interaction between amphiphilic beta-sheets and lipids in complex molecular systems, like LDL particles, in which the large apoB-100 is the main protein component.

Polyunsaturation in lipid membranes: dynamic properties and lateral pressure profiles.

We elucidate the influence of unsaturation on single-component membrane properties, focusing on their dynamical aspects and lateral pressure profiles across the membrane. To this end, we employ atomistic molecular dynamics simulations to study five different membrane systems with varying degrees of unsaturation, starting from saturated membranes and systematically increasing the level of unsaturation, ending up with a bilayer of phospholipids containing the docosahexaenoic acid. For an increasing level of unsaturation, we find considerable effects on dynamical properties, such as accelerated dynamics of the phosphocholine head groups and glycerol backbones and speeded up rotational dynamics of the lipid molecules. The lateral pressure profile is found to be altered by the degree of unsaturation. For an increasing number of double bonds, the peak in the middle of the bilayer decreases. This is compensated for by changes in the membrane-water interface region in terms of increasing peak heights of the lateral pressure profile. Implications of the findings are briefly discussed.

Assessing the nature of lipid raft membranes.

The paradigm of biological membranes has recently gone through a major update. Instead of being fluid and homogeneous, recent studies suggest that membranes are characterized by transient domains with varying fluidity. In particular, a number of experimental studies have revealed the existence of highly ordered lateral domains rich in sphingomyelin and cholesterol (CHOL). These domains, called functional lipid rafts, have been suggested to take part in a variety of dynamic cellular processes such as membrane trafficking, signal transduction, and regulation of the activity of membrane proteins. However, despite the proposed importance of these domains, their properties, and even the precise nature of the lipid phases, have remained open issues mainly because the associated short time and length scales have posed a major challenge to experiments. In this work, we employ extensive atom-scale simulations to elucidate the properties of ternary raft mixtures with CHOL, palmitoylsphingomyelin (PSM), and palmitoyloleoylphosphatidylcholine. We simulate two bilayers of 1,024 lipids for 100 ns in the liquid-ordered phase and one system of the same size in the liquid-disordered phase. The studies provide evidence that the presence of PSM and CHOL in raft-like membranes leads to strongly packed and rigid bilayers. We also find that the simulated raft bilayers are characterized by nanoscale lateral heterogeneity, though the slow lateral diffusion renders the interpretation of the observed lateral heterogeneity more difficult. The findings reveal aspects of the role of favored (specific) lipid-lipid interactions within rafts and clarify the prominent role of CHOL in altering the properties of the membrane locally in its neighborhood. Also, we show that the presence of PSM and CHOL in rafts leads to intriguing lateral pressure profiles that are distinctly different from corresponding profiles in nonraft-like membranes. The results propose that the functioning of certain classes of membrane proteins is regulated by changes in the lateral pressure profile, which can be altered by a change in lipid content.

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.

Atomistic simulation studies of cholesteryl oleates: model for the core of lipoprotein particles.

We have conducted molecular dynamics simulations to gain insight into the atomic-scale properties of an isotropic system of cholesteryl oleate (CO) molecules. Cholesteryl esters are major constituents of low density lipoprotein particles, the key players in the formation of atherosclerosis, as well as the storage form of cholesterol. Here the aim is to clarify structural and dynamical properties of CO molecules under conditions, which are suggestive of those in the core of low density lipoprotein particles. The simulations in the fluid phase indicate that the system of CO molecules is characterized by an absence of translational order, as expected, while the orientational order between distinct CO molecules is significant at short distances, persisting over a molecular size. As for intramolecular properties, the bonds along the oleate chain are observed to be weakly ordered with respect to the sterol structure, unlike the bonds along the short hydrocarbon chain of cholesterol where the ordering is significant. The orientational distribution of the oleate chain as a whole with respect to the sterol moiety is of broad nature, having a major amount of extended and a less considerable proportion of bended structures. Distinct transient peaks at specific angles also appear. The diffusion of CO molecules is found to be a slow process and characterized by a diffusion coefficient of the order of 2x10(-9) cm2/s. This is considerably slower than diffusion, e.g., in ordered domains of lipid membranes rich in sphingomyelin and cholesterol. Analysis of the rotational diffusion rates and trans-to-gauche transition rates yield results consistent with experiments.

Influence of chain length and unsaturation on sphingomyelin bilayers.

Sphingomyelins (SMs) are among the most common phospholipid components of plasma membranes, usually constituting a mixture of several molecular species with various fatty acyl chain moieties. In this work, we utilize atomistic molecular dynamics simulations to study the differences in structural and dynamical properties of bilayers comprised of the most common natural SM species. Keeping the sphingosine moiety unchanged, we vary the amide bonded acyl chain from 16 to 24 carbons in length and examine the effect of unsaturation by comparing lipids with saturated and monounsaturated chains. As for structural properties, we find a slight decrease in average area per lipid and a clear linear increase in bilayer thickness with increasing acyl chain length both in saturated and unsaturated systems. Increasing the acyl chain length is found to further the interdigitation across the bilayer center. This is related to the dynamics of SM molecules, as the lateral diffusion rates decrease slightly for an increasing acyl chain length. Interdigitation also plays a role in interleaflet friction, which is stronger for unsaturated chains. The effect of the cis double bond is most significant on the local order parameters and rotation rates of the chains, though unsaturation shows global effects on overall lipid packing and dynamics as well. Regarding hydrogen bonding or properties related to the lipid/water interface region, no significant effects were observed due to varying chain length or unsaturation. The significance of the findings presented is discussed.

Response to comment by Almeida et al.: free area theories for lipid bilayers--predictive or not?

Free area theories for lateral diffusion in lipid bilayers are reviewed and discussed. It has been suggested by Almeida et al. that free area theories yield quantitative predictions for lateral diffusion coefficients of lipids. We investigate the plausibility of this suggestion by first sketching what is to be expected of a quantitative theory with predictive power, and subsequently examining whether existing free area theories comply with these expectations. Our conclusion is that current free area theories for lipid bilayers are not quantitative theories with predictive power. They involve a number of adjustable parameters, all of which are not estimated independently, but derived from fitting the theory to the very data whose behavior the theory is supposed to predict. Further, the interpretation and behavior of some of the parameters are ambiguous. The best example is the so-called activation barrier, whose qualitative behavior with the cholesterol concentration in a DMPC bilayer varies depending on the experimental method used to generate the input data and the exact assumptions made to formulate the theory. Independent determination of the activation barrier from numerical simulations or experiments appears to be very difficult.

Molecular dynamics simulations of unsaturated lipid bilayers: effects of varying the numbers of double bonds.

The importance of unsaturated, and especially polyunsaturated phosphatidylcholine molecules for the functional properties of biological membranes is widely accepted. Here, the effects of unsaturation on the nanosecond-scale structural and dynamic properties of the phosphatidylcholine bilayer were elucidated by performance of multinanosecond molecular dynamics simulations of all-atom bilayer models. Bilayers of dipalmitoylphosphatidylcholine and its mono-, di-, and tetraunsaturated counterparts were simulated, containing, respectively, oleoyl, linoleoyl, or arachidonoyl chains in the sn-2 position. Analysis of the simulations focused on comparison of the structural properties, especially the ordering of the chains in the membranes. Although the results suggest some problems in the CHARMM force field of the lipids when applied in a constant pressure ensemble, the features appearing in the ordering of the unsaturated chains are consistent with the behaviour known from (2)H NMR experiments. The rigidity of the double bonds is compensated by the flexibility of skew state single bonds juxtaposed with double bonds. The presence of double bonds in the sn-2 chains considerably reduces the order parameters of the CH bonds. Moreover, the double bond region of tetraunsaturated chains is shown to span all the way from the bilayer centre to the head group region.

Impact of cholesterol on voids in phospholipid membranes.

Free volume pockets or voids are important to many biological processes in cell membranes. Free volume fluctuations are a prerequisite for diffusion of lipids and other macromolecules in lipid bilayers. Permeation of small solutes across a membrane, as well as diffusion of solutes in the membrane interior are further examples of phenomena where voids and their properties play a central role. Cholesterol has been suggested to change the structure and function of membranes by altering their free volume properties. We study the effect of cholesterol on the properties of voids in dipalmitoylphosphatidylcholine (DPPC) bilayers by means of atomistic molecular dynamics simulations. We find that an increasing cholesterol concentration reduces the total amount of free volume in a bilayer. The effect of cholesterol on individual voids is most prominent in the region where the steroid ring structures of cholesterol molecules are located. Here a growing cholesterol content reduces the number of voids, completely removing voids of the size of a cholesterol molecule. The voids also become more elongated. The broad orientational distribution of voids observed in pure DPPC is, with a 30% molar concentration of cholesterol, replaced by a distribution where orientation along the bilayer normal is favored. Our results suggest that instead of being uniformly distributed to the whole bilayer, these effects are localized to the close vicinity of cholesterol molecules.

Lessons of slicing membranes: interplay of packing, free area, and lateral diffusion in phospholipid/cholesterol bilayers.

We employ 100-ns molecular dynamics simulations to study the influence of cholesterol on structural and dynamic properties of dipalmitoylphosphatidylcholine bilayers in the fluid phase. The effects of the cholesterol content on the bilayer structure are considered by varying the cholesterol concentration between 0 and 50%. We concentrate on the free area in the membrane and investigate quantities that are likely to be affected by changes in the free area and free volume properties. It is found that cholesterol has a strong impact on the free area properties of the bilayer. The changes in the amount of free area are shown to be intimately related to alterations in molecular packing, ordering of phospholipid tails, and behavior of compressibility moduli. Also the behavior of the lateral diffusion of both dipalmitoylphosphatidylcholine and cholesterol molecules with an increasing amount of cholesterol can in part be understood in terms of free area. Summarizing, our results highlight the central role of free area in comprehending the structural and dynamic properties of membranes containing cholesterol.

Structure and dynamics of sphingomyelin bilayer: insight gained through systematic comparison to phosphatidylcholine.

Sphingomyelin, one of the main lipid components of biological membranes, is actively involved in various cellular processes such as protein trafficking and signal transduction. In particular, specific lateral domains enriched in sphingomyelin and cholesterol have been proposed to play an important functional role in biomembranes, although their precise characteristics have remained unclear. A thorough understanding of the functional role of membranes requires detailed knowledge of their individual lipid components. Here, we employ molecular dynamics simulations to conduct a systematic comparison of a palmitoylsphingomyelin (PSM, 16:0-SM) bilayer with a membrane that comprises dipalmitoylphosphatidylcholine (DPPC) above the main phase transition temperature. We clarify atomic-scale properties that are specific to sphingomyelin due to its sphingosine moiety, and further discuss their implications for SM-rich membranes. We find that PSM bilayers, and in particular the dynamics of PSM systems, are distinctly different from those of a DPPC bilayer. When compared with DPPC, the strong hydrogen bonding properties characteristic to PSM are observed to lead to considerable structural changes in the polar headgroup and interface regions. The strong ordering of PSM acyl chains and specific ordering effects in the vicinity of a PSM-water interface reflect this issue clearly. The sphingosine moiety and related hydrogen bonding further play a crucial role in the dynamics of PSM bilayers, as most dynamic properties, such as lateral and rotational diffusion, are strongly suppressed. This is most evident in the rotational motion characterized by spin-lattice relaxation times and the decay of hydrogen bond autocorrelation functions that are expected to be important in complexation of SM with other lipids in many-component bilayers. A thorough understanding of SM bilayers would greatly benefit from nuclear magnetic resonance experiments for acyl chain ordering and dynamics, allowing full comparison of these simulations to experiments.