We also performed MD simulations of ternary phase separating lipid mixtures based on the simulations of Thallmair et al. . Systems consisting of 1,276 dipalmitoyl-phosphatidylcholine (DPPC), 912 cholesterol, and 950 dilinoleoyl-phosphatidyl-choline (DLiPC) molecules corresponding to a molar ratio of 0.42/0.28/0.30 were built using the program insane.py in a random distribution. The bilayers were solvated in a 0.15 M NaCl solution. This resulted in an overall number of 38,917 CG water beads (10% of which were anti-freeze particles) and 428 NaCl ion pairs. All simulations were performed with Gromacs 18 version 2020.1. Energy minimization of the initial structures using a steepest descent algorithm was followed by two 500 ps pre-equilibration runs with time steps of 1 and 10 fs, respectively. The equilibration was concluded by a run of 1 µs, which is long enough for phase separation to occur. During production, we simulated 750,000,000 steps corresponding to 7.5 µs with ∆t = 10 fs, 15 µs with ∆t = 20 fs, and 22.5 µs with ∆t = 30 fs. The MD parameters and settings corresponded to the "New-RF" values. During the simulations, the pressure was maintained at 1 bar using a Parrinello-Rahman barostat with semi-isotropic coupling. We employed a coupling constant of τ p = 12.0 ps and a compressibility of β = 3 × 10 -4 bar -1 . The temperature was kept constant at 310 K using a velocity rescaling thermostat with a coupling constant of τ T = 1.0 ps. One thermostat was used for the solvent beads (water, antifreeze, and ion beads) and a second independent thermostat for the lipids.