PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of springeropenLink to Publisher's site
Journal of Computer-Aided Molecular Design
 
J Comput Aided Mol Des. 2010 March; 24(3): 173–182.
Published online 2010 March 13. doi:  10.1007/s10822-009-9316-x
PMCID: PMC2858276

Spatial chemical distance based on atomic property fields

Abstract

Similarity of compound chemical structures often leads to close pharmacological profiles, including binding to the same protein targets. The opposite, however, is not always true, as distinct chemical scaffolds can exhibit similar pharmacology as well. Therefore, relying on chemical similarity to known binders in search for novel chemicals targeting the same protein artificially narrows down the results and makes lead hopping impossible. In this study we attempt to design a compound similarity/distance measure that better captures structural aspects of their pharmacology and molecular interactions. The measure is based on our recently published method for compound spatial alignment with atomic property fields as a generalized 3D pharmacophoric potential. We optimized contributions of different atomic properties for better discrimination of compound pairs with the same pharmacology from those with different pharmacology using Partial Least Squares regression. Our proposed similarity measure was then tested for its ability to discriminate pharmacologically similar pairs from decoys on a large diverse dataset of 115 protein–ligand complexes. Compared to 2D Tanimoto and Shape Tanimoto approaches, our new approach led to improvement in the area under the receiver operating characteristic curve values in 66 and 58% of domains respectively. The improvement was particularly high for the previously problematic cases (weak performance of the 2D Tanimoto and Shape Tanimoto measures) with original AUC values below 0.8. In fact for these cases we obtained improvement in 86% of domains compare to 2D Tanimoto measure and 85% compare to Shape Tanimoto measure. The proposed spatial chemical distance measure can be used in virtual ligand screening.

Electronic supplementary material

The online version of this article (doi:10.1007/s10822-009-9316-x) contains supplementary material, which is available to authorized users.

Keywords: Chemical distance, Chemical similarity, Quantitative structure-property relationship, Internal coordinates mechanics (ICM), Atomic property fields (APF), Spatial alignment

Introduction

Ligand-based approaches to protein family profiling has been widely studied and used for in silico pharmacology [1]. Similarity of compound chemical structures often leads to close pharmacological profiles, including binding to the same protein targets. By this reason, chemical similarity criterion is widely used for identification of novel lead molecules in the development of pharmaceuticals. A variety of chemical similar measures has been proposed. However, in many cases compounds with similar pharmacology escape correct recognition as they appear to be dissimilar by any existing measure.

In order to navigate in ligand space, one need to represent the compound using appropriate properties (descriptors) and then use a master equation to measure a distance between two compounds.

Descriptors are usually classified according to their dimensionality ranging from one-dimensional (1-D) to three-dimensional (3D) properties [2, 3, 10]. Easy and fast to compute 1-D descriptors describe global properties which can be derived from chemical formula and classify compounds or ligands from various target families [35, 10]. To perform fast comparison 1-D linear representations of compounds are often used. The most popular of this kind of simplified string is the ‘Simplified Molecular Input Line Entry System’ or SMILES [3, 6, 10].

To improve discrimination, 2D topological descriptors are used. Graph-based methods, such as maximum common subgraph (MCS) [3, 7, 10] and fingerprint-based methods [3, 8, 10] are popular for substructure clustering chemical compounds into subfamilies. Subgraph isomorphism in large molecular databases is quite often time consuming to perform on large numbers of structures and it was for this reason that substructure screening was developed as a rapid method of filtering out those molecules that definitely do not contain the substructure of interest [10, 46]. The similarity between two molecules represented by 2D binary fingerprints is most frequently quantified using the Tanimoto coefficient, which gives a measure of the number of fragments in common between the two molecules [3, 9, 10].

It is well known that molecular recognition depends on the 3D structure and properties of molecule rather than the underlying substructure(s) [10]. 3D methods are computationally more expensive than 2D descriptor based methods, because they require consideration of conformational space of the molecule. These methods can be divided into methods that are alignment-independent and methods that require the molecules to be aligned in 3D space before similarity function is used [10].

Some computationally expensive alignment-independent methods use 3D geometrical descriptors represent them in a binary fingerprint and then use with the Tanimoto coefficient exactly as for 2D fingerprints [10, 11]. Other methods are 3D equivalent of the MCS [10, 12, 13]. Many 3D approaches are based on the use of distances matrices where the value of each element (i, j) equals the interatomic distance between atoms i and j [10, 14]. Also there are approaches where the pharmacophore points are used for similarity comparisons [10, 1517].

Consideration of conformational flexibility of the molecules as well as their relative orientation is required for alignment dependent methods [10]. These methods devised to align the compared structures via maximization of the similarity function that is used [10, 45]. Many different ways have been developed to represent molecules and calculate similarity based on molecular shape and/or field [1831, 45]. For reviews of molecular similarity methods, see refs [2, 10, 3236].

The aim of this study is to design a spatial distance measure between two chemicals that optimizes recognition of their pharmacological similarity by using their 3D conformational ensembles and properties pertaining to molecular interactions. We recently introduced a novel spatial alignment method based on atomic property fields (APF) as a generalized 3D pharmacophoric potential [37]. APF is the representation of the ligand by a multi-component (vector) 3D potential, with the components corresponding to various physico-chemical atomic properties. In the present study, the APF alignment is used to measure spatial chemical similarity/distance between ligands.

A diverse benchmark of 99 proteins (see Supplementary Table 1 for details) and ligands co-crystallized with these proteins (with 6 ligands per protein on average) was used to train APF parameters for better discrimination of pharmacologically similar pairs from dissimilar ones. All possible combinations of pairs of ligands from the same receptors as well as for ligands co-crystallized with certain protein all possible combinations of pairs with ligands co-crystallized with 20 different randomly chosen from benchmark other proteins, were taken and APF representation of larger ligand was used as a reference to superimpose them. Distances between all superimposed pairs of ligands have been evaluated. Performance and results of proposed approach are reported and systematically compared to those obtained with standard 2D fingerprints (using 2D Tanimoto equation [38]) and 3D shape (using shape Tanimoto similarity measure [38]) based approaches (see Supplementary Table 2 for details).

Materials and methods

Atomic property field

Atomic property field (APF) [37] is the representation of the ligand by multi-component (vector) 3D grid potential Pi (equation M1), with the components i corresponding to various physico-chemical properties. Each property component of the APF determines whether the presence at any specific point equation M2 in space of an atom with that particular property is favorable or unfavorable. Pseudo-energy (score) of an atom j in this field is a dot product of its property vector equation M3 and the APF potential at its position equation M4

equation M5

The minima for this pseudo-energy for an atom with any specific property vector equation M6 will be in the areas of space with similar APF potential components.

Assignment of property vectors equation M7 for various atom types is carried out as described in [37], according to the general knowledge of their empiric physico-chemical behavior. Seven property field components were introduced: hydrogen bond donor, hydrogen bond acceptor, sp2 hybridized, lipophilic, size, charged, and electronegative/electropositive. Five of them (hydrogen bond donor, hydrogen bond acceptor, sp2 hybridized, lipophilic and charged) are classic pharmacophoric types. Other two (size and electronegative/electropositive) are extension to make APF vectors differentiate certain atom types that are indistinguishable by the first five components: for example, aliphatic carbon (as in a methyl group) and large halogens are differentiated by electronegativity. For more details, see [37].

APF-based local ligand superimposition

Suboptimal ligand superimpositions were obtained for all pairs of compounds binding to the same protein by superimposing the protein binding pockets. These suboptimal superimpositions were improved by locally minimizing APF-correspondence between the two ligands. APF representation was generated from the larger ligand. The second ligand was flexibly minimized in the obtained field in a search for the local minimum of its APF-score calculated as described above. The local minimization procedure used the Newton method [39] if the number free variables of a given molecule are less than 100 and switched to the conjugate gradient method [39] if number of free variables exceeded 100. The maximal number of function evaluations was set to N2, where N is number of free variables of a molecule. This local minimization procedure is implemented in the ICM software [40].

APF-based global ligand superimposition

For all ligand pairs that do not bind to the same protein, global ligand superimposition was performed by Monte Carlo search of the minimum of the APF-score of the smaller ligand in the APF of the larger. Similarly, the maximal number of function evaluations was set to N2, where N is number of free variables of a given molecule. The global minimization procedure is implemented in the ICM software [40].

Dataset of protein–ligand complexes

The proposed approach to similarity measure was validated on a comprehensive benchmark of protein–ligand complexes from RCSB Protein Data Bank with high resolution and drug-like ligands (Kufareva et al. manuscript in preparation). For this validation, all ligands were paired up with one another to form the so-called correct and decoy pairs. A pair of ligands was considered correct if these ligands exhibited similar pharmacology exemplified by binding to the same protein target, as found in PDB. If two ligands were not found in co-crystal complexes with any single protein, they formed an incorrect, or decoy pair. As the number of decoy pairs was several order of magnitude larger than the number of correct pairs, we compressed the decoy set by pairing each ligand up with all ligands from 20 randomly chosen proteins (rather than all other proteins), excluding the cognate receptor.

The dataset was randomly split into a training set and a test set. The training set consisted of 2,538 correct and 9,819 decoy pairs, so that their numbers are in the same order of magnitude. The test set consisted of 1898 correct and 317,807 decoy pairs. For the training data set, correct ligand pairs were superimposed locally starting from their crystallographic poses as described above, while the decoy pairs were superimposed globally. For the test set all pairs were superimposed using global APF superimposition described above.

Partial least square regression

Partial least square (PLS) regression [41, 42] was used to design the proposed similarity measure. PLS is a recent technique that generalizes and combines features from principal component analysis and multiple regression. The purpose of PLS is to approximate a dependent (target) variable Y with a linear combination of independent variables Xi (descriptors). The output of the learning algorithm then consists of a linear coefficient ci for each of the descriptors Xi and a free term c0.

For each ligand pair in the training set, we assigned the value of Y equal to 1 for correct pairs and 0 for decoy pairs. PLS regression was trained to approximate these values as a linear combination of the seven descriptors equation M8 The optimal vector of coefficients, ci, was used to access pairwise similarity if the compounds in the test set. The relative contributions wi of seven descriptors, further referred as weights, were evaluated as equation M9. A recent implementation of PLS in ICM software [40] was used.

Design of novel spatial chemical distance measure

To design a novel and optimal spatial chemical distance measure we undertook the following approach (see Fig. 1). For each ligand pair, we performed the following steps. First, the larger ligand (greater number of heavy atoms) was locally relaxed from its bioactive conformation in its receptor and its APF representation was built. Second, the smaller ligand was flexibly superimposed onto the larger ligand using local and global APF-based superimposition for correct and decoy pairs, respectively. Third, the pseudo-energies, equation M10 and equation M11 of the larger and the smaller ligands in the APF potential grid of the larger ligand were calculated.

Fig. 1
Schematic outline of approach to design a novel spatial chemical distance measure

Further on, the seven ratios equation M12 were recorded as descriptors of the ligand pair, and PLS regression was trained to approximate the correct/decoy (1/0) value as a weighted sum

equation M13

The protocol is built on top of a well established and implemented in ICM APF concept and ICM gradient minimization.

The weights wi (i = 1, …,7) obtained by the PLS regression were used in the definition of the compound chemical distance as:

equation M14

2D tanimoto distance calculation

Compounds chemical similarity was evaluated as Tanimoto distance between their molecular fingerprints as implemented in ICM [40]. Briefly, given a molecule, all linear and non-linear fragments of different size were enumerated and hashed into a bit string called a fingerprint. The Tanimoto coefficient, T, for two fingerprints was calculated as the number of bits in which they differ divided by the number of non-zero bits they have in common. The Tanimoto distance was defined as 1 − T.

3D shape tanimoto distances

For each compound, the volume of its molecular envelope was calculated using ICM. For every pair of compounds superimposed by the APF-based algorithm, the volume of the smallest envelope enclosing both compounds, V, was calculated. The overlapping volume, V12, was found by V12 = V1 + V2 − V. The 3D shape Tanimoto distance [38] between the two compounds was calculated as

equation M15

where V12 is the volume overlap between two compounds, V1 and V2 are volumes of compounds.

Distance distribution analysis

Distance probability distribution curves were built by normalization (division by number of observations) of histogram plots of 2D Tanimoto, 3D Shape Tanimoto, and our newly designed distances for all correct and decoy pairs from the test set. Global discrimination abilities of 2D Tanimoto, 3D Shape Tanimoto, and our newly designed distance measure are evaluated by calculating an overlapping area between distance probability distribution curves for all correct and decoy pairs from the test set.

ROC AUC values

The new measure was evaluated for its ability to distinguish correct compound pairs from decoy pairs on the large test set. The performance of the measure was evaluated as the area under the ROC (receiver operating characteristic) curve, or ROC AUC. ROC curve analysis [43] describes the ability of a screening method to avoid false positives and false negatives. The ideal screening device demonstrates the ROC AUC value of 1, while a random selection performance corresponds to the ROC AUC value of 0.5.

ROC AUC analysis was performed for 2D Tanimoto, 3D Shape Tanimoto, and our newly designed chemical distance measure for the entire test set. Also ROC AUC analysis was performed for 2D Tanimoto, 3D Shape Tanimoto, and our newly designed chemical distance measure for each of the 115 protein ensembles in the testing set.

Statistical analysis

To evaluate the statistical significance of the obtained ROC AUC values, we assumed that they are distributed according to Gaussian distribution with mean 0.5. The uncertainty σ is calculated by performing 20 random statistical experiments in which the ranks are reshuffled, calculating individual ROC AUCs and calculating the root mean square deviation of those ROC AUCs. For each ROC AUC values obtained on a given subset of correct and decoy pairs, its P-value was calculated as probability to obtain the same AUC by random coincidence:

equation M16

where erf is the error function [44] and σ is the standard deviation. Small P-values below 0.05 mean high statistical significance, while cases with large P-values are statistically less significant.

Software and hardware

Ligand preparation, superimpositions, APF generations, evaluation of pseudo-energies and PLS regression analysis were carried out with ICM 3.6 (Molsoft LLC, La Jolla, CA).

The hardware facility employed in the present study was Intel Core 2 Duo workstation (2.4 GHz with 2 GB of RAM memory).

Results and discussion

APF-based compounds chemical distance measure

In this study, we propose an algorithm for evaluation of compound spatial chemical distance that consists of five consecutive steps (Fig. 2): (1) the ligand with greater number of heavy atoms is locally relaxed from its bioactive coordinates (found in the co-crystal structure with its target protein); (2) APF potential grid maps are built from this ligand; (3) the smaller ligand is flexibly superimposed onto the larger ligand by minimizing its APF-score in the field of the larger ligand; (4) ratios equation M17 are calculated, where equation M18 are the pseudo-energies of the larger ligand in its own APF maps and equation M19 are the pseudo-energies of the smaller ligand in APF representation of the larger one; (5) the similarity score is calculated by (2) and the chemical distance between the compounds is found by (3).

Fig. 2
Schematic flowchart of calculation of spatial chemical distance between two compounds

The weights wi in Eq. (2) were derived by Partial Least Squares regression optimizing the pairwise correct/decoy (1/0) discrimination value as a linear combination of the equation M20 ratios as described in “Materials and Methods”. The results of regression are shown in the Table 1. The APF potential component corresponding to the atomic size was found to contribute the most to the compound similarity. This result indicates the fact that shape complementarity between compounds has a fundamental role, so that compounds with similar shape often have similar pharmacological profiles.

Table 1
Relative and absolute contribution of components of APF potential to spatial chemical distance measure

According to these results, the spatial chemical distance between two compounds can be measured by

equation M21

Screening performance of novel spatial chemical distance measure

The screening performance of the proposed spatial chemical distance measure was tested on the diverse set of 1898 correct (from 115 proteins) and 317,807 decoy pairs (see Supplementary Table 2). Global discrimination abilities of different methods used in this study are presented in the Fig. 3. In this plot, the relative frequencies of three types of compound distances are shown for correct and decoy compound pairs. As expected, the distribution peaks are shifted with respect to each other in all three cases, reflecting the fact that on average, compounds in correct pairs tend to be closer to each other than compounds in decoy pairs. However, it is clear that the newly proposed chemical distance measure provides better separation of the peaks than either 2D fingerprint Tanimoto distance or 3D shape Tanimoto distance. Indeed, by our new measure, only 52% of the correct pairs appear in the distance region of decoy pairs (Fig. 3c). The corresponding number for the other two distance measures appear to be 64 and 56% (Fig. 3a, b). Discrimination abilities of different methods used in this study also can be seen in the Fig. 4, where ROC curves and overall correct/decoy discrimination ROC AUC values for the test set are presented. Our new proposed chemical distance measure gives ROC AUC value of 0.81, in contrast to 0.79 and 0.74 obtained by using 3D shape Tanimoto distance and 2D fingerprint Tanimoto distance. These results clearly shows an advantage of our new approach in lead/decoy discrimination on the benchmark compared to 2D Tanimoto and 3D Shape Tanimoto approaches.

Fig. 3
Probability distribution for correct and decoy pairs: a 2D Tanimoto approach. Overlapped area is 64%; b Shape Tanimoto approach. Overlapped area is 56%; c New approach based on APF. Overlapped area is 52%
Fig. 4
ROC curves and overall correct/decoy discrimination ROC AUC values evaluated by using our novel chemical distance, 3D Shape Tanimoto distance and 2D Tanimoto fingerprint

It also can be seen in the Fig. 5 where the performance of the proposed measure was individually evaluated for each of the 115 protein targets from the test set. By using the novel spatial compounds chemical distance measure, the correct/decoy discrimination ROC AUC value was improved for 66% of the test set when compared to 2D Tanimoto, and 58% when compared to 3D Shape Tanimoto. The improvement was particularly high for the previously problematic cases (weak performance of the 2D Tanimoto and Shape Tanimoto measures) with original AUC values below 0.8. For these cases we obtained improvement in 86% of domains compared to 2D Tanimoto measure and 85% compared to 3D Shape Tanimoto measure. The average ROC AUC value improvement for these cases equals 0.12 and 0.24 with standard deviation 0.15 and 0.24, respectively (Fig. 5a, b).

Fig. 5
Color-coded dependence of achieved improvement in ROC AUC values on original ROC AUC values obtained by 2D Tanimoto (a) and Shape Tanimoto (b) measures. Each dot represents one domain. The dots are colored in different colors spread evenly among P-values ...

Compared to 2D Tanimoto approach, the best improvement (up to 0.92 in AUC units) was achieved for compounds with perfectly matching pharmacophores but no common substructure. An example is given in the Fig. 6, where the two ligands were taken from Geranyltranstransferase and globally superimposed as described above. The newly proposed distance measure returns the value of 0.28 for these two compounds, in contrast to 0.88 obtained by 2D Tanimoto measure. Consequently, using our measure for correct/decoy discrimination for this protein resulted in the ROC AUC value of 0.99, in contrast to 0.07 obtained by 2D Tanimoto approach. As another extreme, there are cases where 2D Tanimoto approach provides better recognition (Fig. 5a, AUCs based on 2D Tanimoto distance greater than 0.8). In general these are the cases where the two compounds have a common substructure but other parts of them are significantly different. However, even in these cases, our new approach performs reasonably accurate and provides systematically high AUCs with 0.85 on average and 0.11 standard deviation.

Fig. 6
Indicative example of weak performance of approach based on 2D Tanimoto similarity measure: a 2D structures of two ligands, which have taken from Geranyltranstransferase; b same ligands in their bioactive conformations; c ligands superimposed by approach ...

There is no doubt that ligand shape, and its complimentary to the shape of the protein interaction site are of fundamental importance for their bioactivity. However, comparison of ligands by shape causes compounds of dissimilar dimensions to be down-ranked during virtual screening. Such ligands, however, may have high affinity to a target as they may form additional interactions with the protein. The results presented in the Fig. 5b), namely up to 0.54 improvement in AUC, show that our new spatial chemical distance measure is able to discover bioactive compounds that are down-ranked by the 3D-shape based measure.

Conclusions

In this study, we developed a spatial chemical distance measure between two chemicals by using their 3D conformational ensembles and properties pertaining to molecular interactions. The parameters for our novel measure were obtained by training a PLS regression to distinguish the correct ligand pairs from decoys on a large and sufficiently diverse set of correct pairs of ligands, which were taken from the same receptor and decoy pairs of ligands, which were taken from different receptors. The screening performance of the proposed spatial chemical distance measure was tested on a diverse and pharmaceutically relevant test set of correct and decoy pairs. Compared to 2D Tanimoto and Shape Tanimoto approaches, our new approach led to improvement in the area under the receiver operating characteristic curve values in 66 and 58% of domains respectively. The improvement was particularly high for the previously problematic cases (weak performance of the 2D Tanimoto and Shape Tanimoto measures) with original AUC values below 0.8. In fact for these cases we obtained improvement in 86% of domains compare to 2D Tanimoto measure and 85% compare to Shape Tanimoto measure.

The presented results suggest that our new spatial chemical distance measure can be successfully used in virtual ligand screening for novel chemical scaffolds targeting the existing proteins of therapeutic relevance.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Acknowledgments

This work was supported by NIH grants 5-R01-GM071872 and 1-R01-GM074832. We are indebted to Manuel Rueda for helpful discussions.

Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

References

1. Mestres J, Martin-Couce L, Gregori-Puigjane E, Cases M, Boyer S. Ligand-based approach to in silico pharmacology: nuclear Receptor Profiling. J Chem Inf Model. 2004;46:2725–2736. doi: 10.1021/ci600300k. [PubMed] [Cross Ref]
2. Bender A, Glen RC. Molecular similarity: a key technique in molecular information. Org Biomol Chem. 2004;2:3204–3218. doi: 10.1039/b409813g. [PubMed] [Cross Ref]
3. Rognan D. Chemogenomic approaches to rational drug design. Br J Pharmacol. 2007;152:38–52. doi: 10.1038/sj.bjp.0707307. [PMC free article] [PubMed] [Cross Ref]
4. Sadowski J, Kubinyi HA. scoring scheme for discriminating between drugs and nondrugs. J Med Chem. 1998;41:3325–3329. doi: 10.1021/jm9706776. [PubMed] [Cross Ref]
5. Morphy R. The influence of target family and functional activity on the physicochemical properties of pre-clinical compounds. J Med Chem. 2006;49:2969–2978. doi: 10.1021/jm0512185. [PubMed] [Cross Ref]
6. Weininger D. SMILES 1. Introduction and encoding rules. J Chem Inf Comput Sci. 1988;28:31–36.
7. Raymond JW, Blankley CJ, Willett P. Comparison of chemical clustering methods using graph- and fingerprint-based similarity measures. J Mol Graph Model. 2003;21:421–433. doi: 10.1016/S1093-3263(02)00188-2. [PubMed] [Cross Ref]
8. Willett P. Similarity-based virtual screening using 2D fingerprints. Drug Discov Today. 2006;11:1046–1053. doi: 10.1016/j.drudis.2006.10.005. [PubMed] [Cross Ref]
9. Willett P, Winterman V, Bawden D. Implementation of nearest neighbor searching in an online chemical structure search system. J Chem Inf Comput Sci. 1986;38:983–996.
10. Leach AR, Gillet VJ. An introduction to cheminformatics. Dordrecht: Springer; 2007.
11. Fisanick W, Lipkus AH, Rusinko A., III Similarity searching on CAS registry substances. 1. Global molecular property and generic atom triangle geometry searching. J Chem Inf Comput Sci. 1992;32:664–674.
12. Moon JB, Howe WJ. 3D database searching and de novo construction methods in molecular design. Tetrahedron Comput Methodol. 1990;3:697–711. doi: 10.1016/0898-5529(90)90168-8. [Cross Ref]
13. Pepperrell CA, Willett P. Techniques for the calculation of three-dimensional structural similarity using inter-atomic distances. J Comput-Aided Mol Des. 1991;5:455–474. doi: 10.1007/BF00125665. [PubMed] [Cross Ref]
14. Pepperrell CA, Willett P, Taylor R. Implementation and use of an atom-mapping procedure for similarity searching in databases of 3D chemical structures. Tetrahedron Comput Methodol. 1990;3:575–593. doi: 10.1016/0898-5529(90)90160-A. [Cross Ref]
15. Good A, Lewis RA. New methodology for profiling combinatorial libraries and screening sets: cleaning up the design with HARPick. J Med Chem. 1997;40:3926–3936. doi: 10.1021/jm970403i. [PubMed] [Cross Ref]
16. Mason JS, Morize I, Menard PR, Cheney DL, Hulme C, Lebaudiniere RF. New 4-point pharmacophore method for molecular similarity and diversity applications: overview of the method and applications, including a novel approach to the design of combinatorial libraries containing privileged substructures. J Med Chem. 1999;42:3251–3264. doi: 10.1021/jm9806998. [PubMed] [Cross Ref]
17. Good AC, Mason JS, Green DVS, Leach AR. Pharmacophore-based approaches to combinatorial library design. In: Ghose AK, Viswanadhan VN, editors. Combinatorial library design and evaluation. Principles, software tools and applications in drug discovery. New York: Marcel Dekker; 2001. pp. 399–428.
18. Kearsley SK, Smith GM. An alternative method for the alignment of molecular structures: maximizing electrostatic and steric overlap. Tetrahedron Comput Methodol. 1990;3:615–633. doi: 10.1016/0898-5529(90)90162-2. [Cross Ref]
19. Bemis GW, Kuntz ID. A fast and efficient method for 2D and 3D molecular shape description. J Comput-Aided Mol Des. 1992;6:607–628. doi: 10.1007/BF00126218. [PubMed] [Cross Ref]
20. Wild D, Willett P. Similarity searching in files of three-dimensional chemical structures. Alignment of molecular electrostatic potential fields with a genetic algorithm. J Chem Inf Comput Sci. 1996;36:159–167.
21. Thorner D, Wild D, Willett P, Wright P. Similarity searching in files of three-dimensional chemical structures: flexible field-based searching of molecular electrostatic potentials. J Chem Inf Comput Sci. 1996;36:900–908.
22. Hahn M. Three-dimensional shape-based searching of conformationally flexible compounds. J Chem Inf Comput Sci. 1997;37:80–86.
23. Rush T, Grant J, Mosyak L, Nicholls A. A shape-based 3-D scaffold hopping method and its application to a bacterial protein-protein interaction. J Med Chem. 2005;48:1489–1495. doi: 10.1021/jm040163o. [PubMed] [Cross Ref]
24. Labute P, Williams C, Feher M, Sourial E, Schmidt J. Flexible alignment of small molecules. J Med Chem. 2001;44:1483–1490. doi: 10.1021/jm0002634. [PubMed] [Cross Ref]
25. Pitman MC, Huber WK, Horn H, Kramer A, Rice JE, Swope WC. FLASHFLOOD: a 3D field-based similarity search and alignment method for flexible molecules. J Comput-Aided Mol Des. 2001;15:587–612. doi: 10.1023/A:1011921423829. [PubMed] [Cross Ref]
26. Putta S, Landrum G, Penzotti J. Conformation mining: an algorithm for finding biologically relevant conformations. J Med Chem. 2005;48:3313–3318. doi: 10.1021/jm049066l. [PubMed] [Cross Ref]
27. Tervo A, Ronkko T, Nyronen T, Poso A. BRUTUS: optimization of a grid-based similarity function for rigid-body molecular superposition. 1. Alignment and virtual screening applications. J Med Chem. 2005;48:4076–4086. doi: 10.1021/jm049123a. [PubMed] [Cross Ref]
28. Rönkkö T, Tervo AJ, Parkkinen J, Poso A. BRUTUS: optimization of a grid-based similarity function for rigid-body molecular superposition. II. Description and characterization. J Comput-Aided Mol Des. 2006;20:227–236. doi: 10.1007/s10822-006-9052-4. [PubMed] [Cross Ref]
29. Cheeseright T, Mackey M, Rose S, Vinter A. Molecular field extrema as descriptors of biological activity: definition and Valida-tion. J Chem Inf Model. 2006;46:665–676. doi: 10.1021/ci050357s. [PubMed] [Cross Ref]
30. Good AC. Novel DOCK clique driven 3D similarity database search tools for molecule shape matching and beyond: adding flexibility to the search for ligand kin. J Mol Graph Model. 2007;26:656–666. doi: 10.1016/j.jmgm.2007.03.016. [PubMed] [Cross Ref]
31. Kirchmair J, Distinto S, Markt P, Schuster D, Spitzer GM, Liedl KR, Wolber G. How to optimize shape-based virtual screening: choosing the right query and including chemical information. J Chem Inf Model. 2009;49:678–692. doi: 10.1021/ci8004226. [PubMed] [Cross Ref]
32. Lemmen C, Lengauer T. Computational methods for the structural alignment of molecules. J Comput-Aided MolDes. 2000;14:215–232. doi: 10.1023/A:1008194019144. [PubMed] [Cross Ref]
33. Eckert H, Bajorath J. Molecular similarity analysis in virtual screening: foundations, limitations and novel approaches. Drug Discov Today. 2007;12:225–233. doi: 10.1016/j.drudis.2007.01.011. [PubMed] [Cross Ref]
34. Willett P, Barnard JM, Downs GM. Chemical similarity searching. J Chem Inf Comput Sci. 1998;38:983–996.
35. Maldonado AG, Doucet JP, Petitjean M, Fan B-T. Molecular similarity and diversity in chemoinformatics: from theory to applications. Mol Divers. 2006;10:39–79. doi: 10.1007/s11030-006-8697-1. [PubMed] [Cross Ref]
36. Good AC, Richards WG (1998) Explicit calculation of 3D molecular similarity. Perspect Drug Discov Design 9/10/11, 321–338
37. Totrov M. Atomic property fields: generalized 3D pharmacophoric potential for automated ligand superposition, pharmacophore elucidation and 3D QSAR. Chem Biol Drug Des. 2008;71:15–27. [PubMed]
38. Fontaine F, Bolton E, Borodina Y, Bryant SH. Fast 3D shape screening of large chemical databases through alignment-recycling. Chem Cent J. 2007;1:12. doi: 10.1186/1752-153X-1-12. [PMC free article] [PubMed] [Cross Ref]
39. Leach AR. Molecular modeling. Principles and applications. Edinburgh: Principles and Applications; 2001.
40. Abagyan R (2008) ICM Manual v. 3.6
41. Geladi P, Kowalski B. Partial least squares regression: a tutorial. Anal Chim Acta. 1986;185:1–17. doi: 10.1016/0003-2670(86)80028-9. [Cross Ref]
42. Hand D, Mannila H, Smyth P. Principles of data mining. Cambridge: The MIT Press; 2001.
43. Hanley JA, McNeil BJ. The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. 1982;143:29–36. [PubMed]
44. Milton Abramowitz M, Stegun IA, editors. Handbook of mathematical functions with formulas, graphs, and mathematical tables. New York: Dover; 1972.
45. Vainio MJ, Puranen JS, Johnson MS. ShaEP: molecular overlay based on shape and electrostatic potential. J Chem Inf Model. 2009;49:492–502. doi: 10.1021/ci800315d. [PubMed] [Cross Ref]
46. Brown N. Cheminformtics–an introduction for computer scientists. ACM Comput Surv. 2009;41:1–38. doi: 10.1145/1459352.1459353. [Cross Ref]

Articles from Springer Open Choice are provided here courtesy of Springer