Atomistic Model for Nearly Quantitative Simulations of Langmuir Monolayers.

Lung surfactant and a tear film lipid layer are examples of biologically relevant macromolecular structures found at the air-water interface. Because of their complexity, they are often studied in terms of simplified lipid layers, the simplest example being a Langmuir monolayer. Given the profound biological significance of these lipid assemblies, there is a need to understand their structure and dynamics on the nanoscale, yet there are not many techniques able to provide this information. Atomistic molecular dynamics simulations would be a tool fit for this purpose; however, the simulation models suggested until now have been qualitative instead of quantitative. This limitation has mainly stemmed from the challenge to correctly describe the surface tension of water with simulation parameters compatible with other biomolecules. In this work, we show that this limitation can be overcome by using the recently introduced four-point OPC water model, whose surface tension for water is demonstrated to be quantitatively consistent with experimental data and which is also shown to be compatible with the commonly employed lipid models. We further establish that the approach of combining the OPC four-point water model with the CHARMM36 lipid force field provides nearly quantitative agreement with experiments for the surface pressure-area isotherm for POPC and DPPC monolayers, also including the experimentally observed phase coexistence in a DPPC monolayer. The simulation models reported in this work pave the way for nearly quantitative atomistic studies of lipid-rich biological structures at air-water interfaces.

Freezing of stressed bilayers and vesicles.

We show using a minimalistic theoretical framework that phase transition decoupling in lipid bilayers is caused by a nonuniform stress profile due to an asymmetric distribution of lipids between the two leaflets. Applying this framework to vesicles, we demonstrate that their anomalous freezing is also caused by a stress asymmetry, but that this is due to lipid tail extension on freezing. Finally, we predict a previously unknown dependence of surface tension on temperature, find the phenomenon also in coarse grained molecular dynamics simulations, and suggest that it might have relevance in thermosensitive protein gating.

Interfacial tension and surface pressure of high density lipoprotein, low density lipoprotein, and related lipid droplets.

Lipid droplets play a central role in energy storage and metabolism on a cellular scale. Their core is comprised of hydrophobic lipids covered by a surface region consisting of amphiphilic lipids and proteins. For example, high and low density lipoproteins (HDL and LDL, respectively) are essentially lipid droplets surrounded by specific proteins, their main function being to transport cholesterol. Interfacial tension and surface pressure of these particles are of great interest because they are related to the shape and the stability of the droplets and to protein adsorption at the interface. Here we use coarse-grained molecular-dynamics simulations to consider a number of related issues by calculating the interfacial tension in protein-free lipid droplets, and in HDL and LDL particles mimicking physiological conditions. First, our results suggest that the curvature dependence of interfacial tension becomes significant for particles with a radius of ∼5 nm, when the area per molecule in the surface region is <1.4 nm(2). Further, interfacial tensions in the used HDL and LDL models are essentially unaffected by single apo-proteins at the surface. Finally, interfacial tensions of lipoproteins are higher than in thermodynamically stable droplets, suggesting that HDL and LDL are kinetically trapped into a metastable state.

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.

Coarse-grained computational studies of supported bilayers: current problems and their root causes.

The supported bilayer is an important experimental method for probing the features of lipid bilayers, yet relatively few computational studies have taken up its modeling. Using coarse-grained molecular dynamics simulations based on the MARTINI force field, we show that some of the few previous attempts employing similar methodologies are problematic, and have indeed led to erroneous conclusions. We further build a theoretical framework to see where exactly the failures of parametrization lie and suggest how these earlier results should be interpreted. The lessons learned may also prove useful in areas other than the supported bilayer as there is a current trend in modifying MARTINI, and other force fields, to suit a wide variety of specific applications. Finally, we analyze the polarizable water version of MARTINI and, although finding it imperfect in some of the same regards as the standard MARTINI, suggest using it in conjunction with supported bilayers.