21.12 Induced Fit (Flexibility) Docking Tutorials
[ Multiple Receptor | Explicit Group ]
21.12.1 Multiple Receptor Conformation Ensemble Docking Example
Here we will demonstrate multiple receptor conformation docking in ICM using Aldose Reductase as an example. Aldose Reductase has a flexible loop in the ligand binding pocket vicinity which enables a variety of inhibitors to bind therefore in order to identify these ligands via docking it is necessary to sample the conformations of this loop and dock to an ensemble of structures. This example demonstrates that a single conformation of the Aldose Reductase receptor cannot account for the binding of all known aldose reductase inhibitors. To account for the flexibility we model the flexible loop and then dock using ICM Multiple Receptor Conformation (4D) docking. You can read more about 4D docking and other methods for incorporating induced fit here.
| Step 1:
- Load PDB file 1pwm. You can search for a PDB file using the PDB Search tab.
- Delete the chlorine atom. This can be done by right clicking on the chlorine in the ICM Workspace and selecting delete.
- Move the ligand (fid) out of the receptor. You can do this by right clicking on the ligand (fid) in the ICM Workspace and select Move From Object.
| Step 2: Convert 1pwm into an ICM object. * Rename the small molecule (fid) object "ligand" and the receptor object "receptor". You can rename an object by right clicking on the name of the object in the ICM Workspace and select rename.|
| Step 3: Setup the docking project using the ligand to guide the binding site. Double click on the ligand in the ICM Workspace and then select the Define Site Around Selected Ligand button (you should see green crosses surrounding the ligand on the receptor).
| Step 4: The probe represents the docking starting position. Generally you can keep it in the default position. If you want to move it use the middle mouse button and hold the SHIFT button for global rotation.
| Step 5: The purple box represents the region in which maps will be generated. To change the box size you can use the left mouse button with the cursor at any corner of the purple box to change it. Make sure the purple box is as close to the pocket as possible.
| Step 6:
In the next part of this exercise we will try to determine if other inhibitors to this receptor could be correctly docked using this one crystal structure or if more than one representation of the receptor is needed. We have already setup the docking project and made all the required maps (Steps 1 to 3) so now let us dock a set of inhibitors to the receptor. We will dock a database called ALDR_ligs.sdf which should be present in your ICM distribution if you cannot find this example sdf file then please E mail firstname.lastname@example.org and we will E mail it to you.
| Step 7: Dock the chemical table. |
| Step 8:
Browse the docked compounds and notice the positioning of the loop (res# 298:302) in the vicinity of the pocket (red). The docking demonstrates that a single conformation of the Aldose Reductase receptor cannot account for the binding of all known aldose reductase inhibitors.
| Step 9: |
- Expand the sequence display for the ALDR_rec object in the ICM Workspace and select the loop region "CALLSC" (it is located towards the C-terminal). See screenshot below. You can select the region by clicking and dragging over the sequence in the ICM Workspace.
| Step 10: Model the loop using MolMechanics/Loop/Sampling-Modeling |
- If the option is greyed out and you cannot click on it this means the selection you have made is not in the current object. You can fix this by right clicking on the object in the ICM Workspace (in this example the object is called ALDR_rec) and select "Set to Current".
- Choose the interactive tab.
- Choose the option Loop Dbase Search.
- Choose the option Store in object.
| Step 11:
- Once the loop sampling has finished you will see the embedded stack in the receptor in the ICM workspace (see image below). The stack is ranked by energy and you can display a conformation by clicking on it in the ICM workspace
| Step 12:
- Delete all but the top 4 conformations in the stack. This can be done by clicking on conformation #5 and this will highlight the row in blue and then whilst holding the Shift key click on the the last member of the stack (see image above). Next right click on the selection and select Delete.
| Step 13: Build the maps for each loop conformation using Choose the menu option Docking/Flexible Receptor/Setup 4D Grid.
| Step 14: Dock the ALDR_ligs table again to the multiple receptors using Docking/Dock Chemical Table. |
| Step 15: |
- Browse the hitlist by clicking on the hitlist table and watch the receptor conformation switching along wih the ligand and note the additional column rec.conformation in the list. This column indicates which receptor in the stack the ligand is docked to.
- Load PDB entry 1iei and compare the pose of the first ligand in the hit list and the corresponding receptor conformation to the experimental data.
- You will see that the flexibility in the loop region is accounted for by multiple receptor conformation docking and therefore inhibitors which scored poorly in Step 4 now score well and the corresponding receptor loop conformation for a ligand is identified.
21.12.2 Explicit Group Docking
This example demonstrates how to keep certain residues as explicit during grid docking. Full side-chains or hydroxyls and water can be treated as flexible. In this example we look at some crystal structures of the tyrosine kinase domain of Fibroblast Growth Factor where we observe that inhibitors cannot be docked into the apo-form without a Lysine (K514) residue in the pocket being treated as flexible.
| Step 1: Download and open the prepared ICM file. The file can be downloaded here ftp://ftp.molsoft.com/pub/workshop/ (1. cut and paste the address into a web browser. 2. Right click on hybrid_docking_kinase.icb and choose "Save Link As". 3. Open the .icb in ICM File/Open). Notice that the inhibitor from PDB 2fgi cannot dock to the apo structure (PDB 1fgk) because of the position of K514).
| Step 2: Setup the docking project using the ligand from PDB 2fgi to define the docking site.
| Step 3: In this example you do not need to change the position of the probe or size of the box. Press GO.
| Step 4: Dock the chemical table called "chem". This table contains two FGFR1 inhibitors from PDB 4RWL and 2FGI.
| Step 5: Browse the hitlist and you will see that the inhibitors cannot dock correctly and the scores are poor because of the position of K514.
| Step 6: Select K514 in the docking project receptor object.
| Step 7: Go to Docking/Flexible Receptor/Explicit Groups and choose Full Side-Chain.
| Step 8: Remake the receptor maps (Docking/(Re) Make Receptor Maps).
| Step 9: Dock the chemical table (chem) again but this time uncheck the box "Dock Immediately".
| Step 10: Prepare the table for docking by removing the check from the "docked" box in the chem table.
| Step 11: Change the thoroughness (effort) to 5 or more and click dock. This option increases the length of the simulation to account for the flexibile residue.
| Step 12: Once the docking has finished observe the improved scores and poses and see that the Lysine residue has moved to accomodate the ligand.