[WIP] Analysis tools for identifying ligand binding modes from MD or BLUES simulations.#113
[WIP] Analysis tools for identifying ligand binding modes from MD or BLUES simulations.#113nathanmlim wants to merge 31 commits intomasterfrom
Conversation
… pose sampling plot.
|
This will need some changes with upcoming changes to our NCMC reporter. |
|
This PR looks quite useful—it adds a modular analysis toolkit (blues.analysis) for identifying metastable ligand binding modes using PyEMMA. It’s well-scoped and does not alter the core functionality of BLUES, making it a safe and valuable addition for post-simulation analysis. For future development, it may be worth refactoring to support deeptime ( https://deeptime-ml.github.io/latest/index.html ) , given the deprecation of PyEMMA. I did encounter some merge conflicts with the following files: These are all straightforward and can be resolved by merging both sets of changes. |
This PR adds in an analysis module I've been working on that automatically identify metastable ligand binding modes using tools from PyEMMA. Included is a juypter notebook that provides a brief walk through usage of the analysis modules, located in
notebooks/example-id_bmodes.ipynbThe analysis modules can be accessed like
blues.analysis import msm,cluster. Helper functions for these modules are located intools.pymsmsubmodule contains theConstructMSMclass:clustersubmodule contains theFindBindingModesclass:populationsubmodule contains theBindingModeOccupancyclass: