Difference between revisions of "2014 DOCK tutorial with HIV Protease"
(→IV. Generating Box and Grid)
|Line 249:||Line 249:|
command or input file
command or input file
Revision as of 15:40, 3 March 2014
For additional Rizzo Lab tutorials see DOCK Tutorials. Use this link Wiki Formatting as a reference for editing the wiki. This tutorial was developed collaboratively by the AMS 536 class of 2013, using DOCK v6.6.
- 1 I. Introduction
- 2 II. Preparing the Receptor and Ligand
- 3 Downloading the PDB Structure (1HVR)
- 4 Preparing the ligand and receptor in Chimera
- 5 III. Generating Receptor Surface and Spheres
- 6 IV. Generating Box and Grid
- 7 V. Docking a Single Molecule for Pose Reproduction
- 8 VI. Virtual Screening
- 9 VII. Running DOCK in Parallel on Seawulf
- 10 VIII. Frequently Encountered Problems
DOCK is a molecular docking program used in drug discovery. It was developed by Irwin D. Kuntz, Jr. and colleagues at UCSF (see UCSF DOCK). This program, given a protein binding site and a small molecule, tries to predict the correct binding mode of the small molecule in the binding site, and the associated binding energy. Small molecules with highly favorable binding energies could be new drug leads. This makes DOCK a valuable drug discovery tool. DOCK is typically used to screen massive libraries of millions of compounds against a protein to isolate potential drug leads. These leads are then further studied, and could eventually result in a new, marketable drug. DOCK works well as a screening procedure for generating leads, but is not currently as useful for optimization of those leads.
DOCK 6 uses an incremental construction algorithm called anchor and grow. It is described by a three-step process:
- Rigid portion of ligand (anchor) is docked by geometric methods.
- Non-rigid segments added in layers; energy minimized.
- The resulting configurations are 'pruned' and energy re-minimized, yielding the docked configurations.
HIV protease (HIVPR) is a protein involved in viral maturation during the life cycle of HIV. HIVPR is an approximately 22 kDa homodimer with 99 residues per chain. Inhibition of this protein has been shown to be an effective form of treatment of HIV. Currently-available HIVPR inhibitors generally take the form of a symmetric cyclic urea compound. For more information, see Lam et al.
While performing docking, it is convenient to adopt a standard directory structure / naming scheme, so that files are easy to find / identify. For this tutorial, we will use something similar to the following:
~username/AMS536/dock-tutorial/00.files/ /01.dockprep/ /02.surface-spheres/ /03.box-grid/ /04.dock/ /05.mini-virtual-screen/ /06.virtual-screen/
In addition, most of the important files that are derived from the original crystal structure will be given a prefix that is the same as the PDB code, '1HVR'. The following sections in this tutorial will adhere to this directory structure / naming scheme.
II. Preparing the Receptor and Ligand
Downloading the PDB Structure (1HVR)
Go to PDB (Protein Data Bank) website (http://www.rcsb.org/pdb/home/home.do) enter the protein ID (1HVR), search for the PDB file and download it as a text form.
Preparing the ligand and receptor in Chimera
Put the 1HVR PDB file in 00.file/folder. If you are in the 00.files/directory, then tap the command:
cp ~/Downloads/1HVR.pdb ./
When you are preparing you PDB files, you have to make some modifications on your original file. For example: we changed the atom name form "CSO" to "CYS" and deleted to lines "OD" and "HD". When you finish the modifications, save it as "1HVR.modified.pdb" in the 00.files/. And then we will create 4 files in 01.dockprep/ directory:
1HVR.dockprep.mol2 1HVR.ligand.mol2 1HVR.receptor.mol2 1HVR.receptor.noH.pdb
For the "1HVR.dockprep.mol2" file: open the 1HVR.modified.pdb in Chimera; delete the water molecules; delete the original hydrogen atoms; add the charge and add the hydrogen atoms manually. Or you can do all of the above by clicking Tools -> Structure Editing -> Dock Prep. Note when adding the charge to the ligand, you can choose AMBER ff99SB as the charge model and chose gasteiger as the charge method. In this 1HVR case, we set Net Charge to 0. (You may have to consider the chemistry of the ligand when assigning a charge state).
III. Generating Receptor Surface and Spheres
Generating the Receptor Surface
Check to make sure 02.surface-spheres directory exists under dock-tutorial. If not then make the following directory:
mkdir 02.surface-sphgen cd 02.surface-sphgen
The following steps will be carried out to generate the receptor surface using Chimera:
Open Chimera by simply typing chimera into the terminal window
| Go File -> Open and choose the PDB file of the protein containing no hydrogens (1HVR.receptor.noH.pdb) from 01.dockprep
| Further, Actions -> Surface -> Show
| Go Tools -> Structure Editing -> Write DMS in order to obtain a dms file, which we will need to place spheres
| In the new window save the surface as 1HVR.receptor.dms
We will be using SPHGEN to generate spheres: see the DOCK online owners manual for additional information:
The following steps will be used to place the spheres on the receptor surface:
1. Create a file called INSPH and fill it out as follows, then save it. This input file tells SPHGEN what to do, details of each line are below:
1HVR.receptor.dms R X 0.0 4.0 1.4 1HVR.receptor.sph
Input File Details:
1HVR.receptor.dms - surface file from the previous step R - tells SPHGEN to place spheres either outside of the surface (R) or inside the surface (L) X - tells SPHGEN the subset of surface points to be used (X=all points) 0.0 - prevents generation of large spheres with close surface contacts(defalut=0.0) 4.0 - maximum sphere radius in angstroms (default=4.0) 1.4 - minimum sphere radius in angstroms (default=radius of probe) 1HVR.receptor.sph - clustered spheres file that we want to generate
2. Run the sphgen program from the terminal:
sphgen -i INSPH -o OUTSPH
-i tells sphgen where the input file INSPH is INSPH tells sphgen what to do -o tells sphgen what to call the oputput file OUTSPH is the output file containing the sphere information
3. (optional) To look at the spheres generated, you need to put them into PDB format.
Run showsphere, by typing the follwoing into the terminal:
You will be prompted with the following questions:
Enter name of sphere cluster file: 1HVR.receptor.sph Enter cluster number to process (<0 = all): -1 Generate surfaces as well as pdb files (<N>/Y)? N Enter name for output file prefix: output_spheres Process cluster 0 (contains ALL spheres) (<N>/Y)? N
You can then open the receptor file in Chimera as well as the output_spheres.pdb file and should see many spheres placed all over the receptor surface.
4. Run the sphere_selector program from in the terminal to select spheres of interest (note: not all of the spheres in the previous image are in the active site so we want to eliminate them):
sphere_selector 1HVR.receptor.sph ../01.dockprep/1HVR.ligand.mol2 8.0
This program goes to select the spheres within a user-defined radius (8.0 here) of a target molecule from a previously obtained file: 1HVR.receptor.sph. In turn, a new file selected_spheres.sph will be generated.
5. Run showsphere to visualize the spheres:
When prompted on the command line, answering the questions as follows:
Enter name of sphere cluster file: selected_spheres.sph Enter cluster number to process (<0 = all): -1 Generate surfaces as well as pdb files (<N>/Y)? N Enter name for output file prefix: output_spheres_selected Process cluster 0 (contains ALL spheres) (<N>/Y)? N
- Launch Chimera, choose File -> Open, choose 1HVR.receptor.noH.pdb
- Go File -> Open, choosing output_spheres_selected.pdb
- Go Select -> Residue -> SPH
- Finally, Actions -> Atoms/Bonds -> sphere
The final image should look similar to the example below:
IV. Generating Box and Grid
1.) Make a new directory and name it: 03.box-grid/
2.) Make a new file in this directory and name it showbox.in
3.) This will automatically open the file showbox.in. Edit the file showbox.in as follows:
Y #Yes, generate a box 8.0 #Size of the box in Angstroms ../02.surfaces-spheres/selected_spheres.sph #Sphere.sph file 1 #Cluster number 1HVR.box.pdb #Name of the output file
4.) Save the file by the inputting the command:
5.) Run the command:
showbox < showbox.in
V. Docking a Single Molecule for Pose Reproduction
Jess Junjie Kai
VI. Virtual Screening
Virtual Screening Introduction
Virtual screening is a method used to predict most favorable ligand binding to a target protein within a ligand database. It also allows for comparison of both qualitative (e.g. position in binding site) and quantitative (e.g. grid scores, internal energy) data pertaining to the each screened ligand with the originally docked molecule.
To perform virtual screening, we use HIVPR.ligands.005.mol2, a mol2 file which contains 5 small molecules to be the virtual library. This is a reasonable computational cost for a quick search, so we can conduct it on own computer . After which, we may able to conduct virtual screening within a larger database HIVPR.ligands.100.mol2. Since the computational cost of virtual screening is much higher, it is better to run it on Seawulf.
Running Virtual Screen
mkdir 05.mini-virtual-screen cd 05.mini-virtual-screen cp ~wjallen/AMS536/multi-mol2/files/HIVPR.ligands.005.mol2 ./ cp ../04.dock/dock.in/ ./ mv dock.in mini-virtual-screen.in
VII. Running DOCK in Parallel on Seawulf
VIII. Frequently Encountered Problems
Write some text here..
command or input file