可电离光响应光开关eDHFR抑制剂质子化耦合光响应结合的多尺度自由能方法

《Journal of Chemical Information and Modeling》：Multiscale Free-Energy Methods for Protonation-Coupled Light-Responsive Binding of Ionizable Photoswitchable eDHFR Inhibitors

　　

光开关配体通过以异构体依赖、光响应的方式结合靶标，实现对生物分子活性的光控。近期可电离光开关配体的发展拓展了应用，但也带来主要设计挑战：光响应结合可依赖于异构体形式、化学取代基以及结合诱导的质子化平衡偏移，这些效应紧密耦合、量级细微且难以预测。因此，少有计算方法被开发并系统基准化以定量预测上述行为。研究人员建立一种多尺度自由能方法，并以近期开发的一系列大肠杆菌二氢叶酸还原酶（Escherichia coli dihydrofolate reductase, eDHFR）可电离光开关抑制剂为模型系统予以基准测试。恒pH副本交换分子动力学（constant pH replica-exchange molecular dynamics, pH-REMD）与量子力学/分子力学（quantum mechanics/molecular mechanics, QM/MM）伞形采样模拟定量刻画了配体结合eDHFR活性位点时质子化态的变化。研究人员评估了采用不同炼金路径、热力学循环与质子化态指派的热力学积分（thermodynamic integration, TI）模拟，用以预测光响应亲和度差值及取代基效应。对每种配体直接进行cis-to-trans转化并在显式处理环境依赖质子化态时，最佳复现实验趋势；化合物间转化路径可靠性较低，因为力场误差引入难以修正的pKa偏差，且质子化/去质子化过程隐式进入热力学循环。忽略结合诱导质子化态变化的TI模拟无法一致复现实验趋势。蛋白-配体与配体-水相互作用分析进一步揭示了异构体依赖结合的能量与结构起源。该研究提出了面向靶向eDHFR可电离光开关的多尺度自由能框架，并强调预测光响应结合行为时显式考虑结合诱导质子化态变化的重要性。

研究背景方面，分子光开关为以光可逆调控生物分子系统提供强有力工具，光药理学中光开关配体可实现非侵入、远程、选择性调控。合理设计光开关配体需准确预测三量：给定靶标的活性异构体、活性与惰性异构体间光响应结合自由能差、化学取代对该差值的影响。异构体结构差异常局域于异构化键附近，结合自由能差通常低于几kcal/mol，标准分子对接与端点自由能法精度不足。对于含可滴定基团的可电离光开关，其优势质子化态在体溶液与蛋白口袋间可能因局部静电环境而偏移，使光响应结合同时依赖几何差与蛋白-配体相互作用对质子化平衡的扰动。已有实验合成TMP（trimethoprim）衍生偶氮苯光开关抑制剂化合物6、11、15靶向eDHFR，其中羧基与氯取代显著改变绝对抑制力与cis-on效应，但缺乏定量自由能框架解释取代基如何通过质子化耦合影响光响应亲和差。

研究人员以eDHFR为模型，建立多尺度自由能方法并基准化。体系源于PDB 3DAU晶体结构，配体cis/trans异构体经AutoDock Vina对接入底物口袋，关键二氨基嘧啶环与Asp27相互作用被保留。力场采用AMBER ff14SB与GAFF，偶氮键周边三个二面角经FOMO-hh-TDA-DFT多参考电子结构方法重参数化以修正GAFF过低异构化势垒。

关键技术方法上，研究人员依次采用恒pH副本交换分子动力学（pH-REMD）判定水溶液与蛋白结合态优势质子化态及pK_a偏移；以B3LYP-D3/MM伞形采样（umbrella sampling, US）验证Asp27-配体质子转移自由能面；以热力学积分（TI）沿两条炼金路径（Method A单配体cis→trans、Method B同异构体化合物间转换）结合不同质子化态协议计算光响应结合差与取代基效应；以MM-GB(PB)SA端点法作对照；常规MD与轨迹分析用于蛋白-配体及配体-水相互作用分解。

研究结果部分，Overview and Rationale for the Multiscale Computational Framework小节指出多尺度框架必要，因需同时描述质子化平衡、蛋白-配体结合与异构自由能变化。System Setup小节中，eDHFR（PDB 3DAU）对接 poses 保留二氨基嘧啶-Asp27关键氢键，力场二面角经多参考QM重参。Force-Field-Based MD Equilibration Simulations小节通过能最小化与NVT/NPT平衡获得150 ns无约束轨迹用于相互作用分析。pH-REMD Simulations小节对化合物6显示：水溶液中配体pK_a低，结合eDHFR后pK_a显著升高，Asp27全程去质子化，表明结合诱导N1质子化。QM/MM US小节以δ=d(O–H)–d(N–H)为反应坐标，PMF最低谷位于质子在配体（δ≈0.5 ?），即质子化配体/去质子化Asp27优于反向，验证pH-REMD结论。TI基准小节表明Method A（显式cis→trans、环境依赖质子化态）最佳复现实验cis-on趋势；Method B化合物间路径受力场pK_a误差污染；忽略结合诱导质子化的TI协议一致失败。MM-GB/PBSA不能稳定复现趋势。Interaction Analysis小节中，cis-6相对trans-6在eDHFR内具更强静电作用（二氨基嘧啶-Asp27盐桥），vdW贡献异构体依赖；配体-水相互作用在cis态减弱，利于结合。Figure 7 poses显示异构体改变偶氮取向从而重排口袋接触。

讨论与结论部分，研究人员总结：离子化光开关的光响应亲和差不能仅由几何异构解释，必须显式计入结合诱导质子化态变化；直接cis-to-trans TI配合pH-REMD/QM-MM界定质子化态是最可靠定量路径；化合物间TI因力场pK_a误差放大而不稳；端点法不适于细微光响应差预测。该多尺度自由能框架发表于《Journal of Chemical Information and Modeling》，为光药学可电离光开关的理性设计提供了首个经实验基准化的质子化耦合预测方案，明确了Asp27与二氨基嘧啶N1质子共享是eDHFR系统中cis-on效应保留或削弱的结构能量根源，对后续TMP类光开关取代基优化具有直接指导价值。