In all steps, the LINCS algorithm was used to constrain all bonds containing hydrogen atoms and a time step of 2 fs was used as the integration time step.