Docking and MD study of histamine H4R based on the crystal structure of H1R

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.

Evolutionary Couplings and Molecular Dynamic Simulations Highlight Details of GPCRs Heterodimers' Interfaces

A growing body of evidence suggests that only a few amino acids ("hot-spots") at the interface contribute most of the binding energy in transient protein-protein interactions. However, experimental protocols to identify these hot-spots are highly labor-intensive and expensive. Computational methods, including evolutionary couplings, have been proposed to predict the hot-spots, but they generally fail to provide details of the interacting amino acids. Here we showed that unbiased evolutionary methods followed by biased molecular dynamic simulations could achieve this goal and reveal critical elements of protein complexes. We applied the methodology to selected G-protein coupled receptors (GPCRs), known for their therapeutic properties. We used the structure-prior-assisted direct coupling analysis (SP-DCA) to predict the binding interfaces of A2aR/D2R, CB1R/D2R, A2aR/CB1R, 5HT2AR/D2R, and 5-HT2AR/mGluR2 receptor heterodimers, which all agreed with published data. In order to highlight details of the interactions, we performed molecular dynamic (MD) simulations using the newly developed AWSEM energy model. We found that these receptors interact primarily through critical residues at the C and N terminal domains and the third intracellular loop (ICL3). The MD simulations showed that these residues are energetically necessary for dimerization and revealed their native conformational state. We subsequently applied the methodology to the 5-HT2AR/5-HTR4R heterodimer, given its implication in drug addiction and neurodegenerative pathologies such as Alzheimer's disease (AD). Further, the SP-DCA analysis showed that 5-HT2AR and 5-HTR4R heterodimerize through the C-terminal domain of 5-HT2AR and ICL3 of 5-HT4R. However, elucidating the details of GPCR interactions would accelerate the discovery of druggable sites and improve our knowledge of the etiology of common diseases, including AD.

Effect of Ions and Sequence Variants on the Antagonist Binding Properties of the Histamine H1 Receptor

The histamine H₁ receptor (H₁R) is a G protein-coupled receptor (GPCR) and represents a main target in the treatment of allergic reactions as well as inflammatory reactions and depressions. Although the overall effect of antagonists on H₁ function has been extensively investigated, rather little is known about the potential modulatory effect of ions or sequence variants on antagonist binding. We investigated the dynamics of a phosphate ion present in the crystal structure and of a sodium ion, for which we determined the position in the allosteric pocket by metadynamics simulations. Both types of ions exhibit significant dynamics within their binding site; however, some key contacts remain stable over the simulation time, which might be exploited to develop more potent drugs targeting these sites. The dynamics of the ions is almost unaffected by the presence or absence of doxepin, as also reflected in their small effect (less than 1 kcal·mol^-1) on doxepin binding affinity. We also examined the effect of four H₁R sequence variants observed in the human population on doxepin binding. These variants cause a reduction in doxepin affinity of up to 2.5 kcal·mol^-1, indicating that personalized medical treatments that take into account individual mutation patterns could increase precision in the dosage of GPCR-targeting drugs.

Use multiscale simulation to explore the effects of the homodimerizations between different conformation states on the activation and allosteric pathway for the μ-opioid receptor

Recently, oligomers of G-protein coupled receptors (GPCRs) have been an important topic in the GPCR fields. However, knowledge about their structures and activation mechanisms is very limited due to the absence of crystal structures reported. In this work, we used multiscale simulations to study the effects of homodimerization between different conformation states on their activation, dynamic behaviors, and allosteric communication pathways for μ-OR. The results indicated that the dimerization of one inactive monomer with either one inactive monomer or one active one could enhance its constitutive activation. However, the conformation state of the other protomer (e.g., active or inactive) can influence the activated extent. The dimerization between the two inactive protomers leads to a negative cooperativity for their activation, which should contribute to the asymmetric activation of GPCR dimers observed in some experiments. On the other hand, for the active monomer, its dimerization with one inactive receptor could alleviate its deactivation, whereby negative and positive cooperativities can be observed between the two subunits of the dimer, depending on the different regions. Observations from protein structure network (PSN) analysis indicated that the dimerization of one inactive monomer with one active one would cause a significant drop in the number of main pathways from the ligand binding pocket to the G-protein coupled region for the inactive protomer, while the impact is minor for the active protomer. But, for the active monomer or the inactive one, its dimerization with one inactive monomer would significantly change the types of residues participating in the pathway with the highest frequency.

Design and synthesis of histamine H3/H4 receptor ligands with a cyclopropane scaffold

We previously designed and synthesized a series of histamine analogues with an imidazolylcyclopropane scaffold and identified potent non-selective antagonists for histamine H₃ and H₄ receptor subtypes. In this study, to develop H₄ selective ligands, we newly designed and synthesized cyclopropane-based derivatives having an indole, benzimidazole, or piperazine structure, which are components of representative H₄ selective antagonists such as JNJ7777120 and JNJ10191584. Among the synthesized derivatives, imidazolylcyclopropanes 12 and 13 conjugated with a benzimidazole showed binding affinity to the H₃ and H₄ receptors comparable to that of a well-known non-selective H₃/H₄ antagonist, thioperamide. These results suggest that the binding modes of the cyclopropane-based H₃/H₄ ligands in the H₄ receptor can be different from those of the indole/benzimidazole-piperazine derivatives.

Discovery of nanomolar ligands with novel scaffolds for the histamine H4 receptor by virtual screening

The involvement of histamine H4 receptor (H4R) in immune cells chemotaxis and mediator release makes it an attractive target for the treatment of inflammation disorders. A decade of medicinal chemistry efforts has led to several promising ligands, although the chemical structures described so far possesses a singular limited diversity. We report here the discovery of novel structures, belonging to completely different scaffolds. The virtual screening was planed as a two-steps process. First, using a "scout screening" methodology, we have experimentally probed the H4R ligand binding site using a small size chemical library with very diverse structures, and identified a hit that further assist us in refining a raw 3D homology model. Second, the refined 3D model was used to conduct a widened virtual screening. This two-steps strategy proved to be very successful, both in terms of structural diversity and hit rate (23%). Moreover, the hits have high affinity for the H4R, with most potent ligands in the nanomolar range.

Computational study on the interaction between CCR5 and HIV-1 entry inhibitor maraviroc: insight from accelerated molecular dynamics simulation and free energy calculation

C-C chemokine receptor type 5 (CCR5) is the co-receptor of human immunodeficiency virus type 1 (HIV-1) and plays an important role in HIV-1 virus infection. Maraviroc has been proved to be effective for anti-HIV-1 by targeting CCR5. Understanding the detailed interaction mechanism between CCR5 and Maraviroc will be of great help to the rational design of a more potential inverse agonist to block HIV-1 infection. Here, we performed molecular dynamics (MD) simulation and accelerated MD simulation (aMD) to study the interaction mechanism between CCR5 and Maraviroc based on a recently reported crystal structure. The results of MD simulation demonstrate that Maraviroc can form stable hydrogen bonds with residues Tyr37(1.39), Tyr251(6.51) and Glu283(7.39). The results of aMD simulation indicate that the carboxamide moiety is more flexible than the tropane group of Maraviroc in the pocket of CCR5. The electrostatic potential analysis proves that Maraviroc can escape from the pocket of CCR5 along the negative electrostatic potential pathway during the dissociation process. The free energy calculation illustrates that there exist three binding pockets during the dissociation process of Maraviroc. Our results will be useful for understanding the interaction mechanism between CCR5 and Maraviroc as well as for the rational design of a more potent inverse agonist.

Computational Advances for the Development of Allosteric Modulators and Bitopic Ligands in G Protein-Coupled Receptors

Allosteric modulators of G protein-coupled receptors (GPCRs), which target at allosteric sites, have significant advantages against the corresponding orthosteric compounds including higher selectivity, improved chemical tractability or physicochemical properties, and reduced risk of receptor oversensitization. Bitopic ligands of GPCRs target both orthosteric and allosteric sites. Bitopic ligands can improve binding affinity, enhance subtype selectivity, stabilize receptors, and reduce side effects. Discovering allosteric modulators or bitopic ligands for GPCRs has become an emerging research area, in which the design of allosteric modulators is a key step in the detection of bitopic ligands. Radioligand binding and functional assays ([(35)S]GTPγS and ERK1/2 phosphorylation) are used to test the effects for potential modulators or bitopic ligands. High-throughput screening (HTS) in combination with disulfide trapping and fragment-based screening are used to aid the discovery of the allosteric modulators or bitopic ligands of GPCRs. When used alone, these methods are costly and can often result in too many potential drug targets, including false positives. Alternatively, low-cost and efficient computational approaches are useful in drug discovery of novel allosteric modulators and bitopic ligands to help refine the number of targets and reduce the false-positive rates. This review summarizes the state-of-the-art computational methods for the discovery of modulators and bitopic ligands. The challenges and opportunities for future drug discovery are also discussed.

Allosteric Binding Site and Activation Mechanism of Class C G-Protein Coupled Receptors: Metabotropic Glutamate Receptor Family

Metabotropic glutamate receptors (mGluR) are mainly expressed in the central nervous system (CNS) and contain eight receptor subtypes, named mGluR1 to mGluR8. The crystal structures of mGluR1 and mGluR5 that are bound with the negative allosteric modulator (NAM) were reported recently. These structures provide a basic model for all class C of G-protein coupled receptors (GPCRs) and may aid in the design of new allosteric modulators for the treatment of CNS disorders. However, these structures are only combined with NAMs in the previous reports. The conformations that are bound with positive allosteric modulator (PAM) or agonist of mGluR1/5 remain unknown. Moreover, the structural information of the other six mGluRs and the comparisons of the mGluRs family have not been explored in terms of their binding pockets, the binding modes of different compounds, and important binding residues. With these crystal structures as the starting point, we built 3D structural models for six mGluRs by using homology modeling and molecular dynamics (MD) simulations. We systematically compared their allosteric binding sites/pockets, the important residues, and the selective residues by using a series of comparable dockings with both the NAM and the PAM. Our results show that several residues played important roles for the receptors' selectivity. The observations of detailed interactions between compounds and their correspondent receptors are congruent with the specificity and potency of derivatives or compounds bioassayed in vitro. We then carried out 100 ns MD simulations of mGluR5 (residue 26-832, formed by Venus Flytrap domain, a so-called cysteine-rich domain, and 7 trans-membrane domains) bound with antagonist/NAM and with agonist/PAM. Our results show that both the NAM and the PAM seemed stable in class C GPCRs during the MD. However, the movements of "ionic lock," of trans-membrane domains, and of some activation-related residues in 7 trans-membrane domains of mGluR5 were congruent with the findings in class A GPCRs. Finally, we selected nine representative bound structures to perform 30 ns MD simulations for validating the stabilities of interactions, respectively. All these bound structures kept stable during the MD simulations, indicating that the binding poses in this present work are reasonable. We provided new insight into better understanding of the structural and functional roles of the mGluRs family and facilitated the future structure-based design of novel ligands of mGluRs family with therapeutic potential.

Histamine H4 receptor ligands: future applications and state of art

Histamine is a chemical transmitter found practically in whole organism and exerts its effects through the interaction with H1 to H4 histaminergic receptors. Specifically, H4 receptors are found mainly in immune cells and blood-forming tissues, thus are involved in inflammatory and immune processes, as well as some actions in central nervous system. Therefore, H4 receptor ligands can have applications in the treatment of chronic inflammatory and immune diseases and may be novel therapeutic option in these conditions. Several H4 receptor ligands have been described from early 2000's until nowadays, being imidazole, indolecarboxamide, 2-aminopyrimidine, quinazoline, and quinoxaline scaffolds the most explored and discussed in this review. Moreover, several studies of molecular modeling using homology models of H4 receptor and QSAR data of the ligands are summarized. The increasing and promising therapeutic applications are leading these compounds to clinical trials, which probably will be part of the next generation of blockbuster drugs.

Sequential application of ligand and structure based modeling approaches to index chemicals for their hH4R antagonism

The human histamine H4 receptor (hH4R), a member of the G-protein coupled receptors (GPCR) family, is an increasingly attractive drug target. It plays a key role in many cell pathways and many hH4R ligands are studied for the treatment of several inflammatory, allergic and autoimmune disorders, as well as for analgesic activity. Due to the challenging difficulties in the experimental elucidation of hH4R structure, virtual screening campaigns are normally run on homology based models. However, a wealth of information about the chemical properties of GPCR ligands has also accumulated over the last few years and an appropriate combination of these ligand-based knowledge with structure-based molecular modeling studies emerges as a promising strategy for computer-assisted drug design. Here, two chemoinformatics techniques, the Intelligent Learning Engine (ILE) and Iterative Stochastic Elimination (ISE) approach, were used to index chemicals for their hH4R bioactivity. An application of the prediction model on external test set composed of more than 160 hH4R antagonists picked from the chEMBL database gave enrichment factor of 16.4. A virtual high throughput screening on ZINC database was carried out, picking ∼ 4000 chemicals highly indexed as H4R antagonists' candidates. Next, a series of 3D models of hH4R were generated by molecular modeling and molecular dynamics simulations performed in fully atomistic lipid membranes. The efficacy of the hH4R 3D models in discrimination between actives and non-actives were checked and the 3D model with the best performance was chosen for further docking studies performed on the focused library. The output of these docking studies was a consensus library of 11 highly active scored drug candidates. Our findings suggest that a sequential combination of ligand-based chemoinformatics approaches with structure-based ones has the potential to improve the success rate in discovering new biologically active GPCR drugs and increase the enrichment factors in a synergistic manner.

Histamine H4 receptor in oral lichen planus

OBJECTIVES

Oral lichen planus (OLP) is an autoimmune disease characterized by a band-like T-cell infiltrate below the apoptotic epithelial cells and degenerated basement membrane. We tested the hypothesis that the high-affinity histamine H4 receptors (H4 Rs) are downregulated in OLP by high histamine concentrations and proinflammatory T-cell cytokines.

MATERIALS AND METHODS

Immunohistochemistry and immunofluorescence staining, image analysis and quantitative real-time polymerase chain reaction of tissue samples and cytokine-stimulated cultured SCC-25 and primary human oral keratinocytes.

RESULTS

H4 R immunoreactivity was weak in OLP and characterized by mast cell (MC) hyperplasia and degranulation. In contrast to controls, H4 R immunostaining and MC counts were negatively correlated in OLP (P = 0.003). H4 R agonist at nanomolar levels led to a rapid internalization of H4 Rs, whereas high histamine concentration and interferon-γ decreased HRH4 -gene transcripts.

CONCLUSION

Healthy oral epithelial cells are equipped with H4 R, which displays a uniform staining pattern in a MC-independent fashion. In contrast, in OLP, increased numbers of activated MCs associate with increasing loss of epithelial H4 R. Cell culture experiments suggest a rapid H4 R stimulation-dependent receptor internalization and a slow cytokine-driven decrease in H4 R synthesis. H4 R may be involved in the maintenance of healthy oral mucosa. In OLP, this maintenance might be impaired by MC degranulation and inflammatory cytokines.