Lesson · 40 min · Free
Molecular Docking Lab with AutoDock Vina
Molecular Docking Lab with AutoDock Vina Molecular Docking Lab with AutoDock Vina Welcome to the practical laboratory session for Molecular Docking with AutoDock Vina! In this lab, we will walk through the essential step
Molecular Docking Lab with AutoDock Vina
Welcome to the practical laboratory session for Molecular Docking with AutoDock Vina! In this lab, we will walk through the essential steps of preparing your ligand and receptor files, setting up the AutoDock Vina calculation, and interpreting the results. This hands-on experience is crucial for understanding the theoretical concepts we've covered and for gaining practical skills in computational drug discovery. Before we begin, ensure you have AutoDock Tools (ADT) and AutoDock Vina installed on your system. We will primarily use ADT for preparing our input files and Vina for performing the docking simulation. The goal is to dock a known ligand into a protein target and analyze the predicted binding poses and affinities.
Setting Up Your Docking Experiment
The first critical step is to prepare your macromolecule (receptor) and ligand files in the correct format. AutoDock Vina requires input files in PDBQT format. AutoDock Tools provides a user-friendly graphical interface to convert standard PDB files into PDBQT, add necessary hydrogen atoms, assign Gasteiger charges, and define rotatable bonds for the ligand.
1. Receptor Preparation using AutoDock Tools
Open AutoDock Tools (ADT). Go to File > Read Molecule and load your protein PDB file (e.g., 1iep.pdb ). Once loaded, you need to prepare it for docking: Remove Water Molecules: Go to Edit > Delete > Water . Add Hydrogens: Go to Edit > Hydrogens > Add . Add Gasteiger Charges: Go to Edit > Charges > Add Gasteiger . Merge Non-polar Hydrogens: Go to Edit > Hydrogens > Merge Non-polar . Set Torsions (Optional, but good practice): Go to Ligand > Torsion Tree > Detect Root and Set Up Torsions . This step is more critical for ligands but doesn't hurt for receptors. Save as PDBQT: Go to File > Save > Write PDBQT . Save it as receptor.pdbqt .
2. Ligand Preparation using AutoDock Tools
Similarly, prepare your ligand. If your ligand is in a 2D or 3D format other than PDB (e.g., SDF, MOL2), you might first need to convert it to PDB using other software (like OpenBabel) or directly load it if ADT supports the format. Load Ligand: Go to Ligand > Input > Open . Load your ligand PDB file (e.g., ligand.pdb ). Add Hydrogens: Go to Ligand > Hydrogens > Add . Add Gasteiger Charges: Go to Ligand > Charges > Compute Gasteiger . Detect Torsions: Go to Ligand > Torsion Tree > Detect Root and Set Up Torsions . This step is crucial for defining the flexible degrees of freedom for the ligand during docking. Save as PDBQT: Go to Ligand > Output > Save as PDBQT . Save it as ligand.pdbqt .
3. Defining the Search Space (Grid Box)
The grid box defines the region where AutoDock Vina will search for optimal ligand poses. It's crucial to center this box around your active site or a region of interest. In ADT, after loading your receptor (and optionally, a reference ligand if you know the binding site): Go to Grid > Macromolecule > Choose and select your prepared receptor.pdbqt . Go to Grid > Grid Box . A green box will appear. Adjust the Center X, Y, Z coordinates and Size X, Y, Z dimensions to encompass your active site. You can manually type values or drag the box in the viewer. A common practice is to center it on a known ligand or a set of active site residues. Note down these coordinates and dimensions. You will need them for the Vina configuration file.
4. Creating the Configuration File (conf.txt)
AutoDock Vina uses a simple text file to specify all docking parameters. This file includes the paths to your receptor and ligand, the grid box dimensions, and optional parameters like exhaustiveness. receptor = receptor.pdbqt ligand = ligand.pdbqt center_x = 15.0 center_y = 25.0 center_z = 30.0 size_x = 20.0 size_y = 20.0 size_z = 20.0 out = docking_results.pdbqt log = vina_log.txt exhaustiveness = 8 num_modes = 9 energy_range = 3 Explanation of parameters: receptor : Path to your prepared receptor PDBQT file. ligand : Path to your prepared ligand PDBQT file. center_x, center_y, center_z : Coordinates of the center of your grid box (obtained from ADT). size_x, size_y, size_z : Dimensions of your grid box in Ångstroms (obtained from ADT). out : Output file name for the docked poses (will be a PDBQT file containing multiple poses). log : Output file name for the Vina log, including binding affinities. exhaustiveness : Controls the thoroughness of the search. Higher values (e.g., 16-32) lead to more accurate results but take longer. Default is 8. num_modes : The maximum number of binding modes to generate. Default is 9. energy_range : The maximum energy difference between the best binding mode and the worst binding mode displayed. Default is 3 kcal/mol.
5. Running AutoDock Vina
Navigate to the directory containing your receptor.pdbqt , ligand.pdbqt , and conf.txt files using your terminal or command prompt. Then execute AutoDock Vina: vina --config conf.txt Vina will run, displaying progress updates in the terminal. Once finished, it will output the docking_results.pdbqt and vina_log.txt files.
6. Analyzing Results
Open vina_log.txt to see the predicted binding affinities (in kcal/mol) for each generated pose. Lower (more negative) values indicate stronger predicted binding. To visualize the docked poses, open docking_results.pdbqt in AutoDock Tools ( File > Read Molecule ). You can then load your receptor ( File > Read Molecule ) and visualize the ligand poses within the active site. ADT allows you to cycle through the different poses and analyze their interactions (hydrogen bonds, hydrophobic interactions, etc.) with the receptor.
Key Takeaways:
Proper preparation of receptor and ligand PDBQT files is paramount for successful docking. The grid box must accurately cover the active site to ensure relevant binding poses are found. The configuration file ( conf.txt ) allows fine-tuning of Vina's parameters. Binding affinities are reported in kcal/mol, with more negative values indicating stronger binding. Visualization tools like AutoDock Tools are essential for analyzing predicted poses and interactions.
Practice Exercise:
Choose a protein target from the PDB database (e.g., a common drug target like DHFR or HIV-1 Protease) and download its PDB file. Identify a known ligand bound to this protein. Prepare both the protein and the ligand for AutoDock Vina using AutoDock Tools. Define an appropriate grid box centered on the known ligand. Create a conf.txt file and run AutoDock Vina. Finally, analyze the generated poses and compare the predicted binding affinity to any experimental data you can find for this ligand-protein complex. Document your steps and observations, including screenshots of key preparation steps and the final docked poses.
Watch the full lesson — free
This topic is part of Molecular Docking with AutoDock Vina, a complete AI-narrated video course. Press play once and watch the entire lecture like a movie.
Start the course free →