The first annual NMRlipids workshop

Twelve contributors participated in the first annual NMRlipids workshop in Berlin from 15th to 17th of May 2019. The flexible schedule included six presentations and a significant fraction of time was devoted to discussions and working towards the goals of the workshop.

The workshop was opened by Samuli Ollila presenting the current status and open issues in the NMRlipids project (slides available here). The presentation was followed by a discussion to determine the main goals of the workshop:

1. Define a quantitative quality measure for the structural quality of a lipid bilayer simulation.
2. Create instructions on how to contribute data into the NMRlipids databank.
3. Create a universal code to visualize lipid structures in simulations (for example, see this figure).

The second day of the workshop was opened by Ivan Gushchin presenting a principal component analysis of lipid structures (slides available from here). This was followed by Patrick Fuchs presenting the status of analysis code to calculate order parameters from united atom simulations, and Tiago Ferreira presenting order parameter experiments using natural abundance ¹³C NMR. At the end of the second day,  Stephanie Dawson presented publishing possibilities offered by www.scienceopen.com. The third day of the workshop was opened by Ivan Gushchin presenting the relation between lipid bilayer densities and scattering form factors (slides available from here).

The main outcomes of the workshop are listed here (topics 1-3 correspond to the goals listed above and 4-7 are additional outcomes):

1. A quantitative measure for the structural quality of lipid bilayer simulations. The aim was to define a simple machine and human readable measure for the structural quality of lipid bilayer in simulations using scattering form factors and order parameters. The main applications would be the NMRlipids III project (see post and post and issue) and the lipid databank. For order parameters, the proposal was to measure the distance from experimental values using the size of the error bar as a unit (see the first page in the document). For the form factor, the task turned out to be more complicated and universal quality measure was not yet defined (see the pages 1-2 in the document).
2. Instructions to contribute data into the NMRlipids databank. The aim was to generate straightforward instructions on how to contribute data into the NMRlipids project in the way to enable automatic analysis and indexing. The suggestion was to update the current indexing system and read the metadata directly from hashtags given in Zenodo descriptions. Preliminary script to pull data from Zenodo was created, but it is not publicly available yet. I have opened an issue for further discussion.
3. Visualize structural differences between force fields and lipids. The aim was to create a code which automatically visualizes lipid structure by generating figures like this, and analyze the differences between PC, PS, PG and PE headgroup structures as suggested in the NMRlipids IVb project. The code for this was generated, but it is not yet publicly available.
4. The analysis code to calculate order parameters from united atom models by Patrick Fuchs is progressing, but it is not yet publicly available. The merging with the current order parameters calculation code was superficially discussed. See also the discussion in the blog after the workshop.
5. In the discussion following the presentation by Tiago Ferreira, we concluded that the sn-1 chain order parameters from natural abundance ¹³C NMR may be less accurate than 0.02 (the previously used quantitative accuracy of order parameters) due to the spectral overlap. The order parameters from specifically deuterated samples are probably more accurate for this region. This is relevant for the quantitative quality measure and for the NMRlipids III project. Further discussion is in the issue.  
6. In the discussion following the presentation by Ivan Gushchin about lipid bilayer density profiles and scattering form factors, we concluded that some of the complications in the NMRlipids III project could be potentially resolved by separately analyzing the effect of POPC and cholesterol densities to the form factors from POPC/cholesterol mixture. 
7. The possibilities of using Markov state modeling presented by Ivan Gushchin (slide 5) to understand lipid structures and transition rates were discussed. 

The workshop was overall successful and the second NMRlipids workshop will be most likely organized at least partially with the same concept.

Explanation how to measure order parameters with ¹³C NMR