In the optimization of the cholesterol model, we ensured that the energetic and dynamic properties of the original model are fully maintained. The four additional beads in the constraint topology incur a computational cost. On the other hand, stricter LINCS settings required for the original Martini 2 cholesterol model also impact the computational cost.