Oxadiazoles are five-membered heteroaromatic rings containing two carbons, two nitrogens, and one oxygen atom, and they exist in different regioisomeric forms. Oxadiazoles are frequently occurring motifs in druglike molecules, and they are often used with the intention of being bioisosteric replacements for ester and amide functionalities. The current study presents a systematic comparison of 1,2,4- and 1,3,4-oxadiazole matched pairs in the AstraZeneca compound collection. In virtually all cases, the 1,3,4-oxadiazole isomer shows an order of magnitude lower lipophilicity (log D), as compared to its isomeric partner. Significant differences are also observed with respect to metabolic stability, hERG inhibition, and aqueous solubility, favoring the 1,3,4-oxadiazole isomers. The difference in profile between the 1,2,4 and 1,3,4 regioisomers can be rationalized by their intrinsically different charge distributions (e.g., dipole moments). To facilitate the use of these heteroaromatic rings, novel synthetic routes for ready access of a broad spectrum of 1,3,4-oxadiazoles, under mild conditions, are described.

No abstract

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.

Steroid receptor drugs have been available for more than half a century, but details of the ligand binding mechanism have remained elusive. We solved X-ray structures of the glucocorticoid and mineralocorticoid receptors to identify a conserved plasticity at the helix 6-7 region that extends the ligand binding pocket toward the receptor surface. Since none of the endogenous ligands exploit this region, we hypothesized that it constitutes an integral part of the binding event. Extensive all-atom unbiased ligand exit and entrance simulations corroborate a ligand binding pathway that gives the observed structural plasticity a key functional role. Kinetic measurements reveal that the receptor residence time correlates with structural rearrangements observed in both structures and simulations. Ultimately, our findings reveal why nature has conserved the capacity to open up this region, and highlight how differences in the details of the ligand entry process result in differential evolutionary constraints across the steroid receptors.

Two X-ray structures of the GluR2 ligand-binding core in complex with (S)-2-amino-3-(5-tert-butyl-3-hydroxy-4-isoxazolyl)propionic acid ((S)-ATPA) have been determined with and without Zn(2+) ions. (S)-ATPA induces a domain closure of ca. 21 degrees compared to the apo form. The tert-butyl moiety of (S)-ATPA is buried in a partially hydrophobic pocket and forces the ligand into the glutamate-like binding mode. The structures provide new insight into the molecular basis of agonist selectivity between AMPA and kainate receptors.

A novel series of melanin concentrating hormone receptor 1 (MCHr1) antagonists were the starting point for a drug discovery program that culminated in the discovery of 103 (AZD1979). The lead optimization program was conducted with a focus on reducing lipophilicity and understanding the physicochemical properties governing CNS exposure and undesired off-target pharmacology such as hERG interactions. An integrated approach was taken where the key assay was ex vivo receptor occupancy in mice. The candidate compound 103 displayed appropriate lipophilicity for a CNS indication and showed excellent permeability with no efflux. Preclinical GLP toxicology and safety pharmacology studies were without major findings and 103 was taken into clinical trials.

In this study, we performed an extensive exploration of the ligand entry mechanism for members of the steroid nuclear hormone receptor family (androgen receptor, estrogen receptor α, glucocorticoid receptor, mineralocorticoid receptor, and progesterone receptor) and their endogenous ligands. The exploration revealed a shared entry path through the helix 3, 7, and 11 regions. Examination of the x-ray structures of the receptor-ligand complexes further showed two distinct folds of the helix 6-7 region, classified as "open" and "closed", which could potentially affect ligand binding. To improve sampling of the helix 6-7 loop, we incorporated motion modes based on principal component analysis of existing crystal structures of the receptors and applied them to the protein-ligand sampling. A detailed comparison with the anisotropic network model (an elastic network model) highlights the importance of flexibility in the entrance region. While the binding (interaction) energy of individual simulations can be used to score different ligands, extensive sampling further allows us to predict absolute binding free energies and analyze reaction kinetics using Markov state models and Perron-cluster cluster analysis, respectively. The predicted relative binding free energies for three ligands binding to the progesterone receptor are in very good agreement with experimental results and the Perron-cluster cluster analysis highlighted the importance of a peripheral binding site. Our analysis revealed that the flexibility of the helix 3, 7, and 11 regions represents the most important factor for ligand binding. Furthermore, the hydrophobicity of the ligand influences the transition between the peripheral and the active binding site.

In this study, we present a fully automatic platform based on our Monte Carlo algorithm, the Protein Energy Landscape Exploration method (PELE), for the estimation of absolute protein–ligand binding free energies, one of the most significant challenges in computer aided drug design. Based on a ligand pathway approach, an initial short enhanced sampling simulation is performed to identify reasonable starting positions for more extended sampling. This stepwise approach allows for a significant faster convergence of the free energy estimation using the Markov State Model (MSM) technique. PELE-MSM was applied on four diverse protein and ligand systems, successfully ranking compounds for two systems. Based on the results, current limitations and challenges with physics-based methods in computational structural biology are discussed. Overall, PELE-MSM constitutes a promising step toward computing absolute binding free energies and in their application into drug discovery projects.