Effect of cholesterol on structural and mechanical properties of membranes depends on lipid chain saturation

Molecular simulation of the effect of cholesterol on lipid-mediated protein-protein interactions

Experiments and molecular simulations have shown that the hydrophobic mismatch between proteins and membranes contributes significantly to lipid-mediated protein-protein interactions. In this article, we discuss the effect of cholesterol on lipid-mediated protein-protein interactions as function of hydrophobic mismatch, protein diameter and protein cluster size, lipid tail length, and temperature. To do so, we study a mesoscopic model of a hydrated bilayer containing lipids and cholesterol in which proteins are embedded, with a hybrid dissipative particle dynamics-Monte Carlo method. We propose a mechanism by which cholesterol affects protein interactions: protein-induced, cholesterol-enriched, or cholesterol-depleted lipid shells surrounding the proteins affect the lipid-mediated protein-protein interactions. Our calculations of the potential of mean force between proteins and protein clusters show that the addition of cholesterol dramatically reduces repulsive lipid-mediated interactions between proteins (protein clusters) with positive mismatch, but does not affect attractive interactions between proteins with negative mismatch. Cholesterol has only a modest effect on the repulsive interactions between proteins with different mismatch.

Self-consistent mean-field model for palmitoyloleoylphosphatidylcholine-palmitoyl sphingomyelin-cholesterol lipid bilayers

The connection between membrane inhomogeneity and the structural basis of lipid rafts has sparked interest in the lateral organization of model lipid bilayers of two and three components. In an effort to investigate anisotropic lipid distribution in mixed bilayers, a self-consistent mean-field theoretical model is applied to palmitoyloleoylphosphatidylcholine (POPC)--palmitoyl sphingomyelin (PSM)--cholesterol mixtures. The compositional dependence of lateral organization in these mixtures is mapped onto a ternary plot. The model utilizes molecular dynamics simulations to estimate interaction parameters and to construct chain conformation libraries. We find that at some concentration ratios the bilayers separate spatially into regions of higher and lower chain order coinciding with areas enriched with PSM and POPC, respectively. To examine the effect of the asymmetric chain structure of POPC on bilayer lateral inhomogeneity, we consider POPC-lipid interactions with and without angular dependence. Results are compared with experimental data and with results from a similar model for mixtures of dioleoylphosphatidylcholine, steroyl sphingomyelin, and cholesterol.

Lipid metabolism in myelinating glial cells: lessons from human inherited disorders and mouse models

The integrity of central and peripheral nervous system myelin is affected in numerous lipid metabolism disorders. This vulnerability was so far mostly attributed to the extraordinarily high level of lipid synthesis that is required for the formation of myelin, and to the relative autonomy in lipid synthesis of myelinating glial cells because of blood barriers shielding the nervous system from circulating lipids. Recent insights from analysis of inherited lipid disorders, especially those with prevailing lipid depletion and from mouse models with glia-specific disruption of lipid metabolism, shed new light on this issue. The particular lipid composition of myelin, the transport of lipid-associated myelin proteins, and the necessity for timely assembly of the myelin sheath all contribute to the observed vulnerability of myelin to perturbed lipid metabolism. Furthermore, the uptake of external lipids may also play a role in the formation of myelin membranes. In addition to an improved understanding of basic myelin biology, these data provide a foundation for future therapeutic interventions aiming at preserving glial cell integrity in metabolic disorders.

Lipid chain branching at the iso- and anteiso-positions in complex Chlamydia membranes: a molecular dynamics study

Membranes in the intracellular eubacterial parasite Chlamydia trachomatis consist of the elementary body (EB) and reticular body (RB), and contain methyl branches at the iso- and anteiso-positions for some phospholipid chains. Acyl chain branching is the focus of this study. Molecular dynamics simulations were used to study bilayers of 1-13-methylpentadecanoyl-2-palmitoyl-phosphatidylcholine (13-MpPPC), 1-14-methylpentadecanoyl-2-palmitoyl-phosphatidylcholine (14-MpPPC), and diphytanoylphosphatidylcholine (DPhPC). These three membranes were simulated at 323K and simulations of DPhPC at 298K were also performed for better comparison to existing experimental data. Two simulations of representative EB and RB membranes of C. trachomatis composed of nine different lipid components were performed at 310.15K, to accurately reflect compositions determined by experiment and physiological conditions. Based on nearly 0.5μs of simulation data, we report that branching increases average lipid surface area, area elastic moduli, and lipid axial relaxation times, while decreasing lipid chain order. Branching also has a distinct effect on electron density profiles. Due to their high cholesterol concentrations, the EB and RB membranes were found to have relatively high area elastic moduli, which may have important biological implications.

HIV fusion peptide penetrates, disorders, and softens T-cell membrane mimics

This work investigates the interaction of N-terminal gp41 fusion peptide (FP) of human immunodeficiency virus type 1 (HIV-1) with model membranes in order to elucidate how FP leads to fusion of HIV and T-cell membranes. FP constructs were (i) wild-type FP23 (23 N-terminal amino acids of gp41), (ii) water-soluble monomeric FP that adds six lysines on the C-terminus of FP23 (FPwsm), and (iii) the C-terminus covalently linked trimeric version (FPtri) of FPwsm. Model membranes were (i) LM3 (a T-cell mimic), (ii) 1,2-dioleoyl-sn-glycero-3-phosphocholine, (iii) 1,2-dioleoyl-sn-glycero-3-phosphocholine/30 mol% cholesterol, (iv) 1,2-dierucoyl-sn-glycero-3-phosphocholine, and (v) 1,2-dierucoyl-sn-glycero-3-phosphocholine/30 mol% cholesterol. Diffuse synchrotron low-angle x-ray scattering from fully hydrated samples, supplemented by volumetric data, showed that FP23 and FPtri penetrate into the hydrocarbon region and cause membranes to thin. Depth of penetration appears to depend upon a complex combination of factors including bilayer thickness, presence of cholesterol, and electrostatics. X-ray data showed an increase in curvature in hexagonal phase 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine, which further indicates that FP23 penetrates into the hydrocarbon region rather than residing in the interfacial headgroup region. Low-angle x-ray scattering data also yielded the bending modulus K(C), a measure of membrane stiffness, and wide-angle x-ray scattering yielded the S(xray) orientational order parameter. Both FP23 and FPtri decreased K(C) and S(xray) considerably, while the weak effect of FPwsm suggests that it did not partition strongly into LM3 model membranes. Our results are consistent with the HIV FP disordering and softening the T-cell membrane, thereby lowering the activation energy for viral membrane fusion.

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.

Lateral organization of complex lipid mixtures from multiscale modeling

The organizational properties of complex lipid mixtures can give rise to functionally important structures in cell membranes. In model membranes, ternary lipid-cholesterol (CHOL) mixtures are often used as representative systems to investigate the formation and stabilization of localized structural domains ("rafts"). In this work, we describe a self-consistent mean-field model that builds on molecular dynamics simulations to incorporate multiple lipid components and to investigate the lateral organization of such mixtures. The model predictions reveal regions of bimodal order on ternary plots that are in good agreement with experiment. Specifically, we have applied the model to ternary mixtures composed of dioleoylphosphatidylcholine:18:0 sphingomyelin:CHOL. This work provides insight into the specific intermolecular interactions that drive the formation of localized domains in these mixtures. The model makes use of molecular dynamics simulations to extract interaction parameters and to provide chain configuration order parameter libraries.

Alamethicin aggregation in lipid membranes

X-ray scattering features induced by aggregates of alamethicin (Alm) were obtained in oriented stacks of model membranes of DOPC(diC18:1PC) and diC22:1PC. The first feature obtained near full hydration was Bragg rod in-plane scattering near 0.11 A(-1) in DOPC and near 0.08 A(-1) in diC22:1PC at a 1:10 Alm:lipid ratio. This feature is interpreted as bundles consisting of n Alm monomers in a barrel-stave configuration surrounding a water pore. Fitting the scattering data to previously published molecular dynamics simulations indicates that the number of peptides per bundle is n = 6 in DOPC and n >or= 9 in diC22:1PC. The larger bundle size in diC22:1PC is explained by hydrophobic mismatch of Alm with the thicker bilayer. A second diffuse scattering peak located at q(r) approximately 0.7 A(-1) is obtained for both DOPC and diC22:1PC at several peptide concentrations. Theoretical calculations indicate that this peak cannot be caused by the Alm bundle structure. Instead, we interpret it as being due to two-dimensional hexagonally packed clusters in equilibrium with Alm bundles. As the relative humidity was reduced, interactions between Alm in neighboring bilayers produced more peaks with three-dimensional crystallographic character that do not index with the conventional hexagonal space groups.