Despite the vast amount of experimental and theoretical studies on the binding affinity of cations – especially the biologically relevant Na⁺ and Ca²⁺ – for phospholipid bilayers, there is no consensus in the literature. Here we show that by interpreting changes in the choline headgroup order parameters according to the ‘molecular electrometer’ concept [Seelig et al., Biochemistry, 1987, 26, 7535], one can directly compare the ion binding affinities between simulations and experiments. Our findings strongly support the view that in contrast to Ca²⁺ and other multivalent ions, Na⁺ and other monovalent ions (except Li⁺) do not specifically bind to phosphatidylcholine lipid bilayers at sub-molar concentrations. However, the Na⁺ binding affinity was overestimated by several molecular dynamics simulation models, resulting in artificially positively charged bilayers and exaggerated structural effects in the lipid headgroups. While qualitatively correct headgroup order parameter response was observed with Ca²⁺ binding in all the tested models, no model had sufficient quantitative accuracy to interpret the Ca²⁺:lipid stoichiometry or the induced atomistic resolution structural changes. All scientific contributions to this open collaboration work were made publicly, using nmrlipids.blogspot.fi as the main communication platform.

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. Cardiovascular diseases are the primary cause of death in western countries. One of the main causes is lipid accumulation and plaque formation on arterial walls, called atherosclerosis. The risk of being exposed to this condition is reduced by high levels of high density lipoprotein (HDL). The functionality of HDL has remained elusive, and even its structure is not well understood. Through extensive coarse-grained simulations, we have clarified the structure of the lipid droplet in HDL and elucidated its interactions with the apolipoprotein A-I (apoA-I) that surrounds the droplet. We have found that the structural and dynamic properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). As for apoA-I, we have observed it alter the overall structure of the lipid droplet close to the HDL-water interface, with prominent changes taking place for cholesterol and other polar lipids. The nature of lipid-protein interactions is most favorable for cholesterol. Our results reveal that not only hydrophobicity but also conformational entropy are the driving forces in the formation of HDL structure, suggesting how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.

Conversion of discoidal phospholipid (PL)-rich high density lipoprotein (HDL) to spheroidal cholesteryl ester-rich HDL is a central step in reverse cholesterol transport. A detailed understanding of this process and the atheroprotective role of apolipoprotein A-I (apoA-I) requires knowledge of the structure and dynamics of these various particles. This study, combining computation with experimentation, illuminates structural features of apoA-I allowing it to incorporate varying amounts of PL. Molecular dynamics simulated annealing of PL-rich HDL models containing unesterified cholesterol results in double belt structures with the same general saddle-shaped conformation of both our previous molecular dynamics simulations at 310 K and the x-ray structure of lipid-free apoA-I. Conversion from a discoidal to a saddle-shaped particle involves loss of helicity and formation of loops in opposing antiparallel parts of the double belt. During surface expansion caused by the temperature-jump step, the curved palmitoyloleoylphosphatidylcholine bilayer surfaces approach planarity. Relaxation back into saddle-shaped structures after cool down and equilibration further supports the saddle-shaped particle model. Our kinetic analyses of reconstituted particles demonstrate that PL-rich particles exist in discrete sizes corresponding to local energetic minima. Agreement of experimental and computational determinations of particle size/shape and apoA-I helicity provide additional support for the saddle-shaped particle model. Truncation experiments combined with simulations suggest that the N-terminal proline-rich domain of apoA-I influences the stability of PL-rich HDL particles. We propose that apoA-I incorporates increasing PL in the form of minimal surface bilayers through the incremental unwinding of an initially twisted saddle-shaped apoA-I double belt structure.

For several decades, the standard model for high density lipoprotein (HDL) particles reconstituted from apolipoprotein A-I (apoA-I) and phospholipid (apoA-I/HDL) has been a discoidal particle ∼100 Å in diameter and the thickness of a phospholipid bilayer. Recently, Wu et al. (Wu, Z., Gogonea, V., Lee, X., Wagner, M. A., Li, X. M., Huang, Y., Undurti, A., May, R. P., Haertlein, M., Moulin, M., Gutsche, I., Zaccai, G., Didonato, J. A., and Hazen, S. L. (2009) J. Biol. Chem. 284, 36605–36619) used small angle neutron scattering to develop a new model they termed double superhelix (DSH) apoA-I that is dramatically different from the standard model. Their model possesses an open helical shape that wraps around a prolate ellipsoidal type I hexagonal lyotropic liquid crystalline phase. Here, we used three independent approaches, molecular dynamics, EM tomography, and fluorescence resonance energy transfer spectroscopy (FRET) to assess the validity of the DSH model. (i) By using molecular dynamics, two different approaches, all-atom simulated annealing and coarse-grained simulation, show that initial ellipsoidal DSH particles rapidly collapse to discoidal bilayer structures. These results suggest that, compatible with current knowledge of lipid phase diagrams, apoA-I cannot stabilize hexagonal I phase particles of phospholipid. (ii) By using EM, two different approaches, negative stain and cryo-EM tomography, show that reconstituted apoA-I/HDL particles are discoidal in shape. (iii) By using FRET, reconstituted apoA-I/HDL particles show a 28–34-Å intermolecular separation between terminal domain residues 40 and 240, a distance that is incompatible with the dimensions of the DSH model. Therefore, we suggest that, although novel, the DSH model is energetically unfavorable and not likely to be correct. Rather, we conclude that all evidence supports the likelihood that reconstituted apoA-I/HDL particles, in general, are discoidal in shape.