[gmx-users] ACE parameter error in ffG53a5-6?
Mauricio Pablo Sica
msica at unq.edu.ar
Wed Jul 25 15:52:26 CEST 2007
Dear all
sorry for being here the 3rd time with this same (maybe trivial)
question.
I have found what appears to be an error in the files ffG53a(5 and 6).rtp
related to the ACE group topology. Maybe it is a trivial error, or a
change due to something I am missing.
I was constructing the topology for N-acetylgalactosamine ffG53a6
parameters. In general, my strategy was to acetylate a 2-amino galactose.
According to data in ffG53a6bon.itp, “my” acetyl group bonded to the
sugar amino should have:
[ bonds ]
CA C gb_27
C O gb_5
C N gb_10 ; parameter for the peptide (amide) bond.
I choose gb_10 according to ffG53a6bon.itp, where:
#define gb_10 0.1330 1.1800e+07
; C, CR1 - N, NR, CR1, C (peptide, 5-ring) 1000
The problem arose when I compared what I did to the ACE group in the
ffG53a6.rtp. There I found that:
[ bonds ]
C CA gb_27
C O gb_5
C +N gb_19; <<<< (Here is the discrepancy)
According to ffG53a6bon.itp, gb_19 is a CHn-O bond, as you can see:
#define gb_19 0.1430 9.2100e+06
; CHn - OM 900
Could it be an error or the choice is due to something I am missing?
Additionally, in the force fields 43a1, 43a2, 43b1 and 45a3, the
parameter for the said C-N bond is gb_9, which in all cases is defined as:
define gb_9 0.1330 1.1800e+07
; C, CR1 - N, NR, CR1, C (peptide, 5-ring) 1000
Also, in ffG53a6.rtp, "C +N" bonds are gb_10 for all residues, except for
ACE.
The question, again, is: is gb_19 correct? If yes, please tell me why.
Sorry for the length of the mail and thanks for your attention.
________________________________________________________________
Dr. Mauricio P. Sica LEPP
Phone: (54011)4365-7100 ext. 169 Laboratorio de Expresión
R. Saenz Peña 354 y Plegado Proteico
(B1876BXF) Bernal, Buenos Aires Universidad de Quilmes
Argentina
________________________________________________________________
More information about the gromacs.org_gmx-users
mailing list