On the value of homology models for virtual screening: discovering hCXCR3 antagonists by pharmacophore-based and structure-based approaches

Human chemokine receptor CXCR3 (hCXCR3) antagonists have potential therapeutic applications as antivirus, antitumor, and anti-inflammatory agents. A novel virtual screening protocol, which combines pharmacophore-based and structure-based approaches, was proposed. A three-dimensional QSAR pharmacophore model and a structure-based docking model were built to virtually screen for hCXCR3 antagonists. The hCXCR3 antagonist binding site was constructed by homology modeling and molecular dynamics (MD) simulation. By combining the structure-based and ligand-based screenings results, 95% of the compounds satisfied either pharmacophore or docking score criteria and would be chosen as hits if the union of the two searches was taken. The false negative rates were 15% for the pharmacophore model, 14% for the homology model, and 5% for the combined model. Therefore, the consistency of the pharmacophore model and the structural binding model is 219/273 = 80%. The hit rate for the virtual screening protocol is 273/286 = 95%. This work demonstrated that the quality of both the pharmacophore model and homology model can be measured by the consistency of the two models, and the false negatives in virtual screening can be reduced by combining two virtual screening approaches.     

A virtual screen for diverse ligands: discovery of selective G protein-coupled receptor antagonists

Virtual screening has become a major focus of bioactive small molecule lead identification, and reports of agonists and antagonists discovered via virtual methods are becoming more frequent. G protein-coupled receptors (GPCRs) are the one class of protein targets for which success with this approach has been limited. This is likely due to the paucity of detailed experimental information describing GPCR structure and the intrinsic function-associated structural flexibility of GPCRs which present major challenges in the application of receptor-based virtual screening. Here we describe an in silico methodology that diminishes the effects of structural uncertainty, allowing for more inclusive representation of a potential docking interaction with exogenous ligands. Using this approach, we screened one million compounds from a virtual database, and a diverse subgroup of 100 compounds was selected, leading to experimental identification of five structurally diverse antagonists of the thyrotropin-releasing hormone receptors (TRH-R1 and TRH-R2). The chirality of the most potent chemotype was demonstrated to be important in its binding affinity to TRH receptors; the most potent stereoisomer was noted to have a 13-fold selectivity for TRH-R1 over TRH-R2. A comprehensive mutational analysis of key amino acid residues that form the putative binding pocket of TRH receptors further verified the binding modality of these small molecule antagonists. The described virtual screening approach may prove applicable in the search for novel small molecule agonists and antagonists of other GPCRs.     

Binding mode similarity measures for ranking of docking poses: a case study on the adenosine A<Subscript>2A</Subscript> receptor

We report an investigation designed to explore alternative approaches for ranking of docking poses in the search for antagonists of the adenosine A_2A receptor, an attractive target for structure-based virtual screening. Calculation of 3D similarity of docking poses to crystallographic ligand(s) as well as similarity of receptor–ligand interaction patterns was consistently superior to conventional scoring functions for prioritizing antagonists over decoys. Moreover, the use of crystallographic antagonists and agonists, a core fragment of an antagonist, and a model of an agonist placed into the binding site of an antagonist-bound form of the receptor resulted in a significant early enrichment of antagonists in compound rankings. Taken together, these findings showed that the use of binding modes of agonists and/or antagonists, even if they were only approximate, for similarity assessment of docking poses or comparison of interaction patterns increased the odds of identifying new active compounds over conventional scoring.     

Identification of New Lead Molecules Against UBE2NL Enzyme for Cancer Therapy

Cancer is characterized by abnormal growth of cells. Targeting ubiquitin proteins in the discovery of new anticancer therapeutics is an attractive strategy. The present study uses the structure-based drug discovery methods to identify new lead structures, which are selective to the putative ubiquitin-conjugating enzyme E2N-like (UBE2NL). The 3D structure of the UBE2NL was evaluated using homology modeling techniques. The model was validated using standard in silico methods. The hydrophobic pocket of UBE2NL that aids in binding with its natural receptor ubiquitin-conjugating enzyme E2 variant (UBE2V) was identified through protein-protein docking study. The binding site region of the UBE2NL was identified using active site prediction tools. The binding site of UBE2NL which is responsible for cancer cell progression is considered for docking study. Virtual screening study with the small molecular structural database was carried out against the active site of UBE2NL. The ligand molecules that have shown affinity towards UBE2NL were considered for ADME prediction studies. The ligand molecules that obey the Lipinski’s rule of five and Jorgensen’s rule of three pharmacokinetic properties like human oral absorption etc. are prioritized. The resultant ligand molecules can be considered for the development of potent UBE2NL enzyme inhibitors for cancer therapy.     

A_1腺苷受体的同源模建及其结构验证

采用同源模建的方法构建了A1腺苷受体的三维结构,并与拮抗剂分子DPCPX对接,将得到的复合物结构进行5 ns的分子动力学模拟,以最后2 ns的平均结构和平衡后抽取的11帧构象共12个蛋白结构为研究对象,用包含52个活性分子和1000个诱饵分子的测试库,分别通过DOCK、VINA和GOLD三种对接软件进行评价,最终得出合理的蛋白质模型.根据top10%的富集因子(EF)和ROC曲线下面积(AU-ROC)的计算结果,我们认为GOLD是最适合A1腺苷受体的对接软件,而12个蛋白质结构中F5和Favg的三维结构模型比较合理,可以作为进一步大规模虚拟筛选的模型.     

Identification of potential Tpx inhibitors against pathogen-host interactions

Yersinia organisms cause many infectious diseases by invading human cells and delivering their virulence factors via the type three secretion system (T3SS). One alternative strategy in the fight against these pathogenic organisms is to interfere with their T3SS. Previous studies demonstrated that thiol peroxidase, Tpx is functional in the assembly of T3SS and its inhibition by salicylidene acylhydrazides prevents the secretion of pathogenic effectors. In this study, the aim was to identify potential inhibitors of Tpx using an integrated approach starting with high throughput virtual screening and ending with molecular dynamics simulations of selected ligands. Virtual screening of ZINC database of 500,000 compounds via ligand-based and structure-based pharmacophore models retrieved 10,000 hits. The structure-based pharmacophore model was validated using high-throughput virtual screening (HTVS). After multistep docking (SP and XP), common scaffolds were used to find common substructures and the ligand binding poses were optimized using induced fit docking. The stability of the protein–ligand complex was examined with molecular dynamics simulations and the binding free energy of the complex was calculated. As a final outcome eight compounds with different chemotypes were proposed as potential inhibitors for Tpx. The eight ligands identified by a detailed virtual screening protocol can serve as leads in future drug design efforts against the destructive actions of pathogenic bacteria.     

An automated system for the analysis of G protein-coupled receptor transmembrane binding pockets: alignment, receptor-based pharmacophores, and their application

G protein-coupled receptors (GPCRs) share a common architecture consisting of seven transmembrane (TM) domains. Various lines of evidence suggest that this fold provides a generic binding pocket within the TM region for hosting agonists, antagonists, and allosteric modulators. Here, a comprehensive and automated method allowing fast analysis and comparison of these putative binding pockets across the entire GPCR family is presented. The method relies on a robust alignment algorithm based on conservation indices, focusing on pharmacophore-like relationships between amino acids. Analysis of conservation patterns across the GPCR family and alignment to the rhodopsin X-ray structure allows the extraction of the amino acids lining the TM binding pocket in a so-called ligand binding pocket vector (LPV). In a second step, LPVs are translated to simple 3D receptor pharmacophore models, where each amino acid is represented by a single spherical pharmacophore feature and all atomic detail is omitted. Applications of the method include the assessment of selectivity issues, support of mutagenesis studies, and the derivation of rules for focused screening to identify chemical starting points in early drug discovery projects. Because of the coarseness of this 3D receptor pharmacophore model, however, meaningful scoring and ranking procedures of large sets of molecules are not justified. The LPV analysis of the trace amine-associated receptor family and its experimental validation is discussed as an example. The value of the 3D receptor model is demonstrated for a class C GPCR family, the metabotropic glutamate receptors.     

Structure-based virtual screening with supervised consensus scoring: evaluation of pose prediction and enrichment factors

Since the evaluation of ligand conformations is a crucial aspect of structure-based virtual screening, scoring functions play significant roles in it. However, it is known that a scoring function does not always work well for all target proteins. When one cannot know which scoring function works best against a target protein a priori, there is no standard scoring method to know it even if 3D structure of a target protein-ligand complex is available. Therefore, development of the method to achieve high enrichments from given scoring functions and 3D structure of protein-ligand complex is a crucial and challenging task. To address this problem, we applied SCS (supervised consensus scoring), which employs a rough linear correlation between the binding free energy and the root-mean-square deviation (rmsd) of a native ligand conformations and incorporates protein-ligand binding process with docked ligand conformations using supervised learning, to virtual screening. We evaluated both the docking poses and enrichments of SCS and five scoring functions (F-Score, G-Score, D-Score, ChemScore, and PMF) for three different target proteins: thymidine kinase (TK), thrombin (thrombin), and peroxisome proliferator-activated receptor gamma (PPARgamma). Our enrichment studies show that SCS is competitive or superior to a best single scoring function at the top ranks of screened database. We found that the enrichments of SCS could be limited by a best scoring function, because SCS is obtained on the basis of the five individual scoring functions. Therefore, it is concluded that SCS works very successfully from our results. Moreover, from docking pose analysis, we revealed the connection between enrichment and average centroid distance of top-scored docking poses. Since SCS requires only one 3D structure of protein-ligand complex, SCS will be useful for identifying new ligands.     

Pharmacophore-based virtual screening,biological evaluation and binding mode analysis of a novel protease-activated receptor 2 antagonist

Protease-activated receptor 2 (PAR₂) is a G protein-coupled receptor, mediating inflammation and pain signaling in neurons, thus it is considered to be a potential therapeutic target for inflammatory diseases. In this study, we performed a ligand-based virtual screening of 1.6 million compounds by employing a common-feature pharmacophore model and two-dimensional similarity search to identify a new PAR₂ antagonist. The common-feature pharmacophore model was established based on the biological screening results of our in-house library. The initial virtual screening yielded a total number of 47 hits, and additional biological activity tests including PAR₂ antagonism and anti-inflammatory effects resulted in a promising candidate, compound 43, which demonstrated an IC₅₀ value of 8.22 µM against PAR₂. In next step, a PAR₂ homology model was constructed using the crystal structure of the PAR₁ as a template to explore the binding mode of the identified ligands. A molecular docking method was optimized by comparing the binding modes of a known PAR₂ agonist GB110 and antagonist GB83, and applied to predict the binding mode of our hit compound 43. In-depth docking analyses revealed that the hydrophobic interaction with Phe243^5.39 is crucial for PAR₂ ligands to exert antagonistic activity. MD simulation results supported the predicted docking poses that PAR₂ antagonist blocked a conformational rearrangement of Na⁺ allosteric site in contrast to PAR₂ agonist that showed Na⁺ relocation upon GPCR activation. In conclusion, we identified new a PAR₂ antagonist together with its binding mode, which provides useful insights for the design and development of PAR₂ ligands.     

Identifying conformational changes of the β2 adrenoceptor that enable accurate prediction of ligand/receptor interactions and screening for GPCR modulators

The new β₂ Adrenoceptor (β₂AR) crystal structures provide a high-resolution snapshot of receptor interactions with two particular partial inverse agonists, (−)-carazolol and timolol. However, both experimental and computational studies of GPCR structure are significantly complicated by the existence of multiple conformational states coupled to ligand type and receptor activity. Agonists and antagonists induce or stabilize distinct changes in receptor structure that mediate a range of pharmacological activities. In this work, we (1) established that the existing β₂AR crystallographic conformers can be extended to describe ligand/receptor interactions for additional antagonist types, (2) generated agonist-bound receptor conformations, and (3) validated these models for agonist and antagonist virtual ligand screening (VLS). Using a ligand directed refinement protocol, we derived a single agonist-bound receptor conformation that selectively retrieved a diverse set of full and partial β₂AR agonists in VLS trials. Additionally, the impact of extracellular loop two conformation on VLS was assessed by docking studies with rhodopsin-based β₂AR homology models, and loop-deleted receptor models. A general strategy for constructing and selecting agonist-bound receptor pocket conformations is presented, which may prove broadly useful in creating agonist and antagonist bound models for other GPCRs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structure-based virtual screening of small-molecule antagonists of platelet integrin αIIbβ3 that do not prime the receptor to bind ligand

Integrin αIIbβ3 has emerged as an important therapeutic target for thrombotic vascular diseases owing to its pivotal role in mediating platelet aggregation through interaction with adhesive ligands. In the search for effective anti-thrombotic agents that can be administered orally without inducing the high-affinity ligand binding state, we recently discovered via high-throughput screening of 33,264 compounds a novel, αIIbβ3-selective inhibitor (RUC-1) of adenosine-5′-diphosphate (ADP) -induced platelet aggregation that exhibits a different chemical scaffold and mode of binding with respect to classical Arg-Gly-Asp (RGD)-mimicking αIIbβ3 antagonists. Most importantly, RUC-1 and its higher-affinity derivative, RUC-2, do not induce major conformational changes in the protein β3 subunit or prime the receptor to bind ligand. To identify additional αIIbβ3-selective chemotypes that inhibit platelet aggregation through similar mechanisms, we screened in silico over 2.5 million commercially available, ‘lead-like’ small molecules based on complementarity to the predicted binding mode of RUC-2 into the RUC-1-αIIbβ3 crystal structure. This first reported structure-based virtual screening application to the αIIbβ3 integrin led to the identification of 2 αIIbβ3-selective antagonists out of 4 tested, which compares favorably with the 0.003?% “hit rate” of our previous high-throughput chemical screening study. The newly identified compounds, like RUC-1 and RUC-2, showed specificity for αIIbβ3 compared to αVβ3 and did not prime the receptor to bind ligand. They thus may hold promise as αIIbβ3 antagonist therapeutic scaffolds.     

Representing receptor flexibility in ligand docking through relevant normal modes

Inspired by the current representation of the ligand-receptor binding process, a normal-mode-based methodology is presented to incorporate receptor flexibility in ligand docking and virtual screening. However, the systematic representation of the deformation space grows geometrically with the number of modes, and furthermore, midscale loop rearrangements like those found in protein kinase binding pockets cannot be accounted for with the first lowest-frequency modes. We thus introduced a measure of relevance of normal modes on a given region of interest and showed that only very few modes in the low-frequency range are necessary and sufficient to describe loop flexibility in cAMP-dependent protein kinase. We used this approach to generate an ensemble of representative receptor backbone conformations by perturbing the structure along a combination of relevant modes. Each ensemble conformation is complexed with known non-native binders to optimize the position of the binding-pocket side chains through a full flexible docking procedure. The multiple receptor conformations thus obtained are used in a small-scale virtual screening using receptor ensemble docking. We evaluated this algorithm on holo and apo structures of cAMP-dependent protein kinase that exhibit backbone rearrangements on two independent loop regions close to the binding pocket. Docking accuracy is improved, since the ligands considered in the virtual screening docked within 1.5 A to at least one of the structures. The discrimination between binders and nonbinders is also enhanced, as shown by the improvement of the enrichment factor. This constitutes a new step toward the systematic integration of flexible ligand-flexible receptor docking tools in structure-based drug discovery.     

Ligand-guided optimization of CXCR4 homology models for virtual screening using a multiple chemotype approach

CXCR4 is a G-protein coupled receptor for CXCL12 that plays an important role in human immunodeficiency virus infection, cancer growth and metastasization, immune cell trafficking and WHIM syndrome. In the absence of an X-ray crystal structure, theoretical modeling of the CXCR4 receptor remains an important tool for structure–function analysis and to guide the discovery of new antagonists with potential clinical use. In this study, the combination of experimental data and molecular modeling approaches allowed the development of optimized ligand-receptor models useful for elucidation of the molecular determinants of small molecule binding and functional antagonism. The ligand-guided homology modeling approach used in this study explicitly re-shaped the CXCR4 binding pocket in order to improve discrimination between known CXCR4 antagonists and random decoys. Refinement based on multiple test-sets with small compounds from single chemotypes provided the best early enrichment performance. These results provide an important tool for structure-based drug design and virtual ligand screening of new CXCR4 antagonists.     

A promising in silico protocol to develop novel PPARγ antagonists as potential anticancer agents: Design,synthesis and experimental validation via PPARγ protein activity and competitive binding assay

Peroxisome proliferator-activated receptor gamma (PPARγ), a member of the nuclear receptor superfamily is an excellent example of targets that orchestrates cancer, inflammation, lipid and glucose metabolism. We report a protocol for the development of novel PPARγ antagonists by employing 3D QSAR based virtual screening for the identification of ligands with anticancer properties. The models are generated based on a large and diverse set of PPARγ antagonist ligands by the HYPOGEN algorithm using Discovery Studio 2019 drug design software. Among the 10 hypotheses generated, Hypotheses 2 showed the highest correlation coefficient values of 0.95 with less RMS deviation of 1.193. Validation of the developed pharmacophore model was performed by Fischer’s randomization and screening against test and decoy set. The GH score or goodness score was found to be 0.81 indicating moderate to a good model. The selected pharmacophore model Hypo 2 was used as a query model for further screening of 11,145 compounds from the PubChem, sc-PDB structure database, and designed novel ligands. Based on fit values and ADMET filter, the final 10 compounds with the predicated activity of ≤ 3 nM were subjected for docking analysis. Docking analysis revealed the unique binding mode with hydrophobic amino acid that can cause destabilization of the H12 which is an important molecular mechanism to prove its antagonist action. Based on high CDocker scores, Cpd31 was synthesized, purified, analyzed and screened for PPARγ competitive binding by TR-FRET assay. The biochemical protein binding results matched the predicted results. Further, Cpd31 was screened against cancer cells and validated the results.     

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.     

Hot-spots-guided receptor-based pharmacophores (HS-Pharm): a knowledge-based approach to identify ligand-anchoring atoms in protein cavities and prioritize structure-based pharmacophores

The design of biologically active compounds from ligand-free protein structures using a structure-based approach is still a major challenge. In this paper, we present a fast knowledge-based approach (HS-Pharm) that allows the prioritization of cavity atoms that should be targeted for ligand binding, by training machine learning algorithms with atom-based fingerprints of known ligand-binding pockets. The knowledge of hot spots for ligand binding is here used for focusing structure-based pharmacophore models. Three targets of pharmacological interest (neuraminidase, beta2 adrenergic receptor, and cyclooxygenase-2) were used to test the evaluated methodology, and the derived structure-based pharmacophores were used in retrospective virtual screening studies. The current study shows that structure-based pharmacophore screening is a powerful technique for the fast identification of potential hits in a chemical library, and that it is a valid alternative to virtual screening by molecular docking.     

Identification of Novel TRPC5 Inhibitors by Pharmacophore-Based and Structure-Based Approaches

Canonical transient receptor potential-5 (TRPC5), which belongs to the subfamily of transient receptor potential (TRP) channels, is a non-selective cation channel mainly expressed in the central nervous system and shows more restricted expression in the periphery. TRPC5 plays a crucial role in human physiology and pathology, for instance, anxiety, depression, epilepsy, pain, memory and chronic kidney disease (CKD). However, due to lack of the effective and selective inhibitors, its physiological and pathological mechanism remains so far unknown. It is therefore pivotal to identify potential TRPC5 inhibitors. We have applied ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS) methods. The pharmacophore models of TRPC5 antagonists generated by using the HypoGen and HipHop algorithms were used as a query model for the screening of potential inhibitors against the Specs database. The resultant hits from LBVS were further screened by SBVS. SBVS was carried out based on the homology model generation of human TRPC5, binding site identification, molecular dynamics optimization and molecular docking studies. In our systematic screening approaches, we have identified 7 hits compounds with comparable dock score after Lipinski and Veber rules, ADMET, PAINS analysis, cluster analysis, and similarity analysis. In conclusion, the current research provides novel backbones for the new-generation of TRPC5 inhibitors.     

Virtual screening of organic molecule databases. Design of focused libraries of potential ligands of NMDA and AMPA receptors

A system of virtual screening of organic molecule databases is designed, which permits preprocessing of databases, molecular docking to a three-dimensional model of receptor, and post-processing of the results obtained. Using this screening system, it is possible to reproduce positions of the known ligands in the glutamate sites of the NMDA and AMPA receptors and in the glycine site of the NMDA receptor, to substantially enrich the database with potentially active compounds, and to distinguish between the agonistic and antagonistic character of the action of these compounds in the case of docking to the open and closed forms of the binding sites. Based on the results of screening of a database of low-molecular-weight organic compounds (total of 135,000 structures) using models of the open and closed forms of the glutamate and glycine sites of the NMDA receptor and of the glutamate site of the AMPA receptor, focused libraries of potential agonists and antagonists of these sites were designed.