别构抑制剂G6PDi-1与人G6PD结合全景图的绘制

《The Journal of Physical Chemistry Letters》：Mapping the Binding Landscape of Allosteric Inhibitor G6PDi-1 on Human G6PD

　　

研究人员通过抑制葡萄糖-6-磷酸脱氢酶（G6PD)靶向氧化型磷酸戊糖途径(PPP)是一种有前景的抗癌策略，但目前尚无临床可用的抑制剂。关键局限在于缺乏关于别构抑制剂和非竞争性抑制剂如何扰动酶活性的分子水平认识，从而限制了理性的优化。本研究中，研究人员结合生化实验、分子动力学(MD)模拟、马尔可夫状态模型(MSM)和MM/PBSA结合自由能计算方法，对一种强效、可逆、非竞争性G6PD抑制剂G6PDi-1的作用机制进行了探究。实验上，G6PDi-1降低了人肝母细胞瘤HepG2细胞中活性二聚体的浓度，并增加了无活性单体形式的占比，将抑制作用与寡聚化的破坏联系起来。计算研究方面，研究人员在酶表面鉴定了多个结合位点，其中G6PDi-1优先结合于二聚体界面，通过空间位阻效应阻断寡聚化。重要的是，额外的远端位点表现出增强的残基间相关性，提示存在二次别构效应，促使G6PD向不利于寡聚化的单体状态转变。这些发现为基于结构的G6PD抑制策略的改进开发提供了分子基础，以适用于癌症治疗。

本研究围绕人葡萄糖-6-磷酸脱氢酶(G6PD)的别构抑制剂G6PDi-1展开，旨在揭示其分子作用机制，为癌症治疗中G6PD靶向药物的开发提供理论基础。G6PD是氧化型磷酸戊糖途径(PPP)的限速酶，该途径生成核糖-5-磷酸以供核苷酸合成，并产生NADPH用于还原性生物合成及抗氧化防御。G6PD的表达与活性在多种肿瘤中上调，可促进细胞增殖、转移、代谢重编程等癌症特征。尽管靶向PPP具有治疗潜力，目前尚无G6PD抑制剂进入临床使用，开发高质量抑制剂仍面临挑战。人G6PD以二聚体或四聚体形式存在时具有催化活性，单体形式则无活性。此前虽有多种甾体和非甾体小分子被报道为G6PD抑制剂，但均因靶点特异性不足而难以临床转化。Rabinowitz及其同事发现的G6PDi-1是一种可逆、非竞争性抑制剂，虽知其具有别构作用且能破坏G6PD寡聚化，但具体分子机制不明。因此，开展本研究以深入解析G6PDi-1的作用机制，对于理性设计新型G6PD抑制剂具有重要意义。

本研究用到的主要技术方法包括：在人肝母细胞瘤HepG2细胞系中进行生化实验检测G6PDi-1对G6PD寡聚状态的影响；利用Schrodinger软件的SiteMap工具进行结合口袋预测，Glide工具进行分子对接，MM/GBSA方法进行结合自由能重打分；开展全原子分子动力学(MD)模拟（总计约120微秒模拟时长）；构建隐马尔可夫状态模型(HMSM)进行结合动力学分析；采用MM/PBSA方法进行结合自由能计算；运用归一化线性互信息(nLMI)分析残基间相关性变化；通过PLUMED包计算距离均方根偏差(DRMSD)评估二聚体界面构象分布。

研究背景与生化实验验证

G6PDi-1对G6PD寡聚化的影响：研究人员首先在人肝母细胞瘤HepG2细胞中验证了G6PDi-1的生物学效应。结果显示，G6PDi-1处理降低了G6PD二聚体形式的占比，同时增加了单体形式的占比。该实验将G6PDi-1的作用效应拓展至更具疾病相关性的肝细胞背景下，但无法区分二聚体减少是由于二聚体组装受损、预先形成的二聚体解离，还是两者兼有。这一观察凸显了阐明G6PDi-1抑制分子基础的必要性，以支持类似物的优化和未来抑制剂的设计。

结合位点的预测与分子动力学模拟

基于分子对接和MD模拟的结合全景分析：为探究G6PDi-1介导的抑制机制，研究人员首先使用SiteMap工具基于几何和物理化学标准鉴定了G6PD上的潜在空腔，将G6PDi-1对接到六个鉴定出的空腔中，并使用MM/GBSA方法对顶级构象进行重打分。随后对六个对接复合物进行全原子MD模拟，每个复合物运行四个独立重复，每段约5微秒。模拟揭示G6PDi-1以高度动态的方式与G6PD相互作用，在数十至数百纳秒内发生配体解离事件，并可重新结合至相同或蛋白质表面的其他区域。基于残基与配体的接触频率，定义了七个不同的空腔(C1-C7)：C1和C2分别对应催化性NADP⁺结合位点(NADP⁺-C)和结构性NADP⁺结合位点(NADP⁺-S)附近区域；C3代表位于二聚体界面对侧的远端口袋；C4标记二聚体界面结合区域；C5、C6和C7为间歇性位点。然而，位点间的频繁转换阻碍了仅基于接触分析或轨迹目视检查来识别单一主导结合口袋。

马尔可夫状态模型分析结合动力学

HMSM揭示的结合动力学特征：由于空腔间的广泛转换，研究人员使用PyEMMA包构建了隐马尔可夫状态模型(HMSM)。模型验证通过隐式时间尺度(ITS)分析和Chapman-Kolmogorov(CK)检验完成。最终构建的八态HMSM显示，从各空腔解离的平均首次通过时间(MFPT)约为20微秒，时间尺度相对相似；而结合动力学则显示位点依赖性，其中二聚体界面空腔C4表现出最快的结合速率(MFPT约136微秒)，其次是C1(约172微秒)和C3(约251微秒)，而结构性NADP⁺位点的结合则慢得多(约3毫秒)。稳态种群分析表明，C3和C4是两个主要的配体结合态，种群分别约为11.7%和9.8%，其次是C1(7.2%)。基于MFPT和稳态种群的综合分析，选择C1(催化性NADP⁺位点)、C3(远端口袋)和C4(二聚体界面)作为动力学上优先的结合区域进行进一步分析。

结合亲和力与相互作用特征

residence时间和结合自由能分析：研究人员测定了G6PDi-1在三个优先位点的平均 residence时间。结果显示，G6PDi-1在催化性NADP⁺位点(C1)和二聚体界面(C4)的停留时间分别约为27和29纳秒，而远端口袋(C3)的占据时间稍长(约33纳秒)。随后对结合帧进行聚类分析，发现G6PDi-1在C1和C4呈现异质性构象，与局部蛋白质柔性和溶剂暴露性质一致；而C3的结合则与较窄的配体位姿分布相关，表明局部环境更为受限。MM/PBSA结合自由能计算显示，尽管C4的配体波动较大且停留时间较短，但二聚体界面口袋(C4)对G6PDi-1具有最有利的结合能，其次是C1，而C3的相互作用能相对较弱。残基相互作用分析表明，疏水接触和间歇性氢键(H-键)在所有三个位点稳定配体：C1位点中Tyr249和Lys171参与氢键形成；C3位点主要由脂肪族残基(如V303、Leu468、Pro477、I480)形成的疏水空腔构成，并与Leu305主链形成一致的氢键；C4位点的结合则由疏水相互作用(如Ile220、Phe373、Leu390、Met404、Leu420)和氢键(涉及Asn218、Asn388、Asp389、Thr402和Ser418)共同驱动，产生了最有利的总体结合能。

别构效应的分子机制探究

配体结合对蛋白内部通讯和寡聚化能力的影响：为探究配体结合可能的别构调节机制，研究人员比较了配体结合态(C1、C3和C4)与活性四聚体态(含底物G6P和NADP⁺)的构象动力学差异，并进行了apo单体G6PD的MD模拟(四个重复，各1微秒)。活性形式的G6PD显示出比单体状态更窄的界面RMSD分布和稳定的亚基间接触面积。原生接触分析鉴定了稳定界面的亚基间盐桥，包括D228-R219、D375-R219、E206-K407和E419-R427。这些残基与单体状态中显示增强柔性的区域重叠，表明寡聚化通过特定的亚基间相互作用稳定界面。与之相应，G6PDi-1结合系统在形成亚基间接触和二聚体界面的残基处表现出RMSF增加，反映了由于亚基间堆积破坏所导致的增强柔性。

归一化线性互信息(nLMI)分析显示，C1和C3结合态相比于活性状态表现出更广泛的非对角耦合，差分图以正值为主，表明配体结合增强了残基间耦合；而C4结合态则呈近似中性轮廓，仅显示有限和局部的变化，与界面阻塞的直接机制一致。通过计算界面形成残基的DRMSD分布，发现apo单体约27%的帧采样了类二聚体界面构象，而这一比例在配体结合态降至4-7%，支持了G6PDi-1减少二聚体胜任性单体构象种群的机制假设。选择性高相关的长程nLMI残基对的结构映射进一步提示了从C1和C3配体相互作用残基(A217、N218、F381、V400、D421和Y424)到二聚体界面的可能别构通路，表明这些位点的扰动可能传播至参与二聚体界面堆积的区域。因此，研究结果支持位点间的功能区分：C4位点的配体结合直接干扰寡聚化，而C1和C3位点的配体结合则与增强的残基间相关性相关，促使构象集合偏离二聚体胜任性状态。

讨论与结论

研究人员在讨论部分指出，本发现提供了单体G6PD上配体可及状态的机制性描述，而非G6PDi-1依赖性寡聚状态再分布的完整热力学模型，因此在解释结果时需考虑若干方法学限制。使用单体形式虽然实现了对蛋白质表面G6PDi-1的广泛采样和动力学及能量上优先相互作用区域的识别，但无法直接定量细胞中配体诱导的单体-寡聚体平衡偏移，也不能解答G6PDi-1是否能结合并破坏已存在的G6PD二聚体。此外，停留时间和MFPT应作为比较性描述符解释，因其基于距离定义的边界态和简化拟合程序推导而来，且较慢的事件可能在可及模拟时间尺度内无法被定量解析。MM/PBSA和MM/GBSA等端点能量方法对所选择的构象集合敏感，且未明确包含与配体扩散、大规模构象变化或寡聚化平衡相关的熵贡献。

研究结论总结如下：研究人员结合生化数据、MD模拟和结合自由能计算，绘制了别构G6PD抑制剂G6PDi-1的结合全景图并定义了其可能的作用机制。G6PDi-1在HepG2细胞中减少活性二聚体形式的实验表型，在计算上得到了以下支持：优先与二聚体界面区域相互作用、界面形成残基的柔性增加、以及二聚体兼容性单体构象种群的减少。值得注意的是，G6PDi-1在远端位点的结合重塑了残基间耦合网络，揭示了额外的别构调节模式。这些发现共同揭示了G6PDi-1抑制的结构和机制基础，为新型G6PD抑制剂的基于结构的开发及有效的抗肿瘤治疗奠定了 foundation。