AMBER Archive (2005)

Subject: RE: AMBER: An atom name with 4 characters is not compatible with Leap ?

From: Arnoult, Eric [PRDFR] (EArnoult_at_prdfr.JNJ.com)
Date: Wed Jul 27 2005 - 10:59:43 CDT


Dear Dr. Case,

I construct the pdb file with the following command lines :

${MOL}=test
${CHARGE}=0
${ATMTYPE}=gaff

${AMBHOME}/antechamber -i ${PREPDIR}/${MOL}.mol2 -fi mol2 -o ${MOL}.ac -fo
ac -j 5 -c bcc -nc ${CHARGE} -at ${ATMTYPE} -s 2 -rn ${MOL}
${AMBHOME}/antechamber -i ${MOL}.ac -fi ac -o ${MOL}_prp_resp.in -fo prepi
${AMBHOME}/antechamber -i ${MOL}.ac -fi ac -o ${MOL}.mpdb -fo mpdb -j 5 -s 2
-rn ${MOL}
${AMBHOME}/antechamber -i ${MOL}.ac -fi ac -c wc -cf ${MOL}.crg
${AMBHOME}/parmchk -i ${MOL}_prp_resp.in -f prepi -o frcmod${MOL}.par

May be I do a mistake in these command lines, I don't think so

thanks

éric

-----Original Message-----
From: owner-amber_at_scripps.edu [mailto:owner-amber_at_scripps.edu]On Behalf
Of David A. Case
Sent: Wednesday, July 27, 2005 5:08 PM
To: amber_at_scripps.edu
Subject: Re: AMBER: An atom name with 4 characters is not compatible
with Leap ?

On Wed, Jul 27, 2005, Arnoult, Eric [PRDFR] wrote:
>
> Can you try this leap.in ?
> "source leaprc.ff03
> source leaprc.gaff
> loadAmberPrep test_prp_resp.in
> lig = loadPdb test.mpdb
> saveAmberParm lig test.top test.crd
> quit"
>

The problem is with the pdb file, where the atom name for at 25 is in the
wrong columns:

ATOM 25 Cl25tes 1 9.155 ...

should be:

ATOM 25 Cl25 tes 1 9.155 ...

(Atom names are always in columns 13-16 of a pdb record.)

Notes:

1. After doing the loadAmberPrep step, you can type "desc tes.1", which will
show you all the atoms in the unit, with the four letter name correctly
there.

2. If you do a "savePdb tes foo.pdb", you will get a pdb file where the atom
names are in the correct column.

So the real question is: how did you construct the pdb file? There may be
some place that is not correctly handling the four letter atom name in the
pdb file, but I can't say more until I know how you created this file.

...good luck...dac

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