《Frontiers in Bioinformatics》：Targeting IDH-driven Oncometabolism in glioblastoma: a computational approach for the identification of potential inhibitors
编辑推荐：
中文摘要翻译：
背景：胶质母细胞瘤（GBM）是一种恶性脑肿瘤，常由异柠檬酸脱氢酶（IDH1和IDH2）基因突变驱动，这些突变促进肿瘤代谢物2-羟基戊二酸（2-HG）的促新合成，进而导致代谢功能障碍。尽管这些突变与不同的临床结局相关，但由于肿瘤异质性、细胞可塑性
中文摘要翻译：
背景：胶质母细胞瘤（GBM）是一种恶性脑肿瘤，常由异柠檬酸脱氢酶（IDH1和IDH2）基因突变驱动，这些突变促进肿瘤代谢物2-羟基戊二酸（2-HG）的促新合成，进而导致代谢功能障碍。尽管这些突变与不同的临床结局相关，但由于肿瘤异质性、细胞可塑性、血脑屏障通透性受限、药物外排机制和有效的DNA修复通路，治疗干预仍具有挑战性。因此，本研究旨在通过克服这些挑战来识别突变IDH蛋白的选择性抑制剂。方法：使用来自ChemFaces和MedChemExpress库的10,309种化合物进行大规模虚拟筛选和多级ADMET过滤，筛选出19种具有良好类药性的先导化合物。进一步进行了分子对接、300 ns分子动力学模拟、PCA分析和MM-PBSA计算以探索构象动力学。结果：分子对接显示，PubChem ID 91457（I1Hit1）对IDH1的结合亲和力为?7.57 kcal/mol，PubChem ID 4946（I1Hit2）为?6.82 kcal/mol，而PubChem ID 71550939（I2Hit1）对IDH2的结合亲和力为?11.23 kcal/mol，PubChem ID 4946（I2Hit2）为?7.87 kcal/mol，均优于标准对照（?6.17 kcal/mol）。分子动力学模拟、PCA分析和MM-PBSA计算确认了复合物的稳定性和抑制潜力。值得注意的是，I1Hit1和I2Hit1分别对IDH1和IDH2显示出增强的结合自由能?20.67 ± 2.32 kcal/mol和?28.78 ± 2.67 kcal/mol，优于标准对照，且计算药效团评估进一步支持了这一结果。结论：综合而言，这些发现强调了I1Hit1和I2Hit1作为新型治疗化合物，能够有效抑制IDH靶点以应对GBM中的表观遗传修饰，这些化合物需要进一步的实验验证以证明其作为IDH驱动代谢的潜在抑制剂在GBM中的作用。
论文解读文章：
**研究背景与问题**
胶质母细胞瘤（GBM）是中枢神经系统中最具侵袭性的原发性脑肿瘤，占所有胶质瘤的约45%和恶性脑肿瘤的近80%，其年发病率在美国约为3.19/10万，且男性略多（男女比约1.6:1）。GBM患者预后极差，即使接受标准治疗（手术切除、放疗和替莫唑胺化疗），中位生存期仅约14个月，5年生存率低于5%。近年来研究发现，异柠檬酸脱氢酶（IDH）基因突变在包括胶质瘤在内的多种实体瘤中频繁出现。IDH1和IDH2的体细胞功能获得性突变导致肿瘤代谢物2-羟基戊二酸（2-HG）异常积累，2-HG竞争性抑制α-酮戊二酸依赖性双加氧酶（如TET甲基胞嘧啶双加氧酶和JmjC结构域组蛋白去甲基化酶），引发DNA和组蛋白高甲基化表型，导致表观遗传重编程和代谢失调，从而驱动肿瘤发生发展。IDH1突变主要发生在R132残基（以R132H最为常见，约占89%），IDH2突变常见于R140和R172残基。目前已有IDH抑制剂如ivosidenib和vorasidenib在临床试验中显示出初步疗效，但GBM治疗仍面临多重挑战：肿瘤异质性、细胞可塑性、血脑屏障（BBB）通透性限制、药物外排泵活性、胶质瘤干细胞耐药性及DNA损伤修复机制等。此外，IDH突变常与TP53、ATRX、EGFR、MGMT等基因改变共存，这些共突变对治疗反应有独立影响。因此，迫切需要开发能够有效穿越血脑屏障、选择性抑制突变IDH酶且具有良好安全性的新型治疗策略。本研究旨在通过整合计算生物学方法，从大规模化合物库中筛选和鉴定针对突变IDH1和IDH2的潜在抑制剂。
**研究方法概述**
研究人员从ChemFaces（970种化合物）和MedChemExpress（9339种化合物）两个抗癌化合物库中获取共10,309种化合物，经过多级ADMET（吸收、分布、代谢、排泄和毒性）筛选，包括Lipinski五规则、血脑屏障通透性、PAINS、生物利用度、GI吸收、致突变性、致癌性、毒效团、肝毒性（H-HT）和药物性肝损伤（DILI）等标准，最终获得19种先导化合物。随后对排名靠前的化合物进行分子对接（使用AutoDockTools-1.5.7，采用Lamarckian遗传算法），针对IDH1（PDB:6VEI）和IDH2（PDB:6VFZ）的活性位点进行对接。接着，对最优复合物进行300 ns的分子动力学（MD）模拟（使用GROMACS 2023.1，CHARMM27力场，TIP3P水模型），并分析RMSD、RMSF、Rg、SASA和氢键。随后进行主成分分析（PCA）和自由能景观（FEL）分析，以及MM-PBSA结合自由能计算。最后，使用PASS、CLC-Pred 2.0、Pred-hERG 5.0和GUSAR等工具进行毒性预测（包括生物活性、细胞毒性、心脏毒性和大鼠急性毒性）。
**研究结果**
**3.1 ADMET分析**：从ChemFaces库中，经过ADME筛选，最终17个化合物通过所有毒性过滤；从MedChemExpress库中，经过多级筛选，最终2个化合物通过所有毒性过滤。共19个化合物进入对接研究。标准化合物vorasidenib和所有命中化合物均通过BBB渗透性预测（BOILED-egg分析显示位于卵黄区域），且均为P-糖蛋白非底物（PGP-），表明其具有中枢神经系统靶向递送的理论适用性。
**3.2 分子对接**：对接结果显示，对IDH1，I1Hit1（β-Eudesmol）结合能为?7.57 kcal/mol（Ki=2.84 μM），I1Hit2（Propranolol）为?6.82 kcal/mol（Ki=10.08 μM），均优于标准vorasidenib的?6.17 kcal/mol（Ki=30.15 μM）。对IDH2，I2Hit1（Isocoronarin D）结合能为?11.23 kcal/mol（Ki=5.82 nM），I2Hit2（Propranolol）为?7.87 kcal/mol（Ki=1.71 μM），均优于标准。关键相互作用残基包括IDH1的GLN277、SER280、TRP124、VAL255、MET259等，以及IDH2的VAL315、LEU320、TRP164、ILE168等。
**3.3 分子动力学模拟**：对IDH1复合物，I1Hit1和I1Hit2的平均RMSD值（0.40 nm和0.44 nm）低于标准（0.55 nm），RMSF值（0.16 nm和0.17 nm）也低于标准，Rg值（2.28-2.33 nm）显示良好紧凑性，SASA值（208.7和208.5 nm2）低于标准（211.9 nm2），氢键分析显示I1Hit2形成最多5个氢键。对IDH2复合物，I2Hit1平均RMSD为0.41 nm，RMSF为0.18 nm，Rg为2.31 nm，SASA为209.76 nm2，氢键分析显示I2Hit1在272-300 ns期间保持2个氢键。总体而言，命中化合物表现出更好的结构稳定性和动态特性。
**3.4 主成分分析（PCA）**：PCA结果显示I1Hit1和I1Hit2的复合物特征值低于标准，分布更集中，表明更小的构象变化和更高的稳定性。FEL分析显示，I1Hit1和I1Hit2的自由能分别为16.9和16.4 kJ/mol，与标准（16.4 kJ/mol）相当，但具有更窄的能垒盆地，表明形成更稳定的构象。对IDH2，I2Hit1和I2Hit2的自由能分别为16.3和15.9 kJ/mol，I2Hit1显示单一狭窄的能量最小值，表明构象稳定性高。
**3.5 结合自由能计算**：MM-PBSA计算结果显示，对IDH1，I1Hit1的总结合自由能为?20.67 ± 2.32 kcal/mol，显著优于标准的?14.96 ± 4.22 kcal/mol；I1Hit2为?11.77 ± 4.84 kcal/mol。对IDH2，I2Hit1的总结合自由能为?28.78 ± 2.67 kcal/mol，优于标准的?24.87 ± 3.55 kcal/mol；I2Hit2为?8.91 ± 4.74 kcal/mol。I1Hit1和I2Hit1均表现出增强的结合亲和力。
**3.6 毒性预测**：PASS预测显示I1Hit1和I2Hit1具有高抗肿瘤活性（Pa=0.792和0.960，Pi分别为0.013和0.004）和TP53表达增强潜力，均优于标准。CLC-Pred 2.0预测所有命中化合物对脑肿瘤细胞系（如U-87MG、C6、SF-539等）具有较高细胞毒性，不变精度（IAP）为0.8-0.9。Pred-hERG 5.0预测I2Hit1为非hERG阻断剂，I1Hit1为弱阻断剂，I1Hit2/I2Hit2为中等阻断剂。GUSAR预测所有命中化合物在大鼠急性毒性中均处于适用域内，毒性等级为3-5级（LD50值范围较广），显示低急性毒性潜力。
**讨论与结论**
讨论部分指出，IDH突变不仅通过2-HG产生直接影响肿瘤代谢，还通过表观遗传重编程、免疫调节和代谢可塑性等机制影响胶质瘤生物学。2-HG积累导致DNA和组蛋白高甲基化，沉默干扰素信号，降低PD-L1表达，同时IDH突变与免疫抑制微环境（髓源性抑制细胞、肿瘤相关巨噬细胞等）相互作用。此外，IDH突变常与TP53、ATRX、EGFR、MGMT等共突变，这些共突变影响治疗反应。尽管现有IDH抑制剂如ivosidenib和vorasidenib在临床试验中显示疗效（如INDIGO试验改善无进展生存期），但肿瘤异质性、耐药机制和BBB通透性等问题依然存在。本研究的计算流程（虚拟筛选、ADMET过滤、对接、MD模拟、PCA、MM-PBSA）识别出I1Hit1（β-Eudesmol）和I2Hit1（Isocoronarin D）作为针对突变IDH1和IDH2的潜在抑制剂，它们表现出优于标准vorasidenib的结合亲和力、稳定性和安全性预测。值得注意的是，I1Hit2/I2Hit2（Propranolol）是一种已获批的心血管药物，其作为IDH抑制剂的多靶点潜力具有药物重定位意义。研究结论强调，I1Hit1和I2Hit1是经过计算优先排序的候选药物，但需要进一步的实验验证（如酶抑制活性测定、2-HG定量、细胞系评估）来确认其治疗潜力并推进临床转化。
**结论翻译**：本研究采用整合计算方法识别针对GBM中突变IDH1和IDH2酶的潜在抑制剂。经过多轮虚拟筛选、ADMET过滤和多种分子分析，PubChem ID 91457（I1Hit1）针对IDH1和PubChem ID 71550939（I2Hit1）针对IDH2被确定为单独的潜在命中化合物。与标准参比化合物相比，这些化合物表现出更好的结合亲和力、良好的类药性和增强的结构稳定性。此外，分子对接、分子动力学模拟、PCA和MM-PBSA分析支持了这些复合物的稳定性和相互作用强度，表明它们具有抑制胞质和线粒体IDH亚型的潜力。所识别的命中化合物可能调节IDH相关的代谢改变，包括与GBM进展相关的氧化还原失衡、NADPH耗竭、氧化应激和表观遗传失调。总之，这些发现强调I1Hit1和I2Hit1是IDH靶向治疗GBM的计算优先候选药物，本研究所概述的整合计算框架为推进新型IDH靶向抑制剂走向临床前研究提供了可复制的平台。由于本研究纯属计算性质，未来的实验验证，包括酶抑制测定、2-HG定量和IDH突变GBM细胞系中的细胞评估，对于确认这些候选药物的治疗潜力并将其推进至临床转化至关重要。
