Lesson · 40 min · Free
Molecular Docking Lab
Molecular Docking Lab body { font-family: sans-serif; line-height: 1.6; margin: 20px; } h1 { color: #2c3e50; } h2 { color: #34495e; border-bottom: 2px solid #ccc; padding-bottom: 5px; margin-top: 30px; } p { margin-botto
Molecular Docking Lab
Welcome to the "Molecular Docking Lab" lesson, part of our "Molecular Docking with AutoDock Vina" course. In this practical session, we will apply the theoretical knowledge gained in previous modules to perform a complete molecular docking experiment using AutoDock Vina. Our goal is to predict the binding affinity and pose of a small molecule ligand to a target protein, a critical step in drug discovery and development. Molecular docking is a computational method that predicts the preferred orientation of one molecule to a second when bound to each other to form a stable complex. It plays a pivotal role in lead optimization, virtual screening, and understanding molecular recognition. For pharmacy and biotech students, mastering this technique provides a powerful tool for rational drug design. Before we begin, ensure you have AutoDock Vina and MGLTools (specifically AutoDockTools, ADT) installed on your system. We will use ADT for preparing our protein and ligand files into the necessary PDBQT format, and Vina for performing the docking simulation itself. A text editor will also be useful for creating and modifying configuration files.
Setting Up Your Docking Experiment
The first step in any docking experiment is preparing your input files. This involves obtaining the 3D structures of your target protein and ligand, and then converting them into the PDBQT format. PDBQT files contain atomic coordinates, atom types, and partial charges, which are essential for Vina's calculations. For our exercise, let's assume you have downloaded a protein structure (e.g., from the PDB database) and a ligand structure (e.g., from PubChem or ChemSpider).
Protein Preparation with AutoDockTools (ADT)
Open ADT. Load your protein (File > Read Molecule...). Remove water molecules (Edit > Delete > Water). Add polar hydrogens (Edit > Hydrogens > Add > Polar Only). Compute Gasteiger charges (Edit > Charges > Compute Gasteiger). Finally, save the protein in PDBQT format (File > Save > Write PDBQT). This process ensures the protein is correctly prepared for docking, with appropriate charges and hydrogen atoms for interaction calculations.
Ligand Preparation with AutoDockTools (ADT)
Similarly, load your ligand into ADT (Ligand > Input > Open...). Detect root and torsions (Ligand > Torsion Tree > Detect Root and Torsions). This step identifies rotatable bonds in the ligand, allowing Vina to explore different conformations. Set the number of active torsions if necessary. Compute Gasteiger charges (Ligand > Charges > Compute Gasteiger). Save the ligand as a PDBQT file (Ligand > Output > Save as PDBQT...).
Creating the Configuration File
The configuration file tells AutoDock Vina where to perform the docking (the search space or "grid box") and other parameters. It's a simple text file. Here's an example: receptor = protein.pdbqt ligand = ligand.pdbqt center_x = 10.0 center_y = 20.0 center_z = 30.0 size_x = 20.0 size_y = 20.0 size_z = 20.0 out = docking_results.pdbqt log = docking_log.txt cpu = 8 exhaustiveness = 10 The center_x , center_y , center_z parameters define the center of your search box in Ångstroms, and size_x , size_y , size_z define its dimensions. You can determine these coordinates by examining your protein's active site in a visualization tool like PyMOL or ADT. The out parameter specifies the output file for docked poses, and log for the docking log. cpu sets the number of processor cores to use, and exhaustiveness controls the thoroughness of the search (higher values mean a more extensive search, but take longer).
Running AutoDock Vina
Once your protein PDBQT, ligand PDBQT, and configuration files are ready, you can run Vina from your terminal or command prompt. Navigate to the directory containing your files and execute the following command: vina --config conf.txt Replace conf.txt with the actual name of your configuration file. Vina will then perform the docking and output the results to the specified .pdbqt and .txt files.
Analyzing Results
The docking_results.pdbqt file will contain multiple predicted poses of the ligand within the binding site, ordered by their estimated binding affinity (lowest energy is best). You can visualize these poses in ADT or PyMOL to examine the interactions between the ligand and protein. The docking_log.txt file provides detailed information about the docking process, including the binding affinities for each pose.
Key Takeaways:
Molecular docking predicts ligand-protein binding poses and affinities. AutoDockTools (ADT) is essential for preparing protein and ligand PDBQT files. The configuration file defines the docking parameters, including the search space. AutoDock Vina is executed via the command line with a configuration file. Results are analyzed by visualizing docked poses and interpreting binding affinities.
Practice Exercise: Docking an Analgesic to a Cyclooxygenase Enzyme
For this exercise, you will dock Ibuprofen (ligand) to Cyclooxygenase-2 (COX-2) enzyme (receptor). Download the structure of COX-2 (PDB ID: 6COX) from the RCSB PDB database. Download the 3D structure of Ibuprofen from PubChem (CID: 3672). Save it as an SDF or MOL2 file. Use AutoDockTools to prepare both 6COX and Ibuprofen into their respective PDBQT formats. Using a visualization tool (e.g., PyMOL or ADT), identify the active site of COX-2 (e.g., where a co-crystallized ligand might be, or based on known literature). Determine appropriate center_x , center_y , center_z , and size_x , size_y , size_z values to create a docking box that encompasses this active site. Create a configuration file (e.g., ibuprofen_cox2.txt ) with the prepared PDBQT files and the determined grid box parameters. Run AutoDock Vina using your configuration file. Analyze the output: What is the estimated binding affinity (in kcal/mol) of the top-ranked pose? Visualize the top pose of Ibuprofen in the COX-2 active site. Describe any key interactions (e.g., hydrogen bonds, hydrophobic interactions) you observe.
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 →