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
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.
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).
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.
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…)
[Tutorial] How to perform docking of zinc metalloproteins using Autodock Vina?
Proteins containing zinc atoms are docked in a different way than that of the normal simple proteins. Zinc atoms must be considered by a force field during the docking process. In this article, we are going to dock zinc metalloprotein with a ligand using Autodock Vina . (more…)