Molecular dynamics trajectory was then
run using the minimized structure as the input, and periodic
boundary conditions and particle mesh Ewald electrostatics
(Essmann et al., 1995) were used in the simulation. SHAKE was
employed to keep all bonds involving hydrogen atoms rigid. The
time step of the simulations was 2.0 fs with a cutoff of 12 Å for
the non-bonded interaction. The General Amber Force Field
(GAFF) parameters were assigned to the ligand. The partial charges
were determined using the AM1-BCC method of the Antechamber
suite of AMBER 11.