both Coulomb and Lennard-Jones potentials with particle-mesh Ewald (PME) correction for long-range electrostatics and a shift potential applied between 0.9 and 1.1 nm to make the Lennard-Jones potential converges smoothly to zero at the cut-off. Also, a relative dielectric constant ε r = 15 was used to attenuate the Coulomb potential, which is done in Martini force field to compensate the absence of molecular dipoles which would attenuate the interactions between ions. An integration time step of 0.02 ps was used in all simulations.