CN108399315B - 一种Bcr-Abl蛋白激酶抑制剂的筛选方法 - Google Patents

一种Bcr-Abl蛋白激酶抑制剂的筛选方法 Download PDF

Info

Publication number
CN108399315B
CN108399315B CN201810173005.2A CN201810173005A CN108399315B CN 108399315 B CN108399315 B CN 108399315B CN 201810173005 A CN201810173005 A CN 201810173005A CN 108399315 B CN108399315 B CN 108399315B
Authority
CN
China
Prior art keywords
formula
specificity
protein kinase
bcr
docking
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810173005.2A
Other languages
English (en)
Other versions
CN108399315A (zh
Inventor
晏致强
汪劲
夏艳
赵银苹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun Institute of Applied Chemistry of CAS
Original Assignee
Changchun Institute of Applied Chemistry of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changchun Institute of Applied Chemistry of CAS filed Critical Changchun Institute of Applied Chemistry of CAS
Priority to CN201810173005.2A priority Critical patent/CN108399315B/zh
Publication of CN108399315A publication Critical patent/CN108399315A/zh
Application granted granted Critical
Publication of CN108399315B publication Critical patent/CN108399315B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/50—Molecular design, e.g. of drugs

Abstract

本发明提供了一种Bcr‑Abl蛋白激酶抑制剂的筛选方法，首先通过谷本相似性搜索得到的与已知Bcr‑Abl蛋白激酶抑制剂具有高相似度的小分子和对接靶标进行对接，然后采用打分函数SPA对蛋白质‑小分子对接构象进行分子间相互作用打分，得到亲和性值E；并将得到的亲和性值E，结合热力学特异性和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值，最后得到具有式(I)结构的Bcr‑Abl蛋白激酶抑制剂；结果表明筛选方法得到的抑制剂对CML细胞有很好的抑制作用，在很大浓度范围内对正常细胞无毒性，而且具有跟普纳替尼相媲美的抑制CML变异体肿瘤的作用。

Description

一种Bcr-Abl蛋白激酶抑制剂的筛选方法
技术领域
本发明涉及药物筛选领域，尤其涉及一种Bcr-Abl蛋白激酶抑制剂的筛选方法。
背景技术
Bcr-Abl酪氨酸激酶是蛋白激酶的一种，它是由bcr-abl致癌基因所表达。Bcr-Abl致癌基因是由染色体9上的abelson(abl)基因和染色体22上的break point cluster(bcr)基因相互重排而形成。Bcr-Abl酪氨酸激酶的持续活性是导致慢性粒细胞白血病(chronicmyeloid leukaemia，CML)的原因。增强的激酶活性承担细胞内的磷酸化和多种信号分子的激活，这些功能加速了细胞的增殖同时抑制了细胞的凋亡。因此，对Bcr-Abl酪氨酸激酶活性的抑制是靶向治疗CML的有效途径。
伊马替尼是2001年FDA首次批准治疗CML的药物，也是第一个激酶抑制剂。Bcr-Abl酪氨酸激酶是典型的会对CML病人产生耐药性的激酶。超过100个位点突变被发现对一代药物伊马替尼(imatinib)的治疗具有耐药性。由于耐药性的产生，使得针对Bcr-Abl酪氨酸激酶药物开发进入后伊马替尼时代(post-imatinib era)以对抗耐药性。在后伊马替尼时代，两种第二代药物尼罗替尼(nilotinib)和达沙替尼(dasatinib)具有更强的药效，更重要的是，这两种药物能够对Bcr-Abl耐药性突变体保持抑制作用。然而，这两种药物都不能抑制耐药T315I激酶突变体。这使得第三代药物普纳替尼(ponatinib)的开发与产生，并在2012年12月获得美国FDA批准上市。普纳替尼是第一种泛Bcr-Abl激酶抑制剂(pan-Bcr-Ablkinase inhibitor)，它能抑制所有已知的抗伊马替尼的突变体，包括T315I突变体。然而，近些年研究表明，普纳替尼具有严重的副作用，因此，临床上急需比普纳替尼具有更好药效和更高特异性的药物。
尽管针对Abl-Bcr酪氨酸激酶的药物开发取得了很大的进展，例如第三代药物普纳替尼的开发，但新的耐药性位点突变的可能性依然存在。这促使必须持续地开发更加有效的抑制剂。同时，每一种抑制剂会产生它自身的耐药性突变谱，CML病人同时使用多个抑制剂能够获得更好的疗效。另外，目前治疗CML的靶向药物价格昂贵，这就需要我们能够设计出比普纳替尼具有更好的药效和抗耐药性并且价格较低的替代药物。因此，开发全新的抑制剂以改善和提高现有的疗效对于CML的治疗依然非常紧迫，尤其是开发能够抑制T315I突变体的第三代泛Abl-Bcr激酶抑制剂。
发现能够与靶蛋白有潜在识别能力和相互作用的先导化合物(lead compound)是前期药物开发(drug discovery)的关键，研究和理解蛋白质与小分子的识别和相互作用机制对新药开发与设计有着直接的应用。小分子高效且特异地结合其靶蛋白取决于两个方面，一是结合亲和力(affinity)，二是结合特异性(specificity)。亲和力是指蛋白质-小分子结合复合物的稳定性，特异性是指小分子偏好于与靶蛋白结合而不与其他蛋白结合。高的亲和力并不能保证高的特异性，而特异性是避免药物分子副作用的关键。在筛选和设计一个新药时，预期的效果是新药能有效地与靶蛋白结合，同时又尽量少地与其它蛋白发生作用，从而避免副作用的产生。因此量化蛋白质与小分子结合的特异性对新药的开发与设计至关重要。目前，对蛋白质-小分子相互作用的理论计算与实验测定主要是量化其亲和力的大小，而忽略了对特异性的考虑。这是因为特异性的量化非常困难。特异性包括热力学特异性(thermodynamic specificity)和动力学特异性(kinetic specificity)。热力学特异性的量化需要比较一个小分子与所有蛋白之间的结合亲和力，这在现有的数据信息和实验条件下，是难以实现的。动力学特异性是指小分子在蛋白质上的停留时间(residencetime，RT)。虽然实验技术可以批量测定小分子停留时间以量化动力学特异性，但在资金和小分子测定数量上还是有限。在理论计算上，利用分子动力学模拟计算大量的蛋白质-小分子复合物的动力学停留时间仍然非常具有挑战性。因此，提供更为有效的筛选Abl-Bcr激酶抑制剂的方法以期节约成本，且得到活性较高的药物是目前需要解决的问题。
发明内容
有鉴于此，本发明所要解决的技术问题在于提供一种Bcr-Abl蛋白激酶抑制剂的筛选方法，本发明提供的筛选方法不仅筛选效率高，而且得到的抑制剂的活性高。
本发明提供了一种Bcr-Abl蛋白激酶抑制剂的筛选方法，包括：
1)将通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子和对接靶标进行对接，然后采用打分函数SPA对蛋白质-小分子对接构象进行分子间相互作用打分，得到靶蛋白与小分子间的总能量，即亲和性值E；
2)利用步骤1)得到的亲和性值E，结合热力学特异性公式和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值，
其中，所述热力学特性公式为：
Figure BDA0001585997560000031
δE为最低能量与平均能量差，即δE＝EN-<E>，
ΔE为能量涨落，即(<E2>-<E>2)1/2；
所述动力学特性公式为：
Figure BDA0001585997560000032
F(r)为不同结合位置r上的亲和性值，
KB，T以及D(r)设为常数；
3)对得到的小分子的亲和性值E、热力学特异性ISR值和动力学特异性RT值进行归一化处理得到小分子的PE，PISR，PRT，然后对三者加和得到小分子多维参数S，将得到的小分子多维参数S与普纳替尼的多维参数进行比较，得到具有式(I)结构的Bcr-Abl蛋白激酶抑制剂，
Figure BDA0001585997560000041
其中，R1为式(R1-1)、式(R1-2)或式(R1-3)，
Figure BDA0001585997560000042
R2为C5～C8的含氮芳香基。
优选的，所述通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子的方法具体为：
利用谷本相似性搜索对PubChem数据库里大约40,000,000个小分子进行初步筛选，选出与伊马替尼，尼罗替尼，达沙替尼和普纳替尼中的任意一个结构相似度大于谷本系数90％库中小分子，即为与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子。
优选的，所述对接靶标为活化野生型蛋白激酶，活化T315I突变体蛋白激酶，非活化野生型蛋白激酶和非活化T315I突变体蛋白激酶。
优选的，所述对接采用的软件为AutoDock4.2。
优选的，所述R2为式(R2-1)、或式(R2-2)，
Figure BDA0001585997560000043
优选的，所述Bcr-Abl蛋白激酶抑制剂为式(I-1)、式(I-2)或式(I-3)，
Figure BDA0001585997560000044
与现有技术相比，本发明提供了一种Bcr-Abl蛋白激酶抑制剂的筛选方法，首先将通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子和对接靶标进行对接，然后采用打分函数SPA对蛋白质-小分子对接构象进行分子间相互作用打分，得到靶蛋白与小分子间的总能量，即亲和性值E；并将得到的亲和性值E，结合热力学特异性和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值，并对上述值进行归一化处理得到小分子的PE，PISR，PRT，然后对三者加和得到小分子多维参数S，将得到的小分子多维参数与普纳替尼的多维参数进行比较，得到具有式(I)结构的Bcr-Abl蛋白激酶抑制剂；实验结果表明，本发明提供的筛选方法得到的抑制剂对CML细胞有很好的抑制作用，在很大浓度范围内对正常细胞无毒性，而且具有跟普纳替尼相媲美的抑制CML变异体肿瘤的作用。
附图说明
图1为非活化态和活化态两种构象的Bcr-Abl蛋白的结构示意图；
图2为不同浓度的PA对三种白血病细胞和两种正常细胞的增殖抑制作用结果；
图3为不同浓度的PA诱导三种白血病细胞凋亡结果；
图4为不同浓度的PA对三种白血病细胞周期的影响；
图5为PA对细胞线粒体膜电位的影响；
图6为检测不同浓度PA处理24h后三种白血病细胞内活性氧(ROS)水平结果；
图7为不同浓度的PA对三种白血病细胞中融合蛋白Bcr-Abl以及其磷酸化蛋白表达水平的影响；
图8为PA有效抑制动物体内移植肿瘤细胞的生长的结果。
具体实施方式
本发明提供了一种Bcr-Abl蛋白激酶抑制剂的筛选方法，包括：
1)将通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子与对接靶标进行对接，然后采用打分函数SPA对蛋白质-小分子对接构象进行分子间相互作用打分，得到靶蛋白与小分子间的总能量，即亲和性值E；
2)利用步骤1)得到的亲和性值E，结合热力学特异性和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值，
其中，所述热力学特性公式为：
Figure BDA0001585997560000061
δE为最低能量与平均能量差，即δE＝EN-<E>，
ΔE为能量涨落，即(<E2>-<E>2)1/2；
所述动力学特性公式为：
Figure BDA0001585997560000062
F(r)为不同结合位置r上的亲和性值，
KB，T以及D(r)设为常数；
3)对得到的小分子的亲和性值E、热力学特异性ISR值和动力学特异性RT值进行归一化处理得到小分子的PE，PISR，PRT，然后对三者加和得到小分子多维参数S，将得到的小分子多维参数与普纳替尼的多维参数进行比较，得到具有式(I)结构的Bcr-Abl蛋白激酶抑制剂，
Figure BDA0001585997560000063
其中，R1为式(R1-1)、式(R1-2)或式(R1-3)，
Figure BDA0001585997560000064
R2为C5～C8的含氮芳香基。
按照本发明，本发明将通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子与对接靶标进行对接，然后采用打分函数SPA对蛋白质-小分子对接构象进行分子间相互作用打分，得到靶蛋白与小分子间的总能量，即亲和性值E；其中，所述对接靶标优选为活化野生型蛋白激酶(PDB ID 2GQG)，活化T315I突变体蛋白激酶(PDBID 2V7A)，非活化野生型蛋白激酶(PDB ID 1E1P)和非活化T315I突变体蛋白激酶(PDB ID3IK3)；所述述对接采用的软件优选为AutoDock4.2；本发明对对蛋白质-小分子对接构象进行分子间相互作用打分的打分函数没有特殊要求，本领域公知的打分函数均可，如可以按照“Yan，Z.and Wang，J.(2012)Specificity quantification of biomolecularrecognition and its implication for drug discovery.Sci.Rep.，2，1--7”的文章中公开的打分函数SPA进行打分。
本发明中，所述通过谷本相似性搜索得到的与现有药物分子具有高相似度的小分子的方法优选为：
利用谷本相似性搜索对PubChem数据库里大约40,000,000个小分子进行初步筛选，选出与伊马替尼，尼罗替尼，达沙替尼和普纳替尼中的任意一个结构相似度大于谷本系数90％库中小分子，即为与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子。
按照本发明，利用步骤1)得到的亲和性值E，结合热力学特异性和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值，
其中，所述热力学特性公式为：
Figure BDA0001585997560000071
δE为最低能量与平均能量差，即δE＝EN-<E>，
ΔE为能量涨落，即(<E2>-<E>2)1/2；
所述动力学特性公式为：
Figure BDA0001585997560000072
F(r)为不同结合位置r上的亲和性值，也即结合能量值。
在计算时，RT公式采用离散化计算求和，KB，T以及D(r)设为常数。
按照本发明，本发明对得到的小分子的亲和性值E、热力学特异性ISR值和动力学特异性RT值进行归一化处理得到小分子的PE，PISR，PRT，然后对三者加和得到小分子多维参数S，将得到的小分子多维参数与普纳替尼的多维参数进行比较，得到具有式(I)结构的Bcr-Abl蛋白激酶抑制剂；其中，本发明对归一化处理的方法没有特殊要求，本领域公知的归一化处理方法均可；所述式(I)结构的Bcr-Abl蛋白激酶抑制剂中，所述R2优选为式(R2-1)、或式(R2-2)，
Figure BDA0001585997560000081
更具体的，所述Bcr-Abl蛋白激酶抑制剂为式(I-1)、式(I-2)或式(I-3)，
Figure BDA0001585997560000082
本发明提供了一种Bcr-Abl蛋白激酶抑制剂的筛选方法，首先将通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子与对接靶标进行对接，然后采用打分函数SPA对蛋白质-小分子对接构象进行分子间相互作用打分，得到靶蛋白与小分子间的总能量，即亲和性值E；并将得到的亲和性值E，结合热力学特异性和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值，并对上述值进行归一化处理得到小分子的PE，PISR，PRT，然后对三者加和得到小分子多维参数S，将得到的小分子多维参数与普纳替尼的多维参数进行比较，得到具有式(I)结构的Bcr-Abl蛋白激酶抑制剂；使得本发明得到的抑制剂对CML细胞有很好的抑制作用，而且在很大浓度范围内对正常细胞无毒性，具有跟普纳替尼相媲美的抑制CML变异体肿瘤的作用。
下面将结合本发明实施例的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例
1)相似度筛选
通过利用谷本相似性搜索(Tanimoto similarity search)对PubChem数据库里大约40,000,000个小分子进行初步筛选。库中小分子与伊马替尼，尼罗替尼，达沙替尼和普纳替尼中任何一个的结构相似度大于谷本系数90％的被保留进入下一步的分子对接。通过这步筛选，得到数量相对较少但结构与药物分子非常相似的一组小分子化合物，结果见表1，另外，选取美国/国家癌症机构(National Cancer Institute，NCI)多样性小分子组(diversity set)作为对比组。
表1.pubchem数据库中与现有药物分子具有高相似度的小分子化合物的个数
Figure BDA0001585997560000091
2)对接与能量计算
以上得到的小分子分别与四个对接靶标(活化野生型(PDB ID 2GQG)，活化T315I突变体(PDB ID 2V7A)，非活化野生型(PDB ID 1E1P)和非活化T315I突变体(PDB ID3IK3))进行分子对接，采用的对接软件是AutoDock4.2。对接完成后，利用Yan，Z.and Wang，J.(2012)Specificity quantification of biomolecular recognition and itsimplication for drug discovery.Sci.Rep.，2，1-7文章中报道的打分函数SPA重新对蛋白质-小分子对接构象进行分子间相互作用打分，所得到的能量值用于热力学特异性和动力学特异性的计算。SPA是一套原子间相互作用能量值，利用SPA的值可计算靶蛋白与小分子间的总能量，即亲和性。其中，所述非活化态和活化态两种构象的Bcr-Abl蛋白的结构示意图见图1，图1为非活化态和活化态两种构象的Bcr-Abl蛋白的结构示意图。
3)特异性计算
利用得到的对接分子构象和SPA打分所得到的亲和性值E，结合热力学特异性和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值。ISR计算公式由漏斗能量地貌理论推导所得。RT计算公式由扩散方程推导所得。两者公式分别如下：
Figure BDA0001585997560000101
Figure BDA0001585997560000102
δE为最低能量与平均能量差，即δE＝EN-<E>。ΔE为能量涨落，即(<E2>-<E>2)1/2。F(r)为不同结合位置r上的亲和性值，也即结合能量值。在计算时，RT公式采用离散化计算求和，KB，T以及D(r)设为常数。
4)多维筛选
通过以上步骤，得到每个分子与对接靶标的亲和性值E，热力学特异性ISR值和动力学特异性RT值。由于每个参数的量纲以及所代表的意义不一样，分别对每个量纲进行归一化处理得到PE，PISR，PRT，然后对三者加和得到参数S，即多维参数。利用多维参数S的值，对所有小分子进行排序。通过比较，得到具有式(I-1)、式(I-2)或式(I-3)的Bcr-Abl蛋白激酶抑制剂。这三个小分子化合物可作为先导化合物用于前期实验测试和可能的临床测试。以具有式(I-3)的小分子化合物(命名为PA，CID 24826974)，其多维参数值明显优于现有药物普纳替尼，结果见表2，此小分子化合物作为先导化合物用于一系列的实验测试。
Figure BDA0001585997560000103
表2先导化合物(PA)与药物普纳替尼的比较。
Figure BDA0001585997560000111
5)实验验证
针对筛选出的具有式(I-3)结构的先导小分子化合物普纳替尼衍生物PA，我们科研组通过一系列体外和体内实验对这些小分子化合物进行生化测试，包括增殖实验(proliferation assay)即MTS方法测定先导化合物对癌细胞系以及正常细胞系的增殖影响；利用DAPI细胞染色、流式细胞实验测定先导化合物对三种CML细胞凋亡的影响；进一步测定细胞内线粒体膜电位、活性氧ROS以及蛋白表达确定先导化合物PA引起细胞凋亡的机制；通过肿瘤移植动物实验观测先导化合物对肿瘤的抑制作用。
具体的，通过对三种CML细胞以及正常细胞的毒性实验，且所有实验都是利用MTS方法酶标仪测试；结果见图2，图2为不同浓度的PA对三种白血病细胞和两种正常细胞的增殖抑制作用结果(n＝3)；其中，图中a为PA化合物的结构和对三种白血病细胞的毒性实验；b为PA对小鼠成纤维细胞(L929)和人正常结肠上皮细胞(NCM460)的毒性。从图中可以看出，PA化合物对野生型Bcr-Abl和变异体T315I表达的CML细胞都具有纳摩级(nM)的IC50值，而且抑制作用随着浓度的增加呈浓度依赖关系，且和普纳替尼相比PA对正常细胞没有毒性；
采用不同浓度的PA诱导三种白血病细胞凋亡，通过DAPI染色对结果进行处理，结果见图3，图3为不同浓度的PA诱导三种白血病细胞凋亡结果，其中，(a)PA处理24h后荧光显微镜下细胞形态变化(比例尺：50μm)；(b)Annexin V-FITC/PI双染流式细胞术检测细胞凋亡率。从图中可以看出，对照组细胞的细胞核为圆形，染色质密度均匀一致。用药组细胞的细胞核边缘不规则，染色体浓集，着色较重，并伴有细胞核固缩。因此，从形态学上说明PA能够诱导细胞发生凋亡。
通过流式细胞术检测，结果见图4，图4为不同浓度的PA对三种白血病细胞周期的影响；从图中可以看出，流式细胞实验检测结构显示与对照组相比，随着浓度的增加，细胞凋亡率显著提高。
凋亡机制可能与线粒体膜电位和活性氧有关，线粒体膜电位下降是细胞凋亡早期的标志性事件，通过检测JC-1单体和JC-1聚合物，使用荧光显微镜观察JC-1从红色荧光到绿色荧光的转变来判断细胞是否发生凋亡。结果见图5，图5为PA对细胞线粒体膜电位的影响；其中，JC-1是一种用于检测线粒体膜电位ΔΨm的荧光探针。(a)JC-1单体(绿色荧光)；(b)JC-1二聚体(红色荧光)；(c)混合；从图5可以看出，随着PA的加入，细胞中红色荧光减弱，绿色荧光增强，说明线粒体膜电位下降，细胞发生凋亡。
通过流式细胞仪测定DCFH-DA荧光染料的荧光强度变化定量检测细胞内ROS水平，结果见图6，图6为检测不同浓度PA处理24h后三种白血病细胞内活性氧(ROS)水平结果。
通过蛋白印迹法研究细胞凋亡和融合蛋白Bcr-Abl的表达水平之间的关系，结果见图7，图7为不同浓度的PA对三种白血病细胞中融合蛋白Bcr-Abl以及其磷酸化蛋白表达水平的影响，从图中可以看出，随着PA浓度的增加，蛋白表达降低。同时移植的动物(小鼠)肿瘤随着给药时间的增加，肿瘤大小明显减小，说明此先导化合物对变异体肿瘤的治疗效果良好。
H&E免疫组织染色结果见图8，图8为PA有效抑制动物体内移植肿瘤细胞的生长的结果，其中，(a-b)PA对变异体肿瘤生长的抑制作用；(c)用药后肿瘤体积的变化；(d)用药后老鼠体重的变化；(e)H&E对五脏组织免疫染色和抗体肿瘤组织染色。从图中可以看出，PA对五脏没有毒性，而对c-abl蛋白具有很好的抑制作用。
以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。

Claims (4)

1.一种Bcr-Abl蛋白激酶抑制剂的筛选方法，包括：
1)将通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子和对接靶标进行对接，然后采用打分函数SPA对蛋白质-小分子对接构象进行分子间相互作用打分，得到靶蛋白与小分子间的总能量，即亲和性值E；
所述通过谷本相似性搜索得到的与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子的方法具体为：
利用谷本相似性搜索对PubChem数据库里40,000,000个小分子进行初步筛选，选出与伊马替尼,尼罗替尼，达沙替尼和普纳替尼中的任意一个结构相似度大于谷本系数90％的PubChem数据库中的小分子，即为与已知Bcr-Abl蛋白激酶抑制剂具有高相似度的小分子；
所述对接靶标为活化野生型蛋白激酶，活化T315I突变体蛋白激酶，非活化野生型蛋白激酶和非活化T315I突变体蛋白激酶；
2)利用步骤1)得到的亲和性值E，结合热力学特异性公式和动力学特异性公式，计算出每个小分子与对接靶标的热力学特异性ISR值和动力学特异性RT值，
其中，所述热力学特性公式为：
Figure FDA0003438742740000011
δE为最低能量与平均能量差，即δE＝EN-<E>，
ΔE为能量涨落，即(<E2>-<E>2)1/2；
所述动力学特性公式为：
Figure FDA0003438742740000012
F(r)为不同结合位置r上的亲和性值，
KB,T以及D(r)设为常数；
3)对得到的小分子的亲和性值E、热力学特异性ISR值和动力学特异性RT值进行归一化处理得到小分子的PE,PISR,PRT，然后对三者加和得到小分子多维参数S，将得到的小分子多维参数S与普纳替尼的多维参数进行比较，得到具有式(I)结构的Bcr-Abl蛋白激酶抑制剂，
Figure FDA0003438742740000021
其中，R1为式(R1-1)、式(R1-2)或式(R1-3)，
Figure FDA0003438742740000022
R2为C5～C8的含氮芳香基。
2.根据权利要求1所述的筛选方法，其特征在于，所述对接采用的软件为AutoDock4.2。
3.根据权利要求1所述的筛选方法，其特征在于，所述R2为式(R2-1)、或式(R2-2)，
Figure FDA0003438742740000023
4.根据权利要求1所述的筛选方法，其特征在于，所述Bcr-Abl蛋白激酶抑制剂为式(I-1)、式(I-2)或式(I-3)，
Figure FDA0003438742740000024
CN201810173005.2A 2018-03-01 2018-03-01 一种Bcr-Abl蛋白激酶抑制剂的筛选方法 Active CN108399315B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810173005.2A CN108399315B (zh) 2018-03-01 2018-03-01 一种Bcr-Abl蛋白激酶抑制剂的筛选方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810173005.2A CN108399315B (zh) 2018-03-01 2018-03-01 一种Bcr-Abl蛋白激酶抑制剂的筛选方法

Publications (2)

Publication Number Publication Date
CN108399315A CN108399315A (zh) 2018-08-14
CN108399315B true CN108399315B (zh) 2022-02-22

Family

ID=63091574

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810173005.2A Active CN108399315B (zh) 2018-03-01 2018-03-01 一种Bcr-Abl蛋白激酶抑制剂的筛选方法

Country Status (1)

Country Link
CN (1) CN108399315B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112053742A (zh) * 2020-07-23 2020-12-08 中南大学湘雅医院 分子靶标蛋白的筛选方法、装置、计算机设备和存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101381320A (zh) * 2007-09-04 2009-03-11 西南大学 一类具有曼尼希碱结构的化合物及在制备治疗和/或预防白血病药物上的应用
CN102798708A (zh) * 2012-08-23 2012-11-28 中国科学院长春应用化学研究所 一种检测配体与靶标结合特异性的方法及药物筛选方法
CN103063657A (zh) * 2012-12-24 2013-04-24 中国科学院长春应用化学研究所 基于化学发光法检测蛋白酶活性及其抑制剂活性的方法
CN103582478A (zh) * 2011-04-07 2014-02-12 阿里亚德医药股份有限公司 治疗神经退行性疾病的方法和组合物

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2007380B1 (en) * 2006-03-30 2011-09-14 Nerviano Medical Sciences S.r.l. Use of a kinase inhibitor for the treatment of particular resistant tumors

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101381320A (zh) * 2007-09-04 2009-03-11 西南大学 一类具有曼尼希碱结构的化合物及在制备治疗和/或预防白血病药物上的应用
CN103582478A (zh) * 2011-04-07 2014-02-12 阿里亚德医药股份有限公司 治疗神经退行性疾病的方法和组合物
CN102798708A (zh) * 2012-08-23 2012-11-28 中国科学院长春应用化学研究所 一种检测配体与靶标结合特异性的方法及药物筛选方法
CN103063657A (zh) * 2012-12-24 2013-04-24 中国科学院长春应用化学研究所 基于化学发光法检测蛋白酶活性及其抑制剂活性的方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
BCR-ABL 激酶抑制剂的研究进展;李琳毅等;《中国现代应用药学》;20160731;第33卷(第7期);第954-962页 *
Ponatinib抑制野生型和T315I突变型BCR-ABL1激酶的分子模拟;陈清等;《化工进展》;20180228;第37卷(第02期);第702-707页 *
SPA-LN: a scoring function of ligand–nucleic acid interactions via optimizing both specificity and affinity;Zhiqiang Yan等;《Nucleic Acids Research》;20170707;第45卷(第12期);第1-8页 *
Substrate Activity Screening (SAS) and Related Approaches in Medicinal Chemistry;Rafaela Gladysz等;《ChemMedChem》;20160204;第11卷(第5期);第467-476页 *
Thermodynamic and kinetic specificities of ligand binding;Zhiqiang Yan等;《Chemical Science》;20130430;第2387-2395页 *
酪氨酸激酶抑制剂筛选方法的建立及姜黄素衍生物的筛选;许娇红;《中国优秀硕士论文全文数据库医药卫生科技辑》;20110215;E079-51 *

Also Published As

Publication number Publication date
CN108399315A (zh) 2018-08-14

Similar Documents

Publication Publication Date Title
Rzymski et al. SEL120-34A is a novel CDK8 inhibitor active in AML cells with high levels of serine phosphorylation of STAT1 and STAT5 transactivation domains
Kim et al. Rapamycin-insensitive companion of mTOR (RICTOR) amplification defines a subset of advanced gastric cancer and is sensitive to AZD2014-mediated mTORC1/2 inhibition
Lawal et al. mTOR/EGFR/iNOS/MAP2K1/FGFR/TGFB1 are druggable candidates for N-(2, 4-difluorophenyl)-2′, 4′-difluoro-4-hydroxybiphenyl-3-carboxamide (NSC765598), with consequent anticancer implications
Su et al. Targeting of multiple oncogenic signaling pathways by Hsp90 inhibitor alone or in combination with berberine for treatment of colorectal cancer
Conde et al. Allosteric antagonist modulation of TRPV2 by piperlongumine impairs glioblastoma progression
CN102216775A (zh) 对hsp90-抑制剂的易感性
Yan et al. Fascaplysin exert anti-tumor effects through apoptotic and anti-angiogenesis pathways in sarcoma mice model
Kim et al. Identification of novel synergistic targets for rational drug combinations with PI3 kinase inhibitors using siRNA synthetic lethality screening against GBM
Liang et al. CD271+ cells are diagnostic and prognostic and exhibit elevated MAPK activity in SHH medulloblastoma
Rathi et al. Therapeutic targeting of PIM KINASE signaling in cancer therapy: Structural and clinical prospects
Zhu et al. Discovery of novel selective PI3Kγ inhibitors through combining machine learning-based virtual screening with multiple protein structures and bio-evaluation
CN108399315B (zh) 一种Bcr-Abl蛋白激酶抑制剂的筛选方法
Meng et al. Raddeanin A induces apoptosis and cycle arrest in human HCT116 cells through PI3K/AKT pathway regulation in vitro and in vivo
Li et al. Upregulation of CCNB2 and its perspective mechanisms in cerebral ischemic stroke and all subtypes of lung cancer: a comprehensive study
Bhole et al. Unlocking the potential of PROTACs: A comprehensive review of protein degradation strategies in disease therapy
Zhang et al. Molecular docking and in vitro experiments verified that kaempferol induced apoptosis and inhibited human HepG2 cell proliferation by targeting BAX, CDK1, and JUN
Ahmed et al. Dual inhibition of glycolysis and glutaminolysis for synergistic therapy of rheumatoid arthritis
WO2014056627A1 (en) Novel approaches for individualized therapy of pancreatic ductal adenocarcinoma
Somadder et al. Drug repurposing and systems biology approaches of Enzastaurin can target potential biomarkers and critical pathways in colorectal Cancer
Kelly et al. Target-agnostic p-glycoprotein assessment yields strategies to evade efflux, leading to a BRAF inhibitor with intracranial efficacy
Gu et al. Discovery of 2-(5-(quinolin-6-yl)-1, 3, 4-oxadiazol-2-yl) acetamide derivatives as novel PI3Kα inhibitors via docking-based virtual screening
Zhang et al. Discovery of novel antagonists of human neurotensin receptor 1 on the basis of ligand and protein structure
Daniel et al. Anticancer potential of mebendazole against chronic myeloid leukemia: In silico and in vitro studies revealed new insights about the mechanism of action
Nair et al. Application of comprehensive bioinformatics approaches to reconnoiter crucial genes and pathways underpinning hepatocellular carcinoma: a drug repurposing endeavor
Qayoom et al. Decoding the molecular mechanism of stypoldione against breast cancer through network pharmacology and experimental validation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant