Studies of H4R antagonists using 3D-QSAR, molecular docking and molecular dynamics

Three-dimensional quantitative structure–activity relationship studies were performed on a series of 88 histamine receptor 4 (H4R) antagonists in an attempt to elucidate the 3D structural features required for activity. Several in silico modeling approaches, including comparative molecular field analysis (CoMFA), comparative similarity indices analysis (CoMSIA), molecular docking, and molecular dynamics (MD), were carried out. The results show that both the ligand-based CoMFA model (Q2 = 0.548, Rncv2 = 0.870, Rpre2 = 0.879, SEE = 0.410, SEP = 0.386) and the CoMSIA model (Q2 = 0.526, Rncv2 =0.866, Rpre2 = 0.848, SEE = 0.416, SEP = 0.413) are acceptable, as they show good predictive capabilities. Furthermore, a combined analysis incorporating CoMFA, CoMSIA contour maps and MD results shows that (1) compounds with bulky or hydrophobic substituents at positions 4–6 in ring A (R2 substituent), positively charged or hydrogen-bonding (HB) donor groups in the R1 substituent, and hydrophilic or HB acceptor groups in ring C show enhanced biological activities, and (2) the key amino acids in the binding pocket are TRP67, LEU71, ASP94, TYR95, PHE263 and GLN266. To our best knowledge, this work is the first to report the 3D-QSAR modeling of these H4R antagonists. The conclusions of this work may lead to a better understanding of the mechanism of antagonism and aid in the design of new, more potent H4R antagonists.