Thomas,
Thanks a lot. Based on the curve of EPTOT vs. time and DENSITY vs.
time, my system looks well equilibrated. Do I have to run production MD
under NVT condition or I still can use NPT for my production run. Is there any
standard such as rmsd for how well the system is equilibrated?
Thanks a lot for help!
Nan
Thomas Steinbrecher <steinbrt.scripps.edu> wrote:Hi Nan,
>a little surprise to me since the PBCAL values for the same receptor in
>those two systems( through separate MD and mm/pbsa) have about 40 Kcal
>difference(std=22), and the ligand also has 1.5 difference(std=1). What
>could be the possible reasons to make the PB values of two exactly same
>receptors so different?
I have noted that even small differences in starting geometries
can result in quite different deltaG results from MM-PBSA.
This is probably only more true if you study two similar ligands instead
of two conformations of one ligand.
But still, 40kcal/mol difference sounds like something else is going
on. Are you sure your simulation is well equilibrated? You mentioned that
your crystal structure does not have a ligand bound to it, so maybe the
geometry of the binding site is far from equilibrium with the ligand
bound.
Regards,
Thomas
---------------------------------
Do you Yahoo!?
New and Improved Yahoo! Mail - 100MB free storage!
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jun 30 2004 - 01:53:01 PDT