学术研究报告：多巴胺D3受体拮抗剂揭示胺能GPCRs中的隐秘口袋
一、 研究团队、发表期刊与时间
本研究由来自Universitat Pompeu Fabra、Acellera、Institució Catalana de Recerca i Estudis Avançats (ICREA)以及Pfizer Worldwide Research and Development多个研究机构的科研人员共同完成。主要作者包括Noelia Ferruz, Stefan Doerr, Michelle A. Vanase-Frawley, Yaozhong Zou, Xiaomin Chen, Eric S. Marr, Robin T. Nelson, Bethany L. Kormos, Travis T. Wager, Xinjun Hou, Anabella Villalobos, Simone Sciabola与Gianni De Fabritiis。该研究于2018年1月17日在线发表在期刊 Scientific Reports 上（卷8，文章号897）。
二、 学术背景与研究目的
本研究隶属于结构生物学与计算药物设计交叉领域，聚焦于G蛋白偶联受体（GPCRs）。GPCRs是最大的膜蛋白家族之一，是约34%上市药物的作用靶点，在细胞信号转导中扮演核心角色。胺能GPCRs是其中一类重要的亚家族，多巴胺D3受体（D3R）因其在药物成瘾等物质使用障碍中的潜在治疗作用而备受关注。
传统的基于结构的药物设计（SBDD）严重依赖静态的X射线晶体结构。然而，GPCRs是高度动态的蛋白质，其功能依赖于固有的柔性。使用静态结构来理解配体结合，尤其是在构象灵活性高的系统中，可能会产生误导。本研究旨在解决一个具体矛盾：将一种名为PF-4363467（文中化合物1）的D3R双重拮抗剂，对接至已知的、与拮抗剂eticlopride（化合物2）共结晶的D3R静态晶体结构（PDB ID: 3PBL）时，所得的对接构象与定点突变实验获得的数据不一致。
因此，本研究的主要目标是：1）通过超越静态对接的方法，阐明PF-4363467与D3R的真实结合模式；2）探究这种结合模式是否揭示了D3R或更广泛的胺能GPCRs中未被认知的构象状态，从而为药物设计提供新思路。
三、 详细研究流程与方法
本研究是一个结合了实验生物化学、计算化学与生物物理学模拟的综合性工作，流程严谨，环环相扣。
1. 突变结合实验与刚性分子对接： 首先，研究团队对D3R进行了系统的定点突变。他们在D3R的正构结合位点（OBS）及其附近选择了9个关键残基，构建了12个不同的D3R突变体。随后，使用放射性配体结合实验，精确测定了四种D3R拮抗剂（PF-4363467 (1)、eticlopride (2)、haloperidol (3) 和GSK598809 (4)）对野生型（WT）及各个突变体D3R的抑制常数（Ki）。通过计算突变体与野生型Ki的比值（Fold-shift），量化了每个突变对配体亲和力的影响。
同时，研究人员使用分子对接软件Glide，将四种化合物刚性对接至eticlopride结合的D3R晶体结构（3PBL）中。对接网格覆盖了整个正构结合口袋。对于eticlopride (2)，对接成功重现了其晶体结构中的结合模式。对于haloperidol (3) 和GSK598809 (4)，对接产生的构象也与已知的文献报道及本次突变数据基本吻合。然而，对于PF-4363467 (1)，对接产生的最佳构象显示其吗啉环上的碱性氮原子与保守残基D110³.³²形成盐桥，但其4-异丙基苯基伸向了由胞外环构成的区域（胞外前庭），而OBS内部由螺旋III、V和VI围成的大量疏水区域未被占据。这一对接构象无法解释两个关键的突变实验现象：I183ECL2F突变导致亲和力下降14倍，V189⁵.³⁹A突变导致亲和力下降8倍。在对接构象中，化合物1距离这两个突变位点分别超过6Å和9Å，其结合不应受到如此显著的影响。这一矛盾表明，静态对接可能未能捕捉到PF-4363467的真实结合模式。
2. 大规模高通量分子动力学模拟与马尔可夫状态模型分析： 为了突破静态结构的局限，研究团队采用了大规模、无偏的分子动力学（MD）模拟来研究PF-4363467与D3R的结合过程。这是一个关键且具有创新性的方法学应用。 * 系统构建与模拟设置： 模拟体系基于D3R晶体结构（3PBL）构建，并进行了必要的修正（如将L119突变回野生型的Trp）。配体PF-4363467被随机放置在受体蛋白的胞外区域。体系使用CHARMM36力场进行参数化，配体参数由基于HTMD软件的工具生成。体系被置于POPC脂质双分子层中，并用水和离子进行溶剂化。 * 模拟执行与自适应采样： 模拟使用ACEMD软件在分布式计算项目GPUGRID.net上进行。研究采用了无监督的自适应采样协议，这是一种高效的构象空间探索方法。模拟分多个“纪元”进行：每个纪元结束后，利用马尔可夫状态模型（MSM）对累积的模拟轨迹进行自动分析，识别出未充分采样或具有动力学重要性的构象，并将其作为下一个纪元模拟的起始点。这种方法能够智能地引导采样，更有效地探索复杂的结合路径和构象变化。 * 数据分析与结合路径重构： 总共累积了约700微秒的模拟数据（其中220微秒用于最终分析）。研究人员使用时间结构独立成分分析（tICA）和聚类方法，将高维模拟数据投影并离散化为多个微观状态。然后，通过PCCA+算法将这些微观状态归并为几个宏观的动力学状态，这些状态代表了结合路径上的不同阶段（如：配体在溶液中的“本体”状态、在胞外前庭的初始结合状态、在OBS内的中间状态、最终结合状态等）。通过构建马尔可夫状态模型，可以定量计算状态间的转换速率、平均首次通过时间（MFPT）以及结合速率常数（kon）。
3. 结构分析与比较： 从MSM确定的最终结合状态（State 4）中，选取了最具代表性的构象进行详细的结构分析。研究人员将该预测构象与D3R/eticlopride晶体结构以及其他已公开的胺能GPCR晶体结构进行了系统的比较，重点分析了关键残基的构象变化、保守相互作用的维持情况以及可能的新颖结构特征。
四、 主要研究结果
1. 突变与对接结果凸显矛盾： 如流程部分所述，对化合物2、3、4的突变数据均能用其刚性对接构象合理解释。唯独化合物1的对接构象与I183ECL2F和V189⁵.³⁹A两个位点的突变数据严重不符，这强烈暗示存在一种未被晶体结构捕捉的替代结合模式。
2. MD模拟揭示全新的结合模式与隐秘口袋： 大规模MD模拟成功地观察到了PF-4363467多次自发结合到D3R的过程。MSM分析揭示了其结合路径：配体首先识别胞外前庭（State 2），然后经历一个限速步骤（约7.4微秒），即其带正电的胺基与D110³.³²形成盐桥，并伴随构象重排，最终到达稳定的结合状态（State 4）。 * 最终结合构象细节： 在MSM确定的最终结合构象中，PF-4363467的结合方向与对接预测截然不同。其吗啉环仍与D110³.³²形成保守的盐桥，但分子的其余部分深深地插入到OBS的疏水核心中。关键的发现是，其4-异丙基苯基部分推开了位于螺旋V和VI之间的两个高度保守的芳香族残基：F197⁵.⁴⁷和F346⁶.⁵²。特别是F346⁶.⁵²，其侧链二面角χ1旋转到了约-174.4°，这是一个在已解析的39个胺能GPCR晶体结构中从未观察到的构象（已知构象集中在-90.8°和70.2°附近）。这一位移导致在F338⁶.⁴⁴、W342⁶.⁴⁸、L343⁶.⁴⁹、F345⁶.⁵¹和F346⁶.⁵²之间形成了一个隐秘口袋，该口袋被PF-4363467的4-异丙基苯基所占据。 * 合理解释突变数据： 这一新结合构象完美解释了突变实验结果。在新的构象中，V189⁵.³⁹的侧链与化合物1的磺酰胺基团紧密接触，将其突变为侧链更小的丙氨酸（V189A）会削弱有利的疏水相互作用，导致亲和力下降。同时，I183ECL2的侧链与化合物1的2-甲基苯基以及V189⁵.³⁹的侧链紧密堆积。将其突变为体积更大的苯丙氨酸（I183F）会产生空间位阻，与配体或邻近的V189⁵.³⁹发生冲突，从而显著降低结合亲和力。这与实验观察到的14倍亲和力下降完全一致。
3. 结合路径与动力学参数： MSM分析计算出PF-4363467结合到D3R的平均首次通过时间约为7.0微秒，对应的结合速率常数kon约为12 × 10⁶ M⁻¹ s⁻¹，与其他GPCR药物结合的研究结果相符。模拟揭示了两种主要的结合路径，其中约70%的结合事件通过一个两步路径（胞外识别→限速重排）完成。
4. 受体整体构象保持非活性状态： 尽管局部形成了隐秘口袋，但预测的D3R-PF-4363467复合物整体构象仍保持非活性状态。关键的“离子锁”（R³.⁵⁰与E⁶.³⁰之间的盐桥）以及NPxxY基序的构象均与D3R/eticlopride晶体结构相似，表明PF-4363467是一个拮抗剂。
五、 研究结论与意义
本研究得出结论：PF-4363467通过与静态晶体结构不同的模式结合D3R，其结合诱导了一个由F346⁶.⁵²构象变化形成的隐秘口袋。这一结合模式得到了突变实验数据的验证，并只能通过纳入蛋白质动态信息的大规模分子动力学模拟来发现。
本研究具有重要的科学价值与应用价值： 1. 强调了GPCR动态性在药物设计中的关键作用： 研究证明，仅依赖静态晶体结构进行药物设计可能遗漏重要的结合模式，尤其是对于那些能够诱导受体产生非经典构象的配体。蛋白质的固有柔性必须被纳入考虑。 2. 揭示了胺能GPCR中一个新的可靶向构象状态： 发现的由F6.52位移形成的隐秘口袋，为设计新型、高选择性的D3R拮抗剂（乃至其他胺能GPCR配体）提供了全新的思路和结构基础。这个口袋在已知晶体结构中不存在，但通过配体结合可以被诱导出来。 3. 展示了计算方法的强大能力： 研究成功地将大规模MD模拟与MSM分析相结合，不仅重现了已知配体的结合，还解决了一个实验与静态模型之间的矛盾，预测并验证了一个全新的结合模式。这为在难以获得晶体结构的GPCR药物发现项目中，利用计算模拟指导先导化合物优化和虚拟筛选提供了强有力的范例。
六、 研究亮点
七、 其他有价值的内容
研究团队还以已知结合模式的eticlopride (2) 作为对照，验证了其MD模拟流程能够重现实验确定的结合构象，增强了整个方法学的可信度。此外，文章在讨论部分提到，在另一项关于氯氮平和氟哌啶醇结合D2R/D3R的大规模MD研究（Thomas et al., 2016）的补充材料中，发现了与本研究F346⁶.⁵²构象相似的氯氮平结合模式，这为本文发现的隐秘口袋构象提供了独立的佐证，表明这种构象可能在胺能GPCR配体结合中具有一定的普适性。