-
针对PD-1/PD-L1免疫检查点的单克隆抗体抑制剂逐渐进入市场并在多种类型的肿瘤治疗中取得一定的积极效果. 然而, 随着应用范围的不断扩展, 抗体药物的局限性以及过多同质化研究等问题逐渐显现出来, 小分子化合物抑制剂成为了研究者们关注的新焦点. 本文旨在利用基于配体和基于结构的结合活性预测方法实现针对PD-L1靶点的小分子化合物虚拟筛选, 从而帮助加速小分子药物的开发. 通过从相关研究文献及专利收集PD-L1小分子抑制活性数据集, 根据不同分子表征方法和算法构建机器学习活性判定分类模型和活性强度预测回归模型, 两类模型从大型类药小分子库(ZINC15)中筛选获得68种高PD-L1抑制活性候选化合物. 其中10种化合物不仅具备良好的药物相似性和药代动力学, 还在分子对接中与已报道的热点化合物表现出同等水平的结合强度和相似的作用机制, 这一现象在后续分子动力学模拟和结合自由能估计中得到进一步验证. 本文提出了一个融合基于配体方法和基于结构方法的计算机辅助药物研发工作流程, 其在大型化合物数据库中有效筛选出有潜力的PD-L1小分子抑制剂, 有望助力加速肿瘤免疫治疗的应用.
-
关键词:
- PD-1/PD-L1/
- 虚拟筛选/
- 机器学习/
- 分子动力学模拟
Monoclonal antibody inhibitors targeting PD-1/PD-L1 immune checkpoints are gradually entering the market and have achieved certain positive effects in the treatments of various types of tumors. However, with the expansion of application, the limitations of antibody drugs and problems such as excessive homogenization of research gradually appear, making small-molecule inhibitors the new focus of researchers. This study aims to use ligand-based and structure-based binding activity prediction methods to achieve virtual screening of small-molecule inhibitors targeting PD-L1, thereby helping to accelerate the development of small molecule drugs. A dataset of PD-L1 small-molecule inhibitory activity from relevant research literature and patents is collected and activity judgment classification models with intensity prediction regression models are constructed based on different molecular featurization methods and machine learning algorithms. The two types of models filter 68 candidate compounds with high PD-L1 inhibitory activity from a large drug-like small molecule screening pool (ZINC15). Ten of these compounds not only have good drug similarities and pharmacokinetics, but also exhibit comparable binding affinities and similar mechanisms of action with previous reported hotspot compounds in molecular docking. This phenomenon is further verified in subsequent molecular dynamics simulation and the estimation of binding free energy. In this study, a virtual screening workflow integrating ligand-based method and structure-based method is developed, and potential PD-L1 small-molecule inhibitors are effectively screened from large compound databases, which is expected to help accelerate the application and expansion of tumor immunotherapy.[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] -
序号 化合物 化学结构式 对接分数/(kcal·mol–1) Hit 1 ZINC000021723762 –11.8 Hit 2 ZINC000021874692 –10.7 Hit 3 ZINC000175468610 –10.7 Hit 4 ZINC000019770413 –10.6 Hit 5 ZINC000021874694 –10.6 Hit 6 ZINC000952973550 –10.5 Hit 7 ZINC000064987401 –10.5 Hit 8 ZINC000003908573 –10.4 Hit 9 ZINC000020538424 –10.4 Hit 10 ZINC000004063088 –10.3 Ctr 1 BMS202 –11.1 Ctr 2 INCB086550 –10.0 化合物 $ {E}_{{\mathrm{v}}{\mathrm{d}}{\mathrm{w}}} $/(kcal·mol–1) $ {E}_{{\mathrm{e}}{\mathrm{l}}{\mathrm{e}}} $/(kcal·mol–1) $ {E}_{{\mathrm{P}}{\mathrm{B}}} $/(kcal·mol–1) $ {E}_{{\mathrm{S}}{\mathrm{A}}} $/(kcal·mol–1) $ \Delta G $/(kcal·mol–1) ZINC000021723762 –58.86$ \pm 0.12 $ –8.76$ \pm $0.10 35.53$ \pm $0.09 –4.03$ \pm $0.00 –36.12$ \pm $0.12 ZINC000021874692 –50.42$ \pm $0.11 –6.50$ \pm $0.14 35.65$ \pm $0.16 –4.50$ \pm $0.01 –25.78$ \pm 0.12 $ ZINC000175468610 –42.46$ \pm 0.09 $ –14.34$ \pm 0.08 $ 36.20$ \pm 0.13 $ –4.02$ \pm $0.00 –24.61$ \pm 0.10 $ ZINC000019770413 –47.31$ \pm 0.08 $ –3.45$ \pm 0.08 $ 26.32$ \pm 0.07 $ –3.48$ \pm $0.00 –27.91$ \pm 0.08 $ ZINC000021874694 –55.17$ \pm 0.09 $ –17.87$ \pm 0.13 $ 45.67$ \pm 0.14 $ –4.59$ \pm $0.00 –31.97$ \pm 0.12 $ ZINC000952973550 –55.82$ \pm 0.08 $ –13.47$ \pm 0.11 $ 45.83$ \pm 0.11 $ –4.19$ \pm $0.00 –27.65$ \pm 0.11 $ ZINC000064987401 –48.50$ \pm 0.08 $ –14.71$ \pm 0.12 $ 38.59$ \pm 0.11 $ –4.06$ \pm $0.00 –28.68$ \pm 0.12 $ ZINC000003908573 –44.75$ \pm 0.07 $ 4.59$ \pm 0.11 $ 20.20$ \pm 0.13 $ –3.74$ \pm $0.00 –23.70$ \pm 0.12 $ ZINC000020538424 –49.50$ \pm 0.09 $ –8.50$ \pm 0.13 $ 32.61$ \pm 0.16 $ –4.17$ \pm $0.01 –29.57$ \pm 0.10 $ ZINC000004063088 –64.14$ \pm 0.09 $ –5.38$ \pm 0.11 $ 36.75$ \pm 0.10 $ –4.21$ \pm $0.00 –36.98$ \pm 0.10 $ BMS202 –40.23$ \pm 0.08 $ –3.57$ \pm 0.07 $ 22.35$ \pm 0.10 $ –4.31$ \pm $0.01 –25.75$ \pm 0.09 $ INCB086550 –50.65$ \pm 0.14 $ –9.87$ \pm 0.17 $ 44.44$ \pm 0.27 $ –5.13$ \pm $0.02 –21.22$ \pm $0.15 -
[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67]
计量
- 文章访问数:2601
- PDF下载量:137
- 被引次数:0