The histamine H₄ receptor, belonging to the family of G-protein coupled receptors, is an increasingly attractive drug target. It plays an indispensable role in many cellular pathways, and numerous H₄R ligands are being studied for the treatment of several inflammatory, allergic, and autoimmune disorders, including pulmonary fibrosis. Activation of H₄R is involved in cytokine production and mediates mast cell activation and eosinophil chemotaxis. The importance of this receptor has also been shown in inflammatory models: peritonitis, respiratory tract inflammation, colitis, osteoarthritis, and rheumatoid arthritis. Recent studies suggest that H₄R acts as a modulator in cancer, neuropathic pain, vestibular disorders, and type-2 diabetes, however, its role is still not fully understood.

The recent developments of fast reliable docking, virtual screening and other algorithms gave rise to discovery of many novel ligands of histamine receptors that could be used for treatment of allergic inflammatory disorders, central nervous system pathologies, pain, cancer and obesity. Furthermore, the pharmacological profiles of ligands clearly indicate that these receptors may be considered as targets not only for selective but also for multi-target drugs that could be used for treatment of complex disorders such as Alzheimer’s disease. Therefore, analysis of protein-ligand recognition in the binding site of histamine receptors and also other molecular targets has become a valuable tool in drug design toolkit. This review covers the period 2014–2020 in the field of theoretical investigations of histamine receptors mostly based on molecular modeling as well as the experimental characterization of novel ligands of these receptors.

Structure-based screening approach targeting mGlu2 receptor was carried out to identify good chemical starting points for anti-epileptic therapy. Interactive modes of final 12 compounds identified on the basis of screening of Asinex library, binding energy analysis, ADME profiling with special emphasis for CNS ranges, and toxicity analysis were studied and showed good binding modes in the mGluR2-active site. Enrichment studies for validating screening protocol were carried out which gave ROC values 0.98 (AUC = 0.96) for SP, 0.97 (AUC = 0.95) for XP with BEDROC analysis. Our results indicate that all the 12 hits showed good CNS drug-like properties, have better binding free energy and ADME profile as compared to co-crystallized ligand with the best ligand hit retaining conserved hydrogen bond interactions with Ala-166, Thr-168, Ser-145, and Arg-61 residues in bilobatevenus fly-trap domain of mGluR2 receptor. Molecular dynamics simulations proved that the two potential hits, qualifying all screening parameters, are stable in the receptor active site pocket, confirming the potential of the identified hits as a specific target for mGluR2. Because the newly discovered mGluR2 agonists are structurally different with T_c values ranging from 0.57 to 0.92, all of them can be considered for further de novo design methods.

Biological mechanism attributing mutations in KCNQ2/Q3 results in benign familial neonatal epilepsy (BFNE), a rare form of epilepsy and thus neglected. It offers a potential target for antiepileptic drug discovery. In the present work, a pharmacophore-based 3D-QSAR model was generated for a series of N-pyridyl and pyrimidine benzamides possessing KCNQ2/Q3 opening activity. The pharmacophore model generated contains one hydrogen bond donor (D), one hydrophobic (H), and two aromatic rings (R). They are the crucial molecular write-up detailing predicted binding efficacy of high affinity and low affinity ligands for KCNQ2/Q3 opening activity. Furthermore, it has been validated by using a biological correlation between pharmacophore hypothesis-based 3D-QSAR variables and functional fingerprints of openers responsible for the receptor binding and also by docking of these benzamides into the validated homology model. Excellent statistical computational tools of QSAR model such as good correlation coefficient (R²>0.80), higher F value (F > 39), and excellent predictive power (Q² > 0.7) with low standard deviation (SD <0.3) strongly suggest that the developed model could be used for prediction of antiepileptic activity of newer analogs. A preliminary pharmacokinetic profile of these derivatives was also performed on the basis of QikProp predictions.

Alzheimer’s disease (AD), a progressive neurodegenerative disorder is the most common cause of dementia among elderly people. To date, the successful therapeutic strategy to treat AD is maintaining the levels of acetylcholine via inhibiting acetylcholinesterase (AChE). The present study involves identification of newer AChE inhibitors by dual approach of e-pharmacophore and structure-based virtual screening of Asinex library. Robustness of docking protocol was validated by enrichment calculation with ROC value .71 and BEDROC value .028. Among 11 selected hits, ZINC72338524 with best MM-GBSA dG binding shows optimal range of CNS properties and ligand–AChE complex stability. Further, molecular dynamics study revealed its molecular interactions with Trp86, Phe338, and Tyr341 amino acid residues of catalytic anionic site and Tyr124, Ser125, and Trp286 amino acid residues of peripheral anionic site. Physicochemical properties and ADMET risk prediction indicates their potential in druggability and safety.