Dear All,
I've been trying to calculate the binding free energy with MM-BPSA using
the mm_pbsa.pl script
implemented in amber 10. For the single trajectory all looks fine, the
binding energy looks
reasonable, however, for a multiple trajectory I've encountered a
strange problem.
I've run separately the MD for the ligand and receptor, and calculate
the MM-PBSA energies
in three different cases by considering the complex, the ligand, and
receptor only as a receptor,
accordingly to the method we have to follow.
The ligand and the complex gives similar energy components to the one
calculated for the single
trajectory case, however, the receptor has a very different values,
except for the Eint which is very
similar to the single traj case, that should indicate the geometry is ok
(all the STDs look fine too).
I've check few times and all seems consistent, the method the parameters
and the ff used.
I've tried two different systems and I got the same problem, even if I
use amber 9.
I would really appreciate if someone could point me to where a potential
error might be.
I'm sure I'm missing something.
thanks a lot
sincerely
claudio
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Sun Aug 31 2008 - 06:07:09 PDT
