[AMBER] peptide leaving box, image command not working

From: Omair A. Khan <oak3.psu.edu>
Date: Mon, 14 Jun 2010 13:51:08 -0400

Dear All,

I'm having a very frustrating problem. I ran a 250 ns trajectory of a
protein and peptide complex with PMEMD in Amber 10. When I generated PDB
files from the restrt (saved as .xyz) files, I saw that my peptide had
drifted extremely far away in some of the structures. It's likely that it
left the primary unit cell. However, in my Sander input file, I made sure
that iwrap was set to 1. (The exact same input and production settings used
for this trajectory have worked successfully on other ones.)

Unsure as to why iwrap didn't work in the first place (it always has for my
other simulations), I tried using every possible combination of *image*, *
center*, and *rms *to correct the mdcrd output coordinate files (saved as
.crd) so that I wouldn't see these wild fluctuations. The link below
includes my topology file, one of the output crd and xyz files, and an
example ptraj script I tried using to correct the crd file. I need to
generate corrected coordinate files because they feed directly into my
scripts for further analysis.

One thing I noticed when running the ptraj script was that although most of
the coordinates in the new crd file seemed acceptable, there were a few that
had values of more than 300. It's unlikely that the peptide or the protein
has moved such a great distance in one picosecond.

The files can be downloaded here:
http://www.omairadilkhan.com/public/amber_community.tar.gz (~200 MB)

Thanks in advance,

AMBER mailing list
Received on Mon Jun 14 2010 - 11:00:05 PDT
Custom Search