Ligand-based combinatorial design of selective purinergic receptor (A2A) antagonists using self-organizing maps

A virtual screening procedure based on a topological pharmacophore similarity metric and self-organizing maps (SOM) was developed and applied to optimizing combinatorial products functioning as P(1) purinergic receptor antagonists. The target was the human A(2A) receptor. A SOM was developed using a set of biologically tested molecules to establish a preliminary structure-activity relationship. A combinatorial library design was performed by projecting virtually assembled new molecules onto the SOM. A small focused library of 17 selected combinatorial products was synthesized and tested. On average, the designed structures yielded a 3-fold smaller binding constant ( approximately 33 vs approximately 100 nM) and 3.5-fold higher selectivity (50 vs 14) than the initial library. The most selective compound obtained revealed a 121-fold relative selectivity for A(2A) with K(i) (A(2A)) = 2.4 nM, and K(i) (A(1)) = 292 nM. This result demonstrates that it was possible to design a small, activity-enriched focused library with an improved property profile using the SOM virtual screening approach. The strategy might be particularly useful in projects in which structure-based design cannot be applied because of a lack of receptor structure information, for example, in the many projects aiming at finding new GPCR modulators.     

Design of combinatorial libraries for the exploration of virtual hits from fragment space searches with LoFT

A case study is presented illustrating the design of a focused CDK2 library. The scaffold of the library was detected by a feature trees search in a fragment space based on reactions from combinatorial chemistry. For the design the software LoFT (Library optimizer using Feature Trees) was used. The special feature called FTMatch was applied to restrict the parts of the queries where the reagents are permitted to match. This way a 3D scoring function could be simulated. Results were compared with alternative designs by GOLD docking and ROCS 3D alignments.     

Development and evaluation of a generic evolutionary method for protein-ligand docking

We have developed a generic evolutionary method with an empirical scoring function for the protein-ligand docking, which is a problem of paramount importance in structure-based drug design. This approach, referred to as the GEMDOCK (Generic Evolutionary Method for molecular DOCKing), combines both continuous and discrete search mechanisms. We tested our approach on seven protein-ligand complexes, and the docked lowest energy structures have root-mean-square derivations ranging from 0.32 to 0.99 A with respect to the corresponding crystal ligand structures. In addition, we evaluated GEMDOCK on crossdocking experiments, in which some complexes with an identical protein used for docking all crystallized ligands of these complexes. GEMDOCK yielded 98% docked structures with RMSD below 2.0 A when the ligands were docked into foreign protein structures. We have reported the validation and analysis of our approach on various search spaces and scoring functions. Experimental results show that our approach is robust, and the empirical scoring function is simple and fast to recognize compounds. We found that if GEMDOCK used the RMSD scoring function, then the prediction accuracy was 100% and the docked structures had RMSD below 0.1 A for each test system. These results suggest that GEMDOCK is a useful tool, and may systematically improve the forms and parameters of a scoring function, which is one of major bottlenecks for molecular recognition.     

Maximum common binding modes (MCBM): consensus docking scoring using multiple ligand information and interaction fingerprints

Improving the scoring functions for small molecule-protein docking is a highly challenging task in current computational drug design. Here we present a novel consensus scoring concept for the prediction of binding modes for multiple known active ligands. Similar ligands are generally believed to bind to their receptor in a similar fashion. The presumption of our approach was that the true binding modes of similar ligands should be more similar to each other compared to false positive binding modes. The number of conserved (consensus) interactions between similar ligands was used as a docking score. Patterns of interactions were modeled using ligand receptor interaction fingerprints. Our approach was evaluated for four different data sets of known cocrystal structures (CDK-2, dihydrofolate reductase, HIV-1 protease, and thrombin). Docking poses were generated with FlexX and rescored by our approach. For comparison the CScore scoring functions from Sybyl were used, and consensus scores were calculated thereof. Our approach performed better than individual scoring functions and was comparable to consensus scoring. Analysis of the distribution of docking poses by self-organizing maps (SOM) and interaction fingerprints confirmed that clusters of docking poses composed of multiple ligands were preferentially observed near the native binding mode. Being conceptually unrelated to commonly used docking scoring functions our approach provides a powerful method to complement and improve computational docking experiments.     

Library Design Guided by In-house Developed Computer Tools

In recent years, combinatorial library synthesis for drug discovery begins to migrate from library synthesis solely dictated by chemistry availability to design and synthesis of libraries with more drug-like properties. Lipinski's rule of five has been used to evaluate drug-like properties of individual compound; recently LibProTM, a new computation program has been developed at Pharmacopeia to evaluate durg-like properties of libraries. By using LibPrpTM, chemists at Pharmacopeia are able to obtain information of molecular weight and ClogP distribution of a library, and percentage of library members that violate Lipinski's rule after input structures of synthons for each combinatorial step. Currently, a "virtual library design” approach that is to calculate properties of a library at conceptual phase of the library design has been used to predetermine the value of the library. Also a new computer program used to predict "Absorption” of compounds will also be discussed.     

CombiDOCK: Structure-based combinatorial docking and library design

We have developed a strategy for efficiently docking a large combinatorial library into a target receptor. For each scaffold orientation, all potential fragments are attached to the scaffold, their interactions with the receptor are individually scored and factorial combinations of fragments are constructed. To test its effectiveness, this approach is compared to two simple control algorithms. Our method is more efficient than the controls at selecting best scoring molecules and at selecting fragments for the construction of an exhaustive combinatorial library. We also carried out a retrospective analysis of the experimental results of a 10×10×10 exhaustive combinatorial library. An enrichment factor of approximately 4 was found for identifying the compounds in the library that are active at 330 nM.     

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.     

Homology model-based virtual screening for GPCR ligands using docking and target-biased scoring

The current study investigates the combination of two recently reported techniques for the improvement of homology model-based virtual screening for G-protein coupled receptor (GPCR) ligands. First, ligand-supported homology modeling was used to generate receptor models that were in agreement with mutagenesis data and structure-activity relationship information of the ligands. Second, interaction patterns from known ligands to the receptor were applied for scoring and rank ordering compounds from a virtual library using ligand-receptor interaction fingerprint-based similarity (IFS). Our approach was evaluated in retrospective virtual screening experiments for antagonists of the metabotropic glutamate receptor (mGluR) subtype 5. The results of our approach were compared to the results obtained by conventional scoring functions (Dock-Score, PMF-Score, Gold-Score, ChemScore, and FlexX-Score). The IFS lead to significantly higher enrichment rates, relative to the competing scoring functions. Though using a target-biased scoring approach, the results were not biased toward the chemical classes of the reference structures. Our results indicate that the presented approach has the potential to serve as a general setup for successful structure-based GPCR virtual screening.     

Application of self-organizing maps in compounds pattern recognition and combinatorial library design

In the computer-aided drug design, in order to find some new leads from a large library of compounds, the pattern recognition study of the diversity and similarity assessment of the chemical compounds is required; meanwhile in the combinatorial library design, more attention is given to design target focusing library along with diversity and drug-likeness criteria. This review presents the current state-of-art applications of Kohonen self-organizing maps (SOM) for studying the compounds pattern recognition, comparing the property of molecular surfaces, distinguishing drug-like and nondrug-like molecules, splitting a dataset into the proper training and test sets before constructing a QSAR (Quantitative Structural-Activity Relationship) model, and also for the combinatorial libraries comparison and the combinatorial library design. The Kohonen self-organizing map will continue to play an important role in drug discovery and library design.     

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.     

针对蛋白质复合物Other类型的打分函数

在不同类型复合物结合界面的物理化学特征不同的基础上, 针对较难预测的Other 类型复合物设计出特异性打分函数, 用于在对接过程中挑选出有效结构. 该函数由原子接触能(EACE)、范德华和静电相互作用能组成,通过多元线性回归方法获得各项的权重系数. 对来自CAPRI benchmark1 中17 个Other 类复合物例子进行打分测试. 结果表明,组合打分能够刻画出Other 类型复合物单体间相互作用的特征, 反映出复合物形成前后的能量变化, 具备一定的从众多样本中筛选出有效结构的能力. 相对于残基成对势(RP), 该组合打分获得了更高的打分成功率. 对CAPRI 第八轮竞赛中两个结构预测模型进行打分排序, 该组合打分也体现出强于RP 的鉴别有效结合模式潜力.     

Strategy for discovering chemical inhibitors of human cyclophilin a: focused library design, virtual screening, chemical synthesis and bioassay

The discovery of cyclophilin A (CypA) inhibitor is now of special interest in the treatment of immunological disorders. In this work, using a strategy integrating focused combinatorial library design, virtual screening, chemical synthesis, and bioassay, a series of novel small molecular CypA inhibitors have been discovered. First, using the fragments taken from our previously discovered CypA inhibitors (Bioorg. Med. Chem. 2006, 14, 2209-2224) as building blocks, we designed a focused combinatorial library containing 255 molecules employing the LD1.0 program (J. Comb. Chem. 2005, 7, 398-406) developed by us. Sixteen compounds (1a-e, 2a-b, 3a-b, and 4a-g) were selected by using virtual screening against the X-ray crystal structure of CypA as well as druglike analysis for further synthesis and bioassay. All these sixteen molecules are CypA binders with binding affinities (K(D) values) ranging from 0.076 to 41.0 microM, and five of them (4a, 4c, and 4e-g) are potent CypA inhibitors with PPIase inhibitory activities (IC(50) values) of 0.25-6.43 microM. The hit rates for binders and inhibitors are as high as 100% and 31.25%, respectively. Remarkably, both the binding affinity and inhibitory activity of the most potent compound increase approximately 10 times than that of the most active compound discovered previously. The high hit rate and the high potency of the new CypA inhibitors demonstrated the efficiency of the strategy for focused library design and screening. In addition, the novel chemical entities reported in this study could be leads for discovering new therapies against the CypA pathway.     

Pharmacophore identification, in silico screening, and virtual library design for inhibitors of the human factor Xa

Factor Xa inhibitors are innovative anticoagulant agents that provide a better safety/efficacy profile compared to other anticoagulative drugs. A chemical feature-based modeling approach was applied to identify crucial pharmacophore patterns from 3D crystal structures of inhibitors bound to human factor Xa (Pdb entries 1fjs, 1kns, 1eqz) using the software LIGANDSCOUT and CATALYST. The complex structures were selected regarding the criteria of high inhibitory potency (i.e. all ligands show K(i) values against factor Xa in the subnanomolar range) and good resolution (i.e. at least 2.2 A) in order to generate selective and high quality pharmacophore models. The resulting chemical-feature based hypotheses were used for virtual screening of commercial molecular databases such as the WDI database. Furthermore, a ligand-based molecular modeling approach was performed to obtain common-feature hypotheses that represent the relevant chemical interactions between 10 bioactive factor Xa inhibitors and the protein, respectively. In a next step a virtual combinatorial library was designed in order to generate new compounds with similar chemical and spatial properties as known inhibitors. The software tool ILIB DIVERSE was used for this procedure in order to provide new scaffolds of this group of anticoagulants. Finally we present the combination of these two techniques, hence virtual screening was performed with selective pharmacophore models in a focused virtual combinatorial database. De novo derived molecular scaffolds that were able to adequately satisfy the pharmacophore criteria are revealed and are promising templates for candidates for further development.     

A virtual high throughput screen for high affinity cytochrome P450cam substrates. Implications for in silico prediction of drug metabolism

Structure-based virtual screening techniques require reliable scoring functions to discriminate potential substrates effectively. In this study we compared the performance of GOLD, PMF, DOCK and FlexX scoring functions in FlexX flexible docking to cytochrome P450cam binding site. Crystal structures of protein-substrate complexes were most effectively reproduced by the FlexX/PMF method. On the other hand, the FlexX/GOLD approach provided the best correlation between experimental binding constants and predicted scores. Binding modes selected by the FlexX/PMF approach were rescored by GOLD to obtain a reliable measure of binding energetics. The effectiveness of the FlexX/PMF/GOLD method was demonstrated by the correct classification of 32 out of the 33 experimentally studied compounds and also in a virtual HTS test on a library of 10,000 compounds. Although almost all the available functions were developed to be general, our study on cytochrome P450cam substrates suggests that careful selection or even tailoring the scoring function might increase the prediction power of virtual screens significantly. The FlexX/PMF/GOLD methodology was tested on cytochrome P450 3A4 substrates and inhibitors. This preliminary study revealed that the combined function was able to recognise 334 out of the 345 compounds bound to 3A4.     

组合化学法在筛选真空紫外荧光材料中的应用

综述了近年来组合化学法在筛选新型真空紫外荧光材料研究方面取得的进展, 包括组合材料库的并行合成和高通量表征技术, 并重点介绍了组合材料库的设计艺术, 最后列举了组合化学法筛选真空紫外荧光材料体系研究的实例及结果.     

Toward designing drug-like libraries: a novel computational approach for prediction of drug feasibility of compounds

Prediction of the degree of drug-like character in small molecules is of great industrial interest. The major barrier, however, is the lack of a definition for drug-like character. We used the concept of the multilevel chemical compatibility (MLCC) between a compound and a drug library as a measure of the drug-like character of a compound. The rationale is that the local chemical environment of each atom or group of atoms in a compound largely contributes to the stability, toxicity, and metabolism in vivo. A systematic comparison of the local environments within a compound and those within the existing drugs provides a basis for determining whether and how much a compound is drug-like. We applied the MLCC calculations to four test sets: top selling drugs, compounds under biological testing prior to the preclinical test, anticancer drugs, and compounds known to have poor drug-like character. The following conclusions were obtained: (1) A convergent number of unique local structure types were found in the analysis of the library of the existing drugs. It suggests that the current drug library contains about 80% of all the viable types; therefore, discovery of a drug with new local structures is only an event of relatively small probability. (2) The method is highly selective in discerning drug-like compounds: most of the top drugs are predicted to be drug-like, about one-quarter of the biological testing compounds are drug-like, and about one-fifth of the anticancer drugs are drug-like. (3) The method also correctly predicted that none of the known problematic compounds are drug-like. (4) The method is fast enough for computational screening of virtual combinatorial chemistry libraries and databases of available compounds.     

The Application of Combinatorial Chemistry in Catalysis

In recent few years combinatorial methodology has been extensively used in material science research. Based on the desired properties of materials, various high throughput synthesizing and screening technologies were developed. These high throughput technologies can increase our speed to more than hundred folds for finding and optimizing materials. One of the most active areas is catalysis. Scientists are developing novel high throughput technologies to screen catalyst libraries to find and optimize new catalysts for chemical industry. In this area die key is combinatorial catalytic reactor design, catalyst library synthesis, and product detection. Systematic technologies for catalyst library synthesis and characterization were developed in our laboratory. In this work, catalyst in situ synthesis, parallel reactor design, and detection methods will be introduced. Combining with the powerful combinatorial methodology, good chemistry design will make our work even more efficient. Hence, as an example of combining combinatorial technologies with chemistry design, a successful catalyst design is also introduced.