- Open Access
PharmDock: a pharmacophore-based docking program
© Hu and Lill; licensee Chemistry Central Ltd. 2014
- Received: 20 January 2014
- Accepted: 4 April 2014
- Published: 16 April 2014
Protein-based pharmacophore models are enriched with the information of potential interactions between ligands and the protein target. We have shown in a previous study that protein-based pharmacophore models can be applied for ligand pose prediction and pose ranking. In this publication, we present a new pharmacophore-based docking program PharmDock that combines pose sampling and ranking based on optimized protein-based pharmacophore models with local optimization using an empirical scoring function.
Tests of PharmDock on ligand pose prediction, binding affinity estimation, compound ranking and virtual screening yielded comparable or better performance to existing and widely used docking programs. The docking program comes with an easy-to-use GUI within PyMOL. Two features have been incorporated in the program suite that allow for user-defined guidance of the docking process based on previous experimental data. Docking with those features demonstrated superior performance compared to unbiased docking.
A protein pharmacophore-based docking program, PharmDock, has been made available with a PyMOL plugin. PharmDock and the PyMOL plugin are freely available from http://people.pharmacy.purdue.edu/~mlill/software/pharmdock.
- Protein pharmacophores
- Biased docking
- Constraint docking
- Confined docking
Pharmacophore models aim to comprise the features of ligand-protein interactions that are most crucial for binding and biological activity. These models are used for virtual screening (VS) to identify potential new actives or for generating ligand alignments for subsequent QSAR simulations. Pharmacophore models are typically derived from structural features common to biologically active ligands that are hypothesized to be important for biological activity [1–5]. Such ligand-based pharmacophore models are dependent on the chemical features present in the known actives. Physicochemical features that are absent in the particular set of actives, but are important for the binding of structurally different ligands, will likely be neglected in the pharmacophore model. Alternatively, the binding site of the target protein can be used to generate a protein-based pharmacophore model without the inclusion of ligand information. These protein-based pharmacophore models are advantageous because a priori knowledge of active ligands is not required and the models are not biased by the chemical space of previously identified actives. Several approaches [6–9] have been developed to derive protein-based pharmacophore models from ligand-free proteins and apply the pharmacophore models in virtual screening.
Protein-based pharmacophore models are, by definition, enriched with the information of potential interactions between ligands and the protein target. Consequently, a direct application of the protein-based pharmacophore models is to use them for ligand pose prediction and pose ranking. In a recent study, we have explored the potential of protein-based pharmacophore models in ligand pose prediction and ranking . We carefully optimized the pharmacophore-generation process to reproduce native contacts for a large number of experimentally-determined protein-ligand complexes. We then developed a fast pharmacophore-based matching and scoring scheme and tested it on the PDBbind  core set. When the native ligand conformations, i.e. the ligand conformations from the crystal structures, were used as input, our pharmacophore-based scheme was able to identify native-like poses (with RMSD to the X-ray pose ≤ 2 Å) within the top-100 ranked poses for 94% of the tested entries. When the low-energy conformations generated by OpenEye Omega [12–14] were used as input, we were still able to maintain a success rate of 71% for predicting native-like binding pose within the top-100 ranked poses. These results were comparable in quality to several widely used docking programs.
Inspired by the promising results in ligand pose prediction and pose ranking using protein-based pharmacophore models, we extended the pharmacophore-based matching and scoring scheme into a docking program, named PharmDock (Pharmacophore-based Docking). The docking program further optimizes the top ranked binding poses predicted from the pharmacophore-based scheme and re-scores the optimized binding poses with a widely used empirical scoring function. We report here PharmDock’s performance in binding pose prediction and free energy of binding estimation tested on the PDBbind core set [11, 15], as well as its performance in virtual screening on 29 targets from the dictionary of useful decoys (DUD) dataset . We also present an open-source graphical user interface (GUI) adapted to PyMOL [17, 18] that we have developed for PharmDock for ease use of the docking software by the scientific community. In addition, we developed two new features within the PyMOL GUI allowing the users to guide the docking process towards specific residues identified from previous experimental data.
The docking program PharmDock
Generation of protein-based pharmacophores
Protein-based pharmacophores refer to the potential interaction sites for the ligand to favorably interact with the protein atoms in the binding site. They can be viewed as the negative or complementary image of the topology and physicochemical properties of the protein binding site. Four types of protein-based pharmacophores are defined in our program: hydrogen-bond donor/acceptor, hydrophobic, aromatic and ionic pharmacophores. In addition, the exclusive volume of the protein is also represented by the so-called forbidden pharmacophores, representing the binding site residues that would sterically overlap with ligand atoms placed at this respective position.
The detailed method in generating protein-based pharmacophores has been described in our previous paper . Briefly, the binding site of the protein is defined based on its known active ligand with a minimum of 3 Å to any ligand heavy atom. A 3D grid with 0.4 Å spacing between grid points was placed in the binding site for each protein structure. The interaction potentials (hydrogen-bond donor/acceptor, hydrophobic, aromatic and ionic) between the protein atoms and probes representing hypothetical ligand atoms were computed on each grid point. The interaction potentials for hydrogen-bonding and hydrophobic probes placed at the grid points were computed using a continuous form of the ChemScore [19, 20] scoring function. The aromatic and ionic interactions were calculated using a functional form similar to ChemScore. The detailed equations can be found in the previous publication . The pharmacophores were generated using the computed interaction energies with the probes on the 3D grid points. The hydrophobic pharmacophores were computed by a k-means clustering over all grid points with favorable hydrophobic scores. For each cluster, the hydrophobic pharmacophore element was then defined as the energy-weighted geometric center of all grid points of a particular cluster. The number of clusters, k, was adjusted until the minimum distance between a cluster center i and any other cluster center was on average smaller than a certain distance cutoff. K-means clustering to generate hydrogen-bond, aromatic and ionic pharmacophores was performed over the grid points associated with the same nearest functional group. For example, in generating a hydrogen-bond donor pharmacophore, the program iterates through all protein acceptors, and groups the grid points closest to the same acceptor into one patch. K-means clustering was then performed within this patch. In our previous study , we have investigated the influence of clustering distance cutoff of each pharmacophore type on the ligand pose sampling accuracy and efficiency. We found that pharmacophore models comprised by only hydrophobic and hydrogen bond elements, which are generated using a distance cutoff of 1.5 Å and 2.0 Å respectively, provide the best compromise between pose sampling accuracy and efficiency. These values will be used for the pharmacophore-based pose sampling process described below. For the pose-ranking process, a more detailed pharmacophore model using a 1 Å cluster distance cutoff for all pharmacophore types was adopted. The rationale is that the densest pharmacophore model provides the best description of the potential protein-ligand interactions and consequently should provide the largest amount of information for scoring.
Generation of ligand conformation and pharmacophores
PharmDock uses the low-energy conformers for a ligand generated by Openeye Omega [12–14] as docking input. For each ligand, a maximum of 100 conformations are generated with the calculated internal energy no more than 15 kcal/mol above the energy of the ligand conformation with the lowest internal energy. Duplicate conformers are removed using a 0.2 Å root-mean-square deviation (RMSD) cutoff for ligands with zero to three rotatable bonds, a 0.3 Å cutoff for ligands with four to six rotatable bonds, and a 0.4 Å cutoff for all ligands with more than six rotatable bonds. The in-house program clusterconformer is then used to generate the pharmacophore elements for each ligand conformation. Four types of pharmacophores are defined for each ligand: hydrogen-bond donor/acceptor, hydrophobic, aromatic and ionic pharmacophores. Hydrogen-bond pharmacophores are placed at the position of potential donor and acceptor groups of the ligand: Hydrogen-bond donors are polar hydrogen atoms bonded to oxygen, nitrogen and sulfur atoms, acceptors are oxygen, nitrogen and sulfur atoms with at least one lone pair. Ligand atoms (excluding hydrogen atoms) are defined to be hydrophobic if they were not hydrogen-bond donors or acceptors or directly bonded to a ligand’s donor or acceptor atoms. The hydrophobic atoms from each ligand conformation are clustered using hierarchical clustering with a minimum distance between cluster centers of 2.0 Å. Clustering is performed to reduce the number of hydrophobic ligand pharmacophores. This significantly reduces the cost of clique detection and consequently increases the efficiency of the docking process. Aromatic pharmacophores are defined as centers of aromatic rings. Ionic groups included functional groups that are formally charged positive or negative, e.g. protonated amines or deprotonated carboxylic acids, and are placed at the centroid of the functional group.
Pharmacophore-based pose sampling and ranking
The binding pose sampling and ranking process of PharmDock has been described and discussed in our previous publication . To provide the best compromise between accuracy and efficiency, only hydrophobic and hydrogen bond pharmacophore elements were used for the pose sampling process. This can be substantiated by our observation that on average only one aromatic interaction and less than one ionic interaction per protein - ligand complex are present in the 190 protein - ligand complexes we examined . This is in contrast to an average of four H-bond and ten hydrophobic interactions that were observed in the same dataset. A detailed evaluation of the pose sampling process can be found in our previous paper . Briefly, the pose sampling is based on a modified Bron-Kerbosch clique detection algorithm [21, 22] that enumerates all possible multi-points (> = 3) matches of ligand and protein-based pharmacophores. First, the length of the edge between each pair of ligand pharmacophores is determined. The edge length is also determined for each protein-based pharmacophore pair. All ligand pharmacophore edges that match the protein-based pharmacophore edges, based on the pharmacophore types (hydrogen bond donor/acceptor and hydrophobic) of their vertices and edge lengths, are identified. Throughout the matching process, a tolerance of 0.3 Å for the edge lengths is allowed. The matching process can be represented by a graph in which each node represents a matching ligand-protein pharmacophore pair. The clique detection algorithm then identifies all the completely connected subgraphs from this graph. The Kabsch algorithm  is then used to spatially align the ligand pharmacophore elements to the matching protein-based pharmacophores in each clique, thus placing the ligand into the protein binding site. To avoid steric clashes between ligand and protein atoms, the number of heavy atoms of the ligand that are located within 1.3 Å to any of the forbidden pharmacophores is counted for each ligand pose. If more than 10% of the ligand’s heavy atoms overlapped with forbidden pharmacophores, the pose is rejected.
r is the distance between the ligand pharmacophore and its closest matching protein-based pharmacophore of the same type. It is noteworthy that equation 2 calculates the score of a pose based on all the ligand pharmacophores rather than only those involved in forming the matching cliques. Also, the detailed pharmacophore models with 1 Å clustering distance cutoff were used for the ranking process.
Local optimization of the ligand binding poses
The pharmacophore-based pose sampling and ranking scheme was deployed to efficiently filter out ligand binding poses that are unlikely to be the native binding pose. As shown in our previous study , this procedure was indeed very effective to enrich the ranking list of binding poses with native-like poses within the top-100 positions. To further optimize the predicted ligand binding poses and estimate the binding energy of the ligand to the target, we perform a local Metropolis Monte Carlo (MC)  optimization for the top-100 ranked poses with the Piecewise Linear Pairwise (PLP) scoring function .
A detailed description of the PLP function can be found in Gehlhaar et al.’s  original work. Briefly, each ligand and protein atom is categorized into four types: H-bond donor, H-bond acceptor, H-bond donor/acceptor, nonpolar. Each pair of ligand-protein atoms is then assigned with one and only one of the interaction types: H-bond or steric. Both EH - bond and E steric are computed using similar piecewise linear functions but with different parameters . Despite the simplicity of PLP, it was shown to be one of the best scoring functions in identifying the true ligand binding poses among decoy poses and ranking the different ligands that bind to the same protein according to their binding affinities . Therefore, we finally settled on PLP scoring function for our local optimization and final scoring of the ligand binding poses.
where R is the ideal gas constant (1.986 cal/mol-K) and T the temperature (300 K). If this factor is larger than a random number generated from a uniform distribution between 0 and 1, the current pose is accepted.
Where ϵ is set to 0.15 and N max is 8. N atoms is the total number of heavy atoms attached to the smaller branch of this torsion. Equation 5 determines a scaling factor ‘ratio’ with which the maximum allowed torsional change is modified. For example, if the maximum torsional change is set to max = 60°, this value will be scaled by ratio = 0.5 + 0.5ϵ = 0.575 (yielding max = 34.5°) and ratio = 0.15 (max = 9°) for torsions with 4 and 8 heavy atoms attached to the smaller branch of the torsion, respectively.
Development of the PyMOL plugin
The docking calculations are performed in the background after submitting the job. The location of the docking results is stored in a monitoring file that allows the user to check the progress of the docking runs and re-import the docking results into PyMOL. All docking poses will be automatically displayed to the user with a separate dialog displaying the docking solutions and associated docking scores.
Tests of PharmDock’s docking performance
The “core set” of the PDBbind [11, 15] database (version 2007) was used to test PharmDock for its performance on ligand binding pose prediction and binding energy estimation. The PDBbind “core set” provides 210 protein-ligand complexes non-redundantly sampled from 1,300 protein-ligand complexes . It covers 70 different proteins, each of which contains three protein-ligand complexes with different binding affinities. All the protein-ligand complexes in the PDBbind core set were pre-processed with hydrogen atoms added and were therefore used directly without additional preparations. Due to various reasons, 20 protein-ligand complexes were excluded from the pose prediction and ranking study as described in our previous study .
The performance of PharmDock in ligand binding pose prediction was evaluated by its ability to reproduce the native-like pose for each protein-ligand complex in the PDBbind core set at the top-1 position. The top pose RMSD, i.e. the RMSD between the top-1 ranked pose to the native binding pose denoted as RMSDtop, was calculated. The average RMSDtop over the studied protein-ligand complexes was reported to assess the overall sampling performance. In addition, the percentages of complexes that were successfully predicted with RMSDtop within 1 Å, 2 Å and 3 Å to the native conformation were also used to evaluate the overall performance of pose generation and ranking. The correlation between PharmDock’s predicted binding energies with the experimentally measured binding affinities was used to evaluate PharmDock’s performance in binding energy estimation. And finally, the success rate of PharmDock in ranking three ligands bound to the same protein according to their binding affinities was also calculated.
The dictionary of useful decoys (DUD)  dataset was used to perform virtual screening (VS) studies. The DUD dataset contains 40 protein targets and a set of active and decoy ligands corresponding to each target. In the current version of PharmDock, the parameters of ions and cofactors were not included. Therefore, the four metalloenzymes, two folate enzymes and five other enzymes (aldose reductase, enoyl ACP reductase, glycogen phosphorylase β, purine nucleoside phosphorylase and S-adenosyl-homocysteine hydrolase) were excluded in our VS experiment. For each protein structure in DUD, the side-chain conformations of ASN, GLN and HIS, and tautomers and protonation states of HIS were adjusted using the Reduce program . The hydrogens were added to the protein using the tleap module of Amber 10 . The protein-based pharmacophores and ligand conformations and pharmacophores were generated for each target following the methods described in “Generation of Protein-based Pharmacophores” and “Generation of Ligand Conformation and Pharmacophores”.
To analyze the VS results, the ligands for each protein system were ranked based on their predicted binding energies. The Receiver Operating Characteristic (ROC) curve displaying the fraction of ranked actives (true positive rate) at a given fraction of ranked decoys (false positive rate) was plotted for each VS run. The area-under-the-curve (AUC) was calculated for each ROC curve and used to assess the overall enrichment quality.
Prediction of binding poses
A detailed analysis of the distribution of the RMSD values between the top-1 ranked binding pose and the X-ray binding pose (RMSDtop) for PharmDock with MC optimization is shown in Figure 6B. When the native ligand conformations were provided together with the low-energy conformations, PharmDock was able to predict binding poses with RMSDtop ≤ 1 Å for 45% of the tested protein-ligand complexes. When the native ligand conformations were excluded from the docking inputs, this number dropped to 20%. Obviously, the starting ligand conformations are critical for obtaining the correct docking solutions in PharmDock. As demonstrated in our previous study , Omega was not always able to sample ligand conformations within 1 Å RMSD to the native ligand conformation. Under our current setting (see “Generation of Ligand Conformation and Pharmacophore” in Methods section), 67% ligands had at least one conformer generated within 1 Å RMSD to its native conformation. This lack of generating native ligand conformations with Omega for a significant fraction of the docked ligands was the major reason for the dramatic drop in the success rate after excluding the native conformation from the docking inputs.
Comparison of PharmDock with docking programs evaluated by Plewczynski et al .
Native conformation as docking input1
Omega conformations as docking input2
%entries ≤ 2 Å
%entries ≤ 2 Å
Comparison of PharmDock with docking programs evaluated by Li et al .
Native conformation as docking input1
Low energy conformations as docking input2
Prediction of binding affinities
The PDBbind “core set” contains three protein-ligand complexes for each protein target. We studied if PharmDock is able to correctly rank the three different ligands bound to the same protein according to their binding affinities. There are 70 different proteins in the original PDBbind core set. After removal of the 20 complexes, as described before, 55 proteins remained that still had three different complexes available for our ranking study. For 42% of these 55 protein targets in both native-seeded and omega-only cases PharmDock was able to correctly rank the three ligands according to their binding affinities. According to the results reported by Cheng et al., when the original crystal complex structures are used for estimating the binding energies, the highest success rate of ranking the ligands for the same protein target was 58.5%. Our 42% success rate ranks at the 9th position compared to the 16 scoring functions tested by Cheng et al.. However, we need to mention again, that the top ranked docking poses rather than the X-ray conformations were used in our case, where the former is more complicated than the latter approach, as additional uncertainties are introduced throughout the pose sampling stage. Therefore a firm conclusion cannot be drawn directly from this comparison. However, combined with the comparable high correlation coefficients between the PharmDock predicted binding energies and the experimentally measured binding affinities, this result suggests that a good prediction in the overall correlation does not guarantee a good ranking power for the ligands bound with the same protein. A scoring function specifically developed for ranking the ligands bound to the same protein might be necessary for improving the performance of PharmDock in such studies.
Virtual screening experiment
Confined and constraint docking with PharmDock
To guarantee easy use of PharmDock, a PyMOL GUI was developed as described in the “Implementation” section. Besides the basic docking settings, the plugin also provides the user with the option to restrict docking to selections of the protein-based pharmacophores in order to concentrate the docking search volume to certain regions of the binding site. We termed this docking protocol as “confined docking”. Alternatively, the user can choose to utilize the “constraint docking” option, where at least one ligand pharmacophore of an accepted docking pose has to match with the user-selected pharmacophore elements. This option allows to enforce certain interactions that are critical for ligand binding to the target protein of interest. Typically, all pharmacophore elements of a property such as hydrogen-bonding are treated equally in the scoring function. The physicochemical environment of the hydrogen-bonding group or desolvation effects, however, can significantly alter the strength of a hydrogen bond . Pharmacophores derived from an analysis of existing protein–ligand structures or previously identified active compounds can be used to differentiate between important and weak hydrogen bonds. The pharmacophores can be incorporated as constraints into docking, thus binding modes that do not match any of the selected pharmacophores will be discarded. To demonstrate how these two novel docking features can benefit the docking performance, examples of using the confined and constraint docking are shown in this section.
Starting from pre-generated ligand and protein-based pharmacophores, we extended our pharmacophore-based pose sampling and ranking into a docking program, named PharmDock. PharmDock’s performance in ligand binding pose prediction and binding energy estimation was tested using the PDBbind core set. We found that the presence of the native ligand conformation greatly influenced PharmDock’s performance in pose prediction: when native ligand conformation was provided together with the low energy conformations in docking, PharmDock was able to reproduce the native binding pose within an RMSDtop ≤ 2 Å for 56% of the protein-ligand complexes. This success rate dropped to 37% when the native ligand conformation is excluded from the docking inputs. Interestingly, the influence of the native ligand conformation was not that dramatic for binding energy estimation. Linear correlations between the predicted binding scores and the experimentally measured binding affinities were observed. The Pearson correlation coefficient of 0.580 and 0.567 were reached with and without the presence of native ligand conformations in the docking input.
We also measured PharmDock’s ability in ranking different ligands bound to the same protein according to their binding affinities. This test is directly related to the virtual screening study where the rankings rather than the absolute binding energies are critical to distinguish potentially active ligands from decoys. PharmDock was able to correctly rank the ligands bound to the same protein according to their binding affinities for 42% of the tested proteins. Subsequent test of PharmDock in virtual screening on 29 targets of DUD dataset yielded an average AUC of 0.61. Dependent on the size of protein binding site and ligand molecule, the required computation time needed to dock one ligand using PharmDock varies between 0.5 and 7 min for the DUD dataset on a single core of a 2.5 GHz Quad Core AMD2380 computer. The average time required for docking was about 2.7 min per ligand.
Many docking programs have been developed since the establishment of the computer-aided molecular design field. While we understand that due to many influencing factors a fair comparison between different docking programs is quite difficult , we chose to compare PharmDock with several widely used docking programs evaluated by other studies so we can better understand PharmDock’s weaknesses and strengths. Using a simple empirical scoring function, PLP, PharmDock’s performances in binding pose prediction and binding energy estimation were comparable to or better than many widely used docking programs. Whereas PharmDock provided a high Pearson correlation coefficient between the estimated binding energy and experimentally measured binding affinities compared to other programs, the overall binding affinity prediction is still not accurate enough for drug optimization purposes. This difficulty in predicting binding affinities is also reflected in its mediocre performances in ranking ligands bound to the same protein and of virtual screening experiments, although comparable to other available docking methods. A more sophisticated scoring function is still necessary for improving PharmDock’s performance for separating physically similar ligands with distinct biological affinity. Furthermore, the comparison between the Native-Seeded and Omega-Only results also indicates that a better sampling of the ligand conformations, in particular for large, flexible ligands would improve PharmDock’s performance for this class of compounds. Last but not least, it has also been shown that the inclusion of protein flexibility and dynamics using the “ligand model concept” recently developed in our group can significantly improve the enrichment of virtual screening experiment . Future version of PharmDock will be combined with this “ligand model concept” to take protein flexibility and dynamics into consideration.
To make PharmDock accessible to any researcher in the field of biological and medicinal chemistry, we developed an open access PyMOl GUI for PharmDock. Both PharmDock and the PyMOl GUI can be downloaded from http://people.pharmacy.purdue.edu/~mlill/software/pharmdock. Two new features, confined docking and constraint docking, were built into the GUI. The idea is to provide users with the flexibility to include their expert knowledge about the target protein into the docking simulations. Two examples of these features demonstrated their usefulness in binding pose prediction and virtual screening.
Project name: PharmDock
Project home page: http://people.pharmacy.purdue.edu/~mlill/software/pharmdock
Operating system(s): Linux
Programming language: C, Python
Other requirements: OpenEye Omega (version 2.2.0 or newer), Python (version 2.4 or newer), Scipy, Numpy, OpenBabel (with Python bindings), PyMOL (tested on version 1.1r1 and 1.5.0)
License: GNU GPL v3
Any restrictions to use by non-academics: None
BH developed PharmDock and PharmDock plugin, prepared the data sets, tested PharmDock, and drafted the manuscript. MAL participated in the development of PharmDock and PharmDock plugin, the discussion of the results, and helped to draft the manuscript. Both authors read and approved the final manuscript.
M.A.L. thanks the National Institutes of Health (GM092855) for funding the present research.
- Martin Y: Distance comparisons (DISCO): a new strategy for examining 3D structure-activity relationships. Classical and Three-Dimensional QSAR in Agrochemistry. Edited by: Hansch C, Fujita T. 1995, Washington, DC: American Chemical Society, 318-329.View ArticleGoogle Scholar
- Barnum D, Greene J, Smellie A, Sprague P: Identification of common functional configurations among molecules. J Chem Inf Comput Sci. 1996, 36: 563-571. 10.1021/ci950273r.View ArticleGoogle Scholar
- Dixon SL, Smondyrev AM, Knoll EH, Rao SN, Shaw DE, Friesner RA: PHASE: a new engine for pharmacophore perception, 3D QSAR model development, and 3D database screening: 1. Methodology and preliminary results. J Comput Aided Mol Des. 2006, 20: 647-671. 10.1007/s10822-006-9087-6.View ArticleGoogle Scholar
- Richmond NJ, Abrams CA, Wolohan PRN, Abrahamian E, Willett P, Clark RD: GALAHAD: 1. Pharmacophore identification by hypermolecular alignment of ligands in 3D. J Comput Aided Mol Des. 2006, 20: 567-587. 10.1007/s10822-006-9082-y.View ArticleGoogle Scholar
- Chen X, Rusinko A, Tropsha A, Young SS: Automated pharmacophore identification for large chemical data sets 1. J Chem Inf Comput Sci. 1999, 39: 887-896. 10.1021/ci990327n.View ArticleGoogle Scholar
- Kirchhoff PD, Brown R, Kahn S, Waldman M, Venkatachalam C: Application of structure‒based focusing to the estrogen receptor. J Comput Chem. 2001, 22: 993-1003. 10.1002/jcc.1060.View ArticleGoogle Scholar
- Barillari C, Marcou G, Rognan D: Hot-spots-guided receptor-based pharmacophores (HS-Pharm): a knowledge-based approach to identify ligand-anchoring atoms in protein cavities and prioritize structure-based pharmacophores. J Chem Inf Model. 2008, 48: 1396-1410. 10.1021/ci800064z.View ArticleGoogle Scholar
- Tintori C, Corradi V, Magnani M, Manetti F, Botta M: Targets looking for drugs: A multistep computational protocol for the development of structure-based pharmacophores and their applications for hit discovery. J Chem Inf Model. 2008, 48: 2166-2179. 10.1021/ci800105p.View ArticleGoogle Scholar
- Cross S, Cruciani G: Grid-derived structure-based 3D pharmacophores and their performance compared to docking. Drug Discov Today Technol. 2011, 7: e213-e219.View ArticleGoogle Scholar
- Hu B, Lill MA: Exploring the potential of protein-based pharmacophore models in ligand pose prediction and ranking. J Chem Inf Model. 2013, 53: 1179-1190. 10.1021/ci400143r.View ArticleGoogle Scholar
- Wang R, Fang X, Lu Y, Yang CY, Wang S: The PDBbind database: methodologies and updates. J Med Chem. 2005, 48: 4111-4119. 10.1021/jm048957q.View ArticleGoogle Scholar
- OpenEye Scientific Software. Santa Fe N, USA, http://www.eyesopen.com OMEGA: version 2.2.0
- Hawkins PCD, Skillman AG, Warren GL, Ellingson BA, Stahl MT: Conformer generation with OMEGA: algorithm and validation using high quality structures from the Protein Databank and Cambridge Structural Database. J Chem Inf Model. 2010, 50: 572-584. 10.1021/ci100031x.View ArticleGoogle Scholar
- Hawkins PCD, Nicholls A: Conformer generation with OMEGA: learning from the dataset and the analysis of failures. J Chem Inf Model. 2012, 52: 2919-2936. 10.1021/ci300314k.View ArticleGoogle Scholar
- Wang R, Fang X, Lu Y, Wang S: The PDBbind database: collection of binding affinities for protein-ligand complexes with known three-dimensional structures. J Med Chem. 2004, 47: 2977-2980. 10.1021/jm030580l.View ArticleGoogle Scholar
- Huang N, Shoichet BK, Irwin JJ: Benchmarking sets for molecular docking. J Med Chem. 2006, 49: 6789-6801. 10.1021/jm0608356.View ArticleGoogle Scholar
- DeLano WL: The PyMOL Molecular Graphics System. 2002, San Carlos CA, USA: DeLano ScientificGoogle Scholar
- The PyMOL Molecular Graphics System. Schrödinger: LLC, Version 1.1r1: http://www.pymol.org/citing,
- Eldridge MD, Murray CW, Auton TR, Paolini GV, Mee RP: Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes. J Comput Aided Mol Des. 1997, 11: 425-445. 10.1023/A:1007996124545.View ArticleGoogle Scholar
- Baxter CA, Murray CW, Clark DE, Westhead DR, Eldridge MD: Flexible docking using tabu search and an empirical estimate of binding affinity. Proteins: Struct, Funct, Bioinf. 1998, 33: 367-382. 10.1002/(SICI)1097-0134(19981115)33:3<367::AID-PROT6>3.0.CO;2-W.View ArticleGoogle Scholar
- Bron C, Kerbosch J: Algorithm 457: finding all cliques of an undirected graph. Commun ACM. 1973, 16: 575-577. 10.1145/362342.362367.View ArticleGoogle Scholar
- Harley ER: Graph algorithms for assembling integrated genome maps. 2003, University of TorontoGoogle Scholar
- Kabsch W: A solution for the best rotation to relate two sets of vectors. Acta Crystallogr A. 1976, 32: 922-923. 10.1107/S0567739476001873.View ArticleGoogle Scholar
- Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH, Teller E: Equation of state calculations by fast computing machines. J Chem Phys. 1953, 21: 1087-10.1063/1.1699114.View ArticleGoogle Scholar
- Gehlhaar DK, Verkhivker GM, Rejto PA, Sherman CJ, Fogel DB, Fogel LJ, Freer ST: Molecular recognition of the inhibitor AG-1343 by HIV-1 protease: conformationally flexible docking by evolutionary programming. Chem Biol. 1995, 2: 317-324. 10.1016/1074-5521(95)90050-0.View ArticleGoogle Scholar
- Cheng T, Li X, Li Y, Liu Z, Wang R: Comparative assessment of scoring functions on a diverse test set. J Chem Inf Model. 2009, 49: 1079-1093. 10.1021/ci9000053.View ArticleGoogle Scholar
- Word J, Lovell S, Richardson J, Richardson D: Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation1. J Mol Biol. 1999, 285: 1735-1747. 10.1006/jmbi.1998.2401.View ArticleGoogle Scholar
- Case DA, Darden TA, Cheatham TE, Simmerling CL, Wang J, Duke RE, Luo R, Crowley M, Walker RC, Zhang W, Merz KM, Wang B, Hayik S, Roitberg A, Seabra G, Kolossváry I, Wong KF, Paesani F, Vanicek J, Wu X, Brozell SR, Steinbrecher T, Gohlke H, Yang L, Tan C, Mongan J, Hornak V, Cui G, Mathews DH, Seetin MG, et al: AMBER 10. 2008, San Francisco: University of CaliforniaGoogle Scholar
- Li X, Li Y, Cheng T, Liu Z, Wang R: Evaluation of the performance of four molecular docking programs on a diverse set of protein-ligand complexes. J Comput Chem. 2010, 31: 2109-2125. 10.1002/jcc.21498.View ArticleGoogle Scholar
- Kellenberger E, Rodrigo J, Muller P, Rognan D: Comparative evaluation of eight docking tools for docking and virtual screening accuracy. Proteins. 2004, 57: 225-242. 10.1002/prot.20149.View ArticleGoogle Scholar
- Plewczynski D, Łaźniewski M, Augustyniak R, Ginalski K: Can we trust docking results? Evaluation of seven commonly used programs on PDBbind database. J Comput Chem. 2011, 32: 742-755. 10.1002/jcc.21643.View ArticleGoogle Scholar
- Cross JB, Thompson DC, Rai BK, Baber JC, Fan KY, Hu Y, Humblet C: Comparison of several molecular docking programs: pose prediction and virtual screening accuracy. J Chem Inf Model. 2009, 49: 1455-1474. 10.1021/ci900056c.View ArticleGoogle Scholar
- Abel R, Salam NK, Shelley J, Farid R, Friesner RA, Sherman W: Contribution of explicit solvent effects to the binding affinity of small‒molecule inhibitors in blood coagulation factor serine proteases. Chemmedchem. 2011, 6: 1049-1066. 10.1002/cmdc.201000533.View ArticleGoogle Scholar
- Smith R, Hubbard RE, Gschwend DA, Leach AR, Good AC: Analysis and optimization of structure-based virtual screening protocols: (3). New methods and old problems in scoring function design. J Mol Graph Model. 2003, 22: 41-53. 10.1016/S1093-3263(03)00125-6.View ArticleGoogle Scholar
- Presnell SR, Patil GS, Mura C, Jude KM, Conley JM, Bertrand JA, Kam C-M, Powers JC, Williams LD: Oxyanion-mediated inhibition of serine proteases. Biochemistry (Mosc). 1998, 37: 17068-17081. 10.1021/bi981636u.View ArticleGoogle Scholar
- Sherawat M, Kaur P, Perbandt M, Betzel C, Slusarchyk WA, Bisacchi GS, Chang C, Jacobson BL, Einspahr HM, Singh TP: Structure of the complex of trypsin with a highly potent synthetic inhibitor at 0.97 A resolution. Acta Crystallogr Sect D: Biol Crystallogr. 2007, 63: 500-507. 10.1107/S090744490700697X.View ArticleGoogle Scholar
- Cole JC, Murray CW, Nissink JWM, Taylor RD, Taylor R: Comparing protein–ligand docking programs is difficult. Proteins: Struct, Funct, Bioinf. 2005, 60: 325-332. 10.1002/prot.20497.View ArticleGoogle Scholar
- Xu M, Lill MA: Significant enhancement of docking sensitivity using implicit ligand sampling. J Chem Inf Model. 2011, 51: 693-706. 10.1021/ci100457t.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.