A structural chemogenomics analysis of aminergic GPCRs: lessons for histamine receptor ligand design

Computational insights into diverse binding modes of the allosteric modulator and their regulation on dopamine D1 receptor

Allosteric drugs hold the promise of addressing many challenges in the current drug development of GPCRs. However, the molecular mechanism underlying their allosteric modulations remain largely elusive. The dopamine D1 receptor (DRD1), a member of Class A GPCRs, is critical for treating psychiatric disorders, and LY3154207 serves as its promising positive allosteric modulator (PAM). In the work, we utilized extensive Gaussian-accelerated molecular dynamics simulations (a total of 41μs) for the first time probe the diverse binding modes of the allosteric modulator and their regulation effects, based on the DRD1 and LY3154207 as representative. Our simulations identify four binding modes of LY3154207 (one boat mode, two metastable vertical modes and a novel cleft-anchored mode), in which the boat mode is the most stable while there three modes are similar in the stability. However, it is interesting to observed that the most stable boat mode inversely exhibits the weakest positive allosteric effect on influencing the orthosteric ligand binding and maintaining the activity of the transducer binding site. It should result from its induced weaker correlation between the allosteric site and the orthosteric site, and between the orthosteric site and the transducer binding site than the other three binding modes, as well as its weakened interaction between a crucial activation-related residue (S2025.46) and the orthosteric ligand (dopamine). Overall, the work offers atomic-level information to advance our understanding of the complex allosteric regulation on GPCRs, which is beneficial to the allosteric modulator design and development.

Pharmacological Characterization of the Zebrafish (Danio Rerio) Histamine H1 Receptor Reveals the Involvement of the Second Extracellular Loop in the Binding of Histamine

The zebrafish (Danio rerio) histamine H1 receptor gene (zfH1R) was cloned in 2007 and reported to be involved in fish locomotion. Yet, no detailed characterization of its pharmacology and signaling properties have so far been reported. In this study, we pharmacologically characterized the zfH1R expressed in HEK-293T cells by means of [3H]-mepyramine binding and G protein-signaling assays. The zfH1R [dissociation constant (KD), 0.7 nM] displayed similar affinity for the antagonist [3H]-mepyramine as the human histamine H1 receptor (hH1R) (KD, 1.5 nM), whereas the affinity for histamine is 100-fold higher than for the human H1R. The zfH1R couples to Gαq/11 proteins and activates several reporter genes, i.e., NFAT, NFϰB, CRE, VEGF, COX-2, SRE, and AP-1, and zfH1R-mediated signaling is prevented by the Gαq/11 inhibitor YM-254890 and the antagonist mepyramine. Molecular modeling of the zfH1R and human H1R shows that the binding pockets are identical, implying that variations along the ligand binding pathway could underly the differences in histamine affinity instead. Targeting differentially charged residues in extracellular loop 2 (ECL2) using site-directed mutagenesis revealed that Arg21045x55 is most likely involved in the binding process of histamine in zfH1R. This study aids the understanding of the pharmacological differences between H1R orthologs and the role of ECL2 in histamine binding and provides fundamental information for the understanding of the histaminergic system in the zebrafish. SIGNIFICANCE STATEMENT: The use of the zebrafish as in vivo models in neuroscience is growing exponentially, which asks for detailed characterization of the aminergic neurotransmitter systems in this model. This study is the first to pharmacologically characterize the zebrafish histamine H1 receptor after expression in HEK-293T cells. The results show a high pharmacological and functional resemblance with the human ortholog but also reveal interesting structural differences and unveils an important role of the second extracellular loop in histamine binding.

Structural insights into the agonists binding and receptor selectivity of human histamine H4 receptor

Histamine is a biogenic amine that participates in allergic and inflammatory processes by stimulating histamine receptors. The histamine H4 receptor (H4R) is a potential therapeutic target for chronic inflammatory diseases such as asthma and atopic dermatitis. Here, we show the cryo-electron microscopy structures of the H4R-Gq complex bound with an endogenous agonist histamine or the selective agonist imetit bound in the orthosteric binding pocket. The structures demonstrate binding mode of histamine agonists and that the subtype-selective agonist binding causes conformational changes in Phe3447.39, which, in turn, form the "aromatic slot". The results provide insights into the molecular underpinnings of the agonism of H4R and subtype selectivity of histamine receptors, and show that the H4R structures may be valuable in rational drug design of drugs targeting the H4R.

Xanomeline displays concomitant orthosteric and allosteric binding modes at the M4 mAChR

The M4 muscarinic acetylcholine receptor (M4 mAChR) has emerged as a drug target of high therapeutic interest due to its expression in regions of the brain involved in the regulation of psychosis, cognition, and addiction. The mAChR agonist, xanomeline, has provided significant improvement in the Positive and Negative Symptom Scale (PANSS) scores in a Phase II clinical trial for the treatment of patients suffering from schizophrenia. Here we report the active state cryo-EM structure of xanomeline bound to the human M4 mAChR in complex with the heterotrimeric Gi1 transducer protein. Unexpectedly, two molecules of xanomeline were found to concomitantly bind to the monomeric M4 mAChR, with one molecule bound in the orthosteric (acetylcholine-binding) site and a second molecule in an extracellular vestibular allosteric site. Molecular dynamic simulations supports the structural findings, and pharmacological validation confirmed that xanomeline acts as a dual orthosteric and allosteric ligand at the human M4 mAChR. These findings provide a basis for further understanding xanomeline's complex pharmacology and highlight the myriad of ways through which clinically relevant ligands can bind to and regulate GPCRs.

Role of adrenergic receptors in shock

Shock is a severe, life-threatening medical condition with a high mortality rate worldwide. All four major categories of shock (along with their various subtypes)-hypovolemic, distributive, cardiogenic, and obstructive, involve a dramatic mismatch between oxygen supply and demand, and share standard features of decreased cardiac output, reduced blood pressure, and overall hypoperfusion. Immediate and appropriate intervention is required regardless of shock type, as a delay can result in cellular dysfunction, irreversible multiple organ failure, and death. Studies have shown that dysfunction and downregulation of adrenergic receptors (ARs) are often implicated in these shock conditions; for example, their density is shown to be decreased in hypovolemic and cardiogenic shock, while their reduced signaling in the brain and vasculature decrease blood perfusion and oxygen supply. There are two main categories of ARs, α, and β, each with its subtypes and distributions. Our group has demonstrated that a dose of .02 mg/kg body wt of centhaquine (CQ) specifically activates α2B ARs on venous circulation along with the central α2A ARs after hypovolemic/hemorrhagic shock. Activating these receptors by CQ increases cardiac output (CO) and reduces systemic vascular resistance (SVR), with a net increase in blood pressure and tissue perfusion. The clinical trials of CQ conducted by Pharmazz Inc. in India have demonstrated significantly improved survival in shock patients. CQ improved blood pressure and shock index, indicating better blood circulation, and reduced lactate levels in the blood compared to in-use standard resuscitative agents. After successful clinical trials, CQ is being marketed as a drug (Lyfaquin^®) for hypovolemic/hemorrhagic shock in India, and United States FDA has approved the phase III IND application. It is anticipated that the phase III trial in the United States will begin in 2023. Thus, we have demonstrated that α2 ARs could be suitable targets for treating or managing hypovolemic/hemorrhagic shock. Further understanding of ARs in shock would help find new potential pharmacological targets.

Mind the Gap—Deciphering GPCR Pharmacology Using 3D Pharmacophores and Artificial Intelligence

G protein-coupled receptors (GPCRs) are amongst the most pharmaceutically relevant and well-studied protein targets, yet unanswered questions in the field leave significant gaps in our understanding of their nuanced structure and function. Three-dimensional pharmacophore models are powerful computational tools in in silico drug discovery, presenting myriad opportunities for the integration of GPCR structural biology and cheminformatics. This review highlights success stories in the application of 3D pharmacophore modeling to de novo drug design, the discovery of biased and allosteric ligands, scaffold hopping, QSAR analysis, hit-to-lead optimization, GPCR de-orphanization, mechanistic understanding of GPCR pharmacology and the elucidation of ligand–receptor interactions. Furthermore, advances in the incorporation of dynamics and machine learning are highlighted. The review will analyze challenges in the field of GPCR drug discovery, detailing how 3D pharmacophore modeling can be used to address them. Finally, we will present opportunities afforded by 3D pharmacophore modeling in the advancement of our understanding and targeting of GPCRs.

Binding and Activation of Serotonergic G-Protein Coupled Receptors by the Multimodal Antidepressant Vortioxetine

G-protein coupled receptors (GPCRs) are important pharmacological targets. Despite substantial progress, important questions still remain concerning the details of activation: how can a ligand act as an agonist in one receptor but as an antagonist in a homologous receptor, and how can agonists activate a receptor despite lacking polar functional groups able to interact with helix 5 as is the case for the related adrenergic receptors? Studying vortioxetine (VXT), an important multimodal antidepressant drug, may elucidate both questions. Herein, we present a thorough in silico analysis of VXT binding to 5-HT_1A, 5-HT_1B, and 5-HT₇ receptors and compare it with available experimental data. We are able to rationalize the differential mode of action of VXT at different receptors, but also, in the case of the 5-HT_1A receptor, we observe the initial steps of activation that inform about an activation mechanism that does not involve polar interaction with helix 5. The results extend our current understanding of agonist and antagonist action at aminergic GPCRs.

Crystal structure of the α1B-adrenergic receptor reveals molecular determinants of selective ligand recognition

α-adrenergic receptors (αARs) are G protein-coupled receptors that regulate vital functions of the cardiovascular and nervous systems. The therapeutic potential of αARs, however, is largely unexploited and hampered by the scarcity of subtype-selective ligands. Moreover, several aminergic drugs either show off-target binding to αARs or fail to interact with the desired subtype. Here, we report the crystal structure of human α_1BAR bound to the inverse agonist (+)-cyclazosin, enabled by the fusion to a DARPin crystallization chaperone. The α_1BAR structure allows the identification of two unique secondary binding pockets. By structural comparison of α_1BAR with α₂ARs, and by constructing α_1BAR-α_2CAR chimeras, we identify residues 3.29 and 6.55 as key determinants of ligand selectivity. Our findings provide a basis for discovery of α_1BAR-selective ligands and may guide the optimization of aminergic drugs to prevent off-target binding to αARs, or to elicit a selective interaction with the desired subtype.

This study reports the X-ray structure of the α_1B-adrenergic G protein-coupled receptor bound to an inverse agonist, and unveils key determinants of subtype-selective ligand binding that may help the design of aminergic drugs with fewer side-effects.

Modulating physicochemical properties of tetrahydropyridine-2-amine BACE1 inhibitors with electron-withdrawing groups: A systematic study

A common challenge for medicinal chemists is to reduce the pK_a of strongly basic groups' conjugate acids into a range that preserves the desired effects, usually potency and/or solubility, but avoids undesired effects like high volume of distribution (V_d), limited membrane permeation, and off-target binding to, notably, the hERG channel and monoamine receptors. We faced this challenge with a 3,4,5,6-tetrahydropyridine-2-amine scaffold harboring an amidine, a key structural component of potential inhibitors of BACE1, the rate-limiting enzyme in the production of Aβ species that make up amyloid plaques in Alzheimer's disease. In our endeavor to balance potency with desirable properties to achieve brain penetration, we introduced a diverse set of groups in beta position of the amidine that modulate logD, PSA and pK_a. Given the synthetic challenge to prepare these highly functionalized warheads, we first developed a design flow including predicted physicochemical parameters which allowed us to select only the most promising candidates for synthesis. For this we evaluated a set of commercial packages to predict physicochemical properties, which can guide medicinal chemists in their endeavors to modulate pK_a values of amidine and amine bases.

New Chemical Biology Tools for the Histamine Receptor Family

The histamine research community has in the last decade been very active and generated a number of exciting new chemical biology tools for the study of histamine receptors, their ligands, and their pharmacology. In this paper we describe the development of histamine receptor structural biology, the use of receptor conformational biosensors, and the development of new ligands for covalent or fluorescent labeling or for photopharmacological approaches (photocaging and photoswitching). These new tools allow new approaches to study histamine receptors and hopefully will lead to better insights in the molecular aspects of histamine receptors and their ligands.

The Structural Basis of Peptide Binding at Class A G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) represent the largest membrane protein family and a significant target class for therapeutics. Receptors from GPCRs’ largest class, class A, influence virtually every aspect of human physiology. About 45% of the members of this family endogenously bind flexible peptides or peptides segments within larger protein ligands. While many of these peptides have been structurally characterized in their solution state, the few studies of peptides in their receptor-bound state suggest that these peptides interact with a shared set of residues and undergo significant conformational changes. For the purpose of understanding binding dynamics and the development of peptidomimetic drug compounds, further studies should investigate the peptide ligands that are complexed to their cognate receptor.

PyPLIF HIPPOS-Assisted Prediction of Molecular Determinants of Ligand Binding to Receptors

Identification of molecular determinants of receptor-ligand binding could significantly increase the quality of structure-based virtual screening protocols. In turn, drug design process, especially the fragment-based approaches, could benefit from the knowledge. Retrospective virtual screening campaigns by employing AutoDock Vina followed by protein-ligand interaction fingerprinting (PLIF) identification by using recently published PyPLIF HIPPOS were the main techniques used here. The ligands and decoys datasets from the enhanced version of the database of useful decoys (DUDE) targeting human G protein-coupled receptors (GPCRs) were employed in this research since the mutation data are available and could be used to retrospectively verify the prediction. The results show that the method presented in this article could pinpoint some retrospectively verified molecular determinants. The method is therefore suggested to be employed as a routine in drug design and discovery.

Agonist Binding and G Protein Coupling in Histamine H2 Receptor: A Molecular Dynamics Study

The histamine H2 receptor (H₂R) plays an important role in the regulation of gastric acid secretion. Therefore, it is a main drug target for the treatment of gastroesophageal reflux or peptic ulcer disease. However, there is as of yet no 3D-structural information available hampering a mechanistic understanding of H₂R. Therefore, we created a model of the histamine-H₂R-G_s complex based on the structure of the ternary complex of the β2-adrenoceptor and investigated the conformational stability of this active GPCR conformation. Since the physiologically relevant motions with respect to ligand binding and conformational changes of GPCRs can only partly be assessed on the timescale of conventional MD (cMD) simulations, we also applied metadynamics and Gaussian accelerated molecular dynamics (GaMD) simulations. A multiple walker metadynamics simulation in combination with cMD was applied for the determination of the histamine binding mode. The preferential binding pose detected is in good agreement with previous data from site directed mutagenesis and provides a basis for rational ligand design. Inspection of the H₂R-G_s interface reveals a network of polar interactions that may contribute to H₂R coupling selectivity. The cMD and GaMD simulations demonstrate that the active conformation is retained on a μs-timescale in the ternary histamine-H₂R-G_s complex and in a truncated complex that contains only G_s helix α5 instead of the entire G protein. In contrast, histamine alone is unable to stabilize the active conformation, which is in line with previous studies of other GPCRs.

PyPLIF HIPPOS: A Molecular Interaction Fingerprinting Tool for Docking Results of AutoDock Vina and PLANTS

We describe here our tool named PyPLIF HIPPOS, which was newly developed to analyze the docking results of AutoDock Vina and PLANTS. Its predecessor, PyPLIF (https://github.com/radifar/pyplif), is a molecular interaction fingerprinting tool for the docking results of PLANTS, exclusively. Unlike its predecessor, PyPLIF HIPPOS speeds up the computational times by separating the reference generation and docking analysis. PyPLIF HIPPOS also offers more options compared to PyPLIF. PyPLIF HIPPOS for Linux is stored as the Supporting Information in this application note and can be accessed in GitHub (https://github.com/radifar/PyPLIF-HIPPOS). Additionally, we present here the application of the tool in a retrospective structure-based virtual screening campaign targeting neuraminidase.

Structural Insight into the Mechanism of 4-Aminoquinolines Selectivity for the alpha2A-Adrenoceptor

Background

α_2A-adrenoceptor (AR) is a potential target for the treatment of degenerative diseases of the central nervous system, and α_2A-AR agonists are effective drugs for this condition. However, the lack of high selectivity for α_2A-AR subtype of traditional drugs greatly limits their clinic usage.

Methods

A series of homobivalent 4-aminoquinolines conjugated by two 4-aminoquinoline moieties via varying alkane linker length (C2-C12) were characterized for their affinities for each α₂-AR subtype. Subsequently, docking, molecular dynamics and mutagenesis were applied to uncover the molecular mechanism.

Results

Most 4-aminoquinolines (4-aminoquinoline monomer, C2-C6, C8-C10) were selective for α_2A-AR over α_2B- and α_2C-ARs. Besides, the affinities are of similar linker length-dependence for each α₂-AR subtype. Among all the compounds tested, C10 has the highest affinity for α_2A-AR (pKi=−7.45±0.62), which is 12-fold and 60-fold selective over α_2B-AR and α_2C-AR, respectively. Docking and molecular dynamics suggest that C10 simultaneously interacts with orthosteric and “allosteric” sites of the α_2A-AR. The mutation of F205 decreases the affinity by 2-fold. The potential allosteric residues include S90, N93, E94 and W99.

Conclusion

The specificity of C10 for the α_2A-AR and the potential orthosteric and allosteric binding sites proposed in this study provide valuable guidance for the development of novel α_2A-AR subtype selective compounds.

Distinct binding of cetirizine enantiomers to human serum albumin and the human histamine receptor H1

Cetirizine, a major metabolite of hydroxyzine, became a marketed second-generation H₁ antihistamine that is orally active and has a rapid onset of action, long duration of effects and a very good safety record at recommended doses. The approved drug is a racemic mixture of (S)-cetirizine and (R)-cetirizine, the latter being the levorotary enantiomer that also exists in the market as a third-generation, non-sedating and highly selective antihistamine. Both enantiomers bind tightly to the human histamine H₁ receptor (hH₁R) and behave as inverse agonists but the affinity and residence time of (R)-cetirizine are greater than those of (S)-cetirizine. In blood plasma, cetirizine exists in the zwitterionic form and more than 90% of the circulating drug is bound to human serum albumin (HSA), which acts as an inactive reservoir. Independent X-ray crystallographic work has solved the structure of the hH₁R:doxepin complex and has identified two drug-binding sites for cetirizine on equine serum albumin (ESA). Given this background, we decided to model a membrane-embedded hH₁R in complex with either (R)- or (S)-cetirizine and also the complexes of both ESA and HSA with these two enantiomeric drugs to analyze possible differences in binding modes between enantiomers and also among targets. The ensuing molecular dynamics simulations in explicit solvent and additional computational chemistry calculations provided structural and energetic information about all of these complexes that is normally beyond current experimental possibilities. Overall, we found very good agreement between our binding energy estimates and extant biochemical and pharmacological evidence. A much higher degree of solvent exposure in the cetirizine-binding site(s) of HSA and ESA relative to the more occluded orthosteric binding site in hH₁R is translated into larger positional fluctuations and considerably lower affinities for these two nonspecific targets. Whereas it is demonstrated that the two known pockets in ESA provide enough stability for cetirizine binding, only one such site does so in HSA due to a number of amino acid replacements. At the histamine-binding site in hH₁R, the distinct interactions established between the phenyl and chlorophenyl moieties of the two enantiomers with the amino acids lining up the pocket and between their free carboxylates and Lys179 in the second extracellular loop account for the improved pharmacological profile of (R)-cetirizine.

Phote-HrTH (Phormia terraenovae Hypertrehalosaemic Hormone), the Metabolic Hormone of the Fruit Fly: Solution Structure and Receptor Binding Model

Fruit flies are a widely distributed pest insect that pose a significant threat to food security. Flight is essential for the dispersal of the adult flies to find new food sources and ideal breeding spots. The supply of metabolic fuel to power the flight muscles of insects is regulated by adipokinetic hormones (AKHs). The fruit fly, Drosophila melanogaster, has the same AKH that is present in the blowfly, Phormia terraenovae; this AKH has the code-name Phote-HrTH. Binding of the AKH to the extra-cellular binding site of a G protein-coupled receptor causes its activation. In this paper, the structure of Phote-HrTH in sodium dodecyl sulfate (SDS) micelle solution was determined using NMR restrained molecular dynamics. The peptide was found to bind to the micelle and be fairly rigid, with an S2 order parameter of 0.96. The translated protein sequence of the AKH receptor from the fruit fly, D. melanogaster, Drome-AKHR, was used to construct two models of the receptor. It is proposed that these two models represent the active and inactive state of the receptor. The model based on the crystal structure of the β-2 adrenergic receptor was found to bind Phote-HrTH with a binding constant of −102kJmol−1, while the other model, based on the crystal structure of rhodopsin, did not bind the peptide. Under molecular dynamic simulation, in a palmitoyloleoylphosphatidylcholine (POPC) membrane, the receptor complex changed from an inactive to an active state. The identification and characterisation of the ligand binding site of Drome-AKHR provide novel information of ligand–receptor interaction, which could lead to the development of species-specific control substances to use discriminately against the fruit fly.

4-(3-Aminoazetidin-1-yl)pyrimidin-2-amines as High-Affinity Non-imidazole Histamine H3 Receptor Agonists with in Vivo Central Nervous System Activity

Despite the high diversity of histamine H3 receptor (H3R) antagonist/inverse agonist structures, partial or full H3R agonists have typically been imidazole derivatives. An in-house screening campaign intriguingly afforded the non-imidazole 4-(3-azetidin-1-yl)pyrimidin-2-amine 11b as a partial H3R agonist. Here, the design, synthesis, and structure-activity relationships of 11b analogues are described. This series yields several non-imidazole full agonists with potencies varying with the alkyl substitution pattern on the basic amine following the in vitro evaluation of H3R agonism using a cyclic adenosine monophosphate response element-luciferase reporter gene assay. The key compound VUF16839 (14d) combines nanomolar on-target activity (pKi = 8.5, pEC50 = 9.5) with weak activity on cytochrome P450 enzymes and good metabolic stability. The proposed H3R binding mode of 14d indicates key interactions similar to those attained by histamine. In vivo evaluation of 14d in a social recognition test in mice revealed an amnesic effect at 5 mg/kg intraperitoneally. The excellent in vitro and in vivo pharmacological profiles and the non-imidazole structure of 14d make it a promising tool compound in H3R research.

Covalent Inhibition of the Histamine H3 Receptor

Covalent binding of G protein-coupled receptors by small molecules is a useful approach for better understanding of the structure and function of these proteins. We designed, synthesized and characterized a series of 6 potential covalent ligands for the histamine H₃ receptor (H₃R). Starting from a 2-amino-pyrimidine scaffold, optimization of anchor moiety and warhead followed by fine-tuning of the required reactivity via scaffold hopping resulted in the isothiocyanate H₃R ligand 44. It shows high reactivity toward glutathione combined with appropriate stability in water and reacts selectively with the cysteine sidechain in a model nonapeptide equipped with nucleophilic residues. The covalent interaction of 44 with H₃R was validated with washout experiments and leads to inverse agonism on H₃R. Irreversible binder 44 (VUF15662) may serve as a useful tool compound to stabilize the inactive H₃R conformation and to study the consequences of prolonged inhibition of the H₃R.

Virtual screening-driven discovery of dual 5-HT6/5-HT2A receptor ligands with pro-cognitive properties

A virtual screening campaign aimed at finding structurally new compounds active at 5-HT₆R provided a set of candidates. Among those, one structure, 4-(5-{[(2-{5-fluoro-1H-pyrrolo[2,3-b]pyridin-3-yl}ethyl)amino]methyl}furan-2-yl)phenol (1, 5-HT₆R K_i = 91 nM), was selected as a hit for further optimization. As expected, the chemical scaffold of selected compound was significantly different from all the serotonin receptor ligands published to date. Synthetic efforts, supported by molecular modelling, provided 43 compounds representing different substitution patterns. The derivative 42, 4-(5-{[(2-{5-fluoro-1H-pyrrolo[2,3-b]pyridin-3-yl}ethyl)amino]methyl}furan-2-yl)phenol (5-HT₆R K_i = 25, 5-HT_2AR K_i = 32 nM), was selected as a lead and showed a good brain/plasma concentration profile, and it reversed phencyclidine-induced memory impairment. Considering the unique activity profile, the obtained series might be a good starting point for the development of a novel antipsychotic or antidepressant with pro-cognitive properties.

The adipokinetic hormones and their cognate receptor from the desert locust, Schistocerca gregaria: solution structure of endogenous peptides and models of their binding to the receptor

Background

Neuropeptides exert their activity through binding to G protein-coupled receptors (GPCRs). GPCRs are well-known drug targets in the pharmaceutical industry and are currently discussed as targets to control pest insects. Here, we investigate the neuropeptide adipokinetic hormone (AKH) system of the desert locust Schistocerca gregaria. The desert locust is known for its high reproduction, and for forming devastating swarms consisting of billions of individual insects. It is also known that S. gregaria produces three different AKHs as ligands but has only one AKH receptor (AKHR). The AKH system is known to be essential for metabolic regulation, which is necessary for reproduction and flight activity.

Methods

Nuclear magnetic resonance techniques (NMR) in a dodecylphosphocholin (DPC) micelle solution were used to determine the structure of the three AKHs. The primary sequence of the S. gregaria AKHR was used to construct a 3D molecular model. Next, the three AKHs were individually docked to the receptor, and dynamic simulation of the whole ligand–receptor complex in a model membrane was performed.

Results

Although the three endogenous AKHs of S. gregaria have quite different amino acids sequences and chain length (two octa- and one decapeptide), NMR experiments assigned a turn structure in DPC micelle solution for all. The GPCR-ModSim program identified human kappa opioid receptor to be the best template after which the S. gregaria AKHR was modeled. All three AKHs were found to have the same binding site on this receptor, interact with similar residues of the receptor and have comparable binding constants. Molecular switches were also identified; the movement of the receptor could be visually shown when ligands (AKHs) were docked and the receptor was activated.

Conclusions

The study proposes a model of binding of the three endogenous ligands to the one existing AKHR in the desert locust and paves the way to use such a model for the design of peptide analogs and finally, peptide mimetics, in the search for novel species-specific insecticides based on receptor–ligand interaction.

Dual 5-HT6 and D3 Receptor Antagonists in a Group of 1H-Pyrrolo[3,2-c]quinolines with Neuroprotective and Procognitive Activity

In light of the multifactorial origin of neurodegenerative disorders and some body of evidence indicating that pharmacological blockade of serotonin 5-HT₆ and dopamine D₃ receptors might be beneficial for cognitive decline, we envisioned (S)-1-[(3-chlorophenyl)sulfonyl]-4-(pyrrolidine-3-yl-amino)-1H-pyrrolo[3,2-c]quinoline (CPPQ), a neutral antagonist of 5-HT₆R, as a chemical template for designing dual antagonists of 5-HT₆/D₃ receptors. As shown by in vitro experiments, supported by quantum chemical calculations and molecular dynamic simulations, introducing alkyl substituents at the pyrrolidine nitrogen of CPPQ, fulfilled structural requirements for simultaneous modulation of 5-HT₆ and D₃ receptors. The study identified compound 19 ((S)-1-((3-chlorophenyl)sulfonyl)-N-(1-isobutylpyrrolidin-3-yl)-1H-pyrrolo[3,2-c]quinolin-4-amine), which was classified as a dual 5-HT₆/D₃R antagonist (K_i(5-HT6) = 27 nM, K_i(D3) = 7 nM). Compound 19 behaved as a neutral antagonist at G_s signaling and had no influence on receptor-operated, cyclin-dependent kinase 5 (Cdk5)-dependent neurite growth. In contrast to the well characterized 5-HT₆R antagonist intepirdine, compound 19 displayed neuroprotective properties against astrocyte damage induced by doxorubicin, as shown using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium (MTT) staining to assess cell metabolic activity and lactate dehydrogenase (LDH) release as an index of cell membrane disruption. This feature is of particular importance considering the involvement of loss of homeostatic function of glial cells in the progress of neurodegeneration. Biological results obtained for 19 in in vitro tests, translated into procognitive properties in phencyclidine (PCP)-induced memory decline in the novel object recognition (NOR) task in rats.

Route to Prolonged Residence Time at the Histamine H1 Receptor: Growing from Desloratadine to Rupatadine

Drug–target binding kinetics are an important predictor of in vivo drug efficacy, yet the relationship between ligand structures and their binding kinetics is often poorly understood. We show that both rupatadine (1) and desloratadine (2) have a long residence time at the histamine H₁ receptor (H₁R). Through development of a [³H]levocetirizine radiolabel, we find that the residence time of 1 exceeds that of 2 more than 10-fold. This was further explored with 22 synthesized rupatadine and desloratadine analogues. Methylene-linked cycloaliphatic or β-branched substitutions of desloratadine increase the residence time at the H₁R, conveying a longer duration of receptor antagonism. However, cycloaliphatic substituents directly attached to the piperidine amine (i.e., lacking the spacer) have decreased binding affinity and residence time compared to their methylene-linked structural analogues. Guided by docking studies, steric constraints within the binding pocket are hypothesized to explain the observed differences in affinity and binding kinetics between analogues.

Ligand-guided homology modeling drives identification of novel histamine H3 receptor ligands

In this study, we report a ligand-guided homology modeling approach allowing the analysis of relevant binding site residue conformations and the identification of two novel histamine H3 receptor ligands with binding affinity in the nanomolar range. The newly developed method is based on exploiting an essential charge interaction characteristic for aminergic G-protein coupled receptors for ranking 3D receptor models appropriate for the discovery of novel compounds through virtual screening.

Probe dependency in the determination of ligand binding kinetics at a prototypical G protein-coupled receptor

Drug-target binding kinetics are suggested to be important parameters for the prediction of in vivo drug-efficacy. For G protein-coupled receptors (GPCRs), the binding kinetics of ligands are typically determined using association binding experiments in competition with radiolabelled probes, followed by analysis with the widely used competitive binding kinetics theory developed by Motulsky and Mahan. Despite this, the influence of the radioligand binding kinetics on the kinetic parameters derived for the ligands tested is often overlooked. To address this, binding rate constants for a series of histamine H₁ receptor (H₁R) antagonists were determined using radioligands with either slow (low k_off) or fast (high k_off) dissociation characteristics. A correlation was observed between the probe-specific datasets for the kinetic binding affinities, association rate constants and dissociation rate constants. However, the magnitude and accuracy of the binding rate constant-values was highly dependent on the used radioligand probe. Further analysis using recently developed fluorescent binding methods corroborates the finding that the Motulsky-Mahan methodology is limited by the employed assay conditions. The presented data suggest that kinetic parameters of GPCR ligands depend largely on the characteristics of the probe used and results should therefore be viewed within the experimental context and limitations of the applied methodology.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

Structure-based exploration and pharmacological evaluation of N-substituted piperidin-4-yl-methanamine CXCR4 chemokine receptor antagonists

Using the available structural information of the chemokine receptor CXCR4, we present hit finding and hit exploration studies that make use of virtual fragment screening, design, synthesis and structure-activity relationship (SAR) studies. Fragment 2 was identified as virtual screening hit and used as a starting point for the exploration of 31 N-substituted piperidin-4-yl-methanamine derivatives to investigate and improve the interactions with the CXCR4 binding site. Additionally, subtle structural ligand changes lead to distinct interactions with CXCR4 resulting in a full to partial displacement of CXCL12 binding and competitive and/or non-competitive antagonism. Three-dimensional quantitative structure-activity relationship (3D-QSAR) and binding model studies were used to identify important hydrophobic interactions that determine binding affinity and indicate key ligand-receptor interactions.

Computer-Aided Studies for Novel Arylhydantoin 1,3,5-Triazine Derivatives as 5-HT₆ Serotonin Receptor Ligands with Antidepressive-Like, Anxiolytic and Antiobesity Action In Vivo

This study focuses on the design, synthesis, biological evaluation, and computer-aided structure-activity relationship (SAR) analysis for a novel group of aromatic triazine-methylpiperazines, with an hydantoin spacer between 1,3,5-traizine and the aromatic fragment. New compounds were synthesized and their affinities for serotonin 5-HT₆, 5-HT_1A, 5-HT_2A, 5-HT₇, and dopamine D₂ receptors were evaluated. The induced-fit docking (IFD) procedure was performed to explore the 5-HT₆ receptor conformation space employing two lead structures. It resulted in a consistent binding mode with the activity data. For the most active compounds found in each modification line, anti-obesity and anti-depressive-like activity in vivo, as well as "druglikeness" in vitro, were examined. Two 2-naphthyl compounds (18 and 26) were identified as the most active 5-HT₆R agents within each lead modification line, respectively. The 5-(2-naphthyl)hydantoin derivative 26, the most active one in the series (5-HT₆R: K_i = 87 nM), displayed also significant selectivity towards competitive G-protein coupled receptors (6⁻197-fold). Docking studies indicated that the hydantoin ring is stabilized by hydrogen bonding, but due to its different orientation, the hydrogen bonds form with S5.44 and N6.55 or Q6.58 for 18 and 26, respectively. Compound 26 exerted anxiolytic-like and antidepressant-like activities. Importantly, it demonstrated anti-obesity properties in animals fed palatable feed, and did not show toxic effects in vitro.

Discovery of a Novel Highly Selective Histamine H4 Receptor Antagonist for the Treatment of Atopic Dermatitis

The histamine H4 receptor (H4R), a member of the G-protein coupled receptor family, has been considered as a potential therapeutic target for treating atopic dermatitis (AD). A large number of H4R antagonists have been disclosed, but no efficient agents controlling both pruritus and inflammation in AD have been developed yet. Here, we have discovered a novel class of orally available H4R antagonists showing strong anti-itching and anti-inflammation activity as well as excellent selectivity against off-targets. A pharmacophore-based virtual screening system constructed in-house successfully identified initial hit compound 9, and the subsequent homology model-guided optimization efficiently led us to discover pyrido[2,3- e]tetrazolo[1,5- a]pyrazine analogue 48 as a novel chemotype of a potent and highly selective H4R antagonist. Importantly, orally administered compound 48 exhibits remarkable efficacy on antipruritus and anti-inflammation with a favorable pharmacokinetic (PK) profile in several mouse models of AD. Thus, these data strongly suggest that our compound 48 is a promising clinical candidate for treatment of AD.

Chemical Diversity in the G Protein-Coupled Receptor Superfamily

G protein-coupled receptors (GPCRs) are the largest family of cell signaling transmembrane proteins that can be modulated by a plethora of chemical compounds. Systematic cheminformatics analysis of structurally and pharmacologically characterized GPCR ligands shows that cocrystallized GPCR ligands cover a significant part of chemical ligand space, despite their limited number. Many GPCR ligands and substructures interact with multiple receptors, providing a basis for polypharmacological ligand design. Experimentally determined GPCR structures represent a variety of binding sites and receptor-ligand interactions that can be translated to chemically similar ligands for which structural data are lacking. This integration of structural, pharmacological, and chemical information on GPCR-ligand interactions enables the extension of the structural GPCR-ligand interactome and the structure-based design of novel modulators of GPCR function.

3D-e-Chem: Structural Cheminformatics Workflows for Computer-Aided Drug Discovery

eScience technologies are needed to process the information available in many heterogeneous types of protein-ligand interaction data and to capture these data into models that enable the design of efficacious and safe medicines. Here we present scientific KNIME tools and workflows that enable the integration of chemical, pharmacological, and structural information for: i) structure-based bioactivity data mapping, ii) structure-based identification of scaffold replacement strategies for ligand design, iii) ligand-based target prediction, iv) protein sequence-based binding site identification and ligand repurposing, and v) structure-based pharmacophore comparison for ligand repurposing across protein families. The modular setup of the workflows and the use of well-established standards allows the re-use of these protocols and facilitates the design of customized computer-aided drug discovery workflows.

Synthesis and Characterization of a Bidirectional Photoswitchable Antagonist Toolbox for Real-Time GPCR Photopharmacology

Noninvasive methods to modulate G protein-coupled receptors (GPCRs) with temporal and spatial precision are in great demand. Photopharmacology uses photons to control in situ the biological properties of photoswitchable small-molecule ligands, which bodes well for chemical biological precision approaches. Integrating the light-switchable configurational properties of an azobenzene into the ligand core, we developed a bidirectional antagonist toolbox for an archetypical family A GPCR, the histamine H₃ receptor (H₃R). From 16 newly synthesized photoswitchable compounds, VUF14738 (28) and VUF14862 (33) were selected as they swiftly and reversibly photoisomerize and show over 10-fold increased or decreased H₃R binding affinities, respectively, upon illumination at 360 nm. Both ligands combine long thermal half-lives with fast and high photochemical trans-/cis conversion, allowing their use in real-time electrophysiology experiments with oocytes to confirm dynamic photomodulation of H₃R activation in repeated second-scale cycles. VUF14738 and VUF14862 are robust and fatigue-resistant photoswitchable GPCR antagonists suitable for spatiotemporal studies of H₃R signaling.

Comparative molecular field analysis and molecular dynamics studies of the dopamine D2 receptor antagonists without a protonatable nitrogen atom

The dopaminergic hypothesis of schizophrenia is the main concept explaining the direct reasons of schizophrenia and the effectiveness of current antipsychotics. All antipsychotics present on the market are potent dopamine D₂ receptor antagonists or partial agonists. In this work we investigate a series of dopamine D₂ receptor antagonists which do not fulfill the criteria of the classical pharmacophore model as they do not possess a protonatable nitrogen atom necessary to interact with the conserved Asp(3.32). Such compounds are interesting, inter alia, due to possible better pharmacokinetic profile when compared to basic, ionizable molecules. By means of homology modeling, molecular docking and molecular dynamics we determined that the compounds investigated interact with Asp(3.32) via their amide nitrogen atom. It was found that the studied compounds stabilize the receptor inactive conformation through the effect on the ionic lock, which is typical for GPCR antagonists. We constructed a CoMFA model for the studied compounds with the following statistics: R² = 0.95, Q² = 0.63. The quality of the CoMFA model was confirmed by high value of R² of the test set, equal 0.96. The CoMFA model indicated two regions where bulky substituents are favored and two regions where bulky substituents are not beneficial. Two red contour regions near carbonyl groups were identified meaning that negative charge would be favored here. Furthermore, the S-oxide group is connected with blue contour region meaning that positive charge is favored in this position. These findings may be applied for further optimization of the studied compound series.

Identifying human diamine sensors for death related putrescine and cadaverine molecules

Pungent chemical compounds originating from decaying tissue are strong drivers of animal behavior. Two of the best-characterized death smell components are putrescine (PUT) and cadaverine (CAD), foul-smelling molecules produced by decarboxylation of amino acids during decomposition. These volatile polyamines act as ‘necromones’, triggering avoidance or attractive responses, which are fundamental for the survival of a wide range of species. The few studies that have attempted to identify the cognate receptors for these molecules have suggested the involvement of the seven-helix trace amine-associated receptors (TAARs), localized in the olfactory epithelium. However, very little is known about the precise chemosensory receptors that sense these compounds in the majority of organisms and the molecular basis of their interactions. In this work, we have used computational strategies to characterize the binding between PUT and CAD with the TAAR6 and TAAR8 human receptors. Sequence analysis, homology modeling, docking and molecular dynamics studies suggest a tandem of negatively charged aspartates in the binding pocket of these receptors which are likely to be involved in the recognition of these small biogenic diamines.

The distinctive dead smell comes largely from molecules like cadaverine and putrescine that are produced during decomposition of organic tissues. These volatile compounds act as powerful chemical signals important for the survival of a wide range of species. Previous studies have identified the trace amine-associated receptor 13c (or TAAR13c) in zebrafish as the cognate receptor of cadaverine in bony fishes. In this work, we employed computational strategies to disclose the human TAAR6 and TAAR8 receptors as sensors of the putrescine and cadaverine molecules. Our results indicate that several negatively charged residues in the ligand binding pocket of these receptors constitute the molecular basis for recognition of these necromones in humans.

A Structural Framework for GPCR Chemogenomics: What's In a Residue Number?

The recent surge of crystal structures of G protein-coupled receptors (GPCRs), as well as comprehensive collections of sequence, structural, ligand bioactivity, and mutation data, has enabled the development of integrated chemogenomics workflows for this important target family. This chapter will focus on cross-family and cross-class studies of GPCRs that have pinpointed the need for, and the implementation of, a generic numbering scheme for referring to specific structural elements of GPCRs. Sequence- and structure-based numbering schemes for different receptor classes will be introduced and the remaining caveats will be discussed. The use of these numbering schemes has facilitated many chemogenomics studies such as consensus binding site definition, binding site comparison, ligand repurposing (e.g. for orphan receptors), sequence-based pharmacophore generation for homology modeling or virtual screening, and class-wide chemogenomics studies of GPCRs.

From the Promiscuous Asenapine to Potent Fluorescent Ligands Acting at a Series of Aminergic G-Protein-Coupled Receptors

Monoamine neurotransmitters such as serotonin, dopamine, histamine, and noradrenaline have important and varied physiological functions and similar chemical structures. Representing important pharmaceutical drug targets, the corresponding G-protein-coupled receptors (termed aminergic GPCRs) belong to the class of cell membrane receptors and share many levels of similarity as well. Given their pharmacological and structural closeness, one could hypothesize the possibility to derivatize a ubiquitous ligand to afford rapidly fluorescent probes for a large set of GPCRs to be used for instance in FRET-based binding assays. Here we report fluorescent derivatives of the nonselective agent asenapine which were designed, synthesized, and evaluated as ligands of 34 serotonin, dopamine, histamine, melatonin, acetylcholine, and adrenergic receptors. It appears that this strategy led rapidly to the discovery and development of nanomolar affinity fluorescent probes for 14 aminergic GPCRs. Selected probes were tested in competition binding assays with unlabeled competitors in order to demonstrate their suitability for drug discovery purposes.

The Target Residence Time of Antihistamines Determines Their Antagonism of the G Protein-Coupled Histamine H1 Receptor

The pharmacodynamics of drug-candidates is often optimized by metrics that describe target binding (K_d or K_i value) or target modulation (IC₅₀). However, these metrics are determined at equilibrium conditions, and consequently information regarding the onset and offset of target engagement and modulation is lost. Drug-target residence time is a measure for the lifetime of the drug-target complex, which has recently been receiving considerable interest, as target residence time is shown to have prognostic value for the in vivo efficacy of several drugs. In this study, we have investigated the relation between the increased residence time of antihistamines at the histamine H₁ receptor (H₁R) and the duration of effective target-inhibition by these antagonists. Hela cells, endogenously expressing low levels of the H₁R, were incubated with a series of antihistamines and dissociation was initiated by washing away the unbound antihistamines. Using a calcium-sensitive fluorescent dye and a label free, dynamic mass redistribution based assay, functional recovery of the H₁R responsiveness was measured by stimulating the cells with histamine over time, and the recovery was quantified as the receptor recovery time. Using these assays, we determined that the receptor recovery time for a set of antihistamines differed more than 40-fold and was highly correlated to their H₁R residence times, as determined with competitive radioligand binding experiments to the H₁R in a cell homogenate. Thus, the receptor recovery time is proposed as a cell-based and physiologically relevant metric for the lead optimization of G protein-coupled receptor antagonists, like the H₁R antagonists. Both, label-free or real-time, classical signaling assays allow an efficient and physiologically relevant determination of kinetic properties of drug molecules.

Interaction of the red pigment-concentrating hormone of the crustacean Daphnia pulex, with its cognate receptor, Dappu-RPCHR: A nuclear magnetic resonance and modeling study

The primary sequence of the red pigment-concentrating hormone (RPCH) receptor of the water flea, Daphnia pulex, was used in homology modeling to construct the first 3D model of a crustacean G-protein coupled receptor, Dappu-RPCHR. This receptor was found to belong to the class A subfamily of GPCRs with a disulfide bridge between Cys⁷² and Cys¹⁵⁰ and an ionic lock between Arg⁹⁷ and Thr²²⁴ and Thr²²⁰. NMR restrained molecular dynamics was used to determine the structure of an agonist, Dappu-RPCH, in a membrane-mimicking environment. The agonist was found to be flexible but has two main conformations in solution, both having β-turns. Docking of the predominant structure was used to find a binding pocket on the receptor. The pocket's spatial location was similar to that of the AKH receptor of Anopheles gambiae. The binding affinity was -69kcalmol^-1 with the N-terminus of Dappu-RPCH inserted between helices 4 and 6, and the C-terminus interacting with extra-cellular loop, ECL2. Upon binding, H-bonding to the peptide may activate the receptor. This development of the first Dappu-RPCH/Dappu-RPCHR model could be useful for understanding ligand-receptor interactions in crustaceans.

Structure and Function of Peptide-Binding G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors and are important human drug targets. Of the 826 human GPCRs, 118 of them recognize endogenous peptide or protein ligands, and 30 of the 118 are targeted by approved drug molecules, including the very high-profile class B glucagon-like peptide 1 receptor. In this review, we analyze the 21 experimentally determined three-dimensional structures of the known peptide-binding GPCRs in relation to the endogenous peptides and drug molecules that modulate their cell signaling processes. Our integrated analyses reveal that half of the marketed drugs and most of the drugs in clinical trials that interact with peptide GPCRs are small molecules with a wide range of binding modes distinct from those of large peptide ligands. As we continue to collect additional data on these receptors from orthogonal approaches, including nuclear magnetic resonance and electron microscopy, we are beginning to understand how these receptors interact with their ligands at the molecular level and how improving the pharmacology of GPCR signal transduction requires us to study these receptors using multiple biophysical techniques.

3D-e-Chem-VM: Structural Cheminformatics Research Infrastructure in a Freely Available Virtual Machine

3D-e-Chem-VM is an open source, freely available Virtual Machine (http://3d-e-chem.github.io/3D-e-Chem-VM/) that integrates cheminformatics and bioinformatics tools for the analysis of protein–ligand interaction data. 3D-e-Chem-VM consists of software libraries, and database and workflow tools that can analyze and combine small molecule and protein structural information in a graphical programming environment. New chemical and biological data analytics tools and workflows have been developed for the efficient exploitation of structural and pharmacological protein–ligand interaction data from proteomewide databases (e.g., ChEMBLdb and PDB), as well as customized information systems focused on, e.g., G protein-coupled receptors (GPCRdb) and protein kinases (KLIFS). The integrated structural cheminformatics research infrastructure compiled in the 3D-e-Chem-VM enables the design of new approaches in virtual ligand screening (Chemdb4VS), ligand-based metabolism prediction (SyGMa), and structure-based protein binding site comparison and bioisosteric replacement for ligand design (KRIPOdb).

Molecular Modelling Approaches for the Analysis of Histamine Receptors and Their Interaction with Ligands

Several experimental techniques to analyse histamine receptors are available, e.g. pharmacological characterisation of known or new compounds by different types of assays or mutagenesis studies. To obtain insights into the histamine receptors on a molecular and structural level, crystal structures have to be determined and molecular modelling studies have to be performed. It is widely accepted to generate homology models of the receptor of interest based on an appropriate crystal structure as a template and to refine the resulting models by molecular dynamic simulations. A lot of modelling techniques, e.g. docking, QSAR or interaction fingerprint methods, are used to predict binding modes of ligands and pharmacological data, e.g. affinity or even efficacy. However, within the last years, molecular dynamic simulations got more and more important: First of all, molecular dynamic simulations are very helpful to refine the binding mode of a ligand to a histamine receptor, obtained by docking studies. Furthermore, with increasing computational performance it got possible to simulate complete binding pathways of ions or ligands from the aqueous extracellular phase into the allosteric or orthosteric binding pocket of histamine receptors.

Molecular interaction fingerprint approaches for GPCR drug discovery

Protein-ligand interaction fingerprints (IFPs) are binary 1D representations of the 3D structure of protein-ligand complexes encoding the presence or absence of specific interactions between the binding pocket amino acids and the ligand. Various implementations of IFPs have been developed and successfully applied for post-processing molecular docking results for G Protein-Coupled Receptor (GPCR) ligand binding mode prediction and virtual ligand screening. Novel interaction fingerprint methods enable structural chemogenomics and polypharmacology predictions by complementing the increasing amount of GPCR structural data. Machine learning methods are increasingly used to derive relationships between bioactivity data and fingerprint descriptors of chemical and structural information of binding sites, ligands, and protein-ligand interactions. Factors that influence the application of IFPs include structure preparation, binding site definition, fingerprint similarity assessment, and data processing and these factors pose challenges as well possibilities to optimize interaction fingerprint methods for GPCR drug discovery.

Function-specific virtual screening for GPCR ligands using a combined scoring method

The ability of scoring functions to correctly select and rank docking poses of small molecules in protein binding sites is highly target dependent, which presents a challenge for structure-based drug discovery. Here we describe a virtual screening method that combines an energy-based docking scoring function with a molecular interaction fingerprint (IFP) to identify new ligands based on G protein-coupled receptor (GPCR) crystal structures. The consensus scoring method is prospectively evaluated by: 1) the discovery of chemically novel, fragment-like, high affinity histamine H₁ receptor (H₁R) antagonists/inverse agonists, 2) the selective structure-based identification of ß₂-adrenoceptor (ß₂R) agonists, and 3) the experimental validation and comparison of the combined and individual scoring approaches. Systematic retrospective virtual screening simulations allowed the definition of scoring cut-offs for the identification of H₁R and ß₂R ligands and the selection of an optimal ß-adrenoceptor crystal structure for the discrimination between ß₂R agonists and antagonists. The consensus approach resulted in the experimental validation of 53% of the ß₂R and 73% of the H₁R virtual screening hits with up to nanomolar affinities and potencies. The selective identification of ß₂R agonists shows the possibilities of structure-based prediction of GPCR ligand function by integrating protein-ligand binding mode information.

Polypharmacology: in silico methods of ligand design and development

How to design a ligand to bind multiple targets, rather than to a single target, is the focus of this review. Rational polypharmacology draws on knowledge that is both broad ranging and hierarchical. Computer-aided multitarget ligand design methods are described according to their nested knowledge level. Ligand-only and then receptor-ligand strategies are first described; followed by the metabolic network viewpoint. Subsequently strategies that view infectious diseases as multigenomic targets are discussed, and finally the disease level interpretation of medicinal therapy is considered. As yet there is no consensus on how best to proceed in designing a multitarget ligand. The current methodologies are bought together in an attempt to give a practical overview of how polypharmacology design might be best initiated.

International Union of Basic and Clinical Pharmacology. XCVIII. Histamine Receptors

Histamine is a developmentally highly conserved autacoid found in most vertebrate tissues. Its physiological functions are mediated by four 7-transmembrane G protein-coupled receptors (H1R, H2R, H3R, H4R) that are all targets of pharmacological intervention. The receptors display molecular heterogeneity and constitutive activity. H1R antagonists are long known antiallergic and sedating drugs, whereas the H2R was identified in the 1970s and led to the development of H2R-antagonists that revolutionized stomach ulcer treatment. The crystal structure of ligand-bound H1R has rendered it possible to design new ligands with novel properties. The H3R is an autoreceptor and heteroreceptor providing negative feedback on histaminergic and inhibition on other neurons. A block of these actions promotes waking. The H4R occurs on immuncompetent cells and the development of anti-inflammatory drugs is anticipated.