Hello, Amber users
javascript:SetCmd(cmdSend);
I'm trying to calculate relative binding free energies of two drug molecules using the TI facilities in Amber10. The perturbation is an O-atom to a NH-group. I've followed the tutorial "small molecule binding to T4-lysozyme" and therefore breaked up the perturbation in 3 steps: 1) remove charge of O-atom, 2) transformation of neutral O-atom to neutral NH-group with soft-core potential and 3) addition of charge to NH-group. I have used 9 lambda values from 0.1 to 0.9. The experimental difference in binding energy is very small, in the order of 1 kcal/mole. The 1:st and the 3:rd steps seemed to work okey since the deltaG are on the order of 0.1 kcal/mole. However the problem is in the 2:nd step where I get unreasonable high dv/dl values, both in the protein and in the water phase. Here the dv/dl are about 15-30 kcal/mole and deltaG evaluates to approx. -20kcal/mol. Here is an example of the input file for the production run:
Production
&cntrl
irest=1,ntx=5,
nstlim=1000000,dt=0.002,
temp0=300.0,ntt=3,gamma_ln=2.0,
ntc=2,ntf=1,
cut=8.0,
ntpr=5000,ntwx=5000,ntwv=0,ntwe=5000,iwrap=1,
ntb=1,
ipol=0,igb=0,
scnb=2.0,scee=1.2,
ntr=0,
icfe=1, clambda = 0.5,
ifsc=1,
crgmask=':CBB.032',
scmask=':CBB.O32',
&end
I hope someone can shed some light on my problem, and maybe point to possible sources of error.
Best regards,
Samuel
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Sun Jun 29 2008 - 06:07:07 PDT
