基于生物信息学与网络药理学表征呫诺美林对槟榔依赖的治疗干预潜力

《Technology and Health Care》：Bioinformatics and network pharmacological characterization of xanomeline for therapeutic intervention in areca nut dependence

　　

背景：槟榔作为一种已知致癌物和成瘾物质，仍在全球范围内被消费。目前针对槟榔成瘾的戒断策略存在空白。呫诺美林（Xanomeline）作为槟榔碱（arecoline）的合成衍生物，近年来在治疗慢性神经疾病方面受到关注。目的：开展一项槟榔碱与呫诺美林之间的计算机模拟（in-silico）比较研究，以识别可能揭示未来干预策略开发的相互作用。方法：研究人员采用了一种整合性计算机模拟方法，包括理化性质分析、毒性预测、分子对接（MD）、分子动力学模拟（MDS）以及网络药理学，对槟榔碱和呫诺美林进行比较。对接对象为与胆碱能信号传导和成瘾相关的毒蕈碱受体（CHRM1–5）和烟碱受体（CHRNA3/4/7, CHRNB2/4）亚型。MDS用于评估蛋白质-配体相互作用。结果：呫诺美林在大多数胆碱能靶点上表现出更优的结合亲和力，特别是对CHRM1（–6.2?kcal/mol），而槟榔碱为（–4.4?kcal/mol）。两种配体均表现出稳定的RMSD轨迹，尽管呫诺美林复合物显示出相对较低的RMSF波动。呫诺美林的血浆蛋白结合率（PPB）（61.72%）高于槟榔碱（12.740%），血脑屏障（BBB）穿透概率（0.992?cm/s）也高于槟榔碱（0.89?cm/s）。AMES毒性预测显示槟榔碱具有中等预测致突变潜能（0.625），而呫诺美林则具有低预测致突变潜能（0.024）。结论：在这项计算机模拟研究的局限性内，呫诺美林在与成瘾相关的胆碱能受体上表现出更强的结合亲和力、在MDS中更高的稳定性、较低的预测致突变潜能以及优于槟榔碱的中枢神经系统分布特性。这些发现支持了呫诺美林调节与槟榔依赖和成瘾相关的胆碱能通路的潜力，并值得进一步开展实验验证，以将其作为一种有前景的药理学干预策略。

研究背景与问题

槟榔是全球仅次于咖啡因、烟草和酒精的最常消费的精神活性物质之一，尤其在亚太地区流行。尽管国际癌症研究机构（IARC）将其列为I类致癌物，但由于文化信仰、社会接受度、感知的健康益处以及无意识形成的成瘾/依赖性，其使用仍在持续。长期使用影响多个器官系统，显著增加口腔潜在恶性病变和口咽癌的风险。近年来，多样化产品（如果味槟榔糖）的出现使儿童群体中的使用率上升，甚至在印度孕妇中也存在高达16%的使用率，导致死产、低出生体重、早产、流产和新生儿戒断综合征等不良胎儿结局。然而，与受世界卫生组织《烟草控制框架公约》（FCTC）监管的烟草不同，槟榔缺乏全球性的监管措施，公众对其危害认识不足。目前，针对槟榔成瘾的戒断策略研究十分匮乏，尚无获批的特异性药物疗法。槟榔的主要活性成分槟榔碱（arecoline）是一种非选择性毒蕈碱乙酰胆碱受体（mAChRs）的部分激动剂，也能作用于与成瘾相关的烟碱乙酰胆碱受体（nAChRs）。研究表明，M₅受体抑制可能有助于治疗药物滥用，而M₁/M₄偏向性激动剂呫诺美林（xanomeline）可减弱可卡因的强化和辨别效应。因此，研究人员提出，作为槟榔碱的结构类似物且具有中枢活性的神经精神药物呫诺美林，可能成为治疗槟榔成瘾/依赖的干预策略。

研究方法概述

本研究采用整合性计算机模拟（in-silico）方法，主要包括：利用ADMETlab2.0平台进行理化性质和ADMET预测；使用GUSAR数据库进行急性毒性预测；通过网络药理学（SwissTargetPrediction、SuperPred、GeneCards、STRING数据库）挖掘共同靶点；运用SWISS-MODEL进行蛋白质同源建模并通过PROCHECK和Ramachandran图验证结构质量；采用AutoDock 4.2进行分子对接（MD）；利用Schr?dinger套件中的Desmond软件包进行100 ns的分子动力学模拟（MDS）。

研究结果

3.1 理化性质与ADMET分析

两种化合物均符合Lipinski“五规则”。呫诺美林的氢键受体数（4）多于槟榔碱（3），总极性表面积（TPSA）为38.25 ?2，高于槟榔碱的29.54 ?2。呫诺美林的logP值（4.0）显著高于槟榔碱（0.151），logS值（-4.088）低于槟榔碱（0.623），表明其亲脂性更强。ADMET预测显示，呫诺美林的血浆蛋白结合率（PPB）（61.72%）远高于槟榔碱（12.74%），血脑屏障（BBB）穿透概率（0.992 cm/s）也高于槟榔碱（0.89 cm/s），体积分布（VD）（2.473 L/kg）更大。代谢预测表明两者均为CYP3A4和CYP2D6底物，但呫诺美林对CYP2C19（0.836）和CYP2D6（0.938）的抑制潜力更高。毒性预测中，槟榔碱AMES致突变性评分较高（0.625），而呫诺美林较低（0.024）；但呫诺美林的人肝毒性评分（0.901）高于槟榔碱（0.224）。

3.2 急性毒性预测

在所有给药途径中，槟榔碱的半数致死剂量（LD₅₀）均低于呫诺美林。例如，静脉注射途径槟榔碱LD₅₀为19.5 mg/kg，呫诺美林为45.58 mg/kg；口服途径槟榔碱为361.4 mg/kg，呫诺美林为393.3 mg/kg。两种化合物均在静脉注射途径显示出最低的LD₅₀值，表明其急性毒性最高。

3.3 数据挖掘、网络药理学与靶点选择

从GeneCards数据库获得1597个与“成瘾”和“毒性遗传”相关的基因。通过SwissTargetPrediction和SuperPred预测，槟榔碱有51个靶基因，其中24个与GeneCards数据库重叠，包括所有毒蕈碱受体（CHRM1-5）和选定的烟碱受体（CHRNA7, CHRNA3, CHRNA4, CHRNB2, CHRNB4）。呫诺美林有65个预测靶基因，其中28个与GeneCards数据库重叠，但不包含上述烟碱受体。值得注意的是，呫诺美林的靶点包括大麻素受体1（CB1r）和2（CB2r）。

3.4 蛋白质同源建模

通过SWISS-MODEL生成的同源模型经PROCHECK和Ramachandran图验证，大部分残基位于最适区域，表明结构稳定性和准确性良好，其中CHRM5在最适区域的残基比例最高（95.1%）。

3.5 分子对接

槟榔碱对所有靶标的结合能约为-4 kcal/mol，其中对CHRNA4的结合能最优（-4.7 kcal/mol），但未形成氢键；对CHRM1的结合能为-4.4 kcal/mol，形成了四个氢键。呫诺美林对CHRM1展现出最佳结合能（-6.2 kcal/mol），并形成了三个氢键（Tyr B:82, Tyr B:85, Glu B:397），表明其与该靶点的相互作用最强且最稳定。

3.6 分子动力学模拟

100 ns的MDS结果显示，槟榔碱-CHRM1和呫诺美林-CHRM1复合物的RMSD值均在2–3 ?范围内，并在75 ns后趋于稳定。RMSF分析显示，呫诺美林-CHRM1复合物的波动较小。氢键分析表明，呫诺美林通过残基Tyr82在整个模拟期间约22%的时间内维持了与CHRM1的相互作用，而槟榔碱的相关作用较弱。

讨论与结论

讨论部分指出，本研究的计算机模拟分析首次比较了槟榔碱与其合成类似物呫诺美林在槟榔成瘾方面的潜力。呫诺美林更高的亲脂性（logP 4.0）有利于其穿过细胞膜和血脑屏障，这与其中枢神经系统活性一致。其较高的血浆蛋白结合率和更大的体积分布预示着更广泛的组织分布和较长的循环时间。尽管呫诺美林在人肝毒性预测上得分较高，但其较低的致突变潜力和更高的LD₅₀值显示出比槟榔碱更有利的安全性特征。分子对接和动力学模拟共同证实，呫诺美林对CHRM1具有更强的结合亲和力和更稳定的相互作用。此外，呫诺美林靶向大麻素受体（尤其是CB2r）的特性可能为其在物质使用障碍中的作用提供了额外机制通路。研究人员假设，呫诺美林的中枢毒蕈碱受体激动活性，特别是对M₁/M₄受体的偏好性，可能通过与槟榔碱重叠的受体靶向来调节奖赏相关的胆碱能信号传导，从而减轻槟榔依赖。研究结论为：尽管两种化合物均具有良好的类药物性质，但在药代动力学和毒理学特征上存在差异。预测显示槟榔碱具有较低的LD₅₀值和中等致突变潜力，而呫诺美林具有更高的亲脂性、血浆蛋白结合率、血脑屏障穿透概率和体积分布。呫诺美林与CHRM1相对更好的受体结合能力和稳定性，以及其对大麻素受体的作用，表明其在槟榔依赖中可能具有治疗相关性。