Re: [AMBER] RMSD with cpptraj

From: Rui, Huan <hrui.amgen.com>
Date: Fri, 9 Feb 2018 19:50:37 +0000

A correction. The CA and backbone atoms do have the same order but not the sidechain heavy atoms.

-----Original Message-----
From: Rui, Huan
Sent: Friday, February 09, 2018 2:46 PM
To: AMBER Mailing List <amber.ambermd.org>
Subject: RE: [AMBER] RMSD with cpptraj

Thanks for the reply, Carlos. The heavy atoms do have the same order but the reference structure does not have hydrogens. The reference structure is also 2 residues shorter even though the missing residues are not included in the RMS calculation command. Could any of the two be the issue?

-----Original Message-----
From: Carlos Simmerling [mailto:carlos.simmerling.gmail.com]
Sent: Friday, February 09, 2018 2:41 PM
To: AMBER Mailing List <amber.ambermd.org>
Subject: Re: [AMBER] RMSD with cpptraj

are the atoms in the same order in the pdb file as in the trajectory file?

On Fri, Feb 9, 2018 at 2:38 PM, Rui, Huan <hrui.amgen.com> wrote:

> Hi all,
>
>
>
> I'm trying to use cpptraj to compute the RMSD between my trajectory
> and a crystal structure. The script works fine, but the RMSD values
> seems incorrect. My reference is this clean_update.pdb. it has the
> same residues as the protein in my trajectory. The first rms command
> produced RMSD at around 6 Ang for all the snapshots, but when I look
> at the final structure overlay produced by my "trajout" command, the
> two structures do not look like they are 6 Ang apart. Also after the
> alignment from the first rms command, I wrote down the distance
> between two CA atoms from resid 100, the output rmsd_test2.plo show
> that they are 4 Ang apart, though in the superposed pdb files they are
> only 0.5 Ang away. Here's how my cpptraj script look like:
>
>
>
> trajin step4_1.nc 1 500 1
>
> parm clean_update.pdb
>
> reference clean_update.pdb parm clean_update.pdb
>
> autoimage :1-170
>
> rms reference :3-165.CA :3-165.CA out rmsd_test.plo
>
> rms reference :100.CA :100.CA nofit out rmsd_test2.plo
>
> trajout testrmsd.pdb onlyframes 500
>
> run
>
> quit
>
>
>
> Could anybody give me some insights on why this is?
>
>
>
> Thank you very much in advance.
>
>
>
> Thanks,
>
> Huan
>
>
>
>
>
>
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.ambermd.org_
> mailman_listinfo_amber&d=DwICAg&c=Sexio4usKrYWFsrnxgjbcQ&r=btLuLrQZVp2
> xq6t8VXeaNQ&m=BqOrzECEElfhW80JOQ8f9Vjn9nBJnd8npQPEq6c2O2k&s=gPsYcM-WXb
> pUD7MlCeULDHulY3XugJCOjl1tSdT0-pM&e=
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.ambermd.org_mailman_listinfo_amber&d=DwICAg&c=Sexio4usKrYWFsrnxgjbcQ&r=btLuLrQZVp2xq6t8VXeaNQ&m=BqOrzECEElfhW80JOQ8f9Vjn9nBJnd8npQPEq6c2O2k&s=gPsYcM-WXbpUD7MlCeULDHulY3XugJCOjl1tSdT0-pM&e=

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Feb 09 2018 - 12:00:03 PST
Custom Search