An inverse agonist of orphan receptor GPR61 acts by a G protein-competitive allosteric mechanism.

GPR61 is an orphan GPCR related to biogenic amine receptors. Its association with phenotypes relating to appetite makes it of interest as a druggable target to treat disorders of metabolism and body weight, such as obesity and cachexia. To date, the lack of structural information or a known biological ligand or tool compound has hindered comprehensive efforts to study GPR61 structure and function. Here, we report a structural characterization of GPR61, in both its active-like complex with heterotrimeric G protein and in its inactive state. Moreover, we report the discovery of a potent and selective small-molecule inverse agonist against GPR61 and structural elucidation of its allosteric binding site and mode of action. These findings offer mechanistic insights into an orphan GPCR while providing both a structural framework and tool compound to support further studies of GPR61 function and modulation.

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism.

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B' and F-G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.

Critical evaluation of methods to incorporate entropy loss upon binding in high-throughput docking.

Proper accounting of the positional/orientational/conformational entropy loss associated with protein-ligand binding is important to obtain reliable predictions of binding affinity. Herein, we critically examine two simplified statistical mechanics-based approaches, namely a constant penalty per rotor method, and a more rigorous method, referred to here as the partition function-based scoring (PFS) method, to account for such entropy losses in high-throughput docking calculations. Our results on the estrogen receptor beta and dihydrofolate reductase proteins demonstrate that, while the constant penalty method over-penalizes molecules for their conformational flexibility, the PFS method behaves in a more "DeltaG-like" manner by penalizing different rotors differently depending on their residual entropy in the bound state. Furthermore, in contrast to no entropic penalty or the constant penalty approximation, the PFS method does not exhibit any bias towards either rigid or flexible molecules in the hit list. Preliminary enrichment studies using a lead-like random molecular database suggest that an accurate representation of the "true" energy landscape of the protein-ligand complex is critical for reliable predictions of relative binding affinities by the PFS method.

Understanding the selectivity of genistein for human estrogen receptor-beta using X-ray crystallography and computational methods.

We present X-ray crystallographic and molecular modeling studies of estrogen receptors-alpha and -beta complexed with the estrogen receptor-beta-selective phytoestrogen genistein, and coactivator-derived NR box peptides containing an LXXLL motif. We demonstrate that the ligand binding mode is essentially identical when genistein is bound to both isoforms, despite the considerably weaker affinity of this ligand for estrogen receptor-alpha. In addition, we examine subtle differences between binding site residues, providing an explanation for why genistein is modestly selective for the beta isoform. To this end, we also present the results of quantum chemical studies and thermodynamic arguments that yield insight to the nature of the interactions leading to estrogen receptor-beta selectivity. The importance of our analysis to structure-based drug design is discussed.

Molecular determinants of ER alpha and ER beta involved in selectivity of 16 alpha-iodo-7 beta estradiol.

The two known estrogen receptors, ER alpha and ER beta, are hormone inducible transcription factors that have distinct roles in regulating cell proliferation and differentiation. The natural ligand, 17 beta-estradiol (E2), binds with high affinity to both ER alpha and ER beta. However, a close analogue, 16 alpha-iodo-17 beta-estradiol (16 alpha IE2) showed about 10-fold selectivity for ER alpha over ER beta. From X-ray studies, it has been shown that the ligand-binding domains (LBD) of the two receptors are strikingly similar, and that only two changes fall within the binding cavity (ER alpha Leu384 to ER beta Met336, and ER alpha Met421 to ER beta Ile373). To understand the molecular basis for the ER alpha selectivity of 16 alpha IE2, mutants and chimeras of ER alpha and ER beta were generated, and ligand-binding and transactivation functions were studied. The ER alpha Leu384 Met mutant behaved like ER alpha WT in the presence of 16 alpha IE2; whereas the profile of the ER alpha Met421 Ile mutant was similar to that of ER beta WT. The ER beta mutant Ile373 Met behaved like ER alpha with 16 alpha IE2. The results clearly demonstrate the role of ER alpha Met421 in the ER alpha selectivity of 16 alpha IE2.

Identification of phenylsulfone-substituted quinoxaline (WYE-672) as a tissue selective liver X-receptor (LXR) agonist.

A series of phenyl sulfone substituted quinoxaline were prepared and the lead compound 13 (WYE-672) was shown to be a tissue selective LXR Agonist. Compound 13 demonstrated partial agonism for LXRbeta in kidney HEK-293 cells but did not activate Gal4 LXRbeta fusion proteins in huh-7 liver cells. Although 13 showed potent binding affinity to LXRbeta (IC(50) = 53 nM), it had little binding affinity for LXRalpha (IC(50) > 1.0 microM) and did not recruit any coactivator/corepressor peptides in the LXRalpha multiplex assay. However, compound 13 showed good agonism in THP-1 cells with respect to increasing ABCA1 gene expression and good potency on cholesterol efflux in THP-1 foam cells. In an eight-week lesion study in LDLR -/- mice, compound 13 showed reduction of aortic arch lesion progression and no plasma or hepatic triglyceride increase. These results suggest quinoxaline 13 may have an improved biological profile for potential use as a therapeutic agent.

Improvement of physiochemical properties of the tetrahydroazepinoindole series of farnesoid X receptor (FXR) agonists: beneficial modulation of lipids in primates.

In an effort to develop orally active farnesoid X receptor (FXR) agonists, a series of tetrahydroazepinoindoles with appended solubilizing amine functionalities were synthesized. The crystal structure of the previously disclosed FXR agonist, 1 (FXR-450), aided in the design of compounds with tethered solubilizing functionalities designed to reach the solvent cavity around the hFXR receptor. These compounds were soluble in 0.5% methylcellulose/2% Tween-80 in water (MC/T) for oral administration. In vitro and in vivo optimization led to the identification of 14dd and 14cc, which in a dose-dependent fashion regulated low density lipoprotein cholesterol (LDLc) in low density lipoprotein receptor knockout (LDLR(-/-)) mice. Compound 14cc was dosed in female rhesus monkeys for 4 weeks at 60 mg/kg daily in MC/T vehicle. After 7 days, triglyceride (TG) levels and very low density lipoprotein cholesterol (VLDLc) levels were significantly decreased and LDLc was decreased 63%. These data are the first to demonstrate the dramatic lowering of serum LDLc levels by a FXR agonist in primates and supports the potential utility of 14cc in treating dyslipidemia in humans beyond just TG lowering.

Estrogen receptor beta ligands: design and synthesis of new 2-phenyl-isoindole-1,3-diones.

The design, synthesis, and biological evaluation of the 2-phenyl-isoindole-1,3-diones will be discussed. Detailed modeling studies with X-ray support were used to understand the ligand binding orientation and observed selectivity.

ERbeta ligands. Part 6: 6H-Chromeno[4,3-b]quinolines as a new series of estrogen receptor beta-selective ligands.

A new class of estrogen receptor beta (ERbeta) ligands based on the 6H-chromeno[4,3-b]quinoline scaffold has been prepared. Several C7-substituted analogues displayed high affinity and modest selectivity for ERbeta.

Estrogen receptor ligands: design and synthesis of new 2-arylindene-1-ones.

The syntheses of a series of 2-arylindene-1-ones as potent ligands of ERbeta and ERalpha are described. Several compounds exhibited high potency and moderate selectivity for the ERbeta receptor. X-ray and modeling studies were used to understand ligand binding orientation and observed affinity.

Molecular and pharmacological properties of a potent and selective novel nonsteroidal progesterone receptor agonist tanaproget.

Progesterone receptor (PR) agonists have several important applications in women's health, such as in oral contraception and post-menopausal hormone therapy. Currently, all PR agonists used clinically are steroids. Because of their interactions with other steroid receptors, steroid-metabolizing enzymes, or other steroid-signaling pathways, these drugs can pose significant side effects in some women. Efforts to discover novel nonsteroidal PR agonists with improved biological properties led to the discovery of tanaproget (TNPR). TNPR binds to the PR from various species with a higher relative affinity than reference steroidal progestins. In T47D cells, TNPR induces alkaline phosphatase activity with an EC(50) value of 0.1 nm, comparable with potent steroidal progestins such as medroxyprogesterone acetate (MPA) and trimegestone (TMG), albeit with a reduced efficacy ( approximately 60%). In a mammalian two-hybrid assay to measure PR agonist-induced interaction between steroid receptor co-activator-1 and PR, TNPR showed similar potency (EC(50) value of 0.02 nm) and efficacy to MPA and TMG. Importantly, in key animal models such as the rat ovulation inhibition assay, TNPR demonstrates full efficacy and an enhanced progestational potency (30-fold) when compared with MPA and TMG. Furthermore, TNPR has relatively weak interactions with other steroid receptors and binding proteins and little effect on cytochrome P450 metabolic pathways. Finally, the three-dimensional crystal structure of the PR ligand binding domain with TNPR has been delineated to demonstrate how this nonsteroidal ligand achieves its high binding affinity. Therefore, TNPR is a structurally novel and very selective PR agonist with an improved preclinical pharmacological profile.

7-Substituted 2-phenyl-benzofurans as ER beta selective ligands.

A series of 2-(4-hydroxy-phenyl)-benzofuran-5-ols with relatively lipophilic groups in the 7-position of the benzofuran was prepared and the affinity and selectivity for ER beta was measured. Many of the analogues were found to be potent and selective ER beta ligands. Additional modifications at the benzofuran 4-position as well as at the 3'-position of the 2-phenyl group were found to further increase selectivity. Such modifications led to compounds with <10 nM potency and >100-fold selectivity for ER beta.

Structure-based design of estrogen receptor-beta selective ligands.

We present the structure-based optimization of a series of estrogen receptor-beta (ERbeta) selective ligands. X-ray cocrystal structures of these ligands complexed to both ERalpha and ERbeta are described. We also discuss how molecular modeling was used to take advantage of subtle differences between the two binding cavities in order to optimize selectivity for ERbeta over ERalpha. Quantum chemical calculations are utilized to gain insight into the mechanism of selectivity enhancement. Despite only two relatively conservative residue substitutions in the ligand binding pocket, the most selective compounds have greater than 100-fold selectivity for ERbeta relative to ERalpha when measured using a competitive radioligand binding assay.