2021 DOCK tutorial 4 with PDBID 1EFY
DOCK is a member of a large suite of molecular docking packages. Notably, it was the first of its kind, and it was created by Irwin Kuntz’s group in the 1980s. The original authors of DOCK include Brian K. Shoichet, David A. Case, and Robert C.Rizzo (link respective group pages). Molecular docking plays a key role in computer aided drug design (CADD). Generally, molecular docking allows us to explore the interactions between small molecules, or ligands, and proteins by evaluating the shape and energetics of the ligand. Add more about the specific algorithm. We will be using the 6.9 release of DOCK, the latest as of November 2018, for this tutorial. New features of DOCK6.9 include DOCK_DN, a new ligand searching method for ligand de novo design and fragment library generation, which allows users to generate their own fragment libraries from mol2 files. For a comprehensive list of what’s new in DOCK6.9, see: http://dock.compbio.ucsf.edu/DOCK_6/new_in_6.9.txt (link).
UCSF Chimera is a standalone program that allows us to visualize and analyze macromolecules. General features of Chimera include automatic identification of atom, hydrogen addition and partial charge assignment, measurements of distances, angles, surface area, volume, and many more. For a complete description of everything Chimera can do, see https://www.rbvi.ucsf.edu/chimera/.
In this tutorial, we will be working with the crystal structure of a member of the Poly (ADP-Ribose) polymerase (PARP) family. Specifically, we will be working with the catalytic fragment of PARP along with a benzimidazole inhibitor. PARP family proteins are involved in DNA repair. It has been shown that PARP may allow cancer cells to repair their DNA after exposure to chemotherapy, which allows for resistance to develop. As such, it is of particular interest to develop novel inhibitors for PARP. Goal for this tutorial: The goal for this tutorial is to demonstrate an end to end example of performing a virtual screen
Before we begin our molecular docking and virtual screen, we should setup our directory structure:
mkdir INSERT DIRECTORY NAMES
Preparing the Ligand and Receptor
Visit https://www.rcsb.org/ and search 1EFY. Click on Download files > PDB Format This file contains the crystal structure of the catalytic fragment of Poly (ADP-ribose) Polymerase complexed with a benzimidazole inhibitor. We will be visualizing the structure in Chimera. Open Chimera and click File > Open Select the PDB file that you have just downloaded
In order to prepare our system for DOCK, we must perform some preprocessing steps. Broadly speaking, we will be creating two new mol2 files out of the original PDB file: one for the receptor and one for the ligand. Both of these new structures will be “clean,” that is, they will not contain any of the extra ions that were present during the experiment but not biologically relevant. We will first create the mol2 file for the receptor. Using the original PDB structure, we will remove everything but the receptor. There are many ways to do this, but the easiest way is the following: Click Select > Residue > all nonstandard The benzimidazole inhibitor and the two water molecules should now be highlighted like so
Now select Actions > Atoms/Bonds > delete
This leaves us with a “cleaned” version of the receptor which we can now save as a mol2 file. Click File > Save Mol2… and save the file as “1EFY_rec_woH.mol2” where “woH” indicates that this is the receptor without hydrogens. Next we must add hydrogens and charges so that … // add why it’s important Select Tools > Structure Editing > AddH A window should pop up like this:
Select “OK.” Now each residue in the receptor should be properly protonated. However, it is important that we validate the hydrogens that Chimera has added, as they may not be correct. Similarly, we must add partial charges to each atom: Select Tools > Structure Editing > Add Charge. Make sure that AM1-BCC is selected in the window that pops up and select “OK”: Now each atom of the receptor should have the appropriate partial charge. As with the protonation states, it is important to doublecheck Chimera’s work. This can be done by // This receptor, which has now been prepared for DOCK, can be saved as a mol2 file named “1EFY_rec_dockprep.mol2”
Just as we preprocessed the receptor, we will now prepare the ligand. Open the original PDB structure in Chimera Now, we will be deleting all of the constituents in the structure except for the benzimidazole inhibitor. Click Select > Residue > BZC. This should highlight the inhibitor. Click Select > Invert (all models). This should highlight everything except the inhibitor. Click Actions > Atoms/Bonds > delete Only the ligand should now be present:
We will save this mol2 file File > save as mol2 > 1EFY_ligand_noH.mol2” where noH indicates no hydrogens Just as before, we must now add the hydrogens and charge.
Add charge: All four of the files should now be scp'd to the 001 directory.