Friday, November 15, 2019

NMRlipids IVb: Preparing the manuscript about PE & PG results

Thanks to the various contributors, significant progress in the NMRlipids IVb project concerning PE and PG lipids has been made after the previous post.

First of all, as a response to the first ToDo issue in the previous post, two different codes to calculate the order parameters from united atom simulations have been contributed and also thoroughly discussed (see discussion in the previous post and buildH program). This enabled the analysis of previously contributed united atom simulations and the results are now added into the manuscript. Also, significant amount of new data has been contributed and we have now a good collection of force fields for both PE and PG lipids.

The data is still suggesting that the CHARMM36 simulations capture the essential differences between PC, PE, PS, and PG headgroups. A figure of structural ensembles of different headgroups from CHARMM36 simulations is now contributed.
https://raw.githubusercontent.com/NMRLipids/NMRlipidsIVPEandPG/master/Figs/struct.png
Figure 1: Overlayed snapshots of a glycerol backbone and headgroup conformations of different headgroups from CHARMM36 simulations. Figure is made by Pavel Buslaev.
However, we need also the quantitative analysis of dihedral distributions to finalize the discussion. I think that we should also compare the P-N vector angle with respect to membrane normal between lipids with different headgroups using CHARMM36 simulations. Also, further ideas to characterize the structural ensembles of lipid headgroups in a comprehensible way are welcomed. I have opened an issue in GitHub for further discussions about this topic.

For calcium-binding affinity to PG containing membranes we still need more data. The current data is not fully converged, and I think that we need longer CHARMM36 simulations with the recent NBfix correction for calcium and counterions. I have opened an issue in GitHub with more details for further discussion.

I have opened also other issues in GitHub. All the issues are important, but the most critical for the progress of the manuscript are the above mentioned analysis of structural differences between different lipid headgroups and additional data for PG-calcium interactions. In addition to the issues listed in GitHub, there are also Todo points in the manuscript (tex file) and SI (tex file). For example, we need citations for the force fields in tables I-III, and to finish the method sections for MD simulation and NMR experiments.