top_hits.py to Analyse AutoDock-Vina-GPU results to screen ligand library. #10
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
To analyse the data when number of ligands are higher. Please use this python script.
First while running the
AutoDock-Vina-GPU-2.1 by following Command, ./AutoDock-Vina-GPU-2-1 --config ./hnod2_input/config.txt 2>&1 | tee result.txt
Here,
2>&1 | tee result.txt
will dump all the data in result.txt file.Later, using this python script we can analyses all the data in result.txt file to list top hits. specify the path to result.txt file in the top_hits.py file (The name can be anything you provided in last step).
Run this script using command,
python3 top_hits.py
It will ask user to give number of hits wanted (lets say we want 20 top hits) So give 20 to next line in terminal.
It will list all the ligands.
**My ligands were named from 1.pdbqt to 10000.pdbqt. **I got 32 GB RAM hence was able to do 10000 ligands at time for approx 1.26 hours.