[gmx-users] Re: converting pdb to gmx successfully - how?

[Ken Rotondi]

David,

Unfortunately computers are not crazy, just very inflexible when 
executing commands.

It sounds like your hydrogen data base file (ffXXXX.hdb) is not 
complete. Are the warnings (without -missing, I really wouldn't use 
that unless you are sure that you want to use incomplete molecules) of 
the type:

WARNING: atom H62 is missing in residue ADE XXX in the pdb file

If that is the case you need to modify your .hdb file to include the 
DNA hydrogens.

Good luck,

Ken

If it doesn't specify
On Feb 2, 2006, at 5:36 PM, David Mathog wrote:

> More progress, sort of.
>
> By changing C7 -> C5M in the DTHY records I finally managed
> to get pdb2gmx -ignh -missing to run to completion.
>
> That said, it does not appear to have worked correctly.
>
> 1.  None of the DNA is protonated.
> 2.  The protein is protonated
>
> which is craziness, it should be one way or the other.
>
> 3.  There were 3 long bond warnings right after "Making bonds...".
> However I cannot figure out why they are there since the
> corresponding atom positions in the gmx output file don't correspond
> to atoms that make bonds in the rtp file.
>
> What do the position numbers in the long bond message refer to?
>
> Thanks,
>
> David Mathog
> mathog at caltech.edu
> Manager, Sequence Analysis Facility, Biology Division, Caltech
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> 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