Re: Question on making 'reversible' solvated simulation

From: Randy J. Zauhar (r.zauhar_at_usciences.edu)
Date: Sun Feb 11 2018 - 12:52:18 CST

Jerome, thanks - I see the issue of site de-solvation not being properly included. ;-(

One thing I will check is my protein preparation - it is possible that I removed crystallographic active site waters, and the ones that crept in are in fact occupying similar positions. In that case I could simply include the xtal waters. That would at least improve things.

Another would be to retain the ‘solvent wall’, and introduce some active site waters in place of the ligand. Could then do a separate simulation to annihilate the active site waters, which would also provide some measure of de-solvation by the ligand.

Ultimately we want to compare ligands for the same target, and the desolvation of the site should be similar for each.

Randy

P.S. I will take a look at the work by Yang, but do not contemplate moving to CHARMm over this!

On 11Feb, 2018, at 5:56 AM, Jérôme Hénin <jerome.henin_at_ibpc.fr<mailto:jerome.henin_at_ibpc.fr>> wrote:

Hi Randy,

While adding a barrier to entry for water molecules may well make the simulation more reversible, it will also make it biased and give an overestimated binding affinity, by killing the binding site desolvation contribution.

It's a tough problem, but ideally this solvation/desolvation process needs to be sampled.

Improving the sampling of that kind of "orthogonal" process was the focus of some work by Wei Yang:
http://www.pnas.org/content/105/51/20227.full
https://pubs.acs.org/doi/abs/10.1021/ct200726v

But that's only implemented in CHARMM afaik. The theory is somewhat involved.

So I'm sorry to say, if it's a common problem, there is no universal, simple and consensual approach to it.

Best,
Jerome

On 10 February 2018 at 19:15, Randy J. Zauhar <r.zauhar_at_usciences.edu<mailto:r.zauhar_at_usciences.edu>> wrote:
Hi,

I am running a ligand annihilation in a solvated system, and I would like to make it truly reversible if possible.

Issue is that a couple water molecules sneak into the active site when the ligand disappears. Then when I reverse and start creating the ligand, they are trapped.

A simple thing I will try is to constraint the water molecules close to the active site, hopefully creating a barrier to entrance.

However. this has to be a pretty common problem, is there a recommended approach to address it?

Thanks!
Randy

Randy J. Zauhar, PhD

Prof. of Biochemistry

Dept. of Chemistry & Biochemistry
University of the Sciences in Philadelphia
600 S. 43rd Street
Philadelphia, PA 19104

Phone: (215)596-8691
FAX: (215)596-8543
E-mail: r.zauhar_at_usciences.edu<mailto:r.zauhar_at_usciences.edu><mailto:r.zauhar_at_usciences.edu<mailto:r.zauhar_at_usciences.edu>>

“Yeah the night is gonna fall, and the vultures will surround you /
And when you’re lookin’ in the mirror what you see is gon’ astound you"

  — Death Cab for Cutie, “Monday Morning"

Randy J. Zauhar, PhD

Prof. of Biochemistry

Dept. of Chemistry & Biochemistry
University of the Sciences in Philadelphia
600 S. 43rd Street
Philadelphia, PA 19104

Phone: (215)596-8691
FAX: (215)596-8543
E-mail: r.zauhar_at_usciences.edu<mailto:r.zauhar_at_usciences.edu>

“Yeah the night is gonna fall, and the vultures will surround you /
And when you’re lookin’ in the mirror what you see is gon’ astound you"

  — Death Cab for Cutie, “Monday Morning"

This archive was generated by hypermail 2.1.6 : Tue Dec 31 2019 - 23:19:41 CST