Version 1 of the repository for the Unbinding-MLTSA, this code corresponds to the bioRxiv preprint: https://www.biorxiv.org/content/10.1101/2021.09.08.459492v1
Combined free energy calculation and machine learning methods for understanding ligand unbinding kinetics. Magd Badaoui, Pedro J Buigues, Dénes Berta, Gaurav M. Mandana, Hankang Gu, Callum J Dickson, Viktor Hornak, Mitsunori Kato, Carla Molteni, Simon Parsons, Edina Rosta bioRxiv 2021.09.08.459492; doi: https://doi.org/10.1101/2021.09.08.459492
Under "MLTSA_examples/OneD_pot/" the file "OneD_sklearn_example.ipynb" is a Jupyter Notebook going through the basics of generating the one dimensional data and subsequent features for the ML. It contains code for the training, testing and global mean calculation as well as the Accuracy Drop for feature selection.
You can run this example using
-
"MLTSA_sklearn/MLTSA_sk.py" : Contains the main code for calculating the Global Means and subsequent Accuracy Drop calculation as well as a wrapping function for plotting nicely the outcome of the approach.
-
"MLTSA_sklearn/models.py" : Contains a wrapping function for an easy train/test of almost any Sklearn classifier.
-
"MLTSA_examples/OneD_pot/OneD_pot_data.py" : Contains the main code for generating simulations on 1D potentials (Single-well and Double-well) and prepare a dataset with the desired number of features. It has functions to generate data on demand based on the
GIF files for each system are available under the "Animated_Trajectories" folder.
Otherwise see here: