In silico studies of diarylpyridine derivatives as novel HIV-1 NNRTIs using docking-based 3D-QSAR, molecular dynamics, and pharmacophore modeling approaches

A series of novel diarylpyridine derivatives has recently been identified as HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs), and most of them exhibited potent activities against HIV-1 strains, with EC50 values in the low nanomolar range. However, the three-dimensional quantitative structure-activity relationships (3D-QSARs) and pharmacophore characteristics of these compounds remain to be studied. In the present study, forty-two diarylpyridine derivatives were firstly docked into HIV-1 reverse transcriptase, and molecular dynamics (10 ns) simulations were further performed to validate the reliability of the docking results, which indicated that residues Lys101, Tyr181, Tyr188, Phe227, and Trp229 might play important roles in binding with these diarylpyridines. The U-shaped docking conformations of all compounds were then used to construct 3D-QSAR and pharmacophore models. The satisfactory statistical parameters of CoMFA (q(loo)(2) = 0.665, r(ncv)(2) = 0.989, r(pred)(2) = 0.962, etc.) and CoMSIA (q(loo)(2) = 0.727, r(ncv)(2) = 0.988, r(pred)(2) = 0.912, etc.) models demonstrated that both constructed models had excellent predictability, and their contour maps gave insights into the structural requirements of the diarylpyridines for the anti-HIV-1 activity. A docking-conformation-based pharmacophore model, containing three hydrophobic centers, three hydrogen-bond acceptors, and three hydrogen-bond donors, was also established. The observations in this study might provide important information for the rational design and development of novel HIV-1 NNRTIs.