[gmx-users] Trjconv PDB files define solvent as "ATOM"?
tsjerkw at gmail.com
Tue Mar 20 06:24:01 CET 2012
Maybe it was laziness, maybe it was the realization that the distinction
between atom and hetatm is a bit awkward that caused the developers not to
bother with it. It also doesn't relate to your problem, probably. There are
non-standard residues and covalently bound ligands that are happily part of
chains, and probably parsed correctly. Chain identifiers and ter statements
should be used to mark the polymers. Solvent should not have a chain
identifier. You can see if changing the tag helps, though:
sed '/SOL/s/^ATOM /HETATM/' in.pdb
Hope it helps,
On Mar 20, 2012 3:50 AM, "John Ladasky" <blind.watchmaker at yahoo.com> wrote:
Following up to myself:
Looking back through the archives, I have learned that the issue of solvent
atoms being exported from trjconv into PDB files as "ATOM" rather than as
"HETATM" is an issue that was raised fully six years ago.
I can hack my way around this problem. But I would still like to
understand whether, and why, I should do so.
gmx-users mailing list gmx-users at gromacs.org
Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-users