Difference between revisions of "2024 DOCK tutorial 2 with PDBID 1NDV"

From Rizzo_Lab
Jump to: navigation, search
(Introduction)
(Gridbox)
Line 101: Line 101:
 
=== sphere selector ===
 
=== sphere selector ===
 
== Gridbox ==
 
== Gridbox ==
=== showbox ===
+
While it is possible to calculate interactions between the ligand and the entire protein surface, this is computationally expensive. In this tutorial we will define a box of interest in the active site of the protein, and DOCK will generate a grid to calculate interaction energies within this box.
=== grid ===
+
=== Generating the Box ===
 +
To generate the box we use a DOCK program called showbox. This program is accessed through the command line. First navigate to your 003.gridbox directory and create a new file called showbox.in. Place the following commands into the new file:
 +
    Y
 +
    8.0
 +
    ../002.surface_spheres/selected_spheres.sph
 +
    1
 +
    1ndv.box.pdb
 +
 
 +
Remember to change "1ndv" to your protein name. This file may need to be modified for your system if you need a larger box. This file currently draws a box of 8.0 angstroms around the selected spheres. To change this, simply change the "8.0" in the input file.
 +
 
 +
To run this file type:
 +
    showbox < showbox.in
 +
If showbox ran successfully you will now have a file named 1ndv.box.pdb in the directory.
 +
 
 +
=== Generating the Grid===
 +
 
 
== Energy Minimization for Ligand ==
 
== Energy Minimization for Ligand ==
  

Revision as of 15:31, 16 March 2024

Introduction

DOCK is a powerful computational tool used in drug design. "Docking" refers to the identification of low energy binding conformations of a small molecule, or ligand, in an active site. There are three methods in which DOCK performs docking: rigid, fixed anchor, and flexible. DOCK also has the ability to perform a virtual screen, which is analogous to a high throughput screen often performed in wet labs. This method allows for thousands of molecules to be docked and ranked so that the highest scoring molecules can move to the next stage of testing. A virtual screen can save both time and money when compared to a traditional high throughput screen. This tutorial will guide you through the process of performing each method of docking, a virtual screen, and associated analysis using PDB 1NDV on a HPC cluster (Seawulf). Replace 1NDV with your target protein in this tutorial.

Required Software and Skills

In order to complete this tutorial UCSF Chimera is required. Chimera is a powerful visualization tool that will aid in the preparation and analysis of the virtual screen. Please note that the newer version, Chimera X, can not be used for this tutorial as it lacks certain integrations with DOCK. While this tutorial will walk you through all steps necessary to use Chimera, it may be helpful to view this tutorial for a more in depth guide on the program.

It is also recommended that you have experience with Unix and vi. A Unix tutorial can be found on the site here and a vi tutorial can be found here.

Set-up

Before we begin, we need to set up our directories and download our file from the Protein Data Bank (PDB).

Creating Directories

Create a parent directory which will contain all subdirectors for each docking step. For this tutorial we will name the parent directory 1NDV_DOCK_VS. Now create the subdirectories using the following scheme.


1NDV Chart.png


Downloading PDB File

Visit the RCSB Protein Data Bank (PDB) website and use the search bar to search for your target protein. As stated before, we will be using 1NDV for this tutorial.


PDB SearchBar.png


The page for the protein will look like this. Click the "Download Files" dropdown menu, and download the "PDB Format".


1NDV-Download.png


With the PDB file saved, we are now ready to begin generating files for DOCK!

Preparing Structure Files

Protein

Analyzing the Protein Structure:

Open your protein.pdb in Chimera.

Assess the protein structure:

1. Are there any missing loops?

2. Are there any metal coordination atoms forming key interactions with the protein?

The .pdb for 1ndv from the RCSB Protein Data Bank is shown below.

1ndv raw pdb.jpg

Assessment:

1. No missing loops.

2. Zinc coordination with the protein and a water molecule. This zinc atom must be kept in the protein prep.


1ndv coordination.jpg

Ligand

We must first create a .pdb file with only the ligand. To do this we must select the ligand in Chimera. There are two ways to do this:

1.Select an atom on the ligand, and use the up arrow until the entire ligand is selected

-or-

2. Use the toolbar and click "Select" → "Structure" → "Ligand"

Once the ligand is selected use "Select" → "Invert" to select all other atoms present in the session.

Then use "Actions" → "Atoms/Bonds" → "Delete" to delete everything except the ligand.

It should look like this.

1ndv-ligand-only.png

You can now save only the ligand in a separate .pbd file. We will name this file 1ndv_ligand_only.pdb.

We will also save the ligand in this state as a .mol2 file. We will name this file 1ndv_no_hydrogens_no_charges.mol2

We now need to add hydrogens and charges like we did previously for the protein. It is important to be careful and to cross reference the structure of the ligand in the PDB to ensure that hydrogens are being added in the right places, and that the overall charge of the ligand remains the same. After adding hydrogens, your structure should look similar to this:

1ndv-ligand-only-H.png

We must now look at the structure in the PDB to ensure Chimera is adding the hydrogens in the proper place. To do this, scroll down on the PDB page of the protein until you see the table labeled "Small Molecules".

1ndv-small-mol.png

Now click on the "2D Diagram" of the ligand to bring up a larger picture.

1ndv-small-mol-big.png

It appears that Chimera has added an extra Hydrogen to N26. We need to remove that by simply selecting and deleting that hydrogen. The ligand should now look like this:

1ndv-ligand-with-H.png

We can now add charges to the ligand. When the dialogue box appears, be sure to set the overall charge to 0 for this ligand. With this step complete we can now save a .pdb file and .mol2 with the name 1ndv_ligand_hydrogens_charges.

Upload all four of the .mol2 files created during the protein and ligand preparation process to Seawulf for use with DOCK.

Surface Spheres

DMS file

sphgen

sphere selector

Gridbox

While it is possible to calculate interactions between the ligand and the entire protein surface, this is computationally expensive. In this tutorial we will define a box of interest in the active site of the protein, and DOCK will generate a grid to calculate interaction energies within this box.

Generating the Box

To generate the box we use a DOCK program called showbox. This program is accessed through the command line. First navigate to your 003.gridbox directory and create a new file called showbox.in. Place the following commands into the new file:

    Y
    8.0 
    ../002.surface_spheres/selected_spheres.sph
    1 
    1ndv.box.pdb

Remember to change "1ndv" to your protein name. This file may need to be modified for your system if you need a larger box. This file currently draws a box of 8.0 angstroms around the selected spheres. To change this, simply change the "8.0" in the input file.

To run this file type:

    showbox < showbox.in

If showbox ran successfully you will now have a file named 1ndv.box.pdb in the directory.

Generating the Grid

Energy Minimization for Ligand

Footprint Generation

Rigid Docking

Fixed Anchor Docking

Flex Docking

Virtual Screening

Slurm

Cartesian Minimization

reScore