[AMBER] Problem with Decomposition of Binding Free Energy

From: loue chen <lauimi.883.gmail.com>
Date: Sat, 17 May 2014 00:08:46 +0500

*Dear Amber Experts,*

I am running MD simulation for a protein ligand complex to determine
selectivity. When I calculate the binding free energy, the MMGBS value is
according to the published data. But when I perform energy decomposition
using MMGBSA.py (Amber Tools 13) to calculate per residue contribution, all
of the energy shift difference shifts toward ligand, and the contribution
of amino acids become opposite.

I am using following input command for MMGBSA and decomposition. Please
check it weather there is some problem in my input command or I am missing
something. Or please if possible send me input command files for MMGBSA,
and Decomposition. Mine are given below:



*MMGBSA.in*

Input file for running PB and GB

&general

   verbose=1, interval=5,

   strip_mask=":WAT:Cl-:CIO:Cs+:IB:K+:Li+:Na+:Rb+",

/

&gb

  igb=5, saltcon=0.100

/

&pb

  istrng=0.100,

/
..........................................................
*DECOMP.in*
Per-residue GB and PB decomposition
&general
   verbose=1,
   #search_path=1,
   interval=10,
   strip_mask=":WAT:Cl-:CIO:Cs+:IB:K+:Li+:Na+:Rb+"
/
&gb
  igb=5, saltcon=0.100,
/
&pb
  istrng=0.100, inp=1, radiopt=0,
/
&decomp
  idecomp=1,
  dec_verbose=3,
/
-- 
*Regards,*
*Tahir Ali Chohan*
*B.Pharm., M.Phil (Pharm. Chem.)*
*PhD Scholar*
*College of Pharmaceutical Sciences*
*Zhejiang University, Hangzhou, China.*
*Cell # 0086-13018996850*
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri May 16 2014 - 12:30:02 PDT
Custom Search