Despite the vast amount of experimental and theoretical studies on the binding affinity of cations -especially the biologically relevant Na + and Ca 2+ -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 2+ 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 2+ binding in all the tested models, no model had sufficient quantitative accuracy to interpret the Ca 2+ :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.

ApoA-I is a uniquely flexible lipid-scavenging protein capable of incorporating phospholipids into stable particles. Here we report molecular dynamics simulations on a series of progressively smaller discoidal high density lipoprotein particles produced by incremental removal of palmitoyloleoylphosphatidylcholine via four different pathways. The starting model contained 160 palmitoyloleoylphosphatidylcholines and a belt of two antiparallel amphipathic helical lipid-associating domains of apolipoprotein (apo) A-I. The results are particularly compelling. After a few nanoseconds of molecular dynamics simulation, independent of the starting particle and method of size reduction, all simulated double belts of the four lipidated apoA-I particles have helical domains that impressively approximate the x-ray crystal structure of lipid-free apoA-I, particularly between residues 88 and 186. These results provide atomic resolution models for two of the particles produced by in vitro reconstitution of nascent high density lipoprotein particles. These particles, measuring 95 angstroms and 78 angstroms by nondenaturing gradient gel electrophoresis, correspond in composition and in size/shape (by negative stain electron microscopy) to the simulated particles with molar ratios of 100:2 and 50:2, respectively. The lipids of the 100:2 particle family form minimal surfaces at their monolayer-monolayer interface, whereas the 50:2 particle family displays a lipid pocket capable of binding a dynamic range of phospholipid molecules.

Spheroidal high-density lipoprotein (HDL) particles circulating in the blood are formed through an enzymatic process activated by apoA-I, leading to the esterification of cholesterol, which creates a hydrophobic core of cholesteryl ester molecules in the middle of the discoidal phospholipid bilayer. In this study, we investigated the conformation of apoA-I in model spheroidal HDL (ms-HDL) particles using both atomistic and coarse-grained molecular dynamics simulations, which are found to provide consistent results for all HDL properties we studied. The observed small contribution of cholesteryl oleate molecules to the solvent-accessible surface area of the entire ms-HDL particle indicates that palmitoyloleoylphosphatidylcholines and apoA-I molecules cover the hydrophobic core comprised of cholesteryl esters particularly well. The ms-HDL particles are found to form a prolate ellipsoidal shape, with sizes consistent with experimental results. Large rigid domains and low mobility of the protein are seen in all the simulations. Additionally, the average number of contacts of cholesteryl ester molecules with apoA-I residues indicates that cholesteryl esters interact with protein residues mainly through their cholesterol moiety. We propose that the interaction of annular cholesteryl oleate molecules contributes to apoA-I rigidity stabilizing and regulating the structure and function of the ms-HDL particle.