- Research article
- Open Access
Combining docking with pharmacophore filtering for improved virtual screening
© Peach and Nicklaus; licensee BioMed Central Ltd. 2009
- Received: 13 May 2009
- Accepted: 20 May 2009
- Published: 20 May 2009
Virtual screening is used to distinguish potential leads from inactive compounds in a database of chemical samples. One method for accomplishing this is by docking compounds into the structure of a receptor binding site in order to rank-order compounds by the quality of the interactions they form with the receptor. It is generally established that docking can be reasonably successful at generating good poses of a ligand in an active site. However, the scoring functions that are used with docking are typically not successful at correctly ranking ligands according to binding affinity or even distinguishing correct poses of a given ligand from incorrect ones.
We have developed a simple method for reducing the number of false positives in a virtual screen, meaning ligands which are scored highly by the docking program but do not bind well in reality. This method uses a docking program for pose generation without regard to scoring, followed by filtering with receptor-based pharmacophore searches. We applied it to three test-case targets: neuraminidase A, cyclin-dependent kinase 2, and the C1 domain of protein kinase C.
The pharmacophore filtering method can perform better than more traditional docking + scoring methods, and allows the advantages of both docking-based and pharmacophore-based approaches to virtual screening to be fully realized.
- Virtual Screening
- Hydrogen Bond Donor
- Hydrogen Bond Acceptor
- Pharmacophore Model
The goal of virtual screening is to select, relatively rapidly and cheaply, a small subset of compounds predicted to have activity against a given biological target out of a large database of compounds. While it is possible to screen large databases in their entirety using automated high-throughput screening methods, this is expensive and requires a substantial investment in infrastructure and assay development. The idea of virtual screening is to test compounds computationally in order to reduce the number of compounds to be screened experimentally, with the additional advantage that the number of compounds in the final set can easily be adjusted according to the resources available for assaying.
The database used for virtual screening can be a collection of commercially available compounds, such as ZINC  or the ChemNavigator iResearch Library , both of which are meta-collections of supplier catalogs. Pharmaceutical companies typically have an in-house database of previously synthesized molecules. A publicly available alternative to this is the open NCI database , a collection of compounds that have been tested over the past few decades in the National Cancer Institute's screens for anticancer activity. Small samples from a subset of this collection are available for research purposes upon request .
A variety of computational methods can be used for virtual screening depending on the desired size of the final subset and on the amount of information known about the target, its natural ligands, and any known inhibitors. Here we focus on the method of receptor-ligand docking and scoring, which can be used when a three-dimensional structure of the target is available. This method can be divided into two parts: first docking to position ligand structures into the target binding site, generating a set of poses for each ligand; and secondly scoring to evaluate and rank-order poses and ligands according to how well each pose for each ligand fits into the binding site and the quality of the interactions it forms with the target. Generally, in these methods the ligand has conformational flexibility while the receptor remains essentially rigid. The output from the docking program is thus a set of poses saved for each ligand, with a numerical score for each pose.
The main problem with virtual screening is that many, and in some cases the vast majority, of the compounds that are predicted to be active are in fact not active when screened experimentally. There are two theories found in the literature on the reasons for this phenomenon. Some researchers argue that docking can usually generate good poses of a ligand in an active site, however scoring functions are generally not successful at correctly ranking either ligands or poses . Others have pointed out that docking programs do not always generate correct poses, and that the highest ranked pose for a given ligand is often incorrect. Scoring functions would therefore function much better at ranking if docking programs did not produce so many incorrect poses for each compound [6, 7].
Regardless of whether it is the docking programs or the scoring functions that are at fault, the issue is that virtual screening can generate an enormous number of false positives – compounds that are scored highly in silico but do not actually bind to the target in vivo or in vitro. These false positives can also be blamed to some extent for false negatives, in the case where true positives are scored relatively poorly (and perhaps even eliminated) because of spuriously high-scoring false positives that are ranked ahead of them. A method of eliminating at least some of these false positives at some stage in the virtual screening procedure would therefore be very useful.
Here we present such a method, which we call pharmacophore filtering. It is a means of post-processing docking results to rapidly eliminate poses and molecules that are not fully chemically compatible with the binding site. This includes, for example, poses that do not completely fill the site, or that leave unpaired buried hydrogen bond donors or acceptors. This method could be viewed as an enforcement of the basic principle of structure-based drug design, namely that good-binding ligands must be chemically complementary to their receptors.
The advantages of including information about the target, such as specific, required hydrogen bonds into docking simulations is already well-appreciated . One advantage of our method over others is that it allows multiple such target-ligand interactions to be quickly tested, compared and re-adjusted without re-running the entire docking calculation. Thus it is a useful addition to the virtual screening arsenal, and as we attempt to demonstrate here, it is broadly applicable to a variety of targets and ligand design goals.
Pharmacophore Filtering Method Description
The pharmacophore filtering method begins by using a docking program in the ordinary way, for pose generation and alignment in the binding site. All the poses output by the docking program are saved to a file, while their scores and rankings are ignored. Next, the docked poses are searched and filtered using a series of pharmacophore query models. The pharmacophores are elucidated based on available crystal structures of bound ligands and/or simple examination of the binding site, and are used to filter the saved poses output from the docking program, and remove any that are incompatible with the model. This method thus adds an element of ligand-based drug design to what is essentially a structure-based drug design method of docking and scoring.
Clearly this method works best when at least one co-crystal structure with the natural ligand or an inhibitor is available, but it would also in theory be possible to use an apo structure and examine the interactions made with water. There are no special restrictions on what software to use, though it is probably better if the docking program has a stochastic component and can generate diverse poses rather than converging on a "best solution." In the test cases for this method, described below, we used the well-known docking programs GOLD  and Glide .
Rather than using a full-blown pharmacophore generation and search program, it is possible to simply filter the output file of docked poses according to whether or not they fulfill certain receptor-ligand contacts or interactions. For example, with Glide, Schrödinger provides a "Pose-Filter" Python script to perform this type of filtering, and for GOLD, such an analysis can be done with its companion programs Hermes and GoldMine. However, a dedicated pharmacophore program will provide more options and greater flexibility in defining the filters to be used.
The pharmacophore program must allow the import of a set of pre-generated conformers that are pre-aligned to the pharmacophore model. The site points in the pharmacophore model can be defined based on either the positions of ligand atoms (ligand-sided) or the positions of receptor atoms (protein-sided), or both, depending on the nature of the interaction to be captured with the filter. The radius of the site points can also be adjusted to alter the sensitivity of the model. A smaller radius gives a tighter, more selective filter, whereas a larger radius can capture some flexibility of ligands within the binding site or account for some variability in binding modes. For our test cases the pharmacophore models were generated in MOE  by simple visual inspection of crystal structure binding sites and co-crystallized ligands, along with information on the binding modes of other known ligands from the literature.
Alternately, there are several published methods and software programs available for automatically or semi-automatically generating structure-based pharmacophore models. For example, the program LUDI  (currently commercially available as part of Discovery Studio from Accelrys) calculates an interaction map of locations in the receptor binding site where an atom from a bound ligand would be in position to form a favorable hydrogen-bonding or hydrophobic interaction. This map can then be converted into a pharmacophore model. Another such program is LigandScout , which is a fully automated method of generating a pharmacophore model from a set of protein-ligand complexes.
We evaluated our methodology with three test cases: protein targets of pharmaceutical relevance with a variety of binding site characteristics.
Influenza neuraminidase A is a surface glycoprotein of the influenza virus whose function is to cleave the linkage between sialic acid and an adjacent sugar in glycoconjugates on the surface of cells targeted for infection . The sialic acid binding site is small, deep, and highly polar. For this target, there existed a literature reference examining which of the many available crystal structures of neuraminidase is best suited for docking a variety of ligands . Based on this study's conclusions, we chose the 1MWE crystal structure , which has good (1.7 Å) resolution and could accommodate all members of a set of co-crystallized ligands . To prepare the structure for docking, we deleted all bound ligands and crystallographic waters and flipped the orientation of residue Asn 294 in the binding site (the sidechain amide nitrogen and oxygen had been assigned incorrectly).
The screening database for this test case was constructed by first downloading a set of 245 neuraminidase ligands from the BindingDB [17, 18]. Ligands with Ki or IC50 values of 10 μM or lower (195 compounds) were considered actives. The inactive ligands were added to a set of 4775 decoys obtained from both the MDL Drug Data Report database (MDDR)  and the ChemNavigator iResearch Library (iRL)  of commercially available compounds. The decoys were filtered using Pipeline Pilot  to choose compounds whose physicochemical properties fit into the ranges seen with the known neuraminidase ligands. This was to ensure that the docking and scoring protocols were not biased in distinguishing hits from decoys due to differences in their property distributions. In this case, decoys were required to have a molecular weight between 190 and 500, a logP between -7.5 and 4.0, a polar surface area between 95 and 250 Å2, fewer than 10 rotatable bonds, at least 3 hydrogen bond donors and at least 2 hydrogen bond acceptors. Thus the final set of 4775 decoys were all as small and highly polar as the set of known neuraminidase binders.
The screening database was docked into the binding site using two docking programs, GOLD 3.0.1  and Glide 4.0 . In GOLD, the binding site was defined as a sphere with a radius of 10 Å centered at the position of atom C6 in the bound sialic acid ligand. We used the 7–8x speedup rather than the library screening settings for the genetic algorithm, for some sacrifice in speed but an improvement in the quality of the generated poses, and the Goldscore scoring function. In Glide, the protein structure was prepared using the Protein Preparation and Grid Generation modules, with all default settings. For docking we used the HTVS precision mode, with all other parameters left in their default settings. With both programs, ten poses were saved for each compound, then all the saved poses from both programs were submitted to pharmacophore filtering (see below). For comparison purposes, all the saved poses from both programs were also re-scored with a total of four different scoring functions: GScore in Glide, Goldscore and Chemscore as implemented in GOLD, and an "Affinity" score. This latter score is a reformulation of two of the terms in the Goldscore scoring function, defined as hbond.external + 1.375*vdw.external, and has been shown to correlate better with experimental binding affinities than Goldscore itself, which was optimized to give good poses .
Another sub-pocket of the binding site holds an acetyl group via a hydrogen bond from the carbonyl oxygen to Arg 152, and a hydrophobic interaction of the methyl group with residues Ile 222 and Trp 178 . We constructed a second filter for this interaction with a ligand-sided hydrogen bond acceptor or anionic site point and a ligand-sided hydrophobic site point (Figure 2B).
Finally, the third filter took into consideration the three highly conserved glutamate residues with unpaired acceptor atoms across the center of the binding site (Glu 119, 227, and 276) . The superimposed co-crystal ligands showed a variety of methods of interacting with these glutamates, so we built receptor-sided site points (meaning that they were centered on oxygen atoms in the glutamate residues of the receptor) for hydrogen bond donor or positively charged atoms, with a radius of about 3 Å to allow for the hydrogen bonding distance to ligand atoms (Figure 2C). The filter required an interaction with at least one of the glutamates.
Cyclin-dependent kinase 2 (CDK2)
The family of cyclin-dependent kinases are regulators of the cell cycle. In particular CDK2, in complex with cyclins E and A, is involved in the G1-S transition and the progression through S phase, although it has recently been shown not to be essential for mitotic cell division in mice . Like all other kinases, CDK2 has an ATP binding site located in a cleft between the N-terminal and C-terminal lobes of the protein. Here again, as with neuraminidase, we consulted the literature to determine which of the available crystal structures would be best suited for virtual screening. There is a published study examining a diverse subset of 20 of the available CDK2 crystal structures, and we chose the one that performed best at docking the most structurally diverse group of CDK2 inhibitors . This crystal structure, 1OIT, is in an active, open conformation and has a high crystallographic resolution of 1.6 Å . Several loop regions in the structure were disordered and coordinates for these were not present. We capped all such protein chain ends, including the N- and C-termini, with hydrogens to form neutral NH2 and COOH groups. We also built coordinates for the missing sidechain of residue Lys 9, which is near the binding site, in a conformation pointing out into solvent with a χ1 angle of -60°.
Along the same lines as with neuraminidase, the screening database was constructed by first downloading a set of 1278 CDK2 ligands from the BindingDB [17, 18]. Ligands with Ki or IC50 values of 10 μM or lower (1063 compounds) were considered active, and the inactive ligands were added to a set of 26792 decoys combined from the MDDR  and ChemNavigator iRL  databases. These decoys were filtered using Pipeline Pilot  to choose compounds physicochemically similar to known kinase inhibitors, with a molecular weight between 200 and 550, at least two aromatic rings, a logP between 0.5 and 6.0, a polar surface area of less than 150 Å2, fewer than 10 rotatable bonds, less than 4 hydrogen bond donors and between 2 and 8 hydrogen bond acceptors.
The screening database was docked into the binding site using both GOLD  and Glide , with the same settings used as before with neuraminidase. In GOLD, the binding site was defined as a sphere with a radius of 10 Å centered at the position of atom C14 in the bound ligand. Ten docked poses were saved for each compound from each program, and all the poses were again submitted to pharmacophore filtering. All the saved poses were also re-scored as before with the four scoring functions: GScore, Goldscore, Chemscore, and GOLD Affinity.
Protein kinase C (PKC) C1 domain
Protein kinase C isozymes are a family of serine/threonine kinases which are centrally involved in cell signaling. The regulatory C1 domain is a 50-residue zinc-finger-like structure that responds to the second messenger diacylglycerol (DAG) by translocating PKC to the cell membrane, where the whole protein undergoes a conformational change and the kinase domain becomes activated . This test case was a more difficult target. The natural activator is a lipid, and other known ligands are all also lipophilic natural products  and derivatives thereof . Furthermore, the binding site in the C1 domain is a half-site in that ligands make interactions with the C1 domain itself, but interactions with lipids in the cell membrane are also important for binding . There is only one holo C1 domain crystal structure (1PTR) available in the PDB, with phorbol ester (a natural tumor promoter and DAG mimetic) in the binding site .
The screening database consisted of a small collection of 27 known C1 domain ligands [29, 30, 33], along with a set of 1096 natural product decoys from several databases: the Natural Products subset of the Open NCI Database , a database of compounds used in traditional Chinese medicine , and a database of natural products isolated from marine species . As with the other two test cases, we selected decoys using Pipeline Pilot  that were physicochemically similar to the known ligands. Decoys were required to have a molecular weight greater than 250, a logP between 0.5 and 6, at least three hydrogen bond acceptors and at least one hydrogen bond donor, less than 250 Å2 of polar surface area, and more than 220 Å2 of non-polar surface area. High-quality three-dimensional structures for these often-complex molecules were generated using CORINA .
The screening database was docked, as before, with GOLD  and Glide . Due to the fact that the natural products tended to be large compounds with many rotatable bonds, we used standard precision docking in Glide and the default genetic algorithm settings in GOLD rather than high-throughput or speeded-up screening settings, and we increased the sampling by saving 20 poses for each compound. In GOLD, the binding site was defined as a sphere with a radius of 10 Å centered at the position of the Nε atom in residue Gln 257. All saved poses were submitted to pharmacophore filtering, and also re-scored as before with the four scoring functions: GScore, Goldscore, Chemscore, and GOLD Affinity.
In all three of the test cases at least one crystal structure of the target and a collection of known inhibitors were available, allowing us to use properties calculated for known inhibitors to make rational decisions on the ranges of values for physicochemical properties that are reasonable for new potential inhibitors of the target. We used the structural characteristics of the binding site and the observed interactions between target and ligands to choose the parameters used in docking and scoring in order to reproduce known inhibitor binding modes, and for analysis of the docked poses with the pharmacophore filtering method. We have shown here that, compared to simply ranking the poses output by the docking program according to a single docking score, an initial filtering of the docked poses with a pharmacophore query to remove poses that do not form certain essential interactions with the target binding site greatly improves the quality of the results. Typically we have seen the ability to eliminate well over 90% of the decoy molecules in the starting database while retaining about 80% of the true positives.
This pharmacophore filtering method is in many ways similar to other approaches that attempt to add additional information to improve the results from docking and scoring. We will briefly delineate them, and then discuss what we believe are the advantages of the method presented here. One alternative way to incorporate the use of pharmacophores into a virtual screening protocol is to use them for pre-screening, rather than post-screening, the docking database. Such a strategy has been used successfully with several targets [39, 40], and has the advantage of reducing the size of the database to be docked. A variation on this method is called pharmacophore-based molecular docking, or PhDOCK , as it is a variation of the DOCK program. A database of conformers is overlaid according to 3D pharmacophores and the pharmacophore site points are then matched with docking site points in the binding site. This allows the simultaneous docking of many molecules at once. The docking site points can either be generated from standard DOCK "spheres" or can be based on the positions of atoms in crystallographic ligands. The advantage of these pre-screening methods is that they can significantly shorten the amount of computer time required for docking, though this becomes less and less relevant as computer speed increases and computer clusters become more common.
A second method of incorporating pharmacophores into virtual screening is to post-process docking results with various kinds of ligand interaction-based filters. One such method is the structural interaction fingerprint , or molecular interaction fingerprint , which encodes into a binary bit string the set of binding interactions made by a docked ligand pose. These strings can be clustered and compared using Tanimoto similarities, and used to rapidly filter poses for the presence or absence of specific interactions, providing an improvement in virtual screening hit list enrichment over conventional scoring functions. Similarly, an interactions-based accuracy classification method  has been shown to provide a superior assessment of docking pose quality compared to simple RMS deviations from crystal structure poses. Rather than applying pharmacophores to filter docked poses directly, the AutoShim method  uses point pharmacophore interaction features in the binding site to weight or "shim" the docking scoring function according to experimental IC50 data for the target receptor, and can produce large improvements in the ability of the scoring function to predict binding affinities. In contrast, the pharmacophore filtering method presented here requires no additional new or proprietary software unlike these other methods for post-processing docking results – only a docking program and a pharmacophore generation program, both of which are highly likely to be in use by any modeling group.
Finally, a third method of incorporating pharmacophores into virtual screening is to add constraints to the docking run to ensure that certain interactions between the ligands and the target binding site are formed. Such pharmacophoric constraints can be implemented in GOLD  as well as in many other modern docking programs, and are the basis for the screening program FlexX-Pharm , where ligands are incrementally constructed into the active site in a manner incorporating "look-ahead checks" to ensure that desired interactions are formed. These types of constraints are particularly useful with kinases, in which hydrogen bonds between ligands and the hinge region of the binding site are almost universally present . Adding this kind of additional information to the docking has been shown to improve enrichment, though this may come at the expense of a reduction in ability to identify novel chemotypes .
Although we have not compared our method directly with these constrained docking methods, we would anticipate that the enrichments seen would be quite similar (if not better with constrained docking because the conformational search done by the docking program is biased in favor of solutions that fit the constraints). However, the great advantage of the pharmacophore filtering method is that it allows the modeler to change his or her mind about which features are most important for a productive binding interaction without having to go back and redo the docking runs. A typical pharmacophore filtering run on a file of 100,000 docked poses takes only a few minutes, whereas the docking runs themselves can take days or weeks. It is therefore possible to quickly test multiple filtering combinations and variations, and to develop refined hit sets focusing on different regions of the binding site.
It is also possible to adjust the stringency of the filters to tune the number of compounds that are output from the pharmacophore filtering, to be sent on to the next step in the screening. If the hit set is to be assayed experimentally, the number of compounds can be scaled according to the resources available and the throughput level of the assay. If more detailed calculations are to be done, such as MM/PBSA scoring or other free energy estimations of binding affinity, the size of the hit set can be adjusted according to the available computational resources and time.
A final advantage of our pharmacophore filtering method is the ability to easily go back and look for compounds occupying specific sub-pockets in the binding site, as for example with CDK2 and residue Lys 89 on the top edge of the binding site. It is also possible to explore novel binding modes, to look for compounds that interact with the target receptor in ways that are not seen with existing known ligands.
The three test cases presented here are of course not sufficient to provide an estimate of the statistical significance of the improvement in performance of this method over traditional docking + scoring methods, and in fact for two of the targets (neuraminidase and CDK2) the improvement seen with the pharmacophore filtering method versus Glide docking with compounds ranked by Gscore is very slight. However, there is no way to know a priori whether or not a given docking + scoring method will perform well for a given target. Many published studies comparing docking programs to one another have shown widely varying results for different targets [5, 48, 49]. In this study, although Glide performs very well with neuraminidase and CDK-2, GOLD is more successful with the third target, the PKC C1 domain.
We also believe that human oversight and intervention in any computational modeling work is essential. In virtual screening, time spent thinking carefully about the receptor binding site and interactions made by ligands binding in it would pay off in a greater understanding of the structural characteristics of the system – knowledge that could be useful for subsequent refinement into leads of the hits out of the screening database. Our rapid and resource-wise undemanding method facilitates this thoughtful approach to one of the challenges in drug design.
In summary, combining docking with pharmacophore searches as a post-processing filter allows the advantages of both methods to be fully exploited. The docking program is used to fit the ligands into the binding site in as wide a variety as possible of reasonable ways. Put another way, we are using it as a conformation generator for conformations that fit into the active site. Compared to a traditional pharmacophore search, where conformations for each molecule in the database are pre-generated to be dispersed over all low-energy conformational space, the docking program allows us to focus on a set of conformations for each molecule that all fit into the binding site. Additionally, the docking program performs the alignment of all the conformations to one another, which is much simpler than dealing with the combinatorial explosion of ways in which the active compounds can be aligned and the number of features to include. Post-processing the docking results with pharmacophore filtering allows us to bypass to a large extent the difficulty with scoring functions, which is that while they are generally good at producing reasonable docked poses of a molecule in a binding site, they are not necessarily good at discriminating between good binders and poor binders, probably due to the non-linear nature of ligand-receptor recognition. The pharmacophore filtering method greatly increases the likelihood that the best and/or most correct pose is selected from the set of docked poses, regardless of the numerical value of its docking score. This method has been used successfully in several virtual screening projects in our laboratory, one for inhibitors of Met tyrosine kinase , and others which will be reported separately.
This research was supported by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research. This project has been funded with federal funds from the National Cancer Institute, National Institutes of Health, under contracts N01-CO-12400 and HHSN261-2008-00001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.
- Irwin JJ, Shoichet BK: ZINC – A free database of commercially available compounds for virtual screening. J Chem Inf Model. 2005, 45: 177-182. 10.1021/ci049714+.View ArticleGoogle Scholar
- ChemNavigator iResearch Library. [http://www.chemnavigator.com/cnc/products/iRL.asp]
- Milne GWA, Nicklaus MC, Driscoll JS, Wang S, Zaharevitz D: National Cancer Institute Drug Information System 3D database. J Chem Inf Comput Sci. 1994, 34: 1219-1224.View ArticleGoogle Scholar
- The NCI/DTP Open Chemical Repository Collection. [http://dtp.nci.nih.gov/branches/dscb/repo_open.html]
- Warren GL, Andrews CW, Capelli AM, Clarke B, LaLonde J, Lambert MH, Lindvall M, Nevins N, Semus SF, Senger S, Tedesco G, Wall ID, Woolven JM, Peishoff CE, Head MS: A critical assessment of docking programs and scoring functions. J Med Chem. 2006, 49: 5912-5931. 10.1021/jm050362n.View ArticleGoogle Scholar
- Kontoyianni M, Sokol GS, McClellan LM: Evaluation of library ranking efficacy in virtual screening. J Comput Chem. 2005, 26 (1): 11-22. 10.1002/jcc.20141.View ArticleGoogle Scholar
- Perola E: Minimizing false positives in kinase virtual screens. Proteins. 2006, 64: 422-435. 10.1002/prot.21002.View ArticleGoogle Scholar
- Jansen JM, Martin EJ: Target-biased scoring approaches and expert systems in structure-based virtual screening. Curr Opin Chem Biol. 2004, 8: 359-364. 10.1016/j.cbpa.2004.06.002.View ArticleGoogle Scholar
- Jones G, Willett P, Glen RC, Leach AR, Taylor R: Development and validation of a genetic algorithm for flexible docking. J Mol Biol. 1997, 267: 727-748. 10.1006/jmbi.1996.0897.View ArticleGoogle Scholar
- Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, Repasky MP, Knoll EH, Shelley M, Perry JK, Shaw DE, Francis P, Shenkin PS: Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem. 2004, 47: 1739-1749. 10.1021/jm0306430.View ArticleGoogle Scholar
- MOE: Molecular Operating Environment version 2005.06. 2005, Chemical Computing Group, Inc.: Montreal, CanadaGoogle Scholar
- Böhm HJ: The computer program LUDI: A new method for the de novo design of enzyme inhibitors. J Comput Aided Mol Des. 1992, 6: 61-78. 10.1007/BF00124387.View ArticleGoogle Scholar
- Wolber G, Langer T: LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters. J Chem Inf Model. 2005, 45: 160-169. 10.1021/ci049885e.View ArticleGoogle Scholar
- Colman PM: Influenza virus neuraminidase: structure, antibodies, and inhibitors. Protein Sci. 1994, 3: 1687-1696. 10.1002/pro.5560031007.View ArticleGoogle Scholar
- Birch L, Murray CW, Hartshorn MJ, Tickle IJ, Verdonk ML: Sensitivity of molecular docking to induced fit effects in influenza virus neuraminidase. J Comput Aided Mol Des. 2002, 16: 855-869. 10.1023/A:1023844626572.View ArticleGoogle Scholar
- Varghese JN, Colman PM, van Donkelaar A, Blick TJ, Sahasrabudhe A, McKimm-Breschkin JL: Structural evidence for a second sialic acid binding site in avian influenza virus neuraminidases. Proc Natl Acad Sci USA. 1997, 94: 11808-11812. 10.1073/pnas.94.22.11808.View ArticleGoogle Scholar
- The Binding Database. [http://www.bindingdb.org]
- Liu T, Lin Y, Wen X, Jorissen RN, Gilson MK: BindingDB: a web-accessible database of experimentally determined protein-ligand binding affinities. Nucleic Acids Res. 2007, 35: D198-201. 10.1093/nar/gkl999.View ArticleGoogle Scholar
- Symyx Databases: Current bioactivity findings for newly launched and developmental drugs. [http://www.symyx.com/products/databases/bioactivity/mddr/index.jsp]
- Pipeline Pilot version 5.0. 2005, SciTegic, Inc.: San Diego, CAGoogle Scholar
- GOLD Support: Scientific FAQs. [http://www.ccdc.cam.ac.uk/products/life_sciences/gold/faqs/scientific_faq.php#bindingaffinity]
- Delano WL: The PyMOL Molecular Graphics System version 0.99. 2002, DeLano Scientific: Palo Alto, CAGoogle Scholar
- Malumbres M, Barbacid M: Mammalian cyclin-dependent kinases. Trends Biochem Sci. 2005, 30: 630-641. 10.1016/j.tibs.2005.09.005.View ArticleGoogle Scholar
- Thomas MP, McInnes C, Fischer PM: Protein structures in virtual screening: a case study with CDK2. J Med Chem. 2006, 49: 92-104. 10.1021/jm050554i.View ArticleGoogle Scholar
- Anderson M, Beattie JF, Breault GA, Breed J, Byth KF, Culshaw JD, Ellston RP, Green S, Minshull CA, Norman RA, Pauptit RA, Stanway J, Thomas AP, Jewsbury PJ: Imidazo[1,2-a]pyridines: a potent and selective class of cyclin-dependent kinase inhibitors identified through structure-based hybridisation. Bioorg Med Chem Lett. 2003, 13: 3021-3026. 10.1016/S0960-894X(03)00638-3.View ArticleGoogle Scholar
- Bramson HN, Corona J, Davis ST, Dickerson SH, Edelstein M, Frye SV, Gampe RT, Harris PA, Hassell A, Holmes WD, Hunter RN, Lackey KE, Lovejoy B, Luzzio MJ, Montana V, Rocque WJ, Rusnak D, Shewchuk L, Veal JM, Walker DH, Kuyper LF: Oxindole-based inhibitors of cyclin-dependent kinase 2 (CDK2): design, synthesis, enzymatic activities, and X-ray crystallographic analysis. J Med Chem. 2001, 44: 4339-4358. 10.1021/jm010117d.View ArticleGoogle Scholar
- Byth KF, Cooper N, Culshaw JD, Heaton DW, Oakes SE, Minshull CA, Norman RA, Pauptit RA, Tucker JA, Breed J, Pannifer A, Rowsell S, Stanway JJ, Valentine AL, Thomas AP: Imidazo[1,2-b]pyridazines: a potent and selective class of cyclin-dependent kinase inhibitors. Bioorg Med Chem Lett. 2004, 14: 2249-2252. 10.1016/j.bmcl.2004.02.008.View ArticleGoogle Scholar
- Newton AC: Protein kinase C: Structural and spatial recognition by phosphorylation, cofactors, and macromolecular interactions. Chem Rev. 2001, 101: 2353-2364. 10.1021/cr0002801.View ArticleGoogle Scholar
- Kishi Y, Rando RR: Structural basis of protein kinase C activation by tumor promoters. Acc Chem Res. 1998, 31: 163-172. 10.1021/ar9600751.View ArticleGoogle Scholar
- Marquez VE, Blumberg PM: Synthetic diacylglycerols (DAG) and DAG-lactones as activators of protein kinase C (PK-C). Acc Chem Res. 2003, 36: 434-443. 10.1021/ar020124b.View ArticleGoogle Scholar
- Kang JH, Peach ML, Pu Y, Lewin NE, Nicklaus MC, Blumberg PM, Marquez VE: Conformationally constrained analogues of diacylglycerol (DAG). 25. Exploration of the sn-1 and sn-2 carbonyl functionality reveals the essential role of the sn-1 carbonyl at the lipid interface in the binding of DAG-lactones to protein kinase C. J Med Chem. 2005, 48: 5738-5748. 10.1021/jm050352m.View ArticleGoogle Scholar
- Zhang G, Kazanietz MG, Blumberg P, Hurley JH: Crystal structure of the cys2 activator-binding domain of protein kinase Cδ in complex with phorbol ester. Cell. 1995, 81: 917-924. 10.1016/0092-8674(95)90011-X.View ArticleGoogle Scholar
- Shao L, Lewin NE, Lorenzo PS, Hu Z, Enyedy IJ, Garfield SH, Stone JC, Marner FJ, Blumberg PM, Wang S: Iridals are a novel class of ligands for phorbol ester receptors with modest selectivity for the RasGRP receptor subfamily. J Med Chem. 2001, 44: 3872-3880. 10.1021/jm010258f.View ArticleGoogle Scholar
- DTP – Natural Products Set Information. [http://dtp.nci.nih.gov/branches/dscb/natprod_explanation.html]
- Fang X, Shao L, Zhang H, Wang S: CHMIS-C: A comprehensive herbal medicine information system for cancer. J Med Chem. 2005, 48: 1481-1488. 10.1021/jm049838d.View ArticleGoogle Scholar
- Lei J, Zhou J: A marine natural product database. J Chem Inf Comput Sci. 2002, 42: 742-748.View ArticleGoogle Scholar
- Gasteiger J, Rudolph C, Sadowski J: Automatic generation of 3D-atomic coordinates for organic molecules. Tetrahedron Comput Methodol. 1990, 3: 547-547. 10.1016/0898-5529(90)90156-3.View ArticleGoogle Scholar
- Pak Y, Enyedy IJ, Varady J, Kung JW, Lorenzo PS, Blumberg PM, Wang S: Structural basis of binding of high-affinity ligands to protein kinase C: Prediction of the binding modes through a new molecular dynamics method and evaluation by site-directed mutagenesis. J Med Chem. 2001, 44: 1690-1701. 10.1021/jm000488e.View ArticleGoogle Scholar
- Lyne PD, Kenny PW, Cosgrove DA, Deng C, Zabludoff S, Wendoloski JJ, Ashwell S: Identification of compounds with nanomolar binding affinity for checkpoint kinase-1 using knowledge-based virtual screening. J Med Chem. 2004, 47: 1962-1968. 10.1021/jm030504i.View ArticleGoogle Scholar
- Brenk R, Naerum L, Grädler U, Gerber HD, Garcia GA, Reuter K, Stubbs MT, Klebe G: Virtual screening for submicromolar leads of tRNA-guanine transglycosylase based on a new unexpected binding mode detected by crystal structure analysis. J Med Chem. 2003, 46: 1133-1143. 10.1021/jm0209937.View ArticleGoogle Scholar
- Joseph-McCarthy D, Thomas BE, Belmarsh M, Moustakas D, Alvarez JC: Pharmacophore-based molecular docking to account for ligand flexibility. Proteins. 2003, 51 (2): 172-188. 10.1002/prot.10266.View ArticleGoogle Scholar
- Deng Z, Chuaqui C, Singh J: Structural interaction fingerprint (SIFt): a novel method for analyzing three-dimensional protein-ligand binding interactions. J Med Chem. 2004, 47: 337-344. 10.1021/jm030331x.View ArticleGoogle Scholar
- Marcou G, Rognan D: Optimizing fragment and scaffold docking by use of molecular interaction fingerprints. J Chem Inf Model. 2007, 47: 195-207. 10.1021/ci600342e.View ArticleGoogle Scholar
- Kroemer RT, Vulpetti A, McDonald JJ, Rohrer DC, Trosset JY, Giordanetto F, Cotesta S, McMartin C, Kihlén M, Stouten PFW: Assessment of docking poses: interactions-based accuracy classification (IBAC) versus crystal structure deviations. J Chem Inf Comput Sci. 2004, 44: 871-881.View ArticleGoogle Scholar
- Martin EJ, Sullivan DC: AutoShim: Empirically corrected scoring functions for quantitative docking with a crystal structure and IC50 training data. J Chem Inf Model. 2008, 48: 861-872. 10.1021/ci7004548.View ArticleGoogle Scholar
- Verdonk ML, Berdini V, Hartshorn MJ, Mooij WTM, Murray CW, Taylor RD, Watson P: Virtual screening using protein-ligand docking: avoiding artificial enrichment. J Chem Inf Comput Sci. 2004, 44: 793-806.View ArticleGoogle Scholar
- Hindle SA, Rarey M, Buning C, Lengauer T: Flexible docking under pharmacophore type constraints. J Comput Aided Mol Des. 2002, 16: 129-149. 10.1023/A:1016399411208.View ArticleGoogle Scholar
- Kontoyianni M, McClellan LM, Sokol GS: Evaluation of docking performance: comparative data on docking algorithms. J Med Chem. 2004, 47: 558-565. 10.1021/jm0302997.View ArticleGoogle Scholar
- Kellenberger E, Rodrigo J, Rognan D: Comparative evaluation of eight docking tools for docking and virtual screening accuracy. Proteins. 2004, 57: 225-242. 10.1002/prot.20149.View ArticleGoogle Scholar
- Peach ML, Tan N, Choyke SJ, Giubellino A, Athauda G, Burke TR, Nicklaus MC, Bottaro DP: Directed discovery of agents targeting the Met tyrosine kinase domain by virtual screening. J Med Chem. 2009, 52: 943-951. 10.1021/jm800791f.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.