Benchmarking different docking protocols for predicting the binding poses of ligands complexed with cyclooxygenase
Sara Shamsian1,2, Babak Sokouti3, Siavoush Dastmalchi2,3,4
1Student Research Committee, Tabriz University of Medical Sciences, Tabriz, 5165665931, Iran.
Introduction:
Non-steroidal anti-inflammatory drugs (NSAIDs) constitute an important class of pharmaceuticals acting on cyclooxygenase COX-1 and COX-2 enzymes. Due to their numerous severe side effects, it is necessary to search for new selective, safe, and effective anti-inflammatory drugs. In silico design of novel therapeutics plays an important role in nowadays drug discovery pipelines. In most cases, the design strategies require the use of molecular docking calculations. The docking procedure may require case-specific condition for a successful result. Additionally, many different docking programs are available, which highlights the importance of identifying the most proper docking method and condition for a given problem.
Methods:
In the current work, the performances of five popular molecular docking programs, namely, GOLD, AutoDock, FlexX, Molegro Virtual Docker (MVD) and Glide to predict the binding mode of co- crystallized inhibitors in the structures of known complexes available for cyclooxygenases were evaluated. Furthermore, the best performers, Glide, AutoDock, GOLD and FlexX, were further evaluated in docking-based virtual screening of libraries consisted of active ligands and decoy molecules for cyclooxygenase enzymes and the obtained docking scores were assessed by receiver operating characteristics (ROC) analysis.
Results:
The results of docking experiments indicated that Glide program outperformed other docking programs by correctly predicting the binding poses (RMSD less than 2 Å) of all studied co-crystallized ligands of COX-1 and COX-2 enzymes (i.e., the performance was 100%). However, the performances of the other studied docking methods for correctly predicting the binding poses of the ligands were between 59% to 82%. Virtual screening results treated by ROC analysis revealed that all tested methods are useful tools for classification and enrichment of molecules targeting COX enzymes. The obtained AUCs range between 0.61-0.92 with enrichment factors of 8 - 40 folds.
Conclusion:
The obtained results support the importance of choosing appropriate docking method for predicting ligand-receptor binding modes, and provide specific information about docking calculations on COXs ligands.
More Related Videos
10:21Author Spotlight: Streamlining Protein Target Prediction and Validation via Molecular Docking and CETSA
Published on: February 23, 2024
11:06Network Pharmacology Prediction and Metabolomics Validation of the Mechanism of Fructus Phyllanthi against Hyperlipidemia
Published on: April 7, 2023
Related Concept Videos
02:18The Equilibrium Binding Constant and Binding Strength
02:40Ligand Binding Sites
Protein-ligand interactions are quite specific; even though numerous potential ligands surround a cellular protein at any given time, only a particular ligand can bind to that protein. Moreover, a ligand binds only to a dedicated area on the surface of the protein, known as the...
00:49Ligand Binding and Linkage
01:49Conserved Binding Sites
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
02:04Protein-protein Interfaces
01:22Protein-Drug Binding: Determination Methods
Indirect methods involve isolating the bound drug from its free form in biological samples such as blood, serum, or plasma. These techniques aim to measure the percentage of drugs bound to proteins. Equilibrium dialysis is a commonly used method where the free drug concentration at equilibrium is measured by separating the bound...