Generation of a homology model of the human histamine H3 receptor for ligand docking and pharmacophore-based screening

Efficiency of Homology Modeling Assisted Molecular Docking in G-protein Coupled Receptors

BACKGROUND

Molecular docking is in regular practice to assess ligand affinity on a target protein crystal structure. In the absence of protein crystal structure, the homology modeling or comparative modeling is the best alternative to elucidate the relationship details between a ligand and protein at the molecular level. The development of accurate homology modeling (HM) and its integration with molecular docking (MD) is essential for successful, rational drug discovery.

OBJECTIVE

The G-protein coupled receptors (GPCRs) are attractive therapeutic targets due to their immense role in human pharmacology. The GPCRs are membrane-bound proteins with the complex constitution, and the understanding of their activation and inactivation mechanisms is quite challenging. Over the past decade, there has been a rapid expansion in the number of solved G-protein-coupled receptor (GPCR) crystal structures; however, the majority of the GPCR structures remain unsolved. In this context, HM guided MD has been widely used for structure-based drug design (SBDD) of GPCRs.

METHODS

The focus of this review is on the recent (i) developments on HM supported GPCR drug discovery in the absence of GPCR crystal structures and (ii) application of HM in understanding the ligand interactions at the binding site, virtual screening, determining receptor subtype selectivity and receptor behaviour in comparison with GPCR crystal structures.

RESULTS

The HM in GPCRs has been extremely challenging due to the scarcity in template structures. In such a scenario, it is difficult to get accurate HM that can facilitate understanding of the ligand-receptor interactions. This problem has been alleviated to some extent by developing refined HM based on incorporating active /inactive ligand information and inducing protein flexibility. In some cases, HM proteins were found to outscore crystal structures.

CONCLUSION

The developments in HM have been highly operative to gain insights about the ligand interaction at the binding site and receptor functioning at the molecular level. Thus, HM guided molecular docking may be useful for rational drug discovery for the GPCRs mediated diseases.

Structural modifications in the distal, regulatory region of histamine H3 receptor antagonists leading to the identification of a potent anti-obesity agent

A series of 4-pyridylpiperazine derivatives with varying regulatory region substituents proved to be potent histamine H3 receptor (H3R) ligands in the nanomolar concentration range. The most influential modification that affected the affinity toward the H3R appeared by introducing electron-withdrawing moieties into the distal aromatic ring. In order to finally discuss the influence of the characteristic 4-pyridylpiperazine moiety on H3R affinity, two Ciproxifan analogues 2 and 3 with a slight modification in their basic part were obtained. The replacement of piperazine in 3 with piperidine in compound 2, led to slightly reduced affinity towards the H3R (Ki = 3.17 and 7.70 nM, respectively). In fact, 3 showed the highest antagonistic properties among all compounds in this series, hence affirming our previous assumptions, that the 4-pyridylpiperazine moiety is the key element for suitable interaction with the human histamine H3 receptor. While its structural replacement to piperidine is also tolerated for H3R binding, the heteroaromatic 4-pyridyl moiety seems to be essential for proper ligand-receptor interaction. The putative protein-ligand interactions responsible for their high affinity were demonstrated using molecular modeling techniques. Furthermore, selectivity, intrinsic activity at the H3R, as well as drug-like properties of ligands were evaluated using in vitro methods. Moreover, pharmacological in vivo test results of compound 9 (structural analogue of Abbott's A-331440) clearly indicate that it may affect the amount of calories consumed, thus act as an anorectic compound.

Novel potent (dihydro)benzofuranyl piperazines as human histamine receptor ligands - Functional characterization and modeling studies on H3 and H4 receptors

Histamine acts through four different receptors (H₁R-H₄R), the H₃R and H₄R being the most explored in the last years as drug targets. The H₃R is a potential target to treat narcolepsy, Parkinson's disease, epilepsy, schizophrenia and several other CNS-related conditions, while H₄R blockade leads to anti-inflammatory and immunomodulatory effects. Our group has been exploring the dihydrobenzofuranyl-piperazines (LINS01 series) as human H₃R/H₄R ligands as potential drug candidates. In the present study, a set of 12 compounds were synthesized from adequate (dihydro)benzofuran synthons through simple reactions with corresponding piperazines, giving moderate to high yields. Four compounds (1b, 1f, 1g and 1h) showed high hH₃R affinity (pK_i > 7), compound 1h being the most potent (pK_i 8.4), and compound 1f showed the best efficiency (pKi 8.2, LE 0.53, LLE 5.85). BRET-based assays monitoring Gα_i activity indicated that the compounds are potent antagonists. Only one compound (2c, pK_i 7.1) presented high affinity for hH₄R. In contrast to what was observed for hH₃R, it showed partial agonist activity. Docking experiments indicated that bulky substituents occupy a hydrophobic pocket in hH₃R, while the N-allyl group forms favorable interactions with hydrophobic residues in the TM2, 3 and 7, increasing the selectivity towards hH₃R. Additionally, the importance of the indole NH in the interaction with Glu5.46 from hH₄R was confirmed by the modeling results, explaining the affinity and agonistic activity of compound 2c. The data reported in this work represent important findings for the rational design of future compounds for hH₃R and hH₄R.

Molecular Modelling of Potential Candidates for the Treatment of Depression

A lot of research initiatives in the last decades have been focused on the search of new strategies to treat depression. However, despite the availability of various antidepressants, current treatment is still far from ideal. Unwanted side effects, modest response rates and the slow onset of action are the main shortcomings. As a strategy to improve symptomatic relief and response rates, the dual modulation of the serotonin transporter and the histamine H₃ receptor by a single chemical entity has been proposed in the literature. Accordingly, this work aims to elucidate key structural features responsible for the dual inhibitory activity of the hexahydro-pyrrolo-isoquinoline derivatives. For this purpose, two approaches were employed, four-dimensional quantitative structure-activity relationship (4D-QSAR) and molecular docking. The 4D-QSAR models for both receptors allowed the identification of the pharmacophore groups critical for the modelled biological activity, whereas the binding mode of this class of compounds to the human serotonin transporter was assessed by molecular docking. The findings can be applicable to design new antidepressants.

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.

Design, synthesis, and biological evaluation of novel oxadiazole- and thiazole-based histamine H3R ligands

Histamine H₃ receptor (H₃R) is largely expressed in the CNS and modulation of the H₃R function can affect histamine synthesis and liberation, and modulate the release of many other neurotransmitters. Targeting H₃R with antagonists/inverse agonists may have therapeutic applications in neurodegenerative disorders, gastrointestinal and inflammatory diseases. This prompted us to design and synthesize azole-based H₃R ligands, i.e. having oxadiazole- or thiazole-based core structures. While ligands of oxadiazole scaffold were almost inactive, thiazole-based ligands were very potent and several exhibited binding affinities in a nanomolar concentration range. Ligands combining 4-cyanophenyl moiety as arbitrary region, para-xylene or piperidine carbamoyl linkers, and/or pyrrolidine or piperidine basic heads were found to be the most active within this series of thiazole-based H₃R ligands. The most active ligands were in silico screened for ADMET properties and drug-likeness. They fulfilled Lipinski's and Veber's rules and exhibited potential activities for oral administration, blood-brain barrier penetration, low hepatotoxicity, combined with an overall good toxicity profile.

Synthesis and biological activity of novel tert-amylphenoxyalkyl (homo)piperidine derivatives as histamine H3R ligands

As a continuation of our search for novel histamine H₃ receptor ligands a series of twenty new tert-amyl phenoxyalkylamine derivatives (2-21) was synthesized. Compounds of four to eight carbon atoms spacer alkyl chain were evaluated on their binding properties at human histamine H₃ receptor (hH₃R). The highest affinities were observed for pentyl derivatives 6-8 (K_i=8.8-23.4nM range) and among them piperidine derivative 6 with K_i=8.8nM. Structures 6, 7 were also classified as antagonists in cAMP accumulation assay (with EC₅₀=157 and 164nM, respectively). Moreover, new compounds were also evaluated for anticonvulsant activity in Antiepileptic Screening Program (ASP) at National Institute of Neurological Disorders and Stroke (USA). Seven compounds (2-4, 9, 11, 12 and 20) showed anticonvulsant activity at maximal electroshock (MES) test in the dose of 30mg/kg at 0.5h. In the subcutaneous pentetrazole (scMET) test compound 4 showed protection at 100 and 300mg/kg dose at mice, however compounds showed high neurotoxicity in rotarod test at used doses. Also, molecular modeling studies were undertaken, to explain affinity of compounds at hH₃R (taking into the consideration X-ray analysis of compound 18). In order to estimate "drug-likeness" of selected compounds in silico and experimental evaluation of lipophilicity, metabolic stability and cytotoxicity was performed.

Integrated in Silico Methods for the Design and Optimization of Novel Drug Candidates

Although almost fully automated, the discovery of novel, effective, and safe drugs is still a long-term and highly expensive process. Consequently, the need for fleet, rational, and cost-efficient development of novel drugs is crucial, and nowadays the advanced in silico drug design methodologies seem to effectively meet these issues. The aim of this chapter is to provide a comprehensive overview of some of the current trends and advances in the in silico design of novel drug candidates with a special emphasis on 6-fluoroquinolone (6-FQ) antibacterials as potential novel Mycobacterium tuberculosis DNA gyrase inhibitors. In particular, the chapter covers some of the recent aspects of a wide range of in silico drug discovery approaches including multidimensional machine-learning methods, ligand-based and structure-based methodologies, as well as their proficient combination and integration into an intelligent virtual screening protocol for design and optimization of novel 6-FQ analogs.

Cholinesterase inhibitory activity of chlorophenoxy derivatives-Histamine H3 receptor ligands

In recent years, multitarget-directed ligands have become an interesting strategy in a search for a new treatment of Alzheimer's disease. Combination of both: a histamine H3 receptor antagonist/inverse agonist and a cholinesterases inhibitor in one molecule could provide a new therapeutic opportunity. Here, we present biological evaluation of histamine H3 receptor ligands-chlorophenoxyalkylamine derivatives against cholinesterases: acetyl- and butyrylcholinesterase. The target compounds showed cholinesterase inhibitory activity in a low micromolar range. The most potent in this group was 1-(7-(4-chlorophenoxy)heptyl)homopiperidine (18) inhibiting the both enzymes (EeAChE IC50=1.93μM and EqBuChE IC50=1.64μM). Molecular modeling studies were performed to explain the binding mode of 18 with histamine H3 receptor as well as with cholinesterases.

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.

Chlorophenoxy aminoalkyl derivatives as histamine H(3)R ligands and antiseizure agents

A series of twenty new chlorophenoxyalkylamine derivatives (9-28) was synthesized and evaluated on their binding properties at the human histamine H3 receptor (hH3R). The spacer alkyl chain contained five to seven carbon atoms. The highest affinities have shown the 4-chloro substituted derivatives 10 and 25 (Ki=133 and 128 nM, respectively) classified as antagonists in cAMP accumulation assay (EC50=72 and 75 nM, respectively). Synthesized compounds were also evaluated for anticonvulsant activity in Antiepileptic Screening Program (ASP) at National Institute of Neurological Disorders and Stroke (USA). Two compounds (4-chloro substituted derivatives: 20 and 26) were the most promising and showed in the MES seizure model in rats (after ip administration) ED50 values of 14 mg/kg and 13.18 mg/kg, respectively. Protective indexes (PI=TD50/ED50) were 3.2 for 20 and 3.8 for 26. Moreover, molecular modeling and docking studies were undertaken to explain affinity at hH3R of target compounds, and the experimentally and in silico estimation of properties like lipophilicity and metabolism was performed. Antiproliferative effects have been also investigated in vitro for selected compounds (10 and 25). These compounds neither possessed significant antiproliferative and antitumor activity, nor modulated CYP3A4 activity up to concentration of 10 μM.

Cherry-picked ligands at histamine receptor subtypes

Histamine, a biogenic amine, is considered as a principle mediator of multiple physiological effects through binding to its H1, H2, H3, and H4 receptors (H1-H4Rs). Currently, the HRs have gained attention as important targets for the treatment of several diseases and disorders ranging from allergy to Alzheimer's disease and immune deficiency. Accordingly, medicinal chemistry studies exploring histamine-like molecules and their physicochemical properties by binding and interacting with the four HRs has led to the development of a diversity of agonists and antagonists that display selectivity for each HR subtype. An overview on H1-R4Rs and developed ligands representing some key steps in development is provided here combined with a short description of structure-activity relationships for each class. Main chemical diversities, pharmacophores, and pharmacological profiles of most innovative H1-H4R agonists and antagonists are highlighted. Therefore, this overview should support the rational choice for the optimal ligand selection based on affinity, selectivity and efficacy data in biochemical and pharmacological studies. This article is part of the Special Issue entitled 'Histamine Receptors'.

Ligand orientation in a membrane-embedded receptor site revealed by solid-state NMR with paramagnetic relaxation enhancement

Paramagnetic relaxation-enhanced solid-state NMR reveals a ouabain analogue with an inverted orientation in the Na,K-ATPase inhibitory site.

Scaffold variations in amine warhead of histamine H₃ receptor antagonists

The histamine H₃ receptor (H₃R) is involved in numerous regulatory neurotransmission processes and there-fore, is a prominent target for centrally occurring disease with some promising clinical candidates. Previous research resulted in the identification of a core pharmacophore blueprint for H₃R antagonists/inverse agonists, which when inserted in a molecule, mostly ensures acceptable affinity. Nevertheless, variations of scaffold and peripheral areas can increase potency and pharmacokinetic profile of drug candidates. The variations in amine scaffolds of antagonists for this aminergic GPCR are of special importance.

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

BACKGROUND AND PURPOSE

Chemogenomics focuses on the discovery of new connections between chemical and biological space leading to the discovery of new protein targets and biologically active molecules. G-protein coupled receptors (GPCRs) are a particularly interesting protein family for chemogenomics studies because there is an overwhelming amount of ligand binding affinity data available. The increasing number of aminergic GPCR crystal structures now for the first time allows the integration of chemogenomics studies with high-resolution structural analyses of GPCR-ligand complexes.

EXPERIMENTAL APPROACH

In this study, we have combined ligand affinity data, receptor mutagenesis studies, and amino acid sequence analyses to high-resolution structural analyses of (hist)aminergic GPCR-ligand interactions. This integrated structural chemogenomics analysis is used to more accurately describe the molecular and structural determinants of ligand affinity and selectivity in different key binding regions of the crystallized aminergic GPCRs, and histamine receptors in particular.

KEY RESULTS

Our investigations highlight interesting correlations and differences between ligand similarity and ligand binding site similarity of different aminergic receptors. Apparent discrepancies can be explained by combining detailed analysis of crystallized or predicted protein-ligand binding modes, receptor mutation studies, and ligand structure-selectivity relationships that identify local differences in essential pharmacophore features in the ligand binding sites of different receptors.

CONCLUSIONS AND IMPLICATIONS

We have performed structural chemogenomics studies that identify links between (hist)aminergic receptor ligands and their binding sites and binding modes. This knowledge can be used to identify structure-selectivity relationships that increase our understanding of ligand binding to (hist)aminergic receptors and hence can be used in future GPCR ligand discovery and design.

Mathematical analysis of the sodium sensitivity of the human histamine H3 receptor

Purpose

It was shown by several experimental studies that some G protein coupled receptors (GPCR) are sensitive to sodium ions. Furthermore, mutagenesis studies or the determination of crystal structures of the adenosine A_2A or δ-opioid receptor revealed an allosteric Na⁺ binding pocket near to the highly conserved Asp^2.50. Within a previous study, the influence of NaCl concentration onto the steady-state GTPase activity at the human histamine H₃ receptor (hH₃R) in presence of the endogenous histamine or the inverse agonist thioperamide was analyzed. The purpose of the present study was to examine and quantify the Na⁺-sensitivity of hH₃R on a molecular level.

Methods

To achieve this, we developed a set of equations, describing constitutive activity and the different ligand-receptor equilibria in absence or presence of sodium ions. Furthermore, in order to gain a better understanding of the ligand- and Na⁺-binding to hH₃R on molecular level, we performed molecular dynamic (MD) simulations.

Results

The analysis of the previously determined experimental steady-state GTPase data with the set of equations presented within this study, reveals that thioperamide binds into the orthosteric binding pocket of the hH₃R in absence or presence of a Na⁺ in its allosteric binding site. However, the data suggest that thioperamide binds preferentially into the hH₃R in absence of a sodium ion in its allosteric site. These experimental results were supported by MD simulations of thioperamide in the binding pocket of the inactive hH₃R. Furthermore, the MD simulations revealed two different binding modes for thioperamide in presence or absence of a Na⁺ in its allosteric site.

Conclusion

The mathematical model presented within this study describes the experimental data regarding the Na⁺-sensitivity of hH₃R in an excellent manner. Although the present study is focused onto the Na⁺-sensitivity of the hH₃R, the resulting equations, describing Na⁺- and ligand-binding to a GPCR, can be used for all other ion-sensitive GPCRs.

Procognitive properties of drugs with single and multitargeting H3 receptor antagonist activities

The histamine H3 receptor (H3 R) is an important modulator of numerous central control mechanisms. Novel lead optimizations for H3 R antagonists/inverse agonists involved studies of structure-activity relationships, cross-affinities, and pharmacokinetic properties of promising ligands. Blockade of inhibitory histamine H3 autoreceptors reinforces histaminergic transmission, while antagonism of H3 heteroreceptors accelerates the corticolimbic liberation of acetylcholine, norepinephrine, glutamate, dopamine, serotonin and gamma-aminobutyric acid (GABA). The H3 R positioned at numerous neurotransmission crossroads indicates therapeutic applications of small-molecule H3 R modulators in a number of psychiatric and neurodegenerative diseases with various clinical candidates available. Dual target drugs displaying H3 R antagonism/inverse agonism with inhibition of acetylcholine esterase (AChE), histamine N-methyltransferase (HMT), or serotonin transporter (SERT) are novel class of procognitive agents. Main chemical diversities, pharmacophores, and pharmacological profiles of procognitive agents acting as H3 R antagonists/inverse agonists and dual H3 R antagonists/inverse agonists with inhibiting activity on AChE, HMT, or SERT are highlighted here.

Synthesis, biological evaluation, and computational studies of Tri- and tetracyclic nitrogen-bridgehead compounds as potent dual-acting AChE inhibitors and hH3 receptor antagonists

Combination of AChE inhibiting and histamine H3 receptor antagonizing properties in a single molecule might show synergistic effects to improve cognitive deficits in Alzheimer's disease, since both pharmacological actions are able to enhance cholinergic neurotransmission in the cortex. However, whereas AChE inhibitors prevent hydrolysis of acetylcholine also peripherally, histamine H3 antagonists will raise acetylcholine levels mostly in the brain due to predominant occurrence of the receptor in the central nervous system. In this work, we designed and synthesized two novel classes of tri- and tetracyclic nitrogen-bridgehead compounds acting as dual AChE inhibitors and histamine H3 antagonists by combining the nitrogen-bridgehead moiety of novel AChE inhibitors with a second N-basic fragment based on the piperidinylpropoxy pharmacophore with different spacer lengths. Intensive structure-activity relationships (SARs) with regard to both biological targets led to compound 41 which showed balanced affinities as hAChE inhibitor with IC50 = 33.9 nM, and hH3R antagonism with Ki = 76.2 nM with greater than 200-fold selectivity over the other histamine receptor subtypes. Molecular docking studies were performed to explain the potent AChE inhibition of the target compounds and molecular dynamics studies to explain high affinity at the hH3R.

Molecular dynamics simulations of the adenosine A2a receptor in POPC and POPE lipid bilayers: effects of membrane on protein behavior

Analysis of 300 ns (ns) molecular dynamics (MD) simulations of an adenosine A2a receptor (A2a AR) model, conducted in triplicate, in 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) bilayers reveals significantly different protein dynamical behavior. Principal component analysis (PCA) shows that the dissimilarities stem from interhelical rather than intrahelical motions. The difference in the hydrophobic thicknesses of these simulated lipid bilayers is potentially a significant reason for the observed difference in results. The distinct lipid headgroups might also lead to different molecular interactions and hence different protein loop motions. Overall, the A2a AR shows higher mobility and flexibility in POPC as compared to POPE.

Molecular dynamics simulations of the adenosine A2a receptor: structural stability, sampling, and convergence

Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.

Current progress in Structure-Based Rational Drug Design marks a new mindset in drug discovery

The past decade has witnessed a paradigm shift in preclinical drug discovery with structure-based drug design (SBDD) making a comeback while high-throughput screening (HTS) methods have continued to generate disappointing results. There is a deficit of information between identified hits and the many criteria that must be fulfilled in parallel to convert them into preclinical candidates that have a real chance to become a drug. This gap can be bridged by investigating the interactions between the ligands and their receptors. Accurate calculations of the free energy of binding are still elusive; however progresses were made with respect to how one may deal with the versatile role of water. A corpus of knowledge combining X-ray structures, bioinformatics and molecular modeling techniques now allows drug designers to routinely produce receptor homology models of increasing quality. These models serve as a basis to establish and validate efficient rationales used to tailor and/or screen virtual libraries with enhanced chances of obtaining hits. Many case reports of successful SBDD show how synergy can be gained from the combined use of several techniques. The role of SBDD with respect to two different classes of widely investigated pharmaceutical targets: (a) protein kinases (PK) and (b) G-protein coupled receptors (GPCR) is discussed. Throughout these examples prototypical situations covering the current possibilities and limitations of SBDD are presented.

Drug design for ever, from hype to hope

In its first 25 years JCAMD has been disseminating a large number of techniques aimed at finding better medicines faster. These include genetic algorithms, COMFA, QSAR, structure based techniques, homology modelling, high throughput screening, combichem, and dozens more that were a hype in their time and that now are just a useful addition to the drug-designers toolbox. Despite massive efforts throughout academic and industrial drug design research departments, the number of FDA-approved new molecular entities per year stagnates, and the pharmaceutical industry is reorganising accordingly. The recent spate of industrial consolidations and the concomitant move towards outsourcing of research activities requires better integration of all activities along the chain from bench to bedside. The next 25 years will undoubtedly show a series of translational science activities that are aimed at a better communication between all parties involved, from quantum chemistry to bedside and from academia to industry. This will above all include understanding the underlying biological problem and optimal use of all available data.

Molecular determinants of ligand binding modes in the histamine H(4) receptor: linking ligand-based three-dimensional quantitative structure-activity relationship (3D-QSAR) models to in silico guided receptor mutagenesis studies

The histamine H(4) receptor (H(4)R) is a G protein-coupled receptor (GPCR) that plays an important role in inflammation. Similar to the homologous histamine H(3) receptor (H(3)R), two acidic residues in the H(4)R binding pocket, D(3.32) and E(5.46), act as essential hydrogen bond acceptors of positively ionizable hydrogen bond donors in H(4)R ligands. Given the symmetric distribution of these complementary pharmacophore features in H(4)R and its ligands, different alternative ligand binding mode hypotheses have been proposed. The current study focuses on the elucidation of the molecular determinants of H(4)R-ligand binding modes by combining (3D) quantitative structure-activity relationship (QSAR), protein homology modeling, molecular dynamics simulations, and site-directed mutagenesis studies. We have designed and synthesized a series of clobenpropit (N-(4-chlorobenzyl)-S-[3-(4(5)-imidazolyl)propyl]isothiourea) derivatives to investigate H(4)R-ligand interactions and ligand binding orientations. Interestingly, our studies indicate that clobenpropit (2) itself can bind to H(4)R in two distinct binding modes, while the addition of a cyclohexyl group to the clobenpropit isothiourea moiety allows VUF5228 (5) to adopt only one specific binding mode in the H(4)R binding pocket. Our ligand-steered, experimentally supported protein modeling method gives new insights into ligand recognition by H(4)R and can be used as a general approach to elucidate the structure of protein-ligand complexes.

The histamine H₄ receptor: targeting inflammatory disorders

The discovery of the histamine H(4) receptor has added a new chapter to the century of extensive biogenic amine research. The human histamine H(4) receptor is mainly expressed in cells of the human immune system (e.g. mast cells, eosinophils, monocytes, dendritic cells, T cells) and mediates several effects on chemotaxis with numerous cell types. The distinct expression pattern and the immunomodulatory role highlight its physiological relevance in inflammatory and immunological processes. Inflammatory conditions, e.g. allergy, asthma and autoimmune diseases, were for a long time thought to be mainly mediated by activation of the histamine H(1) receptor subtype. However, in the treatment of diseases as chronic pruritus, asthma and allergic rhinitis the use of histamine H(1) receptor antagonists is unsatisfying. Selective H(4) receptor ligands and/or synergism of histamine H(1) and H(4) receptor modulation may be more effective in such pathophysiological conditions. Promising preclinical studies underline its role as an attractive target in the treatment of inflammatory and autoimmune disorders. Meanwhile, first histamine H(4) receptor antagonist has reached clinical phases for the treatment of respiratory diseases.

Triazole ligands reveal distinct molecular features that induce histamine H4 receptor affinity and subtly govern H4/H3 subtype selectivity

The histamine H(3) (H(3)R) and H(4) (H(4)R) receptors attract considerable interest from the medicinal chemistry community. Given their relatively high homology yet widely differing therapeutic promises, ligand selectivity for the two receptors is crucial. We interrogated H(4)R/H(3)R selectivities using ligands with a [1,2,3]triazole core. Cu(I)-assisted "click chemistry" was used to assemble diverse [1,2,3]triazole compounds (6a-w and 7a-f), many containing a peripheral imidazole group. The imidazole ring posed some problems in the click chemistry putatively due to Cu(II) coordination, but Boc protection of the imidazole and removal of oxygen from the reaction mixture provided effective strategies. Pharmacological studies revealed two monosubstituted imidazoles (6h,p) with <10 nM H(4)R affinities and >10-fold H(4)R/H(3)R selectivity. Both compounds possess a cycloalkylmethyl group and appear to target a lipophilic pocket in H(4)R with high steric precision. The use of the [1,2,3]triazole scaffold is further demonstrated by the notion that simple changes in spacer length or peripheral groups can reverse the selectivity toward H(3)R. Computational evidence is provided to account for two key selectivity switches and to pinpoint a lipophilic pocket as an important handle for H(4)R over H(3)R selectivity.

Kojic acid derivatives as histamine H(3) receptor ligands

The histamine H(3) receptor (H(3)R) is a promising target in the development of new compounds for the treatment of mainly centrally occurring diseases. However, emerging novel therapeutic concepts have been introduced and some indications in the H(3)R field, e.g. migraine, pain or allergic rhinitis, might take advantage of peripherally acting ligands. In this work, kojic acid-derived structural elements were inserted into a well established H(3)R antagonist/inverse agonist scaffold to investigate the bioisosteric potential of γ-pyranones with respect to the different moieties of the H(3)R pharmacophore. The most affine compounds showed receptor binding in the low nanomolar concentration range. Evaluation and comparison of kojic acid-containing ligands and their corresponding phenyl analogues (3-7) revealed that the newly integrated scaffold greatly influences chemical properties (S Log P, topological polar surface area (tPSA)) and hence, potentially modifies the pharmacokinetic profile of the different derivatives. Benzyl-1-(4-(3-(piperidin-1-yl)propoxy)phenyl)methanamine ligands 3 and 4 belong to the centrally acting diamine-based class of H(3)R antagonist/inverse agonist, whereas kojic acid analogues 6 and 7 might act peripherally. The latter compounds state promising lead structures in the development of H(3)R ligands with a modified profile of action.

Comparison of the pharmacological properties of human and rat histamine H(3)-receptors

Ligand pharmacology of histamine H(3)-receptors is species-dependent. In previous studies, two amino acids in transmembrane domain 3 (TM III) were shown to play a significant role. In this study, we characterized human and rat histamine H(3)-receptors (hH(3)R and rH(3)R, respectively), co-expressed with mammalian G proteins in Sf9 insect cell membranes. We compared a series of imidazole-containing H(3)R ligands in radioligand binding and steady-state GTPase assays. H(3)Rs similarly coupled to Gα(i/o)-proteins. Affinities and potencies of the agonists histamine, N(α)-methylhistamine and R-(α)-methylhistamine were in the same range. Imetit was only a partial agonist. The pharmacology of imetit and proxifan was similar at both species. However, impentamine was more potent and efficacious at rH(3)R. The inverse agonists ciproxifan and thioperamide showed higher potency but lower efficacy at rH(3)R. Clobenpropit was not species-selective. Strikingly, imoproxifan was almost full agonist at hH(3)R, but an inverse agonist at rH(3)R. Imoproxifan was docked into the binding pocket of inactive and active hH(3)R- and rH(3)R-models and molecular dynamic simulations were performed. Imoproxifan bound to hH(3)R and rH(3)R in E-configuration, which represents the trans-isomer of the oxime-moiety as determined in crystallization studies, and stabilized active hH(3)R-, but inactive rH(3)R-conformations. Large differences in electrostatic surfaces between TM III and TM V cause differential orientation of the oxime-moiety of imoproxifan, which then differently interacts with the rotamer toggle switch Trp(6.48) in TM VI. Collectively, the substantial species differences at H(3)Rs are explained at a molecular level by the use of novel H(3)R active-state models.

Homology modelling of the human adenosine A2B receptor based on X-ray structures of bovine rhodopsin, the beta2-adrenergic receptor and the human adenosine A2A receptor

A three-dimensional model of the human adenosine A2B receptor was generated by means of homology modelling, using the crystal structures of bovine rhodopsin, the beta2-adrenergic receptor, and the human adenosine A2A receptor as templates. In order to compare the three resulting models, the binding modes of the adenosine A2B receptor antagonists theophylline, ZM241385, MRS1706, and PSB601 were investigated. The A2A-based model was much better able to stabilize the ligands in the binding site than the other models reflecting the high degree of similarity between A2A and A2B receptors: while the A2B receptor shares about 21% of the residues with rhodopsin, and 31% with the beta2-adrenergic receptor, it is 56% identical to the adenosine A2A receptor. The A2A-based model was used for further studies. The model included the transmembrane domains, the extracellular and the intracellular hydrophilic loops as well as the terminal domains. In order to validate the usefulness of this model, a docking analysis of several selective and nonselective agonists and antagonists was carried out including a study of binding affinities and selectivities of these ligands with respect to the adenosine A2A and A2B receptors. A common binding site is proposed for antagonists and agonists based on homology modelling combined with site-directed mutagenesis and a comparison between experimental and calculated affinity data. The new, validated A2B receptor model may serve as a basis for developing more potent and selective drugs.

Acidic elements in histamine H(3) receptor antagonists

Antagonists of the human histamine H(3) receptor (hH(3)R) often contain a second basic moiety, which is well known to boost affinity on this histamine receptor subtype. Here, we prepared compounds with acidic moieties of different pK(a) values to figure out that the hH(3)R tolerates these functionalities when added to a common pharmacophore blueprint. Depending on the acidic, electronic and steric features the designed ligands showed hH(3)R affinities in the nanomolar concentration range. Additionally, selected ligands were tested but failed as dual acting hH(3)R/hPPAR (human peroxisome proliferator-activated receptor) ligands.

Physical binding pocket induction for affinity prediction

Computational methods for predicting ligand affinity where no protein structure is known generally take the form of regression analysis based on molecular features that have only a tangential relationship to a protein/ligand binding event. Such methods have limited utility when structural variation moves beyond congeneric series. We present a novel approach based on the multiple-instance learning method of Compass, where a physical model of a binding site is induced from ligands and their corresponding activity data. The model consists of molecular fragments that can account for multiple positions of literal protein residues. We demonstrate the method on 5HT1a ligands by training on a series with limited scaffold variation and testing on numerous ligands with variant scaffolds. Predictive error was between 0.5 and 1.0 log units (0.7-1.4 kcal/mol), with statistically significant rank correlations. Accurate activity predictions of novel ligands were demonstrated using a validation approach where a small number of ligands of limited structural variation known at a fixed time point were used to make predictions on a blind test set of widely varying molecules, some discovered at a much later time point.

Docking and chemoinformatic screens for new ligands and targets

Computer-based docking screens are now widely used to discover new ligands for targets of known structure; in the last two years alone, the discovery of ligands for more than 20 proteins has been reported. Recently, investigators have also turned to predicting new substrates for enzymes of unknown function, taking docking in a wholly new direction. Increasingly, the hit rates, the true-positives, and the false-positives from the docking screens are being compared to those from empirical, high-throughput screens, revealing the strengths, weaknesses, and complementarities of both techniques. The recent efflorescence of GPCR structures has made these quintessential drug targets available to structure-based approaches. Consistent with their 'druggability', the docking screens have returned high hit rates and potent molecules. Finally, in the last several years, an approach almost exactly opposite to docking has also appeared; this pharmacological network approach begins not with the structure of the target but rather those of drug molecules and asks, given a pattern of chemistry in the ligands, what targets may a particular drug bind to? This method, which returns to an older, pharmacology logic, has been surprisingly successful in predicting new 'off-targets' for established drugs.

New vistas in GPCR 3D structure prediction

Human G-protein coupled receptors (hGPCRs) comprise the most prominent family of validated drug targets. More than 50% of approved drugs reveal their therapeutic effects by targeting this family. Accurate models would greatly facilitate the process of drug discovery and development. However, 3-D structure prediction of GPCRs remains a challenge due to limited availability of resolved structure. The X-ray structures have been solved for only four such proteins. The identity between hGPCRs and the potential templates is mostly less than 30%, well below the level at which sequence alignment can be done regularly. In this study, we analyze a large database of human G-protein coupled receptors that are members of family A in order to optimize usage of the available crystal structures for molecular modeling of hGPCRs. On the basis of our findings in this study, we propose to regard specific parts from the trans-membrane domains of the reference receptor helices as appropriate template for constructing models of other GPCRs, while other residues require other techniques for their remodeling and refinement. The proposed hypothesis in the current study has been tested by modeling human beta2-adrenergic receptor based on crystal structures of bovine rhodopsin (1F88) and human A2A adenosine receptor (3EML). The results have shown some improvement in the quality of the predicted models compared to Modeller software.