On Tue, Jul 25, 2006, Lwin, ThuZar wrote:
>
> I am learning how to do thermodynamic integration using amber8. I would
> like to mutate Alanine to Glycine by mutating the side chain CH3 to H.
> Here is what I did in tleap. Please let me know what I am doing wrong.
>
> a = sequence{ALA}
> set a.CB pertName HA3
The syntax is unit-name.residue.atom-name, e.g. "a.1.CB", not "a.CB".
However, I really recommend that you use xleap, visually select all atoms
in your structure, and use the "edit selected atoms" command. That will
give you a spreadsheet that shows which things need to be modified.
And, if you are thinking long-range, it would be good to update to Amber 9,
since the way in which thermodynamic integration is implemented is quite
different there.
...regards...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jul 30 2006 - 06:07:04 PDT
