Re: AMBER: how to get the potential energy of a protein from the system

From: David A. Case <case.scripps.edu>
Date: Thu, 31 Jul 2008 07:45:47 -0700

On Thu, Jul 31, 2008, TaoPaul wrote:
>
> I am a freshman in MD. As far as I know the resulting potential energy (
> "EPtot" in out file ) was that of the whole system rather than the protein
> if ran a MD of protien in explicit water. So I wonder to know how to
> extract the potential energy of the protein itself. Could anybody give some
> suggestions?

Most modern energy routines (including PME and generalized Born) are not
"pairwise decomposable": so they cannot be written as just a "protein part"
plus a "solvent part" plus an "interaction part".

There *are* some sorts of analyses where it makes sense to extract the
coordinates of the protein, and compute just its potential energy. The
MM-PBSA analysis option in Amber automates this process for you. But be aware
that this involves some additional approximations.

....dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo.scripps.edu
Received on Sun Aug 03 2008 - 06:07:31 PDT
Custom Search