Connect with us

Bioinformatics Programming

sminalog_analysis.py – A new Python script to fetch top binding affinities from SMINA log file

Dr. Muniba Faiza

Published

on

sminalog_analysis.py – A new Python script to fetch top binding affinities from SMINA log file

In one of our previous posts, we provided a Python script for the virtual screening analysis of Autodock Vina. This script analyzes all log files obtained from docking of multiple ligands to a receptor and provides the binding affinities for top poses from each file. In this article, we are publishing a new Python script for the virtual screening analysis of SMINA [1].

SMINA is a fork of Autodock Vina that improves scoring and minimization. It provides a single log file as output which contains all the binding affinities of the docked ligands. However, it does include the name or identity of the docked ligands. Therefore, it is difficult to identify which binding affinity belongs to which ligand. For this purpose, an SDF file as output does the trick. You can use sdsorter to fetch the binding affinities including the name or identity of the associated ligands. Here is a tutorial to install sdsorter on Ubuntu.

SMINA log file can be useful in the conditions where you need only binding affinities of the ligands regardless of their identities. One such example is the AUC-ROC analysis, calculation enrichment factor (EF) values, and so on.

sminalog_analysis.py script fetches all binding affinities of top poses of each ligand giving a serial number to each one of them from the SMINA output log file.

Usage

Download the file and save it in the directory where you have kept all log files. If you will save it elsewhere, don’t forget to provide the full path to this file. It will run on Linux as well as on Windows.

$ python sminalog_analysis.py

After running the above command, you will see the directory path where you have kept the file along with its name. A few seconds later, it will display “Done! The result is provided in the output.txt file”. A new file will be created namely, ‘output.txt’, in the same directory. This file will contain the final results.

Example:

Say, you have docked 3000 ligands from a database and want to get the binding affinities of the top poses of each of these ligands. Then run the above command and wait for the output file. The output will look like this:

sminalog_analysis.py

Availability

The script is freely available to download from the Bitbucket account of Bioinformatics Review.

NOTE:

  • This script screens for the log files containing the word ‘log’ in their filenames, e.g., “log_docking.txt” or “log_3YPI.txt”, or “2yp1_log_docking.txt” etc.

How to cite this script:

Faiza M., Abdullah T. (2021). sminalog_analysis.py – A Python script to fetch top binding affinities from SMINA log file. 8(6):page 16-18. The article is available at https://bioinformaticsreview.com/20210611/sminalog_analysis-py-a-new-python-script-to-fetch-top-binding-affinities-from-smina-log-file

For any queries, contact me at [email protected].


References

  1. Koes, D. R., Baumgartner, M. P., & Camacho, C. J. (2013). Lessons learned in empirical scoring with smina from the CSAR 2011 benchmarking exercise. Journal of chemical information and modeling53(8), 1893-1904.

Further Reading

vs_Analysis.py: A Python Script to Analyze Virtual Screening Results of Autodock Vina

How to install sdsorter on Ubuntu (Linux)?

Virtual Screening Methodology for Structure-based Drug Designing

Dr. Muniba is a Bioinformatician based in New Delhi, India. She has completed her PhD in Bioinformatics from South China University of Technology, Guangzhou, China. She has cutting edge knowledge of bioinformatics tools, algorithms, and drug designing. When she is not reading she is found enjoying with the family. Know more about Muniba

Bioinformatics Programming

Installing Pycharm on Ubuntu (Linux)

Tariq Abdullah

Published

on

Installing pycharm on Ubuntu

Pycharm [1] is an integrated development environment (IDE) for developers. It combines Python developer tools and provides an easy graphical user interface. In this article, we are going to install Pycharm on Ubuntu. (more…)

Continue Reading

Algorithms

vs_Analysis.py: A Python Script to Analyze Virtual Screening Results of Autodock Vina

Dr. Muniba Faiza

Published

on

VS-Analysis: A Python Script to Analyze Virtual Screening Results of Autodock Vina

The output files obtained as a result of virtual screening (VS) using Autodock Vina may be large in number. It is difficult or quite impossible to analyze them manually. Therefore, we are providing a Python script to fetch top results (i.e., compounds showing low binding affinities). (more…)

Continue Reading

Bioinformatics Programming

How to execute Unix/shell commands in PHP?

Dr. Muniba Faiza

Published

on

programming in php

Besides, Python & Perl, PHP is another language mostly used in bioinformatics programming. In this article, we provide a simple command to execute shell commands in a PHP script. (more…)

Continue Reading

LATEST ISSUE

ADVERT