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 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 () has become the most extensive publicly available collection of
molecular dynamics simulation trajectories of lipid bilayers.