Headgroup Structure and Cation Binding in Phosphatidylserine Lipid Bilayers.

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, |SCH|, for pure PS and mixed PS:PC (phosphatidylcholine) lipid bilayers and augment this data set by measuring the signs of SCH 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 not captured correctly by any of the PS force fields-an observation that is in line with our previous results for PC lipids. Moreover, the simulated response of the PS headgroup to bound ions can differ from experiments even qualitatively. The collected experimental data set 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 ).

Shape and Phase Transitions in a PEGylated Phospholipid System.

Poly(ethylene glycol) (PEG) polymers and PEG-conjugated lipids are widely used in bioengineering and drug transport applications. A PEG layer in a drug carrier increases hydrophilic repulsion, inhibits membrane fusion and serum opsonin interactions, and prolongs the storage and circulation time. It can also change the carrier shape and have an influence on many properties related to the content release of the carrier. In this paper, we focus on the physicochemical effects of PEGylation in the lipid bilayer. We introduce laurdanC as a fluorophore for shape recognition and phase transition detection. Together with laurdanC, cryogenic transmission electron microscopy, differential scanning calorimetry, molecular dynamics simulations, and small-angle X-ray scattering/wide-angle X-ray scattering, we acquire information of the particle/bilayer morphology and phase behavior in systems containing 1,2-dipalmitoyl- sn-glycero-3-phosphocholine:1,2-distearoyl- sn-glycero-3-phosphoethanolamine-PEG(2000) with different fractions. We find that PEGylation leads to two important and potentially usable features of the system. (1) Spherical vesicles present a window of elevated chain-melting temperatures and (2) lipid packing shape-controlled liposome-to-bicelle transition. The first finding is significant for targets requiring multiple release sequences and the second enables tuning the release by composition and the PEG polymer length. Besides drug delivery systems, the findings can be used in other smart soft materials with trigger-polymers as well.

Simulations Study of Single-Component and Mixed n-Alkyl-PEG Micelles.

Here, we study one-component and mixed n-alkyl-poly(ethylene glycol) (C mE n) micelles with varying poly(ethylene glycol) (PEG) chain lengths n using coarse-grained molecular simulations. These nonionic alkyl-PEG surfactants and their aggregates are widely used in bio and chemical technology. As expected, the simulations show that increasing the PEG chain length decreases the alkyl-PEG micelle core diameter and the aggregation number but also enhances PEG chain penetration to the core region and spreads the micelle corona. Both the core and corona density are heavily dependent on the PEG chain length and decrease with increasing PEG length. Furthermore, we find that the alkyl-PEG surfactants exhibit two distinct micellization modes: surfactants with short PEG chains as their hydrophilic heads aggregate with the PEG heads relatively extended. Their aggregation number and the PEG corona density are dictated by the core carbon density. For longer PEG chains, the PEG sterics, that is, the volume occupied by the PEG head group, becomes the critical factor limiting the aggregation. Finally, simulations of binary mixtures of alkyl-PEGs of two different PEG chain lengths show that even in the absence of core-freezing, the surfactants prefer the aggregate size of their single-component solutions with the segregation propelled via enthalpic contributions. The findings, especially as they provide a handle on the density and the density profile of the aggregates, raise attention to effective packing shape as a design factor of micellar systems, for example, drug transport, solubilization, or partitioning.

Molecular electrometer and binding of cations to phospholipid bilayers.

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

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions.

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C-H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P-N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files ( https://zenodo.org/collection/user-nmrlipids ) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Controlling carbon-nanotube-phospholipid solubility by curvature-dependent self-assembly.

Control of aqueous dispersion is central in the processing and usage of nanoscale hydrophobic objects. However, selecting dispersive agents based on the size and form of the hydrophobic object and the role of coating morphology in dispersion efficiency remain important open questions. Here, the effect of the substrate and the dispersing molecule curvature, as well as, the influence of dispersant concentration on the adsorption morphology are examined by molecular simulations of graphene and carbon nanotube (CNT) substrates with phospholipids of varying curvature as the dispersing agents. Lipid spontaneous curvature is increased from close to zero (effectively cylindrical lipid) to highly positive (effectively conical lipid) by studying double tailed dipalmitoylphosphadidylcholine (DPPC) and single tailed lysophosphadidylcholine (LPC) which differ in the number of acyl chains but have identical headgroup. We find that lipids are good dispersion agents for both planar and curved nanoparticles and induce a dispersive barrier nonsize selectively. Differences in dispersion efficiency arise from lipid headgroup density and their extension from the hydrophobic substrate in the adsorption morphology. We map the packing morphology contributing factors and report that the aggregate morphologies depend on the competition of interactions rising from (1) hydrophobicity driven maximization of lipid-substrate contacts and lipid self-adhesion, (2) tail bending energy cost, (3) preferential alignment along the graphitic substrate principal axes, and (4) lipid headgroup preferential packing. Curved substrates adjust the morphology by changing the balance between the interaction strengths. Jointly, the findings show substrate curvature and dimensions are a way to tune lipid adsorption to desired, self-assembling patterns. Besides engineering dispersion efficiency, the findings could bear significance in designing materials with defined molecular scale, molecular coatings for orientation specific CNT assembly or lipid-based molecular masks and patterning on graphene.

Intrinsic α-helical and ß-sheet conformational preferences: a computational case study of alanine.

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, ß-sheet, or other backbone dihedral angle ( ϕ-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that ß-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) ß-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and ß-sheet structures only occur when the backbone bond angle τ (NCα C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the "bridge" region of ϕ-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between ϕ-ψ conformations and τ.