Docking
[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 [1].
For docking zinc metalloprotein, we will be using the Autodock4Zn force field. For that, you need to install the Autodock-Vina-develop package [2,3] as shown below. Download the package from GitHub. For detailed steps on its installation, read this article.
Now, let’s start with the docking of metalloprotein with a ligand. We will use ‘Human protein farnesyltransferase’ complexed with L-778,123 and FPP (PDB ID: 1s63). Download this protein and then we will prepare it using python scripts.
Preparing Protein
This protein contains two chains A & B. You can remove one chain and other HETATMS such as HOH. Do not remove Zn from the PDB file. Now prepare the PDBQT file of protein using a Python script, namely, ‘prepare_receptor.py’. You can find this script here. You can also use Autodock tools for the same.
$ prepare_receptor.py -r 1s63.pdb -v -o protein.pdbqt
It will output the pdbqt file, namely, ‘protein.pdbqt‘. After that run the following command to add the tetrahedral zinc pseudo atoms.
$ python /home/user/Downloads/AutoDock-Vina-develop/example/autodock_scripts/zinc_pseudo.py -r protein.pdbqt -o protein_tz.pdbqt
Replace user
to your username in the above command. After running this command, it will display something like this:
$ Wrote 1 TZ atoms on protein_tz.pdbqt.
Your protein has been prepared, let’s prepare the ligand.
Preparing Ligand
Here, the ligand is 778 as given in 1s63 PDB structure. To download the ligand, click on ‘Download Instance Coordinates‘ below the ligand name under ‘Small Molecule‘ section. Click on SDF format. Download it in the same directory and we are going to save it as ‘ligand.sdf‘. Convert it into pdb format using Pymol. For more information, read here.
Now, let’s prepare the pdbqt file of ligand using Python script, namely, ‘prepare_ligand.py’. You can download this file from here.
$ prepare_ligand.py -l ligand.pdb -v -o ligand.pdbqt
Generating Affinity Maps
To generate affinity maps, we will use the ‘prepare_gpf4zn.py‘ script. You can find this script inside the AutoDock-Vina-develop/example/autodock_scripts directory. In this command, we will define the grid box dimensions and Autodock4Zn forcefield.
$ python /home/user/Downloads/AutoDock-Vina-develop/example/autodock_scripts/prepare_gpf4zn.py -l ligand.pdbqt -r protein_tz.pdbqt -o protein_tz.gpf -p npts=40,30,50 -p gridcenter=18,134,-1 –p parameter_file=AD4Zn.dat
Now, run autogrid4 to generate affinity maps.
$ autogrid4 -p protein_tz.gpf -o protein_tz.glg
After this command, we have successfully created protein and ligand files in pdbqt format and all information regarding forcefield is incorporated into the affinity maps. Now, we can run docking.
Docking using Autodock Vina
We will use --scoring ad4
argument to define Autodock4 forcefield as shown below. Run docking command:
$ vina --ligand 1s63_ligand.pdbqt --maps protein_tz --scoring ad4 --exhaustiveness 32 --out 1s63_ligand_ad4_out.pdbqt
It will take the usual time to finish docking. Soon, you will get docking results.
References
- Trott, O., & Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. Journal of computational chemistry, 31(2), 455-461.
- J. Eberhardt, D. Santos-Martins, A. F. Tillack, and S. Forli. (2021). AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings. Journal of Chemical Information and Modeling.
- O. Trott and A. J. Olson. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. Journal of computational chemistry, 31(2), 455-461.
Docking
[Tutorial] Performing docking using DockingPie plugin in PyMOL.
DockingPie [1] is a PyMOL plugin to perform computational docking within PyMOL [2]. In this article, we will perform simple docking using DockingPie1.2.
Docking
How to install the DockingPie plugin on PyMOL?
DockingPie [1] is a plugin of PyMOL [2] made to fulfill the purpose of docking within the PyMOL interface. This plugin will allow you to dock using four different algorithms, namely, Vina, RxDock, SMINA, and ADFR. It will also allow you to perform flexible docking. Though the installation procedure is the same for all OSs, in this article, we are installing this plugin on Ubuntu (Linux).
Docking
[Tutorial] Installing Pyrx on Windows.
Pyrx [1] is another virtual screening software that also offers to perform docking using Autodock Vina. In this article, we will install Pyrx on Windows. (more…)
Docking
How to install Autodock4 on Ubuntu?
Autodock suite is used for docking small molecules [1]. Recently, Autodock-GPU [2] 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.
Docking
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…)
Bioinformatics Programming
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…)
Docking
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…)
Docking
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 [1]. (more…)
Docking
How to install Autodock on Ubuntu (Linux) with CUDA GPU support?
Autodock [1] is most widely used for docking. To accelerate the docking process, especially, during virtual screening, Autodock-GPU [2] provides great help. In this article, we are going to install Autodock-GPU on Ubuntu. (more…)
Docking
Autodock Vina: Uses & Applications
Autodock Vina [1] 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…)
Docking
How to perform virtual screening using Pyrx?
Pyrx is a bioinformatics tool to perform virtual screening [1]. 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…)
Docking
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…)
Docking
How to install idock on Ubuntu?
idock [1] 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…)
Docking
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 [1]. In this article, we are going to demonstrate the HADDOCK results analysis using a Pymol script generated from the PRODIGY server [2]. (more…)
Docking
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…)
Docking
Protein-protein docking using HADDOCK2.4 web server
HADDOCK2.4 web server (https://bianca.science.uu.nl/haddock2.4/) [1] 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…)
Docking
How to perform site-specific docking using Pyrx?
Pyrx [1] 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. (more…)
Software
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…)
Docking
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…)
Docking
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…)
You must be logged in to post a comment Login