Structural basis for ligand recognition and activation of the prostanoid receptors

Prostaglandin F2α (PGF2α) and thromboxane A2 (TXA2) are endogenous arachidonic acid metabolites, modulating diverse physiological processes including inflammation and cardiovascular homeostasis through activating PGF2α receptor (FP) and TXA2 receptor (TP). Ligands targeting FP and TP have demonstrated efficacy in treating conditions like glaucoma and cardiovascular diseases in humans, as well as reproductive-related diseases in animals. Here, we present five cryoelectron microscopy structures illustrating FP and TP in complex with Gq and bound to PGF2α (endogenous ligand), latanoprost acid (a clinical drug), and two other synthetic agonists. Combined with mutational and functional studies, these structures reveal not only structural features for the specific recognition of endogenous ligands and attainment of receptor selectivity of FP and TP but also the common mechanisms of receptor activation and Gq protein coupling. The findings may enrich our knowledge of ligand recognition and signal transduction of the prostanoid receptor family and facilitate rational ligand design toward these two receptors.

Probing ligand selectivity in pathogens

Why does protein kinase A respond to purine nucleosides in certain pathogens, but not to the cyclic nucleotides that activate this kinase in most other organisms?

Molecular recognition of niacin and lipid-lowering drugs by the human hydroxycarboxylic acid receptor 2

Niacin, an age-old lipid-lowering drug, acts through the hydroxycarboxylic acid receptor 2 (HCAR2), a G-protein-coupled receptor (GPCR). Yet, its use is hindered by side effects like skin flushing. To address this, specific HCAR2 agonists, like MK-6892 and GSK256073, with fewer adverse effects have been created. However, the activation mechanism of HCAR2 by niacin and these new agonists is not well understood. Here, we present three cryoelectron microscopy structures of Gi-coupled HCAR2 bound to niacin, MK-6892, and GSK256073. Our findings show that different ligands induce varying binding pockets in HCAR2, influenced by aromatic amino acid clusters (W91ECL1, H1614.59, W1885.38, H1895.39, and F1935.43) from receptors ECL1, TM4, and TM5. Additionally, conserved residues R1113.36 and Y2847.43, unique to the HCA receptor family, likely initiate activation signal propagation in HCAR2. This study provides insights into ligand recognition, receptor activation, and G protein coupling mediated by HCAR2, laying the groundwork for developing HCAR2-targeted drugs.

Structural and biochemical insights into PsEst3, a new GHSR-type esterase obtained from Paenibacillus sp. R4

PsEst3, a psychrophilic esterase obtained from Paenibacillus sp. R4, which was isolated from the permafrost of Alaska, exhibits relatively high activity at low temperatures. Here, crystal structures of PsEst3 complexed with various ligands were generated and studied at atomic resolution, and biochemical studies were performed to analyze the structure-function relationship of PsEst3. Certain unique characteristics of PsEst3 distinct from those of other classes of lipases/esterases were identified. Firstly, PsEst3 contains a conserved GHSRA/G pentapeptide sequence in the GxSxG motif around the nucleophilic serine. Additionally, it contains a conserved HGFR/K consensus sequence in the oxyanion hole, which is distinct from that in other lipase/esterase families, as well as a specific domain composition (for example a helix-turn-helix motif) and a degenerative lid domain that exposes the active site to the solvent. Secondly, the electrostatic potential of the active site in PsEst3 is positive, which may cause unintended binding of negatively charged chemicals in the active site. Thirdly, the last residue of the oxyanion hole-forming sequence, Arg44, separates the active site from the solvent by sealing the acyl-binding pocket, suggesting that PsEst3 is an enzyme that is customized to sense an unidentified substrate that is distinct from those of classical lipases/esterases. Collectively, this evidence strongly suggests that PsEst3 belongs to a distinct family of esterases.

Structures of the entire human opioid receptor family

Opioids are effective analgesics, but their use is beset by serious side effects, including addiction and respiratory depression, which contribute to the ongoing opioid crisis. The human opioid system contains four opioid receptors (μOR, δOR, κOR, and NOPR) and a set of related endogenous opioid peptides (EOPs), which show distinct selectivity toward their respective opioid receptors (ORs). Despite being key to the development of safer analgesics, the mechanisms of molecular recognition and selectivity of EOPs to ORs remain unclear. Here, we systematically characterize the binding of EOPs to ORs and present five structures of EOP-OR-Gi complexes, including β-endorphin- and endomorphin-bound μOR, deltorphin-bound δOR, dynorphin-bound κOR, and nociceptin-bound NOPR. These structures, supported by biochemical results, uncover the specific recognition and selectivity of opioid peptides and the conserved mechanism of opioid receptor activation. These results provide a structural framework to facilitate rational design of safer opioid drugs for pain relief.

The Third Extracellular Loop of Mammalian Odorant Receptors Is Involved in Ligand Binding

Mammals recognize chemicals in the air via G protein-coupled odorant receptors (ORs). In addition to their orthosteric binding site, other segments of these receptors modulate ligand recognition. Focusing on human hOR1A1, which is considered prototypical of class II ORs, we used a combination of molecular modeling, site-directed mutagenesis, and in vitro functional assays. We showed that the third extracellular loop of ORs (ECL3) contributes to ligand recognition and receptor activation. Indeed, site-directed mutations in ECL3 showed differential effects on the potency and efficacy of both carvones, citronellol, and 2-nonanone.

Soluble guanylyl cyclase: Molecular basis for ligand selectivity and action in vitro and in vivo

Nitric oxide (NO), carbon monoxide (CO), oxygen (O₂), hydrogen sulfide (H₂S) are gaseous molecules that play important roles in the physiology and pathophysiology of eukaryotes. Tissue concentrations of these physiologically relevant gases vary remarkable from nM range for NO to high μM range of O₂. Various hemoproteins play a significant role in sensing and transducing cellular signals encoded by gaseous molecules or in transporting them. Soluble guanylyl cyclase (sGC) is a hemoprotein that plays vital roles in a wide range of physiological functions and combines the functions of gaseous sensor and signal transducer. sGC uniquely evolved to sense low non-toxic levels of NO and respond to elevated NO levels by increasing its catalytic ability to generate the secondary signaling messenger cyclic guanosine monophosphate (cGMP). This review discusses sGC's gaseous ligand selectivity and the molecular basis for sGC function as high-affinity and selectivity NO receptor. The effects of other gaseous molecules and small molecules of cellular origin on sGC's function are also discussed.

Structural perspectives on the mechanism of signal activation, ligand selectivity and allosteric modulation in angiotensin receptors: IUPHAR Review 34

Functional advances have guided our knowledge of physiological and fatal pathological mechanisms of the hormone angiotensin II (AngII) and its antagonists. Such studies revealed that tissue response to a given dose of the hormone or its antagonist depends on receptors that engage the ligand. Thus, we need to know much more about the structures of receptor-ligand complexes at high resolution. Recently, X-ray structures of both AngII receptors (AT1 and AT2 receptors) bound to peptide and non-peptide ligands have been elucidated, providing new opportunities to examine the dynamic fluxes in the 3D architecture of the receptors, as the basis of ligand selectivity, efficacy, and regulation of the molecular functions of the receptors. Constituent structural motifs cooperatively transform ligand selectivity into specific functions, thus conceptualizing the primacy of the 3D structure over individual motifs of receptors. This review covers the new data elucidating the structural dynamics of AngII receptors and how structural knowledge can be transformative in understanding the mechanisms underlying the physiology of AngII.

GPCRs steer Gi and Gs selectivity via TM5-TM6 switches as revealed by structures of serotonin receptors

Serotonin (or 5-hydroxytryptamine, 5-HT) is an important neurotransmitter that activates 12 different G protein-coupled receptors (GPCRs) through selective coupling of G_s, G_i, or G_q proteins. The structural basis for G protein subtype selectivity by these GPCRs remains elusive. Here, we report the structures of the serotonin receptors 5-HT₄, 5-HT₆, and 5-HT₇ with G_s, and 5-HT₄ with G_i1. The structures reveal that transmembrane helices TM5 and TM6 alternate lengths as a macro-switch to determine receptor's selectivity for G_s and G_i, respectively. We find that the macro-switch by the TM5-TM6 length is shared by class A GPCR-G protein structures. Furthermore, we discover specific residues within TM5 and TM6 that function as micro-switches to form specific interactions with G_s or G_i. Together, these results present a common mechanism of G_s versus G_i protein coupling selectivity or promiscuity by class A GPCRs and extend the basis of ligand recognition at serotonin receptors.

Cryo-EM structure of the β3-adrenergic receptor reveals the molecular basis of subtype selectivity

The β₃-adrenergic receptor (β₃AR) is predominantly expressed in adipose tissue and urinary bladder and has emerged as an attractive drug target for the treatment of type 2 diabetes, obesity, and overactive bladder (OAB). Here, we report the cryogenic electron microscopy structure of the β₃AR-G_s signaling complex with the selective agonist mirabegron, a first-in-class drug for OAB. Comparison of this structure with the previously reported β₁AR and β₂AR structures reveals a receptor activation mechanism upon mirabegron binding to the orthosteric site. Notably, the narrower exosite in β₃AR creates a perpendicular pocket for mirabegron. Mutational analyses suggest that a combination of both the exosite shape and the amino-acid-residue substitutions defines the drug selectivity of the βAR agonists. Our findings provide a molecular basis for βAR subtype selectivity, allowing the design of more-selective agents with fewer adverse effects.

Selectivity in the Ligand Functionalization of Photocatalytic Metal Oxide Nanoparticles for Phase Transfer and Self-Assembly Applications

The functionalization of photocatalytic metal oxide nanoparticles of TiO₂ , ZnO, WO₃ and CuO with amine-terminated (oleylamine) and thiol-terminated (dodecane-1-thiol) alkyl-chain ligands was studied under ambient conditions. A high selectivity was observed in the binding specificity of a ligand towards nanoparticles of these different oxides. It was observed that oleylamine binds stably to only TiO₂ and WO₃ , whereas dodecane-1-thiol binds stably only to ZnO and CuO. Similarly, polar-to-nonpolar solvent phase transfer of TiO₂ and WO₃ nanoparticles could be achieved by using oleylamine, but not dodecane-1-thiol, whereas the opposite holds for ZnO and CuO. The surface chemistry of ligand-functionalized nanoparticles was probed by attenuated total reflectance (ATR)-FTIR spectroscopy, which enabled the occupation of the ligands at the active sites to be elucidated. The photostability of the ligands on the nanoparticle surface was determined by the photocatalytic self-cleaning properties of the material. Although TiO₂ and WO₃ degrade the ligands within 24 h under both UV and visible light, ligands on ZnO and CuO remain unaffected. The gathered insights are also highly relevant from an application point of view. As an example, because the ligand-functionalized nanoparticles are hydrophobic in nature, they can be self-assembled at the air-water interface to give nanoparticle films with demonstrated photocatalytic as well as anti-fogging properties.

Molecular Basis of Class B GPCR Selectivity for the Neuropeptides PACAP and VIP

The related neuropeptides PACAP and VIP, and their shared PAC1, VPAC1 and VPAC2 receptors, regulate a large array of physiological activities in the central and peripheral nervous systems. However, the lack of comparative and molecular mechanistic investigations hinder further understanding of their preferred binding selectivity and function. PACAP and VIP have comparable affinity at the VPAC1 and VPAC2 receptor, but PACAP is 400-1,000 fold more potent than VIP at the PAC1 receptor. A molecular understanding of the differing neuropeptide-receptor interactions and the details underlying the receptor transitions leading to receptor activation are much needed for the rational design of selective ligands. To these ends, we have combined structural information and advanced simulation techniques to study PACAP/VIP binding selectivity, full-length receptor conformation ensembles and transitions of the PACAP/VIP receptor variants and subtypes, and a few key interactions in the orthosteric-binding pocket. Our results reveal differential peptide-receptor interactions (at the atomistic detail) important for PAC1, VPAC1 and VPAC2 receptor ligand selectivity. Using microsecond-long molecular dynamics simulations and the Markov State Models, we have also identified diverse receptor conformational ensembles and microstate transition paths for each receptor, the potential mechanisms underlying receptor open and closed states, and the interactions and dynamics at the transmembrane orthosteric pocket for receptor activation. These analyses reveal important features in class B GPCR structure-dynamics-function relationships, which provide novel insights for structure-based drug discovery.

Structural insights into the human D1 and D2 dopamine receptor signaling complexes

The D1- and D2-dopamine receptors (D1R and D2R), which signal through Gs and Gi, respectively, represent the principal stimulatory and inhibitory dopamine receptors in the central nervous system. D1R and D2R also represent the main therapeutic targets for Parkinson's disease, schizophrenia, and many other neuropsychiatric disorders, and insight into their signaling is essential for understanding both therapeutic and side effects of dopaminergic drugs. Here, we report four cryoelectron microscopy (cryo-EM) structures of D1R-Gs and D2R-Gi signaling complexes with selective and non-selective dopamine agonists, including two currently used anti-Parkinson's disease drugs, apomorphine and bromocriptine. These structures, together with mutagenesis studies, reveal the conserved binding mode of dopamine agonists, the unique pocket topology underlying ligand selectivity, the conformational changes in receptor activation, and potential structural determinants for G protein-coupling selectivity. These results provide both a molecular understanding of dopamine signaling and multiple structural templates for drug design targeting the dopaminergic system.

Structures of the human dopamine D3 receptor-Gi complexes

The dopamine system, including five dopamine receptors (D1R-D5R), plays essential roles in the central nervous system (CNS), and ligands that activate dopamine receptors have been used to treat many neuropsychiatric disorders. Here, we report two cryo-EM structures of human D3R in complex with an inhibitory G protein and bound to the D3R-selective agonists PD128907 and pramipexole, the latter of which is used to treat patients with Parkinson's disease. The structures reveal agonist binding modes distinct from the antagonist-bound D3R structure and conformational signatures for ligand-induced receptor activation. Mutagenesis and homology modeling illuminate determinants of ligand specificity across dopamine receptors and the mechanisms for G_i protein coupling. Collectively our work reveals the basis of agonist binding and ligand-induced receptor activation and provides structural templates for designing specific ligands to treat CNS diseases targeting the dopaminergic system.

Aryl Hydrocarbon Receptor (AHR) Ligands as Selective AHR Modulators (SAhRMs)

The aryl hydrocarbon receptor (AhR) was first identified as the intracellular protein that bound and mediated the toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD, dioxin) and dioxin-like compounds (DLCs). Subsequent studies show that the AhR plays an important role in maintaining cellular homeostasis and in pathophysiology, and there is increasing evidence that the AhR is an important drug target. The AhR binds structurally diverse compounds, including pharmaceuticals, phytochemicals and endogenous biochemicals, some of which may serve as endogenous ligands. Classification of DLCs and non-DLCs based on their persistence (metabolism), toxicities, binding to wild-type/mutant AhR and structural similarities have been reported. This review provides data suggesting that ligands for the AhR are selective AhR modulators (SAhRMs) that exhibit tissue/cell-specific AhR agonist and antagonist activities, and that their functional diversity is similar to selective receptor modulators that target steroid hormone and other nuclear receptors.

The Molecular Basis for Purine Binding Selectivity in the Bacterial ATP Synthase ϵ Subunit

The ϵ subunit of ATP synthases has been proposed to regulate ATP hydrolysis in bacteria. Prevailing evidence supports the notion that when the ATP concentration falls below a certain threshold, the ϵ subunit changes its conformation from a non-inhibitory down-state to an extended up-state that then inhibits enzymatic ATP hydrolysis by binding to the catalytic domain. It has been demonstrated that the ϵ subunit from Bacillus PS3 is selective for ATP over other nucleotides, including GTP. In this study, the purine triphosphate selectivity is rationalized by using results from MD simulations and free energy calculations for the R103A/R115A mutant of the ϵ subunit from Bacillus PS3, which binds ATP more strongly than the wild-type protein. Our results are in good agreement with experimental data, and the elucidated molecular basis for selectivity could help to guide the design of novel GTP sensors.

Custom G4 Microarrays Reveal Selective G-Quadruplex Recognition of Small Molecule BMVC: A Large-Scale Assessment of Ligand Binding Selectivity

G-quadruplexes (G4) are considered new drug targets for human diseases such as cancer. More than 10,000 G4s have been discovered in human chromatin, posing challenges for assessing the selectivity of a G4-interactive ligand. 3,6-bis(1-Methyl-4-vinylpyridinium) carbazole diiodide (BMVC) is the first fluorescent small molecule for G4 detection in vivo. Our previous structural study shows that BMVC binds to the MYC promoter G4 (MycG4) with high specificity. Here, we utilize high-throughput, large-scale custom DNA G4 microarrays to analyze the G4-binding selectivity of BMVC. BMVC preferentially binds to the parallel MycG4 and selectively recognizes flanking sequences of parallel G4s, especially the 3′-flanking thymine. Importantly, the microarray results are confirmed by orthogonal NMR and fluorescence binding analyses. Our study demonstrates the potential of custom G4 microarrays as a platform to broadly and unbiasedly assess the binding selectivity of G4-interactive ligands, and to help understand the properties that govern molecular recognition.

Comprehensive Exploration of Target-specific Ligands Using a Graph Convolution Neural Network

Machine learning approaches are widely used to evaluate ligand activities of chemical compounds toward potential target proteins. Especially, exploration of highly selective ligands is important for the development of new drugs with higher safety. One difficulty in constructing well-performing model predicting such a ligand activity is the absence of data on true negative ligand-protein interactions. In other words, in many cases we can access to plenty of information on ligands that bind to specific protein, but less or almost no information showing that compounds don't bind to proteins of interest. In this paper, we suggested an approach to comprehensively explore candidates for ligands specifically targeting toward proteins without using information on the true negative interaction. The approach consists of 4 steps: 1) constructing a model that distinguishes ligands for the target proteins of interest from those targeting proteins that cause off-target effects, by using graph convolution neural network (GCNN); 2) extracting feature vectors after convolution/pooling processes and mapping their principal components in two dimensions; 3) specifying regions with higher density for two ligand groups through kernel density estimation; and 4) investigating the distribution of compounds for exploration on the density map using the same classifier and decomposer. If compounds for exploration are located in higher-density regions of ligand compounds, these compounds can be regarded as having relatively high binding affinity to the major target or off-target proteins compared with other compounds. We applied the approach to the exploration of ligands for β-site amyloid precursor protein [APP]-cleaving enzyme 1 (BACE1), a major target for Alzheimer Disease (AD), with less off-target effect toward cathepsin D. We demonstrated that the density region of BACE1 and cathepsin D ligands are well-divided, and a group of natural compounds as a target for exploration of new drug candidates also has significantly different distribution on the density map.

N 6-modified cAMP derivatives that activate protein kinase A also act as full agonists of murine HCN2 channels

cAMP acts as a second messenger in many cellular processes. Three protein types mainly mediate cAMP-induced effects: PKA, exchange protein directly activated by cAMP (Epac), and cyclic nucleotide-modulated channels (cyclic nucleotide-gated or hyperpolarization-activated and cyclic nucleotide-modulated (HCN) channels). Discrimination among these cAMP signaling pathways requires specific targeting of only one protein. Previously, cAMP modifications at position N ⁶ of the adenine ring (PKA) and position 2'-OH of the ribose (Epac) have been used to produce target-selective compounds. However, cyclic nucleotide-modulated ion channels were usually outside of the scope of these previous studies. These channels are widely distributed, so possible channel cross-activation by PKA- or Epac-selective agonists warrants serious consideration. Here we demonstrate the agonistic effects of three PKA-selective cAMP derivatives, N ⁶-phenyladenosine-3',5'-cyclic monophosphate (N ⁶-Phe-cAMP), N ⁶-benzyladenosine-3',5'-cyclic monophosphate (N ⁶-Bn-cAMP), and N ⁶-benzoyl-adenosine-3',5'-cyclic monophosphate (N ⁶-Bnz-cAMP), on murine HCN2 pacemaker channels. Electrophysiological characterization in Xenopus oocytes revealed that these derivatives differ in apparent affinities depending on the modification type but that their efficacy and effects on HCN2 activation kinetics are similar to those of cAMP. Docking experiments suggested a pivotal role of Arg-635 at the entrance of the binding pocket in HCN2, either causing stabilizing cation-π interactions with the aromatic ring in N ⁶-Phe-cAMP or N ⁶-Bn-cAMP or a steric clash with the aromatic ring in N ⁶-Bnz-cAMP. A reduced apparent affinity of N ⁶-Phe-cAMP toward the variants R635A and R635E strengthened that notion. We conclude that some PKA activators also effectively activate HCN2 channels. Hence, when studying PKA-mediated cAMP signaling with cAMP derivatives in a native environment, activation of HCN channels should be considered.

Molecular Determinants for Ligand Selectivity of the Cell-Free Synthesized Human Endothelin B Receptor

Extracellular domains of G-protein-coupled receptors act as initial molecular selectivity filters for subtype specific ligands and drugs. Chimeras of the human endothelin-B receptor containing structural units from the extracellular domains of the endothelin-A receptor were analyzed after their co-translational insertion into preformed nanodiscs. A short β-strand and a linker region in the second extracellular loop as well as parts of the extracellular N-terminal domain were identified as molecular discrimination sites for the endothelin-B receptor-selective agonists IRL1620, sarafotoxin 6c, 4Ala-ET-1 and ET-3, but not for the non-selective agonist ET-1 recognized by both endothelin receptors. A proposed second disulfide bridge in the endothelin-B receptor tethering the N-terminal domain with the third extracellular loop was not essential for ET-1 recognition and binding, but increased the receptor thermostability. We further demonstrate an experimental approach with cell-free synthesized engineered agonists to analyze the differential discrimination of peptide ligand topologies by the two endothelin receptors. The study is based on the engineering and cell-free insertion of G-protein-coupled receptors into defined membranes and may become interesting also for other targets as an alternative platform to reveal molecular details of ligand selectivity and ligand binding mechanisms.

Binding, Thermodynamics, and Selectivity of a Non-peptide Antagonist to the Melanocortin-4 Receptor

The melanocortin-4 receptor (MC4R) is a potential drug target for treatment of obesity, anxiety, depression, and sexual dysfunction. Crystal structures for MC4R are not yet available, which has hindered successful structure-based drug design. Using microsecond-scale molecular-dynamics simulations, we have investigated selective binding of the non-peptide antagonist MCL0129 to a homology model of human MC4R (hMC4R). This approach revealed that, at the end of a multi-step binding process, MCL0129 spontaneously adopts a binding mode in which it blocks the agonistic-binding site. This binding mode was confirmed in subsequent metadynamics simulations, which gave an affinity for human hMC4R that matches the experimentally determined value. Extending our simulations of MCL0129 binding to hMC1R and hMC3R, we find that receptor subtype selectivity for hMC4R depends on few amino acids located in various structural elements of the receptor. These insights may support rational drug design targeting the melanocortin systems.

The small-molecule inhibitor selectivity between IKKα and IKKβ kinases in NF-κB signaling pathway

The enzyme complex IκB kinase (IKK) is an essential activator of NF-κB signaling pathway involved in propagating the cellular response to inflammation. The complex contains two functional subunits IKKα and IKKβ, which are structurally conserved kinases and selective inhibition of them would result in distinct biological effects. However, most existing IKK inhibitors show moderate or high promiscuity for the two homologous kinases. Understanding of the molecular mechanism and biological implication underlying the specific interactions in IKK-ligand recognition is thus fundamentally important for the rational design of selective IKK inhibitors. In the current work, we integrated molecular docking, quantum mechanics/molecular mechanics calculation and Poisson-Boltzmann/surface area analysis to investigate the structural basis and energetic property of the selective binding of small-molecule ligands to IKKα and IKKβ. It was found that the selectivity is primarily determined by the size and topology difference in ATP-binding pocket of IKKα and IKKβ kinase domains; bulky inhibitor molecules commonly have, respectively, low and appropriate affinities towards IKKα and IKKβ, and thus exhibit relatively high selectivity for IKKβ over IKKα, whereas small ligands can only bind weakly to both the two kinases with low selectivity. In addition, the conformation, arrangement and distribution of residues in IKK pockets are also responsible for constituting the exquisite specificity of ligand binding to KKα and IKKβ. Next, a novel quantitative structure-selectivity relationship model was developed to characterize the relative contribution of each kinase residue to inhibitor selectivity and to predict the selectivity and specificity for a number of known IKK inhibitors. Results showed that the active-site residues contribute significantly to the selectivity by directly interacting with inhibitor ligands, while those protein portions far away from the kinase active sites may also play an important role in determining the selectivity through long-range non-bonded forces and indirect allosteric effect.