Berendsen barostat with a 1ps pressure relaxation time, and Langevin temperature control with a collision frequency of 1 ps -1 . The NPT simulations were again performed using a 1 fs simulation time step, and for all simulations the SHAKE algorithm 81 was applied to restrain all bonds to hydrogen atoms.