Virtual Fragment Screening: Discovery of Histamine H3 Receptor Ligands Using Ligand-Based and Protein-Based Molecular Fingerprints

VirtuousPocketome: a computational tool for screening protein-ligand complexes to identify similar binding sites

Protein residues within binding pockets play a critical role in determining the range of ligands that can interact with a protein, influencing its structure and function. Identifying structural similarities in proteins offers valuable insights into their function and activation mechanisms, aiding in predicting protein-ligand interactions, anticipating off-target effects, and facilitating the development of therapeutic agents. Numerous computational methods assessing global or local similarity in protein cavities have emerged, but their utilization is impeded by complexity, impractical automation for amino acid pattern searches, and an inability to evaluate the dynamics of scrutinized protein-ligand systems. Here, we present a general, automatic and unbiased computational pipeline, named VirtuousPocketome, aimed at screening huge databases of proteins for similar binding pockets starting from an interested protein-ligand complex. We demonstrate the pipeline's potential by exploring a recently-solved human bitter taste receptor, i.e. the TAS2R46, complexed with strychnine. We pinpointed 145 proteins sharing similar binding sites compared to the analysed bitter taste receptor and the enrichment analysis highlighted the related biological processes, molecular functions and cellular components. This work represents the foundation for future studies aimed at understanding the effective role of tastants outside the gustatory system: this could pave the way towards the rationalization of the diet as a supplement to standard pharmacological treatments and the design of novel tastants-inspired compounds to target other proteins involved in specific diseases or disorders. The proposed pipeline is publicly accessible, can be applied to any protein-ligand complex, and could be expanded to screen any database of protein structures.

PyPLIF HIPPOS-Assisted Prediction of Molecular Determinants of Ligand Binding to Receptors

Identification of molecular determinants of receptor-ligand binding could significantly increase the quality of structure-based virtual screening protocols. In turn, drug design process, especially the fragment-based approaches, could benefit from the knowledge. Retrospective virtual screening campaigns by employing AutoDock Vina followed by protein-ligand interaction fingerprinting (PLIF) identification by using recently published PyPLIF HIPPOS were the main techniques used here. The ligands and decoys datasets from the enhanced version of the database of useful decoys (DUDE) targeting human G protein-coupled receptors (GPCRs) were employed in this research since the mutation data are available and could be used to retrospectively verify the prediction. The results show that the method presented in this article could pinpoint some retrospectively verified molecular determinants. The method is therefore suggested to be employed as a routine in drug design and discovery.

FragRep: A Web Server for Structure-Based Drug Design by Fragment Replacement

The design of efficient computational tools for structure-guided ligand design is essential for the drug discovery process. We hereby present FragRep, a new web server for structure-based ligand design by fragment replacement. The input is a protein and a ligand structure, either from protein data bank or from molecular docking. Users can choose specific substructures they want to modify. The server tries to find suitable fragments that not only meet the geometric requirements of the remaining part of the ligand but also fit well with local protein environments. FragRep is a powerful computational tool for the rapid generation of ligand design ideas; either in scaffold hopping or bioisosteric replacing. The FragRep Server is freely available to researchers and can be accessed at http://xundrug.cn/fragrep.

Novel potent (dihydro)benzofuranyl piperazines as human histamine receptor ligands - Functional characterization and modeling studies on H3 and H4 receptors

Histamine acts through four different receptors (H₁R-H₄R), the H₃R and H₄R being the most explored in the last years as drug targets. The H₃R is a potential target to treat narcolepsy, Parkinson's disease, epilepsy, schizophrenia and several other CNS-related conditions, while H₄R blockade leads to anti-inflammatory and immunomodulatory effects. Our group has been exploring the dihydrobenzofuranyl-piperazines (LINS01 series) as human H₃R/H₄R ligands as potential drug candidates. In the present study, a set of 12 compounds were synthesized from adequate (dihydro)benzofuran synthons through simple reactions with corresponding piperazines, giving moderate to high yields. Four compounds (1b, 1f, 1g and 1h) showed high hH₃R affinity (pK_i > 7), compound 1h being the most potent (pK_i 8.4), and compound 1f showed the best efficiency (pKi 8.2, LE 0.53, LLE 5.85). BRET-based assays monitoring Gα_i activity indicated that the compounds are potent antagonists. Only one compound (2c, pK_i 7.1) presented high affinity for hH₄R. In contrast to what was observed for hH₃R, it showed partial agonist activity. Docking experiments indicated that bulky substituents occupy a hydrophobic pocket in hH₃R, while the N-allyl group forms favorable interactions with hydrophobic residues in the TM2, 3 and 7, increasing the selectivity towards hH₃R. Additionally, the importance of the indole NH in the interaction with Glu5.46 from hH₄R was confirmed by the modeling results, explaining the affinity and agonistic activity of compound 2c. The data reported in this work represent important findings for the rational design of future compounds for hH₃R and hH₄R.

PyPLIF HIPPOS: A Molecular Interaction Fingerprinting Tool for Docking Results of AutoDock Vina and PLANTS

We describe here our tool named PyPLIF HIPPOS, which was newly developed to analyze the docking results of AutoDock Vina and PLANTS. Its predecessor, PyPLIF (https://github.com/radifar/pyplif), is a molecular interaction fingerprinting tool for the docking results of PLANTS, exclusively. Unlike its predecessor, PyPLIF HIPPOS speeds up the computational times by separating the reference generation and docking analysis. PyPLIF HIPPOS also offers more options compared to PyPLIF. PyPLIF HIPPOS for Linux is stored as the Supporting Information in this application note and can be accessed in GitHub (https://github.com/radifar/PyPLIF-HIPPOS). Additionally, we present here the application of the tool in a retrospective structure-based virtual screening campaign targeting neuraminidase.

Performance of virtual screening against GPCR homology models: Impact of template selection and treatment of binding site plasticity

Rational drug design for G protein-coupled receptors (GPCRs) is limited by the small number of available atomic resolution structures. We assessed the use of homology modeling to predict the structures of two therapeutically relevant GPCRs and strategies to improve the performance of virtual screening against modeled binding sites. Homology models of the D2 dopamine (D2R) and serotonin 5-HT2A receptors (5-HT2AR) were generated based on crystal structures of 16 different GPCRs. Comparison of the homology models to D2R and 5-HT2AR crystal structures showed that accurate predictions could be obtained, but not necessarily using the most closely related template. Assessment of virtual screening performance was based on molecular docking of ligands and decoys. The results demonstrated that several templates and multiple models based on each of these must be evaluated to identify the optimal binding site structure. Models based on aminergic GPCRs showed substantial ligand enrichment and there was a trend toward improved virtual screening performance with increasing binding site accuracy. The best models even yielded ligand enrichment comparable to or better than that of the D2R and 5-HT2AR crystal structures. Methods to consider binding site plasticity were explored to further improve predictions. Molecular docking to ensembles of structures did not outperform the best individual binding site models, but could increase the diversity of hits from virtual screens and be advantageous for GPCR targets with few known ligands. Molecular dynamics refinement resulted in moderate improvements of structural accuracy and the virtual screening performance of snapshots was either comparable to or worse than that of the raw homology models. These results provide guidelines for successful application of structure-based ligand discovery using GPCR homology models.

4-(3-Aminoazetidin-1-yl)pyrimidin-2-amines as High-Affinity Non-imidazole Histamine H3 Receptor Agonists with in Vivo Central Nervous System Activity

Despite the high diversity of histamine H3 receptor (H3R) antagonist/inverse agonist structures, partial or full H3R agonists have typically been imidazole derivatives. An in-house screening campaign intriguingly afforded the non-imidazole 4-(3-azetidin-1-yl)pyrimidin-2-amine 11b as a partial H3R agonist. Here, the design, synthesis, and structure-activity relationships of 11b analogues are described. This series yields several non-imidazole full agonists with potencies varying with the alkyl substitution pattern on the basic amine following the in vitro evaluation of H3R agonism using a cyclic adenosine monophosphate response element-luciferase reporter gene assay. The key compound VUF16839 (14d) combines nanomolar on-target activity (pKi = 8.5, pEC50 = 9.5) with weak activity on cytochrome P450 enzymes and good metabolic stability. The proposed H3R binding mode of 14d indicates key interactions similar to those attained by histamine. In vivo evaluation of 14d in a social recognition test in mice revealed an amnesic effect at 5 mg/kg intraperitoneally. The excellent in vitro and in vivo pharmacological profiles and the non-imidazole structure of 14d make it a promising tool compound in H3R research.

Comparing structural and transcriptional drug networks reveals signatures of drug activity and toxicity in transcriptional responses

We performed an integrated analysis of drug chemical structures and drug-induced transcriptional responses. We demonstrated that a network representing three-dimensional structural similarities among 5452 compounds can be used to automatically group together drugs with similar scaffolds, physicochemical parameters and mode-of-action. We compared the structural network to a network representing transcriptional similarities among a subset of 1309 drugs for which transcriptional response were available in the Connectivity Map data set. Analysis of structurally similar, but transcriptionally different drugs sharing the same MOA enabled us to detect and remove weak and noisy transcriptional responses, greatly enhancing the reliability of transcription-based approaches to drug discovery and drug repositioning. Cardiac glycosides exhibited the strongest transcriptional responses with a significant induction of pathways related to epigenetic regulation, which suggests an epigenetic mechanism of action for these drugs. Drug classes with the weakest transcriptional responses tended to induce expression of cytochrome P450 enzymes, hinting at drug-induced drug resistance. Analysis of transcriptionally similar, but structurally different drugs with unrelated MOA, led us to the identification of a 'toxic' transcriptional signature indicative of lysosomal stress (lysosomotropism) and lipid accumulation (phospholipidosis) partially masking the target-specific transcriptional effects of these drugs. We found that this transcriptional signature is shared by 258 compounds and it is associated to the activation of the transcription factor TFEB, a master regulator of lysosomal biogenesis and autophagy. Finally, we built a predictive Random Forest model of these 258 compounds based on 128 physicochemical parameters, which should help in the early identification of potentially toxic drug candidates.

Transcriptional responses to drug treatment can reveal mechanism of action and off-target effects thus enabling drug repositioning, but only if measured in the appropriate cells at clinically relevant concentrations. A team led by Diego di Bernardo and Diego Medina generated a network representing structural similarities among compounds to automatically group together drugs with similar scaffolds and MOA. By comparing the structural drug network with a transcriptional drug network based on similarities in transcriptional response, the team observed broad differences between the two. This observation led to the identification of a transcriptional signature related lysosomal stress and phospholipidosis, and a physicochemical model to identify such compounds. These results provide general guidelines to prevent erroneous conclusion when using transcriptional responses of small molecules for drug discovery and drug repositioning

Identification of Histamine H3 Receptor Ligands Using a New Crystal Structure Fragment-based Method

Virtual screening offers an efficient alternative to high-throughput screening in the identification of pharmacological tools and lead compounds. Virtual screening is typically based on the matching of target structures or ligand pharmacophores to commercial or in-house compound catalogues. This study provides the first proof-of-concept for our recently reported method where pharmacophores are instead constructed based on the inference of residue-ligand fragments from crystal structures. We demonstrate its unique utility for G protein-coupled receptors, which represent the largest families of human membrane proteins and drug targets. We identified five neutral antagonists and one inverse agonist for the histamine H₃ receptor with potencies of 0.7-8.5 μM in a recombinant receptor cell-based inositol phosphate accumulation assay and validated their activity using a radioligand competition binding assay. H₃ receptor antagonism is of large therapeutic value and our ligands could serve as starting points for further lead optimisation. The six ligands exhibit four chemical scaffolds, whereof three have high novelty in comparison to the known H₃ receptor ligands in the ChEMBL database. The complete pharmacophore fragment library is freely available through the GPCR database, GPCRdb, allowing the successful application herein to be repeated for most of the 285 class A GPCR targets. The method could also easily be adapted to other protein families.

Do Fragments and Crystallization Additives Bind Similarly to Drug-like Ligands?

The success of fragment-based drug design (FBDD) hinges upon the optimization of low-molecular-weight compounds (MW < 300 Da) with weak binding affinities to lead compounds with high affinity and selectivity. Usually, structural information from fragment-protein complexes is used to develop ideas about the binding mode of similar but drug-like molecules. In this regard, crystallization additives such as cryoprotectants or buffer components, which are highly abundant in crystal structures, are frequently ignored. Thus, the aim of this study was to investigate the information present in protein complexes with fragments as well as those with additives and how they relate to the binding modes of their drug-like counterparts. We present a thorough analysis of the binding modes of crystallographic additives, fragments, and drug-like ligands bound to four diverse targets of wide interest in drug discovery and highly represented in the Protein Data Bank: cyclin-dependent kinase 2, β-secretase 1, carbonic anhydrase 2, and trypsin. We identified a total of 630 unique molecules bound to the catalytic binding sites, among them 31 additives, 222 fragments, and 377 drug-like ligands. In general, we observed that, independent of the target, protein-fragment interaction patterns are highly similar to those of drug-like ligands and mostly cover the residues crucial for binding. Crystallographic additives are also able to show conserved binding modes and recover the residues important for binding in some of the cases. Moreover, we show evidence that the information from fragments and drug-like ligands can be applied to rescore docking poses in order to improve the prediction of binding modes.

Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors

The growing number of studies on G protein-coupled receptors (GPCRs) family are a source of noticeable improvement in our understanding of the functioning of these proteins. GPCRs are responsible for a vast part of signaling in vertebrates and, as such, invariably remain in the spotlight of medicinal chemistry. A deeper insight into the underlying mechanisms of interesting phenomena observed in GPCRs, such as biased signaling or allosteric modulation, can be gained with experimental and computational studies. The latter play an important role in this process, since they allow for observations on scales inaccessible for most other methods. One of the key steps in such studies is proper computational reconstruction of actual ligand-receptor or protein-protein interactions, a process called molecular docking. A number of improvements and innovative applications of this method were documented recently. In this review, we focus particularly on innovations in docking to GPCRs.

Molecular Modelling Approaches for the Analysis of Histamine Receptors and Their Interaction with Ligands

Several experimental techniques to analyse histamine receptors are available, e.g. pharmacological characterisation of known or new compounds by different types of assays or mutagenesis studies. To obtain insights into the histamine receptors on a molecular and structural level, crystal structures have to be determined and molecular modelling studies have to be performed. It is widely accepted to generate homology models of the receptor of interest based on an appropriate crystal structure as a template and to refine the resulting models by molecular dynamic simulations. A lot of modelling techniques, e.g. docking, QSAR or interaction fingerprint methods, are used to predict binding modes of ligands and pharmacological data, e.g. affinity or even efficacy. However, within the last years, molecular dynamic simulations got more and more important: First of all, molecular dynamic simulations are very helpful to refine the binding mode of a ligand to a histamine receptor, obtained by docking studies. Furthermore, with increasing computational performance it got possible to simulate complete binding pathways of ions or ligands from the aqueous extracellular phase into the allosteric or orthosteric binding pocket of histamine receptors.

Improved Potency of Indole-Based NorA Efflux Pump Inhibitors: From Serendipity toward Rational Design and Development

The NorA efflux pump is a potential drug target for reversal of resistance to selected antibacterial agents, and recently we described indole-based inhibitor candidates. Herein we report a second class of inhibitors derived from them but with significant differences in shape and size. In particular, compounds 13 and 14 are very potent inhibitors in that they demonstrated the lowest IC₅₀ values (2 μM) ever observed among all indole-based compounds we have evaluated.

Molecular interaction fingerprint approaches for GPCR drug discovery

Protein-ligand interaction fingerprints (IFPs) are binary 1D representations of the 3D structure of protein-ligand complexes encoding the presence or absence of specific interactions between the binding pocket amino acids and the ligand. Various implementations of IFPs have been developed and successfully applied for post-processing molecular docking results for G Protein-Coupled Receptor (GPCR) ligand binding mode prediction and virtual ligand screening. Novel interaction fingerprint methods enable structural chemogenomics and polypharmacology predictions by complementing the increasing amount of GPCR structural data. Machine learning methods are increasingly used to derive relationships between bioactivity data and fingerprint descriptors of chemical and structural information of binding sites, ligands, and protein-ligand interactions. Factors that influence the application of IFPs include structure preparation, binding site definition, fingerprint similarity assessment, and data processing and these factors pose challenges as well possibilities to optimize interaction fingerprint methods for GPCR drug discovery.

Drug Design for CNS Diseases: Polypharmacological Profiling of Compounds Using Cheminformatic, 3D-QSAR and Virtual Screening Methodologies

HIGHLIGHTS

Many CNS targets are being explored for multi-target drug design

New databases and cheminformatic methods enable prediction of primary pharmaceutical target and off-targets of compounds

QSAR, virtual screening and docking methods increase the potential of rational drug design

The diverse cerebral mechanisms implicated in Central Nervous System (CNS) diseases together with the heterogeneous and overlapping nature of phenotypes indicated that multitarget strategies may be appropriate for the improved treatment of complex brain diseases. Understanding how the neurotransmitter systems interact is also important in optimizing therapeutic strategies. Pharmacological intervention on one target will often influence another one, such as the well-established serotonin-dopamine interaction or the dopamine-glutamate interaction. It is now accepted that drug action can involve plural targets and that polypharmacological interaction with multiple targets, to address disease in more subtle and effective ways, is a key concept for development of novel drug candidates against complex CNS diseases. A multi-target therapeutic strategy for Alzheimer‘s disease resulted in the development of very effective Multi-Target Designed Ligands (MTDL) that act on both the cholinergic and monoaminergic systems, and also retard the progression of neurodegeneration by inhibiting amyloid aggregation. Many compounds already in databases have been investigated as ligands for multiple targets in drug-discovery programs. A probabilistic method, the Parzen-Rosenblatt Window approach, was used to build a “predictor” model using data collected from the ChEMBL database. The model can be used to predict both the primary pharmaceutical target and off-targets of a compound based on its structure. Several multi-target ligands were selected for further study, as compounds with possible additional beneficial pharmacological activities. Based on all these findings, it is concluded that multipotent ligands targeting AChE/MAO-A/MAO-B and also D₁-R/D₂-R/5-HT_2A-R/H₃-R are promising novel drug candidates with improved efficacy and beneficial neuroleptic and procognitive activities in treatment of Alzheimer's and related neurodegenerative diseases. Structural information for drug targets permits docking and virtual screening and exploration of the molecular determinants of binding, hence facilitating the design of multi-targeted drugs. The crystal structures and models of enzymes of the monoaminergic and cholinergic systems have been used to investigate the structural origins of target selectivity and to identify molecular determinants, in order to design MTDLs.

A combined ligand- and structure-based approach for the identification of rilmenidine-derived compounds which synergize the antitumor effects of doxorubicin

The clonidine-like central antihypertensive agent rilmenidine, which has high affinity for I1-type imidazoline receptors (I1-IR) was recently found to have cytotoxic effects on cultured cancer cell lines. However, due to its pharmacological effects resulting also from α2-adrenoceptor activation, rilmenidine cannot be considered a suitable anticancer drug candidate. Here, we report the identification of novel rilmenidine-derived compounds with anticancer potential and devoid of α2-adrenoceptor effects by means of ligand- and structure-based drug design approaches. Starting from a large virtual library, eleven compounds were selected, synthesized and submitted to biological evaluation. The most active compound 5 exhibited a cytotoxic profile similar to that of rilmenidine, but without appreciable affinity to α2-adrenoceptors. In addition, compound 5 significantly enhanced the apoptotic response to doxorubicin, and may thus represent an important tool for the development of better adjuvant chemotherapeutic strategies for doxorubicin-insensitive cancers.

Combinatorial Consensus Scoring for Ligand-Based Virtual Fragment Screening: A Comparative Case Study for Serotonin 5-HT(3)A, Histamine H(1), and Histamine H(4) Receptors

In the current study we have evaluated the applicability of ligand-based virtual screening (LBVS) methods for the identification of small fragment-like biologically active molecules using different similarity descriptors and different consensus scoring approaches. For this purpose, we have evaluated the performance of 14 chemical similarity descriptors in retrospective virtual screening studies to discriminate fragment-like ligands of three membrane-bound receptors from fragments that are experimentally determined to have no affinity for these proteins (true inactives). We used a complete fragment affinity data set of experimentally determined ligands and inactives for two G protein-coupled receptors (GPCRs), the histamine H1 receptor (H1R) and the histamine H4 receptor (H4R), and one ligand-gated ion channel (LGIC), the serotonin receptor (5-HT3AR), to validate our retrospective virtual screening studies. We have exhaustively tested consensus scoring strategies that combine the results of multiple actives (group fusion) or combine different similarity descriptors (similarity fusion), and for the first time systematically evaluated different combinations of group fusion and similarity fusion approaches. Our studies show that for these three case study protein targets both consensus scoring approaches can increase virtual screening enrichments compared to single chemical similarity search methods. Our cheminformatics analyses recommend to use a combination of both group fusion and similarity fusion for prospective ligand-based virtual fragment screening.

Improvements, trends, and new ideas in molecular docking: 2012-2013 in review

Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher-level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods.

Identifying modulators of CXC receptors 3 and 4 with tailored selectivity using multi-target docking

The G protein-coupled receptors of the C-X-C subfamily form a group among the chemokine receptors whose endogenous ligands are peptides with a common Cys-X-Cys motif. The CXC chemokine receptors 3 and 4 (CXCR3, CXCR4), which are investigated in this study, are linked to severe diseases such as cancer, multiple sclerosis, and HIV infections. Of particular interest, this receptor pair potentially forms a target for a polypharmacological drug treatment. Considering known ligands from public databases, such dual binders have not been identified yet. We therefore applied large-scale docking to the structure of CXCR4 and a homology model of CXCR3 with the goal to predict such dual binders, as well as compounds selective for either one of the receptors. Using signaling and biochemical assays, we showed that more than 50% of these predictions were correct in each category, yielding ligands with excellent binding efficiencies. These results highlight that docking is a suitable tool for the identification of ligands with tailored binding profiles to GPCRs, even when using homology models. More importantly, we present novel CXCR3-CXCR4 dual modulators that might pave the road to understanding the mechanisms of polypharmacological inhibition of these receptors.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

Computational methods for fragment-based ligand design: growing and linking

Fragment-based drug design has proved itself as a powerful technique for increasing the sampling and diversity of chemical space and enabling the design of novel leads and compounds. Computational techniques for identifying fragments, binding sites and particularly for linking, growing, and evolving fragments play a significant role in the process. Information from ADME studies and clustering property information in the form of toxicophores and chemotypes can play a significant role in aiding the design of novel, selective fragments with good activity profiles.

Structure-based virtual screening for fragment-like ligands of the G protein-coupled histamine H4 receptor

Structure-based virtual screening using H₁R- and β₂R-based histamine H₄R homology models identified 9 fragments with an affinity ranging from 0.14 to 6.3 μm for H₄R.

Sequential application of ligand and structure based modeling approaches to index chemicals for their hH4R antagonism

The human histamine H4 receptor (hH4R), a member of the G-protein coupled receptors (GPCR) family, is an increasingly attractive drug target. It plays a key role in many cell pathways and many hH4R ligands are studied for the treatment of several inflammatory, allergic and autoimmune disorders, as well as for analgesic activity. Due to the challenging difficulties in the experimental elucidation of hH4R structure, virtual screening campaigns are normally run on homology based models. However, a wealth of information about the chemical properties of GPCR ligands has also accumulated over the last few years and an appropriate combination of these ligand-based knowledge with structure-based molecular modeling studies emerges as a promising strategy for computer-assisted drug design. Here, two chemoinformatics techniques, the Intelligent Learning Engine (ILE) and Iterative Stochastic Elimination (ISE) approach, were used to index chemicals for their hH4R bioactivity. An application of the prediction model on external test set composed of more than 160 hH4R antagonists picked from the chEMBL database gave enrichment factor of 16.4. A virtual high throughput screening on ZINC database was carried out, picking ∼ 4000 chemicals highly indexed as H4R antagonists' candidates. Next, a series of 3D models of hH4R were generated by molecular modeling and molecular dynamics simulations performed in fully atomistic lipid membranes. The efficacy of the hH4R 3D models in discrimination between actives and non-actives were checked and the 3D model with the best performance was chosen for further docking studies performed on the focused library. The output of these docking studies was a consensus library of 11 highly active scored drug candidates. Our findings suggest that a sequential combination of ligand-based chemoinformatics approaches with structure-based ones has the potential to improve the success rate in discovering new biologically active GPCR drugs and increase the enrichment factors in a synergistic manner.

Multiple fragment docking and linking in primary and secondary pockets of dopamine receptors

A sequential docking methodology was applied to computationally predict starting points for fragment linking using the human dopamine D3 receptor crystal structure and a human dopamine D2 receptor homology model. Two focused fragment libraries were docked in the primary and secondary binding sites, and best fragment combinations were enumerated. Similar top scoring fragments were found for the primary site, while secondary site fragments were predicted to convey selectivity. Three linked compounds were synthesized that had 9-, 39-, and 55-fold selectivity in favor of D3 and the subtype selectivity of the compounds was assessed on a structural basis.

A structural chemogenomics analysis of aminergic GPCRs: lessons for histamine receptor ligand design

BACKGROUND AND PURPOSE

Chemogenomics focuses on the discovery of new connections between chemical and biological space leading to the discovery of new protein targets and biologically active molecules. G-protein coupled receptors (GPCRs) are a particularly interesting protein family for chemogenomics studies because there is an overwhelming amount of ligand binding affinity data available. The increasing number of aminergic GPCR crystal structures now for the first time allows the integration of chemogenomics studies with high-resolution structural analyses of GPCR-ligand complexes.

EXPERIMENTAL APPROACH

In this study, we have combined ligand affinity data, receptor mutagenesis studies, and amino acid sequence analyses to high-resolution structural analyses of (hist)aminergic GPCR-ligand interactions. This integrated structural chemogenomics analysis is used to more accurately describe the molecular and structural determinants of ligand affinity and selectivity in different key binding regions of the crystallized aminergic GPCRs, and histamine receptors in particular.

KEY RESULTS

Our investigations highlight interesting correlations and differences between ligand similarity and ligand binding site similarity of different aminergic receptors. Apparent discrepancies can be explained by combining detailed analysis of crystallized or predicted protein-ligand binding modes, receptor mutation studies, and ligand structure-selectivity relationships that identify local differences in essential pharmacophore features in the ligand binding sites of different receptors.

CONCLUSIONS AND IMPLICATIONS

We have performed structural chemogenomics studies that identify links between (hist)aminergic receptor ligands and their binding sites and binding modes. This knowledge can be used to identify structure-selectivity relationships that increase our understanding of ligand binding to (hist)aminergic receptors and hence can be used in future GPCR ligand discovery and design.

Detailed analysis of biased histamine H₄ receptor signalling by JNJ 7777120 analogues

BACKGROUND AND PURPOSE

The histamine H₄ receptor, originally thought to signal merely through Gαi proteins, has recently been shown to also recruit and signal via β-arrestin2. Following the discovery that the reference antagonist indolecarboxamide JNJ 7777120 appears to be a partial agonist in β-arrestin2 recruitment, we have identified additional biased hH₄R ligands that preferentially couple to Gαi or β-arrestin2 proteins. In this study, we explored ligand and receptor regions that are important for biased hH₄R signalling.

EXPERIMENTAL APPROACH

We evaluated a series of 48 indolecarboxamides with subtle structural differences for their ability to induce hH₄R-mediated Gαi protein signalling or β-arrestin2 recruitment. Subsequently, a Fingerprints for Ligands and Proteins three-dimensional quantitative structure-activity relationship analysis correlated intrinsic activity values with structural ligand requirements. Moreover, a hH₄R homology model was used to identify receptor regions important for biased hH₄R signalling.

KEY RESULTS

One indolecarboxamide (75) with a nitro substituent on position R7 of the aromatic ring displayed an equal preference for the Gαi and β-arrestin2 pathway and was classified as unbiased hH₄R ligand. The other 47 indolecarboxamides were β-arrestin2-biased agonists. Intrinsic activities of the unbiased as well as β-arrestin2-biased indolecarboxamides to induce β-arrestin2 recruitment could be correlated with different ligand features and hH₄R regions.

CONCLUSION AND IMPLICATIONS

Small structural modifications resulted in diverse intrinsic activities for unbiased (75) and β-arrestin2-biased indolecarboxamides. Analysis of ligand and receptor features revealed efficacy hotspots responsible for biased-β-arrestin2 recruitment. This knowledge is useful for the design of hH₄R ligands with biased intrinsic activities and aids our understanding of the mechanism of H₄R activation.

The receptor concept in 3D: from hypothesis and metaphor to GPCR-ligand structures

The first mentioning of the word "receptor" for the structure with which a bioactive compound should react for obtaining its specific influence on a physiological system goes back to the years around 1900. The receptor concept was adapted from the lock and key theory for the enzyme substrate and blockers interactions. Through the years the concept, in the beginning rather being a metaphor, not a model, was refined and became reality in recent years. Not only the structures of receptors were elucidated, also the receptor machineries were unraveled. Following a brief historical review we will describe how the recent breakthroughs in the experimental determination of G protein-coupled receptor (GPCR) crystal structures can be complemented by computational modeling, medicinal chemistry, biochemical, and molecular pharmacological studies to obtain new insights into the molecular determinants of GPCR-ligand binding and activation. We will furthermore discuss how this information can be used for structure-based discovery of novel GPCR ligands that bind specific (allosteric) binding sites with desired effects on GPCR functional activity.

Structure-based discovery of selective serotonin 5-HT(1B) receptor ligands

The development of safe and effective drugs relies on the discovery of selective ligands. Serotonin (5-hydroxytryptamine [5-HT]) G protein-coupled receptors are therapeutic targets for CNS disorders but are also associated with adverse drug effects. The determination of crystal structures for the 5-HT1B and 5-HT2B receptors provided an opportunity to identify subtype selective ligands using structure-based methods. From docking screens of 1.3 million compounds, 22 molecules were predicted to be selective for the 5-HT1B receptor over the 5-HT2B subtype, a requirement for safe serotonergic drugs. Nine compounds were experimentally verified as 5-HT1B-selective ligands, with up to 300-fold higher affinities for this subtype. Three of the ligands were agonists of the G protein pathway. Analysis of state-of-the-art homology models of the two 5-HT receptors revealed that the crystal structures were critical for predicting selective ligands. Our results demonstrate that structure-based screening can guide the discovery of ligands with specific selectivity profiles.

Virtual fragment screening on GPCRs: a case study on dopamine D3 and histamine H4 receptors

Prospective structure based virtual fragment screening methodologies on two GPCR targets namely the dopamine D3 and the histamine H4 receptors with a library of 12,905 fragments were evaluated. Fragments were docked to the X-ray structure and the homology model of the D3 and H4 receptors, respectively. Representative receptor conformations for ensemble docking were obtained from molecular dynamics trajectories. In vitro confirmed hit rates ranged from 16% to 32%. Hits had high ligand efficiency (LE) values in the range of 0.31-0.74 and also acceptable lipophilic efficiency. The X-ray structure, the homology model and structural ensembles were all found suitable for docking based virtual screening of fragments against these GPCRs. However, there was little overlap among different hit sets and methodologies were thus complementary to each other.

Isozyme-specific ligands for O-acetylserine sulfhydrylase, a novel antibiotic target

The last step of cysteine biosynthesis in bacteria and plants is catalyzed by O-acetylserine sulfhydrylase. In bacteria, two isozymes, O-acetylserine sulfhydrylase-A and O-acetylserine sulfhydrylase-B, have been identified that share similar binding sites, although the respective specific functions are still debated. O-acetylserine sulfhydrylase plays a key role in the adaptation of bacteria to the host environment, in the defense mechanisms to oxidative stress and in antibiotic resistance. Because mammals synthesize cysteine from methionine and lack O-acetylserine sulfhydrylase, the enzyme is a potential target for antimicrobials. With this aim, we first identified potential inhibitors of the two isozymes via a ligand- and structure-based in silico screening of a subset of the ZINC library using FLAP. The binding affinities of the most promising candidates were measured in vitro on purified O-acetylserine sulfhydrylase-A and O-acetylserine sulfhydrylase-B from Salmonella typhimurium by a direct method that exploits the change in the cofactor fluorescence. Two molecules were identified with dissociation constants of 3.7 and 33 µM for O-acetylserine sulfhydrylase-A and O-acetylserine sulfhydrylase-B, respectively. Because GRID analysis of the two isoenzymes indicates the presence of a few common pharmacophoric features, cross binding titrations were carried out. It was found that the best binder for O-acetylserine sulfhydrylase-B exhibits a dissociation constant of 29 µM for O-acetylserine sulfhydrylase-A, thus displaying a limited selectivity, whereas the best binder for O-acetylserine sulfhydrylase-A exhibits a dissociation constant of 50 µM for O-acetylserine sulfhydrylase-B and is thus 8-fold selective towards the former isozyme. Therefore, isoform-specific and isoform-independent ligands allow to either selectively target the isozyme that predominantly supports bacteria during infection and long-term survival or to completely block bacterial cysteine biosynthesis.

Complementarity between in silico and biophysical screening approaches in fragment-based lead discovery against the A(2A) adenosine receptor

Fragment-based lead discovery (FBLD) is becoming an increasingly important method in drug development. We have explored the potential to complement NMR-based biophysical screening of chemical libraries with molecular docking in FBLD against the A(2A) adenosine receptor (A(2A)AR), a drug target for inflammation and Parkinson's disease. Prior to an NMR-based screen of a fragment library against the A(2A)AR, molecular docking against a crystal structure was used to rank the same set of molecules by their predicted affinities. Molecular docking was able to predict four out of the five orthosteric ligands discovered by NMR among the top 5% of the ranked library, suggesting that structure-based methods could be used to prioritize among primary hits from biophysical screens. In addition, three fragments that were top-ranked by molecular docking, but had not been picked up by the NMR-based method, were demonstrated to be A(2A)AR ligands. While biophysical approaches for fragment screening are typically limited to a few thousand compounds, the docking screen was extended to include 328,000 commercially available fragments. Twenty-two top-ranked compounds were tested in radioligand binding assays, and 14 of these were A(2A)AR ligands with K(i) values ranging from 2 to 240 μM. Optimization of fragments was guided by molecular dynamics simulations and free energy calculations. The results illuminate strengths and weaknesses of molecular docking and demonstrate that this method can serve as a valuable complementary tool to biophysical screening in FBLD.

Fragment-based lead discovery on G-protein-coupled receptors

INTRODUCTION

G-protein-coupled receptors (GPCRs) form one of the largest groups of potential targets for novel medications. Low druggability of many GPCR targets and inefficient sampling of chemical space in high-throughput screening expertise however often hinder discovery of drug discovery leads for GPCRs. Fragment-based drug discovery is an alternative approach to the conventional strategy and has proven its efficiency on several enzyme targets. Based on developments in biophysical screening techniques, receptor stabilization and in vitro assays, virtual and experimental fragment screening and fragment-based lead discovery recently became applicable for GPCR targets.

AREAS COVERED

This article provides a review of the biophysical as well as biological detection techniques suitable to study GPCRs together with their applications to screen fragment libraries and identify fragment-size ligands of cell surface receptors. The article presents several recent examples including both virtual and experimental protocols for fragment hit discovery and early hit to lead progress.

EXPERT OPINION

With the recent progress in biophysical detection techniques, the advantages of fragment-based drug discovery could be exploited for GPCR targets. Structural information on GPCRs will be more abundantly available for early stages of drug discovery projects, providing information on the binding process and efficiently supporting the progression of fragment hit to lead. In silico approaches in combination with biological assays can be used to address structurally challenging GPCRs and confirm biological relevance of interaction early in the drug discovery project.

Small and colorful stones make beautiful mosaics: fragment-based chemogenomics

Smaller stones with a wide variety of colors make a higher resolution mosaic. In much the same way, smaller chemical entities that are structurally diverse are better able to interrogate protein binding sites. This feature article describes the construction of a diverse fragment library and an analysis of the screening of six representative protein targets belonging to three diverse target classes (G protein-coupled receptors ADRB2, H1R, H3R, and H4R, the ligand-gated ion channel 5-HT3R, and the kinase PKA) using chemogenomics approaches. The integration of experimentally determined bioaffinity profiles across related and unrelated protein targets and chemogenomics analysis of fragment binding and protein structure allow the identification of: (i) unexpected similarities and differences in ligand binding properties, and (ii) subtle ligand affinity and selectivity cliffs. With a wealth of fragment screening data being generated in industry and academia, such approaches will contribute to a more detailed structural understanding of ligand-protein interactions.