In the systems containing lipid bilayers, we computed the temperature of the different lipids using the Gromacs 18 tool gmx traj. For the different cholesterol models, the temperatures calculated from the kinetic energies were corrected for the respective numbers of free and constrained degrees of freedom. The standard Martini 2 cholesterol consisting of eight beads has 3N = 24 degrees of freedom. However, the three mass-less virtual sites do not contribute to the kinetic energy. Together with the five constraints imposed on the structure, they leave only 10 degrees of freedom. Hence, a correction of 24/10 was applied to the temperature values. The same reasoning results in a correction factor of 36/10 for our optimized model.