- Contemporary messages sorted: [ by date ] [ by thread ] [ by subject ] [ by author ] [ by messages with attachments ]

From: Jason Swails <jason.swails.gmail.com>

Date: Wed, 9 Feb 2011 21:46:25 -0500

Hi Yokota,

As the GB equations consist of pairwise, analytical terms, you're right that

they're much easier to decompose than PB energies.

Since I haven't looked at the PBSA/decomp code, I can't say any of this for

sure. If I'm mistaken here, I'm hoping that a PBSA developer will step in

and correct me.

PBSA works by placing your system in a grid and solving the PB equation via

a finite difference method. It assigns each grid point a specific charge,

dielectric constant, and *starting* electrostatic potential (and ionic

strength if you specified one), which it then solves iteratively to get the

converged electrostatic potential. Thus, at the end of the calculation,

each gridpoint has associated with it a charge, potential, and dielectric.

You can use geometrical arguments to partition the grid points into

different residues (i.e. each grid point belongs to its *closest* residue).

Once you've created this partition, you simply toss it into the equation

DELTA G = 1/2 integral(charge density * reaction field potential)

Since it's a discrete problem, the integral becomes a sum, and for each

residue you just sum over that residue's grid points. I agree that it's not

as straightforward as GB, but it seems reasonable to me.

Hope this helps (and that it's right :) ),

Jason

On Mon, Feb 7, 2011 at 2:54 AM, <Yokota_Akihiro.takeda.co.jp> wrote:

*> Dear AMBER users,
*

*>
*

*> I have a question about "Polar Solvation" energy for each residue in
*

*> MMPBSA.py tutorial section 3.6.
*

*> http://ambermd.org/tutorials/advanced/tutorial3/py_script/section6.htm
*

*>
*

*> In the middle of the file below, there is a comment "Energy Decomposition
*

*> Analysis (All units kcal/mol):
*

*> Poisson Boltzmann solvent", I think this means that "Polar Solvation(PB)"
*

*> of "each residue" can be calculated.
*

*>
*

*> http://ambermd.org/tutorials/advanced/tutorial3/py_script/files/FINAL_DECOMP_MMPBSA_perres.dat
*

*>
*

*> However PB energy cannot be decomposed straightly by the definition (On the
*

*> other hand, GB energies can be decomposed).
*

*> Does anyone teach me how the PB energies are decomposed into each residue?
*

*>
*

*> Thanks in advance,
*

*> Akihiro Yokota
*

*>
*

*> _______________________________________________
*

*> AMBER mailing list
*

*> AMBER.ambermd.org
*

*> http://lists.ambermd.org/mailman/listinfo/amber
*

*>
*

Date: Wed, 9 Feb 2011 21:46:25 -0500

Hi Yokota,

As the GB equations consist of pairwise, analytical terms, you're right that

they're much easier to decompose than PB energies.

Since I haven't looked at the PBSA/decomp code, I can't say any of this for

sure. If I'm mistaken here, I'm hoping that a PBSA developer will step in

and correct me.

PBSA works by placing your system in a grid and solving the PB equation via

a finite difference method. It assigns each grid point a specific charge,

dielectric constant, and *starting* electrostatic potential (and ionic

strength if you specified one), which it then solves iteratively to get the

converged electrostatic potential. Thus, at the end of the calculation,

each gridpoint has associated with it a charge, potential, and dielectric.

You can use geometrical arguments to partition the grid points into

different residues (i.e. each grid point belongs to its *closest* residue).

Once you've created this partition, you simply toss it into the equation

DELTA G = 1/2 integral(charge density * reaction field potential)

Since it's a discrete problem, the integral becomes a sum, and for each

residue you just sum over that residue's grid points. I agree that it's not

as straightforward as GB, but it seems reasonable to me.

Hope this helps (and that it's right :) ),

Jason

On Mon, Feb 7, 2011 at 2:54 AM, <Yokota_Akihiro.takeda.co.jp> wrote:

-- Jason M. Swails Quantum Theory Project, University of Florida Ph.D. Graduate Student 352-392-4032 _______________________________________________ AMBER mailing list AMBER.ambermd.org http://lists.ambermd.org/mailman/listinfo/amberReceived on Wed Feb 09 2011 - 19:00:02 PST

Custom Search