- Research article
- Open Access
A physicochemical descriptor-based scoring scheme for effective and rapid filtering of kinase-like chemical space
Journal of Cheminformaticsvolume 4, Article number: 4 (2012)
The current chemical space of known small molecules is estimated to exceed 1060 structures. Though the largest physical compound repositories contain only a few tens of millions of unique compounds, virtual screening of databases of this size is still difficult. In recent years, the application of physicochemical descriptor-based profiling, such as Lipinski's rule-of-five for drug-likeness and Oprea's criteria of lead-likeness, as early stage filters in drug discovery has gained widespread acceptance. In the current study, we outline a kinase-likeness scoring function based on known kinase inhibitors.
The method employs a collection of 22,615 known kinase inhibitors from the ChEMBL database. A kinase-likeness score is computed using statistical analysis of nine key physicochemical descriptors for these inhibitors. Based on this score, the kinase-likeness of four publicly and commercially available databases, i.e., National Cancer Institute database (NCI), the Natural Products database (NPD), the National Institute of Health's Molecular Libraries Small Molecule Repository (MLSMR), and the World Drug Index (WDI) database, is analyzed. Three of these databases, i.e., NCI, NPD, and MLSMR are frequently used in the virtual screening of kinase inhibitors, while the fourth WDI database is for comparison since it covers a wide range of known chemical space. Based on the kinase-likeness score, a kinase-focused library is also developed and tested against three different kinase targets selected from three different branches of the human kinome tree.
Our proposed methodology is one of the first that explores how the narrow chemical space of kinase inhibitors and its relevant physicochemical information can be utilized to build kinase-focused libraries and prioritize pre-existing compound databases for screening. We have shown that focused libraries generated by filtering compounds using the kinase-likeness score have, on average, better docking scores than an equivalent number of randomly selected compounds. Beyond library design, our findings also impact the broader efforts to identify kinase inhibitors by screening pre-existing compound libraries. Currently, the NCI library is the most commonly used database for screening kinase inhibitors. Our research suggests that other libraries, such as MLSMR, are more kinase-like and should be given priority in kinase screenings.
Chemical space can be defined as "the total descriptor space covered by all the known and possible small organic compounds" . Chemical space is thus so vast that it prompted Lipinsky and Hopkins to compare it to the total number of stars in the cosmos . Estimates of the total number of possible small molecules vary from 108 to 10200 depending upon the criteria used. For example, Bohacek et al.  estimated it to be 1060, when based on a maximum number of 30 C, N, O, and S atoms; Ertl  estimated a total of 1020-1024 possible small molecules, based on current synthetic methods; and Ogata et al.  estimated a range of 108-1019 possible small molecules, based on combinations of known Protein Data Bank (PDB) ligands.
The CAS registry  is the largest collection of disclosed chemical substance information and currently contains more than 55 million organic and inorganic compounds. Other notable collections of compounds include the Chemical Structure Lookup Service (CSLS) , with around 46 million unique compounds, PubChem  and Chemspider , with around 20 million compounds each, and ZINC  with around 13 million compounds, along with hundreds of other public or private collections ranging from a few thousands to a few millions of compounds. Even though such vast collections only constitute a small fraction of possible chemical space, it is still very difficult to apply a typical biological screen to all molecules in a collection when seeking novel hits on targets of interest . Along with database size, another concern is that very few compounds in these databases are biologically relevant; in other words, the sub-regions of chemical space that are relevant to biology is small [1, 12]. Since not every region of chemical space defined by a compound database is biologically relevant, screening the entire database for a particular target is a waste of resources. In recent years, the focus has shifted away from screening large compound libraries to screening smaller, more target-focused libraries that are generated using all relevant information about the target and its known active compounds [13–17].
The design of focused libraries using physicochemical-based descriptors is known as chemography. The underlying principle of this technique is that structurally similar compounds are likely to have similar interactions with associated targets, along with having similar physicochemical property ranges [18–22]. Such profiling of compounds based on physicochemical descriptors has been in use since the late 1990's and many excellent research articles on this concept exist [23–34]. The most popular methods are the rules defining drug-likeness proposed by Lipinski et al.  and more recently by Veber et al.  and Oprea et al. [11, 29, 36, 37]. These rules are based on simple physicochemical descriptors such as molecular weight, number of hydrogen bond donors and acceptors, logP, polar surface area, and number of rotatable bonds. Since their publication, these rules have been extensively used to differentiate between drugs, lead-like compounds, and other compounds, and have also been used as filters to reduce the size of screening databases. Ideally, these rules must be based on individual target-based known small molecule exemplars. Previously such rules have been applied in a few target classes like G protein-coupled receptors (GPCRs) and kinases [38, 39]. The ultimate goal of such applications is to maximize the diversity of hits while minimizing the size of screening libraries to cover only the biologically relevant space for a particular drug target.
Our focus in this work is to apply physicochemical descriptor-based filters to make kinase-likeness rules. Previously, to make a kinase-focused library, two main types of approaches have been used: 1) scaffold-based design, [40–45] or 2) 2D and 3D (pharmacophore) similarity-based design [46–48]. Both approaches utilize information from known kinase active and inactive compounds. While these approaches have been successfully used in test cases, their utility in finding novel active compounds is limited due to fact that the new compounds selected are often structurally similar to the compounds that are initially used as templates. In comparison, our method, which is based on physicochemical descriptors of known active compounds, has the advantage of being able to select novel structural scaffolds since the physicochemical properties are not directly tied to any specific chemical scaffold. The method is also easy to implement and can handle very large compound data sets.
Protein kinases, which are important signaling enzymes that catalyze a phosphoryl transfer from ATP to a given protein or peptide substrate, are the second most popular drug targets after GPCRs and account for roughly one third of all drug discovery projects . There are currently more than 518 distinct human protein kinase genes spread across seven kinase families . Each kinase contains a conserved catalytic domain consisting of N-terminal and C-terminal lobes connected by a hinge region . The ATP-binding site is located in a deep cleft between these lobes. Numerous structural studies of kinases have shown that, although there is substantial diversity in sequences, the structure of the ATP-binding site  is relatively conserved. Over the years, thousands of inhibitors (actives) have been identified for a variety of kinases. Most of these actives compete with ATP for the ATP-binding site and are classified as ATP-mimetics. The conserved nature of the binding site architecture among kinases suggests that kinase actives occupy a common, biologically relevant region of chemical space that is compatible with the structural and chemical constraints of the ATP-binding site topology. This region can be roughly described by the common steric and electronic features, i.e., physicochemical properties, among known ATP-mimetic kinase actives [39, 52]. A rigorous, high-throughput, statistical approach for identifying these features presents the opportunity to design kinase-focused libraries which consist of compounds that contain complementary features to the common ATP-mimetic actives.
To this end, we collected a set of 22,615 known kinase actives from the ChEMBL database . This database was assembled from more than 1,200 published research articles. We then tested a range of physicochemical descriptors and identified those that significantly distinguish between kinase and non-kinase active compounds. The measured distributions of these properties from kinase actives was then used to train a kinase-likeness scoring function. We used this scoring function to rank order the following small molecule databases for their kinase-likeness: The National Cancer Institute database (NCI),  the Natural Products database (NPD),  and the National Institute of Health's Molecular Libraries Small Molecule Repository (MLSMR),  all of which are commonly used in screening for novel kinase actives. We further validated our results by rank ordering the same databases using another correlation matrix-based scoring function originally developed by Ohno et al. .
Finally, we used the kinase-likeness scoring function to generate a kinase-focused library. In order to demonstrate the utility of this approach, we employed a docking-based screen of this focused library against three diverse and well-known kinase targets to see if the screening results are enriched for tighter binding compounds relative to standard screening libraries. We found that the top compounds from the kinase-focused library indeed made better contacts with the kinase binding pockets than the top compounds from standard screening libraries.
Results show that our kinase scoring function, based on physicochemical descriptors, successfully generates kinase-focused libraries. As a result, application of this scoring function can be a useful first step in kinase actives lead discovery workflows.
Materials and methods
Experimentally validated kinase actives were compiled from the ChEMBL compound database, hereafter referred to as the kinase binding database (KBD). A compound was defined as a kinase active only if its Ki or IC50 value was less than 10 μM. To test kinase-likeness of the compounds that are frequently used in virtual screening of kinase inhibitors, the following publicly available small-molecule databases were also collected: the National Cancer Institute database (NCI), the Natural Products database (NPD), the National Institute of Health's Molecular Libraries Small Molecule Repository (MLSMR). Additionally, for comparison, compounds from the World Drug Index (WDI) database were also collected since they cover a wide range of known chemical space. Both the NCI and the NPD databases were downloaded from the ZINC UCSF collection  of compounds, while the MLSMR compounds were downloaded from the PubChem database,  and the WDI was obtained from Thomson Reuters . The data preprocessing and descriptor calculations were done with Chemical Computing Group's Molecular Operating Environment (MOE), v. 2010.10,  and implemented through Accelrys Pipeline Pilot, v. 8.0. . The preprocessing steps included removing counterions and solvent molecules, setting protonation states and removing metal-based structures. Molecules not in the weight range of 100-700 daltons were also removed from the databases. The final prepared databases numbered 22,615 (KBD), 264,554 (NCI), 103,668 (NPD), 343,605 (MLSMR), and 61,368 (WDI) compounds.
Properties previously used by Oprea et al.  in formulating lead-likeness were used as the starting point. These properties included molecular weight (Wt), octanol/water partition coefficient (SlogP), topological polar surface area (TPSA), and number of rotatable bonds (RB) (bonds were considered rotatable if they a) had a bond order of 1, b) were not a ring and c) had at least two heavy neighbors). Additionally, the number of rings (Rings), hydrogen bond acceptors (HBA), hydrogen bond donors (HBD), number of nitrogen atoms (nN), and number of oxygen atoms (nO) were also included.
Kinase-likeness Scoring Functions
For each KBD physicochemical property, there is a trend towards a normal distribution (details in Results section). Therefore, we utilized a multi-parameter Gaussian scoring function, called the Kinase-Like Score (KLS), for "kinase-likeness". This function is given by:
The values X i , with i = 1 to 9, represent the physicochemical properties of Wt, RB, HBA, HBD, SlogP, TPSA, Rings, nN, and nO, respectively. The parameters μi and σi represent the average and standard deviation of each physicochemical property, X i . The averages and standard deviations were calculated using property distributions obtained from the KBD. In addition, for any compound that fell outside the P10-P90 range (the descriptor range that covers 80% of the compounds) in any of its descriptors, a value of zero was assigned for that descriptor's contribution to the overall score. The KLS score ranges from 0 to 9, and compounds with scores toward the upper end of this range are more kinase-like than compounds with scores toward the lower end of this range. This scoring scheme shares similarities to that published by Wager et al. . In their work, multi-parameter scoring functions were derived for central nervous system compounds using the properties ClogP, ClogD, Wt, TPSA, HBD, and pKa. More elaborate schemes are found in papers such as that by Oashi et al.  that describe a multivariate Gaussian probability distribution to ascertain bioavailability and drug likeness of small molecule compounds.
We calculated the KLS for all of the compounds in the NCI, NPD, MLSMR, and WDI databases. High KLS values of compounds in these libraries indicate that these libraries are enriched with compounds having similar physicochemical properties as those of known kinase actives and ultimately may yield a higher enrichment of kinase actives in screening.
Correlation Coefficient-based Scoring
In addition to the KLS, we applied another orthogonal statistical approach based on a correlation coefficient-based methodology to rank various databases for kinase-likeness. A similar approach has already been applied by Ohno et al. to rank order various test databases for drug-likeness . In this methodology we first computed the correlation coefficient between every pair of descriptors in our descriptor set using Pipeline Pilot software. This was done for the NPD, NCI, MLSMR, WDI, and KBD databases. Nine descriptor values were considered, therefore a total of 36 correlation coefficient values were determined for each database. In theory, if the correlation matrices of two compound libraries are similar, compounds in each library would be similar too. Therefore, a compound library with a correlation matrix similar to the KBD should include kinase-like molecules.
In this study, a scoring function of kinase-likeness based on a comparison of correlation matrices (KLSC) was defined as follows :
In this equation, C i, j (test database) and C i, j (kinase database) represent the correlation coefficient values between two descriptors i and j in the test database and kinase database, respectively. Finally, we calculated a KLSC value for each of the databases, NPD, NCI, and MLSMR with respect to the KBD.
Receiver Operating Characteristic (ROC) Plots
To determine the kinase-likeness of any given database (test database) of compounds, we randomly chose 1% of the known kinase compounds from the KBD and seeded them into the NCI, NPD, MLMSR, and WDI databases. The KLS was computed for each compound in the seeded databases. The compounds were then sorted in decreasing order by their KLS value, and a ROC plot was constructed based on this ordering. The step of mixing a random selection of kinase actives from the KBD was repeated 10 times, and the average, maximum, and minimum ROC curves and area under the curve (AUC) were computed. For validation, a random 1% of kinase actives were also extracted from the KBD and mixed back with the KBD database to make ROC plots. This strategy should result in a ROC curve corresponding to a random selection of compounds from the database. To determine if a different percentage of seeded kinase compounds from the KBD gives different results, we tried other percentages as well (2%, 5%, 10%, 50%, and 100%). We found that the AUC value remains more or less independent of percentage of seeded kinase actives. As a result, only a selection of 1% was used.
Focused Library Construction, Target Selection, and Glide Docking
We applied the KLS (Equation 1) to make a target-focused library for kinases. For this, we initially pooled the major databases (NCI, NPD, and MLSMR) together. From this pool of 711,827 compounds, all the known kinase actives and duplicates were removed and a kinase-focused library was constructed by selecting the top 10,000 compounds based on their kinase-likeness score (KLS). This set is called the focused set. Another test set of 10,000 randomly selected compounds was also selected and called the random set. To test enrichment, docking-based virtual screening (VS) was performed against three kinase targets.
For docking calculations, the receptor and ligand structures were prepared using the Protein Preparation Workflow implemented in the Maestro graphical user interface of the Schrodinger program suite. The various steps involved in this process included removing water molecules, protonating receptor and ligand, assigning optimal tautomeric states of histidine residues, flipping side chains of glutamine and asparagine residues if required, optimizing hydrogen bonds, and energy minimization. The native ligands, along with compounds from the focused set and random set, were processed using LigPrep  with Epik  to expand the protonation and tautomeric states in the range of 7.0 ± 2.0 pH units. Computational docking calculations were done using GLIDE in standard precision (SP) mode  for flexible docking. The binding affinities were estimated as a Docking score (Glide score + Epik state penalties). The docking protocol was validated by re-docking the native ligands, each of which reproduced the binding mode observed in the respective crystal structures. Both the focused-set and random-set compounds were docked, and the top hits were selected using a docking score cutoff value equal to the docking score of the native ligand in each crystal structure. The number of kinase ligands with respect to the total number of ligands in the resulting hit list is an indication of library enrichment.
Additionally, to show that the compounds in the focused library constructed using the physicochemical property based filtering are more diverse than scaffold-based or similarity-based focused libraries, we conducted a test. We first selected the Aldrich Market Select (CNC-AMS) library available through ChemNavigator's iResearch library . This library contains ~ 7 million commercial compounds, all of which are readily purchasable. To test for diversity of focused libraries, we used an example of an already approved kinase drug Sunitinib. Sunitinib is a small-molecule targeting receptor tyrosine kinase and is approved by the FDA for the treatment of both renal cell carcinoma and the imatinib-resistant gastrointestinal stromal tumor . Using Sunitinib's chemical structure, we created three types of focused libraries. First, we created a scaffold-based library by collecting all the compounds from CNC-AMS that contain the same cyclic scaffold as Sunitinib. Second, we created a similarity-based library by collecting the top 100 most similar compounds to drug Sunitinib from CNC-AMS based on Tanimoto similarity score. Third, we created a physiochemical property -based focused library by applying our KLS score and selecting the top 100 KLS scoring compounds from the CNC-AMS library. To evaluate the diversity of each library, we calculated the average of distance, where distance is defined as '1 - similarity', for every pair of collected compounds in all the three libraries. Extended-connectivity fingerprints_6 (ECFP_6) were used to calculate the similarity . For comparison, we also calculated the average distance of sets of 100 randomly selected molecules from the CNC-AMS database.
Studies based on Lipinski's rule-of-five  and Oprea's lead-like criteria  have established that specific sets of properties, such as size, solubility, flexibility, polarity, and hydrogen bond capacities, are sufficient to capture significant molecular features and, in turn, the biological properties of small molecules. In this study, we used nine physicochemical descriptors that influence the above-described properties to effectively quantify the kinase-likeness of a given small molecule. Out of these nine descriptors, six (Wt, RB, HBA, HBD, TPSA, and SlogP) were used by Oprea . In the case of kinase targets, we have prior knowledge of binding site architecture, so three additional descriptors (Rings, nN and nO) were added. These were included because planar ring scaffolds are known to fit well in the ATP binding pocket, and nitrogen and oxygen atoms play a large role in making hydrogen bond contacts.
We compared the distribution of these descriptors for kinase actives (KBD database) with those of molecules from the NPD, NCI, MLSMR, and WDI databases. Comparison was done using a two-sided Student's t test  to determine if the difference in descriptor distributions between two databases was significant. The analysis showed that, with the exception of the number of RB between the KBD and MLSMR databases (a p value of 0.528), all other descriptors were significantly different (p < 0.05). Hence, these nine different descriptors are suitable for formulating a scoring function that can differentiate between small molecules that are likely to bind kinases (actives that have kinase-likeness) and those that are not. The distribution plots of the nine descriptors, along with the outliers, average values, standard deviation (SD) is shown in Figure 1 and Table 1. A complete table with included P10-P90 (range covering 80% of the compounds) values is provided in the Additional File 1.
The statistical analysis of all descriptors shows that, with the exception of nO, the average values of all the other descriptors, i.e., (Wt, RB, HBA, HBD, TPSA, SlogP, Rings, and nN) were higher for kinase active compounds than for the compounds contained in the WDI, NPD, NCI, and MLSMR databases. Of note, the average RB value of the MLSMR compounds and average TPSA value of the WDI compounds were similar to that of the kinase active compounds. Kinase actives had the lowest average value for nO while compounds from the WDI database had the highest nO value. This analysis showed that, on average, the kinase active compounds were heavier (more weight), more flexible (more RB), more polar (more nN, HBA, and HBD) and more hydrophobic (more SlogP and rings) than the compounds in the other databases.
We utilized the distinctive properties of kinase active compounds to develop a kinase-likeness scoring function. Based on the distribution of the descriptor values of kinase actives from the KBD, we calculated a kinase-likeness score, KLS, for each small molecule in the KBD, WDI, NPD, NCI, and MLSMR databases (Equation 1). This score reflects how similar each compound is with respect to the physicochemical properties of known kinase actives. To calculate how kinase-like these databases were, we generated an ROC plot (as described in the Methods section) based on the sorted KLS. This is shown in Figure 2. If a small random set of kinase actives is seeded into the KBD, and VS is performed using the KLS to calculate an ROC AUC, this should produce an AUC value near 0.5. Our results showed this to be 0.506 (upper left plot in Figure 2). This means that the KLS is unable to distinguish the 1% mixed kinase active compounds from other kinase actives in the KBD; in other words, the compounds in the KBD are kinase-like. This is an obvious result but provides a good way to test if databases other than the KBD contain a large number of kinase-like compounds. If similar screening experiments with databases other than the KBD result in computed AUC values near 0.5, these databases would be deemed kinase-like. This in turn, indicates that these databases are composed of compounds that primarily have physicochemical properties similar to that of kinase active compounds. Using these ROC plots and associated AUC values, we can determine which databases are kinase-like. Here we focused only on the NCI, NPD, and MLSMR databases because of their common usage in screening for kinase targets (Figure 2). The WDI database was used as a comparison based on the diversity of compounds within it.
Table 2 shows the p values computed using a two-sided Student's t test for comparison of AUC values between the KBD and the other databases (WDI, NPD, NCI, and MLSMR). Note that all p values were < 0.05, indicating that the AUC distributions for the NCI, WDI, NPD, and MLMSR databases were significantly different than the AUC distribution for the KBD.
From our ROC analysis, based on AUC values, we divided these databases into high, medium, and low kinase-like categories. High kinase-likeness databases have AUC values near 0.5. This is the case for the MLSMR database with an AUC value of 0.57. Low kinase-likeness databases have AUC values closer to 1.0. This is the case for the NCI (AUC = 0.82) and the WDI (AUC = 0.81) databases. Medium kinase-likeness databases have AUC values near 0.75. This is the case for the NPD database (AUC = 0.70). The high AUC result for WDI is expected since it contains marketed and in-development small molecule drugs that bind to a wide variety of protein targets. We believe that the lower AUC value for the MLSMR and the comparatively higher AUC value for the NCI database is a consequence of the more recent assembly of the MLSMR database compared to the NCI database. The MLSMR initiative started in 2005, whereas the NCI has been operational since the late 1950's. Research on kinase inhibitors has only ramped up in recent decades and consequently MLSMR, being the newer of the two databases, contains more kinase active compounds than the NCI. The NPD sits right in the middle with an AUC of 0.7. The hypothesis is that although NPD represents a hugely diverse set of natural product compounds (due to the large size range of the compounds), some portions of NPD have been optimized for kinase-likeness. Therefore, the NPD AUC value is not quite as high as the AUC values of the WDI or MLSMR databases (because of the kinase-like portions), but it is not as low as the KBD either (because of the diversity).
Using in silico rules prior to screening can help focus efforts on databases appropriate to protein targets of interest. The results presented here suggest that, when screening for potential kinase actives, it would be most efficient to probe the databases in the following order: MLSMR > NPD > WDI ≈ NCI.
To support our AUC-based results for the ranking of VS databases, we further analyzed them using the KLSC score (Equation 2), which is based on the similarity between the correlation matrix of a test database and the KBD. The calculated correlation matrix of nine physicochemical descriptors is shown in Table 3 for KBD compounds.
The correlation matrices of the other databases (WDI, NPD, NCI, and MLSMR) are provided in the Additional File 2. Ohno et al. has recently applied a similar scoring function, called the Drug Like Score (DLS), to compare commercially available compound libraries, clinical candidate libraries, and several other chemical libraries with a launched drug library to calculate drug likeness based on six physicochemical descriptors .
The calculated KLSC values for the NCI, WDI, NPD, and MLSMR databases are 0.114, 0.161, 0.079, and 0.066, respectively. A lower KLSC value signifies higher kinase-likeness for a database. The KLSC values shown above reiterate what we found before, i.e., MLSMR contain the most kinase-like active compounds of all the databases studied. Overall, the ranking of kinase-likeness based on KLSC values is MLSMR > NPD > WDI ≈ NCI. This is similar to that observed in ROC-based AUC analysis described in the previous section and shown in Table 2. An additional statistical test based on Kolmogorov-Smirnov methodology (KS test)  is shown in Additional File 3 along with the distribution of KLS values for all the studied databases. The k value obtained by KS test further confirmed the ranking of databases obtained by both the AUC based ranking and the KLSC value based ranking.
Application in Docking-based Virtual Screening
With the aim of demonstrating the practicalities of the kinase-likeness scoring function, two datasets were created: a kinase-focused set and a random set. Both these datasets contain 10,000 small molecules (details in Methods section). To test which of these two datasets was more enriched in kinase-binding compounds, we performed VS of each against three different kinase targets selected from three different branches (families) of the human kinome tree. The criteria for selecting these kinase drug targets was that, 1) they should have been used in previous structure-based VS studies, 2) they should have been crystallized with reasonable resolution (≤ 3.0Å) with bound ligand, and 3) each should be from a different kinase family. With these criteria, we selected B-Raf to represent the tyrosine kinase-like (TKL) family, mitogen-activated protein kinase-activated protein kinase 2 (MK2) to represent the Ca2+/calmodulin-dependent protein kinases (CaMK) family, and cyclin-dependent kinase 2 (CDK2) to represent the CMGC kinase family. The PDB codes for these proteins are 2FB8 , 2PZY , and 1OIT  respectively.
A docking-based VS, as explained in the Methods section, was performed using both the kinase-focused compound set and the random compound set. The outcomes were analyzed in terms of number of "hit" compounds, defined as compounds that scored better than the docked native ligand in each respective crystal structure of the selected targets. The results are shown in Figure 3. In each case, three to four times more hits were obtained from the focused compound set than from the random compound set. Of note, many of the top hit compounds from the focused set yielded a much better docking score than the top hits from the random set. These results demonstrate that the compounds from VS of the focused set are not only more enriched for kinase-like compounds than the random screen, but also interact better with the kinase binding pockets.
To show that physicochemical property-based focused libraries are more diverse than structure-based focused libraries, we created three focused libraries based on the kinase drug Sunitinib: scaffold-based, similarity-based, and physicochemical properties-based. All three focused libraries were extracted from the CNC-AMS database of ~ 7 million purchasable compounds. In the first case, only the cyclic scaffold of drug Sunitinib was taken into account when performing scaffold screening of CNC-AMS. A total of 44 compounds were found to contain similar cyclic scaffold. In the second case, the full structure of Sunitinib was used for screening and the top 100 similar compounds, based on Tanimoto similarity value, were selected. In the third case, the KLS scoring function was applied on all the compounds of CNC-AMS and top 100 scoring compounds were selected. The average distance (1-similarity) between each pair of molecules was 0.69, 0.74, and 0.85, for scaffold-base, similarity-based, and the physicochemical property-based (KLS score) focused libraries, respectively. For comparison a collection of 100 randomly select compounds from CNC-AMS shows average distance of 0.88. This proves that the physicochemical property-based focused library is more diverse than structure-based focused libraries. In fact, the average distance value of 0.85 for the physicochemical property-based focused library generated from CNC-AMS is as diverse as a collection of random selected compounds from CNC-AMS.
Kinases represent one of the most important and therapeutically relevant drug target families involved in diseases like cancer and various other neurodegenerative, infectious and inflammatory diseases. In the past few years, various computational methodologies have been employed to facilitate the kinase drug discovery effort through the screening of small compound databases. Many of these methodologies focus on steps that reduce the size of screening databases. For example, the most common types of pre-screening filters are based on selective and preferential kinase scaffold-based filtering and 2D- and 3D-based similarity filtering using the templates of known kinase active compounds. In the work presented here, we demonstrated another easily calculable and amenable statistical methodology for the generation of virtual kinase-focused libraries.
This methodology is based on the observation that various physicochemical properties of known kinase active compounds such as molecular weight, polar surface area, and rotatable bonds are unique for this structural target class and can be used to differentiate kinase active compounds from molecules of other structural classes. This method calculates a kinase-likeness score based on property distribution ranges of nine physicochemical descriptors derived from known kinase active compounds. These physicochemical properties are molecular weight, SlogP, TPSA, rotatable bonds, hydrogen bond acceptors and donors, number of rings, number of nitrogen and oxygen atoms. Previous studies of making kinase-focused libraries use scaffold-based and similarity-based filters, which are primarily inclusive filters. These methodologies obtain enrichment by including molecules containing scaffolds identified in known actives. Such methods are popular, but the collected compounds lack structural diversity, since they are all based on either the same or similar known scaffolds. By comparison, our physicochemical filtering scheme achieves enrichment by including molecules that exhibit certain desirable physicochemical properties. Since this method relies solely on the physicochemical properties, any selection done using it becomes independent of known scaffolds. Hence the possibility of finding novel scaffolds that are structurally unrelated to known kinase actives is greatly increased. Using an example of kinase drug Sunitinib, we showed that focused libraries based on physicochemical property screening are more diverse than those generated by structure based techniques.
Our method can be used to prioritize small molecule databases for kinase-based virtual screening. Subsequent to database prioritization, the method can be used to constrain the size of screening databases and create focused libraries. The method can even be used to guide the de novo design of molecules prior to synthesis and testing. Extension of this method to improve the quality of other focused libraries is easily possible; the scoring function is readily modifiable since new physicochemical parameters may be added and less useful ones may be deleted.
Using our kinase-likeness score and a validated correlation matrix-based score for comparison, we showed that the MLSMR database compounds are much closer to kinase actives, in terms of physicochemical properties, than compounds in the NCI and NPD databases. We also developed and tested a kinase-focused library based on our kinase-likeness score and showed its superior enrichment for three specific kinase targets from different branches of the human kinome, relative to random screening of a collection of compounds. In each case, the focused library showed three to four times more enrichment compared to random screening of compounds. Given the demonstrated performance of our kinase-like scoring function for generating enriched libraries, we suggest that it may be useful as a first step in kinase-targeted virtual screening workflows.
Dobson CM: Chemical space and biology. Nature. 2004, 432: 824-828. 10.1038/nature03192.
Lipinski C, Hopkins A: Navigating chemical space for biology and medicine. Nature. 2004, 432: 855-861. 10.1038/nature03193.
Bohacek RS, McMartin C, Guida WC: The art and practice of structure-based drug design: a molecular modeling perspective. Med Res Rev. 1996, 16: 3-50. 10.1002/(SICI)1098-1128(199601)16:1<3::AID-MED1>3.0.CO;2-6.
Ertl P: Cheminformatics analysis of organic substituents: identification of the most common substituents, calculation of substituent properties, and automatic identification of drug-like bioisosteric groups. J Chem Inf Comput Sci. 2003, 43: 374-380. 10.1021/ci0255782.
Ogata K, Isomura T, Yamashita H, Kubodera H: A Quantitative Approach to the Estimation of Chemical Space from a Given Geometry by the Combination of Atomic Species. QSAR & Combinatorial Science. 2007, 26: 596-607. 10.1002/qsar.200630037.
Chemical Abstracts Service. [http://www.cas.org]
Chemical Structure Lookup Service. [http://cactus.nci.nih.gov/cgi-bin/lookup/search]
Hann MM, Oprea TI: Pursuing the leadlikeness concept in pharmaceutical research. Curr Opin Chem Biol. 2004, 8: 255-263. 10.1016/j.cbpa.2004.04.003.
Koch MA, Schuffenhauer A, Scheck M, Wetzel S, Casaulta M, Odermatt A, Ertl P, Waldmann H: Charting biologically relevant chemical space: a structural classification of natural products (SCONP). Proc Natl Acad Sci USA. 2005, 102: 17272-17277. 10.1073/pnas.0503647102.
Chen G, Zheng S, Luo X, Shen J, Zhu W, Liu H, Gui C, Zhang J, Zheng M, Puah CM, et al: Focused combinatorial library design based on structural diversity, druglikeness and binding affinity score. J Comb Chem. 2005, 7: 398-406. 10.1021/cc049866h.
Mok NY, Brenk R: Mining the ChEMBL Database: An Efficient Chemoinformatics Workflow for Assembling an Ion Channel-Focused Screening Library. J Chem Inf Model. 2011, 51: 2449-2454. 10.1021/ci200260t.
Deng Z, Chuaqui C, Singh J: Knowledge-based design of target-focused libraries using protein-ligand interaction constraints. J Med Chem. 2006, 49: 490-500. 10.1021/jm050381x.
Schnur D, Beno BR, Good A, Tebben A: Approaches to target class combinatorial library design. Methods Mol Biol. 2004, 275: 355-378. 10.1385/1-59259-802-1:355.
Schreyer SK, Parker CN, Maggiora GM: Data shaving: a focused screening approach. J Chem Inf Comput Sci. 2004, 44: 470-479. 10.1021/ci030025s.
Di L, Kerns EH: Profiling drug-like properties in discovery research. Curr Opin Chem Biol. 2003, 7: 402-408. 10.1016/S1367-5931(03)00055-3.
Horio K, Muta H, Goto J, Hirayama N: A simple method to improve the odds in finding 'lead-like' compounds from chemical libraries. Chem Pharm Bull (Tokyo). 2007, 55: 980-984. 10.1248/cpb.55.980.
Khanna V, Ranganathan S: Physiochemical property space distribution among human metabolites, drugs and toxins. BMC Bioinformatics. 2009, 10 (Suppl 15): S10-10.1186/1471-2105-10-S15-S10.
Raevsky OA, Trepalin SV, Trepalina HP, Gerasimenko VA, Raevskaja OE: SLIPPER-2001 -- software for predicting molecular properties on the basis of physicochemical descriptors and structural similarity. J Chem Inf Comput Sci. 2002, 42: 540-549. 10.1021/ci010097o.
van de Waterbeemd H, Gifford E: ADMET in silico modelling: towards prediction paradise?. Nat Rev Drug Discov. 2003, 2: 192-204. 10.1038/nrd1032.
Ajay A, Walters WP, Murcko MA: Can we learn to distinguish between "drug-like" and "nondrug-like" molecules?. J Med Chem. 1998, 41: 3314-3324. 10.1021/jm970666c.
Egan WJ, Walters WP, Murcko MA: Guiding molecules towards drug-likeness. Curr Opin Drug Discov Devel. 2002, 5: 540-549.
Lajiness MS, Vieth M, Erickson J: Molecular properties that influence oral drug-like behavior. Curr Opin Drug Discov Devel. 2004, 7: 470-477.
Lipinski CA: Drug-like properties and the causes of poor solubility and poor permeability. J Pharmacol Toxicol Methods. 2000, 44: 235-249. 10.1016/S1056-8719(00)00107-6.
Muegge I: Selection criteria for drug-like compounds. Med Res Rev. 2003, 23: 302-321. 10.1002/med.10041.
Muegge I, Heald SL, Brittelli D: Simple selection criteria for drug-like chemical matter. J Med Chem. 2001, 44: 1841-1846. 10.1021/jm015507e.
Oprea TI, Gottfries J: Chemography: the art of navigating in chemical space. J Comb Chem. 2001, 3: 157-166. 10.1021/cc0000388.
Proudfoot JR: Drugs, leads, and drug-likeness: an analysis of some recently launched drugs. Bioorg Med Chem Lett. 2002, 12: 1647-1650. 10.1016/S0960-894X(02)00244-5.
Veber DF, Johnson SR, Cheng HY, Smith BR, Ward KW, Kopple KD: Molecular properties that influence the oral bioavailability of drug candidates. J Med Chem. 2002, 45: 2615-2623. 10.1021/jm020017n.
Walters WP, Murcko MA: Recognizing molecules with drug-like properties. Curr Opin Chem Biol. 1999, 3: 384-387. 10.1016/S1367-5931(99)80058-1.
Walters WP, Murcko MA: Prediction of 'drug-likeness'. Adv Drug Deliv Rev. 2002, 54: 255-271. 10.1016/S0169-409X(02)00003-0.
Wang J, Ramnarayan K: Toward designing drug-like libraries: a novel computational approach for prediction of drug feasibility of compounds. J Comb Chem. 1999, 1: 524-533. 10.1021/cc990032m.
Lipinski CA, Lombardo F, Dominy BW, Feeney PJ: Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Advanced Drug Delivery Reviews. 1997, 23: 3-25. 10.1016/S0169-409X(96)00423-1.
Oprea TI: Current trends in lead discovery: are we looking for the appropriate properties?. Mol Divers. 2002, 5: 199-208.
Oprea TI, Davis AM, Teague SJ, Leeson PD: Is there a difference between leads and drugs? A historical perspective. J Chem Inf Comput Sci. 2001, 41: 1308-1315. 10.1021/ci010366a.
Balakin KV, Tkachenko SE, Lang SA, Okun I, Ivashchenko AA, Savchuk NP: Property-based design of GPCR-targeted library. J Chem Inf Comput Sci. 2002, 42: 1332-1342. 10.1021/ci025538y.
Zhang J, Yang PL, Gray NS: Targeting cancer with small molecule kinase inhibitors. Nat Rev Cancer. 2009, 9: 28-39. 10.1038/nrc2559.
Aronov AM, McClain B, Moody CS, Murcko MA: Kinase-likeness and kinase-privileged fragments: toward virtual polypharmacology. J Med Chem. 2008, 51: 1214-1222. 10.1021/jm701021b.
Brenk R, Schipani A, James D, Krasowski A, Gilbert IH, Frearson J, Wyatt PG: Lessons learnt from assembling screening libraries for drug discovery for neglected diseases. ChemMedChem. 2008, 3: 435-444. 10.1002/cmdc.200700139.
Keri G, Szekelyhidi Z, Banhegyi P, Varga Z, Hegymegi-Barakonyi B, Szantai-Kis C, Hafenbradl D, Klebl B, Muller G, Ullrich A, et al: Drug discovery in the kinase inhibitory field using the Nested Chemical Library technology. Assay Drug Dev Technol. 2005, 3: 543-551. 10.1089/adt.2005.3.543.
Lowrie JF, Delisle RK, Hobbs DW, Diller DJ: The different strategies for designing GPCR and kinase targeted libraries. Comb Chem High Throughput Screen. 2004, 7: 495-510.
Xi H, Lunney EA: The design, annotation, and application of a kinase-targeted library. Methods Mol Biol. 2011, 685: 279-291. 10.1007/978-1-60761-931-4_14.
Zhang C, Bollag G: Scaffold-based design of kinase inhibitors for cancer therapy. Curr Opin Genet Dev. 2010, 20: 79-86. 10.1016/j.gde.2009.10.010.
Deanda F, Stewart EL, Reno MJ, Drewry DH: Kinase-targeted library design through the application of the PharmPrint methodology. J Chem Inf Model. 2008, 48: 2395-2403. 10.1021/ci800276t.
Decornez H, Gulyas-Forro A, Papp A, Szabo M, Sarmay G, Hajdu I, Cseh S, Dorman G, Kitchen DB: Design, selection, and evaluation of a general kinase-focused library. ChemMedChem. 2009, 4: 1273-1278. 10.1002/cmdc.200900164.
Sun D, Chuaqui C, Deng Z, Bowes S, Chin D, Singh J, Cullen P, Hankins G, Lee WC, Donnelly J, et al: A kinase-focused compound collection: compilation and screening strategy. Chem Biol Drug Des. 2006, 67: 385-394. 10.1111/j.1747-0285.2006.00396.x.
Cohen P: Protein kinases--the major drug targets of the twenty-first century?. Nat Rev Drug Discov. 2002, 1: 309-315. 10.1038/nrd773.
Manning G, Whyte DB, Martinez R, Hunter T, Sudarsanam S: The protein kinase complement of the human genome. Science. 2002, 298: 1912-1934. 10.1126/science.1075762.
Bossemeyer D: Protein kinases--structure and function. FEBS Lett. 1995, 369: 57-61. 10.1016/0014-5793(95)00580-3.
Milletti F, Vulpetti A: Predicting polypharmacology by binding site similarity: from kinases to the protein universe. J Chem Inf Model. 2010, 50: 1418-1431. 10.1021/ci1001263.
National Cancer Institute. [http://dtp.nci.nih.gov/webdata.html]
Austin CP, Brady LS, Insel TR, Collins FS: NIH Molecular Libraries Initiative. Science. 2004, 306: 1138-1139. 10.1126/science.1105511.
Ohno K, Nagahara Y, Tsunoyama K, Orita M: Are there differences between launched drugs, clinical candidates, and commercially available compounds?. J Chem Inf Model. 2010, 50: 815-821. 10.1021/ci100023s.
Irwin JJ, Shoichet BK: ZINC--a free database of commercially available compounds for virtual screening. J Chem Inf Model. 2005, 45: 177-182. 10.1021/ci049714+.
Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Bryant SH: PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009, 37: W623-633. 10.1093/nar/gkp456.
World Drug Index. 1999, Dervent Information: London, UK
The Molecular Operating Environment Version 2010.10, Chemical Computing Group Inc., 1010 Sherbrooke Street West, Suite 910, Montreal, Canada H3A 2R7, 2010. [http://www.chemcomp.com]
Pipeline Pilot, version 8.0. 2010, Accelrys San Diego, California, USA
Wager TT, Chandrasekaran RY, Hou XT, Verhoest PR, Villalobos A, Will Y: Defining Desirable Central Nervous System Drug Space through the Alignment of Molecular Properties, in Vitro ADME, and Safety Attributes. ACS Chemical Neuroscience. 2010, 1: 420-434. 10.1021/cn100007x.
Oashi T, Ringer AL, Raman EP, Mackerell AD: Automated selection of compounds with physicochemical properties to maximize bioavailability and druglikeness. J Chem Inf Model. 2011, 51: 148-158. 10.1021/ci100359a.
LigPrep, version 2.4.107. 2010, Schrodinger, LLC, New York, NY
Epik, version 2.1.209. 2010, Schrodinger, LLC, New York, NY
Glide, version 5.6.109. 2010, Schrodinger, LLC, New York, NY
O'Farrell AM, Abrams TJ, Yuen HA, Ngai TJ, Louie SG, Yee KW, Wong LM, Hong W, Lee LB, Town A, et al: SU11248 is a novel FLT3 tyrosine kinase inhibitor with potent activity in vitro and in vivo. Blood. 2003, 101: 3597-3605. 10.1182/blood-2002-07-2307.
Rogers D, Hahn M: Extended-connectivity fingerprints. J Chem Inf Model. 2010, 50: 742-754. 10.1021/ci100050t.
Bland JM: An Introduction To Medical Statistics. 1995, Oxford: Oxford University Press, 2
Corder GW, Foreman DI: Nonparametric statistics for non-statisticians: a step-by-step approach. 2009, New York: Wiley
King AJ, Patrick DR, Batorsky RS, Ho ML, Do HT, Zhang SY, Kumar R, Rusnak DW, Takle AK, Wilson DM, et al: Demonstration of a genetic therapeutic index for tumors expressing oncogenic BRAF by the kinase inhibitor SB-590885. Cancer Res. 2006, 66: 11100-11105. 10.1158/0008-5472.CAN-06-2554.
Wu JP, Wang J, Abeywardane A, Andersen D, Emmanuel M, Gautschi E, Goldberg DR, Kashem MA, Lukas S, Mao W, et al: The discovery of carboline analogs as potent MAPKAP-K2 inhibitors. Bioorg Med Chem Lett. 2007, 17: 4664-4669. 10.1016/j.bmcl.2007.05.101.
Anderson M, Beattie JF, Breault GA, Breed J, Byth KF, Culshaw JD, Ellston RP, Green S, Minshull CA, Norman RA, et al: Imidazo[1,2-a]pyridines: a potent and selective class of cyclin-dependent kinase inhibitors identified through structure-based hybridisation. Bioorg Med Chem Lett. 2003, 13: 3021-3026. 10.1016/S0960-894X(03)00638-3.
Funding of this research was provided by the U.S. Department of Defense Threat Reduction Agency Grant TMTI0004_09_BH_T. The opinions or assertions contained herein are the private views of the authors and are not to be construed as official or as reflecting the views of the U.S. Army or of the U.S. Department of Defense.
The authors declare that they have no competing interests.
NS implemented the methods and conducted the experiments with assistance from HS, SC, MDMA, AW and GT. All five participated in study's design, coordination and manuscript drafting, and gave the final approval of the version to be submitted.