[gmx-users] no residue type

Nick Johans johans.nick93 at gmail.com
Tue Jan 30 19:26:25 CET 2018


Hi, i,m a beginner. i us pdb2gmx or generating topology files.this is
aminoacid pdb file;

ATOM      1  N   HIS     1       0.102  -0.143   0.000  1.00  0.00
 N
ATOM      2  CA  HIS     1       1.558  -0.029   0.000  1.00  0.00
 C
ATOM      3  C   HIS     1       1.986   1.419   0.000  1.00  0.00
 C
ATOM      4  O   HIS     1       1.165   2.348   0.000  1.00  0.00
 O
ATOM      5  CB  HIS     1       2.047  -0.732  -1.277  1.00  0.00
 C
ATOM      6  CG  HIS     1       2.952  -1.892  -0.978  1.00  0.00
 C
ATOM      7  ND1 HIS     1       3.335  -2.320   0.291  1.00  0.00
 N
ATOM      8  CD2 HIS     1       3.517  -2.678  -1.972  1.00  0.00
 C
ATOM      9  CE1 HIS     1       4.118  -3.363  -0.046  1.00  0.00
 C
ATOM     10  NE2 HIS     1       4.282  -3.641  -1.368  1.00  0.00
 N
ATOM     11 1H   HIS     1      -0.540   0.725   0.000  1.00  0.00
 H
ATOM     12 2H   HIS     1      -0.397  -1.101   0.000  1.00  0.00
 H
ATOM     13  HA  HIS     1       1.928  -0.468   0.941  1.00  0.00
 H
ATOM     14 1HB  HIS     1       1.206  -1.113  -1.893  1.00  0.00
 H
ATOM     15 2HB  HIS     1       2.598  -0.031  -1.935  1.00  0.00
 H
ATOM     16  HD2 HIS     1       3.362  -2.530  -3.037  1.00  0.00
 H
ATOM     17  HE1 HIS     1       4.579  -3.923   0.758  1.00  0.00
 H
ATOM     18  HE2 HIS     1       4.837  -4.388  -1.810  1.00  0.00
 H
ATOM     19  OC  HIS     1       3.396   1.655   0.000  1.00  0.00
 O
ATOM     20  HC  HIS     1       3.896   0.754   0.000  1.00  0.00
 H
TER

i know i shouldn't  use AA as an independent structure so i try to add
termini by -ter then i choose amber99  and TIP3P water model but i get this
error;

In the chosen force field there is no residue type for 'HIS' as a starting
terminus

what should i do?
------------------------------
i used ambertools or generating .prmtop and .prcpd file. then
https://www.dropbox.com/s/rfjpz45we1f5h0z/amb2gmx_dihe_old.pl to convert
them to gromacs input file.

>> perl amb2gmx.pl --prmtop threonine.prmtop --crd threonine.inpcrd
--outname THRR


but thereis nothing in files as it says;
[roja at localhost cmd]$ perl amb2gmx.pl --prmtop threonine.prmtop --crd
threonine.inpcrd --outname THRR
sh: ambpdb: command not found
sh: rdparm: command not found
Parsing PDB file THRR.pdb...
found:
    0 protein or solute residues
    0 sodium ions
    0 chloride ions
    0 magnesium ions
    0 waters

    0 TOTAL residues read
Reading coordinates from CRD file...Reading prmtop file...
0 atom type indices read.

found:
        0 atoms
        0 bonds
        0 angles
        0 torsions
BAD BOX OR BOX NOT DETECTED; USING DUMMY BOX SIZE IN GRO FILEBAD BOX OR BOX
NOT DETECTED; USING DUMMY BOX SIZE IN GRO FILE
Qtot = 0
Qmax = 0 at atom -1

Constructing list of unique nonbonded atomtypes...
0 unique atomtypes found.
Writing atom types...
Done.
------------------------
is there sth i didn't consider?

Best regards
-Nick


More information about the gromacs.org_gmx-users mailing list