Hi, I am trying to perturb a molecule to nothing in a protein binding
pocket. By looking at the pdb files at the end of each lambda run, the
molecule has migrated out of the protein when clambda is equal to
0.43738 or bigger. One explanation would be that the small molecule
with shrinking vdw radii do not have favorable interactions with
protein any more comparing to water molecules and expelled therefore.
Therefore I wonder if a weak restraint is needed in order to evaluate
the dV/dL ~ L integral. Thanks.
Eric
&cntrl
ntr=0,
nstlim =20000, nscm=2000, ntave=5000,
ntx=1, irest=0, ntb=1, ntpr=100, tempi=300.0, ig=974651,
ntp=0, taup=1.0,
dt=0.001, nrespa=1,
ntt=1, temp0 = 300., tautp=2.0,
ntc=2, ntf=2, tol=0.000001,
ntwr = 100, ntwx=100,
icfe=1, klambda=6, clambda=0.99078,
cut=9.0,
/
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Mar 18 2005 - 23:53:00 PST
