Currently, the manuscript is mainly a collection of relevant experimental data. For example, Fig. 1 compares the experimental headgroup and glycerol backbone order parameters between PC, PE, PG and PS lipids.
|Fig.1 Absolute values of order parameters for headgroup and glycerol|
backbone with different headgroups from experiments. For references and other details see the manuscript.
Some preliminary comparison between experiments and simulations with CHARMM GUI parameters are shown in Figs. 2 and 3, suggesting that the model has some difficulties to reproduce the experimental order parameters for PS and PG headgroups. More detailed conclusions are difficult to draw only from these data, because experimentally the signs of order parameters for PS and PG are not available (as far as I know). However, the results from other models might help to draw some connections between order parameters and structural details, as was done in NMRlipids I for PC lipids.
|Fig 2. Order parameters for POPS headgroup and glycerol backbone|
from simulations and experiments. For references and details see the manuscript. Absolute
values are shown for experimental data, because signs are not
known. Simulations values are -SCH
|Fig 3. Order parameters for PG headgroup and glycerol backbone|
from simulations and experiments. For references and details see the manuscript. Absolute values are shown, because signs are not known for experimental data.
|Fig. 6 Changes in the PC headgroup order parameter as a function of CaCl2 concentration in bilayers containing various amounts of negatively charged lipids. For references and details see the manuscript.|
"(i) Ca2+ binds to neutral lipids (phosphatidylcholine, phosphatidylethanolamine) and negatively charged lipids (phosphatidylglycerol) with approximately the same binding constant of K = 10-20 M-1;
(ii) the free Ca2+ concentration at the membrane interface is distinctly enhanced if the membrane carries a negative surface charge, either due to protein or to lipid;
(iii) increased interfacial Ca2+ also means increased amounts of bound Ca2+ at neutral and charged lipids;
(iv) the actual binding step can be described by a Langmuir adsorption isotherm with a 1_lipid_:_1_Ca2+ stoichiometry, provided the interfacial concentration CM is used to describe the chemical binding equilibrium."
I believe that an MD simulation model correctly reproducing the cation binding in negatively charged lipids could further sharpen this interpretation.
The goal of this project will be to test if currently available models can be used for an such interpretation. This should also help the model development (if needed), however the actual improvement of force fields is beyond the scope of NMRlipids IV.
As in all NMRlipids projects, all types of contributions (data, comments, criticism, etc.) are welcomed from everyone. The authorship of the publication will be offered to all contributors and the final acceptance is based on self-assessment according to NMRlipids rules. The following contributions would be especially relevant at this stage:
- Results from different simulation models. Simulations of bilayers containing PE, PG or PS almost under any conditions would be currently useful to map the behavior of different models. Direct delivery of calculated order parameters through GitHub or blog comments, or by making the simulation trajectories accessible (for example, through Zenodo) would be ideal ways of contributing.
- Order parameter signs for PE, PG and PS. The order parameter signs are very important for the structural interpretation. However, I am not aware of the order parameter sign measurements for other than PC lipids. If such data would be somehow available, this would be highly useful contribution.