Re: [AMBER] tleap adds in "missing" terminal hydrogens - how to stop?

From: David A Case <case.biomaps.rutgers.edu>
Date: Wed, 27 Jun 2012 15:51:22 -0400

On Wed, Jun 27, 2012, Brandon Sim wrote:
>
> A PDB file I am working on has several missing residues, so I added them in
> manually by copying and pasting the residue info in from another similar
> PDB file that does have those residues. However, as expected, the bond
> lengths between these residues is a bit off (~4 angstroms).
>
> I was hoping that sander would just snap these bonds back into place during
> minimization and would fix the problem. However, when I am preparing the
> inpcrd and parmtop files through tleap, it keeps thinking that because of
> the long distance between the atoms mentioned before, the two belong to
> separate molecules; hence, tleap adds in HXT to the end of the preceding
> residue.

Something else must be going on: tleap doesn't take such actions based on
bond distances, as far as I know. In particular, Amber libraries don't
have any atoms named "HXT". So we would need to know details about exactly
what you are doing--can you make a small example that exhibits the behavior
you describe?

....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jun 27 2012 - 13:00:02 PDT
Custom Search