The systems were simulated with the Martini 3 force field . All systems were subjected to the same relaxation protocol. First a minimization of 500 steps using steepest descent Velocity-rescale thermostat with 𝜏 = 1 ps, coupling lipids separately from the solvent. The pressure was maintained with a semiisotropic pressure coupling during relaxation and production run using respectively the Berendsen [104] and Parrinello-Rahman [106] pressure coupled with 𝜏 = 12 ps and the compressibility at 3 × 10 -4 bar -1 . The verlet buffer tolerance was set to 0.002 kJ mol -1 ps -1 . The reaction field method was used to treat electrostatics and van der Waals interactions were truncated after 1.1 nm ("New-RF" Martini settings [110]) using the potential shift Verlet scheme [108].