On 3/1/13 3:49 PM, Villarealed wrote:
Dear Gromacs-Users,
I am using Gromos96 43a1 FF extended to include phosphorylated residues.
I have a peptide-capped with 17 aminoacids.
When I tried to add ions to the virtual box
I obtained these error
ERROR 1 [file topol.top, line 566]:
No default G96Angle types
ERROR 2 [file topol.top, line 801]:
No default Proper Dih. types
The second error was resolved by modified the archive residuetypes.dat.
I only add this line SEP Protein
To be honest I don't know how to resolve the ERROR 1.
If someone could help me I will appreciate.
Eduardo
grompp -f ions.mdp -c newbox_solv.gro -p topol.top -o ions.tpr
ERROR 1 [file topol.top, line 566]:
No default G96Angle types
topol.top
[ angles ]
; aiajak functc0c1c2
c3
1 2 3 2
newbox_solc.gro
1ACE CA1 2.759 1.040 1.830
1ACE C2 2.881 1.130 1.839
1ACE O3 2.988 1.083 1.878
2SER N4 2.872 1.258 1.805
aminoacids.rtp
[ ACE ]
[ atoms ]
CA CH3 0.000 0
C C 0.380 1
O O -0.380 1
[ bonds ]
CCA gb_26
C O gb_4
C+N gb_9
[ angles ]
CA C+Nga_18
O C+Nga_32
[ impropers ]
CCA+N Ogi_1
It looks like you are using an old version of Gromacs. This was a bug that was
fixed some time ago. There is a missing angle here:
CA C Oga_30
ERROR 2 [file topol.top, line 801]:
No default Proper Dih. types
topol.top
[ dihedrals ]
; aiajakal functc0c1c2
c3c4c5
72717374 1
H N CACB
There should not be a dihedral assigned here. It is not specified in the .rtp,
nor do the Gromos force fields assign dihedrals for H-N-X-X. The -N-CA- torsion
is defined by heavy atoms only. Are you sure the numbering is correct? The
numbers listed in topol.top correspond to the [moleculetype] numbering, which is
not necessarily the same numbering in the coordinate file, if (for example) the
residue/atom numbering does not start at 1 in the coordinate file.
-Justin
newbox_sol.gro
10SEP N 71 2.348 2.562 2.883
10SEP H 72 2.446 2.556 2.863
10SEP CA 73 2.260 2.594 2.768
10SEP CB 74 2.186 2.467 2.720
10SEP OG 75 2.110 2.424 2.827
10SEP P 76 2.001 2.314 2.797
10SEPO1P 77 2.066 2.200 2.732
10SEPO2P 78 1.949 2.275 2.928
10SEPO3P 79 1.895 2.372 2.710
10SEP C 80 2.164 2.714 2.789
10SEP O 81 2.051 2.714 2.741
aminoacids.rtp
[ SEP ]tar -zcvf prog-1-jan-2005.tar.gz /home/jerry/prog
[ atoms ]
N N-0.28000 0
H H 0.28000 0
CA CH1 0.0 1
CB CH2 0.0 2
OGOA-0.36000 2
P P 0.94000 2
O1POP-0.86000 2
O2POP-0.86000 2
O3POP-0.86000 2
C C 0.380 3
O O -0.380 3
[ bonds ]
N Hgb_2
NCAgb_20
CA Cgb_26
C Ogb_4
C+Ngb_9
CACBgb_26
CBOGgb_17
OG Pgb_27
P O1Pgb_23
P O2Pgb_23
P O3Pgb_23
[ angles ]
; aiajak gromos type
-C N H ga_31
H NCA ga_17
-C NCA ga_30
NCA C ga_12
CA C+N ga_18
CA C O ga_29
O C+N ga_32
NCACB ga_12
CCACB ga_12
CACBOG ga_12
CBOG P ga_25
OG P O1P ga_13
OG P O2P ga_13
OG P O3P ga_13
O1P P O2P ga_28
O1P P O3P ga_28
O2P P O3P ga_28
[ impropers ]
; aiajakal gromos type
N-CCA H gi_1
CCA+N O gi_1
CA N CCB gi_2
[ dihedrals ]
; aiajakal gromos type
-CA-C NCA gd_4
-C NCA C gd_19
NCA C+N gd_20
NCACBOG gd_17
CACBOG P gd_14
CBOG P O1P gd_11
-
Eduardo Villarreal Ramírez
Postdoctoral Research Fellow
Mineralized Tissue Laboratory,
Hospital for Special Surgery.
--
View this message in context:
http://gromacs.5086.n6.nabble.com/ERROR-No-default-G96Angle-types-tp5006044.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
--
Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin