Flexible Ligand Docking Research Articles

A Molecular Basis for Agonist and Antagonist Actions at GABAC Receptors

We modelled the N-terminal ligand-binding domain of the rho1 GABA(C) receptor based on the Lymnaea stagnalis acetylcholine-binding protein (L-AChBP) crystal structure using comparative modelling and validated using flexible docking guided by known mutagenesis studies. A range of known rho1 GABA(C) receptor ligands comprising seven full agonists, 10 partial agonists, 43 antagonists and 12 inactive molecules were used to evaluate and validate the models. Of the 50 models identified, six models that allowed flexible ligand docking in accordance with the experimental data were selected and used to study detailed receptor-ligand interactions. The most refined model to accommodate all known active ligands featured a cavity comprising of a volume of 488 A(3). A detailed analysis of the interaction between the rho1 GABA(C) receptor model and the docked ligands revealed possible H-bonds and cation-pi interactions between the different ligands and binding site residues. Based on quantum mechanical/molecular mechanical (QM/MM) calculations, the model showed distinctive conformations of loop C that provided a molecular basis for agonist and antagonist actions. Agonists elicit loop C closure, while a more open loop C was observed upon antagonist binding. The model differentiates the role for key residues known to be involved in either binding and/or gating.

Conformational Analysis of TMC114, a Novel HIV-1 Protease Inhibitor

TMC114, a potent novel HIV-1 protease inhibitor, remains active against a broad spectrum of mutant viruses. In order to bind to a variety of mutants, the compound needs to make strong, preferably backbone, interactions and have enough conformational flexibility to adapt to the changing geometry of the active site. The conformational analysis of TMC114 in the gas phase yielded 43 conformers in which five types of intramolecular H-bond interactions could be observed. All 43 conformers were subject to both rigid and flexible ligand docking in the wild-type and a triple mutant (L63P/V82T/I84V) of HIV-1 protease. The largest binding energy was calculated for the conformations that are close to the conformation observed in the X-ray complexes of TMC114 and HIV-1 protease.

Flexible protein–flexible ligand docking with disrupted velocity simulated annealing

By docking flexible balanol to a rigid model of protein kinase A (PKA), we found that a new simulated annealing protocol termed disrupted velocity simulated annealing (DIVE-SA) outperformed the replica-exchange method and the traditional simulated annealing method in identifying the correct docking pose. In this protocol, the atomic velocities were reassigned periodically to encourage the system to sample a large conformational space. We also found that scaling potential energy surface to reduce structural transition barriers could further facilitate docking. The DIVE-SA method was then evaluated on its ability to perform flexible ligand-flexible protein docking of three ligands (balanol, a balanol analog, and ATP) to PKA. To reduce computational time and to avoid possible unphysical structural changes resulting from the use of nonoptimal force fields, a soft restrain was applied to keep the root-mean-square-deviation (RMSD) between instantaneous protein structures and a chosen reference structure small. Because the restrain was applied to the overall RMSD rather than to individual atoms, a protein could still experience relatively large conformational changes during docking. To examine the impact of applying such a restrain on docking, we constructed two semi-flexible protein models by choosing two different crystal structures as reference. Both the balanol analog and ATP were able to dock to either one of these semi-flexible protein models. On the other hand, balanol could only dock well to one of them. Further analysis indicated that the restrain on the glycine-rich loop was too strong, preventing it to adjust its structure to accommodate balanol in the binding pocket of PKA. Removing the restrain on the glycine-rich loop resulted in much better docking poses. This finding demonstrates the important role that the flexibility of the glycine-rich loop play in accepting different ligands and should profitably not be restrained in molecular docking so that more diverse ligands can be studied.

FlexE ensemble docking approach to virtual screening for CDK2 inhibitors

In spite of a proven potential and effectiveness of FlexE in docking flexible ligands into an ensemble of protein structures, FlexE has rarely been successful in virtual screening situations. In this study, we constructed cyclin-dependent kinase 2 (CDK2) ensemble structures which have exactly the same backbone conformations as 1AQ1 but differ only at the side chain torsion angles of the key amino acid (Lys33, Phe80, Lys89 and Asp145) residues: the torsion angles observed in the 17 CDK2 crystal structures were adapted to represent conformational flexibility. FlexE ensemble docking protocol then completely samples the full conformational fields generated by combination of torsions of the four amino acids in the ATP binding site. Virtual screening for CDK2 inhibitors by using the FlexE ensemble docking of a database composed of 48,703 inactives and 82 actives showed significant enrichment factor (EF = 18.5), and successfully identified 71 actives among the top 132 ligands (53.8%) ranked by total energy scores. Moreover, total energy scoring followed by visual inspection filtered-off non-specific binders among the highly-ranked ligands to increase the ratio of actives-to-inactives to 71:13 at the top 5% of the virtual screening solutions.

GlamDock: Development and Validation of a New Docking Tool on Several Thousand Protein−Ligand Complexes

In this study, we present GlamDock, a new docking tool for flexible ligand docking. GlamDock (version 1.0) is based on a simple Monte Carlo with minimization procedure. The main features of the method are the energy function, which is a continuously differentiable empirical potential, and the definition of the search space, which combines internal coordinates for the conformation of the ligand, with a mapping-based description of the rigid body translation and rotation. First, we validate GlamDock on a standard benchmark, a set of 100 protein-ligand complexes, which allows comparative evaluation to existing docking tools. The results on this benchmark show that GlamDock is at least comparable in efficiency and accuracy to the best existing docking tools. The main focus of this work is the validation on the scPDB database of protein-ligand complexes. The size of this data set allows a thorough analysis of the dependencies of docking accuracy on features of the protein-ligand system. In particular, it allows a two-dimensional analysis of the results, which identifies a number of interesting dependencies that are generally lost or even misinterpreted in the one-dimensional approach. The overall result that GlamDock correctly predicts the complex structure in practically half of the cases in the scPDB is important not only for screening ligands against a particular protein but even more so for inverse screening, that is, the identification of the correct targets for a particular ligand.

Flexible Ligand Docking with Glide

Glide is a ligand docking program for predicting protein-ligand binding modes and ranking ligands via high-throughput virtual screening. Glide utilizes two different scoring functions, SP and XP GlideScore, to rank-order compounds. Three modes of sampling ligand conformational and positional degrees of freedom are available to determine the optimal ligand orientation relative to a rigid protein receptor geometry. This unit presents protocols for flexible ligand docking with Glide, optionally including ligand constraints or ligand molecular similarities.

FLIPDock: Docking flexible ligands into flexible receptors

Conformational changes of biological macromolecules when binding with ligands have long been observed and remain a challenge for automated docking methods. Here we present a novel protein-ligand docking software called FLIPDock (Flexible LIgand-Protein Docking) allowing the automated docking of flexible ligand molecules into active sites of flexible receptor molecules. In FLIPDock, conformational spaces of molecules are encoded using a data structure that we have developed recently called the Flexibility Tree (FT). While the FT can represent fully flexible ligands, it was initially designed as a hierarchical and multiresolution data structure for the selective encoding of conformational subspaces of large biological macromolecules. These conformational subspaces can be built to span a range of conformations important for the biological activity of a protein. A variety of motions can be combined, ranging from domains moving as rigid bodies or backbone atoms undergoing normal mode-based deformations, to side chains assuming rotameric conformations. In addition, these conformational subspaces are parameterized by a small number of variables which can be searched during the docking process, thus effectively modeling the conformational changes in a flexible receptor. FLIPDock searches the variables using genetic algorithm-based search techniques and evaluates putative docking complexes with a scoring function based on the AutoDock3.05 force-field. In this paper, we describe the concepts behind FLIPDock and the overall architecture of the program. We demonstrate FLIPDock's ability to solve docking problems in which the assumption of a rigid receptor previously prevented the successful docking of known ligands. In particular, we repeat an earlier cross docking experiment and demonstrate an increased success rate of 93.5%, compared to original 72% success rate achieved by AutoDock over the 400 cross-docking calculations. We also demonstrate FLIPDock's ability to handle conformational changes involving backbone motion by docking balanol to an adenosine-binding pocket of protein kinase A.

Prediction of HIV-1 entry inhibitors neomycin–arginine conjugates interaction with the CD4-gp120 binding site by molecular modeling and multistep docking procedure

Developing of multi-target HIV-1 entry inhibitors represents an important avenue of drug therapy. Two such inhibitors are hexa-arginine-neomycin-conjugate (NeoR6) and nona- d-arginine-neomycin-conjugate (Neo-r9). Our findings that NeoR6-resistant mutations appear in the gp120 constant regions; and NeoR6 is not CCR5 antagonist, but inhibits CXCR4 and CCR5 HIV-1 using isolates, led us to suggest that NeoR6 may inhibit HIV-1 entry by interfering with the CD4-gp120 binding. To support this notion, we constructed a homology model of unliganded HIV-1 IIIB gp120 and docked NeoR6 and Neo-r9 to it, using a multistep docking procedure: geometric–electrostatic docking by MolFit; flexible ligand docking by Autodock3 and final refinement of the obtained complexes by Discover3. Binding free energies were calculated by MM-PBSA methodology. The model predicts competitive inhibition of CD4-gp120 binding by NeoR6 and Neo-r9. We determined plausible binding sites between constructed CD4-bound gp120 trimer and homology modeled membranal CXCR4, and tested NeoR6 and Neo-r9 interfering with this interaction. These models support our notion that another mechanism of anti-HIV-1 activity of NeoR6 is inhibition of gp120-CXCR4 binding. These structural models and interaction of NeoR6 and Neo-r9 with gp120 and CXCR4 provide a powerful approach for structural based drug design for selective targeting of HIV-1 entry and/or for inhibition of other retroviruses with similar mechanism of entry.

Alternative to Consensus ScoringA New Approach Toward the Qualitative Combination of Docking Algorithms

Since the development of the first docking algorithm in the early 1980s a variety of different docking approaches and tools has been created in order to solve the docking problem. Subsequent studies have shown that the docking performance of most tools strongly depends on the considered target. Thus it is hard to choose the best algorithm in the situation at hand. The docking tools FlexX and AutoDock are among the most popular programs for docking flexible ligands into target proteins. Their analysis, comparison, and combination are the topics of this study. In contrast to standard consensus scoring techniques which integrate different scoring algorithms usually only by their rank, we focus on a more general approach. Our new combined docking workflow-AutoxX-unifies the interaction models of AutoDock and FlexX rather than combining the scores afterward which allows interpretability of the results. The performance of FlexX, AutoDock, and the combined algorithm AutoxX was evaluated on the basis of a test set of 204 structures from the Protein Data Bank (PDB). AutoDock and FlexX show a highly diverse redocking accuracy at the different complexes which assures again the usefulness of taking several docking algorithms into account. With the combined docking the number of complexes reproduced below an rmsd of 2.5 A could be raised by 10. AutoxX had a strong positive effect on several targets. The highest performance increase could be found when redocking 20 protein-ligand complexes of alpha-thrombin, plasmepsin, neuraminidase, and d-xylose isomerase. A decrease was found for gamma-chymotrypsin. The results show that--applied to the right target-AutoxX can improve the docking performance compared to AutoDock and FlexX alone.

Docking Ligands into Flexible and Solvated Macromolecules. 1. Development and Validation of FITTED 1.0

We report the development and validation of a novel suite of programs, FITTED 1.0, for the docking of flexible ligands into flexible proteins. This docking tool is unique in that it can deal with both the flexibility of macromolecules (side chains and main chains) and the presence of bridging water molecules while treating protein/ligand complexes as realistically dynamic systems. This software relies on a genetic algorithm to account for the flexibility of the two molecules as well as the location of bridging water molecules. In addition, FITTED 1.0 features a novel application of a switching function to retain or displace key water molecules from the protein-ligand complexes. Two independent modules, ProCESS and SMART, were developed to set up the proteins and the ligands prior to the docking stage. Validation of the accuracy of the software was achieved via the application of FITTED 1.0 to the docking of inhibitors of HIV-1 protease, thymidine kinase, trypsin, factor Xa, and MMP to their respective proteins.

Development and validation of a modular, extensible docking program: DOCK 5

We report on the development and validation of a new version of DOCK. The algorithm has been rewritten in a modular format, which allows for easy implementation of new scoring functions, sampling methods and analysis tools. We validated the sampling algorithm with a test set of 114 protein-ligand complexes. Using an optimized parameter set, we are able to reproduce the crystal ligand pose to within 2 A of the crystal structure for 79% of the test cases using our rigid ligand docking algorithm with an average run time of 1 min per complex and for 72% of the test cases using our flexible ligand docking algorithm with an average run time of 5 min per complex. Finally, we perform an analysis of the docking failures in the test set and determine that the sampling algorithm is generally sufficient for the binding pose prediction problem for up to 7 rotatable bonds; i.e. 99% of the rigid ligand docking cases and 95% of the flexible ligand docking cases are sampled successfully. We point out that success rates could be improved through more advanced modeling of the receptor prior to docking and through improvement of the force field parameters, particularly for structures containing metal-based cofactors.

Activity–structure correlations in divergent lectin evolution: fine specificity of chicken galectin CG-14 and computational analysis of flexible ligand docking for CG-14 and the closely related CG-16

Gene duplication and sequence divergence are driving forces toward establishing protein families. To examine how sequence changes affect carbohydrate specificity, the two closely related proto-type chicken galectins CG-14 and CG-16 were selected as models. Binding properties were analyzed using a highly sensitive solid-phase assay. We tested 56 free saccharides and 34 well-defined glycoproteins. The two galectins share preference for the II (Galbeta1-4GlcNAc) versus I (Galbeta1-3GlcNAc) version of beta-galactosides. A pronounced difference is found owing to the reactivity of CG-14 with histo-blood group ABH active oligosaccharides and A/B active glycoproteins. These experimental results prompted to determine activity-structure correlations by modeling. Computational analysis included consideration of the flexibility of binding partners and the presence of water molecules. It provided a comparative description of complete carbohydrate recognition domains, which had so far not been characterized in animal galectins. The structural models assigned II, I selectivity to a region downstream of the central Trp moiety. Docking revealed that the tetrasaccharides can be accommodated in their free-state low-energy conformations. CG-14's preference for A versus B epitopes could be attributed to a contact between His124 and the N-acetyl group of GalNAc. Regarding intergalectin comparison, the Ala53/Cys51 exchange affects the interaction potential of His54/His52. Close inspection of simulated dynamic interplay revealed reorientation of His124 at the site of the His124/Glu123 substitution, with potential impact on ligand dissociation. In summary, this study identifies activity differences and provides information on their relation to structural divergence, epitomizing the value of this combined approach beyond galectins.

Evaluation of monoquaternary pyridinium oximes potency to reactivate tabun-inhibited human acetylcholinesterase

Monoquaternary N-benzyl-4-hydroxyiminomethylpyridinium bromide (Py-4-H) and its analogous with diverse substituents introduced into the phenyl ring (Py-4-CH 3, Py-4-Br, Py-4-Cl and Py-4-NO 2) were synthesized in order to examine their potency as reactivators of tabun-inhibited human erythrocyte acetylcholinesterase (AChE; EC 3.1.1.7). Within 24 h, the reactivation of tabun-inhibited AChE reached 80% with Py-4-CH 3, Py-4-Br and Py-4-Cl, 40% with Py-4-NO 2, and 30% with Py-4-H. The overall reactivation rate constants were up to 5.0 min −1 M −1. All oximes inhibited human AChE reversibly, and the inhibition potency increased in the following order Py-4-Br < Py-4-Cl < Py-4-CH 3 < Py-4-H < Py-4-NO 2. Although oximes Py-4-H and Py-4-NO 2 did not show significant reactivation ability, these oximes might be of interest as pre-treatment drugs due to their high affinity for the native AChE. Docking studies were carried out to elucidate the differences in oximes potency. The orientations of all studied oximes in the active site of human AChE have been proposed by flexible ligand docking with AutoDock 3.0. Analyses of the obtained complexes revealed the presence of numerous hydrogen bonds and close contacts between the oximes and the residues in the active site. Final docked energies predicted correctly the relative order of the inhibition potency of compounds (except in the case of Py-4-CH 3) as well as the most probable orientation of the best reactivator, Py-4-Br, which can result in an attack on the phosphorus atom of the tabun-phosphorylated human AChE.

EHiTS: A new fast, exhaustive flexible ligand docking system

The flexible ligand docking problem is divided into two subproblems: pose/conformation search and scoring function. For successful virtual screening the search algorithm must be fast and able to find the optimal binding pose and conformation of the ligand. Statistical analysis of experimental data of bound ligand conformations is presented with conclusions about the sampling requirements for docking algorithms. eHiTS is an exhaustive flexible-docking method that systematically covers the part of the conformational and positional search space that avoids severe steric clashes, producing highly accurate docking poses at a speed practical for virtual high-throughput screening. The customizable scoring function of eHiTS combines novel terms (based on local surface point contact evaluation) with traditional empirical and statistical approaches. Validation results of eHiTS are presented and compared to three other docking software on a set of 91 PDB structures that are common to the validation sets published for the other programs.

MolDock: A New Technique for High-Accuracy Molecular Docking

In this article we introduce a molecular docking algorithm called MolDock. MolDock is based on a new heuristic search algorithm that combines differential evolution with a cavity prediction algorithm. The docking scoring function of MolDock is an extension of the piecewise linear potential (PLP) including new hydrogen bonding and electrostatic terms. To further improve docking accuracy, a re-ranking scoring function is introduced, which identifies the most promising docking solution from the solutions obtained by the docking algorithm. The docking accuracy of MolDock has been evaluated by docking flexible ligands to 77 protein targets. MolDock was able to identify the correct binding mode of 87% of the complexes. In comparison, the accuracy of Glide and Surflex is 82% and 75%, respectively. FlexX obtained 58% and GOLD 78% on subsets containing 76 and 55 cases, respectively.

A theory of mode of action of azolylalkylquinolines as DNA binding agents using automated flexible ligand docking

Azolylalkylquinolines (AAQs) are a family of quinolines with varying degrees of cytotoxic activity (comparable or moderately superior to adriamycin in some cases) developed in the past decade in our group where their exact mode of action is still unclear. In this study the most probable DNA binding mode of AAQs was investigated employing a novel flexible ligand docking approach by using AutoDock 3.0. Forty-nine AAQs with known experimental inhibitory activity were docked onto d(CGCAAATTTGCG) 2, d(CGATCG) 2 and d(CGCG) 2 oligonucleotides retrieved from the Protein Data Bank (PDB IDs: 102D, 1D12 and 1D32, respectively) as the representatives of the three plausible models of interactions between chemotherapeutic agents and DNA (groove binding, groove binding plus intercalation and bisintercalation, respectively). Good correlation ( r 2 = 0.64) between calculated binding energies and experimental inhibitory activities was obtained using groove binding plus intercalation model for phenyl-azolylalkylquinoline (PAAQ) series. Our findings show that the most probable mode of action of PAAQs as DNA binding agents is via intercalation of quinolinic moiety between CG base pairs with linker chain and azole moiety binding to the minor groove.

A Novel Search Engine for Virtual Screening of Very Large Databases

Virtual screening of large chemical databases using the structure of the receptor can be computationally very demanding. We present a novel strategy that combines exhaustive similarity searches directly in SMILES format with the docking of flexible ligands, whose 3D structure is generated on the fly from the SMILES representation. Our strategy makes use of the recently developed LINGO tools to extract implicit chemical information from SMILES strings and integrates LINGO similarities into a pseudo-evolutionary algorithm. The algorithm represents a combination of a fast target-independent similarity method with a slower but information richer target-focused method. A virtual search of FactorXa ligands provided 62% of the potential hits after docking only 6.5% of a database of nearly 1 million molecules. The set of solutions showed good diversity, indicating that the method shows good scaffold hopping capabilities.

PSI‐DOCK: Towards highly efficient and accurate flexible ligand docking

We have developed a new docking method, Pose-Sensitive Inclined (PSI)-DOCK, for flexible ligand docking. An improved SCORE function has been developed and used in PSI-DOCK for binding free energy evaluation. The improved SCORE function was able to reproduce the absolute binding free energies of a training set of 200 protein-ligand complexes with a correlation coefficient of 0.788 and a standard error of 8.13 kJ/mol. For ligand binding pose exploration, a unique searching strategy was designed in PSI-DOCK. In the first step, a tabu-enhanced genetic algorithm with a rapid shape-complementary scoring function is used to roughly explore and store potential binding poses of the ligand. Then, these predicted binding poses are optimized and compete against each other by using a genetic algorithm with the accurate SCORE function to determine the binding pose with the lowest docking energy. The PSI-DOCK 1.0 program is highly efficient in identifying the experimental binding pose. For a test dataset of 194 complexes, PSI-DOCK 1.0 achieved a 67% success rate (RMSD < 2.0 A) for only one run and a 74% success rate for 10 runs. PSI-DOCK can also predict the docking binding free energy with high accuracy. For a test set of 64 complexes, the correlation between the experimentally observed binding free energies and the docking binding free energies for 64 complexes is r = 0.777 with a standard deviation of 7.96 kJ/mol. Moreover, compared with other docking methods, PSI-DOCK 1.0 is extremely easy to use and requires minimum docking preparations. There is no requirement for the users to add hydrogen atoms to proteins because all protein hydrogen atoms and the flexibility of the terminal protein atoms are intrinsically taken into account in PSI-DOCK. There is also no requirement for the users to calculate partial atomic charges because PSI-DOCK does not calculate an electrostatic energy term. These features are not only convenient for the users but also help to avoid the influence of different preparation methods.

The β-Carboline Analog Mana-Hox Causes Mitotic Aberration by Interacting with DNA

Mana-Hox, an analog of beta-carbolines with anticancer activity, induces aberrant mitosis and delays mitotic exit. However, the cellular target is not known. In this study, we visualized the intracellular localization of Mana-Hox. Mana-Hox rapidly penetrated into cells (within 1 min) and concentrated on disorganized metaphase chromosomes after 13 hr of exposure. We demonstrated that Mana-Hox is a noncovalent DNA binder that can interact with DNA through intercalation and/or through minor groove binding. Furthermore, Mana-Hox also inhibits topoisomerase II relaxation activity in vitro, suggesting that Mana-Hox could perturb mitotic chromosome decatenation. Overall, Mana-Hox binding to DNA plays a critical role in the induction of aberrant mitosis and contributes to its anticancer activity.

A Method for Induced-Fit Docking, Scoring, and Ranking of Flexible Ligands. Application to Peptidic and Pseudopeptidic β-secretase (BACE 1) Inhibitors

Inhibition of beta-secretase (BACE 1) has recently been investigated as a promising therapeutic approach in the treatment of Alzheimer's disease, and a growing number of BACE 1 inhibitors and crystal structures of BACE 1/inhibitors complexes have been reported. We report herein a predictive computational method and its application to potential BACE 1 inhibitors. Using a training set of 50 known highly flexible inhibitors, we developed a docking method that accounts for the flexibility of both the protein and the inhibitors. Protein flexibility is accounted for using a specifically designed genetic algorithm. We next developed a scoring function consisting of force field evaluation of the inhibitor/protein interactions and two additional terms for hydrogen bonding and entropy change upon binding. Discarding three outliers from the training set, our protocol was found to perform well with an rmsd of 1.19 kcal/mol. Evaluation of the predictive power was next carried out by virtual screening of 80 synthetic compounds. The significant enrichment at the top of the ranking list in active compounds demonstrated the ability of the docking and scoring protocol to rank the compounds relative to their activities.