Hi,
I am trying to calculate RMSd with amber trajectory "crd" file using
"carnal". (It is peptide with about 35 residues).
I need to calculate the RMSD for every crd time frame(about 1000 frames)
based on backbone(or CA atom).
The problem is that I need to compare with reference which is given
in "pdb" format. I know I can set reference for RMSd in "carnal" but
it seems it doesn't recognize pbd format.
How can I calculate RMSD with this ?
Here is my script with running results.
Thanks.
PS. I converted crd file into pdb file using VMD.
I know 'VMD' can calculate RMSD also. Is there any method I can do in
VMD ?
protege[3] 14> more rmsd.in
FILES_IN
PARM p1 vii.parm;
STREAM s1
md_200_rmsd.pdb;
STATIC ref_set viiu.pdb;
FILES_OUT
TABLE tbl sm.rms;
DECLARE
GROUP grp1 ( ATOM NAME CA );
RMS fit1 FIT grp1 s1 ref_set;
OUTPUT
TABLE tbl fit1 ;
END
protege[3] 15> carnal<rmsd.in
- CARNAL -
AMBER 6.0
COORDINATE ANALYSIS
input stdin
> FILES_IN
> PARM p1 vii.parm;
Reading parm file (vii.parm)
parm: opening vii.parm
> STREAM s1
> md_200_rmsd.pdb;
Using default parm (vii.parm) for STREAM s1
stream: opening md_200_rmsd.pdb
> STATIC ref_set viiu.pdb;
Using default parm vii.parm for STAT ref_set
stat: opening viiu.pdb
Error: viiu.pdb: full pdb format not supported
Only ATOM records in amber order allowed
protege[3] 16>
Received on Fri Oct 04 2002 - 00:02:37 PDT
