Hi everyone,
I have a protein-ligand docked complex and I'm trying to apply energy
minimization to take care of the steric clashes and hopefully optimize
ligand position and geometry. Currently I'm running the minimization
in gas phase and the input looks like:
&cntrl
ntxo = 1,
ntf = 1, ntb = 0,
dielc = 4,
cut = 12.0, nsnb = 99999,
scnb = 2.0, scee = 1.2,
imin = 1, maxcyc = 1500,
ncyc = 500, drms = 0.0001
ibelly = 1, bellymask=':ligname | (:ligname < :12.0)'
&end
&ewald
eedmeth= 5,
&end
My question is if it's better off to use the GB model (for example
igb=1) than gas phase for that purpose?
I also want to run MM-GBSA (w/ igb=1) on this complex, so would GB
minimization offer better consistency and results if the same
parameters (dielc, intdiel, extdiel, saltcon etc.) are used?
Could anyone kindly give me an example GB minimization input with
restraint that can be used to substitute the above?
Thanks a lot in advance,
-DX
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Wed Jul 30 2008 - 06:07:35 PDT
