AMBER Archive (2007)Subject: AMBER: Reg restraint value obtained during umbrella sampling
From: D.Usharani (usha_at_ipc.iisc.ernet.in)
Date: Thu May 31 2007 - 12:22:21 CDT
Dear amber users,
I am a new subscriber for amber. I am trying to do an umbrella sampling
for a 12 base pairs two DNAs separated by certain distance using amber8.
The reaction coordinate is along the two center of mass of DNA. The
distance restraint file I used is according to amber8 manual (pg 173) and
is as follows:
&rst
iat= -1,-1, nstep1=1, nstep2=6000,
iresid=1,irstyp=0,ifvari=1,ninc=0,imulti=0,ir6=0,ifntyp=0,
r1=0.00000E+00, r2=112.0000, r3=112.0000, r4=99.0000,
rk2=5.0000, rk3=5.0000,
r1a=0.00000E+00, r2a=110.0000, r3=110.0000, r4=99.0000,
rk2a=5.0000, rk3a=5.0000,
igr1 = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18,
19, 20, 21, 22, 23, 24, 0,
grnam1(1) ='P', grnam1(2)= 'P', grnam1(3)='P', grnam1(4)='P',
grnam1(5)='P', grnam1(6)='P', grnam1(7)='P', grnam1(8)='P',
grnam1(9)='P', grnam1(10)='P', grnam1(11)='P', grnam1(12)='P',
grnam1(13)='P', grnam1(14)='P', grnam1(15)='P', grnam1(16)='P',
grnam1(17)='P', grnam1(18)='P', grnam1(19)='P', grnam1(20)='P',
grnam1(21)='P', grnam1(22)='P', grnam1(23)='P', grnam1(24)='P',
igr2 =
25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,0,
grnam2(1) ='P', grnam2(2)= 'P', grnam2(3)='P', grnam2(4)='P',
grnam2(5)='P', grnam2(6)='P', grnam2(7)='P', grnam2(8)='P',
grnam2(9)='P', grnam2(10)='P', grnam2(11)='P', grnam2(12)='P',
grnam2(13)='P', grnam2(14)='P', grnam2(15)='P', grnam2(16)='P',
grnam2(17)='P', grnam2(18)='P', grnam2(19)='P', grnam2(20)='P',
grnam2(21)='P', grnam2(22)='P', grnam2(23)='P', grnam2(24)='P',
/
my sander input file is :
&cntrl
nmropt = 1,
ntx = 7, irest = 1, ntrx = 1, ntxo = 1,
ntpr = 100, ntwx = 500, ntwv = 0, ntwe = 0,
ntwprt = 0, ntwr = 500,
ntf = 2, ntb = 1, dielc = 1.0,
cut = 9.0, nsnb = 10,
scnb = 2.0, scee = 1.2,
ipol = 0,
ibelly = 0,
imin = 0,
nstlim = 100000,
nscm = 1000,
t = 0.0, dt = 0.001,
temp0 = 300.0, tempi = 300.0,
ig = 71277,
ntt = 1,
tautp = 1.0,
vlimit = 20.0,
ntp = 0, pres0 = 1.0, comp = 44.6,
taup = 0.5,
ntc = 2, tol = 0.0005,
/
&ewald
/
&wt type='DUMPFREQ', istep1=10 /
&wt type='END' /
DISANG= distcom112.RST
DUMPAVE=umbcom.112
The output file and trajectory files didn't have any errors written, but
the umbcom.112 file that has been generated during the run should have
restraint value of distance which was not printed, instead it is like this
Restraints on step 0:
0.000
Restraints on step 10:
0.000
Even with varying the force constant (rk2, rk3 ) values of 5-25kcal/mol
also it didn't print any value.
Can any one tell me where is the mistake? Is that i cant define the
residue in the center of mass restraint? If so how should be the input
distance restraint file for umbrella sampling?
Thanks in Advance
usha
Inorganic and physical chemistry
Indian Institute of Science
Bangalore -12
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|