The lateral pressure proiles (LPPs) were calculated using the GROMACS-LS tool. The calculation of the LPPs was based on the last 200 ns of the simulation data (Table ), with the coordinates and velocities saved every 10 ps. An electrostatic cutoff of 2.0 nm was used for the LPP calculation. The proiles were calculated with about 140 slabs, corresponding to an approximate slab width of 1 Å. The outputs of GROMACS-LS were the components of the stress tensor, σ, as a function of the z coordinate (bilayer normal). LPP was calculated using the formula