Re: [AMBER] deltaG binding obtained from AmberTools 12 and 13 are different

From: Xioling Chuang <xioling_c.hotmail.com>
Date: Tue, 3 Sep 2013 11:04:50 +0200

Dear Jason,

Sorry for incomplete information. I used MMPBSA.py and I checked the energy components. I found that the energy difference obtained from AT12 and AT13 is vibrational energy (from nmode), not GB energy. The vibrational energy were
----------------------------------------
            AT12 AT13
----------------------------------------
cpx 540.8 524.6
receptor 303.3 295.5
ligand 243.9 247.9
delta -34.4 -46.9
----------------------------------------

The input is

&general
   startframe=1, endframe=1,
   interval=1, keep_files=0,
   verbose=1,
/
&gb
  igb=2, saltcon=0.0,
/
&nmode
   nmode_igb=1, nmode_istrng=0.0,
   drms=0.1, maxcyc=1000,
/

Any suggestion?

Thank you,
Xioling


> Date: Mon, 2 Sep 2013 12:12:23 -0400
> From: jason.swails.gmail.com
> To: amber.ambermd.org
> Subject: Re: [AMBER] deltaG binding obtained from AmberTools 12 and 13 are different
>
> On Mon, Sep 2, 2013 at 9:18 AM, Xioling Chuang <xioling_c.hotmail.com>wrote:
>
> > Dear All,
> >
> > I have used MM/GBSA to calculate the delta G binding of DNA with Amber12
> > and AmberTools12. Next, I installed AmberTools13 and re-compiled Amber12.
> > After re-calculating the delta G binding using the same input file (igb=2,
> > other parameters were default), the obtained result (including NMODE) are
> > different
> > AmberTools12: -70 kcal/mol
> > AmberTools13: -58 kcal/mol
> >
> > Why did two versions provide different result? Which one is reasonable?
> >
>
> We don't have enough information. I'm assuming you used MMPBSA.py. In any
> case, the GB code has remained almost identical, but there have been some
> small changes to the standard SA calculation using LCPO. (You can use
> molsurf in AmberTools 13 for GB/SA, which is a better option, IMO).
>
> It would be instructive to look at exactly which terms change between
> AmberTools 13 and AmberTools 12 to properly diagnose the problem. Without
> this information I can't help isolate the problem.
>
> HTH,
> Jason
>
> --
> Jason M. Swails
> BioMaPS,
> Rutgers University
> Postdoctoral Researcher
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
                                               
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Sep 03 2013 - 02:30:02 PDT
Custom Search