Re: AMBER: PBCAL in mmpbsa

From: nlxc <nlxc2000.yahoo.com>
Date: Tue, 29 Jun 2004 17:29:51 -0700 (PDT)

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
Custom Search