Discovery of high affinity ligands for β2-adrenergic receptor through pharmacophore-based high-throughput virtual screening and docking

Novel high affinity compounds for human β₂-adrenergic receptor (β₂-AR) were searched among the clean drug-like subset of ZINC database consisting of 9,928,465 molecules that satisfy the Lipinski's rule of five. The screening protocol consisted of a high-throughput pharmacophore screening followed by an extensive amount of docking and rescoring. The pharmacophore model was composed of key features shared by all five inactive states of β₂-AR in complex with inverse agonists and antagonists. To test the discriminatory power of the pharmacophore model, a small-scale screening was initially performed on a database consisting of 117 compounds of which 53 antagonists were taken as active inhibitors and 64 agonists as inactive inhibitors. Accordingly, 7.3% of the ZINC database subset (729,413 compounds) satisfied the pharmacophore requirements, along with 44 antagonists and 17 agonists. Afterwards, all these hit compounds were docked to the inactive apo form of the receptor using various docking and scoring protocols. Following each docking experiment, the best pose was further evaluated based on the existence of key residues for antagonist binding in its vicinity. After final evaluations based on the human intestinal absorption (HIA) and the blood brain barrier (BBB) penetration properties, 62 hit compounds have been clustered based on their structural similarity and as a result four scaffolds were revealed. Two of these scaffolds were also observed in three high affinity compounds with experimentally known K_i values. Moreover, novel chemical compounds with distinct structures have been determined as potential β₂-AR drug candidates.     

Pharmacophore modeling,virtual screening,docking and in silico ADMET analysis of protein kinase B (PKB β) inhibitors

Protein kinase B (PKB) is a key mediator of proliferation and survival pathways that are critical for cancer growth. Therefore, inhibitors of PKB are useful agents for the treatment of cancer. Herein, we describe pharmacophore-based virtual screening combined with docking study as a rational strategy for identification of novel hits or leads. Pharmacophore models of PKB β inhibitors were established using the DISCOtech and refined with GASP from compounds with IC₅₀ values ranging from 2.2 to 246 nM. The best pharmacophore model consists of one hydrogen bond acceptor (HBA), one hydrogen bond donor (HBD) site and two hydrophobic (HY) features. The pharmacophore models were validated through receiver operating characteristic (ROC) and Güner-Henry (GH) scoring methods indicated that the model-3 was statistically valuable and reliable in identifying PKB β inhibitors. Pharmacophore model as a 3D search query was searched against NCI database. Several compounds with different structures (scaffolds) were retrieved as hits. Molecules with a Q_fit value of more than 95 and three other known inhibitors were docked in the active site of PKB to further explore the binding mode of these compounds. Finally in silico pharmacokinetic and toxicities were predicted for active hit molecules. The hits reported here showed good potential to be PKB β inhibitors.     

In silico screening and study of novel ERK2 inhibitors using 3D QSAR,docking and molecular dynamics

ERK2 is a dual specificity protein kinase, part of the Ras/Raf/MEK/ERK signal transduction cascade. It forms an interesting target for inhibition based on its relationship with cell proliferation and oncogenesis. A 3D QSAR pharmacophore model (Hypo1) with high correlation (r = 0.938) was developed for ERK2 ATP site on the basis of experimentally known inhibitors. The model included three hydrogen bonds, and one hydrophobic site. Assessment of Hypo1 through Fisher randomization, cost analysis, leave one out method and decoy test suggested that the model can reliably detect ERK2 inhibitors. Hypo1 has been used for virtual screening of potential inhibitors from ZINC, Drug Bank, NCI, Maybridge and Chembank databases. Using Hypo1 as a query, databases have been interrogated for compounds who meet the pharmacophore features. The resulting hit compounds were subject to docking and analysis. Docking and molecular dynamics analysis showed that in order to achieve a higher potency compounds have to interact with catalytic site, glycine rich loop, Hinge region, Gatekeeper region and ATP site entrance residues. We also identified catalytic site and Glycine rich loop as important regions to bind by molecules for better potency and selectivity.     

A systematic methodology for large scale compound screening: A case study on the discovery of novel S1PL inhibitors

Decrease in sphingosine 1-phosphate (S1P) concentration induces migration of pathogenic T cells to the blood stream, disrupts the CNS and it is implicated in multiple sclerosis (MS), a progressive inflammatory disorder of the central nervous system (CNS), and Alzheimer’s disease (AD). A promising treatment alternative for MS and AD is inhibition of the activity of the microsomal enzyme sphingosine 1-phosphate lyase (S1PL), which degrades intracellular S1P. This report describes an integrated systematic approach comprising virtual screening, molecular docking, substructure search and molecular dynamics simulation to discover novel S1PL inhibitors. Virtual screening of the ZINC database via ligand-based and structure-based pharmacophore models yielded 10000 hits. After molecular docking, common substructures of the top ranking hits were identified. The ligand binding poses were optimized by induced fit docking. MD simulations were performed on the complex structures to determine the stability of the S1PL-ligand complex and to calculate the binding free energy. Selectivity of the selected molecules was examined by docking them to hERG and cytochrome P450 receptors. As a final outcome, 15 compounds from different chemotypes were proposed as potential S1PL inhibitors. These molecules may guide future medicinal chemistry efforts in the discovery of new compounds against the destructive action of pathogenic T cells.     

Strategies for generating less toxic P-selectin inhibitors: pharmacophore modeling, virtual screening and counter pharmacophore screening to remove toxic hits

This paper describes the generation of ligand-based as well as structure-based models and virtual screening of less toxic P-selectin receptor inhibitors. Ligand-based model, 3D-pharmacophore was generated using 27 quinoline salicylic acid compounds and is used to retrieve the actives of P-selectin. This model contains three hydrogen bond acceptors (HBA), two ring aromatics (RA) and one hydrophobic feature (HY). To remove the toxic hits from the screened molecules, a counter pharmacophore model was generated using inhibitors of dihydrooratate dehydrogenase (DHOD), an important enzyme involved in nucleic acid synthesis, whose inhibition leads to toxic effects. Structure-based models were generated by docking and scoring of inhibitors against P-selectin receptor, to remove the false positives committed by pharmacophore screening. The combination of these ligand-based and structure-based virtual screening models were used to screen a commercial database containing 538,000 compounds.     

Combined in silico approaches for the identification of novel inhibitors of human islet amyloid polypeptide (hIAPP) fibrillation

Human islet amyloid polypeptide (hIAPP) is a natively unfolded polypeptide hormone of glucose metabolism, which is co-secreted with insulin by the β-cells of the pancreas. In patients with type 2 diabetes, IAPP forms amyloid fibrils because of diabetes-associated β-cells dysfunction and increasing fibrillation, in turn, lead to failure of secretory function of β-cells. This provides a target for the discovery of small organic molecules against protein aggregation diseases. However, the binding mechanism of these molecules with monomers, oligomers and fibrils to inhibit fibrillation is still an open question. In this work, ligand and structure-based in silico approaches were used to identify novel fibrillation inhibitors and/or fibril binding compounds. The best pharmacophore model was used as a 3D search query for virtual screening of a compound database to identify novel molecules having the potential to be therapeutic agents against protein aggregation diseases. Docking and molecular dynamics simulation studies were used to explore the interaction pattern and mechanism of the identified novel small molecules with predicted hIAPP structure, its aggregation prone conformation and fibril forming segments. We show that catechins with galloyl group and molecules having two to three planar apolar rings bind to hIAPP structures and fibril forming segments with greater affinity. The differences in binding affinities of different compounds against several fibril forming segments of the peptide suggest that a mixture of active compounds may be required for treatment of aggregation diseases.     

Pharmacophore identification and virtual screening for methionyl-tRNA synthetase inhibitors

Aminoacyl-tRNA synthetases (aaRSs) are essential enzymes involved in protein biosynthesis in all living organisms and are an unexploited antibacterial targets, as many strains of bacteria have become resistant to all established classes of antibiotics. Therefore, the main aim of this study is to discover new lead molecules which would be useful as anti-bacterial compounds. Pharmacophore models were developed by using CATALYST HypoGen with a training set of 29 diverse methionyl-tRNA synthetase (MetRS) inhibitors. The best quantitative pharmacophore hypothesis (Hypo1) obtained a correlation coefficient of 0.975, root mean square deviation (RMSD) of 0.55 and cost difference (null cost-total cost) of 70.32. This Hypo1 was validated by two methods, first by using 104 test set molecules which resulted a correlation of 0.926 between HypoGen estimated activities versus experimental activities and secondly by Cat-Scramble validation method. This validated pharmacophore model was further used for screening databases for discovery of new MetRS inhibitors. The new lead compounds were further analyzed for drug-like properties. Homology modeled structure of Staphylococcus aureus MetRS was built and molecular docking studies were performed with many inhibitors using the newly built protein structure. Finally, it was found that the new leads exhibited good estimated inhibitory activity, calculated binding properties similar to experimentally proven compounds and also favorable drug-like properties.     

流感病毒神经氨酸酶抑制剂的理性筛选

介绍了理性筛选流感病毒神经氨酸酶抑制剂的全过程，共分4个阶段：1)化合物数据库类药性处理；2)建立神经氨酸酶抑制剂三维药效团并对目标数据库进行构象搜索；3)分子对接及对接后分析；4)神经氨酸酶抑制模型的建立及待测化合物的活性检测。活性检测后发现4个活性化合物，其中Ic。为10。M的化合物1个．Ic，。为10^-6M的化合物2个，IC50为10^-7M的化合物1个。应用理性筛选方法，从化合物数据库中挑选出部分化合物进行神经氨酸酶抑制活性的筛选，减少了药物筛选的盲目性，提高了药物发现的机率。     

表皮生长因子受体酪氨酸激酶抑制剂的药效团研究

根据一系列表皮生长因子受体酪氨酸激酶抑制剂的三维定量构效关系研究,得到了该类抑制剂的药效团,研究结果与Novartis的药效团模型相当类似。药效团包括一个氢键受体,一个氢键给体,一个疏水区和一个带有氯或溴原子的苯环。该药效团对于研究表皮生长因子受体酪氨酸激酶抑制剂结构与活性的关系具有重要的意义。通过三维数据库搜索可能会得到新的先导化合物。     

In silico identification of novel inhibitors against Plasmodium falciparum dihydroorate dehydrogenase

Plasmodium falciparum causes the most fatal form of malaria and accounts for over 1 million deaths annually, yet currently used drug therapies are compromised by resistance. The malaria parasite cannot salvage pyrimidines and relies on de novo biosynthesis for survival. The enzyme dihydrooratate dehydrogenase (DHODH), a mitochondrial flavoenzyme, catalyzes the rate-limiting step of this pathway and is therefore an attractive anti-malarial chemotherapeutic target. In an effort to design new and potential anti-malarials, structure-based pharmacophore mapping, molecular docking, binding energy calculations and binding affinity predictions were employed in a virtual screening strategy to design new and potent P. falciparum dihydrooratate dehydrogenase (PfDHODH) inhibitors. A structure-based pharmacophore model was generated which consist of important interactions as observed in co-crystal of PfDHODH enzyme. The developed model was used to retrieve molecules from ChemBridge database, a freely available commercial database. A total of 87 molecules mapped on the modeled pharmacophore from the database. The retrieved hits were further screened by docking simulation, binding energy calculations and biding affinity predictions using genetic optimization for ligand docking (GOLD) and MOE. Based on these results, finally 26 chemo-types molecules were predicted as new, potential and structurally diverse PfDHODH inhibitors.     

Towards more accurate pharmacophore modeling: Multicomplex-based comprehensive pharmacophore map and most-frequent-feature pharmacophore model of CDK2

Pharmacophore modeling, including ligand- and structure-based approaches, has become an important tool in drug discovery. However, the ligand-based method often strongly depends on the training set selection, and the structure-based pharmacophore model is usually created based on apo structures or a single protein-ligand complex, which might miss some important information. In this study, multicomplex-based method has been suggested to generate a comprehensive pharmacophore map of cyclin-dependent kinase 2 (CDK2) based on a collection of 124 crystal structures of human CDK2-inhibitor complex. Our multicomplex-based comprehensive pharmacophore map contains almost all the chemical features important for CDK2-inhibitor interactions. A comparison with previously reported ligand-based pharmacophores has revealed that the ligand-based models are just a subset of our comprehensive map. Furthermore, one most-frequent-feature pharmacophore model consisting of the most frequent pharmacophore features was constructed based on the statistical frequency information provided by the comprehensive map. Validations to the most-frequent-feature model show that it can not only successfully discriminate between known CDK2 inhibitors and the molecules of focused inactive dataset, but also is capable of correctly predicting the activities of a wide variety of CDK2 inhibitors in an external active dataset. Obviously, this investigation provides some new ideas about how to develop a multicomplex-based pharmacophore model that can be used in virtual screening to discover novel potential lead compounds.     

乙酰胆碱酯酶抑制剂虚拟筛选方法研究

在虚拟筛选过程中，虚拟筛选策略和方法是获取结果的基础，但需要通过实验数据来检验。获得虚拟筛选合理的限制因素，将为大规模虚拟筛选提供方法依据。通过建立乙酰胆碱酯酶抑制剂的活性检测方法，检测了4000多个化合物的生物活性，并发现了34个活性较高的化合物，同时将设计的6个不同的虚拟筛选方法分别用于2个乙酰胆碱醋酶抑制剂虚拟筛选模型。将所预测的可能活性化合物及两模型共有的可能活性化合物分别与实验所得的活性化合物对照，综合分析讨论虚拟筛选乙酰胆碱酯酶抑制剂的重要因素和进一步富集活性化合物的方法，为大规模的虚拟筛选乙酰胆碱醋酶抑制剂提供可靠依据，并为其他基于蛋白酶结构的虚拟筛选提供参考。     

Identifying potential selective fluorescent probes for cancer-associated protein carbonic anhydrase IX using a computational approach

Carbonic anhydrase IX (CAIX) is a biomarker for tumor hypoxia. Fluorescent inhibitors of CAIX have been used to study hypoxic tumor cell lines. However, these inhibitor-based fluorescent probes may have a therapeutic effect that is not appropriate for monitoring treatment efficacy. In the search for novel fluorescent probes that are not based on known inhibitors, a database of 20,860 fluorescent compounds was virtually screened against CAIX using hierarchical virtual ligand screening (HierVLS). The screening database contained 14,862 compounds tagged with the ATTO680 fluorophore plus an additional 5998 intrinsically fluorescent compounds. Overall ranking of compounds to identify hit molecular probe candidates utilized a principal component analysis (PCA) approach. Four potential binding sites, including the catalytic site, were identified within the structure of the protein and targeted for virtual screening. Available sequence information for 23 carbonic anhydrase isoforms was used to prioritize the four sites based on the estimated “uniqueness” of each site in CAIX relative to the other isoforms. A database of 32 known inhibitors and 478 decoy compounds was used to validate the methodology. A receiver–operating characteristic (ROC) analysis using the first principal component (PC1) as predictive score for the validation database yielded an area under the curve (AUC) of 0.92. AUC is interpreted as the probability that a binder will have a better score than a non-binder. The use of first component analysis of binding energies for multiple sites is a novel approach for hit selection. The very high prediction power for this approach increases confidence in the outcome from the fluorescent library screening. Ten of the top scoring candidates for isoform-selective putative binding sites are suggested for future testing as fluorescent molecular probe candidates.     

First universal pharmacophore model for hERG1 K+ channel activators: acthER

The intra-cavitary drug blockade of hERG1 channel has been extensively studied, both experimentally and theoretically. Structurally diverse ligands inadvertently block the hERG1 K⁺ channel currents lead to drug induced Long QT Syndrome (LQTS). Accordingly, designing either hERG1 channel openers or current activators, with the potential to target other binding pockets of the channel, has been introduced as a viable approach in modern anti-arrhythmia drug development. However, reports and investigations on the molecular mechanisms underlying activators binding to the hERG1 channel remain sparse and the overall molecular design principles are largely unknown. Most of the hERG1 activators were discovered during mandatory screening for hERG1 blockade. To fill this apparent deficit, the first universal pharmacophore model for hERG1 K⁺ channel activators was developed using PHASE. 3D structures of 18 hERG1 K⁺ channel activators and their corresponding measured binding affinity values were used in the development of pharmacophore models. These compounds spanned a range of structurally different chemotypes with moderate variation in binding affinity. A five sites AAHRR (A, hydrogen-bond accepting, H, hydrophobic, R, aromatic) pharmacophore model has shown reasonable high statistical results compared to the other developed more than 1000 hypotheses. This model was used to construct steric and electrostatic contour maps. The predictive power of the model was tested with 3 external test set compounds as true unknowns. Finally, the pharmacophore model was combined with the previously developed receptor-based model of hERG1 K⁺ channel to develop and screen novel activators. The results are quite striking and it suggests a greater future role for pharmacophore modeling and virtual drug screening simulations in deciphering complex patterns of molecular mechanisms of hERG1 channel openers at the target sites. The developed model is available upon request and it may serve as basis for the synthesis of novel therapeutic hERG1 activators.     

The use of a pharmacophore model for identification of novel ligands for the benzodiazepine binding site of the GABAA receptor

A Catalyst pharmacophore model has been developed for the benzodiazepine site within the GABA(A) receptor complex. The model is based on a pharmacophore model originally proposed by Cook and co-workers (Drug Des. Discovery 1995, 12, 193-248) and further developed by Kahnberg et al. (J. Med. Chem. 2002, 45, 4188-4201). The Catalyst pharmacophore model has been validated by using a series of flavonoids with varying affinities for the benzodiazepine receptor and has then been used as a search query in database searching with the aim of finding novel structures which have the possibility to be modified into novel lead compounds. Five of the hits from the database searching were purchased and their affinities for the benzodiazepine site of the GABA(A) receptor were determined. Two of the compounds displayed K(i) values below 10 microM. The substance showing highest potency in-vitro displayed an affinity of 121 nM making it an interesting compound for optimization. The false positive compounds (K(i) values >10 microM affinities) have been analysed in terms of conformational energy penalties and possibilities for hydrogen bond interactions. The analysis clearly demonstrates the need for post processing of Catalyst hits.     

Structure based virtual screening-driven identification of monastrol as a potent urease inhibitor

Virtual screening uses computer based methods to discover new ligands on the basis of biological structures. Among all virtual screening methods structure based docking has received considerable attention. In an attempt to identify new ligands as urease inhibitors, structure-based virtual screening (SBVS) of an in-house database of 10,000 organic compounds was carried out. The X-ray crystallographic structure of Bacillus pasteurii (BP) in complex with acetohydroxamic acid (PDB Code 4UBP) was used as a protein structure. As a starting point, ~10,000 compounds of our in-house database were analyzed to check redundancy and the compounds found repeated were removed from the database. Finally 6993 compounds were docked into the active site of BP urease using GOLD and MOE-Dock software. A remarkable feature of this study was the identification of monastrol, a well-known KSP inhibitor already in clinical trials, as a novel urease inhibitor. The hits identified were further evaluated by molecular docking and on examination of the affinity predictions, twenty-seven analogs of monastrol were synthesized by a multicomponent Biginelli reaction followed by their in vitro screening as urease inhibitors. Finally twelve compounds were identified as new urease inhibitors. The excellent in vitro activity suggested that these compounds may serve as viable lead compounds for the treatment of urease related problems.     

Identification of novel PAD4 inhibitors based on a pharmacophore model derived from transition state coordinates

1.4 Protein arginine deiminases 4 (PAD4) is an attractive target for the development of novel and selective inhibitors of Rheumatoid Arthritis (RA). F-amidine is known as mechanism-based inhibitor targeting PAD4 and used as inactivators by covalently modifying the active site Cys645. To identify novel structural inhibitors of PAD4, we investigated the flexibility of protein on basis of the transition state geometry of PAD4 inhibited by F-amidine from our previous QM/MM calculation. And a pharmacophore model was generated containing four features (ADHH) using five representative structures from molecular dynamic (MD) simulation on basis of the transition state geometry of PAD4 inhibited by F-amidine. We performed virtual screening using the pharmacophore model and molecular docking methods, resulting in the discovery of two molecules with K_D (dissociation equilibrium constant) values of 112 μM and 218 μΜ against PAD4 through Surface Plasmon Resonance (SPR) experiments. These two molecules could potentially serve as PAD4 inhibitors.     

Integration of ligand and structure based approaches for identification of novel MbtI inhibitors in Mycobacterium tuberculosis and molecular dynamics simulation studies

Mycobacterium tuberculosis is an obligate pathogen of mammals and is responsible for more than two million deaths annually. The ability to acquire iron from the extracellular environment is a key determinant of pathogenicity in mycobacteria. M. tuberculosis acquires iron exclusively through the siderophores. Several lines of evidence suggest that siderophores have a critical role in bacterial growth and virulence. Hence, in the present study, we have used a combined ligand and structure-based drug design approach for identification of novel inhibitors against salicylate synthase MbtI, a unique and essential enzyme for the biosynthesis of siderophores in M. tuberculosis. We have generated the ligand based and structure based pharmacophores and validated exhaustively. From the validation results it was found that GH (Goodness of Hit) scores for the selected ligand based and structure based pharmacophore models were 0.89 and 0.97, respectively, which indicate that the quality of the pharmacophore models are acceptable as GH value is >0.7. The validated pharmacophores were used for screening the ZINC database. A total of 73 hits, obtained through various insilico screening techniques, were further enriched to 17 hits using docking studies. Molecular dynamics simulations were carried out to compare the binding mode and stability of complexes of MbtI bound with substrate, known inhibitors, and three top ranked hits. The results obtained in this study gave assurance about the identified hits as prospective inhibitors of MbtI.     

A support vector machines approach for virtual screening of active compounds of single and multiple mechanisms from large libraries at an improved hit-rate and enrichment factor

Support vector machines (SVM) and other machine-learning (ML) methods have been explored as ligand-based virtual screening (VS) tools for facilitating lead discovery. While exhibiting good hit selection performance, in screening large compound libraries, these methods tend to produce lower hit-rate than those of the best performing VS tools, partly because their training-sets contain limited spectrum of inactive compounds. We tested whether the performance of SVM can be improved by using training-sets of diverse inactive compounds. In retrospective database screening of active compounds of single mechanism (HIV protease inhibitors, DHFR inhibitors, dopamine antagonists) and multiple mechanisms (CNS active agents) from large libraries of 2.986 million compounds, the yields, hit-rates, and enrichment factors of our SVM models are 52.4–78.0%, 4.7–73.8%, and 214–10,543, respectively, compared to those of 62–95%, 0.65–35%, and 20–1200 by structure-based VS and 55–81%, 0.2–0.7%, and 110–795 by other ligand-based VS tools in screening libraries of ≥1 million compounds. The hit-rates are comparable and the enrichment factors are substantially better than the best results of other VS tools. 24.3–87.6% of the predicted hits are outside the known hit families. SVM appears to be potentially useful for facilitating lead discovery in VS of large compound libraries.     

Identification of potent inhibitors of DNA methyltransferase 1 (DNMT1) through a pharmacophore-based virtual screening approach

DNA methylation is an epigenetic change that results in the addition of a methyl group at the carbon-5 position of cytosine residues. DNA methyltransferase (DNMT) inhibitors can suppress tumour growth and have significant therapeutic value. However, the established inhibitors are limited in their application due to their substantial cytotoxicity. Additionally, the standard drugs for DNMT inhibition are non-selective cytosine analogues with considerable cytotoxic side-effects. In the present study, we have designed a workflow by integrating various ligand-based and structure-based approaches to discover new agents active against DNMT1. We have derived a pharmacophore model with the help of available DNMT1 inhibitors. Utilising this model, we performed the virtual screening of Maybridge chemical library and the identified hits were then subsequently filtered based on the Naïve Bayesian classification model. The molecules that have returned from this classification model were subjected to ensemble based docking. We have selected 10 molecules for the biological assay by inspecting the interactions portrayed by these molecules. Three out of the ten tested compounds have shown DNMT1 inhibitory activity. These compounds were also found to demonstrate potential inhibition of cellular proliferation in human breast cancer MDA-MB-231 cells. In the present study, we have utilized a multi-step virtual screening protocol to identify inhibitors of DNMT1, which offers a starting point to develop more potent DNMT1 inhibitors as anti-cancer agents.