[gmx-users] Issues with my PDB file
Justin A. Lemkul
jalemkul at vt.edu
Thu Apr 2 21:13:19 CEST 2009
Jacob Harvey wrote:
> Dear all,
>
> I appear to be having problems with my pdb file for gromacs. If someone
> could help me out that would be awesome. I enter the following code to
> run pdb2gmx:
>
> pdb2gmx -f test.pdb
>
> and then I select a force field. I have selected OPLS-AA/L all-atom
> force field (2001 aminoacid dihedrals). Once I select the force field I
> get the error
>
> "Fatal error:
> Residue '2' not found in residue topology database"
>
That one's pretty self-explanatory. You can't expect pdb2gmx to be magic.
Parameters have to exist in the force field .rtp file for each building block
you want to use. A residue named "2" is meaningless to pdb2gmx.
-Justin
> This is what my pdb file looks like (I've just posted the first 10-15
> lines or so). Does anyone see where I am going wrong? Thank you very
> much in advance.
>
> HETATM 1 C 2 2 -4.438 3.334 2.988
> HETATM 2 C 2 2 0.588 3.313 0.517
> HETATM 3 C 2 2 5.476 3.310 -2.216
> HETATM 4 H 2 2 -4.134 2.590 3.731
> HETATM 5 H 2 2 -5.525 3.438 3.047
> HETATM 6 H 2 2 -3.996 4.291 3.276
> HETATM 7 H 2 2 1.480 3.111 1.117
> HETATM 8 H 2 2 -0.266 3.349 1.200
> HETATM 9 H 2 2 0.702 4.307 0.076
> HETATM 10 H 2 2 6.375 3.525 -1.633
> HETATM 11 H 2 2 4.648 3.862 -1.761
> HETATM 12 H 2 2 5.628 3.709 -3.223
> HETATM 13 C 2 2 -3.996 2.922 1.581
> HETATM 14 H 2 2 -2.902 2.862 1.545
> HETATM 15 H 2 2 -4.281 3.703 0.865
> HETATM 16 C 2 2 0.394 2.242 -0.557
> HETATM 17 H 2 2 -0.485 2.486 -1.167
> .
> .
> .
> .
> .
>
> Jacob Harvey
>
> --
> --
> Jacob Harvey
>
> Graduate Student
>
> University of Massachusetts Amherst
>
> j.harv8 at gmail.com <mailto:j.harv8 at gmail.com>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> gmx-users mailing list gmx-users at 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 gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
--
========================================
Justin A. Lemkul
Graduate Research Assistant
ICTAS Doctoral Scholar
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
More information about the gromacs.org_gmx-users
mailing list