|
|||||||||||||||||||||||||||||||||
AMBER Archive (2008)Subject: AMBER: about NMR refinement in explicit solvent
From: Rogelio Hernández (rahl_emc_at_yahoo.com.mx)
Dear Amber Users
I am using one of the protocols described by Xia et. al.
The protocol is:
&cntrl
ntwe = 5, ntwx = 5, ntpr = 5,
scnb = 2.0, scee = 1.2, nsnb = 25, dielc = 1, cut = 15.0,
ntb = 0, igb = 0,
maxcyc = 1000, ntmin = 1, ncyc = 100, dx0 = 0.01, drms = 0.01,
ntr = 0,
imin = 1, nmropt = 0,
/
then 2 cycles of 30 ps MD simulated annealing in vacuum, from
&cntrl
nstlim=30000, pencut=-0.001,
ntpr=200, ntt=1, ntwx=200,
cut=15.0, ntb=0, vlimit=10,
/
&ewald
eedmeth=5,
/
&wt type='TEMP0',
value2=0.0, /
value2=1.0, /
&wt type='REST',
value2=1.0, /
&wt type='END' /
LISTOUT=POUT
DISANG=RST
and 2000 steps energy minimization with NMR constraints in
&cntrl
ntwe = 5, ntwx = 5, ntpr = 5,
scnb = 2.0, scee = 1.2, nsnb = 25,
ntb = 0,
maxcyc =2000, ntmin = 1, ncyc = 500,
ntr = 0,
imin = 1, nmropt = 1,
/
&wt type='REST',
value2=1.0, /
&wt type='END' /
LISTOUT=POUT
DISANG=RST
REFINEMET IN EXPLICIT SOLVENT
&cntrl
nstlim = 20000, dt = 0.002
ntwe = 500, ntwx = 500, ntpr = 50,
scnb = 2.0, scee = 1.2, nsnb = 25, dielc = 1, cut = 10.0,
ntb = 1,
ntr = 1,
ntc = 2, ntf = 2,
ntx = 1, irest = 0,
imin = 0, nmropt = 0,
ntt=1, tempi=0.0, temp0=50.0, tautp=1.0,
/
Keep protein fixed with weak restraints
5.0
RES 1 36
END
END
40 ps of MD with constant volume, gradually heating up from 50
nstlim = 20000, dt = 0.002
ntwe = 500, ntwx = 500, ntpr = 50,
scnb = 2.0, scee = 1.2, nsnb = 25,
ntb = 1,
ntr = 1,
ntc = 2, ntf = 2,
ntx = 1, irest = 0,
imin = 0, nmropt = 0,
ntt=1, tempi=50.0, temp0=100.0,
/
Keep protein fixed with weak
5.0
RES 1 36
END
END
40 ps of MD with constant volume, gradually heating up from 100
nstlim = 20000, dt = 0.002
ntwe = 500, ntwx = 500, ntpr = 50,
scnb = 2.0, scee = 1.2, nsnb = 25, dielc = 1, cut = 10.0,
ntb = 1,
ntr = 1,
ntc = 2, ntf = 2,
ntx = 1, irest = 0,
imin = 0, nmropt = 0,
ntt=1, tempi=100.0, temp0=300.0, tautp=1.0,
/
Keep protein fixed with weak restraints
5.0
RES 1 36
END
END
This was followed by 2 cycles of 40 ps constant pressure MD
nstlim = 20000, dt = 0.002
ntwe = 500, ntwx = 500, ntpr = 50,
scnb = 2.0, scee = 1.2, nsnb = 25, dielc = 1, cut = 10.0,
ntb = 2,
ntr = 1,
ntc = 2, ntf = 2,
ntx = 1, irest = 0,
imin = 0, nmropt = 0,
ntt=1, tempi=300.0, temp0=300.0, tautp=1.0,
ntp=1, taup=0.2,
/
Keep protein fixed with weak restraints
5.0
RES 1 36
END
END
Then another 40 ps constant volume MD with 5 kcal/(mol A2
nstlim = 20000, dt = 0.002
ntwe = 500, ntwx = 500, ntpr = 50,
scnb = 2.0, scee = 1.2, nsnb = 25, dielc = 1, cut = 10.0,
ntb = 1,
ntr = 1,
ntc = 2, ntf = 2,
ntx = 1, irest = 0,
imin = 0, nmropt = 0,
ntt=1, tempi=300.0, temp0=300.0, tautp=1.0,
/
Keep protein fixed with weak restraints
5.0
RES 1 36
END
END
and 60 ps constant volume MD with 0.5 kcal/(mol A2 )
&cntrl
nstlim = 30000, dt = 0.002
ntwe = 500, ntwx = 500, ntpr = 50,
scnb = 2.0, scee = 1.2, nsnb = 25, dielc = 1, cut = 10.0,
ntb = 1,
ntr = 1,
ntc = 2, ntf = 2,
ntx = 1, irest = 0,
imin = 0, nmropt = 0,
ntt=1, tempi=300.0, temp0=300.0, tautp=1.0,
/
Keep protein fixed with weak restraints
0.5
RES 1 36
END
END
Finally, a 100 ps MD simulation without harmonic constraints
nstlim = 100000, dt = 0.001
ntwe = 500, ntwx = 500, ntpr = 50,
scnb = 2.0, scee = 1.2, nsnb = 25, dielc = 1, cut = 10.0,
ntb = 1,
ntr = 0,
ntc = 2, ntf = 2,
ntx = 1, irest = 0,
imin = 0, nmropt = 1,
ntt=1, tempi=300.0, temp0=300.0, tautp=1.0,
/
&wt type='REST', istep1=0,istep2=100000,value1=1.0,
value2=1.0, /
&wt type='END' /
LISTOUT=POUT
DISANG=RST
After the final MD cycle, water molecules were removed from the
igb = 2, saltcon=0.2, gbsa=1, cut=300.0,
intdiel=1.0, extdiel=78.5,
maxcyc = 2000, ntmin = 1, ncyc = 100, dx0 = 0.01, drms = 0.0001,
ntr = 0, ntpr=10, ntx=1, ntb=0,
imin = 1, nmropt = 1,
/
&wt type='REST', istep1=0,istep2=2000,value1=1.0,
value2=1.0, /
&wt type='END' /
LISTOUT=POUT
DISANG=RST
The trouble was I expected improve the RMSD between energy minimized structures respect the CYANA structures but this was not the outcome, maybe the parameters that I entered were not good values.
On the other side I want
Thanks for your
-----------------------------------------------------------------------
| |||||||||||||||||||||||||||||||||
|