- Research article
- Open Access
A theoretical entropy score as a single value to express inhibitor selectivity
BMC Bioinformatics volume 12, Article number: 94 (2011)
Designing maximally selective ligands that act on individual targets is the dominant paradigm in drug discovery. Poor selectivity can underlie toxicity and side effects in the clinic, and for this reason compound selectivity is increasingly monitored from very early on in the drug discovery process. To make sense of large amounts of profiling data, and to determine when a compound is sufficiently selective, there is a need for a proper quantitative measure of selectivity.
Here we propose a new theoretical entropy score that can be calculated from a set of IC50 data. In contrast to previous measures such as the 'selectivity score', Gini score, or partition index, the entropy score is non-arbitary, fully exploits IC50 data, and is not dependent on a reference enzyme. In addition, the entropy score gives the most robust values with data from different sources, because it is less sensitive to errors. We apply the new score to kinase and nuclear receptor profiling data, and to high-throughput screening data. In addition, through analyzing profiles of clinical compounds, we show quantitatively that a more selective kinase inhibitor is not necessarily more drug-like.
For quantifying selectivity from panel profiling, a theoretical entropy score is the best method. It is valuable for studying the molecular mechanisms of selectivity, and to steer compound progression in drug discovery programs.
In recent years, the kinase field has developed the practice of monitoring inhibitor selectivity through profiling on panels of biochemical assays [1–7], and other fields are following this example [8, 9]. Such profiling means that scientists are faced with increasing amounts of data that need to be distilled into human sense. It would be powerful to have a good single selectivity value for quantitatively steering the drug discovery process, for measuring progress of series within a program, for computational drug design [10–12], and for establishing when a compound is sufficiently selective. However, in contrast to, for instance, lipophilicity and potency, where values such as logP or binding constant (K d ) are guiding, quantitative measures for selectivity are still under debate. Often graphic methods are used to give insight, for example dotting a kinome tree [13, 14], heat maps [4, 6], or a radius plot, but such methods only allow qualitative comparison of a limited set of compounds at a time.
To make quantitative selectivity comparisons, three notable methods have been proposed (Figure 1). The first is the 'selectivity score' , which simply divides the number of kinases hit at an arbitrary K d or IC50 value (e.g. 3 μM) by the number of kinases tested (S(3 μM), Figure 1a). A related score is S(10x), which divides the number of kinases hit at 10 times the K d of the target by the number of kinases tested . The disadvantage of both methods is that 3 μM, or the factor 10, is an arbitrary cut-off value. For example, take two inhibitors, one that binds to two kinases with K d s of 1 nM and 1 μM, and another with K d s of 1 nM and 1 nM. Both are ranked equally specific by both S(3 μM) and S(10x), whereas the first compound is clearly more specific.
A less arbitrary parameter for selectivity is the Gini score . This uses %-inhibition data at a single inhibitor concentration. These data are rank-ordered, summed and normalized (Figure 1b) to arrive at a cumulative fraction inhibition plot, after which the score is calculated by the relative area outside the curve (Figure 1b). Though this solves the problem with the selectivity score, it leaves other disadvantages. One is that the Gini score has no conceptual or thermodynamic meaning such as a K d value has. Another is that it performs suboptimally with smaller profiling panels . In addition, the use of %-inhibition data makes the value more dependent on experimental conditions than a K d -based score . For instance, profiling with 1 μM inhibitor concentration results in higher percentages inhibition than using 0.1 μM of inhibitor. The 1 μM test therefore yields a more promiscuous Gini value, requiring the arbitrary 1 μM to be mentioned when calculating Gini scores. The same goes for concentrations of ATP or other co-factors. This is confusing and limits comparisons across profiles.
A recently proposed method is the partition index . This selects a reference kinase (usually the most potently hit one), and calculates the fraction of inhibitor molecules that would bind this kinase, in an imaginary pool of all panel kinases (Figure 1c). The partition index (Pmax) is a K d -based score with a thermodynamical underpinning, and performs well when test panels are smaller . However, this score is still not ideal, since it doesn't characterize the complete inhibitor distribution in the imaginary kinase mixture, but just the fraction bound to the reference enzyme. Consider two inhibitors: A binds to 11 kinases, one with a K d of 1 nM and ten others at 10 nM. Inhibitor B binds to 2 kinases, both with K d s of 1 nM. The partition index would score both inhibitors as equally specific (Pmax = 0.5), whereas the second is intuitively more specific. Another downside is the necessary choice of a reference kinase. If an inhibitor is relevant in two projects, it can have two different Pmax values. Moreover, because the score is relative to a particular kinase, the error on the K d of this reference kinase dominates the error in the partition index. Ideally, in panel profiling, the errors on all K d s are equally weighted.
Here we propose a novel selectivity metric without these disadvantages. Our method is based on the principle that, when confronted with multiple kinases, inhibitor molecules will assume a Boltzmann distribution over the various targets (Figure 1d). The broadness of this distribution can be assessed through a theoretical entropy calculation (it is not actually measuring entropy). We show the advantages of this method and some applications. Because it can be used with any activity profiling dataset, it is a universal parameter for expressing selectivity.
Results and discussion
Imagine a theoretical mixture of all protein targets on which selectivity was assessed. No competing factors are present such as ATP. To this mixture we add a small amount of inhibitor, in such a way that approximately all inhibitor molecules are bound by targets, and no particular binding site gets saturated. A selective inhibitor will bind to one target almost exclusively and have a narrow distribution (low entropy, Figure 1d). A promiscuous inhibitor will bind to many targets and have a broad distribution (high entropy, Figure 1d). The broadness of the inhibitor distribution on the target mixture reflects the selectivity of the compound.
The binding of one inhibitor molecule to a particular protein can be seen as a thermodynamical state with an energy level determined by K d (through ΔG = RTlnK d ). For simplicity we use the term K d to represent both K d and K i . The distribution of molecules over these energy states is given by the Boltzmann law. As the broadness of a Boltzmann distribution is measured by entropy, the selectivity implied in the distributions of Figure 1d can be captured in an entropy.
A similar insight is given by information theory. It is well-established that information can be quantified using entropy . A selective kinase inhibitor can be seen as containing more information about which active site to bind than a promiscuous inhibitor. The selectivity difference between the inhibitors can therefore be quantified by information entropy.
The distribution of a compound across energy states is given by the Boltzmann formula :
Where ϕ1 is the fraction of molecules occupying state 1, and ΔG1 is the free energy of occupying state 1 when the inhibitor comes from solution. In order to arrive at a fraction, the denominator in equation (1) contains the summation of occupancies of all states, which are labelled i, with free energies ΔG i .
In general, entropy can be calculated from fractions of all l states using the Gibbs formula :
Ssel is shorthand for selectivity entropy. Compared to the original Gibbs formulation, equation (2) contains a minus sign on the right hand to ensure that Ssel is a positive value. Now, we need to evaluate equation (2) from a set of measurements. For this we need
Where K a,i is the association constant of the inhibitor to target (or state) i, which is the inverse of the binding constant K d,i (which is a dissociation constant). In short: K a,i = 1/K d,i . If we express the free energy in units of 'per molecule' rather than 'per mole', equation (3) becomes
and equation (1) can be rewritten as
Using this result in equation (2) gives
Simplifying notation gives
Equation (7) defines how a selectivity entropy can be calculated from a collection of association constants K a . Here ΣK is the sum of all association constants.
It is most simple to apply equation (7) to directly measured binding constants or inhibition constants. Also IC50s can be used, but this is only really meaningful if they are related to K d . Fortunately, for kinases it is standard to measure IC50 values at [ATP] = KM,ATP. Ideally, such IC50s equal 2 times K d , according to the Cheng-Prusoff equation [19, 20]. The factor 2 will drop out in equation (7), and we therefore can use data of the format IC50-at-KM, ATP directly as if they were K d .
Protocol for calculating a selectivity entropy
From the above, it follows that a selectivity entropy can be quickly calculated from a set of profiling data with the following protocol:
Generate K a values by taking 1/K d or 1/IC50
Add all K a values to obtain ΣK
For every K a , calculate K a /ΣK
For every K a , evaluate (K a /ΣK) ln (K a /ΣK)
Sum all terms and multiply by -1
This process can be easily automated for use with large datasets  or internal databases.
The selectivity entropy is based on calculating the entropy of the hypothetical inhibitor distribution in a protein mixture. To give more insights into the properties of this metric, some examples are useful.
An inhibitor that only binds to a single kinase with a K d of 1 nM (K a = 109 M-1) has K a /ΣK a = 1. Then Ssel = -[1 ln 1]= 0, which is the lowest possibly entropy.
An inhibitor that binds to two kinases (X and Y) with a K d of 1 nM has K x /ΣK a = K y /ΣK a = 0.5 and a selectivity entropy of -[0.5 ln 0.5 + 0.5 ln 0.5] = 0.69. Thus lower selectivity results in higher entropy.
If we modify the compound such that it still inhibits kinase X with a K d of 1 nM, but inhibits less strongly kinase Y with a K d of 1 μM, then the new inhibitor is more specific. Now K x /ΣK a = 109/(109+106) and K y /ΣK a = 106/(109+106), resulting in Ssel = -[0.999 ln 0.999 + 0.001 ln 0.001] = 0.0079. This is less than 0.69. This shows that the selectivity entropy can distinguish in the case where the selectivity scores S(3 μM) and S(10x) cannot (see above).
A less selective inhibitor that binds three targets with K d s of 1 nM, has Ssel = -3·[0.3 ln 0.3 ] = 1.08, and an even more promiscuous inhibitor that binds 5 targets, of which 3 at 1 nM, and 2 at 1 μM, has ΣK = 3·109+ 2·106 = 3.002·109 and Ssel = -3·[1·109/3·109 ln 1·109/3·109 ] + 2·[1·106/3·109 ln 1·106/3·109 ] = 3.07. Thus Ssel gradually increases when more targets are more potently hit.
If we take the inhibitors A and B that were mentioned earlier, then A (with an inhibition profile of 1 nM, and ten times 10 nM), has ΣK = 1·109+ 10·108 = 2·109 and Ssel = - [1·109/2·109 ln 1·109/2·109 ] + 10·[1·108/2·109 ln 1·108/2·109 ] = 1.84. This is a more aselective value than inhibitor B with an inhibition profile of twice 1 nM, which has Ssel = 0.69 (see above). Thus the selectivity entropy can distinguish in a case where the partition coefficient Pmax cannot.
Comparison to other methods
Having defined the entropy, we next investigated its performance relative to the most widely-used methods, on a public profiling dataset of 38 inhibitors on 290 non-mutant kinases  (Table 1 and Additional file 1). The values for Gini score, S(3 μM), S(10x) and partition coefficient, were taken from earlier work . To this we added a K a -Gini value and the selectivity entropy. The K a -Gini is a Gini score directly calculated on K a s, without reverting to %-inhibition values (see below). From each of these scores we determined an inhibitor selectivity ranking, and a rank order difference compared to the entropy method (Uitdehaag_S1). In addition, to get an overview of the profiling raw data , we appended an activity-based heat map (Uitdehaag_S1).
From the rankings it is apparent that each of the earlier methods such as the classic Gini score, S(3 μM) and S(10x) generate considerable ranking differences compared to all other methods. This was observed earlier . For the Gini score, this is related to the conversion from IC50 to %-inhibition, because the K a -Gini gives more consistent rankings. For the S(3 μM) and the S(10x), the use of a cut-off is likely too coarse an approach. For instance in the case of S(10x), there are six inhibitors with a score of 0, making it impossible to distinguish between those highly specific compounds.
The newer methods such as Pmax, K a -Gini, and the selectivity entropy, give a more consistent ranking between them. For example, all three methods have PI-103, CI-1033, GW2580, VX-745 and gefitinib in their selectivity top five. There are differences however, most strikingly illustrated by the inhibitor SB-431542. This is ranked by Pmax as 31st most selective, but by K a -Gini and the selectivity entropy as 15th and 14th (Uitdehaag_S1). Also S(3 μM) ranks this ALK5 inhibitor  as selective. However, SB-431542 hits four kinases with very similar IC50s between 100-300 nM, which leads to a broad partitioning over these kinases, resulting in a very promiscuous Pmax of 0.14. The partition coefficient therefore ranks SB-431542 as almost equally selective to sunitinib (Pmax = 0.11, rank 33). Nevertheless, sunitinib inhibits 181 kinases below 3 μM, and SB-431542 only 5. Therefore we think that K a -Gini and the selectivity entropy are a better 'general' measure of selectivity in this case.
Another inhibitor scored differently is MLN-518 , which ranks 26st by Pmax, but 14th and 15th by K a -Gini and the selectivity entropy (Table 1 and Uitdehaag_S1). Again, these differences arise because this inhibitor hits 4 kinases with roughly equal potencies between 2-10 nM, leading to a promiscuous Pmax (0.26). However, MLN-518 only hits 10 kinases below 3 μM, making it intuitively more selective than e.g. ZD-6474  (Pmax = 0.28, ranked 25th by Pmax), which hits 79 kinases below 3 μM. These cases illustrate the earlier point that Pmax underscores inhibitors that only hit a few kinases at comparable potencies. The Gini score and selectivity entropy assign a higher selectivity to these cases.
Finally, any selectivity score should be in line with the visual ranking from a heat map. The Additional file 1 shows that, generally, compounds with a higher entropy indeed have a busier heat map. A few exceptions stand out, which by eye appear more promiscuous than their entropy ranking indicates, for instance SU-14813, sunitinib and staurosporin. However, these compounds have extreme low K d s on selected targets (SU-14813: 0.29 nM on PDGFRβ, sunitinib: 0.075 nM on PDGFRβ, staurosporin: 0.037 nM on LOK and 0.024 nM on SLK). Therefore they are relatively selective over activities in the 1-100 nM range, whereas these activities still fall within the highlighted ranges in Uitdehaag_S1. In a sense, the large dynamic range of the data limits visual assessment through a heat map.
Consistency across profiling methods
As a next step we selected 16 compounds from the public profile (Ambit) , and measured activity data on these using a different profiling service (Millipore, data available as Additional file 2). The 16 compounds represent a diversity of molecular scaffolds, promiscuity and target classes (Table 2). Also for these new data, we calculated the selectivity metrics (Uitdehaag_S2). In the ideal case, the selectivity values are similar irrespective of profiling technology (in the same way that a K d value is ideally independent of laboratory and assay format). The data of both methods are plotted in Figure 2.
All metrics except the entropy and Pmax tend to be quite unevenly distributed. For instance all K a -Gini scores fall between 0.93 and 1.00, where they can theoretically range from 0 to 1. If we nevertheless calculate the correlation statistics between both datasets, the R-square from linear regression and the correlation indicate that the selectivity entropy, S(3 μM) and K a -Gini are the most robust methods (Figure 2).
It would be ideal if the absolute value of the metrics could also be compared between datasets. This means that a specificity of e.g. 1.2 in the first profile, would also score 1.2 in the second profile. To get insight in this, we calculated the best fit to a 1:1 correlation (the diagonal line in Figure 2), using normalized data. The K a -Gini score was rescaled to its useful range of 0.93-1.00 (see legend to Figure 2), and then fitted. The S(3 μM) and the selectivity entropy have the best fit. The fact that here the K a -Gini performs poorer is probably caused by the use of cumulative inhibition values (Figure 1b), which leads to the accumulation of errors (as pointed out in ref. 16).
In all fits, the Pmax and S(10x) scores show worse fits and more scatter, indicating that these methods generate more error in their final value. For S(10x) and for Pmax, this is because both methods make use of a reference value, usually the most potent IC50, and errors in this reference value propagate more than errors in other IC50s. Ideally, for S(10x) and Pmax, the reference value specifically would have to be more accurately established.
If all analyses are taken together, the selectivity entropy avoids many pitfalls of the other methods (see above), shows consistent compound ranking (Table 1, Uitdehaag_S1), and is among the most robust methods across profiling datasets (Figure 2). For this reason, we propose the entropy method as the best metric for general selectivity.
Defining average selectivity
Quantification of selectivity helps to define when a compound is selective or promiscuous. Because of its consistency, the entropy method is ideally suited for benchmarking selectivity values. In the 290-kinase profiling dataset, the entropies are monomodally distributed, with an average of 1.8 (median of 1.9) and a standard deviation (σ) of 1.0 (not shown). Based on the correlation in Figure 2, it is expected that these statistics will be conserved in other profiling sets. Therefore, in general, a kinase compound with an entropy less than about 2 can be called selective, and more than 2 promiscuous. This provides a first quantitative definition of kinase selectivity.
Selectivity of allosteric inhibitors
It is generally thought that allosteric kinase inhibitors (known as type II, type III, or DFG-out inhibitors) are more selective [25, 26]. The selectivity entropy now allows quantitative testing of this idea. We identified, from literature, which inhibitors in the profiling datasets are type II and III, based on X-ray structures. Sorafenib induces the kinase DFG-out conformation in B-RAF , nilotinib and gleevec in Abl , GW-2580 in Fms  and BIRB-796 in p38α . Lapatinib induces a C-helix shift in EGFR . PD-0325901  and AZD-6244 induce a C-helix shift in MEK1 . All other kinase inhibitors in the profile were labelled type I. Comparing the entropy distributions in both samples shows that type II/III inhibitors have significantly lower entropies (Figure 3a). Although other factors, such as the time at which a compound was developed, could influence the entropy differences, the correlation between low entropy and allostery strongly supports the focus on allostery for developing specific inhibitors [25, 26].
Among the specific inhibitors in the type I category, 3D-structures of PI-103, CI-1033 and VX-745 bound to their targets have not been determined. Therefore, potentially, these inhibitors could also derive their specificity from a form of undiscovered induced fit. Indeed, VX-745-related compounds induce a peptide flip near Met109/Gly110 in P38α . Of the five most selective compounds in Table 1, only gefitinib so far is undoubtedly a type I inhibitor , making this EGFR inhibitor an interesting model for the structural biology of non-allosteric specificity.
Use of selectivity measures in nuclear receptor profiling
Selectivity profiling is most advanced in the kinase field, but is emerging in other fields. To illustrate that selectivity metrics such as the entropy can also be used with other target families, we investigated a long-standing question in the nuclear receptor field: are non-steroidal ligands more selective than steroidals? . For this, we calculated the entropies of a published profile of 35 antagonists on a panel of 6 steroid receptors  (the androgen receptor, estrogen receptor α, estrogen receptor β, mineralocorticoid receptor, glucocorticoid receptor, and progesterone receptor). This shows that there are no statistically significant selectivity differences between steroidals and non-steroidals (Figure 3b). A more important determinant for selectivity could be, in parallel to kinase inhibitors, if a ligand induces a conformational change. Indeed, many nuclear receptor agonists are known to induce a transformation from a flexible receptor to a rigid agonistic form [36–40], or a heterodimer form [41, 42]. In contrast, antagonists are know to displace helix 12 specifically from the agonistic form . Thus, the large role of induced fit in ligand binding to nuclear receptors might explain the relative high selectivity of these ligands [9, 36, 43, 44].
Use in hit prioritization
Aside from solving questions in the structure-function area, the selectivity entropy can be used during drug discovery. Previously it has been shown that selectivity metrics can be used in lead optimization projects to classify compounds, set targets, and rationalize improvement . In addition, metrics such as the entropy are useful in evaluating screening data, especially now screening larger compound collections in parallel assays is increasingly popular.
We downloaded PubChem data of 59 compounds tested in a panel of four assays for regulators of G protein signalling (RGS) . These data were selected because they were publicly available and were neither a kinase nor a nuclear receptor panel. In addition the data were dose-response, were all in a similar assay format, and were ran in the same lab with the same compound set.
We calculated the compound entropies across the RGS panel, and used them for ranking, which immediately distinguishes the scaffolds that are specific (Figure 3c). The best are ID 24785302, a pyrazole-phenoxy derivative, and ID 24834029, a bicyclo-octane derivative, which are likely to be better lead optimization starting points than more promiscuous scaffolds. Triaging compounds by entropy is a far more time-efficient and unbiased way than manual evaluation of four parallel columns of data. Indeed, listing of the selectivity entropy in public databases of screening data would provide users with immediate information on scaffold promiscuity.
Selectivity and clinical outcome
Finally, the selectivity entropy can be used to study clinical success. Selective compounds are generated because they are thought to be less toxic and therefore better doseable to effective ranges . To test the hypothesis that clinically approved inhibitors are more selective, we binned the compounds in the public kinase profile  according to their clinical history, and calculated their average entropies (Figure 3d, Additional file 3). Compared to the average discontinued compound, the average marketed kinase inhibitor is not more selective, and the average Phase III compound is even significantly more aselective. To exclude therapy area effects, we also performed the analysis for compounds in the oncology area, which is the only therapeutic area with a statistically significant amount of projects. This leads to a similar conclusion (Figure 3d). To exclude effects of time from this analysis (more recently invented kinase inhibitors might be more selective, because of advances in the kinase field), we repeated the analysis for compounds that entered clinical phase I before 2005. This shows even more clearly that more succesful compounds are, if anything, more broadly selective (Figure 3d).
Behind such statistics lies the success of, for instance, the spectrum selective drugs dasatinib, sorafenib and sunitinib (an average entropy of 3.13), and the failure of the highly selective MEK-targeted drugs PD-0325901 and CI-1040 (an average entropy of 0.32). Because 66-100% of the analysed compounds in each clinical bin are (or were) developed for oncology, our conclusion is primarily valid for oncology, until more kinase inhibitors enter the clinic for other indications. Nevertheless, the finding that a selective kinase inhibitor has fewer chances of surviving early clinical trials fuels the notion that polypharmacology is sometimes required to achieve effect (in oncology) [45–47].
In order to quantify compound selectivity as a single value, based on data from profiling in parallel assays, we have presented a selectivity entropy method, and compared this to other existing methods. The best method should avoid artifacts that obscure compound ranking, and show consistent values across profiling methods. Based on these criteria, the selectivity entropy is the best method.
A few cautionary notes are in order. First, the method is labelled an entropy in the sense of information theory , which is different to entropy in the sense of vibrational modes in enzyme active sites. Whereas these vibrations can form a physical basis for selectivity [39, 48, 49], our method is a computational metric to condense large datasets.
Secondly, any selectivity metric that produces a general value does not take into account the specific importance of individual targets. Therefore, the entropy is useful for generally characterizing tool compounds and drug candidates, but if particular targets need to be hit, or avoided, the K d s on these individual targets need to be monitored. It is possible to calculate an entropy on any particular panel of all-important targets, or to assign a weighing factor to every kinase, as suggested for Pmax and calculate a weighted entropy. However, the practicality of this needs to be assessed.
Next, it is good custom to perform profiling in biochemical assays at [ATP] = KM-ATP, because this generates IC50s that are directly related to the ATP-independent K d value. However, in a cellular environment, there is a constant high (~5 mM) ATP concentration and therefore a biochemically selective inhibitor will act with different specificity in a cell. If the inhibitor has a specificity for a target with a KM,ATP above the panel average, then that inhibitor will act even more specifically in a cell and vice versa (KM,ATP values can generally be found on websites of profiling research organizations). Selectivity inside the cell is also determined by factors such as cellular penetration, compartimentalization and metabolic activity . Therefore, selectivity from biochemical panel profiling is only a first step in developing selective inhibitors.
Another point is that any selectivity metric is always associated with the assay panel used, and the entropy value will change if an inhibited protein is added to the panel. Adding a protein that does not bind inhibitor will not affect the entropy value. In this way the discovery of new inhibitor targets by e.g. pulldown experiments, can change the idea of inhibitor selectivity, and also the entropy value. A good example is PI-103, the most selective inhibitor in Table 1, which in the literature is known as a dual PI3-kinase/mTOR inhibitor , and which appears specific in Table 1 because PI3-kinase is not incorporated in the profiling panel.
In addition, an inhibitor that hits 2 kinases at 1 nM from a panel of 10 has the same selectivity entropy as an inhibitor that inhibits 2 kinases at 1 nM in a panel of 100. However, intuitively, the second inhibitor is more specific (the 'selectivity score' differentiates in this case). This illustrates that it is important to compare entropy scores on similar panels. At the same time, when results from different panels are weighed, as in the example, it should not be assumed for the first inhibitor, that it is inactive against all 90 other kinases in the second panel. It would be better to assign an average K d where measurements are missing. In that case the first inhibitor would score a more promiscuous entropy compared to the second inhibitor.
Finally it must be stressed that the selectivity entropy could be applied in many more fields. It could, for instance, be a useful metric in the computational studies that attempt to link compound in vitro safety profiles to compound characteristics [51–53]. Currently, that field uses various forms of 'promiscuity scores' which bear similarity to the selectivity score. A more robust and non-arbitrary metric such as the selectivity entropy could be of help in building more detailed pharmacological models of compound activity-selectivity relationships [51–53].
In summary, the selectivity entropy is a very useful tool for making sense of large arrays of profiling data. We have demonstrated its use in characterizing tool compounds and drug candidates. Many more applications are imaginable in fields where an array of data is available and the selectivity of a response needs to be assessed. In that sense, the selectivity entropy is a general aid in the study of selectivity.
Calculation of other selectivity scores
For comparisons between currently used methods, we calculated the selectivity scores S(3 μM) and S(10x) as outlined above and in ref. 5. The partition coefficient Pmax was calculated as originally proposed , by taking the K a value of the most potently hit kinase, and dividing it by Σ K a . It is worth to note that the partition coefficient is the same as ϕl in our entropy equation (eq. 2).
The Gini score was calculated from data on %-inhibition . In Figure 1b, these data were extracted from K d values using the Hill expression: %-inhibition = 100/(1+10-(pKd - pconc)), where pK d = -log (K d ) and pconc = -log (inhibitor concentration evaluated). In addition, to work more directly with K d s, we also introduce a K a -Gini score, in which association constants are used for rank-ordering the kinase profile. From this K a -rank ordering, a cumulative effect is calculated and normalized, after which the areas are determined, in the same way as for the original Gini score . All calculations were done in Microsoft Excel.
Sources of existing and new data
For our comparative rank-ordering (Table 1, Uitdehaag_S1) we used the publicly available dataset released by Ambit http://www.ambitbio.com, which contains binding data (K d s) of 38 inhibitors on 290 kinases (excluding mutants), and which is currently the largest single profiling set available .
For comparing profiles across methods (Figure 2), we selected 16 kinase inhibitors of the Ambit profile (Table 2) and submitted these to the kinase profiling service from Millipore (http://www.millipore.com/drugdiscovery/svp3/kpservices, data available as Additional file 2). Both profiling methods are described earlier [3, 5, 14] and differ (among other variations) in the following way: Ambit uses a competitive binding setup in absence of ATP on kinases from T7 or HEK293 expression systems . Millipore uses a radioactive filter binding activity assay, with kinases purified from Escherichia coli or baculovirus expression systems . All Millipore profiling was done on 222 human kinases at [ATP] = KM,ATP.
For comparing inhibitors with an allosteric (actually: induced fit) profile (Figure 3a), we used data from the Ambit profile , supplemented with Millipore profiling data on nilotinib, PD-0325901 and AZD6244, because these important inhibitors were lacking in the Ambit dataset (data available in Additional file 2).
For comparing nuclear receptor data (Figure 3b), we used the published profiling dataset of 35 inhibitors on a panel consisting of all six steroid hormone receptors  The data we used were EC50s in cell-based assays.
For evaluation of a screening dataset (Figure 3c), we selected data from the PubChem initiative, determined at the University of New Mexico on regulators of G protein signalling (isoforms 4, 19, 7 and 16. Assay identifiers: 1872, 1884, 1888 and 1869) .
For evaluating clinical success (Figure 3d), we tracked the clinical status of each compound in the Ambit profile using the Thompson Pharma® database (status February 2011, analysis availabe as Additional file 3).
Davies SP, Reddy H, Caivano M, Cohen P: Specificity and mechanism of some commonly used protein kinase inhibitors. Biochem J 2000, 351(Pt1):95–105. 10.1042/0264-6021:3510095
Vieth M, Higgs RE, Robertson DH, Shapiro M, Gragg EA, Hemmerle H: Kinomics - structural biology and chemogenomics of kinase inhibitors and targets. Biochem Biophys Acta 2004, 1697(1–2):243–257.
Bain J, Plater L, Elliott M, Shapiro N, Hastie CJ, McLauchlan H, Klevernic I, Arthur JSC, Alessi DR, Cohen P: The selectivity of protein kinase inhibitors: a further update. Biochem J 2007, 408(3):297–315. 10.1042/BJ20070797
Fedorov O, Marsden B, Pogagic V, Rellos P, Müller S, Bullock AN, Schwaller J, Sundström M, Knapp S: A systematic interaction map of validated kinase inhibitors with Ser/Thr kinases. Proc Natl Acad Sci USA 2007, 104(51):20523–20528. 10.1073/pnas.0708800104
Karaman MW, Herrgard S, Treiber DK, Gallant P, Atteridge CE, Campbell BT, Chan KW, Ciceri P, Davis MI, Edeen PT, Faraoni R, Floyd M, Hunt JP, Lockhart DJ, Milanov ZV, Morrison MJ, Pallares G, Patel HK, Pritchard S, Wodicka LM, Zarrinkar PP: A quantitative analysis of kinase inhibitor selectivity. Nature Biotechnol 2008, 26(1):127–132. 10.1038/nbt1358
Bamborough P, Drewry D, Harper G, Smith GK, Schneider K: Assessment of chemical coverage of kinome space and its implications for kinase drug discovery. J Med Chem 2008, 51(24):7898–7914. 10.1021/jm8011036
Smyth LA, Collins I: Measuring and interpreting the selectivity of protein kinase inhibitors. J Chem Biol 2009, 2(3):131–151. 10.1007/s12154-009-0023-9
Heilker R, Wolff M, Tautermann CS, Bieler M: G-protein-coupled receptor focused drug discovery using a target class platform approach. Drug Discov Today 2009, 14(5–6):231–240. 10.1016/j.drudis.2008.11.011
Wilkinson JM, Hayes S, Thompson D, Whitney P, Bi K: Compound profiling using a panel of steroid hormone receptor cell-based assays. J Biomol Screen 2008, 13(8):755–765. 10.1177/1087057108322155
Sciabola S, Stanton RV, Wittkop S, Wildman S, Moshinsky D, Potluri S, Xi H: Predicting kinase selectivity profiles using free-wilson QSAR analysis. J Chem Inf Model 2008, 48(9):1851–1867. 10.1021/ci800138n
Sheridan RP, Nam K, Maiorov VN, McMasters DR, Cornell WD: QSAR models for predicting the similarity in binding profiles for pairs of protein kinases and the variation of models between experimental data sets. J Chem Inf Model 2009, 49(8):1974–1985. 10.1021/ci900176y
Brandt P, Jensen AJ, Nilsson J: Small kinase assay panels can provide a measure of selectivity. BioOrg Med Chem Letters 2009, 19(20):5861–5863. 10.1016/j.bmcl.2009.08.083
Manning G, Whyte DB, Martinez R, Hunter T, Sudarsanam S: The protein kinase complement of the human genome. Science 2002, 298(5600):1912–1934. 10.1126/science.1075762
Fabian MA, Biggs WH III, Treiber DK, Atteridge CE, Azimioara MD, Benedetti MG, Carter TA, Ciceri P, Edeen PT, Floyd M, Ford JM, Galvin M, Gerlach JL, Grotzfeld RM, Herrgard S, Insko DE, Insko MA, Lai AG, Lélias JM, Mehta SA, Milanov ZV, Velasco AM, Wodicka LM, Patel HK, Zarrinkar PP, Lockhart DJ: A small molecule-kinase interaction map for clinical kinase inhibitors. Nature Biotechnol 2005, 23(3):329–336.
Graczyk P: Gini coefficient: a new way to express kinase selectivity against a family of kinases. J Med Chem 2007, 50(23):5773–5779. 10.1021/jm070562u
Cheng AC, Eksterowicz J, Geuns-Meyer S, Sun Y: Analysis of kinase inhibitor selectivity using a thermodynamics-based partition index. J Med Chem 2010, 53(11):4502–4510. 10.1021/jm100301x
Shannon CE: A mathematical theory of communication. The Bell Systems Technical J 1948, 27: 379–423.
Atkins P, de Paula J: Atkins' Physical Chemistry. Oxford University Press: Oxford; 1970.
Cheng Y, Prusoff WH: Relationship between the inhibition constant (Ki) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction. Biochem Pharmacol 1973, 22(23):3099–3108. 10.1016/0006-2952(73)90196-2
Knight ZA, Shokat KM: Features of selective kinase inhibitors. Chem Biol 2005, 12(6):621–637. 10.1016/j.chembiol.2005.04.011
Roman DL, Talbot JN, Roof RA, Sunahara RK, Traynor JR, Neubig RR: Identification of small-molecule inhibitors of RGS4 using a high-troughput flow cytometry protein interaction assay. Mol Pharmacol 2006, 71(1):169–175. 10.1124/mol.106.028670
Inman GJ, Nicolás FJ, Callahan JF, Harling JD, Gaster LM, Reith AD, Laping NJ, Hill CS: SB-431542 is a potent and specific inhibitor of transforming growth factor-β superfamily type I activin receptor-like kinase (ALK) receptors ALK4, ALK5, and ALK7. Mol Pharmacol 2002, 62(1):65–74. 10.1124/mol.62.1.65
Kelly LM, Yu JC, Boulton CL, Apatira M, Li J, Sullivan CM, Williams I, Amaral SM, Curley DP, Duclos N, Neuberg D, Scarborough RM, Pandey A, Hollenbach S, Abe K, Lokker NA, Gilliland DG, Giese NA: CT53518, a novel selective Flt3 antagonist for the treatment of acute myelogenous leukemia (AML). Cancer Cell 2002, 1(5):421–432. 10.1016/S1535-6108(02)00070-3
Hennequin LF, Stokes ES, Thomas AP, Johnstone C, Plé PA, Ogilvie DJ, Dukes M, Wedge SR, Kendrew J, Curwen JO: Novel 4-anilinoquinazolines with C-7 basic side chains: design and structure activity relationship of a series of potent, orally active, VEGF receptor tyrosine kinase inhibitors. J Med Chem 2002, 45(6):1300–1312. 10.1021/jm011022e
Liu Y, Gray NS: Rational design of inhibitors that bind to inactive kinase conformations. Nature Chem Biol 2006, 2(7):358–364. 10.1038/nchembio799
Simard JR, Klüter S, Grütter C, Getlik M, Rabiller M, Rode HB, Rauh D: A new screening assay for allosteric inhibitors of Src. Nature Chem Biol 2009, 5(6):395–396. 10.1038/nchembio.162
Wan PTC, Garnett MJ, Roe SM, Lee S, Niculescu-Duvaz D, Good VM, Cancer Genome Project, Jones CM, Marshall CJ, Springer CJ, Barford D, Marais R: Mechanism of activation of the RAF-ERK signaling pathway by oncogenic mutations of B-RAF. Cell 2004, 116(6):855–867. 10.1016/S0092-8674(04)00215-6
Vajpai N, Strauss A, Fendrich G, Cowan-Jacob SW, Manley PW, Grzesiek S, Jahnke W: Solution conformations and dynamics of ABL kinase-inhibitor complexes determined by NMR substantiate the different binding modes of imatinib/nilotinib and dasatinib. J Biol Chem 2008, 283(26):18292–18302. 10.1074/jbc.M801337200
Meyers MJ, Pelc M, Kamtekar S, Day J, Poda GI, Hall MK, Michener ML, Reit BA, Mathis KJ, Pierce BS, Parikh MD, Mischke DA, Long SA, Parlow JJ, Anderson DR, Thorarensen A: Structure-based drug design enables conversion of a DFG-in binding CSF-1R kinase inhibitor to a DFG-out binding mode. BioOrg Med Chem Letters 2010, 20(5):1543–1547. 10.1016/j.bmcl.2010.01.078
Pargellis C, Tong L, Churchill L, Cirillo PF, Gilmore T, Graham AG, Grob PM, Hickey ER, Moss N, Pav S, Regan J: Inhibition of p38 MAP kinase by utilizing a novel allosteric binding site. Nature Struct Biol 2002, 9(4):268–272. 10.1038/nsb770
Wood ER, Truesdale AT, McDonald OB, Yuan D, Hassell A, Dickerson SH, Ellis B, Pennisi C, Horne E, Lackey K, Alligood KJ, Rusnak DW, Gilmer TM, Shewchuk L: A unique structure of epidermal growth factor receptor bound to GW572016 (lapatinib): relationships among protein conformation, inhibitor off-rate, and receptor activity in tumor cells. Cancer Res 2004, 64(18):6652–6659. 10.1158/0008-5472.CAN-04-1168
Fischmann TO, Smith CK, Mayhood TW, Meyers JE Jr, Reichert P, Mannarino A, Carr D, Zhu H, Wong J, Yang RS, Le HV, Madison VS: Crystal structures of MEK1 binary and ternary complexes with nucleotides and inhibitors. Biochemistry 2009, 48(12):2661–2674. 10.1021/bi801898e
Fitzgerald CE, Patel SB, Becker JW, Cameron PM, Zaller D, Pikounis VB, O'Keefe SJ, Scapin G: Structural basis for p38α MAP kinase quinazolinone and pyridol-pyrimidine inhibitor specificity. Nature Struct Biol 2003, 10(9):764–769. 10.1038/nsb949
Johnson LN: Protein kinase inhibitors: contributions from structure to clinical compounds. Q Rev Biophysics 2009, 42(1):1–40. 10.1017/S0033583508004745
Hermkens PH, Kamp S, Lusher S, Veeneman GH: Non-steroidal steroid receptor modulators. IDrugs 2006, 9(7):488–494.
Egea PF, Klaholz BP, Moras D: Ligand-protein interaction in nuclear receptors of hormones. FEBS Letters 2000, 476(1–2):62–67. 10.1016/S0014-5793(00)01672-0
Nettles KW, Sun J, Radek JT, Sheng S, Rodriguez L, Katzenellenbogen JA, Katzenellenbogen BS, Greene GL: Allosteric control of ligand selectivity between estrogen receptors α and β. Mol Cell 2004, 13(3):317–327. 10.1016/S1097-2765(04)00054-1
Raaijmakers HC, Versteegh JE, Uitdehaag JCM: The X-ray structure of RU486 bound to the progesterone receptor in a destabilized agonistic conformation. J Biol Chem 2009, 284(9):19572–19579. 10.1074/jbc.M109.007872
Martinez L, Nascimento AS, Nunes FM, Phillips K, Aparicio R, Dias SM, Figueira AC, Lin JH, Nguyen P, Apriletti JW, Neves FAR, Baxter JD, Webb P, Skaf MS, Polikarpov I: Gaining ligand selectivity in thyroid hormone receptors via entropy. Proc Natl Acad Sci USA 2009, 106(49):20717–20722. 10.1073/pnas.0911024106
Togashi M, Borngraeber S, Sandler B, Fletterick RJ, Webb P, Baxter JD: Conformational adaptation of nuclear receptor ligand binding domains to agonists: potential for novel approaches to ligand design. J Steroid Biochem Mol Biol 2005, 93(2–5):127–137. 10.1016/j.jsbmb.2005.01.004
Bourguet W, Vivat V, Wurtz JM, Chambon P, Gronemeyer H, Moras D: Crystal structure of a heterodimeric complex of RAR and RXR ligand-binding domains. Mol Cell 2000, 5(2):289–298. 10.1016/S1097-2765(00)80424-4
Fradera X, Vu D, Nimz O, Skene R, Hosfield D, Wynands R, Cooke AJ, Haunsø A, King A, Bennett DJ, McGuire R, Uitdehaag JCM: X-ray structures of the LXRα LBD in its homodimeric form and implication for heterodimer signaling. J Mol Biol 2010, 399(1):120–132. 10.1016/j.jmb.2010.04.005
Cornell W, Nam K: Steroid hormone binding receptors: application of homology modeling, induced fit docking, and molecular dynamics to study structure function relationships. Curr Top Med Chem 2009, 9(9):844–853. 10.2174/156802609789207109
Nabuurs SB, Wagener M, de Vlieg J: A flexible approach to induced fit docking. J Med Chem 2007, 50(26):6507–6518. 10.1021/jm070593p
Knight ZA, Lin H, Shokat KM: Targeting the cancer kinome through polypharmacology. Nature Rev Cancer 2010, 10(2):130–137. 10.1038/nrc2787
Morphy R, Rankovic Z: Designing multiple ligands - medicinal chemistry strategies and challenges. Curr Pharm Design 2009, 15(6):587–600. 10.2174/138161209787315594
Hopkins A: Network pharmacology: the next paradigm in drug discovery. Nature Chem Biol 2008, 4(11):682–690. 10.1038/nchembio.118
Robinson D, Sherman W, Farid R: Understanding kinase selectivity through energetic analysis of binding site waters. Chem Med Chem 2010, 5(4):618–627.
Scapin G: Protein kinase inhibition: different approaches to selective inhibitor design. Curr Drug Targets 2006, 7(11):1443–1454.
Raynaud FI, Eccles S, Clarke PA, Hayes A, Nutley B, Alix S, Henley A, Di-Stefano F, Ahmad Z, Guillard S, Bjerke LM, Kelland L, Valenti M, Patterson L, Gowan S, de Haven Brandon A, Hayakawa M, Kaizawa H, Koizumi T, Ohishi T, Patel S, Saghir N, Parker P, Waterfield M, Workman P: Pharmacologic characterization of a potent inhibitor of class I phosphatidylinositide 3-kinases. Cancer Res 2007, 67(12):5840–5850. 10.1158/0008-5472.CAN-06-4615
Yang Y, Chen H, Nilsson I, Muresan S, Engkvist O: Investigation of the relationship between topology and selectivity for druglike molecules. J Med Chem 2010, 53(21):7709–7714. 10.1021/jm1008456
Peters JU, Schnider P, Mattei P, Kansy M: Pharmacological promiscuity: dependence on compound properties and target specificity in a set of recent Roche compounds. Chem Med Chem 2009, 4(4):680–686.
Azzaoui K, Hamon J, Faller B, Whitebread S, Jacoby E, Bender A, Jenkins JL, Urban L: Modelling promiscuity based on in vitro safety pharmacology data. Chem Med Chem 2007, 2(6):874–880.
We thank our colleagues Rogier Buijsman, Husam Alwan and Koen Dechering for discussion and critical reading of the manuscript. We thank Jennifer Wilkinson (Invitrogen) for sharing nuclear receptor data.
JU conceived the entropy principle and drafted the manuscript. GZ organized the kinase profiling data and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Selectivity metrics, heat maps, and selectivity rank ordering for the Ambit profiling dataset (an extension of Table
Additional file 1:1). (TIFF 4 MB)
Additional file 2:. (PDF 364 KB)
Additional file 3:. (PDF 129 KB)
About this article
Cite this article
Uitdehaag, J.C., Zaman, G.J. A theoretical entropy score as a single value to express inhibitor selectivity. BMC Bioinformatics 12, 94 (2011). https://doi.org/10.1186/1471-2105-12-94