AMBER Archive (2007)

Subject: Re: AMBER: Problem using Ptraj

From: Thomas Cheatham III (tec3_at_utah.edu)
Date: Mon Apr 23 2007 - 12:02:17 CDT


> I am having troubles to prepare the snapshots for mm_pbsa
> analysis. I have the trajectory file and the corresponding prmtop file, wich I
> can visualize in VMD without problems. The point is that after striping the
> waters and ions with ptraj, and obtaining the new desolvated and deionized
> prmtop files, when I try to check the new trajectory, I find that the ligand
> atoms are all missplaced, with corrupted binding distances, angles and lenghts.

Usually this is a problem with the box information, either still being
present when not expected or not being present when expected. The first
frame/snapshot will look good and all subsequent will be offset.

ptraj by default (and stupidly) automatically outputs box information;
MM-PBSA can read box information or not, depending on the "BOX" keyword.

If you are stripping, generally you want "nobox" appended to the trajout
line,

strip :WAT
trajout mdcrd.strip nobox
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu