Molecular Docking using different Tools

 

Soham Pawar*, Chaitrali Kulkarni, Puja Gadade, Supriya Pujari, Surajkumar Kakade,

S.H. Rohane, V.K. Redasani

Department of Pharmacy, Yashoda Technical Campus, Satara, 415003.

*Corresponding Author E-mail: sohamsp105@gmail.com

 

ABSTRACT:

Molecular docking is a powerful and effective tool in modern era of computer aided drug design. The purpose of ligand-protein docking is to predict the predominant binding mode(s) of a ligand with a protein of known three-dimensional structure. Successful docking methods search high-dimensional spaces effectively and use a scoring function that correctly ranks candidate dockings. Docking can be used to perform virtual screening of wide variety of compounds, rank the results, and propose structural hypotheses of how the ligands inhibit the target, which is invaluable in lead optimization. This review represents the overview of advanced molecular docking using various software. It analyses the various approaches regarding molecular docking like the use of machine learning algorithms in molecular docking, different physicochemical aspects related to ligand-protein complex. These recent improvements in modern technology affect the whole area of healthcare and welfare science.

 

KEYWORDS: Molecular Docking, Drugdesign, Autodock, Pymol, Discovery Studio.

 

 


INTRODUCTION:

The drug discovery process originated in the 19th century by John Langley in 1905, when he proposed the theory of receptive substances. Ehrlich was awarded by Nobel Prize in 1908. In 1960, Hansch and Fujita introduced the concept of QSAR (quantitative structure–activity relationship). In the drug discovery process, the first step is to identify an appropriate ‘druggable’ target, which can be a biomolecule or a protein receptor that is explicitly associated with a disease condition or pathology1.

 

QSAR Quantitative structure-activity relationships (QSAR):

Quantitative structure-activity relationshipsattempt to identify and quantify the physicochemical properties of a drug and to see whether any of these properties have an effect on the drugs biological activity.

 

 

It helps to find consistent relationship between biological activity and molecular properties, so that these “rules” can be used to evaluate the activity of new compounds. It helps in understanding QSAR regarding electronic effects, steric effects and lipophilicity. In 3D QSAR, 3D properties of a molecule are considered as whole rather than considering individual substituent2,3. This method involves the analysis of the quantitative relationship between the biological activity of a set of compounds and their three-dimensional properties using statistical correlation methods. It revolves around the important features of a molecule, its overall size and shape, and its electronic properties4,5.

 

Physicochemical parameters:

In QSAR modelling the prediction contains different Phsicochemical properties and theoretical molecular descriptors of various chemicals.

 

Various parameters used in QSAR studies are:

Hydrophobicity: partition coefficient, π- substitution constant

 

 

Electronic parameter: Hammet constant, dipole moment.

Stearic parameter: Taft’s constant, Ver loop steric parameter.

 

Molecular docking:

Molecular docking is a key tool in structural molecular biology and computer-assisted drug design. The goal of ligand-protein docking is to predict the predominant binding mode (s) of a ligand with a protein of known three-dimensional structure6,7.

 

Three steps of molecular docking

(1) Definition of the structure of the target molecule

(2) Location of the binding site

(3) Determination of the binding mode

 

The molecular docking is used for demonstration of feasibility of any biological reaction as it is carried out before experimental part of any investigation. This method is one of the most frequently used method in structure-based drug design due to its ability to predict the binding conformation of small molecule ligands to the appropriate target binding sites.

 

Auto-Dock:

AutoDock is an automated suite of protein-ligand docking tools. AutoDock tools are abbreviated as ADT. It is designed to predict the protein interactions with small molecules such as drug molecule and substrate. The application of this tool is immense, ranging from structure-based drug design, lead molecule optimization, protein-ligand docking, protein-protein docking, analysis and validation of mechanism of action of drug molecules. AutoDock4 has two key programs to be executed, i.e., Autogrid4 and AutoDock4. Autogrid4 prepares a grid map of the amino acids presents within the Grid Box defined by the user. AutoDock4 then analyzes the interactions of those amino acids with the ligand molecule. AutoDock4 predicts the free binding energy with a scoring function based on the AMBER force field and linear regression analysis, additionally by referring to a large set of library data of known protein ligand interactions with their inhibition constants that were used in AutoDock3. It is found that the standard error in free binding energy in AutoDock4 is approximately 2.5kcal/mol8.

 

This protocol is presented to aid students and research scholars those who are interested in using AutoDock4 for learning and or research purpose. This protocol is considered to be adequate and effective by the authors, to use AutoDock4 for protein-ligand docking analysis. The discussed protocol can be used to study the interactions of selected ligand molecule (drug molecule) with chosen protein targets (9–11). Users, however, have to do a little preparative work on the protein molecule pdb file. Downloaded 3D structure of proteins from RCSB website has to be edited before docking in AutoDock4. Although a small discussion on the preparation of protein molecules is presented in this tutorial, users are advised to do further reading in the preparation of protein structures12-14.

 

Installation:

Download the installation files “autodocksuite-4.2.5.1-i86Windows.exe and mgltools_win32_1.5.6_Setup.exe” from Auto Dock website (http://autodock.scripps.edu/).

Run “mgltools_win32_1.5.6_Setup.exe” file and leave al the settings default and the installation folder should not be altered. Finish the installation with the default settings. Once the installation is finished, “AutoDockTools-1.5.6.exe” will be executed by default. Close the application and start the “AutoDockTools-1.5.6.exe” once again. Once the application is opened for the second time, a new folder “. mgltools” will be created automatically in the users folder (C:\Users\Guest\.mgltools). To successfully run the AutoDock software, the protein and ligand structures in their “.pdb” format have to be present within the “.mgltools” folder. Next, run the “autodocksuite-4.2.5.1-i86Windows.exe” file and select a custom folder for this. At the end of the installation, there will be two execution files created in the selected custom folder (i.e., AutoDock4. exe and autogrid4.exe). At this stage, the software is ready to be executed.

 

Protein Molecule Preparation:

Protein molecules can be downloaded from Protein Data Bank website (www.rcsb.org). The downloaded protein structure in their “.pdb” format has to be edited to remove the non-amino acid residues, such as water molecules, ions, ligands that are in the complex. These can be removed using either PyMol software or WordPad15-17.

 

Executing AutoDock:

AutoDock software calculates and predicts the interaction between the ligand molecule and protein molecule based on predefined parameters. To be precise, the interactions between the molecules will be calculated at a user specified region in the protein. This region can be defined by users, using the Grip map option. Ultimately, the software predicts the interaction and binding energy of the ligand molecule and the amino acids present within the Grid Box only. Thus setting, the Grid Box at the binding site or active site or other essential regions of the protein is important.

 

Analysis in AutoDock can be divided into following categories; (a) Initializing molecules; (b) Running Auto Grid; (c) Running AutoDock; (d) Analyzing Interaction energy.

 

 

(a) Initializing molecules - Initializing the molecule mainly includes addition of hydrogen atoms and addition of Kolman charge to the protein molecule. While for the ligand molecule, addition of Gasteiger charge, identifying aromatic carbons, detecting rotatable bonds, and setting TORSDOF value. The protein has to be initialized manually, while the ligand is automatically initialized when opened in the tool. Once the protein molecule is opened, it is important to change the view of the protein. It is essential that the protein is in “Ribbon view,” failing to do so, might result in errors in the later steps.

 

(b) Running Auto Grid - AutoGrid has to be executed, to define the region/area in the protein to be analyzed for the interaction with the ligand molecule. In general, the region of interaction could be identified using prediction tools such as Q-Site finder and MetaPocket to identify the binding pockets on the surface of the proteins. Then, the GridBox is set in AutoDock to cover the identified binding sites. AutoDock only analyzes the interactions of ligand molecule and the amino acids that are present within the GridBox. So, setting up, the GridBox is a crucial step.

 

(c)Running AutoDock - Once Auto Grid is successfully completed, AutoDock can be executed. AutoDock calculates the interactions of the amino acid within the GridBox and ligand molecule.

 

AutoDock calculations are performed in several steps:

 

Step 1—Coordinate File Preparation:

AutoDock4.2 is parameterized to use a model of the protein and ligand that includes polar hydrogen atoms, but not hydrogen atoms bonded to carbon atoms. An extended PDB format, termed PDBQT, is used for coordinate files, which includes atomic partial charges and atom types. The current AutoDock force field uses several atom types for the most common atoms, including separate types for aliphatic and aromatic carbon atoms, and separate types for polar atoms that form hydrogen bonds and those that do not. PDBQT files also include information on the torsional degrees of freedom. In cases where specific sidechains in the protein are treated as flexible, a separate PDBQT file is also created for the sidechain coordinates. AutoDockTools, the Graphical User Interface for AutoDock, may be used for creating PDBQT files from traditional PDB files.

 

Step2—AutoGrid Calculation:

Rapid energy evaluation is achieved by recalculating atomic affinity potentials for each atom type in the ligand molecule being docked. In the AutoGrid procedure the protein is embedded in a three-dimensional grid and a probe atom is placed at each grid point. The energy of interaction of this single atom with the protein is assigned to the grid point. AutoGrid affinity grids are calculated for each type of atom in the ligand, typically carbon, oxygen, nitrogen and hydrogen, as well as grids of electrostatic and desolvation potentials. Then, during the AutoDock calculation, the energetics of a particular ligand configuration is evaluated using the values from the grids.

 

Step 3—docking using AutoDock:

Docking is carried out using one of several search methods. The most efficient method is a Lamarckian genetic algorithm (LGA), but traditional genetic algorithms and simulated annealing are also available. For typical systems, AutoDock is run several times to give several docked conformations, and analysis of the predicted energy and the consistency of results is combined to identify the best solution.

 

Step 4—Analysis using Auto Dock Tools:

AutoDockTools includes a number of methods for analysing the results of docking simulations, including tools for clustering results by conformational similarity, visualizing conformations, visualizing interactions between ligands and proteins, and visualizing the affinity potentials created by AutoGrid.

 

Analysing Interaction Energy:

Once AutoDock4.exe is successfully executed. The result will be given as the ten best confirmations. These can be viewed in the analyze options. The confirmations can be viewed in the order of their free energy binding, by choosing the “Play, ranked by energy” option. The ten conformations can be viewed by changing the conformations number. The interaction energy of the given conformation can also be viewed. The number of hydrogen bonds formed between the ligand and protein can be viewed.

 

Visualisation by Pymol:

Open pymol>> select protein file (“.pdb”file ) set the protein  view to ribbon and surface

Remove water molecules>>Press Ctrl L>>remove unwanted ligand.

 

File>>export molecule>>save>>rename as protein1 and use PDB (*.pdb*.pdb.gz) extension to save.

 

Pymol-

Open pymol>> add protein1>>add ligand1>>first docked molecule>>file>>export molecule>>save rename (complex1)>>use PDB (*.pdb*.pdb.gz) extension to save.(Fig)

 

2D Study by Discovery studio:

Open discovery studio>>add complex1 file>>receptor-ligand interaction>>ligand interactions>> show distance (bond length), show type (bond angle).

 

To see the names of amino acid:

Change the visibility of receptor and ligand >>interacting atom>> right click>>label >> choose amino acid in object section>> select full name in attribute section >> apply >> ok

 

To see the theoretical part:

Click on ^ >> Click on non-bond>>it will show various aspects of bonds

 

To see ligand non bond monitor:

Click on >situated on left side of 3D diagram >>click on > near ligand non bond monitor

 

To see 2D diagram:

Select show 2D diagram

 

CONCLUSION:

Molecular docking is an inexpensive, safe and easy to use tool which helps in investigating, Interpreting, explaining and identification of molecular properties using three-dimensional and two-dimensional structures as well. It is advance techniquewhich allow us to characterize the behaviour of small molecules in the binding site of target proteins as well as to elucidate fundamental biochemical processes.

 

Since different models yield different results, it is necessary to have a small number of standard models which are applicable to very large systems. Molecular docking is used to predict the structural intermolecular complexes formed between two or more constituting molecules. These techniques are used in the field of computational chemistry, computerizedBiology and material used for molecular system ranges from small molecules to large biological molecules and material assembly. It is open-source software for computational drug discovery that can be used to libraries of computer against potential drug targets. By using Pymol and Discovery studio like software a new genre of molecular research opened for world.

 

REFERENCES:

1.      Jejurikar BL, Rohane SH. Drug Designing in Discovery Studio. 2021;14(April):135–8.

2.      Aware D, Rohane S. Review Article a Role of Herbal Drug as an Immunity Booster during Covid-19 Pandemic. 2021; 11(3):206–11.

3.      Galande AK, Rohane SH. Insilico Molecular docking analysis in Maestro Software. 2021;14: 97–100.

4.      Patil NS, Rohane SH. Organization of Swiss Dock : In study of Computational and Molecular Docking Study. 2021;14(April):145–8.

5.      Nangare SA, Rohane SH. Review on Guassion, the General Purpose in Computational Chemistry for Medicinal Chemistry. 2021; 14 :89–91.

6.      Satpute UM, Rohane SH. Efficiency of AUTODOCK: Insilico study of Pharmaceutical Drug Molecules. Asian J Res Chem. 2021; 14(1): 1–5.

7.      Rohane SH, Makwana AG. In silico study for the prediction of multiple pharmacological activities of novel hydrazone derivatives. Indian J Chem -Section B. 2019; 58(3): 387–402.

8.      Hemalatha K, Selvin J, Girija K. Synthesis, In silico Molecular Docking Study and Anti-bacterial Evaluation of some Novel 4-Anilino Quinazolines. 2018; 8(3):125–32.

9.      Jadhav R, More K, Khandekar D, Yamgar R, Sawant S. An Efficient Synthesis of Rabeprazole Sodium. 2011;1(1):13–6.

10.   Otuokere IE, Amaku FJ, Alisa CO. In Silico Geometry Optimization, Excited – State Properties of (Belinostat) and its Molecular Docking Studies with Ebola Virus Glycoprotein. 2015; 5(3): 131–7.

11.   Nandhini S, Radha R, Vadivu R. Docking of Hematoporphyrin on Various Anticancer Drugs Targeting Enzymes. 2016; 6(3):123–30.

12.   Akula G, Nirmala R, Reddy SR. Quinazolinones and evaluation of their Anti-inflammatory activity. 2017; 7(2): 60–2.

13.   Bhalerao R, Rishipathak D, Udavant P, Kshirsagar S. Synthesis and Evaluation of Antidiabetic Activity of Few 5- (4- (5-. 2017; 7(4): 261–4.

14.   Kamala G, Srinivasan N, Shankar KR, Suresh R. Synthesis, Characterization and Antimicrobial Evaluation of N-Mannich Bases of (2- Substituted Phenyl) Benzimidazole Derivatives. 2018; 8(2): 87–93.

15.   Shahare H V, Amrutkar RD. Synthesis, Characterization and Antimicrobial Activity of Diphenylamino Isoxazoline Derivatives. 2018; 8(3): 148–50.

16.   Synthesis, Characterization, Anti Microbial, Analgesic activities of 3- (3-Chlorophenyl) -5-Phenyl-4, 5-Dihydro-1, 2-Oxazole Derivatives. 2017; 7(3): 171–4.

17.   Kumar SS, Melchias G, Ravikumar P, Chandrasekar R, Kumaravel P. Bioinspired synthesis of silver nanoparticles using Euphorbia hirta leaf extracts and their antibacterial activity. 2014; 4(1): 39–43.

 

 

 

 

 

Received on 27.12.2022         Modified on 22.03.2023

Accepted on 19.06.2023   ©Asian Pharma Press All Right Reserved

Asian J. Pharm. Res. 2023; 13(4):292-296.

DOI: 10.52711/2231-5691.2023.00053