The analysis of Autodock Vina  results is a bit tricky in the sense of viewing all interactions and selecting the best pose. In our last video tutorial, we explained how to analyze docking results obtained from Vina using PyMol. This article is the written guide for the same.
We need a PDB file of protein and vina output file in .pdbqt format. Here, our protein file is 2bxa.pdb and the vina output file is cmpf.pdbqt. Let’s get started!
Open PyMol then go to 'File' --> 'Open' --> then select the PDB file of your protein (here, 2bxa.pdb) --> click 'Open'.It will display the protein.
- Now open the vina output file. Again
go to 'File' --> 'Open' --> then select the pdbqt file (here, cmpf.pdbqt) --> click 'Open'.
You can now see the protein and the ligand but there are no bondings/interactions between the two. Now, in order to see the interactions, follow these steps:
Look at the middle right side in PyMol window, there you can see some tabs named as ‘all’, ‘2bxa.pdb’, and ‘cmpf.pdbqt’, and some letters such as ‘A’ (Action), ‘S’ (Show), ‘H’ (Hide), and ‘L’ (Label) right in front of these tabs.
Go to the tab 'all' --> click 'A' (it will show you a drop-down menu) --> 'preset' --> 'ligand sites' --> 'cartoon'.
It will re-center your complex and show the interactions the protein is making with the ligand.
You can also measure the bond lengths following these steps:
Go to 'Wizard' --> 'Measurement'.
It will ask to select the atoms.
Select two atoms from both the protein and the ligand showing interaction. It will display the length.
Similarly, you can label the residues as well by going to
'Wizard' --> 'Label'
Now, there comes the point where you need to choose the best docking pose. Well, the best pose is supposed to be the first docking pose generated by Vina which shows ‘zero’ RMSD value and best binding affinity. You may also need to make sure whether your ligand is interacting to all the binding/catalytic residues you defined for docking. Sometimes, your second pose may also show the required maximum interactions within the binding pocket.
This is the basic introductory tutorial, you can always explore more to represent the complex more beautifully. There are several ways to represent the protein and the ligand such as ribbon, surface, cartoon, and more. You can also show only binding pockets of the protein as per your requirements.
- 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.
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…)