Re: AMBER: ptraj Hbond analysis

From: Thomas Cheatham <cheatham.chpc.utah.edu>
Date: Sun, 20 Nov 2005 16:13:44 -0700 (Mountain Standard Time)

> I have used the imaged trajectory file to save a
> trajectory of closest 200 waters with noimage option.
> When I viewed the trajectory in VMD the rna molecule
> is rotating with some spead in random directions. I

The RNA will rotate around (as it is expected to on a ns time scale),
ideally hopefully not in total random directions...

The reason why you see random rotation is due to the RMS command.

(1) to see if your RNA is rotating and how fast, look at the raw
trajectory in a movie (likely stripped of water and box information) with
NO rms fitting. This will show you how the RNA moves in the box (i.e. is
it slowly rotating,etc.) If you have not set NSCM to remove overall
center of mass motion of the whole box, the whole box can move over
time...

(2) Do the RMS on a particular set of residues. In your script, you have:

  trajin ade1l_9-30sets.traj
  trajout ade1l_close200.traj
  rms first
  closest 200 :1-17 oxygen noimage

If you look at the ptraj output, you will see that you are RMS fitting to
the first frame using all atoms (INCLUDING the water). As the water
exchanges very rapidly, this will lead to apparent random rotation of the
RNA.

If you want to RMS fit only to the RNA atoms,

rms first mass out rms.dat :1-17

 or

rms first :1-17

Since you are doing noimage on the closest, you might want to image the
waters prior to make sure periodically wrapped waters are included into
the closest command...

  center :1-17 mass origin
  image origin center
  rms first :1-17
  closest 200 :1-17 oxygen noimage


--tom




-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Nov 20 2005 - 23:53:01 PST
Custom Search