In silico study of 3-hydroxypyrimidine-2,4-diones as inhibitors of HIV RT-associated RNase H using molecular docking, molecular dynamics, 3D-QSAR, and pharmacophore models

• Published in: New journal of chemistry Vol. 43; no. 43; pp. 17004 - 17017
• Main Authors: Gao, Ya, Chen, Yanming, Tian, Yafeng, Zhao, Yilan, Wu, Fengshou, Luo, Xiaogang, Ju, Xiulian, Liu, Genyan
• Format: Journal Article
• Language: English
• Published: Cambridge Royal Society of Chemistry 2019
• Subjects: Computer simulation; Dynamic structural analysis; Inhibitors; Molecular docking; Molecular dynamics; Molecular structure; Multiplication; Pharmacology; Three dimensional models
• ISSN: 1144-0546; 1369-9261
• DOI: 10.1039/C9NJ03353J

HIV reverse transcriptase (RT)-associated ribonuclease H (RNase H) plays an important role in HIV multiplication and represents a novel anti-HIV target. Recently, a novel series of 3-hydroxypyrimidine-2,4-dione ( HPD ) derivatives have been reported as potent inhibitors of HIV RT-associated RNase H; they also exhibit antiviral activities, with EC 50 values in the low micromolar range. To better understand their structure–activity relationships and mechanisms of action, an integrated computational study, including molecular docking, molecular dynamics (MD), three-dimensional quantitative structure–activity relationship (3D-QSAR), and pharmacophore modeling, was performed on these HPD s. Ninety-three HPD s were firstly docked into the RNase H active site; then, MD simulations were performed to validate the accuracy of the docking results. The comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods were used to generate 3D-QSAR models. Structure-based CoMFA ( q 2 = 0.908, R 2 = 0.978, r pred 2 = 0.949) and CoMSIA ( q 2 = 0.908, R 2 = 0.960, r pred 2 = 0.919) models and docking-based CoMFA ( q 2 = 0.832, R 2 = 0.985, r pred 2 = 0.967) and CoMSIA ( q 2 = 0.927, R 2 = 0.990, r pred 2 = 0.977) models were constructed and exhibited excellent predictive ability. The generated pharmacophore model provided deep insight into the pharmacological structural characteristics of the HPD s. Nine virtually screened compounds and six newly designed compounds based on the pharmacophore and 3D-QSAR models are potential leads of HIV RNase H inhibitors. These results may provide important information for the design and development of potent and novel HIV RNase H inhibitors.