CHARMM force field (Vanommeslaeghe et al., 2010[14]) and Gasteiger charges were used.