On 8/12/2011 9:40 PM, Matthias Ernst wrote:
Hi,
I try to simulate DNA or better parts of a long DNA double-helix.
Therefore, the system I use is not "correctly" terminated at the ends
from chemical point of view with OH-groups but rather cut out the
longer helix. Nevertheless, all residues are named correctly, all
atoms needed in the residues are present and corresponding parameters
exist in the Amber-Forcefield which I intent to use.
OK, but the valences still have to be satisifed. It's up to you to
decide how.
When I run pdb2gmx, it always tries to recognize start end ending
termini:
Identified residue DG1 as a starting terminus.
Identified residue DC34 as a ending terminus.
and then later on, it fails because the P-Atoms are not found in the
rtp as it searches e.g. in DG5 instead of DG:
Fatal error:
Atom P in residue DG 1 was not found in rtp entry DG5 with 31 atoms
while sorting atoms.
How can I avoid this error? I tried -noter option but this does not help.
Is there a way to just treat the residues as the occur in the pdb-file
without automatically detection of starting and ending terminus?
pdb2gmx can only put pieces together that are defined in the forcefield
.rtp and .tdb. It knows that termini have to be different from mid-chain
residues. Either it needs to find a named capping residue in the
coordinate file that matches the .rtp, or the -ter mechanism needs to be
used to generate a terminus from the .tdb on a normal residue in the
.rtp identified from the coordinate file.
Mark
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
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-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists