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.

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.

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.

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.

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.

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.

Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes.

A novel scoring function to estimate protein-ligand binding affinities has been developed and implemented as the Glide 4.0 XP scoring function and docking protocol. In addition to unique water desolvation energy terms, protein-ligand structural motifs leading to enhanced binding affinity are included: (1) hydrophobic enclosure where groups of lipophilic ligand atoms are enclosed on opposite faces by lipophilic protein atoms, (2) neutral-neutral single or correlated hydrogen bonds in a hydrophobically enclosed environment, and (3) five categories of charged-charged hydrogen bonds. The XP scoring function and docking protocol have been developed to reproduce experimental binding affinities for a set of 198 complexes (RMSDs of 2.26 and 1.73 kcal/mol over all and well-docked ligands, respectively) and to yield quality enrichments for a set of fifteen screens of pharmaceutical importance. Enrichment results demonstrate the importance of the novel XP molecular recognition and water scoring in separating active and inactive ligands and avoiding false positives.

Multiconformation, Density Functional Theory-Based pKa Prediction in Application to Large, Flexible Organic Molecules with Diverse Functional Groups.

We consider the conformational flexibility of molecules and its implications for micro- and macro-pKa. The corresponding formulas are derived and discussed against the background of a comprehensive scientific and algorithmic description of the latest version of our computer program Jaguar pKa, a density functional theory-based pKa predictor, which is now capable of acting on multiple conformations explicitly. Jaguar pKa is essentially a complex computational workflow incorporating research and technologies from the fields of cheminformatics, molecular mechanics, quantum mechanics, and implicit solvation models. The workflow also makes use of automatically applied empirical corrections which account for the systematic errors resulting from the neglect of explicit solvent interactions in the algorithm's implicit solvent model. Applications of our program to large, flexible organic molecules representing several classes of functional groups are shown, with a particular emphasis in illustrations laid on drug-like molecules. It is demonstrated that a combination of aggressive conformational search and an explicit consideration of multiple conformations nearly eliminates the dependence of results on the initially chosen conformation. In certain cases this leads to unprecedented accuracy, which is sufficient for distinguishing stereoisomers that have slightly different pKa values. An application of Jaguar pKa to proton sponges, the pKa of which are strongly influenced by steric effects, showcases the advantages that pKa predictors based on quantum mechanical calculations have over similar empirical programs.

WScore: A Flexible and Accurate Treatment of Explicit Water Molecules in Ligand-Receptor Docking.

We have developed a new methodology for protein-ligand docking and scoring, WScore, incorporating a flexible description of explicit water molecules. The locations and thermodynamics of the waters are derived from a WaterMap molecular dynamics simulation. The water structure is employed to provide an atomic level description of ligand and protein desolvation. WScore also contains a detailed model for localized ligand and protein strain energy and integrates an MM-GBSA scoring component with these terms to assess delocalized strain of the complex. Ensemble docking is used to take into account induced fit effects on the receptor conformation, and protein reorganization free energies are assigned via fitting to experimental data. The performance of the method is evaluated for pose prediction, rank ordering of self-docked complexes, and enrichment in virtual screening, using a large data set of PDB complexes and compared with the Glide SP and Glide XP models; significant improvements are obtained.

Exploring the GluR2 ligand-binding core in complex with the bicyclical AMPA analogue (S)-4-AHCP.

The X-ray structure of the ionotropic GluR2 ligand-binding core (GluR2-S1S2J) in complex with the bicyclical AMPA analogue (S)-2-amino-3-(3-hydroxy-7,8-dihydro-6H-cyclohepta[d]-4-isoxazolyl)propionic acid [(S)-4-AHCP] has been determined, as well as the binding pharmacology of this construct and of the full-length GluR2 receptor. (S)-4-AHCP binds with a glutamate-like binding mode and the ligand adopts two different conformations. The K(i) of (S)-4-AHCP at GluR2-S1S2J was determined to be 185 +/- 29 nM and at full-length GluR2(R)o it was 175 +/- 8 nM. (S)-4-AHCP appears to elicit partial agonism at GluR2 by inducing an intermediate degree of domain closure (17 degrees). Also, functionally (S)-4-AHCP has an efficacy of 0.38 at GluR2(Q)i, relative to (S)-glutamate. The proximity of bound (S)-4-AHCP to domain D2 prevents full D1-D2 domain closure, which is limited by steric repulsion, especially between Leu704 and the ligand.

Convergent synthesis and pharmacology of substituted tetrazolyl-2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid analogues.

The synthesis and pharmacological characterization of 1- and 2-alkyltetrazolyl analogues of (RS)-2-amino-3-[3-hydroxy-5-(2-methyl-2H-5-tetrazolyl)-4-isoxazolyl]propionic acid (2-Me-Tet-AMPA), a highly potent and selective agonist at AMPA receptors, are presented. A shorter and more convergent synthetic route than previously described, employing a new method for introducing the amino acid moiety, was developed for these derivatives. The 2-substituted isomers were selective agonists, and their activity correlated inversely with the size of the substituent. Structural explanations of the structure-activity relationship are provided.

Design, synthesis, and pharmacology of a highly subtype-selective GluR1/2 agonist, (RS)-2-amino-3-(4-chloro-3-hydroxy-5-isoxazolyl)propionic acid (Cl-HIBO).

On the basis of structural studies, chloro-homoibotenic acid (Cl-HIBO) was designed and synthesized. Cl-HIBO was characterized in binding and electrophysiology experiments on native and cloned subtypes of GluRs. Electrophysiological selectivities ranged from 275 to 1600 for GluR1/2 over GluR3/4. The potent AMPA receptor activity was strongly desensitizing and the neurotoxicity similar to AMPA. Thus, Cl-HIBO is the most subtype selective agonist reported to date on GluR1/2, and offers a new standard for selectively studying subtypes of AMPA receptors.

Selective agonists at group II metabotropic glutamate receptors: synthesis, stereochemistry, and molecular pharmacology of (S)- and (R)-2-amino-4-(4-hydroxy[1,2,5]thiadiazol-3-yl)butyric acid.

Homologation of analogues of the central excitatory neurotransmitter glutamic acid (Glu), in which the distal carboxy group has been bioisosterically replaced by acidic heterocyclic units, has previously provided subtype selective ligands for metabotropic Glu receptors (mGluRs). The (S)-form of the 1,2,5-thiadiazol-3-ol Glu analogue, 2-amino-3-(4-hydroxy[1,2,5]thiadiazol-3-yl)propionic acid (TDPA, 6), is an 2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid (AMPA) receptor agonist, which in addition stereospecifically activates group I mGluRs. We have now synthesized the (S)- and (R)-forms of 2-amino-4-(4-hydroxy[1,2,5]thiadiazol-3-yl)butyric acid (homo-TDPA, 7) and shown that whereas neither enantiomer interacts with AMPA receptors, (S)- and (R)-7 appear to be selective and equipotent agonists at group II mGluRs as represented by the mGluR2 subtype. The activities of (S)- and (R)-7 are rationalized by conformational analysis, comparison with the potent and specific group II mGluR agonist (-)-LY379268 [(-)-12], and docking to a homology model of mGluR2.

Competitive antagonism of AMPA receptors by ligands of different classes: crystal structure of ATPO bound to the GluR2 ligand-binding core, in comparison with DNQX.

Ionotropic glutamate receptors (iGluRs) constitute a family of ligand-gated ion channels that are essential for mediating fast synaptic transmission in the central nervous system. This study presents a high-resolution X-ray structure of the competitive antagonist (S)-2-amino-3-[5-tert-butyl-3-(phosphonomethoxy)-4-isoxazolyl]propionic acid (ATPO) in complex with the ligand-binding core of the receptor. Comparison with the only previous structure of the ligand-binding core in complex with an antagonist, 6,7-dinitro-2,3-quinoxalinedione (DNQX) (Armstrong, N.; Gouaux, E. Neuron 2000, 28, 165-181), reveals that ATPO and DNQX stabilize an open form of the ligand-binding core by different sets of interactions. Computational techniques are used to quantify the differences between these two ligands and to map the binding site. The isoxazole moiety of ATPO acts primarily as a spacer, and other scaffolds could potentially be used. Whereas agonists induce substantial domain closures compared to the apo structure, ATPO only induces minor conformational changes. These results are consistent with the hypothesis that domain closure is related to receptor activation. To facilitate the design of novel AMPA receptor antagonists, we present a modified model of the binding site that includes key residues involved in ligand recognition.

(S)-2-Amino-3-(3-hydroxy-7,8-dihydro-6H-cyclohepta[d]isoxazol-4-yl)propionic acid, a potent and selective agonist at the GluR5 subtype of ionotropic glutamate receptors. Synthesis, modeling, and molecular pharmacology.

We have previously described (RS)-2-amino-3-(3-hydroxy-7,8-dihydro-6H-cyclohepta[d]isoxazol-4-yl)propionic acid (4-AHCP) as a highly effective agonist at non-N-methyl-d-aspartate (non-NMDA) glutamate (Glu) receptors in vivo, which is more potent than (RS)-2-amino-3-(3-hydroxy-5-methylisoxazol-4-yl)propionic acid (AMPA) but inactive at NMDA receptors. However, 4-AHCP was found to be much weaker than AMPA as an inhibitor of [(3)H]AMPA binding and to have limited effect in a [(3)H]kainic acid binding assay using rat cortical membranes. To shed light on the mechanism(s) underlying this quite enigmatic pharmacological profile of 4-AHCP, we have now developed a synthesis of (S)-4-AHCP (6) and (R)-4-AHCP (7). At cloned metabotropic Glu receptors mGluR1alpha (group I), mGluR2 (group II), and mGluR4a (group III), neither 6 nor 7 showed significant agonist or antagonist effects. The stereoisomer 6, but not 7, activated cloned AMPA receptor subunits GluR1o, GluR3o, and GluR4o with EC(50) values in the range 4.5-15 microM and the coexpressed kainate-preferring subunits GluR6 + KA2 (EC(50) = 6.4 microM). Compound 6, but not 7, proved to be a very potent agonist (EC(50) = 0.13 microM) at the kainate-preferring GluR5 subunit, equipotent with (S)-2-amino-3-(5-tert-butyl-3-hydroxyisothiazol-4-yl)propionic acid [(S)-Thio-ATPA, 4] and almost 4 times more potent than (S)-2-amino-3-(5-tert-butyl-3-hydroxyisoxazol-4-yl)propionic acid [(S)-ATPA, 3]. Compound 6 thus represents a new structural class of GluR5 agonists. Molecular modeling and docking to a crystal structure of the extracellular binding domain of the AMPA subunit GluR2 has enabled identification of the probable active conformation and binding mode of 6. We are able to rationalize the observed selectivities by comparing the docking of 4 and 6 to subtype constructs, i.e., a crystal structure of the extracellular binding domain of GluR2 and a homology model of GluR5.

Novel 1-hydroxyazole bioisosteres of glutamic acid. Synthesis, protolytic properties, and pharmacology.

A number of 1-hydroxyazole derivatives were synthesized as bioisosteres of (S)-glutamic acid (Glu) and as analogues of the AMPA receptor agonist (R,S)-2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid (AMPA, 3b). All compounds were subjected to in vitro pharmacological studies, including a series of Glu receptor binding assays, uptake studies on native as well as cloned Glu uptake systems, and the electrophysiological rat cortical slice model. Compounds 7a,b, analogues of AMPA bearing a 1-hydroxy-5-pyrazolyl moiety as the distal carboxylic functionality, showed only moderate affinity for [3H]AMPA receptor binding sites (IC(50) = 2.7 +/- 0.4 microM and IC(50) = 2.6 +/- 0.6 microM, respectively), correlating with electrophysiological data from the rat cortical wedge model (EC(50) = 280 +/- 48 microM and EC(50) = 586 +/- 41 microM, respectively). 1-Hydroxy-1,2,3-triazol-5-yl analogues of AMPA, compounds 8a,b, showed high affinity for [3H]AMPA receptor binding sites (IC(50) = 0.15 +/- 0.03 microM and IC(50) = 0.13 +/- 0.02 microM, respectively). Electrophysiological data showed that compound 8a was devoid of activity in the rat cortical wedge model (EC(50) > 1000 microM), whereas the corresponding 4-methyl analogue 8b was a potent AMPA receptor agonist (EC(50) = 15 +/- 2 microM). In accordance with this disparity, compound 8a was found to inhibit synaptosomal [3H]D-aspartic acid uptake (IC(50) = 93 +/- 25 microM), as well as excitatory amino acid transporters (EAATs) EAAT1 (IC(50) = 100 +/- 30 microM) and EAAT2 (IC(50) = 300 +/- 80 microM). By contrast, compound 8b showed no appreciable affinity for Glu uptake sites, neither synaptosomal nor cloned. Compounds 9a-c and 10a,b, possessing 1-hydroxyimidazole as the terminal acidic function, were devoid of activity in all of the systems tested. Protolytic properties of compounds 7a,b, 8b, and 9b were determined by titration, and a correlation between the pK(a) values and the activity at AMPA receptors was apparent. Optimized structures of all the synthesized ligands were fitted to the known crystal structure of an AMPA-GluR2 construct. Where substantial reduction or abolition of affinity at AMPA receptors was observed, this could be rationalized on the basis of the ability of the ligand to fit the construct. The results presented in this article point to the utility of 1-hydroxypyrazole and 1,2,3-hydroxytriazole as bioisosteres of carboxylic acids at Glu receptors and transporters. None of the compounds showed significant activity at metabotropic Glu receptors.