汪　瀅，唐國榮，劉澍楠

(中國人民解放軍第476醫(yī)院，福州 350002)

結(jié)合分子相似性、藥效團(tuán)和分子對(duì)接篩選新的HIV-1蛋白酶抑制劑

汪瀅，唐國榮，劉澍楠*

(中國人民解放軍第476醫(yī)院，福州 350002)

摘要：結(jié)合分子相似性、藥效團(tuán)和分子對(duì)接建立兼顧計(jì)算效率和預(yù)測(cè)準(zhǔn)確度的HIV-1蛋白酶抑制劑篩選方法。首先通過對(duì)現(xiàn)有HIV-1蛋白酶抑制劑分子進(jìn)行相似性分析，選取代表性的HIV-1蛋白酶抑制劑作為模板分子，構(gòu)建和優(yōu)化藥效團(tuán)模型，并從1萬個(gè)化合物中優(yōu)先篩選出500個(gè)化合物。而后采用分子對(duì)接方法進(jìn)一步考察化合物與HIV-1 蛋白酶結(jié)合情況，得到 4個(gè)新的活性候選化合物，并進(jìn)行其結(jié)合自由能計(jì)算和抗突變性分析。結(jié)果表明新候選化合物ST025723和HIV-1蛋白酶表現(xiàn)出較好的相互作用和抗突變性，具有深入研究的價(jià)值，同時(shí)也證明分子相似性、藥效團(tuán)和分子對(duì)接相結(jié)合能夠快速有效地發(fā)現(xiàn)新穎活性候選化合物。

關(guān)鍵詞：HIV-1蛋白酶抑制劑；分子相似性；藥效團(tuán)；分子對(duì)接

獲得性免疫缺陷綜合征(AIDS)是由人類免疫缺陷病毒(HIV)感染引起的傳染病，對(duì)公共健康構(gòu)成嚴(yán)重威脅[1]。HIV蛋白前體由結(jié)構(gòu)基因gag、pol和env編碼而成[2]。HIV-1蛋白酶是HIV病毒粒子組裝、成熟過程的關(guān)鍵酶，其主要功能是特異性地裂解gag和gag-pol多聚蛋白前體形成具有活性的病毒結(jié)構(gòu)蛋白和酶[3,4]。HIV-1蛋白酶以二聚體形式發(fā)揮其活性，干擾二聚體形成或破壞已形成的二聚體能夠有效地抑制其活性[5]。因此，HIV-1蛋白酶可作為抗艾滋病治療的重要靶點(diǎn)，其抑制劑也成為目前抗HIV藥物復(fù)合療法的重要組成部分[6,7]。近二十年，F(xiàn)DA批準(zhǔn)了9個(gè)HIV-1蛋白酶抑制劑藥物，其中包括最新的Darunavir。但是HIV病毒變異率高，耐藥病毒株不斷出現(xiàn)，而且耐藥毒株的耐藥性越來越強(qiáng)[8,9]。Darunavir是至今研發(fā)出的活性最高的HIV-1蛋白酶抑制劑，也是唯一一個(gè)被美國FDA 批準(zhǔn)作為治療耐藥型HIV-1 的蛋白酶抑制劑。但是，Darunavir具有肝臟毒性，特別是對(duì)伴有乙型和丙型肝炎病毒感染的HIV患者[10]。另外，Darunavir還可誘發(fā)嚴(yán)重的過敏反應(yīng)[11]。因此，開發(fā)新型的、安全性高的、抗耐藥性的HIV-1蛋白酶抑制劑具有重要意義。

目前，虛擬篩選包括分子對(duì)接、藥效團(tuán)、三維定量構(gòu)效關(guān)系以及機(jī)器學(xué)習(xí)，已經(jīng)成功應(yīng)用于先導(dǎo)化合物的發(fā)現(xiàn)[12-15]。雖然單獨(dú)的分子對(duì)接以及單個(gè)分子為模板的藥效團(tuán)模型已被用于HIV-1 蛋白酶抑制劑的發(fā)現(xiàn)[16,17]，但這并不能滿足海量化合物庫篩選和新穎活性骨架捕獲的需要。DataWarrior是全新的集分子描述符計(jì)算和數(shù)據(jù)分析于一體的化學(xué)信息學(xué)軟件，能夠利用自帶的分子描述符對(duì)小分子化合物進(jìn)行快速的相似性分析和聚類[18]。LigandScout是基于3D-藥效團(tuán)模型進(jìn)行精確虛擬篩選的集成平臺(tái)。支持基于配體或基于結(jié)構(gòu)的藥效團(tuán)建模，還包含基于藥效團(tuán)的疊合，共有藥效團(tuán)特征的創(chuàng)建，以及自動(dòng)生成ROC曲線以便進(jìn)行性能評(píng)估。其篩選速度快，特別適合于大化合物庫的篩選[19]。

因此，本研究采用分子相似性，藥效團(tuán)和分子對(duì)接相結(jié)合的策略。首先，對(duì)已知的HIV-1蛋白酶抑制劑進(jìn)行相似性分析，選出代表性化合物作為模板建立藥效團(tuán)篩選模型，根據(jù)Darunavir和HIV-1蛋白酶間的相互作用指導(dǎo)藥效團(tuán)模型的藥效團(tuán)特征優(yōu)化，并利用已報(bào)道的HIV-1蛋白酶抑制劑作為訓(xùn)練集對(duì)模型篩選能力進(jìn)行評(píng)價(jià)。最后，我們將優(yōu)化好的模型應(yīng)用于化合物庫篩選，對(duì)打分較高的化合物進(jìn)行分子對(duì)接，觀察化合物與HIV-1 蛋白酶的結(jié)合情況，選取活性候選化合物進(jìn)行結(jié)合自由能計(jì)算和抗突變性分析。

1材料與方法

1.1　軟件和材料

采用軟件DataWarrior (開源軟件)，LigandScout 3.12 (Inte: Ligand)，PyMOL (開源軟件)，Gold 軟件 (CCDC Software Ltd)，MOE (Chemical Computing Group Inc.) 和Discovery Studio Visualization (Accelrys)。蛋白質(zhì)晶體結(jié)構(gòu)下載于PDB (Protein Data Bank)，收集文獻(xiàn)已報(bào)道的HIV-1蛋白酶抑制劑組成訓(xùn)練集(Training Set)[20-23]， 類藥性偽活性化合物庫(Drug-Like Ligand Decoys Set)來自Schr?dinger，化合物篩選庫采用TimTec的多樣化合物篩選庫(Diversity-screening-set)。

1.2　藥效團(tuán)模型的建立、優(yōu)化及驗(yàn)證

化合物包括Training Set、Decoys Set 和Screening Set均采用LigandScout 中的OMEGA生成ldb 格式的化合物構(gòu)象數(shù)據(jù)庫，單個(gè)化合物最多構(gòu)象限制為300。Darunavir的HIV-1蛋白晶體復(fù)合物(PDB ID: 4LL3)從PDB獲得，并利用PyMOL對(duì)Darunavir和HIV-1蛋白酶的相互作用模式進(jìn)行分析。采用DataWarrior對(duì)141個(gè)已知的HIV-1蛋白酶抑制劑進(jìn)行相似性分析，選取7個(gè)代表性的化合物進(jìn)行藥效團(tuán)模型構(gòu)建，并根據(jù)Darunavir和HIV-1蛋白酶的相互作用分析指導(dǎo)藥效團(tuán)特征的優(yōu)化。在模型建立和優(yōu)化過程中，我們利用訓(xùn)練集和類藥性偽活性化合物庫對(duì)模型進(jìn)行評(píng)價(jià)，考察模型對(duì)已知HIV-1蛋白酶抑制劑的富集能力。

1.3　基于藥效團(tuán)模型的化合物篩選和分子對(duì)接

將優(yōu)化好的藥效團(tuán)模型作為提問結(jié)構(gòu)，應(yīng)用LigandScout分子模擬軟件包對(duì)TimTec公司的多樣化合物篩選庫(含10 000個(gè)化合物)進(jìn)行試篩，得到與所建藥效團(tuán)模型配備較高的前500個(gè)化合物。再利用Gold軟件對(duì)這些化合物進(jìn)行分子對(duì)接，以HIV-1蛋白酶晶體結(jié)構(gòu)(PDB ID: 4LL3)作為受體篩選模型，定義復(fù)合物中Darunavir結(jié)合區(qū)域作為活性位點(diǎn)，綜合考慮受體—配體之間的極性和非極性相互作用評(píng)價(jià)對(duì)接結(jié)果。最后選取打分高于或相當(dāng)于Darunavir的化合物作為活性候選物。

1.4　結(jié)合自由能計(jì)算和抗突變性分析

現(xiàn)有的分子對(duì)接軟件很少或很難考慮受體大分子的柔性, 采用MOE軟件包中的MMFF94力場(chǎng)進(jìn)行能量優(yōu)化可以同時(shí)考慮受體和配體的柔性, 得到更為合理的結(jié)合自由能[24]。在HIV-1 蛋白酶抑制劑研發(fā)過程中，藥物學(xué)家已經(jīng)發(fā)現(xiàn)了至少50種HIV蛋白酶抗性突變，Darunavir 的抗耐藥性源于它能與蛋白酶活性中心的氨基酸骨架形成廣泛而強(qiáng)的氫鍵作用，而且大部分作用殘基并不位于突變位點(diǎn)，能有效抵抗蛋白酶局部突變所致的親和力下降[9]。因此分析活性候選物的氫鍵形成位點(diǎn)和蛋白酶突變位點(diǎn)，有助于了解化合物的抗突變性。

2結(jié)果與討論

2.1　蛋白晶體復(fù)合物及相互作用分析

從蛋白數(shù)據(jù)庫(PDB)下載Darunavir的HIV-1蛋白酶晶體結(jié)構(gòu)(PDB ID: 4LL3)。通過Discovery Studio Visualization對(duì)結(jié)構(gòu)進(jìn)行分析，可以看到晶體復(fù)合物中HIV-1蛋白酶由AB兩條鏈共同構(gòu)成了抑制劑Darunavir的結(jié)合口袋(見圖1a)。采用PyMOL軟件對(duì)結(jié)合口袋進(jìn)行相互作用分析，發(fā)現(xiàn)Darunavir與活性口袋匹配較好，能夠同HIV-1蛋白酶上的ASP-25，GLY-27，ASP-29以及ASP30氨基酸殘基產(chǎn)生強(qiáng)的氫鍵相互作用(見圖1b)，另外Darunavir 還可與HIV 蛋白酶的Leu23，Gly49，Ile50，Pro81，Vla82以及Ile84等殘基有較強(qiáng)的范德華相互作用，所有這些特性對(duì)于它的高活性和抗突變型蛋白酶的能力來說都是至關(guān)重要的。

圖1　Darunavir和HIV-1蛋白酶的晶體復(fù)合物結(jié)構(gòu)及其相互作用Fig.1　The crystal structure of Darunavir bound to HIV-1 protease and their interactions

2.2　藥效團(tuán)的構(gòu)建、驗(yàn)證和優(yōu)化

利用DataWarrior對(duì)已知的HIV-1蛋白酶抑制劑進(jìn)行相似性分析(見圖2a)，發(fā)現(xiàn)其大致可分為7類，每類的中心分子被選作模板分子。利用LigandScout軟件構(gòu)建藥效團(tuán)模型并進(jìn)行優(yōu)化和驗(yàn)證，優(yōu)化后的藥效團(tuán)模型如圖2b所示，其包括3個(gè)芳香環(huán)，4個(gè)氫鍵受體和2個(gè)氫鍵供體。采用訓(xùn)練集(由141個(gè)已知的HIV-1蛋白酶抑制劑組成) 和類藥性偽活性化合物庫 (999個(gè)來自Schr?dinger的類藥分子構(gòu)成) 對(duì)模型進(jìn)行評(píng)價(jià)，考察藥效團(tuán)模型將活性分子和偽活性分子(Decoys)分開的能力。通過計(jì)算ROC曲線下面積(Area Under Curves, AUC)和1%，5%，10%，100% 四個(gè)水平的富集因子(Enrichment Factors, EFs)作為模型驗(yàn)證的基準(zhǔn)。如圖3所示，最后優(yōu)化好的模型AUC為0.80，四個(gè)水平的EFs分別為：8.1，7.8， 5.2和2.8。

圖2　HIV-1蛋白酶抑制劑的相似性分析和以代表性分子為模板構(gòu)建的藥效團(tuán)模型Fig.2　The similarity analysis of known HIV-1 protease inhibitors using DataWarrior and the pharmacophorederived from seven representative HIV-1 protease inhibitors using LigandScout

圖3　藥效團(tuán)模型篩選能力的評(píng)價(jià)(AUC=0.80)Fig.3　The screening performance of the optimizedpharmacophore model (AUC=0.80)

2.3　藥效團(tuán)虛擬篩選和分子對(duì)接

首先利用優(yōu)化好的藥效團(tuán)模型對(duì)TimTec的Diversity-screening-set (10 000個(gè)化合物)進(jìn)行篩選，保留打分較高的前500個(gè)化合物，然后利用Gold軟件進(jìn)行分子對(duì)接。為了檢驗(yàn)對(duì)接方法的可靠性，選擇將原有配體Darunavir重新對(duì)接到HIV-1蛋白酶受體中，然后比較配體對(duì)接后預(yù)測(cè)構(gòu)象與晶體構(gòu)象的差異。結(jié)果預(yù)測(cè)構(gòu)象可以重現(xiàn)原晶體中的構(gòu)象，表明分子對(duì)接參數(shù)能夠有效用于化合物的進(jìn)一步篩選。最后從分子對(duì)接結(jié)果中選取打分高于或相當(dāng)于Darunavir的4個(gè)候選化合物做進(jìn)一步分析(見表1)。如圖4所示，對(duì)獲得的4個(gè)候選化合物與HIV-1蛋白酶間的相互作用進(jìn)行分析，發(fā)現(xiàn)候選化合物ST070845可與HIV-1蛋白酶的ARG-8，ASP-25，ASP-30以及GLY-48形成多個(gè)氫鍵，而且分子中的苯并吲哚母核可以很好地嵌入疏水口袋?；衔颯T088084中的三氮唑母核能夠同時(shí)和蛋白酶A、B兩條鏈的ILE50殘基形成強(qiáng)的氫鍵作用，另外其末端氨基亦可與蛋白酶中的GLY-48殘基形成氫鍵?；衔颯T025723與蛋白酶的ASP-25和GLY-27殘基產(chǎn)生氫鍵相互作用，這兩個(gè)殘基也是Darunavir的作用殘基?；衔颯T025679的酰胺結(jié)構(gòu)可與ILE-50形成氫鍵，同時(shí)其甲氧基也可與ASP-30殘基有氫鍵作用。

表1　化合物的分子對(duì)接打分，結(jié)合自由能和氫鍵作用

2.4　結(jié)合自由能計(jì)算、抗突變性分析和結(jié)構(gòu)新穎性評(píng)價(jià)

分子對(duì)接可以快速地對(duì)配體小分子進(jìn)行一定空間的構(gòu)象搜索，以保證在合理的時(shí)間內(nèi)從大量數(shù)據(jù)庫中搜索到類藥的化合物。但是，現(xiàn)有的分子對(duì)接軟件很少或很難考慮受體大分子的柔性，難以在對(duì)接時(shí)適時(shí)調(diào)整受體的構(gòu)象。對(duì)受體—配體復(fù)合物進(jìn)行能量優(yōu)化和結(jié)合自由能計(jì)算可以同時(shí)考慮受體和配體的柔性，受體活性中心的構(gòu)象可以隨著配體構(gòu)象的改變而改變, 真正做到“誘導(dǎo)契合”， 得到更為合理的結(jié)合自由能[25]。通過結(jié)合自由能計(jì)算，發(fā)現(xiàn)這些化合物的結(jié)合自由能要高于Darunavir，但是化合物ST025679的結(jié)合自由能非常接近于Darunavir，另外化合物ST088084和ST025723也顯示出較低的結(jié)合自由能，而化合物ST070845雖然對(duì)接打分較高，但也呈現(xiàn)出較高的結(jié)合自由能。耐藥性是艾滋病治療過程中的一個(gè)重要問題，主要原因是HIV-1蛋白酶序列中的某些氨基酸發(fā)生了抗性突變，按照它們?cè)诘鞍酌附Y(jié)構(gòu)中的位置，又可以分為活性位點(diǎn)突變和非活性位點(diǎn)突變，其中活性位點(diǎn)突變占了較大比重。根據(jù)HIV蛋白酶常見突變位點(diǎn)和4個(gè)新活性候選化合物的作用殘基，其中化合物ST088084和ST025679可能會(huì)因I50V突變導(dǎo)致其活性下降[26,27]。最后，采用FCFP_4和MACCS兩種分子指紋對(duì)化合物ST025723與FDA批準(zhǔn)上市的9個(gè)HIV-1蛋白酶抑制藥物進(jìn)行相似性分析，發(fā)現(xiàn)ST025723與藥物Atazanavir最為相似，其相似系數(shù)分別為0.34和0.45，說明ST025723結(jié)構(gòu)新穎，不同于已有的HIV-1蛋白酶抑制藥物。因此，綜合考慮結(jié)合自由能、抗突變性和結(jié)構(gòu)新穎性，化合物ST025723具有更潛在的深入研究價(jià)值。

圖4 　HIV-1蛋白酶抑制劑候選化合物和HIV-1蛋白酶的相互作用模式Fig.4　The interactions between candidate inhibitors and HIV-1 protease

3結(jié)論

HIV-1蛋白酶是HIV復(fù)制過程中的關(guān)鍵酶，將其作為藥物靶點(diǎn)并開發(fā)其抑制劑具有廣闊的前景。本研究采用分子相似性、藥效團(tuán)和分子對(duì)接相結(jié)合的虛擬篩選策略，使各自的優(yōu)點(diǎn)得到了很好的融合。由于藥效團(tuán)篩選速度快，適合于對(duì)大化合物庫進(jìn)行快速篩選。而后采用分子對(duì)接，對(duì)藥效團(tuán)篩選結(jié)果進(jìn)行進(jìn)一步的篩選，而且還能對(duì)獲得的候選化合物進(jìn)行相互作用模式分析，最后采用結(jié)合自由能計(jì)算，同時(shí)考慮受體和配體的柔性，使活性評(píng)價(jià)更為合理和準(zhǔn)確。此外，通過化合物結(jié)合位點(diǎn)和HIV-1蛋白酶常見突變位點(diǎn)進(jìn)行抗突變性分析，進(jìn)一步從篩選到的活性候選化合物中挑選出具有抗耐藥性的活性化合物。該策略確保了篩選結(jié)果的質(zhì)量，達(dá)到在更短時(shí)間內(nèi)更高效地發(fā)現(xiàn)抗耐藥活性候選化合物的目的。根據(jù)對(duì)TimTec多樣化合物庫篩選結(jié)果，4個(gè)活性最好的候選化合物，雖然結(jié)構(gòu)和現(xiàn)有HIV-1蛋白酶抑制劑Darunavir存在較大差異，但其作用模式卻類似，這也充分展現(xiàn)了組合藥效團(tuán)模型在捕獲新穎活性骨架方面的優(yōu)勢(shì)。

參考文獻(xiàn)

[1]MOORE C,CIARALD E.Quality of care and service expansion for HIV care and treatment[J].Current HIV/AIDS Reports, 2015, 12(2): 223-230.

[2]MCGETTIGAN J, NAPER K, ORENSTEIN J, et al. Functional human immunodeficiency virus type 1 (HIV-1) Gag-Pol or HIV-1 Gag-Pol and env expressed from a single rhabdovirus-based vaccine vector genome[J].Journal of Virology, 2003, 77(20): 10889-10899.

[3]ADAMSON C. Protease-mediated maturation of HIV:inhibitors of protease and the maturation process[J].Molecular Biology International,2012,604261: 1-13.

[4]袁嘯塵, 鈕冰, 尹京苑. 基于二層特征篩選的HIV-1蛋白酶特異位點(diǎn)預(yù)測(cè)[J]. 生物信息學(xué), 2012, 10(4): 274-279.

CHEN Yuanxiao, NIU Bing, YIN Jingyuan. HIV-1 protease cleavage site prediction based on two stage feature selection method[J]. Chinese Journal of Bioinformatics, 2012, 10(4): 274-279.

[5]HAYASHI H,TAKAMUNE N,NIRASAWA T,et al. Dimerization of HIV-1 protease occurs through two steps relating to the mechanism of protease dimerization inhibition by darunavir[J].Proceedings of the National Academy of Sciences,2014,111(33):12234-12239.

[6]DEBOUCK C. The HIV-1 protease as a therapeutic target for AIDS[J].AIDS Research and Human Retroviruses, 1992, 8(2):153-164.

[7]POTEMPA M, LEE S, WOLFENDEN R,et al. The triple threat of HIV-1 protease inhibitors[J].Current Topics in Microbiology and Immunology, 2015, 389: 203-241.

[8]HOU L, WANG H, DUAN S, et al. The prevalence and determinants of drug-resistance-associated mutations in the HIV-1infected MSM population of Henan Province in China[J].Archives of Virology,2015,160(8):2051-2061.

[9]郝格非,楊光富. 新型非擬肽磺胺類HIV蛋白酶抑制劑Darunavir 的研究進(jìn)展[J].有機(jī)化學(xué),2008,28(9):1545-1552.

HAO Gefei,YANG Guangfu. Progress in the new generation nonpeptide sulfonamide anti-HIV protease inhibitor darunavir[J].Chinese Journal of Organic Chemistry, 2008, 28(9): 1545-1552.

[10]VISPO E.Warning on hepatotoxicity of darunavir[J].AIDS Reviews, 2008, 10(1): 62-63.

[11]LIN K,CHENG C,YANG C,et al.Skin rash related to once-daily boosted darunavir-containing antiretroviral therapy in HIV-infected Taiwanese:incidence and associated factor[J].Journal of Infection and Chemotherapy,2014,20(8):465-470.

[12]LAVECCHIA A.Machine-learning approaches in drug discovery: methods and applications[J].Drug Discovery Today,2015,20(3):318-331.

[13]VUORINEN A,SCHUSTER D.Methods for generating and applying pharmacophore models as virtual screening filters and for bioactivity profiling[J].Methods,2014,71: 113-134.

[14]宋新蕊, 李達(dá), 陳潔, 等.計(jì)算機(jī)輔助藥物篩選平臺(tái)及應(yīng)用[J].生物信息學(xué),2014,12(4):300-304.

SONG Xinrui,LI Da,CHEN Jie, et al.Computer aided drug screening platform and its application[J]. Chinese Journal of Bioinformatics, 2014, 12(4): 300-304.

[15]LAGARDE N,ZAGURY J,MONTES M.Benchmarking data sets for the evaluation of virtual ligand screening methods: review and perspectives[J].Journal of Chemical Information and Modeling,2015,55(7):1297-1307.

[16]谷萬港,張麗,張旋. 基于結(jié)構(gòu)的HIV-1蛋白酶抑制劑的虛擬篩選[J].昆明醫(yī)科大學(xué)學(xué)報(bào),2013,34(8):19-22.

GU Wangang, ZHANG Li,ZHANG Xuan.Discovery of protease inhibitors of HIV-1 through structure based virtual screening[J]. Journal of Kunming Medical University,2013,34(8):19-22.

[17]繆有盼, 李愛秀, 劉 濤,等. 基于“底物包膜”假說篩選新型HIV-1蛋白酶抑制劑[J].化學(xué)學(xué)報(bào),2009, 67(10): 1098-1102.

MIAO Youpan,LI Aixiu,LIU Tao,et al.Screening novel HIV-1 protease inhibitors based on the substrate envelope hypothesis[J].Acta Chimica Sinica,2009,67(10):1098-1102.

[18]SANDER T,FREYSS J,VON KORFF M, et al.DataWarrior: an open-source program for chemistry aware data visualization and analysis[J].Journal of Chemical Information and Modeling,2015,55(2):460-473.

[19]WOLBER G,LANGER T.LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters[J].Journal of Chemical Information and Modeling,2005,45(1):160-169.

[20]WEI Y,LI J,CHEN Z, et al.Multistage virtual screening and identification of novel HIV-1 protease inhibitors by integrating SVM,shape,pharmacophore and docking methods[J]. European Journal of Medicinal Chemistry,2015,101:409-418.

[21]PARAI M,HUGGINS D,CAO H,et al.Design,synthesis,and biological and structural evaluations of novel HIV-1 protease inhibitors to combat drug resistance[J]. Journal of Medicinal Chemistry, 2012,55(14):6328-6341.

[22]ISLAM M,PILLAY T. Exploration of the structural requirements of HIV-protease inhibitors using pharmacophore, virtual screening and molecular docking approaches for lead identification[J].Journal of Molecular Graphics and Modelling,2015, 56: 20-30.

[23]ALI A,REDDY G,NALAM M, et al.Structure-based design,synthesis and structure-activity relationship studies of HIV-1 protease inhibitors incorporating phenyloxazolidinones[J].Journal of Medicinal Chemistry, 2010, 53(21): 7699-7708.

[24]DEETH R, FEY N, WILLIAMS-HUBBARD B.DommiMOE:an implementation of ligand field molecular mechanics in the molecular operating environment[J].Journal of Computational Chemistry,2005,26(2): 123-130.

[25]曹 冉, 李 偉, 孫漢資, 等. 計(jì)算化學(xué)方法在基于受體結(jié)構(gòu)的藥物分子設(shè)計(jì)中的基礎(chǔ)理論及應(yīng)用[J]. 藥學(xué)學(xué)報(bào), 2013, 48 (7): 1041-1052.

CAO Ran,LI Wei,SUN Hanzi,et al.Computational chemistry in structure-based drug design[J].Acta Pharmaceutica Sinica,2013,48 (7):1041-1052.

[26]PARTALEDIS J,YAMAGUCHI K,TISDALE M,et al. In vitro selection and characterization of human immunodeficiency virus type 1 (HIV-1) isolates with reduced sensitivity to hydroxyethyl amino sulfonamide inhibitors of HIV-1 aspartyl protease[J].Journal of Virology, 1995, 69(9):5228-5235.

[27]MITTAL S,BANDARANAYAKE R,KING N,et al.Structural and thermodynamic basis of amprenavir/darunavir and atazanavir resistance in HIV-1 protease with mutations at residue 50[J].Journal of Virology,2013,87(8):4176-4184.

Discovery of new HIV-1 protease inhibitors by integrating molecular

similarity,pharmacophore and docking methods

WANG Ying, TANG Guorong, LIU Shunan*

(DepartmentofPharmacy,The476thHospitalofPLA,Fuzhou350002,China)

Abstract：An efficient and accurate method was established for HIV-1 protease inhibitor screening by integrating molecular similarity, pharmacophore and docking methods. First, the similarity of known HIV-1 protease inhibitors was analyzed, and representative inhibitors were chosen as templates to build pharmacophore model for database screening. Five hundred compounds were picked out from 10 000 compounds by pharmacophore screening, which were further submitted to docking experiment. Four new candidate compounds were obtained, and their binding free energies and antimutagenic abilities were further evaluated. The candidate compound ST025723 showed good interactions with HIV-1 protease and antimutagenic ability,and deserved to be further studied.Therefore, the combination of molecular similarity, pharmacophore, and docking could be an efficient strategy to discover novel active candidates.

Keywords：HIV-1 protease inhibitor;Molecular similarity;Pharmacophore;Docking

中圖分類號(hào)：R91

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1672-5565(2015)04-244-07

doi:10.3969/j.issn.1672-5565.2015.04.07

作者簡(jiǎn)介：汪瀅，女，藥師，研究方向：化合物數(shù)據(jù)庫的構(gòu)建與虛擬篩選；E-mail：yingwang912@163.com;*通信作者：劉澍楠，男，主任藥師，研究方向：創(chuàng)新藥物研究與藥學(xué)信息學(xué)；E-mail：liusnfj@gmail.com.

收稿日期：2015-10-14;修回日期：2015-11-17.