Difference between revisions of "2020 DOCK tutorial 4 with PDBID 6UZW"
Stonybrook (talk | contribs) (→Generating the box and grid) |
Stonybrook (talk | contribs) (→Selecting Spheres) |
||
| Line 133: | Line 133: | ||
===Selecting Spheres=== | ===Selecting Spheres=== | ||
| − | Here we will be selecting the spheres which | + | Here we will be selecting the spheres which define the binding pocket of the ligand because we are trying to direct the ligand towards that binding site rather than all over the receptor. To select the spheres type the following command. |
sphere_selector 6UZW_rec.sph ../01_dockprep/6UZW_lig_withH.mol2 10.0 | sphere_selector 6UZW_rec.sph ../01_dockprep/6UZW_lig_withH.mol2 10.0 | ||
| − | This command will select all of the spheres within 10.0 angstroms of the ligand and output them to selected_spheres.sph. Visualize the selected spheres using Chimera to make sure the correct spheres are selected. Notice that, spheres around the ligand binding site are kept and all the other spheres are deleted in the image below. | + | This command will select all of the spheres within 10.0 angstroms of the ligand and output them to selected_spheres.sph. Visualize the selected spheres using Chimera to make sure the correct spheres are selected. Notice that, spheres around the ligand-binding site are kept and all the other spheres are deleted in the image below. |
[[File:selected_spheres_6UZW.jpeg|thumb|center|500px|Figure . Selected spheres files of 6UZW]] | [[File:selected_spheres_6UZW.jpeg|thumb|center|500px|Figure . Selected spheres files of 6UZW]] | ||
| Line 154: | Line 154: | ||
Step 4. Add charge. | Step 4. Add charge. | ||
*Add charges Tools > Structure Editing > Add charge Typically use AM1-BCC | *Add charges Tools > Structure Editing > Add charge Typically use AM1-BCC | ||
| − | *You need to be careful about selecting charge because this will affect predictions for interactions (i.e. Coulomb's Law). Read paper corresponding to your PDB | + | *You need to be careful about selecting charge because this will affect predictions for interactions (i.e. Coulomb's Law). Read paper corresponding to your PDB structure to determine conditions of purification and proper protonation. In the case of 6UZW our paper of reference is: https://doi.org/10.1038/s41467-020-14321-0 |
| − | Pose | + | Pose Reproduction |
Flexible, Fixed-anchor, rigid | Flexible, Fixed-anchor, rigid | ||
Revision as of 20:27, 15 February 2020
This tutorial is created by AMS 536 Spring 2020 Group 4. Group 4 members include Steven Pak, Caitlyn Cardetti, MiaoMiao He and Chuanzhou (Joey) Zhu.
How to for selecting protein: Pick a protein with few substrates and rotational bonds?
The PDB code chosen is 6UZW which is for the crystal structure of GLUN1/GLUN2A ligand-binding domain in complex with glycine and UBP791.
Dock Prep Program used: Chimera
Contents
I. Introduction
DOCK
DOCK 6.9 is a molecular modeling program that identifies interactions between chemical compounds and the receptors. DOCK has many features that can help with drug discovery( ex: Virtual Screening and Denovo). This tutorial will show us how to use Virtual Screening.
6UZW
The tutorial will be based on the PDB file 6UZE downloaded from the PDB Database. 6UZW is the crystal structure for GLUN1/GLUN2A complexed with (2S,3R)-1-[7-(2-carboxyethyl)phenanthrene-2-carbonyl]piperazine-2,3-dicarboxylic acid and glycine.
Organization of Directories
We set up the files in our project space as such. It will be helpful to have these folders ready ahead of time.
00_files
01_dockprep
02_surface_spheres
03_gridbox
04_dock
05_footprint
06_virtual_screen
07_virtual_screen_mpi
08_cartesianmin
09_rescore
II. Preparation of the ligand and receptor (preparation before the usage of DOCK6.9
- Download 6UZW to local computer off of the PDB database - Save the .pdb file into your 00_files directory.
Checking the Structure and Preparing the Complex without Hydrogens
- Open Chimera, open the downloaded .pdb file, and check the structure for missing residues, gaps, heme groups, missing loops, and size. We recommend using proteins with no heme groups and proteins of a relatively smaller size. - Save the .pdb file (the ligand and the receptor) as a mol2 file (6UZW_complex_noH.mol2) - When preparing both the ligand the receptor, you can open this complex mol2 file and delete the component that you won't need (eg. for preparing the receptor, open the complex, delete the ligand, and save).
This complex (receptor and ligand) has been prepared without Hydrogens.
Preparation of Ligand without Hydrogens
- Open 6UZW_complex_noH.mol2 through Chimera again - Isolate the ligand by deleting the receptor. Select on the receptor by clicking on (Select --> Residue --> QM4 --> Select --> Inverted (selected models). The receptor should be surrounded by a green film. - Delete it (Actions -> Atoms/Bonds -> Delete). - The ligand should be by itself without any amino acids from the receptors and any water molecules. - Save the isolated ligand as a mol2 file (File -> Save mol2 -> 6UZW_lig_noH_mol2) into 01_dockprep
This ligand has been prepared without Hydrogens.
Preparation of Receptor without Hydrogens
- Open 6UZW_complex_noH.mol2 through Chimera again - Isolate the ligand by deleting the receptor. Select on the receptor by clicking on (Select --> Residue --> Standard Amino acids --> Select --> Inverted (selected models). The ligand and water molecules should be surrounded by a green film. - Delete it (Actions -> Atoms/Bonds -> Delete). - The the receptor should be by itself without the ligand - Save the isolated receptor as a mol2 file (File -> Save mol2 -> 6UZW_rec_noH_mol2) into 01_dockprep
This ligand has been prepared without Hydrogens.
Preparation of ligand with Hydrogens
-Open the 6UZW_lig_noH.mol2 -Tools -> Surface Binding Analysis
Tools -> Structure Editing -> Add H (To add Hydrogen atoms)
Tools -> Structure Editing -> Add Charge (To add the charge use the latest AMBER force field available for standard residues. Here we used AMBER ff14SB)
Do not click or edit anything while going through this process.
Save as a mol2 file. (6UZW_lig_withH.mol2) into 01_dockprep
- If you follow the step below all the above stated steps will automatically appear one after the other to prepare the receptor.
Tools -> Structure/Binding Analysis -> DockPrep
Preparation of receptor with Hydrogens
-Open the 6UZW_rec_noH.mol2 -Tools -> Surface Binding Analysis
Tools -> Structure Editing -> Add H (To add Hydrogen atoms)
Tools -> Structure Editing -> Add Charge (To add the charge use the latest AMBER force field available for standard residues. Here we used AMBER ff14SB)
Do not click or edit anything while going through this process.
Save as a mol2 file. (6UZW_rec_withH.mol2) into 01_dockprep
- If you follow the step below all the above stated steps will automatically appear one after the other to prepare the receptor.
Tools -> Structure/Binding Analysis -> DockPrep
III. Generating receptor surface and spheres
Preparation of DMS file
- Open 6UZW_rec_noH.mol2 using chimera. - Action -> Surface -> Show - Tools -> Structure Editing -> Write DMS - Save the 6UZW_rec_noH.dms into 02_surface_spheres folder
Reopen the file and make sure the surface was generated.
Transfer all the folders created so far to seawulf cluster to be used in DOCK using the Secure Copy command on unix from the local terminal that you have saved your files. NOT from the seawulf terminal
scp -r [FROM] [TO] scp -r PATH/folder/ PATH/seawulf_folder
Generating spheres
- Go to 02_surface_spheres folder - Create a new input file to create spheres by typing vim INSPH and type the following lines inside the file.
6UZW_rec_noH.dms R X 0.0 4.0 1.4 6UZW_rec.sph
The first line 6UZW_rec_noH.dms specifies the input file. R indicates that spheres generated will be outside of the receptor surface. X specifies all the points will be used. 0.0 is the distance in angstroms and it will avoid steric clashes. 4.0 is the maximum surface radius of the spheres and 1.4 is the minimum radius in angstroms.The last line 6UZW_spheres.sph creates a sphere file that contains clustered spheres.
Once the INSPH file is ready, type the following command to generate the spheres.
sphgen -i INSPH -o OUTSPH
Once sphgen command is successful, 6UZW_rec.sph file will be created. Open it up using Chimera along with 6UZW_rec_noH.mol2 file. You should get a similar output like the image below.
Selecting Spheres
Here we will be selecting the spheres which define the binding pocket of the ligand because we are trying to direct the ligand towards that binding site rather than all over the receptor. To select the spheres type the following command.
sphere_selector 6UZW_rec.sph ../01_dockprep/6UZW_lig_withH.mol2 10.0
This command will select all of the spheres within 10.0 angstroms of the ligand and output them to selected_spheres.sph. Visualize the selected spheres using Chimera to make sure the correct spheres are selected. Notice that, spheres around the ligand-binding site are kept and all the other spheres are deleted in the image below.
Step 1. Separate protein and ligand
Step 2. Delete H2O molecules in Chimera to reduce computing time.
- To delete: Select > Residues > H O H - they should be highlighted Click Actions > Atoms/Bonds > delete H2O molecules should now be deleted.
Step 3. Add H's.
- PDB codes do not contain H's due to low electron density.
- To add H's Tools > Structure Editing > AddH
- H's MUST be added before you add charge or they will not be accounted for and will affect predictions for interactions.
Step 4. Add charge.
- Add charges Tools > Structure Editing > Add charge Typically use AM1-BCC
- You need to be careful about selecting charge because this will affect predictions for interactions (i.e. Coulomb's Law). Read paper corresponding to your PDB structure to determine conditions of purification and proper protonation. In the case of 6UZW our paper of reference is: https://doi.org/10.1038/s41467-020-14321-0
Pose Reproduction Flexible, Fixed-anchor, rigid
Generating the box and grid
For this part, we will be generating a box and grid for docking the ligand. A space-defined box will be searching efficiently, and a grid-based energy system is beneficial for a rapid and thoroughly process.
Step 1. Create a Box.
- Go to 03_gridbox folder - Create a new input file to create a box by typing vi showbox.in and type the following lines inside the file.
Y 8.0 ../02_surface_spheres/selected_spheres.sph 1 6UZW.box.pdb
Those are basically parameters/answers for questions specified in the program. The first line Y represents "Yes, we want to create a box around the sphere". 8.0 indicates that the box will have a length of 8 Angstroms, in all six directions. The third line specifies the sphere file should be used. 1 is the cluster number for creating the box. The last line 6UZW_box.pdb creates a file that contains the box.
Once the showbox.in file is ready, type the following command to create the box.
showbox < showbox.in
Once this command is successful, 6UZW.box.pdb file will be created. Open it up using Chimera and you should be able to visualize the structure like the image below.
Step 2. Generating the Grid
-In the same directory, create a new file by typing vi grid.in and type the following lines inside the file.
compute_grids yes grid_spacing 0.4 output_molecule no contact_score no energy_score yes energy_cutoff_distance 9999 atom_model a attractive_exponent 6 repulsive_exponent 9 distance_dielectric yes dielectric_factor 4 bump_filter yes bump_overlap 0.75 receptor_file ../01_dockprep/6UZW_rec_prepped.mol2 box_file 6UZW.box.pdb vdw_definition_file /gpfs/projects/AMS536/zzz.programs/dock6(check again)/parameters/vdw_AMBER_parm99.defn score_grid_prefix grid
Once the file is ready, type the following command to generate the grid. More information can be found in http://dock.compbio.ucsf.edu/DOCK_6/tutorials/grid_generation/generating_grid.html
grid -i grid.in -o grid.out
The grid calculation might takes a while to complete. But once this command is successful, four more files will be in 03_gridbox directory:
6uzw.box.pdb grid.in gridinfo.out grid.bmp grid.nrg showbox.in
Open it up using Chimera and make sure everything looks reasonable. If all is good, you should be able to visualize the structure like the image below.
