[AMBER] NMR restraint ?

From: xuemeiwang1103 <xuemeiwang1103.163.com>
Date: Wed, 5 May 2010 15:28:44 +0800 (CST)

Dear all:
      I want to constraint the distance between four atoms(GLU355 OE1-WAT588 H1,GLU355 OE2-WAT637 H1),so I use the NMR restraint in my mdin file , After minimal the system,I began to heat it to 300K over 6 stages.The strange thing is that the NMR restraint seems only work at first two stages ,when the Tem reach to 300K,the distance is still far.Could someone tell me the reason? Thank you !


My input file:

Stage 1 heating 50 to 100K

&cntrl

 imin=0,irest=1,ntx=5,

 nstlim=5000,dt=0.001,

 ntc=2,ntf=2,

 ntt=3,gamma_ln=1.0,

 tempi=50.0,temp0=100.0,

 ntb=1,cut=8,

 nmropt=1,

 ntpr=250,ntwx=500,ntwr=500,

/

&wt type='END' ,

/

LISTIN=POUT

LISTOUT=RST.out

DISANG=dis.RST

The dis.RST

#DISTANTENT.FILE

 &rst

  ixpk= 0, nxpk= 0, iat= 9299, 5551, r1= 0.84, r2= 1.34, r3= 2.34, r4= 2.84, rk2= 25.0, rk3= 25.0, ir6=1,ialtd= 0,

/

 &rst

  iat= 69425, 5550, r1= 0.84,r2= 1.34, r3= 2.34, r4=2.84,rk2= 25.0,rk3= 25.0,

/

Outputfile:

3. ATOMIC COORDINATES AND VELOCITIES

-------------------------------------------------------------------------------

 begin time read from input coords = 5.000 ps

          Begin reading energy term weight changes/NMR restraints

 WEIGHT CHANGES:

                         ** No weight changes given **

 RESTRAINTS:

 Requested file redirections:

  LISTIN = POUT

  LISTOUT = RST.out

  DISANG = dis.RST

 Restraints will be read from file: dis.RST

Here are comments from the DISANG input file:

#DISTANTENT.FILE

******

 H1 ( 9299)-OE2 ( 5551) NSTEP1= 0 NSTEP2= 0

R1 = 0.840 R2 = 1.340 R3 = 2.340 R4 = 2.840 RK2 = 25.000 RK3 = 25.000

 Rcurr: 2.272 Rcurr-(R2+R3)/2: 0.432 MIN(Rcurr-R2,Rcurr-R3): 0.000

******

 H1 (69425)-OE1 ( 5550) NSTEP1= 0 NSTEP2= 0

R1 = 0.840 R2 = 1.340 R3 = 2.340 R4 = 2.840 RK2 = 25.000 RK3 = 25.000

 Rcurr: 2.123 Rcurr-(R2+R3)/2: 0.283 MIN(Rcurr-R2,Rcurr-R3): 0.000

                       Number of restraints read = 2

                  Done reading weight changes/NMR restraints

 Input file:

Stage 1 heating 150 to 200K

&cntrl

 imin=0,irest=1,ntx=5,

 nstlim=5000,dt=0.001,

 ntc=2,ntf=2,

 ntt=3,gamma_ln=1.0,

 tempi=150.0,temp0=200.0,

 ntb=1,cut=8,

 nmropt=1,

 ntpr=250,ntwx=500,ntwr=500,

/

&wt type='END' ,

/

lISTIN=POUT

LISTOUT=RST.out

DISANG=dis.RST

Outputfile:

3. ATOMIC COORDINATES AND VELOCITIES

-------------------------------------------------------------------------------

 begin time read from input coords = 10.000 ps

           Begin reading energy term weight changes/NMR restraints

 WEIGHT CHANGES:

                         ** No weight changes given **

 RESTRAINTS:

  No valid redirection requests found

                          ** No restraint defined **

                 Done reading weight changes/NMR restraints






_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed May 05 2010 - 00:30:03 PDT
Custom Search