Structure, function and pharmacology of human itch GPCRs.

The MRGPRX family of receptors (MRGPRX1-4) is a family of mas-related G-protein-coupled receptors that have evolved relatively recently1. Of these, MRGPRX2 and MRGPRX4 are key physiological and pathological mediators of itch and related mast cell-mediated hypersensitivity reactions2-5. MRGPRX2 couples to both Gi and Gq in mast cells6. Here we describe agonist-stabilized structures of MRGPRX2 coupled to Gi1 and Gq in ternary complexes with the endogenous peptide cortistatin-14 and with a synthetic agonist probe, respectively, and the development of potent antagonist probes for MRGPRX2. We also describe a specific MRGPRX4 agonist and the structure of this agonist in a complex with MRGPRX4 and Gq. Together, these findings should accelerate the structure-guided discovery of therapeutic agents for pain, itch and mast cell-mediated hypersensitivity.

Selectivity Challenges in Docking Screens for GPCR Targets and Antitargets.

To investigate large library docking's ability to find molecules with joint activity against on-targets and selectivity versus antitargets, the dopamine D2 and serotonin 5-HT2A receptors were targeted, seeking selectivity against the histamine H1 receptor. In a second campaign, κ-opioid receptor ligands were sought with selectivity versus the µ-opioid receptor. While hit rates ranged from 40% to 63% against the on-targets, they were just as good against the antitargets, even though the molecules were selected for their putative lack of binding to the off-targets. Affinities, too, were often as good or better for the off-targets. Even though it was occasionally possible to find selective molecules, such as a mid-nanomolar D2/5-HT2A ligand with 21-fold selectivity versus the H1 receptor, this was the exception. Whereas false-negatives are tolerable in docking screens against on-targets, they are intolerable against antitargets; addressing this problem may demand new strategies in the field.

Author Correction: Novel scaffolds for inhibition of Cruzipain identified from high-throughput screening of anti-kinetoplastid chemical boxes.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Structure-inspired design of ß-arrestin-biased ligands for aminergic GPCRs.

Development of biased ligands targeting G protein-coupled receptors (GPCRs) is a promising approach for current drug discovery. Although structure-based drug design of biased agonists remains challenging even with an abundance of GPCR crystal structures, we present an approach for translating GPCR structural data into ß-arrestin-biased ligands for aminergic GPCRs. We identified specific amino acid-ligand contacts at transmembrane helix 5 (TM5) and extracellular loop 2 (EL2) responsible for Gi/o and ß-arrestin signaling, respectively, and targeted those residues to develop biased ligands. For these ligands, we found that bias is conserved at other aminergic GPCRs that retain similar residues at TM5 and EL2. Our approach provides a template for generating arrestin-biased ligands by modifying predicted ligand interactions that block TM5 interactions and promote EL2 interactions. This strategy may facilitate the structure-guided design of arrestin-biased ligands at other GPCRs, including polypharmacological biased ligands.

Novel scaffolds for inhibition of Cruzipain identified from high-throughput screening of anti-kinetoplastid chemical boxes.

American Trypanosomiasis or Chagas disease is a prevalent, neglected and serious debilitating illness caused by the kinetoplastid protozoan parasite Trypanosoma cruzi. The current chemotherapy is limited only to nifurtimox and benznidazole, two drugs that have poor efficacy in the chronic phase and are rather toxic. In this scenario, more efficacious and safer drugs, preferentially acting through a different mechanism of action and directed against novel targets, are particularly welcome. Cruzipain, the main papain-like cysteine peptidase of T. cruzi, is an important virulence factor and a chemotherapeutic target with excellent pre-clinical validation evidence. Here, we present the identification of new Cruzipain inhibitory scaffolds within the GlaxoSmithKline HAT (Human African Trypanosomiasis) and Chagas chemical boxes, two collections grouping 404 non-cytotoxic compounds with high antiparasitic potency, drug-likeness, structural diversity and scientific novelty. We have adapted a continuous enzymatic assay to a medium-throughput format and carried out a primary screening of both collections, followed by construction and analysis of dose-response curves of the most promising hits. Using the identified compounds as a starting point a substructure directed search against CHEMBL Database revealed plausible common scaffolds while docking experiments predicted binding poses and specific interactions between Cruzipain and the novel inhibitors.

In silico design of novel probes for the atypical opioid receptor MRGPRX2.

The primate-exclusive MRGPRX2 G protein-coupled receptor (GPCR) has been suggested to modulate pain and itch. Despite putative peptide and small-molecule MRGPRX2 agonists, selective nanomolar-potency probes have not yet been reported. To identify a MRGPRX2 probe, we first screened 5,695 small molecules and found that many opioid compounds activated MRGPRX2, including (-)- and (+)-morphine, hydrocodone, sinomenine, dextromethorphan, and the prodynorphin-derived peptides dynorphin A, dynorphin B, and α- and ß-neoendorphin. We used these to select for mutagenesis-validated homology models and docked almost 4 million small molecules. From this docking, we predicted ZINC-3573-a potent MRGPRX2-selective agonist, showing little activity against 315 other GPCRs and 97 representative kinases-along with an essentially inactive enantiomer. ZINC-3573 activates endogenous MRGPRX2 in a human mast cell line, inducing degranulation and calcium release. MRGPRX2 is a unique atypical opioid-like receptor important for modulating mast cell degranulation, which can now be specifically modulated with ZINC-3573.

Allosteric ligands for the pharmacologically dark receptors GPR68 and GPR65.

At least 120 non-olfactory G-protein-coupled receptors in the human genome are 'orphans' for which endogenous ligands are unknown, and many have no selective ligands, hindering the determination of their biological functions and clinical relevance. Among these is GPR68, a proton receptor that lacks small molecule modulators for probing its biology. Using yeast-based screens against GPR68, here we identify the benzodiazepine drug lorazepam as a non-selective GPR68 positive allosteric modulator. More than 3,000 GPR68 homology models were refined to recognize lorazepam in a putative allosteric site. Docking 3.1 million molecules predicted new GPR68 modulators, many of which were confirmed in functional assays. One potent GPR68 modulator, ogerin, suppressed recall in fear conditioning in wild-type but not in GPR68-knockout mice. The same approach led to the discovery of allosteric agonists and negative allosteric modulators for GPR65. Combining physical and structure-based screening may be broadly useful for ligand discovery for understudied and orphan GPCRs.

Small-Molecule Allosteric Modulators of the Protein Kinase PDK1 from Structure-Based Docking.

Finding small molecules that target allosteric sites remains a grand challenge for ligand discovery. In the protein kinase field, only a handful of highly selective allosteric modulators have been found. Thus, more general methods are needed to discover allosteric modulators for additional kinases. Here, we use virtual screening against an ensemble of both crystal structures and comparative models to identify ligands for an allosteric peptide-binding site on the protein kinase PDK1 (the PIF pocket). We optimized these ligands through an analog-by-catalog search that yielded compound 4, which binds to PDK1 with 8 µM affinity. We confirmed the docking poses by determining a crystal structure of PDK1 in complex with 4. Because the PIF pocket appears to be a recurring structural feature of the kinase fold, known generally as the helix αC patch, this approach may enable the discovery of allosteric modulators for other kinases.

Quantification of the transferability of a designed protein specificity switch reveals extensive epistasis in molecular recognition.

Reengineering protein-protein recognition is an important route to dissecting and controlling complex interaction networks. Experimental approaches have used the strategy of "second-site suppressors," where a functional interaction is inferred between two proteins if a mutation in one protein can be compensated by a mutation in the second. Mimicking this strategy, computational design has been applied successfully to change protein recognition specificity by predicting such sets of compensatory mutations in protein-protein interfaces. To extend this approach, it would be advantageous to be able to "transplant" existing engineered and experimentally validated specificity changes to other homologous protein-protein complexes. Here, we test this strategy by designing a pair of mutations that modulates peptide recognition specificity in the Syntrophin PDZ domain, confirming the designed interaction biochemically and structurally, and then transplanting the mutations into the context of five related PDZ domain-peptide complexes. We find a wide range of energetic effects of identical mutations in structurally similar positions, revealing a dramatic context dependence (epistasis) of designed mutations in homologous protein-protein interactions. To better understand the structural basis of this context dependence, we apply a structure-based computational model that recapitulates these energetic effects and we use this model to make and validate forward predictions. Although the context dependence of these mutations is captured by computational predictions, our results both highlight the considerable difficulties in designing protein-protein interactions and provide challenging benchmark cases for the development of improved protein modeling and design methods that accurately account for the context.

Actin is required for IFT regulation in Chlamydomonas reinhardtii.

Assembly of cilia and flagella requires intraflagellar transport (IFT), a highly regulated kinesin-based transport system that moves cargo from the basal body to the tip of flagella [1]. The recruitment of IFT components to basal bodies is a function of flagellar length, with increased recruitment in rapidly growing short flagella [2]. The molecular pathways regulating IFT are largely a mystery. Because actin network disruption leads to changes in ciliary length and number, actin has been proposed to have a role in ciliary assembly. However, the mechanisms involved are unknown. In Chlamydomonas reinhardtii, conventional actin is found in both the cell body and the inner dynein arm complexes within flagella [3, 4]. Previous work showed that treating Chlamydomonas cells with the actin-depolymerizing compound cytochalasin D resulted in reversible flagellar shortening [5], but how actin is related to flagellar length or assembly remains unknown. Here we utilize small-molecule inhibitors and genetic mutants to analyze the role of actin dynamics in flagellar assembly in Chlamydomonas reinhardtii. We demonstrate that actin plays a role in IFT recruitment to basal bodies during flagellar elongation and that when actin is perturbed, the normal dependence of IFT recruitment on flagellar length is lost. We also find that actin is required for sufficient entry of IFT material into flagella during assembly. These same effects are recapitulated with a myosin inhibitor, suggesting that actin may act via myosin in a pathway by which flagellar assembly is regulated by flagellar length.

Structure-based ligand discovery for the protein-protein interface of chemokine receptor CXCR4.

G-protein-coupled receptors (GPCRs) are key signaling molecules and are intensely studied. Whereas GPCRs recognizing small-molecules have been successfully targeted for drug discovery, protein-recognizing GPCRs, such as the chemokine receptors, claim few drugs or even useful small molecule reagents. This reflects both the difficulties that attend protein-protein interface inhibitor discovery, and the lack of structures for these targets. Imminent structure determination of chemokine receptor CXCR4 motivated docking screens for new ligands against a homology model and subsequently the crystal structure. More than 3 million molecules were docked against the model and then against the crystal structure; 24 and 23 high-scoring compounds from the respective screens were tested experimentally. Docking against the model yielded only one antagonist, which resembled known ligands and lacked specificity, whereas the crystal structure docking yielded four that were dissimilar to previously known scaffolds and apparently specific. Intriguingly, several were potent and relatively small, with IC(50) values as low as 306 nM, ligand efficiencies as high as 0.36, and with efficacy in cellular chemotaxis. The potency and efficiency of these molecules has few precedents among protein-protein interface inhibitors, and supports structure-based efforts to discover leads for chemokine GPCRs.

In silico analysis of the binding of agonists and blockers to the ß2-adrenergic receptor.

Activation of G protein-coupled receptors (GPCRs) is a complex phenomenon. Here, we applied Induced Fit Docking (IFD) in tandem with linear discriminant analysis (LDA) to generate hypotheses on the conformational changes induced to the ß(2)-adrenergic receptor by agonist binding, preliminary to the sequence of events that characterize activation of the receptor. This analysis, corroborated by a follow-up molecular dynamics study, suggested that agonists induce subtle movements to the fifth transmembrane domain (TM5) of the receptor. Furthermore, molecular dynamics also highlighted a correlation between movements of TM5 and the second extracellular loop (EL2), suggesting that freedom of motion of EL2 is required for the agonist-induced TM5 displacement. Importantly, we also showed that the IFD/LDA procedure can be used as a computational means to distinguish agonists from blockers on the basis of the differential conformational changes induced to the receptor. In particular, the two most predictive models obtained are based on the RMSD induced to Ser207 and on the counterclockwise rotation induced to TM5.

Design, synthesis, crystal structures, and antimicrobial activity of sulfonamide boronic acids as ß-lactamase inhibitors.

We investigated a series of sulfonamide boronic acids that resulted from the merging of two unrelated AmpC ß-lactamase inhibitor series. The new boronic acids differed in the replacement of the canonical carboxamide, found in all penicillin and cephalosporin antibiotics, with a sulfonamide. Surprisingly, these sulfonamides had a highly distinct structure-activity relationship from the previously explored carboxamides, high ligand efficiencies (up to 0.91), and K(i) values down to 25 nM and up to 23 times better for smaller analogues. Conversely, K(i) values were 10-20 times worse for larger molecules than in the carboxamide congener series. X-ray crystal structures (1.6-1.8 Å) of AmpC with three of the new sulfonamides suggest that this altered structure-activity relationship results from the different geometry and polarity of the sulfonamide versus the carboxamide. The most potent inhibitor reversed ß-lactamase-mediated resistance to third generation cephalosporins, lowering their minimum inhibitory concentrations up to 32-fold in cell culture.

Structural basis of G protein-coupled receptor-G protein interactions.

The interaction of G protein-coupled receptors (GPCRs) with heterotrimeric G proteins represents one of the most fundamental biological processes. However, the molecular architecture of the GPCR-G protein complex remains poorly defined. In the present study, we applied a comprehensive GPCR-G protein alpha subunit (Galpha) chemical cross-linking strategy to map a receptor-Galpha interface, both before and after agonist-induced receptor activation. Using the M(3) muscarinic acetylcholine receptor (M3R)-Galpha(q) system as a model system, we examined the ability of approximately 250 combinations of cysteine-substituted M3R and Galpha(q) proteins to undergo cross-link formation. We identified many specific M3R-Galpha(q) contact sites, in both the inactive and active receptor conformations, allowing us to draw conclusions regarding the basic architecture of the M3R-Galpha(q) interface and the nature of the conformational changes following receptor activation. As heterotrimeric G proteins as well as most GPCRs share a high degree of structural homology, our findings should be of broad general relevance.

Ligand and structure-based models for the prediction of ligand-receptor affinities and virtual screenings: Development and application to the beta(2)-adrenergic receptor.

In this study, we evaluated the applicability of ligand-based and structure-based models to quantitative affinity predictions and virtual screenings for ligands of the beta(2)-adrenergic receptor, a G protein-coupled receptor (GPCR). We also devised and evaluated a number of consensus models obtained through partial least square regressions, to combine the strengths of the individual components. In all cases, the bioactive conformation of each ligand was derived from molecular docking at the crystal structure of the receptor. We identified the most effective models applicable to the different scenarios, in the presence or in the absence of a training set. For ranking the affinity of closely related analogs when a training set is available, a ligand-based consensus model (LI-CM) seems to be the best choice, while the structure-based MM-GBSA score seems the best alternative in the absence of a training set. For virtual screening purposes, the structure-based MM-GBSA score was found to be the method of choice. Consensus models consistently had performances superior or close to those of the best individual components, and were endowed with a significantly increased robustness. Given multiple models with no a priori knowledge of their predictive capabilities, constructing a consensus model ensures results very close to those that the best model alone would have yielded.

Structural basis of the selectivity of the beta(2)-adrenergic receptor for fluorinated catecholamines.

The important and diverse biological functions of adrenergic receptors, a subclass of G protein-coupled receptors (GPCRs), have made the search for compounds that selectively stimulate or inhibit the activity of different adrenergic receptor subtypes an important area of medicinal chemistry. We previously synthesized 2-, 5-, and 6-fluoronorepinehprine (FNE) and 2-, 5-, and 6-fluoroepinephrine (FEPI) and found that 2FNE and 2FEPI were selective beta-adrenergic agonists and that 6FNE and 6FEPI were selective alpha-adrenergic agonists, while 5FNE and 5FEPI were unselective. Agonist potencies correlated well with receptor binding affinities. Here, through a combination of molecular modeling and site-directed mutagenesis, we have identified N293 in the beta(2)-adrenergic receptor as a crucial residue for the selectivity of the receptor for catecholamines fluorinated at different positions.