The concentration of cholesterol in cell membranes affects membrane fluidity and thickness, and might regulate different processes such as the formation of lipid rafts. Since interpreting experimental data from biological membranes is rather intricate, investigations on simple models with biological relevance are necessary to understand the natural systems. We study the effect of cholesterol on the molecular structure of multi-lamellar vesicles (MLVs) composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a phospholipid ubiquitous in cell membranes, with compositions in the range 0-60 mol% cholesterol. Order parameters, |S(CH)|, are experimentally determined by using (1)H-(13)C solid-state nuclear magnetic resonance (NMR) spectroscopy with segmental detail for all parts of both the cholesterol and POPC molecules, namely the ring system and alkyl chain of the sterol, as well as the glycerol backbone, choline headgroup and the sn-1 and sn-2 acyl chains of POPC. With increasing cholesterol concentration the acyl chains gradually adopt a more extended conformation while the orientation and dynamics of the polar groups are rather unaffected. Additionally, we perform classical molecular dynamics simulations on virtual bilayers mimicking the POPC-cholesterol MLVs investigated by NMR. Good agreement between experiments and simulations is found for the cholesterol alignment in the bilayer and for the |S(CH)| profiles of acyl chains below 15 mol% cholesterol. Deviations occur for the choline headgroup and glycerol backbone parts of POPC, as well as for the phospholipid and cholesterol alkyl chains at higher cholesterol concentrations. The unprecedented detail of the NMR data enables a more complete comparison between simulations and experiments on POPC-cholesterol bilayers and may aid in developing more realistic model descriptions of biological membranes.

Phosphatidylserine (PS) is a negatively charged lipid type commonly found in eukaryotic membranes, where it interacts with proteins via nonspecific electrostatic interactions as well as via specific binding. Moreover, in the presence of calcium ions, PS lipids can induce membrane fusion and phase separation. Molecular details of these phenomena remain poorly understood, partly because accurate models to interpret the experimental data have not been available. Here we gather a set of previously published experimental NMR data of C-H bond order parameter magnitudes, |S CH |, for pure PS and mixed PS:PC (phosphatidylcholine) lipid bilayers, and augment this data set by measuring the signs of S CH in the PS headgroup using S-DROSS solid-state NMR spectroscopy. The augmented data set is then used to assess the accuracy of the PS headgroup structures in, and the cation binding to, PS-containing membranes in the most commonly used classical molecular dynamics (MD) force fields including CHARMM36, Lipid17, MacRog, Slipids, GROMOS-CKP, Berger, and variants. We show large discrepancies between different force fields, and that none of them reproduces the NMR data within experimental accuracy. However, the best MD models can detect the most essential differences between PC and PS headgroup structures. The cation binding affinity is, in line with our previous results for PC lipids, not captured correctly by any of the PS force fields. Moreover, the simulated response of PS headgroup to bound ions can differ from experiments even qualitatively. The collected experimental dataset and simulation results will pave the way for development of lipid force fields that correctly describe the biologically relevant negatively charged membranes and their interactions with ions. This work is part of the NMRlipids open collaboration project (nmrlipids.blogspot.fi).

Interest in lipid interactions with proteins and other biomolecules is emerging not only in fundamental biochemistry but also in the field of nanobiotechnology where lipids are commonly used, for example, in carriers of mRNA vaccines. The outward-facing components of cellular membranes and lipid nanoparticles, the lipid headgroups, regulate membrane interactions with approaching substances, such as proteins, drugs, RNA, or viruses. Because lipid headgroup conformational ensembles have not been experimentally determined in physiologically relevant conditions, an essential question about their interactions with other biomolecules remains unanswered: Do headgroups exchange between a few rigid structures, or fluctuate freely across a practically continuous spectrum of conformations? Here, we combine solid-state NMR experiments and molecular dynamics simulations from the NMRlipids Project to resolve the conformational ensembles of headgroups of four key lipid types in various biologically relevant conditions. We find that lipid headgroups sample a wide range of overlapping conformations in both neutral and charged cellular membranes, and that differences in the headgroup chemistry manifest only in probability distributions of conformations. Furthermore, the analysis of 894 protein-bound lipid structures from the Protein Data Bank suggests that lipids can bind to proteins in a wide range of conformations, which are not limited by the headgroup chemistry. We propose that lipids can select a suitable headgroup conformation from the wide range available to them to fit the various binding sites in proteins. The proposed inverse conformational selection model will extend also to lipid binding to targets other than proteins, such as drugs, RNA, and viruses.

Phosphatidylserine (PS) lipids are important signaling molecules and the most common negatively charged lipids in eukaryotic membranes. The signaling can be often regulated by calcium, but its interactions with PS headgroups are not fully understood. Classical molecular dynamics (MD) simulations can potentially give detailed description of lipid-ion interactions, but the results strongly depend on the used force eld. Here, we apply the electronic continuum correction (ECC) to the Amber Lipid17 parameters of 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine (POPS) lipid to improve its interactions with K + , Na + and Ca 2+ ions. The partial charges of headgroup, 1 glycerol backbone and carbonyls of POPS, bearing a unit negative charge, were scaled with a factor of 0.75, derived for monovalent ions and the Lennard-Jones σ parameters of the same segments were scaled with a factor of 0.89. The resulting ECC-POPS model gives more realistic interactions with Na + and Ca 2+ cations than the original Amber Lipid17 parameters, when validated using headgroup order parameters and the electrometer concept. In ECC-lipids simulations, populations of complexes of Ca 2+ cations with more than two PS lipids are negligible, and interactions of Ca 2+ cations with only carboxylate groups is twice more likely than with only phosphate groups, while interactions with carbonyls almost entirely involve also other groups. Our results pave the way for more realistic MD simulations of biomolecular systems with anionic membranes allowing to elucidate signaling processes involving PS and Ca 2+ .

Molecular dynamics (MD) simulations give atomically detailed information on structure and dynamics in amphiphilic bilayer systems on timescales up to about 1 μs. The reorientational dynamics of the C-H bonds is conventionally verified by measurements of (13)C or (2)H nuclear magnetic resonance (NMR) longitudinal relaxation rates R1, which are more sensitive to motional processes with correlation times close to the inverse Larmor frequency, typically around 1-10 ns on standard NMR instrumentation, and are thus less sensitive to the 10-1000 ns timescale motion that can be observed in the MD simulations. We propose an experimental procedure for atomically resolved model-free estimation of the C-H bond effective reorientational correlation time τe, which includes contributions from the entire range of all-atom MD timescales and that can be calculated directly from the MD trajectories. The approach is based on measurements of (13)C R1 and R1ρ relaxation rates, as well as (1)H-(13)C dipolar couplings, and is applicable to anisotropic liquid crystalline lipid or surfactant systems using a conventional solid-state NMR spectrometer and samples with natural isotopic composition. The procedure is demonstrated on a fully hydrated lamellar phase of 1-palmitoyl-2-oleoyl-phosphatidylcholine, yielding values of τe from 0.1 ns for the methyl groups in the choline moiety and at the end of the acyl chains to 3 ns for the g1 methylene group of the glycerol backbone. MD simulations performed with a widely used united-atom force-field reproduce the τe-profile of the major part of the acyl chains but underestimate the dynamics of the glycerol backbone and adjacent molecular segments. The measurement of experimental τe-profiles can be used to study subtle effects on C-H bond reorientational motions in anisotropic liquid crystals, as well as to validate the C-H bond reorientation dynamics predicted in MD simulations of amphiphilic bilayers such as lipid membranes.