From: David Huggins (djh210_at_cam.ac.uk)
Date: Fri May 27 2016 - 15:21:01 CDT
I am trying to implement a NAMD version of the separated topologies
method from "Separated topologies—A method for relative binding free
energy calculations using orientational restraints".
I have used alchemify to switch off interactions between two ligands
in a protein binding site and the protocol requires me to run dynamics
where both ligands are fully coupled to the system (though not to each
However, the resulting simulations always exit early with a "Low
global exclusion count". I have tried gentle heating and reduced
timesteps but neither works.
Does anyone have any experience of running simulations where two fully
coupled species are overlapping and have no mutual interactions?
Dr David J. Huggins
MRC New Investigator
Theory of Condensed Matter Group
JJ Thomson Avenue
Cambridge, CB3 0HE
Phone: +44 (0)1223 764164
Fax: +44 (0)1223 337356
This archive was generated by hypermail 2.1.6 : Tue Dec 27 2016 - 23:22:12 CST