本文发表于 *ACS Catalysis*，作者为 Rakesh K. Roy、Dimitri Antoniou 与 Steven D. Schwartz，研究机构为亚利桑那大学化学与生物化学系。该研究于 2026 年 3 月 19 日投稿，2026 年 4 月 27 日修订，2026 年 4 月 28 日接收，最终于 2026 年 5 月 19 日发表。

本研究的核心科学问题是酶催化中的快速运动如何影响反应动力学。类固醇 5α-还原酶 2 型（steroid 5α-reductase type 2, SRD5A2）催化 NADPH 依赖的睾酮（testosterone）还原为二氢睾酮（dihydrotestosterone, DHT），该反应通过连续的氢负转移（hydride transfer）和质子转移（proton transfer）完成。在此前的研究中，人们已经知道酶的静态结构对催化至关重要，但蛋白质中从飞秒到毫秒尺度的动态运动同样可以直接影响反应速率。尤其是速率促进振动（rate-promoting vibrations, RPVs）被认为能够通过压缩反应中心、增强静电稳定或预组织催化残基来降低能垒。然而，SRD5A2 的活化能垒、速率促进动力学以及静电贡献此前均未被表征。更重要的是，此前对 RPVs 的研究主要集中于可溶性酶，对于脂双层膜结合蛋白中这类效应是否仍然存在、是否会被膜环境削弱尚不清楚。本研究旨在利用量子力学/分子力学（QM/MM）模拟与转变路径采样（transition path sampling, TPS）方法，揭示 SRD5A2 催化过程中的分子运动如何影响反应动力学，并通过与疾病相关突变体 L224P 的比较，验证动态-静电耦合在催化中的重要作用。

在方法学上，本研究首先从 RCSB 蛋白质数据库获取 SRD5A2 的晶体结构（PDB ID: 7BW1），将其中的抑制剂非那雄胺（finasteride）替换为睾酮，构建酶-底物复合物体系。该体系被嵌入 POPC 脂双层中，并在 TIP3P 水盒子中溶剂化，同时加入 0.15 M KCl 以模拟生理盐浓度。QM 区域包含睾酮、Tyr91、Glu57 以及 NADH 的烟酰胺环，共 100 个原子，使用 PM3 半经验方法处理；其余部分使用 CHARMM36 力场描述。体系依次经过能量最小化、逐步加热、1 ns 平衡和 40 ns 无偏平衡。L224P 突变体体系采用相同流程制备。为了生成无偏反应轨迹，研究者首先利用定向谐振力构造了一条有偏反应轨迹，随后以 Metropolis Monte Carlo 方式在轨迹空间中进行随机行走，生成 280 条无偏反应轨迹。在此基础上，研究者计算了氢负转移和质子转移步骤的自由能曲线，反应坐标定义为供体-受体距离的反对称组合。自由能曲线通过 12 个重叠窗口的射击算法采样获得，每个窗口累积 1000 条短轨迹。电场计算则通过叠加 QM 区域原子和 12 Å 截断范围内经典蛋白质原子的静电贡献，获得投影在睾酮羰基 C=O 键偶极上的电场强度。此外，研究者通过系统的承诺子分析（committor analysis），在约束 QM 区域的同时逐步加入不同的蛋白质残基，寻找能够使承诺子概率分布在 p_B = 0.5 附近达到最优的最小残基集合。为了进一步确定残基间的协同关系，研究者对 8 个候选残基的所有 70 种四残基组合进行了组合承诺子分析，并计算了所有 28 个残基对的成对重要性评分。整个分析流程中，研究者使用了自行开发的代码进行自由能计算和窗口拼接，同时利用自举法（bootstrap）进行误差估计。

研究结果首先体现在自由能垒方面。野生型酶中，氢负转移步骤的自由能垒为 7.8 ± 0.24 kcal/mol，质子转移步骤为 5.4 ± 0.12 kcal/mol。相比之下，L224P 突变体的氢负转移能垒升高至 11.1 ± 0.28 kcal/mol，质子转移能垒升高至 8.8 ± 0.26 kcal/mol，分别增加了 3.3 和 3.4 kcal/mol。与此同时，野生型中两个化学步骤之间的时间间隔为 35–45 fs，而突变体中这一间隔延长至 115–130 fs，超过野生型的两倍。这些结果表明 L224P 突变同时破坏了氢负转移和质子转移的过渡态稳定性。在电场分析中，野生型酶表现出与压缩事件精确相关的电场增强现象。反应物态时电场基线约为 −170 至 −185 mV/cm，在约 55 fs 时 Tyr33 对烟酰胺环的压缩使电场增强至 −209 mV/cm；在约 95 fs 时，Leu224 的推动达到最大，电场也达到峰值 −221 mV/cm。此后电场维持在约 −210 mV/cm 直至过渡态。这种约 50 mV/cm 的电场增强根据 δΔG‡ ≈ −μ·δE 的关系，可以对应 2–4 kcal/mol 的能垒降低，与突变体中观察到的能垒升高幅度高度一致。突变体的电场则几乎保持不变，维持在 −185 至 −195 mV/cm 之间，没有出现与压缩相对应的动态增强。此外，残基水平的电场分解显示，突变体中 Arg94 的稳定贡献降低 6.1 mV/cm，Trp53 和 Cys119 从主要贡献者中消失，而 Asp164、His90、Arg114 和 Arg227 则以相反的电场方向出现，进一步破坏氢负转移过渡态的静电稳定。

承诺子分析揭示了更为深层的反应坐标网络。对于氢负转移步骤，仅约束 QM 区域时，野生型的承诺子分布严重偏向产物侧，表明仅靠化学坐标无法定义过渡态。当额外约束 Tyr33、Trp53、Arg94、Cys119、Glu197、Phe223、Leu224 和 Arg227 八个残基后，分布显著拓宽并围绕 p_B = 0.5 中心化，说明这些残基的蛋白质运动是反应坐标的必要组成部分。在组合承诺子分析中，最优四残基组合为 Leu224–Tyr33–Phe223–Glu197，其中心区比例达 42%，平均承诺子概率为 0.496，几乎与理想值 0.5 完全一致。成对重要性分析显示 Glu197–Phe223 是最重要的残基对，综合重要性评分达 0.96，其次是 Tyr33–Phe223（0.92）、Tyr33–Arg227（0.88）和 Trp53–Leu224（0.85）。Leu224 虽然不属于过渡态核心残基，但通过与多个核心残基形成高重要性连接而成为一个次级枢纽。突变体中，同一四残基组合的承诺子分布严重偏向反应物侧，中心区比例仅为 6%，平均承诺子概率为 0.312，表明该组合已无法稳定过渡态。突变体中重新寻找最优组合后，得到组合 55（Cys119–Phe223–Arg94–Glu197）和组合 52（Cys119–Arg227–Phe223–Arg94），其中 Cys119、Phe223 和 Arg94 是突变体过渡态网络的不变核心。值得注意的是，Leu224 完全从突变体最优组合中消失，Tyr33 也被排除出过渡态核心，而 Cys119 从野生型中的弱贡献者转变为突变体中的主要稳定残基。

对于质子转移步骤，野生型中仅约束 QM 区域同样不足以定义过渡态，需要额外约束 Trp53、Cys119、Leu224 和 Arg227 才能获得良好的承诺子分布。然而突变体中仅约束 QM 区域就出现了接近 p_B = 0.5 的分布，说明突变体的质子转移过渡态对蛋白质坐标的依赖性降低，过渡态区域更加弥散。结构分析显示，L224P 突变导致 Arg227 向 NADPH 磷酸骨架移动，破坏了 Tyr33–Arg227 氢键，使 Tyr33 与 NADPH 的 CZ–NC4 距离从 3.24 Å 增大至 5.5 Å。同时，Trp53 侧链从野生型中与 Glu57 近乎平行的闭合构象转变为突变体中近乎垂直的开放构象。角度分布分析显示，野生型中 Trp53 吲哚环与 Glu57 羧酸根的矢量夹角最可几值约为 9°，而突变体中该角度变化至约 108°，定量证实了 Trp53 的闭合-开放构象转变。

本研究的核心结论是，SRD5A2 的催化效率依赖于压缩动力学、残基网络和电场调节三者的整合。轨迹分析观察到的两次顺序压缩事件——Tyr33 在约 55 fs 对烟酰胺环的压缩和 Leu224 在 80–120 fs 的推动——并非偶然运动，而是反应坐标的必要组成部分。承诺子分析证明，缺少这些残基时过渡态无法被正确定义。同时，电场计算显示这些压缩运动在降低供体-受体距离的同时，也动态增强了投影在底物羰基上的电场。这一电场增强约为 50 mV/cm，对应 2–4 kcal/mol 的能垒降低，与 L224P 突变体中因失去压缩而升高的 3.3–3.4 kcal/mol 能垒高度吻合。该研究还发现，与乳酸脱氢酶等可溶性酶相比，SRD5A2 中供体-受体压缩幅度较小，电场增强也相对较弱，这可能与脂双层环境的阻尼作用有关。但即便在膜环境中，快速促进振动和静电预组织仍然在催化中发挥重要作用。

本研究的价值在于将速率促进振动的概念首次系统性地拓展到膜结合酶体系。此前对于酶动态学的研究多集中于水溶性酶，而 SRD5A2 作为位于内质网膜中的整合膜蛋白，为理解膜环境对酶动态-静电耦合的影响提供了独特模型。此外，L224P 是临床上已知的导致 5α-还原酶缺乏症的致病突变，本研究从分子层面阐明了该突变如何通过破坏 Leu224 压缩、消除电场增强并重塑过渡态稳定网络来废除酶活性。这一发现不仅具有基础科学意义，也为理解 SRD5A2 相关疾病的分子病理机制以及设计下一代 5α-还原酶抑制剂提供了重要的理论依据。