PRODIGY: A Contact-based Predictor of Binding Affinity in Protein-protein Complexes

Equipment

1. A computer with internet access

Software

1. A web browser (the PRODIGY server has been tested successfully on Chrome, Firefox and Safari)
   PRODIGY web server address: http://milou.science.uu.nl/services/PRODIGY
   
2. Software repositories for running a local version (not described in this protocol) under a Linux or MacOSX operating system:
   
1. PRODIGY repository (https://github.com/haddocking/binding_affinity)
   
2. freeSASA (http://freesasa.github.io)
   

Procedure

1. The software
   1. Technical description
      PRODIGY is made freely available to the scientific community either as standalone software (https://github.com/haddocking/binding_affinity), which can be used locally on a desktop computer, or more conveniently as an online web-server, for which the usage is explained in this protocol. The PRODIGY software consists of a collection of Python scripts, a few Perl scripts to handle the online submission and the open-source tool freeSASA (Mitternacht, 2016) used to calculate the solvent accessible surface area, using default NACCESS (Hubbard and Thornton, 1993) parameters for atomic radii (http://freesasa.github.io).
      
   2. Data requirement
      
   1. Input file – mandatory
      The main input required to perform binding affinity prediction is a text file containing the atomic coordinates describing the 3D structure of the protein-protein complex (or ensemble of complexes). They can either be experimental structures solved e.g., by X-ray crystallography or NMR spectroscopy, which can be obtained from the worldwide Protein Data Bank (wwPDB) (http://www.wwpdb.org) (Berman et al., 2003), or structures modeled through computational approaches, e.g., by homology modelling or docking approaches. The 3D coordinates should be provided to the server in PDB or mmCIF format.
      The input structure/structures can be provided in different ways:
      
      1. By uploading a PDB or mmCIF file.
         
      2. By providing a PDB code for automatic retrieval form the wwPDB (http://www.wwpdb.org).
         
      3. By uploading an archive file (.tar, .tgz, .zip, .bz2 or .tar.gz) containing multiple PDB/mmCIF files. This option allows the submission of a unique file when many structures have to be analyzed (e.g., models derived from docking simulations).
         
   2. Chains – mandatory
      It is necessary to specify chains identifiers for the molecules involved in the interaction. If one (or both) interacting molecule is made of multiple chains at the interface, they all have to be provided separated by comma.
      
   3. Temperature – optional
      The user can specify at which temperature to perform the calculation of the dissociation constant (K_d). If nothing is specified, PRODIGY will use 25 °C by default.
      
   4. Job name – optional
      If provided, the job name will be used to identify your run, otherwise a random label will be assigned.
      
   5. Email – optional
      If an email is provided, a link to the results will be sent when the job has completed.
      
      
2. How to use PRODIGY web-server
   1. Submitting a prediction
      Here we describe the process of submitting a prediction run to the PRODIGY web-server (http://milou.science.uu.nl/services/PRODIGY). As example, we will use the protein-protein complex between an antibody (FAB) and HIV-1 capsid protein p24, that is present in the Protein Data Back (PDB) with the access code ‘1E6J’.
      
   1. Open an internet browser and go to http://milou.science.uu.nl/services/PRODIGY.
      
   2. Fill in the PRODIGY input page (Figure 1):
      
      1. Insert the PDB code ‘1E6J’ into the ‘Structure’ box for automatic retrieval from the wwPDB.
         
      2. The complex 1E6J is made of 3 chains: P (corresponding to HIV-1 capsid protein p24) and L and H (corresponding to the FAB). Considering that in this protocol we want to investigate the binding affinity at the interface between the antibody (chains L + H) and the antigen (chain P), you will need to insert:
         Interactor 1 ID_chain(s): P
         Interactor 2 IC_chain(s): L, H
         
         
         Figure 1. Example view of an input page of the PRODIGY web-server. (http://milou.science.uu.nl/services/PRODIGY)
         
         
   3. Personalize your job by defining some (optional) parameters if needed:
      
      1. In the box related to the temperature, change the 25 °C default if needed (note that this only affects the calculation of the dissociation constant and not the binding affinity ΔG).
         
      2. Give a name to your prediction run. No spaces or special characters other then ‘-‘ or ‘_’ are allowed. For this example we will name our run ‘1E6J_prediction’.
         
      3. Add your email address to be notified when your job is done and receive the link to the results page.
         
   4. We are now ready to send the prediction to PRODIGY: Click on the Submit button at the bottom of the page.
      
   5. A prediction usually does not take more than a few minutes. After this time, you will be redirected to the result page. If an email has been provided (see the above step B1c.iii), you will be notified when the prediction is complete and receive a link to the results page. Please note that the results are only stored for 2 weeks.
      
   2. The result page
      The result page is organized in three sections, reporting different information:
      
      1. Binding affinity and K_d prediction
         The name identifiers of your complex, which contains the PDB code of the retrieved file (or the name of the input you upload) is reported, together with the predicted ΔG (in kcal mol^-1) and K_d (in M) values at the given temperature. In this example, -9.1 kcal mol^-1 has been predicted for ΔG, corresponding to a K_d of 2.1e-07 M at 25 °C.
         
      2. Prediction details
         
         1. Number of ICs calculated within a threshold of 5.5 Å and % NIS classified according to the charged/polar/apolar character of the amino acids are reported. In this case, for example, there are 7 ICs between charged and polar residues and the % NIS charged atoms is 20.48.
            
         2. Further, the full table (format .txt) of ICs is provided and can be viewed by clicking on the link reported under ‘Table of the ICs at the interface’. The format of the table is the following:
            #chain1 #aa1 #res_num1 #chain2 #aa2 #res_num2
            H → THR → 33 P THR 210
            In which chain ID, residue type and residue number are reported for both residues interacting in Protein 1 and Protein 2.
            
      3. Download outputs
         In this foldable menu, it is possible to download a ready-to-run Pymol script (.pml) (http://www.pymol.org) that will highlight the interaction interface by displaying and coloring the interacting residues, see Figure 2. Further, it is possible to download a compressed file (.tgz) with all the result files. 
         
         
         Figure 2. A three-dimensional representation of the complex 1E6J with the color-coding of the PRODIGY script (.pml). This script can be downloaded from the PRODIGY output page. Interactor 1 is shown in light pink (chains L and H in this example) and Interactor 2 in light blue (chain P), respectively. The interacting residues are represented in sticks in blue and dark pink for Interactor 1 and Interactor 2, respectively.
         
         
3. Useful information
   
   1. Make sure to check and input the correct chain_IDs for the PDB file that you are uploading/retrieving: chain_IDs have to be present in the file, and correspond to the chains that are interacting. In this example, the FAB has two chains labeled as L and H, and both of them are interacting with the HIV1 capsid protein, which is labeled as chain P.
      
   2. PRODIGY can deal with files consisting of an ensemble of structures (e.g., as is typical for NMR structures). In the current implementation, only the first model will be considered for prediction. If you wish to analyze every models present in such an ensemble you should split the PDB file into single-model PDB files and submit them all as an archive file. A collection of useful Python scripts for the manipulation of PDB files, such as splitting of ensemble file, residue renumbering, changing chain ID and so on, can be found in our freely available pdb-tools GitHub repository available online at https://github.com/haddocking/pdb-tools.
      
   3. The PRODIGY web-server currently only supports the 20 canonical amino acids.
      
   4. Information about the web-server input/output, the prediction method and its performance, and the dataset used for training/testing the method can be found online under the Manual/Method/Dataset PROGIDY pages respectively. These are reachable through the corresponding tabs located at the beginning of each page. 
      
      
4. Distribution/Software download
   The PRODIGY web-server is made freely available to the scientific community at: http://milou.science.uu.nl/services/PRODIGY. The prediction scripts are also available from our GitHub repository for local setup and usage at: https://github.com/haddocking/binding_affinity.
   The collection of software developed by the HADDOCK group can be found at: http://www.bonvinlab.org/software.
   The freeSASA software (Mitternacht, 2016) used to calculate the solvent accessible surface area can be downloaded from http://freesasa.github.io.

Notes

To run the ready-to-run Pymol script (.pml) provided by PRODIGY (see step B2c), open a Pymol session with the PDB code that you submitted to PRODIGY and follow one of the possible options:

1. From the bar menu of Pymol, choose File  Run and navigate in the directory where the PRODIGY Pymol script has been saved. Then select the .pml file clicking on ‘Open’.
   
2. In the Pymol terminal bar, type @ followed by the .pml file. Please note, if the Pymol session is not open in that folder, the user will need to type the full path. For example: @home/my_path/prodigy_pymol_script.pml
   

Acknowledgments

This protocol has been adapted from: Vangone and Bonvin (2015) and Xue et al. (2016). Anna Vangone was supported by H2020 Marie-Skłodowska-Curie Individual Fellowship MCSA-IF-2015 [BAP-659025].

References

1. Berman, H., Henrick, K. and Nakamura, H. (2003). Announcing the worldwide Protein Data Bank. Nat Struct Biol 10(12): 980.
   
2. Hubbard, S. J. and Thornton, J. M. (1993). Naccess. Computer Program.
   
3. Kastritis, P. L. and Bonvin, A. M. (2010). Are scoring functions in protein-protein docking ready to predict interactomes? Clues from a novel binding affinity benchmark. J Proteome Res 9(5): 2216-2225.
   
4. Kastritis, P. L., Moal, I. H., Hwang, H., Weng, Z., Bates, P. A., Bonvin, A. M. and Janin, J. (2011). A structure-based benchmark for protein-protein binding affinity. Protein Sci 20(3): 482-491.
   
5. Mitternacht, S. (2016). FreeSASA: An open source C library for solvent accessible surface area calculations. F1000Res 5: 189.
   
6. Vangone, A., and Bonvin, A. M. J. J. (2015). Contacts-based prediction of binding affinity in protein-protein complexes. eLife 4: 291.
   
7. Vreven, T., Moal, I. H., Vangone, A., Pierce, B. G., Kastritis, P. L., Torchala, M., Chaleil, R., Jimenez-Garcia, B., Bates, P. A., Fernandez-Recio, J., Bonvin, A. M. and Weng, Z. (2015). Updates to the integrated protein-protein interaction benchmarks: Docking benchmark version 5 and affinity benchmark version 2. J Mol Biol 427(19): 3031-3041.
   
8. Xue, L. C., Rodrigues, J. P., Kastritis, P. L., Bonvin, A. M. and Vangone, A. (2016). PRODIGY: a web server for predicting the binding affinity of protein-protein complexes. Bioinformatics.