Combinatorial library design using a multiobjective genetic algorithm

Early results from screening combinatorial libraries have been disappointing with libraries either failing to deliver the improved hit rates that were expected or resulting in hits with characteristics that make them undesirable as lead compounds. Consequently, the focus in library design has shifted toward designing libraries that are optimized on multiple properties simultaneously, for example, diversity and "druglike" physicochemical properties. Here we describe the program MoSELECT that is based on a multiobjective genetic algorithm and which is able to suggest a family of solutions to multiobjective library design where all the solutions are equally valid and each represents a different compromise between the objectives. MoSELECT also allows the relationships between the different objectives to be explored with competing objectives easily identified. The library designer can then make an informed choice on which solution(s) to explore. Various performance characteristics of MoSELECT are reported based on a number of different combinatorial libraries.     

Use of catalyst pharmacophore models for screening of large combinatorial libraries

Using a data set comprised of literature compounds and structure-activity data for cyclin dependent kinase 2, several pharmacophore hypotheses were generated using Catalyst and evaluated using several criteria. The two best were used in retrospective searches of 10 three-dimensional databases containing over 1,000,000 proprietary compounds. The results were then analyzed for the efficiency with which the hypotheses performed in the areas of compound prioritization, library prioritization, and library design. First as a test of their compound prioritization capabilities, the pharmacophore models were used to search combinatorial libraries that were known to contain CDK active compounds to see if the pharmacophore models could selectively choose the active compounds over the inactive compounds. Second as a test of their utility in library design again the pharmacophore models were used to search the active combinatorial libraries to see if the key synthons were over represented in the hits from the pharmacophore searches. Finally as a test of their ability to prioritize combinatorial libraries, several inactive libraries were searched in addition to the active libraries in order to see if the active libraries produced significantly more hits than the inactive libraries. For this study the pharmacophore models showed potential in all three areas. For compound prioritization, one of the models selected active compounds at a rate nearly 11 times that of random compound selection though in other cases models missed the active compounds entirely. For library design, most of the key fragments were over represented in the hits from at least one of the searches though again some key fragments were missed. Finally, for library prioritization, the two active libraries both produced a significant number of hits with both pharmacophore models, whereas none of the eight inactive libraries produced a significant number of hits for both models.     

High Speed Development and Synthesis of Novel Small Molecule Libraries

Combinatorial chemistry has produced libraries of millions of compounds in the last decade. Screening of those compounds, unfortunately, has not yet yielded as many new drug candidates as initially expected. Among a number of possible reasons, one is that many libraries combinatorial chemistry produced in the early periods are of the nature of linear, flat, and flexible molecules such as peptides and oligonucleotides, which do not have the desired properties to selectively interact with their targets to yield high quality hits and leads. In order to increase the number of quality hits and leads, rigid, structural featurerich and drug-like compound libraries are highly desirable. Design and development of structural features-rich and natural product-like combinatorial libraries, as well as high speed library production using modern solution and solid phase synthesis techniques such as IRORI's Directed Sorting technology, will be discussed.     

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.     

Simulated molecular evolution in a full combinatorial library

BACKGROUND: The Darwinian concept of 'survival of the fittest' has inspired the development of evolutionary optimization methods to find molecules with desired properties in iterative feedback cycles of synthesis and testing. These methods have recently been applied to the computer-guided heuristic selection of molecules that bind with high affinity to a given biological target. We describe the optimization behavior and performance of genetic algorithms (GAs) that select molecules from a combinatorial library of potential thrombin inhibitors in 'artificial molecular evolution' experiments, on the basis of biological screening results. RESULTS: A full combinatorial library of 15,360 members structurally biased towards the serine protease thrombin was synthesized, and all were tested for their ability to inhibit the protease activity of thrombin. Using the resulting large structure-activity landscape, we simulated the evolutionary selection of potent thrombin inhibitors from this library using GAs. Optimal parameter sets were found (encoding strategy, population size, mutation and cross-over rate) for this artificial molecular evolution. CONCLUSIONS: A GA-based evolutionary selection is a valuable combinatorial optimization strategy to discover compounds with desired properties without needing to synthesize and test all possible combinations (i.e. all molecules). GAs are especially powerful when dealing with very large combinatorial libraries for which synthesis and screening of all members is not possible and/or when only a small number of compounds compared with the library size can be synthesized or tested. The optimization gradient or 'learning' per individual increases when using smaller population sizes and decreases for higher mutation rates.     

Constructing virtual combinatorial fragment libraries based upon MDL Drug Data Report database

Structural analysis of known drugs or drug-like compounds provides important information for drug design. The 142553 drug molecules in the MDL Drug Data Report database were analyzed, and then the common structural features were extracted. According to the common structural features, drug molecules were segmented into 32017 fragments, including 13642 ring fragments, 10076 linker fragments, and 8299 side chain fragments. These fragments were further used to establish three types of virtual combinatorial fragment libraries: a basic framework library containing 13574 rings; a linker library of 8051 linkers and a pharmacophore library of 34244 fragments combined by rings and side chains. After energy minimization, all fragments in the above three libraries maintain reasonable geometrical features and spatial conformations, and would be useful for building a virtual combinatorial database and de novo drug design.     

Polymer-supported synthesis of pyridone-focused libraries as inhibitors of anaplastic lymphoma kinase

The optimization of screening hits on a promising new target for therapy of certain cancers involving anaplastic lymphoma kinase (ALK) inspired the development of this efficient solid-phase chemistry. A series of novel pyridones have been recently discovered as inhibitors of ALK, which led to the design of focused libraries around the pyridone scaffold. A stepwise process involving iterative template modification based on both medicinal chemistry insights and computational ranking of virtual libraries was employed in the design. The unique solid-phase chemistry has addressed the need for rapid optimization of this "early lead" series. Herein the methodology and scope of the chemistry, as well as its application for library synthesis, are discussed.     

Where are the GaPs? A rational approach to monomer acquisition and selection

Gridding and partitioning (GaP) is a computational method for the classification and selection of monomers for combinatorial libraries. The molecules are described in terms of the pharmacophoric groups they contain and where those pharmacophoric groups can be located in three-dimensional space. The approach involves a detailed conformational analysis of each molecule. This conformational analysis is done within a common coordinate frame, thus enabling the monomers to be compared. The use of a partitioned space is central to this particular application as it facilitates the identification of regions of space which are not well represented by existing compounds. Several ways to extend the use of partitioned pharmacophore spaces are described. Applications of the approach in monomer acquisition and in library design are outlined.     

Implementation of a system for reagent selection and library enumeration, profiling, and design

We describe an integrated suite of computational tools which are used to assist in the selection of compounds for biological assays and the design of combinatorial libraries. These functions are delivered in a platform-independent manner via a corporate intranet and are used by computational experts and nonexperts alike. While the system was primarily designed to be used prior to synthesis, it can also be used to provide structural information for library registration and for decoding beads in tagged libraries. We describe a simple statistical method for monomer selection and compare it to computationally more demanding approaches.     

Scalable methods for the construction and analysis of virtual combinatorial libraries

One can distinguish between two kinds of virtual combinatorial libraries: viable and accessible . Viable libraries are relatively small in size, are assembled from readily available reagents that have been filtered by the medicinal chemist, and often have a physical counterpart. Conversely, accessible libraries can encompass millions or billions of structures, typically include all possible reagents that are in principle compatible with a particular reaction scheme, and they can never be physically synthesized in their entirety. Although the analysis of viable virtual libraries is relatively straightforward, the handling of large accessible libraries requires methods that scale well with respect to library size. In this work, we present novel, efficient and scalable techniques for the construction, analysis, and in silico screening of massive virtual combinatorial libraries.     

Integrating virtual screening and combinatorial chemistry for accelerated drug discovery

Virtual screening is increasingly being used in drug discovery programs with a growing number of successful applications. Experimental methodologies developed to speed up the drug discovery processes include high-throughput screening and combinatorial chemistry. The complementarities between computational and experimental screenings have been recognized and reviewed in the literature. Computational methods have also been used in the combinatorial chemistry field, in particular in library design. However, the integration of computational and combinatorial chemistry screenings has been attempted only recently. Combinatorial libraries (experimental or virtual) represent a notable source of chemically related compounds. Advances in combinatorial chemistry and deconvolution strategies, have enabled the rapid exploration of novel and dense regions in the chemical space. The present review is focused on the integration of virtual and experimental screening of combinatorial libraries. Applications of virtual screening to discover novel anticancer agents and our ongoing efforts towards the integration of virtual screening and combinatorial chemistry are also discussed.     

A probabilistic approach to high throughput drug discovery

A methodology is presented in which high throughput screening experimental data are used to construct a probabilistic QSAR model which is subsequently used to select building blocks for a virtual combinatorial library. The methodology is based upon statistical probability estimation and not regression. The methodology is applied to the construction of two focused virtual combinatorial libraries: one for cyclic GMP phosphodiesterase type V inhibitors and one for acyl-CoA:cholesterol O-acyltransferase inhibitors. The results suggest that the methodology is capable of selecting combinatorial substituents that lead to active compounds starting with binary (pass/fail) activity measurements.     

A Cleavable Scaffold Strategy for the Synthesis of One-Bead One-Compound Cyclic Peptoid Libraries That Can Be Sequenced By Tandem Mass Spectrometry

Many macrocyclic depsipeptides or related compounds have interesting medicinal properties and often display more favorable pharmacokinetic properties than linear analogues. Therefore, there is considerable interest in the development of large combinatorial libraries of macrocyclic peptidomimetic compounds. However, such molecules cannot be easily sequenced by tandem mass spectrometry, making it difficult to identify hits isolated from library screens using one bead one compound libraries. Here we report a strategy to solve this problem by placing a methionine in both the linker connecting the cyclic molecule to the bead as well as within the cycle itself. Treatment with CNBr both linearizes the molecule at a specific position and releases the molecule from the bead, making its characterization by tandem MALDI mass spectrometry straightforward.     

Fragment-based drug design: computational & experimental state of the art

Fragment-based screening is an emerging technology which is used as an alternative to high-throughput screening (HTS), and often in parallel. Fragment screening focuses on very small compounds. Because of their small size and simplicity, fragments exhibit a low to medium binding affinity (mM to μM) and must therefore be screened at high concentration in order to detect binding events. Since some issues are associated with high-concentration screening in biochemical assays, biophysical methods are generally employed in fragment screening campaigns. Moreover, these techniques are very sensitive and some of them can give precise information about the binding mode of fragments, which facilitates the mandatory hit-to-lead optimization. One of the main advantages of fragment-based screening is that fragment hits generally exhibit a strong binding with respect to their size, and their subsequent optimization should lead to compounds with better pharmacokinetic properties compared to molecules evolved from HTS hits. In other words, fragments are interesting starting points for drug discovery projects. Besides, the chemical space of low-complexity compounds is very limited in comparison to that of drug-like molecules, and thus easier to explore with a screening library of limited size. Furthermore, the "combinatorial explosion" effect ensures that the resulting combinations of interlinked binding fragments may cover a significant part of "drug-like" chemical space. In parallel to experimental screening, virtual screening techniques, dedicated to fragments or wider compounds, are gaining momentum in order to further reduce the number of compounds to test. This article is a review of the latest news in both experimental and in silico virtual screening in the fragment-based discovery field. Given the specificity of this journal, special attention will be given to fragment library design.     

Immunoassays: biological tools for high throughput screening and characterisation of combinatorial libraries

In the demanding field of proteomics, there is an urgent need for affinity-catcher molecules to implement effective and high throughput methods for analysing the human proteome or parts of it. Antibodies have an essential role in this endeavour, and selection, isolation and characterisation of specific antibodies represent a key issue to meet success. Alternatively, it is expected that new, well-characterised affinity reagents generated in rapid and cost-effective manners will also be used to facilitate the deciphering of the function, location and interactions of the high number of encoded protein products. Combinatorial approaches combined with high throughput screening (HTS) technologies have become essential for the generation and identification of robust affinity reagents from biological combinatorial libraries and the lead discovery of active/mimic molecules in large chemical libraries. Phage and yeast display provide the means for engineering a multitude of antibody-like molecules against any desired antigen. The construction of peptide libraries is commonly used for the identification and characterisation of ligand-receptor specific interactions, and the search for novel ligands for protein purification. Further improvement of chemical and biological resistance of affinity ligands encouraged the "intelligent" design and synthesis of chemical libraries of low-molecular-weight bio-inspired mimic compounds. No matter what the ligand source, selection and characterisation of leads is a most relevant task. Immunological assays, in microtiter plates, biosensors or microarrays, are a biological tool of inestimable value for the iterative screening of combinatorial ligand libraries for tailored specificities, and improved affinities. Particularly, enzyme-linked immunosorbent assays are frequently the method of choice in a large number of screening strategies, for both biological and chemical libraries.     

Scaffold architecture and pharmacophoric properties of natural products and trade drugs: application in the design of natural product-based combinatorial libraries

Natural products were analyzed to determine whether they contain appealing novel scaffold architectures for potential use in combinatorial chemistry. Ring systems were extracted and clustered on the basis of structural similarity. Several such potential scaffolds for combinatorial chemistry were identified that are not present in current trade drugs. For one of these scaffolds a virtual combinatorial library was generated. Pharmacophoric properties of natural products, trade drugs, and the virtual combinatorial library were assessed using a self-organizing map. Obviously, current trade drugs and natural products have several topological pharmacophore patterns in common. These features can be systematically explored with selected combinatorial libraries based on a combination of natural product-derived and synthetic molecular building blocks.     

Automatic tailoring and transplanting: a practical method that makes virtual screening more useful

Docking-based virtual screening of large compound libraries has been widely applied to lead discovery in structure-based drug design. However, subsequent lead optimizations often rely on other types of computational methods, such as de novo design methods. We have developed an automatic method, namely automatic tailoring and transplanting (AutoT&T), which can effectively utilize the outcomes of virtual screening in lead optimization. This method detects suitable fragments on virtual screening hits and then transplants them onto a lead compound to generate new ligand molecules. Binding affinities, synthetic feasibilities, and drug-likeness properties are considered in the selection of final designs. In this study, our AutoT&T program was tested on three different target proteins, including p38 MAP kinase, PPAR-α, and Mcl-1. In the first two cases, AutoT&T was able to produce molecules identical or similar to known inhibitors with better potency than the given lead compound. In the third case, we demonstrated how to apply AutoT&T to design novel ligand molecules from scratch. Compared to the solutions generated by other two de novo design methods, i.e., LUDI and EA-Inventor, the solutions generated by AutoT&T were structurally more diverse and more promising in terms of binding scores in all three cases. AutoT&T also completed the assigned jobs more efficiently than LUDI and EA-Inventor by several folds. Our AutoT&T method has certain technical advantages over de novo design methods. Importantly, it expands the application of virtual screening from lead discovery to lead optimization and thus may serve as a valuable tool for many researchers.