hi

Thanks for your reply. I tried -ignh and it doesn't help. I understand -missing is bad but I made sure the only missing atoms the program complains are the first residue in each three chains. I am using a crystal structure. However the warning messages are only for the first residue in each chain.

Zhong

On Jan 14, 2009, at 1:59 PM, Justin A. Lemkul wrote:



Zhong Zheng wrote:
Hi all
I am running Gromacs on a three-chain protein. The program complains about "atom H is missing in residue ALA 1 in the pdb file, You might need to add atom H to the hydrogen database of residue ALA in the file ff???.hdb (see the manual)". It's a warning message. I can use the option -missing to ignore it. However I would like to understand this. It's always the first residue in the chain. So I am thinking this is because they are the n-terminus. Then why no complain about c-terminus? How can I fix this?

Use -ignh to allow the relevant .hdb file to add the appropriate hydrogens. The problem is that, without -ignh, pdb2gmx expects all atoms (including H) to be present in the .pdb file. If you are using a crystal structure, this requirement will not be satisfied.

Likely pdb2gmx does not complain about the C-terminus because it is unprotonated, or otherwise the program exits after finding a missing N-terminal hydrogen.

As an aside, using -missing is generally a very bad idea, unless you have a very specific (and anticipated!) reason to use it.

-Justin

Thanks a lot.
Zhong Zheng
Ph.D candidate
Physics Dept., City College-CUNY
Marshak Science Building J-419 , 160 Convent Ave, New York, NY 10031 USA Tel: (212)650-6079 MCCE wiki: scratchpad.wikia.com/wiki/MCCE <http://scratchpad.wikia.com/wiki/MCCE >
------------------------------------------------------------------------
_______________________________________________
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 gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

--
========================================

Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
_______________________________________________
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 gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

_______________________________________________
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 gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to