|
|||||||||||||||||||||||||||||||||
AMBER Archive (2005)Subject: AMBER: PBSA residual and pairwise decompostion
From: Jiten (jiten_at_postech.ac.kr)
Dear all,
I tried to calculated the residue and pairwise decomposition using DCTYPE =2 and 4 options. I have notice some irregularity in the result. I tried to find the free energy contribution for the drug binding to the protein active site (It is drug-258, only ligand residue in the complex as in the enclosed excel file).
1. In pairwise data all the B-terms (backbone) are displayed as 0 (zero), while it shows some number in the case of residue decompostion data. Despite this the TVDW, TELE and TGAS values are the same in both the cases.
2. The most drastic differences are in GB calculations. The residue decompposition calculation seems to be alright as the total of the sum of the energy values of the residue is same with of the free energy binding calculated using MM-PBGB.
Therefore, I wonder if the energy values reported for pairwiase decomposition may be right or wrong.
Your valuable suggestions/corrections on this regard will be highly appreciated.
Sincerely,
N. Jiten Singh
-----------------------------------------------------------------------
| |||||||||||||||||||||||||||||||||
|