AMBER Archive (2004)

Subject: AMBER: Possible amber8 antechamber PDB parsing bug

From: Mark Williamson (mark.williamson_at_imperial.ac.uk)
Date: Thu May 06 2004 - 10:16:19 CDT

Dear amber users,

I have noticed a possible bug with amber8's antechamber. It manifests
itself in antechamber's PDB output when charge fitting using g03 output.
The gaussian output file (2.3 Meg) used in this example can be found at
http://dumb.ch.ic.ac.uk/~mjw99/share/mcg_b3lyp_6-31G*_opt.log

The command line used is:

$AMBERHOME/exe/antechamber -i mcg_b3lyp_6-31G\*_opt.log -fi gout -o
mcg.prep -c resp -fo prepi -nc 1

Upon examining the PDB output (NEWPDB.PDB):

With amber7:
-----------------

ATOM 1 C8 MOL 1 -2.429 0.672 -0.560 -0.11095
ATOM 2 H6 MOL 1 -2.362 1.640 -1.044 0.12412
ATOM 3 C9 MOL 1 -3.632 0.008 -0.576 -0.14637

With amber8:
-----------------

ATOM 1 C8 MOL 1 -2.429 0.672 -0.560 -0.11095
ATOM 2 H6 MOL 1 -2.362 1.640 -1.044 0.12412
ATOM 3 C9 MOL 1 -3.632 0.008 -0.576 -0.14637

The problems is with spacing (this may be lost here due to email
formatting), the 5 column is moved right by one space (26 --> 27) . This
causes viewers like VMD to fail. According to the PDB file format (
http://www.rcsb.org/pdb/docs/format/pdbguide2.2/part_62.html ):

23 - 26 Integer resSeq Residue sequence number.

27 AChar iCode Code for insertion of
residues.

This movement is significant since it is at a column boundary; resSeq is
being expressed as iCode. This has also been attempted amber8 with
patches 1-8 applied, with the same result.
Can anyone confirm the above problem?

regards,

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