On Wed, May 27, 2009, Tomasio, Susana wrote:
>
> I noticed that the N-cap residue Acetyl (ACE) in the amoeba force field
> is in two files (amoeba_amino.off and amoeba_aminont.off), but has
> different atom labels (2nd column): One is labelled with numbers and the
> other one is labelled with letters.
> In the amoeba_amino.off:
> "HH31" "181" 0 1 131072 1 1 0.076010
> "CH3" "180" 0 1 131072 2 6 -0.190264
> "HH32" "181" 0 1 131072 3 1 0.076011
> "HH33" "181" 0 1 131072 4 1 0.076010
> "C" "182" 0 1 131072 5 6 0.512403
> "O" "183" 0 1 131072 6 8 -0.550170
>
> And in the amoeba_aminont.off:
> "HH31" "HC" 0 1 131072 1 1 0.112300
> "CH3" "CT" 0 1 131072 2 6 -0.366200
> "HH32" "HC" 0 1 131072 3 1 0.112300
> "HH33" "HC" 0 1 131072 4 1 0.112300
> "C" "C" 0 1 131072 5 6 0.597200
> "O" "O" 0 1 131072 6 8 -0.567900
>
> When I tried to load the ACE residue it didn't work, but by changing the
> 2nd column in the amoeba_aminont.off to include the labels that are in
> the the amoeba_amino.off file, it did work.
Thanks for the note. This certainly looks like a bug. I'm cc-ing this to Wei
Zhang, who is our (only?) amoeba expert, but to me, what you did looks
correct: the other entries in amoeba_aminont.off us numbers for atom types.
Did you ever solve the other problems you were having with amoeba simulations?
I've been leaving that up to Wei as well, but we can probably recruit others
to look at this if needed. If you are having problems, are there any
differences between pmemd.amba and sander? That might help pin down the
problem.
...thanks!...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu May 28 2009 - 01:08:40 PDT
