[AMBER] Specifying improper dihedrals

From: Ben Roberts <roberts.qtp.ufl.edu>
Date: Wed, 17 Jun 2009 16:24:14 +0100

Hi all,

Thanks to all those who responded regarding specifying improper
dihedrals manually. In the end, I've chosen to edit the prmtop by hand
(since I only have four impropers to add, and generating a DISANG file
seemed the more complicated of the two approaches).

So, I changed the MPHIA and NPHIA numbers at the top of the PRMTOP
file, and added the new dihedrals in the DIHEDRALS_WITHOUT_HYDROGEN
section of the PRMTOP. Specifically:

<snip>
    34623 34647 -34620 -34611 43
%FLAG EXCLUDED_ATOMS_LIST
<snip>

becomes

<snip>
    34623 34647 -34620 -34611 43 20088 20094 20091
-0 44
    21228 21237 21234 -0 44 15294 15303 15300
-3 44
    15402 15411 15408 -3 44
%FLAG EXCLUDED_ATOMS_LIST
<snip>

Which is all well and good, but when I ran "rdparm dihedrals" on this
file, it marked the new dihedrals - which is fine - but only the
second two showed up as impropers (i.e., with an I before them).

Output from rdparm:

<snip>
    45137: 1.000 3.14 2.0 :448.CG :448.CE1 :448.ND1 :1.NI
(6697,6699,6698,1)
    45138: 1.000 3.14 2.0 :474.CE1 :474.CD2 :474.NE2 :1.NI
(7077,7080,7079,1)
I 45139: 1.000 3.14 2.0 :336.CE1 :336.CD2 :336.NE2 :2.NI
(5099,5102,5101,2)
I 45140: 1.000 3.14 2.0 :338.CE1 :338.CD2 :338.NE2 :2.NI
(5135,5138,5137,2)
<snip>

Now I know that to specify an improper, one makes the fourth number
negative - as you can see I've done. But it looks to me as though,
because that fourth number is zero in this case, it's interpreted it
as a non-negative number!


Is there any way around this?

Cheers,
Ben

-- 
Benjamin P. Roberts
Postdoctoral Research Associate
Quantum Theory Project
University of Florida
2301 New Physics Building #92
PO Box 118435
Gainesville FL 32611-8435
USA
Phone: +1 352 392 6712
Cell:  +1 352 222 3677
Member of the Royal Australian Chemical Institute
and of the American Chemical Society
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 06 2009 - 09:54:46 PDT
Custom Search