IMP Tutorial  for IMP version 2.2.1
Docking of PCSK9 with idock

Introduction

This example demonstrates the integrative docking protocol by application to the PCSK9 antigen–J16 Fab antibody complex. It uses a command line tool called idock.py. For full help on this tool, run from a command line:

idock.py -h

Setup

First, obtain the input files used in this example and put them in the current directory, by typing:

cp <imp_example_path>/saxs/pcsk9/* .

(On a Windows machine, use 'copy' rather than 'cp'.) Here, <imp_example_path> is the directory containing the IMP example files. The full path to the files can be determined by running in a Python interpreter 'import IMP.saxs; print IMP.saxs.get_example_path('pcsk9')'.

The input files include the Protein Data Bank (PDB) structures of the isolated J16 Fab antibody and PCSK9 antigen, antibody_cut.pdb and 2p4e.pdb, respectively. The protocol can also include one or more additional sources of experimental information. In this case, a set of cross links between the antibody and antigen are provided in the file cxms.txt. The cross links file has a simple format; each line lists the residue name and chain ID of the antibody and the antigen, followed by the minimum and maximum CA-CA distance in angstroms.

Calculation

All steps of the integrative docking procedure are carried out by simply running:

idock.py --complex_type AA --cxms=cxms.txt --precision=1 antibody_cut.pdb 2p4e.pdb

The application makes use of the PatchDock and FiberDock programs for docking and refinement, which must be obtained separately from the PatchDock website, and the reduce program for adding hydrogens to PDB files. It expects to find 'reduce' in the system path. It must be told explicitly where to find PatchDock and FiberDock (either by giving the –patch_dock and –fiber_dock options to idock.py, or by setting the PATCH_DOCK_HOME and FIBER_DOCK_HOME environment variables).

The –precision argument above controls how many candidate dockings are sampled by PatchDock, and can range from 1 (normal) to 3 (high precision). Higher values result in longer run times. The –complex_type argument is passed to PatchDock and FiberDock, and identifies the type of complex (AA means antibody-antigen in this case). It alters how the docking is performed (for more information, see the PatchDock help).

Results

Once the docking procedure has finished, the primary output file generated is results_cxms.txt, the first few lines of which look similar to:

     # |  Score  | filt| ZScore |cxms    | Zscore | Energy | Zscore | Transformation
     1 |  -1.177 |  +  | -1.695 |  0.016 | -0.454 |  7.310 | -0.723 |   -1.5638 -1.2242 1.9240 2.4535 -2.4522 1.4710
     2 |   0.276 |  +  | 0.398 |  0.004 |  0.999 |  0.000 | -0.723 |   -0.9194 -1.0958 -2.0672 -2.7843 -1.1249 -3.0764

Each line corresponds to one model; the models are ranked by total score, best first. The scores and z-score pairs for each experimental method (just cxms in this case) are also shown (only docking solutions that were not filtered out by any of the data sources – i.e. they scored well against every source – are included in this file). The last column is a transformation (3 rotation angles and a translation vector) that transforms the antibody relative to the antigen (the antigen is not transformed).