Potent and selective TYK2-JH1 inhibitors highly efficacious in rodent model of psoriasis.

TYK2 is a member of the JAK family of kinases and a key mediator of IL-12, IL-23, and type I interferon signaling. These cytokines have been implicated in the pathogenesis of multiple inflammatory and autoimmune diseases such as psoriasis, rheumatoid arthritis, lupus, and inflammatory bowel diseases. Supported by compelling data from human genetic association studies, TYK2 inhibition is an attractive therapeutic strategy for these diseases. Herein, we report the discovery of a series of highly selective catalytic site TYK2 inhibitors designed using FEP+ and structurally enabled design starting from a virtual screen hit. We highlight the structure-based optimization to identify a lead candidate 30, a potent cellular TYK2 inhibitor with excellent selectivity, pharmacokinetic properties, and in vivo efficacy in a mouse psoriasis model.

Acetyl-CoA carboxylase inhibition by ND-630 reduces hepatic steatosis, improves insulin sensitivity, and modulates dyslipidemia in rats.

Simultaneous inhibition of the acetyl-CoA carboxylase (ACC) isozymes ACC1 and ACC2 results in concomitant inhibition of fatty acid synthesis and stimulation of fatty acid oxidation and may favorably affect the morbidity and mortality associated with obesity, diabetes, and fatty liver disease. Using structure-based drug design, we have identified a series of potent allosteric protein-protein interaction inhibitors, exemplified by ND-630, that interact within the ACC phosphopeptide acceptor and dimerization site to prevent dimerization and inhibit the enzymatic activity of both ACC isozymes, reduce fatty acid synthesis and stimulate fatty acid oxidation in cultured cells and in animals, and exhibit favorable drug-like properties. When administered chronically to rats with diet-induced obesity, ND-630 reduces hepatic steatosis, improves insulin sensitivity, reduces weight gain without affecting food intake, and favorably affects dyslipidemia. When administered chronically to Zucker diabetic fatty rats, ND-630 reduces hepatic steatosis, improves glucose-stimulated insulin secretion, and reduces hemoglobin A1c (0.9% reduction). Together, these data suggest that ACC inhibition by representatives of this series may be useful in treating a variety of metabolic disorders, including metabolic syndrome, type 2 diabetes mellitus, and fatty liver disease.

Anti-leukaemic activity of the TYK2 selective inhibitor NDI-031301 in T-cell acute lymphoblastic leukaemia.

Activation of tyrosine kinase 2 (TYK2) contributes to the aberrant survival of T-cell acute lymphoblastic leukaemia (T-ALL) cells. Here we demonstrate the anti-leukaemic activity of a novel TYK2 inhibitor, NDI-031301. NDI-031301 is a potent and selective inhibitor of TYK2 that induced robust growth inhibition of human T-ALL cell lines. NDI-031301 treatment of human T-ALL cell lines resulted in induction of apoptosis that was not observed with the JAK inhibitors tofacitinib and baricitinib. Further investigation revealed that NDI-031301 treatment uniquely leads to activation of three mitogen-activated protein kinases (MAPKs), resulting in phosphorylation of ERK, SAPK/JNK and p38 MAPK coincident with PARP cleavage. Activation of p38 MAPK occurred within 1 h of NDI-031301 treatment and was responsible for NDI-031301-induced T-ALL cell death, as pharmacological inhibition of p38 MAPK partially rescued apoptosis induced by TYK2 inhibitor. Finally, daily oral administration of NDI-031301 at 100 mg/kg bid to immunodeficient mice engrafted with KOPT-K1 T-ALL cells was well tolerated, and led to decreased tumour burden and a significant survival benefit. These results support selective inhibition of TYK2 as a promising potential therapeutic strategy for T-ALL.

Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field.

Designing tight-binding ligands is a primary objective of small-molecule drug discovery. Over the past few decades, free-energy calculations have benefited from improved force fields and sampling algorithms, as well as the advent of low-cost parallel computing. However, it has proven to be challenging to reliably achieve the level of accuracy that would be needed to guide lead optimization (∼5× in binding affinity) for a wide range of ligands and protein targets. Not surprisingly, widespread commercial application of free-energy simulations has been limited due to the lack of large-scale validation coupled with the technical challenges traditionally associated with running these types of calculations. Here, we report an approach that achieves an unprecedented level of accuracy across a broad range of target classes and ligands, with retrospective results encompassing 200 ligands and a wide variety of chemical perturbations, many of which involve significant changes in ligand chemical structures. In addition, we have applied the method in prospective drug discovery projects and found a significant improvement in the quality of the compounds synthesized that have been predicted to be potent. Compounds predicted to be potent by this approach have a substantial reduction in false positives relative to compounds synthesized on the basis of other computational or medicinal chemistry approaches. Furthermore, the results are consistent with those obtained from our retrospective studies, demonstrating the robustness and broad range of applicability of this approach, which can be used to drive decisions in lead optimization.

Crystal structure and pharmacological characterization of a novel N-methyl-D-aspartate (NMDA) receptor antagonist at the GluN1 glycine binding site.

NMDA receptors are ligand-gated ion channels that mediate excitatory neurotransmission in the brain. They are tetrameric complexes composed of glycine-binding GluN1 and GluN3 subunits together with glutamate-binding GluN2 subunits. Subunit-selective antagonists that discriminate between the glycine sites of GluN1 and GluN3 subunits would be valuable pharmacological tools for studies on the function and physiological roles of NMDA receptor subtypes. In a virtual screening for antagonists that exploit differences in the orthosteric binding site of GluN1 and GluN3 subunits, we identified a novel glycine site antagonist, 1-thioxo-1,2-dihydro-[1,2,4]triazolo[4,3-a]quinoxalin-4(5H)-one (TK40). Here, we show by Schild analysis that TK40 is a potent competitive antagonist with Kb values of 21-63 nM at the GluN1 glycine-binding site of the four recombinant GluN1/N2A-D receptors. In addition, TK40 displayed >100-fold selectivity for GluN1/N2 NMDA receptors over GluN3A- and GluN3B-containing NMDA receptors and no appreciable effects at AMPA receptors. Binding experiments on rat brain membranes and the purified GluN1 ligand-binding domain using glycine site GluN1 radioligands further confirmed the competitive interaction and high potency. To delineate the binding mechanism, we have solved the crystal structure of the GluN1 ligand-binding domain in complex with TK40 and show that TK40 binds to the orthosteric binding site of the GluN1 subunit with a binding mode that was also predicted by virtual screening. Furthermore, the structure reveals that the imino acetamido group of TK40 acts as an α-amino acid bioisostere, which could be of importance in bioisosteric replacement strategies for future ligand design.

Docking covalent inhibitors: a parameter free approach to pose prediction and scoring.

Although many popular docking programs include a facility to account for covalent ligands, large-scale systematic docking validation studies of covalent inhibitors have been sparse. In this paper, we present the development and validation of a novel approach for docking and scoring covalent inhibitors, which consists of conventional noncovalent docking, heuristic formation of the covalent attachment point, and structural refinement of the protein-ligand complex. This approach combines the strengths of the docking program Glide and the protein structure modeling program Prime and does not require any parameter fitting for the study of additional covalent reaction types. We first test this method by predicting the native binding geometry of 38 covalently bound complexes. The average RMSD of the predicted poses is 1.52 Å, and 76% of test set inhibitors have an RMSD of less than 2.0 Å. In addition, the apparent affinity score constructed herein is tested on a virtual screening study and the characterization of the SAR properties of two different series of congeneric compounds with satisfactory success.

Discovery of a Potent and Selective Tyrosine Kinase 2 Inhibitor: TAK-279.

TYK2 is a key mediator of IL12, IL23, and type I interferon signaling, and these cytokines have been implicated in the pathogenesis of multiple inflammatory and autoimmune diseases such as psoriasis, rheumatoid arthritis, lupus, and inflammatory bowel diseases. Supported by compelling data from human genome-wide association studies and clinical results, TYK2 inhibition through small molecules is an attractive therapeutic strategy to treat these diseases. Herein, we report the discovery of a series of highly selective pseudokinase (Janus homology 2, JH2) domain inhibitors of TYK2 enzymatic activity. A computationally enabled design strategy, including the use of FEP+, was instrumental in identifying a pyrazolo-pyrimidine core. We highlight the utility of computational physics-based predictions used to optimize this series of molecules to identify the development candidate 30, a potent, exquisitely selective cellular TYK2 inhibitor that is currently in Phase 2 clinical trials for the treatment of psoriasis and psoriatic arthritis.

Docking performance of the glide program as evaluated on the Astex and DUD datasets: a complete set of glide SP results and selected results for a new scoring function integrating WaterMap and glide.

Glide SP mode enrichment results for two preparations of the DUD dataset and native ligand docking RMSDs for two preparations of the Astex dataset are presented. Following a best-practices preparation scheme, an average RMSD of 1.140 Å for native ligand docking with Glide SP is computed. Following the same best-practices preparation scheme for the DUD dataset an average area under the ROC curve (AUC) of 0.80 and average early enrichment via the ROC (0.1 %) metric of 0.12 were observed. 74 and 56 % of the 39 best-practices prepared targets showed AUC over 0.7 and 0.8, respectively. Average AUC was greater than 0.7 for all best-practices protein families demonstrating consistent enrichment performance across a broad range of proteins and ligand chemotypes. In both Astex and DUD datasets, docking performance is significantly improved employing a best-practices preparation scheme over using minimally-prepared structures from the PDB. Enrichment results for WScore, a new scoring function and sampling methodology integrating WaterMap and Glide, are presented for four DUD targets, hivrt, hsp90, cdk2, and fxa. WScore performance in early enrichment is consistently strong and all systems examined show AUC > 0.9 and superior early enrichment to DUD best-practices Glide SP results.

Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution.

Generating the appropriate protonation states of drug-like molecules in solution is important for success in both ligand- and structure-based virtual screening. Screening collections of millions of compounds requires a method for determining tautomers and their energies that is sufficiently rapid, accurate, and comprehensive. To maximise enrichment, the lowest energy tautomers must be determined from heterogeneous input, without over-enumerating unfavourable states. While computationally expensive, the density functional theory (DFT) method M06-2X/aug-cc-pVTZ(-f) [PB-SCRF] provides accurate energies for enumerated model tautomeric systems. The empirical Hammett-Taft methodology can very rapidly extrapolate substituent effects from model systems to drug-like molecules via the relationship between pK(T) and pK(a). Combining the two complementary approaches transforms the tautomer problem from a scientific challenge to one of engineering scale-up, and avoids issues that arise due to the very limited number of measured pK(T) values, especially for the complicated heterocycles often favoured by medicinal chemists for their novelty and versatility. Several hundreds of pre-calculated tautomer energies and substituent pK(a) effects are tabulated in databases for use in structural adjustment by the program Epik, which treats tautomers as a subset of the larger problem of the protonation states in aqueous ensembles and their energy penalties. Accuracy and coverage is continually improved and expanded by parameterizing new systems of interest using DFT and experimental data. Recommendations are made for how to best incorporate tautomers in molecular design and virtual screening workflows.

Inhibition of Acetyl-CoA Carboxylase by Phosphorylation or the Inhibitor ND-654 Suppresses Lipogenesis and Hepatocellular Carcinoma.

The incidence of hepatocellular carcinoma (HCC) is rapidly increasing due to the prevalence of obesity and non-alcoholic fatty liver disease, but the molecular triggers that initiate disease development are not fully understood. We demonstrate that mice with targeted loss-of-function point mutations within the AMP-activated protein kinase (AMPK) phosphorylation sites on acetyl-CoA carboxylase 1 (ACC1 Ser79Ala) and ACC2 (ACC2 Ser212Ala) have increased liver de novo lipogenesis (DNL) and liver lesions. The same mutation in ACC1 also increases DNL and proliferation in human liver cancer cells. Consistent with these findings, a novel, liver-specific ACC inhibitor (ND-654) that mimics the effects of ACC phosphorylation inhibits hepatic DNL and the development of HCC, improving survival of tumor-bearing rats when used alone and in combination with the multi-kinase inhibitor sorafenib. These studies highlight the importance of DNL and dysregulation of AMPK-mediated ACC phosphorylation in accelerating HCC and the potential of ACC inhibitors for treatment.

Structures of the ligand-binding core of iGluR2 in complex with the agonists (R)- and (S)-2-amino-3-(4-hydroxy-1,2,5-thiadiazol-3-yl)propionic acid explain their unusual equipotency.

AMPA-type ionotropic glutamate receptors generally display high stereoselectivity in agonist binding. However, the stereoisomers of 2-amino-3-(4-hydroxy-1,2,5-thiadiazol-3-yl)propionic acid (TDPA) have similar enantiopharmacology. To understand this observation, we have determined the X-ray structures of ( R)-TDPA and ( S)-TDPA in complex with the ligand-binding core of iGluR2 and investigated the binding pharmacology at AMPA and kainate receptors. Both enantiomers induce full domain closure in iGluR2 but adopt different conformations when binding to the receptor, which may explain the similar enantiopharmacology.

Spatial scale, abundance and the species-energy relationship in British birds.

1. The spatial scale of analysis may influence the nature, strength and underlying drivers of macroecological patterns, one of the most frequently discussed of which is the relationship between species richness and environmental energy availability. 2. It has been suggested that species-energy relationships are hump-shaped at fine spatial grains and consistently positive at larger regional grains. The exact nature of this scale dependency is, however, the subject of much debate as relatively few studies have investigated species-energy relationships for the same assemblage across a range of spatial grains. Here, we contrast species-energy relationships for the British breeding avifauna at spatial grains of 1 km x 1 km, 2 km x 2 km and 10 km x 10 km plots, while maintaining a constant spatial extent. 3. Analyses were principally conducted using data on observed species richness. While survey work may fail to detect some species, observed species richness and that estimated using nonparametric techniques were strongly positively correlated with each other, and thus exhibit very similar spatial patterns. Moreover, the forms of species-energy relationships using observed and estimated species richness were statistically indistinguishable from each other. 4. Positive decelerating species-energy relationships arise at all three spatial grains. There is little evidence that the explanatory power of these relationships varies with spatial scale. However, ratios of regional (large-scale) to local (small-scale) species richness decrease with increasing energy availability, indicating that local richness responds to energy with a steeper gradient than does regional richness. Local assemblages thus sample a greater proportion of regional richness at higher energy levels, suggesting that spatial turnover of species richness is lower in high-energy regions. Similarly, a crude measure of temporal turnover, the ratio of cumulative species richness over a 4-year period to species richness in a single year, is lower in high-energy regions. These negative relationships between turnover and energy appear to be causal as both total and mean occupancy per species increases with energy. 5. While total density in 1 km x 1 km plots correlates positively with energy availability, such relationships are very weak for mean density per species. This suggests that the observed association between total abundance and species richness may not be mediated by population extinction rates, as predicted by the more individuals hypothesis. 6. The sampling mechanism suggests that species-energy relationships arise as high-energy areas support a greater number of individuals, and that random allocation of these individuals to local areas from a regional assemblage will generate species-energy relationships. While randomized local species-energy relationships are linear and positive, predicted richness is consistently greater than that observed. The mismatch between the observed and randomized species-energy relationships probably arises as a consequence of the aggregated nature of species distributions. The sampling mechanism, together with species spatial aggregation driven by limited habitat availability, may thus explain the species-energy relationship observed at this spatial scale.

The structure of a mixed GluR2 ligand-binding core dimer in complex with (S)-glutamate and the antagonist (S)-NS1209.

Ionotropic glutamate receptors (iGluRs) mediate fast synaptic transmission between cells of the central nervous system and are involved in various aspects of normal brain function. iGluRs are implicated in several brain disorders, e.g. in the high-frequency discharge of impulses during an epileptic seizure. (RS)-NS1209 functions as a competitive antagonist at 2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionate receptors, and shows robust preclinical anticonvulsant and neuroprotective effects. This study explores 2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionate receptor binding and selectivity of this novel class of antagonists. We present here the first X-ray structure of a mixed GluR2 ligand-binding core dimer, with the high-affinity antagonist (S)-8-methyl-5-(4-(N,N-dimethylsulfamoyl)phenyl)-6,7,8,9,-tetrahydro-1H-pyrrolo[3,2-h]-isoquinoline-2,3-dione-3-O-(4-hydroxybutyrate-2-yl)oxime [(S)-NS1209] in one protomer and the endogenous ligand (S)-glutamate in the other. (S)-NS1209 stabilises an even more open conformation of the D1 and D2 domains of the ligand-binding core than that of the apo structure due to steric hindrance. This is the first time ligand-induced hyperextension of the binding domains has been observed. (S)-NS1209 adopts a novel binding mode, including hydrogen bonding to Tyr450 and Gly451 of D1. Parts of (S)-NS1209 occupy new areas of the GluR2 ligand-binding cleft, and bind near residues that are not conserved among receptor subtypes. The affinities of (RS)-NS1209 at the GluR2 ligand-binding core as well as at GluR1-6 and mutated GluR1 and GluR3 receptors have been measured. Two distinct binding affinities were observed at the GluR3 and GluR4 receptors. In a functional in vitro assay, no difference in potency was observed between GluR2(Q)(o) and GluR3(o) receptors. The thermodynamics of binding of the antagonists (S)-NS1209, DNQX and (S)-ATPO to the GluR2 ligand-binding core have been determined by displacement isothermal titration calorimetry. The displacement of (S)-glutamate by all antagonists was shown to be driven by enthalpy.

Functional characterization of Tet-AMPA [tetrazolyl-2-amino-3-(3-hydroxy-5-methyl- 4-isoxazolyl)propionic acid] analogues at ionotropic glutamate receptors GluR1-GluR4. The molecular basis for the functional selectivity profile of 2-Bn-Tet-AMPA.

Four 2-substituted Tet-AMPA [Tet = tetrazolyl, AMPA = 2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid] analogues were characterized functionally at the homomeric AMPA receptors GluR1i, GluR2Qi, GluR3i, and GluR4i in a Fluo-4/Ca2+ assay. Whereas 2-Et-Tet-AMPA, 2-Pr-Tet-AMPA, and 2-iPr-Tet-AMPA were nonselective GluR agonists, 2-Bn-Tet-AMPA exhibited a 40-fold higher potency at GluR4i than at GluR1i. Examination of homology models of the S1-S2 domains of GluR1 and GluR4 containing 2-Bn-Tet-AMPA suggested four nonconserved residues in a region adjacent to the orthosteric site as possible determinants of the GluR4i/GluR1i selectivity. In a mutagenesis study, doubly mutating M686V/I687A in GluR1i in combination with either D399S or E683A increased both the potency and the maximal response of 2-Bn-Tet-AMPA at this receptor to levels similar to those elicited by the agonist at GluR4i. The dependence of the novel selectivity profile of 2-Bn-Tet-AMPA upon residues located outside of the orthosteric site underlines the potential for developing GluR subtype selective ligands by designing compounds with substituents that protrude beyond the (S)-Glu binding pocket.

A tetrazolyl-substituted subtype-selective AMPA receptor agonist.

Replacement of the methyl group of the AMPA receptor agonist 2-amino-3-[3-hydroxy-5-(2-methyl-2H-5-tetrazolyl)-4-isoxazolyl]propionic acid (2-Me-Tet-AMPA) with a benzyl group provided the first AMPA receptor agonist, compound 7, capable of discriminating GluR2-4 from GluR1 by its more than 10-fold preference for the former receptor subtypes. An X-ray crystallographic analysis of this new analogue in complex with the GluR2-S1S2J construct shows that accommodation of the benzyl group creates a previously unobserved pocket in the receptor, which may explain the remarkable pharmacological profile of compound 7.

An experiment on the impact of a neonicotinoid pesticide on honeybees: the value of a formal analysis of the data.

BACKGROUND: We assess the analysis of the data resulting from a field experiment conducted by Pilling et al. (PLoS ONE. doi: 10.1371/journal.pone.0077193, 5) on the potential effects of thiamethoxam on honeybees. The experiment had low levels of replication, so Pilling et al. concluded that formal statistical analysis would be misleading. This would be true if such an analysis merely comprised tests of statistical significance and if the investigators concluded that lack of significance meant little or no effect. However, an analysis that includes estimation of the size of any effects-with confidence limits-allows one to reach conclusions that are not misleading and that produce useful insights. MAIN BODY: For the data of Pilling et al., we use straightforward statistical analysis to show that the confidence limits are generally so wide that any effects of thiamethoxam could have been large without being statistically significant. Instead of formal analysis, Pilling et al. simply inspected the data and concluded that they provided no evidence of detrimental effects and from this that thiamethoxam poses a "low risk" to bees. CONCLUSIONS: Conclusions derived from the inspection of the data were not just misleading in this case but also are unacceptable in principle, for if data are inadequate for a formal analysis (or only good enough to provide estimates with wide confidence intervals), then they are bound to be inadequate as a basis for reaching any sound conclusions. Given that the data in this case are largely uninformative with respect to the treatment effect, any conclusions reached from such informal approaches can do little more than reflect the prior beliefs of those involved.

Accelerating drug discovery through tight integration of expert molecular design and predictive scoring.

Modeling protein-ligand interactions has been a central goal of computational chemistry for many years. We here review recent progress toward this goal, and highlight the role free energy calculation methods and computational solvent analysis techniques are now having in drug discovery. We further describe recent use of these methodologies to advance two separate drug discovery programs targeting acetyl-CoA carboxylase and tyrosine kinase 2. These examples suggest that tight integration of sophisticated chemistry teams with state-of-the-art computational methods can dramatically improve the efficiency of small molecule drug discovery.

A novel inhibitor of active protein kinase G attenuates chronic inflammatory and osteoarthritic pain.

Activating PKG-1α induces a long-term hyperexcitability (LTH) in nociceptive neurons. Since the LTH correlates directly with chronic pain in many animal models, we tested the hypothesis that inhibiting PKG-1α would attenuate LTH-mediated pain. We first synthesized and characterized compound N46 (N-((3R,4R)-4-(4-(2-fluoro-3-methoxy-6-propoxybenzoyl)benzamido)pyrrolidin-3-yl)-1H-indazole-5-carboxamide). N46 inhibits PKG-1α with an IC50 of 7.5 nmol, was highly selective when tested against a panel of 274 kinases, and tissue distribution studies indicate that it does not enter the CNS. To evaluate its antinociceptive potential, we used 2 animal models in which the pain involves both activated PKG-1α and LTH. Injecting complete Freund's adjuvant (CFA) into the rat hind paw causes a thermal hyperalgesia that was significantly attenuated 24 hours after a single intravenous injection of N46. Next, we used a rat model of osteoarthritic knee joint pain and found that a single intra-articular injection of N46 alleviated the pain 14 days after the pain was established and the relief lasted for 7 days. Thermal hyperalgesia and osteoarthritic pain are also associated with the activation of the capsaicin-activated transient receptor protein vanilloid-1 (TRPV1) channel. We show that capsaicin activates PKG-1α in nerves and that a subcutaneous delivery of N46 attenuated the mechanical and thermal hypersensitivity elicited by exposure to capsaicin. Thus, PKG-1α appears to be downstream of the transient receptor protein vanilloid-1. Our studies provide proof of concept in animal models that a PKG-1α antagonist has a powerful antinociceptive effect on persistent, already existing inflammatory pain. They further suggest that N46 is a valid chemotype for the further development of such antagonists.

3-hydroxypyridazine 1-oxides as carboxylate bioisosteres: a new series of subtype-selective AMPA receptor agonists.

Three positional isomers (compounds 1, 2, and 3) of 1-uracilylalanine (willardiine) based on a 3-hydroxypyridazine 1-oxide scaffold with an alanine side-chain at positions 4 (1), 5 (2) or 6 (3) were tested for binding to recombinant homomeric AMPA receptor (AMPA-R) subtypes GluR1-4, as well for excitatory activity on the rat cortical wedge preparation. 1 had approximately 30 times higher affinity than willardiine while showing a similar selectivity profile, i.e. 22-fold selectivity for GluR1/2 over GluR3/4. The GluR1-4 affinities of 3 were similar to 1, however, its 31-fold selectivity for GluR1/2 over GluR3/4 is the highest yet observed among azine-based glutamate analogues. The non-isosteric congener 2 showed weaker binding to AMPA-Rs. In the cortical wedge, 1 evokes similar responses to AMPA, while 3 and 2 are 10- and 100-fold weaker, respectively. Dose-response curves on Xenopus laevis oocytes expressing GluR1-4(flip) confirmed that 1 and 3 are potent GluR1/2 receptor agonists (EC(50)s from 0.26 to 1.7microM) but are 10- to 160-fold less potent at GluR3/4. The structures, potencies and selectivities of this new class of AMPA agonists are compared with those of willardiine, 5-fluorowillardiine and azawillardiine, referring to the binding mode observed in the crystal structure of willardiine bound to GluR2-S1S2. The results indicate that the 3-hydroxypyridazine 1-oxide moiety can function as an outstanding carboxylate mimic at AMPA-Rs, leading the way to further fine-tuning of subtype selectivity. This little-explored molecular motif may find wider application in medicinal chemistry.