Week 11 - Ligand Binding, Docking and Drug Discovery¶
How to find a ligand or drug¶
There are many databases which contain coordinate files for ligands, market drugs, and drug-like molecules that potentially target enzymes or receptors involved in disease.
It is possible to search similar molecules for a drug or ligand of interest via these databases.
Ligand Expo¶
Ligand Expo provides chemical and structural information about small molecules within the structure entries of the Protein Data Bank.
Example 1: Find the ligand for c-Abl kinase (PDB ID: 1IEP) using PDB and Ligand Expo¶
- Go to the PDB website and search
1IEP. - See that one of its ligands is imatinib, an anti-cancer drug.
- As you notice, it has a 3-letter code:
STI. - Go to Ligand Expo:
- http://ligand-expo.rcsb.org
- From the menu bar, click the Search icon.

- Search for the ligand using its 3-letter code (
STI) by typing it into the search box and clicking Go.

- In this way you can learn molecular weight, charge, molecular shape, and formula of the ligand.
- You can also download the SDF (structure data file) or PDB file of this molecule to further analyze it.
- Click Downloads from the sub-menu bar.

- Choose either SDF or PDB format to download the structure of this ligand.
- You can analyze this structure in Jmol, UCSF Chimera, or SPDBV.
How to find a target protein for a chosen drug or ligand¶
Ligand Expo has another search tool which makes finding proteins that a ligand/drug targets possible.
Example 2: Find the PDB structures that STI (imatinib) binds¶
- Go to http://ligand-expo.rcsb.org
- Go to Search option from the menu bar.
- At the very bottom of the page, there is another search option:
- SEARCH FOR INSTANCES OF CHEMICAL COMPONENTS BY 3-LETTER ID CODE

- Click Go.
- You will see the list of PDB structures that contain imatinib.
How to find alternative drugs for a target protein¶
BindingDB¶
BindingDB is a public, web-accessible database of measured binding affinities, focusing chiefly on interactions of drug-target proteins with small, drug-like molecules.
Example 3: Find alternative drugs that may bind to Abl kinase (besides imatinib)¶
- Go to http://www.bindingdb.org
- Choose
ABLfrom the Quick Search → Choose target option and click Go.

- This page displays binding characteristics/parameters like ΔG°, Ki, Kd, etc.
- Select PD153035 (Novartis Pharmaceuticals®).
- You will see links beside the drug name (B.MOAD, DrugBank, Google Scholar, PDB, etc.).
- To learn more, go to DrugBank via the link provided.
- After clicking the link, a list of DrugBank IDs will be displayed. Click the uppermost ID:
DB02567. - You can find:
- chemical formula
- 3-letter code
- molecular weight
- water solubility
- SDF and PDB files of the ligand
At the same time, you can find information about the target protein (sequence, molecular weight, function, cellular processes, etc.).
ZINC Database¶
Example 4: Find similar drugs of zoledronic acid (ZOL) using ZINC¶
Another anti-cancer drug is zoledronic acid (ZOL), which targets farnesyl pyrophosphate synthetase.
- Find the structure of
ZOLfrom Ligand Expo. - Go to http://zinc.docking.org/
- Click Search and Browse from the link list on the left.
- You will see a box to draw a molecule.
- Draw zoledronic acid using the atoms and bonds around the drawing box.

- To draw a bond or atom, click on it first, then click on the drawing box.
- Set similarity to 80% from the similarity option on the left.
- Run Query.
- A list of 80% similar drugs to
ZOLis displayed.
You can, for example, use one of these drugs in cell culture or animal studies and move onto phase I, II, III studies. If it works, you may design a more efficient drug to treat cancer.
Docking¶
Docking predicts the orientation of a ligand and its binding site on the target molecule by calculating binding affinities for each binding state.
There are several docking software options:
- Some treat the ligand as movable and the target as stable
- Some treat the ligand as stable and the target as movable
- Very few consider both ligand and target as movable
AutoDock Vina¶
AutoDock Vina is an open-source program for drug discovery, molecular docking and virtual screening. It offers multi-core capability, high performance, and enhanced accuracy and ease of use.
AutoDock Vina has been designed and implemented by Dr. Oleg Trott in the Molecular Graphics Lab at The Scripps Research Institute.
Example 5: Predict the binding site and orientation of STI (imatinib) on Abl kinase (1IEP) using AutoDock Vina¶
For docking:
- First step is to create PDBQT files, which contain positions of atoms and their charges.
- Then, a region of the protein is specified as the site where the ligand may bind.
- Then the program is run to calculate binding affinities of possible ligand conformations.
Notes:
- The first model is always the one with the best binding affinity, so it is the most possible binding mode.
- Binding modes can be visualized with AutoDock Tools.
Please watch the tutorial for AutoDock Vina:
Links to learn more and download AutoDock Vina and AutoDock Tools:
- http://vina.scripps.edu/
- http://vina.scripps.edu/download.html
- http://mgltools.scripps.edu/
- http://mgltools.scripps.edu/downloads
AutoDock Vina Steps¶
- Open AutoDock.
- Open protein PDB.
- Click Edit → Hydrogens → Add → select polar only.
- Click Grid → Macromolecule → Choose → select the protein you want to dock on.
- Save protein with hydrogens in PDBQT format.
- Select the binding pocket using Grid → Gridbox.
- Save coordinates of the box.
- Open ligand PDB (Ligand → Input → Open →
ligand.pdb). - Pick rotatable bonds using Ligand → Torsion Tree → Choose torsion.
- Save ligand as PDBQT.
- Now close all windows and open a command prompt and go to the folder containing your PDBQT files.
- Use command prompt to run Vina, or prepare a configuration file and let Vina read preferences.
- Prepare a configuration file: create a text file
conf.txtunder the same folder. - Go to the Vina folder containing
conf.txtand PDBQT files in a terminal. - Type the command:
vina --config conf.txt –log log.txt