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.     

Development of a system to evaluate compound identity, purity, and concentration in a single experiment and its application in quality assessment of combinatorial libraries and screening hits

The development and use of a new assay system for the simultaneous determination of identity, purity, and concentration of sample components from combinatorial libraries produced by parallel synthesis are described. The system makes use of high-performance liquid chromatography with UV/vis photodiode array (PDA), evaporative light scattering (ELSD), chemiluminescent nitrogen (CLND), and time-of-flight mass spectrometer (TOFMS) detectors (HPLC-PDA-ELSD-CLND-TOFMS). Although these detectors have previously been utilized separately for the analysis of combinatorial chemistry libraries, the use of TOFMS along with CLND provides a synergistic combination enabling target and side-product structures to be identified and their concentrations and purities determined in a single experiment from a solution containing microgram levels of material. The CLND was found to give a linear response based on the number of moles of nitrogen present. Therefore, if the number of nitrogens per molecule is known, the concentration of each nitrogen-containing sample component may be determined utilizing an unrelated co-injected standard. A molecular formula for an impurity may often be calculated from the exact mass determined by the TOFMS and knowledge of the chemistry involved. Thus, if the sample components contain nitrogen, the concentration of every identified HPLC peak may be determined even in the absence of primary standards. This combination of detectors enabled the characterization of both target compounds and byproducts in combinatorial libraries, allowing the optimization of library synthetic procedures. This system was also used to survey the quality of libraries, enabling the selection of the best libraries for screening. This method also facilitated the characterization of samples from combinatorial libraries found as hits in high-throughput screening to establish the potency of the leads based on their actual concentration. In addition, concentrations and potencies of impurities were determined after identification of their structures, utilizing exact mass data, determination of charge states, and knowledge of the synthetic chemistry.     

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.     

How large does a compound screening collection need to be?

Increasingly, chemical libraries are being produced which are focused on a biological target or group of related targets, rather than simply being constructed in a combinatorial fashion. A screening collection compiled from such libraries will contain multiple analogues of a number of discrete series of compounds. The question arises as to how many analogues are necessary to represent each series in order to ensure that an active series will be identified. Based on a simple probabilistic argument and supported by in-house screening data, guidelines are given for the number of compounds necessary to achieve a "hit", or series of hits, at various levels of certainty. Obtaining more than one hit from the same series is useful since this gives early acquisition of SAR (structure-activity relationship) and confirms a hit is not a singleton. We show that screening collections composed of only small numbers of analogues of each series are sub-optimal for SAR acquisition. Based on these studies, we recommend a minimum series size of about 200 compounds. This gives a high probability of confirmatory SAR (i.e. at least two hits from the same series). More substantial early SAR (at least 5 hits from the same series) can be gained by using series of about 650 compounds each. With this level of information being generated, more accurate assessment of the likely success of the series in hit-to-lead and later stage development becomes possible.     

From protein domains to drug candidates-natural products as guiding principles in the design and synthesis of compound libraries

In the continuing effort to find small molecules that alter protein function and ultimately might lead to new drugs, combinatorial chemistry has emerged as a very powerful tool. Contrary to original expectations that large libraries would result in the discovery of many hit and lead structures, it has been recognized that the biological relevance, design, and diversity of the library are more important. As the universe of conceivable compounds is almost infinite, the question arises: where is a biologically validated starting point from which to build a combinatorial library? Nature itself might provide an answer: natural products have been evolved to bind to proteins. Recent results in structural biology and bioinformatics indicate that the number of distinct protein families and folds is fairly limited. Often the same structural domain is used by many proteins in a more or less modified form created by divergent evolution. Recent progress in solid-phase organic synthesis has enabled the synthesis of combinatorial libraries based on the structure of complex natural products. It can be envisioned that natural-product-based combinatorial synthesis may permit hit or lead compounds to be found with enhanced probability and quality.     

Mass spectrometry and combinatorial chemistry: a short outline.

The rapid evolution of combinatorial chemistry in recent years has led to a dramatic improvement in synthetic capabilities. The goal is to accelerate the discovery of molecules showing affinity against a target, such as an enzyme or a receptor, through the simultaneous synthesis of a great number of structurally diverse compounds. This is done by generating combinatorial libraries containing as many as hundreds or thousands of compounds. The need to test all these compounds led to the development of high-throughput screening (HTS) techniques, and also high-throughput analytical techniques capable of assessing the occurrence, structure and purity of the products. In order to be applied effectively to the characterization of combinatorial libraries, an analytical technique must be adequately sensitive (to analyse samples which are typically produced in nanomole amounts or less), fast, affordable and easy to automate (to minimize analysis time and operator intervention). Although no method alone can meet all the analytical challenges underlying this task, the recent progress in mass spectrometric (MS) instrumentation renders this technique an essential tool for scientists working in this area. We describe here relevant aspects of the use of MS in combinatorial technologies, such as current methods of characterization, purification and screening of libraries. Some examples from our laboratory deal with the analysis of pooled oligomeric libraries containing n x 324(n = 1, 2) compounds, using both on-line high-performance liquid chromatography/MS with an ion trap mass spectrometer, and direct infusion into a triple quadrupole instrument. In the first approach, MS and product ion MS/MS with automatic selection of the precursor were performed in one run, allowing library confirmation and structural elucidation of unexpected by-products. The second approach used MS scans to characterize the entire library and also precursor ion and neutral loss scans to detect selectively components with given structural characteristics.     

Preparative high-performance liquid chromatography-mass spectrometry for the high-throughput purification of combinatorial libraries

Accurate results for the testing of combinatorial libraries necessitates high purity of the library members. Therefore, combinatorial libraries derived from a combinatorial solution or solid-phase synthesis often require the purification of compounds that do not achieve a certain purity threshold. This study describes that preparative high-performance liquid chromatography (HPLC)-mass spectrometry (MS) is the method of choice for the purification of large arrays of diverse compounds. The adoption of this technology to the workflow of a solution phase combinatorial chemistry laboratory producing more than 20,000 compounds per year is described. Furthermore, the setup and logistics are discussed as well as the purity achievable for large libraries. Efficiency, speed, quality, and universality of preparative HPLC-MS are presented in detail for a library of 140 compounds, including data logistics and downstream processes as well.     

Molecular libraries in liquid phase via UGI-MCR

Standard chemistry prescribes the conversion of one or two compounds into their products. In contrast, Eintopf (one-pot) multicomponent reactions (MCRs) involve at least three different compounds. One-pot MCRs are a useful tool in combinatorial chemistry: From a mixture of educts a large number of products can be simultaneously formed in liquid phase, called a soluble molecular library. The member compounds of such libraries are investigated simultaneously for desired properties, e.g. antibiotic activity. The main constraint is, that the underlying chemistry must not produce unknown side reactions and must lead to a broad spectrum of stable products with high yields. Isocyanide multicomponent chemistry allows the generation of soluble libraries of very different sizes, which are easy to screen for biological or pharmaceutical efficacy using the algorithms presented. Products can easily be enumerated and the kinetics of the isocyanide chemistry is simple to investigate. Combinatorial chemistry is capable of generating and optimizing leads faster and with fewer resources than by conventional means. Combinatorial chemistry based on isocyanide chemistry is by far the most important and most impressive technique in use today to reducing time and costs associated with lead generation and optimization during the drug discovery process. The simplicity of the reaction conditions involved means that the generation and screening of libraries can be automated.     

Quality assessment of combinatorial pooled libraries using mass spectrometry

Parallel synthesis techniques aim to prepare collections of single compounds which, once tested, can easily be identified by their sole location in the synthesic array. On the other hand, true combinatorial chemistry produces libraries of compounds as mixtures of variable size which require a deconvolution procedure for identification of the active hits or leads. In the latter case, analytical methods are crucial for the success of the strategy and mass spectrometry plays a major role. If the goal is to identify all the library components, including expected products as well as by-products, various mass spectrometric techniques may be necessary. Library components can be separated according to their mass by increasing mass resolution or by their elution time by coupling liquid chromatography and mass spectrometry. The efficiency of such separation techniques are discussed as a function of the size and the degeneracy of the library. Library members possess common structural features which impart similar fragmentation patterns after ionization in the gas phase. This feature can be exploited by tandem mass spectrometry to specifically detect subfamilies of products. Examples of precursor ion scans, product ion scans and constant neutral loss scans will be shown that facilitate partial characterization of libraries. To solve the difficult problem of the quantitative analysis of libraries, i.e., to evaluate their equimolarity, the use of an evaporative light scattering detector (ELSD) or a chemiluminescent nitrogen detector (CLND) is suggested as more appropriate.     

Combinatorial Chemistry—Challenge and Chance for the Development of New Catalysts and Materials

One should not underestimate the capability of the combinatorial method in solid-state chemistry; this is the opinion of the author. Combinatorial chemistry can provide a large number of new compounds, but once the components that are interesting for a certain application have been successfully selected, the techniques of conventional catalysis and materials research are required. The strengths of conventional chemistry lie in the optimization, systematic modification, and improvement of new lead structures. In contrast, discovery is the potential strength of combinatorial chemistry. Careful design is most important for the synthesis of useful libraries, since the diversity of the periodic table is much too large to be accessed comprehensively or systematically by such large libraries.     

Recent developments in the combinatorial synthesis of nitrogen heterocycles using solid phase technology

Recognizing the potential of combinatorial chemistry to accelerate drug discovery and development, most pharmaceutical and related industries are seriously looking toward combinatorial synthesis of compounds in order to facilitate the identification of 'lead' molecules. In particular, solid phase synthesis is the core technology for combinatorial chemistry and is widely used for generating libraries of structurally related compounds. Since many drugs contain the nitrogen heterocyclic component and since heterocycles possess a high order of structural diversity, a precise overview of recent progress in the combinatorial synthesis of nitrogen heterocycles using solid phase methodology would be useful. Since the progress in solid phase synthesis of organic molecules has been reviewed regularly from 1992 to 1998, only the development of solid phase combinatorial synthetic approaches of small nitrogen heterocycles since 1999 will be reviewed here. This review describes the solid phase synthesis of azepanes, benzodiazepines, benzimidazoles, benzothiazepines, cinnolines, indolizines, beta lactams, oxazepins, oxazoles including benzisooxazoles, hydantoins, piperidines, pyrimidines, pyrazolones, quinolones, trizolopyridazines and thiazoles.     

Novel affinity ligands for chromatography using combinatorial chemistry

Spatially addressable combinatorial libraries were synthesized by solution phase chemistry and screened for binding to human serum albumin. Members of arylidene diamide libraries were among the best hits found, having submicromolar binding affinities. The results were analyzed by the frequency with which particular substituents appeared among the most potent compounds. After immobilization of the ligands either through the oxazolone or the amine substituent, characterization by surface plasmon resonance showed that ibuprofen affected the binding kinetics, but phenylbutazone did not. It is therefore likely that these compounds bind to Site 2 in sub domain IIIA of human serum albumin (HSA).     

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.     

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.     

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.     

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.     

组合化学、分子库与新药研究

组合化学是进入90 年代以来寻找及优化新药先导化合物的主要研究方法, 其特点是改变了传统的逐一合成、逐一纯化、逐一筛选的模式, 而是以合成和筛选化学库的形式完成寻找及优化药物先导化合物, 极大地加快了药物先导化合物出现的速度。本文就目前有关组合化学研究的基本理论、基本方法、发展趋势、研究成果以及我国应当采取的措施进行了综述。     

Streamlined approach to high-quality purification and identification of compound series using high-resolution MS and NMR

Automated medicinal chemistry (parallel chemistry) has become an integral part of the drug-discovery process in almost every large pharmaceutical company. Parallel array synthesis of individual organic compounds has been used extensively to generate diverse structural libraries to support different phases of the drug-discovery process, such as hit-to-lead, lead finding, or lead optimization. In order to guarantee effective project support, efficiency in the production of compound libraries has been maximized. As a consequence, also throughput in chromatographic purification and analysis has been adapted. As a recent trend, more laboratories are preparing smaller, yet more focused libraries with even increasing demands towards quality, i.e. optimal purity and unambiguous confirmation of identity. This paper presents an automated approach how to combine effective purification and structural conformation of a lead optimization library created by microwave-assisted organic synthesis. The results of complementary analytical techniques such as UHPLC-HRMS and NMR are not only regarded but even merged for fast and easy decision making, providing optimal quality of compound stock. In comparison with the previous procedures, throughput times are at least four times faster, while compound consumption could be decreased more than threefold.     

Quality control in combinatorial chemistry: determination of the quantity,purity, and quantitative purity of compounds in combinatorial libraries

The quality of combinatorial libraries determines the success of biological screening in drug discovery programs. In this paper, we evaluate and compare various methods for measuring identity, purity, and quantity (yield) of combinatorial libraries. Determination of quantitative purity reveals the true library quality and often indicates potential quality problems before full-scale library production. The relative purity can be determined for every member in a large library in a high-throughput mode, but must be cautiously interpreted. In particular, many impurities are not observable by relative purity measurements using detectors such as UV(214), UV(254), and evaporative light-scattering detection. These "invisible" impurities may constitute a significant portion of the sample weight. We found that TFA, plastic extracts, inorganic compounds, and resin washout are among these impurities. With compelling evidence, we reach a conclusion that purification is the only way to remove "invisible" impurities and improve the quantitative purity of any compound even though some compounds may have a high relative purity before purification.