基于分子动力学模拟揭示多巴胺转运体与底物、可卡因及治疗药物的差异性结合机制

《International Journal of Biological Macromolecules》：Insight into the differential binding mechanisms of substrate, cocaine, and therapeutic drugs with dopamine transporter based on molecular dynamics simulations

　　

摘要 多巴胺转运体（DAT）是一种关键的膜蛋白，负责调节多巴胺能信号传导，同时也是成瘾药物和治疗性抑制剂的主要靶点。然而，内源性底物、成瘾物质和治疗性抑制剂差异化调控DAT构象的原子水平机制仍未完全阐明。本研究采用全原子分子动力学（MD）模拟，结合MM/PBSA计算、主成分分析（PCA）和自由能景观（FEL）分析，比较了多巴胺（DA）、可卡因、苯扎托品（BZT）和GBR12909（GBR）与DAT的结合模式。研究识别出三种不同的结合模式：DA主要以静电作用为驱动，但这种有利的贡献大部分被极性溶剂化罚分所抵消；可卡因和GBR依赖于范德华力和疏水作用；而BZT表现出双重驱动模式，同时具有强静电和范德华力贡献。DA通过Asp12介导的盐桥和多点氢键网络实现识别，而可卡因和GBR通过疏水/芳香残基嵌入结合，BZT则通过协同的离子锚定和疏水包裹结合。值得注意的是，Asp12发生功能角色转换——在DA和BZT结合时有利，而在可卡因和GBR结合时不利。门控距离、PCA和FEL分析进一步揭示了配体依赖的门控调控：BZT施加最强的构象约束，可卡因稳定外开状态，GBR封闭细胞外入口并伴随细胞内门控转换，而DA保留构象可调性。这些发现为理解DAT的差异化调控提供了系统的原子水平机制框架，并为新型竞争性抑制剂的开发提供了指导。引言 多巴胺（DA）是大脑中最丰富的儿茶酚胺能神经递质，也是一种被广泛研究的感觉神经递质[1], [2]。作为调节中枢神经系统（CNS）功能的核心分子，它不仅参与愉悦感和满足感的诱导，还在认知、运动、情绪和奖赏通路的信号传导中发挥关键的调节作用[3]。DA主要在多巴胺能神经元的末端合成并释放到突触间隙，在那里激活突触前和突触后膜上的多巴胺受体[4], [5]。通过这些受体，DA参与多种关键生理过程的信息调节，从而维持脑生理功能的稳态平衡。多巴胺转运体（DAT）[6], [7], [8]由SLC6A3基因编码，属于神经递质钠同向转运体（NSS）家族，是神经系统中高度保守的膜蛋白（图1）[9], [10], [11]。位于突触前膜上的DAT负责将DA从突触间隙重新摄取回突触前神经元，从而终止多巴胺能信号传导并维持DA稳态[12]。因此，DAT是多巴胺能神经传递的关键调节器。其功能依赖于钠离子和氯离子梯度来驱动底物的跨膜转运。异常的DAT水平与多种精神疾病密切相关，包括抑郁症、双相情感障碍、帕金森病（PD）和注意缺陷多动障碍（ADHD）[13], [14], [15]。由于其在调节DA神经传递时空动态中的关键作用，DAT已成为这些疾病的核心治疗靶点。DAT也是各种精神活性物质的重要靶点，包括可卡因、甲基苯丙胺和D-苯丙胺等非法精神兴奋剂，以及用于治疗相关疾病的治疗药物[16], [17], [18]。所有这些物质都会抑制DAT的DA转运活性，导致突触间隙DA浓度迅速升高。其中，非法精神兴奋剂可以诱发强烈的欣快感，通过奖赏反馈和强化效应促进物质滥用和成瘾，构成了一个全球性的公共卫生问题。作为一种典型的成瘾物质，长期使用可卡因不仅严重损害个体健康，还会加剧精神分裂症和焦虑障碍等疾病的发展。在机制上，可卡因特异性结合DAT底物结合口袋，从而阻断DA的重摄取。这导致突触间隙DA异常积累，持续刺激大脑奖赏回路，最终产生正性强化。然而，可卡因诱导DAT功能障碍的原子水平动态机制仍未得到充分认识。这一知识空白阻碍了有效药物治疗方法的开发。迄今为止，尚无FDA批准的可卡因使用障碍药物治疗方案，这仍然是一个紧迫的全球公共卫生挑战。近年来，研究人员致力于开发可卡因成瘾的治疗药物，并发现GBR12909（GBR）、苯扎托品（BZT）及其衍生物能够减轻可卡因诱导的运动刺激和奖赏效应，减少自我给药和强化学习行为，从而对可卡因成瘾具有潜在的治疗价值[12], [19], [20], [21]。本研究聚焦于底物DA、成瘾物质可卡因以及潜在的抗可卡因成瘾药物GBR和BZT。目前仍有两个关键问题尚未解决。首先，尽管GBR、BZT和可卡因都以DAT为靶点，但它们诱导了截然不同的构象锁定模式——一种"同靶点-不同效应"现象，缺乏原子水平的机制解释。其次，TM1a和TM6b的协同运动如何驱动DA跨膜转运仍然 elusive，限制了对DAT交替访问机制的完整理解。鉴于DAT的生理重要性和治疗相关性，已开展了大量实验研究[12], [22], [23], [24], [25], [26], [27], [28], [29]。一些研究相继解析了果蝇DAT（dDAT）[30]和人类DAT[12]的静态结构，鉴定了底物结合位点、离子结合位点和配体结合区域，揭示了DAT在转运过程中的三种构象（外开、封闭和内开），并阐明了DA重摄取过程中DAT的底物识别和构象转换特征。同时，不同类型DAT靶向药物的结合模式差异也得到了阐明，为相关药物的开发提供了基础依据。然而，静态结构只能捕获DAT转运周期中的中间状态，无法完全描绘底物转运、可卡因作用以及GBR和BZT抑制过程的动态全貌或关键细节。这种静态结构与动态功能之间的脱节已成为阐明可卡因作用机制和合理设计DAT靶向药物的关键瓶颈。分子动力学（MD）模拟作为生物物理和药物发现领域的核心计算工具，可以在原子水平捕获生物分子系统的动态行为，有效弥补静态结构研究的局限性，弥合静态结构与动态功能之间的差距[31]。与静态结构技术不同，MD模拟可以实时追踪蛋白质构象波动、分子间相互作用（如氢键和疏水接触）的形成与断裂，以及脂质膜环境对蛋白-配体相互作用的影响，从而克服传统实验方法在捕获微秒级构象转换方面的固有局限。近年来，一些研究已应用MD模拟探索DAT构象动力学[32], [33], [34], [35], [36]。然而，大多数研究仅关注单一配体与DAT的相互作用，缺乏对DA、GBR、BZT和可卡因结合机制的系统比较。此外，尚缺乏解释底物识别特异性和抑制剂诱导构象锁定的原子水平动态机制的系统理论框架。在静态结构知识库的基础上，本研究采用全原子MD模拟结合结合自由能计算、主成分分析和自由能景观分析，系统研究了DAT的构象调控机制。以DA（底物）、可卡因（成瘾物质）以及GBR和BZT（潜在治疗药物）作为比较体系，本研究旨在：（i）阐明这些配体对DAT构象动力学的差异化影响，包括细胞内和细胞外门的开合特征；（ii）解析可卡因干扰DAT功能的关键动态机制；（iii）阐明GBR、BZT和可卡因之间"同靶点-不同效应"现象的分子基础。本研究的结果将为合理设计DAT靶向的低成瘾性治疗药物提供原子水平理论支持，为开发有效的抗成瘾干预措施提供科学依据，并促进神经精神疾病治疗领域的进展。需要注意的是，本研究所使用的五个冷冻电镜结构作为起始模型，是最近由Li等（2024）[12]和Nielsen等（2024）[37]报道的。仅从那些静态结构出发，那些作者建立了配体依赖的构象偏好定性分类：DA处于封闭状态，可卡因处于外开状态，GBR12909和苯扎托品处于内开状态。本基于MD的研究与该分类完全一致。本研究并未重复该分类，而是提供了静态冷冻电镜结构无法提供的四项贡献：（i）定量结合自由能分解——明确分解为范德华力（ΔEvdW）、静电（ΔEele）、极性溶剂化（ΔGPB）和非极性溶剂化（ΔGSA）组分；（ii）残基水平的Asp12角色转换——在DATDA和DATBZT中通过盐桥有利，在DATCOC和DATGBR中通过羧基-配体竞争不利；（iii）时间分辨门控距离动力学——D12-F253（细胞外）和F9-F259（细胞内），揭示了四种配体之间的不同门控模式；（iv）自由能景观采样——沿PC1/PC2和Rg/RMSD重建连续自由能面，揭示了超越五个离散结构点的亚盆连通性和相对深度。章节片段 结构准备 本研究所使用的多巴胺转运体（DAT）三维结构均来自蛋白质数据银行（PDB），包括五种不同的结合状态：无底物状态（DATapo，PDB ID: 8Y2C）、多巴胺结合状态（DATDA，PDB ID: 8Y2D）、可卡因结合状态（DATCOC，PDB ID: 9EO4）、苯扎托品结合状态（DATBZT，PDB ID: 8Y2E）和GBR12909结合状态（DATGBR，PDB ID: 8Y2F），这些结构作为后续MD模拟的初始基础[12], [37]。值得注意的是，不同配体结合对DAT整体构象稳定性的影响 不同配体状态下DAT的构象稳定性通过计算400 ns MD模拟期间Cα原子的均方根偏差（RMSD）进行评估（图2；重复数据见附图S2）。如图2A所示，DATapo系统在前50 ns内迅速稳定，随后维持较低的RMSD水平。相比之下，DATDA的RMSD在模拟初期急剧上升，经历了一个短暂的过渡期（50–100 ns），然后达到平台期。结论 在本研究中，采用全原子分子动力学模拟结合MM/PBSA计算、主成分分析和自由能景观分析，系统比较了内源性底物（DA）、成瘾药物（可卡因）和两种治疗性抑制剂（GBR和BZT）如何差异化调控DAT构象。综合结果描绘了一个连贯的机制框架，其中四种配体沿单一谱系占据四个不同的位置。作者贡献声明 吴建华：论文撰写——初稿，可视化，软件，方法学，调查研究，形式分析，数据整理，概念化。王波：论文撰写——初稿，验证，监督，项目管理，概念化。李宗航：论文撰写——审阅与编辑，资金获取，形式分析，数据整理。刘涵钰：论文撰写——初稿，验证，形式分析，数据整理。张婷：论文撰写——审阅与编辑，可视化，调查研究，资金声明 利益冲突声明 作者声明他们没有已知的竞争性经济利益或个人关系，这些利益或关系可能影响本文所报告的工作。致谢 本研究得到辽宁省科学技术厅博士启动基金计划（2025BS0257、2025BS0255和2025BS0250）、中央高校基本科研业务费专项资金（D2026056）、辽宁省科学技术厅重点研发计划（民生科技）（2024JH2/102500083）以及校级研究生教育教学改革项目（2025YJG005）的支持。吴建华 | 王波 | 李宗航 | 刘涵钰 | 张婷 | 郭嘉明