In MD simulations of phase separating ternary lipid mixtures with LINCS settings chosen to ensure similarly small temperature gradients, we achieved comparable performance with the original and optimized cholesterol model in terms of simulated time per wall-clock time (ns/day). A slight performance advantage of the optimized model increased further to ≈ 30% in presence of the membrane protein β 2 AR because the increase in the lincs order required for the original cholesterol model applies also to constraints in the membrane protein.