More FEP using separated topologies

From: David Huggins (djh210_at_cam.ac.uk)
Date: Fri May 27 2016 - 15:21:01 CDT

Dear All,

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
other).

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?

Best wishes,

Dave

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Dr David J. Huggins
MRC New Investigator
Theory of Condensed Matter Group
Cavendish Laboratory
JJ Thomson Avenue
Cambridge, CB3 0HE
United Kingdom

Phone: +44 (0)1223 764164
Fax: +44 (0)1223 337356
Email: djh210_at_cam.ac.uk
Web: http://www.tcm.phy.cam.ac.uk/~djh210/

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

This archive was generated by hypermail 2.1.6 : Tue Dec 27 2016 - 23:22:12 CST