How to perform site-specific docking using Pyrx?
Pyrx  is a virtual screening software that allows the docking of multiple ligands with a target protein. Previously, we provided an article on the installation of Pyrx in Ubuntu. In this article, we are going to perform protein-ligand docking using Pyrx.
Here, we will dock one ligand only for demonstration. For virtual screening of multiple ligands, follow the same steps as explained in this tutorial except you will have to load all ligands.
We are going to use the same protein as we have used in our previous docking tutorials, i.e., human serum albumin (HSA). The PDB ID is 2BXA. It is complexed with a compound called 3-carboxy-4-methyl-5-propyl-2- furanpropanoic acid (CMPF).
- Download the structure of HSA from PDB (www.rcsb.org).
- You can either download the PDB format as usual or download the PDBx/mmCIF format if you face any problems regarding the loading of the structure.
- Save this file in a folder.
We are going to dock Sodium Octanoate (SO) with the HSA.
- Go to Pubchem (www.pubchem.ncbi.nlm.nih.gov).
- Search for the compound “Sodium Octanoate”.
- Download the 3D conformer of the same in SDF format.
- You can download the structure from another database as well such as ZINC.
Preparing the protein
To prepare the protein for docking, you need to know the binding site/ binding residues. They can be present in a single chain of a protein or might involve two. Read the associated paper published for that structure. It consists of all structural details about it. If there is no literature available on the binding site of that protein, then go for the prediction. Read this article for more information.
In the present case, you can also get all binding residues by visualizing the HSA-CMPF complex in Pymol and note down all interacting residues. But always try to find the binding sites in published papers.
Since we already know the binding sites of HSA, we can proceed to prepare our protein.
- Open the PDB file and remove HETATOMS.
The structure we are using is a crystal structure complexed with ligand(s), therefore, to dock the desired ligand with the protein in that particular position we need to remove the bound ligand by removing hetatoms from the PDB file. If we will dock our ligand without removing the already complexed ligand, then we will not get the correct results. We can also easily remove ligands by visualizing the protein in PyMol.
Read this article to know how to remove hetatoms.
- Remove extra chains
We will keep only one of the four chains (here, Chain A was taken) and remove the rest of the three chains. The chains are removed from the protein structure just to avoid complexity. Read this article to know how to remove chains.
- Save the file
Now save the file as “protein.pdb”.
We have prepared our protein structure to proceed further for docking. Now let’s prepare our ligand.
Preparing the ligand
Since we need the protein and the ligand in a .pdb format, therefore, we have to convert the ligand from SDF to pdb format. We will use PyMol for this purpose and never use online converters because they may ruin your ligand file. You can also use Openbabel.
- Open PyMol and open the downloaded ligand.
- Click on “File” –> “Save Molecule” –> select the molecule –> click “OK”.
- Save it in the same folder as ‘ligand.pdb’ or with any other name.
Performing docking using Pyrx
We will be using Autodock Vina for docking. It is provided in Pyrx that you can dock using the Vina algorithm. Run Pyrx GUI and follow the steps explained below.
Loading protein and ligand
- Go to ‘File’ –> ‘Load Molecule’ or just click on the very first icon given on top left corner.
- Select the downloaded protein structure. Here, named as ‘protein.pdb’.
- Again go to ‘File’ –> ‘Load Molecule’ or just click on the very first icon given on top left corner.
- Select the converted ligand structure. Here, named as ‘ligand.pdb’.
Defining protein and ligand
You will see the loaded protein and ligand under the ‘Molecules’ tab. Now we have to define which one is protein and which one is the ligand. To do that
- Right-click on protein –> ‘Autodock’ –> ‘Make Macromolecule’.
- Right-click on ligand –> ‘Autodock’ –> ‘Make Ligand’.
After that you will notice under the ‘Autodock’ tab, it has prepared their PDBQT files itself.
Defining grid box
Look at the ‘Controls’ panel given at the bottom under the ‘Select molecule’ tab. You will see “1 ligand selected” and the path of the selected protein.
- Now, click the ‘Forward’ button present at the bottom right corner. You will see a grid box on the display screen.
- Go back to the ‘Molecules’ tab on the right side. Click on the “+” sign of the loaded protein. You will see all the residues present in the chain.
- Select the binding residues by right-clicking on the residue –> Display –> Label –> Atoms. You will see the atoms appearing on the protein.
- Now adjust the grid box accordingly such that all the selected residues come inside the box. You don’t have to enclose the ligand inside the grid box.
Running Vina for docking
You can set the exhaustiveness by entering an appropriate number in the box given on the bottom left. After adjusting everything, click on the ‘Forward’ button. It will start docking and will display the processing. After the docking is finished, you will see the poses and their binding affinities in the bottom panel. It will show all nine poses along with other values such as RMSD.
Now to analyze these results, you will have to open the PDB/PDBQT file of the protein and the vina output file in Pymol. For analyzing results, read this article. You can also analyze these results in Discovery Studio Visualizer, read here to analyze docking results in DS Visualizer.
vs_Analysis.py: A Python Script to Analyze Virtual Screening Results of Autodock Vina
How to install Autodock4 on Ubuntu?
Autodock suite is used for docking small molecules . Recently, Autodock-GPU  is developed to accelerate the docking process. Its installation is described in this article. In this tutorial, we will install Autodock 4.2.6 on Ubuntu.
What values are considered as good or bad in computational docking?
After performing computational docking, a question that comes to mind most is “what docking score is considered good or bad”. In this article, we will discuss this in detail. (more…)
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?
A configuration file is one of the required files for docking using Autodock Vina. In this article, we are going to generate a config file using Autodock Tools GUI . (more…)
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?
Pyrx is a bioinformatics tool to perform virtual screening . We previously provided an article on performing simple protein-ligand site-specific docking using Pyrx. In this article, we are going to perform virtual screening using Pyrx. (more…)
How to perform metal ion-protein docking using idock?
Previously, we provided a tutorial on the installation of idock on Ubuntu (Linux). In this article, we are going to demonstrate the docking of a metal ion (such as Zn, Mg, Fe, etc.,) with a protein using idock. (more…)
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…)
How to analyze HADDOCK results using Pymol script generated from PRODIGY?
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 . (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 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?
This is a demonstration video of our previously mentioned Python script (vs_analysis.py) for virtual screening results analysis of Autodock Vina. (more…)
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…)
Blind Docking vs Site-specific Docking
We have provided several tutorials on docking including both types: blind docking & site-specific docking using Autodock Vina. In this article, we are going to compare both methods of docking. (more…)
Most widely used software for docking results image generation
Computational docking is an important technique in bioinformatics. The docking result interpretation and representation are important as well. There are various software available for docking result analysis and offer to generate high-quality images for publications. In this article, we will mention such a few software for the same purpose. (more…)
You must be logged in to post a comment Login