Three-dimensional pharmacophore modeling of liver-X receptor agonists

High cholesterol levels contribute to hyperlipidemia. Liver X receptors (LXRs) are the drug targets. LXRs regulate the cholesterol absorption, biosynthesis, transportation, and metabolism. Novel agonists of LXR, especially LXRβ, are attractive solutions for treating hyperlipidemia. In order to discover novel LXRβ agonists, a three-dimensional pharmacophore model was built based upon known LXRβ agonists. The model was validated with a test set, a virtual screening experiment, and the FlexX docking approach. Results show that the model is capable of predicting a LXRβ agonist activity. Ligand-based virtual screening results can be refined by cross-linking by structure-based approaches. This is because two ligands that are mapped in the same way to the same pharmacophore model may have significantly different binding behaviors in the receptor's binding pocket. This paper reports our approach to identify reliable pharmacophore models through combining both ligand- and structure-based approaches.     

Hierarchical virtual screening: identification of potential high-affinity and selective β(3)-adrenergic receptor agonists

The hierarchical virtual screening (HVS) study, consisting of pharmacophore modelling, docking and VS of the generated focussed virtual library, has been carried out to identify novel high-affinity and selective β(3)-adrenergic receptor (β-AR) agonists. The best pharmacophore model, comprising one H-bond donor, two hydrophobes, one positive ionizable and one negative ionizable feature, was developed based on a training set of 51 β(3)-AR agonists using the pharmacophore generation protocol implemented in Discovery Studio. The model was further validated with the test set, external set and ability of the pharmacophoric features to complement the active site amino acids of the homology modelled β(3)-AR developed using MODELLER software. The focussed virtual library was generated using the structure-based insights gained from our earlier reported comprehensive study focussing on the structural basis of β-AR subtype selectivity of representative agonists and antagonists. The HVS with the sequential use of the best pharmacophore model and homology modelled β(3)-AR in the screening of the generated focussed library has led to the identification of potential virtual leads as novel high-affinity and selective β(3)-AR agonists.     

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.     

GPCR structure-based virtual screening approach for CB2 antagonist search

The potential for therapeutic specificity in regulating diseases has made cannabinoid (CB) receptors one of the most important G-protein-coupled receptor (GPCR) targets in search for new drugs. Considering the lack of related 3D experimental structures, we have established a structure-based virtual screening protocol to search for CB2 bioactive antagonists based on the 3D CB2 homology structure model. However, the existing homology-predicted 3D models often deviate from the native structure and therefore may incorrectly bias the in silico design. To overcome this problem, we have developed a 3D testing database query algorithm to examine the constructed 3D CB2 receptor structure model as well as the predicted binding pocket. In the present study, an antagonist-bound CB2 receptor complex model was initially generated using flexible docking simulation and then further optimized by molecular dynamic and mechanical (MD/MM) calculations. The refined 3D structural model of the CB2-ligand complex was then inspected by exploring the interactions between the receptor and ligands in order to predict the potential CB2 binding pocket for its antagonist. The ligand-receptor complex model and the predicted antagonist binding pockets were further processed and validated by FlexX-Pharm docking against a testing compound database that contains known antagonists. Furthermore, a consensus scoring (CScore) function algorithm was established to rank the binding interaction modes of a ligand on the CB2 receptor. Our results indicated that the known antagonists seeded in the testing database can be distinguished from a significant amount of randomly chosen molecules. Our studies demonstrated that the established GPCR structure-based virtual screening approach provided a new strategy with a high potential for in silico identifying novel CB2 antagonist leads based on the homology-generated 3D CB2 structure model.     

Effective virtual screening protocol for CYP2C9 ligands using a screening site constructed from flurbiprofen and S-warfarin pockets

An effective virtual screening protocol was developed against an extended active site of CYP2C9, which was derived from X-ray structures complexed with flubiprofen and S-warfarin. Virtual screening has been effectively supported by our structure-based pharmacophore model. Importance of hot residues identified by mutation data and structural analysis was first estimated in an enrichment study. Key role of Arg108 and Phe114 in ligand binding was also underlined. Our screening protocol successfully identified 76% of known CYP2C9 ligands in the top 1% of the ranked database resulting 76-fold enrichment relative to random situation. Relevance of the protocol was further confirmed in selectivity studies, when 89% of CYP2C9 ligands were retrieved from a mixture of CYP2C9 and CYP2C8 ligands, while only 22% of CYP2C8 ligands were found applying the structure-based pharmacophore constraints. Moderate discrimination of CYP2C9 ligands from CYP2C18 and CYP2C19 ligands could also be achieved extending the application domain of our virtual screening protocol for the entire CYP2C family. Our findings further demonstrate the existence of an active site comprising of at least two binding pockets and strengthens the need of involvement of protein flexibility in 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.     

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.     

Integrating pharmacophore mapping,virtual screening,density functional theory,molecular simulation towards the discovery of novel apolipoprotein (apoE ε4) inhibitors

AimAn integrated protocol of virtual screening involving molecular docking, pharmacophore probing, and simulations was established to identify small novel molecules targeting crucial residues involved in the variant apoE ε4 to mimic its behavior as apoE2 thereby eliminating the amyloid plaque accumulation and facilitating its clearance.Materials and MethodsAn excellent ligand-based and structure-based approach was made to identify common pharmacophoric features involving structure-based docking with respect to apoE ε4 leading to the development of apoE ε4 inhibitors possessing new scaffolds. An effort was made to design multiple-substituted triazine derivatives series bearing a novel scaffold. A structure-based pharmacophore mapping was developed to explore the binding sites of apoE ε4 which was taken into consideration. Subsequently, virtual screening, ADMET, DFT searches were at work to narrow down the proposed hits to be forwarded as a potential drug likes candidates. Further, the binding patterns of the best-proposed hits were studied and were forwarded for molecular dynamic simulations of 10 ns for its structural optimization.ResultsSelectivity profile for the most promising candidates was studied, revealing significantly C13 and C15 to be the most potent compounds. The proposed hits can be forwarded for further study against apoE ε4 involved in neurological disorder Alzheimer’s.     

Ensemble pharmacophore meets ensemble docking: a novel screening strategy for the identification of RIPK1 inhibitors

Programmed cell death has been a fascinating area of research since it throws new challenges and questions in spite of the tremendous ongoing research in this field. Recently, necroptosis, a programmed form of necrotic cell death, has been implicated in many diseases including neurological disorders. Receptor interacting serine/threonine protein kinase 1 (RIPK1) is an important regulatory protein involved in the necroptosis and inhibition of this protein is essential to stop necroptotic process and eventually cell death. Current structure-based virtual screening methods involve a wide range of strategies and recently, considering the multiple protein structures for pharmacophore extraction has been emphasized as a way to improve the outcome. However, using the pharmacophoric information completely during docking is very important. Further, in such methods, using the appropriate protein structures for docking is desirable. If not, potential compound hits, obtained through pharmacophore-based screening, may not have correct ranks and scores after docking. Therefore, a comprehensive integration of different ensemble methods is essential, which may provide better virtual screening results. In this study, dual ensemble screening, a novel computational strategy was used to identify diverse and potent inhibitors against RIPK1. All the pharmacophore features present in the binding site were captured using both the apo and holo protein structures and an ensemble pharmacophore was built by combining these features. This ensemble pharmacophore was employed in pharmacophore-based screening of ZINC database. The compound hits, thus obtained, were subjected to ensemble docking. The leads acquired through docking were further validated through feature evaluation and molecular dynamics simulation.     

The discovery of Kv1.5 blockers as a case study for the application of virtual screening approaches

Different virtual screening techniques are available as alternatives to high throughput screening. These different techniques have been rarely used together on the same target. We had the opportunity to do so in order to discover novel blockers of the voltage-dependent potassium channel Kv1.5, a potential target for the treatment of atrial fibrillation. Our corporate database was searched, using a protein-based pharmacophore, derived from a homology model, as query. As a result, 244 molecules were screened in vitro, 19 of them (7.8%) were found to be active. Five of them, belonging to five different chemical classes, exhibited IC50 values under 10 microM. The performance of this structure-based virtual screening protocol has been compared with those of similarity and ligand-based pharmacophore searches. The analysis of the results supports the conventional wisdom of using as many virtual screening techniques as possible in order to maximize the chance of finding as many chemotypes as possible.     

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

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

Development of purely structure-based pharmacophores for the topoisomerase I-DNA-ligand binding pocket

Purely structure-based pharmacophores (SBPs) are an alternative method to ligand-based approaches and have the advantage of describing the entire interaction capability of a binding pocket. Here, we present the development of SBPs for topoisomerase I, an anticancer target with an unusual ligand binding pocket consisting of protein and DNA atoms. Different approaches to cluster and select pharmacophore features are investigated, including hierarchical clustering and energy calculations. In addition, the performance of SBPs is evaluated retrospectively and compared to the performance of ligand- and complex-based pharmacophores. SBPs emerge as a valid method in virtual screening and a complementary approach to ligand-focussed methods. The study further reveals that the choice of pharmacophore feature clustering and selection methods has a large impact on the virtual screening hit lists. A prospective application of the SBPs in virtual screening reveals that they can be used successfully to identify novel topoisomerase inhibitors.     

Novel approach to structure-based pharmacophore search using computational geometry and shape matching techniques

Computationally efficient structure-based virtual screening methods have recently been reported that seek to find effective means to utilize experimental structure information without employing detailed molecular docking calculations. These tools can be coupled with efficient experimental screening technologies to improve the probability of identifying hits and leads for drug discovery research. Commercial software ROCS (rapid overlay of chemical structures) from Open Eye Scientific is such an example, which is a shape-based virtual screening method using the 3D structure of a ligand, typically from a bound X-ray costructure, as the query. We report here the development of a new structure-based pharmacophore search method (called Shape4) for virtual screening. This method adopts a variant of the ROCS shape technology and expands its use to work with an empty crystal structure. It employs a rigorous computational geometry method and a deterministic geometric casting algorithm to derive the negative image (i.e., pseudoligand) of a target binding site. Once the negative image (or pseudoligand) is generated, an efficient shape comparison algorithm in the commercial OE SHAPE Toolkit is adopted to compare and match small organic molecules with the shape of the pseudoligand. We report the detailed computational protocol and its computational validation using known biologically active compounds extracted from the WOMBAT database. Models derived for five selected targets were used to perform the virtual screening experiments to obtain the enrichment data for various virtual screening methods. It was found that our approach afforded similar or better enrichment ratios than other related methods, often with better diversity among the top ranking computational hits.     

人类腺苷受体A3亚型拮抗剂的构效关系分析

构建人类腺苷受体A3亚型药效团模型和三维蛋白结构模型用于作用模式研究.以18个来源于文献具有腺苷受体A3亚型拮抗活性的化合物作为训练集,使用HypoGen方法构建药效团模型.通过同源模建和分子动力学模拟构建了人类腺苷受体A3亚型的三维蛋白模型,并利用PROCHECK方法评估该模型的合理性,对所得的结构使用分子对接程序进行作用模式分析,药效团模型和同源模建结果相互匹配较好.使用新药效团模型对MDL药物数据库(MDDR)中包含的约120000个化合物进行虚拟筛选,得到了8个候选化合物,用于进一步的生物学评价和活性测定.本工作对于人类腺苷受体A3亚型拮抗剂的设计和抗哮喘药物的研发具有一定的理论指导和应用价值.     

Pharmacophore model-aided virtual screening combined with comparative molecular docking and molecular dynamics for identification of marine natural products as SARS-CoV-2 papain-like protease inhibitors

Targeting SARS-CoV-2 papain-like protease using inhibitors is a suitable approach for inhibition of virus replication and dysregulation of host anti-viral immunity. Engaging all five binding sites far from the catalytic site of PLpro is essential for developing a potent inhibitor. We developed and validated a structure-based pharmacophore model with 9 features of a potent PLpro inhibitor. The pharmacophore model-aided virtual screening of the comprehensive marine natural product database predicted 66 initial hits. This hit library was downsized by filtration through a molecular weight filter of ≤ 500 g/mol. The 50 resultant hits were screened by comparative molecular docking using AutoDock and AutoDock Vina. Comparative molecular docking enables benchmarking docking and relieves the disparities in the search and scoring functions of docking engines. Both docking engines retrieved 3 same compounds at different positions in the top 1 % rank, hence consensus scoring was applied, through which CMNPD28766, aspergillipeptide F emerged as the best PLpro inhibitor. Aspergillipeptide F topped the 50-hit library with a pharmacophore-fit score of 75.916. Favorable binding interactions were predicted between aspergillipeptide F and PLpro similar to the native ligand XR8-24. Aspergillipeptide F was able to engage all the 5 binding sites including the newly discovered BL2 groove, site V. Molecular dynamics for quantification of Cα-atom movements of PLpro after ligand binding indicated that it exhibits highly correlated domain movements contributing to the low free energy of binding and a stable conformation. Thus, aspergillipeptide F is a promising candidate for pharmaceutical and clinical development as a potent SARS-CoV-2 PLpro inhibitor.     

Energetic analysis of fragment docking and application to structure-based pharmacophore hypothesis generation

We have developed a method that uses energetic analysis of structure-based fragment docking to elucidate key features for molecular recognition. This hybrid ligand- and structure-based methodology uses an atomic breakdown of the energy terms from the Glide XP scoring function to locate key pharmacophoric features from the docked fragments. First, we show that Glide accurately docks fragments, producing a root mean squared deviation (RMSD) of <1.0 Å for the top scoring pose to the native crystal structure. We then describe fragment-specific docking settings developed to generate poses that explore every pocket of a binding site while maintaining the docking accuracy of the top scoring pose. Next, we describe how the energy terms from the Glide XP scoring function are mapped onto pharmacophore sites from the docked fragments in order to rank their importance for binding. Using this energetic analysis we show that the most energetically favorable pharmacophore sites are consistent with features from known tight binding compounds. Finally, we describe a method to use the energetically selected sites from fragment docking to develop a pharmacophore hypothesis that can be used in virtual database screening to retrieve diverse compounds. We find that this method produces viable hypotheses that are consistent with known active compounds. In addition to retrieving diverse compounds that are not biased by the co-crystallized ligand, the method is able to recover known active compounds from a database screen, with an average enrichment of 8.1 in the top 1% of the database.