: 糖胺聚糖如何介导组织蛋白酶K原的成熟：来自计算研究的新见解

《Biophysical Chemistry》：How glycosaminoglycans mediate procathepsin K maturation: New insights from computational studies

　　

组织蛋白酶主要是在溶酶体和细胞外基质（ECM）中发挥功能的半胱氨酸蛋白酶，它们调节着重要的蛋白水解过程。它们被合成为无活性的酶原（组织蛋白酶原），其中N端肽段阻挡了活性位点的通路并在成熟过程中被移除。已知糖胺聚糖是一类由重复二糖单位组成的线性硫酸化多糖，能够调节组织蛋白酶的活性和酶原加工。在此，研究人员使用分子建模来阐明GAG在组织蛋白酶K原成熟过程中的作用。分子对接、分子动力学（MD）模拟和结合自由能计算识别了组织蛋白酶原表面上介导GAG识别的推测别构位点。对最稳定复合物进行的微秒级MD模拟，通过主成分分析（PCA）进行分析，揭示了GAG依赖性的酶原构象景观变化。为了探究环境效应，研究人员进一步在溶酶体条件（pH 4）下模拟了apo和结合GAG的组织蛋白酶K原。这些模拟表明，酸化和GAG结合在促进激活方面存在协同作用。低pH值破坏了肽段的二级结构，而GAG在一个独特的别构热点结合则放大了这种效应，促进了螺旋解旋并有助于肽段解离。总体而言，研究结果为GAG辅助的组织蛋白酶K原成熟提供了分子水平的框架，突出了酶原激活的协同环境和别构调节。

研究背景与意义
组织蛋白酶主要是在溶酶体和细胞外基质（ECM）中发挥功能的半胱氨酸蛋白酶，参与细胞稳态和细胞外基质重塑等生理过程。其结构采用典型的木瓜蛋白酶样折叠，由L和R两个亚结构域组成V形活性位点裂隙，催化三联体Cys¹²⁴、His²⁶¹和Asn²⁸¹位于结构域界面。组织蛋白酶K（Cathepsin K, CatK）是骨吸收的核心酶，负责降解I型和II型三股螺旋胶原蛋白。其胶原酶活性高度依赖于与糖胺聚糖（GAGs，一类高度阴离子化的线性硫酸化多糖）形成高分子量复合物。由于CatK功能失调与骨硬化症、骨质疏松症、类风湿性关节炎等多种病理密切相关，靶向其别构外位点以破坏CatK-GAG复合物形成已成为避免活性位点抑制脱靶效应的替代治疗策略。与其他半胱氨酸组织蛋白酶一样，CatK被合成为无活性的酶原（Zymogens，即procatK），其N端肽段空间阻断活性位点。成熟过程通过酸性内溶酶体条件下的自催化处理或其他内吞蛋白酶的蛋白水解裂解发生。尽管已知GAG能够调节组织蛋白酶活性和酶原加工，但GAG介导的组织蛋白酶K原成熟的分子机制尚不清楚，尤其缺乏计算层面的深入理解。为此，研究人员系统地筛选了不同链长、电荷和硫酸化模式的GAG，以确定稳定的procatK复合物及其结合位点，并在溶酶体条件下模拟评估酸化与GAG结合的相互作用。该研究发表在《Biophysical Chemistry》，揭示了GAG结合和溶酶体酸化在促进procatK成熟中的协同作用，为酶原激活的协同环境和别构调节提供了分子水平的框架，对开发针对CatK相关疾病的别构药物具有重要意义。

关键技术方法
研究人员从蛋白质数据库（PDB）获取了procatK（PDB ID: 1BY8）和CatK（PDB ID: 4X6H）的晶体结构作为研究基础。主要关键技术方法包括：采用泊松-玻尔兹曼表面积（PBSA）计算静电势等值面以预测GAG结合区域；使用AutoDock 3进行分子对接确定结合模式；运用AMBER21在中性及酸性（pH 4）条件下进行分子动力学（MD）模拟及微秒级模拟以解析构象动态；通过分子力学广义Born表面积（MM-GBSA）方法计算结合自由能评估复合物稳定性；利用主成分分析（PCA）、均方根涨落（RMSF）、DSSP分析及断开图等手段解析变构调节和构象变化。

研究结果
3.1. 静电势表面描述
研究人员通过PBSA分析识别了proCatK上可能介导与阴离子GAG结合的正电荷区域。结果显示，与成熟酶相比，proCatK具有更明显的带正电荷表面斑块，最大的区域位于R亚结构域。此外，在肽段内部及附近也存在较小的正电区，表明GAG结合并不局限于成熟结构域。

3.2. proCatK-GAG结构的分子描述
通过分子对接，研究人员发现了复发的GAG结合区域。位点I位于R亚结构域背面，是主导结合区域，与实验晶体结构中观察到的结合区一致。研究还发现GAG链长和净电荷与结合位点偏好相关，较短的dp₄链优先采样多个位点，而延伸至dp₆则使结合转向R亚结构域背面（位点I），位点I占据了超过80%的dp₆簇。

3.3. procatK-GAG复合物稳定性评估
通过结合自由能计算和RMSD分析评估了结合位点的稳定性。位点I是最主要的结合位点，高度硫酸化的GAG（如HP和C462S/C463S）在此形成最稳定的复合物，其ΔG_bind值低于-100 kcal mol^-1。稳定性分析表明，GAG链所带负电荷越多，结合自由能越有利，证实了静电互补性在procatK-GAG相互作用中的主导作用。

3.4. 选定procatK-GAG复合物的变构调节分析
对最稳定复合物的微秒级MD轨迹进行PCA分析表明，在中性pH下未观察到肽段在活性位点区域的显著重排，但揭示了远端肽段区域存在GAG依赖性的构象变化，如α-螺旋局部解旋。在pH 4的溶酶体条件下，酸化显著增加了第一主成分的方差，表明肽段集体运动增强。RMSF分析显示，酸化和C4S结合协同增强了肽段N端区域的灵活性。DSSP分析进一步证实，酸性条件下肽段α-螺旋的二级结构发生显著丢失。

3.5. 成熟机制见解
通过分析C4S结合的procatK复合物构象变化及断开图，研究发现中等长度α-螺旋发生明显重排，包括部分解旋和横向移动，同时肽段从核心酶上脱离，增加溶剂暴露以促进裂解。关键的是，肽段构象变化的启动与C4S在成熟结构域表面特定别构热点（涉及Gly²¹¹和Asn²⁹⁷残基）的结合密切相关。GAG在此位点的结合充当变构开关，触发信号显著增加肽段回转半径（R_g）并引发展开。

讨论与结论总结
讨论部分指出，本研究识别出的结合位点与先前实验晶体结构及计算预测的别构位点高度重叠，构成更广泛的GAG结合景观。模拟结果深入阐明了proCatK成熟的机制：低pH单独破坏二级结构并增加肽段移动性，而C4-S结合显著放大这一效应，通过肽段α-螺旋解旋及其从成熟结构域脱离来推进成熟，从而增加肽段对溶酶体蛋白酶的可及性并加速激活。结论部分总结如下：组织蛋白酶在生物系统中发挥多样且关键的功能，别构调节因能避免直接活性位点抑制的并发症而成为极具吸引力的治疗策略。GAG被广泛认为是组织蛋白酶活性的别构调节剂。本研究探讨了负责骨吸收的主要酶CatK的酶原形式与各种GAG的复合物。能量和几何分析确定了最有利的结合位点，高硫酸化的GAG链形成最稳定的复合物。模拟提供了由酸性和GAG结合协同效应驱动的proCatK成熟的机制见解。低pH单独破坏二级结构，而C4-S结合显著放大了这种效应，成熟过程通过肽段中α-螺旋的解旋和随后从成熟结构域的脱离来进行。这增加了肽段对溶酶体蛋白酶的可及性，从而促进裂解并加速激活。研究人员在成熟结构域表面识别出一个关键别构热点（涉及Gly²¹¹和Asn²⁹⁷残基）。GAG结合到这个口袋充当变构开关，触发使肽段回转半径显著增加并引发展开的信号。研究人员提出proCatK成熟可能发生在更长的时间尺度上，可能需要超过1μs的模拟或增强采样方法。未来的工作应该研究负电荷较少的链对成熟的影响，假设较不稳定的复合物可能赋予糖链更大的构象灵活性，从而增加定位到别构口袋的概率。