On Mon, Sep 9, 2013 at 5:49 AM, 千语千寻 <qianyuqianxun88.gmail.com> wrote:
> Hi,everyone
>
> I have used MMPBSA.py to calculate the binding erergy of a complex which
> includes a protein and a ligand .For a same complex,I use the same
> inputfilein sander and
> pmemd programme at the equilibrium stage.Then I use the mdcrd file to begin
> MMPBSA. But the results of GB are quite different. The difference value is
> 22 kcal/mol.
> Have you met this question? I have several questions
> 1, the process which i use is right?
> 2, what is the reason of the D-value?
>
You may not have run long enough to obtain converged results. The
trajectories of 2 simulations will never be identical with different
programs, but that doesn't make one of them 'wrong', just 'different'
(i.e., the two simulations could be exploring completely different regions
of phase space).
However, there is so little information here that we cannot possibly
provide specific help or advice. pmemd and sander have been so thoroughly
evaluated and checked that I find it highly unlikely that 'standard' MD in
either program is 'wrong'.
Also, you should not be analyzing the 'equilibration' dynamics, since the
equilibration stage is typically non-equilibrium. (It is probably better
termed the 'relaxation' stage) You should run at least 5 to 10 ns of
simulation AFTER a multi-nanosecond equilibration/relaxation stage before
attempting MM/PBSA.
Good luck,
Jason
--
Jason M. Swails
BioMaPS,
Rutgers University
Postdoctoral Researcher
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Sep 09 2013 - 10:00:02 PDT
