All MD simulations were performed using the opensource package GROMACS All bonds were constrained using the LINCS algorithm. The SPC/E model was employed to describe intermolecular interactions between water molecules, while intramolecular and intermolecular interactions involving SDS and C 12 E 6 were modeled using a combination of the TRaPPE-UA and OPLS force fields described in Ref. 34 A short-range cut-off of 1.0 nm was employed for both the (shifted) Lennard-Jones and Coulomb interactions. Long-range electrostatics were modeled using the Particle Mesh Ewald-Switch method, and a long-range dispersion correction was used correct both the energy and pressure for the cutoff in the van der Waals interaction. A time step of 2 fs was used for the Verlet integrator. All simulations were performed in the canonical ensemble (NVT) at a fixed temperature at 298.15 K using the stochastic velocity rescaling of Bussi et. al. with a time constant of 0.1 ps.