“…The molecular dynamic simulations were performed with gromacs5 packages (Abraham et al, 2015), with the charmm36 force field (Klauda et al, 2010) and the TIP3P water model (Jorgensen, Chandrasekhar, Madura, Impey, & Klein, 1983) used. The topology files, charges, and force field parameters of RIV were generated by online software charmm General Force Field (CGenFF; Vanommeslaeghe, Raman, & Mackerell, 2012;Vanommeslaeghe et al, 2009Vanommeslaeghe et al, , 2011. Each system was put in a SPC216 water box (Berendsen, Postma, van Gunsteren, & Hermans, 1981) where the distance from the edge of the box to the surface of the complex was at least 1 nm.…”