计算机辅助CDK6抑制剂先导化合物初步设计及化合物活性的测定
【关键词】 机辅助药物设计
Computerassisted primary design of lead compounds of CDK6 inhibitors and determination of their bioactivity
【Abstract】 AIM: To primarily design the lead compounds of CDK16 inhibitors based on CDK6 structure and to test their bioactivity. METHODS: Using LigBuilderv1.2 ［a kind of computerassisted drug design (CADD) program］, we tried to design new lead compunds of CDK6 inhibitors by method of de novo design and according to Lipinski rule, and test their IC50 by MTT assay. RESULTS: We got the structures of 50 new chemical compounds as the lead compounds for further drug filtering, and 3 of 8 compunds we synthesized successfully had a bioactivity between 220 and 310 μmol/L. CONCLUSION: The new compunds have some inhibitory activity and the efficiency of finding lead compounds can be improved with the help of CADD.
【Keywords】 CADD; CDK6 inhibitor; IC50
【摘要】 目的： 初步设计CDK6抑制剂先导化合物，并测定化合物活性. 方法：利用计算机辅助药物设计（CADD）程序LigBuilderv1.2，依据Lipinski法则，采用从头设计、片断生长的方法，在Linux操作系统下，设计全新CDK6抑制剂化合物；MTT法测定化合物IC50. 结果：得到50种化合物结构并成功合成8种，经测定其中3种有活性（IC50在220~310 μmol/L间）. 结论：新型CDK6抑制剂化合物具有一定抑制活性，CADD可以大大提高先导化合物发现的效率，加快新药研究的步代.
【关键词】 计算机辅助药物设计; CDK6抑制剂; IC50
0引言
根据肿瘤干细胞理论［1-3］，干细胞相关基因及其编码蛋白是对于肿瘤的发生、维持具有重要意义. 近年来的研究表明［2-4］，BMI1信号通路不仅对于维持干细胞的生物特性具有重要意义，而且在肿瘤发生过程中也扮演重要的角色. 本研究试图从BMI1信号通路的重要组分，细胞周期素依赖性蛋白激酶（CDK6）入手，利用计算机辅助进行药物设计，以期找到可以对该信号通路进行干预的良好的先导化合物，为进一步药物筛选打下良好的基础.
1材料和方法
1.1CADD（计算机辅助药物设计）获得CDK6活性位点： 建立参数文件cdk6pocket. index，进入LibBuilder程序pocket子程序，运行pocket,得到CDK6的活性位点信息. 用C自行设计程序取得CDK6活性口袋的坐标信息，并查找配体p16ink4d位于活性口袋的原子坐标信息. 选取其中的原子N， CA， C，0结构信息作为种子结构,用Hyperchem正确加氢后存储为mo12文件.
获得配体： 根据Lipinski法则编辑参数文件cdk6grow. index,进入LigBuilder子程序grow,运行grow，选取50个生成结构，在Linux Red Hat9.0操作系统下，CPU为AMD Athlon1800+的PC上运行时间为1.5 h. 编辑process参数文件，运行process，得到结果文件.
1.2细胞株的培养急性白血病细胞系K562（本室保存）细胞在含100 mL/L热灭活小牛血清的RPMI 1640培养液中常规培养传代（37℃, 50 mL/L CO2）. 10种化合物分别溶于三蒸水中，配成5×10-3 mol/L储存液储存在-20℃，对照组培养中仅含相应量的DMSO（0.1%）.
1.3细胞毒试验采用常规MTT法测定药物体外抗肿瘤活性,简言之,取对数生长期细胞,用含有100 mL/L小牛血清的RPMI 1640培养液配成1×108/L，接种于96孔微培养板（细胞数为1×105， 180 μL/孔），37℃， 50 mL/L CO2条件下培养过夜. 将化合物对倍稀释成6个浓度，每孔加药液20 μL，以STI571为阳性对照，阴性对照加等体积生理盐水，使终体积为200 μL/孔. 药物作用72 h后，每孔加MTT 20 μL（5 g/L）,37℃继续培养4 h,离心(1000 r/min, 10 min)弃上清,每孔加入DMSO 150 μL,振荡至沉淀完全溶解；在酶标仪上检测432 nm处光密度A值. 细胞生长抑制率按以下公式计算：抑制率（%）=（对照组A值-加药组A值）/对照组A值×100%，以同一药物的不同浓度对肿瘤细胞生长抑制率作图可得到剂量反应曲线，根据回归方程求出该药物的半数抑制浓度(IC50)，即细胞存活率减少一半时的药物剂量.
2结果
生成50个候选化合物，除去对受体配体结合不利的原子共有5种骨架结构，合成得到8种结构进行生物活性测定，化合物（图1）对K562肿瘤细胞生长抑制实验结果表明，8个化合物中有3个具有抑制活性，IC50介于220~310 μmol/L之间（表1）. 化合物主要和CDK6活性口袋中的LEU166， ARG168上的疏水集团，ARG168， GLN149上的氢键供体氨基形成稳定结合（图2~4）.
图1CDK6抑制化合物母核结构式（略）
表1CDK6抑制剂对K562细胞系的抑制作用
3讨论
BMI1是维持干细胞生物学特性的重要蛋白,近年来研究发现该蛋白同样对于肿瘤的维持具有重要作用. BMI1发挥作用的途径称为BMI1信号通路，该通路有两个分支： 一条分支是通过抑制p19arf，最终抑制p53发挥作用. 该通路中有一个重要蛋白MDM2，由于该蛋白可以抑制p53蛋白发挥抗癌作用，利用CADD方法设计了MDM2小分子抑制剂，成功地在人类肿瘤长期移植的裸鼠身上抑制了肿瘤的生长；在该通路的一个分支中，BMI1通过抑制p16ink4a发挥作用［5］. P16ink4a也是一种肿瘤抑制蛋白，它可以抑制细胞周期素依赖性蛋白激酶（CDK6）和细胞周期素D(cyclin D)的结合，该蛋白受抑制时，CDK6可以和cyclin D结合激活，从而进一步磷酸化视网膜母细胞瘤肿瘤抑制蛋白（pRB）使E2F依赖性的转录过程发生，细胞内DNA合成，进入周期循环. 干细胞（肿瘤细胞）依赖此机制不停地进行细胞循环和DNA合成以维持自身. 根据这一特性，如果能设计一种可以和CDK6稳定结合的抑制剂发挥类似p16ink4a的作用，将对肿瘤产生稳定的抑制作用，本研究得到的50种化合物均可以和靶蛋白稳定结合（pkd>8.10）.
图2化合物1 C24H2803和CDK6相互作用图 pkd=8.1（略）
图3化合物2 C28H3803和CDK6相互作用图 pkd=9.2（略）
由于软件接受的受体蛋白信息是基于p16ink4a/CDK6的复合物，所以得到的与CDK6相结合的活性位点是和p16ink4a对应的. CADD已成为药物设计中的常规方法之一，利用CADD，新药开发的周期可以大大缩短，经费开支大大减少. 利用CADD已成功设计了多种新药，如陈凯先等［6］用CoMRA方法对喹诺酮类抗菌药物的定量构效关系（QSAR）的研究发现了3种生物活性强于氧氟沙星的新化合物；Gschwend等［7］用分子对接方法设计的二氢叶酸还原酶（DHFR）抑制剂. 本研究采用从头设计法（de novo design）不依赖于已知三维结构的先导化合物，得出的分子是全新的化合物，但同时也可能出现化合物合成困难的问题. 本研究得到的50种化合物只有8种可以合成.
图4化合物3 C25H3202和CDK6相互作用图 pkd=8.6（略）
LigBuilder是北京大学来鲁华教授课题组开发的基于分子结构的CADD程序［8-9］，可运行于除Windows外的多数机操作系统，包括Unix, Linux等. 利用该软件可采用从头设计的方法设计先导化合物，该软件结合了国外从头设计软件LUDI和grid的优点，使利用一种软件就可以完成包括种子结构的设计，片段生长（grow），片段连接（link）在内的各种工作成为可能. 该程序包还为使用者提供了极大的灵活性，如种子结构的生成就需要使用者灵活运用各种手段以得到合适的种子结构，而种子结构的选择对于化合物生成具有决定性的意义，一般是从配体位于靶蛋白的活性口袋内的原子中选取. 本研究自行设计了寻找种子结构的程序，但不提供优化功能. LigBuilder程序采用遗传算法最终得到的先导化合物是成千上万符合筛选条件的化合物中的佼佼者，因此可以最大程度地满足使用者的要求.
本研究利用CADD方法成功设计了CDK6抑制剂的先导化合物，并初步证明部分化合物具有一定的生物活性，以后的研究拟对先导化合物结构做进一步优化，期望最终得到具有良好抑制活性的CDK6抑制剂.
【】
［1］ Kopper L, Hajdu M. Tumor stem cells ［J］. Pathol Oncol Res, 2004,10:69-73.
［2］ Max J. Mutant stem cells may seed cancer ［J］. Science, 2003,301:1308-1310.
［3］ Muhammad AH, Michael WB, Max W. Therapeutic implications of cancer stem cells ［J］. Curr Opin Genet Dev, 2004,14:43-47.
［4］ Park IK, Morrison SJ, Clarke MF. Bmil, stem cells, and senescence regulation ［J］. J Clin Invest, 2004,113(2):175-179.
［5］ Vassilev LT, Vu BT, Graves B, et al. In vivo activation of the p53 pathway by smallmolecule antagonists of MDM2 ［J］. Science, 2004,303(5659):844-848.
［6］ 陈凯先,蒋华良,嵇汝运. 计算机辅助药物设计原理、方法及应用［M］. 上海： 上海技术出版社，2000：12-13.
［7］ Gschwend DA, Sirawaraporn W, Santi DV, et al. Specificity in structurebased drug design: Identification of a novel, selective inhibitor of Pneumocystis carinii dihydrofolate reductase ［J］. Proteins, 1997,29:59-67.
［8］ Wang RX, Gao Y, Lai LH. LigBuilder: A multipurpose program for structurebased drug design ［J］. J Mol Model, 2000,6:498-516.
［9］ Wang RX, Liu L, Lai LH, et al. SCORE: A new empirical method for estimating the binding affinity of a proteinligand complex ［J］. J Mol Model, 1998,4(12):379-394.
