Re: [AMBER] How to interpret MMGBSA ligand binding free energy decomposition results

From: Carlos Simmerling <carlos.simmerling.gmail.com>
Date: Wed, 9 Mar 2016 06:13:24 -0500

All residues have van der waals interactions. Polar ones may also have
stronger electrostatics.
On Mar 9, 2016 3:26 AM, "Zhang Marc" <marczhang_md.zoho.com> wrote:

> Hi Jason
>
>
>
> Thanks a lot for your explanation.
>
>
>
> The non-polar interaction I refereed is the value of van der waals
> component.
>
>
>
> The result show van der waals= -1.21 that could be difficult to
> understand, as Threonine is a polar residue.
>
>
>
> I do not know how it contributes to the van der waals interaction.
>
>
>
> Cheers
>
> Marc​
>
>
>
>
>
> ---- On Wed, 09 Mar 2016 03:44:13 +0100 Jason Swails &
> lt;jason.swails.gmail.com&gt;wrote ----
>
>
>
>
> On Tue, Mar 8, 2016 at 12:08 PM, Zhang Marc &lt;marczhang_md.zoho.com&gt;
> wrote:
>
>
>
> &gt;
>
> &gt;
>
> &gt; Dear experts in Amber users group
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; I performed a 20 ns MD simulations of a ligand-protein complex system.
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; The ligand binding free energy was calculated using MMGBSA.
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; The energy was decomposed onto residues in the binding pocket for
>
> &gt; identifying residues important for ligand-protein interactions.
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; Now I have difficulties to interpret the decomposition result, as it
>
> &gt; suggested Threonine 119 is critical for both polar &amp;amp; non-polar
>
> &gt; interactions.
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; Here is the part of the result:
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; Resides van der Waals Electrostatic Polar Solvation
>
> &gt; Non-Polar Solvation TOTAL
>
> &gt;
>
> &gt; Thr119 -1.21 -2.37 0.27
>
> &gt; -0.93 -4.23
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; But Threonine is a polar residue and it is seldomly contributes to
>
> &gt; hydrophobic interactions, as I was taught during biochemistry course.
>
> &gt;
>
> &gt;
>
> &gt;
>
> &gt; Could you please help me figure out where the problem is?
>
> &gt;
>
>
>
> The nonpolar solvation free energy term is a very simple SASA based term.
>
> Basically, you take your solvent accessible surface area, multiply it by a
>
> constant, and that gives you the "nonpolar" part of solvation.
>
>
>
> It's unclear here whether these results are the *binding* results (i.e.,
>
> the bound - unbound energies), or if they are the bound or unbound values.
>
> But assuming this is the binding free energy, what a negative non-polar
>
> solvation energy here means is that the surface area exposed to solvent for
>
> this threonine in the bound state is smaller than that in the unbound
>
> state. If this is near the ligand, that makes sense. In the unbound state
>
> it's interacting with solvent (large SASA). In the bound state, it's
>
> interacting with the ligand (small SASA).
>
>
>
> HTH,
>
> Jason
>
> _______________________________________________
>
> 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
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 09 2016 - 03:30:21 PST
Custom Search