AMBER: GAFF atom typing for flavin

From: Hans Martin Senn <>
Date: Sat, 4 Aug 2007 17:15:30 +0100

Dear AMBER community

I am using antechamber/leap to generate topologies for largish
organic molecules (up to ~300 atoms) containing a flavin (FAD/FMN/
isoalloxazine) moiety. While this works perfectly for most parts of
the molecules, the automatic atom-type assignment for the
isoalloxazine ring system is completely unreasonable. In particular,
antechamber fails to recognize that the three rings together form an
aromatic system.
So I have been trying to assign the GAFF atom types manually, using
the definitions and examples in the GAFF and antechamber papers.
However, being new to this, I would like to ask for your comments on
whether my choices make sense. Specifically: Should I use cc/cd or
rather ca for the C atoms in the fused bonds? Is na OK for the imide N?

                    || |
                    c nd ca
                   / \ // \ / \\
               --na cc cd ca--
                  | | || |
               o==c cc cc ca--
                   \ // \ / \ //
                    nd na ca
                           | |

(Sorry for this somewhat awkward drawing... you must use a monospaced
font to view it. Dangling bonds connect to H or alkyl substituents.)

A follow-up question then is: Do I need to adjust the IMPROPER
section in the prepin file manually as well because I turned, e.g.,
pyramidal nh into planar na?

Thanks for any hints!

Best wishes

Dr. Hans Martin Senn
Lord Kelvin/Adam Smith Research Fellow
University of Glasgow
Department of Chemistry
Joseph Black Building, University Avenue, Glasgow G12 8QQ, Scotland/UK
E-mail:, phone: +44 141 330 6574

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Sun Aug 05 2007 - 06:08:10 PDT
Custom Search