A knowledge-guided strategy for improving the accuracy of scoring functions in binding affinity prediction
© Cheng et al; licensee BioMed Central Ltd. 2010
Received: 19 September 2009
Accepted: 17 April 2010
Published: 17 April 2010
Current scoring functions are not very successful in protein-ligand binding affinity prediction albeit their popularity in structure-based drug designs. Here, we propose a general knowledge-guided scoring (KGS) strategy to tackle this problem. Our KGS strategy computes the binding constant of a given protein-ligand complex based on the known binding constant of an appropriate reference complex. A good training set that includes a sufficient number of protein-ligand complexes with known binding data needs to be supplied for finding the reference complex. The reference complex is required to share a similar pattern of key protein-ligand interactions to that of the complex of interest. Thus, some uncertain factors in protein-ligand binding may cancel out, resulting in a more accurate prediction of absolute binding constants.
In our study, an automatic algorithm was developed for summarizing key protein-ligand interactions as a pharmacophore model and identifying the reference complex with a maximal similarity to the query complex. Our KGS strategy was evaluated in combination with two scoring functions (X-Score and PLP) on three test sets, containing 112 HIV protease complexes, 44 carbonic anhydrase complexes, and 73 trypsin complexes, respectively. Our results obtained on crystal structures as well as computer-generated docking poses indicated that application of the KGS strategy produced more accurate predictions especially when X-Score or PLP alone did not perform well.
Compared to other targeted scoring functions, our KGS strategy does not require any re-parameterization or modification on current scoring methods, and its application is not tied to certain systems. The effectiveness of our KGS strategy is in theory proportional to the ever-increasing knowledge of experimental protein-ligand binding data. Our KGS strategy may serve as a more practical remedy for current scoring functions to improve their accuracy in binding affinity prediction.
Molecular recognition plays an important role in many fundamental processes in biological systems [1, 2]. The basic concept of molecular recognition was first narrated by Emil Fischer more than 100 years ago. His "lock-and-key" theory , i.e. "... enzyme and glycoside must fit together like a key and a lock in order to initiate a chemical action upon each other...", has long been regarded as the basis for studying the binding between a ligand molecule to its biological receptor. Molecular docking, as a computational simulation of the ligand-receptor binding process, is widely applied in many research areas, such as structure-based drug design. For example, docking-based virtual screening [4–7] has become a complementary approach to high-throughput screening for the discovery of novel lead compounds, which is popular among academic research groups as well as pharmaceutical companies. In such a process, a library of small molecules are fit into the binding pocket of a given target protein through molecular docking, aiming at achieving an optimal complementarity of steric and physicochemical properties. Then, a computational method, which is often referred to as "scoring function", is used to evaluate the fitness between the ligand and the protein. All of the molecules are subsequently ranked by their binding scores, and only the most promising ones will be examined later in experiments. Obviously, a quantitative prediction of protein-ligand binding affinities is the key to the success of such studies.
As demonstrated in many previous studies [8, 9], today's molecular docking programs, such as DOCK , AutoDock [11–13], FlexX , Surflex [15, 16], LigandFit , GOLD [18, 19], and Glide [20, 21], are able to identify the correct binding pose of a flexible ligand to its receptor with a reasonable accuracy. However, binding affinity prediction (the "scoring problem") is still the Achilles' heel of molecular docking because in many cases the binding scores produced by scoring functions do not correlate well with true binding affinities, sometimes even cannot rank a set of compounds correctly [22–27]. The relatively disappointing performance of scoring functions in this aspect may be the result of a compromise between accuracy and speed since simplifications have to be made regarding solvation effects, conformational flexibility, and other factors in protein-ligand binding. A number of computationally more expensive methods have also been developed for binding affinity prediction in the past two decades or so. These methods typically conduct conformational sampling of the protein-ligand complex of interest through extensive molecular dynamics simulation in explicit solvent. Thus, they are able to address solvation effect and conformational flexibility in theory. Free energy pathway methods, such as free energy perturbation  and thermodynamics integration , sometimes can reproduce protein-ligand binding affinities within 1~2 kcal/mol. Nevertheless, they are normally used for the computation of the relative binding affinities of closely resembled ligand molecules, and thus have rather limited applications. "End-point" methods, such as the linear interaction energy (LIE) approximation [30, 31] and the MM-PB/SA method , avoid the integration of free energy pathway in order to save computation cost. They are normally applied to the modeling of a congeneric set of compounds binding to the same target protein. A number of successful applications of these methods have already been reported in literature. However, their success seems to rely on well-selected systems, and their robustness still needs to be validated more extensively. Inaccuracy in force field or inadequate sampling may account for their possible failures. In fact, some comparative tests  indicate that such methods are not necessarily more accurate than scoring functions although they definitely consume more computational resources. Technically, these methods are by far too computationally expensive for high-throughput tasks. It is also complicated to set up a job with such methods. Due to these concerns, such methods are not likely to be integrated into molecular docking programs for practical uses.
Considering the balance between accuracy, efficiency, and applicability, scoring functions are still the best choice to tackle the scoring problem for molecular docking and some other tasks in structure-based drug design. Thus, improving the general performance of scoring functions is undoubtedly a worthwhile aim. A good number of scoring functions have already been reported in literature since 1990s. They can be classified roughly into three categories: (i) Force field-based methods [10–13, 18, 19] rely on established force fields to compute the non-covalent interactions between protein and ligand, including van der Waals and electrostatic interactions. They are often augmented by GB/SA or PB/SA terms in order to consider solvation effect. (ii) Empirical scoring functions [20, 21, 34–43] decompose the protein-ligand binding free energy into some basic terms, such as hydrogen bonding, hydrophobic effect and so on. Each term is computed with an intuitive algorithm, and the weight factors of each term are typically derived from a regression analysis on a set of protein-ligand complexes with known binding affinities. Hence, empirical scoring functions are also referred to as regression-based methods. (iii) Methods based on potentials of mean force [44–52] compute protein-ligand interactions as a sum of distance-dependent pairwise potentials. A technical advantage of these methods is that deduction of potentials of mean force only requires the knowledge of protein-ligand complex structures.
Most of today's scoring functions are developed as all-purpose models, which are presumably applicable to all sorts of protein-ligand complexes. A recent comparative assessment of 16 popular scoring functions conducted by us  revealed that an accurate prediction of the binding affinities across a variety of protein-ligand complexes is still a major challenge for them. Nevertheless, it was also noticed that on certain classes of complexes, some scoring functions actually produced very promising results. This observation indicates that it is perhaps more practical to develop specific scoring functions applicable to certain classes of protein-ligand complexes in order to improve the accuracy in binding affinity prediction. In fact, this idea has been practiced by some researchers. For instance, Teramoto et al. introduced the supervised scoring modes as well as an optimized consensus scoring scheme with feature selection to enhance the enrichment factor in structure-based virtual screening [53, 54]. Recently, Seifert and co-workers published a review article on target-specific scoring functions (or "targeted scoring functions") . Current all-purpose scoring functions can be re-calibrated on certain classes of protein-ligand complexes to become target-specific, which is probably the most straightforward approach. For example, the DrugScore-RNA  adopts the same framework as DrugScore [45, 46], a scoring function based on potentials of mean force, but is derived from 670 crystal structures of nucleic acid-ligand and nucleic acid-protein complexes. Antes proposed the POEM approach (Parameter Optimization using Ensemble Methods)  and applied it to the optimization of two scoring functions (FlexX and ScreenScore) on kinases and ATPases. Seifert described a statistical method (ProPose) for improving the signal-to-noise ratio of scoring functions in molecular docking and successfully customized the Böhm scoring function on three selected proteins: cyclin dependent kinase 2, estrogen receptor, and cyclooxygenase 2 .
The targeted scoring functions mentioned above typically require re-parameterization or special customization on current scoring functions to become suitable for specific systems. In this study, we have proposed an alternative strategy, namely knowledge-guided scoring (KGS). This strategy requires no re-parameterization in prior, and in principle can be applied in combination with any scoring functions to any classes of protein-ligand complexes. The key idea is that the unknown binding affinity of a given protein-ligand complex can be estimated more reliably based on the known binding affinity of an appropriate reference complex. The reference complex is required to share a similar pattern of key protein-ligand interactions with the given protein-ligand complex. For this purpose, our KGS strategy utilizes a sufficient number of relevant protein-ligand complexes with known structures and binding affinities as a knowledge set. The key protein-ligand interactions in each complex are summarized as a pharmacophore model, which can be elucidated by an automatic algorithm implemented by us. The knowledge set is then searched through for the appropriate reference complex for any given protein-ligand complex. Thus, our KGS strategy can take full advantage of known knowledge, resulting in an improved accuracy in binding affinity prediction. In our study, the KGS strategy was tested in combination with two all-purpose scoring functions, i.e. X-Score  and PLP [37, 38] on three sets of protein-ligand complexes. An improved average accuracy for both X-Score and PLP was indeed obtained. Detailed descriptions are given in the following sections.
Equation 4 indicates how the binding affinity of a given protein-ligand complex is computed using the known binding affinity of a proper reference complex as a starting point. For the convenience of narration, this scoring strategy will be referred to as the KGS strategy, i.e. Knowledge-Guided Scoring, throughout this article. In principle, any scoring method may be employed to compute the required binding scores of both the reference complex and the query complex in Equation 4. Nevertheless, it is certainly more reasonable in reality to choose a capable scoring method for this purpose. The reference complex can be selected among a database of protein-ligand complexes with reliable structures and binding data. The constant k in Equation 4 can be derived through a regression analysis between the experimental binding data and the computed binding scores by the employed scoring method on the same database. It is introduced to scale the outcomes of scoring functions, which could be in arbitrary units, to a realistic range comparable to the experimental binding data of the reference complex.
The only input needed by our algorithm is the three-dimensional structure of a given protein-ligand complex, which can be an experimentally resolved structure or a computer model generated by molecular docking or other methods. In either case, the ligand molecule is required to be inside the binding site on the protein with the desired binding pose. The very first step of our algorithm is to identify the amino acid residues on the protein which are in direct contact with the bound ligand molecule. An amino acid residue is considered to be in direct contact with the ligand if any heavy atom on it is within a distance of 5.0 Å from any heavy atom on the ligand. This distance cutoff is adjustable to the users. A box large enough is then used to enclose all of these residues as well as the ligand molecule (Figure 1). Next, evenly spaced grids are created inside the box with a spacing of 0.5 Å by default. Each grid is checked for its accessibility by placing a hydrogen atom on it. If the hydrogen probe bumps with any protein atom, i.e. when the inter-atom distance of the two participating atoms is shorter than the sum of their van der Waals radii minus 0.5 Å, the grid under consideration will be labeled as "ignored". The van der Waals atomic radii used in our algorithm are cited from the Tripos force field. If a grid is more than 5.0 Å away from any atom on the protein, it will also be labeled as "ignored", as there will be no direct interaction between the protein and any ligand atom placed on this grid. All "ignored" grids will be removed later on to speed up the following processes. As result, the remaining grids actually define the binding pocket for subsequent analyses.
For each remaining grid, three different types of probes are placed on it and the binding score between each probe and the protein are evaluated thereby. These probes include (1) a positively charged sp3 nitrogen atom (ammonium cation), representing a hydrogen bond donor; (2) a negatively charged sp2 oxygen atom (as in a carboxyl group), representing a hydrogen bond acceptor; and (3) a sp3 carbon atom (methane), representing a hydrophobic group. The binding scores between each probe and the protein are calculated by using the corresponding algorithms in the empirical scoring function X-Score (version 1.2) . Final classification of each grid will be determined according to the particular probe which produces the highest binding score, either as "donor", "acceptor", or "hydrophobic". For example, a grid labeled as "donor" indicates that a hydrogen bond donor is mostly preferred on this particular grid. All of the scored grids are further refined and clustered by a two-step process. At the first step, the average score is calculated over all donor grids. The donor grids whose scores are lower than the average score are re-labeled as "ignored". The same process is also repeated on the acceptor grids and the hydrophobic grids. Consequently, only the grids with significant contributions to protein-ligand interactions will survive (Figure 1 & Figure 2B). At the second step, our algorithm checks each remaining donor grid and counts the total number of its "neighbors", i.e. the remaining grids of the same type within a range of 2.0 Å. The average number of neighbors for all donor grids is calculated. Those grids with a total number of neighbors below the average will be re-labeled as "ignored" and filtered out. The same process is also repeated on all acceptor grids and the hydrophobic grids. After this process, only the grids in aggregation will still survive, which represent the key interaction sites inside the binding pocket more clearly (Figure 1 & Figure 2C).
Here, PS i is the pharmacophore score of the i th grid under consideration, S j is the binding score of the j th neighboring grid of this grid, and r ij is the distance between the i th grid and the j th grid. According to this algorithm, the pharmacophore score of a certain grid combines the contributions from all its neighboring grids within 2.0 Å. The spatial distribution of the neighboring grids also has an impact: a dense group of grids will be associated with a higher pharmacophore score; whereas the pharmacophore score for a sparse group of grids will be somewhat lower. Besides, a correction is introduced in Equation 5 to favor the pharmacophore features overlapping with the atoms on the ligand: The d i in Equation 5 is the distance between the i th grid and a close non-hydrogen atom on the given ligand. Thus, if a pharmacophore feature overlaps exactly with a certain atom on the ligand, its pharmacophore score will retain by 100%; otherwise its pharmacophore score will receive a certain distance-dependent discount. In order to avoid the generation of too many pharmacophore features, our algorithm by default sets the minimal distance between two features to 3.5 Å, approximately the average van der Waals distance between two non-hydrogen atoms. The above process is repeated until all groups of grids are attributed to certain pharmacophore features. Occasionally, too many pharmacophore features are deduced by our algorithm when the binding pocket is really large. Thus, an upper limit of 15 features in a pharmacophore model is set for the sake of subsequent similarity searching. By using our program, the average computation time consumed by pharmacophore elucidation for a typical protein-ligand complex (with ~30 residues at the binding site) was merely 170 ms on a low-end laptop with one 2.13 GHz CPU and 2GB memory inside.
Here, N P and N q are the numbers of pharmacophore features in models P and q, respectively; while N Pq is the maximal number of common features found between P and q.
For any given protein-ligand complex, an external database of known protein-ligand complexes supplied by the user will be examined with the algorithm described above. The complex with the highest similarity index, if higher than a user-set cutoff, is chosen as the reference in Equation 4 for calculating the binding affinity of the query complex. By using our program, the similarity search is fairly fast, which took less than 1 ms to map a pair of pharmacophore models on a low-end laptop with one 2.13 GHz CPU and 2GB memory inside.
Preparation of test sets
The complex structure is resolved through X-ray crystal diffraction with an overall resolution better or equal to 2.5 angstroms. Neither the protein nor the ligand has any missing fragment in the crystal structure.
Binding data of the complex is experimentally measured as either dissociation constant (K d ) or inhibition constant (K i ). Both the protein and the ligand used in binding assay match exactly the ones used in structure determination.
The complex is formed by one protein molecule and one ligand molecule in a binary manner, and the binding is non-covalent in nature. The ligand molecule does not contain any uncommon elements, such as Be, B, Si and metal atoms, and its molecular weight does not exceed 1000.
PDB codes of the protein-ligand complexes in the three test sets
HIV protease complexes (N= 112)
1GNM, 1GNN, 1GNO, 1A30, 1A9M, 1AAQ, 1AJV, 1AJX, 1B6J, 1B6K, 1B6L, 1B6M, 1BDQ, 1BV7, 1BV9, 1BWA, 1BWB, 1C70, 1D4K, 1D4L, 1D4Y, 1DIF, 1DMP, 1G2K, 1G35, 1HBV, 1HEG, 1HIH, 1HII, 1HOS, 1HPO, 1HPS, 1HPV, 1HPX, 1HSH, 1HVH, 1HVI, 1HVJ, 1HVK, 1HVL, 1HVR, 1HVS, 1HWR, 1HXB, 1HXW, 1IIQ, 1IZH, 1IZI, 1LZQ, 1MES, 1MET, 1MEU, 1MRW, 1MRX, 1MSM, 1MSN, 1MTR, 1NH0, 1ODY, 1OHR, 1PRO, 1QBR, 1QBS, 1QBT, 1QBU, 1SBG, 1SDT, 1SDU, 1SDV, 1SGU, 1SH9, 1T7J, 1W5V, 1W5W, 1W5X, 1W5Y, 1Z1H, 1Z1R, 1ZP8, 1ZPA, 1ZSF, 1ZSR, 2AOC, 2AOD, 2AOE, 2AQU, 2AVM, 2AVO, 2AVQ, 2AVS, 2AVV, 2BPV, 2BPY, 2BQV, 2F80, 2F81, 2F8G, 2FGU, 2FGV, 2HB3, 2I0A, 2I0D, 7HVP, 7UPJ, 2HS2, 2AOG, 2HS1, 1A94, 1AID, 1KZK, 1TCX, 3AID
Trypsin complexes ( N = 73)
1C1R, 1C5P, 1C5Q, 1C5S, 1C5T, 1CE5, 1F0T, 1F0U, 1K1I, 1K1J, 1K1L, 1K1M, 1K1N, 1OSS, 1PPC, 1PPH, 1QB1, 1QB6, 1QB9, 1QBN, 1QBO, 1TNG, 1TNH, 1TNI, 1TNJ, 1TNK, 1TNL, 1V2J, 1V2K, 1V2L, 1V2N, 1V2Q, 1V2R, 1V2S, 1V2T, 1V2U, 1V2W, 2A31, 2BZA, 2FX6, 1BRA, 1G3B, 1G3C, 1G3D, 1G3E, 1GHZ, 1GI1, 1GI4, 1GI6, 1GJ6, 1J16, 1J17, 1O2H, 1O2J, 1O2K, 1O2N, 1O2O, 1O2Q, 1O2S, 1O2W, 1O2X, 1O2Z, 1O30, 1O33, 1O36, 1O38, 1O3D, 1O3F, 1O3H, 1O3I, 1O3J, 1O3K, 1V2O
Carbonic anhydrase complexes ( N = 44)
1BN1, 1BN3, 1BN4, 1BNN, 1BNQ, 1BNT, 1BNU, 1BNV, 1BNW, 1A42, 1AVN, 1BCD, 1CIL, 1CIM, 1CIN, 1CNW, 1CNX, 1CNY, 1G1D, 1G45, 1G46, 1G48, 1G4J, 1G4O, 1G52, 1G53, 1G54, 1I9L, 1I9M, 1I9N, 1I9O, 1I9P, 1I9Q, 1IF7, 1IF8, 1OKL, 1TTM, 1XPZ, 1XQ0, 1YDA, 1YDB, 1YDD, 2EZ7, 2H4N
In addition, a set of putative binding poses were prepared for the ligand molecule in each complex in all three test sets. These binding poses were generated with the GOLD software (version 4.1) by docking the native ligand into its co-crystallized protein target. The parameter "No. of GA operations" was set to 10000 to get docking poses as diverse as possible. The ChemScore scoring function implemented in GOLD was chosen as the scoring engine. All other parameters were assigned the default values. A total of 100 top-ranked docking poses were retained for each ligand. The root-mean-square deviations (RMSD) from the native ligand pose of each docking pose was calculated using the "rms_analysis" utility in GOLD.
Regarding the HIV protease test set, at least one docking pose with RMSD < 2 Å from the native binding pose are found among all GOLD-generated docking poses for 90 ligands out of the total 112. It indicates that near-native poses were successfully generated for most ligands with our method. The average RMSD of the docking poses for these 90 ligands was 5.2 Å, indicating that those docking poses were also diverse. Another 19 ligands had at least one docking pose with RMSD falling in the range of 2~3 Å. Our method failed to produce docking poses with RMSD < 3 Å only for three ligands in the HIV protease test set. These three cases were all associated with large and flexible ligand molecules. As for the carbonic anhydrase and the trypsin test sets, 39 out of 44 and 64 out of 73 ligands, respectively, had at least one docking pose with RMSD < 2 Å from the native ligand pose among all GOLD-generated docking poses. Finally, the corresponding native binding pose observed in crystal structure was added to the ensemble of GOLD-generated docking poses for each ligand to ensure that the most important point in the conformational space of each ligand was sampled.
Evaluation of scoring methods
Our recent study  revealed that two scoring functions, i.e. X-Score  and PLP [37, 38], have relatively better performance in binding affinity prediction than other scoring functions. They were thus chosen in this study to test the KGS strategy. Note that X-Score has three built-in options, i.e. HPScore, HMScore, and HSScore. In many cases, the difference in the outcomes of these three options is marginal. For the sake of convenience, only the average value of these three options in X-Score was considered in our study. Similarly, PLP also has two variations, i.e. PLP1 and PLP2. Only PLP1 was considered in our study. The X-Score program was obtained from its original authors. PLP1 was implemented by us according to the descriptions given in the original references [37, 38] as well as the information given in the user manual of the Discovery Studio software (version 2.0) . We compared the results produced by our in-house implementation of PLP1 and those produced by the one implemented in Discovery Studio on the entire PDBbind refined set. These two sets of results were found to be almost identical (data not shown), indicating that our own implementation of PLP1 was correct.
For each test set, both X-Score and PLP1 were used to compute the binding scores of all member protein-ligand complexes. For each given complex, the pharmacophore mapping algorithm described in a previous section was applied to identify a proper reference complex among the other complexes included in the same test set. If such a reference complex was found, an adjusted binding score of the given protein-ligand complex was computed with Equation 4 based on the known binding constant of the reference complex. If not, the binding score of the given protein-ligand complex was unchanged. This process was repeated until all protein-ligand complexes in each test set had been processed. Then, the Pearson correlation coefficient (R p ) between the experimentally determined binding constants and the final binding scores was calculated. The standard deviations (SD) in fitting the computed binding scores to the experimental binding constants were used as a quantitative measurement of the accuracy of each scoring method for comparison. The smaller is the standard deviation, the better is the accuracy.
The computations described above were all based on experimentally resolved protein-ligand complex structures, and the corresponding results are referred to as Set I throughout this article. In the reality of binding affinity prediction, however, one needs to rely on predicted structures in most cases. Thus, it is necessary to test our KGS strategy in such scenarios as well. For this purpose, a total number of 100 putative docking poses were prepared for each ligand in all three test sets. For each complex, all of these docking poses as well as the native binding pose were computed by X-Score and PLP, respectively, using the same procedure described in the previous paragraph. The binding score of each docking pose was adjusted using Equation 4 whenever applicable. Then, the best binding score obtained across all docking poses and the native binding pose for each protein-ligand complex was considered in the correlation analysis with the experimental data. In other words, the docking poses generated by the GOLD software were actually rescored by X-Score and PLP in combination with our KGS strategy. The corresponding results are referred to as Set II throughout this article.
Results and Discussion
Performance on HIV protease complexes
Statistical results produced by two scoring functions alone on the three test sets
HIV protease complexes
(N = 112)
Carbonic anhydrase complexes
(N = 44)
(N = 73)
Information on the five pairs of HIV protease complexes with the highest similarities
The query complex
The reference complex
-log K d (exp)
-log K d (exp)
The results produced by PLP in combination with the KGS strategy (Figure 4A) basically reveal the same trend: at lower similarity cutoffs, the improvements exhibited by PLP+KGS over PLP alone are not obvious; while at higher similarity cutoffs, considerable improvements are observed, demonstrating the success of the KGS strategy again. An exception is that PLP alone produced acceptable results on the five complexes listed in Table 3, but this can be well interpreted as a coincidence. Interestingly, considerable improvements are also observed above the similarity cutoff of 0.40, exactly the same as the case of X-Score+KGS. In addition, the corresponding results produced by X-Score+KGS and PLP+KGS are really close on all subsets of complexes when the similarity cutoff > 0.40 although the performance of two scoring functions alone can be very different. These findings further support our statement that the success of the KGS strategy is in principle independent from the scoring method employed in computation.
When computer-generated docking poses are considered in scoring for instead, i.e. Set II, the standard deviations produced by X-Score and PLP as a function of the cutoffs used in similarity search are plotted in Figure 4B. Detailed statistical results can be found in the Additional file 1, Table S2. The same trend as in Set I results are observed for both X-Score and PLP: the performance of both scoring functions is improved by the KGS strategy when the similarity cutoff used in reference searching is higher than 0.40. The standard deviations of both scoring functions are lowered by 0.2 units or even more. This observation indicates that our KGS strategy can be applied not only to the complexes with experimentally determined structures but also predicted structures, e.g. the docking poses generated by a molecular docking program. This feature may make KGS a valuable strategy for "real" drug design tasks, such as structure-based virtual screening. It should be mentioned that the standard deviations at different similarity cutoffs in Set II results are somewhat larger than the counterparts in Set I results. It is understandable since for some complexes in this test set, the native binding pose of the ligand molecule is not necessarily the best-scored binding pose selected by scoring function, which is a well-known phenomenon [22, 26, 27]. This defective aspect of scoring function of course introduces extra noises in binding affinity prediction.
Performance on carbonic anhydrase complexes
When only the crystal structures are considered in scoring (Set I), one can see that the standard deviations produced by X-Score+KGS at various similarity cutoffs are consistently lower than or comparable to those produced by X-Score alone (Figure 6A). The same trend retains when the scoring method is switched to PLP. When the docking poses are considered in scoring (Set II), the standard deviations produced by X-Score+KGS and PLP+KGS are basically lower than those produced by X-Score and PLP alone when the similarity cutoff is above 0.40 (Figure 6B). In both cases, however, the improvement after the application of the KGS strategy is not as significant as the one observed on the HIV protease test set. It is understandable since as mentioned above, the performance of X-Score and PLP alone is already good on this test set, leaving not much room for improvement. Considering the intrinsic accuracy of the employed scoring functions as well as the uncertainties in experimental binding data, there is certainly a limit on the average accuracy of binding affinity prediction. Note that after the similarity cutoff in defining the reference complexes is raised above 0.40, the standard deviations produced by X-Score+KGS and PLP+KGS are consistently below 0.60 log units (corresponding to ~0.8 kcal/mol in binding free energy at room temperature). We believe that this level of accuracy, if having not reached the limit of scoring functions, should be rather close to it.
Another trend observed in Figure 6 is that X-Score+KGS and PLP+KGS produced comparable statistical results especially when the similarity cutoffs are relatively high. In contrast, there is noticeable difference in the statistical results produced by X-Score and PLP alone under the same circumstances. In fact, exactly the same trend can be observed on the HIV protease test set as well (Figure 4). Each scoring function has its own strength and weakness, and thus one would expect that different scoring functions produce different results on the given systems. Therefore, the users have to test on their selected targets a number of scoring functions or even combinations of scoring functions before any prediction can be made. Our results indicate that once scoring functions are combined with the KGS strategy, they tend to produce converged results since the difference between their outcomes are largely leveled off by the use of a reference. This feature will bring great convenience to the users in practice.
Performance on trypsin complexes
Information on the six pairs of trypsin complexes with the highest similarities
The query complex
The reference complex
-log K d (exp)
-log K d (exp)
The above analysis suggests that failure of the KGS strategy in this case is caused directly by the intrinsic inaccuracy of the employed scoring function. In theory, our KGS strategy can effectively reduce the systematic errors in the predicted absolute binding constants. It thus becomes critical for the employed scoring function to produce correct rankings for the given complexes, a feature termed as "ranking power" in our recent study . In that study, we demonstrated that today's best scoring functions are able to provide correct rankings for only 50-60% of the protein-ligand complex families under consideration. A much improved "ranking power" perhaps should be the primary aim for future scoring functions. Technically, developers of scoring functions may want to examine the ranking coefficient, such as the Spearman coefficient, more closely than the conventional Pearson coefficient between experimental binding constants and computed binding scores.
Another possible reason for the failure of the KGS strategy in this case lies in our algorithm for pharmacophore elucidation. In our study, a pharmacophore model actually represents a set of key protein-ligand interactions, and our KGS strategy relies on pharmacophore models for defining reference complexes. Such a pharmacophore model is dependent on the compositions as well as the conformations of both the protein and the ligand. Our algorithm only considers three types of features, i.e. hydrogen bond donor, hydrogen bond acceptor, and hydrophobic center. Other basic protein-ligand interactions, such as cation-pi interaction, are ignored for convenience. Besides, a pharmacophore has to contain a limited number of features to be practical. Therefore, some aspects in protein-ligand interactions may be missing. Simplification made in a pharmacophore model may lead to the choice of an inappropriate reference complex. Given these flaws, we have observed that our algorithm in some cases produced results controversial to commonsense. For example, complexes 1O2J and 1O2K are formed by an identical protein (trypsin) and an identical ligand (Figure 8), but the similarity between their pharmacophore models produced by our algorithm is only 0.58. Although these two complexes do exhibit some conformational difference in their binding pockets, one certainly expects a higher similarity score in this case. Our algorithm for pharmacophore elucidation certainly can be improved further.
Comparison with other targeted scoring methods
As mentioned in the Introduction section, current targeted scoring functions [55–58] are developed typically through re-parameterization or modification on existing all-purpose scoring functions, which often relies on some sophisticated statistical procedures. In contrast, our KGS strategy works in a different fashion. There are two basic modules in the framework of KGS: one is a scoring method, and the other is a method for defining the reference complexes. As for the scoring method, we have tested X-Score and PLP in this study. But it can be any well-validated scoring function or other approach. Once the scoring method is chosen, it can be applied as is without re-parameterization or modification. This is important since many programs are available to the end-users as black boxes. As for the method for defining reference complexes, our current algorithm is based on comparison of structure-based pharmacophore models of relevant protein-ligand complexes. Other algorithms of course may be considered as well, such as the protein-ligand interaction fingerprints [71, 72]. Thus, these two modules can be chosen independently, and it is in principle flexible to combine them. This is a notable technical advantage of our KGS strategy, and that is why we refer to it as a strategy rather than a particular method.
Compared to most targeted scoring functions, another technical advantage of our KGS strategy is that its application is in principle not limited to certain classes of targets. In order to apply the KGS strategy, an external database of protein-ligand complexes with known three-dimensional structures and experimental binding data needs to be supplied as a knowledge set. If one needs to consider protein-ligand complexes of various types, this knowledge set may consist of a sufficient number of protein-ligand complexes of various types, such as the PDBbind database. Otherwise, if one's study focuses on a particular class of protein-ligand complexes, one may want to supply a knowledge set only consisting of relevant protein-ligand complexes, such as the test sets used in this study. In practice, most researchers study a certain congeneric class of ligands bound to a common target protein, and thus the latter approach is more suitable for this purpose. Note that in such a case, one can even employ a targeted scoring function, if available, as the internal scoring method for the KGS strategy to obtain more accurate results. In this sense, our KGS strategy is fully compatible with targeted scoring functions.
It should be mentioned that our KGS strategy is similar to the AutoShim method proposed by Martin and Sullivan  in certain aspects. AutoShim does not require re-parameterization of scoring function either, and in principle does not tie to any specific scoring functions. According to AutoShim, point-pharmacophore like "shims" are generated in the binding pocket on the target protein. These "shims" are then weighted by partial least squares (PLS) regression to adjust the outcomes of the Flo+ scoring function in order to better reproduce known binding data. Nevertheless, AutoShim is basically a 3D-QSAR model, which integrates the Flo+ score as well as several hundred of descriptors ("shims"). Such a model normally needs to be carefully validated to avoid over-fitting since it relies on so many parameters. As a matter of fact, Martin et al mentioned that their PLS model produced comparable results with or without the Flo+ score in the entire descriptor set (R2 = 0.60 vs. R2 = 0.56). In contrast, the internal scoring function (e.g. X-Score or PLP) in our KGS strategy plays an indispensable role. Another common aspect between AutoShim and our KGS strategy is that both methods require a training set of protein-ligand complexes with known experimental binding data. As a QSAR-like method, AutoShim relies much on this training set for model calibration. For our KGS strategy, the training set is also used for deriving the scaling parameter k in Equation 4, but it serves primarily as a knowledge set for defining appropriate reference complexes. Using an appropriate reference contributes most to the improved accuracy of our KGS strategy in binding affinity prediction.
Finally, as mentioned repeatedly in this article, our KGS strategy computes the binding constant of a given protein-ligand complex based on the known binding constant of a reference complex. It is essentially an interpolation method. We all know that more accurate results can be obtained through interpolation if more known data exist in the problem space. The knowledge of protein-ligand binding data is certainly in a constant increase. For example, binding data included in the PDBbind database increase by approximately 25% each year. Thus, application of the KGS strategy will hopefully produce more and more accurate results with the increase in such knowledge. In contrast, other targeted scoring functions benefit indirectly from the increasing knowledge of binding data since re-calibration on a larger data set does not necessarily produce a more accurate model. Besides, the majority of experimental binding data as well as crystal structures are actually owned by pharmaceutical companies, which are not available to the public for understandable reasons. Researchers from pharmaceutical companies can supply their in-house data as the external knowledge set required by the KGS strategy so that their own data can be effectively utilized in binding affinity prediction as well.
We have developed a general strategy, i.e. knowledge-guided scoring (KGS), for improving the accuracy of scoring functions in binding affinity prediction. Our KGS strategy computes the binding constant of a given protein-ligand complex based on the known binding constant of an appropriate reference complex. The reference complex is required to share a similar pattern of protein-ligand interactions to that of the complex of interest. Thus, some uncertain factors in protein-ligand binding, which are difficult to be accurately considered by scoring functions, may cancel out in computation, resulting in more accurate prediction of absolute binding affinities.
Our KGS strategy was evaluated in combination with X-Score and PLP on three sets of protein-ligand complexes. As for the HIV protease complexes, X-Score and PLP alone failed to provide acceptable prediction of binding constants; while both X-Score+KGS and PLP+KGS demonstrated notably improved performance especially when the similarity cutoff used in reference searching was relatively high. As for the carbonic anhydrase complexes and trypsin complexes, both X-Score and PLP were able to provide reasonable results by themselves. Application of the KGS strategy in these two cases only produced marginally better or comparable results due to limited remaining room for improvement. An interesting observation is that X-Score+KGS and PLP+KGS produced converged results despite the difference in the outcomes of X-Score and PLP alone. This prompts that application of the KGS strategy may save the end-users, at least to some extents, from the troublesome evaluation and shopping among different scoring functions. Besides, the standard deviations between the experimental and computed binding constants produced by X-Score+KGS and PLP+KGS are below 1.0 logK d units (corresponding to 1.36 kcal/mol in binding free energy at room temperature) on all three test sets when the similarity cutoff is over 0.40. This level of accuracy is certainly acceptable for high-throughput tasks in structure-based drug design. It however remains to be verified if application of KGS is able to achieve this level accuracy consistently on other classes of protein-ligand complexes.
In principle, KGS can be applied in combination with any existing scoring functions. Unlike many other targeted scoring functions, KGS does not require re-parameterization of the given scoring function, and its application is not limited to certain classes of protein-ligand complexes. In addition, KGS is essentially an interpolation method, and thus its effectiveness is in theory proportional to the ever-increasing knowledge of experimental binding data. In-house collections of binding data can be effectively utilized by KGS in computation as well. Compared to other targeted scoring functions, these features make our KGS strategy a more practical remedy for current scoring functions to improve their accuracy in binding affinity prediction. Nevertheless, we did not attempt to test the KGS strategy in this study to see if it also helps with binding mode prediction, or achieves higher success rates in cross-docking or hit rates in virtual screening. For those purposes, we believe that our current algorithms for applying KGS may need certain adjustments.
The authors are grateful to the financial supports from the Chinese National Natural Science Foundation (Grants No.20772149 & No.90813006), the Chinese Ministry of Science and Technology (Grants No.2006AA02Z337 & 2009ZX09501-002), and the Science and Technology Commission of Shanghai Municipality (Grant No.074319113). The authors also thank Dr. Xun Li, Dr. Fu Lin, Dr. Yan Li, Jie Li and Jiajie Li for their contributions to the development and maintenance of the PDBbind database.
- Lybrand TP: Ligand-Protein Docking and Rational Drug Design. Curr Opin Struct Biol 1995, 5(2):224–228. 10.1016/0959-440X(95)80080-8View ArticlePubMedGoogle Scholar
- Kollman PA: Theory of Macromolecule-Ligand Interactions. Curr Opin Struct Biol 1994, 4(2):240–245. 10.1016/S0959-440X(94)90315-8View ArticleGoogle Scholar
- Fischer E: Einfluss der Configuration auf die Wirkung der Enzyme. Ber Dtsch Chem Ges 1894, 27: 2985–2993. 10.1002/cber.18940270364View ArticleGoogle Scholar
- McInnes C: Virtual Screening Strategies in Drug Discovery. Curr Opin Chem Biol 2007, 11(5):494–502. 10.1016/j.cbpa.2007.08.033View ArticlePubMedGoogle Scholar
- Shoichet BK: Virtual screening of chemical libraries. Nature 2004, 432(7019):862–865. 10.1038/nature03197View ArticlePubMedPubMed CentralGoogle Scholar
- Lyne PD: Structure-Based Virtual Screening: an Overview. Drug Discov Today 2002, 7(20):1047–1055. 10.1016/S1359-6446(02)02483-2View ArticlePubMedGoogle Scholar
- Walters WP, Stahl MT, Murcko MA: Virtual Screening - an Overview. Drug Discov Today 1998, 3: 160–178. 10.1016/S1359-6446(97)01163-XView ArticleGoogle Scholar
- Warren GL, Andrews CW, Capelli AM, Clarke B, LaLonde J, Lambert MH, Lindvall M, Nevins N, Semus SF, Senger S, et al.: A Critical Assessment of Docking Programs and Scoring Functions. J Med Chem 2006, 49(20):5912–5931. 10.1021/jm050362nView ArticlePubMedGoogle Scholar
- Kellenberger E, Rodrigo J, Muller P, Rognan D: Comparative evaluation of eight docking tools for docking and virtual screening accuracy. Proteins: Struct, Funct, Bioinf 2004, 57(2):225–242. 10.1002/prot.20149View ArticleGoogle Scholar
- Ewing TJA, Makino S, Skillman AG, Kuntz ID: DOCK 4.0: Search Strategies for Automated Molecular Docking of Flexible Molecule Databases. J Comput-Aided Mol Des 2001, 15(5):411–428. 10.1023/A:1011115820450View ArticlePubMedGoogle Scholar
- Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ: Automated Docking Using a Lamarckian Genetic Algorithm and an Empirical Binding Free Energy Function. J Comput Chem 1998, 19(14):1639–1662. 10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-BView ArticleGoogle Scholar
- Morris GM, Goodsell DS, Huey R, Olson AJ: Distributed Automated Docking of Flexible Ligands to Proteins: Parallel Applications of AutoDock 2.4. J Comput-Aided Mol Des 1996, 10(4):293–304. 10.1007/BF00124499View ArticlePubMedGoogle Scholar
- Goodsell DS, Olson AJ: Automated Docking of Substrates to Proteins by Simulated Annealing. Proteins: Struct, Funct, Genet 1990, 8(3):195–202. 10.1002/prot.340080302View ArticleGoogle Scholar
- Rarey M, Kramer B, Lengauer T, Klebe G: A Fast Flexible Docking Method Using an Incremental Construction Algorithm. J Mol Biol 1996, 261(3):470–489. 10.1006/jmbi.1996.0477View ArticlePubMedGoogle Scholar
- Jain A: Surflex-Dock 2.1: Robust Performance from Ligand Energetic Modeling, Ring Flexibility, and Knowledge-Based Search. J Comput-Aided Mol Des 2007, 21(5):281–306. 10.1007/s10822-007-9114-2View ArticlePubMedGoogle Scholar
- Jain AN: Surflex: Fully Automatic Flexible Molecular Docking Using a Molecular Similarity-Based Search Engine. J Med Chem 2003, 46(4):499–511. 10.1021/jm020406hView ArticlePubMedGoogle Scholar
- Venkatachalam CM, Jiang X, Oldfield T, Waldman M: LigandFit: A Novel Method for the Shape-Directed Rapid Docking of Ligands to Protein Active Sites. J Mol Graphics Modell 2003, 21(4):289–307. 10.1016/S1093-3263(02)00164-XView ArticleGoogle Scholar
- Jones G, Willett P, Glen RC, Leach AR, Taylor R: Development and Validation of a Genetic Algorithm for Flexible Docking. J Mol Biol 1997, 267(3):727–748. 10.1006/jmbi.1996.0897View ArticlePubMedGoogle Scholar
- Jones G, Willett P, Glen RC: Molecular Recognition of Receptor Sites Using a Genetic Algorithm with a Description of Desolvation. J Mol Biol 1995, 245(1):43–53. 10.1016/S0022-2836(95)80037-9View ArticlePubMedGoogle Scholar
- Friesner RA, Murphy RB, Repasky MP, Frye LL, Greenwood JR, Halgren TA, Sanschagrin PC, Mainz DT: Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes. J Med Chem 2006, 49(21):6177–6196. 10.1021/jm051256oView ArticlePubMedGoogle Scholar
- Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, Repasky MP, Knoll EH, Shelley M, Perry JK, et al.: Glide: A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy. J Med Chem 2004, 47(7):1739–1749. 10.1021/jm0306430View ArticlePubMedGoogle 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(4):1079–1093. 10.1021/ci9000053View ArticlePubMedGoogle Scholar
- Brewerton SC: The use of protein-ligand interaction fingerprints in docking. Curr Opin Drug Discovery Dev 2008, 11(3):356–364.Google Scholar
- Ferrara P, Gohlke H, Price DJ, Klebe G, Brooks CL: Assessing Scoring Functions for Protein-Ligand Interactions. J Med Chem 2004, 47(12):3032–3047. 10.1021/jm030489hView ArticlePubMedGoogle Scholar
- Marsden PM, Puvanendrampillai D, Mitchell JBO, Glen RC: Predicting protein-ligand binding affinities: a low scoring game? Org Biomol Chem 2004, 2: 3267–3273. 10.1039/b409570gView ArticlePubMedGoogle Scholar
- Wang R, Lu Y, Fang X, Wang S: An Extensive Test of 14 Scoring Functions Using the PDBbind Refined Set of 800 Protein-Ligand Complexes. J Chem Inf Comput Sci 2004, 44(6):2114–2125.View ArticlePubMedGoogle Scholar
- Wang R, Lu Y, Wang S: Comparative Evaluation of 11 Scoring Functions for Molecular Docking. J Med Chem 2003, 46(12):2287–2303. 10.1021/jm0203783View ArticlePubMedGoogle Scholar
- Kollman P: Free Energy Calculations: Applications to Chemical and Biochemical Phenomena. Chem Rev 1993, 93(7):2395–2417. 10.1021/cr00023a004View ArticleGoogle Scholar
- Jorgensen WL: Free Energy Calculations: A Breakthrough for Modeling Organic Chemistry in Solution. Acc Chem Res 1989, 22(5):184–189. 10.1021/ar00161a004View ArticleGoogle Scholar
- Carlson HA, Jorgensen WL: An Extended Linear Response Method for Determining Free Energies of Hydration. J Phys Chem 1995, 99(26):10667–10673. 10.1021/j100026a034View ArticleGoogle Scholar
- Aqvist J, Medina C, Samuelsson J-E: A New Method for Predicting Binding Affinity in Computer-Aided Drug Design. Protein Eng 1994, 7(3):385–391. 10.1093/protein/7.3.385View ArticlePubMedGoogle Scholar
- Kollman PA, Massova I, Reyes C, Kuhn B, Huo S, Chong L, Lee M, Lee T, Duan Y, Wang W, et al.: Calculating Structures and Free Energies of Complex Molecules: Combining Molecular Mechanics and Continuum Models. Acc Chem Res 2000, 33(12):889–897. 10.1021/ar000033jView ArticlePubMedGoogle Scholar
- Zhang X, Li X, Wang R: Interpretation of the Binding Affinities of PTP1B Inhibitors with the MM-GB/SA Method and the X-Score Scoring Function. J Chem Inf Model 2009, 49(4):1033–1048. 10.1021/ci8004429View ArticlePubMedGoogle Scholar
- Krammer A, Kirchhoff PD, Jiang X, Venkatachalam CM, Waldman M: LigScore: A Novel Scoring Function for Predicting Binding Affinities. J Mol Graphics Modell 2005, 23(5):395–407. 10.1016/j.jmgm.2004.11.007View ArticleGoogle Scholar
- Wang R, Lai L, Wang S: Further Development and Validation of Empirical Scoring Functions for Structure-Based Binding Affinity Prediction. J Comput-Aided Mol Des 2002, 16(1):11–26. 10.1023/A:1016357811882View ArticlePubMedGoogle Scholar
- Wang R, Liu L, Lai L, Tang Y: SCORE: A New Empirical Method for Estimating the Binding Affinity of a Protein-Ligand Complex. J Mol Model 1998, 4(12):379–394. 10.1007/s008940050096View ArticleGoogle Scholar
- Gehlhaar DK, Bouzida D, Rejto PA, eds: Rational Drug Design: Novel Methodology and Practical Applications. American Chemical Society: Washington DC; 1999.
- Gehlhaar DK, Verkhivker GM, Rejto PA, Sherman CJ, Fogel DR, 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(5):317–324. 10.1016/1074-5521(95)90050-0View ArticlePubMedGoogle 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, Genet 1998, 33(3):367–382. 10.1002/(SICI)1097-0134(19981115)33:3<367::AID-PROT6>3.0.CO;2-WView ArticleGoogle Scholar
- 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(5):425–445. 10.1023/A:1007996124545View ArticlePubMedGoogle Scholar
- Jain AN: Scoring Noncovalent Protein-Ligand Interactions: A Continuous Differentiable Function Tuned to Compute Binding Affinities. J Comput-Aided Mol Des 1996, 10(5):427–440. 10.1007/BF00124474View ArticlePubMedGoogle Scholar
- Böhm H-J: Prediction of Binding Constants of Protein Ligands: A Fast Method for the Prioritization of Hits Obtained from de novo Design or 3D Database Search Programs. J Comput-Aided Mol Des 1998, 12(4):309–323. 10.1023/A:1007999920146View ArticlePubMedGoogle Scholar
- Böhm H-J: The Development of a Simple Empirical Scoring Function to Estimate the Binding Constant for a Protein-Ligand Complex of Known Three-Dimensional Structure. J Comput-Aided Mol Des 1994, 8(3):243–256. 10.1007/BF00126743View ArticlePubMedGoogle Scholar
- Mooij WTM, Verdonk ML: General and Targeted Statistical Potentials for Protein-Ligand Interactions. Proteins: Struct, Funct, Bioinf 2005, 61(2):272–287. 10.1002/prot.20588View ArticleGoogle Scholar
- Velec HFG, Gohlke H, Klebe G: DrugScoreCSD - Knowledge-Based Scoring Function Derived from Small Molecule Crystal Data with Superior Recognition Rate of Near-Native Ligand Poses and Better Affinity Prediction. J Med Chem 2005, 48(20):6296–6303. 10.1021/jm050436vView ArticlePubMedGoogle Scholar
- Gohlke H, Hendlich M, Klebe G: Knowledge-Based Scoring Function to Predict Protein-Ligand Interactions. J Mol Biol 2000, 295(2):337–356. 10.1006/jmbi.1999.3371View ArticlePubMedGoogle Scholar
- Muegge I: PMF Scoring Revisited. J Med Chem 2006, 49(20):5895–5902. 10.1021/jm050038sView ArticlePubMedGoogle Scholar
- Muegge I: Effect of Ligand Volume Correction on PMF Scoring. J Comput Chem 2001, 22(4):418–425. 10.1002/1096-987X(200103)22:4<418::AID-JCC1012>3.0.CO;2-3View ArticleGoogle Scholar
- Muegge I: A Knowledge-Based Scoring Function for Protein-Ligand Interactions: Probing the Reference State. Perspect Drug Discov Des 2000, 20: 99–114. 10.1023/A:1008729005958View ArticleGoogle Scholar
- Muegge I, Martin YC: A General and Fast Scoring Function for Protein-Ligand Interactions: A Simplified Potential Approach. J Med Chem 1999, 42(5):791–804. 10.1021/jm980536jView ArticlePubMedGoogle Scholar
- Huang S-Y, Zou X: An iterative knowledge-based scoring function to predict protein-ligand interactions: II. Validation of the scoring function. J Comput Chem 2006, 27(15):1876–1882. 10.1002/jcc.20505View ArticlePubMedGoogle Scholar
- Huang S-Y, Zou X: An iterative knowledge-based scoring function to predict protein-ligand interactions: I. Derivation of interaction potentials. J Comput Chem 2006, 27(15):1866–1875. 10.1002/jcc.20504View ArticlePubMedGoogle Scholar
- Teramoto R, Fukunishi H: Consensus Scoring with Feature Selection for Structure-Based Virtual Screening. J Chem Inf Model 2008, 48(2):288–295. 10.1021/ci700239tView ArticlePubMedGoogle Scholar
- Teramoto R, Fukunishi H: Supervised Scoring Models with Docked Ligand Conformations for Structure-Based Virtual Screening. J Chem Inf Model 2007, 47(5):1858–1867. 10.1021/ci700116zView ArticlePubMedGoogle Scholar
- Seifert MHJ: Targeted scoring functions for virtual screening. Drug Discov Today 2009, 14(11–12):562–569. 10.1016/j.drudis.2009.03.013View ArticlePubMedGoogle Scholar
- Pfeffer P, Gohlke H: DrugScoreRNA-Knowledge-Based Scoring Function To Predict RNA-Ligand Interactions. J Chem Inf Model 2007, 47(5):1868–1876. 10.1021/ci700134pView ArticlePubMedGoogle Scholar
- Antes I, Merkwirth C, Lengauer T: POEM: Parameter Optimization Using Ensemble Methods: Application to Target Specific Scoring Functions. J Chem Inf Model 2005, 45(5):1291–1302. 10.1021/ci050036gView ArticlePubMedGoogle Scholar
- Seifert MHJ: Optimizing the Signal-to-Noise Ratio of Scoring Functions for Protein-Ligand Docking. J Chem Inf Model 2008, 48(3):602–612. 10.1021/ci700345nView ArticlePubMedGoogle Scholar
- Cheng T, Zhao Y, Li X, Lin F, Xu Y, Zhang X, Li Y, Wang R, Lai L: Computation of Octanol-Water Partition Coefficients by Guiding an Additive Model with Knowledge. J Chem Inf Model 2007, 47(6):2140–2148. 10.1021/ci700257yView ArticlePubMedGoogle Scholar
- Kuhne R, Ebert R-U, Schuurmann G: Model Selection Based on Structural Similarity-Method Description and Application to Water Solubility Prediction. J Chem Inf Model 2006, 46(2):636–641. 10.1021/ci0503762View ArticlePubMedGoogle Scholar
- Chen J, Lai L: Pocket v.2: Further Developments on Receptor-Based Pharmacophore Modeling. J Chem Inf Model 2006, 46(6):2684–2691. 10.1021/ci600246sView ArticlePubMedGoogle Scholar
- Wang R, Gao Y, Lai L: LigBuilder: A Multi-Purpose Program for Structure-Based Drug Design. J Mol Model 2000, 6(7):498–516. 10.1007/s0089400060498View ArticleGoogle Scholar
- Bron C, Kerbosch J: Algorithm 457: Finding All Cliques of an Undirected Graph. Commun ACM 1973, 16(9):575–577. 10.1145/362342.362367View ArticleGoogle Scholar
- Willett P, Barnard JM, Downs GM: Chemical Similarity Searching. J Chem Inf Comput Sci 1998, 38(6):983–996.View ArticleGoogle Scholar
- Wang R, Fang X, Lu Y, Yang CY, Wang S: The PDBbind Database: Methodologies and Updates. J Med Chem 2005, 48(12):4111–4119. 10.1021/jm048957qView ArticlePubMedGoogle 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(12):2977–2980. 10.1021/jm030580lView ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucl Acids Res 2000, 28(1):235–242. 10.1093/nar/28.1.235View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao Y, Cheng T, Wang R: Automatic Perception of Organic Molecules Based on Essential Structural Information. J Chem Inf Model 2007, 47(4):1379–1385. 10.1021/ci700028wView ArticlePubMedGoogle Scholar
- The Sybyl Software (version 7.2) Tripos Inc.: St. Louis, MO; 2006.
- The Discovery Studio software (version 2.0) Accelrys Software Inc.: San Diego, CA 92121 USA; 2001.
- Baroni M, Cruciani G, Sciabola S, Perruccio F, Mason JS: A Common Reference Framework for Analyzing/Comparing Proteins and Ligands. Fingerprints for Ligands And Proteins (FLAP): Theory and Application. J Chem Inf Model 2007, 47(2):279–294. 10.1021/ci600253eView ArticlePubMedGoogle Scholar
- Marcou G, Rognan D: Optimizing Fragment and Scaffold Docking by Use of Molecular Interaction Fingerprints. J Chem Inf Model 2007, 47(1):195–207. 10.1021/ci600342eView ArticlePubMedGoogle Scholar
- Martin EJ, Sullivan DC: AutoShim: Empirically Corrected Scoring Functions for Quantitative Docking with a Crystal Structure and IC50 Training Data. J Chem Inf Model 2008, 48(4):861–872. 10.1021/ci7004548View ArticlePubMedGoogle 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 cited.