Molecular docking and pharmacophore filtering in the discovery of dual-inhibitors for human leukotriene A4 hydrolase and leukotriene C4 synthase

Combination of drugs for multiple targets has been a standard treatment in treating various diseases. A single chemical entity that acts upon multiple targets is emerging nowadays because of their predictable pharmacokinetic and pharmacodynamic properties. We have employed a computer-aided methodology combining molecular docking and pharmacophore filtering to identify chemical compounds that can simultaneously inhibit the human leukotriene hydrolase (hLTA4H) and the human leukotriene C4 synthase (hLTC4S) enzymes. These enzymes are the members of arachidonic acid pathway and act upon the same substrate, LTA4, producing different inflammatory products. A huge set of 4966 druglike compounds from the Maybridge database were docked into the active site of hLTA4H using the GOLD program. Common feature pharmacophore models were developed from the known inhibitors of both the targets using Accelrys Discovery Studio 2.5. The hits from the hLTA4H docking were filtered to match the chemical features of both the pharmacophore models. The compounds that resulted from the pharmacophore filtering were docked into the active site of hLTC4S and the hits those bind well at both the active sites and matched the pharmacophore models were identified as possible dual inhibitors for hLTA4H and hLTC4S enzymes. Reverse validation was performed to ensure the results of the study.     

Epitopes based drug design for dengue virus envelope protein: A computational approach

Dengue virus (DENV) has emerged as a rapidly spreading epidemic throughout the tropical and subtropical regions around the globe. No suitable drug has been designed yet to fight against DENV, therefore, the need for safe and effective antiviral drug has become imperative. The envelope protein of DENV is responsible for mediating the fusion process between viral and host membranes. This work reports an in silico approach to target B and T cell epitopes for dengue envelope protein inhibition. A conserved region “QHGTI” in B and T cell epitopes of dengue envelope glycoprotein was confirmed to be valid for targeting by visualizing its interactions with the host cell membrane TIM-1 protein which acts as a receptor for serotype 2 and 3. A reverse pharmacophore mapping approach was used to generate a seven featured pharmacophore model on the basis of predicted epitope. This pharmacophore model as a 3D query was used to virtually screen a chemical compounds dataset “Chembridge”. A total of 1010 compounds mapped on the developed pharmacophore model. These retrieved hits were subjected to filtering via Lipinski’s rule of five, as a result 442 molecules were shortlisted for further assessment using molecular docking. Finally, 14 hits of different structural properties having interactions with the active site residues of dengue envelope glycoprotein were selected as lead candidates. These structurally diverse lead candidates have strong likelihood to act as further starting structures in the development of novel and potential drugs for the treatment of dengue fever.     

Discovery and design of new PI3K inhibitors through pharmacophore-based virtual screening,molecular docking,and binding free energy analysis

Phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT)/the mammalian target of rapamycin (mTOR) signaling pathway is one of the sought after therapeutic target for treating human cancers. This pathway is often hyper activated in cancers. In the present study, pharmacophore-based virtual screening, molecular docking, and binding free energy calculations were performed on a series of quinoline derivatives which were reported to be effective against PI3Kα. A five-point pharmacophore hypothesis with one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic group (H), and two aromatic rings (R) was developed with acceptable R² and Q² values of 0.93 and 0.60 respectively. Eventually, common pharmacophore hypothesis-based screening was conducted against TOSLab, CPP, and ASINEX macrocylce databases, and potential hits were identified which were further subjected to rigorous docking process in order to screen out drug like molecules having crucial interactions with the target PI3Kα. Finally, binding free energy analysis was carried out for the top hits obtained from docking process. We also designed new 1, 3, 4-oxadiazole-based cyclic peptides by incorporating the structural features of the hits obtained from the above databases. Among the designed cyclic peptides, the cyclic peptide with tryptophan moiety showed good interactions and free binding energy values. On the whole, this study helped us in identifying new promising molecules as PI3Kα inhibitors which can be explored further to generate greater number of compounds with better pharmacokinetic properties.     

Development of energetic pharmacophore for the designing of 1,2,3,4-tetrahydropyrimidine derivatives as selective cyclooxygenase-2 inhibitors

We present here the Energetic pharmacophore model representing complementary features of the 1,2,3,4-tetrahydropyrimidine for selective cyclooxygenase-2 (COX-2) inhibition. For the development of pharmacophore hypothesis, a total of 43 previously reported compounds were docked on active site of COX-2 enzyme. The generated pharmacophore features were ranked using energetic terms of Glide XP docking for 1,2,3,4-tetrahydropyrimidine scaffold to optimize its structure requirement for COX-2 inhibition. The thirty new 4,5,6-triphenyl-1,2,3,4-tetrahydropyrimidine derivatives were synthesized and assessed for selective COX-2 inhibitory activity. Two compounds 4B1 and 4B11 were found to be potent and selective COX-2 inhibitors. The molecular docking studies revealed that the newly synthesized compounds can be docked into COX-2 binding site and also provide the molecular basis for their activity.     

Structure-based virtual screening of dipeptidyl peptidase 4 inhibitors and their in vitro analysis

Type 2 diabetes mellitus (T2DM) is one of the most widely prevalent metabolic disorders with no cure to date thus remains the most challenging task in the current drug discovery. Therefore, the only strategy to control diabetes prevalence is to develop novel efficacious therapeutics. Dipeptidyl Peptidase 4 (DPP-4) inhibitors are currently used as anti-diabetic drugs for the inhibition of incretins. This study aims to construct the chemical feature based on pharmacophore models for dipeptidyl peptidase IV. The structure-based pharmacophore modeling has been employed to evaluate new inhibitors of DPP-4. A four-featured pharmacophore model was developed from crystal structure of DPP-4 enzyme with 4-(2-aminoethyl) benzenesulfonyl fluoride in its active site via pharmacophore constructing tool of Molecular Operating Environment (MOE) consisting F1 Hyd (hydrophobic region), F2 Hyd|Cat|Don (hydrophobic cationic and donor region), F3 Acc (acceptor region) and F4 Hyd (hydrophobic region). The generated pharmacophore model was used for virtual screening of in-house compound library (the available compounds which were used for initial screening to get the few compounds for the current studies). The resultant selected compounds, after virtual screening were further validated using in vitro assay. Furthermore, structure-activity relationship was carried out for the compounds possessing significant inhibition potential after docking studies. The binding free energy of analogs was evaluated via molecular mechanics generalized Born surface area (MM-GBSA) and Poisson-Boltzmann surface area (MM-PBSA) methods using AMBER 16 as a molecular dynamics (MD) simulation package. Based on potential findings, we report that selected candidates are more likely to be used as DPP-4 inhibitors or as starting leads for the development of novel and potent DPP-4 inhibitors.     

Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.     

Identification of potential Tpx inhibitors against pathogen-host interactions

Yersinia organisms cause many infectious diseases by invading human cells and delivering their virulence factors via the type three secretion system (T3SS). One alternative strategy in the fight against these pathogenic organisms is to interfere with their T3SS. Previous studies demonstrated that thiol peroxidase, Tpx is functional in the assembly of T3SS and its inhibition by salicylidene acylhydrazides prevents the secretion of pathogenic effectors. In this study, the aim was to identify potential inhibitors of Tpx using an integrated approach starting with high throughput virtual screening and ending with molecular dynamics simulations of selected ligands. Virtual screening of ZINC database of 500,000 compounds via ligand-based and structure-based pharmacophore models retrieved 10,000 hits. The structure-based pharmacophore model was validated using high-throughput virtual screening (HTVS). After multistep docking (SP and XP), common scaffolds were used to find common substructures and the ligand binding poses were optimized using induced fit docking. The stability of the protein–ligand complex was examined with molecular dynamics simulations and the binding free energy of the complex was calculated. As a final outcome eight compounds with different chemotypes were proposed as potential inhibitors for Tpx. The eight ligands identified by a detailed virtual screening protocol can serve as leads in future drug design efforts against the destructive actions of pathogenic bacteria.     

Pharmacophore model-aided virtual screening combined with comparative molecular docking and molecular dynamics for identification of marine natural products as SARS-CoV-2 papain-like protease inhibitors

Targeting SARS-CoV-2 papain-like protease using inhibitors is a suitable approach for inhibition of virus replication and dysregulation of host anti-viral immunity. Engaging all five binding sites far from the catalytic site of PLpro is essential for developing a potent inhibitor. We developed and validated a structure-based pharmacophore model with 9 features of a potent PLpro inhibitor. The pharmacophore model-aided virtual screening of the comprehensive marine natural product database predicted 66 initial hits. This hit library was downsized by filtration through a molecular weight filter of ≤ 500 g/mol. The 50 resultant hits were screened by comparative molecular docking using AutoDock and AutoDock Vina. Comparative molecular docking enables benchmarking docking and relieves the disparities in the search and scoring functions of docking engines. Both docking engines retrieved 3 same compounds at different positions in the top 1 % rank, hence consensus scoring was applied, through which CMNPD28766, aspergillipeptide F emerged as the best PLpro inhibitor. Aspergillipeptide F topped the 50-hit library with a pharmacophore-fit score of 75.916. Favorable binding interactions were predicted between aspergillipeptide F and PLpro similar to the native ligand XR8-24. Aspergillipeptide F was able to engage all the 5 binding sites including the newly discovered BL2 groove, site V. Molecular dynamics for quantification of Cα-atom movements of PLpro after ligand binding indicated that it exhibits highly correlated domain movements contributing to the low free energy of binding and a stable conformation. Thus, aspergillipeptide F is a promising candidate for pharmaceutical and clinical development as a potent SARS-CoV-2 PLpro inhibitor.     

In Silico Strategy for Targeting the mTOR Kinase at Rapamycin Binding Site by Small Molecules

Computer aided drug-design methods proved to be powerful tools for the identification of new therapeutic agents. We employed a structure-based workflow to identify new inhibitors targeting mTOR kinase at rapamycin binding site. By combining molecular dynamics (MD) simulation and pharmacophore modelling, a simplified structure-based pharmacophore hypothesis was built starting from the FKBP12-rapamycin-FRB ternary complex retrieved from RCSB Protein Data Bank (PDB code 1FAP). Then, the obtained model was used as filter to screen the ZINC biogenic compounds library, containing molecules derived from natural sources or natural-inspired compounds. The resulting hits were clustered according to their similarity; moreover, compounds showing the highest pharmacophore fit-score were chosen from each cluster. The selected molecules were subjected to docking studies to clarify their putative binding mode. The binding free energy of the obtained complexes was calculated by MM/GBSA method and the hits characterized by the lowest ΔG_bind values were identified as potential mTOR inhibitors. Furthermore, the stability of the resulting complexes was studied by means of MD simulation which revealed that the selected compounds were able to form a stable ternary complex with FKBP12 and FRB domain, thus underlining their potential ability to inhibit mTOR with a rapamycin-like mechanism.     

On the value of homology models for virtual screening: discovering hCXCR3 antagonists by pharmacophore-based and structure-based approaches

Human chemokine receptor CXCR3 (hCXCR3) antagonists have potential therapeutic applications as antivirus, antitumor, and anti-inflammatory agents. A novel virtual screening protocol, which combines pharmacophore-based and structure-based approaches, was proposed. A three-dimensional QSAR pharmacophore model and a structure-based docking model were built to virtually screen for hCXCR3 antagonists. The hCXCR3 antagonist binding site was constructed by homology modeling and molecular dynamics (MD) simulation. By combining the structure-based and ligand-based screenings results, 95% of the compounds satisfied either pharmacophore or docking score criteria and would be chosen as hits if the union of the two searches was taken. The false negative rates were 15% for the pharmacophore model, 14% for the homology model, and 5% for the combined model. Therefore, the consistency of the pharmacophore model and the structural binding model is 219/273 = 80%. The hit rate for the virtual screening protocol is 273/286 = 95%. This work demonstrated that the quality of both the pharmacophore model and homology model can be measured by the consistency of the two models, and the false negatives in virtual screening can be reduced by combining two virtual screening approaches.     

Integrating pharmacophore mapping,virtual screening,density functional theory,molecular simulation towards the discovery of novel apolipoprotein (apoE ε4) inhibitors

AimAn integrated protocol of virtual screening involving molecular docking, pharmacophore probing, and simulations was established to identify small novel molecules targeting crucial residues involved in the variant apoE ε4 to mimic its behavior as apoE2 thereby eliminating the amyloid plaque accumulation and facilitating its clearance.Materials and MethodsAn excellent ligand-based and structure-based approach was made to identify common pharmacophoric features involving structure-based docking with respect to apoE ε4 leading to the development of apoE ε4 inhibitors possessing new scaffolds. An effort was made to design multiple-substituted triazine derivatives series bearing a novel scaffold. A structure-based pharmacophore mapping was developed to explore the binding sites of apoE ε4 which was taken into consideration. Subsequently, virtual screening, ADMET, DFT searches were at work to narrow down the proposed hits to be forwarded as a potential drug likes candidates. Further, the binding patterns of the best-proposed hits were studied and were forwarded for molecular dynamic simulations of 10 ns for its structural optimization.ResultsSelectivity profile for the most promising candidates was studied, revealing significantly C13 and C15 to be the most potent compounds. The proposed hits can be forwarded for further study against apoE ε4 involved in neurological disorder Alzheimer’s.     

Evaluation of potential inhibitors of squalene synthase based on virtual screening and in vitro studies

Squalene synthase (SQS) is a potential target for hyperlipidemia treatment. To identify novel chemical scaffolds of SQS inhibitors, we generated 3D-QSAR pharmacophore models using HypoGen. The best quantitative pharmacophore model, Hypo 1, was selected for virtual screening using two chemical databases, Specs and Traditional Chinese Medicine database (TCM). The best-mapped hit compounds were then subjected to filtering by Lipinski^’s rule of five and docking studies to refine the hits. Finally, five compounds were selected from the top-ranked hit compounds for SQS inhibitory assay in vitro. Three of these compounds could inhibit SQS in vitro, and should be further evaluated pre-clinically as a treatment for hyperlipidemia.     

Molecular dynamics simulation study of PTP1B with allosteric inhibitor and its application in receptor based pharmacophore modeling

Allosteric inhibition of protein tyrosine phosphatase 1B (PTP1B), has paved a new path to design specific inhibitors for PTP1B, which is an important drug target for the treatment of type II diabetes and obesity. The PTP1B_1–282-allosteric inhibitor complex crystal structure lacks α7 (287–298) and moreover there is no available 3D structure of PTP1B_1–298 in open form. As the interaction between α7 and α6–α3 helices plays a crucial role in allosteric inhibition, α7 was modeled to the PTP1B_1–282 in open form complexed with an allosteric inhibitor (compound-2) and a 5 ns MD simulation was performed to investigate the relative orientation of the α7–α6–α3 helices. The simulation conformational space was statistically sampled by clustering analyses. This approach was helpful to reveal certain clues on PTP1B allosteric inhibition. The simulation was also utilized in the generation of receptor based pharmacophore models to include the conformational flexibility of the protein-inhibitor complex. Three cluster representative structures of the highly populated clusters were selected for pharmacophore model generation. The three pharmacophore models were subsequently utilized for screening databases to retrieve molecules containing the features that complement the allosteric site. The retrieved hits were filtered based on certain drug-like properties and molecular docking simulations were performed in two different conformations of protein. Thus, performing MD simulation with α7 to investigate the changes at the allosteric site, then developing receptor based pharmacophore models and finally docking the retrieved hits into two distinct conformations will be a reliable methodology in identifying PTP1B allosteric inhibitors. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Pharmacophore generation,atom-based 3D-QSAR and molecular dynamics simulation analyses of pyridine-3-carboxamide-6-yl-urea analogues as potential gyrase B inhibitors

DNA gyrase subunit B (GyrB) is an attractive drug target for the development of antibacterial agents with therapeutic potential. In the present study, computational studies based on pharmacophore modelling, atom-based QSAR, molecular docking, free binding energy calculation and dynamics simulation were performed on a series of pyridine-3-carboxamide-6-yl-urea derivatives. A pharmacophore model using 49 molecules revealed structural and chemical features necessary for these molecules to inhibit GyrB. The best fitted model AADDR.13 was generated with a coefficient of determination (r²) of 0.918. This model was validated using test set molecules and had a good r² of 0.78. 3D contour maps generated by the 3D atom-based QSAR revealed the key structural features responsible for the GyrB inhibitory activity. Extra precision molecular docking showed hydrogen bond interactions with key amino acid residues of ATP-binding pocket, important for inhibitor binding. Further, binding free energy was calculated by the MM-GBSA rescoring approach to validate the binding affinity. A 10 ns MD simulation of inhibitor #47 showed the stability of the predicted binding conformations. We identified 10 virtual hits by in silico high-throughput screening. A few new molecules were also designed as potent GyrB inhibitors. The information obtained from these methodologies may be helpful to design novel inhibitors of GyrB.     

Combination of pharmacophore modeling and 3D-QSAR analysis of potential glyoxalase-I inhibitors as anticancer agents

Glyoxalase system is an ubiquitous system in human cells which has been examined thoroughly for its role in different diseases. It comprises two enzymes; Glyoxalase I (Glo-I) and Glyoxalase II (Glo-II) which perform detoxifying endogenous harmful metabolites, mainly methylglyoxal (MG) into non-toxic bystanders. In silico computer Aided Drug Design approaches were used and ninety two diverse pharmacophore models were generated from eighteen Glyoxalase I crystallographic complexes. Subsequent QSAR modeling followed by ROC evaluation identified a single pharmacophore model which was able to predict the expected Glyoxalase I inhibition. Screening of the National Cancer Institute (NCI) database using the optimal pharmacophore Hypo(3VW9) identified several promising hits. Thirty eight hits were successfully predicted then ordered and evaluated in vitro. Seven hits out of the thirty eight tested compounds showed more than 50% inhibition with low micromolar IC₅₀.     

E-pharmacophore-based screening of mGluR5 negative allosteric modulators for central nervous system disorder

Glutamate, a major neurotransmitter in the central nervous system of human, plays a crucial role in various neurological pathways by activating the ligand-gated ion channels such as mGluR and iGluR. Dysfunction of mGluR 5 can cause Alzheimer’s disease, Parkinson’s disease, epilepsy, depression, anxiety, etc. In the current study, we have developed the energetically optimized pharmacophore model to screen the eMolecules database having more than 6 million compounds with the help of reported cocrystal structure with 3-chloro-5-[6-(5-fluoropyridin-2 yl)pyrimidin-4-yl]benzonitrile (PDB ID: 5CGD). The obtained hits were docked into the allosteric site of the target and further validation of E-pharmacophore was done by enrichment calculations followed by the molecular dynamics simulations to analyze the specific amino acid interactions with the compound present in the allosteric site of the receptor.     

Pharmacophore Modeling and Virtual Screening of Novel Inhibitors for c-Kit Kinase

The stem cell factor receptor (c‐Kit) has been known to play critical roles in regulating numerous aspects of cellular behavior including cell growth, differentiation, migration and metabolism. In this investigation, a three‐dimensional pharmacophore model of c‐Kit inhibitors has been established by using the HypoGen algorithms implemented in the catalyst software package. The best quantitative pharmacophore model, hypothesis 1, which has the highest correlation coefficient (0.989), consists of one hydrogen bond acceptor, two hydrogen bond donors and one hydrophobic feature. To our knowledge, this is the first report on the pharmacophore modeling study of c‐Kit inhibitors. The best hypothesis, hypothesis 1, was used to screen molecular structural databases, including Specs and China Natural Products Database for potential lead compounds. The hit compounds were subsequently subjected to filtering by Lipinski's rules and docking study to refine the retrieved hits and as a result to reduce the rate of false positive. Finally 28 compounds were purchased or synthesized for further in vitro assay against several human tumour cell lines including A549, MCF‐7, HepG2 and PC‐3, in which c‐Kit is overexpressed. Two compounds show very low micromolar inhibition potency against the PC‐3 and HepG2 cell lines respectively. And they were selected for further modification and testing.     

Identification of novel PfDHODH inhibitors as antimalarial agents via pharmacophore-based virtual screening followed by molecular docking and in vivo antimalarial activity

Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) catalyses the fourth reaction of de novo pyrimidine biosynthesis in parasites, and represents an important target for the treatment of malaria. In this study, we describe pharmacophore-based virtual screening combined with docking study and biological evaluation as a rational strategy for identification of novel hits as antimalarial agents. Pharmacophore models were established from known PfDHODH inhibitors using the GALAHAD module with IC₅₀ values ranging from 0.033 μM to 142 μM. The best pharmacophore model consisted of three hydrogen bond acceptor, one hydrogen bond donor and one hydrophobic features. The pharmacophore models were validated through receiver operating characteristic and Günere–Henry scoring methods. The best pharmacophore model as a 3D search query was searched against the IBS database. Several compounds with different structures (scaffolds) were retrieved as hit molecules. Among these compounds, those with a Q_FIT value of more than 81 were docked in the PfDHODH enzyme to further explore the binding modes of these compounds. In silico pharmacokinetic and toxicities were predicted for the best docked molecules. Finally, the identified hits were evaluated in vivo for their antimalarial activity in a parasite inhibition assay. The hits reported here showed good potential to become novel antimalarial agents.     

Receptor-Based Pharmacophore Modeling in the Search for Natural Products for COVID-19 Mpro

Considering the urgency of the COVID-19 pandemic, we developed a receptor-based pharmacophore model for identifying FDA-approved drugs and hits from natural products. The COVID-19 main protease (M^pro) was selected for the development of the pharmacophore model. The model consisted of a hydrogen bond acceptor, donor, and hydrophobic features. These features demonstrated good corroboration with a previously reported model that was used to validate the present model, showing an RMSD value of 0.32. The virtual screening was carried out using the ZINC database. A set of 208,000 hits was extracted and filtered using the ligand pharmacophore mapping, applying the lead-like properties. Lipinski’s filter and the fit value filter were used to minimize hits to the top 2000. Simultaneous docking was carried out for 200 hits for natural drugs belonging to the FDA-approved drug database. The top 28 hits from these experiments, with promising predicted pharmacodynamic and pharmacokinetic properties, are reported here. To optimize these hits as M^pro inhibitors and potential treatment options for COVID-19, bench work investigations are needed.     

Surrogate docking: structure-based virtual screening at high throughput speed

Summary Structure-based screening using fully flexible docking is still too slow for large molecular libraries. High quality docking of a million molecule library can take days even on a cluster with hundreds of CPUs. This performance issue prohibits the use of fully flexible docking in the design of large combinatorial libraries. We have developed a fast structure-based screening method, which utilizes docking of a limited number of compounds to build a 2D QSAR model used to rapidly score the rest of the database. We compare here a model based on radial basis functions and a Bayesian categorization model. The number of compounds that need to be actually docked depends on the number of docking hits found. In our case studies reasonable quality models are built after docking of the number of molecules containing 50 docking hits. The rest of the library is screened by the QSAR model. Optionally a fraction of the QSAR-prioritized library can be docked in order to find the true docking hits. The quality of the model only depends on the training set size – not on the size of the library to be screened. Therefore, for larger libraries the method yields higher gain in speed no change in performance. Prioritizing a large library with these models provides a significant enrichment with docking hits: it attains the values of 13 and 35 at the beginning of the score-sorted libraries in our two case studies: screening of the NCI collection and a combinatorial libraries on CDK2 kinase structure. With such enrichments, only a fraction of the database must actually be docked to find many of the true hits. The throughput of the method allows its use in screening of large compound collections and in the design of large combinatorial libraries. The strategy proposed has an important effect on efficiency but does not affect retrieval of actives, the latter being determined by the quality of the docking method itself. Electronic supplementary material is available at http://dx.doi.org/10.1007/s10822-005-9002-6.