This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.
Perspective Cite This: ACS Omega 2019, 4, 2758−2765
http://pubs.acs.org/journal/acsodf
Recent Progress in Structure-Based Evaluation of Compound Promiscuity Erik Gilberg†,‡ and Jürgen Bajorath*,† †
Downloaded via 146.185.201.209 on February 8, 2019 at 06:05:40 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.
Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Endenicher Allee 19c, D-53115 Bonn, Germany ‡ Pharmaceutical Institute, Rheinische Friedrich-Wilhelms-Universität, An der Immenburg 4, D-53121 Bonn, Germany ABSTRACT: In chemistry and drug discovery, compound promiscuity is a controversial issue and often viewed differently. On the one hand, promiscuity has a clear-cut negative connotation, as promiscuous behavior of small molecules is often associated with nonspecific binding or assay artifacts. On the other hand, it is also well-established that compounds frequently interact with multiple targets including distantly related or unrelated proteins. In drug discovery, multitarget activity of small molecules receives increasing attention because it is therapeutically relevant, giving rise to desired or undesired pharmacological effects. Exploration of compound promiscuity is far from being simple because true and artificial activities are often difficult to discern. In light of these complications, X-ray structures of ligand−target complexes provide a wealth of information about molecular promiscuity, which is just beginning to be recognized and explored. Systematic analysis of structurally confirmed binding events involving artifact-prone compounds or multitarget ligands eliminates some of the uncertainties associated with promiscuity analysis and puts it on a new level, enabling the study of promiscuity-conferring interactions at the molecular level of detail. Herein, we discuss recent progress made in structure-based promiscuity analysis and put key findings into scientific context. inhibitor drug development in oncology19,20 and is also emerging in other therapeutic areas.16−18 However, polypharmacology is also critically viewed in pharmaceutical research because it represents a major departure from the single-target specificity paradigm,21,22 which has long governed discovery efforts. The extent to which drugs and other bioactive compounds might be promiscuous also continues to be debated and views often differ here as well. Several estimates and expectation values have been put forward. For example, on the basis of drug−target network analysis, it has been proposed that drugs might interact with on average 3−13 targets,23,24 depending on drug classes and data sets used. Data incompleteness certainly affects promiscuity estimates23 because not all drugs have been tested in vitro against all possible targets (and likely never will be). Hence, not surprisingly, it has been presumed that drugs should generally be promiscuous.14,23 However, currently available activity data do not support expectation values claiming the presence of widespread and strong promiscuity among drugs and other bioactive compounds.25−27 There are several lines of evidence to balance such expectations. For example, even most extensively tested screening compounds are on average only active against two to three targets (with a median value of two targets) and the majority of compounds
1. INTRODUCTION Apparent multitarget activity of small molecules, also referred to as promiscuity,1 is a controversial issue in chemistry and pharmaceutical research for a variety of reasons. First of all, it is often difficult to draw a line between “true” promiscuity and apparent promiscuity resulting from assay artifacts.2 Such artifacts might be caused by aggregators3−5 or other molecules with assay interference potential such as pan-assay interference compounds (PAINS),6−8 which may have a variety of undesired reactivities under assay conditions. PAINS typically occur as substructures in larger compounds. Assay interference effects are often difficult to recognize. To further complicate matters, even notorious PAINS do not display consistent activity patterns. Compounds containing such PAINS substructures might be frequently active in screening assays, rarely active, or consistently inactive, depending on how these substructures are embedded.9,10 Moreover, some well-known assay interference compounds also have a long history as privileged structures in medicinal chemistry.11,12 Hence, distinguishing target-based activity and promiscuous ligand-binding events from artificial activity signals is far from being a trivial task.13 On the other hand, true multitarget activity of compounds is of increasing relevance in drug discovery because it provides the basis of polypharmacology.14−18 The concept of polypharmacology implies that drugs often elicit their therapeutic effects through in vivo interactions with multiple targets. Polypharmacology has become a major theme for kinase © 2019 American Chemical Society
Received: December 27, 2018 Accepted: January 29, 2019 Published: February 6, 2019 2758
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765
ACS Omega
Perspective
are consistently inactive.25 Moreover, comprehensive analysis of high-confidence activity data has revealed that the majority of bioactive compounds from medicinal chemistry including drugs are only annotated with a single target and that promiscuous compounds are rare.26,27 Hence, further largescale experimental assessment of single- vs multitarget activity will be required to assess the influence of data incompleteness, advance our understanding of molecular promiscuity, and evaluate preconceived notions of polypharmacological drug actions that cannot be generalized at present. It is confirmed, however, that confined subsets of highly promiscuous compounds and drugs exist19,26 and that multitarget engagement is highly relevant for a number of therapeutic applications.18,19 In addition, multitarget activity also plays a pivotal role for repurposing of existing drugs.28,29 Given the topical nature of molecular promiscuity and the uncertainties associated with its assessment, a variety of computational approaches have been introduced to study promiscuity patterns and/or predict promiscuous compounds. For example, promiscuity was analyzed computationally through systematic compound data mining9,10,25−27,30 and ligand−target network modeling.14,24 In addition, statistically balanced ligand-based similarity searching was employed to infer new targets of active compounds,29,31 algorithms were introduced to detect promiscuity patterns at the level of molecular core structures,32 and machine learning models were developed to systematically distinguish between promiscuous and inactive interference compounds.33 Furthermore, promiscuity was also studied at the target structure level. In particular, protein binding site similarity was analyzed to explain multitarget activity of bound ligands.34,35 In addition, drugs with multitarget activity were assigned to proteins having similar functions on the basis of X-ray structures.36 Hence, promiscuity as the basis of polypharmacology has been, and continues to be, investigated from different perspectives but rationalizing multitarget activity of small molecules is far from being simplistic.
ity analysis was to determine whether or not a sufficient number of X-ray structures containing promiscuous ligands were available.
3. SEARCHING FOR X-RAY STRUCTURES WITH PROMISCUOUS LIGANDS Therefore, we searched the RCSB Protein Data Bank (PDB)37 and its Ligand Expo section38 for complex structures with bound promiscuous compounds. This was required to identify ligands that were complexed with different proteins and classify these targets. Activity measurements for ligands were retrieved from PDBbind39 and ChEMBL40 databases. Crystallographic target proteins were assigned to families on the basis of UniProt identifiers41 and the ChEMBL target classification scheme. In addition, for qualifying crystallographic ligands, a computational search for structural analogues42 was carried out in ChEMBL and in the PDB. To prioritize template compounds from medicinal chemistry for the design of ligands with multitarget activity (see below), X-ray ligands were mapped to ChEMBL and structural analogues available in ChEMBL were identified. We initially searched for X-ray complexes that contained compounds with PAINS substructures. In this case, the initial goal was to determine how frequently PAINS that tend to be reactive under assay conditions might form specific ligand− target interactions. Then, PAINS-containing complex structures were analyzed in detail. Our systematic search identified an unexpectedly large number of 2874 X-ray structures containing 1107 compounds with 70 different PAINS substructures, including many wellknown assay interference motifs.43 Thus, a large number of PAINS motifs were found in ligand−target complexes, revealing their ability to engage in specific interactions, despite assay liabilities. These structures provided a basis for investigating binding of PAINS vs potential reactivity in a variety of protein environments,43 as further discussed below. Subsequently, we identified crystallographic ligands bound to multiple targets (termed multitarget ligands) and targets from different families (multifamily ligands). Such ligands were considered to represent prime examples of promiscuous compounds with confirmed target-based interactions. Multifamily ligands were of particular interest because they bound to distantly related or unrelated targets and hence provided test cases for exploring how specific interactions were formed in distinct binding sites. Again, an unexpectedly large number of 1418 crystallographic multitarget ligands (with molecular weight >300 Da) were identified, which contained 702 multifamily ligands.44 These multifamily ligands were active against a median of three targets from two families. For 355 ligands available in ChEMBL, the median target number increased to 17 (but the family median remained at 2). For 168 of the 702 multifamily ligands, series of analogues (comprising a total of 4829 compounds covering 190 additional targets) were identified in ChEMBL, which yielded 133 different analogue series-based (ASB) scaffolds45 as templates for polypharmacology-directed compound design. “Polypharmacology by design” represents an extension of the polypharmacology concept and another emerging theme in drug discovery.46−48 Multifamily ligands we identified also included a subset of 192 compounds for which one or more potency values of at least 10 μM or lower (pIC50, pKi, or pKd ≥ 5) were reported in
2. PROMISCUOUS COMPOUNDS AND THE QUEST FOR CONFIRMATORY EVIDENCE Given the uncertainties associated with interpreting compound promiscuity on the basis of assay data, the question becomes how to best evaluate multitarget binding events. In our view, the most conservative and least error-prone approach to exploring compound promiscuity is the analysis of X-ray structures of ligand−target complexes with promiscuous compounds. Of course, X-ray data also have intrinsic limitations. For example, details of ligand−target interactions might not always be fully discernible due to limited resolution or other crystallographic ambiguities. Moreover, X-ray structures of ligand−target complexes provide a time-average static view of interactions that might well be dynamically driven and as such difficult to reconcile. However, in light of possible assay artifacts, a major advantage of focusing promiscuity assessment on X-ray structures is that these structures confirm binding of a chemical entity to a given target. Simply put, we can “see” electron density of a bound ligand and know a compound is active against the target. In addition, complex X-ray structures enable the study of multitarget activity at the molecular level of detail. A major limitation of structure-based promiscuity analysis is the sparseness of X-ray structures relative to biological assay data. Thus, an essential prerequisite for meaningful promiscu2759
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765
ACS Omega
Perspective
compound databases. Hence, in these instances, activity measurements corroborated structural evidence of specific binding events. This subset of multifamily ligands was further analyzed to rationalize and compare multifamily interactions in detail.49 The 192 ligands were represented by a total of 3378 X-ray structures that involved proteins from 2 to 16 different families per ligand. The multifamily ligands were found to contain 40 endogenous compounds and 10 others that preferentially bound to metabolizing enzymes or serum proteins. For each of the remaining 142 ligands, the targets were compared in detail to identify ligands that interacted with structurally or functionally related proteins assigned to different families, similar binding domains in different proteins, or distinct targets. On the basis of this analysis, a subset of 91 multifamily ligands was identified that exclusively bound to unrelated proteins. This subset was assigned highest priority for promiscuity analysis and their complexes were studied in detail.49 Taken together, analyses of complex X-ray structures yielded an unexpectedly large knowledge base of multitarget and multifamily ligands, which provided a sound basis for structurebased promiscuity analysis. In the following, key findings from our studies43,44,49 are discussed using representative examples.
4. STRUCTURE-BASED COMPOUND PROMISCUITY ANALYSIS 4.1. Assay Interference Compounds: Specific Interactions vs Artifacts. Figure 1 shows different complexes formed by assay interference compounds, which illustrate their dual ability to engage in specific ligand−target interactions or cause artifacts in assays. Figure 1a shows quercetin in complex with two unrelated targets including a kinase and a Niquercetinase, an oxidoreductase. Quercetin contains a catechol substructure, a prominent PAINS motif known to cause assay artifacts due to undesired metal chelation, redox activity, reactivity against biological nucleophiles in its oxidized form, or membrane perturbation.5,7,50,51 Catechol substructures are not only found in synthetic compounds but also in a variety of natural products. Quercetin’s catechol moiety is part of a polyphenolic flavonoid scaffold. In drug discovery, flavonoids have for long been considered privileged structures52 that preferentially interact with kinases,53 DNA,54 or oxidoreductases.55 The ability of quercetin and related polyphenolic flavonoids to bind to these targets was confirmed by many Xray structures. Quercetin is a relatively small and rigid ligand, and its reactive catechol moiety participates in the formation of specific interactions. For example, in the kinase complex, it contributes to a hydrogen-bond network involving water molecules and residue Glu64 (Figure 1a, left).53 In contrast, in the active site of Ni-quercetinase, the catechol forms π-stacking interactions with Phe78 and a backbone hydrogen bond with residue Val77 (Figure 1a, right).55 Clearly, despite their reactivity and biological membrane perturbation potential, quercetin and other polyphenolic flavonoids are capable of engaging in specific target−ligand interactions, consistent with their proposed privileged structure character. It follows that apparent activity of such catecholcontaining compounds must be carefully evaluated on a caseby-case basis.14,43 Importantly, neither artificial reactivity nor specific activity can be generally assumed. Figure 1b shows a contrasting example revealing consequences of an interference mechanism. Unsaturated fivemembered heterocycles such as rhodanines represent another
Figure 1. X-ray structures with assay interference compounds. Bound crystallographic ligands are shown (carbon atoms are colored green) that contain PAINS substructures implicated in assay interference mechanisms. Protein backbone ribbons and carbon atoms of selected residues are colored in silver or orange. Other ligand and protein atoms are shown using standard atom coloring. Residues discussed in 2760
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765
ACS Omega
Perspective
Figure 1. continued the text are labeled. In molecular graphs, PAINS substructures are colored red. Activity values are reported, if available. (a) Quercetin (PAINS code: catechol_A) in complex with death-associated protein kinase 1 (left, PDB ID: 5AUW) and Ni-quercetinase (right, PDB ID: 5FLJ). In the X-ray structures, the catechol moiety is colored dark red. (b) Inhibitory reaction product in complex with metallo-β-lactamase type 2 (PDB ID: 4PVO) resulting from hydrolysis of an unsaturated rhodanine (PAINS code: ene_rhod_A). The figure was adopted from ref 43 and modified. (c) Mannich base-type inhibitor (PAINS code: mannich_A) in the active site of polyketide synthase 13 (PDB ID: 5V3X). The figure was adopted from ref 43 and modified. All structural representations were generated with Molecular Operating Environment.68
class of prominent PAINS but have also been considered as privileged scaffolds in drug discovery,11,56 similar to flavonoids. Rhodanines are prone to ring-opening reactions and have a strong tendency to react with nucleophilic thiol groups through a Michael-type addition at the exocyclic double bond.57,58 As shown in Figure 1b (bottom left), β-lactamase was cocrystallized with an arylidene rhodanine inhibitor but this inhibitor was not found in the X-ray structure. Rather, the enzyme formed a complex with a smaller hydrolyzed thioenolate that was produced by the rhodanine ring-opening reaction (Figure 1b, bottom right).59 This degradation product specifically bound to β-lactamase with well-defined interactions. The thiol and hydroxyl groups of the inhibitor formed an interaction network with the zinc cations in the active site involving residues Asp120, Cys221, and His263. Hence, the Xray structure identified a reaction product of unsaturated rhodanine as the active component and illustrated the outcome of an interference mechanism. The structurally confirmed “pseudo prodrug” character of rhodanines should be taken into consideration when judging potential activities of these compounds and might also be further explored and exploited in compound design. Figure 1c shows another example of a well-known PAINS motif, a Mannich base. In this case, however, the structural context in which the Mannich base was presented restricted its reactivity and led to specific inhibition. The assay interference potential of Mannich bases primarily depends on their ability to form reactive quinone methides.60 Notably, this structural motif also frequently occurs in natural products.60 However, as shown in Figure 1c, the Mannich base inhibitor complexed with polyketide synthase 1361 contained a tertiary amine as part of an aliphatic piperidine ring system that prevented the formation of the reactive tautomer. Moreover, in the X-ray structure, the aliphatic amine formed a charge-assisted hydrogen bond with the side chain of Asn1640, whereas the phenolic core of the inhibitor was involved in aromatic stacking interactions with Phe1670. Hence, in this case, the structural context of the PAINS motif impeded reactivity and the stabilizing chemical modification led to a potent inhibitor of polyketide synthase. 4.2. Crystallographic Multitarget Compounds and Expansion into Analogue Space. The many crystallographically characterized multitarget and multifamily compounds we identified yielded extensive target coverage and displayed a variety of binding patterns. Figure 2 shows a representative example, the matrix metalloproteinase inhibitor galardin62 in a complex with the zinc metalloproteinase domain of anthrax toxin.63 Galardin was found in six X-ray
Figure 2. Multitarget inhibitor and promiscuous analogues. Shown is galardin bound to the zinc metalloproteinase domain of anthrax toxin (PDB ID: 4PKW). The inhibitor was present in six X-ray structures with three different targets from the same family. The ASB scaffold45 representing a galardin analogue series with modifications at two sites (R1 and R2) is shown on the lower left. Analogues were reported to be active against a total of 22 target proteins from 10 different families. “#” means “number of”. The figure was adopted from ref 44 and modified.
structures bound to three different proteases. Thus, following our classification scheme, galardin was a crystallographic multitarget ligand. A computational search for structural analogues of galardin in ChEMBL identified 22 compounds with substitutions at two sites. The ASB scaffold representing the analogue series is also shown in Figure 2 (bottom left). Taken together, for the analogues, activities against a total of 22 targets belonging to 10 different families were reported, hence suggesting additional target hypotheses for galardin and closely related compounds and indicating high potential for polypharmacology. Thus, galardin was among the structurally confirmed promiscuous compounds for which analogue series were identified. These series provided more than 100 ASB scaffolds as templates for polypharmacological ligand design.44 4.3. Multifamily Compounds and Varying Interaction Hotspots. Given the unexpectedly large number of crystallographic multifamily compounds we identified, a particularly interesting question was how these compounds facilitated binding to distantly or unrelated targets.49 Figure 3 shows a representative example. Indomethacin is a drug with known polypharmacology and a multifamily ligand for which a variety of crystallographic complexes with unrelated targets were available including cyclooxygenases,64 peroxisome proliferatoractivated receptor γ (PPARγ),65 phospholipases,66 and serum albumin.67 As shown in Figure 3, indomethacin bound with similar conformations to cyclooxygenase-2 (top left) and PPARγ (top right), with an all-atom root-mean-square deviation of 1.2 Å. However, in these complexes, polar and van der Waals (vdW) interactions were centered on different regions of indomethacin such that different atoms became target-dependent interaction hotspots. In cyclooxygenase-2, πstacking interactions were formed with the central indole core 2761
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765
ACS Omega
Perspective
family. Nonetheless, the quest for specificity has been the guiding principle in compound design and development for many years. Prior to the availability of profiling assays and compound arrays covering many different targets, specificity was typically evaluated by testing active compounds against their primary (intended) target and a few related ones. Hence, it is not surprising that a rather narrow view of in vitro specificity has dominated drug development for long time. Even the first clinical kinase inhibitors such as imatinib, approved for cancer treatment 15 years ago, were originally thought to be kinase-specific, before it was discovered that they were promiscuous and that polypharmacology was responsible for their efficacy. Kinase inhibitors were only the beginning. Over the past two decades, increasing insights into the complex multifactorial nature and cellular network dependence of many diseases were complemented with mounting evidence that drugs were often active against multiple targets. Consequently, polypharmacology became a new paradigm in drug discovery. Evidence of multitarget activity helped rationalize desired therapeutic as well as undesired side effects. Furthermore, several existing drugs were identified to interact with unrelated targets having diverse functions, which provided the conceptual basis for drug repurposing. Given the increasing popularity of the polypharmacology paradigm, it was frequently proposed that drugs would generally act on multiple targets. In contrast, examples of drugs requiring target specificity were emphasized such as agents used for the treatment of chronic diseases. Such assumptions and controversial views prompted us and others to have a careful look at large collections of bioactive compounds and activity data and determine the frequency and magnitude of promiscuity. Large-scale analyses of activity data, including data from high-throughput screening, did not support the assumption that drugs and other bioactive compounds would generally be promiscuous. Confined subsets of highly promiscuous drugs and other bioactive compounds were identified, including clinical kinase inhibitors. However, the majority of kinase inhibitors reported to date in the medicinal chemistry literature are only annotated with a single kinase. This also applies to other active compounds covering the current spectrum of therapeutic targets. In this context, data incompleteness must be taken into consideration, giving rise to principally low data-driven promiscuity estimates. On the other hand, currently available volumes of compounds and activity data are already so large that statistical trends deduced from these data cannot be disregarded, indicating that data sparseness alone can hardly be responsible for overall low promiscuity detected among bioactive compounds. These observations do not call the principles of polypharmacology into question; neither does polypharmacology per se invalidate the compound specificity paradigm. Rather, it is likely that there is a balance between polypharmacological and target-specific drug actions and ensuing in vivo effects. Regardless, increasing attention has also been paid to the generation of compounds with activity against two or more therapeutic targets of interest. Going forward, multitarget compound design will greatly benefit from a deeper and more refined understanding of multitarget activity at the molecular level of detail. The discussion of structure-based promiscuity analysis presented herein should be viewed in light of the above considerations. Moreover, in this context, another critically important aspect must be taken into account. Promiscuity
Figure 3. Crystallographic interaction hotspots of multifamily ligands. Shown is indomethacin, an exemplary multifamily ligand, bound to cyclooxygenase-2 (left, orange, PDB ID: 4COX) and peroxisome proliferator-activated receptor γ (right, blue, PDB ID: 3ADS). Dashed lines represent polar ligand−target interactions. In each complex, polar and van der Waals (vdW) interactions were determined for each ligand atom using Molecular Operating Environment68 and mapped onto the molecular graph of the ligand.69 Atom positions were colorcoded according to vdW or polar interaction counts using a continuous color spectrum (light, minimal (min) number of respective interactions; dark, maximal (max) number).69 Polar and vdW interactions were separately counted for each participating protein and ligand atom (e.g., if a ligand atom was in vdW contact with three protein atoms, the interaction count was 3). The conformation of indomethacin was similar in both binding sites. However, the two-dimensional ligand depictions reveal the formation of different interaction hotspots in these binding sites.
of the ligand and hydrogen bonds with the polarized methoxy group, whereas the 4-chlorophenyl substituent and an aliphatic methyl residue were involved in vdW contacts. In contrast, in PPARγ, the 4-chlorophenyl moiety participated in stacking interactions and both the methoxy group and the indole core were involved in vdW contacts. In the presence of similar binding conformations, the formation of different targetdependent interaction hotspots provided a molecular rationale for multifamily binding of indomethacin. Equivalent observations were made for many multifamily ligands that interacted with distinct proteins.49
5. PERSPECTIVE Generating active compounds that are target-specific has long been a paradigm in pharmaceutical research. It is well known that specificity is often difficult to achieve, especially among closely related targets such as members of the same enzyme 2762
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765
ACS Omega
Perspective
proteins, confirmed by manual curation of structural data. These compounds provide an extensive knowledge base for rationalizing promiscuity at the molecular level of detail. A key question is how small molecules interact with diverse targets. Are there common patterns? Are the interactions specific that these compounds form in different binding sites? Previous studies have largely focused on exploring binding site similarity to rationalize promiscuous binding events, mostly on a case-bycase basis. Binding site similarity encompasses geometry and shape as well as chemical features, and it is intuitive that similar binding sites should be likely to accommodate the same or similar ligands. This usually applies to members of the same protein family. Our extensive knowledge base of multifamily ligands has enabled us to go further and systematically explore binding conformations of ligands in different binding sites, including unrelated targets, and quantitatively analyze their interaction patterns. Many multifamily ligands interacting with distinct proteins were found to display surprisingly similar binding conformations in different active sites. Many of these active sites had only little or no detectable similarity. Similar binding conformations were not only observed for small and fairly rigid ligands but also for increasingly flexible compounds. In other cases, such ligands bound in different conformations and formed different types of interactions. However, similar binding conformations in different active sites generally resulted in the formation of distinct interaction hotspots. In our studies, the formation of different interaction hotspots on the basis of preferred binding conformations has emerged as a promiscuity-conferring mechanism, which helps us to better understand how small molecules are capable of interacting with unrelated targets. In conclusion, in this work, we have discussed the potential of structure-based promiscuity analysis by focusing on crystallographic ligands in complexes with different targets. Our analyses have relied on the systematic identification of such ligands, which yielded a large, and previously unconsidered, knowledge base of crystallographically confirmed multitarget and multifamily compounds. The structures revealed a variety of binding characteristics and helped to rationalize promiscuity at the molecular level of detail. Promiscuous X-ray ligands we have identified and associated activity information have been made freely available as a part of the original studies.43,44,49
analysis is often influenced by artifacts, which further complicate matters. Importantly, just around the time when the first promiscuous kinase inhibitors were in final stages of clinical development, evidence was mounting that artificial activity readouts were affecting biological assays on a large scale due to various interference mechanisms. Compound aggregation as a cause of activity artifacts first emerged as a common threat to biological screening and medicinal chemistry. A few years later, different facets of chemical reactivity were established as recurrent assay interference mechanisms. These caveats also complicate the study of compound promiscuity and provide another compelling reason to investigate promiscuity at the level of X-ray structures confirming binding events. As we have shown, currently available X-ray data yield many crystallographic multitarget and multifamily ligands, more so than we originally anticipated, which establish a sound basis for structure-based promiscuity exploration. Structurally characterized assay interference compounds such as PAINS provide prototypes for rationalizing promiscuity at different levels. Representative examples presented herein delineate a wide spectrum of possible effects. Even prominent PAINS with undeniable assay interference potential such as catechols or rhodanines have the structurally confirmed ability to interact with a variety of targets. PAINS-containing ligands can be truly promiscuous, although they might also cause assay artifacts. Accordingly, the privileged structure character of, for example, flavonoids or rhodanines, which has long been considered in medicinal chemistry, cannot simply be disregarded, despite their interference potential. However, we have also shown that degradation products of interference compounds might be active and form specific ligand−target interactions. Moreover, chemical modifications of PAINS can prevent undesired reactions to occur and contribute to targetspecific binding. These observations emphasize the critically important role of the structural context in which PAINS or other interference motifs occur. As discussed, there is firm crystallographic evidence for different binding characteristics. Taken together, these findings strongly suggest that activity of interference compounds must be carefully evaluated and judged on a case-by-case basis. Disregarding the ability of such compounds to specifically interact with diverse targets would be as careless scientifically as disregarding their interference potential. We have also shown that taking (noncrystallographic) analogues of multitarget and multifamily compounds into account provides additional information. Compounds forming an analogue series are likely to have similar binding characteristics and share targets. Thus, by considering analogues of crystallographic ligands, additional targets might be inferred, as long as the underlying activity data are carefully curated and reliable. The larger the analogue series containing crystallographic ligands, the more the opportunities to evaluate the consistency of activity annotations and identify “outliers”. Hence, complementing X-ray structures of multitarget or multifamily ligands with analogue series and medicinal chemistry data is a promising approach, as long as the accumulated data are not overinterpreted. We have identified more than 700 crystallographic multifamily ligands in complexes with distantly related or unrelated targets, and for nearly 200 of these ligands, well-defined activity measurements were available. A subset of nearly 100 of these multifamily ligands exclusively interacted with unrelated
■
AUTHOR INFORMATION
Corresponding Author
*E-mail:
[email protected]. Phone: +49-228-7369100. ORCID
Jürgen Bajorath: 0000-0002-0557-5714 Author Contributions
The study was carried out and the manuscript was written with contributions of all authors. All authors have approved the final version of the manuscript. Notes
The authors declare no competing financial interest.
■
ACKNOWLEDGMENTS We thank OpenEye and Chemical Computing Group for academic software licenses. 2763
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765
ACS Omega
■
Perspective
(24) Mestres, J.; Gregori-Puigjané, E.; Valverde, S.; Solé, R. V. The Topology of Drug-Target Interaction Networks: Implicit Dependence on Drug Properties and Target Families. Mol. Biosyst. 2009, 5, 1051− 1057. (25) Jasial, S.; Hu, Y.; Bajorath, J. Determining the Degree of Promiscuity of Extensively Assayed Compounds. PLoS One 2016, 11, No. e0153873. (26) Hu, Y.; Bajorath, J. Entering the ‘Big Data’ Era in Medicinal Chemistry: Molecular Promiscuity Analysis Revisited. Future Sci. OA 2017, 3, No. FSO179. (27) Stumpfe, D.; Tinivella, A.; Rastelli, G.; Bajorath, J. Promiscuity of Inhibitors of Human Protein Kinases at Varying Data Confidence Levels and Test Frequencies. RSC Adv. 2017, 7, 41265−41271. (28) Chong, C. R.; Sullivan, D. J. New Uses for Old Drugs. Nature 2007, 448, 645−646. (29) Keiser, M. J.; Setola, V.; Irwin, J. J.; Laggner, C.; Abbas, A. I.; Hufeisen, S. J.; Jensen, N. H.; Kuijer, M. B.; Matos, R. C.; Tran, T. B.; Whaley, R.; Glennon, R. A.; Hert, J.; Thomas, K. L. H.; Edwards, D. D.; Shoichet, B. K.; Roth, B. L. Predicting New Molecular Targets for Known Drugs. Nature 2009, 462, 175−181. (30) Canny, S. A.; Cruz, Y.; Southern, M. R.; Griffin, P. R. PubChem Promiscuity: A Web Resource for Gathering Compound Promiscuity Data from PubChem. Bioinformatics 2012, 28, 140−141. (31) Keiser, M. J.; Roth, B. L.; Armbruster, B. N.; Ernsberger, P.; Irwin, J. J.; Shoichet, B. K. Relating Protein Pharmacology by Ligand Chemistry. Nat. Biotechnol. 2007, 25, 197−206. (32) Yang, J. J.; Ursu, O.; Lipinski, C. A.; Sklar, L. A.; Oprea, T. I.; Bologa, C. G. Badapple: Promiscuity Patterns from Noisy Evidence. J. Cheminf. 2016, 8, No. e29. (33) Jasial, S.; Gilberg, E.; Blaschke, T.; Bajorath, J. Machine Learning Distinguishes with High Accuracy Between Pan-Assay Interference Compounds that are Promiscuous or Represent Dark Chemical Matter. J. Med. Chem. 2018, 61, 10255−10264. (34) De Franchi, E.; Schalon, C.; Messa, M.; Onofri, F.; Benfenati, F.; Rognan, D. Binding of Protein Kinase Inhibitors to Synapsin I Inferred from Pairwise Binding Site Similarity measurements. PLoS One 2010, 5, No. e12214. (35) Haupt, J. V.; Daminelli, S.; Schroeder, M. Drug Promiscuity in PDB: Protein Binding Site Similarity Is Key. Plos One 2013, 8, No. e65894. (36) Moya-García, A.; Adeyelu, T.; Kruger, F. A.; Dawson, N. L.; Lees, J. G.; Overington, J. P.; Orengo, C.; Ranea, J. A. G. Structural and Functional View of Polypharmacology. Sci. Rep. 2017, 7, No. e10102. (37) Berman, H. M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T. N.; Weissig, H.; Shindyalov, I. N.; Bourne, P. E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235−242. (38) Feng, Z.; Chen, L.; Maddula, H.; Akcan, O.; Oughtred, R.; Berman, H. M.; Westbrook, J. Ligand Depot: A Data Warehouse for Ligands Bound to Macromolecules. Bioinformatics 2004, 20, 2153− 2155. (39) Wang, R.; Fang, X.; Lu, Y.; Wang, S. The PDBbind Database: Collection of Binding Affinities for Protein-Ligand Complexes with Known Three-Dimensional Structures. J. Med. Chem. 2004, 47, 2977−2980. (40) Gaulton, A.; Bellis, L. J.; Bento, A. P.; Chambers, J.; Davies, M.; Hersey, A.; Light, Y.; McGlinchey, S.; Michalovich, D.; Al-Lazikani, B.; Overington, J. P. ChEMBL: A Large-Scale Bioactivity Database for Drug Discovery. Nucleic Acids Res. 2012, 40, D1100−D1107. (41) The UniProt Consortium UniProt. The Universal Protein Knowledgebase. Nucleic Acids Res. 2017, 45, D158−D169. (42) Stumpfe, D.; Dimova, D.; Bajorath, J. Computational Method for the Systematic Identification of Analog Series and Key Compounds Representing Series and Their Biological Activity Profiles. J. Med. Chem. 2016, 59, 7667−7676. (43) Gilberg, E.; Gütschow, M.; Bajorath, J. X-Ray Structures of Target-Ligand Complexes Containing Compounds with Assay Interference Potential. J. Med. Chem. 2018, 61, 1276−1284.
REFERENCES
(1) Hu, Y.; Bajorath, J. Compound Promiscuity - What Can We Learn From Current Data. Drug Discovery Today 2013, 18, 644−650. (2) Gilberg, E.; Jasial, S.; Stumpfe, D.; Dimova, D.; Bajorath, J. Highly Promiscuous Small Molecules from Biological Screening Assays Include Many Pan-Assay Interference Compounds but Also Candidates for Polypharmacology. J. Med. Chem. 2016, 59, 10285− 10290. (3) McGovern, S. L.; Caselli, E.; Griorieff, N.; Shoichet, B. K. A Common Mechanism Underlying Promiscuous Inhibitors from Virtual and High-Throughput Screening. J. Med. Chem. 2002, 45, 1712−1722. (4) Shoichet, B. K. Screening in a Spirit Haunted World. Drug Discovery Today 2006, 11, 607−615. (5) Irwin, J. J.; Duan, D.; Torosyan, H.; Doak, A. K.; Ziebart, K. T.; Sterling, T.; Tumanian, G.; Shoichet, B. K. An Aggregation Advisor for Ligand Discovery. J. Med. Chem. 2015, 58, 1712−1722. (6) Baell, J. B.; Holloway, G. A. New Substructure Filters for Removal of Pan Assay Interference Compounds (PAINS) from Screening Libraries and for Their Exclusion in Bioassays. J. Med. Chem. 2010, 53, 2719−2740. (7) Baell, J. B.; Walters, M. A. Chemistry: Chemical Con Artists Foil Drug Discovery. Nature 2014, 513, 481−483. (8) Baell, J. B.; Nissink, J. W. M. Seven Year Itch: Pan-Assay Interference Compounds (PAINS) in 2017 − Utility and Limitations. ACS Chem. Biol. 2018, 13, 36−44. (9) Capuzzi, S. J.; Muratov, E. N.; Tropsha, A. Phantom PAINS: Problems with the Utility of Alerts for Pan-Assay INterference CompoundS. J. Chem. Inf. Model. 2017, 57, 417−427. (10) Jasial, S.; Hu, Y.; Bajorath, J. How Frequently Are Pan Assay Interference Compounds Active? Large-Scale Analysis of Screening Data Reveals Diverse Activity Profiles, Low Global Hit Frequency, and Many Consistently Inactive Compounds. J. Med. Chem. 2017, 60, 3879−3886. (11) Mendgen, T.; Steuer, C.; Klein, C. D. Privileged Scaffolds or Promiscuous Binders: A Comparative Study on Rhodanines and Related Heterocycles in Medicinal Chemistry. J. Med. Chem. 2012, 55, 743−753. (12) Nelson, K. M.; Dahlin, J. L.; Bisson, J.; Graham, J.; Pauli, G. F.; Walters, M. A. The Essential Medicinal Chemistry of Curcumin. J. Med. Chem. 2017, 60, 1620−1637. (13) Aldrich, C.; Bertozzi, C.; Georg, G. I.; Kiessling, L.; Lindsley, C.; Liotta, D.; Merz, K. M., Jr.; Schepartz, A.; Wang, S. The Ecstasy and Agony of Assay Interference Compounds. J. Med. Chem. 2017, 60, 2165−2168. (14) Hopkins, A. L. Network Pharmacology: The Next Paradigm in Drug Discovery. Nat. Chem. Biol. 2008, 4, 682−690. (15) Peters, J. U. Polypharmacology − Foe or Friend? J. Med. Chem. 2013, 56, 8955−8971. (16) Anighoro, A.; Bajorath, J.; Rastelli, G. Polypharmacology: Challenges and Opportunities in Drug Discovery. J. Med. Chem. 2014, 57, 7874−7887. (17) Rosini, M. Polypharmacology: The Rise of Multitarget Drugs over Combination Therapies. Future Med. Chem. 2014, 6, 485−487. (18) Bolognesi, M. L.; Cavalli, A. Multitarget Drug Discovery and Polypharmacology. ChemMedChem. 2016, 11, 1190−1192. (19) Knight, Z. A.; Lin, H.; Shokat, K. M. Targeting the Cancer Kinome through Polypharmacology. Nat. Rev. Cancer 2010, 10, 130− 137. (20) Gross, S.; Rahal, R.; Stransky, N.; Lengauer, C.; Hoeflich, K. P. Targeting Cancer with Kinase Inhibitors. J. Clin. Invest. 2015, 125, 1780−1789. (21) Nurse, P. The Ends of Understanding. Nature 1997, 387, 657. (22) Roukos, D. H. Networks Medicine: From Reductionism to Evidence of Complex Dynamic Biomolecular Interactions. Pharmacogenomics 2011, 12, 695−698. (23) Mestres, J.; Gregori-Puigjane, E.; Valverde, S.; Sole, R. V. Data Completeness − The Achilles Heel of Drug-Target Networks. Nat. Biotechnol. 2008, 26, 983−984. 2764
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765
ACS Omega
Perspective
(44) Gilberg, E.; Stumpfe, D.; Bajorath, J. X-Ray Structure Based Identification of Compounds with Activity against Targets from Different Families and Generation of Templates for Multi-Target Ligand Design. ACS Omega 2018, 3, 106−111. (45) Dimova, D.; Stumpfe, D.; Bajorath, J. Computational Design of New Molecular Scaffolds for Medicinal Chemistry, Part II: Generalization of ASB Scaffolds. Future Sci. OA 2018, 4, No. FSO267. (46) Morphy, R.; Rankovic, Z. Designed Multiple Ligands. An Emerging Drug Discovery Paradigm. J. Med. Chem. 2005, 48, 6523− 6543. (47) Hopkins, A. L.; Mason, J. S.; Overington, J. P. Can We Rationally Design Promiscuous Drugs? Curr. Opin. Struct. Biol. 2006, 16, 127−136. (48) Proschak, E.; Stark, H.; Merk, D. Polypharmacology by Design: A Medicinal Chemist’s Perspective on Multitargeting Compounds. J. Med. Chem. 2019, 62, 420−444. (49) Gilberg, E.; Gütschow, M.; Bajorath, J. Promiscuous Ligands from Experimentally Determined Structures, Binding Conformations, and Protein Family Dependent Interaction Hotspots. ACS Omega 2019, 4, 1729−1737. (50) Dahlin, J. L.; Nissink, J. W. M.; Strasser, J. M.; Francis, S.; Higgins, L.; Zhou, H.; Zhang, Z.; Walters, M. A. PAINS in the Assay: Chemical Mechanisms of Assay Interference and Promiscuous Enzymatic Inhibition Observed During a Sulfhydryl-Scavenging HTS. J. Med. Chem. 2015, 58, 2091−2113. (51) Ingólfsson, H. I.; Thakur, P.; Herold, K. F.; Hobart, E. A.; Ramsey, N. B.; Periole, X.; de Jong, D. H.; Zwama, M.; Yilmaz, D.; Hall, K.; Maretzky, T.; Hemmings, H. C., Jr.; Blobel, C.; Marrink, S. J.; Kocer, A.; Sack, J. T.; Andersen, O. S. Phytochemicals Perturb Membranes and Promiscuously alter Protein Function. ACS Chem. Biol. 2014, 9, 1788−1798. (52) Reis, J.; Gaspar, A.; Milhazes, N.; Borges, F. Chromone as a Privileged Scaffold in Drug Discovery: Recent Advances. J. Med. Chem. 2017, 60, 7941−7957. (53) Yokoyama, T.; Kosaka, Y.; Mizuguchi, M. Structural Insight into the Interactions between Death-Associated Protein Kinase 1 and Natural Flavonoids. J. Med. Chem. 2015, 58, 7400−7408. (54) Srivastava, S.; Somasagara, R. R.; Hegde, M.; Nishana, M.; Tadi, S. K.; Srivastava, M.; Choudhary, B.; Raghavan, S. C. Quercetin, a Natural Flavonoid Interacts with DNA, Arrests Cell Cycle and Causes Tumor Regression by Activating Mitochondrial Pathway of Apoptosis. Sci. Rep. 2016, 6, No. e24049. (55) Jeoung, J. H.; Nianios, D.; Fetzner, S.; Dobbek, H. Quercetin 2,4-Dioxygenase Activates Dioxygen in a Side-on O2-Ni Complex. Angew. Chem., Int. Ed. 2016, 55, 3281−3284. (56) Tomasíc, T.; Masic, L. P. Rhodanine as a Privileged Scaffold in Drug Discovery. Curr. Med. Chem. 2009, 16, 1596−1629. (57) Metz, J. T.; Huth, J. R.; Hajduk, P. J. Enhancement of Chemical Rules for Predicting Compound Reactivity Towards Protein Thiol Groups. J. Comput.-Aided Mol. Des. 2007, 21, 139−144. (58) Carlson, E. E.; May, J. F.; Kiessling, L. L. Chemical Probes of UDP-Galactopyranose Mutase. Chem. Biol. 2006, 13, 825−837. (59) Brem, J.; van Berkel, S. S.; Aik, W.; Rydzik, A. M.; Avison, M. B.; Pettinati, I.; Umland, K.; Kawamura, A.; Spencer, J.; Claridge, T. D. W.; McDonough, M. A.; Schofield, C. J. Rhodanine Hydrolysis Leads to Potent Thioenolate Mediated Metallo-β-Lactamase Inhibition. Nat. Chem. 2014, 6, 1084−1090. (60) Baell, J. B. Feeling Naturés PAINS: Natural Products, Natural Product Drugs, and Pan Assay Interference Compounds (PAINS). J. Nat. Prod. 2016, 79, 616−628. (61) Aggarwal, A.; Parai, M. K.; Shetty, N.; Wallis, D.; Woolhiser, L.; Hastings, C.; Dutta, N. K.; Galaviz, S.; Dhakal, R. C.; Shrestha, R.; Wakabayashi, S.; Walpole, C.; Matthews, D.; Floyd, D.; Scullion, P.; Riley, J.; Epemolu, O.; Norval, S.; Snavely, T.; Robertson, G. T.; Rubin, E. J.; Ioerger, T. R.; Sirgel, F. A.; van der Merwe, R.; van Helden, P. D.; Keller, P.; Bottger, E. C.; Karakousis, P. C.; Lenaerts, A. J.; Sacchettini, J. C. Development of a Novel Lead that Targets M. tuberculosis Polyketide Synthase 13. Cell 2017, 170, 249−259.
(62) Santiskulvong, C.; Rozengurt, E. Galardin (GM 6001), a BroadSpectrum Matrix Metalloproteinase Inhibitor, blocks Bombesin- and LPA-induced EGF Receptor Transactivation and DNA Synthesis in Rat-1 Cells. Exp. Cell Res. 2003, 290, 437−46. (63) Maize, K. M.; Kurbanov, E. K.; De La Mora-Rey, T.; Geders, T. W.; Hwang, D. J.; Walters, M. A.; Johnson, R. L.; Amin, E. A.; Finzel, B. C. Anthrax Toxin Lethal Factor Domain 3 is Highly Mobile and Responsive to Ligand Binding. Acta Crystallogr., Sect. D: Struct. Biol. 2014, 70, 2813−2822. (64) Kurumbail, R. G.; Stevens, A. M.; Gierse, J. K.; McDonald, J. J.; Stegeman, R. A.; Pak, J. Y.; Gildehaus, D.; Miyashiro, J. M.; Penning, T. D.; Seibert, K.; Isakson, P. C.; Stallings, W. C. Structural Basis for Selective Inhibition of Cyclooxygenase-2 by Anti-Inflammatory Agents. Nature 1996, 384, 644−648. (65) Waku, T.; Shiraki, T.; Oyama, T.; Maebara, K.; Nakamori, R.; Morikawa, K. The Nuclear Receptor PPARγ Individually Responds to Serotonin- and Fatty Acid-Metabolites. EMBO J. 2010, 29, 3395− 3407. (66) Singh, N.; Kumar, R. P.; Kumar, S.; Sharma, S.; Mir, R.; Kaur, P.; Srinivasan, A.; Singh, T. P. Simultaneous Inhibition of AntiCoagulation and Inflammation: Crystal Structure of Phospholipase A2 Complexed with Indomethacin at 1.4 Å Resolution Reveals the Presence of the New Common Ligand-Binding Site. J. Mol. Recognit. 2009, 22, 437−445. (67) Bogdan, M.; Pirnau, A.; Floare, C.; Bugeac, C. Binding Interaction of Indomethacin with Human Serum Albumin. J. Pharm. Biomed. Anal. 2008, 47, 981−984. (68) Chemical Computing Group. Molecular Operating Environment (MOE), 2018.01; Chemical Computing Group ULC: 1010 Sherbooke St. West, Suite #910, Montreal, QC, Canada, H3A 2R7, 2018. (69) Clark, A. M.; Labute, P. 2D Depiction of Protein-Ligand Complexes. J. Chem. Inf. Model. 2007, 47, 1933−1944.
2765
DOI: 10.1021/acsomega.8b03639 ACS Omega 2019, 4, 2758−2765