I am trying to compute a free energy of binding of a ligand to a large
protein using mmpbsa and nmode analysis with Amber 9.
My question is that i see in many publications that a 10-15A region of
the system is cut out around the ligand, and the subsystem used to
obtain an estimate of dS. How is this done using amber? Is there a way
to tell the leap program not to add hydrogens to the termini when
generating a .top file? There will be many segments when the region is
cut, and I do not want amber to patch the termini with hydrogens as the
added atoms will not match the coordinates of the trajectory file I have
previously created.
Any suggestions on generating the new .top file for the subsystem
without introducing extra atoms?
Thanks,
Matt Danielson
Purdue University
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jun 17 2007 - 06:07:37 PDT
