劉炳旺，馬英，李彬寒，王潤玲

（1.天津醫(yī)科大學藥學院，天津市臨床藥物關鍵技術重點實驗室，天津300070；2.天津醫(yī)科大學腫瘤醫(yī)院分子影像及核醫(yī)學診療科，國家腫瘤臨床醫(yī)學研究中心，天津市“腫瘤防治”重點實驗室，天津 300060）

論著

膽固醇酯轉(zhuǎn)運蛋白抑制劑的分子設計

劉炳旺1，2，馬英1，李彬寒1，王潤玲1

（1.天津醫(yī)科大學藥學院，天津市臨床藥物關鍵技術重點實驗室，天津300070；2.天津醫(yī)科大學腫瘤醫(yī)院分子影像及核醫(yī)學診療科，國家腫瘤臨床醫(yī)學研究中心，天津市“腫瘤防治”重點實驗室，天津 300060）

目的：設計針對膽固醇酯轉(zhuǎn)運蛋白(CETP)的抑制劑，提高其降低血脂水平，為以后相關疾病的治療提供科學的方法和依據(jù)。方法：以現(xiàn)有torcetrapib抑制劑為模板，利用Schrodinger Suite 2009軟件中的“Core Hopping”模塊進行修飾，對修飾后的結(jié)構(gòu)用“Core Hopping”模塊改造結(jié)構(gòu)；利用Schrodinger Suite2009中的Glide模塊對druglike數(shù)據(jù)庫中10萬個化合物進行高通量虛擬篩選，得到結(jié)構(gòu)進行修飾。利用Glide模塊對所建立的化合物庫進行虛擬篩選。應用Qikprop模塊做ADME（吸收、分布、代謝、排泄）預測來推測這些化合物的成藥可能性。結(jié)果:對torcetrapib結(jié)構(gòu)進行修飾，得到8個較好的化合物；通過對druglike數(shù)據(jù)庫進行高通量虛擬篩選得到1個化合物-ZINC26608950，修飾后得到6個化合物。用分子對接方法分析了新抑制劑和CETP的相互作用機制，與現(xiàn)有的抑制劑torcetrapib相比有更好的結(jié)合能力。通過ADME預測得出設計出的化合物均符合類藥5原則。結(jié)論：通過兩種方法得到14個CETP抑制劑，分子對接初步證明其可以抑制CETP，此為CETP抑制劑的結(jié)構(gòu)改造及活性測定奠定了基礎。

膽固醇酯轉(zhuǎn)運蛋白；抑制劑；藥物設計；分子對接

心血管疾病具有“發(fā)病率高、致殘率高、死亡率高、復發(fā)率高、并發(fā)癥多”的特點[1-2]，嚴重危害到人類健康，已位于人類死亡病因之首，而高血壓、冠狀動脈粥樣硬化是其病因之一。研究治療冠狀動脈粥樣硬化藥物對預防和治療心血管疾病有著重要意義[3]。研究發(fā)現(xiàn)膽固醇酯轉(zhuǎn)運蛋白（CETP）可作為治療動脈硬化的新靶點，一個基因突變致CETP缺乏的家系，高密度脂蛋白（HDL）水平升高，同時冠心病發(fā)病率明顯下降，并且當抑制CETP的活性時，能有效提高血漿內(nèi)HDL含量水平[4-6]，所以CETP抑制劑被認為是冠狀動脈粥樣硬化(AS)最有治療潛力的藥物之一[7-12]。本研究旨在通過兩種方法-高通量虛擬篩選和結(jié)構(gòu)修飾的方法設計出新的CETP抑制劑。

1 材料與方法

1.1 材料 軟件及數(shù)據(jù)庫在DELL T5500計算工作站上進行操作，涉及的軟件為Schrodinger 2009，模塊包括 ligPrep,QikProp,Combiglide,Glide。在ChemBioOffice2010上進行繪圖，片段分子來自已有碎片庫。受體來自PDBbank數(shù)據(jù)庫。

1.2 方法

1.2.1 受體結(jié)構(gòu)的準備 蛋白結(jié)構(gòu)預處理由Schrodinger Suite 2009軟件中的Protein preparation wizard模塊完成。首先對CETP蛋白質(zhì)結(jié)構(gòu)（PDB ID:4EWS）[13]進行修正，包括修正化學鍵，加氫，處理金屬離子，修正N端和C端的殘基結(jié)構(gòu)錯誤，除去蛋白質(zhì)晶體結(jié)構(gòu)中的水分子和雜原子，只保留蛋白分子與原始配體。再對蛋白質(zhì)氫鍵進行優(yōu)化，優(yōu)化參數(shù)為Sample water orientations，優(yōu)化力場選擇OPLS 2005[14]，以RMSD 0.5 ?為收斂標準。

1.2.2 盒子的產(chǎn)生 以原始配體小分子torcetrapib為中心，大小選擇為11?生成對接盒子文件。為確保盒子文件能夠適合模擬體系，對盒子文件進行驗證。本實驗通過將小分子torcetrapib對接到CETP（4EWS）的配體結(jié)合區(qū)，再重新與晶體結(jié)構(gòu)中的配體構(gòu)象疊合比對來驗證對接方法和盒子文件的選擇是否合理。

1.2.3 配體準備 對分子進行配體優(yōu)化處理，此步在 LigPre模塊中完成。小分子優(yōu)化的力場為OPLS_2005，產(chǎn)生pH值7.0±2.0下的離子化狀態(tài)，去鹽，產(chǎn)生各種可能的互變異構(gòu)體。

1.2.4 分子對接 應用Glide[15]模塊進行分子對接。以原始配體小分子為盒子中心，大小為11×11×11?自動生成對接盒子文件。虛擬篩選中對接精度為HTVS（高通量篩選）。對接精度為SP（標準精度），采用柔性對接，柔性對接時要求考慮各種環(huán)的柔性情況，要求對接函數(shù)在打分時對非平面型的酰胺鍵構(gòu)象進行罰分，其它參數(shù)為默認值。

1.2.5 結(jié)構(gòu)修飾 “Core Hopping”模塊具有片段替換和分子剛性對接兩大功能，可以發(fā)現(xiàn)全新的拓撲結(jié)構(gòu)骨架和活性分子。在“Core Hopping”過程中，第一步是定義模板化合物，在Define Combinations step中定義需要替換的片段。第二步將定義好的受體盒子導入Receptor grid file中。第三步導入分子片段，分子片段來自于實驗室已有片段數(shù)據(jù)庫和軟件中自帶的片段數(shù)據(jù)庫。小分子片段會自動連接到母核上生成新的化合物數(shù)據(jù)庫，得到的化合物數(shù)據(jù)庫與受體進行剛性對接。依據(jù)對接得分進行下一步篩選，得分的絕對值越大，配體與受體的親和力越強，結(jié)合力越好。

1.2.6 ADME預測 QikProp模塊可以預測化合物，可以對整個分子的3D結(jié)構(gòu)預測，預測性質(zhì)包括：辛醇/水和水/氣體的log Ps,log S,log BB；整個CNS活性；Caco-2和MDCK細胞滲透性；血清白蛋白結(jié)合活性log Khsa；HERG的K離子通道阻斷性log IC50。本實驗對logPo/w、PSA、logS和PMDCK進行了預測。

2 結(jié)果

2.1 對接方法驗證 對接后的結(jié)構(gòu)與晶體中配體原始構(gòu)象的均方根偏差（RMSD）值為0.8 ?。說明盒子文件能夠適用模擬體系。圖 1為 CETP與torcetrapib的結(jié)合模式圖。

圖1 CETP與compound1、torcetrapib的結(jié)合模式圖Fig 1 The conformation obtained by docking torcetrapib and compound1 to CETP

2.2 CETP抑制劑

2.2.1 對torcetrapib結(jié)構(gòu)進行修飾 本實驗選取torcetrapib作為結(jié)構(gòu)修飾的母核。具體修飾過程如圖2所示：首先對A片段進行替換旨在找到結(jié)合更強的片段，找出2個優(yōu)選片段A1、A2。然后，對B片段進行替換以期找到疏水作用更強的片段，找出4個優(yōu)選片段B1、B2、B3、B4。這樣就產(chǎn)生了2×4=8個化合物。根據(jù)這個分子修飾的原則，表1中列出了修飾后的8個torcetrapib的類似物分子，包括torcetrapib本身在內(nèi)的對接結(jié)果。

圖2 以torcetrapib為母核的結(jié)構(gòu)修飾流程圖Fig 2 A flowchart to show the agonists design process of torcetrapib

2.2.2 高通量虛擬篩選 對druglike[16]數(shù)據(jù)庫中10 000個化合物進行篩選，選出能與受體CETP的4個結(jié)合口袋結(jié)合的化合物最后得到1個化合物ZINC26608950。此化合物能結(jié)合到第2,3,4口袋。由此，選ZINC26608950作為進一步修改的候選化合物。ZINC26608950的缺點是沒有向CETP受體的第一結(jié)合空腔延伸，沒能與CETP受體結(jié)合區(qū)形成較好的匹配，所以可以嘗試設計出能夠占據(jù)CETP受體4個結(jié)合區(qū)的高效抑制劑。對此結(jié)構(gòu)應用“Core Hopping”進行結(jié)構(gòu)改造得出6個較好的結(jié)構(gòu),見表1。具體修飾過程如圖3所示：在環(huán)上添加分子片段以提高與受體的疏水區(qū)的匹配性，找出6個優(yōu)選片段R1、R2、R3、R4、R5和R6。對接結(jié)果顯示：分子結(jié)構(gòu)與CETP具有較好的匹配性。圖4為torcetrapib、comp-vs#1與CETP結(jié)合模式圖。

圖3 對ZINC26608950的結(jié)構(gòu)修飾流程圖Fig 3 A flowchart to show the agonists design process of ZINC26608950

圖4 Torcetrapib、comp-vs#1與CETP結(jié)合模式圖Fig 4 The conformation obtained by docking torcetrapib and comp-vs#1 to CETP

2.3 ADME預測結(jié)果 應用Schrodinger2009軟件包中的QikProp模塊對與藥效相關的5個性質(zhì)進行預測。結(jié)果如表1所示。新設計的化合物的logPo/w、PSA、logS、PMDCK均在合理的范圍內(nèi)。

表1 ADME預測和對接得分Tab 1 Physiochemical descriptors calculated by QP and docking score

3 討論

3.1 CETP受體抑制劑的設計 靶點的選擇對化合物的設計起到至關重要的作用。CETP靶點蛋白晶體結(jié)構(gòu)從布魯克海文蛋白數(shù)據(jù)庫下載，遵循以下條件：所選靶點晶體結(jié)構(gòu)為人源的；晶體結(jié)構(gòu)中氨基酸殘基沒有突變；晶體結(jié)構(gòu)分辨率小于2.6?[17]。

Torcetrapib[18]是由美國輝瑞公司研發(fā)的第一個CETP抑制劑，也是第一個進行大規(guī)模終點事件評價臨床試驗的CETP抑制劑。torcetrapib使收縮壓升高約5 mmHg，同時升高血鈉、碳酸氫鹽及醛固酮水平。似乎這些脫靶效應導致了死亡率的增加，后續(xù)研究顯示這些不良作用與CETP抑制機制無關，即其他CETP抑制劑不一定有這樣的作用。CETP抑制劑torcetrapib具有高脂溶性（cLogP=8.395），極低水溶性，這就使化合物產(chǎn)生了脫靶現(xiàn)象。torcetrapib有兩點不符合成藥五規(guī)則：分子質(zhì)量超過500和高脂溶性。所以在設計化合物的時候要避免此現(xiàn)象產(chǎn)生?，F(xiàn)今，一個成功用作口服藥物的化合物必須具有足夠的水溶性，以及和其臨床藥效相關的通透性。準確的ADME預測對進一步合成起到指導作用。本實驗設計的化合物 logPo/w、PSA、logS、PMDCK均在合理的范圍內(nèi)。盡管logPo/w越高，結(jié)合能力越強，但不能過度地增強logPo/w水平，這樣會導致藥物在體內(nèi)分布較差。這些化合物的創(chuàng)新性在于不僅均符合成藥五規(guī)則還能與CETP有很好的相互作用，在理論上這些化合物的結(jié)合能力超過torcetrapib。

分子對接初步證明其可以抑制CETP，具有理論活性，表明基于CETP治療冠狀動脈粥樣硬化具有一定可行性，為治療冠狀動脈粥樣硬化的先導化合物的結(jié)構(gòu)改造及活性測定奠定了基礎。

[1] Zhang L,Yan F,Zhang S,et al.Structural basis of transfer between lipoproteins by cholesteryl ester transfer protein[J].Nat Chem Biol, 2012,8(4):342

[2] Bishop B M.Systematic review of CETP inhibitors for increasing high-density lipoprotein cholesterol:where do these agents stand in the approval process[J]Am J Ther,2015,22(2):147

[3] Chapman M J,Wilfried L G,Guerin M,et al.Cholesteryl ester transfer protein:at the heart of the action of lipid-modulating therapy with statins,fibrates,niacin,and cholesteryl ester transfer protein inhibitors[J].Eur Heart J,2010,31(2):149

[4] Barter P J,Brewer H B,Chapman M J,et al.Cholesteryl ester transfer protein:a novel target for raising HDL and inhibiting atherosclerosis[J].Arterioscler Thromb Vasc Biol,2003,23(2):160

[5] Ding Y Y,Zhang W,Zhang M Q,et al.Functional and association studies of the cholesteryl ester transfer protein (CETP)gene in a Wannan Black pig model[J].Anim Genet,2015,46(6):702

[6] Morton R E,Izem L.Modification of CETP function by changing its substrate preference:a new paradigm for CETP drug design[J].J Lipid Res,2015,56(3):612

[7]楊征,邱敏.丹參酮IIA的心血管作用及機制研究進[J].中國動脈硬化雜志,2011，19（4）:372

[8] Collins H L,Drazul-Schrader D,Sulpizio A,et al.L-Carnitine intake and high trimethylamine N-oxide plasma levels correlate with low aortic lesions in ApoE transgenic mice expressing CETP[J]. Atherosclerosis,2015,244:29

[9] Miyosawa K,Watanabe Y,Murakami K,et al.New CETP inhibitor K-312 reduces PCSK9 expression:a potential effect on LDL cholesterol metabolism[J].Am J Physiol Endocrinol Metab,2015, 309(2):177

[10]Jeenduang N,Porntadavity S,Nuinoon M,et al.Studies of the CETP TaqIB and ApoE polymorphisms in southern thai subjects with the metabolic syndrome[J].Biochem Genet,2015,53（7/8）:184

[11]Hitchcock E,Patankar J V,Tyson C,et al.A novel microdeletion affecting the CETP gene raises HDL-associated cholesterol levels [J].Clin Genet,2015,1

[12]Greene D J,Izem L,Morton R E.Defective triglyceride biosynthesis in CETP-deficient SW872 cells[J].J Lipid Res,2015,56（9）:1669 [13]Liu S P,Mistry A,Jennifer M,et al.Crystal structures of cholesteryl ester transfer protein incomplex with inhibitors[J].J Biol chem,2015, 287（44）:37321

[14]Oostenbrink C,Soares T A,Vander Vegt N F,et al.Validation of the 53A6 GROMOS force field[J].Eur Biophys J,2005,34(4):273

[15]Halgren T A,Murphy R B,Friesner R A,et al.Glide:a new approach for rapid,accurate docking and scoring.2.Enrichment factors in database screening[J].J Med Chem,2004,47(7):1750

[16]Irwin J J,Shoichet B K.ZINC--a free database of commercially availablecompoundsforvirtualscreening[J].JChemInf Model,2005, 45（1）:177

[17]洪淵．基質(zhì)金屬蛋白酶及其抑制劑的分子對接研究[D]．吉林:吉林大學，2007

[18]Morehouse L A,Sugarman E D,Bourassa P A,et al.Inhibition of CETP activity by torcetrapib reduces susceptibility to diet-induced atherosclerosis in NZW rabbits[J].J Lipid Res,2007,48（6）:1263

（2015-10-09收稿）

Designing novel inhibitors by targeting cholesteryl ester transfer protein

LIU Bing-wang1,2,MA Ying1,LI Bin-han1,WANG Run-ling1
（1.School of Pharmacy,Tianjin Medical University,Tianjin Key Laboratory on Technologies Enabling Development of Clinical Therapeutics and Diagnostics,Tianjin 300070,China；2.Department of Molecular Imaging and Nuclear Medicine,Cancer Institute and Hospital,Tianjin Medical University,National Clinical Research Center for Cancer,Key Laboratory of Cancer Prevention and Therapy, Tianjin 300060,China)

Objective：To develop highly active cholesteryl ester transfer protein（CETP）novel inhibitors to facilitate further research.Methods：In this study,torcetrapib was modified by means of“core hopping”;the small molecule drug-like database was screened by Glide of Schrodinger Suite 2009 and then was modified.It was further validated by the outcomes of their ADME (absorption,distribution, metabolism,and excretion)predictions that the new agonists could have high potential to become drug candidates.Results：Eight compounds were discovered to inhibit CETP based on torcetrapib by means of“core hopping”.One CETP inhibitor-ZINC26608950 was discovered by high throughput virtual screening,and 6 compounds were discovered to inhibit CETP based on ZINC26608950.It was observed by docking that these novel inhibitors had more favorable conformation for binding to the receptor than torcetrapib.It was predicted that the values for these novel candidates were all within the reasonable ranges.Conclusion：The binding affinity of the designed molecules is theoretically higher than the torcetrapib in CETP.And the new inhibitors may become potential drug hits.

cholesteryl ester transfer protein;inhibitors;computer aid design drug;molecular docking

R914.2

A

1006－8147（2016）02-0122-04

劉炳旺(1981-)，男，碩士在讀，研究方向：計算機輔助藥物設計；通信作者：王潤玲，E-mail:wangrunling@tijmu.edu.cn。