PLUMS: a program for the rapid optimization of focused libraries

PLUMS is a new method to perform rational monomer selection for combinatorial chemistry libraries. The algorithm has been developed to optimize focused libraries with specific two-dimensional and/or three-dimensional properties. A preliminary step is the identification of those molecules in the initial virtual library which satisfy the imposed property constraints; we define these molecules as the virtual hits. From the virtual hits, PLUMS generates a starting library, which is the true combinatorial library that includes all the virtual hits. Monomers are then removed in an iterative fashion, thus reducing the size of the library. At each iteration, the worst monomer is removed. Each sublibrary is selected using a global scoring function, which balances effectiveness and efficiency. The iterative process continues until one is left with a library that consists entirely of virtual hits. The optimal library, which is the best compromise between effectiveness and efficiency, can then be selected according to the score. During the iterative process, equivalent solutions may well occur and are taken into account by the algorithm, according to a user-defined parameter. The number of monomers for each substitution site and the size of the library are parameters that can be either optimized or used to constrain the selection. The results obtained on two test libraries are presented. PLUMS was compared with genetic algorithms (GA) and monomer frequency analysis (MFA), which are widely used for monomer selection. For the two test libraries, PLUMS and GA gave equivalent results. MFA is the fastest method, but it can give misleading solutions. Possible advantages and disadvantages of the different methods are discussed.     

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.     

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.     

Combinatorial library-based design with Basis Products

Uncovering useful lead compounds from a vast virtual library of synthesizable compounds continues to be of tremendous interest to pharmaceutical researchers. Here we present the concept of Basis Products (BPs), a new and broadly applicable method for achieving efficient selections from a combinatorial library. By definition, Basis Products are a strategically selected subset of compounds from a potentially very large combinatorial library, and any compound in a combinatorial library can represented by its BPs. In this article we will show how to use BP docking scores to find the top compounds of a combinatorial library. Compared with the brute-force docking of an entire virtual library, docking with BPs are much more efficient because of the substantial size reduction, saving both time and resources. We will also demonstrate how BPs can be used for property-based combinatorial library designs. Furthermore, BPs can also be considered as fragments carrying chemistry knowledge, hence they can potentially be used in combination with any fragment-based design method. Therefore, BPs can be used to integrate combinatorial design with structure-based design and/or fragment-based design. Other potential applications of BPs include lead hopping and consensus core building, which we will describe briefly as well in this report.     

OptiDock: virtual HTS of combinatorial libraries by efficient sampling of binding modes in product space

Products from combinatorial libraries generally share a common core structure that can be exploited to improve the efficiency of virtual high-throughput screening (vHTS). In general, it is more efficient to find a method that scales with the total number of reagents (Sigma growth) rather with the number of products (Pi growth). The OptiDock methodology described herein entails selecting a diverse but representative subset of compounds that span the structural space encompassed by the full library. These compounds are docked individually using the FlexX program (Rarey, M.; Kramer, B.; Lengauer, T.; Klebe, G. J. Mol. Biol. 1995, 251, 470-489) to define distinct docking modes in terms of reference placements for combinatorial core atoms. Thereafter, substituents in R-cores (consisting of the core structure substituted at a single variation site) are docked, keeping the core atoms fixed at the coordinates dictated by each reference placement. Interaction energies are calculated for each docked R-core with respect to the target protein, and energies for whole compounds are calculated by finding the reference core placement for which the sum of corresponding R-core energies is most negative. The use of diverse whole compounds to define binding modes is a key advantage of the protocol over other combinatorial docking programs. As a result, OptiDock returns better-scoring conformers than does serially applied FlexX. OptiDock is also better able to find a viable docked pose for each library member than are other combinatorial approaches.     

Design of libraries to explore receptor sites

Despite rapid progress in both combinatorial chemistry and high-throughput screening, the number of molecules that could potentially be made and tested for biological activity still far exceeds the capacity for synthesis or screening. Consequently, it is potentially valuable to select and synthesize sublibraries that contain rationally selected subsets. When the structure of the protein receptor site is known, this may be used to impose restrictions of the selection on molecules. This paper describes a method for rapid analysis of large virtual libraries to select a subset that can exhibit at least one conformer which wil interact strongly with the receptor and fit within the receptor site.     

Recent developments in the HPLC enantiomeric separation using chiral selectors identified by a combinatorial strategy

Two kinds of chiral stationary phases (CSPs) can be distinguished: (i) the "conventional" CSPs for which the selectivity is not pre-determined and (ii) the CSPs which are characterized by a predictable elution order, depending on the target enantiomer used for the selection of the chiral selector. At the present time, three general methodologies have been described to create chiral selectors specifically designed against the racemate to resolve: the molecular imprinting technology, the production of antibodies and the combinatorial approach. The latter methodology involves two categories of procedures to develop CSPs: an approach from a small library of low-molecular weight selectors and an approach from a very large library of single-stranded oligonucleotides (DNA and RNA aptamers). In this review, the recent advances in the HPLC applications of the chiral selectors identified through these two combinatorial procedures are addressed.     

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.     

Comparative molecular surface analysis (CoMSA) for virtual combinatorial library screening of styrylquinoline HIV-1 blocking agents

We used comparative molecular surface analysis to design molecules for the synthesis as part of the search for new HIV-1 integrase inhibitors. We analyzed the virtual combinatorial library (VCL) constituted from various moieties of styrylquinoline and styrylquinazoline inhibitors. Since imines can be applied in a strategy of dynamic combinatorial chemistry (DCC), we also tested similar compounds in which the -C=N- or -N=C- linker connected the heteroaromatic and aromatic moieties. We then used principal component analysis (PCA) or self-organizing maps (SOM), namely, the Kohonen neural networks to obtain a clustering plot analyzing the diversity of the VCL formed. Previously synthesized compounds of known activity, used as molecular probes, were projected onto this plot, which provided a set of promising virtual drugs. Moreover, we further modified the above mentioned VCL to include the single bond linker -C-N- or -N-C-. This allowed increasing compound stability but expanded also the diversity between the available molecular probes and virtual targets. The application of the CoMSA with SOM indicated important differences between such compounds and active molecular probes. We synthesized such compounds to verify the computational predictions.     

Virtual darwinian drug design: QSAR inverse problem, virtual combinatorial chemistry, and computational screening

The generation of diversity and its further selection by an external system is a common mechanism for the evolution of the living species and for the current drug design methods. This assumption allows us to label the methods based on generation and selection of molecular diversity as "Darwinian" ones, and to distinguish them from the structure-based, structure-modulation approaches. An example of a Darwinian method is the inverse QSAR. It consists of the computational generation of candidate chemical structures and their selection according to a previously established QSAR model. New trends in the field of combinatorial chemical syntheses comprise the concepts of virtual combinatorial synthesis and virtual or computational screening. Virtual combinatorial synthesis, closely related to inverse QSAR, can be defined as the computational simulation of the generation of new chemical structures by using a combinatorial strategy to generate a virtual library. Virtual screening is the selection of chemical structures having potential desirable properties from a database or virtual library in order to be synthesized and assayed. This review is mainly focused on graph theoretical drug design approaches, but a survey with key references is provided that covers other simulation methods.     

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.     

Library fingerprints: a novel approach to the screening of virtual libraries

We propose a novel method to prioritize libraries for combinatorial synthesis and high-throughput screening that assesses the viability of a particular library on the basis of the aggregate physical-chemical properties of the compounds using a na?ve Bayesian classifier. This approach prioritizes collections of related compounds according to the aggregate values of their physical-chemical parameters in contrast to single-compound screening. The method is also shown to be useful in screening existing noncombinatorial libraries when the compounds in these libraries have been previously clustered according to their molecular graphs. We show that the method used here is comparable or superior to the single-compound virtual screening of combinatorial libraries and noncombinatorial libraries and is superior to the pairwise Tanimoto similarity searching of a collection of combinatorial libraries.     

A fast exchange algorithm for designing focused libraries in lead optimization

Combinatorial chemistry is widely used in drug discovery. Once a lead compound has been identified, a series of R-groups and reagents can be selected and combined to generate new potential drugs. The combinatorial nature of this problem leads to chemical libraries containing usually a very large number of virtual compounds, far too large to permit their chemical synthesis. Therefore, one often wants to select a subset of "good" reagents for each R-group of reagents and synthesize all their possible combinations. In this research, one encounters some difficulties. First, the selection of reagents has to be done such that the compounds of the resulting sublibrary simultaneously optimize a series of chemical properties. For each compound, a desirability index, a concept proposed by Harrington,(20) is used to summarize those properties in one fitness value. Then a loss function is used as objective criteria to globally quantify the quality of a sublibrary. Second, there are a huge number of possible sublibraries, and the solutions space has to be explored as fast as possible. The WEALD algorithm proposed in this paper starts with a random solution and iterates by applying exchanges, a simple method proposed by Fedorov(13) and often used in the generation of optimal designs. Those exchanges are guided by a weighting of the reagents adapted recursively as the solutions space is explored. The algorithm is applied on a real database and reveals to converge rapidly. It is compared to results given by two other algorithms presented in the combinatorial chemistry literature: the Ultrafast algorithm of D. Agrafiotis and V. Lobanov and the Piccolo algorithm of W. Zheng et al.     

Web-based cheminformatics and molecular property prediction tools supporting drug design and development at Novartis

Web-based tools offer many advantages for processing chemical information, most notably ease of use and high interactivity. Therefore more and more pharmaceutical companies are using web technology to deliver sophisticated molecular processing tools directly to the desks of their chemists, to assist them in the process of designing and developing new drugs. In this paper, the web-based cheminformatics system developed at Novartis and currently used by more than thousand users is described. The system allows various molecular modeling and molecular processing tasks, including the calculation of molecular and substituent properties, property-based virtual screening, visualization of molecules, bioisosteric design, diversity analysis, and support of combinatorial chemistry. The methodology to calculate various molecular properties relevant to drug design is described, including the prediction of intestinal absorption, blood-brain barrier penetration, efflux, and water solubility. Information about the web technology used is also provided.     

Successful screening of large encoded combinatorial libraries leading to the discovery of novel p38 MAP kinase inhibitors

Screening of more than 2 million compounds comprising 41 distinct encoded combinatorial libraries revealed a novel structural class of p38 mitogen-activated protein (MAP) kinase inhibitors. The methodology used for screening large encoded combinatorial libraries combined with the statistical interpretation of screening results is described. A strong preference for a particular triaminotriazine aniline amide was discovered based on biological activity observed in the screening campaign. Additional screening of a focused follow-up combinatorial library yielded data expanding the unique combinatorial SAR and emphasizing an extraordinary preference for this particular building block and structural class. The preference is further highlighted when the p38 inhibitor data set is compared to data obtained for a panel of other kinases.     

RNA aptamers as pathway-specific MAP kinase inhibitors

BACKGROUND: In eukaryotic cells, many intracellular signaling pathways have closely related mitogen activated protein kinase (MAPK) paralogs as central components. Although MAPKs are therefore obvious targets to control the cellular responses resulting from the activation of these signaling pathways, the development of inhibitors which target specific cell signaling pathways involving MAPKs has proven difficult. RESULTS: We used an RNA combinatorial approach to isolate RNAs that inhibit the in vitro phosphorylation activity of extracellular regulated kinase 2 (ERK2). These inhibitors block phosphorylation by ERK1 and ERK2, but do not inhibit Jun N-terminal kinase or p38 MAPKs. Kinetic analysis indicates these inhibitors function at high picomolar concentrations through the steric exclusion of substrate and ATP binding. In one case, we identified a compact RNA structural domain responsible for inhibition. CONCLUSIONS: RNA reagents can selectively recognize and inhibit MAPKs involved in a single signal transduction pathway. The methodology described here is readily generalizable, and can be used to develop inhibitors of MAPKs involved in other signal transduction pathways. Such reagents may be valuable tools to analyze and distinguish homologous effectors which regulate distinct signaling responses.