Large-Scale Modeling of Sparse Protein Kinase Activity Data.

Protein kinases are a protein family that plays an important role in several complex diseases such as cancer and cardiovascular and immunological diseases. Protein kinases have conserved ATP binding sites, which when targeted can lead to similar activities of inhibitors against different kinases. This can be exploited to create multitarget drugs. On the other hand, selectivity (lack of similar activities) is desirable in order to avoid toxicity issues. There is a vast amount of protein kinase activity data in the public domain, which can be used in many different ways. Multitask machine learning models are expected to excel for these kinds of data sets because they can learn from implicit correlations between tasks (in this case activities against a variety of kinases). However, multitask modeling of sparse data poses two major challenges: (i) creating a balanced train-test split without data leakage and (ii) handling missing data. In this work, we construct a protein kinase benchmark set composed of two balanced splits without data leakage, using random and dissimilarity-driven cluster-based mechanisms, respectively. This data set can be used for benchmarking and developing protein kinase activity prediction models. Overall, the performance on the dissimilarity-driven cluster-based split is lower than on random split-based sets for all models, indicating poor generalizability of models. Nevertheless, we show that multitask deep learning models, on this very sparse data set, outperform single-task deep learning and tree-based models. Finally, we demonstrate that data imputation does not improve the performance of (multitask) models on this benchmark set.

Collaborative SAR Modeling and Prospective In Vitro Validation of Oxidative Stress Activation in Human HepG2 Cells.

Oxidative stress is the consequence of an abnormal increase of reactive oxygen species (ROS). ROS are generated mainly during the metabolism in both normal and pathological conditions as well as from exposure to xenobiotics. Xenobiotics can, on the one hand, disrupt molecular machinery involved in redox processes and, on the other hand, reduce the effectiveness of the antioxidant activity. Such dysregulation may lead to oxidative damage when combined with oxidative stress overpassing the cell capacity to detoxify ROS. In this work, a green fluorescent protein (GFP)-tagged nuclear factor erythroid 2-related factor 2 (NRF2)-regulated sulfiredoxin reporter (Srxn1-GFP) was used to measure the antioxidant response of HepG2 cells to a large series of drug and drug-like compounds (2230 compounds). These compounds were then classified as positive or negative depending on cellular response and distributed among different modeling groups to establish structure-activity relationship (SAR) models. A selection of models was used to prospectively predict oxidative stress induced by a new set of compounds subsequently experimentally tested to validate the model predictions. Altogether, this exercise exemplifies the different challenges of developing SAR models of a phenotypic cellular readout, model combination, chemical space selection, and results interpretation.

Enhancing Hit Discovery in Virtual Screening through Absolute Protein-Ligand Binding Free-Energy Calculations.

In the hit identification stage of drug discovery, a diverse chemical space needs to be explored to identify initial hits. Contrary to empirical scoring functions, absolute protein-ligand binding free-energy perturbation (ABFEP) provides a theoretically more rigorous and accurate description of protein-ligand binding thermodynamics and could, in principle, greatly improve the hit rates in virtual screening. In this work, we describe an implementation of an accurate and reliable ABFEP method in FEP+. We validated the ABFEP method on eight congeneric compound series binding to eight protein receptors including both neutral and charged ligands. For ligands with net charges, the alchemical ion approach is adopted to avoid artifacts in electrostatic potential energy calculations. The calculated binding free energies correlate with experimental results with a weighted average of R2 = 0.55 for the entire dataset. We also observe an overall root-mean-square error (RMSE) of 1.1 kcal/mol after shifting the zero-point of the simulation data to match the average experimental values. Through ABFEP calculations using apo versus holo protein structures, we demonstrated that the protein conformational and protonation state changes between the apo and holo proteins are the main physical factors contributing to the protein reorganization free energy manifested by the overestimation of raw ABFEP calculated binding free energies using the holo structures of the proteins. Furthermore, we performed ABFEP calculations in three virtual screening applications for hit enrichment. ABFEP greatly improves the hit rates as compared to docking scores or other methods like metadynamics. The good performance of ABFEP in rank ordering compounds demonstrated in this work confirms it as a useful tool to improve the hit rates in virtual screening, thus facilitating hit discovery.

Multitask machine learning models for predicting lipophilicity (logP) in the SAMPL7 challenge.

Accurate prediction of lipophilicity-logP-based on molecular structures is a well-established field. Predictions of logP are often used to drive forward drug discovery projects. Driven by the SAMPL7 challenge, in this manuscript we describe the steps that were taken to construct a novel machine learning model that can predict and generalize well. This model is based on the recently described Directed-Message Passing Neural Networks (D-MPNNs). Further enhancements included: both the inclusion of additional datasets from ChEMBL (RMSE improvement of 0.03), and the addition of helper tasks (RMSE improvement of 0.04). To the best of our knowledge, the concept of adding predictions from other models (Simulations Plus logP and logD@pH7.4, respectively) as helper tasks is novel and could be applied in a broader context. The final model that we constructed and used to participate in the challenge ranked 2/17 ranked submissions with an RMSE of 0.66, and an MAE of 0.48 (submission: Chemprop). On other datasets the model also works well, especially retrospectively applied to the SAMPL6 challenge where it would have ranked number one out of all submissions (RMSE of 0.35). Despite the fact that our model works well, we conclude with suggestions that are expected to improve the model even further.

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation.

The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand-protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (R2) compared with experiment in the range of 0.6-0.9 and mean unsigned errors compared to experiment of 0.6-0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.

Successive Statistical and Structure-Based Modeling to Identify Chemically Novel Kinase Inhibitors.

Kinases are frequently studied in the context of anticancer drugs. Their involvement in cell responses, such as proliferation, differentiation, and apoptosis, makes them interesting subjects in multitarget drug design. In this study, a workflow is presented that models the bioactivity spectra for two panels of kinases: (1) inhibition of RET, BRAF, SRC, and S6K, while avoiding inhibition of MKNK1, TTK, ERK8, PDK1, and PAK3, and (2) inhibition of AURKA, PAK1, FGFR1, and LKB1, while avoiding inhibition of PAK3, TAK1, and PIK3CA. Both statistical and structure-based models were included, which were thoroughly benchmarked and optimized. A virtual screening was performed to test the workflow for one of the main targets, RET kinase. This resulted in 5 novel and chemically dissimilar RET inhibitors with remaining RET activity of <60% (at a concentration of 10 µM) and similarities with known RET inhibitors from 0.18 to 0.29 (Tanimoto, ECFP6). The four more potent inhibitors were assessed in a concentration range and proved to be modestly active with a pIC50 value of 5.1 for the most active compound. The experimental validation of inhibitors for RET strongly indicates that the multitarget workflow is able to detect novel inhibitors for kinases, and hence, this workflow can potentially be applied in polypharmacology modeling. We conclude that this approach can identify new chemical matter for existing targets. Moreover, this workflow can easily be applied to other targets as well.

Drug Discovery Maps, a Machine Learning Model That Visualizes and Predicts Kinome-Inhibitor Interaction Landscapes.

The interpretation of high-dimensional structure-activity data sets in drug discovery to predict ligand-protein interaction landscapes is a challenging task. Here we present Drug Discovery Maps (DDM), a machine learning model that maps the activity profile of compounds across an entire protein family, as illustrated here for the kinase family. DDM is based on the t-distributed stochastic neighbor embedding (t-SNE) algorithm to generate a visualization of molecular and biological similarity. DDM maps chemical and target space and predicts the activities of novel kinase inhibitors across the kinome. The model was validated using independent data sets and in a prospective experimental setting, where DDM predicted new inhibitors for FMS-like tyrosine kinase 3 (FLT3), a therapeutic target for the treatment of acute myeloid leukemia. Compounds were resynthesized, yielding highly potent, cellularly active FLT3 inhibitors. Biochemical assays confirmed most of the predicted off-targets. DDM is further unique in that it is completely open-source and available as a ready-to-use executable to facilitate broad and easy adoption.

Comprehensive structure-activity-relationship of azaindoles as highly potent FLT3 inhibitors.

Acute myeloid leukemia (AML) is characterized by fast progression and low survival rates, in which Fms-like tyrosine kinase 3 (FLT3) receptor mutations have been identified as a driver mutation in cancer progression in a subgroup of AML patients. Clinical trials have shown emergence of drug resistant mutants, emphasizing the ongoing need for new chemical matter to enable the treatment of this disease. Here, we present the discovery and topological structure-activity relationship (SAR) study of analogs of isoquinolinesulfonamide H-89, a well-known PKA inhibitor, as FLT3 inhibitors. Surprisingly, we found that the SAR was not consistent with the observed binding mode of H-89 in PKA. Matched molecular pair analysis resulted in the identification of highly active sub-nanomolar azaindoles as novel FLT3-inhibitors. Structure based modelling using the FLT3 crystal structure suggested an alternative, flipped binding orientation of the new inhibitors.

Getting from A to B-exploring the activation motifs of the class B adhesion G protein-coupled receptor subfamily G member 4/GPR112.

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.

A covalent antagonist for the human adenosine A2A receptor.

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.

Relative Binding Free Energy Calculations Applied to Protein Homology Models.

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.

Interacting with GPCRs: Using Interaction Fingerprints for Virtual Screening.

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.

In search of novel ligands using a structure-based approach: a case study on the adenosine A2A receptor.

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.

Sodium ion binding pocket mutations and adenosine A2A receptor function.

Recently we identified a sodium ion binding pocket in a high-resolution structure of the human adenosine A2A receptor. In the present study we explored this binding site through site-directed mutagenesis and molecular dynamics simulations. Amino acids in the pocket were mutated to alanine, and their influence on agonist and antagonist affinity, allosterism by sodium ions and amilorides, and receptor functionality was explored. Mutation of the polar residues in the Na(+) pocket were shown to either abrogate (D52A(2.50) and N284A(7.49)) or reduce (S91A(3.39), W246A(6.48), and N280A(7.45)) the negative allosteric effect of sodium ions on agonist binding. Mutations D52A(2.50) and N284A(7.49) completely abolished receptor signaling, whereas mutations S91A(3.39) and N280A(7.45) elevated basal activity and mutations S91A(3.39), W246A(6.48), and N280A(7.45) decreased agonist-stimulated receptor signaling. In molecular dynamics simulations D52A(2.50) directly affected the mobility of sodium ions, which readily migrated to another pocket formed by Glu13(1.39) and His278(7.43). The D52A(2.50) mutation also decreased the potency of amiloride with respect to ligand displacement but did not change orthosteric ligand affinity. In contrast, W246A(6.48) increased some of the allosteric effects of sodium ions and amiloride, whereas orthosteric ligand binding was decreased. These new findings suggest that the sodium ion in the allosteric binding pocket not only impacts ligand affinity but also plays a vital role in receptor signaling. Because the sodium ion binding pocket is highly conserved in other class A G protein-coupled receptors, our findings may have a general relevance for these receptors and may guide the design of novel synthetic allosteric modulators or bitopic ligands.

Discovery and mapping of an intracellular antagonist binding site at the chemokine receptor CCR2.

The chemokine receptor CCR2 is a G protein-coupled receptor that is involved in many diseases characterized by chronic inflammation, and therefore a large variety of CCR2 small molecule antagonists has been developed. On the basis of their chemical structures these antagonists can roughly be divided into two groups with most likely two topographically distinct binding sites. The aim of the current study was to identify the binding site of one such group of ligands, exemplified by three allosteric antagonists, CCR2-RA-[R], JNJ-27141491, and SD-24. We first used a chimeric CCR2/CCR5 receptor approach to obtain insight into the binding site of the allosteric antagonists and additionally introduced eight single point mutations in CCR2 to further characterize the putative binding pocket. All constructs were studied in radioligand binding and/or functional IP turnover assays, providing evidence for an intracellular binding site for CCR2-RA-[R], JNJ-27141491, and SD-24. For CCR2-RA-[R] the most important residues for binding were found to be the highly conserved tyrosine Y(7.53) and phenylalanine F(8.50) of the NPxxYx(5,6)F motif, as well as V(6.36) at the bottom of TM-VI and K(8.49) in helix-VIII. These findings demonstrate for the first time the presence of an allosteric intracellular binding site for CCR2 antagonists. This contributes to an increased understanding of the interactions of diverse ligands at CCR2 and may allow for a more rational design of future allosteric antagonists.

A yeast screening method to decipher the interaction between the adenosine A2B receptor and the C-terminus of different G protein α-subunits.

The expression of human G protein-coupled receptors (GPCRs) in Saccharomyces cerevisiae containing chimeric yeast/mammalian Gα subunits provides a useful tool for the study of GPCR activation. In this study, we used a one-GPCR-one-G protein yeast screening method in combination with molecular modeling and mutagenesis studies to decipher the interaction between GPCRs and the C-terminus of different α-subunits of G proteins. We chose the human adenosine A2B receptor (hA2BR) as a paradigm, a typical class A GPCR that shows promiscuous behavior in G protein coupling in this yeast system. The wild-type hA2BR and five mutant receptors were expressed in 8 yeast strains with different humanized G proteins, covering the four major classes: Gαi, Gαs, Gαq, and Gα12. Our experiments showed that a tyrosine residue (Y) at the C-terminus of the Gα subunit plays an important role in controlling the activation of GPCRs. Receptor residues R103(3.50) and I107(3.54) are vital too in G protein-coupling and the activation of the hA2BR, whereas L213(IL3) is more important in G protein inactivation. Substitution of S235(6.36) to alanine provided the most divergent G protein-coupling profile. Finally, L236(6.37) substitution decreased receptor activation in all G protein pathways, although to a different extent. In conclusion, our findings shed light on the selectivity of receptor/G protein coupling, which may help in further understanding GPCR signaling.

Selecting an optimal number of binding site waters to improve virtual screening enrichments against the adenosine A2A receptor.

A major challenge in structure-based virtual screening (VS) involves the treatment of explicit water molecules during docking in order to improve the enrichment of active compounds over decoys. Here we have investigated this in the context of the adenosine A2A receptor, where water molecules have previously been shown to be important for achieving high enrichment rates with docking, and where the positions of some binding site waters are known from a high-resolution crystal structure. The effect of these waters (both their presence and orientations) on VS enrichment was assessed using a carefully curated set of 299 high affinity A2A antagonists and 17,337 decoys. We show that including certain crystal waters greatly improves VS enrichment and that optimization of water hydrogen positions is needed in order to achieve the best results. We also show that waters derived from a molecular dynamics simulation - without any knowledge of crystallographic waters - can improve enrichments to a similar degree as the crystallographic waters, which makes this strategy applicable to structures without experimental knowledge of water positions. Finally, we used decision trees to select an ensemble of structures with different water molecule positions and orientations that outperforms any single structure with water molecules. The approach presented here is validated against independent test sets of A2A receptor antagonists and decoys from the literature. In general, this water optimization strategy could be applied to any target with waters-mediated protein-ligand interactions.

Chemical genetics strategy to profile kinase target engagement reveals role of FES in neutrophil phagocytosis.

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.

Synthesis and Pharmacological Evaluation of Triazolopyrimidinone Derivatives as Noncompetitive, Intracellular Antagonists for CC Chemokine Receptors 2 and 5.

CC chemokine receptors 2 (CCR2) and 5 (CCR5) are involved in many inflammatory diseases; however, most CCR2 and CCR5 clinical candidates have been unsuccessful. (Pre)clinical evidence suggests that dual CCR2/CCR5 inhibition might be more effective in the treatment of such multifactorial diseases. In this regard, the highly conserved intracellular binding site in chemokine receptors provides a new avenue for the design of multitarget ligands. In this study, we synthesized and evaluated the biological activity of a series of triazolopyrimidinone derivatives in CCR2 and CCR5. Radioligand binding assays first showed that they bind to the intracellular site of CCR2, and in combination with functional assays on CCR5, we explored structure-affinity/activity relationships in both receptors. Although most compounds were CCR2-selective, 39 and 43 inhibited ß-arrestin recruitment in CCR5 with high potency. Moreover, these compounds displayed an insurmountable mechanism of inhibition in both receptors, which holds promise for improved efficacy in inflammatory diseases.

Development of Covalent Ligands for G Protein-Coupled Receptors: A Case for the Human Adenosine A3 Receptor.

The development of covalent ligands for G protein-coupled receptors (GPCRs) is not a trivial process. Here, we report a streamlined workflow thereto from synthesis to validation, exemplified by the discovery of a covalent antagonist for the human adenosine A3 receptor (hA3AR). Based on the 1 H,3 H-pyrido[2,1- f]purine-2,4-dione scaffold, a series of ligands bearing a fluorosulfonyl warhead and a varying linker was synthesized. This series was subjected to an affinity screen, revealing compound 17b as the most potent antagonist. In addition, a nonreactive methylsulfonyl derivative 19 was developed as a reversible control compound. A series of assays, comprising time-dependent affinity determination, washout experiments, and [35S]GTPγS binding assays, then validated 17b as the covalent antagonist. A combined in silico hA3AR-homology model and site-directed mutagenesis study was performed to demonstrate that amino acid residue Y2657.36 was the unique anchor point of the covalent interaction. This workflow might be applied to other GPCRs to guide the discovery of covalent ligands.