Dear amber users,
I performed NPT simulation (NTT=1) for a dimmer solvated in a truncated octahedral box, with a minimum distance of 10 angstrom from the atoms of protein to the edge of the solvent box. After about 5ns NPT simulation, the protein protruded from the unit cell. It seems that the protein didn't rotate around its center (geometry or center of mass), which caused the problem. Can the problem be avoided by translating and orienting center of the initial dimmer model to overlap with the origin of the coordinate set and then solvating it in the OCT box? Or this problem is due to the parameters used (DT, NSCM etc.)?
Another surprising phenomena is that the ligand in one protein of the dimmer run away from the binding pocket while the other one kept in the original position of the other protein in the dimmer.
Here is the script I used in PMEMD:
---------------------
NPT production
&cntrl
imin=0,
ntx=5, irest=1,
ntpr=100, ntwx=1000,ntwr=5000,
dt=0.002,nstlim=1000000, nscm=5000,
cut=9.,dielc=1.0,
ntc=2, ntf=2,
ntb=2,
ntp=1,pres0=1.0,taup=3.0,
ntt=1, tautp=3.0,temp0=300.0,
ioutfm=1
/
-----------------------
Any suggestion is greatly appreciated.
Thanks.
----
Jeffrey
嗃i?櫒ド╤∈&
Received on Fri Mar 20 2009 - 01:07:14 PDT
