Molecular dynamics simulations were performed using the LAMMPS software . For these simulations we used the Universal Force Field (UFF) for valence interactions (bond, angle, and dihedral terms) combined with the RexPoN universal nonbond potentials (UNB) to describe van der Waals interactions, and the polarizable charge equilibrium (PQEq) scheme for electrostatics.