Combining fragment docking with graph theory to improve ligand docking for homology model structures

Abstract

Computational protein-ligand docking is well-known to be prone to inaccuracies in input receptor structures, and it is challenging to obtain good docking results with computationally predicted receptor structures (e.g. through homology modeling). Here we introduce a fragment-based docking method and test if it reduces requirements on the accuracy of an input receptor structures relative to non-fragment docking approaches. In this method, small rigid fragments are docked first using AutoDock Vina to generate a large number of favorably docked poses spanning the receptor binding pocket. Then a graph theory maximum clique algorithm is applied to find combined sets of docked poses of different fragment types onto which the complete ligand can be properly aligned. On the basis of these alignments, possible binding poses of complete ligand are determined. This docking method is first tested for bound docking on a series of Cytochrome P450 (CYP450) enzyme-substrate complexes, in which experimentally determined receptor structures are used. For all complexes tested, ligand poses of less than 1 Å root mean square deviations (RMSD) from the actual binding positions can be recovered. Then the method is tested for unbound docking with modeled receptor structures for a number of protein-ligand complexes from different families including the very recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) protease. For all complexes, poses with RMSD less than 3 Å from actual binding positions can be recovered. Our results suggest that for docking with approximately modeled receptor structures, fragment-based methods can be more effective than common complete ligand docking approaches.

Iron(III), copper(II), cadmium(II), and mercury(II) complexes of isatin carbohydrazone Schiff base ligand (H3L): Synthesis, characterization, X‐ray diffraction, cyclic voltammetry, fluorescence, density functional theory, biological activity, and molecular docking studies

“…On the other hand, it does not take into account the experimental factors such as physicochemical properties of the drug (size, concentration, and pKa), properties of the membrane (thickness, area, solution pH, and permeability for the drug), and mechanism of the drug transport. [ 75,76 ] It was also remarked the fine binding affinity in interaction with the binding site residues implying van der Waals interactions, conventional H‐bonding and π‐interactions, π‐alkyl, and π‐sulfur.…”

In silico approach through molecular docking and study ADME on imine derivative compounds as a potential antibacterial agent against Staphylococcus aureus

“…Prediction of physicochemical properties was carried out using SwissADME, which is based on Lipinski's rule of five and requires several physicochemical parameters. Drug-likeness or drug similarity is used to estimate the probability of a molecule becoming an oral drug by considering the qualitative bioavailability of the compound (Lipinski et al, 2001;Sarfaraz et al, 2020).…”

Cleaning of Wastewater Using Crosslinked Poly(Acrylamide-co-Acrylic Acid) Hydrogels: Analysis of Rotatable Bonds, Binding Energy and Hydrogen Bonding

“…The docking method was applied to analyze different interactions, with the receptor and ligand representing the polymer matrix and dye, respectively. This simulation method has the advantage of enabling us to predict the preferred orientation of one molecule to a second one when bound to each other to form a stable complex [ 38 , 39 , 40 ]. Interestingly, this helps us to economize cost and time of research work.…”