立体化学引导的不可逆到慢结合组织蛋白酶L抑制的切换机制研究

《ACS Catalysis》：Stereochemistry-Guided Switch from Irreversible to Slow-Binding Cathepsin L Inhibition

为调控组织蛋白酶L（CatL）的抑制模式，研究人员通过结合计算与实验方法，系统探究了环氧酮类抑制剂立体构型与叠氮基功能化对其抑制机制的影响，成功实现了从不可逆到慢结合可逆抑制的精准调控，为开发可调控驻留时间的共价抑制剂提供了新策略。

　　

在生命体的复杂精密工厂里，蛋白质的合成与降解维持着微妙的平衡。组织蛋白酶L（Cathepsin L， CatL）就像是工厂里的一个关键“分解师”，作为一种溶酶体半胱氨酸蛋白酶，它负责清理废旧蛋白、重塑细胞外基质，并在抗原呈递等免疫调节中扮演重要角色。然而，一旦这位“分解师”失控，就可能引发一系列严重的健康问题。研究表明，CatL的异常活跃与癌症转移、病毒（如SARS-CoV-2）入侵细胞、自身免疫性疾病（如类风湿性关节炎、1型糖尿病）乃至肾脏疾病（如蛋白尿）等多种病理过程密切相关。因此，精准地“管制”CatL，开发高效特异的抑制剂，已成为对抗这些疾病的重要策略。

在众多抑制剂中，共价抑制剂因其能与靶点形成稳定、持久的复合物，具有作用时间长、不易产生耐药性等优势，备受关注。其中，环氧酮（epoxyketone）类化合物是一类经典的共价抑制剂“弹头”，它能与蛋白酶活性位点的亲核性半胱氨酸残基形成稳定的共价键。然而，传统的共价抑制剂往往“一锤定音”，形成不可逆的复合物，这在带来持久疗效的同时，也可能增加脱靶毒性的风险。那么，能否设计一种“智能开关”，让抑制剂的结合时间变得可调，从“永久关闭”变为“可逆的长时间暂停”呢？这正是发表在《ACS Catalysis》上的这项研究试图回答的核心科学问题。

研究人员设想，通过对抑制剂“弹头”的立体化学（即原子在空间中的排列方式）和化学结构进行精妙改造，或许能实现对抑制动力学（kinetics）的精细调控。为了验证这一设想，他们开展了一项结合计算化学与多维度实验验证的系统性研究。研究团队设计、合成了两类抑制剂：一对具有不同立体构型的二肽基环氧酮类化合物（D1和D2），以及两个相应的叠氮基衍生物（AZA和AMK）。通过经典的分子动力学（Molecular Dynamics， MD）模拟、杂化量子力学/分子力学（QM/MM）自由能计算，并结合酶动力学分析、电喷雾电离质谱（ESI-MS）、饱和转移差（STD）核磁共振（NMR）以及纳米差示扫描荧光（nanoDSF）等多种技术，他们深入揭示了立体化学和叠氮基功能化如何像“分子旋钮”一样，精确调节CatL的抑制模式。

研究者主要采用了分子动力学模拟分析抑制剂-酶复合物的构象稳定性与相互作用，利用QM/MM方法计算了共价键形成路径的自由能垒以阐明反应机理。通过酶动力学实验测定了抑制常数（K_i）和二级速率常数（k_2nd），评估了抑制效力与模式。利用电喷雾电离质谱直接检测并验证了共价酶-抑制剂加合物的形成。借助饱和转移差核磁共振技术，实时监测了非共价相互作用随时间的变化，揭示了慢结合过程。此外，纳米差示扫描荧光技术被用于评估抑制剂结合对酶热稳定性的影响。

3.1. 合成和两种二肽基环氧酮的立体化学阐明

研究首先完成了目标抑制剂D1、D2、AZA和AMK的合成，并通过实验手段明确了D1和D2的立体构型，为后续的构效关系研究奠定了基础。

3.2. 二肽基环氧酮非共价E:I反应复合物的计算研究

通过500纳秒的分子动力学模拟，比较了D1和D2与CatL活性位点的结合模式。结果显示，尽管结构相似，但二者的行为截然不同：(R)-构型的D1的“弹头”部分在活性位点中具有更高的构象柔性，能够采样到适合催化三联体（Cys25?His163?Asn187）中Cys25亲核攻击的反应性构象。相反，(S)-构型的D2则主要稳定在一个非反应性构象中，其环氧酮与Trp189形成氢键，使其远离攻击位点。有趣的是，模拟发现D2可以翻转形成一个较不稳定的“翻转”构象（D2-flipped），此时其叠氮乙酰基部分靠近Cys25，提示了另一种潜在的反应路径。结合能分析也表明，D2的非反应性构象结合更强，而D1和D2-flipped构象的结合相对较弱。

3.3. 环氧酮抑制CatL化学步骤的QM/MM研究

为了从能量角度理解反应可行性，研究者采用了QM/MM方法计算了共价键形成的自由能面。对于D1，无论亲核攻击发生在环氧环的哪个碳上，其速决步骤（共价键形成）的自由能垒均在20-23 kcal·mol^-1范围内，且形成的最终共价加合物非常稳定（ΔG低至约-40 kcal·mol^-1），反向能垒极高，这从理论上证实了D1是一种高效的不可逆共价抑制剂。对于D2-flipped构象，计算发现其通过叠氮基相邻碳受到攻击的S_N2型反应能垒高达38.8 kcal·mol^-1，预示这是一个极其缓慢的过程。

3.4. E-I复合物的质谱分析

实验数据完美支持了计算预测。将CatL与过量抑制剂孵育1小时后，电喷雾电离质谱分析显示，D1几乎完全与酶形成了共价加合物，未见游离酶。而D2则大部分以游离酶形式存在，仅检测到微量的完整D2加合物和丢失了叠氮基的加合物（CatL+D2-N₃），这与计算预测的D2具有极低的反应性但可通过翻转构象缓慢形成加合物的“双功能”特性相符。

3.5. 饱和度转移差核磁共振评估环氧酮与CatL反应的时间依赖性

STD-NMR实验进一步从动力学上验证了这一点。D2的叠氮乙酰基部分显示出随时间变化的信号，表明其与酶的结合是一个缓慢的、时间依赖的共价过程。而D1则未观察到类似的慢结合特征，符合其快速不可逆的反应机制。

3.6. CatL的动力学实验

酶动力学实验定量揭示了四种抑制剂的效力差异。D1表现出典型的时间依赖性不可逆抑制，其二级速率常数k_2nd为3.9 x 10⁴M^-1s^-1，活性很高。D2、AZA和AMK则未显示时间依赖性，其抑制曲线符合可逆抑制模型，测得了相应的抑制常数K_i。值得注意的是，D2的K_i为0.134 μM，表明其非共价结合亲和力很强，而将其环氧酮替换为酰胺（AZA）后，K_i骤降至126 μM，说明(S)-环氧酮对于维持高亲和力结合至关重要。将环氧酮替换为叠氮甲基酮（AMK）后，K_i为2.6 μM，活性得以部分保留。

3.7. CatL与AMK非共价E:I反应复合物的计算研究 & 3.8. AMK抑制CatL化学步骤的计算研究

对AMK的分子动力学模拟显示，其结合模式与D2类似，但结合稍弱。QM/MM计算预测其通过叠氮基的S_N2型反应能垒为28.2 kcal·mol^-1，虽比D2-flipped的能垒低了约10 kcal·mol^-1，但仍属于较高能垒，预示着慢反应特性。

3.9. CatL与AMK形成的E-I产物复合物的质谱和STD-NMR分析 & 3.10. CatL+抑制剂非共价复合物结合能分析

质谱实验证实，在高浓度下（酶:抑制剂=1:60），AMK能与CatL形成约50%的共价加合物，而在等摩尔浓度下则几乎不形成，确证了其慢结合、浓度依赖的特性。相比之下，其前体氯甲基酮（ClMK）在相同条件下则完全形成加合物，凸显了叠氮基替换带来的动力学调节作用。结合能计算（MM-GBSA/MM-PBSA）的结果定性地支持了各抑制剂观测到的结合亲和力趋势。

综上所述，这项研究通过精巧的分子设计和多学科交叉的研究手段，成功演示了如何通过“立体化学”和“弹头工程”这两个“分子旋钮”，对组织蛋白酶L（CatL）的抑制模式进行精准编程。研究得出的核心结论是：(R)-构型的环氧酮（D1）是一种高效、快速的不可逆抑制剂，其反应能垒适中，产物极其稳定；(S)-构型的环氧酮（D2）则展现出独特的“双功能”特性：其优势构象具有高亲和力但无反应性，而一个次要的“翻转”构象使其叠氮基部分能作为弱亲电“弹头”，以极高的能垒（极慢的速度）发生反应，从而表现为一种慢结合、时间依赖的共价抑制剂。用叠氮甲基酮（AMK）直接替换环氧酮，可以模拟这种慢结合机制，且反应性高于D2的翻转构象，而简单地用酰胺（AZA）替换环氧酮则会丧失大部分活性。

这项工作的意义重大。首先，它从原子层面深入揭示了抑制剂立体化学如何通过影响其在酶活性口袋中的构象、氢键网络和非键相互作用，最终决定其反应性和抑制动力学，为基于结构的理性药物设计提供了关键见解。其次，它成功地将CatL的抑制模式从经典的不可逆共价抑制，拓展到了慢结合可逆/共价抑制，实现了对抑制剂-靶标复合物“驻留时间”的可控调节。这种可调性对于平衡疗效与潜在脱靶毒性至关重要，尤其在需要长期给药的慢性病治疗中。最后，研究验证了叠氮基作为一种可调控的弱亲电“弹头”的潜力，为开发新型“可点击化学”探针和具有优化药代动力学特性的治疗剂开辟了新途径。这项研究不仅为针对CatL的相关疾病（如癌症、病毒感染、自身免疫病）的药物开发提供了新的设计框架和先导化合物，其揭示的“结构-动力学”关系原理也可望应用于其他靶点的共价抑制剂设计，具有广泛的启发意义。