Pharmacophore Modeling and in Silico Screening Studies to Design Potential KDR Kinase Inhibitors

A novel ligand‐based pharmacophore model for KDR kinase was generated on the basis of chemical features of 30 KDR kinase inhibitors. This pharmacophore model consists of one hydrogen‐bond acceptor, one hydrogen‐bond donor and two hydrophobic groups. Several methods have been used to validate the model, suggesting that it can serve as a reliable tool for virtual screening to facilitate the discovery of novel KDR inhibitors. The model was then used as database search query from the National Cancer Institute (NCI) database for the rational design to identify new hit compound.     

Rapid Identification of Potential Drug Candidates from Multi-Million Compounds’ Repositories. Combination of 2D Similarity Search with 3D Ligand/Structure Based Methods and In Vitro Screening

Rapid in silico selection of target focused libraries from commercial repositories is an attractive and cost-effective approach in early drug discovery. If structures of active compounds are available, rapid 2D similarity search can be performed on multimillion compounds’ databases. This approach can be combined with physico-chemical parameter and diversity filtering, bioisosteric replacements, and fragment-based approaches for performing a first round biological screening. Our objectives were to investigate the combination of 2D similarity search with various 3D ligand and structure-based methods for hit expansion and validation, in order to increase the hit rate and novelty. In the present account, six case studies are described and the efficiency of mixing is evaluated. While sequentially combined 2D/3D similarity approach increases the hit rate significantly, sequential combination of 2D similarity with pharmacophore model or 3D docking enriched the resulting focused library with novel chemotypes. Parallel integrated approaches allowed the comparison of the various 2D and 3D methods and revealed that 2D similarity-based and 3D ligand and structure-based techniques are often complementary, and their combinations represent a powerful synergy. Finally, the lessons we learnt including the advantages and pitfalls of the described approaches are discussed.     

New serotonin 5-HT(6) ligands from common feature pharmacophore hypotheses

Serotonin 5-HT6 receptor antagonists are thought to play an important role in the treatment of psychiatry, Alzheimer's disease, and probably obesity. To find novel and potent 5-HT6 antagonists and to provide a new idea for drug design, we used a ligand-based pharmacophore to perform the virtual screening of a commercially available database. A three-dimensional common feature pharmacophore model was developed by using the HipHop program provided in Catalyst software and was used as a query for screening the database. A recursive partitioning (RP) model which can separate active and inactive compounds was used as a filtering system. Finally a sequential virtual screening procedure (SQSP) was conducted, wherein both the common feature pharmacophore and the RP model were used in succession to improve the results. Some of the hits were selected based on druglikeness, ADME properties, structural diversity, and synthetic accessibility for real biological evaluation. The best hit compound showed a significant IC50 value of 9.6 nM and can be used as a lead for further drug development.     

Integrating medicinal chemistry, organic/combinatorial chemistry, and computational chemistry for the discovery of selective estrogen receptor modulators with Forecaster, a novel platform for drug discovery

As part of a large medicinal chemistry program, we wish to develop novel selective estrogen receptor modulators (SERMs) as potential breast cancer treatments using a combination of experimental and computational approaches. However, one of the remaining difficulties nowadays is to fully integrate computational (i.e., virtual, theoretical) and medicinal (i.e., experimental, intuitive) chemistry to take advantage of the full potential of both. For this purpose, we have developed a Web-based platform, Forecaster, and a number of programs (e.g., Prepare, React, Select) with the aim of combining computational chemistry and medicinal chemistry expertise to facilitate drug discovery and development and more specifically to integrate synthesis into computer-aided drug design. In our quest for potent SERMs, this platform was used to build virtual combinatorial libraries, filter and extract a highly diverse library from the NCI database, and dock them to the estrogen receptor (ER), with all of these steps being fully automated by computational chemists for use by medicinal chemists. As a result, virtual screening of a diverse library seeded with active compounds followed by a search for analogs yielded an enrichment factor of 129, with 98% of the seeded active compounds recovered, while the screening of a designed virtual combinatorial library including known actives yielded an area under the receiver operating characteristic (AU-ROC) of 0.78. The lead optimization proved less successful, further demonstrating the challenge to simulate structure activity relationship studies.     

Development of predictive pharmacophore model for in silico screening,and 3D QSAR CoMFA and CoMSIA studies for lead optimization,for designing of potent tumor necrosis factor alpha converting enzyme inhibitors

A chemical feature-based pharmacophore model was developed for Tumor Necrosis Factor-α converting enzyme (TACE) inhibitors. A five point pharmacophore model having two hydrogen bond acceptors (A), one hydrogen bond donor (D) and two aromatic rings (R) with discrete geometries as pharmacophoric features was developed. The pharmacophore model so generated was then utilized for in silico screening of a database. The pharmacophore model so developed was validated by using four compounds having proven TACE inhibitory activity which were grafted into the database. These compounds mapped well onto the five listed pharmacophoric features. This validated pharmacophore model was also used for alignment of molecules in CoMFA and CoMSIA analysis. The contour maps of the CoMFA/CoMSIA models were utilized to provide structural insight for activity improvement of potential novel TACE inhibitors. The pharmacophore model so developed could be used for in silico screening of any commercial/in house database for identification of TACE inhibiting lead compounds, and the leads so identified could be optimized using the developed CoMSIA model. The present work highlights the tremendous potential of the two mutually complementary ligand-based drug designing techniques (i.e. pharmacophore mapping and 3D-QSAR analysis) using TACE inhibitors as prototype biologically active molecules.     

Evaluation of potential inhibitors of squalene synthase based on virtual screening and in vitro studies

Squalene synthase (SQS) is a potential target for hyperlipidemia treatment. To identify novel chemical scaffolds of SQS inhibitors, we generated 3D-QSAR pharmacophore models using HypoGen. The best quantitative pharmacophore model, Hypo 1, was selected for virtual screening using two chemical databases, Specs and Traditional Chinese Medicine database (TCM). The best-mapped hit compounds were then subjected to filtering by Lipinski^’s rule of five and docking studies to refine the hits. Finally, five compounds were selected from the top-ranked hit compounds for SQS inhibitory assay in vitro. Three of these compounds could inhibit SQS in vitro, and should be further evaluated pre-clinically as a treatment for hyperlipidemia.     

Pharmacophore features distributions in different classes of compounds

A pharmacophore analysis approach was used to investigate and compare different classes of compounds relevant to the drug discovery process (specifically, drug molecules, compounds in high throughput screening libraries, combinatorial chemistry building blocks and nondrug molecules). The distributions for a set of pharmacophore features including hydrogen bond acceptors, hydrogen bond donors, negatively ionizable centers, positively ionizable centers and hydrophobic points, were generated and examined. Significant differences were observed between the pharmacophore profiles obtained for the drug molecules and those obtained for the high-throughput screening compounds, which appear to be closely related to the nondrug pharmacophore distribution. It is suggested that the analysis of pharmacophore profiles could be used as an additional tool for the property-based optimization of compound selection and library design processes, thus improving the odds of success in lead discovery projects.     

Computational validation of the importance of absolute stereochemistry in virtual screening

Consideration of stereochemistry early in the identification and optimization of lead compounds can improve the efficiency and efficacy of the drug discovery process and reduce the time spent on subsequent drug development. These improvements can result by focusing on specific enantiomers that have the desired potential therapeutic effect (eutomers), while removing from consideration enantiomers that may have no, or even undesirable, effects (distomers). A virtual screening campaign that correctly takes stereochemical information into account can, in theory, be utilized to provide information about the relative binding affinities of enantiomers. Thus, the proper enumeration of the relevant stereoisomers in general, and enantiomeric pairs in particular, of chiral compounds is crucial if one is to use virtual screening as an effective drug discovery tool. As is obvious, in cases where no stereochemical information is provided for chiral compounds in a 2D chemical database, then each possible stereoisomer should be generated for construction of the subsequent 3D database to be used for virtual screening. However, acute problems can arise in 3D database construction when relative stereochemistry is encoded in a 2D database for a chiral compound containing multiple stereogenic atoms but absolute stereochemistry is not implied. In this case, we report that generation of enantiomeric pairs is imperative in database development if one is to obtain accurate docking results. A study is described on the impact of the neglect of enantiomeric pairs on virtual screening using the human homolog of murine double minute 2 (MDM2) protein, the product of a proto-oncogene, as the target. Docking in MDM2 with GLIDE 4.0 was performed using the NCI Diversity Set 3D database and, for comparison, a set of enantiomers we created corresponding to mirror image structures of the single enantiomers of chiral compounds present in the NCI Diversity Set. Our results demonstrate that potential lead candidates may be overlooked when databases contain 3D structures representing only a single enantiomer of racemic chiral compounds.     

SHAFTS: a hybrid approach for 3D molecular similarity calculation. 1. Method and assessment of virtual screening

We developed a novel approach called SHAFTS (SHApe-FeaTure Similarity) for 3D molecular similarity calculation and ligand-based virtual screening. SHAFTS adopts a hybrid similarity metric combined with molecular shape and colored (labeled) chemistry groups annotated by pharmacophore features for 3D similarity calculation and ranking, which is designed to integrate the strength of pharmacophore matching and volumetric overlay approaches. A feature triplet hashing method is used for fast molecular alignment poses enumeration, and the optimal superposition between the target and the query molecules can be prioritized by calculating corresponding "hybrid similarities". SHAFTS is suitable for large-scale virtual screening with single or multiple bioactive compounds as the query "templates" regardless of whether corresponding experimentally determined conformations are available. Two public test sets (DUD and Jain's sets) including active and decoy molecules from a panel of useful drug targets were adopted to evaluate the virtual screening performance. SHAFTS outperformed several other widely used virtual screening methods in terms of enrichment of known active compounds as well as novel chemotypes, thereby indicating its robustness in hit compounds identification and potential of scaffold hopping in virtual screening.     

Pharmacophore alignment search tool: influence of scoring systems on text-based similarity searching

The text-based similarity searching method Pharmacophore Alignment Search Tool is grounded on pairwise comparisons of potential pharmacophoric points between a query and screening compounds. The underlying scoring matrix is of critical importance for successful virtual screening and hit retrieval from large compound libraries. Here, we compare three conceptually different computational methods for systematic deduction of scoring matrices: assignment-based, alignment-based, and stochastic optimization. All three methods resulted in optimized pharmacophore scoring matrices with significantly superior retrospective performance in comparison with simplistic scoring schemes. Computer-generated similarity matrices of pharmacophoric features turned out to agree well with a manually constructed matrix. We introduce the concept of position-specific scoring to text-based similarity searching so that knowledge about specific ligand-receptor binding patterns can be included and demonstrate its benefit for hit retrieval. The approach was also used for automated pharmacophore elucidation in agonists of peroxisome proliferator activated receptor gamma, successfully identifying key interactions for receptor activation.     

SARS-CoV protease inhibitors design using virtual screening method from natural products libraries

Two natural products databases, the marine natural products database (MNPD) and the traditional Chinese medicines database (TCMD), were used to find novel structures of potent SARS-CoV protease inhibitors through virtual screening. Before the procedure, the databases were filtered by Lipinski's ROF and Xu's extension rules. The results were analyzed by statistic methods to eliminate the bias in target-based database screening toward higher molecular weight compounds for enhancing the hit rate. Eighteen lead compounds were recommended by the screening procedure. They were useful for experimental scientists in prioritizing drug candidates and studying the interaction mechanism. The binding mechanism was also analyzed between the best screening compound and the SARS protein.     

A pseudo-ligand approach to virtual screening

A virtual screening method is presented that is grounded on a receptor-derived pharmacophore model termed "virtual ligand" or "pseudo-ligand". The model represents an idealized constellation of potential ligand sites that interact with residues of the binding pocket. For rapid virtual screening of compound libraries the potential pharmacophore points of the virtual ligand are encoded as an alignment-free correlation vector, avoiding spatial alignment of pharmacophore features between the pharmacophore query (i.e., the virtual ligand) and the candidate molecule. The method was successfully applied to retrieving factor Xa inhibitors from a Ugi three-component combinatorial library, and yielded high enrichment of actives in a retrospective search for cyclooxygenase-2 (COX-2) inhibitors. The approach provides a concept for "de-orphanizing" potential drug targets and identifying ligands for hitherto unexplored or allosteric binding pockets.     

基于虚拟筛选策略发现新型潜在JAK3小分子抑制剂

利用已知活性的分子采用基于配体的策略构建药效团模型,通过基于类药规则、药效团模型、多种精度的分子对接算法、MM/GBSA结合能预测以及ADMET筛选手段对含约250万个分子的数据库进行虚拟筛选。发现5种JAK3抑制剂的新型骨架,其中6个以1-苯基咪唑烷-2-酮为骨架的分子在与JAK3激酶的结合能以及分子的ADMET性质评价方面均表现优异,具备高JAK3抑制剂潜力,被认为是虚拟筛选的命中分子。     

Identification of novel malarial cysteine protease inhibitors using structure-based virtual screening of a focused cysteine protease inhibitor library

Malaria, in particular that caused by Plasmodium falciparum , is prevalent across the tropics, and its medicinal control is limited by widespread drug resistance. Cysteine proteases of P. falciparum , falcipain-2 (FP-2) and falcipain-3 (FP-3), are major hemoglobinases, validated as potential antimalarial drug targets. Structure-based virtual screening of a focused cysteine protease inhibitor library built with soft rather than hard electrophiles was performed against an X-ray crystal structure of FP-2 using the Glide docking program. An enrichment study was performed to select a suitable scoring function and to retrieve potential candidates against FP-2 from a large chemical database. Biological evaluation of 50 selected compounds identified 21 diverse nonpeptidic inhibitors of FP-2 with a hit rate of 42%. Atomic Fukui indices were used to predict the most electrophilic center and its electrophilicity in the identified hits. Comparison of predicted electrophilicity of electrophiles in identified hits with those in known irreversible inhibitors suggested the soft-nature of electrophiles in the selected target compounds. The present study highlights the importance of focused libraries and enrichment studies in structure-based virtual screening. In addition, few compounds were screened against homologous human cysteine proteases for selectivity analysis. Further evaluation of structure-activity relationships around these nonpeptidic scaffolds could help in the development of selective leads for antimalarial chemotherapy.     

Pharmacophore modeling coupled with molecular dynamic simulation approach to identify new leads for meprin-β metalloprotease

Human meprin beta metalloprotease, a small subgroup of the astacin family, is a potent drug target for the treatment of several disorders such as fibrosis, neurodegenerative disease in particular Alzheimer and inflammatory bowel diseases. In this study, a ligand-based pharmacophore approach has been used for the selection of potentially active compounds to understand the inhibitory activities of meprin-β by using the sulfonamide scaffold based inhibitors. Using this dataset, a pharmacophore model (Hypo1) was selected on the basis of a highest correlation coefficient (0.959), lowest total cost (105.89) and lowest root mean square deviation (1.31 Å) values. All the pharmacophore hypotheses generated from the candidate inhibitors comprised four features: two hydrogen-bond acceptor, one hydrogen-bond donor and one zinc binder feature. The best validated pharmacophore model (Hypo1) was used for virtual screening of compounds from several databases. The selective hit compounds were filtered by drug likeness property, acceptable ADMET profile, molecular docking and DFT study. Molecular dynamic simulations with the final 10 hit compounds revealed that a large number of non-covalent interactions were formed with the active site and specificity sub-pockets of the meprin beta metalloprotease. This study assists in the development of the new lead molecules as well as gives a better understanding of their interaction with meprin-β.     

Targets looking for drugs: a multistep computational protocol for the development of structure-based pharmacophores and their applications for hit discovery

Pharmacophoresthree-dimensional (3D) arrangements of essential features enabling a molecule to exert a particular biological effectconstitute a very useful tool in drug design both in hit discovery and hit-to-lead optimization process. Two basic approaches for pharmacophoric model generation can be used by chemists, depending on the availability or not of the target 3D structure. In view of the rapidly growing number of protein structures that are now available, receptor-based pharmacophore generation methods are becoming more and more used. Since most of them require the knowledge of the 3D structure of the ligand-target complex, they cannot be applied when no compounds targeting the binding site of interest are known. Here, a GRID-based procedure for the generation of receptor-based pharmacophores starting from the knowledge of the sole protein structure is described and successfully applied to address three different tasks in the field of medicinal chemistry.     

Comparative analysis of pharmacophore screening tools

The pharmacophore concept is of central importance in computer-aided drug design (CADD) mainly because of its successful application in medicinal chemistry and, in particular, high-throughput virtual screening (HTVS). The simplicity of the pharmacophore definition enables the complexity of molecular interactions between ligand and receptor to be reduced to a handful set of features. With many pharmacophore screening softwares available, it is of the utmost interest to explore the behavior of these tools when applied to different biological systems. In this work, we present a comparative analysis of eight pharmacophore screening algorithms (Catalyst, Unity, LigandScout, Phase, Pharao, MOE, Pharmer, and POT) for their use in typical HTVS campaigns against four different biological targets by using default settings. The results herein presented show how the performance of each pharmacophore screening tool might be specifically related to factors such as the characteristics of the binding pocket, the use of specific pharmacophore features, and the use of these techniques in specific steps/contexts of the drug discovery pipeline. Algorithms with rmsd-based scoring functions are able to predict more compound poses correctly as overlay-based scoring functions. However, the ratio of correctly predicted compound poses versus incorrectly predicted poses is better for overlay-based scoring functions that also ensure better performances in compound library enrichments. While the ensemble of these observations can be used to choose the most appropriate class of algorithm for specific virtual screening projects, we remarked that pharmacophore algorithms are often equally good, and in this respect, we also analyzed how pharmacophore algorithms can be combined together in order to increase the success of hit compound identification. This study provides a valuable benchmark set for further developments in the field of pharmacophore search algorithms, e.g., by using pose predictions and compound library enrichment criteria.     

A fast algorithm for searching for molecules containing a pharmacophore in very large virtual combinatorial libraries.

We present a new algorithm for identifying molecules that display a pharmacophore, or in general a structural motif, by efficiently constructing and screening huge virtual combinatorial libraries of diverse compounds. The uniqueness of this algorithm is its ability to build and screen libraries of ca. 10(18) 3D molecular conformations within a reasonable time scale, thereby increasing the chemical space that can be virtually screened by many orders of magnitude. The algorithm may be used to design new molecules that display a desired pharmacophore on predefined sets of chemical scaffolds. This is demonstrated herein by screening a library of backbone cyclic peptides to find candidate peptido- and proteinomimetics.     

In silico drug repositioning of FDA-approved drugs to predict new inhibitors for alpha-synuclein aggregation

One of the hallmarks of Parkinson’s disease (PD), a long-term neurodegenerative syndrome, is the accumulation of alpha-synuclein (α-syn) fibrils. Despite numerous studies and efforts, inhibition of α-syn protein aggregation is still a challenge. To overcome this issue, we propose an in silico pharmacophore-based repositioning strategy, to find a pharmaceutical drug that, in addition to their defined role, can be used to prevent aggregation of the α-syn protein. Ligand-based pharmacophore modeling was developed and the best model was selected with validation parameters including 72 % sensitivity, 98 % specificity and goodness score about 0.7. The optimal model has three groups of hydrogen bond donor (HBD), three groups of hydrogen bond acceptor (HBA), and two aromatic rings (AR). The FDA-Approved reports in the ZINC15 database were screened with the pharmacophore model taken from inhibitor compounds. The model identified 22 hits, as promising candidate drugs for Parkinson's therapy. It is noteworthy that among these, 10 drugs have been reported to inhibition of α-syn aggregation or treat/reduce Parkinson's pathogenesis. This model was used to virtual screen ZINC, NCI databases, and natural products from the pomegranate. The results of this screen were filtered for their inability to cross the blood-brain barrier, poor oral bioavailability, etc. Finally, the selected compounds of two ZINC and NCI databases were combined and structurally clustered. Remained compounds were clustered in 28 different clusters, and the 17 compounds were introduced as final candidates.