How to perform virtual screening using Autodock Vina?
Virtual screening is used to identify small molecules that are most likely to bind to a target protein. There is various software available for virtual screening including GOLD  and GLIDE . Autodock Vina  is a freely accessible software and provides good results can also be used for screening various ligands. Recently, Autodock Vina has provided plugins  to facilitate virtual screening. In previous articles, we have provided tutorials for the installation of Raccoon and Raccoon2 plugins on Ubuntu.
These plugins provide a graphical user interface and there is no need to prepare pdbqt files of protein and ligands. You can simply load your target protein and ligand pdb files. In this article, we will use a Perl script or a bash script to screen several ligands for a target protein.
Preparing input files
The docking procedure will be the same as mentioned in previous articles on site-specific docking and blind docking. (Video tutorials: site-specific docking and blind docking).
Prepare pdbqt files of the target protein and all ligands. You can also use OpenBabel software to convert them from sdf into pdbqt format. If you have downloaded smiles and want to convert them into 3D conformers then again use OpenBabel for that. For more details, read this article. Save them into the same directory.
For Ubuntu/Linux, download the following files from our Bitbucket account of Bioinformatics-Review or my personal account MunibaFaiza:
Perl script (vina_vs_linux.pl) and bash script (vina_vs.sh).
If you are using Perl script, then list the names of prepared pdbqt files of ligands in a text file.
If you are using the bash script, save the ligands pdbqt files as “ligand_name.pdbqt”, where the name could be any letter. For example, “ligand_3.pdbqt” or “ligand_hydrocl.pdbqt”, and so on.
Further, the configuration file will remain the same as mentioned in previous articles except for the second line, i.e., “Ligand = ligand.pdbqt”. You can download the configuration file from here (change the coordinates as per your grid dimensions).
For Windows, download the following file from here or here:
Perl script (vina_vs_win.pl).
Running Vina for virtual screening
Open the terminal (Ctrl+Alt+T) and type the following command:
$ perl vina_vs_linux.pl
It will prompt for the name of a file containing all the names of ligands. Enter the name of the file and continue.
Open command prompt (type cmd in the search tab) and navigate to the folder where you have kept all files including vina.exe. Type the following command:
It will also prompt for the name of the file containing all the names of ligands. Enter the name of the file and continue.
> bash vina_vs.sh
Make sure you have already installed and added Vina in your path on Ubuntu. If not, then follow this tutorial.
- Jones, G., Willett, P., Glen, R. C., Leach, A. R., & Taylor, R. (1997). Development and validation of a genetic algorithm for flexible docking. Journal of molecular biology, 267(3), 727-748.
- Friesner, R. A., Banks, J. L., Murphy, R. B., Halgren, T. A., Klicic, J. J., Mainz, D. T., … & Shaw, D. E. (2004). Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. Journal of medicinal chemistry, 47(7), 1739-1749.
- Trott, O., & Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. Journal of computational chemistry, 31(2), 455-461.
- Forli, S., Huey, R., Pique, M. E., Sanner, M. F., Goodsell, D. S., & Olson, A. J. (2016). Computational protein-ligand docking and virtual drug screening with the AutoDock suite. Nature Protocols, 11(5), 905-919.
What values are considered as good or bad in computational docking?
After performing computational docking, a question that comes to mind most is “what docking score is considered good or bad”. In this article, we will discuss this in detail. (more…)
How to sort binding affinities based on a cutoff using vs_analysis.py script?
Previously, we have provided a Python script (vs_analysis.py) to analyze the virtual screening (VS) results of Autodock Vina. Now, we have updated this script to sort binding affinities based on user inputted cutoff value. (more…)
Virtual Screening using Autodock Vina: Frequently Asked Questions & Answers for Starters
Virtual Screening (VS) is one of the important techniques in bioinformatics. It can be easily performed using Autodock Vina. We have provided detailed articles on this topic. In this article, we are trying to answer some FAQs for beginners. (more…)
July 5, 2021 at 12:08 pm
How can I install Raccoon2 plugins on windows?
Dr. Muniba Faiza
July 12, 2021 at 3:33 pm
Please visit the following link: