《RSC Advances》：Computational insights into the electrostatic origins of mechanistic selectivity: human topoisomerase IIα enforces ATP hydrolysis via substrate-assisted catalysisOpen Access
编辑推荐：
人类DNA拓扑异构酶IIα（topo IIα）是一种关键的分子马达，也是重要的抗癌靶点。尽管ATP水解对其功能至关重要，但底物辅助催化（substrate-assisted catalysis, SAC）相对另一种被提出的通用碱催化（general base
人类DNA拓扑异构酶IIα（topo IIα）是一种关键的分子马达，也是重要的抗癌靶点。尽管ATP水解对其功能至关重要，但底物辅助催化（substrate-assisted catalysis, SAC）相对另一种被提出的通用碱催化（general base catalysis, GBC）通路为何具有机制上的优先性，其根源仍是一个悬而未决的问题。研究人员采用密度泛函理论（Density Functional Theory, DFT）结合点电荷静电嵌入方法，对先前由QM/MM模拟推导出的反应通路进行了研究，从而明确分离出酶预组织的静电效应。该分析揭示了SAC通路相对于GBC通路存在高达约40 kcal mol?1的静电偏好，这一现象被称为机理选择性。逐残基静电分解表明，仅16个氨基酸就贡献了该选择性的绝大部分。此外，库仑相互作用和自然键轨道（Natural Bond Orbital, NBO）分析表明，这一静电网络通过调节反应中心的电荷分布和特定轨道占据数来执行机理选择性。最终，识别出在静电上促进ATP水解的氨基酸，为未来开发能够绕过传统拓扑异构酶毒剂遗传毒性的新一代催化抑制剂提供了机理上的洞见。
人类DNA拓扑异构酶IIα（human DNA topoisomerase IIα, topo IIα）是IIA型拓扑异构酶家族成员，在增殖细胞中高表达，长期以来被用作抗癌化疗的靶点。以依托泊苷（etoposide）和多柔比星（doxorubicin）为代表的拓扑异构酶II毒剂通过稳定DNA与topo II之间的瞬时共价复合物，将酶转化为细胞毒素，但伴随严重的心毒性及继发性恶性肿瘤等副作用。因此，研究人员致力于开发能够在不可起DNA损伤的前提下抑制酶活性的替代化合物——催化抑制剂。靶向ATP结合位点是其中的重要策略之一。IIA型拓扑异构酶为同源二聚体，含N端ATPase结构域、中央DNA-gate和C-gate三个功能区域，其中ATP结合稳定N-gate二聚化并促进T-段捕获，而ATP水解对于链通过、构象重置和催化循环完成至关重要。
研究人员此前通过多尺度QM/MM计算（replica path method, RPATh）对topo IIα中ATP水解机制进行了探索，发现该反应经由解离机制进行，优先遵循底物辅助催化（substrate-assisted catalysis, SAC）通路，即ATP的γ-磷酸作为裂解水的质子受体；该机制优于涉及QTK环Lys378的通用碱催化（general base catalysis, GBC）模型，突变研究亦支持Lys378主要起稳定而非催化作用。SAC机制的活化势垒约为22 kcal mol
?1，而GBC高达52 kcal mol
?1，后者在生命过程中极不典型。然而，酶对SAC机制强烈偏好的物理根源尚不清楚。鉴于酶催化能力主要源于预组织静电的假设日益获得支持，研究人员采用聚焦的QM/点电荷方法，严格分离酶环境的静电贡献，以探究topo IIα静电环境是否及如何在机制层面决定这一偏好。该工作不仅有助于理解酶催化的基本物理基础，也为设计靶向ATP结合位点、避免DNA损伤的新型催化抑制剂提供亚原子层面的指导。
在技术方法上，研究人员以实验解析的topo IIα与AMP-PNP复合物结构（PDB: 1ZXM）为起点，通过RPATh方法分别生成SAC和GBC通路的16点最小能量路径。反应中心包含ATP的甲基三磷酸片段、Glu87和Lys378侧链、裂解水、Mg
2+阳离子及其配位水分子；酶环境（除Glu87和Lys378外的全部氨基酸）以CHARMM26力场原子点电荷表示。对每个反应坐标点，在有无点电荷两种条件下进行M06-2X/6-31+G(d,p)泛函的单点DFT计算，通过能量差ΔE量化酶静电对反应中心势能的影响；进一步对SAC和GBC的R与TS结构计算活化能变化ΔE
a，并以双差量ΔΔE
a量化机制间偏好。此外，采用库仑定律计算反应中心与酶环境间的成对静电相互作用，结合自然键轨道（NBO）分析提取原子电荷和键序，将机理选择性追踪至亚原子层面。
3.1 Topo IIα静电仅对SAC机制起催化作用：将酶静电嵌入前后DFT能量剖面进行比较，发现topo IIα点电荷稳定SAC通路每个反应坐标点，尤其对RC点9–16产生超过?35 kcal mol
?1的强稳定化；而无点电荷时SAC和GBC的活化能分别为57.1和45.6 kcal mol
?1，加入酶静电后SAC的E
a降低逾37 kcal mol
?1，而GBC反而升高近2 kcal mol
?1。这表明topo IIα的静电环境对SAC机制是强催化的，对GBC机制则轻微反催化，产生约40 kcal mol
?1的机制间能量判别。
3.2 SAC与GBC机制催化效应随距离的累积模式截然不同：通过逐步扩展距反应中心的截断半径并累加所包含氨基酸的ΔE
a值，发现SAC机制中强催化效应在内部15 ?内稳定累积，远端贡献基本抵消；而GBC的累积剖面高度无序、围绕零值反复波动，类似无序溶剂效应，表明topo IIα在静电上仅对SAC机制有进化适应。
3.3 少数氨基酸贡献绝大部分静电催化效应：按残基分解发现，|ΔE
a| ≥ 1.6 kcal mol
?1的18个氨基酸贡献了SAC机制?38.1 kcal mol
?1总催化效应的97%以上；其中Lys168贡献最大（ΔE
a = ?8.2 kcal mol
?1），实验显示Lys168Ala突变使催化活性降低7倍。GBC机制总净效应为+3.6 kcal mol
?1，催化与反催化残基相互抵消。
3.4 16个特定残基驱动机理选择性：通过ΔΔE
a分析识别出16个机理分化氨基酸（|ΔΔE
a| ≥ 3.0 kcal mol
?1），其贡献超过99%的SAC偏好。其中Lys83、Asp153和Arg242更利于GBC，而Asp86、Ser148、Asn150、Lys157、Gly161、Arg162、Gly164、Tyr165、Gly166、Lys168、Glu379、Lys37*和Arg46*更利于SAC。10个选择性残基集中于高度保守的ATP盖子（142–171）区域，Lys37*和Arg46*来自另一单体，凸显跨二聚体界面的静电协同。
3.5 关键残基在R→TS转化中结构变化极小：16个机理分化氨基酸在SAC和GBC中的RMSD分别仅为0.063 ?和0.119 ?，表明机理选择性并非源于酶结构的几何重排。
3.6 机理选择性主要由电荷重分布而非几何变化驱动：库仑分解显示，16个选择性氨基酸与反应中心的相互作用在SAC中为强催化（ΔE
Coulomb = ?32.3 kcal mol
?1），在GBC中为反催化（+14.3 kcal mol
?1）。ΔΔE
Coulomb与ΔΔQ呈现较强相关性（R2 = 0.73），表明原子在TS中变得越负，其所在机制越受静电偏好，证明酶静电通过匹配SAC特有的电荷重排模式而非几何变化来执行机理选择性。
3.7 Topo IIα静电通过调节自然键轨道选择性执行SAC通路：NBO键序分析显示，酶静电在SAC的TS中增强裂解水氧与γ-磷之间即将形成的键（+1.3%）以及质子转移相关键（+0.1%），同时削弱将被断裂的β-γ磷酸键（?0.1%）和O40–H41键（?0.3%），总键序提高5.5%；而GBC中总键序仅提高0.2%，且酶静电显著阻碍了亲核攻击所需的关键键形成（P29–O40为?3.9%），并主动抑制另一关键断裂过程。这表明机理选择性可直接追溯到酶静电对特定自然键轨道的选择性调控。
在讨论部分，研究人员指出，预组织静电不仅降低反应总体活化势垒，更充当反应通路本身的仲裁者。16个残基中ATP盖子区域的集中分布表明该结构并非简单的机械性盖子，而是通过其静电特性主动选择SAC通路。尽管存在固定点电荷、静态模型忽略熵贡献、固定质子化状态等局限，约40 kcal mol
?1的巨大能量判别使定性结论保持稳健。研究结论部分翻译如下：研究人员证实，人类topo IIα预组织的静电环境不仅仅降低ATP水解的整体活化势垒，更是反应通路本身的主要仲裁者。通过计算分离酶预组织静电，揭示了约40 kcal mol
?1的巨大能量判别，决定性地偏向SAC通路。这种被定义为机理选择性的效应主要源于亚原子层面电荷重新分布，而非几何重排；16个特定残基（以ATP盖子区域的Lys168为代表，并包含来自相邻单体的残基）主导了这一选择性，表明topo IIα的四级结构通过二聚体界面的主动协调来执行机理选择。该研究为未来设计靶向topo IIα ATP结合位点的新型催化抑制剂提供了亚原子层面的机理指导，并提出了酶静电机理选择性是否在不同酶家族中普遍存在的科学问题。
