基于药效团筛选与分子模拟的BRD4抑制剂计算发现及其在神经母细胞瘤治疗中的应用

《Scientific Reports》：Computational discovery of BRD4 inhibitors for neuroblastoma therapy using pharmacophore screening and molecular simulations

本研究针对神经母细胞瘤治疗中BRD4靶点抑制剂的开发难题，采用整合药效信息学策略，通过药效团虚拟筛选、分子对接、ADMET分析和分子动力学模拟，从五大数据库中筛选出10个高亲和力候选化合物（结合能-9.623至-8.894 kcal/mol）。其中CSC101378549和MCULE-2,532,030,345展现出优异的结合稳定性和药代动力学特性，为BRD4靶向治疗提供了新型先导化合物。

　　

神经母细胞瘤作为儿童最常见的颅外实体肿瘤，其治疗面临巨大挑战。这种起源于未成熟神经细胞的恶性肿瘤具有高度异质性，从可自行消退的惰性肿瘤到快速转移的侵袭性类型，临床表现千差万别。近年来，表观遗传调控因子BRD4（Bromodomain-containing protein 4）被证实在该疾病发生发展中起关键作用，成为极具潜力的治疗靶点。

BRD4作为BET蛋白家族成员，通过识别组蛋白乙酰化标记调控基因表达，参与细胞周期进程、凋亡抵抗等致癌过程。特别值得注意的是，BRD4在神经母细胞瘤中常出现异常活化，驱动肿瘤细胞增殖和存活。然而，现有BET抑制剂如(+)-JQ1存在选择性差、毒副作用大等问题，严重限制其临床应用。开发高选择性、强效的BRD4特异性抑制剂成为当前研究的热点难点。

针对这一需求，研究人员开展了一项整合计算药物发现研究。该研究创新性地采用多层次计算生物学方法，旨在发现新型BRD4抑制剂候选化合物。研究团队首先基于BRD4蛋白晶体结构（PDB ID: 4BJX）中共结晶配体73B的相互作用特征，构建包含疏水基团、芳香环和氢键受体等关键药效特征的查询模型。

通过药效团虚拟筛选技术，研究人员对CHEMBL、Enamine、MCULE、Chemspace和ChemDiv五大化合物数据库进行系统性筛选，重点关注符合Lipinski规则（氢键受体<10，分子量<500，氢键供体<5，logP<5）的类药分子。从初始筛选获得的1089个命中化合物中，采用分子对接技术进一步评估其与BRD4活性位点的结合能力。

关键技术方法包括：基于Pharmit服务器的药效团建模、Schrodinger Maestro平台的分子对接（Glide SP模式）、ADMET性质预测（QikProp工具）、结合自由能计算（MM-GBSA方法）以及200纳秒分子动力学模拟（Desmond软件）。所有计算均采用OPLS_2005力场，并在生理条件（pH 7.0，300 K，1 atm）下进行验证。

Pharmacophore modelling and virtual screening

研究基于共结晶配体73B的相互作用特征，构建了包含2个疏水基团、2个芳香环和1个氢键受体的药效团模型。虚拟筛选从五大数据库中获得1089个符合标准的候选化合物，其中MCULE数据库贡献最多命中分子。

Structure-based molecular docking

分子对接结果显示，前十名化合物与BRD4的结合能介于-9.623至-8.894 kcal/mol之间，显著优于已知抑制剂RVX-208和MS436。关键相互作用分析表明，这些化合物主要通过氢键与Pro82、Gln85、Tyr97、Asn140等残基结合，同时形成大量疏水相互作用稳定结合构象。

Prediction of biological activity,availability,and drug-likeness of top ten hits

生物活性预测显示，所有候选化合物在GPCR配体、酶抑制、激酶抑制等关键指标均表现良好（活性评分0.68至-0.48）。其中MCULE-2,532,030,345和CSC101378549表现出最优的类药性，生物利用度评分达0.55，具备口服给药潜力。

ADME and toxicological properties of top hits

ADMET分析表明，CSC101378549和MCULE-2,532,030,345完全符合药物开发标准：QPlogPo/w（-2.0-6.5）反映适宜脂溶性，QPlogHERG（<-5）提示无心脏毒性风险，QPPCaco（>500）显示良好肠道渗透性，QPlogBB（-3.0-1.2）表明可控血脑屏障穿透能力。

Binding free energy calculations of docked complexes

MM-GBSA结合自由能计算证实，CSC101378549和MCULE-2,532,030,345与BRD4的结合能分别为-82.851 kcal/mol和-86.300 kcal/mol，主要贡献来自疏水相互作用（ΔG_bindLipo）和范德华力（ΔG_bindvdW），表明结合过程受疏水作用和形状互补性主导。

Binding pose analysis of potent complexes

结合模式比对显示，两个先导化合物与共结晶配体73B在活性位点呈现高度重叠，验证了筛选策略的准确性。这种相似的空间布局暗示它们可能通过类似机制抑制BRD4功能。

Structure stability analysis of potential complexes

200纳秒分子动力学模拟证实了复合物的稳定性。CSC101378549-BRD4体系的RMSD值稳定在2.1-3.2 ?范围内，而MCULE-2,532,030,345复合物RMSD波动更小（1.6-2.8 ?）。RMSF分析显示活性位点残基柔性变化可控，二级结构元素保持稳定。

氢键分析表明，CSC101378549与Pro82、Gln85、Tyr97、Asn135维持2-4个稳定氢键，MCULE-2,532,030,345与Gln85、Tyr97、Asn135等残基形成2-3个氢键。溶剂可及表面积（SASA）和回转半径（Rg）分析进一步验证了蛋白质结构的紧凑性和稳定性。

Principal component analysis

主成分分析显示，两个复合物在前五个特征向量中的累积方差分别达66%和77.4%，其中PC1主导运动（21.61%和39.18%）。自由能面分析表明体系主要处于局部能量极小值区域，证实结合构象的热力学稳定性。

Maps of positive and negative correlation motions

动态交叉相关矩阵（DCCM）揭示残基20-45和80-110区域存在显著协同运动，活性位点残基表现出高度相关性，进一步支持了配体结合诱导的构象稳定效应。

本研究成功建立了一套完整的BRD4抑制剂计算发现流程。通过多层次虚拟筛选策略，从大规模化合物库中识别出具有高结合亲和力和良好类药性的先导化合物。特别值得注意的是，CSC101378549和MCULE-2,532,030,345展现出优于已知抑制剂的结合特性，其与BRD4活性位点的相互作用模式与共结晶配体高度一致。

分子动力学模拟结果表明，这两个先导化合物能够与BRD4形成稳定的复合物结构，关键相互作用在200纳秒模拟时间内保持稳定。结合自由能计算进一步证实了结合过程的热力学可行性。ADMET性质预测显示这些化合物具备良好的药代动力学特征和可接受的安全性风险。

该研究的创新性在于将药效团筛选、分子对接、ADMET预测和分子动力学模拟有机整合，形成了一套高效的计算药物发现范式。这不仅为神经母细胞瘤的靶向治疗提供了新的候选化合物，也为表观遗传靶点的药物开发提供了重要方法论参考。未来通过实验验证这些先导化合物的生物学活性，将有望推动新型BRD4抑制剂的临床前开发进程。