From: JC Gumbart (gumbart_at_ks.uiuc.edu)
Date: Mon Jun 17 2013 - 11:37:10 CDT
Take a look at this paper and see if addresses your questions: http://pubs.acs.org/doi/abs/10.1021/ct3008099
On Jun 17, 2013, at 10:43 AM, Felipe Merino wrote:
> Dear all,
> I know this is a little bit off topic but i think somebody could help us here. We have been doing some umbrella sampling simulations to calculate the binding free energy of a protein-DNA complex. We are using the minimal interatomic distance as reaction coordinate. The thing is that the protein has two domains and we are only pulling them separately, so in the end we have always an endpoint with the other domain still bound (that was planned). In the end, the binding free energies are much higher than expected (around 25 kcal/mol). The point is that it sounds to me that there should be a correction on the PMF to account for the "confinement" of being bound to the other (still attached) domain (very much line in the case where you restrain the ligand when doing FEP annihilation) but i am not sure if this is correct. Any insights for this will be highly appreciated.
This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:23:21 CST