Discovery of novel human histamine H4 receptor ligands by large-scale structure-based virtual screening

Docking of T6361 Analogues as Potential Inhibitors of E.coli MurA Followed by ADME-Toxicity Study

BACKGROUND

Developing more potent antibacterial agents is one of the most important tasks of scientists in the health field due to the problem of antibiotic resistance. Among the antibiotic targets, we mention MurA (UDP-N-Acetylglucosamine Enolpyruvyl Transferase), which is a key enzyme of peptidoglycan biosynthesis of the bacterial cell wall.

OBJECTIVE

Our objective was to search for new inhibitors of the bacterial enzyme MurA by docking the analogues of its inhibitor T6361, a derivative of 5-sulfonoxy-anthranilic acid.

METHODS

990 analogues of T6361 were docked in the first binding site of E.coli MurA (open form) using the FlexX program, and the ADME-Toxicity profile of the best ones was evaluated by SwissADME and PreADMET web servers. .

RESULTS

Docking results revealed two T6361 analogues to provide better energy scores than T6361, and have similar interactions with the binding site of E.coli MurA namely,3-{[2-(piperidine-1-carbonyl) phenyl]sulfamoyl}benzoic acid and 3-{[2-(pyrrolidine-1 carbonyl)phenyl]sulfamoyl}benzoic acid. Moreover, the two molecules were found to possess good pharmacokinetics and low toxicity.

CONCLUSION

We propose two analogues of T6361 as new potential inhibitors of MurA enzyme. Their good ADME-Toxicity profile qualifies them to reach in vitro and in vivo assays as future lead molecules.

Computational Chemistry for the Identification of Lead Compounds for Radiotracer Development

The use of computer-aided drug design (CADD) for the identification of lead compounds in radiotracer development is steadily increasing. Traditional CADD methods, such as structure-based and ligand-based virtual screening and optimization, have been successfully utilized in many drug discovery programs and are highlighted throughout this review. First, we discuss the use of virtual screening for hit identification at the beginning of drug discovery programs. This is followed by an analysis of how the hits derived from virtual screening can be filtered and culled to highly probable candidates to test in in vitro assays. We then illustrate how CADD can be used to optimize the potency of experimentally validated hit compounds from virtual screening for use in positron emission tomography (PET). Finally, we conclude with a survey of the newest techniques in CADD employing machine learning (ML).

FSDscore: An Effective Target-focused Scoring Criterion for Virtual Screening

Improving screening efficiency is one of the most challenging tasks of virtual screening (VS). In this work, we propose an effective target-focused scoring criterion for VS and apply it to the screening of a specific target scaffold replacement library constructed by enumeration of suitable substitution fragments and R-groups of known ligands. This criterion is based on both ligand- and structure-based scoring methods, which includes feature maps, 3D shape similarity, and the pairwise distance information between proteins and ligands (FSDscore). It is precisely due to the hybrid advantages of ligand- and structure-based approaches that FSDscore performs far better on the validation dataset than other scoring methods. We apply FSDscore to the VS of different kinase targets, MERTK (Mer tyrosine kinase) and ABL1 (tyrosine-protein kinase ABL1) in order to avoid occasionality. Finally, a VS case study shows the potential and effectiveness of our scoring criterion in drug discovery and molecular dynamics simulation further verifies its powerful ability.

Structure-based discovery of nonopioid analgesics acting through the α2A-adrenergic receptor

Because nonopioid analgesics are much sought after, we computationally docked more than 301 million virtual molecules against a validated pain target, the α2A-adrenergic receptor (α2AAR), seeking new α2AAR agonists chemotypes that lack the sedation conferred by known α2AAR drugs, such as dexmedetomidine. We identified 17 ligands with potencies as low as 12 nanomolar, many with partial agonism and preferential Gi and Go signaling. Experimental structures of α2AAR complexed with two of these agonists confirmed the docking predictions and templated further optimization. Several compounds, including the initial docking hit '9087 [mean effective concentration (EC50) of 52 nanomolar] and two analogs, '7075 and PS75 (EC50 4.1 and 4.8 nanomolar), exerted on-target analgesic activity in multiple in vivo pain models without sedation. These newly discovered agonists are interesting as therapeutic leads that lack the liabilities of opioids and the sedation of dexmedetomidine.

Property-Unmatched Decoys in Docking Benchmarks

Enrichment of ligands versus property-matched decoys is widely used to test and optimize docking library screens. However, the unconstrained optimization of enrichment alone can mislead, leading to false confidence in prospective performance. This can arise by over-optimizing for enrichment against property-matched decoys, without considering the full spectrum of molecules to be found in a true large library screen. Adding decoys representing charge extrema helps mitigate over-optimizing for electrostatic interactions. Adding decoys that represent the overall characteristics of the library to be docked allows one to sample molecules not represented by ligands and property-matched decoys but that one will encounter in a prospective screen. An optimized version of the DUD-E set (DUDE-Z), as well as Extrema and sets representing broad features of the library (Goldilocks), is developed here. We also explore the variability that one can encounter in enrichment calculations and how that can temper one's confidence in small enrichment differences. The new tools and new decoy sets are freely available at http://tldr.docking.org and http://dudez.docking.org.

Improving Protein-Ligand Docking Results with High-Throughput Molecular Dynamics Simulations

Structure-based virtual screening relies on classical scoring functions that often fail to reliably discriminate binders from nonbinders. In this work, we present a high-throughput protein-ligand complex molecular dynamics (MD) simulation that uses the output from AutoDock Vina to improve docking results in distinguishing active from decoy ligands in a directory of useful decoy-enhanced (DUD-E) dataset. MD trajectories are processed by evaluating ligand-binding stability using root-mean-square deviations. We select 56 protein targets (of 7 different protein classes) and 560 ligands (280 actives, 280 decoys) and show 22% improvement in ROC AUC (area under the curve, receiver operating characteristics curve), from an initial value of 0.68 (AutoDock Vina) to a final value of 0.83. The MD simulation demonstrates a robust performance across all seven different protein classes. In addition, some predicted ligand-binding modes are moderately refined during MD simulations. These results systematically validate the reliability of a physics-based approach to evaluate protein-ligand binding interactions.

Accurate Hit Estimation for Iterative Screening Using Venn-ABERS Predictors

Iterative screening has emerged as a promising approach to increase the efficiency of high-throughput screening (HTS) campaigns in drug discovery. By learning from a subset of the compound library, inferences on what compounds to screen next can be made by predictive models. One of the challenges of iterative screening is to decide how many iterations to perform. This is mainly related to difficulties in estimating the prospective hit rate in any given iteration. In this article, a novel method based on Venn-ABERS predictors is proposed. The method provides accurate estimates of the number of hits retrieved in any given iteration during an HTS campaign. The estimates provide the necessary information to support the decision on the number of iterations needed to maximize the screening outcome. Thus, this method offers a prospective screening strategy for early-stage drug discovery.

Ultra-large library docking for discovering new chemotypes

Despite intense interest in expanding chemical space, libraries containing hundreds-of-millions to billions of diverse molecules have remained inaccessible. Here we investigate structure-based docking of 170 million make-on-demand compounds from 130 well-characterized reactions. The resulting library is diverse, representing over 10.7 million scaffolds that are otherwise unavailable. For each compound in the library, docking against AmpC β-lactamase (AmpC) and the D4 dopamine receptor were simulated. From the top-ranking molecules, 44 and 549 compounds were synthesized and tested for interactions with AmpC and the D4 dopamine receptor, respectively. We found a phenolate inhibitor of AmpC, which revealed a group of inhibitors without known precedent. This molecule was optimized to 77 nM, which places it among the most potent non-covalent AmpC inhibitors known. Crystal structures of this and other AmpC inhibitors confirmed the docking predictions. Against the D4 dopamine receptor, hit rates fell almost monotonically with docking score, and a hit-rate versus score curve predicted that the library contained 453,000 ligands for the D4 dopamine receptor. Of 81 new chemotypes discovered, 30 showed submicromolar activity, including a 180-pM subtype-selective agonist of the D4 dopamine receptor.

Strategies of Virtual Screening in Medicinal Chemistry

Virtual screening represents an effective computational strategy to rise-up the chances of finding new bioactive compounds by accelerating the time needed to move from an initial intuition to market. Classically, the most pursued approaches rely on ligand- and structure-based studies, the former employed when structural data information about the target is missing while the latter employed when X-ray/NMR solved or homology models are instead available for the target. The authors will focus on the most advanced techniques applied in this area. In particular, they will survey the key concepts of virtual screening by discussing how to properly select chemical libraries, how to make database curation, how to applying and- and structure-based techniques, how to wisely use post-processing methods. Emphasis will be also given to the most meaningful databases used in VS protocols. For the ease of discussion several examples will be presented.

Anti-inflammatory effect of a new piperazine derivative: (4-methylpiperazin-1-yl)(1-phenyl-1H-pyrazol-4-yl)methanone

AIMS

This study investigates the anti-nociceptive and anti-inflammatory effects of new piperazine compound (LQFM182) as well as the toxicity acute in vitro.

MAIN METHODS

To evaluate the anti-nociceptive activity, the acetic acid-induced abdominal writhing test, tail flick test and formalin-induced pain test were used. The anti-inflammatory activity was evaluated using the models of paw oedema and pleurisy induced by carrageenan and some inflammatory parameters were evaluated, including cell migration, myeloperoxidase enzyme activity and the levels of TNF-α and IL-1β cytokines in pleural exudate. The acute oral systemic toxicity of LQFM182 in mice was evaluated through the neutral red uptake (nru) assay.

KEY FINDINGS

LQFM182 (50, 100 or 200 mg/kg, p.o.) decreased the number of writhings induced by acetic acid in a dose-dependent manner, and an intermediate dose (100 mg/kg, p.o.) reduced the paw licking time of animals in the second phase of the formalin test. Furthermore, LQFM182 (100 mg/kg, p.o.) reduced oedema formation at all hours of the paw oedema induced by carrageenan test and in pleurisy test reduced cell migration from the reduction of polymorphonuclear cells, myeloperoxidase enzyme activity and the levels of pro-inflammatory cytokines IL-1β and TNF-α. Therefore, it was classified in GHS category 300 < LD₅₀ < 2000 mg/kg.

SIGNIFICANCE

Reduction of the TNF-α and IL-1β levels.

Performance of machine-learning scoring functions in structure-based virtual screening

Classical scoring functions have reached a plateau in their performance in virtual screening and binding affinity prediction. Recently, machine-learning scoring functions trained on protein-ligand complexes have shown great promise in small tailored studies. They have also raised controversy, specifically concerning model overfitting and applicability to novel targets. Here we provide a new ready-to-use scoring function (RF-Score-VS) trained on 15 426 active and 893 897 inactive molecules docked to a set of 102 targets. We use the full DUD-E data sets along with three docking tools, five classical and three machine-learning scoring functions for model building and performance assessment. Our results show RF-Score-VS can substantially improve virtual screening performance: RF-Score-VS top 1% provides 55.6% hit rate, whereas that of Vina only 16.2% (for smaller percent the difference is even more encouraging: RF-Score-VS top 0.1% achieves 88.6% hit rate for 27.5% using Vina). In addition, RF-Score-VS provides much better prediction of measured binding affinity than Vina (Pearson correlation of 0.56 and −0.18, respectively). Lastly, we test RF-Score-VS on an independent test set from the DEKOIS benchmark and observed comparable results. We provide full data sets to facilitate further research in this area (http://github.com/oddt/rfscorevs) as well as ready-to-use RF-Score-VS (http://github.com/oddt/rfscorevs_binary).

Identifying ligands at orphan GPCRs: current status using structure-based approaches

GPCRs are the most successful pharmaceutical targets in history. Nevertheless, the pharmacology of many GPCRs remains inaccessible as their endogenous or exogenous modulators have not been discovered. Tools that explore the physiological functions and pharmacological potential of these 'orphan' GPCRs, whether they are endogenous and/or surrogate ligands, are therefore of paramount importance. Rates of receptor deorphanization determined by traditional reverse pharmacology methods have slowed, indicating a need for the development of more sophisticated and efficient ligand screening approaches. Here, we discuss the use of structure-based ligand discovery approaches to identify small molecule modulators for exploring the function of orphan GPCRs. These studies have been buoyed by the growing number of GPCR crystal structures solved in the past decade, providing a broad range of template structures for homology modelling of orphans. This review discusses the methods used to establish the appropriate signalling assays to test orphan receptor activity and provides current examples of structure-based methods used to identify ligands of orphan GPCRs. Linked Articles This article is part of a themed section on Molecular Pharmacology of G Protein-Coupled Receptors. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v173.20/issuetoc.

Discovery of nanomolar ligands with novel scaffolds for the histamine H4 receptor by virtual screening

The involvement of histamine H4 receptor (H4R) in immune cells chemotaxis and mediator release makes it an attractive target for the treatment of inflammation disorders. A decade of medicinal chemistry efforts has led to several promising ligands, although the chemical structures described so far possesses a singular limited diversity. We report here the discovery of novel structures, belonging to completely different scaffolds. The virtual screening was planed as a two-steps process. First, using a "scout screening" methodology, we have experimentally probed the H4R ligand binding site using a small size chemical library with very diverse structures, and identified a hit that further assist us in refining a raw 3D homology model. Second, the refined 3D model was used to conduct a widened virtual screening. This two-steps strategy proved to be very successful, both in terms of structural diversity and hit rate (23%). Moreover, the hits have high affinity for the H4R, with most potent ligands in the nanomolar range.

Role of computer-aided drug design in modern drug discovery

Drug discovery utilizes chemical biology and computational drug design approaches for the efficient identification and optimization of lead compounds. Chemical biology is mostly involved in the elucidation of the biological function of a target and the mechanism of action of a chemical modulator. On the other hand, computer-aided drug design makes use of the structural knowledge of either the target (structure-based) or known ligands with bioactivity (ligand-based) to facilitate the determination of promising candidate drugs. Various virtual screening techniques are now being used by both pharmaceutical companies and academic research groups to reduce the cost and time required for the discovery of a potent drug. Despite the rapid advances in these methods, continuous improvements are critical for future drug discovery tools. Advantages presented by structure-based and ligand-based drug design suggest that their complementary use, as well as their integration with experimental routines, has a powerful impact on rational drug design. In this article, we give an overview of the current computational drug design and their application in integrated rational drug development to aid in the progress of drug discovery research.

What can crystal structures of aminergic receptors tell us about designing subtype-selective ligands?

G protein-coupled receptors (GPCRs) are integral membrane proteins that represent an important class of drug targets. In particular, aminergic GPCRs interact with a significant portion of drugs currently on the market. However, most drugs that target these receptors are associated with undesirable side effects, which are due in part to promiscuous interactions with close homologs of the intended target receptors. Here, based on a systematic analysis of all 37 of the currently available high-resolution crystal structures of aminergic GPCRs, we review structural elements that contribute to and can be exploited for designing subtype-selective compounds. We describe the roles of secondary binding pockets (SBPs), as well as differences in ligand entry pathways to the orthosteric binding site, in determining selectivity. In addition, using the available crystal structures, we have identified conformational changes in the SBPs that are associated with receptor activation and explore the implications of these changes for the rational development of selective ligands with tailored efficacy.

Applying DEKOIS 2.0 in structure-based virtual screening to probe the impact of preparation procedures and score normalization

Background

Structure-based virtual screening techniques can help to identify new lead structures and complement other screening approaches in drug discovery. Prior to docking, the data (protein crystal structures and ligands) should be prepared with great attention to molecular and chemical details.

Results

Using a subset of 18 diverse targets from the recently introduced DEKOIS 2.0 benchmark set library, we found differences in the virtual screening performance of two popular docking tools (GOLD and Glide) when employing two different commercial packages (e.g. MOE and Maestro) for preparing input data. We systematically investigated the possible factors that can be responsible for the found differences in selected sets. For the Angiotensin-I-converting enzyme dataset, preparation of the bioactive molecules clearly exerted the highest influence on VS performance compared to preparation of the decoys or the target structure. The major contributing factors were different protonation states, molecular flexibility, and differences in the input conformation (particularly for cyclic moieties) of bioactives. In addition, score normalization strategies eliminated the biased docking scores shown by GOLD (ChemPLP) for the larger bioactives and produced a better performance. Generalizing these normalization strategies on the 18 DEKOIS 2.0 sets, improved the performances for the majority of GOLD (ChemPLP) docking, while it showed detrimental performances for the majority of Glide (SP) docking.

Conclusions

In conclusion, we exemplify herein possible issues particularly during the preparation stage of molecular data and demonstrate to which extent these issues can cause perturbations in the virtual screening performance. We provide insights into what problems can occur and should be avoided, when generating benchmarks to characterize the virtual screening performance. Particularly, careful selection of an appropriate molecular preparation setup for the bioactive set and the use of score normalization for docking with GOLD (ChemPLP) appear to have a great importance for the screening performance. For virtual screening campaigns, we recommend to invest time and effort into including alternative preparation workflows into the generation of the master library, even at the cost of including multiple representations of each molecule.

Graphical Abstract

Advancements in therapeutically targeting orphan GPCRs

G-protein coupled receptors (GPCRs) are popular biological targets for drug discovery and development. To date there are more than 140 orphan GPCRs, i.e., receptors whose endogenous ligands are unknown. Traditionally orphan GPCRs have been difficult to study and the development of therapeutic compounds targeting these receptors has been extremely slow although these GPCRs are considered important targets based on their distribution and behavioral phenotype as revealed by animals lacking the receptor. Recent advances in several methods used to study orphan receptors, including protein crystallography and homology modeling are likely to be useful in the identification of therapeutics targeting these receptors. In the past 13 years, over a dozen different Class A GPCRs have been crystallized; this trend is exciting, since homology modeling of GPCRs has previously been limited by the availability of solved structures. As the number of solved GPCR structures continues to grow so does the number of templates that can be used to generate increasingly accurate models of phylogenetically related orphan GPCRs. The availability of solved structures along with the advances in using multiple templates to build models (in combination with molecular dynamics simulations that reveal structural information not provided by crystallographic data and methods for modeling hard-to-predict flexible loop regions) have improved the quality of GPCR homology models. This, in turn, has improved the success rates of virtual ligand screens that use homology models to identify potential receptor binding compounds. Experimental testing of the predicted hits and validation using traditional GPCR pharmacological approaches can be used to drive ligand-based efforts to probe orphan receptor biology as well as to define the chemotypes and chemical scaffolds important for binding. As a result of these advances, orphan GPCRs are emerging from relative obscurity as a new class of drug targets.

Hierarchical virtual screening approaches in small molecule drug discovery

Virtual screening has played a significant role in the discovery of small molecule inhibitors of therapeutic targets in last two decades. Various ligand and structure-based virtual screening approaches are employed to identify small molecule ligands for proteins of interest. These approaches are often combined in either hierarchical or parallel manner to take advantage of the strength and avoid the limitations associated with individual methods. Hierarchical combination of ligand and structure-based virtual screening approaches has received noteworthy success in numerous drug discovery campaigns. In hierarchical virtual screening, several filters using ligand and structure-based approaches are sequentially applied to reduce a large screening library to a number small enough for experimental testing. In this review, we focus on different hierarchical virtual screening strategies and their application in the discovery of small molecule modulators of important drug targets. Several virtual screening studies are discussed to demonstrate the successful application of hierarchical virtual screening in small molecule drug discovery.

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.

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.

Histamine H4 receptor ligands: future applications and state of art

Histamine is a chemical transmitter found practically in whole organism and exerts its effects through the interaction with H1 to H4 histaminergic receptors. Specifically, H4 receptors are found mainly in immune cells and blood-forming tissues, thus are involved in inflammatory and immune processes, as well as some actions in central nervous system. Therefore, H4 receptor ligands can have applications in the treatment of chronic inflammatory and immune diseases and may be novel therapeutic option in these conditions. Several H4 receptor ligands have been described from early 2000's until nowadays, being imidazole, indolecarboxamide, 2-aminopyrimidine, quinazoline, and quinoxaline scaffolds the most explored and discussed in this review. Moreover, several studies of molecular modeling using homology models of H4 receptor and QSAR data of the ligands are summarized. The increasing and promising therapeutic applications are leading these compounds to clinical trials, which probably will be part of the next generation of blockbuster drugs.

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.

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.

Identification of Novel Histamine H4 Ligands by Virtual Screening on Molecular Dynamics Ensembles

We report the identification of novel histamine H4 receptor ligands by ensemble docking on homology model conformers derived from molecular dynamics simulations. Selected receptor models from the trajectories demonstrated superior virtual screening performance compared to the initial models. The ensemble of the best models was able to retrieve a diverse set of known H4 ligands. Prospective virtual screening against these models and subsequent in vitro experimental validation identified novel H4 ligands. Compound 3 showing highest affinity and ligand efficiency represents an interesting scaffold for further medicinal chemistry exploration.

Comparative Modeling of Drug Target Proteins☆

In this perspective, we begin by describing the comparative protein structure modeling technique and the accuracy of the corresponding models. We then discuss the significant role that comparative prediction plays in drug discovery. We focus on virtual ligand screening against comparative models and illustrate the state-of-the-art by a number of specific examples.

The impact of molecular dynamics sampling on the performance of virtual screening against GPCRs

The formation of ligand-protein complexes requires simultaneous adaptation of the binding partners. In structure based virtual screening, high throughput docking approaches typically consider the ligand flexibility, but the conformational freedom of the protein is usually taken into account in a limited way. The goal of this study is to elaborate a methodology for incorporating protein flexibility to improve the virtual screening enrichments on GPCRs. Explicit-solvated molecular dynamics simulations (MD) were carried out in lipid bilayers to generate an ensemble of protein conformations for the X-ray structures and homology models of both aminergic and peptidergic GPCRs including the chemokine CXCR4, dopamine D3, histamine H4, and serotonin 5HT6 holo receptor complexes. The quality of the receptor models was assessed by enrichment studies to compare X-ray structures, homology models, and snapshots from the MD trajectory. According to our results, selected frames from the MD trajectory can outperform X-ray structures and homology models in terms of enrichment factor and AUC values. Significant changes were observed considering EF1% values: comparing the original CXCR4, D3, and H4 targets and the additional 5HT6 initial models to that of the best MD frame resulted in 0 to 6.7, 0.32 to 3.5 (10×), 13.3 to 26.7 (2×), and 0 to 14.1 improvements, respectively. It is worth noting that rank-average based ensemble evaluation calculated for different ensemble sizes could not improve the results further. We propose here that MD simulation can capture protein conformations representing the key interacting points of the receptor but less biased toward one specific chemotype. These conformations are useful for the identification of a "consensus" binding site with improved performance in virtual screening.

Ligand pose and orientational sampling in molecular docking

Molecular docking remains an important tool for structure-based screening to find new ligands and chemical probes. As docking ambitions grow to include new scoring function terms, and to address ever more targets, the reliability and extendability of the orientation sampling, and the throughput of the method, become pressing. Here we explore sampling techniques that eliminate stochastic behavior in DOCK3.6, allowing us to optimize the method for regularly variable sampling of orientations. This also enabled a focused effort to optimize the code for efficiency, with a three-fold increase in the speed of the program. This, in turn, facilitated extensive testing of the method on the 102 targets, 22,805 ligands and 1,411,214 decoys of the Directory of Useful Decoys - Enhanced (DUD-E) benchmarking set, at multiple levels of sampling. Encouragingly, we observe that as sampling increases from 50 to 500 to 2000 to 5000 to 20000 molecular orientations in the binding site (and so from about 1×10¹⁰ to 4×10¹⁰ to 1×10¹¹ to 2×10¹¹ to 5×10¹¹ mean atoms scored per target, since multiple conformations are sampled per orientation), the enrichment of ligands over decoys monotonically increases for most DUD-E targets. Meanwhile, including internal electrostatics in the evaluation ligand conformational energies, and restricting aromatic hydroxyls to low energy rotamers, further improved enrichment values. Several of the strategies used here to improve the efficiency of the code are broadly applicable in the field.

Optimized method of G-protein-coupled receptor homology modeling: its application to the discovery of novel CXCR7 ligands

Homology modeling of G-protein-coupled seven-transmembrane receptors (GPCRs) remains a challenge despite the increasing number of released GPCR crystal structures. This challenge can be attributed to the low sequence identity and structural diversity of the ligand-binding pocket of GPCRs. We have developed an optimized GPCR structure modeling method based on multiple GPCR crystal structures. This method was designed to be applicable to distantly related receptors of known structural templates. CXC chemokine receptor (CXCR7) is a potential drug target for cancer chemotherapy. Homology modeling, docking, and virtual screening for CXCR7 were carried out using our method. The predicted docking poses of the known antagonists were different from the crystal structure of human CXCR4 with the small-molecule antagonist IT1t. Furthermore, 21 novel CXCR7 ligands with IC50 values of 1.29-11.4 μM with various scaffolds were identified by structure-based virtual screening.

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.

Mapping histamine H4 receptor–ligand binding modes

Computational prediction of ligand binding modes in G protein-coupled receptors (GPCRs) remains a challenging task. Systematic consideration of different protein modelling templates, ligand binding poses, and ligand protonation states in extensive molecular dynamics (MD) simulation studies enabled the prediction of ligand-specific mutation effects in the histamine H₄ receptor, a key player in inflammation.

Structural modelling and dynamics of proteins for insights into drug interactions

Proteins are the workhorses of biomolecules and their function is affected by their structure and their structural rearrangements during ligand entry, ligand binding and protein-protein interactions. Hence, the knowledge of protein structure and, importantly, the dynamic behaviour of the structure are critical for understanding how the protein performs its function. The predictions of the structure and the dynamic behaviour can be performed by combinations of structure modelling and molecular dynamics simulations. The simulations also need to be sensitive to the constraints of the environment in which the protein resides. Standard computational methods now exist in this field to support the experimental effort of solving protein structures. This review presents a comprehensive overview of the basis of the calculations and the well-established computational methods used to generate and understand protein structure and function and the study of their dynamic behaviour with the reference to lung-related targets.

Virtual ligand screening against comparative protein structure models

Virtual ligand screening uses computation to discover new ligands of a protein by screening one or more of its structural models against a database of potential ligands. Comparative protein structure modeling extends the applicability of virtual screening beyond the atomic structures determined by X-ray crystallography or NMR spectroscopy. Here, we describe an integrated modeling and docking protocol, combining comparative modeling by MODELLER and virtual ligand screening by DOCK.

Discovery of novel histamine H4 and serotonin transporter ligands using the topological feature tree descriptor

Ligand-based approaches are particularly important in the hit identification process of drug discovery when no structural information on the target is available. Pharmacophore descriptors that use a topological representation of the ligands are usually fast enough to screen large compound libraries effectively when seeking novel lead candidates. One example of this kind is the Feature Tree descriptor, a reduced graph representation implemented in the FTrees software. In this study, we tested the screening efficiency of FTrees by both retrospective and prospective screens using known histamine H4 antagonists and serotonin transporter (SERT) inhibitors as query molecules. Our results demonstrate that FTrees can effectively find actives. Particularly when combined with a subsequent 2D fingerprint-based diversity selection, FTrees was found to be extremely effective at discovering a diverse set of scaffolds. Prospective screening of our in-house compound deck provided several novel H4 and SERT ligands that could serve as suitable starting points for further optimization.

Structure-based prediction of subtype selectivity of histamine H3 receptor selective antagonists in clinical trials

Histamine receptors (HRs) are excellent drug targets for the treatment of diseases, such as schizophrenia, psychosis, depression, migraine, allergies, asthma, ulcers, and hypertension. Among them, the human H(3) histamine receptor (hH(3)HR) antagonists have been proposed for specific therapeutic applications, including treatment of Alzheimer's disease, attention deficit hyperactivity disorder (ADHD), epilepsy, and obesity. However, many of these drug candidates cause undesired side effects through the cross-reactivity with other histamine receptor subtypes. In order to develop improved selectivity and activity for such treatments, it would be useful to have the three-dimensional structures for all four HRs. We report here the predicted structures of four HR subtypes (H(1), H(2), H(3), and H(4)) using the GEnSeMBLE (GPCR ensemble of structures in membrane bilayer environment) Monte Carlo protocol, sampling ∼35 million combinations of helix packings to predict the 10 most stable packings for each of the four subtypes. Then we used these 10 best protein structures with the DarwinDock Monte Carlo protocol to sample ∼50 000 × 10(20) poses to predict the optimum ligand-protein structures for various agonists and antagonists. We find that E206(5.46) contributes most in binding H(3) selective agonists (5, 6, 7) in agreement with experimental mutation studies. We also find that conserved E5.46/S5.43 in both of hH(3)HR and hH(4)HR are involved in H(3)/ H(4) subtype selectivity. In addition, we find that M378(6.55) in hH(3)HR provides additional hydrophobic interactions different from hH(4)HR (the corresponding amino acid of T323(6.55) in hH(4)HR) to provide additional subtype bias. From these studies, we developed a pharmacophore model based on our predictions for known hH(3)HR selective antagonists in clinical study [ABT-239 1, GSK-189,254 2, PF-3654746 3, and BF2.649 (tiprolisant) 4] that suggests critical selectivity directing elements are: the basic proton interacting with D114(3.32), the spacer, the aromatic ring substituted with the hydrophilic or lipophilic groups interacting with lipophilic pockets in transmembranes (TMs) 3-5-6 and the aliphatic ring located in TMs 2-3-7. These 3D structures for all four HRs should help guide the rational design of novel drugs for the subtype selective antagonists and agonists with reduced side effects.

Crystal structure-based virtual screening for fragment-like ligands of the human histamine H(1) receptor

The recent crystal structure determinations of druggable class A G protein-coupled receptors (GPCRs) have opened up excellent opportunities in structure-based ligand discovery for this pharmaceutically important protein family. We have developed and validated a customized structure-based virtual fragment screening protocol against the recently determined human histamine H(1) receptor (H(1)R) crystal structure. The method combines molecular docking simulations with a protein-ligand interaction fingerprint (IFP) scoring method. The optimized in silico screening approach was successfully applied to identify a chemically diverse set of novel fragment-like (≤22 heavy atoms) H(1)R ligands with an exceptionally high hit rate of 73%. Of the 26 tested fragments, 19 compounds had affinities ranging from 10 μM to 6 nM. The current study shows the potential of in silico screening against GPCR crystal structures to explore novel, fragment-like GPCR ligand space.