Skip to main content

AutoDock Vina protocol: For Molecular Docking

There are many softwares available for carrying out molecular docking calculations. Many of these softwares are paid, while some of them are free. AutoDock is one such freely available software. Earlier, there was AutoDock 4.0 version available. One of the problem with the software was the space it would occupy, if you attempt to do multi docking.
Recently, there is another algorithm made available under autodock. It is called AutoDock Vina.
This is more efficient, robust and fast algorithm.
There are many tutorials available to learn AutoDock Vina including:
http://autodock.scripps.edu/faqs-help/tutorial/using-autodock-4-with-autodocktools/2012_ADTtut.pdf

and also a video tutorial
http://vina.scripps.edu/tutorial.html 

 But for some reason, few changes or explanation in little detail is always required to make it understandable for a beginner to follow. 
Here is a protocol designed in our group, which our Master's student follow to get acquainted to Molecular Docking studies.
This protocol can be pretty much followed for any molecular docking study:
----------------------------------------------

*You need to download files for ligand and target form the link as suggested in the tutorial
Or else you may directly attempt to work with your set of ligands and targets. 

Step I: Preparation of Ligand PDBQT files:


1)       Open AutoDockTools.
2)      On the Menu bar Click on File --- Read Molecule --- Desktop---Docking_Files---ind.pdb.
3)      Press Open
Below Menu bar you will find an option named Ligand.
4)      Click on Ligand --- Input --- Choose --- ind
5)      Press  Select Molecule for AutoDock 4.
A widget appears on the screen describing the chemical nature of the Drug Molecule. Press OK
6)      Click on Ligand --- Output --- Save as PDBQT.
7)      Save the PDBQT file on Desktop in Docking_Files folder.
In the Dashboard Click on the red circle (fig.1) beside the name of the drug molecule ind.
Drug molecule will disappear from the screen.


 


Step II: Preparation of  Protein PDBQT files.
1)      On the Menu bar Click on File --- Read Molecule --- Desktop---Docking_Files---1hsg.pdb.
2)      Press Open
3)      On the Menu Bar, click on Select --- Select From String.
     A Dialogue Box appears on the screen  with the name Select From String
4)      Type HOH* within the Residue Box.
5)      Type * within the Atom Box (fig.2).


1)       Click Add and then Press Dismiss.
2)       On the Menu Bar click on Edit --- Delete Water
This step deletes all water molecules from the PDB files.
Below the Menu bar you will find an Option named Grid.
3)      Click on Grid --- Macromolecule --- Choose ---- hsg1.
4)       Press Select Molecule.
A WARNING widget appears on the screen.
5)       Click on OK.
6)      A Dialogue Box appears on screen for saving the PDBQT.
7)      Save this file within the folder named Docking_File on Desktop as hsg1.pdbqt.
8)      Click on Grid --- Grid Box
Change the X,Y and Z dimension to 26.
Enter the values 3.18, -0.294 and -3.784 in x, y and z center respectively.
9)      Increase Spacing from 0.375 to 1 (fig.3).



Before Closing this widget, open notepad and enter the following.
receptor = hsg1.pdbqt
ligand = ind.pdbqt


center_x = 3.18
center_y = -0.294
center_z = -3.784

size_x = 26
size_y = 26
size_z = 26

1)      Save this file as conf.txt (fig.4).

 

1)       In Grid Option widget, Click on File--- Close Saving Current
Close AutodockTools.
Before starting with Docking check whether  hsg1.pdbqt, ind.pdbqt and conf.txt files are all saved in the Docking_Files folder.

Step III: Running Autodock Vina.
1)      Press start button and click on Command Prompt
2)      In the command prompt type cd, press space and drag the folder Docking_Files over the command prompt, Press enter .
3)      Open Docking_Files folder drag the file named Vina over the command prompt, press space and type --config, again press space then drag the conf.txt file over the command prompt.
4)      Press Space and type --log log.txt. Press enter.
Command Prompt should be similar to fig.5

 

Wait for the Docking Process to complete. A file with the name ind_out.pdbqt will be generated in the folder Docking_Files folder

Step IV: Analysing Vina Result.
1)      Open AutodockTools.
2)      Click on Analyze --- Docking --- Open AutodockVina Result --- Desktop --- Docking_Files --- ind_out.pdbqt.
Click on the option Single molecule with muliple conformations and press OK.
3)      Click on Analyze --- Macromolecule --- Open --- Desktop --- Docking_Files --- hsg1.pdbqt.
4)      Click in Analyze --- Docking --- Show Interaction.
Study the interactions between the protein and ligand. Check for hydrogen bond interactions in different ligand conformers. Use left and right arrow keys on the Keyboard to change conformations.
Close AutodockTools.
5)      Open DiscoveryStudio from Desktop Icon.
6)       Drag hsg1.pdbqt , then drag ind_out.pdbqt.
7)      Click on Scripts --- Ligand Interactions --- Show ligand binding site atoms.
8)      Below the Menu Bar Click on Receptor-Ligand Interactions.
9)      In the Tools explorer widget click on Hydrophobic.
10)  Right click on Black window Labels--- Object --- from the Drop Down Menu select Amino Acids and press OK.

Repeat the Steps I, III and IV for other ligand file. Remember to  change the ligand name in the conf file accordingly.


Comments

Popular posts from this blog

Material Studio - selective molecule movement

So, working with Material Studio is awesome. You get to handle bigger, complicated and even unimaginable system. But I am facing problems in very minute functioning. I recently faced two problems, but managed to resolve on my own. Somehow, tutorials doesn't help in resolving them and you end up doing some exercise on your own to get the solution. Among these, the problem of bringing two molecular structures in one window and then adjusting one molecule by moving it around, while other structure remains untouched, as it is. Actually, this should not be a problem to do. There is a simple way of doing this. Select the molecule to be moved and then there is a symbol (button) to select, which allows you to selectively move that molecule only. This is the same way it can be done in Material studio too.  But among all the dozens of symbols, the new person handling MS, kind of gets lost. So, what he does is check the tutorials, for any help. But the tutorials don't talk about suc...

How would COVID bahave in August - few projections

Entering August with COVID by side   Entering in the month of August, COVID is still here by our side. Situation seems grim, worse than earlier, as number of active cases as well as the count of deceased patients are increasing day by day. Yesterday we crossed the mark of 50,000 cases in a day. Closely competing with Brazil with US counts also in range now. On the other side, few other numbers like reduced mortality rate to 2.18%, recovery rate of 64.6% and slightly increased doubling rate of 21 days can reduce the anxiety. My projections of 16.75 lakhs for India might be crossed by coupe of thousands of more cases after today’s count is added. This number would reach to a huge count of 28.5 lakhs by August 15. Maharashtra just doesn’t like competition. Just like what Delhi tried few days back, Andhra Pradesh also managed to cross Maharashtra’s per day count for2 consecutive days. Seems like not a good idea. Maharashtra immediately went to the next zone. The state which w...

Rigorous method to create lib file

Sometimes you totally fail to go ahead with the prepin and frcmod files you have generated from the newly generated PDB file even after you follow the procedure given in previous two posts. And you keep wondering? In such a case you need a more detailed procedure. First take a look at the pdb file. For example: Here is one such pdb file generated from the gabedit that troubled me for some time. --------------- HEADER    PROTEIN COMPND    UNNAMED AUTHOR    GENERATED BY GABEDIT 2.3.8 at Tue Jul 26 16:53:48 2011 ATOM  1      C   C     1       -2.907  -0.997  -0.385  1.00  300.00          C 0.0000   ATOM  2      N   N     1       -2.060  -0.981  -1.189  1.00  300.00 ...