I read and revised my whole workflow, i don't think I had missed something. My protein is metalloprotein containing two zn atoms in the a active siite. The docked ligands show positive mm_gbsa values with OPLS3e forcefield and the vaues became negative only when i use Opls2005 forcefeild. I couldn't understand what was going on.@zhaoxi Sun