AMBER Archive (2004)Subject: Re[4]: AMBER: rmsd over residues using ptraj
From: Qiang Lu (qiangl_at_uci.edu)
Date: Wed May 19 2004 - 14:19:33 CDT
Hi Bill,
Thanks for your help. I finished it. The input of carnal is as the following:
FILES_IN PARM p1 1tsr_b_ww_nc_nd_lb_sa_5.parm7;
STREAM s1 1tsr_b_ww_nc_nd_lb_sa_5.10ns.crd;
STATIC ref_set 1tsr_b_ww_nc_nd_lb_sa_5.00.rst;
FILES_OUT
TABLE tb1 sm1.rms;
TABLE tb2 sm2.rms;
TABLE tb3 sm3.rms;
TABLE tb4 sm4.rms;
TABLE tb5 sm5.rms;
DECLARE
GROUP grp ((ATOM NAME CA C O N)&(RES 1-194)); RMS fit FIT grp s1 ref_set;
GROUP g1 ((RES 1)&(ATOM NAME CA C O N)); RMS fit1 g1 fit;
GROUP g2 ((RES 2)&(ATOM NAME CA C O N)); RMS fit2 g2 fit;
GROUP g3 ((RES 3)&(ATOM NAME CA C O N)); RMS fit3 g3 fit;
GROUP g4 ((RES 4)&(ATOM NAME CA C O N)); RMS fit4 g4 fit;
GROUP g5 ((RES 5)&(ATOM NAME CA C O N)); RMS fit5 g5 fit;
OUTPUT
TABLE tb1 fit1;
TABLE tb2 fit2;
TABLE tb3 fit3;
TABLE tb4 fit4;
TABLE tb5 fit5;
END
--
Best regards,
Qiang mailto:qiangl_at_uci.edu
=================Original message text===============
> That will align the proteins only based on the marked residue
> but not the whole protein. That will cause error.
The carnal version allows simply "rmsid%residues" to get the
per-residue rms in a table, where the protein can be positioned
to fit to any subset of atoms in the reference frame including
(default) all atoms.
Bill
-----------------------------------------------------------------------
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|