OptiDock: Virtual HTS of Combinatorial Libraries by Efficient Sampling of Binding Modes in Product Space

Methods for combinatorial and parallel library design

Diversity has historically played a critical role in design of combinatorial libraries, screening sets and corporate collections for lead discovery. Large library design dominated the field in the 1990s with methods ranging anywhere from purely arbitrary through property based reagent selection to product based approaches. In recent years, however, there has been a downward trend in library size. This was due to increased information about the desirable targets gleaned from the genomics revolution and to the ever growing availability of target protein structures from crystallography and homology modeling. Creation of libraries directed toward families of receptors such as GPCRs, kinases, nuclear hormone receptors, proteases, etc., replaced the generation of libraries based primarily on diversity while single target focused library design has remained an important objective. Concurrently, computing grids and cpu clusters have facilitated the development of structure based tools that screen hundreds of thousands of molecules. Smaller "smarter" combinatorial and focused parallel libraries replaced those early un-focused large libraries in the twenty-first century drug design paradigm. While diversity still plays a role in lead discovery, the focus of current library design methods has shifted to receptor based methods, scaffold hopping/bio-isostere searching, and a much needed emphasis on synthetic feasibility. Methods such as "privileged substructures based design" and pharmacophore based design still are important methods for parallel and small combinatorial library design. This chapter discusses some of the possible design methods and presents examples where they are available.

Docking methods for structure-based library design

The drug discovery process mainly relies on the experimental high-throughput screening of huge compound libraries in their pursuit of new active compounds. However, spiraling research and development costs and unimpressive success rates have driven the development of more rational, efficient, and cost-effective methods. With the increasing availability of protein structural information, advancement in computational algorithms, and faster computing resources, in silico docking-based methods are increasingly used to design smaller and focused compound libraries in order to reduce screening efforts and costs and at the same time identify active compounds with a better chance of progressing through the optimization stages. This chapter is a primer on the various docking-based methods developed for the purpose of structure-based library design. Our aim is to elucidate some basic terms related to the docking technique and explain the methodology behind several docking-based library design methods. This chapter also aims to guide the novice computational practitioner by laying out the general steps involved for such an exercise. Selected successful case studies conclude this chapter.

Structure-directed combinatorial library design

In recent years pharmaceutical companies have utilized structure-based drug design and combinatorial library design techniques to speed up their drug discovery efforts. Both approaches are routinely used in the lead discovery and lead optimization stages of the drug discovery process. Fragment-based drug design, a new power tool in the drug design toolbox, is also gaining acceptance across the pharmaceutical industry. This review will focus on the interplay between these three design techniques and recent developments in computational methodologies that enhance their integration. Examples of successful synergistic applications of these three techniques will be highlighted. Opinion regarding possible future directions of the field will be given.

eHiTS-to-VMD interface application. The search for tyrosine-tRNA ligase inhibitors

Owing to the recent development of virtual high-throughput screening (vHTS) and a vast number of compounds subjected to vHTS analyses, it has been essential to automate the processing of computational data required for the analysis and visualization of research results. Using the search for tyrosine-tRNA ligase inhibitors as an example, we present a computer application, an interface between eHiTS software for virtual high-throughput screening and VMD graphic software used to visualize calculation results.

Balancing focused combinatorial libraries based on multiple GPCR ligands

G-Protein coupled receptors (GPCRs) are important targets for drug discovery, and combinatorial chemistry is an important tool for pharmaceutical development. The absence of detailed structural information, however, limits the kinds of combinatorial design techniques that can be applied to GPCR targets. This is particularly problematic given the current emphasis on focused combinatorial libraries. By linking an incremental construction method (OptDesign) to the very fast shape-matching capability of ChemSpace, we have created an efficient method for designing targeted sublibraries that are topomerically similar to known actives. Multi-objective scoring allows consideration of multiple queries (actives) simultaneously. This can lead to a distribution of products skewed towards one particular query structure, however, particularly when the ligands of interest are quite dissimilar to one another. A novel pivoting technique is described which makes it possible to generate promising designs even under those circumstances. The approach is illustrated by application to some serotonergic agonists and chemokine antagonists.

4SCan/vADME: intelligent library screening as a shortcut from hits to lead compounds

Managing to solve the first step in drug discovery - the hit finding - can be a quite elaborate task, but it is only the initial step to the final goal; hit-to-lead optimisation still lies ahead and consumes even more time and resources. The solution is rather simple, that is, to take only the most promising compounds into account; but who is going to decide which ones are the most promising among a list of tens of millions of compounds in a virtual combinatorial library? 4SCan/vADME helps by bridging the gap between virtual (combinatorial) libraries designed by chemists and the in silico methods, docking and alignment, for screening databases. After choosing a random starting set, the implemented learning and prediction algorithm iteratively considers only combinations of fragments that have shown to result in more suitable interactions by the chosen method. ADME properties of the final list are then calculated via several in silico methods, resulting in a combined evaluation of the individual compound's target-specific, as well as ADME, properties. Based on the latter list of evaluated compounds, medicinal chemists can then decide which compounds might be the best ones to synthesise first and to serve as possible lead candidates. Following a brief introduction to virtual high-throughput screening techniques, the 4SCan/vADME method is outlined and discussed in this paper, using an example coming out of the 4SC pipeline.

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

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 approximately 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 approximately 13 and approximately 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.