AMBER Archive (2005)Subject: Re: AMBER: new tutorial available
From: Ilyas Yildirim (yildirim_at_pas.rochester.edu)
Date: Wed Jan 26 2005  05:59:16 CST
Dear Mr. Case,
Thanks for the new tutorial. I have a couple of questions to ask.
1. The transformation is from toluene > nothing, and in order to do that,
u do a 2step MD simulation; one called as 'charging' and the other as
'disappearing'. For charging, the free energy difference is calculated as
2.35 kcal/mol, and for disappearing, the free energy difference is
calculated as (for klambda=6) 2.7 kcal/mol. So, the total free energy
difference of toluene > nothing transformation is 2.35 + 2.7 = 5.05
kcal/mol, am I right?
2. For the charging case, u are doing a 3point Gaussian Integration,
while for disappearing case, u are doing a 14point Gaussian Integration.
Is there any optimum number of lambda values for different structures at
interest? For instance if I do a simple perturbation on the amino group of
an RNA molecule, how can I decide how many lambda values I should use?
3. For the disappearing case with klambda=6, you have chosen lambda values
of 0.90, 0.92, 0.94, 0.96 and 0.98. Any reason for that?
4. In the frcmod file which describes the parameters for the dummy atoms,
you did not use zero values for harmonic force constants of the bonds and
angles, and the barrier heights (PK). If the atoms are disappearing, why
not choosing them as zero?
Thanks in advance. Best,
On Thu, 20 Jan 2005, David A. Case wrote:
> I have (finally!) assembled a tutorial on thermodynamic integration using
> sander. Follow the "Tutorials (English)" link at the Amber web site.
>
> Note that this is an "advanced" application: you will need to be familiar with
> running other Amber simulations before attempting this one.
>
> Comments, corrections and suggestions are welcome.
>
> ...enjoy...dac
>
> 
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
>

Ilyas Yildirim

 Department of Chemisty  
 University of Rochester  
 Hutchison Hall, # B10  
 Rochester, NY 146270216  Ph.:(585) 275 67 66 (Office) 
 http://www.pas.rochester.edu/~yildirim/ 


The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
