Fully automated flexible docking of ligands into flexible synthetic receptors using forward and inverse docking strategies

The prediction of the structure of host-guest complexes is one of the most challenging problems in supramolecular chemistry. Usual procedures for docking of ligands into receptors do not take full conformational freedom of the host molecule into account. We describe and apply a new docking approach which performs a conformational sampling of the host and then sequentially docks the ligand into all receptor conformers using the incremental construction technique of the FlexX software platform. The applicability of this approach is validated on a set of host-guest complexes with known crystal structure. Moreover, we demonstrate that due to the interchangeability of the roles of host and guest, the docking process can be inverted. In this inverse docking mode, the receptor molecule is docked around its ligand. For all investigated test cases, the predicted structures are in good agreement with the experiment for both normal (forward) and inverse docking. Since the ligand is often smaller than the receptor and, thus, its conformational space is more restricted, the inverse docking approach leads in most cases to considerable speed-up. By having the choice between two alternative docking directions, the application range of the method is significantly extended. Finally, an important result of this study is the suitability of the simple energy function used here for structure prediction of complexes in organic media.     

Flexible ligand docking using a genetic algorithm

Summary Two computational techniques have been developed to explore the orientational and conformational space of a flexible ligand within an enzyme. Both methods use the Genetic Algorithm (GA) to generate conformationally flexible ligands in conjunction with algorithms from the DOCK suite of programs to characterize the receptor site. The methods are applied to three enzyme-ligand complexes: dihydrofolate reductase-methotrexate, thymidylate synthase-phenolpthalein and HIV protease-thioketal haloperidol. Conformations and orientations close to the crystallographically determined structures are obtained, as well as alternative structures with low energy. The potential for the GA method to screen a database of compounds is also examined. A collection of ligands is evaluated simultaneously, rather than docking the ligands individually into the enzyme.Abbreviations GA genetic algorithm; dhfr, dihydrofolate reductase - mtx methotrexate - ts thymidylate synthase - fen phenolphalein - HIV human immune deficiency virus - hivp HIV protease - thk thioketal haloperidol     

Flexible ligand docking using a robust evolutionary algorithm

A flexible ligand docking protocol based on evolutionary algorithms is investigated. The proposed approach incorporates family competition and adaptive rules to integrate decreasing‐based mutations and self‐adaptive mutations to act as global and local search strategies, respectively. The method is applied to a dihydrofolate reductase enzyme with the anticancer drug methotrexate and two analogues of antibacterial drug trimethoprim. Conformations and orientations closed to the crystallographically determined structures are obtained, as well as alternative structures with low energy. Numerical results indicate that the new approach is very robust. The docked lowest‐energy structures have root‐mean‐square derivations ranging from 0.67 to 1.96 Å with respect to the corresponding crystal structures. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 988–998, 2000     

Comparison of two implementations of the incremental construction algorithm in flexible docking of thrombin inhibitors

A set of 32 known thrombin inhibitors representing different chemical classes has been used to evaluate the performance of two implementations of incremental construction algorithms for flexible molecular docking: DOCK 4.0 and FlexX 1.5. Both docking tools are able to dock 10–35% of our test set within 2 Å of their known, bound conformations using default sampling and scoring parameters. Although flexible docking with DOCK or FlexX is not able to reconstruct all native complexes, it does offer a significant improvement over rigid body docking of single, rule-based conformations, which is still often used for docking of large databases. Docking of sets of multiple conformers of each inhibitor, obtained with a novel protocol for diverse conformer generation and selection, yielded results comparable to those obtained by flexible docking. Chemical scoring, which is an empirically modified force field scoring method implemented in DOCK 4.0, outperforms both interaction energy scoring by DOCK and the Böhm scoring function used by FlexX in rigid and flexible docking of thrombin inhibitors. Our results indicate that for reliable docking of flexible ligands the selection of anchor fragments, conformational sampling and currently available scoring methods still require improvement.     

Flexible docking in solution using metadynamics

We apply our recently developed metadynamics method to the docking of ligands on flexible receptors in water solution. This method mimics the real dynamics of a ligand exiting or entering an enzyme and in so doing reconstructs the free energy surface. We apply it to four docking cases: beta-trypsin/benzamidine, beta-trypsin/chlorobenzamidine, immunoglobulin McPC-603/phosphocholine, and cyclin-dependent kinase 2/staurosporine. In every case studied, the method is able to predict the docked geometry and the free energy of docking. Its added value with respect to many other available methods is that it reconstructs the complete free energy surface, including all the relevant minima and the barriers between them.     

Conformational searching using a population-based incremental learning algorithm

A new population-based incremental learning algorithm for conformational searching of molecules is presented. This algorithm is particularly effective at determining, by relatively small number of energy minimizations, global energy minima of large flexible molecules. The algorithm is also able to find a large set of low energy conformations of more rigid small molecules. The performance of the algorithm is relation to other algorithm is examined via the test molecules: C(18) H(38) , C(39)H(80) , cycloheptadecane and a set of five drug-like molecules.     

Complexing photolabile cholinergic ligands with synthetic and biological receptors: a dynamic survey

Molecular recognition processes between synthetic or biological receptors and small photolabile cholinergic ligands are described for adaptive host–guest complexes with calixarene derivatives and for the photoregulation of cholinergic enzymes, respectively.     

On-virus construction of polyvalent glycan ligands for cell-surface receptors

Glycans arrayed on the exterior of virus particles were used as substrates for glycosyltransferase reactions to build di- and trisaccharides from the virus surface. The resulting particles exhibited tight and specific associations with cognate receptors on beads and cells, in one example defeating in cis cell-surface interactions in a manner characteristic of polyvalent binding. Combined with the ability of viruses to provide structurally well-defined attachment points, the methodology provides a convenient and powerful way to prepare complex carbohydrate ligands for clustered receptors.     

An anchor-dependent molecular docking process for docking small flexible molecules into rigid protein receptors

A molecular docking method designated as ADDock, anchor- dependent molecular docking process for docking small flexible molecules into rigid protein receptors, is presented in this article. ADDock makes the bond connection lists for atoms based on anchors chosen for building molecular structures for docking small flexible molecules or ligands into rigid active sites of protein receptors. ADDock employs an extended version of piecewise linear potential for scoring the docked structures. Since no translational motion for small molecules is implemented during the docking process, ADDock searches the best docking result by systematically changing the anchors chosen, which are usually the single-edge connected nodes or terminal hydrogen atoms of ligands. ADDock takes intact ligand structures generated during the docking process for computing the docked scores; therefore, no energy minimization is required in the evaluation phase of docking. The docking accuracy by ADDock for 92 receptor-ligand complexes docked is 91.3%. All these complexes have been docked by other groups using other docking methods. The receptor-ligand steric interaction energies computed by ADDock for some sets of active and inactive compounds selected and docked into the same receptor active sites are apparently separated. These results show that based on the steric interaction energies computed between the docked structures and receptor active sites, ADDock is able to separate active from inactive compounds for both being docked into the same receptor.     

Flexible protein‐ligand docking using the Fleksy protocol

Considering protein plasticity is important in accurately predicting the three‐dimensional geometry of protein‐ligand complexes. Here, we present the first public release of our flexible docking tool Fleksy, which is able to consider both ligand and protein flexibility in the docking process. We describe the workflow and different features of the software and present its performance on two cross‐docking benchmark datasets. © 2012 Wiley Periodicals, Inc.     

ConCept: de novo design of synthetic receptors for targeted ligands

Low-molecular-weight receptors that bind targeted guest molecules have a wide range of potential applications but are difficult to design. This paper describes an evolutionary method for computer-aided design of such receptors that works by linking together chemical components from a user-defined library around a stable conformation of the targeted ligand. The software can operate in three modes: de novo design, in which it builds a wide variety of receptors from small components; macrocycle design, in which it builds homopolymeric macrocycles around the ligand; and elaboration of an existing receptor structure. The top candidates generated by the automatic construction process are further studied with detailed affinity calculations whose validity is supported by prior studies of experimentally characterized host-guest systems. All three modes of operation are illustrated here through the design of novel adenine receptors.     

Protein farnesyltransferase: Flexible docking studies on inhibitors using computational modeling

Summary Using MacroModel, peptide, peptidomimetic and non-peptidomimetic inhibitors of the zinc metalloenzyme, farnesyltransferase (FTase), were docked into the enzyme binding site. Inhibitor flexibility, farnesyl pyrophosphate substrate flexibility, and partial protein flexibility were taken into account in these docking studies. In addition to CVFM and CVIM, as well as our own inhibitors FTI-276 and FTI-2148, we have docked other farnesyltransferase inhibitors (FTIs) including Zarnestra, which presently is in advanced clinical trials. The AMBER* force field was employed, augmented with parameters that were derived for zinc. A single binding site model that was derived from the crystal structure of CVFM complexed with farnesyltransferase and farnesylpyrophosphate was used for these studies. The docking results using the lowest energy structure from the simulation, or one of the lowest energy structures, were generally in excellent agreement with the X-ray structures. One of the most important findings of this study is that numerous alternative conformations for the methionine side chain can be accommodated by the enzyme suggesting that the methionine pocket can tolerate groups larger than methionine at the C-terminus of the tetrapeptide and suggesting alternative locations for the placement of side chains that may improve potency.     

Multiple automatic base selection: Protein–ligand docking based on incremental construction without manual intervention

A possible way of tackling the molecular docking problem arising in computer- aided drug design is the use of the incremental construction method. This method consists of three steps: the selection of a part of a molecule, a so- called base fragment, the placement of the base fragment into the active site of a protein, and the subsequent reconstruction of the complete drug molecule. Assuming that a part of a drug molecule is known, which is specific enough to be a good base fragment, the method is proven to be successful for a large set of docking examples. In addition, it leads to the fastest algorithms for flexible docking published so far. In most real-world applications of docking, large sets of ligands have to be tested for affinity to a given protein. Thus, manual selection of a base fragment is not practical. On the other hand, the selection of a base fragment is critical in that only few selections lead to a low-energy structure. We overcome this limitation by selecting a representative set of base fragments instead of a single one. In this paper, we present a set of rules and algorithms to automate this selection. In addition, we extend the incremental construction method to deal with multiple fragmentations of the drug molecule. Our results show that with multiple automated base selection, the quality of the docking predictions is almost as good as with one manually preselected base fragment. In addition, the set of solutions is more diverse and alternative binding modes with low scores are found. Although the run time of the overall algorithm increases, the method remains fast enough to search through large ligand data sets.     

Biotin-specific synthetic receptors prepared using molecular imprinting

The composition of new molecularly imprinted polymers (MIPs) specific for biotin was optimised using molecular modelling software. Three functional monomers: methacrylic acid (MAA), 2-(trifluoromethyl)acrylic acid (TFAA) and 2-acrylamido-2-methylpropanesulfonic acid (AMPSA), which demonstrated the highest binding scores with biotin, were tested on their ability to generate specific binding sites. The imprinted polymers were photografted to the surface of polystyrene microspheres in water. The affinity of the synthetic “receptor” sites was evaluated in binding experiments using horseradish peroxidase-labelled biotin. Good correlation was found between the modelling results and the performance of the materials in the template re-binding study. The dissociation constants for all MIPs were 1.4-16.8 nM, which is sufficient for most analytical applications where biotin is used as a label.     

Massive docking of flexible ligands using environmental niches in parallelized genetic algorithms

Virtual screening of large libraries of small compounds requires fast and reliable automatic docking methods. In this article we present a parallel implementation of a genetic algorithm (GA) and the implementation of an enhanced genetic algorithm (EGA) with niching that lead to remarkable speedups compared to the original version AutoDock 3.0. The niching concept is introduced naturally by sharing genetic information between evolutions of subpopulations that run independently, each on one CPU. A unique set of additionally introduced search parameters that control this information flow has been obtained for drug‐like molecules based on the detailed study of three test cases of different complexity. The average docking time for one compound is of 8.6 s using eight R10,000 processors running at 200 MHz in an Origin 2000 computer. Different genetic algorithms with and without local search (LS) have been compared on an equal workload basis showing EGA/LS to be superior over all alternatives because it finds lower energy solutions faster and more often, particularly for high dimensionality problems. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1971–1982, 2001     

Flexible docking simulations: Scaled collective variable Monte Carlo minimization approach using Bezier splines,and comparison with a standard Monte Carlo algorithm

An algorithm for docking a flexible ligand onto a flexible or rigid receptor, using the scaled‐collective‐variables Monte Carlo with energy minimization approach, is presented. Energy minimization is shown to be one of the best techniques for distinguishing between native‐ and nonnative‐generated conformations. Incorporation of this technique into a Monte Carlo procedure enables one to distinguish the native conformation directly during the conformational search. It avoids the generation of a large number of ligand conformers for which more sophisticated energy evaluation tools would have had to be applied to identify the nativelike conformations. The efficiency of the Monte Carlo minimization was greatly improved by incorporating a new grid‐based energy evaluation technique using Bezier splines for which the energy function, as well as all of its derivatives, can be deduced from the values at the grid points. Comparison between our ECEPP/3‐based algorithm and the Monte Carlo algorithm presented elsewhere (Hart, T. N.; Read, R. J. Prot Struct Funct Genet 1992, 13, 206–222) has been made for docking NH₂ D Phe Pro Arg COOH, the noncovalent analog of NH₂ D Phe Pro Arg chloromethylketone (PPACK), onto the active site of human α‐thrombin. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 244–252, 1999     

Flexible protein‐protein docking based on Best‐First search algorithm

We developed a new high resolution protein‐protein docking method based on Best‐First search algorithm that loosely imitates protein‐protein associations. The method operates in two stages: first, we perform a rigid search on the unbound proteins. Second, we search alternately on rigid and flexible degrees of freedom starting from multiple configurations from the rigid search. Both stages use heuristics added to the energy function, which causes the proteins to rapidly approach each other and remain adjacent, while optimizing on the energy. The method deals with backbone flexibility explicitly by searching over ensembles of conformations generated before docking. We ran the rigid docking stage on 66 complexes and grouped the results into four classes according to evaluation criteria used in Critical Assessment of Predicted Interactions (CAPRI; “high,” “medium,” “acceptable,” and “incorrect”). Our method found medium binding conformations for 26% of the complexes and acceptable for additional 44% among the top 10 configurations. Considering all the configurations, we found medium binding conformations for 55% of the complexes and acceptable for additional 39% of the complexes. Introducing side‐chains flexibility in the second stage improves the best found binding conformation but harms the ranking. However, introducing side‐chains and backbone flexibility improve both the best found binding conformation and the best found conformation in the top 10. Our approach is a basis for incorporating multiple flexible motions into protein‐protein docking and is of interest even with the current use of a simple energy function. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Homology model-based virtual screening for GPCR ligands using docking and target-biased scoring

The current study investigates the combination of two recently reported techniques for the improvement of homology model-based virtual screening for G-protein coupled receptor (GPCR) ligands. First, ligand-supported homology modeling was used to generate receptor models that were in agreement with mutagenesis data and structure-activity relationship information of the ligands. Second, interaction patterns from known ligands to the receptor were applied for scoring and rank ordering compounds from a virtual library using ligand-receptor interaction fingerprint-based similarity (IFS). Our approach was evaluated in retrospective virtual screening experiments for antagonists of the metabotropic glutamate receptor (mGluR) subtype 5. The results of our approach were compared to the results obtained by conventional scoring functions (Dock-Score, PMF-Score, Gold-Score, ChemScore, and FlexX-Score). The IFS lead to significantly higher enrichment rates, relative to the competing scoring functions. Though using a target-biased scoring approach, the results were not biased toward the chemical classes of the reference structures. Our results indicate that the presented approach has the potential to serve as a general setup for successful structure-based GPCR virtual screening.     

Unsupervised guided docking of covalently bound ligands

Summary An approach for docking covalently bound ligands in protein enzymes or receptors was implemented in MacDOCK, a similarity-driven docking program based on DOCK 4.0. This approach was tested with a small number of covalent ligand–protein structures, using both native and non-native protein structures. In all cases, MacDOCK was able to generate orientations consistent with the known covalent binding mode of these complexes, with a performance similar to that of other docking programs. This method was also applied to search for known covalent thrombin inhibitors in a medium-sized molecular database (ca. 11,000 compounds). Detection of functional groups suitable for covalent docking was carried out automatically. A significant enrichment in known active molecules in the first 5% of the database was obtained, showing that MacDOCK can be used efficiently for the virtual screening of covalently bound ligands.     

FPDock: Protein–protein docking using flower pollination algorithm

Proteins play their vital role in biological systems through interaction and complex formation with other biological molecules. Indeed, abnormalities in the interaction patterns affect the proteins’ structure and have detrimental effects on living organisms. Research in structure prediction gains its gravity as the functions of proteins depend on their structures. Protein–protein docking is one of the computational methods devised to understand the interaction between proteins. Metaheuristic algorithms are promising to use owing to the hardness of the structure prediction problem. In this paper, a variant of the Flower Pollination Algorithm (FPA) is applied to get an accurate protein–protein complex structure. The algorithm begins execution from a randomly generated initial population, which gets flourished in different isolated islands, trying to find their local optimum. The abiotic and biotic pollination applied in different generations brings diversity and intensity to the solutions. Each round of pollination applies an energy-based scoring function whose value influences the choice to accept a new solution. Analysis of final predictions based on CAPRI quality criteria shows that the proposed method has a success rate of 58% in top10 ranks, which in comparison with other methods like SwarmDock, pyDock, ZDOCK is better. Source code of the work is available at: https://github.com/Sharon1989Sunny/_FPDock_.