In one of our previously published articles, we demonstrated protein-protein docking using HADDOCK2.4 . In this article, we are going to demonstrate the HADDOCK results analysis using a Pymol script generated from the PRODIGY server .
PRODIGY (Protein Binding energy prediction) webserver provides a variety of functions for the binding affinity prediction of biological complexes as well as the identification of biological interfaces from crystallographic ones . It can also be used to closely analyze the docking output of the HADDOCK web server.
After you get results from the HADDOCK2.4 web server, download them as a tar file and later extract them. You will find several structures in PDB format. They would be named as “cluster1_1”, “cluster1_2”, “cluster2_1”, and so and so forth. The best cluster structure is provided by the HADDOCK itself. You can select that structure for further analysis.
Submitting structure to PRODIGY
- First, select the best structure as provided by HADDOCK.
- Go to PRODIGY webserver.
- Select an appropriate tab (Protein-Protein, or Protein-small molecule, or biological/crystallographic). For this tutorial, we will select protein-protein because we docked a protein to another protein in our last tutorial on HADDOCK2.4.
- After that, upload your selected structure, let’s say cluster1_1.
- Write the names of the chains for both proteins. Here, it is chain A and chain B.
- Set the appropriate temperature, we are going to use the default one (25 degrees C).
- Provide a name to your job, it’s optional though.
- Now, submit it.
Within a few minutes, your job will be finished and a result page will appear. The result page will display the binding affinity and Kd values, along with the other prediction details such as the Number of Interfacial Contacts, Non-Interacting Surface, and a table of the interfacial contacts at the interface. Below this, you will find a Pymol script namely, “Pymol ICs-based representation script (.pml)” under the section “Download Outputs”.
Analysis using Pymol script
Now, to analyze your output,
- open the Cluster1_1.pdb file in Pymol.
- open the downloaded .pml script in Pymol.
- As soon as you open the Pymol script, you will see your proteins in different colors and the interacting residues with white backgrounds. The interacting residues would appear as selected.
- After that, you can represent your docked complex as per your requirements. For example, as surfaces, cartoon models, or highlight interacting residues, as so on.
You can also read the table for interacting residues in the docked complexes. That’s how you analyze HADDOCK results using the Pymol script generated by the PRODIGY web server.
- De Vries, S. J., Van Dijk, M., & Bonvin, A. M. (2010). The HADDOCK web server for data-driven biomolecular docking. Nature protocols, 5(5), 883-897.
- Vangone A. and Bonvin A.M.J.J. “Contact-based prediction of binding affinity in protein-protein complexes”, eLife, 4, e07454 (2015).
- Xue L., Rodrigues J., Kastritis P., Bonvin A.M.J.J.*, Vangone A.*, “PRODIGY: a web-server for predicting the binding affinity in protein-protein complexes”, Bioinformatics, doi:10.1093/bioinformatics/btw514 (2016).
[Tutorial] Installing Pyrx on Windows.
How to install Autodock4 on Ubuntu?
What values are considered as good or bad in computational docking?
How to sort binding affinities based on a cutoff using vs_analysis.py script?
Previously, we have provided a Python script (vs_analysis.py) to analyze the virtual screening (VS) results of Autodock Vina. Now, we have updated this script to sort binding affinities based on user inputted cutoff value. (more…)
Virtual Screening using Autodock Vina: Frequently Asked Questions & Answers for Starters
Virtual Screening (VS) is one of the important techniques in bioinformatics. It can be easily performed using Autodock Vina. We have provided detailed articles on this topic. In this article, we are trying to answer some FAQs for beginners. (more…)
[Tutorial] How to perform docking of zinc metalloproteins using Autodock Vina?
Proteins containing zinc atoms are docked in a different way than that of the normal simple proteins. Zinc atoms must be considered by a force field during the docking process. In this article, we are going to dock zinc metalloprotein with a ligand using Autodock Vina . (more…)
How to generate config file for docking using Autodock Tools?
How to install Autodock on Ubuntu (Linux) with CUDA GPU support?
Autodock  is most widely used for docking. To accelerate the docking process, especially, during virtual screening, Autodock-GPU  provides great help. In this article, we are going to install Autodock-GPU on Ubuntu. (more…)
Autodock Vina: Uses & Applications
Autodock Vina  is one of the most widely used bioinformatics software for computational docking. For beginners and those who are new to the field of bioinformatics, such software may appear confusing without having an initial idea of what this software actually used for. In this article, we have discussed some of the primary uses and applications of Autodock Vina software. (more…)
How to perform virtual screening using Pyrx?
How to perform metal ion-protein docking using idock?
How to install idock on Ubuntu?
idock  is a multithreaded software based on Autodock Vina. It is a virtual screening tool for flexible ligand docking. It also supports 27 different chemical elements including zinc, magnesium, iron, calcium, etc. In this article, we are going to install idock on Ubuntu. (more…)
Installing PatchDock on Ubuntu (Linux)
PatchDock is a docking algorithm that is based on the shape complementarity principle [1,2]. It performs molecular docking of any two types of molecules including proteins, DNA, drugs, and peptides. In this article, we are going to install PatchDock on Ubuntu (Linux). (more…)
Protein-protein docking using HADDOCK2.4 web server
HADDOCK2.4 web server (https://bianca.science.uu.nl/haddock2.4/)  is used for protein-protein docking and their modeling. Not only proteins, but it can also process peptides, small molecules, and nucleic acids. In this article, we are going to perform protein-protein docking. (more…)
How to perform site-specific docking using Pyrx?
How to install sdsorter on Ubuntu (Linux)?
sdsorter is a tool that helps in the easy manipulation of SDF files. It takes a single file as input and also provides a single output file. In this article, we are going to install sdsorter on Ubuntu (Linux). (more…)
Video Tutorial: How to use vs_analysis Python script for virtual screening result analysis of Autodock Vina?
Basic Methodology to Predict Antigen-Antibody Interactions in silico
Antigen-antibody interactions play an important role in protecting our body from foreign molecules. This is applied in vaccine development against a disease. In this article, we are going to mention a simple methodology used in studying/ identifying interactions between antigen-antibody complex. (more…)
What is Virtual Screening in Bioinformatics?
If you are a beginner in the field of bioinformatics, you must have heard of the term “virtual screening” (VS). In this article, we will discuss the basic meaning of VS and its significance in bioinformatics. (more…)