AMBER: TI-FEP for ALA --> GLY

From: Jiten <jiten.postech.ac.kr>
Date: Sat, 25 Jun 2005 14:07:03 +0900

Dear Amber users

For my protein-drug system : while doing the free energy peturbation using TI method using 12 windows, and using the perl script attached here, I got almost zero difference in free energy. I used the following atom types, where DH is dummy atom with zero charge and the charges of the alaline is replaced by glycine charges in the pert.charge options


  H1 HC H0 DH
    | | | |
--CT---CT--HC ---> --CT --H0--DH
             | |
            HC DH

My frcmod.DH is as follows,

         ALA to GLY
MASS
DH 0.000

BOND
H0-DH 340.0 1.090 !CT-HC

ANGLE
DH - H0 - DH 35.0 109.50 !HC-CT-HC
CT - H0 - DH 35.0 109.50 !HC-CT-HC

DIHE
- The default values as given by leap 
--
--
NONB
DH     0.000    0.00
Am I doing something stupid here ? Do I need to set the Bond and angles : K and R values as zero ?
However, I calculate the absolute free energy of the two syatems using MM-PBSA from the last 1ns trajectory of the 2ns MD simulations (well converged - rmsd, temp and energy) and I got that there is clear differences between the two protein systems.
      ELE    -8.99 -8.96 
      VDW    -41.125 -37.68 
      GAS    -50.115 -46.64 
      PBSUR  -3.025 -3.07 
      PBCAL  27.105 24.16 
      PBSOL  24.09 21.09 
      PBELE  18.12 15.2 
      PBTOT  -26.025 -25.55 
      TSTOT 12.675 15.79 
      PBTOT-S  -13.35 -9.76 
Thanks for any suggestions,
Sincerely,
N. Jiten Singh
C/O Prof. Kwang S. Kim
Department of Chemistry
Pohang University of Science and Technology
San 31, Hyojadong, Namgu
Pohang 790-784, Korea
Phone : 82-54-279-5853 ( Lab ) / 279-4138 ( Appt )
Fax : 82-54-279-8137 (or +82-54-279-3399) 
Web : http://csm.postech.ac.kr/ and http://www.postech.ac.kr/e
Home Page : http://www.geocities.com/njs_19

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu

Received on Sat Jun 25 2005 - 06:53:00 PDT
Custom Search