[gmx-users] a pdb2gmx problem for NON aminoacids

Markus O Kaukonen mokaukon at cc.helsinki.fi
Thu Nov 29 09:02:22 CET 2001


Hi,
I generated by hand a rather poor starting guess for CH3CH2OH
in xyz format:
---------------
9
ch3ch2oh to be done
C       0.      0.      0.
H       0.      0.      -1.0
H       -.7.    -.7     0.0
H       .7      .7      0.0
C       0.      0.      1.5
H       0.      1.      1.5
H       0.      -1.     1.5
O       0.      0.      3.0
H       0.      0.      4.0
--------------------------

then translated it to pdb format:
----------
ATOM      1  C   UNK     1       0.000   0.000   0.000  1.00  0.00
ATOM      2  H   UNK     1       0.000   0.000  -1.000  1.00  0.00
ATOM      3  H   UNK     1      -0.700  -0.700   0.000  1.00  0.00
ATOM      4  H   UNK     1       0.700   0.700   0.000  1.00  0.00
ATOM      5  C   UNK     1       0.000   0.000   1.500  1.00  0.00
ATOM      6  H   UNK     1       0.000   1.000   1.500  1.00  0.00
ATOM      7  H   UNK     1       0.000  -1.000   1.500  1.00  0.00
ATOM      8  O   UNK     1       0.000   0.000   3.000  1.00  0.00
ATOM      9  H   UNK     1       0.000   0.000   4.000  1.00  0.00
-------------------

then tried to change it to gromacs format end get the topology file:
pdb2gmx -f c-c-o.pdb
but got
Fatal error: Residue 'UNK' not found in residue topology database
after giving 0 for force field.

Must all the atoms in the system really belong to some aminoacid?

I also tried residue id 'SOL' but got error:
Fatal error: Atom C in residue HOH 1 not found in database
             while sorting atoms

What is the simple workaround?
I just want to optimise the structure of this molecule.

Cheers, Markus
-- 
++++++ OFFICE ++++++++++++++++++++++++++++++++ HOME ++++++++++++
Markus.Kaukonen at helsinki.fi   \          Verajakallionkatu 6 a 5
Helsinki Bioenergetics Group   \_______  02600 Espoo, Finland
Institute of Biotechnology	       \ tel:+358-(0)9-5127122/h
Biocenter 2, P.O.BOX 56 (Viikinkaari 5) \____
FIN-00014 University of Helsinki              \______ H+ ______
tel: +358-(0)9-19159749/o, fax: +358-(0)9-19158001             \
--- Rikos ei kannata, eika maatalous Suomessa. (Paimio 1998) ---





More information about the gromacs.org_gmx-users mailing list