基于特征的深度学习预测化合物-蛋白质相互作用的研究进展
Research progress of feature-based deep learning for predicting compound-protein interaction
-
摘要: 药物研发过程中，化合物-蛋白质相互作用（compound-protein interaction，CPI）预测是发现苗头化合物、药物重定位等研究的关键技术手段。近年来，深度学习被广泛应用于CPI研究，加速了药物发现中CPI预测的发展。本文重点讨论基于特征的 CPI 预测模型，首先，介绍了CPI预测中常见的数据库、化合物和蛋白质的典型特征表示方法。根据建模中的关键问题，从多模态和注意力机制两个方面，对基于特征的CPI预测模型展开论述。在此基础上，选取其中12个模型，在3个经典数据集上评估了模型在分类任务和回归任务中的性能。本文总结当前该领域面临的挑战，对未来的发展方向进行展望，为CPI预测方法进一步研究提供思路。
-
关键词:
- 深度学习 /
- 化合物-蛋白质相互作用 /
- 药物重定位 /
- 苗头化合物筛选
Abstract: The prediction of compound-protein interaction (CPI) is a critical technological tool for discovering lead compounds and drug repurposing during the process of drug development.In recent years, deep learning has been widely used in CPI research, which has accelerated the development of CPI prediction in drug discovery.This review focuses on feature-based CPI prediction models.First, we described the datasets, as well as typical feature representation methods commonly used for compounds and proteins in CPI prediction.Based on the critical problems in modeling, we discussed models for CPI prediction from two perspectives: multimodal features and attention mechanisms.Then, the performance of 12 selected models was evaluated on 3 benchmark datasets for both classification and regression tasks.Finally, the review summarizes the existing challenges in this field and prospects for future directions.We believe that this investigation will provide some reference and insight for further research on CPI prediction. -
肝癌复发率高、转移性强，预后差， 是全球最常见的恶性肿瘤之一[1−2]。因此，对其发生发展机制的研究及后续药物的开发迫在眉睫。肿瘤免疫微环境（tumour immune microenvironment, TIME）是由巨噬细胞、淋巴细胞、树突状细胞等免疫细胞共同构成[3]。肿瘤相关巨噬细胞（tumor-associated macrophages, TAMs）是TIME中占比最多的免疫细胞[4]。TAMs浸润的数量与肝癌的恶性程度、不良预后等密切相关[5]。TAMs可通过代谢重编程等方式影响肿瘤恶性进展[6]。因此，对TAMs进行靶向调控或改造，可为肝癌的治疗提供新的策略。
脂质代谢作为肝癌细胞代谢重编程的主要标志之一，对脂质代谢进行调控可成为一种临床肝癌治疗的候选方案[7]。例如，脂肪酸合酶（fatty acid synthase, FASN）的小分子抑制剂（TVB-2640）目前已进入肝细胞癌Ⅱ期临床试验[8]。此外，TAMs在肿瘤细胞脂质合成中也发挥着重要作用，可能通过分泌促进脂质合成途径的细胞因子和生长因子，如白细胞介素-6（interleukin-6, IL-6）、肿瘤坏死因子-α（tumor necrosis factor α, TNF-α）、TGF-β，激活肿瘤细胞内的脂质合成途径或抑制脂质代谢途径，影响乙酰辅酶 A 羧化酶（acetyl-CoA carboxylase alpha, ACC1）等关键酶的表达和活性，从而促进脂质的合成和积累[9−11]，进而为肿瘤细胞的恶性增殖提供物质基础[12]。然而，关于TAMs是否能够介导肝癌细胞的脂质代谢重编程鲜有报道。因此，深入研究TAMs如何介导肝癌细胞脂质积累的机制将有助于揭示肝癌发展的新机制，为肝癌肿瘤治疗提供更多的靶点和策略。
为此，本研究首先建立了TAMs上清诱导Hepa1-6肝癌细胞的脂滴形成的共培养体系，在此基础上，探究了TAMs与Hepa1-6互作的关键因子。进一步探究了IL-10与ACC1在肝癌细胞恶性生物学中的作用。通过对肝癌细胞脂滴形成的作用机制进行初步研究，以期为肝癌的免疫治疗提供理论依据。
1. 材 料
1.1 试 剂
高糖DMEM培养基（上海源培生物科技有限公司）；牛血清（以色列Biological Industries公司）；cell counting kit-8（CCK-8）试剂盒、乳酸脱氢酶细胞毒性检测试剂盒、油红O（上海碧云天生物科技有限公司）；细胞凋亡检测试剂盒（南京诺唯赞生物科技有限公司）；脂滴绿色荧光检测试剂盒BODIPY 493/503（美国Thermo Fisher公司）；乙酰辅酶A羧化酶-α抑制剂TOFA（美国APExBIO公司）；脂滴荧光染料ND640（美国MedChemExpress公司）；乙酰辅酶A羧化酶Rabbit pAb（成都正能生物技术有限公司）；GAPDH-HRP conjugated（美国Bioworld Technology公司）；HRP标记的山羊抗兔IgG（H+L）（上海碧云天生物科技有限公司）；羧基荧光素琥珀酰亚胺酯荧光染料（CFSE，上海翌圣生物科技股份有限公司）；碘化丙啶荧光染料（PI，上海翌圣生物科技股份有限公司） ；甘油三酯（TG）测定试剂盒（南京建成生物工程研究所）； 乳酸脱氢酶（lactate dehydrogenase，LDH）细胞毒性试剂盒（上海碧云天生物技术股份有限公司）；细胞凋亡检测试剂盒（南京诺唯赞生物科技股份有限公司）。
1.2 细胞株与质粒
小鼠肝癌细胞株Hepa1-6、小鼠永生化骨髓源性巨噬细胞iBMDM（上海中科院细胞所）；ACC1-siRNA（苏州吉玛基因股份有限公司）；ACC1 CRISPR-cas9质粒实验室自己合成，并测序验证（上海生工生物工程有限公司）。
1.3 仪 器
超净工作台、CO2细胞培养箱（美国Thermo Fisher公司）；高速冷冻离心机（德国Eppendorf公司）；酶标仪（杭州奥盛仪器有限公司）；转膜仪、实时荧光定量PCR仪（美国Bio-Rad公司）；Novo Cyte流式细胞仪（美国艾森生物公司）；数码倒置相差荧光显微镜（日本Nikon公司）。
2. 方 法
2.1 细胞培养
用含10% FBS的DMEM培养基培养Hepa1-6与iBMDM细胞。将细胞放入37 ℃、5% CO2的细胞培养箱中培养。
2.2 细胞模型
以1∶1的比例分别在培养板和Transwell小室中接种iBMDM细胞和Hepa1-6细胞，培养3 d后得到TAMs。用TAMs上清液刺激Hepa1-6细胞，研究TAMs浸润对该细胞的影响。实验中，TOFA和ND-630终浓度分别为10 μmol/L和10 nmol/L。
2.3 脂滴检测
BODIPY：用2 mmol/L BODIPY493/503染料（1∶
1000 ）对Hepa1-6细胞染色，室温避光孵育20 min，PBS洗涤，进行流式分析。油红染色：参照油红O染色试剂盒说明书。
2.4 免疫印迹试验
细胞中加入蛋白裂解液，用BCA法进行定量，蛋白变性后，用8% SDS-PAGE凝胶进行电泳，转移至PVDF膜上，用5% BSA室温封闭1 h，Acetyl Coenzyme A Carboxylase Rabbit pAb、GAPDH-HRP conjugated过夜封闭，TBST洗脱6次共30 min。HRP标记的山羊抗兔IgG（H+L）室温封闭1 h，TBST洗脱，ECL化学发光法显影，化学发光仪拍照。
2.5 甘油三酯含量检测
使用试剂盒检测甘油三酯含量，详细实验步骤参考该说明书。
2.6 siRNA干扰
将Hepa1-6按每孔1×106个的密度接种于12孔细胞培养板中，待细胞贴壁后进行转染。
2.7 ACC1 CRISPR-cas9质粒构建
使用在线工具（https://benchling.com）设计靶向鼠源ACC1不同外显子区域的gRNA，序列如下：
Table 1. Sequence for guide RNA (gRNA) about acetyl-CoA carboxylase alpha (ACC1)gRNA Forward primer Reverse primer gRNA-1 CACCGAATGCATGCGATCTATCCGT AAACACGGATAGATCGCATGCATTC gRNA-2 CACCGAAGTGTATCTGAGCTGACGG AAACCCGTCAGCTCAGATACACTTC gRNA-3 CACCGCAAACGTGAATGCTTGACCA AAACTGGTCAAGCATTCACGTTTGC 将Lenti-CRPISPRv2载体用BsmBⅠ酶酶切后回收。连接，转化，挑选单克隆进行扩大培养，提质粒后测序验证。后续进行慢病毒包装，侵染Hepa1-6细胞。
2.8 细胞活力检测（CCK-8法）
在96孔细胞培养板中，每孔加入200个Hepa1-6细胞，总体积为100 µL；然后按1∶10的比例加入CCK-8溶液；37 ℃孵育2~4 h；用酶标仪检测样品的吸收度（450 nm）。
2.9 细胞增殖（CFSE法）
首先用PBS清洗细胞，在
1000 r/min条件下离心4~5 min后收集细胞；向细胞中按1∶500的比例加入CFSE，然后置培养箱中孵育25~30 min；用PBS洗去未结合的CFSE并取一部分细胞置4 ℃冰箱中作为阳性对照；将实验组细胞继续避光培养48 h；最后收集细胞，用流式细胞仪对其进行检测。2.10 细胞周期检测
用75%的乙醇重悬细胞，−20 ℃过夜，离心弃去废液，PBS洗涤后重悬细胞，加入RNA裂解酶，37 ℃水浴15 min；然后加入PI染液，避光静置30 min；PBS洗涤后用流式细胞仪分析。
2.11 细胞毒性检测
使用乳酸脱氢酶（lactate dehydrogenase，LDH）细胞毒性试剂盒检测 ，实验步骤参考说明书。
2.12 细胞凋亡
使用细胞凋亡检测试剂盒进行检测，实验步骤参考说明书（A211-01）。
2.13 统计学分析
用GraphPadPrism8软件统计分析。数据以算数$\bar{x}\pm s $形式表示，组间两两比较采用Student’s t-test检验进行显著性检验。3组以上比较采用One-Way ANOVA方差分析进行显著性检验。P <0.05认为是存在显著性差异。
3. 结 果
3.1 TAMs浸润诱导Hepa1-6细胞积累LDs
将鼠肝癌细胞Hepa1-6和iBMDM进行共培养（TAMs组），使用q-PCR与流式检测未共培养或共培养后的iBMDM的表型，实验发现M2相关标志物精氨酸酶-1（arginase-1, ARG1）、血管内皮生长因子A（vascular endothelial growth factor A, VEGFA）、IL-10的mRNA水平显著升高，细胞程序性死亡配体1（programmed cell death ligand 1，PD-L1）的表达也有升高趋势（P < 0. 01，图1-A、图1-B），结果表明共培养后iBMDM偏向M2型且PD-L1表达升高。
![]()
Figure 1. tumor-associated macrophage (TAMs) model establishmentA: q-PCR for CD206, ARG1, VEGFα, IL-10, inducible nitric oxide synthase (iNOS), interleukin-12 subunit alpha (IL12A), TNFα mRNA expression levels; B: flow cytometry was used to detect the expression of PD-L1 in TAMs*P < 0. 05,**P < 0. 01,***P < 0. 001使用浓度分别为10%、20%、30%、40%的TAMs上清液（即TAMs上清液占完全培养基的比例）刺激Hepa1-6细胞。流式结果显示，随着TAMs上清液浓度的增加，Hepa1-6细胞LDs积累量逐渐增高，且TAMs上清液含量为30%和40%时并未出现显著性差异，即使用30%的TAMs上清液即可增强Hepa1-6细胞LDs的积累（P < 0. 05，图2-A、图2-B）。当用30% TAMs上清液刺激Hepa1-6细胞1，2，3 d后，LDs积累量呈现逐渐上升趋势（P <0.05，图2-C、图2-D）。因此，构建了一个TAMs浸润诱导肝癌LDs积累细胞模型，该模型用30% TAMs上清液与Hepa1-6细胞共同培养3 d，在此过程中，Hepa1-6细胞LDs积累量增加，这些LDs积累量增加的Hepa1-6细胞被标记为LDhigh Hepa1-6。油红染色实验结果表明，TAMs上清液刺激后，Hepa1-6细胞LDs增加（图2-E）。采用q-PCR检测细胞内围脂滴蛋白家族的mRNA水平表达量，结果表明脂蛋白1（Plin1）、脂蛋白4（Plin4）、脂蛋白5（Plin5）基因表达量显著性提高（P <0.05，图2-F）。以上结果表明，TAMs浸润可以诱导Hepa1-6细胞LDs的积累。
![]()
Figure 2. LDhigh Hepa1-6 cell model constructs and TAMs infiltration induced accumulation of LDs in hepatocellular carcinoma cellsA and B: Hepa1-6 was stimulated with different concentrations of TAMs supernatants (10%, 20%, 30%, 40%), stained with BODIPY 493/503, and LDs accumulation in Hepa1-6 was detected by flow cytometry; C and D: 30% TAMs supernatants were stimulated for 1, 2 and 3 days, stained with BODIPY 493/503 and flow cytometry was used to detect the accumulation of LDs in Hepa1-6; E: Oil Red O labelled LDs; F: q-PCR for Plin 1, Plin 2, Plin 3, Plin 4, Plin5 mRNA expression levels*P < 0. 05, **P < 0. 01,***P < 0. 001 vs control group3.2 TAMs通过分泌IL-10影响肝癌细胞积累LDs
为了深入了解TAMs如何诱导Hepa1-6细胞累积LDs，使用q-PCR对IL-10、TGF-β、白细胞介素12（IL-12）和TNF-α的mRNA表达水平进行了检测。结果显示，TAMs组IL-10和TGF-β的表达水平均较高，其中IL-10的表达水平尤其显著（P <0.001，图3-A）。为进一步探讨IL-10对Hepa1-6细胞LDs积累的影响，采用不同质量浓度（1、5和10 ng/mL）的IL-10刺激Hepa1-6细胞24 h，用BODIPY493/503标记的脂滴，流式结果显示，当IL-10质量浓度为10 ng/mL时，Hepa1-6细胞的LDs含量显著增加（P <0.001，图3-B、图3-C）。以上结果表明，肿瘤相关巨噬细胞分泌的IL-10可能是促进Hepa1-6细胞积累脂滴的重要机制之一。
![]()
Figure 3. IL-10 secretion by TAMs induces accumulation of LDs in hepatocellular carcinoma cells A: q-PCR for IL-10, TGFβ, IL-12, TNFα mRNA expression levels; B and C: Stimulation with different concentrations of IL-10, BODIPY 493/503 staining, and flow cytometric detection of LDs accumulation in Hepa1-6, MFI: Mean fluorescence intensity*P < 0. 05,**P < 0. 01,***P < 0. 001为了进一步说明IL-10在TAMs中的作用，采用IL-10封闭抗体处理LDhigh Hepa1-6，与未处理的LDhigh Hepa1-6相比，其脂滴生成明显减少（图4-A），q-PCR结果显示，IL-10封闭后的LDhigh Hepa1-6 ACC1、Plin4、Plin5表达量有降低趋势（图4-B，P <0.05）。流式结果显示，IL-10封闭后抑制了LDhigh Hepa1-6 细胞增殖能力，促进了LDhigh Hepa1-6细胞凋亡（图4-C~图4-E）。以上结果表明，TAMs可以通过分泌IL-10促进肝癌细胞LDs的积累。
![]()
Figure 4. IL-10 deletion affects lipid droplet formation, cell proliferation and apoptosisA: Oil Red O labelled LDs; B: q-PCR for ACC1, Plin4, Plin5 mRNA expression levels; C: CFSE staining labelling to detect the effect of IL-10 deletion on the proliferation of Hepa1-6 cells; D and E: Effects of IL-10 deletion on apoptosis in Hepa1-6 cells*P < 0. 05 ,**P < 0. 01, ***P < 0. 001 vs control group3.3 肝癌细胞中ACC1高表达且与肝癌不良预后密切相关
为了进一步探究肝癌细胞积累LDs的机制, q-PCR检测了脂质合成相关基因ACC1，FASN，酰基辅酶A合成酶长链家族成员（acyl coenzyme A synthetase long chain family, ACSL），二酯酰甘油酰基转移酶1（diacylglycerol-O-acyltransferase homolog 1, DGAT1）与酯酰基转移酶2（diacylglycerol-O-acyltransferase homolog 2, DGAT2）的表达水平。结果显示，TAMs上清液刺激后, LDhigh Hepa1-6细胞中ACC1的表达显著增高（P <0.01，图5-A）。同时，Western blot结果显示LDhigh Hepa1-6细胞中ACC1的表达水平也显著提升（图5-B）。
![]()
Figure 5. High expression of ACC1 in LDhigh Hepa1-6 cells and implications for cancerA: q-PCR for ACC1, FASN, ACSL, DGAT1, DGAT2 mRNA expression levels; B: Western blotting for ACC1 protein levels; C: ACC1 expression levels in clinical samples of hepatocellular carcinoma compared with normal tissues; D: ACC1 expression levels in different hepatocellular carcinoma stages (stage1, stage2, stage3, stage4); E: ACC1 expression levels in different hepatocellular carcinoma grades (grade1, grade2, grade3, grade4); F: Survival curves of patients with ACC1 high expression and F: Survival curves of patients with high ACC1 expression and low ACC1 expression; G: Survival curves of patients with high ACC1 expression and low ACC1 expression in different liver cancer grades (grade1, grade2, grade3, grade4) (*P < 0. 05,**P < 0. 01,***P < 0. 001)基于以上结果，本研究利用TCGA数据库对LIHC肝癌临床样本中ACC1表达水平进行分析，结果显示，与正常组织相比，肝癌组织中ACC1的表达水平显著上调，且ACC1的表达随着HCC的恶化有增高趋势（P <0.01，图5-C~图5-E）。因此推测ACC1可能调控肝癌细胞LDs的积累。
接下来，利用TCGA数据库，对ACC1高表达与HCC患者生存期之间的关系进行研究，结果显示，ACC1表达水平与HCC患者生存呈负相关性，表明ACC1高表达不利于HCC患者生存（P <0.001，图5-F-G）。
3.4 ACC1介导LDhigh Hepa1-6细胞LDs积累
为了探究ACC1是否调控Hepa1-6细胞LDs积累，Western blot结果显示，ACC1抑制剂TOFA和ND-630能抑制ACC1表达（图6-A），与其他人报道一致。分别用BODIPY493/503和油红O标记LDs，结果显示，TOFA和ND-630能显著减少LDhigh Hepa1-6细胞中LDs的积累（P < 0. 001，图6-B~图6-D）。对Hepa1-6细胞中甘油三酯（TG）含量进行检测，结果显示，TOFA和ND-630处理后，LDhigh Hepa1-6细胞中的TG含量显著减少（P <0.05，图6-E）。
![]()
Figure 6. ACC1 mediates the accumulation of LDs in LDhigh Hepa1-6 cellsA: Western blotting to detect ACC1 protein level; B and C: BODIPY 493/503 staining and flow cytometry to detect the accumulation of LDs in TOFA and ND-630-treated Hepa1-6 cells, MFI: Mean fluorescence intensity; D: Oil Red O labelling of LDs; E: TG content of the cells; F: q-PCR to detect the expression level of ACC1 mRNA and assess the effect of ACC1-siRNA interference; G: Western blotting to detect ACC1 expression level and assess the effect of ACC1-siRNA interference; H and I: BODIPY 493/503 staining and flow cytometry to detect the effect of siRNA interference on the accumulation of LDs in Hepa1-6, MFI: Mean fluorescence intensity (*P < 0. 05,**P < 0. 01,***P < 0. 001 vs control group)利用siRNA技术设计了3条siRNA（ACC1-siRNA-563，ACC1-siRNA-1722，ACC1-siRNA-4838），进行ACC1敲减细胞模型构建，并通过q-PCR和Western blot检测了干扰效果。结果表明，3条siRNA均能降低Hepa1-6细胞中ACC1的表达，其中ACC1-siRNA-4838效果最好，ACC1敲减效率为68%（P <0.001，图6-F~图6-G）。接下来，对LDs的积累进行检测，结果显示，ACC1-siRNA-4838干扰后，LDs积累减少（P <0.05，图6-H~图6-I）。
此外，构建了慢病毒包装的ACC1 CRISPR-cas9质粒，侵染Hepa1-6细胞敲除ACC1，Western blot验证敲除效果，结果显示敲除后ACC1明显减少（标记为ACC1-KO Hepa1-6）（图7-A）。用BODIPY 493/503标记LDs，结果显示TAMs上清刺激后，ACC1-KO Hepa1-6细胞LDs显著减少（P <0.01，图7-B、图7-C）。通过检测TG含量发现，相比于LDhigh Hepa1-6细胞，TAMs上清液刺激后，ACC1-KO Hepa1-6细胞中TG含量显著减少（P <0.01，图7-D）。
![]()
Figure 7. ACC1 knockdown inhibits the accumulation of LDs in LDhigh Hepa1-6 cellsA: Western blotting to detect the expression level of ACC1 and evaluate the effect of ACC1 knockdown; B and C: BODIPY 493/503 staining and flow cytometry to detect the effect of ACC1 knockdown on the accumulation of LDs in Hepa1-6; D: triglyceride (TG) content of the cells, MFI: Mean fluorescence intensity*P < 0. 05,**P < 0. 01,***P < 0. 001通过以上3种方式干预ACC1（抑制剂、小干扰RNA和敲除）来研究其对Hepa1-6细胞LDs积累的影响。结果表明，抑制或敲除ACC1可以显著减少脂肪积累，表明ACC1可能是调控细胞内脂肪积累的关键分子。
3.5 ACC1调控LDhigh Hepa1-6细胞增殖
为了探究ACC1介导的LDs积累是否会促进Hepa1-6细胞恶性生物学行为，首先研究了ACC1对Hepa1-6细胞毒性以及增殖的影响。TOFA和ND-630处理24 h后，Hepa1-6细胞LDH释放显著增加，表明抑制ACC1会促进Hepa1-6细胞死亡（P <0.001，图8-A）。CCK-8结果显示，使用TOFA和ND-630处理后，LDhigh Hepa1-6细胞活力显著降低，说明抑制ACC1会降低Hepa1-6细胞增殖（P <0.001，图8-B）。相比于LDhigh Hepa1-6细胞，TAMs上清刺激后，ACC1-KO Hepa1-6细胞活力降低，表明敲除ACC1会使Hepa1-6细胞增殖能力降低（P <0.001，图8-C）。
![]()
Figure 8. ACC1 affects Hepa1-6 cell viabilityA: Detection of LDH release to assess the toxicity of TOFA and ND-630 on Hepa1-6 cells; B: CCK-8 to detect the effect of TOFA and ND-630 on the viability of Hepa1-6 cells; C: CCK-8 to detect the effect of ACC1 knockdown on the viability of Hepa1-6 cells***P < 0. 001, **P < 0. 01, *P < 0. 05为了进一步评估肿瘤细胞增殖能力，用CFSE对Hepa1-6细胞进行荧光标记。流式结果表明TOFA和ND-630会抑制LDhigh Hepa1-6细胞的增殖（图9-A）。同时，敲除ACC1的Hepa1-6细胞增殖能力显著降低（图9-B）。细胞周期检测结果显示，TOFA和ND-630处理后，S期细胞比例增加，即抑制ACC1会阻滞细胞周期在S期（图9-C~9-D）。此外，对ACC1-KO Hepa1-6细胞周期分析显示，其S期细胞比例增加，进一步验证了抑制ACC1会阻滞细胞周期在S期（图9-E~9-F）。因此，以上结果表明ACC1调控Hepa1-6细胞的增殖能力，抑制ACC1可以将细胞周期阻滞在S期，导致肿瘤细胞增殖能力降低。
![]()
Figure 9. ACC1 regulates Hepa1-6 cell proliferationA: CFSE staining labeling to detect the effect of TOFA and ND-630 on the proliferation of Hepa1-6 cells; B: CFSE staining labeling to detect the effect of ACC1 knockdown on the proliferation of Hepa1-6 cells; C and D: Effect of TOFA and ND-630 on the cell cycle of Hepa1-6 cells; E and F: The effect of ACC1 knockdown on the cell cycle of Hepa1-6 cells3.6 ACC1影响LDhigh Hepa1-6细胞凋亡
用AnnexinV-FITC和PI对细胞进行标记，流式分析结果显示，TOFA处理后，LDhigh Hepa1-6早期凋亡无显著差异，晚期凋亡显著增强，用ND-630处理后，细胞早期凋亡和晚期凋亡都显著增强（图10-A、图10-B）。并且与LDhigh Hepa1-6细胞相比，ACC1-KO Hepa1-6细胞早期凋亡和晚期凋亡均显著增强（图10-C、图10-D）。实验结果表明，ACC1调控Hepa1-6细胞凋亡，抑制ACC1会促进肿瘤细胞凋亡。
![]()
Figure 10. ACC1 regulates apoptosis in Hepa1-6 cellsA and B: Effects of TOFA and ND-630 on apoptosis in Hepa1-6 cells; C and D: Effects of ACC1 knockdown on apoptosis in Hepa1-6 cells*P < 0. 05,**P < 0. 01,***P < 0. 0014. 讨 论
脂质代谢作为肝癌代谢重编程的主要标志之一，在肝癌细胞中产生更多的脂滴，以此为肝癌细胞的恶性生长与转移提供能量[13−14]。因此，对肝癌中脂滴形成的机制进行探究，有利于为肝癌的潜在靶点提供理论基础。本研究建立了iBMDM上清液诱导Hepa1-6产生脂滴的共培养模型，即iBMDM与Hepa1-6共培养3 d后，得到TAMs。TAMs上清与Hepa1-6共培养3 d后便可有效诱导肿瘤细胞的脂滴生成。并在此基础上，本研究探究了脂滴形成促进肝癌细胞恶性生物学行为的机制，旨在对肝癌的免疫治疗提供新的思路。
TAMs是大多数肿瘤中浸润最多的一类免疫细胞，通过M2表型、促进血管生成、产生基质金属蛋白酶和产生抑制性受体等机制在肿瘤发生、发展和转移中发挥关键性作用。基于清除TAMs、抑制TAMs募集和对TAMs进行重编程等手段，靶向TAMs已成为肿瘤免疫治疗的新型策略[15]。Zhang等[16]发现TAMs分泌白细胞介素6（IL-6）能够作用于胶质瘤细胞，促进胶质瘤细胞中3-磷酸肌醇蛋白激酶1（PDPK1）介导的磷酸化，进而促进了胶质瘤细胞的有氧糖酵解和肿瘤生长。使用IL-6中和抗体可以抑制经TAMs促进的胶质瘤细胞的有氧糖酵解和肿瘤生长。另外研究发现，TAMs通过摄取较高的葡萄糖，以此来增强糖胺生物合成，并伴随分泌更多的组织蛋白酶B到TME中，进而促进肿瘤转移和化疗耐药性[17]。然而，关于TAMs如何调控肿瘤细胞脂质合成的报道较少。本研究发现，TAMs通过释放IL-10可促进Hepa1-6细胞的脂滴囤积，进而促进肿瘤的生长。IL-10可通过IL-10-DDIT4-mTOR通路影响细胞的糖脂代谢，但巨噬细胞分泌的IL-10如何影响肝癌细胞的脂滴形成，仍需做进一步研究。
肝癌细胞中，脂肪酸主要通过FA合成产生。ACC1是脂肪酸从头合成第一步的限速酶，是脂质生物合成的关键代谢酶[18−19]。本研究发现LDhigh Hepa1-6细胞中的ACC1表达水平显著升高。同时，TCGA数据表明，肝癌患者肿瘤组织中ACC1的表达水平也显著提高，且ACC1的高表达与HCC患者的不良预后有密切关系。在本研究中，抑制ACC1表达后，LDhigh Hepa1-6细胞中LDs积累量显著减少，表明ACC1是调控LDs积累的关键分子。
LDs在维持脂质稳态、调节细胞应激、蛋白质处理等方面发挥着至关重要的作用。本研究采用了3种模型：ACC1抑制剂TOFA、ND-630，以及ACC1 siRNA干扰和ACC1 CRISPR-cas9敲除，对ACC1进行抑制，旨在探讨其在肝癌细胞恶性行为中的影响。实验结果显示，抑制ACC1可有效减少LDs的积累，抑制LDhigh Hepa1-6细胞的增殖，并促进细胞凋亡的发生。这些结果表明，ACC1是调控LDs积累的关键因子。研究表明，ACC1在肝癌细胞中的表达水平通常较高，与肝癌的发生和发展密切相关[20]。ACC1在肝癌细胞中扮演着重要角色，一是通过合成脂肪酸满足能量需求，二是调节脂滴形成，维持肝癌细胞生存和增殖[21]。此外，在肝癌前期，由于细胞受到损伤或炎症刺激，ACC1的表达水平可能会逐渐上调[22]。在肝癌发展过程中，ACC1的表达和活性可能会进一步增加，以满足肿瘤细胞快速增殖和生长的需求[23]。然而，在肝癌晚期，肿瘤细胞可能会经历代谢异常，导致脂质代谢紊乱和脂滴形成受到抑制[24]。在这种情况下，ACC1的表达和活性可能会下降，因为肿瘤细胞对脂肪酸的需求减少。这种变化可能与肝癌细胞的代谢适应和生长状态密切相关。因此，ACC1在不同肝癌病理阶段的表达和活性变化可能反映了肝癌细胞代谢适应性的调节和肿瘤生长状态的变化。总之，ACC1在肝癌细胞中的作用机制及其潜在的临床意义值得深入研究，有望为肝癌的治疗提供新的思路和方法。
综上所述，本研究揭示了TAMs能促进肝癌细胞的脂滴形成，且TAMs释放的IL-10可能会通过调控ACC1的表达在此过程中发挥作用。肝癌细胞通过ACC1介导的脂质从头合成在其增殖和凋亡等恶性生物学中扮演了关键作用。本研究确认，ACC1是肝癌细胞脂质代谢和恶性行为调控的关键因子，其抑制不仅可以减少TAMs促进的肝癌细胞脂滴积累，还可以抑制肝癌细胞的增殖并诱导凋亡，显示了其作为治疗靶点的潜力。
本刊编委徐寒梅教授团队在Signal Transduction and Targeted Therapy上对核酸药物的最新研究及
未来发展进行总结与展望近日，中国药科大学江苏省合成多肽药物发现与评价工程研究中心徐寒梅教授团队在Signal Transduction and Targeted Therapy（IF=40.8）杂志发表题为“Nucleic acid drugs: recent progress and future perspectives”的综述文章。本文以孙小艺、Sarra Setrerrahmane、李臣诚、胡加亮为主要作者，徐寒梅教授为通信作者。中国药科大学为该论文第一通信单位。 全文链接：https://doi.org/10.1038/s41392-024-02035-4 -
[1] Fu MY, Zhu YY, Wu CY, et al. Prediction of plasma protein binding rate based on machine learning[J]. J China Pharm Univ (中国药科大学学报), 2021, 52(6): 699-706. [2] Zhong FS. Studies on drug-target interactions based on graph neural network(基于图神经网络的药物—靶标作用研究)[D]. Shanghai: Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 2022. [3] Kim S, Chen J, Cheng TJ, et al. PubChem in 2021: new data content and improved web interfaces[J]. Nucleic Acids Res, 2021, 49(D1): D1388-D1395. [4] Bento AP, Gaulton A, Hersey A, et al. The ChEMBL bioactivity database: an update[J]. Nucleic Acids Res, 2014, 42(Database issue): D1083-D1090. [5] Szklarczyk D, Santos A, von Mering C, et al. STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data[J]. Nucleic Acids Res, 2016, 44(D1): D380-D384. [6] Liu TQ, Lin Y, Wen X, et al. BindingDB: a web-accessible database of experimentally determined protein-ligand binding affinities[J]. Nucleic Acids Res, 2007, 35(Database issue): D198-D201. [7] Su MY, Yang QF, Du Y, et al. Comparative assessment of scoring functions: the CASF-2016 update[J]. J Chem Inf Model, 2019, 59(2): 895-913. [8] Mysinger MM, Carchia M, Irwin JJ, et al. Directory of useful decoys, enhanced (DUD-E): better ligands and decoys for better benchmarking[J]. J Med Chem, 2012, 55(14): 6582-6594. [9] Wishart DS, Feunang YD, Guo AC, et al. DrugBank 5.0: a major update to the DrugBank database for 2018[J]. Nucleic Acids Res, 2018, 46(D1): D1074-D1082. [10] Davis MI, Hunt JP, Herrgard S, et al. Comprehensive analysis of kinase inhibitor selectivity[J]. Nat Biotechnol, 2011, 29(11): 1046-1051. [11] Tang J, Szwajda A, Shakyawar S, et al. Making sense of large-scale kinase inhibitor bioactivity data sets: a comparative and integrative analysis[J]. J Chem Inf Model, 2014, 54(3): 735-743. [12] Du PF, Wang X, Xu C, et al. PseAAC-Builder: a cross-platform stand-alone program for generating various special Chou''s pseudo-amino acid compositions[J]. Anal Biochem, 2012, 425(2): 117-119. [13] Chen Z, Zhou Y, Song JN, et al. hCKSAAP_UbSite: improved prediction of human ubiquitination sites by exploiting amino acid pattern and properties[J]. Biochim Biophys Acta, 2013, 1834(8): 1461-1467. [14] Wan XZ, Wu XL, Wang DY, et al. An inductive graph neural network model for compound-protein interaction prediction based on a homogeneous graph[J]. Brief Bioinform, 2022, 23(3): bbac073 .[15] Chen Y, Wu H, Ge WH, et al. Research on entity relation extraction of Chinese adverse drug reaction reports based on deep learning method[J]. J China Pharm Univ (中国药科大学学报), 2019, 50(6): 753-759. [16] ?ztürk H, ?zgür A, OzkirimLi E. DeepDTA: deep drug-target binding affinity prediction[J]. Bioinformatics, 2018, 34(17): i821-i829. [17] ?ztürk H, OzkirimLi E, ?zgür A. WideDTA: prediction of drug-target binding affinity[J].arXiv,2019:1902. 04166 v1.[18] Yang JB, Cai YY, Zhao KR, et al. Concepts and applications of chemical fingerprint for hit and lead screening[J]. Drug Discov Today, 2022, 27(11): 103356. [19] Lee I, Keum J, Nam H. DeepConv-DTI: prediction of drug-target interactions via deep learning with convolution on protein sequences[J]. PLoS Comput Biol, 2019, 15(6): e1007129 .[20] Rifaioglu AS, Cetin Atalay R, Cansen Kahraman D, et al. MDeePred: novel multi-channel protein featurization for deep learning-based binding affinity prediction in drug discovery[J]. Bioinformatics, 2021, 37(5): 693-704. [21] Wang SD, Du ZZ, Ding M, et al. LDCNN-DTI: a novel light deep convolutional neural network for drug-target interaction predictions[C]// 2020 IEEE International Conference on Bioinformatics and Biomedicine (BIBM) . Seoul:IEEE, 2021: 1132-1136.[22] Wang KL, Zhou RY, Li YH, et al. DeepDTAF: a deep learning method to predict protein-ligand binding affinity[J]. Brief Bioinform, 2021, 22(5): bbab072 .[23] Wang S, Jiang MJ, Zhang SG, et al. MCN-CPI: multiscale convolutional network for compound-protein interaction prediction[J]. Biomolecules, 2021, 11(8): 1119. [24] Lin X. DeepGS: Deep Representation Learning of graphs and sequences for drug-target binding affinity prediction[J]. arXiv,2020: 2003.13902v2. [25] Wei LS, Long WT, Wei LY. MDL-CPI: multi-view deep learning model for compound-protein interaction prediction[J]. Methods, 2022, 204: 418-427. [26] Nguyen T, Le H, Quinn TP, et al. GraphDTA: predicting drug-target binding affinity with graph neural networks[J]. Bioinformatics, 2021, 37(8): 1140-1147. [27] Chen S, Sun Z, Lin LH, et al. To improve protein sequence profile prediction through image captioning on pairwise residue distance map[J]. J Chem Inf Model, 2020, 60(1): 391-399. [28] Wang S, Sun SQ, Li Z, et al. Accurate de novo prediction of protein contact map by ultra-deep learning model[J]. PLoS Comput Biol, 2017, 13(1): e1005324 .[29] Jiang MJ, Li Z, Zhang SG, et al. Drug-target affinity prediction using graph neural network and contact maps[J]. RSC Adv, 2020, 10(35): 20701-20712. [30] Michel M, Menéndez Hurtado D, Elofsson A. PconsC4: fast, accurate and hassle-free contact predictions[J]. Bioinformatics, 2019, 35(15): 2677-2679. [31] Guo BJ, Zheng HY, Jiang HH, et al. Enhanced compound-protein binding affinity prediction by representing protein multimodal information via a coevolutionary strategy[J]. Brief Bioinform, 2023, 24(2): bbac628 .[32] Li SY, Wan FP, Shu HT, et al. MONN: a multi-objective neural network for predicting compound-protein interactions and affinities[J]. Cell Syst, 2020, 10(4): 308-322.e11. [33] Yang ZD, Zhong WH, Zhao L, et al. MGraphDTA: deep multiscale graph neural network for explainable drug-target binding affinity prediction[J]. Chem Sci, 2022, 13(3): 816-833. [34] Gao KY, Fokoue A, Luo H, et al. Interpretable drug target prediction using deep neural representation[C]// Proceedings of the Twenty-Seventh International Joint Conference on Artificial Intelligence . California: International Joint Conferences on Artificial Intelligence Organization, 2018: 3371-3377.[35] Zhao Q, Xiao F, Yang M, et al. AttentionDTA: prediction of drug-target binding affinity using attention model[C]// IEEE International Conference on Bioinformatics and Biomedicine (BIBM) . San Diego: IEEE, 2019: 64-69.[36] Chen LF, Tan XQ, Wang DY, et al. TransformerCPI: improving compound-protein interaction prediction by sequence-based deep learning with self-attention mechanism and label reversal experiments[J]. Bioinformatics, 2020, 36(16): 4406-4414. [37] Nguyen NQ, Jang G, Kim H, et al. Perceiver CPI: a nested cross-attention network for compound-protein interaction prediction[J]. Bioinformatics, 2023, 39(1): btac731 .[38] Huang KX, Xiao C, Glass LM, et al. MolTrans: molecular Interaction Transformer for drug-target interaction prediction[J]. Bioinformatics, 2021, 37(6): 830-836. [39] Wang J, Dokholyan NV. Yuel: improving the generalizability of structure-free compound-protein interaction prediction[J]. J Chem Inf Model, 2022, 62(3): 463-471. [40] Bai PZ, Miljkovi? F, John B, et al. Interpretable bilinear attention network with domain adaptation improves drug-target prediction[J]. Nat Mach Intell, 2023, 5(2): 126-136. [41] Zhao QC, Yang MY, Cheng ZJ, et al. Biomedical data and deep learning computational models for predicting compound-protein relations[J]. IEEE/ACM Trans Comput Biol Bioinform, 2022, 19(4): 2092-2110. [42] Wang SD, Du ZZ, Ding M, et al. KG-DTI: a knowledge graph based deep learning method for drug-target interaction predictions and Alzheimer''s disease drug repositions[J]. Appl Intell, 2022, 52(1): 846-857. [43] Li YF, Sun C, Wei JM, et al. Drug-Protein interaction prediction by correcting the effect of incomplete information in heterogeneous information[J]. Bioinformatics, 2022, 38(22): 5073-5080. [44] Ma WJ, Zhang SG, Li Z, et al. Predicting drug-target affinity by learning protein knowledge from biological networks[J]. IEEE J Biomed Health Inform, 2023, 27(4): 2128-2137. [45] Wu F, Jin ST, Jiang YH, et al. Pre-training of equivariant graph matching networks with conformation flexibility for drug binding[J]. Adv Sci, 2022, 9(33): e2203796 .[46] Wan FP, Zhu Y, Hu HL, et al. DeepCPI: a deep learning-based framework for large-scale in silico drug screening[J]. Genomics Proteomics Bioinformatics, 2019, 17(5): 478-495. [47] Nguyen TM, Nguyen T, Tran T. Mitigating cold-start problems in drug-target affinity prediction with interaction knowledge transferring[J]. Brief Bioinform, 2022, 23(4): bbac269 .
下载:
计量
- 文章访问数: 342
- HTML全文浏览量: 31
- PDF下载量: 357
