Journal
JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 15, Issue 11, Pages 5817-5828Publisher
AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.9b00262
Keywords
-
Ask authors/readers for more resources
Molecular dynamics (MD) simulations can be a powerful tool for modeling complex dissociative processes such as sed ligand unbinding. However, many biologically relevant dissociative processes occur on timescales that far exceed the timescales of typical MD simulations. Here, we implement and apply an enhanced sampling method in which specific energy terms in the potential energy function are selectively scaled to accelerate dissociative events during simulations. Using ligand unbinding as an example of a complex dissociative process, we selectively scaled up ligand-water interactions in an attempt to increase the rate of ligand unbinding. After applying our selectively scaled MD (ssMD) approach to several cyclin-dependent kinase-inhibitor complexes, we discovered that we could accelerate ligand unbinding, thereby allowing, in some cases, unbinding events to occur within as little as 2 ns. Moreover, we found that we could make realistic estimates of the initial unbinding times (tau(sim)(unbind)) as well as the accompanying change in free energy (Delta G(sim)) of the inhibitors from our ssMD simulation data. To accomplish this, we employed a previously described Kramers'-based rate extrapolation method and a newly described free energy extrapolation method. Because our ssMD approach is general, it should find utility as an easy-to-deploy, enhanced sampling method for modeling other dissociative processes.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available