Re: AMBER: RMSD calculations in ptraj and Carnal

From: Joseph Nachman <nachman.hera.med.utoronto.ca>
Date: Mon, 21 Jun 2004 11:07:57 -0400 (EDT)

I tried mass-weighting, and rhe results were identical to those woithout
mass weighting, which is as expected, since I'm calculating RMSDs only for
CAs.

What puzzles me is that the is that the RMSDs calculated by ptraj seem to
have no relation whatsoever to those calculated by hand (and thesame holds
true for Carnal):

RMSD(manual):

    1 2.05
    2 3.20
    3 3.63
    4 5.22
    5 4.04

RMSD(ptraj):

    1.00 2.04581
    2.00 2.15204
    3.00 2.50675
    4.00 3.22698
    5.00 0.81087

RMSD(carnal)


L0 8.693087e-03
L0 3.986790e-03
L0 2.426629e-03
L0 1.264287e-02
L0 1.153996e-02

The inputs I used are:

ptraj:

trajin 1.xyz
trajin 2.xyz
trajin 3.xyz
trajin 4.xyz
trajin 5.xyz
reference minhm.xyz
rms reference out test.rms :106-110.CA

Carnal:

FILES_IN
  PARM p1 1spshm.top;
  STREAM s1 1.xyz 2.xyz 3.xyz 4.xyz 5.xyz;
  STATIC ref_set minhm.xyz;
FILES_OUT
  TABLE t1 table.rms;
DECLARE
  GROUP grpid ((ATOM TYPE CA) & (RES 106-110));
  RMS fit1 grpid s1 ref_set;
OUTPUT
  TABLE t1 fit1;
END

All the structures are already fitted, since most of the structure was
kept fixed (ibelly option in sander-classic, Amber6), and only a few
residues allowed to move.

Any help and ideas welcome!

Thanks in advance,

Joseph

 On Tue, 15 Jun 2004, Thomas E. Cheatham, III wrote:

>
> By default, ptraj RMS fits *without* mass weighting. To do mass
> weighting, simply add the "mass" keyword. Back a long time ago when testing, carnal
> and ptraj give equivalent RMSd results (~1995) and the routines have not
> changed significantly since then.
>
> rms first mass out rms.dat .CA
>
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
>

-- 
--------------------------------------------------------------------
Joseph Nachman				Department of Biochemistry
nachman.hera.med.utoronto.ca		University of Toronto
					Medical Sciences Building
tel: +1 416 978-5510			Toronto, Ontario M5S 1A8
fax: +1 416 978-8548			Canada
--------------------------------------------------------------------
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Mon Jun 21 2004 - 16:53:01 PDT
Custom Search