Phospholipid bilayers of pure composition were prepared using the CHARMM-GUI The CHARMM36 force field was used for all reference phospholipid simulations . Each bilayer consists of 200 lipid molecules with a water layer of 2.25 nm on top and bottom of the bilayer. For the anionic PG and PS lipid systems, 200 sodium counter ions were added to the solution but otherwise, the systems were simulated without salt. The systems were first minimized and equilibrated in a series of small equilibrations, with gradually weaker restraints on protein and lipids using generic CHARMM-GUI simulation settings. Production simulations were 500 ns long at 310 K (POPC, POPE, POPG, POPS, PAPC, PDPC, PLC18), 323 K (PSM, SSM) or 333 K (DMPCE, DHPCE) and 1 bar using the v-rescale thermostat (τT = 1 ps) and Parrinello-Rahman barostat (τp = 5 ps) using semi-isotropic pressure coupling (Bussi, 2007). The non-bonded settings were identical to the hydrocarbon simulations. Simulations were performed using GROMACS v2019.4.