I'm trying to calculate binding energy of md simulation trajectory by g_mmpbsa. But van der Waals forces are very high. What should i do to lower the vander waal forces?
Mainly it is caused by either large conformational changes or broken molecule across the periodic box. If MD simulation with poorly optimized ligands, the calculated energies may not be matched with traditional values. Also check the grid box coordination, the vicinity of catalytic residues to backbone atoms of your ligand, etc.