Coordinates after the MM minimization step were used as starting to prevent a potential drift manifested due to the incorrect protonation and forcefield limitations. Then, 500 ps of NVT and NPT equilibration steps were performed while restraining heavy atoms including water oxygens by 100 000 kJ/mol/nm 2 . Water was treated as flexible for finer position adjustment.