[AMBER] distance calculation by reading in the restraint file

From: xiaoqin huang <xqhuang1018.msn.com>
Date: Fri, 6 Feb 2009 12:52:56 -0500

Hi, all
anyone can tell me what is wrong here.
 
I used the DISANG=dus1.f in the input file to read in the restraint file dus1.f
 
and the content of the dus1.f is:
-------------------------------------
 &rst IRESID=1, iat(1)=197,ATNAM(1)=CB,ATNAM(1)=OG, iat(2)=532, ATNAM(2)=O31,ATNAM(2)=C32,ATNAM(2)=O33,ATNAM(2)=C34,ATNAM(2)=C35,ATNAM(2)=C37,ATNAM(2)=C39,ATNAM(2)=C41,ATNAM(2)=C43, r1=0.00,r2=4.650,r3=4.650,r4=20.00,rk2=80,rk3=80, &end------------------------------------------
 
it is supposed to restrain the distance between the group1 atoms and the group 2 atoms, and according to the amber8 manual, such distance is calcuated as the DISTANCE BETWEEN THE AVERAGEED COORDINATES FOR GROUP 1 ATOMS AND THE AVERAGED COORDINATES OF GROUP 2 ATOM.
 
SO, IN THE MD OUTPUT FILE, I see:
 
-------------------------------------------------
 Here are comments from the DISANG input file: Number of restraints read = 1
                  Done reading weight changes/NMR restraints--------------------------------------------------------
 
and the fist number of the dumped file is 4.450.
 
THEN, I calculated the distance according to the definition in the manual, i.e. coordinates-averaged, using the input coordinates, it is 3.927.
 
it is so different, where is the error?
 
and how to get more reaing-in information about the restratints from the MD output file?
it is better to output a piece of information like " matches ??? atoms" in the output file.
 
 
thanks
 
xiaoqing
 
 
 
 
 
 
 
 
 
 
 
 
 
 

_________________________________________________________________
Windows Live™: Keep your life in sync.
http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_022009_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Feb 08 2009 - 01:11:39 PST
Custom Search