Re: [AMBER] Interpreting energy decomposition in Amber thermodynamic integration

From: Stefan Ivanov <stefan.ivanov.postgrad.manchester.ac.uk>
Date: Sun, 22 Jan 2017 10:31:05 +0000

Hello Prof Case,

Once again, thank you very much for your reply.

What I meant was "What are the reported per-residue internal, VDW, and electrostatic energies?" Are they also averages of dV/dlambda for the given lambda value? If so, I would need to integrate the per-residue values between 0 and 1 to get meaningful numbers, correct?

Sincerely,

Stefan

________________________________________
From: David Case [david.case.rutgers.edu]
Sent: Sunday, January 22, 2017 3:05 AM
To: AMBER Mailing List
Subject: Re: [AMBER] Interpreting energy decomposition in Amber thermodynamic integration

On Fri, Jan 20, 2017, Stefan Ivanov wrote:
>
> I'm transforming a template ligand bound to a protein to a modeled
> ligand bound to the same protein with 9 windows (for starters). Say I
> look at the computed per-residue energies at = 0.1 and at = 0.9. I
> want to know what the reported energies mean at the specified values
> and how I can compare them, say for residue X at = 0.1 and at = 0.9.

You ask "what the reported energies mean": they are the average values
of dV/dlambda at the given lambda point. Only the integral of dV/dlambda has
any usual "meaning": it is a free energy.
>
>
> As I said, I'm doing decomposition with idecomp = 1 and only the ligands
> have any significant internal energy contributions and I was wondering
> why that is.

Ah: the word "significant" is important here. You are changing the ligand, so
the protein energies will only change to the extent that the protein
configuration differs between the two ligands. Unless residue 129 (say) is in
a different environment with ligand A vs ligand B, it will not contribute
to the free energy change. (One reason that I am not much of fan of TI energy
decomposition: it often tells you no more than what you could have figured
out with a careful visual study of the protein binding pocket.)

I encourage you to experiment with some very simple TI transformations to
get a clearer understanding of how everything works.

...dac



_______________________________________________
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 Sun Jan 22 2017 - 03:00:02 PST
Custom Search