- Open Access
Genome-wide characterisation of the binding repertoire of small molecule drugs
Human Genomicsvolume 1, Article number: 41 (2003)
Most, if not all, drugs interact with multiple proteins. One or more of these interactions are responsible for carrying out the primary therapeutic effects of the drug. Others are involved in the transport or metabolic processing of the drug or in the mediation of side effects. Still others may be responsible for activities that correspond to alternate therapeutic applications. The potential clinical impact of a drug and its cost of development are affected by the sum of all these interactions. The drug development process includes the identification and characterisation of a drug's clinically relevant interactions. This characterisation is presently accomplished by a combination of experimental laboratory techniques and clinical trials, with increasing numbers of patient participants. Efficient methods for the identification of all the molecular targets of a drug prior to clinical trials could greatly expedite the drug development process. Combinatorial peptide and cDNA phage display have the potential for achieving a complete characterisation of the binding repertoire of a small molecule. This paper will discuss the current state of phage display technology, as applied to the identification of novel receptors for small molecules, using a successful application with the drug Taxol™ as an example of the technical and theoretical benefits and pitfalls of this method.
The concept of a receptor as the component of an organism with which a chemical agent interacts was based originally on the independent work of Paul Ehrlich and J. N. Langley during the late nineteenth century. Pharmaceutical agents are generally believed to exert their therapeutic effect by binding to a particular protein or nucleic acid target . Both natural product and synthetic chemical drugs, however, can and do have multiple macromolecule-binding partners (Table 1) [2–48], such as metabolic breakdown enzymes (ie cytochrome P450 enzymes), side effect-mediating targets and, in some cases, secondary therapeutic targets (see aspirin and Taxol™ below as specific examples). Recognition of the binding promiscuity of drug molecules is reflected in the recent assembly of a publicly accessible web-based therapeutic target database .
For a drug to reach the market, an average input of 14 years and US$359 million is required. The assessment of toxicity, side effects and metabolic breakdown pathways consumes a large percentage of this investment[51–53]. The speed with which a drug passes through the various stages of the drug discovery process can be a critical factor for success. Novel strategies that provide information about the potential target proteins of a small molecule are badly needed to augment the current drug development process. A complete accounting of the repertoire of targets for a particular drug provides important information towards a complete understanding of the mode of action of a drug, as well as the prediction of potential complications from its use in the clinic. Detection of potential toxicity and side effect properties early on during drug development can translate into large financial savings and lowered patient risk, as efforts can be redirected towards drug leads that are more likely to be safe for patient use. Unfortunately, identifying macromolecules that bind to small molecule drugs can be a slow and laborious process, particularly if the target(s) are expressed at low levels in their respective tissue(s). Although computer methods for toxicity prediction offer a fast, attractive solution to this problem [54–56], even the best of this predictive software generates hypotheses based upon crystal structure-based rigid protein docking studies. In many cases, small molecules (ie non-peptide molecules of less than about 1,000 Daltons) appear to interact with their receptors through regions that are disordered prior to the interaction . This suggests that computer-based docking strategies may be limited in their success with a subset of potential drug/protein interactions.
In this paper, we outline a strategy for identifying the targets of small molecules based on the screening of protein or combinatorial peptide libraries displayed on phage particles. The ultimate goal is the complete delineation of all the targets for a small molecule. Since known receptors are mainly proteins, the discussion will be limited to the identification of protein targets. These strategies have been validated using libraries of cDNA expression products[59–61] and libraries of combinatorial peptides. The latter have been used to identify a novel target for Taxol™ , and carry out an in vitro approach to the identification of ATP-binding proteins within entire proteomes (Makowski et al., in preparation).
The number of targets of a drug
Most, if not all, drugs interact with multiple macromolecules. Some of these molecules conform to the concept of a 'receptor', in that they are responsible for mediating the primary thera peutic effects of the drug. Other target molecules would not classically be referred to as receptors because they are involved in transport, chemical processing or the side effects of the drug. Some of these interactions are non-specific, such as the inter actions between hydrophobic drugs and serum proteins. These interactions are, however, critical elements in defining the clinical effect of the drug. Table 1 lists the known 'specific' macromolecular targets for a few representative drugs. Some of the best drugs may be more effective than related molecules because they bind two or more targets relevant to their clinical action. Aspirin, the most potent anti-inflammatory agent identified to date, has at least three targets that mediate anti-inflammatory activity. The anti-cancer activity of Taxol™ is greater than that of other mitotic spindle poisons, presumably because it interacts with the anti-apoptotic protein Bcl-2, as well as with β-tubulin.
A small molecule that binds to no macromolecules (if such a thing exists) will be biologically inert. Conversely, a small molecule that binds to hundreds of proteins may possess numerous therapeutic applications but will have a very high chance of being toxic, as any one of these interactions might disrupt an essential molecular function. A good drug, there fore, should be expected to bind to somewhere between one and a dozen proteins. One, or possibly a few, of these inter-actions will mediate the primary activity of the drug. The necessity of specific action demands that small molecule drugs bind with high affinity and specificity to a relatively limited number of protein targets.
Between 3,000 to 10,000 new disease-related targets are expected to emerge out of the Human Genome Project and related sequencing projects such as the Microbial Genome Program of the US Department of Energy . Parallel synthesis or combinatorial chemistry is similarly expected to generate a coordinate build-up of lead compounds. This concurrent accumulation of targets and drug candidates may not automatically translate into an increased introduction of Food and Drug Administration-approved pharmaceuticals into the medical marketplace. Although high-throughput combinatorial chemistry can create staggering molecular diversity, this diversity may be irrelevant or only semi relevant to the drug target population.
Natural products and their derivatives have evolved to bind to a relatively limited number of proteins in order to be non-toxic to the organism from which they derive. The only exceptions to this rule are a few widely used metabolites, such as ATP, that bind to many different proteins. Synthetic molecules are not so limited. No evolutionary pressure has restricted the number of their biological targets on the basis of inactivity or toxicity. It is possible that the members of many classes of synthetic molecules will bind to many different macromolecules, thus increasing the percentage that will fail toxicity screening. This effect may be partially responsible for the difficulties that combinatorial chemistry has had in superseding natural product leads in the drug development pipeline. There has been a growing trend towards the synthesis of focused, natural product-like libraries which 'partake of the "wisdom" of natural products' . This focused library approach may limit the chemical diversity to a biologically relevant subpopulation.
Methods for identification of the molecular targets of a small molecule
Conceptually, it should be possible to develop methods that are capable of identifying all the targets of a small molecule drug or drug candidate. One could conceive of an ideal tool -- a biological interaction chip containing an indexed set of all the molecular surfaces available to a drug molecule within the human body (and within human pathogens and animal models for human pathologies). A simple screen, allowing for very high sensitivity detection of a drug candidate bound to each of the proteins on this chip, could then be used to identify the dozen or so targets expected for a successful drug candidate and rule out promiscuous, potentially toxic drug candidates. An ideal, universal system of this kind is far beyond the relatively simple protein or peptide chips currently under development . What is available now?
The most obvious path to a full accounting of the targets of a small molecule is to create an affinity column containing an immobilised form of the small molecule (preferably attached via tethers to at least two different locations so that every portion of the small molecule is free for binding in at least one of the tagged forms). A cell extract poured through the column will allow those molecules with high affinity to the drug to be bound to the column. Identification of those targets, either by gel electrophoresis and/or mass spectrometry, provides an accounting of the molecules with high affinity for the immobilised small molecule. This strategy has been successful for the identification of numerous molecular targets and will continue to be used. It may be difficult to identify poorly expressed targets, targets that are membrane bound, or large complexes via this technique, as well as to distinguish true positives from false positives.
Bacteriophage display libraries
An efficient way of improving the signal-to-noise ratio in an affinity selection process is to physically link a phenotype to a biologically 'amplifiable' genotype. Several such 'surface display systems' are now in use (such as ribosome display and cell surface display)[67–69] and some of these have been validated for novel drug-protein interaction identification (see McPherson et al., where a tagged version of the drug FK506 was used as a probe molecule against libraries of mRNA-protein fusion molecules constructed from human liver, kidney and bone marrow transcripts). The simplest and most widespread cost-effective systems for carrying out this linkage, however, are bacteriophage display libraries . In these libraries, the DNA encoding a protein or peptide is inserted into the gene of a structural protein of the phage in such a way that its expression product is displayed on the surface of the phage as a fusion to the structural protein, but without disrupting an essential function of that structural protein. A library of phages, each displaying a different expression product, is screened for members that exhibit affinity to a target (using an affinity column, as described above, or a related technique referred to as 'biopanning'). Even if this affinity selection results in only a handful of positive 'hits', because of the linkage of phenotype to genotype, it is possible to characterise these 'hits' because the affinity-selected phage can be grown and the sequence of the inserted DNA ascertained. This process circumvents the need to carry out highly sensitive chemical analysis of small amounts of material separated via binding to the drug target. The majority of bacteriophage display systems have been based upon the filamentous phage M13 , but the inherent limitations of a biological system based upon membrane extrusion for replication has led to the development of phage display systems based upon suchlytic phages as T7 and λ .
cDNA expression product libraries
It is now possible to construct libraries of cDNA expression products displayed on the surface of phage particles and to screen them for members exhibiting affinity to small molecules. cDNA sequences are inserted into the gene for a structural protein in such a way that the expression product is displayed on the surface of a virion. Even if only a single phage particle is selected for affinity to a small molecule, it is possible to identify the gene product responsible because the virus can be grown in quantities sufficient for sequencing the viral DNA which contains an insert coding for the gene product displayed on the surface. A similar approach using genomic libraries is also possible, whereby whole genomes are chemically fragmented and cloned into the phage genome [75–79]. A principal shortcoming of this method is that many gene products do not fold properly on a phage surface, particularly large proteins and membrane proteins. This draw back can be alleviated by the use of genomic DNA fragmented into epitope-sized coding pieces as mentioned above. A second problem with full-length protein display is that most phage display systems append the foreign protein onto the carboxy terminus of the phage structural protein, due to the structural requirements of the viral particle (such as a number of commercially available T7 cDNA display libraries). This results in the translation of many cDNA sequences in the incorrect reading frame, as well as the inclusion of 'nonsense' peptides translated from 3' non-coding regions. A number of systems are available which attach the foreign protein to the amino terminus of the phage structural protein, thus weeding out DNA insert sequences which do not fully translate in-frame to the beginning of the phage structural protein gene (for example, M13 gene III-based libraries) .
Several successful uses of this strategy have been reported. Sche et al., used a cDNA phage display library to isolate clones of the FKBP12 protein from a human brain cDNA library using a biotinylated FK506 probe molecule [59, 60]. During multiple rounds of affinity selection with the immobilised drug molecule, the FKBP12 gene emerged as the dominant library member with detection possible via polymerase chain reaction amplification after the second round of affinity selection. Six rounds of affinity selection with the modified drug were required for the FKBP gene to be the predominant clone within the selected population. A similar phage-based approach has been used to identify the protein hNopp140 as a binding partner for the drug doxorubicin . DNA-binding proteins were successfully affinity-selected from a yeast genomic library displayed on a λ phage display vector . Carbohydrate-binding proteins have also been identified by similar techniques .
These representative results demonstrate that, in cases where the targets can be successfully displayed in their proper three-dimensional configuration on the surface of phage particles, the use of cDNA expression product display libraries can provide a powerful approach to the identification of the molecular targets of small molecule drugs. Many targets, such as membrane proteins, however, may require other strategies, such as the fragmented epitope approaches mentioned above or the use of combinatorial peptide libraries outlined below.
Combinatorial peptide libraries
An alternative to the use of cDNA-expression product display is combinatorial peptide display. In this method, random synthetic oligonucleotides are inserted into the gene for a structural protein at a site known to be presented on the surface of the phage particle; the resulting population is screened for members exhibiting affinity for the small molecule probe using multiple rounds of biopanning. The sequences of the affinity-selected peptides are then compared with the sequences of naturally occurring proteins and those proteins with segments having the greatest sequence similarity to the affinity-selected peptides are identified as potential targets for the small molecule. This process has been successfully used to identify novel targets for small molecules, as well as to map the binding sites of small molecules on proteins [3, 58].
The principal limitation of this technique is that, for a successful outcome, the binding properties of peptides on the surface of the phage must mimic the binding properties of similar sequences in naturally occurring proteins. This will not always be the case. Our understanding of how a protein recognises a small molecule has changed little since 1894, when Emil Fischer first proposed that '. . .enzyme and glucoside have to fit to each other like a lock and key. . .'. This metaphor of protein structure conveys an image of rigidity that belies the complexity of the process of molecular recognition. It is becoming increasingly clear that many binding sites exhibit significant disorder prior to interaction with a ligand, particularly small molecule-binding sites, with the protein-ligand complex moving through a series of progressively lower energy states towards the bottom of a 'binding funnel' (similar to the proposed 'folding funnel' model for protein folding) [82, 83]. The initial steps in the process of molecular recognition of a small molecule by a protein might well involve the interaction of a disordered segment of peptide with the small molecule, suggesting that the structural context of the interaction is far less important than implied by the 'lock and key' vision . In this case, peptides on the surface of a phage might well reflect the binding properties of peptide segments of naturally occurring proteins, at least during the initial 'handshake' stage of the binding process.
Elucidation of three-dimensional protein structures via X-ray crystallography and nuclear magnetic resonance spectroscopy has provided extensive evidence that small molecule-binding sites are made of short stretches of peptide, usually arising from different discontinuous portions of the polypeptide chain. This would suggest that the sequence of any one portion of the protein would make an energetically trivial contribution to the total binding energy. But a survey of the structures of protein-small molecule interactions has shown that virtually all ligands larger than 300 Daltons interact with at least five contiguous amino acids in at least one peptide segment , implying that there is an energetic 'hot spot' in the ligand-binding site similar to those identified in protein-protein interactions via mutagenesis studies [84–86]. A peptide segment, five amino acids in length, is adequate to define a binding motif identifiable on a genome-wide scan.
It may be possible to categorise small molecule-binding sites by their level of flexibility prior to ligand binding. Those binding sites that act most like the classical 'lock and key' are rigid prior to binding. Although these may be detected through the screening of cDNA display libraries, it is unlikely that they will be identified through affinity screening of peptide libraries. Detection of these receptors may be more amenable to cDNA display approaches as outlined above, or even to in silico docking approaches. Other sites may have rigid bases and flexible over lying peptides, as in examples presented in a survey by Freire . Portions of these binding sites appear to be detectable using peptide phage display (see Taxol™ discussion below). Finally, there appear to be proteins that are significantly disordered prior to binding a ligand, such as the bacterial quorum-sensing pheromone from Agrobacterium tumefaciens . These highly disordered binding sites may well be the most easily mimicked by the binding properties of peptides on the surface of phage particles and appear to be detectable through the use of peptide phage display (again, see Taxol™ discussion below). Initial successes with this approach suggest an important role for random peptide display as a drug target identification tool.
The case of Taxol™
Taxol™ is a natural product with broad cytotoxic activity against a range of human tumours. It was first isolated from the bark of Pacific yew trees, and is now produced in a two-step semi-synthesis from precursors extracted from the needles of farmed yew bushes . The primary activity of Taxol™ is its cytotoxicity towards rapidly dividing cells, apparently brought on through p53-mediated apoptosis initiated by binding to β-tubulin, a process that induces microtubule bundling and halts mitosis.
Isolation of Taxol™-selected peptides
Libraries of random 12-mer and disulfide-bond cyclised random 7-mer peptides displayed at the amino terminus of the mature pIII of M13 phage were selected for members that exhibited affinity to a biotinylated derivative of Taxol™ which had been immobilised on streptavidin-coated plates [3, 58]. The sequences of Taxol™-selected peptides were analysed in three ways: (i) The sequences of proteins known to bind to Taxol™ were compared with the Taxol™-selected sequences as an internal control; (ii) the sequences of Taxol™-selected peptides were analysed to identify motifs that occurred in multiple peptides and these motifs were used to search the human genome database for novel Taxol™-binding proteins; and (iii) the full-length sequences of the Taxol™-selected peptides were compared with the sequences of proteins on a genome-wide scale to identify those proteins with regions exhibiting sequences with high similarity to the Taxol™-selected peptide population.
Only weak consensus sequences were found within either of the two selected peptide populations; however, the physi-cochemical properties of the selected phage populations were statistically evaluated and found to be substantially different from that of an unselected population . This observation suggested that, in spite of the absence of a strong consensus sequence, the selection process had been successful. A comparison of the Taxol™-binding peptide sequences to known Taxol™-binding proteins was carried out to determine if any peptides were mimicking known Taxol™-binding motifs.
Comparison of the Taxol™-selected peptides to the sequence of β-tubulin
Initial comparison of sequence motifs from the Taxol™-selected populations with the sequence of human β-tubulin suggested no detectable similarity . The development of algorithms which made it possible to compare the complete sequences of all the selected peptides with the sequences of naturally occurring proteins allowed for the detection of similarities between the sequences of the Taxol™-binding peptides and the sequence of β-tubulin, however . Figure 1 is a depiction of the three-dimensional structure of β-tubulin, in which the colour table has been chosen to reflect the similarity of the tubulin sequence to that of the random 12-mers selected for affinity to Taxol™. Examination of this Figure indicates that the region of highest similarity is a surface loop in direct contact with the bound Taxol™. This may be an example of the rigid base/flexible overlying loop scenario discussed above. Figure 1 demonstrates that the sequence of at least one segment of β-tubulin involved in Taxol™ binding is a weak motif within the sequences of the Taxol™-selected peptides. Nevertheless, the representation of this putative Taxol™-binding motif is relatively weak in the Taxol™-selected peptide population. Although the surface loop in contact with the Taxol™ molecule exhibits the highest similarity to the sequences of the Taxol™-selected peptides in β-tubulin, this similarity is not exceptionally high in comparison with that found when scanning the entire genome.
A similar comparison of the affinity-selected confor-mationally constrained 7-mer population with that of β-tubulin indicated that the protein segment of highest similarity to the peptides was remote from the Taxol™-binding site, but that the second-highest similarity peak corresponded to a protein segment in direct contact with Taxol™ (results not shown). Again, the observed level of similarity was not exceptional in comparison with whole genome scanning.
Other Taxol™-binding proteins
A comparison of the sequence of the human multi-drug trans porter MDR1 (P-glycoprotein) to the sequences of Taxol™-selected peptides resulted in the prediction of two putative Taxol™-binding regions within the P-glycoprotein sequence . Both of these regions fall within the putative drug-binding regions of the protein that have been identified previously using photo-affinity labelling [4, 5], demonstrating that the data are consistent with those peptides possessing Taxol™-binding motifs.
Using biotin-labelled Taxol™, Byrd et al. identified two Taxol™ targets from both mouse macrophages and brain as heat shock proteins of the 70- and 90-kDa families . The structure of Hsp90 bound to geldamycin, a specific inhibitor of the Hsp90 family, has been elucidated . The Hsp90 sequence has three segments exhibiting similarity to the sequences of the Taxol™-selected peptides. Two of these are in alpha-helical structures that flank the geldamycin binding site; the third is at the end of a beta strand beneath the geldamycin-binding site, as shown in Figure 2. These results, when compared with the previous analyses, predict that Taxol™ will bind Hsp90 at the same site as geldamycin.
Experimental evidence from affinity chromatography suggests that Taxol™ interacts with murine CD18 in a manner mimicking that of lipopolysaccharide to produce an inflammatory response . Although the three-dimensional structure of one subunit of CD18 has been elucidated and the Taxol™-selected peptides have been used to predict regions with high probability of binding Taxol™, there are no biochemical data available to evaluate the Taxol™-binding site on CD18 predicted by this analysis.
Genome-wide survey of Taxol™-binding motifs
An analysis of the Taxol™-selected peptides identified several sequence motifs which appeared multiple times in the peptides. Sequences of all proteins in the human genome database were searched to identify proteins containing these sequences. This process led to the identification of the anti-apoptotic protein Bcl-2 as a candidate Taxol™-binding protein containing the amino acid sequence motif HTPHP. Enzyme-linked immunosorbent assay binding studies and circular dichroism spectroscopy were used to show that the affinity of Taxol™ for Bcl-2 is in the 200 nM range . Calculation of the similarity of the Taxol™-selected peptides to the sequence of Bcl-2 indicated that the most likely site for the binding of Taxol™ to Bcl-2 was within a 50 amino acid unstructured loop near the amino terminus of the protein. A model of the three-dimensional structure of Bcl-2 is shown in Figure 3. This model is based upon the three-dimensional structure of Bcl-XL and the more recent structures of Bcl-2, in which the unstructured loop was excised and replaced by the corresponding sequence from Bcl-XL .
The significance of the Taxol™-Bcl-2 interaction
Is the interaction of Taxol™ with Bcl-2 relevant to the clinical activity of Taxol™? It has long been known that Bcl-2 is phos-phorylated and inactivated in cells treated with Taxol™ [93, 94]. This inactivation blocks the anti-apoptotic activity of Bcl-2, leading to apoptosis, a major mechanism of the cytotoxic activity of Taxol™. It has been noted, however, that Taxol™ is a much more potent anti-cancer drug than other tubulin-binding mitotic spindle poisons. Experiments have pointed to the existence of a non-p53-dependent pathway towards apoptosis subsequent to the addition of Taxol™ to the cell culture [95–97].
The direct interaction of Taxol™ with Bcl-2 may mediate a predicted alternative pathway to apoptosis, independent of mitotic block and the subsequent activity of p53. This hypothesis was reinforced when it was shown that application of Taxol™ to isolated mitochondria induced the permeability transition that results in the release of cytochrome c, an obligate step in the apoptotic pathway. Bcl-2 is anchored in the outer membrane of mitochondria through a C-terminal tail, making it a potential mitochondrial target for mediating this activity. Although tubulin has also been shown to be a component of mitochondrial membranes , making it possible that interaction with tubulin also induces this activity, the functional connection between tubulin and the molecular events involved in the mitochondrial permeability transition has been shown to involve Bcl-2. Removal of the flexible loop of Bcl-2 blocks the apoptotic action of Taxol™ [100, 101], as does phosphorylation . It is a reasonable hypothesis that when this flexible regulatory loop is bound to Taxol™, Bcl-2 is incapable of preventing the progression to apoptosis.
Whole-genome similarity calculations
In the genome-wide scan for potential Taxol™-binding proteins, the authors identified a single protein, Bcl-2, that has since been shown to exhibit significant affinity for Taxol™ . None of the previously known Taxol™-targets were identified in this scan. There are several possible reasons for this, some technical and some fundamental. From a technical perspective, only 69 peptides selected for binding to Taxol™ were used for this novel ligand identification. In their more recent work, aimed at the identification of ATP-binding proteins (Makowski et al., 2003; in preparation), the authors routinely use the sequences of 500 affinity-selected peptides for genome-wide scans. While other Taxol™-binding motifs were present in the Taxol™-selected peptides, such as those matching the β-tubulin and P-glycoprotein sequences, these were less strongly represented than the motif found in Bcl-2. Sequencing of a larger number of peptide sequences makes it possible to observe multiple copies of the motifs that exhibit affinity for the target small molecule, greatly enhancing the likely success of a genome-wide search.
A more fundamental consideration is the fact that Taxol™ binding to Bcl-2 involves a 50 amino acid flexible loop, a structural context readily mimicked by the phage-displayed peptides. The Taxol™-binding sites in the previously identified target structures (eg β-tubulin) appear to be more rigidly structured, involving several discontinuous segments of peptide. Any one peptide segment involved in a relatively rigid binding site may correspond to a motif that exhibits relatively low affinity on its own, leading to a relatively lower abundance in the affinity-selected peptides. These short motifs are observable -- these loops were picked out in a scan of the individual proteins as the most likely site within the protein to be involved in binding of Taxol™ -- but their observation in a genome-wide scan presents a formidable signal-to-noise problem. Careful analysis of the results of the authors' ongoing study of ATP-binding peptides, utilising the extensive body of knowledge available for ATP-binding proteins, is being used to provide a basis for further development of genome-wide methods.
The results from the Taxol™ affinity screening suggest that the utilisation of combinatorial peptide phage display to identify novel drug targets will be most successful for targets with binding sites that are highly flexible prior to ligand binding. Considerable technical improvements in the current methods are possible, both experimental and algorithmic, and the greatest potential is for identification of disordered binding sites that are difficult to find by other methods. cDNA and other display methods may be more appropriate for finding binding sites that involve relatively rigid constructs or multiple punctate interactions with short peptide segments. In silico docking methods for the prediction of binding sites are also more likely to be successful in the identification of rigid binding sites in target proteins. A complete accounting of the binding repertoire of a small molecule drug will require a combination of methods, each best suited to a subset of the potential targets.
The results outlined here indicate that utilisation of phage display of cDNA expression products and random peptides holds significant promise as a strategy for identification of all the protein targets of a small molecule, as well as identification of specific regions in those proteins that may be involved in the binding process.
Drew J: 'Drug discovery: A historical perspective'. Science. 2000, 287: 1960-1964. 10.1126/science.287.5460.1960.
Schiff PB, Fant J, Horwitz SB: 'Promotion of microtubule assembly in vitro by taxol'. Nature. 1979, 277: 665-667. 10.1038/277665a0.
Rodi DJ, Janes RW, Sanganee HJ, et al: 'Screening a library of phage-displayed peptides identifies human Bcl-2 as a Taxol™-binding protein'. J Mol Biol. 1999, 285: 197-204. 10.1006/jmbi.1998.2303.
Wu Q, Bounard PY, Kuduk SD, et al: 'Identification of the domains of photoincorporation of the 3- and 7-benzophenone analogues of Taxol™ in the carboxyl-terminal half of murine mdr1b P-glycoprotein'. Biochemistry. 1998, 37: 11272-11279. 10.1021/bi980598b.
Ambudkar SV, Dey S, Hrycyna CA, et al: 'Biochemical, cellular and pharmacological aspects of the multi-drug transporter'. Ann Rev Pharmacol Toxicol. 1999, 39: 361-398. 10.1146/annurev.pharmtox.39.1.361.
Bhat N, Perera PY, Carboni JM, et al: 'Use of a photoacti-vatable Taxol™ analogue to identify unique cellular targets in murine macrophages: Identification of murine CD18 as a major Taxol™-binding protein and a role for Mac-1 in Taxol™-induced gene expression'. J Immunol. 1999, 162: 7335-7342.
Dhawan V, Swaffar DS: 'A unique paclitaxel-mediated modulation of the catalytic activity of topoisomerase IIα'. Anticancer Drugs. 1999, 10: 397-404. 10.1097/00001813-199904000-00008.
Byrd CA, Bornmann W, Erdjument-Bromage H, et al: 'Heat shock protein 90 mediates macrophage activation by Taxol™ and bacterial lipopolysacchride'. Proc Natl Acad Sci USA. 1999, 96: 5645-5650. 10.1073/pnas.96.10.5645.
Vane JR: 'Inhibition of prostaglandin synthesis as a mechanism of action for aspirin-like drugs'. Nat New Biol. 1971, 231: 232-235.
Yin MJ, Yamamoto Y, Gaynor RB: 'The anti-inflam matory agents aspirin and salicylate inhibit the activity of IκB kinase β'. Nature. 1998, 396: 77-80. 10.1038/23948.
Puscas I, Coltau M, Pasca R: 'Nonsteroidal anti-inflammatory drugs activate carbonic anhydrase by a direct mechanism of action'. J Pharmacol Exp Ther. 1996, 277: 1464-1466.
Villanueva GB, Allen N: 'Acetylation of antithrombin III by aspirin'. Semin Thromb Hemost. 1986, 12: 213-215. 10.1055/s-2007-1003553.
Glasgow JF, Middleton B, Moore R, et al: 'The mechanism of inhibition of beta-oxidation by aspirin metabolites in skin fibroblasts from Reyes syndrome patients and controls'. Biochim Biophys Acta. 1999, 1454: 115-125. 10.1016/S0925-4439(99)00025-3.
Kawashima Y, Nakagawa S, Tachibana Y, Kozuka H: 'Effects of peroxisome proliferators on fatty acid-binding protein in rat liver'. Biochim Biophys Acta. 1983, 754: 21-27. 10.1016/0005-2760(83)90077-2.
Woodcock J, Moazed D, Cannon M, et al: 'Interaction of antibiotics with A- and P-site-specific bases in 16S ribosomal RNA'. EMBO J. 1991, 10: 3099-3103.
Rothlin CV, Katz E, Verbitsky M, et al: 'Block of the α9 nicotinic receptor by ototoxic aminoglycosides'. Neuropharmacology. 2000, 39: 2525-2532. 10.1016/S0028-3908(00)00056-3.
Ylikoski J, Xing-Qun L, Virkkala J, Pirvola U: 'Blockade of c-Jun N-terminal kinase pathway attenuates gentamicin-induced cochlear and vestibular hair cell death'. Hear Res. 2002, 163: 71-81. 10.1016/S0378-5955(01)00380-X.
Bryla J, Lietz T, Jarzyna R, et al: 'Differential in vivo and in vitro effect of gentamicin on glutamate synthesis and glutamate deamination in rabbit kidney-cortex tubules and mitochondria'. Pharmacol Res. 1992, 26: 367-375. 10.1016/1043-6618(92)90235-4.
Ghosh S, Dasgupta S: 'Gentamicin induced inhibition of steroidogenic enzymes in rat testis'. Indian J Physiol Pharmacol. 1999, 43: 247-250.
Chavez ML, DeKorte CJ: 'Valdecoxib: A review'. Clin Ther. 2003, 25: 817-851. 10.1016/S0149-2918(03)80110-8.
Bamforth KJ, Dalgliesh K, Coughtrie MW: 'Inhibition of human liver steroid sulfotransferase activities by drugs: A novel mechanism of drug toxicity?'. Eur J Pharmacol. 1992, 228: 15-21.
Goodwin JS: 'Mechanism of action of nonsteroidal anti-inflammatory agents'. Am J Med. 1984, 77: 57-64. 10.1016/S0002-9343(84)80020-0.
Martin JA, Mobberley MA, Redshaw S, et al: 'The inhibitory activity of a peptide derivative against the growth of simian immunode ficiency virus in C8166 cells'. Biochem Biophys Res Commun. 1991, 176: 180-188. 10.1016/0006-291X(91)90906-N.
Lenhard JM, Furfine ES, Jain RG, et al: 'HIV protease inhibitors block adipogenesis and increase lipolysis in vitro'. Antiviral Res. 2000, 47: 121-129. 10.1016/S0166-3542(00)00102-9.
Fourmy D, Recht MI, Blanchard SC, Puglisi JD: 'Structure of the A site of Escherichia coli 16S ribosomal RNA complexed with an aminoglycoside antibiotic'. Science. 1996, 274: 1367-1371. 10.1126/science.274.5291.1367.
von Ahsen U, Davies J, Schroeder R: 'Antibiotic inhibition of group I ribozyme function'. Nature. 1991, 353: 368-370. 10.1038/353368a0.
Mikkelsen NE, Brännvall M, Virtanen A, Kirsebom LA: 'Inhibition of RNase P RNA cleavage by aminoglycosides'. Proc Natl Acad Sci USA. 1999, 96: 6155-6160. 10.1073/pnas.96.11.6155.
Gabev E, Kasianowicz J, Abbott T, McLaughlin S: 'Binding of neomycin to phosphatidylinositol 4,5-bisphosphate (PIP2)'. Biochim Biophys Acta. 1989, 979: 105-112. 10.1016/0005-2736(89)90529-4.
Wang Y, Hamasaki K, Rando RR: 'Specificity of aminoglycoside binding to RNA constructs derived from the 16S rRNA decoding region and the HIV-RRE activator region'. Biochemistry. 1997, 36: 768-779. 10.1021/bi962095g.
Polidoro G, Del Boccio G, Di Ilio C, et al: 'In vitro interaction of penicillins and cephalosporins with human placenta GSH S-transferase'. Res Commun Chem Pathol Pharmacol. 1984, 46: 411-423.
Wright AJ: 'The penicillins'. Mayo Clin Proc. 1999, 74: 290-307. 10.4065/74.3.290.
Padovan E, von Greyerz S, Pichler WJ, Weltzien HU: 'Antigen-dependent and -independent IFN-γ modulation by penicillins'. J Immunol. 1999, 162: 1171-1177.
Ogawara H, Akiyama T, Ishida J, et al: 'A specific inhibitor for tyrosine protein kinase from Pseudomonas'. J Antibiot (Tokyo). 1986, 39: 606-608. 10.7164/antibiotics.39.606.
Akiyama T, Ishida J, Nakagawa S, et al: 'Genistein, a specific inhibitor of tyrosine-specific protein kinases'. J Biol Chem. 1987, 262: 5592-5595.
Traxler P, Green J, Mett H, et al: 'Use of a pharmaco-phore model for the design of EGFR tyrosine kinase inhibitors: isoflavones and 3-phenyl-4(1H)-quinolones'. J Med Chem. 1999, 42: 1018-1026. 10.1021/jm980551o.
Okura A, Arakawa H, Oka H, et al: 'Effect of genistein on topoisomerase activity and on the growth of [Val 12]Ha-ras-transformed NIH 3T3 cells'. Biochem Biophys Res Commun. 1988, 157: 183-189. 10.1016/S0006-291X(88)80030-5.
Markovits J, Linassier C, Fosse P, et al: 'Inhibitor effects of the tyrosine kinase inhibitor genistein on mammalian DNA topoisomerase II'. Cancer Res. 1989, 49: 5111-5117.
Pike ACW, Brzozowski AM, Hubbard RE, et al: 'Structure of the ligand-binding domain of oestrogen receptor β in the presence of a partial agonist and a full agonist'. EMBO J. 1999, 18: 4608-4618. 10.1093/emboj/18.17.4608.
Margeat E, Bourdoncle A, Margueron R, et al: 'Ligands differentially modulate the protein interactions of the human estrogen receptors α and β'. J Mol Biol. 2003, 326: 77-92. 10.1016/S0022-2836(02)01355-4.
Wolfson JS, Hooper DC: 'Fluoroquinolone antimicrobial agents'. Clin Microbiol Rev. 1989, 2: 378-424.
Heisig P: 'Genetic evidence for a role of parC mutations in development of high-level fluoroquinolone resistance in Escherichia coli'. Antimicrobial Agents Chemother. 1996, 40: 879-885.
Varon E, Janoir C, Kitzis M-D, Gutmann L: 'ParC and GyrA may be interchangeable initial targets of some fluoroquinolones in Streptococcus pneumoniae'. Antimicrobial Agents Chemother. 1999, 43: 302-306.
Hamy F, Felder ER, Heizmann G, et al: 'An inhibitor of the Tat/TAR RNA interaction that effectively suppresses HIV-1 replication'. Proc Natl Acad Sci USA. 1997, 94: 3548-3553. 10.1073/pnas.94.8.3548.
Daelemans D, Schols D, Witvrouw M, et al: 'A second target for the peptoid Tat/transactivation response element inhibitor CGP64222: inhibition of human immunodeficiency virus replication by blocking CXC0Chemokine receptor 4-mediated virus entry'. Mol Pharmacol. 2000, 57: 116-124.
Sedlacek HH: 'Mechanisms of action of flavopiridol'. Crit Rev Oncol Hematol. 2001, 38: 139-170. 10.1016/S1040-8428(00)00124-4.
Bible KC, Bible RH, Kottke TJ, et al: 'Flavopiridol binds to duplex DNA'. Cancer Res. 2000, 60: 2419-2428.
Maxwell A: 'DNA gyrase as a drug target'. Trends Microbiol. 1997, 5: 102-109. 10.1016/S0966-842X(96)10085-8.
Kranz RG, Beckman DL, Foster-Hartnett D: 'DNA gyrase activities from Rhodobacter capsulatus: Analysis of target(s) of coumarins and cloning of the gyrB locus'. FEMS Microbiol Lett. 1982, 72: 25-32.
Chen X, Ji ZL, Chen YZ: 'TTD: therapeutic target database'. Nucl Acids Res. 2002, 30: 412-415. 10.1093/nar/30.1.412.
McNeal J: 'Efforts in combinatorial chemistry'. Todays Chemist at Work. 1998, 7: 34-36.
DiMasi JA, Bryant NR, Lasagna L: 'New drug development in the United States from 1963 to 1990'. Clin Pharmacol Ther. 1991, 50: 471-486. 10.1038/clpt.1991.172.
Drews J: 'Strategic choices facing the pharmaceutical industry: A case for innovation'. Drug Discov Today. 1997, 2: 72-78. 10.1016/S1359-6446(96)10051-9.
Zivin JA: 'Understanding clinical trials'. Sci Am. 2000, 282: 69-75. 10.1038/scientificamerican0400-69.
Greene N: 'Computer systems for the prediction of toxicity: An update'. Adv Drug Deliv Rev. 2002, 54: 417-431. 10.1016/S0169-409X(02)00012-1.
Cronin MT: 'Computer-aided prediction of drug toxicity and metabolism'. EXS. 2003, 93: 259-278.
Wilson AG, White AC, Mueller RA: 'Role of predictive metabolism and toxicity modeling in drug discovery -- a summary of some recent advancements'. Curr Opin Drug Discov Devel. 2003, 6: 123-128.
Chen YZ, Ung CY: 'Prediction of potential toxicity and side effect protein targets of a small molecule by a ligand-protein inverse docking approach'. J Mol Graph Model. 2001, 20: 199-218. 10.1016/S1093-3263(01)00109-7.
Rodi DJ, Agoston GE, Manon R, et al: 'Identification of small molecule binding sites within proteins using phage display technology'. Comb Chem High Throughput Screen. 2001, 4: 553-572.
Sche PP, McKenzie KM, White JD, Austin DJ: 'Display cloning: Functional identification of natural product receptors using cDNA-phage display'. Chem Biol. 1999, 6: 707-716. 10.1016/S1074-5521(00)80018-6.
Sche PP, McKenzie KM, White JD, Austin DJ: 'Corrigendum to: 'Display cloning: functional identification of natural product receptors using cDNA-phage display' [Chem. Biol., Vol. 6, pp. 707-716]'. Chem Biol. 2001, 8: 399-400. 10.1016/S1074-5521(01)00012-6.
Jin Y, Yu J, Yu YG: 'Identification of hNopp140 as a binding partner for doxorubicin with a phage display cloning method'. Chem Biol. 2002, 9: 157-162. 10.1016/S1074-5521(02)00096-0.
Darvas F, Dorman G, Urge L, et al: 'Combinatorial chemistry. Facing the challenge of chemical genomics'. Pure Appl Chem. 2001, 73: 1487-1498. 10.1351/pac200173091487.
Newman DJ, Cragg GM, Snader KM: 'Natural products as sources of new drugs over the period 1981-2002'. J Nat Prod. 2003, 66: 1022-1037. 10.1021/np030096l.
Borman S: 'Organic lab sparks drug discovery'. Chem Eng News. 2002, 80: 23-24.
Phizicky E, Bastiaens PI, Zhu H, et al: 'Protein analysis on proteomic scale'. Nature. 2003, 422: 208-215. 10.1038/nature01512.
Graves PR, Kwiek JJ, Fadden P, et al: 'Discovery of novel targets of quinoline drugs in the human purine binding proteome'. Mol Pharmacol. 2002, 62: 1364-1372. 10.1124/mol.62.6.1364.
Takahishi TT, Austin RJ, Roberts RW: 'mRNA display: Ligand discovery, interaction analysis and beyond'. Trends Biochem Sci. 2003, 28: 159-165. 10.1016/S0968-0004(03)00036-7.
Lee SY, Choi JH, Xu Z: 'Microbial cell-surface display'. Trends Biotechnol. 2003, 21: 45-52. 10.1016/S0167-7799(02)00006-9.
Williams C: 'Biotechnology match making: Screening orphan ligands and receptors'. Curr Opin Biotechnol. 2003, 11: 42-46.
McPherson M, Yang YF, Hammond PW, Kreider BL: 'Drug receptor identification from multiple tissues using cellular-derived mRNA display libraries'. Chem Biol. 2002, 9: 691-698. 10.1016/S1074-5521(02)00148-5.
Scott JK, Barbas CF: 'Phage-display vectors'. Phage Display: A Laboratory Manual. Edited by: Barbas CF. 2001, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, 2.1-2.19.
Scott JK, Smith GP: 'Searching for peptide ligands with an epitope library'. Science. 1990, 249: 386-390. 10.1126/science.1696028.
Rodi DJ, Soares A, Makowski L: 'Quantitative assessment of peptide sequence diversity in M13 combinatorial peptide phage display libraries'. J Mol Biol. 2002, 322: 1039-1052. 10.1016/S0022-2836(02)00844-6.
Castagnoli L, Zucconi A, Quondam M, et al: 'Alternative bacteriophage display systems'. Comb Chem High Throughput Screen. 2001, 4: 121-133.
Santantelo GM, Tornow J, Moldave K: 'Cloning of open reading frames and promoters from the Saccharomyces cervisiae genome: Construction of genomic libraries of random small fragments'. Gene. 1986, 46: 181-186. 10.1016/0378-1119(86)90402-6.
Bluthner M, Bautz EKF, Bautz FA: 'Mapping of epitopes recognized by PM/Scl autoantibodies with gene-fragmented phage display libraries'. J Immunol Methods. 1996, 198: 187-198. 10.1016/S0022-1759(96)00160-3.
Jacobsson K, Frykberg L: 'Gene VIII-based phage-display vectors for selection against complex mixtures of ligands'. BioTechniques. 1998, 24: 294-301.
Cochrane D, Webster C, Masih G, McCafferty J: 'Identification of natural ligands for SH2 domains from a phage display cDNA library'. J Mol Biol. 2000, 297: 89-97. 10.1006/jmbi.2000.3561.
Robben J, Hertveldt K, Bosmans E, Volckaert G: 'Selection and identification of dense granule antigen GRA3 by Toxoplasma gondii whole genome phage display'. J Biol Chem. 2002, 277: 17544-17547. 10.1074/jbc.M110275200.
Hagiwara H, Kunihiro S, Nakajima K, et al: 'Affinity selection of DNA-binding proteins from yeast genomic DNA libraries by improved lambda phage display vector'. J Biochem (Tokyo). 2002, 132: 975-982. 10.1093/oxfordjournals.jbchem.a003312.
Yamamoto M, Kominato Y, Yamamoto F: 'Phage display cDNA cloning of protein with carbohydrate affinity'. Biochem Biophys Res Commun. 1999, 16: 194-199.
Tsai C-J, Ma B, Nussinov R: 'Folding and binding cascades: Shifts in energy landscapes'. Proc Natl Acad Sci USA. 1999, 96: 9970-9972. 10.1073/pnas.96.18.9970.
Tsai C-J, Kumar S, Ma B, Nussinov R: 'Folding funnels, binding funnels and protein function'. Protein Sci. 1999, 8: 1181-1190. 10.1110/ps.8.6.1181.
Clackson T, Wells JA: 'A hot spot of binding energy in a hormone-receptor interface'. Science. 1995, 267: 383-386. 10.1126/science.7529940.
Dwyer JJ, Dwyer MA, Kossiakoff AA: 'High affinity RNase S-peptide variants obtained by phage display have a novel hot-spot of binding energy'. Biochemistry. 2001, 40: 13491-13500. 10.1021/bi011703b.
Scott JK, Huang SF, Gangadhar BP, et al: 'Evidence that a protein-protein interaction hot spot on heterotrimeric G protein γ sub- units is used for recognition of a subclass of effectors'. EMBO J. 2001, 20: 767-776. 10.1093/emboj/20.4.767.
Freire E: 'The propagation of binding interactions to remote sites in proteins: Analysis of the binding of the monoclonal antibody D1.3 to lysozyme'. Proc Natl Acad Sci USA. 1999, 96: 10118-10122. 10.1073/pnas.96.18.10118.
Zhang RG, Pappas T, Brace JL, et al: 'Structure of a bacterial quorum-sensing transcription factor complexed with pheromone and DNA'. Nature. 2002, 417: 971-974. 10.1038/nature00833.
'Taxane anticancer agents. Basic science and current status'. ACS Symp Ser. Edited by: Georg GI, Chen TT, Ojima I, Vyas DM. 1995, Washington, DC, 583:
Stebbins CE, Russo AA, Schneider C, et al: 'Crystal structure of an Hsp90-geldamycin complex; targeting of a protein chaperone by an anti-tumor agent'. Cell. 1997, 89: 239-252. 10.1016/S0092-8674(00)80203-2.
Muchmore SW, Sattler M, Liang H, et al: 'X-ray and NMR structure of human Bcl-xL, an inhibitor of programmed cell death'. Nature. 1996, 381: 335-339. 10.1038/381335a0.
Petros AM, Medek A, Nettesheim DG, et al: 'Solution structure of the anti-apoptotic protein Bcl-2'. Proc Natl Acad Sci USA. 2001, 98: 3012-3016. 10.1073/pnas.041619798.
Maundrell K, Antonsson B, Magnenat E, et al: 'Bcl-2 undergoes phosphorylation by c-Jun N-terminal kinase/stress-activated protein kinases in the presence of the constitutively active GTP-binding protein Rac1'. J Biol Chem. 1997, 272: 25238-25242. 10.1074/jbc.272.40.25238.
Ito T, Deng X, Carr B, May WS: 'Bcl-2 phosphory-lation required for anti-apoptosis function'. J Biol Chem. 1997, 272: 11671-11673. 10.1074/jbc.272.18.11671.
Woods CM, Zhu J, McQueney PA, et al: 'Taxol-induced mitotic block triggers rapid onset of a p53-independent apoptotic pathway'. Mol Med. 1995, 1: 506-526.
Jordan MA, Wendell K, Gardiner S, et al: 'Mitotic block induced in HeLa cells by low concentrations of paclitaxel (Taxol) results in abnormal mitotic exit and apoptotic cell death'. Cancer Res. 1996, 56: 816-825.
Lieu CH, Yang YN, Lai YK: 'Dual cytotoxic mechanisms of submicromolar taxol on human leukemia HL-60 cells'. Biochem Pharmacol. 1997, 53: 1587-1596. 10.1016/S0006-2952(97)82450-1.
Andre N, Carre M, Brasseur G, et al: 'Paclitaxel targets mitochondria upstream of caspase activation in intact human neuroblastoma cells'. FEBS Lett. 2002, 532: 256-260. 10.1016/S0014-5793(02)03691-8.
Carre M, Andre N, Carles G, et al: 'Tubulin is an inherent component of mitochondrial membranes that interacts with the volta ge-dependent anion channel'. J Biol Chem. 2002, 277: 33664-33669. 10.1074/jbc.M203834200.
Srinivastava R, Mi QS, Hardwick JM, Longo DL: 'Deletion of the loop region of Bcl-2 completely clocks paclitaxel-induced apoptosis'. Proc Natl Acad Sci USA. 1999, 96: 3775-3780. 10.1073/pnas.96.7.3775.
Wang S, Wang Z, Biose L, et al: 'Loss of the Bcl-2 phos-phorylation loop domain increases resistance of human leukemia cells to paclitaxel-mediated mitochondrial dysfunction and apoptosis'. Biochem Biophys Res Commun. 1999, 259: 67-72. 10.1006/bbrc.1999.0669.
Haldar S, Basu A, Croce CM: 'Serine-70 is one of the critical sites for drug-induced Bcl2 phosphorylation in cancer cells'. Cancer Res. 1998, 58: 1609-1615.
Muchmore SW, Sattler M, Liang H, et al: 'X-ray and NMR structure of human Bcl-xL, an inhibitor of programmed cell death'. Nature. 1996, 381: 335-339. 10.1038/381335a0.
Petros AM, Medek A, Nettesheim DG, et al: 'Solution structure of the anti-apoptotic protein Bcl-2'. Proc Natl Acad Sci USA. 2001, 98: 3012-3016. 10.1073/pnas.041619798.
We thank B. Kay for comments on the manuscript and S. Mandava for continuing assistance with the similarity calculations. This work was funded by a grant from the Lucille P. Markey Foundation to L.M. Continued support for development of these methods comes from the Office of Biological and Environmental Research, Department of Energy under Contract No. W-31-109-Eng-38 to D.J.R.