Search

Chemical tools to monitor drug-target engagement of endogenously expressed protein kinases are highly desirable for preclinical target validation in drug discovery. Here, we describe a chemical genetics strategy to selectively study target engagement of endogenous kinases. By substituting a serine residue into cysteine at the DFG-1 position in the ATP-binding pocket, we sensitize the non-receptor tyrosine kinase FES towards covalent labeling by a complementary fluorescent chemical probe. This mutation is introduced in the endogenous FES gene of HL-60 cells using CRISPR/Cas9 gene editing. Leveraging the temporal and acute control offered by our strategy, we show that FES activity is dispensable for differentiation of HL-60 cells towards macrophages. Instead, FES plays a key role in neutrophil phagocytosis via SYK kinase activation. This chemical genetics strategy holds promise as a target validation method for kinases., Chemical tools to monitor drug-target engagement of endogenous enzymes are essential for preclinical target validation. Here, the authors present a chemical genetics strategy to study target engagement of endogenous kinases, achieving specific labeling and inactivation of FES kinase to provide insights into FES’ role in neutrophil phagocytosis.

The recent crystal structures of CC chemokine receptors 2 and 9 (CCR2 and CCR9) have provided structural evidence for an allosteric, intracellular binding site. The high conservation of residues involved in this site suggests its presence in most chemokine receptors, including the close homologue CCR1. By using [H-3]CCR2-RA-[R], a high-affinity, CCR2 intracellular ligand, we report an intracellular binding site in CCR1, where this radioligand also binds with high affinity. In addition, we report the synthesis and biological characterization of a series of pyrrolone derivatives for CCR1 and CCR2, which allowed us to identify several high-affinity intracellular ligands, including selective and potential multitarget antagonists. Evaluation of selected compounds in a functional [S-35]GTP gamma S assay revealed that they act as inverse agonists in CCR1, providing a new manner of pharmacological modulation. Thus, this intracellular binding site enables the design of selective and multitarget inhibitors as a novel therapeutic approach.

In this work, we present a case study to explore the challenges associated with finding novel molecules for a receptor that has been studied in depth and has a wealth of chemical information available. Specifically, we apply a previously described protocol that incorporates explicit water molecules in the ligand binding site to prospectively screen over 2.5 million drug-like and lead-like compounds from the commercially available eMolecules database in search of novel binders to the adenosine A2A receptor (A2AAR). A total of seventy-one compounds were selected for purchase and biochemical assaying based on high ligand efficiency and high novelty (Tanimoto coefficient ≤0.25 to any A2AAR tested compound). These molecules were then tested for their affinity to the adenosine A2A receptor in a radioligand binding assay. We identified two hits that fulfilled the criterion of ~50 % radioligand displacement at a concentration of 10 μM. Next we selected an additional eight novel molecules that were predicted to make a bidentate interaction with Asn2536.55, a key interacting residue in the binding pocket of the A2AAR. None of these eight molecules were found to be active. Based on these results we discuss the advantages of structure-based methods and the challenges associated with finding chemically novel molecules for well-explored targets. Electronic supplementary material The online version of this article (doi:10.1007/s10822-016-9963-7) contains supplementary material, which is available to authorized users.

The rapid growth of structural information for G-protein-coupled receptors (GPCRs) has led to a greater understanding of their structure, function, selectivity, and ligand binding. Although novel ligands have been identified using methods such as virtual screening, computationally driven lead optimization has been possible only in isolated cases because of challenges associated with predicting binding free energies for related compounds. Here, we provide a systematic characterization of the performance of free-energy perturbation (FEP) calculations to predict relative binding free energies of congeneric ligands binding to GPCR targets using a consistent protocol and no adjustable parameters. Using the FEP+ package, first we validated the protocol, which includes a full lipid bilayer and explicit solvent, by predicting the binding affinity for a total of 45 different ligands across four different GPCRs (adenosine A2AAR, β1 adrenergic, CXCR4 chemokine, and δ opioid receptors). Comparison with experimental binding affinity measurements revealed a highly predictive ranking correlation (average spearman ρ = 0.55) and low root-mean-square error (0.80 kcal/mol). Next, we applied FEP+ in a prospective project, where we predicted the affinity of novel, potent adenosine A2A receptor (A2AR) antagonists. Four novel compounds were synthesized and tested in a radioligand displacement assay, yielding affinity values in the nanomolar range. The affinity of two out of the four novel ligands (plus three previously reported compounds) was correctly predicted (within 1 kcal/mol), including one compound with approximately a tenfold increase in affinity compared to the starting compound. Detailed analyses of the simulations underlying the predictions provided insights into the structural basis for the two cases where the affinity was overpredicted. Taken together, these results establish a protocol for systematically applying FEP+ to GPCRs and provide guidelines for identifying potent molecules in drug discovery lead optimization projects.

The sodium ion site is an allosteric site conserved among many G protein-coupled receptors (GPCRs). Amiloride 1 and 5-(N,N-hexamethylene)amiloride 2 (HMA) supposedly bind in this sodium ion site and can influence orthosteric ligand binding. The availability of a high-resolution X-ray crystal structure of the human adenosine A2A receptor (hA2AAR), in which the allosteric sodium ion site was elucidated, makes it an appropriate model receptor for investigating the allosteric site. In this study, we report the synthesis and evaluation of novel 5'-substituted amiloride derivatives as hA2AAR allosteric antagonists. The potency of the amiloride derivatives was assessed by their ability to displace orthosteric radioligand [(3)H]4-(2-((7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a]-[1,3,5]triazin-5-yl)amino)ethyl)phenol ([(3)H]ZM-241,385) from both the wild-type and sodium ion site W246A mutant hA2AAR. 4-Ethoxyphenethyl-substituted amiloride 12l was found to be more potent than both amiloride and HMA, and the shift in potency between the wild-type and mutated receptor confirmed its likely binding to the sodium ion site.

Recent crystal structures of multiple G protein-coupled receptors (GPCRs) have revealed a highly conserved intracellular pocket that can be used to modulate these receptors from the inside. This novel intracellular site partially overlaps with the G protein and β-arrestin binding site, providing a new manner of pharmacological intervention. Here we provide an update of the architecture and function of the intracellular region of GPCRs, until now portrayed as the signaling domain. We review the available evidence on the presence of intracellular binding sites among chemokine receptors and other class A GPCRs, as well as different strategies to target it, including small molecules, pepducins, and nanobodies. Finally, the potential advantages of intracellular (allosteric) ligands over orthosteric ligands are also discussed.

We report on the synthesis and biological evaluation of a series of 1,2-diarylimidazol-4-carboxamide derivatives developed as CB1 receptor antagonists. These were evaluated in a radioligand displacement binding assay, a [35S]GTPγS binding assay, and in a competition association assay that enables the relatively fast kinetic screening of multiple compounds. The compounds show high affinities and a diverse range of kinetic profiles at the CB1 receptor and their structure–kinetic relationships (SKRs) were established. Using the recently resolved hCB1 receptor crystal structures, we also performed a modeling study that sheds light on the crucial interactions for both the affinity and dissociation kinetics of this family of ligands. We provide evidence that, next to affinity, additional knowledge of binding kinetics is useful for selecting new hCB1 receptor antagonists in the early phases of drug discovery.

We expanded on a series of pyrido[2,1-f]purine-2,4-dione derivatives as human adenosine A3 receptor (hA3R) antagonists to determine their kinetic profiles and affinities. Many compounds showed high affinities and a diverse range of kinetic profiles. We found hA3R antagonists with very short residence time (RT) at the receptor (2.2 min for 5) and much longer RTs (e.g., 376 min for 27 or 391 min for 31). Two representative antagonists (5 and 27) were tested in [35S]GTPγS binding assays, and their RTs appeared correlated to their (in)surmountable antagonism. From a kon-koff-KD kinetic map, we divided the antagonists into three subgroups, providing a possible direction for the further development of hA3R antagonists. Additionally, we performed a computational modeling study that sheds light on the crucial receptor interactions, dictating the compounds' binding kinetics. Knowledge of target binding kinetics appears useful for developing and triaging new hA3R antagonists in the early phase of drug discovery.

The increase of publicly available bioactivity data in recent years has fueled and catalyzed research in chemogenomics, data mining, and modeling approaches. As a direct result, over the past few years a multitude of different methods have been reported and evaluated, such as target fishing, nearest neighbor similarity-based methods, and Quantitative Structure Activity Relationship (QSAR)-based protocols. However, such studies are typically conducted on different datasets, using different validation strategies, and different metrics. In this study, different methods were compared using one single standardized dataset obtained from ChEMBL, which is made available to the public, using standardized metrics (BEDROC and Matthews Correlation Coefficient). Specifically, the performance of Naïve Bayes, Random Forests, Support Vector Machines, Logistic Regression, and Deep Neural Networks was assessed using QSAR and proteochemometric (PCM) methods. All methods were validated using both a random split validation and a temporal validation, with the latter being a more realistic benchmark of expected prospective execution. Deep Neural Networks are the top performing classifiers, highlighting the added value of Deep Neural Networks over other more conventional methods. Moreover, the best method (‘DNN_PCM’) performed significantly better at almost one standard deviation higher than the mean performance. Furthermore, Multi-task and PCM implementations were shown to improve performance over single task Deep Neural Networks. Conversely, target prediction performed almost two standard deviations under the mean performance. Random Forests, Support Vector Machines, and Logistic Regression performed around mean performance. Finally, using an ensemble of DNNs, alongside additional tuning, enhanced the relative performance by another 27% (compared with unoptimized ‘DNN_PCM’). Here, a standardized set to test and evaluate different machine learning algorithms in the context of multi-task learning is offered by providing the data and the protocols.Graphical Abstract. Electronic supplementary material The online version of this article (doi:10.1186/s13321-017-0232-0) contains supplementary material, which is available to authorized users.

The structure of the human A2A adenosine receptor has been elucidated by X-ray crystallography with a high affinity non-xanthine antagonist, ZM241385, bound to it. This template molecule served as a starting point for the incorporation of reactive moieties that cause the ligand to covalently bind to the receptor. In particular, we incorporated a fluorosulfonyl moiety onto ZM241385, which yielded LUF7445 (4-((3-((7-amino-2-(furan-2-yl)-[1, 2, 4]triazolo[1,5-a][1, 3, 5]triazin-5-yl)amino)propyl)carbamoyl)benzene sulfonyl fluoride). In a radioligand binding assay, LUF7445 acted as a potent antagonist, with an apparent affinity for the hA2A receptor in the nanomolar range. Its apparent affinity increased with longer incubation time, suggesting an increasing level of covalent binding over time. An in silico A2A-structure-based docking model was used to study the binding mode of LUF7445. This led us to perform site-directed mutagenesis of the A2A receptor to probe and validate the target lysine amino acid K153 for covalent binding. Meanwhile, a functional assay combined with wash-out experiments was set up to investigate the efficacy of covalent binding of LUF7445. All these experiments led us to conclude LUF7445 is a valuable molecular tool for further investigating covalent interactions at this receptor. It may also serve as a prototype for a therapeutic approach in which a covalent antagonist may be needed to counteract prolonged and persistent presence of the endogenous ligand adenosine.KEYWORDS: A2A adenosine receptor; Adenosine; Covalent antagonist; G protein-coupled receptors; Radioligand binding

A significant challenge and potential high-value application of computer-aided drug design is the accurate prediction of protein–ligand binding affinities. Free energy perturbation (FEP) using molecular dynamics (MD) sampling is among the most suitable approaches to achieve accurate binding free energy predictions, due to the rigorous statistical framework of the methodology, correct representation of the energetics, and thorough treatment of the important degrees of freedom in the system (including explicit waters). Recent advances in sampling methods and force fields coupled with vast increases in computational resources have made FEP a viable technology to drive hit-to-lead and lead optimization, allowing for more efficient cycles of medicinal chemistry and the possibility to explore much larger chemical spaces. However, previous FEP applications have focused on systems with high-resolution crystal structures of the target as starting points—something that is not always available in drug discovery projects. As such, the ability to apply FEP on homology models would greatly expand the domain of applicability of FEP in drug discovery. In this work we apply a particular implementation of FEP, called FEP+, on congeneric ligand series binding to four diverse targets: a kinase (Tyk2), an epigenetic bromodomain (BRD4), a transmembrane GPCR (A2A), and a protein–protein interaction interface (BCL-2 family protein MCL-1). We apply FEP+ using both crystal structures and homology models as starting points and find that the performance using homology models is generally on a par with the results when using crystal structures. The robustness of the calculations to structural variations in the input models can likely be attributed to the conformational sampling in the molecular dynamics simulations, which allows the modeled receptor to adapt to the “real” conformation for each ligand in the series. This work exemplifies the advantages of using all-atom simulation methods with full system flexibility and offers promise for the general application of FEP to homology models, although additional validation studies should be performed to further understand the limitations of the method and the scenarios where FEP will work best.

The expanding number of crystal structures of G protein-coupled receptors (GPCRs) has increased the knowledge on receptor function and their ability to recognize ligands. Although structure-based virtual screening has been quite successful on GPCRs, scores obtained by docking are typically not indicative for ligand affinity. Methods capturing interactions between protein and ligand in a more explicit manner, such as interaction fingerprints (IFPs), have been applied as an addition or alternative to docking. Originally IFPs captured the interactions of amino acid residues with ligands with specific definitions for the various interaction types. More complex IFPs now capture atom-atom interactions, such as in SYBYL, or fragment-fragment co-occurrences such as in SPLIF. Overall, most of the IFPs have been studied in comparison with docking in retrospective studies. For GPCRs it remains unclear which IFP should be used, if at all, and in what manner. Thus, the performance between five different IFPs was compared on five different representative GPCRs, including several extensions of the original implementations,. Results show that the more detailed IFPs, SYBYL and SPLIF, perform better than the other IFPs (Deng, Credo, and Elements). SPLIF was further tuned based on the number of poses, fingerprint similarity coefficient, and using an ensemble of structures. Enrichments were obtained that were significantly higher than initial enrichments and those obtained by 2D-similarity. With the increase in available crystal structures for GPCRs, and given that IFPs such as SPLIF enhance enrichment in virtual screens, it is anticipated that IFPs will be used in conjunction with docking, especially for GPCRs with a large binding pocket.

The adhesion G protein-coupled receptors [ADGRs/class B2 G protein-coupled receptors (GPCRs)] constitute an ancient family of GPCRs that have recently been demonstrated to play important roles in cellular and developmental processes. Here, we describe a first insight into the structure-function relationship of ADGRs using the family member ADGR subfamily G member 4 (ADGRG4)/GPR112 as a model receptor. In a bioinformatics approach, we compared conserved, functional elements of the well-characterized class A and class B1 secretin-like GPCRs with the ADGRs. We identified several potential equivalent motifs and subjected those to mutational analysis. The importance of the mutated residues was evaluated by examining their effect on the high constitutive activity of the N-terminally truncated ADGRG4/GPR112 in a 1-receptor-1-G protein Saccharomyces cerevisiae screening system and was further confirmed in a transfected mammalian human embryonic kidney 293 cell line. We evaluated the results in light of the crystal structures of the class A adenosine A2A receptor and the class B1 corticotropin-releasing factor receptor 1. ADGRG4 proved to have functionally important motifs resembling class A, class B, and combined elements, but also a unique highly conserved ADGR motif (H3.33). Given the high conservation of these motifs and residues across the adhesion GPCR family, it can be assumed that these are general elements of ADGR function.-Peeters, M. C., Mos, I., Lenselink, E. B., Lucchesi, M., IJzerman, A. P., Schwartz, T. W. Getting from A to B-exploring the activation motifs of the class B adhesion G protein-coupled receptor subfamily G member 4/GPR112.