thanks for the hints but unfortunately i still can not get the molecule running. i tried with both prodrg server and x2top but with limited success:

in the first case i generate the gro.pdb and .top file via the webinterface [1] of prodrg by pasting my coordinate file into the textbox. then i run the following commands:

editconf -f sedoheptulose.gro.pdb -o sedoheptulose.cube.pdb -bt cubic -d 2.0
genbox -cp sedoheptulose.cube.pdb -cs spc216.gro -o sedoheptulose.water.pdb -p sedoheptulose.top -nice 0 grompp -np 1 -f em.mdp -c sedoheptulose.water.pdb -p sedoheptulose.top -o sedoheptulose.em.tpr

but in the grompp commands the following error occurs:

-------------------------------------------------------
Program grompp, VERSION 3.3.1
Source code file: topio.c, line: 388

Fatal error:
Invalid order for directive moleculetype, file ""sedoheptulose.top"", line 17
-------------------------------------------------------

but in line 17 only "[ moleculetype ]" is written?



in the second case I do not really understand the error message:

x2top -f sedoheptulose.pdb -o sedoheptulose.top

leads to:

-------------------------------------------------------
Program x2top, VERSION 3.3.1
Source code file: futil.c, line: 537

Fatal error:
Library file ffG43a1.n2t not found in current dir nor in default directories.
(You can set the directories to search with the GMXLIB path variable)
-------------------------------------------------------


also in the faq i could not find a real guide how to simulate such a molecule ...



Bernhard




[1] http://davapc1.bioch.dundee.ac.uk/prodrg/index.html



[EMAIL PROTECTED] wrote:

Message: 1
Date: Wed, 19 Nov 2008 13:25:07 +0000
From: "Kukol, Andreas" <[EMAIL PROTECTED]>
Subject: RE: [gmx-users] sedoheptulose
To: Discussion list for GROMACS users <gmx-users@gromacs.org>
Message-ID:
        <[EMAIL PROTECTED]>
Content-Type: text/plain; charset="iso-8859-1"

Yes, it is possible by deveoping your own topology. Pdb2gmx is not suitable for 
that purpose.

It would be helpful to check out the documentation resources and manual, e.g.

http://wiki.gromacs.org/index.php/FAQs

x2top and the prodgr web-server might be some good starting points.

Andreas

From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Bernhard Knapp
Sent: 19 November 2008 12:54
To: gmx-users@gromacs.org
Subject: [gmx-users] sedoheptulose

Dear users

I would like to perform a MD simulation of Sedoheptulose (C7H14O7) via Gromacs. 
But already during pdb2gmx it comes up with:
"Residue 'UNK' not found in residue topology database" which is obviously clear 
since sedoheptulose is not an amino acid. Is it somehow possible to simulate this 
molecule?

the sedoheptulose data to start from is:
ATOM      1  C   UNK     0      -3.442   2.742   0.000  0.00  0.00           C+0
ATOM      2  C   UNK     0      -2.108   3.512   0.000  0.00  0.00           C+0
ATOM      3  C   UNK     0      -0.774   2.742   0.000  0.00  0.00           C+0
ATOM      4  C   UNK     0       0.559   3.512   0.000  0.00  0.00           C+0
ATOM      5  C   UNK     0       1.893   2.742   0.000  0.00  0.00           C+0
ATOM      6  C   UNK     0       3.227   3.512   0.000  0.00  0.00           C+0
ATOM      7  O   UNK     0      -0.774   1.202   0.000  0.00  0.00           O+0
ATOM      8  O   UNK     0       0.559   5.052   0.000  0.00  0.00           O+0
ATOM      9  O   UNK     0       1.893   1.202   0.000  0.00  0.00           O+0
ATOM     10  O   UNK     0      -4.775   3.512   0.000  0.00  0.00           O+0
ATOM     11  O   UNK     0      -2.108   5.052   0.000  0.00  0.00           O+0
ATOM     12  H   UNK     0      -6.109   2.742   0.000  0.00  0.00           H+0
ATOM     13  H   UNK     0      -0.774   5.822   0.000  0.00  0.00           H+0
ATOM     14  H   UNK     0       1.893   5.822   0.000  0.00  0.00           H+0
ATOM     15  H   UNK     0       0.559   0.432   0.000  0.00  0.00           H+0
ATOM     16  H   UNK     0       3.227   0.432   0.000  0.00  0.00           H+0
ATOM     17  C   UNK     0       4.560   2.742   0.000  0.00  0.00           C+0
ATOM     18  O   UNK     0       5.894   3.512   0.000  0.00  0.00           O+0
ATOM     19  H   UNK     0       7.228   2.742   0.000  0.00  0.00           H+0
ATOM     20  O   UNK     0       3.227   5.052   0.000  0.00  0.00           O+0

cheers
bernhard


_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to