We primarily used the Martini 3 force field to run CG simulations. However, a number of parameters were adjusted to better capture atomistic behaviors of the system (see SI for the details of these changes). We adopted a timestep of 20 fs for the production run, which is appropriate since atomic bond vibrations are absent in these simulations. Some initial equilibration simulations required a shorter time step for stability. Periodic boundary conditions were implemented similar to all-atom MD. A cutoff of 1.0 nm with a potential-shift-Verlet modifier was applied for both van der Waals and Coulomb interactions. The Coulomb potential was modified using the reaction-field method, which assumes a constant dielectric environment beyond the cutoff distance. Temperature coupling used the velocity rescale thermostat 58 with tau-t = 1 ps and pressure coupling used the Berendsen algorithm with tau-p = 1 ps.