Alexandre Bonvin bio photo

Computational Structural Biology group focusing on dissecting, understanding and predicting biomolecular interactions at the molecular level.

Email Twitter Github Youtube Subscribe


Supported by:




HADDOCK2.4 manual - Pseudo Contact Shifts

Pseudo contact shifts (PCS) can provide useful information on both the distances and the orientation of the molecules to be docked. They can be used directly as restraints in HADDOCK using the XPCS energy term in CNS. Note that for this the standard CNS version needs to be recompiled with the files provided in the cns1.3 directory in the HADDOCK distribution.

For each set of PCS restraints used, a tensor must be defined. Its axial (D) and rhombic (R) components must be defined.

Note that the proper units for use in HADDOCK should be: 10-28 / (12*pi) m3 (which gives a scaling factor of 265.26 compared to values expressed in 10-32 m3)

For pseudo contact shifts, R and D this can be obtained using for example the Numbat software (Schmitz C, Stanton-Cook MJ, Su XC, Otting G, Huber T (2008) Numbat: an interactive software tool for fitting delta chi-tensors to molecular coordinates using pseudocontact shifts. J Biomol NMR 41:179-189). You can request by email the source code for Linux (recommended) or a Windows executable by writing to numbat.pcs@gmail.com .

The position of the tensor also needs to be defined by specifying distance restraints with respect to the protein to which it is attached. These restraints should be defined by the TENSOR_FILE in the run.param file used to setup a run. The entry in run.param looks like:

  TENSOR_FILE=tensor_basic.tbl

and this tensor restraint file contains distance restraints. For example, if the tensor is attached to a metal ion in the case of a metallo-protein (as in the protein-protein-pcs example in HADDOCK2.4) it looks like:

assi (resid  999 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  998 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  997 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  996 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  995 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  994 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  993 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  992 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  991 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0
assi (resid  990 and resn XAN and name  OO ) (resid 190 and segid A ) 0.0 0.0 0.0

The distances are thus all set to 0 to keep the tensor on top of the metal ion. 10 such distances are defined, one for each PCS set and its associated tensor (numbered from 990 to 999).

The proper format for PCS restraints is the following:

assign ( resid 999  and name OO )
       ( resid 999  and name Z )
       ( resid 999  and name X )
       ( resid 999  and name Y )
       ( resid 7 and name N and segid A ) -0.119906 0.15

where the last two numbers are the PCS value and its associated error.

The 2.4 version of HADDOCK supports up to 10 different PCS restraints sets. Each can have a separate tensor. The tensor residue number should be in the range 999-990. For each class you can specify the first and last stage at which the various RDCs will be used:

  • 0: rigid body EM (it0)
  • 1: semi-flexible simulated annealing (SA) (it1)
  • 2: explicit solvent refinement (water)

The PCS restraints should be defined by the PCSX_FILE parameter in the run.param file used to setup a run, with X being the number of the restraint set. An example entry in run.param would be:

  PCS1_FILE=./DATA/Dy_eps_HOT_0.15_.tbl
  PCS2_FILE=./DATA/Er_eps_HOT_0.15_.tbl
  PCS3_FILE=./DATA/Tb_eps_HOT_0.15_.tbl
  TENSOR_FILE=./DATA/tensor_basic.tbl

Example taken from the protein-protein-pcs example provided with HADDOCK2.4.