Overall Simulation analyses: The MD trajectories were analyzed with tools included in the GROMACS 2020.2 package and by in-house Python scripts. Molecular clusters of lipids were calculated with the gmx clustsize tool of GROMACS using a cut-off of 0.24nm (unless mentioned otherwise in the text), the largest distance to be considered in a cluster. The cut-off choice was based on the database of atomic van der Waals radii used by GROMACS . The equilibrium minimum pairwise distance between two lipid molecules is the equilibrium nucleus-nucleus distance of two hydrogen atoms i.e. 0.24nm (twice the van der Waals radius of a hydrogen atom). Visualization of molecular conformations was made with PyMOL (The PyMOL Molecular Graphics System, Version 2.3.0 Schrödinger, LLC.) and VMD, version 1.9.4a43 (June 11, 2020) . Further analyses and graphs were generated using Python.