AMBER Archive (2004)

Subject: RE: AMBER: (no subject)

From: FyD (fyd_at_u-picardie.fr)
Date: Wed Jul 14 2004 - 14:30:04 CDT


Hi,

> I am going through the plastocynin tutorial. Can you tell me
> how I can use
> GAFF of antechamber for this purpose, since i could not
> figure out how to
> use the resp charges generated by RED as there is no relevent
> file type which can be used.

R.E.D. generates a Tripos "Mol_o.mol2" file with the RESP/ESP charges in it.

To load this file into leap/amber8 (only), run:
xleap -f leaprc.ff99

M = loadmol2 Mol_o.mol2
saveamberparm M M.top M.crd

Before be sure to manually add the atom types needed according to the Cornell
et
al. paper, (JACS 117 5179. See Table 1 "List of atom types"), i.e. modify the
6th column of this .mol2 file:
cf http://www.tripos.com/custResources/mol2Files/.
Personnally, I prefer to "manually" add the atom types because this step is
crucial and determines your FF accuracy. As said Ross, you have to be aware of
what you do...

Using "saveamberparm M M.top M.crd" leap (1st run), it is going to tell you,
which FF params are missing. Then, you can create your personnal FRCMOD file
(using GAFF or by fitting ab-initio data) which contains 'your' NEW FF params
that can be load into leap (in a 2nd run) i.e. using:
xleap -f leaprc.ff99

M = loadmol2 Mol_o.mol2
FRCMOD = loadamberparam frcmod.Anshul
saveamberparm M M.top M.crd

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