Association of lipidome remodeling in the adipocyte membrane with acquired obesity in humans.

Identification of early mechanisms that may lead from obesity towards complications such as metabolic syndrome is of great interest. Here we performed lipidomic analyses of adipose tissue in twin pairs discordant for obesity but still metabolically compensated. In parallel we studied more evolved states of obesity by investigating a separated set of individuals considered to be morbidly obese. Despite lower dietary polyunsaturated fatty acid intake, the obese twin individuals had increased proportions of palmitoleic and arachidonic acids in their adipose tissue, including increased levels of ethanolamine plasmalogens containing arachidonic acid. Information gathered from these experimental groups was used for molecular dynamics simulations of lipid bilayers combined with dependency network analysis of combined clinical, lipidomics, and gene expression data. The simulations suggested that the observed lipid remodeling maintains the biophysical properties of lipid membranes, at the price, however, of increasing their vulnerability to inflammation. Conversely, in morbidly obese subjects, the proportion of plasmalogens containing arachidonic acid in the adipose tissue was markedly decreased. We also show by in vitro Elovl6 knockdown that the lipid network regulating the observed remodeling may be amenable to genetic modulation. Together, our novel approach suggests a physiological mechanism by which adaptation of adipocyte membranes to adipose tissue expansion associates with positive energy balance, potentially leading to higher vulnerability to inflammation in acquired obesity. Further studies will be needed to determine the cause of this effect.

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 proteins diffuse as dynamic complexes with lipids.

We describe how membrane proteins diffuse laterally in the membrane plane together with the lipids surrounding them. We find a number of intriguing phenomena. The lateral displacements of the protein and the lipids are strongly correlated, as the protein and the neighboring lipids form a dynamical protein-lipid complex, consisting of approximately 50-100 lipids. The diffusion of the lipids in the complex is much slower compared to the rest of the lipids. We also find a strong directional correlation between the movements of the protein and the lipids in its vicinity. The results imply that in crowded membrane environments there are no "free" lipids, as they are all influenced by the protein structure and dynamics. Our results indicate that, in studies of cell membranes, protein and lipid dynamics have to be considered together.

Composition and lipid spatial distribution of HDL particles in subjects with low and high HDL-cholesterol.

A low level of high density lipoprotein cholesterol (HDL-C) is a powerful risk factor for cardiovascular disease. However, despite the reported key role of apolipo-proteins, specifically, apoA-I, in HDL metabolism, lipid molecular composition of HDL particles in subjects with high and low HDL-C levels is currently unknown. Here lipidomics was used to study HDL derived from well-characterized high and low HDL-C subjects. Low HDL-C subjects had elevated triacylglycerols and diminished lysophosphatidylcholines and sphingomyelins. Using information about the lipid composition of HDL particles in these two groups, we reconstituted HDL particles in silico by performing large-scale molecular dynamics simulations. In addition to confirming the measured change in particle size, we found that the changes in lipid composition also induced specific spatial distributions of lipids within the HDL particles, including a higher amount of triacylglycerols at the surface of HDL particles in low HDL-C subjects. Our findings have important implications for understanding HDL metabolism and function. For the first time we demonstrate the power of combining molecular profiling of lipoproteins with dynamic modeling of lipoprotein structure.

Bioinformatics and computational methods for lipidomics.

Large amounts of lipidomics data are rapidly becoming available. However, there is a lack of tools capable of taking the full advantage of the wealth of new information. Lipid bioinformatics is thus an emerging need as well as challenge for lipid research. Lipid concentration changes in biological systems reflect regulation at multiple spatial and dynamic scales, e.g., biochemical reactions in the cells, intercellular lipid trafficking, changes in cell membrane composition, systemic lipid metabolism or lipid oxidation. In order to address the complexity of lipids and their regulation, four areas of bioinformatics need to be developed: (1) data processing and lipid identification, (2) statistical data analysis, (3) pathway analysis, and (4) lipid modeling in systems and biophysical contexts. In this paper we overview the current state of the lipid bioinformatics field as well as suggest few potential new areas of research.

Conformational changes and slow dynamics through microsecond polarized atomistic molecular simulation of an integral Kv1.2 ion channel.

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins-not to mention numerous applications in drug design. Here, we present a full 1 micros atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120 degrees rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation ( approximately 35 degrees ) of the extracellular end of all S4 segments is present also in a reference 0.5 micros simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 3(10) helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4-lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5-1 micros). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations.

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.

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.

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.