We used the GROMACS insert-molecule function to add glycerol molecules. The solvation process was done iteratively as GROMACS insert-molecule function cannot insert all the glycerol molecules in a single run. Each iteration is followed by energy minimization and short NVT equilibrations (50 ps). Once no more glycerol molecules could (or need to) be added, we ran a longer NVT equilibration (1 ns) and a short NPT equilibration (50 ps) on the system. If the solvation level was below the target we attempted to add more glycerol molecules by repeating the iteration but with short NPT (instead of NVT) equilibration at the end of each iteration. Finally, we ran a longer NPT (5 ns) equilibration on the system.