I'm using protocols like "Abalone\Scripts\standard protocols\optimization by dynamics (B).js" (in the Abalone program) for the local optimization and the replica exchange for the global optimization. Amber99SB or Amber03 force field, always with a Generalized Born method, and often with reduced integer charges.