A pseudo-ligand approach to virtual screening

A virtual screening method is presented that is grounded on a receptor-derived pharmacophore model termed "virtual ligand" or "pseudo-ligand". The model represents an idealized constellation of potential ligand sites that interact with residues of the binding pocket. For rapid virtual screening of compound libraries the potential pharmacophore points of the virtual ligand are encoded as an alignment-free correlation vector, avoiding spatial alignment of pharmacophore features between the pharmacophore query (i.e., the virtual ligand) and the candidate molecule. The method was successfully applied to retrieving factor Xa inhibitors from a Ugi three-component combinatorial library, and yielded high enrichment of actives in a retrospective search for cyclooxygenase-2 (COX-2) inhibitors. The approach provides a concept for "de-orphanizing" potential drug targets and identifying ligands for hitherto unexplored or allosteric binding pockets.     

Selection of DNA-Encoded Dynamic Chemical Libraries for Direct Inhibitor Discovery

Dynamic combinatorial libraries (DCLs) is a powerful tool for ligand discovery in biomedical research; however, the application of DCLs has been hampered by their low diversity. Recently, the concept of DNA encoding has been employed in DCLs to create DNA-encoded dynamic libraries (DEDLs); however, all current DEDLs are limited to fragment identification, and a challenging process of fragment linking is required after selection. We report an anchor-directed DEDL approach that can identify full ligand structures from large-scale DEDLs. This method is also able to convert unbiased libraries into focused ones targeting specific protein classes. We demonstrated this method by selecting DEDLs against five proteins, and novel inhibitors were identified for all targets. Notably, several selective BD1/BD2 inhibitors were identified from the selections against bromodomain 4 (BRD4), an important anti-cancer drug target. This work may provide a broadly applicable method for inhibitor discovery.     

Selection of DNA‐Encoded Dynamic Chemical Libraries for Direct Inhibitor Discovery

Dynamic combinatorial libraries (DCLs) is a powerful tool for ligand discovery in biomedical research; however, the application of DCLs has been hampered by their low diversity. Recently, the concept of DNA encoding has been employed in DCLs to create DNA‐encoded dynamic libraries (DEDLs); however, all current DEDLs are limited to fragment identification, and a challenging process of fragment linking is required after selection. We report an anchor‐directed DEDL approach that can identify full ligand structures from large‐scale DEDLs. This method is also able to convert unbiased libraries into focused ones targeting specific protein classes. We demonstrated this method by selecting DEDLs against five proteins, and novel inhibitors were identified for all targets. Notably, several selective BD1/BD2 inhibitors were identified from the selections against bromodomain 4 (BRD4), an important anti‐cancer drug target. This work may provide a broadly applicable method for inhibitor discovery.     

Dynamic Amino Acid Side-Chains Grafting on Folded Peptide Backbone**

An efficient strategy for the synthesis of large libraries of conformationally defined peptides is reported, using dynamic combinatorial chemistry as a tool to graft amino acid side chains on a well-ordered 3D (3-dimension) peptide backbone. Combining rationally designed scaffolds with combinatorial side chains selection represents an alternative method to access peptide libraries for structures that are not genetically encodable. This method would allow a breakthrough for the discovery of protein mimetic for unconventional targets for which little is known.     

Virtual screening using protein-ligand docking: avoiding artificial enrichment

This study addresses a number of topical issues around the use of protein-ligand docking in virtual screening. We show that, for the validation of such methods, it is key to use focused libraries (containing compounds with one-dimensional properties, similar to the actives), rather than "random" or "drug-like" libraries to test the actives against. We also show that, to obtain good enrichments, the docking program needs to produce reliable binding modes. We demonstrate how pharmacophores can be used to guide the dockings and improve enrichments, and we compare the performance of three consensus-ranking protocols against ranking based on individual scoring functions. Finally, we show that protein-ligand docking can be an effective aid in the screening for weak, fragment-like binders, which has rapidly become a popular strategy for hit identification. All results presented are based on carefully constructed virtual screening experiments against four targets, using the protein-ligand docking program GOLD.     

Evaluating virtual screening methods: good and bad metrics for the "early recognition" problem

Many metrics are currently used to evaluate the performance of ranking methods in virtual screening (VS), for instance, the area under the receiver operating characteristic curve (ROC), the area under the accumulation curve (AUAC), the average rank of actives, the enrichment factor (EF), and the robust initial enhancement (RIE) proposed by Sheridan et al. In this work, we show that the ROC, the AUAC, and the average rank metrics have the same inappropriate behaviors that make them poor metrics for comparing VS methods whose purpose is to rank actives early in an ordered list (the "early recognition problem"). In doing so, we derive mathematical formulas that relate those metrics together. Moreover, we show that the EF metric is not sensitive to ranking performance before and after the cutoff. Instead, we formally generalize the ROC metric to the early recognition problem which leads us to propose a novel metric called the Boltzmann-enhanced discrimination of receiver operating characteristic that turns out to contain the discrimination power of the RIE metric but incorporates the statistical significance from ROC and its well-behaved boundaries. Finally, two major sources of errors, namely, the statistical error and the "saturation effects", are examined. This leads to practical recommendations for the number of actives, the number of inactives, and the "early recognition" importance parameter that one should use when comparing ranking methods. Although this work is applied specifically to VS, it is general and can be used to analyze any method that needs to segregate actives toward the front of a rank-ordered list.     

Enhanced encapsulation of actives in self-sealing microcapsules by precipitation in capsule shells

Microcapsules with core-shell structures are excellent vehicles for the encapsulation of active ingredients; however, the actives often leak out of these structures over time, without observable damage to them. We present a novel approach to enhancing the encapsulation of active ingredients inside microcapsules. We use two components that can form solid precipitates upon mixing and add one each to the microcapsule core and to the continuous phase. The components diffuse through the shell in the same manner as the actives, but upon meeting, they precipitate to form solid particles within the shell; this significantly reduces leakage through the shell of the microcapsules. We show that the reduction in the leakage of actives is due to the blockage of channels or pores that exist in the shell of the capsules by the solid precipitates.     

A general method for designing combinatorial peptide libraries decodable by amino acid analysis

Herein we describe an algorithm for designing combinatorial peptide libraries for split-and-mix synthesis on solid support that are decodable by amino acid analysis (AAA) of the beads. AAA is a standard service analysis available in most biochemical laboratories, and it allows one to control the quality of the peptide on each bead, an important feature that is missing from most library decoding protocols. In the algorithm, each AA is assigned to two variable positions in the sequence grouped in a "unique pair". This arrangement limits sequence design because both the number of unique pairs U (setting the maximum number of variable AA) and the maximum number S of different AA per variable position depend on the peptide length N (U=N(N-1)/2), S=N-1). The method is therefore only suitable for focused libraries. An application example is shown for the selection of peptides with N-terminal proline or hydroxyproline catalyzing an aldol reaction from a combinatorial library of 65536 octapeptides. A simple enumeration program is available to help design combinatorial libraries decodable by amino acid analysis. The method applies to linear and cyclic peptides, can be used for nonnatural building blocks, including beta-amino acids, and should help to explore the vast chemistry of linear and cyclic peptide for catalysis and bioactivity.     

Retrospective analysis of an experimental high-throughput screening data set by recursive partitioning

With the emergence of combinatorial chemistry, whether based on parallel, mixture, solution, or solid phase chemistry, it is now possible to generate large numbers of diverse or focused compound libraries. In this paper we aim to demonstrate that it is possible to design targeted libraries by applying nonparametric statistical methods, recursive partitioning in particular, to large data sets containing thousands of compounds and their associated biological data. Moreover, when applied to an experimental high-throughput screening (HTS) data set, our data strongly suggest that this method can improve the hit rate of our primary screens (about 4- to 5-fold) while increasing screening efficiency: less than one-fifth of the complete selection needs to be screened in order to identify about 75% of all actives present.     

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.     

Ligand design by a combinatorial approach based on modeling and experiment: application to HLA-DR4

Combinatorial synthesis and large scale screening methods are being used increasingly in drug discovery, particularly for finding novel lead compounds. Although these "random" methods sample larger areas of chemical space than traditional synthetic approaches, only a relatively small percentage of all possible compounds are practically accessible. It is therefore helpful to select regions of chemical space that have greater likelihood of yielding useful leads. When three-dimensional structural data are available for the target molecule this can be achieved by applying structure-based computational design methods to focus the combinatorial library. This is advantageous over the standard usage of computational methods to design a small number of specific novel ligands, because here computation is employed as part of the combinatorial design process and so is required only to determine a propensity for binding of certain chemical moieties in regions of the target molecule. This paper describes the application of the Multiple Copy Simultaneous Search (MCSS) method, an active site mapping and de novo structure-based design tool, to design a focused combinatorial library for the class II MHC protein HLA-DR4. Methods for the synthesizing and screening the computationally designed library are presented; evidence is provided to show that binding was achieved. Although the structure of the protein-ligand complex could not be determined, experimental results including cross-exclusion of a known HLA-DR4 peptide ligand (HA) by a compound from the library. Computational model building suggest that at least one of the ligands designed and identified by the methods described binds in a mode similar to that of native peptides.     

Rapid shape-based ligand alignment and virtual screening method based on atom/feature-pair similarities and volume overlap scoring

Shape-based methods for aligning and scoring ligands have proven to be valuable in the field of computer-aided drug design. Here, we describe a new shape-based flexible ligand superposition and virtual screening method, Phase Shape, which is shown to rapidly produce accurate 3D ligand alignments and efficiently enrich actives in virtual screening. We describe the methodology, which is based on the principle of atom distribution triplets to rapidly define trial alignments, followed by refinement of top alignments to maximize the volume overlap. The method can be run in a shape-only mode or it can include atom types or pharmacophore feature encoding, the latter consistently producing the best results for database screening. We apply Phase Shape to flexibly align molecules that bind to the same target and show that the method consistently produces correct alignments when compared with crystal structures. We then illustrate the effectiveness of the method for identifying active compounds in virtual screening of eleven diverse targets. Multiple parameters are explored, including atom typing, query structure conformation, and the database conformer generation protocol. We show that Phase Shape performs well in database screening calculations when compared with other shape-based methods using a common set of actives and decoys from the literature.     

Combining combinatorial chemistry and affinity chromatography: highly selective inhibitors of human betaine: homocysteine S-methyltransferase

A new method to find novel protein targets for ligands of interest is proposed. The principle of this approach is based on affinity chromatography and combinatorial chemistry. The proteins within a crude rat liver homogenate were allowed to interact with a combinatorial library of phosphinic pseudopeptides immobilized on affinity columns. Betaine: homocysteine S-methyltransferase (BHMT) was one of the proteins that was retained and subsequently eluted from these supports. The phosphinic pseudopeptides, which served as immobilized ligands for the isolation of rat BHMT, were then tested for their ability to inhibit human recombinant BHMT in solution. The most potent inhibitor also behaved as a selective ligand for the affinity purification of BHMT from a complex media. Further optimization uncovered Val-Phe-psi[PO(2-)-CH(2)]-Leu-His-NH(2) as a potent BHMT inhibitor that has an IC(50) of about 1 microM.     

Comprehensive comparison of ligand-based virtual screening tools against the DUD data set reveals limitations of current 3D methods

In recent years, many virtual screening (VS) tools have been developed that employ different molecular representations and have different speed and accuracy characteristics. In this paper, we compare ten popular ligand-based VS tools using the publicly available Directory of Useful Decoys (DUD) data set comprising over 100?000 compounds distributed across 40 protein targets. The DUD was developed initially to evaluate docking algorithms, but our results from an operational correlation analysis show that it is also well suited for comparing ligand-based VS tools. Although it is conventional wisdom that 3D molecular shape is an important determinant of biological activity, our results based on permutational significance tests of several commonly used VS metrics show that the 2D fingerprint-based methods generally give better VS performance than the 3D shape-based approaches for surprisingly many of the DUD targets. To help understand this finding, we have analyzed the nature of the scoring functions used and the composition of the DUD data set itself. We propose that to improve the VS performance of current 3D methods, it will be necessary to devise screening queries that can represent multiple possible conformations and which can exploit knowledge of known actives that span multiple scaffold families.     

Multiobjective optimization of combinatorial libraries

Combinatorial chemistry and high-throughput screening have caused a fundamental shift in the way chemists contemplate experiments. Designing a combinatorial library is a controversial art that involves a heterogeneous mix of chemistry, mathematics, economics, experience, and intuition. Although there seems to be little agreement as to what constitutes an ideal library, one thing is certain: only one property or measure seldom defines the quality of the design. In most real-world applications, a good experiment requires the simultaneous optimization of several, often conflicting, design objectives, some of which may be vague and uncertain. In this paper, we discuss a class of algorithms for subset selection rooted in the principles of multiobjective optimization. Our approach is to employ an objective function that encodes all of the desired selection criteria, and then use a simulated annealing or evolutionary approach to identify the optimal (or a nearly optimal) subset from among the vast number of possibilities. Many design criteria can be accommodated, including diversity, similarity to known actives, predicted activity and/or selectivity determined by quantitative structure-activity relationship (QSAR) models or receptor binding models, enforcement of certain property distributions, reagent cost and availability, and many others. The method is robust, convergent, and extensible, offers the user full control over the relative significance of the various objectives in the final design, and permits the simultaneous selection of compounds from multiple libraries in full- or sparse-array format.     

Achieving High Affinity and Selectivity for Asymmetric Dimethylarginine by Putting a Lid on a Box

The methylation states of Lys and Arg represent a particularly challenging set of targets to distinguish selectively in water using synthetic receptors. To date, trimethyllysine (Kme3) is the only post translational modification (PTM) of the eight possible methylation states of Lys and Arg that can be recognized selectively. Here, we report the first synthetic receptor capable of selectively recognizing asymmetric dimethylarginine (Rme2a). This was achieved by using a biased dynamic combinatorial chemistry (DCC) library to generate a receptor mimicking the 5‐sided box‐like shape of Rme2 reader proteins, a feature that has been hypothesized to impart selectivity. Additionally, we synthesized a thioether‐linked analogue of the resulting receptor to provide a novel scaffold with maintained selectivity but greater stability. This work introduces strategies that can be applied towards achieving selectivity based on subtle differences in hydrophilic guests in aqueous solutions.