Re: [AMBER] In MMPBSY.py , how to set interval value when large conformation change ?

From: windy <chxp_moon.163.com>
Date: Wed, 22 Apr 2015 09:58:22 +0800 (CST)

Hi, Jason:

     Thanks for you reply . I finaly using the same frames to calculate the free energy.

     What you reply that we are using MMPBSA.py to compute a binding free energy for a noncovalent protein-ligand complex, and that you only care about that binding free energy.

      But my system is a convalent protein-ligand complex, may I use the MMPBSA.py to calculate the free energy ? Is the result reasonable enough to make others trust on?

      What I need to be care when using MMPBSA.py to calculate a convalent protein-ligand complex ?

      Excuse me for ask so many problem , haha....








在2015年04月21 02时47分, "Jason Swails"<jason.swails.gmail.com>写道:

On Mon, Apr 20, 2015 at 2:20 PM, George Tzotzos <gtzotzos.me.com> wrote:

> Jason,
>
> As I’ve come across similar situations, I’d like to chip in.
>
> In situations like this, would it help to do (i) clustering of the
> trajectory, and (ii) running MMPBSY.py on each cluster?
>

Again, it depends on what you hope to learn, and how you "run MMPBSA.py on
each cluster". The points I made about over-sampling lesser-populated
clusters by, e.g., taking the same number of snapshots from each cluster
remains (depending, of course, on the clustering algorithm you used).

Random sampling (or periodic sampling of a random, stationary distribution)
should give you the correct ensemble statistics (to a confidence determined
by the size of the subsample). Trying to cleverly use the ensemble to
determine how to sample that ensemble is just asking to introduce bias into
your binding free energy results, in my opinion.

What I'm saying corresponds to an assumption that you are using MMPBSA.py
to compute a binding free energy for a noncovalent protein-ligand complex,
and that you only care about that binding free energy. If you are instead
interested in the relative stabilities between two conformations (and you
don't feel that they are properly converged so you can't simply use their
populations), then some of these clustering-based approaches may be
appropriate.

All the best,
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 Apr 21 2015 - 19:00:03 PDT
Custom Search